A data-driven computational model enables integrative and mechanistic characterization of dynamic macrophage polarization

A mechanistic systems biology model of brain microvascular endothelial cell signaling reveals dynamic pathway-based therapeutic targets for brain ischemia

Ischemic stroke is a significant threat to human health. Currently, there is a lack of effective treatments for stroke, and progress in new neuron-centered drug target development is relatively slow. On the other hand, studies have demonstrated that brain microvascular endothelial cells (BMECs) are crucial components of the neurovascular unit and play pivotal roles in ischemic stroke progression. To better understand the complex multifaceted roles of BMECs in the regulation of ischemic stroke pathophysiology and facilitate BMEC-based drug target discovery, we utilized a transcriptomics-informed systems biology modeling approach and constructed a mechanism-based computational multipathway model to systematically investigate BMEC function and its modulatory potential. Extensive multilevel data regarding complex BMEC pathway signal transduction and biomarker expression under various pathophysiological conditions were used for quantitative model calibration and validation, and we generated dynamic BMEC phenotype maps in response to various stroke-related stimuli to identify potential determinants of BMEC fate under stress conditions. Through high-throughput model sensitivity analyses and virtual target perturbations in model-based single cells, our model predicted that targeting succinate could effectively reverse the detrimental cell phenotype of BMECs under oxygen and glucose deprivation/reoxygenation, a condition that mimics stroke pathogenesis, and we experimentally validated the utility of this new target in terms of regulating inflammatory factor production, free radical generation and tight junction protection in vitro and in vivo. Our work is the first that complementarily couples transcriptomic analysis with mechanistic systems-level pathway modeling in the study of BMEC function and endothelium-based therapeutic targets in ischemic stroke.

Linking signal input, cell state, and spatial context to inflammatory responses

Signal integration is central to a causal understanding of appropriately scaled inflammatory responses. Here, we discuss recent progress in our understanding of the stimulus-response linkages downstream of pro-inflammatory inputs, with special attention to (1) the impact of cell state on the specificity of evoked gene expression and (2) the critical role of the spatial context of stimulus exposure. Advances in these directions are emerging from new tools for inferring cell-cell interactions and the activities of cytokines and transcription factors in complex microenvironments, enabling analysis of signal integration in tissue settings. Building on data-driven elucidation of factors driving inflammatory outcomes, mechanistic modeling can then contribute to a quantitative understanding of regulatory events that balance protective versus pathological inflammation.

Mathematical Modeling and Inference of Epidermal Growth Factor-Induced Mitogen-Activated Protein Kinase Cell Signaling Pathways

The mitogen-activated protein kinase (MAPK) pathway is an important intracellular signaling cascade that plays a key role in various cellular processes. Understanding the regulatory mechanisms of this pathway is essential for developing effective interventions and targeted therapies for related diseases. Recent advances in single-cell proteomic technologies have provided unprecedented opportunities to investigate the heterogeneity and noise within complex, multi-signaling networks across diverse cells and cell types. Mathematical modeling has become a powerful interdisciplinary tool that bridges mathematics and experimental biology, providing valuable insights into these intricate cellular processes. In addition, statistical methods have been developed to infer pathway topologies and estimate unknown parameters within dynamic models. This review presents a comprehensive analysis of how mathematical modeling of the MAPK pathway deepens our understanding of its regulatory mechanisms, enhances the prediction of system behavior, and informs experimental research, with a particular focus on recent advances in modeling and inference using single-cell proteomic data.

From transcriptomics to digital twins of organ function

