Role of Cholesterol on Binding of Amyloid Fibrils to Lipid Bilayers

Pore Formation by Amyloid-like Peptides: Effects of the Nonpolar-Polar Sequence Pattern

One of the mechanisms accounting for the toxicity of amyloid peptides in diseases like Alzheimer's and Parkinson's is the formation of pores on the plasma membrane of neurons. Here, we perform unbiased all-atom simulations of the full membrane damaging pathway, which includes adsorption, aggregation, and perforation of the lipid bilayer accounting for pore-like structures. Simulations are performed using four peptides made with the same amino acids. Differences in the nonpolar-polar sequence pattern of these peptides prompt them to adsorb into the membrane with the extended conformations oriented either parallel [peptide labeled F1, Ac-(FKFE)2-NH2], perpendicular (F4, Ac-FFFFKKEE-NH2), or with an intermediate orientation (F2, Ac-FFKKFFEE-NH2, and F3, Ac-FFFKFEKE-NH2) in regard to the membrane surface. At the water-lipid interface, only F1 fully self-assembles into β-sheets, and F2 peptides partially fold into an α-helical structure. The β-sheets of F1 emerge as electrostatic interactions attract neighboring peptides to intermediate distances where nonpolar side chains can interact within the dry core of the bilayer. This complex interplay between electrostatic and nonpolar interactions is not observed for the other peptides. Although β-sheets of F1 peptides are mostly parallel to the membrane, some of their edges penetrate deep inside the bilayer, dragging water molecules with them. This precedes pore formation, which starts with the flow of two water layers through the membrane that expand into a stable cylindrical pore delimited by polar faces of β-sheets spanning both leaflets of the bilayer.

Divalent cations promote huntingtin fibril formation on endoplasmic reticulum derived and model membranes

Huntington's Disease (HD) is caused by an abnormal expansion of the polyglutamine (polyQ) domain within the first exon of the huntingtin protein (htt). This expansion promotes disease-related htt aggregation into amyloid fibrils and the formation of proteinaceous inclusion bodies within neurons. Fibril formation is a complex heterogenous process involving an array of aggregate species such as oligomers, protofibrils, and fibrils. In HD, structural abnormalities of membranes of several organelles develop. In particular, the accumulation of htt fibrils near the endoplasmic reticulum (ER) impinges upon the membrane, resulting in ER damage, altered dynamics, and leakage of Ca2+. Here, the aggregation of htt at a bilayer interface assembled from ER-derived liposomes was investigated, and fibril formation directly on these membranes was enhanced. Based on these observations, simplified model systems were used to investigate mechanisms associated with htt aggregation on ER membranes. As the ER-derived liposome fractions contained residual Ca2+, the role of divalent cations was also investigated. In the absence of lipids, divalent cations had minimal impact on htt structure and aggregation. However, the presence of Ca2+ or Mg2+ played a key role in promoting fibril formation on lipid membranes despite reduced htt insertion into and association with lipid interfaces, suggesting that the ability of divalent cations to promote fibril formation on membranes is mediated by induced changes to the lipid membrane physicochemical properties. With enhanced concentrations of intracellular calcium being a hallmark of HD, the ability of divalent cations to influence htt aggregation at lipid membranes may play a role in aggregation events that lead to organelle abnormalities associated with disease.

Insights into the mechanism of peptide fibril growth on gold surface

Understanding the formation of β-fibrils over the gold surface is of paramount interest in nano-bio-medicinal Chemistry. The intricate mechanism of self-assembly of neurofibrillogenic peptides and their growth over the gold surface remains elusive, as experiments are limited in unveiling the microscopic dynamic details, in particular, at the early stage of the peptide aggregation. In this work, we carried out equilibrium molecular dynamics and enhanced sampling simulations to elucidate the underlying mechanism of the growth of an amyloid-forming sequence of tau fragments over the gold surface. Our results disclose that the collective intermolecular interactions between the peptide chains and peptides with the gold surface facilitate the peptide adsorption, followed by integration, finally leading to the fibril formation.

Recent Computational Advances Regarding Amyloid-β and Tau Membrane Interactions in Alzheimer's Disease

