Effect of lipid saturation on amyloid-beta peptide partitioning and aggregation in neuronal membranes: molecular dynamics simulations

Development of Self-Assembled Biomimetic Nanoscale Collagen-like Peptide-Based Scaffolds for Tissue Engineering: An In Silico and Laboratory Study

Development of biocomposite scaffolds has gained tremendous attention due to their potential for tissue regeneration. However, most scaffolds often contain animal-derived collagen that may elicit an immunological response, necessitating the development of new biomaterials. Herein, we developed a new collagen-like peptide,(Pro-Ala-His)10 (PAH)10, and explored its ability to be utilized as a functional biomaterial by incorporating it with a newly synthesized peptide-based self-assembled gel. The gel was prepared by conjugating a pectin derivative, galataric acid, with a pro-angiogenic peptide (LHYQDLLQLQY) and further functionalized with a cortistatin-derived peptide, (Phe-Trp-Lys-Thr)4 (FWKT)4, and the bio-ionic liquid choline acetate. The self-assembly of (PAH)10 and its interactions with the galactarate-peptide conjugates were examined using replica exchange molecular dynamics (REMD) simulations. Results revealed the formation of a multi-layered scaffold, with enhanced stability at higher temperatures. We then synthesized the scaffold and examined its physicochemical properties and its ability to integrate with aortic smooth muscle cells. The scaffold was further utilized as a bioink for bioprinting to form three-dimensional cell-scaffold matrices. Furthermore, the formation of actin filaments and elongated cell morphology was observed. These results indicate that the (PAH)10 hybrid scaffold provides a suitable environment for cell adhesion, proliferation and growth, making it a potentially valuable biomaterial for tissue engineering.

Aggregation and partitioning of amyloid peptide fragments in the presence of a lipid bilayer: A coarse grained molecular dynamics study

Amyloidogenicity and toxicity of amyloid peptides have been linked to the peptide aggregation and interactions with lipid bilayers. In this work we used the coarse grained MARTINI model to study the aggregation and partitioning of amyloid peptide fragments Aβ(1-28) and Aβ(25-35) in the presence of a dipalmitoylphosphatidylcholine bilayer. We explored the peptide aggregation starting from three initial spatial arrangements where free monomers were placed in solution outside the membrane, at the membrane-solution interface, or in the membrane. We found that Aβ(1-28) and Aβ(25-35) interact with the bilayer quite differently. The Aβ(1-28) fragments show strong peptide-peptide and peptide-lipid interactions leading to irreversible aggregation where the aggregates stay confined to their initial spatial location. The Aβ(25-35) fragments show weaker peptide-peptide and peptide-lipid interaction leading to reversible aggregation and accumulation at the membrane-solution interface irrespective of their initial spatial arrangement. Those findings can be explained in terms of the shape of the potential of mean force for the single-peptide translocation across the membrane.

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.

Destabilization of Aβ fibrils by omega-3 polyunsaturated fatty acids: a molecular dynamics study

The senile plaques of neurotoxic aggregates of Aβ protein, deposited extraneuronally, mark the pathological hallmark of Alzheimer's disease (AD). The natural compounds such as omega-3 (ω-3) polyunsaturated fatty acids (PUFAs), which can access blood-brain barrier, are believed to be potential disruptors of preformed Aβ fibrils to cure AD with unknown mechanism. Herein, we present the destabilization potential of three ω-3 PUFAs, viz. Eicosapentaenoic acid (EPA), Docosahexaenoic acid (HXA), and α-linolenic acid (LNL) by molecular dynamics simulation. After an initial testing of 300 ns, EPA and HXA have been considered further for extended production run time, 500 ns. The increased value of root mean square deviation (RMSD), radius of gyration, and solvent-accessible surface area (SASA), the reduced number of H-bonds and β-sheet content, and disruption of salt bridges and hydrophobic contacts establish the binding of these ligands to Aβ fibril leading to destabilization. The polar head was found to interact with positively charged lysine (K28) residue in the fibril. However, the hydrophobicity of the long aliphatic tail competes with the intrinsic hydrophobic interactions of Aβ fibril. This amphiphilic nature of EPA and HXA led to the breaking of inherent hydrophobic contacts and formation of new bonds between the tail of PUFA and hydrophobic residues of Aβ fibril, leading to the destabilization of fibril. The Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) results explain the binding of EPA and HXA to Aβ fibril by interacting with different residues. The destabilization potential of EPA and HXA establishes them as promising drug leads to cure AD, and encourages prospecting of other fatty acids for therapeutic intervention in AD.Communicated by Ramaswamy H. Sarma.

Design of short peptides and peptide amphiphiles as collagen mimics and an investigation of their interactions with collagen using molecular dynamics simulations and docking studies

Short peptide sequences and bolaamphiphiles derived from natural proteins are gaining importance due to their ability to form unique nanoscale architectures for a variety of biological applications. In this work, we have designed six short peptides (triplet or monomeric forms) and two peptide bolaamphiphiles that either incorporate the bioactive collagen motif (Gly-X-Y) or sequences where Gly, Pro, or hydroxyproline (Hyp) are replaced by Ala or His. For the bolaamphiphiles, a malate moiety was used as the aliphatic linker for connecting His with Hyp to create collagen mimics. Stability of the assemblies was assessed through molecular dynamics simulations and results indicated that (Pro-Ala-His)₃ and (Ala-His-Hyp)₃ formed the most stable structures, while the amphiphiles and the monomers showed some disintegration over the course of the 200 ns simulation, though most regained structural integrity and formed fibrillar structures, and micelles by the end of the simulation, likely due to the formation of more thermodynamically stable conformations. Multiple replica simulations (REMD) were also conducted where the sequences were simulated at different temperatures. Our results showed excellent convergence in most cases compared to constant temperature molecular dynamics simulation. Furthermore, molecular docking and MD simulations of the sequences bound to collagen triple helix structure revealed that several of the sequences had a high binding affinity and formed stable complexes, particularly (Pro-Ala-His)₃ and (Ala-His-Hyp)₃. Thus, we have designed new hybrid-peptide-based sequences which may be developed for potential applications as biomaterials for tissue engineering or drug delivery.

Two statins and cromolyn as possible drugs against the cytotoxicity of Aβ(31-35) and Aβ(25-35) peptides: a comparative study by advanced computer simulation methods

In this work, possible effective mechanisms of cromolyn, atorvastatin and lovastatin on the cytotoxicity of Aβ(31-35) and Aβ(25-35) peptides were investigated by classical molecular dynamics and well-tempered metadynamics simulations. The results demonstrate that all the drugs affect the behavior of the peptides, such as their ability to aggregate, and alter their secondary structures and their affinity to a particular drug. Our findings from the computed properties suggest that the best drug candidate is lovastatin. This medicine inhibits peptide aggregation, adsorbs the peptides on the surface of the drug clusters, changes the secondary structure and binds to MET₃₅, which has been seen as the reason for the toxicity of the studied peptide sequences. Moreover, lovastatin is the drug which previously has demonstrated the strongest ability to penetrate the blood-brain barrier and makes lovastatin the most promising medicine among the three investigated drugs. Atorvastatin is also seen as a potential candidate if its penetration through the blood-brain barrier could be improved. Otherwise, its properties are even better than the ones demonstrated by lovastatin. Cromolyn appears to be less interesting as an anti-aggregant from the computational data, in comparison to the two statins.

Accelerated Molecular Dynamics to Explore the Binding of Transition Metals to Amyloid-β

We report the accelerated molecular dynamics (aMD) simulation of amyloid-β (Aβ) peptides of four different lengths (16, 28, 40, and 42 residues) and their complexes when bound to Cu(II), Fe(II), or Zn(II). 600 ns equilibrated trajectory data were analyzed for each structure from three independent 200 ns aMD simulations, generating 16 aMD trajectories. We show that the presence of a metal ion leads to reduced size and decreased mobility relative to the free peptide due to the anchoring effect of the ions. The reduced mobility was shown largely to be due to the restricted movement in N-terminal residues, most notably Asp1 and His6 that are involved in the metal-ion coordination in all cases. Significant disruption of the secondary structure and patterns of salt bridge interactions arise on the coordination of metal ions. In this regard, similarities were noted between results for Zn(II) and Fe(II), whereas results for Cu(II) are more comparable to that of the free peptides. Reweighting of free energy surfaces was carried out from aMD data to identify the properties and descriptions of local minima structures.

