Computational investigation of copper-mediated conformational changes in α-synuclein dimer.

We report molecular dynamics simulation of dimers of α-synuclein, the peptide closely associated with onset of Parkinson's disease, both as metal-free dimer and with inter-chain bridging provided by Cu(II) ions. Our investigation reveals that the presence of copper-induced inter-chain bridging not only stabilizes α-synuclein dimers, but also leads to enhanced ß-sheet formation at critical regions within the N-terminal and NAC regions of the protein. These contacts are larger and longer-lived in the presence of copper, and as a result each peptide chain is more extended and less flexible than in the metal-free dimer. The persistence of these inter-peptide contacts underscores their significance in stabilising the dimers, potentially influencing the aggregation pathway. Moreover, the increased flexibility in the two termini, as well as the absence of persistent contacts in the metal-free dimer, correlates with the presence of amorphous aggregates. This phenomenon is known to mitigate fibrillation, while their absence in the metal-bound dimer suggests an increased propensity to form fibrils in the presence of copper ions.

Exploring the impact of mutation and post-translational modification on α-Synuclein: Insights from molecular dynamics simulations with and without copper.

We report molecular dynamics simulations of two modifications to α-Synuclein, namely A53T mutation and phosphorylation at Ser129, which have been observed in Parkinson's disease patients. Both modifications are close to known metal binding sites, so as well as each modified peptide we also study Cu(II) bound to N-terminal and C-terminal residues. We show that A53T is predicted to cause increased ß-sheet content of the peptide, with a persistent ß-hairpin between residues 35-55 particularly notable. Phosphorylation has less effect on secondary structure but is predicted to significantly increase the size of the peptide, especially when bound to Cu(II), which is ascribed to reduced interaction of C-terminal sequence with central non-amyloid component. In addition, estimate of binding free energy to Cu(II) indicates A53T has little effect on metal-ion affinity, whereas phosphorylation markedly enhances the strength of binding. We suggest that the predicted changes in spatial extent and secondary structure of α-Synuclein may have implications for aggregation into Lewy bodies.

Evaluation of implicit solvent models in molecular dynamics simulation of α-Synuclein.

We report conventional and accelerated molecular dynamics simulations of α-Synuclein, designed to assess performance of using different starting conformation, solvation environment and force field combination. Backbone and sidechain chemical shifts, radius of gyration, presence of ß-hairpin structures in KTK(E/Q)GV repeats and secondary structure percentages were used to evaluate how variations in forcefield, solvation model and simulation protocol provide results that correlate with experimental findings. We show that with suitable choice of forcefield and solvent, ff03ws and OBC implicit model, respectively, acceptable reproduction of experimental data on size and secondary structure is obtained by both conventional and accelerated MD. In contrast to the implicit solvent model, simulations in explicit TIP4P/2005 solvent do not properly represent size or secondary structure of α-Synuclein.Communicated by Ramaswamy H. Sarma.

How Cu(II) binding affects structure and dynamics of α-synuclein revealed by molecular dynamics simulations.

We report accelerated molecular dynamics simulations of α-Synuclein and its complex with two Cu(II) ions bound to experimentally determined binding sites. Adding two Cu(II) ions, one bound to the N-terminal region and one to the C-terminus, decreases size and flexibility of the peptide while introducing significant new contacts within and between N-terminus and non-Aß component (NAC). Cu(II) ions also alter the pattern of secondary structure within the peptide, inducing more and longer-lasting elements of secondary structure such as ß-strands and hairpins. Free energy surfaces, obtained from reweighting the accelerated molecular dynamics boost potential, further demonstrate the restriction on size and flexibility that results from binding of copper ions.

Forcefield evaluation and accelerated molecular dynamics simulation of Zn(II) binding to N-terminus of amyloid-ß.

We report conventional and accelerated molecular dynamics simulation of Zn(II) bound to the N-terminus of amyloid-ß. By comparison against NMR data for the experimentally determined binding mode, we find that certain combinations of forcefield and solvent model perform acceptably in describing the size, shape and secondary structure, and that there is no appreciable difference between implicit and explicit solvent models. We therefore used the combination of ff14SB forcefield and GBSA solvent model to compare the result of different binding modes of Zn(II) to the same peptide, using accelerated MD to enhance sampling and comparing the free peptide simulated in the same way. We show that Zn(II) imparts significant rigidity to the peptide, disrupts the secondary structure and pattern of salt bridges seen in the free peptide, and induces closer contact between residues. Free energy surfaces in 1 or 2 dimensions further highlight the effect of metal coordination on peptide's spatial extent. We also provide evidence that accelerated MD provides improved sampling over conventional MD by visiting as many or more configurations in much shorter simulation times.