Mathematical analysis and dynamic active subspaces for a long term model of HIV

  • Received: 04 November 2015 Accepted: 23 October 2016 Published: 01 June 2017
  • MSC : Primary: 92B99; Secondary: 34A34

  • Recently, a long-term model of HIV infection dynamics [8] was developed to describe the entire time course of the disease. It consists of a large system of ODEs with many parameters, and is expensive to simulate. In the current paper, this model is analyzed by determining all infection-free steady states and studying the local stability properties of the unique biologically-relevant equilibrium. Active subspace methods are then used to perform a global sensitivity analysis and study the dependence of an infected individual's T-cell count on the parameter space. Building on these results, a global-in-time approximation of the T-cell count is created by constructing dynamic active subspaces and reduced order models are generated, thereby allowing for inexpensive computation.

    Citation: Tyson Loudon, Stephen Pankavich. Mathematical analysis and dynamic active subspaces for a long term model of HIV[J]. Mathematical Biosciences and Engineering, 2017, 14(3): 709-733. doi: 10.3934/mbe.2017040

    Related Papers:

    [1] Miguel Atencia, Esther García-Garaluz, Gonzalo Joya . The ratio of hidden HIV infection in Cuba. Mathematical Biosciences and Engineering, 2013, 10(4): 959-977. doi: 10.3934/mbe.2013.10.959
    [2] Yicang Zhou, Yiming Shao, Yuhua Ruan, Jianqing Xu, Zhien Ma, Changlin Mei, Jianhong Wu . Modeling and prediction of HIV in China: transmission rates structured by infection ages. Mathematical Biosciences and Engineering, 2008, 5(2): 403-418. doi: 10.3934/mbe.2008.5.403
    [3] Titin Siswantining, Alhadi Bustamam, Devvi Sarwinda, Saskya Mary Soemartojo, Moh. Abdul Latief, Elke Annisa Octaria, Anggrainy Togi Marito Siregar, Oon Septa, Herley Shaori Al-Ash, Noval Saputra . Triclustering method for finding biomarkers in human immunodeficiency virus-1 gene expression data. Mathematical Biosciences and Engineering, 2022, 19(7): 6743-6763. doi: 10.3934/mbe.2022318
    [4] S. J. Gutowska, K. A. Hoffman, K. F. Gurski . The effect of PrEP uptake and adherence on the spread of HIV in the presence of casual and long-term partnerships. Mathematical Biosciences and Engineering, 2022, 19(12): 11903-11934. doi: 10.3934/mbe.2022555
    [5] A. M. Elaiw, N. H. AlShamrani . Stability of HTLV/HIV dual infection model with mitosis and latency. Mathematical Biosciences and Engineering, 2021, 18(2): 1077-1120. doi: 10.3934/mbe.2021059
    [6] Divine Wanduku . A nonlinear multi-population behavioral model to assess the roles of education campaigns, random supply of aids, and delayed ART treatment in HIV/AIDS epidemics. Mathematical Biosciences and Engineering, 2020, 17(6): 6791-6837. doi: 10.3934/mbe.2020354
    [7] Lih-Ing W. Roeger, Z. Feng, Carlos Castillo-Chávez . Modeling TB and HIV co-infections. Mathematical Biosciences and Engineering, 2009, 6(4): 815-837. doi: 10.3934/mbe.2009.6.815
    [8] Arni S. R. Srinivasa Rao, Kurien Thomas, Kurapati Sudhakar, Philip K. Maini . HIV/AIDS epidemic in India and predicting the impact of the national response: Mathematical modeling and analysis. Mathematical Biosciences and Engineering, 2009, 6(4): 779-813. doi: 10.3934/mbe.2009.6.779
    [9] Sophia Y. Rong, Ting Guo, J. Tyler Smith, Xia Wang . The role of cell-to-cell transmission in HIV infection: insights from a mathematical modeling approach. Mathematical Biosciences and Engineering, 2023, 20(7): 12093-12117. doi: 10.3934/mbe.2023538
    [10] Andrew Omame, Sarafa A. Iyaniwura, Qing Han, Adeniyi Ebenezer, Nicola L. Bragazzi, Xiaoying Wang, Woldegebriel A. Woldegerima, Jude D. Kong . Dynamics of Mpox in an HIV endemic community: A mathematical modelling approach. Mathematical Biosciences and Engineering, 2025, 22(2): 225-259. doi: 10.3934/mbe.2025010
  • Recently, a long-term model of HIV infection dynamics [8] was developed to describe the entire time course of the disease. It consists of a large system of ODEs with many parameters, and is expensive to simulate. In the current paper, this model is analyzed by determining all infection-free steady states and studying the local stability properties of the unique biologically-relevant equilibrium. Active subspace methods are then used to perform a global sensitivity analysis and study the dependence of an infected individual's T-cell count on the parameter space. Building on these results, a global-in-time approximation of the T-cell count is created by constructing dynamic active subspaces and reduced order models are generated, thereby allowing for inexpensive computation.


    1. Introduction

    The Human Immunodeficiency Virus (HIV) disables many components of the body's immune system and, without antiretroviral treatment, leads to the onset of Acquired Immune Deficiency Syndrome (AIDS). Despite the vast amount of study devoted to understanding viral pathogenesis and developing new therapeutics, no procedure or medication currently exists to reliably eliminate the virus from a host. However, new advances in long-term treatment strategies and insight into disease dynamics have stemmed from mathematical and computational modeling approaches, in addition to clinical experimentation.

    A variety of mathematical models have been proposed to describe HIV infection and disease dynamics [1,12,13,14,15,16,17,19,23,24,26]. Unfortunately, using a model to capture the entire time course of infection within the body can be troublesome as many oversimplify the biological dynamics of the disease in an effort to gain mathematical and rudimentary biological insight, and fail to capture all stages of infection. The majority of models accurately capture only the first stage(s) of infection [1,12,13,16,17,18,21] with the T-cell count and viral load asymptotically approaching a nonzero limit -the latter referred to as the viral set point. One recent description has been able to provide a holistic understanding of disease dynamics by accurately capturing all three stages of infection. This model, proposed in [8], is comprised of a system of seven nonlinear autonomous differential equations that are fully coupled and augmented by twenty-seven distinct parameters. In the current study we investigate the dynamical properties of the model and establish a result concerning its large-time asymptotic behavior. Further, we analyze this model, utilizing mathematical and statistical methods to elucidate the contribution of the parameters to an infected individual's T-cell count. In particular, we conduct a global sensitivity analysis in order to determine which of the twenty-seven parameters possess the greatest influence over the time course of the T-cell population. With this information, a reduced model (in the parameter space) can be constructed and solutions can be approximated as a function of time.

    The paper proceeds as follows. Within the next section, all infection-free steady states of the model are determined and the local asymptotic stability properties of the biologically-relevant equilibrium are studied. The associated theorems have interesting implications for the model's predictive nature, especially upon the introduction of antiretroviral therapy. In Section 3, we utilize active subspace methods to perform a global sensitivity analysis of the model with respect to its parameter space and further investigate the dependence of an infected individual's T-cell count on system parameters. With this information, a global-in-time approximation of the T-cell count is constructed using dynamic active subspaces. Additionally, various reduced models are constructed to represent different stages of the disease, and these are discussed in detail. In general, active subspace methods are a useful tool to perform global sensitivity analysis of a given parameter space, provide a clear picture of the most important activity in a model arising from parameter variation, construct dimensionally-reduced approximations to complex, dynamical models, and execute inexpensive numerical approximations from models that require computationally-intensive simulation. This is of particular importance in the current context as parameters for the original long-term HIV model are individual-dependent, and therefore must be determined for each new patient. Hence, the construction of a significantly less expensive computational approximation will allow one to utilize the model for large sets of patient data. Finally, we note that all of the MATLAB scripts and functions used to generate our results are provided free, open source, and available to the public at

    http://inside.mines.edu/~pankavic/activeHIV

    while an interactive iPython notebook is freely accessible at

    http://nbviewer.jupyter.org/github/paulcon/as-data-sets/blob/master/HIV/HIV.ipynb.

    Additional examples of biological and physical models for which active subspace methods have been implemented are also accessible at this website.


    2. Model description, parameters, and analysis

    To begin, we consider the following long term model of HIV disease dynamics within a host, as recently formulated in [8]:

    dTdt=s1+p1C1+VTVδ1T(K1V+K2MI)TdTIdt=ψ(K1V+K2MI)T+α1TLδ2TIK3TICTLdTLdt=(1ψ)(K1V+K2MI)Tα1TLδ3TLdMdt=s2+K4MVK5MVδ4MdMIdt=K5MVδ5MIK6MICTLdCTLdt=s3+(K7TI+K8MI)CTLδ6CTLdVdt=K9TI+K10MIK11TV(K12+K13)MVδ7V.} (1)

    The population of CD4+ T-cells, denoted here by T(t), is one of the most critical components in determining the body's response to HIV infection, and the first equation represents its time evolution. In (1) the T-cell population is increased by a standard source, s1, which represents the constant supply rate of immunocompetent T-cells from the thymus, and a nonlinear generation term p1C1+VTV, which accounts for the appearance of new T-cells due to the immune system's response to the infection. In contrast, because T-cells have a finite lifespan, a natural death term δ1T is also included. The last two terms in the first equation model the infection of T-cells by either virions, the population of which is denoted by V(t), or infected macrophages, MI(t), at rates K1 and K2, respectively. The latter term is introduced in (1) as studies have shown that infected macrophages likely play a vital role in the progression of infection [7,10] by producing large amounts of virus even after the T-cell population has been depleted.

    The second equation describes changes within the actively-infected T-cell population (TI), the first term of which represents the addition to this compartment due to new infections. However, not all interactions between virions and healthy T-cells produce actively-infected T-cells. Some proportion ψ of these new infections contribute to the TI population, while the remainder (1ψ) lack necessary host factors, resulting in the absence of viral protein expression, and become latently-infected T-cells (TL), contributing to the third equation. The infected T-cell population adds, in addition to the number created by the virions and macrophages, a supply of newly-activated latent T-cells, at a per-capita rate α1. In turn, this portion of the TL population is lost within the third equation. Infected T-cells are lost due to natural death, which is represented by δ2TI, while a more interesting term, K3TICTL, represents the loss of infected T-cells due to cytotoxic lymphocytes (CTL), one of the many attacker cells the immune system employs. Latently-infected T-cells also possess a natural death rate denoted by δ3.

    Similar to T-cells, the macrophage population (M) possess a natural birth rate, s2, as well as a natural death rate, δ4. As discussed in [17], macrophages divide and become more aggressive in seeking out pathogens once they are alerted of a viral presence by CD4+ T-cells. Hence, the macrophage population is increased in response to HIV infection at a rate of K4. Macrophages also attempt to eliminate virions, but may also become infected, adding to the infected macrophage population (MI). Infected macrophages naturally expire at a certain rate δ5, and are also destroyed by cytotoxic lymphocytes at a rate K6. Once infected, macrophages produce virions at a rate given by K10. Infected macrophages may also infect healthy T-cells, and the rate at which this occurs is denoted by K2.

    The main defender of the body against infected cells is the cytotoxic lymphocyte, the population of which is denoted by CTL. These cells seek to destroy renegade T-cells and macrophages that have been infected and altered by HIV. As for the other components of the immune system, we assume that cytotoxic lymphocytes are produced at a constant rate s3 by the bone marrow. Additionally, new lymphocytes are produced proportionate to the other aspects of the body's immune response. Hence, the term (K7TI+K8MI)CTL, occurs within the sixth equation to represent the recruitment of new lymphocytes.

    Lastly, the growth of the virion population V depends on a variety of parameters. New virions are continually produced by infected T-cells and infected macrophages at rates K9 and K10, respectively. In addition, the virion population is decreased through a variety of means. Due to infection, virions are lost at the rates K11 and K12, proportional to the interaction of T-cells and macrophages with virions, respectively. Also, macrophages, responding to antibodies produced by a host, ingest and destroy virions at a rate K13. Finally, virus particles are cleared from the host by other factors, including the innate immune response, at a rate δ7.

    Note that all parameter values in (1) are positive. Typical values and ranges for the parameters taken from [8] can be found within Table 1. The time course of infection predicted by this model with the established parameter values was shown to agree well with clinical data from [5,6,22], and representative simulations of (1) with realistic parameter values are displayed in Figure 1.

    Table 1. Parameter values and ranges.
    Parameter Value Range Value taken from: Units
    s1 10 5 -36 [13] mm3d1
    s2 0.15 0.03 -0.15 [13] mm3d1
    s3 5 - [8] mm3d1
    p1 0.2 0.01 -0.5 [8] d1
    C1 55.6 1 -188 [8] mm3
    K1 3.87 x 103 108 -102 [8] mm3d1
    K2 106 106 [13] mm3d1
    K3 4.5 x 104 104 -1 [8] mm3d1
    K4 7.45 x 104 - [8] mm3d1
    K5 5.22 x 104 4.7 x 109 -103 [8] mm3d1
    K6 3 x 106 - [8] mm3d1
    K7 3.3 x 104 106 -103 [8] mm3d1
    K8 6 x 109 - [8] mm3d1
    K9 0.537 0.24 -500 [8] d1
    K10 0.285 0.005 -300 [8] d1
    K11 7.79 x 106 - [8] mm3d1
    K12 106 - [8] mm3d1
    K13 4 x 105 - [8] mm3d1
    δ1 0.01 0.01 -0.02 [8] d1
    δ2 0.28 0.24 -0.7 [8] d1
    δ3 0.05 0.02 -0.069 [8] d1
    δ4 0.005 0.005 [13] d1
    δ5 0.005 0.005 [13] d1
    δ6 0.015 0.015 -0.05 [27] d1
    δ7 2.39 2.39 -13 [13] d1
    α1 3 x 104 - [8] d1
    ψ 0.97 0.93 -0.98 [8] -
     | Show Table
    DownLoad: CSV
    Figure 1. Ten simulations of (1) with representative parameter values.

    2.1. Infection-free steady states

    Though the system (1) possesses a large number of steady states -the authors have discovered at least ten using standard parameter values and a computational root finder -one is often most interested in understanding the dynamical properties of the disease-free equilibrium. In this section, we identify such equilibria and investigate the stability of the biologically-relevant state. Our first result demonstrates that only one such equilibrium state exists when all populations of (1) are positive.

    Theorem 2.1. The model (1) possesses exactly two virus-free (i.e. V0) steady states. One of these states, namely

    E:=(s1δ1ωK2K9,ωK10,ωK10ξK6(α1+δ3ψ),s2δ4,ωK9,δ5K6,0)

    achieves negative values, where

    ω=s3K6+δ5δ6δ5(K7K10K8K9)andξ=(1ψ)(δ2K6δ5K3).

    The only nonnegative (i.e. biologically relevant) steady state of (1) satisfying V0 is

    ENI:=(s1δ1,0,0,s2δ4,0,s3δ6,0).

    Hence, the only guarantee of viral clearance as t occurs when actively and latently infected populations are also eradicated, resulting in healthy T-cell and macrophage populations tending asymptotically to background values. The proof of Theorem 2.1 is contained in Appendix A. Utilizing standard parameter values for this model from Table 1, we find the following equilibrium populations for E:

    T=1010.39mm3,TI=54.57mm3TL=15.77mm3,M=30mm3MI=102.83mm3,CTL=1666.67mm3,V=0mm3

    Since the parameter values in (1) are positive, the steady state E given in Theorem 2.1 must have a negative cytotoxic T-lymphocyte population, namely δ5K6. So, under no parameter regime will E be biologically relevant. With the unique infection-free steady state identified, we turn to its stability properties.


    2.2. Stability analysis

    Next, we provide necessary and sufficient conditions which guarantee the local asymptotic stability of the disease-free equilibrium ENI.

    Theorem 2.2. The equilibrium state ENI is locally asymptotically stable if and only if R01, where

    R0=max{R1,R2,R3}

    and

    R1=K1K9δ2K11,R2=K5K10(K12+K13)δ5
    R3=K2K5K9s1s2δ1δ2δ4δ5δ7+δ4δ5K1K9s1+δ1δ2K5K10s2.

    The proof of Theorem 2.2 is also contained within Appendix A. Computing the basic reproduction number of Theorem 2.2 by using the standard parameter values given in Table 1, we find that R0=R1=953>>1. Hence, as expected, the non-infective steady state ENI is not locally asymptotically stable.

    Introducing antiretroviral therapy, or ART, into the system provides additional insight into this result. Two specific classes of ART drugs, namely Reverse Transcriptase Inhibitors (RTIs) and Protease Inhibitors (PIs), serve to reduce the amount of new virus produced by either reducing the ability of virions to replicate through reverse transcription or disabling the capability of newly-produced virions to mature, thereby rendering them uninfective. The efficacies of these classes of drug, denoted ϵRTI,ϵPI[0,1], can be incorporated using the transformations

    K1K1(1ϵRTI),  K5K5(1ϵRTI),  K9K9(1ϵPI),  K10K10(1ϵPI).

    As these constants appear within the stability result and decrease each of the ratios R1, R2, and R3, it follows that a large enough efficacy will force the system to tend towards ENI as t. Therefore, one merely needs to perform this transformation on the result of Theorem 2.2 in order to determine what efficacy is needed to guarantee viral clearance. Using the previously-determined parameter values and defining ϵ[0,1] by

    1ϵ:=(1ϵRTI)(1ϵPI),

    we find ϵ>0.998 in order to force R0<1. Thus, a combined drug efficacy greater than 99.8% would be needed to asymptotically drive the system to clearance. Of course, lesser drug efficacies could still give rise to viral loads that are effectively negligible rather than tending to zero, and therefore correspond to viral clearance.

    Notice that the asymptotic stability result in Theorem 2.2 depends only upon a relatively few number (15 of 27) of the parameters. Therefore, the majority of the pertinent dynamics takes place on a lower-dimensional subspace of the entire parameter space. Hence, in the next section we explore a dynamic tool to better understand the contribution of the parameter space to populations within the model. In particular, this will lead to reducing the dimension of the parameter space with minimal loss of information using an active subspace decomposition.


    3. Dynamic active subspaces, sensitivity, and reduced models

    In this section we will use active subspace methods to approximate the T-cell count at a specific time given the parameter values in (1). An active subspace is a low-dimensional linear subspace of the set of parameters, in which input perturbations along these directions alter the model's predictions more, on average, than perturbations which are orthogonal to the subspace. These subspaces allow for a global measurement of sensitivity of output variables with respect to parameters, and often the construction of reduced-order models that greatly decrease the dimension of the parameter space.


    3.1. Active subspace methods

    The general structure of an active subspace decomposition begins by letting mN be given and defining the space X=[1,1]m. Also given is a differentiable function f:XR and an associated probability density ρ:XR+ satisfying

    ρ(x) dx=1.

    Here, the space X represents a normalized set of parameter values. With these quantities in place, consider the matrix C defined by

    C=(xf)(xf)Tρdx. (2)

    For any smooth f, the matrix C represents an average derivative functional which weights input values according to the density ρ. In general terms, f(x) represents the quantity of interest in a given model, while gradients of f are taken with respect to normalized model parameters xX, and ρ(x) is the probability density associated to the values of these parameters. We note here that a single normalized parameter is a random variable taking values in [1,1], which when appropriately scaled represents a parameter in the original model (1). Since the dimension of the parameter space in this model is 27, we take m=27 throughout. The matrix C is the average of the outer product of the gradient of f with itself and has some useful properties that will allow us to deduce information about how f is altered by perturbations in its arguments.

    Considering each entry of the matrix

    Cij=fxifxjρ(x)dx

    we note that C is symmetric, and thus permits the spectral eigendecompostion

    C=WΛWT,whereΛ=diag(λ1,,λm),λ1λm0. (3)

    Here, W is an orthogonal matrix whose columns wi,(i=1,,m) are the orthonormal eigenvectors of C.

    From (3) we can further solve for the eigenvalues of C, which are given by

    λi=((xf)Twi)2ρ(x)dx,i=1,,m. (4)

    From (4) we see that the eigenvalues of the C matrix are the mean squared directional derivatives of f, in the direction of the corresponding eigenvector. Thus, the eigenvalues of C provide useful information about the quantity of interest. For instance, if a particular eigenvalue is small then (4) tells us that, on average, f does not change significantly in the direction of the corresponding eigenvector. Conversely, if the eigenvalue under consideration is large, then we may deduce that f changes considerably in the direction of the corresponding eigenvector. Therefore, it will be of interest to further investigate the behavior of the function in this direction.

    Once the eigendecomposition (3) has been determined, the eigenvalues and eigenvectors can be separated in the following way:

    Λ=[Λ100Λ2],W=[W1W2]. (5)

    where Λ1 contains the "large" eigenvalues of C, Λ2 contains the "small" eigenvalues, and Wk contains the eigenvectors associated with each Λk, for k=1,2. An easy way to differentiate between the "large" and "small" eigenvalues is to list them on a log plot from greatest to least and determine a spectral gap. This gap will correspond to differences of at least an order of magnitude, and thus allow one to compartmentalize large eigenvalues within Λ1 and the remaining smaller eigenvalues in Λ2. A more systematic method of choosing how many eigenvalues to store within Λ1 will be presented in Section 3.2.

    With the decomposition (5), we can represent any element x of the parameter space by

    x=WWTIx=W1WT1xy+W2WT2xz=W1y+W2z. (6)

    Thus, evaluating the quantity of interest at x is equivalent to doing so at the point W1y+W2z, i.e.

    f(x)=f(W1y+W2z).

    By the definition of W1 and W2 it's clear that small perturbations in z will not, on average, alter the values of f. However, small perturbations in y will, on average, change f significantly. For this reason we define the range of W1 to be the active subspace of the model and the range of W2 to be the corresponding inactive subspace. The linear combinations that generate these subspaces will then represent the contributions of differing parameters in the model and describe the sensitivity of the quantity of interest with respect to parameter variations.

    In general, the eigenvalues and eigenvectors of C defined by (2) can be well-approximated, using a random sampling algorithm. We will briefly outline the method, but full details can be found in [2] (Algorithm 3.1) and [3]. The algorithm can be described concisely as follows:

    1. Draw N samples {xj}Nj=1 independently according to the density ρ.

    2. For each parameter sample xj, approximate the gradient xfj=xf(xj) using the finite difference, i.e.

    xif(xj)f(xj+hi)f(xj)|hi|

    where

    (hi)k={δif i=k0if ik.

    represents a vector perturbation from the sampled parameter values and δ>0 can be taken arbitrarily small.

    3. Approximate the matrix C by

    CˆC=1NNj=1(xfj)(xfj)T

    4. Compute the eigendecompositions ˆC=ˆWˆΛˆWT.

    We note that the last step is equivalent to computing the singular value decomposition of the matrix

    1N[xf1xfN]=ˆWˆΛˆV, (7)

    where it can be shown that the singular values are the square roots of the eigenvalues of ˆC and the left singular vectors are the eigenvectors of ˆC. The singular value decomposition method of approximating ˆC was developed first in [25].


    3.2. An illustrative example -approximating T(1700)

    Next, we will demonstrate the active subspace method at one specific point in time by applying the aforementioned algorithm to the HIV model (1). A suitable quantity of interest would likely involve the T-cell or virus population evaluated at a fixed time. For this example we select T(1700), the T-cell count 1700 days after initial infection. This quantity was chosen because the T-cell count is a strong indicator of a patient's overall health and the most important factor in the progression of HIV infection. The time of 1700 days after initial infection was chosen because regardless of parameter values, preliminary simulations have shown that the patient's T-cell count within (1) will not have decreased to zero by this time. However, if a later time is chosen, the patient's T-cell count could vanish before the final time is reached.

    In order to compute C, one must be able to construct the gradient of the quantity of interest with respect to the normalized parameter space as required in Step 3 of the algorithm. Since we do not have an explicit representation for T(1700) as a function of system parameters, we cannot explicitly compute gradients. Instead we approximate them using the aforementioned finite difference scheme with a step size of δ=106. Additionally, each sample is chosen so that every normalized parameter is uniformly distributed between 1 and 1, i.e. xj(U[1,1])27. In order to map this normalized parameter space onto the biologically-relevant range of parameter values, we use the linear mapping

    p=12(diag(xuxl)xj+(xu+xl)), (8)

    for each of the random samples xj(U[1,1])27, where xu and xl are vectors containing the upper and lower bounds on the parameters, respectively. Thus, the resulting vector p represents the actual parameter values input within the model. Next, a stiff differential equations solver (MATLAB's ode23s function) is used to compute the T-cell count after 1700 days. Then, each of the 27 parameters is perturbed by δ=106, and again the corresponding T-cell count after 1700 days is computed. With these two values, the aforementioned finite difference approximation is used to calculate the gradient of f(x)=T(1700;x) with respect to the normalized model parameters. Regarding the upper and lower limits, xu and xl are taken to be 2.5% above and below the standard values given in Table 1. Figure 2 displays the approximate eigenvalues of the corresponding C matrix. Clearly a spectral gap exists between the first and second eigenvalues. In order to automatically determine the optimal decomposition of Λ we can use the relative measure of separation given by

    Figure 2. Approximation of eigenvalues of C using 1000 random samples.
    ˆλk=λkλk+1λ1,k=1,2,,26. (9)

    Then, the dimension of the active subspace, i.e. the number of eigenvalues stored within Λ1, will be given by

    dim=argmaxk=1,,26ˆλk. (10)

    While the index of the largest value of ˆλk describes the location of the largest spectral gap, it is often convenient to consider only the first two values ˆλ1 and ˆλ2. Doing so limits the dimension of the active subspace to one and two respectively, which allows for easy visualization of the quantity of interest as a function of the active subspace and allows one to fit a curve or surface to the data. Plotting the values of ˆλk results in Figure 4.

    Figure 4. Measure of separation for the eigenvalues of C.

    Clearly, with this measure of separation, the optimal choice for the dimension of the active subspace is merely one. Consequently, we store λ1 in the matrix Λ1 and the remaining eigenvalues λi, i=2,,27, along the diagonal of Λ2. The active subspace is then generated by linear combinations of the entries of w, the first eigenvector. Figure 3 displays the eigenvector corresponding to the maximal eigenvalue shown in Figure 2, and it can be seen that three parameters possess associated weights greater than 0.3. These are the 9th, 22nd, and 25th parameters as ordered in Table 1. From the table these parameters can be identified as K4, δ4, and δ7, which represent the increase in macrophage population due to the immune system, the death rate of the macrophage population, and the death rate of the virus population, respectively. Hence, small perturbations in these parameters will significantly alter the value of T(1700), as they are the most heavily weighted. Contrastingly, changes within the remaining parameters, whose weights are near zero, will not have an appreciable affect on T(1700).

    Figure 3. Approximation of the 1st eigenvector of C using 1000 random samples. This is referred to as the first active variable vector and denoted by w.

    Plotting T(1700) along the active subspace results in Figure 5. Here, the horizontal axis is represented by values of the first active variable y=xw, which represents a linear combination of the normalized parameters x with weights given by entries of the first active variable vector w. Adopting the terminology of [2], we will refer to plots of the quantity of interest along the active subspace as sufficient summary plots. Figure 5 (left) shows a clear trend, namely that T(1700;y) is a decreasing function of the active variable. Additionally, the data points within the sufficient summary plot can be fit to a particular function. In this case, a four parameter arctangent function is fit to the data to approximate the quantity of interest. This was performed using the MATLAB function lsqcurvefit which minimizes the residual (in the least squares sense) of the difference in the data and the approximation. Hence, we have determined that the data is best fit by

    Figure 5. Sufficient summary plot after 1700 days (left). Approximation to the T-cell count after 1700 days (right).
    T(1700;y)=79.2532492.5680 tan1(0.8933y1.9069), (11)

    and the resulting curve can be found in the plot on the right side of Figure 5. In order to test the accuracy of the nonlinear approximation, 100 simulations were run and the relative error computed. The results are displayed in Figure 6 and indicate that the approximation is typically within 5% from the value computed by computationally solving (1) to determine T(1700). Stated another way, over 90% of the test simulations returned approximations whose error was 5% or less.

    Figure 6. Relative errors in the approximation of T (1700).

    Of course, not every eigenvalue decomposition will necessarily result in a one-dimensional active subspace of the model parameters. For instance, should one wish to approximate the T-cell count 2000 days after initial infection, the same order reduction method can be employed and the eigenvalues of the C matrix can be computed. The plot of the eigenvalues of C can be seen in Figure 7 (left). In this case the largest spectral gap exists between the second and third eigenvalues rather than the first and second, and therefore a more descriptive choice for the dimension of the active subspace is two instead of one.

    Figure 7. Eigenvalues of the matrix C after 2000 days (left). Dimension of the active subspace for each time (right).

    The two weight vectors corresponding to these dominant eigenvalues are then computed and sufficient summary plots are obtained. Figure 8 shows the one and two dimensional sufficient summary plots, respectively, for the approximation of T(2000). Considering the two dimensional sufficient summary plot in Figure 8 (right), it can be seen that very little variation occurs within the T-cell count in the vertical direction, which represents the second active variable. All of the variation in the T-cell count appears to transpire mostly in the horizontal direction. Also, in Figure 8 (left) we can see that the one dimensional sufficient summary plot clearly shows a distinct trend. For these reasons, the one-dimensional active subspace representation provides a sufficient description of the dynamics in the parameter space, while the second active variable doesn't appear to contribute as greatly. Hence, even though the largest spectral gap appears between the second and third eigenvalues, a single linear combination of the parameters captures the overwhelming majority of the dynamics, and only the one-dimensional representation is utilized. Figure 7 (right) displays the dimension of the active subspace at each time using the measure of separation and (10).

    Figure 8. Sufficient summary plots after 2000 days, displaying the one-dimensional (left) and two-dimensional (right) active subspace representations.

    Since this algorithm can be used to approximate the T-cell count at any fixed time, we can further compute active subspaces to reduce complexity in the parameter space and obtain a simple time course for T(t). We refer to such a model as a dynamic active subspace approximation. Within the next section we will demonstrate precisely how this reduced description can be constructed.


    3.3. Dynamic active subspaces

    Now that the approximation for T(1700) has been computed, this method can be repeated for a predetermined discretization of the time domain in order to create a global-in-time approximation for the T-cell count, say T(t;y). First, it is necessary to choose the mesh on which T(t) will be approximated. In this case the time interval [0,3400] is divided into 86 non-uniform subintervals and the active subspace decomposition, along with the associated eigenvalues and eigenvectors, is used to determine a functional approximation at each time step. Next, one must orient the eigenvectors to point in approximately the same direction so that they transition smoothly from one time step to the next. By this we mean that the magnitude of the components of the consecutive weight vectors differ only slightly, but because the normalized eigenvector decomposition is only unique up to a sign, weight vectors at different time steps must be oriented so that they do not point in opposing directions. Upon correctly orienting the weight vectors, nonlinear curve-fitting is employed to compute the functions that best fit the associated sufficient summary plots.

    Once the approximations have been computed at each time step, the final task is to dynamically assemble them using linear basis functions, namely

    T(t;x)86i=1Ti(xw(t))ϕi(t) (12)

    where w(t)=86i=1wiϕi(t) for all t[0,3400] is a linear interpolation of the first active variable vector, wi, within the ith time interval, y(t)=xw(t) is the corresponding active variable, Ti(y) is the approximation to the T-cell count within the ith time interval, and ϕi(t) is the hat function on the interval [ti1,ti+1] given by

    ϕi(t)={tti1titi1,ti1ttiti+1tti+1ti,titti+10,t[ti1,ti+1].

    Recall that the original parameter values p and the normalized parameters x can be interchanged using (8). Here we use the notation T(t;x) instead of T(t;y) because within the global-in-time approximation, the active variable vector, w and hence the active variable y, changes within differing time intervals. Of course, a higher-order approximation can be obtained by expanding ϕi(t) in a polynomial basis of greater degree.

    The results of computing the active subspace and sufficient summary plot for 18 of the 86 time steps can be found in Appendix B. With this we see that trends in the sufficient summary plot transition smoothly from one time step to the next and at each fixed time t the T-cell count can be represented by one of three distinct functional forms. In particular, curves that are fit to the data generated within sufficient summary plots transition from (a) linear to (b) arctangent trends at or around 5565 days after initial transmission, while for times greater than or equal to 1800 days, the trend resembles (c) an arctangent function multiplied by a heaviside step function. The last of these aproximations occurs during later periods of the time course when the T-cell count may vanish for certain values of the active variable. For example, Figure 9 displays a sufficient summary plot of the T-cell count at 2600 days after infection occurs. For large values of the first active variable, say y0.5, we see that T(2600;y)=0 since a proportion of sample runs result in a patient's T-cell count tending to zero prior to 2600 days. Contrastingly, for y<0.5, the T-cell count trend resembles an arctangent function. These three specific functional approximations and their transitions comprehensibly display the biological aspects of all stages of infection dynamics. Before discussing this further, a minor discussion of HIV disease pathogenesis is needed.

    Figure 9. Sufficient summary plot after 2600 days using 1000 trials (left). Same plot with function approximation (right).

    The time-course of HIV infection is characterized by three distinct stages: acute (or primary) infection, chronic infection (also referred to as clinical latency), and the transition to AIDS. The first of these phases takes place within 1020 weeks of initial introduction of the virus within a host and is characterized by a rapid fluctuation in the T-cell and virion population. With respect to the T-cell population there is initially a rapid decrease from the introduction of the virus, and then a rapid rebound arising from the body's immune response. Symptoms during this phase of the infection include fever, swollen glands, fatigue, rash, and sore throat. The next stage, chronic infection, ranges from a number of years to over a decade without treatment. During the chronic period the T-cell and virion populations remain at relatively constant levels, with the T-cell population decreasing at a particularly slow rate and the number of virions increasing steadily. During the last stage of infection, the transition to AIDS occurs as the T-cell population reaches a density lower than 200 cells per mm3. The progression to AIDS is typically associated with a sharp decrease in the T-cell population within a year or so.

    From the sufficient summary plots in Figure 10 and Appendix B, the three distinct stages of the disease become clear. The three functional forms arising from sufficient summary plots precisely separate the three distinct stages of HIV disease progression within an infected individual -the initial arctangent function representing the acute phase, followed by a slow linear decline that denotes the asymptomatic or chronic phase, and finally the heaviside arctangent function detailing the decline of the T-cell count as a patient develops AIDS. This description of the T-cell population provides a detailed visual account of the three stages and the transitions between them. For clarity, a representative graph of each trend in the data is provided in Figure 10. Therefore, the abrupt changes in the T-cell count arising within the one-dimensional active subspace of parameters completely categorize the stages of disease pathogenesis.

    Figure 10. Sufficient summary plots representing the three stages of infection -Acute (left), Chronic (center), AIDS (right).

    In addition to separating these stages of the disease, the active subspace method allows for the creation of three different types of approximate models. First, a visual representation of the quantity of interest -in this case, the T-cell count -as a function of the most pertinent linear combination of parameters within the original model is provided. Additionally, the method allows for the construction of an explicit, analytic model by combining nonlinear function approximations such as (11) over an interval of time. In this direction a second approximation method is available instead of utilizing basis functions as in (12). Namely, one may prescribe the functional form during a particular stage of the disease and fit time-dependent coefficients to the transitions within sufficient summary plots. For instance, one may express the arctangent approximation over the timespan t[0,55] by

    T(t;y)=a(t)+b(t) tan1(c(t)yd(t))

    for well-behaved functions a,b,c, and d that are fit to the changes in the data. This type of explicit approximation would be easiest to implement during the asymptomatic phase for which

    T(t;y)=m(t)y+b(t)

    and the slope and T-intercept functions, m(t) and b(t) respectively, are given by the piecewise-defined or smoothed functions approximated in Figure 11.

    Figure 11. Slope (left) and T-intercept (right) functions, m(t) and b(t), respectively for t ∈ [55,1300].

    Lastly, the nonlinear fits arising from the sufficient summary data can be easily stored and supply the basis for a low-cost computational approximation without the need to simulate the original model over long time periods. In contrast to simulating the full system of ODEs for each new set of parameter values, this computational model need only be precomputed once and can easily describe which of the parameters are most important and during which stages they significantly alter the biological quantities of interest. As parameters within the original HIV model (1) vary amongst patients, they must be recomputed for each individual, and the computational savings provided by the reduced model is vital.

    Utilizing the dynamical algorithm represented by (12), this computational model can be constructed and compared with a representative simulation of the full dynamical system. Figure 12 displays both the simulated T-cell count and its global-in-time approximation using dynamic active subspaces. We note that because the active subspace method is global with respect to the parameter space, the precise values of parameters in Table 1 do not necessarily influence the structure of the active subspace model, as long as a feasible range of parameter values is available.

    Figure 12. Global-in-time approximation of the T-cell count.

    In order to test the accuracy of the active subspace approximation to solutions of (1), 100 independent simulations were performed and the relative error was computed. Within these error calculations a uniformly-distributed random time was selected along with a uniformly-drawn selection of the parameter values within their respective ranges. The result is shown in Figure 13 , and demonstrates that for 96% of simulations the analytic approximation varied less than 5% from the value given by the stiff differential equation solver. Hence, one may conclude that the global-in-time active subspace model well-approximates solutions to the original system of ODEs given in (1).

    Figure 13. Relative error in the global approximation of the T-Student Version of MATLAB cell count.

    3.4. Dimension reduction in the parameter space

    As determined in the previous section, not all parameters are required to describe the behavior of solutions within the model and those that are needed may not be important during each stage. Therefore, it makes sense to investigate the construction of reduced models, namely those which eliminate the contributions of certain parameters. Though the active subspace method has reduced the dimension of the parameter space upon which the T-cell count depends, all of the parameter values are still needed in order to compute the approximate solution, as the reduced parameter space has been expressed merely as a linear combination of the original parameter values, i.e. y=xw. However, the weights (or coefficients) within this linear combination, given by w, should provide a clear method to reduce the original parameter space as well. For instance, assume that a parameter, say x3, possesses a corresponding weight entry w3 that is very small in comparison to the other weights. Then, variations in x3 will have little to no effect on the reduced parameter space y, and hence will not appreciably influence the output variable T(t;y). Thus, one can merely eliminate x3 within the model or set x30 rather than considering it as a parameter whose value could possibly vary depending upon the patient.

    While it would be most beneficial to determine those weight vectors that remain small throughout the entire course of infection and remove the corresponding parameters, this is an unlikely scenario as different parameters will typically influence different stages in some substantial manner. For this reason, a true global-in-time parameter reduction is highly unlikely. Indeed, the computed weight vectors for this study show that at most 3 of the 27 parameters -namely, K11, K12, and K13 -could be eliminated without introducing enormous variations in the behavior of the model during some stage of infection. Such a limited reduction in complexity would be marginally beneficial to the expense of computations. Instead, we may separate the dynamics into the three distinct stages of the disease and use parameter reduction to create new models for each stage separately. By eliminating the associated interactions from (1) we may derive a simpler system that still accurately predicts each of the three phases of infection.

    As an illustrative example, we consider the weight vectors arising within the first 40 days of initial infection and remove parameters whose weights remain below a fixed threshold throughout this interval of time. Using a relatively strict threshold, in this case 0.032, a number of parameters are removed from (1), namely s2, s3, K2, K4, K5, K6, K7, K8, K10, K11, K12, K13, δ3, δ4, δ5, δ6, and α1. With these parameters eliminated, the TL, M, MI, and CTL populations decouple from the remaining equations. This implies that latently-infected T-cells, macrophages (both healthy and infected), and cytotoxic lymphocytes do not play an important role in the early behavior of the disease, i.e. within the first five to six weeks of introduction of the virus within the body. Because these populations decouple, it's necessary to remove K3 (whose maximum weight during the first 40 days is 0.1256) and set ψ = 1, thereby yielding the reduced system for the acute stage, namely

    dTdt=s1+p1C1+VTVδ1TK1TVdTIdt=K1TVδ2TIdVdt=K9TIδ7V.} (13)

    This model is an augmented form of the well-known three-component model (see [11]) with an additional Michaelis-Menten term within the T-cell population, which accounts for the homeostatic proliferation of such cells upon depletion of this compartment due to interactions with virions. This model was recently analyzed in detail in [20] and found to possess many desirable properties, including bistable equilibria which explain the dependence of infection dynamics on the initial T-cell count and viral load, as well as the existence of a Hopf bifurcation which describe oscillations within the system. Using the fitted parameter values given in [20] for the model (13) and plotting against the full model (1) with the standard parameter values given in Table 1 results in Figure 14.

    Figure 14. Full HIV model versus reduced HIV model for the first 100 days. Parameter values within the reduced model are s1 = 10, p1 = 0.2, C1 = 55.6, δ1 = 0.01, K1 = 4.72 × 10−3, δ2 = 0.69, K9 = 5.37 × 10−1, and δ7 = 2.39.

    Hence, for the acute stage we have reduced (1) with 27 independent parameters, to (13), which features only 8 parameters. Following the same procedure for the other two stages would further reduce the dimension of the parameter space in the model (1) at later stages of the infection. These reduced models could then be used to analyze the disease separately within each distinct stage of infection and with much less computational cost. If a global model is still preferred in comparison to separating models by stage, one can perform this parameter reduction for all three stages subject to the constraint that values of T(t) within time-adjacent models be equal and a new global-in-time dynamical model arises by construction. Such a reduced, long-term model would be comparable to that of Figure 12.


    4. Conclusion

    The current study concerns an analysis of the system (1), which is one of the only mathematical models to accurately represent all three stages of HIV infection within a host. A unique, biologically-relevant virus-free equilibrium was shown to exist, and conditions were determined that guarantee the local asymptotic stability of this state. Then, using dynamic active subspaces, the system of seven ODEs with 27 parameters was approximated by algebraic and computational models while retaining the majority of pertinent information and system behavior. This method enabled the discovery of a dominant subspace of parameters, within which the T-cell count is most sensitive to perturbations, and allowed us to perform thorough parameter studies in this direction, rather than the complete 27-dimensional parameter space. In general, active subspace methods also provide for a deeper visualization of the dependence of solutions on the parameter space, as well as, an analytic model (12) to describe the quantity of interest and a computational model for which solutions are less expensive to construct. The efficacy of the model (12) was investigated by calculating the relative error compared to solutions of (1) solved with a stiff ODE solver.

    As with any study, a number of future questions arise from our investigation. First, as many parameter values within (1) are not specifically known and vary greatly in the literature, one would like to quantify the uncertainty inherent in choosing a particular value for each. While this can be investigated using the approximations established herein, it was not the focus of the study and more complex mathematical tools are needed to do so. Another direction to consider is the explicit quantification of error within the lower-dimensional approximations provided by active subspaces. Some preliminary results appear in [2], but exact bounds and convergence theorems for dynamic, rather than time-independent, active subspaces are currently unavailable within the literature. Finally, we note that these methods can be used for other in-host models of HIV (or other physical and biological models [4]) that possess high-dimensional parameter spaces. For instance, a recent refinement of (1) was proposed in [9], and appears to display less sensitivity to variations in parameter values. Hence, a similar technique would likely be useful to conduct a parameter study or construct a reduced-order model for this system.


    Appendix A. Proofs of theorems

    In the first appendix, we outline the proofs of the theorems stated in Section 2.

    Proof of Theorem 2.1. Beginning with (1), we search for steady states by assuming that all time derivatives are zero within the equations, and attempt to solve for the constant states (T,TI,TL,M,MI,CTL,V). Assuming V=0 within the system of ODEs provides a significant reduction in the complexity of the system. The M equation implies M=s2δ4. Using this within the equation for MI implies that either MI=0 or CTL=δ5K6.

    Consider the latter case first. Multiplying the TI equation by (1ψ) and the TL equation by ψ and adding gives

    0=(α1+ψδ3)TL+((1ψ)K3δ5K6(1ψ)δ2)TI. (14)

    Creating a linear system with (14), the CTL equation, and the V equation then solving for TI, TL, and MI gives

    TI=K10ω,TL=K10ξωK6(α1+δ3ψ),MI=K9ω

    where

    ω=s3K6+δ5δ6δ5(K7K10K8K9)andξ=(1ψ)(δ2K6δ5K3).

    Lastly, inserting the value of MI into the T equation and solving for T yields

    T=s1δ1ωK2K9

    Finally, consider the former case. Then, it follows from the equation for V that TI=0 as well. Collecting these terms in the CTL differential equation implies that CTL=s3δ6. The equations for TI and TL together imply TL=0, and finally, with the remaining populations determined, the first equation implies T=s1δ1. Hence, we find the steady state

    ENI:=(s1δ1,0,0,s2δ4,0,s3δ6,0).

    We note that assuming all parameter values are strictly positive implies that the only non-infective steady state of biological significance is ENI.

    Finally, we sketch the proof of the asymptotic stability result, which utilizes standard methods from the theory of dynamical systems (i.e. the Hartman-Grobman and Routh-Hurtwitz theorems) to determine the qualitative behavior of the ENI steady state from (1).

    Proof of Theorem 2.2. As the model is not positivity-preserving a simple technique like the next-generation method does not apply. Hence, we will utilize the Hartman-Grobman Theorem to arrive at the stated result and omit some of the more technical details. We first compute the Jacobian of (1) evaluated at the steady states ENI.

    J(ENI)=(δ1000K2s1δ10(p1c1K1)s1c1δ10δ2δ6+K3s3δ6a10pK2s1δ10pK1s1δ100a1δ30(p1)K2s1δ10(p1)K1s1δ1000δ400(K4K5)s2δ40000δ5δ6+K6s3δ60K5s2δ40K7s3δ600K8s3δ6δ600K900K100δ7K11s1δ1(K12+K13)s2δ4).

    From this, we can see that three eigenvalues are certainly real and negative

    λ1=δ1,λ2=δ4,λ3=δ6.

    The remaining four eigenvalues are more difficult to identify as they are determined by the quartic equation

    a4λ4+a3λ3+a2λ2+a1λ+a0=0

    where

    a4=δ1δ4δ26>0
    a3=α1δ1δ4δ26+δ4δ26K11s1+δ1δ26K12s2+δ1δ26K13s2+δ1δ4δ6K3s3+δ1δ4δ6K6s3+δ1δ2δ4δ26+δ1δ3δ4δ26+δ1δ4δ5δ26+δ1δ4δ7δ26>0

    and a0, a1, and a2 are given by much longer expressions and not necessarily positive.

    Instead, we must impose conditions on each term to guarantee positivity, which is needed for the roots of the quartic to possess negative real part by the Routh-Hurwitz criteria. In particular, the two negative terms in a2 are dominated by the remaining positive terms if and only if

    K1K9δ2K11
    K5K10(K12+K13)δ5.

    The same conditions imply the positivity of a1. For a0, the negative terms are dominated by positive terms if and only if the two previous conditions hold and

    K2K5K9s1s2δ1δ2δ4δ5δ7+δ4δ5K1K4s1+δ1δ2K5K10s2.

    The final inequalities of the Routh-Hurwitz criteria are also implied by these conditions. Hence, defining

    R1=K1K9δ2K11,R2=K5K10(K12+K13)δ5,R3=K2K5K9s1s2δ1δ2δ4δ5δ7+δ4δ5K1K9s1+δ1δ2K5K10s2,

    and

    R0=max{R1,R2,R3}

    we see that the equilibrium is locally asymptotically stable if and only if all three conditions are satisfied, and thus R01.


    Appendix B. Visual representation of three stages of infection

    In Figures 15, 16, and 17, we provide a more temporally refined representation of the three stages of infection displayed by the model.

    Figure 15. Sufficient summary plots throughout the Acute stage.
    Figure 16. Sufficient summary plots throughout the Chronic stage.
    Figure 17. Sufficient summary plots during the progression to AIDS.

    Acknowledgments

    The authors would like to thank Prof. Paul Constantine for helpful discussions and advice, and Ryan Howard for constructing an openly-available iPython notebook that details the work in our computational study.


    [1] [ D. Callaway,A. Perelson, HIV-1 infection and low steady state viral loads, Bull.Math.Biol., 64 (2002): 29-64.
    [2] [ P. Constantine, Active Subspaces: Emerging Ideas for Dimension Reduction in Parameter Studies SIAM, 2015.
    [3] [ P. Constantine and D. Gleich, Computing active subspaces with monte carlo, arXiv: 1408.0545
    [4] [ P. Constantine,B. Zaharatos,M. Campanelli, Discovering an active subspace in a single-diode solar cell model, Statistical Analysis and Data Mining: The ASA Data Science Journal, 8 (2015): 264-273.
    [5] [ A. S. Fauci,G. Pantaleo,S. Stanley, Immunopathogenic mechanisms of HIV infection, Annals of Internal Medicine, 124 (1996): 654-663.
    [6] [ T. C. Greenough,D. B. Brettler,F. Kirchhoff, Long-term non-progressive infection with Human Immunodeficiency Virus in a Hemophilia cohort, J Infect Dis, 180 (1999): 1790-1802.
    [7] [ A. B. Gumel,P. N. Shivakumar,B. M. Sahai, A mathematical model for the dynamics of HIV-1 during the typical course of infection, Nonlinear Analysis, 47 (2001): 1773-1783.
    [8] [ M. Hadjiandreou,R. Conejeros,V. S. Vassiliadis, Towards a long-term model construction for the dynamic simulation of HIV infection, Mathematical Biosciences and Engineering, 4 (2007): 489-504.
    [9] [ E. Hernandez-Vargas,R. Middleton, Modeling the three stages in HIV infection, J TheorBiol., 320 (2013): 33-40.
    [10] [ T. Igarashi,C. R. Brown,Y. Endo, Macrophages are the principal reservoir and sustain high virus loads in Rhesus Macaques following the depletion of CD4+ T-cells by a highly pathogenic SIV: Implications for HIV-1 infections of man, Proc Natl Acad Sci., 98 (2001): 658-663.
    [11] [ E. Jones and P. Roemer (sponsors: S. Pankavich and M. Raghupathi), Analysis and simulation of the three-component model of HIV dynamics, SIAM Undergraduate Research Online, 7 (2014), 89–106
    [12] [ D. Kirschner, Using mathematics to understand HIV immunodynamics, Am. Math. Soc., 43 (1996): 191-202.
    [13] [ D. E. Kirschner and A. S. Perelson, A model for the immune response to HIV: AZT treatment studies, Mathematical Population Dynamics: Analysis of Heterogeneity, Volume One: Theory of Epidemics Eds. O. Arino, D. Axelrod, M. Kimmel, and M. Langlais, Wuerz Publishing Ltd., Winnipeg, Canada, (1993), 295–310.
    [14] [ D. Kirschner,G. F. Webb, Immunotherapy of HIV-1 infection, J Biological Systems, 6 (1998): 71-83.
    [15] [ D. Kirschner,G. F. Webb,M. Cloyd, A model of HIV-1 disease progression based on virus-induced lymph node homing-induced apoptosis of CD4+ lymphocytes, J Acquir Immune Dec Syndr, 24 (2000): 352-362.
    [16] [ J. M. Murray,G. Kaufmann,A. D. Kelleher, A model of primary HIV-1 infection, Math Biosci, 154 (1998): 57-85.
    [17] [ M. Nowak and R. May, Virus Dynamics: Mathematical Principles of Immunology and Virology Oxford University Press, NewYork, 2000.
    [18] [ S. Pankavich, The effects of latent infection on the dynamics of HIV, Differential Equations and Dynamical Systems, 24 (2016): 281-303.
    [19] [ S. Pankavich,D. Shutt, An in-host model of HIV incorporating latent infection and viral mutation, Dynamical Systems, Differential Equations, and Applications, AIMS Proceedings, null (2015): 913-922.
    [20] [ S. Pankavich, N. Neri and D. Shutt, Bistable dynamics and Hopf bifurcation in a refined model of the acute stage of HIV infection, submitted, (2015).
    [21] [ S. Pankavich,C. Parkinson, Mathematical analysis of an in-host model of viral dynamics with spatial heterogeneity, Discrete and Continuous Dynamical Systems B, 21 (2016): 1237-1257.
    [22] [ E. Pennisi and J. Cohen, Eradicating HIV from a patient: Not just a dream?, Science, 272 (1996), 1884.
    [23] [ A. S. Perelson, Modeling the Interaction of the Immune System with HIV, Lecture Notes in Biomath. Berlin: Springer, 1989.
    [24] [ A. Perelson,P. Nelson, Mathematical analysis of HIV-1 dynamics in vivo, SIAM Rev., 41 (1999): 3-44.
    [25] [ T. M. Russi, Uncertainty Quantification with Experimental data and Complex System Models, Ph. D. thesis, UC Berkeley, 2010.
    [26] [ W. Y. Tan,H. Wu, Stochastic modeing of the dynamics of CD4+ T-cell infection by HIV and some monte carlo studies, Math Biosci, 147 (1997): 173-205.
    [27] [ E. Vergu,A. Mallet,J. Golmard, A modeling approach to the impact of HIV mutations on the immune system, Comput Biol Med., 35 (2005): 1-24.
  • This article has been cited by:

    1. Hyunjung Lee, Elaine T. Spiller, Susan E. Minkoff, 2019, Dimension reduction and global sensitivity metrics using active subspaces for coupled flow and deformation modeling, 3240, 10.1190/segam2019-3215234.1
    2. C. A. Pena Fernandez, 2020, Fractional Complex-order Hopfield Neural Networks to Analyze the Effect of Drug-resistance in the HIV Infection, 978-1-7281-4248-7, 8, 10.1109/ICACI49185.2020.9177513
    3. Jeffrey M. Hokanson, Paul G. Constantine, Data-Driven Polynomial Ridge Approximation Using Variable Projection, 2018, 40, 1064-8275, A1566, 10.1137/17M1117690
    4. Tongqian Zhang, Junling Wang, Yi Song, Zhichao Jiang, Dynamical Analysis of a Delayed HIV Virus Dynamic Model with Cell-to-Cell Transmission and Apoptosis of Bystander Cells, 2020, 2020, 1076-2787, 1, 10.1155/2020/2313102
    5. James C. Gross, Pranay Seshadri, Geoff Parks, 2020, Optimisation with Intrinsic Dimension Reduction: A Ridge Informed Trust-Region Method, 978-1-62410-595-1, 10.2514/6.2020-0157
    6. Charlotte Lew, 2022, Neural Network Modeling of HIV Acute and Chronic Phases With and Without Antiretroviral Intervention, 978-1-6654-7184-8, 47, 10.1109/TransAI54797.2022.00014
    7. Wouter Edeling, On the Deep Active-Subspace Method, 2023, 11, 2166-2525, 62, 10.1137/21M1463240
    8. Zixiao Ma, Bai Cui, Zhaoyu Wang, Dongbo Zhao, Parameter Reduction of Composite Load Model Using Active Subspace Method, 2021, 36, 0885-8950, 5441, 10.1109/TPWRS.2021.3078671
    9. Qiang-hui Xu, Ji-cai Huang, Yue-ping Dong, Yasuhiro Takeuchi, A Delayed HIV Infection Model with the Homeostatic Proliferation of CD4+ T Cells, 2022, 38, 0168-9673, 441, 10.1007/s10255-022-1088-2
    10. Xifu Sun, Barry Croke, Anthony Jakeman, Stephen Roberts, Benchmarking Active Subspace methods of global sensitivity analysis against variance-based Sobol' and Morris methods with established test functions, 2022, 149, 13648152, 105310, 10.1016/j.envsoft.2022.105310
    11. Cameron Clarke, Stephen Pankavich, Three-stage modeling of HIV infection and implications for antiretroviral therapy, 2024, 88, 0303-6812, 10.1007/s00285-024-02056-1
    12. Soraya Terrab, Stephen Pankavich, Global sensitivity analysis of plasma instabilities via active subspaces, 2024, 134, 10075704, 107994, 10.1016/j.cnsns.2024.107994
    13. Wouter Edeling, Maxime Vassaux, Yiming Yang, Shunzhou Wan, Serge Guillas, Peter V. Coveney, Global ranking of the sensitivity of interaction potential contributions within classical molecular dynamics force fields, 2024, 10, 2057-3960, 10.1038/s41524-024-01272-z
    14. Valentin Breaz, Richard Wilkinson, 2024, Randomized Maximum Likelihood Via High-Dimensional Bayesian Optimization, 979-8-3503-4485-1, 5300, 10.1109/ICASSP48485.2024.10446755
    15. Sanket Jantre, Nathan M. Urban, Xiaoning Qian, Byung-Jun Yoon, 2024, Learning Active Subspaces for Effective and Scalable Uncertainty Quantification in Deep Neural Networks, 979-8-3503-4485-1, 5330, 10.1109/ICASSP48485.2024.10448265
    16. Hui Duan, Giray Okten, DERIVATIVE-BASED SHAPLEY VALUE FOR GLOBAL SENSITIVITY ANALYSIS AND MACHINE LEARNING EXPLAINABILITY , 2025, 15, 2152-5080, 1, 10.1615/Int.J.UncertaintyQuantification.2024051548
    17. Yueping Dong, Jicai Huang, Yasuhiro Takeuchi, Qianghui Xu, 2024, Chapter 2, 978-981-97-7849-2, 11, 10.1007/978-981-97-7850-8_2
    18. Xia Wang, Yue Wang, Yueping Dong, Libin Rong, Dynamics of a delayed HIV infection model with cell-to-cell transmission and homeostatic proliferation, 2024, 139, 2190-5444, 10.1140/epjp/s13360-024-05845-1
    19. Jing Cai, Jun Zhang, Kai Wang, Zhixiang Dai, Zhiliang Hu, Yueping Dong, Zhihang Peng, Evaluating the long-term effects of combination antiretroviral therapy of HIV infection: a modeling study, 2025, 90, 0303-6812, 10.1007/s00285-025-02196-y
  • Reader Comments
  • © 2017 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(4237) PDF downloads(568) Cited by(19)

Article outline

Figures and Tables

Figures(17)  /  Tables(1)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog