Genome-Scale Imaging of the 3D Organization and Transcriptional Activity of Chromatin.

The 3D organization of chromatin regulates many genome functions. Our understanding of 3D genome organization requires tools to directly visualize chromatin conformation in its native context. Here we report an imaging technology for visualizing chromatin organization across multiple scales in single cells with high genomic throughput. First we demonstrate multiplexed imaging of hundreds of genomic loci by sequential hybridization, which allows high-resolution conformation tracing of whole chromosomes. Next we report a multiplexed error-robust fluorescence in situ hybridization (MERFISH)-based method for genome-scale chromatin tracing and demonstrate simultaneous imaging of more than 1,000 genomic loci and nascent transcripts of more than 1,000 genes together with landmark nuclear structures. Using this technology, we characterize chromatin domains, compartments, and trans-chromosomal interactions and their relationship to transcription in single cells. We envision broad application of this high-throughput, multi-scale, and multi-modal imaging technology, which provides an integrated view of chromatin organization in its native structural and functional context.

Genome-Scale Identification of Essential Metabolic Processes for Targeting the Plasmodium Liver Stage.

Plasmodium gene functions in mosquito and liver stages remain poorly characterized due to limitations in the throughput of phenotyping at these stages. To fill this gap, we followed more than 1,300 barcoded P. berghei mutants through the life cycle. We discover 461 genes required for efficient parasite transmission to mosquitoes through the liver stage and back into the bloodstream of mice. We analyze the screen in the context of genomic, transcriptomic, and metabolomic data by building a thermodynamic model of P. berghei liver-stage metabolism, which shows a major reprogramming of parasite metabolism to achieve rapid growth in the liver. We identify seven metabolic subsystems that become essential at the liver stages compared with asexual blood stages: type II fatty acid synthesis and elongation (FAE), tricarboxylic acid, amino sugar, heme, lipoate, and shikimate metabolism. Selected predictions from the model are individually validated in single mutants to provide future targets for drug development.

Systems Biology of Metabolism.

Metabolism is highly complex and involves thousands of different connected reactions; it is therefore necessary to use mathematical models for holistic studies. The use of mathematical models in biology is referred to as systems biology. In this review, the principles of systems biology are described, and two different types of mathematical models used for studying metabolism are discussed: kinetic models and genome-scale metabolic models. The use of different omics technologies, including transcriptomics, proteomics, metabolomics, and fluxomics, for studying metabolism is presented. Finally, the application of systems biology for analyzing global regulatory structures, engineering the metabolism of cell factories, and analyzing human diseases is discussed.

A Multiplexed Single-Cell CRISPR Screening Platform Enables Systematic Dissection of the Unfolded Protein Response.

Functional genomics efforts face tradeoffs between number of perturbations examined and complexity of phenotypes measured. We bridge this gap with Perturb-seq, which combines droplet-based single-cell RNA-seq with a strategy for barcoding CRISPR-mediated perturbations, allowing many perturbations to be profiled in pooled format. We applied Perturb-seq to dissect the mammalian unfolded protein response (UPR) using single and combinatorial CRISPR perturbations. Two genome-scale CRISPR interference (CRISPRi) screens identified genes whose repression perturbs ER homeostasis. Subjecting ∼100 hits to Perturb-seq enabled high-precision functional clustering of genes. Single-cell analyses decoupled the three UPR branches, revealed bifurcated UPR branch activation among cells subject to the same perturbation, and uncovered differential activation of the branches across hits, including an isolated feedback loop between the translocon and IRE1α. These studies provide insight into how the three sensors of ER homeostasis monitor distinct types of stress and highlight the ability of Perturb-seq to dissect complex cellular responses.

Brain energy metabolism is optimized to minimize the cost of enzyme synthesis and transport.

