Chronic hepatitis B virus (HBV) infection is strongly associated with increased risk of liver cancer and cirrhosis. While existing treatments effectively inhibit the HBV life cycle, viral rebound frequently occurs following treatment interruption. Consequently, functional cure rates of chronic HBV infection remain low and there is increased interest in a novel treatment modality, capsid assembly modulators (CAMs). Here, we develop a multiscale mathematical model of CAM treatment in chronic HBV infection. By fitting the model to participant data from a phase I trial of the first-generation CAM vebicorvir, we estimate the drug's dose-dependent effectiveness and identify the physiological mechanisms that drive the observed biphasic decline in HBV DNA and RNA, and mechanistic differences between HBeAg-positive and negative infection. Finally, we demonstrate analytically and numerically that the relative change of HBV RNA more accurately reflects the antiviral effectiveness of a CAM than the relative change in HBV DNA.

Quantification of intrahepatic covalently closed circular DNA (cccDNA) is a key for evaluating an elimination of hepatitis B virus (HBV) in infected patients. However, quantifying cccDNA requires invasive methods such as a liver biopsy, which makes it impractical to access the dynamics of cccDNA in patients. Although HBV RNA and HBV core-related antigens (HBcrAg) have been proposed as surrogate markers for evaluating cccDNA activity, they do not necessarily estimate the amount of cccDNA. Here, we employed a recently developed multiscale mathematical model describing intra- and intercellular viral propagation and applied it in HBV-infected patients under treatment. We developed a model that can predict intracellular HBV dynamics by use of extracellular viral markers, including HBsAg, HBV DNA, and HBcrAg in peripheral blood. Importantly, the model prediction of the amount of cccDNA in patients over time was confirmed to be well correlated with the data for quantified cccDNA by paired liver biopsy. Thus, our method combining classic and emerging surrogate markers enables us to predict the decay dynamics of cccDNA in patients undergoing treatment.

PGT121 is a broadly neutralizing antibody in clinical development for the treatment and prevention of HIV-1 infection via passive administration. PGT121 targets the HIV-1 V3-glycan and demonstrated potent antiviral activity in a phase I clinical trial. Resistance to PGT121 monotherapy rapidly occurred in the majority of participants in this trial with the sampled rebound viruses being entirely resistant to PGT121 mediated neutralization. However, two individuals experienced long-term ART-free viral suppression following antibody infusion and retained sensitivity to PGT121 upon viral rebound. Here, we develop mathematical models of the HIV-1 dynamics during this phase I clinical trial. We utilize these models to understand the dynamics leading to PGT121 resistance and to identify the mechanisms driving the observed long-term viral control. Our modeling highlights the importance of the relative fitness difference between PGT121 sensitive and resistant subpopulations prior to treatment. Specifically, by fitting our models to data, we identify the treatment-induced competitive advantage of previously existing or newly generated resistant population as a primary driver of resistance. Finally, our modeling emphasizes the high neutralization ability of PGT121 in both participants who exhibited long-term viral control.

Chronic infection with hepatitis B virus (HBV) is caused by the persistence of closed circular DNA (cccDNA) in the nucleus of infected hepatocytes. Despite available therapeutic anti-HBV agents, eliminating the cccDNA remains challenging. Thus, quantifying and understanding the dynamics of cccDNA are essential for developing effective treatment strategies and new drugs. However, such study requires repeated liver biopsy to measure the intrahepatic cccDNA, which is basically not accepted because liver biopsy is potentially morbid and not common during hepatitis B treatment. We here aimed to develop a noninvasive method for quantifying cccDNA in the liver using surrogate markers in peripheral blood. We constructed a multiscale mathematical model that explicitly incorporates both intracellular and intercellular HBV infection processes. The model, based on age-structured partial differential equations, integrates experimental data from in vitro and in vivo investigations. By applying this model, we roughly predicted the amount and dynamics of intrahepatic cccDNA within a certain range using specific viral markers in serum samples, including HBV DNA, HBsAg, HBeAg, and HBcrAg. Our study represents a significant step towards advancing the understanding of chronic HBV infection. The noninvasive quantification of cccDNA using our proposed method holds promise for improving clinical analyses and treatment strategies. By comprehensively describing the interactions of all components involved in HBV infection, our multiscale mathematical model provides a valuable framework for further research and the development of targeted interventions.

Mathematical models of viral infection have been developed, fitted to data, and provide insight into disease pathogenesis for multiple agents that cause chronic infection, including HIV, hepatitis C, and B virus. However, for agents that cause acute infections or during the acute stage of agents that cause chronic infections, viral load data are often collected after symptoms develop, usually around or after the peak viral load. Consequently, we frequently lack data in the initial phase of viral growth, i.e., when pre-symptomatic transmission events occur. Missing data may make estimating the time of infection, the infectious period, and parameters in viral dynamic models, such as the cell infection rate, difficult. However, having extra information, such as the average time to peak viral load, may improve the robustness of the estimation. Here, we evaluated the robustness of estimates of key model parameters when viral load data prior to the viral load peak is missing, when we know the values of some parameters and/or the time from infection to peak viral load. Although estimates of the time of infection are sensitive to the quality and amount of available data, particularly pre-peak, other parameters important in understanding disease pathogenesis, such as the loss rate of infected cells, are less sensitive. Viral infectivity and the viral production rate are key parameters affecting the robustness of data fits. Fixing their values to literature values can help estimate the remaining model parameters when pre-peak data is missing or limited. We find a lack of data in the pre-peak growth phase underestimates the time to peak viral load by several days, leading to a shorter predicted growth phase. On the other hand, knowing the time of infection (e.g., from epidemiological data) and fixing it results in good estimates of dynamical parameters even in the absence of early data. While we provide ways to approximate model parameters in the absence of early viral load data, our results also suggest that these data, when available, are needed to estimate model parameters more precisely.

