The cortical microtubules of Toxoplasma gondii underlie the helicity of parasite movement.

Motility is essential for apicomplexan parasites to infect their hosts. In a three-dimensional (3D) environment, the apicomplexan parasite Toxoplasma gondii moves along a helical path. The cortical microtubules, which are ultra-stable and spirally arranged, have been considered to be a structure that guides the long-distance movement of the parasite. Here, we address the role of the cortical microtubules in parasite motility, invasion and egress by utilizing a previously generated mutant (dubbed 'TKO') in which these microtubules are destabilized in mature parasites. We found that the cortical microtubules in ∼80% of the non-dividing (i.e. daughter-free) TKO parasites are much shorter than normal. The extent of depolymerization was further exacerbated upon commencement of daughter formation or cold treatment, but parasite replication was not affected. In a 3D Matrigel matrix, the TKO mutant moved directionally over long distances, but along trajectories that were significantly more linear (i.e. less helical) than those of wild-type parasites. Interestingly, this change in trajectory did not impact either movement speed in the matrix or the speed and behavior of the parasite during entry into and egress from the host cell.

ENTRAIN: integrating trajectory inference and gene regulatory networks with spatial data to co-localize the receptor-ligand interactions that specify cell fate.

MOTIVATION: Cell fate is commonly studied by profiling the gene expression of single cells to infer developmental trajectories based on expression similarity, RNA velocity, or statistical mechanical properties. However, current approaches do not recover microenvironmental signals from the cellular niche that drive a differentiation trajectory. RESULTS: We resolve this with environment-aware trajectory inference (ENTRAIN), a computational method that integrates trajectory inference methods with ligand-receptor pair gene regulatory networks to identify extracellular signals and evaluate their relative contribution towards a differentiation trajectory. The output from ENTRAIN can be superimposed on spatial data to co-localize cells and molecules in space and time to map cell fate potentials to cell-cell interactions. We validate and benchmark our approach on single-cell bone marrow and spatially resolved embryonic neurogenesis datasets to identify known and novel environmental drivers of cellular differentiation. AVAILABILITY AND IMPLEMENTATION: ENTRAIN is available as a public package at https://github.com/theimagelab/entrain and can be used on both single-cell and spatially resolved datasets.

An apical protein, Pcr2, is required for persistent movement by the human parasite Toxoplasma gondii.

The phylum Apicomplexa includes thousands of species of unicellular parasites that cause a wide range of human and animal diseases such as malaria and toxoplasmosis. To infect, the parasite must first initiate active movement to disseminate through tissue and invade into a host cell, and then cease moving once inside. The parasite moves by gliding on a surface, propelled by an internal cortical actomyosin-based motility apparatus. One of the most effective invaders in Apicomplexa is Toxoplasma gondii, which can infect any nucleated cell and any warm-blooded animal. During invasion, the parasite first makes contact with the host cell "head-on" with the apical complex, which features an elaborate cytoskeletal apparatus and associated structures. Here we report the identification and characterization of a new component of the apical complex, Preconoidal region protein 2 (Pcr2). Pcr2 knockout parasites replicate normally, but they are severely diminished in their capacity for host tissue destruction due to significantly impaired invasion and egress, two vital steps in the lytic cycle. When stimulated for calcium-induced egress, Pcr2 knockout parasites become active, and secrete effectors to lyse the host cell. Calcium-induced secretion of the major adhesin, MIC2, also appears to be normal. However, the movement of the Pcr2 knockout parasite is spasmodic, which drastically compromises egress. In addition to faulty motility, the ability of the Pcr2 knockout parasite to assemble the moving junction is impaired. Both defects likely contribute to the poor efficiency of invasion. Interestingly, actomyosin activity, as indicated by the motion of mEmerald tagged actin chromobody, appears to be largely unperturbed by the loss of Pcr2, raising the possibility that Pcr2 may act downstream of or in parallel with the actomyosin machinery.

Dynamics of latent HIV under clonal expansion.

