Rapid QM/MM approach for biomolecular systems under periodic boundary conditions: Combination of the density‐functional tight‐binding theory and particle mesh Ewald method

A quantum mechanical/molecular mechanical (QM/MM) approach based on the density‐functional tight‐binding (DFTB) theory is a useful tool for analyzing chemical reaction systems in detail. In this study, an efficient QM/MM method is developed by the combination of the DFTB/MM and particle mesh Ewald (PME) methods. Because the Fock matrix, which is required in the DFTB calculation, is analytically obtained by the PME method, the Coulomb energy is accurately and rapidly computed. For assessing the performance of this method, DFTB/MM calculations and molecular dynamics simulation are conducted for a system consisting of two amyloid‐β(1‐16) peptides and a zinc ion in explicit water under periodic boundary conditions. As compared with that of the conventional Ewald summation method, the computational cost of the Coulomb energy by utilizing the present approach is drastically reduced, i.e., 166.5 times faster. Furthermore, the deviation of the electronic energy is less than 10−6 Eh . © 2016 Wiley Periodicals, Inc.

[1]  Seifert,et al.  Construction of tight-binding-like potentials on the basis of density-functional theory: Application to carbon. , 1995, Physical review. B, Condensed matter.

[2]  Nathalie Reuter,et al.  Frontier Bonds in QM/MM Methods: A Comparison of Different Approaches , 2000 .

[3]  H. Naiki,et al.  Ultrasonication-dependent production and breakdown lead to minimum-sized amyloid fibrils , 2009, Proceedings of the National Academy of Sciences.

[4]  G. Seifert,et al.  Calculations of molecules, clusters, and solids with a simplified LCAO-DFT-LDA scheme , 1996 .

[5]  Michael Gaus,et al.  Parametrization of DFTB3/3OB for Magnesium and Zinc for Chemical and Biological Applications , 2014, The journal of physical chemistry. B.

[6]  C. Masters,et al.  Alzheimer's Disease Amyloid-β Binds Copper and Zinc to Generate an Allosterically Ordered Membrane-penetrating Structure Containing Superoxide Dismutase-like Subunits* , 2001, The Journal of Biological Chemistry.

[7]  J. Herbert,et al.  Periodic boundary conditions for QM/MM calculations: Ewald summation for extended Gaussian basis sets. , 2013, The Journal of chemical physics.

[8]  J. Sipe,et al.  Review: history of the amyloid fibril. , 2000, Journal of structural biology.

[9]  Yuko Okamoto,et al.  Thermodynamic perspective on the dock-lock growth mechanism of amyloid fibrils. , 2009, The journal of physical chemistry. B.

[10]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[11]  Michael Gaus,et al.  DFTB3: Extension of the self-consistent-charge density-functional tight-binding method (SCC-DFTB). , 2011, Journal of chemical theory and computation.

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

[13]  M. Shiga,et al.  Ab initio quantum mechanical/molecular mechanical molecular dynamics using multiple-time-scale approach and perturbation theory , 2007 .

[14]  S. Ten-no,et al.  A combined quantum mechanical and molecular mechanical method using modified generalized hybrid orbitals: implementation for electronic excited states. , 2011, Physical chemistry chemical physics : PCCP.

[15]  Yuko Okamoto,et al.  Explicit symplectic integrators of molecular dynamics algorithms for rigid-body molecules in the canonical, isobaric-isothermal, and related ensembles. , 2007, The Journal of chemical physics.

[16]  H. Naiki,et al.  Ultrasonication-induced Amyloid Fibril Formation of β2-Microglobulin* , 2005, Journal of Biological Chemistry.

[17]  A Mitsutake,et al.  Generalized-ensemble algorithms for molecular simulations of biopolymers. , 2000, Biopolymers.

[18]  P. P. Ewald.  Die Berechnung optischer und elektrostatischer Gitterpotentiale , 1921 .

[19]  J. Lee,et al.  Zn(2+) effect on structure and residual hydrophobicity of amyloid β-peptide monomers. , 2014, The journal of physical chemistry. B.

[20]  C. Dobson,et al.  Protein misfolding, functional amyloid, and human disease. , 2006, Annual review of biochemistry.

[21]  Omar Al-Ubaydli,et al.  The Causal Effect of Market Priming on Trust: An Experimental Investigation Using Randomized Control , 2013, PloS one.

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