Plus-strand RNA viruses are the largest group of viruses. Many are human pathogens that inflict a socio-economic burden. Interestingly, plus-strand RNA viruses share remarkable similarities in their replication. A hallmark of plus-strand RNA viruses is the remodeling of intracellular membranes to establish replication organelles (so-called "replication factories"), which provide a protected environment for the replicase complex, consisting of the viral genome and proteins necessary for viral RNA synthesis. In the current study, we investigate pan-viral similarities and virus-specific differences in the life cycle of this highly relevant group of viruses. We first measured the kinetics of viral RNA, viral protein, and infectious virus particle production of hepatitis C virus (HCV), dengue virus (DENV), and coxsackievirus B3 (CVB3) in the immuno-compromised Huh7 cell line and thus without perturbations by an intrinsic immune response. Based on these measurements, we developed a detailed mathematical model of the replication of HCV, DENV, and CVB3 and showed that only small virus-specific changes in the model were necessary to describe the in vitro dynamics of the different viruses. Our model correctly predicted virus-specific mechanisms such as host cell translation shut off and different kinetics of replication organelles. Further, our model suggests that the ability to suppress or shut down host cell mRNA translation may be a key factor for in vitro replication efficiency, which may determine acute self-limited or chronic infection. We further analyzed potential broad-spectrum antiviral treatment options in silico and found that targeting viral RNA translation, such as polyprotein cleavage and viral RNA synthesis, may be the most promising drug targets for all plus-strand RNA viruses. Moreover, we found that targeting only the formation of replicase complexes did not stop the in vitro viral replication early in infection, while inhibiting intracellular trafficking processes may even lead to amplified viral growth.

In HIV-1-infected individuals, transmitted/founder (TF) virus contributes to establish new infection and expands during the acute phase of infection, while chronic control (CC) virus emerges during the chronic phase of infection. TF viruses are more resistant to interferon-alpha (IFN-α)-mediated antiviral effects than CC virus, however, its virological relevance in infected individuals remains unclear. Here we perform an experimental-mathematical investigation and reveal that IFN-α strongly inhibits cell-to-cell infection by CC virus but only weakly affects that by TF virus. Surprisingly, IFN-α enhances cell-free infection of HIV-1, particularly that of CC virus, in a virus-cell density-dependent manner. We further demonstrate that LY6E, an IFN-stimulated gene, can contribute to the density-dependent enhancement of cell-free HIV-1 infection. Altogether, our findings suggest that the major difference between TF and CC viruses can be explained by their resistance to IFN-α-mediated inhibition of cell-to-cell infection and their sensitivity to IFN-α-mediated enhancement of cell-free infection.

Experimental Zika virus infection in non-human primates results in acute viral load dynamics that can be well-described by mathematical models. The inoculum dose that would be received in a natural infection setting is likely lower than the experimental infections and how this difference affects the viral dynamics and immune response is unclear. Here we study a dataset of experimental infection of non-human primates with a range of doses of Zika virus. We develop new models of infection incorporating both an innate immune response and viral interference with that response. We find that such a model explains the data better than models with no interaction between virus and the immune response. We also find that larger inoculum doses lead to faster dynamics of infection, but approximately the same total amount of viral production.

Repurposed drugs that are safe and immediately available constitute a first line of defense against new viral infections. Despite limited antiviral activity against SARS-CoV-2, several drugs are being tested as medication or as prophylaxis to prevent infection. Using a stochastic model of early phase infection, we evaluate the success of prophylactic treatment with different drug types to prevent viral infection. We find that there exists a critical efficacy that a treatment must reach in order to block viral establishment. Treatment by a combination of drugs reduces the critical efficacy, most effectively by the combination of a drug blocking viral entry into cells and a drug increasing viral clearance. Below the critical efficacy, the risk of infection can nonetheless be reduced. Drugs blocking viral entry into cells or enhancing viral clearance reduce the risk of infection more than drugs that reduce viral production in infected cells. The larger the initial inoculum of infectious virus, the less likely is prevention of an infection. In our model, we find that as long as the viral inoculum is smaller than 10 infectious virus particles, viral infection can be prevented almost certainly with drugs of 90% efficacy (or more). Even when a viral infection cannot be prevented, antivirals delay the time to detectable viral loads. The largest delay of viral infection is achieved by drugs reducing viral production in infected cells. A delay of virus infection flattens the within-host viral dynamic curve, possibly reducing transmission and symptom severity. Thus, antiviral prophylaxis, even with reduced efficacy, could be efficiently used to prevent or alleviate infection in people at high risk.

