Exploring α-synuclein stability under the external electrostatic fields: Effect of repeat unit.

Parkinson's disease (PD) is a category of neurodegenerative disorders (ND) that currently lack comprehensive and definitive treatment strategies. The etiology of PD can be attributed to the presence and aggregation of a protein known as α-synuclein. Researchers have observed that the application of an external electrostatic field holds the potential to induce the separation of the fibrous structures into peptides. To comprehend this phenomenon, our investigation involved simulations conducted on the α-synuclein peptides through the application of Molecular Dynamics (MD) simulation techniques under the influence of a 0.1 V/nm electric field. The results obtained from the MD simulations revealed that in the presence of external electric field, the monomer and oligomeric forms of α-synuclein are experienced significant conformational changes which could prevent them from further aggregation. However, as the number of peptide units in the model system increases, forming trimers and tetramers, the stability against the electric field also increases. This enhanced stability in larger aggregates indicates a critical threshold in α-synuclein assembly where the electric field's effectiveness in disrupting the aggregation diminishes. Therefore, our findings suggest that early diagnosis and intervention could be crucial in preventing PD progression. When α-synuclein predominantly exists in its monomeric or dimeric form, applying even a lower electric field could effectively disrupt the initial aggregation process. Inhibition of α-synuclein fibril formation at early stages might serve as a viable solution to combat PD by halting the formation of more stable and pathogenic α-synuclein fibrils.

In silico study of the impact of oxidation on pyruvate transmission across the hVDAC1 protein channel.

The overexpression of voltage dependent anion channels (VDACs), particularly VDAC1, in cancer cells compared to normal cells, plays a crucial role in cancer cell metabolism, apoptosis regulation, and energy homeostasis. In this study, we used molecular dynamics (MD) simulations to investigate the effect of a low level of VDAC1 oxidation (induced e.g., by cold atmospheric plasma (CAP)) on the pyruvate (Pyr) uptake by VDAC1. Inhibiting Pyr uptake through VDAC1 can suppress cancer cell proliferation. Our primary target was to study the translocation of Pyr across the native and oxidized forms of hVDAC1, the human VDAC1. Specifically, we employed MD simulations to analyze the hVDAC1 structure by modifying certain cysteine residues to cysteic acids and methionine residues to methionine sulfoxides, which allowed us to investigate the effect of oxidation. Our results showed that the free energy barrier for Pyr translocation through the native and oxidized channel was approximately 4.3 ± 0.7 kJ mol-1 and 10.8 ± 1.8 kJ mol-1, respectively. An increase in barrier results in a decrease in rate of Pyr permeation through the oxidized channel. Thus, our results indicate that low levels of CAP oxidation reduce Pyr translocation, resulting in decreased cancer cell proliferation. Therefore, low levels of oxidation are likely sufficient to treat cancer cells given the inhibition of Pyr uptake.

An exploration of the binding prediction of anatoxin-a and atropine to acetylcholinesterase enzyme using multi-level computer simulations.

Acetylcholinesterase (AChE) is crucial for the breakdown of acetylcholine to acetate and choline, while the inhibition of AChE by anatoxin-a (ATX-a) results in severe health complications. This study explores the structural characteristics of ATX-a and its interactions with AChE, comparing to the reference molecule atropine for binding mechanisms. Molecular docking simulations reveal strong binding affinity of both ATX-a and atropine to AChE, interacting effectively with specific amino acids in the binding site as potential inhibitors. Quantitative assessment using the MM-PBSA method demonstrates a significantly negative binding free energy of -81.659 kJ mol-1for ATX-a, indicating robust binding, while atropine exhibits a stronger binding affinity with a free energy of -127.565 kJ mol-1. Umbrella sampling calculates the ΔGbindvalues to evaluate binding free energies, showing a favorable ΔGbindof -36.432 kJ mol-1for ATX-a and a slightly lower value of -30.12 kJ mol-1for atropine. This study reveals the dual functionality of ATX-a, acting as both a nicotinic acetylcholine receptor agonist and an AChE inhibitor. Remarkably, stable complexes form between ATX-a and atropine with AChE at its active site, exhibiting remarkable binding free energies. These findings provide valuable insights into the potential use of ATX-a and atropine as promising candidates for modulating AChE activity.

Effect of Electric Field on α-Synuclein Fibrils: Revealed by Molecular Dynamics Simulations.

