Extension of the Effective Fragment Potential Method to Macromolecules.

The effective fragment potential (EFP) approach, which can be described as a nonempirical polarizable force field, affords an accurate first-principles treatment of noncovalent interactions in extended systems. EFP can also describe the effect of the environment on the electronic properties (e.g., electronic excitation energies and ionization and electron-attachment energies) of a subsystem via the QM/EFP (quantum mechanics/EFP) polarizable embedding scheme. The original formulation of the method assumes that the system can be separated, without breaking covalent bonds, into closed-shell fragments, such as solvent and solute molecules. Here, we present an extension of the EFP method to macromolecules (mEFP). Several schemes for breaking a large molecule into small fragments described by EFP are presented and benchmarked. We focus on the electronic properties of molecules embedded into a protein environment and consider ionization, electron-attachment, and excitation energies (single-point calculations only). The model systems include chromophores of green and red fluorescent proteins surrounded by several nearby amino acid residues and phenolate bound to the T4 lysozyme. All mEFP schemes show robust performance and accurately reproduce the reference full QM calculations. For further applications of mEFP, we recommend either the scheme in which the peptide is cut along the Cα-C bond, giving rise to one fragment per amino acid, or the scheme with two cuts per amino acid, along the Cα-C and Cα-N bonds. While using these fragmentation schemes, the errors in solvatochromic shifts in electronic energy differences (excitation, ionization, electron detachment, or electron-attachment) do not exceed 0.1 eV. The largest error of QM/mEFP against QM/EFP (no fragmentation of the EFP part) is 0.06 eV (in most cases, the errors are 0.01-0.02 eV). The errors in the QM/molecular mechanics calculations with standard point charges can be as large as 0.3 eV.

[1]  A. Krylov,et al.  First-principle protocol for calculating ionization energies and redox potentials of solvated molecules and ions: theory and application to aqueous phenol and phenolate. , 2012, The journal of physical chemistry. B.

[2]  A. Lever,et al.  Vertical Ionization Energies and Electron Affinities of Ions in Solution from Outer-Sphere Charge Transfer Transition Energies. , 1998, Inorganic chemistry.

[3]  N. Reuter,et al.  Green fluorescent proteins: Empirical force field for the neutral and deprotonated forms of the chromophore. Molecular dynamics simulation's of the wild type and S65T mutant , 2002 .

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

[5]  X. H. Chen,et al.  Molecular caps for full quantum mechanical computation of peptide–water interaction energy , 2003, J. Comput. Chem..

[6]  B. Grigorenko,et al.  Modeling of biomolecular systems with the quantum mechanical and molecular mechanical method based on the effective fragment potential technique: Proposal of flexible fragments , 2002 .

[7]  Joost VandeVondele,et al.  Calculation of redox properties: understanding short- and long-range effects in rubredoxin. , 2007, The journal of physical chemistry. B.

[8]  Mark S Gordon,et al.  Fully Integrated Effective Fragment Molecular Orbital Method. , 2013, Journal of chemical theory and computation.

[9]  Mark S. Gordon,et al.  Damping functions in the effective fragment potential method , 2009 .

[10]  C. Cramer,et al.  Computational electrochemistry: prediction of liquid-phase reduction potentials. , 2014, Physical chemistry chemical physics : PCCP.

[11]  Hai Lin,et al.  QM/MM: What Have We Learned, Where Are We, and Where Do We Go from Here? , 2007 .

[12]  Benjamin T. Miller,et al.  A parallel implementation of the analytic nuclear gradient for time-dependent density functional theory within the Tamm–Dancoff approximation , 1999 .

[13]  Donald G Truhlar,et al.  Perspective on Foundations of Solvation Modeling: The Electrostatic Contribution to the Free Energy of Solvation. , 2008, Journal of chemical theory and computation.

[14]  Kenneth M. Merz,et al.  Combined Quantum Mechanical/Molecular Mechanical Methodologies Applied to Biomolecular Systems , 1999 .

[15]  Anna I Krylov,et al.  Quantum chemistry behind bioimaging: insights from ab initio studies of fluorescent proteins and their chromophores. , 2012, Accounts of chemical research.

[16]  Andreas Dreuw,et al.  Single-reference ab initio methods for the calculation of excited states of large molecules. , 2005, Chemical reviews.

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

[18]  John M Herbert,et al.  Aiming for benchmark accuracy with the many-body expansion. , 2014, Accounts of chemical research.

