Search

Single-cell RNA sequencing (scRNA-seq) data has been widely used for cell trajectory inference, with the assumption that cells with similar expression profiles share the same differentiation state. However, the inferred trajectory may not reveal clonal differentiation heterogeneity among T cell clones. Single-cell T cell receptor sequencing (scTCR-seq) data provides invaluable insights into the clonal relationship among cells, yet it lacks functional characteristics. Therefore, scRNA-seq and scTCR-seq data complement each other in improving trajectory inference, where a reliable computational tool is still missing. We developed LRT, a computational framework for the integrative analysis of scTCR-seq and scRNA-seq data to explore clonal differentiation trajectory heterogeneity. Specifically, LRT uses the transcriptomics information from scRNA-seq data to construct overall cell trajectories and then utilizes both the TCR sequence information and phenotype information to identify clonotype clusters with distinct differentiation biasedness. LRT provides a comprehensive analysis workflow, including preprocessing, cell trajectory inference, clonotype clustering, trajectory biasedness evaluation, and clonotype cluster characterization. We illustrated its utility using scRNA-seq and scTCR-seq data of CD8+ T cells and CD4+ T cells with acute lymphocytic choriomeningitis virus infection. These analyses identified several clonotype clusters with distinct skewed distribution along the differentiation path, which cannot be revealed solely based on scRNA-seq data. Clones from different clonotype clusters exhibited diverse expansion capability, V-J gene usage pattern and CDR3 motifs. The LRT framework was implemented as an R package 'LRT', and it is now publicly accessible at https://github.com/JuanXie19/LRT. In addition, it provides two Shiny apps 'shinyClone' and 'shinyClust' that allow users to interactively explore distributions of clonotypes, conduct repertoire analysis, implement clustering of clonotypes, trajectory biasedness evaluation, and clonotype cluster characterization. Author summary: Understanding the dynamic changes behind biological processes is important for determining molecular mechanisms underlying normal tissue formulation, developmental disorders and pathologies. Usually, a biological process can be characterized by identifying a trajectory, a path that goes through the various cellular states associated with the process. Since cells in different states may express different sets of genes, researchers often infer cell trajectory via capturing transcriptomics changes. Dozens of methods have been developed for cell trajectory inference, and scRNA-seq data is predominantly utilized. However, cells from different clones may exist distinct differentiation patterns, while methods based only on scRNA-seq data cannot capture this heterogeneity. T cells play a key role in the immune system, and their high antigen recognition specificity is largely determined by their TCR sequences. Thanks to the advent of scTCR-seq technology, we can now identify the group of cells coming from the same clone. This paper describes our novel computational framework, namely LRT, and demonstrates that by complementing scRNA-seq data with the clonal information from scTCR-seq data using LRT, we are able to examine the clonal differentiation heterogeneity that cannot be revealed solely based on scRNA-seq data. [ABSTRACT FROM AUTHOR]

