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
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
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
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
Ru(II) "piano-stool" complexes belong to group of biologically active metallocomplexes with promising anticancer activity. In this study, we investigate the reaction mechanism of [(η(6)-benzene)Ru(II)(en)(H(2)O)](2+) (en = ethylenediamine) complex binding to DNA by hybrid QM/MM computational techniques. The reaction when the Ru(II) complex is coordinated on N7-guanine from major groove is explored. Two reaction pathways, direct binding to N7 position and two-step mechanism passing through O6 position, are considered. It was found that the reaction is exothermic and the direct binding process is preferred kinetically. In analogy to cisplatin, we also explored the possibility of intrastrand cross-link formation where the Ru(II) complex makes a bridge between two adjacent guanines. Two different pathways were found, leading to a final structure with released benzene ligand. This process is exothermic; however, one pathway is blocked by relatively high initial activation barrier. Geometries, energies, and electronic properties analyzed by atoms in molecules and natural population analysis methods are discussed.
- MeSH
- hydroxid vápenatý * farmakologie terapeutické užití MeSH
- lidé MeSH
- zubní cementy * MeSH
- Check Tag
- lidé 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
The ISWI family protein SMARCA5 contains the ATP-binding pocket that coordinates the catalytic Mg2+ ion and water molecules for ATP hydrolysis. In this study, we demonstrate that SMARCA5 can also possess an alternative metal-binding ability. First, we isolated SMARCA5 on the cobalt column (IMAC) to near homogeneity. Examination of the interactions of SMARCA5 with metal-chelating supports showed that, apart from Co2+, it binds to Cu2+, Zn2+ and Ni2+. The efficiency of the binding to the last-listed metal was influenced by the chelating ligand, resulting in a strong preference for Ni-NTA over the Ni-CM-Asp equivalent. To gain insight in the preferential affinity for the Ni-NTA ligand, QM calculations were performed on model systems and metal-ligand complexes with a limited protein fragment of SMARCA5 containing the double-histidine (dHis) motif. The calculations correlated the observed affinity with the relative stability of the d-block metals to tetradentate ligand coordination over tridentate, as well as their overall octahedral coordination capacity. Likewise, binding free energies derived from model imidazole complexes mirrored the observed Ni-NTA/Ni-CM-Asp preferential affinity. Finally, similar calculations on complexes with a SMARCA5 peptide fragment derived from the AlphaFold structural prediction, captured almost accurately the expected relative stability of the TM complexes, and produced a large energetic separation (~10 kcal∙mol-1) between Ni-NTA and Ni-CM-Asp in favour of the former.
- MeSH
- adenosintrifosfatasy MeSH
- chromozomální proteiny, nehistonové metabolismus chemie MeSH
- kovy chemie metabolismus MeSH
- kvantová teorie MeSH
- lidé MeSH
- ligandy MeSH
- molekulární modely MeSH
- restrukturace chromatinu MeSH
- vazba proteinů * MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
