Evaluation of single drug treatment against SARS-CoV-2 in Vero E6 cells
Figure 1 and Table 1 show the anti-SARS-CoV-2 activities and cytotoxicity of the repurposed drugs in Vero E6 cells. The plaque assay was used to determine the viral production and is expressed as the percent inhibition relative to the viral titer of DMSO-treated cells. The one-step qRT-PCR was used to quantitate the viral RNA in virus supernatants and is also expressed as the percent inhibition relative to the DMSO-treated cells. The IC50 values calculated from the dose-response determined by plaque assay for Niclosamide, Ivermectin, and Chloroquine were 0.049, 1.23, 0.046 and 0.83 µM, respectively. The IC50 values calculated from the dose-response determined by one-step qRT-PCR for Niclosamide, Ivermectin, and Chloroquine were 0.043, 1.27, and 0.89 µM, respectively. Both methods used for viral quantification resulted in similar IC50 values. Thus the viral RNA quantification by the one-step qRT-PCR accurately determined the infectious virus output in these experiments, and could be used for the further two-drug combination experiments for the high throughput screening.
Table 1
Single drug treatment against SARS-CoV-2 in vitro.
Drug candidates | Drug class | Drug indication | CC50 (µM) | IC50 (µM) Plaque assay | IC50 (µM) qRT-PCR |
Niclosamide | Anthelminthic agents | Treatment of tapeworm and intestinal fluke infections [13] | 0.29 | 0.049 | 0.043 |
Ivermectin | Anti-parasitic agents | Treatment of onchocerciasis, and other worm infestations [14] | 10.55 | 1.23 | 1.27 |
Chloroquine | Anti-malarial agents | Treatment of malaria, rheumatic diseases and Zika virus infection [50] | 118.20 | 0.83 | 0.89 |
Evaluation of two-drug combination treatments against SARS-CoV-2 in Vero E6 cells
Firstly, the antiviral activities of two-drug combinations were assessed in vitro in Vero E6 cells. The cells were treated with 16 different pairwise combinations of two drugs, including, Niclosamide-Ivermectin, Niclosamide-Chloroquine and Ivermectin-chloroquine.
Niclosamide-ivermectin Combination
The presence of Ivermectin induced a shift in the dose-response curve of Niclosamide, with 10.75-fold reduction of Niclosamide IC50 value in the presence of 2.4 and 1.2 µM Ivermectin and approximately 2-fold reduction of Niclosamide IC50 value in the presence of 0.6 and 0.3 µM Ivermectin (Fig. 2A, Table 2). In a similar way, the presence of Niclosamide also induced a shift in Ivermectin dose-response curve, with 26.46-fold reduction of Ivermectin IC50 value in the presence of 0.09 µM Niclosamide (Fig. 2B, Table 2). The presence of 0.045 µM, 0.0225 µM and 0.01125 µM Niclosamide resulted in 7.18, 4.06 and 1.92-fold reduction of Ivermectin IC50 value, respectively. The dose-response matrix of Niclosamide and Ivermectin combination showed the obvious increasing inhibitory effects (Fig. 2C). A synergy landscape plot showed high positive Loewe synergy scores in combinations with Ivermectin concentration higher than 0.6 µM with a peak score of 22.76 indicating a synergistic effect. The scores were low positive to slightly negative in the other part of the plot with lower Ivermectin concentration indicating only additive effect at these low concentrations (Fig. 2D). The mean Loewe synergy score is 6.60. The ZIP, Bliss independence and HSA reference models were also used, the results showed the synergy scores of 12.64, 12.77 and 19.03, respectively, which accounted for the synergistic effect between Niclosamide and Ivermectin in Vero E6. No significant cytotoxicity in all 16 pairwise combinations (Fig. 2A, B).
Table 2
Antiviral activity of two-drug combinations treatment against SARS-CoV-2 in Vero E6 cells.
Drug treatment | IC50 (µM) qRT-PCR | Fold reduction of IC50 (single/combined) |
Niclosamide-Ivermectin | Niclosamide | 0.043 | |
Niclosamide + Ivermectin 2.4 µM | 0.004 | 10.750 |
Niclosamide + Ivermectin 1.2 µM | 0.004 | 10.750 |
Niclosamide + Ivermectin 0.6 µM | 0.018 | 2.399 |
Niclosamide + Ivermectin 0.3 µM | 0.022 | 1.955 |
Ivermectin | 1.27 | |
Ivermectin + Niclosamide 0.09 µM | 0.048 | 26.46 |
Ivermectin + Niclosamide 0.0045 µM | 0.177 | 7.18 |
Ivermectin + Niclosamide 0.0225 µM | 0.313 | 4.06 |
Ivermectin + Niclosamide 0.01125 µM | 0.660 | 1.92 |
Niclosamide-Chloroquine | Niclosamide | 0.043 | |
Niclosamide + Chloroquine 1.7 µM | 0.003 | 14.333 |
Niclosamide + Chloroquine 0.85 µM | 0.009 | 4.778 |
Niclosamide + Chloroquine 0.425 µM | 0.013 | 3.308 |
Niclosamide + Chloroquine 0.2125 µM | 0.029 | 1.483 |
Chloroquine | 0.89 | |
Chloroquine + Niclosamide 0.09 µM | 0.028 | 31.78 |
Chloroquine + Niclosamide 0.0045 µM | 0.193 | 4.61 |
Chloroquine + Niclosamide 0.0225 µM | 0.249 | 3.57 |
Chloroquine + Niclosamide 0.01125 µM | 0.531 | 1.68 |
Ivermectin-Chloroquine | Ivermectin | 1.27 | |
Ivermectin + Chloroquine 1.7 µM | 0.023 | 55.22 |
Ivermectin + Chloroquine 0.85 µM | 0.122 | 10.41 |
Ivermectin + Chloroquine 0.425 µM | 0.515 | 2.47 |
Ivermectin + Chloroquine 0.2125 µM | 0.821 | 1.55 |
Chloroquine | 0.89 | |
Chloroquine + Ivermectin 2.4 µM | 0.014 | 63.57 |
Chloroquine + Ivermectin 1.2 µM | 0.221 | 4.03 |
Chloroquine + Ivermectin 0.6 µM | 0.315 | 2.83 |
Chloroquine + Ivermectin 0.3 µM | 0.514 | 1.73 |
Niclosamide-chloroquine Combination
By using the same approach, it was found that the presence of Chloroquine induced a shift in Niclosamide dose-response curve, with 14.333, 4.778, 3.308 and 1.483-fold reduction of Niclosamide IC50 value in the presence of 1.7, 0.85, 0.425, and 0.2125 µM Chloroquine, respectively (Fig. 3A, Table 2). A similar trend was observed for the Chloroquine dose-response curve in the presence of Niclosamide, with 31.78, 4.61, 3.57 and 1.68-fold reduction of Chloroquine IC50 value in the presence of 0.09, 0.045, 0.0225 and 0.01125 µM Niclosamide, respectively (Fig. 3B, Table 2). The dose-response matrix shows increasing inhibitory effect of the combination with higher concentrations of Niclosamide and Chloroquine (Fig. 3C). The synergy map shows positive synergy scores at high concentrations of both drugs, while the lower concentrations gave zero and negative synergy scores with a peak positive score of 18.57, indicating a synergistic effect. (Fig. 3D). As most parts of the surface had Loewe synergy scores between -10 and 10, except for the highest concentrations of both drugs, with a mean score of 0.073, it suggests an additive effect between Niclosamide and Chloroquine. Additionally, the synergy scores calculated using ZIP and Bliss independence reference models gave the values of 3.86 and 3.67, respectively, which similarly indicated the additive effect. The HSA model resulted in the synergy score of 11.41, which accounted for the small level in synergistic effect. No significant cytotoxicity in all 16 pairwise combinations (Fig. 3A, B).
Ivermectin-chloroquine Combination
The results showed that the presence of Chloroquine induced a shift in Ivermectin dose-response curve, with 55.22, 10.41, 2.47, 1.55-fold reduction of Ivermectin IC50 value in the presence of 1.7, 0.85, 0.425, and 0.2125 µM Chloroquine, respectively (Fig. 4A, Table 2). Similarly, the presence of Ivermectin also induced a shift in Chloroquine dose-response curve, with 63.57, 4.03, 2.83 and 1.73-fold reduction of Chloroquine IC50 value in the presence of 2.4, 1.2, 0.6 and 0.3 µM Ivermectin, respectively (Fig. 4B, Table 2). The dose-response matrix shows increasing inhibitory effect with higher concentrations of Ivermectin and Chloroquine (Fig. 4C). Most parts of the synergy surface show negative synergy scores except for a small positive area with a peak positive score of 7.43 at the highest concentration of Chloroquine (Fig. 4D). The peak negative score of the surface is -8.16. As all of the surface had Loewe synergy scores between -10 and 10 with a mean score of -1.812, it suggests an additive effect between Ivermectin and Chloroquine. Moreover, both ZIP, Bliss independence and HSA reference models showed the synergy scores of 1.97, 1.98 and 9.63 which indicated the additive effect. No significant cytotoxicity in all 16 pairwise combinations (Fig. 4A, B).
Evaluation of single drug treatment against SARS-CoV-2 in Calu-3 cells
The best antiviral activity and calculated synergy scores demonstrated in the treatment with Niclosamide-Ivermectin combination in Vero E6 cells. Therefore, this two-drug combination was selected for the further evaluation in the human lung cancer cell line, Calu-3. The antiviral activities of single Niclosamide and Ivermectin treatments were assessed in Calu-3 cells (Fig. 5). The IC50 values of both drugs were 0.2 µM in Calu-3 cells. The CC50 values of Niclosamide and Ivermectin were 5.62 µM and 3.10 µM, respectively. The SI values of Niclosamide and Ivermectin were 28.1 and 15.5, respectively.
Evaluation of Niclosamide-Ivermectin combination treatment against SARS-CoV-2 in Calu-3 cells
The strong shifts were observed in the dose-response curves of Niclosamide combined with 0.4 and 0.2 µM Ivermectin (Fig. 6A, Table 3). The pairwise combinations of four different concentrations of Niclosamide with 0.4 and 0.2 µM Ivermectin resulted in a similar percent inhibition, thus, it was unable to calculate accurate IC50 values with the least curve fit. The presence of 0.1 and 0.05 µM Ivermectin also induced a shift in the dose-response curve of Niclosamide, in a similar level of 2.38 and 2.33-fold reduction of Niclosamide IC50 values, respectively (Fig. 6A, Table 3). In a similar way, the presence of Niclosamide induced a shift in Ivermectin dose-response curve with 8.69, 4.88, 3.64, and 2.41-fold reduction of Ivermectin IC50 value in the presence of 0.4, 0.2, 0.1 and 0.05 µM Niclosamide, respectively (Fig. 6B, Table 3). The dose–response matrix shows the increasing antiviral activity compared to the single drug treatments (Fig. 6C). The combination synergy analysis showed the mean Loewe synergy score of 2.897, which accounted for the additive effect between Niclosamide and Ivermectin in Calu-3 cells (Fig. 6D). Additionally, the peak Loewe synergy score was 13.19. The synergy score obtained from ZIP and Bliss independence reference models were 0.886 and 0.954, respectively, which also accounted for the additive effect. The synergy score calculated using HSA models was 10.795, which indicated a small synergistic effect between Niclosamide and Ivermectin. All 16 pairwise combinations showed no significant cytotoxicity (Fig. 6A, B)
Table 3
Evaluation of Niclosamide-Ivermectin combination treatments against SARS-CoV-2 in Calu-3 cells
Drug treatment | IC50 (µM) | Fold reduction of IC50 (single/combined) |
Favipiravir-Ivermectin | Favipiravir | 0.20 | |
Niclosamide + Ivermectin 0.4 µM | ND | ND |
Niclosamide + Ivermectin 0.2 µM | ND | ND |
Niclosamide + Ivermectin 0.1 µM | 0.084 | 2.38 |
Niclosamide + Ivermectin 0.05 µM | 0.086 | 2.33 |
Ivermectin | 0.20 | |
Ivermectin + Niclosamide 0.4 µM | 0.023 | 8.69 |
Ivermectin + Niclosamide 0.2 µM | 0.041 | 4.88 |
Ivermectin + Niclosamide 0.1 µM | 0.055 | 3.64 |
Ivermectin + Niclosamide 0.05 µM | 0.083 | 2.41 |
ND = not determined, cannot calculate IC50 with the least curve fit of the data sets. |