Comparing normal modes across different models and scales: Hessian reduction versus coarse‐graining

Dimension reduction is often necessary when attempting to reach longer length and time scales in molecular simulations. It is realized by constraining degrees of freedom or by coarse‐graining the system. When evaluating the accuracy of a dimensional reduction, there is a practical challenge: the models yield vectors with different lengths, making a comparison by calculating their dot product impossible. This article investigates mapping procedures for normal mode analysis. We first review a horizontal mapping procedure for the reduced Hessian techniques, which projects out degrees of freedom. We then design a vertical mapping procedure for the “implosion” of the all‐atom (AA) Hessian to a coarse‐grained scale that is based upon vibrational subsystem analysis. This latter method derives both effective force constants and an effective kinetic tensor. Next, a series of metrics is presented for comparison across different scales, where special attention is given to proper mass‐weighting. The dimension‐dependent metrics, which require prior mapping for proper evaluation, are frequencies, overlap of normal mode vectors, probability similarity, Hessian similarity, collectivity of modes, and thermal fluctuations. The dimension‐independent metrics are shape derivatives, elastic modulus, vibrational free energy differences, heat capacity, and projection on a predefined basis set. The power of these metrics to distinguish between reasonable and unreasonable models is tested on a toy alpha helix system and a globular protein; both are represented at several scales: the AA scale, a Gō‐like model, a canonical elastic network model, and a network model with intentionally unphysical force constants. Published 2012 Wiley Periodicals, Inc.

[1]  L. Mouawad,et al.  Diagonalization in a mixed basis: A method to compute low‐frequency normal modes for large macromolecules , 1993 .

[2]  George N Phillips,et al.  Application of elastic network models to proteins in the crystalline state. , 2009, Biophysical journal.

[3]  L. Björck,et al.  Protein PAB, a mosaic albumin-binding bacterial protein representing the first contemporary example of module shuffling. , 1994, The Journal of biological chemistry.

[4]  A. Atilgan,et al.  Direct evaluation of thermal fluctuations in proteins using a single-parameter harmonic potential. , 1997, Folding & design.

[5]  Andrew C. Simmonett,et al.  Capture of hydroxymethylene and its fast disappearance through tunnelling , 2008, Nature.

[6]  Sebastian Doniach,et al.  Toward the mechanism of dynamical couplings and translocation in hepatitis C virus NS3 helicase using elastic network model , 2007, Proteins.

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

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

[9]  R Brüschweiler,et al.  Characterization of collective and anisotropic reorientational protein dynamics. , 2000, Physical review letters.

[10]  Michel Waroquier,et al.  Experimental and theoretical IR study of methanol and ethanol conversion over H-SAPO-34 , 2011 .

[11]  Feng Gai,et al.  Folding of a three-helix bundle at the folding speed limit , 2004 .

[12]  Wenjun Zheng,et al.  A unification of the elastic network model and the Gaussian network model for optimal description of protein conformational motions and fluctuations. , 2008, Biophysical journal.

[13]  M. Karplus,et al.  Normal modes for specific motions of macromolecules: application to the hinge-bending mode of lysozyme. , 1985, Proceedings of the National Academy of Sciences of the United States of America.

[14]  R. M. Stratt,et al.  THE INSTANTANEOUS NORMAL MODES OF LIQUIDS , 1995 .

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

[16]  R. Jernigan,et al.  Anisotropy of fluctuation dynamics of proteins with an elastic network model. , 2001, Biophysical journal.

[17]  H. Berendsen,et al.  ALGORITHMS FOR MACROMOLECULAR DYNAMICS AND CONSTRAINT DYNAMICS , 1977 .

[18]  Gerhard Hummer,et al.  Coordinate-dependent diffusion in protein folding , 2009, Proceedings of the National Academy of Sciences.

[19]  Michel Waroquier,et al.  Normal Mode Analysis in Zeolites: Toward an Efficient Calculation of Adsorption Entropies. , 2011, Journal of chemical theory and computation.

