Coarse-grained MD simulations reveal beta-amyloid fibrils of various sizes bind to interfacial liquid-ordered and liquid-disordered regions in phase separated lipid rafts with diverse membrane-bound conformational states

Multiscale Modeling of Macromolecular Interactions between Tau-Amylin Oligomers and Asymmetric Lipid Nanodomains That Link Alzheimer's and Diabetic Diseases

The molecular events of protein misfolding and self-aggregation of tau and amylin are associated with the progression of Alzheimer's and diabetes, respectively. Recent studies suggest that tau and amylin can form hetero-tau-amylin oligomers. Those hetero-oligomers are more neurotoxic than homo-tau oligomers. So far, the detailed interactions between the hetero-oligomers and the neuronal membrane are unknown. Using multiscale MD simulations, the lipid binding and protein folding behaviors of hetero-oligomers on asymmetric lipid nanodomains or raft membranes were examined. Our raft membranes contain phase-separated phosphatidylcholine (PC), cholesterol, and anionic phosphatidylserine (PS) or ganglioside (GM1) in one leaflet of the lipid bilayer. The hetero-oligomers bound more strongly to the PS and GM1 than other lipids via the hydrophobic and hydrophilic interactions, respectively, in the raft membranes. The hetero-tetramer disrupted the acyl chain orders of both PC and PS in the PS-containing raft membrane, but only the GM1 in the GM1-containing raft membrane as effectively as the homo-tau-tetramer. We discovered that the alpha-helical content in the heterodimer was greater than the sum of alpha-helical contents from isolated tau and amylin monomers on both raft membranes, indicative of a synergetic effect of tau-amylin interactions in surface-induced protein folding. Our results provide new molecular insights into understanding the cross-talk between Alzheimer's and diabetes.

Pragmatic Coarse-Graining of Proteins: Models and Applications

The molecular details involved in the folding, dynamics, organization, and interaction of proteins with other molecules are often difficult to assess by experimental techniques. Consequently, computational models play an ever-increasing role in the field. However, biological processes involving large-scale protein assemblies or long time scale dynamics are still computationally expensive to study in atomistic detail. For these applications, employing coarse-grained (CG) modeling approaches has become a key strategy. In this Review, we provide an overview of what we call pragmatic CG protein models, which are strategies combining, at least in part, a physics-based implementation and a top-down experimental approach to their parametrization. In particular, we focus on CG models in which most protein residues are represented by at least two beads, allowing these models to retain some degree of chemical specificity. A description of the main modern pragmatic protein CG models is provided, including a review of the most recent applications and an outlook on future perspectives in the field.

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.

Extension of the iSoLF implicit-solvent coarse-grained model for multicomponent lipid bilayers

iSoLF is a coarse-grained (CG) model for lipid molecules with the implicit-solvent approximation used in molecular dynamics (MD) simulations of biological membranes. Using the original iSoLF (iSoLFv1), MD simulations of lipid bilayers consisting of either POPC or DPPC and these bilayers, including membrane proteins, can be performed. Here, we improve the original model, explicitly treating the electrostatic interactions between different lipid molecules and adding CG particle types. As a result, the available lipid types increase to 30. To parameterize the potential functions of the new model, we performed all-atom MD simulations of each lipid at three different temperatures using the CHARMM36 force field and the modified TIP3P model. Then, we parameterized both the bonded and non-bonded interactions to fit the area per lipid and the membrane thickness of each lipid bilayer by using the multistate Boltzmann Inversion method. The final model reproduces the area per lipid and the membrane thickness of each lipid bilayer at the three temperatures. We also examined the applicability of the new model, iSoLFv2, to simulate the phase behaviors of mixtures of DOPC and DPPC at different concentrations. The simulation results with iSoLFv2 are consistent with those using Dry Martini and Martini 3, although iSoLFv2 requires much fewer computations. iSoLFv2 has been implemented in the GENESIS MD software and is publicly available.

Lipid domain boundary triggers membrane damage and protein folding of human islet amyloid polypeptide in the early pathogenesis of amyloid diseases

