Molecular Dynamics Study of the Interaction between the N-terminal of α-Synuclein and a Lipid Bilayer Mimicking Synaptic Vesicles

Folding of N-terminally acetylated α-synuclein upon interaction with lipid membranes

α-Synuclein (α-syn) is an abundant presynaptic neuronal protein whose aggregation is strongly associated with Parkinson's disease. It has been proposed that lipid membranes significantly affect α-syn's aggregation process. Extensive studies have been conducted to understand the interactions between α-syn and lipid membranes and have demonstrated that the N-terminus plays a critical role. However, the dynamics of the interactions and the conformational transitions of the N-terminus of α-syn at the atomistic scale details are still highly desired. In this study, we performed extensive enhanced sampling molecular dynamics simulations to quantify the folding and interactions of wild-type and N-terminally acetylated α-syn when interacting with lipid structures. We found that N-terminal acetylation significantly increases the helicity of the first few residues in solution or when interacting with lipid membranes. The observations in simulations showed that the binding of α-syn with lipid membranes mainly follows the induced-fit model, where the disordered α-syn binds with the lipid membrane through the electrostatic interactions and hydrophobic contacts with the packing defects; after stable insertion, N-terminal acetylation promotes the helical folding of the N-terminus to enhance the anchoring, thus increasing the binding affinity. We have shown the critical role of the first N-terminal residue methionine for recognition and anchoring to the negatively charged membrane. Although N-terminal acetylation neutralizes the positive charge of Met1 that may affect the electrostatic interactions of α-syn with membranes, the increase in helicity of the N-terminus should compensate for the binding affinity. This study provides detailed insight into the folding dynamics of α-syn's N-terminus with or without acetylation in solution and upon interaction with lipids, which clarifies how the N-terminal acetylation regulates the affinity of α-syn binding to lipid membranes. It also shows how packing defects and electrostatic effects coregulate the N-terminus of α-syn folding and its interaction with membranes.

Atomistic Insight into the Lipid Nanodomains of Synaptic Vesicles

Membrane curvature, once regarded as a passive consequence of membrane composition and cellular architecture, has been shown to actively modulate various properties of the cellular membrane. These changes could also lead to segregation of the constituents of the membrane, generating nanodomains with precise biological properties. Proteins often linked with neurodegeneration (e.g., tau, alpha-synuclein) exhibit an unintuitive affinity for synaptic vesicles in neurons, which are reported to lack distinct, ordered nanodomains based on their composition. In this study, all-atom molecular dynamics simulations are used to study a full-scale synaptic vesicle of realistic Gaussian curvature and its effect on the membrane dynamics and lipid nanodomain organization. Compelling indicators of nanodomain formation, from the perspective of composition, surface areas per lipid, order parameter, and domain lifetime, are identified in the vesicle membrane, which are absent in a flat bilayer of the same lipid composition. Therefore, our study supports the idea that curvature may induce phase separation in an otherwise fluid, disordered membrane.

Interaction between Permeation Enhancers and Lipid Bilayers

Permeation enhancers (PEs) are a class of molecules that interact with the epithelial membrane and transiently increase its transcellular permeability. Although there have been few clinical trials of PE coformulated drugs, the mechanism of action of PEs remains elusive. In this paper, the interaction between two archetypes of PEs [salcaprozate sodium (SNAC) and sodium caprate (C10)] and membranes is investigated with extensive all-atom molecular dynamics simulations. The simulations show that (1) the association between the neutral PEs and membranes is favored in free energy, (2) the propensity of neutral PE aggregation is larger in aqueous solution than in lipid bilayers, (3) the equilibrium distribution of neutral PEs in membranes is fast, e.g., accessible with unbiased MD simulations, and (4) the micelle of neutral PEs formed in aqueous solution does not rupture the membranes (e.g., not forming pores or breaking up the membrane) under simulation conditions. All results combined, this study indicates that PEs insert into the membranes in an equilibrium or near equilibrium process. This study lays the foundation for future investigations of how PEs impact the free energy of permeation for small molecules.

Impact of the Unstirred Water Layer on the Permeation of Small-Molecule Drugs

Over the last two decades, numerous molecular dynamics (MD) simulation-based investigations have attempted to predict the membrane permeability to small-molecule drugs as indicators of their bioavailability, a majority of which utilize the inhomogeneous solubility diffusion (ISD) model. However, MD-based membrane permeability is routinely 3-4 orders of magnitude larger than the values measured with the intestinal perfusion technique. There have been contentious discussions on the sources of the large discrepancies, and the two indisputable, potentially dominant ones are the fixed protonation state of the permeant and the neglect of the unstirred water layer (UWL). Employing six small-molecule drugs of different biopharmaceutical classification system classes, the current MD study relies on the ISD model but introduces the (de)protonation of the permeant by characterizing the permeation free energy of both neutral and charged states. In addition, the role of the UWL as a potential resistance against permeation is explored. The new MD protocol closely mimics the nature of small-molecule permeation and yields estimates that agree well with in vivo intestinal permeability.

Transmembrane β-Barrel Models of α-Synuclein Oligomers

The aggregation of α-synuclein is implicated in a number of neurodegenerative diseases, such as Parkinson's and Multiple System Atrophy, but the role of these aggregates in disease development is not clear. One possible mechanism of cytotoxicity is the disturbance or permeabilization of cell membranes by certain types of oligomers. However, no high-resolution structure of such membrane-embedded complexes has ever been determined. Here we construct and evaluate putative transmembrane β-barrels formed by this protein. Examination of the α-synuclein sequence reveals two regions that could form membrane-embedded β-hairpins: 64-92 (the NAC), and 35-56, which harbors many familial Parkinson's mutations. The stability of β-barrels formed by these hairpins is examined first in implicit membrane pores and then by multimicrosecond all-atom simulations. We find that a NAC region barrel remains stably inserted and hydrated for at least 10 μs. A 35-56 barrel remains stably inserted in the membrane but dehydrates and collapses if all His50 are neutral or if His50 is replaced by Q. If half of the His50 are doubly protonated, the barrel takes an oval shape but remains hydrated for at least 10 μs. Possible implications of these findings for α-synuclein pathology are discussed.

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.

Polyphenols-Based Nanosheets of Propolis Modulate Cytotoxic Amyloid Fibril Assembly of α-Synuclein

Natural compounds with anti-aggregation capacity are increasingly recognized as viable candidates against neurodegenerative diseases. Recently, the polyphenolic fraction of propolis (PFP), a complex bee product, has been shown to inhibit amyloid aggregation of a model protein especially in the nanosheet form. Here, we examine the aggregation-modulating effects of the PFP nanosheets on α-synuclein (α-syn), an intrinsically disordered protein involved in the pathogenesis of Parkinson's disease. Based on a range of biophysical data including intrinsic and extrinsic fluorescence, circular dichroism (CD) data, and nuclear magnetic resonance spectroscopy, we propose a model for the interaction of α-syn with PFP nanosheets, where the positively charged N-terminal and the middle non-amyloid component regions of α-syn act as the main binding sites with the negatively charged PFP nanosheets. The Thioflavin T (ThT) fluorescence, Congo red absorbance, and CD data reveal a prominent dose-dependent inhibitory effect of PFP nanosheets on α-syn amyloid aggregation, and the microscopy images and MTT assay data suggest that the PFP nanosheets redirect α-syn aggregation toward nontoxic off-pathway oligomers. When preformed α-syn amyloid fibrils are present, fluorescence images show co-localization of PFP nanosheets and ThT, further confirming the binding of PFP nanosheets with α-syn amyloid fibrils. Taken together, our results demonstrate the binding and anti-aggregation activity of PFP nanosheets in a disease-related protein system and propose them as potential nature-based tools for probing and targeting pathological protein aggregates in neurodegenerative diseases.

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.

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".