The energy metabolism of the brain is poorly understood partly due to the complex morphology of neurons and fluctuations in ATP demand over time. To investigate this, we used metabolic models that estimate enzyme usage per pathway, enzyme utilization over time, and enzyme transportation to evaluate how these parameters and processes affect ATP costs for enzyme synthesis and transportation. Our models show that the total enzyme maintenance energy expenditure of the human body depends on how glycolysis and mitochondrial respiration are distributed both across and within cell types in the brain. We suggest that brain metabolism is optimized to minimize the ATP maintenance cost by distributing the different ATP generation pathways in an advantageous way across cell types and potentially also across synapses within the same cell. Our models support this hypothesis by predicting export of lactate from both neurons and astrocytes during peak ATP demand, reproducing results from experimental measurements reported in the literature. Furthermore, our models provide potential explanation for parts of the astrocyte-neuron lactate shuttle theory, which is recapitulated under some conditions in the brain, while contradicting other aspects of the theory. We conclude that enzyme usage per pathway, enzyme utilization over time, and enzyme transportation are important factors for defining the optimal distribution of ATP production pathways, opening a broad avenue to explore in brain metabolism.

Review and meta-analysis of the genetic Minimal Cut Set approach for gene essentiality prediction in cancer metabolism.

Cancer metabolism is a marvellously complex topic, in part, due to the reprogramming of its pathways to self-sustain the malignant phenotype in the disease, to the detriment of its healthy counterpart. Understanding these adjustments can provide novel targeted therapies that could disrupt and impair proliferation of cancerous cells. For this very purpose, genome-scale metabolic models (GEMs) have been developed, with Human1 being the most recent reconstruction of the human metabolism. Based on GEMs, we introduced the genetic Minimal Cut Set (gMCS) approach, an uncontextualized methodology that exploits the concepts of synthetic lethality to predict metabolic vulnerabilities in cancer. gMCSs define a set of genes whose knockout would render the cell unviable by disrupting an essential metabolic task in GEMs, thus, making cellular proliferation impossible. Here, we summarize the gMCS approach and review the current state of the methodology by performing a systematic meta-analysis based on two datasets of gene essentiality in cancer. First, we assess several thresholds and distinct methodologies for discerning highly and lowly expressed genes. Then, we address the premise that gMCSs of distinct length should have the same predictive power. Finally, we question the importance of a gene partaking in multiple gMCSs and analyze the importance of all the essential metabolic tasks defined in Human1. Our meta-analysis resulted in parameter evaluation to increase the predictive power for the gMCS approach, as well as a significant reduction of computation times by only selecting the crucial gMCS lengths, proposing the pertinency of particular parameters for the peak processing of gMCS.

Systems-level temporal immune-metabolic profile in Crimean-Congo hemorrhagic fever virus infection.

Crimean-Congo hemorrhagic fever (CCHF) caused by CCHF virus (CCHFV) is one of the epidemic-prone diseases prioritized by the World Health Organisation as public health emergency with an urgent need for accelerated research. The trajectory of host response against CCHFV is multifarious and remains unknown. Here, we reported the temporal spectrum of pathogenesis following the CCHFV infection using genome-wide blood transcriptomics analysis followed by advanced systems biology analysis, temporal immune-pathogenic alterations, and context-specific progressive and postinfection genome-scale metabolic models (GSMM) on samples collected during the acute (T0), early convalescent (T1), and convalescent-phase (T2). The interplay between the retinoic acid-inducible gene-I-like/nucleotide-binding oligomerization domain-like receptor and tumor necrosis factor signaling governed the trajectory of antiviral immune responses. The rearrangement of intracellular metabolic fluxes toward the amino acid metabolism and metabolic shift toward oxidative phosphorylation and fatty acid oxidation during acute CCHFV infection determine the pathogenicity. The upregulation of the tricarboxylic acid cycle during CCHFV infection, compared to the noninfected healthy control and between the severity groups, indicated an increased energy demand and cellular stress. The upregulation of glycolysis and pyruvate metabolism potentiated energy generation through alternative pathways associated with the severity of the infection. The downregulation of metabolic processes at the convalescent phase identified by blood cell transcriptomics and single-cell type proteomics of five immune cells (CD4+ and CD8+ T cells, CD14+ monocytes, B cells, and NK cells) potentially leads to metabolic rewiring through the recovery due to hyperactivity during the acute phase leading to post-viral fatigue syndrome.

DLTKcat: deep learning-based prediction of temperature-dependent enzyme turnover rates.