The HIV latent reservoir exhibits slow decay on antiretroviral therapy (ART), impacted by homeostatic proliferation and activation. How these processes contribute to the total dynamic while also producing the observed profile of sampled latent clone sizes is unclear. An agent-based model was developed that tracks individual latent clones, incorporating homeostatic proliferation of cells and activation of clones. The model was calibrated to produce observed latent reservoir dynamics as well as observed clonal size profiles. Simulations were compared to previously published latent HIV integration data from 5 adults and 3 children. The model simulations reproduced reservoir dynamics as well as generating residual plasma viremia levels (pVL) consistent with observations on ART. Over 382 Latin Hypercube Sample simulations, the median latent reservoir grew by only 0.3 log10 over the 10 years prior to ART initiation, after which time it decreased with a half-life of 15 years, despite number of clones decreasing at a faster rate. Activation produced a maximum size of genetically intact clones of around one million cells. The individual simulation that best reproduced the sampled clone profile, produced a reservoir that decayed with a 13.9 year half-life and where pVL, produced mainly from proliferation, decayed with a half-life of 10.8 years. These slow decay rates were achieved with mean cell life-spans of only 14.2 months, due to expansion of the reservoir through proliferation and activation. Although the reservoir decayed on ART, a number of clones increased in size more than 4,000-fold. While small sampled clones may have expanded through proliferation, the large sizes exclusively arose from activation. Simulations where homeostatic proliferation contributed more to pVL than activation, produced pVL that was less variable over time and exhibited fewer viral blips. While homeostatic proliferation adds to the latent reservoir, activation can both add and remove latent cells. Latent activation can produce large clones, where these may have been seeded much earlier than when first sampled. Elimination of the reservoir is complicated by expanding clones whose dynamic differ considerably to that of the entire reservoir.

NGlyAlign: an automated library building tool to align highly divergent HIV envelope sequences.

BACKGROUND: The high variability in envelope regions of some viruses such as HIV allow the virus to establish infection and to escape subsequent immune surveillance. This variability, as well as increasing incorporation of N-linked glycosylation sites, is fundamental to this evasion. It also creates difficulties for multiple sequence alignment methods (MSA) that provide the first step in their analysis. Existing MSA tools often fail to properly align highly variable HIV envelope sequences requiring extensive manual editing that is impractical with even a moderate number of these variable sequences. RESULTS: We developed an automated library building tool NGlyAlign, that organizes similar N-linked glycosylation sites as block constraints and statistically conserved global sites as single site constraints to automatically enforce partial columns in consistency-based MSA methods such as Dialign. This combined method accurately aligns variable HIV-1 envelope sequences. We tested the method on two datasets: a set of 156 founder and chronic gp160 HIV-1 subtype B sequences as well as a set of reference sequences of gp120 in the highly variable region 1. On measures such as entropy scores, sum of pair scores, column score, and similarity heat maps, NGlyAlign+Dialign proved superior against methods such as T-Coffee, ClustalOmega, ClustalW, Praline, HIValign and Muscle. The method is scalable to large sequence sets producing accurate alignments without requiring manual editing. As well as this application to HIV, our method can be used for other highly variable glycoproteins such as hepatitis C virus envelope. CONCLUSIONS: NGlyAlign is an automated tool for mapping and building glycosylation motif libraries to accurately align highly variable regions in HIV sequences. It can provide the basis for many studies reliant on single robust alignments. NGlyAlign has been developed as an open-source tool and is freely available at https://github.com/UNSW-Mathematical-Biology/NGlyAlign_v1.0 .

Elimination of cervical cancer in Tanzania: Modelled analysis of elimination in the context of endemic HIV infection and active HIV control.

