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.

Distribution of lipid aldehydes in phase-separated membranes: A molecular dynamics study.

It is well established that lipid aldehydes (LAs) are able to increase the permeability of cell membranes and induce their rupture. However, it is not yet clear how LAs are distributed in phase-separated membranes (PSMs), which are responsible for the transport of selected molecules and intracellular signaling. Thus, we investigate here the distribution of LAs in a PSM by coarse-grained molecular dynamics simulations. Our results reveal that LAs derived from mono-unsaturated lipids tend to accumulate at the interface between the liquid-ordered/liquid-disordered domains, whereas those derived from poly-unsaturated lipids remain in the liquid-disordered domain. These results are important for understanding the effects caused by oxidized lipids in membrane structure, properties and organization.

Lipid Oxidation: Role of Membrane Phase-Separated Domains.

Lipid oxidation is associated with several inflammatory and neurodegenerative diseases, but many questions to unravel its effects on biomembranes are still open due to the complexity of the topic. For instance, recent studies indicated that phase-separated domains can have a significant effect on membrane function. It is reported that domain interfaces are "hot spots" for pore formation, but the underlying mechanisms and the effect of oxidation-induced phase separation on membranes remain elusive. Thus, to evaluate the permeability of the membrane coexisting of liquid-ordered (Lo) and liquid-disordered (Ld) domains, we performed atomistic molecular dynamics simulations. Specifically, we studied the membrane permeability of nonoxidized or oxidized homogeneous membranes (single-phase) and at the Lo/Ld domain interfaces of heterogeneous membranes, where the Ld domain is composed of either oxidized or nonoxidized lipids. Our simulation results reveal that the addition of only 1.5% of lipid aldehyde molecules at the Lo/Ld domain interfaces of heterogeneous membranes increases the membrane permeability, whereas their addition at homogeneous membranes does not have any effect. This study is of interest for a better understanding of cancer treatment methods based on oxidative stress (causing among others lipid oxidation), such as plasma medicine and photodynamic therapy.

How do nitrated lipids affect the properties of phospholipid membranes?

Biological membranes are under constant attack of free radicals, which may lead to lipid nitro-oxidation, producing a complex mixture of nitro-oxidized lipids that are responsible for structural and dynamic changes on the membrane. Despite the latter, nitro-oxidized lipids are also associated with several inflammatory and neurodegenerative diseases, the underlying mechanisms of which remain elusive. We perform atomistic molecular dynamics simulations using several isomers of nitro-oxidized lipids to study their effect on the structure and permeability of the membrane, as well as the interaction between the mixture of these products in the phospholipid membrane environment. Our results show that the stereo- and positional isomers have a stronger effect on the properties of the membrane composed of oxidized lipids compared to that containing nitrated lipids. Nevertheless, nitrated lipids lead to three-fold increase in water permeability compared to oxidized lipids. In addition, we show that in a membrane consisting of combined nitro-oxidized lipid products, the presence of oxidized lipids protects the membrane from transient pores. Is well stablished that plasma application and photodynamic therapy produces a number of oxidative species used to kill cancer cells, through membrane damage induced by nitro-oxidative stress. This study is important to elucidate the mechanisms and the molecular level properties involving the reactive species produced during that cancer therapies.

Photosensitized Membrane Permeabilization Requires Contact-Dependent Reactions between Photosensitizer and Lipids.

Although the general mechanisms of lipid oxidation are known, the chemical steps through which photosensitizers and light permeabilize lipid membranes are still poorly understood. Herein we characterized the products of lipid photooxidation and their effects on lipid bilayers, also giving insight into their formation pathways. Our experimental system was designed to allow two phenothiazinium-based photosensitizers (methylene blue, MB, and DO15) to deliver the same amount of singlet oxygen molecules per second to 1-palmitoyl-2-oleoyl- sn-glycero-3-phosphocholine liposome membranes, but with a substantial difference in terms of the extent of direct physical contact with lipid double bonds; that is, DO15 has a 27-times higher colocalization with ω-9 lipid double bonds than MB. Under this condition, DO15 permeabilizes membranes at least 1 order of magnitude more efficiently than MB, a result that was also valid for liposomes made of polyunsaturated lipids. Quantification of reaction products uncovered a mixture of phospholipid hydroperoxides, alcohols, ketones, and aldehydes. Although both photosensitizers allowed the formation of hydroperoxides, the oxidized products that require direct reactions between photosensitizer and lipids were more prevalent in liposomes oxidized by DO15. Membrane permeabilization was always connected with the presence of lipid aldehydes, which cause a substantial decrease in the Gibbs free energy barrier for water permeation. Processes depending on direct contact between photosensitizers and lipids were revealed to be essential for the progress of lipid oxidation and consequently for aldehyde formation, providing a molecular-level explanation of why membrane binding correlates so well with the cell-killing efficiency of photosensitizers.

