scholarly journals Thermodynamics and kinetics of the amyloid-β peptide revealed by Markov state models based on MD data in agreement with experiment

2021 ◽  
Author(s):  
Arghadwip Paul ◽  
Suman Samantray ◽  
Marco Anteghini ◽  
Mohammed Khaled ◽  
Birgit Strodel
Keyword(s):  
Amyloid Β ◽  
Md Simulations ◽  
Amyloid Β Peptide ◽  
Intrinsically Disordered ◽  
Thermodynamics And Kinetics ◽  
Markov State ◽  
Markov State Models ◽  
State Models ◽  
Kinetics Of

The convergence of MD simulations is tested using varying measures for the intrinsically disordered amyloid-β peptide (Aβ). Markov state models show that 20–30 μs of MD is needed to reliably reproduce the thermodynamics and kinetics of Aβ.

scholarly journals Thermodynamics and kinetics of the amyloid-β peptide revealed by Markov state models based on MD data in agreement with experiment

2020 ◽  
Author(s):  
Arghadwip Paul ◽  
Suman Samantray ◽  
Marco Anteghini ◽  
Birgit Strodel
Keyword(s):  
Alzheimer’S Disease ◽  
Force Fields ◽  
Amyloid Β ◽  
Md Simulations ◽  
Amyloid Β Peptide ◽  
Thermodynamics And Kinetics ◽  
Markov State ◽  
Markov State Models ◽  
State Models ◽  
Kinetics Of

AbstractThe amlyoid-β peptide (Aβ) is closely linked to the development of Alzheimer’s disease. Molecular dynamics (MD) simulations have become an indispensable tool for studying the behavior of this peptide at the (sub)molecular level, thereby providing insight into the molecular basis of Alzheimer’s disease. General key aspects of MD simulations are the force field used for modeling the peptide or protein and its environment, which is important for accurate modeling of the system of interest, and the length of the simulations, which determines whether or not equilibrium is reached. In this study we address these points by analyzing 30-µs MD simulations acquired for Aβ40 using seven different force fields. We assess the convergence of these simulations based on the convergence of various structural properties and of NMR and fluorescence spectroscopic observables. Moreover, we calculate Markov state models for each of the seven MD simulations, which provide an unprecedented view of the thermodynamics and kinetics of the amyloid-β peptide. This further allows us to provide answers for pertinent questions, like: Which force fields are suitable for modeling Aβ? (a99SB-UCB and a99SB-ILDN/TIP4P-D); What does Aβ peptide really look like? (mostly extended and disordered) and; How long does it take MD simulations of Aβ to attain equilibrium? (20–30 µs). We believe the analyses presented in this study will provide a useful reference guide for important questions relating to the structure and dynamics of Aβin particular, and by extension other similar disordered peptides.

scholarly journals Predicting the Kinetics of RNA Oligonucleotides Using Markov State Models

2017 ◽  
Vol 13 (2) ◽  
pp. 926-934 ◽  
Author(s):  
Giovanni Pinamonti ◽  
Jianbo Zhao ◽  
David E. Condon ◽  
Fabian Paul ◽  
Frank Noè ◽  
...  
Keyword(s):  
Markov State ◽  
Markov State Models ◽  
State Models ◽  
Kinetics Of

scholarly journals Markov models for the elucidation of allosteric regulation

2018 ◽  
Vol 373 (1749) ◽  
pp. 20170178 ◽  
Author(s):  
Ushnish Sengupta ◽  
Birgit Strodel
Keyword(s):  
Markov Models ◽  
Allosteric Regulation ◽  
Md Simulations ◽  
Markov Modelling ◽  
Conformational Ensembles ◽  
Markov State ◽  
Protein Allostery ◽  
Markov State Models ◽  
State Models ◽  
Allosteric Mechanisms

Allosteric regulation refers to the process where the effect of binding of a ligand at one site of a protein is transmitted to another, often distant, functional site. In recent years, it has been demonstrated that allosteric mechanisms can be understood by the conformational ensembles of a protein. Molecular dynamics (MD) simulations are often used for the study of protein allostery as they provide an atomistic view of the dynamics of a protein. However, given the wealth of detailed information hidden in MD data, one has to apply a method that allows extraction of the conformational ensembles underlying allosteric regulation from these data. Markov state models are one of the most promising methods for this purpose. We provide a short introduction to the theory of Markov state models and review their application to various examples of protein allostery studied by MD simulations. We also include a discussion of studies where Markov modelling has been employed to analyse experimental data on allosteric regulation. We conclude our review by advertising the wider application of Markov state models to elucidate allosteric mechanisms, especially since in recent years it has become straightforward to construct such models thanks to software programs like PyEMMA and MSMBuilder. This article is part of a discussion meeting issue ‘Allostery and molecular machines’.

MasterMSM: A Package for Constructing Master Equation Models of Molecular Dynamics

2019 ◽  
Author(s):  
David De Sancho ◽  
Anne Aguirre
Keyword(s):  
Molecular Dynamics ◽  
Master Equation ◽  
Md Simulations ◽  
Distinctive Features ◽  
Markov State Models ◽  
Software Packages ◽  
Long Time ◽  
State Models ◽  
Definition Of ◽  
Kinetics Of

