Electrostatically Embedded Molecular Tailoring Approach and Validation for Peptides.

We add higher-order electronic polarization effects to the molecular tailoring approach (MTA) by embedding each fragment in background charges as in combined quantum mechanical and molecular mechanical (QM/MM) methods; the resulting method considered here is called electrostatically embedded MTA (EE-MTA). We compare EE-MTA to MTA for a test peptide, Ace-(Ala)20-NMe, and we find that including background charges (embedding charges) greatly improves the performance. The fragmentation is performed on the basis of amino acids as monomers, and several sizes of fragment are tested. The fragments are capped with either hydrogen cap atoms or tuned fluorine cap atoms. The effective core potential of the tuned fluorine cap atom is optimized so as to reproduce the proton affinity for seven types of tetrapeptide, and the proton affinity calculated with the tuned cap atom shows a lower mean unsigned error than that obtained by using a hydrogen cap atom. In the application to the test peptide, these generically tuned cap atoms show better performance compared with the hydrogen cap atom for both the electronic energy and the energy difference between an α helix and a β sheet (in the latter case, 1.0% vs 2.7% when averaged over three sizes of fragments and two locations for cut bonds). Also, we compared the accuracy of several charge redistribution schemes, and we find that the results are not particularly sensitive to these choices for the Ace-(Ala)20-NMe peptide. We also illustrate the dependence on the choice of charge model for the embedding charges, including both fixed embedding charges and embedding charges that depend on conformation.

[1]  Ross C. Walker,et al.  The implementation of a fast and accurate QM/MM potential method in Amber , 2008, J. Comput. Chem..

[2]  P. Kollman,et al.  Application of RESP charges to calculate conformational energies, hydrogen bond energies, and free energies of solvation , 1993 .

[3]  R. S. Mulliken.  Electronic Population Analysis on LCAO–MO Molecular Wave Functions. I , 1955 .

[4]  Peter A. Kollman,et al.  Calculation of molecular geometries, relative conformational energies, dipole moments, and molecular electrostatic potential fitted charges of small organic molecules of biochemical interest by density functional theory , 1995, J. Comput. Chem..

[5]  Jiali Gao,et al.  Toward a Molecular Orbital Derived Empirical Potential for Liquid Simulations , 1997 .

[6]  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.

[7]  Donald G Truhlar,et al.  Electrostatically Embedded Many-Body Approximation for Systems of Water, Ammonia, and Sulfuric Acid and the Dependence of Its Performance on Embedding Charges. , 2009, Journal of chemical theory and computation.

[8]  K. Kitaura,et al.  Multilayer formulation of the fragment molecular orbital method (FMO). , 2005, The journal of physical chemistry. A.

[9]  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.

[10]  Nan Jiang,et al.  Electrostatic field-adapted molecular fractionation with conjugated caps for energy calculations of charged biomolecules. , 2006, The Journal of chemical physics.

[11]  U. Singh,et al.  A combined ab initio quantum mechanical and molecular mechanical method for carrying out simulations on complex molecular systems: Applications to the CH3Cl + Cl− exchange reaction and gas phase protonation of polyethers , 1986 .

[12]  Kazuo Kitaura,et al.  Structural and interaction analysis of helical heparin oligosaccharides with the fragment molecular orbital method , 2009 .

[13]  P. Kollman,et al.  How well does a restrained electrostatic potential (RESP) model perform in calculating conformational energies of organic and biological molecules? , 2000 .

[14]  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..

[15]  D. Suárez,et al.  Thermochemical Fragment Energy Method for Biomolecules: Application to a Collagen Model Peptide. , 2009, Journal of chemical theory and computation.

[16]  Margaret M. Hurley,et al.  Simple one-electron quantum capping potentials for use in hybrid QM/MM studies of biological molecules , 2002 .

[17]  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.

[18]  Marat Valiev,et al.  Fast electron correlation methods for molecular clusters in the ground and excited states , 2005 .

[19]  Michael J. Frisch,et al.  Achieving Linear Scaling for the Electronic Quantum Coulomb Problem , 1996, Science.

