MSCALE: A General Utility for Multiscale Modeling.

The combination of theoretical models of macromolecules that exist at different spatial and temporal scales has become increasingly important for addressing complex biochemical problems. This work describes the extension of concurrent multiscale approaches, introduces a general framework for carrying out calculations, and describes its implementation into the CHARMM macromolecular modeling package. This functionality, termed MSCALE, generalizes both the additive and subtractive multiscale scheme (e.g. QM/MM ONIOM-type), and extends its support to classical force fields, coarse grained modeling (e.g. ENM, GNM, etc.), and a mixture of them all. The MSCALE scheme is completely parallelized with each subsystem running as an independent, but connected calculation. One of the most attractive features of MSCALE is the relative ease of implementation using the standard MPI communication protocol. This allows external access to the framework and facilitates the combination of functionality previously isolated in separate programs. This new facility is fully integrated with free energy perturbation methods, Hessian based methods, and the use of periodicity and symmetry, which allows the calculation of accurate pressures. We demonstrate the utility of this new technique with four examples; (1) subtractive QM/MM and QM/QM calculations; (2) multi-force field alchemical free energy perturbation; (3) integration with the SANDER module of AMBER and the TINKER package to gain access to potentials not available in CHARMM; and (4) mixed resolution (i.e. coarse grain / all-atom) normal mode analysis. The potential of this new tool is clearly established and in conclusion an interesting mathematical problem is highlighted and future improvements are proposed.

[1]  Ivet Bahar,et al.  Elastic network models for understanding biomolecular machinery: from enzymes to supramolecular assemblies , 2005, Physical biology.

[2]  Martin Karplus,et al.  Large amplitude conformational change in proteins explored with a plastic network model: adenylate kinase. , 2005, Journal of molecular biology.

[3]  T. Straatsma,et al.  Free energy of ionic hydration: Analysis of a thermodynamic integration technique to evaluate free energy differences by molecular dynamics simulations , 1988 .

[4]  P. S. Kim,et al.  A switch between two-, three-, and four-stranded coiled coils in GCN4 leucine zipper mutants. , 1993, Science.

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

[6]  Bernard R. Brooks,et al.  Exploring the quantum mechanical/molecular mechanical replica path method: a pathway optimization of the chorismate to prephenate Claisen rearrangement catalyzed by chorismate mutase , 2003 .

[7]  Bernard R Brooks,et al.  Normal-modes-based prediction of protein conformational changes guided by distance constraints. , 2005, Biophysical journal.

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

[9]  B. Brooks,et al.  Multiscale methods for macromolecular simulations. , 2008, Current opinion in structural biology.

[10]  I. Bahar,et al.  Coarse-grained normal mode analysis in structural biology. , 2005, Current opinion in structural biology.

[11]  B. Brooks,et al.  A super-linear minimization scheme for the nudged elastic band method , 2003 .

[12]  Pengyu Y. Ren,et al.  Polarizable Atomic Multipole Water Model for Molecular Mechanics Simulation , 2003 .

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

[14]  H. C. Andersen Rattle: A “velocity” version of the shake algorithm for molecular dynamics calculations , 1983 .

[15]  Valentina Tozzini,et al.  Multiscale modeling of proteins. , 2010, Accounts of chemical research.

[16]  Gregory A. Voth,et al.  The multiscale coarse-graining method. I. A rigorous bridge between atomistic and coarse-grained models. , 2008, The Journal of chemical physics.

[17]  Bernard R Brooks,et al.  Vibrational subsystem analysis: A method for probing free energies and correlations in the harmonic limit. , 2008, The Journal of chemical physics.

[18]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[19]  Stefan Bruckner,et al.  Unorthodox uses of Bennett's acceptance ratio method , 2009, J. Comput. Chem..

[20]  Shawn T. Brown,et al.  Advances in methods and algorithms in a modern quantum chemistry program package. , 2006, Physical chemistry chemical physics : PCCP.

[21]  D. Case,et al.  Exploring protein native states and large‐scale conformational changes with a modified generalized born model , 2004, Proteins.

[22]  Jay W. Ponder,et al.  Algorithms for calculating excluded volume and its derivatives as a function of molecular conformation and their use in energy minimization , 1991 .

[23]  Robert A Latour,et al.  The relationship between platelet adhesion on surfaces and the structure versus the amount of adsorbed fibrinogen. , 2010, Biomaterials.

[24]  Matthew L. Leininger,et al.  PSI3: An open‐source Ab Initio electronic structure package , 2007, J. Comput. Chem..

[25]  Ian J. Bush,et al.  The GAMESS-UK electronic structure package: algorithms, developments and applications , 2005 .

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

[27]  R. Larson,et al.  The MARTINI Coarse-Grained Force Field: Extension to Proteins. , 2008, Journal of chemical theory and computation.

[28]  Bernard R. Brooks,et al.  Interfacing Q‐Chem and CHARMM to perform QM/MM reaction path calculations , 2007, J. Comput. Chem..

[29]  Michel Waroquier,et al.  Efficient Calculation of QM/MM Frequencies with the Mobile Block Hessian. , 2011, Journal of chemical theory and computation.

[30]  Charles H. Bennett,et al.  Efficient estimation of free energy differences from Monte Carlo data , 1976 .

[31]  R. Jernigan,et al.  Structure-derived potentials and protein simulations. , 1996, Current opinion in structural biology.