Molecular simulations of the effects of phospholipid and cholesterol peroxidation on lipid membrane properties.

Non-enzymatic lipid peroxidation may change biomembrane structure and function. Here, we employed molecular dynamics simulations to study the effects of either phospholipid or cholesterol peroxidation individually, as well as the combined peroxidation of both components. When lipids were peroxidized, the generated OOH groups migrated to the membrane surface and engaged in H-bonds with each other and the phospholipid carbonyl ester groups. It caused the sn-2 acyl chains of phospholipid hydroperoxides to bend and the whole sterol backbone of cholesterol hydroperoxides to tilt. When phospholipids were kept intact, peroxidation of the sterol backbone led to a partial degradation of its condensing and ordering properties, independently of the position and isomerism of the OOH substitution. However, even in massively peroxidized membranes in which all phospholipids and cholesterol were peroxidized, the condensing and ordering properties of the sterol backbone were still significant. The possible implications for the formation of membrane lateral domains were discussed. Cholesterol peroxyl radicals were also investigated and we found that the OO groups did not migrate to the headgroups region.

Molecular dynamics simulations of the transport of reactive oxygen species by mammalian and plant aquaporins.

BACKGROUND: Aquaporins are responsible for water transport across lipid membranes. They are also able to transport reactive oxygen species, playing an important role in redox signaling. Certain plant aquaporins have even the ability to be regulated by oxidative stress. However, the underlying mechanisms are still not fully understood. METHODS: Here, molecular dynamics simulations were employed to determine the activation free energies related to the transport of reactive oxygen species through both mammalian and plant aquaporin models. RESULTS AND CONCLUSIONS: Both aquaporins may transport hydrogen peroxide (H2O2) and the protonated form of superoxide radicals (HO2). The solution-to-pore transfer free energies were low for small oxy-radicals, suggesting that even highly reactive hydroxyl radicals (HO) might have access to the pore interior and oxidize amino acids responsible for channel selectivity. In the plant aquaporin, no significant change in water permeability was observed upon oxidation of the solvent-exposed disulfide bonds at the extracellular region. During the simulated time scale, the existence of a direct oxidative gating mechanism involving these disulfide bonds could not be demonstrated. GENERAL SIGNIFICANCE: Simulation results may improve the understanding of redox signaling mechanisms and help in the interpretation of protein oxidative labeling experiments.

Reactive oxygen species at phospholipid bilayers: distribution, mobility and permeation.

Reactive oxygen species (ROS) are involved in biochemical processes such as redox signaling, aging, carcinogenesis and neurodegeneration. Although biomembranes are targets for reactive oxygen species attack, little is known about the role of their specific interactions. Here, molecular dynamics simulations were employed to determine the distribution, mobility and residence times of various reactive oxygen species at the membrane-water interface. Simulations showed that molecular oxygen (O2) accumulated at the membrane interior. The applicability of this result to singlet oxygen ((1)O2) was discussed. Conversely, superoxide (O2(-)) radicals and hydrogen peroxide (H2O2) remained at the aqueous phase. Both hydroxyl (HO) and hydroperoxyl (HO2) radicals were able to penetrate deep into the lipid headgroups region. Due to membrane fluidity and disorder, these radicals had access to potential peroxidation sites along the lipid hydrocarbon chains, without having to overcome the permeation free energy barrier. Strikingly, HO2 radicals were an order of magnitude more concentrated in the headgroups region than in water, implying a large shift in the acid-base equilibrium between HO2 and O2(-). In comparison with O2, both HO and HO2 radicals had lower lateral mobility at the membrane. Simulations revealed that there were intermittent interruptions in the H-bond network around the HO radicals at the headgroups region. This effect is expected to be unfavorable for the H-transfer mechanism involved in HO diffusion. The implications for lipid peroxidation and for the effectiveness of membrane antioxidants were evaluated.

The pro- and anti-tumoral properties of gap junctions in cancer and their role in therapeutic strategies.