Tuberculosis (TB) remains a global public health threat. Understanding the dynamics of host-pathogen interactions within TB granulomas will assist in identifying what leads to the successful elimination of infection. In vitro TB models provide a controllable environment to study these granuloma dynamics. Previously we developed a biomimetic 3D spheroid granuloma model that controls bacteria better than a traditional monolayer culture counterpart. We used agent-based simulations to predict the mechanistic reason for this difference. Our calibrated simulations were able to predict heterogeneous bacterial dynamics that are consistent with experimental data. In one group of simulations, spheroids are found to have higher macrophage activation than their traditional counterparts, leading to better bacterial control. This higher macrophage activation in the spheroids was not due to higher counts of activated T cells, instead fewer activated T cells were able to activate more macrophages due to the proximity of these cells to each other within the spheroid. In a second group of simulations, spheroids again have more macrophage activation but also more T cell activation, specifically CD8+ T cells. This higher level of CD8+ T cell activation is predicted to be due to the proximity of these cells to the cells that activate them. Multiple mechanisms of control were predicted. Simulations removing individual mechanisms show that one group of simulations has a CD4+ T cell dominant response, while the other has a mixed/CD8+ T cell dominant response. Lastly, we demonstrated that in spheroids the initial structure and movement rules work synergistically to reduce bacterial load. These findings provide valuable insights into how the structural complexity of in vitro models impacts immune responses. Moreover, our study has implications for engineering more physiologically relevant in vitro models and advancing our understanding of TB pathogenesis and potential therapeutic interventions. Author summary: Tuberculosis (TB) is a respiratory infection that is associated with granulomas, organized clusters of immune cells that form around the bacteria. Recent work has focused on making in vitro models that better approximate these granulomas, such as our spheroid granuloma model. We found that this spheroid was better able to control bacteria than its traditional culture counterpart, so, in this work, we use agent-based modeling to explore the causes of this differential control. The simulated spheroids have increased macrophage activation. Surprisingly, this is not due to larger numbers of activated T cells, which activate these macrophages. Rather, the increase in activated macrophages is due to their proximity to fewer activated T cells. We also see a subset of simulations where the spheroids have increased cytotoxic CD8+ T cell activation. The simulation is further used to simulate changes in the structure of the in vitro models to explore its impact on immune response and bacterial control. These findings reveal the importance of structural complexity in in vitro models and offer insights into TB pathogenesis. Ultimately, our work paves the way for more realistic in vitro models and a deeper understanding of how our bodies combat TB. [ABSTRACT FROM AUTHOR]

The advancements in next-generation sequencing have made it possible to effectively detect somatic mutations, which has led to the development of personalized neoantigen cancer vaccines that are tailored to the unique variants found in a patient's cancer. These vaccines can provide significant clinical benefit by leveraging the patient's immune response to eliminate malignant cells. However, determining the optimal vaccine dose for each patient is a challenge due to the heterogeneity of tumors. To address this challenge, we formulate a mathematical dose optimization problem based on a previous mathematical model that encompasses the immune response cascade produced by the vaccine in a patient. We propose an optimization approach to identify the optimal personalized vaccine doses, considering a fixed vaccination schedule, while simultaneously minimizing the overall number of tumor and activated T cells. To validate our approach, we perform in silico experiments on six real-world clinical trial patients with advanced melanoma. We compare the results of applying an optimal vaccine dose to those of a suboptimal dose (the dose used in the clinical trial and its deviations). Our simulations reveal that an optimal vaccine regimen of higher initial doses and lower final doses may lead to a reduction in tumor size for certain patients. Our mathematical dose optimization offers a promising approach to determining an optimal vaccine dose for each patient and improving clinical outcomes. Author summary: The development of neoantigen cancer vaccines have rapidly increased over the past decade with the advancement of next-generation sequencing technologies to determine immunogenic peptides from patient's somatic mutations. However, traditional methods to determine the cancer vaccine dose often produce suboptimal clinical outcomes. This work use our previous mathematical model to represent the immunological cascade at the cellular and subcellular levels elicited by the vaccine dose, and focuses on developing a mathematical optimization approach to identify the optimal vaccine dose to minimize two objective functionals, (i) minimize the amount of peptide dose and tumor size, and (ii) minimize the number of activated T cells in addition to the objective functional i. The optimization approach allows the identification of optimal vaccine doses among a set of tested doses for a higher clinical benefit in tumor reduction. To demonstrate and validate our optimization approach, we perform in silico experiments on six patients with melanoma from a clinical trial study. The results show that the predicted optimal vaccine doses can provide higher clinical benefit in tumor reduction when compared to the clinical trial doses for some patients. [ABSTRACT FROM AUTHOR]

