Extension of the Slipids Force Field to Polyunsaturated Lipids

Molecular Modeling of the Adsorption of an Egg Yolk Protein on a Water-Oil Interface

Egg yolk contains several molecular species with emulsifying properties, such as proteins and phospholipids. In particular, these molecules have both polar and nonpolar parts and thus can act as surfactants. One of the most surface-active proteins from egg yolk low-density lipoproteins is the so-called Apovitellenin-1. Experimental studies have been hindered by difficulties in isolating individual species from egg yolk lipoproteins. The purpose of this work was to assess the emulsifying properties of Apovitellenin-1 and any potential cooperative or competitive behavior in the presence of phospholipids. To do so, molecular simulations were carried out in a liquid-liquid interfacial system consisting of water and soybean oil, with varying concentrations of phospholipids and for different spatial configurations. To evaluate the conformational stability of the protein at the water-oil interface, the Gibbs free energy was computed from Metadynamics simulations as a function of the distance from the interface and of the radius of gyration. Moreover, a detailed analysis was also performed to determine which peptide residues were responsible for the protein adsorption at the oil-water interface as well as the lowering of the interfacial tension. Lastly, we combined the simulation results with a thermodynamic model to predict the interfacial tension behavior at increasing protein bulk concentration, which cannot be measured experimentally.

Molecular dynamics study on micelle-small molecule interactions: developing a strategy for an extensive comparison

Theoretical predictions of the solubilizing capacity of micelles and vesicles present in intestinal fluid are important for the development of new delivery techniques and bioavailability improvement. A balance between accuracy and computational cost is a key factor for an extensive study of numerous compounds in diverse environments. In this study, we aimed to determine an optimal molecular dynamics (MD) protocol to evaluate small-molecule interactions with micelles composed of bile salts and phospholipids. MD simulations were used to produce free energy profiles for three drug molecules (danazol, probucol, and prednisolone) and one surfactant molecule (sodium caprate) as a function of the distance from the colloid center of mass. To address the challenges associated with such tasks, we compared different simulation setups, including freely assembled colloids versus pre-organized spherical micelles, full free energy profiles versus only a few points of interest, and a coarse-grained model versus an all-atom model. Our findings demonstrate that combining these techniques is advantageous for achieving optimal performance and accuracy when evaluating the solubilization capacity of micelles. All-atom (AA) and coarse-grained (CG) umbrella sampling (US) simulations and point-wise free energy (FE) calculations were compared to their efficiency to computationally analyze the solubilization of active pharmaceutical ingredients in intestinal fluid colloids.

Alchemical Free Energy Calculations on Membrane-Associated Proteins

Membrane proteins have diverse functions within cells and are well-established drug targets. The advances in membrane protein structural biology have revealed drug and lipid binding sites on membrane proteins, while computational methods such as molecular simulations can resolve the thermodynamic basis of these interactions. Particularly, alchemical free energy calculations have shown promise in the calculation of reliable and reproducible binding free energies of protein-ligand and protein-lipid complexes in membrane-associated systems. In this review, we present an overview of representative alchemical free energy studies on G-protein-coupled receptors, ion channels, transporters as well as protein-lipid interactions, with emphasis on best practices and critical aspects of running these simulations. Additionally, we analyze challenges and successes when running alchemical free energy calculations on membrane-associated proteins. Finally, we highlight the value of alchemical free energy calculations calculations in drug discovery and their applicability in the pharmaceutical industry.

Puncturing lipid membranes: onset of pore formation and the role of hydrogen bonding in the presence of flavonoids

Products of lipid peroxidation induce detrimental structural changes in cell membranes, such as the formation of water pores, which occur in the presence of lipids with partially oxidized chains. However, the influence of another class of products, dicarboxylic acids, is still unclear. These products have greater mobility in the lipid bilayer, which enables their aggregation and the formation of favorable sites for the appearance of pores. Therefore, dodecanedioic acid (DDA) was selected as a model product. Additionally, the influence of several structurally different flavonoids on DDA aggregation via formation of hydrogen bonds with carboxyl groups was investigated. The molecular dynamics of DDA in DOPC lipid bilayer revealed the formation of aggregates extending over the hydrophobic region of the bilayer and increasing its polarity. Consequently, water penetration and the appearance of water wires was observed, representing a new step in the mechanism of pore formation. Furthermore, DDA molecules were found to interact with lipid polar groups, causing them to be buried in the bilayer. The addition of flavonoids to the system disrupted aggregate formation, resulting in the displacement of DDA molecules from the center of the bilayer. The placement of DDA and flavonoids in the lipid bilayer was confirmed by small-angle X-ray scattering. Atomic force microscopy and electron paramagnetic resonance were used to characterize the structural properties. The presence of DDA increased bilayer roughness and decreased the ordering of lipid chains, confirming its detrimental effects on the membrane surface, while flavonoids were found to reduce or reverse these changes.

