The Alzheimer's peptides Abeta40 and 42 adopt distinct conformations in water: a combined MD / NMR study

Molecular Simulations Reveal Terminal Group Mediated Stabilization of Helical Conformers in Both Amyloid-β42 and α-Synuclein

The presence of partially structured helices in natively unfolded amyloid-β42 (Aβ42) and α-synuclein (αS) has been shown to accelerate fibrillation in the onset of Alzheimer's and Parkinson's disease, respectively. At the other extreme, folded stable helical conformers have also been reported to resist amyloid formation. Recent studies indicate that amyloidogenic aggregation can be impeded using small molecules that stabilize the α-helical monomers and switch off the neurotoxic pathway. We predict a common intrapeptide route to stabilization based on the plasticity of helical conformations of Aβ42 and αS as assessed through extensive atomistic molecular dynamics (MD) computer simulations (∼36 μs) across ten distinct protein force field and water model combinations. Computed free energies and interaction maps (not obtainable from experiments alone) show that flexible terminal groups (N-terminus of Aβ42 and C-terminus of αS) show a tendency to stabilize folded helical conformations in both peptides via primary hydrophobic interactions with central hydrophobic domains, and secondary salt bridges with other domains. These interactions confer aggregation resistance by decreasing the population of partially structured helices and are absent in control simulations of complete unfolding. Computed helical stability is also significantly reduced in terminal-deleted variants. The models suggest new strategies to tackle neurodegeneration by rationally re-engineering terminal groups to optimize their predicted ability to deactivate helical monomers.

Femtosecond Hydration Map of Intrinsically Disordered α-Synuclein

Protein hydration water plays a fundamentally important role in protein folding, binding, assembly, and function. Little is known about the hydration water in intrinsically disordered proteins that challenge the conventional sequence-structure-function paradigm. Here, by combining experiments and simulations, we show the existence of dynamical heterogeneity of hydration water in an intrinsically disordered presynaptic protein, namely α-synuclein, implicated in Parkinson's disease. We took advantage of nonoccurrence of cysteine in the sequence and incorporated a number of cysteine residues at the N-terminal segment, the central amyloidogenic nonamyloid-β component (NAC) domain, and the C-terminal end of α-synuclein. We then labeled these cysteine variants using environment-sensitive thiol-active fluorophore and monitored the solvation dynamics using femtosecond time-resolved fluorescence. The site-specific femtosecond time-resolved experiments allowed us to construct the hydration map of α-synuclein. Our results show the presence of three dynamically distinct types of water: bulk, hydration, and confined water. The amyloidogenic NAC domain contains dynamically restrained water molecules that are strikingly different from the water molecules present in the other two domains. Atomistic molecular dynamics simulations revealed longer residence times for water molecules near the NAC domain and supported our experimental observations. Additionally, our simulations allowed us to decipher the molecular origin of the dynamical heterogeneity of water in α-synuclein. These simulations captured the quasi-bound water molecules within the NAC domain originating from a complex interplay between the local chain compaction and the sequence composition. Our findings from this synergistic experimental simulation approach suggest longer trapping of interfacial water molecules near the amyloidogenic hotspot that triggers the pathological conversion into amyloids via chain sequestration, chain desolvation, and entropic liberation of ordered water molecules.

The role of the protein-RNA recognition code in neurodegeneration

MicroRNAs are small endogenous RNAs that pair and bind to sites on mRNAs to direct post-transcriptional repression. However, there is a possibility that microRNAs directly influence protein structure and activity, and this influence can be termed post-translational riboregulation. This conceptual review explores the literature on neurodegenerative disorders. Research on the association between neurodegeneration and RNA-repeat toxicity provides data that support a protein-RNA recognition code. For example, this code explains why hnRNP H and SFPQ proteins, which are involved in amyotrophic lateral sclerosis, are sequestered by the (GGGGCC)n repeat sequence. Similarly, it explains why MNBL proteins and (CTG)n repeats in RNA, which are involved in myotonic dystrophy, are sequestered into RNA foci. Using this code, proteins involved in diseases can be identified. A simple protein BLAST search of the human genome for amino acid repeats that correspond to the nucleotide repeats reveals new proteins among already known proteins that are involved in diseases. For example, the (CAG)n repeat sequence, when transcribed into possible peptide sequences, leads to the identification of PTCD3, Rem2, MESP2, SYPL2, WDR33, COL23A1, and others. After confirming this approach on RNA repeats, in the next step, the code was used in the opposite manner. Proteins that are involved in diseases were compared with microRNAs involved in those diseases. For example, a reasonable correspondence of microRNA 9 and 107 with amyloid-β-peptide (Aβ42) was identified. In the last step, a miRBase search for micro-nucleotides, obtained by transcription of a prion amino acid sequence, revealed new microRNAs and microRNAs that have previously been identified as involved in prion diseases. This concept provides a useful key for designing RNA or peptide probes.

Trafficking and proteolytic processing of amyloid precursor protein and secretases in Alzheimer's disease development: An up-to-date review

Alzheimer's disease (AD), which is predicted to affect 1 in 85 persons worldwide by 2050, results in progressive loss of neuronal functions, leading to impairments in memory and cognitive abilities. As being one of the major neuropathological hallmarks of AD, senile plaques mainly consist of amyloid-β (Aβ) peptides, which are derived from amyloid precursor protein (APP) via the sequential cleavage by β- and γ-secretases. Although the overproduction and accumulation of Aβ peptides are at the center of AD research, the new discoveries point out to the complexity of the disease development. In this respect, it is crucial to understand the processing and the trafficking of APP, the enzymes involved in its processing, the cleavage products and their therapeutic potentials. This review summarizes the salient features of APP processing focusing on APP, the canonical secretases as well as the novel secretases and the cleavage products with an update of the recent developments. We also discussed the intracellular trafficking of APP and secretases in addition to their potential in AD therapy.

O-GlcNAc Modification Protects against Protein Misfolding and Aggregation in Neurodegenerative Disease

Post-translational modifications (PTMs) of proteins are becoming the focus of intense research due to their implications in a broad spectrum of neurodegenerative diseases. Various PTMs have been identified to alter the toxic profiles of proteins which play critical roles in disease etiology. In Alzheimer's disease (AD), dysregulated phosphorylation is reported to promote pathogenic processing of the microtubule-associated tau protein. Among the PTMs, the enzymatic addition of N-acetyl-d-glucosamine (GlcNAc) residues to Ser/Thr residues is reported to deliver protective effects against the pathogenic processing of both amyloid precursor protein (APP) and tau. Modification of tau with as few as one single O-GlcNAc residue inhibits its toxic self-assembly. This modification also has the same effect on the assembly of the Parkinson's disease (PD) associated α-synuclein (ASyn) protein. In fact, O-GlcNAcylation ( O-linked GlcNAc modification) affects the processing of numerous proteins implicated in AD, PD, amyotrophic lateral sclerosis (ALS), and Huntington's disease (HD) in a similar manner. As such, manipulation of a protein's O-GlcNAcylation status has been proposed to offer therapeutic routes toward addressing multiple neurodegenerative pathologies. Here we review the various effects that O-GlcNAc modification, and its modulated expression, have on pathogenically significant proteins involved in neurodegenerative disease.

