Abstract: 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 (H₂O₂and OH) and hydrophobic (NO₂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 H₂O₂and NO₂, 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 Cys₁₉₁) on the transport of the above-mentioned RONS across AQP1. Both lipid and protein oxidation seem to slightly increase the free energy barrier for H₂O₂and NO₂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.

Abstract: 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.

Abstract: Biological membranes are under constant attack of free radicals, which may lead to lipid nitro-oxidation, pro­ ducing 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 neuro­ degenerative 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 phos­pholipid 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.

Abstract: 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.

Abstract: Natural deep eutectic solvents (NADES) represent a green alternative to common organic solvents in the biochemical industry due to their benign behavior and tailorable properties, in particular as media for enzymatic reactions. However, to fully exploit their potential in enzymatic reactions, there is a need for a more fundamental understanding of how these neoteric solvents influence the course of these reac-tions. Thus, the aim of this study is to investigate the influence of NADES with various molar composi-tions on the stability and structure of enzymes, applying molecular dynamics simulations. This can help to better understand the effect of individual compounds of NADES, in addition to eutectic mixtures. More specifically, we simulate the behavior of Candida antarctica lipase B (CALB) enzyme in NADES com-posed of choline chloride with either urea, ethylene glycol or glycerol. Hereto, we monitor the NADES microstructure, the general stability of the enzyme and changes in the structure of its active sites and sur-face residues. Our simulations show that none of the studied NADES systems significantly disrupt the microstructure of the solvent or the stability of the CALB enzyme within the time scales of the simula-tions. The enzyme preserves its initial structure, size and intra-chain hydrogen bonds in all investigated compositions and, for the first time reported, also in NADES with increased hydrogen bond donating com-pound ratios. As the main novelty, our results indicate that, in addition to the composition, the molar ratio can be an additional variable to fine-tune the physicochemical properties of NADES without altering the enzyme characteristics. These findings could facilitate the development and application of task -tailored NADES media for biocatalytic processes. (c) 2022 Elsevier B.V. All rights reserved.

Abstract: The adsorption of CO2 on surfaces of graphene and carbon nanotubes (CNTs), decorated with Fe atoms, are investigated using the self-consistent-charge density-functional tight-binding (SCC-DFTB) method, neglecting the heat effects. Fe ad-atoms are more stable when they are dispersed on hollow sites. They introduce a large density of states at the Fermi level (N-F); where keeping such density low would help in gas sensing. Furthermore, the Fe ad-atom can weaken the C=O double bonds of the chemisorbed CO2 molecule, paving the way for oxygen atoms to drain more charges from Fe. Consequently, chemisorption of CO2 molecules reduces both N-F and the conductance while it enhances the sensitivity with the increasing gas dose. Conducting armchair CNTs (ac-CNTs) have higher sensitivity than graphene and semiconducting zigzag CNTs (zz-CNT5). Comparative study of sensitivity of ac-CNT-Fe composite towards various gases (e.g., O-2, N-2, H-2, H2O, CO and CO2) has shown high sensitivity and selectivity towards CO, CO2 and H2O gases. (C) 2016 Elsevier B.V. All rights reserved.

Abstract: Reactive oxygen and nitrogen species (RONS) are involved in many biochemical processes, including nitrooxidative 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 perme­ation 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 hy­drophilic 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. Hy­drophobic RONS like molecular oxygen (O2), nitrogen dioxide (NO2) and nitric oxide (NO) even have insignif­icant 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.

Abstract: 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.

Abstract: Nonthermal atmospheric pressure plasmas are gaining increasing attention for biomedical applications. However, very little fundamental information on the interaction mechanisms between the plasma species and biological cells is currently available. We investigate the interaction of important plasma species, such as OH, H2O2, O, O3, as well as O2 and H2O, with bacterial peptidoglycan by means of reactive molecular dynamics simulations, aiming for a better understanding of plasma disinfection. Our results show that OH, O, O3, and H2O2 can break structurally important bonds of peptidoglycan (i.e., CO, CN, or CC bonds), which consequently leads to the destruction of the bacterial cell wall. The mechanisms behind these breakups are, however, dependent on the impinging plasma species, and this also determines the effectiveness of the cell wall destruction.