The interactions of amyloid proteins with membranes have been subject to many experimental and computational studies, as these interactions contribute in part to neurodegenerative diseases. In this review, we report on recent simulations that have focused on the adsorption and insertion modes of amyloid-β and tau proteins in membranes. The atomistic-resolution characterization of the conformational changes of these amyloid proteins upon lipid cell membrane and free lipid interactions is of interest to rationally design drugs targeting transient oligomers in Alzheimer's disease.

Integrative modeling of diverse protein-peptide systems using CABS-dock

The CABS model can be applied to a wide range of protein-protein and protein-peptide molecular modeling tasks, such as simulating folding pathways, predicting structures, docking, and analyzing the structural dynamics of molecular complexes. In this work, we use the CABS-dock tool in two diverse modeling tasks: 1) predicting the structures of amyloid protofilaments and 2) identifying cleavage sites in the peptide substrates of proteolytic enzymes. In the first case, simulations of the simultaneous docking of amyloidogenic peptides indicated that the CABS model can accurately predict the structures of amyloid protofilaments which have an in-register parallel architecture. Scoring based on a combination of symmetry criteria and estimated interaction energy values for bound monomers enables the identification of protofilament models that closely match their experimental structures for 5 out of 6 analyzed systems. For the second task, it has been shown that CABS-dock coarse-grained docking simulations can be used to identify the positions of cleavage sites in the peptide substrates of proteolytic enzymes. The cleavage site position was correctly identified for 12 out of 15 analyzed peptides. When combined with sequence-based methods, these docking simulations may lead to an efficient way of predicting cleavage sites in degraded proteins. The method also provides the atomic structures of enzyme-substrate complexes, which can give insights into enzyme-substrate interactions that are crucial for the design of new potent inhibitors.

Interaction of the tau fibrils with the neuronal membrane

Tau proteins are recently gaining a lot of interest due to their active role in causing a range of tauopathies. Molecular mechanisms underlying the tau interaction with the neuronal membrane are hitherto unknown and difficult to characterize using experimental methods. Using the cryo-EM structure of the tau-fibrils we have used atomistic molecular dynamics simulation to model the tau fibril and neuronal membrane interaction using explicit solvation. The dynamics and structural characteristics of the tau fibril with the neuronal membrane are compared to the tau fibril in the aqueous phase to corroborate the effect of the neuronal membrane in the tau structure. Tau fibrils have been modelled using CHARMM-36m force field and the six component neuronal membrane composition is taken from the earlier simulation results. The timescale conceivable in our molecular dynamics simulations is of the order of microseconds which captures the onset of the interaction of the tau fibrils with the neuronal membrane. This interaction is found to impact the tau pathogenesis that finally causes neuronal toxicity. Our study initiates the understanding of tau conformational ensemble in the presence of neuronal membrane and sheds the light on the significant tau-membrane interactions.

Lipid oxidation controls peptide self-assembly near membranes through a surface attraction mechanism

The self-assembly of peptides into supramolecular structures has been linked to neurodegenerative diseases but has also been observed in functional roles. Peptides are physiologically exposed to crowded environments of biomacromolecules, and particularly cellular membrane lipids. Previous research has shown that membranes can both accelerate and inhibit peptide self-assembly. Here, we studied the impact of membrane models that mimic cellular oxidative stress and compared this to mammalian and bacterial membranes. Using molecular dynamics simulations and experiments, we propose a model that explains how changes in peptide-membrane binding, electrostatics, and peptide secondary structure stabilization determine the nature of peptide self-assembly. We explored the influence of zwitterionic (POPC), anionic (POPG) and oxidized (PazePC) phospholipids, as well as cholesterol, and mixtures thereof, on the self-assembly kinetics of the amyloid β (1-40) peptide (Aβ40), linked to Alzheimer's disease, and the amyloid-forming antimicrobial peptide uperin 3.5 (U3.5). We show that the presence of an oxidized lipid had similar effects on peptide self-assembly as the bacterial mimetic membrane. While Aβ40 fibril formation was accelerated, U3.5 aggregation was inhibited by the same lipids at the same peptide-to-lipid ratio. We attribute these findings and peptide-specific effects to differences in peptide-membrane adsorption with U3.5 being more strongly bound to the membrane surface and stabilized in an α-helical conformation compared to Aβ40. Different peptide-to-lipid ratios resulted in different effects. We found that electrostatic interactions are a primary driving force for peptide-membrane interaction, enabling us to propose a model for predicting how cellular changes might impact peptide self-assembly in vivo.

