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: 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: 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: Recent advances in biomedical research in cancer immunotherapy have identified the use of an oxidative stress-based approach to treat cancers, which works by inducing immunogenic cell death (ICD) in cancer cells. Since the anti-cancer effects of non-thermal plasma (NTP) are largely attributed to the reactive oxygen and nitrogen species that are delivered to and generated inside the target cancer cells, it is reasonable to postulate that NTP would be an effective modality for ICD induction. NTP treatment of tumors has been shown to destroy cancer cells rapidly and, under specific treatment regimens, this leads to systemic tumorspecific immunity. The translational benefit of NTP for treatment of cancer relies on its ability to enhance the interactions between NTP-exposed minor cells and local immune cells which initiates subsequent protective immune responses. This review discusses results from recent investigations of NTP application to induce ICD in cancer cells. With further optimization of clinical devices and treatment protocols, NTP can become an essential part of the therapeutic armament against cancer.

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 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 catalysis is gaining increasing interest for various gas conversion applications, such as CO2 conversion into value-added chemicals and fuels, CH4 activation into hydrogen, higher hydrocarbons or oxygenates, and NH3 synthesis. Other applications are already more established, such as for air pollution control, e.g. volatile organic compound remediation, particulate matter and NOx removal. In addition, plasma is also very promising for catalyst synthesis and treatment. Plasma catalysis clearly has benefits over ‘conventional’ catalysis, as outlined in the Introduction. However, a better insight into the underlying physical and chemical processes is crucial. This can be obtained by experiments applying diagnostics, studying both the chemical processes at the catalyst surface and the physicochemical mechanisms of plasma-catalyst interactions, as well as by computer modeling. The key challenge is to design cost-effective, highly active and stable catalysts tailored to the plasma environment. Therefore, insight from thermal catalysis as well as electro- and photocatalysis is crucial. All these aspects are covered in this Roadmap paper, written by specialists in their field, presenting the state-of-the-art, the current and future challenges, as well as the advances in science and technology needed to meet these challenges.

Abstract: Atomistic simulation methods for the quantification of free energies are in wide use. These methods operate by sampling the probability density of a system along a small set of suitable collective variables (CVs), which is, in turn, expressed in the form of a free energy surface (FES). This definition of the FES can capture the relative stability of metastable states but not that of the transition state because the barrier height is not invariant to the choice of CVs. Free energy barriers therefore cannot be consistently computed from the FES. Here, we present a simple approach to calculate the gauge correction necessary to eliminate this inconsistency. Using our procedure, the standard FES as well as its gauge-corrected counterpart can be obtained by reweighing the same simulated trajectory at little additional cost. We apply the method to a number of systems—a particle solvated in a Lennard-Jones fluid, a Diels–Alder reaction, and crystallization of liquid sodium—to demonstrate its ability to produce consistent free energy barriers that correctly capture the kinetics of chemical or physical transformations, and discuss the additional demands it puts on the chosen CVs. Because the FES can be converged at relatively short (sub-ns) time scales, a free energy-based description of reaction kinetics is a particularly attractive option to study chemical processes at more expensive quantum mechanical levels of theory.

Abstract: Antibiotic resistance is one of the world’s most urgent public health problems. Due to its antibacterial properties, cold atmospheric plasma (CAP) may serve as an alternative method to antibiotics. It is claimed that oxidative stress caused by CAP is the main reason of bacteria inactivation. In this work, we computationally investigated the effect of plasma-induced oxidation on various glycolysis metabolites, by monitoring the production of the biomass. We observed that in addition to the significant reduction in biomass production, the rate of some re­actions has increased. These reactions produce anti-oxidant products, showing the bacterial defense mechanism to escape the oxidative damage. Nevertheless, the simulations show that the plasma-induced oxidation effect is much stronger than the defense mechanism, causing killing of the bacteria.

Abstract: Hydroxyapatite (HAP) is frequently used as biocompatible coating on Ti-based implants. In this context, the HAP-Ti adhesion is of crucial importance. Here, we report ab initio calculations to investigate the influence of Si incorporation into the amorphous calcium-phosphate (a-HAP) structure on the interfacial bonding mechanism between the a-HAP coating and an amorphous titanium dioxide (a-TiO2) substrate, contrasting two different density functionals: PBE-GGA, and DFT-D3, which are capable of describing the influence of the van der Waals (vdW) interactions. In particular, we discuss the effect of dispersion on the work of adhesion (W-ad), equilibrium geometries, and charge density difference (CDD). We find that replacement of P by Si in a-HAP (a-Si-HAP) with the creation of OH vacancies as charge compensation results in a significant increase in the bond strength between the coating and substrate in the case of using the PBE-GGA functional. However, including the vdW interactions shows that these forces considerably contribute to the W-ad. We show that the difference (W-ad – W-ad(vdW)) is on average more than 1.1 J m(-2) and 0.5 J m(-2) for a-HAP/a-TiO2 and a-Si-HAP/a-TiO2, respectively. These results reveal that including vdW interactions is essential for accurately describing the chemical bonding at the a-HAP/a-TiO2 interface.

