Effect of the Tottori familial disease mutation (D7N) on the monomers and dimers of Aβ40 and Aβ42

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.

Effect of Terahertz Waves on the Structure of the Aβ42 Monomer, Dimer, and Protofibril: Insights from Molecular Dynamics Simulations

Amyloid-β (Aβ) and its assemblies play important roles in the pathogenesis of Alzheimer's disease (AD). Recent studies conducted by experimental and computational researchers have extensively explored the structure, assembly, and influence of biomolecules and cell membranes on Aβ. However, the impact of terahertz waves on the structures of Aβ monomers and aggregates remains largely unexplored. In this study, we systematically investigate the molecular mechanisms by which terahertz waves affect the structure of the Aβ42 monomer, dimer, and tetramer through all-atom molecular dynamics (MD) simulations. Our findings indicate that terahertz waves at a specific frequency (42.55 THz) can enhance intramolecular and intermolecular interactions in the Aβ42 monomer and dimer, respectively, by resonating with the symmetric stretching mode of the -COO- groups and the symmetric bending/stretching mode of -CH3 groups. Consequently, the β-structure content of the Aβ42 monomer is greatly increased, and the binding energy between the monomers in the Aβ42 dimer is significantly enhanced. Additionally, our observations suggest that terahertz waves can mildly stabilize the structure of tetrameric protofibrils by enhancing the interactions among peripheral peptides. Furthermore, we also investigated the effect of the frequency of terahertz waves on the structure of Aβ42. The present study contributes to a better understanding of the impact of external fields on the biobehavior of Aβ42 peptides and may shed some light on the potential risks associated with electromagnetic field radiation.

N-terminal Domain of Amyloid-β Impacts Fibrillation and Neurotoxicity

Alzheimer's disease is characterized by the presence of distinct amyloid-β peptide (Aβ) assemblies with diverse sizes, shapes, and toxicity. However, the primary determinants of Aβ aggregation and neurotoxicity remain unknown. Here, the N-terminal amino acid residues of Aβ42 that distinguished between humans and rats were substituted. The effects of these modifications on the ability of Aβ to aggregate and its neurotoxicity were investigated using biochemical, biophysical, and cellular techniques. The Aβ-derived diffusible ligand, protofibrils, and fibrils formed by the N-terminal mutational peptides, including Aβ42(R5G), Aβ42(Y10F), and rat Aβ42, were indistinguishable by conventional techniques such as size-exclusion chromatography, negative-staining transmission electron microscopy and silver staining, whereas the amyloid fibrillation detected by thioflavin T assay was greatly inhibited in vitro. Using circular dichroism spectroscopy, we discovered that both Aβ42 and Aβ42(Y10F) generated protofibrils and fibrils with a high proportion of parallel β-sheet structures. Furthermore, protofibrils formed by other mutant Aβ peptides and N-terminally shortened peptides were incapable of inducing neuronal death, with the exception of Aβ42 and Aβ42(Y10F). Our findings indicate that the N-terminus of Aβ is important for its fibrillation and neurotoxicity.

Atomistic Insights into A315E Mutation-Enhanced Pathogenicity of TDP-43 Core Fibrils

The aggregation of TAR DNA-binding protein of 43 kDa (TDP-43) into fibrillary deposits is implicated in amyotrophic lateral sclerosis (ALS), and some hereditary mutations localized in the low complexity domain (LCD) facilitate the formation of pathogenic TDP-43 fibrils. A recent cryo-EM study reported the atomic-level structures of the A315E TDP-43 LCD (residues 288-319, TDP-43_288-319) core fibril in which the protofilaments have R-shaped structures and hypothesized that A315E U-shaped protofilaments can readily convert to R-shaped protofilaments compared to the wild-type (WT) ones. There are no atomic structures of WT protofilaments available yet. Herein, we performed extensive all-atom explicit-solvent molecular dynamics simulations on A315E and WT protofilaments starting from both the cryo-EM-determined R-shaped and our constructed U-shaped structures. Our simulations show that WT protofilaments also adopt the R-shaped structures but are less stable than their A315E counterparts. Except for R293-E315 salt bridges, N312-F316 hydrophobic interactions and F316-F316 π-π stacking interactions are also crucial for the stabilization of the neck region of the R-shaped A315E protofilaments. The loss of R293-E315 salt bridges and the weakened interactions of N312-F316 and F316-F316 result in the reduced stability of the R-shaped WT protofilaments. Simulations starting from U-shaped folds reveal that A315E protofilaments can spontaneously convert to the cryo-EM-derived R-shaped protofilaments, whereas WT protofilaments convert to R-shape-like structures with remodeled neck regions. The R-shape-like WT protofilaments could act as intermediate states slowing down the U-to-R transition. This study reveals that A315E mutation can not only enhance the structural stability of the R-shaped TDP-43_288-319 protofilaments but also promote the U-to-R transition, which provides atomistic insights into the A315E mutation-enhanced TDP-43 pathogenicity in ALS.

Effect of Cholesterol Molecules on Aβ1-42 Wild-Type and Mutants Trimers

