
MBW:Estimation of HIV/AIDS parametersFrom MathBioContentsExecutive SummaryThis paper determines HIV/AIDS disease parameters for a basic HIV model from measurements of healthy CD4+T concentrations in plasma and the viral load in plasma. First, the minimum number of necessary measurements is obtained. Next, the stages of HIV progression when parameter estimation is impossible are determined. After that, using system control techniques, an estimation algorithm of HIV parameters is created. The use of adaptive observers helps filter out the noise in measurement data and predict all six parameters reasonably well. Project SummaryThis project covers aspects of mathematical biology and biological systems, specifically within Immunology. This is an HIV/AIDS specific mathematical model. The model uses concentration of certain T cells as well as the concentration of the virus. The healthy CD4+ T cells are created at a constant rate, die depending on their concentration, and are infected at a rate proportional to both their concentration and virus concentration. The model used to effectively capture that information is:
The infected T cells reproduce and die at a concentrationproportional rate. The virus is produced by the infected T cells and also dies proportionally to the concentration. Background InformationHIV, human immunodeficiency virus, is a virus which causes the immune system to fail, making the host vulnerable to opportunistic infections. It infects approximately 0.6% of the world's population. Since its discovery in 1981, it has taken approximately 25 million lives. Infection StagesThe infection happens in 4 stages: incubation, acute infection, latency, and AIDS. The initial incubation period lasts twofour weeks and generally has no symptoms. During the acute infection stage, however, the patients show symptoms such as fever, lymphadenopathy, rashes, pharyngitis, and sores. The symptoms are difficult to link specifically to HIV because they are very scattered. The acute stage lasts approximately 4 weeks. During this phase, the virus replicates itself quickly. The latency stage occurs next; it is a stage when there are few to no symptoms showing. It can last anywhere between a couple weeks and 20 years. During this stage, the virus count stays relatively low. This can give the impression that the disease is under control. However, following the latency stage comes AIDS, when the T cell count drops below 200 cells/microliter. At that point, the immune system can't fight off opportunistic infections. Virus Reproduction
The virus enters into the host cell by fusing with the cell membrane and releasing the contents of the virus capsule into the interior of the cell. Following entering the cell, the virus takes over the protein producing machinery of the cell, and uses the cell's capabilities to replicate. See the image for some of the details of replication.
The virus replicates itself very fast, to virions per day. The mutation rate is approximately per nucleotide base, due to the errors in the reverse transcriptase process (when the RNA of the virus is used to create cDNA). The errorprone process is responsible for the high number of virus mutations created in a single host during one day. Since a single T cell may also be infected with several mutants of HIV, it becomes difficult to determine a proper course of treatment very quickly. The drug resistance caused by the frequent mutations is one of the reasons why HIV is such a difficult disease to approach. For more information on Virus Replication please see: [1] Mathematical Models BackgroundMathematical models of HIV dynamics describe the interactions between the immune system and the virus at different stage. These models give insights to the dynamics of the disease. For more background on the dynamics of this disease please see MBW: HIV Dynamics
Mathematical ModelThe basic model of HIV/AIDS is as follows:
In the model, is the concentration of healthy CD4+ T cells, is the concentration of infected CD4+ T cells, and is the concentration of virus in plasma. s is the rate at which new CD4+ T cells are created, d is the death rate of T cells, is the infection rate, is the death rate of infected T cells, is the death rate of virions, and k is the rate of creation of virions. It has been previously shown that the system is observable and identifiable. Since only the concentration of uninfected T cells and concentration of virus can be measured, the outputs are redefined as . Determining Conditions for Parameter IdentificationThe system can be redefined as:
Where values define a onetoone map. The conditions are that and . From clinical data, and . Therefore,
With the measurements labeled so that the superscript indicates the number of the measurement, and the parameter d indicating the time difference, the following system can be set up to calculate the first three theta values based on the measured concentrations.
The parameters can be similarly derived. The method will fail to predict the parameters when the viral load remains constant, so during the latency stage and during short periods after treatment, the data will not provide for useful results. The main difficulty with this method though is that it fails with noisy data. Clinical data is usually noisy, with variance of 20 for the T cell concentration and log variance of 0.2 for virion concentration. The results section shows how this system performs with parameter prediction. Estimation of Parameters Using Adaptive ObserversTo see a general process on how to estimate parameters refer to: [2] In order to analyze noisy data, the system can be transformed into observer form by setting z to be the following:
The observer form is as shown:
The filtered transformation can be defined as well.
In those expressions, b has to be greater than zero, so that the system can be transformed into the adaptive observer form and then an adaptive observer can be designed.
The are symmetric positive definite matrices, </math>k_1</math> is a negative constant, . SimulationThe simulation was set up in Matlab Simbiology package, in order to obtain data for parameter estimation. The parameters were as follows: . The time is on the x axis and is in days. The HIV progression can be seen for all three variables both from the paper and the Simbiology run. They are very similar. ResultsThe parameter estimation done on s, d and shows that without the adaptive observer technique, the results tend to be very noisy. When the adaptive observer technique is used, the parameters s, d and look much smoother. The estimations are still not very smooth because the virus and infected cells have a much higher turnover (a factor of 10 compared to d for the infected cell death and a factor of over 100 for the virus death). DiscussionThe parameters can be determined using the adaptive observer method. They are most easily determined during the early infection stage when there is sufficient change in viral load. Therefore the main question to be asked is how likely is it that the disease can be monitored at the early stage? (See MBW:When to Initiate HIV Therapy for more information on HIV at an early stage) In any case, by disturbing the steady state of the disease with drugs, the parameters can be calculated as well. By monitoring the parameters for the disease in one patient, the drug resistance development could be potentially monitored. Potentially, a study of how the parameters and k change with resistance could be done. ConclusionThe determination of parameters for the HIV/AIDS model was studied. It was determined that at least 4 measurements of CD4+ T cell load and 5 measurements of viral load need to be obtained in order to estimate model parameters. The model does not allow for estimation of HIV/AIDS parameters during the latency stage. Using adaptive observers, the parameters can be calculated while the data is obtained, and all 6 parameters needed converge at close to the correct value when the model is simulated. For a more detailed description of what each these parameters mean see MBW:HIV_Dynamics. Recent PapersOver the years there have been about 41 different papers citing the original paper, "Estimation of HIV/AIDS Parameters". One specific paper touches the topic of the estimation of parameters and extends upon the research that was already done. In this paper the concept of algebraic identifiability of a three dimensional HIV/AIDS dynamic model. Identifiability is defined as a condition that a model must satisfy in order to be statistically inferred. The researchers proposed a new method of finding the parameters of the model by using the multiple time point method. Previously, the highorder derivative method was used in order to find the missing parameters. The paper uses data and the new method to compute the new parameters and presents the new method. ReferencesXiaohua Xia. "Estimation of HIV/AIDS Parameters." Automatica 39 (2003) 19831988. Xia, X. Moog, C.H. (2003) "Identifiability of nonlinear systems with application to HIV/AIDS models." IEEE Transactions on Automatic Control, 48, 330336. Wu, H. Zhu, H. Miao, H. Perelson, A. (2008) "Parameter Identifiability and Estimation of HIV/AIDS Dynamic Models." Bulletin of Mathematical Biology, 70, 785799. [3] 