On the Application of Molecular-Dynamics Based Markov State Models to Functional Proteins.
Ontology highlight
ABSTRACT:
SUBMITTER: Malmstrom RD
PROVIDER: S-EPMC4248791 | biostudies-other | 2014 Jul
REPOSITORIES: biostudies-other
Similar Datasets
Project description:Ecological systems can often be characterised by changes among a finite set of underlying states pertaining to individuals, populations, communities or entire ecosystems through time. Owing to the inherent difficulty of empirical field studies, ecological state dynamics operating at any level of this hierarchy can often be unobservable or 'hidden'. Ecologists must therefore often contend with incomplete or indirect observations that are somehow related to these underlying processes. By formally disentangling state and observation processes based on simple yet powerful mathematical properties that can be used to describe many ecological phenomena, hidden Markov models (HMMs) can facilitate inferences about complex system state dynamics that might otherwise be intractable. However, HMMs have only recently begun to gain traction within the broader ecological community. We provide a gentle introduction to HMMs, establish some common terminology, review the immense scope of HMMs for applied ecological research and provide a tutorial on implementation and interpretation. By illustrating how practitioners can use a simple conceptual template to customise HMMs for their specific systems of interest, revealing methodological links between existing applications, and highlighting some practical considerations and limitations of these approaches, our goal is to help establish HMMs as a fundamental inferential tool for ecologists.
| S-EPMC7702077 | biostudies-literature
Project description:Mutations in leucine-rich repeat kinase 2 (LRRK2) are recognized as the most frequent cause of Parkinson's disease (PD). As a multidomain ROCO protein, LRRK2 is characterized by the presence of both a Ras-of-complex (ROC) GTPase domain and a kinase domain connected through the C-terminal of an ROC domain (COR). The bienzymatic ROC-COR-kinase catalytic triad indicated the potential role of GTPase domain in regulating kinase activity. However, as a functional GTPase, the detailed intrinsic regulation of the ROC activation cycle remains poorly understood. Here, combining extensive molecular dynamics simulations and Markov state models, we disclosed the dynamic structural rearrangement of ROC's homodimer during nucleotide turnover. Our study revealed the coupling between dimerization extent and nucleotide-binding state, indicating a nucleotide-dependent dimerization-based activation scheme adopted by ROC GTPase. Furthermore, inspired by the well-known R1441C/G/H PD-relevant mutations within the ROC domain, we illuminated the potential allosteric molecular mechanism for its pathogenetic effects through enabling faster interconversion between inactive and active states, thus trapping ROC in a prolonged activated state, while the implicated allostery could provide further guidance for identification of regulatory allosteric pockets on the ROC complex. Our investigations illuminated the thermodynamics and kinetics of ROC homodimer during nucleotide-dependent activation for the first time and provided guidance for further exploiting ROC as therapeutic targets for controlling LRRK2 functionality in PD treatment.
| S-EPMC8467336 | biostudies-literature
Project description:Understanding protein conformational variability remains a challenge in drug discovery. The issue arises in protein kinases, whose multiple conformational states can affect the binding of small-molecule inhibitors. To overcome this challenge, we propose a comprehensive computational framework based on Markov state models (MSMs). Our framework integrates the information from explicit-solvent molecular dynamics simulations to accurately rank-order the accessible conformational variants of a target protein. We tested the methodology using Abl kinase with a reference and blind-test set. Only half of the Abl conformational variants discovered by our approach are present in the disclosed X-ray structures. The approach successfully identified a protein conformational state not previously observed in public structures but evident in a retrospective analysis of Lilly in-house structures: the X-ray structure of Abl with WHI-P154. Using a MSM-derived model, the free energy landscape and kinetic profile of Abl was analyzed in detail highlighting opportunities for targeting the unique metastable states.
| S-EPMC6317529 | biostudies-literature
Project description:α-Synuclein (αS) is a presynaptic protein that binds to cell membranes and is linked to Parkinson's disease (PD). Binding of αS to membranes is a likely first step in the molecular pathophysiology of PD. The αS molecule can adopt multiple conformations, being largely disordered in water, adopting a β-sheet conformation when present in amyloid fibrils, and forming a dynamic multiplicity of α-helical conformations when bound to lipid bilayers and related membrane-mimetic surfaces. Multiscale molecular dynamics simulations in conjunction with nuclear magnetic resonance (NMR) and cross-linking mass spectrometry (XLMS) measurements are used to explore the interactions of αS with an anionic lipid bilayer. The simulations and NMR measurements together reveal a break in the helical structure of the central non-amyloid-β component (NAC) region of αS in the vicinity of residues 65-70, which may facilitate subsequent oligomer formation. Coarse-grained simulations of αS starting from the structure of αS when bound to a detergent micelle reveal the overall pattern of protein contacts to anionic lipid bilayers, while subsequent all-atom simulations provide details of conformational changes upon membrane binding. In particular, simulations and NMR data for liposome-bound αS indicate incipient β-strand formation in the NAC region, which is supported by intramolecular contacts seen via XLMS and simulations. Markov state models based on the all-atom simulations suggest a mechanism of conformational change of membrane-bound αS via a dynamic helix break in the region of residue 65 in the NAC region. The emergent dynamic model of membrane-interacting αS advances our understanding of the mechanism of PD, potentially aiding the design of novel therapeutic approaches.
| S-EPMC8006134 | biostudies-literature
Project description:An enhanced bioinformatics tool incorporating the participation of molecular structure as well as sequence in protein DNA recognition is proposed and tested. Boltzmann probability models of sequence-dependent DNA structure from all-atom molecular dynamics simulations were obtained and incorporated into hidden Markov models (HMMs) that can recognize molecular structural signals as well as sequence in protein-DNA binding sites on a genome. The binding of catabolite activator protein (CAP) to cognate DNA sequences was used as a prototype case for implementation and testing of the method. The results indicate that even HMMs based on probabilistic roll/tilt dinucleotide models of sequence-dependent DNA structure have some capability to discriminate between known CAP binding and nonbinding sites and to predict putative CAP binding sites in unknowns. Restricting HMMs to sequence only in regions of strong consensus in which the protein makes base specific contacts with the cognate DNA further improved the discriminatory capabilities of the HMMs. Comparison of results with controls based on sequence only indicates that extending the definition of consensus from sequence to structure improves the transferability of the HMMs, and provides further supportive evidence of a role for dynamical molecular structure as well as sequence in genomic regulatory mechanisms.
| S-EPMC124344 | biostudies-literature
Project description:Atomistic molecular dynamics simulations of the TZ1 beta-hairpin peptide have been carried out using an implicit model for the solvent. The trajectories have been analyzed using a Markov state model defined on the projections along two significant observables and a kinetic network approach. The Markov state model allowed for an unbiased identification of the metastable states of the system, and provided the basis for commitment probability calculations performed on the kinetic network. The kinetic network analysis served to extract the main transition state for folding of the peptide and to validate the results from the Markov state analysis. The combination of the two techniques allowed for a consistent and concise characterization of the dynamics of the peptide. The slowest relaxation process identified is the exchange between variably folded and denatured species, and the second slowest process is the exchange between two different subsets of the denatured state which could not be otherwise identified by simple inspection of the projected trajectory. The third slowest process is the exchange between a fully native and a partially folded intermediate state characterized by a native turn with a proximal backbone H-bond, and frayed side-chain packing and termini. The transition state for the main folding reaction is similar to the intermediate state, although a more native like side-chain packing is observed.
| S-EPMC3106446 | biostudies-literature
Project description:Phosphorylation of tyrosine 32 in K-Ras has been shown to influence function by disrupting the GTPase cycle. To shed light on the underlying mechanism and atomic basis of this process, we carried out a comparative investigation of the oncogenic G12D K-Ras mutant and its phosphorylated variant (pTyr32) using all-atom molecular dynamics simulations and Markov state models. We show that, despite sharing a number of common features, G12D and pTyr32-G12D K-Ras exhibit some distinct conformational states and fluctuations. In addition to notable differences in conformation and dynamics of residues surrounding the GTP binding site, nonlocal changes were observed at a number of loops. Switch I is more flexible in pTyr32-G12D K-Ras while switch II is more flexible in G12D K-Ras. We also used time-lagged independent component analysis and k-means clustering to identify five metastable states for each system. We utilized transition path theory to calculate the transition probabilities for each state to build a Markov state model for each system. These models and other close inspections suggest that the phosphorylation of Tyr32 strongly affects protein dynamics and the active site conformation, especially with regards to the canonical switch conformations and dynamics.
| S-EPMC7020251 | biostudies-literature
Project description:The extent to which glass-like kinetics govern dynamics in protein folding has been heavily debated. Here, we address the subject with an application of space-time perturbation theory to the dynamics of protein folding Markov state models. Borrowing techniques from the s-ensemble method, we argue that distinct active and inactive phases exist for protein folding dynamics, and that kinetics for specific systems can fall into either dynamical regime. We do not, however, observe a true glass transition in any system studied. We go on to discuss how these inactive and active phases might relate to general protein folding properties.
| S-EPMC3677858 | biostudies-literature
Project description:The study of functional brain networks has grown rapidly over the past decade. While most functional connectivity (FC) analyses estimate one static network structure for the entire length of the functional magnetic resonance imaging (fMRI) time series, recently there has been increased interest in studying time-varying changes in FC. Hidden Markov models (HMMs) have proven to be a useful modeling approach for discovering repeating graphs of interacting brain regions (brain states). However, a limitation lies in HMMs assuming that the sojourn time, the number of consecutive time points in a state, is geometrically distributed. This may encourage inaccurate estimation of the time spent in a state before switching to another state. We propose a hidden semi-Markov model (HSMM) approach for inferring time-varying brain networks from fMRI data, which explicitly models the sojourn distribution. Specifically, we propose using HSMMs to find each subject's most probable series of network states and the graphs associated with each state, while properly estimating and modeling the sojourn distribution for each state. We perform a simulation study, as well as an analysis on both task-based fMRI data from an anxiety-inducing experiment and resting-state fMRI data from the Human Connectome Project. Our results demonstrate the importance of model choice when estimating sojourn times and reveal their potential for understanding healthy and diseased brain mechanisms.
| S-EPMC6504179 | biostudies-literature
Project description:Recent studies have proposed that one can summarize brain activity into dynamics among a relatively small number of hidden states and that such an approach is a promising tool for revealing brain function. Hidden Markov models (HMMs) are a prevalent approach to inferring such neural dynamics among discrete brain states. However, the impact of assuming Markovian structure in neural time series data has not been sufficiently examined. Here, to address this situation and examine the performance of the HMM, we compare the model with the Gaussian mixture model (GMM), which is with no temporal regularization and thus a statistically simpler model than the HMM, by applying both models to synthetic time series generated from empirical resting-state functional magnetic resonance imaging (fMRI) data. We compared the GMM and HMM for various sampling frequencies, lengths of recording per participant, numbers of participants and numbers of independent component signals. We find that the HMM attains a better accuracy of estimating the hidden state than the GMM in a majority of cases. However, we also find that the accuracy of the GMM is comparable to that of the HMM under the condition that the sampling frequency is reasonably low (e.g., TR = 2.88 or 3.60 s) or the data are relatively short. These results suggest that the GMM can be a viable alternative to the HMM for investigating hidden-state dynamics under this condition.
| S-EPMC9291560 | biostudies-literature