Stability of Aβ11-40 Trimers with Parallel and Antiparallel β-Sheet Organizations in a Membrane-Mimicking Environment by Replica Exchange Molecular Dynamics Simulation

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.

An S-Shaped Aβ42 Cross-β Hexamer Embedded into a Lipid Bilayer Reveals Membrane Disruption and Permeability

The interactions of amyloid oligomers with membranes are known to contribute to cellular toxicity. Numerous in vitro experimental studies reported on the insertion of oligomers of different sizes that can induce cell membrane disruption, extract lipids, and form ion-permeable transmembrane pores. The current repertoire of amyloid-beta (Aβ) membrane-inserted folds that was subject to high-resolution structure NMR spectroscopy and computer simulations is devoid of any cross-β fibrillar structure. In this study, we explored the dynamics of an S-shaped Aβ42 cross-β hexamer model inserted into a lipid bilayer membrane by two atomistic molecular dynamics simulations. The initial model is characterized by the hydrophobic residues at the central hydrophobic core (residues 17-21, CHC) and the C-terminus (residues 30-42) embedded into the membrane. We observed major structural secondary, tertiary, and quaternary rearrangements leading to two distinct species, hexamer and two trimers, accompanied by membrane disruption and water permeation. The simulations show that some configurations, but not the majority, have the CHC and C-terminus hydrophobic residues exposed to the solvent. Overall, our computational results offer new perspectives to understand the relationship between Aβ42 assemblies and membrane permeability.

Concentration-Dependent Effects of Cholesterol on the Dimerization of Amyloid-β Peptides in Lipid Bilayers

Membrane disruption mediated by the accumulation of amyloid-β (Aβ) on cell membranes is central to the pathogenesis of Alzheimer's disease (AD). Cholesterol, an important component of membranes, is well-recognized as a risk factor in AD. It can affect the aggregation and pore formation of Aβ on membranes whereas the specific effects are rather complex, particularly regarding the non-linear response to cholesterol concentrations. Yet, the mechanistic understanding of the role of cholesterol in Aβ-membrane interactions remains incomplete. Herein, we employed microsecond-scale molecular dynamics simulations to investigate the effects of cholesterol on Aβ dimerization in a lipid bilayer containing different molar ratios of cholesterol (0, 20, and 40 mol %). Cholesterol reduces the time required for the formation of stable dimers and exerts dual effects on Aβ-membrane interactions. First, cholesterol promotes the extraction of the C-terminal region from the membrane to water. Consequently, at the ratios of 0 and 20 mol %, peptides are anchored at the membrane-water interface, but they are repelled to water at a ratio of 40 mol % with high structural flexibility. Second, cholesterol weakens Aβ-membrane interactions, thereby enhancing inter-peptide interactions. The former is favorable for dimerization while the latter is not. The balance between two factors eventually leads to a non-monotonic effect on the degree of dimerization, whereby the number of inter-peptide contacts is the largest at a cholesterol ratio of 20 mol %. These results provide atomistic insights into the regulation mechanism of Aβ42 aggregation by cholesterol and help to understand the pathological link between cholesterol and AD.

501Y.V2 spike protein resists the neutralizing antibody in atomistic simulations

SARS-CoV-2 outbreaks worldwide caused COVID-19 pandemic, which is related to several million deaths. In particular, SARS-CoV-2 Spike (S) protein is a major biological target for COVID-19 vaccine design. Unfortunately, recent reports indicated that Spike (S) protein mutations can lead to antibody resistance. However, understanding the process is limited, especially at the atomic scale. The structural change of S protein and neutralizing antibody fragment (FAb) complexes was thus probed using molecular dynamics (MD) simulations. In particular, the backbone RMSD of the 501Y.V2 complex was significantly larger than that of the wild-type one implying a large structural change of the mutation system. Moreover, the mean of Rg, CCS, and SASA are almost the same when compared two complexes, but the distributions of these values are absolutely different. Furthermore, the free energy landscape of the complexes was significantly changed when the 501Y.V2 variant was induced. The binding pose between S protein and FAb was thus altered. The FAb-binding affinity to S protein was thus reduced due to revealing over steered-MD (SMD) simulations. The observation is in good agreement with the respective experiment that the 501Y.V2 SARS-CoV-2 variant can escape from neutralizing antibody (NAb).

Molecular insights into the primary nucleation of polymorphic amyloid β dimers in DOPC lipid bilayer membrane

