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.

Abstract: This study focuses on a magnetically stabilized gliding arc (MGA) plasma. Two fully coupled flow-plasma models (in 3D and 2D) are presented. The 3D model is applied to compare the arc dynamics of the MGA with a traditional gas-driven gliding arc. The 2D model is used for a detailed parametric study on the effect of the external magnetic field. The results show that the relative velocity between the plasma and feed gas is generated due to the Lorentz force, which can increase the plasma-treated gas fraction. The magnetic field also helps to decrease the gas temperature by enhancing heat transfer and to increase the electron number density. This work shows the potential of an external magnetic field to control the gliding arc behavior, for enhanced gas conversion at low gas flow rates.

Abstract: This work proposes to use core-shell structured spheres to evaluate whether it allows to individually optimize bulk and surface effects of a packing material, in order to optimize conversion and energy efficiency. Different core-shell materials have been prepared by spray coating, using dense spheres (as core) and powders (as shell) of SiO2, Al2O3, and BaTiO3. The materials are investigated for their performance in CO2 dissociation and compared against a benchmark consisting of a packed-bed reactor with the pure dense spheres, as well as an empty reactor. The results in terms of CO2 conversion and energy efficiency show various interactions between the core and shell material, depending on their combination. Al2O3 was found as the best core material under the applied conditions here, followed by BaTiO3 and SiO2, in agreement with their behaviour for the pure spheres. Applying a thin shell layer on the cores showed equal performance between the different shell materials. Increasing the layer thickness shifts this behaviour, and strong combination effects were observed depending on the specific material. Therefore, this method of core-shell spheres has the potential to allow tuning of the packing properties more closely to the application by designing an optimal combination of core and shell.

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 technology provides a sustainable, fossil-free method for N2 fixation, i.e., the conversion of inert atmospheric N2 into valuable substances, such as NOx or ammonia. In this work, we present a novel gliding arc plasmatron at atmospheric pressure for NOx production at different N2/O2 gas feed ratios, offering a promising NOx yield of 1.5% with an energy cost of 3.6 MJ/mol NOx produced. To explain the underlying mechanisms, we present a chemical kinetics model, validated by experiments, which provides insight into the NOx formation pathways and into the ambivalent role of the vibrational kinetics. This allows us to pinpoint the factors limiting the yield and energy cost, which can help to further improve the process.

Abstract: We explore the consequences of nonthermal plasma-activation on product yields in catalytic ammonia synthesis, a reaction that is equilibrium-limited at elevated temperatures. We employ a minimal microkinetic model that incorporates the influence of plasma-activation on N2 dissociation rates to predict NH3 yields into and across the equilibrium-limited regime. NH3 yields are predicted to exceed bulk thermodynamic equilibrium limits on materials that are thermal-rate-limited by N2 dissociation. In all cases, yields revert to bulk equilibrium at temperatures at which thermal reaction rates exceed plasma-activated ones. Beyond-equilibrium NH3 yields are observed in a packed bed dielectric barrier discharge reactor and exhibit sensitivity to catalytic material choice in a way consistent with model predictions. The approach and results highlight the opportunity to exploit synergies between nonthermal plasmas and catalysts to affect transformations at conditions inaccessible through thermal routes.

Abstract: Plasma technology is gaining increasing interest for CO2 conversion. The gas temperature in (and after) the plasma reactor largely affects the performance. Therefore, we examine the effect of cooling/quenching, during and after the plasma, on the CO2 conversion and energy efficiency, for typical “warm” plasmas, by means of chemical kinetics modeling. For plasmas at low specific energy input (SEI ∼ 0.5 eV/molecule), it is best to quench at the plasma end, while for high-SEI plasmas (SEI ∼ 4 eV/molecule), quenching at maximum conversion is better. For low-SEI plasmas, quenching can even increase the conversion beyond the dissociation in the plasma, known as superideal quenching. To better understand the effects of quenching at different plasma conditions, we study the dissociation and recombination rates, as well as the vibrational distribution functions (VDFs) of CO2, CO, and O2. When a high vibrational−translational (VT) nonequilibrium exists at the moment of quenching, the dissociation and recombination reaction rates both increase. Depending on the conversion degree at the moment of quenching, this can lead to a net increase or decrease of CO2 conversion. In general, however, and certainly for equilibrium plasmas at high temperature, quenching after the plasma helps prevent recombination reactions and clearly enhances the final CO2 conversion. We also investigate the effect of different quenching cooling rates on the CO2 conversion and energy efficiency. Finally, we compare plasma-based conversion to purely thermal conversion. For warm plasmas with typical temperatures of 3000−4000 K, the conversion is roughly thermal.

Abstract: We studied hydrogen sulfide (H2S) decomposition into hydrogen (H2) and sulfur (S2) in a gliding arc plasmatron (GAP) and microwave (MW) plasma by a combination of 0D and 2D models. The conversion, energy efficiency, and plasma distribution are examined for different discharge conditions, and validated with available experiments from literature. Furthermore, a comparison is made between GAP and MW plasma. The GAP operates at atmospheric pressure, while the MW plasma experiments to which comparison is made were performed at reduced pressure. Indeed, the MW discharge region becomes very much contracted near atmospheric pressure, at the conditions under study, as revealed by our 2D model. The models predict that thermal reactions play the most important role in H2S decomposition in both plasma types. The GAP has a higher energy efficiency but lower conversion than the MW plasma at their typical conditions. When compared at the same conversion, the GAP exhibits a higher energy efficiency and lower energy cost than the MW plasma.