Molecular dynamics simulation of amyloid beta dimer formation

Oligomer Formation by Physiologically Relevant C-Terminal Isoforms of Amyloid β-Protein

Alzheimer's disease (AD) is a neurological disorder associated with amyloid β-protein (Aβ) assembly into toxic oligomers. In addition to the two predominant alloforms, Aβ1-40 and Aβ1-42, other C-terminally truncated Aβ peptides, including Aβ1-38 and Aβ1-43, are produced in the brain. Here, we use discrete molecular dynamics (DMD) and a four-bead protein model with amino acid-specific hydropathic interactions, DMD4B-HYDRA, to examine oligomer formation of Aβ1-38, Aβ1-40, Aβ1-42, and Aβ1-43. Self-assembly of 32 unstructured monomer peptides into oligomers is examined using 32 replica DMD trajectories for each of the four peptides. In a quasi-steady state, Aβ1-38 and Aβ1-40 adopt similar unimodal oligomer size distributions with a maximum at trimers, whereas Aβ1-42 and Aβ1-43 oligomer size distributions are multimodal with the dominant maximum at trimers or tetramers, and additional maxima at hexamers and unidecamers (for Aβ1-42) or octamers and pentadecamers (for Aβ1-43). The free energy landscapes reveal isoform- and oligomer-order specific structural and morphological features of oligomer ensembles. Our results show that oligomers of each of the four isoforms have unique features, with Aβ1-42 alone resulting in oligomers with disordered and solvent-exposed N-termini. Our findings help unravel the structure-function paradigm governing oligomers formed by various Aβ isoforms.

Advanced computational approaches to understand protein aggregation

Protein aggregation is a widespread phenomenon implicated in debilitating diseases like Alzheimer's, Parkinson's, and cataracts, presenting complex hurdles for the field of molecular biology. In this review, we explore the evolving realm of computational methods and bioinformatics tools that have revolutionized our comprehension of protein aggregation. Beginning with a discussion of the multifaceted challenges associated with understanding this process and emphasizing the critical need for precise predictive tools, we highlight how computational techniques have become indispensable for understanding protein aggregation. We focus on molecular simulations, notably molecular dynamics (MD) simulations, spanning from atomistic to coarse-grained levels, which have emerged as pivotal tools in unraveling the complex dynamics governing protein aggregation in diseases such as cataracts, Alzheimer's, and Parkinson's. MD simulations provide microscopic insights into protein interactions and the subtleties of aggregation pathways, with advanced techniques like replica exchange molecular dynamics, Metadynamics (MetaD), and umbrella sampling enhancing our understanding by probing intricate energy landscapes and transition states. We delve into specific applications of MD simulations, elucidating the chaperone mechanism underlying cataract formation using Markov state modeling and the intricate pathways and interactions driving the toxic aggregate formation in Alzheimer's and Parkinson's disease. Transitioning we highlight how computational techniques, including bioinformatics, sequence analysis, structural data, machine learning algorithms, and artificial intelligence have become indispensable for predicting protein aggregation propensity and locating aggregation-prone regions within protein sequences. Throughout our exploration, we underscore the symbiotic relationship between computational approaches and empirical data, which has paved the way for potential therapeutic strategies against protein aggregation-related diseases. In conclusion, this review offers a comprehensive overview of advanced computational methodologies and bioinformatics tools that have catalyzed breakthroughs in unraveling the molecular basis of protein aggregation, with significant implications for clinical interventions, standing at the intersection of computational biology and experimental research.

A turn for the worse: Aβ β-hairpins in Alzheimer's disease

Amyloid-β (Aβ) oligomers are a cause of neurodegeneration in Alzheimer's disease (AD). These soluble aggregates of the Aβ peptide have proven difficult to study due to their inherent metastability and heterogeneity. Strategies to isolate and stabilize homogenous Aβ oligomer populations have emerged such as mutations, covalent cross-linking, and protein fusions. These strategies along with molecular dynamics simulations have provided a variety of proposed structures of Aβ oligomers, many of which consist of molecules of Aβ in β-hairpin conformations. β-Hairpins are intramolecular antiparallel β-sheets composed of two β-strands connected by a loop or turn. Three decades of research suggests that Aβ peptides form several different β-hairpin conformations, some of which are building blocks of toxic Aβ oligomers. The insights from these studies are currently being used to design anti-Aβ antibodies and vaccines to treat AD. Research suggests that antibody therapies designed to target oligomeric Aβ may be more successful at treating AD than antibodies designed to target linear epitopes of Aβ or fibrillar Aβ. Aβ β-hairpins are good epitopes to use in antibody development to selectively target oligomeric Aβ. This review summarizes the research on β-hairpins in Aβ peptides and discusses the relevance of this conformation in AD pathogenesis and drug development.

Edge Substitution Effects of Histidine Tautomerization Behaviors on the Structural Properties and Aggregation Properties of Aβ(1-42) Mature Fibril

Histidine behaviors play critical roles in folding and misfolding processes due to the changes in net charge and the various N/N-H orientations on imidazole rings. However, the effect of histidine tautomerization (HIE (Nε-H, ε) and HID (Nδ-H, δ) states) behaviors on the edge chain of Aβ mature fibrils remains inadequately understood, which is critical for finding a strategy to disturb fibril elongation and growth. In the current study, eight independent molecular dynamics simulations were conducted to investigate such impacts on the structural and aggregation properties. Our results from three different binding models revealed that the binding contributions of edge substitution effects are primarily located between chains 1 and 2. Histidine states significantly influence the secondary structure of each domain. Further analysis confirmed that the C1_H6//C1_E11 intrachain interaction is essential in maintaining the internal stability of chain 1, while the C1_H13//C2_H13 and C1_H14//C2_H13 interchain interactions are critical in maintaining the interchain stability of the fibril structure. Our subsequent analysis revealed that the current edge substitution leads to the loss of the C1_H13//C1_E11 intrachain and C1_H13//C2_H14 interchain interactions. The N-terminal regularity was significantly directly influenced by histidine states, particularly by the residue of C1_H13. Our study provides valuable insights into the effect of histidine behaviors on the edge chain of Aβ mature fibril, advancing our understanding of the histidine behavior hypothesis in misfolding diseases.

Aquaporin-4 Mediated Aggregation of Alzheimer's Amyloid β-Peptide

