Ligand-Receptor Interactions of Lamivudine: A View from Charge Density Study and QM/MM Calculations

The nature and strength of interactions for an anti-HIV drug, Lamivudine, were studied in a pure crystal form of the drug and the ligand–receptor complexes. High-resolution single-crystal X-ray diffraction studies of the tetragonal polymorph allowed the drug’s experimental charge density distribution in the solid state to be obtained. The QM/MM calculations were performed for a simplified model of the Lamivudine complex with deoxycytidine kinase (two complexes with different binding modes) to reconstruct the theoretical charge density distribution. The peculiarities of intramolecular interactions were compared with previously reported data for an isolated molecule. Intermolecular interactions were revealed within the quantum theory of ‘Atoms in Molecules’, and their contributions to the total crystal energy or ligand–receptor binding energy were evaluated. It was demonstrated that the crystal field effect weakened the intramolecular interactions. Overall, the energies of intermolecular interactions in ligand–receptor complexes (320.1–394.8 kJ/mol) were higher than the energies of interactions in the crystal (276.9 kJ/mol) due to the larger number of hydrophilic interactions. In contrast, the sum of the energies of hydrophobic interactions was found to be unchanged. It was demonstrated by means of the Voronoi tessellation that molecular volume remained constant for different molecular conformations (250(13) Å(3)) and increased up to 399 Å(3) and 521(30) Å(3) for the Lamivudine phosphate and triphosphate.