Treating macaques with an anti-α4β7 antibody under the umbrella of combination antiretroviral therapy (cART) during early SIV infection can lead to viral remission, with viral loads maintained at < 50 SIV RNA copies/ml after removal of all treatment in a subset of animals. Depletion of CD8+ lymphocytes in controllers resulted in transient recrudescence of viremia, suggesting that the combination of cART and anti-α4β7 antibody treatment led to a state where ongoing immune responses kept the virus undetectable in the absence of treatment. A previous mathematical model of HIV infection and cART incorporates immune effector cell responses and exhibits the property of two different viral load set-points. While the lower set-point could correspond to the attainment of long-term viral remission, attaining the higher set-point may be the result of viral rebound. Here we expand that model to include possible mechanisms of action of an anti-α4β7 antibody operating in these treated animals. We show that the model can fit the longitudinal viral load data from both IgG control and anti-α4β7 antibody treated macaques, suggesting explanations for the viral control associated with cART and an anti-α4β7 antibody treatment. This effective perturbation to the virus-host interaction can also explain observations in other nonhuman primate experiments in which cART and immunotherapy have led to post-treatment control or resetting of the viral load set-point. Interestingly, because the viral kinetics in the various treated animals differed-some animals exhibited large fluctuations in viral load after cART cessation-the model suggests that anti-α4β7 treatment could act by different primary mechanisms in different animals and still lead to post-treatment viral control. This outcome is nonetheless in accordance with a model with two stable viral load set-points, in which therapy can perturb the system from one set-point to a lower one through different biological mechanisms.

In order to assess the efficacy of novel HIV-1 treatments leading to a functional cure, the time to viral rebound is frequently used as a surrogate endpoint. The longer the time to viral rebound, the more efficacious the therapy. In support of such an approach, mathematical models serve as a connection between the size of the latent reservoir and the time to HIV-1 rebound after treatment interruption. The simplest of such models assumes that a single successful latent cell reactivation event leads to observable viremia after a period of exponential viral growth. Here we consider a generalization developed by Pinkevych et al. and Hill et al. of this simple model in which multiple reactivation events can occur, each contributing to the exponential growth of the viral load. We formalize and improve the previous derivation of the dynamics predicted by this model, and use the model to estimate relevant biological parameters from SIV rebound data. We confirm a previously described effect of very early antiretroviral therapy (ART) initiation on the rate of recrudescence and the viral load growth rate after treatment interruption. We find that every day ART initiation is delayed results in a 39% increase in the recrudescence rate (95% credible interval: [18%, 62%]), and a 11% decrease of the viral growth rate (95% credible interval: [4%, 20%]). We show that when viral rebound occurs early relative to the viral load doubling time, a model with multiple successful reactivation events fits the data better than a model with only a single successful reactivation event.

Mathematical modelling has successfully been used to provide quantitative descriptions of many viral infections, but for the Ebola virus, which requires biosafety level 4 facilities for experimentation, modelling can play a crucial role. Ebola virus modelling efforts have primarily focused on in vivo virus kinetics, e.g., in animal models, to aid the development of antivirals and vaccines. But, thus far, these studies have not yielded a detailed specification of the infection cycle, which could provide a foundational description of the virus kinetics and thus a deeper understanding of their clinical manifestation. Here, we obtain a diverse experimental data set of the Ebola virus infection in vitro, and then make use of Bayesian inference methods to fully identify parameters in a mathematical model of the infection. Our results provide insights into the distribution of time an infected cell spends in the eclipse phase (the period between infection and the start of virus production), as well as the rate at which infectious virions lose infectivity. We suggest how these results can be used in future models to describe co-infection with defective interfering particles, which are an emerging alternative therapeutic.

Hepatitis C virus (HCV) causes acute hepatitis C and can lead to life-threatening complications if it becomes chronic. The HCV genome is a single plus strand of RNA. Its intracellular replication is a spatiotemporally coordinated process of RNA translation upon cell infection, RNA synthesis within a replication compartment, and virus particle production. While HCV is mainly transmitted via mature infectious virus particles, it has also been suggested that HCV-infected cells can secrete HCV RNA carrying exosomes that can infect cells in a receptor independent manner. In order to gain insight into these two routes of transmission, we developed a series of intracellular HCV replication models that include HCV RNA secretion and/or virus assembly and release. Fitting our models to in vitro data, in which cells were infected with HCV, suggests that initially most secreted HCV RNA derives from intracellular cytosolic plus-strand RNA, but subsequently secreted HCV RNA derives equally from the cytoplasm and the replication compartments. Furthermore, our model fits to the data suggest that the rate of virus assembly and release is limited by host cell resources. Including the effects of direct acting antivirals in our models, we found that in spite of decreasing intracellular HCV RNA and extracellular virus concentration, low level HCV RNA secretion may continue as long as intracellular RNA is available. This may possibly explain the presence of detectable levels of plasma HCV RNA at the end of treatment even in patients that ultimately attain a sustained virologic response.

Antiretroviral therapy (ART) effectively controls HIV infection, suppressing HIV viral loads. Suspension of therapy is followed by rebound of viral loads to high, pre-therapy levels. However, there is significant heterogeneity in speed of rebound, with some rebounds occurring within days, weeks, or sometimes years. We present a stochastic mathematical model to gain insight into these post-treatment dynamics, specifically characterizing the dynamics of short term viral rebounds (≤ 60 days). Li et al. (2016) report that the size of the expressed HIV reservoir, i.e., cell-associated HIV RNA levels, and drug regimen correlate with the time between ART suspension and viral rebound to detectable levels. We incorporate this information and viral rebound times to parametrize our model. We then investigate insights offered by our model into the underlying dynamics of the latent reservoir. In particular, we refine previous estimates of viral recrudescence after ART interruption by accounting for heterogeneity in infection rebound dynamics, and determine a recrudescence rate of once every 2-4 days. Our parametrized model can be used to aid in design of clinical trials to study viral dynamics following analytic treatment interruption. We show how to derive informative personalized testing frequencies from our model and offer a proof-of-concept example. Our results represent first steps towards a model that can make predictions on a person living with HIV (PLWH)'s rebound time distribution based on biomarkers, and help identify PLWH with long viral rebound delays.

