Modeling Favipiravir Antiviral Efficacy Against Emerging Viruses: From Animal Studies to Clinical Trials

In 2014, our research network was involved in the evaluation of favipiravir, an anti-influenza polymerase inhibitor, against Ebola virus. In this review, we discuss how mathematical modeling was used, first to propose a relevant dosing regimen in humans, and then to optimize its antiviral efficacy in a nonhuman primate (NHP) model. The data collected in NHPs were finally used to develop a model of Ebola pathogenesis integrating the interactions among the virus, the innate and adaptive immune response, and the action of favipiravir. We conclude the review of this work by discussing how these results are of relevance for future human studies in the context of Ebola virus, but also for other emerging viral diseases for which no therapeutics are available.

Antiviral efficacy of favipiravir against Ebola virus: A translational study in cynomolgus macaques

BACKGROUND

Despite repeated outbreaks, in particular the devastating 2014-2016 epidemic, there is no effective treatment validated for patients with Ebola virus disease (EVD). Among the drug candidates is the broad-spectrum polymerase inhibitor favipiravir, which showed a good tolerance profile in patients with EVD (JIKI trial) but did not demonstrate a strong antiviral efficacy. In order to gain new insights into the antiviral efficacy of favipiravir and improve preparedness and public health management of future outbreaks, we assess the efficacy achieved by ascending doses of favipiravir in Ebola-virus-infected nonhuman primates (NHPs).

METHODS AND FINDINGS

A total of 26 animals (Macaca fascicularis) were challenged intramuscularly at day 0 with 1,000 focus-forming units of Ebola virus Gabon 2001 strain and followed for 21 days (study termination). This included 13 animals left untreated and 13 treated with doses of 100, 150, and 180 mg/kg (N = 3, 5, and 5, respectively) favipiravir administered intravenously twice a day for 14 days, starting 2 days before infection. All animals left untreated or treated with 100 mg/kg died within 10 days post-infection, while animals receiving 150 and 180 mg/kg had extended survival (P < 0.001 and 0.001, respectively, compared to untreated animals), leading to a survival rate of 40% (2/5) and 60% (3/5), respectively, at day 21. Favipiravir inhibited viral replication (molecular and infectious viral loads) in a drug-concentration-dependent manner (P values < 0.001), and genomic deep sequencing analyses showed an increase in virus mutagenesis over time. These results allowed us to identify that plasma trough favipiravir concentrations greater than 70-80 μg/ml were associated with reduced viral loads, lower virus infectivity, and extended survival. These levels are higher than those found in the JIKI trial, where patients had median trough drug concentrations equal to 46 and 26 μg/ml at day 2 and day 4 post-treatment, respectively, and suggest that the dosing regimen in the JIKI trial was suboptimal. The environment of a biosafety level 4 laboratory introduces a number of limitations, in particular the difficulty of conducting blind studies and performing detailed pharmacological assessments. Further, the extrapolation of the results to patients with EVD is limited by the fact that the model is fully lethal and that treatment initiation in patients with EVD is most often initiated several days after infection, when symptoms and high levels of viral replication are already present.

CONCLUSIONS

Our results suggest that favipiravir may be an effective antiviral drug against Ebola virus that relies on RNA chain termination and possibly error catastrophe. These results, together with previous data collected on tolerance and pharmacokinetics in both NHPs and humans, support a potential role for high doses of favipiravir for future human interventions.

Viral kinetics analysis and virological characterization of treatment failures in patients with chronic hepatitis C treated with sofosbuvir and an NS5A inhibitor

BACKGROUND

The combination of sofosbuvir (SOF) plus an NS5A inhibitor for 12 weeks is highly efficacious in patients with chronic hepatitis C. As the costs of generic production of sofosbuvir and NS5A inhibitor are rapidly decreasing, the combination of these DAAs will be the standard treatment in most low- to middle-income countries in the future.

AIM

To identify key predictors of response that can be used to tailor treatment decisions.

METHODS

A cohort of 216 consecutive patients infected with HCV genotype 1 (1a: n = 57; 1b: n = 77), 2 (n = 4), 3 (n = 33) or 4 (n = 44) were treated with sofosbuvir (SOF) + daclatasvir (n = 176) or SOF + ledipasvir (n = 40) for 12 weeks. The viral kinetics was analysed using the biphasic model and the cure boundary was used to predict time to clear HCV.

RESULTS

The overall SVR rate was high (94.4%; n = 204), regardless of the time to viral suppression or low-level viraemia at the end of treatment. The model-based predicted HCV RNA levels at the end of treatment could not differentiate patients who did from those who did not achieve SVR. The presence of NS5A resistance-associated substitutions [position 28 (OR = 70.3, P<.001) and/or 31 (OR = 61.6, P = .002)] at baseline was predictive of virological failure in cirrhotic patients but was not associated with on-treatment viral kinetics.

CONCLUSION

