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

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.

Mechanism of Cationic Lipid Induced DNA Condensation: Lipid-DNA Coordination and Divalent Cation Charge Fluctuations

The condensation of nucleic acids by lipids is a widespread phenomenon in biology with crucial implications for drug delivery. However, the mechanisms of DNA assembly in lipid bilayers remain insufficiently understood due to challenges in measuring and assessing each component's contribution in the lipid-DNA-cation system. This study uses all-atom molecular dynamics simulations to investigate DNA condensation in cationic lipid bilayers. Our exhaustive exploration of the thermodynamic factors reveals unique roles for phospholipid head groups and cations. We observed that bridging cations between lipid and DNA drastically reduce charges, while mobile magnesium cations "ping-ponging" between double strands create charge fluctuations. While the first factor stabilizes the DNA-lipid complex, the latter creates attractive forces to induce the spontaneous condensation of DNAs. This novel mechanism not only sheds light on the current data regarding cationic lipid-induced DNA condensation but also provides potential design strategies for creating efficient gene delivery vectors for drug delivery.

Data-mining unveils structure-property-activity correlation of viral infectivity enhancing self-assembling peptides

Gene therapy via retroviral vectors holds great promise for treating a variety of serious diseases. It requires the use of additives to boost infectivity. Amyloid-like peptide nanofibers (PNFs) were shown to efficiently enhance retroviral gene transfer. However, the underlying mode of action of these peptides remains largely unknown. Data-mining is an efficient method to systematically study structure-function relationship and unveil patterns in a database. This data-mining study elucidates the multi-scale structure-property-activity relationship of transduction enhancing peptides for retroviral gene transfer. In contrast to previous reports, we find that not the amyloid fibrils themselves, but rather µm-sized β-sheet rich aggregates enhance infectivity. Specifically, microscopic aggregation of β-sheet rich amyloid structures with a hydrophobic surface pattern and positive surface charge are identified as key material properties. We validate the reliability of the amphiphilic sequence pattern and the general applicability of the key properties by rationally creating new active sequences and identifying short amyloidal peptides from various pathogenic and functional origin. Data-mining-even for small datasets-enables the development of new efficient retroviral transduction enhancers and provides important insights into the diverse bioactivity of the functional material class of amyloids.

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.

Influence of Cholesterol on the Membrane Binding and Conformation of α-Synuclein

The α-Synuclein (α-Syn) plays an important role in the pathology of Parkinson's disease (PD), and its oligomers and fibrils are toxic to the nervous system. As organisms age, the cholesterol content in biological membranes increases, which is a potential cause of PD. Cholesterol may affect the membrane binding of α-Syn and its abnormal aggregation, but the mechanism remains unclear. Here, we present our molecular dynamics simulation studies on the interaction between α-Syn and lipid membranes, with or without cholesterol. It is demonstrated that cholesterol provides additional hydrogen bond interaction with α-Syn; however, the coulomb interaction and hydrophobic interaction between α-Syn and lipid membranes could be weakened by cholesterol. In addition, cholesterol leads to the shrinking of lipid packing defects and the decrease of lipid fluidity, thereby shortening the membrane binding region of α-Syn. Under these multifaceted effects of cholesterol, membrane-bound α-Syn shows signs of forming a β-sheet structure, which may further induce the formation of abnormal α-Syn fibrils. These results provide important information for the understanding of membrane binding of α-Syn, and they are expected to promote the bridging between cholesterol and the pathological aggregation of α-Syn.

“Interplay of lipid-head group and packing defects in driving Amyloid-beta mediated myelin-like model membrane deformation”