The combined force field-sampling problem in simulations of disordered amyloid-β peptides

Molecular dynamics simulations of intrinsically disordered proteins (IDPs) can provide high resolution structural ensembles if the force field is accurate enough and if the simulation sufficiently samples the conformational space of the IDP with the correct weighting of sub-populations. Here, we investigate the combined force field-sampling problem by testing a standard force field as well as newer fixed charge force fields, the latter specifically motivated for better description of unfolded states and IDPs, and comparing them with a standard temperature replica exchange (TREx) protocol and a non-equilibrium Temperature Cool Walking (TCW) sampling algorithm. The force field and sampling combinations are used to characterize the structural ensembles of the amyloid-beta peptides Aβ42 and Aβ43, which both should be random coils as shown recently by experimental nuclear magnetic resonance (NMR) and 2D Förster resonance energy transfer (FRET) experiments. The results illustrate the key importance of the sampling algorithm: while the standard force field using TREx is in poor agreement with the NMR J-coupling and nuclear Overhauser effect and 2D FRET data, when using the TCW method, the standard and optimized protein-water force field combinations are in very good agreement with the same experimental data since the TCW sampling method produces qualitatively different ensembles than TREx. We also discuss the relative merit of the 2D FRET data when validating structural ensembles using the different force fields and sampling protocols investigated in this work for small IDPs such as the Aβ42 and Aβ43 peptides.

Separation of β-Amyloid Tryptic Peptide Species with Isomerized and Racemized l-Aspartic Residues with Ion Mobility in Structures for Lossless Ion Manipulations

Accumulation of β-amyloid (Aβ) is one of the hallmarks of Alzheimer's disease. The deposition of β-amyloid plaques is likely to start years in advance of manifestation of clinical symptoms, although the exact timing is unknown. Over the years, Aβ peptides undergo both post-translational modification and stereoisomerization. Analysis of the resulting stereoisomers is particularly challenging because of their identical elemental composition and similar physicochemical properties. Herein, we have utilized our recently developed structures for lossless ion manipulations ion mobility-mass spectrometry platform (SLIM IM-MS), in conjunction with serpentine ultralong path with extended routing (SUPER), to baseline resolve four distinct sets of Aβ17-28 tryptic peptide epimers on a rapid (∼1 s) time scale. We discovered that sodium adduct ions, [M + H + Na]2+, allowed baseline SLIM SUPER IM resolution for all Aβ epimer sets assessed, while such baseline separations were unachievable for their [M + 2H]2+ doubly protonated ions.

General Purpose Water Model Can Improve Atomistic Simulations of Intrinsically Disordered Proteins

Unconstrained atomistic simulations of intrinsically disordered proteins and peptides (IDP) remain a challenge: widely used, "general purpose" water models tend to favor overly compact structures relative to experiment. Here we have performed a total of 93 μs of unrestrained MD simulations to explore, in the context of IDPs, a recently developed "general-purpose" 4-point rigid water model OPC, which describes liquid state of water close to experiment. We demonstrate that OPC, together with a popular AMBER force field ff99SB, offers a noticeable improvement over TIP3P in producing more realistic structural ensembles of three common IDPs benchmarks: 55-residue apo N-terminal zinc-binding domain of HIV-1 integrase ("protein IN"), amyloid β-peptide (Aβ₄₂) (residues 1-42), and 26-reside H4 histone tail. As a negative control, computed folding profile of a regular globular miniprotein (CLN025) in OPC water is in appreciably better agreement with experiment than that obtained in TIP3P, which tends to overstabilize the compact native state relative to the extended conformations. We employed Aβ₄₂ peptide to investigate the possible influence of the solvent box size on simulation outcomes. We advocate a cautious approach for simulations of IDPs: we suggest that the solvent box size should be at least four times the radius of gyration of the random coil corresponding to the IDP. The computed free energy landscape of protein IN in OPC resembles a shallow "tub" - conformations with substantially different degrees of compactness that are within 2 k_B T of each other. Conformations with very different secondary structure content coexist within 1 k_B T of the global free energy minimum. States with higher free energy tend to have less secondary structure. Computed low helical content of the protein has virtually no correlation with its degree of compactness, which calls into question the possibility of using the helicity as a metric for assessing performance of water models for IDPs, when the helicity is low. Predicted radius of gyration ( R _g) of H4 histone tail in OPC water falls in-between that of a typical globular protein and a fully denatured protein of the same size; the predicted R _g is consistent with two independent predictions. In contrast, H4 tail in TIP3P water is as compact as the corresponding globular protein. The computed free energy landscape of H4 tail in OPC is relatively flat over a significant range of compactness, which, we argue, is consistent with its biological function as facilitator of internucleosome interactions.

Molecular Dynamics Simulations Combined with Nuclear Magnetic Resonance and/or Small-Angle X-ray Scattering Data for Characterizing Intrinsically Disordered Protein Conformational Ensembles

The concept of intrinsically disordered proteins (IDPs) has emerged relatively slowly, but over the past 20 years, it has become an intense research area in structural biology. Indeed, because of their considerable flexibility and structural heterogeneity, the determination of IDP conformational ensemble is particularly challenging and often requires a combination of experimental measurements and computational approaches. With the improved accuracy of all-atom force fields and the increasing computing performances, molecular dynamics (MD) simulations have become more and more reliable to generate realistic conformational ensembles. And the combination of MD simulations with experimental approaches, such as nuclear magnetic resonance (NMR) and/or small-angle X-ray scattering (SAXS) allows one to converge toward a more accurate and exhaustive description of IDP structures. In this Review, we discuss the state of the art of MD simulations of IDP conformational ensembles, with a special focus on studies that back-calculated and directly compared theoretical and experimental NMR or SAXS observables, such as chemical shifts (CS), ³J-couplings (³Jc), residual dipolar couplings (RDC), or SAXS intensities. We organize the review in three parts. In the first section, we discuss the studies which used NMR and/or SAXS data to test and validate the development of force fields or enhanced sampling techniques. In the second part, we explore different methods for the refinement of MD-derived structural ensembles, such as NMR or SAXS data-restrained MD simulations or ensemble reweighting to better fit experiments. Finally, we survey some recent studies combining MD simulations with NMR and/or SAXS measurements to investigate the relationship between IDP conformational ensemble and biological activity, as well as their implication in human diseases. From this review, we noticed that quite a few studies compared MD-generated conformational ensembles with both NMR and SAXS measurements to validate IDP structures at both local and global levels. Yet, beside the IDP propensity to form local secondary structures, their dynamic extension or compactness also appears important for their activity. Thus, we believe that a close synergy between MD simulations, NMR, and SAXS experiments would be greatly appropriate to address the challenges of characterizing the disordered structures of proteins and their complexes, relative to their biological functions.

Solvent Composition Effects on the Structural Properties of the Aβ42 Monomer from the 3D-RISM-KH Molecular Theory of Solvation

