In silico study of the impact of oxidation on pyruvate transmission across the hVDAC1 protein channel.

The overexpression of voltage dependent anion channels (VDACs), particularly VDAC1, in cancer cells compared to normal cells, plays a crucial role in cancer cell metabolism, apoptosis regulation, and energy homeostasis. In this study, we used molecular dynamics (MD) simulations to investigate the effect of a low level of VDAC1 oxidation (induced e.g., by cold atmospheric plasma (CAP)) on the pyruvate (Pyr) uptake by VDAC1. Inhibiting Pyr uptake through VDAC1 can suppress cancer cell proliferation. Our primary target was to study the translocation of Pyr across the native and oxidized forms of hVDAC1, the human VDAC1. Specifically, we employed MD simulations to analyze the hVDAC1 structure by modifying certain cysteine residues to cysteic acids and methionine residues to methionine sulfoxides, which allowed us to investigate the effect of oxidation. Our results showed that the free energy barrier for Pyr translocation through the native and oxidized channel was approximately 4.3 ± 0.7 kJ mol-1 and 10.8 ± 1.8 kJ mol-1, respectively. An increase in barrier results in a decrease in rate of Pyr permeation through the oxidized channel. Thus, our results indicate that low levels of CAP oxidation reduce Pyr translocation, resulting in decreased cancer cell proliferation. Therefore, low levels of oxidation are likely sufficient to treat cancer cells given the inhibition of Pyr uptake.

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.

Unraveling the Transport Properties of RONS across Nitro-Oxidized Membranes.

The potential of cold atmospheric plasma (CAP) in biomedical applications has received significant interest, due to its ability to generate reactive oxygen and nitrogen species (RONS). Upon exposure to living cells, CAP triggers alterations in various cellular components, such as the cell membrane. However, the permeation of RONS across nitrated and oxidized membranes remains understudied. To address this gap, we conducted molecular dynamics simulations, to investigate the permeation capabilities of RONS across modified cell membranes. This computational study investigated the translocation processes of less hydrophilic and hydrophilic RONS across the phospholipid bilayer (PLB), with various degrees of oxidation and nitration, and elucidated the impact of RONS on PLB permeability. The simulation results showed that less hydrophilic species, i.e., NO, NO2, N2O4, and O3, have a higher penetration ability through nitro-oxidized PLB compared to hydrophilic RONS, i.e., HNO3, s-cis-HONO, s-trans-HONO, H2O2, HO2, and OH. In particular, nitro-oxidation of PLB, induced by, e.g., cold atmospheric plasma, has minimal impact on the penetration of free energy barriers of less hydrophilic species, while it lowers these barriers for hydrophilic RONS, thereby enhancing their translocation across nitro-oxidized PLB. This research contributes to a better understanding of the translocation abilities of RONS in the field of plasma biomedical applications and highlights the need for further analysis of their role in intracellular signaling pathways.

Does non-thermal plasma modify biopolymers in solution? A chemical and mechanistic study for alginate.

In the last decades, non-thermal plasma has been extensively investigated as a relevant tool for various biomedical applications, ranging from tissue decontamination to regeneration and from skin treatment to tumor therapies. This high versatility is due to the different kinds and amount of reactive oxygen and nitrogen species that can be generated during a plasma treatment and put in contact with the biological target. Some recent studies report that solutions of biopolymers with the ability to generate hydrogels, when treated with plasma, can enhance the generation of reactive species and influence their stability, resulting thus in the ideal media for indirect treatments of biological targets. The direct effects of the plasma treatment on the structure of biopolymers in water solution, as well as the chemical mechanisms responsible for the enhanced generation of RONS, are not yet fully understood. In this study, we aim at filling this gap by investigating, on the one hand, the nature and extent of the modifications induced by plasma treatment in alginate solutions, and, on the other hand, at using this information to explain the mechanisms responsible for the enhanced generation of reactive species as a consequence of the treatment. The approach we use is twofold: (i) investigating the effects of plasma treatment on alginate solutions, by size exclusion chromatography, rheology and scanning electron microscopy and (ii) study of a molecular model (glucuronate) sharing its chemical structure, by chromatography coupled with mass spectrometry and by molecular dynamics simulations. Our results point out the active role of the biopolymer chemistry during direct plasma treatment. Short-lived reactive species, such as OH radicals and O atoms, can modify the polymer structure, affecting its functional groups and causing partial fragmentation. Some of these chemical modifications, like the generation of organic peroxide, are likely responsible for the secondary generation of long-lived reactive species such as hydrogen peroxide and nitrite ions. This is relevant in view of using biocompatible hydrogels as vehicles for storage and delivery reactive species for targeted therapies.

