QM/MM model study on properties and structure of some antibiotics in gas phase: Comparison of energy and NMR chemical shift

The combination of Quantum Mechanics (QM) and Molecular Mechanics (MM) methods has become an alternative tool for many applications for which pure QM and MM are not suitable. The QM/MM method has been used for different types of problems, for example: structural biology, surface phenomena, and liquid phase. In this paper, we have used these methods for antibiotics and then we compared results. The calculations were done by the full ab initio method (HF/3-21G) and the (HF/STO-3G) and QM/MM (ONIOM) method with HF (3-21G)/AM1/UFF and HF (STO-3G)/AM1/UFF. We found the geometry obtained by the QM/MM method to be very accurate, and we can use this rapid method in place of time consuming ab initio methods for large molecules. Comparison of energy values in the QM/MM and QM methods is given. In the present work, we compare chemical shifts and conclude that the QM/MM method is a perturbed full QM method. The work has been done on penicillin, streptomycin, benzyl penicillin, neomycin, kanamycin, gentamicin, and amoxicillin.     

A computational study of a host-guest complex

The geometry and energetics of a complex involving pyrazine and an acridine diacid cleft-like host designed by Rebek were investigated at several levels of theory. Molecular mechanics (using the Tripos and CHARMm force fields), semiempirical quantum chemical approaches (with the AM1 and PM3 methods), and an ab initio quantum chemical method (RHF/STO-3G) were used in the complete relaxation of the complex. The geometry of the complex optimized by the RHF/STO-3G method is in excellent agreement with a published X-ray structure; upon superposition, the rms deviation between the corresponding cleft heavy atoms is only 0.17 Å and the pyrazine molecules are superimposable. In addition, ab initio quantum chemical techniques were used to study the complex when the cleft is modeled by a pair of acetic acid molecules. All the calculations presented herein support a two-point interaction mechanism. The similarities found in the results for the full complex and the truncated model are consistent with a purely structural role for the acridine linker of the host. © 1997 John Wiley & Sons, Ltd.     

QM/MM study of the active site of free papain and of the NMA-papain complex.

Hybrid quantum mechanical/molecular mechanical (QM/MM) calculations using restricted and unrestricted Hartree-Fock and B3LYP ab initio (QM) and Amber force field (MM), respectively, have been applied to study the catalytic site of papain in both free and substrate bonded forms. Ab initio geometry optimizations have been performed for the active site of papain and the N-methyl-acetamide (NMA)-papain complex within the molecular mechanical treatment of the protein environment. A covalent tetrahedral intermediate structure could be obtained only when the amide N atom of the substrate molecule was protonated through a proton transfer from the His-159 in the catalytic site. Our results support the previous assumption that a proton transfer from His-159 to the amide N atom of the substrate occurs prior to or concerted with the nucleophilic attack of the Cys-25 sulfur atom to the carbonyl group of the substrate. The electron correlation effect will reduce the proton transfer barrier. Therefore, this proton transfer can be easily observed in the B3LYP/6-31G* calculations. The HF/6-31G* method overestimates the reaction barrier against this proton transfer. The sulfur atom of Cys-25 and the imidazole ring of His-159 are found to be coplanar in the free form of the enzyme. However, the rotation of the imidazole ring of His-159 was observed during the formation of the tetrahedral intermediate. Without the papain environment, the coplanar thiolate-imidazolium ion pair RS-...ImH+ is much less stable than the neutral form of RSH....Im. Within the protein environment, however, the thiolate-imidazolium ion pair becomes more stable than its neutral form by 4.1 and 0.4 kcal/mol in HF/6-31G* and B3LYP/6-31G* calculations, respectively. The barrier of proton transfer from S-H group of Cys-25 to the imidazole ring of His-159 was reduced from 22.0 kcal/mol to 15.2 kcal/mol by the protein environment in HF/6-31G* calculations. This barrier is found to be much smaller (2.5 kcal/mol) in B3LYP/6-31G* calculations.     