The misfolding and self-aggregation of human Islet Amyloid Polypeptide (hIAPP) are linked to the onset of type 2 diabetes (T2D). However, the mechanism of how the disordered hIAPP aggregates trigger membrane damage leading to the loss of Islet cells in T2D is unknown. Using coarse-grained (CG) and all-atom (AA) molecular dynamics simulations, we have investigated the membrane-disruption behaviors of hIAPP oligomers on the phase-separated lipid nanodomains that mimic the highly heterogeneous lipid raft structures of cell membranes. Our results revealed that hIAPP oligomers preferentially bind to the liquid-ordered and liquid-disordered domain boundary around two hydrophobic residues at L16 and I26, and lipid acyl chain order disruption and beta-sheet formation occur upon hIAPP binding to the membrane surface. We propose that the lipid order disruption and surface-induced beta-sheet formation on the lipid domain boundary represent the early molecular events of membrane damage associated with the early pathogenesis of T2D.

Amphipathic peptide-phospholipid nanofibers: Kinetics of fiber formation and molecular transfer between assemblies

Understanding the kinetics of nano-assembly formation is important to elucidate the biological processes involved and develop novel nanomaterials with biological functions. In the present study, we report the kinetic mechanisms of nanofiber formation from a mixture of phospholipids and the amphipathic peptide 18A[A11C], carrying cysteine substitution of the apolipoprotein A-I-derived peptide 18A at residue 11. 18A[A11C] with acetylated N-terminus and amidated C-terminus can associate with phosphatidylcholine to form fibrous aggregates at neutral pH and lipid-to-peptide molar ratio of ∼1, although the reaction pathways of self-assembly remain unclear. Here, the peptide was added to giant 1-palmitoyl-2-oleoyl phosphatidylcholine vesicles to monitor nanofiber formation under fluorescence microscopy. The peptide initially solubilized the lipid vesicles into particles smaller than the resolution of optical microscope, and fibrous aggregates appeared subsequently. Transmission electron microscopy and dynamic light scattering analyses revealed that the vesicle-solubilized particles were spherical or circular, measuring ∼10-20 nm in diameter. The rate of nanofiber formation of 18A with 1,2-dipalmitoyl phosphatidylcholine from the particles was proportional to the square of lipid-peptide concentration in the system, suggesting that the association of particles, accompanied by conformational changes, was the rate-limiting step. Moreover, molecules in the nanofibers could be transferred between aggregates faster than those in the lipid vesicles. These findings provide useful information for the development and control of nano-assembling structures using peptides and phospholipids.

Exploring Membrane Binding Targets of Disordered Human Tau Aggregates on Lipid Rafts Using Multiscale Molecular Dynamics Simulations

The self-aggregation of tau, a microtubule-binding protein, has been linked to the onset of Alzheimer's Disease. Recent studies indicate that the disordered tau aggregates, or oligomers, are more toxic than the ordered fibrils found in the intracellular neurofibrillary tangles of tau. At present, details of tau oligomer interactions with lipid rafts, a model of neuronal membranes, are not known. Using molecular dynamics simulations, the lipid-binding events, membrane-damage, and protein folding of tau oligomers on various lipid raft surfaces were investigated. Tau oligomers preferred to bind to the boundary domains (Lod) created by the coexisting liquid-ordered (Lo) and liquid-disordered (Ld) domains in the lipid rafts. Additionally, stronger binding of tau oligomers to the ganglioside (GM1) and phosphatidylserine (PS) domains, and subsequent protein-induced lipid chain order disruption and beta-sheet formation were detected. Our results suggest that GM1 and PS domains, located exclusively in the outer and inner leaflets, respectively, of the neuronal membranes, are specific membrane domain targets, whereas the Lod domains are non-specific targets, of tau oligomers binding to neurons. The molecular details of these specific and non-specific tau bindings to lipid rafts may provide new insights into understanding membrane-associated tauopathies leading to Alzheimer's Disease.

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.

Simulation study of domain formation in a model bacterial membrane

Recent experimental studies revealed that functional membrane microdomains (FMMs) are formed in prokaryotic cells which are structurally and functionally similar to the lipid rafts formed in eukaryotic cells. In this study, we employ coarse-grained molecular dynamics simulations to investigate the mechanism of domain formation and its physiochemical properties in a model methicillin-resistant staphylococcus aureus (MRSA) cell membrane. We find that domains are formed through lateral segregation of staphyloxanthin (STX), a carotenoid which shields the bacteria from the host's immune because of its antioxidant nature. Simulation results suggest that membrane integrity increases with the size of the domain, which is assessed by computing bond order parameter of the lipid tails, membrane expansion modulus and water permeability across the membrane. Various membrane domain proteins such as flotillin-like protein floA and penicillin binding protein (PBP2a) preferentially bind with the STX and accumulate in the membrane domain which is consistent with the recent experimental results.