Bioactive Nonthermal Biocompatible Plasma Enhances Migration on Human Gingival Fibroblasts.

This study hypothesizes that the application of low-dose nonthermal biocompatible dielectric barrier discharge plasma (DBD-NBP) to human gingival fibroblasts (HGFs) will inhibit colony formation but not cell death and induce matrix metalloproteinase (MMP) expression, extracellular matrix (ECM) degradation, and subsequent cell migration, which can result in enhanced wound healing. HGFs treated with plasma for 3 min migrate to each other across the gap faster than those in the control and 5-min treatment groups on days 1 and 3. The plasma-treated HGFs show significantly high expression levels of the cell cycle arrest-related p21 gene and enhanced MMP activity. Focal adhesion kinase (FAK) mediated attenuation of wound healing or actin cytoskeleton rearrangement, and plasma-mediated reversal of this attenuation support the migratory effect of DBD-NBP. Further, this work performs computer simulations to investigate the effect of oxidation on the stability and conformation of the catalytic kinase domain (KD) of FAK. It is found that the oxidation of highly reactive amino acids (AAs) Cys427, Met442, Cys559, Met571, Met617, and Met643 changes the conformation and increases the structural flexibility of the FAK protein and thus modulates its function and activity. Low-dose DBD-NBP-induces host cell cycle arrest, ECM breakdown, and subsequent migration, thus contributing to the enhanced wound healing process.

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.

Effect of Cysteine Oxidation in SARS-CoV-2 Receptor-Binding Domain on Its Interaction with Two Cell Receptors: Insights from Atomistic Simulations.

Binding of the SARS-CoV-2 S-glycoprotein to cell receptors is vital for the entry of the virus into cells and subsequent infection. ACE2 is the main cell receptor for SARS-CoV-2, which can attach to the C-terminal receptor-binding domain (RBD) of the SARS-CoV-2 S-glycoprotein. The GRP78 receptor plays an anchoring role, which attaches to the RBD and increases the chance of other RBDs binding to ACE2. Although high levels of reactive oxygen and nitrogen species (RONS) are produced during viral infections, it is not clear how they affect the RBD structure and its binding to ACE2 and GRP78. In this research, we apply molecular dynamics simulations to study the effect of oxidation of the highly reactive cysteine (Cys) amino acids of the RBD on its binding to ACE2 and GRP78. The interaction energy of both ACE2 and GRP78 with the whole RBD, as well as with the RBD main regions, is compared in both the native and oxidized RBDs. Our results show that the interaction energy between the oxidized RBD and ACE2 is strengthened by 155 kJ/mol, increasing the binding of the RBD to ACE2 after oxidation. In addition, the interaction energy between the RBD and GRP78 is slightly increased by 8 kJ/mol after oxidation, but this difference is not significant. Overall, these findings highlight the role of RONS in the binding of the SARS-CoV-2 S-glycoprotein to host cell receptors and suggest an alternative mechanism by which RONS could modulate the entrance of viral particles into the cells.

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.

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.

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.

Oxidation of Innate Immune Checkpoint CD47 on Cancer Cells with Non-Thermal Plasma.