Peptide-Membrane Binding: Effects of the Amino Acid Sequence

An understanding of how the amino acid sequence affects the interaction of peptides with lipid membranes remains mostly unknown. This type of knowledge is required to rationalize membrane-induced toxicity of amyloid peptides and to design peptides that can interact with lipid bilayers. Here, we perform a systematic study of how variations in the sequence of the amphipathic Ac-(FKFE)₂-NH₂ peptide affect its interaction with zwitterionic lipid bilayers using extensive all-atom molecular dynamics simulations in explicit solvent. Our results show that peptides with a net positive charge bind more frequently to the lipid bilayer than neutral or negatively charged sequences. Moreover, neutral amphipathic peptides made with the same numbers of phenylalanine (F), lysine (K), and glutamic (E) amino acids at different positions in the sequence differ significantly in their frequency of binding to the membrane. We find that peptides bind with a higher frequency to the membrane if their positive lysine side chains are more exposed to the solvent, which occurs if they are located at the extremity (as opposed to the middle) of the sequence. Non-polar residues play an important role in accounting for the adsorption of peptides onto the membrane. In particular, peptides made with less hydrophobic non-polar residues (e.g., valine and alanine) are significantly less adsorbed to the membrane compared to peptides made with phenylalanine. We also find that sequences where phenylalanine residues are located at the extremities of the peptide have a higher tendency to be adsorbed.

A coarse-grained molecular dynamics investigation on spontaneous binding of Aβ1-40 fibrils with cholesterol-mixed DPPC bilayers

Alzheimer's disease is the most common form of dementia. Its aetiology is characterized by the misfolding and aggregation of amyloid-β (Aβ) peptides into β-sheet-rich Aβ oligomers/fibrils. Although multiple experimental studies have suggested that Aβ oligomers/fibrils interact with the cell membranes and perturb their structures and dynamics, the molecular mechanism of this interaction is still not fully understood. In the present work, we have performed a total of 120 μs-long simulations to investigate the interaction between trimeric or hexameric Aβ_1-40 fibrils with either a 100% DPPC bilayer, a 70% DPPC-30% cholesterol bilayer or a 50% DPPC-50% cholesterol bilayer. Our simulation data capture the spontaneous binding of the aqueous Aβ_1-40 fibrils with the membranes and show that the central hydrophobic amino acid cluster, the lysine residue adjacent to it and the C-terminal hydrophobic residues are all involved in the process. Moreover, our data show that while the Aβ_1-40 fibril does not bind to the 100% DPPC bilayer, its binding affinity for the membrane increases with the amount of cholesterol. Overall, our data suggest that two clusters of hydrophobic residues and one lysine help Aβ_1-40 fibrils establish stable interactions with a cholesterol-rich DPPC bilayer. These residues are likely to represent potential target regions for the design of inhibitors, thus opening new avenues in structure-based drug design against Aβ oligomer/fibril-membrane interaction.

Exploring the binding kinetics and behaviors of self-aggregated beta-amyloid oligomers to phase-separated lipid rafts with or without ganglioside-clusters

Lipid binding kinetics and energetics of self-aggregated and disordered beta-amyloid oligomers of various sizes, from solution to lipid raft surfaces, were investigated using MD simulations. Our systems include small (monomers to tetramers) and larger (octamers and dodecamers) oligomers. Our lipid rafts contain saturated and unsaturated phosphatidylcholine (PC), cholesterol, and with or without asymmetrically distributed monosialotetrahexosylganglioside (GM1). All rafts exhibited dynamic and structurally diversified domains including liquid-ordered (Lo), liquid-disordered (Ld), and interfacial Lod domains. For rafts without GM1, all oligomers bound to the Lod domain. For GM1-containing rafts, all small oligomers and most larger oligomers bound specifically to the GM1-clusters embedded in the Lo domain. Lipid-protein binding energies followed an order of GM1 >> unsaturated PC > saturated PC > cholesterol for all rafts. In addition, protein-induced membrane structural disruption increased progressively with the size of the oligomer for the annular lipids surrounding the membrane-bound protein in non-GM1-containing rafts. We propose that the tight binding of beta-amyloid oligomers to the GM1-clusters and the structural perturbation of lipids surrounding the membrane-bound proteins at the Lod domain are early molecular events of the beta-amyloid aggregation process on neuronal membrane surfaces that trigger the onset of Alzheimer's.

