Abstract: Reactive oxygen and nitrogen species released by cold physical plasma are being proposed as effectors in various clinical conditions connected to inflammatory processes. As these plasmas can be tailored in a wide range, models to compare and control their biochemical footprint are desired to infer on the molecular mechanisms underlying the observed effects and to enable the discrimination between different plasma sources. Here, an improved model to trace short-lived reactive species is presented. Using FTIR, high-resolution mass spectrometry, and molecular dynamics computational simulation, covalent modifications of cysteine treated with different plasmas were deciphered and the respective product pattern used to generate a fingerprint of each plasma source. Such, our experimental model allows a fast and reliable grading of the chemical potential of plasmas used for medical purposes. Major reaction products were identified to be cysteine sulfonic acid, cystine, and cysteine fragments. Less abundant products, such as oxidized cystine derivatives or S-nitrosylated cysteines, were unique to different plasma sources or operating conditions. The data collected point at hydroxyl radicals, atomic O, and singlet oxygen as major contributing species that enable an impact on cellular thiol groups when applying cold plasma in vitro or in vivo.

Abstract: Cold atmospheric plasma (CAP) has great potential for sterilization in the food industry, by deactivation of thermophilic bacteria, but the underlying mechanisms are largely unknown. Therefore, we investigate here whether CAP is able to denature/modify protein from thermophilic bacteria. We focus on MTH1880 (MTH) from Methanobacterium thermoautotrophicum as model protein, which we treated with dielectric barrier discharge (DBD) plasma operating in air for 10, 15 and 20 mins. We analysed the structural changes of MTH using circular dichroism, fluorescence and NMR spectroscopy, as well as the thermal and chemical denaturation, upon CAP treatment. Additionally, we performed molecular dynamics (MD) simulations to determine the stability, flexibility and solvent accessible surface area (SASA) of both the native and oxidised protein.

Abstract: Because of its unique properties, plasma technology has gained much prominence in the
microelectronics industry. Recently, environmental and energy applications of plasmas have gained a lot of attention. In this area, the focus is on converting CO 2 and reforming hydrocarbons, with the goal of developing an efficient single-step ‘gas-to-liquid’ (GTL) process. Here we show that applying tri-reforming principles to plasma—further called ‘plasma-based multi-reforming’—allows us to better control the plasma chemistry and thus the formed products. To demonstrate this, we used chemical kinetics calculations supported by experiments and reveal that better control of the plasma chemistry can be achieved by adding O 2 or H 2 O to a mixture containing CH 4 and CO 2 (diluted in N 2 ). Moreover, by adding O 2 and H 2 O simultaneously, we can tune the plasma chemistry even further, improving the conversions, thermal efficiency and methanol yield. Unlike thermocatalytic reforming, plasma-based reforming is capable of producing methanol in a single step; and compared with traditional plasma-based dry reforming, plasma-based multi-reforming increases the methanol yield by more than seven times and the thermal efficiency by 49%, as revealed by our model calculations. Thus, we believe that by using plasma-based multi-reforming, ‘gas-to-liquid’ conversion may be made efficient and scalable.

Abstract: A growing body of literature has recognized the non-thermal effect of pulsed microwave radiation (PMR) on bacterial systems. However, its mode of action in deactivating bacteria has not yet been extensively investigated. Nevertheless, it is highly important to advance the applications of PMR from simple to complex biological systems. In this study, we first optimized the conditions of the PMR device and we assessed the results by simulations, using ANSYS HFSS (High Frequency Structure Simulator) and a 3D particle-in-cell code for the electron behavior, to provide a better overview of the bacterial cell exposure to microwave radiation. To determine the sensitivity of PMR,<italic>Escherichia coli</italic> and<italic>Staphylococcus aureus</italic>cultures were exposed to PMR (pulse duration: 60 ns, peak frequency: 3.5 GHz) with power density of 17 kW/cm²at the free space of sample position, which would induce electric field of 8.0 kV/cm inside the PBS solution of falcon tube in this experiment at 25 °C. At various discharges (D) of microwaves, the colony forming unit curves were analyzed. The highest ratios of viable count reductions were observed when the doses were increased from 20D to 80D, which resulted in an approximate 6 log reduction in <italic>E. coli</italic>and 4 log reduction in<italic>S. aureus.</italic>Moreover, scanning electron microscopy also revealed surface damage in both bacterial strains after PMR exposure. The bacterial inactivation was attributed to the deactivation of oxidation-regulating genes and DNA damage.