Ionizable lipids penetrate phospholipid bilayers with high phase transition temperatures: perspectives from free energy calculations

The efficacies of modern gene-therapies strongly depend on their contents. At the same time the most potent formulations might not contain the best compounds. In this work we investigated the effect of phospholipids and their saturation on the binding ability of (6Z,9Z,28Z,31Z)-heptatriacont-6,9,28,31-tetraene-19-yl 4-(dimethylamino) butanoate (DLin-MC3-DMA) to model membranes at the neutral pH. We discovered that DLin-MC3-DMA has affinity to the most saturated monocomponent lipid bilayer 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and an aversion to the unsaturated one 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC). The preference to a certain membrane was also well-correlated to the phase transition temperatures of phospholipid bilayers, and to their structural and dynamical properties. Additionally, in the case of the presence of DLin-MC3-DMA in the membrane with DOPC the ionizable lipid penetrated it, which indicates possible synergistic effects. Comparisons with other ionizable lipids were performed using a model lipid bilayer of 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC). Particularly, the lipids heptadecan-9-yl 8-[2-hydroxyethyl-(6-oxo-6-undecoxyhexyl)amino]octanoate (SM-102) and [(4-Hydroxybutyl) azanediyl] di(hexane-6,1-diyl) bis(2-hexyldecanoate) (ALC-0315) from modern mRNA-vaccines against COVID-19 were investigated and force fields parameters were derived for those new lipids. It was discovered that ALC-0315 binds strongest to the membrane, while DLin-MC3-DMA is not able to reside in the bilayer center. The ability to penetrate the membrane POPC by SM-102 and ALC-0315 can be related to their saturation, comparing to DLin-MC3-DMA.

Unraveling membrane properties at the organelle-level with LipidDyn

•

Lipidomics of organelles could be used to design models for molecular simulations.

•

The bottleneck is the analysis and rationalization of the data from simulations.

•

LipidDyn is an automated pipeline to streamline the analyses of lipid bilayers.

•

LipidDyn allows to collect analysis in a non-time-consuming and reproducible manner.

•

We applied LipidDyn to different case studies to illustrate its potential.

Cellular membranes are formed from different lipids in various amounts and proportions depending on the subcellular localization. The lipid composition of membranes is sensitive to changes in the cellular environment, and its alterations are linked to several diseases. Lipids not only form lipid-lipid interactions but also interact with other biomolecules, including proteins.

Molecular dynamics (MD) simulations are a powerful tool to study the properties of cellular membranes and membrane-protein interactions on different timescales and resolutions. Over the last few years, software and hardware for biomolecular simulations have been optimized to routinely run long simulations of large and complex biological systems. On the other hand, high-throughput techniques based on lipidomics provide accurate estimates of the composition of cellular membranes at the level of subcellular compartments. Lipidomic data can be analyzed to design biologically relevant models of membranes for MD simulations. Similar applications easily result in a massive amount of simulation data where the bottleneck becomes the analysis of the data. In this context, we developed LipidDyn, a Python-based pipeline to streamline the analyses of MD simulations of membranes of different compositions. Once the simulations are collected, LipidDyn provides average properties and time series for several membrane properties such as area per lipid, thickness, order parameters, diffusion motions, lipid density, and lipid enrichment/depletion. The calculations exploit parallelization, and the pipeline includes graphical outputs in a publication-ready form. We applied LipidDyn to different case studies to illustrate its potential, including membranes from cellular compartments and transmembrane protein domains. LipidDyn is available free of charge under the GNU General Public License from https://github.com/ELELAB/LipidDyn.

Effect of encapsulated protein on the dynamics of lipid sponge phase: a neutron spin echo and molecular dynamics simulation study

