Constraint-Based Reconstruction and Analyses of Metabolic Models: Open-Source Python Tools and Applications to Cancer.

The influence of metabolism on signaling, epigenetic markers, and transcription is highly complex yet important for understanding cancer physiology. Despite the development of high-resolution multi-omics technologies, it is difficult to infer metabolic activity from these indirect measurements. Fortunately, genome-scale metabolic models and constraint-based modeling provide a systems biology framework to investigate the metabolic states and define the genotype-phenotype associations by integrations of multi-omics data. Constraint-Based Reconstruction and Analysis (COBRA) methods are used to build and simulate metabolic networks using mathematical representations of biochemical reactions, gene-protein reaction associations, and physiological and biochemical constraints. These methods have led to advancements in metabolic reconstruction, network analysis, perturbation studies as well as prediction of metabolic state. Most computational tools for performing these analyses are written for MATLAB, a proprietary software. In order to increase accessibility and handle more complex datasets and models, community efforts have started to develop similar open-source tools in Python. To date there is a comprehensive set of tools in Python to perform various flux analyses and visualizations; however, there are still missing algorithms in some key areas. This review summarizes the availability of Python software for several components of COBRA methods and their applications in cancer metabolism. These tools are evolving rapidly and should offer a readily accessible, versatile way to model the intricacies of cancer metabolism for identifying cancer-specific metabolic features that constitute potential drug targets.

Integrated analysis of plasma and single immune cells uncovers metabolic changes in individuals with COVID-19.

A better understanding of the metabolic alterations in immune cells during severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection may elucidate the wide diversity of clinical symptoms experienced by individuals with coronavirus disease 2019 (COVID-19). Here, we report the metabolic changes associated with the peripheral immune response of 198 individuals with COVID-19 through an integrated analysis of plasma metabolite and protein levels as well as single-cell multiomics analyses from serial blood draws collected during the first week after clinical diagnosis. We document the emergence of rare but metabolically dominant T cell subpopulations and find that increasing disease severity correlates with a bifurcation of monocytes into two metabolically distinct subsets. This integrated analysis reveals a robust interplay between plasma metabolites and cell-type-specific metabolic reprogramming networks that is associated with disease severity and could predict survival.

Multi-omic single-cell snapshots reveal multiple independent trajectories to drug tolerance in a melanoma cell line.

The determination of individual cell trajectories through a high-dimensional cell-state space is an outstanding challenge for understanding biological changes ranging from cellular differentiation to epigenetic responses of diseased cells upon drugging. We integrate experiments and theory to determine the trajectories that single BRAFV600E mutant melanoma cancer cells take between drug-naive and drug-tolerant states. Although single-cell omics tools can yield snapshots of the cell-state landscape, the determination of individual cell trajectories through that space can be confounded by stochastic cell-state switching. We assayed for a panel of signaling, phenotypic, and metabolic regulators at points across 5 days of drug treatment to uncover a cell-state landscape with two paths connecting drug-naive and drug-tolerant states. The trajectory a given cell takes depends upon the drug-naive level of a lineage-restricted transcription factor. Each trajectory exhibits unique druggable susceptibilities, thus updating the paradigm of adaptive resistance development in an isogenic cell population.

Deciphering UV-induced DNA Damage Responses to Prevent and Treat Skin Cancer.

Ultraviolet (UV) radiation is among the most prevalent environmental factors that influence human health and disease. Even 1 h of UV irradiation extensively damages the genome. To cope with resulting deleterious DNA lesions, cells activate a multitude of DNA damage response pathways, including DNA repair. Strikingly, UV-induced DNA damage formation and repair are affected by chromatin state. When cells enter S phase with these lesions, a distinct mutation signature is created via error-prone translesion synthesis. Chronic UV exposure leads to high mutation burden in skin and consequently the development of skin cancer, the most common cancer in the United States. Intriguingly, UV-induced oxidative stress has opposing effects on carcinogenesis. Elucidating the molecular mechanisms of UV-induced DNA damage responses will be useful for preventing and treating skin cancer with greater precision. Excitingly, recent studies have uncovered substantial depth of novel findings regarding the molecular and cellular consequences of UV irradiation. In this review, we will discuss updated mechanisms of UV-induced DNA damage responses including the ATR pathway, which maintains genome integrity following UV irradiation. We will also present current strategies for preventing and treating nonmelanoma skin cancer, including ATR pathway inhibition for prevention and photodynamic therapy for treatment.

Single-cell analysis resolves the cell state transition and signaling dynamics associated with melanoma drug-induced resistance.

Continuous BRAF inhibition of BRAF mutant melanomas triggers a series of cell state changes that lead to therapy resistance and escape from immune control before establishing acquired resistance genetically. We used genome-wide transcriptomics and single-cell phenotyping to explore the response kinetics to BRAF inhibition for a panel of patient-derived BRAFV600 -mutant melanoma cell lines. A subset of plastic cell lines, which followed a trajectory covering multiple known cell state transitions, provided models for more detailed biophysical investigations. Markov modeling revealed that the cell state transitions were reversible and mediated by both Lamarckian induction and nongenetic Darwinian selection of drug-tolerant states. Single-cell functional proteomics revealed activation of certain signaling networks shortly after BRAF inhibition, and before the appearance of drug-resistant phenotypes. Drug targeting those networks, in combination with BRAF inhibition, halted the adaptive transition and led to prolonged growth inhibition in multiple patient-derived cell lines.