Original Research ARTICLE
Investigating Different Mechanisms of Action in Combination Therapy for Influenza
- 1Physics Department, East Carolina University, Greenville, NC, United States
- 2Department of Physics and Astronomy, Texas Christian University, Fort Worth, TX, United States
Combination therapy for influenza can have several benefits, from reducing the emergence of drug resistant virus strains to decreasing the cost of antivirals. However, there are currently only two classes of antivirals approved for use against influenza, limiting the possible combinations that can be considered for treatment. However, new antivirals are being developed that target different parts of the viral replication cycle, and their potential for use in combination therapy should be considered. The role of antiviral mechanism of action in the effectiveness of combination therapy has not yet been systematically investigated to determine whether certain antiviral mechanisms of action pair well in combination. Here, we use a mathematical model of influenza to model combination treatment with antivirals having different mechanisms of action to measure peak viral load, infection duration, and synergy of different drug combinations. We find that antivirals that lower the infection rate and antivirals that increase the duration of the eclipse phase perform poorly in combination with other antivirals.
Influenza, more commonly known as the flu, is a seasonal illness that has symptoms of runny nose, cough, fever and an aching body (Khandaker et al., 2011). While most people recover fairly quickly, influenza can be lethal, especially in children and the elderly (Pop-Vicas and Gravenstein, 2011; Ruf and Knuf, 2014). There are vaccines to prevent influenza, though unfortunately, they are strain-specific (Jang and Seong, 2014), change annually, and don't always match the circulating strain (Dos Santos et al., 2016). This sometimes leaves antivirals as our only defense against influenza. However, the usefulness of antivirals is limited by the fast mutation rate of influenza (Sanjuán and Domingo-Calap, 2016; Villa and Lässig, 2017) and its ability to quickly develop resistance to antivirals (Perelson et al., 2012; Dobrovolny and Beauchemin, 2017). Of the two classes of drugs currently approved for use against influenza, one (adamantanes) is rarely used anymore due to high resistance in circulating strains (Bright et al., 2006; Dong et al., 2015). The class of drugs most broadly used against influenza is the neuraminidase inhibitors (NAIs), such as oseltamivir, zanamivir, and peramivir. While most strains of influenza are currently susceptible to NAIS, resistance to oseltamivir rose rapidly to 98–100% of circulating strains in the 2008–2009 flu season (Dharan et al., 2009; Zaraket et al., 2010), showing that resistance to these antivirals can suddenly emerge. A number of polymerase complex inhibitors (pimodivir, faviparivir, and S-033188) are also in development (Shaw, 2017; Yuan et al., 2018). Since these antivirals are fairly new, little is yet known about influenza's ability to develop resistance, although many think that since these domains are highly conserved, resistance is less likely to develop (Zhang et al., 2018). As new antivirals are developed (Naesens et al., 2016), physicians will need to develop treatment strategies that will limit the emergence of drug resistance.
Combination therapy, the use of two or more antivirals, is one possible strategy for limiting drug resistance. While influenza can very quickly develop resistance to a single antiviral, it takes longer to develop resistance to two drugs simultaneously (Perelson et al., 2012). In using two drugs then, we ensure that nearly all virus produced during the infection are susceptible to at least one of the drugs, so they will not replicate and the infection will be suppressed. There have been many experimental studies, both in vivo and in vitro, investigating the effect of different combinations of influenza antivirals (Hayden et al., 1980; Smee et al., 2002, 2009, 2010a,b; Govorkova et al., 2004; Ilyushina et al., 2006, 2007, 2008; Masihi et al., 2007; Bantia et al., 2010; Duval et al., 2010; Kim et al., 2010; Chen et al., 2011; Haasbach et al., 2013; Seo et al., 2013; Tarbet et al., 2014; Belardo et al., 2015; Morokutti-Kurz et al., 2015; Marathe et al., 2016; Beigel et al., 2017; de Mello et al., 2018), with some examining not just combinations of two drugs, but even examining triple combinations (Nguyen et al., 2009, 2010, 2012; Hoopes et al., 2011; Kim et al., 2011; Lu et al., 2015; Pavlova et al., 2018). These experimental studies are largely limited to combinations of antivirals from the two widely available classes of influenza antivirals, neuraminidase inhibitors (NAIs) and adamantanes.
The mechanisms of action of NAIs and adamantanes are at different points in the viral replication cycle with adamantanes acting to block viral replication (Abed et al., 2005) and NAIs acting to block viral release (Gubareva et al., 2000; Abed et al., 2002). However, new antivirals with different mechanisms of action are being developed (Koszalka et al., 2017; White et al., 2018; Zabrodskaya et al., 2018), and will potentially be used in combination therapy (Loregian et al., 2014; Koszalka et al., 2017). With the development of these new antivirals, there has been some interest in exploring how the mechanism of action of antivirals involved in combination therapy affects the effectiveness of combination therapy (Dunning et al., 2014; Popov et al., 2018). This is most often measured by assessing, via in vitro experiments, the synergy or antagony of the drug combination. Due to interactions between the drugs, the effect of their combination can be larger (synergy) or smaller (antagony) than the sum of their individual effects (Bliss, 1939; Loewe, 1953; Berenbaum, 1989). Unfortunately, experimental examination of a wide range of combination therapy doses is costly and time-consuming.
Mathematical models can help in the effort to find optimal combination therapy doses. Within host mathematical models of influenza have previously been used to study many aspects of antiviral treatment including extracting of drug efficacy parameters (Beauchemin et al., 2008; Brown et al., 2011; Beggs and Dobrovolny, 2015; Liao et al., 2017), treatment of severe influenza (Dobrovolny et al., 2010, 2011; Deecke and Dobrovolny, 2018), emergence of drug resistance (Handel et al., 2007; Perelson et al., 2012; Hur et al., 2013; Canini et al., 2014; Dobrovolny and Beauchemin, 2017; Deecke and Dobrovolny, 2018), and to optimize antiviral treatments (Perelson et al., 2012; Heldt et al., 2013; Hur et al., 2013; Canini et al., 2014). While there are some mathematical models that attempt to model infections in patients by including an immune response (Dobrovolny et al., 2013; Cao and McCaw, 2015; Cao et al., 2015; Price et al., 2015; Zarnitsyna et al., 2016; Yan et al., 2017), the lack of appropriate human data for parameterizing and validating these models limits their predictive ability (Dobrovolny et al., 2013; Boianelli et al., 2015). However, simpler mathematical models can successfully reproduce in vitro dynamics (Beauchemin and Handel, 2011; Pinilla et al., 2012), and since mathematical models can quickly and efficiently simulate hundreds of combinations of doses, they are ideally suited as preliminary studies to ascertain whether combination therapy is effective and, if so, which combinations of doses produce the best results.
In this paper, we use an in vitro mathematical model of influenza infection to study combination therapy of influenza antivirals with different mechanisms of action. We measure the peak viral load, infection duration, and synergy/antagony of the various drug combinations to determine whether mechanisms of action pair better in combination therapy. We find that drugs that lengthen the eclipse phase and drugs that decrease the infection rate perform poorly, using all three measures, in combination with all other drugs.
2. Materials and Methods
We use a simple model of viral infection (Pinilla et al., 2012) given by the equations
In this model, virus V infects healthy target cells T at a rate β. Once infected, the cells move into the eclipse phase Ej during which there is internal replication of viral proteins and RNA, but no external production of virus. After some average time τE, the cells move into the infectious phase Ij where they are actively producing virus at rate p. Virus is cleared from the system at a clearance rate c. After an average time τI, infectious cells die. Both the eclipse and infectious phases are modeled as having Erlang distributions represented as the multiple (nE eclipse and nI infectious) compartments in each phase. Recent work has suggested that this distribution more faithfully reproduces all aspects of viral dynamics (Holder and Beauchemin, 2011; Kakizoe et al., 2015; Beauchemin et al., 2017).
In order to perform simulations with our model, we need estimates of the parameter values. We used parameter values determined by Pinilla et al. (2012) through fits of this model to data from in vitro infections of influenza A/Québec/144147/09 in MDCK cells. Parameter values are given in Table 1.
2.2. Modeling the Drug Effect
To model the effectiveness of influenza antivirals, we use the efficacy, ε, a parameter that varies from 0 to 1. An efficacy of 0 means the drug has no effect while an efficacy of 1 means the drug is completely effective. The efficacy of a drug is related to the drug dose through the Emax model (Holford and Sheiner, 1981),
where D is the dose of the antiviral, εmax is the maximum possible drug efficacy, IC50 is the drug dose needed to achieve half the maximum effect, and γ is the Hill coefficient. Biologically, γ is determined by the number of binding reactions needed for the drug to function (Weiss, 1997), which is assumed to be 1 for influenza antivirals (Beauchemin et al., 2008; Beggs and Dobrovolny, 2015). We also assume that εmax =1 so that we can explore the full range of possible behavior. Finally, we set IC50 =1 which amounts to expressing drug doses relative to the IC50, i.e., D → D/IC50.
In this study, we are interested in examining how antiviral mechanism of action affects the efficacy of combination therapy of influenza. Therefore, we examine not only currently available antivirals, but also hypothetical antivirals with other mechanisms of action. We model different mechanisms of action by applying drug efficacy to different parameters of the model. It is important to note that the model is quite general in that most biological processes are not explicitly represented. This means that it is not always clear which processes are captured by each parameter and some processes may affect more than one parameter. For example, the boundary between cell infection, represented by β, and the intracellular processes that are part of the eclipse phase, represented by τE, is not well-defined; fusion and endocytosis fall into this fuzzy boundary region and might well affect both model parameters.
The mechanisms of action for antiviral drugs that we modeled are:
• Reducing infection rate. This is modeled by applying the efficacy to the parameter β → (1−ε)β for β in both the target cell and eclipse cell equations. This represents an antiviral that blocks entry into the cell.
• Protecting target cells. To do this we place the efficacy on the parameter β → (1 − ε)β, but only for the β that appears in the eclipse equation. This represents an antiviral that blocks intracellular processes, but not entry into the cell. Cells containing virus are removed from the target cell class, but cannot be infected and so can be thought of as being “protected” from infection. This was shown to be the best mathematical model for replicating the effects of the influenza antiviral amantadine (Beauchemin et al., 2008). This mechanism will be denoted by β2.
• Reducing virion production rate. For this mechanism, we put the efficacy on p → (1 − ε)p. This is often used to model the action of neuraminidase inhibitors (Baccam et al., 2006; Handel et al., 2007; Dobrovolny et al., 2011).
• Increasing the rate of viral clearance. In this case, the efficacy is applied to c → c/(1 − ε). In patients, this most likely represents a drug that stimulates the adaptive immune response, particularly antibodies, to enhance clearance of the virus (Taylor and Dimmock, 1985a,b), but could also represent a drug that inactivates virus (Fujimori et al., 2012).
• Increasing the length of the eclipse phase. A drug effect on parameter τE → τE/(1 − ε) represents an antiviral that delays assembly of the virions. There could be several possible mechanisms for this (Heldt et al., 2013), such as delay of production of proteins or RNA.
• Decreasing the lifespan of infectious cells. A drug effect on parameter τI → (1 − ε)τI could also represent a stimulant of the immune response in patients, although in this case increasing the cytotoxic T lymphocytes which are responsible for killing infected cells (Zweerink et al., 1977; Mbawuike et al., 2007), but could also represent a drug that stimulates autophagy of infected cells (Feizi et al., 2017).
2.3. Measuring the Effect of Combination Therapy
We use our model to determine which combination of antivirals are the best at treating influenza. We use two measurements of the viral titer to assess the efficacy of treatment, the maximum viral titer Vmax and infection duration, Tinf. The maximum viral titer is indicative of the viral burden in patients and is thought to be a measure of the transmissibility of the infection (Handel et al., 2009). The infection duration is defined as the time the viral titer remains above , as described in Dobrovolny et al. (2010).
A quantity that is often used to characterize combination therapy is synergy or antagony. If the effect of a combination of antivirals is better than expected based on the individual effects, the combination is said to be synergistic; if the effect of the combination is worse than expected, the combination is said to be antagonistic. If we assume a multiplicative effect for drugs, known as Bliss synergy (Bliss, 1939), this is calculated via
where W is the observed percent inhibition of viral titer due to the combination of two drugs, and M and N are the observed percent inhibition of each drug individually (Koizumi and Iwami, 2014). A positive value of S represents a synergistic combination, and a negative value represents antagony.
An alternative assumption, known as Loewe additivity (Loewe, 1953), assumes that the effect of drug combinations is additive. Synergy in this case is measured using the combination index (CI)
where D1 and D2 are the doses of the two drugs that result in a particular effect during monotherapy, and and are the doses of the two drugs in combination that result in the same effect. If the CI is greater than 1, the drugs are considered synergistic, while if the CI is less than 1, the drugs are considered antagonistic. Note that for these calculations, viral titer is measured at a specific time; we used 72 h, a time that is commonly used in experimental measurements of synergy (Smee et al., 2002, 2009; Ilyushina et al., 2006, 2008; Nguyen et al., 2009, 2010).
3.1. Curing the Infection
For ordinary differential equation models, whether or not an infection progresses is determined by the basic reproductive number, R0, which is defined as the number of secondary infections produced by a single infected cell in a heterogeneous population. R0 for this model is given by
If R0 > 1, we will have an infection and if R0 < 1, infection is suppressed. Modeling the application of drug as described in the Methods section results in modification of R0, such that R0, treated = (1 − ε1)(1 − ε2)R0, where εi are the efficacies of each of the two drugs being used in combination. Applying the condition that R0, treated = 1, and keeping the assumption that the Hill coefficient is 1 (γ=1), we can find the boundary between combinations of doses that prevent infection and those that don't.
where i, j = 1, 2 indicating each of the two drugs in the combination, and f(i) is given by
Using Equation (2), with our assumptions of εmax=1 and IC50=1, to write the condition in terms of drug dose, we have
Note that this analysis does not capture the effect of a drug effect applied to τE since τE does not appear in R0.
The mathematical model (Equation 1) can be computationally solved, giving a picture of how viral load changes as different drugs are applied. In Figure 1, we show the model predictions of treatment with an antiviral that reduces virion production rate (red line), an antiviral that reduces infection rate (blue line), and their combination (green line). The three panels show antiviral doses of 10 × IC50 (left), 100 × IC50 (center), and 1,000 × IC50. We see that a drug that reduces infection has little effect on the viral titer, shifting the curve slightly to the right at high doses, but not affecting the peak viral load much. The drug that reduces virion production has a larger effect. At a dose of 10 × IC50, the combination regimen is no different than monotherapy with a drug that reduces production of virions. At higher doses, the combination treatment reflects the reduced peak viral load of the drug reducing virion production and the rightward shift of a drug that reduces infection rate.
Figure 1. Model simulations of antiviral treatment. The untreated infection is given by the black dashed line. Monotherapy with an antiviral reducing virion production is shown in red. Monotherapy with an antiviral reducing the infection rate is shown in blue. The combination of the two (at the same doses) is shown in green. The left panel shows infections when the antiviral is treated with doses of 10 × IC50; the center panel shows doses of 100 × IC50; and the right panel shows doses of 1,000 × IC50.
To get a more complete view of drug effect, we investigate the effect of combination therapy on measures of disease burden, Vmax and Tinf. We examine possible combinations of theoretical antivirals (14 pairs), excluding combinations of drugs modeled as acting on β and β2. For each dose combination, we simulate an infection and measure Vmax and Tinf. Results for Vmax are shown in Figure 2 and results for Tinf are shown in Figure 3. The boundary between infection and no infection, given by Equation (8), is indicated by the black line on graphs in Figure 2. Note that while we investigate the effect of antivirals up to doses 10, such high doses are not achievable for many antivirals due to toxicity.
Figure 2. Maximum viral titer for a range of doses of various combinations of influenza antivirals. Lowest peak viral titer is in blue with highest peak values indicated in red. The black line indicates the theoretical curve given by Equation (8) that delineates regions where the infection is cured.
Figure 3. Infection duration for a range of doses of various combinations of influenza antivirals. Regions where the infection has been cured, (duration=0) are indicated by dark blue, while regions where the infection is longer than 10 days are indicated by red.
While the boundary between infection and no infection is the same for most drug combinations, we can see from Figure 2 that different mechanisms of action of the antivirals affect how quickly viral load is reduced as we approach the boundary. For example, a drug that reduces infection rate (reducing β) does not assist in lowering viral load no matter which other drug is used in combination. This is seen in the first four graphs of Figure 2 which show no change in viral load as dose on β (y-axis) is varied. We can also see that drugs that increase the duration of the eclipse phase appear to reduce the viral load to zero beyond a certain dose (about 10 IC50), although this is just a computational effect. A drug that increases the duration of the eclipse phase simply delays the time of the peak viral titer (González-Parra and Dobrovolny, 2018), in this case pushing the peak beyond the duration of the simulation. An infection that grows that slowly, however, is likely to be suppressed by the patient's immune response (Beauchemin and Handel, 2011). The remaining drug combinations, those not including a drug reducing infection rate or a drug increasing eclipse duration, all have a similar effect on peak viral load.
Since the effect of different drug combinations on viral load is so similar, we look for other measures of infection severity that might help differentiate certain drug combinations. The infection duration is shown in Figure 3. We again see boundaries, for all drug combinations, beyond which Tinf falls to 0 and there is no infection. For combinations that include a drug that reduces infection rate, Tinf increases quite drastically if the dose of either drug is just below this boundary. We see a similar increasing Tinf when a drug that increases the eclipse phase is part of the combination due to the delay of the replication cycle as the eclipse phase increases. The remaining drug combinations show gradual decrease of Tinf as drug dose of either drug in the combination is increased. We can, however, see some slight differences in how quickly Tinf decreases. We find that the combination of a drug increasing the clearance rate and a drug decreasing the infectious lifespan reduces Tinf is the most effective at reducing the duration of the infection since it requires lower doses than the other combinations to achieve short infection durations.
Of particular interest is whether a drug combination improves outcomes over treatment with a single drug. There are two models for drug interaction; one assumes a multiplicative effect (Bliss synergy) (Bliss, 1939), and one assumes an additive effect (Loewe, 1953). Since we are not investigating specific antivirals, we cannot make assumptions about how the antivirals might interact, so we include information about both types of synergy here.
The Bliss synergy/antagony of the different drug combinations is shown in Figure 4. In these figures, blue indicates antagony, red indicates synergy, and green indicates no enhancement or impediment of the antiviral effect. Drug combinations that include an antiviral that decreases infection rate show antagony over a broad range of doses. Drug combinations that include a drug that increases the eclipse duration also have large regions of antagony. The remaining drug combinations are largely neutral, showing little synergy or antagony over most combinations of doses. Particularly interesting is the combination of an antiviral that reduces infection rate and one that prolongs the eclipse phase (second row, left panel) where there is high synergy at low doses of the drug that prolongs the eclipse phase and antagony at high doses of this antiviral. Remember that synergy is a measure of whether the drug combination had a larger effect than expected. At low doses, the drugs are not expected to have much effect, so as long as the combination outperforms that low expectation, synergy will be high. Both of these antivirals shift the viral titer curve to the right. On their own, they will not shift the curve much (even up to very high doses for a drug that reduces infection rate, as seen in Figure 1), but combined they shift the curve just enough to have the peak move past the 72 h measurement time resulting in a large drop in the viral titer measurement.
Figure 4. Synergy for a range of doses of various combinations of influenza antivirals. Antagony of the antivirals is indicated by regions of blue, while regions of synergy are indicated by red. Note that most combinations result in neither synergy nor antagony (pale green).
The combination index of the different drug combinations is shown in Figure 5. Like the previous figure, blue indicates antagony, while red indicates synergy. In this case, no effect is given by CI = 1. If we assume Loewe additivity for our antivirals, we see many more regions of synergy than for Bliss synergy. There are only five drug combinations (β1/p, β1/c, β2/p, β2/c, p/c) where the majority of dose combinations show antagony or no effect. For the remaining drug combinations, the CI varies rapidly, going from anatagony to synergy over a very small range of doses.
Figure 5. Combination index for a range of doses of various combinations of influenza antivirals. Antagony of the antivirals is indicated by regions of blue, while regions of synergy are indicated by red. A combination index of one (no effect) is pale green.
In this paper, we computationally evaluated combinations of theoretical influenza antivirals to assess the effect of antiviral mechanism of action on drug interactions. We found that combination therapies that include drugs that reduce the infection rate or combination therapies that include a drug that lengthens the eclipse phase show markedly different effects on peak viral load, infection duration, and synergy than other combination therapies. The decreased efficacy in reducing viral load, dose combinations that lead to increased infection duration, and large number of dose combinations that produce antagonistic reactions suggest that antivirals with these two mechanisms of action are not good candidates for combination therapy. Combination therapies using antivirals with other mechanisms of action have similar effects on peak viral load, infection duration, and synergy, although a combination of an antiviral that increases viral clearance and an antiviral that decreases infectious cell lifespan reduces infection duration at lower doses than other drug combinations, perhaps making it the best combination. Note that both of these antivirals decrease the duration of the viral decay phase of influenza. Since this phase is typically much longer than the viral growth phase (Smith et al., 2010; Beauchemin and Handel, 2011), a fractional reduction of the duration of this phase results in a greater shortening of the infection than a similar fractional reduction in the duration of the growth phase.
While we have identified a “best” antiviral combination by using infection duration, other measurements could lead to a different conclusion. Other features of the viral titer curve, besides peak viral load or some measure of infection duration, are sometimes used to assess the effectiveness of antiviral treatment (Beggs and Dobrovolny, 2015). These include the area under the viral titer curve (AUC) (Ryan et al., 1994; Hayden et al., 1999; Heldens and van den Hoven, 2002; Rayner et al., 2013; Beggs and Dobrovolny, 2015) or the viral growth and decay rates (Beggs and Dobrovolny, 2015). Measurement of these features might allow further differentiation of the role of antiviral mechanism of action in combination therapy, although the value of these measures in terms of clinical benefit to the patient are unclear. The most common method of assessing efficacy of an antiviral (or combination of antivirals) is to measure the viral load at a specific time (Noah et al., 2007), although the efficacy estimated in this manner is known to depend on the measurement time (Stresser et al., 2014).
Besides looking at features of the viral titer curve, other factors could play a role in deciding which antiviral combinations are most beneficial to patients. Combination therapy has been proposed as a strategy to decrease the occurrence of drug resistance (Dunning et al., 2014). A recent modeling study suggests that mechanism of action of an antiviral plays a role in how quickly drug resistant mutants emerge during an influenza infection (Dobrovolny and Beauchemin, 2017) during monotherapy. This effect could carry over to combination therapies such that certain antiviral combinations will be more effective in blocking the appearance of drug resistant mutants, a factor that should be considered when determining optimal drug combinations for treating influenza. The toxicity/side-effects of antivirals also needs to be considered when determining optimal treatment strategies. Other factors that could be considered when designing drug combinations are the cost or cost-benefit ratio of the medications (Burch et al., 2009) and toxicity or side-effects of the medications.
The mathematical model used in this study is fairly simple and does not include a full description of all biological processes involved in influenza replication. We do not include an explicit immune response since there is still no consensus on the correct mathematical formulation of immune responses (Dobrovolny et al., 2013), although this will be needed to properly assess combination therapy in vivo. A recent modeling study has indicated that inclusion of an immune response alters the predicted effect of antivirals (Cao and McCaw, 2015). Inclusion of an immune response would also allow for more accurate representation of immune-stimulating antivirals (Zweerink et al., 1977; Taylor and Dimmock, 1985a,b; Mbawuike et al., 2007). More broadly, due to the generality of the model, we also do not fully capture the full range of possible mechanisms of action of influenza antivirals (Heldt et al., 2013; Liao et al., 2017), although this could be corrected by using more detailed models if combination therapy of specific antivirals needs to be investigated. While this study might not be detailed enough to make predictions about specific antiviral combinations, the model used here is sufficient to show that antiviral mechanisms of action affect the effectiveness of combination therapy. Clinicians and drug developers should consider interactions between different mechanisms of action when developing combination therapies.
Data Availability Statement
The datasets generated for this study can be found in github [https://github.com/hdobrovo/combination_therapy.git].
HD conceived the experiments. KM and HD conducted the experiments. KM, TR, and HD analyzed the results, KM and HD wrote the manuscript. All authors reviewed and approved the manuscript.
This work was supported by a grant from the National Science Foundation Research Experience for Undergraduate program PHY-1358770.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Abed, Y., Bourgault, A.-M., Fenton, R. J., Morley, P. J., Gower, D., Owens, I. J., et al. (2002). Characterization of 2 influenza A(H3N2) clinical isolates with reduced susceptibility to neuraminidase inhibitors due to mutations in the hemagglutinin gene. J. Infect. Dis. 186, 1074–1080. doi: 10.1086/344237
Abed, Y., Goyette, N., and Boivin, G. (2005). Generation and characterization of recombinant influenza A (H1N1) viruses harboring amantadine resistance mutations. Antimicrob. Agents Chemother. 49, 556–559. doi: 10.1128/AAC.49.2.556-559.2005
Bantia, S., Kellogg, D., Parker, C. D., and Babu, Y. (2010). Combination of peramivir and rimantadine demonstrate synergistic antiviral effects in sub-lethal influenza A (H3N2) virus mouse model. Antiviral Res. 88, 276–280. doi: 10.1016/j.antiviral.2010.09.020
Beauchemin, C. A., and Handel, A. (2011). A review of mathematical models of influenza A infections within a host or cell culture: lessons learned and challenges ahead. BMC Public Health 11(Suppl. 1):S7. doi: 10.1186/1471-2458-11-S1-S7
Beauchemin, C. A., McSharry, J. J., Drusano, G. L., Nguyen, J. T., Went, G. T., Ribeiro, R. M., et al. (2008). Modeling amantadine treatment of influenza A virus in vitro. J. Theor. Biol. 254, 439–451. doi: 10.1016/j.jtbi.2008.05.031
Beauchemin, C. A., Miura, T., and Iwami, S. (2017). Duration of SHIV production by infected cells is not exponentially distributed: implications for estimates of infection parameters and antiviral efficacy. Sci. Rep. 7:42765. doi: 10.1038/srep42765
Beigel, J. H., Bao, Y., Beeler, J., Manosuthi, W., Slandzicki, A., Dar, S. M., et al. (2017). Oseltamivir, amantadine, and ribavirin combination antiviral therapy versus oseltamivir monotherapy for the treatment of influenza: a multicentre, double-blind, randomised phase 2 trial. Lancet Infect. Dis. 17, 1255–1265. doi: 10.1016/S1473-3099(17)30476-0
Belardo, G., Cenciarelli, O., La Frazia, S., Rossignol, J. F., and Santoro, M. G. (2015). Synergistic effect of nitazoxanide with neuraminidase inhibitors against influenza A viruses in vitro. Antimicrob. Agents Chemother. 59, 1061–1069. doi: 10.1128/AAC.03947-14
Boianelli, A., Nguyen, V. K., Ebensen, T., Schulze, K., Wilk, E., Sharma, N., et al. (2015). Modeling influenza virus infection: a roadmap for influenza research. Viruses 7, 5274–5304. doi: 10.3390/v7102875
Bright, R. A., Shay, D. K., Shu, B., Cox, N. J., and Klimov, A. I. (2006). Adamantane resistance among influenza A viruses isolated early during the 2005-2006 influenza season in the United States. JAMA 295, 891–894. doi: 10.1001/jama.295.8.joc60020
Brown, A. N., Bulitta, J. B., McSharry, J. J., Weng, Q., Adams, J. R., Kulawy, R., et al. (2011). Effect of half-life on the pharmacodynamic index of zanamivir against influenza virus delineated by a mathematical model. Antimicrob. Agents Chemother. 55, 1747–1753. doi: 10.1128/AAC.01629-10
Burch, J., Paulden, M., Conti, S., Stock, C., Corbett, M., Welton, N., et al. (2009). Antiviral drugs for the treatment of influenza: a systematic review and economic evaluation. Health Tech. Assess. 13, 1–265. doi: 10.3310/hta13580
Canini, L., Conway, J. M., Perelson, A. S., and Carrat, F. (2014). Impact of different oseltamivir regimens on treating influenza A virus infection and resistance emergence: insights from a modelling study. PLoS Comput. Biol. 10:e1003568. doi: 10.1371/journal.pcbi.1003568
Cao, P., and McCaw, J. M. (2015). The mechanisms for within-host influenza virus control affect model-based assessment and prediction of antiviral treatment. Viruses - Basel 9:197. doi: 10.3390/v9080197
Cao, P., Yan, A. W., Heffernan, J. M., Petrie, S., Moss, R. G., Carolan, L. A., et al. (2015). Innate immunity and the inter-exposure interval determine the dynamics of secondary influenza virus infection and explain observed viral hierarchies. PLoS Comput. Biol. 11:e1004334. doi: 10.1371/journal.pcbi.1004334
Chen, L., Dou, J., Su, Z., Zhou, H., Wang, H., Zhou, W., et al. (2011). Synergistic activity of baicalen with ribavirin against influenza A (H1N1) virus infections in cell culture and mice. Antivir. Res. 91, 314–320. doi: 10.1016/j.antiviral.2011.07.008
de Mello, C. P., Drusano, G. L., Adams, J. R., Shudt, M., Kulawy, R., and Brown, A. N. (2018). Oseltamivir-zanamivir combination therapy suppresses drug-resistant H1N1 influenza A viruses in the hollow fiber infection model (HFIM) system. Eur. J. Pharmaceut. Sci. 111, 443–449. doi: 10.1016/j.ejps.2017.10.027
Dharan, N. J., Gubareva, L. V., Meyer, J. J., Okomo-Adhiambo, M., McClinton, R. C., Marshall, S. A., et al. (2009). Local and systemic cytokine responses during experimental human influenza A virus infection. JAMA 301, 1034–1041. doi: 10.1001/jama.2009.294
Dobrovolny, H. M., Baron, M. J., Gieschke, R., Davies, B. E., Jumbe, N. L., and Beauchemin, C. A. A. (2010). Exploring cell tropism as a possible contributor to influenza infection severity. PLoS ONE 5:e13811. doi: 10.1371/journal.pone.0013811
Dobrovolny, H. M., and Beauchemin, C. A. (2017). Modelling the emergence of influenza drug resistance: the roles of surface proteins, the immune response and antiviral mechanisms. PLoS ONE 12:e0180582. doi: 10.1371/journal.pone.0180582
Dobrovolny, H. M., Gieschke, R., Davies, B. E., Jumbe, N. L., and Beauchemin, C. A. A. (2011). Neuraminidase inhibitors for treatment of human and avian strain influenza: a comparative study. J. Theor. Biol. 269, 234–244. doi: 10.1016/j.jtbi.2010.10.017
Dobrovolny, H. M., Reddy, M. B., Kamal, M. A., Rayner, C. R., and Beauchemin, C. A. (2013). Assessing mathematical models of influenza infections using features of the immune response. PLoS ONE 8:e57088. doi: 10.1371/journal.pone.0057088
Dong, G., Peng, C., Luo, J., Wang, C., Han, L., Wu, B., et al. (2015). Adamantane-resistant influenza A viruses in the world (1902-2013): Frequency and distribution of M2 gene mutations. PLoS ONE 10:e0119115. doi: 10.1371/journal.pone.0119115
Duval, X., van der Werf, S., Blanchon, T., Mosnier, A., Bouscambert-Duchamp, M., Tibi, A., et al. (2010). Efficacy of oseltamivir-zanamivir combination compared to each monotherapy for seasonal influenza: a randomized placebo-controlled trial. PLoS Med. 7:e1000362. doi: 10.1371/journal.pmed.1000362
Feizi, N., Mehrbod, P., Romani, B., Soleimanjahi, H., Bamdad, T., Feizi, A., et al. (2017). Autophagy induction regulates influenza virus replication in a time-dependent manner. J. Med. Microbiol. 66, 536–541. doi: 10.1099/jmm.0.000455
Fujimori, Y., Sato, T., Hayata, T., Nagao, T., Nakayama, M., Nakayama, T., et al. (2012). Novel antiviral characteristics of nanosized copper(I) iodide particles showing inactivation activity against 2009 pandemic H1N1 influenza virus. Appl. Env. Microbiol. 78, 951–955. doi: 10.1128/AEM.06284-11
Govorkova, E. A., Fang, H.-B., Tan, M., and Webster, R. G. (2004). Neuraminidase inhibitor-rimantadine combinations exert additive and synergistic anti-influenza virus effects in MDCK cells. Antimicrob. Agents Chemother. 48, 4855–4863. doi: 10.1128/AAC.48.12.4855-4863.2004
Haasbach, E., Hartmayer, C., and Planz, O. (2013). Combination of MEK inhibitors and oseltamivir leads to synergistic antiviral effects after influenza A virus infection in vitro. Antiviral Res. 98, 319–324. doi: 10.1016/j.antiviral.2013.03.006
Handel, A., Longini, I. M., and Antia, R. (2009). Antiviral resistance and the control of pandemic influenza: the roles of stochasticity, evolution and model details. J. Theor. Biol. 256, 117–125. doi: 10.1016/j.jtbi.2008.09.021
Hayden, F. G., Cote, K. M., and Gordon Douglas, J. (1980). Plaque inhibition assay for drug susceptibility testing of influenza viruses. Antimicrob. Agents Chemother. 17, 865–870. doi: 10.1128/AAC.17.5.865
Hayden, F. G., Treanor, J. J., Fritz, R. S., Lobo, M., Betts, R. F., Miller, M., et al. (1999). Use of the oral neuraminidase inhibitor oseltamivir in experimental human influenza: randomized controlled trials for prevention and treatment. JAMA 282, 1240–1246. doi: 10.1001/jama.282.13.1240
Heldens, J., and van den Hoven, M. W. R. (2002). Area under the curve calculations as a tool to compare the efficacy of equine influenza vaccines — a retrospective analysis of three independent field trials. J. Immunol. Meth. 254, 11–17. doi: 10.1016/S0022-1759(01)00571-3
Heldt, F. S., Frensing, T., Pflugmacher, A., Gropler, R., Peschel, B., and Reichl, U. (2013). Multiscale modeling of influenza A virus infection supports the development of direct-acting antivirals. PLoS Comput. Biol. 9:e1003372. doi: 10.1371/journal.pcbi.1003372
Holder, B. P., and Beauchemin, C. A. (2011). Exploring the effect of biological delays in kinetic models of influenza within a host or cell culture. BMC Public Health 11:S10. doi: 10.1186/1471-2458-11-S1-S10
Holford, N., and Sheiner, L. (1981). Understanding the dose-effect relationship: clinical application of pharmacokinetic-pharmacodynamic models. Clin. Pharmacokinet. 6, 429–453. doi: 10.2165/00003088-198106060-00002
Hoopes, J. D., Driebe, E. M., Kelley, E., Engelthaler, D. M., Keim, P. S., Perelson, A. S., et al. (2011). Triple combination antiviral drug (TCAD) composed of amantadine, oseltamivir, and ribavirin impedes the selection of drug-resistant influenza A virus. PLoS ONE 6:e29778. doi: 10.1371/journal.pone.0029778
Hur, K.-Y., Moon, J.-Y., Kim, S.-H., and Yoo, J.-Y. (2013). Model-based simulation and prediction of an antiviral strategy against influenza A infection. PLoS ONE 8:e68235. doi: 10.1371/journal.pone.0068235
Ilyushina, N. A., Bovin, N. V., Webster, R. G., and Govorkova, E. A. (2006). Combination chemotherapy, a potential strategy for reducing the emergence of drug-resistant influenza A variants. Antiviral Res. 70, 121–131. doi: 10.1016/j.antiviral.2006.01.012
Ilyushina, N. A., Hay, A., Yilmaz, N., Boon, A. C. M., Webster, R. G., and Govorkova, E. A. (2008). Oseltamivir-ribavirin combination therapy for highly pathogenic H5N1 influenza virus infection in mice. Antimicrob. Agents Chemother. 52, 3889–3897. doi: 10.1128/AAC.01579-07
Ilyushina, N. A., Hoffmann, E., Salomon, R., Webster, R. G., and Govorkova, E. A. (2007). Amantadine-oseltamivir combination therapy for H5N1 influenza virus infection in mice. Antiviral Ther. 12, 363–370.
Kakizoe, Y., Nakaoka, S., Beauchemin, C. A., Morita, S., Mori, H., Igarashi, T., et al. (2015). A method to determine the duration of the eclipse phase for in vitro infection with a highly pathogenic SHIV strain. Sci. Rep. 5:10371. doi: 10.1038/srep10371
Khandaker, G., Dierig, A., Rashid, H., King, C., Heron, L., and Booy, R. (2011). Systematic review of clinical and epidemiological features of the pandemic influenza A (H1N1) 2009. Influenza Other Respir Viruses 5, 148–156. doi: 10.1111/j.1750-2659.2011.00199.x
Kim, W.-Y., Suh, G. Y., Huh, J. W., Kim, S.-H., ju Kim, M., Kim, Y. S., et al. (2011). Triple-combination antiviral drug for pandemic H1N1 influenza virus infection in critically ill patients on mechanical ventilation. Antimicrob. Agents Chemother. 55, 5703–5709. doi: 10.1128/AAC.05529-11
Koizumi, Y., and Iwami, S. (2014). Mathematical modeling of multi-drugs therapy: a challenge for determining the optimal combinations of antiviral drugs. Theor. Biol. Med. Model. 11:41. doi: 10.1186/1742-4682-11-41
Koszalka, P., Tilmanis, D., and Hurt, A. C. (2017). Synthesis and in vitro study of novel neuraminidase inhibitors against avian influenza virus. Influenza Other Respir. Vir. 11, 240–246. doi: 10.1111/irv.12446
Liao, L. E., Kowal, S., Cardenas, D. A., and Beauchemin, C. A. A. (2017). Exploring virus release as a bottleneck for the spread of influenza a virus infection in vitro and the implications for antiviral therapy with neuraminidase inhibitors. PLoS ONE 12:e0183621. doi: 10.1371/journal.pone.0183621
Loregian, A., Mercorelli, B., Nannetti, G., Compagnin, C., and Palu, G. (2014). Antiviral strategies against influenza virus: towards new therapeutic approaches. Cell. Mol. Life Sci. 71, 3659–3683. doi: 10.1007/s00018-014-1615-2
Lu, Y., Hardes, K., Dahms, S. O., Bottcher-Friebertshauser, E., Steinmetzer, T., Than, M. E., et al. (2015). Peptidomimetic furin inhibitor MI-701 in combination with oseltamivir and ribavirin efficiently blocks propagation of highly pathogenic avian influenza viruses and delays high level oseltamivir resistance in MDCK cells. Antiviral Res. 120, 89–100. doi: 10.1016/j.antiviral.2015.05.006
Marathe, B. M., Wong, S.-S., Vogel, P., Garcia-Alcalde, F., Webster, R. G., Webby, R. J., et al. (2016). Combinations of oseltamivir and T-705 extend the treatment window for highly pathogenic influenza A(H5N1) virus infection in mice. Sci. Rep. 6:26742. doi: 10.1038/srep26742
Masihi, K., Schweiger, B., Finsterbusch, T., and Hengel, H. (2007). Low dose oral combination chemoprophylaxis with oseltamivir and amantadine for influenza A virus infections in mice. J. Chemother. 19, 295–303. doi: 10.1179/joc.2007.19.3.295
Mbawuike, I. N., Zhang, Y., and Couch, R. B. (2007). Control of mucosal virus infection by influenza nucleoprotein-specific CD8+ cytotoxic T lymphocytes. Resp. Res. 8, 44–52. doi: 10.1186/1465-9921-8-44
Morokutti-Kurz, M., Konig-Schuster, M., Koller, C., Graf, C., Graf, P., Kirchoff, N., et al. (2015). The intranasal application of zanamivir and carrageenan is synergistically active against influenza A virus in the murine model. PLoS ONE 10:e0128794. doi: 10.1371/journal.pone.0128794
Nguyen, J. T., Hoopes, J. D., Le, M. H., Smee, D. F., Patick, A. K., Faix, D. J., et al. (2010). Triple combination of amantadine, ribavirin, and oseltamivir is highly active and synergistic against drug resistant influenza virus strains in vitro. PLoS ONE 5:e9332. doi: 10.1371/journal.pone.0009332
Nguyen, J. T., Hoopes, J. D., Smee, D. F., Prichard, M. N., Driebe, E. M., Engelthaler, D. M., et al. (2009). Triple combination of oseltamivir, amantadine, and ribavirin displays synergistic activity against multiple influenza virus strains in vitro. Antimicrob. Agents Chemother. 53, 4115–4126. doi: 10.1128/AAC.00476-09
Nguyen, J. T., Smee, D. F., Barnard, D. L., Julander, J. G., Gross, M., de Jong, M. D., et al. (2012). Efficacy of combined therapy with amantadine, oseltamivir, and ribavirin in vivo against susceptible and amantadine-resistant influenza A viruses. PLoS ONE 7:e31006. doi: 10.1371/journal.pone.0031006
Noah, J. W., Severson, W., Noah, D. L., Rasmussen, L., White, E. L., and Jonsson, C. B. (2007). A cell-based luminescence assay is effective for high-throughput screening of potential influenza antivirals. Antivir. Res. 73, 50–59. doi: 10.1016/j.antiviral.2006.07.006
Pavlova, E. L., Simeonova, L. S., and Gegova, G. A. (2018). Combined efficacy of oseltamivir, isoprinosine and ellagic acid in influenza A(H3N2)-infected mice. Biomed. Pharmacother. 98, 29–35. doi: 10.1016/j.biopha.2017.12.014
Perelson, A. S., Rong, L., and Hayden, F. G. (2012). Combination antiviral therapy for influenza: predictions from modeling of human infections. J. Infect. Dis. 205, 1642–1645. doi: 10.1093/infdis/jis265
Pinilla, L. T., Holder, B. P., Abed, Y., Boivin, G., and Beauchemin, C. A. A. (2012). The H275Y neuraminidase mutation of the pandemic A/H1N1 influenza virus lengthens the eclipse phase and reduces viral output of infected cells, potentially compromising fitness in ferrets. J. Virol. 86, 10651–10660. doi: 10.1128/JVI.07244-11
Popov, A. F., Shchelkanov, M. Y., Dmitrenko, K. A., and Simakova, A. I. (2018). Combined therapy of influenza with antiviral drugs with a different mechanism of action in comparison with monotherapy. J. Pharm. Sci. Res. 10, 357–360.
Price, I., Mochan-Keef, E. D., Swigon, D., Ermentrout, B. G., Lukens, S., Toapanta, F. R., et al. (2015). The inflammatory response to influenza A virus (H1N1): an experimental and mathematical study. J. Theor. Biol. 374, 83–93. doi: 10.1016/j.jtbi.2015.03.017
Rayner, C. R., Bulik, C. C., Kamal, M. A., Reynolds, D. K., Toovey, S., Hammel, J. P., et al. (2013). Pharmacokinetic-pharmacodynamic determinants of oseltamivir efficacy using data from phase 2 inoculation studies. Antimicrob. Agents Chemother. 57, 3478–3487. doi: 10.1128/AAC.02440-12
Ryan, D. M., Ticehurst, J., Dempsey, M. H., and Penn, C. R. (1994). Inhibition of influenza virus replication in mice by GG167 (4-guanidino-2,4-dideoxy-2,3-dehydro-N-acetylneuraminic acid) is consistent with extracellular activity of viral neuraminidase (sialidase). Antimicrob. Agents Chemother. 38, 2270–2275. doi: 10.1128/AAC.38.10.2270
Seo, S., Englund, J. A., Nguyen, J. T., Pukrittayakamee, S., Lindegardh, N., Tarning, J., et al. (2013). Combination therapy with amantadine, oseltamivir and ribavirin for influenza a infection: safety and pharmacokinetics. Antivir. Ther. 18, 377–386. doi: 10.3851/IMP2475
Smee, D., Hurst, B., Wong, M. H., Bailey, K., Tarbet, E. B., Morrey, J. D., et al. (2010a). Effects of the combination of favipiravir (T-705) and oseltamivir on influenza A virus infection in mice. Antimicrob. Agents Chemother. 54, 126–133. doi: 10.1128/AAC.00933-09
Smee, D., Hurst, B., Wong, M. H., Tarbet, E., Babu, Y., Klumpp, K., et al. (2010b). Combinations of oseltamivir and peramivir for the treatment of influenza A (H1N1) virus infections in cell culture and in mice. Antivir. Res. 88, 38–44. doi: 10.1016/j.antiviral.2010.07.003
Smee, D. F., Bailey, K. W., Morrison, A. C., and Sidwell, R. W. (2002). Combination treatment of influenza A virus infections in cell culture and in mice with the cyclopentane neuraminidase inhibitor RWJ-270201 and ribavarin. Chemotherepy 48, 88–93. doi: 10.1159/000057668
Smee, D. F., Hurst, B. L., Wong, M.-H., Bailey, K. W., and Morrey, J. D. (2009). Effects of double combinations of amantadine, oseltamivir, and ribavirin on influenza A (H5N1) virus infections in cell culture and in mice. Antimicrob. Agents Chemother. 53, 2120–2128. doi: 10.1128/AAC.01012-08
Stresser, D. M., Mao, J., Kenny, J. R., Jones, B. C., and Grime, K. (2014). Exploring concepts of in vitro time-dependent CYP inhibition assays. Exp. Opin. Drug Metabol. Toxicol. 10, 157–174. doi: 10.1517/17425255.2014.856882
Tarbet, E. B., Vollmer, A. H., Hurst, B. L., Barnard, D. L., Furuta, Y., and Smee, D. F. (2014). In vitro activity of favipiravir and neuraminidase inhibitor combinations against oseltamivir-sensitive and oseltamivir-resistant pandemic influenza A (H1N1) virus. Arch. Virol. 159, 1279–1291. doi: 10.1007/s00705-013-1922-1
Taylor, H. P., and Dimmock, N. J. (1985a). Mechanism of neutralization of influenza virus by secretory IgA is different from that of monomeric IgA or IgG. J. Exp. Med. 161, 198–209. doi: 10.1084/jem.161.1.198
White, K. M., Abreu, P., Wang, H., De Jesus, P. D., Manicassamy, B., Garcia-Sastre, A., et al. (2018). Broad spectrum inhibitor of influenza A and B viruses targeting the viral nucleoprotein. ACS Infect. Dis. 4, 146–157. doi: 10.1021/acsinfecdis.7b00120
Yan, A. W., Cao, P., Heffernan, J. M., McVernon, J., Quinn, K. M., La Gruta, N. L., et al. (2017). Modelling cross-reactivity and memory in the cellular adaptive immune response to influenza infection in the host. J. Theor. Biol. 413, 34–49. doi: 10.1016/j.jtbi.2016.11.008
Zabrodskaya, Y. A., Lebedev, D. V., Egorova, M. A., Shaldzhyan, A. A., Shvetsov, A. V., Kuklin, A. I., et al. (2018). The amyloidogenicity of the influenza virus PB1-derived peptide sheds light on its antiviral activity. Biophys. J. 234, 16–23. doi: 10.1016/j.bpc.2018.01.001
Zaraket, H., Saito, R., Suzuki, Y., Baranovich, T., Dapat, C., Caperig-Dapat, I., et al. (2010). Genetic makeup of amantadine-resistant and oseltamivir-resistant human influenza A/H1N1 viruses. J. Clin. Microbiol. 48, 1085–1092. doi: 10.1128/JCM.01532-09
Zarnitsyna, V. I., Handel, A., McMaster, S. R., Hayward, S. L., Kohlmeier, J. E., and Antia, R. (2016). Mathematical model reveals the role of memory CD8 T cell populations in recall responses to influenza. Front. Immunol. 7:165. doi: 10.3389/fimmu.2016.00165
Zhang, J., Hu, Y., Musharrafieh, R. G., Yin, H., and Wang, J. (2018). Focusing on the influenza virus polymerase complex: recent progress in drug discovery and assay development. Curr. Med. Chem. doi: 10.2174/0929867325666180706112940. [Epub ahead of print].
Keywords: influenza, antiviral, mathematical modeling, combination therapy, mechanism of action
Citation: Melville K, Rodriguez T and Dobrovolny HM (2018) Investigating Different Mechanisms of Action in Combination Therapy for Influenza. Front. Pharmacol. 9:1207. doi: 10.3389/fphar.2018.01207
Received: 29 June 2018; Accepted: 03 October 2018;
Published: 23 October 2018.
Edited by:Mauro Martins Teixeira, Universidade Federal de Minas Gerais, Brazil
Reviewed by:Joana Rocha-Pereira, KU Leuven, Belgium
Tero Aittokallio, Institute for Molecular Medicine Finland, Finland
Copyright © 2018 Melville, Rodriguez and Dobrovolny. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hana M. Dobrovolny, email@example.com