Lipid membranes are highly mobile systems with hierarchical, time and length scale dependent, collective motions including thickness fluctuations, undulations, and topological membrane changes, which play an important role in membrane interactions. In this work we have characterised the effect of encapsulating two industrially important enzymes, β-galactosidase and aspartic protease, in lipid sponge phase nanoparticles on the dynamics of the lipid membrane using neutron spin echo (NSE) spectroscopy and molecular dynamics (MD) simulations. From NSE, reduced membrane dynamics were observed upon enzyme encapsulation, which were dependent on the enzyme concentration and type. By fitting the intermediate scattering functions (ISFs) with a modified Zilman and Granek model including nanoparticle diffusion, an increase in membrane bending rigidity was observed, with a larger effect for β-galactosidase than aspartic protease at the same concentration. MD simulations for the system with and without aspartic protease showed that the lipids relax more slowly in the system with protein due to the replacement of the lipid carbonyl-water hydrogen bonds with lipid-protein hydrogen bonds. This indicates that the most likely cause of the increase in membrane rigidity observed in the NSE measurements was dehydration of the lipid head groups. The dynamics of the protein itself were also studied, which showed a stable secondary structure of protein over the simulation, indicating no unfolding events occurred.

Membrane Mixer: A Toolkit for Efficient Shuffling of Lipids in Heterogeneous Biological Membranes

Molecular dynamics (MD) simulations of biological membranes have achieved such levels of sophistication that are commonly used to predict unresolved structures and various properties of lipids and to substantiate experimental data. While achieving sufficient sampling of lipid dynamics remains a major challenge, a commonly used method to improve lipid sampling, e.g., in terms of specific interactions with membrane-associated proteins, is to randomize the initial arrangement of lipid constituents in multiple replicas of simulations, without changing the overall lipid composition of the membrane of interest. Here, we introduce a method that can rapidly generate multiple replicas of lipid bilayers with different spatial and conformational configurations for any given lipid composition. The underlying algorithm, which allows one to shuffle lipids at any desired level, relies on the application of an external potential, here referred to as the "carving potential", that removes clashes/entanglements before lipid positions are exchanged (shuffled), thereby minimizing the energy penalty due to abrupt lipid repositioning. The method is implemented as "Membrane Mixer Plugin (MMP) 1.0" in VMD, with a convenient graphical user interface that guides the user in setting various options and parameters. The plugin is fully automated and generates new membrane replicas more rapidly and conveniently than other analogous tools. The plugin and its capabilities introduced here can be extended to include additional features in future versions.

Molecular Organization of a Raft-like Domain in a Polyunsaturated Phospholipid Bilayer: A Supervised Machine Learning Analysis of Molecular Dynamics Simulations

Numerous health benefits are associated with omega-3 polyunsaturated fatty acids (n-3 PUFA) consumed in fish oils. An understanding of the mechanism remains elusive. The plasma membrane as a site of action is the focus in this study. With large-scale all-atom MD simulations run on a model membrane (1050 lipid molecules), we observed the evolution over time (6 μs) of a circular (raft-like) domain composed of N-palmitoylsphingomyelin (PSM) and cholesterol embedded into a surrounding (non-raft) patch composed of polyunsaturated 1-palmitoyl-2-docosahexaenoylphosphatylcholine (PDPC) (1:1:1 mol). A supervised machine learning algorithm was developed to characterize the migration of each lipid based on molecular conformation and the local environment. PDPC molecules were seen to infiltrate the ordered raft-like domain in a small amount, while a small concentration of PSM and cholesterol molecules was seen to migrate into the disordered non-raft region. Enclosing the raft-like domain, a narrow (∼2 nm in width) interfacial zone composed of PDPC, PSM, and cholesterol that buffers the substantial difference in order (ΔS_CD ≈ 0.12) between raft-like and non-raft environments was seen to form. Our results suggest that n-3 PUFA regulate the architecture of lipid rafts enriched in sphingolipids and cholesterol with a minimal effect on order within their interior in membranes.

Site of Cholesterol Oxidation Impacts Its Localization and Domain Formation in the Neuronal Plasma Membrane

