skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Stochastic investigation of HIV infection and the emergence of drug resistance
Drug-resistant HIV-1 has caused a growing concern in clinic and public health. Although combination antiretroviral therapy can contribute massively to the suppression of viral loads in patients with HIV-1, it cannot lead to viral eradication. Continuing viral replication during sub-optimal therapy (due to poor adherence or other reasons) may lead to the accumulation of drug resistance mutations, resulting in an increased risk of disease progression. Many studies also suggest that events occurring during the early stage of HIV-1 infection (i.e., the first few hours to days following HIV exposure) may determine whether the infection can be successfully established. However, the numbers of infected cells and viruses during the early stage are extremely low and stochasticity may play a critical role in dictating the fate of infection. In this paper, we use stochastic models to investigate viral infection and the emergence of drug resistance of HIV-1. The stochastic model is formulated by a continuous-time Markov chain (CTMC), which is derived based on an ordinary differential equation model proposed by Kitayimbwa et al. that includes both forward and backward mutations. An analytic estimate of the probability of the clearance of HIV infection of the CTMC model near the infection-free equilibrium is obtained by a multitype branching process approximation. The analytical predictions are validated by numerical simulations. Unlike the deterministic dynamics where the basic reproduction number $$ \mathcal{R}_0 $$ serves as a sharp threshold parameter (i.e., the disease dies out if $$ \mathcal{R}_0 < 1 $$ and persists if $$ \mathcal{R}_0 > 1 $$), the stochastic models indicate that there is always a positive probability for HIV infection to be eradicated in patients. In the presence of antiretroviral therapy, our results show that the chance of clearance of the infection tends to increase although drug resistance is likely to emerge.  more » « less
Award ID(s):
1950254
PAR ID:
10350818
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Mathematical Biosciences and Engineering
Volume:
19
Issue:
2
ISSN:
1551-0018
Page Range / eLocation ID:
1174 to 1194
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; (, Mathematical Biosciences and Engineering)
    In this paper we develop a four compartment within-host model of nutrition and HIV. We show that the model has two equilibria: an infection-free equilibrium and infection equilibrium. The infection free equilibrium is locally asymptotically stable when the basic reproduction number $$ \mathcal{R}_0 < 1 $$, and unstable when $$ \mathcal{R}_0 > 1 $$. The infection equilibrium is locally asymptotically stable if $$ \mathcal{R}_0 > 1 $$ and an additional condition holds. We show that the within-host model of HIV and nutrition is structured to reveal its parameters from the observations of viral load, CD4 cell count and total protein data. We then estimate the model parameters for these 3 data sets. We have also studied the practical identifiability of the model parameters by performing Monte Carlo simulations, and found that the rate of clearance of the virus by immunoglobulins is practically unidentifiable, and that the rest of the model parameters are only weakly identifiable given the experimental data. Furthermore, we have studied how the data frequency impacts the practical identifiability of model parameters. 
    more » « less
  2. ; ; ; (, International Journal of Bifurcation and Chaos)
    Individuals infected by human immunodeficiency virus (HIV) are under oxidative stress due to the imbalance between reactive oxygen species (ROS) production and elimination. This paper presents a mathematical model with the cytotoxic T lymphocytes (CTL) immune response to examine the role of ROS in the dynamics of HIV infection. We classify the equilibria of the model and study the stability of these equilibria. Numerical simulations show that incorporating ROS and CTL immune response into the model leads to very rich dynamics, including bistable phenomena and periodic solutions. Although the current antiretroviral therapy can suppress viral load to the undetectable level, it cannot eradicate the virus. A high level of ROS may be a factor for HIV persistence in patients despite suppressive therapy. These results suggest that oxidative damage and anti-oxidant therapy should be considered in the study of HIV infection and treatment. 
    more » « less
  3. ; ; ; ; ; ; (, Viruses)
    null (Ed.)
    Multidrug resistance continues to be a barrier to the effectiveness of highly active antiretroviral therapy in the treatment of human immunodeficiency virus 1 (HIV-1) infection. Darunavir (DRV) is a highly potent protease inhibitor (PI) that is oftentimes effective when drug resistance has emerged against first-generation inhibitors. Resistance to darunavir does evolve and requires 10–20 amino acid substitutions. The conformational landscapes of six highly characterized HIV-1 protease (PR) constructs that harbor up to 19 DRV-associated mutations were characterized by distance measurements with pulsed electron double resonance (PELDOR) paramagnetic resonance spectroscopy, namely double electron–electron resonance (DEER). The results show that the accumulated substitutions alter the conformational landscape compared to PI-naïve protease where the semi-open conformation is destabilized as the dominant population with open-like states becoming prevalent in many cases. A linear correlation is found between values of the DRV inhibition parameter Ki and the open-like to closed-state population ratio determined from DEER. The nearly 50% decrease in occupancy of the semi-open conformation is associated with reduced enzymatic activity, characterized previously in the literature. 
    more » « less
  4. ; ; (, Mathematical Biosciences and Engineering)
    In this paper, we introduce a novel multi-scale network model of two epidemics: HIV infection and opioid addiction. The HIV infection dynamics is modeled on a complex network. We determine the basic reproduction number of HIV infection, $$ \mathcal{R}_{v} $$, and the basic reproduction number of opioid addiction, $$ \mathcal{R}_{u} $$. We show that the model has a unique disease-free equilibrium which is locally asymptotically stable when both $$ \mathcal{R}_{u} $$ and $$ \mathcal{R}_{v} $$ are less than one. If $$ \mathcal{R}_{u} > 1 $$ or $$ \mathcal{R}_{v} > 1 $$, then the disease-free equilibrium is unstable and there exists a unique semi-trivial equilibrium corresponding to each disease. The unique opioid only equilibrium exist when the basic reproduction number of opioid addiction is greater than one and it is locally asymptotically stable when the invasion number of HIV infection, $$ \mathcal{R}^{1}_{v_i} $$ is less than one. Similarly, the unique HIV only equilibrium exist when the basic reproduction number of HIV is greater than one and it is locally asymptotically stable when the invasion number of opioid addiction, $$ \mathcal{R}^{2}_{u_i} $$ is less than one. Existence and stability of co-existence equilibria remains an open problem. We performed numerical simulations to better understand the impact of three epidemiologically important parameters that are at the intersection of two epidemics: $$ q_v $$ the likelihood of an opioid user being infected with HIV, $$ q_u $$ the likelihood of an HIV-infected individual becoming addicted to opioids, and $$ \delta $$ recovery from opioid addiction. Simulations suggest that as the recovery from opioid use increases, the prevalence of co-affected individuals, those who are addicted to opioids and are infected with HIV, increase significantly. We demonstrate that the dependence of the co-affected population on $$ q_u $$ and $$ q_v $$ are not monotone. 
    more » « less
  5. ; (, PLOS Computational Biology)
    Regoes, Roland R. (Ed.)
    While highly active antiretroviral therapy (HAART) is successful in controlling the replication of Human Immunodeficiency Virus (HIV-1) in many patients, currently there is no cure for HIV-1, presumably due to the presence of reservoirs of the virus. One of the least studied viral reservoirs is the brain, which the virus enters by crossing the blood-brain barrier (BBB) via macrophages, which are considered as conduits between the blood and the brain. The presence of HIV-1 in the brain often leads to HIV associated neurocognitive disorders (HAND), such as encephalitis and early-onset dementia. In this study we develop a novel mathematical model that describes HIV-1 infection in the brain and in the plasma coupled via the BBB. The model predictions are consistent with data from macaques infected with a mixture of simian immunodeficiency virus (SIV) and simian-human immunodeficiency virus (SHIV). Using our model, we estimate the rate of virus transport across the BBB as well as viral replication inside the brain, and we compute the basic reproduction number. We also carry out thorough sensitivity analysis to define the robustness of the model predictions on virus dynamics inside the brain. Our model provides useful insight into virus replication within the brain and suggests that the brain can be an important reservoir causing long-term viral persistence. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email