Determinants of reactivity and selectivity in soluble epoxide hydrolase from quantum mechanics/molecular mechanics modeling

Soluble epoxide hydrolase (sEH) is an enzyme involved in drug metabolism that catalyzes the hydrolysis of epoxides to form their corresponding diols. sEH has a broad substrate range and shows high regio- and enantioselectivity for nucleophilic ring opening by Asp333. Epoxide hydrolases therefore have potential synthetic applications. We have used combined quantum mechanics/molecular mechanics (QM/MM) umbrella sampling molecular dynamics (MD) simulations (at the AM1/CHARMM22 level) and high-level ab initio (SCS-MP2) QM/MM calculations to analyze the reactions, and determinants of selectivity, for two substrates: trans-stilbene oxide (t-SO) and trans-diphenylpropene oxide (t-DPPO). The calculated free energy barriers from the QM/MM (AM1/CHARMM22) umbrella sampling MD simulations show a lower barrier for phenyl attack in t-DPPO, compared with that for benzylic attack, in agreement with experiment. Activation barriers in agreement with experimental rate constants are obtained only with the highest level of QM theory (SCS-MP2) used. Our results show that the selectivity of the ring-opening reaction is influenced by several factors, including proximity to the nucleophile, electronic stabilization of the transition state, and hydrogen bonding to two active site tyrosine residues. The protonation state of His523 during nucleophilic attack has also been investigated, and our results show that the protonated form is most consistent with experimental findings. The work presented here illustrates how determinants of selectivity can be identified from QM/MM simulations. These insights may also provide useful information for the design of novel catalysts for use in the synthesis of enantiopure compounds.     

A QM/MM analysis of the conformations of crystalline sucrose moieties

Both ab initio quantum mechanics (QM) and molecular mechanics (MM) were used to produce a hybrid energy surface for sucrose that simultaneously provides low energies for conformations that are observed in crystal structures and high energies for most unobserved structures. HF/6-31G* QM energies were calculated for an analogue based on tetrahydropyran (THP) and tetrahydrofuran (THF). Remaining contributions to the potential energy of sucrose were calculated with MM. To do this, the MM surface for the analogue was subtracted from the MM surface for the disaccharide, and the QM surface for the analogue was added. Prediction of the distribution of observable geometries was enhanced by reducing the strength of the hydrogen bonding. Reduced hydrogen-bonding strength is probably useful because many crystalline sucrose moieties do not have intramolecular hydrogen bonds between the fructose and glucose residues. Therefore, hydrogen bonding does not play a large role in determining the molecular conformation. On the hybrid energy surface that was constructed with a dielectric constant of 3.5, the average potential energy of 23 sucrose moieties from crystal structures is 1.16 kcal/mol, and the population of observed structures drops off exponentially as the energy increases.     

Study of nonplanarity of peptide bond using theoretical calculations

The conformational dependence of nonplanarity of the peptide bond of formylglycinamide has been studied using ab initio and density functional theory methods. Hartree-Fock self-consistent field theory (HF), M?ller-Plesset perturbation theory (MP2) of ab initio and B3LYP level of theory of dft method have been used employing 6-31++G** basis set. The MP2 method predicts better results than HF and B3LYP levels of theory for conformational stability dependence of nonplanarity. Systematic dependence of planarity deviation has been observed in MP2 theory. The chemical hardness values successfully predict the conformational region, but fail to obey maximum hardness principle. It is concluded that the most reliable dft method could not successfully predict the planarity of peptide bond in comparison with electron correlated method of ab initio method.     

Theoretical study of antagonists and inhibitors of mammalian adenosine deaminase. II. Isomeric aza-deazaanalogues of adenosine