Alzheimer’s disease displays aggregates of the amyloid-beta (Aβ) peptide in the brain, and there is increasing evidence that cholesterol may contribute to the pathogenesis of the disease. Though many experimental and theoretical studies have focused on the interactions of Aβ oligomers with membrane models containing cholesterol, an understanding of the effect of free cholesterol on small Aβ42 oligomers is not fully established. To address this question, we report on replica exchange with a solute tempering simulation of an Aβ42 trimer with cholesterol and compare it with a previous replica exchange molecular dynamics simulation. We show that the binding hot spots of cholesterol are rather complex, involving hydrophobic residues L17–F20 and L30–M35 with a non-negligible contribution of loop residues D22–K28 and N-terminus residues. We also examine the effects of cholesterol on the trimers of the disease-causing A21G and disease-protective A2T mutations by molecular dynamics simulations. We show that these two mutations moderately impact cholesterol-binding modes. In our REST2 simulations, we find that cholesterol is rarely inserted into aggregates but rather attached as dimers and trimers at the surface of Aβ42 oligomers. We propose that cholesterol acts as a glue to speed up the formation of larger aggregates; this provides a mechanistic link between cholesterol and Alzheimer’s disease.

Computer Simulations Aimed at Exploring Protein Aggregation and Dissociation

Protein aggregation can lead to well-defined structures that are functional, but is also the cause of the death of neuron cells in many neurodegenerative diseases. The complexity of the molecular events involved in the aggregation kinetics of amyloid proteins and the transient and heterogeneous characters of all oligomers prevent high-resolution structural experiments. As a result, computer simulations have been used to determine the atomic structures of amyloid proteins at different association stages as well as to understand fibril dissociation. In this chapter, we first review the current computer simulation methods used for aggregation with some atomistic and coarse-grained results aimed at better characterizing the early formed oligomers and amyloid fibril formation. Then we present the applications of non-equilibrium molecular dynamics simulations to comprehend the dissociation of protein assemblies.

Dynamics of Amyloid Formation from Simplified Representation to Atomistic Simulations

Amyloid fibril formation is an intrinsic property of short peptides, non-disease proteins, and proteins associated with neurodegenerative diseases. Aggregates of the Aβ and tau proteins, the α-synuclein protein, and the prion protein are observed in the brain of Alzheimer's, Parkinson's, and prion disease patients, respectively. Due to the transient short-range and long-range interactions of all species and their high aggregation propensities, the conformational ensemble of these devastating proteins, the exception being for the monomeric prion protein, remains elusive by standard structural biology methods in bulk solution and in lipid membranes. To overcome these limitations, an increasing number of simulations using different sampling methods and protein models have been performed. In this chapter, we first review our main contributions to the field of amyloid protein simulations aimed at understanding the early aggregation steps of short linear amyloid peptides, the conformational ensemble of the Aβ40/42 dimers in bulk solution, and the stability of Aβ aggregates in lipid membrane models. Then we focus on our studies on the interactions of amyloid peptides/inhibitors to prevent aggregation, and long amyloid sequences, including new results on a monomeric tau construct.

The cellular modifier MOAG-4/SERF drives amyloid formation through charge complementation

While aggregation-prone proteins are known to accelerate aging and cause age-related diseases, the cellular mechanisms that drive their cytotoxicity remain unresolved. The orthologous proteins MOAG-4, SERF1A, and SERF2 have recently been identified as cellular modifiers of such proteotoxicity. Using a peptide array screening approach on human amyloidogenic proteins, we found that SERF2 interacted with protein segments enriched in negatively charged and hydrophobic, aromatic amino acids. The absence of such segments, or the neutralization of the positive charge in SERF2, prevented these interactions and abolished the amyloid-promoting activity of SERF2. In protein aggregation models in the nematode worm Caenorhabditis elegans, protein aggregation and toxicity were suppressed by mutating the endogenous locus of MOAG-4 to neutralize charge. Our data indicate that MOAG-4 and SERF2 drive protein aggregation and toxicity by interactions with negatively charged segments in aggregation-prone proteins. Such charge interactions might accelerate primary nucleation of amyloid by initiating structural changes and by decreasing colloidal stability. Our study points at charge interactions between cellular modifiers and amyloidogenic proteins as potential targets for interventions to reduce age-related protein toxicity.

Molecular mechanisms of resveratrol and EGCG in the inhibition of Aβ42 aggregation and disruption of Aβ42 protofibril: similarities and differences

The aggregation of amyloid-β protein (Aβ) into fibrillary deposits is implicated in Alzheimer's disease (AD), and inhibiting Aβ aggregation and clearing Aβ fibrils are considered as promising strategies to treat AD. It has been reported that resveratrol (RSV) and epigallocatechin-3-gallate (EGCG), two of the most extensively studied natural polyphenols, are able to inhibit Aβ fibrillization and remodel the preformed fibrillary aggregates into amorphous, non-toxic species. However, the mechanisms by which RSV inhibits Aβ₄₂ aggregation and disrupts Aβ₄₂ protofibril, as well as the inhibitory/disruptive mechanistic similarities and differences between RSV and EGCG, remain mostly elusive. Herein, we performed extensive all-atom molecular dynamics (MD) simulations on Aβ₄₂ dimers (the early aggregation state of Aβ₄₂) and protofibrils (the intermediate of Aβ₄₂ fibril formation and elongation) in the absence/presence of RSV or EGCG molecules. Our simulations show that both RSV and EGCG can bind with Aβ₄₂ monomers and inhibit the dimerization of Aβ₄₂. The binding of RSV with Aβ₄₂ peptide is mostly viaπ-π stacking interactions, while the binding of EGCG with Aβ₄₂ is mainly through hydrophobic, π-π stacking, and hydrogen-bonding interactions. Moreover, both RSV and EGCG disrupt the β-sheet structure and K28-A42 salt bridges, leading to a disruption of Aβ₄₂ protofibril structure. RSV mainly binds with residues whose side-chains point inwards from the surface of the protofibril, while EGCG mostly binds with residues whose side-chains point outwards from the surface of the protofibril. Furthermore, RSV interacts with Aβ₄₂ protofibrils mostly viaπ-π stacking interactions, while EGCG interacts with Aβ₄₂ protofibrils mainly via hydrogen-bonding and hydrophobic interactions. For comparison, we also explore the effects of RSV/EGCG molecules on the aggregation inhibition and protofibril disruption of the Iowa mutant (D23N) Aβ. Our findings may pave the way for the design of more effective drug candidates as well as the utilization of cocktail therapy using RSV and EGCG for the treatment of AD.

Small-Molecule Targeted Aβ42 Aggregate Degradation: Negatively Charged Small Molecules Are More Promising than the Neutral Ones

Heavy evidence has confirmed that Aβ42 oligomers are the most neurotoxic aggregates and play a critical role in the occurrence and development of Alzheimer's disease by causing functional neuron death, cognitive damage, and dementia. Disordered Aβ42 oligomers are challenging therapeutic targets, and no drug is currently in clinical use that modifies the properties of their monomeric states. Here, a negatively charged molecule (ER), rather than the neutral TS1 one, is identified by a molecular dynamics simulation method to be more capable of binding and sequestering the intrinsically disordered amyloid-β peptide Aβ42 in its soluble pentameric state as well as its monomeric components. Results reveal that the ERs interact with Aβ and inhibit the primary nucleation pathways in its aggregation process in entropic expansion mechanism for both Aβ42 and Aβ40 oligomers but with opposite characteristics of hydrophobic surface area (HSA). The interaction between Aβ42 oligomer and either charged ER or neutral TS1/TS0 characterizes decreased HSA, and the decrease in ER-involved case is highly visible, consistent with the observations from in silico and in vitro studies. By contrast, the presence of these inhibitors causes the HSA of Aβ40 oligomer to change undetectably and there is even a bit of increase in the histidine isomerized Aβ40 oligomer. The HSA distinction between Aβ42 and Aβ40 oligomer is possibly derived from the different effects of M35-inhibitor interaction, which is analogous to the effect of M35 oxidation. In comparison with the neutral TS1/TS0 inhibitors, ER is more prone to bind the residues located in the central (β1) and C-terminal (β2) regions of Aβ42 peptide, two key nucleation regions for Aβ intramolecular folding, intermolecular aggregation, and assembly. Notably, ER can strongly bind the charged residues, such as K16, K28, D23, to greatly disturb the potential stabilizer (e.g., salt-bridge, etc.) in metastable Aβ42 oligomers and protofibrils. These results illustrate the strategy of overcoming Alzheimer's disease from inhibiting its early stage Aβ aggregation with two kinds of small molecules to alter their behavior for therapeutic purposes and strongly recommend paying more attention to the engineering and development of negatively charged inhibitors, the long-term underappreciated ones, targeting the early stage Aβ aggregates.

Impact of Mutations on the Conformational Transition from α-Helix to β-Sheet Structures in Arctic-Type Aβ40: Insights from Molecular Dynamics Simulations

The amyloid-β (Aβ) protein aggregation into toxic oligomers and fibrils has been recognized as a key player in the pathogenesis of Alzheimer's disease. Recent experiments reported that a double alanine mutation (L17A/F19A) in the central hydrophobic core (CHC) region of [G22]Aβ₄₀ (familial Arctic mutation) diminished the self-assembly propensity of [G22]Aβ₄₀. However, the molecular mechanism behind the decreased aggregation tendency of [A17/A19/G22]Aβ₄₀ is not well understood. Herein, we carried out molecular dynamics simulations to elucidate the structure and dynamics of [G22]Aβ₄₀ and [A17/A19/G22]Aβ₄₀. The results for the secondary structure analysis reveal a significantly increased amount of the helical content in the CHC and C-terminal region of [A17/A19/G22]Aβ₄₀ as compared to [G22]Aβ₄₀. The bending free-energy analysis of D23-K28 salt bridge suggests that the double alanine mutation in the CHC region of [G22]Aβ₄₀ has the potential to reduce the fibril formation rate by 0.57 times of [G22]Aβ₄₀. Unlike [G22]Aβ₄₀, [A17/A19/G22]Aβ₄₀ largely sampled helical conformation, as determined by the minimum energy conformations extracted from the free-energy landscape. The present study provided atomic level details into the experimentally observed diminished aggregation tendency of [A17/A19/G22]Aβ₄₀ as compared to [G22]Aβ₄₀.

Structures of the intrinsically disordered Aβ, tau and α-synuclein proteins in aqueous solution from computer simulations

Intrinsically disordered proteins (IDPs) play many biological roles in the human proteome ranging from vesicular transport, signal transduction to neurodegenerative diseases. The Aβ and tau proteins, and the α-synuclein protein, key players in Alzheimer's and Parkinson's diseases, respectively are fully disordered at the monomer level. The structural heterogeneity of the monomeric and oligomeric states and the high self-assembly propensity of these three IDPs have precluded experimental structural determination. Simulations have been used to determine the atomic structures of these IDPs. In this article, we review recent computer models to capture the equilibrium ensemble of Aβ, tau and α-synuclein proteins at different association steps in aqueous solution and present new results of the PEP-FOLD framework on α-synuclein monomer.

Molecular dynamics simulations of copper binding to N-terminus mutants of amyloid-β

We report results of molecular dynamic (MD) simulations on N-terminus mutants of the copper-bound, amyloid-β (Aβ) peptide. Eight structures of Aβ were modelled, including seven mutant peptides in addition to the unaltered wild-type (WT). Trajectories analysed for each individual system were all approximately 1.4 μs in length, yielding a total of over 11 μs in total. The impact of these mutations are marked and varied compared to the wild-type peptide, including effects on secondary structure, stability and conformational changes. Each system showed differing levels of stability with some showing consistent, compact conformations whereas others displayed more flexible structures. Contrasts between comparable mutations at similar sites, such as A2T/A2V and D7H/D7N, show the location as well as the type of mutation have effects on protein structure observed in Ramachandran plots. We also report notable changes in peptide structure at residues remote to the site of substitution showing these mutations influence the entirety of Aβ. Salt-bridge profiles show this most clearly: addition or removal of charged residues affecting all salt-bridge interactions present in WT, even those remote from the site of mutation. Effects on secondary structure differ between mutations, most notably a change in incidence of β-strand, which has been linked to enhanced aggregational properties for the peptide. GFN2-xTB semi-empirical calculations show clear differences in binding energies of the copper-centre for each system.Communicated by Ramaswamy H. Sarma.

Effect of Piedmont mutation (L34V) on the structure, dynamics, and aggregation of Alzheimer's Aβ40 peptide

The amyloid-β (Aβ) aggregation in the brain has been associated with the development of Alzheimer's disease (AD). The previous studies have reported that Piedmont mutation (L34V) increases the rate of Aβ₄₀ aggregation. However, the underlying molecular mechanism of the effect of L34V mutation on Aβ₄₀ structure, dynamics, and aggregation remains largely unclear. In the present study, molecular dynamics (MD) simulations were performed to elucidate the effect of L34V mutation on the structural changes and conformational dynamics of Aβ₄₀. The secondary structure analysis highlight that L34V mutation enhances Aβ₄₀ self-assembly due to the formation of aggregation-prone β-sheet structure at the C-terminus of Aβ₄₀ monomeric structure. The higher probability of Asp23-Lys28 salt bridge interaction in Aβ₄₀(L34V) leads to aggregation prone β-sheet conformations, which has the potential to increase the fibril formation rate. The free energy landscape (FEL) analysis depict a sampling of coil conformation in the free energy minima of Aβ₄₀, whereas the aggregation-prone β-sheet conformation was observed at the C-terminal region of Aβ₄₀(L34V) in the minimum energy conformations extracted from FEL of Aβ₄₀(L34V). MD simulations, in agreement with experiment, highlight that L34V mutation increases Aβ₄₀ aggregation as the sampling of the aggregation-prone β-sheet conformation substantially increased. Overall, MD simulations provided atomic level details into the increased fibril formation tendency upon L34V mutation and physical insights into the L34V-mediated conformational as well as structural changes in Aβ₄₀.

Combining molecular dynamics simulations and experimental analyses in protein misfolding

The fold of a protein determines its function and its misfolding can result in loss-of-function defects. In addition, for certain proteins their misfolding can lead to gain-of-function toxicities resulting in protein misfolding diseases such as Alzheimer's, Parkinson's, or the prion diseases. In all of these diseases one or more proteins misfold and aggregate into disease-specific assemblies, often in the form of fibrillar amyloid deposits. Most, if not all, protein misfolding diseases share a fundamental molecular mechanism that governs the misfolding and subsequent aggregation. A wide variety of experimental methods have contributed to our knowledge about misfolded protein aggregates, some of which are briefly described in this review. The misfolding mechanism itself is difficult to investigate, as the necessary timescale and resolution of the misfolding events often lie outside of the observable parameter space. Molecular dynamics simulations fill this gap by virtue of their intrinsic, molecular perspective and the step-by-step iterative process that forms the basis of the simulations. This review focuses on molecular dynamics simulations and how they combine with experimental analyses to provide detailed insights into protein misfolding and the ensuing diseases.

Folding dynamics of Aβ42 monomer at pH 4.0–7.5 with and without physiological salt conditions – does the β1 or β2 region fold first?

The seeding region of Aβ42 monomer is jointly affected by the solution acidity, ionic distribution of the salt, and charged residues.

Sequence specificity despite intrinsic disorder: How a disease-associated Val/Met polymorphism rearranges tertiary interactions in a long disordered protein

The role of electrostatic interactions and mutations that change charge states in intrinsically disordered proteins (IDPs) is well-established, but many disease-associated mutations in IDPs are charge-neutral. The Val66Met single nucleotide polymorphism (SNP) in precursor brain-derived neurotrophic factor (BDNF) is one of the earliest SNPs to be associated with neuropsychiatric disorders, and the underlying molecular mechanism is unknown. Here we report on over 250 μs of fully-atomistic, explicit solvent, temperature replica-exchange molecular dynamics (MD) simulations of the 91 residue BDNF prodomain, for both the V66 and M66 sequence. The simulations were able to correctly reproduce the location of both local and non-local secondary structure changes due to the Val66Met mutation, when compared with NMR spectroscopy. We find that the change in local structure is mediated via entropic and sequence specific effects. We developed a hierarchical sequence-based framework for analysis and conceptualization, which first identifies “blobs” of 4-15 residues representing local globular regions or linkers. We use this framework within a novel test for enrichment of higher-order (tertiary) structure in disordered proteins; the size and shape of each blob is extracted from MD simulation of the real protein (RP), and used to parameterize a self-avoiding heterogenous polymer (SAHP). The SAHP version of the BDNF prodomain suggested a protein segmented into three regions, with a central long, highly disordered polyampholyte linker separating two globular regions. This effective segmentation was also observed in full simulations of the RP, but the Val66Met substitution significantly increased interactions across the linker, as well as the number of participating residues. The Val66Met substitution replaces β-bridging between V66 and V94 (on either side of the linker) with specific side-chain interactions between M66 and M95. The protein backbone in the vicinity of M95 is then free to form β-bridges with residues 31-41 near the N-terminus, which condenses the protein. A significant role for Met/Met interactions is consistent with previously-observed non-local effects of the Val66Met SNP, as well as established interactions between the Met66 sequence and a Met-rich receptor that initiates neuronal growth cone retraction.

Intrinsically disordered proteins are proteins that have no well-defined structure in at least one functional form. Mutations in one amino acid may still affect their function significantly, especially in subtle ways with cumulative adverse effects on health. Here we report on molecular dynamics simulations of a protein that is critical for neuronal health throughout adulthood (brain-derived neurotrophic factor). We investigate the effects of a mutation carried by 30% of human population, which has been widely studied for its association with aging-related and stress-related disorders, reduced volume of the hippocampus, and variations in episodic memory. We identify a molecular mechanism in which the mutation may change the global conformations of the protein and its ability to bind to receptors.

Unusual Segregation of APP Mutations in Monogenic Alzheimer Disease

APP gene mutations causing Alzheimer disease (AD) segregate in an autosomal dominant pattern. We report on a 40-year-old woman with a severe cognitive decline starting at 36 years, while her affected relatives presented symptoms onset in the 6th decade. The proband carried an APP missense variant in homozygous state (NM_000484.4: c.2032G>A; NP_000475.1: p.Asp678Asn; rs63750064) and showed a more severe clinical picture than the other AD relatives, as regards the age of onset and the rate of disease progression. This mutation behaves as a semi-dominant trait. The very rare chance of studying APP mutations in the homozygous state demonstrates they are not always dominant and other segregation models are possible.

Amyloid-β(29-42) Dimeric Conformations in Membranes Rich in Omega-3 and Omega-6 Polyunsaturated Fatty Acids

The omega-3 and omega-6 polyunsaturated fatty acids are two important components of cell membranes in human brains. When incorporated into phospholipids, omega-3 slows the progression of Alzheimer's disease (AD), whereas omega-6 is linked to increased risk of AD. Little is known on the amyloid-β (Aβ) conformations in membranes rich in omega-3 and omega-6 phospholipids. Herein, the structural properties of the Aβ_29-42 dimer embedded in both fatty acid membranes were comparatively studied to a 1-palmitoyl-2-oleoyl- sn-glycero-3-phosphocholine (POPC) bilayer using all-atom molecular dynamics (MD) simulations. Starting from α-helix, both omega-6 and omega-3 membranes promote new orientations and conformations of the dimer, in agreement with the observed dependence of Aβ production upon addition of these two fatty acids. This conformational result is corroborated by atomistic MD simulations of the dimer of the 99 amino acid C-terminal fragment of amyloid precursor protein spanning the residues 15-55. Starting from β-sheet, omega-6 membrane promotes helical and disordered structures of Aβ_29-42 dimer, whereas omega-3 membrane preserves the β-sheet structures differing however from those observed in POPC. Remarkably, the mixture of the two fatty acids and POPC depicts another conformational ensemble of the Aβ_29-42 dimer. This finding demonstrates that variation in the abundance of the molecular phospholipids, which changes with age, modulates membrane-embedded Aβ oligomerization.

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.

Compound CID 9998128 Is a Potential Multitarget Drug for Alzheimer's Disease

We have probed small molecule compound CID 9998128 as a potential multitarget drug for the Alzheimer's disease (AD) using in silico and in vitro experiments. By all-atom simulation and molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method, we have demonstrated that this compound strongly binds to both amyloid β42 (Aβ₄₂) fibrils and β-secretase, and the van der Waals interaction dominates over the electrostatic interaction in binding affinity. A detailed analysis at the atomic level revealed that indazole in CID 99998128 structure made a major contribution to instability of all studied complexes. In vitro experiments have shown that CID 9998128 inhibits the Aβ₄₂ amyloid fibrillization and is capable to clear Aβ₄₂ fibrils. Moreover, the compound dose-dependently decreases β-site amyloid precursor protein cleaving enzyme (BACE-1) activity with EC₅₀ value in micromolar range. Thus, our study has revealed that CID 9998128 is a good candidate for AD treatment through preventing production of Aβ peptides and degrading their aggregates. For drug design, we predict that the chemical structure of potent AD multitarget inhibitors should not contain indazole.

Effect of pH on Aβ42 Monomer and Fibril-like Oligomers-Decoding in Silico of the Roles of pK Values of Charged Residues

Amyloid beta-peptide (Aβ) is the key to developing Alzheimer's disease. Experiments have confirmed that different acidity influences directly not only the structural morphology and population of Aβ oligomers, but also the toxicity. The atomic-level association between the pH, charged residues, and Aβ properties remains obscure. Herein, conformational changes of Aβ₄₂ monomer, fibril-like trimer, and pentamer in the medium pH range of 4.0-7.5 are studied. The results reveal that, as the pH changes from 7.5 to the isoelectric pH, His6, His13, and His14 are protonated in turn, successively approach the center of mass of folded Aβ monomer, trigger ionic interactions and changes of neighboring turns (Asp7-Ser8, His14-Lys16) and even a distant one (Leu34-Met35), as well as concomitant changes of secondary structure, and promote the conformation transition from unfolded to folded. This observation discloses that protonation can convert these charged residues from originally hydrophilic to "hydrophobic-like". For fibril-like oligomers, the pH susceptibility essentially stems from the pK values of charged residues in the context of the Aβ fibril, and in turn one can predict the dynamic behavior of these residues in the processes of dissociation or stabilization of a fibril by comparing the pK values of residues involved in salt bridges in the normal state with those in the current context. This idea is justified by two fibril models and appears to be applicable to other peptides and their fibril systems.

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.

A chemically engineered, stable oligomer mimic of amyloid β42 containing an oxime switch for fibril formation

A stable Aβ oligomer mimic that is transformed into fibrils by a chemical stimulus, i.e., an oxime exchange reaction, is disclosed.

Dihydrochalcone molecules destabilize Alzheimer's amyloid-β protofibrils through binding to the protofibril cavity

Dihydrochalcone molecules destabilize Aβ_17–42protofibrils by disrupting the N-terminal β1 region and the turn region through binding to the protofibril cavity.

APP/Aβ structural diversity and Alzheimer's disease pathogenesis

The amyloid cascade hypothesis of Alzheimer's disease (AD) proposes amyloid- β (Aβ) is a chief pathological element of dementia. AD therapies have targeted monomeric and oligomeric Aβ 1-40 and 1-42 peptides. However, alternative APP proteolytic processing produces a complex roster of Aβ species. In addition, Aβ peptides are subject to extensive posttranslational modification (PTM). We propose that amplified production of some APP/Aβ species, perhaps exacerbated by differential gene expression and reduced peptide degradation, creates a diverse spectrum of modified species which disrupt brain homeostasis and accelerate AD neurodegeneration. We surveyed the literature to catalog Aβ PTM including species with isoAsp at positions 7 and 23 which may phenocopy the Tottori and Iowa Aβ mutations that result in early onset AD. We speculate that accumulation of these alterations induce changes in secondary and tertiary structure of Aβ that favor increased toxicity, and seeding and propagation in sporadic AD. Additionally, amyloid-β peptides with a pyroglutamate modification at position 3 and oxidation of Met35 make up a substantial portion of sporadic AD amyloid deposits. The intrinsic physical properties of these species, including resistance to degradation, an enhanced aggregation rate, increased neurotoxicity, and association with behavioral deficits, suggest their emergence is linked to dementia. The generation of specific 3D-molecular conformations of Aβ impart unique biophysical properties and a capacity to seed the prion-like global transmission of amyloid through the brain. The accumulation of rogue Aβ ultimately contributes to the destruction of vascular walls, neurons and glial cells culminating in dementia. A systematic examination of Aβ PTM and the analysis of the toxicity that they induced may help create essential biomarkers to more precisely stage AD pathology, design countermeasures and gauge the impacts of interventions.

Bexarotene Does Not Clear Amyloid Beta Plaques but Delays Fibril Growth: Molecular Mechanisms

In 2012, it was reported that anticancer drug bexarotene reduced amyloid plaque and improved mental functioning in a small sample of mice engineered to exhibit Alzheimer's like symptoms. It has been suggested that bexarotene stimulates expression of apolipoprotein E (ApoE) leading to intracellular clearance of amyloid beta (Aβ). However, the effect of bexarotene on clearance of plaques has not been seen in some mouse models. Two interesting questions include whether bexarotene can destroy Aβ fibrils via direct interaction with them and how this compound impacts the lag phase in the fibril growth process. By the Thioflavin T fluorescence assay and atomic force microscopy, we have shown that bexarotene prolongs the lag phase, but it does not degrade Aβ fibrils. The impotence of bexarotene in destroying fibrils means that this compound is weakly bound to Aβ. On the other hand, the weak binding would prevent bexarotene from prolonging the lag phase. Thus, our two main in vitro observations seem to contradict each other. In order to settle this problem at the atomic level, we have performed all-atom molecular dynamics simulations in explicit water. We have demonstrated that bexarotene is not capable to reduce amyloid deposits due to weak binding to Aβ fibrils. However, it delays the self-assembly through reduction of the β-content of Aβ monomers at high enough ligand concentrations. Bexarotene is the first compound which displays such an unusual behavior. We have also shown that bexarotene has a low binding propensity to Aβ monomer and dimer.

EGCG inhibits the oligomerization of amyloid beta (16-22) hexamer: Theoretical studies

An extensive replica exchange molecular dynamics (REMD) simulation was performed to investigate the progress patterns of the inhibition of (-)-epigallocatechin-3-gallate (EGCG) on the Aβ_16-22 hexamer. Structural variations of the oligomers without and with EGCG were monitored and analyzed in detail. It has been found that EGCG prevents the formation of Aβ oligomer through two different ways by either accelerating the Aβ oligomerization or reducing the β-content of the hexamer. It also decreases the potential "highly toxic" conformations of Aβ oligomer, which is related to the conformations having high order β-sheet sizes. Both electrostatic and van der Waals interaction energies are found to be involved to the binding process. Computed results using quantum chemical methods show that the π-π stacking is a critical factor of the interaction between EGCG and the peptides. As a result, the binding free energy of the EGCG to the Aβ peptides is slightly larger than that of the curcumin.

High-Resolution Structures of the Amyloid-β 1-42 Dimers from the Comparison of Four Atomistic Force Fields

The dimer of the amyloid-β peptide Aβ of 42 residues is the smallest toxic species in Alzheimer's disease, but its equilibrium structures are unknown. Here we determined the equilibrium ensembles generated by the four atomistic OPLS-AA, CHARMM22*, AMBER99sb-ildn, and AMBERsb14 force fields with the TIP3P water model. On the basis of 144 μs replica exchange molecular dynamics simulations (with 750 ns per replica), we find that the four force fields lead to random coil ensembles with calculated cross-collision sections, hydrodynamics properties, and small-angle X-ray scattering profiles independent of the force field. There are, however, marked differences in secondary structure, with the AMBERsb14 and CHARMM22* ensembles overestimating the CD-derived helix content, and the OPLS-AA and AMBER99sb-ildn secondary structure contents in agreement with CD data. Also the intramolecular beta-hairpin content spanning residues 17-21 and 30-36 varies between 1.5% and 13%. Overall, there are significant differences in tertiary and quaternary conformations among all force fields, and the key finding, irrespective of the force field, is that the dimer is stabilized by nonspecific interactions, explaining therefore its possible transient binding to multiple cellular partners and, in part, its toxicity.

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.

Conformational Ensembles of the Wild-Type and S8C Aβ1-42 Dimers

We characterized the dimer of the amyloid-β wild-type (WT) peptide, Aβ, of 42 residues and its disulfide-bond-locked double mutant (S8C) by replica exchange molecular dynamics simulations. Aβ dimers are known to be the smallest toxic species in Alzheimer's disease, and the S8C mutant has been shown experimentally to form an exclusive homogeneous and neurotoxic dimer. Our 50 μs all-atom simulations reveal similar secondary structures and collision cross-sections but very different intramolecular and intermolecular conformations upon double S8C mutation. Both dimers are very dynamic with hundreds of free-energy minima that differ from the U-shape and S-shape conformations of the peptides in the fibrils. The only common structural feature, shared by both species with a probability of 4% in WT and 12% in S8C-S8C, is a three-stranded β-sheet spanning the 17-23, 29-36, and 39-41 residues, which does not exist in the Aβ40 WT dimers.

Replica exchange molecular dynamics study of the amyloid beta (11–40) trimer penetrating a membrane

The transmembrane Aβ_11–40 trimer is investigated for the first time using REMD and FEP.

Replica exchange molecular dynamics study of the truncated amyloid beta (11–40) trimer in solution

The structure of the 3Aβ_11–40 oligomer is determined for the first time using T-REMD simulations.

Impact of Cu(II) Binding on Structures and Dynamics of Aβ42 Monomer and Dimer: Molecular Dynamics Study