Cholesterol is integral to the structure of mammalian cell membranes. Oxidation of cholesterol alters how it behaves in the membrane and influences the membrane biophysical properties. Elevated levels of oxidized cholesterol are associated with neurodegenerative diseases such as Alzheimer's disease, Parkinson's disease, multiple sclerosis, and Huntington's disease. Previous work has investigated the impact of oxidized cholesterol in the context of simple model membrane systems. However, there is a growing body of literature that shows that complex membranes possessing physiological phospholipid distributions have different properties from those of binary or trinary model membranes. In the current work, the impact of oxidized cholesterol on the biophysical properties of a complex neuronal plasma membrane is investigated using coarse-grained Martini molecular dynamics simulations. Comparison of the native neuronal membrane to neuronal membranes containing 10% tail-oxidized or 10% head-oxidized cholesterol shows that the site of oxidization changes the behavior of the oxidized cholesterol in the membrane. Furthermore, species-specific domain formation is observed between each oxidized cholesterol and minor lipid classes. Although both tail-oxidized and head-oxidized cholesterols modulate the biophysical properties of the membrane, smaller changes are observed in the complex neuronal membrane than seen in the previous work on simple binary or trinary model membranes. This work highlights the presence of compensatory effects of lipid diversity in the complex neuronal membrane. Overall, this study improves our molecular-level understanding of the effects of oxidized cholesterol on the properties of neuronal tissue and emphasizes the importance of studying membranes with realistic lipid compositions.

Amyloid-β peptide dimers undergo a random coil to β-sheet transition in the aqueous phase but not at the neuronal membrane

Mounting evidence suggests that the neuronal cell membrane is the main site of oligomer-mediated neuronal toxicity of amyloid-β peptides in Alzheimer's disease. To gain a detailed understanding of the mutual interference of amyloid-β oligomers and the neuronal membrane, we carried out microseconds of all-atom molecular dynamics (MD) simulations on the dimerization of amyloid-β (Aβ)42 in the aqueous phase and in the presence of a lipid bilayer mimicking the in vivo composition of neuronal membranes. The dimerization in solution is characterized by a random coil to β-sheet transition that seems on pathway to amyloid aggregation, while the interactions with the neuronal membrane decrease the order of the Aβ42 dimer by attenuating its propensity to form a β-sheet structure. The main lipid interaction partners of Aβ42 are the surface-exposed sugar groups of the gangliosides GM1. As the neurotoxic activity of amyloid oligomers increases with oligomer order, these results suggest that GM1 is neuroprotective against Aβ-mediated toxicity.

Investigating the lipid fingerprint of SLC6 neurotransmitter transporters: a comparison of dDAT, hDAT, hSERT, and GlyT2

The local lipid annulus, or "fingerprint", of four SLC6 transporters (dDAT, hDAT, hSERT, and GlyT2) embedded in a complex neuronal membrane were compared and characterised using molecular dynamics. Our analysis included the development of new tools to improve membrane leaflet detection and the analysis of leaflet-dependent properties. Overall, the lipid fingerprints of the four transporters are comprised of similar lipids when grouped by headgroup or tail saturation. The enrichment and depletion of specific lipids, including sites of cholesterol contacts, varies between transporters. The subtle differences in lipid fingerprints results in varying membrane biophysical properties near the transporter. Our results highlight that the lipid-fingerprint of SLC6 transporters in complex membranes is highly dependent on membrane composition. Our results further characterize how the presence and identity of membrane proteins affects the complex interplay of lipid-protein interactions, influencing the local lipid environment and membrane biophysical properties.

DOPC versus DOPE as a helper lipid for gene-therapies: molecular dynamics simulations with DLin-MC3-DMA

