Abstract: This paper gives an overview of the modeling work developed in our group in the last 25 years for various plasmas used in analytical spectrochemistry, i.e., glow discharges (GDs), inductively coupled plasmas (ICPs), and laser ablation (LA) for sample introduction in the ICP and for laser-induced breakdown spectroscopy (LIBS). The modeling approaches are briefly presented, which are different for each case, and some characteristic results are illustrated. These plasmas are used not only in analytical chemistry but also in other applications, and the insights obtained in these other fields were quite helpful for us to develop models for the analytical plasmas. Likewise, there is now a huge interest in plasma–liquid interaction, atmospheric pressure glow discharges (APGDs), and dielectric barrier discharges (DBDs) for environmental, medical, and materials applications of plasmas. The insights obtained in these fields are also very relevant for ambient desorption/ionization sources and for liquid sampling, which are nowadays very popular in analytical chemistry, and they could be very helpful in developing models for these sources as well.

Abstract: Nonthermal plasma is a promising alternative for ammonia synthesis at gentle conditions. Metal meshes of Fe, Cu, Pd, Ag, and Au were employed as catalysts in radio frequency plasma for ammonia synthesis. The energy yield for all these transition metal catalysts ranged between 0.12 and 0.19 g-NH3/kWh at 300 W and, thus, needs further improvement. In addition, a semimetal, pure gallium, was used for the first time as catalyst for ammonia synthesis, with energy yield of 0.22 g-NH3/kWh and with a maximum yield of ∼10% at 150 W. The emission spectra, as well as computer simulations, revealed hydrogen recombination as a primary governing parameter, which depends on the concentration or flux of H atoms in the plasma and on the catalyst surface. The simulations helped to elucidate the underlying mechanism, implicating the dominance of surface reactions and surface adsorbed species. The rate limiting step appears to be NH2 formation on the surface of the reactor wall and on the catalyst surface, which is different from classical catalysis.

Abstract: The aggregation of insoluble amyloid beta (Aβ) peptides in the brain is known to trigger the onset of neurodegenerative diseases, such as Alzheimer’s disease. In spite of the massive number of investigations, the underlying mechanisms to destabilize the Aβ aggregates are still poorly understood. Some studies indicate the importance of oxidation to destabilize the Aβ aggregates. In particular, oxidation induced by cold atmospheric plasma (CAP) has demonstrated promising results in eliminating these toxic aggregates. In this paper, we investigate the effect of oxidation on the stability of an Aβ pentamer. By means of molecular dynamics simulations and umbrella sampling, we elucidate the conformational changes of Aβ pentamer in the presence of oxidized residues, and we estimate the dissociation free energy of the terminal peptide out of the pentamer form. The calculated dissociation free energy of the terminal peptide is also found to decrease with increasing oxidation. This indicates that Aβ pentamer aggregation becomes less favorable upon oxidation. Our study contributes to a better insight in one of the potential mechanisms for inhibition of toxic Aβ peptide aggregation, which is considered to be the main culprit to Alzheimer’s disease.

Abstract: The modelling of a gliding arc discharge (GAD) is studied by means of the quasineutral (QN) plasma modelling approach. The model is first evaluated for reliability and proper description of a gliding arc discharge at atmospheric pressure, by comparing with a more elaborate non-quasineutral (NQN) plasma model in two different geometries – a 2D axisymmetric and a Cartesian geometry. The NQN model is considered as a reference, since it provides a continuous self-consistent plasma description, including the near electrode regions. In general, the results of the QN model agree very well with those obtained from the NQN model. The small differences between both models are attributed to the approximations in the derivation of the QN model. The use of the QN model provides a substantial reduction of the computation time compared to the NQN model, which is crucial for the development of more complex models in three dimensions or with complicated chemistries. The latter is illustrated for (i) a reverse vortex flow(RVF) GAD in argon, and (ii) a GAD in CO2. The RVF discharge is modelled in three dimensions and the effect of the turbulent heat transport on the plasma and gas characteristics is
discussed. The GAD model in CO2 is in a 1D geometry with axial symmetry and provides results for the time evolution of the electron, gas and vibrational temperature of CO2, as well as for the molar fractions of the different species.

Abstract: We performed computer simulations to investigate the effect of oxidation on the extracellular cystine (CYC) uptake by the xC− antiporter. The latter is important for killing of cancer cells. Specifically, applying molecular dynamics (MD) simulations we studied the transport of CYC across xCT, i.e., the light subunit of the xC− antiporter, in charge of bidirectional transport of CYC and glutamate. We considered the outward facing (OF) configuration of xCT, and to study the effect of oxidation, we modified the Cys327 residue, located in the vicinity of the extracellular milieu, to cysteic acid (CYO327). Our computational results showed that oxidation of Cys327 results in a free energy barrier for CYC translocation, thereby blocking the access of CYC to the substrate binding site of the OF system. The formation of the energy barrier was found to be due to the conformational changes in the channel. Analysis of the MD trajectories revealed that the reorganization of the side chains of the Tyr244 and CYO327 residues play a critical role in the OF channel blocking. Indeed, the calculated distance between Tyr244 and either Cys327 or CYO327 showed a narrowing of the channel after oxidation. The obtained free energy barrier for CYC translocation was found to be 33.9kJmol−1, indicating that oxidation of Cys327, by e.g., cold atmospheric plasma, is more effective in inhibiting the xC− antiporter than in the mutation of this amino acid to Ala (yielding a barrier of 32.4kJmol−1). The inhibition of the xC− antiporter may lead to Cys starvation in some cancer cells, eventually resulting in cancer cell death.