[20]  Gregory J O Beran,et al.  Approximating quantum many-body intermolecular interactions in molecular clusters using classical polarizable force fields. , 2009, The Journal of chemical physics.

[21]  Weitao Yang,et al.  A density‐matrix divide‐and‐conquer approach for electronic structure calculations of large molecules , 1995 .

[22]  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.

[23]  Donald G Truhlar,et al.  Evaluation of the Electrostatically Embedded Many-Body Expansion and the Electrostatically Embedded Many-Body Expansion of the Correlation Energy by Application to Low-Lying Water Hexamers. , 2008, Journal of chemical theory and computation.

[24]  Frank Weinhold,et al.  Natural hybrid orbitals , 1980 .

[25]  Kazuo Kitaura,et al.  Cluster hydration model for binding energy calculations of protein-ligand complexes. , 2009, The journal of physical chemistry. B.

[26]  Shridhar R Gadre,et al.  On the applicability of fragmentation methods to conjugated pi systems within density functional framework. , 2010, The Journal of chemical physics.

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

[28]  A. E. Hirsh,et al.  The application of statistical physics to evolutionary biology. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[29]  Wei Li.  Linear scaling explicitly correlated MP2-F12 and ONIOM methods for the long-range interactions of the nanoscale clusters in methanol aqueous solutions. , 2013, The Journal of chemical physics.

[30]  Notker Rösch,et al.  Elastic Polarizable Environment Cluster Embedding Approach for Covalent Oxides and Zeolites Based on a Density Functional Method , 2003 .

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

[32]  G Narahari Sastry,et al.  Ab initio investigation of benzene clusters: molecular tailoring approach. , 2010, The Journal of chemical physics.

[33]  J. Karle,et al.  Kernel energy method: Basis functions and quantum methods† , 2006 .

[34]  Lu,et al.  Density-functional calculations of the structure and stability of C240. , 1994, Physical review. B, Condensed matter.

[35]  N. Burton,et al.  Catalytic Mechanism of the Enzyme Papain: Predictions with a Hybrid Quantum Mechanical/Molecular Mechanical Potential , 1997 .

[36]  Yingkai Zhang,et al.  Improved pseudobonds for combined ab initio quantum mechanical/molecular mechanical methods. , 2005, The Journal of chemical physics.

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

[38]  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.

[39]  Yuto Komeiji,et al.  Fragment molecular orbital method: analytical energy gradients , 2001 .

[40]  Jan H. Jensen,et al.  Covalent bond fragmentation suitable to describe solids in the fragment molecular orbital method. , 2008, The journal of physical chemistry. A.

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

[42]  Anuja P. Rahalkar,et al.  Enabling ab initio Hessian and frequency calculations of large molecules. , 2008, The Journal of chemical physics.

[43]  D. Truhlar,et al.  QM/MM: what have we learned, where are we, and where do we go from here? , 2007 .

[44]  D. Truhlar,et al.  Combined Quantum Mechanical and Molecular Mechanical Methods for Calculating Potential Energy Surfaces: Tuned and Balanced Redistributed-Charge Algorithm. , 2010, Journal of chemical theory and computation.

[45]  Wei Zhang,et al.  A point‐charge force field for molecular mechanics simulations of proteins based on condensed‐phase quantum mechanical calculations , 2003, J. Comput. Chem..

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

[47]  Wei Li,et al.  Generalized energy-based fragmentation approach for computing the ground-state energies and properties of large molecules. , 2007, The journal of physical chemistry. A.

[48]  D. Truhlar,et al.  Geometry optimization using tuned and balanced redistributed charge schemes for combined quantum mechanical and molecular mechanical calculations. , 2011, Physical chemistry chemical physics : PCCP.

[49]  Tai-Sung Lee,et al.  A pseudobond approach to combining quantum mechanical and molecular mechanical methods , 1999 .

[50]  Nicholas J Mayhall,et al.  Molecules-in-Molecules: An Extrapolated Fragment-Based Approach for Accurate Calculations on Large Molecules and Materials. , 2011, Journal of chemical theory and computation.

[51]  Lu,et al.  Shape of large single- and multiple-shell fullerenes. , 1994, Physical review. B, Condensed matter.

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

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

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

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