MOTIVATION: Illustrating how HIV-1 is transmitted and how it evolves in the following weeks is an important step for developing effective vaccination and prevention strategies. It is currently possible through DNA sequencing to account for the diverse array of viral strains within an infected individual. This provides an unprecedented opportunity to pinpoint when each patient was infected and which viruses were transmitted. RESULTS: Here we develop a mathematical tool for early HIV-1 evolution within a subject whose infection originates either from a single or multiple viral variants. The shifted Poisson mixture model (SPMM) provides a quantitative guideline for segregating viral lineages, which in turn enables us to assess when a subject was infected. The infection duration estimated by SPMM showed a statistically significant linear relationship with that by Fiebig laboratory staging (P = 0.00059) among 37 acutely infected subjects. Our tool provides a functional approach to understanding early genetic diversity, one of the most important parameters for deciphering HIV-1 transmission and predicting the rate of disease progression. AVAILABILITY AND IMPLEMENTATION: SPMM, webserver, is available at http://www.hayounlee.org/web-tools.html. CONTACT: [email protected] SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Antiretroviral therapy (ART) effectively controls HIV infection, suppressing HIV viral loads. However, some residual virus remains, below the level of detection, in HIV-infected patients on ART. The source of this viremia is an area of debate: does it derive primarily from activation of infected cells in the latent reservoir, or from ongoing viral replication? Observations seem to be contradictory: there is evidence of short term evolution, implying that there must be ongoing viral replication, and viral strains should thus evolve. However, phylogenetic analyses, and rare emergent drug resistance, suggest no long-term viral evolution, implying that virus derived from activated latent cells must dominate. We use simple deterministic and stochastic models to gain insight into residual viremia dynamics in HIV-infected patients. Our modeling relies on two underlying assumptions for patients on suppressive ART: that latent cell activation drives viral dynamics and that the reproductive ratio of treated infection is less than 1. Nonetheless, the contribution of viral replication to residual viremia in patients on ART may be non-negligible. However, even if the portion of viremia attributable to viral replication is significant, our model predicts (1) that latent reservoir re-seeding remains negligible, and (2) some short-term viral evolution is permitted, but long-term evolution can still be limited: stochastic analysis of our model shows that de novo emergence of drug resistance is rare. Thus, our simple models reconcile the seemingly contradictory observations on residual viremia and, with relatively few parameters, recapitulates HIV viral dynamics observed in patients on suppressive therapy.

Complications of HIV-1 infection in individuals who utilize drugs of abuse is a significant problem, because these drugs have been associated with higher virus replication and accelerated disease progression as well as severe neuropathogenesis. To gain further insight it is important to quantify the effects of drugs of abuse on HIV-1 infection dynamics. Here, we develop a mathematical model that incorporates experimentally observed effects of morphine on inducing HIV-1 co-receptor expression. For comparison we also considered viral dynamic models with cytolytic or noncytolytic effector cell responses. Based on the small sample size Akaike information criterion, these models were inferior to the new model based on changes in co-receptor expression. The model with morphine affecting co-receptor expression agrees well with the experimental data from simian immunodeficiency virus infections in morphine-addicted macaques. Our results show that morphine promotes a target cell subpopulation switch from a lower level of susceptibility to a state that is about 2-orders of magnitude higher in susceptibility to SIV infection. As a result, the proportion of target cells with higher susceptibility remains extremely high in morphine conditioning. Such a morphine-induced population switch not only has adverse effects on the replication rate, but also results in a higher steady state viral load and larger CD4 count drops. Moreover, morphine conditioning may pose extra obstacles to controlling viral load during antiretroviral therapy, such as pre-exposure prophylaxis and post infection treatments. This study provides, for the first time, a viral dynamics model, viral dynamics parameters, and related analytical and simulation results for SIV dynamics under drugs of abuse.

HIV-1 is subject to immune pressure exerted by the host, giving variants that escape the immune response an advantage. Virus released from activated latent cells competes against variants that have continually evolved and adapted to host immune pressure. Nevertheless, there is increasing evidence that virus displaying a signal of latency survives in patient plasma despite having reduced fitness due to long-term immune memory. We investigated the survival of virus with latent envelope genomic fragments by simulating within-host HIV-1 sequence evolution and the cycling of viral lineages in and out of the latent reservoir. Our model incorporates a detailed mutation process including nucleotide substitution, recombination, latent reservoir dynamics, diversifying selection pressure driven by the immune response, and purifying selection pressure asserted by deleterious mutations. We evaluated the ability of our model to capture sequence evolution in vivo by comparing our simulated sequences to HIV-1 envelope sequence data from 16 HIV-infected untreated patients. Empirical sequence divergence and diversity measures were qualitatively and quantitatively similar to those of our simulated HIV-1 populations, suggesting that our model invokes realistic trends of HIV-1 genetic evolution. Moreover, reconstructed phylogenies of simulated and patient HIV-1 populations showed similar topological structures. Our simulation results suggest that recombination is a key mechanism facilitating the persistence of virus with latent envelope genomic fragments in the productively infected cell population. Recombination increased the survival probability of latent virus forms approximately 13-fold. Prevalence of virus with latent fragments in productively infected cells was observed in only 2% of simulations when we ignored recombination, while the proportion increased to 27% of simulations when we allowed recombination. We also found that the selection pressures exerted by different fitness landscapes influenced the shape of phylogenies, diversity trends, and survival of virus with latent genomic fragments. Our model predicts that the persistence of latent genomic fragments from multiple different ancestral origins increases sequence diversity in plasma for reasonable fitness landscapes.