The World Health Organisation (WHO) has launched a strategic initiative for cervical cancer (CC) elimination which involves scaling up three interventions: human papillomavirus (HPV) vaccination, twice-lifetime HPV-screening screening and pre-cancer/cancer treatment by 2030. CC is challenging to control in countries with endemic human immunodeficiency virus (HIV), as women living with HIV (WLHIV) are at elevated risk of HPV infection, persistence and progression. This analysis estimated the impact of the elimination interventions on CC incidence and mortality but additionally considered more intensive screening for WLHIV, using Tanzania as an example. A dynamic HIV/HPV model was used to simulate the elimination strategy for vaccination, screening and pre-cancer/cancer treatment, with 3-yearly HPV-screening in WLHIV starting at age 25 years, in the context of sustained HIV control in Tanzania from 2020 to 2119. Without vaccination or HPV screening, CC incidence rates per 100 000 women are predicted to fall from 58.0 in 2020 to 41.6 (range: 39.1-44.7) in 2119, due to existing HIV control. HPV vaccination and twice-lifetime HPV-screening for the general population and 3-yearly screening for WLHIV, would reduce CC incidence to 1.3 (range: 1.3-2.5) by 2119, with elimination (<4/100 000) in 2076 (range: 2076-2092). CC mortality rates per 100 000 women are predicted to reach 1.1 (range: 1.1-2.1) with further reductions contingent on increased CC treatment access. Vaccination and 3-yearly HPV-screening for WLHIV is predicted to achieve elimination in the subgroup of WLHIV potentially as early as 2061 (range: 2061-2078), with a 2119 CC incidence rate of 1.7 (range: 1.7-3.3). Scaling-up vaccination and HPV-screening will substantially reduce CC incidence in Tanzania, with elimination predicted within a century. Three-yearly HPV-screening and HPV vaccination, at high coverage rates, would facilitate CC elimination among WLHIV, and thus accelerate elimination in the overall population.

Variability in molecular characteristics of Hepatitis E virus quasispecies could modify viral surface properties and transmission.

Hepatitis E virus (HEV) usually causes self-limited liver diseases but can also result in severe cases. Genotypes 1 (G1) and 2 circulate in developing countries are human-restricted and waterborne, while zoonotic G3 and G4 circulating in industrialized countries preferentially infect human through consumption of contaminated meat. Our aims were to identify amino acid patterns in HEV variants that could be involved in pathogenicity or in transmission modes, related to their impact on antigenicity and viral surface hydrophobicity. HEV sequences from human (n = 37) and environmental origins (wild boar [n = 3], pig slaughterhouse effluent [n = 6] and urban wastewater [n = 2]) were collected for the characterization of quasispecies using ultra-deep sequencing (ORF2/ORF3 overlap). Predictive and functional assays were carried out to investigate viral particle antigenicity and hydrophobicity. Most quasispecies showed a major variant while a mixture was observed in urban wastewater and in one chronically infected patient. Amino acid signatures were identified, as a rabbit-linked HEV pattern in two infected patients, or the S68L (ORF2) / H81C (ORF3) residue mostly identified in wild boars. By comparison with environmental strains, molecular patterns less likely represented in humans were identified. Patterns impacting viral hydrophobicity and/or antigenicity were also observed, and the higher hydrophobicity of HEV naked particles compared with the enveloped forms was demonstrated. HEV variants isolated from human and environment present molecular patterns that could impact their surface properties as well as their transmission. These molecular patterns may concern only one minor variant of a quasispecies and could emerge under selective pressure.

High rates of cirrhosis and severe clinical events in patients with HBV/HDV co-infection: longitudinal analysis of a German cohort.

