Identification of high sugar diet-induced dysregulated metabolic pathways in muscle using tissue-specific metabolic models in Drosophila.

Individual tissues perform highly specialized metabolic functions to maintain whole-body homeostasis. Although Drosophila serves as a powerful model for studying human metabolic diseases, a lack of tissue-specific metabolic models makes it challenging to quantitatively assess the metabolic processes of individual tissues and disease models in this organism. To address this issue, we reconstructed 32 tissue-specific genome-scale metabolic models (GEMs) using pseudo-bulk single cell transcriptomics data, revealing distinct metabolic network structures across tissues. Leveraging enzyme kinetics and flux analyses, we predicted tissue-dependent metabolic pathway activities, recapitulating known tissue functions and identifying tissue-specific metabolic signatures, as supported by metabolite profiling. Moreover, to demonstrate the utility of tissue-specific GEMs in a disease context, we examined the effect of a high sugar diet (HSD) on muscle metabolism. Together with 13C-glucose isotopic tracer studies, we identified glyceraldehyde 3-phosphate dehydrogenase (GAPDH) as a rate-limiting enzyme in response to HSD. Mechanistically, the decreased GAPDH activity was linked to elevated NADH/NAD+ ratio, caused by disturbed NAD+ regeneration rates, and oxidation of GAPDH. Furthermore, we introduced a pathway flux index to predict and validate additionally perturbed pathways, including fructose and butanoate metabolism. Altogether, our results represent a significant advance in generating quantitative tissue-specific GEMs and flux analyses in Drosophila, highlighting their use for identifying dysregulated metabolic pathways and their regulation in a human disease model.

NADPH composite index analysis quantifies the relationship between compartmentalized NADPH dynamics and growth rates in cancer cells.

NADPH, a highly compartmentalized electron donor in mammalian cells, plays essential roles in cell metabolism. However, little is known about how cytosolic and mitochondrial NADPH dynamics relate to cancer cell growth rates in response to varying nutrient conditions. To address this issue, we present NADPH composite index analysis, which quantifies the relationship between compartmentalized NADPH dynamics and growth rates using genetically encoded NADPH sensors, automated image analysis pipeline, and correlation analysis. Through this analysis, we demonstrated that compartmentalized NADPH dynamics patterns were cancer cell-type dependent. Specifically, cytosolic and mitochondrial NADPH dynamics of MDA-MB-231 decreased in response to serine deprivation, while those of HCT-116 increased in response to serine or glutamine deprivation. Furthermore, by introducing a fractional contribution parameter, we correlated cytosolic and mitochondrial NADPH dynamics to growth rates. Using this parameter, we identified cancer cell lines whose growth rates were selectively inhibited by targeting cytosolic or mitochondrial NADPH metabolism. Mechanistically, we identified citrate transporter as a key mitochondrial transporter that maintains compartmentalized NADPH dynamics and growth rates. Altogether, our results present a significant advance in quantifying the relationship between compartmentalized NADPH dynamics and cancer cell growth rates, highlighting a potential of targeting compartmentalized NADPH metabolism for selective cancer cell growth inhibitions.

Screening compound libraries for H2O2-mediated cancer therapeutics using a peroxiredoxin-based sensor.

Compounds that modulate H2O2 reaction networks have applications as targeted cancer therapeutics, as a subset of cancers exhibit sensitivity to this redox signal. Previous studies to identify therapeutics that induce oxidants have relied upon probes that respond to many different oxidants in cells, and thus do not report on only H2O2, a redox signal that selectively oxidizes proteins. Here we use a genetically encoded fluorescent probe for human peroxiredoxin-2 (Prx2) oxidation in screens for small-molecule compounds that modulate H2O2 pathways. We further characterize cellular responses to several compounds selected from the screen. Our results reveal that some, but not all, of the compounds enact H2O2-mediated toxicity in cells. Among them, SMER3, an antifungal, has not been reported as an oxidant-inducing drug. Several drugs, including cisplatin, that previously have been shown to induce reactive oxygen species (ROS) do not appear to oxidize Prx2, suggesting H2O2 is not among the ROS induced by those drugs.

Nanosensor Chemical Cytometry for Characterizing the Efflux Heterogeneity of Nitric Oxide from Macrophages.