The enzyme turnover rate, ${k}_{cat}$, quantifies enzyme kinetics by indicating the maximum efficiency of enzyme catalysis. Despite its importance, ${k}_{cat}$ values remain scarce in databases for most organisms, primarily because of the cost of experimental measurements. To predict ${k}_{cat}$ and account for its strong temperature dependence, DLTKcat was developed in this study and demonstrated superior performance (log10-scale root mean squared error = 0.88, R-squared = 0.66) than previously published models. Through two case studies, DLTKcat showed its ability to predict the effects of protein sequence mutations and temperature changes on ${k}_{cat}$ values. Although its quantitative accuracy is not high enough yet to model the responses of cellular metabolism to temperature changes, DLTKcat has the potential to eventually become a computational tool to describe the temperature dependence of biological systems.

Applications of genome-scale metabolic models to investigate microbial metabolic adaptations in response to genetic or environmental perturbations.

Environmental perturbations are encountered by microorganisms regularly and will require metabolic adaptations to ensure an organism can survive in the newly presenting conditions. In order to study the mechanisms of metabolic adaptation in such conditions, various experimental and computational approaches have been used. Genome-scale metabolic models (GEMs) are one of the most powerful approaches to study metabolism, providing a platform to study the systems level adaptations of an organism to different environments which could otherwise be infeasible experimentally. In this review, we are describing the application of GEMs in understanding how microbes reprogram their metabolic system as a result of environmental variation. In particular, we provide the details of metabolic model reconstruction approaches, various algorithms and tools for model simulation, consequences of genetic perturbations, integration of '-omics' datasets for creating context-specific models and their application in studying metabolic adaptation due to the change in environmental conditions.

Genome-scale modeling drives 70-fold improvement of intracellular heme production in Saccharomyces cerevisiae.

Heme is an oxygen carrier and a cofactor of both industrial enzymes and food additives. The intracellular level of free heme is low, which limits the synthesis of heme proteins. Therefore, increasing heme synthesis allows an increased production of heme proteins. Using the genome-scale metabolic model (GEM) Yeast8 for the yeast Saccharomyces cerevisiae, we identified fluxes potentially important to heme synthesis. With this model, in silico simulations highlighted 84 gene targets for balancing biomass and increasing heme production. Of those identified, 76 genes were individually deleted or overexpressed in experiments. Empirically, 40 genes individually increased heme production (up to threefold). Heme was increased by modifying target genes, which not only included the genes involved in heme biosynthesis, but also those involved in glycolysis, pyruvate, Fe-S clusters, glycine, and succinyl-coenzyme A (CoA) metabolism. Next, we developed an algorithmic method for predicting an optimal combination of these genes by using the enzyme-constrained extension of the Yeast8 model, ecYeast8. The computationally identified combination for enhanced heme production was evaluated using the heme ligand-binding biosensor (Heme-LBB). The positive targets were combined using CRISPR-Cas9 in the yeast strain (IMX581-HEM15-HEM14-HEM3-Δshm1-HEM2-Δhmx1-FET4-Δgcv2-HEM1-Δgcv1-HEM13), which produces 70-fold-higher levels of intracellular heme.

Metastatic triple negative breast cancer adapts its metabolism to destination tissues while retaining key metabolic signatures.

Triple negative breast cancer (TNBC) metastases are assumed to exhibit similar functions in different organs as in the original primary tumor. However, studies of metastasis are often limited to a comparison of metastatic tumors with primary tumors of their origin, and little is known about the adaptation to the local environment of the metastatic sites. We therefore used transcriptomic data and metabolic network analyses to investigate whether metastatic tumors adapt their metabolism to the metastatic site and found that metastatic tumors adopt a metabolic signature with some similarity to primary tumors of their destinations. The extent of adaptation, however, varies across different organs, and metastatic tumors retain metabolic signatures associated with TNBC. Our findings suggest that a combination of anti-metastatic approaches and metabolic inhibitors selected specifically for different metastatic sites, rather than solely targeting TNBC primary tumors, may constitute a more effective treatment approach.

LooplessFluxSampler: an efficient toolbox for sampling the loopless flux solution space of metabolic models.