HIV-1-infected cells in peripheral blood can be grouped into different transcriptional subclasses. Quantifying the turnover of these cellular subclasses can provide important insights into the viral life cycle and the generation and maintenance of latently infected cells. We used previously published data from five patients chronically infected with HIV-1 that initiated combination antiretroviral therapy (cART). Patient-matched PCR for unspliced and multiply spliced viral RNAs combined with limiting dilution analysis provided measurements of transcriptional profiles at the single cell level. Furthermore, measurement of intracellular transcripts and extracellular virion-enclosed HIV-1 RNA allowed us to distinguish productive from non-productive cells. We developed a mathematical model describing the dynamics of plasma virus and the transcriptional subclasses of HIV-1-infected cells. Fitting the model to the data allowed us to better understand the phenotype of different transcriptional subclasses and their contribution to the overall turnover of HIV-1 before and during cART. The average number of virus-producing cells in peripheral blood is small during chronic infection. We find that a substantial fraction of cells can become defectively infected. Assuming that the infection is homogenous throughout the body, we estimate an average in vivo viral burst size on the order of 104 virions per cell. Our study provides novel quantitative insights into the turnover and development of different subclasses of HIV-1-infected cells, and indicates that cells containing solely unspliced viral RNA are a good marker for viral latency. The model illustrates how the pool of latently infected cells becomes rapidly established during the first months of acute infection and continues to increase slowly during the first years of chronic infection. Having a detailed understanding of this process will be useful for the evaluation of viral eradication strategies that aim to deplete the latent reservoir of HIV-1.

Several studies have proven oseltamivir to be efficient in reducing influenza viral titer and symptom intensity. However, the usefulness of oseltamivir can be compromised by the emergence and spread of drug-resistant virus. The selective pressure exerted by different oseltamivir therapy regimens have received little attention. Combining models of drug pharmacokinetics, pharmacodynamics, viral kinetics and symptom dynamics, we explored the efficacy of oseltamivir in reducing both symptoms (symptom efficacy) and viral load (virological efficacy). We simulated samples of 1000 subjects using previously estimated between-subject variability in viral and symptom dynamic parameters to describe the observed heterogeneity in a patient population. We simulated random mutations conferring resistance to oseltamivir. We explored the effect of therapy initiation time, dose, intake frequency and therapy duration on influenza infection, illness dynamics, and emergence of viral resistance. Symptom and virological efficacies were strongly associated with therapy initiation time. The proportion of subjects shedding resistant virus was 27-fold higher when prophylaxis was initiated during the incubation period compared with no treatment. It fell to below 1% when treatment was initiated after symptom onset for twice-a-day intakes. Lower doses and prophylaxis regimens led to lower efficacies and increased risk of resistance emergence. We conclude that prophylaxis initiated during the incubation period is the main factor leading to resistance emergence.

Hepatitis B is a DNA virus that infects liver cells and can cause both acute and chronic disease. It is believed that both viral and host factors are responsible for determining whether the infection is cleared or becomes chronic. Here we investigate the mechanism of protection by developing a mathematical model of the antibody response following hepatitis B virus (HBV) infection. We fitted the model to data from seven infected adults identified during acute infection and determined the ability of the virus to escape neutralization through overproduction of non-infectious subviral particles, which have HBs proteins on their surface, but do not contain nucleocapsid protein and viral nucleic acids. We showed that viral clearance can be achieved for high anti-HBV antibody levels, as in vaccinated individuals, when: (1) the rate of synthesis of hepatitis B subviral particles is slow; (2) the rate of synthesis of hepatitis B subviral particles is high but either anti-HBV antibody production is fast, the antibody affinity is high, or the levels of pre-existent HBV-specific antibody at the time of infection are high, as could be attained by vaccination. We further showed that viral clearance can be achieved for low equilibrium anti-HBV antibody levels, as in unvaccinated individuals, when a strong cellular immune response controls early infection.

Simple models of therapy for viral diseases such as hepatitis C virus (HCV) or human immunodeficiency virus assume that, once therapy is started, the drug has a constant effectiveness. More realistic models have assumed either that the drug effectiveness depends on the drug concentration or that the effectiveness varies over time. Here a previously introduced varying-effectiveness (VE) model is studied mathematically in the context of HCV infection. We show that while the model is linear, it has no closed-form solution due to the time-varying nature of the effectiveness. We then show that the model can be transformed into a Bessel equation and derive an analytic solution in terms of modified Bessel functions, which are defined as infinite series, with time-varying arguments. Fitting the solution to data from HCV infected patients under therapy has yielded values for the parameters in the model. We show that for biologically realistic parameters, the predicted viral decay on therapy is generally biphasic and resembles that predicted by constant-effectiveness (CE) models. We introduce a general method for determining the time at which the transition between decay phases occurs based on calculating the point of maximum curvature of the viral decay curve. For the parameter regimes of interest, we also find approximate solutions for the VE model and establish the asymptotic behavior of the system. We show that the rate of second phase decay is determined by the death rate of infected cells multiplied by the maximum effectiveness of therapy, whereas the rate of first phase decline depends on multiple parameters including the rate of increase of drug effectiveness with time.