BACKGROUND: Chronic hepatitis delta virus (HDV) infection causes severe liver disease which often leads to cirrhosis and hepatocellular carcinoma (HCC). Aim of this study was to establish the disease severity and prognostic factors for disease outcome by analysing frequencies of clinical events and their correlation with baseline virological and biochemical parameters as well as interferon and nucleos(t)ide analogue treatment choice. METHODS: We studied a single-centre cohort of 49 anti-HDAg-positive patients with HBsAg persistence for at least 6 months. Virological and biochemical parameters, interferon and nucleos(t)ide analogue treatment choice as well as clinical events during follow-up were analysed by retrospective chart review (mean follow-up time 3 years, range 0.25-7.67 years). RESULTS: Severe clinical events occurred in 11/49 hepatitis D patients, including HCC (8/49), death (8/49) or liver transplantation (2/49). HCCs only occurred secondary to liver cirrhosis and their event rates in this cohort of hepatitis D patients did not differ from a matched HBV mono-infected cohort with comparable frequency of liver cirrhosis. A stepwise multivariate logistic regression revealed low platelet count (p = 0. 0290) and older age (p = 0.0337) correlating most strongly with overall clinical events, while serum HDV RNA positivity at baseline did not correlate with any clinical outcome. Interferon-free but not nucleos(t)ide analogue-free patient care correlated with the occurrence of HCC at logistic regression, although only 3/18 interferon-treated patients demonstrated repeatedly negative HDV PCR results post therapy. CONCLUSIONS: Our data indicate that progressive liver disease at baseline plays a major role as predictive factor for overall clinical outcome of hepatitis D patients. In particular, HCC risk may not be underestimated in hepatitis D virus RNA negative hepatitis D patients with advanced liver fibrosis.

Analysis and dissociation of anti-HIV effects of shRNA to CCR5 and the fusion inhibitor C46.

BACKGROUND: The gene therapeutic Cal-1 comprises the anti-HIV agents: (i) sh5, a short hairpin RNA to CCR5 that down-regulates CCR5 expression and (ii) maC46 (C46), a peptide that inhibits viral fusion with the cell membrane. These constructs were assessed for inhibition of viral replication and selective cell expansion in a number of settings. METHODS: HIV replication, selective outgrowth and cell surface viral binding were analysed with a single cycle infection assay of six pseudotyped HIV strains and a static and longitudinal passaging of MOLT4/CCR5 cells with HIV. Pronase digestion of surface virus and fluorescence microscopy assessed interactions between HIV virions and transduced cells. RESULTS: Cal-1 reduced CCR5 expression in peripheral blood mononuclear cells to CCR5Δ32 heterozygote levels. Even low level transduction resulted in significant preferential expansion in MOLT4/CCR5 gene-containing cells over a 3-week HIV challenge regardless of viral suppression [12.5% to 47.0% (C46), 46.7% (sh5), 62.2% (Dual), respectively]. The sh5 and Dual constructs at > 95% transduction also significantly suppressed virus to day 12 in the passage assay and all constructs, at varying percentage transduction inhibited virus in static culture. No escape mutations were present through 9 weeks of challenge. The Dual construct significantly suppressed infection by a panel of CCR5-using viruses, with its efficacy being independently determined from the single constructs. Dual and sh5 inhibited virion internalisation, as determined via pronase digestion of surface bound virus, by 70% compared to 13% for C46. CONCLUSIONS: The use of two anti-HIV genes allows optimal preferential survival and inhibition of HIV replication, with the impact on viral load being dependent on the percentage of gene marked cells.

Roadmap to control HBV and HDV epidemics in China.

PURPOSE: Hepatitis B virus (HBV) is endemic in China. Almost 10% of HBV infected individuals are also infected with hepatitis D virus (HDV) which has a 5-10 times higher mortality rate than HBV mono-infection. The aim of this manuscript is to devise strategies that can not only control HBV infections but also HDV infections in China under the current health care budget in an optimal manner. METHODS: Using a mathematical model, an annual budget of $10billion was optimally allocated among five interventions namely, testing and HBV adult vaccination, treatment for mono-infected and dually-infected individuals, second line treatment for HBV mono-infections, and awareness programs. RESULTS: We determine that the optimal strategy is to test and treat both infections as early as possible while applying awareness programs at full intensity. Under this strategy, an additional 19.8million HBV, 1.9million HDV infections and 0.25million lives will be saved over the next 10years at a cost-savings of $79billion than performing no intervention. Introduction of second line treatment does not add a significant economic burden yet prevents 1.4million new HBV infections and 15,000 new HDV infections. CONCLUSION: Test and treatment programs are highly efficient in reducing HBV and HDV prevalence in the population. Under the current health budget in China, not only test and treat programs but awareness programs and second line treatment can also be implemented that minimizes prevalence and mortality, and maximizes economic benefits.