Cell level functions underlie tissue and organ physiology. Gene expression patterns offer extensive views of the pathways and processes within and between cells. Single cell transcriptomics provides detailed information on gene expression within cells, cell types, subtypes and their relative proportions in organs. Functional pathways can be scalably connected to physiological functions at the cell and organ levels. Integrating experimentally obtained gene expression patterns with prior knowledge of pathway interactions enables identification of networks underlying whole cell functions such as growth, contractility, and secretion. These pathways can be computationally modeled using differential equations to simulate cell and organ physiological dynamics regulated by gene expression changes. Such computational systems can be thought of as parts of digital twins of organs. Digital twins, at the core, need computational models that represent in detail and simulate how dynamics of pathways and networks give rise to whole cell level physiological functions. Integration of transcriptomic responses and numerical simulations could simulate and predict whole cell functional outputs from transcriptomic data. We developed a computational pipeline that integrates gene expression timelines and systems of coupled differential equations to generate cell-type selective dynamical models. We tested our integrative algorithm on the eicosanoid biosynthesis network in macrophages. Converting transcriptomic changes to a dynamical model allowed us to predict dynamics of prostaglandin and thromboxane synthesis and secretion by macrophages that matched published lipidomics data obtained in the same experiments. Integration of cell-level system biology simulations with genomic and clinical data using a knowledge graph framework will allow us to create explicit predictive models that mechanistically link genomic determinants to organ function. Such integration requires a multi-domain ontological framework to connect genomic determinants to gene expression and cell pathways and functions to organ level phenotypes in healthy and diseased states. These integrated scalable models of tissues and organs as accurate digital twins predict health and disease states for precision medicine.

Endothelial cells signaling and patterning under hypoxia: a mechanistic integrative computational model including the Notch-Dll4 pathway

Introduction: Several signaling pathways are activated during hypoxia to promote angiogenesis, leading to endothelial cell patterning, interaction, and downstream signaling. Understanding the mechanistic signaling differences between endothelial cells under normoxia and hypoxia and their response to different stimuli can guide therapies to modulate angiogenesis. We present a novel mechanistic model of interacting endothelial cells, including the main pathways involved in angiogenesis. Methods: We calibrate and fit the model parameters based on well-established modeling techniques that include structural and practical parameter identifiability, uncertainty quantification, and global sensitivity. Results: Our results indicate that the main pathways involved in patterning tip and stalk endothelial cells under hypoxia differ, and the time under hypoxia interferes with how different stimuli affect patterning. Additionally, our simulations indicate that Notch signaling might regulate vascular permeability and establish different Nitric Oxide release patterns for tip/stalk cells. Following simulations with various stimuli, our model suggests that factors such as time under hypoxia and oxygen availability must be considered for EC pattern control. Discussion: This project provides insights into the signaling and patterning of endothelial cells under various oxygen levels and stimulation by VEGFA and is our first integrative approach toward achieving EC control as a method for improving angiogenesis. Overall, our model provides a computational framework that can be built on to test angiogenesis-related therapies by modulation of different pathways, such as the Notch pathway.

Getting everyone to agree on gene signatures for murine macrophage polarization in vitro

Macrophages, key players in the innate immune system, showcase remarkable adaptability. Derived from monocytes, these phagocytic cells excel in engulfing and digesting pathogens and foreign substances as well as contributing to antigen presentation, initiating and regulating adaptive immunity. Macrophages are highly plastic, and the microenvironment can shaper their phenotype leading to numerous distinct polarized subsets, exemplified by the two ends of the spectrum: M1 (classical activation, inflammatory) and M2 (alternative activation, anti-inflammatory). RNA sequencing (RNA-Seq) has revolutionized molecular biology, offering a comprehensive view of transcriptomes. Unlike microarrays, RNA-Seq detects known and novel transcripts, alternative splicing, and rare transcripts, providing a deeper understanding of genome complexity. Despite the decreasing costs of RNA-Seq, data consolidation remains limited, hindering noise reduction and the identification of authentic signatures. Macrophages polarization is routinely ascertained by qPCR to evaluate those genes known to be characteristic of M1 or M2 skewing. Yet, the choice of these genes is literature- and experience-based, lacking therefore a systematic approach. This manuscript builds on the significant increase in deposited RNA-Seq datasets to determine an unbiased and robust murine M1 and M2 polarization profile. We now provide a consolidated list of global M1 differentially expressed genes (i.e. robustly modulated by IFN-γ, LPS, and LPS+ IFN-γ) as well as consolidated lists of genes modulated by each stimulus (IFN-γ, LPS, LPS+ IFN-γ, and IL-4).