Abstract: 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.

Abstract: The safe and efficient use of hydrogen energy, which is in high demand worldwide today, requires efficient hydrogen storage. Despite significant advances in hydrogen storage using carbon-based nanomaterials, including carbon nanotubes (CNTs), efforts to substantially increase the storage capacity remain less effective. In this work, we demonstrate the effect of endohedral transition metal atoms on the hydrogen storage capacity of CNTs using reactive molecular dynamics simulations. We find that an increase in the volume fraction of endohedral nickel atoms leads to an increase in the concentration of physisorbed hydrogen molecules around single-walled CNTs (SWNTs) by approximately 1.6 times compared to pure SWNTs. The obtained results provide insight into the underlying mechanisms of how endohedral transition metal atoms enhance the hydrogen storage ability of SWNTs under nearly ambient conditions.

Abstract: Cold atmospheric plasmas (CAPs) have attracted significant interest for their potential benefits in medical applications, including cancer therapy. The therapeutic effects of CAPs are related to reactive oxygen and nitrogen species (ROS and RNS) present in the plasma. The impact of ROS has been extensively studied, but the role of RNS in CAP-treatment remains poorly understood at the molecular level. Here, we investigate the permeation of RNS and ROS across native and oxidized phospholipid bilayers (PLBs) by means of computer simulations. The results reveal significantly lower free energy barriers for RNS (i.e. NO, NO2, N2O4) and O3 compared to hydrophilic ROS, such as OH, HO2 and H2O2. This suggests that the investigated RNS and O3 can permeate more easily through both native and oxidized PLBs in comparison to hydrophilic ROS, indicating their potentially important role in plasma medicine.

Abstract: While plasma treatment of skin diseases and wound healing has been proven highly effective, the underlying mechanisms, and more generally the effect of plasma radicals on skin tissue, are not yet completely understood. In this paper, we perform ReaxFF-based reactive molecular dynamics simulations to investigate the interaction of plasma generated OH radicals with a model system composed of free fatty acids, ceramides, and cholesterol molecules. This model system is an approximation of the upper layer of the skin (stratum corneum). All interaction mechanisms observed in our simulations are initiated by H-abstraction from one of the ceramides. This reaction, in turn, often starts a cascade of other reactions, which eventually lead to the formation of aldehydes, the dissociation of ceramides or the elimination of formaldehyde, and thus eventually to the degradation of the skin barrier function.

Abstract: The application of atmospheric pressure plasmas in medicine is increasingly gaining attention in recent years, although very little is currently known about the plasma-induced processes occurring on the surface of living organisms. It is known that most bio-organisms, including bacteria, are coated by a liquid film surrounding them, and there might be many interactions between plasma species and the liquid layer before the plasma species reach the surface of the bio-organisms. Therefore, it is essential to study the behaviour of the reactive species in a liquid film, in order to determine whether these species can travel through this layer and reach the biomolecules, or whether new species are formed along the way. In this work, we investigate the interaction of reactive oxygen species (i.e. O, OH, HO2 and H2O2) with water, which is assumed as a simple model system for the liquid layer surrounding biomolecules. Our computational investigations show that OH, HO2 and H2O2 can travel deep into the liquid layer and are hence in principle able to reach the bio-organism. Furthermore, O, OH and HO2 radicals react with water molecules through hydrogen-abstraction reactions, whereas no H-abstraction reaction takes place in the case of H2O2. This study is important to gain insight into the fundamental operating mechanisms in plasma medicine, in general, and the interaction mechanisms of plasma species with a liquid film, in particular.