Dynamic neutron scattering from conformational dynamics. II. Application using molecular dynamics simulation and Markov modeling.

Neutron scattering experiments directly probe the dynamics of complex molecules on the sub pico- to microsecond time scales. However, the assignment of the relaxations seen experimentally to specific structural rearrangements is difficult, since many of the underlying dynamical processes may exist on similar timescales. In an accompanying article, we present a theoretical approach to the analysis of molecular dynamics simulations with a Markov State Model (MSM) that permits the direct identification of structural transitions leading to each contributing relaxation process. Here, we demonstrate the use of the method by applying it to the configurational dynamics of the well-characterized alanine dipeptide. A practical procedure for deriving the MSM from an MD is introduced. The result is a 9-state MSM in the space of the backbone dihedral angles and the side-chain methyl group. The agreement between the quasielastic spectrum calculated directly from the atomic trajectories and that derived from the Markov state model is excellent. The dependence on the wavevector of the individual Markov processes is described. The procedure means that it is now practicable to interpret quasielastic scattering spectra in terms of well-defined intramolecular transitions with minimal a priori assumptions as to the nature of the dynamics taking place.

[1]  Bettina Keller,et al.  An Analysis of the Validity of Markov State Models for Emulating the Dynamics of Classical Molecular Systems and Ensembles. , 2011, Journal of chemical theory and computation.

[2]  Xuhui Huang,et al.  Using generalized ensemble simulations and Markov state models to identify conformational states. , 2009, Methods.

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

[4]  C. Schütte,et al.  Supplementary Information for “ Constructing the Equilibrium Ensemble of Folding Pathways from Short Off-Equilibrium Simulations ” , 2009 .

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

[6]  W. Kabsch,et al.  Dictionary of protein secondary structure: Pattern recognition of hydrogen‐bonded and geometrical features , 1983, Biopolymers.

[7]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[8]  D. Oesterhelt,et al.  Neutron scattering reveals extremely slow cell water in a Dead Sea organism , 2007, Proceedings of the National Academy of Sciences.

[9]  V. F. Sears Neutron scattering lengths and cross sections , 1992 .

[10]  F. Noé,et al.  Transition networks for modeling the kinetics of conformational change in macromolecules. , 2008, Current opinion in structural biology.

[11]  Frank Noé,et al.  Markov models and dynamical fingerprints: Unraveling the complexity of molecular kinetics , 2012 .

[12]  Jeremy C. Smith,et al.  Three classes of motion in the dynamic neutron-scattering susceptibility of a globular protein. , 2011, Physical review letters.

[13]  V. Pande,et al.  Error analysis and efficient sampling in Markovian state models for molecular dynamics. , 2005, The Journal of chemical physics.

[14]  E. Mamontov,et al.  Apparent Decoupling of the Dynamics of a Protein from the Dynamics of its Aqueous Solvent. , 2012, The journal of physical chemistry letters.

[15]  Jeremy C. Smith,et al.  Hierarchical analysis of conformational dynamics in biomolecules: transition networks of metastable states. , 2007, The Journal of chemical physics.

[16]  Jeremy C. Smith,et al.  Dynamical fingerprints for probing individual relaxation processes in biomolecular dynamics with simulations and kinetic experiments , 2011, Proceedings of the National Academy of Sciences.

[17]  Frank Noé,et al.  On the Approximation Quality of Markov State Models , 2010, Multiscale Model. Simul..

[18]  Vijay S. Pande,et al.  Markov State Models of Millisecond Folder ACBP Reveals New Views of the Folding Reaction , 2011 .

[19]  L. Hove Correlations in Space and Time and Born Approximation Scattering in Systems of Interacting Particles , 1954 .

[20]  K. Dill,et al.  Automatic discovery of metastable states for the construction of Markov models of macromolecular conformational dynamics. , 2007, The Journal of chemical physics.

[21]  Frank Noé,et al.  Dynamic neutron scattering from conformational dynamics. I. Theory and Markov models. , 2013, The Journal of chemical physics.

[22]  Thomas J Lane,et al.  MSMBuilder2: Modeling Conformational Dynamics at the Picosecond to Millisecond Scale. , 2011, Journal of chemical theory and computation.

[23]  P. Deuflhard,et al.  A Direct Approach to Conformational Dynamics Based on Hybrid Monte Carlo , 1999 .

[24]  Jeremy C. Smith,et al.  Sassena - X-ray and neutron scattering calculated from molecular dynamics trajectories using massively parallel computers , 2012, Comput. Phys. Commun..

[25]  William Swope,et al.  Describing Protein Folding Kinetics by Molecular Dynamics Simulations. 1. Theory , 2004 .

[26]  Frank Noé,et al.  Probing molecular kinetics with Markov models: metastable states, transition pathways and spectroscopic observables. , 2011, Physical chemistry chemical physics : PCCP.

[27]  J C Smith,et al.  Liquid-like side-chain dynamics in myoglobin. , 1994, Journal of molecular biology.

[28]  Vincent A. Voelz,et al.  Atomistic folding simulations of the five-helix bundle protein λ(6−85). , 2011, Journal of the American Chemical Society.

[29]  J. Smith,et al.  Protein dynamics: comparison of simulations with inelastic neutron scattering experiments , 1991, Quarterly Reviews of Biophysics.

[30]  Vijay S Pande,et al.  Progress and challenges in the automated construction of Markov state models for full protein systems. , 2009, The Journal of chemical physics.

[31]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

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

[33]  Frank Noé,et al.  Markov models of molecular kinetics: generation and validation. , 2011, The Journal of chemical physics.

[34]  A. Sokolov,et al.  Influence of hydration on the dynamics of lysozyme. , 2006, Biophysical journal.

[35]  Frank Noé,et al.  A Variational Approach to Modeling Slow Processes in Stochastic Dynamical Systems , 2012, Multiscale Model. Simul..

[36]  J. Swenson,et al.  Dynamics of a protein and its surrounding environment: a quasielastic neutron scattering study of myoglobin in water and glycerol mixtures. , 2009, The Journal of chemical physics.

[37]  Frank Noé,et al.  Kinetic characterization of the critical step in HIV-1 protease maturation , 2012, Proceedings of the National Academy of Sciences.

[38]  Vijay S. Pande,et al.  Everything you wanted to know about Markov State Models but were afraid to ask. , 2010, Methods.

[39]  Philip M. Long,et al.  Performance guarantees for hierarchical clustering , 2002, J. Comput. Syst. Sci..

[40]  Jeremy C. Smith,et al.  Picosecond timescale rigid‐helix and side‐chain motions in deoxymyoglobin , 1993, Proteins.

[41]  Frank Noé,et al.  EMMA: A Software Package for Markov Model Building and Analysis. , 2012, Journal of chemical theory and computation.