{"title":"Protein-based QM-CGMM","authors":"Yi Sun","doi":"10.1016/j.jmgm.2025.109081","DOIUrl":null,"url":null,"abstract":"<div><div>We present an approach to Coarse Grain (CG) the Molecular Mechanics (MM) region in a thermodynamically consistent approach when the region is a protein backbone. This approach is achieved, via using 4 beads on each amino acid residue to replicate the electrostatic field that it generates, followed by connecting them via a heterogeneous Elastic Network Model (heteroENM) force-field. We subsequently treat the boundary between the QM and CGMM regions using appropriate Lennard-Jones parameters. We apply this procedure to the classic Chorismate Mutase system, where the enzyme backbone in the MM region catalyzes the Claisen rearrangement in the QM region using Transition-Tempered metadynamics (TTMetaD) to obtain a converged Potential of Mean Force (PMF). The CGMM reveals a free energy barrier similar to the experimental value. This work demonstrates the potential of utilizing a CGMM forcefield to capture all-atom thermodynamic properties at a reduced cost.</div></div>","PeriodicalId":16361,"journal":{"name":"Journal of molecular graphics & modelling","volume":"140 ","pages":"Article 109081"},"PeriodicalIF":2.7000,"publicationDate":"2025-06-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of molecular graphics & modelling","FirstCategoryId":"99","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S109332632500141X","RegionNum":4,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"BIOCHEMICAL RESEARCH METHODS","Score":null,"Total":0}
引用次数: 0
Abstract
We present an approach to Coarse Grain (CG) the Molecular Mechanics (MM) region in a thermodynamically consistent approach when the region is a protein backbone. This approach is achieved, via using 4 beads on each amino acid residue to replicate the electrostatic field that it generates, followed by connecting them via a heterogeneous Elastic Network Model (heteroENM) force-field. We subsequently treat the boundary between the QM and CGMM regions using appropriate Lennard-Jones parameters. We apply this procedure to the classic Chorismate Mutase system, where the enzyme backbone in the MM region catalyzes the Claisen rearrangement in the QM region using Transition-Tempered metadynamics (TTMetaD) to obtain a converged Potential of Mean Force (PMF). The CGMM reveals a free energy barrier similar to the experimental value. This work demonstrates the potential of utilizing a CGMM forcefield to capture all-atom thermodynamic properties at a reduced cost.
期刊介绍:
The Journal of Molecular Graphics and Modelling is devoted to the publication of papers on the uses of computers in theoretical investigations of molecular structure, function, interaction, and design. The scope of the journal includes all aspects of molecular modeling and computational chemistry, including, for instance, the study of molecular shape and properties, molecular simulations, protein and polymer engineering, drug design, materials design, structure-activity and structure-property relationships, database mining, and compound library design.
As a primary research journal, JMGM seeks to bring new knowledge to the attention of our readers. As such, submissions to the journal need to not only report results, but must draw conclusions and explore implications of the work presented. Authors are strongly encouraged to bear this in mind when preparing manuscripts. Routine applications of standard modelling approaches, providing only very limited new scientific insight, will not meet our criteria for publication. Reproducibility of reported calculations is an important issue. Wherever possible, we urge authors to enhance their papers with Supplementary Data, for example, in QSAR studies machine-readable versions of molecular datasets or in the development of new force-field parameters versions of the topology and force field parameter files. Routine applications of existing methods that do not lead to genuinely new insight will not be considered.