Clearance of Alzheimer's amyloid oligomers from the brain is crucial for preventing cell toxicity. Dementia complications arise as a result of apoptosis, which is caused by peptide plaques on the lipid surface of cells. Here, we employed all-atom and coarse-grained molecular dynamics simulations to investigate the aggregation of amyloid peptides at the lipid surface and the role of aquaporin-4 (AQP4) in facilitating peptide clearance from astrocytes. The network of protein-protein interactions through text mining revealed that the expression of AQP4 and amyloid aggregation were strongly correlated. It has also been revealed that the role of aquaporins in the etiology of Alzheimer's disease involves several interconnected proteins and pathways. The nature of aggregation at the surface of the 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid bilayer was revealed by the interaction of amyloid oligomers. The membrane-bound pore region of AQP4 interacts with the peptide and slows its aggregation. This interaction maintains the helical content of the peptide while lowering its toxicity at the lipid surface. The hydrophobicity of the peptide also decreased because of these interactions, which may help in the removal of the peptide from astrocytes. Long-term coarse-grained MD simulations demonstrated different features of oligomer aggregation at the surface and strong oligomer attraction to AQP4, which inhibited aggregation. Additionally, the water dynamics of aquaporins demonstrate how the selectivity filter is broken to disrupt water flow. Our findings also provide insight into the physiological alterations in brain tissue associated with Alzheimer's disease, including water retention and increased water flow in the CSF. Furthermore, in vitro thioflavin fluorescence spectroscopy revealed a slower aggregation of the peptide in the presence of AQP4.

Big versus small: The impact of aggregate size in disease

Protein aggregation results in an array of different size soluble oligomers and larger insoluble fibrils. Insoluble fibrils were originally thought to cause neuronal cell deaths in neurodegenerative diseases due to their prevalence in tissue samples and disease models. Despite recent studies demonstrating the toxicity associated with soluble oligomers, many therapeutic strategies still focus on fibrils or consider all types of aggregates as one group. Oligomers and fibrils require different modeling and therapeutic strategies, targeting the toxic species is crucial for successful study and therapeutic development. Here, we review the role of different-size aggregates in disease, and how factors contributing to aggregation (mutations, metals, post-translational modifications, and lipid interactions) may promote oligomers opposed to fibrils. We review two different computational modeling strategies (molecular dynamics and kinetic modeling) and how they are used to model both oligomers and fibrils. Finally, we outline the current therapeutic strategies targeting aggregating proteins and their strengths and weaknesses for targeting oligomers versus fibrils. Altogether, we aim to highlight the importance of distinguishing the difference between oligomers and fibrils and determining which species is toxic when modeling and creating therapeutics for protein aggregation in disease.

Inhibition Effect and Molecular Mechanisms of Quercetin on the Aβ42 Dimer: A Molecular Dynamics Simulation Study

Amyloid-β (Aβ) dimer as the smallest oligomer has recently been drawing attention due to its neurotoxicity, transient nature, and heterogeneity. The inhibition of Aβ dimer's aggregation is the key to primary intervention of Alzheimer's disease. Previous experimental studies have reported that quercetin, the widespread polyphenolic constituent of multiple fruits and vegetables, can hamper the formation of Aβ protofibrils and disaggregate Aβ fibrils. However, the molecular mechanisms of quercetin in the suppression of the Aβ(1-42) dimer's conformational changes still remain elusive. In this work, to investigate the inhibitory mechanisms of quercetin molecules on the Aβ(1-42) dimer, an Aβ(1-42) dimer based on monomeric the Aβ(1-42) peptide with enriched coil structures is constructed. The early molecular mechanisms of quercetin molecules on inhibiting the Aβ(1-42) dimer at two different Aβ42-to-quercetin molar ratios (1:5 and 1:10) are explored via all-atom molecular dynamics simulations. The results indicate that quercetin molecules can impede the configurational change of the Aβ(1-42) dimer. The interactions and the binding affinity between the Aβ(1-42) dimer and quercetin molecules in the Aβ42 dimer + 20 quercetin system are stronger in comparison with that in the Aβ42 dimer + 10 quercetin system. Our work may be helpful in developing new drug candidates for preventing the conformational transition and further aggregation of the Aβ dimer.

Modulation of Mechanical Properties of Short Bioinspired Peptide Materials by Single Amino-Acid Mutations

The occurrence of modular peptide repeats in load-bearing (structural) proteins is common in nature, with distinctive peptide sequences that often remain conserved across different phylogenetic lineages. These highly conserved peptide sequences endow specific mechanical properties to the material, such as toughness or elasticity. Here, using bioinformatic tools and phylogenetic analysis, we have identified the GX8 peptide with the sequence GLYGGYGX (where X can be any residue) in a wide range of organisms. By simple mutation of the X residue, we demonstrate that GX8 can be self-assembled into various supramolecular structures, exhibiting vastly different physicochemical and viscoelastic properties, from liquid-like coacervate microdroplets to hydrogels to stiff solid materials. A combination of spectroscopic, electron microscopy, mechanical, and molecular dynamics studies is employed to obtain insights into molecular scale interactions driving self-assembly of GX8 peptides, underscoring that π-π stacking and hydrophobic interactions are the drivers of peptide self-assembly, whereas the X residue determines the extent of hydrogen bonding that regulates the macroscopic mechanical response. This study highlights the ability of single amino-acid polymorphism to tune the supramolecular assembly and bulk material properties of GX8 peptides, enabling us to cover a broad range of potential biomedical applications such as hydrogels for tissue engineering or coacervates for drug delivery.

Structural properties of Aβ (1-40) peptide in protonation stage of one, two, and three: New insights from the histidine protonation behaviors

Structural properties and aggregation tendency can be significantly influenced by histidine behaviors (histidine on N^δ-H state is defined as δ, likewise, N^ε-H: ε and both N^δ-H and N^ε-H: p). In current study, we investigated structural properties of Aβ(1-40) peptide during protonation evolution stage of one, two, and three by total 19 independent replica exchange molecular dynamics simulations using implicit solvent. Our results show that any kind of protonated state will promote β-sheet structure formation in comparison with deprotonated (εεε). With increase in number of protonation, the lowest β-sheet content increased. The highest averaged β-sheet structure content was detected in (δpδ) (46.0 %), (εpp) (36.8 %), and (ppp) (16.0 %) in each protonation stage. With three β-strand structures, (δpδ) shows more stable features and high hydrophobic properties. Further analysis confirmed that H13 and H14 are more important than H6. Specifically, H13 and H14 have a synergistic effect for structural formations by controlling H-bond networks in H13(p) with V39/V40 and H14(p/δ) with G37/G38. Finally, the Pearson correlation coefficient results confirmed that experimental result (ref. 44) is corresponding to our (εpp) system. Our current study will be conducive to understanding the effects of the histidine behaviors, it provides new insights for exploration protein folding and misfolding processes.

Design and Study of In Silico Binding Dynamics of Certain Isoxazole Bearing Leads Against Aβ-42 and BACE-1 Loop in Protein Fibrillation

Aims:

Design isoxazole bearing leads as dual inhibitors against Amyloid β and BACE-1 loop in protein fibrillation.

Background:

Protein fibrillation is one of the key reasons for several diseases, namely Alzheimer’s, Parkinson’s, and many others. One of the key strategies of preventing protein fibrillation is destabilizing the protein fibrils themselves or inhibiting the amyloid fibril-forming pathway in the initial stage.

Introduction:

Attempts have been taken to design newer leads to inhibit protein fibrillation by targeting the β-amyloidogenesis pathway in the brain. To exploit interfenestration between Amyloid β -42 protein and BACE-1 (β-site amyloid precursor protein cleaving enzyme) for amyloidogenesis, studies are undertaken to design dual inhibitors against the same.

Method:

In vitro binding interactions were found using docking, de novo ligand design, and MD simulation study.

Results:

Three compounds bearing an isoxazole heterocyclic nucleus were designed which could successfully bind to the hydrophobic raft and salt bridge residues Asp 23-Lys-26 of Amyloid β, destabilizing the growing fibril. Additionally, one of our candidate compounds exhibited force of interaction with Thr232 at the S3 pocket of BACE-1, interacted with key residue Asp228, Tyr71, and Thr72 of the β-hairpin flap and hydrogen bonding with Gly11 at loop 10s.

Conclusion:

Protein flexibility dynamics of the Aβ-42 protein revealed that there is a considerable conformational change of the same with or without ligand binding. The lower RMSF of the bound region and reprogramming residual contacts within the Aβ-42 protein suggested successful binding of the ligand with the protein, lowering the access for further β-β dimerization.

Amyloid-β peptide dimers undergo a random coil to β-sheet transition in the aqueous phase but not at the neuronal membrane

Mounting evidence suggests that the neuronal cell membrane is the main site of oligomer-mediated neuronal toxicity of amyloid-β peptides in Alzheimer's disease. To gain a detailed understanding of the mutual interference of amyloid-β oligomers and the neuronal membrane, we carried out microseconds of all-atom molecular dynamics (MD) simulations on the dimerization of amyloid-β (Aβ)42 in the aqueous phase and in the presence of a lipid bilayer mimicking the in vivo composition of neuronal membranes. The dimerization in solution is characterized by a random coil to β-sheet transition that seems on pathway to amyloid aggregation, while the interactions with the neuronal membrane decrease the order of the Aβ42 dimer by attenuating its propensity to form a β-sheet structure. The main lipid interaction partners of Aβ42 are the surface-exposed sugar groups of the gangliosides GM1. As the neurotoxic activity of amyloid oligomers increases with oligomer order, these results suggest that GM1 is neuroprotective against Aβ-mediated toxicity.

Effects of Aβ-derived peptide fragments on fibrillogenesis of Aβ

Amyloid β (Aβ) peptide aggregation plays a central role in Alzheimer's disease (AD) etiology. AD drug candidates have included small molecules or peptides directed towards inhibition of Aβ fibrillogenesis. Although some Aβ-derived peptide fragments suppress Aβ fibril growth, comprehensive analysis of inhibitory potencies of peptide fragments along the whole Aβ sequence has not been reported. The aim of this work is (a) to identify the region(s) of Aβ with highest propensities for aggregation and (b) to use those fragments to inhibit Aβ fibrillogenesis. Structural and aggregation properties of the parent Aβ_1-42 peptide and seven overlapping peptide fragments have been studied, i.e. Aβ_1-10 (P1), Aβ_6-15 (P2), Aβ_11-20 (P3), Aβ_16-25 (P4), Aβ_21-30 (P5), Aβ_26-36 (P6), and Aβ_31-42 (P7). Structural transitions of the peptides in aqueous buffer have been monitored by circular dichroism and Fourier transform infrared spectroscopy. Aggregation and fibrillogenesis were analyzed by light scattering and thioflavin-T fluorescence. The mode of peptide-peptide interactions was characterized by fluorescence resonance energy transfer. Three peptide fragments, P3, P6, and P7, exhibited exceptionally high propensity for β-sheet formation and aggregation. Remarkably, only P3 and P6 exerted strong inhibitory effect on the aggregation of Aβ_1-42, whereas P7 and P2 displayed moderate inhibitory potency. It is proposed that P3 and P6 intercalate between Aβ_1-42 molecules and thereby inhibit Aβ_1-42 aggregation. These findings may facilitate therapeutic strategies of inhibition of Aβ fibrillogenesis by Aβ-derived peptides.

Effects of applied surface-tension on membrane-assisted Aβ aggregation

Accumulation of protein-based (Aβ) aggregates on cellular membranes with varying structural properties is commonly recognized as the key step in Alzheimer's pathogenesis. But experimental and computational challenges have made this biophysical characterization difficult. In particular, studies connecting biological membrane organization and Aβ aggregation are limited. While experiments have suggested that an increased membrane curvature results in faster Aβ peptide aggregation in the context of Alzheimer's disease, a mechanistic explanation for this relation is missing. In this work, we are leveraging molecular simulations with a physics-based coarse grained model to address and understand the relationships between curved cellular membranes and aggregation of a model template peptide Aβ 16-22. In agreement with experimental results, our simulations also suggest a positive correlation between increased peptide aggregation and membrane curvature. More curved membranes have higher lipid packing defects that engage peptide hydrophobic groups and promote faster diffusion leading to peptide fibrillar structures. In addition, we curated the effects of peptide aggregation on the membrane's structure and organization. Interfacial peptide aggregation results in heterogeneous headgroup-peptide interactions and an induced crowding effect at the lipid headgroup region, leading to a more ordered headgroup region and disordered lipid-tails at the membrane core. This work presents a mechanistic and morphological overview of the relationships between the biomembrane local structure and organization, and Aβ peptide aggregation.

Effect of charged mutation on aggregation of a pentapeptide: Insights from molecular dynamics simulations

Aggregation of therapeutic monoclonal antibodies (mAbs) can negatively affect their chemistry, manufacturing, and control attributes and lead to undesirable immune responses in patients. Therefore, optimization of lead mAb drug candidates during discovery stages to mitigate aggregation is increasingly becoming an integral part of their developability assessments. The disruption of short sequence motifs called aggregation prone regions (APRs) found in amino acid sequences of mAb candidates can potentially mitigate their aggregation. In this work, we have performed molecular dynamics simulations to study the aggregation of an APR (VLVIY) found in λ light chains of human antibodies and its single point mutant KLVIY. Eighteen different multicopy peptide simulation systems of "VLVIY" and "KLVIY" were constructed by varying their concentrations, temperatures, termini capping, and flanking gate-keeper regions. Within 20 ns of the simulation, peptide "VLVIY" formed an aggregate of 100 peptides at ~0.1 M concentration with a 60% reduction in solvent accessible surface area (SASA). Furthermore, analysis of the SASA change, peptide cluster distribution, and water residence time demonstrated how Val ➔ Lys mutation resists aggregation and improves solubility. Presence of Lys slows down aggregation kinetics via charge-charge repulsions and by raising the kinetic barrier to formation of large oligomers. However, the effect of the Val ➔ Lys mutation is dependent on sequence and structural contexts around the APR. This mutation also alters the solvation shell around the peptide by favoring solute-solvent interactions, thereby increasing its solubility. This work has provided a detailed mechanistic explanation of how APR disruption can mitigate aggregation in biotherapeutics and improve their developability.