Structural characterization of amyloid (A)β peptides implicated in Alzheimer's disease is a challenging problem due to their intrinsically disordered nature and their high propensity for aggregation. Only limited information is currently available from experiments on conformational properties and aggregation pathways of the peptides in cellular environments. In silico modeling complements experimental information, providing atomistic insight into structure and dynamics of different Aβ species. All-atom explicit solvent molecular dynamics (MD) simulations with a properly selected force field can deliver reliable structural and dynamic information. In the case of intrinsically disordered Aβ peptides, enhanced sampling simulations beyond the nanosecond time scale are required to obtain statistically meaningful results even for simple solvent conditions. To overcome the challenges of conformational sampling in crowded cellular environments, alternative approaches have to be used, including postprocessing of MD data. In this study, we employ the statistical-mechanical, three-dimensional reference interaction site model with the Kovalenko-Hirata closure integral equation molecular theory of solvation to describe solvent composition effects on the conformational equilibrium in a structural ensemble of the Aβ42 (covering residues 1-42) monomer based on a statistical reweighting technique. The methodology enables a computationally efficient prediction on how different factors in the cellular environment, such as solvent composition, nonpolar solvation, and macromolecular crowding, affect the structural properties of the monomer. Similarities have been identified between changes in the structural ensemble caused by nonpolar solvation and crowded environments modeled by ionic solution with large negative ions. In particular, both solvent conditions reduce the random coil content and enhance the helical structure content of the monomer. In contrast to the previous studies, which reported increased α-helical content of peptides in crowded environments, this work attributes these structural features to the difference in solvent exposure of hydrophilic residues of the monomer for different secondary structure elements, rather than to (entropic) excluded volume effects.

Highly Disordered Amyloid-β Monomer Probed by Single-Molecule FRET and MD Simulation

Monomers of amyloid-β (Aβ) protein are known to be disordered, but there is considerable controversy over the existence of residual or transient conformations that can potentially promote oligomerization and fibril formation. We employed single-molecule Förster resonance energy transfer (FRET) spectroscopy with site-specific dye labeling using an unnatural amino acid and molecular dynamics simulations to investigate conformations and dynamics of Aβ isoforms with 40 (Aβ40) and 42 residues (Aβ42). The FRET efficiency distributions of both proteins measured in phosphate-buffered saline at room temperature show a single peak with very similar FRET efficiencies, indicating there is apparently only one state. 2D FRET efficiency-donor lifetime analysis reveals, however, that there is a broad distribution of rapidly interconverting conformations. Using nanosecond fluorescence correlation spectroscopy, we measured the timescale of the fluctuations between these conformations to be ∼35 ns, similar to that of disordered proteins. These results suggest that both Aβ40 and Aβ42 populate an ensemble of rapidly reconfiguring unfolded states, with no long-lived conformational state distinguishable from that of the disordered ensemble. To gain molecular-level insights into these observations, we performed molecular dynamics simulations with a force field optimized to describe disordered proteins. We find, as in experiments, that both peptides populate configurations consistent with random polymer chains, with the vast majority of conformations lacking significant secondary structure, giving rise to very similar ensemble-averaged FRET efficiencies.

Binding Modes of Phthalocyanines to Amyloid β Peptide and Their Effects on Amyloid Fibril Formation

The inherent tendency of proteins to convert from their native states into amyloid aggregates is associated with a range of human disorders, including Alzheimer's and Parkinson's diseases. In that sense, the use of small molecules as probes for the structural and toxic mechanism related to amyloid aggregation has become an active area of research. Compared with other compounds, the structural and molecular basis behind the inhibitory interaction of phthalocyanine tetrasulfonate (PcTS) with proteins such as αS and tau has been well established, contributing to a better understanding of the amyloid aggregation process in these proteins. We present here the structural characterization of the binding of PcTS and its Cu(II) and Zn(II)-loaded forms to the amyloid β-peptide (Aβ) and the impact of these interactions on the peptide amyloid fibril assembly. Elucidation of the PcTS binding modes to Aβ₄₀ revealed the involvement of specific aromatic and hydrophobic interactions in the formation of the Aβ₄₀-PcTS complex, ascribed to a binding mode in which the planarity and hydrophobicity of the aromatic ring system in the phthalocyanine act as main structural determinants for the interaction. Our results demonstrated that formation of the Aβ₄₀-PcTS complex does not interfere with the progression of the peptide toward the formation of amyloid fibrils. On the other hand, conjugation of Zn(II) but not Cu(II) at the center of the PcTS macrocyclic ring modified substantially the binding profile of this phthalocyanine to Aβ₄₀ and became crucial to reverse the effects of metal-free PcTS on the fibril assembly of the peptide. Overall, our results provide a firm basis to understand the structural rules directing phthalocyanine-protein interactions and their implications on the amyloid fibril assembly of the target proteins; in particular, our results contradict the hypothesis that PcTS might have similar mechanisms of action in slowing the formation of a variety of pathological aggregates.

Solid-state NMR reveals a comprehensive view of the dynamics of the flexible, disordered N-terminal domain of amyloid-β fibrils

Amyloid fibril deposits observed in Alzheimer's disease comprise amyloid-β (Aβ) protein possessing a structured hydrophobic core and a disordered N-terminal domain (residues 1-16). The internal flexibility of the disordered domain is likely essential for Aβ aggregation. Here, we used ²H static solid-state NMR methods to probe the dynamics of selected side chains of the N-terminal domain of Aβ_1-40 fibrils. Line shape and relaxation data suggested a two-state model in which the domain's free state undergoes a diffusive motion that is quenched in the bound state, likely because of transient interactions with the structured C-terminal domain. At 37 °C, we observed freezing of the dynamics progressively along the Aβ sequence, with the fraction of the bound state increasing and the rate of diffusion decreasing. We also found that without solvation, the diffusive motion is quenched. The solvent acted as a plasticizer reminiscent of its role in the onset of global dynamics in globular proteins. As the temperature was lowered, the fraction of the bound state exhibited sigmoidal behavior. The midpoint of the freezing curve coincided with the bulk solvent freezing for the N-terminal residues and increased further along the sequence. Using ²H R _1ρ measurements, we determined the conformational exchange rate constant between the free and bound states under physiological conditions. Zinc-induced aggregation leads to the enhancement of the dynamics, manifested by the faster conformational exchange, faster diffusion, and lower freezing-curve midpoints.

Characterizing the structural and thermodynamic properties of Aβ42 and Aβ40

The self-assembly of amyloid-beta (Aβ) proteins in aqueous extracellular environments is implicated in Alzheimer's disease. Among several alloforms of Aβ proteins differing in sequence length, the 42- and 40-residue forms (Aβ42 and Aβ40) are the most abundant ones in the human body. Although the only difference is the additional I₄₁A₄₂ residues in the C-terminus, Aβ42 exhibits more aggregation tendency and stronger neurotoxicity than Aβ40. Here, we investigate the molecular factors that confer more aggregation potential to Aβ42 than to Aβ40 based on molecular dynamics simulations combined with solvation thermodynamic analyses. It is observed that the most salient structural feature of Aβ42 relative to Aβ40 is the more enhanced β-sheet forming tendency, in particular in the C-terminal region. While such a structural characteristic of Aβ42 will certainly serve to facilitate the formation of aggregate species rich in β-sheet structure, we also detect its interesting thermodynamic consequence. Indeed, we find from the decomposition analysis that the C-terminal region substantially increases the solvation free energy (i.e., overall "hydrophobicity") of Aβ42, which is caused by the dehydration of the backbone moieties showing the enhanced tendency of forming the β-structure. Together with the two additional hydrophobic residues (I₄₁A₄₂), this leads to the higher solvation free energy of Aβ42, implying the larger water-mediated attraction toward the self-assembly. Thus, our computational results provide structural and thermodynamic grounds on why Aβ42 has more aggregation propensity than Aβ40 in aqueous environments.

Pathways of amyloid-beta absorption and aggregation in a membranous environment

Aggregation of misfolded oligomeric amyloid-beta (Aβ) peptides on lipid membranes has been identified as a primary event in Alzheimer's pathogenesis. However, the structural and dynamical features of this membrane assisted Aβ aggregation have not been well characterized. The microscopic characterization of dynamic molecular-level interactions in peptide aggregation pathways has been challenging both computationally and experimentally. In this work, we explore differential patterns of membrane-induced Aβ 16-22 (K-L-V-F-F-A-E) aggregation from the microscopic perspective of molecular interactions. Physics-based coarse-grained molecular dynamics (CG-MD) simulations were employed to investigate the effect of lipid headgroup charge - zwitterionic (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine: POPC) and anionic (1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-l-serine: POPS) - on Aβ 16-22 peptide aggregation. Our analyses present an extensive overview of multiple pathways for peptide absorption and biomechanical forces governing peptide folding and aggregation. In agreement with experimental observations, anionic POPS molecules promote extended configurations in Aβ peptides that contribute towards faster emergence of ordered β-sheet-rich peptide assemblies compared to POPC, suggesting faster fibrillation. In addition, lower cumulative rates of peptide aggregation in POPS due to higher peptide-lipid interactions and slower lipid diffusion result in multiple distinct ordered peptide aggregates that can serve as nucleation seeds for subsequent Aβ aggregation. This study provides an in-silico assessment of experimentally observed aggregation patterns, presents new morphological insights and highlights the importance of lipid headgroup chemistry in modulating the peptide absorption and aggregation process.

Significantly different contact patterns between Aβ40 and Aβ42 monomers involving the N-terminal region

Aβ42 peptide, with two additional residues at C-terminus, aggregates much faster than Aβ40. We performed equilibrium replica-exchange molecular dynamics simulations of their monomers using our residue-specific force field. Simulated ³ J_{HNH α} -coupling constants agree excellently with experimental data. Aβ40 and Aβ42 have very similar local conformational features, with considerable β-strand structures in the segments: A2-H6 (A), L17-A21 (B), A30-V36 (C) of both peptides and V39-I41 (D) of Aβ42. Both peptides have abundant A-B and B-C contacts, but Aβ40 has much more contacts between A and C than Aβ42, which may retard its aggregation. Only Aβ42 has considerable A-B-C-D topology. Decreased probability of A-C contact in Aβ42 relates to the competition from C-D contact. Increased A-C contact probability may also explain the slower aggregation of A2T and A2V mutants of Aβ42.

On the Applicability of Force Fields To Study the Aggregation of Amyloidogenic Peptides Using Molecular Dynamics Simulations

Molecular dynamics simulations play an essential role in understanding biomolecular processes such as protein aggregation at temporal and spatial resolutions which are not attainable by experimental methods. For a correct modeling of protein aggregation, force fields must accurately represent molecular interactions. Here, we study the effect of five different force fields on the oligomer formation of Alzheimer's Aβ_16-22 peptide and two of its mutants: Aβ_16-22(F19V,F20V), which does not form fibrils, and Aβ_16-22(F19L) which forms fibrils faster than the wild type. We observe that while oligomer formation kinetics depends strongly on the force field, structural properties, such as the most relevant protein-protein contacts, are similar between them. The oligomer formation kinetics obtained with different force fields differ more from each other than the kinetics between aggregating and nonaggregating peptides simulated with a single force field. We discuss the difficulties in comparing atomistic simulations of amyloid oligomer formation with experimental observables.

Nonproductive Binding Modes as a Prominent Feature of Aβ40 Fiber Elongation: Insights from Molecular Dynamics Simulation

The formation of amyloid fibers has been implicated in a number of neurodegenerative diseases. The growth of amyloid fibers is strongly thermodynamically favorable, but kinetic traps exist where the incoming monomer binds in an incompatible conformation that blocks further elongation. Unfortunately, this process is difficult to follow experimentally at the atomic level. It is also too complex to simulate in full detail and to date has been explored either through coarse-grained simulations, which may miss many important interactions, or full atomic simulations, in which the incoming peptide is constrained to be near the ideal fiber geometry. Here we use an alternate approach starting from a docked complex in which the monomer is from an experimental NMR structure of one of the major conformations in the unbound ensemble, a largely unstructured peptide with the central hydrophobic region in a 3₁₀ helix. A 1000 ns full atomic simulation in explicit solvent shows the formation of a metastable intermediate by sequential, concerted movements of both the fiber and the monomer. A Markov state model shows that the unfolded monomer is trapped at the end of the fiber in a set of interconverting antiparallel β-hairpin conformations. The simulation here may serve as a model for the binding of other non-β-sheet conformations to amyloid fibers.

The effect of amyloid-β peptide on synaptic plasticity and memory is influenced by different isoforms, concentrations, and aggregation status

The increase of oligomeric amyloid-beta (oAβ) has been related to synaptic dysfunction, thought to be the earliest event in Alzheimer's disease pathophysiology. Conversely, the suppression of endogenous Aβ impaired synaptic plasticity and memory, suggesting that the peptide is needed in the healthy brain. However, different species, aggregation forms and concentrations of Aβ might differently influence synaptic function/dysfunction. Here, we have tested the contribution of monomeric and oligomeric Aβ42 and Aβ40 at 200 nM and 200 pM concentrations on hippocampal long-term potentiation and spatial memory. We found that, when at 200 nM, oAβ40, oAβ42, and monomeric Aβ42 impaired long-term potentiation and memory, whereas only oAβ42 200 pM enhanced synaptic plasticity and memory and rescued the detrimental effect due to depletion of endogenous Aβ. Interestingly, quantification of monomer-like and oligomer-like species carried out by transmission electron microscopy revealed an increase of the monomer/oligomer ratio in the oAβ42 200 pM preparation, suggesting that the content of monomers and oligomers depends on the final concentration of the solution.

Peptides, Peptidomimetics, and Carbohydrate-Peptide Conjugates as Amyloidogenic Aggregation Inhibitors for Alzheimer's Disease

Alzheimer's disease (AD) is a progressive neurodegenerative disorder accounting for 60-80% of dementia cases. For many years, AD causality was attributed to amyloid-β (Aβ) aggregated species. Recently, multiple therapies that target Aβ aggregation have failed in clinical trials, since Aβ aggregation is found in AD and healthy patients. Attention has therefore shifted toward the aggregation of the tau protein as a major driver of AD. Numerous inhibitors of tau-based pathology have recently been developed. Diagnosis of AD has shifted from measuring late stage senile plaques to early stage biomarkers, amyloid-β and tau monomers and oligomeric assemblies. Synthetic peptides and some derivative structures are being explored for use as theranostic tools as they possess the capacity both to bind the biomarkers and to inhibit their pathological self-assembly. Several studies have demonstrated that O-linked glycoside addition can significantly alter amyloid aggregation kinetics. Furthermore, natural O-glycosylation of amyloid-forming proteins, including amyloid precursor protein (APP), tau, and α-synuclein, promotes alternative nonamyloidogenic processing pathways. As such, glycopeptides and related peptidomimetics are being investigated within the AD field. Here we review advancements made in the last 5 years, as well as the arrival of sugar-based derivatives.

High pressure NMR reveals conformational perturbations by disease-causing mutations in amyloid β-peptide

Here we present the high pressure NMR characterization of Aβ42 and two Aβ40 variants with Alzheimer-causing mutations E22G and D23N. While chemical shifts only identified localized changes at ambient pressure compared with Aβ40, high pressure NMR revealed a common site with heightened pressure sensitivity at Q15, K16 and L17 in all three variants, which correlates to higher β-propensity at central hydrophobic cluster (CHC) and faster aggregation.

Weaker N-Terminal Interactions for the Protective over the Causative Aβ Peptide Dimer Mutants

Knowing that abeta amyloid peptide (Aβ₄₂) dimers are the smallest and most abundant neurotoxic oligomers for Alzheimer's disease (AD), we used molecular simulations with advanced sampling methods (replica-exchange) to characterize and compare interactions between the N-termini (residues 1-16) of wild type (WT-WT) and five mutant dimers under constrained and unconstrained conditions. The number of contacts and distances between the N-termini, and contact maps of their conformational landscape illustrate substantial differences for a single residue change. The N-terminal contacts are significantly diminished for the dimers containing the monomers that protect against (WT-A2T) as compared with those that predispose toward (A2V-A2V) AD and for the control WT-WT dimers. The reduced number of N-terminal contacts not only occurs at or near the second residue mutations but also is distributed through to the 10th residue. These findings provide added support to the accumulating evidence for the "N-terminal hypothesis of AD" and offer an alternate mechanism for the cause of protection from the A2T mutant.

Developing a molecular dynamics force field for both folded and disordered protein states

Many proteins that perform important biological functions are completely or partially disordered under physiological conditions. Molecular dynamics simulations could be a powerful tool for the structural characterization of such proteins, but it has been unclear whether the physical models (force fields) used in simulations are sufficiently accurate. Here, we systematically compare the accuracy of a number of different force fields in simulations of both ordered and disordered proteins, finding that each force field has strengths and limitations. We then describe a force field that substantially improves on the state-of-the-art accuracy for simulations of disordered proteins without sacrificing accuracy for folded proteins, thus broadening the range of biological systems amenable to molecular dynamics simulations.

Molecular dynamics (MD) simulation is a valuable tool for characterizing the structural dynamics of folded proteins and should be similarly applicable to disordered proteins and proteins with both folded and disordered regions. It has been unclear, however, whether any physical model (force field) used in MD simulations accurately describes both folded and disordered proteins. Here, we select a benchmark set of 21 systems, including folded and disordered proteins, simulate these systems with six state-of-the-art force fields, and compare the results to over 9,000 available experimental data points. We find that none of the tested force fields simultaneously provided accurate descriptions of folded proteins, of the dimensions of disordered proteins, and of the secondary structure propensities of disordered proteins. Guided by simulation results on a subset of our benchmark, however, we modified parameters of one force field, achieving excellent agreement with experiment for disordered proteins, while maintaining state-of-the-art accuracy for folded proteins. The resulting force field, a99SB-disp, should thus greatly expand the range of biological systems amenable to MD simulation. A similar approach could be taken to improve other force fields.

Impact of membrane curvature on amyloid aggregation

The misfolding, amyloid aggregation, and fibril formation of intrinsically disordered proteins/peptides (or amyloid proteins) have been shown to cause a number of disorders. The underlying mechanisms of amyloid fibrillation and structural properties of amyloidogenic precursors, intermediates, and amyloid fibrils have been elucidated in detail; however, in-depth examinations on physiologically relevant contributing factors that induce amyloidogenesis and lead to cell death remain challenging. A large number of studies have attempted to characterize the roles of biomembranes on protein aggregation and membrane-mediated cell death by designing various membrane components, such as gangliosides, cholesterol, and other lipid compositions, and by using various membrane mimetics, including liposomes, bicelles, and different types of lipid-nanodiscs. We herein review the dynamic effects of membrane curvature on amyloid generation and the inhibition of amyloidogenic proteins and peptides, and also discuss how amyloid formation affects membrane curvature and integrity, which are key for understanding relationships with cell death. Small unilamellar vesicles with high curvature and large unilamellar vesicles with low curvature have been demonstrated to exhibit different capabilities to induce the nucleation, amyloid formation, and inhibition of amyloid-β peptides and α-synuclein. Polymorphic amyloidogenesis in small unilamellar vesicles was revealed and may be viewed as one of the generic properties of interprotein interaction-dominated amyloid formation. Several mechanical models and phase diagrams are comprehensively shown to better explain experimental findings. The negative membrane curvature-mediated mechanisms responsible for the toxicity of pancreatic β cells by the amyloid aggregation of human islet amyloid polypeptide (IAPP) and binding of the precursors of the semen-derived enhancer of viral infection (SEVI) are also described. The curvature-dependent binding modes of several types of islet amyloid polypeptides with high-resolution NMR structures are also discussed.

Time- and dose-related effects of amyloid beta1-40 on retina and optic nerve morphology in rats

PURPOSE

Amyloid beta (Aβ) is known to contribute to the pathophysiology of retinal neurodegenerative diseases such as glaucoma. Effects of intravitreal Aβ(1-42) on retinal and optic nerve morphology in animal models have widely been studied but not those of Aβ(1-40). Hence, we evaluated the time- and dose-related effects of intravitreal Aβ(1-40) on retinal and optic nerve morphology. Since oxidative stress and brain derived neurotrophic factor (BDNF) are associated with Aβ-induced neuronal damage, we also studied dose and time-related effects of Aβ(1-40) on retinal oxidative stress and BDNF levels.

MATERIALS AND METHODS

Five groups of rats were intravitreally administered with vehicle or Aβ(1-40) in doses of 1.0, 2.5, 5 and 10 nmol. Animals were sacrificed and eyes were enucleated at weeks 1, 2 and 4 post-injection. The retinae were subjected to morphometric analysis and TUNEL staining. Optic nerve sections were stained with toluidine blue and were graded for neurodegenerative effects. The estimation of BDNF and markers of oxidative stress in retina were done using ELISA technique.

RESULTS AND CONCLUSIONS

It was observed that intravitreal Aβ(1-40) causes significant retinal and optic nerve damage up to day 14 post-injection and there was increasing damage with increase in dose. However, on day 30 post-injection both the retinal and optic nerve morphology showed a trend towards normalization. The observations made for retinal cell apoptosis, retinal glutathione, superoxide dismutase activity and BDNF were in accordance with those of morphological changes with deterioration till day 14 and recovery by day 30 post-injection. The findings of this study may provide a guide for selection of appropriate experimental conditions for future studies.

