Abstract: The peculiarities of sputtering processes at 0.5-5 keV Ne grazing ion bombardment of Si(001) and SiC(001) surfaces and their possible application for the surface modification have been studied by computer simulation. Sputtering yields in the primary knock-on recoil atoms regime versus the initial energy of incident ions (E(0) = 0.5-5 keV) and angle of incidence (psi = 0-30 degrees) counted from a target surface have been calculated. Comparative studies of layer-by-layer sputtering for Si(001) and SiC(001) surfaces versus the initial energy of incident ions as well as an effective sputtering and sputtering threshold are discussed.

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

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

Abstract: Plasma is gaining increasing interest for cancer
treatment, but the underlying mechanisms are not yet fully
understood. Using computer simulations at the molecular
level, we try to gain better insight in how plasma-generated
reactive oxygen and nitrogen species (RONS) can
penetrate through the cell membrane. Specifically, we
compare the permeability of various (hydrophilic and
hydrophobic) RONS across both oxidized and nonoxidized cell membranes. We also study pore formation,
and how it is hampered by higher concentrations of
cholesterol in the cell membrane, and we illustrate the
much higher permeability of H2O2 through aquaporin
channels. Both mechanisms may explain the selective
cytotoxic effect of plasma towards cancer cells. Finally, we
also discuss the synergistic effect of plasma-induced
oxidation and electric fields towards pore formation.
Keywords plasma medicine, cancer treatment, computer
modelling, cell membrane, reactive oxygen and nitrogen
species

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

Abstract: Modelling at the molecular or atomic scale can be very useful for obtaining a better insight in plasma medicine. This paper gives an overview of different atomic/molecular scale modelling approaches that can be used to study the direct interaction of plasma species with biomolecules or the consequences of these interactions for the biomolecules on a somewhat longer time-scale. These approaches include density functional theory (DFT), density functional based tight binding (DFTB), classical reactive and non-reactive molecular dynamics (MD) and united-atom or coarse-grained MD, as well as hybrid quantum mechanics/molecular mechanics (QM/MM) methods. Specific examples will be given for three important types of biomolecules, present in human cells, i.e. proteins, DNA and phospholipids found in the cell membrane. The results show that each of these modelling approaches has its specific strengths and limitations, and is particularly useful for certain applications. A multi-level approach is therefore most suitable for obtaining a global picture of the plasma–biomolecule interactions.

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: 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: 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: 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: 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: The unique physical properties of carbyne, a novel carbon nanostructure, have attracted considerable interest in modern nanotechnology. While carbyne synthesis has been accomplished successfully using diverse techniques, the underlying mechanisms governing the carbon monoxide-dependent catalytic synthesis of endohedral carbyne remain poorly understood. In this simulation-based study, we investigate the synthesis of endohedral carbyne from carbon and carbon monoxide radicals in the presence of a nickel catalyst inside double-walled carbon nanotubes with a (5,5)@(10,10) structure. The outcome of our investigation demonstrates that the incorporation of the carbon atom within the Ni-n@(5,5)@(10,10) model system initiates the formation of an elongated carbon chain. In contrast, upon the introduction of carbon monoxide radicals, the growth of the carbyne chain is inhibited as a result of the oxidation of endohedral nickel clusters by oxygen atoms after the initial steps of nucleation. Our findings align with prior theoretical, simulation, and experimental investigations, reinforcing their consistency and providing valuable insights into the synthesis of carbyne-based nanodevices that hold promising potential for future advancements in nanotechnology.