Next Article in Journal
Inflammation and Epstein–Barr Virus at the Crossroads of Multiple Sclerosis and Post-Acute Sequelae of COVID-19 Infection
Previous Article in Journal
Use of Remdesivir in Patients with SARS-CoV-2 Pneumonia in a Real-Life Setting during the Second and Third COVID-19 Epidemic Waves
Previous Article in Special Issue
Antiviral RNAi Response in Culex quinquefasciatus-Derived HSU Cells
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 15
Issue 4
10.3390/v15040948
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Antiviral Activity of Acetylsalicylic Acid against Bunyamwera Virus in Cell Culture

by
Sara Yolanda Fernández-Sánchez
1,
José P. Cerón-Carrasco
2,
Cristina Risco
1,* and
Isabel Fernández de Castro
1,*
1
Cell Structure Laboratory, Centro Nacional de Biotecnología, CSIC, Campus de Cantoblanco, 28049 Madrid, Spain
2
Centro Universitario de la Defensa, Universidad Politécnica de Cartagena, C/Coronel López Peña s/n, Base Aérea de San Javier, Santiago de la Ribera, 30720 Murcia, Spain
*
Authors to whom correspondence should be addressed.
Viruses 2023, 15(4), 948; https://doi.org/10.3390/v15040948
Submission received: 2 March 2023 / Revised: 5 April 2023 / Accepted: 6 April 2023 / Published: 11 April 2023
(This article belongs to the Special Issue Bunyavirus, Volume II)
Browse Figures
Review Reports Versions Notes

Abstract

:
The Bunyavirales order is a large group of RNA viruses that includes important pathogens for humans, animals and plants. With high-throughput screening of clinically tested compounds we have looked for potential inhibitors of the endonuclease domain of a bunyavirus RNA polymerase. From a list of fifteen top candidates, five compounds were selected and their antiviral properties studied with Bunyamwera virus (BUNV), a prototypic bunyavirus widely used for studies about the biology of this group of viruses and to test antivirals. Four compounds (silibinin A, myricetin, L-phenylalanine and p-aminohippuric acid) showed no antiviral activity in BUNV-infected Vero cells. On the contrary, acetylsalicylic acid (ASA) efficiently inhibited BUNV infection with a half maximal inhibitory concentration (IC50) of 2.02 mM. In cell culture supernatants, ASA reduced viral titer up to three logarithmic units. A significant dose-dependent reduction of the expression levels of Gc and N viral proteins was also measured. Immunofluorescence and confocal microscopy showed that ASA protects the Golgi complex from the characteristic BUNV-induced fragmentation in Vero cells. Electron microscopy showed that ASA inhibits the assembly of Golgi-associated BUNV spherules that are the replication organelles of bunyaviruses. As a consequence, the assembly of new viral particles is also significantly reduced. Considering its availability and low cost, the potential usability of ASA to treat bunyavirus infections deserves further investigation.

Graphical Abstract

1. Introduction