This real-world study confirms the excellent results of clinical trials with therapies based on a combination of SOF plus an NS5A inhibitor. It suggests that a personalized approach including baseline NS5A inhibitor resistance testing may inform treatment decisions in cirrhotic patients.

Persistence and clearance of Ebola virus RNA from seminal fluid of Ebola virus disease survivors: a longitudinal analysis and modelling study

BACKGROUND

By January, 2016, all known transmission chains of the Ebola virus disease (EVD) outbreak in west Africa had been stopped. However, there is concern about persistence of Ebola virus in the reproductive tract of men who have survived EVD. We aimed to use biostatistical modelling to describe the dynamics of Ebola virus RNA load in seminal fluid, including clearance parameters.

METHODS

In this longitudinal study, we recruited men who had been discharged from three Ebola treatment units in Guinea between January and July, 2015. Participants provided samples of seminal fluid at follow-up every 3-6 weeks, which we tested for Ebola virus RNA using quantitative real-time RT-PCR. Representative specimens from eight participants were then inoculated into immunodeficient mice to test for infectivity. We used a linear mixed-effect model to analyse the dynamics of virus persistence in seminal fluid over time.

FINDINGS

We enrolled 26 participants and tested 130 seminal fluid specimens; median follow up was 197 days (IQR 187-209 days) after enrolment, which corresponded to 255 days (228-287) after disease onset. Ebola virus RNA was detected in 86 semen specimens from 19 (73%) participants. Median duration of Ebola virus RNA detection was 158 days after onset (73-181; maximum 407 days at end of follow-up). Mathematical modelling of the quantitative time-series data showed a mean clearance rate of Ebola virus RNA from seminal fluid of -0·58 log units per month, although the clearance kinetic varied greatly between participants. Using our biostatistical model, we predict that 50% and 90% of male survivors clear Ebola virus RNA from seminal fluid at 115 days (90% prediction interval 72-160) and 294 days (212-399) after disease onset, respectively. We also predicted that the number of men positive for Ebola virus RNA in affected countries would decrease from about 50 in January 2016, to fewer than 1 person by July, 2016. Infectious virus was detected in 15 of 26 (58%) specimens tested in mice.

INTERPRETATION

Time to clearance of Ebola virus RNA from seminal fluid varies greatly between individuals and could be more than 13 months. Our predictions will assist in decision-making about surveillance and preventive measures in EVD outbreaks.

FUNDING

This study was funded by European Union's Horizon 2020 research and innovation programme, Directorate-General for International Cooperation and Development of the European Commission, Institut national de la santé et de la recherche médicale (INSERM), German Research Foundation (DFG), and Innovative Medicines Initiative 2 Joint Undertaking.

Determining Ribavirin's mechanism of action against Lassa virus infection

Ribavirin is a broad spectrum antiviral which inhibits Lassa virus (LASV) replication in vitro but exhibits a minor effect on viremia in vivo. However, ribavirin significantly improves the disease outcome when administered in combination with sub-optimal doses of favipiravir, a strong antiviral drug. The mechanisms explaining these conflicting findings have not been determined, so far. Here, we used an interdisciplinary approach combining mathematical models and experimental data in LASV-infected mice that were treated with ribavirin alone or in combination with the drug favipiravir to explore different putative mechanisms of action for ribavirin. We test four different hypotheses that have been previously suggested for ribavirin’s mode of action: (i) acting as a mutagen, thereby limiting the infectivity of new virions; (ii) reducing viremia by impairing viral production; (iii) modulating cell damage, i.e., by reducing inflammation, and (iv) enhancing antiviral immunity. Our analysis indicates that enhancement of antiviral immunity, as well as effects on viral production or transmission are unlikely to be ribavirin’s main mechanism mediating its antiviral effectiveness against LASV infection. Instead, the modeled viral kinetics suggest that the main mode of action of ribavirin is to protect infected cells from dying, possibly reducing the inflammatory response.

The paradox of highly effective sofosbuvir-based combination therapy despite slow viral decline: can we still rely on viral kinetics?

High sustained virologic response (SVR) rates have been observed after 6 weeks of anti-HCV treatment using sofosbuvir, ledipasvir and a non-nucleoside polymerase-inhibitor (GS-9669) or a protease-inhibitor (GS-9451) and after 12 weeks with sofosbuvir + ledipasvir. Here we analyze the viral kinetics observed during these treatments to decipher the origin of the rapid cure and to evaluate the possibility of further reducing treatment duration. We found that viral kinetics were surprisingly slow in all treatment groups and could not reproduce the high SVR rates observed. Based on experimental results suggesting that NS5A- or protease-inhibitors can generate non-infectious virus, we incorporated this effect into a mathematical model. We found that to predict observed SVR rates it was necessary to assume that ledipasvir, GS-9669 and GS-9451 rapidly reduce virus infectivity. We predicted with this model that 4 weeks of triple therapy could be sufficient to achieve SVR in patients with undetectable viremia at week 1, but would be suboptimal in general. In conclusion, the rapid cure rate achieved with these combinations is largely disconnected from viral loads measured during treatment. A model assuming that rapid cure is due to a drug effect of generating non-infectious virus could be a basis for future response guided therapy.

