Markov State Models Provide Insights into Dynamic Modulation of Protein Function

Abstract

ConspectusProtein function is inextricably linked to protein dynamics. As we move from a static structural picture to a dynamic ensemble view of protein structure and function, novel computational paradigms are required for observing and understanding conformational dynamics of proteins and its functional implications. In principle, molecular dynamics simulations can provide the time evolution of atomistic models of proteins, but the long time scales associated with functional dynamics make it difficult to observe rare dynamical transitions. The issue of extracting essential functional components of protein dynamics from noisy simulation data presents another set of challenges in obtaining an unbiased understanding of protein motions. Therefore, a methodology that provides a statistical framework for efficient sampling and a human-readable view of the key aspects of functional dynamics from data analysis is required. The Markov state model (MSM), which has recently become popular worldwide for studying protein dynamics, is an example of such a framework.In this Account, we review the use of Markov state models for efficient sampling of the hierarchy of time scales associated with protein dynamics, automatic identification of key conformational states, and the degrees of freedom associated with slow dynamical processes. Applications of MSMs for studying long time scale phenomena such as activation mechanisms of cellular signaling proteins has yielded novel insights into protein function. In particular, from MSMs built using large-scale simulations of GPCRs and kinases, we have shown that complex conformational changes in proteins can be described in terms of structural changes in key structural motifs or “molecular switches” within the protein, the transitions between functionally active and inactive states of proteins proceed via multiple pathways, and ligand or substrate binding modulates the flux through these pathways. Finally, MSMs also provide a theoretical toolbox for studying the effect of nonequilibrium perturbations on conformational dynamics. Considering that protein dynamics in vivo occur under nonequilibrium conditions, MSMs coupled with nonequilibrium statistical mechanics provide a way to connect cellular components to their functional environments. Nonequilibrium perturbations of protein folding MSMs reveal the presence of dynamically frozen glass-like states in their conformational landscape. These frozen states are also observed to be rich in β-sheets, which indicates their possible role in the nucleation of β-sheet rich aggregates such as those observed in amyloid-fibril formation. Finally, we describe how MSMs have been used to understand the dynamical behavior of intrinsically disordered proteins such as amyloid-β, human islet amyloid polypeptide, and p53. While certainly not a panacea for studying functional dynamics, MSMs provide a rigorous theoretical foundation for understanding complex entropically dominated processes and a convenient lens for viewing protein motions.

The effects of temperature on aerobic metabolism: towards a mechanistic understanding of the responses of ectotherms to a changing environment

“…However, both the equilibrium and Sharpe-Schoolfield models may represent over-simplifications of the highly dynamic structures of proteins (Shukla et al, 2015). Rather than occurring in a small number of conformations (e.g.…”

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

“…[12][13][14]16,[19][20][21][24][25][26][27][28][29][30][31][32][33][34][35][36][37][38] Building a MSM involves decomposing the phase space sampled by one or more MD trajectories into a set of discrete states and estimating the (conditional) transition probabilities between each pair of states. However, there are many ways to perform the state decomposition, and the choice of MSM building protocol can introduce subjectivity into the analysis.…”

Recognition Dynamics of p53 and MDM2: Implications for Peptide Design

“…Any kinetic information has been gleaned from experimental techniques such as SPR which provide the rates of association and dissociation. However, details of these processes have only been available from a few Brownian dynamics simulations and more recently long timescale atomistic simulations coupled to markov chain analyses 32,33 . In a recent Brownian dynamics study, based on kinetic arguments, we showed that residues on the MDM2 surface outside the p53 binding pocket could stabilize the p53-MDM2 encounter complex, a state that is formed prior to the formation of the final bound complex, by increasing p53 residence times and allowing for allosteric interactions with MDM2 34 .…”