Full QM Calculation of RNA Energy Using Electrostatically Embedded Generalized Molecular Fractionation with Conjugate Caps Method.

In this study, the electrostatically embedded generalized molecular fractionation with conjugate caps (concaps) method (EE-GMFCC) was employed for efficient linear-scaling quantum mechanical (QM) calculation of total energies of RNAs. In the EE-GMFCC approach, the total energy of RNA is calculated by taking a proper combination of the QM energy of each nucleotide-centric fragment with large caps or small caps (termed EE-GMFCC-LC and EE-GMFCC-SC, respectively) deducted by the energies of concaps. The two-body QM interaction energy between non-neighboring ribonucleotides which are spatially in close contact are also taken into account for the energy calculation. Numerical studies were carried out to calculate the total energies of a number of RNAs using the EE-GMFCC-LC and EE-GMFCC-SC methods at levels of the Hartree-Fock (HF) method, density functional theory (DFT), and second-order many-body perturbation theory (MP2), respectively. The results show that the efficiency of the EE-GMFCC-SC method is about 3 times faster than the EE-GMFCC-LC method with minimal accuracy sacrifice. The EE-GMFCC-SC method is also applied for relative energy calculations of 20 different conformers of two RNA systems using HF and DFT, respectively. Both single-point and relative energy calculations demonstrate that the EE-GMFCC method has deviations from the full system results of only a few kcal/mol.

[1]  Anuja P. Rahalkar,et al.  Structure, energetics, and reactivity of boric acid nanotubes: a molecular tailoring approach. , 2008, The journal of physical chemistry. A.

[2]  Kaori Fukuzawa,et al.  Fragment molecular orbital method: use of approximate electrostatic potential , 2002 .

[3]  Xiao He,et al.  A systematic study on RNA NMR chemical shift calculation based on the automated fragmentation QM/MM approach , 2016 .

[4]  Donald G Truhlar,et al.  Electrostatically Embedded Many-Body Correlation Energy, with Applications to the Calculation of Accurate Second-Order Møller-Plesset Perturbation Theory Energies for Large Water Clusters. , 2007, Journal of chemical theory and computation.

[5]  J. P. Loria,et al.  MoD-QM/MM Structural Refinement Method: Characterization of Hydrogen Bonding in the Oxytricha nova G-Quadruplex. , 2014, Journal of chemical theory and computation.

[6]  Xiao He,et al.  Automated Fragmentation QM/MM Calculation of Amide Proton Chemical Shifts in Proteins with Explicit Solvent Model. , 2013, Journal of chemical theory and computation.

[7]  M. Frank-Kamenetskii,et al.  Base-stacking and base-pairing contributions into thermal stability of the DNA double helix , 2006, Nucleic acids research.

[8]  Michael A Collins,et al.  Approximate ab initio energies by systematic molecular fragmentation. , 2005, The Journal of chemical physics.

[9]  Xiao He,et al.  Fragment quantum mechanical calculation of proteins and its applications. , 2014, Accounts of chemical research.

[10]  Ryan P A Bettens,et al.  The combined fragmentation and systematic molecular fragmentation methods. , 2014, Accounts of chemical research.

[11]  B. Brooks,et al.  An analysis of the accuracy of Langevin and molecular dynamics algorithms , 1988 .

[12]  Lou Massa,et al.  Kernel energy method: application to DNA. , 2005, Biochemistry.

[13]  Shinichiro Nakamura,et al.  Analytic second derivatives of the energy in the fragment molecular orbital method. , 2013, The Journal of chemical physics.

[14]  John F. Ouyang,et al.  Many-Body Basis Set Superposition Effect. , 2015, Journal of chemical theory and computation.

[15]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[16]  John Z. H. Zhang,et al.  Calculation of protein–ligand binding affinities based on a fragment quantum mechanical method , 2015 .

[17]  Shridhar R. Gadre,et al.  Molecular Tailoring Approach for Simulation of Electrostatic Properties , 1994 .

[18]  Wei Li,et al.  Geometry optimizations and vibrational spectra of large molecules from a generalized energy-based fragmentation approach. , 2008, The journal of physical chemistry. A.

[19]  Donald G Truhlar,et al.  Electrostatically Embedded Many-Body Expansion for Large Systems, with Applications to Water Clusters. , 2007, Journal of chemical theory and computation.

[20]  Ryan P A Bettens,et al.  Energy-Based Molecular Fragmentation Methods. , 2015, Chemical reviews.