BACKGROUND: Uniform random sampling of mass-balanced flux solutions offers an unbiased appraisal of the capabilities of metabolic networks. Unfortunately, it is impossible to avoid thermodynamically infeasible loops in flux samples when using convex samplers on large metabolic models. Current strategies for randomly sampling the non-convex loopless flux space display limited efficiency and lack theoretical guarantees. RESULTS: Here, we present LooplessFluxSampler, an efficient algorithm for exploring the loopless mass-balanced flux solution space of metabolic models, based on an Adaptive Directions Sampling on a Box (ADSB) algorithm. ADSB is rooted in the general Adaptive Direction Sampling (ADS) framework, specifically the Parallel ADS, for which theoretical convergence and irreducibility results are available for sampling from arbitrary distributions. By sampling directions that adapt to the target distribution, ADSB traverses more efficiently the sample space achieving faster mixing than other methods. Importantly, the presented algorithm is guaranteed to target the uniform distribution over convex regions, and it provably converges on the latter distribution over more general (non-convex) regions provided the sample can have full support. CONCLUSIONS: LooplessFluxSampler enables scalable statistical inference of the loopless mass-balanced solution space of large metabolic models. Grounded in a theoretically sound framework, this toolbox provides not only efficient but also reliable results for exploring the properties of the almost surely non-convex loopless flux space. Finally, LooplessFluxSampler includes a Markov Chain diagnostics suite for assessing the quality of the final sample and the performance of the algorithm.

Flux sampling in genome-scale metabolic modeling of microbial communities.

BACKGROUND: Microbial communities play a crucial role in ecosystem function through metabolic interactions. Genome-scale modeling is a promising method to understand these interactions and identify strategies to optimize the community. Flux balance analysis (FBA) is most often used to predict the flux through all reactions in a genome-scale model; however, the fluxes predicted by FBA depend on a user-defined cellular objective. Flux sampling is an alternative to FBA, as it provides the range of fluxes possible within a microbial community. Furthermore, flux sampling can capture additional heterogeneity across a population, especially when cells exhibit sub-maximal growth rates. RESULTS: In this study, we simulate the metabolism of microbial communities and compare the metabolic characteristics found with FBA and flux sampling. With sampling, we find significant differences in the predicted metabolism, including an increase in cooperative interactions and pathway-specific changes in predicted flux. CONCLUSIONS: Our results suggest the importance of sampling-based approaches to evaluate metabolic interactions. Furthermore, we emphasize the utility of flux sampling in quantitatively studying interactions between cells and organisms.

Inferring secretory and metabolic pathway activity from omic data with secCellFie.

Understanding protein secretion has considerable importance in biotechnology and important implications in a broad range of normal and pathological conditions including development, immunology, and tissue function. While great progress has been made in studying individual proteins in the secretory pathway, measuring and quantifying mechanistic changes in the pathway's activity remains challenging due to the complexity of the biomolecular systems involved. Systems biology has begun to address this issue with the development of algorithmic tools for analyzing biological pathways; however most of these tools remain accessible only to experts in systems biology with extensive computational experience. Here, we expand upon the user-friendly CellFie tool which quantifies metabolic activity from omic data to include secretory pathway functions, allowing any scientist to infer properties of protein secretion from omic data. We demonstrate how the secretory expansion of CellFie (secCellFie) can help predict metabolic and secretory functions across diverse immune cells, hepatokine secretion in a cell model of NAFLD, and antibody production in Chinese Hamster Ovary cells.

Evaluation of enzyme-constrained genome-scale model through metabolic engineering of anaerobic co-production of 2,3-butanediol and glycerol by Saccharomyces cerevisiae.

