High Content Image-Based Screening of a Protease Inhibitor Library Reveals Compounds Broadly Active against Rift Valley Fever Virus and Other Highly Pathogenic RNA Viruses

High content image-based screening was developed as an approach to test a protease inhibitor small molecule library for antiviral activity against Rift Valley fever virus (RVFV) and to determine their mechanism of action. RVFV is the causative agent of severe disease of humans and animals throughout Africa and the Arabian Peninsula. Of the 849 compounds screened, 34 compounds exhibited ≥50% inhibition against RVFV. All of the hit compounds could be classified into 4 distinct groups based on their unique chemical backbone. Some of the compounds also showed broad antiviral activity against several highly pathogenic RNA viruses including Ebola, Marburg, Venezuela equine encephalitis, and Lassa viruses. Four hit compounds (C795-0925, D011-2120, F694-1532 and G202-0362), which were most active against RVFV and showed broad-spectrum antiviral activity, were selected for further evaluation for their cytotoxicity, dose response profile, and mode of action using classical virological methods and high-content imaging analysis. Time-of-addition assays in RVFV infections suggested that D011-2120 and G202-0362 targeted virus egress, while C795-0925 and F694-1532 inhibited virus replication. We showed that D011-2120 exhibited its antiviral effects by blocking microtubule polymerization, thereby disrupting the Golgi complex and inhibiting viral trafficking to the plasma membrane during virus egress. While G202-0362 also affected virus egress, it appears to do so by a different mechanism, namely by blocking virus budding from the trans Golgi. F694-1532 inhibited viral replication, but also appeared to inhibit overall cellular gene expression. However, G202-0362 and C795-0925 did not alter any of the morphological features that we examined and thus may prove to be good candidates for antiviral drug development. Overall this work demonstrates that high-content image analysis can be used to screen chemical libraries for new antivirals and to determine their mechanism of action and any possible deleterious effects on host cellular biology.


Introduction
Many RNA viruses are highly pathogenic to humans and can cause hemorrhagic fever and/or encephalitis. Among these, Rift Valley fever virus (RVFV), a member of the genus Phlebovirus (family Bunyaviridae), causes Rift Valley fever (RVF) that severely affects human and livestock throughout Africa and the Arabian Peninsula [1,2]. In humans, RVF is characterized by a flu-like illness that occasionally progresses to hemorrhagic fever, encephalitis, or ocular disease with a significant death rate [3]. In ruminants, RVF is characterized by a high rate of abortion in pregnant females (80-100%), as well as high mortality in younger animals [4]. RVF is a mosquito-borne disease transmitted by the bite of various species of mosquitoes or through direct contact with infected animals [5,6]. Similar to RVFV, Lake Victoria marburgvirus (MARV), Zaire ebolavirus (ZEBOV), Lassa virus (LASV) and Venezuelan equine encephalomyelitis virus (VEEV) are highly pathogenic viruses and classified as National Institute of Allergy and Infectious Disease (NIAID) category A priority pathogens. Marburgviruses and ebolaviruses, members of the family Filoviridae, as well as LASV, members of the family Arenaviridae, cause severe hemorrhagic fevers in humans and nonhuman primates with extraordinarily high case-fatality rates (reviewed in [7,8]). VEEV, an alpha virus and member of Togaviridae family, causes severe encephalitis in horses and humans (reviewed in [9]). RVFV, as with other highly pathogenic RNA viruses, including EBOV, MARV, VEEV and LASV, cause severe disease in many developing countries that already suffer from fragile economies and health care infrastructures. There is currently no U.S. Food and Drug Administration (FDA) approved therapeutic or prophylactic treatments for any of these agents, thus there is an urgent need for research to develop effective new drugs and vaccines to combat these diseases.
Recent advancements in high content image (HCI)-based screening (HCS) technologies have contributed greatly to increasing the efficiency of the drug discovery process. HCS utilizes automated high-speed, high-resolution microscopy and image analysis to measure morphological changes in the cells in a quantitative and high-throughput manner [10]. Most importantly, HCI-based analysis enables simultaneous measurement of multiple features of cellular biology that are relevant to therapeutic and cytotoxic characteristics of potential antiviral compounds. As a result, HCI-based analyses not only allows for rapid screening of compounds, but can provide early insights into their cytotoxicity and mode of action, thereby facilitating the decision-making processes that govern the progression from a candidate compound to a successful antiviral drug.
RVFV is an enveloped spherical virus with containing a has a tri-segmented, single-stranded RNA genome, which encodes for the RNA-dependent RNA polymerase (RdRp), envelope glycoproteins (Gn/Gc), nucelocapsid proteins (N) and non-structural proteins including NSs and NSm [11]. Virus entry into cells is mediated by the binding of the envelope glycoproteins (Gn/Gc) to an unknown cell surface receptor which mediates virus endocytosis. Acidification of the virus-containing endocytotic vesicle promotes virus-host membrane fusion and results in the release of the encapsidated genome and RdRp into the cytoplasm, where transcription and replication of the viral genome occurs [12]. The glycoproteins Gn and Gc form a heteromeric complex and localize in steady-state to the Golgi apparatus due to a Golgi localization signal on Gn [6,10,13]. Multiple interactions between the glycoporteins, the encapsidated viral genome, and remaining structural proteins are believed to cause a change in membrane curvature leading to virus budding into the Golgi lumen [14][15][16]. Although the exact mechanism is not known, based on knowledge from other bunyaviruses, it is possible that the release of the virusfilled vesicles from the Golgi and its subsequent fusion with cell plasma membrane, releases mature virions into extracellular medium.
During the past two years, small molecule libraries containing compounds related to currently marketed drugs with known targets; mostly targeting kinases, phosphatases, or G protein coupled receptors; have been used with high throughput or HCIbased screening for antiviral drug discovery [17][18][19]. However these studies in which HCI analysis was used primarily focused on assay development and optimization specifically for determining the antiviral activity of the compounds and were not applied for the determination of the compound's mechanism of action [18,19]. Proteases constitute one of the largest drug target enzyme families, and their promise has been demonstrated by several successful drugs on the market, including the angiotensin converting enzyme inhibitors and anti-HIV protease inhibitors (reviewed in [20]). Many other drugs targeting key proteases involved in various diseases are predicted to reach the market in near future. Despite their success as drugs in the marketplace, to our knowledge, no studies using a protease inhibitor library and HCI-based analysis to screen for their antiviral activities have been previously reported.
Here we describe the development and application of HCS to test a protease inhibitor library of 849 small molecules for antiviral activity against multiple pathogenic RNA viruses and to determine their mechanism of action for RVFV.