Isomeric aza-deazaanalogues of adenosine and their N1-protonated forms (except for that of 8-aza-1-deazaadenosine) were studied by computer modeling to find a relationship between their molecular structures and the properties as substrates for the mammalian adenosine deaminase. The atomic charge distribution and maps of the electrostatic potential around their van der Waals molecular surface were calculated using the ab initio STO-3G method. The conformational studies were carried out by the MM+ method of molecular mechanics. The previously proposed mechanism of the substrate acceptance in the active site of mammalian adenosine deaminase was refined, and the potential substrate properties were predicted for two previously unstudied adenosine analogues, 5-aza-9-deazaadenosine and 8-aza-3-deazaadenosine.     

Theoretical investigation of the triphosphate forms of azidothymidine and thymidine

In this paper we investigate (using AM1 semi-empirical as well as HF methods at the STO-3G, 3-21G, 6-31G, 6-31G* and 6-31+G** level) the conformations, geometrical parameters, Mulliken charges, and solvation effects of the triphosphate form of AZT (AZTTP), as well as the thymidine nucleotide (dTTP) structure. Our calculated geometrical parameters and Mulliken charges, with and without solvation effects, are correlated with recent experimental results.     

An AB initio investigation of molecules with a disulfide bond: (HS)2, (CH3S)2 and cystine

ab initio Calculations at the Gaussian-70 STO-3G and 4-31G basis levels have been carried out for (HS)2 and (CH3S)2. Cystine was investigated at the STO-3G level. The STO-3G energy minimized geometry agrees well with experiments for (HS)2 and (CH3S)2. The barriers to internal rotation are predicted to be (at the 4.31G level): (HS)2, cis 8.5 kcal, trans 3.03 kcal; (CH3S)2, cis 18.47 kcal, trans 6.04 kcal.     

Theoretical study of antagonists and inhibitors of mammalian adenosine deaminase. I. Adenosine and its aza- and deaza-analogs

Aza- and deazaanalogues of adenosine, including their 1-protonated forms (except for that of 1-deazaadenosine), were studied by computer computation to find a relationship between their molecular structures and substrate properties for the mammalian adenosine deaminase. The atomic charge distribution and maps of the electrostatic potential around their van der Waals molecular surface were calculated for these compounds using the ab initio STO-3G method. The conformational studies were carried out by the MM+ method of molecular mechanics. The mechanism that determines the substrate selectivity of mammalian adenosine deaminase is discussed. The English version of the paper: Russian Journal of Bioorganic Chemistry, 2002, vol. 28, no. 4; see also http://www.maik.ru.     

The Hydrolysis Mechanism of Formamide Revisited: Comparison Between ab initio, Semiempirical and DFT Results

The hydrolysis of amides is a model reaction to study peptide hydrolysis. This process has been previously considered in the literature at the ab initio level. In this work, we revisit different reaction mechanisms (water-assisted, non-assisted, neutral and acid-catalyzed) with various theoretical methods : semiempirical, ab initio and Density Functional. The ab initio calculations are carried out at a computational level which is substantially higher than in previous studies. We describe the structure of the transition states and discuss the influence of the catalyst. We also compute the activation free energies for these processes at the Density Functional Theory level. Comparison of the methods allows to outline the main trends of these theoretical approaches which may be useful to design new computational strategies for investigating biological reaction mechanisms through the use of combined Quantum Mechanics/Molecular Mechanics methods.     

Lanthanide complex coordination polyhedron geometry prediction accuracies of ab initio effective core potential calculations