Chronic liver infection by hepatitis C virus (HCV) is a major public health concern. Despite partly successful treatment options, several aspects of intrahepatic HCV infection dynamics are still poorly understood, including the preferred mode of viral propagation, as well as the proportion of infected hepatocytes. Answers to these questions have important implications for the development of therapeutic interventions. In this study, we present methods to analyze the spatial distribution of infected hepatocytes obtained by single cell laser capture microdissection from liver biopsy samples of patients chronically infected with HCV. By characterizing the internal structure of clusters of infected cells, we are able to evaluate hypotheses about intrahepatic infection dynamics. We found that individual clusters on biopsy samples range in size from 4-50 infected cells. In addition, the HCV RNA content in a cluster declines from the cell that presumably founded the cluster to cells at the maximal cluster extension. These observations support the idea that HCV infection in the liver is seeded randomly (e.g. from the blood) and then spreads locally. Assuming that the amount of intracellular HCV RNA is a proxy for how long a cell has been infected, we estimate based on models of intracellular HCV RNA replication and accumulation that cells in clusters have been infected on average for less than a week. Further, we do not find a relationship between the cluster size and the estimated cluster expansion time. Our method represents a novel approach to make inferences about infection dynamics in solid tissues from static spatial data.

The current paradigm for studying hepatitis C virus (HCV) dynamics in patients utilizes a standard viral dynamic model that keeps track of uninfected (target) cells, infected cells, and virus. The model does not account for the dynamics of intracellular viral replication, which is the major target of direct-acting antiviral agents (DAAs). Here we describe and study a recently developed multiscale age-structured model that explicitly considers the potential effects of DAAs on intracellular viral RNA production, degradation, and secretion as virus into the circulation. We show that when therapy significantly blocks both intracellular viral RNA production and virus secretion, the serum viral load decline has three phases, with slopes reflecting the rate of serum viral clearance, the rate of loss of intracellular viral RNA, and the rate of loss of intracellular replication templates and infected cells, respectively. We also derive analytical approximations of the multiscale model and use one of them to analyze data from patients treated for 14 days with the HCV protease inhibitor danoprevir. Analysis suggests that danoprevir significantly blocks intracellular viral production (with mean effectiveness 99.2%), enhances intracellular viral RNA degradation about 5-fold, and moderately inhibits viral secretion (with mean effectiveness 56%). The multiscale model can be used to study viral dynamics in patients treated with other DAAs and explore their mechanisms of action in treatment of hepatitis C.

Influenza virus infection remains a public health problem worldwide. The mechanisms underlying viral control during an uncomplicated influenza virus infection are not fully understood. Here, we developed a mathematical model including both innate and adaptive immune responses to study the within-host dynamics of equine influenza virus infection in horses. By comparing modeling predictions with both interferon and viral kinetic data, we examined the relative roles of target cell availability, and innate and adaptive immune responses in controlling the virus. Our results show that the rapid and substantial viral decline (about 2 to 4 logs within 1 day) after the peak can be explained by the killing of infected cells mediated by interferon activated cells, such as natural killer cells, during the innate immune response. After the viral load declines to a lower level, the loss of interferon-induced antiviral effect and an increased availability of target cells due to loss of the antiviral state can explain the observed short phase of viral plateau in which the viral level remains unchanged or even experiences a minor second peak in some animals. An adaptive immune response is needed in our model to explain the eventual viral clearance. This study provides a quantitative understanding of the biological factors that can explain the viral and interferon kinetics during a typical influenza virus infection.

Viral production from infected cells can occur continuously or in a burst that generally kills the cell. For HIV infection, both modes of production have been suggested. Standard viral dynamic models formulated as sets of ordinary differential equations can not distinguish between these two modes of viral production, as the predicted dynamics is identical as long as infected cells produce the same total number of virions over their lifespan. Here we show that in stochastic models of viral infection the two modes of viral production yield different early term dynamics. Further, we analytically determine the probability that infections initiated with any number of virions and infected cells reach extinction, the state when both the population of virions and infected cells vanish, and show this too has different solutions for continuous and burst production. We also compute the distributions of times to establish infection as well as the distribution of times to extinction starting from both a single virion as well as from a single infected cell for both modes of virion production.

