Genome-scale modeling predicts metabolic differences between macrophage subtypes in colorectal cancer

Dissecting metabolic landscape of alveolar macrophage

The highly plastic nature of Alveolar Macrophage (AM) plays a crucial role in the defense against inhaled particulates and pathogens in the lungs. Depending on the signal, AM acquires either the classically activated M1 phenotype or the alternatively activated M2 phenotype. In this study, we investigate the metabolic shift in the activated phases of AM (M1 and M2 phases) by reconstructing context specific Genome-Scale Metabolic (GSM) models. Metabolic pathways such as pyruvate metabolism, arachidonic acid metabolism, chondroitin/heparan sulfate biosynthesis, and heparan sulfate degradation are found to be important driving forces in the development of the M1/M2 phenotypes. Additionally, we formulated a bilevel optimization framework named MetaShiftOptimizer to identify minimal modifications that shift one activated state (M1/M2) to the other. The identified reactions involve metabolites such as glycogenin, L-carnitine, 5-hydroperoxy eicosatetraenoic acid, and leukotriene B4, which show potential to be further investigated as significant factors for developing efficient therapy targets for severe respiratory disorders in the future. Overall, our study contributes to the understanding of the metabolic capabilities of the M1 and M2 phenotype of AM and identifies pathways and reactions that can be potential targets for polarization shift and also be used as therapeutic strategies against respiratory diseases.

Integrating mechanism-based T cell phenotypes into a model of tumor-immune cell interactions

Interactions between cancer cells and immune cells in the tumor microenvironment influence tumor growth and can contribute to the response to cancer immunotherapies. It is difficult to gain mechanistic insights into the effects of cell-cell interactions in tumors using a purely experimental approach. However, computational modeling enables quantitative investigation of the tumor microenvironment, and agent-based modeling, in particular, provides relevant biological insights into the spatial and temporal evolution of tumors. Here, we develop a novel agent-based model (ABM) to predict the consequences of intercellular interactions. Furthermore, we leverage our prior work that predicts the transitions of CD8+ T cells from a naïve state to a terminally differentiated state using Boolean modeling. Given the details incorporated to predict T cell state, we apply the integrated Boolean-ABM framework to study how the properties of CD8+ T cells influence the composition and spatial organization of tumors and the efficacy of an immune checkpoint blockade. Overall, we present a mechanistic understanding of tumor evolution that can be leveraged to study targeted immunotherapeutic strategies.

Single-cell omics analysis with genome-scale metabolic modeling

Single-cell technologies have been widely used in biological studies and generated a plethora of single-cell data to be interpreted. Due to the inclusion of the priori metabolic network knowledge as well as gene-protein-reaction associations, genome-scale metabolic models (GEMs) have been a powerful tool to integrate and thereby interpret various omics data mostly from bulk samples. Here, we first review two common ways to leverage bulk omics data with GEMs and then discuss advances on integrative analysis of single-cell omics data with GEMs. We end by presenting our views on current challenges and perspectives in this field.