Glycation induces conformational changes in the amyloid-β peptide and enhances its aggregation propensity: molecular insights

The cytotoxicity of the amyloid beta (Aβ) peptide, implicated in the pathogenesis of Alzheimer's disease (AD), can be enhanced by its post-translational glycation, a series of non-enzymatic reactions with reducing sugars and reactive dicarbonyls. However, little is known about the underlying mechanisms that potentially enhance the cytotoxicity of the advanced glycation modified Aβ. In this work, fully atomistic molecular dynamics (MD) simulations are exploited to obtain direct molecular insights into the process of early Aβ self-assembly in the presence and absence of glycated lysine residues. Analyses of data exceeding cumulative timescales of 1 microsecond for each system reveal that glycation results in a stronger enthalpy of association between Aβ monomers and lower conformational entropy, in addition to a sharp overall increase in the beta-sheet content. Further analyses reveal that the enhanced interactions originate, in large part, due to markedly stronger, as well as new, inter-monomer salt bridging propensities in the glycated variety. Interestingly, these conformational and energetic effects are broadly reflected in preformed protofibrillar forms of Aβ small oligomers modified with glycation. Our combined results imply that glycation consolidates Aβ self-assembly regardless of its point of occurrence in the pathway. They provide a basis for further mechanistic studies and therapeutic endeavors that could potentially result in novel ways of combating AGE related AD progression.

Critical aggregation concentration for the formation of early Amyloid-β (1-42) oligomers

The oligomers formed during the early steps of amyloid aggregation are thought to be responsible for the neurotoxic damage associated with Alzheimer’s disease. It is therefore of great interest to characterize this early aggregation process and the aggregates formed, especially for the most significant peptide in amyloid fibrils, Amyloid-β(1–42) (Aβ42). For this purpose, we directly monitored the changes in size and concentration of initially monomeric Aβ42 samples, using Fluorescence Correlation Spectroscopy. We found that Aβ42 undergoes aggregation only when the amount of amyloid monomers exceeds the critical aggregation concentration (cac) of about 90 nM. This spontaneous, cooperative process resembles surfactants self-assembly and yields stable micelle-like oligomers whose size (≈50 monomers, R_h ≈ 7–11 nm) and elongated shape are independent of incubation time and peptide concentration. These findings reveal essential features of in vitro amyloid aggregation, which may illuminate the complex in vivo process.

Insights into the Molecular Mechanisms of Alzheimer's and Parkinson's Diseases with Molecular Simulations: Understanding the Roles of Artificial and Pathological Missense Mutations in Intrinsically Disordered Proteins Related to Pathology

Amyloid-β and α-synuclein are intrinsically disordered proteins (IDPs), which are at the center of Alzheimer's and Parkinson's disease pathologies, respectively. These IDPs are extremely flexible and do not adopt stable structures. Furthermore, both amyloid-β and α-synuclein can form toxic oligomers, amyloid fibrils and other type of aggregates in Alzheimer's and Parkinson's diseases. Experimentalists face challenges in investigating the structures and thermodynamic properties of these IDPs in their monomeric and oligomeric forms due to the rapid conformational changes, fast aggregation processes and strong solvent effects. Classical molecular dynamics simulations complement experiments and provide structural information at the atomic level with dynamics without facing the same experimental limitations. Artificial missense mutations are employed experimentally and computationally for providing insights into the structure-function relationships of amyloid-β and α-synuclein in relation to the pathologies of Alzheimer's and Parkinson's diseases. Furthermore, there are several natural genetic variations that play a role in the pathogenesis of familial cases of Alzheimer's and Parkinson's diseases, which are related to specific genetic defects inherited in dominant or recessive patterns. The present review summarizes the current understanding of monomeric and oligomeric forms of amyloid-β and α-synuclein, as well as the impacts of artificial and pathological missense mutations on the structural ensembles of these IDPs using molecular dynamics simulations. We also emphasize the recent investigations on residual secondary structure formation in dynamic conformational ensembles of amyloid-β and α-synuclein, such as β-structure linked to the oligomerization and fibrillation mechanisms related to the pathologies of Alzheimer's and Parkinson's diseases. This information represents an important foundation for the successful and efficient drug design studies.

Investigation of the interaction of amyloid β peptide (11–42) oligomers with a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) membrane using molecular dynamics simulation

The mechanism of pore formation in model neural cell membranes by β amyloid (Aβ) peptides was investigated using molecular dynamics simulation which indicated that Aβ oligomers of size equal or greater than 3 has a higher tendency for pore formation than monomers and that cholesterol tends to retard Aβ binding and insertion into the membrane.

Effect of a Paramagnetic Spin Label on the Intrinsically Disordered Peptide Ensemble of Amyloid-β

Paramagnetic relaxation enhancement is an NMR technique that has yielded important insight into the structure of folded proteins, although the perturbation introduced by the large spin probe might be thought to diminish its usefulness when applied to characterizing the structural ensembles of intrinsically disordered proteins (IDPs). We compare the computationally generated structural ensembles of the IDP amyloid-β42 (Aβ42) to an alternative sequence in which a nitroxide spin label attached to cysteine has been introduced at its N-terminus. Based on this internally consistent computational comparison, we find that the spin label does not perturb the signature population of the β-hairpin formed by residues 16-21 and 29-36 that is dominant in the Aβ42 reference ensemble. However, the presence of the tag induces a strong population shift in a subset of the original Aβ42 structural sub-populations, including a sevenfold enhancement of the β-hairpin formed by residues 27-31 and 33-38. Through back-calculation of NMR observables from the computational structural ensembles, we show that the structural differences between the labeled and unlabeled peptide would be evident in local residual dipolar couplings, and possibly differences in homonuclear ¹H-¹H nuclear Overhauser effects (NOEs) and heteronuclear ¹H-¹⁵N NOEs if the paramagnetic contribution to the longitudinal relaxation does not suppress the NOE intensities in the real experiment. This work shows that molecular simulation provides a complementary approach to resolving the potential structural perturbations introduced by reporter tags that can aid in the interpretation of paramagnetic relaxation enhancement, double electron-electron resonance, and fluorescence resonance energy transfer experiments applied to IDPs.

Comparing the Aggregation Free Energy Landscapes of Amyloid Beta(1-42) and Amyloid Beta(1-40)

Using a predictive coarse-grained protein force field, we compute and compare the free energy landscapes and relative stabilities of amyloid-β protein (1-42) and amyloid-β protein (1-40) in their monomeric and oligomeric forms up to the octamer. At the same concentration, the aggregation free energy profile of Aβ42 is more downhill, with a computed solubility that is about 10 times smaller than that of Aβ40. At a concentration of 40 μM, the clear free energy barrier between the pre-fibrillar tetramer form and the fibrillar pentamer in the Aβ40 aggregation landscape disappears for Aβ42, suggesting that the Aβ42 tetramer has a more diverse structural range. To further compare the landscapes, we develop a cluster analysis based on the structural similarity between configurations and use it to construct an oligomerization map that captures the paths of easy interconversion between different but structurally similar states of oligomers for both species. A taxonomy of the oligomer species based on β-sheet stacking topologies is proposed. The comparison of the two oligomerization maps highlights several key differences in the landscapes that can be attributed to the two additional C-terminal residues that Aβ40 lacks. In general, the two terminal residues strongly stabilize the oligomeric structures for Aβ42 relative to Aβ40, and greatly facilitate the conversion from pre-fibrillar trimers to fibrillar tetramers.

