Modelling hepatitis D virus RNA and HBsAg dynamics during nucleic acid polymer monotherapy suggest rapid turnover of HBsAg

Hepatitis D virus (HDV) requires hepatitis B surface antigen (HBsAg) for its assembly and release. Current HBV treatments are only marginally effective against HDV because they fail to inhibit HBsAg production/secretion. However, monotherapy with the nucleic acid polymer REP 2139-Ca is accompanied by rapid declines in both HBsAg and HDV RNA. We used mathematical modeling to estimate HDV-HBsAg-host parameters and to elucidate the mode of action and efficacy of REP 2139-Ca against HDV in 12 treatment-naive HBV/HDV co-infected patients. The model accurately reproduced the observed decline of HBsAg and HDV, which was simultaneous. Median serum HBsAg half-life (t1/2) was estimated as 1.3 [0.9–1.8] days corresponding to a pretreatment production and clearance of ~108 [107.7–108.3] IU/day. The HDV-infected cell loss was estimated to be 0.052 [0.035–0.074] days−1 corresponding to an infected cell t1/2 = 13.3 days. The efficacy of blocking HBsAg and HDV production were 98.2 [94.5–99.9]% and 99.7 [96.0–99.8]%, respectively. In conclusion, both HBsAg production and HDV replication are effectively inhibited by REP 2139-Ca. Modeling HBsAg kinetics during REP 2139-Ca monotherapy indicates a short HBsAg half-life (1.3 days) suggesting a rapid turnover of HBsAg in HBV/HDV co-infection.