Cells and Viruses
HeLa, Vero and Vero E6 cells were obtained from the American Type Culture Collection (ATCC) and maintained in Dulbecco's Modified Eagle Medium (Life Technologies, USA) supplemented with 10% fetal bovine serum (FBS, Life Technologies, USA). Human small airway epithelial cells (HSAECs) (Lonza, USA) were maintained under humidified conditions at 37uC, 5% CO 2 in Ham's F12 medium supplemented with nonessential amino acids, pyruvate, b-mercaptoethanol and 10% fetal calf serum (FCS). RVFV strains ZH548/MP-12 (MP-12) and ZH501 were obtained from The Salk Institute's Government Services Division and from Dr. Michael Turell (USAMRIID), respectively. RVFV strain MP12 was derived from a virulent strain of RVFV (ZH548) and is highly attenuated in its virulence due to several nucleotide mutations in its genome, however, it retains its sensitivity to IFN-a [16,21] and can be handled safety under biosafety level (BSL)-2 laboratory conditions. ZH501 is the wild-type strain of RVFV and is fully virulent, requiring to be handled under BSL-3 laboratory conditions. VEEV (1CSH3), LASV (Josiah), EBOV (Zaire) and MARV (Ci67) were obtained from the USAMRIID collection and were propagated in Vero cells. Virus-containing supernatant was clarified by centrifugation at 12,0006g for 30 min prior to storage at 280uC. All virus stock titers were determined by plaque assay on Vero E6 cells as previously described [22].

Virus Infections
Work with infectious filoviruses or arenaviruses were performed under BSL-4 laboratory conditions. RVFV ZH501 and VEEV infections were performed in a BSL-3 lab, and RVFV MP12 infections were performed under BSL-2 conditions. HeLa cells were infected with virus at the indicated MOIs. Inocula were removed after 1 h, unless stated otherwise, washed one time with 16 PBS, and replaced with the same amount of fresh medium. Infection was allowed to proceed for a specific duration of time, as indicated for each experiment. At the end of the incubation time, culture supernatants were collect to determine the virus yield by plaque assay. Alternatively, virus-infected cells were fixed in 10% neutral buffered Formalin (Sigma) for 3 days for BSL-3/4 viruses or 15 min for BSL-2 viruses.

Compound Library and Reference Reagents
ChemDiv was the source of the Serine/Cysteine Protease (SCP) inhibitor library used in this study. ChemDiv built the library using a variety of computational programs to qualify a small

Author Summary
Rift Valley fever (RVF) is an arthropod-borne viral zoonosis that occurs in large parts of sub-Saharan and North Africa and in 2000 emerged outside the African continent for the first time, raising concerns that it could further expand its geographical range. The disease in humans can result in encephalitis or hemorrhagic fever and in ruminants often results in abortion in pregnant females. Due to the lack of a licensed and commercially available vaccine, efforts to discover effective antiviral drugs are underway. Drug discovery using high content image-based screening is an effective tool that has been successfully used to identify new drugs. In this study, we developed an image-based assay to identify compounds active against RVF virus and other highly pathogenic human viruses. We demonstrated the usefulness of our image-based high content assay in identifying potential RVF antivirals by screening a small subset of chemical compounds for inhibition of RVF virus in a human cell line (HeLa) and partially characterized their mechanism of action within infected cells. The methods we developed in this study will be useful in discovering new effective drugs to combat Rift Valley fever.
molecule as a potential SCP ligand (reviewed in [23]). Basically, their method applied a set of pre-selected descriptors for encoding the molecular structures, and a trained neural network to qualify the molecules as potential SCP ligands. The molecular requirements were profiled by using available databases of active SCP and non-SCP-active agents. The library was composed of 849 compounds in a 96-well plate format and was used in both primary and secondary screening assays. For subsequent mechanism of action studies, compounds were reordered (ChemDiv, USA) at least three different times and tested independently to ensure consistency across batches. Reference reagents included human interferon-a2a (IFN-a) (PBL Interferon Source, USA), ribavirin (Sigma), Brefeldin A (cat # B-7450, Life Technologies) and Nocodazole (cat # M1404, Sigma). All compounds from the ChemDiv protease inhibitor library and ribavirin were prepared in 100% dimethyl sulfoxide (DMSO, Sigma-Aldrich, USA) at 8 mM stock solution and stored at 220uC. IFN-a was prepared in Dulbecco's phosphate-buffered saline (DPBS, Sigma-Aldrich, USA) containing 5% FBS.

Primary Screening of the ChemDiv Protease Inhibitor Library
The ChemDiv protease inhibitor library was screened against RVFV strain MP12 at 20 mM in 0.5% (v/v) DMSO. Screening was performed in triplicates wells of a 96-well plate. HeLa cells were treated with compound for 2 h prior to virus infection. After virus infection, cells were washed with PBS 16, and media containing the compound was added back to the cells and remained for the duration of the infection (24 h). To test the potential broad-spectrum antiviral activity of the primary hit compounds, HeLa cells were treated with the compounds in the same way as in the primary screen and infected with multiple RNA viruses as described below. The multiplicity of infection (MOI) and duration of infection for each virus were optimized to achieve an infection rate of 40-60% and were as follows: RVFV ZH501 (MOI = 1 for 24 h), EBOV (MOI = 3 for 48 h), MARV (MOI = 3 for 48 h), LASV (MO1 = 1 for 24 h), and VEEV (MOI = 0.5 for 20 h). Mock-infected HeLa cells were considered as negative controls, and wells in which the cells were infected with virus, but treated with 0.5% DMSO were considered as positive controls. Compound activity based on percent inhibition and cell toxicity was assessed by IFA and image analysis as described below. Data were plotted in Excel to generate a heat map using ''color scales'' in the conditional format menu. Scatter plot distribution was generated using TIBCO Spotfire 4.5.0 (TIBCO Software, USA).

Image Analysis
Confocal images were collected using a Leica TCS-SP5 confocal/multiphoton microscope. High-content quantitative imaging data were acquired and analyzed on an Opera confocal reader (model 3842 [Quadruple Excitation High Sensitivity] or model 5025; PerkinElmer) at two exposures using a 106 air objective or 406 water objective. Analyses of the images were accomplished within the Opera or Columbus environment using standard Acapella scripts.

Assay Validation and Data Normalization
The statistical reliability of the RVFV MP12 infection was evaluated by calculating for the Z9-factor using the formula: 12 [(3s p +3s n )/|m p 2m n |)], where m p , s p , and m n , s n are the mean (m) and standard deviations (s) of both positive (p) and negative (n) controls [20]. The mock-infected and RVFV-infected HeLa cells were used as negative and positive controls, respectively. Percentages of infected cells were normalized with mock-infected and RVFV-infected cells which were considered as 100%, and the values obtained were subtracted from 100 to determine the percent infection inhibition. The same formula was used for all percent infection inhibition calculations unless specified differently. For primary screening, Z9-factor $0.5 was used to validate the results of the assay. The hit (i.e., a compound that demonstrates inhibition of virus infection) selection criteria for the primary screening was set at $50% inhibition of RVFV infection.