[20]  N. Go,et al.  Dynamics of a small globular protein in terms of low-frequency vibrational modes. , 1983, Proceedings of the National Academy of Sciences of the United States of America.

[21]  Bernard R. Brooks,et al.  CHARMMing: A New, Flexible Web Portal for CHARMM , 2008, J. Chem. Inf. Model..

[22]  Ramon Carbo,et al.  How similar is a molecule to another? An electron density measure of similarity between two molecular structures , 1980 .

[23]  H. Berendsen,et al.  Essential dynamics of proteins , 1993, Proteins.

[24]  Y. Sanejouand,et al.  Conformational change of proteins arising from normal mode calculations. , 2001, Protein engineering.

[25]  Markus Reiher,et al.  Localizing normal modes in large molecules. , 2009, The Journal of chemical physics.

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

[27]  Dusanka Janezic,et al.  Harmonic analysis of large systems. II. Comparison of different protein models , 1995, J. Comput. Chem..

[28]  Gregory A Voth,et al.  A multiscale coarse-graining method for biomolecular systems. , 2005, The journal of physical chemistry. B.

[29]  John D. Head,et al.  COMPUTATION OF VIBRATIONAL FREQUENCIES FOR ADSORBATES ON SURFACES , 1997 .

[30]  Sungsoo Na,et al.  Coarse‐graining of protein structures for the normal mode studies , 2007, J. Comput. Chem..

[31]  V Van Speybroeck,et al.  Normal modes for large molecules with arbitrary link constraints in the mobile block Hessian approach. , 2009, The Journal of chemical physics.

[32]  Changbong Hyeon,et al.  Theoretical perspectives on protein folding. , 2010, Annual review of biophysics.

[33]  B. Brooks,et al.  Probing the local dynamics of nucleotide-binding pocket coupled to the global dynamics: myosin versus kinesin. , 2005, Biophysical journal.

[34]  Jörg Grunenberg,et al.  Efficient computation of compliance matrices in redundant internal coordinates from Cartesian Hessians for nonstationary points , 2010 .

[35]  M. Karplus,et al.  Harmonic dynamics of proteins: normal modes and fluctuations in bovine pancreatic trypsin inhibitor. , 1983, Proceedings of the National Academy of Sciences of the United States of America.

[36]  Y. Sanejouand,et al.  Building‐block approach for determining low‐frequency normal modes of macromolecules , 2000, Proteins.

[37]  Alexander D. MacKerell,et al.  Extending the treatment of backbone energetics in protein force fields: Limitations of gas‐phase quantum mechanics in reproducing protein conformational distributions in molecular dynamics simulations , 2004, J. Comput. Chem..

[38]  Bernard R Brooks,et al.  Effects of denaturants and osmolytes on proteins are accurately predicted by the molecular transfer model , 2008, Proceedings of the National Academy of Sciences.

[39]  K. Hinsen,et al.  Harmonicity in slow protein dynamics , 2000 .

[40]  Michel Waroquier,et al.  First principle kinetic studies of zeolite-catalyzed methylation reactions. , 2011, Journal of the American Chemical Society.

[41]  Holger Gohlke,et al.  A natural coarse graining for simulating large biomolecular motion. , 2006, Biophysical journal.

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

[43]  T Verstraelen,et al.  Calculating Reaction Rates with Partial Hessians: Validation of the Mobile Block Hessian Approach. , 2008, Journal of chemical theory and computation.

[44]  N. Go,et al.  Studies on protein folding, unfolding, and fluctuations by computer simulation. II. A. Three‐dimensional lattice model of lysozyme , 1978 .

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

[46]  Gregory A Voth,et al.  Effective force fields for condensed phase systems from ab initio molecular dynamics simulation: a new method for force-matching. , 2004, The Journal of chemical physics.

[47]  Mai Suan Li,et al.  Biomolecules under mechanical force , 2010 .

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

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