The self-association of amylogenic proteins to the fibril form is considered a pivotal factor in the pathogenesis of neurodegenerative diseases, including Parkinson's disease (PD). PD causes unintended or uncontrollable movements in its common symptoms. α-synuclein is the major cause of PD development and thus has been the main target of numerous studies to suppress and sequester its expression or effectively degrade it. Nonetheless, to date, there are no efficient and proven ways to prevent pathological protein aggregation. Recent investigations proposed applying an external electric field to interrupt the fibrils. This method is a non-invasive approach that has a certain benefit over others. We performed molecular dynamics (MD) simulations by applying an electric field on highly toxic fibrils of α-synuclein to gain a molecular-level insight into fibril disruption mechanisms. The results revealed that the applied external electric field induces substantial changes in the conformation of the α-synuclein fibrils. Furthermore, we show the threshold value for electric field strength required to completely disrupt the α-synuclein fibrils by opening the hydrophobic core of the fibril. Thus, our findings might serve as a valuable foundation to better understand molecular-level mechanisms of the α-synuclein fibrils disaggregation process under an applied external electric field.

Mechanistic Insight into Permeation of Plasma-Generated Species from Vacuum into Water Bulk.

Due to their potential benefits, cold atmospheric plasmas (CAPs), as biotechnological tools, have been used for various purposes, especially in medical and agricultural applications. The main effect of CAP is associated with reactive oxygen and nitrogen species (RONS). In order to deliver these RONS to the target, direct or indirect treatment approaches have been employed. The indirect method is put into practice via plasma-activated water (PAW). Despite many studies being available in the field, the permeation mechanisms of RONS into water at the molecular level still remain elusive. Here, we performed molecular dynamics simulations to study the permeation of RONS from vacuum into the water interface and bulk. The calculated free energy profiles unravel the most favourable accumulation positions of RONS. Our results, therefore, provide fundamental insights into PAW and RONS chemistry to increase the efficiency of PAW in biological applications.

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.

Hypoglycemic and Anti-Inflammatory Effects of Triterpene Glycoside Fractions from Aeculus hippocastanum Seeds.

Horse chestnut (Aesculus hippocastanum L.)-derived drugs have shown their potential in biomedical applications. The seed of A. hippocastanum contains various kinds of chemical compounds including phenolics, flavonoids, coumarins, and triterpene saponins. Here, we investigated the chemical components in A. hippocastanum L. grown in Uzbekistan, which has not yet been studied in detail. We identified 30 kinds of triterpene saponins in an extract of A. hippocastanum L. Classifying extracted saponins into eight fractions, we next studied the hypoglycemic and the anti-inflammatory activities of escin and its derivatives through in vivo experiments. We came by data indicating the highest (SF-1 and SF-2) and the lowest (SF-5 and SF-8) antidiabetic and anti-inflammatory effects of those eight fractions. These results imply the prospective use of A. hippocastanum L. grown in Uzbekistan in the production of pharmaceutical drugs to treat diabetes and inflammation.

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.

Nanoscale insight into silk-like protein self-assembly: effect of design and number of repeat units.

By means of replica exchange molecular dynamics simulations we investigate how the length of a silk-like, alternating diblock oligopeptide influences its secondary and quaternary structure. We carry out simulations for two protein sizes consisting of three and five blocks, and study the stability of a single protein, a dimer, a trimer and a tetramer. Initial configurations of our simulations are ß-roll and ß-sheet structures. We find that for the triblock the secondary and quaternary structures upto and including the tetramer are unstable: the proteins melt into random coil structures and the aggregates disassemble either completely or partially. We attribute this to the competition between conformational entropy of the proteins and the formation of hydrogen bonds and hydrophobic interactions between proteins. This is confirmed by our simulations on the pentablock proteins, where we find that, as the number of monomers in the aggregate increases, individual monomers form more hydrogen bonds whereas their solvent accessible surface area decreases. For the pentablock ß-sheet protein, the monomer and the dimer melt as well, although for the ß-roll protein only the monomer melts. For both trimers and tetramers remain stable. Apparently, for these the entropy loss of forming ß-rolls and ß-sheets is compensated for in the free-energy gain due to the hydrogen-bonding and hydrophobic interactions. We also find that the middle monomers in the trimers and tetramers are conformationally much more stable than the ones on the top and the bottom. Interestingly, the latter are more stable on the tetramer than on the trimer, suggesting that as the number of monomers increases protein-protein interactions cooperatively stabilize the assembly. According to our simulations, the ß-roll and ß-sheet aggregates must be approximately equally stable.