[19]  Ilya Kaliman,et al.  Hybrid MPI/OpenMP parallelization of the effective fragment potential method in the libefp software library , 2015, J. Comput. Chem..

[20]  Jan H. Jensen,et al.  Mapping Enzymatic Catalysis Using the Effective Fragment Molecular Orbital Method: Towards all ab initio Biochemistry , 2012, PloS one.

[21]  Anna I Krylov,et al.  On the Nature of an Extended Stokes Shift in the mPlum Fluorescent Protein. , 2015, The journal of physical chemistry. B.

[22]  A. Kolomeisky,et al.  Turning On and Off Photoinduced Electron Transfer in Fluorescent Proteins by π-Stacking, Halide Binding, and Tyr145 Mutations. , 2016, Journal of the American Chemical Society.

[23]  M. Gordon,et al.  Modeling Solvent Effects on Electronic Excited States , 2011 .

[24]  Trygve Helgaker,et al.  Molecular Electronic-Structure Theory: Helgaker/Molecular Electronic-Structure Theory , 2000 .

[25]  Mark S. Gordon,et al.  Electrostatic energy in the effective fragment potential method: Theory and application to benzene dimer , 2007, J. Comput. Chem..

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

[27]  Evgeny Epifanovsky,et al.  New implementation of high‐level correlated methods using a general block tensor library for high‐performance electronic structure calculations , 2013, J. Comput. Chem..

[28]  Yihan Shao,et al.  YinYang atom: a simple combined ab initio quantum mechanical molecular mechanical model. , 2007, The journal of physical chemistry. A.

[29]  the Carbonium Ion in the , 1976 .

[30]  Michael D. Hands,et al.  Interactions between halide anions and a molecular hydrophobic interface. , 2013, Faraday discussions.

[31]  Alán Aspuru-Guzik,et al.  Advances in molecular quantum chemistry contained in the Q-Chem 4 program package , 2014, Molecular Physics.

[32]  Jacopo Tomasi,et al.  Quantum Mechanical Continuum Solvation Models , 2005 .

[33]  Richard A. Friesner,et al.  A mixed quantum mechanics/molecular mechanics (QM/MM) method for large-scale modeling of chemistry in protein environments , 2000, J. Comput. Chem..

[34]  Jacob Kongsted,et al.  Scrutinizing the effects of polarization in QM/MM excited state calculations. , 2011, Physical chemistry chemical physics : PCCP.

[35]  L. Slipchenko.  Effective Fragment Potential Method , 2016 .

[36]  Shridhar R Gadre,et al.  Molecular tailoring approach: a route for ab initio treatment of large clusters. , 2014, Accounts of chemical research.

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

[38]  Mark S. Gordon,et al.  General atomic and molecular electronic structure system , 1993, J. Comput. Chem..

[39]  A. Krylov,et al.  Toward understanding the redox properties of model chromophores from the green fluorescent protein family: an interplay between conjugation, resonance stabilization, and solvent effects. , 2012, The journal of physical chemistry. B.

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

[41]  Jan H. Jensen,et al.  Accurate Intraprotein Electrostatics Derived from First Principles: An Effective Fragment Potential Method Study of the Proton Affinities of Lysine 55 and Tyrosine 20 in Turkey Ovomucoid Third Domain , 2001 .

[42]  T. Crawford,et al.  Electronic structure of the two isomers of the anionic form of p-coumaric acid chromophore. , 2011, The Journal of chemical physics.

[43]  Ilya Kaliman,et al.  LIBEFP: A new parallel implementation of the effective fragment potential method as a portable software library , 2013, J. Comput. Chem..

[44]  Jacob Kongsted,et al.  The polarizable embedding coupled cluster method. , 2011, Journal of Chemical Physics.

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

[46]  M. Gordon,et al.  Accurate first principles model potentials for intermolecular interactions. , 2013, Annual review of physical chemistry.

[47]  Evgeny Epifanovsky,et al.  General implementation of the resolution-of-the-identity and Cholesky representations of electron repulsion integrals within coupled-cluster and equation-of-motion methods: theory and benchmarks. , 2013, The Journal of chemical physics.

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

[49]  Martin Head-Gordon,et al.  Scaled second-order perturbation corrections to configuration interaction singles: efficient and reliable excitation energy methods. , 2007, The journal of physical chemistry. A.