Atomic resolution conformational dynamics of intrinsically disordered proteins from NMR spin relaxation

Nuclear magnetic resonance (NMR) spectroscopy is one of the most powerful experimental approaches for investigating the conformational behaviour of intrinsically disordered proteins (IDPs). IDPs represent a significant fraction of all proteomes, and, despite their importance for understanding fundamental biological processes, the molecular basis of their activity still remains largely unknown. The functional mechanisms exploited by IDPs in their interactions with other biomolecules are defined by their intrinsic dynamic modes and associated timescales, justifying the considerable interest over recent years in the development of technologies adapted to measure and describe this behaviour. NMR spin relaxation delivers information-rich, site-specific data reporting on conformational fluctuations occurring throughout the molecule. Here we review recent progress in the use of ¹⁵N relaxation to identify local backbone dynamics and long-range chain-like motions in unfolded proteins.

How accurate are your simulations? Effects of confined aqueous volume and AMBER FF99SB and CHARMM22/CMAP force field parameters on structural ensembles of intrinsically disordered proteins: Amyloid-β42 in water

Amyloid-β₄₂ (Aβ₄₂) is an intrinsically disordered peptide intimately related to the pathogenesis of several neurodegenerative diseases. Molecular dynamics (MD) simulations are extensively utilized in the characterization of the structures and conformational dynamics of intrinsically disordered proteins (IDPs) including Aβ₄₂, with AMBER and CHARMM parameters being commonly used in these studies. Recently, comparison of the effects of force field parameters on the Aβ₄₂ structures has started to gain significant attention. In this study, the structures of Aβ₄₂ are simulated using AMBER FF99SB and CHARMM22/CMAP parameters via replica exchange MD simulations utilizing a widely used clustering algorithm. These analyses show that the structural properties (extent and positioning of the elements of secondary and tertiary structure), radius of gyration values, number and position of salt bridges are extremely dependent on the chosen force field parameters notably with the usage of clustering algorithms. For example, predicted secondary structure elements, which are of the great importance for better understanding of the molecular mechanisms of neurodegenerative diseases, deviate enormously in models generated using currently available force field parameters for proteins. Based on the derived models, chemical shift values are calculated and compared to the experimentally determined data. This comparison revealed that although both force field parameters yield results in agreement with experiments, the obtained structural properties were rather different using a clustering algorithm. In other words, these results show that the predicted structures depend heavily on the force field parameters. Importantly, since none of the force field parameters currently utilized in MD studies were developed specifically taking into account the disordered nature of IDPs, these findings clearly indicate that new force field parameters have to be developed for IDPs considering their rapid flexibility and dynamics with high amplitude. Furthermore, molecular simulations of IDPs are typically conducted using one water volume. We show that the confined aqueous volume impacts the predicted structural properties of Aβ₄₂ in water. Although up to date, confined aqueous volume effects have been ignored in the MD simulations of IDPs in water, our data indicate that these effects have to be taken into account in predicting the structural and thermodynamic properties of disordered proteins in solution.

Amyloid β Fibril Elongation by Monomers Involves Disorder at the Tip

The growth of amyloid fibrils from Aβ_1-42 peptide, one of the key pathogenic players in Alzheimer's disease, is believed to follow a nucleation-elongation mechanism. Fibril elongation is often described as a "dock-lock" procedure, where a disordered monomer adsorbs to an existing fibril in a relatively fast process (docking), followed by a slower conformational transition toward the ordered state of the template (locking). Here, we use molecular dynamics simulations of an ordered pentamer of Aβ42 at fully atomistic resolution, which includes solvent, to characterize the elongation process. We construct a Markov state model from an ensemble of short trajectories generated by an advanced sampling algorithm that efficiently diversifies a subset of the system without any bias forces. This subset corresponds to selected dihedral angles of the peptide chain at the fibril tip favored to be the fast growing one experimentally. From the network model, we extract distinct locking pathways covering time scales in the high microsecond regime. Slow steps are associated with the exchange of hydrophobic contacts, between nonnative and native intermolecular contacts as well as between intra- and intermolecular ones. The N-terminal segments, which are disordered in fibrils and typically considered inert, are able to shield the lateral interfaces of the pentamer. We conclude by discussing our findings in the context of a refined dock-lock model of Aβ fibril elongation, which involves structural disorder for more than one monomer at the growing tip.

Amyloid beta: structure, biology and structure-based therapeutic development

Amyloid beta peptide (Aβ) is produced through the proteolytic processing of a transmembrane protein, amyloid precursor protein (APP), by β- and γ-secretases. Aβ accumulation in the brain is proposed to be an early toxic event in the pathogenesis of Alzheimer's disease, which is the most common form of dementia associated with plaques and tangles in the brain. Currently, it is unclear what the physiological and pathological forms of Aβ are and by what mechanism Aβ causes dementia. Moreover, there are no efficient drugs to stop or reverse the progression of Alzheimer's disease. In this paper, we review the structures, biological functions, and neurotoxicity role of Aβ. We also discuss the potential receptors that interact with Aβ and mediate Aβ intake, clearance, and metabolism. Additionally, we summarize the therapeutic developments and recent advances of different strategies for treating Alzheimer's disease. Finally, we will report on the progress in searching for novel, potentially effective agents as well as selected promising strategies for the treatment of Alzheimer's disease. These prospects include agents acting on Aβ, its receptors and tau protein, such as small molecules, vaccines and antibodies against Aβ; inhibitors or modulators of β- and γ-secretase; Aβ-degrading proteases; tau protein inhibitors and vaccines; amyloid dyes and microRNAs.

1,2,3,4,6-penta-O-galloyl-β-D-glucopyranose Binds to the N-terminal Metal Binding Region to Inhibit Amyloid β-protein Oligomer and Fibril Formation

The early oligomerization of amyloid β-protein (Aβ) is a crucial step in the etiology of Alzheimer's disease (AD), in which soluble and highly neurotoxic oligomers are produced and accumulated inside neurons. In search of therapeutic solutions for AD treatment and prevention, potent inhibitors that remodel Aβ assembly and prevent neurotoxic oligomer formation offer a promising approach. In particular, several polyphenolic compounds have shown anti-aggregation properties and good efficacy on inhibiting oligomeric amyloid formation. 1,2,3,4,6-penta-O-galloyl-β-D-glucopyranose is a large polyphenol that has been shown to be effective at inhibiting aggregation of full-length Aβ_1-40 and Aβ_1-42, but has the opposite effect on the C-terminal fragment Aβ_25-35. Here, we use a combination of ion mobility coupled to mass spectrometry (IMS-MS), transmission electron microscopy (TEM) and molecular dynamics (MD) simulations to elucidate the inhibitory effect of PGG on aggregation of full-length Aβ_1-40 and Aβ_1-42. We show that PGG interacts strongly with these two peptides, especially in their N-terminal metal binding regions, and suppresses the formation of Aβ_1-40 tetramer and Aβ_1-42 dodecamer. By exploring multiple facets of polyphenol-amyloid interactions, we provide a molecular basis for the opposing effects of PGG on full-length Aβ and its C-terminal fragments.

