In vitro assays to measure virus replication have been greatly improved by the development of recombinant RNA viruses expressing luciferase or other enzymes capable of bioluminescence. Here we detail a high-throughput screening pipeline that combines such recombinant strains of measles and chikungunya viruses to isolate broad-spectrum antivirals from chemical libraries.
RNA viruses are responsible for major human diseases such as flu, bronchitis, dengue, Hepatitis C or measles. They also represent an emerging threat because of increased worldwide exchanges and human populations penetrating more and more natural ecosystems. A good example of such an emerging situation is chikungunya virus epidemics of 2005-2006 in the Indian Ocean. Recent progresses in our understanding of cellular pathways controlling viral replication suggest that compounds targeting host cell functions, rather than the virus itself, could inhibit a large panel of RNA viruses. Some broad-spectrum antiviral compounds have been identified with host target-oriented assays. However, measuring the inhibition of viral replication in cell cultures using reduction of cytopathic effects as a readout still represents a paramount screening strategy. Such functional screens have been greatly improved by the development of recombinant viruses expressing reporter enzymes capable of bioluminescence such as luciferase. In the present report, we detail a high-throughput screening pipeline, which combines recombinant measles and chikungunya viruses with cellular viability assays, to identify compounds with a broad-spectrum antiviral profile.
RNA viruses are responsible for a large variety of human infections, and have a huge impact on populations worldwide both in terms of public health and economical cost. Efficient vaccines have been developed against several human RNA viruses, and are widely used as prophylactic treatments. However, there is still a critical lack of therapeutic drugs against RNA virus infections. Indeed, efficient vaccines are not available against major human pathogens such as dengue virus, Hepatitis C virus or human respiratory syncytial virus (hRSV). Besides, RNA viruses are responsible for a majority of emerging diseases, which have increased in frequency because of global exchanges and human impact on ecological systems. Against this threat that represent RNA viruses, our therapeutic arsenal is extremely limited and relatively inefficient1-3. Current therapies are essentially based on recombinant type I interferons (IFN-α/β) to stimulate innate immunity, or the administration of ribavirin. Although the mode of action of this ribonucleoside analog is controversial and probably relies on various mechanisms, the inhibition of cellular IMPDH (inosine monophosphate dehydrogenase), which depletes intracellular GTP pools, is clearly essential4. Ribavirin, in combination with pegylated IFN-α, is the main treatment against Hepatitis C virus. However, IFN-α/β and ribavirin treatments are of relatively poor efficacy in vivo against most RNA viruses as they efficiently blunt IFN-α/β signaling through expression of virulence factors5 and often escape ribavirin3. This added to the fact that ribavirin treatment is raising important toxicity issues, although it was recently approved against severe hRSV disease with controversial benefits6. More recently, some virus-specific treatments have been marketed, in particular against influenza virus with the development of neuraminidase inhibitors3. However, the large diversity and permanent emergence of RNA viruses precludes the development of specific treatments against each one of them in a relatively close future. Altogether, this stresses the need for efficient strategies to identify and develop potent antiviral molecules in the near future.
It is trivial to say that a broad-spectrum inhibitor active against a large panel of RNA viruses would solve this problem. Although such a molecule is still a virologist's dream, our better understanding of cellular defense mechanisms and innate immune system suggest that some possibilities exist7,8. Several academic and industrial laboratories are now seeking molecules that stimulate specific facets of cellular defense mechanisms or metabolic pathways to blunt viral replication. Although such compounds will probably show significant side effects, treatments against acute viral infections will be administered for a relatively short time, making them acceptable despite some potential toxicity on the long term. Various strategies have been developed to identify such broad-spectrum antiviral molecules. Some research programs aim at finding molecules that target specific defense or metabolic pathways. This includes, for example, pathogen recognition receptors to elicit antiviral gene expression9 and activate antiviral factors such as RNaseL10, autophagy machinery to promote virus degradation11, nucleoside synthesis pathways12,13, or apoptotic cascades to precipitate death of virus-infected cells14. Other groups have developed phenotypic screens that are not target-based13,15-17. In that case, antiviral molecules are simply identified by their capacity to block viral replication in a given cellular system. The general assumption is that a compound inhibiting 2-3 unrelated RNA viruses would have a suitable profile for a broad-spectrum antiviral molecule. The mode of action of hit compounds selected with such an empirical approach is only determined in a second time and eventually, may lead to the identification of novel cellular targets for antivirals. Interestingly, a retrospective analysis of new drugs approved by the US Food and Drug Administration between 1999 and 2008 has shown that in general, such phenotypic screenings tend to perform better than target-based approaches to discover first-in-class small-molecule drugs18.
Viral replication in high-throughput cell-based assays is usually determined from virus cytopathic effects. Cells are infected and cultured in 96- or 384-well plates in the presence of tested compounds. After few days, cellular layers are fixed and stained with dyes such as crystal violet. Finally, absorbance is determined with a plate reader and compounds inhibiting viral replication are identified by their capacity to preserve cellular layers from virus-induced cytopathic effect. Alternatively, virally-mediated cytopathic effects are assessed using standard viability assays such as MTS reduction. Such assays are highly tractable and cost-effective, but suffer from three major limitations. First, they require a virus-cell combination where viral replication is cytopathic in only few days but this is not always possible, thus calling for alternative approaches19. Second, they are poorly quantitative since they are based on an indirect measure of viral replication. Finally, toxic compounds can be scored as positive hits, and therefore must be eliminated with a counter screen measuring cellular viability. To overcome some of these hurdles, recombinant viruses or replicons have been engineered by reverse genetics to express reporter proteins, such as EGFP or luciferase, from an additional transcription unit or in frame with viral protein genes (few examples are 20-23). When these viruses replicate, reporter proteins are produced together with viral proteins themselves. This provides a very quantitative assay to measure viral replication and evaluate the inhibitory activity of candidate molecules. This is particularly true for recombinant viruses expressing luciferase (or other enzymes capable of bioluminescence) since this reporter system exhibits a wide dynamic range with a high sensitivity and virtually no background. Furthermore, there is no excitation light source, thus preventing interference with compound fluorescence24.
Here, we detail a high-throughput protocol to screen chemical libraries for broad-spectrum inhibitors of RNA viruses. Compounds are tested first on human cells infected with a recombinant measles virus (MV) expressing firefly luciferase25 (rMV2/Luc, Figure 1, primary screen). MV belongs to Mononegavirales order, and is often considered as a prototypical member of negative-strand RNA viruses. As such, MV genome is used as a template by the viral polymerase to synthesize mRNA molecules encoding for viral proteins. In the recombinant MV strain called rMV2/Luc, luciferase expression is expressed from an additional transcription unit inserted between P and M genes (Figure 2A). In parallel, compounds are tested for their toxicity on human cells using a commercial luciferase-based reagent that evaluates, by ATP quantification, the number of metabolically active cells in culture (Figure 1, primary screen). Entire chemical libraries can be easily screened with these two assays in order to select compounds that are not toxic and efficiently block MV replication. Then, hits are retested for dose-response inhibition of MV replication, the lack of toxicity as well as for their capacity to impair chikungunya virus (CHIKV) replication (Figure 1, secondary screen). CHIKV is a member of Togaviridae family and its genome is a positive single-strand RNA molecule. As such, it is completely unrelated to measles and compounds inhibiting both MV and CHIKV stand a great chance to inhibit a large panel of RNA viruses. CHIKV nonstructural proteins are directly translated from the viral genome, whereas structural proteins are encoded by transcription and translation of a subgenomic mRNA molecule. Our in vitro replication assay for CHIKV is based on a recombinant strain called CHIKV/Ren, which expresses Renilla luciferase enzyme as a cleaved part of the nonstructural polyprotein through an insertion of the reporter gene between nsP3 and nsP4 sequences26 (Figure 2B). The measure of Renilla luciferase activity allows the monitoring of viral replication at the early stage of CHIKV life cycle.
This high-throughput protocol was used to quickly identify compounds with a suitable profile for broad-spectrum antivirals in a commercial library of 10,000 molecules enriched for chemical diversity. Compounds were essentially following Lipinski's rule of five, with molecular weights ranging from 250 to 600 daltons, and log D values below 5. Most of these molecules were new chemical entities not available in other commercial libraries.
1. Preparation of 96-well Daughter Plates with Compounds
2. Preparation of Cell Cultures to Determine Compound Toxicity
3. Preparation of Cell Cultures Infected by rMV2/Luc
4. Luciferase Activity Measures and Data Analysis
5. Secondary Screen for Toxicity and Broad-spectrum Antiviral Activity
This screening pipeline relies first on the selection of compounds that inhibit MV replication and do not show any significant cellular toxicity (Figure 1, primary screen). It takes advantage of luciferase-based assays to determine cellular toxicity and the inhibition of viral replication. All pipetting steps and data acquisition can be performed in a high-throughput setting with a robotic platform.
Cellular toxicity of compounds is assessed using a luciferase-based viability assay that evaluates the number of living cells in culture based on quantification of the ATP present, which corresponds to metabolically active cells. Because luciferin oxidation by luciferase is ATP dependent, luminescence signal is directly proportional to cellular ATP provided by living cells. Figure 3A is showing results for one representative screening plate. As expected, very distinct luciferase activity values were measured in control wells that contained either living cells or dead cells killed with IGEPAL. Threshold is empirically set at µ+/2 to filter out compounds showing some cellular toxicity and in the present case, 21 compounds from this plate were discarded. In parallel, antiviral activity is measured using rMV2/Luc, a recombinant strain of MV expressing luciferase (Figure 2). Human HEK-293T cells are used because they are highly sensitive to MV infection. As a result, luciferase activity in infected culture wells is very high (from 200 to 300 x 103 luciferase activity units in control wells). This provides a large dynamic range that facilitates the selection of MV replication inhibitors. Figure 3B is showing representative results obtained for one screening plate. Four compounds inhibited viral replication by more than 75% in this example and were selected.
Finally, compounds that are not toxic and scored positive in the viral replication assay are selected (Figure 3A and 3B). Their half maximal inhibitory concentration (IC50) is determined on both MV and CHIKV, which are two unrelated RNA viruses (Figure 1, secondary screen). Antiviral activity of hit compounds is determined using recombinant virus strains expressing either firefly (rMV2/Luc) or Renilla (CHIKV/Ren) luciferase. Figure 4A and 4B show representative inhibition curves for MV and CHIKV replication, respectively. These curves are used to determine IC50s of hit compounds. In parallel, the lack of cellular toxicity of selected compounds is confirmed in a dose-response experiment using the luciferase-based viability assay described above (Figure 4C). Table 1 is showing IC50 values for MV and CHIKV replication on a set of 13 non-toxic compounds identified from a chemical library of 10,000 molecules. These compounds are novel both in terms of chemical structure and biological activity. Some of these molecules shared structural similarities, and can be gathered into three distinct chemical families. As a reference, IC50 values obtained with brequinar are also shown in Table 1. Brequinar is an inhibitor of dihydroorotate dehydrogenase, the fourth enzyme of pyrimidine biosynthesis pathway, with a potent antiviral activity in vitro30,31. When looking at IC50 values, less than one order of magnitude separated this reference molecule from most active hits identified by screening. This established the strong antiviral activity of selected compounds.
Figure 1. Summary of the screening pipeline. Left panel is showing the primary screen where compounds are selected for the lack of toxicity and capacity to block MV replication. Right panel is showing the secondary screen where the antiviral activity of selected compounds is confirmed, and their capacity to block CHIKV replication is determined. Compounds that efficiently block MV and CHIKV without any significant toxicity are considered as potential broad-spectrum inhibitors of RNA viruses. Click here to view larger image.
Figure 2. Genomic organization of recombinant RNA viruses carrying luciferase gene. (A) rMV2/Luc: the MV negative-sense RNA genome is displayed with its 3' end on the left, with the six genes indicated in capital letters and depicted as white rectangles. The additional transcriptional units encoding for the reporter gene is inserted between P and M genes and depicted as a yellow rectangle. (B) CHIKV/Ren: the CHIKV positive-sense RNA genome is displayed with its 5'-capped end on the left, with the open reading frame encoding the four non-structural proteins (nsP1-4). Click here to view larger image.
Figure 3. Representative results for one 96-well screening plate. (A) Luciferase values obtained with the luciferase-based viability assay. Columns 1 and 12 correspond to alternate controls with DMSO alone ("+"; no toxicity) or IGEPAL ("-"; high toxicity). Blue color is highlighting compounds that are considered not toxic according to ATP concentrations detected in culture wells, whereas white wells correspond to toxic compounds (luciferase values < 12,203 x 103). (B) Luciferase values obtained with rMV2/Luc-infected cells. Columns 1 and 12 correspond to alternate controls, i.e. non-infected HEK-293T cells ("-") or rMV2/Luc-infected cells treated with DMSO alone ("+"). For each tested compound, the figure is showing both the luciferase luminescence signal and the percentage of inhibition relative to control wells. Compounds inhibiting luminescence signal by more than 75% are highlighted in yellow or green depending if they were considered toxic or not as determined in (A). Thus, green color is showing hit compounds that inhibit MV replication and do not affect cellular viability. Click here to view larger image.
Figure 4. Dose-response antiviral activity and toxicity of compounds C864, C648 and C062. (A) HEK-293T cells were infected with rMV2/Luc strain of MV expressing luciferase (MOI=0.1), and incubated with increasing doses of C864, C648, C062 or DMSO alone. After 24 hr, luciferase expression was determined. * indicates that observed luciferase inhibitions were statistically significant with all three compounds (p-values < 0.05). (B) HEK-293T cells were infected with CHIKV/Ren strain of CHIKV expressing Renilla luciferase (MOI=0.2), and incubated with increasing doses of C864, C648, C062 or DMSO alone. After 24 hr, Renilla luciferase expression was determined. * indicates statistically significant inhibitions with all three compounds (p-values < 0.05). (C) HEK-293T cells were incubated with increasing doses of C864, C648, C062 or DMSO alone. As a control for toxicity, cell cultures were supplemented with 2.5 μl of a 0.5% IGEPAL solution. After 24 hr, the number of living cells was determined using the luciferase-based viability assay. Click here to view larger image.
Chemical Family | Compound | MV (IC50 μM) | CHIKV (IC50 μM) | CC50 (μM) |
1 | C864 | 0.6 | 0.6 | > 5 |
1 | C877 | 0.2 | 0.2 | > 5 |
1 | C957 | 1.2 | 1.2 | > 5 |
1 | C963 | 1.2 | 0.9 | > 5 |
1 | C967 | 1.7 | 1.6 | > 5 |
1 | C348 | 0.25 | 0.3 | > 5 |
1 | C265 | 0.25 | 0.3 | > 5 |
1 | C270 | 0.5 | 0.5 | > 5 |
1 | C350 | 0.4 | 0.5 | > 5 |
2 | C646 | 0.16 | 0.15 | > 5 |
2 | C814 | 1.9 | 1.1 | > 5 |
2 | C648 | 0.7 | 0.4 | > 5 |
3 | C062 | 1 | 1.2 | > 5 |
Brequinar | 0.04 | 0.04 | > 5 |
Table 1. Antiviral activity and cytotoxicity of the 13 compounds selected from the original chemical library. Results are expressed as IC50 values for the inhibition of MV or CHIKV replication. CC50 values were found >5 μM for all compounds in dose-response experiments from the secondary screen. But the scoring filter threshold for cytoxicity from the primary screen suggest that CC50 values are even >20 μM.
The screening pipeline described here aims at selecting compounds with a suitable profile as broad-spectrum inhibitors of RNA viruses. When a library of 10,000 compounds was screened with this protocol and aforementioned filtering criteria were applied, i.e. inhibition of MV replication superior to 75%, 40 compounds (0.4%) scored positive. Besides, about half of them showed some toxicity in the luciferase-based viability assay, and were disregarded for this reason. Finally, a dozen of hits readily available in larger amount were retested. This small set was easily retested for both cellular toxicity and the inhibition of viral replication against MV and CHIKV in dose-response experiments (Figure 4). The antiviral activity of all selected compounds was confirmed, and the overall identification hit rate for this screening pipeline was 1.3‰. This score is very similar to results obtained by other groups with different high-throughput replication-based screening assays for antivirals27,32, and falls in the expected range for cell-based phenotypic assays in general (1-5‰33). However, a hit rate is highly dependent on both the compound library that is screened and the scoring filter thresholds, which can be adjusted depending on the profile of collected data. Here, selected thresholds for the inhibition of MV replication and cellular toxicity were empirically determined, and trade off was set to obtain a tractable number of hits that can be easily retested in the secondary screen against CHIKV.
In the primary screen, compounds were tested for both their capacity to block MV replication and the lack of cellular toxicity. Antiviral activity was determined at a concentration of 2 μM on rMV2/Luc-infected cells whereas cellular toxicity was determined at 20 μM. These experimental settings were found valuable to directly select for compounds with a high antiviral activity relative to cellular toxicity. In addition, the primary screen was performed in the presence of uridine. Indeed, recent publications suggest that compounds targeting early steps of pyrimidine biosynthesis pathway are frequently isolated when looking for inhibitors of RNA viruses13,16,27,28. In particular, research groups looking for antiviral molecules have found inhibitors of dihydroorotate dehydrogenase, the fourth enzyme of this metabolic pathway31. Uridine in culture medium efficiently trans-complements for the inhibition of this enzyme, and thus restores viral replication. In this way, compounds that prevent MV replication through inhibition of pyrimidine biosynthesis pathway are directly excluded from the antiviral hits. Finally, the Z’-factor was found to be systematically above 0.5, thus demonstrating the high robustness and quality of the assay.
In parallel, cellular toxicity was determined using a high-throughput luciferase-based assay that determines the number of metabolically active cells in culture based on quantification of the ATP. In each screening plate, IGEPAL was added in control wells for toxicity. This non-denaturing, non-ionic detergent has the advantage of killing efficiently all cell types non-specifically without affecting luciferase enzyme activity. Z’-factor was systematically above 0.5 over 125 plates, thus demonstrating the high robustness and quality of the assay. One limitation of this cellular toxicity assay is its relatively high cost. As an alternative assay to determine cellular viability, some groups have used cell lines constitutively expressing luciferase from a stable transgene34. In such assays, toxic compounds are expected to decrease or even suppress luciferase expression when affecting cellular viability. We recently developed a stable cell line expressing luciferase upon IFN-β stimulation35, and this cell line was used to test a training set of 80 compounds for cellular toxicity. Surprisingly, the correlation coefficient (CC) between this assay and the luciferase-based viability assay described above was relatively low (CC = 0.67). Indeed, several compounds impaired luciferase expression from the transgene, whereas ATP metabolism was unaffected. Some broad-spectrum antiviral compounds are likely to affect cellular functions such as cell signaling, gene transcription or protein translation, which will also impair cellular gene expression. Although there is no perfect screening system for cellular viability, it is at risk that filtering toxic compounds based on cellular gene expression will potentially lead to exclude bona fide antiviral molecules.
The secondary screen aims at establishing the antiviral activity of selected compounds against two completely unrelated RNA viruses, while determining precisely their IC50s. Interestingly, it should be noticed that luciferase enzyme inhibitors could be wrongly scored as positives in the rMV2/Luc assay but such false-positives will be filtered-out by the secondary screen when tested for the inhibition of CHIKV/Ren. Indeed, rMV2/Luc and CHIKV/Ren express firefly and Renilla luciferase, respectively, and these two enzymes are unrelated and catalyze different chemical reactions. Compounds that inhibit firefly luciferase reporter from rMV2/Luc in the primary screen will not show any activity against CHIKV/Ren and will be discarded, thus preventing their erroneous selection as antivirals. Unexpectedly, all 13 primary hits blocking MV replication also inhibited CHIKV in this secondary screen. This is not always the case as few MV-specific inhibitors were isolated when screening other libraries (data not shown). Furthermore, selected compounds share structural similarities and can be grouped into only three chemical families (Table 1). Compounds within a chemical family probably have the same mode of action, and therefore can be considered as analogs of the same chemical entity. This puts into perspective the fact that all 13 selected compounds inhibited both MV and CHIKV replication. Nevertheless, these results could also suggest that broad-spectrum antivirals are more frequently identified than virus-specific inhibitors. The rational behind this hypothesis is that, despite obvious specificities, all RNA viruses rely on the same basic cellular functions to replicate, such as ribosomal machinery, cytoskeleton or mitochondria as a source of energy. These functional modules are complex molecular systems, which provide a relatively large panel of potential targets for broad-spectrum antivirals. In contrast, there is only a limited set of viral or cellular factors that represent appropriate targets to raise virus-specific inhibitors. In the future, new panels of antivirals identified by phenotypic screening should provide necessary data to confirm or invalidate this hypothesis.
As aforementioned, isolated compounds had a suitable profile for broad-spectrum inhibitors of RNA viruses with several showing IC50 values below 1 μM in both MV and CHIKV replication assays. This will need to be confirmed with more conventional viral replication assays, including measures of viral progeny in culture supernatants. Antiviral activity of selected compounds will also have to be confirmed on a much larger panel of RNA viruses. In conclusion, these compounds represent excellent starting points for their pharmaceutical development, with the demonstration of their antiviral activity in animal infection models as a primary goal.
The authors have nothing to disclose.
We thank Dr. Yves L. Janin for his fruitful comments and suggestions. We would like to thank H.H. Gad for CHIK/Ren and C. Combredet for her technical support. This work was supported by the Institut Pasteur, the Centre National de la Recherche Scientifique (CNRS), the Institut National de la Santé Et de la Recherche Médicale (INSERM), the Institut Carnot – Pasteur Maladies Infectieuses (Programme STING to P.O.V and H.M.-L), the Agence Nationale pour la Recherche (ANR-RPIB, Programme STING 2.0 to P.O.V), and the "Conseil Régional d’Ile-de-France" (Chemical Library Project, grants n° I 06-222/R and I 09-1739/R to H.M.-L.). The work on CHIKV/Ren was supported by the project ArbOAS (ANR grant 2010-INTB-1601-02).
Freedom EVO platform | TECAN | Robotic platform | |
96-Well Polystyrene Cell Culture Microplates, White | Greiner Bio One | 655083 | |
CellTiter-Glo Luminescent Cell Viability Assay | Promega | G7570 | Luciferase-based viability assay |
Bright-Glo Luciferase Assay System | Promega | E2610 | Reagent containing firefly luciferase substrate |
Renilla-Glo Luciferase Assay System | Promega | E2710 | Reagent containing Renilla luciferase substrate |
Britelite plus Reporter Gene Assay System | Perkin-Elmer | 6016761 | Reagent containing firefly luciferase substrate. Can be used as an alternative to Brigh-Glo reagent to determine luciferase activity |