Gap junctions (GJs), essential structures for cell-cell communication, are made of two hemichannels (commonly called connexons), one on each adjacent cell. Found in almost all cells, GJs play a pivotal role in many physiological and cellular processes, and have even been linked to the progression of diseases, such as cancer. Modulation of GJs is under investigation as a therapeutic strategy to kill tumor cells. Furthermore, GJs have also been studied for their key role in activating anti-cancer immunity and propagating radiation- and oxidative stress-induced cell death to neighboring cells, a process known as the bystander effect. While, gap junction (GJ)-based therapeutic strategies are being developed, one major challenge has been the paradoxical role of GJs in both tumor progression and suppression, based on GJ composition, cancer factors, and tumoral context. Therefore, understanding the mechanisms of action, regulation, and the dual characteristics of GJs in cancer is critical for developing effective therapeutics. In this review, we provide an overview of the current understanding of GJs structure, function, and paradoxical pro- and anti-tumoral role in cancer. We also discuss the treatment strategies to target these GJs properties for anti-cancer responses, via modulation of GJ function.

Unraveling the permeation of reactive species across nitrated membranes by computer simulations.

Reactive oxygen and nitrogen species (RONS) are involved in many biochemical processes, including nitro-oxidative stress that causes cancer cell death, observed in cancer therapies such as photodynamic therapy and cold atmospheric plasma. However, their mechanisms of action and selectivity still remain elusive due to the complexity of biological cells. For example, it is not well known how RONS generated by cancer therapies permeate the cell membrane to cause nitro-oxidative damage. There are many studies dedicated to the permeation of RONS across native and oxidized membranes, but not across nitrated membranes, another lipid product also generated during nitro-oxidative stress. Herein, we performed molecular dynamics (MD) simulations to calculate the free energy barrier of RONS permeation across nitrated membranes. Our results show that hydrophilic RONS, such as hydroperoxyl radical (HO2) and peroxynitrous acid (ONOOH), have relatively low barriers compared to hydrogen peroxide (H2O2) and hydroxyl radical (HO), and are more prone to permeate the membrane than for the native or peroxidized membranes, and similar to aldehyde-oxidized membranes. Hydrophobic RONS like molecular oxygen (O2), nitrogen dioxide (NO2) and nitric oxide (NO) even have insignificant barriers for permeation. Compared to native and peroxidized membranes, nitrated membranes are more permeable, suggesting that we must not only consider oxidized membranes during nitro-oxidative stress, but also nitrated membranes, and their role in cancer therapies.

Oxidative damage to hyaluronan-CD44 interactions as an underlying mechanism of action of oxidative stress-inducing cancer therapy.

Multiple cancer therapies nowadays rely on oxidative stress to damage cancer cells. Here we investigated the biological and molecular effect of oxidative stress on the interaction between CD44 and hyaluronan (HA), as interrupting their binding can hinder cancer progression. Our experiments demonstrated that the oxidation of HA decreased its recognition by CD44, which was further enhanced when both CD44 and HA were oxidized. The reduction of CD44-HA binding negatively affected the proliferative state of cancer cells. Our multi-level atomistic simulations revealed that the binding free energy of HA to CD44 decreased upon oxidation. The effect of HA and CD44 oxidation on CD44-HA binding was similar, but when both HA and CD44 were oxidized, the effect was much larger, in agreement with our experiments. Hence, our experiments and computations support our hypothesis on the role of oxidation in the disturbance of CD44-HA interaction, which can lead to the inhibition of proliferative signaling pathways inside the tumor cell to induce cell death.

Parametrization and Molecular Dynamics Simulations of Nitrogen Oxyanions and Oxyacids for Applications in Atmospheric and Biomolecular Sciences.

Nitrogen oxyanions and oxyacids are important agents in atmospheric chemistry and medical biology. Although their chemical behavior in solution is relatively well understood, they may behave very differently at the water/air interface of atmospheric aerosols or at the membrane/water interface of cells. Here, we developed a fully classical model for molecular dynamics simulations of NO3-, NO2-, HNO3, and HNO2 in the framework of the GROMOS 53A6 and 54A7 force field versions. The model successfully accounted for the poorly structured solvation shell and ion pairing tendency of NO3-. Accurate pure-liquid properties and hydration free energies were obtained for the oxyacids. Simulations at the water/air interface showed a local enrichment of HNO3 and depletion of NO3-. The effect was discussed in light of earlier spectroscopic data and ab initio calculations, suggesting that HNO3 behaves as a weaker acid at the surface of water. Our model will hopefully allow for efficient and accurate simulations of nitrogen oxyanions and oxyacids in solution and at microheterogeneous interface environments.