Hit Confirmation by Dose-Response Curve Analysis
Dose response curve analysis was used to confirm hit compound activity against RVFV, VEEV, EBOV, MARV and LASV. Briefly, HeLa cells were seeded at a concentration of 2610 4 cells per well in a 96-well plate. After an overnight incubation at 37uC, hit compounds and references compounds, except for IFN-a, were tested in a 10-point dose-response curve (2-fold serial dilution from 200 mM) assay. IFN-a, was used as reference only in RVFV infections. The assay was performed in the same manner as the primary screening assay described above. Cells subjected to IFN-a were incubated with media containing serial dilutions (2-fold from 3000 units/mL) of IFN-a for 16 h (or overnight) prior to the start of virus infection. Two hours post-infection, cells were washed (one time with 16 PBS) and fresh medium was added back. At the end of virus incubation, cells were fixed and subjected to IFA to visualize the corresponding viral antigen expression as described in the secondary screen above. Each concentration of the hit compound was tested in triplicate. Percent infection inhibition was determined normalizing the percentage of infected cells in compound treated cells with mock treated (0.5% DMSO) and virus infected cells. Data generated from the image analysis of the dose response curve assay were plotted and analyzed using the non-linear regression formula: log (inhibitor) vs. response -variable slope (4 parameters) in GraphPad Prism 6. The EC 50 , defined as the effective concentration resulting in a 50% inhibition of virus infection, was used to evaluate compound activity. Compound toxicity was determined by normalizing cell number of compound treated+virus-infected cells with mock treated (0.5% DMSO)+ virus-infected cells, which were considered as 100%. The CC 50 value, defined as the compound concentration resulting in a 50% reduction in cell viability (based on normalized data) compared with mock infection was used to evaluate cell toxicity. The relative effectiveness of the compound is defined in terms of its selectivity index (SI), a value that indicates the relationship between the compound's effective and toxic concentrations, and is calculated as: SI = CC 50 /EC 50 . It is therefore desirable for a compound to have a high SI value, indicating maximum antiviral activity and minimal cell toxicity.

Time-of-Addition (ToA) Assay
HeLa cells at a concentration of 2610 4 cells per well in a 96well plate were either mock-infected or treated with compound 2 h prior to infection (pre), concurrent with virus infection (0 h), or at various time points post-infection as indicated in the experiment. Cells were infected with RVFV strain MP12 at an MOI of 10 for 1 h at 4uC. After the 1 h incubation, cells were washed one time with 16 PBS to remove unattached virus and incubated at 37uC with complete media for the duration of the experiment. After 13 h, cells were fixed in 10% Formalin and immunostained to detect viral G protein expression. Experiments were performed in triplicate and the average (6 standard deviation) of three independent experiments is shown.

Quantitative RT-PCR
RVFV RNA yields were determined by Quantitative RT-PCR (qRT-PCR) using TaqMan-based probe sequences as previously described [24]. In brief, total RNA from supernatants or cells were prepared using a MagMax 96 RNA extraction kit (Life Technologies) or RNeasy Plus minikit (Qiagen) respectively following manufacturer's instructions. Purified RNA was then measured with a Quant-IT RiboGreen RNA assay kit (Life Technologies) of which fifty nanograms of RNA was used in the assays. qRT-PCR assays were then performed on an ABI Prism 7900HT sequence detection system with an RNA UltraSense onestep kit (Life Technologies) and TaqMan probes (Applied Biosystems) according to the manufacturers' instructions. Serial 10-fold dilutions of the assayed virus (10 2 to 10 7 copies) were used as standards. The viral RNA levels from cell extracts was normalized with the PPIB housekeeping gene. Four wells of a 96 well plate were used to generate one data point and the experiment was done in triplicates. Relative expression levels were determined using the comparative threshold cycle (C T ) method [25].

Production of Stable Cell Line Expressing RVFV Gn-Gc Polyprotein
A HeLa cell line stably expressing RVFV Gn-Gc polyprotein, referred as HeLa-G was made by transfecting HeLa cells with the expression vector pCAGGS containing the open reading frame of RVFV Gn and Gc genes (kindly provided by Dr. Shinji Makino [26]) and an empty pcDNA3.1 vector expressing a neomycin selection marker at a 8:1 ratio followed by selection using Geneticin (G418; Life Technologies, USA), an analog of neomycin sulfate, at a concentration of 400 mg/mL. Isolated clones expressing RVFV Gn protein were verified by IFA and western blots analysis.
High Content-Based Image Analysis to Quantify Changes in Golgi Phenotype and to Evaluate Dose Response Curve Analysis RVFV-Gn localized to Golgi complex and serves as a Golgi marker. Thus HeLa-G cells were used to quantify changes in Golgi phenotype. HeLa-G cells were seeded at 1610 4 the day before the experiment. Cells were mock treated or treated with hit compounds for 12 h, with the exception of IFN-a2a, in which cells were treated overnight prior to start of the experiment. Following treatment, cells were fixed in 10% Formalin, permeabilized with 0.1% Triton-X-100 and immunostained to detect Gn protein expression. Subsequently cells were also stained with Hoechst for nuclear detection. Images were acquired on Opera using 206 air objective, and the images were analyzed with Columbus High Content Imaging software. Individual cells were segmented based on the Hoechst nuclear stain using the ''Find Nuclei'' building block and using ''Select Cell Region'' building block, a ring region of 15 pixel width was then created around the individual nuclei that was further masked to classify them as cells. A common boundary was drawn connecting the individual Golgi stacks appearing as spots and spaced less than 2 pixels apart to create a ''Golgi mask.'' Subsequently, the area of the Golgi mask was evaluated and only those with surface area $100 pixels were considered as ''intact Golgi.'' Percent Golgi was determined by normalizing ''Golgi number'' per field with the ''number of nuclei'' in the same field, and the resulting values were then normalized with vehicle control (0.5% DMSO) treated cells. The values obtained were plotted and analyzed using the non-linear regression formula: log (inhibitor) vs. response-variable slope (4 parameters) in GraphPad Prism 6. The GC 50 was defined as the effective concentration resulting in a 50% reduction in Golgi numbers.