Relatively little is known about the viral factors contributing to the lethality of the 1918 pandemic, although its unparalleled virulence was likely due in part to the newly discovered PB1-F2 protein. This protein, while unnecessary for replication, increases apoptosis in monocytes, alters viral polymerase activity in vitro, enhances inflammation and increases secondary pneumonia in vivo. However, the effects the PB1-F2 protein have in vivo remain unclear. To address the mechanisms involved, we intranasally infected groups of mice with either influenza A virus PR8 or a genetically engineered virus that expresses the 1918 PB1-F2 protein on a PR8 background, PR8-PB1-F2(1918). Mice inoculated with PR8 had viral concentrations peaking at 72 hours, while those infected with PR8-PB1-F2(1918) reached peak concentrations earlier, 48 hours. Mice given PR8-PB1-F2(1918) also showed a faster decline in viral loads. We fit a mathematical model to these data to estimate parameter values. The model supports a higher viral production rate per cell and a higher infected cell death rate with the PR8-PB1-F2(1918) virus. We discuss the implications these mechanisms have during an infection with a virus expressing a virulent PB1-F2 on the possibility of a pandemic and on the importance of antiviral treatments.

BACKGROUND: The occurrence of a genetic bottleneck in HIV sexual or mother-to-infant transmission has been well documented. This results in a majority of new infections being homogeneous, i.e., initiated by a single genetic strain. Early after infection, prior to the onset of the host immune response, the viral population grows exponentially. In this simple setting, an approach for estimating evolutionary and demographic parameters based on comparison of diversity measures is a feasible alternative to the existing Bayesian methods (e.g., BEAST), which are instead based on the simulation of genealogies. RESULTS: We have devised a web tool that analyzes genetic diversity in acutely infected HIV-1 patients by comparing it to a model of neutral growth. More specifically, we consider a homogeneous infection (i.e., initiated by a unique genetic strain) prior to the onset of host-induced selection, where we can assume a random accumulation of mutations. Previously, we have shown that such a model successfully describes about 80% of sexual HIV-1 transmissions provided the samples are drawn early enough in the infection. Violation of the model is an indicator of either heterogeneous infections or the initiation of selection. CONCLUSIONS: When the underlying assumptions of our model (homogeneous infection prior to selection and fast exponential growth) are met, we are under a very particular scenario for which we can use a forward approach (instead of backwards in time as provided by coalescent methods). This allows for more computationally efficient methods to derive the time since the most recent common ancestor. Furthermore, the tool performs statistical tests on the Hamming distance frequency distribution, and outputs summary statistics (mean of the best fitting Poisson distribution, goodness of fit p-value, etc). The tool runs within minutes and can readily accommodate the tens of thousands of sequences generated through new ultradeep pyrosequencing technologies. The tool is available on the LANL website.

It has recently been estimated that a single HIV-1 infected cell produces between and more than viral particles over its life span. Since body-wide estimates of the ratio of free virus to productively infected cells are smaller than and much smaller than , individual virions must be cleared rapidly. This seems difficult to reconcile with the fact that most of the total body virus is trapped on follicular dendritic cells where it can survive for many months. It has also been difficult to reconcile the vast difference in the rates at which the virus is cleared from the blood in rhesus macaques and in chronically infected patients. Here we attempt to reconcile these seemingly contradictory observations by considering the virion clearance rate in various organs and the virion exchange rates between them. The main results are that the per capita clearance rate of free virus in lymphoid tissue should be fast, the virion exchange rate between lymphoid tissue and the blood should be slow, and the comparatively slow previous estimates for the virion clearance rate from the blood correspond to the rate of virion efflux from the blood to other organs where the virus is ultimately cleared.

For most HIV-infected patients, antiretroviral therapy controls viral replication. However, in some patients drug resistance can cause therapy to fail. Nonetheless, continued therapy with a failing regimen can preserve or even lead to increases in CD4+ T cell counts. To understand the biological basis of these observations, we used mathematical models to explain observations made in patients with drug-resistant HIV treated with enfuvirtide (ENF/T-20), an HIV-1 fusion inhibitor. Due to resistance emergence, ENF was removed from the drug regimen, drug-sensitive virus regrown, and ENF was re-administered. We used our model to study the dynamics of plasma-viral RNA and CD4+ T cell levels, and the competition between drug-sensitive and resistant viruses during therapy interruption and re-administration. Focusing on resistant viruses carrying the V38A mutation in gp41, we found ENF-resistant virus to be 17±3% less fit than ENF-sensitive virus in the absence of the drug, and that the loss of resistant virus during therapy interruption was primarily due to this fitness cost. Using viral dynamic parameters estimated from these patients, we show that although re-administration of ENF cannot suppress viral load, it can, in the presence of resistant virus, increase CD4+ T cell counts, which should yield clinical benefits. This study provides a framework to investigate HIV and T cell dynamics in patients who develop drug resistance to other antiretroviral agents and may help to develop more effective strategies for treatment.

Although potent combination therapy is usually able to suppress plasma viral loads in HIV-1 patients to below the detection limit of conventional clinical assays, a low level of viremia frequently can be detected in plasma by more sensitive assays. Additionally, many patients experience transient episodes of viremia above the detection limit, termed viral blips, even after being on highly suppressive therapy for many years. An obstacle to viral eradication is the persistence of a latent reservoir for HIV-1 in resting memory CD4(+) T cells. The mechanisms underlying low viral load persistence, slow decay of the latent reservoir, and intermittent viral blips are not fully characterized. The quantitative contributions of residual viral replication to viral and the latent reservoir persistence remain unclear. In this paper, we probe these issues by developing a mathematical model that considers latently infected cell activation in response to stochastic antigenic stimulation. We demonstrate that programmed expansion and contraction of latently infected cells upon immune activation can generate both low-level persistent viremia and intermittent viral blips. Also, a small fraction of activated T cells revert to latency, providing a potential to replenish the latent reservoir. By this means, occasional activation of latently infected cells can explain the variable decay characteristics of the latent reservoir observed in different clinical studies. Finally, we propose a phenomenological model that includes a logistic term representing homeostatic proliferation of latently infected cells. The model is simple but can robustly generate the multiphasic viral decline seen after initiation of therapy, as well as low-level persistent viremia and intermittent HIV-1 blips. Using these models, we provide a quantitative and integrated prospective into the long-term dynamics of HIV-1 and the latent reservoir in the setting of potent antiretroviral therapy.

