QM/MM calculations
Dotaz
Zobrazit nápovědu
The density functional code deMon2k employs a fitted density throughout (Auxiliary Density Functional Theory), which offers a great speed advantage without sacrificing necessary accuracy. Powerful Quantum Mechanical/Molecular Mechanical (QM/MM) approaches are reviewed. Following an overview of the basic features of deMon2k that make it efficient while retaining accuracy, three QM/MM implementations are compared and contrasted. In the first, deMon2k is interfaced with the CHARMM MM code (CHARMM-deMon2k); in the second MM is coded directly within the deMon2k software; and in the third the Chemistry in Ruby (Cuby) wrapper is used to drive the calculations. Cuby is also used in the context of constrained-DFT/MM calculations. Each of these implementations is described briefly; pros and cons are discussed and a few recent applications are described briefly. Applications include solvated ions and biomolecules, polyglutamine peptides important in polyQ neurodegenerative diseases, copper monooxygenases and ultra-rapid electron transfer in cryptochromes.
- MeSH
- kvantová teorie MeSH
- lidé MeSH
- molekulární modely MeSH
- peptidy chemie MeSH
- simulace molekulární dynamiky MeSH
- software * MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
- přehledy MeSH
Hybrid QM/MM methods combine the rigor of quantum mechanical (QM) calculations with the low computational cost of empirical molecular mechanical (MM) treatment allowing to capture dynamic properties to probe critical atomistic details of enzyme reactions. Catalysis by RNA enzymes (ribozymes) has only recently begun to be addressed with QM/MM approaches and is thus still a field under development. This review surveys methodology as well as recent advances in QM/MM applications to RNA mechanisms, including those of the HDV, hairpin, and hammerhead ribozymes, as well as the ribosome. We compare and correlate QM/MM results with those from QM and/or molecular dynamics (MD) simulations, and discuss scope and limitations with a critical eye on current shortcomings in available methodologies and computer resources. We thus hope to foster mutual appreciation and facilitate collaboration between experimentalists and theorists to jointly advance our understanding of RNA catalysis at an atomistic level.
- MeSH
- biofyzika metody MeSH
- fosfáty chemie MeSH
- fosforylace MeSH
- hořčík chemie MeSH
- katalýza MeSH
- konformace nukleové kyseliny MeSH
- kvantová teorie MeSH
- lidé MeSH
- molekulární modely MeSH
- počítačová simulace MeSH
- ribozomy chemie MeSH
- RNA katalytická chemie MeSH
- RNA virová chemie MeSH
- RNA chemie MeSH
- software MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- práce podpořená grantem MeSH
- přehledy MeSH
- Research Support, N.I.H., Extramural MeSH
One of the most common methods to treat the electrostatic effect of the environment in QM/MM calculations is to include the MM atoms as point charges in the QM Hamiltonian. In this case, a microiterative geometry optimization ignoring the QM contributions to the forces in the relaxation of the environment cannot yield exact stationary points. One solution that has been suggested in the literature is based on using a constant additive correction to the MM gradient during the microiterations, determined in the preceding macroiteration. Here, we analyze the convergence properties of the gradient correction method and point out that a smooth relaxation is not ensured if the curvature of the approximate, MM-based description of the potential energy surface of the environment is too small in comparison with the exact one. We suggest a computationally cheap second-order correction that uses an estimated Hessian from the Davidon-Fletcher-Powell method to tackle the problems caused by the too small curvature. Test calculations on four metalloenzymatic systems (∼100 QM atoms, ∼2000 relaxed MM atoms, ∼20,000 atoms in total) show that our approach efficiently restores the convergence where gradient correction alone would lead to oscillations.
- MeSH
- kvantová teorie MeSH
- simulace molekulární dynamiky MeSH
- software MeSH
- výpočetní biologie metody MeSH
- Publikační typ
- časopisecké články MeSH
The glycosylation of cell surface proteins plays a crucial role in a multitude of biological processes, such as cell adhesion and recognition. To understand the process of protein glycosylation, the reaction mechanisms of the participating enzymes need to be known. However, the reaction mechanism of retaining glycosyltransferases has not yet been sufficiently explained. Here we investigated the catalytic mechanism of human isoform 2 of the retaining glycosyltransferase polypeptide UDP-GalNAc transferase by coupling two different QM/MM-based approaches, namely a potential energy surface scan in two distance difference dimensions and a minimum energy reaction path optimisation using the Nudged Elastic Band method. Potential energy scan studies often suffer from inadequate sampling of reactive processes due to a predefined scan coordinate system. At the same time, path optimisation methods enable the sampling of a virtually unlimited number of dimensions, but their results cannot be unambiguously interpreted without knowledge of the potential energy surface. By combining these methods, we have been able to eliminate the most significant sources of potential errors inherent to each of these approaches. The structural model is based on the crystal structure of human isoform 2. In the QM/MM method, the QM region consists of 275 atoms, the remaining 5776 atoms were in the MM region. We found that ppGalNAcT2 catalyzes a same-face nucleophilic substitution with internal return (SNi). The optimized transition state for the reaction is 13.8 kcal/mol higher in energy than the reactant while the energy of the product complex is 6.7 kcal/mol lower. During the process of nucleophilic attack, a proton is synchronously transferred to the leaving phosphate. The presence of a short-lived metastable oxocarbenium intermediate is likely, as indicated by the reaction energy profiles obtained using high-level density functionals.
- MeSH
- algoritmy MeSH
- chemické modely * MeSH
- glykosylace MeSH
- glykosyltransferasy chemie ultrastruktura MeSH
- katalýza MeSH
- kinetika MeSH
- konformace proteinů MeSH
- molekulární modely * MeSH
- počítačová simulace MeSH
- polysacharidy chemie ultrastruktura MeSH
- vazba proteinů MeSH
- vazebná místa MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
Deltahedral metallacarborane compounds have recently been discovered as potent, specific, stable, and nontoxic inhibitors of HIV-1 protease (PR), the major target for AIDS therapy. The 2.15 A-resolution X-ray structure has exhibited a nonsymmetrical binding of the parental compound [Co(3+)-(C2B9H11)2](-) (GB-18) into PR dimer and a symmetrical arrangement in the crystal of two PR dimer complexes into a tetramer. In order to explore structural and energetic details of the inhibitor binding, quantum mechanics coupled with molecular mechanics approach was utilized. Realizing the close positioning of anionic inhibitors in the active site cavity, the possibility of an exchange of structural water molecules Wat50 and Wat128 by Na+ counterions was studied. The energy profiles for the rotation of the GB-18 molecules along their longitudinal axes in complex with PR were calculated. The results show that two Na+ counterions are present in the active site cavity and provide energetically favorable and unfavorable positions for carbon atoms within the carborane cages. Eighty-one rotamer combinations of four molecules of GB-18 bound to PR out of 4 x 10(5) are predicted to be highly populated. These results lay ground for further calculations of interaction energies between GB-18 and amino acids of PR active site and will make it possible to interpret computationally the binding of similar metallacarborane molecules to PR as well as to resistant PR variants. Moreover, this computational tool will allow the design of new, more potent metallacarborane-based HIV-1 protease inhibitors.
Organophosphorus compounds (OP) nerve agents are among the most toxic chemical substances known. Their toxicity is due to their ability to bind to acetylcholinesterase. Currently, some enzymes, such as phosphotriesterase, human serum paraoxonase 1 and diisopropyl fluorophosphatase, capable of degrading OP, have been characterized. Regarding the importance of bioremediation methods for detoxication of OP, this work aims to study the interaction modes between the human human deoxyuridine triphosphate nucleotidohydrolase (dUTPase) and Sarin and VX, considering their Rp and Sp enantiomers, to evaluate the asymmetric catalysis of those compounds. In previous work, this enzyme has shown good potential to degrade phosphotriesters, and based on this characteristic, we have applied the human dUTPase to the OP degradation. Molecular docking, chemometrics and mixed quantum and molecular mechanics calculations have been employed, showing a good interaction between dUTPase and OP. Two possible reaction mechanisms were tested, and according to our theoretical results, the catalytic degradation of OP by dUTPase can take place via both mechanisms, beyond being stereoselective, that is, dUTPase cleaves one enantiomer preferentially in relation to other. Chemometric techniques provided excellent assistance for performing this theoretical investigation. The dUTPase study shows importance by the fact of it being a human enzyme. Communicated by Ramaswamy H. Sarma.
- MeSH
- analýza hlavních komponent MeSH
- biodegradace MeSH
- katalytická doména MeSH
- kvantová teorie * MeSH
- lidé MeSH
- nervová bojová látka chemie metabolismus MeSH
- organofosforové sloučeniny chemie metabolismus MeSH
- organothiofosforové sloučeniny chemie metabolismus MeSH
- pyrofosfatasy metabolismus MeSH
- sarin chemie metabolismus MeSH
- simulace molekulového dockingu * MeSH
- vodíková vazba MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
Mechanistic studies on the hydrolytic dehalogenation catalyzed by haloalkane dehalogenases are of importance for environmental and industrial applications. Here, Car-Parrinello (CP) and ONIOM hybrid quantum-mechanical/molecular mechanics (QM/MM) are used investigate the second reaction step of the catalytic cycle, which comprises a general base-catalyzed hydrolysis of an ester intermediate (EI) to alcohol and free enzyme. We focus on the enzyme LinB from Sphingomonas paucimobilis UT26, for which the X-ray structure at atomic resolution is available. In agreement with previous proposals, our calculations suggest that a histidine residue (His272), polarized by glutamate (Glu132), acts as a base, accepting a proton from the catalytic water molecule and transferring it to an alcoholate ion. The reaction proceeds through a metastable tetrahedral intermediate, which shows an easily reversed reaction to the EI. In the formation of the products, the protonated aspartic acid (Asp108) can easily adopt conformation of the relaxed state found in the free enzyme. The overall free energy barrier of the reaction calculated by potential of the mean force integration using CP-QM/MM calculations is equal to 19.5 +/- 2 kcal . mol(-1). The lowering of the energy barrier of catalyzed reaction with respect to the water reaction is caused by strong stabilization of the reaction intermediate and transition state and their preorganization by electrostatic field of the enzyme. (c) 2007 Wiley-Liss, Inc.
- MeSH
- bakteriální proteiny chemie metabolismus MeSH
- financování organizované MeSH
- halogenace MeSH
- hydrolasy chemie metabolismus MeSH
- hydrolýza MeSH
- katalýza MeSH
- kinetika MeSH
- konformace proteinů MeSH
- kvantová teorie MeSH
- molekulární modely MeSH
- Sphingomonas enzymologie metabolismus MeSH
- termodynamika MeSH
- vazebná místa MeSH
The hepatitis delta virus (HDV) ribozyme is an RNA motif embedded in human pathogenic HDV RNA. Previous experimental studies have established that the active-site nucleotide C75 is essential for self-cleavage of the ribozyme, although its exact catalytic role in the process remains debated. Structural data from X-ray crystallography generally indicate that C75 acts as the general base that initiates catalysis by deprotonating the 2'-OH nucleophile at the cleavage site, while a hydrated magnesium ion likely protonates the 5'-oxygen leaving group. In contrast, some mechanistic studies support the role of C75 acting as general acid and thus being protonated before the reaction. We report combined quantum chemical/molecular mechanical calculations for the C75 general base pathway, utilizing the available structural data for the wild type HDV genomic ribozyme as a starting point. Several starting configurations differing in magnesium ion placement were considered and both one-dimensional and two-dimensional potential energy surface scans were used to explore plausible reaction paths. Our calculations show that C75 is readily capable of acting as the general base, in concert with the hydrated magnesium ion as the general acid. We identify a most likely position for the magnesium ion, which also suggests it acts as a Lewis acid. The calculated energy barrier of the proposed mechanism, approximately 20 kcal/mol, would lower the reaction barrier by approximately 15 kcal/mol compared with the uncatalyzed reaction and is in good agreement with experimental data.
In numerous Gram-positive bacteria, the glmS ribozyme or catalytic riboswitch regulates the expression of glucosamine-6-phosphate (GlcN6P) synthase via site-specific cleavage of its sugar-phosphate backbone in response to GlcN6P ligand binding. Biochemical data have suggested a crucial catalytic role for an active site guanine (G40 in Thermoanaerobacter tengcongensis, G33 in Bacillus anthracis). We used hybrid quantum chemical/molecular mechanical (QM/MM) calculations to probe the mechanism where G40 is deprotonated and acts as a general base. The calculations suggest that the deprotonated guanine G40(-) is sufficiently reactive to overcome the thermodynamic penalty arising from its rare protonation state, and thus is able to activate the A-1(2'-OH) group toward nucleophilic attack on the adjacent backbone. Furthermore, deprotonation of A-1(2'-OH) and nucleophilic attack are predicted to occur as separate steps, where activation of A-1(2'-OH) precedes nucleophilic attack. Conversely, the transition state associated with the rate-determining step corresponds to concurrent nucleophilic attack and protonation of the G1(O5') leaving group by the ammonium moiety of the GlcN6P cofactor. Overall, our calculations help to explain the crucial roles of G40 (as a general base) and GlcN6P (as a general acid) during glmS ribozyme self-cleavage. In addition, we show that the QM/MM description of the glmS ribozyme self-cleavage reaction is significantly more sensitive to the size of the QM region and the quality of the QM-MM coupling than that of other small ribozymes.
The excision of 8-oxoguanine (oxoG) by the human 8-oxoguanine DNA glycosylase 1 (hOGG1) base-excision repair enzyme was studied by using the QM/MM (M06-2X/6-31G(d,p):OPLS2005) calculation method and nuclear magnetic resonance (NMR) spectroscopy. The calculated glycosylase reaction included excision of the oxoG base, formation of Lys249-ribose enzyme-substrate covalent adduct and formation of a Schiff base. The formation of a Schiff base with ΔG# = 17.7 kcal/mol was the rate-limiting step of the reaction. The excision of the oxoG base with ΔG# = 16.1 kcal/mol proceeded via substitution of the C1΄-N9 N-glycosidic bond with an H-N9 bond where the negative charge on the oxoG base and the positive charge on the ribose were compensated in a concerted manner by NH3+(Lys249) and CO2-(Asp268), respectively. The effect of Asp268 on the oxoG excision was demonstrated with 1H NMR for WT hOGG1 and the hOGG1(D268N) mutant: the excision of oxoG was notably suppressed when Asp268 was mutated to Asn. The loss of the base-excision function was rationalized with QM/MM calculations and Asp268 was confirmed as the electrostatic stabilizer of ribose oxocarbenium through the initial base-excision step of DNA repair. The NMR experiments and QM/MM calculations consistently illustrated the base-excision reaction operated by hOGG1.
- MeSH
- biokatalýza MeSH
- DNA-glykosylasy chemie metabolismus MeSH
- guanin analogy a deriváty metabolismus MeSH
- kyselina asparagová metabolismus MeSH
- lidé MeSH
- lysin metabolismus MeSH
- molekulární modely MeSH
- mutantní proteiny chemie metabolismus MeSH
- oprava DNA * MeSH
- protonová magnetická rezonanční spektroskopie MeSH
- termodynamika MeSH
- vztahy mezi strukturou a aktivitou MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