Accumulating evidence suggests that amyloid plaque associated myelin lipid loss as a result of elevated amyloid burden might also contribute to Alzheimer's disease. The amyloid fibrils though closely associated with lipids under physiological conditions, however, the progression of membrane remodeling events leading to lipid-fibril assembly remains unknown. Here we first reconstitute the interaction of Aβ-40 with myelin-like model membrane and show that the binding of Aβ-40 induces extensive tubulation. To look into the mechanism of membrane tubulation we chose a set of membrane conditions varying in lipid packing density and net charge that allows us to dissect the contribution of lipid specificity of Aβ-40 binding, aggregation kinetics, and subsequent changes in membrane parameters such as fluidity, diffusion, and compressibility modulus. We show that the binding of Aβ-40 depends predominantly on the lipid packing defect densities and electrostatic interactions and results in rigidification of the myelin-like model membrane during the early phase of amyloid aggregation. Furthermore, elongation of Aβ-40 into higher oligomeric and fibrillar species leads to eventual fluidization of the model membrane followed by extensive lipid membrane tubulation observed in the late phase. Taken together, our results capture mechanistic insights into snapshots of temporal dynamics of Aβ-40 - myelin-like model membrane interaction and demonstrate how short timescale, local phenomena of binding, and fibril-mediated load generation results in the consequent association of lipids with growing amyloid fibrils.

Probing the interactions between amyloidogenic proteins and bio-membranes

Protein misfolding diseases (PMDs) in humans are characterized by the deposition of protein aggregates in tissues, including Alzheimer's disease, Parkinson's disease, type 2 diabetes, and amyotrophic lateral sclerosis. Misfolding and aggregation of amyloidogenic proteins play a central role in the onset and progression of PMDs, and these processes are regulated by multiple factors, especially the interaction between proteins and bio-membranes. Bio-membranes induce conformational changes in amyloidogenic proteins and affect their aggregation; on the other hand, the aggregates of amyloidogenic proteins may cause membrane damage or dysfunction leading to cytotoxicity. In this review, we summarize the factors that affect the binding of amyloidogenic proteins and membranes, the effects of bio-membranes on the aggregation of amyloidogenic proteins, mechanisms of membrane disruption by amyloidogenic aggregates, technical approaches for detecting these interactions, and finally therapeutic strategies targeting membrane damage caused by amyloidogenic proteins.

The Mechanism of Action of SAAP-148 Antimicrobial Peptide as Studied with NMR and Molecular Dynamics Simulations

Background: SAAP-148 is an antimicrobial peptide derived from LL-37. It exhibits excellent activity against drug-resistant bacteria and biofilms while resisting degradation in physiological conditions. Despite its optimal pharmacological properties, its mechanism of action at the molecular level has not been explored. Methods: The structural properties of SAAP-148 and its interaction with phospholipid membranes mimicking mammalian and bacterial cells were studied using liquid and solid-state NMR spectroscopy as well as molecular dynamics simulations. Results: SAAP-148 is partially structured in solution and stabilizes its helical conformation when interacting with DPC micelles. The orientation of the helix within the micelles was defined by paramagnetic relaxation enhancements and found similar to that obtained using solid-state NMR, where the tilt and pitch angles were determined based on 15N chemical shift in oriented models of bacterial membranes (POPE/POPG). Molecular dynamic simulations revealed that SAAP-148 approaches the bacterial membrane by forming salt bridges between lysine and arginine residues and lipid phosphate groups while interacting minimally with mammalian models containing POPC and cholesterol. Conclusions: SAAP-148 stabilizes its helical fold onto bacterial-like membranes, placing its helix axis almost perpendicular to the surface normal, thus probably acting by a carpet-like mechanism on the bacterial membrane rather than forming well-defined pores.

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.

The Increased Amyloidogenicity of Spike RBD and pH-Dependent Binding to ACE2 May Contribute to the Transmissibility and Pathogenic Properties of SARS-CoV-2 Omicron as Suggested by In Silico Study