Macrophages are a critical part of the human immune response, and their collective heterogeneity is implicated in disease progression and prevention. A nondestructive, label-free tool does not currently exist for profiling the dynamic, antigenic responses of single macrophages in a collection to correlate with specific molecular expression and correlated biophysical properties at the cellular level, despite the potential for diagnosis and therapeutics. Herein, we develop a nanosensor chemical cytometry (NCC) that can profile the heterogeneity of inducible nitric oxide synthase (iNOS) responses from macrophage populations. By integrating a near-infrared (nIR) fluorescent nanosensor array and collagen layer with microfluidics, the cellular lensing effect of the macrophage was utilized to characterize both nitric oxide (NO) efflux and refractive index (RI) changes at a single-cell level. Using a parallel, multichannel approach, distinct iNOS heterogeneities of macrophages can be monitored at an attomolar (10-18 mol) sensitivity in a nondestructive and real-time manner with a throughput of exceeding the 200 cells/frame. We demonstrate that estimated mean NO efflux rates of macrophage populations are elevated from 342 (σ = 199) to 464 (σ = 206) attomol/cell·hr with a 3% larger increase in the heterogeneity, and estimated RI of macrophage decrease from 1.366 (σ = 0.015) to 1.359 (σ = 0.009) with trimodal subpopulations under lipopolysaccharide (LPS) activation. These measured values are also in good agreement with Griess assay results and previously reported measurements. This work provides an efficient strategy for single-cell analysis of macrophage populations for cellular manufacturing and biopharmaceutical engineering.

Cellular lensing and near infrared fluorescent nanosensor arrays to enable chemical efflux cytometry.

Nanosensors have proven to be powerful tools to monitor single cells, achieving spatiotemporal precision even at molecular level. However, there has not been way of extending this approach to statistically relevant numbers of living cells. Herein, we design and fabricate nanosensor array in microfluidics that addresses this limitation, creating a Nanosensor Chemical Cytometry (NCC). nIR fluorescent carbon nanotube array is integrated along microfluidic channel through which flowing cells is guided. We can utilize the flowing cell itself as highly informative Gaussian lenses projecting nIR profiles and extract rich information. This unique biophotonic waveguide allows for quantified cross-correlation of biomolecular information with various physical properties and creates label-free chemical cytometer for cellular heterogeneity measurement. As an example, the NCC can profile the immune heterogeneities of human monocyte populations at attomolar sensitivity in completely non-destructive and real-time manner with rate of ~600 cells/hr, highest range demonstrated to date for state-of-the-art chemical cytometry.

Differential substrate use in EGF- and oncogenic KRAS-stimulated human mammary epithelial cells.

Many metabolic phenotypes in cancer cells are also characteristic of proliferating nontransformed mammalian cells, and attempts to distinguish between phenotypes resulting from oncogenic perturbation from those associated with increased proliferation are limited. Here, we examined the extent to which metabolic changes corresponding to oncogenic KRAS expression differed from those corresponding to epidermal growth factor (EGF)-driven proliferation in human mammary epithelial cells (HMECs). Removal of EGF from culture medium reduced growth rates and glucose/glutamine consumption in control HMECs despite limited changes in respiration and fatty acid synthesis, while the relative contribution of branched-chain amino acids to the TCA cycle and lipogenesis increased in the near-quiescent conditions. Most metabolic phenotypes measured in HMECs expressing mutant KRAS were similar to those observed in EGF-stimulated control HMECs that were growing at comparable rates. However, glucose and glutamine consumption as well as lactate and glutamate production were lower in KRAS-expressing cells cultured in media without added EGF, and these changes correlated with reduced sensitivity to GLUT1 inhibitor and phenformin treatment. Our results demonstrate the strong dependence of metabolic behavior on growth rate and provide a model to distinguish the metabolic influences of oncogenic mutations and nononcogenic growth.

Oxidative pentose phosphate pathway and glucose anaplerosis support maintenance of mitochondrial NADPH pool under mitochondrial oxidative stress.

Mitochondrial NADPH protects cells against mitochondrial oxidative stress by serving as an electron donor to antioxidant defense systems. However, due to technical challenges, it still remains unknown as to the pool size of mitochondrial NADPH, its dynamics, and NADPH/NADP+ ratio. Here, we have systemically modulated production rates of H2O2 in mitochondria and assessed mitochondrial NADPH metabolism using iNap sensors, 13C glucose isotopic tracers, and a mathematical model. Using sensors, we observed decreases in mitochondrial NADPH caused by excessive generation of mitochondrial H2O2, whereas the cytosolic NADPH was maintained upon perturbation. We further quantified the extent of mitochondrial NADPH/NADP+ based on the mathematical analysis. Utilizing 13C glucose isotopic tracers, we found increased activity in the pentose phosphate pathway (PPP) accompanied small decreases in the mitochondrial NADPH pool, whereas larger decreases induced both PPP activity and glucose anaplerosis. Thus, our integrative and quantitative approach provides insight into mitochondrial NADPH metabolism during mitochondrial oxidative stress.

Kinetic modeling of H2O2 dynamics in the mitochondria of HeLa cells.