Quantifying the dynamics of intrahost HIV-1 sequence evolution is one means of uncovering information about the interaction between HIV-1 and the host immune system. In the chronic phase of infection, common dynamics of sequence divergence and diversity have been reported. We developed an HIV-1 sequence evolution model that simulated the effects of mutation and fitness of sequence variants. The amount of evolution was described by the distance from the founder strain, and fitness was described by the number of offspring a parent sequence produces. Analysis of the model suggested that the previously observed saturation of divergence and decrease of diversity in later stages of infection can be explained by a decrease in the proportion of offspring that are mutants as the distance from the founder strain increases rather than due to an increase of viral fitness. The prediction of the model was examined by performing phylogenetic analysis to estimate the change in the rate of evolution during infection. In agreement with our modeling, in 13 out of 15 patients (followed for 3-12 years) we found that the rate of intrahost HIV-1 evolution was not constant but rather slowed down at a rate correlated with the rate of CD4+ T-cell decline. The correlation between the dynamics of the evolutionary rate and the rate of CD4+ T-cell decline, coupled with our HIV-1 sequence evolution model, explains previously conflicting observations of the relationships between the rate of HIV-1 quasispecies evolution and disease progression.

Despite many years of potent antiretroviral therapy, latently infected cells and low levels of plasma virus have been found to persist in HIV-infected patients. The factors influencing this persistence and their relative contributions have not been fully elucidated and remain controversial. Here, we address these issues by developing and employing a simple, but mechanistic viral dynamics model. The model has two novel features. First, it assumes that latently infected T cells can undergo bystander proliferation without transitioning into active viral production. Second, it assumes that the rate of latent cell activation decreases with time on antiretroviral therapy due to the activation and subsequent loss of latently infected cells specific for common antigens, leaving behind cells that are successively less frequently activated. Using the model, we examined the quantitative contributions of T cell bystander proliferation, latent cell activation, and ongoing viral replication to the stability of the latent reservoir and persisting low-level viremia. Not surprisingly, proliferation of latently infected cells helped maintain the latent reservoir in spite of loss of latent infected cells through activation and death, and affected viral dynamics to an extent that depended on the magnitude of latent cell activation. In the limit of zero latent cell activation, the latent cell pool and viral load became uncoupled. However, as the activation rate increased, the plasma viral load could be maintained without depleting the latent reservoir, even in the absence of viral replication. The influence of ongoing viral replication on the latent reservoir remained insignificant for drug efficacies above the "critical efficacy" irrespective of the activation rate. However, for lower drug efficacies viral replication enabled the stable maintenance of both the latent reservoir and the virus. Our model and analysis methods provide a quantitative and qualitative framework for probing how different viral and host factors contribute to the dynamics of the latent reservoir and the virus, offering new insights into the principal determinants of their persistence.

The problem of protecting computer systems can be viewed generally as the problem of learning to distinguish self from other. The authors describe a method for change detection which is based on the generation of T cells in the immune system. Mathematical analysis reveals computational costs of the system, and preliminary experiments illustrate how the method might be applied to the problem of computer viruses.>

This paper describes an immune system model based on binary strings. The purpose of the model is to study the pattern-recognition processes and learning that take place at both the individual and species levels in the immune system. The genetic algorithm (GA) is a central component of the model. The paper reports simulation experiments on two pattern-recognition problems that are relevant to natural immune systems. Finally, it reviews the relation between the model and explicit fitness-sharing techniques for genetic algorithms, showing that the immune system model implements a form of implicit fitness sharing.

In typical applications, genetic algorithms (GAs) process populations of potential problem solutions to evolve a single population member that specifies an ‘optimized’ solution. The majority of GA analysis has focused on these optimization applications. In other applications (notably learning classifier systems and certain connectionist learning systems), a GA searches for a population of cooperative structures that jointly perform a computational task. This paper presents an analysis of this type of GA problem. The analysis considers a simplified genetics-based machine learning system: a model of an immune system. In this model, a GA must discover a set of pattern-matching antibodies that effectively match a set of antigen patterns. Analysis shows how a GA can automatically evolve and sustain a diverse, cooperative population. The cooperation emerges as a natural part of the antigen-antibody matching procedure. This emergent effect is shown to be similar to fitness sharing, an explicit technique for multimodal GA optimization. Further analysis shows how the GA population can adapt to express various degrees of generalization. The results show how GAs can automatically and simultaneously discover effective groups of cooperative computational structures.

Multiport Systems and Bond Graphs Basic Component Models System Models State-Space Equations Analysis of Linear Systems Multiport Fields and Junction Structures Transducers, Amplifiers, and Instruments Rigid-Body Mechanics Distributed Parameter Systems Magnetic Circuits and Devices Thermofluid Systems Nonlinear Simulation Index.