Merging Metabolic Modeling and Imaging for Screening Therapeutic Targets in Colorectal Cancer.

Cancer-associated fibroblasts (CAFs) play a key role in metabolic reprogramming and are well-established contributors to drug resistance in colorectal cancer (CRC). To exploit this metabolic crosstalk, we integrated a systems biology approach that identified key metabolic targets in a data-driven method and validated them experimentally. This process involved high-throughput computational screening to investigate the effects of enzyme perturbations predicted by a computational model of CRC metabolism to understand system-wide effects efficiently. Our results highlighted hexokinase (HK) as one of the crucial targets, which subsequently became our focus for experimental validation using patient-derived tumor organoids (PDTOs). Through metabolic imaging and viability assays, we found that PDTOs cultured in CAF conditioned media exhibited increased sensitivity to HK inhibition. Our approach emphasizes the critical role of integrating computational and experimental techniques in exploring and exploiting CRC-CAF crosstalk.

Modeling the heterogeneous apoptotic response of caspase-mediated signaling in tumor cells.

Resisting apoptosis is a hallmark of cancer. For this reason, it may be possible to force cancer cells to die by targeting components along the apoptotic signaling pathway. However, apoptosis signaling is challenging to understand due to dynamic and complex behaviors of ligands, receptors, and intracellular signaling components in response to cancer therapy. In this work, we forecast the apoptotic response based on the combined impact of these features. We expanded a previously established mathematical model of caspase-mediated apoptosis to include extracellular activation and receptor dynamics. In addition, three potential threshold values of caspase-3 necessary for the activation of apoptosis were selected to forecast which cells become apoptotic over time. We first vary ligand and receptor levels with the number of intracellular signaling proteins remaining consistent. Then, we vary the intracellular protein molecules in each simulated tumor cell to forecast the response of a heterogeneous population. By leveraging the benefits of computational modeling, we investigate the combined effect of several factors on the onset of apoptosis. This work provides quantitative insights for how the apoptotic signaling response can be forecasted, and precisely triggered, amongst heterogeneous cells via extracellular activation.

Integrating Diversity, Equity, and Inclusion into Preclinical, Clinical, and Public Health Mathematical Models.

Mathematical modelling applied to preclinical, clinical, and public health research is critical for our understanding of a multitude of biological principles. Biology is fundamentally heterogeneous, and mathematical modelling must meet the challenge of variability head on to ensure the principles of diversity, equity, and inclusion (DEI) are integrated into quantitative analyses. Here we provide a follow-up perspective on the DEI plenary session held at the 2023 Society for Mathematical Biology Annual Meeting to discuss key issues for the increased integration of DEI in mathematical modelling in biology.

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.

Information-Theoretic Analysis of a Model of CAR-4-1BB-Mediated NFκB Activation.

Systems biology utilizes computational approaches to examine an array of biological processes, such as cell signaling, metabolomics and pharmacology. This includes mathematical modeling of CAR T cells, a modality of cancer therapy by which genetically engineered immune cells recognize and combat a cancerous target. While successful against hematologic malignancies, CAR T cells have shown limited success against other cancer types. Thus, more research is needed to understand their mechanisms of action and leverage their full potential. In our work, we set out to apply information theory on a mathematical model of NFκB signaling initiated by the CAR following antigen encounter. First, we estimated channel capacity for CAR-4-1BB-mediated NFκB signal transduction. Next, we evaluated the pathway's ability to distinguish contrasting "low" and "high" antigen concentration levels, depending on the amount of variability in protein concentrations. Finally, we assessed the fidelity by which NFκB activation reflects the encountered antigen concentration, depending on the prevalence of antigen-positive targets in tumor population. We found that in most scenarios, fold change in the nuclear concentration of NFκB carries a higher channel capacity for the pathway than NFκB's absolute response. Additionally, we found that most errors in transducing the antigen signal through the pathway skew towards underestimating the concentration of encountered antigen. Finally, we found that disabling IKKß deactivation could increase signaling fidelity against targets with antigen-negative cells. Our information-theoretic analysis of signal transduction can provide novel perspectives on biological signaling, as well as enable a more informed path to cell engineering.Kindly check and confirm whether the corresponding affiliation is correctly identified.this is correct.

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