The ability to efficiently and accurately predict solid-state geometries of lanthanide coordination compounds efficiently and accurately is central for the design of new ligands capable of forming stable and highly luminescent complexes. Accordingly, we present in this paper a report on the capability of various ab initio effective core potential calculations in reproducing the coordination polyhedron geometries of lanthanide complexes. Starting with all combinations of HF, B3LYP and MP2(Full) with STO-3G, 3-21G, 6-31G, 6-31G* and 6-31+G basis sets for [Eu(H₂O)₉]³⁺ and closing with more manageable calculations for the larger complexes, we computed the fully predicted ab initio geometries for a total of 80 calculations on 52 complexes of Sm(III), Eu(III), Gd(III), Tb(III), Dy(III), Ho(III), Er(III) and Tm(III), the largest containing 164 atoms. Our results indicate that RHF/STO-3G/ECP appears to be the most efficient model chemistry in terms of coordination polyhedron crystallographic geometry predictions from isolated lanthanide complex ion calculations. Moreover, both augmenting the basis set and/or including electron correlation generally enlarged the deviations and aggravated the quality of the predicted coordination polyhedron crystallographic geometry. Our results further indicate that Cosentino et al.’s suggestion of using RHF/3-21G/ECP geometries appears to be indeed a more robust, but not necessarily, more accurate recommendation to be adopted for the general lanthanide complex case. Figure Graphical visualization of unsigned mean errors, UME(Eu-L)s, involving only the interatomic distances between the europium central ion and the oxygen atoms of the coordination polyhedron of the cation nona-aqua-europium(III) for various model chemistries, all compared to the “Cambridge Structural Database 2004” crystallographic geometry     

Azaguanine: A Theoretical Study of Its Tautomerism and Protonation in the Gas Phase and Aqueous Solution

The tautomerism and protonation of 8-azaguanine (azaG) have been studied by means of ab initio methods, both in the gas phase and in aqueous solution. An elimination procedure to choose the most stable tautomeric forms, based on AM1 and HF/6–31G* energies, has been applied. Tautomers azaG(1,9), azaG(1,7) and azaG(9,15) have been selected and their energies calculated at MP2/6–311++G**//HF/6–31G* level. Self-consistent reaction field IPCM, based on polarizable continuum model (PCM), has been applied to study the solvent effects. The stability order in the gas phase is azaG(1,7) = azag(1,9) > azaG(9,15), whereas in solution the order becomes azaG(1,7) > azaG(1,9) > azaG(9,15), the latter being just 0.76 kcal/mol over azaG(1,7). The calculations of proton affinities allowed to unambiguously determine the preferred sites of protonation of these species.     

Computational study of Au_4 cluster on a carbon nanotube with and without defects using QM/MM methodology

We use ONIOM (QM/MM) methodology to carry out geometry calculations in a 4-atom nanocluster supported by an (8, 8) armchair carbon nanotube with and without defects employing LSDA/SDD for the QM system and UFF for MM. In two particular cases, defects were added in the carbon nanotube wall. These defects are a double oxygenated vacancy (Vac(2)O(2)) and a double vacancy but without oxygen which creates two pentagons and an octagon. Our results show how geometries using QM/MM and energies calculations carried out with QM, change on both the gold nanocluster and the carbon nanotube. In addition, an application of ONIOM methodology in a comparative study to predict behavior of structures as hybrid materials based in carbon nanotubes combined with gold nanoclusters is shown. In this work we examine geometry changes on both the gold nanocluster and the carbon nanotube. A comparison is made with the binding energy resulting values as well as with the orbital energies such as the frontier orbitals HOMO and LUMO.     

Modelling flavin and substrate substituent effects on the activation barrier and rate of oxygen transfer by p-hydroxybenzoate hydroxylase

The simulation of enzymatic reactions, using computer models, is becoming a powerful tool in the most fundamental challenge in biochemistry: to relate the catalytic activity of enzymes to their structure. In the present study, various computed parameters were correlated with the natural logarithm of experimental rate constants for the hydroxylation of various substrate derivatives catalysed by wild-type para-hydroxybenzoate hydroxylase (PHBH) as well as for the hydroxylation of the native substrate (p-hydroxybenzoate) by PHBH reconstituted with a series of 8-substituted flavins. The following relative parameters have been calculated and tested: (a) energy barriers from combined quantum mechanical/molecular mechanical (QM/MM) (AM1/CHARMM) reaction pathway calculations, (b) gas-phase reaction enthalpies (AM1) and (c) differences between the HOMO and LUMO energies of the isolated substrate and cofactor molecules (AM1 and B3LYP/6-31+G(d)). The gas-phase approaches yielded good correlations, as long as similarly charged species are involved. The QM/MM approach resulted in a good correlation, even including differently charged species. This indicates that the QM/MM model accounts quite well for the solvation effects of the active site surroundings, which vary for differently charged species. The correlations obtained demonstrate quantitative structure activity relationships for an enzyme-catalysed reaction including, for the first time, substitutions on both substrate and cofactor.     