Alzheimer's disease (AD) pathology is characterized by loss of memory cognitive and behavioral deterioration. One of the hallmarks of AD is amyloid β (Aβ) plaques in the brain that consists of Aβ oligomers and fibrils. It is accepted that oligomers, particularly dimers, are toxic species that are produced extracellularly and intracellularly in membranes. It is believed that the disruption of membranes by polymorphic Aβ oligomers is the key for the pathology of AD. This is a first study that investigate the effect of polymorphic “α‐helix/random coil” and “fibril‐like” Aβ dimers on 1,2‐dioleoyl‐sn‐glycero‐3‐phosphocholine (DOPC) membrane. It has been found that the DOPC membrane promotes Aβ_1–42 “fibril‐like” dimers and impedes Aβ_1–42 “α‐helix/random coil” dimers. The N‐termini domains within Aβ_1–42 dimers play a role in Aβ aggregation in membrane milieus. In addition, the aromatic π–π interactions (involving residues F19 and F20 in Aβ_1–42) are the driving forces for the hydrophobic interactions that initiate the primary nucleation of polymorphic Aβ_1–42 dimers within DOPC membrane. Finally, the DOPC bilayer membrane thickness is locally decreased, and it is disrupted by an embedded distinct Aβ_1–42 dimer, due to relatively large contacts between Aβ_1–42 monomers and the DOPC membrane. This study reveals insights into the molecular mechanisms by which polymorphic early‐stage Aβ_1–42 dimers have distinct impacts on DOPC membrane.

Computational Models for the Study of Protein Aggregation

Protein aggregation has been studied by many groups around the world for many years because it can be the cause of a number of neurodegenerative diseases that have no effective treatment. Obtaining the structure of related fibrils and toxic oligomers, as well as describing the pathways and main factors that govern the self-organization process, is of paramount importance, but it is also very difficult. To solve this problem, experimental and computational methods are often combined to get the most out of each method. The effectiveness of the computational approach largely depends on the construction of a reasonable molecular model. Here we discussed different versions of the four most popular all-atom force fields AMBER, CHARMM, GROMOS, and OPLS, which have been developed for folded and intrinsically disordered proteins, or both. Continuous and discrete coarse-grained models, which were mainly used to study the kinetics of aggregation, are also summarized.

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.

Effect of Lauric Acid on the Stability of Aβ42 Oligomers

While Alzheimer's disease is correlated with the presence of Aβ fibrils in patient brains, the more likely agents are their precursors, soluble oligomers that may form pores or otherwise distort cell membranes. Using all-atom molecular dynamics simulation, we study how the presence of fatty acids such as lauric acid changes the stability of pore-forming oligomers built from three-stranded Aβ₄₂ chains. Such a change would alter the distribution of amyloids in the fatty acid-rich brain environment and therefore could explain the lower polymorphism observed in Aβ fibrils derived from brains of patients with Alzheimer's disease. We find that lauric acid stabilizes both ring-like and barrel-shaped models, with the effect being stronger for barrel-like models than for ring-like oligomers.

The F19W mutation reduces the binding affinity of the transmembrane Aβ11-40 trimer to the membrane bilayer

Alzheimer's disease is linked to the aggregation of the amyloid-β protein (Aβ) of 40 or 42 amino acids. Lipid membranes are known to modulate the rate and mechanisms of the Aβ aggregation. Point mutations in Aβ can alter these rates and mechanisms. In particular, experiments show that F19 mutations influence the aggregation rate, but maintain the fibril structures. Here, we used molecular dynamics simulations to examine the effect of the F19W mutation in the 3Aβ11-40 trimer immersed in DPPC lipid bilayers submerged in aqueous solution. Substituting Phe by its closest (non-polar) aromatic amino acid Trp has a dramatic reduction in binding affinity to the phospholipid membrane (measured with respect to the solvated protein) compared to the wild type: the binding free energy of the protein-DPPC lipid bilayer increases by 40-50 kcal mol-1 over the wild-type. This is accompanied by conformational changes and loss of salt bridges, as well as a more complex free energy surface, all indicative of a more flexible and less stable mutated trimer. These results suggest that the impact of mutations can be assessed, at least partially, by evaluating the interaction of the mutated peptides with the lipid membranes.

Emergence of Barrel Motif in Amyloid-β Trimer: A Computational Study

Amyloid-β (Aβ) peptides form assemblies that are pathological hallmarks of Alzheimer's disease. Aβ oligomers are soluble, mobile, and toxic forms of the peptide that act in the extracellular space before assembling into protofibrils and fibrils. Therefore, oligomers play an important role in the mechanism of Alzheimer's disease. Since it is difficult to determine by experiment the atomic structures of oligomers, which accumulate fast and are polymorphic, computer simulation is a useful tool to investigate elusive oligomers' structures. In this work, we report extended all-atom molecular dynamics simulations, both canonical and replica exchange, of Aβ(1-42) trimer starting from two different initial conformations: (i) the pose produced by the best docking of a monomer aside of a dimer (simulation 1), representing oligomers freshly formed by assembling monomers, and (ii) a configuration extracted from an experimental mature fibril structure (simulation 2), representing settled oligomers in equilibrium with extended fibrils. We showed that in simulation 1, regions with small β-barrels are populated, indicating the chance of spontaneous formation of domains resembling channel-like structures. These structural domains are alternative to those more representative of mature fibrils (simulation 2), the latter showing a stable bundle of C-termini that is not sampled in simulation 1. Moreover, trimer of Aβ(1-42) can form internal pores that are large enough to be accessed by water molecules and Ca2+ ions.

High Speed AFM and NanoInfrared Spectroscopy Investigation of Aβ1-42 Peptide Variants and Their Interaction With POPC/SM/Chol/GM1 Model Membranes

Due to an aging population, neurodegenerative diseases such as Alzheimer's disease (AD) have become a major health issue. In the case of AD, Aβ_{1 - 42} peptides have been identified as one of the markers of the disease with the formation of senile plaques via their aggregation, and could play a role in memory impairment and other tragic syndromes associated with the disease. Many studies have shown that not only the morphology and structure of Aβ_{1 - 42} peptide assembly are playing an important role in the formation of amyloid plaques, but also the interactions between Aβ_{1 - 42} and the cellular membrane are crucial regarding the aggregation processes and toxicity of the amyloid peptides. Despite the increasing amount of information on AD associated amyloids and their toxicity, the molecular mechanisms involved still remain unclear and require in-depth investigation at the local scale to clearly decipher the role of the sequence of the amyloid peptides, of their secondary structures, of their oligomeric states, and of their interactions with lipid membranes. In this original study, through the use of Atomic Force Microscopy (AFM) related-techniques, high-speed AFM and nanoInfrared AFM, we tried to unravel at the nanoscale the link between aggregation state, structure and interaction with membranes in the amyloid/membrane interaction. Using three mutants of Aβ peptides, L34T, oG37C, and WT Aβ_{1 - 42} peptides, with differences in morphology, structure and assembly process, as well as model lipidic membranes whose composition and structure allow interactions with the peptides, our AFM study coupling high spatial and temporal resolution and nanoscale structure information clearly evidences a local correlation between the secondary structure of the peptides, their fibrillization kinetics and their interactions with model membranes. Membrane disruption is associated to small transient oligomeric entities in the early stages of aggregation that strongly interact with the membrane, and present an antiparallel β-sheet secondary structure. The strong effect on membrane integrity that exists when these oligomeric Aβ_{1 - 42} peptides interact with membranes of a particular composition could be a lead for therapeutic studies.

Tau R3-R4 Domain Dimer of the Wild Type and Phosphorylated Ser356 Sequences. I. In Solution by Atomistic Simulations

In Alzheimer's disease, neurofibrillary lesions correlate with cognitive deficits and consist of inclusions of tau protein with cross-β structure. A stable dimeric form of soluble tau has been evidenced in the cells, but its high-resolution structure is missing in solution. We know, however, that cryo-electron microscopy (c-EM) of full-length tau in the brain of an individual with AD displays a core of eight β-sheets with a C-shaped architecture spanning the R3-R4 repeat domain, while the rest of the protein is very flexible. To address the conformational ensemble of the dimer, we performed atomistic replica exchange molecular dynamics simulations on the tau R3-R4 domain starting from the c-EM configuration. We find that the wild type tau R3-R4 dimer explores elongated, U-shaped, V-shaped, and globular forms rather than the C-shape. Phosphorylation of Ser356, pSer356, is known to block the interaction between the tau protein and the amyloid-β42 peptide. Standard molecular dynamics simulations of this phosphorylated sequence for a total of 5 μs compared to its wild type counterpart show a modulation of the population of β-helices and accessible topologies and a decrease of intermediates near the fibril-like conformers.