Results
Viral-host kinetics during REP 2139-Ca monotherapy. The kinetics of HBV DNA, anti-HBs, ALT, HBsAg and HDV for individual participants are described in Fig. 1. The pre-treatment HBV DNA was <3 log 10 IU/mL in all 12 participants and it was <LLoQ in 5 cases (Table 1). HBV DNA levels varied by <1 log 10 over the course of treatment and did not correlate with observed responses in HBsAg and HDV RNA levels.
Pre-treatment ALT had a median value of 97 U/L [interquartile range:85. . Half of the cases experienced transient ALT increases and the highest ALT flare was 281 U/L. The median pre-treatment HBsAg titer was 4.15 [3.96-4.31] log 10 IU/mL. Two patients (Pt 20 and Pt 22) experienced HBsAg decline <1 log 10 IU/mL from baseline and were classified as non-responders. A third patient, Pt 9, was a borderline responder, who had a minimal HBsAg decrease and HDV decreased only towards the end of treatment. Nonetheless, due to the magnitude of HDV decline, we chose to include Pt 9 in our modeling analysis. After a delay of 3-7 weeks during which HBsAg remained at pre-treatment levels, all responding patients exhibited a biphasic decrease in HBsAg, except for Pt 9 who had a monophasic HBsAg decline (Fig. 1). The median anti-HBs titers at week 15 post initiation of monotherapy in the 6 patients who experienced seroconversion ( Fig. 1) was 25 [IQR: [23][24][25][26][27][28][29][30]] mIU/mL.
The median pre-treatment HDV level was 6.7 [6.1-7.1] log 10 U/mL. In the 2 HBsAg non-responders, HDV decline was either 0.83 log 10 U/mL from baseline (Pt20) or transiently increased but eventually declined to 2.48 (Pt22) log 10 U/mL from baseline. After a 3-7 week delay during which HDV remained at pre-treatment levels, 5 of 10 responding patients experienced a rapid monophasic decline in HDV either reaching LLoQ or TND and 5 had a biphasic pattern with a slower 2 nd phase HDV decline (Fig. 1). By the end of therapy, the 10 responding patients experienced median HDV declines of 4.4 [3.5-5.8] log 10 U/mL. HDV was undetectable in 4 cases and below LLoQ in 3 additional patients at week 15 ( Fig. 1).
Following monotherapy, patients underwent combination therapy with REP 2139-Ca and pegIFN. A total of 9 patients were HDV RNA negative at the end of the full treatment course and 7 of these 9 patients remained HDV RNA negative 2.5 years later 8,12 .
Overall, the kinetic analysis suggests that: (i) there was no association between changes in HBV DNA and HBsAg levels or between HBV DNA and HDV RNA levels, (ii) anti-HBs sero-conversion occurred in half of responding patients and took place several weeks after HBsAg and HDV started to decline from pre-treatment www.nature.com/scientificreports www.nature.com/scientificreports/ levels and typically occurred towards the end of treatment, (iii) reminiscent of our previous findings in mono-infected HBV patients 10 , ALT elevations (with the exception of patient 1) did not correlate with changes in HBsAg or HDV levels, and (iv) after apparently similar delays HBsAg and HDV kinetics were highly correlated. Thus, we chose to model only the kinetics of HDV RNA and HBsAg in Eq. 1 (Fig. 2).
Modeling results. The model ( Fig. 2 and Eq. 1) reproduces well the HDV RNA and HBsAg kinetics in the 10 responding patients (Fig. 3) and provides estimates of unknown model parameters ( Table 2). Modeling calibration with measured data estimates median baseline HDV RNA V 0 of 6. The median estimated loss rate of infected cells δ, was 0.052 [0.035-0.074] cells/day. We note that model fits for Pt 9, who was a border-line responder were included in our analysis of the median and IQR, however as these are nonparametric measures, removing this Pt 9 from the analysis had a minimal effect on these values.
We did not find any association between the individual fit parameter values t b , ε H , ε V , and δ and baseline characteristics including duration of infection, ALT, gender, liver stiffness, V 0 , and H 0 . We found that p > 0.1 for all tested associations. In particular, the lack of association with ALT further justifies excluding ALT dynamics from our model.

Discussion
In the current study we estimated viral and host kinetic parameters under REP 2139-Ca monotherapy using a model modified from analysis of HDV replication under pegIFN monotherapy 13 . The pre-treatment median and interquartile range for HDV RNA (V 0 = 6.6 [6.1-7.0] U/mL) were similar to previously studied pegIFN-treated patients (V 0 = 7. The reason for this difference in loss rates is unclear, particularly as REP 2139-Ca has not previously been shown to increase loss of infected cells. This difference could be the result of an immune response, although, thus far REP 2139 has not been shown to augment the host immune response to HDV 8,10,11 . Unfortunately, the effects of NAPs observed in humans were not reproduced in rodent models 14 , so alternative methods are needed to determine whether NAPs do affect the immune system and then models could be developed to incorporate such a response. The estimation of HBsAg turnover based on pegIFN inhibition kinetics is confounded by the multiple antiviral mechanisms of pegIFN, which not only affects viral replication but alters immune function 15,16 . Additionally, we previously showed 13 that under pegIFN it is not feasible to estimate HBsAg turnover because decreases in HBsAg occurred only during the second phase decrease of HDV RNA, which corresponds to cell death/loss. In contrast, REP 2139-Ca targets HBsAg SVP assembly and secretion 17 , which is the source of almost all circulating HBsAg, allowing for a direct assessment of HBsAg turnover. The estimated half-life of HBsAg (1.3 days) under REP 2139-Ca is strikingly short compared to the half-life of 38 days estimated under lamivudine 18 and approximately 7-fold shorter than estimated with deuterated HBsAg 19 suggesting that the turnover of HBV SVP may be more rapid in these patients than estimated in previous studies. However, both studies only examined a Figure 2. A schematic description of the model (Eq.1). Target cells, T 0 , are infected with rate β, and become infectious cells I. Infectious cells loss at a rate , and produce virions, V, at rate p. After treatment the production rate of virions is reduced by a factor (1-ε v ). Virions are cleared at rate c. Infectious cells also produce HBsAg, H, at rate P H . This rate is reduced by a factor (1-ε H ) after treatment. HBsAg is cleared at rate c H . As was done previously 13 , we assume that T 0 was constant during the 15 weeks of treatment at its pre-treatment steady-state value. (2020) 10:7837 | https://doi.org/10.1038/s41598-020-64122-0 www.nature.com/scientificreports www.nature.com/scientificreports/ small subset of patients and the estimates are within an order of magnitude. The decline in circulating HBsAg is likely driven at least in part by a reduction in HBsAg production/secretion 20,21 . Nonetheless it is also possible that NAPs may increase HBsAg clearance through some yet unknown mechanism. The previous studies investigated HBV mono-infections, whereas the current patients were HBV/HDV coinfected, raising the possibility that co-infections behave differently, and this could perhaps lead to the observed differences.
We also note that our model does not include intracellular dynamics. The identified mechanism of action of REP 2139-Ca consisting of blocking SVP assembly and production with an expected corresponding blockage of HDV production (which is dependent on SVP morphogenesis) is consistent with the patient data. However, the data do not eliminate the possibility of additional modes of action. Preliminary analysis of HDV protein interaction suggests that REP 2139 binds to the small and large forms of HDAg in vitro 22 similar to other sequence independent oligonucleotide interactions with HDAg 23 . These interactions could potentially result in a direct antiviral effect against HDV. However, more frequent sampling in additional patients is needed to separate the onset of blocking HBsAg and HDV production and investigate this possible mechanism using a multiscale modeling approach. Further studies are needed to validate the model presented here and to determine how different concentrations of NAPs relate to the mode of action and effectivity of treatment. Such studies also should investigate changes in HBsAg isoforms during NAP therapy and the effects of NAPs on the release of HBV SVP, Dane particles / HBV filaments and HDV virions.
In conclusion, by simultaneously modeling HBsAg and HDV kinetic data under REP 2139-Ca monotherapy, we estimated for the first time the efficacy of REP 2139-Ca in blocking HBsAg and HDV production, as well as the serum HBsAg half-life and the loss/death rate of HDV-infected cells. These analyses demonstrate a potent effect of REP 2139-Ca against HBsAg and HDV and also seem to indicate that the turnover of HBsAg SVP in the serum may be faster than previously estimated.

Methods
Patients. The study population consisted of 12 treatment-naive, HBeAg-negative, HDV RNA positive participants with serum HBsAg titers >1000 IU/mL who were treated with REP 2139-Ca (the calcium chelate complex formulation of REP 2139) for 15 weeks in the phase IIA REP 301 clinical trial 8 . Baseline participant characteristics were as previously described 8 (Table 1). All methods were carried out in accordance with the Helsinki declaration HDV RNA, HBV DNA, anti-HBs, HBsAg, and ALT measurements. Serum HDV RNA (Robogene MK I), HBV DNA (Abbott Realtime) and HBsAg levels (Abbott Architect quantitative) were measured every two weeks. For HDV, the lower limit of quantification (LLoQ) was 3.26 log U/mL, for HBsAg, LLoQ = −1.3 log IU/mL and for HBV DNA, LLoQ=1 log IU/mL. For anti-HBs the architect LLoQ is 2mIU/mL and the current convention, which we followed, is to define seroconversion as a titer> 10mIU/mL 24 . The upper limit of normal for ALT was 50 U/L.

Modeling.
Previous in vitro studies demonstrated that REP 2139 blocks the assembly of HBV subviral particles (SVPs), simultaneously lowering intracellular HBsAg and blocking HBsAg secretion from SVPs 17 . This effect is driven by a post entry mechanism as REP 2139 does not block entry of HBV or HDV into the host cell 25 . Therefore, we modified our previous dual-mathematical model of HDV RNA and HBsAg dynamics under pegIFN therapy 13 by adding a possible effect of REP 2139 in blocking both HDV and HBsAg production (Eq. 1 and Fig. 2). The modified model includes the following parameters: I = productively HDV infected cells,   www.nature.com/scientificreports www.nature.com/scientificreports/ V = HDV RNA, H = HBsAg, and T 0 = the number of target/susceptible cells (i.e., HBsAg-producing cells). V, infects T 0 with constant rate β, generating infected cells, I. Parameter δ is the loss rate of HDV-infected cells, p is the production rate of virions, c is the clearance rate constant of virions, p H is the production rate of HBsAg, and c H is the clearance rate constant of HBsAg. Treatment is assumed to begin blocking HDV and HBsAg production after time t b , with efficacy ε V and ε H , respectively. Model equations are: HBsAg production rate. Having stable levels of HBsAg implies production and clearance are in balance before treatment. Therefore, from Eq. 1 the production rate of serum HBsAg before treatment initiation must equal the HBsAg clearance rate c H H 0 . Assuming that the total body fluid volume, F, was 13,360 mL for body weight of 70 kg, as in our previous study 13 we estimated the total HBsAg production in each patient before treatment initiation by the product F·c H ·H 0 .

HBsAg production ** (IU/day) V 0 (log(U) /mL)
Statistical analysis. Non-parametric Spearman and Mann-Whitney U Tests were performed using Python version 3.7 and Scipy version 1.3. In all cases P ≤ 0.05 was considered significant.