An overview of molecular dynamics simulations of oxidized lipid systems, with a comparison of ELBA and MARTINI force fields for coarse grained lipid simulations

CHARMM GUI Membrane Builder for oxidized phospholipid membrane modeling and simulation

Oxidative stress leads to the production of oxidized phospholipids (oxPLs) that modulate the biophysical properties of phospholipid monolayers and bilayers. As many immune cells are responsible for surveilling cells and tissues for the presence of oxPLs, oxPL-dependent mechanisms have been suggested as targets for treating chronic kidney disease, atherosclerosis, diabetes, and cancer metastasis. This review details recent experimental and computational studies that characterize oxPLs' behaviors in various monolayers and bilayers. These studies investigate how the tail length and polar functional groups of OxPLs impact membrane properties, how oxidized membranes can be stabilized, and how membrane integrity is generally affected by oxidized lipids. In addition, for oxPL-containing membrane modeling and simulation, CHARMM-GUI Membrane Builder has been extended to support a variety of oxPLs, accelerating the simulation system building process for these biologically relevant lipid bilayers.

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.

Effect of lipid oxidation on the channel properties of Cx26 hemichannels: A molecular dynamics study

Intercellular communication plays a crucial role in cancer, as well as other diseases, such as inflammation, tissue degeneration, and neurological disorders. One of the proteins responsible for this, are connexins (Cxs), which come together to form a hemichannel. When two hemichannels of opposite cells interact with each other, they form a gap junction (GJ) channel, connecting the intracellular space of these cells. They allow the passage of ions, reactive oxygen and nitrogen species (RONS), and signaling molecules from the interior of one cell to another cell, thus playing an essential role in cell growth, differentiation, and homeostasis. The importance of GJs for disease induction and therapy development is becoming more appreciated, especially in the context of oncology. Studies have shown that one of the mechanisms to control the formation and disruption of GJs is mediated by lipid oxidation pathways, but the underlying mechanisms are not well understood. In this study, we performed atomistic molecular dynamics simulations to evaluate how lipid oxidation influences the channel properties of Cx26 hemichannels, such as channel gating and permeability. Our results demonstrate that the Cx26 hemichannel is more compact in the presence of oxidized lipids, decreasing its pore diameter at the extracellular side and increasing it at the amino terminus domains, respectively. The permeability of the Cx26 hemichannel for water and RONS molecules is higher in the presence of oxidized lipids. The latter may facilitate the intracellular accumulation of RONS, possibly increasing oxidative stress in cells. A better understanding of this process will help to enhance the efficacy of oxidative stress-based cancer treatments.

Effects of Nitro-Oxidative Stress on Biomolecules: Part 1-Non-Reactive Molecular Dynamics Simulations

Plasma medicine, or the biomedical application of cold atmospheric plasma (CAP), is an expanding field within plasma research. CAP has demonstrated remarkable versatility in diverse biological applications, including cancer treatment, wound healing, microorganism inactivation, and skin disease therapy. However, the precise mechanisms underlying the effects of CAP remain incompletely understood. The therapeutic effects of CAP are largely attributed to the generation of reactive oxygen and nitrogen species (RONS), which play a crucial role in the biological responses induced by CAP. Specifically, RONS produced during CAP treatment have the ability to chemically modify cell membranes and membrane proteins, causing nitro-oxidative stress, thereby leading to changes in membrane permeability and disruption of cellular processes. To gain atomic-level insights into these interactions, non-reactive molecular dynamics (MD) simulations have emerged as a valuable tool. These simulations facilitate the examination of larger-scale system dynamics, including protein-protein and protein-membrane interactions. In this comprehensive review, we focus on the applications of non-reactive MD simulations in studying the effects of CAP on cellular components and interactions at the atomic level, providing a detailed overview of the potential of CAP in medicine. We also review the results of other MD studies that are not related to plasma medicine but explore the effects of nitro-oxidative stress on cellular components and are therefore important for a broader understanding of the underlying processes.

A theoretical model of dietary lipid variance as the origin of primary ciliary dysfunction in preeclampsia