Prediction of the structure of a silk-like protein in oligomeric states using explicit and implicit solvent models.

We perform Replica Exchange Molecular Dynamics (REMD) simulations on a silk-like protein design with amino-acid sequence [(Gly-Ala)3-Gly-Glu]5 to investigate the stability of a single protein, a dimer, a trimer and a tetramer made up of these proteins starting from ß-roll and ß-sheet structures in both explicit (TIP3P) and implicit (GBSA) solvent models. Our simulation results for the implicit solvent model agree with those for the explicit solvent model for simulation times up to the longest tested, being 30 ns per replica. From this we infer that the implicit solvent model that we use is reliable, allowing us to reach much longer time scales (up to 200 ns per replica). We find that the self-assembly of fibers of these proteins in solution must be a nucleated process, involving nuclei made up of at least three monomers. We also find that the conformation of the protein changes upon assembly, i.e., there is a transition from a disordered globular state to an ordered ß-sheet structure in the self-assembled state of aggregates containing more than two monomers. This indicates that autosteric effects must be important in the polymerization of this protein, reminiscent of what is observed for ß-amyloids. Our findings are consistent with recent experimental results on a protein with an amino acid sequence similar to that of the protein we study.

Structural and thermal analyses in semiconducting and metallic zigzag single-walled carbon nanotubes using molecular dynamics simulations.

Equilibrium molecular dynamics (EMD) simulations have been performed to investigate the structural analysis and thermal conductivity (λ) of semiconducting (8,0) and metallic (12,0) zigzag single-walled carbon nanotubes (SWCNTs) for varying ±Î³(%) strains. For the first time, the present outcomes provide valuable insights into the relationship between the structural properties of zigzag SWCNTs and corresponding thermal behavior, which is essential for the development of high-performance nanocomposites. The radial distribution function (RDF) has been employed to assess the buckling and deformation understandings of the (8,0) and (12,0) SWCNTs for a wide range of temperature T(K) and varying ±Î³(%) strains. The visualization of SWCNTs shows that the earlier buckling and deformation processes are observed for semiconducting SWCNTs as compared to metallic SWCNTs for high T(K) and it also evident through an abrupt increase in RDF peaks. The RDF and visualization analyses demonstrate that the (8,0) SWCNTs can more tunable under compressive than tensile strains, however, the (12,0) zigzag SWCNTs indicate an opposite trend and may tolerate more tensile than compressive strains. Investigations show that the tunable domain of ±Î³(%) strains decreases from (-10%≤ γ ≤+19%) to (-5%≤ γ ≤+10%) for (8,0) SWCNTs and the buckling process shifts to lower ±Î³(%) for (12,0) SWCNTs with increasing T(K). For intermediate-high T(K), the λ(T) of (12,0) SWCNTs is high but the (8,0) SWCNTs show certainly high λ(T) for low T(K). The present λ(T, ±Î³) data are in reasonable agreement with parts of previous NEMD, GK-HNEMD data and experimental investigations with simulation results generally under predicting the λ(T, ±Î³) by the ∼1% to ∼20%, regardless of the ±Î³(%) strains, depending on T(K). Our simulation data significantly expand the strain range to -10% ≤ γ ≤ +19% for both zigzag SWCNTs, depending on temperature T(K). This extension of the range aims to establish a tunable regime and delve into the intrinsic characteristics of zigzag SWCNTs, building upon previous work.

Unveiling the Borohydride Ion through Force-Field Development.

The borohydride ion, BH4-, is an essential reducing agent in many technological processes, yet its full understanding has been elusive, because of at least two significant challenges. One challenge arises from its marginal stability in aqueous solutions outside of basic pH conditions, which considerably limits the experimental thermodynamic data. The other challenge comes from its unique and atypical hydration shell, stemming from the negative excess charge on its hydrogen atoms, which complicates the accurate modeling in classical atomistic simulations. In this study, we combine experimental and computer simulation techniques to devise a classical force field for NaBH4 and deepen our understanding of its characteristics. We report the first measurement of the ion's activity coefficient and extrapolate it to neutral pH conditions. Given the difficulties in directly measuring its solvation free energies, owing to its instability, we resort to quantum chemistry calculations. This combined strategy allows us to derive a set of nonpolarizable force-field parameters for the borohydride ion for classical molecular dynamics simulations. The derived force field simultaneously captures the solvation free energy, the hydration structure, as well as the activity coefficient of NaBH4 salt across a broad concentration range. The obtained insights into the hydration shell of the BH4- ion are crucial for accurately modeling and understanding its interactions with other molecules, ions, materials, and interfaces.