Molecular insight into the early stage of amyloid-β(1-42) Homodimers aggregation influenced by histidine tautomerism

Aggregated amyloid β-peptide (Aβ) in small oligomeric forms inside the brain causes synaptic function disruption and the development of Alzheimer's disease (AD). Histidine is an important amino acid that may lead to structural changes. Aβ42 monomer chain includes 3 histidine residues that considering two ε and δ tautomers 8 isomers, including (εεε) and (εδδ) could be formed. Molecular dynamics simulation on homodimerization of (εεε) (the most common type of tautomers) and (εδδ) tautomers with different initial configurations using monomer chains from our previous work were performed to uncover the tautomeric behavior of histidine on Aβ42 aggregation in a physiological pH which is still largely unknown and impossible to observe experimentally. We found a higher propensity of forming β-sheet in (εδδ) homodimers and specifically in a greater amount from Aβ42 than from Aβ40. A smaller amount of β-sheet formation was observed for (εεε) homodimers compared with (εδδ). Additionally, interactions in (εδδ) homodimers may indicate the importance of the hydrophobic core and C-/N-terminals during oligomerization. Our findings indicate the important role of the tautomeric effect of histidine and (εδδ) homodimers at the early stage of Aβ aggregation.

Defining the Neuropathological Aggresome across in Silico, in Vitro, and ex Vivo Experiments

The loss of proteostasis over the life course is associated with a wide range of debilitating degenerative diseases and is a central hallmark of human aging. When left unchecked, proteins that are intrinsically disordered can pathologically aggregate into highly ordered fibrils, plaques, and tangles (termed amyloids), which are associated with countless disorders such as Alzheimer's disease, Parkinson's disease, type II diabetes, cancer, and even certain viral infections. However, despite significant advances in protein folding and solution biophysics techniques, determining the molecular cause of these conditions in humans has remained elusive. This has been due, in part, to recent discoveries showing that soluble protein oligomers, not insoluble fibrils or plaques, drive the majority of pathological processes. This has subsequently led researchers to focus instead on heterogeneous and often promiscuous protein oligomers. Unfortunately, significant gaps remain in how to prepare, model, experimentally corroborate, and extract amyloid oligomers relevant to human disease in a systematic manner. This Review will report on each of these techniques and their successes and shortcomings in an attempt to standardize comparisons between protein oligomers across disciplines, especially in the context of neurodegeneration. By standardizing multiple techniques and identifying their common overlap, a clearer picture of the soluble neuropathological aggresome can be constructed and used as a baseline for studying human disease and aging.

Mutual structural effects of unmodified and pyroglutamylated amyloid β peptides during aggregation

Amyloid β (Aβ) peptide aggregates are linked to Alzheimer's disease (AD). Posttranslationally pyroglutamylated Aβ (pEAβ) occurs in AD brains in significant quantities and is hypertoxic, but the underlying structural and aggregation properties remain poorly understood. Here, the structure and aggregation of Aβ_1-40 and pEAβ_3-40 are analyzed separately and in equimolar combination. Circular dichroism data show that Aβ_1-40 , pEAβ_3-40 , and their combination assume α-helical structure in dry state and transition to unordered structure in aqueous buffer. Aβ_1-40 and the 1:1 combination gradually acquire β-sheet structure while pEAβ_3-40 adopts an α-helix/β-sheet conformation. Thioflavin-T fluorescence studies suggest that the two peptides mutually inhibit fibrillogenesis. Fourier transform infrared (FTIR) spectroscopy identifies the presence of β-turn and α-helical structures in addition to β-sheet structure in peptides in aqueous buffer. The kinetics of transitions from the initial α-helical structure to β-sheet structure were resolved by slow hydration of dry peptides by D₂ O vapor, coupled with isotope-edited FTIR. These data confirmed the mutual suppression of β-sheet formation by the two peptides. Remarkably, pEAβ_3-40 maintained a significant fraction of α-helical structure in the combined sample, implying a reduced β-sheet propensity of pEAβ_3-40 . Altogether, the data imply that the combination of unmodified and pyroglutamylated Aβ peptides resists fibrillogenesis and favors the prefibrillar state, which may underlie hypertoxicity of pEAβ.

The protective mutation A673T in amyloid precursor protein gene decreases Aβ peptides production for 14 forms of Familial Alzheimer's Disease in SH-SY5Y cells

The deposition of Aβ plaques in the brain leads to the onset and development of Alzheimer’s disease. The Amyloid precursor protein (APP) is cleaved by α-secretase (non-amyloidogenic processing of APP), however increased cleavage by β-secretase (BACE1) leads to the accumulation of Aβ peptides, which forms plaques. APP mutations mapping to exons 16 and 17 favor plaque accumulation and cause Familial Alzheimer Disease (FAD). However, a variant of the APP gene (A673T) originally found in an Icelandic population reduces BACE1 cleavage by 40%. A series of plasmids containing the APP gene, each with one of 29 different FAD mutations mapping to exon 16 and exon 17 was created. These plasmids were then replicated with the addition of the A673T mutation. Combined these formed the library of plasmids that was used in this study. The plasmids were transfected in neuroblastomas to assess the effect of this mutation on Aβ peptide production. The production of Aβ peptides was decreased for some FAD mutations due to the presence of the co-dominant A673T mutation. The reduction of Aβ peptide concentrations for the London mutation (V717I) even reached the same level as for A673T control in SH-SY5Y cells. These preliminary results suggest that the insertion of A673T in APP genes containing FAD mutations might confer a clinical benefit in preventing or delaying the onset of some FADs.

A lowly populated, transient β-sheet structure in monomeric Aβ1-42 identified by multinuclear NMR of chemical denaturation

Chemical denaturation is a well-established approach for probing the equilibrium between folded and unfolded states of proteins. We demonstrate applicability of this method to the detection of a small population of a transiently folded structural element in a system that is often considered to be intrinsically fully disordered. The ¹H^N, ¹⁵N, ¹³C^α, and ¹³C' chemical shifts of Aβ^1-40 and Aβ^1-42 peptides and their M35-oxidized variants were monitored as a function of urea concentration and compared to analogous urea titrations of synthetic pentapeptides of homologous sequence. Fitting of the chemical shift titrations yields a 10 ± 1% population for a structured element at the C-terminus of Aβ^1-42 that folds with a cooperativity of m = 0.06 kcal/mol·M. The fit also yields the chemical shifts of the folded state and, using a database search, for Aβ^1-42 these shifts identified an antiparallel intramolecular β-sheet for residues I32-A42, linked by a type I' β-turn at G37 and G38. The structure is destabilized by oxidation of M35. Paramagnetic relaxation rates and two previously reported weak, medium-range NOE interactions are consistent with this transient β-sheet. Introduction of the requisite A42C mutation and tagging with MTSL resulted in a small stabilization of this β-sheet. Chemical shift analysis suggests a C-terminal β-sheet may be present in Aβ^1-40 too, but the turn type at G37 is not type I'. The approach to derive Transient Structure from chemical Denaturation by NMR (TSD-NMR), demonstrated here for Aβ peptides, provides a sensitive tool for identifying the presence of lowly populated, transiently ordered elements in proteins that are considered to be intrinsically disordered, and permits extraction of structural data for such elements.