Abstract: n/a

Abstract: Plasma degradation is one of the most effective methods for the abatement of greenhouse gas sulfur hexafluoride
(SF6). To evaluate the potential of MgO as a catalyst in plasma degradation, we investigate the catalytic properties
of MgO on SF6 adsorption and activation by density functional theory (DFT) where the O-defective and
hydroxylated surfaces are considered as two typical plasma-generated surfaces. Our results show that perfect
MgO (001) and (111) surfaces cannot interact with SF6 and only physical adsorption happens. In case of Odefective
MgO surfaces, the O vacancy is the most stable adsorption site. SF6 undergoes a decomposition to SF5
and F over the O-defective MgO (001) surface and undergoes an elongation of the bottom S-F bond over the Odefective
(111) surface. Besides, SF6 shows a physically adsorption at the stepsite of the MgO (001) surface,
accompanied by small changes in its bond angle and length. Furthermore, SF6 is found to be physically and
chemically adsorbed over 0.5 and 1.0 ML (monolayer) H-covered O-terminated MgO (111) surfaces, respectively.
The SF6 molecule undergoes a self-decomposition on the 1.0 ML hydroxylated surface via a surface bonding
process. This study shows that defective and hydroxylated MgO surfaces have the surface capacities for SF6
activation, which shows that MgO has potential as packing material in SF6 waste treatment in packed-bed
plasmas.

Abstract: Non-small cell lung cancer (NSCLC) is the most common type of lung cancer, with 85% of all lung cancer reported as NSCLC. Moreover, there are no effective treatments in advanced NSCLC. This study shows for the first time that oral administration of plasma-treated water (PTW) can cure advanced NSCLC. The cold plasma in water generates a cocktail of reactive species, and oral administration of this cocktail to mice showed no toxicities even at the highest dose of PTW, after a single dose and repeated doses for 28 d in mice. In vivo studies reveal that PTW showed favorable anticancer effects on chemo-resistant lung cancer, similarly to gefitinib treatment as a reference drug in a chemo-resistant NSCLC model. The anticancer activities of PTW seem to be involved in inhibiting proliferation and angiogenesis and enhancing apoptosis in the cancer cells. Interestingly, the PTW contributes to enhanced immune response and improved cachexia in the model.

Abstract: This study hypothesizes that the application of low-dose nonthermal biocompatible dielectric barrier discharge plasma (DBD-NBP) to human gingival fibroblasts (HGFs) will inhibit colony formation but not cell death and induce matrix metalloproteinase (MMP) expression, extracellular matrix (ECM) degradation, and subsequent cell migration, which can result in enhanced wound healing. HGFs treated with plasma for 3 min migrate to each other across the gap faster than those in the control and 5-min treatment groups on days 1 and 3. The plasma-treated HGFs show significantly high expression levels of the cell cycle arrest-related p21 gene and enhanced MMP activity. Focal adhesion kinase (FAK) mediated attenuation of wound healing or actin cytoskeleton rearrangement, and plasma-mediated reversal of this attenuation support the migratory effect of DBD-NBP. Further, this work performs computer simulations to investigate the effect of oxidation on the stability and conformation of the catalytic kinase domain (KD) of FAK. It is found that the oxidation of highly reactive amino acids (AAs) Cys427, Met442, Cys559, Met571, Met617, and Met643 changes the conformation and increases the structural flexibility of the FAK protein and thus modulates its function and activity. Low-dose DBD-NBP-induces host cell cycle arrest, ECM breakdown, and subsequent migration, thus contributing to the enhanced wound healing process.

Abstract: Molecular dynamics simulations are essential for a better understanding of dissociative chemisorption on metal surfaces, which is often the rate-controlling step in heterogeneous and plasma catalysis. The workhorse quasi-classical trajectory approach ubiquitous in molecular dynamics is able to accurately predict reactivity only for high translational and low vibrational energies. In contrast, catalytically relevant conditions generally involve low translational and elevated vibrational energies. Existing quantum dynamics approaches are intractable or approximate as a result of the large number of degrees of freedom present in molecule−metal surface reactions. Here, we extend a ring polymer molecular dynamics approach to fully include, for the first time, the degrees of freedom of a moving metal surface. With this approach, experimental sticking probabilities for the dissociative chemisorption of methane on Pt(111) are reproduced for a large range of translational and vibrational energies by including nuclear quantum effects and employing full-dimensional simulations.

