Abstract: Xenobiotic nucleic acids (XNA) are nucleic acid analogues not present in nature that can be used for the storage of genetic information. In vivo XNA applications could be developed into novel biocontainment strategies, but are currently limited by the challenge of developing XNA processing enzymes such as polymerases, ligases and nucleases. Here, we present a structure-guided modelling-based strategy for the rational design of those enzymes essential for the development of XNA molecular biology. Docking of protein domains to unbound double-stranded nucleic acids is used to generate a first approximation of the extensive interaction of nucleic acid processing enzymes with their substrate. Molecular dynamics is used to optimise that prediction allowing, for the first time, the accurate prediction of how proteins that form toroidal complexes with nucleic acids interact with their substrate. Using the Chlorella virus DNA ligase as a proof of principle, we recapitulate the ligase's substrate specificity and successfully predict how to convert it into an XNA-templated XNA ligase.

Abstract: Plasma catalysis is promising for various environmental, energy and chemical synthesis applications, but the underlying mechanisms are far from understood. Modeling can help to obtain a better insight in these mechanisms. Some burning questions relate to the plasma behavior inside packed bed reactors and whether plasma can penetrate into catalyst pores. In this paper, we try to provide answers to these questions, by means of both fluid modeling and particle-in-cell/Monte Carlo collision simulations. We present a short overview of recent findings obtained in our group by means of modeling, i.e., the enhanced electric field near the contact points and the streamer propagation through the packing in packed bed reactors, as well as the plasma behavior in catalyst pores, to determine the minimum pore size in which plasma streamers can penetrate.

Abstract: Physical plasmas generate unique mixes of reactive oxygen and nitrogen species (RONS or ROS). Only a bit more than a decade ago, these plasmas, operating at body temperature, started to be considered for medical therapy with considerably little mechanistic redox chemistry or biomedical research existing on that topic at that time. Today, a vast body of evidence is available on physical plasma-derived ROS, from their spatiotemporal resolution in the plasma gas phase to sophisticated chemical and biochemical analysis of these species once dissolved in liquids. Data from<italic>in silico</italic>analysis dissected potential reaction pathways of plasma-derived reactive species with biological membranes, and<italic>in vitro</italic>and<italic>in vivo</italic>experiments in cell and animal disease models identified molecular mechanisms and potential therapeutic benefits of physical plasmas. In 2013, the first medical plasma systems entered the European market as class IIa devices and have proven to be a valuable resource in dermatology, especially for supporting the healing of chronic wounds. The first results in cancer patients treated with plasma are promising, too. Due to the many potentials of this blooming new field ahead, there is a need to highlight the main concepts distilled from plasma research in chemistry and biology that serve as a mechanistic link between plasma physics (how and which plasma-derived ROS are produced) and therapy (what is the medical benefit). This inevitably puts cellular membranes in focus, as these are the natural interphase between ROS produced by plasmas and translation of their chemical reactivity into distinct biological responses.

Abstract: Pancreatic ductal adenocarcinoma (PDAC) is one of the most aggressive cancers with a low response to treatment and a five-year survival rate below 5%. The ineffectiveness of treatment is partly because of an immunosuppressive tumor microenvironment, which comprises tumor-supportive pancreatic stellate cells (PSCs). Therefore, new therapeutic strategies are needed to tackle both the immunosuppressive PSC and pancreatic cancer cells (PCCs). Recently, physical cold atmospheric plasma consisting of reactive oxygen and nitrogen species has emerged as a novel treatment option for cancer. In this study, we investigated the cytotoxicity of plasma-treated phosphate-buffered saline (pPBS) using three PSC lines and four PCC lines and examined the immunogenicity of the induced cell death. We observed a decrease in the viability of PSC and PCC after pPBS treatment, with a higher efficacy in the latter. Two PCC lines expressed and released damage-associated molecular patterns characteristic of the induction of immunogenic cell death (ICD). In addition, pPBS-treated PCC were highly phagocytosed by dendritic cells (DCs), resulting in the maturation of DC. This indicates the high potential of pPBS to trigger ICD. In contrast, pPBS induced no ICD in PSC. In general, pPBS treatment of PCCs and PSCs created a more immunostimulatory secretion profile (higher TNF-α and IFN-γ, lower TGF-β) in coculture with DC. Altogether, these data show that plasma treatment via pPBS has the potential to induce ICD in PCCs and to reduce the immunosuppressive tumor microenvironment created by PSCs. Therefore, these data provide a strong experimental basis for further in vivo validation, which might potentially open the way for more successful combination strategies with immunotherapy for PDAC.

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: 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: The combination of catalysis and nonthermal plasma holds promise for enabling difficult chemical conversions. The possible synergy between both depends strongly on the nature of the reactive plasma species and the catalyst material. In this paper, we show how vibrationally excited species and plasma-generated radicals interact with transition metal catalysts and how changing the catalyst material can improve the conversion rates and product selectivity. We developed a microkinetic model to investigate the impact of vibrational excitations and plasma-generated radicals on the nonoxidative coupling of methane over transition metal surfaces. We predict a significant increase in ethylene formation for vibrationally excited methane. Plasma-generated radicals have a stronger impact on the turnover frequencies with high selectivity toward ethylene on noble catalysts and mixed selectivity on non-noble catalysts. In general, we show how the optimal catalyst material depends on the desired products as well as the plasma conditions.

Abstract: Mcl1 is a primary member of the Bcl–2 family—anti–apoptotic proteins (AAP)—that is overexpressed in several cancer pathologies. The apoptotic regulation is mediated through the binding of pro-apoptotic peptides (PAPs) (e.g., Bak and Bid) at the canonical hydrophobic binding groove (CBG) of Mcl1. Although all PAPs form amphipathic α-helices, their amino acid sequences vary to different degree. This sequence variation exhibits a central role in the binding partner selectivity towards different AAPs. Thus, constructing a novel peptide or small organic molecule with the ability to mimic the natural regulatory process of PAP is essential to inhibit various AAPs. Previously reported experimental binding free energies (BFEs) were utilized in the current investigation aimed to understand the mechanistic basis of different PAPs targeted to mMcl1. Molecular dynamics (MD) simulations used to estimate BFEs between mMcl1—PAP complexes using Molecular Mechanics-Generalized Born Solvent Accessible (MMGBSA) approach with multiple parameters. Predicted BFE values showed an excellent agreement with the experiment (R2 = 0.92). The van–der Waals (ΔGvdw) and electrostatic (ΔGele) energy terms found to be the main energy components that drive heterodimerization of mMcl1—PAP complexes. Finally, the dynamic network analysis predicted the allosteric signal transmission pathway involves more favorable energy contributing residues. In total, the results obtained from the current investigation may provide valuable insights for the synthesis of a novel peptide or small organic inhibitor targeting Mcl1.

Abstract: The combination of catalysis with non-thermal plasma is a promising alternative to thermal catalysis. A dielectric-barrier discharge reactor was used to study plasma-catalytic dry reforming of methane at ambient pressure and temperature and a fixed plasma power of 45 W. The effect of different catalytic packing materials was evaluated in terms of conversion, product selectivity, and energy efficiency. The conversion of CO2 (~22%) and CH4 (~33%) were found to be similar in plasma-only and when introducing packing materials in plasma. The main reason is the shorter residence time of the gas due to packing geometry, when compared at identical flow rates. H2, CO, C2-C4 hydrocarbons, and oxygenates were identified in the product gas. High selectivity towards H2 and CO were found for all catalysts and plasma-only, with a H2/CO molar ratio of ~0.9. The lowest syngas selectivity was obtained with Cu/Al2O3 (~66%), which instead, had the highest alcohol selectivity (~3.6%).

Abstract: This paper gives an overview of the modeling work developed in our group in the last 25 years for various plasmas used in analytical spectrochemistry, i.e., glow discharges (GDs), inductively coupled plasmas (ICPs), and laser ablation (LA) for sample introduction in the ICP and for laser-induced breakdown spectroscopy (LIBS). The modeling approaches are briefly presented, which are different for each case, and some characteristic results are illustrated. These plasmas are used not only in analytical chemistry but also in other applications, and the insights obtained in these other fields were quite helpful for us to develop models for the analytical plasmas. Likewise, there is now a huge interest in plasma–liquid interaction, atmospheric pressure glow discharges (APGDs), and dielectric barrier discharges (DBDs) for environmental, medical, and materials applications of plasmas. The insights obtained in these fields are also very relevant for ambient desorption/ionization sources and for liquid sampling, which are nowadays very popular in analytical chemistry, and they could be very helpful in developing models for these sources as well.

Abstract: Pancreatic ductal adenocarcinoma (PDAC) is a lethal disease with low five-year survival rates of 8% by conventional treatment methods, e.g., chemotherapy, radiotherapy, and surgery. PDAC shows high resistance towards chemo- and radiotherapy and only 15–20% of all patients can have surgery. This disease is predicted to become the third global leading cause of cancer death due to its significant rise in incidence. Therefore, the development of an alternative or combinational method is necessary to improve current approaches. Cold atmospheric plasma (CAP) treatments could offer multiple advantages to this emerging situation. The plasma-derived reactive species can induce oxidative damage and a cascade of intracellular signaling pathways, which could lead to cell death. Previous reports have shown that CAP treatment also influences cells in the tumor microenvironment, such as the pancreatic stellate cells (PSCs). These PSCs, when activated, play a crucial role in the propagation, growth and survival of PDAC tumors. However, the effect of CAP on PSCs is not yet fully understood. This review focuses on the application of CAP for PDAC treatment and the importance of PSCs in the response to treatment.

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.