Prerpints.org logo
Preprint
Article

In Vitro Combinatorial Activity of Direct Acting Antivirals and Monoclonal Antibodies against the Ancestral B.1 and BQ.1.1 SARS-CoV-2 Viral Variants

Altmetrics

Downloads

97

Views

32

Comments

0

A peer-reviewed article of this preprint also exists.

supplementary.pdf (592.33KB )

Submitted:

29 December 2023

Posted:

03 January 2024

You are already at the latest version

Alerts
Abstract
Combination antiviral therapy may be helpful in the treatment of SARS-CoV-2 infection, however no clinical trial data are available and combined use of direct acting antivirals (DAA) and monoclonal antibodies (mAb) has been reported only anecdotally. To assess the cooperative effects of dual drug combinations in vitro, we used a VERO E6 cell based in vitro system with the ancestral B.1 or the highly divergent BQ.1.1 virus to test pairwise combinations of the licensed DAA, including nirmatrelvir (NRM), remdesivir (RDV) and the active metabolite of molnupiravir (EIDD-1931) as well the combination of RDV with four licensed mAbs (sotrovimab, bebtelovimab, cilgavimab, tixagevimab; tested only with the susceptible B.1 virus). According to SynergyFinder 3.0 summary and weighted scores, all the combinations had an additive effect. Within DAA/DAA combinations, paired scores with the B.1 and BQ.1.1 variants were comparable. In the post-hoc analysis weighting synergy by concentrations, several cases of highly synergistic scores were detected at specific drug concentrations, both for DAA/DAA and for RDV/mAb combinations. This was supported by in vitro confirmation experiments showing a more than linear shift of a drug effective concentration (IC50) at increasing concentrations of the companion drug, although the effect was prominent with DAA/DAA combinations but minimal or null with RDV/mAb combinations. These results support a cooperative effects of dual drug combinations in vitro, which should be further investigated in animal models before introduction into the clinic.
Keywords: 
Subject: Biology and Life Sciences  -   Virology