Multipole electrostatic potential derived atomic charges in NDDO-methods with <Emphasis Type="Italic">spd</Emphasis>-basis sets

The recently introduced multipole approach for computing the molecular electrostatic potential (MEP) within the semiempirical neglect of diatomic differential overlap (NDDO) framework [Horn AHC, Lin Jr-H., Clark T (2005) Theor Chem Acc 114:159–168] has been used to obtain atomic charges of nearly ab initio quality by scaling the semiempirical MEP. The parameterization set comprised a total of 797 compounds and included not only the newly parameterized AM1* elements Al, Si, P, S, Cl, Ti, Zr, and Mo but also the standard AM1 elements H, C, N, O and F. For comparison, the ZDO-approximated MEP was also calculated analytically in the spd-basis. For the AM1*-optimized structures, single-point calculations at the B3LYP, HF and MP2 levels with the 6-31G(d) and LanL2DZP basis sets were performed to obtain the MEP. The regression analysis of all 12 combinations of semiempirical and ab initio MEP data yielded correlation coefficients of at least 0.99 in all cases. Scaling the analytical and multipole-derived semiempirical MEP by the regression coefficients yielded mean unsigned errors below 2.6 and 1.9 kcal mol⁻¹, respectively. Subsequently, for 22 drug molecules from the World Drug Index, atomic charges were computed according to the RESP procedure using XX/6-31G(d) (XX=B3LYP, HF, MP2) and scaled AM1* multipole MEP; the correlation coefficients obtained are 0.83, 0.85 and 0.83, respectively. Figure: Schematic representation of the atomic charge generation: The molecular electrostatic potential (MEP) is calculated using the AM1* Hamiltonian; then the semiempirical MEP is scaled to DFT or ab initio level, and atomic charges are generated subsequently by the restraint electrostatic potential (RESP) fit method. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible to authorized users. Proceedings of “Modeling Interactions in Biomolecules II”, Prague, September 5th–9th, 2005.     

A comparative study of trypsin specificity based on QM/MM molecular dynamics simulation and QM/MM GBSA calculation

Hydrogen bonding and polar interactions play a key role in identification of protein-inhibitor binding specificity. Quantum mechanics/molecular mechanics molecular dynamics (QM/MM MD) simulations combined with DFT and semi-empirical Hamiltonian (AM1d, RM1, PM3, and PM6) methods were performed to study the hydrogen bonding and polar interactions of two inhibitors BEN and BEN1 with trypsin. The results show that the accuracy of treating the hydrogen bonding and polar interactions using QM/MM MD simulation of PM6 can reach the one obtained by the DFT QM/MM MD simulation. Quantum mechanics/molecular mechanics generalized Born surface area (QM/MM-GBSA) method was applied to calculate binding affinities of inhibitors to trypsin and the results suggest that the accuracy of binding affinity prediction can be significantly affected by the accurate treatment of the hydrogen bonding and polar interactions. In addition, the calculated results also reveal the binding specificity of trypsin: (1) the amidinium groups of two inhibitors generate favorable salt bridge interaction with Asp189 and form hydrogen bonding interactions with Ser190 and Gly214, (2) the phenyl of inhibitors can produce favorable van der Waals interactions with the residues His58, Cys191, Gln192, Trp211, Gly212, and Cys215. This systematic and comparative study can provide guidance for the choice of QM/MM MD methods and the designs of new potent inhibitors targeting trypsin.