Signaling-biophysical modeling unravels mechanistic control of red blood cell phagocytosis by macrophages in sickle cell disease

Red blood cell (RBC) aging manifests through progressive changes in cell morphology, rigidity, and expression of membrane proteins. To maintain the quality of circulating blood, splenic macrophages detect the biochemical signals and biophysical changes of RBCs and selectively clear them through erythrophagocytosis. In sickle cell disease (SCD), RBCs display alterations affecting their interaction with macrophages, leading to aberrant phagocytosis that may cause life-threatening spleen sequestration crises. To illuminate the mechanistic control of RBC engulfment by macrophages in SCD, we integrate a system biology model of RBC-macrophage signaling interactions with a biophysical model of macrophage engulfment, as well as in vitro phagocytosis experiments using the spleen-on-a-chip technology. Our modeling framework accurately predicts the phagocytosis dynamics of RBCs under different disease conditions, reveals patterns distinguishing normal and sickle RBCs, and identifies molecular targets including Src homology 2 domain-containing protein tyrosine phosphatase-1 (SHP1) and cluster of differentiation 47 (CD47)/signal regulatory protein α (SIRPα) as therapeutic targets to facilitate the controlled clearance of sickle RBCs in the spleen.

A comparative analysis of 2D and 3D experimental data for the identification of the parameters of computational models

Computational models are becoming an increasingly valuable tool in biomedical research. Their accuracy and effectiveness, however, rely on the identification of suitable parameters and on appropriate validation of the in-silico framework. Both these steps are highly dependent on the experimental model used as a reference to acquire the data. Selecting the most appropriate experimental framework thus becomes key, together with the analysis of the effect of combining results from different experimental models, a common practice often necessary due to limited data availability. In this work, the same in-silico model of ovarian cancer cell growth and metastasis, was calibrated with datasets acquired from traditional 2D monolayers, 3D cell culture models or a combination of the two. The comparison between the parameters sets obtained in the different conditions, together with the corresponding simulated behaviours, is presented. It provides a framework for the study of the effect of the different experimental models on the development of computational systems. This work also provides a set of general guidelines for the comparative testing and selection of experimental models and protocols to be used for parameter optimization in computational models.

Experiment-based computational model predicts that IL-6 classic and trans-signaling exhibit similar potency in inducing downstream signaling in endothelial cells

Inflammatory cytokine mediated responses are important in the development of many diseases that are associated with angiogenesis. Targeting angiogenesis as a prominent strategy has shown limited effects in many contexts such as cardiovascular diseases and cancer. One potential reason for the unsuccessful outcome is the mutual dependent role between inflammation and angiogenesis. Inflammation-based therapies primarily target inflammatory cytokines such as interleukin-6 (IL-6) in T cells, macrophages, cancer cells, and muscle cells, and there is a limited understanding of how these cytokines act on endothelial cells. Thus, we focus on one of the major inflammatory cytokines, IL-6, mediated intracellular signaling in endothelial cells by developing a detailed computational model. Our model quantitatively characterized the effects of IL-6 classic and trans-signaling in activating the signal transducer and activator of transcription 3 (STAT3), phosphatidylinositol 3-kinase/protein kinase B (PI3K/Akt), and mitogen-activated protein kinase (MAPK) signaling to phosphorylate STAT3, extracellular regulated kinase (ERK) and Akt, respectively. We applied the trained and validated experiment-based computational model to characterize the dynamics of phosphorylated STAT3 (pSTAT3), Akt (pAkt), and ERK (pERK) in response to IL-6 classic and/or trans-signaling. The model predicts that IL-6 classic and trans-signaling induced responses are IL-6 and soluble IL-6 receptor (sIL-6R) dose-dependent. Also, IL-6 classic and trans-signaling showed similar potency in inducing downstream signaling; however, trans-signaling induces stronger downstream responses and plays a dominant role in the overall effects from IL-6 due to the in vitro experimental setting of abundant sIL-6R. In addition, both IL-6 and sIL-6R levels regulate signaling strength. Moreover, our model identifies the influential species and kinetic parameters that specifically modulate the downstream inflammatory and/or angiogenic signals, pSTAT3, pAkt, and pERK responses. Overall, the model predicts the effects of IL-6 classic and/or trans-signaling stimulation quantitatively and provides a framework for analyzing and integrating experimental data. More broadly, this model can be utilized to identify potential targets that influence IL-6 mediated signaling in endothelial cells and to study their effects quantitatively in modulating STAT3, Akt, and ERK activation.