1. Introduction

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), which causes coronavirus disease 2019 (COVID-19), has expanded worldwide since its emergence in China in December 2019, resulting in a global pandemic with more than 700 million confirmed cases and about 7 million deaths documented as of September 2023 (https://covid19.who.int/). Massive vaccination and natural SARS-CoV-2 infections and reinfections have progressively expanded immunity to SARS-CoV-2, drastically reducing the incidence of COVID-19 and its severe outcomes in the immunocompetent population. As a result, the end to COVID-19 as a global health emergency was declared in May 2023 by the World Health Organization (https://www.who.int/europe/emergencies/situations/covid-19). However, SARS-CoV-2 continues to circulate endemically in the global population, causing serious illness in the most fragile categories such as elderly people, comorbid and immunocompromised patients (Semenzato et al., 2022; Turtle et al., 2023). The persistent circulation of SARS-CoV-2 in the human population has been driven by the ability to evade the host immune response, particularly since the emergence of the Omicron variant and the plethora of its sublineages, which continue to evolve rapidly and challenge population immunity (Qu et al., 2023). In addition, waning immunity combines with decreased vaccine acceptance to pave the way to novel epidemic waves (Menegale et al., 2023). Thus, early treatment of symptomatic patients with available antiviral medications remains crucial to protect fragile individuals and limit virus spread and evolution.
Currently available anti-SARS-CoV-2 agents include the direct acting antivirals (DAA) remdesivir (RDV), molnupiravir (MNP) and ritonavir boosted nirmatrelvir (NRM) as well as anti-spike monoclonal antibodies (mAb) (https://www.ema.europa.eu/en/human-regulatory/overview/public-health-threats/coronavirus-disease-covid-19/covid-19-medicines#authorised-covid-19-treatments-section). Unfortunately, continued virus evolution has progressively decreased or abolished mAb neutralization activity, with few cases of restored activity against the currently dominating virus lineages (Touret et al., 2023). On the other hand, DAA have maintained their activity across multiple variants, due to the high conservation of the targeted viral enzymes (Cho et al., 2023). However, current DAA may be unable to fully prevent severe disease progression, and their use is limited by a number of drug-drug interactions (NRM), inconvenient administration (RDV) or mutagenicity concerns (MNP) (Del Borgo et al., 2023). Indeed, the development of novel drugs or different therapeutic approaches remains a key priority (Li et al., 2023). While waiting for novel agents, one option recently gaining interest is to combine available drugs into a drug cocktail, based on the long-standing success of combination therapy in HIV and HCV infection (Shyr et al., 2021). The key advantage of this approach is that each drug used in the combination can act at a different step in the virus replication cycle, increasing the effectiveness of the treatment while reducing the emergence of drug-resistant viral variants. Interestingly, if the drugs interact synergistically, the active dose can be reduced with respect to monotherapy, potentially lowering toxicity. Although no combination therapy is approved for the treatment of COVID-19, some studies have recently evaluated the combined effect of licensed or candidate anti-SARS-CoV-2 compounds in vitro (Gidari et al., 2022; Nguyenla et al., 2022; Wagoner et al., 2022; Wang et al., 2022) and in animal models (Abdelnabi et al., 2022; Jeong et al., 2022) and anecdotal cases of combined administration of two anti-SARS-CoV-2 treatments in clinical practice have been reported (Gidari et al., 2023; Lanzafame et al., 2023; Marangoni et al., 2023; Mikulska et al., 2023).
In this work, we evaluated the in vitro interactions between approved anti-SARS-CoV-2 compounds in a live virus cell-based assay using the prototype SARS-CoV-2 B.1 virus and a more recent omicron lineage (BQ.1.1). We included all the three DAA (RDV, MNP, NRM) plus four mAb with activity against the prototype virus and some of the later variants.

2. Materials and methods

2.1. Cells and viral stocks

The VERO E6 monkey cell line (ATCC® CRL-1586) was used to propagate and titrate the viral stocks and to determine the antiviral activity of the antiviral compounds, alone or in combination. VERO E6 cells were propagated in DMEM High Glucose medium (Euroclone) supplemented with 10% Fetal Bovine Serum (FBS) and 1% Streptomycin/Penicillin (PS) (Euroclone) in a humidified incubator at 37°C with 5% CO2. The same medium with a lower FBS concentration (1%) was used for viral infection experiments. For all DAA activity experiments, VERO E6 cells were pre-treated with 0.5 µM CP-100356 hydrochloride inhibitor to reduce the efflux activity of P-glycoprotein, which is overexpressed by this cell line (Fiaschi et al., 2022; Zhu et al., 2022).
The SARS-CoV-2 stocks used for the experiments included the wild type B.1 (GISAID code EPI_ISL_2472896) and the Omicron BQ.1.1 (EPI_ISL_17257516) variant, kindly provided by the Department of Biomedical and Clinical Sciences Luigi Sacco, University of Milan, Italy. Viral stocks were expanded in VERO E6 cells and titrated as previously described (Vicenti et al., 2021).

2.2. Drugs and cytotoxicity assay

The P-glycoprotein inhibitor CP-100356 (cat. HY-108347) and the DAA, including RDV (cat. HY-104077), NRM (cat. HY-138687) and EIDD-1931 (the active metabolite of MNP; cat. HY-125033), were supplied by MedChemExpress (Monmouth Junction, NJ, USA) as powder and dissolved in 100% dimethyl sulfoxide (DMSO). Biosimilar mAb including bebtelovimab (BEB; ref. PXTA1750-100), sotrovimab (SOT; ref. PXTA1637-100), tixagevimab (TIX; ref. PXTA1032-100) and cilgavimab (CIL; ref. PXTA1033-100) were supplied by ProteoGenix SAS (Schiltigheim, France) as powder and dissolved in phosphate buffered saline (PBS).
The cytotoxicity of all investigated compounds was assayed in VERO E6 cells as previously described [Fiaschi 2022]. Briefly, cell viability was calculated using the CellTiter Glo 2.0 kit (Promega) and luminescence values were measured through the GloMax® Discover Multimode Microplate Reader (Promega) and elaborated with the GraphPad PRISM software version 9.0 (La Jolla). The half-maximal cytotoxic concentration (CC50) and the compound concentration allowing 90% cell viability (CC10) were calculated using a non-linear regression analysis of the dose-response curves and the ECanything GraphPad function. The CC10 was used for each compound as the starting concentration in the antiviral activity assay.

2.3. Antiviral activity assay

To determine the antiviral activity of each compound against the B.1 and BQ.1.1 SARS-CoV-2 variants, a direct yield reduction assay, based on the infection of cells in the presence of serial drug dilutions, was performed as previously described (Fiaschi et al., 2022). Briefly, semiconfluent pre-seeded VERO E6 cells were treated with 4-fold decreasing concentrations of each antiviral and, after 1 h incubation at 37°C with 5% CO2, the cells were infected at 0.005 multiplicity of infection (MOI). To test mAb activity, the viral stocks were pre-incubated at MOI 0.005 with mAb dilutions for 1h at 37°C and then the virus/mAb mixture was transferred onto semiconfluent pre-seeded VERO E6 cells. In both cases, after additional 72 h incubation, cell viability was measured with CellTiter Glo 2.0. The half-maximal inhibitory concentration (IC50) was calculated through a non-linear regression analysis of the dose-response curves generated with GraphPad PRISM. Infected and uninfected cells without drugs were used to calculate the 0% and 100% cell viability, respectively, and used to normalize the results.
The different DAA combinations (NRM/EIDD-1931, NRM/RDV, RDV/EIDD-1931) were evaluated against the B.1 and BQ.1.1 variants. In addition, RDV only was tested in combination with the four mAb (RDV/SOT, RDV/BEB, RDV/CIL, RDV/TIX) based on convenience in the perspectives of clinical use since both RDV and mAb must be administered parenterally. These combinations could only be tested against the B.1 prototype virus because none of the mAb considered was active against the BQ.1.1 variant. The protocol previously described to determine the antiviral activity of individual compounds was adapted to the dual compound combinations. To generate the dual combinations, six 4-fold serial dilutions of the first compound were mixed with six 4-fold serial dilutions of the second compound in 96-well plates to obtain 36 different combos in a two-dimensional matrix. Briefly, for the DAA combos, the drug combination matrix was incubated for 1h at 37°C on pre-seeded semiconfluent VERO E6 cells and the cultures were then infected with B.1 and BQ.1.1 viral stock at 0.005 MOI. For the mAb/antiviral combos, the B.1 viral stock was pre-incubated for 1h at 37°C with mAb dilutions and then the virus/mAb mixture was transferred at 0.005 MOI onto VERO-E6 cells pretreated with serial DAA dilutions. For both conditions, infected cells were incubated for additional 72 h at 37°C with 5% CO2 and cell viability was calculated with the Cell Titer-Glo 2.0 protocol as described above. Each plate also included the dose-response curves for each individual compound, a mock infection control and an infection control without antivirals to normalize the results. Two independent experiments were performed for each combination.

2.4. Dose-effect relationships with compound combinations

Antiviral activity data derived from compound combinations were analyzed by the online SynergyFinder 3.0 tool (https://synergyfinder.fimm.fi/). This version of the software was recently implemented with post-analysis options enabling to better explore and interpret the dose-effect relationships between antiviral agents at any tested concentration (Ianevski et al., 2022). The default Zero interaction potency (ZIP) model was used to determine the type of drug-drug interaction by comparing the change in the potency of the dose-response curves between individual drugs and their combinations. The ZIP model combines the advantages of both the Loewe additivity and the Bliss independence models, aiming at a systematic assessment of various types of drug interactions patterns that may arise in a high-throughput drug combination screening (Yadav et al., 2015). For each dual compound combination, the results were expressed as a single summary Synergy Score (SS), which is averaged over all the dose combination measurements and is interpreted as the percent change in response due to drug interactions with respect to the effect of the individual drugs (e.g., a synergy score of 20 indicates 20% increased activity beyond expectation). The interaction is visualized by 2D and 3D synergy maps that highlight the synergistic and antagonistic dose areas in red and green, respectively. As per the indications of the system, an antagonistic, additive, and synergistic effect of the drug combination is defined by a summary SS value below or equal to -10, from -10 to 10, and above 10, respectively (Ianevski et al., 2022). In addition, SynergyFinder 3.0 implements a “Weighting synergy by concentrations” post-hoc analysis returning a weighted SS metric as a measure of synergy at lower drug dose windows where toxicity is less likely and clinical application is favored (Ianevski et al., 2022).
To confirm the results obtained by SynergyFinder 3.0, we selected drug concentrations that were repeatedly detected as synergistic (SS >10) in the drug-drug matrix and directly measured the IC50 shift of the drug combination in experiments where Compound 1 at three fixed concentrations (IC50 plus two 2-fold dilutions) was mixed with six 4-fold decreasing concentrations of Compound 2. The three Compound 2 IC50 values obtained in the presence of the fixed concentrations of Compound 1 were compared with the IC50 of Compound 2 alone and three corresponding fold shift values in Compound 2 IC50 were calculated as: IC50 [Compound 2 alone] / IC50 [Compound 1 + Compound 2].

2.5. Statistical analysis

SynergyFinder 3.0 does not perform any statistical comparison among summary or weighted SS values obtained with different compound combinations. Since the summary SS is the average of all the individual SS obtained with all the concentration pairs tested, the overall combinatorial effects of the three DAA pairs as well as those of the three RDV/mAb groups were compared by the Kruskal-Wallis test followed by Mann-Whitney pairwise comparisons between groups. SS values obtained with the same DAA/DAA concentrations against the B.1 and BQ.1.1 viral variants were compared by the Wilcoxon signed rank test.

3. Results

The IC50 of EIDD-1931, RDV and NRM, when tested alone, were 2.40±0.40/1.59±0.44, 0.06±0.03/0.03±0.01 and 0.10±0.03/0.10±0.01 μM against B.1/BQ.1.1, respectively. The mAb tested, namely CIL, TIX, BEB and SOT were active only against the wild type B.1 SARS-CoV-2 variant (IC50 values: 0.20±0.13, 0.07±0.04, 0.03±0.01, 0.81±0.32 µg/ml, respectively). All the IC50, CC50 and CC90 values are reported in Supplementary Table 1.
The DAA combos tested against the wild type B.1 and BQ.1.1 strains were NRM/EIDD-1931, NRM/RDV, RDV/EIDD-1931 while all mAb were tested only against the wild type B.1 strain in combination with RDV, the only injectable DAA (RDV/SOT, RDV/BEB, RDV/CIL, RDV/TIX). To evaluate DAA/DAA or RDV/mAb interactions, we analyzed the data generated from combination experiments using the SynergyFinder 3.0 software and applying the ZIP multiple synergy algorithm that combines Bliss and Loewe reference models. An example of the output generated by analyzing the 36-cell matrix obtained with the NRM/RDV combo against the B.1 SARS-CoV-2 strain is shown in Figure 1. The weighted SS for all the DAA combinations indicated additivity against B.1 (2.1±0.4 for NRM/EIDD-1931, 1.9±1.4 for NRM/RDV and 0.9±0.7 for RDV/EIDD-1931) and BQ.1.1 (0.0±0.7 for NRM/EIDD-1931, 2.1±2.1 for NRM/RDV and 1.7±2.5 for RDV/EIDD-1931). Additivity was also observed for all RDV/mAb combos against B.1 (6.1±2.9 for RDV/SOT, 4.2±3.3 for RDV/BEB, 0.8±0.1 for RDV/CIL and 2.6±2.7 for RDV/TIX). The more comprehensive summary score metric confirmed additivity for all the combinations tested. Summary and weighted SS are indicated in Table 1, while the 2-D Synergy plots generated by each experiment are shown in Supplementary Figure 1.
Despite the overall additive effect shown by summary and weighted SS for all the combinations, there were specific drug combo concentrations, below to the IC90 of the individual drugs, showing synergy. Table 2 shows these cases, stratified by their SS (10-29.9, 30-50, and above 50). The highest SS were obtained for several RDV/mAb concentrations, particularly with SOT, BEB and TIX, which was also reflected by the highest overall summary and weighted SS (Table 1). Among DAA combinations, only NRM/RDV reached the highest SS stratum but only at one concentration while hits at lower synergistic strata were detected for NRM/EIDD-1931, RDV/EIDD-1931 and again NRM/RDV. The number of hits showing synergy for NRM/EIDD-1931, NRM/RDV and RDV/EIDD-1931 was 5, 6 and 2, respectively. When comparing SS values obtained with all the DAA/DAA and DAA/mAb combinations and the B.1 virus, the Kruskal-Wallis test revealed significant differences (p < 0.001) and post-hoc pairwise analysis detected significantly higher synergy for RDV/SOT compared with the other RDV/mAb combinations (RDV/BEB, p = 0.006; RDV/CIL, p < 0.001; RDV/TIX, p = 0.006) and with two DAA/DAA combinations (NRM/EIDD-1931, p < 0.001; NRM/RDV, p = 0.010; RDV/EIDD-1931, p = 0.004). However, SS values indicating synergistic effects often were outliers in the overall distribution of data points for most dual combinations (Supplementary Figure 2). SS values remained significantly higher for the RDV/SOT combo when comparing only the four DAA/mAb combinations (p < 0.001 compared with RDV/CIL and p = 0.002 compared with RDV/BEB and RDV/TIX). Within the DAA/DAA combination subset a comparable number of cases showing synergy was detected with the B.1 and BQ.1.1 viruses (8 each) and there was no statistically significant difference in the SS values when comparing B.1 to BQ.1.1 at paired DAA/DAA concentrations (p = 0.153).
As a proof of concept, the IC50 synergistic potency shift was measured in infected cells treated with three fixed drug concentrations of the first compound plus scalar concentrations of the second compound. As fixed drug, we selected NRM and RDV because they repeatedly showed synergy at specific concentrations (0.05 and 0.02 μM, respectively) close to their IC50 (Table 2). As shown in Figure 2 and in Supplementary Table 2, the IC50 of NRM with the addition of RDV 0.06 μM and 0.03 μM was reduced by 33- and 8-fold against B.1 and by >140- and 14-fold against BQ.1.1, respectively. A synergistic potency shift in EIDD-1931 IC50 was also induced by the addition of NRM at 0.1 μM and 0.05 μM against B.1 (88- and 11-fold, respectively) and against BQ.1.1 (28- and 7-fold, respectively) and by the addition of RDV at 0.06 μM and 0.03 μM against B.1 (>26- and 26- fold, respectively) and against BQ.1.1 (>30- and 30-fold, respectively). As shown in Figure 2 and in Supplementary Table 2, lower synergistic potency shifts were observed when fixed RDV concentrations were tested in combination with mAb. Indeed, the addition of RDV 0.06 μM reduced SOT, BEB and TIX IC50 by 4-, 2- and 3- fold respectively, while the addition of RDV 0.03 μM reduced SOT and TIX IC50 by 2-fold and did not affect BEB IC50. At the lowest concentration tested (0.015 μM), RDV induced a 2-fold reduction only in SOT IC50. In conclusion, the results obtained in the proof-of-concept experiments suggested a synergistic interaction with the DAA/DAA combinations, however a cooperative effect between RDV and mAb was measurable only for RDV/SOT.

4. Discussion

While antiviral combination therapy has met with extraordinary success in the treatment of chronic infections, including HIV and HCV, expanding the concept to acute infections such as influenza and COVID-19 has just started to be considered (Koszalka et al., 2022; Yan and Yan, 2023). A strong argument for combining drugs with different mechanisms of actions is the need to limit the emergence of drug resistance which is clearly a major issue with persistent viruses. However, there have been reports of rapid selection of drug resistance also with exposure of SARS-CoV-2 to mAb or DAA, particularly in immunocompromised subjects requiring prolonged treatment (Ahmadi et al., 2023; Hogan et al., 2023; Yang et al., 2023; Zuckerman et al., 2023). Combination therapy with SARS-CoV-2 has been so far limited to mAb cocktails, including casirivimab/imdevimab, bamlanivimab/etesevimab and CIL/TIX. However, these combinations target the same virus replication step.
Although combined administration of two anti-SARS-CoV-2 treatments in the clinical setting has been anecdotally reported (Gidari et al., 2023; Lanzafame et al., 2023; Marangoni et al., 2023; Mikulska et al., 2023), it is currently unlikely that clinical trials comparing drug combinations to single-drug anti-SARS-CoV-2 therapy are designed. Thus, in vitro assessment of the potential for cooperative effects with multiple drugs plays a key role at this time. Of note, privileging in vitro over difficult-to-obtain in vivo data is not novel in the SARS-CoV-2 treatment arena, where authorization of mAb has been often revoked based on in vitro inefficacy data, rather than on clinical trials, with newly emerging viral variants (Wu et al., 2022). Assessment of the dose-effect relationship with drug combinations involves measuring inhibition of viral replication in a standardized virus-cell system where different concentrations of the drugs of interest are combined. Synergistic interactions are formally defined by activity scores which are higher than the sum of the activities of the individual drugs. Several tools are available online which are based on different mathematical models to define combinatorial effects as synergy, additivity, or antagonism (e.g., MacSynergy II™, CompuSyn, Combenefit, SynergyFinder)(Foucquier and Guedj, 2015). The available tools range from simple spreadsheet-based systems to software packages assisted by a graphical user interface and based on machine learning algorithms (Fuentealba-Manosalva et al., 2023). Since there is not a consensus on how to define synergy (Wooten et al., 2021), it is important to include post-hoc analysis and discover the specific drug concentrations resulting in unequivocal synergy. Indeed, by using the latest release of the popular SynergyFinder, we detected highly scored synergistic potency shifts at specific concentrations of the drug combinations tested although summary weighted SS indicated additive effects only and, notably, we could largely support this indication by confirmatory laboratory tests.
Our results can be compared with few published data of this kind. The NRM/EIDD-1931 combo was scored as synergistic in two previous studies, one using SynergyFinder 2.0 with the HSA algorithm to process data from a VERO E6 cell line-based system with the B.1 virus (Gidari et al., 2022), the other using SynergyFinder 3.0 with the Bliss algorithm, Calu-3 cells and the Delta virus (Wagoner et al., 2022). In addition, NRM and MNP were shown to increase the survival rate by >2-fold in the K18-hACE2 transgenic mouse model compared with either drug alone (Jeong et al., 2022). Another study reported an additive effect of the RDV/EIDD-1931 combination by applying the original SynergyFinder version with the ZIP algorithm in A549-ACE2-TMPRSS2 cells infected with the Omicron BA.5 virus (Abdelnabi et al., 2022). Thus, previous reports and our study consistently show a cooperative interaction with combinations of licensed DAA. The difference between additive and synergistic effects is not straightforward to define, due to the diversity of methods employed and the absence of a gold standard assay. We used the most updated version of SynergyFinder which, in addition to the summary SS metric, returns also a weighted SS computed at lower drug concentrations and allows post-hoc analysis at specific drug concentrations. Such analysis indeed revealed highly synergistic interactions which may have driven an overall synergistic result with earlier SynergyFinder versions (Gidari et al., 2022). Synergistic hits were scored for all the three DAA/DAA combination and for RDV with three of the four mAb considered. Indeed, our work is the first to comprehensively test all the three licensed DAA and perform a first-time analysis of RDV plus mAb. Other strengths of our experimental design include the use of the CP-100356 hydrochloride drug efflux inhibitor to better mimic in vivo DAA activity and the assessment of the DAA combinatorial effects against two highly divergent viruses such as the ancestral B.1 and the later emerging Omicron BQ.1.1 variant. Comparable effects were showed with the two viral strains, as expected from the known variant-independent activity of SARS-CoV-2 polymerase and protease inhibitors (Cho et al., 2023).
Finally, we corroborated SynergyFinder data with in vitro experiments to measure the IC50 shift for representative drug combinations where three 2-fold dilutions of one compound, downward from its IC50, were used together with a complete dilution series of a second compound. By selecting cases with repeated synergistic scores at concentrations close to the IC50, we showed that a 2-fold increase in the concentration of the first compound decreased the IC50 of the second compound by much more than 2-fold, supporting the synergistic effect indicated by SynergyFinder. Of note, while this was evident with DAA combinations, the effect was smaller or null with RDV/mAb combinations advising for further studies to better define the interaction between these drugs acting at extracellular and intracellular levels. This highlights the need to support data obtained from SynergyFinder or similar tools with in vitro confirmation experiments, particularly in light of the multiplicity and frequent update of the algorithms. Finally, in vivo experiments should be designed in a suitable animal model before moving forward to a wide introduction of combination therapy into the clinic.

Funding

This work received funding from the Ministero dell’Università e della Ricerca (PRIN 202022GZEHE_01); EU funding within the NextGenerationEU-MUR PNRR Extended Partnership initiative on Emerging Infectious Diseases (Project no. PE00000007, INF-ACT) and by the EuCARE Project funded by the European Union’s Horizon Europe Research and Innovation Program, Grant Agreement No. 101046016.

Declaration of competing interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Abdelnabi, R.; Maes, P.; de Jonghe, S.; Weynand, B.; Neyts, J. Combination of the parent analogue of remdesivir (GS-441524) and molnupiravir results in a markedly potent antiviral effect in SARS-CoV-2 infected Syrian hamsters. Front. Pharmacol. 2022, 13, 1072202. [Google Scholar] [CrossRef] [PubMed]
  2. Ahmadi, A.S.; Zadheidar, S.; Sadeghi, K.; Nejati, A.; Salimi, V.; Hajiabdolbaghi, M.; Mokhtari-Azad, T.; Yavarian, J. SARS-CoV-2 intrahost evolution in immunocompromised patients in comparison with immunocompetent populations after treatment. J. Med Virol. 2023, 95, e28877. [Google Scholar] [CrossRef]
  3. Cho, J.; Shin, Y.; Yang, J.-S.; Kim, J.W.; Kim, K.-C.; Lee, J.-Y. Evaluation of antiviral drugs against newly emerged SARS-CoV-2 Omicron subvariants. Antivir. Res. 2023, 214, 105609. [Google Scholar] [CrossRef] [PubMed]
  4. Del Borgo, C.; Garattini, S.; Bortignon, C.; Carraro, A.; Di Trento, D.; Gasperin, A.; Grimaldi, A.; De Maria, S.G.; Corazza, S.; Tieghi, T.; et al. Effectiveness, Tolerability and Prescribing Choice of Antiviral Molecules Molnupiravir, Remdesivir and Nirmatrelvir/r: A Real-World Comparison in the First Ten Months of Use. Viruses 2023, 15, 1025. [Google Scholar] [CrossRef] [PubMed]
  5. Fiaschi, L.; Dragoni, F.; Schiaroli, E.; Bergna, A.; Rossetti, B.; Giammarino, F.; Biba, C.; Gidari, A.; Lai, A.; Nencioni, C.; et al. Efficacy of Licensed Monoclonal Antibodies and Antiviral Agents against the SARS-CoV-2 Omicron Sublineages BA.1 and BA.2. Viruses 2022, 14, 1374. [Google Scholar] [CrossRef] [PubMed]
  6. Foucquier, J.; Guedj, M. Analysis of drug combinations: current methodological landscape. Pharmacol. Res. Perspect. 2015, 3, e00149. [Google Scholar] [CrossRef] [PubMed]
  7. Fuentealba-Manosalva, O.; Mansilla, M.; Buelvas, N.; Martin-Martin, A.; Torres, C.G.; López-Muñoz, R.A. Mind the Curve: Dose–Response Fitting Biases the Synergy Scores across Software Used for Chemotherapy Combination Studies. Int. J. Mol. Sci. 2023, 24, 9705. [Google Scholar] [CrossRef]
  8. Gidari, A.; Sabbatini, S.; Schiaroli, E.; Bastianelli, S.; Pierucci, S.; Busti, C.; Comez, L.; Libera, V.; Macchiarulo, A.; Paciaroni, A.; et al. The Combination of Molnupiravir with Nirmatrelvir or GC376 Has a Synergic Role in the Inhibition of SARS-CoV-2 Replication In Vitro. Microorganisms 2022, 10, 1475. [Google Scholar] [CrossRef] [PubMed]
  9. Gidari, A.; Sabbatini, S.; Schiaroli, E.; Bastianelli, S.; Pierucci, S.; Busti, C.; Saraca, L.M.; Capogrossi, L.; Pasticci, M.B.; Francisci, D. Synergistic Activity of Remdesivir–Nirmatrelvir Combination on a SARS-CoV-2 In Vitro Model and a Case Report. Viruses 2023, 15, 1577. [Google Scholar] [CrossRef]
  10. I Hogan, J.; Duerr, R.; Dimartino, D.; Marier, C.; E Hochman, S.; Mehta, S.; Wang, G.; Heguy, A. Remdesivir Resistance in Transplant Recipients With Persistent Coronavirus Disease 2019. Clin. Infect. Dis. 2022, 76, 342–345. [Google Scholar] [CrossRef]
  11. Ianevski, A.; Giri, A.K.; Aittokallio, T. SynergyFinder 3.0: an interactive analysis and consensus interpretation of multi-drug synergies across multiple samples. Nucleic Acids Res. 2022, 50, W739–W743. [Google Scholar] [CrossRef] [PubMed]
  12. Jeong, J.H.; Chokkakula, S.; Min, S.C.; Kim, B.K.; Choi, W.-S.; Oh, S.; Yun, Y.S.; Kang, D.H.; Lee, O.-J.; Kim, E.-G.; et al. Combination therapy with nirmatrelvir and molnupiravir improves the survival of SARS-CoV-2 infected mice. Antivir. Res. 2022, 208, 105430. [Google Scholar] [CrossRef] [PubMed]
  13. Koszalka, P.; Subbarao, K.; Baz, M. Preclinical and clinical developments for combination treatment of influenza. PLOS Pathog. 2022, 18, e1010481. [Google Scholar] [CrossRef] [PubMed]
  14. Lanzafame, M.; Gottardi, M.; Guella, L.; Collini, L.; Costa, G.; Guella, A.; Vento, S. Successful treatment of persistent SARS-CoV-2 infection with nirmatrelvir/ritonavir plus sotrovimab in four immunocompromised patients. J. Chemother. 2023, 35, 623–626. [Google Scholar] [CrossRef]
  15. Li, G.; Hilgenfeld, R.; Whitley, R.; De Clercq, E. Therapeutic strategies for COVID-19: progress and lessons learned. Nat. Rev. Drug Discov. 2023, 22, 449–475. [Google Scholar] [CrossRef] [PubMed]
  16. Marangoni, D.; Antonello, R.M.; Coppi, M.; Palazzo, M.; Nassi, L.; Streva, N.; Povolo, L.; Malentacchi, F.; Zammarchi, L.; Rossolini, G.M.; et al. Combination regimen of nirmatrelvir/ritonavir and molnupiravir for the treatment of persistent SARS-CoV-2 infection: A case report and a scoping review of the literature. Int. J. Infect. Dis. 2023, 133, 53–56. [Google Scholar] [CrossRef] [PubMed]
  17. Menegale, F.; Manica, M.; Zardini, A.; Guzzetta, G.; Marziano, V.; D'Andrea, V.; Trentini, F.; Ajelli, M.; Poletti, P.; Merler, S. Evaluation of Waning of SARS-CoV-2 Vaccine–Induced Immunity. JAMA Netw. Open 2023, 6, e2310650. [Google Scholar] [CrossRef] [PubMed]
  18. Mikulska, M.; Sepulcri, C.; Dentone, C.; Magne, F.; Balletto, E.; Baldi, F.; Labate, L.; Russo, C.; Mirabella, M.; Magnasco, L.; et al. Triple Combination Therapy With 2 Antivirals and Monoclonal Antibodies for Persistent or Relapsed Severe Acute Respiratory Syndrome Coronavirus 2 Infection in Immunocompromised Patients. Clin. Infect. Dis. 2023, 77, 280–286. [Google Scholar] [CrossRef]
  19. Nguyenla, X.; Wehri, E.; Van Dis, E.; Biering, S.B.; Yamashiro, L.H.; Zhu, C.; Stroumza, J.; Dugast-Darzacq, C.; Graham, T.G.W.; Wang, X.; et al. Discovery of SARS-CoV-2 antiviral synergy between remdesivir and approved drugs in human lung cells. Sci. Rep. 2022, 12, 1–17. [Google Scholar] [CrossRef]
  20. Qu, P.; Faraone, J.N.; Evans, J.P.; Zheng, Y.-M.; Carlin, C.; Anghelina, M.; Stevens, P.; Fernandez, S.; Jones, D.; Panchal, A.R.; et al. Enhanced evasion of neutralizing antibody response by Omicron XBB.1.5, CH.1.1, and CA.3.1 variants. Cell Rep. 2023, 42, 112443. [Google Scholar] [CrossRef]
  21. Semenzato, L.; Botton, J.; Drouin, J.; Baricault, B.; Bertrand, M.; Jabagi, M.-J.; Cuenot, F.; Le Vu, S.; Dray-Spira, R.; Weill, A.; et al. Characteristics associated with the residual risk of severe COVID-19 after a complete vaccination schedule: A cohort study of 28 million people in France. Lancet Reg. Heal. - Eur. 2022, 19, 100441. [Google Scholar] [CrossRef]
  22. Shyr, Z.A.; Cheng, Y.-S.; Lo, D.C.; Zheng, W. Drug combination therapy for emerging viral diseases. Drug Discov. Today 2021, 26, 2367–2376. [Google Scholar] [CrossRef] [PubMed]
  23. Touret, F.; Giraud, E.; Bourret, J.; Donati, F.; Tran-Rajau, J.; Chiaravalli, J.; Lemoine, F.; Agou, F.; Simon-Loriere, E.; van der Werf, S.; et al. Enhanced neutralization escape to therapeutic monoclonal antibodies by SARS-CoV-2 omicron sub-lineages. iScience 2023, 26, 106413. [Google Scholar] [CrossRef] [PubMed]
  24. Turtle, L.; Thorpe, M.; Drake, T.M.; Swets, M.; Palmieri, C.; Russell, C.D.; Ho, A.; Aston, S.; Wootton, D.G.; Richter, A.; et al. Outcome of COVID-19 in hospitalised immunocompromised patients: An analysis of the WHO ISARIC CCP-UK prospective cohort study. PLOS Med. 2023, 20, e1004086. [Google Scholar] [CrossRef] [PubMed]
  25. Vicenti, I.; Martina, M.G.; Boccuto, A.; De Angelis, M.; Giavarini, G.; Dragoni, F.; Marchi, S.; Trombetta, C.M.; Crespan, E.; Maga, G. System-oriented optimization of multi-target 2,6-diaminopurine derivatives: Easily accessible broad-spectrum antivirals active against flaviviruses, influenza virus and SARS-CoV-2. Eur. J. Med. Chem. 2021, 224, 113683. [Google Scholar] [CrossRef]
  26. Wagoner, J.; Herring, S.; Hsiang, T.-Y.; Ianevski, A.; Biering, S.B.; Xu, S.; Hoffmann, M.; Pöhlmann, S.; Gale, M.; Aittokallio, T.; et al. Combinations of Host- and Virus-Targeting Antiviral Drugs Confer Synergistic Suppression of SARS-CoV-2. Microbiol. Spectr. 2022, 10, e0333122. [Google Scholar] [CrossRef] [PubMed]
  27. Wang, X.; Sacramento, C.Q.; Jockusch, S.; Chaves, O.A.; Tao, C.; Fintelman-Rodrigues, N.; Chien, M.; Temerozo, J.R.; Li, X.; Kumar, S.; et al. Combination of antiviral drugs inhibits SARS-CoV-2 polymerase and exonuclease and demonstrates COVID-19 therapeutic potential in viral cell culture. Commun. Biol. 2022, 5, 1–14. [Google Scholar] [CrossRef] [PubMed]
  28. Wooten, D.J.; Meyer, C.T.; Lubbock, A.L.R.; Quaranta, V.; Lopez, C.F. MuSyC is a consensus framework that unifies multi-drug synergy metrics for combinatorial drug discovery. Nat. Commun. 2021, 12, 1–16. [Google Scholar] [CrossRef] [PubMed]
  29. Wu, M.Y.; Carr, E.J.; Harvey, R.; Mears, H.V.; Kjaer, S.; Townsley, H.; Hobbs, A.; Ragno, M.; Herman, L.S.; Adams, L.; et al. WHO's Therapeutics and COVID-19 Living Guideline on mAbs needs to be reassessed. Lancet 2022, 400, 2193–2196. [Google Scholar] [CrossRef]
  30. Yadav, B.; Wennerberg, K.; Aittokallio, T.; Tang, J. Searching for Drug Synergy in Complex Dose–Response Landscapes Using an Interaction Potency Model. Comput. Struct. Biotechnol. J. 2015, 13, 504–513. [Google Scholar] [CrossRef]
  31. Yan, D.; Yan, B. Viral target and metabolism-based rationale for combined use of recently authorized small molecule COVID-19 medicines: Molnupiravir, nirmatrelvir, and remdesivir. Fundam. Clin. Pharmacol. 2023, 37, 726–738. [Google Scholar] [CrossRef] [PubMed]
  32. Yang, S.; Multani, A.; Garrigues, J.M.; Oh, M.S.; Hemarajata, P.; Burleson, T.; Green, N.M.; Oliai, C.; Gaynor, P.T.; Beaird, O.E.; et al. Transient SARS-CoV-2 RNA-Dependent RNA Polymerase Mutations after Remdesivir Treatment for Chronic COVID-19 in Two Transplant Recipients: Case Report and Intra-Host Viral Genomic Investigation. Microorganisms 2023, 11, 2096. [Google Scholar] [CrossRef] [PubMed]
  33. Zhu, Y.; Binder, J.; Yurgelonis, I.; Rai, D.K.; Lazarro, S.; Costales, C.; Kobylarz, K.; McMonagle, P.; Steppan, C.M.; Aschenbrenner, L.; et al. Generation of a VeroE6 Pgp gene knock out cell line and its use in SARS-CoV-2 antiviral study. Antivir. Res. 2022, 208, 105429. [Google Scholar] [CrossRef]
  34. Zuckerman, N.S.; Bucris, E.; Keidar-Friedman, D.; Amsalem, M.; Brosh-Nissimov, T. Nirmatrelvir Resistance—de Novo E166V/L50V Mutations in an Immunocompromised Patient Treated With Prolonged Nirmatrelvir/Ritonavir Monotherapy Leading to Clinical and Virological Treatment Failure—a Case Report. Clin. Infect. Dis. 2023, 78, 352–355. [Google Scholar] [CrossRef] [PubMed]
Preprints 94871 g001
Figure 1. Dose-response matrix, two-dimensional (2D) and three-dimensional (3D) synergy maps generated by SinergyFinder 3.0 by applying the ZIP model on the Nirmatrelvir (NRM)/Remdesivir (RDV) combination tested against the ancestral B.1 SARS-COV-2 strain.
Figure 1. Dose-response matrix, two-dimensional (2D) and three-dimensional (3D) synergy maps generated by SinergyFinder 3.0 by applying the ZIP model on the Nirmatrelvir (NRM)/Remdesivir (RDV) combination tested against the ancestral B.1 SARS-COV-2 strain.
Preprints 94871 g001
Preprints 94871 g002
Figure 2. IC50 synergistic potency shift experiments in B.1 or BQ.1.1 virus infected VERO E6 cells treated with three fixed drug concentrations of Compound 1 (as indicated in the graph title) plus scalar concentrations of Compound 2 (as indicated in the legend). The three concentrations of Compound 1 were chosen as the IC50 of Compound 1 alone plus two 2-fold dilutions. The scalar concentration of Compound 2 is the same as that used in the SynergyFinder matrix. The concentration of Compound 2 on the x-axis is expressed as micromolar for DAA and as ng/ml for mAb while the y-axis shows the percentage of cell viability. DAA/DAA combinations were tested against the ancestral B.1 SARS-COV-2 strain and the BQ.1.1 Omicron variant while RDV/mAb combinations were tested only against the ancestral B.1 SARS-COV-2 strain. Graphs were generated by GraphPad PRISM software version 9 (La Jolla, CA, USA).
Figure 2. IC50 synergistic potency shift experiments in B.1 or BQ.1.1 virus infected VERO E6 cells treated with three fixed drug concentrations of Compound 1 (as indicated in the graph title) plus scalar concentrations of Compound 2 (as indicated in the legend). The three concentrations of Compound 1 were chosen as the IC50 of Compound 1 alone plus two 2-fold dilutions. The scalar concentration of Compound 2 is the same as that used in the SynergyFinder matrix. The concentration of Compound 2 on the x-axis is expressed as micromolar for DAA and as ng/ml for mAb while the y-axis shows the percentage of cell viability. DAA/DAA combinations were tested against the ancestral B.1 SARS-COV-2 strain and the BQ.1.1 Omicron variant while RDV/mAb combinations were tested only against the ancestral B.1 SARS-COV-2 strain. Graphs were generated by GraphPad PRISM software version 9 (La Jolla, CA, USA).
Preprints 94871 g002
Table 1. Summary synergy scores (SS) and weighted SS of all the tested antiviral combinations, generated by the Synergy Finder 3.0 online software (https://synergyfinder.fimm.fi/). The antiviral pairwise combinations of nirmatrelvir (NRM), EIDD-1931 (the active form of molnupiravir) and remdesivir (RDV) were tested against the ancestral B.1 and the Omicron BQ.1.1 SARS-CoV-2 variant. The combinations of RDV with each of the four monoclonal antibodies, including sotrovimab (SOT), bebtelovimab (BEB), cilgavimab (CIL) and tixagevimab (TIX), were tested only against the ancestral B.1 SARS-CoV-2 variant. NP: Not performed.
Table 1. Summary synergy scores (SS) and weighted SS of all the tested antiviral combinations, generated by the Synergy Finder 3.0 online software (https://synergyfinder.fimm.fi/). The antiviral pairwise combinations of nirmatrelvir (NRM), EIDD-1931 (the active form of molnupiravir) and remdesivir (RDV) were tested against the ancestral B.1 and the Omicron BQ.1.1 SARS-CoV-2 variant. The combinations of RDV with each of the four monoclonal antibodies, including sotrovimab (SOT), bebtelovimab (BEB), cilgavimab (CIL) and tixagevimab (TIX), were tested only against the ancestral B.1 SARS-CoV-2 variant. NP: Not performed.
Summary SS against B.1
Mean ± SD		 Weighted SS against B.1 Mean ± SD Summary SS against BQ.1.1 Mean ± SD Weighted SS against BQ.1.1 Mean ± SD
NRM/EIDD-1931 -3.1 ± 3.1 2.1 ± 0.4 -2.7 ± 1.3 0.0 ± 0.7
RDV/EIDD-1931 0.4 ± 6.2 0.9 ± 0.7 2.6 ± 6.5 1.7 ± 2.5
NRM/RDV 2.0 ± 2.0 1.9 ± 1.4 -0.2 ± 3.7 2.1 ± 2.1
RDV/SOT 8.8 ± 10.3 6.1 ± 2.9 NP NP
RDV/BEB 5.7 ± 2.6 4.2 ± 3.3 NP NP
RDV/TIX 0.6 ± 6.4 2.6 ± 2.7 NP NP
RDV/CIL -2.0 ± 0.5 0.8 ± 0.1 NP NP
Table 2. Concentrations of pairwise drug combinations below each drug IC90 and showing synergism in the SynergyFinder 3.0 matrix, stratified by synergy score values (10-30, 30-50, and above 50).
Table 2. Concentrations of pairwise drug combinations below each drug IC90 and showing synergism in the SynergyFinder 3.0 matrix, stratified by synergy score values (10-30, 30-50, and above 50).
Synergy score stratum Drug combination Concentrationsa Viral variant
>50 RDV/NRM 0.02/0.05 B.1 and BQ1.1
RDV/SOT 0.04/8.8
0.04/35.2
0.04/140.6		 B.1
RDV/BEB 0.04/0.6
0.04/2.2
0.04/8.9		 B.1
RDV/TIX 0.04/2.2
0.04/8.8		 B.1
30-50 NRM/EIDD-1931 0.05/0.75
0.05/1.50		 B.1
B.1 and BQ.1.1
RDV/SOT 0.01/562.5
0.04/562.5		 B.1
B.1
RDV/TIX 0.04/35 B.1
10-29.9 NRM/EIDD-1931 0.003/0.07
0.01/0.08
0.05/0.09		 B.1
B.1
B.1
RDV/EIDD-1931 0.02/1.5
0.03/0.75		 B.1 and BQ.1.1
BQ.1.1
RDV/NRM 0.02/0.05
0.004/0.05
0.006/0.05
0.02/0.03
0.02/0.01		 B.1
B.1 and BQ.1.1
BQ.1.1
BQ.1.1
BQ.1.1
RDV/BEB 0.01/8.9 B.1
RDV/TIX 0.01/35 B.1
aConcentrations are expressed as µM for DAA and ng/ml for mAb.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

facebook logotwitter logolinkedin logo
MDPI Initiatives
Important Links
Subscribe

Disclaimer

Terms of Use

Privacy Policy

© 2024 MDPI (Basel, Switzerland) unless otherwise stated