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β.
Protein disorder and aggregation play significant roles in the pathogenesis of numerous neurodegenerative diseases, such as Alzheimer's and Parkinson's disease. The end products of the aggregation process in these diseases are highly structured amyloid fibrils. Though in most cases small, soluble oligomers formed during amyloid aggregation are the toxic species. A full understanding of the physicochemical forces that drive protein aggregation is thus required if one aims for the rational design of drugs targeting the formation of amyloid oligomers. 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.
The 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.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.