Transport of Reactive Oxygen and Nitrogen Species across Aquaporin: A Molecular Level Picture.

Aquaporins (AQPs) are transmembrane proteins that conduct not only water molecules across the cell membrane but also other solutes, such as reactive oxygen and nitrogen species (RONS), produced (among others) by cold atmospheric plasma (CAP). These RONS may induce oxidative stress in the cell interior, which plays a role in cancer treatment. The underlying mechanisms of the transport of RONS across AQPs, however, still remain obscure. We apply molecular dynamics simulations to investigate the permeation of both hydrophilic (H2O2 and OH) and hydrophobic (NO2 and NO) RONS through AQP1. Our simulations show that these RONS can all penetrate across the pores of AQP1. The permeation free energy barrier of OH and NO is lower than that of H2O2 and NO2, indicating that these radicals may have easier access to the pore interior and interact with the amino acid residues of AQP1. We also study the effect of RONS-induced oxidation of both the phospholipids and AQP1 (i.e., sulfenylation of Cys191) on the transport of the above-mentioned RONS across AQP1. Both lipid and protein oxidation seem to slightly increase the free energy barrier for H2O2 and NO2 permeation, while for OH and NO, we do not observe a strong effect of oxidation. The simulation results help to gain insight in the underlying mechanisms of the noticeable rise of CAP-induced RONS in cancer cells, thereby improving our understanding on the role of AQPs in the selective anticancer capacity of CAP.

Molecular dynamics simulations of mechanical stress on oxidized membranes.

Biomembranes are under constant attack of free radicals that may lead to lipid oxidation in conditions of oxidative stress. The products generated during lipid oxidation are responsible for structural and dynamical changes which may jeopardize the membrane function. For instance, the local rearrangements of oxidized lipid molecules may induce membrane rupture. In this study, we investigated the effects of mechanical stress on oxidized phospholipid bilayers (PLBs). Model bilayers were stretched until pore formation (or poration) using non-equilibrium molecular dynamics simulations. We studied single-component homogeneous membranes composed of lipid oxidation products, as well as two-component heterogeneous membranes with coexisting native and oxidized domains. In homogeneous membranes, the oxidation products with -OH and -OOH groups reduced the areal strain required for pore formation, whereas the oxidation product with O group behaved similarly to the native membrane. In heterogeneous membranes composed of oxidized and non-oxidized domains, we tested the hypothesis according to which poration may be facilitated at the domain interface region. However, results were inconclusive due to their large statistical variance and sensitivity to simulation setup parameters. We pointed out important technical issues that need to be considered in future simulations of mechanically-induced poration of heterogeneous membranes. This research is of interest for photodynamic therapy and plasma medicine, because ruptured and intact plasma membranes are experimentally considered hallmarks of necrotic and apoptotic cell death.

Oxygen distribution in the fluid/gel phases of lipid membranes.

The interactions between oxygen and lipid membranes play fundamental roles in basic biological processes (e.g., cellular respiration). Obviously, membrane oxidation is expected to be critically dependent on the distribution and concentration of oxygen in the membrane. Here, we combined theoretical and experimental methods to investigate oxygen partition and distribution in lipid membranes of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) and 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) in a temperature range between 298 and 323 K, specifically focusing on the changes caused by the lipid phase and phase transition. Even though oxygen is known to be more concentrated in the center of fluid phase membranes than on the headgroup regions, the distribution profile of oxygen inside gel-phase bilayers remained to be determined. Molecular dynamics simulations now show that the distribution of oxygen inside DPPC bilayers dramatically changes upon crossing the main transition temperature, with oxygen being nearly depleted halfway from the headgroups to the membrane center below the transition temperature. In a parallel approach, singlet oxygen luminescence emission measurements employing the photosensitizer Pheophorbide-a (Pheo) confirmed the differences in oxygen distribution and concentration profiles between gel- and fluid-phase membranes, revealing changes in the microenvironment of the embedded photosensitizer. Our results also reveal that excited triplet state lifetime, as it can be determined from the singlet oxygen luminescence kinetics, is a useful probe to assess oxygen distribution in lipid membranes with distinct lipid compositions.

Molecular Structure and Permeability at the Interface between Phase-Separated Membrane Domains.