An icosahedral virus as a fluorescent calibration standard: a method for counting protein molecules in cells by fluorescence microscopy.

The ability to replace genes coding for cellular proteins with DNA that codes for fluorescent protein-tagged versions opens the way to counting the number of molecules of each protein component of macromolecular assemblies in vivo by measuring fluorescence microscopically. Converting fluorescence to absolute numbers of molecules requires a fluorescent standard whose molecular composition is known precisely. In this report, the construction, properties and mode of using a set of fluorescence calibration standards are described. The standards are based on an icosahedral virus engineered to contain exactly 240 copies of one of seven different fluorescent proteins. Two applications of the fluorescent standards to counting molecules in the human parasite Toxoplasma gondii are described. Methods for improving the preciseness of the measurements and minimizing potential inaccuracies are emphasized.

Recognizing the impact of endemic hepatitis D virus on hepatitis B virus eradication.

BACKGROUND: Hepatitis delta virus (HDV) in conjunction with hepatitis B virus (HBV) increases adult morbidity and mortality. A number of studies have performed cost-benefit analyses for HBV interventions, but they have ignored the impact of HDV on these outcomes. METHODS: Using a mathematical model of HBV-HDV epidemiology, we compare health benefits and cost outcomes of four interventions: testing with HBV adult vaccination (diagnosis), diagnosis with antiviral treatment for HBV infections (mono-infections), diagnosis with antiviral treatment for HBV-HDV infections (dual-infections), and awareness programs. The relationship between optimal levels and outcomes of each of these interventions and HDV prevalence in HBV infected individuals ranging from 0 to 50% is determined. RESULTS: Over a 50 year period under no intervention, HBV prevalence, per capita total cost and death toll increase by 2.25%, -$11 and 2.6-fold respectively in moderate HDV endemic regions compared to mono-infected regions; the corresponding values for high HDV endemic regions are 4.2%, -$21 and 3.9-fold. Optimal interventions can be strategized similarly in mono and dually endemic regions. Only implementation of all four interventions achieves a very low HBV prevalence of around 1.5% in a moderate HDV endemic region such as China, with 2.8 million fewer deaths compared to no intervention. Although the policy of implementation of all four interventions costs additional $382 billion compared to no intervention, it still remains cost-effective with an incremental cost-effectiveness ratio of $1400/QALY. Very high efficacy awareness programs achieve less prevalence with fewer deaths at a lower cost compared to treatment and/or vaccination programs. CONCLUSION: HDV substantially affects the performance of any HBV-related intervention. Its exclusion results in over-estimation of the effectiveness of HBV interventions.

Dynamics of in vivo hepatitis D virus infection.

UNLABELLED: Hepatitis-D virus (HDV) is a satellite virus of hepatitis-B virus (HBV) whose intracellular products are required for the completion of the HDV life cycle. HDV can replicate in a cell without the presence of HBV but needs hepatitis B surface antigen (HBsAg) to complete virus assembly and packaging. In order to better understand HDV dynamics, we developed a mathematical model and successfully simulated HBV and HDV data under a range of scenarios. Compared to HBV mono-infection, dual HDV infection resulted in lower chronic HBV DNA levels, with more marked suppression for coinfection (1 logs HBV DNA copies/ml lower) compared to superinfection (0.6 logs HBV DNA copies/ml). Although they have no effect on HBV, prenylation inhibitors may provide the best therapy for reducing HDV viremia irrespective of the stage in which they are commenced. We found that highly effective long term pegylated interferon (IFN) therapy (99.99%) eliminates HBV and HDV viremia while less effective long term IFN therapy (99%) will only produce approximately 2.03 logs and no decrease in HBV and HDV viremia respectively in both coinfection and superinfection settings. Our study also suggests that there is a substantial difference in the outcome of therapies depending upon the time of commencement. CONCLUSION: Mathematical modeling of HDV infection can describe the complex interplay between this virus and HBV. Simulations suggest that HDV impacts on the feedback mechanisms that maintain cccDNA levels and that targeting these mechanisms may result in new therapeutic agents for both viruses.