[50]  L. Björck,et al.  Solution structure of the albumin-binding GA module: a versatile bacterial protein domain. , 1997, Journal of molecular biology.

[51]  S. Takada Protein folding simulation with solvent‐induced force field: Folding pathway ensemble of three‐helix‐bundle proteins , 2001, Proteins.

[52]  Asim Okur,et al.  MSCALE: A General Utility for Multiscale Modeling. , 2011, Journal of chemical theory and computation.

[53]  V Van Speybroeck,et al.  Mobile Block Hessian Approach with Adjoined Blocks: An Efficient Approach for the Calculation of Frequencies in Macromolecules. , 2009, Journal of chemical theory and computation.

[54]  Y. Sanejouand,et al.  Hinge‐bending motion in citrate synthase arising from normal mode calculations , 1995, Proteins.

[55]  T Verstraelen,et al.  Vibrational modes in partially optimized molecular systems. , 2007, The Journal of chemical physics.

[56]  Toon Verstraelen,et al.  TAMkin: A Versatile Package for Vibrational Analysis and Chemical Kinetics , 2010, J. Chem. Inf. Model..

[57]  B. Hess Convergence of sampling in protein simulations. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[58]  Hui Li,et al.  Partial Hessian vibrational analysis: the localization of the molecular vibrational energy and entropy , 2002 .

[59]  Michel Waroquier,et al.  Comparative study of various normal mode analysis techniques based on partial Hessians , 2009, J. Comput. Chem..

[60]  R. Jernigan,et al.  Residue-residue potentials with a favorable contact pair term and an unfavorable high packing density term, for simulation and threading. , 1996, Journal of molecular biology.

[61]  Steven A. Siegelbaum,et al.  Effects of Surface Water on Protein Dynamics Studied by a Novel Coarse-Grained Normal Mode Approach , 2008, Biophysical journal.

[62]  Wenjun Zheng,et al.  Approximate normal mode analysis based on vibrational subsystem analysis with high accuracy and efficiency. , 2009, The Journal of chemical physics.

[63]  Bernard R Brooks,et al.  Modeling protein conformational changes by iterative fitting of distance constraints using reoriented normal modes. , 2006, Biophysical journal.

[64]  Nicholas A Besley,et al.  Computation of the amide I band of polypeptides and proteins using a partial Hessian approach. , 2007, The Journal of chemical physics.

[65]  Martin Karplus,et al.  Molecular properties from combined QM/MM methods. I. Analytical second derivative and vibrational calculations , 2000 .

[66]  K. Hinsen Analysis of domain motions by approximate normal mode calculations , 1998, Proteins.

[67]  Markus Reiher,et al.  Convergence characteristics and efficiency of mode-tracking calculations on pre-selected molecular vibrations , 2004 .

[68]  I. Bahar,et al.  Normal mode analysis : theory and applications to biological and chemical systems , 2005 .

[69]  Dusanka Janezic,et al.  Harmonic analysis of large systems. I. Methodology , 1995, J. Comput. Chem..

[70]  Dusanka Janezic,et al.  Harmonic analysis of large systems. III. Comparison with molecular dynamics , 1995, J. Comput. Chem..

[71]  D. ben-Avraham,et al.  Dynamic and elastic properties of F-actin: a normal-modes analysis. , 1995, Biophysical journal.

[72]  Y. Sanejouand,et al.  A new approach for determining low‐frequency normal modes in macromolecules , 1994 .

[73]  M Waroquier,et al.  Cartesian formulation of the mobile block Hessian approach to vibrational analysis in partially optimized systems. , 2007, The Journal of chemical physics.

[74]  D. Thirumalai,et al.  LINKING RATES OF FOLDING IN LATTICE MODELS OF PROTEINS WITH UNDERLYING THERMODYNAMIC CHARACTERISTICS , 1998, cond-mat/9805061.

[75]  Valentina Tozzini,et al.  Coarse-grained models for proteins. , 2005, Current opinion in structural biology.

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