[21]  Xiao He,et al.  Fragment quantum chemical approach to geometry optimization and vibrational spectrum calculation of proteins. , 2016, Physical chemistry chemical physics : PCCP.

[22]  Xiao He,et al.  Electrostatically embedded generalized molecular fractionation with conjugate caps method for full quantum mechanical calculation of protein energy. , 2013, The journal of physical chemistry. A.

[23]  Michael W. Schmidt,et al.  Analytic second derivative of the energy for density functional theory based on the three-body fragment molecular orbital method. , 2015, The Journal of chemical physics.

[24]  Holger Gohlke,et al.  The Amber biomolecular simulation programs , 2005, J. Comput. Chem..

[25]  Lou Massa,et al.  Kernel energy method illustrated with peptides , 2005 .

[26]  Ryan P A Bettens,et al.  A new algorithm for molecular fragmentation in quantum chemical calculations. , 2006, The journal of physical chemistry. A.

[27]  Xiao He,et al.  Fragment density functional theory calculation of NMR chemical shifts for proteins with implicit solvation. , 2012, Physical chemistry chemical physics : PCCP.

[28]  W. L. Jorgensen,et al.  Temperature dependence of TIP3P, SPC, and TIP4P water from NPT Monte Carlo simulations: Seeking temperatures of maximum density , 1998 .

[29]  John M Herbert,et al.  A generalized many-body expansion and a unified view of fragment-based methods in electronic structure theory. , 2012, The Journal of chemical physics.

[30]  Donald G Truhlar,et al.  The variational explicit polarization potential and analytical first derivative of energy: Towards a next generation force field. , 2008, The Journal of chemical physics.

[31]  Yingkai Zhang,et al.  An efficient linear scaling method for ab initio calculation of electron density of proteins , 2004 .

[32]  R. Spizzo,et al.  Long non-coding RNAs and cancer: a new frontier of translational research? , 2012, Oncogene.

[33]  Lin He,et al.  MicroRNAs: small RNAs with a big role in gene regulation , 2004, Nature reviews genetics.

[34]  Michael A Collins.  Molecular forces, geometries, and frequencies by systematic molecular fragmentation including embedded charges. , 2014, The Journal of chemical physics.

[35]  Jianing Song,et al.  The critical effect of polarization on the dynamical structure of guanine quadruplex DNA. , 2013, Physical chemistry chemical physics : PCCP.

[36]  John Z H Zhang,et al.  Theoretical method for full ab initio calculation of DNA/RNA-ligand interaction energy. , 2004, The Journal of chemical physics.

[37]  Lou Massa,et al.  The kernel energy method of quantum mechanical approximation carried to fourth-order terms , 2008, Proceedings of the National Academy of Sciences.

[38]  U. Ryde,et al.  Ligand-Binding Affinity Estimates Supported by Quantum-Mechanical Methods. , 2016, Chemical reviews.

[39]  V Ganesh,et al.  Molecular tailoring approach for geometry optimization of large molecules: energy evaluation and parallelization strategies. , 2006, The Journal of chemical physics.

[40]  Xiao He,et al.  A new method for direct calculation of total energy of protein. , 2005, The Journal of chemical physics.

[41]  R. P. Bettens,et al.  Ab initio NMR chemical-shift calculations based on the combined fragmentation method. , 2013, Physical chemistry chemical physics : PCCP.

[42]  Michal Otyepka,et al.  Performance of Molecular Mechanics Force Fields for RNA Simulations: Stability of UUCG and GNRA Hairpins. , 2010, Journal of chemical theory and computation.

[43]  Shuhua Li,et al.  An efficient implementation of the generalized energy-based fragmentation approach for general large molecules. , 2010, The journal of physical chemistry. A.

[44]  John M Herbert,et al.  Many-Body Expansion with Overlapping Fragments: Analysis of Two Approaches. , 2013, Journal of chemical theory and computation.

[45]  Shridhar R Gadre,et al.  Vibrational infrared and Raman spectra of polypeptides: Fragments-in-fragments within molecular tailoring approach. , 2016, The Journal of chemical physics.

[46]  K. Merz,et al.  Protein NMR chemical shift calculations based on the automated fragmentation QM/MM approach. , 2009, The journal of physical chemistry. B.

[47]  Kazuo Kitaura,et al.  Extending the power of quantum chemistry to large systems with the fragment molecular orbital method. , 2007, The journal of physical chemistry. A.

[48]  V. Hornak,et al.  Comparison of multiple Amber force fields and development of improved protein backbone parameters , 2006, Proteins.