Depth-Dependent Segmental Melting of the Sphingomyelin Alkyl Chain in Lipid Bilayers

The chain melting of lipid bilayers has often been investigated in detail using calorimetric methods, such as differential scanning calorimetry (DSC), and the resultant main transition temperature is regarded as one of the most important parameters in model membrane experiments. However, it is not always clear whether the hydrocarbon chains of lipids are gradually melting along the depth of the lipid bilayer or whether they all melt concurrently in a very narrow temperature range, as implied by DSC. In this study, we focused on stearoyl-d-sphingomyelin (SSM) as an example of raft-forming lipids. We synthesized deuterium-labeled SSMs at the 4', 10', and 16' positions, and their depth-dependent melting was measured using solid-state deuterium NMR by changing the temperature by 1.0 °C, and comparing with that observed from a saturated lipid, palmitoylstearoylphosphatidylcholine (PSPC). The results showed that SSM exhibited a characteristic depth-dependent melting, which was not observed for PSPC. The strong intermolecular hydrogen bonds between the sphingomyelin amide moiety probably caused the chain melting to start from the chain terminus through the middle part and end in the upper part. This depth-dependent melting implies that the small gel-like domains of SSM remain at temperatures slightly above the main transition temperature. These sphingomyelin features may be responsible for the biological properties of SM-based lipid rafts.

All-atom molecular dynamics simulations of the combined effects of different phospholipids and cholesterol content on electroporation

In this paper, we applied all-atom molecular dynamics (MD) simulations to study the effects of phospholipids and cholesterol content on bilayer membrane electroporation.

Interaction of membrane vesicles with the Pseudomonas functional amyloid protein FapC facilitates amyloid formation

Functional amyloids (FA) are proteins which are evolutionarily optimized to form highly stable fibrillar structures that strengthen the bacterial biofilm matrix. FA such as CsgA (E. coli) and FapC (Pseudomonas) are secreted to the bacterial surface where they integrate into growing fibril structures projecting from the outer membrane. FA are exposed to membrane surfaces in this process, but it remains unclear how membranes can interact with FA and potentially affect the self-assembly. Here we report the effect of different vesicles (DOPG, DMPG, DOPS, DOPC and DMPC) on the kinetics and structural endpoints of FapC fibrillation using various biophysical techniques. Particularly anionic lipids such as DMPG trigger FapC fibrillation, and the protein's second repeat sequence (R2) appears to be important for this interaction. Vesicles formed from phospholipids extracted from three different Pseudomonas strains (Δfap, ΔFapC and pfap) induce FapC fibrillation by accelerating nucleation. The general aggregation inhibitor epigallocatechin gallate (EGCG) inhibits FapC fibrillation by blocking interactions between FapC and vesicles and redirecting FapC monomers to oligomer structures. Our work indicates that biological membranes can contribute significantly to the fibrillation of functional amyloid.

Characterization of a Conformation-Restricted Amyloid β Peptide and Immunoreactivity of Its Antibody in Human AD brain

Characterization of amyloid β (Aβ) oligomers, the transition species present prior to the formation of Aβ fibrils and that have cytotoxicity, has become one of the major topics in the investigations of Alzheimer's disease (AD) pathogenesis. However, studying pathophysiological properties of Aβ oligomers is challenging due to the instability of these protein complexes in vitro. Here, we report that conformation-restricted Aβ42 with an intramolecular disulfide bond at positions 17 and 28 (SS-Aβ42) formed stable Aβ oligomers in vitro. Thioflavin T binding assays, nondenaturing gel electrophoresis, and morphological analyses revealed that SS-Aβ42 maintained oligomeric structure, whereas wild-type Aβ42 and the highly aggregative Aβ42 mutant with E22P substitution (E22P-Aβ42) formed Aβ fibrils. In agreement with these observations, SS-Aβ42 was more cytotoxic compared to the wild-type and E22P-Aβ42 in cell cultures. Furthermore, we developed a monoclonal antibody, designated TxCo-1, using the toxic conformation of SS-Aβ42 as immunogen. X-ray crystallography of the TxCo-1/SS-Aβ42 complex, enzyme immunoassay, and immunohistochemical studies confirmed the recognition site and specificity of TxCo-1 to SS-Aβ42. Immunohistochemistry with TxCo-1 antibody identified structures resembling senile plaques and vascular Aβ in brain samples of AD subjects. However, TxCo-1 immunoreactivity did not colocalize extensively with Aβ plaques identified with conventional Aβ antibodies. Together, these findings indicate that Aβ with a turn at positions 22 and 23, which is prone to form Aβ oligomers, could show strong cytotoxicity and accumulated in brains of AD subjects. The SS-Aβ42 and TxCo-1 antibody should facilitate understanding of the pathological role of Aβ with toxic conformation in AD.