Atomistic molecular dynamics simulations of tubulin heterodimers explain the motion of a microtubule

Microtubules are essential parts of the cytoskeleton that are built by polymerization of tubulin heterodimers into a hollow tube. Regardless that their structures and functions have been comprehensively investigated in a modern soft matter, it is unclear how properties of tubulin heterodimer influence and promote the self-assembly. A detailed knowledge of such structural mechanisms would be helpful in drug design against neurodegenerative diseases, cancer, diabetes etc. In this work atomistic molecular dynamics simulations were used to investigate the fundamental dynamics of tubulin heterodimers in a sheet and a short microtubule utilizing well-equilibrated structures. The breathing motions of the tubulin heterodimers during assembly show that the movement at the lateral interface between heterodimers (wobbling) dominates in the lattice. The simulations of the protofilament curvature agrees well with recently published experimental data, showing curved protofilaments at polymerization of the microtubule plus end. The tubulin heterodimers exposed at the microtubule minus end were less curved and displayed altered interactions at the site of sheet closure around the outmost heterodimers, which may slow heterodimer binding and polymerization, providing a potential explanation for the limited dynamics observed at the minus end.

Potential Anti-Alzheimer Agents from Guanidinyl Tryptophan Derivatives with Activities of Membrane Adhesion and Conformational Transition Inhibitions

Guanidinyl tryptophan derivatives TGN1, TGN2, TGN3, and TGN4 were synthesized, and these compounds were shown to possess in vitro inhibitory activity for amyloid aggregation in a previous study. Nevertheless, the influence of the TGN series of compounds on the binding and permeation behaviors of an Aβ monomer to the cell membranes was not elucidated. In this study, we investigated the effect of compounds in the TGN series on the behavior of an Aβ monomer regarding its toxicity toward the bilayer lipid membrane using molecular dynamics (MD) simulation. MD simulations suggest that TGN4 is a potential agent that can interfere with the movement of the Aβ monomer into the membrane. The MM-GBSA result demonstrated that TGN4 exhibits the highest affinity to the Aβ_1–42 monomer but has the lowest affinity to the bilayer. Moreover, TGN4 also contributes to a decrease in the binding affinity between the Aβ_1–42 monomer and the POPC membrane. Regarding the results of the binding mode and conformational analyses, a high number of amino-acid residues were shown to provide the binding interactions between TGN4 and the Aβ_1–42 monomer. TGN4 also reduces the conformational transition of the Aβ_1–42 monomer by means of interacting with the monomer. The present study presents molecular-level insights into how the TGN series of compounds affect the membrane adsorption and the conformational transition of the Aβ_1–42 monomer, which could be valuable for the further development of new anti-Alzheimer agents.

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

Molecular dynamics simulation for mechanism revelation of the safety and nutrition of lipids and derivatives in food: State of the art

Molecular dynamics (MD) simulation has proved to be a powerful tool in the study of proteins, nucleic acids, lipids, and carbohydrates et al. in fields of health, nutrition, and food science. In particular, MD simulation has been employed in the investigation of various lipid systems such as triglycerides, phospholipid membranes, etc. Due to the continuous updating of computing resources and the development of new MD simulation methods and force field parameters, the simulation's time and size scale of lipids system has increased by several orders of magnitude. However, MD simulation cannot be used for systems invovle chemical reactions. These greatly limit its further application in the field of lipid research. This paper reviews the progress and development of MD simulation, especially for the application of MD simulation in different lipid systems. In this paper, MD simulation and its general workflow was briefly introduced firstly. Subsequently, the application of MD simulation in various lipid systems was reviewed in-depth. Finally, the limitation and future prospects of MD simulation in lipid research were also discussed. This review provided new insights into the investigation of MD simulation, and a novel thought for lipid study. We believe that MD simulation will exhibit more and more great advantages in the investigation of lipids in the future due to the development of novlel methods.