Serving as the cell's key interface in communicating with the outside world, primary cilia have emerged as an area of multidisciplinary research interest over the last 2 decades. Although the term "ciliopathy" was first used to describe abnormal cilia caused by gene mutations, recent studies focus on abnormalities of cilia that are found in diseases without clear genetic antecedents, such as obesity, diabetes, cancer, and cardiovascular disease. Preeclampsia, a hypertensive disease of pregnancy, is intensely studied as a model for cardiovascular disease partially due to many shared pathophysiologic elements, but also because changes that develop over decades in cardiovascular disease arise in days with preeclampsia yet resolve rapidly after delivery, thus providing a time-lapse view of the development of cardiovascular pathology. As with genetic primary ciliopathies, preeclampsia affects multiple organ systems. While aspirin delays the onset of preeclampsia, there is no cure other than delivery. The primary etiology of preeclampsia is unknown; however, recent reviews emphasize the fundamental role of abnormal placentation. During normal embryonic development, trophoblastic cells, which arise from the outer layer of the 4-day-old blastocyst, invade the maternal endometrium and establish extensive placental vascular connections between mother and fetus. In primary cilia of trophoblasts, Hedgehog and Wnt/catenin signaling operate upstream of vascular endothelial growth factor to advance placental angiogenesis in a process that is promoted by accessible membrane cholesterol. In preeclampsia, impaired proangiogenic signaling combined with an increase in apoptotic signaling results in shallow invasion and inadequate placental function. Recent studies show primary cilia in preeclampsia to be fewer in number and shortened with functional signaling abnormalities. Presented here is a model that integrates preeclampsia lipidomics and physiology with the molecular mechanisms of liquid-liquid phase separation in model membrane studies and the known changes in human dietary lipids over the last century to explain how changes in dietary lipids might reduce accessible membrane cholesterol and give rise to shortened cilia and defects in angiogenic signaling, which underlie placental dysfunction of preeclampsia. This model offers a possible mechanism for non-genetic dysfunction in cilia and proposes a proof-of-concept study to treat preeclampsia with dietary lipids.

CADD, AI and ML in drug discovery: A comprehensive review

Computer-aided drug design (CADD) is an emerging field that has drawn a lot of interest because of its potential to expedite and lower the cost of the drug development process. Drug discovery research is expensive and time-consuming, and it frequently took 10-15 years for a drug to be commercially available. CADD has significantly impacted this area of research. Further, the combination of CADD with Artificial Intelligence (AI), Machine Learning (ML), and Deep Learning (DL) technologies to handle enormous amounts of biological data has reduced the time and cost associated with the drug development process. This review will discuss how CADD, AI, ML, and DL approaches help identify drug candidates and various other steps of the drug discovery process. It will also provide a detailed overview of the different in silico tools used and how these approaches interact.

Directly imaging emergence of phase separation in peroxidized lipid membranes

Lipid peroxidation is a process which is key in cell signaling and disease, it is exploited in cancer therapy in the form of photodynamic therapy. The appearance of hydrophilic moieties within the bilayer's hydrocarbon core will dramatically alter the structure and mechanical behavior of membranes. Here, we combine viscosity sensitive fluorophores, advanced microscopy, and X-ray diffraction and molecular simulations to directly and quantitatively measure the bilayer's structural and viscoelastic properties, and correlate these with atomistic molecular modelling. Our results indicate an increase in microviscosity and a decrease in the bending rigidity upon peroxidation of the membranes, contrary to the trend observed with non-oxidized lipids. Fluorescence lifetime imaging microscopy and MD simulations give evidence for the presence of membrane regions of different local order in the oxidized membranes. We hypothesize that oxidation promotes stronger lipid-lipid interactions, which lead to an increase in the lateral heterogeneity within the bilayer and the creation of lipid clusters of higher order.

Effect of Disulfide Bridge on the Binding of SARS-CoV-2 Fusion Peptide to Cell Membrane: A Coarse-Grained Study

In this paper, we present the parameterization of the CAVS coarse-grained (CG) force field for 20 amino acids, and our CG simulations show that the CAVS force field could accurately predict the amino acid tendency of the secondary structure. Then, we used the CAVS force field to investigate the binding of a severe acute respiratory syndrome-associated coronavirus fusion peptide (SARS-CoV-2 FP) to a phospholipid bilayer: a long FP (FP-L) containing 40 amino acids and a short FP (FP-S) containing 26 amino acids. Our CAVS CG simulations displayed that the binding affinity of the FP-L to the bilayer is higher than that of the FP-S. We found that the FP-L interacted more strongly with membrane cholesterol than the FP-S, which should be attributed to the stable helical structure of the FP-L at the C-terminus. In addition, we discovered that the FP-S had one major and two minor membrane-bound states, in agreement with previous all-atom molecular dynamics (MD) studies. However, we found that both the C-terminal and N-terminal amino acid residues of the FP-L can strongly interact with the bilayer membrane. Furthermore, we found that the disulfide bond formed between Cys840 and Cys851 stabilized the helices of the FP-L at the C-terminus, enhancing the interaction between the FP-L and the bilayer membrane. Our work indicates that the stable helical structure is crucial for binding the SARS-CoV-2 FP to cell membranes. In particular, the helical stability of FP should have a significant influence on the FP-membrane binding.