Enzyme-constrained genome-scale models (ecGEMs) have potential to predict phenotypes in a variety of conditions, such as growth rates or carbon sources. This study investigated if ecGEMs can guide metabolic engineering efforts to swap anaerobic redox-neutral ATP-providing pathways in yeast from alcoholic fermentation to equimolar co-production of 2,3-butanediol and glycerol. With proven pathways and low product toxicity, the ecGEM solution space aligned well with observed phenotypes. Since this catabolic pathway provides only one-third of the ATP of alcoholic fermentation (2/3 versus 2 ATP per glucose), the ecGEM predicted a growth decrease from 0.36 h-1 in the reference to 0.175 h-1 in the engineered strain. However, this <3-fold decrease would require the specific glucose consumption rate to increase. Surprisingly, after the pathway swap the engineered strain immediately grew at 0.15 h-1 with a glucose consumption rate of 29 mmol (g CDW)-1 h-1, which was indeed higher than reference (23 mmol (g CDW)-1 h-1) and one of the highest reported for S. cerevisiae. The accompanying 2,3-butanediol- (15.8 mmol (g CDW)-1 h-1) and glycerol (19.6 mmol (g CDW)-1 h-1) production rates were close to predicted values. Proteomics confirmed that this increased consumption rate was facilitated by enzyme reallocation from especially ribosomes (from 25.5 to 18.5 %) towards glycolysis (from 28.7 to 43.5 %). Subsequently, 200 generations of sequential transfer did not improve growth of the engineered strain, showing the use of ecGEMs in predicting opportunity space for laboratory evolution. The observations in this study illustrate both the current potential, as well as future improvements, of ecGEMs as a tool for both metabolic engineering and laboratory evolution.

Machine learning predicts system-wide metabolic flux control in cyanobacteria.

Metabolic fluxes and their control mechanisms are fundamental in cellular metabolism, offering insights for the study of biological systems and biotechnological applications. However, quantitative and predictive understanding of controlling biochemical reactions in microbial cell factories, especially at the system level, is limited. In this work, we present ARCTICA, a computational framework that integrates constraint-based modelling with machine learning tools to address this challenge. Using the model cyanobacterium Synechocystis sp. PCC 6803 as chassis, we demonstrate that ARCTICA effectively simulates global-scale metabolic flux control. Key findings are that (i) the photosynthetic bioproduction is mainly governed by enzymes within the Calvin-Benson-Bassham (CBB) cycle, rather than by those involve in the biosynthesis of the end-product, (ii) the catalytic capacity of the CBB cycle limits the photosynthetic activity and downstream pathways and (iii) ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO) is a major, but not the most, limiting step within the CBB cycle. Predicted metabolic reactions qualitatively align with prior experimental observations, validating our modelling approach. ARCTICA serves as a valuable pipeline for understanding cellular physiology and predicting rate-limiting steps in genome-scale metabolic networks, and thus provides guidance for bioengineering of cyanobacteria.

Self-controlled in silico gene knockdown strategies to enhance the sustainable production of heterologous terpenoid by Saccharomyces cerevisiae.

Microbial bioengineering is a growing field for producing plant natural products (PNPs) in recent decades, using heterologous metabolic pathways in host cells. Once heterologous metabolic pathways have been introduced into host cells, traditional metabolic engineering techniques are employed to enhance the productivity and yield of PNP biosynthetic routes, as well as to manage competing pathways. The advent of computational biology has marked the beginning of a novel epoch in strain design through in silico methods. These methods utilize genome-scale metabolic models (GEMs) and flux optimization algorithms to facilitate rational design across the entire cellular metabolic network. However, the implementation of in silico strategies can often result in an uneven distribution of metabolic fluxes due to the rigid knocking out of endogenous genes, which can impede cell growth and ultimately impact the accumulation of target products. In this study, we creatively utilized synthetic biology to refine in silico strain design for efficient PNPs production. OptKnock simulation was performed on the GEM of Saccharomyces cerevisiae OA07, an engineered strain for oleanolic acid (OA) bioproduction that has been reported previously. The simulation predicted that the single deletion of fol1, fol2, fol3, abz1, and abz2, or a combined knockout of hfd1, ald2 and ald3 could improve its OA production. Consequently, strains EK1∼EK7 were constructed and cultivated. EK3 (OA07△fol3), EK5 (OA07△abz1), and EK6 (OA07△abz2) had significantly higher OA titers in a batch cultivation compared to the original strain OA07. However, these increases were less pronounced in the fed-batch mode, indicating that gene deletion did not support sustainable OA production. To address this, we designed a negative feedback circuit regulated by malonyl-CoA, a growth-associated intermediate whose synthesis served as a bypass to OA synthesis, at fol3, abz1, abz2, and at acetyl-CoA carboxylase-encoding gene acc1, to dynamically and autonomously regulate the expression of these genes in OA07. The constructed strains R_3A, R_5A and R_6A had significantly higher OA titers than the initial strain and the responding gene-knockout mutants in either batch or fed-batch culture modes. Among them, strain R_3A stand out with the highest OA titer reported to date. Its OA titer doubled that of the initial strain in the flask-level fed-batch cultivation, and achieved at 1.23 ± 0.04 g L-1 in 96 h in the fermenter-level fed-batch mode. This indicated that the integration of optimization algorithm and synthetic biology approaches was efficiently rational for PNP-producing strain design.