Emerging and re-emerging viruses are a major health concern, and the validation of broad-spectrum antivirals that can be quickly applied to the clinic is, therefore, a priority. There are many viruses that can cause disease in humans. In particular, the World Health Organization’s (WHO) current list of prioritized diseases [1] includes three diseases caused by coronaviruses (SARS-CoV, MERS-CoV and SARS-CoV-2) and three more caused by members of the order Bunyavirales: Rift Valley Fever virus (RVFV), Crimean–Congo Hemorrhagic Fever virus (CCHFV) and Lassa Fever Virus (LFV). As a consequence of the intense research developed after the emergence of SARS-CoV-2, three antivirals have been approved to treat the associated pathology known as COVID-19: remdesivir, molnupiravir and paxlovid, the latter being a mixture of ritonavir and nirmatrelvir [2,3,4,5], and several more are currently being evaluated in clinical trials (https://clinicaltrials.gov/, accessed on 5 April 2023). However, there are no vaccines used worldwide, and no specific antivirals have been approved so far to treat bunyavirus infections. Ribavirin and Favipiravir, two inhibitors of viral RNA polymerases, are currently being studied in ongoing clinical trials for Lassa fever and Crimean–Congo hemorrhagic fever (https://clinicaltrials.gov/, accessed on 5 April 2023).
The Bunyavirales order comprises more than 500 viruses that are generally defined as bunyaviruses. With the exception of Hantaviruses and Arenaviruses, all viruses in the Bunyavirales order are transmitted by arthropods [6]. Like many other RNA viruses, bunyaviruses replicate their genome in association with intracellular membranes [7,8]. BUNV viral replication complexes (VRC) assemble in Golgi-derived single-membrane vesicles or spherules, which constitute the viral replication organelle (RO) and which are often connected with tubular structures [7,9,10]. Immature viral particles also assemble in Golgi membranes and undergo two maturation steps, the first one in the trans-Golgi compartment and the second one during virus egress that renders fully infectious virions [11]. All these steps are susceptible to being blocked by antiviral compounds. Antivirals are molecules that interfere with a step of the virus life cycle. There are basically two types: (i) direct acting antivirals (DAAs) that target viral components, and (ii) drugs against cell factors and pathways. DAAs are very specific and efficient for a particular virus or group of related viruses. Their main disadvantage is that virus variants might escape. On the other hand, inhibitors of cell components used by viruses should be efficient against emerging variants. However, a blockade of cell pathways can produce more severe secondary effects [12].
The incidence of bunyaviruses is growing due to several factors, such as climate change and globalization which favor spreading of arthropod vectors [13,14,15]. Thus, reports of novel bunyaviruses that cause severe disease are frequent [16,17,18,19,20]. For these reasons, we are looking for affordable compounds to treat bunyavirus infections. Viral RNA polymerases (vRNApol) are adequate targets for DAAs, because these proteins are essential for viral genome replication and share sequence as well as structural similarities across different families of RNA viruses [20,21,22]. Structural studies of a known influenza endonuclease inhibitor (2,4-dioxo-4-phenylbutanoic acid or DPBA) bound to the endonuclease domain of La Crosse bunyavirus (LACV) vRNApol, suggested that structure-based studies can identify compounds with antiviral activity against different RNA viruses [23]. DPBA is a known inhibitor of the influenza virus vRNPpol endonuclease domain that performs the cap-snatching process [24]. Bunyaviruses also perform cap-snatching to prime viral transcription. Basically, the vRNApol synthesizes viral messenger RNAs using short capped primers derived from cellular transcripts after endonuclease-mediated cleavage [25]. Structural and functional similarities between the endonucleases of influenza virus and bunyaviruses suggest that endonuclease inhibitors designed against the influenza virus endonuclease could be repurposed to target bunyaviruses [23,25,26].
To identify potential inhibitors of BUNV vRNApol, we looked for molecules already used to treat other diseases and performed an in silico study to identify compounds similar to DPBA that could bind to the vRNApol endonuclease domain of Hantaan virus, a bunyavirus belonging to hantavirus genus [27]. From a list of fifteen top candidates, we selected five compounds for functional studies in cell culture. One of them, acetylsalicylic acid, efficiently inhibited BUNV infection by blocking the biogenesis of the viral replication organelle.

2. Materials and Methods

2.1. Computational Models

The discovery of novel antivirals has initially been guided by two computational strategies. A ligand-based virtual screening was first performed by using BRUSELAS code [28]. In that approach, the structure of DPBA, a well-known inhibitor of the virus endonuclease reaction, is used as a chemical template for the search of similar compounds in the DrugBank (DB) database [29]. For the records, that similarity scheme is based on a fingerprint derived from a wide panel of numeric descriptors including molecular weight, number of rotatable bonds, number of hydrogen acceptor/donor atoms, topological polar surface area and octanol-water partition coefficient [28]. In a second stage, docking simulations with LeadFinder [30] as implemented in Flare [31] are conducted to screen DB against the crystal structure of the cap-snatching endonuclease of Hantaan virus, which is deposited in the Protein Data Bank with code 5IZE [27]. Docking allows assessment of the protein-ligand interaction by generating a series of ‘poses’ of the drug in the target binding site [32,33,34]. DB was selected as the working library because it compiles both Food and Drug Administration (FDA) approved drugs and experimental molecules and has been successfully used in drug repurposing [35,36].

2.2. Cells, Virus and Drugs

Vero cells (CCL-81) were obtained from the American Type Culture Collection and grown in DMEM (Sigma-Aldrich, St. Louis, MO, USA; D6429) supplemented with 10% fetal bovine serum (FBS, Reactiva S.A., Barcelona, Spain), 1% MEM nonessential amino acids (Sigma-Aldrich; M7145), 1% L-glutamine (Merck, Rahway, NJ, USA; 59202C) and 1% penicillin/streptomycin (Sigma-Aldrich; P4333).
Bunyamwera virus (BUNV, VR-87) was obtained from the American Type Culture Collection (Virginia, WV, USA). Propagation and production of BUNV stocks was performed in Vero cells as previously described [9,37]. Virus titers measured as plaque forming units, or PFU per milliliter, were determined by standard plaque assay in Vero cells using an agar overlay [37].
Acetylsalicylic acid (A5376), myricetin (M6760), L-phenylalanine (5482), p-aminohippuric acid (A1422) and silibinin A (S0417) were purchased from Sigma-Aldrich, dissolved in ethanol (acetylsalicylic acid and myricetin), water (L-phenylalanine and p-aminohippuric acid) or dimethyl sulfoxide (DMSO) (silibinin A), as recommended by the manufacturer, to prepare stock solutions and stored at −20 °C.

2.3. Cytotoxicity Assays

Cytotoxicity of the drugs was analyzed with the MTT (3-[4,5-dimethylthiazol-2-yl]-2,5 diphenyl tetrazolium bromide) assay, which measures the mitochondrial dehydrogenase activity of living cells [38]. Vero cells were seeded in a 96-well plate at a density of 10,000 cells per well and cultured until an 80% confluency was reached. Cells were incubated with serial dilutions of the drugs in DMEM for 24 h and then treated with MTT reagent (Sigma-Aldrich; M5655) at 0.5 mg/mL for 4 h. The formazan crystals were dissolved by incubation with the MTT solvent (10% SDS, 0.01 M HCl in 85% isopropanol) for 30 min in an orbital shaker protected from light. The amount of MTT formazan formed was determined spectroscopically at 570 nm in an ELISA plate reader. The background was also measured and subtracted at 690 nm. The cytotoxic effect of the drugs was measured with the CC50 calculated from dose–response curves. Three independent replicas were analyzed.

2.4. Viral Infection and Antiviral Activity

Vero cells were absorbed with BUNV at a multiplicity of infection (MOI) of 1 PFU per cell (viral stock diluted in DMEM culture medium supplemented to contain 2% FBS) at 37 °C for 1 h, gently moving the plate side to side every 15 min. After the inoculum was removed, fresh DMEM supplemented with 2% FBS was added, and cells were incubated for 8 h at 37 °C.
The effect of the compounds on the viral infection was determined by immunofluorescence and fluorescence microscopy. Vero cells were seeded in 96-well flat-bottom plates and incubated at 37 °C to reach 80–90% confluency. Cell monolayers were infected with BUNV at an MOI of 1 PFU/cell for 1 h at 37 °C. The viral inoculum was then removed and a serial dilution of acetylsalicylic acid (ASA), myricetin, L-phenylalanine, p-aminohippuric acid and silibinin A diluted in DMEM supplemented with 2% FBS was added. Infected cells not treated with drugs were also analyzed as controls. After 8 h, cells were fixed with 4% paraformaldehyde (PFA, TAAB Laboratories, Berks, England) for 20 min and washed three times with PBS. The percentage of infected cells was determined by immunofluorescence with an antibody specific for the viral N protein and fluorescence microscopy (see below). To analyze the recovery of BUNV after removing ASA from cell cultures, Vero cells were absorbed with BUNV at an MOI of 1 PFU/cell for 1 h and treated with ASA 10 mM for 8 h. The medium was removed and replaced by new DMEM with 2% FBS. After incubation at different time points (2, 4, 6 and 8 h), cells were fixed for 20 min with 4% PFA, washed three times with PBS and processed for immunofluorescence with an antibody specific for the viral N protein and studied by fluorescence microscopy. Images were acquired using a Leica DMi8 S widefield epifluorescence microscope and analyzed with the Image J software. The IC50 was calculated from the percentage of infected cells. Inhibition data were plotted as dose–effect curves fitted to a nonlinear regression model using the GraphPad Prism 9.5.0 software. All experiments were replicated three times.

2.5. Western Blot

Vero cells were infected with BUNV at 1 PFU/cell and treated with ASA at different doses (5, 10 and 15 mM) in 6-well plates. Non-infected cells were also processed as controls. Eight hours after infection, cells were harvested with CellStripper™ reagent (Corning®, Somerville, MA, USA; 25-056-CI) and incubated with 50 μL of RIPA lysis buffer (Thermo Scientific™, Waltham, MA, USA; 89900) with protease inhibitors for 30 min on ice. The amount of protein in cell lysates was quantified to ensure equal loading (10 μg) in Western blot gels using Pierce™ BCA Protein Assay Kit (Thermo Scientific™; 23227) according to the manufacturer’s specifications. The samples were analyzed by SDS-PAGE and transferred to polyvinylidene fluoride (PVDF) membranes by standard blotting procedures. Membranes were blocked overnight at 4 °C with PBS containing non-fat dry milk and 0.05% Tween 20, incubated for 1 h at room temperature with a rabbit anti-BUNV polyclonal antiserum, kindly provided by Prof. Richard Elliott [37] diluted 1:2000 in blocking buffer and a mouse anti-α-tubulin monoclonal antibody (Sigma-Aldrich; T5168) diluted 1:1000 in blocking buffer. After washing with this buffer, membranes were incubated with horseradish peroxidase (HRP)-conjugated secondary antibodies for 1 h at room temperature (RT). The target proteins were detected by ECL solutions (SuperSignal™ West Dura Extended Duration Substrate, Thermo Scientific™) with a ChemiDoc Imager (BioRad). The band intensity of Gc and N viral proteins was calculated using Image Lab software (BioRad) and normalized to α-tubulin. Three independent replicas were analyzed.

2.6. Immunofluorescence and Confocal Microscopy

Cells were grown on sterile glass coverslips, fixed with 4% PFA, and washed three times with PBS. Fixed cells were permeabilized and blocked by incubation for 40 min with PBS containing 0.25% saponin and 2% FBS. Cells were incubated for 1 h with primary antibodies diluted 1/200 in PBS with saponin and FBS. The following antibodies were used: a rabbit polyclonal antibody specific for the BUNV nucleocapsid (N) protein [7], a mouse monoclonal antibody specific for BUNV Gc glycoprotein, kindly provided by Prof. Richard M. Elliott [39] and a rabbit anti-Giantin polyclonal antibody that is a marker of the Golgi complex (BioLegend, San Diego, CA, United States; 924302). Alexa Fluor–conjugated antibodies (Invitrogen) were used as secondary antibodies diluted 1/500 in PBS with 0.25% saponin and 2% FBS. Cell nuclei were stained with 4′,6-diamidino-2-phenylindole (DAPI) (Sigma-Aldrich) diluted 1/200 in PBS with 0.25% saponin and 2% FBS for 20 min. All incubations were conducted at RT. Coverslips were mounted using Prolong-Gold (Life Technologies, Carlsbad, CA, USA). Confocal microscopy images were acquired using a Leica STELLARIS 5 confocal multispectral microscope equipped with an HCX PL APO 63.0 X/1.4 NA oil objective and LAS X Life Science software (Leica Microsystems).

2.7. Transmission Electron Microscopy (TEM)

Cells were grown on sterile glass coverslips in 6-well plates at 37 °C to reach a 90% confluency, before infection and drug treatment. Vero cells were mock-infected or infected with BUNV at an MOI of 1 PFU/cell for 1 h at 37 °C, and incubated in the absence or presence of three different concentrations (5, 10 and 15 mM) of ASA and fixed at 8 hpi with 1% glutaraldehyde in 0.2 M Hepes buffer, pH 7.4, at RT for 2 h. Cells were washed twice with Hepes 0.4 M, pH 7.4, postfixed by incubation at 4 °C for 1 h with a mixture of 1% osmium tetroxide and 0.8% potassium ferricyanide in water and washed three times with Hepes 0.4 M. Cells were then dehydrated in 5-min steps with increasing concentrations of ethanol (25, 50, 70 and 95%) in 10-min steps, and twice in 100% ethanol at 4 °C. Samples were processed for flat embedding in the epoxy resin EML-812 (TAAB Laboratories) by incubating for 2 h with a 1:1 mixture of ethanol–resin at RT. Cells were infiltrated for 16 h in pure resin and polymerized at 60 °C for 48 h. Ultrathin (∼50–70 nm) oriented serial sections were obtained using an UC6 ultramicrotome (Leica Microsystems), collected on uncoated 300-mesh copper grids (TAAB Laboratories), stained with saturated uranyl acetate and lead citrate, and imaged by TEM. Images were acquired using a Jeol JEM 1400 Flash electron microscope operating at 120 kV and equipped with a Gatan CMOS Oneview camera. At least 100 cells per condition were studied by TEM.

2.8. Quantification and Statistical Analysis

Experiments were conducted with a minimum of three independent replicates. Data presented in the text are expressed as the mean ± SD. Statistical significance was determined using two-sample unequal variance student’s t-test with two-tailed distribution (α = 0.05). Graphs were prepared and statistical analyses were conducted using GraphPad Prism 9.5.0 software.

3. Results

3.1. Antiviral Selection by Molecular Modeling

As described above, our computational models combine both ligand-based and structure-based molecular models. The former is based on the DPBA inhibitor, while the latter uses the crystal structure of the endonuclease of Hantaan virus resolved by Reguera, Cusack and co-workers [27]. These authors demonstrated the essential role played not only by the catalytic His36, but also by other critical residues, e.g., Glu54, Asp97 and Lys124, as they are also involved in Mn2+ ion binding [27]. Virtual Screening with docking was consequently focused on that specific region. The performed simulations lead to a series of five hits: ASA, L-phenylalanine and p-aminohippuric acid were selected by similarity; myricetin and silibinin A are associated with the more negative (stable) binding poses in the catalytic pocket. Figure 1 summarizes predicted interactions energies and illustrates the binding mode at the endonuclease catalytic site.

3.2. Antiviral Activity in Cells Infected with BUNV

To evaluate the antiviral activity of the drugs, Vero cells were infected with BUNV for 1 h and treated in parallel with increasing concentrations of the selected compounds. Drug cytotoxicity was first measured using the MTT test, and safe drug doses were used to determine the percentage of infected cells by immunofluorescence, using an antibody specific for the BUNV N protein. Myricetin, L-phenylalanine, p-aminohippuric acid and silibinin A showed no antiviral activity. However, ASA presented antiviral activity against BUNV (Figure 2). The percentage of cells with an N viral protein fluorescence signal decreased in infected Vero cells treated with increasing concentrations of ASA (Figure 2A–E). By regression analysis, we calculated the concentration of compound required to inhibit 50% of virus infection (IC50), and the concentration for 50% cytotoxic effect (CC50). The graph in Figure 2F shows that ASA efficiently inhibited BUNV infection at non-toxic concentrations, with an IC50 of 2.02 mM. To study ASA antiviral capacity in detail, further assays were performed with this compound at 5, 10 and 15 mM.

3.3. Effect of ASA on Viral Titer and Viral Protein Expression

To study the action of ASA in the production of infectious viral particles, Vero cells were infected at an MOI of 1 PFU/cell for 1 h at 37 °C, and incubated in the absence or presence of 5, 10 or 15 mM of ASA. At 8 hpi, supernatants were collected and titrated by standard plaque assay [37]. ASA treatment dramatically reduced the number of lysis plaques (Figure 3A). The amount of infectious viral particles in cell culture supernatants was reduced by 10 to 1000-fold (Figure 3B).
The effect of ASA in the expression of BUNV proteins was measured by Western blot (WB) using a rabbit anti-BUNV polyclonal antiserum that mainly recognizes the nucleocapsid N protein and the large Gc glycoprotein of BUNV (Figure 4). The WB image on the left shows that expression of N and Gc viral proteins in Vero cells was significantly reduced with the ASA treatment. The table of relative frequencies of Figure 4 shows a dose-dependent reduction of N and Gc expression levels in cells treated with ASA compared to untreated cells. These results suggest that the decrease in the amount of infectious viral particles in cell supernatants from ASA-treated cells shown in Figure 3, is due to the blockade of an early step of BUNV infection.

3.4. Distribution of Viral Gc and N Proteins in Cells Treated with ASA

BUNV Gc glycoprotein is a marker of viral factories that assemble in the Golgi complex [9]. BUNV N is a multifunctional protein that is synthesized in large amounts early in infection and accumulates in the cytosol and in association with cell membranes [7,40]. To study the effects of ASA in Gc and N localization in cells, we processed infected cells, treated or untreated with ASA by immunofluorescence using antibodies specific for Gc, N and the Golgi marker Giantin. As mentioned before, the number of infected cells decreased with ASA treatment. The percentage of cells with N positive signal was 48%, 22% and 8% of cells treated with ASA 5, 10 and 15 mM, respectively. A Gc positive signal was presented in 22%, 14% and 1% of cells treated with ASA 5, 10 and 15 mM, respectively. Confocal microscopy showed that the fluorescence signal of both viral proteins decreased in a dose-dependent manner (Figure 2A–D and Figure 5A–D). Moreover, ASA inhibited the characteristic BUNV-induced Golgi fragmentation and dispersion. For comparison, the normal morphology and intracellular localization of the Golgi complex in uninfected Vero cells non-treated with ASA are shown in Figure S1. With ASA, Gc was retained in perinuclear areas where it colocalized with Golgi (Figure 5A–D). The number of cells with N fluorescent signal as well as the intensity of the signal decreased with increasing concentrations of ASA, but the cytosolic distribution of N protein did not change in the presence of ASA (Figure 2A–D). Three-dimensional reconstructions (Figure 5E,F) showed that in non-treated cells, Gc is dispersed in the cytosol where it associates to fragments of Golgi (Figure 5E) and to virions inside secretory vesicles [9]. With ASA, both Gc and Giantin signals grouped in the perinuclear region (arrows in Figure 5F). With 5 mM ASA, no significant differences were observed in the number of cells with N and Gc protein signals; however, there was a significant reduction with the 10 and 15 mM treatments (Figure 5G). Quantification of Gc distribution in infected cells not treated and treated with 5, 10 or 15 mM of ASA, is shown in Figure 5H that shows two phenotypes corresponding to scattered or perinuclear distribution of Gc in the cytosol (Figure 5H). This analysis revealed significant differences in the localization of Gc signal in infected cells without ASA or in cells treated with 5 mM ASA, compared to cells treated with higher ASA doses.

3.5. Effect of ASA on the Morphology of Viral Structures and Cellular Organelles

To further investigate the mechanism of action of ASA, we studied the ultrastructure of infected Vero cells not treated or treated with ASA (5, 10 and 15 mM) by transmission electron microscopy (Figure 6). Infected cells not treated with the drug showed a recruitment of cellular organelles to particular areas of the cytoplasm forming the viral factory (VF), structure that contains elements of the Golgi complex, endoplasmic reticulum (ER) and mitochondria (Figure 6A,B). Extracellular virions were also seen close to the plasma membrane of infected cells (inset in Figure 6A). High magnification images of VFs show characteristic spherules, tubes and viral particles (Figure 6C).
In infected cells treated with 5 mM ASA, cellular organelles were not apparently recruited to a specific area into the cytoplasm forming a factory (Figure 6D). However, some infected and ASA-treated cells showed the BUNV factory (Figure 6E) with spherules in Golgi membranes (Figure 6F). In infected cells treated with 10 mM ASA, the cellular organelles were randomly distributed through the cytoplasm and no characteristic VFs were seen (Figure 6G). Higher magnification images showed altered Golgi complexes with swollen sacculi, and neither spherules nor viral particles were seen in the majority of cells (Figure 6H,I). Similar results were obtained with infected cells treated with 15 mM ASA (Figure 6J–L). Quantification of viral spherules, intracellular and extracellular viral particles is shown in Figure 7. Tables of frequencies and graphs show that ASA treatment remarkably reduced the number of spherules and viral particles, being practically undetectable with the treatment of 15 mM ASA. We also analyzed the effect of the ASA treatment in the ultrastructure of uninfected Vero cells (Figure S2). The distribution of cellular organelles was not modified in cells treated with lower or higher concentrations of ASA in comparison to untreated cells (Figure S2, images on the left). Closer views of the Golgi complex, ER and mitochondria did not show visible alterations of those cellular organelles with any of the ASA treatments (Figure S2, images in the middle). High magnification views of the Golgi complex showed that Golgi sacculi were not swollen (Figure S2, images on the right), a feature detected only in BUNV-infected and ASA-treated cells. These results suggest that ASA treatment interferes with BUNV infection by blocking the formation of viral spherules.

3.6. Recovery of Cell Morphology and BUNV Infection after Elimination of ASA

To evaluate the capacity of BUNV to recover infection after removing ASA from cell cultures, Vero cells were infected and treated with 10 mM ASA, and 8 h later the medium was replaced by fresh DMEM with 2% FBS. Cells were then incubated for 0, 2, 4, 6 or 8 h without ASA and processed for immunofluorescence and fluorescence microscopy, titration by plaque assay and electron microscopy. The percentage of infected cells increased after removing the drug and reached values around 90–100% of infected cells after 8 h of incubation without ASA (Figure 8A). In addition, infectious viral particles were detected by plaque assay after incubation without drug (Figure 8B). Electron microscopy of infected cells treated with ASA showed Golgi complexes with swollen sacculi and no viral structures (Figure 8C). After incubation without drug for 2, 4, 6 or 8 h, the Golgi complex had already recovered its normal morphology and showed viral spherules and viral particles (Figure 8D–H). Therefore, after the elimination of ASA, the Golgi complex recovered and BUNV was able to induce a productive infection.

4. Discussion

In our search for antivirals to combat bunyavirus infections, we have discovered that ASA inhibits the infection of a bunyavirus by interfering with the biogenesis of its replication organelle. ASA is on the WHO model list of essential medicines [41] that includes the most efficacious, safe and cost–effective medicines for priority conditions. In addition, ASA is one of the most widely used medications in the world to treat pain, rheumatic fever, arthritis and inflammation [42]. A potential use of ASA to reduce the long-term risk of several cancers has been proposed [42,43,44,45] and ASA’s antiviral properties have also been reported. ASA inhibits the replication of flaviviruses, such as Japanese encephalitis virus (JEV), dengue virus (DENV) and hepatitis C virus (HCV) [46,47]. ASA also blocks the propagation of respiratory viruses such as influenza A virus and the human rhinovirus (HRV) [48].
Structural and functional similarities between the endonucleases of the influenza virus and bunyaviruses suggested that endonuclease inhibitors designed against influenza virus endonuclease could be repurposed to target bunyaviruses [23,25,26]. That earlier evidence motivated the present screening of the DrugBank database, a library that contains ca. 10,000 compounds previously used to treat humans, including both approved and experimental drugs [29]. Our hypothesis was that a similar drug to DPBA, a known inhibitor of influenza virus endonuclease [27], might also exhort a biological activity against bunyaviruses. That search was not manually conducted but guided by computational methods that used fingerprints to systematically assess similarities. Virtual screening also confirms the compatibility with the catalytic site, a prerequisite to halt the endonuclease machinery.
The antiviral properties of diketoacid inhibitors such as DPBA have been tested in cells infected with different viruses such as the influenza virus [49], the human immunodeficiency virus 1 (HIV-1) [50,51] and bunyaviruses [52,53]. To our knowledge, no clinical trials for diketoacid inhibitors to treat virus disease have been initiated.
Here we show that ASA induced a decrease of up to three logarithmic units in viral titers of culture supernatants from BUNV-infected Vero cells. It also produced a significant dose-dependent reduction of the expression levels of Gc and N viral proteins and inhibited the assembly of Golgi-associated BUNV spherules and new viral particles. Light and electron microscopy showed that ASA protected the Golgi complex from the characteristic BUNV-induced fragmentation and caused swelling of Golgi sacculi in BUNV-infected cells only, while in mock-infected cells treated with ASA the Golgi complex looked normal. These results suggest that the way vRNApol molecules interact with Golgi membranes changes when bound to ASA, something that probably impairs the structural role of vRNApol in the biogenesis of viral replication organelles and VRCs. Similar results were obtained in Vero cells infected with BUNV and treated with the typical vRNApol inhibitor Ribavirin where viral spherules were absent or disrupted [54].
The elimination of ASA from the cell culture medium led to the recovery of the normal morphology of Golgi and of all viral structures, such as spherules (ROs), intracellular and extracellular viral particles. A positive aspect of these results is that the secondary effects in cell morphology (swollen Golgi) did not affect cell viability (Figure 2) and were reversible (Figure 8). However, the reversibility of the antiviral effect has to be also taken into consideration for stablishing the most adequate treatment conditions in experiments with animals, such as number of doses over time. When using therapeutic doses, plasma levels of ASA can range from 0.2 to 0.5 mM and the toxic level is stablished above 8 mM [55], which is clearly higher than the concentration of 2 mM (IC50) that efficiently inhibited BUNV replication in cell culture. However, due to its anti-coagulant properties, ASA could not be used to treat patients with hemorrhagic fever. Examples of bunyavirus infections that could be treated with ASA and others that should not be treated with ASA are shown in Table 1. ASA is a cheap and safe drug with the potential to become a broad-spectrum antiviral.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/v15040948/s1, Figure S1: Morphology and localization of the Golgi complex in mock-infected Vero cells in the absence of ASA; Figure S2: Transmission electron microscopy of uninfected cells untreated or treated with ASA.

Author Contributions

Conceptualization, S.Y.F.-S., J.P.C.-C., I.F.d.C. and C.R.; methodology, J.P.C.-C.; software, J.P.C.-C.; validation, S.Y.F.-S., I.F.d.C. and J.P.C.-C.; formal analysis, S.Y.F.-S., J.P.C.-C. and I.F.d.C.; investigation, S.Y.F.-S., J.P.C.-C. and I.F.d.C.; resources, J.P.C.-C. and C.R.; writing—original draft preparation, I.F.d.C., C.R.; writing—review and editing, S.Y.F.-S., J.P.C.-C., I.F.d.C. and C.R.; visualization, S.Y.F.-S., J.P.C.-C., I.F.d.C. and C.R.; supervision, J.P.C.-C. and C.R.; project administration, J.P.C.-C. and C.R.; funding acquisition, C.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Spanish Ministry of Science and Innovation [RTI2018-094445-B-I00. MCI/AEI/FEDER, UE].

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data will be available from the corresponding author at reasonable request.

Acknowledgments

We acknowledge that the results of this research have been achieved using the DECI-17 resource NIIF High Performance Computing supercomputer at University of Debrecen, Hungary, with support from the PRACE.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. WHO. Prioritizing Diseases for Research and Development in Emergency Contexts. Available online: https://www.who.int/activities/prioritizing-diseases-for-research-and-development-in-emergency-contexts (accessed on 5 April 2023).
  2. Buxeraud, J.; Faure, S.; Fougere, E. Nirmatrelvir/ritonavir (Paxlovid(R)), a treatment for COVID-19. Actual. Pharm. 2022, 61, 10–12. [Google Scholar] [PubMed]
  3. Jayk Bernal, A.; Gomes da Silva, M.M.; Musungaie, D.B.; Kovalchuk, E.; Gonzalez, A.; Delos Reyes, V.; Martin-Quiros, A.; Caraco, Y.; Williams-Diaz, A.; Brown, M.L.; et al. Molnupiravir for Oral Treatment of COVID-19 in Nonhospitalized Patients. N. Engl. J. Med. 2022, 386, 509–520. [Google Scholar] [CrossRef] [PubMed]
  4. Lamontagne, F.; Agarwal, A.; Rochwerg, B.; Siemieniuk, R.A.; Agoritsas, T.; Askie, L.; Lytvyn, L.; Leo, Y.S.; Macdonald, H.; Zeng, L.; et al. A living WHO guideline on drugs for COVID-19. BMJ 2020, 370, m3379. [Google Scholar] [CrossRef] [PubMed]
  5. Schafer, A.; Martinez, D.R.; Won, J.J.; Meganck, R.M.; Moreira, F.R.; Brown, A.J.; Gully, K.L.; Zweigart, M.R.; Conrad, W.S.; May, S.R.; et al. Therapeutic treatment with an oral prodrug of the remdesivir parental nucleoside is protective against SARS-CoV-2 pathogenesis in mice. Sci. Transl. Med. 2022, 14, eabm3410. [Google Scholar] [CrossRef]
  6. Horne, K.M.; Vanlandingham, D.L. Bunyavirus-vector interactions. Viruses 2014, 6, 4373–4397. [Google Scholar] [CrossRef] [Green Version]
  7. Fontana, J.; Lopez-Montero, N.; Elliott, R.M.; Fernandez, J.J.; Risco, C. The unique architecture of Bunyamwera virus factories around the Golgi complex. Cell Microbiol. 2008, 10, 2012–2028. [Google Scholar] [CrossRef]
  8. Mackenzie, J. Wrapping things up about virus RNA replication. Traffic 2005, 6, 967–977. [Google Scholar] [CrossRef]
  9. Salanueva, I.J.; Novoa, R.R.; Cabezas, P.; Lopez-Iglesias, C.; Carrascosa, J.L.; Elliott, R.M.; Risco, C. Polymorphism and structural maturation of bunyamwera virus in Golgi and post-Golgi compartments. J. Virol. 2003, 77, 1368–1381. [Google Scholar] [CrossRef] [Green Version]
  10. Sun, Y.; Li, J.; Gao, G.F.; Tien, P.; Liu, W. Bunyavirales ribonucleoproteins: The viral replication and transcription machinery. Crit. Rev. Microbiol. 2018, 44, 522–540. [Google Scholar] [CrossRef]
  11. Novoa, R.R.; Calderita, G.; Arranz, R.; Fontana, J.; Granzow, H.; Risco, C. Virus factories: Associations of cell organelles for viral replication and morphogenesis. Biol. Cell 2005, 97, 147–172. [Google Scholar] [CrossRef]
  12. Garcia-Serradilla, M.; Risco, C.; Pacheco, B. Drug repurposing for new, efficient, broad spectrum antivirals. Virus Res. 2019, 264, 22–31. [Google Scholar] [CrossRef]
  13. Buhler, K.J.; Dibernardo, A.; Pilfold, N.W.; Harms, N.J.; Fenton, H.; Carriere, S.; Kelly, A.; Schwantje, H.; Aguilar, X.F.; Leclerc, L.M.; et al. Widespread Exposure to Mosquitoborne California Serogroup Viruses in Caribou, Arctic Fox, Red Fox, and Polar Bears, Canada. Emerg. Infect. Dis. 2023, 29, 54–63. [Google Scholar] [CrossRef]
  14. Ciota, A.T.; Keyel, A.C. The Role of Temperature in Transmission of Zoonotic Arboviruses. Viruses 2019, 11, 1013. [Google Scholar] [CrossRef] [Green Version]
  15. Dai, S.; Deng, F.; Wang, H.; Ning, Y. Crimean-Congo Hemorrhagic Fever Virus: Current Advances and Future Prospects of Antiviral Strategies. Viruses 2021, 13, 1195. [Google Scholar] [CrossRef]
  16. Dutuze, M.F.; Nzayirambaho, M.; Mores, C.N.; Christofferson, R.C. A Review of Bunyamwera, Batai, and Ngari Viruses: Understudied Orthobunyaviruses With Potential One Health Implications. Front. Vet. Sci. 2018, 5, 69. [Google Scholar] [CrossRef] [Green Version]
  17. Edridge, A.W.D.; van der Hoek, L. Emerging orthobunyaviruses associated with CNS disease. PLoS Negl. Trop. Dis. 2020, 14, e0008856. [Google Scholar] [CrossRef] [PubMed]
  18. Evans, A.B.; Peterson, K.E. Throw out the Map: Neuropathogenesis of the Globally Expanding California Serogroup of Orthobunyaviruses. Viruses 2019, 11, 794. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Li, J.; Li, S.; Yang, L.; Cao, P.; Lu, J. Severe fever with thrombocytopenia syndrome virus: A highly lethal bunyavirus. Crit. Rev. Microbiol. 2021, 47, 112–125. [Google Scholar] [CrossRef]
  20. Ter Horst, S.; Conceicao-Neto, N.; Neyts, J.; Rocha-Pereira, J. Structural and functional similarities in bunyaviruses: Perspectives for pan-bunya antivirals. Rev. Med. Virol. 2019, 29, e2039. [Google Scholar] [CrossRef]
  21. Bivacqua, R.; Barreca, M.; Spano, V.; Raimondi, M.V.; Romeo, I.; Alcaro, S.; Andrei, G.; Barraja, P.; Montalbano, A. Insight into non-nucleoside triazole-based systems as viral polymerases inhibitors. Eur. J. Med. Chem. 2023, 249, 115136. [Google Scholar] [CrossRef] [PubMed]
  22. Ferrer-Orta, C.; Arias, A.; Escarmis, C.; Verdaguer, N. A comparison of viral RNA-dependent RNA polymerases. Curr. Opin. Struct. Biol. 2006, 16, 27–34. [Google Scholar] [CrossRef] [PubMed]
  23. Reguera, J.; Weber, F.; Cusack, S. Bunyaviridae RNA polymerases (L-protein) have an N-terminal, influenza-like endonuclease domain, essential for viral cap-dependent transcription. PLoS Pathog. 2010, 6, e1001101. [Google Scholar] [CrossRef]
  24. Dias, A.; Bouvier, D.; Crepin, T.; McCarthy, A.A.; Hart, D.J.; Baudin, F.; Cusack, S.; Ruigrok, R.W. The cap-snatching endonuclease of influenza virus polymerase resides in the PA subunit. Nature 2009, 458, 914–918. [Google Scholar] [CrossRef]
  25. Olschewski, S.; Cusack, S.; Rosenthal, M. The Cap-Snatching Mechanism of Bunyaviruses. Trends Microbiol. 2020, 28, 293–303. [Google Scholar] [CrossRef] [Green Version]
  26. Jones, R.; Lessoued, S.; Meier, K.; Devignot, S.; Barata-Garcia, S.; Mate, M.; Bragagnolo, G.; Weber, F.; Rosenthal, M.; Reguera, J. Structure and function of the Toscana virus cap-snatching endonuclease. Nucleic Acids Res. 2019, 47, 10914–10930. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Reguera, J.; Gerlach, P.; Rosenthal, M.; Gaudon, S.; Coscia, F.; Gunther, S.; Cusack, S. Comparative Structural and Functional Analysis of Bunyavirus and Arenavirus Cap-Snatching Endonucleases. PLoS Pathog. 2016, 12, e1005636. [Google Scholar] [CrossRef] [Green Version]
  28. Banegas-Luna, A.J.; Ceron-Carrasco, J.P.; Puertas-Martin, S.; Perez-Sanchez, H. BRUSELAS: HPC Generic and Customizable Software Architecture for 3D Ligand-Based Virtual Screening of Large Molecular Databases. J. Chem. Inf. Model. 2019, 59, 2805–2817. [Google Scholar] [CrossRef]
  29. Wishart, D.S.; Knox, C.; Guo, A.C.; Cheng, D.; Shrivastava, S.; Tzur, D.; Gautam, B.; Hassanali, M. DrugBank: A knowledgebase for drugs, drug actions and drug targets. Nucleic Acids Res. 2008, 36, D901–D906. [Google Scholar] [CrossRef] [PubMed]
  30. BioMolTech®. Lead Finder™, 6.0.1.; BioMolTech: Toronto, ON, Canada. Available online: http://www.cresset-group.com/lead-finder/(accessed on 5 April 2023).
  31. Cresset®. Flare™, Litlington: Cambridgeshire, UK. Available online: http://www.cresset-group.com/flare/(accessed on 5 April 2023).
  32. Bauer, M.R.; Mackey, M.D. Electrostatic Complementarity as a Fast and Effective Tool to Optimize Binding and Selectivity of Protein-Ligand Complexes. J. Med. Chem. 2019, 62, 3036–3050. [Google Scholar] [CrossRef]
  33. Cheeseright, T.; Mackey, M.; Rose, S.; Vinter, A. Molecular field extrema as descriptors of biological activity: Definition and validation. J. Chem. Inf. Model. 2006, 46, 665–676. [Google Scholar] [CrossRef]
  34. Kuhn, M.; Firth-Clark, S.; Tosco, P.; Mey, A.; Mackey, M.; Michel, J. Assessment of Binding Affinity via Alchemical Free-Energy Calculations. J. Chem. Inf. Model. 2020, 60, 3120–3130. [Google Scholar] [CrossRef]
  35. Diallo, B.N.; Swart, T.; Hoppe, H.C.; Tastan Bishop, O.; Lobb, K. Potential repurposing of four FDA approved compounds with antiplasmodial activity identified through pr.roteome scale computational drug discovery and in vitro assay. Sci. Rep. 2021, 11, 1413. [Google Scholar] [CrossRef] [PubMed]
  36. Yuan, S.; Chan, J.F.; den-Haan, H.; Chik, K.K.; Zhang, A.J.; Chan, C.C.; Poon, V.K.; Yip, C.C.; Mak, W.W.; Zhu, Z.; et al. Structure-based discovery of clinically approved drugs as Zika virus NS2B-NS3 protease inhibitors that potently inhibit Zika virus infection in vitro and in vivo. Antiviral. Res. 2017, 145, 33–43. [Google Scholar] [CrossRef]
  37. Watret, G.E.; Pringle, C.R.; Elliott, R.M. Synthesis of bunyavirus-specific proteins in a continuous cell line (XTC-2) derived from Xenopus laevis. J. Gen. Virol. 1985, 66 Pt 3, 473–482. [Google Scholar] [CrossRef]
  38. Tolosa, L.; Donato, M.T.; Gomez-Lechon, M.J. General Cytotoxicity Assessment by Means of the MTT Assay. Methods Mol. Biol. 2015, 1250, 333–348. [Google Scholar]
  39. Lappin, D.F.; Nakitare, G.W.; Palfreyman, J.W.; Elliott, R.M. Localization of Bunyamwera bunyavirus G1 glycoprotein to the Golgi requires association with G2 but not with NSm. J. Gen. Virol. 1994, 75 Pt 12, 3441–3451. [Google Scholar] [CrossRef] [PubMed]
  40. Eifan, S.A.; Elliott, R.M. Mutational analysis of the Bunyamwera orthobunyavirus nucleocapsid protein gene. J. Virol. 2009, 83, 11307–11317. [Google Scholar] [CrossRef] [Green Version]
  41. WHO. WHO Model List of Essential Medicines-22nd List. 2021. Available online: https://www.who.int/publications/i/item/WHO-MHP-HPS-EML-2021.02 (accessed on 5 April 2023).
  42. Desborough, M.J.R.; Keeling, D.M. The aspirin story-from willow to wonder drug. Br. J. Haematol. 2017, 177, 674–683. [Google Scholar] [CrossRef] [Green Version]
  43. Algra, A.M.; Rothwell, P.M. Effects of regular aspirin on long-term cancer incidence and metastasis: A systematic comparison of evidence from observational studies versus randomised trials. Lancet Oncol. 2012, 13, 518–527. [Google Scholar] [CrossRef]
  44. Li, X.; Wu, S.; Yu, Y. Aspirin Use and the Incidence of Hepatocellular Carcinoma in Patients With Hepatitis B Virus or Hepatitis C Virus Infection: A Meta-Analysis of Cohort Studies. Front. Med. 2020, 7, 569759. [Google Scholar] [CrossRef]
  45. Rothwell, P.M.; Wilson, M.; Elwin, C.E.; Norrving, B.; Algra, A.; Warlow, C.P.; Meade, T.W. Long-term effect of aspirin on colorectal cancer incidence and mortality: 20-year follow-up of five randomised trials. Lancet 2010, 376, 1741–1750. [Google Scholar] [CrossRef]
  46. Liao, C.L.; Lin, Y.L.; Wu, B.C.; Tsao, C.H.; Wang, M.C.; Liu, C.I.; Huang, Y.L.; Chen, J.H.; Wang, J.P.; Chen, L.K. Salicylates inhibit flavivirus replication independently of blocking nuclear factor kappa B activation. J. Virol. 2001, 75, 7828–7839. [Google Scholar] [CrossRef] [Green Version]
  47. Trujillo-Murillo, K.; Rincon-Sanchez, A.R.; Martinez-Rodriguez, H.; Bosques-Padilla, F.; Ramos-Jimenez, J.; Barrera-Saldana, H.A.; Rojkind, M.; Rivas-Estilla, A.M. Acetylsalicylic acid inhibits hepatitis C virus RNA and protein expression through cyclooxygenase 2 signaling pathways. Hepatology 2008, 47, 1462–1472. [Google Scholar] [CrossRef]
  48. Glatthaar-Saalmuller, B.; Mair, K.H.; Saalmuller, A. Antiviral activity of aspirin against RNA viruses of the respiratory tract-an in vitro study. Influenza Other Respir. Viruses 2017, 11, 85–92. [Google Scholar] [CrossRef]
  49. Song, M.S.; Kumar, G.; Shadrick, W.R.; Zhou, W.; Jeevan, T.; Li, Z.; Slavish, P.J.; Fabrizio, T.P.; Yoon, S.W.; Webb, T.R.; et al. Identification and characte…erization of influenza variants resistant to a viral endonuclease inhibitor. Proc. Natl. Acad. Sci. USA 2016, 113, 3669–3674. [Google Scholar] [CrossRef] [Green Version]
  50. Alves, C.N.; Marti, S.; Castillo, R.; Andres, J.; Moliner, V.; Tunon, I.; Silla, E. A quantum mechanic/molecular mechanic study of the wild-type and N155S mutant HIV-1 integrase complexed with diketo acid. Biophys. J. 2008, 94, 2443–2451. [Google Scholar] [CrossRef] [Green Version]
  51. Corona, A.; di Leva, F.S.; Rigogliuso, G.; Pescatori, L.; Madia, V.N.; Subra, F.; Delelis, O.; Esposito, F.; Cadeddu, M.; Costi, R.; et al. New insights into the interaction between pyrrolyl diketoacids and HIV-1 integrase active site and comparison with RNase H. Antivir. Res. 2016, 134, 236–243. [Google Scholar] [CrossRef]
  52. Fernandez-Garcia, Y.; Horst, S.T.; Bassetto, M.; Brancale, A.; Neyts, J.; Rogolino, D.; Sechi, M.; Carcelli, M.; Gunther, S.; Rocha-Pereira, J. Diketo acids inhibit the cap-snatching endonuclease of several Bunyavirales. Antivir. Res. 2020, 183, 104947. [Google Scholar] [CrossRef]
  53. Kuang, W.; Zhang, H.; Cai, Y.; Zhang, G.; Deng, F.; Li, H.; Zhou, Y.; Wang, M.; Gong, P.; Guo, Y.; et al. Structural and Biochemical Basis for Development of Diketo Acid Inhibitors Targeting the Cap-Snatching Endonuclease of the Ebinur Lake Virus (Order: Bunyavirales). J. Virol. 2022, 96, e0217321. [Google Scholar] [CrossRef] [PubMed]
  54. Garcia-Serradilla, M.; Risco, C. Light and electron microscopy imaging unveils new aspects of the antiviral capacity of silver nanoparticles in bunyavirus-infected cells. Virus Res. 2021, 302, 198444. [Google Scholar] [CrossRef]
  55. Arif, H.; Aggarwal, S. Salicylic Acid (Aspirin). In StatPearls; StatPearls Publishing: Treasure Island, FL, USA, 2022. [Google Scholar]
Viruses 15 00948 g001 550
Figure 1. Molecular modeling and summary of the numeric outcomes. (Left panel): poses adopted by the selected compounds at the endonuclease catalytic site (drugs are represented as colored wireframes and target is sketched in green cartoons). Main residues are indicated by arrows. (Right panel): Computed free energies for the drug–endonuclease complexes (in kcal/mol).
Figure 1. Molecular modeling and summary of the numeric outcomes. (Left panel): poses adopted by the selected compounds at the endonuclease catalytic site (drugs are represented as colored wireframes and target is sketched in green cartoons). Main residues are indicated by arrows. (Right panel): Computed free energies for the drug–endonuclease complexes (in kcal/mol).
Viruses 15 00948 g001
Viruses 15 00948 g002 550
Figure 2. Antiviral activity of ASA against BUNV. Vero cells were absorbed 1 h with BUNV at an MOI of 1 PFU/cell, exposed to increasing concentrations of ASA for 8 h and processed by immunofluorescence. (A–D) Fluorescence microscopy of infected cells untreated (A) and treated with 5, 10 and 15 mM ASA (B–D). Nuclei are labeled with DAPI (blue) and viral infection with an antibody specific for BUNV nucleoprotein (N) (green). (E) Percentage of infected cells measured by immunofluorescence in the absence of ASA treatment and with three different concentrations of ASA. Data are shown as mean ± s.e.m. of three independent experiments. (F) Dose–response curve (red line) of ASA was determined by nonlinear regression. Cytotoxic effect on Vero cells exposed to increasing concentrations of the drug in the absence of virus is also shown (black line). IC50 and CC50 were determined from dose–response curves based on treatment with seven different concentrations. IC50 was calculated from the percentage of infected cells, and CC50 was determined by the MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide) cell viability assay. ND, no drug. Scale bars: 50 μm.
Figure 2. Antiviral activity of ASA against BUNV. Vero cells were absorbed 1 h with BUNV at an MOI of 1 PFU/cell, exposed to increasing concentrations of ASA for 8 h and processed by immunofluorescence. (A–D) Fluorescence microscopy of infected cells untreated (A) and treated with 5, 10 and 15 mM ASA (B–D). Nuclei are labeled with DAPI (blue) and viral infection with an antibody specific for BUNV nucleoprotein (N) (green). (E) Percentage of infected cells measured by immunofluorescence in the absence of ASA treatment and with three different concentrations of ASA. Data are shown as mean ± s.e.m. of three independent experiments. (F) Dose–response curve (red line) of ASA was determined by nonlinear regression. Cytotoxic effect on Vero cells exposed to increasing concentrations of the drug in the absence of virus is also shown (black line). IC50 and CC50 were determined from dose–response curves based on treatment with seven different concentrations. IC50 was calculated from the percentage of infected cells, and CC50 was determined by the MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide) cell viability assay. ND, no drug. Scale bars: 50 μm.
Viruses 15 00948 g002
Viruses 15 00948 g003 550
Figure 3. Study of ASA antiviral effect by plaque reduction assay. BUNV was propagated in Vero cells in the absence or presence of different concentrations of ASA (5, 10 and 15 mM) for 8 h, the supernatants were collected and the virus titer was tested by plaque assay. (A) A representative picture showing the plaque reduction assay with the drug. (B) BUNV titers from cell supernatants without ASA and with three different concentrations of ASA. Each bar shows the mean value ± s.e.m. Significance was determined by unpaired two-tailed student’s t-test; ns indicates “not significant”, ** p < 0.01 and *** p < 0.001.
Figure 3. Study of ASA antiviral effect by plaque reduction assay. BUNV was propagated in Vero cells in the absence or presence of different concentrations of ASA (5, 10 and 15 mM) for 8 h, the supernatants were collected and the virus titer was tested by plaque assay. (A) A representative picture showing the plaque reduction assay with the drug. (B) BUNV titers from cell supernatants without ASA and with three different concentrations of ASA. Each bar shows the mean value ± s.e.m. Significance was determined by unpaired two-tailed student’s t-test; ns indicates “not significant”, ** p < 0.01 and *** p < 0.001.
Viruses 15 00948 g003
Viruses 15 00948 g004 550
Figure 4. Effect of ASA on viral protein expression. Non-infected and BUNV-infected Vero cells at an MOI of 1 PFU/cell for 1 h and incubated in the presence or absence of ASA (5, 10 and 15 mM) for 8 h, were processed for Western blot analysis. An amount of 10 micrograms of the cell lysates were loaded and samples were probed with an antiserum against BUNV. The tubulin signal was used as loading control. NI, non-infected; ND, no drug. One experiment representative of 3 is shown. The tables show the average relative quantity of the Gc protein and N protein ± standard deviation (SD) for all the experiments. The p-values (one- and two-tailed, respectively) were calculated with the student’s t-test.
Figure 4. Effect of ASA on viral protein expression. Non-infected and BUNV-infected Vero cells at an MOI of 1 PFU/cell for 1 h and incubated in the presence or absence of ASA (5, 10 and 15 mM) for 8 h, were processed for Western blot analysis. An amount of 10 micrograms of the cell lysates were loaded and samples were probed with an antiserum against BUNV. The tubulin signal was used as loading control. NI, non-infected; ND, no drug. One experiment representative of 3 is shown. The tables show the average relative quantity of the Gc protein and N protein ± standard deviation (SD) for all the experiments. The p-values (one- and two-tailed, respectively) were calculated with the student’s t-test.
Viruses 15 00948 g004
Viruses 15 00948 g005 550
Figure 5. ASA treatment induces the accumulation of BUNV Gc glycoprotein in the cell perinuclear region. Vero cells were absorbed 1 h with BUNV at an MOI of 1 PFU/cell, incubated in the presence or absence of ASA (5, 10 and 15 mM) for 8 h and processed by immunofluorescence. (A–F) Confocal fluorescence microscopy of infected cells untreated and treated with ASA. The nucleus is labeled with DAPI (blue), the BUNV Gc glycoprotein (green) with a specific antibody, and the Golgi complex with the anti-Giantin antibody (red). In untreated cells (A), Gc signal is distributed throughout the cytoplasm, while Gc signal in ASA-treated cells (B–D) is retained in a perinuclear region associated with the Golgi complex. (E,F) Three-dimensional reconstructions from a series of confocal images. In the absence of the drug (E), Gc signal is dispersed, while in cells treated with 15 mM ASA (F), Gc signal accumulates near the nucleus in the Golgi complex (white arrows). (G) Average of the number of cells with both N and Gc viral proteins-associated fluorescent signals in the absence and presence of ASA measured by immunofluorescence. This average was determined by counting the number of cells with both N and Gc protein signals. (H) Quantification of the distribution of Gc IF signal shows that ASA induces a retention of Gc signal in the perinuclear region. Quantification was performed by counting the number of cells with Gc protein signal in the cytosol or in the perinuclear region. A total of 50 cells per condition were studied. The results in panels G and H are presented as the mean of the number of cells with viral proteins-associated fluorescence signals ± s.e.m. of three independent experiments. Significance was determined by unpaired two-tailed student’s t-test; ns indicates “not significant”, ** p < 0.01 and *** p < 0.001. ND, no drug. Scale bars: 25 μm in (A–D); 10 μm in (E,F).
Figure 5. ASA treatment induces the accumulation of BUNV Gc glycoprotein in the cell perinuclear region. Vero cells were absorbed 1 h with BUNV at an MOI of 1 PFU/cell, incubated in the presence or absence of ASA (5, 10 and 15 mM) for 8 h and processed by immunofluorescence. (A–F) Confocal fluorescence microscopy of infected cells untreated and treated with ASA. The nucleus is labeled with DAPI (blue), the BUNV Gc glycoprotein (green) with a specific antibody, and the Golgi complex with the anti-Giantin antibody (red). In untreated cells (A), Gc signal is distributed throughout the cytoplasm, while Gc signal in ASA-treated cells (B–D) is retained in a perinuclear region associated with the Golgi complex. (E,F) Three-dimensional reconstructions from a series of confocal images. In the absence of the drug (E), Gc signal is dispersed, while in cells treated with 15 mM ASA (F), Gc signal accumulates near the nucleus in the Golgi complex (white arrows). (G) Average of the number of cells with both N and Gc viral proteins-associated fluorescent signals in the absence and presence of ASA measured by immunofluorescence. This average was determined by counting the number of cells with both N and Gc protein signals. (H) Quantification of the distribution of Gc IF signal shows that ASA induces a retention of Gc signal in the perinuclear region. Quantification was performed by counting the number of cells with Gc protein signal in the cytosol or in the perinuclear region. A total of 50 cells per condition were studied. The results in panels G and H are presented as the mean of the number of cells with viral proteins-associated fluorescence signals ± s.e.m. of three independent experiments. Significance was determined by unpaired two-tailed student’s t-test; ns indicates “not significant”, ** p < 0.01 and *** p < 0.001. ND, no drug. Scale bars: 25 μm in (A–D); 10 μm in (E,F).
Viruses 15 00948 g005
Viruses 15 00948 g006 550
Figure 6. Transmission electron microscopy of BUNV-infected cells untreated or treated with ASA. Vero cells were absorbed 1 h with BUNV at an MOI of 1 PFU/cell, incubated in the presence or absence of ASA (5, 10 and 15 mM) for 8 h and processed by transmission electron microscopy. (A–C) BUNV-infected Vero cells not treated with ASA (no drug or ND). (A) Recruitment of cellular organelles to a discrete area of the cytoplasm forming the viral factory (VF). The inset shows extracellular virions. (B) Higher magnification view of the VF shows the Golgi complex (G) surrounded by mitochondria (mi) and endoplasmic reticulum (ER) membranes. Mitochondria have swollen cristae. (C) Bunyavirus replication organelles (spherules and tubules, arrowheads) in Golgi membranes. The inset shows an intracellular virus in the Golgi complex. (D–F) BUNV-infected Vero cells treated with ASA 5 mM. The cellular organelles are not recruited to a specific area into the cytoplasm (D). Golgi complex stacks surrounded by endoplasmic reticulum membranes and modified mitochondria (E). The arrowhead in (F) points to a viral spherule in Golgi membranes. (G–L) BUNV-infected Vero cells treated with ASA 10 or 15 mM. Low magnification images (G,J) show cellular organelles distributed through the cytoplasm not forming a characteristic VF. Golgi membranes are swollen without bunyavirus spherules and mitochondria have normal morphology and shape (H,K). Details of modified Golgi complex (asterisks) are shown in (I,L). Neither BUNV particles nor spherules were detected. N, nucleus. Scale bars: 2 μm in (A,D,G,J); 500 nm in (B,E,H,K); 200 nm in the inset of (A,C,F,I,L); 100 nm in the inset of (C).
Figure 6. Transmission electron microscopy of BUNV-infected cells untreated or treated with ASA. Vero cells were absorbed 1 h with BUNV at an MOI of 1 PFU/cell, incubated in the presence or absence of ASA (5, 10 and 15 mM) for 8 h and processed by transmission electron microscopy. (A–C) BUNV-infected Vero cells not treated with ASA (no drug or ND). (A) Recruitment of cellular organelles to a discrete area of the cytoplasm forming the viral factory (VF). The inset shows extracellular virions. (B) Higher magnification view of the VF shows the Golgi complex (G) surrounded by mitochondria (mi) and endoplasmic reticulum (ER) membranes. Mitochondria have swollen cristae. (C) Bunyavirus replication organelles (spherules and tubules, arrowheads) in Golgi membranes. The inset shows an intracellular virus in the Golgi complex. (D–F) BUNV-infected Vero cells treated with ASA 5 mM. The cellular organelles are not recruited to a specific area into the cytoplasm (D). Golgi complex stacks surrounded by endoplasmic reticulum membranes and modified mitochondria (E). The arrowhead in (F) points to a viral spherule in Golgi membranes. (G–L) BUNV-infected Vero cells treated with ASA 10 or 15 mM. Low magnification images (G,J) show cellular organelles distributed through the cytoplasm not forming a characteristic VF. Golgi membranes are swollen without bunyavirus spherules and mitochondria have normal morphology and shape (H,K). Details of modified Golgi complex (asterisks) are shown in (I,L). Neither BUNV particles nor spherules were detected. N, nucleus. Scale bars: 2 μm in (A,D,G,J); 500 nm in (B,E,H,K); 200 nm in the inset of (A,C,F,I,L); 100 nm in the inset of (C).
Viruses 15 00948 g006
Viruses 15 00948 g007 550
Figure 7. Quantification of spherules, intracellular and extracellular viruses in BUNV-infected cells, with or without ASA using transmission electron microscopy. Vero cells were absorbed 1 h with BUNV at an MOI of 1 PFU/cell, incubated in the presence or absence of ASA (5, 10 and 15 mM) for 8 h and processed by embedding in epoxy resin, ultramicrotomy and electron microscopy. The table shows the absolute and relative frequencies of spherules, intracellular and extracellular viruses counted in BUNV-infected cells in each condition: BUNV ND (infected with no ASA), BUNV + ASA 5 mM, BUNV + ASA 10 mM and BUNV + ASA 15 mM. ASA treatment produces a significant reduction of the number of viral structures. The graphs showing the means of spherules, intracellular and extracellular viruses, counted in a total of 20 cells per condition, are shown near each table, on the right.
Figure 7. Quantification of spherules, intracellular and extracellular viruses in BUNV-infected cells, with or without ASA using transmission electron microscopy. Vero cells were absorbed 1 h with BUNV at an MOI of 1 PFU/cell, incubated in the presence or absence of ASA (5, 10 and 15 mM) for 8 h and processed by embedding in epoxy resin, ultramicrotomy and electron microscopy. The table shows the absolute and relative frequencies of spherules, intracellular and extracellular viruses counted in BUNV-infected cells in each condition: BUNV ND (infected with no ASA), BUNV + ASA 5 mM, BUNV + ASA 10 mM and BUNV + ASA 15 mM. ASA treatment produces a significant reduction of the number of viral structures. The graphs showing the means of spherules, intracellular and extracellular viruses, counted in a total of 20 cells per condition, are shown near each table, on the right.
Viruses 15 00948 g007
Viruses 15 00948 g008 550
Figure 8. Recovery of a productive BUNV infection after elimination of ASA from cell cultures. Vero cells were absorbed 1 h with BUNV at an MOI of 1 PFU/cell and treated with ASA 10 mM for 8 h. The medium was removed and replaced by new DMEM with 2% FBS. After incubation at different time points (2, 4, 6 and 8 h), cells were processed for immunofluorescence, titration by plaque assay and electron microscopy. (A) The percentage of infected cells per condition was determined by fluorescence microscopy using an antibody specific for BUNV N protein. After removing the drug from the cell culture medium, the number of infected cells increased over time. (B) Graph showing the viral titer at different time points. Each bar shows the mean value ± s.e.m. Significance was determined by unpaired two-tailed student’s t-test, where ** indicates p < 0.01. (C–H) Electron microscopy study. After 8 h of ASA treatment (C), the Golgi membranes were swollen. Neither BUNV spherules nor viral particles were detected in these infected and ASA-treated cells. (D–G) After ASA elimination, BUNV particles (arrows) and spherules (arrowheads) are seen associated with Golgi membranes. (H) High magnification image of a Golgi complex with a viral spherule (arrowhead). N, nucleus; G, Golgi complex. Scale bars: 500 nm in (C–G); 200 nm in (H).
Figure 8. Recovery of a productive BUNV infection after elimination of ASA from cell cultures. Vero cells were absorbed 1 h with BUNV at an MOI of 1 PFU/cell and treated with ASA 10 mM for 8 h. The medium was removed and replaced by new DMEM with 2% FBS. After incubation at different time points (2, 4, 6 and 8 h), cells were processed for immunofluorescence, titration by plaque assay and electron microscopy. (A) The percentage of infected cells per condition was determined by fluorescence microscopy using an antibody specific for BUNV N protein. After removing the drug from the cell culture medium, the number of infected cells increased over time. (B) Graph showing the viral titer at different time points. Each bar shows the mean value ± s.e.m. Significance was determined by unpaired two-tailed student’s t-test, where ** indicates p < 0.01. (C–H) Electron microscopy study. After 8 h of ASA treatment (C), the Golgi membranes were swollen. Neither BUNV spherules nor viral particles were detected in these infected and ASA-treated cells. (D–G) After ASA elimination, BUNV particles (arrows) and spherules (arrowheads) are seen associated with Golgi membranes. (H) High magnification image of a Golgi complex with a viral spherule (arrowhead). N, nucleus; G, Golgi complex. Scale bars: 500 nm in (C–G); 200 nm in (H).
Viruses 15 00948 g008
Table
Table 1. Examples of members of the Bunyavirales order and associated diseases that could potentially be treated with ASA, assuming that studies with animal models and clinical trials will give good results, and those that could not be treated with ASA because they produce hemorrhagic fever. Distribution, reservoirs and vectors are included. N.I., not identified. Sources: CDC, NIH and ECDC.
Table 1. Examples of members of the Bunyavirales order and associated diseases that could potentially be treated with ASA, assuming that studies with animal models and clinical trials will give good results, and those that could not be treated with ASA because they produce hemorrhagic fever. Distribution, reservoirs and vectors are included. N.I., not identified. Sources: CDC, NIH and ECDC.
VirusDiseaseDistributionReservoirVectorASA
La CrosseChildhood encephalitisNorth, Central and South AmericaSmall mammals (chipmunks and squirrels)MosquitoesYes
Cache Valley
virus		Fever, encephalitis, meningitis North and
Central America		Deer, cattle, horses and sheepN.I.Yes
Toscana virusMeningitis, encephalitisMediterranean region of Western Europe (Portugal, Spain, France, Greece and Croatia), Cyprus and TurkeyPossibly migratory birds and domestic animalsSandfliesYes
Uukuniemi virusFever, headache, muscle and joint painScandinaviaN.I.TicksYes
SchmallenbergCongenital malformations and abortion in livestockBelgium, Germany, France, Italy, Luxemburg,
the Netherlands, Spain and United Kingdom		N.I.Possibly mosquitoes and midgesYes
RVFVAcute and febrile illness with few severe cases (neurological disorders, partial or complete blindness, hemorrhagic fever or thrombosis)
High mortality and abortion in domestic ruminants		Africa,
Saudi Arabia and Yemen		Livestock (cattle, sheep, goats, buffalos and camels)MosquitoesNo
CCHFVHemorrhagic feverEastern and Southern Europe, Northwestern China, Africa, the Middle East and the Indian subcontinentWild and domestic
animals		TicksNo
Lassa Fever virusHemorrhagic feverEndemic in parts of West Africa including Sierra Leone, Liberia, Guinea and NigeriaMultimammate ratMultimammate ratNo
SFTSV
(Dabie bandavirus)		Hemorrhagic fever
Severe fever with thrombocytopenia		Northwest and central China, Japan, South Korea, Vietnam and TaiwanN.I.TicksNo
HantaanHemorrhagic fever with renal syndrome (Old World)
Hantavirus pulmonary syndrome (New World)		Europe and Asia
North and South America		RodentsRodentsNo
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Fernández-Sánchez, S.Y.; Cerón-Carrasco, J.P.; Risco, C.; Fernández de Castro, I. Antiviral Activity of Acetylsalicylic Acid against Bunyamwera Virus in Cell Culture. Viruses 2023, 15, 948. https://doi.org/10.3390/v15040948

AMA Style

Fernández-Sánchez SY, Cerón-Carrasco JP, Risco C, Fernández de Castro I. Antiviral Activity of Acetylsalicylic Acid against Bunyamwera Virus in Cell Culture. Viruses. 2023; 15(4):948. https://doi.org/10.3390/v15040948

Chicago/Turabian Style

Fernández-Sánchez, Sara Yolanda, José P. Cerón-Carrasco, Cristina Risco, and Isabel Fernández de Castro. 2023. "Antiviral Activity of Acetylsalicylic Acid against Bunyamwera Virus in Cell Culture" Viruses 15, no. 4: 948. https://doi.org/10.3390/v15040948

APA Style

Fernández-Sánchez, S. Y., Cerón-Carrasco, J. P., Risco, C., & Fernández de Castro, I. (2023). Antiviral Activity of Acetylsalicylic Acid against Bunyamwera Virus in Cell Culture. Viruses, 15(4), 948. https://doi.org/10.3390/v15040948

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop