Influence of ester-modified lipids on bilayer structure

Salt effects on mixed composition membranes containing an antioxidant lipophilic edaravone derivative: a computational-experimental study

The protection of lipid membranes against oxidation avoids diseases associated with oxidative stress. As a strategy to contrast it, functionalized lipids with antioxidant activity are used to become part of membranes thus protecting them. For this purpose, a lipophilic edaravone derivative has been synthesized, adding a C18 saturated chain to the original structure. The antioxidant activity of C18-Edv has been demonstrated in our previous work. In this study, molecular dynamics simulations have been performed to define the effects of NaCl, MgCl₂, KCl, and CaCl₂ salts on a palmitoyl-oleoyl-sn-glycero-phosphocholine (POPC) lipid bilayer encapsulating C18-Edv. The results showed how different salts influence POPC lateral diffusion, and the movements of C18-Edv heads, which are antioxidant moieties, were correlated to the ability of C18-Edv molecules to protect membranes. MgCl₂ showed a negative impact leading to C18-Edv clusterization and membrane stretching, while KCl and NaCl showed a moderate influence on the mixed lipid membrane structure. CaCl₂ increased the exposure of the C18-Edv heads to the lipid-water interface, resulting in the salt with a higher propensity to guarantee protection against radicals in the aqueous phase. Finally, C18-Edv-POPC liposomes have been prepared following the simulation conditions, and then an experimental Oxygen Radical Absorbance Capacity (ORAC) assay has been performed to confirm the in silico predicted results.

Influence of a lipophilic edaravone on physical state and activity of antioxidant liposomes: An experimental and in silico study

The influence of a lipophilic derivative of Edaravone (C18Edv) on a POPC liposomal bilayer has been investigated by a combined computational-experimental approach. The order and hydration degree of three different systems composed by 10%, 20% and 40% in w/w percentage of C18Edv respect to POPC were investigated through Molecular Dynamics (MD) simulations and fluorescence spectroscopy experiments. Dynamic Light Scattering measurements showed how the presence of different amounts of C18EdV determines differences on liposome size and stability. The survey revealed that the content of lipophilic antioxidant tunes liposome rigidity and influences cellular uptake and antioxidant activity which are maximized for formulation containing 20% of C18Edv.

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.

Liposomal Formulations for an Efficient Encapsulation of Epigallocatechin-3-gallate: An in- Silico/Experimental Approach

As a part of research project aimed to optimize antioxidant delivery, here we studied the influence of both salts and lipid matrix composition on the interaction of epigallocatechin-3-gallate (EGCG) with bilayer leaflets. Thus, we combined in silico and experimental methods to study the ability of neutral and anionic vesicles to encapsulate EGCG in the presence of Ca²⁺ and Mg²⁺ divalent salts. Experimental and in silico results show a very high correlation, thus confirming the efficiency of the developed methodology. In particular, we found out that the presence of calcium ions hinders the insertion of EGCG in the liposome bilayer in both neutral and anionic systems. On the contrary, the presence of MgCl₂ improves the insertion degree of EGCG molecules respect to the liposomes without divalent salts. The best and most efficient salt concentration is that corresponding to a 5:1 molar ratio between Mg²⁺ and EGCG, in both neutral and anionic vesicles. Concerning the lipid matrix composition, the anionic one results in better promotion of the catechin insertion within the bilayer since experimentally we achieved 100% EGCG encapsulation in the lipid carrier in the presence of a 5:1 molar ratio of magnesium. Thus, the combination of this anionic liposomal formulation with magnesium chloride, avoids time-consuming separation steps of unentrapped active principle and appears particularly suitable for EGCG delivery applications.

Salts Influence Cathechins and Flavonoids Encapsulation in Liposomes: A Molecular Dynamics Investigation

