Fully analytic energy gradient in the fragment molecular orbital method.

The Z-vector equations are derived and implemented for solving the response term due to the external electrostatic potentials, and the corresponding contribution is added to the energy gradients in the framework of the fragment molecular orbital (FMO) method. To practically solve the equations for large molecules like proteins, the equations are decoupled by taking advantage of the local nature of fragments in the FMO method and establishing the self-consistent Z-vector method. The resulting gradients are compared with numerical gradients for the test molecular systems: (H(2)O)(64), alanine decamer, hydrated chignolin with the protein data bank (PDB) ID of 1UAO, and a Trp-cage miniprotein construct (PDB ID: 1L2Y). The computation time for calculating the response contribution is comparable to or less than that of the FMO self-consistent charge calculation. It is also shown that the energy gradients for the electrostatic dimer approximation are fully analytic, which significantly reduces the computational costs. The fully analytic FMO gradient is parallelized with an efficiency of about 98% on 32 nodes.

[1]  Kazuo Kitaura,et al.  On the accuracy of the 3-body fragment molecular orbital method (FMO) applied to density functional theory , 2004 .

[2]  Kazuo Kitaura,et al.  Second order Møller-Plesset perturbation theory based upon the fragment molecular orbital method. , 2004, The Journal of chemical physics.

[3]  U. Nagashima,et al.  Development of an ab initio MO-MD program based on fragment MO method – an attempt to analyze the fluctuation of protein , 2004 .

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

[5]  Jan H. Jensen,et al.  Analytic gradient for the adaptive frozen orbital bond detachment in the fragment molecular orbital method , 2009 .

[6]  Kazuo Kitaura,et al.  The three-body fragment molecular orbital method for accurate calculations of large systems , 2006 .

[7]  Yuichi Inadomi,et al.  Fragment molecular orbital method: application to molecular dynamics simulation, ‘ab initio FMO-MD’ , 2003 .

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

[9]  Shigenori Tanaka,et al.  Ab initio Path Integral Molecular Dynamics Based on Fragment Molecular Orbital Method , 2009 .

[10]  A. Imamura,et al.  Electronic structures and nonlinear optical properties of supramolecular associations of benzo-2,1,3-chalcogendiazoles by the elongation method , 2010 .

[11]  Takeshi Ishikawa,et al.  Fragment Molecular Orbital method‐based Molecular Dynamics (FMO‐MD) as a simulator for chemical reactions in explicit solvation , 2009, J. Comput. Chem..

[12]  Mark S. Gordon,et al.  Accurate methods for large molecular systems. , 2009, The journal of physical chemistry. B.

[13]  Takeshi Ishikawa,et al.  How does an S(N)2 reaction take place in solution? Full ab initio MD simulations for the hydrolysis of the methyl diazonium ion. , 2008, Journal of the American Chemical Society.

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

[15]  Thom Vreven,et al.  QM:QM electronic embedding using Mulliken atomic charges: energies and analytic gradients in an ONIOM framework. , 2008, The Journal of chemical physics.

[16]  Kenneth M. Merz,et al.  Importance of dispersion and electron correlation in ab initio protein folding. , 2009, The journal of physical chemistry. B.

[17]  Alexander D. MacKerell.  Empirical force fields for biological macromolecules: Overview and issues , 2004, J. Comput. Chem..

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

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

[20]  Kazuo Kitaura,et al.  A combined effective fragment potential-fragment molecular orbital method. II. Analytic gradient and application to the geometry optimization of solvated tetraglycine and chignolin. , 2011, The Journal of chemical physics.

[21]  Noriyuki Kurita,et al.  Fragment molecular orbital method with density functional theory and DIIS convergence acceleration , 2003 .

[22]  K. Ohno,et al.  A self-consistent molecular field theory for aggregates of neutral molecules. I , 1972 .

[23]  Peter A. Kollman,et al.  AMBER, a package of computer programs for applying molecular mechanics, normal mode analysis, molecular dynamics and free energy calculations to simulate the structural and energetic properties of molecules , 1995 .

[24]  Masato Kobayashi,et al.  Time-dependent Hartree–Fock frequency-dependent polarizability calculation applied to divide-and-conquer electronic structure method , 2010 .

[25]  Donald G Truhlar,et al.  X-Pol Potential: An Electronic Structure-Based Force Field for Molecular Dynamics Simulation of a Solvated Protein in Water. , 2009, Journal of chemical theory and computation.

[26]  Y. Mochizuki,et al.  Theoretical study of hydration models of trivalent rare-earth ions using model core potentials , 2010 .

[27]  K. Kitaura,et al.  Fragment-Molecular-Orbital-Method-Based ab Initio NMR Chemical-Shift Calculations for Large Molecular Systems , 2010 .

[28]  Yuichi Inadomi,et al.  Fragment molecular orbital study of the electronic excitations in the photosynthetic reaction center of Blastochloris viridis , 2009, J. Comput. Chem..

[29]  Kazuo Kitaura,et al.  Polarizable continuum model with the fragment molecular orbital‐based time‐dependent density functional theory , 2008, J. Comput. Chem..

