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: Ammonia is a crucial nutrient used for plant growth and as a building block in pharmaceutical and chemical industry, produced via nitrogen fixation of the ubiquitous atmospheric N2. Current industrial ammonia production relies heavily on fossil resources, but a lot of work is put into developing non-fossil based pathways. Among these is the use of nonequilibrium plasma. In this work, we investigated water vapor as H source for nitrogen fixation into NH3 by non-equilibrium plasma. The highest selectivity towards NH3 was observed with low amounts of added H2O vapor, but the highest production rate was reached at high H2O vapor.

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: We present a zero-dimensional plasma kinetics model, including both surface and gas phase kinetics, to determine the role of vibrationally excited states in plasma-catalytic ammonia synthesis. We defined a new method to systematically capture the conditions of dielectric barrier discharges (DBDs), including those found in packed bed DBDs. We included the spatial and temporal nature of such discharges by special consideration of the number of micro-discharges in the model. We introduce a parameter that assigns only a part of the plasma power to the microdischarges, to scale the model conditions from filamentary to uniform plasma. Because of the spatial and temporal behaviour of the micro-discharges, not all micro-discharges occurring in the plasma reactor during a certain gas residence time are affecting the molecules. The fraction of power considered in the model ranges from 0.005 %, for filamentary plasma, to 100 %, for uniform plasma. If vibrational excitation is included in the plasma chemistry, these different conditions, however, yield an ammonia density that is only varying within one order of magnitude. At only 0.05 % of the power put into the uniform plasma component, a model neglecting vibrational excitation clearly does not result in adequate amounts of ammonia. Thus, our new model, which accounts for the concept in which not all the power is deposited by the micro-discharges, but some part may also be distributed in between them, suggests that vibrational kinetic processes are really important in (packed bed) DBDs. Indeed, vibrational excitation takes place in both the uniform plasma between the micro-discharges and in the strong micro-discharges, and is responsible for an increased N2 dissociation rate. This is shown here for plasma-catalytic ammonia synthesis, but might also be valid for other gas conversion applications.

Abstract: In this work, we use density functional theory calculations to study the combined effect of external electric fields, surface morphology, and surface charge on CO2 activation over Cu(111), Cu(211), Cu(110), and Cu(001) surfaces. We observe that the binding energy of the CO2 molecule on Cu surfaces increases significantly upon increasing the applied electric field strength. In addition, rougher surfaces respond more effectively to the presence of the external electric field toward facilitating the formation of a carbonate-like CO2 structure and the transformation of the most stable adsorption mode from physisorption to chemisorption. The presence of surface charges further strengthens the electric field effect and consequently causes an improved bending of the CO2 molecule and C−O bond length elongation. On the other hand, a net charge in the absence of an externally applied electric field shows only a marginal effect on CO2 binding. The chemisorbed CO2 is more stable and further activated when the effects of an external electric field, rough surface, and surface charge are combined. These results can help to elucidate the underlying factors that control CO2 activation in heterogeneous and plasma catalysis, as well as in electrochemical processes.

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: 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: 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: 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.

Abstract: Cancer is the second leading cause of death worldwide, and while science has advanced significantly to improve the treatment outcome and quality of life in cancer patients, there are still many issues with the current therapies, such as toxicity and the development of resistance to treatment [...]

Abstract: NADPH oxidases 1 (NOX1) derived reactive oxygen species (ROS) play an important role in the progression of cancer through signaling pathways. Therefore, in this paper, we demonstrate the effect of cold atmospheric plasma (CAP) on the structural changes of Noxa1 SH3 protein, one of the regulatory subunits of NOX1. For this purpose, firstly we purified the Noxa1 SH3 protein and analyzed the structure using X-ray crystallography, and subsequently, we treated the protein with two types of CAP reactors such as pulsed dielectric barrier discharge (DBD) and Soft Jet for different time intervals. The structural deformation of Noxa1 SH3 protein was analyzed by various experimental methods (circular dichroism, fluorescence, and NMR spectroscopy) and by MD simulations. Additionally, we demonstrate the effect of CAP (DBD and Soft Jet) on the viability and expression of NOX1 in A375 cancer cells. Our results are useful to understand the structural modification/oxidation occur in protein due to reactive oxygen and nitrogen (RONS) species generated by CAP.

Abstract: The requirements for new technologies to serve as anticancer agents go far beyond their toxicity potential. Novel applications also need to be safe on a molecular and patient level. In a broader sense, this also relates to cancer metastasis and inflammation. In a previous study, the toxicity of an atmospheric pressure argon plasma jet in four human pancreatic cancer cell lines was confirmed and plasma treatment did not promote metastasis in vitro and in ovo. Here, these results are extended by additional types of analysis and new models to validate and define on a molecular level the changes related to metastatic processes in pancreatic cancer cells following plasma treatment in vitro and in ovo. In solid tumors that were grown on the chorion-allantois membrane of fertilized chicken eggs (TUM-CAM), plasma treatment induced modest to profound apoptosis in the tissues. This, however, was not associated with a change in the expression levels of adhesion molecules, as shown using immunofluorescence of ultrathin tissue sections. Culturing of the cells detached from these solid tumors for 6d revealed a similar or smaller total growth area and expression of ZEB1, a transcription factor associated with cancer metastasis, in the plasma-treated pancreatic cancer tissues. Analysis of in vitro and in ovo supernatants of 13 different cytokines and chemokines revealed cell line-specific effects of the plasma treatment but a noticeable increase of, e.g., growth-promoting interleukin 10 was not observed. Moreover, markers of epithelial-to-mesenchymal transition (EMT), a metastasis-promoting cellular program, were investigated. Plasma-treated pancreatic cancer cells did not present an EMT-profile. Finally, a realistic 3D tumor spheroid co-culture model with pancreatic stellate cells was employed, and the invasive properties in a gel-like cellular matrix were investigated. Tumor outgrowth and spread was similar or decreased in the plasma conditions. Altogether, these results provide valuable insights into the effect of plasma treatment on metastasis-related properties of cancer cells and did not suggest EMT-promoting effects of this novel cancer therapy.

Abstract: Ammonia is one of the most produced chemicals, mainly synthesized from fossil fuels for fertilizer applications. Furthermore, ammonia may be one of the energy carriers of the future, when it is produced from renewable electricity. This has spurred research on alternative technologies for green ammonia production. Research on plasma-driven ammonia synthesis has recently gained traction in academic literature. In the current review, we summarize the literature on plasma-driven ammonia synthesis. We distinguish between mechanisms for ammonia synthesis in the presence of a plasma, with and without a catalyst, for different plasma conditions. Strategies for catalyst design are discussed, as well as the current understanding regarding the potential plasma-catalyst synergies as function of the plasma conditions and their implications on energy efficiency. Finally, we discuss the limitations in currently reported models and experiments, as an outlook for research opportunities for further unravelling the complexities of plasma-catalytic ammonia synthesis, in order to bridge the gap between the currently reported models and experimental results.

Abstract: Plasma has been studied for several years to convert CO2 into value-added products. If CO2 could be converted in the presence of H2O as a cheap H-source for making syngas and oxygenates, it would mimic natural photosynthesis. However, CO2/H2O plasmas have not yet been extensively studied, not by experiments, and certainly not computationally. Therefore, we present here a kinetic modelling study to obtain a greater understanding of the vibrational kinetics of a CO2/H2O microwave plasma. For this purpose, we first created an electron impact cross section set for H2O, using a swarm-derived method. We added the new cross section set and CO2/H2O-related chemistry to a pure CO2 model. While it was expected that H2O addition mainly causes quenching of the CO2 asymmetric mode vibrational levels due to the additional CO2/H2O vibrational-translational relaxation, our model shows that the modifications in the vibrational kinetics are mainly induced by the strong electron dissociative attachment to H2O molecules, causing a reduction in electron density, and the corresponding changes in the input of energy into the CO2 vibrational levels by electron impact processes.