Ionizable lipids are important compounds of modern therapeutic lipid nano-particles (LNPs). One of the most promising ionizable lipids (or amine lipids) is DLin-MC3-DMA. Depending on their pharmaceutical application these LNPs can also contain various helper lipids, such as phospho- and pegylated lipids, cholesterol and nucleic acids as a cargo. Due to their complex compositions the structures of these therapeutics have not been refined properly. Therefore, the role of each lipid in the pharmacological properties of LNPs has not been determined. In this work an atomistic model for the neutral form of DLin-MC3-DMA was derived and all-atom molecular dynamics (MD) simulations were carried out in order to investigate the effect of the phospholipid headgroup on the possible properties of the shell-membranes of LNPs. Bilayers containing either DOPC or DOPE lipids at two different ratios of DLin-MC3-DMA (5 mol% and 15 mol%) were constructed and simulated at neutral pH 7.4. The results from the analysis of MD trajectories revealed that DOPE lipid headgroups associated strongly with lipid tails and carbonyl oxygens of DLin-MC3-DMA, while for DOPC lipid headgroups no significant associations were observed. Furthermore, the strong associations between DOPE and DLin-MC3-DMA result in the positioning of DLin-MC3-DMA at the surface of the membrane. Such an interplay between the lipids slows down the lateral diffusion of all simulated bilayers, where a more dramatic decrease of the diffusion rate is observed in membranes with DOPE. This can explain the low water penetration of lipid bilayers with phosphatidylethanolamines and, probably, can relate to the bad transfection properties of LNPs with DOPE and DLin-MC3-DMA.

Flavonol clustering in model lipid membranes: DSC, AFM, force spectroscopy and MD simulations study

We here report on flavonols (myricetin (MCE) and its glycoside myricitrin (MCI)) - 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) membrane interactions focusing on the effects of flavonol clustering on the membrane thermotropic and nanomechanical properties. Atomic force microscopy (AFM), force spectroscopy (FS) and differential scanning calorimetry (DSC) together with molecular dynamics (MD) simulations provided a consistent picture of flavonol - DMPC membrane interactions. DMPC membrane as a supported lipid bilayer preserved its integrity even at higher flavonol molar fraction x. When present at x = 0.1 - 0.3, MCE and MCI both slightly improve DMPC bilayer fluidity which is evidenced by the decrease in the main phase transition temperature T_m. MCE is found within the interior of the bilayer, while MCI incorporates in the head group-water interface region. AFM and FS confirmed clusters as protrusions with an average height of 0.012 μm and average diameters of 0.60 and 0.24 μm for MCE and MCI clusters, respectively. The average membrane thickness in DMPC fluid phase decreases for 7% at x_MCE = 0.30, while only 4% at x_MCI = 0.27. The induced membrane changes are dependent on the chemical and physical properties of inserted flavonols. The hypothesis regarding the tendency of flavonol to clustering in membranes by increasing flavonol molar fraction has been confirmed.

Structure of Supramers Formed by the Amphiphile Biotin-CMG-DOPE

The synthetic function-spacer-lipid (FSL) amphiphile biotin-CMG-DOPE is widely used for delicate ligation of living cells with biotin residues under physiological conditions. Since this molecule has an "apolar-polar-hydrophobic" gemini structure, the supramolecular organization is expected to differ significantly from the classical micelle. Its organization is investigated with experimental methods and molecular dynamics simulations (MDS). Although the linear length of a single biotin-CMG-DOPE molecule is 9.5 nm, the size of the dominant supramer globule is only 14.6 nm. Investigations found that while the DOPE tails form a hydrophobic core, the polar CMG spacer folds back upon itself and predominantly places the biotin reside inside the globule or planar layer. MDS demonstrates that <10 % of biotin residues on the highly water dispersible globules and only 1 % of biotin residues in layer coatings are in an linear conformation and exposing biotin into the aqueous medium. This explains why in biotin-CMG-DOPE apolar biotin residues both in water dispersible globules and coatings on solid surfaces are still capable of interacting with streptavidin.

Contributions and Limitations of Biophysical Approaches to Study of the Interactions between Amphiphilic Molecules and the Plant Plasma Membrane

Some amphiphilic molecules are able to interact with the lipid matrix of plant plasma membranes and trigger the immune response in plants. This original mode of perception is not yet fully understood and biophysical approaches could help to obtain molecular insights. In this review, we focus on such membrane-interacting molecules, and present biophysically grounded methods that are used and are particularly interesting in the investigation of this mode of perception. Rather than going into overly technical details, the aim of this review was to provide to readers with a plant biochemistry background a good overview of how biophysics can help to study molecular interactions between bioactive amphiphilic molecules and plant lipid membranes. In particular, we present the biomimetic membrane models typically used, solid-state nuclear magnetic resonance, molecular modeling, and fluorescence approaches, because they are especially suitable for this field of research. For each technique, we provide a brief description, a few case studies, and the inherent limitations, so non-specialists can gain a good grasp on how they could extend their toolbox and/or could apply new techniques to study amphiphilic bioactive compound and lipid interactions.