Abstract: While spectroscopic data on small hydrocarbons in interstellar media in combination with crossed molecular beam (CMB) experiments have provided a wealth of information on astrochemically relevant species, much of the underlying mechanistic pathways of their formation remain elusive. Therefore, in this work, the chemical reaction mechanisms of C(³P_J) + C₆H₆and C⁺(²P) + C₆H₆systems using the quantum mechanical molecular dynamics (QMMD) technique at the PBE0-D3(BJ) level of theory is investigated, mimicking a CMB experiment. Both the dynamics of the reactions as well as the electronic structure for the purpose of the reaction network are evaluated. The method is validated for the first reaction by comparison to the available experimental data. The reaction scheme for the C(³P_J) + C₆H₆system covers the literature data,<italic>e.g.</italic>the major products are the 1,2-didehydrocycloheptatrienyl radical (C₇H₅) and benzocyclopropenyl radical (C₆H₅–CH), and it reveals the existence of less common pathways for the first time. The chemistry of the C⁺(²P_J) + C₆H₆system is found to be much richer, and we have found that this is because of more exothermic reactions in this system in comparison to those in the C(³P_J) + C₆H₆system. Moreover, using the QMMD simulation, a number of reaction paths have been revealed that produce three distinct classes of reaction products with different ring sizes. All in all, at all the collision energies and orientations, the major product is the heptagon molecular ion for the ionic system. It is also revealed that the collision orientation has a dominant effect on the reaction products in both systems, while the collision energy mostly affects the charged system. These simulations both prove the applicability of this approach to simulate crossed molecular beams, and provide fundamental information on reactions relevant for the interstellar medium.

Abstract: It is well established that lipid aldehydes (LAs) are able to increase the permeability of cell membranes and induce their rupture. However, it is not yet clear how LAs are distributed in phase-separated membranes (PSMs), which are responsible for the transport of selected molecules and intracellular signaling. Thus, we investigate here the distribution of LAs in a PSM by coarse-grained molecular dynamics simulations. Our results reveal that LAs derived from mono-unsaturated lipids tend to accumulate at the interface between the liquid-ordered/liquiddisordered domains, whereas those derived from poly-unsaturated lipids remain in the liquid-disordered domain. These results are important for understanding the effects caused by oxidized lipids in membrane structure, properties and organization.

Abstract: The 2022 Roadmap is the next update in the series of Plasma Roadmaps published by<italic>Journal of Physics</italic>D with the intent to identify important outstanding challenges in the field of low-temperature plasma (LTP) physics and technology. The format of the Roadmap is the same as the previous Roadmaps representing the visions of 41 leading experts representing 21 countries and five continents in the various sub-fields of LTP science and technology. In recognition of the evolution in the field, several new topics have been introduced or given more prominence. These new topics and emphasis highlight increased interests in plasma-enabled additive manufacturing, soft materials, electrification of chemical conversions, plasma propulsion, extreme plasma regimes, plasmas in hypersonics, data-driven plasma science and technology and the contribution of LTP to combat COVID-19. In the last few decades, LTP science and technology has made a tremendously positive impact on our society. It is our hope that this roadmap will help continue this excellent track record over the next 5–10 years.

Abstract: Utilizing ammonia as a storage medium for hydrogen is currently receiving increased attention. A possible method to retrieve the hydrogen is by plasma-catalytic decomposition. In this work, we combined an experimental study, using a dielectric barrier discharge plasma reactor, with a plasma kinetic model, to get insights into the decomposition mechanism. The experimental results revealed a similar effect on the ammonia conversion when changing the flow rate and power, where increasing the specific energy input (higher power or lower flow rate) gave an increased conversion. A conversion as high as 82 % was achieved at a specific energy input of 18 kJ/Nl. Furthermore, when changing the discharge volume from 31 to 10 cm3, a change in the plasma distribution factor from 0.2 to 0.1 was needed in the model to best describe the conversions of the experiments. This means that a smaller plasma volume caused a higher transfer of energy through micro-discharges (non-uniform plasma), which was found to promote the decomposition of ammonia. These results indicate that it is the collisions between NH3 and the high-energy electrons that initiate the decomposition. Moreover, the rate of ammonia destruction was found by the model to be in the order of 1022 molecules/(cm3 s) during the micro-discharges, which is 5 to 6 orders of magnitude higher than in the afterglows. A considerable re-formation of ammonia was found to take place in the afterglows, limiting the overall conversion. In addition, the model revealed that implementation of packing material in the plasma introduced high concentrations of surface-bound hydrogen atoms, which introduced an additional ammonia re-formation pathway through an Eley-Rideal reaction with gas phase NH2. Furthermore, a more uniform plasma is predicted in the presence of MgAl2O4, which leads to a lower average electron energy during micro-discharges and a lower conversion (37 %) at a comparable residence time for the plasma alone (51 %).

Abstract: We present a zero-dimensional kinetic model to characterise specifically the gas-phase dynamics of methane
conversion in a nanosecond pulsed discharge (NPD) plasma reactor. The model includes a systematic approach to
capture the nanoscale power discharges and the rapid ensuing changes in electric field, gas and electron temperature,
as well as species densities. The effects of gas temperature and reactor pressure on gas conversion and
product selectivity are extensively investigated and validated against experimental work. We discuss the
important reaction pathways and provide an analysis of the dynamics of the heating and cooling mechanisms. H
radicals are found to be the most populous plasma species and they participate in hydrogenation and dehydrogenation
reactions, which are the dominant recombination reactions leading to C2H4 and C2H2 as main
products (depending on the pressure).