Cathechins and flavonoids are responsible of numerous health benefits. Two of the most representatives' compounds for their antioxidant and therapeutic effects are Epigallocatechin 3-Gallate (EGCG), from green tea extracts, and morelloflavone (MF), from Garcinia dulcis. Here we explore, by atomistic Molecular Dynamics simulations, how EGCG and MF interact with lipid bilayers and we show the salts' influence on their encapsulation degree in neutral liposomes. As a result, we found out that EGCGs naturally bind to the hydrophilic regions of phospholipids, positioning themselves mostly at the interface between water and lipid phases. The presence of a salt clearly influences the EGCG molecules' absorption and the total effect depends strongly on the salt nature and concentration. Beside, for MF, we observed a high stability of the intermolecular MFs aggregates in water that strongly penalizes the flavonoid's interaction with the lipid polar heads. However, salts can influence MF's liposomal penetration, even if they are not able to promote completely its absorption inside the bilayer. For both compounds, the increase of penetration is more marked in presence of magnesium chloride, whilst calcium chloride showed the opposite effect.

A molecular insight into the electro-transfer of small molecules through electropores driven by electric fields

The transport of chemical compounds across the plasma membrane into the cell is relevant for several biological and medical applications. One of the most efficient techniques to enhance this uptake is reversible electroporation. Nevertheless, the detailed molecular mechanism of transport of chemical species (dyes, drugs, genetic materials, …) following the application of electric pulses is not yet fully elucidated. In the past decade, molecular dynamics (MD) simulations have been conducted to model the effect of pulsed electric fields on membranes, describing several aspects of this phenomenon. Here, we first present a comprehensive review of the results obtained so far modeling the electroporation of lipid membranes, then we extend these findings to study the electrotransfer across lipid bilayers subject to microsecond pulsed electric fields of Tat11, a small hydrophilic charged peptide, and of siRNA. We use in particular a MD simulation protocol that allows to characterize the transport of charged species through stable pores. Unexpectedly, our results show that for an electroporated bilayer subject to transmembrane voltages in the order of 500mV, i.e. consistent with experimental conditions, both Tat11 and siRNA can translocate through nanoelectropores within tens of ns. We discuss these results in comparison to experiments in order to rationalize the mechanism of drug uptake by cells. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

CHARMM-GUI Input Generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM Simulations Using the CHARMM36 Additive Force Field

Proper treatment of nonbonded interactions is essential for the accuracy of molecular dynamics (MD) simulations, especially in studies of lipid bilayers. The use of the CHARMM36 force field (C36 FF) in different MD simulation programs can result in disagreements with published simulations performed with CHARMM due to differences in the protocols used to treat the long-range and 1-4 nonbonded interactions. In this study, we systematically test the use of the C36 lipid FF in NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM. A wide range of Lennard-Jones (LJ) cutoff schemes and integrator algorithms were tested to find the optimal simulation protocol to best match bilayer properties of six lipids with varying acyl chain saturation and head groups. MD simulations of a 1,2-dipalmitoyl-sn-phosphatidylcholine (DPPC) bilayer were used to obtain the optimal protocol for each program. MD simulations with all programs were found to reasonably match the DPPC bilayer properties (surface area per lipid, chain order parameters, and area compressibility modulus) obtained using the standard protocol used in CHARMM as well as from experiments. The optimal simulation protocol was then applied to the other five lipid simulations and resulted in excellent agreement between results from most simulation programs as well as with experimental data. AMBER compared least favorably with the expected membrane properties, which appears to be due to its use of the hard-truncation in the LJ potential versus a force-based switching function used to smooth the LJ potential as it approaches the cutoff distance. The optimal simulation protocol for each program has been implemented in CHARMM-GUI. This protocol is expected to be applicable to the remainder of the additive C36 FF including the proteins, nucleic acids, carbohydrates, and small molecules.

Membrane pore formation in atomistic and coarse-grained simulations

Biological cells and their organelles are protected by ultra thin membranes. These membranes accomplish a broad variety of important tasks like separating the cell content from the outer environment, they are the site for cell-cell interactions and many enzymatic reactions, and control the in- and efflux of metabolites. For certain physiological functions e.g. in the fusion of membranes and also in a number of biotechnological applications like gene transfection the membrane integrity needs to be compromised to allow for instance for the exchange of polar molecules across the membrane barrier. Mechanisms enabling the transport of molecules across the membrane involve membrane proteins that form specific pores or act as transporters, but also so-called lipid pores induced by external fields, stress, or peptides. Recent progress in the simulation field enabled to closely mimic pore formation as supposed to occur in vivo or in vitro. Here, we review different simulation-based approaches in the study of membrane pores with a focus on lipid pore properties such as their size and energetics, poration mechanisms based on the application of external fields, charge imbalances, or surface tension, and on pores that are induced by small molecules, peptides, and lipids. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Bilayer Deformation, Pores, and Micellation Induced by Oxidized Lipids