[49]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[50]  John F. Ouyang,et al.  Combined Fragmentation Method: A Simple Method for Fragmentation of Large Molecules. , 2012, Journal of chemical theory and computation.

[51]  John F. Ouyang,et al.  Trouble with the Many-Body Expansion. , 2014, Journal of chemical theory and computation.

[52]  Nicholas J Mayhall,et al.  Many-Overlapping-Body (MOB) Expansion: A Generalized Many Body Expansion for Nondisjoint Monomers in Molecular Fragmentation Calculations of Covalent Molecules. , 2012, Journal of chemical theory and computation.

[53]  Lou Massa,et al.  The Kernel Energy Method: Application to a tRNA , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[54]  Xiao He,et al.  The generalized molecular fractionation with conjugate caps/molecular mechanics method for direct calculation of protein energy. , 2006, The Journal of chemical physics.

[55]  Tai-Sung Lee,et al.  Linear-scaling quantum mechanical calculations of biological molecules: The divide-and-conquer approach , 1998 .

[56]  John Z H Zhang,et al.  Quantum mechanical map for protein-ligand binding with application to beta-trypsin/benzamidine complex. , 2004, The Journal of chemical physics.

[57]  John Z H Zhang,et al.  An efficient approach for ab initio energy calculation of biopolymers. , 2005, The Journal of chemical physics.

[58]  Ross C. Walker,et al.  An overview of the Amber biomolecular simulation package , 2013 .

[59]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[60]  Krishnan Raghavachari,et al.  Accurate Composite and Fragment-Based Quantum Chemical Models for Large Molecules. , 2015, Chemical reviews.

[61]  Jiali Gao,et al.  The Design of a Next Generation Force Field: The X-POL Potential. , 2007, Journal of chemical theory and computation.

[62]  Mark S Gordon,et al.  Systematic fragmentation method and the effective fragment potential: an efficient method for capturing molecular energies. , 2009, The journal of physical chemistry. A.

[63]  R. P. Bettens,et al.  First principles NMR calculations by fragmentation. , 2007, The journal of physical chemistry. A.

[64]  X. Chen,et al.  Molecular fractionation with conjugated caps density matrix with pairwise interaction correction for protein energy calculation. , 2006, The Journal of chemical physics.

[65]  John Z. H. Zhang,et al.  Molecular fractionation with conjugate caps for full quantum mechanical calculation of protein-molecule interaction energy , 2003 .

[66]  Xiao He,et al.  Quantum Fragment Based ab Initio Molecular Dynamics for Proteins. , 2015, Journal of chemical theory and computation.

[67]  J. Šponer,et al.  Refinement of the AMBER Force Field for Nucleic Acids: Improving the Description of α/γ Conformers , 2007 .

[68]  K. Kitaura,et al.  Fragment molecular orbital method: an approximate computational method for large molecules , 1999 .

[69]  Paul G. Mezey,et al.  Ab Initio-Quality Electrostatic Potentials for Proteins: An Application of the ADMA Approach , 2002 .

[70]  C. Croce,et al.  A microRNA expression signature of human solid tumors defines cancer gene targets , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[71]  P. Kollman,et al.  A well-behaved electrostatic potential-based method using charge restraints for deriving atomic char , 1993 .

[72]  Donny D. Licatalosi,et al.  RNA processing and its regulation: global insights into biological networks , 2010, Nature Reviews Genetics.

[73]  Michael A Collins,et al.  Accuracy and efficiency of electronic energies from systematic molecular fragmentation. , 2006, The Journal of chemical physics.

[74]  Spencer R Pruitt,et al.  Fragmentation methods: a route to accurate calculations on large systems. , 2012, Chemical reviews.

[75]  Marat Valiev,et al.  Fast electron correlation methods for molecular clusters without basis set superposition errors. , 2008, The Journal of chemical physics.

[76]  John M Herbert,et al.  Pair-Pair Approximation to the Generalized Many-Body Expansion: An Alternative to the Four-Body Expansion for ab Initio Prediction of Protein Energetics via Molecular Fragmentation. , 2016, Journal of chemical theory and computation.

[77]  J. Šponer,et al.  Refinement of the Cornell et al. Nucleic Acids Force Field Based on Reference Quantum Chemical Calculations of Glycosidic Torsion Profiles , 2011, Journal of chemical theory and computation.

[78]  Shridhar R. Gadre,et al.  Ab initio quality one‐electron properties of large molecules: Development and testing of molecular tailoring approach , 2003, J. Comput. Chem..