Sara Sadat Aghamiri
Bhanwar Lal Puniya
Rada Amin
Tomáš Helikar- Department of Biochemistry, University of Nebraska-Lincoln, Lincoln, NE, United States
Dendritic cells (DCs) are professional antigen-presenting cells (APCs) with the unique ability to mediate inflammatory responses of the immune system. Given the critical role of DCs in shaping immunity, they present an attractive avenue as a therapeutic target to program the immune system and reverse immune disease disorders. To ensure appropriate immune response, DCs utilize intricate and complex molecular and cellular interactions that converge into a seamless phenotype. Computational models open novel frontiers in research by integrating large-scale interaction to interrogate the influence of complex biological behavior across scales. The ability to model large biological networks will likely pave the way to understanding any complex system in more approachable ways. We developed a logical and predictive model of DC function that integrates the heterogeneity of DCs population, APC function, and cell-cell interaction, spanning molecular to population levels. Our logical model consists of 281 components that connect environmental stimuli with various layers of the cell compartments, including the plasma membrane, cytoplasm, and nucleus to represent the dynamic processes within and outside the DC, such as signaling pathways and cell-cell interactions. We also provided three sample use cases to apply the model in the context of studying cell dynamics and disease environments. First, we characterized the DC response to Sars-CoV-2 and influenza co-infection by in-silico experiments and analyzed the activity level of 107 molecules that play a role in this co-infection. The second example presents simulations to predict the crosstalk between DCs and T cells in a cancer microenvironment. Finally, for the third example, we used the Kyoto Encyclopedia of Genes and Genomes enrichment analysis against the model’s components to identify 45 diseases and 24 molecular pathways that the DC model can address. This study presents a resource to decode the complex dynamics underlying DC-derived APC communication and provides a platform for researchers to perform in-silico experiments on human DC for vaccine design, drug discovery, and immunotherapies.
Highlights
● The predictive model of the human Dendritic Cell (DC) bridges the gap between experimental data and the in-silico simulations.
● Constructing the first large-scale logical model of Dendritic Cell
● Applications of the DC model in human immunology
Introduction
Dendritic cells (DCs) comprise a diverse set of antigen-presenting cells that are responsible for the recognition of foreign and self-antigens and the subsequent regulation and initiation of specialized adaptive and innate immune responses (1, 2). Via pattern recognition receptors, DCs can sense a wide range of epitopes expressed by pathogens and damaged cells (3). The sophisticated ontogeny of DCs enables them to maintain tolerance in the presence of foreign and self-antigens or to initiate an inflammatory response (4–6). Striking the right balance to antigen response puts DCs in a critical pathway for disease management (7, 8). An insufficient immune response to an antigen can suppress downstream cell differentiation leading to an increased risk of infection and malignancy (9, 10). An over-reactive or chronic immune response, however, can lead to auto-immune diseases, allergies, and chronic inflammation (11, 12).
DCs mediate adaptive responses through cell-cell interactions (e.g., antigen presentation via the major histocompatibility complex (MHC) classes), the increase of co-stimulatory immune checkpoint ligands/receptors, and through the secretion of pro-and anti-inflammatory interleukins, growth factors, and chemokines (13, 14). Antigen recognition triggers a cascade of signaling pathways that switch the DC cellular state from tolerant (immature) to inflammatory (mature) (15, 16). DCs comprise three major subtypes with distinct immunogenicity and plasticity: conventional DCs (cDC1 and cDC2), plasmacytoid DCs (pDCs), and monocyte-derived DCs (MoDCs). As mature DCs, they can prime effectors and suppressors cells (e.g., lymphocyte T and B cells) to stimulate a wide range of immune responses (17, 18).
The significance of DCs in identifying and initiating an adaptive response to foreign and self-antigens has stimulated interest in isolating DCs as a potential therapeutic tool to program specific immune responses to pathogens and malignant cells (19–21). For example, in 2010, a DC-based vaccine was approved for the prevention of prostate cancer (22). DC-based vaccine development for other diseases has not been as successful; achieving full maturation of DCs and a limited ability for DCs to activate T cells are some of the challenges that have been encountered (23). Improved methods for characterizing and perturbing the complex mechanisms underlying DC maturation in the context of the broader immune system may help translate biological knowledge to clinical applications. Computational models, for example, have been gaining traction as a means to study the dynamics of immune responses in the context of homeostasis and diseased states (24–26) by utilizing a variety of mathematical frameworks (27–29) to represent multiple levels of biological regulation (e.g., genome-scale metabolic network regulation, signal transduction, cell-to-cell communication, etc.). Despite previous modeling efforts of immune-related biological diseases, a large-scale model of major DC functions and its communication with other immune cells is still lacking.
Multiscale models have the potential to uncover the underlying mechanisms behind emergent behaviors at various scales such as intracellular, cellular, and systemic levels. Such models can consider various temporal and organizational scales including signal transduction, gene regulation, metabolism, cellular behaviors, and cytokine transport (25). Different multiscale models have been developed to study the dynamic response of DC under different extracellular environments. For example, Klinke II DJ. developed a multiscale model to investigate the impact of the lung microenvironment on the education of DC for optimal T cell polarization. The model considered DC trafficking and education in the lung while taking into consideration the time, maturation, spatial distribution and IL12 response (30). Mei Y. and colleagues created a multiscale platform, the ENteric Immune Simulator (ENISI), to study the mucosal immune response during colonic inflammation. The multiscale tool has the advantage of connecting three different scales - intracellular, cellular, and tissular - using different mathematical languages (31). Lai et al. developed a multiscale model of DC-based vaccine by considering the signaling pathways underlying DC maturation, the bio-distribution of DCs in multiple organs, and the DC-T-cell response to identify optimal targets for enhancing anti-cancer DC vaccination in the context of melanoma (28). Although those multiscale models included a tissue scale to study the dynamic distribution of DC in different organs, so far, these models have only focused on DC functionality under specific disease conditions or specific signaling pathways. However, DC functions involve complex intracellular and cellular networks that are critical for regulating cell activation and initiating immune responses.
Logical modeling formalism has emerged as a particularly effective approach to modeling large-scale biological systems due to scalability and independence of kinetic parameters that are largely unknown (32–35). Logical models of different scales and complexity (a few to hundreds of components) have been applied to study various biological and translational questions (36), such as studying cellular crosstalk (37) and predicting cellular phenotypes (24, 38) and drug targets (39).
Here, we present a multiscale mechanistic model of human DCs that captures the complex interplay of intracellular molecular signaling to intercellular cell-cell communications. DC model enables researchers to easily modify, expand and test new hypotheses of the immune system. Our aim is to provide the researchers with computational tools to gain insight into DC mechanisms and disease pathology. The mechanistic model uses the logical mathematical framework (40) and focuses on signal transduction networks responsible for regulating DCs’ antigen-presenting cellular function, cellular interactions, maturation process, and immune cell population dynamics. It captures the dynamic biological events in response to diverse stimuli (pathogens, malignancy) and the downstream biological coordination between surface molecules (receptors, integrins, lectins), signal transduction (kinases, enzymes, transcription factors), and secretory factors (cytokines, chemokines). Two diseases are highlighted to demonstrate the utility of the model under diverse conditions. Lastly, receptor-ligand interactions between DCs and four immune cell types that DCs commonly interact with (T cells, B cells, natural killer (NK) cells, and neutrophils) have also been represented. The results of in-silico simulations of the model under various environmental conditions and network perturbations were validated using peer-reviewed published literature.
Methods
Model construction
The computational model is a mechanistic, logic-based model. Each component of the model can assume an active (1) or inactive (0) state at any time t. The activity state of the model’s internal components is determined by the regulatory mechanisms of other directly interacting components. These regulatory mechanisms are described with Boolean functions comprised of AND, OR, NOT operators (40).
To gain a comprehensive understanding of the molecular pathways involved in dendritic cells and antigen-presenting cells, we conducted a systematic search of the literature using PubMed. Our search was specifically focused on exploring the molecular pathways involved in each DC subtype: pDC, cDC1, cDC2, and MoDC. To limit the search results, we utilized a combination of search terms, including: “dendritic cells AND antigen-presenting cell AND MoDC AND molecular pathway,” “dendritic cells AND antigen-presenting cell AND pDC AND molecular pathway,” “dendritic cells AND antigen-presenting cell AND cDC1 AND molecular pathway,” and “dendritic cells AND antigen-presenting cell AND cDC2 AND molecular pathway.” This comprehensive search allowed us to obtain a wealth of information related to the molecular pathways involved in each DC subtype, providing a foundation for our investigation into the function and activation of these important immune cells. In the development of our model, we followed strict data inclusion criteria, limiting our selection to original research articles focused on healthy human subjects. Studies using mice and clinical trials were excluded from our manual literature mining process. The first draft of the model was constructed using the information obtained from the manual literature mining of the original studies. Upon reviewing the initial draft, we conducted a supplementary search of the literature utilizing both review and original studies to obtain well-established biological information related to regulators of the unconnected components. This thorough and systematic approach allowed us to develop a comprehensive model using 92 publications (83 original and 9 reviews) that represents the molecular pathways involved in dendritic cells. We defined subtype-specific markers to differentiate between pDC, cDC1, cDC2, MoDC (Results and Supplemental Figure S1A).
To validate the model, we collected literature reporting specific DC response to different extracellular conditions and simulated emergent behaviors that were not directly programmed into the model (41). Because logical models are of qualitative nature, model validations focus on the ability of the model to reproduce qualitative behaviors seen in wet-lab experiments (e.g., change in activity level of a component(s) under specific extracellular conditions) - a standard process for logical models (24, 41–43).
From the publications, we retrieved information related to DC-specific stimuli, the effect of the studied environment, and comprehensive signaling pathways (receptors, kinases, transcription factors). The model consists of 281 components. These components are categorized into various classifications and compartments. There are 178 proteins, 87 RNAs, and 16 components representing phenotypes and cells. The proteins are organized in cell membrane (64 components), cytoplasm (40 components), nucleus (22 components), extracellular space (52 components). Figure 1 Created with BioRender.com.
To standardize the naming convention of the components in the model, we used protein and gene names from the HUGO Gene Nomenclature Committee (HGNC) (44). The model was built in the web-based modeling and analysis platform, Cell Collective, and manually curated using the aforementioned literature (45). All components used to build the regulatory mechanisms have been annotated in Cell Collective with the exact quote from the reference literature. The model is publicly available in Cell Collective (under Published Models) where it can be simulated as well as downloaded (and other logical models published by the community) in several file formats (such as SBML-qual, text file of logical functions, and truth tables) (46, 47).
Model simulations and analyses
Cell Collective was used to perform all computational simulations and analyses of the model. Cell Collective uses discrete mathematics to construct the model, but the simulated output values are semi-continuous, ranging from 0 to 100% activity levels (48, 49). The activity levels of external components are unitless and defined as a percent chance (probability * 100) of the component being active in a specific time t (24). Depending on the desired experiment, the activity levels of external components can be set by the user to specific values, or they can be set to ranges from which values during each simulation are selected randomly (e.g., to simulate dose-response experiments).
We used Cell Collective for two types of analyses: real-time and dose-response using asynchronous updates such that all genes take different times to make a transition, which is closer to biological phenomena (50). The initial condition of the model was set to immature cellular phenotype as 1 (active) and all other components were set to zero since DCs are considered immature under the resting condition and before stimuli activation (51). The immature DCs are recruited to the inflamed site by pathogen signals, capture foreign antigens and undergo maturation to DC subtypes (52). While simulating the model in Cell Collective, the user can define the activity levels of external components to a specific point or provide ranges (e.g., varying between 0% to 100%). When a range is defined for external components, their activity levels are selected randomly in each simulation. In the real-time simulation, we showed the activity of components at different times (steps), which was presented using the mean activity level of multiple simulations. For dose-response analysis, we conducted each simulation consisting of 800 steps. The activity levels output components are fractions of ones over the last 300 iterations (500 to 800 steps) that describe the model’s steady behavior as described by (48, 49). Under each environment set for a biological scenario, we used 1,000 simulations.
Global sensitivity analysis
We used sensitivity analysis in Cell Collective to determine the association between external components (e.g., in-vitro inducers) and internal components (such as TLRs, cells, cytokines, and phenotypes). We used probabilistic global sensitivity analysis based on standardized regression coefficient (SRC) using the “sensitivity” package in R (24, 53) on input and output data of Cell Collective. In a single-input setting, we used SRC, which measures the strength of association between dependent and independent variables (53). We performed Cell Collective simulations under input activity levels ranging from 0 to 100%. The activity levels of inputs and outputs were independent and dependent variables in the statistical model. A higher SRC value represents a higher strength of association between input and output variables. We used SRC and k-means clustering algorithms (900 samples, specified three clusters for low, medium, and high activity levels) methods to visualize the simulation results.
Kyoto Encyclopedia of Genes and Genomes pathways analysis
The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway (54) enrichment analysis was conducted using the DAVID Bioinformatics Resources (2021 Update) (55) to explore the model components at the functional level. DAVID is a gene functional classification tool in which we used a p-value <0.05 with a false discovery rate (< 5%) as the cutoff criterion for KEGG pathway enrichment. We used the ggplot2 R package to visualize the fold enrichment and P-values of the top 20 enriched KEGG pathways.
Results
A large-scale multicellular, mechanistic model of signal transduction regulation of dendritic cell immune responses
We constructed a mechanistic multiscale model of signal transduction networks governing the proper function of human DCs, spanning biological scales from molecular to cell-to-cell communication. The model comprises 281 components and 702 interactions between these components that regulate the DC responses. The multiscale nature of our model is based on the intercellular and intracellular communication dynamics between DCs and four other immune cell types, serving as a bridge between innate and adaptive immunity (56).
To facilitate the model’s utility, its architecture across various biological scales is first depicted in a biological illustration of the pathways, communication molecules, cell markers, and receptors involved in regulating DCs from immature cells to mature phenotypes (Figure 1).
Figure 1 Schematic representation summarizing the main components of the model and the connection between each biological layer of the system.
Further, the model’s architecture is depicted in Figure 2. Herein, the organization follows communication from the DC’s extracellular space to the plasma membrane (ligand-receptors, markers), to the cytoplasm (kinases and signaling cascades), to the nucleus (transcription factors and gene regulations) and to the secretory compartments that communicate with various DC phenotypes and other immune cell types that interact with DCs (Figure 2; Table S1 ). The model is freely available in the Cell Collective modeling software and repository (45, 57) (DC model Link and workflow to Cell Collective environment in Supplemental Figure S2). Each component and interaction of the model have been fully annotated to facilitate data transparency and reusability. In Cell Collective, the community can also directly simulate and analyze the model, further improve it, or download it as an SBML file (46).
Figure 2 Visualization of the large-scale DC model in Cell Collective. The network view of the model. Dots represent signaling molecules; edges represent interactions between the model components. Red edges represent the inhibitions, while green and gray edges show activatory interactions. External (stimuli) and internal components are colored yellow and gray, respectively.
The model includes main signaling pathways, immune checkpoints, cytokines, and DC response mechanisms, and is able to represent the DC antigen presentation and maturation functions. The intra-cellular scale includes receptors that can sense the cellular environment and the downstream pathways that regulate DC responses that result in different cellular phenotypes. A diverse group of pattern recognition receptors (PRRs) is included in the model, including toll-like receptors (TLRs; TLR1, 2, 3, 4, 7, 8, and 9), C-type lectin receptors (CLRs; CLEC4C, CLEC7A, CLEC9A, CLEC10A, CD209), nucleotide-binding domain/leucine-rich repeat-containing receptors (NLRs; NOD2). The model also contains regulatory proteins (NF-κB, PYCARD, NLRP3, SYK, LILRA4, ISG20, ADAR, BST2, DDIT3, ATF4, PPP1R15A), maturation molecules (CD80, CD83, CD86, CD40, HLA-DR), and cytokines (IL6, TNF, IL12, IFNA1, IL12A, IL1B, IL12B, IL23A, IL10, IFNB1) (58).
In our study, we differentiated between various DC subtypes using a combination of subtype-specific markers. For pDCs, we employed CLEC4C (C-type lectin domain family 4 member C), and NRP1 (neuropilin 1) (59–61). For cDC1, we utilized CLEC9A (Dendritic cell C-type lectin receptor 9A), XCR1 (XC chemokine receptor 1), and THBD (thrombomodulin) (59, 60, 62–64). cDC2 was characterized using CLEC10A (C-type lectin domain containing 10A), and CD163 (CD163 molecule) (59, 60, 65). For MoDCs, we applied MRC1 (mannose receptor C-type 1), and CD1A (CD1a molecule) (59, 60, 66).
The intracellular molecular cascades stimulate interaction with other immune cells, including effector and exhausted T cells, B cells, NK cells, and neutrophils. For example, DCs capture and display antigen protein fragments on their plasma membrane through the antigen presentation process (67) then antigenic peptides are bound to appropriate molecules of the MHC, also known in humans as the human leukocyte antigen (HLA). T cells can recognize the antigens at the T cell-APC interface. DCs’ highly stimulatory and versatile APC function produces cytokines, interferons (IFNB), and tumor necrosis factor superfamily (TNF) to stimulate naive T cells to differentiate into effector subsets (68). As such, cytokines (IL1, IL6, IL10, IL12, IL23) and IFNB are also included in the model. DCs increase the expression of the MHC, the adhesion molecules, and the co-stimulators upon maturation, further stimulating T-cell proliferation and cytokine release (69). DC immune checkpoints are also included in the model (ICOS-LG, TNFSF9, TNFSF4, CD70, PVR, Nectin-2, BTLA, and PD-L1); the checkpoints regulate stimulatory and inhibitory pathways capable of maintaining self-tolerance and facilitating the immune response (70).
In-silico model validation
To validate the mechanistic model of DC functions, we collected experimental data from 30 different studies (Supplemental Data, Table S2, Figure S1) and reproduced them via in-silico experiments. Below we describe four in-silico experiments that showcase how the model was used to validate well-known (extensively published) in-vitro experiments spanning intracellular communication dynamics, intercellular communication dynamics, and both inter and intracellular communication dynamics.
The first experiment assessed the model’s ability to reproduce the behavior of TLRs. Namely, Grandclaudon (26) studied a range of DC molecular states expressing various patterns of communication signals. The authors present that DCs were treated for 24 hours with lipopolysaccharide (LPS), which activated TLR4 signaling pathways and induced DC communication molecules, including IL1β, IL6, TNF-α, and IL12 cytokines. Similar studies (58, 71, 72) present the LPS-induced secretion of DC cytokines as a non-trivial test to investigate the TLR4 cooperation in response to infections. A better understanding of the mechanisms of host resistance can provide a basis for the development of more effective adjuvants and immunotherapeutic regimens.
To validate the ability of the model to reproduce TLR behavior, we first validated that the model contained all components and regulatory pathways to support this experiment. Next, we ran 1,000 dose-response simulations while defining three activity levels of LPS (0, 50, and 100). We then evaluated and compared the secretion of inflammatory cytokines – IL1β, IL6, IL12, IFN-α, and TNF – upon activation by LPS at each of the three different dose responses against well-established cytokine responses. Figure 3A displays the secretion level for each of the five cytokines at LPS doses 0, 50, and 100. As expected, there was no cytokine secretion at dose 0 and subsequent secretion and elevated secretion at doses 50 and 100, respectively.
Figure 3 In-silico model validations. (A) Inflammatory cytokine activity level in response to LPS using dose-response analysis. (B) Standard regression coefficient of main factors activated in the HIV infection environment. (C) Activity level of IL2, IFN-α and NK under CpG-containing oligonucleotides and poly(I:C) stimulation. (D) Markers of maturation at different time points in the presence of neutrophils. (E) Time course distribution of DC immature (im-DC) and mature (mDC) states with (+) or without (-) neutrophils.
The second experiment assessed the model’s ability to mount an appropriate immune response to the presence (and initiation) of the human immunodeficiency virus (HIV) infection. HIV initiates viral transcription through TLR8 and promotes the maturation of DCs (from immature (imDC) to mature plasmacytoid (pDC)) as defined by the expression of CD83 and CCR7 surface markers and the production of IFN-α and TNF (73, 74). We used Cell Collective’s global sensitivity analysis (refer to the “Methods” section, “Global sensitivity analysis”) as a method to determine the association between HIV and each of the internal components. Figure 3B displays the correlation of activity for imDC, pDC, CD83, CCR7, IFN-α, and TNF in the presence of HIV infection. As expected, immature DCs (imDC) exhibit a negative correlation in the presence of HIV, which shows HIV-bearing immature DCs can differentiate into mature DCs in response to the infection, presenting HIV antigens to T cells and initiating viral immune responses. Further, mature DCs (pDC), as well as surface markers CD83 and CCR7, and IFN-α and TNF exhibit a positive standard regression coefficient (SRC), which means an increase in HIV load results in increased activity of these components. The simulation results are consistent with biological experiments that describe pDCs exposed to HIV strongly upregulating the expression of CD83 and functional CCR7 maturation markers, IFN-α, and TNF cytokines (73).
The third experiment assessed the model’s ability to simulate known intercellular dynamic crosstalk between DCs and other immune cells in tandem with intracellular communication dynamics. Gerosa (75) showed that human peripheral pDC and MoDCs are necessary to induce NK cell function depending on the type of microbial stimulus. In this experiment, pDCs and MoDCs were stimulated in response to CpG-containing oligonucleotides (CpG) and poly(I:C), and evaluated the mean activity level for NK cells, IL2, IFN-α, as a result of (CpG)/poly(I:C)-induced release of IL2 and IFN-α and subsequent activation of NK cells. Figure 3C displays the expected activity of IL2 and IFN-α as well as NK cells when CpG/poly(I:C) is inactive compared to an active state.
Last, we validated intercellular communication dynamics between imDCs and neutrophils. Neutrophils stimulate imDCs to become competent antigen-presenting cells. This maturation phenotype is characterized by the expression of specific surface markers (e.g., HLA-DR, CD86, and CD40) and the secretion of IL12 in response to DC-neutrophil interactions (76–78). Figure 3D displays the activity of model components IL12, CD209, CD40, CD86, and HLA-DR in response to the presence of neutrophils over time, demonstrating the pathways responsible for neutrophil-induced DC maturation. Figure 3E shows the activity level of immature and mature DCs in the presence and absence of neutrophils. On the left, when neutrophils are absent, immature DCs continue to increase in activity over time, whereas mature DCs do not become active. On the right, as neutrophils become present, immature DC activity tapers, and mature DC activity increases.
The aforementioned experiments illustrate the ability of the model to reproduce major experiments spanning complex inter- and intracellular communication dynamics.
Case studies
To aid researchers in identifying how to use this model, we showcase three case studies by presenting a brief application background, the method we used to apply the model in this context, and model results.
Case 1: Intracellular communication dynamics. Characterization of DC response to a combinatorial COVID-19 and Influenza infection environment.
In this case study, we integrated Covid-19 and Influenza pathogens into the model to characterize the molecular response of DC under single and co-infection conditions. Coronavirus disease 2019 (COVID-19) and Influenza respiratory disease, caused by Sars-CoV-2 and influenza virus, respectively, share similarities in seasonal manifestations, viral transmission method, symptoms, and immunopathogenesis (79, 80). Co-infection with Sars-CoV-2 and influenza virus increases disease severity and impairs neutralizing antibody and CD4+ T cell responses (81). Patients can develop both infections, and in some cases, co-infection leads to a poor prognosis (82–84). Despite the comprehensive investigation of DC behavior in single infections with Sars-CoV-2 or influenza, the comparative understanding of DC programming under co-infection is not fully explored due to limited patient cohorts and case studies (85–88).
Thus, the purpose of this model-based study was to investigate the molecular behavior of DCs in three infectious states: infection with i) Influenza type A virus (IAV), ii) Sars-CoV-2, and ii) co-infection with IAV and Sars-CoV-2. We ran 900 dose-response simulations and identified three molecular patterns (Figure 4A, Table S3) based on similarity in activity levels of DCs molecular components. Each pattern presents a list of molecules that play a role in the co-infection. We reported the mean activity level of DC model molecules ranging from low activity (green) to fully activated (red) molecular state in single infection and co-infection cellular environments.
Figure 4 In-silico predictions of molecular activity across the whole DC model comparing Sars-CoV-2 and Influenza A virus (IAV) co-infection to the single infection. (A) The differential molecule expressions are grouped into three main patterns in response to each environmental setting. The first pattern grouped molecules that are regulated similarly in all three conditions. The second pattern is related to similar regulation between Influenza and co-infection, and the third one grouped similar behavior between Sars-CoV-2 and co-infection conditions. The scale represented the activity level ranging between 0 to 100%, 100 being the highest activity level. (B) The molecular signatures in pattern 1 in single and co-infection cellular environments. (C) An example of the second pattern shows that under different activity levels of Sars-CoV-2 (green, low; yellow, medium; purple, high), ICOSLG is inactive while it is upregulated in both co-infection and single IAV. (D) CLEC9A is categorized as the third pattern, and for both Sars-CoV-2 and co-infection, it has low activity levels compared to the IAV single infection.
In the first pattern (Figure 4A – Pattern 1 and 4B), we identified molecular signatures with similar activities in single infection and co-infection. Figure 4B presents molecular signatures following this pattern, including markers of DC differentiation (CD86, CD1A, CD40, CD83, ITGAM), PRRs (TLR8), immune checkpoint molecules (PVR, Nectin2), chemokines/chemokine receptor (CCR7, CCL19, CXCL8), cytokines (IL6, TNF, IL12, IL1B, IL12A, IL12B, IL23, IL10, IFNA1, IFNB1), signaling molecules (NF-κB, PYCARD, NLRP3, SYK, LILRA4, ISG20, ADAR, BST2, DDIT3, ATF4, PPP1R15A), and CLRs (CLEC4C, CLEC10A).
Several of these signatures are expressed during the single infection studies on human samples (89–91). For example, separate studies on Sars-CoV-2 and influenza virus infections show expression of IL1B, IL10, TNF, CD86, CCR7, IL6, CXCL8, IFN (79, 89, 92, 93).
In the second pattern (Figure 4A – Pattern 2), the molecular signature characterizes the similarity between IAV single infection and co-infection. Previous studies indicated that immune checkpoints are increased in influenza single infection (94) but not in Sars-CoV-2 single infection (89, 95, 96). Thus this experiment focuses on the significance of the immune checkpoint signatures. The immune checkpoints (TNFSF4, CD70, ICOSLG, PDCD1LG2), followed by cytokines (IL2, IFNL2, CXCL10), markers of DC differentiation (CD80, CD86), and signaling signature (SEMAD4), are upregulated in both co-infection and single IAV but downregulated under Sars-CoV-2 infection. As an example, Figure 4C shows the activity level of the ICOSLG immune checkpoint in 300 simulations per each infection condition (single and co-infection), which is higher in the presence of both viruses.
In the third pattern (Figure 4A – Pattern 3), the molecular signatures of Sars-CoV-2 and co-infection were similar. The major signature includes a decreased expression of signaling and decreased expression of transcription factors in both Sars-CoV-2 and co-infection, suggesting a disruption of the signaling network associated with Sars-CoV-2 infection. Neuropilin-1 (NRP1), the only signaling protein to be highly expressed in the third pattern, was previously shown to facilitate Sars-CoV-2 entry by interacting with spike protein (97, 98). Additional signatures are related to the decrease of pathogen sensors and maturation marker expressions, such as TLRs (TLR1, TLR7, MYD88), CLRs (CLEC9A, CLEC7A, CD209), and MHC class signatures (HLA-DQA, HLADPB1, HLA-DM, HLA-DRB1, HLA-DR, HLA-DQB1), suggesting the loss of DC function to sense and present antigen to other immune cells properly. For example, Figure 4D presents the simulation results of the C-type lectin domain containing 9A (CLEC9A) with a low activity level in co-infection.
Several studies indicated that DCs displayed a defect in maturation and are depleted in COVID-19 patients, and as our in-silico simulations predicted, one of the mechanisms might be due to the defect of the signaling compartment and pathogen sensors (85, 86, 99). Nevertheless, further experimental investigations are needed to explore these hypotheses.
Case 2: Intercellular communication dynamics. Crosstalk between DCs and T cells in a cancer microenvironment.
DCs play a crucial role in initiating a protective anti-tumoral response by presenting tumor antigens and providing co-stimulatory immune checkpoint to T cells (100). However, tumor microenvironments sustain DCs in an immature/tolerant phenotype, thereby altering antigen presentation, co-stimulatory signals, and thus the ability to effectively activate T cells. Therefore, T-cells become exhausted due to continuous exposure to antigens and increase in multiple inhibitory immune checkpoints that further benefit the mechanism of resistance to immunotherapies (101).
Several factors with immunoregulatory properties are involved in DC-T cell interplay. For example, the cytokine HMGB1 released by cancer cells contributes to cancer development by promoting tolerogenic DC differentiation and the suppression of anti-tumoral T cells (102–104). Moreover, a study conducted in-vivo reported the role of HMGB1 in promoting T-cell exhaustion in the condition of trauma (105). However, the role of cancer-derived HMGB1 in promoting exhaustion through the modulation of immune checkpoint expression has not been investigated. Modern immunotherapy approaches aim to reverse T cell exhaustion by blocking inhibitory immune checkpoint receptors (106). Combinatorial treatments using approved inhibitors of PD-L1 immune checkpoint and two receptors PD-1 and CTLA-4, showed promising results. Additional immune checkpoint inhibitors are under clinical trial investigations (107). However, not all cancer types respond equally, and patients can acquire resistance to immune checkpoint inhibitors (ICI) (108, 109). Because many experimental studies have investigated the role of immune checkpoints individually, a computational approach can help to better understand the dynamic distribution of inhibitory and stimulatory immune checkpoints that can aid in identifying ideal checkpoint candidates and facilitate combinatorial therapeutic strategies.
In this case study, we examined the impact of cancer-derived HMGB1 on the DC-T cell synapse interaction. We included the cancerous cytokine HMGB1 as environmental (cancer) stimulus. The DC model includes two groups of immune checkpoints: stimulatory and inhibitory ligands/receptors that are enable us to study the impact of HMGB1 on DC-T immune checkpoints. The ligands are expressed on DCs, while receptors are particularly expressed by T cells. In the model’s plasma membrane compartment (Figure 5A, Table S1) are included six stimulatory ligands (e.g., ICOS-LG, TNFSF9, TNFSF4, CD70), two inhibitory molecules (PD-L1, BTLA) and three molecules with a dual function depending on the receptors they are binding (CD80-CD86, PVR, and Nectin-2). From the T cell side, we included six stimulatory receptors (CD28, ICOS, TNFRSF9, TNFRSF4, CD27, CD226) to define effector T cells and three inhibitory receptors (PD-1, TNFRSF14, CTLA-4) that define exhausted T cells. Figure 5C showed the interaction between ligands with their respective receptors. Of note, CD80-CD86, used as main maturation markers, binds two different immune checkpoint receptors with opposite functions (the stimulatory receptor CD28 and the inhibitory receptor CTLA4), and the two ligands PVR and Nectin-2 share the same stimulatory receptor CD226 (Figure 5A).
Figure 5 Investigation of DCs-T cells crosstalk under an HMGB1 tumor environment. (A) The table indicates the immune checkpoint pairing between ligands and associated receptors along with the type of functions (stimulatory or inhibitory). (B) Assessment of IL6, IL8, IL12, and TNF cytokine expression under HMGB1 environment. (C) The activity level of MHC classes I and II are between 0 and 10, with the time 100 added at the final expression of both classes. (D) Standardized regression coefficient (SRC) between ligands and associated receptors on DC and T cells with HMGB1 environment. High and low SRC presented with yellow and purple, respectively. Arrows link the ligands to their respective receptors. (E) Time course of effector T and exhausted T cells activity level expression at a probability of activation at time steps 0, 1, 10, and 100 in HMGB1 simulation.
We compared the model’s simulation results under the HMGB1/cancer environment with published experimental data (110). Messmer D. et al. showed that HMGB1 promotes the secretion of inflammatory cytokines (IL6, IL12, IL8, TNFα) and our in-silico simulation is consistent with the experimental data (Figure 5B). Then, in the HMGB1 environment, we evaluated: i) the activation of MHC Class I and II, ii) the distribution of ligands/receptors for co-stimulatory and inhibitory immune checkpoints expressed at the DC-T cell interface, and iii) the dynamic distribution of effector and exhausted T cells.
We evaluated the distribution of the mean percentage activity level of MHC class I and II in five in-silico experiments at time steps 0, 1, 10, and 100. The mean activity level started from 0, and at time step 10 reaches 32.74% for class I and 61.8% for class II. MHC class I and II increase to maximum activity level at step 100 (Figure 5C). Our simulation indicated that MHC class I and II expression increased in response to HMGB1 simulation. Because our model does not include specific tumor antigens that can be presented by MHC to TCR, we cannot conclude that the increase of MHC expression is due to antigen overload, however, our simulation indicates an increase of MHC classes under HMGB1 simulation.
Next, we investigated the dynamics of immune checkpoint pairs under HMGB1 environmental stimulation using dose-response and sensitivity analyses (Figure 5D). Our in-silico results showed that the stimulatory molecules such as CD27 (receptor for CD70), ICOS (receptor for ICOS-LG), TNFRSF4 (receptor for TNFSF4), and TNFRSF9 (receptor of TNFSF9) showed no significant correlation with HMGB1 stimulation and shared similar distribution with their receptors expressed by T cells. The PVR and Nectin-2 displayed high correlation as well as their receptor CD226. CD80-CD86 showed no significant correlation in response to HMGB1 stimulation (Figure 5D). Regarding the dual receptors of CD80-CD86, the stimulatory CD28 receptor, and the inhibitory receptor CTLA-4 showed no correlation under HMGB1 stimuli. The inhibitory pairing PD-L1 ligand with its receptor PD-1, the main target for immune checkpoint inhibitors, is highly represented. Moreover, the additional inhibitory receptor TNFRSF14 and its ligand BTLA ligand don’t show a significant distribution in response to HMGB1 stimulation (Figure 5D).
We simulated the model under HMGB1 environment, and we evaluated the mean activity level of effector and exhausted T cells at time 0’, 1’, 10’, and 100’ from five in-silico experiments (Figure 5E). The exhausted T cell activity level is faster than the effector T cells at 10’ (Figure 5E). At the maximum time of the simulation, both T cell phenotypes do reach maximum accumulation (activity). The results demonstrate HMGB1 promotes both effectors and exhausted T cells and exhausted phenotype accumulated faster than effector.
Immune checkpoint immunomodulatory functions are initiated by ligand-receptor interaction that can either promote or suppress T cell function (111). CD226 is important in generating an anti-tumor response. While CD226 expression is required as a co-stimulatory factor for T cells during antigen presentation by APCs, the loss of CD226 can lead to impaired effector T activation and increased susceptibility to tumor development (112–114). Hence, in our model, the activity of CD226-PVR/Nectin2 contributed to the increase of effector T cells and is associated with the MHC expression. Among inhibitor pairs, only PD-L1-PD-1 displays a high correlation with an HMGB1 simulation. The interaction between PD-L1 and PD-1 drives T-cell dysfunction and exhaustion to prevent an efficient anti-tumor T-cell response (115, 116). Previous studies indicated that HMGB1 increases PD-L1 expression in cancer cells; however, the modulation of PD-L1-PD-1 by HMGB1 in immune cells remains unknown. Our in-silico simulation suggests that HMGB1 can also promote PD-L1-PD-1 expression at the DC-T cell interface, thereby explaining the increase of exhausted T cells.
In summary, using the example of analysis of multiple ligand/receptor-mediated cellular programming at the time, our in-silico experiments illustrated the capacity of the model to provide complex and dynamic insight into biological processes at the molecular and cellular scales.
Case 3: The scope of the DC models offers potential applications in several immune-related diseases.
The crosstalk between the disease environment and DCs highly contributes to the organization of the immune response (11, 12, 117–119). Because each disease environment is unique and complex, a multiscale model can be an effective tool to investigate the complexities underlying multiscale, systemic diseases.
Given the DCs’ role in initiating both innate and adaptive immune responses, we sought to explore the links between the DC model’s core disease pathways we identified for IAV, Sars-CoV-2, and tumor microenvironment and additional diseases to identify the extensibility of our model. To do this, we performed a Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis against the model’s components (refer to “Methods” section, “Kyoto Encyclopedia of Genes and Genomes KEGG pathways analysis”). Using a cut-off p-value of <0.05, we identified 69 enriched pathways (Table S4) for the top 20 diseases (Figure 6). The fold enrichment analysis of the top 20 diseases revealed multiple categories, such as autoimmune disease, infection, and transplantation.
Figure 6 The top 20 human diseases and signaling pathways associated with the DC model. P-values and fold enrichment of KEGG pathways in DC. -log 10 of P-values was used for visualization. Thus, larger dot sizes correspond to lower P-values. The cutoffs P-value < 0.05 and a false discovery rate (FDR) < 5% were set for significant enrichment.
The two highest scores are represented by inflammatory bowel disease and tuberculosis infection. The enrichment for those two diseases can be explained by the presence of TLRs, lectins, and cytokines in supporting pathology (120). For example, in the context of tuberculosis, TLRs and lectins can recognize different motifs of Mycobacterium tuberculosis, which in turn can trigger pro- or anti-inflammatory cytokine response (121).
The DC model can be utilized to identify pathways related to many diseases, including infections, cancer, and autoimmune diseases. The high fidelity of the model predictions will depend on the extension of the pathways related to the diseases of study.
Discussion
We have developed a mechanistic multiscale model of human DCs that spans biological scales from molecular interactions to cell-cell communication. We included biological events that occur between DCs’ environmental stimuli and their receptors, followed by activation of signal transduction in response to each signal. Moreover, we constructed the molecular network that links the downstream signal transduction of kinases and transcription factors to secreted cytokines/chemokines and growth factors. We extended the model further by integrating a cellular compartment that includes the communication between DCs and several innate/adaptive immune cells through direct (ligand-receptor) and indirect (cytokine, chemokines releases) interactions. Our model can be used to study DC maturation, differentiation to each subset, APCs function, and the bidirectional crosstalk between DCs and other immune cells. Because the model incorporates pathways that regulate and facilitate many key functions of DCs, it can be applied to study several diseases as well as the basic mechanism of DCs’ functions.
The presented DC model leverages the widely used logical modeling formalism (40). The advantages of this modeling approach include its scalability (efficient simulations) as evidenced by the fact that some of the largest computational models have been constructed using this formalism (e.g (122–125). Another advantage is that logical models do not rely on kinetic parameters that are mostly unknown (40, 126). On the other hand, if one is interested in modeling relatively small and well-studied pathways (with known parameters), a kinetic modeling approach may be more appropriate. The model is limited by the missing data in the literature about any unknown interactions. Our model includes major pathways involved in DC immunobiology. Nevertheless, the model is limited in scope as it does not include all known DC signaling and cell-cell communication. The model is being provided in a readily exchangeable format (SBML) and easy-to-use modeling software (Cell Collective), making it relatively easy for the community to build on the model and continue to expand as needed by different applications. For example, to specifically investigate DC-T cell communication, T cell subsets such as CD4 and CD8 can be integrated by adding molecular and cellular components of the immunological synapse. We previously published a logical model of signal transduction networks governing CD4+ T cell differentiation in response to various cytokines (24). Those same cytokines are also included in our DC model, creating the possibility of integrating both systems to study how DCs might influence CD4+ T cell fate and plasticity.
As another example, HMGB1 interacts with several TLRs (e.g., TLR2, TLR4, and TLR9), which have been included in the DC model. HMGB1 also interacts with RAGE - a receptor for advanced glycation end-products - that is not currently included in our model. Adding RAGE to the system would increase the complex interplay between receptors and signaling pathways to mediate cytokine release and immune response (127–129). The model would then be able to simulate the different molecular intersections during single or multiple TLRs/RAGE activation and predict the multiple environmental conditions for efficient DCs maturation without compromising the adaptive response (e.g., T and B cells). Therefore, the multiscale model could be further used to characterize APC function in response to a stochastic tumor micro-environment with multiple components simultaneously.
In our cancer in-silico simulations, our model-generated hypotheses suggested a list of potential immune checkpoints to explore for studying the effect of single and multiple combinatorial ICI on DC-T cell interaction outcome (Figure 5D). We showed the dynamics of immune checkpoint pairs under a tumor HMGB1 environment. Recent therapeutic approaches include the optimization of DC-based therapies by combining DC vaccines with immune checkpoint inhibitors (ICI), such as anti-CTLA-4 and anti-PD1/PDL1 (130, 131), or by silencing immune checkpoint signaling pathways (132). Despite being in early clinical phases, combinatorial therapy holds a potential to balance toxicity, safety, and clinical outcomes (130, 131). Additional ICI to restore T cell or APC activation is currently under investigation to expand therapeutic options and optimize the efficacy of the immune checkpoint targeting strategy (107, 133). Nevertheless, the complexity of immune checkpoint ligands resides in their capacity to bind several different receptors with opposite functions, therefore switching between stimulatory and inhibitory signals. As the model prediction suggested, PVR and Nectin-2 showed a high activity similar to their receptor, CD226. Of note, PVR and Nectin-2 can trigger opposite signals whether they bind the stimulatory receptor (CD226) or the inhibitory receptors (TIGIT and CD96, not included in the model) (134). Moreover, the optimal combination can depend on the ligands/receptors’ availability and the balance between stimulatory and inhibitory expression. Our model simulations suggest that the inhibitory receptor CTLA-4 has no activity under HMGB1 stimuli. At the same time, PD-1 and PD-L1 are highly correlated, suggesting that the use of anti-CTLA-4 might not be as effective as the use of anti-PD-1 or anti-PD-L1 to restore DC-T cell function in a cancerous HMGB1 environment (135).
The development of computational models that recapitulate complex human disease behavior can be a resource for scientists and clinicians to simulate thousands of possibilities for studying the complex biological process at multiple scales. The disease enrichment analysis highlighted the potential of our model to incorporate additional pathological events as some disease modules are already implemented. For example, Type I diabetes (T1DM), an auto-immune disease characterized by immune-mediated destruction of insulin-producing beta cells, is enriched in our model (136). The loss of tolerance to self-antigens and the increase of autoreactive T cells instead of immunosuppressive T cells are the main cause of insulin deficiency. Several studies indicated that DCs presented self-antigen generated from degraded b-islet to prime autoreactive T cells via dysfunctional NF-κB and MAPK pathways (137, 138). Current therapies focus on generating tolerant DCs and immunosuppressive T cells to target the auto-immune disease and restore the imbalance of tolerance (139). To address those mechanisms, incorporating tolerogenic DCs and immunosuppressive T cell phenotype components under the stimulation of a self-antigen input could predict molecular conditions by which immunosuppressive cells are amplified to respond to disease pathology (140).
In summary, we have demonstrated the potential for a multiscale DC model to investigate the immunobiology of DCs and identify potential targets for improving the effectiveness of DC-based cell therapies. Lastly, the model can be further expanded to support additional mechanistic and therapeutic questioning related to DC ontogeny.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.
Author contributions
SSA, RA and TH conceived the study. SSA, RA and TH designed the study. SSA performed literature mining and collected the data. SSA constructed the models. SSA, and RA performed refinement of the constructed models. SSA, RA and BLP analyzed the data, performed the experimental work and analyzed the experimental results. SSA, BLP, RA and TH wrote the manuscript. RA and TH supervised the study. All authors contributed to the article and approved the submitted version.
Funding
The work was supported by an NIH grant R35GM119770 to TH.
Conflict of interest
TH is the majority stakeholder in Discovery Collective, Inc. with proprietary rights to Cell Collective. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2023.1112985/full#supplementary-material
References
1. Joffre O, Nolte MA, Spörri R, Reis e Sousa C. Inflammatory signals in dendritic cell activation and the induction of adaptive immunity. Immunol Rev (2009) 227:234–47. doi: 10.1111/j.1600-065X.2008.00718.x
2. Pasare C, Medzhitov R. Toll-like receptors: linking innate and adaptive immunity. Microbes Infect (2004) 6:1382–7. doi: 10.1016/j.micinf.2004.08.018
3. Pandey S, Singh S, Anang V, Bhatt AN, Natarajan K, Dwarakanath BS. Pattern recognition receptors in cancer progression and metastasis. Cancer Growth Metastasis (2015) 8:25–34. doi: 10.4137/CGM.S24314
4. de Jong JMH, Schuurhuis DH, Ioan-Facsinay A, Welling MM, Camps MGM, van der Voort EIH, et al. Dendritic cells, but not macrophages or b cells, activate major histocompatibility complex class II-restricted CD4+ T cells upon immune-complex uptake. vivo Immunol (2006) 119:499–506. doi: 10.1111/j.1365-2567.2006.02464.x
5. Kambayashi T, Laufer TM. Atypical MHC class II-expressing antigen-presenting cells: can anything replace a dendritic cell? Nat Rev Immunol (2014) 14:719–30. doi: 10.1038/nri3754
6. Figdor CG, de Vries IJM, Lesterhuis WJ, Melief CJM. Dendritic cell immunotherapy: mapping the way. Nat Med (2004) 10:475–80. doi: 10.1038/nm1039
7. Fricke I, Gabrilovich DI. Dendritic cells and tumor microenvironment: a dangerous liaison. Immunol Invest (2006) 35:459–83. doi: 10.1080/08820130600803429
8. Sesti-Costa R, de Moraes-Vieira PMM, Cervantes-Barragan L. Dendritic cells: immune response in infectious diseases and autoimmunity. Mediators Inflammation (2020) 2020:2948525. doi: 10.1155/2020/2948525
9. Audiger C, Rahman MJ, Yun TJ, Tarbell KV, Lesage S. The importance of dendritic cells in maintaining immune tolerance. J Immunol (2017) 198:2223–31. doi: 10.4049/jimmunol.1601629
10. Domogalla MP, Rostan PV, Raker VK, Steinbrink K. Tolerance through education: How tolerogenic dendritic cells shape immunity. Front Immunol (2017) 8:1764. doi: 10.3389/fimmu.2017.01764
11. Ganguly D, Haak S, Sisirak V, Reizis B. The role of dendritic cells in autoimmunity. Nat Rev Immunol (2013) 13:566–77. doi: 10.1038/nri3477
12. Chow TG, Gill MA. Regulation of allergic inflammation by dendritic cells. Curr Opin Allergy Clin Immunol (2020) 20:56–63. doi: 10.1097/ACI.0000000000000603
13. Lanzavecchia A, Sallusto F. Regulation of T cell immunity by dendritic cells. Cell (2001) 106:263–6. doi: 10.1016/s0092-8674(01)00455-x
14. Thaiss CA, Semmling V, Franken L, Wagner H, Kurts C. Chemokines: a new dendritic cell signal for T cell activation. Front Immunol (2011) 2:31. doi: 10.3389/fimmu.2011.00031
15. Wieczorek M, Abualrous ET, Sticht J, Álvaro-Benito M, Stolzenberg S, Noé F, et al. Major histocompatibility complex (MHC) class I and MHC class II proteins: conformational plasticity in antigen presentation. Front Immunol (2017) 8:292. doi: 10.3389/fimmu.2017.00292
16. Reis e Sousa C. Activation of dendritic cells: translating innate into adaptive immunity. Curr Opin Immunol (2004) 16:21–5. doi: 10.1016/j.coi.2003.11.007
17. Geginat J, Nizzoli G, Paroni M, Maglie S, Larghi P, Pascolo S, et al. Immunity to pathogens taught by specialized human dendritic cell subsets. Front Immunol (2015) 6:527. doi: 10.3389/fimmu.2015.00527
18. Eisenbarth SC. Dendritic cell subsets in T cell programming: location dictates function. Nat Rev Immunol (2019) 19:89–103. doi: 10.1038/s41577-018-0088-1
19. Banchereau J, Paczesny S, Blanco P, Bennett L, Pascual V, Fay J, et al. Dendritic cells: controllers of the immune system and a new promise for immunotherapy. Ann N Y Acad Sci (2003) 987:180–7. doi: 10.1111/j.1749-6632.2003.tb06047.x
20. Mastelic-Gavillet B, Balint K, Boudousquie C, Gannon PO, Kandalaft LE. Personalized dendritic cell vaccines-recent breakthroughs and encouraging clinical results. Front Immunol (2019) 10:766. doi: 10.3389/fimmu.2019.00766
21. Yu J, Sun H, Cao W, Song Y, Jiang Z. Research progress on dendritic cell vaccines in cancer immunotherapy. Exp Hematol Oncol (2022) 11:3. doi: 10.1186/s40164-022-00257-2
22. Jähnisch H, Füssel S, Kiessling A, Wehner R, Zastrow S, Bachmann M, et al. Dendritic cell-based immunotherapy for prostate cancer. Clin Dev Immunol (2010) 2010:517493. doi: 10.1155/2010/517493
23. Cintolo JA, Datta J, Mathew SJ, Czerniecki BJ. Dendritic cell-based vaccines: barriers and opportunities. Future Oncol (2012) 8:1273–99. doi: 10.2217/fon.12.125
24. Puniya BL, Todd RG, Mohammed A, Brown DM, Barberis M, Helikar T. A mechanistic computational model reveals that plasticity of CD4+ T cell differentiation is a function of cytokine composition and dosage. Front Physiol (2018) 9:878. doi: 10.3389/fphys.2018.00878
25. Wertheim KY, Puniya BL, La Fleur A, Shah AR, Barberis M, Helikar T. A multi-approach and multi-scale platform to model CD4+ T cells responding to infections. PloS Comput Biol (2021) 17:e1009209. doi: 10.1371/journal.pcbi.1009209
26. Grandclaudon M, Perrot-Dockès M, Trichot C, Karpf L, Abouzid O, Chauvin C, et al. A quantitative multivariate model of human dendritic cell-T helper cell communication. Cell (2019) 179:432–447.e21. doi: 10.1016/j.cell.2019.09.012
27. Verma M, Bassaganya-Riera J, Leber A, Tubau-Juni N, Hoops S, Abedi V, et al. High-resolution computational modeling of immune responses in the gut. Gigascience (2019) 8:giz062. doi: 10.1093/gigascience/giz062
28. Lai X, Keller C, Santos G, Schaft N, Dörrie J, Vera J. Multi-level computational modeling of anti-cancer dendritic cell vaccination utilized to select molecular targets for therapy optimization. Front Cell Dev Biol (2021) 9:746359. doi: 10.3389/fcell.2021.746359
29. Qiao L, Phipps-Yonas H, Hartmann B, Moran TM, Sealfon SC, Hayot F. Immune response modeling of interferon beta-pretreated influenza virus-infected human dendritic cells. Biophys J (2010) 98:505–14. doi: 10.1016/j.bpj.2009.10.049
30. Klinke DJ. A multi-scale model of dendritic cell education and trafficking in the lung: implications for T cell polarization. Ann BioMed Eng (2007) 35:937–55. doi: 10.1007/s10439-007-9318-6
31. Mei Y, Abedi V, Carbo A, Zhang X, Lu P, Philipson C, et al. Multiscale modeling of mucosal immune responses. BMC Bioinf (2015) 16(Suppl 12):S2. doi: 10.1186/1471-2105-16-S12-S2
32. Schwab JD, Kühlwein SD, Ikonomi N, Kühl M, Kestler HA. Concepts in boolean network modeling: What do they all mean? Comput Struct Biotechnol J (2020) 18:571–82. doi: 10.1016/j.csbj.2020.03.001
33. Samaga R, Klamt S. Modeling approaches for qualitative and semi-quantitative analysis of cellular signaling networks. Cell Commun Signal (2013) 11:43. doi: 10.1186/1478-811X-11-43
34. Albert R, Thakar J. Boolean modeling: a logic-based dynamic approach for understanding signaling and regulatory networks and for making useful predictions. Wiley Interdiscip Rev Syst Biol Med (2014) 6:353–69. doi: 10.1002/wsbm.1273
35. Le Novère N. Quantitative and logic modelling of molecular and gene networks. Nat Rev Genet (2015) 16:146–58. doi: 10.1038/nrg3885
36. Hemedan AA, Niarakis A, Schneider R, Ostaszewski M. Boolean modelling as a logic-based dynamic approach in systems medicine. Comput Struct Biotechnol J (2022) 20:3161–72. doi: 10.1016/j.csbj.2022.06.035
37. Siegle L, Schwab JD, Kühlwein SD, Lausser L, Tümpel S, Pfister AS, et al. A boolean network of the crosstalk between IGF and wnt signaling in aging satellite cells. PloS One (2018) 13:e0195126. doi: 10.1371/journal.pone.0195126
38. Béal J, Pantolini L, Noël V, Barillot E, Calzone L. Personalized logical models to investigate cancer response to BRAF treatments in melanomas and colorectal cancers. PloS Comput Biol (2021) 17:e1007900. doi: 10.1371/journal.pcbi.1007900
39. Puniya BL, Allen L, Hochfelder C, Majumder M, Helikar T. Systems perturbation analysis of a Large-scale signal transduction model reveals potentially influential candidates for cancer therapeutics. Front Bioeng Biotechnol (2016) 4:10. doi: 10.3389/fbioe.2016.00010
40. Abou-Jaoudé W, Traynard P, Monteiro PT, Saez-Rodriguez J, Helikar T, Thieffry D, et al. Logical modeling and dynamical analysis of cellular networks. Front Genet (2016) 7:94. doi: 10.3389/fgene.2016.00094
41. Niarakis A, Helikar T. A practical guide to mechanistic systems modeling in biology using a logic-based approach. Brief Bioinf (2021) 22:bbaa236. doi: 10.1093/bib/bbaa236
42. Niarakis A, Bounab Y, Grieco L, Roncagalli R, Hesse A-M, Garin J, et al. Computational modeling of the main signaling pathways involved in mast cell activation. Curr Top Microbiol Immunol (2014) 382:69–93. doi: 10.1007/978-3-319-07911-0_4
43. Biane C, Delaplace F, Klaudel H. Networks and games for precision medicine. BioSystems (2016) 150:52–60. doi: 10.1016/j.biosystems.2016.08.006
44. Tweedie S, Braschi B, Gray K, Jones TEM, Seal RL, Yates B, et al. Genenames.org: the HGNC and VGNC resources in 2021. Nucleic Acids Res (2021) 49:D939–46. doi: 10.1093/nar/gkaa980
45. Helikar T, Kowal B, McClenathan S, Bruckner M, Rowley T, Madrahimov A, et al. The cell collective: toward an open and collaborative approach to systems biology. BMC Syst Biol (2012) 6:96. doi: 10.1186/1752-0509-6-96
46. Chaouiya C, Bérenguier D, Keating SM, Naldi A, van Iersel MP, Rodriguez N, et al. SBML qualitative models: a model representation format and infrastructure to foster interactions between qualitative modelling formalisms and tools. BMC Syst Biol (2013) 7:135. doi: 10.1186/1752-0509-7-135
47. Keating SM, Waltemath D, König M, Zhang F, Dräger A, Chaouiya C, et al. SBML level 3: an extensible format for the exchange and reuse of biological models. Mol Syst Biol (2020) 16:e9110. doi: 10.15252/msb.20199110
48. Helikar T, Konvalina J, Heidel J, Rogers JA. Emergent decision-making in biological signal transduction networks. Proc Natl Acad Sci USA (2008) 105:1913–8. doi: 10.1073/pnas.0705088105
49. Helikar T, Rogers JA. ChemChains: a platform for simulation and analysis of biochemical networks aimed to laboratory scientists. BMC Syst Biol (2009) 3:58. doi: 10.1186/1752-0509-3-58
50. Garg A, Di Cara A, Xenarios I, Mendoza L, De Micheli G. Synchronous versus asynchronous modeling of gene regulatory networks. Bioinformatics (2008) 24:1917–25. doi: 10.1093/bioinformatics/btn336
51. Steinman RM, Hawiger D, Liu K, Bonifaz L, Bonnyay D, Mahnke K, et al. Dendritic cell function in vivo during the steady state: a role in peripheral tolerance. Ann N Y Acad Sci (2003) 987:15–25. doi: 10.1111/j.1749-6632.2003.tb06029.x
52. Kim R, Emi M, Tanabe K. Functional roles of immature dendritic cells in impaired immunity of solid tumour and their targeted strategies for provoking tumour immunity. Clin Exp Immunol (2006) 146:189–96. doi: 10.1111/j.1365-2249.2006.03215.x
53. CRAN - package sensitivity. Available at: https://cran.r-project.org/web/packages/sensitivity/index.html (Accessed February 9, 2023).
54. Kanehisa M, Furumichi M, Tanabe M, Sato Y, Morishima K. KEGG: new perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res (2017) 45:D353–61. doi: 10.1093/nar/gkw1092
55. Huang DW, Sherman BT, Tan Q, Kir J, Liu D, Bryant D, et al. DAVID bioinformatics resources: expanded annotation database and novel algorithms to better extract biology from large gene lists. Nucleic Acids Res (2007) 35:W169–75. doi: 10.1093/nar/gkm415
56. Palucka K, Banchereau J. Dendritic cells: a link between innate and adaptive immunity. J Clin Immunol (1999) 19:12–25. doi: 10.1023/a:1020558317162
57. Helikar T, Kowal B, Rogers JA. A cell simulator platform: the cell collective. Clin Pharmacol Ther (2013) 93:393–5. doi: 10.1038/clpt.2013.41
58. Trinchieri G, Sher A. Cooperation of toll-like receptor signals in innate immune defence. Nat Rev Immunol (2007) 7:179–90. doi: 10.1038/nri2038
59. Schlitzer A, Zhang W, Song M, Ma X. Recent advances in understanding dendritic cell development, classification, and phenotype. [version 1; peer review: 2 approved]. F1000Res (2018) 7:F1000 Faculty Rev–1558. doi: 10.12688/f1000research.14793.1
60. Collin M, Bigley V. Human dendritic cell subsets: an update. Immunology (2018) 154:3–20. doi: 10.1111/imm.12888
61. Alculumbre S, Raieli S, Hoffmann C, Chelbi R, Danlos F-X, Soumelis V. Plasmacytoid pre-dendritic cells (pDC): from molecular pathways to function and disease association. Semin Cell Dev Biol (2019) 86:24–35. doi: 10.1016/j.semcdb.2018.02.014
62. Schraml BU, Reis e Sousa C. Defining dendritic cells. Curr Opin Immunol (2015) 32:13–20. doi: 10.1016/j.coi.2014.11.001
63. Swartz AM, Nair SK. The In vitro differentiation of human CD141+CLEC9A+ dendritic cells from mobilized peripheral blood CD34+ hematopoietic stem cells. Curr Protoc (2022) 2:e410. doi: 10.1002/cpz1.410
64. Sasaki I, Kato T, Hemmi H, Fukuda-Ohta Y, Wakaki-Nishiyama N, Yamamoto A, et al. Conventional type 1 dendritic cells in intestinal immune homeostasis. Front Immunol (2022) 13:857954. doi: 10.3389/fimmu.2022.857954
65. Rhodes JW, Tong O, Harman AN, Turville SG. Human dendritic cell subsets, ontogeny, and impact on HIV infection. Front Immunol (2019) 10:1088. doi: 10.3389/fimmu.2019.01088
66. Patente TA, Pinho MP, Oliveira AA, Evangelista GCM, Bergami-Santos PC, Barbuto JAM. Human dendritic cells: their heterogeneity and clinical application potential in cancer immunotherapy. Front Immunol (2018) 9:3176. doi: 10.3389/fimmu.2018.03176
67. Heath WR, Carbone FR. Cross-presentation in viral immunity and self-tolerance. Nat Rev Immunol (2001) 1:126–34. doi: 10.1038/35100512
68. Pulendran B, Tang H, Manicassamy S. Programming dendritic cells to induce T(H)2 and tolerogenic responses. Nat Immunol (2010) 11:647–55. doi: 10.1038/ni.1894
69. Cella M, Scheidegger D, Palmer-Lehmann K, Lane P, Lanzavecchia A, Alber G. Ligation of CD40 on dendritic cells triggers production of high levels of interleukin-12 and enhances T cell stimulatory capacity: T-T help via APC activation. J Exp Med (1996) 184:747–52. doi: 10.1084/jem.184.2.747
70. Kubo T, Hirohashi Y, Torigoe T. [Dendritic cell and cancer immune checkpoint]. Nihon Rinsho Meneki Gakkai Kaishi (2016) 39:468–72. doi: 10.2177/jsci.39.468
71. Stamatos NM, Carubelli I, van de Vlekkert D, Bonten EJ, Papini N, Feng C, et al. LPS-induced cytokine production in human dendritic cells is regulated by sialidase activity. J Leukoc Biol (2010) 88:1227–39. doi: 10.1189/jlb.1209776
72. Yu L, Wang L, Chen S. Endogenous toll-like receptor ligands and their biological significance. J Cell Mol Med (2010) 14:2592–603. doi: 10.1111/j.1582-4934.2010.01127.x
73. Fonteneau J-F, Larsson M, Beignon A-S, McKenna K, Dasilva I, Amara A, et al. Human immunodeficiency virus type 1 activates plasmacytoid dendritic cells and concomitantly induces the bystander maturation of myeloid dendritic cells. J Virol (2004) 78:5223–32. doi: 10.1128/jvi.78.10.5223-5232.2004
74. Gringhuis SI, van der Vlist M, van den Berg LM, den Dunnen J, Litjens M, Geijtenbeek TBH. HIV-1 exploits innate signaling by TLR8 and DC-SIGN for productive infection of dendritic cells. Nat Immunol (2010) 11:419–26. doi: 10.1038/ni.1858
75. Gerosa F, Gobbi A, Zorzi P, Burg S, Briere F, Carra G, et al. The reciprocal interaction of NK cells with plasmacytoid or myeloid dendritic cells profoundly affects innate resistance functions. J Immunol (2005) 174:727–34. doi: 10.4049/jimmunol.174.2.727
76. Schuster S, Hurrell B, Tacchini-Cottier F. Crosstalk between neutrophils and dendritic cells: a context-dependent process. J Leukoc Biol (2013) 94:671–5. doi: 10.1189/jlb.1012540
77. Megiovanni AM, Sanchez F, Robledo-Sarmiento M, Morel C, Gluckman JC, Boudaly S. Polymorphonuclear neutrophils deliver activation signals and antigenic molecules to dendritic cells: a new link between leukocytes upstream of T lymphocytes. J Leukoc Biol (2006) 79:977–88. doi: 10.1189/jlb.0905526
78. van Gisbergen KPJM, Sanchez-Hernandez M, Geijtenbeek TBH, van Kooyk Y. Neutrophils mediate immune modulation of dendritic cells through glycosylation-dependent interactions between mac-1 and DC-SIGN. J Exp Med (2005) 201:1281–92. doi: 10.1084/jem.20041276
79. Khorramdelazad H, Kazemi MH, Najafi A, Keykhaee M, Zolfaghari Emameh R, Falak R. Immunopathological similarities between COVID-19 and influenza: Investigating the consequences of Co-infection. Microb Pathog (2021) 152:104554. doi: 10.1016/j.micpath.2020.104554
80. Maltezou HC, Theodoridou K, Poland G. Influenza immunization and COVID-19. Vaccine (2020) 38:6078–9. doi: 10.1016/j.vaccine.2020.07.058
81. Kim E-H, Nguyen T-Q, Casel MAB, Rollon R, Kim S-M, Kim Y-I, et al. Coinfection with Sars-CoV-2 and influenza a virus increases disease severity and impairs neutralizing antibody and CD4+ T cell responses. J Virol (2022) 96:e0187321. doi: 10.1128/jvi.01873-21
82. Dadashi M, Khaleghnejad S, Abedi Elkhichi P, Goudarzi M, Goudarzi H, Taghavi A, et al. COVID-19 and influenza Co-infection: A systematic review and meta-analysis. Front Med (Lausanne) (2021) 8:681469. doi: 10.3389/fmed.2021.681469
83. Swets MC, Russell CD, Harrison EM, Docherty AB, Lone N, Girvan M, et al. Sars-CoV-2 co-infection with influenza viruses, respiratory syncytial virus, or adenoviruses. Lancet (2022) 399(10334):1463–4. doi: 10.1016/S0140-6736(22)00383-X
84. Alosaimi B, Naeem A, Hamed ME, Alkadi HS, Alanazi T, Al Rehily SS, et al. Influenza co-infection associated with severity and mortality in COVID-19 patients. Virol J (2021) 18:127. doi: 10.1186/s12985-021-01594-0
85. Borcherding L, Teksen AS, Grosser B, Schaller T, Hirschbühl K, Claus R, et al. Impaired dendritic cell homing in COVID-19. Front Med (Lausanne) (2021) 8:761372. doi: 10.3389/fmed.2021.761372
86. Winheim E, Rinke L, Lutz K, Reischer A, Leutbecher A, Wolfram L, et al. Impaired function and delayed regeneration of dendritic cells in COVID-19. PloS Pathog (2021) 17:e1009742. doi: 10.1371/journal.ppat.1009742
87. Ghanem MH, Shih AJ, Khalili H, Werth EG, Chakrabarty JK, Brown LM, et al. Proteomic and single-cell transcriptomic dissection of human plasmacytoid dendritic cell response to influenza virus. Front Immunol (2022) 13:814627. doi: 10.3389/fimmu.2022.814627
88. Wimmers F, Donato M, Kuo A, Ashuach T, Gupta S, Li C, et al. The single-cell epigenomic and transcriptional landscape of immunity to influenza vaccination. Cell (2021) 184:3915–3935.e21. doi: 10.1016/j.cell.2021.05.039
89. Onodi F, Bonnet-Madin L, Meertens L, Karpf L, Poirot J, Zhang S-Y, et al. Sars-CoV-2 induces human plasmacytoid predendritic cell diversification via UNC93B and IRAK4. J Exp Med (2021) 218:e20201387. doi: 10.1084/jem.20201387
90. Lee JS, Park S, Jeong HW, Ahn JY, Choi SJ, Lee H, et al. Immunophenotyping of COVID-19 and influenza highlights the role of type I interferons in development of severe COVID-19. Sci Immunol (2020) 5:eabd1554. doi: 10.1126/sciimmunol.abd1554
91. Zhu L, Yang P, Zhao Y, Zhuang Z, Wang Z, Song R, et al. Single-cell sequencing of peripheral mononuclear cells reveals distinct immune response landscapes of COVID-19 and influenza patients. Immunity (2020) 53:685–696.e3. doi: 10.1016/j.immuni.2020.07.009
92. Zolfaghari Emameh R, Nosrati H, Eftekhari M, Falak R, Khoshmirsafa M. Expansion of single cell transcriptomics data of SARS-CoV infection in human bronchial epithelial cells to COVID-19. Biol Proced Online (2020) 22:16. doi: 10.1186/s12575-020-00127-3
93. Ding Q, Lu P, Fan Y, Xia Y, Liu M. The clinical characteristics of pneumonia patients coinfected with 2019 novel coronavirus and influenza virus in wuhan, China. J Med Virol (2020) 92:1549–55. doi: 10.1002/jmv.25781
94. Alculumbre SG, Saint-André V, Di Domizio J, Vargas P, Sirven P, Bost P, et al. Diversification of human plasmacytoid predendritic cells in response to a single stimulus. Nat Immunol (2018) 19:63–75. doi: 10.1038/s41590-017-0012-z
95. Kong Y, Wang Y, Wu X, Han J, Li G, Hua M, et al. Storm of soluble immune checkpoints associated with disease severity of COVID-19. Signal Transduct Target Ther (2020) 5:192. doi: 10.1038/s41392-020-00308-2
96. Vanderbeke L, Van Mol P, Van Herck Y, De Smet F, Humblet-Baron S, Martinod K, et al. Monocyte-driven atypical cytokine storm and aberrant neutrophil activation as key mediators of COVID-19 disease severity. Nat Commun (2021) 12:4117. doi: 10.1038/s41467-021-24360-w
97. Cantuti-Castelvetri L, Ojha R, Pedro LD, Djannatian M, Franz J, Kuivanen S, et al. Neuropilin-1 facilitates Sars-CoV-2 cell entry and infectivity. Science (2020) 370:856–60. doi: 10.1126/science.abd2985
98. Gudowska-Sawczuk M, Mroczko B. The role of neuropilin-1 (NRP-1) in Sars-CoV-2 infection: Review. J Clin Med (2021) 10:2772. doi: 10.3390/jcm10132772
99. Saichi M, Ladjemi MZ, Korniotis S, Rousseau C, Ait Hamou Z, Massenet-Regad L, et al. Single-cell RNA sequencing of blood antigen-presenting cells in severe COVID-19 reveals multi-process defects in antiviral immunity. Nat Cell Biol (2021) 23:538–51. doi: 10.1038/s41556-021-00681-2
100. Alfei F, Ho P-C, Lo W-L. DCision-making in tumors governs T cell anti-tumor immunity. Oncogene (2021) 40:5253–61. doi: 10.1038/s41388-021-01946-8
101. Waldman AD, Fritz JM, Lenardo MJ. A guide to cancer immunotherapy: from T cell basic science to clinical practice. Nat Rev Immunol (2020) 20:651–68. doi: 10.1038/s41577-020-0306-5
102. Kang R, Zhang Q, Zeh HJ, Lotze MT, Tang D. HMGB1 in cancer: good, bad, or both? Clin Cancer Res (2013) 19:4046–57. doi: 10.1158/1078-0432.CCR-13-0495
103. Demoulin S, Herfs M, Somja J, Roncarati P, Delvenne P, Hubert P. HMGB1 secretion during cervical carcinogenesis promotes the acquisition of a tolerogenic functionality by plasmacytoid dendritic cells. Int J Cancer (2015) 137:345–58. doi: 10.1002/ijc.29389
104. Tripathi A, Shrinet K, Kumar A. HMGB1 protein as a novel target for cancer. Toxicol Rep (2019) 6:253–61. doi: 10.1016/j.toxrep.2019.03.002
105. Muire PJ, Schwacha MG, Wenke JC. Systemic T cell exhaustion dynamics is linked to early high mobility group box protein 1 (HMGB1) driven hyper-inflammation in a polytrauma rat model. Cells (2021) 10:1646. doi: 10.3390/cells10071646
106. Zarour HM. Reversing T-cell dysfunction and exhaustion in cancer. Clin Cancer Res (2016) 22:1856–64. doi: 10.1158/1078-0432.CCR-15-1849
107. Marin-Acevedo JA, Kimbrough EO, Lou Y. Next generation of immune checkpoint inhibitors and beyond. J Hematol Oncol (2021) 14:45. doi: 10.1186/s13045-021-01056-8
108. Li B, Chan HL, Chen P. Immune checkpoint inhibitors: basics and challenges. Curr Med Chem (2019) 26:3009–25. doi: 10.2174/0929867324666170804143706
109. Robert C. A decade of immune-checkpoint inhibitors in cancer therapy. Nat Commun (2020) 11:3801. doi: 10.1038/s41467-020-17670-y
110. Messmer D, Yang H, Telusma G, Knoll F, Li J, Messmer B, et al. High mobility group box protein 1: an endogenous signal for dendritic cell maturation and Th1 polarization. J Immunol (2004) 173:307–13. doi: 10.4049/jimmunol.173.1.307
111. Zhang Y, Zheng J. Functions of immune checkpoint molecules beyond immune evasion. Adv Exp Med Biol (2020) 1248:201–26. doi: 10.1007/978-981-15-3266-5_9
112. Iguchi-Manaka A, Kai H, Yamashita Y, Shibata K, Tahara-Hanaoka S, Honda S, et al. Accelerated tumor growth in mice deficient in DNAM-1 receptor. J Exp Med (2008) 205:2959–64. doi: 10.1084/jem.20081611
113. Braun M, Aguilera AR, Sundarrajan A, Corvino D, Stannard K, Krumeich S, et al. CD155 on tumor cells drives resistance to immunotherapy by inducing the degradation of the activating receptor CD226 in CD8+ t cells. Immunity (2020) 53:805–823.e15. doi: 10.1016/j.immuni.2020.09.010
114. Weulersse M, Asrir A, Pichler AC, Lemaitre L, Braun M, Carrié N, et al. Eomes-dependent loss of the Co-activating receptor CD226 restrains CD8+ T cell anti-tumor functions and limits the efficacy of cancer immunotherapy. Immunity (2020) 53:824–839.e10. doi: 10.1016/j.immuni.2020.09.006
115. Simon S, Labarriere N. PD-1 expression on tumor-specific T cells: Friend or foe for immunotherapy? Oncoimmunology (2017) 7:e1364828. doi: 10.1080/2162402X.2017.1364828
116. Shi L, Chen S, Yang L, Li Y. The role of PD-1 and PD-L1 in T-cell immune suppression in patients with hematological malignancies. J Hematol Oncol (2013) 6:74. doi: 10.1186/1756-8722-6-74
117. Bigley V, Barge D, Collin M. Dendritic cell analysis in primary immunodeficiency. Curr Opin Allergy Clin Immunol (2016) 16:530–40. doi: 10.1097/ACI.0000000000000322
118. Bandola-Simon J, Roche PA. Dysfunction of antigen processing and presentation by dendritic cells in cancer. Mol Immunol (2019) 113:31–7. doi: 10.1016/j.molimm.2018.03.025
119. Bieber K, Autenrieth SE. Dendritic cell development in infection. Mol Immunol (2020) 121:111–7. doi: 10.1016/j.molimm.2020.02.015
120. Bates J, Diehl L. Dendritic cells in IBD pathogenesis: an area of therapeutic opportunity? J Pathol (2014) 232:112–20. doi: 10.1002/path.4277
121. Sasindran SJ, Torrelles JB. Mycobacterium tuberculosis infection and inflammation: what is beneficial for the host and for the bacterium? Front Microbiol (2011) 2:2. doi: 10.3389/fmicb.2011.00002
122. Helikar T, Kochi N, Kowal B, Dimri M, Naramura M, Raja SM, et al. A comprehensive, multi-scale dynamical model of ErbB receptor signal transduction in human mammary epithelial cells. PloS One (2013) 8:e61757. doi: 10.1371/journal.pone.0061757
123. Hernandez C, Thomas-Chollier M, Naldi A, Thieffry D. Computational verification of Large logical models-application to the prediction of T cell response to checkpoint inhibitors. Front Physiol (2020) 11:558606. doi: 10.3389/fphys.2020.558606
124. Saez-Rodriguez J, Simeoni L, Lindquist JA, Hemenway R, Bommhardt U, Arndt B, et al. A logical model provides insights into T cell receptor signaling. PloS Comput Biol (2007) 3:e163. doi: 10.1371/journal.pcbi.0030163
125. Montagud A, Béal J, Tobalina L, Traynard P, Subramanian V, Szalai B, et al. Patient-specific boolean models of signalling networks guide personalised treatments. eLife (2022) 11:e72626. doi: 10.7554/eLife.72626
126. Wynn ML, Consul N, Merajver SD, Schnell S. Logic-based models in systems biology: a predictive and parameter-free network analysis method. Integr Biol (Camb) (2012) 4:1323–37. doi: 10.1039/c2ib20193c
127. Yang H, Wang H, Andersson U. Targeting inflammation driven by HMGB1. Front Immunol (2020) 11:484. doi: 10.3389/fimmu.2020.00484
128. Lee S-A, Kwak MS, Kim S, Shin J-S. The role of high mobility group box 1 in innate immunity. Yonsei Med J (2014) 55:1165–76. doi: 10.3349/ymj.2014.55.5.1165
129. Yu M, Wang H, Ding A, Golenbock DT, Latz E, Czura CJ, et al. HMGB1 signals through toll-like receptor (TLR) 4 and TLR2. Shock (2006) 26:174–9. doi: 10.1097/01.shk.0000225404.51320.82
130. Mougel A, Terme M, Tanchot C. Therapeutic cancer vaccine and combinations with antiangiogenic therapies and immune checkpoint blockade. Front Immunol (2019) 10:467. doi: 10.3389/fimmu.2019.00467
131. Zhao J, Chen Y, Ding Z-Y, Liu J-Y. Safety and efficacy of therapeutic cancer vaccines alone or in combination with immune checkpoint inhibitors in cancer treatment. Front Pharmacol (2019) 10:1184. doi: 10.3389/fphar.2019.01184
132. Wang D, Huang XF, Hong B, Song X-T, Hu L, Jiang M, et al. Efficacy of intracellular immune checkpoint-silenced DC vaccine. JCI Insight (2018) 3:e98368. doi: 10.1172/jci.insight.98368
133. Kong BY, Bolton H, Kim JW, Silveira PA, Fromm PD, Clark GJ. On the other side: manipulating the immune checkpoint landscape of dendritic cells to enhance cancer immunotherapy. Front Oncol (2019) 9:50. doi: 10.3389/fonc.2019.00050
134. Jin H-S, Park Y. Hitting the complexity of the TIGIT-CD96-CD112R-CD226 axis for next-generation cancer immunotherapy. BMB Rep (2021) 54:2–11. doi: 10.5483/BMBRep.2021.54.1.229
135. Moon J, Oh YM, Ha S-J. Perspectives on immune checkpoint ligands: expression, regulation, and clinical implications. BMB Rep (2021) 54:403–12. doi: 10.5483/BMBRep.2021.54.8.054
136. Clark M, Kroger CJ, Tisch RM. Type 1 diabetes: A chronic anti-Self-Inflammatory response. Front Immunol (2017) 8:1898. doi: 10.3389/fimmu.2017.01898
137. Mollah ZUA, Pai S, Moore C, O’Sullivan BJ, Harrison MJ, Peng J, et al. Abnormal NF-kappa b function characterizes human type 1 diabetes dendritic cells and monocytes. J Immunol (2008) 180:3166–75. doi: 10.4049/jimmunol.180.5.3166
138. Babon JAB, DeNicola ME, Blodgett DM, Crèvecoeur I, Buttrick TS, Maehr R, et al. Analysis of self-antigen specificity of islet-infiltrating T cells from human donors with type 1 diabetes. Nat Med (2016) 22:1482–7. doi: 10.1038/nm.4203
139. Grohová A, Dáňová K, Špíšek R, Palová-Jelínková L. Cell based therapy for type 1 diabetes: should we take hyperglycemia into account? Front Immunol (2019) 10:79. doi: 10.3389/fimmu.2019.00079
Keywords: systems immunology, predictive modeling, In-Silico experiments, antigen-presenting cell, dendritic cell, immunology highlights
Citation: Aghamiri SS, Puniya BL, Amin R and Helikar T (2023) A multiscale mechanistic model of human dendritic cells for in-silico investigation of immune responses and novel therapeutics discovery. Front. Immunol. 14:1112985. doi: 10.3389/fimmu.2023.1112985
Received: 30 November 2022; Accepted: 22 February 2023;
Published: 10 March 2023.
Edited by:
Mohadeseh Zarei Ghobadi, University of Isfahan, IranReviewed by:
Winston Garira, University of Venda, South AfricaMaximilien Grandclaudon, Dana-Farber Cancer Institute, United States
Abazar Arab Ameri, University of Zanjan, Iran
Copyright © 2023 Aghamiri, Puniya, Amin and Helikar. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Rada Amin, rada.aminali@gmail.com; Tomáš Helikar, thelikar2@unl.edu