The classical force field, which is compatible with the Amber force field 99SB, has been obtained for the interaction of Cu(II) with monomer and dimers of amyloid-β peptides using the coordination where Cu(II) is bound to His6, His13 (or His14), and Asp1 with distorted planar geometry. The newly developed force field and molecular dynamics simulation were employed to study the impact of Cu(II) binding on structures and dynamics of Aβ₄₂ monomer and dimers. It was shown that in the presence of Cu(II) the β content of monomer is reduced substantially compared with the wild-type Aβ₄₂ suggesting that, in accord with experiments, metal ions facilitate formation of amorphous aggregates rather than amyloid fibrils with cross-β structures. In addition, one possible mechanism for amorphous assembly is that the Asp23-Lys28 salt bridge, which plays a crucial role in β sheet formation, becomes more flexible upon copper ion binding to the Aβ N-terminus. The simulation of dimers was conducted with the Cu(II)/Aβ stoichiometric ratios of 1:1 and 1:2. For the 1:1 ratio Cu(II) delays the Aβ dimerization process as observed in a number of experiments. The mechanism underlying this phenomenon is associated with slow formation of interchain salt bridges in dimer as well as with decreased hydrophobicity of monomer upon Cu-binding.

Alzheimer's Amyloid-β Sequesters Caspase-3 in Vitro via Its C-Terminal Tail

Amyloid-β (Aβ), the main constituent in senile plaques found in the brain of patients with Alzheimer's disease (AD), is considered as a causative factor in AD pathogenesis. The clinical examination of the brains of patients with AD has demonstrated that caspase-3 colocalizes with senile plaques. Cellular studies have shown that Aβ can induce neuronal apoptosis via caspase-3 activation. Here, we performed biochemical and in silico studies to investigate possible direct effect of Aβ on caspase-3 to understand the molecular mechanism of the interaction between Aβ and caspase-3. We found that Aβ conformers can specifically and directly sequester caspase-3 activity in which freshly prepared Aβ42 is the most potent. The inhibition is noncompetitive, and the C-terminal region of Aβ plays an important role in sequestration. The binding of Aβ to caspase-3 was examined by cross-linking and proteolysis and by docking and all-atom molecular dynamic simulations. Experimental and in silico results revealed that Aβ42 exhibits a higher binding affinity than Aβ40 and the hydrophobic C-terminal region plays a key role in the caspase-Aβ interaction. Overall, our study describes a novel mechanism demonstrating that Aβ sequesters caspase-3 activity via direct interaction and facilitates future therapeutic development in AD.

Kinetics of the Interactions between Copper and Amyloid-β with FAD Mutations and Phosphorylation at the N terminus

Mutations and post‐translational modifications of amyloid‐β (Aβ) peptide in its N terminus have been shown to increase fibril formation, yet the molecular mechanism is not clear. Here we investigated the kinetics of the interactions of copper with two Aβ peptides containing Familial Alzheimer's disease (FAD) mutations (English (H6R) and Tottori (D7N)), as well as with Aβ peptide phosphorylated at serine 8 (pS8). All three peptides bind to copper with a similar rate as the wild‐type (wt). The dissociation rates follow the order pS8>H6R>wt>D7N; the interconversion between the two coordinating species occurs 50 % faster for H6R and pS8, whereas D7N had only a negligible effect. Interestingly, the rate of ternary complex (copper‐bridged heterodimer) formation for the modified peptides was significantly faster than that for wt, thus leading us to propose that FAD and sporadic AD might share a kinetic origin for the enhanced oligomerisation of Aβ.

Fullerenol C60(OH)16 prevents amyloid fibrillization of Aβ40-in vitro and in silico approach

The generation of Aβ amyloid aggregates in the form of senile plaques in the brain is one of the pathological hallmarks of Alzheimer's disease (AD). There is no cure for AD and one of the recent treatment strategies is focused on the inhibition of amyloid fibrillization of Aβ peptide. Fullerene C60 has been proposed as a candidate for destroying Aβ aggregates but it is not soluble in water and its toxicity to cells remains largely ambiguous. To overcome these drawbacks, we synthesized and studied the effect of water-soluble fullerenol C60(OH)16 (fullerene C60 carrying 16 hydroxyl groups) on the amyloid fibrillization of Aβ40 peptide in vitro. Using a Thioflavin T fluorescent assay and atomic force microscopy it was found that C60(OH)16 effectively reduces the formation of amyloid fibrils. The IC50 value is in the low range (μg ml(-1)) suggesting that fullerenol interferes with Aβ40 aggregation at stoichiometric concentrations. The in silico calculations supported the experimental data. It was revealed that fullerenol tightly binds to monomer Aβ40 and polar, negatively charged amino acids play a key role. Electrostatic interactions dominantly contribute to the binding propensity via interaction of the oxygen atoms from the COO(-) groups of side chains of polar, negatively charged amino acids with the OH groups of fullerenol. This stabilizes contact with either the D23 or K28 of the salt bridge. Due to the lack of a well-defined binding pocket fullerenol is also inclined to locate near the central hydrophobic region of Aβ40 and can bind to the hydrophobic C-terminal of the peptide. Upon fullerenol binding the salt bridge becomes flexible, inhibiting Aβ aggregation. In order to assess the toxicity of fullerenol, we found that exposure of neuroblastoma SH-SY5Y cells to fullerenol caused no significant changes in viability after 24 h of treatment. These results suggest that fullerenol C60(OH)16 represents a promising candidate as a therapeutic for Alzheimer's disease.