Abstract: A new, single-step protocol for wrapping individual nanosized β-casein micelles with silica is presented. This biomolecule-friendly synthesis proceeds at low protein concentration at almost neutral pH, and makes use of sodium silicate instead of the common silicon alkoxides. This way, formation of potentially protein-denaturizing alcohols can be avoided. The pH of the citrate-buffered synthesis medium is close to the isoelectric point of β-casein, which favours micelle formation. A limited amount of sodium silicate is added to the protein micelle suspension, to form a thin silica coating around the β-casein micelles. The size distribution of the resulting proteinsilica structures was characterized using DLS and SAXS, as well as 1H NMR DOSY with a dedicated pulsed-field gradient cryo-probehead to cope with the low protein concentration. The degree of silica-condensation was investigated by 29Si MAS NMR, and the nanostructure was revealed by advanced electron microscopy techniques such as ESEM and HAADF-STEM. As indicated by the combined characterization results, a silica shell of 2 nm is formed around individual β-casein micelles giving rise to separate protein coresilica shell nanoparticles of 17 nm diameter. This alcohol-free method at mild temperature and pH is potentially suited for packing protein molecules into bio-compatible silica nanocapsules for a variety of applications in biosensing, therapeutic protein delivery and biocatalysis.

Abstract: Many current anti-cancer therapies rely on increasing the intracellular reactive oxygen and nitrogen species (RONS) contents with the aim to induce irreparable damage, which subsequently results in tumor cell death. A novel tool in cancer therapy is the use of cold atmospheric plasma (CAP), which has been found to be very effective in the treatment of many different cancer cell types in vitro as well as in vivo, mainly through the vast generation of RONS. One of the key determinants of the cell's fate will be the interaction of RONS, generated by CAP, with important proteins, i.e. redox-regulatory proteins. One such protein is cytoglobin (CYGB), a recently discovered globin proposed to be involved in the protection of the cell against oxidative stress. In this study, the effect of plasma-produced RONS on CYGB was investigated through the treatment of CYGB with CAP for different treatment times. Spectroscopic analysis of CYGB showed that although chemical modifications occur, its secondary structure remains intact. Mass spectrometry experiments identified these modifications as oxidations of mainly sulfur-containing and aromatic amino acids. With longer treatment time, the treatment was also found to induce nitration of the heme. Furthermore, the two surface-exposed cysteine residues of CYGB were oxidized upon treatment, leading to the formation of intermolecular disulfide bridges, and potentially also intramolecular disulfide bridges. In addition, molecular dynamics and docking simulations confirmed, and further show, that the formation of an intramolecular disulfide bond, due to oxidative conditions, affects the CYGB 3D structure, thereby opening the access to the heme group, through gate functioning of His117. Altogether, the results obtained in this study (1) show that plasma-produced RONS can extensively oxidize proteins and (2) that the oxidation status of two redox-active cysteines lead to different conformations of CYGB.

Abstract: Melanoma arises from pigment-producing cells called melanocytes located in the basal layers of the epidermis of the skin. Cytoglobin (CYGB) is a ubiquitously expressed hexacoordinated globin that is highly enriched in me­lanocytes and frequently downregulated during melanomagenesis. Previously, we showed that non-thermal plasma (NTP)-produced reactive oxygen and nitrogen species (RONS) lead to the formation of an intra­ molecular disulfide bridge that would allow CYGB to function as a redox-sensitive protein. Here, we investigate the cytotoxic effect of indirect NTP treatment in two melanoma cell lines with divergent endogenous CYGB expression levels, and we explore the role of CYGB in determining treatment outcome. Our findings are consistent with previous studies supporting that NTP cytotoxicity is mediated through the production of RONS and leads to apoptotic cell death in melanoma cells. Furthermore, we show that NTP-treated solutions elicit an antioxidant response through the activation of nuclear factor erythroid 2–related factor 2 (NRF2). The knock­ down and overexpression of CYGB respectively sensitizes and protects melanoma cells from RONS-induced apoptotic cell death. The presence of CYGB enhances heme-oxygenase 1 (HO-1) and NRF2 protein expression levels, whereas the absence impairs their expression. Moreover, analysis of the CYGB-dependent transcriptome demonstrates the tumor suppressor long non-coding RNA maternally expressed 3 (MEG3) as a hitherto unde­ scribed link between CYGB and NRF2. Thus, the presence of CYGB, at least in melanoma cells, seems to play a central role in determining the therapeutic outcome of RONS-inducing anticancer therapies, like NTP-treated solutions, possessing both tumor-suppressive and oncogenic features. Hence, CYGB expression could be of in­ terest either as a biomarker or as a candidate for future targeted therapies in melanoma.

