Screening for drugs that don’t work
In the battle against COVID-19, drugs discovered in repurposing screens are of particular interest because these could be rapidly implemented as treatments. However, Tummino et al. deliver a cautionary tale, finding that many leads from such screens have an antiviral effect in cells through phospholipidosis, a phospholipid storage disorder that can be induced by cationic amphiphilic drugs (see the Perspective by Edwards and Hartung). There is a strong correlation between drug-induced phospholipidosis and inhibition of severe acute respiratory syndrome coronavirus 2 replication in cells. Unfortunately, drugs that have an antiviral effect in cells through phospholipidosis are unlikely to be effective in vivo. Screening out such drugs may allow a focus on drugs with better clinical potential.
Science, abi4708, this issue p. 541; see also abj9488, p. 488
Abstract
Repurposing drugs as treatments for COVID-19, the disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has drawn much attention. Beginning with sigma receptor ligands and expanding to other drugs from screening in the field, we became concerned that phospholipidosis was a shared mechanism underlying the antiviral activity of many repurposed drugs. For all of the 23 cationic amphiphilic drugs we tested, including hydroxychloroquine, azithromycin, amiodarone, and four others already in clinical trials, phospholipidosis was monotonically correlated with antiviral efficacy. Conversely, drugs active against the same targets that did not induce phospholipidosis were not antiviral. Phospholipidosis depends on the physicochemical properties of drugs and does not reflect specific target-based activities—rather, it may be considered a toxic confound in early drug discovery. Early detection of phospholipidosis could eliminate these artifacts, enabling a focus on molecules with therapeutic potential.
The outbreak of COVID-19 has inspired multiple drug repurposing screens to find antiviral therapeutics that can be rapidly brought to the clinic (1). To date, more than 1974 drugs and investigational drugs have been reported to have in vitro activity against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) (1) (Fig. 1). Because almost all of these drugs act against human targets and might be unlikely to be viable against a novel virus (2), the question of mechanism of action arises.
Our interest in this question was motivated by the discovery that human sigma receptors were candidates for modulating SARS-CoV-2 infection (3) and that drugs and reagents like chloroquine, haloperidol, clemastine, and PB28—all with nanomolar affinity against one or both sigma receptors—had cellular antiviral half-maximal inhibitory concentration (IC50) values in the 300-nM to 5-μM range. Subsequently, we investigated more than 50 different molecules with a wide range of affinities at these receptors. Although we found molecules with relatively potent antiviral activity, there was little correlation between receptor potency and antiviral efficacy in cells (figs. S1 to S3 and table S1). Whereas drugs like amiodarone, sertraline, and tamoxifen had mid- to high-nanomolar antiviral IC50s, other potent sigma ligands, such as melperone and ditolylguanidine (DTG), were without measurable antiviral activity. Notably, the antiviral sigma drugs were all cationic at physiological pH and relatively hydrophobic, whereas those that were inactive against the virus were often smaller and more polar. This cationic-amphiphilic character was shared by many of the hits emerging from other phenotypic screens (Fig. 1 and fig. S4), suggesting that it was this physicochemical property that might explain cellular antiviral activity instead of a specific on-target activity (4).
If the cationic-amphiphilic nature of these molecules led to antiviral activity in vitro, rather than their target-based activities, one would expect this physical property to reflect a shared cellular mechanism. In fact, cationic amphiphilic drugs (CADs) can provoke phospholipidosis in cells and organs (5). This side effect is characterized by the formation of vesicle-like structures and “foamy” or “whorled” membranes (5, 6) and is thought to arise by CAD disruption of lipid homeostasis. CADs accumulate in intracellular compartments, such as endosomes and lysosomes, where they can directly or indirectly inhibit lipid processing (5). Modulation of these same lipid processing pathways is critical for viral replication (7), and inhibiting phospholipid production has previously been associated with the inhibition of coronavirus replication (8). CADs have in vitro activity against multiple viruses including severe acute respiratory syndrome, Middle East respiratory syndrome, Ebola, Zika, dengue, and filoviruses (9), though CAD induction of phospholipidosis has only been proposed as an antiviral mechanism for Marburg virus (10). Finally, among the most potent known phospholipidosis inducers are amiodarone (11) and chloroquine (12, 13), which are also potent inhibitors of SARS-CoV-2 replication in vitro (14–16), whereas drugs from SARS-CoV-2 phenotypic screens, such as chlorpromazine (17) and tamoxifen (16), are also known to induce phospholipidosis (18). As an effect that rarely occurs at concentrations lower than 100 nM, which does not appear to translate from in vitro to in vivo antiviral activity and which can result in dose-limiting toxicity (19), phospholipidosis may be a confound to true antiviral drug discovery.
Here, we investigate the association between phospholipidosis and antiviral activity against SARS-CoV-2 in cell culture. This apparently general mechanism may be responsible for many of the drug repurposing hits for SARS-CoV-2 and an extraordinary amount of effort and resources lavished on drug discovery against this disease. We explore the prevalence of this confound in SARS-CoV-2 repurposing studies, how phospholipidosis correlates with inhibition of viral infection, and how to eliminate such hits rapidly so as to focus on drugs with genuine potential against COVID-19 and against pandemics yet to arise.
Results
Correlation of phospholipidosis and antiviral activity
To investigate the role of phospholipidosis in antiviral activity in vitro, we tested 19 drugs for their induction of this effect in A549 cells using the well-established nitrobenzoxadiazole–conjugated phosphoethanolamine (NBD-PE) staining assay (20). Here, the vesicular lipidic bodies characteristic of the effect may be quantified by high-content imaging (Fig. 2A). Three classes of drugs and reagents were initially investigated: (i) sigma-binding antiviral CADs we had discovered, like amiodarone, sertraline, chlorpromazine, and clemastine (nine total), which are predicted or known to induce phospholipidosis; (ii) analogs of these CADs that no longer bound sigma receptors but were still antiviral (four total), which are predicted to induce phospholipidosis despite their lack of sigma binding; and (iii) sigma-binding, nonantiviral drugs, like melperone and DTG, that are more polar than classic CADs (two total), which are predicted not to induce phospholipidosis. Of the nine sigma-binding CADs that were antiviral (the first class)—six of which were also reported in phenotypic screens in the literature as inhibitors of SARS-CoV-2—eight induced phospholipidosis consistent with the hypothesis (Fig. 2, A and B, and figs. S5 and S6). The only non–phospholipidosis-inducing antiviral from this set was elacridar, a promiscuous P-glycoprotein inhibitor; this drug may therefore be active through another mechanism. Notably, analogs of the potent sigma ligand PB28 that had lost their sigma-binding activity but remained CADs (ZZY-10-051 and ZZY-10-061; Fig. 2, B to F, and figs. S5 to S8) did induce phospholipidosis, as did the antipsychotic olanzapine and the antihistamine diphenhydramine, which are weak sigma receptor ligands but are structurally related to potent sigma receptor ligands like chlorpromazine and clemastine. Finally, melperone and DTG, which are potent cationic sigma receptor ligands but are not antiviral, did not induce phospholipidosis (Fig. 2, A and B, and figs. S5 and S6; class iii). These results do not prove phospholipidosis as the antiviral mechanism but are consistent with the phospholipidosis hypothesis.
If phospholipidosis is responsible for antiviral activity, then other molecules known to induce phospholipidosis should be antiviral. We tested three CADs for antiviral activity, including ebastine, ellipticine, and Bix01294, all of which are reported to induce phospholipidosis (21) [Bix01294 and ebastine have also been reported as drug repurposing hits against SARS-CoV-2 (22)]. We further tested azithromycin, which is also reported to induce phospholipidosis (23) but has different physical properties than typical CADs. We first confirmed phospholipidosis-inducing activity for these molecules (Fig. 2B and figs. S5 and S6). All four molecules were next shown to be antiviral using live virus assays (e.g., SARS-CoV-2 strain BetaCoV/France/IDF0372/2020; see materials and methods), with IC50 values in the 400-nM to 3-μM range, overlapping with the activities of other CADs that we and others have identified for SARS-CoV-2 (22) (fig. S6). This too was consistent with the antiviral phospholipidosis hypothesis.
For phospholipidosis to explain antiviral activity, we might expect a correlation between concentration-response curves for phospholipidosis and for antiviral activity. We compared concentrations that induce phospholipidosis with those that inhibit SARS-CoV-2 for each drug individually. Most correlations were high—not only did antiviral activity occur in the same concentration ranges as phospholipidosis, but R2 (coefficient of determination) values, ranging from 0.51 to 0.94, supported a quantitative relationship between the two effects (Fig. 3A). We then fit a sigmoidal model through all of the 107 phospholipidosis versus antiviral activity observations (made up of six concentration measurements for each of the 16 phospholipidosis-inducing drugs) and observed a strong negative correlation [R2 = 0.65; 95% confidence interval (95% CI), 0.52 to 0.76] between induced phospholipidosis and SARS-CoV-2 viral load across all observations for all 16 drugs. Because phospholipidosis and antiviral effects are both saturable, the sigmoidal curve–fit plateaus at the extremes (Fig. 3B).
(A) Correlations between phospholipidosis (DIPL), normalized to amiodarone at 100%, and percent of SARS-CoV-2, normalized to DMSO at 100%, in the reverse transcription quantitative polymerase chain reaction (RT-qPCR) assay in A549-ACE2 cells. Each dot represents the same concentration tested in both assays. A strong negative correlation emerges, with R2 ≥ 0.65 and P ≤ 0.05 for all high and medium phospholipidosis-inducing drugs except ellipticine, which is confounded by its cytotoxicity in both experiments; ebastine; and ZZY-10-61. The latter two examples are marginally significant. (B) The SARS-CoV-2 viral loads and induced phospholipidosis magnitude for each compound and dose in (A) are plotted as sqrt(viral_amount_mean) ~ 10 × inv_logit{hill × 4/10 × [log(DIPL_mean) − logIC50]}. Fitting a sigmoid Bayesian model with weakly informative priors yields parameters and 95% credible intervals of IC50 = 43 (38, 48)%, hill: −5.6 (−7.0, −4.5), and sigma 2.0 (0.14, 1.78). Forty draws from the fit model are shown as blue lines. Salmon-colored points overlaid with the model represent predicted phospholipidosis inducers from the literature (fig. S10).
Concurrent measurement of viral infection and drug-induced phospholipidosis
In the previous experiments, drug-induced phospholipidosis and drug antiviral activity were measured separately. To measure the two effects in the same cells at the same time, we dosed cells with either 1 or 10 μM of five characteristic CADs [amiodarone, sertraline, PB28, hydroxychloroquine (HCQ), and Bix01294], followed by a mock or SARS-CoV-2 infection, and quantified phospholipidosis and the accumulation of viral spike protein (Fig. 4A and fig. S9). Compared with dimethyl sulfoxide (DMSO), drug treatments led to substantial increases in NBD-PE aggregates, indicating increased phospholipidosis (fig. S9). At 1-μM drug concentrations, SARS-CoV-2 spike protein was readily stained, and it was possible to visualize both spike protein and phospholipidosis in the same cells (yellow puncta), which suggests that at this low concentration of drug—often close to the antiviral IC50 value—both phospholipidosis and viral infection co-occur, though even viral staining was reduced relative to that observed in the DMSO-treated controls. As drug concentration rose to 10 μM, viral spike protein staining dropped, whereas staining for phospholipidosis increased (fig. S9); there was nearly complete loss of spike protein signal with a concomitant increase in phospholipidosis (Fig. 4A) for all treatments. In seven-point concentration-response curves for amiodarone, sertraline, and PB28, viral staining monotonically decreased as phospholipidosis increased (Fig. 4, B and C).
(A) Representative images from a costaining experiment measuring phospholipidosis and SARS-CoV-2 spike protein in infected and uninfected A549-ACE2 cells. Five molecules (1 and 10 μM) and DMSO were measured; see fig. S9 for Bix 01294. Blue indicates Hoechst nuclei staining, green indicates NBD-PE phospholipid staining, red indicates SARS-CoV-2 spike protein staining, and yellow indicates coexpression of spike protein and NBD-PE. Scale bar, 20 μm. (B) Concentration-response curves for phospholipidosis induction measured by NBD-PE staining in infected cells for three characteristic CADs. (C) Spike protein in infected cells decreases as phospholipidosis increases. For (B) and (C), data are means ± SEMs from four biological replicates.