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: Non-thermal plasma (NTP), an ionized gas generated at ambient pressure and temperature, has been an emerging technology for medical applications. Through controlled delivery of reactive oxygen and nitrogen species (ROS/RNS), NTP can elicit hormetic cellular responses, thus stimulating broad therapeutic effects. To enable clinical translation of the promising preclinical research into NTP therapy, a deeper understanding of NTP interactions with clinical substrates is profoundly needed. Since NTP-generated ROS/RNS will inevitably interact with blood in several clinical contexts, understanding their stability in this system is crucial. In this study, two medically relevant NTP delivery modalities were used to assess the stability of NTP-generated ROS/RNS in three aqueous solutions with increasing organic complexities: phosphate-buffered saline (PBS), blood plasma (BP), and processed whole blood. NTP-generated RNS collectively (NO2−, ONOO−), H2O2, and ONOO− exclusively were analyzed over time. We demonstrated that NTP-generated RNS and H2O2 were stable in PBS but scavenged by different components of the blood. While RNS remained stable in BP after initial scavenging effects, it was completely reduced in processed whole blood. On the other hand, H2O2 was completely scavenged in both liquids over time. Our previously developed luminescent probe europium(III) was used for precision measurement of ONOO− concentration. NTP-generated ONOO− was detected in all three liquids for up to at least 30 seconds, thus highlighting its therapeutic potential. Based on our results, we discussed the necessary considerations to choose the most optimal NTP modality for delivery of ROS/RNS to and via blood in the clinical context.

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: Hoofd-halskanker vormt een aanzienlijke uitdaging met bijna 900.000 nieuwe diagnoses per jaar, waarbij de jaarlijkse incidentie blijft stijgen. Vaak wordt de diagnose pas in een laat stadium gesteld, wat complexe behandelingen noodzakelijk maakt. Terugval van patiënten is helaas een veelvoorkomend probleem. De gemiddelde overlevingsduur is beperkt tot enkele maanden. Daarom is er een dringende behoefte om nieuwe, veelbelovende behandelingen te ontwikkelen voor patiënten met hoofd-halskanker. Voor het bereiken van deze vooruitgang spelen innovatieve studiemodellen een cruciale rol. Het ontwikkelen van deze nieuwe behandelingen start met laboratoriumonderzoek, waarbij traditionele tweedimensionale celculturen hun beperkingen hebben. Daarom verschuiven onderzoekers hun aandacht meer en meer naar geavanceerdere driedimensionale modellen, met hoofd-halskankerorganoïden als beloftevol nieuw model. Dit model behoudt immers zowel het genetische profiel als de morfologische kenmerken van de originele tumor van de hoofd-halskankerpatiënt. Hoofdhalskankerorganoïden bieden daarom de mogelijkheid om innovatieve behandelingen te testen en kunnen mogelijk zelfs de respons van een patiënt op bepaalde therapieën voorspellen. Hoewel tumororganoïden als ‘patiënt-in-het-lab’ veelbelovend zijn, zijn er uitdagingen te overwinnen, zoals de ontwikkelingstijd en de toepasbaarheid bij alle tumortypes, evenals het ontbreken van immuuncellen en andere micro-omgevingscomponenten. Er is daarom een grote behoefte aan gestandaardiseerde protocollen voor de ontwikkeling van organoïden en verkorting van de ontwikkelingstijd. Concluderend bieden driedimensionale hoofd-halskankerorganoïden een veelbelovend perspectief voor de toekomst van kankerbehandelingen. Ze hebben het potentieel om bij te dragen aan de ontwikkeling van gepersonaliseerde behandelingen en zo de overlevingskansen van kankerpatiënten te verbeteren. Het is echter belangrijk om hun voorspellend vermogen en toepassingsmogelijkheden verder te onderzoeken, voordat ze op grote schaal worden geïmplementeerd.

Abstract: We investigate the theoretical possibility of single layer graphene formation on a nickel surface at different substrate temperatures under far-from-equilibrium high precursor flux conditions, employing state-of-the-art hybrid reactive molecular dynamics/uniform acceptance force bias Monte Carlo simulations. It is predicted that under these conditions, the formation of a single layer graphene-like film may proceed through a combined depositionsegregation mechanism on a nickel substrate, rather than by pure surface segregation as is typically observed for metals with high carbon solubility. At 900 K and above, nearly continuous graphene layers are obtained. These simulations suggest that single layer graphene deposition is theoretically possible on Ni under high flux conditions.

Abstract: The application of coreshell SiSiO2 nanowires as nanoelectronic devices strongly depends on their structure, which is difficult to tune precisely. In this work, we investigate the formation of the coreshell nanowires at the atomic scale, by reactive molecular dynamics simulations. The occurrence of two temperature-dependent oxidation mechanisms of ultra-small diameter Si-NWs is demonstrated. We found that control over the Si-core radius and the SiOx (x ≤ 2) oxide shell is possible by tuning the growth temperature and the initial Si-NW diameter. Two different structures were obtained, i.e., ultrathin SiO2 silica nanowires at high temperature and Si core|ultrathin SiO2 silica nanowires at low temperature. The transition temperature is found to linearly decrease with the nanowire curvature. Finally, the interfacial stress is found to be responsible for self-limiting oxidation, depending on both the initial Si-NW radius and the oxide growth temperature. These novel insights allow us to gain control over the exact morphology and structure of the wires, as is needed for their application in nanoelectronics.

