key: cord-1037367-t6rnpnfc authors: Mebs, Stefan; Lüth, Anja; Luger, Peter title: 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 date: 2010-08-15 journal: Bioorg Med Chem DOI: 10.1016/j.bmc.2010.06.080 sha: c7005ee9025aa72654f99df95eb122cc31254c42 doc_id: 1037367 cord_uid: t6rnpnfc 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. 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. Ó 2010 Elsevier Ltd. All rights reserved. Since the pharmacological action of most drugs is promoted by a recognition process enabling a macromolecular target protein and a low molecular weight substrate to form a complex, the understanding of mutual drug-receptor recognition on an atomic scale is a key challenge in that field. Biological recognition is dominated by intermolecular interactions and it is generally accepted that not only steric but also complementary electronic properties of the involved species play a distinctive role in this process. 1, 2 Whereas steric information can be provided by various methods, for example, by the determination of the molecular structure from a conventional X-ray diffraction experiment using the spherical independent atom model (IAM ), electronic properties can be ob-tained from high-resolution X-ray data at low temperatures followed by an aspherical electron density modelling of the structure. For aspherical atom modelling, the Hansen and Coppens multipole formalism is commonly used that describes the asphericity of an atomic electron density by a sum of spherical harmonic density functions. 3 In addition to intramolecular atomic and bonding properties, experimental electron density (ED) distributions contain information on weak intermolecular interactions in the crystal. For example, the electrostatic potential can be derived, being one of the most important properties in the study of molecular reactivity and the analysis of molecular association, which follows the so called principle of electrostatic complementarity. 1,4-6 ED properties in the crystal can serve in a first approximation as a model for the interactions under biological conditions, since the same types of effects as the ligand-target interactions should exist in this environment (e.g., steric, electrostatic, hydrogen bonding, van-der-Waals effects), and these are expected to be comparable in size. 2, 7 In this respect, experimental ED properties are much better suited to simulate physiological conditions than those from theoretical gas phase calculations on isolated molecules. 2 Several experimental advances have taken place in ED research in the last decade. 8 Since the ED of a chemical structure consists of 0968-0896/$ -see front matter Ó 2010 Elsevier Ltd. All rights reserved. doi:10.1016/j.bmc.2010.06.080 a large spherical and a very small aspherical contribution, located mainly in covalent bonds and lone pair regions, a very precise Xray diffraction experiment at high-resolution (d 6 0.5 Å) and low temperature (T <100 K) has to be carried out to make these small aspherical effects visible. To meet these demanding experimental conditions the appearance of CCD area detectors in the mid-1990s and the availability of an increasing number of highly intense synchrotron beamlines of the third generation have contributed substantially. 9 Nevertheless, there still exist some limitations for a broad application of experimental ED work in the life sciences, where at least two major problems have to be addressed. (i) In medicinal chemistry, especially in drug development, it is common practise to screen over a large number of compounds to identify their potential activity. This makes also ED studies on series of molecules in short time periods necessary. (ii) Important biological molecules are macromolecules (proteins, polynucleotides), of which ED studies are known only in extremely rare cases. [10] [11] [12] [13] The above mentioned experimental advances have led to a dramatic reduction of X-ray exposure time to a few hours or even minutes per data collection so that in combination with methodical and computational developments the examination and quantitative analysis of series of ED distribution seems feasible already at present. 14 For macromolecules, a generally lower crystal quality, ubiquitous disorder, large amount of solvent in the crystal structure, cause altogether limited diffraction properties, so that these compounds are still an experimental challenge. As one approach to overcome these shortcomings, various activities have been initiated in the last few years to establish data bases of aspherical ED fragments obtained either from experimental or theoretical investigations. [15] [16] [17] One concept is to replace the IAM by the invariom formalism, 18 an aspherical scattering model which assigns fixed multipole populations from theoretical calculations (stored in an invariom library 19 ) to each atom of a structure to calculate nonspherical scattering factors. Following the transferability and additivity concept expressed in Bader's Qtaim theory, 20 the data base entries are then applied to the building blocks of a given macromolecule to generate additively the ED of the entire molecule. Making use in parts of this concept, we describe in this study an easy-to-use procedure to derive ED information at least in the active site regions of drug-receptor complexes from a combination of macromolecular X-ray data and theoretical calculations. Applications were considered 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 we analyzed the progesterone receptor complex. A procedure to obtain ED information for drug-receptor complexes was elaborated as follows (see also schematic representation in the flow chart of Fig. 1 ): (1) The X-ray structure atomic coordinates of the considered protein or protein complex were extracted from the protein data bank (PDB). 21 (2) The protein was reduced to an appropriate substructure in a region around the active site, hydrogens were added to free valences, their positions were set according to steric considerations using GaussView. 22 (3) The substructure obtained this way was entered into a single point ab initio density-functional calculation with the B3LYP/6-31G ** basis set, using GAUSSIAN03. 22 (4) The wave functions were entered into the program TONTO, 23 which generates a pseudo periodic structure, so that a theoretical crystallographic structure factor file (hkl-file) could be generated. (5) This hkl-file was used as input into the program XD, 24 which is normally used to refine an aspherical atomic electron density model against experimental structure factors. Aspherical scattering factors (multipole populations) from our invariom library are assigned to the atoms of the protein substructure and, if present, also to the substrate molecule. For atoms not in the invariom library multipoles are refined against the theoretical structure factors, while invariom multipoles, positional and thermal parameters are kept fixed. (6) The obtained electron density of the subfragment can then be analyzed by the corresponding subprograms of the XD program suite. We note that all substrate-in-protein structures are compared to optimized geometries of the corresponding isolated substrate molecules at the same level of theory. The procedure proposed above has a couple of shortcomings. One open question is the accuracy of the protein structure. Proteins in the PDB are frequently not determined up to atomic resolution. It is common practise to apply a restraint refinement with the aim to approach an idealized geometry of the contributing amino acids as introduced by Engh and Huber. 25 Since in our procedure atomic positions are not further refined, a potential inaccuracy remains. In the examples detailed below, the accuracy of bonds in terms of rms deviations from ideality was reported in the range of some hundreds of Å, which is about one order of magnitude larger than for small molecules, where bond length errors are generally smaller than 0.005 Å. A second point is the choice of the substructure in the active site region, that means whether it is appropriate or not. It has to be admitted, that the choice is somewhat ambiguous, however, care was taken to cover this region as complete as necessary on one hand but as small as possible on the other hand to allow affordable computing time in the subsequent calculations. The advantage of this procedure is that it is rather fast and can be applied with reasonable effort, because the considered molecular structure is much smaller than the entire protein or protein receptor complex and that due to the support of invariom parameters only moderate least-squares refinements are necessary. Moreover, the procedure provides an easy access to molecular surfaces and surfaces properties. Since recognition of interaction processes takes place via surface properties, their inspection is considered as a valuable source of information. Supported by corresponding graphical software developed in our lab (Moliso) 26 we concentrate in this study on two types of surfaces: (i) The electrostatic potential (ESP), calculated in the XDPROP routine of XD 24 after the method of Volkov et al. 27 mapped as colour code on an ED iso-surface of 0.0067 eÅ À3 (corresponding to 0.001 au). For the consideration of the reactive behaviour of a chemical system, the three-dimensional distribution of its ESP is very helpful in that negative regions can be regarded as nucleophilic centers, whereas regions with positive electrostatic potential are potential electrophilic sites. Moreover, the electrostatic potential makes the polarization of the electron density visible. (ii) The Hirshfeld surface, [28] [29] [30] which provides a computationally simple and elegant partitioning procedure for molecules in a crystal or another appropriate environment. It is defined by the ratio of the molecular density over the crystal density which should be a constant c, where c = 0.5 is normally chosen. We calculate the Hirshfeld surface according to this condition and map the non-spherical ED on this surface by a colour code. This allows the amount of density at a given site to be visualized. This way the Hirshfeld surface makes ED concentrations and hence sites and strengths of intermolecular interactions visible. The example in Figure 2 shows the experimentally obtained Hirshfeld surface of the recently examined indolylquinazoline derivative LA1810 (1), 31,32 a further potential EGFR tyrosine kinase inhibitor. This compound exists as a dimer in the crystal structure formed by a pair of N-HO hydrogen bonds. Further strong intermolecular contacts are not seen. Consequently the Hirshfeld surface shows only two strongly coloured regions, the N-H donor and the oxygen acceptor sites of the hydrogen bond. As the chemical information provided by the presented procedure is dependent on a proper construction of the three-dimensional ED distribution within the active site substructure of the proteins, it was necessary to evaluate the stability of the results by repeating the procedure for different substructures of the same active site. To meet the two questions raised before (accuracy of the protein structure, choice of the substructure), this quality test was applied to the progesterone receptor as it provides the most complete packing around the substrate molecule. On the one hand, a random generator was used to 'jiggle' the atomic positions. This was done in two steps. First, the amount of displacement in atomic coordinates was chosen to be in the range of the experimental uncertainties of the protein structure. In a second step, the atomic displacements were chosen significantly beyond the experimental errors. On the other hand, the substructure was reduced to half of its size, carrying then 20 amino acids in the protein pocket. For the jiggled structure the multipole populations of the original combined assignment/refinement were applied. For the new smaller substructure, the above presented procedure ( Fig. 1 ) was followed from its second step, which means, a DFT single point calculation was performed, followed by a calculation of theoretical structure factors and subsequent assignment/refinement of multipoles. For these modifications of the progesterone receptor substructures, reactive surfaces were generated and the topology of the intermolecular interactions was derived for the smaller substructure. The results are described in Section 3.4. and in parts in the Supplementary data. The EGFR is a transmembrane glycoprotein with an extracellular ligand-binding domain and an intracellular domain with tyrosine kinase activity for signal transduction. Ligand binding activates the receptor and its signalling pathways leading to the activation or modulation of cellular processes. [33] [34] [35] The receptor is expressed on healthy cells (40-100 EGFR/cell) as well as on malignant tissues (more than 1 000 000 EGFR/cell). 36 The epidermal growth factor (EGF) and the transforming growth factor-a (TGF-a) are the most important stimulatory ligands for EGFR. Overexpression of EGFR has been demonstrated in a wide variety of malignant cells. Thus, therapeutic strategies to inhibit EGFR and EGFR-related pathways have been pursued, including the development of ATP-competitive small molecule inhibitors of the intracellular tyrosine kinase domain of the receptor like gefitinib (Iressa Ò ) (2). The macromolecular structure of an EGFR tyrosine kinase domain complexed with the anilinoquinazoline derivative (2) is available at a resolution of 2.6 Å from the PDB, code 1M17. 37 We have extracted a subfragment consisting of 34 amino acids around the active site and also the structure of the anilinoquinazoline substrate. Making use of the procedure described above, the ESP of the isolated substrate (Fig. 3 above) , the ESP of the substrate in the substructure pocket (Fig. 3 below) and the ESP of the substructure pocket (Fig. 4) were generated. Whereas Figure 3 (above) shows only a moderate polarization of the isolated anilinoquinazoline gefitinib (2), the polarization increases considerably, if the substrate ESP is calculated in the environment of the pocket (Fig. 3 below). Moreover, we note that positively polarized regions dominate the substrate molecular surface. In accordance with these findings, the receptor pocket (Fig. 4) shows an ESP where negative regions dominate so that the principle of electrostatic complementarity is properly realized. Figure 5 shows the Hirshfeld surface of (2) where interactions with the receptor pocket are clearly visible at sites with strong colour code. In contrast to LA1810 (1) (Fig. 2 right) several sites of interaction are observed for the anilinoquinazoline gefitinib (2), indicating its increased potency for forming contacts with the protein (Scheme 1). SARS (severe acute respiratory syndrome) became apparent in fall 2002 and was distributed within a few weeks over China and other countries with the highest number of almost 4000 people infected in May 2003. 38, 39 SARS is caused by a previously unknown coronavirus, SARS-CoV. For the replication of this virus, its main cysteine protease M pro (molecular mass 33.8 kDa per protomer) is essential since it releases the key proteins in viral replication. As a therapeutic approach towards anti-SARS drugs, the inhibition of M pro by aziridine or epoxides was considered. The mechanism is displayed schematically in Figure 6 . A nucleophilic attack from the cysteinic sulphur of M pro takes place at a carbon atom of the epoxide inhibitor. This causes a ring opening and a covalent bond is formed between the cystein S c and the attacked epoxide carbon atom, so that an irreversible inhibition of M pro results. For the application of the procedure described above to the M pro inhibition procedure, we have made use of two corresponding macromolecular structures reported by Lee et al. 39 One is the unbound M pro structure (PDB code 2A5A, resolution 2.08 Å), the second one is a complex of M pro with an azapeptide epoxide (APE) substrate (R 1 = Cbz-Leu-Phe-AGln, R 2 = COOEt, (3), see Fig. 6 ), PDB code 2A5I (resolution 1.88 Å). After extraction of the active site subfragment (41 residues) we have obtained ESP's in these regions, which are shown in Figure 7 (top left) for the case of the unbound protein and in Figure 7 (top right) for the active site of the complex where the interaction with the substrate is taken into account. In addition, Figure 7 (bottom left and right) shows the ESP of the substrate before and after ring opening has taken place. While in the unbound active site pocket a more or less balanced distribution of positive and negative ESP sites is seen, the interior of the pocket is dominated by negative ESP in case of the complex which is in accordance with the mostly positive ESP surface of the substrate. Figure 8 shows the Hirshfeld surface of the APE substrate with the electron density mapped as colour code on this surface. As mentioned before, deeply coloured regions show high electron density concentrations being generally indicative for strong intra or intermolecular interactions. We note that the strongest ED is found just at the site where the covalent bond with the S c has been established. The steroid hormone progesterone is a key component in the complex regulation of normal female reproductive function. The major physiological roles of progesterone in the mammal are (1) in the uterus and ovary: release of mature oocytes, facilitation of implantation, and maintenance of pregnancy, by promotion of uterine growth and suppression of myometrial contractility; (2) in the mammary gland: lobular-alveolar development in preparation for milk secretion and suppression of milk protein synthesis before parturition; and (3) in the brain: mediation of signals required for sexually responsive behaviour. Recent evidence also supports a role for progesterone in modulation of bone mass. Progesterone effects are mediated by its nuclear receptor, which is a dimer composed of two receptor progesterone binding proteins, PR A and PR B. [40] [41] [42] The crystal structure of a progesterone-bound ligand-binding domain of the human progesterone receptor was published at 1.8 Å resolution (PDB code 1A28) by Williams and Sigler. 43 The interactions of the hormone with the ligand-binding domain (LBD) were discussed in details based on steric considerations. Strong interactions were identified from the protein to the 3-keto group and less pronounced to the methyl-ketone substituent at the C17 carbon of the progesterone (Scheme 2). We have generated a 40 residue subfragment around the hormone binding pocket and calculated ESP's for progesterone in the isolated state and bound to the pocket (Fig. 9) . In both cases the hydrocarbon skeleton region exhibits a more or less uniform posi- tive ESP, whereas the 3-keto region and the methyl-ketone substituent are negatively polarized, being less pronounced in the methyl-ketone region. Hence, the interacting protein region should provide the opposite charge at the corresponding sites. Figure 10 shows a selected region of ESP in the hormone binding pocket of the LBD, where a rather strong positive region is next to the 3-keto oxygen, while otherwise mostly negative ESP faces parts of the positive region of the hydrocarbon section of the hormone. In Figure 11 a further Hirshfeld surface representation is given, showing the hormone (4) in the same orientation as in Figure 9 (lower part) in its pocket of the LBD. It can easily be recognized that the most coloured region and hence the strongest interaction is at the 3-keto group followed by the oxygen of the methyl-ketone. However, further less intense sites of interactions are visible, so that a close inspection of this surface can help to identify additional interactions of relevance. The strength of the electron density at a considered site can serve as a quantitative measure for the interaction in question. For evaluation of the stability of the results, three modifications of the progesterone receptor substructure were created and analyzed. In one approach, using a random generator, the original coordinates were displaced ('jiggled') randomly in two steps. In a first step an average displacement of atomic coordinates of 0.05 Å was applied which is slightly above the reported experimental uncertainty of this progesterone receptor complex of $0.01 Å. 43 In a second step, the random displacement was enlarged to 0.5 Å, far beyond the experimental error. For both models the ESP's were generated for the substrate molecule and the protein pocket. as ball-and-stick representation. A water molecule which also interacts with the hormone is displayed in transparent mode. Figure 11 . Hirshfeld representation of the hormone (4) in its pocket of the LBD. For the 0.05 Å randomization the ESP's are indistinguishable to the results for the original substructure, see Figure 10 and upper part of Figure 12 (see also Figs. S1 and S2, Supplementary data). This does no longer hold for the stronger 0.5 Å randomization (Fig. 12 below and Fig. S3 ). Especially the ESP features in the pocket displayed in Figure 12 (below) differ strongly from the original substructure. This finding can be supported quantitatively making use of the parameters of the Politzer analysis 44 which examines the distribution of positive/negative ESP values on a given surface. Table 1 shows the Politzer parameter p, which gives the average deviation from the overall average potential. 45 For both the substrate and the pocket, the p values are identical for the original and the jiggle 0.05 model but change significantly when going to the jiggle 0.5 model. Hence, if the uncertainty of the model is close to or slightly above the experimental error, the procedure described here is unaffected but becomes questionable if the randomization is extended to the large value of 0.5 Å. Alternatively, a smaller substructure of 20 amino acid residues was generated (see Fig. S4 ) for which density-functional calculation, calculation of theoretical structure factors, and subsequent multipole assignment/refinement was performed. Although the results for the ESP's are somewhat different (Fig. 13) , which is expected for a different substructure, the overall agreement with the major features of the original 40 residue fragment is good. Especially the extended negative region neighbouring the positive hydrocarbon section of progesterone is clearly visible (compare Fig. 10 ). The positive region next to the 3-keto oxygen is less pronounced. This is confirmed by the Hirshfeld surface (Fig. 14) where sites of intermolecular interactions are in most cases comparably distributed as for the 40 residue substructure (see Fig. 11 ). The more qualitative discussion on intermolecular interactions between substrates and active sites based on reactive surfaces is complemented by a topological analysis of the ED at the bond critical points (bcp) found between hydrogen bonded atoms. 20 Espinosa and co-workers investigated the dependence of various topological parameters, [46] [47] [48] such as the electron density q(r bcp ), the Laplacian r 2 q(r bcp ) of the ED, the Hessian Eigenvalues k 1 -k 3 of the Laplacian, and the kinetic and potential energy densities, G(r bcp ) and V(r bcp ), respectively, on different geometrical quantities for a characterization of the strength and nature of intermolecular interactions. Amongst others, they found systematical relations between the HÁÁÁO distances and q(r bcp ), r 2 q(r bcp ), and k 3 , respectively, for an analysis of 83 X-HÁÁÁO (X = C, N, O) contacts. 47 The Figure 13 . Above: ESP of progesterone in the 20 residue substructure. There are small differences to the representations in the lower part of Figure 9 . Below: ESP in a smaller section of the 20 residue substructure close to the progesterone substrate. functions obtained by this statistical approach were applied to the 60 hydrogen bonds found in this work for all three protein substructures including the substrate molecules. These hydrogen bonds appear between substrate and protein as well as between different amino acids of the protein site. Figure 15 displays the values for q(r bcp ), r 2 q(r bcp ), and k 3 in relation to the HÁÁÁY (Y = O, N, S) distances. Although the X-H (X = C, N, O, S) distances of the protein structures are not fixed to values obtained by neutron diffraction, an approach which is commonly used in ED studies, the overall agreement between the functions taken from Ref. 47 (details see Fig. 15 ) and the topological properties obtained by the above described procedure is very good. Especially k 3 matches nearly perfectly the expected progression. On one hand, these results confirm the quality of the ED distribution of the whole substructure, on the other hand this procedure enables the detailed analysis of the hydrogen bonding network within these large systems. Hydrogen bonds are one source of weak intra and intermolecular interactions and play an important role in molecular association and in ligand-target interactions in which hydrogen bond donors face hydrogen bond acceptors; hence the topological analysis described above provides a tool for a quantitative study of these important contacts. A corresponding hydrogen bond analysis for the reduced 20 residue progesterone substructure, introduced in Section 3.4 was also carried out. For results, see 'Supplementary data', Figure S5 . For biological processes the understanding of intermolecular recognition has to be understood on an atomic scale. Considering the fact that in a first step the mutual recognition takes place via the surfaces of the involved molecules, we have developed a procedure to generate and display molecular surfaces and selected ED surface properties of biomolecules. The present limitations to obtain ED information for macromolecules are diminished by reducing the macromolecular protein structure to a region around the active site, so that the computational effort becomes affordable. The use of two data bases is needed. Experimental X-ray structural information is provided by the protein data base (PDB), while atomic electron density data are supplied by the Invariom data base which allows then the additive generation of the ED of the entire system following Bader's transferability and additivity concept. 20 The procedure is easy to apply with the presently existing software and provides surface properties that allow the identification of sites and strengths of intermolecular interactions. This was shown for the ESP surfaces in protein active site pockets and corresponding substrates and on Hirshfeld surfaces, on which interactions directly become visible in terms of ED concentrations. We believe that the joint information from the ESP (to indicate nucleophilic/electrophilic sites) and from the Hirshfeld surface (to identify sites and strengths of preferred interactions), which in this combination can only be obtained from the ED, provides various options to enter electronic information into the examination of drug-receptor interactions with acceptable effort. If, for example, docking or fitting techniques for substrate-target systems are carried out based on steric considerations, the obtained geometry can then be entered into the procedure described here to examine also the electronic complementarity. In addition to the possibilities to determine drug-receptor interactions, also potential sites at the receptor can be identified which can preferably be attacked by drugs. These findings can be helpful for the design of chemical structures as new agonists/antagonists on a specific receptor. One shortcoming should not be overlooked. It is the choice and accuracy of the medium molecular weighted substructure which should mimic the active site region. However, the considered examples have shown that even the quantitative analysis of hydrogen bonding in the chosen substructure systems follows the progression developed by Espinosa et al. 47 on the basis of accurate small molecule crystal structures. Furthermore, the analysis of the additional model systems introduced for the progesterone receptor active site (see Section 3.4), the jiggled structure on the one hand and the smaller pocket on the other hand, reveal that the results presented in this study are stable to a high degree. The ESP is virtually identical for the original and the jiggle 0.05 Å pockets and only slightly different for the smaller 20 residue model. The latter one also shows minor changes in the Hirshfeld surface as well as in the topology of the hydrogen bonds (see Supplementary data). So we are confident that as long as fine details of the atomic structure are not considered, the surface properties are robust enough to provide reliable hints for the recognition and interaction of biomolecular species. Atoms in Molecules-A Quantum Theory Protein Data Bank A Fortran Based Object-Oriented System for Quantum Chemistry and Crystallography, User Manual 39, 901. Software available from The authors thank the Deutsche Forschungsgemeinschaft (DFG) for financial support within the special priority program SPP1178, grant Lu 222/29-2. Supplementary data (additional description of the modified substructures of the progesterone receptor complex and illustration of further results in Figures S1-S5) associated with this article can be found, in the online version, at doi:10.1016/j.bmc.2010. 06.080.