Systematic engineering for production of anti-aging sunscreen compound in Pseudomonas putida.

Sunscreen has been used for thousands of years to protect skin from ultraviolet radiation. However, the use of modern commercial sunscreen containing oxybenzone, ZnO, and TiO2 has raised concerns due to their negative effects on human health and the environment. In this study, we aim to establish an efficient microbial platform for production of shinorine, a UV light absorbing compound with anti-aging properties. First, we methodically selected an appropriate host for shinorine production by analyzing central carbon flux distribution data from prior studies alongside predictions from genome-scale metabolic models (GEMs). We enhanced shinorine productivity through CRISPRi-mediated downregulation and utilized shotgun proteomics to pinpoint potential competing pathways. Simultaneously, we improved the shinorine biosynthetic pathway by refining its design, optimizing promoter usage, and altering the strength of ribosome binding sites. Finally, we conducted amino acid feeding experiments under various conditions to identify the key limiting factors in shinorine production. The study combines meta-analysis of 13C-metabolic flux analysis, GEMs, synthetic biology, CRISPRi-mediated gene downregulation, and omics analysis to improve shinorine production, demonstrating the potential of Pseudomonas putida KT2440 as platform for shinorine production.

Elucidation of triacylglycerol catabolism in Yarrowia lipolytica: How cells balance acetyl-CoA and excess reducing equivalents.

Yarrowia lipolytica is an industrial yeast that can convert waste oil to value-added products. However, it is unclear how this yeast metabolizes lipid feedstocks, specifically triacylglycerol (TAG) substrates. This study used 13C-metabolic flux analysis (13C-MFA), genome-scale modeling, and transcriptomics analyses to investigate Y. lipolytica W29 growth with oleic acid, glycerol, and glucose. Transcriptomics data were used to guide 13C-MFA model construction and to validate the 13C-MFA results. The 13C-MFA data were then used to constrain a genome-scale model (GSM), which predicted Y. lipolytica fluxes, cofactor balance, and theoretical yields of terpene products. The three data sources provided new insights into cellular regulation during catabolism of glycerol and fatty acid components of TAG substrates, and how their consumption routes differ from glucose catabolism. We found that (1) over 80% of acetyl-CoA from oleic acid is processed through the glyoxylate shunt, a pathway that generates less CO2 compared to the TCA cycle, (2) the carnitine shuttle is a key regulator of the cytosolic acetyl-CoA pool in oleic acid and glycerol cultures, (3) the oxidative pentose phosphate pathway and mannitol cycle are key routes for NADPH generation, (4) the mannitol cycle and alternative oxidase activity help balance excess NADH generated from ß-oxidation of oleic acid, and (5) asymmetrical gene expressions and GSM simulations of enzyme usage suggest an increased metabolic burden for oleic acid catabolism.

Genome-scale and pathway engineering for the sustainable aviation fuel precursor isoprenol production in Pseudomonas putida.

Sustainable aviation fuel (SAF) will significantly impact global warming in the aviation sector, and important SAF targets are emerging. Isoprenol is a precursor for a promising SAF compound DMCO (1,4-dimethylcyclooctane) and has been produced in several engineered microorganisms. Recently, Pseudomonas putida has gained interest as a future host for isoprenol bioproduction as it can utilize carbon sources from inexpensive plant biomass. Here, we engineer metabolically versatile host P. putida for isoprenol production. We employ two computational modeling approaches (Bilevel optimization and Constrained Minimal Cut Sets) to predict gene knockout targets and optimize the "IPP-bypass" pathway in P. putida to maximize isoprenol production. Altogether, the highest isoprenol production titer from P. putida was achieved at 3.5 g/L under fed-batch conditions. This combination of computational modeling and strain engineering on P. putida for an advanced biofuels production has vital significance in enabling a bioproduction process that can use renewable carbon streams.