Coarse-grained description of monounsaturated peroxidized phospholipidbilayers

Lipid peroxides result from a reaction between cis-unsaturated lipid chains and singlet oxygen molecules leading to the addition of a peroxide OOH side group next to the acyl-chain double bond. It is now established that HP-POPC (hydroperoxidized POPC) molecules form stable, thin and laterally expanded bilayers. The difference in structural organisation arises from the hydrophilic character of the OOH side group that has strong affinity with the water interface region, leading to significant reorganisation of the bilayer. In this Article we describe a coarse-grained (CG) model of POPC and DOPC lipid peroxides within the framework of the Martini CG force-field (v2.2), derived from experimental data. We then discuss extensively the predicted structure and the influence of hydration and show how shifting the position of the unsaturated bonds along the chain changes the structure. Finally, we provide electron and neutron scattering length density profiles of the simulated bilayers.

Molecular basis of the anchoring and stabilization of human islet amyloid polypeptide in lipid hydroperoxidized bilayers

The molecular structure of membrane lipids is formed by mono- or polyunsaturations on their aliphatic tails that make them susceptible to oxidation, facilitating the incorporation of hydroperoxide (R-OOH) functional groups. Such groups promote changes in both composition and complexity of the membrane significantly modifying its physicochemical properties. Human Langerhans islets amyloid polypeptide (hIAPP) is the main component of amyloid deposits found in the pancreas of patients with type-2 diabetes (T2D). hIAPP in the presence of membranes with oxidized lipid species accelerates the formation of amyloid fibrils or the formation of intermediate oligomeric structures. However, the molecular bases at the initial stage of the anchoring and stabilization of the hIAPP in a hydroperoxidized membrane are not yet well understood. To shed some light on this matter, in this contribution, three bilayer models were modeled: neutral (POPC), anionic (POPS), and oxidized (POPC_OOH), and full atom Molecular Dynamics (MD) simulations were performed. Our results show that the POPC_OOH bilayer increases the helicity in hIAPP when compared to POPC or POPS bilayer. The modification in the secondary structure covers the residues of the so-called amyloidogenic core of the hIAPP. Overall, the hydroperoxidation of the neutral lipids modifies both the anchoring and the stabilization of the peptide hIAPP by reducing the random conformations of the peptide and increasing of hydrogen bond population with the hydroperoxidized lipids.

Stiffening Effect of Ceramide on Lipid Membranes Provides Non-Sacrificial Protection against Potent Chemical Damage

Ceramide is a sphingolipid that constitutes only a small fraction of membrane biomolecules but plays a central role in regulating many biological processes. The ceramide level in cell membranes can drastically increase in response to external damage, which has been hypothesized to involve ceramide's biophysical role that increases the membrane packing density and lowers the membrane permeability. However, direct observation of the consequent influence on membrane chemistry resulting from these ceramide-induced physical properties has been absent. Using our unique field-induced droplet ionization mass spectrometry technique combined with molecular dynamics simulations, here we report that the addition of ceramide to POPC monolayer membranes at the air-water interface greatly reduces the chemical damage from potent chemicals, ^•OH radicals, and HCl vapor, by stiffening the membrane packing and lowering the permeability. These results shed new light on the underlying chemoprotective role of ceramide in lipid membranes, which might serve as a previously unknown protection mechanism in response to external stimuli that cause cell stress or death.

The pore-forming activity of sticholysin I is enhanced by the presence of a phospholipid hydroperoxide in membrane

Sticholysin I (StI) is a pore-forming toxin (PFT) belonging to the actinoporin protein family characterized by high permeabilizing activity in membranes. StI readily associates with sphingomyelin (SM)-containing membranes originating pores that can lead to cell death. Binding and pore-formation are critically dependent on the physicochemical properties of membrane. 1-palmitoyl-2-oleoylphosphatidylcholine hydroperoxide (POPC-OOH) is an oxidized phospholipid (OxPL) containing an -OOH moiety in the unsaturated hydrocarbon chain which orientates towards the bilayer interface. This orientation causes an increase in the lipid molecular area, lateral expansion and decrease in bilayer thickness, elastic and bending modulus, as well as modification of lipid packing. Taking advantage of membrane structural changes promoted by POPC-OOH, we investigated its influence on the permeabilizing ability of StI. Here we report the action of StI on Giant Unilamellar Vesicles (GUVs) made of 1-palmitoyl-2-oleoylphosphatidylcholine (POPC) and SM containing increasing amount of POPC-OOH to assess vesicle permeability changes when compared to OxPL-lacking membranes. Inclusion of POPC-OOH in membranes did not promote spontaneous vesicle leaking but resulted in increased membrane permeability due to StI action. StI activity did not modify the fluid-gel phase coexistence boundaries neither in POPC:SM or POPC-OOH:SM membranes. However, the StI insertion mechanism in membrane seems to differ between POPC:SM and POPC-OOH:SM mixtures as suggested by changes in the time course of monolayer surface tension measurements, even though a preferable binding of the toxin to OxPL-containing systems could not be here demonstrated. In summary, modifications in the membrane imposed by lipid hydroperoxidation favor StI permeabilizing activity.

Lipid Hydroperoxide Compromises the Membrane Structure Organization and Softens Bending Rigidity

Lipid hydroperoxides are key mediators of diseases and cell death. In this work, the structural and dynamic perturbations induced by the hydroperoxidized POPC lipid (POPC-OOH) in fluid POPC membranes, at both 23 and 37 °C, were addressed using advanced small-angle X-ray scattering (SAXS) and fluorescence methodologies. Notably, SAXS reveals that the hydroperoxide group decreases the lipid bilayer bending rigidity. This alteration disfavors the bilayer stacking and increases the swelling in-between stacked bilayers. We further investigated the changes in the apolar/polar interface of hydroperoxide-containing membranes through time-resolved fluorescence/anisotropy experiments of the probe TMA-DPH and time-dependent fluorescence shifts of Laurdan. A shorter mean fluorescence lifetime for TMA-DPH was obtained in enriched POPC-OOH membranes, revealing a higher degree of hydration near the membrane interface. Moreover, a higher microviscosity near TMA-DPH and lower order are predicted for these oxidized membranes, at variance with the usual trend of variation of these two parameters. Finally, the complex relaxation process of Laurdan in pure POPC-OOH membranes also indicates a higher membrane hydration and viscosity in the close vicinity of the -OOH moiety. Altogether, our combined approach reveals that the hydroperoxide group promotes alterations in the membrane structure organization, namely, at the level of membrane order, viscosity, and bending rigidity.

Lipid Metabolism Regulates Oxidative Stress and Ferroptosis in RAS-Driven Cancers: A Perspective on Cancer Progression and Therapy

HRAS, NRAS and KRAS, collectively referred to as oncogenic RAS, are the most frequently mutated driver proto-oncogenes in cancer. Oncogenic RAS aberrantly rewires metabolic pathways promoting the generation of intracellular reactive oxygen species (ROS). In particular, lipids have gained increasing attention serving critical biological roles as building blocks for cellular membranes, moieties for post-translational protein modifications, signaling molecules and substrates for ß-oxidation. However, thus far, the understanding of lipid metabolism in cancer has been hampered by the lack of sensitive analytical platforms able to identify and quantify such complex molecules and to assess their metabolic flux in vitro and, even more so, in primary tumors. Similarly, the role of ROS in RAS-driven cancer cells has remained elusive. On the one hand, ROS are beneficial to the development and progression of precancerous lesions, by upregulating survival and growth factor signaling, on the other, they promote accumulation of oxidative by-products that decrease the threshold of cancer cells to undergo ferroptosis. Here, we overview the recent advances in the study of the relation between RAS and lipid metabolism, in the context of different cancer types. In particular, we will focus our attention on how lipids and oxidative stress can either promote or sensitize to ferroptosis RAS driven cancers. Finally, we will explore whether this fine balance could be modulated for therapeutic gain.

Does membrane curvature elastic energy play a role in mediating oxidative stress in lipid membranes?

The effects of oxidative stress on cells are associated with a wide range of pathologies. Oxidative stress is predominantly initiated by the action of reactive oxygen species and/or lipoxygenases on polyunsaturated fatty acid containing lipids. The downstream products are oxidised phospholipids, bioactive aldehydes and a range of Schiff base by-products between aldehydes and lipids, or other biomacromolecules. In this review we assess the impact of oxidative stress on lipid membranes, focusing on the changes that occur to the curvature preference (lipid spontaneous curvature) and elastic properties of membranes, since these biophysical properties modulate phospholipid homeostasis. Studies show that the lipid products of oxidative stress reduce stored curvature elastic energy in membranes. Based upon this observation, we hypothesize that the effects of oxidative stress on lipid membranes will be reduced by compounds that increase stored curvature elastic energy. We find a strong correlation appears across literature studies that we have reviewed, such that many compounds like vitamin E, Curcumin, Coenzyme Q10 and vitamin A show behaviour consistent with this hypothesis. Finally, we consider whether age-related changes in lipid composition represent the homeostatic response of cells to compensate for the accumulation of in vivo lipid oxidation products.