Non-thermal plasma (NTP) therapy has been emerging as a promising cancer treatment strategy, and recently, its ability to locally induce immunogenic cancer cell death is being unraveled. We hypothesized that the chemical species produced by NTP reduce immunosuppressive surface proteins and checkpoints that are overexpressed on cancerous cells. Here, 3D in vitro tumor models, an in vivo mouse model, and molecular dynamics simulations are used to investigate the effect of NTP on CD47, a key innate immune checkpoint. CD47 is immediately modulated after NTP treatment and simulations reveal the potential oxidized salt-bridges responsible for conformational changes. Umbrella sampling simulations of CD47 with its receptor, signal-regulatory protein alpha (SIRPα), demonstrate that the induced-conformational changes reduce its binding affinity. Taken together, this work provides new insight into fundamental, chemical NTP-cancer cell interaction mechanisms and a previously overlooked advantage of present NTP cancer therapy: reducing immunosuppressive signals on the surface of cancer cells.

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.

Influence of osmolytes and ionic liquids on the Bacteriorhodopsin structure in the absence and presence of oxidative stress: A combined experimental and computational study.

Understanding the folding and stability of membrane proteins is of great importance in protein science. Recently, osmolytes and ionic liquids (ILs) are increasingly being used as drug delivery systems in the biopharmaceutical industry. However, the stability of membrane proteins in the presence of osmolytes and ILs is not yet fully understood. Besides, the effect of oxidative stress on membrane proteins with osmolytes or ILs has not been investigated. Therefore, we studied the influence of osmolytes and ILs as co-solvents on the stability of a model membrane protein (i.e., Bacteriorhodopsin in purple membrane of Halobacterium salinarum), using UV-Vis spectroscopy and molecular dynamics (MD) simulations. The MD simulations allowed us to determine the flexibility and solvent accessible surface area (SASA) of Bacteriorhodopsin protein in the presence and/or absence of co-solvents, as well as to carry out principal component analysis (PCA) to identify the most important movements in this protein. In addition, by means of UV-Vis spectroscopy we studied the effect of oxidative stress generated by cold atmospheric plasma on the stability of Bacteriorhodopsin in the presence and/or absence of co-solvents. This study is important for a better understanding of the stability of proteins in the presence of oxidative stress.

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.

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.

Effect of oxidative stress on cystine transportation by xC‾ antiporter.

We performed computer simulations to investigate the effect of oxidation on the extracellular cystine (CYC) uptake by the xC- antiporter. The latter is important for killing of cancer cells. Specifically, applying molecular dynamics (MD) simulations we studied the transport of CYC across xCT, i.e., the light subunit of the xC- antiporter, in charge of bidirectional transport of CYC and glutamate. We considered the outward facing (OF) configuration of xCT, and to study the effect of oxidation, we modified the Cys327 residue, located in the vicinity of the extracellular milieu, to cysteic acid (CYO327). Our computational results showed that oxidation of Cys327 results in a free energy barrier for CYC translocation, thereby blocking the access of CYC to the substrate binding site of the OF system. The formation of the energy barrier was found to be due to the conformational changes in the channel. Analysis of the MD trajectories revealed that the reorganization of the side chains of the Tyr244 and CYO327 residues play a critical role in the OF channel blocking. Indeed, the calculated distance between Tyr244 and either Cys327 or CYO327 showed a narrowing of the channel after oxidation. The obtained free energy barrier for CYC translocation was found to be 33.9kJmol-1, indicating that oxidation of Cys327, by e.g., cold atmospheric plasma, is more effective in inhibiting the xC- antiporter than in the mutation of this amino acid to Ala (yielding a barrier of 32.4kJmol-1). The inhibition of the xC- antiporter may lead to Cys starvation in some cancer cells, eventually resulting in cancer cell death.

How membrane lipids influence plasma delivery of reactive oxygen species into cells and subsequent DNA damage: an experimental and computational study.