Polarization of Melatonin-Modulated Colostrum Macrophages in the Presence of Breast Tumor Cell Lines

Human colostrum and milk contain diverse cells and soluble components that have the potential to act against tumors. In breast cancer, macrophages play a significant role in immune infiltration and contribute to the progression and spread of tumors. However, studies suggest that these cells can be reprogrammed to act as an antitumor immune response. This study aimed to evaluate the levels of melatonin and its receptors, MT1 (melatonin receptor 1) and MT2 (melatonin receptor 2), in colostrum and assess the differentiation and polarization of the colostrum macrophages modulated by melatonin in the presence of breast tumor cells. Colostrum samples were collected from 116 mothers and tested for their melatonin and receptor levels. The colostrum cells were treated with or without melatonin and then cultured for 24 h in the presence or absence of breast tumor cells. The results showed that melatonin treatment increased the expression of MT1 and MT2 in the colostrum cells. Furthermore, melatonin treatment increased the percentage of M1 macrophages and decreased the percentage of M2 macrophages. When the colostrum macrophages were cocultured with breast tumor cells, melatonin reduced the percentage of both macrophage phenotypes and the cytokines tumor necrosis factor-alpha (TNF-α) and interleukin 8 (IL-8). These data suggest that melatonin can regulate the inflammatory process via M1 macrophages in the tumor microenvironment and, simultaneously, the progression of M2 macrophages that favor tumorigenesis.

Signaling is the pathway to macrophage function

Tissue and inflammatory contexts are well appreciated to shape macrophage function to promote health or disease. However, there has been minimal progress towards understanding how these contexts modify signaling-to-transcription networks. Integration of mechanistic modeling and data-driven approaches will be crucial for investigating how cell state impacts macrophage decision-making.

A tractable physical model for the yeast polarity predicts epistasis and fitness

Accurate phenotype prediction based on genetic information has numerous societal applications, such as crop design or cellular factories. Epistasis, when biological components interact, complicates modelling phenotypes from genotypes. Here we show an approach to mitigate this complication for polarity establishment in budding yeast, where mechanistic information is abundant. We coarse-grain molecular interactions into a so-called mesotype, which we combine with gene expression noise into a physical cell cycle model. First, we show with computer simulations that the mesotype allows validation of the most current biochemical polarity models by quantitatively matching doubling times. Second, the mesotype elucidates epistasis emergence as exemplified by evaluating the predicted mutational effect of key polarity protein Bem1p when combined with known interactors or under different growth conditions. This example also illustrates how unlikely evolutionary trajectories can become more accessible. The tractability of our biophysically justifiable approach inspires a road-map towards bottom-up modelling complementary to statistical inferences.

This article is part of the theme issue ‘Interdisciplinary approaches to predicting evolutionary biology’.

Experiment-based Computational Model Predicts that IL-6 Trans-Signaling Plays a Dominant Role in IL-6 mediated signaling in Endothelial Cells