Abstract: In this paper, an overview is given of modeling activities going on in our research group, for describing the plasma chemistry and plasmasurface interactions in reactive plasmas. The plasma chemistry is calculated by a fluid approach or by hybrid Monte Carlo (MC)fluid modeling. An example of both is illustrated in the first part of the paper. The example of fluid modeling is given for a dielectric barrier discharge (DBD) in CH4/O2, to describe the partial oxidation of CH4 into value-added chemicals. The example of hybrid MCfluid modeling concerns an inductively coupled plasma (ICP) etch reactor in Ar/Cl2/O2, including also the description of the etch process. The second part of the paper deals with the treatment of plasmasurface interactions on the atomic level, with molecular dynamics (MD) simulations or a combination of MD and MC simulations.

Abstract: Background
During normal pregnancy, placental oxidative stress (OS) is present during all three trimesters and is necessary to obtain normal cell function. However, if OS reaches a certain level, pregnancy complications might arise. In preeclampsia (PE), a dangerous pregnancy specific hypertensive disorder, OS induced in the ischemic placenta causes a systemic inflammatory response and activates maternal endothelial cells. In this study, we aimed to quantify superoxide concentrations (as a measure of systemic OS) using electron paramagnetic resonance (EPR) and correlate them to markers of systemic inflammation, iron status and vascular function.
Methods
Fifty-nine women with a healthy pregnancy (HP), 10 non-pregnant controls (NP) and 28 PE patients (32±3.3weeks) were included. During HP, blood samples for superoxide, neutrophil to lymphocyte ratio (NLR), mean platelet volume (MPV) and iron status were taken at 10, 25 and 39 weeks. Vascular measurements for arterial stiffness (carotid-femoral pulse wave velocity (CF-PWV), augmentation index (AIx), augmentation Pressure (AP)) and microvascular endothelial function (reactive hyperemia index (RHI)) were performed at 35 weeks. In PE, all measurements were performed at diagnosis. CMH (1-hydroxy-3-methoxycarbonyl-2,2,5,5-tetramethylpyrrolidine) was used as spin probe for EPR, since the formed CM radical
corresponds to the amount of superoxide.
Results
Superoxide concentration remains stable during pregnancy (p = 0.92), but is significantly higher compared to the NP controls (p<0.0001). At 25 weeks, there is a significant positive correlation between superoxide and ferritin concentration. (p = 0.04) In PE, superoxide, systemic inflammation and iron status are much higher compared to HP (all p<0.001). During HP, superoxide concentrations correlate significantly with arterial stiffness (all p<0.04), while in PE superoxide is significantly correlated to microvascular endothelial function (p = 0.03).
Conclusions
During HP there is an increased but stable oxidative environment, which is correlated to ferritin concentration. If superoxide levels increase, there is an augmentation in arterial stiffness. In PE pregnancies, systemic inflammation and superoxide concentrations are higher and result in a deterioration of endothelial function. Together, these findings support the hypothesis that vascular function is directly linked to the amount of OS and that measurement of OS in combination with vascular function tests might be used in the prediction of PE.

Abstract: Since the onset of pattern transfer technologies for chip manufacturing, various strategies have been developed to circumvent or overcome aspect ratio dependent etching (ARDE). These methods have, however, their own limitations in terms of etch non-idealities, throughput or costs. Moreover, they have mainly been optimized for individual in-device features and die-scale patterns, while occasionally ending up with poor patterning of metrology marks, affecting the alignment and overlay in lithography. Obtaining a better understanding of the underlying mechanisms of ARDE and how to mitigate them therefore remains a relevant challenge to date, for both marks and advanced nodes. In this work, we accordingly assessed the neutral transport mechanisms in ARDE by means of experiments and multiscale modeling for SiO₂etching with CHF₃/Ar and CF₄/Ar plasmas. The experiments revealed a local maximum in the etch rate for an aspect ratio around unity, i.e. the simultaneous occurrence of regular and inverse reactive ion etching lag for a given etch condition. We were able to reproduce this ARDE trend in the simulations without taking into account charging effects and the polymer layer thickness, suggesting shadowing and diffuse reflection of neutrals as the primary underlying mechanisms. Subsequently, we explored four methods with the simulations to regulate ARDE, by varying the incident plasma species fluxes, the amount of polymer deposition, the ion energy and angular distribution and the initial hardmask sidewall angle, for which the latter was found to be promising in particular. Although our study focusses on feature dimensions characteristic to metrology marks and back-end-of-the-line integration, the obtained insights have a broader relevance, e.g. to the patterning of advanced nodes. Additionally, this work supports the insight that physisorption may be more important in plasma etching at room temperature than originally thought, in line with other recent studies, a topic on which we recommend further research.

Abstract: In this study we report on a 2D fluid model of a gliding arc discharge in argon. Despite the 3D nature of the discharge, 2D models are found to be capable of providing very useful information about the operation of the discharge. We employ two modelsan axisymmetric and a Cartesian one. We show that for the considered experiment and the conditions of a low current arc (around 30 mA) in argon, there is no significant heating of the cathode surface and the discharge is sustained by field electron emission from the cathode accompanied by the formation of a cathode spot. The obtained discharge power and voltage are relatively sensitive to the surface properties and particularly to the surface roughness, causing effectively an amplification of the normal electric field. The arc body and anode region are not influenced by this and depend mainly on the current value. The gliding of the arc is modelled by means of a 2D Cartesian model. The arcelectrode contact points are analysed and the gliding mechanism along the electrode surface is discussed. Following experimental observations, the cathode spot is simulated as jumping from one point to another. A complete arc cycle is modelled from initial ignition to arc decay. The results show that there is no interaction between the successive gliding arcs.

Abstract: The aim of this work consists of the evaluation of atmospheric pressure dielectric barrier discharges for the conversion of greenhouse gases into useful compounds. Therefore, pure CO2 feed flows are administered to the discharge zone at varying discharge frequency, power input, gas temperature and feed flow rates, aiming at the formation of CO and O2. The discharge obtained in CO2 is characterized as a filamentary mode with a microdischarge zone in each half cycle of the applied voltage. It is shown that the most important parameter affecting the CO2-conversion levels is the gas flow rate. At low flow rates, both the conversion and the CO-yield are significantly higher. In addition, also an increase in the gas temperature and the power input give rise to higher conversion levels, although the effect on the CO-yield is limited. The optimum discharge frequency depends on the power input level and it cannot be unambiguously stated that higher frequencies give rise to increased conversion levels. A maximum CO2 conversion of 30% is achieved at a flow rate of 0.05 L min−1, a power density of 14.75 W cm−3 and a frequency of 60 kHz. The most energy efficient conversions are achieved at a flow rate of 0.2 L min−1, a power density of 11 W cm−3 and a discharge frequency of 30 kHz.

Abstract: The chemical complexity of non-equilibrium plasmas poses a challenge for plasma modeling because of the computational load. This paper presents a dimension reduction method for such chemically complex plasmas based on principal component analysis (PCA). PCA is used to identify a low-dimensional manifold in chemical state space that is described by a small number of parameters: the principal components. Reduction is obtained since continuity equations only need to be solved for these principal components and not for all the species. Application of the presented method to a CO2 plasma model including state-to-state vibrational kinetics of CO2 and CO demonstrates the potential of the PCA method for dimension reduction. A manifold described by only two principal components is able to predict the CO2 to CO conversion at varying ionization degrees very accurately.

Abstract: The electron bounce resonance heating (BRH) in dual-frequency capacitively coupled plasmas operated in oxygen is studied by different experimental methods and a particle-in-cell/Monte Carlo collision (PIC/MCC) simulation, and compared with the electropositive argon discharge. In comparison with argon, the experimental results show that in an oxygen discharge the resonance peaks in positive-ion density and light intensity tend to occur at larger electrode gaps. Moreover, at electrode gaps L > 2.5 cm, the positive-ion (and electron) density and the light emission drop monotonically in the oxygen discharge upon increasing L, whereas they rise (after an initial drop) in the argon case. At resonance gap the electronegativity reaches its maximum due to the BRH. All these experimental observations are explained by PIC/MCC simulations, which show that in the oxygen discharge the bulk electric field becomes quite strong and is out of phase with the sheath field. Therefore, it retards the resonance electrons when traversing the bulk, resulting in a suppressed BRH. Both experiment and simulation results show that this effect becomes more pronounced at lower high-frequency power, when the discharge mode changes from electropositive to electronegative. In a pure oxygen discharge, the BRH is suppressed with increasing pressure and almost diminishes at 12 Pa. Finally, the driving frequency significantly affects the BRH, because it determines the phase relation between bulk electric field and sheath electric field.