Variability in long-term hepatitis B virus dynamics under antiviral therapy.

Hepatitis B virus (HBV) dynamics in treated patients can be complex and differ considerably from other viral infections. We analyse dynamics of liver and serum levels of HBV DNA in 24 chronically HBV-infected individuals undergoing 1 year of combination therapy with pegylated interferon alpha and adefovir dipivoxil (ADV), followed by 2 years of ADV monotherapy. Serum viral dynamics differentiated the patients into four response groups dependent on how quickly viremia became undetectable: quickly suppressed (HBV DNA <100 copies/ml within 8 weeks and staying suppressed, GRP1); quickly suppressed but some rebound (<10,000 copies/ml, GRP2); slow decay (GRP3); virological failures (>10,000 copies/ml, GRP4). These groups did not differ before start of therapy by serum HBV DNA (p=0.2), HBsAg (p=0.1), ALT (p=0.4), total HBV DNA within the liver (p=0.08), or cccDNA (p=0.3). Despite very different serum HBV DNA levels after 3 years, there was no statistical difference in total HBV DNA within the liver (p=0.08), nor in cccDNA levels (p=0.1), but HBsAg levels in serum were significantly lower for GRP1 compared to GRP4 (p=0.02). Efficacy in terms of reduction over the 3 years of serum HBV DNA, liver HBV DNA, cccDNA, and ratios of liver HBV DNA to cccDNA were 99.98%, 99.5%, 98.4%, and 83.2% respectively, exhibiting larger antiviral effects in serum than in liver. Over the course of therapy, HBV DNA viremia exhibited large oscillations for some individuals. Mathematical modelling reproduced the dynamics of these diverse groups by assuming a number of viral clones arose that experienced delayed recognition by the antibody response. Large viremia oscillations under therapy suggest sequential outgrowth of viral clones with delayed recognition by the humoral response.

Predicting first traversal times for virions and nanoparticles in mucus with slowed diffusion.

Particle-tracking experiments focusing on virions or nanoparticles in mucus have measured mean-square displacements and reported diffusion coefficients that are orders of magnitude smaller than the diffusion coefficients of such particles in water. Accurate description of this subdiffusion is important to properly estimate the likelihood of virions traversing the mucus boundary layer and infecting cells in the epithelium. However, there are several candidate models for diffusion that can fit experimental measurements of mean-square displacements. We show that these models yield very different estimates for the time taken for subdiffusive virions to traverse through a mucus layer. We explain why fits of subdiffusive mean-square displacements to standard diffusion models may be misleading. Relevant to human immunodeficiency virus infection, using computational methods for fractional subdiffusion, we show that subdiffusion in normal acidic mucus provides a more effective barrier against infection than previously thought. By contrast, the neutralization of the mucus by alkaline semen, after sexual intercourse, allows virions to cross the mucus layer and reach the epithelium in a short timeframe. The computed barrier protection from fractional subdiffusion is some orders of magnitude greater than that derived by fitting standard models of diffusion to subdiffusive data.

HIV DNA subspecies persist in both activated and resting memory CD4+ T cells during antiretroviral therapy.