Inflammatory cytokine mediated responses are important in the development of many diseases that are associated with angiogenesis. Targeting angiogenesis as a prominent strategy has shown limited effects in many contexts such as peripheral arterial disease (PAD) and cancer. One potential reason for the unsuccessful outcome is the mutual dependent role between inflammation and angiogenesis. Inflammation-based therapies primarily target inflammatory cytokines such as interleukin-6 (IL-6) in T cells, macrophages, cancer cells, muscle cells, and there is a limited understanding of how these cytokines act on endothelial cells. Thus, we focus on one of the major inflammatory cytokines, IL-6, mediated intracellular signaling in endothelial cells by developing a detailed computational model. Our model quantitatively characterized the effects of IL-6 classic and trans-signaling in activating the signal transducer and activator of transcription 3 (STAT3), phosphatidylinositol 3-kinase/protein kinase B (PI3K/Akt), and mitogen-activated protein kinase (MAPK) signaling to phosphorylate STAT3, extracellular regulated kinase (ERK) and Akt, respectively. We applied the trained and validated experiment-based computational model to characterize the dynamics of phosphorylated STAT3 (pSTAT3), Akt (pAkt), and extracellular regulated kinase (pERK) in response to IL-6 classic and/or trans-signaling. The model predicts that IL-6 classic and trans-signaling induced responses are IL-6 and soluble IL-6 receptor (sIL-6R) dose-dependent. Also, IL-6 trans-signaling induces stronger downstream signaling and plays a dominant role in the overall effects from IL-6. In addition, both IL-6 and sIL-6R levels regulate signaling strength. Moreover, our model identifies the influential species and kinetic parameters that specifically modulate the pSTAT3, pAkt, and pERK responses, which represent potential targets for inflammatory cytokine mediated signaling and angiogenesis-based therapies. Overall, the model predicts the effects of IL-6 classic and/or trans-signaling stimulation quantitatively and provides a framework for analyzing and integrating experimental data. More broadly, this model can be utilized to identify targets that influence inflammatory cytokine mediated signaling in endothelial cells and to study the effects of angiogenesis- and inflammation-based therapies.

A novel strategy to dissect multifaceted macrophage function in human diseases

Macrophages are widely distributed immune cells that play central roles in a variety of physiologic and pathologic processes, including obesity and cardiovascular disease (CVD). They are highly plastic cells that execute diverse functions according to a combination of signaling and environmental cues. While macrophages have traditionally been understood to polarize to either proinflammatory M1-like or anti-inflammatory M2-like states, evidence has shown that they exist in a spectrum of states between those 2 phenotypic extremes. In obesity-related disease, M1-like macrophages exacerbate inflammation and promote insulin resistance, while M2-like macrophages reduce inflammation, promoting insulin sensitivity. However, polarization markers are expressed inconsistently in adipose tissue macrophages, and they additionally exhibit phenotypes differing from the M1/M2 paradigm. In atherosclerotic CVD, activated plaque macrophages can also exist in a range of proinflammatory or anti-inflammatory states. Some of these macrophages scavenge lipids, developing into heterogeneous foam cell populations. To better characterize the many actions of macrophages in human disease, we have designed a novel set of computational tools: MacSpectrum and AtheroSpectrum. These tools provide information on the inflammatory polarization status, differentiation, and foaming of macrophages in both human and mouse samples, allowing for better characterization of macrophage subpopulations based on their function. Using these tools, we identified disease-relevant cell states in obesity and CVD, including the novel concept that macrophage-derived foam cell formation can follow homeostatic noninflammatory or pathogenic inflammatory foaming programs.

The osteogenetic activities of mesenchymal stem cells in response to Mg2+ ions and inflammatory cytokines: a numerical approach using fuzzy logic controllers