Effect of calcination temperature induced structural modifications on the photocatalytic efficacy of Fe2O3-ZrO2 nanostructures: mechanochemical synthesis.

Water contamination due to organic pollutants is a challenging issue around the globe, and several attempts have been made to deal with this issue. Out of which, the semiconductor-based photocatalytic process had gained much attention and proved to be an efficient, easy, and economical process for the removal of organic dyes from aqueous solutions. For this purpose, the iron oxide-zirconium dioxide nanocomposite (Fe2O3-ZrO2 NC) was prepared via a simple mechanochemical process using a mortar and pestle, followed by a calcination process at 300, 600, and 900 °C. Different physicochemical analyses were carried out in order to investigate the successful synthesis of Fe2O3-ZrO2 NC and the effect of temperature on the crystallinity, surface area, pore size, phase composition, sample morphology, and particle/crystallite size. The Fe2O3-ZrO2 NCs were subjected to a photocatalytic test under solar light irradiation against fluorescein dye in an aqueous medium, and the photocatalytic performance was examined under the influence of calcination temperatures, pH, catalyst dose, and initial concentration. The stability of the Fe2O3-ZrO2 NCs was also checked by recycling them for five reuse cycles.

Escin's Multifaceted Therapeutic Profile in Treatment and Post-Treatment of Various Cancers: A Comprehensive Review.

Although modern medicine is advancing at an unprecedented rate, basic challenges in cancer treatment and drug resistance remain. Exploiting natural-product-based drugs is a strategy that has been proven over time to provide diverse and efficient approaches in patient care during treatment and post-treatment periods of various diseases, including cancer. Escin-a plant-derived triterpenoid saponin-is one example of natural products with a broad therapeutic scope. Initially, escin was proven to manifest potent anti-inflammatory and anti-oedematous effects. However, in the last two decades, other novel activities of escin relevant to cancer treatment have been reported. Recent studies demonstrated escin's efficacy in compositions with other approved drugs to accomplish synergy and increased bioavailability to broaden their apoptotic, anti-metastasis, and anti-angiogenetic effects. Here, we comprehensively discuss and present an overview of escin's chemistry and bioavailability, and highlight its biological activities against various cancer types. We conclude the review by presenting possible future directions of research involving escin for medical and pharmaceutical applications as well as for basic research.

Unraveling the Transport Properties of RONS across Nitro-Oxidized Membranes.

The potential of cold atmospheric plasma (CAP) in biomedical applications has received significant interest, due to its ability to generate reactive oxygen and nitrogen species (RONS). Upon exposure to living cells, CAP triggers alterations in various cellular components, such as the cell membrane. However, the permeation of RONS across nitrated and oxidized membranes remains understudied. To address this gap, we conducted molecular dynamics simulations, to investigate the permeation capabilities of RONS across modified cell membranes. This computational study investigated the translocation processes of less hydrophilic and hydrophilic RONS across the phospholipid bilayer (PLB), with various degrees of oxidation and nitration, and elucidated the impact of RONS on PLB permeability. The simulation results showed that less hydrophilic species, i.e., NO, NO2, N2O4, and O3, have a higher penetration ability through nitro-oxidized PLB compared to hydrophilic RONS, i.e., HNO3, s-cis-HONO, s-trans-HONO, H2O2, HO2, and OH. In particular, nitro-oxidation of PLB, induced by, e.g., cold atmospheric plasma, has minimal impact on the penetration of free energy barriers of less hydrophilic species, while it lowers these barriers for hydrophilic RONS, thereby enhancing their translocation across nitro-oxidized PLB. This research contributes to a better understanding of the translocation abilities of RONS in the field of plasma biomedical applications and highlights the need for further analysis of their role in intracellular signaling pathways.

Synthesis of nickel nanowires (Ni-NWs) as high ferromagnetic material by electrodeposition technique.

Metallic nanowires (NWs) and their different compounds display incredible prospects for their use in various applications including media storage, sensor and solar cell devices along with the biological drug delivery systems. In this research work, the metallic NWs like nickel nanowires (Ni-NWs) are synthesized successfully by employing electrodeposition process. Anodic aluminum oxide (AAO) templates are employed as a platform with copper metal coating which acts as an active cathode. The synthesized Ni-NWs are examined through various characterization techniques including X-ray diffraction (XRD), scanning electron microscope (SEM) and vibrating sample magnetometer (VSM) to study the crystal structure, surface morphology and magnetic properties, respectively. The XRD analysis shows the development of various diffraction planes like Ni (111), Ni (200), Ni (220) which confirms the formation of polycrystalline nickel NWs. The SEM analysis reveals that the range of diameter and length of nickel NWs are found to be ∼160 to 200 and ∼4 to 11 micron respectively showing high aspect ratio (ranged from ∼200 to 300). The ferromagnetic behavior of Ni-NWs is confirmed by the hysteresis loop carried out for parallel and perpendicular configurations having Hc = 100 and 206 Oe, respectively. The obtained results suggest that the synthesized Ni- NWs may be used for high-density media storage devices.

Bias-force guided simulations combined with experimental validations towards GPR17 modulators identification.

Glioblastoma Multiforme (GBM) is known to be by far the most aggressive brain tumor to affect adults. The median survival rate of GBM patient's is < 15 months, while the GBM cells aggressively develop resistance to chemo- and radiotherapy with their self-renewal capacity which suggests the pressing need to develop novel preventative measures. We have recently proved that GPR17 -an orphan G protein-coupled receptor- is highly expressed on the GBM cell surface and it has a vital role to play in the disease progression. Despite the progress made on GBM downregulation, there still remain difficulties in developing a promising modulator for GPR17, till date. Here, we have performed robust virtual screening combined with biased-force pulling molecular dynamic (MD) simulations to predict high-affinity GPR17 modulators followed by experimental validation. Initially, the database containing 1379 FDA-approved drugs were screened against the orthosteric binding pocket of the GPR17. The external bias-potentials were then applied to the screened hits during the MD simulations which enabled to predict a spectrum of rupture peak force values that were used to select four approved drugs -ZINC000003792417 (Sacubitril), ZINC000014210457 (Victrelis), ZINC000001536109 (Pralatrexate) and ZINC000003925861 (Vorapaxar)- as top hits. The hits selected turns out to demonstrate unique dissociation pathways, interaction pattern, and change in polar network over time. Subsequently the selected hits with GPR17 were measured by inhibiting the forskolin-stimulated cAMP accumulation in GBM cell lines, LN229 and SNB19. The ex vivo validations shows that Sacubitril drug can act as a full agonist, while Vorapaxar functions as a partial agonist for GPR17. The pEC50 of Sacubitril was identified as 4.841 and 4.661 for LN229 and SNB19, respectively. Small interference of the RNA (siRNA)- silenced the GPR17 to further validate the targeted binding of Sacubitril with GPR17. In the current investigation, we have identified new repurposable GPR17 specific drugs which are likely to increase the opportunity to treat orphan deadly diseases.

Formation of Highly Conductive Interfaces in Crystalline Ionic Liquid-Gated Unipolar MoTe2/h-BN Field-Effect Transistor.

2H MoTe2 (molybdenum ditelluride) has generated significant interest because of its superconducting, nonvolatile memory, and semiconducting of new materials, and it has a large range of electrical properties. The combination of transition metal dichalcogenides (TMDCs) and two dimensional (2D) materials like hexagonal boron nitride (h-BN) in lateral heterostructures offers a unique platform for designing and engineering novel electronic devices. We report the fabrication of highly conductive interfaces in crystalline ionic liquid-gated (ILG) field-effect transistors (FETs) consisting of a few layers of MoTe2/h-BN heterojunctions. In our initial exploration of tellurium-based semiconducting TMDs, we directed our attention to MoTe2 crystals with thicknesses exceeding 12 nm. Our primary focus centered on investigating the transport characteristics and quantitatively assessing the surface interface heterostructure. Our transconductance (gm) measurements indicate that the very efficient carrier modulation with an ILG FET is two times larger than standard back gating, and it demonstrates unipolarity of the device. The ILG FET exhibited highly unipolar p-type behavior with a high on/off ratio, and it significantly increased the mobility in MoTe2/h-BN heterochannels, achieving improvement as one of the highest recorded mobility increments. Specifically, we observed hole and electron mobility values ranging from 345 cm2 V-1 s-1 to 285 cm2 V-1 s-1 at 80 K. We predict that our ability to observe the intrinsic, heterointerface conduction in the channels was due to a drastic reduction of the Schottky barriers, and electrostatic gating is suggested as a method for controlling the phase transitions in the few layers of TMDC FETs. Moreover, the simultaneous structural phase transitions throughout the sample, achieved through electrostatic doping control, presents new opportunities for developing phase change devices using atomically thin membranes.