Identifying the leading dynamics of ubiquitin: a comparison between the tICA and the LE4PD slow fluctuations in amino acids' position
arXiv:2112.13171 · doi:10.1063/5.0059688
Abstract
Molecular Dynamics (MD) simulations of proteins implicitly contain the information connecting the atomistic molecular structure and proteins' biologically relevant motion, where large-scale fluctuations are deemed to guide folding and function. In the complex multiscale processes described by MD trajectories it is difficult to identify, separate, and study those large-scale fluctuations. This problem can be formulated as the need to identify a small number of collective variables that guide the slow kinetic processes. Among the methods used to study the slow, leading processes in proteins' dynamics, the time-lagged independent component analysis, or tICA, has been extensively used. Recently, we developed a Langevin coarse-grained approach for the dynamics of proteins, called the Langevin Equation for Protein Dynamics or LE4PD. This approach partitions the protein's MD dynamics into uncorrelated, wavelength-dependent, diffusive modes, and it associates to each mode a free-energy map. In the free energy maps, we measure the spatial extension and the time evolution of the mode-dependent, slow dynamical fluctuations, using the string method and a Markov state model. The theory identifies the slow collective variables in the rescaled LE4PD normal modes. Here, we compare the tICA modes' predictions with the collective LE4PD modes. We observe that the two methods consistently identify the nature and extension of the slowest fluctuation processes. The tICA separates the slow, leading processes in a smaller number of modes than the LE4PD does. However, LE4PD provides time-dependent information and a formal connection to the physics of the kinetic processes missing in the pure statistical analysis of tICA.
arXiv comment20 pages main text; 11 pages SM; 28 total figures; the following article has been accepted by The Journal of Chemical Physics. After it is published, it will be found at https://doi.org/10.1063/5.0059688
References in corpus (16)
- String Method for the Study of Rare Events
- Identification of slow molecular order parameters for Markov model construction
- Induced fit, conformational selection and independent dynamic segments: an extended view of binding events
- Diffusive Dynamics of the Reaction Coordinate for Protein Folding Funnels
- Time-lagged autoencoders: Deep learning of slow collective variables for molecular kinetics
- Caliber based spectral gap optimization of order parameters (SGOOP) for sampling complex molecular systems
- Data-driven model reduction and transfer operator approximation
- Nonlinear Discovery of Slow Molecular Modes using State-Free Reversible VAMPnets
- A variational conformational dynamics approach to the selection of collective variables in metadynamics
- A First Principle Approach to Rescale the Dynamics of Simulated Coarse-Grained Macromolecular Liquids
- Variational Selection of Features for Molecular Kinetics
- Variational approach for learning Markov processes from time series data
- Relaxation Mode Analysis and Markov State Relaxation Mode Analysis for Chignolin in Aqueous Solution near a Transition Temperature
- Capabilities and Limitations of Time-lagged Autoencoders for Slow Mode Discovery in Dynamical Systems
- Predicting protein dynamics from structural ensembles
- Identifying the leading dynamics of ubiquitin: a comparison between the tICA and the LE4PD slow fluctuations in amino acids' position