Cation-Zwitterionic Lipid Interactions Are Affected by the Lateral Area per Lipid

Interactions of the divalent cations Ca²⁺ and Mg²⁺ with the zwitterionic lipid bilayers prepared of a fully saturated dipalmitoylphosphatidylcholine (DPPC) or a di-monounsaturated dioleoylphosphatidylcholine (DOPC) were studied by using the neutron scattering methods and molecular dynamics simulations. The effect on the bilayer structural properties confirms the direct interactions in all cases studied. The changes are observed in the bilayer thickness and lateral area. The extent of these structural changes, moreover, suggests various mechanisms of the cation-lipid interactions. First, we have observed a small difference when studying DPPC bilayers in the gel and fluid phases, with somewhat larger effects in the former case. Second, the hydration proved to be a factor in the case of DOPC bilayers, with the larger effects in the case of less hydrated systems. Most importantly, however, there was a qualitative difference between the results of the fully hydrated DOPC bilayers and the others examined. These observations then prompt us to suggest an interaction model that is plausibly governed by the lateral area of lipid, though affected indirectly also by the hydration level. Namely, when the interlipid distance is small enough to allow for the multiple lipid-ion interactions, the lipid-ion-lipid bridges are formed. The bridges impose strong attractions that increase the order of lipid hydrocarbon chains, resulting in the bilayer thickening. In the other case, when the interlipid distance extends beyond a limiting length corresponding to the area per lipid of ∼65 Ų, Mg²⁺ and Ca²⁺ continue to interact with the lipid groups by forming the separate ion-lipid pairs. As the interactions proposed affect the lipid membrane structure in the lateral direction, they may prove to play their role in other mechanisms lying within the membrane multicomponent systems and regulating for example the lipid-peptide-ion interactions.

Effects of lipid heterogeneity on model human brain lipid membranes

Cell membranes naturally contain a heterogeneous lipid distribution. However, homogeneous bilayers are commonly preferred and utilised in computer simulations due to their relative simplicity, and the availability of lipid force field parameters. Recently, experimental lipidomics data for the human brain cell membranes under healthy and Alzheimer's disease (AD) conditions were investigated, since disruption to the lipid composition has been implicated in neurodegenerative disorders, including AD [R. B. Chan et al., J. Biol. Chem., 2012, 287, 2678-2688]. In order to observe the effects of lipid complexity on the various bilayer properties, molecular dynamics simulations were used to study four membranes with increasing heterogeneity: a pure POPC membrane, a POPC and cholesterol membrane in a 1 : 1 ratio (POPC-CHOL), and to our knowledge, the first realistic models of a healthy brain membrane and an Alzheimer's diseased brain membrane. Numerous structural, interfacial, and dynamical properties, including the area per lipid, interdigitation, dipole potential, and lateral diffusion of the two simple models, POPC and POPC-CHOL, were analysed and compared to those of the complex brain models consisting of 27 lipid components. As the membranes gain heterogeneity, a number of alterations were found in the structural and dynamical properties, and more significant differences were observed in the lateral diffusion. Additionally, we observed snorkeling behaviour of the lipid tails that may play a role in the permeation of small molecules across biological membranes. In this work, atomistic description of realistic brain membrane models is provided, which can add insight towards the permeability and transport pathways of small molecules across these membrane barriers.

Effect of Amyloid-β Monomers on Lipid Membrane Mechanical Parameters-Potential Implications for Mechanically Driven Neurodegeneration in Alzheimer's Disease