Colorectal cancer (CRC) shows high incidence and mortality, partly due to the tumor microenvironment (TME), which is viewed as an active promoter of disease progression. Macrophages are among the most abundant cells in the TME. These immune cells are generally categorized as M1, with inflammatory and anti-cancer properties, or M2, which promote tumor proliferation and survival. Although the M1/M2 subclassification scheme is strongly influenced by metabolism, the metabolic divergence between the subtypes remains poorly understood. Therefore, we generated a suite of computational models that characterize the M1- and M2-specific metabolic states. Our models show key differences between the M1 and M2 metabolic networks and capabilities. We leverage the models to identify metabolic perturbations that cause the metabolic state of M2 macrophages to more closely resemble M1 cells. Overall, this work increases understanding of macrophage metabolism in CRC and elucidates strategies to promote the metabolic state of anti-tumor macrophages.

A data-driven Boolean model explains memory subsets and evolution in CD8+ T cell exhaustion.

T cells play a key role in a variety of immune responses, including infection and cancer. Upon stimulation, naïve CD8+ T cells proliferate and differentiate into a variety of memory and effector cell types; however, failure to clear antigens causes prolonged stimulation of CD8+ T cells, ultimately leading to T cell exhaustion (TCE). The functional and phenotypic changes that occur during CD8+ T cell differentiation are well characterized, but the underlying gene expression state changes are not completely understood. Here, we utilize a previously published data-driven Boolean model of gene regulatory interactions shown to mediate TCE. Our network analysis and modeling reveal the final gene expression states that correspond to TCE, along with the sequence of gene expression patterns that give rise to those final states. With a model that predicts the changes in gene expression that lead to TCE, we could evaluate strategies to inhibit the exhausted state. Overall, we demonstrate that a common pathway model of CD8+ T cell gene regulatory interactions can provide insights into the transcriptional changes underlying the evolution of cell states in TCE.

Information-theoretic analysis of a model of CAR-4-1BB-mediated NFκB activation.

Systems biology utilizes computational approaches to examine an array of biological processes, such as cell signaling, metabolomics and pharmacology. This includes mathematical modeling of CAR T cells, a modality of cancer therapy by which genetically engineered immune cells recognize and combat a cancerous target. While successful against hematologic malignancies, CAR T cells have shown limited success against other cancer types. Thus, more research is needed to understand their mechanisms of action and leverage their full potential. In our work, we set out to apply information theory on a mathematical model of cell signaling of CAR-mediated activation following antigen encounter. First, we estimated channel capacity for CAR-4-1BB-mediated NFκB signal transduction. Next, we evaluated the pathway's ability to distinguish contrasting "low" and "high" antigen concentration levels, depending on the amount of intrinsic noise. Finally, we assessed the fidelity by which NFκB activation reflects the encountered antigen concentration, depending on the prevalence of antigen-positive targets in tumor population. We found that in most scenarios, fold change in the nuclear concentration of NFκB carries a higher channel capacity for the pathway than NFκB's absolute response. Additionally, we found that most errors in transducing the antigen signal through the pathway skew towards underestimating the concentration of encountered antigen. Finally, we found that disabling IKKß deactivation could increase signaling fidelity against targets with antigen-negative cells. Our information-theoretic analysis of signal transduction can provide novel perspectives on biological signaling, as well as enable a more informed path to cell engineering.

A systems and computational biology perspective on advancing CAR therapy.

In the recent decades, chimeric antigen receptor (CAR) therapy signaled a new revolutionary approach to cancer treatment. This method seeks to engineer immune cells expressing an artificially designed receptor, which would endue those cells with the ability to recognize and eliminate tumor cells. While some CAR therapies received FDA approval and others are subject to clinical trials, many aspects of their workings remain elusive. Techniques of systems and computational biology have been frequently employed to explain the operating principles of CAR therapy and suggest further design improvements. In this review, we sought to provide a comprehensive account of those efforts. Specifically, we discuss various computational models of CAR therapy ranging in scale from organismal to molecular. Then, we describe the molecular and functional properties of costimulatory domains frequently incorporated in CAR structure. Finally, we describe the signaling cascades by which those costimulatory domains elicit cellular response against the target. We hope that this comprehensive summary of computational and experimental studies will further motivate the use of systems approaches in advancing CAR therapy.

Flux Sampling in Genome-scale Metabolic Modeling of Microbial Communities.

Microbial communities play a crucial role in ecosystem function through metabolic interactions. Genome-scale modeling is a promising method to understand these interactions. 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 may capture additional heterogeneity across cells, especially when cells exhibit sub-maximal growth rates. In this study, we simulate the metabolism of microbial communities and compare the metabolic characteristics found with FBA and flux sampling. We find significant differences in the predicted metabolism with sampling, including increased cooperative interactions and pathway-specific changes in predicted flux. Our results suggest the importance of sampling-based and objective function-independent approaches to evaluate metabolic interactions and emphasize their utility in quantitatively studying interactions between cells and organisms.