<div>Markov state models (MSMs) have become one of the most important techniques for understanding biomolecular transitions from classical molecular dynamics (MD) simulations. MSMs provide a systematized way of accessing the long time kinetics of the system of interest from the short-timescale microscopic transitions observed in simulation trajectories. At the same time, they provide a consistent description of the equilibrium and dynamical properties of the system of interest, and they are ideally suited for comparisons against experiment. A few software packages exist for building MSMs, which have been widely adopted. Here we introduce MasterMSM, a new Python package that uses the master equation formulation of MSMs and provides a number of new algorithms for building and analyzing these models. We demonstrate some of the most distinctive features of the package, including the estimation of rates, definition of core-sets for transition based assignment of states, the estimation of committors and fluxes, and the sensitivity analysis of the emerging networks. The package is available at https://github.com/daviddesancho/MasterMSM.</div>

scholarly journals Molecular dynamics simulations of protein aggregation: protocols for simulation setup and analysis with Markov state models and transition networks

2020 ◽  
Author(s):  
Suman Samantray ◽  
Wibke Schumann ◽  
Alexander-Maurice Illig ◽  
Martin Carballo-Pacheco ◽  
Arghadwip Paul ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Protein Aggregation ◽  
Amyloid Fibrils ◽  
Md Simulations ◽  
Aggregation Process ◽  
Markov State ◽  
Markov State Models ◽  
State Models ◽  
Dynamics Simulations ◽  
Transition Networks

AbstractProtein disorder and aggregation play significant roles in the pathogenesis of numerous neuro-degenerative diseases, such as Alzheimer’s and Parkinson’s disease. The end products of the aggregation process in these diseases are β-sheet rich amyloid fibrils. Though in most cases small, soluble oligomers formed during amyloid aggregation are the toxic species. A full understanding of the physicochemical forces behind the protein aggregation process is required if one aims to reveal the molecular basis of the various amyloid diseases. Among a multitude of biophysical and biochemical techniques that are employed for studying protein aggregation, molecular dynamics (MD) simulations at the atomic level provide the highest temporal and spatial resolution of this process, capturing key steps during the formation of amyloid oligomers. Here we provide a step-by-step guide for setting up, running, and analyzing MD simulations of aggregating peptides using GROMACS. For the analysis we provide the scripts that were developed in our lab, which allow to determine the oligomer size and inter-peptide contacts that drive the aggregation process. Moreover, we explain and provide the tools to derive Markov state models and transition networks from MD data of peptide aggregation.

MasterMSM: A Package for Constructing Master Equation Models of Molecular Dynamics

2019 ◽  
Author(s):  
David De Sancho ◽  
Anne Aguirre
Keyword(s):  
Molecular Dynamics ◽  
Master Equation ◽  
Md Simulations ◽  
Distinctive Features ◽  
Markov State Models ◽  
Software Packages ◽  
Long Time ◽  
State Models ◽  
Definition Of ◽  
Kinetics Of

<div>Markov state models (MSMs) have become one of the most important techniques for understanding biomolecular transitions from classical molecular dynamics (MD) simulations. MSMs provide a systematized way of accessing the long time kinetics of the system of interest from the short-timescale microscopic transitions observed in simulation trajectories. At the same time, they provide a consistent description of the equilibrium and dynamical properties of the system of interest, and they are ideally suited for comparisons against experiment. A few software packages exist for building MSMs, which have been widely adopted. Here we introduce MasterMSM, a new Python package that uses the master equation formulation of MSMs and provides a number of new algorithms for building and analyzing these models. We demonstrate some of the most distinctive features of the package, including the estimation of rates, definition of core-sets for transition based assignment of states, the estimation of committors and fluxes, and the sensitivity analysis of the emerging networks. The package is available at https://github.com/daviddesancho/MasterMSM.</div>

scholarly journals Improved Kinetics of Molecular Simulations using Biased Markov State Models

2016 ◽  
Vol 110 (3) ◽  
pp. 523a
Author(s):  
Joseph F. Rudzinski ◽  
Kurt Kremer ◽  
Tristan Bereau
Keyword(s):  
Molecular Simulations ◽  
Markov State ◽  
Markov State Models ◽  
State Models ◽  
Kinetics Of

scholarly journals Laplacian score and genetic algorithm based automatic feature selection for Markov State Models in adaptive sampling based molecular dynamics

2020 ◽  
Vol 2 ◽  
pp. e9
Author(s):  
Anu George ◽  
Madhura Purnaprajna ◽  
Prashanth Athri
Keyword(s):  
Molecular Dynamics ◽  
Adaptive Sampling ◽  
Simulated Data ◽  
Md Simulations ◽  
Considerable Effect ◽  
Feature Subset ◽  
Markov State ◽  
Markov State Models ◽  
State Models ◽  
Laplacian Score

Adaptive sampling molecular dynamics based on Markov State Models use short parallel MD simulations to accelerate simulations, and are proven to identify hidden conformers. The accuracy of the predictions provided by it depends on the features extracted from the simulated data that is used to construct it. The identification of the most important features in the trajectories of the simulated system has a considerable effect on the results. Methods In this study, we use a combination of Laplacian scoring and genetic algorithms to obtain an optimized feature subset for the construction of the MSM. The approach is validated on simulations of three protein folding complexes, and two protein ligand binding complexes. Results Our experiments show that this approach produces better results when the number of samples is significantly lesser than the number of features extracted. We also observed that this method mitigates over fitting that occurs due to high dimensionality of large biosystems with shorter simulation times.

Sign in / Sign up

Export Citation Format

Share Document