Magnesium (Mg²⁺) ions are frequently reported to regulate osteogenic activities of mesenchymal stem cells (MSCs). In this study, we propose a numerical model to study the regulatory importance of Mg²⁺ ions on MSCs osteoblastic differentiation in the presence of an inflammatory response. A fuzzy logic controller was formulated to receive the concentrations of Mg²⁺ ions and the inflammatory cytokines of TNF-α, IL-10, IL-1β, and IL-8 as cellular inputs and predict the cells’ early and late differentiation rates. Five sets of empirical data obtained from published cell culture experiments were used to calibrate the model. The model successfully reproduced the empirical data regarding the concentration- and phase-dependent effect of Mg²⁺ ions on the differentiation process. In agreement with the experiments, the model showed the stimulatory role of Mg²⁺ ions on the early differentiation phase, once administered at low concentration, and their inhibitory role on the late differentiation phase. The numerical approach used in this study suggested 6–8 mM as the most effective concentration of Mg²⁺ ions in promoting the early differentiation process. Also, the proposed model sheds light on the fundamental differences in the behavioral properties of cells cultured in different experiments, e.g. differentiation rate and the sensitivity of the cultured cells to stimulatory signals such as Mg²⁺ ions. Thus, it can be used to interpret and compare different empirical findings. Moreover, the model successfully reproduced the nonlinearities in the concentration-dependent role of the inflammatory cytokines in early and late differentiation rates. Overall, the proposed model can be employed in studying the osteogenic properties of Mg-based implants in the presence of an inflammatory response.

Magnesium (Mg) is an attractive material for bone implants as it fully degrades after implantation, saving pain and cost of the second surgery for implant removal. To advance its application in the orthopedic industry, it is paramount to fully understand the biological impact of the degradation products, in particular Mg²⁺ ions. Here, we propose a computer model to study the effects of Mg²⁺ ions on bone regeneration. The model focuses on stem cells and includes both the direct stimulation effects of Mg²⁺ ions on cells and the indirect stimulus through the inflammatory system. The proposed model successfully reproduced the experimental data of five different studies. The model additionally highlighted differences amongst different experiments in terms of the cellular response to Mg²⁺ ions. The proposed system therefore provides an important addition to the field of Mg implant research.

Dynamics of tumor-associated macrophages in a quantitative systems pharmacology model of immunotherapy in triple-negative breast cancer

Quantitative systems pharmacology (QSP) modeling is an emerging mechanistic computational approach that couples drug pharmacokinetics/pharmacodynamics and the course of disease progression. It has begun to play important roles in drug development for complex diseases such as cancer, including triple-negative breast cancer (TNBC). The combination of the anti-PD-L1 antibody atezolizumab and nab-paclitaxel has shown clinical activity in advanced TNBC with PD-L1-positive tumor-infiltrating immune cells. As tumor-associated macrophages (TAMs) serve as major contributors to the immuno-suppressive tumor microenvironment, we incorporated the dynamics of TAMs into our previously published QSP model to investigate their impact on cancer treatment. We show that through proper calibration, the model captures the macrophage heterogeneity in the tumor microenvironment while maintaining its predictive power of the trial results at the population level. Despite its high mechanistic complexity, the modularized QSP platform can be readily reproduced, expanded for new species of interest, and applied in clinical trial simulation.

•

A mechanistic model of quantitative systems pharmacology in immuno-oncology

•

Dynamics of tumor-associated macrophages are integrated into our previous work

•

Conducting in silico clinical trials to predict clinical response to cancer therapy

Artificial neural networks enable genome-scale simulations of intracellular signaling

Mammalian cells adapt their functional state in response to external signals in form of ligands that bind receptors on the cell-surface. Mechanistically, this involves signal-processing through a complex network of molecular interactions that govern transcription factor activity patterns. Computer simulations of the information flow through this network could help predict cellular responses in health and disease. Here we develop a recurrent neural network framework constrained by prior knowledge of the signaling network with ligand-concentrations as input and transcription factor-activity as output. Applied to synthetic data, it predicts unseen test-data (Pearson correlation r = 0.98) and the effects of gene knockouts (r = 0.8). We stimulate macrophages with 59 different ligands, with and without the addition of lipopolysaccharide, and collect transcriptomics data. The framework predicts this data under cross-validation (r = 0.8) and knockout simulations suggest a role for RIPK1 in modulating the lipopolysaccharide response. This work demonstrates the feasibility of genome-scale simulations of intracellular signaling.

