Acquired non-thermal plasma resistance mediates a shift towards aerobic glycolysis and ferroptotic cell death in melanoma.

AIMS: To gain insights into the underlying mechanisms of NTP therapy sensitivity and resistance, using the first-ever NTP-resistant cell line derived from sensitive melanoma cells (A375). METHODS: Melanoma cells were exposed to NTP and re-cultured for 12 consecutive weeks before evaluation against the parental control cells. Whole transcriptome sequencing analysis was performed to identify differentially expressed genes and enriched molecular pathways. Glucose uptake, extracellular lactate, media acidification, and mitochondrial respiration was analyzed to determine metabolic changes. Cell death inhibitors were used to assess the NTP-induced cell death mechanisms, and apoptosis and ferroptosis was further validated via Annexin V, Caspase 3/7, and lipid peroxidation analysis. RESULTS: Cells continuously exposed to NTP became 10 times more resistant to NTP compared to the parental cell line of the same passage, based on their half-maximal inhibitory concentration (IC50). Sequencing and metabolic analysis indicated that NTP-resistant cells had a preference towards aerobic glycolysis, while cell death analysis revealed that NTP-resistant cells exhibited less apoptosis but were more vulnerable to lipid peroxidation and ferroptosis. CONCLUSIONS: A preference towards aerobic glycolysis and ferroptotic cell death are key physiological changes in NTP-resistance cells, which opens new avenues for further, in-depth research into other cancer types.

Effect of lipid oxidation on the channel properties of Cx26 hemichannels: A molecular dynamics study.

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.

Distribution of lipid aldehydes in phase-separated membranes: A molecular dynamics study.

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/liquid-disordered 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.

Effect of Cysteine Oxidation in SARS-CoV-2 Receptor-Binding Domain on Its Interaction with Two Cell Receptors: Insights from Atomistic Simulations.

Binding of the SARS-CoV-2 S-glycoprotein to cell receptors is vital for the entry of the virus into cells and subsequent infection. ACE2 is the main cell receptor for SARS-CoV-2, which can attach to the C-terminal receptor-binding domain (RBD) of the SARS-CoV-2 S-glycoprotein. The GRP78 receptor plays an anchoring role, which attaches to the RBD and increases the chance of other RBDs binding to ACE2. Although high levels of reactive oxygen and nitrogen species (RONS) are produced during viral infections, it is not clear how they affect the RBD structure and its binding to ACE2 and GRP78. In this research, we apply molecular dynamics simulations to study the effect of oxidation of the highly reactive cysteine (Cys) amino acids of the RBD on its binding to ACE2 and GRP78. The interaction energy of both ACE2 and GRP78 with the whole RBD, as well as with the RBD main regions, is compared in both the native and oxidized RBDs. Our results show that the interaction energy between the oxidized RBD and ACE2 is strengthened by 155 kJ/mol, increasing the binding of the RBD to ACE2 after oxidation. In addition, the interaction energy between the RBD and GRP78 is slightly increased by 8 kJ/mol after oxidation, but this difference is not significant. Overall, these findings highlight the role of RONS in the binding of the SARS-CoV-2 S-glycoprotein to host cell receptors and suggest an alternative mechanism by which RONS could modulate the entrance of viral particles into the cells.

Cold Atmospheric Plasma Does Not Affect Stellate Cells Phenotype in Pancreatic Cancer Tissue in Ovo.

Pancreatic ductal adenocarcinoma (PDAC) is a challenging neoplastic disease, mainly due to the development of resistance to radio- and chemotherapy. Cold atmospheric plasma (CAP) is an alternative technology that can eliminate cancer cells through oxidative damage, as shown in vitro, in ovo, and in vivo. However, how CAP affects the pancreatic stellate cells (PSCs), key players in the invasion and metastasis of PDAC, is poorly understood. This study aims to determine the effect of an anti-PDAC CAP treatment on PSCs tissue developed in ovo using mono- and co-cultures of RLT-PSC (PSCs) and Mia PaCa-2 cells (PDAC). We measured tissue reduction upon CAP treatment and mRNA expression of PSC activation markers and extracellular matrix (ECM) remodelling factors via qRT-PCR. Protein expression of selected markers was confirmed via immunohistochemistry. CAP inhibited growth in Mia PaCa-2 and co-cultured tissue, but its effectiveness was reduced in the latter, which correlates with reduced ki67 levels. CAP did not alter the mRNA expression of PSC activation and ECM remodelling markers. No changes in MMP2 and MMP9 expression were observed in RLT-PSCs, but small changes were observed in Mia PaCa-2 cells. Our findings support the ability of CAP to eliminate PDAC cells, without altering the PSCs.

