Project description:Background & aimsBulevirtide (BLV) was approved for the treatment of compensated chronic hepatitis D virus (HDV) infection in Europe in 2020. However, research into the effects of the entry inhibitor BLV on HDV-host dynamics is in its infancy.MethodsEighteen patients with HDV under nucleos(t)ide analogue treatment for hepatitis B, with compensated cirrhosis and clinically significant portal hypertension, received BLV 2 mg/day. HDV RNA, alanine aminotransferase (ALT), and hepatitis B surface antigen (HBsAg) were measured at baseline, weeks 4, 8 and every 8 weeks thereafter. A mathematical model was developed to account for HDV, HBsAg and ALT dynamics during BLV treatment.ResultsMedian baseline HDV RNA, HBsAg, and ALT were 4.9 log IU/ml [IQR: 4.4-5.8], 3.7 log IU/ml [IQR: 3.4-3.9] and 106 U/L [IQR: 81-142], respectively. During therapy, patients fit into four main HDV kinetic patterns: monophasic (n = 2), biphasic (n = 10), flat-partial response (n = 4), and non-responder (n = 2). ALT normalization was achieved in 14 (78%) patients at a median of 8 weeks (range: 4-16). HBsAg remained at pre-treatment levels. Assuming that BLV completely (∼100%) blocks HDV entry, modeling indicated that two HDV-infected cell populations exist: fast HDV clearing (median t1/2 = 13 days) and slow HDV clearing (median t1/2 = 44 days), where the slow HDV-clearing population consisted of ∼1% of total HDV-infected cells, which could explain why most patients exhibited a non-monophasic pattern of HDV decline. Moreover, modeling explained ALT normalization without a change in HBsAg based on a non-cytolytic loss of HDV from infected cells, resulting in HDV-free HBsAg-producing cells that release ALT upon death at a substantially lower rate compared to HDV-infected cells.ConclusionThe entry inhibitor BLV provides a unique opportunity to understand HDV, HBsAg, ALT, and host dynamics.Impact and implicationsMathematical modeling of hepatitis D virus (HDV) treatment with the entry inhibitor bulevirtide (BLV) provides a novel window into the dynamics of HDV RNA and alanine aminotransferase. Kinetic data from patients treated with BLV monotherapy can be explained by hepatocyte populations with different basal HDV clearance rates and non-cytolytic clearance of infected cells. While further studies are needed to test and refine the kinetic characterization described here, this study provides a new perspective on viral dynamics, which could inform evolving treatment strategies for HDV.
Project description:Viruses are major pathogens of agricultural crops. Viral infections often start after the virus enters the outer layer of a tissue, and many successful viruses, after local replication in the infected tissue, are able to spread systemically. Quantitative details of virus dynamics in plants, however, are poorly understood, in part, because of the lack of experimental methods which allow the accurate measurement of the degree of infection in individual plant tissues. Recently, a group of researchers followed the kinetics of infection of individual cells in leaves of Nicotiana tabacum plants using Tobacco etch virus (TEV) expressing either Venus or blue fluorescent protein (BFP). Assuming that viral spread occurs from lower to upper leaves, the authors fitted a simple mathematical model to the frequency of cellular infection by the two viral variants found using flow cytometry. While the original model could accurately describe the kinetics of viral spread locally and systemically, we found that many alternative versions of the model, for example, if viral spread starts at upper leaves and progresses to lower leaves or when virus dissemination is stopped due to an immune response, fit the data with reasonable quality, and yet with different parameter estimates. These results strongly suggest that experimental measurements of the virus infection in individual leaves may not be sufficient to identify the pathways of viral dissemination between different leaves and reasons for viral control. We propose experiments that may allow discrimination between the alternatives. By analyzing the kinetics of coinfection of individual cells by Venus and BFP strains of TEV we found a strong deviation from the random infection model, suggesting cooperation between the two strains when infecting plant cells. Importantly, we showed that many mathematical models on the kinetics of coinfection of cells with two strains could not adequately describe the data, and the best fit model needed to assume (i) different susceptibility of uninfected cells to infection by two viruses locally in the leaf vs. systemically from other leaves, and (ii) decrease in the infection rate depending on the fraction of uninfected cells which could be due to a systemic immune response. Our results thus demonstrate the difficulty in reaching definite conclusions from extensive and yet limited experimental data and provide evidence of potential cooperation between different viral variants infecting individual cells in plants.
Project description:Introduction: Bulevirtide is a first-in-class antiviral drug to treat chronic hepatitis B/D. We investigated the drug-drug interaction potential and pharmacokinetics of high-dose subcutaneous bulevirtide (5 mg twice daily) with organic anion transporting polypeptide 1B1 (OATP1B1) and cytochrome P450 (CYP) 3A4. Methods: This was a single-center, open-label, fixed-sequence drug-drug interaction trial in 19 healthy volunteers. Before and at bulevirtide steady state, participants ingested a single 40 mg dose of pravastatin. A midazolam microdose was applied to quantify CYP3A4 activity. Results: At bulevirtide steady state, pravastatin area under the concentration-time curve (AUC0-∞) increased 1.32-fold (90% CI 1.08-1.61). The 5 mg bulevirtide twice-daily treatment resulted in a mean AUC0-12 of 1210 h*ng/ml (95% CI 1040-1408) and remained essentially unchanged under the influence of pravastatin. CYP3A4 activity did not change to a clinically relevant extent. As expected, total bile acids increased substantially (35-fold) compared to baseline during bulevirtide treatment. All study medication was well tolerated. Discussion: The study demonstrated that high-dose bulevirtide inhibited OATP1B-mediated hepatic uptake of the marker substrate pravastatin but the extent is considered clinically not relevant. Changes in CYP3A4 activity were also not clinically relevant. In conclusion, this study suggests that OATP1B substrate drugs as well as CYP3A4 substrates may safely be used without dose adjustment in patients treated with bulevirtide. However, in patients using high statin doses and where concomitant factors potentially further increase statin exposure, caution may be required when using bulevirtide.
Project description:BackgroundAnalyzing the dynamics of insulin concentration in the blood is necessary for a comprehensive understanding of the effects of insulin in vivo. Insulin removal from the blood has been addressed in many studies. The results are highly variable with respect to insulin clearance and the relative contributions of hepatic and renal insulin degradation.ResultsWe present a dynamic mathematical model of insulin concentration in the blood and of insulin receptor activation in hepatocytes. The model describes renal and hepatic insulin degradation, pancreatic insulin secretion and nonspecific insulin binding in the liver. Hepatic insulin receptor activation by insulin binding, receptor internalization and autophosphorylation is explicitly included in the model. We present a detailed mathematical analysis of insulin degradation and insulin clearance. Stationary model analysis shows that degradation rates, relative contributions of the different tissues to total insulin degradation and insulin clearance highly depend on the insulin concentration.ConclusionThis study provides a detailed dynamic model of insulin concentration in the blood and of insulin receptor activation in hepatocytes. Experimental data sets from literature are used for the model validation. We show that essential dynamic and stationary characteristics of insulin degradation are nonlinear and depend on the actual insulin concentration.
Project description:Hepatitis C virus (HCV) entry inhibitors (EIs) act synergistically with drugs targeting other stages of the HCV lifecycle. The origin of this synergy remains unknown. Here, we argue that the synergy may arise from the complementary activities of the drugs across cell subpopulations expressing different levels of HCV entry receptors. We employ mathematical modeling of viral kinetics in vitro, where cells with a distribution of entry receptor expression levels are exposed to HCV with or without drugs. The drugs act independently in each cell, as expected in the absence of underlying interactions. Yet, at the cell population level our model predicts that the drugs exhibit synergy. EIs effectively block infection of cells with low receptor levels. With high receptor levels, where EIs are compromised, other drugs are potent. This novel mechanism of synergy, arising at the cell population level may facilitate interpretation of drug activity and treatment optimization.
Project description:Tuberculosis (TB), the disease caused by Mycobacterium tuberculosis (Mtb), remains a major health problem with 10.6 million cases of the disease and 1.6 million deaths in 2021. It is well understood that pulmonary TB is due to Mtb growth in the lung but quantitative estimates of rates of Mtb replication and death in lungs of patients or animals such as monkeys or rabbits remain largely unknown. We performed experiments with rabbits infected with a novel, virulent clinical Mtb isolate of the Beijing lineage, HN878, carrying an unstable plasmid pBP10. In our in vitro experiments we found that pBP10 is more stable in HN878 strain than in a more commonly used laboratory-adapted Mtb strain H37Rv (the segregation coefficient being s = 0.10 in HN878 vs. s = 0.18 in H37Rv). Interestingly, the kinetics of plasmid-bearing bacteria in lungs of Mtb-infected rabbits did not follow an expected monotonic decline; the percent of plasmid-bearing cells increased between 28 and 56 days post-infection and remained stable between 84 and 112 days post-infection despite a large increase in bacterial numbers in the lung at late time points. Mathematical modeling suggested that such a non-monotonic change in the percent of plasmid-bearing cells can be explained if the lung Mtb population consists of several (at least 2) sub-populations with different replication/death kinetics: one major population expanding early and being controlled/eliminated, while another, a smaller population expanding at later times causing a counterintuitive increase in the percent of plasmid-bearing cells. Importantly, a model with one kinetically homogeneous Mtb population could not explain the data including when the model was run stochastically. Given that in rabbits HN878 strain forms well circumscribed granulomas, our results suggest independent bacterial dynamics in subsets of such granulomas. Our model predictions can be tested in future experiments in which HN878-pBP10 dynamics in individual granulomas is followed over time. Taken together, our new data and mathematical modeling-based analyses illustrate differences in Mtb dynamics in mice and rabbits confirming a perhaps somewhat obvious observation that "rabbits are not mice".
Project description:The mechanism by which hepatitis C virus (HCV) gains entry into cells is a complex one, involving a broad range of host proteins. Entry is a critical phase of the viral lifecycle, and a potential target for therapeutic or vaccine-mediated intervention. However, the mechanics of HCV entry remain poorly understood. Here we describe a novel computational model of viral entry, encompassing the relationship between HCV and the key host receptors CD81 and SR-B1. We conduct experiments to thoroughly quantify the influence of an increase or decrease in receptor availability upon the extent of viral entry. We use these data to build and parameterise a mathematical model, which we then validate by further experiments. Our results are consistent with sequential HCV-receptor interactions, whereby initial interaction between the HCV E2 glycoprotein and SR-B1 facilitates the accumulation CD81 receptors, leading to viral entry. However, we also demonstrate that a small minority of viruses can achieve entry in the absence of SR-B1. Our model estimates the impact of the different obstacles that viruses must surmount to achieve entry; among virus particles attaching to the cell surface, around one third of viruses accumulate sufficient CD81 receptors, of which 4-8% then complete the subsequent steps to achieve productive infection. Furthermore, we make estimates of receptor stoichiometry; in excess of 10 receptors are likely to be required to achieve viral entry. Our model provides a tool to investigate the entry characteristics of HCV variants and outlines a framework for future quantitative studies of the multi-receptor dynamics of HCV entry.
Project description:Trypanosome RNA polymerase II transcription is polycistronic, individual mRNAs being excised by trans splicing and polyadenylation. In this study, we refined the previously published mathematical model for bloodstream form parasites and extended it to the procyclic form. We used the model, together with known mRNA half-lives, to predict the abundances of individual mRNAs, assuming rapid, unregulated mRNA processing; then we compared the results with measured mRNA abundances. Remarkably, the abundances of most mRNAs in procyclic forms are predicted quite well by the model, being largely explained by variations in mRNA decay rates and length. In bloodstream forms substantially more mRNAs are less abundant than predicted. We list mRNAs that are likely to show particularly slow or inefficient processing, either in both forms or with developmental regulation. We also measured ribosome occupancies of all mRNAs in trypanosomes grown in the same conditions as were used to measure mRNA turnover. In procyclic forms there was a weak positive correlation between ribosome density and mRNA half-life, suggesting cross-talk between translation and mRNA decay; ribosome density was related to the proportion of the mRNA on polysomes, indicating control of translation initiation. Ribosomal protein mRNAs in procyclics appeared to be exceptionally rapidly processed but poorly translated. Through this study, we conclude that lLevels of mRNAs in procyclic form trypanosomes are determined mainly by length and mRNA decay, with some control of precursor processing. In bloodstream forms variations in nuclear events play a larger role in transcriptome regulation, suggesting acquisition of new control mechanisms during adaptation to mammalian parasitism. Ribosome profiling and mRNA libraries were constructed in triplicate from in vitro PCF and in duplicate from in vitro T. brucei Lister427, to understand global differntial gene transcription.