High-Content Image Analysis of Tubulin
HeLa cells were either mock treated (0.5% DMSO) or treated with BFA (50 ng/mL) for 15 min or nocodazole (33 mM) for 5 h or D011-2120 (50 mM) for 4 h and fixed in 10% Formalin, permeabilized with 0.1% Triton-X-100, and immunostained with tubulin antibody. Subsequently, cells were stained with nuclear Hoechst (33342) dye and Deep red (Cell mask). Images were acquired on an Opera System using a 406 air objective and exported to Columbus High Content Imaging and Analysis Software. Individual cells were segmented based on the Hoechst nuclear stain using the ''Find Nuclei'' building block and the cell mask deep red stain using ''Select Cell Region.'' To quantify tubulin phenotype, texture analysis was performed using the Columbus software and the ''SER Features'' method. The ''SER Features'' method analyzes the occurrence of typical patters in the intensity structure of an image. The ''Ridge Value'' was selected as this feature showed the most difference in the intensity pattern of mock-treated vs. nocodazole-treated cells. The ridge values were collected from 4 fields and 3 wells for each condition, and the values were normalized with control treated cells, which were considered as 100%.

High Content Image-Based Assay Development and Validation
HeLa cells were considered the optimal cell line to use for assay development as the cells are permissive to RVFV infection, are the optimal size for image analysis, possess a functional interferon response, and they are relatively easy to grow. Selection criteria for the appropriate cell seeding density included having a sufficiently high number of cells but with enough spatial distribution for proper identification during image acquisition. After testing various seeding densities of HeLa cells, the seeding density of 2610 4 cells per well was selected as minimal variations among replicates were noted without compromising proper segmentation of cells (data not shown). RVFV infection in HeLa cells was visualized by IFA detection of the RVFV Gn protein using the 4D4 mAb and confocal imaging using the Evotec Technologies High-Throughput Cell Analyzer Opera (Fig. 1A). Expression of this protein on the cell surface is indicative of virus egress from the infected cell. For the RVFV infection, a MOI of 1.0 and incubation time of 24 h was used since it achieved an infection rate of ,60% and allowed multiple rounds of virus replication that facilitated the screening of active compounds that target different stages of virus life cycle. Finally, the RVFV infected cells were classified as those cells having G signal within the defined boundaries of a cell and greater than a pre-defined threshold level. Several additional parameters (,50; data not shown) were acquired from the images, some of which include the nuclear and cytoplasmic intensities of the G signal, nuclear size and nuclear intensities. However, for the screening purposes described in this section, only cell number and the infected cell number was used to determine the percentage of infected cells.
The Z9-factor was evaluated to test the robustness of the assay in the 96-well plate format. Figure 1B shows a representation of the validation process done for the RVFV infection in HeLa cells. Cells were dispensed in a 96-well plate and half (48 wells) were mock-infected (0.5% DMSO only) while the other half were infected with RVFV. Mock-infected wells were considered as negative controls, while RVFV-infected wells as positive controls.
After the viral infection period and IFA, image analysis of the acquired images revealed the infection percentage, cell number based on cell nuclei detected, and other pre-defined parameters. The Z9-factor was calculated using average and standard deviations of the percent infection of the positive and negative controls as described in methods section [27]. The experiment was performed in triplicate on three separate days, and the calculated Z9-factor was 0.62. A Z9-factor $0.5 indicates a statistically reliable separation between positive and negative controls.