Flowing Atmospheric Pressure Afterglow for Ambient Ionization: Reaction Pathways Revealed by Modeling.

We describe the plasma chemistry in a helium flowing atmospheric pressure afterglow (FAPA) used for analytical spectrometry, by means of a quasi-one-dimensional (1D) plasma chemical kinetics model. We study the effect of typical impurities present in the feed gas, as well as the afterglow in ambient humid air. The model provides the species density profiles in the discharge and afterglow regions and the chemical pathways. We demonstrate that H, N, and O atoms are formed in the discharge region, while the dominant reactive neutral species in the afterglow are O3 and NO. He* and He2* are responsible for Penning ionization of O2, N2, H2O, H2, and N, and especially O and H atoms. Besides, He2+ also contributes to ionization of N2, O2, H2O, and O through charge transfer reactions. From the pool of ions created in the discharge, NO+ and (H2O)3H+ are the dominant ions in the afterglow. Moreover, negatively charged clusters, such as NO3H2O- and NO2H2O-, are formed and their pathway is discussed as well. Our model predictions are in line with earlier observations in the literature about the important reagent ions and provide a comprehensive overview of the underlying pathways. The model explains in detail why helium provides a high analytical sensitivity because of high reagent ion formation by both Penning ionization and charge transfer. Such insights are very valuable for improving the analytical performance of this (and other) ambient desorption/ionization source(s).

Lipid Oxidation: Role of Membrane Phase-Separated Domains.

Lipid oxidation is associated with several inflammatory and neurodegenerative diseases, but many questions to unravel its effects on biomembranes are still open due to the complexity of the topic. For instance, recent studies indicated that phase-separated domains can have a significant effect on membrane function. It is reported that domain interfaces are "hot spots" for pore formation, but the underlying mechanisms and the effect of oxidation-induced phase separation on membranes remain elusive. Thus, to evaluate the permeability of the membrane coexisting of liquid-ordered (Lo) and liquid-disordered (Ld) domains, we performed atomistic molecular dynamics simulations. Specifically, we studied the membrane permeability of nonoxidized or oxidized homogeneous membranes (single-phase) and at the Lo/Ld domain interfaces of heterogeneous membranes, where the Ld domain is composed of either oxidized or nonoxidized lipids. Our simulation results reveal that the addition of only 1.5% of lipid aldehyde molecules at the Lo/Ld domain interfaces of heterogeneous membranes increases the membrane permeability, whereas their addition at homogeneous membranes does not have any effect. This study is of interest for a better understanding of cancer treatment methods based on oxidative stress (causing among others lipid oxidation), such as plasma medicine and photodynamic therapy.

Rational design of an XNA ligase through docking of unbound nucleic acids to toroidal proteins.

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.

The Quest to Quantify Selective and Synergistic Effects of Plasma for Cancer Treatment: Insights from Mathematical Modeling.

Cold atmospheric plasma (CAP) and plasma-treated liquids (PTLs) have recently become a promising option for cancer treatment, but the underlying mechanisms of the anti-cancer effect are still to a large extent unknown. Although hydrogen peroxide (H2O2) has been recognized as the major anti-cancer agent of PTL and may enable selectivity in a certain concentration regime, the co-existence of nitrite can create a synergistic effect. We develop a mathematical model to describe the key species and features of the cellular response toward PTL. From the numerical solutions, we define a number of dependent variables, which represent feasible measures to quantify cell susceptibility in terms of the H2O2 membrane diffusion rate constant and the intracellular catalase concentration. For each of these dependent variables, we investigate the regimes of selective versus non-selective, and of synergistic versus non-synergistic effect to evaluate their potential role as a measure of cell susceptibility. Our results suggest that the maximal intracellular H2O2 concentration, which in the selective regime is almost four times greater for the most susceptible cells compared to the most resistant cells, could be used to quantify the cell susceptibility toward exogenous H2O2. We believe our theoretical approach brings novelty to the field of plasma oncology, and more broadly, to the field of redox biology, by proposing new ways to quantify the selective and synergistic anti-cancer effect of PTL in terms of inherent cell features.

How do nitrated lipids affect the properties of phospholipid membranes?

Biological membranes are under constant attack of free radicals, which may lead to lipid nitro-oxidation, producing 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 neurodegenerative 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 phospholipid 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.

Modeling plasmas in analytical chemistry-an example of cross-fertilization.

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. Graphical abstract.

Towards Green Ammonia Synthesis through Plasma-Driven Nitrogen Oxidation and Catalytic Reduction.