The influence of different oxidized lipids on lipid bilayers was investigated with 16 individual 1 μs atomistic molecular dynamics (MD) simulations. Binary mixtures of lipid bilayers of 1-palmitoyl-2-linoleoyl-sn-glycero-3-phosphatidylcholine (PLPC) and its peroxide and aldehyde products were performed at different concentrations. In addition, an asymmetrical short chain lipid, 1-palmitoyl-2-decanoyl-sn-glycero-3-phosphatidylcholine (PDPC), was used to compare the effects of polar/apolar groups in the lipid tail on lipid bilayer. Although water defects occurred with both aldehyde and peroxide lipids, full pore formation was observed only for aldehyde lipids. At medium concentrations the pores were stable. At higher concentrations, however, the pores became unstable and micellation occurred. Data analysis shows that aldehyde lipids' propensity for pore formation is due to their shorter and highly mobile tail. The highly polar peroxide lipids are stabilized by strong hydrogen bonds with interfacial water.

CHARMM-GUI Input Generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM Simulations Using the CHARMM36 Additive Force Field

Proper treatment of nonbonded interactions is essential for the accuracy of molecular dynamics (MD) simulations, especially in studies of lipid bilayers. The use of the CHARMM36 force field (C36 FF) in different MD simulation programs can result in disagreements with published simulations performed with CHARMM due to differences in the protocols used to treat the long-range and 1-4 nonbonded interactions. In this study, we systematically test the use of the C36 lipid FF in NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM. A wide range of Lennard-Jones (LJ) cutoff schemes and integrator algorithms were tested to find the optimal simulation protocol to best match bilayer properties of six lipids with varying acyl chain saturation and head groups. MD simulations of a 1,2-dipalmitoyl-sn-phosphatidylcholine (DPPC) bilayer were used to obtain the optimal protocol for each program. MD simulations with all programs were found to reasonably match the DPPC bilayer properties (surface area per lipid, chain order parameters, and area compressibility modulus) obtained using the standard protocol used in CHARMM as well as from experiments. The optimal simulation protocol was then applied to the other five lipid simulations and resulted in excellent agreement between results from most simulation programs as well as with experimental data. AMBER compared least favorably with the expected membrane properties, which appears to be due to its use of the hard-truncation in the LJ potential versus a force-based switching function used to smooth the LJ potential as it approaches the cutoff distance. The optimal simulation protocol for each program has been implemented in CHARMM-GUI. This protocol is expected to be applicable to the remainder of the additive C36 FF including the proteins, nucleic acids, carbohydrates, and small molecules.

Mechanical properties of lipid bilayers from molecular dynamics simulation

Lipid areas (Aℓ), bilayer area compressibilities (KA), bilayer bending constants (KC), and monolayer spontaneous curvatures (c0) from simulations using the CHARMM36 force field are reported for 12 representative homogenous lipid bilayers. Aℓ (or their surrogate, the average deuterium order parameter in the "plateau region" of the chain) agree very well with experiment, as do the KA. Simulated KC are in near quantitative agreement with vesicle flicker experiments, but are somewhat larger than KC from X-ray, pipette aspiration, and neutron spin echo for saturated lipids. Spontaneous curvatures of bilayer leaflets from the simulations are approximately 30% smaller than experimental values of monolayers in the inverse hexagonal phase.

Liposome-based gene delivery systems containing a steroid derivative: computational and small angle X-ray diffraction study

The structural properties and the phase behaviour of mixed composition liposomes containing a functionalized lipid are investigated with the aim to design neutral liposomes able to coordinate metals and to complex DNA.