This work presents a systematic assessment of QM/QM' and QM/MM models with respect to direct QM calculations for the tautomerization (neutral to zwitterion) reactions of amino acids (glycine, alanine, valine, aspartate, and neutral and protonated histidine) solvated in a 160 water cluster. The effect of varying QM region size and choice of embedding potentials, including fixed-charge and polarizable molecular mechanics force fields (TIP3P and EFP) and various semiempirical QM methods (PM7, GFN2-xTB, DFTBA, DFTB3, HF-3c, and PBEh-3c), on the accuracy of the models was examined. A surprising finding was that molecular mechanics force fields outperformed many of the semiempirical methods. Generally, the errors in the QM/QM' and QM/MM models converge slowly with respect to the QM region size, requiring 50 or more waters to be included in the QM region before the error in the model falls below 1 kcal mol of its pure QM result. Different QM region selection schemes were also compared, and it was found that selection based on Natural Population Analysis (NPA) atomic charges significantly reduced the error in the QM/QM' and QM/MM models particularly if a low-quality embedding potential was used. It is envisaged that these results will be useful for the development of future hybrid QM models.
Download full-text PDF |
Source |
|---|---|
| http://dx.doi.org/10.1021/acs.jpcb.1c04876 | DOI Listing |
Publication Analysis
Top Keywords
Similar Publications
Activation and Reactivity of the Deubiquitinylase OTU Cezanne-2 from MD Simulations and QM/MM Calculations.
J Chem Inf Model
January 2025
Molecular Simulations and Design Group, Max Planck Institute for Dynamics of Complex Technical Systems, Sandtorstrasse 1, 39106 Magdeburg, Germany.
Cezanne-2 (Cez2) is a deubiquitinylating (DUB) enzyme involved in the regulation of ubiquitin-driven cellular signaling and selectively targets Lys11-linked polyubiquitin chains. As a representative member of the ovarian tumor (OTU) subfamily DUBs, it performs cysteine proteolytic isopeptide bond cleavage; however, its exact catalytic mechanism is not yet resolved. In this work, we used different computational approaches to get molecular insights into the Cezanne-2 catalytic mechanism.
View Article and Find Full Text PDFDissecting Methionine Oxidation by Hydrogen Peroxide in Proteins: Thermodynamics, Kinetics, and Susceptibility Descriptors.
J Chem Inf Model
January 2025
Department of Chemistry, University of Rome, Sapienza, P.le A. Moro 5, 00185 Rome, Italy.
The oxidation of Met residues in proteins is a complex process, where protein-specific structural and dynamical features play a relevant role in determining the reaction kinetics. Aiming to a full-side perspective, we report here a comprehensive characterization of Met oxidation kinetics by hydrogen peroxide in a leptin protein case study. To do that, we estimated the reaction-free energy profile of the Met oxidation via a QM/MM approach, while the kinetics of the formation of the reactive species were calculated using classical molecular dynamics (MD) simulations.
View Article and Find Full Text PDFModulating Charge Transfer Kinetics along Poly Adenine: Chemical Modifications, Temperature, and Conformational Effects.
J Chem Theory Comput
January 2025
Department of Chemistry, University of Rome, Sapienza, P.le A. Moro 5, 00185 Rome, Italy.
The charge transfer (CT) reactions in nucleic acids are crucial for genome damage and repair and nanoelectronics using DNA as a molecular conductor. Previous experimental and theoretical works underlined the significance of nucleic acid structural dynamics on CT kinetics, requiring models that incorporate the dynamics of the nucleic acid, solvents, and counterions. Here, we investigated hole transfer kinetics in poly adenine single and double strands at various temperatures and the rate enhancement due to adenine-to-7-deazaadenine mutation by means of a QM/MM approach.
View Article and Find Full Text PDFA Closer Look at the FeS Heme Bonds in Azotobacter vinelandii Bacterioferritin: QM/MM and Local Mode Analysis.
J Comput Chem
January 2025
Chemistry Department, Southern Methodist University, Dallas, Texas, USA.
Using the QM/MM methodology and a local mode analysis, we investigated a character and a strength of FeS bonds of heme groups in oxidized and reduced forms of Bacterioferritin from Azotobacter vinelandii. The strength of the FeS bonds was correlated with a bond length, an energy density at a bond critical point, and a charge difference of the F and S atoms. Changing the oxidation state from ferrous to ferric generally makes the FeS bonds weaker, longer, more covalent, and more polar.
View Article and Find Full Text PDFNAC4ED: A high-throughput computational platform for the rational design of enzyme activity and substrate selectivity.
mLife
December 2024
Shanghai Engineering Research Center of Molecular Therapeutics & New Drug Development School of Chemistry and Molecular Engineering, East China Normal University Shanghai China.
In silico computational methods have been widely utilized to study enzyme catalytic mechanisms and design enzyme performance, including molecular docking, molecular dynamics, quantum mechanics, and multiscale QM/MM approaches. However, the manual operation associated with these methods poses challenges for simulating enzymes and enzyme variants in a high-throughput manner. We developed the NAC4ED, a high-throughput enzyme mutagenesis computational platform based on the "near-attack conformation" design strategy for enzyme catalysis substrates.
View Article and Find Full Text PDFWant AI Summaries of new PubMed Abstracts delivered to your In-box?
Enter search terms and have AI summaries delivered each week - change queries or unsubscribe any time!