Ammonia is an industrial large-volume chemical, with its main application in fertilizer production. It also attracts increasing attention as a green-energy vector. Over the past century, ammonia production has been dominated by the Haber-Bosch process, in which a mixture of nitrogen and hydrogen gas is converted to ammonia at high temperatures and pressures. Haber-Bosch processes with natural gas as the source of hydrogen are responsible for a significant share of the global CO2 emissions. Processes involving plasma are currently being investigated as an alternative for decentralized ammonia production powered by renewable energy sources. In this work, we present the PNOCRA process (plasma nitrogen oxidation and catalytic reduction to ammonia), combining plasma-assisted nitrogen oxidation and lean NOx trap technology, adopted from diesel-engine exhaust gas aftertreatment technology. PNOCRA achieves an energy requirement of 4.6 MJ mol-1 NH3 , which is more than four times less than the state-of-the-art plasma-enabled ammonia synthesis from N2 and H2 with reasonable yield (>1 %).

Effect of oxidative stress on cystine transportation by xC‾ antiporter.

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.

Transport of cystine across xC- antiporter.

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.

Reactivity and stability of plasma-generated oxygen and nitrogen species in buffered water solution: a computational study.

Plasma-treated liquids have great potential for biomedical applications. However, insight into the underlying mechanisms and the exact chemistry is still scarce. In this study, we present the combination of a 0D chemical kinetics and a 2D fluid dynamics model to investigate the plasma treatment of a buffered water solution with the kINPen® plasma jet. Using this model, we calculated the gas and liquid flow profiles and the transport and chemistry of all species in the gas and the liquid phase. Moreover, we evaluated the stability of the reactive oxygen and nitrogen species after plasma treatment. We found that of all species, only H2O2, HNO2/NO2-, and HNO3/NO3- are stable in the buffered solution after plasma treatment. This is because both their production and loss processes in the liquid phase are dependent on short-lived radicals (e.g. OH, NO, and NO2). Apart from some discrepancy in the absolute values of the concentrations, which can be explained by the model, all general trends and observations in our model are in qualitative agreement with experimental data and literature.

Reaction of chloride anion with atomic oxygen in aqueous solutions: can cold plasma help in chemistry research?

Cold atmospheric plasma in contact with solutions has many applications, but its chemistry contains many unknowns such as the undescribed reactions with solutes. By combining experiments and modelling, we report the first direct demonstration of the reaction of chloride with oxygen atoms in aqueous solutions exposed to cold plasma.

How membrane lipids influence plasma delivery of reactive oxygen species into cells and subsequent DNA damage: an experimental and computational study.

The mechanisms of plasma in medicine are broadly attributed to plasma-derived reactive oxygen and nitrogen species (RONS). In order to exert any intracellular effects, these plasma-derived RONS must first traverse a major barrier in the cell membrane. The cell membrane lipid composition, and thereby the magnitude of this barrier, is highly variable between cells depending on type and state (e.g. it is widely accepted that healthy and cancerous cells have different membrane lipid compositions). In this study, we investigate how plasma-derived RONS interactions with lipid membrane components can potentially be exploited in the future for treatment of diseases. We couple phospholipid vesicle experiments, used as simple cell models, with molecular dynamics (MD) simulations of the lipid membrane to provide new insights into how the interplay between phospholipids and cholesterol may influence the response of healthy and diseased cell membranes to plasma-derived RONS. We focus on the (i) lipid tail saturation degree, (ii) lipid head group type, and (iii) membrane cholesterol fraction. Using encapsulated molecular probes, we study the influence of the above membrane components on the ingress of RONS into the vesicles, and subsequent DNA damage. Our results indicate that all of the above membrane components can enhance or suppress RONS uptake, depending on their relative concentration within the membrane. Further, we show that higher RONS uptake into the vesicles does not always correlate with increased DNA damage, which is attributed to ROS reactivity and lifetime. The MD simulations indicate the multifactorial chemical and physical processes at play, including (i) lipid oxidation, (ii) lipid packing, and (iii) lipid rafts formation. The methods and findings presented here provide a platform of knowledge that could be leveraged in the development of therapies relying on the action of plasma, in which the cell membrane and oxidative stress response in cells is targeted.

Analysis of Short-Lived Reactive Species in Plasma-Air-Water Systems: The Dos and the Do Nots.

This Feature addresses the analysis of the reactive species generated by nonthermal atmospheric pressure plasmas, which are widely employed in industrial and biomedical research, as well as first clinical applications. We summarize the progress in detection of plasma-generated short-lived reactive oxygen and nitrogen species in aqueous solutions, discuss the potential and limitations of various analytical methods in plasma-liquid systems, and provide an outlook on the possible future research goals in development of short-lived reactive species analysis methods for a general nonspecialist audience.