Effect of packing density of lipid vesicles on the Aβ42 fibril polymorphism

The aggregation of amyloid-β 1-42 (Aβ42) on lipid membranes is closely related to the pathology of Alzheimer's disease (AD). Herein, we demonstrated the effect of the packing density of lipid vesicles on the Aβ42 fibrillation kinetics and fibril morphology. We used three distinct phosphatidylcholine (PC) lipids, containing different numbers of cis-double bonds in acyl chains, and therefore, a different packing density in the lipid vesicles. Our results showed that the fibrillation of Aβ42 was greatly enhanced and the formed fibrils became shorter as the number of double bonds in lipids increased. Due to the low-density characteristics of dioleoyl phosphatidylcholine (DOPC), Aβ42 monomers were able to interact with the hydrophobic acyl chain of lipids exposed to the aqueous phase, thereby inducing rapid fibrillation and short fibril morphologies. Furthermore, the effects of the anionic lipids dioleoyl phosphatidylserine (DOPS) and dioleoyl phosphatidylglycerol (DOPG), and mixed vesicles of DOPC/DOPS and DOPC/DOPG on Aβ42 fibrillations were investigated. The tight binding of Aβ42 to the lipid head groups via electrostatic interactions was able to suppress the modulation of Aβ42 fibrillations compared to accelerated fibrillations on loosely packed membranes. Our proposed mechanism regarding the influence of lipid packing density on Aβ42 fibrillations provides an advanced understanding of lipid-associated amyloid fibrillations.

The role of alpha-helix on the structure-targeting drug design of amyloidogenic proteins

The most accredited hypothesis links the toxicity of amyloid proteins to their harmful effects on membrane integrity through the formation of prefibrillar-transient oligomers able to disrupt cell membranes. However, damage mechanisms necessarily assume a first step in which the amyloidogenic protein transfers from the aqueous phase to the membrane hydrophobic core. This determinant step is still poorly understood. However, according to our lipid-chaperon hypothesis, free lipids in solution play a crucial role in facilitating this footfall. Free phospholipid concentration in the aqueous phase acts as a switch between ion channel-like pore and fibril formation, so that high free lipid concentration in solution promotes pore and repress fibril formation. Conversely, low free lipids in the solution favor fibril and repress pore formation. This behavior is due to the formation of stable lipid-protein complexes. Here, we hypothesize that the helix propensity is a fundamental requirement to fulfill the lipid-chaperon model. The alpha-helix region seems to be responsible for the binding with amphiphilic molecules fostering the proposed mechanism. Indeed, our results show the dependency of protein-lipid binding from the helical structure presence. When the helix content is substantially lower than the wild type, the contact probability decreases. Instead, if the helix is broadening, the contact probability increases. Our findings open a new perspective for in silico screening of secondary structure-targeting drugs of amyloidogenic proteins.

Exploring interactions between lipids and amyloid-forming proteins: A review on applying fluorescence and NMR techniques

A hallmark of Alzheimer's, Parkinson's, and other amyloid diseases is the assembly of amyloid proteins into amyloid aggregates or fibrils. In many cases, the formation and cytotoxicity of amyloid assemblies are associated with their interaction with cell membranes. Despite studied for many years, the characterization of the interaction is challenged for reasons on the multiple aggregation states of amyloid-forming proteins, transient and weak interactions in the complex system. Although several strategies such as computation biology, spectroscopy, and imaging methods have been performed, there is an urgent need to detail the molecular mechanism in different time scales and high resolutions. This review highlighted the recent applications of fluorescence, solution and solid-state NMR in exploring the interactions between amyloid protein and membranes attributing to their advantages of high sensitivity and atomic resolution.

Molecular dynamics simulations reveal the destabilization mechanism of Alzheimer's disease-related tau R3-R4 Protofilament by norepinephrine