Hydrogen peroxide (H2O2) promotes a range of phenotypes depending on its intracellular concentration and dosing kinetics, including cell death. While this qualitative relationship has been well established, the quantitative and mechanistic aspects of H2O2 signaling are still being elucidated. Mitochondria, a putative source of intracellular H2O2, have recently been demonstrated to be particularly vulnerable to localized H2O2 perturbations, eliciting a dramatic cell death response in comparison to similar cytosolic perturbations. We sought to improve our dynamic and mechanistic understanding of the mitochondrial H2O2 reaction network in HeLa cells by creating a kinetic model of this system and using it to explore basal and perturbed conditions. The model uses the most current quantitative proteomic and kinetic data available to predict reaction rates and steady-state concentrations of H2O2 and its reaction partners within individual mitochondria. Time scales ranging from milliseconds to one hour were simulated. We predict that basal, steady-state mitochondrial H2O2 will be in the low nM range (2-4 nM) and will be inversely dependent on the total pool of peroxiredoxin-3 (Prx3). Neglecting efflux of H2O2 to the cytosol, the mitochondrial reaction network is expected to control perturbations well up to H2O2 generation rates ~50 µM/s (0.25 nmol/mg-protein/s), above which point the Prx3 system would be expected to collapse. Comparison of these results with redox Western blots of Prx3 and Prx2 oxidation states demonstrated reasonable trend agreement at short times (≤ 15 min) for a range of experimentally perturbed H2O2 generation rates. At longer times, substantial efflux of H2O2 from the mitochondria to the cytosol was evidenced by peroxiredoxin-2 (Prx2) oxidation, and Prx3 collapse was not observed. A refined model using Monte Carlo parameter sampling was used to explore rates of H2O2 efflux that could reconcile model predictions of Prx3 oxidation states with the experimental observations.

A xenograft and cell line model of SDH-deficient pheochromocytoma derived from Sdhb+/- rats.

Tumors caused by loss-of-function mutations in genes encoding TCA cycle enzymes have been recently discovered and are now of great interest. Mutations in succinate dehydrogenase (SDH) subunits cause pheochromocytoma/paraganglioma (PCPG) and syndromically associated tumors, which differ phenotypically and clinically from more common SDH-intact tumors of the same types. Consequences of SDH deficiency include rewired metabolism, pseudohypoxic signaling and altered redox balance. PCPG with SDHB mutations are particularly aggressive, and development of treatments has been hampered by lack of valid experimental models. Attempts to develop mouse models have been unsuccessful. Using a new strategy, we developed a xenograft and cell line model of SDH-deficient pheochromocytoma from rats with a heterozygous germline Sdhb mutation. The genome, transcriptome and metabolome of this model, called RS0, closely resemble those of SDHB-mutated human PCPGs, making it the most valid model now available. Strategies employed to develop RS0 may be broadly applicable to other SDH-deficient tumors.

Monitoring the action of redox-directed cancer therapeutics using a human peroxiredoxin-2-based probe.

Redox cancer therapeutics target the increased reliance on intracellular antioxidant systems and enhanced susceptibility to oxidant-induced stress of some cancer cells compared to normal cells. Many of these therapeutics are thought to perturb intracellular levels of the oxidant hydrogen peroxide (H2O2), a signaling molecule that modulates a number of different processes in human cells. However, fluorescent probes for this species remain limited in their ability to detect the small perturbations induced during successful treatments. We report a fluorescent sensor based upon human peroxiredoxin-2, which acts as the natural indicator of small H2O2 fluctuations in human cells. The new probe reveals peroxide-induced oxidation in human cells below the detection limit of current probes, as well as peroxiredoxin-2 oxidation caused by two different redox cancer therapeutics in living cells. This capability will be useful in elucidating the mechanism of current redox-based therapeutics and in developing new ones.

Mitochondrial H2O2 Generation Using a Tunable Chemogenetic Tool To Perturb Redox Homeostasis in Human Cells and Induce Cell Death.

Among reactive oxygen species (ROS), H2O2 alone acts as a signaling molecule that promotes diverse phenotypes depending on the intracellular concentration. Mitochondria have been suggested as both sources and sinks of cellular H2O2, and mitochondrial dysfunction has been implicated in diseases such as cancer. A genetically encoded H2O2 generator, d-amino acid oxidase (DAAO), was targeted to the mitochondria of human cells, and its utility in investigating cellular response to a range of H2O2 doses over time was assessed. Organelle-specific peroxiredoxin dimerization and protein S-glutathionylation were measured as indicators of increased H2O2 flux due to the activity of DAAO. Cell death was observed in a concentration- and time-dependent manner, and protein oxidation shifted in localization as the dose increased. This work presents the first systematic study of H2O2-specific perturbation of mitochondria in human cells, and it reveals a marked sensitivity of this organelle to increases in H2O2 in comparison with prior studies that targeted the cytosol.