[56]  Yutaka Akiyama,et al.  Fragment molecular orbital method: application to polypeptides , 2000 .

[57]  Yang,et al.  Direct calculation of electron density in density-functional theory. , 1991, Physical review letters.

[58]  Mark S. Gordon,et al.  The Effective Fragment Potential Method: A QM-Based MM Approach to Modeling Environmental Effects in Chemistry , 2001 .

[59]  Gregory J O Beran,et al.  Practical quantum mechanics-based fragment methods for predicting molecular crystal properties. , 2012, Physical chemistry chemical physics : PCCP.

[60]  Michael A. Collins,et al.  Accurate treatment of nonbonded interactions within systematic molecular fragmentation , 2009 .

[61]  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.

[62]  Weitao Yang,et al.  The divide‐and‐conquer density‐functional approach: Molecular internal rotation and density of states , 1992 .

[63]  Kazuya Ishimura,et al.  Accuracy of the three‐body fragment molecular orbital method applied to Møller–Plesset perturbation theory , 2007, J. Comput. Chem..

[64]  Mark S Gordon,et al.  The fragment molecular orbital and systematic molecular fragmentation methods applied to water clusters. , 2012, Physical chemistry chemical physics : PCCP.

[65]  K. Morokuma,et al.  A simple scheme of estimating substitution or substituent effects in the ab initio MO method based on the shift operator , 1990 .

[66]  Kazuo Kitaura,et al.  The importance of three-body terms in the fragment molecular orbital method. , 2004, The Journal of chemical physics.

[67]  D. Truhlar,et al.  Multilevel X-Pol: a fragment-based method with mixed quantum mechanical representations of different fragments. , 2012, The journal of physical chemistry. B.

[68]  P. Kollman,et al.  Atomic charges derived from semiempirical methods , 1990 .

[69]  Yang,et al.  Direct calculation of electron density in density-functional theory: Implementation for benzene and a tetrapeptide. , 1991, Physical review. A, Atomic, molecular, and optical physics.

[70]  S. R. Gadre,et al.  Molecular cluster building algorithm: electrostatic guidelines and molecular tailoring approach. , 2011, The Journal of chemical physics.

[71]  Jean-Louis Rivail,et al.  Quantum chemical computations on parts of large molecules: the ab initio local self consistent field method , 1996 .

[72]  Shridhar R. Gadre,et al.  Molecular tailoring approach in conjunction with MP2 and Ri‐MP2 codes: A comparison with fragment molecular orbital method , 2010, J. Comput. Chem..

[73]  Benny G. Johnson,et al.  THE CONTINUOUS FAST MULTIPOLE METHOD , 1994 .

[74]  Shuhua Li,et al.  Cooperativity in long α- and 3(10)-helical polyalanines: both electrostatic and van der Waals interactions are essential. , 2011, The journal of physical chemistry. B.

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

[76]  R. P. Bettens,et al.  Accurately reproducing ab initio electrostatic potentials with multipoles and fragmentation. , 2009, The journal of physical chemistry. A.

[77]  Weitao Yang,et al.  Analytical energy gradients and geometry optimization in the divide‐and‐conquer method for large molecules , 1995 .

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

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

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

[81]  P. Löwdin.  On the Non‐Orthogonality Problem Connected with the Use of Atomic Wave Functions in the Theory of Molecules and Crystals , 1950 .

[82]  Donald G Truhlar,et al.  Charge Model 5: An Extension of Hirshfeld Population Analysis for the Accurate Description of Molecular Interactions in Gaseous and Condensed Phases. , 2012, Journal of chemical theory and computation.

[83]  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.

[84]  M. Field,et al.  A Generalized Hybrid Orbital (GHO) Method for the Treatment of Boundary Atoms in Combined QM/MM Calculations , 1998 .

[85]  Donald G Truhlar,et al.  Redistributed charge and dipole schemes for combined quantum mechanical and molecular mechanical calculations. , 2005, The journal of physical chemistry. A.

[86]  Martin Head-Gordon,et al.  Derivation and efficient implementation of the fast multipole method , 1994 .