Effect of Amyloid-β Monomers on Lipid Membrane Mechanical Parameters-Potential Implications for Mechanically Driven Neurodegeneration in Alzheimer's Disease

Alzheimer’s disease (AD) is a neurodegenerative disease that results in memory loss and the impairment of cognitive skills. Several mechanisms of AD’s pathogenesis were proposed, such as the progressive accumulation of amyloid-β (Aβ) and τ pathology. Nevertheless, the exact neurodegenerative mechanism of the Aβ remains complex and not fully understood. This paper proposes an alternative hypothesis of the mechanism based on maintaining the neuron membrane’s mechanical balance. The incorporation of Aβ decreases the lipid membrane’s elastic properties, which eventually leads to the impairment of membrane clustering, disruption of mechanical wave propagation, and change in gamma oscillations. The first two disrupt the neuron’s ability to function correctly while the last one decreases sensory encoding and perception enabling. To begin discussing this mechanical-balance hypothesis, we measured the effect of two selected peptides, Aβ-40 and Aβ-42, as well as their fluorescently labeled modification, on membrane mechanical properties. The decrease of bending rigidity, consistent for all investigated peptides, was observed using molecular dynamic studies and experimental flicker-noise techniques. Additionally, wave propagation was investigated with molecular dynamic studies in membranes with and without incorporated neurodegenerative peptides. A change in membrane behavior was observed in the membrane system with incorporated Aβ.

Edge Strand Dissociation and Conformational Changes in Transthyretin under Amyloidogenic Conditions

During amyloidogenesis, proteins undergo conformational changes that allow them to aggregate and assemble into insoluble, fibrillar structures. Soluble oligomers that form during this process typically contain 2-24 monomeric subunits and are cytotoxic. Before the formation of these soluble oligomers, monomeric species first adopt aggregation-competent conformations. Knowledge of the structures of these intermediate states is invaluable to the development of molecular strategies to arrest pathological amyloid aggregation. However, the highly dynamic and interconverting nature of amyloidogenic species limits biophysical characterization of their structures during amyloidogenesis. Here, we use molecular dynamics simulations to probe conformations sampled by monomeric transthyretin under amyloidogenic conditions. We show that certain β-strands in transthyretin tend to unfold and sample nonnative conformations and that the edge strands in one β-sheet (the DAGH sheet) are particularly susceptible to conformational changes in the monomeric state. We also find that changes in the tertiary structure of transthyretin can be associated with disruptions to the secondary structure. We evaluated the conformations produced by molecular dynamics by calculating how well molecular-dynamics-derived structures reproduced NMR-derived interatomic distances. Finally, we leverage our computational results to produce experimentally testable hypotheses that may aid experimental explorations of pathological conformations of transthyretin.

Differences in the free energies between the excited states of Aβ40 and Aβ42 monomers encode their aggregation propensities

The early events in the aggregation of the intrinsically disordered peptide, amyloid-β (Aβ), involve transitions from the disordered free energy ground state to assembly-competent states. Are the fingerprints of order found in the amyloid fibrils encoded in the conformations that the monomers access at equilibrium? If so, could the enhanced aggregation rate of Aβ42 compared to Aβ40 be rationalized from the sparsely populated high free energy states of the monomers? Here, we answer these questions in the affirmative using coarse-grained simulations of the self-organized polymer-intrinsically disordered protein (SOP-IDP) model of Aβ40 and Aβ42. Although both the peptides have practically identical ensemble-averaged properties, characteristic of random coils (RCs), the conformational ensembles of the two monomers exhibit sequence-specific heterogeneity. Hierarchical clustering of conformations reveals that both the peptides populate high free energy aggregation-prone ([Formula: see text]) states, which resemble the monomers in the fibril structure. The free energy gap between the ground (RC) and the [Formula: see text] states of Aβ42 peptide is smaller than that for Aβ40. By relating the populations of excited states of the two peptides to the fibril formation time scales using an empirical formula, we explain nearly quantitatively the faster aggregation rate of Aβ42 relative to Aβ40. The [Formula: see text] concept accounts for fibril polymorphs, leading to the prediction that the less stable [Formula: see text] state of Aβ42, encoding for the U-bend fibril, should form earlier than the structure with the S-bend topology, which is in accord with Ostwald's rule rationalizing crystal polymorph formation.

Aβ40 Oligomers Promote Survival and Early Neuronal Differentiation of Dentate Gyrus-Isolated Precursor Cells Through Activation of the Akt Signaling Pathway

The amyloid beta-peptide (Aβ) is the low-abundance product of amyloid precursor protein (APP), which is produced lifelong in the healthy brain. The functional properties of Aβ40 and Aβ42 peptides have not been completely elucidated to date. Although, several studies suggest that these peptides have a number of neurotrophic and neurotoxic properties, respectively. Interestingly, low concentrations of Aβ40 and Aβ42 regulate synaptic plasticity and improve cognitive functions, whereas the accumulation of Aβ42, coupled with the effects of age, can cause dysregulation of synaptic function, as is shown in Alzheimer's disease. Additionally, several studies suggest that both peptides, Aβ40 and Aβ42, are involved in neurogenic processes; however, these results are still controversial. Moreover, existing data indicate a direct relationship between the physicochemical characteristics of the peptides and their effects. Herein, we evaluated the effect of Aβ40 oligomers on hippocampal precursor cells isolated from the dentate gyrus of adult female C57Bl6 mice (mADGPCs). To this end, mADGPCs were treated with nanomolar and micromolar range concentrations of oligomeric forms of Aβ40 for 24, 48, and 72 h to evaluate their effects on several events in the neurogenic process in vitro, including viability, proliferation, and early differentiation. The results indicate that Aβ40 favors mADGPC proliferation, survival, and neuronal differentiation following a mechanism that involves activation of the Akt signaling pathway. Thus, this study provides evidence about the positive effects of Aβ40 oligomers on the neurogenic process in adult mouse hippocampal precursor cells in vitro.

Utilisation of the OliveNet™ Library to investigate phenolic compounds using molecular modelling studies in the context of Alzheimer's disease