Cyclopentane rings in hydrophobic chains of a phospholipid enhance the bilayer stability to electric breakdown

Archaeal lipids ensure unprecedented stability of archaea membranes in extreme environments. Here, we incorporate a characteristic structural feature of an archaeal lipid, the cyclopentane ring, into hydrocarbon chains of a short-chain (C12) phosphatidylcholine to explore whether the insertion would allow such a lipid (1,2-di-(3-(3-hexylcyclopentyl)-propanoate)-sn-glycero-3-phosphatidylcholine, diC12cp-PC) to form stable bilayers at room temperature. According to fluorescence-based assays, in water diC12cp-PC formed liquid-crystalline bilayers at room temperature. Liposomes produced from diC12cp-PC retained calcein for over a week when stored at +4 °C. diC12cp-PC could also form model bilayer lipid membranes that were by an order of magnitude more stable to electrical breakdown than egg PC membranes. Molecular dynamics simulation showed that the cyclopentane fragment fixes five carbon atoms (or four C-C bonds), which is compensated by the higher mobility of the rest of the chain. This was found to be the reason for the remarkable stability of the diC12cp-PC bilayer: restricted conformational mobility of a chain segment increases the membrane bending modulus (compared to a normal hydrocarbon chain of the same length). Here, higher stiffness practically does not affect the line tension of a membrane pore edge. Rather it makes it more difficult for diC12cp-PC to rearrange in order to line the edge of a hydrophilic pore; therefore, fewer pores are formed.

Archaeal cyclopentane fragment in a surfactant's hydrophobic tail decreases the Krafft point

Archaea are prokaryotic microorganisms famous for their ability to adapt to extreme environments, including low and high temperatures. Archaeal lipids often are macrocycles with two polar heads and a hydrophobic core that contains methyl groups and in-line cycles. Here we present the design of novel general-purpose surfactants that have inherited features of archaeal lipids. These are C12 and C14 carboxylic acids containing in-line cyclopentanes. The cyclopentanes disturb the chain packing, which results in remarkable expansion of the operational range of the surfactant into the low-temperature region. We report synthesis and properties of these novel archaea-like surfactants and details of their chain packing derived from thermodynamics model predictions, molecular dynamics simulations, and experimental data on CMC and Krafft points.

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.

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

Aggregation of amyloid-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}β (Aβ) peptides, cleaved from the amyloid precursor protein, is known as a precursor of the Alzheimer’s disease (AD). It is also known that Alzheimer’s disease is characterized by a substantial decrease of the amount of polyunsaturated lipids in the neuronal membranes of the frontal gray matter. To get insight into possible interconnection of these phenomena, we have carried out molecular dynamics simulations of two fragments of A\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}β peptide, A\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_{1-28}$$\end{document}1-28 and A\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_{26-40}$$\end{document}26-40, in four different lipid bilayers: two monocomponent ones (14:0-14:0 PC, 18:0-22:6 PC), and two bilayers containing mixtures of 18:0-18:0 PE, 22:6-22:6 PE, 16:0-16:0 PC and 18:1-18:1 PC lipids of composition mimicking neuronal membranes in a “healthy” and “AD” brain. The simulations showed that the presence of lipids with highly unsaturated 22:6cis fatty acids chains strongly affects the interaction of amyloid-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}β peptides with lipid membranes. The polyunsaturated lipids cause stronger adsorption of A\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}β-peptides by the membrane and lead to weaker binding between peptides when the latter form aggregates. This difference in the behaviour observed in monocomponent bilayers is propagated in a similar fashion to the mixed membranes mimicking composition of neuronal membranes in “healthy” and “AD” brains, with “healthy” membrane having higher fraction of polyunsaturated lipids. Our simulations give strong indication that it can be physical–chemical background of the interconnection between amyloid fibrillization causing Alzheimer’s disease, and content of polyunsaturated lipids in the neuronal membranes.

Computational Modeling of Realistic Cell Membranes

Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.

Amyloid-β(29-42) Dimeric Conformations in Membranes Rich in Omega-3 and Omega-6 Polyunsaturated Fatty Acids

