Abstract: A Monte Carlo model has been developed for investigating the electron behavior in a dual-magnetron sputter deposition system. To describe the three-dimensional (3D) geometry, different reference frames, i.e. a local and a global coordinate system, were used. In this study, the influence of both closed and mirror magnetic field configurations on the plasma properties is investigated. In the case of a closed magnetic field configuration, the calculated electron trajectories show that if an electron is emitted in (or near) the center of the cathode, where the influence of the magnetic field is low, it is able to travel from one magnetron to the other. On the other hand, when an electron is created at the race track area, it is more or less trapped in the strong magnetic field and cannot easily escape to the second magnetron region. In the case of a mirror magnetic field configuration, irrespective of where the electron is emitted from the cathode, it cannot travel from one magnetron to the other because the magnetic field lines guide the electron to the substrate. Moreover, the electron density and electron impact ionization rate have been calculated and studied in detail for both configurations.

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: Bacteria can be inactivated through various physical and chemical means, and these have always been the focus of extensive research. To further improve the methodology for these ends, two types of plasma systems were investigated: nano-second pulsed plasma (NPP) as liquid discharge plasma and an Argon gas-feeding dielectric barrier discharge (Ar-DBD) as a form of surface plasma. To understand the sterilizing action of these two different plasma sources, we performed experiments with Staphylococcus aureus (S. aureus) bacteria (wild type) and multidrug resistant bacteria (Penicillum-resistant, Methicillin-resistant and Gentamicin-resistant). We observed that both plasma sources can inactivate both the wild type and multidrug-resistant bacteria to a good extent. Moreover, we observed a change in the surface morphology, gene expression and β-lactamase activity. Furthermore, we used X-ray photoelectron spectroscopy to investigate the variation in functional groups (C-H/C-C, C-OH and C=O) of the peptidoglycan (PG) resulting from exposure to plasma species. To obtain atomic scale insight in the plasma-cell interactions and support our experimental observations, we have performed molecular dynamics simulations to study the effects of plasma species, such as OH, H2O2, O, O3, as well as O2 and H2O, on the dissociation/formation of above mentioned functional groups in PG.

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: Carbon nanotubes (CNTs) have been successfully used in biomedicine, including cancer therapy, due to their unique physico-chemical properties. Because pristine CNTs exhibit hydrophobic behaviour, they can have a cytotoxic effect on cells, which limits their practical use in biomedicine. The toxicity of CNTs can be reduced by adding water-soluble functional radicals to their surface, i.e. by increasing their hydrophilicity. Another possibility for increasing the hydrophilicity of CNTs is probably filling them with endohedral metal atoms, which has not yet been studied. Thus, in this study, we use computer simulations to investigate the combined effect of endohedral nickel atoms and functional groups on the hydrophilicity of CNTs. Our simulation results show that the introduction of endohedral nickel atoms into CNTs increases their binding energy with functional groups. We also find that the addition of functional groups to the surface of CNT, along with filling it with endohedral nickel atoms, leads to an increase in the dipole moment of the CNT as well as its interaction energy with water, thereby increasing the hydrophilicity of the CNT and, consequently, its solubility in water. This, in turn, can lead to a decrease in CNT toxicity.

Abstract: Cold atmospheric plasma (CAP) is a promising technology for several medical applications, including the removal of biofilms from surfaces. However, the molecular mechanisms of CAP treatment are still poorly understood. Here we unravel the possible mechanisms of CAP‐induced oxidation of oligosaccharides, employing reactive molecular dynamics simulations based on the density functional‐tight binding potential. Specifically, we find that the interaction of oxygen atoms (used as CAP‐generated reactive species) with cellotriose (a model system for the oligosaccharides) can break structurally important glycosidic bonds, which subsequently leads to the disruption of the oligosaccharide molecule. The overall results help to shed light on our experimental evidence for cellotriose CAP. This oxidation by study provides atomic‐level insight into the onset of plasma‐induced removal of biofilms, as oligosaccharides are one of the main components of biofilm.

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

Abstract: In this review paper, we present several examples of reactive molecular dynamics simulations, which contribute to a better understanding of the underlying mechanisms in plasma medicine on the atomic scale. This includes the interaction of important reactive oxygen plasma species with the outer cell wall of both gram-positive and gram-negative bacteria, and with lipids present in human skin. Moreover, as most biomolecules are surrounded by a liquid biofilm, the behavior of these plasma species in a liquid (water) layer is presented as well. Finally, a perspective for future atomic scale modeling studies is given, in the field of plasma medicine in general, and for cancer treatment in particular.

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

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

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