SARS-CoV-2 is a rapidly evolving pathogen that has caused a global pandemic characterized by several consecutive waves. Based on epidemiological and NGS data, many different variants of SARS-CoV-2 were described and characterized since the original variant emerged in Wuhan in 2019. Notably, SARS-CoV-2 variants differ in transmissibility and pathogenicity in the human population, although the molecular basis for this difference is still debatable. A significant role is attributed to amino acid changes in the binding surface of the Spike protein to the ACE2 receptor, which may facilitate virus entry into the cell or contribute to immune evasion. We modeled in silico the interaction between Spike RBDs of Wuhan-Hu-1, Delta, and Omicron BA.1 variants and ACE2 at different pHs (pH 5 and pH 7) and showed that the strength of this interaction was higher for the Omicron BA.1 RBD compared to Wuhan-Hu-1 or Delta RBDs and that the effect was more profound at pH 5. This finding is strikingly related to the increased ability of Omicron variants to spread in the population. We also noted that during its spread in the population, SARS-CoV-2 evolved to a more charged, basic composition. We hypothesize that the more basic surface of the Omicron variant may facilitate its spread in the upper respiratory tract but not in the lower respiratory tract, where pH estimates are different. We calculated the amyloidogenic properties of Spike RBDs in different SARS-CoV-2 variants and found eight amyloidogenic regions in the Spike RBDs for each of the variants predicted by the FoldAmyloid program. Although all eight regions were almost identical in the Wuhan to Gamma variants, two of them were significantly longer in both Omicron variants, making the Omicron RBD more amyloidogenic. We discuss how the increased predicted amyloidogenicity of the Omicron variants RBDs may be important for protein stability, influence its interaction with ACE2 and contribute to immune evasion.

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.

Early-Stage β-Amyloid-Membrane Interactions Modulate Lipid Dynamics and Influence Structural Interfaces and Fibrillation

Molecular interactions between β-amyloid (Aβ) peptide and membranes contribute to the neuronal toxicity of Aβ and the pathology of Alzheimer's disease (AD). Neuronal plasma membranes serve as biologically relevant environments for the Aβ aggregation process as well as affect the structural polymorphisms of Aβ aggregates. However, the nature of these interactions is unknown. Here, we utilized solid-state NMR spectroscopy to explore the site-specific interactions between Aβ peptides and lipids in synaptic plasma membranes at the membrane-associated nucleation stage. The key results show that different segments in the hydrophobic sequence of Aβ initiate membrane binding and inter-strand assembling. We demonstrate early-stage Aβ-lipid interactions modulate lipid dynamics, leading to more rapid lipid headgroup motion and reduced lateral diffusive motion. These early events influence the structural polymorphisms of yielded membrane-associated Aβ fibrils with distinct C-terminal quaternary interface structure compared to fibrils grown in aqueous solutions. Based on our results, we propose a schematic mechanism by which Aβ-lipid interactions drive membrane-associated nucleation processes, providing molecular insights into the early events of fibrillation in biological environments.

Physical Mechanism of Amyloid-β Peptide Chain Aggregation on Fluidic Lipid Nanotubules

The question of how peptide chain aggregation is influenced by lipid membranes with varying shapes and structures is crucial for a detailed understanding of the neurotoxicity effect of the peptide chains. Not like the more usual spherical liposomes and planar lipid membranes, herein, we use lipid nanotubules as a model of important neuron synapse nanowire structures and devote particular attention to the effect of nanotubule fluidity on amyloid-β peptide (Aβ) chain aggregation. We apply single-molecule tracking (SMT) to elucidate how Aβ chains diffuse and aggregate on lipid nanotubules with different fluidities. The physical mechanism implies that fluidic lipid nanotubules facilitate the super-diffusion of two-dimensional (2D)-mobile precursor Aβ chains and promote their aggregation. This aggregation mechanism is retarded on less fluidic lipid nanotubules where the super-diffusion of 2D-mobile precursor Aβ chains is restricted by "frozen" lipids with less mobility. This work provides a mechanistic explanation for Aβ chain aggregation on fluidic lipid nanotubules.

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.