T cells rely on their T cell receptors (TCRs) to discern foreign antigens presented by human leukocyte antigen (HLA) proteins. The TCRs of an individual contain a record of this individual's past immune activities, such as immune response to infections or vaccines. Mining the TCR data may recover useful information or biomarkers for immune related diseases or conditions. Some TCRs are observed only in the individuals with certain HLA alleles, and thus characterizing TCRs requires a thorough understanding of TCR-HLA associations. The extensive diversity of HLA alleles and the rareness of some HLA alleles present a formidable challenge for this task. Existing methods either treat HLA as a categorical variable or represent an HLA by its alphanumeric name, and have limited ability to generalize to the HLAs that are not seen in the training process. To address this challenge, we propose a neural network-based method named Deep learning Prediction of TCR-HLA association (DePTH) to predict TCR-HLA associations based on their amino acid sequences. We demonstrate that DePTH is capable of making reasonable predictions for TCR-HLA associations, even when neither the HLA nor the TCR have been included in the training dataset. Furthermore, we establish that DePTH can be used to quantify the functional similarities among HLA alleles, and that these HLA similarities are associated with the survival outcomes of cancer patients who received immune checkpoint blockade treatments. Author summary: T cells are critical components of the human immune system. A T cell receptor (TCR) is a protein complex found on the surface of T cells, and it determines the antigens that a T cell recognizes. The TCR repertoire, which is the collection of all the TCRs within an individual, provides rich information about the history and current activities of the immune system. A TCR-antigen interaction is mediated by a protein called human leukocyte antigen (HLA). HLA genes are highly polymorphic in the human population, contributing to the variation of human immune response to different types of antigens. Studying the associations between TCRs and HLAs can help us identify functional TCRs given the HLAs of an individual. To this end, we develop a deep learning method to predict TCR-HLA associations based on their amino acid sequences. This method allows us to borrow information across different HLA genes. We demonstrate that the predictions of our model can be used to quantify the functional similarities of HLA alleles and such similarities are associated with cancer patients' survival outcome. [ABSTRACT FROM AUTHOR]

T cells form transient cell-to-cell contacts with antigen presenting cells (APCs) to facilitate surface interrogation by membrane bound T cell receptors (TCRs). Upon recognition of molecular signatures (antigen) of pathogen, T cells may initiate an adaptive immune response. The duration of the T cell/APC contact is observed to vary widely, yet it is unclear what constructive role, if any, such variations might play in immune signaling. Modeling efforts describing antigen discrimination often focus on steady-state approximations and do not account for the transient nature of cellular contacts. Within the framework of a kinetic proofreading (KP) mechanism, we develop a stochastic First Receptor Activation Model (FRAM) describing the likelihood that a productive immune signal is produced before the expiry of the contact. Through the use of extreme statistics, we characterize the probability that the first TCR triggering is induced by a rare agonist antigen and not by that of an abundant self-antigen. We show that defining positive immune outcomes as resilience to extreme statistics and sensitivity to rare events mitigates classic tradeoffs associated with KP. By choosing a sufficient number of KP steps, our model is able to yield single agonist sensitivity whilst remaining non-reactive to large populations of self antigen, even when self and agonist antigen are similar in dissociation rate to the TCR but differ largely in expression. Additionally, our model achieves high levels of accuracy even when agonist positive APCs encounters are rare. Finally, we discuss potential biological costs associated with high classification accuracy, particularly in challenging T cell environments. Author summary: Physical contact between the T cell and antigen presenting cell (APC) is essential for productive immune signaling. Wide variations in this contact time have been observed yet little is known of mechanisms controlling this crucial timescale, nor how its duration may impact antigen discrimination. We develop and analyze a probabilistic mathematical model of T cell activation which combines kinetic proofreading (KP) with a finite contact duration. Our model is capable of suppressing large populations of self ligands while remaining sensitive to only a single agonist in T cell/APC cellular contacts. Additionally, we explored two challenging cases, one in which self and agonist antigen are similar and one in which agonist positive APCs are rare. We found that our model could overcome these environmental challenges by increasing the number of kinetic proofreading steps. Finally, we discuss the potential biological costs of achieving such accuracy. Our work demonstrates the extreme effectiveness of kinetic proofreading in a temporal context while also demonstrating the possible challenges in biological implementation of such a model. [ABSTRACT FROM AUTHOR]