The omega-3 and omega-6 polyunsaturated fatty acids are two important components of cell membranes in human brains. When incorporated into phospholipids, omega-3 slows the progression of Alzheimer's disease (AD), whereas omega-6 is linked to increased risk of AD. Little is known on the amyloid-β (Aβ) conformations in membranes rich in omega-3 and omega-6 phospholipids. Herein, the structural properties of the Aβ_29-42 dimer embedded in both fatty acid membranes were comparatively studied to a 1-palmitoyl-2-oleoyl- sn-glycero-3-phosphocholine (POPC) bilayer using all-atom molecular dynamics (MD) simulations. Starting from α-helix, both omega-6 and omega-3 membranes promote new orientations and conformations of the dimer, in agreement with the observed dependence of Aβ production upon addition of these two fatty acids. This conformational result is corroborated by atomistic MD simulations of the dimer of the 99 amino acid C-terminal fragment of amyloid precursor protein spanning the residues 15-55. Starting from β-sheet, omega-6 membrane promotes helical and disordered structures of Aβ_29-42 dimer, whereas omega-3 membrane preserves the β-sheet structures differing however from those observed in POPC. Remarkably, the mixture of the two fatty acids and POPC depicts another conformational ensemble of the Aβ_29-42 dimer. This finding demonstrates that variation in the abundance of the molecular phospholipids, which changes with age, modulates membrane-embedded Aβ oligomerization.

Emerging Diversity in Lipid-Protein Interactions

Membrane lipids interact with proteins in a variety of ways, ranging from providing a stable membrane environment for proteins to being embedded in to detailed roles in complicated and well-regulated protein functions. Experimental and computational advances are converging in a rapidly expanding research area of lipid–protein interactions. Experimentally, the database of high-resolution membrane protein structures is growing, as are capabilities to identify the complex lipid composition of different membranes, to probe the challenging time and length scales of lipid–protein interactions, and to link lipid–protein interactions to protein function in a variety of proteins. Computationally, more accurate membrane models and more powerful computers now enable a detailed look at lipid–protein interactions and increasing overlap with experimental observations for validation and joint interpretation of simulation and experiment. Here we review papers that use computational approaches to study detailed lipid–protein interactions, together with brief experimental and physiological contexts, aiming at comprehensive coverage of simulation papers in the last five years. Overall, a complex picture of lipid–protein interactions emerges, through a range of mechanisms including modulation of the physical properties of the lipid environment, detailed chemical interactions between lipids and proteins, and key functional roles of very specific lipids binding to well-defined binding sites on proteins. Computationally, despite important limitations, molecular dynamics simulations with current computer power and theoretical models are now in an excellent position to answer detailed questions about lipid–protein interactions.

Simulation Best Practices for Lipid Membranes [Article v1.0]

We establish a reliable and robust standardization of settings for practical molecular dynamics (MD) simulations of pure and mixed (single- and multi-component) lipid bilayer membranes. In lipid membranes research, particle-based molecular simulations are a powerful tool alongside continuum theory, lipidomics, and model, in vitro, and in vivo experiments. Molecular simulations can provide precise and reproducible spatiotemporal (atomic- and femtosecond-level) information about membrane structure, mechanics, thermodynamics, kinetics, and dynamics. Yet the simulation of lipid membranes can be a daunting task, given the uniqueness of lipid membranes relative to conventional liquid-liquid and solid-liquid interfaces, the immense and complex thermodynamic and statistical mechanical theory, the diversity of multiscale lipid models, limitations of modern computing power, the difficulty and ambiguity of simulation controls, finite size effects, competitive continuum simulation alternatives, and the desired application, including vesicle experiments and biological membranes. These issues can complicate an essential understanding of the field of lipid membranes, and create major bottlenecks to simulation advancement. In this article, we clarify these issues and present a consistent, thorough, and user-friendly framework for the design of state-of-the-art lipid membrane MD simulations. We hope to allow early-career researchers to quickly overcome common obstacles in the field of lipid membranes and reach maximal impact in their simulations.

Multiscale Modeling and Simulation Approaches to Lipid-Protein Interactions