UNLABELLED: The latent HIV reservoir is a major impediment to curing HIV infection. The contribution of CD4(+) T cell activation status to the establishment and maintenance of the latent reservoir was investigated by enumerating viral DNA components in a cohort of 12 individuals commencing antiretroviral therapy (ART) containing raltegravir, an integrase inhibitor. Prior to ART, the levels of total HIV DNA were similar across HLA-DR(+) and HLA-DR(-) (HLA-DR(±)) CD38(±) memory CD4(+) T cell phenotypes; episomal two-long terminal repeat (2-LTR) HIV DNA levels were higher in resting (HLA-DR(-) CD38(-)) cells, and this phenotype exhibited a significantly higher ratio of 2-LTR to integrated HIV DNA (P = 0.002). After 1 year of ART, there were no significant differences across each of the memory phenotypes of any HIV DNA component. The decay dynamics of integrated HIV DNA were slow within each subset, and integrated HIV DNA in the resting HLA-DR(-) CD38(-) subset per mm(3) of peripheral blood exhibited no significant decay (half-life of 25 years). Episomal 2-LTR HIV DNA decayed relative to integrated HIV DNA in resting cells with a half-life of 134 days. Surprisingly, from week 12 on, the decay rates of both total and episomal HIV DNA were lower in activated CD38(+) cells. By weeks 24 and 52, HIV RNA levels in plasma were most significantly correlated with the numbers of resting cells containing integrated HIV DNA. On the other hand, total HIV DNA levels in all subsets were significantly correlated with the numbers of HLA-DR(+) CD38(-) cells containing integrated HIV DNA. These results provide insights into the interrelatedness of cell activation and reservoir maintenance, with implications for the design of therapeutic strategies targeting HIV persistence. IMPORTANCE: It is generally believed that HIV is not cleared by extensive antiretroviral therapy (ART) due to the difficulty in eradicating the latent reservoir in resting CD4(+) T cells. New therapies that attempt to activate this reservoir so that immune or viral cytopathic mechanisms can remove those infected cells are currently being investigated. However, results obtained in this research indicate that activation, at least on some level, already occurs within this reservoir. Furthermore, we are the first to describe the dynamics of different HIV DNA species in resting and activated memory CD4+ T cell subsets that point to the role different levels of activation play in maintaining the HIV reservoir.

In silico single cell dynamics of hepatitis B virus infection and clearance.

UNLABELLED: The progression of acute hepatitis B virus (HBV) to chronic infection or clearance is highly dependent on the host immune response composed of cytolytic (CTL) and non-cytolytic (non-CTL) effects. Cytolytic processes induce hepatocyte killing while non-CTL processes inhibit intracellular replication. Both effects are widely recognized and accepted. However, there are uncertainties about the assistance provided by either the loss of covalently circular closed DNA (cccDNA) during cell proliferation or the emergence of refractory cells to immune mediated clearance. We developed an agent-based mathematical model and tested the relative roles of different mechanisms of the immune system in the clearance of acute HBV infection. HBV viremia clearance time and hepatocyte turnover (HT) were used as the two major criteria in determining reasonable outcomes. Modelling results in 90% of cells containing between 1 and 17 cccDNA copies and normally distributed at the peak of infection. Variations in p36 levels, responsible for determining export of virions or recirculation to amplify cccDNA numbers, have a much greater impact on mean cccDNA level/cell at peak viremia than virus infectivity and cccDNA half-life. A strong CTL effect alone failed to clear infection with HT ≈ 10. Acute infection clearance was possible with combined CTL and non-CTL effects along with complete loss of intracellular viral components during cell proliferation resulting in the desired range of HT (0.7-1). The emergence of cells refractory to infection can reduce HT by up to 90%. However their impact was less effective than complete loss of intracellular viral components during cell proliferation. CONCLUSION: the existence of refractory cells is not necessary when there is complete loss of intracellular quantities during cell proliferation but is essential with only partial clearance.

A quantitative comparison of anti-HIV gene therapy delivered to hematopoietic stem cells versus CD4+ T cells.

