Abstract: The application of nonthermal atmospheric pressure plasma is emerging as an alternative and efficient technique for the inactivation of bacterial biofilms. In this study, reactive molecular dynamics simulations were used to examine the reaction mechanisms of hydroxyl radicals, as key reactive oxygen plasma species in biological systems, with several organic molecules (i.e., alkane, alcohol, carboxylic acid, and amine), as prototypical components of biomolecules in the biofilm. Our results demonstrate that organic molecules containing hydroxyl and carboxyl groups may act as trapping agents for the OH radicals. Moreover, the impact of OH radicals on N-acetyl-glucosamine, as constituent component of staphylococcus epidermidis biofilms, was investigated. The results show how impacts of OH radicals lead to hydrogen abstraction and subsequent molecular damage. This study thus provides new data on the reaction mechanisms of plasma species, and particularly the OH radicals, with fundamental components of bacterial biofilms.

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

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: 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: Plasma medicine has attracted tremendous interest in a variety of medical conditions, ranging from wound healing to antimicrobial applications, even in cancer treatment, through the interactions of cold atmospheric plasma (CAP) and various biological tissues directly or indirectly. The underlying mechanisms of CAP treatment are still poorly understood although the oxidative effects of CAP with amino acids, peptides, and proteins have been explored experimentally. In this study, machine learning (ML) technology is introduced to efficiently unveil the interaction mechanisms of amino acids and reactive oxygen species (ROS) in seconds based on the data obtained from the reactive molecular dynamics (MD) simulations, which are performed to probe the interaction of five types of amino acids with various ROS on the timescale of hundreds of picoseconds but with the huge computational load of several days. The oxidative reactions typically start with H-abstraction, and the details of the breaking and formation of chemical bonds are revealed; the modification types, such as nitrosylation, hydroxylation, and carbonylation, can be observed. The dose effects of ROS are also investigated by varying the number of ROS in the simulation box, indicating agreement with the experimental observation. To overcome the limits of timescales and the size of molecular systems in reactive MD simulations, a deep neural network (DNN) with five hidden layers is constructed according to the reaction data and employed to predict the type of oxidative modification and the probability of occurrence only in seconds as the dose of ROS varies. The well-trained DNN can effectively and accurately predict the oxidative processes and productions, which greatly improves the computational efficiency by almost ten orders of magnitude compared with the reactive MD simulation. This study shows the great potential of ML technology to efficiently unveil the underpinning mechanisms in plasma medicine based on the data from reactive MD simulations or experimental measurements. In this study, since reactive molecular dynamics simulation can currently only describe interactions between a few hundred atoms in a few hundred picoseconds, deep neural networks (DNN) are introduced to enhance the simulation results by predicting more data efficiently. image

Abstract: Device performance of Si|SiOx core-shell based nanowires critically depends on the exact control over the oxide thickness. Low-temperature plasma oxidation is a highly promising alternative to thermal oxidation allowing for improved control over the oxidation process, in particular for ultrasmall Si nanowires. We here elucidate the room temperature plasma oxidation mechanisms of ultrasmall Si nanowires using hybrid molecular dynamics / force-bias Monte Carlo simulations. We demonstrate how the oxidation and concurrent water formation mechanisms are a function of the oxidizing plasma species and we demonstrate how the resulting core-shell oxide thickness can be controlled through these species. A new mechanism of water formation is discussed in detail. The results provide a detailed atomic level explanation of the oxidation process of highly curved Si surfaces. These results point out a route toward plasma-based formation of ultrathin core-shell Si|SiOx nanowires at room temperature.

Abstract: Background: Strong electric fields are knownto affect cell membrane permeability,which can be applied for therapeutic purposes, e.g., in cancer therapy. A synergistic enhancement of this effect may be accomplished by the presence of reactive oxygen species (ROS), as generated in cold atmospheric plasmas. Little is known about the synergy between lipid oxidation by ROS and the electric field, nor on howthis affects the cell membrane permeability.
Method: We here conduct molecular dynamics simulations to elucidate the dynamics of the permeation process under the influence of combined lipid oxidation and electroporation. A phospholipid bilayer (PLB), consisting of di-oleoyl-phosphatidylcholine molecules covered with water layers, is used as a model system for the plasma membrane.
Results and conclusions:Weshow howoxidation of the lipids in the PLB leads to an increase of the permeability of the bilayer to ROS, although the permeation free energy barriers still remain relatively high. More importantly, oxidation of the lipids results in a drop of the electric field threshold needed for pore formation (i.e., electroporation) in the PLB. The created pores in the membrane facilitate the penetration of reactive plasma species deep into the cell interior, eventually causing oxidative damage.
General significance: This study is of particular interest for plasma medicine, as plasma generates both ROS and electric fields, but it is also of more general interest for applications where strong electric fields and ROS both come into play.

Abstract: Experiments have demonstrated the potential selective anticancer capacity of cold atmospheric plasmas (CAPs), but the underlying mechanisms remain unclear. Using computer simulations, we try to shed light on the mechanism of selectivity, based on aquaporins (AQPs), i.e. transmembrane protein channels transferring external H 2 O 2 and other reactive oxygen species, created e.g. by CAPs, to the cell interior. Specifically, we perform molecular dynamics simulations for the permeation of H 2 O 2 through AQP1 (one of the members of the AQP family) and the palmitoyl-oleoyl-phosphatidylcholine (POPC) phospholipid bilayer (PLB). The free energy barrier of H 2 O 2 across AQP1 is lower than for the POPC PLB, while the permeability coefficient, calculated using the free energy and diffusion rate profiles, is two orders of magnitude higher. This indicates that the delivery of H 2 O 2 into the cell interior should be through AQP. Our study gives a better insight into the role of AQPs in the selectivity of CAPs for treating cancer cells.

Abstract: We perform molecular dynamics simulations to study the flip-flop motion of phosphatidylserine (PS) across the plasma membrane upon increasing oxidation degree of the membrane. Our computational results show that an increase of the oxidation degree in the lipids leads to a decrease of the free energy barrier for translocation of PS through the membrane. In other words, oxidation of the lipids facilitates PS flip-flop motion across the membrane, because in native phospholipid bilayers this is only a “rare event” due to the high energy barriers for the translocation of PS. The present study provides an atomic-scale insight into the mechanisms of the PS flip-flop upon oxidation of lipids, as produced for example by cold atmospheric plasma, in living cells.

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: We perform computer simulations supported by experiments to investigate the oxidation of an important signaling protein, that is, human epidermal growth factor (hEGF), caused by cold atmospheric plasma (CAP) treatment. Specifically, we study the conformational changes of hEGF with different degrees of oxidation, to mimic short and long CAP treatment times. Our results indicate that the oxidized structures become more flexible, due to their conformational changes and breakage of the disulfide bonds, especially at higher oxidation degrees. MM/GBSA calculations reveal that an increasing oxidation level leads to a lower binding free energy of hEGF with its receptor. These results help to understand the fundamentals of the use of CAP for wound healing versus cancer treatment at short and longer treatment times.

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

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.