A simple procedure for the derivation of electron density based surfaces of drug-receptor complexes from a combination of X-ray data and theoretical calculations

Graphical abstract Electrostatic potential on an electron density iso-surface (0.0067 eÅ−3) of an APE substrate in the active site pocket of the SARS main protease Mpro. Geometry data from PDB entry 2A5I.
To contribute to an understanding of biological recognition and interaction, an easy-to-use procedure was developed to generate and display molecular surfaces and selected electron density based surface properties. To overcome the present limitations to derive electron densities of macromolecules, the considered systems were reduced to appropriate substructures around the active centers. The combination of experimental X-ray structural information and aspherical atomic electron density data from theoretical calculations resulted in properties like the electrostatic potential and the Hirshfeld surface which allowed a study of electronic complementarity and the identification of sites and strengths of drug-receptor interactions. Applications were examined for three examples. The anilinoquinazoline gefitinib (Iressa®) belongs to a new class of anticancer drugs that inhibit the tyrosine kinase activity of the epidermal growth factor receptor (EGFR). In the second example, the interaction of epoxide inhibitors with the main protease of the SARS coronavirus was investigated. Furthermore, the progesterone receptor complex was examined. The quantitative analysis of hydrogen bonding in the chosen substructure systems follows a progression elaborated earlier on the basis of accurate small molecule crystal structures. This finding and results from modified substructures suggest that also the surface properties seem robust enough to provide stable information about the recognition of interacting biomolecular species although they are obtained from medium molecular weighted subfragments of macromolecular complexes, which consist of no more than ∼40 residues.