Abstract: Targeting the redox balance of malignant cells via the delivery of high oxidative stress unlocks a potential therapeutic strategy against glioblastoma (GBM). We investigated a novel reactive oxygen species (ROS)-inducing combination treatment strategy, by increasing exogenous ROS via cold atmospheric plasma and inhibiting the endogenous protective antioxidant system via auranofin (AF), a thioredoxin reductase 1 (TrxR) inhibitor. The sequential combination treatment of AF and cold atmospheric plasma-treated PBS (pPBS), or AF and direct plasma application, resulted in a synergistic response in 2D and 3D GBM cell cultures, respectively. Differences in the baseline protein levels related to the antioxidant systems explained the cell-line-dependent sensitivity towards the combination treatment. The highest decrease of TrxR activity and GSH levels was observed after combination treatment of AF and pPBS when compared to AF and pPBS monotherapies. This combination also led to the highest accumulation of intracellular ROS. We confirmed a ROS-mediated response to the combination of AF and pPBS, which was able to induce distinct cell death mechanisms. On the one hand, an increase in caspase-3/7 activity, with an increase in the proportion of annexin V positive cells, indicates the induction of apoptosis in the GBM cells. On the other hand, lipid peroxidation and inhibition of cell death through an iron chelator suggest the involvement of ferroptosis in the GBM cell lines. Both cell death mechanisms induced by the combination of AF and pPBS resulted in a significant increase in danger signals (ecto-calreticulin, ATP and HMGB1) and dendritic cell maturation, indicating a potential increase in immunogenicity, although the phagocytotic capacity of dendritic cells was inhibited by AF. In vivo, sequential combination treatment of AF and cold atmospheric plasma both reduced tumor growth kinetics and prolonged survival in GBM-bearing mice. Thus, our study provides a novel therapeutic strategy for GBM to enhance the efficacy of oxidative stress-inducing therapy through a combination of AF and cold atmospheric plasma.

Abstract: In the last decades, AC powered atmospheric dielectric barrier discharges (DBDs) in air with a liquid electrode have been proposed as a promising plasma technology with versatile applicability in medicine agriculture and water treatment. The fundamental features of the micro-discharge filaments that make up this type of plasma have, however, not been studied yet in sufficient detail. In order to address this need, we investigated a single DBD micro-discharge filament over a water film in a sphere-to-sphere electrode configuration, by means of ICCD imaging and optical emission spectroscopy. When the water film temporarily acts as the cathode, the plasma duration is remarkably long and shows a clear similarity with a resistive barrier discharge, which we attribute to the resistive nature of the water film and the formation of a cathode fall. As another striking difference to DBD with solid electrodes, a constant glow-like plasma is observed at the water surface during the entire duration of the applied voltage cycle, indicating continuous plasma treatment of the liquid. We propose several elementary mechanisms that might underlie the observed unique behavior, based on the specific features of a water electrode.

Abstract: Mcl1 is an important anti-apoptotic member of the Bcl2 family proteins that are upregulated in several cancer malignancies. The canonical binding groove (CBG) located at the surface of Mcl1 exhibits a critical role in binding partners selectively via the BH3-domain of pro-apoptotic Bcl2 family members that trigger the downregulation of Mcl1 function. There are several crystal structures of point-mutated pro-apoptotic Bim peptides in complex with Mcl1. However, the mechanistic effects of such point-mutations towards peptide binding and complex stability still remain unexplored. Here, the effects of the reported point mutations in Bim peptides and their binding mechanisms to Mcl1 were computationally evaluated using atomistic-level steered molecular dynamics (SMD) simulations. A range of external-forces and constant-velocities were applied to the Bim peptides to uncover the mechanistic basis of peptide dissociation from the CBG of Mcl1. Although the peptides showed similarities in their dissociation pathways, the peak rupture forces varied significantly. According to simulations results, the disruption of the conserved polar contacts at the complex interface causes a sequential release of the peptides from the CBG of Mcl1. Overall, the results obtained from the current study may provide valuable insights for the development of novel anti-cancer peptide-inhibitors that can downregulate Mcl1’s function.

Abstract: Based on the current knowledge, a plasma–liquid interface looks and behaves very differently from its counterpart at a solid surface. Local processes characteristic to most liquids include a stronger evaporation, surface deformations, droplet ejection, possibly distinct mechanisms behind secondary electron emission, the formation of an electric double layer, and an ion drift-mediated liquid resistivity. All of them can strongly influence the interfacial charge distribution. Accordingly, the plasma sheath at a liquid surface is most likely unique in its own way, both with respect to its structure and behavior. However, insights into these properties are still rather scarce or uncertain, and more studies are required to further disclose them. In this Perspective, we argue why more research on the plasma sheath is not only recommended but also crucial to an accurate understanding of the plasma–liquid interaction. First, we analyze how the sheath regulates various elementary processes at the plasma–liquid interface, in terms of the electrical coupling, the bidirectional mass transport, and the chemistry between plasma and liquid phase. Next, these three regulatory functions of the sheath are illustrated for concrete applications. Regarding the electrical coupling, a great deal of attention is paid to the penetration of fields into biological systems due to their relevance for plasma medicine, plasma agriculture, and food processing. Furthermore, we illuminate the role of the sheath in nuclear fusion, nanomaterial synthesis, and chemical applications. As such, we hope to motivate the plasma community for more fundamental research on plasma sheaths at liquid surfaces.

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.