Random walks on networks are widely used to model stochastic processes such as search strategies, transportation problems or disease propagation. A prominent example of such process is the dynamics of naive T cells within the lymph node while they are scanning for antigens. The observed T cells trajectories in small sub-volumes of the lymph node are well modeled as a random walk and they have been shown to follow the lymphatic conduit network as substrate for migration. One can then ask how does the connectivity patterns of the lymph node conduit network affect the T cells collective exploration behavior. In particular, does the network display properties that are uniform across the whole volume of the lymph node or can we distinguish some heterogeneities? We propose a workflow to accurately and efficiently define and compute these quantities on large networks, which enables us to characterize heterogeneities within a very large published dataset of Lymph Node Conduit Network. To establish the significance of our results, we compared the results obtained on the lymph node to null models of varying complexity. We identified significantly heterogeneous regions characterized as "remote regions" at the poles and next to the medulla, while a large portion of the network promotes uniform exploration by T cells. Author summary: Lymph nodes are organs in which actors of the immune system meet. In particular, the encounter between the naive T cells and their specific antigens occurs in lymph nodes. This event triggers the adaptive immune response. T cells movement has been shown to be well described as a random walk, at least when they are measured on small sub-volumes of the lymph node. In parallel, it was shown that T-cells migrate following the lymphatic conduit network that span the lymph nodes. In this study, we ask, how does the connectivity pattern of the conduit network on which T cells move influences their collective exploration behavior? Are there regions in the lymph node conduit network which have distinct random walk related properties? The topological reconstruction of the lymph node conduit network was recently made available. The network is very large (about 200 000 nodes) and appears very regular, with most nodes being connected to three neighbours. We propose a workflow to detect heterogeneities in such large and quasi-regular networks, building on random walk on network tools, and the measure of two features which we interpret using a series of generated null models for comparison. We show that the lymph node conduit network displays remotely accessible regions at both poles and near medulla, with however most of the network promoting uniform exploration. [ABSTRACT FROM AUTHOR]

Data clustering plays a significant role in biomedical sciences, particularly in single-cell data analysis. Researchers use clustering algorithms to group individual cells into populations that can be evaluated across different levels of disease progression, drug response, and other clinical statuses. In many cases, multiple sets of clusters must be generated to assess varying levels of cluster specificity. For example, there are many subtypes of leukocytes (e.g. T cells), whose individual preponderance and phenotype must be assessed for statistical/functional significance. In this report, we introduce a novel hierarchical density clustering algorithm (HAL-x) that uses supervised linkage methods to build a cluster hierarchy on raw single-cell data. With this new approach, HAL-x can quickly predict multiple sets of labels for immense datasets, achieving a considerable improvement in computational efficiency on large datasets compared to existing methods. We also show that cell clusters generated by HAL-x yield near-perfect F1-scores when classifying different clinical statuses based on single-cell profiles. Our hierarchical density clustering algorithm achieves high accuracy in single cell classification in a scalable, tunable and rapid manner. Author summary: Modern experimental techniques such as mass cytometry (CyTOF) make it possible to quickly make high-dimensional measurements on upwards of tens of millions of cells with single-cell resolution. An important problem in biology is to use these measurements to group together similar cells to identify biologically meaningful cell types that can be used to study disease progression, drug responses, and other clinical outcomes. However, the size and complexity of experimental data sets makes this problem computationally and theoretically extremely difficult. Here, we present a new algorithm HAL-X that accurately and quickly identifies cell clusters from biological data. Importantly, our algorithm does not require large amounts of memory. This eliminates the need for specialized high-end computing resources, allowing biologists to quickly analyze their data using a standard laptop or desktop computer. [ABSTRACT FROM AUTHOR]