Atomic Insights into Amyloid-Induced Membrane Damage

Amphipathic peptides can cause biological membranes to leak either by dissolving their lipid content via a detergent-like mechanism or by forming pores on the membrane surface. These modes of membrane damage have been related to the toxicity of amyloid peptides and to the activity of antimicrobial peptides. Here, we perform the first all-atom simulations in which membrane-bound amphipathic peptides self-assemble into β-sheets that subsequently either form stable pores inside the bilayer or drag lipids out of the membrane surface. An analysis of these simulations shows that the acyl tail of lipids interact strongly with non-polar side chains of peptides deposited on the membrane. These strong interactions enable lipids to be dragged out of the bilayer by oligomeric structures accounting for detergent-like damage. They also disturb the orientation of lipid tails in the vicinity of peptides. These distortions are minimized around pore structures. We also show that membrane-bound β-sheets become twisted with one of their extremities partially penetrating the lipid bilayer. This allows peptides on opposite leaflets to interact and form a long transmembrane β-sheet, which initiates poration. In simulations, where peptides are deposited on a single leaflet, the twist in β-sheets allows them to penetrate the membrane and form pores. In addition, our simulations show that fibril-like structures produce little damage to lipid membranes, as non-polar side chains in these structures are unavailable to interact with the acyl tail of lipids.

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.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

An Overview of Several Inhibitors for Alzheimer's Disease: Characterization and Failure

Amyloid beta (Aβ) oligomers are the most neurotoxic aggregates causing neuronal death and cognitive damage. A detailed elucidation of the aggregation pathways from oligomers to fibril formation is crucial to develop therapeutic strategies for Alzheimer's disease (AD). Although experimental techniques rely on the measure of time- and space-average properties, they face severe difficulties in the investigation of Aβ peptide aggregation due to their intrinsically disorder character. Computer simulation is a tool that allows tracing the molecular motion of molecules; hence it complements Aβ experiments, as it allows to explore the binding mechanism between metal ions and Aβ oligomers close to the cellular membrane at the atomic resolution. In this context, integrated studies of experiments and computer simulations can assist in mapping the complete pathways of aggregation and toxicity of Aβ peptides. Aβ oligomers are disordered proteins, and due to a rapid exploration of their intrinsic conformational space in real-time, they are challenging therapeutic targets. Therefore, no good drug candidate could have been identified for clinical use. Our previous investigations identified two small molecules, M30 (2-Octahydroisoquinolin-2(1H)-ylethanamine) and Gabapentin, capable of Aβ binding and inhibiting molecular aggregation, synaptotoxicity, intracellular calcium signaling, cellular toxicity and memory losses induced by Aβ. Thus, we recommend these molecules as novel candidates to assist anti-AD drug discovery in the near future. This review discusses the most recent research investigations about the Aβ dynamics in water, close contact with cell membranes, and several therapeutic strategies to remove plaque formation.

Binding Mechanisms of Amyloid-like Peptides to Lipid Bilayers and Effects of Divalent Cations

In several neurodegenerative diseases, cell toxicity can emerge from damage produced by amyloid aggregates to lipid membranes. The details accounting for this damage are poorly understood including how individual amyloid peptides interact with phospholipid membranes before aggregation. Here, we use all-atom molecular dynamics simulations to investigate the molecular mechanisms accounting for amyloid-membrane interactions and the role played by calcium ions in this interaction. Model peptides known to self-assemble into amyloid fibrils and bilayer made from zwitterionic and anionic lipids are used in this study. We find that both electrostatic and hydrophobic interactions contribute to peptide-bilayer binding. In particular, the attraction of peptides to lipid bilayers is dominated by electrostatic interactions between positive residues and negative phosphate moieties of lipid head groups. This attraction is stronger for anionic bilayers than for zwitterionic ones. Hydrophobicity drives the burial of nonpolar residues into the interior of the bilayer producing strong binding in our simulations. Moreover, we observe that the attraction of peptides to the bilayer is significantly reduced in the presence of calcium ions. This is due to the binding of calcium ions to negative phosphate moieties of lipid head groups, which leaves phospholipid bilayers with a net positive charge. Strong binding of the peptide to the membrane occurs less frequently in the presence of calcium ions and involves the formation of a "Ca²⁺ bridge".

