Abstract: The study presents a 2D fluid plasma model of a gliding arc discharge for dissociation of CO2 which allows its subsequent conversion into value-added chemicals. The model is based on the balance equations of charged and neutral particles, the electron energy balance equation, the gas thermal balance equation and the current continuity equation. By choosing the modeling domain to be the plane perpendicular to the arc current, the numerical calculations are significantly simplified. Thus, the model allows us to explore the influence of the gas instabilities (turbulences) on the energy efficiency of CO2 conversion. This paper presents results for plasma parameters at different values of the effective turbulent thermal conductivity leading to enhanced energy transport.

Abstract: We investigate the performance of an atmospheric pressure glow discharge (APGD) reactor for CO2 conversion in three different configurations, through experiments and simulations. The first (basic) configuration utilizes the well-known pin-to-plate design, which offers a limited conversion. The second configuration improves the reactor performance by employing a vortex-flow generator. The third, “confined” configuration is a complete redesign of the reactor, which encloses the discharge in a limited volume, significantly surpassing the conversion rate of the other two designs. The plasma properties are investigated using an advanced plasma model.

Abstract: Melittin (MEL), a small peptide component of bee venom, has been reported to exhibit anti-cancer effects in vitro and in vivo. However, its clinical applicability is disputed because of its non-specific cytotoxicity and haemolytic activity in high treatment doses. Plasma-treated phosphate buffered saline solution (PT-PBS), a solution rich in reactive oxygen and nitrogen species (RONS) can disrupt the cell membrane integrity and induce cancer cell death through oxidative stress-mediated pathways. Thus, PT-PBS could be used in combination with MEL to facilitate its access into cancer cells and to reduce the required therapeutic dose. The aim of our study is to determine the reduction of the effective dose of MEL required to eliminate cancer cells by its combination with PT-PBS. For this purpose, we have optimised the MEL threshold concentration and tested the combined treatment of MEL and PT-PBS on A375 melanoma and MCF7 breast cancer cells, using in vitro, in ovo and in silico approaches. We investigated the cytotoxic effect of MEL and PT-PBS alone and in combination to reveal their synergistic cytological effects. To support the in vitro and in ovo experiments, we showed by computer simulations that plasma-induced oxidation of the phospholipid bilayer leads to a decrease of the free energy barrier for translocation of MEL in comparison with the non-oxidized bilayer, which also suggests a synergistic effect of MEL with plasma induced oxidation. Overall, our findings suggest that MEL in combination with PT-PBS can be a promising combinational therapy to circumvent the non-specific toxicity of MEL, which may help for clinical applicability in the future.

Abstract: Pesticides are emerging contaminants frequently detected in the aquatic environment. In this work, a novel approach combining activated carbon adsorption, oxygen plasma treatment and ozonation was studied for the removal of the persistent chlorinated pesticide alachlor. A comparison was made between the removal efficiency and energy consumption for two different reactor operation modes: batchrecirculation and single-pass mode. The kinetics study revealed that the insufficient removal of alachlor by adsorption was significantly improved in terms of degradation efficiency and energy consumption when combined with the plasma treatment. The best efficiency (ca. 80% removal with an energy cost of 19.4 kWh mÀ3) was found for the single-pass operational mode of the reactor. In the batch-recirculating process, a complete elimination of alachlor by plasma treatment was observed after 30 min of treatment. Analysis of the reactive species induced by plasma in aqueous solutions showed that the decomposition of alachlor mainly occurred through a radical oxidation mechanism, with a minor contribution of long-living oxidants (O3, H2O2). Investigation of the alachlor oxidation pathways revealed six different oxidation mechanisms, including the loss of aromaticity which was never before reported for plasma-assisted degradation of aromatic pesticides. It was revealed that the removal rate and energy cost could be further improved with more than 50% by additional O3 gas bubbling in the solution reservoir.

