Optimized parameter selection reveals trends in Markov state models for protein folding

Abstract

As molecular dynamics simulations access increasingly longer time scales, complementary advances in the analysis of biomolecular time-series data are necessary. Markov state models offer a powerful framework for this analysis by describing a system's states and the transitions between them. A recently established variational theorem for Markov state models now enables modelers to systematically determine the best way to describe a system's dynamics. In the context of the variational theorem, we analyze ultra-long folding simulations for a canonical set of twelve proteins [K. Lindorff-Larsen et al., Science 334, 517 (2011)] by creating and evaluating many types of Markov state models. We present a set of guidelines for constructing Markov state models of protein folding; namely, we recommend the use of cross-validation and a kinetically motivated dimensionality reduction step for improved descriptions of folding dynamics. We also warn that precise kinetics predictions rely on the features chosen to describe the system and pose the description of kinetic uncertainty across ensembles of models as an open issue. Published by AIP Publishing. [http://dx

Specific PIP2 Binding Promotes Calcium Activation of TMEM16A Chloride Channels

“…69,70,71 The simulation trajectories are clustered into micro-states using a hybrid k-centers k-medoids algorithm 68 . The generalized matrix Rayleigh quotient (GMRQ) method 71,72 is applied to optimize the hyperparameters as described in previous studies 73,74 .…”

The Dynamic Conformational Landscapes of the Protein Methyltransferase SETD8

“…We conducted approximately 500×1 µs explicit-solvent MD simulations from each seed and accumulated 5 milliseconds of aggregate data in 10 million conformational snapshots for apo-SETD8 (Figures S16, Table S3). To identify functionally relevant conformational states and their transitions, we built MSMs using a pipeline that employs machine learning and extensive hyperparameter optimization to identify slow degrees of freedom and structural and kinetic criteria to cluster conformational snapshots into discrete conformational states (Figures S17-24, Tables S4, S5) 38 . This approach identified 24 kinetically metastable conformations (macrostates) from an optimized, cross-validated set of 100 microstates (Figures 2c, S25-30, Tables S6, S7).…”

“…To determine the optimal number of microstates, we again used variational scoring [94][95][96][97] combined with cross-validation 38 to evaluate model quality. The full dataset (4.931 ms, 0.5 ns/frame, 9,862,657 frames) was separately featurized with the top-scoring feature sets: 6,567 distances (featurization a above) and 920 dihedral angles (featurization c above).…”

“…The code used for the generation and analysis of the molecular dynamics data is available via a Github repository at https://github.com/choderalab/SETD8-materials.Optimal hyperparameter selection for featurization and tICA. To select the optimal featurization of the data for subsequent Markov state model (MSM) analysis, we used variational scoring[94][95][96][97] combined with cross-validation 38 to evaluate model quality, consistent with modern MSM construction practice38 .…”

Modelling Intrinsically Disordered Protein Dynamics as Networks of Transient Secondary Structure

“…All-atom molecular dynamics (MD) simulations have the potential to offer atomistic-level detail in the dynamics of IDPs; however, modelling the dynamics of a protein using MD simulations requires describing the conformation space of the protein in a way that lends insight into states and processes of interest to the researcher. Significant study has gone into how to partition this phase space in the field of modelling the dynamics of structured proteins (3), and use of methods for featurization and time-lagged dimensionality reduction (6,7) have greatly improved the modelling process (8) for proteins that inhabit a few structured conformations. However, in practice, when we use these same tools to analyze MD simulation data of IDPs, we find that these methods face challenges when used to analyze the conformational landscape of an IDP.…”