Molecular organization in hydroperoxidized POPC bilayers

Lipid hydroperoxides are the primary reaction products of lipid oxidation, a natural outcome of life under oxygen. While playing a major role in cell metabolism, the microscopic origins of the effects of lipid hydroperoxidation on biomembranes remain elusive. Here we probe the polar structure of partially to fully hydroperoxidized bilayers of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) by a combination of environment-sensitive fluorescent probes and coarse-grained Martini numerical simulations. We find that the inserted organic hydroperoxide group -OOH migrates preferentially to the surface for bilayers with small fractions of hydroperoxidized lipids, but populates also significantly the bilayer interior for larger fractions. Our findings suggest that by modifying the intimate polarity of biomembranes, lipid peroxidation will have a significant impact on the activity of transmembrane proteins and on the bio-medical efficiency of membrane active molecules such as cell-penetrating and antimicrobial peptides.

Effect of magnesium sulfate in oxidized lipid bilayers properties by using molecular dynamics

Magnesium sulfate (MgSO₄) has been used as a protector agent for many diseases related to oxidative stress. The effect of MgSO₄ on the oxidized lipid bilayer has not yet been studied using molecular dynamics calculations. In this work, the effects of oxidation were evaluated by using a POPC membrane model at different concentrations of its aldehyde (-CHO) and hydroperoxide (-OOH) derivatives with and without MgSO₄. Several quantitative and qualitative properties were evaluated, such as membrane thickness, area per lipid, area compressibility modulus, snapshots after simulation finish, density distributions, time evolutions of oxidized group positions, and radial distributions of oxidized group concerning Mg. Results indicate that in the absence of MgSO₄ the mobility of oxidized groups, particularly –CHO, toward the surface interface is high. At a low oxidation level of the bilayer there is an increase in the compressibility modulus as compared to the unoxidized bilayer. MgSO₄, at a low oxidation level, tends to lessen the oxidation effects by lowering the dispersion in the distribution of oxidized species toward the membrane surface and the water region. However, MgSO₄ does not change the trends of decreasing membrane thickness and area compressibility modulus and increasing area per lipid upon oxidation. In this regard, MgSO₄ diminishes the electrostatic long-distance attractive interactions between the oxidized groups and the charged headgroups of the interface, owing to the Mg⁺² and SO₄^-2 screening effects and an electrostatic stabilization of the headgroups, preventing the pore formation, which is well-known to occur in oxidized membranes.

•

MgSO₄ in vitro restores oxidized membranes but its molecular mechanism is unclear.

•

MD simulations of oxidized lipid bilayers were performed with and without of MgSO₄.

•

A restriction in the mobility of oxidized groups is produced by MgSO₄.

•

Mg⁺² and SO₄⁼ produce screening effects on the oxidized membranes.

•

MgSO₄ produce a diminution of electrostatic long-distance attractive interactions.

Site-Specific Peroxidation Modulates Lipid Bilayer Mechanics

Peroxidation of plasma membranes, characterized by oxidative attack of lipidic carbon-carbon double-bonds in unsaturated fatty acids, has been identified as an important biochemical event in multiple pathological conditions, including neurodegenerative diseases, atherosclerosis, diabetes, preeclampsia, aging, cancer, etc. Changes to the lipid bilayer structure as a result of lipid peroxidation may lead to lipid membrane malfunction, and consequently initiate further downstream biochemical cascades. However, how lipid peroxidation modulates the mechanical properties of lipid membranes remains largely controversial. In this study, we investigate the peroxidation of lipids with polyunsaturated fatty acid tails using molecular dynamics simulations. By systematically varying the oxidation site, we find that lipid peroxidation alters the biophysical properties of bilayer membrane in a peroxidation site-specific manner. Specifically, our results suggest that peroxidation at sites in the bilayer interior disturbs and softens the membrane, whereas peroxidation at sites near the membrane-water interface results in a more ordered and stiffer membrane. Such a peroxidation site-specific modulation of lipid membrane mechanics provides an explanation for the contradictory results obtained in previous experiments. Our study paves the way for an improved understanding of the initiation of the downstream cellular dysfunction caused by lipid peroxidation.

Cellular compartments challenged by membrane photo-oxidation