Impact of the Particle Diameter on Ion Cloud Formation from Gold Nanoparticles in ICPMS.

The unique capabilities of microsecond dwell time (DT) single-particle inductively coupled plasma mass spectrometry (spICPMS) were utilized to characterize the cloud of ions generated from the introduction of suspensions of gold nanoparticles (AuNPs) into the plasma. A set of narrowly distributed particles with diameters ranging from 15.4 to 100.1 nm was synthesized and characterized according to established protocols. Statistically significant numbers of the short transient spICPMS events were evaluated by using 50 µs DT for their summed intensity, maximum intensity, and duration, of which all three were found to depend on the particle diameter. The summed intensity increases from 10 to 1661 counts and the maximum intensity from 6 to 309 counts for AuNPs with diameters from 15.4 to 83.2 nm. The event duration rises from 322 to 1007 µs upon increasing AuNP diameter. These numbers represent a comprehensive set of key data points of the ion clouds generated in ICPMS from AuNPs. The extension of event duration is of high interest to appoint the maximum possible particle number concentration at which separation of consecutive events in spICPMS can still be achieved. Moreover, the combined evaluation of all above-mentioned ion cloud characteristics can explain the regularly observed prolonged single-particle events. The transport and ionization behavior of AuNPs in the ICP was also computationally modeled to gain insight into the size-dependent signal generation. The simulated data reveals that the plasma temperature, and therefore the point of ionization of the particles, is the same for all diameters. However, the maximum number density of Au+, as well as the extent of the ion cloud, depends on the particle diameter, in agreement with the experimental data, and it provides an adequate explanation for the observed ion cloud characteristics.

Toward the Understanding of Selective Si Nano-Oxidation by Atomic Scale Simulations.

The continuous miniaturization of nanodevices, such as transistors, solar cells, and optical fibers, requires the controlled synthesis of (ultra)thin gate oxides (<10 nm), including Si gate-oxide (SiO2) with high quality at the atomic scale. Traditional thermal growth of SiO2 on planar Si surfaces, however, does not allow one to obtain such ultrathin oxide due to either the high oxygen diffusivity at high temperature or the very low sticking ability of incident oxygen at low temperature. Two recent techniques, both operative at low (room) temperature, have been put forward to overcome these obstacles: (i) hyperthermal oxidation of planar Si surfaces and (ii) thermal or plasma-assisted oxidation of nonplanar Si surfaces, including Si nanowires (SiNWs). These nano-oxidation processes are, however, often difficult to study experimentally, due to the key intermediate processes taking place on the nanosecond time scale. In this Account, these Si nano-oxidation techniques are discussed from a computational point of view and compared to both hyperthermal and thermal oxidation experiments, as well as to well-known models of thermal oxidation, including the Deal-Grove, Cabrera-Mott, and Kao models and several alternative mechanisms. In our studies, we use reactive molecular dynamics (MD) and hybrid MD/Monte Carlo simulation techniques, applying the Reax force field. The incident energy of oxygen species is chosen in the range of 1-5 eV in hyperthermal oxidation of planar Si surfaces in order to prevent energy-induced damage. It turns out that hyperthermal growth allows for two growth modes, where the ultrathin oxide thickness depends on either (1) only the kinetic energy of the incident oxygen species at a growth temperature below Ttrans = 600 K, or (2) both the incident energy and the growth temperature at a growth temperature above Ttrans. These modes are specific to such ultrathin oxides, and are not observed in traditional thermal oxidation, nor theoretically considered by already existing models. In the case of thermal or plasma-assisted oxidation of small Si nanowires, on the other hand, the thickness of the ultrathin oxide is a function of the growth temperature and the nanowire diameter. Below Ttrans, which varies with the nanowire diameter, partially oxidized SiNW are formed, whereas complete oxidation to a SiO2 nanowire occurs only above Ttrans. In both nano-oxidation processes at lower temperature (T < Ttrans), final sandwich c-Si|SiOx|a-SiO2 structures are obtained due to a competition between overcoming the energy barrier to penetrate into Si subsurface layers and the compressive stress (∼2-3 GPa) at the Si crystal/oxide interface. The overall atomic-simulation results strongly indicate that the thickness of the intermediate SiOx (x < 2) region is very limited (∼0.5 nm) and constant irrespective of oxidation parameters. Thus, control over the ultrathin SiO2 thickness with good quality is indeed possible by accurately tuning the oxidant energy, oxidation temperature and surface curvature. In general, we discuss and put in perspective these two oxidation mechanisms for obtaining controllable ultrathin gate-oxide films, offering a new route toward the fabrication of nanodevices via selective nano-oxidation.