Abstract: Cold physical plasma has limited tumor growth in many preclinical models and is, therefore, suggested as a putative therapeutic option against cancer. Yet, studies investigating the cells’ metastatic behavior following plasma treatment are scarce, although being of prime importance to evaluate the safety of this technology. Therefore, we investigated four human pancreatic cancer cell lines for their metastatic behavior in vitro and in chicken embryos (in ovo). Pancreatic cancer was chosen as it is particularly metastatic to the peritoneum and systemically, which is most predictive for outcome. In vitro, treatment with the kINPen plasma jet reduced pancreatic cancer cell activity and viability, along with unchanged or decreased motility. Additionally, the expression of adhesion markers relevant for metastasis was down-regulated, except for increased CD49d. Analysis of 3D tumor spheroid outgrowth showed a lack of plasma-spurred metastatic behavior. Finally, analysis of tumor tissue grown on chicken embryos validated the absence of an increase of metabolically active cells physically or chemically detached with plasma treatment. We conclude that plasma treatment is a safe and promising therapeutic option and that it does not promote metastatic behavior in pancreatic cancer cells in vitro and in ovo.

Abstract: Increasing the selectivity of cancer treatments is attractive, as it has the potential to reduce side-effects of therapy. Cold atmospheric plasma (CAP) is a novel cancer treatment that disrupts the intracellular oxidative balance. Several reports claim CAP treatment to be selective, but retrospective analysis of these studies revealed discrepancies in several biological factors and culturing methods. Before CAP can be conclusively stated as a selective cancer treatment, the importance of these factors must be investigated. In this study, we evaluated the influence of the cell type, cancer type, and cell culture medium on direct and indirect CAP treatment. Comparison of cancerous cells with their non-cancerous counterparts was performed under standardized conditions to determine selectivity of treatment. Analysis of seven human cell lines (cancerous: A549, U87, A375, and Malme-3M; non-cancerous: BEAS-2B, HA, and HEMa) and five different cell culture media (DMEM, RPMI1640, AM, BEGM, and DCBM) revealed that the tested parameters strongly influence indirect CAP treatment, while direct treatment was less affected. Taken together, the results of our study demonstrate that cell type, cancer type, and culturing medium must be taken into account before selectivity of CAP treatment can be claimed and overlooking these parameters can easily result in inaccurate conclusions of selectivity.

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: 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: Tumours are complex systems formed by cellular (malignant, immune, and endothelial cells, fibroblasts) and acellular components (extracellular matrix (ECM) constituents and secreted factors). A close interplay between these factors, collectively called the tumour microenvironment, is required to respond appropriately to external cues and to determine the treatment outcome. Cold plasma (here referred as ‘plasma’) is an emerging anticancer technology that generates a unique cocktail of reactive oxygen and nitrogen species to eliminate cancerous cells via multiple mechanisms of action. While plasma is currently regarded as a local therapy, it can also modulate the mechanisms of cell-to-cell and cell-to-ECM communication, which could facilitate the propagation of its effect in tissue and distant sites. However, it is still largely unknown how the physical interactions occurring between cells and/or the ECM in the tumour microenvironment affect the plasma therapy outcome. In this review, we discuss the effect of plasma on cell-to-cell and cell-to-ECM communication in the context of the tumour microenvironment and suggest new avenues of research to advance our knowledge in the field. Furthermore, we revise the relevant state-of-the-art in three-dimensional in vitro models that could be used to analyse cell-to-cell and cell-to-ECM communication and further strengthen our understanding of the effect of plasma in solid tumours.

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: We present an approach to incorporate the effect of vibrational nonequilibrium in molecular dynamics (MD) simulations. A perturbed canonical ensemble, in which selected modes are excited to higher temperature while all others remain equilibrated at low temperature, is simulated by applying a specifically tailored bias potential. Our method can be readily applied to any (classical or quantum mechanical) MD setup at virtually no additional computational cost and allows the study of reactions of vibrationally excited molecules in nonequilibrium environments such as plasmas. In combination with enhanced sampling methods, the vibrational efficacy and mode selectivity of vibrationally stimulated reactions can then be quantified in terms of chemically relevant observables, such as reaction rates and apparent free energy barriers. We first validate our method for the prototypical hydrogen exchange reaction and then show how it can capture the effect of vibrational excitation on a symmetric SN2 reaction and radical addition on CO2.