Improvement of ALT decay kinetics by all-oral HCV treatment: Role of NS5A inhibitors and differences with IFN-based regimens

BACKGROUND

Intracellular HCV-RNA reduction is a proposed mechanism of action of direct-acting antivirals (DAAs), alternative to hepatocytes elimination by pegylated-interferon plus ribavirin (PR). We modeled ALT and HCV-RNA kinetics in cirrhotic patients treated with currently-used all-DAA combinations to evaluate their mode of action and cytotoxicity compared with telaprevir (TVR)+PR.

STUDY DESIGN

Mathematical modeling of ALT and HCV-RNA kinetics was performed in 111 HCV-1 cirrhotic patients, 81 treated with all-DAA regimens and 30 with TVR+PR. Kinetic-models and Cox-analysis were used to assess determinants of ALT-decay and normalization.

RESULTS

HCV-RNA kinetics was biphasic, reflecting a mean effectiveness in blocking viral production >99.8%. The first-phase of viral-decline was faster in patients receiving NS5A-inhibitors compared to TVR+PR or sofosbuvir+simeprevir (p<0.001), reflecting higher efficacy in blocking assembly/secretion. The second-phase, noted δ and attributed to infected-cell loss, was faster in patients receiving TVR+PR or sofosbuvir+simeprevir compared to NS5A-inhibitors (0.27 vs 0.21 d-1, respectively, p = 0.0012). In contrast the rate of ALT-normalization, noted λ, was slower in patients receiving TVR+PR or sofosbuvir+simeprevir compared to NS5A-inhibitors (0.17 vs 0.27 d-1, respectively, p<0.001). There was no significant association between the second-phase of viral-decline and ALT normalization rate and, for a given level of viral reduction, ALT-normalization was more profound in patients receiving DAA, and NS5A in particular, than TVR+PR.

CONCLUSIONS

Our data support a process of HCV-clearance by all-DAA regimens potentiated by NS5A-inhibitor, and less relying upon hepatocyte death than IFN-containing regimens. This may underline a process of "cell-cure" by DAAs, leading to a fast improvement of liver homeostasis.

Model Evaluation of Continuous Data Pharmacometric Models: Metrics and Graphics

This article represents the first in a series of tutorials on model evaluation in nonlinear mixed effect models (NLMEMs), from the International Society of Pharmacometrics (ISoP) Model Evaluation Group. Numerous tools are available for evaluation of NLMEM, with a particular emphasis on visual assessment. This first basic tutorial focuses on presenting graphical evaluation tools of NLMEM for continuous data. It illustrates graphs for correct or misspecified models, discusses their pros and cons, and recalls the definition of metrics used.

Ebola Virus Infection: Review of the Pharmacokinetic and Pharmacodynamic Properties of Drugs Considered for Testing in Human Efficacy Trials

The 2014-2015 outbreak of Ebola virus disease is the largest epidemic to date in terms of the number of cases, deaths, and affected areas. In October 2015, no antiviral agents had proven antiviral efficacy in patients. However, in September 2014, the World Health Organization inventoried and has since regularly updated a list of potential drug candidates with demonstrated antiviral efficacy in in vitro or animal models. This includes agents belonging to various therapeutic classes, namely direct antiviral agents (favipiravir and BCX4430), a combination of antibodies (ZMapp), type I interferons, RNA interference-based drugs (TKM-Ebola and AVI-7537), and anticoagulant drugs (rNAPc2). Here, we review the pharmacokinetic and pharmacodynamic information presently available for these drugs, using data obtained in healthy volunteers for pharmacokinetics and data obtained in human clinical trials or animal models for pharmacodynamics. Future studies evaluating these drugs in clinical trials are critical to confirm their efficacy in humans, propose appropriate doses, and evaluate the possibility of treatment combinations.

Ebola virus dynamics in mice treated with favipiravir

The polymerase inhibitor favipiravir is a candidate for the treatment of Ebola virus disease. Here, we designed a mathematical model to characterize the viral dynamics in 20 mice experimentally infected with Ebola virus, which were either left untreated or treated with favipiravir at 6 or 8days post infection. This approach provided estimates of kinetic parameters of Ebola virus reproduction, such as the half-life of productively infected cells, of about 6h, and the basic reproductive number which indicates that virus produced by a single infected cell productively infects about 9 new cells. Furthermore, the model predicted that favipiravir efficiently blocks viral production, reaching an antiviral effectiveness of 95% and 99.6% at 2 and 6days after initiation of treatment, respectively. The model could be particularly helpful to guide future studies evaluating favipiravir in larger animals.