Calibrating agent-based models to tumor images using representation learning.

Agent-based models (ABMs) have enabled great advances in the study of tumor development and therapeutic response, allowing researchers to explore the spatiotemporal evolution of the tumor and its microenvironment. However, these models face serious drawbacks in the realm of parameterization - ABM parameters are typically set individually based on various data and literature sources, rather than through a rigorous parameter estimation approach. While ABMs can be fit to simple time-course data (such as tumor volume), that type of data loses the spatial information that is a defining feature of ABMs. While tumor images provide spatial information, it is exceedingly difficult to compare tumor images to ABM simulations beyond a qualitative visual comparison. Without a quantitative method of comparing the similarity of tumor images to ABM simulations, a rigorous parameter fitting is not possible. Here, we present a novel approach that applies neural networks to represent both tumor images and ABM simulations as low dimensional points, with the distance between points acting as a quantitative measure of difference between the two. This enables a quantitative comparison of tumor images and ABM simulations, where the distance between simulated and experimental images can be minimized using standard parameter-fitting algorithms. Here, we describe this method and present two examples to demonstrate the application of the approach to estimate parameters for two distinct ABMs. Overall, we provide a novel method to robustly estimate ABM parameters.

Ensemble-based genome-scale modeling predicts metabolic differences between macrophage subtypes in colorectal cancer.

1Colorectal cancer (CRC) shows high incidence and mortality, partly due to the tumor microenvironment, which is viewed as an active promoter of disease progression. Macrophages are among the most abundant cells in the tumor microenvironment. These immune cells are generally categorized as M1, with inflammatory and anti-cancer properties, or M2, which promote tumor proliferation and survival. Although the M1/M2 subclassification scheme is strongly influenced by metabolism, the metabolic divergence between the subtypes remains poorly understood. Therefore, we generated a suite of computational models that characterize the M1- and M2-specific metabolic states. Our models show key differences between the M1 and M2 metabolic networks and capabilities. We leverage the models to identify metabolic perturbations that cause the metabolic state of M2 macrophages to more closely resemble M1 cells. Overall, this work increases understanding of macrophage metabolism in CRC and elucidates strategies to promote the metabolic state of anti-tumor macrophages.

Proof-of-concept for incorporating mechanistic insights from multi-omics analyses of polymyxin B in combination with chloramphenicol against Klebsiella pneumoniae.

Carbapenemase-resistant Klebsiella pneumoniae (KP) resistant to multiple antibiotic classes necessitates optimized combination therapy. Our objective is to build a workflow leveraging omics and bacterial count data to identify antibiotic mechanisms that can be used to design and optimize combination regimens. For pharmacodynamic (PD) analysis, previously published static time-kill studies (J Antimicrob Chemother 70, 2015, 2589) were used with polymyxin B (PMB) and chloramphenicol (CHL) mono and combination therapy against three KP clinical isolates over 24 h. A mechanism-based model (MBM) was developed using time-kill data in S-ADAPT describing PMB-CHL PD activity against each isolate. Previously published results of PMB (1 mg/L continuous infusion) and CHL (Cmax : 8 mg/L; bolus q6h) mono and combination regimens were evaluated using an in vitro one-compartment dynamic infection model against a KP clinical isolate (108 CFU/ml inoculum) over 24 h to obtain bacterial samples for multi-omics analyses. The differentially expressed genes and metabolites in these bacterial samples served as input to develop a partial least squares regression (PLSR) in R that links PD responses with the multi-omics responses via a multi-omics pathway analysis. PMB efficacy was increased when combined with CHL, and the MBM described the observed PD well for all strains. The PLSR consisted of 29 omics inputs and predicted MBM PD response (R2  = 0.946). Our analysis found that CHL downregulated metabolites and genes pertinent to lipid A, hence limiting the emergence of PMB resistance. Our workflow linked insights from analysis of multi-omics data with MBM to identify biological mechanisms explaining observed PD activity in combination therapy.

Systematic Bayesian posterior analysis guided by Kullback-Leibler divergence facilitates hypothesis formation.