[50]  Mark S. Gordon,et al.  An effective fragment method for modeling solvent effects in quantum mechanical calculations , 1996 .

[51]  Mark S. Gordon,et al.  Solvent-induced frequency shifts: configuration interaction singles combined with the effective fragment potential method. , 2010, The journal of physical chemistry. A.

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

[53]  J. Herbert,et al.  Accurate and efficient quantum chemistry calculations for noncovalent interactions in many-body systems: the XSAPT family of methods. , 2015, The journal of physical chemistry. A.

[54]  C. Cramer,et al.  Reply to Comment on “A Universal Approach to Solvation Modeling” , 2009 .

[55]  J. Herbert,et al.  Experimental benchmark data and systematic evaluation of two a posteriori, polarizable-continuum corrections for vertical excitation energies in solution. , 2015, The journal of physical chemistry. A.

[56]  Hui Li,et al.  Intraprotein electrostatics derived from first principles: Divide‐and‐conquer approaches for QM/MM calculations , 2003, J. Comput. Chem..

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

[58]  Jan H. Jensen,et al.  The Effective Fragment Molecular Orbital Method for Fragments Connected by Covalent Bonds , 2012, PloS one.

[59]  Pär Söderhjelm,et al.  How accurate can a force field become? A polarizable multipole model combined with fragment-wise quantum-mechanical calculations. , 2009, The journal of physical chemistry. A.

[60]  L. Slipchenko,et al.  Solvent effects on the electronic transitions of p-nitroaniline: a QM/EFP study. , 2011, The journal of physical chemistry. A.

[61]  A. Krylov,et al.  Complex-scaled equation-of-motion coupled-cluster method with single and double substitutions for autoionizing excited states: theory, implementation, and examples. , 2013, The Journal of chemical physics.

[62]  Michael W. Schmidt,et al.  Noncovalent interactions in extended systems described by the effective fragment potential method: theory and application to nucleobase oligomers. , 2010, The journal of physical chemistry. A.

[63]  L. Slipchenko,et al.  Solvation of the excited states of chromophores in polarizable environment: orbital relaxation versus polarization. , 2010, The journal of physical chemistry. A.

[64]  Brian W Matthews,et al.  Use of stabilizing mutations to engineer a charged group within a ligand-binding hydrophobic cavity in T4 lysozyme. , 2009, Biochemistry.

[65]  Anna I. Krylov,et al.  Q‐Chem: an engine for innovation , 2013 .

[66]  Evgeny Epifanovsky,et al.  Quantum Chemical Benchmark Studies of the Electronic Properties of the Green Fluorescent Protein Chromophore. 1. Electronically Excited and Ionized States of the Anionic Chromophore in the Gas Phase. , 2009, Journal of chemical theory and computation.

[67]  Anna I. Krylov,et al.  Chromophore Photoreduction in Red Fluorescent Proteins Is Responsible for Bleaching and Phototoxicity , 2014, The journal of physical chemistry. B.

[68]  Rodney J. Bartlett,et al.  Coupled‐cluster theory and its equation‐of‐motion extensions , 2012 .

[69]  Iris Antes,et al.  Adjusted Connection Atoms for Combined Quantum Mechanical and Molecular Mechanical Methods , 1999 .

[70]  I. Polyakov,et al.  Toward Molecular-Level Characterization of Photoinduced Decarboxylation of the Green Fluorescent Protein: Accessibility of the Charge-Transfer States. , 2012, Journal of chemical theory and computation.

[71]  Jacob Kongsted,et al.  Excited States in Solution through Polarizable Embedding , 2010 .

[72]  M. Head‐Gordon,et al.  Long-range corrected hybrid density functionals with damped atom-atom dispersion corrections. , 2008, Physical chemistry chemical physics : PCCP.

[73]  Olexandr Isayev,et al.  Effect of solvation on the vertical ionization energy of thymine: from microhydration to bulk. , 2011, The journal of physical chemistry. A.

[74]  A. Krylov,et al.  Effect of protein environment on electronically excited and ionized states of the green fluorescent protein chromophore. , 2011, The journal of physical chemistry. B.

[75]  R. Ballard,et al.  Photoelectron spectra of anions in solutions and the ionisation energy of the solvated electron , 1983 .

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

[77]  Spencer R Pruitt,et al.  Efficient and accurate fragmentation methods. , 2014, Accounts of chemical research.