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: 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: 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: 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: 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: 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: 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: 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: 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: 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: Cold atmospheric plasma (CAP) has shown its potential in biomedical applications, such as wound healing, cancer treatment and bacterial disinfection. Recent experiments have provided evidence that CAP can also enhance the intracellular uptake of glucose molecules which is important in diabetes therapy. In this respect, it is essential to understand the underlying mechanisms of intracellular glucose uptake induced by CAP, which is still unclear. Hence, in this study we try to elucidate the possible mechanism of glucose uptake by cells by performing computer simulations. Specifically, we study the transport of glucose molecules through native and oxidized membranes. Our simulation results show that the free energy barrier for the permeation of glucose molecules across the membrane decreases upon increasing the degree of oxidized lipids in the membrane. This indicates that the glucose permeation rate into cells increases when the CAP oxidation level in the cell membrane is increased.

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