The lipid composition impacts directly on the structure and function of the cytoplasmic as well as organelle membranes. Depending on the type of membrane, specific lipids are required to accommodate, intercalate, or pack membrane proteins to the proper functioning of the cells/organelles. Rather than being only a physical barrier that separates the inner from the outer spaces, membranes are responsible for many biochemical events such as cell-to-cell communication, protein-lipid interaction, intracellular signaling, and energy storage. Photochemical reactions occur naturally in many biological membranes and are responsible for diverse processes such as photosynthesis and vision/phototaxis. However, excessive exposure to light in the presence of absorbing molecules produces excited states and other oxidant species that may cause cell aging/death, mutations and innumerable diseases including cancer. At the same time, targeting key compartments of diseased cells with light can be a promising strategy to treat many diseases in a clinical procedure called Photodynamic Therapy. Here we analyze the relationships between membrane alterations induced by photo-oxidation and the biochemical responses in mammalian cells. We specifically address the impact of photosensitization reactions in membranes of different organelles such as mitochondria, lysosome, endoplasmic reticulum, and plasma membrane, and the subsequent responses of eukaryotic cells.

Lipid Hydroperoxidation Effect on the Dynamical Evolution of the Conductance Process in Bilayer Lipid Membranes: A Condition Toward Criticality

Cell membranes are one of the main targets of oxidative processes mediated by reactive oxygen species (ROS). These chemical species interact with unsaturated fatty acids of membrane lipids, triggering an autocatalytic chain reaction, producing lipid hydroperoxides (LOOHs) as the first relatively stable product of the ROS-mediated lipid peroxidation (LPO) process. Numerous biophysical and computational studies have been carried out to elucidate the LPO impact on the structure and organization of lipid membranes. However, although LOOHs are the major primary product of LPO of polyunsaturated fatty acids (PUFAs), to the best of our knowledge, there is no experimental evidence on the effects of the accumulation of these LPO byproducts on the electrical properties and the underlying dynamics of lipid membranes. In this work, bilayer lipid membranes (BLMs) containing 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocoline (POPC) with increasing hydroperoxidized POPC (POPC-OOH) molar proportions (BLM_PC/PC-OOH) are used as model membranes to investigate the effect of LOOH-mediated LPO propagation on the electrical behavior of lipid bilayers. Voltage-induced ion current signals are analyzed by applying the fractal method of power spectrum density (PSD) analysis. We experimentally prove that, when certain LOOH concentration and energy threshold are overcome, oxidized membranes evolve toward a critical state characterized by the emergence of non-linear electrical behavior dynamics and the pore-type metastable structures formation. PSD analysis shows that temporal dynamics exhibiting "white" noise (non-time correlations) reflects a linear relationship between the input and output signals, while long-term correlations (β > 0.5) begin to be observed closely to the transition (critical point) from linear (Ohmic) to nonlinear (non-Ohmic) behavior. The generation of lipid pores appears to arise as an optimized energy dissipation mechanism based on the system's ability to self-organize and generate ordered structures capable of dissipating energy gradients more efficiently under stressful oxidative conditions.

Effects of Low-level Lipid Peroxidation on the Permeability of Nitroaromatic Molecules across a Membrane: A Computational Study

