A molecular modeling study on dihydrofolate reductase (DHFR) inhibitors was carried out. By combining molecular dynamics simulations with semiempirical (PM6), ab initio, and density functional theory (DFT) calculations, a simple and generally applicable procedure to evaluate the binding energies of DHFR inhibitors interacting with the human enzyme is reported here, providing a clear picture of the binding interactions of these ligands from both structural and energetic viewpoints. A reduced model for the binding pocket was used. This approach allows us to perform more accurate quantum mechanical calculations as well as to obtain a detailed electronic analysis using the quantum theory of atoms in molecules (QTAIM) technique. Thus, molecular aspects of the binding interactions between inhibitors and the DHFR are discussed in detail. A significant correlation between binding energies obtained from DFT calculations and experimental IC₅₀ values was obtained, predicting with an acceptable qualitative accuracy the potential inhibitor effect of nonsynthesized compounds. Such correlation was experimentally corroborated synthesizing and testing two new inhibitors reported in this paper.

We report here an exhaustive and complete conformational study on the conformational potential energy hypersurface (PEHS) of dopamine (DA) interacting with the dopamine D2 receptor (D2-DR). A reduced 3D model for the binding pocket of the human D2-DR was constructed on the basis of the theoretical model structure of bacteriorhodopsin. In our reduced model system, only 13 amino acids were included to perform the quantum mechanics calculations. To obtain the different complexes of DA/D2-DR, we combined semiempirical (PM6), DFT (B3LYP/6-31G(d)), and QTAIM calculations. The molecular flexibility of DA interacting with the D2-DR was evaluated from potential energy surfaces and potential energy curves. A comparative study between the molecular flexibility of DA in the gas phase and at D2-DR was carried out. In addition, several molecular dynamics simulations were carried out to evaluate the molecular flexibility of the different complexes obtained. Our results allow us to postulate the complexes of type A as the "biologically relevant conformations" of DA. In addition, the theoretical calculations reported here suggested that a mechanistic stepwise process takes place for DA in which the protonated nitrogen group (in any conformation) acts as the anchoring portion, and this process is followed by a rapid rearrangement of the conformation allowing the interaction of the catecholic OH groups.

No abstract

In this paper, we reported the results obtained by charge density analysis of the network of non‐covalent interactions (NCI) established in the binding pocket of a receptor, in relevant conformations of ligand – receptor complexes. Starting with strong and moderate hydrogen bonds, moving on to weaker polar interactions and ending with stacking and T‐shape like interactions between aromatic rings, all of them have been investigated within the framework of the density functional theory and the quantum theory of atoms in molecules. Also, natural bond orbital analysis was carried out, in order to evaluate quantitatively the electronic population of the aromatic rings. The analysis of our “case study” shows that the interactions of the catechol OH groups of the ligand, in the different conformations of the dopamine (DA) – D2 receptor complex, determine the decrease or increase of the electron density on the aromatic ring of DA. In turn, the electronic population of the aromatic ring of DA defines its orientation within the binding site and the type of interactions that it establishes with the aromatic rings of the receptor. Although the approach used here was traditionally applied to the study of NCI in small molecules complexes in gas phase, we show through this work that this methodology is also a very powerful tool for the study of biomolecular complexes, providing a very detailed description of the binding event. Copyright © 2013 John Wiley & Sons, Ltd.

No abstract

We report an exhaustive conformational and electronic study on dopamine (DA) interacting with the D2 dopamine receptor (D2DR). For the first time, the complete surface of the conformational potential energy of the complex DA/D2DR is reported. Such a surface was obtained through the use of QM/MM calculations. A detailed study of the molecular interactions that stabilize and destabilize the different molecular complexes was carried out using two techniques: Quantum Theory of Atoms in Molecules computations and nuclear magnetic shielding constants calculations. A comparative study of the behavior of DA in the gas phase, aqueous solution, and in the active site of D2DR has allowed us to evaluate the degree of deformation suffered by the ligand and, therefore, analyze how rustic are the lock‐key model and the induced fit theory in this case. Our results allow us to propose one of the conformations obtained as the “biologically relevant” conformation of DA when it is interacting with the D2DR.

The multidimensional conformational potential energy hypersurface (PEHS) of cyclic Gly-Gly-Gly (1,4,7-triazonane-2,5,8-trione) was comprehensively investigated at the Hartree-Fock (RHF/6-31G(d)) level of theory. The equilibrium structures, their relative stability, and the transition state (TS) structures involved in the conformational interconversion pathways were analyzed. aug-cc-pVTZ//B3LYP/6-311++G** single point calculations predict a trans-cis-cis conformation as the energetically preferred form for this compound. However, all of the levels of theory employed here predicted that two forms, a trans-cis-cis and a cis-cis-cis (crown), of conformers contribute significantly to the equilibrium mixture at room temperature. The conformational interconversion between the global minimum and the symmetric cis-cis-cis crown form requires 12.49 kcal/mol at the RHF 6-31G(d) level of theory, whereas the conformational interconversion between the cis-cis-cis crown and cis-cis-cis boat form requires 18.70 kcal/mol. An exploratory topological analysis of the PEHS was also carried out. Our results allow us to form a concise idea about the internal intricacies of the PEHSs of these cyclic tripeptides, describing the conformations as well as the conformational interconversion processes in these hypersurfaces.

In this study, we described quantitatively the interactions between two new amino‐2H‐imidazole inhibitors ((R)‐1t and (S)‐1m) and BACE1 using a hybrid quantum mechanics‐molecular mechanical (QM/MM) method together with a quantum theory of atoms In Molecules (QTAIM) analysis. Our computational calculations revealed that the binding affinity of these compounds is mostly related to the amino‐2H‐imidazole core, which interact tightly with the aspartate dyad of the active site. The interactions were stronger when the inhibitors presented a bulky substituent with a hydrogen bond acceptor motif pointing toward Trp76, such as the 3,5‐dimethyl‐4‐methoxyphenyl group of compound (S)‐1m. Furthermore, the QTAIM analysis revealed that many hydrophobic interactions complement cooperatively the hydrogen bond which is not present when compound (R)‐1t is bound to the enzyme. The combined QM/MM‐QTAIM analysis allows identifying the interactions that account for the activity difference between compounds, even at a nanomolar range.