Restoring Pro-healing/remodeling- associated M2a/c Macrophages using ON101 Accelerates Diabetic Wound Healing

Diabetic wounds exhibit chronic inflammation and delayed tissue proliferation or remodeling, mainly owing to prolonged proinflammatory (M1) macrophage activity and defects in transition to prohealing/proremodeling (M2a/M2c; CD206⁺ and/or CD163⁺) macrophages. We found that topical treatment with ON101, a plant-based potential therapeutic for diabetic foot ulcers, increased M2c-like (CD163⁺ and CD206⁺) cells and suppressed M1-like cells, altering the inflammatory gene profile in a diabetic mouse model compared with that in the controls. An in vitro macrophage-polarizing model revealed that ON101 directly suppressed CD80⁺ and CD86⁺ M1-macrophage polarization and M1-associated proinflammatory cytokines at both protein and transcriptional levels. Notably, conditioned medium collected from ON101-treated M1 macrophages reversed the M1-conditioned medium‒mediated suppression of CD206⁺ macrophages. Furthermore, conditioned medium from ON101-treated adipocyte progenitor cells significantly promoted CD206⁺ and CD163⁺ macrophages but strongly inhibited M1-like cells. ON101 treatment also stimulated the expression of GCSF and CXCL3 genes in human adipocyte progenitor cells. Interestingly, treatment with recombinant GCSF protein enhanced both CD206⁺ and CD163⁺ M2 markers, whereas CXCL3 treatment only stimulated CD163⁺ M2 macrophages. Depletion of cutaneous M2 macrophages inhibited ON101-induced diabetic wound healing. Thus, ON101 directly suppressed M1 macrophages and facilitated the GCSF- and CXCL3-mediated transition from M1 to M2 macrophages, lowering inflammation and leading to faster diabetic wound healing.

Dynamic Multiscale Regulation of Perfusion Recovery in Experimental Peripheral Arterial Disease: A Mechanistic Computational Model

In peripheral arterial disease (PAD), the degree of endogenous capacity to modulate revascularization of limb muscle is central to the management of leg ischemia. To characterize the multiscale and multicellular nature of revascularization in PAD, we have developed the first computational systems biology model that mechanistically incorporates intracellular, cellular, and tissue-level features critical for the dynamic reconstitution of perfusion after occlusion-induced ischemia. The computational model was specifically formulated for a preclinical animal model of PAD (mouse hindlimb ischemia [HLI]), and it has gone through multilevel model calibration and validation against a comprehensive set of experimental data so that it accurately captures the complex cellular signaling, cell-cell communication, and function during post-HLI perfusion recovery. As an example, our model simulations generated a highly detailed description of the time-dependent spectrum-like macrophage phenotypes in HLI, and through model sensitivity analysis we identified key cellular processes with potential therapeutic significance in the pathophysiology of PAD. Furthermore, we computationally evaluated the in vivo effects of different targeted interventions on post-HLI tissue perfusion recovery in a model-based, data-driven, virtual mouse population and experimentally confirmed the therapeutic effect of a novel model-predicted intervention in real HLI mice. This novel multiscale model opens up a new avenue to use integrative systems biology modeling to facilitate translational research in PAD.

Systems biology of angiogenesis signaling: Computational models and omics