Alzheimer’s disease (AD) is a neurodegenerative disease that results in memory loss and the impairment of cognitive skills. Several mechanisms of AD’s pathogenesis were proposed, such as the progressive accumulation of amyloid-β (Aβ) and τ pathology. Nevertheless, the exact neurodegenerative mechanism of the Aβ remains complex and not fully understood. This paper proposes an alternative hypothesis of the mechanism based on maintaining the neuron membrane’s mechanical balance. The incorporation of Aβ decreases the lipid membrane’s elastic properties, which eventually leads to the impairment of membrane clustering, disruption of mechanical wave propagation, and change in gamma oscillations. The first two disrupt the neuron’s ability to function correctly while the last one decreases sensory encoding and perception enabling. To begin discussing this mechanical-balance hypothesis, we measured the effect of two selected peptides, Aβ-40 and Aβ-42, as well as their fluorescently labeled modification, on membrane mechanical properties. The decrease of bending rigidity, consistent for all investigated peptides, was observed using molecular dynamic studies and experimental flicker-noise techniques. Additionally, wave propagation was investigated with molecular dynamic studies in membranes with and without incorporated neurodegenerative peptides. A change in membrane behavior was observed in the membrane system with incorporated Aβ.

Lactoferrin-Loaded PEG/PLA Block Copolymer Targeted With Anti-Transferrin Receptor Antibodies for Alzheimer Disease

Last few years, struggles have been reported to develop the nanovesicles for drug delivery via the brain-blood barrier (BBB). Novel drugs, for instance, iAβ5, are efficient to inhibit the aggregates connected to the treatment of Alzheimer disease and are being evaluated, but most of the reports reflect some drawbacks of the drugs to reach the brain in preferred concentrations owing to the less BBB penetrability of the surface dimensions. In this report, we designed and developed a new approach to enhance the transport of drug via BBB, constructed with lactoferrin (Lf)-coated polyethylene glycol-polylactide nanoparticles (Lf-PPN) with superficial monoclonal antibody-functionalized antitransferrin receptor and anti-Aβ to deliver the iAβ5 hooked on the brain. The porcine brain capillary endothelial cells were utilized as BBB typically to examine the framework efficacy and toxicity. The cellular uptake of the immuno-nanoparticles with measured conveyance of the iAβ5 peptide was significantly enhanced and associated with Lf-PPN without monoclonal antibody functionalizations.

Computer simulations of protein-membrane systems

The interactions between proteins and membranes play critical roles in signal transduction, cell motility, and transport, and they are involved in many types of diseases. Molecular dynamics (MD) simulations have greatly contributed to our understanding of protein-membrane interactions, promoted by a dramatic development of MD-related software, increasingly accurate force fields, and available computer power. In this chapter, we present available methods for studying protein-membrane systems with MD simulations, including an overview about the various all-atom and coarse-grained force fields for lipids, and useful software for membrane simulation setup and analysis. A large set of case studies is discussed.

Modelling of interactions between Aβ(25-35) peptide and phospholipid bilayers: effects of cholesterol and lipid saturation

Aggregation of amyloid beta (Aβ) peptides in neuronal membranes is a known promoter of Alzheimer’s disease. To gain insight into the molecular details of Aβ peptide aggregation and its effect on model neuronal membranes, we carried out molecular dynamics simulations of the Aβ(25–35) fragment of the amyloid precursor protein in phospholipid bilayers composed of either fully saturated or highly unsaturated lipids, in the presence or absence of cholesterol. It was found that the peptide does not penetrate through any of the considered membranes, but can reside in the headgroup region and upper part of the lipid tails showing a clear preference to a polyunsaturated cholesterol-free membrane. Due to the ordering and condensing effect upon addition of cholesterol, membranes become more rigid facilitating peptide aggregation on the surface. Except for the case of the cholesterol-free saturated lipid bilayer, the peptides have a small effect on the membrane structure and ordering. It was also found that the most “active” amino-acid for peptide–lipid and peptide–cholesterol interaction is methionine-35, followed by asparagine-27 and serine-26, which form hydrogen bonds between peptides and polar atoms of lipid headgroups. These amino acids are also primarily responsible for peptide aggregation. This work will be relevant for designing strategies to develop drugs to combat Alzheimer’s disease.

Molecular dynamics simulations of Aβ(25–35) peptides in phospholipid bilayers are carried out to investigate the effect of polyunsaturated lipids and cholesterol on aggregation of the peptides.