Abstract: Using reactive molecular dynamics simulations, the melting behavior of nickel-carbon nanoclusters is examined. The phase diagrams of icosahedral and Wulff polyhedron clusters are determined using both the Lindemann index and the potential energy. Formulae are derived for calculating the equilibrium constants and the solid and liquid fractions during a phase transition, allowing more rational determination of the melting temperature with respect to the arbitrary Lindemann value. These results give more insight into the properties of nickel-carbon nanoclusters in general and can specifically be very useful for a better understanding of the synthesis of carbon nanotubes using the catalytic chemical vapor deposition method.

Abstract: We performed molecular dynamics simulations to investigate the effect of lipid peroxidation products on the structural and dynamic properties of the cell membrane. Our simulations predict that the lipid order in a phospholipid bilayer, as a model system for the cell membrane, decreases upon addition of lipid peroxidation products. Eventually, when all phospholipids are oxidized, pore formation can occur. This will allow reactive species, such as reactive oxygen and nitrogen species (RONS), to enter the cell and cause oxidative damage to intracellular macromolecules, such as DNA or proteins. On the other hand, upon increasing the cholesterol fraction of lipid bilayers, the cell membrane order increases, eventually reaching a certain threshold, from which cholesterol is able to protect the membrane against pore formation. This finding is crucial for cancer treatment by plasma technology, producing a large number of RONS, as well as for other cancer treatment methods that cause an increase in the concentration of extracellular RONS. Indeed, cancer cells contain less cholesterol than their healthy counterparts. Thus, they will be more vulnerable to the consequences of lipid peroxidation, eventually enabling the penetration of RONS into the interior of the cell, giving rise to oxidative stress, inducing pro-apoptotic factors. This provides, for the first time, molecular level insight why plasma can selectively treat cancer cells, while leaving their healthy counterparts undamaged, as is indeed experimentally demonstrated.

Abstract: This study reports on the possible effects of OH radical impact on the transmembrane domain 6 of P-glycoprotein, TM6, which plays a crucial role in drug binding in human cells. For the first time, we employ molecular dynamics (MD) simulations based on the self-consistent charge density functional tight binding (SCC-DFTB) method to elucidate the potential sites of fragmentation and mutation in this domain upon impact of OH radicals, and to obtain fundamental information about the underlying reaction mechanisms. Furthermore, we apply non-reactive MD simulations to investigate the long-term effect of this mutation, with possible implications for drug binding. Our simulations indicate that the interaction of OH radicals with TM6 might lead to the breaking of C-C and C-N peptide bonds, which eventually cause fragmentation of TM6. Moreover, according to our simulations, the OH radicals can yield mutation in the aromatic ring of phenylalanine in TM6, which in turn affects its structure. As TM6 plays an important role in the binding of a range of cytotoxic drugs with P-glycoprotein, any changes in its structure are likely to affect the response of the tumor cell in chemotherapy. This is crucial for cancer therapies based on reactive oxygen species, such as plasma treatment.

Abstract: In recent years, the ability of cold atmospheric pressure plasmas (CAPS) to selectively induce cell death in cancer cells has been widely established. This selectivity has been assigned to the reactive oxygen and nitrogen species (RONS) created in CAPs. To provide new insights in the search for an explanation
for the observed selectivity, we calculate the transfer free energy of multiple ROS across membranes containing a varying amount of cholesterol. The cholesterol fraction is investigated as a selectivity parameter because membranes of cancer cells are known to contain lower fractions of cholesterol compared to healthy cells. We find that cholesterol has a significant effect on the permeation of
reactive species across a membrane. Indeed, depending on the specific reactive species, an increasing cholesterol fraction can lead to (i) an increase of the transfer free energy barrier height and width, (ii) the formation of a local free energy minimum in the center of the membrane and (iii) the creation of extra free energy barriers due to the bulky sterol rings. In the context of plasma oncology, these observations suggest that the increased ingress of RONS in cancer cells can be explained by the decreased cholesterol fraction of their cell membrane.

Abstract: We evaluate the anti-cancer capacity of plasma-treated PBS (pPBS), by measuring the concentrations of NO2 − and H2O2 in pPBS, treated with a plasma jet, for different values of gas flow rate, gap and plasma treatment time, as well as the effect of pPBS on cancer cell cytotoxicity, for three different glioblastoma cancer cell lines, at exactly the same plasma treatment conditions. Our experiments reveal that pPBS is cytotoxic for all conditions investigated. A small variation in gap between plasma jet and liquid surface (10 mm vs 15 mm) significantly affects the chemical composition of pPBS and its anti-cancer capacity, attributed to the occurrence of discharges onto the liquid. By correlating the effect of gap, gas flow rate and plasma treatment time on the chemical composition and anti-cancer capacity of pPBS, we may conclude that H2O2 is a more important species for the anti-cancer capacity of pPBS than NO2 −. We also used a 0D model, developed for plasma-liquid interactions, to elucidate the most important mechanisms for the generation of H2O2 and NO2 −. Finally, we found that pPBS might be more suitable for practical applications in a clinical setting than (commonly used) plasma-activated media (PAM), because of its higher stability.

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.