Abstract: Understanding the folding and stability of membrane proteins is of great importance in protein science. Recently, osmolytes and ionic liquids (ILs) are increasingly being used as drug delivery systems in the biopharmaceutical industry. However, the stability of membrane proteins in the presence of osmolytes and ILs is not yet fully understood. Besides, the effect of oxidative stress on membrane proteins with osmolytes or ILs has not been investigated. Therefore, we studied the influence of osmolytes and ILs as co-solvents on the stability of a model membrane protein (i.e., Bacteriorhodopsin in purple membrane of Halobacterium salinarum), using UV–Vis spectroscopy and molecular dynamics (MD) simulations. The MD simulations allowed us to determine the flexibility and solvent accessible surface area (SASA) of Bacteriorhodopsin protein in the presence and/or absence of cosolvents, as well as to carry out principal component analysis (PCA) to identify the most important movements in this protein. In addition, by means of UV–Vis spectroscopy we studied the effect of oxidative stress generated by cold atmospheric plasma on the stability of Bacteriorhodopsin in the presence and/or absence of co-solvents. This study is important for a better understanding of the stability of proteins in the presence of oxidative stress.

Abstract: A gliding arc (GA) plasma has great potential for CO2 conversion into value-added chemicals, because of its high energy efficiency. To improve the application, a 2D/3D fluid model is needed to investigate the CO2 conversion mechanisms in the actual discharge geometry. Therefore, the complex CO2 chemical kinetics description must be reduced due to the huge computational cost associated with 2D/3D models. This paper presents a chemistry reduction method for CO2 plasmas, based on the so-called directed relation graph method. Depending on the defined threshold values, some marginal species are identified. By means of a sensitivity analysis, we can further reduce the chemistry set by removing one by one the marginal species. Based on the socalled flux-sensitivity coupling, we obtain a reduced CO2 kinetics model, consisting of 36 or 15 species (depending on whether the 21 asymmetric mode vibrational states of CO2 are explicitly included or lumped into one group), which is applied to a GA discharge. The results are compared with those predicted with the full chemistry set, and very good agreement is reached. Moreover, the range of validity of the reduced CO2 chemistry set is checked, telling us that this reduced set is suitable for low power GA discharges. Finally, the time and spatial evolution of the CO2 plasma characteristics are presented, based on a 2D model with the reduced kinetics.

Abstract: Nitrogen oxyanions and oxyacids are important agents in atmospheric chemistry and medical biology. Although their chemical behavior in solution is relatively well understood, they may behave very differently at the water/air interface of atmospheric aerosols or at the membrane/water interface of cells. Here, we developed a fully classical model for molecular dynamics simulations of NO3−, NO2−, HNO3, and HNO2 in the framework of the GROMOS 53A6 and 54A7 force field versions. The model successfully accounted for the poorly structured solvation shell and ion pairing tendency of NO3−. Accurate pure-liquid properties and hydration free energies were obtained for the oxyacids. Simulations at the water/air interface showed a local enrichment of HNO3 and depletion of NO3−. The effect was discussed in light of earlier spectroscopic data and ab initio calculations, suggesting that HNO3 behaves as a weaker acid at the surface of water. Our model will hopefully allow for efficient and accurate simulations of nitrogen oxyanions and oxyacids in solution and at microheterogeneous interface environments.

Abstract: Plasma is gaining interest for CH4 conversion into higher hydrocarbons and H2. However, the performance in terms of conversion and selectivity toward different hydrocarbons is different for different plasma types, and the underlying mechanisms are not yet fully understood. Therefore, we study here these mechanisms in different plasma sources, by means of a chemical kinetics model. The model is first validated by comparing the calculated conversions and hydrocarbon/H2 selectivities with experimental results in these different plasma types and over a wide range of specific energy input (SEI) values. Our model predicts that vibrational−translational nonequilibrium is negligible in all CH4 plasmas investigated, and instead, thermal conversion is important. Higher gas temperatures also lead to a more selective production of unsaturated hydrocarbons (mainly C2H2) due to neutral dissociation of CH4 and subsequent dehydrogenation processes, while three-body recombination reactions into saturated hydrocarbons (mainly C2H6, but also higher hydrocarbons) are dominant in low temperature plasmas.