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: 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 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: In recent years there has been growing interest in the use of low-temperature atmospheric pressure plasmas for biomedical applications. Currently, however, there is very little fundamental knowledge regarding the relevant interaction mechanisms of plasma species with living cells. In this paper, we investigate the interaction of important plasma species, such as O3, O2 and O atoms, with bacterial peptidoglycan (or murein) by means of reactive molecular dynamics simulations. Specifically, we use the peptidoglycan structure to model the gram-positive bacterium Staphylococcus aureus murein. Peptidoglycan is the outer protective barrier in bacteria and can therefore interact directly with plasma species. Our results demonstrate that among the species mentioned above, O3 molecules and especially O atoms can break important bonds of the peptidoglycan structure (i.e. CO, CN and CC bonds), which subsequently leads to the destruction of the bacterial cell wall. This study is important for gaining a fundamental insight into the chemical damaging mechanisms of the bacterial peptidoglycan structure on the atomic scale.

Abstract: A multi-species Monte Carlo (MC) model, combined with an analytical surface model, has been developed in order to investigate the general plasma processes occurring during the sputter deposition of complex oxide films in a dual-magnetron sputter deposition system. The important plasma species, such as electrons, Ar+ ions, fast Ar atoms and sputtered metal atoms (i.e. Mg and Al atoms) are described with the so-called multi-species MC model, whereas the deposition of MgxAlyOz films is treated by an analytical surface model. Targetsubstrate distances for both magnetrons in the dual-magnetron setup are varied for the purpose of growing stoichiometric complex oxide thin films. The metal atoms are sputtered from pure metallic targets, whereas the oxygen flux is only directed toward the substrate and is high enough to obtain fully oxidized thin films but low enough to avoid target poisoning. The calculations correspond to typical experimental conditions applied to grow these complex oxide films. In this paper, some calculation results are shown, such as the densities of various plasma species, their fluxes toward the targets and substrate, the deposition rates, as well as the film stoichiometry. Moreover, some results of the combined model are compared with experimental observations. Note that this is the first complete model, which can be applied for large and complicated magnetron reactor geometries, such as dual-magnetron configurations. With this model, we are able to describe all important plasma species as well as the deposition process. It can also be used to predict film stoichiometries of complex oxide films on the substrate.

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: Reactive molecular dynamics simulations are performed to study the interaction of reactive oxygen species, such as OH, HO2 and H2O2, with the endotoxic biomolecule lipid A of the gram-negative bacterium Escherichia coli. It is found that the aforementioned plasma species can destroy the lipid A, which consequently results in reducing its toxic activity. All bond dissociation events are initiated by hydrogen-abstraction reactions. However, the mechanisms behind these dissociations are dependent on the impinging plasma species, i.e. a clear difference is observed in the mechanisms upon impact of HO2 radicals and H2O2 molecules on one hand and OH radicals on the other hand. Our simulation results are in good agreement with experimental observations.

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: The penetration of reactive oxygen and nitrogen species (RONS) across the stratum corneum (SC) is a necessary and crucial process in many skin‐related plasma medical applications. To gain more insights into this penetration behavior, we combined experimental measurements of the permeability of dry and moist SC layers with computer simulations of model lipid membranes. We measured the permeation of relatively stable molecules, which are typically generated by plasma, namely H2O2, NO3−, and NO2−. Furthermore, we calculated the permeation free energy profiles of the major plasma‐generated RONS and their derivatives (i.e., H2O2, OH, HO2, O2, O3, NO, NO2, N2O4, HNO2, HNO3, NO2−, and NO3−) across native and oxidized SC lipid bilayers, to understand the mechanisms of RONS permeation across the SC. Our results indicate that hydrophobic RONS (i.e., NO, NO2, O2, O3, and N2O4) can translocate more easily across the SC lipid bilayer than hydrophilic RONS (i.e., H2O2, OH, HO2, HNO2, and HNO3) and ions (i.e., NO2− and NO3−) that experience much higher permeation barriers. The permeability of RONS through the SC skin lipids is enhanced when the skin is moist and the lipids are oxidized. These findings may help to understand the underlying mechanisms of plasma interaction with a biomaterial and to optimize the environmental parameters in practice in plasma medical applications.

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

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: We report on multi-level atomistic simulations for the interaction of reactive oxygen species (ROS) with the head groups of the phospholipid bilayer, and the subsequent effect of head group and lipid tail oxidation on the structural and dynamic properties of the cell membrane. Our simulations are validated by experiments using a cold atmospheric plasma as external ROS source. We found that plasma treatment leads to a slight initial rise in membrane rigidity, followed by a strong and persistent increase in fluidity, indicating a drop in lipid order. The latter is also revealed by our simulations. This study is important for cancer treatment by therapies producing (extracellular) ROS, such as plasma treatment. These ROS will interact with the cell membrane, first oxidizing the head groups, followed by the lipid tails. A drop in lipid order might allow them to penetrate into the cell interior (e.g., through pores created due to oxidation of the lipid tails) and cause intracellular oxidative damage, eventually leading to cell death. This work in general elucidates the underlying mechanisms of ROS interaction with the cell membrane at the atomic level.

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.