Bayesian inference produces a posterior distribution for the parameters of a mathematical model that can be used to guide the formation of hypotheses; specifically, the posterior may be searched for evidence of alternative model hypotheses, which serves as a starting point for hypothesis formation and model refinement. Previous approaches to search for this evidence are largely qualitative and unsystematic; further, demonstrations of these approaches typically stop at hypothesis formation, leaving the questions they raise unanswered. Here, we introduce a Kullback-Leibler (KL) divergence-based ranking to expedite Bayesian hypothesis formation and investigate the hypotheses it generates, ultimately generating novel, biologically significant insights. Our approach uses KL divergence to rank parameters by how much information they gain from experimental data. Subsequently, rather than searching all model parameters at random, we use this ranking to prioritize examining the posteriors of the parameters that gained the most information from the data for evidence of alternative model hypotheses. We test our approach with two examples, which showcase the ability of our approach to systematically uncover different types of alternative hypothesis evidence. First, we test our KL divergence ranking on an established example of Bayesian hypothesis formation. Our top-ranked parameter matches the one previously identified to produce alternative hypotheses. In the second example, we apply our ranking in a novel study of a computational model of prolactin-induced JAK2-STAT5 signaling, a pathway that mediates beta cell proliferation. Within the top 3 ranked parameters (out of 33), we find a bimodal posterior revealing two possible ranges for the prolactin receptor degradation rate. We go on to refine the model, incorporating new data and determining which degradation rate is most plausible. Overall, while the effectiveness of our approach depends on having a properly formulated prior and on the form of the posterior distribution, we demonstrate that our approach offers a novel and generalizable quantitative framework for Bayesian hypothesis formation and use it to produce a novel, biologically-significant insight into beta cell signaling.

Kinetic and data-driven modeling of pancreatic ß-cell central carbon metabolism and insulin secretion.

Pancreatic ß-cells respond to increased extracellular glucose levels by initiating a metabolic shift. That change in metabolism is part of the process of glucose-stimulated insulin secretion and is of particular interest in the context of diabetes. However, we do not fully understand how the coordinated changes in metabolic pathways and metabolite products influence insulin secretion. In this work, we apply systems biology approaches to develop a detailed kinetic model of the intracellular central carbon metabolic pathways in pancreatic ß-cells upon stimulation with high levels of glucose. The model is calibrated to published metabolomics datasets for the INS1 823/13 cell line, accurately capturing the measured metabolite fold-changes. We first employed the calibrated mechanistic model to estimate the stimulated cell's fluxome. We then used the predicted network fluxes in a data-driven approach to build a partial least squares regression model. By developing the combined kinetic and data-driven modeling framework, we gain insights into the link between ß-cell metabolism and glucose-stimulated insulin secretion. The combined modeling framework was used to predict the effects of common anti-diabetic pharmacological interventions on metabolite levels, flux through the metabolic network, and insulin secretion. Our simulations reveal targets that can be modulated to enhance insulin secretion. The model is a promising tool to contextualize and extend the usefulness of metabolomics data and to predict dynamics and metabolite levels that are difficult to measure in vitro. In addition, the modeling framework can be applied to identify, explain, and assess novel and clinically-relevant interventions that may be particularly valuable in diabetes treatment.

Computational analysis of 4-1BB-induced NFκB signaling suggests improvements to CAR cell design.

BACKGROUND: Chimeric antigen receptor (CAR)-expressing cells are a powerful modality of adoptive cell therapy against cancer. The potency of signaling events initiated upon antigen binding depends on the costimulatory domain within the structure of the CAR. One such costimulatory domain is 4-1BB, which affects cellular response via the NFκB pathway. However, the quantitative aspects of 4-1BB-induced NFκB signaling are not fully understood. METHODS: We developed an ordinary differential equation-based mathematical model representing canonical NFκB signaling activated by CD19scFv-4-1BB. After a global sensitivity analysis on model parameters, we ran Monte Carlo simulations of cell population-wide variability in NFκB signaling and quantified the mutual information between the extracellular signal and different levels of the NFκB signal transduction pathway. RESULTS: In response to a wide range of antigen concentrations, the magnitude of the transient peak in NFκB nuclear concentration varies significantly, while the timing of this peak is relatively consistent. Global sensitivity analysis showed that the model is robust to variations in parameters, and thus, its quantitative predictions would remain applicable to a broad range of parameter values. The model predicts that overexpressing NEMO and disabling IKKß deactivation can increase the mutual information between antigen levels and NFκB activation. CONCLUSIONS: Our modeling predictions provide actionable insights to guide CAR development. Particularly, we propose specific manipulations to the NFκB signal transduction pathway that can fine-tune the response of CD19scFv-4-1BB cells to the antigen concentrations they are likely to encounter. Video Abstract.

Modelling predicts differences in chimeric antigen receptor T-cell signalling due to biological variability.