Screening of the ChemDiv Protease Inhibitor Library
A total of 849 compounds from ChemDiv protease inhibitor library were used for antiviral screening. The library was screened at 20 mM concentration against RVFV and primary hits were selected based on the criteria: $50% inhibition and $50% cell number compared to positive controls. Figure 2 shows a scatter plot of the primary screening data where percent infected cells were normalized with mock-treated RVFV-infected cells to determine the percent infection inhibition (see methods). Thirteen separate master compound plates were used and are represented as separate colors in Figure 2. Many of the hit compounds clustered together on the scatter plot since the compounds were grouped together on the plate by primary chemical structure. Sixty-one compounds were initially selected as primary hits (exhibiting $50% inhibition), although we suspected some of these may be due to edge well effects. Upon repeat screening, the compounds with suspected edge well effects proved not to be active and were eliminated from subsequent analysis. Thus, a total of 34 compounds (4.0%) exhibited $50% inhibition against RVFV (shown above the black horizontal line at 50% infection inhibition in Figure 2).
In order to determine if the primary hits had broad-spectrum activity against other highly pathogenic RNA viruses, they were further screened against multiple viruses, including EBOV, MARV, VEEV, LASV, and the virulent RVFV ZH501 strain, in addition to MP12. Figure 3A shows the generated heat map of percent virus infection inhibition and percent cell viability. As can be seen in the heat map, most of the compounds had varying degrees of activity across the multiple RNA viruses tested. The 34 primary hit compounds were classified into four distinct groups based on their primary chemical structure as shown in brackets in Figure 3A. Four compounds, each comprising a unique primary    (Fig. 3B), were selected based on activity against three or more of the viruses tested (Fig. 3A). All other primary hit compounds were chemical analogs of these four compounds. These compounds were selected for further study to determine their EC 50 , CC 50 , SI, and mechanism of action against RVFV.
Antiviral activities of the four hit compounds against RVFV were determined by dose response curve analysis in HeLa cells. The EC 50 , CC 50 and SI were determined for the four compounds and two reference compounds, ribavirin and IFN-a2a, known inhibitors of RVFV. Figure 3C shows the dose response curve analysis of the four hit compounds and two reference compounds, with percent infection normalized against RVFV-infected HeLa cells. All compounds tested showed inhibition of RVFV-infected HeLa cells in a dose-dependent manner (Fig. 3C). Furthermore, all compounds tested had a SI of greater than 1, indicating that the virus inhibitory effects were apparent before the cytotoxic effects appeared. G202-0362 had the lowest effective dose (EC 50 = 13.6 mM), and significantly, appeared to have no toxic effects on the HeLa cells at the concentration used in the experiment (Table 1). C795-0925 had a similar effective dose (EC 50 = 14.4 mM), albeit the compound was more toxic to cells (CC 50 = 56.7 mM) compared to G202-0362. D011-2120 had the highest effective dose (EC 50 = 61.5 mM) and consequently the lowest SI (1.5) ( Table 1). The effective doses of ribavirin and IFNa inhibited RVFV were consistent with what has been reported in the literature [17,28,29]. In order to examine the potential that the inhibitory effects of these compounds were due to an artifact of using a cancer cell line (HeLa), we repeated the screens of these four compounds and references compounds against RVFVinfected human small airway epithelial (SAEC) cells, a primary cell line derived normal human lung tissue. The effective doses of all the compounds, except for IFN-a2a, were similar to those observed in HeLa cells (Table 1). IFN-a was approximately 10fold less effective against RVFV infection in SAEC compared to in HeLa cells. IFN-a has several targets in mammalian cells and therefore, although highly effective, is not an ideal model to compare with other drugs. Interestingly for C795-0925, while the EC 50 was comparable between the two cell lines, the compound was less toxic to SAEC (CC 50 = 110.6 mM).
The hit compounds' potency against other pathogenic viruses including VEEV, EBOV, MARV and LASV were further evaluated by does response analysis (Table 2) and compared to ribavirin's antiviral activity. In VEEV infections, all 4 compounds were more potent than ribavirin (a broad spectrum antiviral agent) while C795-0925 was most potent among all with EC 50 = 10.2 mM. G202-0362 was most potent in filoviruses EBOV and MARV with EC 50 = 16 mM and 34 mM respectively and more potent than ribavirin. However, in LASV infections, ribavirin was most effective compared to all 4 compounds. Collectively, these data show that these compounds were effective against more than one virus.
Hit Compounds Target Viral Replication and/or Viral Egress, but Not Entry RVFV infection follows an ordered sequence of steps that can be visualized using IFA for the early nucleocapsid (N) and late glycoprotein (Gn) expression in infected cells [28]. Compounds targeting any stage of the viral life cycle should alter either the kinetics or localization of viral protein expression. Thus we used time-of-addition (ToA) assays, adding compounds at various time points before or after virus infection, combined with HCI analysis to determine at which stage of the viral life cycle the compound is acting. To examine the kinetics of N, Gn, and cell surface exposed Gn (referred to as S-Gn) expression, HeLa cells were infected with RVFV MP12 virus (time = 0 h), and the percent of cells expressing N, Gn, and S-Gn were determined by IFA at 4, 6, 8, 10, and 12 h post-infection (PI). N protein expression was apparent at 4 h PI and steadily increased throughout the 12 h infection cycle, at which time approximately 70% of the cells expressed the protein (Fig. 4A). Gn expression was first apparent at 8 h PI and steadily increased throughout the infection cycle to a peak of approximately 40% of the cells expression the protein at 12 h PI. Cell surface exposed Gn (S-Gn) protein expression is indicative of virus egress. Consistent with this, S-Gn was first apparent at 10 h PI and peaked at 12 h PI in approximately 25% of the cells (Fig. 4A) and corresponds to one round of replication of the RVFV life cycle [28]. Release of virus into the extracellular medium was further confirmed by quantitative PCR analysis (Q-PCR) and plaque assay (data not shown). Compounds targeting virus entry or virus and their primary structures are shown (B). (C) Dose response curve analysis of the four lead compounds and the reference compounds ribavirin and IFN-a2a. Percent infection inhibition is shown in blue and the percent cell number is shown in red. doi:10.1371/journal.pntd.0003095.g003 replication will show aberrant kinetics in the cytoplasmic expression of N and Gn. Alternatively, if the compound only targets virus egress then changes in the kinetics of S-Gn expression, and not cytoplasmic Gn or N expression, should be observed. Therefore, the kinetics of viral N, Gn, and S-Gn expression in the presence of compound was examined. The percentage of cells expressing N, Gn or S-Gn upon compound treatment was normalized with mock treated and infected cells, which was considered as 100% (Fig. 4B). Treatment with all compounds significantly reduced the expression of S-Gn by 70-100%, indicative reduction of infectious virus getting to the cell surface. However, treatment with F694-1532 and C795-0925 also significantly reduced expression of N and Gn, suggesting these compounds exhibited their effects earlier in the virus life cycle. While treatment with D011-2120 and G202-0362 reduced S-Gn expression approximately 90%, they had only minor effect (,10% reduction) on N and Gn expression (Fig. 4B). These data suggest that D011-2120 and G202-0362 are targeting virus egress, while F695-1532 and C795-0925 may be targeting steps earlier in the virus life cycle.
To further define the stage of the viral life cycle the compounds are acting, compounds were added either 2 h prior to infection (2 2 h), at the time of infection (0 h), or 2 h post-infection (2 h). Acid conditions (pH,5.0) inside the endocytosed virus-containing vesicles are necessary to promote virus-host membrane fusion [12,28]. Thus, we used NH 4 Cl (20 mM) as a control as it targets the last stage of the virus entry pathway by neutralizing the pH inside the acidic compartments within the cell. Figure 4C shows that NH 4 Cl inhibited RVFV infection when cells were treated at early time points (22 h or 0 h), but not at the later time point (2 h), consistent with its role as an entry inhibitor. Conversely, levels of virus inhibition with each of the four hit compounds added 2 h PI were similar to levels of inhibition when the compounds were added at 22 h or 0 h. (Fig. 4C). These data suggest that the observed inhibitory effects of the four hit compounds were not due to a block in the virus entry pathway.
In order to examine how the compounds affect subsequent stages of the virus life cycle, compounds were added to virusinfected HeLa cells at 2, 4, 6, 8, 9, 10, 11, 12, or 12.5 h PI and the percentage of infected cells were evaluated at 13 h PI using IFA for S-Gn expression, which is indicative of virus egress.
More than 50% of infection was rescued when F694-1532 or C795-0925 were added at 8 h or 9 h PI respectively (Fig. 4D, marked with red arrows), which suggested that compounds mostly targeted viral RNA replication or viral protein expression. Not surprisingly, D011-2120 and G202-0362, which appeared to target virus egress, did not recover restriction until late stages (. 10 h for 50% infection recovery). Collectively these data suggested that C795-0925 and F694-1532 targeted viral RNA transcription or protein synthesis and not virus entry or egress, while D011-2120 and G202-0362 targeted egress. To further confirm these observations, viral RNA levels were evaluated to determine compound mediated effects on RVFV life cycle as described below.
The compound mediated effect on virus replication was determined by measuring the amount of viral RNA that accumulated inside infected cells or in their supernatants at the end of one life cycle. The viral RNA in the supernatants is a  measure of virus egress. Cells were mock-treated or incubated with any of the four hit compounds beginning 2 hours prior to the start of infection. Twelve hours post-infection, the cell extracts or supernatants were collected separately, and the viral RNA levels were evaluated by qRT-PCR. The RNA levels in compoundtreated cells were normalized with mock-treated cells to calculate the relative change in RNA levels (Fig. 4E). Alternatively, in the duplicate experiment, cells were subjected to IFA of Gn and S-Gn protein expression. HCI was then applied to evaluate the percentage of cells expressing Gn or S-Gn, and the relative change was calculated. (Fig. 4E). As expected, D011-2120 and G202-0362 treatments that targeted virus egress resulted in low levels of 9% and 28%, respectively of viral RNA accumulation in the supernatants as compared to 75% and 95%, respectively inside the cells. In contrast, viral RNA levels inside the cells and in the supernatants were similar for cells treated with C795-025 and F694-1532. In addition, these data corroborated with the viral protein expression levels inside cells (Gn) or on the surface of cells (S-Gn). Collectively, these data suggest that C795-0925 and F694-1532 targeted viral RNA transcription or protein synthesis, and not virus entry or egress, while D011-2120 and G202-0362 targeted egress. HCI was further applied to understand the mode of action as described below.