Phase-separated membrane domains, also known as lipid rafts, are believed to play an important role in cell function. Although most rafts are sterol-enriched membrane regions, evidence suggests that living cells may also contain gel-like rafts. Interactions between gel and fluid domains have a large impact on membrane properties, as is the case with permeability. The membrane permeability may reach a peak at the main phase transition temperature, by far exceeding the values recorded at the fluid phase. It has been proposed that gel-fluid interfaces are leaky, but the effect has not yet been demonstrated at the molecular level. Here, we performed atomistic molecular dynamics simulations of phospholipid bilayers with coexisting gel-like and fluid domains. We found that the thickness mismatch between both phases, the membrane elasticity, and the lipid packing acted together to promote the formation of a thickness minimum at the gel-fluid interface. Free energy calculations showed that pore-mediated ionic permeation was strongly facilitated at the constriction region, whereas water permeation by simple diffusion was only marginally affected. Long-lived, peristaltic undulations were recorded at the bulk fluid phase near the main transition temperature. They gave rise to thickness minima that, although shallower than the interface constrictions, could also enhance permeability. Finally, we demonstrated that an interface constriction was also formed at the boundaries of regular, cholesterol-enriched lipid rafts. Our simulation results will hopefully contribute to a better understanding of biological processes such as transport, signaling, and cellular damage promoted by low temperature and dehydration.

Behavior of evaporating droplets at nonsoluble and soluble surfaces: Modeling with molecular resolution.

Liquid droplets in equilibrium with vapor are simulated at solidlike surfaces using the cooperative motion algorithm (CMA). These droplets behave like real droplets, i.e., the densities of the coexistent liquid and vapor phases obey empirical relations such as rho l - rho v proportional, variant (1 - T/Tc)(1/3). Droplet evaporation was studied under various interaction conditions, i.e., nonsoluble and soluble substrates. In the last case, substrate particles migrate toward the liquid-vapor interface to minimize the droplet surface energy. This leads to the formation of a microwell surrounded by a ringlike deposit on the substrate surface. It is shown that the ring formation in the first stages of evaporation results in pinning of the droplet contact area.

Parallel damage in mitochondrial and lysosomal compartments promotes efficient cell death with autophagy: The case of the pentacyclic triterpenoids.

The role of autophagy in cell death is still controversial and a lot of debate has concerned the transition from its pro-survival to its pro-death roles. The similar structure of the triterpenoids Betulinic (BA) and Oleanolic (OA) acids allowed us to prove that this transition involves parallel damage in mitochondria and lysosome. After treating immortalized human skin keratinocytes (HaCaT) with either BA or OA, we evaluated cell viability, proliferation and mechanism of cell death, function and morphology of mitochondria and lysosomes, and the status of the autophagy flux. We also quantified the interactions of BA and OA with membrane mimics, both in-vitro and in-silico. Essentially, OA caused mitochondrial damage that relied on autophagy to rescue cellular homeostasis, which failed upon lysosomal inhibition by Chloroquine or Bafilomycin-A1. BA caused parallel damage on mitochondria and lysosome, turning autophagy into a destructive process. The higher cytotoxicity of BA correlated with its stronger efficiency in damaging membrane mimics. Based on these findings, we underlined the concept that autophagy will turn into a destructive outcome when there is parallel damage in mitochondrial and lysosomal membranes. We trust that this concept will help the development of new drugs against aggressive cancers.

Photodynamic efficiency of cationic meso-porphyrins at lipid bilayers: insights from molecular dynamics simulations.

Porphyrin derivatives have applications as photoactive drugs in photodynamic therapy. However, little is known about their interactions with phospholipid membranes at the molecular level. We employed molecular dynamics simulations to model the binding between a series of cationic meso-(N-methyl-4-pyridinium)phenylporphyrins and anionic phosphatidylglycerol lipid bilayers. This was done in the presence of molecular oxygen within the membrane. The ability of various porphyrins to cause photodamage was quantified in terms of their immersion depth and degree of exposition to a higher oxygen concentration inside the membrane. Simulations showed that the photodynamic efficiency could be improved as the number of hydrophobic phenyl substituents attached to the porphyrinic ring increased. In the specific case of porphyrins containing two hydrophobic and two charged substituents, the cis isomer was significantly more efficient than the trans. These results correlate well with previous experimental observations. They highlight the importance of both the total charge and amphiphilicity of the photosensitizer for its performance in photodynamic therapy.