Lipid membranes play a crucial role in living systems by compartmentalizing biological processes and forming a barrier between these processes and the environment. Naturally, a large apparatus of biomolecules is responsible for construction, maintenance, transport, and degradation of these lipid barriers. Additional classes of biomolecules are tasked with transport of specific substances or transduction of signals from the environment across lipid membranes. In this article, we intend to describe a set of techniques that enable one to build accurate models of lipid systems and their associated proteins, and to simulate their dynamics over a variety of time and length scales. We discuss the methods and challenges that allow us to derive structural, mechanistic, and thermodynamic information from these models. We also show how these models have recently been applied in research to study some of the most complex lipid-protein systems to date, including bacterial and viral envelopes, neuronal membranes, and mammalian signaling systems.

Cholesterol in phospholipid bilayers: positions and orientations inside membranes with different unsaturation degrees

Cholesterol is an essential component of all animal cell membranes and plays an important role in maintaining the membrane structure and physical-chemical properties necessary for correct cell functioning. The presence of cholesterol is believed to be responsible for domain formation (lipid rafts) due to different interactions of cholesterol with saturated and unsaturated lipids. In order to get detailed atomistic insight into the behaviour of cholesterol in bilayers composed of lipids with varying degrees of unsaturation, we have carried out a series of molecular dynamics simulations of saturated and polyunsaturated lipid bilayers with different contents of cholesterol, as well as well-tempered metadynamics simulations with a single cholesterol molecule in these bilayers. From these simulations we have determined distributions of cholesterol across the bilayer, its orientational properties, free energy profiles, and specific interactions of molecular groups able to form hydrogen bonds. Both molecular dynamics and metadynamics simulations showed that the most unsaturated bilayer with 22:6 fatty acid chains shows behaviour which is most different from other lipids. In this bilayer, cholesterol is relatively often found in a "flipped" configuration with the hydroxyl group oriented towards the membrane middle plane. This bilayer has also the highest (least negative) binding free energy among liquid phase bilayers, and the lowest reorientation barrier. Furthermore, cholesterol molecules in this bilayer are often found to form head-to-tail contacts which may lead to specific clustering behaviour. Overall, our simulations support ideas that there can be a subtle interconnection between the contents of highly unsaturated fatty acids and cholesterol, deficiency or excess of each of them is related to many human afflictions and diseases.

Molecular Dynamics Simulations of Membrane Proteins: An Overview

Simulations of membrane proteins have been rising in popularity in the past decade. Advancements in technology and force fields made it possible to simulate behavior of membrane proteins. Membrane protein simulations can now be used as supporting evidence for experimental findings, for elucidating protein mechanisms, and validating protein crystal structures. Unrelated to experimental data, these simulations can also serve to investigate larger scale processes like protein sorting, protein-membrane interactions, and more. In this review, the history as well as the state-of-the-art methodologies in membrane protein simulations will be summarized. An emphasis will be put on how to set up the system and on the current models for the different components of the simulation system. An overview of the available tools for membrane protein simulation will be given, and current limitations and prospects will also be discussed.

Schistosomiasis: Snail-vector control, molecular modelling and dynamic studies of bioactive N-acetylglycoside saponins from Tetrapleura tetraptera

Schistosomiasis, a chronic neglected tropical disease caused by the Schistosoma spp. parasite, is associated with disabling patient symptoms. The new focus of the WHO roadmap on 'transmission control, wherever possible' offers drug development opportunities for intermediate-host control to prevent human-to-snail-to-human parasite transmission. Reports on the analysis of the impact of 'chemical-based mollusciciding' have concluded that constant application of molluscicides may contribute significantly towards the elimination of schistosomiasis in endemic areas. In South-Western Nigeria, Tetrapleura tetraptera is a tree whose fruit has been widely used in snail vector control. The presence of molluscicidal N-acetyl triterpene glycosides in the fruit has been reported. In this study, a bioactivity-directed fractionation of the fruit extract was performed to isolate the most potent molluscicidal saponin from the fruit. In an attempt to provide mechanistic insight into the observed activity, in silico screening was performed, profiling the molluscicidal N-acetyl triterpene glycosides reported from the fruit against two potential therapeutic targets in the mollusk used, NADH-ubiquinone oxidoreductase (NAD1) and retinoid X receptor. The docking predicted binary complexes of the saponins, which were subjected to explicit solvent conformational sampling from which patterns of structural stability were obtained. The binding energies alone did not account for the potency of the saponins indicating the influence of other factor like pharmacokinetic parameters. The study concluded that there is a preferential suitability of ND1's MWFE site for the rational design and development of novel molluscicidal agents.