Azathioprine is an important drug used in the therapy of autoimmune disorders and in preventing graft rejection. Its molecule is composed of two main moieties:  mercaptopurine and imidazole derivative. It is an immunosuppressive agent whose biological activity results from its in vivo mercaptolysis mediated by a nucleophilic attack on the C(5i) atom of imidazole ring of the azathioprine molecule. Solvation model SM5.4 with the PM3 Hamiltonian have been applied to model the reaction of azathioprine with cysteine. The employed quantum mechanical method shed new light on the mechanism of the reaction of azathioprine with cysteine in aqueous solution. The obtained results indicated that the first step in the reaction most likely involves the nucleophilic attack of the COO^- of cysteine on the C(5i) atom of the imidazole ring of azathioprine, followed by a subsequent intramolecular attack of the SH group of the cysteine residue. It was shown that biogenic thiols such as glutathione or cysteine facilitate the first and crucial step of azathioprine metabolism, due to the presence of COO^-, SH, and NH₃⁺ groups in their molecules.