Abstract
Human immunodeficiency virus (HIV) is a fast replicating ribonucleic acid virus, which can easily mutate in order to escape the effects of drug administration. Hence, understanding the basic mechanisms underlying HIV persistence in the body is essential in the development of new therapies that could eradicate HIV infection. Lymphoid tissues are the primary sites of HIV infection. Despite the recent progress in real-time monitoring technology, HIV infection dynamics in a whole body is unknown. Mathematical modeling and simulations provide speculations on global behavior of HIV infection in the lymphatic system. We propose a new mathematical model that describes the spread of HIV infection throughout the lymphoid tissue network. In order to represent the volume difference between lymphoid tissues, we propose the proportionality of several kinetic parameters to the lymphoid tissues’ volume distribution. Under this assumption, we perform extensive numerical computations in order to simulate the spread of HIV infection in the lymphoid tissue network. Numerical computations simulate single drug treatments of an HIV infection. One of the important biological speculations derived from this study is a drug saturation effect generated by lymphoid network connection. This implies that a portion of reservoir lymphoid tissues to which drug is not sufficiently delivered would inhibit HIV eradication despite of extensive drug injection.
Similar content being viewed by others
References
Agosto LM, Zhong P, Munro J, Mothes W (2014) Highly active antiretroviral therapies are effective against HIV-1 cell-to-cell transmission. PLoS Pathog 10:e1003982
Alberts B, Johnson A, Lewis J, Raff M, Roberts K, Walter P (2008) Molecular biology of the cell 5E. Garland Science, New York
Aloisi F, Pujol-Borrell R (2006) Lymphoid neogenesis in chronic inflammatory diseases. Nat Rev Immunol 6:205-217
Arron ST, Ribeiro RM, Gettie A, Bohm R, Blanchard J, Yu J, Perelson AS, Ho DD, Zhang L (2005) Impact of thymectomy on the peripheral T cell pool in rhesus macaques before and after infection with simian immunodeficiency virus. Eur J Immunol 35:46-55
Bajaria SH, Webb G, Cloyd M, Kirschner D (2002) Dynamics of naive and memory CD4+ T lymphocytes in HIV-1 disease progression. J Acquir Immune Defic Syndr 30:41-58
Bajaria SH, Webb G, Kirschner DE (2004) Predicting differential responses to structured treatment interruptions during HAART. Bull Math Biol 66:1093-1118
Bisset LR, Lung TL, Kaelin M, Ludwig E, Dubs RW (2004) Reference values for peripheral blood lymphocyte phenotypes applicable to the healthy adult population in Switzerland. Eur J Haematol 72:203-212
Boyman O, Létourneau S, Krieg C, Sprent J (2009) Homeostatic proliferation and survival of naïve and memory T cells. Eur J Immunol 39:2088-2094
Bronte V, Pittet MJ (2013) The spleen in local and systemic regulation of immunity. Immunity 39:806-818
Casteleyn C, Cornillie P, Van Ginneken C, Simoens P, Van Cruchten S, Vandevelde K, Van den Broeck W (2014) Lymph drainage from the ovine tonsils: an anatomical study of the tonsillar lymph vessels. Anat Histol Embryol 43:482-489
Catalfamo M, Wilhelm C, Tcheung L, Proschan M, Friesen T, Park J-H, Adelsberger J, Baseler M, Maldarelli F, Davey R, Roby G, Rehm C, Lane C (2011) CD4 and CD8 T cell immune activation during chronic HIV infection: roles of homeostasis, HIV, type I IFN, and IL-7. J Immunol 186:2106-2116
Chen HY, Di Mascio M, Perelson AS, Ho DD, Zhang L (2007) Determination of virus burst size in vivo using a single-cycle SIV in rhesus macaques. Proc Natl Acad Sci USA 104:19079-19084
Core Team R (2014) R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna
Costiniuk CT, Jenabian M-A (2014) Cell-to-cell transfer of HIV infection: implications for HIV viral persistence. J Gen Virol 95:2346-2355
Davenport MP, Ribeiro RM, Zhang L, Wilson DP, Perelson AS (2007) Understanding the mechanisms and limitations of immune control of HIV. Immunol Rev 216:164-175
De Boer RJ, Perelson AS (2013) Quantifying T lymphocyte turnover. J Theor Biol 327:45-87
De Boer RJ, Ribeiro RM, Perelson AS (2010) Current estimates for HIV-1 production imply rapid viral clearance in lymphoid tissues. PLoS Comput Biol 6:e1000906
den Braber I, Mugwagwa T, Vrisekoop N, Westera L, Mögling R, de Boer AB, Willems N, Schrijver EHR, Spierenburg G, Gaiser K, Mul E, Otto SA, Ruiter AFC, Ackermans MT, Miedema F, Borghans JAM, de Boer RJ, Tesselaar K (2012) Maintenance of peripheral naive T cells is sustained by thymus output in mice but not humans. Immunity 36:288-297
Diekmann O, Heesterbeek JAP, Metz JAJ (1990) On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases. J Math Biol 35:503-522
Eberl G (2005) Inducible lymphoid tissues in the adult gut: recapitulation of a fetal developmental pathway? Nat Rev Immunol 5:413-420
Fletcher CV, Staskus K, Wietgrefe SW, Rothenberger M, Reilly C, Chipman JG, Beilman GJ, Khoruts A, Thorkelson A, Schmidt TE, Anderson J, Perkey K, Stevenson M, Perelson AS, Douek DC, Haase AT, Schacker TW (2014) Persistent HIV-1 replication is associated with lower antiretroviral drug concentrations in lymphatic tissues. Proc Natl Acad Sci USA 111:2307-2312
Frost SD, Dumaurier MJ, Wain-Hobson S, Brown AJ (2001) Genetic drift and within-host metapopulation dynamics of HIV-1 infection. Proc Natl Acad Sci USA 98:6975-6980
Graw F, Regoes RR (2012) Influence of the fibroblastic reticular network on cell-cell interactions in lymphoid organs. PLoS Comput Biol 8:e1002436
Gyllenberg M, Hanski I (1992) Single-species metapopulation dynamics: a structured model. Theor Popul Biol 42:35-61
Gyllenberg M, Metz JA (2001) On fitness in structured metapopulations. J Math Biol 43:545-560
Hale JS, Fink PJ (2009) Back to the thymus: peripheral T cells come home. Immunol Cell Biol 87:58-64
Hanski I, Gilpin M (1997) Metapopulation biology: ecology, genetics and evolution. Academic Press, London
Harris M, Patenaude P, Cooperberg P, Filipenko D, Thorne A, Raboud J, Rae S, Dailey P, Chernoff D, Todd J, Conway B, Montaner JS (1997) Correlation of virus load in plasma and lymph node tissue in human immunodeficiency virus infection. INCAS Study Group. Italy, Netherlands, Canada, Australia, and (United) States. J Infect Dis 176:1388-1392
Hogue IB, Bajaria SH, Fallert BA, Qin S, Reinhart TA, Kirschner DE (2008) The dual role of dendritic cells in the immune response to human immunodeficiency virus type 1 infection. J Gen Virol 89:2228-2239
Jilek BL, Zarr M, Sampah ME, Rabi SA, Bullen CK, Lai J, Shen L, Siliciano RF (2012) A quantitative basis for antiretroviral therapy for HIV-1 infection. Nat Med 18:446-451
Kirschner D, Webb GF, Cloyd M (2000) Model of HIV-1 disease progression based on virus-induced lymph node homing and homing-induced apoptosis of CD4+ lymphocytes. J Acquir Immune Defic Syndr 24:352-362
Kodera M, Grailer JJ, Karalewitz AP-A, Subramanian H, Steeber DA (2008) T lymphocyte migration to lymph nodes is maintained during homeostatic proliferation. Microsc Microanal 14:211-224
Komarova NL, Anghelina D, Voznesensky I, Trinité B, Levy DN, Wodarz D (2013) Relative contribution of free-virus and synaptic transmission to the spread of HIV-1 through target cell populations. Biol Lett 9:20121049
Markowitz M, Louie M, Hurley A, Sun E, Di Mascio M, Perelson AS, Ho DD (2003) A novel antiviral intervention results in more accurate assessment of human immunodeficiency virus type 1 replication dynamics and T-cell decay in vivo. J Virol 77:5037-5038
Metz JA, Gyllenberg M (2001) How should we define fitness in structured metapopulation models? Including an application to the calculation of evolutionarily stable dispersal strategies. Proc Biol Sci 268:499-508
Mueller SN, Germain RN (2009) Stromal cell contributions to the homeostasis and functionality of the immune system. Nat Rev Immunol 9:618-629
Murooka TT, Deruaz M, Marangoni F, Vrbanac VD, Seung E, von Andrian UH, Tager AM, Luster AD, Mempel TR (2012) HIV-infected T cells are migratory vehicles for viral dissemination. Nature 490:283-287
Murray JM, Kaufmann GR, Hodgkin PD, Lewin SR, Kelleher AD, Davenport MP, Zaunders JJ (2003) Naive T cells are maintained by thymic output in early ages but by proliferation without phenotypic change after age twenty. Immunol Cell Biol 81:487-495
Pabst O (2012) New concepts in the generation and functions of IgA. Nat Rev Immunol 12:821-832
Pellas TC, Weiss L (1990) Deep splenic lymphatic vessels in the mouse: a route of splenic exit for recirculating lymphocytes. Am J Anat 187:347-354
Perelson AS, Nelson PW (1999) Mathematical analysis of HIV-1 dynamics in vivo. SIAM Rev 41:3-44
Perelson AS, Neumann AU, Markowitz M, Leonard JM, Ho DD (1996) HIV-1 dynamics in vivo: virion clearance rate, infected cell life-span, and viral generation time. Science 271:1582-1586
Perelson AS, Essunger P, Cao Y, Vesanen M, Hurley A, Saksela K, Markowitz M, Ho DD (1997) Decay characteristics of HIV-1-infected compartments during combination therapy. Nature 387:188-191
Qatarneh SM, Kiricuta I-C, Brahme A, Tiede U, Lind BK (2006) Three-dimensional atlas of lymph node topography based on the visible human data set. Anat Rec B New Anat 289:98-111
Rabi SA, Laird GM, Durand CM, Laskey S, Shan L, Bailey JR, Chioma S, Moore RD, Siliciano RF (2013) Multi-step inhibition explains HIV-1 protease inhibitor pharmacodynamics and resistance. J Clin Invest 123:3848-3860
Ramratnam B, Bonhoeffer S, Binley J, Hurley A, Zhang L, Mittler JE, Markowitz M, Moore JP, Perelson AS, Ho DD (1999) Rapid production and clearance of HIV-1 and hepatitis C virus assessed by large volume plasma apheresis. Lancet 354:1782-1785
Ribeiro RM, Perelson AS (2007) Determining thymic output quantitatively: using models to interpret experimental T-cell receptor excision circle (TREC) data. Immunol Rev 216:21-34
Ruddle NH (2014) Lymphatic vessels and tertiary lymphoid organs. J Clin Invest 124:953-959
Sathaliyawala T, Kubota M, Yudanin N, Turner D, Camp P, Thome JJC, Bickham KL, Lerner H, Goldstein M, Sykes M, Kato T, Farber DL (2013) Distribution and compartmentalization of human circulating and tissue-resident memory T cell subsets. Immunity 38:187-197
Shen L, Peterson S, Sedaghat AR, McMahon MA, Callender M, Zhang H, Zhou Y, Pitt E, Anderson KS, Acosta EP, Siliciano RF (2008) Dose-response curve slope sets class-specific limits on inhibitory potential of anti-HIV drugs. Nat Med 14:762-766
Sprent J (1973) Circulating T and B lymphocytes of the mouse. I. Migratory properties. Cell Immunol 7:40-59
Sprent J, Basten A (1973) Circulating T and B lymphocytes of the mouse. II. Lifespan. Cell Immunol 7:10-39
Sprent J, Surh CD (2009) Re-entry of mature T cells to the thymus: an epiphenomenon? Immunol Cell Biol 87:46-49
Stacker SA, Williams SP, Karnezis T, Shayan R, Fox SB, Achen MG (2014) Lymphangiogenesis and lymphatic vessel remodelling in cancer. Nat Rev Cancer 14:159-172
Tan JT, Dudl E, LeRoy E, Murray R, Sprent J, Weinberg KI, Surh CD (2001) IL-7 is critical for homeostatic proliferation and survival of naive T cells. Proc Natl Acad Sci USA 98:8732-8737
Textor J, Henrickson SE, Mandl JN, von Andrian UH, Westermann J, de Boer RJ, Beltman JB (2014) Random migration and signal integration promote rapid and robust T cell recruitment. PLoS Comput Biol 10:e1003752
Turley SJ, Fletcher AL, Elpek KG (2010) The stromal and haematopoietic antigen-presenting cells that reside in secondary lymphoid organs. Nat Rev Immunol 10:813-825
van den Driessche P, Watmough J (2002) Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Math Biosci 180:29-48
Vrisekoop N, den Braber I, de Boer AB, Ruiter AFC, Ackermans MT, van der Crabben SN, Schrijver EHR, Spierenburg G, Sauerwein HP, Hazenberg MD, de Boer RJ, Miedema F, Borghans JAM, Tesselaar K (2008) Sparse production but preferential incorporation of recently produced naive T cells in the human peripheral pool. Proc Natl Acad Sci USA 105:6115-6120
Weinreich MA, Hogquist KA (2008) Thymic emigration: when and how T cells leave home. J Immunol 181:2265-2270
Westermann J, Pabst R (1990) Lymphocyte subsets in the blood: a diagnostic window on the lymphoid system? Immunol Today 11:406-410
Acknowledgments
The authors are grateful to the referees for constructive and helpful comments and suggestions, which led to significant improvement of our original manuscript.
Author information
Authors and Affiliations
Corresponding author
Additional information
Dedicated to Mats Gyllenberg on the occasion of his 60th birthday.
This research was partly supported by (i) the Japan Society for the Promotion of Science (JSPS) through the “Grant-in-Aid for Young Scientists B25871132 (to S.N.) and B25800092 (to S.I.)”, and received funding support from the Commissioned Research program of the Ministry of Health, Labour and Welfare, Japan (to S.N., H26-ShinkoJitsuyoka-General-016).
Appendix: Supplementary information
Appendix: Supplementary information
1.1 Quasi-steady state approximation
Let us consider a situation in which a virion’s turnover rate is significantly faster than that of a CD4 T cell. While we do not apply quasi-steady state approximation to the blood component. Then, quasi-steady state approximation can be applied to the third equation in (2.9). That is,
Equation (A.1) for \(j=1,2,\ldots ,N\) are then written:
If \(\det (\mathbf {I}-\frac{\rho _V{{\mathbf {M}}}_{\mathbf {L}}}{c}) \ne 0\), then
Assume that \(\rho _V{{\mathbf {M}}}_{\mathbf {L}}/c \approx {\mathbf {O}}\), where \({\mathbf {O}}\) is the zero matrix. Then we have
Substituting (A.4) into (2.9), system (2.9) is reduced to the following \(2\times N\)-dimensional ordinary differential equations:
Let \(E_{c}=({\bar{T}}_1,\ldots ,{\bar{T}}_N,{\bar{T}}_B,0,\ldots ,0)\) denote an infectious-free equilibrium for (A.5) with (2.8), and let T(t) and I(t) denote the total number of uninfected and infected CD4 T cells in the lymphoid tissue network at time t, respectively. That is,
Note that the total flux (input and output) of migrated cells between the lymphoid tissues is conserved. More precisely, it follows from (2.6) that
Adding the right hand side of system (A.5) and the first equation of (2.8) at this equilibrium, we find
Note that any solution of (A.8) with positive initial condition converges to \(\lambda /d\) as \(t \rightarrow \infty \). The explicit values of the components of \(E_c\) correspond to the stable equilibrium of the following system of ordinary differential equations
with the constraint:
1.2 Derivation of basic reproduction number
We derive a next generation matrix for (A.5), which will determine the possibility of HIV persistence. The next generation matrix is defined for the linearized system in (A.5) around \(E_{c}\). Define matrices \({\mathcal {F}}\) and \({\mathcal {V}}\) to be
and
The basic reproduction number is given by the spectral radius of next generation matrix \({\mathcal {F}}{\mathcal {V}}^{-1}\) (Diekmann et al. 1990; Driessche and Watmough 2002):
The basic reproduction number \(R_{0,\mathrm{{LN}}}\) is numerically calculated using the predefined function \(\mathtt {eigen()}\) in the statistical computation software \(\mathtt {R}\) (R Core Team 2014).
1.3 Two lymphoid tissues model
Numerical computations in Sect. 4 implies existence of the lower limit of \(R_{0,LN}\). To analytically prove existence of the lower limit, we consider a simplified system in which only two lymphoid tissues exist and are connected each other. The corresponding system of differential equations is given by
Let \(E_0=({\bar{T}}_1,\bar{T_2},0,0)\) denote an equilibrium in which no infection occurs. In the similar way in “Derivation of basic reproduction number”, matrices \({\mathcal {F}}_2\) and \({\mathcal {V}}_2\) for system (A.14) are explicitly given by
and
The basic reproduction number defined for \(E_0\) of system (A.14) is given by the spectral radius of next generation matrix \({\mathcal {F}}_2{\mathcal {V}}_{2}^{-1}\):
Note that \(R_{0,2}\) is explicitly given as
Note that \(R_{0,2}\) is monotonically increasing with respect to \(\alpha \):
Consider a situation that drug treatment reduces \(\alpha \). By taking the limit \(\alpha \rightarrow 0\), \(R_{0,2}\) converges to
If migration rates \(m_1\) and \(m_2\) are sufficiently larger than \(\delta \), then
(A.21) indicates that the basic reproduction number is less than 1 but it is close to 1 if cellular/viral migration are sufficiently fast. For instance, we use the following numerical values to calculate the limit of \(R_{0,2}\) with respect to the limit \(\alpha \rightarrow 0\):
Then \(\bar{R}_{0,2} \simeq 0.8\). The threshold value of \(\alpha \) which determines whether \(R_{0,2}>1\) or not is \(\alpha ^* \simeq 0.12\). To simulate drug therapy, numerical values of basic reproduction number \(R_{0,2}\) are calculated with respect to (A.22) and different \(\alpha \). The left panel of Fig. 8 shows a decline curve of \(R_{0,2}\) with respect to the decrease of \(\alpha \). Saturation effect described as the plateau of the curve at low values of \(\alpha \) exists. Moreover, numerical simulation results with \(\alpha =0.5>\alpha ^*\) and \(\alpha =0.1<\alpha ^*\) shown in the right panel of Fig. 8 indicate that \(R_{0,2}\) plays a role as the threshold to determine whether infected cells persist or not. Importantly, although parameter values used in numerical computations of \(R_{0,2}\) have less biological meaning, qualitatively similar characteristics to the network model under single drug treatment therapy in Sect. 4 are reproduced with two lymphoid tissues model (A.14).
1.4 Scaling
The units of variables T, I and V are (cells/ml), (cells/ml) and (copies/ml), respectively (see Table 1 in Appendix). Because the total amount of blood accounts for about 8 % of total body mass, the average man (60 kg) contains approximately 5 l (5 kg) of blood (Alberts et al. 2008, Table 22-1 Blood Cells). We assume that about 2 % of white blood cells are found in the peripheral blood (Westermann and Pabst 1990; Vrisekoop et al. 2008). In addition, because the average number of lymphocytes in the peripheral blood is \(1.0\times 10^{6}\) cells/ml, the total number of lymphocytes in a body is estimated to be
Let x, y and z denote the number of uninfected CD4 T cells, infected CD4 T cells and virions in a body, respectively. Define the conversion ratio \(\gamma \) to be
Then the relationship between any two variables is
Parameters \(\beta \) and \(\omega \) in the (T, I, V)-system are the product of \(\gamma \) and the corresponding parameters in the un-scaled system. The other parameters remain unchanged.
1.5 Parameter values
The default values for the parameters are summarized in Table 1 in Appendix. We use the reference values \(d=0.01\) and \(\delta =0.7\) for the per capita death rates of the uninfected and infected CD4 T cells, respectively. The normal count of CD4 T cells contained in 1 ml of peripheral blood is approximately \(1.0\times 10^6\) in a healthy individual. Since about 2 % of white blood cells are found in blood, we assume that the total number of CD4 T cells is in a steady state \(\lambda /d\) under the absence of infection. Hence, we have \(\lambda =5.0\times 10^5\) and \(T(0)+T_B(0)=\lambda /d\). We need to point out the limitation of the current model formulation. [R2-1] The current model does not distinguish CD4 T cell subsets. In quantitative estimates of human immune cell turnover rates using in vivo kinetic deuterium labeling (Vrisekoop et al. 2008), an estimate for naive T cell turnover rate is in the order of 0.001 [/day]. Hence the choice of the numerical value \(d=0.01\) in Table 1 in Appendix might be inappropriate. Regardless of the choice of \(\lambda \) and d, however, our numerical simulation results under the absence of infection shows consistent results with clinical observations that a T cell level in the peripheral blood is \(1.0 \times 10^6\) per ml. Values \(\beta \) and \(\omega \) vary. Because the HIV infection’s basic reproduction number is estimated to be between 8 and 10, the ranges of \(\beta \) and \(\omega \) are determined so that they cover the estimated value for \(R_{0,\mathrm{{LN}}}\). Numerical values of free and cell-to-cell infection rate \(\beta \) and \(\omega \) are kept in the same order for the following reasons. In Komarova et al. (2013), in vitro cell culture system was used to infect target cells (Jurkat cells expressing CXCR4) by HIV. Experimental data of a single round HIV infection with/without shaking were fitted with their mathematical model which incorporates two sources of infection: free and cell-to-cell transmission. The ratio of the free to cell-to-cell infection rate (\(\omega /\beta \) in our notation) for several experiments indicated that estimated values take 1.0 on average, suggesting that the free and cell-to-cell infection rates are in the same order (Table 1 in Komarova et al. 2013). The number of HIV virions produced from one infected cell during its life expectancy is known as its burst size. Because the life expectancy of infected cells is \(1/\delta \simeq 1\), the burst size roughly corresponds to the per capita production rate of the virions p. According to the estimate in Chen et al. (2007), burst size in vivo is approximately \(1.0 \times 10^4\) per ml. In addition, given \(1.0 \times 10^4\) infected cells per ml, it follows from the burst size estimate that the estimated number of virions is \(1.0 \times 10^8\). This contradicts clinical observations, in which the estimated number of mRNA copies of HIV is \(1.0 \times 10^{5}\) on average. One possible interpretation of this is given in Boer and Perelson (2013), where the clearance of HIV is suggested to be faster than the usual estimates obtained using lymphoid tissue. In this paper, we adjust the value of p in order to produce clinically relevant observations. The decay rate of virions \(c = 30.0\) is also taken from past literature. Because the volume distribution of lymphoid tissues is currently unknown, the exponential distribution was assumed in our study. Additionally, because the total number of infected CD4 T cells and particles of HIV at infection initiation are also unknown, we use the numerical values in Bajaria et al. (2004) as our reference: \(I^0=I(0)=1.0\) (assumed) and \(V^0=V(0)=10\) (Bajaria et al. 2004). To determine the lymphoid tissue structure and migration rate of cells/virions, numerical values of the population migration rate from the blood to lymphoid tissues \(m_{BL}\) recurrence interval \(\tau \), mean neighborhoods \({\bar{n}}\), mean connections \({\bar{e}}\), and mean lymphoid tissue volume size \({\bar{C}}\) must be specified. We set \(m_{BL}=1000\), \(\tau =1.0\), \({\bar{n}}=5\), \({\bar{e}}=10\), \({\bar{C}}=10\) (mm: axial diameter), respectively. These numerical values are written in Table 1 in Appendix and captions of Figs. 3, 4, 5, and 7. In silico lymphoid tissue network is then generated under the rules described in Sects. 2.3 and 2.4.
Rights and permissions
About this article
Cite this article
Nakaoka, S., Iwami, S. & Sato, K. Dynamics of HIV infection in lymphoid tissue network. J. Math. Biol. 72, 909–938 (2016). https://doi.org/10.1007/s00285-015-0940-x
Received:
Revised:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00285-015-0940-x
Keywords
- Mathematical Modeling
- Numerical computation
- HIV infection
- Lymphoid tissue network
- The basic reproduction number
- Combinational drug treatment