Back

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

Samantray, S.; Schumann, W.; Illig, A.-M.; Paul, A.; Barz, B.; Strodel, B.

2020-04-27 biophysics
10.1101/2020.04.25.060269 bioRxiv
Show abstract

Protein disorder and aggregation play significant roles in the pathogenesis of numerous neuro-degenerative diseases, such as Alzheimers and Parkinsons disease. The end products of the aggregation process in these diseases are {beta}-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.

Matching journals

The top 7 journals account for 50% of the predicted probability mass.

50% of probability mass above

"Similar papers" are the closest papers from that journal in the model's embedding space. They show what the match is built on, but the ranking comes mostly from a classifier over the whole training set, not from these examples alone.