The human immune system consists of a highly intelligent network of billions of independent, self-organized cells that interact with each other. Machine learning (ML) is an artificial intelligence (AI) tool that automatically processes huge amounts of image data. Immunotherapies have revolutionized the treatment of blood cancer. Specifically, one such therapy involves engineering immune cells to express chimeric antigen receptors (CAR), which combine tumor antigen specificity with immune cell activation in a single receptor. To improve their efficacy and expand their applicability to solid tumors, scientists optimize different CARs with different modifications. However, predicting and ranking the efficacy of different "off-the-shelf" immune products (e.g., CAR or Bispecific T-cell Engager [BiTE]) and selection of clinical responders are challenging in clinical practice. Meanwhile, identifying the optimal CAR construct for a researcher to further develop a potential clinical application is limited by the current, time-consuming, costly, and labor-intensive conventional tools used to evaluate efficacy. Particularly, more than 30 years of immunological synapse (IS) research data demonstrate that T cell efficacy is not only controlled by the specificity and avidity of the tumor antigen and T cell interaction, but also it depends on a collective process, involving multiple adhesion and regulatory molecules, as well as tumor microenvironment, spatially and temporally organized at the IS formed by cytotoxic T lymphocytes (CTL) and natural killer (NK) cells. The optimal function of cytotoxic lymphocytes (including CTL and NK) depends on IS quality. Recognizing the inadequacy of conventional tools and the importance of IS in immune cell functions, we investigate a new strategy for assessing CAR-T efficacy by quantifying CAR IS quality using the glass-support planar lipid bilayer system combined with ML-based data analysis. Previous studies in our group show that CAR-T IS quality correlates with antitumor activities in vitro and in vivo. However, current manually quantified IS quality data analysis is time-consuming and labor-intensive with low accuracy, reproducibility, and repeatability. In this study, we develop a novel ML-based method to quantify thousands of CAR cell IS images with enhanced accuracy and speed. Specifically, we used artificial neural networks (ANN) to incorporate object detection into segmentation. The proposed ANN model extracts the most useful information to differentiate different IS datasets. The network output is flexible and produces bounding boxes, instance segmentation, contour outlines (borders), intensities of the borders, and segmentations without borders. Based on requirements, one or a combination of this information is used in statistical analysis. The ML-based automated algorithm quantified CAR-T IS data correlates with the clinical responder and non-responder treated with Kappa-CAR-T cells directly from patients. The results suggest that CAR cell IS quality can be used as a potential composite biomarker and correlates with antitumor activities in patients, which is sufficiently discriminative to further test the CAR IS quality as a clinical biomarker to predict response to CAR immunotherapy in cancer. For translational research, the method developed here can also provide guidelines for designing and optimizing numerous CAR constructs for potential clinical development. Trial Registration: ClinicalTrials.gov NCT00881920. Author summary: Adoptive transfer of chimeric antigen receptor (CAR)-modified immune cells (including CAR-T and CAR-NK cells) have revolutionized the treatment of cancer with success in clinical trials treating multiple myeloma, leukemia, sarcoma, and neuroblastoma. However, CAR-modified immune cells (particularly CAR-T cells) must form a functional immunological synapse (IS) with their susceptible tumor cells to be effective in clinics. Currently, there are no effective biomarkers to predict CAR efficacy in vivo. In this study, we develop a state-of-the-art machine learning (ML) detection, and segmentation method to measure the quality of the CAR-T cell IS using CAR-T samples from patients. We automate the IS quality analysis to develop effective prognostic applications of CAR-T therapies for cancer patients. The fast, easy-to-implement Synapse Predicts Efficacy (SPE) assay we propose will streamline CAR development and selection, ultimately optimizing clinical outcome(s) for patients undergoing these rapidly evolving immunotherapies. This technology can lead to development of fast and easy tools to predict CAR-T cell effectiveness in cancer patients. [ABSTRACT FROM AUTHOR]