Probing oligomerization of amyloid beta peptide in silico

Aggregation of amyloid β (Aβ) peptide is implicated in fatal Alzheimer's disease, for which no cure is available. Understanding the mechanisms responsible for this aggregation is required in order for therapies to be developed. In an effort to better understand the molecular mechanisms involved in spontaneous aggregation of Aβ peptide, extensive molecular dynamics simulations are reported, and the results are analyzed through a combination of structural biology tools and a novel essential collective dynamics method. Several model systems composed of ten or twelve Aβ_17-42 chains in water are investigated, and the influence of metal ions is probed. The results suggest that Aβ monomers tend to aggregate into stable globular-like oligomers with 13-23% of β-sheet content. Two stages of oligomer formation have been identified: quick collapse within the first 40 ns of the simulation, characterized by a decrease in inter-chain separation and build-up of β-sheets, and the subsequent slow relaxation of the oligomer structure. The resulting oligomers comprise a stable, coherently moving sub-aggregate of 6-9 strongly inter-correlated chains. Cu²⁺ and Fe²⁺ ions have been found to develop coordination bonds with carboxylate groups of E22, D23 and A42, which remain stable during 200 ns simulations. The presence of Fe²⁺, and particularly Cu²⁺ ions, in negatively charged cavities has been found to cause significant changes in the structure and dynamics of the oligomers. The results indicate, in particular, that formation of non-fibrillar oligomers might be involved in early template-free aggregation of Aβ_17-42 monomers, with charged species such as Cu²⁺ or Fe²⁺ ions playing an important role.

Capturing Invisible Motions in the Transition from Ground to Rare Excited States of T4 Lysozyme L99A

Proteins commonly sample a number of conformational states to carry out their biological function, often requiring transitions from the ground state to higher-energy states. Characterizing the mechanisms that guide these transitions at the atomic level promises to impact our understanding of functional protein dynamics and energy landscapes. The leucine-99-to-alanine (L99A) mutant of T4 lysozyme is a model system that has an experimentally well characterized excited sparsely populated state as well as a ground state. Despite the exhaustive study of L99A protein dynamics, the conformational changes that permit transitioning to the experimentally detected excited state (∼3%, ΔG ∼2 kcal/mol) remain unclear. Here, we describe the transitions from the ground state to this sparsely populated excited state of L99A as observed through a single molecular dynamics (MD) trajectory on the Anton supercomputer. Aside from detailing the ground-to-excited-state transition, the trajectory samples multiple metastates and an intermediate state en route to the excited state. Dynamic motions between these states enable cavity surface openings large enough to admit benzene on timescales congruent with known rates for benzene binding. Thus, these fluctuations between rare protein states provide an atomic description of the concerted motions that illuminate potential path(s) for ligand binding. These results reveal, to our knowledge, a new level of complexity in the dynamics of buried cavities and their role in creating mobile defects that affect protein dynamics and ligand binding.

Influence of Hyperglycemic Conditions on Self-Association of the Alzheimer's Amyloid β (Aβ1-42) Peptide

Clinical studies have identified a correlation between type-2 diabetes mellitus and cognitive decrements en route to the onset of Alzheimer's disease (AD). Recent studies have established that post-translational modifications of the amyloid β (Aβ) peptide occur under hyperglycemic conditions; particularly, the process of glycation exacerbates its neurotoxicity and accelerates AD progression. In view of the assertion that macromolecular crowding has an altering effect on protein self-assembly, it is crucial to characterize the effects of hyperglycemic conditions via crowding on Aβ self-assembly. Toward this purpose, fully atomistic molecular dynamics simulations were performed to study the effects of glucose crowding on Aβ dimerization, which is the smallest known neurotoxic species. The dimers formed in the glucose-crowded environment were found to have weaker associations as compared to that of those formed in water. Binding free energy calculations show that the reduced binding strength of the dimers can be mainly attributed to the overall weakening of the dispersion interactions correlated with substantial loss of interpeptide contacts in the hydrophobic patches of the Aβ units. Analysis to discern the differential solvation pattern in the glucose-crowded and pure water systems revealed that glucose molecules cluster around the protein, at a distance of 5-7 Å, which traps the water molecules in close association with the protein surface. This preferential exclusion of glucose molecules and resulting hydration of the Aβ peptides has a screening effect on the hydrophobic interactions, which in turn diminishes the binding strength of the resulting dimers. Our results imply that physical effects attributed to crowded hyperglycemic environments are incapable of solely promoting Aβ self-assembly, indicating that further mechanistic studies are required to provide insights into the self-assembly of post-translationally modified Aβ peptides, known to possess aggravated toxicity, under these conditions.

Fluorescence correlation spectroscopy reveals a cooperative unfolding of monomeric amyloid-β 42 with a low Gibbs free energy

The amyloid-beta peptide (Aβ) plays a major role in the progression of Alzheimer's disease. Due to its high toxicity, the 42 amino acid long isoform Aβ42 has become of considerable interest. The Aβ42 monomer is prone to aggregation down to the nanomolar range which makes conventional structural methods such as NMR or X-ray crystallography infeasible. Conformational information, however, will be helpful to understand the different aggregation pathways reported in the literature and will allow to identify potential conditions that favour aggregation-incompetent conformations. In this study, we applied fluorescence correlation spectroscopy (FCS) to investigate the unfolding of Alexa Fluor 488 labelled monomeric Aβ42 using guanidine hydrochloride as a denaturant. We show that our Aβ42 pre-treatment and the low-nanomolar concentrations, typically used for FCS measurements, strongly favour the presence of monomers. Our results reveal that there is an unfolding/folding behaviour of monomeric Aβ42. The existence of a cooperative unfolding curve suggests the presence of structural elements with a Gibbs free energy of unfolding of about 2.8 kcal/mol.

Prediction of ligand effects in platinum-amyloid-β coordination

Ligand field molecular mechanics (LFMM) and semi-empirical Parametric Model 7 (PM7) methods are applied to a series of six Pt^II-Ligand systems binding to the N-terminal domain of the amyloid-β (Aβ) peptide. Molecular dynamics using a combined LFMM/Assisted Model Building with Energy Refinement (AMBER) approach is used to explore the conformational freedom of the peptide fragment, and identifies favourable platinum binding modes and peptide conformations for each ligand investigated. Platinum coordination is found to depend on the nature of the ligand, providing evidence that binding mode may be controlled by suitable ligand design. Boltzmann populations at 310K indicate that each Pt-Aβ complex has a small number of thermodynamically accessible states. Ramachandran maps are constructed for the sampled Pt-Aβ conformations and secondary structural analysis of the obtained complex structures is performed and contrasted with the free peptide; coordination of these platinum complexes disrupts existing secondary structure in the Aβ peptide and promotes formation of ligand-specific turn-type secondary structure.