Aggregation of Tau protein into neurofibrillary tangles is associated with the pathogenesis of Alzheimer's disease (AD) which has no cure yet. Clearing neurofibrillary tangles is one of major therapeutic strategies. Experimental studies reported that norepinephrine (NE) has the ability to disrupt Tau filament and cause Tau degradation. However, the underlying mechanism remains elusive. Herein, we performed molecular dynamic simulations to investigate the influence of NE on the C-shaped Tau R3-R4 protofilament. Our simulations show that NE compound destabilizes Tau protofilament by mostly disrupting β6/β8 and altering the β2-β3 and β6-β7 angles. NE binds mainly with aromatic residues Y310/P312/H374/F378 through ππ stacking and charged residues E338/E342/D348/D358/E372 via hydrogen-bonding interactions. Our results, together with the findings that exercise can markedly increase NE level, suggest that exercise might be a potent therapy against AD. This study reveals the disruptive mechanism of Tau protofilament by NE molecules, which may provide new clues for AD drug candidate design.

The Impact of Halogenated Phenylalanine Derivatives on NFGAIL Amyloid Formation

The hexapeptide hIAPP22-27 (NFGAIL) is known as a crucial amyloid core sequence of the human islet amyloid polypeptide (hIAPP) whose aggregates can be used to better understand the wild-type hIAPP's toxicity to β-cell death. In amyloid research, the role of hydrophobic and aromatic-aromatic interactions as potential driving forces during the aggregation process is controversially discussed not only in case of NFGAIL, but also for amyloidogenic peptides in general. We have used halogenation of the aromatic residue as a strategy to modulate hydrophobic and aromatic-aromatic interactions and prepared a library of NFGAIL variants containing fluorinated and iodinated phenylalanine analogues. We used thioflavin T staining, transmission electron microscopy (TEM) and small-angle X-ray scattering (SAXS) to study the impact of side-chain halogenation on NFGAIL amyloid formation kinetics. Our data revealed a synergy between aggregation behavior and hydrophobicity of the phenylalanine residue. This study introduces systematic fluorination as a toolbox to further investigate the nature of the amyloid self-assembly process.

Time-Dependent Lipid Dynamics, Organization and Peptide-Lipid Interaction in Phospholipid Bilayers with Incorporated β-Amyloid Oligomers

Nonfibrillar β-amyloid (Aβ) oligomers are considered as major neurotoxic species in the pathology of Alzheimer's disease. The presence of Aβ oligomers was shown to cause membrane disruptions in a broad range of model systems. However, the molecular basis of such a disruption process remains unknown. We previously demonstrated that membrane-incorporated 40-residue Aβ (Aβ₄₀) oligomers could form coaggregates with phospholipids. This process occurred more rapidly than the fibrillization of Aβ₄₀ and led to more severe membrane disruption. The present study probes the time-dependent changes in lipid dynamics, bilayer structures, and peptide-lipid interactions along the time course of the oligomer-induced membrane disruption, using solid-state NMR spectroscopy. Our results suggest the presence of certain intermediate states with phospholipid molecules entering the C-terminal hydrogen-bonding networks of the Aβ₄₀ oligomeric cores. This work provides insights on the molecular mechanisms of Aβ₄₀-oligomer-induced membrane disruption.

Data showing the lipid conformations and membrane binding behaviors of beta-amyloid fibrils in phase-separated cholesterol-enriched lipid domains with and without glycolipid and oxidized cholesterol from coarse-grained molecular dynamics simulations

The structural conformations of phospholipids and cholesterol in phase-separated lipid domains were determined by surface area, transverse density profile, and lipid acyl chain orientational parameter calculations. Binding kinetics and characterization of membrane-bound states of beta-amyloid fibrils of various sizes (dimer to pentamer), on those lipid domains, were determined using protein residue orientational parameter and fibril-residue-lipid minimum distance analysis methods. The energy of binding and characterization of annular lipid shells surrounding the surface-bound amyloid fibrils were also determined. The calculations described above support the article “Coarse-Grained MD simulations Reveal Diverse Membrane-Bound Conformational States of Beta-Amyloid Fibrils in the Liquid-ordered and Liquid-disordered Regions of Phase-Separated Lipid Rafts Containing Glycolipid, Cholesterol and Oxidized Cholesterol (Cheng et al., 2020 [1])”. The reported data is valuable for the future design and analysis of any protein fibrils binding to phase-separated lipid domains in model multi-component lipids membranes using either atomistic or coarse-grained molecular dynamics simulations. Additionally, this data can guide or validate future single-molecule experiments on fibril/membrane interactions in model or cell membranes.