Molecular dynamics simulations of an α-synuclein NAC domain fragment with a ff14IDPSFF IDP-specific force field suggest β-sheet intermediate states of fibrillation

Abstract

For the discovery of treatments against synucleinopathies, it is necessary to unravel and fully understand the mechanism of fibrillation of proteins involved. Among them, α-synuclein (αS) plays a key role in the development of these diseases through its aggregation into oligomers found in Lewy bodies. However, its structural disorder as an intrinsically disordered protein (IDP) makes its characterization by experimental techniques arduously difficult. Atomistic simulations aim to provide insights into this blank canvas and, fortunately, some studies have already suggested promising mechanisms. Still, it is urgent to consider the IDP features in simulations, so recently a lot of force fields designed to deal with IDPs have been developed. In this study, we have carried out a total of 12 μs simulations of an αS core fragment using a popular ff14SB AMBER force field and the ff14IDPSFF variation that includes a grid-based energy correction map (CMAP) method. The predicted chemical shifts from the simulations and those measured from the αS protein in the NMR solution indicate that ff14IDPSFF reproduces the experimental data more accurately. Moreover, structural analysis exhibits opposite trends between secondary structure propensities. The ff14SB force field preserves the α-helices found in the micelle-bound αS structure, which is used as an initial conformation, while ff14IDPSFF stands out with increased structural disorder and the formation of β-sheets, which suggests that the IDP-specific force field can capture more suitable conformations representing the possible intermediate states of the fibrillation process.