In recent decades, chimeric antigen receptors (CARs) have been successfully used to generate engineered T cells capable of recognizing and eliminating cancer cells. The structure of CARs typically includes costimulatory domains, which enhance the T-cell response upon antigen encounter. However, it is not fully known how those co-stimulatory domains influence cell activation in the presence of biological variability. In this work, we used mathematical modelling to elucidate how the inclusion of one such costimulatory molecule, CD28, impacts the response of a population of CAR T cells under different sources of variability. Particularly, we demonstrate that CD28-bearing CARs mediate a faster and more consistent population response under both target antigen variability and kinetic rate variability. Next, we identify kinetic parameters that have the most impact on cell response time. Finally, based on our findings, we propose that enhancing the catalytic activity of lymphocyte-specific protein tyrosine kinase can result in drastically reduced and more consistent response times among heterogeneous CAR T-cell populations.

Phosphatases are predicted to govern prolactin-mediated JAK-STAT signaling in pancreatic beta cells.

Patients with diabetes are unable to produce a sufficient amount of insulin to properly regulate their blood glucose levels. One potential method of treating diabetes is to increase the number of insulin-secreting beta cells in the pancreas to enhance insulin secretion. It is known that during pregnancy, pancreatic beta cells proliferate in response to the pregnancy hormone, prolactin (PRL). Leveraging this proliferative response to PRL may be a strategy to restore endogenous insulin production for patients with diabetes. To investigate this potential treatment, we previously developed a computational model to represent the PRL-mediated JAK-STAT signaling pathway in pancreatic beta cells. Here, we applied the model to identify the importance of particular signaling proteins in shaping the response of a population of beta cells. We simulated a population of 10 000 heterogeneous cells with varying initial protein concentrations responding to PRL stimulation. We used partial least squares regression to analyze the significance and role of each of the varied protein concentrations in producing the response of the cell. Our regression models predict that the concentrations of the cytosolic and nuclear phosphatases strongly influence the response of the cell. The model also predicts that increasing PRL receptor strengthens negative feedback mediated by the inhibitor suppressor of cytokine signaling. These findings reveal biological targets that can potentially be used to modulate the proliferation of pancreatic beta cells to enhance insulin secretion and beta cell regeneration in the context of diabetes.

Exogenous Lactogenic Signaling Stimulates Beta Cell Replication In Vivo and In Vitro.

As patients recently diagnosed with T1D and patients with T2D have residual beta cell mass, there is considerable effort in beta cell biology to understand the mechanisms that drive beta cell regeneration as a potential cellular therapy for expanding patients' residual beta cell population. Both mouse and human studies have established that beta cell mass expansion occurs rapidly during pregnancy. To investigate the mechanisms of beta cell mass expansion during pregnancy, we developed a novel in vivo and in vitro models of pseudopregnancy. Our models demonstrate that pseudopregnancy promotes beta cell mass expansion in parous mice, and this expansion is driven by beta cell proliferation rather than hypertrophy. Importantly, estrogen, progesterone, and placental lactogen induce STAT5A signaling in the pseudopregnancy model, demonstrating that this model successfully recapitulates pregnancy-induced beta cell replication. We then created an in vitro model of pseudopregnancy and found that the combination of estrogen and placental lactogen induced beta cell replication in human islets and rat insulinoma cells. Therefore, beta cells both in vitro and in vivo increase proliferation when subjected to the pseudopregnancy cocktail compared to groups treated with estradiol or placental lactogen alone. The pseudopregnancy models described here may help inform novel methods of inducing beta cell replication in patients with diabetes.

Elucidating tumor-stromal metabolic crosstalk in colorectal cancer through integration of constraint-based models and LC-MS metabolomics.

Colorectal cancer (CRC) is a major cause of morbidity and mortality in the United States. Tumor-stromal metabolic crosstalk in the tumor microenvironment promotes CRC development and progression, but exactly how stromal cells, in particular cancer-associated fibroblasts (CAFs), affect the metabolism of tumor cells remains unknown. Here we take a data-driven approach to investigate the metabolic interactions between CRC cells and CAFs, integrating constraint-based modeling and metabolomic profiling. Using metabolomics data, we perform unsteady-state parsimonious flux balance analysis to infer flux distributions for central carbon metabolism in CRC cells treated with or without CAF-conditioned media. We find that CAFs reprogram CRC metabolism through stimulation of glycolysis, the oxidative arm of the pentose phosphate pathway (PPP), and glutaminolysis, as well as inhibition of the tricarboxylic acid cycle. To identify potential therapeutic targets, we simulate enzyme knockouts and find that CAF-treated CRC cells are especially sensitive to inhibitions of hexokinase and glucose-6-phosphate, the rate limiting steps of glycolysis and oxidative PPP. Our work gives mechanistic insights into the metabolic interactions between CRC cells and CAFs and provides a framework for testing hypotheses towards CRC-targeted therapies.