Alzheimer's disease (AD) is a debilitating neurodegenerative disease that affects over 47 million people worldwide, and is the most common form of dementia. There is a vast body of literature demonstrating that the disease is caused by an accumulation of toxic extracellular amyloid-β (Aβ) peptides and intracellular neurofibrillary tangles that consist of hyperphosphorylated tau. Adherence to the Mediterranean diet has been shown to reduce the incidence of AD and the phenolic compounds in extra virgin olive oil, including oleocanthal, have gained a significant amount of attention. A large number of these ligands have been described in the pre-existing literature and 222 of these compounds have been characterised in the OliveNet™ database. In this study, molecular docking was used to screen the 222 phenolic compounds from the OliveNet™ database and assess their ability to bind to various forms of the Aβ and tau proteins. The phenolic ligands were found to be binding strongly to the hairpin-turn of the Aβ1-40 and Aβ1-42 monomers, and binding sites were also identified in the tau fibril protein structures. Luteolin-4'-O-rutinoside, oleuricine A, isorhoifolin, luteolin-7-O-rutinoside, cyanidin-3-O-rutinoside and luteolin-7,4-O-diglucoside were predicted to be novel lead compounds. Molecular dynamics (MD) simulations performed using well-known olive ligands bound to Aβ1-42 oligomers highlighted that future work may examine potential anti-aggregating properties of novel compounds in the OliveNet™ database. This may lead to the development and evaluation of new compounds that may have efficacy against Alzheimer's disease.

High-resolution probing of early events in amyloid-β aggregation related to Alzheimer's disease

In Alzheimer's disease (AD), soluble oligomers of amyloid-β (Aβ) are emerging as a crucial entity in driving disease progression as compared to insoluble amyloid deposits. The lacuna in establishing the structure to function relationship for Aβ oligomers prevents the development of an effective treatment for AD. While the transient and heterogeneous properties of Aβ oligomers impose many challenges for structural investigation, an effective use of a combination of NMR techniques has successfully identified and characterized them at atomic-resolution. Here, we review the successful utilization of solution and solid-state NMR techniques to probe the aggregation and structures of small and large oligomers of Aβ. Biophysical studies utilizing the commonly used solution and 19F based NMR experiments to identify the formation of small size early intermediates and to obtain their structures, and dock-lock mechanism of fiber growth at atomic-resolution are discussed. In addition, the use of proton-detected magic angle spinning (MAS) solid-state NMR experiments to obtain high-resolution insights into the aggregation pathways and structures of large oligomers and other aggregates is also presented. We expect these NMR based studies to be valuable for real-time monitoring of the depletion of monomers and the formation of toxic oligomers and high-order aggregates under a variety of conditions, and to solve the high-resolution structures of small and large size oligomers for most amyloid proteins, and therefore to develop inhibitors and drugs.

Acetylation of Aβ42 at Lysine 16 Disrupts Amyloid Formation

The residue lysine 28 (K28) is known to form an important salt bridge that stabilizes the Aβ amyloid structure, and acetylation of lysine 28 (K28Ac) slows the Aβ42 fibrillization rate but does not affect fibril morphology. On the other hand, acetylation of lysine 16 (K16Ac) residue greatly diminishes the fibrillization property of Aβ42 peptide and also affects its toxicity. This is due to the fact that lysine 16 acetylated amyloid beta peptide forms amorphous aggregates instead of amyloid fibrils. This is likely a result of increased hydrophobicity of the K16-A21 region due to K16 acetylation, as confirmed by molecular dynamic simulation studies. The calculated results show that the hydrophobic patches of aggregates from acetylated peptides were different when compared to wild-type (WT) peptide. K16Ac and double acetylated (KKAc) peptide aggregates show significantly higher cytotoxicity compared to the WT or K28Ac peptide aggregates alone. However, the heterogeneous mixture of WT and acetylated Aβ42 peptide aggregates exhibited higher free radical formation as well as cytotoxicity, suggesting dynamic interactions between different species could be a critical contributor to Aβ pathology.

Oligomerization and Conformational Change Turn Monomeric β-Amyloid and Tau Proteins Toxic: Their Role in Alzheimer's Pathogenesis

The structural polymorphism and the physiological and pathophysiological roles of two important proteins, β-amyloid (Aβ) and tau, that play a key role in Alzheimer's disease (AD) are reviewed. Recent results demonstrate that monomeric Aβ has important physiological functions. Toxic oligomeric Aβ assemblies (AβOs) may play a decisive role in AD pathogenesis. The polymorph fibrillar Aβ (fAβ) form has a very ordered cross-β structure and is assumed to be non-toxic. Tau monomers also have several important physiological actions; however, their oligomerization leads to toxic oligomers (TauOs). Further polymerization results in probably non-toxic fibrillar structures, among others neurofibrillary tangles (NFTs). Their structure was determined by cryo-electron microscopy at atomic level. Both AβOs and TauOs may initiate neurodegenerative processes, and their interactions and crosstalk determine the pathophysiological changes in AD. TauOs (perhaps also AβO) have prionoid character, and they may be responsible for cell-to-cell spreading of the disease. Both extra- and intracellular AβOs and TauOs (and not the previously hypothesized amyloid plaques and NFTs) may represent the novel targets of AD drug research.

In Situ Observation of Amyloid Nucleation and Fibrillation by FastScan Atomic Force Microscopy

Amyloidogenic proteins are key components in various amyloid diseases. The aggregation process and the local structural changes of the toxic species from toxic oligomers to protofibrils and subsequently to mature fibrils are crucial for understanding the molecular mechanism of the amyloidgenic process and also for developing a treatment strategy. Exploration on amyloid aggregation dynamics in situ under real liquid condition is feasible for reflection of the whole process with biological correlations. Herein we report the in situ dynamic study and structure exploration of Amylin_1-37 aggregation by FastScan atomic force microscopy. Amylin_1-37 nucleation process was observed in which smaller oligomers or monomers were assimilated by the surrounding big oligomers. Amylin_1-37 protofibril aggregation was positively correlated with monomer concentration, whereas no direct relationship was observed between fibril elongation and monomer concentration. Growing end and passivated end were found during Amylin_1-37 fibrillation. In the assembly process, the growing end kept its structure, and its stiffness was lower than the aggregate body, whereas the passivated end might experience rearrangements of β-structures, which eventually enabled fibril growth from this end. This work is beneficial to the insights of amyloid fibrillation and may shed light on the development of drugs targeting the specific phase of amyloid aggregation.

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.

Structural discordance in HIV-1 Vpu from brain isolate alarms amyloid fibril forming behavior- a computational perspective

