INTRODUCTION
Malaria is one of the leading global health threats, causing 228 million clinical episodes and 405,000 deaths in 2018 (
1). The continued development of drug resistance in
Plasmodium falciparum over the past decades has led to a global adoption of artemisinin-based combination therapies (ACTs) as a first-line treatment for uncomplicated malaria (
2). However, the emergence and spread of artemisinin resistance in
P. falciparum (
3–5) has evoked a need for the development of new antimalarials, especially alternative treatments in case of standard ACTs treatment failure, and new combinations of available drugs that can be immediately deployable and do not rapidly induce resistance (
6–9).
Another important consideration is that malaria incidence in travelers to areas of endemicity has continued to rise steadily (
10,
11). With continued growth in international travel, more than 120 million international travelers arrive annually in countries where malaria is endemic, resulting in an average of nearly 10,000 travel-related malaria cases per year (
12). Despite many decades of intense research and development effort, there is currently no effective malaria vaccine. Antimalarials remain the essential options for malaria prophylaxis, and nonuse or inappropriate use of chemoprophylaxis is a major risk factor of fatal malaria associated with travel (
13). Unfortunately, the use of available antimalarials for chemoprophylaxis is restricted by drug resistance, intolerable side effects or adverse effects, and noncompliance due to daily or complicated dosing regimens (
10,
11,
14). Thus, more effective, safer, and more convenient chemoprophylaxis drugs for malaria are urgently needed.
Azithromycin (AZ) is a widely used macrolide antibiotic with antimalarial activity (
15). It has been known to inhibit ribosomal protein synthesis in the
Plasmodium apicoplast by binding to the apicoplast ribosomal 50S subunit (
16), leading to failure in the production of infective daughter merozoites in the asexual stage, as well as the reduction of gametocyte-ookinete transformation and sporozoite production in the sexual stage (
17–20), possibly by interfering with protein prenylation, which results in defects in cellular trafficking (
21). Given its slow and relatively weak delayed death effect on asexual-stage parasites, AZ monotherapy is not useful for uncomplicated malaria and is only moderately effective for
P. falciparum prophylaxis (
22). However, it has been considered as a potential candidate for combination therapy, especially in pregnant women and children, owing to its multistage targeting action with unique mechanism, attractive safety profile, and excellent pharmacological properties (
20,
23). AZ combination therapy with other antimalarials has been evaluated extensively in clinical trials over the past dozen years, but no reliable evidence has been found on the superiority or equivalence of AZ combination therapy with current partner drugs compared to the current first-line antimalarial combinations (
22,
24–30). A more appropriate partner for AZ in combination remains to be identified.
Naphthoquine (NQ) is a 4-aminoquinoline antimalarial first synthesized in 1981 by Li et al. (
31). A single-dose, fixed coformulation of NQ with artemisinin (
32) has been marketed under the name ARCO for 15 years. Based on the excellent efficacy and safety data of NQ-containing therapies from available studies involving more than 4,000 patients, NQ has been proposed as an emerging candidate for new ACT or triple therapy to respond to the concern on resistance (
33). In addition, NQ has a long half-life up to 23 days in human (
33,
34), and its excellent efficacy for seasonal malaria chemoprophylaxis with monthly single dose has been reported in China (
35,
36). In a posttreatment follow-up study of artemisinin-NQ in Papua New Guinean children, the posttreatment prophylactic effect of NQ even was evident at day 63 or beyond (
37). Based on these advantages, NQ should be a highly promising candidate for further development as a chemoprophylactic antimalarial.
Recently, a coformulated combination of AZ with NQ has been evaluated in a clinical trial for malaria prophylaxis in Southeast Asia, and it shown a good safety profile and efficacy in adults with a monthly dose of 400 or 800 mg (
38). However, no study has yet investigated the interaction between AZ and NQ, and the possible benefits of combining NQ and AZ for malaria treatment or prophylaxis remains to be clarified. In the present study, we evaluated the combination of AZ and NQ in an animal malaria model. The pharmacodynamic interaction and its effect on the blood-stage prophylaxis and drug resistance development were investigated to determine whether the combination therapy with AZ and NQ can provide some benefits for malaria treatment or prophylaxis.
DISCUSSION
Combination therapy to increase efficacy, decrease toxicity, and minimize or slow down the development of drug resistance has become the standard choice for malaria treatment and is an important strategy for the development of new antimalarial drugs. However, accurate evaluation of the drug-drug interactions in a combination is complicated because of the influence of the drug ratios or doses used and potency desired. There have been many hypotheses, models, theories, and approaches, as well as some controversies, regarding drug combination analysis (
40–42). In the present study, the Chou-Talalay method was used because it has been well recognized as the most impactful approach to quantifying synergy and cited over 9,000 times in more than thousand biomedical journals over the past few decades (
40,
43,
44).
Our data provided evidence of synergy at higher inhibition levels but antagonism at lower inhibition levels between NQ and AZ in combinations at various ratios, and a decrease in the AZ proportion will lead to a decrease in the dose range of synergy from lower inhibition levels to higher inhibition levels (
Fig. 1B). A nonsignificant but visible enhancement of synergistic effect at the ED
90 level was also observed with the increase in AZ proportion (see Results and
Table 2). Apart from one report of synergistic interaction between AZ and NQ against
P. falciparum in vitro, in which only one ratio of AZ to NQ (34:40) was tested (
45), no other studies have been reported concerning the interaction between AZ and NQ. However, a previous study on combinations of AZ and CQ, another 4-aminoquinoline antimalarial drug very similar to NQ, reported a comparable finding in an
in vitro test against
P. falciparum. This study demonstrated synergistic effect at the FIC
90 level but additive interaction at the FIC
50 level, and the increase in AZ proportion led to greater synergy (
23). The interpretation proposed by the authors is that the lysosomotropic property of AZ might cause its accumulation in the digestive vacuole of the parasite at a high concentration, which perhaps can indirectly increase the potency of CQ by affecting pH and/or the process of hemoglobin digestion (
23). Although the drug-drug interaction
in vivo might be different from that
in vitro because the
in vivo pharmacokinetics cause more complex changes in drug exposures over time, the considerable concordance between our results and the findings in the previous
in vitro test of AZ and CQ indicated that, in addition to possible pharmacokinetics interaction
in vivo, the interaction between AZ and NQ in our study could also be partially interpreted as similar to the
in vitro interaction between AZ and CQ. However, the interpretation given above is not satisfactory because CQ and NQ are also lysosomotropic weak bases known to accumulate in the digestive vacuole and exert their antimalarial activity by binding to toxic heme released via hemoglobin proteolysis, thus interfering with parasite-mediated heme detoxification (
33,
46,
47). Although the precise mechanism of AZ and CQ or NQ uptake in
Plasmodium parasite remains unknown, it is believed that their accumulation in the digestive vacuole is largely dependent on the transmembrane pH gradient, and competition of intravacuolar accumulation between them will occur, as previously observed between AZ and ammonium chloride in rat alveolar macrophages (
48). Thus, AZ accumulation in digestive vacuoles should diminish the accumulation of NQ in same cellular compartment and vice versa. The former obviously means a decrease in NQ potency, while the later implies an increase in AZ potency, because it means the higher proportion of AZ molecules were retained in the parasite cytoplasm to target the apicoplast (
16,
17). Based on these considerations, we assume that an imbalance or balance between the accumulation of AZ and NQ in digestive vacuoles could be one reason for the different types of interaction that occurred in the combinations with different drug ratios and inhibition levels. When the increase in the total contribution of AZ to parasite inhibition caused by its decreased accumulation in digestive vacuoles exceeds the decrease in NQ contribution caused by the same reason, the interaction is a synergy; otherwise, it is addition or antagonism. The synergistic interaction at higher inhibition levels and the additive or even antagonistic interaction at lower inhibition levels observed in our study and a previous report about AZ and CQ indicated that the contribution of increased AZ potency tended to exceed that of decreased NQ or CQ potency at higher inhibition levels, such as at the ED
90 level and above.
Of note, although greater synergy at the ED
90 level and a wider dose range of synergy were observed with an increased AZ proportion, a further increase beyond a certain ratio does not seem to work better, as indicated by the only nearly additive interaction observed for two combinations with a AZ/NQ ratio f above 26:1 (
Table 1 and
Fig. 1C) and the similar inhibition level for the cut-point of CI = 1 that occurred in the Fa-CI trendlines of the combinations with AZ/NQ ratios of 20:3 and 10:3 (
Fig. 1B). Moreover, an increase in the AZ proportion with the consequent synergy enhancement did not result in a significant difference in the efficacy of three fixed-ratio combinations (see Results and
Table 3). This finding suggested that stronger synergy may not always lead to higher overall efficacy. One interpretation for this could be that at the dosage tested in our study, AZ predominantly exerted a delayed death effect by targeting apicoplast, which only causes growth arrest in daughter parasites and was much weaker and slower than the parasiticidal activity of NQ (
16,
17). Furthermore, a large decrease in the NQ proportion might cause a large overall potency loss in a combination enough to offset the gain from a synergistic interaction.
Since a high level of parasitemia inhibition is always desired for malaria treatment, the synergy observed at the ED90 level or above in a wide range of dose ratios of AZ to NQ should be much more therapeutically relevant than the addition or antagonism observed at lower inhibition levels, and this may provide clearly positive evidence for the value of combination therapy with AZ and NQ.
Considering the mismatch between the terminal half-life of AZ and NQ, we further investigated how long after oral administration the AZ and NQ remaining in the body can still interact synergistically by using the blood-stage prophylaxis models. In a pilot experiment of
P. berghei K173-challenged mouse models, the 1:1 combination provided the most effective and sustained prophylaxis (
Fig. 2). This suggested that the AZ and NQ doses at this ratio maintained the best balance between the contributions of each drug, and the interaction between them to provide the best efficacy lasted the longest
in vivo. In further tests performed with
P. berghei K173-challenged mice and
P. cynomolgi bastianelli-challenged monkeys, replacing half the amount of NQ with AZ did not significantly decrease effectiveness, whereas adding the same amount of AZ to NQ achieved a great improvement for 4 weeks after the first dose, although AZ alone presented only slight or no prophylactic effect and was much weaker than NQ (
Fig. 3 and
4). This suggested that when combined with NQ at an optimal ratio of 1:1, the weak or ineffective AZ could exert an effect against the growth of blood-stage parasite comparable to that of NQ, and this positive interaction could occur even at 4 weeks after the first dose, although AZ has a much shorter half-life than NQ (6.42 h versus 198.6 h in normal mice after a single oral dose) (
49,
50).
NQ or AZ resistance has not been documented in clinical setting, but previous laboratory studies have found that both NQ- and AZ-resistant lines could be rapidly selected, and an NQ-resistant
P. berghei line was highly cross-resistant to CQ (I
90 = 14.5,
P < 0.01) (
16,
51). Considering the widespread CQ resistance in areas where malaria is endemic, the highly CQ-resistant
P. berghei RCQ/K173 was used to determine the synergistic antimalarial efficacy of the drug combination against CQ-resistant parasite. For the 1:1 combination, a slightly increased synergistic effect against the CQ-resistant parasite was observed. Moreover, compared to the significantly decreased sensitivity to NQ, the decrease in sensitivity to combination in the CQ-resistant parasite was slight and not significant (I
90 = 1.9,
P = 0.146097) (
Table 3). The former suggested that the effect of the decrease in NQ potency caused by AZ on a CQ-resistant parasite might be smaller than that on a sensitive parasite, and the latter indicated that combination therapy should be more likely to avoid efficacy diminishment in CQ-resistant malaria treatment than NQ monotherapy. A parallel comparison of drug resistance development in
P. berghei K173 under selective pressure of drugs was performed to further determine the effect of the combination on the development of resistance. The delayed and slower development of resistance under combination selection pressure compared to that under AZ and NQ pressure alone clearly proved that the combination delayed the emergence of resistance and slowed down the development of resistance (
Fig. 5).
Taken together, these findings provide positive evidence of the benefits of combination therapy with AZ and NQ. Combination therapy using drugs with different mechanisms of action is the current state of the art in antimalarial treatment. However, apart from ACTs, there are relatively few other combinations available. In light of increasing concern regarding the emergence and spread of artemisinin resistance, the combination of AZ and NQ should be a viable candidate for alternative treatment development. In addition, the positive interaction between AZ and NQ after oral administration to provide prophylaxis lasting 4 weeks in mice and rhesus monkeys in this study, as well as the excellent prophylactic efficacy with monthly dose regimens reported in a previous clinical study (
38), suggest that this combination may also be an ideal candidate for long-term chemoprophylaxis of malaria.
MATERIALS AND METHODS
Drugs.
Naphthoquine phosphate and AZ were obtained from Shanghai Sixth Pharmaceutical Factory (Shanghai, China) and Shanghai Modern Pudong Pharmaceutical Co., Ltd. (Shanghai, China), respectively. Drug suspension for oral administration were freshly prepared on the day of drug administration by grinding the drug mixed with small amount of Tween 80 and subsequently suspending it in double-distilled water to the desired concentration. The dose of NQ was calculated as the base.
Parasites and experimental animals.
The P. berghei K173, P. berghei RCQ/K173, and P. cynomolgi bastianelli L parasites were maintained as cryopreserved stabilates or by mechanical blood passage.
Kunming mice (4 to 5 weeks old, 18 to 22 g) and rhesus macaques (Macaca mulatta, 2.7 to 3.7 kg) were obtained from the Animal Research Center of Academy of Military and Medicinal Science (Beijing, China). Animals were housed in standard laboratory cages and maintained under 12-h light-dark cycle with ad libitum access to food and water at a constant temperature (25 to 28°C) and humidity (65 to 80%). All animal experiments were carried out in accordance with institutional guidelines for animal care, using protocols approved by the Institutional Animal Care and Use Committee at Beijing Institute of Microbiology and Epidemiology. Animals were infected with parasites by intraperitoneal or intravenous injection of cryopreservation or blood from an infected donor animal. Parasitemia was monitored microscopically in Giemsa-stained blood smears prepared from mouse tail and monkey ear blood. Drugs were administered by oral gavage by using a feeding needle for mice and a stomach tube for monkeys.
Antimalarial activity assay.
Antimalarial activity against
P. berghei K173 infection in mice was assessed by Peters’ 4-day suppressive test (
39). Mice intraperitoneally inoculated with 1 × 10
7 parasite-infected erythrocytes (iRBC) were randomly grouped (
n = 10) and administered 0.2 ml of drug suspension once daily for 4 consecutive days or vehicle only for the nontreated control. On day 4 after infection, parasitemia was examined in each animal under a microscope at ×100 magnification (oil immersion). The parasitemia level was determined by counting the number of iRBC from more than 200 erythrocytes in random fields of the microscope. Average percentage of parasitemia inhibition for each test group was calculated as 100 × (
A −
B)/
A, where
A is the average percentage of parasitemia in the nontreated control group, and
B is the average percentage of parasitemia in the test groups. For the determination of effective dose, each drug was tested in five dose groups. The dose-response curve was fitted by a linear regression model based on the median-effect equation (
40) (equation 1), in which
D is the drug dose,
fa is the percentage of parasitemia inhibited by
D,
fu is the percentage of uninhibited parasitemia (
fu = 100 –
fa),
Dm is the median-effect dose that inhibits parasitemia by 50% (ED
50), and
m is the coefficient signifying the shape of the dose-effect relationship.
Drug combination analysis.
The Chou-Talalay method (
40) was used for combination analysis. Combination index (CI) values were used to determine an additive (CI = 1), antagonistic (CI > 1), or synergistic (CI < 1) interaction between drugs and were calculated using equation 2, in which
Daz and
Dnq, (ED
x)
az, and (ED
x)
nq are the doses of AZ and NQ used in combination or alone, respectively, to produce that same effect (i.e., % parasitemia inhibition).
The Fa-CI plot and normalized isobologram was created to present an effect-oriented or dose-oriented view on drug interactions. Data points above or below the line of additivity indicate antagonism or synergy, respectively. Curve-shift analysis (
41) was also performed to directly compare the dose-response curves and determine the drug interaction, in which the dose-response curves of single drugs or combinations were normalized by transforming the dose of AZ and NQ alone or in combinations to the equivalents of its own ED
90 value (ED
90 eq). A left shift in the dose-response curve indicates synergy.
Prophylactic effect evaluation in blood-stage parasites challenge models.
For study in P. berghei K173 challenge model, grouped mice (n = 10) were administered combinations or single drugs alone once daily for 3 consecutive days and then intraperitoneal challenged once with 1 × 107 iRBC at 1, 2, 3, or 4 weeks after the first dose. Parasitemia was checked ca. 7 to 10 days after challenge. Blood from cardiac puncture in mice confirmed to be parasite-free via microscopy in 100 random fields was inoculated into a healthy mouse, which was then examined for ca. 7 to 10 days to confirm the absence of infection. Prophylactic efficiency was expressed as (1 − PT/PC) × 100%, where PT and PC are the percentages of animals positive for parasitemia in treated and placebo groups, respectively.
For studies using the P. cynomolgi bastianelli challenge model, grouped monkeys (n = 2 to 3) were administered combinations or single drugs alone once daily for 2 consecutive days and intravenously challenged once with 1 × 106 iRBC at 3, 4, or 5 weeks after the first dose. Parasitemia was checked at 5 and 10 days after challenge, followed by every 2 to 3 days until 30 days after challenge. Prophylaxis was deemed successful if the animal remained parasite-free during the 30-day follow-up monitoring period.
Animals administered vehicle only as placebo controls were included for in experiments to validate the challenge.
Selection of drug-resistant P. berghei.
Selection of drug-resistant parasites was performed as previously described (
50). Briefly, five
P. berghei K173-infected mice were treated with the combination or with single drugs alone at the ED
90 dose 3 to 5 days after infection. On day 7 after infection, parasitemia was checked in each mouse, and infected blood from mice with the highest parasitemia (>2%) was used to inoculate the next passage. The parasites were exposed to a gradually increased dose of drugs in subsequent passages. The level of resistance was evaluated after each five passages by calculating the I
90, which is defined as the ratio of the ED
90 of the resistant line to that of the parent line. The degree of resistance was categorized into four levels by the I
90 values as follows: (i) I
90 ≤ 1, susceptible; (ii) 1< I
90 ≤ 10.0, slightly resistant; (iii) 10 < I
90 ≤ 100.0, moderately resistant; and (iv) I
90 > 100.0, highly resistant.
Statistical analysis.
Statistical analysis and data visualization were conducted using R (
52). Dose-response curves were fitted as described above, and the goodness of fit was determined using
r2. Averages and standard deviations of ED
90 values in a fixed-ratio combination test were calculated from three independent experiments. Significant differences were determined using one-way analysis of variance (ANOVA), followed by Tukey’s honestly significant difference (HSD)
post hoc test. Log-linear or logistic regression analysis was performed on the data from the prophylaxis studies to determine the statistically significant relationship between treatments, challenge times after treatment, and the frequency of parasitemia arising in animals, and a likelihood ratio test was performed for goodness-of-fit and comparison analyses (
53). The R package vcd was used to produce mosaic plots with residual-based shadings for log-linear analysis (
54,
55). The effect of drugs on resistance development was determined by exponential regression analysis using I
90 as a response variable against passage and drugs, in which the significance of the total effects or of each effect of passages, drugs, and interactions between them on I
90 was tested by the F statistic, and the significance of coefficients in the fitted model was checked by the t statistic.