- 1Division of Animal Biotechnology, Faculty of Veterinary Sciences and Animal Husbandry, Sher-e- Kashmir University of Agricultural Sciences and Technology of Kashmir, Srinagar, India
- 2Department of Biotechnology, University of Kashmir, Srinagar, India
- 3Department of Biotechnology, School of Engineering and Technology, Sharda University, Greater Noida, UP, India
- 4Faculty of Animal Science and Veterinary Medicine, Sher-e-Bangla Agricultural University, Dhaka, Bangladesh
- 5Division of Computational Biology, Daskdan Innovations Pvt. Ltd., Srinagar, India
- 6College of Veterinary Medicine, Jeonbuk National University, Iksan, Republic of Korea
Introduction: The antiviral activity of different mutagens against single-stranded RNA viruses is well documented; however, their activity on the replication of double-stranded RNA viruses remains unexplored. This study aims to investigate the effect of different antivirals on the replication of a chicken embryo fibroblast-adapted Infectious Bursal Disease virus, FVSKG2. This study further explores the antiviral mechanism utilized by the most effective anti-IBDV agent.
Methods: The cytotoxicity and anti-FVSKG2 activity of different antiviral agents (ribavirin, 5-fluorouracil, 5-azacytidine, and amiloride) were evaluated. The virus was serially passaged in chicken embryo fibroblasts 11 times at sub-cytotoxic concentrations of ribavirin, 5-fluorouracil or amiloride. Further, the possible mutagenic and non-mutagenic mechanisms utilized by the most effective anti-FVSKG2 agent were explored.
Results and Discussion: Ribavirin was the least cytotoxic on chicken embryo fibroblasts, followed by 5-fluorouracil, amiloride and 5-azacytidine. Ribavirin inhibited the replication of FVSKG2 in chicken embryo fibroblasts significantly at concentrations as low as 0.05 mM. The extinction of FVSKG2 was achieved during serial passage of the virus in chicken embryo fibroblasts at ≥0.05 mM ribavirin; however, the emergence of a mutagen-resistant virus was not observed until the eleventh passage. Further, no mutation was observed in 1,898 nucleotides of the FVSKG2 following its five passages in chicken embryo fibroblasts in the presence of 0.025 mM ribavirin. Ribavarin inhibited the FVSKG2 replication in chicken embryo fibroblasts primarily through IMPDH-mediated depletion of the Guanosine Triphosphate pool of cells. However, other mechanisms like ribavirin-mediated cytokine induction or possible inhibition of viral RNA-dependent RNA polymerase through its interaction with the enzyme’s active sites enhance the anti-IBDV effect. Ribavirin inhibits ds- RNA viruses, likely through IMPDH inhibition and not mutagenesis. The inhibitory effect may, however, be augmented by other non-mutagenic mechanisms, like induction of antiviral cytokines in chicken embryo fibroblasts or interaction of ribavirin with the active sites of RNA-dependent RNA polymerase of the virus.
1. Introduction
RNA viruses are turning out to be the major causes of global pandemics in recent times. The growing human-animal contact can lead to the evolution of certain viruses and generate diseases with the potential to spread between species and cause pandemics (1, 2). Numerous mutagenic nucleoside analogues have already been approved for use after thorough investigations of their antiviral properties on single-stranded RNA viruses; however, the assessment of these nucleoside analogues against double-stranded RNA viruses is yet to be explored.
Infectious Bursal Disease Virus (IBDV) is a non-enveloped, double-stranded (ds)-RNA virus belonging to the genus Avibirnavirus and family Birnaviridae, that causes a substantial economic impact on the poultry industries globally (3–10). The virus possesses a bi-segmented genome containing segment-A (3.2 kb) that encodes viral proteins; VP2, VP3, VP4, and VP-5 (11–15) and Segment-B (2.8 kb) encoding VP1, which is an RNA-dependent RNA polymerase (RdRp) (16). VP1 contains important active sites necessary for IBDV replication. The active site of the polymerase is formed by Asp-402, Asp-416, and Asn-403. At the same time, Ser-166 amino acid offers a self-guanylylation activity to the virus, which is essential for protein priming, an alternate mechanism of genomic replication of the virus (17).
Ribavirin (1-β-D-ribofuranosyl-2,2,4-triazole-3-carboxamide) is a classic example of base-analogue mutagens clinically used against various viruses affecting human health. Lethal mutagenesis is a primary antiviral mechanism utilized by ribavirin against ssRNA viruses through an accumulation of mutations in the virus beyond a threshold (18–20). However, few studies have explored various non-mutagenic mechanisms utilized by ribavirin to inhibit viral replication (21–27). Although ribavirin has been used clinically as an antiviral for a long time and has been the subject of several studies, the mechanism underlying its antiviral activity is still not entirely understood (28–31). Moreover, previous studies have explored the possible antiviral mechanisms of ribavirin using ssRNA viruses (18, 20, 32–34); however, their effect on the replication of dsRNA viruses remains largely unexplored. The current study employs various strategies to evaluate the inhibitory potential of ribavirin against a dsRNA virus, IBDV, and further explore its possible mechanisms of action using a combination of in vitro and comprehensive in silico approaches.
2. Materials and methods
2.1. Viruses and antiviral mutagens
FVSKG2 (Accession number: OP161172), a local IBDV isolate, was used in the study. Four mutagens viz.; ribavirin, 5-fluorouracil, 5-azacytidine, and amiloride (Sigma, St. Louis, Missouri, United States) were used in this study. Each mutagen was dissolved separately in Dulbecco’s modified Eagle media (DMEM, Sigma, St. Louis, Missouri, United States) at stock concentrations of 15 mM (Ribavirin), 20 mM (5-fluorouracil, amiloride) and 5 mM (5-azacytidine) which were further sterile-filtered using a 0.22-μm syringe filter (Millex-GV Filter, 0.22 μm Millipore Sigma, Burlington; Massachusetts; United States). The filtrate was aliquoted and stored at −20°C until use, as described in previous studies (35, 36).
2.2. Chicken embryo fibroblast cell culture
The chicken embryo fibroblasts (CEFs) were isolated from 9 days old embryonated chicken eggs and maintained in growth media containing DMEM, 5% Fetal Bovine Serum (Fetal Bovine Serum, Sigma, St. Louis, Missouri, United States) and 1% Penicillin–Streptomycin Solution with 10,000 U penicillin and 10,000 μg streptomycin/mL (Pen Strep, Gibco Waltham, Massachusetts, United States), at 37°C and 5% CO2 in a humidified chamber as described in our previous studies (37, 38).
2.3. Cytotoxicity assay
The cytotoxic effect of mutagens on CEFs was determined using MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide) assay (Thiazolyl Blue Tetrazolium Bromide, Sigma, St. Louis, Missouri, United States). Confluent CEFs were prepared in 96-well plates (Tissue culture plate, 96-well Falcon, Flowery Branch, Georgia, United States). Cells were rinsed and replenished with 100 μL growth medium containing one of the four concentrations of mutagens (0, 0.5, 1.0, and 1.5 mM) in three replicates. To determine cell viability at specific time point, cells were incubated at 37°C and 5% CO2, as described in the previous study (35, 39). At 24 hpt and 48 hpt the growth medium of the cells was replaced with 100 μL of DMEM containing MTT at a concentration of 0.5 mg/mL and incubated for 4 h at 37°C. After the MTT treatment, 100 μL of the solubilizing solution containing 10% SDS in 0.01 M HCL was added and the cells were incubated overnight. Absorbance was read at 570 nm using a multimode microplate reader (Biotek, Cytation™ 3, Winooski, VT, United States) (40, 41). Control wells (cells with 0.0 mM mutagen) and blank wells (without cells) were utilized for calculating the cell viability by using the following formula as described in our earlier studies (37, 38).
2.4. Evaluation of the effect of mutagens on IBDV replication
CEFs were cultured in 25 cm2 flasks prior to virus inoculation. At the time of virus inoculation, cell number was determined and multiplicity of infection (MOI) of 0.1 was used. The anti-IBDV activity of ribavirin was evaluated against IBDV isolate FVSKG2. Confluent cultures of CEFs were infected with FVSKG2 and incubated for 2 h in a humidified incubator at 37°C and 5% CO2 following which the virus inoculum was discarded, and the cells were replenished with 5 mL maintenance media, each containing a specific concentration of mutagens. Different sub-cytotoxic concentrations of mutagens, ribavirin, 5-Azacytidine, and Amiloride at 0.0, 0.05, 0.1, 0.2, and 0.3 mM, 5-Florouracil at 0.0, 0.1, 0.2, 0.3, and 0.5 mM were tested against FVSKG2 based on previous studies (31, 35, 42, 43). The virus-inoculated cells, each treated with a specific concentration of mutagens, were then incubated for four more days under the same culture conditions as described above, during which time 500 μL of cell culture medium was collected and replaced with fresh media from each flask every 24 h. The collected media was centrifuged and the supernatant was stored at −80°C until further analysis.
2.5. Virus titration assay and determination of TCID50
Confluent CEFs in 96-well plates (Tissue culture plate, 96 well, Falcon, Flowery Branch, Georgia, United States) were used to evaluate the titer of FVSKG2 samples. Eight-fold serial dilution (10–1 to 10–8) of viruses in DMEM were prepared and inoculated in triplicates with 100 μL media per well and incubated for 2 h under cultural conditions. After incubation, the inoculum was discarded and cells were replenished with growth media. The cells were then incubated for 96 h and monitored for cytopathic effects (CPE). The wells were scored for CPE after 96 h. Virus titer was expressed as 50% tissue culture infective dose (TCID50/mL) and calculated using the method of Reed and Muench (44).
2.6. Serial passage of FVSKG2
FVSKG2 was serially passaged in CEFs in different concentrations of mutagens, Ribavirin and Amiloride at 0.0, 0.05, 0.1, 0.2, and 0.3 mM, and 5-Fluorouracil at 0.0, 0.1, 0.2, 0.3, and 0.5 mM for 11 passages. Confluent monolayers of CEFs, prepared in 6-well plates, were pre-treated with mutagens at given concentrations 4hrs prior to inoculation of a virus at 0.01 MOI. After 2 h post incubation, virus inoculum was removed, and cells were replenished with a growth medium containing the same concentrations of mutagens as used in the pre-treatment stage in each well. The infection was then allowed to proceed for 24 h, following which the plates were freeze-thawed thrice. Cell culture fluid was collected from each well and centrifuged at 2,500 rpm for 5 min, and the supernatant was stored at −80°C for virus titration. 200 μL of supernatant from each passage was used as virus inoculum for the next passage. This procedure was repeated for each serial passage of the virus. In addition, the FVSKG2 was serially passaged in the presence of 0.025 mM ribavirin to evaluate for the presence of mutations compared to the control (0 mM Ribavirin).
2.7. Viral RNA isolation and sequencing
According to the manufacturer’s instructions (Virus RNA isolation kit, GeneAll, South Korea), 300 μL of the supernatant of the stock virus of FVSKG2 was used to isolate the RNA, with the following changes. The VL buffer was used for 30 min of incubation with the viral sample. The Takara one-step RT-PCR kit was used to perform the one-step RT-PCR (Primescript One-step RT-PCR kit, Takara, Japan). Amplification was performed according to the manufacturer’s instructions, following some modifications suited for amplifying a dsRNA virus. In a 1:4 ratio, DMSO (Sigma Aldrich, ST. Louis, United States) was added to the RNA template with 0.5 μL RNAse inhibitor (rRNasin, Promega, Wisconsin, United States) and incubated at 99°C for 3 min and snap-chilled on ice following which primers, forward A_124F:CGCAGCGATGACAAACCT; reverse A_1100R: GATCCCCCGCCTGACCACCACTT for 976 bp region of segment A and forward B_1080F: CTGAAAGGTACGACAAAAGCACAT; reverse B_2002R: TACCAACCTCAACGCCTCATACCT for 922 bp region in segment B were added separately. Then 12.5 μL of 2X RT-PCR buffer and 1 μL enzyme mix were added to set up a 25 μL reaction in thermo-cycler (Biometra T Advanced, Analytika Jena, Jena, Germany). PCR products of size 976 bp and 922 bp, respectively, were amplified by using the following thermal cycling conditions 50°C for 30 min, 95°C for 2 min, followed by 35 cycles at 95°C for 30 s, 62°C and 61°C, respectively, for 1.5 min, 72°C for 1.5 min, and final extension of 72°C for 5 min. PCR products were purified by DNA purification kit (Wizard® SV Gel and PCR Clean-Up System, Promega Madison, Wisconsin, United States) and sent for sequencing to Macrogen (Seoul, South Korea).
2.8. Ribavirin-guanosine inhibition study
To evaluate whether the addition of guanosine may compete with ribavirin and rescue viral replication, FVSKG2 was added at MOI of 0.01 in CEFs. Ribavirin and guanosine combinations were then added to achieve final concentrations of 0, 0.05, 0.05 mM and 0, 0, 0.025 mM, respectively. Culture supernatants were harvested at 48 Hours Post Inoculation (HPI) to determine viral titers as described in previous studies (35, 45).
2.9. Inhibition of inosine-5′-monophosphate dehydrogenase
Confluent CEFs in 6-well plates were pre-treated with ribavirin at the concentration of 0, 10 and 40 μM in the presence or absence of a fixed concentration of 40 μM guanosine. Separately, confluent CEFs in 6-well plates were pre-treated with mycophenolic acid at 0, 1, 5, and 10 μM in the presence or absence of fixed concentration of 40 μM guanosine. FVSKG2 at MOI of 0.01 was added to the cells and incubated for 2 h. Inoculum was then removed and the cells were replenished with the same concentrations of the ribavirin, guanosine and mycophenolic acid as in pre-treatment, incubated at 37°C and 5% CO2. The viral titer was determined every 48 h post-infection.
2.10. Evaluation of mRNA expression of cytokines
Ribavirin was added at 0, 0.05 and 0.1 mM concentrations in CEFs cultured in 24-well plates. RNA was isolated from the ribavirin-treated CEFs 24-h post-treatment using an RNA isolation kit (GeneAll Hybrid-RTM kit, GeneAll Biotechnology, Seoul, South Korea) following the manufacturer instructions. RNA was reverse-transcribed into complementary DNA (cDNA) using a high-capacity cDNA reverse transcription kit (Thermofisher Scientific, Maasachusetts, United States) as per the following conditions: 25°C for 10 min, 37°C for 120 min and 85°C for 5 min and holding at 4°C. Real-time PCR was performed on Analytik Jena, qRT-PCR (Jena, Germany) system using various cytokine-specific primers (Table 1), following the manufacturer’s instructions. The qPCR reaction was set up using 10 μL of GoTaq qPCR Master Mix(2X), 0.5 μL of Forward Primer (20X), 0.5 μL of Reverse Primer(20X), 7 μL of Nuclease-Free Water, 2 μL of cDNA template (or water for the no-template control reactions) The cycling conditions for performing the qPCR were as follows: GoTaq® Hot Start Polymerase activation for 1 cycle at 95°C for 2 min, denaturation at 95°C for 15 sec, annealing and data collection (40 cycles) at temperature 60°C for 30 sec. The relative quantities of cytokine mRNA in ribavirin-treated and non-treated CEFs were normalized to GAPDH mRNA, and the amounts were determined using the 2-ΔΔCt method (46).
2.11. Protein and ligand structure preparation
Two proteins, RdRp (Uniprot ID: Q9Q6Q5) and IMPDH2 (Uniprot ID: Q5F4A4), and two ligands ribavirin (PubChem ID:37542) and mycophenolic acid (Pubchem ID: 446541) were used in this study. The Crystal Structure of RdRp (PDB ID: 2PGG) (17) was retrieved and used for docking and simulation studies. For IMPDH2, the protein sequence of gallus gallus of 514 residues was obtained from UniProtKB1 database. The sequence shared 95% similarity with Human IMPDH2 (PDB ID: P12268), and to generate the three-dimensional structure of IMPDH2 of gallus gallus, the structure representing the IMPDH2 with mycophenolic acid (PDB ID: 1JR1) (47) was used as a template for homology modeling. The modeling was carried out using the SWISS-MODEL server2 (48). The generated model was subjected to energy minimization and refinement procedures via the ModRefiner server3 (49). The energy-minimized structure was named IMPDH2_GG and was used throughout the study.
2.12. Molecular docking
AutoDock Tool (50) was used to perform highly extensive molecular docking. Ribavarin was docked into the crystal structure of RdRp (2PGG), and two separate dockings were performed on two different sites. Site 1 represented the vital amino acids viz. ASN402, ASP403, and ASN416 in the catalytic palm of RdRp. Site 2 was built around SER166, a self-guanylation site present in RdRp. For IMPDH2_GG, the potential binding sites were identified based on the structural comparison with human IMPDH2. The crystal structure of IMPDH2 with mycophenolic acid (1JR1) was used as the template to define the binding pocket for IMPDH_GG. Ribavirin and mycophenolic acid were docked into this binding pocket. All the complexes generated were subjected to the molecular dynamics simulation run for further analysis.
2.13. All atom molecular dynamics simulation
GROMACS 2021 series version-2 (51) molecular dynamics package was used to carry out an all-atom molecular dynamics simulation of the four complexes generated from molecular docking. The CHARMM36 (52) force field was used to define the complex of protein, water, and ions in the TIP3P water model. The ligands used, viz. ribavirin, and mycophenolic acid were processed in the CHARMM General Force Field (CGenFF) program4 (53). Energy Minimization and conjugate gradient algorithms of GROMACS were employed for optimizing the final protein-ligand complexes. The protein, ligand, water, and ions systems were equilibrated in NVT and NPT ensembles for 100 ps. Using the Parrinello–Rahman barostat, these ensembles maintained the system at 310 K temperature and 1 bar of pressure. Each production run included three replicas of 100 ns. Trajectories were saved every 2 ps/frame for further analysis.
2.14. Data analysis
The effect of mutagens on CEFs was analysed by Mann–Whitney U test while the effect of mutagen on FVSKG2 replication, serial passages and multi-step growth curve was analyzed by Repeated measures analysis of variance (ANOVA) followed by Dunnet’s post-test using GraphPad Prism 8 (DNASTAR (Inc., Madison, WI, United States). Nucleotide sequences were aligned and analyzed using DNASTAR (Inc., Madison, WI, United States). Cytokine expression was analyzed by the Mann–Whitney U test (SPSS Version 20) (35, 54). This study performed the primary docking using AutoDock Vina (The Scripps Research Institute, United States) and molecular dynamics simulation using the GROMACS-2021. Discovery Studio Biovia 2021 (Dassault Systèmes, USA) and PyMOL (Schrödinger, LLC) were employed for visualizations.
3. Results
3.1. Ribavirin exhibited the least cytotoxicity on CEFs
The cytotoxicity of the four mutagens was evaluated and is summarized in the Figure 1. No significant cytotoxicity was observed after 48 h post-treatment at concentration of 0.25 mM and 0.5 mM, but at 1.0 and 1.5 mM, ribavirin was found to be significantly cytotoxic. 5-fluorouracil, 5-azacytidine and amiloride showed significant cytotoxicity at 0.5 mM concentration and above.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1192583/fvets-10-1192583-HTML-r1/image_m/fvets-10-1192583-g001.jpg)
Figure 1. Cytotoxicity assay of mutagens on CEFs: Cell viability was determined by using MTT assay: Y-axis represents percent cell viability, while X-axis represents hours post-treatment (hpt). Data expressed as mean ± standard error (n = 3). Significant cytotoxicity was observed at and above 1.0 mM of ribavirin, while as in case of 5-fluorouracil, 5-azacytidine and amiloride significant cytotoxicity was observed at 0.5 mM. Asterisks represent significant difference (p < 0.05). The number of asterisks “*” represents the extent of significance. The 0 hpt refers to time period immediately after the treatment.
3.2. Ribavirin inhibited the IBDV replication in CEFs in a dose-dependent manner
The antiviral effect of ribavirin, 5- fluorouracil, 5-azacytidine and amiloride on replication of FVSKG2 in CEFs was evaluated using a multi-step growth curve. This was done by determining the titer (TCID50/mL) of mutagen-treated viral populations at indicated concentrations relative to the control concentration (0 mM) of mutagens, summarized in Figure 2. The replication of FVSKG2 decreased significantly at 0.05, 0.1, 0.2, and 0.3 mM of ribavirin in a dose-dependent manner (p < 0.0001). When FVSKG2 was treated with 0.3 mM ribavirin, there was about a 3.5-log10 reduction in virus titer at 72 h post treatment, followed by a 4.5-log10 reduction at 120 hpt. Although similar activity was observed with 0.2 mM ribavirin, a 4-log10 reduction and 2.5-log10 reduction were observed in the presence of 0.1 mM or 0.05 mM ribavirin, respectively. Significant suppression of FVSKG2 replication was observed at 0.5 mM of 5-fluorouracil with 2-log10 reduction in virus titer (p < 0.0001). However, no significant antiviral activity was observed with 0.1, 0.2, and 0.3 mM 5-fluorouracil. Treatment with 0.2 and 0.3 mM 5-azacytidine significantly achieved up to 2-log10 reduction in FVSKG2 replication (p = 0.0019). Similarly, high concentrations (1 and 2 mM) of amiloride also caused significant suppression of FVSKG2 replication with up to 2.5-log10 reduction in virus titer (p = 0.0266). Further, no significant antiviral activity was measured at low concentrations of 5-azacytidine (<0.1 mM) and amiloride (<0.05 mM).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1192583/fvets-10-1192583-HTML-r1/image_m/fvets-10-1192583-g002.jpg)
Figure 2. Effect of ribavirin, 5- fluorouracil, 5-azacytidine, and amiloride on replication of FVSKG2: The effect of indicated concentrations of the mutagens on the replication of IBDV isolate FVSKG2 in CEFs were evaluated. Results are presented as 50% tissue culture infective dose TCID50 per mL (y-axis) over time (x-axis). Virus titers were evaluated using cell culture fluids collected from pre-treated CEFs every 24 h after being incubated with the FVSKG2. At 0.05, 0.1, 0.2, and 0.3 mM of ribavirin, the replication of FVSKG2 significantly decreased in a dose-dependent manner. At 0.5 mM of 5-fluorouracil, a 2 log10 drop in virus titer was seen along with significant suppression of FVSKG2 replication (p = 0.0001). On the other hand, there was no detectable antiviral action with 0.1, 0.2, and 0.3 mM 5-fluorouracil. Treatment with 5-azacytidine at 0.2 and 0.3 mM significantly reduced FVSKG2 replication by up to 2 log10 (p = 0.0019). Similarly, significant inhibition of FVSKG2 replication was caused by high concentrations (1 and 2 mM) of amiloride. Error bars denote mean ± standard error of the mean (SEM) while Asterisks “*” indicate significant differences in virus titer (TCID50/mL log10) as compared to control (0 mM), p < 0.05. The number of asterisks represents the extent of significance.
3.3. Ribavirin drove an early extinction of IBDV in CEFs during serial passages
Based on the cytotoxicity and antiviral activity of the selected antiviral mutagens, serial passages of FVSKG2 in CEF cells were carried out in the presence of an indicated concentration of each mutagen (Figure 3). We observed that among the mutagens, ribavirin was active against FVSKG2, reducing the viral titer in a concentration-dependent manner. When passaged in the presence of 0.05 mM ribavirin, virus titers dropped for the first two passages causing a significant 4-log10 reduction at p1 and p2. However, the virus disappeared in the third passage and did not emerge until the eleventh passage. Furthermore, FVSKG2 replication was completely suppressed at higher ribavirin concentrations (>0.1 mM). The virus appeared to be less sensitive to 5-fluorouracil and amiloride at indicated concentrations, with no significant change in the virus titer across the 11 passages relative to the control (0 mM). The only exception was the 0.1 mM 5-fluorouracil, which exhibited significant fluctuations in virus titers (p < 0.05). In the presence of 5-fluorouracil at a concentration of 0.1 mM, the virus initially experienced a 1-log10 reduction in titer at passage 2, reaching an equilibrium from passage 4 to passage 6 with 105.5 TCID50/mL virus titer. By passage 8, the virus titer increased again to 105.5 TCID50/mL and reached a new equilibrium by passage 9 with a 1.5 log10 reduction in virus titer. The virus titer persisted at 105 TCID50/mL from passage 9 to passage 11.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1192583/fvets-10-1192583-HTML-r1/image_m/fvets-10-1192583-g003.jpg)
Figure 3. Sequential passage of FVSKG2 in different concentrations of ribavirin, 5-fluorouracil, and amiloride in CEFs: In presence of ribavirin at concentrations of 0.05 and above 0.1 mM, extinction of virus occurred during third and first passage, respectively. In presence of 5-fluorouracil no significant inhibition of FVSKG2 was observed at concentrations at and above 0.1 mM in CEFs. In presence of amiloride, no significant inhibition of FVSKG2 was observed at concentration at and above 0.05 mM in CEFs, Asterisks “*” indicate significant differences.
3.4. Ribavirin did not induce any mutations in FVSKG2 at 0.025 mM
FVSKG2 replication was inhibited to an undetectable level at 0.05 mM of ribavirin in CEFs. However, at 0.025 mM ribavirin, FVSKG2 was replicated at a low titer. No mutation was observed in 1,898 nucleotide bases in FVSKG2 passaged five times in the presence of 0.025 mM ribavirin in CEFs when compared to FVSKG2 passaged for the same number of times in the absence of ribavirin.
3.5. Ribavirin competition assay with guanosine
At 48 h after treatment, ribavirin inhibited the replication of FVSKG2, resulting in a viral titer of 102 TCID50/mL and a 5.5-log10 reduction in viral titer compared to 0 mM ribavirin. However, when 0.025 mM guanosine was added in addition to ribavirin, the viral titer increased to 106 TCID50/mL at 48 h. Therefore, Ribavirin inhibited FVSKG2 replication, while adding guanosine to the culture revived the virus replication (Figure 4).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1192583/fvets-10-1192583-HTML-r1/image_m/fvets-10-1192583-g004.jpg)
Figure 4. Ribavirin competition assay with guanosine: Ribavirin (Rib) competition assay with nucleoside guanosine (Guo) was performed with FVSKG2 at MOI of 0.01 in CEFs, under three different conditions: no drug treatment, 0.05 mM of ribavirin and 0.05 mM ribavirin plus 0.025 mM guanosine. Culture supernatants were harvested at 48 h post-treatment to determine viral titers. The mean viral titers ± s.d. from triplicates derived from one out of two independent experiments are shown.
3.6. Inhibition of IMPDH by ribavirin
The replication of FVSKG2 in CEFs was investigated at indicated concentrations of ribavirin and mycophenolic acid in the presence or absence of 40 mM guanosine. In the non-guanosine treated control, low titers of FVSKG2 have measured at 10 uM ribavirin, reaching a peak titer of almost 103.5 TCID50/mL at 48 h post-infection. However, no virus was detected at 40 μm ribavirin in cell-free supernatant. Guanosine addition increased the titers significantly, with titers reaching >103.5 TCID50/mL at 48 h post-infection. At 10 μm ribavirin, guanosine addition caused a 1-log10 increase in virus titer and completely rescued the virus at 40 μm ribavirin with titer reaching 104.5 TCID50/mL. Consistent with this, similar findings were observed in case of mycophenolic acid, wherein guanosine supplementation repressed the antiviral effect of the mycophenolic acid. Intriguingly, mycophenolic acid inhibited the FVSKG2 replication at the indicated concentration (>1 μM). The addition of guanosine resulted in a 4.5 log10 increase in virus titer at 1 μM mycophenolic acid, which persisted to higher concentrations of mycophenolic acid. The virus achieved a 3.5-log10 increase in titer at 10 μM mycophenolic acid, reaching about 103.5 TCID50/mL at 48 h post-infection compared to non-guanosine treated mycophenolic acid control (Figure 5A).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1192583/fvets-10-1192583-HTML-r1/image_m/fvets-10-1192583-g005.jpg)
Figure 5. (A) The effect of indicated concentrations of ribavirin (Rib) and mycophenolic acid (MPA) in presence or absence of Guanosine (Guo) Supplementation on the replication of FVSKG2 in CEFs. Virus titers are presented as 50% tissue culture infective dose (TCID50/mL, log10) per milliliter (y axis) over drug concentration (x axis). Errors bars denote mean ± standard error of mean (SEM). (B) B-factor diagram indicating binding free-energy contribution of interacting residues. (C) FEL values constructed as a function of PC1 and PC2 eigenvectors. (D) The ligand-IMPDH interaction diagram showing the presence of important interactions.
Molecular docking of ribavirin and mycophenolic Acid with IMPDH was performed using Autodock vina. As a potent IMPDH inhibitor, mycophenolic acid was chosen as the standard reference molecule. Analysis of the molecular interactions summarized in Table 2 showed negative binding energy for both molecules within the IMPDH binding pocket. While Mycophenolic Acid (MPA) exhibited a binding energy of −5.9 Kcal/mol, ribavirin with IMPDH exhibited a binding energy of −8.1 Kcal/mol. The molecular interactions revealed the presence of Van der wall interactions as well as hydrogen bonding, indicating efficient ligand binding (Figure 5D). Overall, these studies suggested stable binding of ribavirin compared to the reference molecule within the IMPDH binding pocket (S1 Movie). Ribavirin can thereby prevent viral replication by acting as a potent inhibitor of IMPDH.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1192583/fvets-10-1192583-HTML-r1/image_m/fvets-10-1192583-t002.jpg)
Table 2. Molecular Docking Presentation of binding affinity and binding pockets of Mycophenolic Acid and Ribavirin with IMPDH.
The docked complexes were subjected to molecular dynamics simulation and MMPBSA analysis using GROMACS and gmx_MMPBSA, respectively. The residual decomposition energy of the interacting residues was calculated at an interval of 20 ps, representing 500 frames over the 100 ns trajectories. The per residue energy decomposition plot showed stable and persistent binding of ribavirin to IMPDH compared to mycophenolic acid (Supplementary Figure S1). In addition, the energy2bfac tool used to map the binding energy contribution of interacting residues illustrated the involvement of a larger number of residues in the binding of ribavirin within the IMPDH binding pocket (Figure 5B). The Free Energy Landscape (FEL) values of Rib-IMPDH and MPA-IMPDH were constructed and plotted (Figure 5C). A comparison of the FEL values for the complexes revealed that Rib-IMPDH FEL spanned significantly larger areas of PC1 and PC2 with more free energy wells in the region of the global free energy minimum region. These findings indicate greater flexibility and conformational diversity of Rib-IMPDH compared to the MPA-IMPDH complex.
3.7. Effect of ribavirin on the expression of cytokines in chicken embryo fibroblasts
Real-time PCR was used to assess the cytokine levels in CEFs following a 24-h treatment with 0.05 or 0.1 mM ribavirin. As compared to mock-treated CEFs, the expression of IFN-α (Figure 6A), TNF-α (Figure 6B), IL-2, IL-12 (Figure 6C), and IL-10 (Figure 6D) was significantly higher in CEFs treated either with 0.05 or 0.1 mM ribavirin. However, there was no significant difference in the expression of IFN-β between mock-treated and ribavirin-treated CEFs. Further, the expression of IL-6 in CEFs treated with 0.05 or 0.1 mM of ribavirin was significantly lower than in the mock-treated CEFs (Figure 6B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1192583/fvets-10-1192583-HTML-r1/image_m/fvets-10-1192583-g006.jpg)
Figure 6. Evaluation of cytokines alteration in CEF cells by ribavirin: CEFs were grown in the presence or absence of indicated concentrations of ribavirin. CEFs were harvested at 24 hpt and subjected to cytokine mRNA expression analysis. Cytokines expression analysis of IL-6, IL-2, IL-12, interferon-alpha (IFN-α), tumor necrosis factor-alpha (TNF-α) and TNF-β mRNAs in CEF cells after 24-h culture. qRT-PCR was used for evaluating mRNA expression. While Relative quantification (RQ) was calculated using the 2-DDCt method. Error bars denote mean ± standard error of the mean (SEM). The bars represent the means, and the error bars represent the standard errors of the mean (SEM). Bars showing different letters represent values that differ significantly from each other (p < 0.05).
3.8. In silico analysis of ribavirin as inhibitor of RdRp of IBDV
Ribavirin shows a brief interaction with the two active sites of RdRp (S2 Movie). While site-1, formed of ASP 402, ASN 403, and ASP 416, constitutes the conserved catalytic active site, site-2, formed by SER 166, constitutes the RdRp self-guanylation site (Figure 7A). Ribavirin shows decent binding energy within both target sites. The analysis of the docking results revealed binding energy of −4.39 kcal/mol at site 2 and − 5.51 kcal/mol at site-1 (Figure 7B). The presence of Van der wall interactions and hydrogen bonding indicate efficient binding of ribavirin at both the target sites (Figure 7C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1192583/fvets-10-1192583-HTML-r1/image_m/fvets-10-1192583-g007.jpg)
Figure 7. Interaction of ribavirin with RdRp: (A) Denotion of the target sites in the IBDV RdRp structure. (B) Autodock vina based Molecular docking analysis of ribavirin at site 1 and site 2. (C) Analysis of interacting residues and nature of interaction between ribavirin and IBDV RdRp target sites.
4. Discussion
In the current study, ribavirin was found to be the least cytotoxic on CEFs among the tested mutagens, which agrees with the previous studies that used different cells (18, 35, 55). Ribavirin significantly inhibited the replication of FVSKG2 at concentrations of 0.05 mM and above which is in agreement with previous studies where the antiviral activity of ribavirin has been demonstrated against many other RNA viruses like influenza virus, poliovirus, human respiratory syncytial virus, hantaan virus and foot and mouth disease virus (FMDV) (18, 20, 32–34).
The mechanism of action of ribavirin on viral replication is not fully understood despite its use as an antiviral (28–30). In the current study, the treatment of FVSKG2 with 0.025 mM ribavirin for series of five passages could not induce any mutations in the 1,898 base pairs of viral genome sequenced, which is in disagreement with previous studies where antiviral activity of ribavirin on ssRNA viruses, viz. FMDV, HCV, polio virus, hepatitis E virus and mumps virus was found to be due to mutagenic action of ribavirin (20, 32, 39, 56, 57). The reasons for no-mutagenic action of ribavirin on FVSKG2 may be due to variable reasons: (a) dsRNA are reported to be less prone to mutations than ssRNA viruses (58), (b) Sequencing only a part of viral genome, and (c) using sanger’s sequencing method which cannot detect the sub-consensus variants that could have likely generated during the treatment of virus with 0.025 mM ribavirin.
In the current study, ribavirin inhibited IBDV replication through guanosine depletion which is in agreement with previous studies where ribavirin was found to deplete guanosine pool of cells (18). Further, in the current study, the ribavirin-mediated depletion of the guanosine pool of cells was demonstrated to occur due to inhibition of IMPDH by ribavirin (Figure 5), which was further validated by the strong interaction of ribavirin with IMPDH shown by the in-silico studies. This is in agreement with previous studies where the antiviral effect of ribavirin on flaviviruses, paramyxoviruses, dengue virus and influenza A virus was found to be due to depletion of cellular guanosine pool via IMPDH inhibition (26, 59–62).
Ribavirin showed brief interaction with two active sites of RdRp of IBDV in the current study (Figure 6) which is in agreement with previous studies where ribavirin was observed to possess a weak inhibitory activity on RdRp of many RNA viruses viz. HCV, Hepatitis E virus, bovine viral diarrhea virus, vesicular stomatitis virus, influenza virus, reovirus, and HIV in vitro or in silico (23, 24, 63).
We observed a significantly higher mRNA expression of IFN-α, TNF- α, IL-2, IL-12, and IL-10 was in CEFs treated with 0.05 or 0.1 mM ribavirin which is partly in agreement with previous studies where ribavirin was also found to up-regulate the interferon-stimulated genes and enhance the effect of interferons (21, 22). In another study, ribavirin was found to augment IL-2, IFN-α, and TNF-α and suppress IL-4 and IL-5 (24). Further, the expression of IL-6 in CEFs treated with 0.05 or 0.1 mM of ribavirin was significantly lower than in mock-treated CEFs which is in agreement with a study where ribavirin was shown to likely suppress IL-6 in alveolar epithelial cells (64). Collectively, ribavirin inhibited the replication of FVSKG2 predominantly through the ribavirin-mediated guanosine depletion of CEFs; however, other mechanisms like induction of antiviral cytokines in CEFs by ribavirin or the interaction of ribavirin with the active sites of viral RdRP could also have an additive effect on the inhibition of the virus replication.
5. Conclusion
The study reveals that dsRNA viruses are likely more resistant to ribavirin-mediated mutagenesis compared to the ssRNA viruses. Ribavirin inhibits dsRNA viruses through non-mutagenic mechanisms, primarily through depletion of the guanosine pool of cells which is likely augmented by its stimulation of cellular cytokines or its inhibition of active sites of viral RdRp.
Data availability statement
The data presented in the study are deposited in the NCBI repository (https://www.ncbi.nlm.nih.gov/), accession number OP161172.1.
Ethics statement
The animal study was reviewed and approved by Institutional Animal Ethics Committee (IAEC), Faculty of Veterinary Sciences and Animal Husbandry, Sher-e-Kashmir University of Agricultural Sciences and Technology of Kashmir (Registration no: 1809/GO/ReBi/15/CPCSEA).
Author contributions
TA and MP performed the in vitro experiments and drafted the manuscript. IG performed the part of in-vitro and in silico study and drafted the manuscript. W-IK and AK conducted thorough reviews and editing of the manuscript. AH performed additional in vitro experiments suggested by the reviewers. RAS, SMA, and NAG contributed to the reading and revision of the manuscript. NS and NC contributed to the conception and design of the study, the data analyses, and manuscript revisions. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by a grant from the Science and Engineering Research Board, Department of Science and Technology, Government of India, under the early career research award scheme (No: ECR/2016/001172).
Conflict of interest
NC is employed by Daskdan Innovations Pvt. Ltd.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2023.1192583/full#supplementary-material
Footnotes
References
1. Carrasco-Hernandez, R, Jácome, R, López Vidal, Y, and Ponce de León, S. Are RNA viruses candidate agents for the next global pandemic? A review. ILAR J. (2017) 58:343–58. doi: 10.1093/ilar/ilx026
2. Joseph, U, Su, YC, Vijaykrishna, D, and Smith, GJ. The ecology and adaptive evolution of influenza a interspecies transmission. Influenza Other Respir Viruses. (2017) 11:74–84. doi: 10.1111/irv.12412
3. Dey, S, Pathak, DC, Ramamurthy, N, Maity, HK, and Chellappa, MM. Infectious bursal disease virus in chickens: prevalence, impact, and management strategies. Vet Med (Auckl). (2019) 10:85–97. doi: 10.2147/VMRR.S185159
4. Sharma, JM, Kim, IJ, Rautenschlein, S, and Yeh, HY. Infectious bursal disease virus of chickens: pathogenesis and immunosuppression. Dev Comp Immunol. (2000) 24:223–35. doi: 10.1016/S0145-305X(99)00074-9
5. Withers, D, Young, J, and Fred, DT. Infectious bursal disease virus-induced immunosuppression in the Chick is associated with the presence of undifferentiated follicles in the recovering Bursa. Viral Immunol. (2005) 18:127–37. doi: 10.1089/vim.2005.18.127
6. Qi, X, Zhang, L, Chen, Y, Gao, L, Wu, G, Qin, L, et al. Mutations of residues 249 and 256 in VP2 are involved in the replication and virulence of infectious bursal disease virus. PLoS One. (2013) 8:e70982. doi: 10.1371/journal.pone.0070982
7. Ingrao, F, Rauw, F, Lambrecht, B, and van den Berg, T. Infectious bursal disease: a complex host-pathogen interaction. Dev Comp Immunol. (2013) 41:429–38. doi: 10.1016/j.dci.2013.03.017
8. Mahgoub, HA, Bailey, M, and Kaiser, P. An overview of infectious bursal disease. Arch Virol. (2012) 157:2047–57. doi: 10.1007/s00705-012-1377-9
9. Müller, H, Islam, MR, and Raue, R. Research on infectious bursal disease—the past, the present and the future. Vet Microbiol. (2003) 97:153–65. doi: 10.1016/j.vetmic.2003.08.005
11. Mundt, E, Beyer, J, and Müller, H. Identification of a novel viral protein in infectious bursal disease virus-infected cells. J Gen Virol. (1995) 76:437–43. doi: 10.1099/0022-1317-76-2-437
12. Muller, H, and Becht, H. Biosynthesis of virus-specific proteins in cells infected with infectious bursal disease virus and their significance as structural elements for infectious virus and incomplete particles. J Virol. (1982) 44:384–92. doi: 10.1128/JVI.44.1.384-392.1982
13. Azad, AA, Barrett, SA, and Fahey, KJ. The characterization and molecular cloning of the double-stranded RNA genome of an Australian strain of infectious bursal disease virus. Virology. (1985) 143:35–44. doi: 10.1016/0042-6822(85)90094-7
14. Hudson, PJ, McKern, NM, Power, BE, and Azad, AA. Genomic structure of tbe large RNA segment of infectious bursal disease virus. Nucleic Acids Res. (1986) 14:5001–12. doi: 10.1093/nar/14.12.5001
15. Kibenge, FS, and Dhama, V. Evidence that virion-associated VP1 of avibirnaviruses contains viral RNA sequences. Arch Virol. (1997) 142:1227–36.
16. Muller, H, and Nitschke, R. The two segments of the infectious bursal disease virus genome are circularized by a 90,000-Da protein. Virology. (1987) 159:174–7.
17. Pan, J, Vakharia, VN, and Tao, YJ. The structure of a birnavirus polymerase reveals a distinct active site topology. Proc Natl Acad Sci U S A. (2007) 104:7385–90. doi: 10.1073/pnas.0611599104
18. Pauly, MD, and Lauring, AS. Effective lethal mutagenesis of influenza virus by three nucleoside analogs. J Virol. (2015) 89:3584–97. doi: 10.1128/JVI.03483-14
19. Anderson, JP, Daifuku, R, and Loeb, LA. Viral error catastrophe by mutagenic nucleosides. Annu Rev Microbiol. (2004) 58:183–205. doi: 10.1146/annurev.micro.58.030603.123649
20. Airaksinen, A, Pariente, N, Menéndez-Arias, L, and Domingo, E. Curing of foot-and-mouth disease virus from persistently infected cells by ribavirin involves enhanced mutagenesis. Virology. (2003) 311:339–49. doi: 10.1016/S0042-6822(03)00144-2
21. Feld, JJ, Lutchman, GA, Heller, T, Hara, K, Pfeiffer, JK, Leff, RD, et al. Ribavirin improves early responses to peginterferon through improved interferon signaling. Gastroenterology. (2010) 139:154–62 e4. doi: 10.1053/j.gastro.2010.03.037
22. Thomas, E, Feld, JJ, Li, Q, Hu, Z, Fried, MW, and Liang, TJ. Ribavirin potentiates interferon action by augmenting interferon-stimulated gene induction in hepatitis C virus cell culture models. Hepatology. (2011) 53:32–41. doi: 10.1002/hep.23985
23. Eriksson, B, Helgstrand, E, Johansson, NG, Larsson, A, Misiorny, A, Noren, JO, et al. Inhibition of influenza virus ribonucleic acid polymerase by ribavirin triphosphate. Antimicrob Agents Chemother. (1977) 11:946–51. doi: 10.1128/AAC.11.6.946
24. Lau, JY, Tam, RC, Liang, TJ, and Hong, Z. Mechanism of action of ribavirin in the combination treatment of chronic HCV infection. Hepatology. (2002) 35:1002–9. doi: 10.1053/jhep.2002.32672
25. Zhou, S, Liu, R, Baroudy, BM, Malcolm, BA, and Reyes, GR. The effect of ribavirin and IMPDH inhibitors on hepatitis C virus subgenomic replicon RNA. Virology. (2003) 310:333–42. doi: 10.1016/S0042-6822(03)00152-1
26. Malinoski, F, and Stollar, V. Inhibitors of IMP dehydrogenase prevent Sindbis virus replication and reduce GTP levels in Aedes albopictus cells. Virology. (1981) 110:281–91. doi: 10.1016/0042-6822(81)90060-X
27. Streeter, DG, Witkowski, JT, Khare, GP, Sidwell, RW, Bauer, RJ, Robins, RK, et al. Mechanism of action of 1- -D-ribofuranosyl-1,2,4-triazole-3-carboxamide (Virazole), a new broad-spectrum antiviral agent. Proc Natl Acad Sci U S A. (1973) 70:1174–8. doi: 10.1073/pnas.70.4.1174
28. Crotty, S, Maag, D, Arnold, JJ, Zhong, W, Lau, JY, Hong, Z, et al. The broad-spectrum antiviral ribonucleoside ribavirin is an RNA virus mutagen. Nat Med. (2000) 6:1375–9. doi: 10.1038/82191
29. Jain, MK, and Zoellner, C. Role of ribavirin in HCV treatment response: now and in the future. Expert Opin Pharmacother. (2010) 11:673–83. doi: 10.1517/14656560903580001
30. Witkowski, JT, Robins, RK, Sidwell, RW, and Simon, LN. Design, synthesis, and broad spectrum antiviral activity of 1- -D-ribofuranosyl-1,2,4-triazole-3-carboxamide and related nucleosides. J Med Chem. (1972) 15:1150–4. doi: 10.1021/jm00281a014
31. Hadj Hassine, I, Ben M’hadheb, M, and Menéndez-Arias, L. Lethal mutagenesis of RNA viruses and approved drugs with antiviral mutagenic activity. Viruses. (2022) 14:841. doi: 10.3390/v14040841
32. Crotty, S, Cameron, CE, and Andino, R. RNA virus error catastrophe: direct molecular test by using ribavirin. PNAS Nexus. (2001) 98:6895–900. doi: 10.1073/pnas.111085598
33. Aljabr, W, Touzelet, O, Pollakis, G, Wu, W, Munday, DC, Hughes, M, et al. Investigating the influence of ribavirin on human respiratory syncytial virus RNA synthesis by using a high-resolution transcriptome sequencing approach. J Virol. (2016) 90:4876–88. doi: 10.1128/JVI.02349-15
34. Mayor, J, Engler, O, and Rothenberger, S. Antiviral efficacy of ribavirin and Favipiravir against Hantaan virus. Microorganisms. (2021) 9:1306. doi: 10.3390/microorganisms9061306
35. Khatun, A, Shabir, N, Yoon, KJ, and Kim, WI. Effects of ribavirin on the replication and genetic stability of porcine reproductive and respiratory syndrome virus. BMC Vet Res. (2015) 11:21. doi: 10.1186/s12917-015-0330-z
36. Sun, Y, Song, M, Niu, L, Bai, X, Sun, N, Zhao, X, et al. Antiviral effects of the constituents derived from Chinese herb medicines on infectious bursal disease virus. Pharm Biol. (2013) 51:1137–43. doi: 10.3109/13880209.2013.781197
37. Sofi, HS, Akram, T, Tamboli, AH, Majeed, A, Shabir, N, and Sheikh, FA. Novel lavender oil and silver nanoparticles simultaneously loaded onto polyurethane nanofibers for wound-healing applications. Int J Pharm. (2019) 569:118590. doi: 10.1016/j.ijpharm.2019.118590
38. Ashraf, R, Sofi, HS, Akram, T, Rather, HA, Abdal-hay, A, Shabir, N, et al. Fabrication of multifunctional cellulose/TiO2/ag composite nanofibers scaffold with antibacterial and bioactivity properties for future tissue engineering applications. J Biomed Mater Res A. (2020) 108:947–62. doi: 10.1002/jbm.a.36872
39. Jurković, M, Slović, A, Forčić, D, Ivančić-Jelečki, J, Košutić-Gulija, T, and Jagušić, M. Influence of ribavirin on mumps virus population diversity. Viruses. (2021) 13:2535. doi: 10.3390/v13122535
40. Musser, JMB, Heatley, JJ, Koinis, AV, Suchodolski, PF, Guo, J, Escandon, P, et al. Ribavirin inhibits parrot Bornavirus 4 replication in cell culture. PLoS One. (2015) 10:e0134080. doi: 10.1371/journal.pone.0134080
41. van Meerloo, J, Kaspers, GJ, and Cloos, J. Cell sensitivity assays: the MTT assay. Methods Mol Biol. (2011) 731:237–45.
42. Bouchard, J, Walker, MC, Leclerc, JM, Lapointe, N, Beaulieu, R, and Thibodeau, L. 5-azacytidine and 5-azadeoxycytidine inhibit human immunodeficiency virus type 1 replication in vitro. ASM J CD. (1990) 34:206–9. doi: 10.1128/AAC.34.2.206
43. Dapp, MJ, Clouser, CL, Patterson, S, and Mansky, LM. 5-Azacytidine can induce lethal mutagenesis in human immunodeficiency virus type 1. J Virol. (2009) 83:11950–8. doi: 10.1128/JVI.01406-09
44. Reed, LJ, and Muench, H. A SIMPLE METHOD OF ESTIMATING FIFTY PER CENT ENDPOINTS12. Am J Epidemiol. (1938) 27:493–7. doi: 10.1093/oxfordjournals.aje.a118408
45. Akram, T, Muneeb, JM, Gul, I, Parveez, M, Shah, RA, Ahmad, SM, et al. Isolation of infectious bursal disease virus from a poultry farm in Srinagar district of Jammu and Kashmir. SKUAST J Res. (2020) 22:86–91.
46. Shabir, N, Khatun, A, Nazki, S, Kim, B, Choi, E-J, Sun, D, et al. Evaluation of the cross-protective efficacy of a chimeric porcine reproductive and respiratory syndrome virus constructed based on two field strains. Viruses. (2016) 8:240. doi: 10.3390/v8080240
47. Sintchak, MD, Fleming, MA, Futer, O, Raybuck, SA, Chambers, SP, Caron, PR, et al. Structure and mechanism of inosine monophosphate dehydrogenase in complex with the immunosuppressant mycophenolic acid. Cells. (1996) 85:921–30. doi: 10.1016/S0092-8674(00)81275-1
48. Waterhouse, A, Bertoni, M, Bienert, S, Studer, G, Tauriello, G, Gumienny, R, et al. SWISS-MODEL: homology modelling of protein structures and complexes. Nucleic Acids Res. (2018) 46:W296–303. doi: 10.1093/nar/gky427
49. Zhang, B, Zhang, X, Pearce, R, Shen, H-B, and Zhang, Y. A new protocol for atomic-level protein structure modeling and refinement using low-to-medium resolution Cryo-EM density maps. J Mol Biol. (2020) 432:5365–77. doi: 10.1016/j.jmb.2020.07.027
50. Morris, GM, Huey, R, Lindstrom, W, Sanner, MF, Belew, RK, Goodsell, DS, et al. AutoDock4 and AutoDockTools4: automated docking with selective receptor flexibility. J Comput Chem. (2009) 30:2785–91. doi: 10.1002/jcc.21256
51. Berendsen, HJ, van der Spoel, D, and van Drunen, R. GROMACS: a message-passing parallel molecular dynamics implementation. Comput Phys Commun. (1995) 91:43–56. doi: 10.1016/0010-4655(95)00042-E
52. Gutiérrez, IS, Lin, F-Y, Vanommeslaeghe, K, Lemkul, JA, Armacost, KA, Brooks, CL III, et al. Parametrization of halogen bonds in the CHARMM general force field: improved treatment of ligand–protein interactions. Bioorg Med Chem. (2016) 24:4812–25. doi: 10.1016/j.bmc.2016.06.034
53. Vanommeslaeghe, K, Raman, EP, and MacKerell, AD Jr. Automation of the CHARMM general force field (CGenFF) II: assignment of bonded parameters and partial atomic charges. J Chem Inf Model. (2012) 52:3155–68. doi: 10.1021/ci3003649
54. Khatun, A, Shabir, N, Seo, BJ, Kim, BS, Yoon, KJ, and Kim, WI. The attenuation phenotype of a ribavirin-resistant porcine reproductive and respiratory syndrome virus is maintained during sequential passages in pigs. J Virol. (2016) 90:4454–68. doi: 10.1128/JVI.02836-15
55. Unal, MA, Bitirim, CV, Summak, GY, Bereketoglu, S, Cevher Zeytin, I, Besbinar, O, et al. Ribavirin shows antiviral activity against SARS-CoV-2 and downregulates the activity of TMPRSS2 and the expression of ACE2 in vitro. Can J Physiol Pharmacol. (2021) 99:449–60. doi: 10.1139/cjpp-2020-0734
56. Ortega-Prieto, AM, Sheldon, J, Grande-Perez, A, Tejero, H, Gregori, J, Quer, J, et al. Extinction of hepatitis C virus by ribavirin in hepatoma cells involves lethal mutagenesis. PLoS One. (2013) 8:e71039. doi: 10.1371/journal.pone.0071039
57. Todt, D, Gisa, A, Radonic, A, Nitsche, A, Behrendt, P, Suneetha, PV, et al. In vivo evidence for ribavirin-induced mutagenesis of the hepatitis E virus genome. Gut. (2016) 65:1733–43. doi: 10.1136/gutjnl-2015-311000
58. Sanjuán, R, and Domingo-Calap, P. Mechanisms of viral mutation. Cell Mol Life Sci. (2016) 73:4433–48. doi: 10.1007/s00018-016-2299-6
59. Streeter, DG, Witkowski, J, Khare, GP, Sidwell, RW, Bauer, RJ, Robins, RK, et al. Mechanism of action of 1-β-D-ribofuranosyl-1, 2, 4-triazole-3-carboxamide (Virazole), a new broad-spectrum antiviral agent. Proc Natl Acad Sci. (1973) 70:1174–8. doi: 10.1073/pnas.70.4.1174
60. Leyssen, P, Balzarini, J, Clercq, ED, and Neyts, J. The predominant mechanism by which ribavirin exerts its antiviral activity in vitro against Flaviviruses and paramyxoviruses is mediated by inhibition of IMP dehydrogenase. J Virol. (2005) 79:1943–7. doi: 10.1128/JVI.79.3.1943-1947.2005
61. Browne, MJ. Mechanism and specificity of action of ribavirin. Antimicrob Agents Chemother. (1979) 15:747–53. doi: 10.1128/AAC.15.6.747
62. Takhampunya, R, Ubol, S, Houng, HS, Cameron, CE, and Padmanabhan, R. Inhibition of dengue virus replication by mycophenolic acid and ribavirin. J Gen Virol. (2006) 87:1947–52. doi: 10.1099/vir.0.81655-0
63. Cancela, F, Rendon-Marin, S, Quintero-Gil, C, Houston, DR, Gumbis, G, Panzera, Y, et al. Modelling of hepatitis E virus RNA-dependent RNA polymerase genotype 3 from a chronic patient and in silico interaction analysis by molecular docking with ribavirin. J Biomol Struct Dyn. (2021) 3:1–17. doi: 10.1080/07391102.2021.2011416
Keywords: infectious bursal disease virus, ribavirin, dsRNA, antiviral, mutagen
Citation: Akram T, Gul I, Parveez Zia M, Hassan A, Khatun A, Shah RA, Ahmad SM, Ganai NA, Chikan NA, Kim W-I and Shabir N (2023) Ribavirin inhibits the replication of infectious bursal disease virus predominantly through depletion of cellular guanosine pool. Front. Vet. Sci. 10:1192583. doi: 10.3389/fvets.2023.1192583
Edited by:
Rajeev Ranjan, ICAR-National Institute on Foot and Mouth Disease, IndiaReviewed by:
Jinyang Zhang, Kunming University of Science and Technology, ChinaBasavaraj S. Mathapati, Indian Council of Medical Research (ICMR), India
Copyright © 2023 Akram, Gul, Parveez Zia, Hassan, Khatun, Shah, Ahmad, Ganai, Chikan, Kim and Shabir. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Won-Il Kim, kwi0621@jbnu.ac.kr; Nadeem Shabir, nadeem.shabir@skuastkashmir.ac.in