HIV-1 being the most widespread type worldwide, its accounts for almost 95% of all infections including HIV associated dementia (HAD) that triggers neurological dysfunction and neurodegeneration in patients. The common features associated with HAD and other neurodegenerative diseases are accumulation of amyloid plaques, neuronal loss and deterioration of cognitive abilities, amongst which amyloid fibrillation is considered to be a hallmark. The success of effective therapeutics lies in the understanding of mechanisms leading to neurotoxicity. Few viral proteins like gp-120 are known to be involved in aggregation and enhancement of viral infectivity while comprehending the neurotoxic role of some other proteins is still underway. In the current study, amyloidogenic potential of HIV-1 Vpu protein from brain isolate is investigated through computational approaches. The aggregation propensity of brain derived HIV-1 Vpu was assessed by several amyloid prediction servers that projected the region 4-35 to be amyloidogenic. The protein structure was modeled and subjected to 70 ns molecular dynamics (MD) simulation to investigate the transformation of α-helical conformation of the predicted aggregate region into β-sheet, proposing the protein's ability to initiate fibril formation that is central to amyloidogenic proteins. The structural features of brain derived HIV-1 Vpu were consistent with the in silico amyloid prediction results that depicts the conformational change in the region 8-28 of which residues Ala8, Ile9, Val10, Ala19, Ile20 and Val21 constitutes β-sheet formation. The α-helix/β-sheet discordance of the predicted region was reflected in the simulation study highlighting the possible structural transition associated with HIV-1 Vpu protein of brain isolate.

Looking at the Disordered Proteins through the Computational Microscope

Intrinsically disordered proteins (IDPs) have attracted wide interest over the past decade due to their surprising prevalence in the proteome and versatile roles in cell physiology and pathology. A large selection of IDPs has been identified as potential targets for therapeutic intervention. Characterizing the structure-function relationship of disordered proteins is therefore an essential but daunting task, as these proteins can adapt transient structure, necessitating a new paradigm for connecting structural disorder to function. Molecular simulation has emerged as a natural complement to experiments for atomic-level characterizations and mechanistic investigations of this intriguing class of proteins. The diverse range of length and time scales involved in IDP function requires performing simulations at multiple levels of resolution. In this Outlook, we focus on summarizing available simulation methods, along with a few interesting example applications. We also provide an outlook on how these simulation methods can be further improved in order to provide a more accurate description of IDP structure, binding, and assembly.

β- N-Methylamino-l-alanine (BMAA) Not Involved in Alzheimer's Disease

β- N-Methylamino-l-alanine (BMAA) is a neurotoxic agent implicated in ALS as well as Parkinson's and Alzheimer's diseases. It is produced by blue-green algae and could find its way via fish and seafood into the human food supply. Isolation from biological samples yields the compound in monomeric and protein-bound form. It has been suggested that the protein-bound fraction may result from genetic misincorporation into proteins in place of serine. Concomitant misfolding of the mutated proteins may be responsible for the neurological diseases. Recent reports that contradict the misincorporation theory leave unresolved the nature of the protein-bound form of BMAA. We have found from quantum mechanical calculations on model systems that it is possible to bind BMAA with high affinity in a noncovalent fashion to proteins. Because of our interest in Alzheimer's disease, molecular dynamics simulations were applied to search for such binding between BMAA and the β-amyloid peptide and to discover whether replacement of either of its two serine residues could affect its aggregation into neurotoxic oligomers. No stable noncovalently bound complex was found, and it was concluded that incorporation of BMAA in place of serine would not alter the conformational dynamics of the β-amyloid peptide.

Bending energy penalty enhances the adhesive strength of functional amyloid curli to surfaces

The functional amyloid curli fiber, a major proteinaceous component of biofilm extracellular matrices, plays an important role in biofilm formation and enterobacteriaceae adhesion. Curli nanofibers exhibit exceptional underwater adhesion to various surfaces, have high rigidity and strong tensile mechanical properties, and thus hold great promise in biomaterials. The mechanisms of how curli fibers strongly attach to surfaces and detach under force remain elusive. To investigate curli fiber adhesion to surfaces, we developed a coarse-grained curli fiber model, in which the protein subunit CsgA (curli specific gene A) self-assembles into the fiber. The coarse-grained model yields physiologically relevant and tunable bending rigidity and persistence length. The force-induced desorption of a single curli fiber is examined using coarse-grained modeling and theoretical analysis. We find that the bending energy penalty arising from high persistence length enhances the resistance of the curli fiber against desorption and thus strengthens the adhesion of the curli fiber to surfaces. The CsgA-surface adhesion energy and the curli fiber bending rigidity both play crucial roles in the resistance of curli fiber against desorption from surfaces. To enable the desorption process, the applied peeling force must overcome both the interfacial adhesion energy and the energy barrier for bending the curli fiber at the peeling front. We show that the energy barrier to desorption increases with the interfacial adhesion energy, however, the bending induced failure of a single curli fiber limits the work of adhesion if the proportion of the CsgA-surface adhesion energy to the CsgA-CsgA cohesive energy becomes large. These results illustrate that the optimal adhesion performance of nanofibers is dictated by the interplay between bending, surface energy and cohesive energy. Our model provides timely insight into enterobacteriaceae adhesion mechanisms as well as future designs of engineered curli fiber based adhesives.

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.

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.

Molecular dynamics simulations of lysozyme-lipid systems: probing the early steps of protein aggregation

Using the molecular dynamics simulation, the role of lipids in the lysozyme transition into the aggregation-competent conformation has been clarified. Analysis of the changes of lysozyme secondary structure upon its interactions with the model bilayer membranes composed of phosphatidylcholine and its mixtures with phosphatidylglycerol (10, 40, and 80 mol%) within the time interval of 100 ns showed that lipid-bound protein is characterized by the increased content of β-structures. Along with this, the formation of protein-lipid complexes was accompanied by the increase in the gyration radius and the decrease in RMSD of polypeptide chain. The results obtained were interpreted in terms of the partial unfolding of lysozyme molecule on the lipid matrix, with the magnitude of this effect being increased with increasing the fraction of anionic lipids. Based on the results of molecular dynamics simulation, a hypothetical model of the nucleation of lysozyme amyloid fibrils in a membrane environment was suggested.

Molecular Dynamics Simulations of Amyloid β-Peptide (1-42): Tetramer Formation and Membrane Interactions

The aggregation cascade and peptide-membrane interactions of the amyloid β-peptide (Aβ) have been implicated as toxic events in the development and progression of Alzheimer's disease. Aβ42 forms oligomers and ultimately plaques, and it has been hypothesized that these oligomeric species are the main toxic species contributing to neuronal cell death. To better understand oligomerization events and subsequent oligomer-membrane interactions of Aβ42, we performed atomistic molecular-dynamics (MD) simulations to characterize both interpeptide interactions and perturbation of model membranes by the peptides. MD simulations were utilized to first show the formation of a tetramer unit by four separate Aβ42 peptides. Aβ42 tetramers adopted an oblate ellipsoid shape and showed a significant increase in β-strand formation in the final tetramer unit relative to the monomers, indicative of on-pathway events for fibril formation. The Aβ42 tetramer unit that formed in the initial simulations was used in subsequent MD simulations in the presence of a pure POPC or cholesterol-rich raft model membrane. Tetramer-membrane simulations resulted in elongation of the tetramer in the presence of both model membranes, with tetramer-raft interactions giving rise to the rearrangement of key hydrophobic regions in the tetramer and the formation of a more rod-like structure indicative of a fibril-seeding aggregate. Membrane perturbation by the tetramer was manifested in the form of more ordered, rigid membranes, with the pure POPC being affected to a greater extent than the raft membrane. These results provide critical atomistic insight into the aggregation pathway of Aβ42 and a putative toxic mechanism in the pathogenesis of Alzheimer's disease.

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.