[30]  Stefan Goedecker,et al.  Linear scaling electronic structure methods in chemistry and physics , 2003, Comput. Sci. Eng..

[31]  D. Truhlar,et al.  Quantum mechanical methods for enzyme kinetics. , 2003, Annual review of physical chemistry.

[32]  T. Nakano,et al.  Does amination of formaldehyde proceed through a zwitterionic intermediate in water? Fragment molecular orbital molecular dynamics simulations by using constraint dynamics. , 2010, Chemistry.

[33]  Takeshi Ishikawa,et al.  A fully quantum mechanical simulation study on the lowest n-π* state of hydrated formaldehyde , 2007 .

[34]  Kwang S. Kim,et al.  Theory and applications of computational chemistry : the first forty years , 2005 .

[35]  P. Kollman,et al.  Calculating structures and free energies of complex molecules: combining molecular mechanics and continuum models. , 2000, Accounts of chemical research.

[36]  Toshio Fujita,et al.  Novel Quantitative Structure-Activity Studies of HIV-1 Protease Inhibitors of the Cyclic Urea Type Using Descriptors Derived from Molecular Dynamics and Molecular Orbital Calculations , 2009 .

[37]  Kaori Fukuzawa,et al.  Large-scale FMO-MP3 calculations on the surface proteins of influenza virus, hemagglutinin (HA) and neuraminidase (NA) , 2010 .

[38]  Ye Mei,et al.  QUANTUM CALCULATION OF PROTEIN SOLVATION AND PROTEIN–LIGAND BINDING FREE ENERGY FOR HIV-1 PROTEASE/WATER COMPLEX , 2009 .

[39]  P. Kollman,et al.  Biomolecular simulations: recent developments in force fields, simulations of enzyme catalysis, protein-ligand, protein-protein, and protein-nucleic acid noncovalent interactions. , 2001, Annual review of biophysics and biomolecular structure.

[40]  G. Scuseria,et al.  Ab initio molecular dynamics: Propagating the density matrix with Gaussian orbitals , 2001 .

[41]  Mark S. Gordon,et al.  A new hierarchical parallelization scheme: Generalized distributed data interface (GDDI), and an application to the fragment molecular orbital method (FMO) , 2004, J. Comput. Chem..

[42]  Takeshi Ishikawa,et al.  Fragment molecular orbital calculations on red fluorescent proteins (DsRed and mFruits). , 2009, The journal of physical chemistry. B.

[43]  Hermann Stoll,et al.  A Simple One-Body Approach to the Calculation of the First Electronic Absorption Band of Water. , 2009, Journal of chemical theory and computation.

[44]  H. Kubota,et al.  Ab initio Study of the Relationship between the Pressure and the Change in the Dipole Moment of a Molecule in a Ferroelectric Crystal , 1994 .

[45]  Jiali Gao,et al.  A molecular-orbital derived polarization potential for liquid water , 1998 .

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

[47]  L. Seijo,et al.  Ab initio model potential embedded cluster calculations including lattice relaxation and polarization: Local distortions on Mn2+‐doped CaF2 , 1995 .

[48]  Keiji Morokuma,et al.  Molecular Orbital Studies of Hydrogen Bonds. III. C=O···H–O Hydrogen Bond in H2CO···H2O and H2CO···2H2O , 1971 .

[49]  Kazuo Kitaura,et al.  Molecular recognition mechanism of FK506 binding protein: An all‐electron fragment molecular orbital study , 2007, Proteins.

[50]  Kazuo Kitaura,et al.  The Fragment Molecular Orbital Method: Practical Applications to Large Molecular Systems , 2009 .

[51]  Lou Massa,et al.  Calculation of strong and weak interactions in TDA1 and RangDP52 by the kernel energy method , 2009, Proceedings of the National Academy of Sciences.

[52]  K. Raghavachari,et al.  ONIOM-based QM:QM electronic embedding method using Löwdin atomic charges: energies and analytic gradients. , 2010, The Journal of chemical physics.

[53]  M. Gordon,et al.  A combined effective fragment potential-fragment molecular orbital method. I. The energy expression and initial applications. , 2009, The Journal of chemical physics.

[54]  Hiroshi Chuman,et al.  Correlation Analyses on Binding Affinity of Substituted Benzenesulfonamides with Carbonic Anhydrase Using ab Initio MO Calculations on Their Complex Structures , 2010, J. Chem. Inf. Model..

[55]  Yuto Komeiji,et al.  Three-body expansion and generalized dynamic fragmentation improve the fragment molecular orbital-based molecular dynamics (FMO-MD)☆ , 2010 .

[56]  L. Seijo,et al.  The ab initio model potential representation of the crystalline environment. Theoretical study of the local distortion on NaCl:Cu+ , 1988 .

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

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

[59]  Masami Uebayasi,et al.  The fragment molecular orbital method for geometry optimizations of polypeptides and proteins. , 2007, The journal of physical chemistry. A.

[60]  Kazuo Kitaura,et al.  Coupled-cluster theory based upon the fragment molecular-orbital method. , 2005, The Journal of chemical physics.