A key question in SARS-CoV-2 infection is why viral loads and patient outcomes vary dramatically across individuals. Because spatial-temporal dynamics of viral spread and immune response are challenging to study in vivo, we developed Spatial Immune Model of Coronavirus (SIMCoV), a scalable computational model that simulates hundreds of millions of lung cells, including respiratory epithelial cells and T cells. SIMCoV replicates viral growth dynamics observed in patients and shows how spatially dispersed infections can lead to increased viral loads. The model also shows how the timing and strength of the T cell response can affect viral persistence, oscillations, and control. By incorporating spatial interactions, SIMCoV provides a parsimonious explanation for the dramatically different viral load trajectories among patients by varying only the number of initial sites of infection and the magnitude and timing of the T cell immune response. When the branching airway structure of the lung is explicitly represented, we find that virus spreads faster than in a 2D layer of epithelial cells, but much more slowly than in an undifferentiated 3D grid or in a well-mixed differential equation model. These results illustrate how realistic, spatially explicit computational models can improve understanding of within-host dynamics of SARS-CoV-2 infection. Author summary: A key question in SARS-CoV-2 infection is why viral loads and patient outcomes are so different across individuals. Because it's difficult to see how the virus spreads in the lungs of infected people, we developed Spatial Immune Model of Coronavirus (SIMCoV), a computational model that simulates hundreds of millions of cells, including lung cells and immune cells. SIMCoV simulates how virus grows and then declines, and the simulations match data observed in patients. SIMCoV shows that when there are more initial infection sites, the virus grows to a higher peak. The model also shows how the timing of the immune response, particularly the T cell response, can affect how long the virus persists and whether it is ultimately cleared from the lungs. SIMCoV shows that the different viral loads in different patients can be explained by how many different places the virus is initially seeded inside their lungs. We explicitly add the branching airway structure of the lung into the model and show that virus spreads slightly faster than it would in a 2D layer of lung cells, but much slower than in traditional mathematical models based on differential equations. These results illustrate how realistic spatial computational models can improve understanding of how SARS-CoV-2 infection spreads in the lung. [ABSTRACT FROM AUTHOR]

MHC-I associated peptides (MAPs) play a central role in the elimination of virus-infected and neoplastic cells by CD8 T cells. However, accurately predicting the MAP repertoire remains difficult, because only a fraction of the transcriptome generates MAPs. In this study, we investigated whether codon arrangement (usage and placement) regulates MAP biogenesis. We developed an artificial neural network called Codon Arrangement MAP Predictor (CAMAP), predicting MAP presentation solely from mRNA sequences flanking the MAP-coding codons (MCCs), while excluding the MCC per se. CAMAP predictions were significantly more accurate when using original codon sequences than shuffled codon sequences which reflect amino acid usage. Furthermore, predictions were independent of mRNA expression and MAP binding affinity to MHC-I molecules and applied to several cell types and species. Combining MAP ligand scores, transcript expression level and CAMAP scores was particularly useful to increase MAP prediction accuracy. Using an in vitro assay, we showed that varying the synonymous codons in the regions flanking the MCCs (without changing the amino acid sequence) resulted in significant modulation of MAP presentation at the cell surface. Taken together, our results demonstrate the role of codon arrangement in the regulation of MAP presentation and support integration of both translational and post-translational events in predictive algorithms to ameliorate modeling of the immunopeptidome. Author summary: MHC-I associated peptides (MAPs) are small fragments of intracellular proteins presented at the surface of cells and used by the immune system to detect and eliminate cancerous or virus-infected cells. While it is theoretically possible to predict which portions of the intracellular proteins will be naturally processed by the cells to ultimately reach the surface, current methodologies have prohibitively high false discovery rates. Here we introduce an artificial neural network called Codon Arrangement MAP Predictor (CAMAP) which integrates information from mRNA-to-protein translation to other factors regulating MAP biogenesis (e.g. MAP ligand score and transcript expression levels) to improve MAP prediction accuracy. While most MAP predictive approaches focus on MAP sequences per se, CAMAP's novelty is to analyze the MAP-flanking mRNA sequences, thereby providing completely independent information for MAP prediction. We show on several datasets that the integration of CAMAP scores with other known factors involved in MAP presentation (i.e. MAP ligand score and mRNA expression) significantly improves MAP prediction accuracy, and further validate CAMAP learned features using an in-vitro assay. These findings may have major implications for the design of vaccines against cancers and viruses, and in times of pandemics could accelerate the identification of relevant MAPs of viral origins. [ABSTRACT FROM AUTHOR]