D011-2120 Inhibits RVFV Glycoprotein Gn Localization by Disrupting Golgi Complex
RVFV envelope glycoproteins localize to the Golgi complex and mediate virion assembly and budding into the Golgi lumen [13,14,30]. Subsequently, based on similarities with other bunyaviruses, it is believed that virus egress occurs by the fusion of the virion-containing Golgi vesicles with the cellular plasma membrane. Because D011-2120 and G202-0362 appeared to target virus egress, HCA was applied to determine if the virus trafficking from the Golgi complex to plasma membrane was affected. Golgi complex was visualized by transducing cells overnight with baculovirus expressing fusion construct of a Golgi marker and green florescent protein (GFP). RVFV G expression in cytoplasm is detectable by 10 h PI (Fig. 4A) and therefore this time point was used to interrogate any potential compound-mediated changes in G localization within the Golgi complex of infected cells. At 10 h PI, cells treated with 0.5% DMSO alone displayed a bright G signal localized to one side of the nucleus, which is consistent with the location of the Golgi complex (Fig. 5A-B). A portion of the G protein co-localized with Golgi as evidenced by the bright yellow staining when the images were overlaid (Fig. 5C). Furthermore, numerous G stained spots (Fig. 5B-C), presumably virus containing Golgi vesicles were found in the cytoplasm (marked with white arrow). G202-0362 treated cells, similar to DMSO treated cells, expressed G that co-localized with Golgi apparatus (Fig. 5D-F). However the numerous G stained spots in the cytoplasm that were observed in DMSO treated cells were not observed. Therefore it is possible that G202-0362 is targeting virus assembly, virus budding into Golgi lumen or the subsequent release of the virus containing Golgi vesicle into the cytoplasm.
On the contrary, D011-2120 treated virus-infected cells resulted in Golgi dispersion (Fig. 5G-I). G remained localized with the disrupted Golgi. The effects of compounds on Golgi and G protein localization were compared to that of two well-known Golgi-disrupting compounds, nocodazole and brefeldin A (BFA). Nocodazole is a microtubule-disrupting agent and induces a slow dispersal of the juxtanuclear Golgi to peripheral sites [31] (Fig. 5J-L), while BFA causes a rapid disassembly and redistribution of the Golgi into the endoplasmic reticulum (ER) [32], (Fig. 5M-O). D011-2120 and nocodazole had same effect on Golgi dispersion (compare Fig. 5I and 5L) suggesting a common mechanism of action. These results suggest that D011-2120, exerts its antiviral activity by disrupting the Golgi complex in a manner similar to nocodazole, but distinct from BFA.

High-Content Image-Based Analysis of Changes in Golgi Phenotype by the Hit Compounds
We further developed HCI analysis to quantify alterations in Golgi structures, to estimate GC 50 values (effective concentration of compound that disrupts 50% of the Golgi complexes) of D011-2120 targeting the Golgi, and to screen other compounds that may also target the Golgi. HCI-based analysis was developed to quantify changes in Golgi phenotype. For this, a HeLa cell line that stably expresses RVFV Gn was generated (HeLa-G). In the absence of viral RNP complex, Gn remains in the Golgi ( [14] and reviewed in [33]) and therefore IFA of Gn served as a Golgi marker. IFA showed specificity in the localization of RVFV Gn with Golgi marker (Fig. 6A). Images of Gn stained HeLa-G that were mock-treated or treated with D011-2120 were used further for the development of HCI-based assay for cellular Golgi distribution. Figure 6B shows the individual steps that allowed successful segmentation of nuclei (I), cells (II) and the Golgi complex (III). The image analysis was applied such that the intact juxtanuclear Golgi appeared as one ''intact Golgi'' (see methods). Figure 6B-IV (see insert) shows a red boundary generated by image analysis to classify ''intact Golgi.'' D011-2120 treated cells showed Golgi staining as numerous spots dispersed throughout the cytoplasm, which was disqualified as 'Intact Golgi'' by image analysis program (Fig. 6B-V, see insert). This image analysis was further applied to determine how the hit compounds modulated Golgi phenotype in HeLa-G cells that were treated with increasing concentrations of the hit compounds or the reference compounds Ribavirin or IFN-a (Fig. 6C). Several features of the Golgi complex were generated, including ''Golgi number'' and ''Golgi area'' per well, and a heat map thus generated shows a clear reduction in Golgi ''number'' or ''area'' in wells with higher concentration of D011-2120. The ''Golgi number'' parameter was further used to determine the GC 50 values by dose response curve analysis (Fig. 6D), which were compared to the compounds' EC 50 as shown in Table 3. The GC 50 value for D011-2120 was less than half that of its EC 50 value. The ratio of GC 50 to EC 50 was further calculated and was determined to be 0.62. A ratio of less than 1 suggests that the compound's antiviral activity is related to its ability to disrupt the Golgi complex. Other compounds had lower GC 50 values, but when compared to their EC 50 , none of their GC50/EC50 ratios were below 1 (Table 3). Surprisingly, F694-1532 treated cells also resulted in a decrease in Golgi numbers. However, further analysis showed that this compound inhibited expression of Gn or any other exogenously introduced gene such as GFP, but did not disrupt Golgi phenotype (data not shown). level or in the supernatants of compound treated cells at 12 h PI. Cells were mock treated or treated with the 4 hit compounds from 2 h prior to the start of infection through the entire duration of infection. In a duplicate experiment the percentage of cells expressing Gn and S-Gn was determined by HCI of the immunostained cells. The relative fold change was determined by normalizing data of compound treated cells with values of corresponding mock treated cells. doi:10.1371/journal.pntd.0003095.g004