[32]  Jianpeng Ma,et al.  CHARMM: The biomolecular simulation program , 2009, J. Comput. Chem..

[33]  Tjerk P. Straatsma,et al.  NWChem: A comprehensive and scalable open-source solution for large scale molecular simulations , 2010, Comput. Phys. Commun..

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

[35]  Markus Reiher,et al.  Atomistic Approaches in Modern Biology , 2007 .

[36]  D. Thirumalai,et al.  Low-frequency normal modes that describe allosteric transitions in biological nanomachines are robust to sequence variations , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[37]  Nadeem A. Vellore,et al.  Assessment of the transferability of a protein force field for the simulation of peptide-surface interactions. , 2010, Langmuir : the ACS journal of surfaces and colloids.

[38]  J. Onuchic,et al.  Protein folding mediated by solvation: Water expulsion and formation of the hydrophobic core occur after the structural collapse , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[39]  M. Karplus,et al.  A combined quantum mechanical and molecular mechanical potential for molecular dynamics simulations , 1990 .

[40]  Bernard R. Brooks,et al.  Optimization of quantum mechanical molecular mechanical partitioning schemes: Gaussian delocalization of molecular mechanical charges and the double link atom method , 2002 .

[41]  Jhih-Wei Chu,et al.  Emerging methods for multiscale simulation of biomolecular systems , 2007 .

[42]  Robert A Latour,et al.  Benchmark experimental data set and assessment of adsorption free energy for peptide-surface interactions. , 2009, Langmuir : the ACS journal of surfaces and colloids.

[43]  Cecilia Clementi,et al.  Balancing energy and entropy: a minimalist model for the characterization of protein folding landscapes. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[44]  Robert A Latour,et al.  Probing the conformation and orientation of adsorbed enzymes using side-chain modification. , 2009, Langmuir : the ACS journal of surfaces and colloids.

[45]  Shayantani Mukherjee,et al.  PRIMO/PRIMONA: A coarse‐grained model for proteins and nucleic acids that preserves near‐atomistic accuracy , 2010, Proteins.

[46]  Tirion,et al.  Large Amplitude Elastic Motions in Proteins from a Single-Parameter, Atomic Analysis. , 1996, Physical review letters.

[47]  William Gropp,et al.  MPI-2: Extending the Message-Passing Interface , 1996, Euro-Par, Vol. I.

[48]  M. Kim,et al.  A connection rule for alpha-carbon coarse-grained elastic network models using chemical bond information. , 2006, Journal of molecular graphics & modelling.

[49]  Eng‐Hui Yap,et al.  A coarse‐grained α‐carbon protein model with anisotropic hydrogen‐bonding , 2008 .

[50]  Thom Vreven,et al.  Combining Quantum Mechanics Methods with Molecular Mechanics Methods in ONIOM. , 2006, Journal of chemical theory and computation.

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

[52]  Bernard R Brooks,et al.  Langevin network model of myosin. , 2008, The journal of physical chemistry. B.

[53]  A. Liwo,et al.  Ab initio simulations of protein-folding pathways by molecular dynamics with the united-residue model of polypeptide chains. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[54]  Evert Jan Baerends,et al.  Calculation of pressure in case of periodic boundary conditions , 2006 .

[55]  Pengyu Y. Ren,et al.  Consistent treatment of inter‐ and intramolecular polarization in molecular mechanics calculations , 2002, J. Comput. Chem..

[56]  Bernard R Brooks,et al.  Exploring SCC-DFTB paths for mapping QM/MM reaction mechanisms. , 2007, The journal of physical chemistry. A.

[57]  M. Dewar,et al.  Ground States of Molecules. 38. The MNDO Method. Approximations and Parameters , 1977 .

[58]  Gregory A Voth,et al.  Coarse-grained free energy functions for studying protein conformational changes: a double-well network model. , 2007, Biophysical journal.

[59]  K. Morokuma,et al.  A NEW ONIOM IMPLEMENTATION IN GAUSSIAN98. PART I. THE CALCULATION OF ENERGIES, GRADIENTS, VIBRATIONAL FREQUENCIES AND ELECTRIC FIELD DERIVATIVES , 1999 .

[60]  J. Ponder,et al.  An efficient newton‐like method for molecular mechanics energy minimization of large molecules , 1987 .

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

[62]  K. Morokuma,et al.  ONIOM: A Multilayered Integrated MO + MM Method for Geometry Optimizations and Single Point Energy Predictions. A Test for Diels−Alder Reactions and Pt(P(t-Bu)3)2 + H2 Oxidative Addition , 1996 .

[63]  I. Bahar,et al.  Gaussian Dynamics of Folded Proteins , 1997 .

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

[65]  Walter Thiel,et al.  QM/MM methods for biomolecular systems. , 2009, Angewandte Chemie.

[66]  Daniel Russel,et al.  The structural dynamics of macromolecular processes. , 2009, Current opinion in cell biology.

[67]  Bernard R Brooks,et al.  Adjacent gauche stabilization in linear alkanes: implications for polymer models and conformational analysis. , 2005, The journal of physical chemistry. B.

[68]  Gregory A Voth,et al.  The multiscale coarse-graining method. II. Numerical implementation for coarse-grained molecular models. , 2008, The Journal of chemical physics.