Lipid peroxidation (LPO) in cellular membranes can cause severe membrane damage and potential cell death. Although oxidized phospholipids have been proved to lead to great changes in the structures and properties of membranes, effects of low-level LPO on membrane permeability have not yet been fully understood. Here, we explored the molecular mechanism of low-level LPO changing the permeability of nitroaromatic molecules across a lipid bilayer by all-atom molecular dynamics simulations. The results reveal that the enhanced passive transport of nitroaromatic molecules lies in the size of defects (i.e., water "finger" and "cone"), which is further dependent on the extent of LPO and the structural feature of solutes. In detail, if the solute can form more hydrogen bonds with water, which stabilizes the water into a large-size cone, there is a greater permeability coefficient (P). Otherwise, a small-size finger only results in a small increase of P. For example, the presence of 15% oxidized lipids could result in an increase of 2,4,6-trinitrotoluene (TNT's) P by more than 2 orders of magnitude (from 1.7 × 10^-2 to 2.39 cm·s^-1). The result suggests that the membrane permeability can be greatly promoted in the physiologically relevant environment with low-level LPO, and more importantly, clarifies the contributions of both the hydrophobicity of the membrane interior and the structural feature of solutes to such enhanced permeability. This work may provide significant insight into the toxic effects of nitroaromatic molecules and the pharmaceutical characteristics of tissues with oxidative damage.

Self-Assembly of Phosphocholine Derivatives Using the ELBA Coarse-Grained Model: Micelles, Bicelles, and Reverse Micelles

The ELBA coarse-grained force field was originally developed for lipids, and its water model is described as a single-site Lennard-Jones particle with electrostatics modeled by an embedded point-dipole, while other molecules in this force field have a three (or four)-to-one mapping scheme. Here, ELBA was applied to investigate the self-assembly processes of dodecyl-phosphocholine (DPC) micelle, 1,2-dipalmitoyl-sn-glycero-3-phosphocholine/1,2-dihexaoyl-sn-glycero-3-phosphocholine (DPPC/DHPC) bicelles, and DPPC/cyclohexane/water reverse micelles through coarse-grained molecular dynamics (MD) simulations. New parameters were obtained using a simplex algorithm-based calibration procedure to determine the Lennard-Jones parameters for cyclohexane, dodecane, and cyclohexane-dodecane cross-interactions. Density, self-diffusion coefficient, surface tension, and mixture excess volume were found to be in fair agreement with experimental data. These new parameters were used in the simulations, and the obtained structures were analyzed for shape, size, volume, and surface area. Except for the shape of DPC micelles, all other properties match well with available experimental data and all-atom simulations. Remarkably, in agreement with experiments the rodlike shape of the DPPC reverse micelle is well described by ELBA, while all-atom data in the literature predicts a disclike shape. To further check the consistency of the force field in reproducing the correct shapes of reverse micelles, additional simulations were performed doubling the system size. Two distinct reverse micelles were obtained both presenting the rodlike shape and correct aggregation number.

Atomic-scale modeling of the effect of lipid peroxidation on the permeability of reactive species

Biomembranes and lipid systems are rich in unsaturated lipid components and are subject to photo-induced lipid peroxidation. The peroxidized lipid products in cellular systems are known to affect the structural organization and function of the biomembrane. We employed molecular dynamics simulations to study the effects of phospholipid peroxidation on membrane properties and the permeability of different reactive species. The results suggest that when the lipids are peroxidized, the peroxide group moves toward the membrane surface, which causes the membrane system to expand laterally and increase in area. The permeability profile revealed that nitrogen species can easily permeate through the native and peroxidized system in comparison to oxygen species, suggesting its importance in plasma-based treatment. Thus, by breaching the energy barrier with lower energy, they can traverse the cell membrane and induce oxidative stress, which leads to apoptosis.Communicated by Ramaswamy H. Sarma.

Understanding Interactions of Curcumin with Lipid Bilayers: A Coarse-Grained Molecular Dynamics Study

The interactions of curcumin with various lipid bilayers (POPC, DOPC, oxidized POPC, and oxidized DOPC) and model biomembranes (symmetric bacteria and yeast plasma membranes, as well as asymmetric mammalian plasma membrane) are investigated. A nonlinear thinning effect of curcumin with respect to its concentration is demonstrated in PC membranes and in the yeast. Curcumin induces asymmetry to the symmetric yeast membranes but reduces the degree of asymmetry of the mammalian plasma membranes when the molecule is placed facing the outer leaflets. The molecule is found to diffuse through oxidized PC bilayers, POPC bilayers at a curcumin to lipid ratio C/L = 1/5, yeast membranes at C/L = 1/100, and the mammalian plasma membranes at C/L = 1/5 and when the molecule placed facing the outer leaflets. The results of this work demonstrate that the lipid type, the lipid distribution, and curcumin amount play a critical role in defining the interactions of curcumin with the lipids and their transport behavior through the bilayers.

Understanding the Role of Lipids in Signaling Through Atomistic and Multiscale Simulations of Cell Membranes

Cell signaling controls essentially all cellular processes. While it is often assumed that proteins are the key architects coordinating cell signaling, recent studies have shown more and more clearly that lipids are also involved in signaling processes in a number of ways. Lipids do, for instance, act as messengers, modulate membrane receptor conformation and dynamics, and control membrane receptor partitioning. Further, through structural modifications such as oxidation, the functions of lipids as part of signaling processes can be modified. In this context, in this article we discuss the understanding recently revealed by atomistic and coarse-grained computer simulations of nanoscale processes and underlying physicochemical principles related to lipids’ functions in cellular signaling.

Solidification to improve the biopharmaceutical performance of SEDDS: Opportunities and challenges

Self-emulsifying drug delivery systems (SEDDS) offer potential for overcoming the inherent slow dissolution and poor oral absorption of hydrophobic drugs by retaining them in a solubilised state during gastrointestinal transit. However, the promising biopharmaceutical benefits of liquid lipid formulations has not translated into widespread commercial success, due to their susceptibility to long term storage and in vivo precipitation issues. One strategy that has emerged to overcome such limitations, is to combine the solubilisation and dissolution enhancing properties of lipids with the stabilising effects of solid carrier materials. The development of intelligent hybrid drug formulations has presented new opportunities to harness the potential of emulsified lipids in optimising oral bioavailability for lipophilic therapeutics. Specific emphasis of this review is placed on the impact of solidification approaches and excipients on the biopharmaceutical performance of self-emulsifying lipids, with findings highlighting the key design considerations that should be implemented when developing hybrid lipid-based formulations.

Physical properties of model biological lipid bilayers: insights from all-atom molecular dynamics simulations

The physical properties of lipid bilayers are sensitive to the specific type and composition of the lipids that make up the many different types of cell membranes. Studying model bilayers of representative heterogeneous compositions can provide key insights into membrane functionality. In this work, we use atomistic molecular dynamics simulations to characterize key properties in a number of bilayer membranes of varying composition. We first examine basic properties, such as lipid area, volume, and bilayer thickness, of simple, homogeneous bilayers comprising several lipid types, which are prevalent in biological membranes. Such lipids are then used in simulations of heterogeneous systems representative of bacterial, mammalian, and cancer membranes. Our analysis is especially focused on depth-dependent, transmembrane profiles; in particular, we calculate lateral pressure and dipole potential profiles, two fundamental properties which play key roles in a large number of biological functions.

Contrasting roles of oxidized lipids in modulating membrane microdomains

Lipid rafts display a lateral heterogeneity forming membrane microdomains that hold a fundamental role on biological membranes and are indispensable to physiological functions of cells. Oxidative stress in cellular environments may cause lipid oxidation, changing membrane composition and organization, thus implying in effects in cell signaling and even loss of homeostasis. The individual contribution of oxidized lipid species to the formation or disruption of lipid rafts in membranes still remains unknown. Here, we investigate the role of different structures of oxidized phospholipids on rafts microdomains by carefully controlling the membrane composition. Our experimental approach based on fluorescence microscopy of giant unilamellar vesicles (GUV) enables the direct visualization of the impact of hydroperoxidized POPC lipid (referred to as POPCOOH) and shortened chain lipid PazePC (1-palmitoyl-2-azelaoyl-sn-glycero-3-phosphocholine) on phase separation. We found that the molecular structure of oxidized lipid is of paramount importance on lipid mixing and/or demixing. The hydrophobic mismatch promoted by POPCOOH coupled to its cylindrical molecular shape favor microdomains formation. In contrast, the conical shape of PazePC causes disarrangement of lipid 2D organized platforms. Our findings contribute to better unraveling how oxidized phospholipids can trigger formation or disruption of lipid rafts. As a consequence, phospholipid oxidation may indirectly affect association or dissociation of key biomolecules in the rafts thus altering cell signaling and homeostasis.

Modeling the effects of lipid peroxidation during ferroptosis on membrane properties

Ferroptosis is a form of regulated cell death characterized by the accumulation of lipid hydroperoxides. There has been significant research on the pathways leading to the accumulation of oxidized lipids, but the downstream effects and how lipid peroxides cause cell death during ferroptosis remain a major puzzle. We evaluated key features of ferroptosis in newly developed molecular dynamics models of lipid membranes to investigate the biophysical consequences of lipid peroxidation, and generated hypotheses about how lipid peroxides contribute to cell death during ferroptosis.

Alpha-tocopherol inhibits pore formation in oxidized bilayers

In biological membranes, alpha-tocopherols (α-toc; vitamin E) protect polyunsaturated lipids from free radicals. Although the interactions of α-toc with non-oxidized lipid bilayers have been studied, their effects on oxidized bilayers remain unknown. In this study, atomistic molecular dynamics (MD) simulations of oxidized lipid bilayers were performed with varying concentrations of α-toc. Bilayers with 1-palmitoyl-2-lauroyl-sn-glycero-3-phosphocholine (PLPC) lipids and their aldehyde derivatives at a 1 : 1 ratio were studied. Our simulations show that oxidized lipids self-assemble into aggregates with a water pore rapidly developing across the bilayer. The free energy of transporting an α-toc molecule in a bilayer suggests that α-tocs can passively adsorb into it. When α-toc molecules were present at low concentrations in bilayers containing oxidized lipids, water pore formation was slowed down. At high α-toc concentrations, no pores were observed. Based on the simulations, we propose that the mechanism of how α-toc inhibits pore formation in bilayers with oxidized lipids is the following: α-tocs trap the polar groups of the oxidized lipids at the membrane-water interface resulting in a decreased probability of the oxidized lipids making contact with the two leaflets and initiating pore formation. This demonstrates that α-toc molecules not only protect the bilayer from oxidation but also help to stabilize the bilayer after lipid peroxidation occurs. These results will help in designing more efficient molecules to protect membranes from oxidative stress.