D011-2120 Acts by Depolymerizing Microtubules Similar to Nocodazole
Because of the similarities observed in Golgi staining patterns in cells treated with D011-2120 and cells treated with nocodazole, we hypothesized that they are both acting by a similar biochemical mechanism. It is well-known that nocodazole exerts its effects by inhibiting polymerization of microtubules [34], which are made up of linear polymers of tubulin that undergo continual assembly and disassembly to form a dynamic filamentous microtubule network within cells. Thus, we reasoned that D011-2120 may also be acting to inhibit this process. To test this, the kinetics of Golgi complex dispersion in HeLa cells exposed to either vehicle control (0.5% DMSO), D011-2120, nocodazole, or BFA were examined by live cell imaging. Our previous data indicated that a 4 h exposure to both D011-2120 and nocodazole was optimal for cells to exhibit the Golgi dispersion phenotype (data not shown). Thus, HeLa cells infected with MP12 virus for one life cycle (12 h) and treated with D011-2120 or Nocodazole during the last 4 h of infection or with BFA during the last 15 min of infection. Tubulin was visualized by IFA. As expected, treatment with either the vehicle control (0.5% DMSO) or BFA showed predominantly filamentous microtubule network (Fig. 7A). In contrast, cells treated with either nocodazole or D011-2120 showed diffused pattern of tubulin staining (Fig. 7A). To further quantify tubulin phenotypic changes, HeLa cells in a 96-well plate were treated with D011-2120 or reference compounds for 4 h and immunostained with tubulin antibody. The images were further subjected to HCA to quantify the changes in the pattern of microtubule organization using the ''Ridge Value'' of the patternanalysis program (Columbus Software). The ''ridge'' feature best distinguished the filamentous pattern in control cells from the diffused staining of tubulin in the compound treated cells (Fig. 7B). The ridge value for cells treated with D011-2120 was similar to cells treated with nocodazole, which were approximately half that of cells treated with BFA or mock-treated (Fig. 7B).
Because the separation of chromosomes requires functional spindle fibers that are made of microtubules, microtubule depolymerizing agents, such as nocodazole, are known to arrest cells in the G2/M phase of the cell cycle [35]. To determine if D011-2120 induces mitotic arrest in treated HeLa cells, cells were treated with various concentrations of the compound and immunostained to detect the expression of phosphorylated histone 3 at serine number 10 (H3-pS10), a known marker of cells undergoing mitosis and that is associated with chromatin condensation [36]. As compared to mock-treated (0.5% DMSO) cells, there was a significant increase in the number of cells expressing H3-pS10 when cells were treated for 24 h with 25 mM of D011-2120 or 3.3 mM nocodazole (Fig. 8A). Furthermore, the effect was dose dependent with the number of cells expressing H3-pS10 increasing with increasing concentration of D011-2120 up to 25 mM (Fig. 8B). Concentrations above this amount were toxic to the cells. Taken together, these data suggest that D011-2120 acts by depolymerizing microtubules resulting in mitotic arrest.