The mechanisms of plasma in medicine are broadly attributed to plasma-derived reactive oxygen and nitrogen species (RONS). In order to exert any intracellular effects, these plasma-derived RONS must first traverse a major barrier in the cell membrane. The cell membrane lipid composition, and thereby the magnitude of this barrier, is highly variable between cells depending on type and state (e.g. it is widely accepted that healthy and cancerous cells have different membrane lipid compositions). In this study, we investigate how plasma-derived RONS interactions with lipid membrane components can potentially be exploited in the future for treatment of diseases. We couple phospholipid vesicle experiments, used as simple cell models, with molecular dynamics (MD) simulations of the lipid membrane to provide new insights into how the interplay between phospholipids and cholesterol may influence the response of healthy and diseased cell membranes to plasma-derived RONS. We focus on the (i) lipid tail saturation degree, (ii) lipid head group type, and (iii) membrane cholesterol fraction. Using encapsulated molecular probes, we study the influence of the above membrane components on the ingress of RONS into the vesicles, and subsequent DNA damage. Our results indicate that all of the above membrane components can enhance or suppress RONS uptake, depending on their relative concentration within the membrane. Further, we show that higher RONS uptake into the vesicles does not always correlate with increased DNA damage, which is attributed to ROS reactivity and lifetime. The MD simulations indicate the multifactorial chemical and physical processes at play, including (i) lipid oxidation, (ii) lipid packing, and (iii) lipid rafts formation. The methods and findings presented here provide a platform of knowledge that could be leveraged in the development of therapies relying on the action of plasma, in which the cell membrane and oxidative stress response in cells is targeted.

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.

Transport of cystine across xC- antiporter.

Extracellular cystine (CYC) uptake by xC- antiporter is important for the cell viability. Especially in cancer cells, the upregulation of xC- activity is observed, which protects these cells from intracellular oxidative stress. Hence, inhibition of the CYC uptake may eventually lead to cancer cell death. Up to now, the molecular level mechanism of the CYC uptake by xC- antiporter has not been studied in detail. In this study, we applied several different simulation techniques to investigate the transport of CYC through xCT, the light subunit of the xC- antiporter, which is responsible for the CYC and glutamate translocation. Specifically, we studied the permeation of CYC across three model systems, i.e., outward facing (OF), occluded (OCC) and inward facing (IF) configurations of xCT. We also investigated the effect of mutation of Cys327 to Ala within xCT, which was also studied experimentally in literature. This allowed us to qualitatively compare our computation results with experimental observations, and thus, to validate our simulations. In summary, our simulations provide a molecular level mechanism of the transport of CYC across the xC- antiporter, more specifically, which amino acid residues in the xC- antiporter play a key role in the uptake, transport and release of CYC.

The effect of reactive oxygen and nitrogen species on the structure of cytoglobin: A potential tumor suppressor.

Many current anti-cancer therapies rely on increasing the intracellular reactive oxygen and nitrogen species (RONS) contents with the aim to induce irreparable damage, which subsequently results in tumor cell death. A novel tool in cancer therapy is the use of cold atmospheric plasma (CAP), which has been found to be very effective in the treatment of many different cancer cell types in vitro as well as in vivo, mainly through the vast generation of RONS. One of the key determinants of the cell's fate will be the interaction of RONS, generated by CAP, with important proteins, i.e. redox-regulatory proteins. One such protein is cytoglobin (CYGB), a recently discovered globin proposed to be involved in the protection of the cell against oxidative stress. In this study, the effect of plasma-produced RONS on CYGB was investigated through the treatment of CYGB with CAP for different treatment times. Spectroscopic analysis of CYGB showed that although chemical modifications occur, its secondary structure remains intact. Mass spectrometry experiments identified these modifications as oxidations of mainly sulfur-containing and aromatic amino acids. With longer treatment time, the treatment was also found to induce nitration of the heme. Furthermore, the two surface-exposed cysteine residues of CYGB were oxidized upon treatment, leading to the formation of intermolecular disulfide bridges, and potentially also intramolecular disulfide bridges. In addition, molecular dynamics and docking simulations confirmed, and further show, that the formation of an intramolecular disulfide bond, due to oxidative conditions, affects the CYGB 3D structure, thereby opening the access to the heme group, through gate functioning of His117. Altogether, the results obtained in this study (1) show that plasma-produced RONS can extensively oxidize proteins and (2) that the oxidation status of two redox-active cysteines lead to different conformations of CYGB.