Understanding how immunological memory lasts a lifetime requires quantifying changes in the number of memory cells as well as how their division and death rates change over time. We address these questions by using a statistically powerful mixed-effects differential equations framework to analyze data from two human studies that follow CD8 T cell responses to the yellow fever vaccine (YFV-17D). Models were first fit to the frequency of YFV-specific memory CD8 T cells and deuterium enrichment in those cells 42 days to 1 year post-vaccination. A different dataset, on the loss of YFV-specific CD8 T cells over three decades, was used to assess out of sample predictions of our models. The commonly used exponential and bi-exponential decline models performed relatively poorly. Models with the cell loss following a power law (exactly or approximately) were most predictive. Notably, using only the first year of data, these models accurately predicted T cell frequencies up to 30 years post-vaccination. Our analyses suggest that division rates of these cells drop and plateau at a low level (0.1% per day, ∼ double the estimated values for naive T cells) within one year following vaccination, whereas death rates continue to decline for much longer. Our results show that power laws can be predictive for T cell memory, a finding that may be useful for vaccine evaluation and epidemiological modeling. Moreover, since power laws asymptotically decline more slowly than any exponential decline, our results help explain the longevity of immune memory phenomenologically. Author summary: Immunological memory, generated in response to infection or vaccination, may provide complete or partial protection from antigenically similar infections for the lifetime. Memory CD8 T cells are important players in protection from secondary viral infections but quantitative understanding of their dynamics in humans is limited. We analyze data from two studies where immunization with the yellow fever virus vaccine (YFV-17D) generates a mild acute infection and long-term memory. We find that: (i) the division rate of YFV-17D-specific CD8 T cells drops and stabilizes at ∼ 0.1% per day during the first year following vaccination whereas the death rate declines more gradually, and (ii) the number of these cells declines approximately in accordance with a power law (∝ time−0.75) for at least several decades following vaccination. [ABSTRACT FROM AUTHOR]

Immune responses rely on a complex adaptive system in which the body and infections interact at multiple scales and in different compartments. We developed a modular model of CD4+ T cells, which uses four modeling approaches to integrate processes at three spatial scales in different tissues. In each cell, signal transduction and gene regulation are described by a logical model, metabolism by constraint-based models. Cell population dynamics are described by an agent-based model and systemic cytokine concentrations by ordinary differential equations. A Monte Carlo simulation algorithm allows information to flow efficiently between the four modules by separating the time scales. Such modularity improves computational performance and versatility and facilitates data integration. We validated our technology by reproducing known experimental results, including differentiation patterns of CD4+ T cells triggered by different combinations of cytokines, metabolic regulation by IL2 in these cells, and their response to influenza infection. In doing so, we added multi-scale insights to single-scale studies and demonstrated its predictive power by discovering switch-like and oscillatory behaviors of CD4+ T cells that arise from nonlinear dynamics interwoven across three scales. We identified the inflamed lymph node's ability to retain naive CD4+ T cells as a key mechanism in generating these emergent behaviors. We envision our model and the generic framework encompassing it to serve as a tool for understanding cellular and molecular immunological problems through the lens of systems immunology. Author summary: CD4+ T cells are a key part of the adaptive immune system. They differentiate into different phenotypes to carry out different functions. They do so by secreting molecules called cytokines to regulate other immune cells. Multi-scale modeling can potentially explain their emergent behaviors by integrating biological phenomena occurring at different spatial (intracellular, cellular, and systemic), temporal, and organizational scales (signal transduction, gene regulation, metabolism, cellular behaviors, and cytokine transport). We built a computational platform by combining disparate modeling frameworks (compartmental ordinary differential equations, agent-based modeling, Boolean network modeling, and constraint-based modeling). We validated the platform's ability to predict CD4+ T cells' emergent behaviors by reproducing their differentiation patterns, metabolic regulation, and population dynamics in response to influenza infection. We then used it to predict and explain novel switch-like and oscillatory behaviors for CD4+ T cells. On the basis of these results, we believe that our multi-approach and multi-scale platform will be a valuable addition to the systems immunology toolkit. In addition to its immediate relevance to CD4+ T cells, it also has the potential to become the foundation of a virtual immune system. [ABSTRACT FROM AUTHOR]