Discussion
RVFV, along with VEEV, MARV, EBOV and LASV all cause serious disease in humans and animals (RVFV and VEEV). As such, they are classified by NIAID as category A priority pathogens with bioterrorism potential. In addition to being highly virulent, there is currently no FDA approved antivirals to treat any of them. In this study, we applied HCI-based analysis as a screening tool to discover novel compounds with antiviral activity against these highly pathogenic viruses and to elucidate their mode of action. Image-based high-throughput screening of potential antiviral compounds is an extremely effective drug discovery tool as it can not only determine antiviral activity, but can also measure many aspects of the viral life cycle and virus interaction with multiple host cellular components. However, many potential antiviral compounds could be targeting a cellular component that is essential for normal functioning of host cells. It is at this juncture that HCI-based screening has promise over traditional highthroughput screening assays. The sophisticated image analysis in HCS offers more information including intensity of expression, localization and pattern of expression (i.e., texture) of the fluorescent target. Since viral infection follows a specific sequence of viral gene expression, any changes in the kinetics of viral protein expression patterns can provide useful insights into the stage of the virus life cycle targeted by the compound. This information can, therefore, be helpful in forming a hypothesis regarding which specific virus-host protein interaction(s) the compounds are acting on.
To discover novel compounds with broad spectrum antiviral activity, the strain MP12 of RVFV served as a model system for primary screening. This attenuated strain was previously shown to faithfully recapitulate the replication properties of wildtype RVFV [15] and unlike the wildtype virus, can be handled safely under BSL-2 laboratory conditions. Thus, we used MP12 virus-infected cells to set the optimal system parameters, which included cell type, cell number, duration and multiplicity of infection, to achieve a statistically reliable high-throughput screening assay with a Z9factor of 0.62 (Fig. 1B). Basic features such as cell number and infected cell number were used to determine the ''percentage of infected cells.'' Several additional parameters such as size of cell and nucleus were also acquired that provided additional information to determine if the compound were toxic to cells was acquired (data not shown) but were not used in this study. Using these optimized conditions, a focused small molecule library of 840 serine/cysteine protease inhibitors from ChemDiv were screened against MP12 RVFV infection in HeLa cells.
Primary screening of the protease inhibitor library yielded 34 primary hits (Fig. 2), which could all be classified as one of four chemical structures (Fig. 3B). Interestingly, none of these four chemical scaffolds had been previously reported to possess antiviral activity. The 34 primary hits were further screened against the fully virulent RVFV ZH501 strain, VEEV, LASV, EBOV and MARV. These viruses represent four different viral families and thus would give us some indication if the hit compounds were broadly active. Many of the hit compounds were broadly active against all the viruses tested, albeit with different efficiencies. For example, compound C795-0925 had the highest activity against RVFV and VEEV, but exhibited moderate antiviral activity against EBOV, MARV, or LASV. Similarly, compound D011-2120 had high activity against EBOV, MARV, and VEEV, but moderate activity against RVFV or LASV (Fig. 3A). One compound from each chemical scaffold that showed broad spectrum antiviral activity (active against at least three of the viruses tested) were chosen for further study using a dose response curve assay against RVFV. Three of the four hit compounds had EC 50 values that were below that of ribavirin, a drug commonly used to treat infections due to RNA viruses. To examine whether the compound's activities were specific to HeLa, we also tested them in the primary cell line SAEC. All of the hit compounds had activities similar in these cells as compared to those seen in HeLa cells; however, compound C795-0925 was less toxic in SAEC (CC 50 56.7 mM in HeLa compared to 110.6 mM in SAEC).
The virus life cycle is an orchestrated series of specific steps that include virus entry, viral RNA replication, transcription, followed by viral protein expression, and finally the assembly of the viral structural proteins and egress out of the cell. To determine the specific stage of virus life cycle targeted by the compounds, we first evaluated the kinetics of RVFV N, Gn and S-Gn expression (Fig. 4A) [28]. Cells treated with C795-0925and F694-1532 exhibited decreased levels of expression of the viral proteins N and Gn. In contrast, cells treated with D011-2120 and G202-0362  Table 3. Golgi disruption properties of the four hit compounds and reference compounds derived from HCI-based analysis of HeLa-G cells. showed normal expression levels of N and Gn, similar to that seen in mock-treated virus infected cells. As expected, cells treated with C795-0925and F694-1532 also exhibited significantly reduced levels of surface expressed Gn (S-Gn), as fewer virions were being made during treatment with these compounds, thus fewer virions were getting to the cell surface. However, with compounds D011-2120 and G202-0362, cells exhibited normal levels of N and Gn, but greatly reduced levels of S-Gn, suggesting the inhibition was at a more downstream point in the viral life cycle, such as virus assembly or egress. Thus, we used time-of-addition assays to determine at which stage these compounds were acting. Since none of the compounds inhibited virus infection when added at the time of infection or 2 h prior to infection, we concluded that the compounds are not having an effect on virus entry (Fig. 4C). When time of compound addition was extended to 12.5 h PI, results suggest that D011-2120 and G202-0362 were affecting virus egress; whereas, C795-0925and F694-1532 were affecting viral RNA replication. We were further interested in using HCI-based analysis to define the mechanism of action used by D011-2120 and G202-0362 to inhibit virus egress. RVFV glycoproteins localize to the Golgi, promote virus assembly, and bud into the Golgi lumen [14,33]. Based on studies with other bunyaviruses, it is speculated that virus egress results from the budding of these virus-containing vesicles from the Golgi, trafficking through the cytoplasm, and their subsequent fusion with the plasma membrane, releasing the virus into the extracellular environment [15,30,33]. This is consistent with our data showing strong Gn staining in the Golgi at early time points (data not shown) and diffuse punctate stain throughout the cytoplasm at later time points (i.e., 10 h PI) (Fig. 5B). Presumably this represents RVFV-containing vesicles released from the Golgi and making their way to the cell surface as has been described for other bunyaviruses [22,33]. These vesicles were not observed in cells stably expressing Gn alone (Fig. 6A), further supporting this idea. Given our observations, we wanted to determine if either G202-0362 or D011-2120 (compounds earlier shown to effect Figure 7. D011-2120 depolymerizes microtubules similar to nocodazole. (A) IFA of HeLa cells that were mock infected of infected with MP12 virus for 12 h and either were mock treated or treated with vehicle control (0.5% DMSO), Brefeldin A (100 ng/ml) nocodazole (33 mM), and D011-2120 (50 mM) for 15 min, 4 h or 10 h respectively prior to stopping infection. Cells were immunostained to detect tubulin (green) or RVFV Gn (red). (B) HeLa cells in a 96 well plate were treated with D011-2120 or reference compounds for 4 h and immunostained with tubulin antibody. Acquired images were subjected to pattern and ridge analysis (Columbus software) to evaluate differences in microtubulin organization within the cells. Each data point is derived from eight replicates. Percent ridge values were obtained by normalizing each data point with 0.5% DMSO treated cells, which was considered 100%. doi:10.1371/journal.pntd.0003095.g007 egress) affected viral trafficking from the Golgi to the cell surface. Interestingly, cells treated with G202-0362 showed Gn localized to the Golgi, but no virus could be seen in vesicles in the cytoplasm, suggesting that this compound was blocking the virus budding from Golgi and subsequent trafficking to the cell surface (Fig. 5). In contrast, D011-2120 appeared to completely disrupt the Golgi in a manner similar to nocodazole (i.e., by depolymerizing microtubules), as evidenced by co-localization of the Golgi marker, 1,4-galactosyltransferase, and Gn in a staining pattern similar to that seen with nocodazole ( Fig. 5G-L). These data suggested that D011-2120, but not G202-0362, may be acting by disrupting cellular microtubules. Because we were interested in further defining the phenotypic alterations in the Golgi caused by treatment with D011-2120, we developed a HCI-based assay for this purpose. An algorithm using Columbus software was built that allowed classification of normal juxtanuclear Golgi complex from stacks of Golgi dispersed throughout the cytoplasm (Fig. 5B). The algorithm was applied to determine the potency of the compounds to determine the concentration of compound that gives 50% reduction in Golgi numbers (termed GC 50 ) ( Table 3). When these values were compared to EC 50 values, it was evident that D011-2021 antiviral activity was related to its Golgi disruptive properties. Its GC 50 was 38.2 mM, which was significantly lower than its EC 50 (61.5 mM). Interestingly, while F694-1532 exhibited a relatively low GC 50 (35.3 mM), it was still higher than its EC 50 (20 mM). Likewise, as expected, F694-1532 and G202-0362 GC 50 values were higher than their corresponding EC 50 values, confirming that they also do not possess Golgi disruptive properties. We also examined the %GN L , which is the lowest percentage of Golgi number achievable at the compound's highest concentration. D011-2120 had the lowest %GN L further suggesting this compound was an active Golgi disruptor. Because the staining pattern of cells treated with D011-2120 was similar to that seen with nocodazole treatment (Fig. 5), we hypothesized that D011-2120 was acting by a mechanism similar to that of nocodazole (i.e., depolymerization of microtubules). Because spindle fibers, which are made of microtubules, are an essential component involved in mitosis, and subsequently, cell division, a consequence of microtubule depolymerization is mitotic arrest. Consistent with the idea that D011-2120 exerts its effect by microtubule depolymerization, we showed that cells treated with this compound undergo mitotic arrest (Fig. 8A).
Collectively, our data suggest that the antiviral action of D011-2120 is by depolymerizing microtubules which disrupted virus trafficking of virions from the Golgi to the cell membrane during virus egress in RVFV infection. Microtubulin plays an important role in the life cycle of many viruses supporting virus entry, RNA replication, assembly and or egress (reviewed in [37])). The specific stage of the virus life cycle that is critically dependent on microtubulin defers between viruses and therefore it is possible that D011-2120 targeted virus entry or replication but not egress of other RNA viruses. Moreover microtubulin depolymerization induces mitotic arrest, and reduces cell viability in cells which can impair virus growth. This is apparent in the SI values of #1.2 derived from D011-2120 treatments in VEEV, EBOV, MARV or LASV infections, which meant that the compounds antiviral action could be due to compromised cell viability.
G202-0362 also affected virus egress, but it appears to do so by a different mechanism to that of D011-2120, namely by blocking virus budding from the trans Golgi. F694-1532 inhibits viral replication by an as yet unknown mechanism. This compound, however, also appeared to inhibit overall cellular gene expression and formed aggregates in the cell. Thus, while D011-2120 and F694-1532 may prove to have utility in treating other diseases, such as certain types of cancer, they are not attractive candidates for antiviral therapeutics. G202-0362 blocked virus budding from the Golgi, but did not disrupt the Golgi and did not alter any of the morphological features that we examined, including changes in cytoskeleton organization, cell or nucleus size, or Golgi or endoplasmic reticulum structure (data not shown). Moreover, this compound had no observable effects on cell division in uninfected cells at the highest concentration tested (200 mM). In addition this compound was effective against filoviruses and VEEV Likewise C795-0925 was highly effective against VEEV with EC 50 = 10.2 mM and SI = 6.2. Thus, these two compounds may prove to be good candidates for antiviral drug development. In this study, we showed how HCI-based analysis could be used to screen for antiviral compounds effective against RVFV or other highly pathogenic RNA viruses and determine their mechanism of action.