Gene therapy represents an alternative and promising anti-HIV modality to highly active antiretroviral therapy. It involves the introduction of a protective gene into a cell, thereby conferring protection against HIV. While clinical trials to date have delivered gene therapy to CD4+T cells or to CD34+ hematopoietic stem cells (HSC), the relative benefits of each of these two cellular targets have not been conclusively determined. In the present analysis, we investigated the relative merits of delivering a dual construct (CCR5 entry inhibitor + C46 fusion inhibitor) to either CD4+T cells or to CD34+ HSC. Using mathematical modelling, we determined the impact of each scenario in terms of total CD4+T cell counts over a 10 year period, and also in terms of inhibition of CCR5 and CXCR4 tropic virus. Our modelling determined that therapy delivery to CD34+ HSC generally resulted in better outcomes than delivery to CD4+T cells. An early one-off therapy delivery to CD34+ HSC, assuming that 20% of CD34+ HSC in the bone marrow were gene-modified (G+), resulted in total CD4+T cell counts ≥ 180 cells/ µL in peripheral blood after 10 years. If the uninfected G+ CD4+T cells (in addition to exhibiting lower likelihood of becoming productively infected) also exhibited reduced levels of bystander apoptosis (92.5% reduction) over non gene-modified (G-) CD4+T cells, then total CD4+T cell counts of ≥ 350 cells/ µL were observed after 10 years, even if initially only 10% of CD34+ HSC in the bone marrow received the protective gene. Taken together our results indicate that: 1.) therapy delivery to CD34+ HSC will result in better outcomes than delivery to CD4+T cells, and 2.) a greater impact of gene therapy will be observed if G+ CD4+T cells exhibit reduced levels of bystander apoptosis over G- CD4+T cells.

Explaining the determinants of first phase HIV decay dynamics through the effects of stage-dependent drug action.

A recent investigation of the effect of different antiretroviral drug classes on first phase dynamics of HIV RNA plasma virus levels has indicated that drugs acting at stages closer to viral production, such as the integrase inhibitor raltegravir, can produce a steeper first phase decay slope that may not be due to drug efficacy. Moreover it was found that for most drug classes the first phase transitions from a faster (phase IA) to a slightly slower decay region (phase IB) before the start of the usual second phase. Neither of these effects has been explained to date. We use a mathematical model that incorporates the different stages of the HIV viral life cycle in CD4+ T cells: viral entry, reverse transcription, integration, and viral production, to investigate the intracellular HIV mechanisms responsible for these complex plasma virus decay dynamics. We find differences in the phase IA slope across drug classes arise from a higher death rate of cells when they enter the productively infected stage post-integration, with a half-life of approximately 8 hours in this stage, whereas cells in earlier stages of the infection cycle have half-lives similar to uninfected cells. This implies any immune clearance is predominantly limited to the productive infection stage. We also show that the slowing of phase IA to phase IB at day 2 to 4 of monotherapy, depending on drug class, is a result of new rounds of infection. The level at which this slowing occurs is a better indicator of drug efficacy than the slope of the initial decay.

Hepatitis C virus envelope glycoprotein signatures are associated with treatment failure and modulation of viral entry and neutralization.

BACKGROUND: A major challenge for antiviral treatment of hepatitis C virus (HCV) infection is viral resistance, potentially resulting from the high variability of HCV envelope glycoproteins and subsequent selection of strains with enhanced infectivity and/or immune escape. METHODS: We used a bioinformatics and functional approach to investigate whether E1/E2 envelope glycoprotein structure and function were associated with treatment failure in 92 patients infected with HCV genotype 1. RESULTS: Bioinformatics analysis identified 1 sustain virological response (R)-related residue in E1 (219T) and 2 non-SVR (NR)-related molecular signatures in E2 (431A and 642V) in HCV genotype 1a. Two of these positions also appeared in minimal networks separating NR patients from R patients. HCV pseudoparticles (HCVpp) expressing 431A and 642V resulted in a decrease in antibody-mediated neutralization by pretreatment sera. 431A/HCVpp entry into Huh7.5 cells increased with overexpression of CD81 and SR-BI. Moreover, an association of envelope glycoprotein signatures with treatment failure was confirmed in an independent cohort (Virahep-C). CONCLUSIONS: Combined in silico and functional analyses demonstrate that envelope glycoprotein signatures associated with treatment failure result in an alteration of host cell entry factor use and escape from neutralizing antibodies, suggesting that virus-host interactions during viral entry contribute to treatment failure.