
MBW:Fundamental Mathematics of the Delay Differential Equations in the HIV1 Virus ModelFrom MathBioThis is an article is written by Sai'd Azzam and John Nelson Kane ContentsExecutive SummaryIn this mathematical model we attempt to better understand the dynamics of the HIV1 infection by incorporating an intracellular delay, which is the time period between an antiviral drug reaching a cell and the effectiveness of the drug. In the past, models have made the unrealistic assumption that the intracellular delay is constant or there is no delay after ingestion of the drug. This model introduces this delay in the form of a continuous distribution with the underlying principle that different cells will be associated with a different delay time. One fundamental equation of this model is known as an integrodifferential equation. The focus of this study will be on the analytic methods used to transform this integrodifferential equation into a set of linear 1st order differential equations in which analytical solutions can be obtained. This will be done with a method known as the linear chain trick. By incorporating this method it will be shown that a time delay equation can be transformed into an ordinary equation of only current time, thus making analytic solutions possible. As a final conclusion we will present the final set of 1st order linear differential equations that will equally describe this model. IntroductionThe mathematics describing this model are motivated by the research and attempted understanding of the Human Immunodeficiency Virus (HIV). HIV is a replicating virus and the cause of acquired immunodeficiency syndrome (AIDS). HIV1 is the most common and pathogenic strain of the virus and is typically associated with three stages. Once a patient is first infected, the amount of virus found in the blood stream is relatively high. This high concentration may last several weeks to a few months at which point the concentration decreases to some stable lower level. Following the beginning of this stable level of virus concentration, a typically long period follows where the CD4+ Tcell count slowly decreases. A CD4+ Tcell is also known as a T helper cell that is a type of white blood cell used to help the activity of other immune cells. This usually long period of a decrease in Tcells appears to be a quasisteady state, because although it appears to be stable, the virus concentration is noticed to change slightly when analyzed on a long time scale. This quasisteady state can last for up to a few years with little change. In the final stage of the disease, the CD4+ Tcell count drops to a point in which the immune system fails and the viral load increases rapidly, usually causing fatality. Much research has been done to better understand this middle period in which the Tcells slowly decrease and the virus concentration stays at a relatively stable value. If one can provide the correct combination of antiretroviral drugs with treatment to better stabilize the white blood cells than one could dramatically prolong life and even essentially make the disease nonfatal to the human body. This model attempts to analytically describe the middle stage, the quasisteady state of the virus. Previous models have attempted to describe this stage with either a constant delay, or the special case of no delay at all. These are not realistic for obvious reasons, although constant delays allow for analytic solutions to be easily obtain for the differential equations in those models are first order and linear. This model will incorporate a gamma distribution for the time delay with the physical interpretation that there is no constant time for which the decay occurs in any given cell. The density function of a gamma distribution is shown below:
Delay / BackgroundThere are two different types of delays that will effect the infectious/treatment process: a pharmacological delay and an intracellular delay. A pharmacological delay is the time between the ingestion of the drug and it’s first appearance within the cells. An intracellular delay is the time between the initial infection of a cell and the release of new viruses. Studies show from Perelson et al. that the pharmacological delay is around two hours while the intracellular delay is of the magnitude of a few days. It is also must be considered that the pharmacological delay only occurs at the commencement of treatment, while intracellular delay occurs continuously. For this reason we can ignore pharmacological delay and we will do so for the following model. The intracellular delay will be in the form of a gamma distribution. As shown in the following methods, this will allow for an analytical solution to be found. The gamma distribution is also considered to be a reasonable delay distribution within a cell from the biological sense. Mathematical Model With Variable/Parameter DefinitionsThe main goal of this paper is to have a biologically plausible answer, and to analytically transform an integrodifferential equation into a system of first order differential equations. The parameters associated with the model are: the density of productively infected cells, I, the concentration of infections virus, V(I), the viral clearance rate constant per day, c, and the concentration of virus has been rendered noninfections by the protease inhibitors, V(NI). The total concentration of virus is V=V(I) + V(NI). We can also assume that the density of the noninfected cells, T, will remain constant. This is a fair assumption since the data are only examined for a short period of time following the initial ingestion of the drug. In the absence of delay, Perelson et al. proposed the following model:
We make note that k is the infection rate constant, p is the rate at which a productively infected cell release viruses. Productivity infected cells die at a rate per cell. Plasma viruses are cleared at rate c per virus. respresents the drug efficiency. If than the drug has no effective and if than the drug is 100% effective. As noted above this model will only account for the intracellular delay, and the pharmacological delay will be ignored. With the addition of the delay the following model is proposed in the form of integrodifferential equations as follows:
We multiplied the gamma distribution, f(t'), by the change of the infectious virus equation, , and integrate to find the probability.
The equations (4b) and (4c) are already in the form of first order differential equations, thus we will not transform them.
Throughout the analysis we will assume that the system being treated is at steady state prior to being treated. We will also make note of at the time the drugs initially takes effect. This allows for the ignorance of the pharmacological delay. Noting that there are zero noninfectious virus cells at we can allow for initial conditions as follows:
The intracellular delay is not always known to an exact, nor is it the same value for every cell undergoing the therapy process. In fact it is unrealistic to represent this delay with a constant thus is will be represented with a gamma distribution known as the delay distribution .
with a mean, , a variance, and a max peak, . This distribution is widely used in mathematical biology because it is a plausible delay and allows for analytic solutions to be obtained. Converting time delay defferential equations into ordinary differential equationFrom 4(a) we wish to eliminate the term , as a first step to approach an analytical solution, thus we make a substitution. We will define two new variables and and
We need to solve for b in terms of and thus,
Here we will solve for k , in term of and
Our integrodifferential equation is of the form:
Substituting for b and k with the expressions above in terms of and , and we get:
As we can see, we cancel the and terms. Thus we arrive at:
Thus our equation is now in terms of and . The equation is now similar to (4a) with the exception of in replacement of b, in replacement of k and no term. We now have a gamma distribution in terms of and is defined as:
We now introduce the function . This term is defined as:
With the gamma distribution expanded we have:
To begin the linear chain trick, we introduce a new variable u = tt' and substitute it into
solving for t', we get:
After substituting u into we get:
Then we differentiate the integrand with respect to the parameter t. We notice that the integral is over du thus we can use Leibniz's rule of differentiating integrals. This derivative with respect to t, applying the chain rule is:
Then we simplify the terms and we get:
Since , we get
If we look closely, we notice that the above terms mean,
This is only valid for thus we have to get a solution in term of . Thus we plug 1 for the equation, and we get:
We differentiate in the same manner with respect to the parameter t. First we manipulate for simplification. We pull the term out of the intergal, which is considered a constant in terms of the variable u.
Using Leibniz's rule of differentiation of an integrand we differentiate with respect to t.
Again, if we look closely is just:
After all of the above manipulations, now the integrodifferential equations given in (4a), (4b) and (4c) are now equivalent to the following system of ordinary 1st order linear differential equations.
From (4b), our initial linear differential equation is:
We will assume 100% effectiveness ( ), thus we arrive at:
From (4c), our initial linear differential equation is:
Assuming 100% effectiveness ( ), we get:
for are auxiliary functions and are related to the delay distribution as follows:
Results and DiscussionWe began with a system of three delay differential equations, one being an integrodifferential equation. Throughout this paper we have shown that through a series of steps and a method known as the linear chain trick, the integrodifferential equation can be rewritten in terms of a system of 1st order differential equations that can be solved for analytically. Our delay is in the form of a gamma distribution. One must account for all possibilities of the delay, thus it is needed to integrate over the whole distribution density function. This is the reason for an integrand appearing in the model to give us an integrodifferential equation. The initial model is shown below with the integrodifferential equation being the first presented.
It has been shown through that with the linear chain trick, the integrodifferential equation can be transformed into a system of 1st order differential equations. Noticing that the original 2nd and 3rd equations in model are already in linear form they will remain in the final model. The only difference is that the drug has been assumed to be 100% effective. It has been shown that the above model presented is equivalent to the system of 1st order linear differential equations shown below:
Now that a first order system has been obtained, one can go further by working towards an analytic result. Noticing that there are now n + 3 linear equations one might find it more useful to solve using a software package with the n that is desired. To further investigate how the mathematics presented here will model a biological system such as HIV1 dynamics and the associated delay, one could study the effects of given parameters with a sensitivity analysis.
External Links1. J. E. Mittler, B. Sulzer, A. U. Neumann, and A. S. Perelson. Influence of delayed viral production on viral dynamics in HIV1 infected patients. Math. Biosciences, 152:143163, 1998. [1] 2. D.D. Ho, A.U. Neumann, A.S. Perelson, W. Chen, J.M. Leonard, M. Markowitz Rapid turnover of plasma virions and CD4 lymphocytes in HIV1 infection Nature, 373 (1995), p. 123. [2] 3. Differentiation under the integral sign [3] 4. S. G. Azzam, J. N. Kane. Delayed Production and Viral Dynamics of the HIV1 Virus [4] 