Abstract: Aquaporins (AQPs) are transmembrane proteins that conduct not only water molecules across the cell membrane but also other solutes, such as reactive oxygen and nitrogen species (RONS), produced (among others) by cold atmospheric plasma (CAP). These RONS may induce oxidative stress in the cell interior, which plays a role in cancer treatment. The underlying mechanisms of the transport of RONS across AQPs, however, still remain obscure. We apply molecular dynamics simulations to investigate the permeation of both hydrophilic (H₂O₂and OH) and hydrophobic (NO₂and NO) RONS through AQP1. Our simulations show that these RONS can all penetrate across the pores of AQP1. The permeation free energy barrier of OH and NO is lower than that of H₂O₂and NO₂, indicating that these radicals may have easier access to the pore interior and interact with the amino acid residues of AQP1. We also study the effect of RONS-induced oxidation of both the phospholipids and AQP1 (i.e., sulfenylation of Cys₁₉₁) on the transport of the above-mentioned RONS across AQP1. Both lipid and protein oxidation seem to slightly increase the free energy barrier for H₂O₂and NO₂permeation, while for OH and NO, we do not observe a strong effect of oxidation. The simulation results help to gain insight in the underlying mechanisms of the noticeable rise of CAP-induced RONS in cancer cells, thereby improving our understanding on the role of AQPs in the selective anticancer capacity of CAP.

Abstract: Background: Strong electric fields are knownto affect cell membrane permeability,which can be applied for therapeutic purposes, e.g., in cancer therapy. A synergistic enhancement of this effect may be accomplished by the presence of reactive oxygen species (ROS), as generated in cold atmospheric plasmas. Little is known about the synergy between lipid oxidation by ROS and the electric field, nor on howthis affects the cell membrane permeability.
Method: We here conduct molecular dynamics simulations to elucidate the dynamics of the permeation process under the influence of combined lipid oxidation and electroporation. A phospholipid bilayer (PLB), consisting of di-oleoyl-phosphatidylcholine molecules covered with water layers, is used as a model system for the plasma membrane.
Results and conclusions:Weshow howoxidation of the lipids in the PLB leads to an increase of the permeability of the bilayer to ROS, although the permeation free energy barriers still remain relatively high. More importantly, oxidation of the lipids results in a drop of the electric field threshold needed for pore formation (i.e., electroporation) in the PLB. The created pores in the membrane facilitate the penetration of reactive plasma species deep into the cell interior, eventually causing oxidative damage.
General significance: This study is of particular interest for plasma medicine, as plasma generates both ROS and electric fields, but it is also of more general interest for applications where strong electric fields and ROS both come into play.

Abstract: The conformational and configurational preferences of Me2PNSNPMe2 (3) and Me2PNSNAsMe2 (4) have been identified using quantum chemical calculations at the DFT/B3LYP/6-311+G* level of theory. An approach in which energetic, structural (geometries and bond orders), electronic (analysis of the electron density) and spectroscopic properties are combined leads to the conclusion that these sulfur-nitrogen-pnictogen chains share many of the properties of their chalcogen-nitrogen analogues but that the through-space intramolecular interactions favouring the Z,Z configuration are even weaker than in these latter compounds. The results of this analysis also lead to an unambiguous assignment of the variable-temperature 31P and 15N NMR spectra of these compounds and their structures both in solution and in the solid state.

Abstract: Intercellular communication plays a crucial role in cancer, as well as other diseases, such as inflammation, tissue degeneration, and neurological disorders. One of the proteins responsible for this, are connexins (Cxs), which come together to form a hemichannel. When two hemichannels of opposite cells interact with each other, they form a gap junction (GJ) channel, connecting the intracellular space of these cells. They allow the passage of ions, reactive oxygen and nitrogen species (RONS), and signaling molecules from the interior of one cell to another cell, thus playing an essential role in cell growth, differentiation, and homeostasis. The importance of GJs for disease induction and therapy development is becoming more appreciated, especially in the context of oncology. Studies have shown that one of the mechanisms to control the formation and disruption of GJs is mediated by lipid oxidation pathways, but the underlying mechanisms are not well understood. In this study, we performed atomistic molecular dynamics simulations to evaluate how lipid oxidation influences the channel properties of Cx26 hemichannels, such as channel gating and permeability. Our results demonstrate that the Cx26 hemichannel is more compact in the presence of oxidized lipids, decreasing its pore diameter at the extracellular side and increasing it at the amino terminus domains, respectively. The permeability of the Cx26 hemichannel for water and RONS molecules is higher in the presence of oxidized lipids. The latter may facilitate the intracellular accumulation of RONS, possibly increasing oxidative stress in cells. A better understanding of this process will help to enhance the efficacy of oxidative stress-based cancer treatments.

Abstract: Biological membranes are under constant attack of free radicals, which may lead to lipid nitro-oxidation, pro­ ducing a complex mixture of nitro-oxidized lipids that are responsible for structural and dynamic changes on the membrane. Despite the latter, nitro-oxidized lipids are also associated with several inflammatory and neuro­ degenerative diseases, the underlying mechanisms of which remain elusive. We perform atomistic molecular dynamics simulations using several isomers of nitro-oxidized lipids to study their effect on the structure and permeability of the membrane, as well as the interaction between the mixture of these products in the phos­pholipid membrane environment. Our results show that the stereo- and positional isomers have a stronger effect on the properties of the membrane composed of oxidized lipids compared to that containing nitrated lipids. Nevertheless, nitrated lipids lead to three-fold increase in water permeability compared to oxidized lipids. In addition, we show that in a membrane consisting of combined nitro-oxidized lipid products, the presence of oxidized lipids protects the membrane from transient pores. Is well stablished that plasma application and photodynamic therapy produces a number of oxidative species used to kill cancer cells, through membrane damage induced by nitro-oxidative stress. This study is important to elucidate the mechanisms and the molecular level properties involving the reactive species produced during that cancer therapies.

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: 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: In this work, the effects of N-doping into the Co-doped single vacancy (Co-SV-G) and di-vacancy graphene flake (Co-dV-G) are investigated and compared toward direct oxidation of methane to methanol (DOMM) employing two different oxidants (N2O and O-2) using density functional theory (DFT) calculation. We found that DOMM on CoN3-G utilizing the N2O molecule as oxygen-donor proceeds via a two-step reaction with low activation energies. In addition, we found that although CoN3-G might be a good catalyst for methane conversion, it can also catalyze the oxidation of methanol to CO2 and H2O due to the required low activation barriers. Moreover, the adsorption behaviors of CHx (x = 0-4) species and dehydrogenation of CHx (x = 1-4) species on CoN3-G are investigated. We concluded that CoN3-G can be used as an efficient catalyst for DOMM and N-2O reduction at ambient conditions which may serve as a guide for fabricating effective C/N catalysts in energy-related devices.

Abstract: A seriesof [TiN/TiAlN]nmultilayer coatingswith different bilayer numbers n=5, 10, 25, 50, and 100 were deposited on stainless steel substrate AISI 304 by a lateral rotating cathode arc technique in a flowing nitrogen atmosphere. The composition and microstructure of the coatings have been analyzed by using energy dispersive X-ray spectroscopy, X-ray diffraction (XRD), and conventional and high-resolution transmission electron microscopy (HRTEM). XRD analysis shows that the preferential orientation growth along the (111) direction is reduced in the multilayer coatings. TEM analysis reveals that the grain size of the coatings decreases with increasing bilayer number. HRTEMimaging of the multilayer coatings shows a high density misfit dislocation between the TiN and TiAlN layers. The cross-plane thermal conductivity of the coatings was measured by a pulsed photothermal reflectance technique. With increasing bilayer number, the multilayer coatings' thermal conductivity decreases gradually. This reduction of thermal conductivity can be ascribed to increased phonon scattering due to the disruption of columnar structure, reduced preferential orientation, decreased grain size of the coatings and present misfit dislocations at the interfaces.

Abstract: The compositional influence of Cr and Y on the microstructure of Mg―Cr―O, and Mg―Y―O films synthesized by reactive magnetron sputtering has been investigated by transmission electron microscopy, X-ray diffraction and molecular dynamics simulations. A decrease in crystallinity is observed in these films as the M (Cr or Y) content is increased. It is found that M forms a solid solution with MgO for metal ratios up to ~ 70% and ~ 50% for Cr and Y respectively. Above ~ 70% Cr metal ratio the Mg―Cr―O films are found to be completely amorphous. The Mg―Y―O films are composed of Mg(Y)O and Y2O3 nano crystallites, up to ~ 50% Y metal ratio. Above this ratio, only Y2O3 nano crystallites are found. The preferential < 111> MgO grain alignment is strongly affected by the increase in M content. For M metal ratios up to ~ 50%, there is a selective promotion of the < 100> MgO grain alignments and a decline in the < 111> grain alignments.

Abstract: Understanding the fundamental dynamics in carbon nanotube (CNT) catalysts is of primary importance to understand CNT nucleation. This Letter reports on calculated bond switching (BS) rates in pure and carbon containing nickel nanoclusters. The rates are analyzed in terms of their temperature dependent spatial distribution and the mobility of the cluster atoms. The BS mechanism is found to change from vibrational to diffusional at around 900 K, with a corresponding strong increase in activation energy. Furthermore, the BS activation energy is observed to decrease as the carbon content in the cluster increases, resulting in an effective liquification of the cluster.

Abstract: Warm plasma offers a promising route for CO2 splitting into valuable CO, yet recombination reactions of CO with oxygen, forming again CO2, have recently emerged as critical limitation. This study combines experiments and fluid dynamics + chemical kinetics modelling to comprehensively analyse the recombination reactions upon CO2 splitting in an atmospheric plasmatron. We introduce an innovative in-situ gas sampling technique, enabling 2D spatial mapping of gas product compositions and temperatures, experimentally confirming for the first time the substantial limiting effect of CO recombination reactions in the afterglow region. Our results show that the CO mole fraction at a 5 L/min flow rate drops significantly from 11.9 % at a vertical distance of z = 20 mm in the afterglow region to 8.6 % at z = 40 mm. We constructed a comprehensive 2D model that allows for spatial reaction rates analysis incorporating crucial reactions, and we validated it to kinetically elucidate this phenomenon. CO2 +M⇌O+CO+M and CO2 +O⇌CO+O2 are the dominant reactions, with the forward reactions prevailing in the plasma region and the backward reactions becoming prominent in the afterglow region. These results allow us to propose an afterglow quenching strategy for performance enhancement, which is further demonstrated through a meticulously developed plasmatron reactor with two-stage cooling. Our approach substantially increases the CO2 conversion (e.g., from 6.6 % to 19.5 % at 3 L/min flow rate) and energy efficiency (from 13.5 % to 28.5 %, again at 3 L/min) and significantly shortens the startup time (from ~ 150 s to 25 s). Our study underscores the critical role of inhibiting recombination reactions in plasma-based CO2 conversion and offers new avenues for performance enhancement.

Abstract: By combining dielectric barrier discharge plasma and external heating, we exploit a two-stage hybrid plasmathermal
system (HPTS), i.e., a plasma stage followed by a thermal stage, for direct non-oxidative coupling of
CH4 to C2H4 and H2, yielding a CH4 conversion of ca. 17 %. In the two-stage HPTS, the plasma first converts CH4
into C2H6 and C3H8, which in the thermal stage leads to a high C2H4 selectivity of ca. 63 % by pyrolysis, with H2
selectivity of ca. 64 %.

Abstract: Given the ecological problems associated to the CO2 emissions of fossil fuels, plasma technology has gained
interest for conversion of CO2 into value-added products. Microwave plasmas operating at atmospheric pressure
have proven to be especially interesting, due to the high gas temperatures inside the reactor (i.e. up to 6000 K)
allowing for efficient thermal dissociation of CO2 into CO and O2. However, the performance of these high
temperature plasmas is limited by recombination of CO back into CO2 once the gas cools down in the afterglow.
In this work, we computationally investigated several quenching nozzles, developed and experimentally tested
by Hecimovic et al., [1] for their ability to quickly cool the gas after the plasma, thereby quenching the CO
recombination reactions. Using a 3D computational fluid dynamics model and a quasi-1D chemical kinetics
model, we reveal that a reactor without nozzle lacks gas mixing between hot gas in the center and cold gas near
the reactor walls. Especially at low flow rates, where there is an inherent lack of convective cooling due to the
low gas flow velocity, the temperature in the afterglow remains high (between 2000 and 3000 K) for a relatively
long time (in the 0.1 s range). As shown by our quasi-1D chemical kinetics model, this results in a important loss
of CO due to recombination reactions. Attaching a nozzle in the effluent of the reactor induces fast gas quenching
right after the plasma. Indeed, it introduces (i) more convective cooling by forcing cool gas near the walls to mix
with hot gas in the center of the reactor, as well as (ii) more conductive cooling through the water-cooled walls of
the nozzle. Our model shows that gas quenching and the suppression of recombination reactions have more
impact at low flow rates, where recombination is the most limiting factor in the conversion process.

Abstract: Plasma-catalytic ammonia synthesis in a dielectric barrier discharge reactor has emerged as a possible route for electrification of nitrogen fixation. In this study, we use a combination of experiments and a plasma kinetic model to investigate the ammonia synthesis from N2 and H2, both with and without a solid packing material in the plasma zone. The effect of plasma power, feed flow rate, N2:H2 feed ratio, gas residence time, temperature, and packing material (MgAl2O4 alone or impregnated with Co or Ru) on the ammonia synthesis rate were examined in the experiments. The kinetic model was employed to improve our understanding of the ammonia formation pathways and identify possible changes in these pathways when altering the N2:H2 feed ratio. A higher NH3 synthesis rate was achieved when increasing the feed flow rate, as well as when increasing the gas tem-perature from 100 to 200 ◦C when a packing material was present in the plasma. At the elevated temperature of 200 ◦C, an optimum in the NH3 synthesis rate was observed at an equimolar feed ratio (N2:H2 =1:1) for the plasma alone and MgAl2O4, while a N2-rich feed was favored for Ru/MgAl2O4 and Co/MgAl2O4. The optimum in the synthesis rate with the N2-rich feed, where high energy electrons are more likely to collide with N2, suggests that the rate-limiting step is the dissociation of N2 in the gas phase. This is supported by the kinetic model when packing material was used. However, for the plasma alone, the model found that the N2 dissociation is only rate limiting in H2-rich feeds, whereas the limited access to H in N2-rich feeds makes the hydrogenation of N species limiting.