- Research
- Open Access
Hybrid approach to model the spatial regulation of T cell responses
- Anass Bouchnita^{1, 2, 3},
- Gennady Bocharov^{4}Email author,
- Andreas Meyerhans^{4, 5, 6} and
- Vitaly Volpert^{1, 4, 7}
https://doi.org/10.1186/s12865-017-0205-0
© Bouchnita et al. 2017
- Published: 21 June 2017
Abstract
Background
Moving from the molecular and cellular level to a multi-scale systems understanding of immune responses requires the development of novel approaches to integrate knowledge and data from different biological levels into mechanism-based integrative mathematical models. The aim of our study is to present a methodology for a hybrid modelling of immunological processes in their spatial context.
Methods
A two-level hybrid mathematical model of immune cell migration and interaction integrating cellular and organ levels of regulation for a 2D spatial consideration of idealized secondary lymphoid organs is developed. It considers the population dynamics of antigen-presenting cells, CD4 ^{+} and CD8 ^{+} T lymphocytes in naive-, proliferation- and differentiated states. Cell division is assumed to be asymmetric and regulated by the extracellular concentration of interleukin-2 (IL-2) and type I interferon (IFN), together controlling the balance between proliferation and differentiation. The cytokine dynamics is described by reaction-diffusion PDEs whereas the intracellular regulation is modelled with a system of ODEs.
Results
The mathematical model has been developed, calibrated and numerically implemented to study various scenarios in the regulation of T cell immune responses to infection, in particular the change in the diffusion coefficient of type I IFN as compared to IL-2. We have shown that a hybrid modelling approach provides an efficient tool to describe and analyze the interplay between spatio-temporal processes in the emergence of abnormal immune response dynamics.
Discussion
Virus persistence in humans is often associated with an exhaustion of T lymphocytes. Many factors can contribute to the development of exhaustion. One of them is associated with a shift from a normal clonal expansion pathway to an altered one characterized by an early terminal differentiation of T cells. We propose that an altered T cell differentiation and proliferation sequence can naturally result from a spatial separation of the signaling events delivered via TCR, IL-2 and type I IFN receptors. Indeed, the spatial overlap of the concentration fields of extracellular IL-2 and IFN in lymph nodes changes dynamically due to different migration patterns of APCs and CD4 ^{+} T cells secreting them.
Conclusions
The proposed hybrid mathematical model of the immune response represents a novel analytical tool to examine challenging issues in the spatio-temporal regulation of cell growth and differentiation, in particular the effect of timing and location of activation signals.
Keywords
- Immune system
- T cell
- Spatial dynamics
- Multi-scale regulation
- Hybrid model
Background
The immune system is regulated by multiple processes at various levels of biological organization including the genetic-, cellular-, tissue-, organ- and the whole organism levels. The resulting structural and functional complexity of the immune system called for a major shift towards information-rich, systems-based approaches in immunological research. High throughput technologies generate vast amounts of data that facilitate dissection of the immunological processes at ever finer resolution. The need to embed immune processes into their spatial context both at the molecular- and cellular level is a hallmark of the systems immunology approach. In fact, there are many examples of how the fate decisions in the immune system depend on the spatial-temporal dynamics of cytokines, e.g. the interleukin-2 (IL-2) [1] and type I interferon (IFN) [2, 3].
Overview of the hybrid and multiscale approaches to model the spatial dynamics of immune responses
Model | Phenomena | Process considered | Types of equations | State variables |
---|---|---|---|---|
Baldazzi et al. [7] | Immune response to antigen in lymph node (500 hrs) | Clonal expansion, 3D: transport, reaction-diffusion | Agent-based for cells, PDEs for molecules antigen, chemokines | DCs, B-cells, CD4 ^{+} T cells, |
Fallahi-Sichani et al. [8] | Immune response in Tuberculosis, Granuloma (200 days) | Clonal expansion, 2D: chemotaxis, cell-to-cell interactions single-cell state regulation | Agent-based for cells, ODEs for cytokines, 2D geometry of lung tissue | Macrophages, CD8 ^{+} T cells, Treg cells, T γ cells, M. tuberculosis, TNF α, TNFR |
Gong et al. [9] | Immune response to antigen in lymph node (550 hrs) | Clonal expansion, 3D: trafficking, cell-to-cell interactions | Agent-based for cells, anatomicaly based 3D geometry of lymph node | 3 states for: DCs, CD4 ^{+} T cells, CD8 ^{+} T cells; Locations for HEVs, FRCs |
Prokopiou et al. [10] | Early CD8 ^{+} T-cell response in lymph node (136 hrs) | Clonal expansion, intracellular regulation, 3D: migration, reaction-diffusion | CPM for cells, PDEs for extracellular cytokines, ODEs for intracellular factors | APCs, T-cells, IL-2, IL-2R, Tbet, Caspase, Fas (activated, non-activated) |
The model developed in this work takes into account: 1) spatial aspects of the immune response in the lymph node by means of cell and concentration distributions, 2) regulation of T lymphocytes in the lymph node including their asymmetric division and their interaction with extracellular cytokine concentrations, 3) the intracellular regulation of T cells depending on IL-2 and type I IFN, 4) the interaction of the tissue level and the whole organism through the infection level.
Conventional models of the immune response are based on ordinary differential equations, and they do not take into account spatial distributions of cells and concentrations in the lymph node. The multi-scale models previously developed (see Table 1) use a similar agent-based cell description. However they do not take into account processes such as asymmetric cell division, their interaction with IFN or the interaction of the tissue level with the organism level which is one of the key features of our multi-scale model.
One of the essential features of our modelling approach is the description of cells as individual objects which can move, divide, differentiate or die by apoptosis. Dead cells are removed from the computation domain. Cells are considered as soft spheres with their motion described by Newton’s second law for their centers. Cell fate is determined by the intracellular regulation. The spatial distribution of cytokines is described by reaction–diffusion equations, while the intracellular regulation and infection dynamics are described by ordinary differential equations conventionally accepted in mathematical immunology. A similar approach was developed and justified in modelling of other physiological systems such as blood diseases and hematopoiesis [11–17]. Here we develop it for an integrative modelling of immunological processes in their spatial context.
Methods
Biophysics of the immune response
To formulate the mathematical model, we consider a part of the lymph node, i.e., the T cell zone, which contains various cell types, mainly the antigen presenting cells (APCs) and subsets of T lymphocytes. Naive T cells and some APCs (such as plasmocytoid Dendritic Cells, pDCs) enter the node with blood flow via the High Endothelial Venules (HEVs) whereas effector and/or memory T cells, and mainly DCs and macrophages home to lymph nodes via afferent lymphatic vessels [18, 19]. Following activation with pathogens, APCs acquire a motile state that allows their translocation to the T cell zone of draining lymph node with the afferent lymph flow [20, 21]. Therefore, we assume that the influx of APCs is proportional to the level of infection in the organism. Differentiation of naive T cells into CD4 ^{+} and CD8 ^{+} T cells occurs in the thymus from progenitor T cells [22]. We suppose that they enter lymph nodes already differentiated and that there is a given influx of each cell type.
Mature CD8 ^{+} T cells (effector cells) leave the lymph node and kill infected cells. Therefore, there is a negative feedback between production of mature CD8 ^{+} T cells and the influx of APCs.
Hybrid model of cell dynamics
In our model of cell dynamics, cells are considered as individual objects that can move, divide, differentiate and die. Their behavior is determined by the surrounding cells, by intracellular regulatory networks described by ordinary differential equations and by various substances in the extracellular matrix whose concentrations are described by partial differential equations. This approach was used to model hematopoiesis and blood diseases [11–17].
- 1.
n _{ APC }(x,t) - the density of APCs in T cell zone;
- 2.
n _{ CD4}(x,t) - the density of CD4 ^{+} T cells in T cell zone (with different levels of maturity);
- 3.
n _{ CD8}(x,t) - the density of CD8 ^{+} T cells in T cell zone (with different levels of maturity);
- 4.
I _{ e }(x,t) - the concentration of IL-2 in T cell zone;
- 5.
C _{ e }(x,t) - the concentration of type I IFN in T cell zone;
- 6.
I _{ i }(t) - the intracellular concentration of IL-2-induced signalling molecules in the ith cell;
- 7.
C _{ i }(t) - the intracellular concentration of type I IFN-induced signalling molecules in the ith cell;
- 8.
N _{ ef }(t) - the total number of effector CD8 ^{+}T cells in the body;
- 9.
N _{ inf }(t) - the total number of infected cells in the body;
Cell displacement. In the model, cells are represented by individual elastic spheres. There are two mechanisms of motion of cells in the lymph node. First of all they move in a random way. This motion allows naive T cells to meet APCs which is necessary for their activation, division and differentiation. Second, each two cells, when they meet, they push each other due to a direct mechanical interaction. We consider this interaction as an elastic force acting on cells and influencing their motion. Let us describe it in more detail.
Cell division and differentiation. APC and naive T cells enter the computational domain with a given frequency if there is available space. Naive T cells move in the computational domain randomly. If they contact APC, they divide asymmetrically (Fig. 3). The distant daughter cell is similar to the mother cell, the proximal daughter cell becomes differentiated.
When the cell reaches the half of its life cycle, it will increase its size. When it divides, two daughter cells appear, the direction of the axis connecting their centers is chosen randomly from 0 to 2π. The duration of the cell cycle is 18 hours with a random perturbation of −3 to 3 hours.
We consider two levels of maturity of CD4 ^{+} T cells and three levels of CD8 ^{+} T cells. If a differentiated cell has enough IL-2 (see the next paragraph), then it divides and gives two more mature cells. Finally differentiated cells leave the lymph node. In the simulations, this means that they are removed from the computational domain.
Here I _{ i } is the intracellular concentration of signalling molecules accumulated as a consequence of IL-2 signals transmitted through transmembrane receptor IL2R downstream the signaling pathway to control the gene expression in the ith cell. The concentrations inside two different cells are in general different from each other. The first term in the right-hand side of this equation shows the cumulative effect of IL-2 signalling. The extracellular concentration I _{ e } is taken at the coordinate x _{ i } of the center of the cell. The second term describes the degradation of IL-2-induced signalling molecules inside the cell. Furthermore, n _{ T } is the number of molecules internalized by T cell receptors.
Here C _{ i } is the intracellular concentration of signalling molecules accumulated as a consequence of IFN signals transmitted through transmembrane receptor IFNR downstream the signaling pathway to control the gene expression in the i-th cell. The concentrations inside two different cells are in general different from each other. The first term in the right-hand side of this equation shows the cumulative effect of IFN signalling. The extracellular concentration C _{ e } is taken at the coordinate x _{ i } of the center of the cell. The second term describes the degradation of IFN-induced signalling molecules inside the cell.
- C1
If the concentration of activation signals induced by type I IFN, C _{ i }, is greater than some critical level \(C_{i}^{*}\) at the beginning of the cell cycle and that of I _{ i }, is smaller than the critical level \(I_{i}^{*}\), then the cell will differentiate resulting in a mature cell.
- C2
If the concentration of activation signals induced by IL-2, I _{ i }, is greater than some critical level \(I_{i}^{*}\) at the end of the cell cycle, then the cell will divide producing two more mature cells.
- C3
If \(C_{i}<C_{i}^{*}\) at the beginning of cell cycle and \(I_{i}<I_{i}^{*}\) at the end of cell cycle, then the cell will die by apoptosis and will be removed from the computational domain.
Stochastic aspects of the model.
As it is discussed above, mechanical interaction of cells results in their displacement described by equation (1) for their centers. In order to describe random motion of cells we add random variables to the cell velocity in the horizontal and vertical directions.
Duration of cell cycle is given as a random variable in the interval [ T−τ,T+τ].
Extracellular dynamics of cytokines.
Here I _{ ex } is the extracellular concentration of IL-2, D is the diffusion coefficient, W _{ IL } is the rate of its production by CD4 ^{+} T cells, and the last term in the right-hand side of this equation describes its consumption and degradation. The production rate W _{ IL } is determined by mature CD4 ^{+} T cells. We consider each such cell as a source term with a constant production rate ρ _{ IL } at the area of the cell. Let us note that we do not take into account explicitly consumption of IL-2 by immature cells in order not to introduce an additional parameter. Implicitly this consumption is taken into account in the degradation term.
Initial and boundary conditions for both concentrations IL-2 and IFN are taken zero. As before, the production rate W _{ IFN } equals ρ _{ IFN } at the area filled by APC cells and zero otherwise.
where T is their number in the lymph nodes. So the first term in the right-hand side of this equation describes production of effector cells in the lymph nodes and the second term their death in the body.
Finally, the influx of APC cells into the lymph nodes is proportional to the number of infected cells N _{ inf }.
This influx is limited by the place available in the lymph node. If there is a free place sufficient to put a cell, the new cells are added. Let us also note that the lymph nodes can increase due to infection in order to produce more effector cells.
Results
We illustrate the model performance by considering two scenarios, reflecting different spatial patterns of IL-2 and type IFN concentration fields. In the first one, both cytokines have the same diffusion coefficient D _{ IL2}=D _{ IFN }, whereas in the second case the diffusion rate of IFN is 10-fold faster. The details of the numerical implementation of the hybrid model and the parameter values used for the simulations are presented in Additional file 1: Appendix. Cell population densities and cytokine concentrations are scaled with respect to some reference values. These are determined by the cell density in the lymph node ∼10^{5}−10^{6} mm ^{−3}, the relative proportions of APCs, CD4 ^{+} T cells and CD8 ^{+} T cells [26–32] and the production rate of the cytokines (described in detail in the appendix, see Additional file 1). The considered cell numbers correspond to a computational domain in the T cell zone of about 100μ m×100μ m×100μ m.
The model presented above contains two compartments, lymph node where effector cells are produced and the body where infection develops. The lymph node is described with the hybrid model while infection development in the organism by ordinary differential equations for infected cells and for effector cells. These two compartments are coupled by means of flux of effector cells from the lymph node to the body and by the flux of APC cells to the lymph node.
The evolution of the total number of CD4 ^{+} and CD8 ^{+} T cells in the lymph node T cell zone is shown in Fig. 5. The dynamics of APC cells in the lymph node T cell zone and the effector cells in the body is shown in Fig. 6. The magnitude of the immune response is not sufficient to eradicate completely the infection. Indeed, the number of infected cells decreases but remains positive (Fig. 7). As virus infection is not cleared, the cell populations fluctuate around some constant values. Overall, the model reproduces the qualitative patterns of long-term persistent infection (experimental infections and in humans) dynamics (e.g., [34–38]). The primary clonal expansion takes about seven days and is followed by an enhanced long-term T cell response to the persistent infection. The increase in the spread of type I IFN changes the relative distributions pattern of IL-2 and IFN, so that the resulting alteration in the cytokine signalling reduces the clonal expansion and increases the overall level of the virus infection.
Cumulative numbers of key variables of the model over 113 days post infection
Model variable | D _{ IFN }=D _{ IL } | D _{ IFN }=10D _{ IL } |
---|---|---|
Number of CD4 ^{+} T cells | 27544 | 27040 |
Number of CD8 ^{+} T cells | 15194 | 14139 |
Number of APCs | 4749 | 5293 |
The infection load | 16.98 | 19.31 |
Number of N _{ ef } | 87849 | 80967 |
The net effect of the increase in the diffusion rate of type I IFN is a reduction in the clonal expansion of the T cells, in particular the effector T cells in the peripheral organs (by ∼10%) and a rise in the infection level (by ∼20%). The changes in the clonal T cell expansion are the consequence of the differences in the cytokine concentration fields, which in turn alter the timing and the sequence of the IL-2 and type I IFN signalling.
Discussion
The activation of T cell division and death is regulated by a set of signals coming from the interactions of TCRs with the MHC class I and II presented peptides, IFN and IL-2 receptors. However, the effect of type I interferon depends on the temporal sequence of the signals obtained by naïve T cells [2]. It can change from a normal activation of T cells followed by their proliferation and differentiation to an already differentiated state followed by apoptosis. We propose that the altered T cell differentiation and proliferation sequence can result from a spatial separation of the signaling events obtained by T cells, i.e. the TCR versus IFN receptor signaling, due to a generally different location of the APCs and the type I IFN concentration field. The hypothesis is formulated using the results of our study of the spatio-temporal dynamics of the T cell response to infection with the developed hybrid mathematical model integrating intracellular, and systemic levels of the immune response regulation.
Virus persistence in humans is often associated with an exhaustion of T lymphocytes. Many factors can contribute to the development of exhaustion. One of them is associated with a shift from a normal clonal expansion pathway to an altered one characterized by an early terminal differentiation of T cells. The proposed hybrid model allows us to investigate the integrative effects of numerous biophysical and biochemical parameters on the immune regulation that are beyond the scope of the existing experimental techniques.
The need for development of multi-scale integrative models in mathematical immunology is well realized [7, 9, 10]. However, implementation of such models embedded into spatial context of immune responses remains a challenge [4, 41]. Our hybrid model has been developed in a modular form to describe a range of specific interactions and regulations of the immune response. This implies that the structure and composition of the model can be adaptively changed to meet the needs of any specific studies. The present study is based upon a simplified description of the virus-target cell interaction. This block of the hybrid model can be easily refined using the existing set of lumped mathematical models. We presented the results of single runs for two different scenarios. Obviously, multiple runs will be used to address the effect of local and global variations in the parameters on the immune response dynamics in future studies.
The presented methodology for developing a hybrid model of immune processes enables integration of data and knowledge across multiple scales. This is essential for understanding the control of the lower level processes by the properties of the higher level processes. We considered here an important example of the regulation of the T cell fate [2] by the spatial structure and material properties of lymphoid organs in which the transport and cell interaction take place [4].
Mathematical modelling of human immune system represents an important challenge. Whereas a macroscopic level population dynamics of the immune response can be followed experimentally and captured with mathematical models under the framework of mono- or two/three compartmental modelling based on ODEs, the embedding of dynamics and interactions of various cells into the spatial organization of real lymphoid organs represents a fundamental challenge. The effects of spatial heterogeneity of the cytokine fields and cells localization in the lymphoid organs on quantitative and phenotypic features of the T cell responses are poorly understood as these types of parameters are difficult to examine experimentally. Indeed, to describe the cytokine and cellular interactions one need to integrate in a unified manner a range of processes including the lymph flow, constrained diffusion, chemotaxis and haptotaxis in soft tissues. A complementary problem is the assimilation of the data on intracellular regulation of immune cells’ states resulting from a parallel signaling via cell surface receptors to antigens, cytokines, chemokines, hormones, etc. The complexity of the issue is well exemplified in a recent review [42]. Following the presented modelling approach, we hope to be able to examine the sensitivity of the immune responses at the macroscopic level to parameters of the cell interaction at the microscopic level. This should assist in fine tuning of the offset dynamics of immune responses using a broad spectrum of modern immunomodulatory drugs.
Conclusions
The proposed hybrid mathematical model of the immune response represents a novel analytical tool to examine challenging issues in the spatio-temporal regulation of cell growth and differentiation, in particular the effect of timing and location of activation signals. It allows us to overcome the limitations of reductionist approach to a single factor analysis of the immune regulation and to proceed to the studies of the structure-function relationships in a genuinely systemic manner. As W.E. Paul stated in [33] “…the behavior of immune cells is highly colored by the cellular/molecular environment in which they exist…It is to the quantitative prediction of the outcome of given perturbations in the immune system that we envisage our mathematical/ modeling colleagues will apply themselves.” In our view, hybrid modelling approach provides the means for a comprehensive analysis and interpretation of content-rich meta-data obtained by a broad range of scale-specific data acquisition techniques, including imaging, flow cytometry, transcriptome sequencing in an anatomically correct and immunologically meaningful way.
The practical implementation of a hybrid approach to multiscale modelling presents a number of challenges ranging from the numerical accuracy and consistency of the different methods being used to compute the system components dynamics on one side to the risk of producing modelling artefacts because of the system complexity and parameter uncertainty on the other side. In fact, one needs to have clear computational methodologies for the development of various mathematical modelling tools including simple single-level resolution phenomenological models, large-scale multi-compartmental models and high-resolution multiscale models. In the end, their appropriateness depends on the questions to study and available data.
Declarations
Acknowledgments
This study was supported by the grant of the Russian Science Foundation, project no. 15-11-00029 (to GB, AM and VV).
Funding
The publication charge for this article was funded by the Russian Science Foundation, project no. 15-11-00029. The funder had no role in study design, analysis, decision to publish, or preparation of the manuscript.
Availability of data and material
Data sharing not applicable to this article as no datasets were generated or analysed during the current study.
Authors’ contributions
All authors of this paper have contributed to the study design. AB and VV developed the hybrid model and conducted the numerical simulations. AB, GB, AM, VV conceived the study design and manuscript preparation. All authors read and approved the final manuscript.
Competing interests
The authors declare that they have no competing interests.
Consent for publication
Not applicable.
Ethics approval and consent to participate
Not applicable.
About this Supplement
This article has been published as part of BMC Immunology Volume 18 Supplement 1, 2017. Systems Immunology & ImmunoInformatics. The full contents of the supplement are available online https://bmcimmunol.biomedcentral.com/articles/supplements/volume-18-supplement-1.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver(http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
Authors’ Affiliations
References
- Khan SH, Martin MD, Starbeck-Miller GR, Xue H-H, Harty JT, Badovinac VP. The Timing of Stimulation and IL-2 Signaling Regulate Secondary CD8 ^{+} T Cell Responses. PLoS Pathog. 2015; 11(10):e1005199.View ArticlePubMedPubMed CentralGoogle Scholar
- Welsh RM, Bahl K, Marshall HD, Urban SL. Type 1 Interferons and Antiviral CD8 ^{+} T-Cell Responses. PLoS Pathog. 2012; 8(1):e1002352.View ArticlePubMedPubMed CentralGoogle Scholar
- Hervas-Stubbs S, Perez-Gracia JL, Rouzaut A, Sanmamed MF, Le Bon A, Melero I. Direct Effects of Type I Interferons on Cells of the Immune System. Clin Cancer Res. 2011; 17:2619–27.View ArticlePubMedGoogle Scholar
- Ludewig B, Stein JV, Sharpe J, Cervantes-Barragan L, Thiel V, Bocharov G. A global “imaging” view on systems approaches in immunology. Eur J Immunol. 2012; 42:3116–25.View ArticlePubMedGoogle Scholar
- Palsson S, Hickling TP, Bradshaw-Pierce EL, Zager M, Jooss K, O’Brien PJ, Spilker ME, Palsson BO, Vicini P. The development of a fully-integrated immune response model (FIRM) simulator of the immune response through integration of multiple subset models. BMC Syst Biol. 2013; 7:95. doi:10.1186/1752-0509-7-95.View ArticlePubMedPubMed CentralGoogle Scholar
- Cilfone NA, Kirschner DE, Linderman JJ. Strategies for efficient numerical implementation of hybrid multi-scale agent-based models to describe biological systems. Cel Mol Bioeng. 2015; 8(1):119–36.View ArticleGoogle Scholar
- Baldazzi V, Paci P, Bernaschi M, Castiglione F. Modeling lymphocyte homing and encounters in lymph nodes. BMC Bioinforma. 2009; 10:387.View ArticleGoogle Scholar
- Fallahi-Sichani M, El-Kebir M, Marino S, Kirschner DE, Linderman JJ. Multi-scale computational modeling reveals a critical role for TNF receptor 1 dynamics in tuberculosis granuloma formation. J Immunol. 2011; 186(6):3472–83.View ArticlePubMedPubMed CentralGoogle Scholar
- Gong C, Mattila JT, Miller M, Flynn JL, Linderman JJ, Kirschner D. Predicting lymph node output efficiency using systems biology. J Theor Biol. 2013; 335:169–84.View ArticlePubMedPubMed CentralGoogle Scholar
- Prokopiou SA, Barbarroux L, Bernard S, Mafille J, Leverrier Y, Arpin C, Marvel J, Gandrillon O, Crauste F. Multiscale Modeling of the Early CD8 ^{+} T T-Cell Immune Response in Lymph Nodes:An Integrative Study. Computation. 2014; 2:159–81.Google Scholar
- Bessonov N, Eymard N, Kurbatova P, Volpert V. Mathematical modelling of erythropoiesis in vivo with multiple erythroblastic islands. Appl Math Lett. 2012; 25:1217–21.View ArticleGoogle Scholar
- Fischer S, Kurbatova P, Bessonov N, Gandrillon O, Volpert V, Crauste F. Modelling erythroblastic islands: using a hybrid model to assess the function of central macrophage. J Theoretical Biol. 2012; 298:92–106.View ArticleGoogle Scholar
- Kurbatova P, Bernard S, Bessonov N, Crauste F, Demin I, Dumontet C, Fischer S, Volpert V. Hybrid model of erythropoiesis and leukemia treatment with cytosine arabinoside. SIAM J Appl Math. 2011; 71(6):2246–68.View ArticleGoogle Scholar
- Volpert V, Bessonov N, Eymard N, Tosenberger A. Modèle multi-échelle de la dynamique cellulaire. In: Glade N, et d’Angelique Stephanou, editeurs. Le vivant discret et continu. Paris: Éditions Matériologiques: 2013. p. 91–111.Google Scholar
- Kurbatova P, Eymard NV, V. Volpert V. Hybrid Model of Erythropoiesis. Acta Biotheoretica. 2013; 61:305–15.View ArticlePubMedGoogle Scholar
- Eymard N, Bessonov N, Gandrillon O, Koury MJ, Volpert V. The role of spatial organization of cells in erythropoiesis. J Math Biol. 2015; 70:71–97.View ArticlePubMedGoogle Scholar
- Stéphanou A, Volpert V. Hybrid modelling in biology: a classification review. Math Model Nat Phenom. 2016; 11(1):37–48.View ArticleGoogle Scholar
- Mueller SN, Germain RN. Stromal cell contributions to the homeostasis and functionality of the immune system. Nat Rev Immunol. 2009; 9(9):618–29.PubMedPubMed CentralGoogle Scholar
- Girard JP, Moussion C, Förster R. HEVs, lymphatics and homeostatic immune cell trafficking in lymph nodes. Nat Rev Immunol. 2012; 12(11):762–73.View ArticlePubMedGoogle Scholar
- Junt T, Scandella E, Ludewig B. Form follows function: lymphoid tissue microarchitecture in antimicrobial immune defence. Nat Rev Immunol. 2008; 8(10):764–75.View ArticlePubMedGoogle Scholar
- Förster R, Braun A, Worbs T. Lymph node homing of T cells and dendritic cells via afferent lymphatics. Trends Immunol. 2012; 33(6):271–80.View ArticlePubMedGoogle Scholar
- Goldsby RA, Kindt TJ, Osborne BA. Kuby Immunology. New York: W H Freeman & Co (Sd); 2000. Ch.10:222.Google Scholar
- Chang JT, Reiner SL. Asymmetric Division and Stem Cell Renewal without a Permanent Niche: Lessons from Lymphocytes. Cold Spring Harb Symp Quant Biol. 2008; 73:73–9.View ArticlePubMedGoogle Scholar
- Broere F, Apasov SG, Sitkovsky MV, van Eden W. T cell subsets and T cell-mediated immunity In: Nijkamp FP, Parnham MJ, editors. Principles of Immunopharmacology: 3rd revised and extended edition. Basel: Springer AG: 2011.Google Scholar
- Nelson BH. IL-2, Regulatory T Cells, and Tolerance. J Immunol. 2004; 172:3983–8.View ArticlePubMedGoogle Scholar
- Ganusov VV, De Boer RJ. Do most lymphocytes in humans really reside in the gut?Trends Immunol. 2007; 28(12):514–8.View ArticlePubMedGoogle Scholar
- Scandella E, Bolinger B, Lattmann E, Miller S, Favre S, Littman DR, Finke D, Luther SA, Junt T, Ludewig B. Restoration of lymphoid organ integrity through the interaction of lymphoid tissue-inducer cells with stroma of the T cell zone. Nat Immunol. 2008; 9(6):667–75.View ArticlePubMedGoogle Scholar
- Kumar V, Scandella E, Danuser R, Onder L, Nitschké M, Fukui Y, Halin C, Ludewig B, Stein JV. Global lymphoid tissue remodeling during a viral infection is orchestrated by a B cell-lymphotoxin-dependent pathway. Blood. 2010; 115(23):4725–33.View ArticlePubMedGoogle Scholar
- Bocharov G, Züst R, Cervantes-Barragan L, Luzyanina T, Chiglintsev E, Chereshnev VA, et al. A Systems Immunology Approach to Plasmacytoid Dendritic Cell Function in Cytopathic Virus Infections. PLoS Pathog. 2010; 6(7):e1001017.View ArticlePubMedPubMed CentralGoogle Scholar
- Bocharov G, Chereshnev V, Gainova I, Bazhan S, Bachmetyev B, Argilaguet J, Martinez J, Meyerhans A. Human Immunodeficiency Virus Infection: from Biological Observations to Mechanistic Mathematical Modelling. Math Model Nat Phenom. 2012; 7(5):78–104.View ArticleGoogle Scholar
- Cremasco V, Woodruff MC, Onder L, Cupovic J, Nieves-Bonilla JM, Schildberg FA, Chang J, Cremasco F, Harvey CJ, Wucherpfennig K, Ludewig B, Carroll MC, Turley SJ. B cell homeostasis and follicle confines are governed by fibroblastic reticular cells. Nat Immunol. 2014; 15(10):973–81.View ArticlePubMedPubMed CentralGoogle Scholar
- Giese C, Marx U. Human immunity in vitro - solving immunogenicity and more. Adv Drug Deliv Rev. 2014; 69–70:103–22.View ArticlePubMedGoogle Scholar
- Paul WE. The Immune System - Complexity Exemplified. Math Model Nat Phenom. 2012; 7(5):4–6.View ArticleGoogle Scholar
- Moskophidis D, Lechner F, Pircher H, Zinkernagel RM. Virus persistence in acutely infected immunocompetent mice by exhaustion of antiviral cytotoxic effector T cells. Nature. 1993; 362(6422):758–61.View ArticlePubMedGoogle Scholar
- Moskophidis D, Battegay M, van den Broek M, Laine E, Hoffmann-Rohrer U. Zinkernagel R. Role of virus and host variables in virus persistence or immunopathological disease caused by a non-cytolytic virus. J. Gen. Virol. 1995; 76(2):381–391.PubMedGoogle Scholar
- McMichael AJ, Borrow P, Tomaras GD, Goonetilleke N, Haynes BF. The immune response during acute HIV-1 infection: clues for vaccine development. Nat Rev Immunol. 2010; 10(1):11–23.View ArticlePubMedGoogle Scholar
- Cervantes-Barragan L, Lewis KL, Firner S, Thiel V, Hugues S, Reith W, Ludewig B, Reizis B. Plasmacytoid dendritic cells control T-cell response to chronic viral infection. Proc Natl Acad Sci USA. 2012; 109(8):3012–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Hansen SG, Piatak M, Ventura AB, Colette M, Hughes CM, Gilbride RM, Ford JC, Oswald K, Shoemaker R, Li Y, Lewis MS, Gilliam AN, Xu G, Whizin N, Burwitz BJ, Planer SL, Turner JM, Legasse AW, Axthelm MK, Nelson JA, Früh K, Sacha JB, Estes JD, Keele BF, Edlefsen PT, Lifson JD, Picker LJ. Immune clearance of highly pathogenic SIV infection. Nature. 2013; 502(7469):100–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Kaech SM, Wherry EJ. Heterogeneity and cell-fate decisions in effector and memory CD8+ T cell differentiation during viral infection. Immunity. 2007; 27(3):393–405.View ArticlePubMedPubMed CentralGoogle Scholar
- Fricke GM, Letendre KA, Moses ME, Cannon JL. Persistence and Adaptation in Immunity: T Cells Balance the Extent and Thoroughness of Search. PLoS Comput Biol. 2016; 12(3):e1004818.View ArticlePubMedPubMed CentralGoogle Scholar
- Cappuccio A, Tieri P, Castiglione F. Multiscale modelling in immunology: a review. Brief Bioinform. 2015; pii:bbv012:1–11.Google Scholar
- Grossman Z, Paul WE. Dynamic tuning of lymphocytes: physiological basis, mechanisms, and function. Annu Rev Immunol. 2015; 33:677–713.View ArticlePubMedGoogle Scholar
- Baker CTH, Bocharov GA, Paul CAH. Mathematical Modelling of the Interleukin-2 T-Cell System: A Comparative Study of Approaches Based on Ordinary and Delay Differential Equation. J Theor Medicine. 1997; 1(2):117–28.View ArticleGoogle Scholar