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

Atomistic molecular simulations of Aβ-Zn conformational ensembles

The amyloid-forming Aβ peptide is able to interact with metal cations to form very stable complexes that influence fibril formation and contribute to the onset of Alzheimer's disease. Multiple structures of peptides derived from Aβ in complex with different metals have been resolved experimentally to provide an atomic-level description of the metal-protein interactions. However, Aβ is intrinsically disordered, and hence more amenable to an ensemble description. Molecular dynamics simulations can now reach the timescales needed to generate ensembles for these type of complexes. However, this requires accurate force fields both for the protein and the protein-metal interactions. Here we use state-of-the-art methods to generate force field parameters for the Zn(II) cations in a set of complexes formed by different Aβ variants and combine them with the Amber99SB*-ILDN optimized force field. Upon comparison of NMR experiments with the simulation results, further optimized with a Bayesian/Maximum entropy approach, we provide an accurate description of the molecular ensembles for most Aβ-metal complexes. We find that the resulting conformational ensembles are more heterogeneous than the NMR models deposited in the Protein Data Bank.

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.

Accelerated Molecular Dynamics to Explore the Binding of Transition Metals to Amyloid-β

We report the accelerated molecular dynamics (aMD) simulation of amyloid-β (Aβ) peptides of four different lengths (16, 28, 40, and 42 residues) and their complexes when bound to Cu(II), Fe(II), or Zn(II). 600 ns equilibrated trajectory data were analyzed for each structure from three independent 200 ns aMD simulations, generating 16 aMD trajectories. We show that the presence of a metal ion leads to reduced size and decreased mobility relative to the free peptide due to the anchoring effect of the ions. The reduced mobility was shown largely to be due to the restricted movement in N-terminal residues, most notably Asp1 and His6 that are involved in the metal-ion coordination in all cases. Significant disruption of the secondary structure and patterns of salt bridge interactions arise on the coordination of metal ions. In this regard, similarities were noted between results for Zn(II) and Fe(II), whereas results for Cu(II) are more comparable to that of the free peptides. Reweighting of free energy surfaces was carried out from aMD data to identify the properties and descriptions of local minima structures.

Zinc Binding to NAP-Type Neuroprotective Peptides: Nuclear Magnetic Resonance Studies and Molecular Modeling

Aggregation of amyloid-β peptides (Aβ) is a hallmark of Alzheimer’s disease (AD), which is affecting an increasing number of people. Hence, there is an urgent need to develop new pharmaceutical treatments which could be used to prevent the AD symptomatology. Activity-dependent neuroprotective protein (ADNP) was found to be deficient in AD, whereas NAP, an 8-amino-acid peptide (¹NAPVSIPQ⁸) derived from ADNP, was shown to enhance cognitive function. The higher tendency of zinc ion to induce Aβ aggregation and formation of amorphous aggregates is also well-known in the scientific literature. Although zinc binding to Aβ peptides was extensively investigated, there is a shortage of knowledge regarding the relationship between NAP peptide and zinc ions. Therefore, here, we investigated the binding of zinc ions to the native NAP peptide and its analog obtained by replacing the serine residue in the NAP sequence with tyrosine (¹NAPVYIPQ⁸) at various molar ratios and pH values by mass spectrometry (MS) and nuclear magnetic resonancespectroscopy (NMR). Matrix-assisted laser desorption/ionization time-of-flight (MALDI ToF) mass spectrometry confirmed the binding of zinc ions to NAP peptides, while the chemical shift of Asp¹, observed in ¹H-NMR spectra, provided direct evidence for the coordinating role of zinc in the N-terminal region. In addition, molecular modeling has also contributed largely to our understanding of Zn binding to NAP peptides.