Membrane Anchored Polymers Modulate Amyloid Fibrillation

The nucleating role of cellular membrane components, such as lipid moieties on amyloid beta (Aβ_1-40 ) fibrillation, has been reported in recent years. The influence of conjugates fabricated from lipid anchors (cholesterol, diacylglycerol) and hydrophilic polymers on Aβ_1-40 fibrillation is reported here, aiming to understand the impact of polymers cloud point temperature (T_cp ) and its hydrophobic tails on the amyloid fibrillation. Novel lipid-polymer conjugates, consisting of poly(oligo(ethylene glycol)_m acrylates) and hydrophobic groups (diacylglyceryl-, cholesteryl-, octyl-, decyl-, hexadecyl-) as anchors are synthesized using reversible addition-fragmentation chain transfer (RAFT) polymerization, allowing to tune the hydrophilic-hydrophobic profile of the conjugates by varying both, the degree of polymerization (n) and number of ethylene glycol units (m) in their side chain. The impact of these conjugates on Aβ_1-40 fibrillation is investigated via in vitro kinetic studies and transmission electron microscopy (TEM). Hydrophobic lipid-anchors are significantly delaying fibrillation (both lag- and half times), observing similar fibrillar structures via TEM when compared to native Aβ_1-40 . Other hydrophobic end groups are also delaying fibrillation of Aβ_1-40 , irrespective of their "n" and "m," whereas more hydrophilic polymers (both with longer ethylene glycol-sidechains, m = 3 for octyl, decyl and m = 5 for cholesterol) are only marginally inhibited fibrillation.

Cholesterol Molecules Alter the Energy Landscape of Small Aβ1-42 Oligomers

Small amyloid-β (Aβ) oligomers are believed to be key pathogenic species in Alzheimer's disease (AD). One suggested toxicity mechanism is the detergent model where oligomers remove lipid molecules from the bilayer. Senile plaques of AD patients also accumulate a 1:1 ratio of cholesterol/Aβ. What are the dominant structures of small Aβ42 oligomers with cholesterol molecules in aqueous solution? Here, we answer this question by performing atomistic replica exchange molecular dynamics simulations of Aβ42 dimers and trimers. Our simulations demonstrate that the interactions with cholesterol molecules change completely the energy landscape of small Aβ42 oligomers. This result shows that simulations in the bulk solution cannot recapitulate aggregation in the brain extracellular space.

Amyloidogenic Intrinsically Disordered Proteins: New Insights into Their Self-Assembly and Their Interaction with Membranes

Aβ, IAPP, α-synuclein, and prion proteins belong to the amyloidogenic intrinsically disordered proteins’ family; indeed, they lack well defined secondary and tertiary structures. It is generally acknowledged that they are involved, respectively, in Alzheimer’s, Type II Diabetes Mellitus, Parkinson’s, and Creutzfeldt–Jakob’s diseases. The molecular mechanism of toxicity is under intense debate, as many hypotheses concerning the involvement of the amyloid and the toxic oligomers have been proposed. However, the main role is represented by the interplay of protein and the cell membrane. Thus, the understanding of the interaction mechanism at the molecular level is crucial to shed light on the dynamics driving this phenomenon. There are plenty of factors influencing the interaction as mentioned above, however, the overall view is made trickier by the apparent irreproducibility and inconsistency of the data reported in the literature. Here, we contextualized this topic in a historical, and even more importantly, in a future perspective. We introduce two novel insights: the chemical equilibrium, always established in the aqueous phase between the free and the membrane phospholipids, as mediators of protein-transport into the core of the bilayer, and the symmetry-breaking of oligomeric aggregates forming an alternating array of partially ordered and disordered monomers.