[23]  Hisashi Okumura,et al.  Dimerization process of amyloid-β(29-42) studied by the Hamiltonian replica-permutation molecular dynamics simulations. , 2014, The journal of physical chemistry. B.

[24]  M. Elstner,et al.  Parametrization and Benchmark of DFTB3 for Organic Molecules. , 2013, Journal of chemical theory and computation.

[25]  Colin L Masters,et al.  Pleomorphic copper coordination by Alzheimer's disease amyloid-beta peptide. , 2009, Journal of the American Chemical Society.

[26]  Christopher M. Dobson,et al.  Molecular recycling within amyloid fibrils , 2005, Nature.

[27]  Yuko Okamoto,et al.  Generalized-ensemble algorithms for molecular dynamics simulations , 2007 .

[28]  Masato Kobayashi,et al.  Three pillars for achieving quantum mechanical molecular dynamics simulations of huge systems: Divide‐and‐conquer, density‐functional tight‐binding, and massively parallel computation , 2016, J. Comput. Chem..

[29]  H. Senn,et al.  QM/MM Methods for Biological Systems , 2006 .

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

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

[32]  Adrian J Mulholland,et al.  Analysis of chorismate mutase catalysis by QM/MM modelling of enzyme-catalysed and uncatalysed reactions. , 2011, Organic & biomolecular chemistry.

[33]  H. Takeuchi,et al.  Metal binding modes of Alzheimer's amyloid beta-peptide in insoluble aggregates and soluble complexes. , 2000, Biochemistry.

[34]  P. Kolandaivel,et al.  Role of zinc and copper metal ions in amyloid β-peptides Aβ1–40 and Aβ1–42 aggregation , 2014 .

[35]  M. Levitt,et al.  Theoretical studies of enzymic reactions: dielectric, electrostatic and steric stabilization of the carbonium ion in the reaction of lysozyme. , 1976, Journal of molecular biology.

[36]  Y. Shigeta.  Hybrid QM/MM studies on energetics of malonaldehyde in condensed phase , 2004 .

[37]  Hisashi Okumura,et al.  Amyloid fibril disruption by ultrasonic cavitation: nonequilibrium molecular dynamics simulations. , 2014, Journal of the American Chemical Society.

[38]  K. Morokuma,et al.  Insights into carbon nanotube and graphene formation mechanisms from molecular simulations: a review , 2015, Reports on progress in physics. Physical Society.

[39]  D Thirumalai,et al.  Monomer adds to preformed structured oligomers of Aβ-peptides by a two-stage dock–lock mechanism , 2007, Proceedings of the National Academy of Sciences.

[40]  L. Ying,et al.  Introduction of a fluorescent probe to amyloid-β to reveal kinetic insights into its interactions with copper(II). , 2015, Angewandte Chemie.

[41]  J. Rumfeldt,et al.  Sonication of proteins causes formation of aggregates that resemble amyloid , 2004, Protein science : a publication of the Protein Society.

[42]  Sándor Suhai,et al.  Self-consistent-charge density-functional tight-binding method for simulations of complex materials properties , 1998 .

[43]  Hisashi Okumura,et al.  Comparison of Replica-Permutation Molecular Dynamics Simulations with and without Detailed Balance Condition , 2015 .

[44]  Jean-Didier Maréchal,et al.  Three dimensional models of Cu(2+)-Aβ(1-16) complexes from computational approaches. , 2011, Journal of the American Chemical Society.

[45]  Kwangho Nam,et al.  Acceleration of Ab Initio QM/MM Calculations under Periodic Boundary Conditions by Multiscale and Multiple Time Step Approaches. , 2014, Journal of chemical theory and computation.

[46]  Fabrizio Chiti,et al.  Amyloid formation by globular proteins under native conditions. , 2009, Nature chemical biology.

[47]  S. Furlan,et al.  Modeling of the Zn2+ binding in the 1-16 region of the amyloid beta peptide involved in Alzheimer's disease. , 2009, Physical chemistry chemical physics : PCCP.

[48]  R. Riek,et al.  3D structure of Alzheimer's amyloid-β(1–42) fibrils , 2005 .

[49]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[50]  Hisashi Okumura,et al.  Temperature and pressure denaturation of chignolin: Folding and unfolding simulation by multibaric‐multithermal molecular dynamics method , 2012, Proteins.