Unmodified and pyroglutamylated amyloid β peptides form hypertoxic hetero-oligomers of unique secondary structure

Amyloid β (Aβ) peptide plays a major role in Alzheimer's disease (AD) and occurs in multiple forms, including pyroglutamylated Aβ (AβpE). Identification and characterization of the most cytotoxic Aβ species is necessary for advancement in AD diagnostics and therapeutics. While in brain tissue multiple Aβ species act in combination, structure/toxicity studies and immunotherapy trials have been focused on individual forms of Aβ. As a result, the molecular composition and the structural features of "toxic Aβ oligomers" have remained unresolved. Here, we have used a novel approach, hydration from gas phase coupled with isotope-edited Fourier transform infrared (FTIR) spectroscopy, to identify the prefibrillar assemblies formed by Aβ and AβpE and to resolve the structures of both peptides in combination. The peptides form unusual β-sheet oligomers stabilized by intramolecular H-bonding as opposed to intermolecular H-bonding in the fibrils. Time-dependent morphological changes in peptide assemblies have been visualized by atomic force microscopy. Aβ/AβpE hetero-oligomers exert unsurpassed cytotoxic effect on PC12 cells as compared to oligomers of individual peptides or fibrils. These findings lead to a novel concept that Aβ/AβpE hetero-oligomers, not just Aβ or AβpE oligomers, constitute the main neurotoxic conformation. The hetero-oligomers thus present a new biomarker that may be targeted for development of more efficient diagnostic and immunotherapeutic strategies to combat AD.

Alzheimer's Protective Cross-Interaction between Wild-Type and A2T Variants Alters Aβ42 Dimer Structure

Whole genome sequencing has recently revealed the protective effect of a single A2T mutation in heterozygous carriers against Alzheimer's disease (AD) and age-related cognitive decline. The impact of the protective cross-interaction between the wild-type (WT) and A2T variants on the dimer structure is therefore of high interest, as the Aβ dimers are the smallest known neurotoxic species. Toward this goal, extensive atomistic replica exchange molecular dynamics simulations of the solvated WT homo- and A2T hetero- Aβ_1-42 dimers have been performed, resulting into a total of 51 μs of sampling for each system. Weakening of a set of transient, intrachain contacts formed between the central and C-terminal hydrophobic residues is observed in the heterodimeric system. The majority of the heterodimers with reduced interaction between central and C-terminal regions lack any significant secondary structure and display a weak interchain interface. Interestingly, the A2T N-terminus, particularly residue F4, is frequently engaged in tertiary and quaternary interactions with central and C-terminal hydrophobic residues in those distinct structures, leading to hydrophobic burial. This atypical involvement of the N-terminus within A2T heterodimer revealed in our simulations implies possible interference on Aβ₄₂ aggregation and toxic oligomer formation, which is consistent with experiments. In conclusion, the present study provides detailed structural insights onto A2T Aβ₄₂ heterodimer, which might provide molecular insights onto the AD protective effect of the A2T mutation in the heterozygous state.

Self-assembly of the full-length amyloid Aβ42 protein in dimers

The self-assembly of amyloid (Aβ) proteins into nano-aggregates is a hallmark of Alzheimer's disease (AD) development, yet the mechanism of how disordered monomers assemble into aggregates remains elusive. Here, we applied long-time molecular dynamics simulations to fully characterize the assembly of Aβ42 monomers into dimers. Monomers undergo conformational changes during their interaction, but the resulting dimer structures do not resemble those found in fibril structures. To identify natural conformations of dimers among a set of simulated ones, validation approaches were developed and applied, and a subset of dimer conformations were characterized. These dimers do not contain long β-strands that are usually found in fibrils. The dimers are stabilized primarily by interactions within the central hydrophobic regions and the C-terminal regions, with a contribution from local hydrogen bonding. The dimers are dynamic, as evidenced by the existence of a set of conformations and by the quantitative analyses of the dimer dissociation process.

Polymorphism of fibrillar structures depending on the size of assembled Aβ17-42 peptides

The size of assembled Aβ_17-42 peptides can determine polymorphism during oligomerization and fibrillization, but the mechanism of this effect is unknown. Starting from separate random monomers, various fibrillar oligomers with distinct structural characteristics were identified using discontinuous molecular dynamics simulations based on a coarse-grained protein model. From the structures observed in the simulations, two characteristic oligomer sizes emerged, trimer and paranuclei, which generated distinct structural patterns during fibrillization. A majority of the simulations for trimers and tetramers formed non-fibrillar oligomers, which primarily progress to off-pathway oligomers. Pentamers and hexamers were significantly converted into U-shape fibrillar structures, meaning that these oligomers, called paranuclei, might be potent on-pathway intermediates in fibril formation. Fibrillar oligomers larger than hexamers generated substantial polymorphism in which hybrid structures were readily formed and homogeneous fibrillar structures appeared infrequently.

Isotope-edited FTIR reveals distinct aggregation and structural behaviors of unmodified and pyroglutamylated amyloid β peptides

Amyloid β peptide (Aβ) is causatively associated with Alzheimer's disease (AD), and N-terminally truncated and pyroglutamylated Aβ peptides (AβpE) exert hypertoxic effect by an unknown mechanism. Recent evidence has identified the prefibrillar oligomers of Aβ, not the fibrils, as the prevalent cytotoxic species. Structural characterization of Aβ and AβpE oligomers is therefore important for better understanding of their toxic effect. Here we have used isotope-edited Fourier transform infrared (FTIR) spectroscopy to identify the conformational changes in Aβ(1-42) and AβpE(3-42) upon aggregation, individually and in 1 : 1 molar combination. During the first two hours of exposure to aqueous buffer, the peptides undergo transition from mostly α-helical to mostly β-sheet structure. Data on peptides (13)C,(15)N-labeled at K(16)L(17)V(18) or V(36)G(37)G(38)V(39) allowed construction of structural models for the monomer and early oligomers. The peptide monomer comprises a β-hairpin that involves residues upstream of the K(16)L(17)V(18) sequence and an N-terminal α-helix. The oligomers form by non-H-bonding interactions between the β-strands of neighboring β-hairpins, in lateral or staggered manner, with the strands running parallel or antiparallel. Relative α-helical and β-sheet propensities of Aβ(1-42) and AβpE(3-42) depend on the ionic strength of the buffer, emphasizing the importance of ionic interactions in Aβ peptide structure and aggregation. It is inferred that N-terminal modification of AβpE(3-42) affects the helix stability and thereby modulates β-sheet oligomer formation. The data thus provide new insight into the molecular mechanism of Aβ oligomerization by emphasizing the role of the N-terminal transient α-helical structure and by identifying structural constraints for molecular organization of the oligomers.