[61]  Kazuo Kitaura,et al.  Derivatives of the approximated electrostatic potentials in the fragment molecular orbital method , 2009 .

[62]  D. Case,et al.  Generalized born models of macromolecular solvation effects. , 2000, Annual review of physical chemistry.

[63]  Kazuo Kitaura,et al.  Importance of the hybrid orbital operator derivative term for the energy gradient in the fragment molecular orbital method , 2010 .

[64]  Nuclear-Electronic Orbital Method within the Fragment Molecular Orbital Approach† , 2010 .

[65]  Junwei Zhang,et al.  VISCANA: Visualized Cluster Analysis of Protein-Ligand Interaction Based on the ab Initio Fragment Molecular Orbital Method for Virtual Ligand Screening , 2006, J. Chem. Inf. Model..

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

[67]  Yuichi Inadomi,et al.  DNA and estrogen receptor interaction revealed by fragment molecular orbital calculations. , 2007, The journal of physical chemistry. B.

[68]  Y. Yamaguchi,et al.  A New Dimension to Quantum Chemistry: Analytic Derivative Methods in AB Initio Molecular Electronic Structure Theory , 1994 .

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

[70]  Kazuo Kitaura,et al.  Time-dependent density functional theory based upon the fragment molecular orbital method. , 2007, The Journal of chemical physics.

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

[72]  K. Kistler,et al.  Solvatochromic shifts of uracil and cytosine using a combined multireference configuration interaction/molecular dynamics approach and the fragment molecular orbital method. , 2009, The journal of physical chemistry. A.

[73]  Spencer R Pruitt,et al.  Open-Shell Formulation of the Fragment Molecular Orbital Method. , 2010, Journal of chemical theory and computation.

[74]  Chakram S. Jayanthi,et al.  Order-/N methodologies and their applications , 2002 .

[75]  T. Nakano,et al.  Parallelized integral-direct CIS(D) calculations with multilayer fragment molecular orbital scheme , 2007 .

[76]  Yuto Komeiji,et al.  Fragment molecular orbital-based molecular dynamics (FMO-MD), a quantum simulation tool for large molecular systems , 2009 .

[77]  Jan H. Jensen,et al.  Effective fragment molecular orbital method: a merger of the effective fragment potential and fragment molecular orbital methods. , 2010, The journal of physical chemistry. A.

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

[79]  Mark S. Gordon,et al.  A parallel distributed data CPHF algorithm for analytic Hessians , 2007, J. Comput. Chem..

[80]  Efthimios Kaxiras,et al.  A QM/MM Implementation of the Self-Consistent Charge Density Functional Tight Binding (SCC-DFTB) Method , 2001 .

[81]  T. Nakano,et al.  Fragment molecular orbital-based molecular dynamics (FMO-MD) simulations on hydrated Zn(II) ion , 2010 .

[82]  Kenneth M Merz,et al.  Divide-and-Conquer Hartree-Fock Calculations on Proteins. , 2010, Journal of chemical theory and computation.

[83]  Yuji Mochizuki,et al.  Configuration interaction singles method with multilayer fragment molecular orbital scheme , 2005 .

[84]  J. Ladik,et al.  Investigation of the interaction between molecules at medium distances: I. SCF LCAO MO supermolecule, perturbational and mutually consistent calculations for two interacting HF and CH2O molecules , 1975 .

[85]  Kazuo Kitaura,et al.  Multiconfiguration self-consistent-field theory based upon the fragment molecular orbital method. , 2005, The Journal of chemical physics.

[86]  Hui Li,et al.  The polarizable continuum model (PCM) interfaced with the fragment molecular orbital method (FMO) , 2006, J. Comput. Chem..

[87]  A. Warshel.  Computer simulations of enzyme catalysis: methods, progress, and insights. , 2003, Annual review of biophysics and biomolecular structure.

[88]  K. Kitaura,et al.  Time-dependent density functional theory with the multilayer fragment molecular orbital method , 2007 .

[89]  Hiroaki Tokiwa,et al.  Theoretical study of intramolecular interaction energies during dynamics simulations of oligopeptides by the fragment molecular orbital-Hamiltonian algorithm method. , 2005, The Journal of chemical physics.

[90]  Jules W. Moskowitz,et al.  Water Molecule Interactions , 1970 .

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

[92]  Kazuo Kitaura,et al.  All electron quantum chemical calculation of the entire enzyme system confirms a collective catalytic device in the chorismate mutase reaction. , 2006, The journal of physical chemistry. B.

[93]  Zhihai Liu,et al.  Evaluation of the performance of four molecular docking programs on a diverse set of protein‐ligand complexes , 2010, J. Comput. Chem..

[94]  Tohru Suzuki,et al.  Ab initio fragment molecular orbital studies of influenza virus hemagglutinin–sialosaccharide complexes toward chemical clarification about the virus host range determination , 2008, Glycoconjugate Journal.

[95]  Car,et al.  Unified approach for molecular dynamics and density-functional theory. , 1985, Physical review letters.