Angiogenesis is a highly regulated multiscale process that involves a plethora of cells, their cellular signal transduction, activation, proliferation, differentiation, as well as their intercellular communication. The coordinated execution and integration of such complex signaling programs is critical for physiological angiogenesis to take place in normal growth, development, exercise, and wound healing, while its dysregulation is critically linked to many major human diseases such as cancer, cardiovascular diseases, and ocular disorders; it is also crucial in regenerative medicine. Although huge efforts have been devoted to drug development for these diseases by investigation of angiogenesis‐targeted therapies, only a few therapeutics and targets have proved effective in humans due to the innate multiscale complexity and nonlinearity in the process of angiogenic signaling. As a promising approach that can help better address this challenge, systems biology modeling allows the integration of knowledge across studies and scales and provides a powerful means to mechanistically elucidate and connect the individual molecular and cellular signaling components that function in concert to regulate angiogenesis. In this review, we summarize and discuss how systems biology modeling studies, at the pathway‐, cell‐, tissue‐, and whole body‐levels, have advanced our understanding of signaling in angiogenesis and thereby delivered new translational insights for human diseases.

This article is categorized under:

Cardiovascular Diseases > Computational Models

Cancer > Computational Models

MECHANISTIC AND DATA-DRIVEN MODELS OF CELL SIGNALING: TOOLS FOR FUNDAMENTAL DISCOVERY AND RATIONAL DESIGN OF THERAPY

A full understanding of cell signaling processes requires knowledge of protein structure/function relationships, protein-protein interactions, and the abilities of pathways to control phenotypes. Computational models offer a valuable framework for integrating that knowledge to predict the effects of system perturbations and interventions in health and disease. Whereas mechanistic models are well suited for understanding the biophysical basis for signal transduction and principles of therapeutic design, data-driven models are particularly suited to distill complex signaling relationships among samples and between multivariate signaling changes and phenotypes. Both approaches have limitations and provide incomplete representations of signaling biology, but their careful implementation and integration can provide new understanding for how manipulating system variables impacts cellular decisions.

Integrating single cell sequencing with a spatial quantitative systems pharmacology model spQSP for personalized prediction of triple-negative breast cancer immunotherapy response

Response to cancer immunotherapies depends on the complex and dynamic interactions between T cell recognition and killing of cancer cells that are counteracted through immunosuppressive pathways in the tumor microenvironment. Therefore, while measurements such as tumor mutational burden provide biomarkers to select patients for immunotherapy, they neither universally predict patient response nor implicate the mechanisms that underlie immunotherapy resistance. Recent advances in single-cell RNA sequencing technology measure cellular heterogeneity within cells of an individual tumor but have yet to realize the promise of predictive oncology. In addition to data, mechanistic multiscale computational models are developed to predict treatment response. Incorporating single-cell data from tumors to parameterize these computational models provides deeper insights into prediction of clinical outcome in individual patients. Here, we integrate whole-exome sequencing and scRNA-seq data from Triple-Negative Breast Cancer patients to model neoantigen burden in tumor cells as input to a spatial Quantitative System Pharmacology model. The model comprises a four-compartmental Quantitative System Pharmacology sub-model to represent a whole patient and a spatial agent-based sub-model to represent tumor volumes at the cellular scale. We use the high-throughput single-cell data to model the role of antigen burden and heterogeneity relative to the tumor microenvironment composition on predicted immunotherapy response. We demonstrate how this integrated modeling and single-cell analysis framework can be used to relate neoantigen heterogeneity to immunotherapy treatment outcomes. Our results demonstrate feasibility of merging single-cell data to initialize cell states in multiscale computational models such as the spQSP for personalized prediction of clinical outcomes to immunotherapy.

Protocol for simulating macrophage signal transduction and phenotype polarization using a large-scale mechanistic computational model

The ability to measure and analyze the complex dynamic multi-marker features of macrophages is critical for the understanding of their diverse phenotypes and functions in health and disease. To that end, we have recently developed a multi-pathway computational model that for the first time enables a systems-level characterization of macrophage signaling and activation from quantitative, temporal, dose-dependent, and single-cell aspects. This protocol includes instructions to utilize this model to computationally explore different biological scenarios with high resolution and efficiency. For complete details on the use and execution of this protocol, please refer to Zhao et al. (2021).