key: cord-1024727-ikn1btuk
authors: Shah, Fahad Hassan; Lim, Kyeong Ho; Kim, Song Ja
title: Do fever-relieving medicines have anti-COVID activity: an in silico insight
date: 2021-03-24
journal: Future Virology
DOI: 10.2217/fvl-2020-0398
sha: f0ac4f9d6eba6fc2b866bc39fe8e86332b5044b6
doc_id: 1024727
cord_uid: ikn1btuk

Aim: The present study was performed to determine the inhibitory interaction of fever-relieving medicines with severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) essential proteins. Materials & methods: Structure-based drug repositioning was performed using PYRX 0.9 and these drugs were directed toward the predicted active site of SARS-CoV-2 spike glycoprotein receptor-binding domain, main protease and RNA-dependent RNA polymerase. Results: Results showed that acetaminophen and naproxen have considerable inhibitory activity and show a high affinity for active residues of these proteins. The prediction of activity spectra for substances (PASS) studies showed that these drugs are anti-inflammatory, antiviral and immunostimulant. Conclusion: Hence, it is proven that these drugs have antiviral activity against SARS-CoV-2 and can stimulate the immune and anti-inflammatory response against this disease.

The pandemic storm led by novel coronaviruses infected billions and has claimed millions of lives worldwide. Infectious waves across different countries have begun at a rapid pace with increased lethality while therapeutics and vaccines are still awaiting approval from higher authorities. Coronaviruses are a group of highly contagious viruses that have their roots in order Nidovirales, family Coronaviridae and genus Beta coronaviruses [1] . Severe acute respiratory syndrome coronavirus-2, also known as 2019 novel coronavirus is the member of beta-coronaviruses, which are positive-sense RNA viruses enclosed within enveloped capsid having a genome of about 30 kilobase pairs [2] .

Upon infection with the host cell, the RNA genome of this virus directs the production of 27 viral proteins in which 16 nonstructural proteins are translated from 5 end. Whereas, the structural region encodes four integral proteins contiguous to one another near 3 end and are responsible for virus-host attachment (spike glycoprotein), capsid enveloping (envelope protein), membrane formation (M-protein) and genome encapsulation (nucleocapsid). The rest of the seven proteins called accessory proteins are translated through the process of ribosomal frameshifting. These viruses primarily affect the respiratory tract, disguising the immune system in order to proliferate viral copies [3] . Upon recognition, the immune system gets engaged in an exaggerated immunological response, also known as a cytokine storm that impairs the normal body physiological homeostasis, leading to death [4] .

These viruses primarily affect the respiratory system by employing spike glycoproteins to interact with ACE receptors on the host cells and injecting their genome to initiate RNA translation and replication. During the initial stages, an individual infected with this virus remains asymptomatic for a definite time period because of immunosuppression facilitated by nonstructural proteins [3] . Upon immunorecognition, individual experiences nonspecific clinical symptoms, in other words, body pain, chills, fever and coughing. As the disease progresses, headache, olfactory dysfunction, diarrhea and later the development of a cytokine storm, which damages pulmonary cells and leads to acute respiratory distress and multiorgan failure.

During the manifestation of these clinically important symptoms, usually, fever-relieving medicines are prescribed and administered to alleviate fever and myalgia. Most patients develop immunity against this disease with these medicines, whereas other patients require extra supportive care to eradicate this infection [5] . So, it was hypothesized that fever medicines not only tackle fever but also perform antiviral activity by improving the immune response and inhibiting SARS-CoV-2 replication [6] . To elucidate this unexplained biological activity of fever medicine, the current research study was designed to determine the anti-Covid activity of fever-relieving medicines using an in silico screening platform.

Critical SARS-CoV-2 proteins responsible for replication (SARS-Cov-2 RNA-dependent RNA polymerase (PDB ID: 6M71), viral-host interaction (SARS-CoV-2 receptor-binding domain (PDB ID: 7JMO) and viral polyprotein processing (SARS-CoV-2 main protease (PDB ID: 6LU7) were selected and procured from the RCSB protein databank (Supplementary Table 1 ). These protein receptors were initially prepared by removing pre-existing ligands and water ions from them through the Discovery Studio program and then minimized with Modrefiner [7] . Fever-alleviating drugs listed in (Table 1) were obtained from PubChem and refined with the PRODRG Server [8] .

In order to predict the probable inhibitory interaction of ligands with corresponding receptors, active site residues play an integral role in substrate processing and inhibition of protein function. The interaction of ligands with these active residues halt substrate processing, which results in protein inhibition. By keeping this strategy in mind, DogSite Scorer [9] and Prank Web [10] were used to predict the active site residues within these receptors to direct these ligands toward these areas in the screening analysis.

Drug repositioning procedure Both minimized and refined versions of receptors and ligands were loaded on PYRX 0.9 software to facilitate the interaction of these drugs with the predicted active site residues of SARS-CoV-2 receptors (Supplementary Table  2 ). The algorithm selected for docking interaction between receptor-ligands was Autodock 4.2, which provides highly accurate results resembling in vitro results.

Screened ligands SMILES (Supplementary Table 3 ) was used to predict important biological activities using the prediction of activity spectra for substances database (PASS) [11] . The pharmacokinetic prediction was performed to evaluate the efficacy and absorption, distribution, metabolism, excretion and toxicity (ADMET) qualities of these ligands through ADMET SAR 2.0 [12] and pKCSM [13] .

Structure-based drug repositioning was performed that uses docking algorithm to target ligands toward a directed area for interaction on the receptor. The selected ligands, (Table 1) were focused on the active interface of SARS-CoV-2 receptors. These receptors were spike glycoprotein receptor binding domain (RBD), main proteases (M PRO ) and RNA-dependent RNA polymerase (RdRp) enzyme implicated in the processing of viral polyproteins, virushost interaction and viral replication. Among seven widely exploited fever-relieving drugs, acetaminophen and naproxen had favorable hydrogen bonding with the targeted active site of these receptors with binding energy of ∼4.5-6.2 kcal/mol, as shown in (Table 2) .

Hydroxyl and amide group of acetaminophen formed hydrogen bonds with MET165, GLU166 and VAL186 (Figure 1 ), whereas MET165 and GLN189 participated in hydrophobic interactions with the benzene group. The carboxylic group of naproxen established interactions with GLY143, SER144 and CYS145, and hydrophobic interaction made by the benzene group with CYS145, and MET165 (Supplementary Figure 1) . 

Acetaminophen established three hydrogen bonds with PHE338, GLY339 and LEU368 whereas benzene ring participated in forming hydrophobic interactions with PHE338 and VAL367 (Supplementary Figure 2) . Naproxen employed carboxylic group to form hydrogen bonds with CYS336, PHE338 and ASP364 and benzene rings and oxy-group established hydrophobic interactions with VAL367, PHE338, LEU368 and PHE374 (Supplementary Figure 3) .

Three hydrogen interactions were formed by three separate functional groups (hydroxyl, amino and carboxyl) of acetaminophen with PHE45, HIS133 and ASN705. Naproxen exploited the carboxylic group to form two hydrogen bonds with GLN468 and ASN705 and benzene rings and oxy-group formed hydrophobic interactions with PHE45, LEU708, TYR129 and ALA34, respectively ( Supplementary Figures 4 & 5) .

The toxic behavior of acetaminophen and naproxen was determined in rodent models by administering various doses through four administration routes. It was observed that acetaminophen incites adverse and organ-damaging effects, when the quantity of the drug is increased to 2,000,000 mg/kg orally, 100,000 mg/kg intravenously and 700,000 mg/kg and 1,100,000 mg/kg for the intraperitoneal and subcutaneous route. For naproxen, acute toxicity recorded for intraperitoneal route was 500,000 mg/kg; 200,000 mg/kg for intravenous route; 900,000 mg/kg oral; and 600,000 mg/kg for subcutaneous. OECD chemical classification results placed acetaminophen in Class 5 chemicals and naproxen in Class 4, respectively. The organ-damaging effects (Supplementary Table 4 ) and adverse effects of these drugs are summarized in (Table 3) .

The activity of acetaminophen and naproxen was analyzed with thousands of biologically active substrates to discover new effects and mechanisms, which was facilitated by PASS online server. The biological activities predicted by PASS algorithm are 90 % synonymous with in vitro experiments [11] . Results showed that acetaminophen is an antiviral agent, active against influenza, picornavirus, adenovirus, coronavirus, HIV and herpes virus. This provides invaluable justification that acetaminophen is involved in subduing viral proliferation activity. Also, they are immunostimulant and inhibit the activity of transcription factor NF-κB. Naproxen is active against rhinoviruses, picornaviruses and influenza and dampens the response of NF-κB and IL-6, which are implicated in the cytokine storm ( Table 4 ).

Both of these drugs have high intestinal absorption equipped with increased blood-brain barrier (BBB) penetration capability and oral bioavailability, except for acetaminophen, which lacks oral bioavailability. Subcellular localization of these drugs are mitochondria and are nonreactive toward P-glycoprotein. Acetaminophen possesses affinity for CYP2C9 substrates and naproxen for CYP2C9 and inhibits CYPE1A2. These drugs are noncarcinogenic and nonmutagenic. Naproxen has low water solubility and high plasma protein binding, whereas acetaminophen has high water solubility and low plasma protein-binding affinity ( 

Asymptomatic disease nature, uncontrollable inflammatory response and delayed immunological reactions are typical manifestations of novel coronavirus infection. These manifestations are caused by nonstructural proteins processed by SARS-CoV-2 main proteases [3] , whereas replication and virus-cell interactions are facilitated by RdRp and spike protein [1, 14] . Currently, these targets are being exploited for vaccine and drug design and development studies. Several different drug candidates have been introduced from various screening, in vitro and clinical studies in the past few months. Drugs approved by the US FDA for the treatment of novel coronavirus in critical patients are remdesivir and barictinib [15] . However, their dosages must be calibrated for an effective response against this virus and prevent deleterious side effects. The mild symptomatic cases are dealt with fever relieving medicines to alleviate fever and other symptoms to buy enough time for the host body to develop an immune response to ward off this infection. Interestingly, 80% of patients recover from this disease by taking these fever relieving medicines [5] . Being analgesic, anti-inflammatory and antipyretic agent, these drugs are somehow acting as an antiviral agent in this disease, which remained elusive until now. To unravel this mystery, we have performed a drug repositioning study to check the inhibitory interaction of these fever-relieving medicines against SARS-CoV-2 essential proteins. It was observed that only two commonly exploited drugs (acetaminophen and naproxen) successfully established hydrogen interaction with the active site residues of these receptors with adequate binding energy, and ligand efficiency. Furthermore, root-mean-square deviation (RMSD) and inhibitory constant values favored these drugs. PASS prediction studies revealed that acetaminophen has moderate antiviral activity against different viruses including coronaviruses and low immunostimulant activity. On the other hand, naproxen has strong anti-inflammatory activity and avert the proliferation of rhinoviruses, picornaviruses and influenza viruses and inhibit NF-κB and IL-6 activity. These activities suggest that, in the early phase of SARS-CoV-2 infection, these drugs will perform antiviral activity by reducing viral entry by binding to the RBD of spike glycoprotein, modulate adaptive immune response and ameliorate viral clearance by dampening the viral polyprotein processing and replication. These drugs will also prevent SARS-CoV-2-induced damage to pulmonary epithelial cells and other organs implicated in the pulmonary phase. Our findings are in consonance with Terrier et al., 2020 [16] preprint, who infected Vero E6 cell and 3D reconstituted human airway epithelia with SARS-CoV-2 to evaluate the antiviral and pulmonary protective activity of naproxen. In their study, naproxen dose-dependently reduced viral titers in both cell lines and induced a protective effect against viral-induced damage. Naproxen also reduces IL-6 and NF-kB levels as evident from the study by Kang et al., 2001 [17] , which can be useful in reducing cytokine storm and organ damage during hyper-inflammation phase. However, these results are specific for naproxen, whereas the role of acetaminophen remains elusive in this disease, but can be explored by employing various in vitro experiments to confirm these in silico predicted activities. Toxicity and pharmacokinetic studies were performed to evaluate the toxicity and pharmacokinetic properties of these drugs, to equip scientists with valuable information about these drugs that will help them devise various approaches to amplify its antiviral and anti-inflammatory response in animal and clinical trials. Acetaminophen has high acute toxicity values with mild adverse effects but possesses low oral bioavailability. Naproxen has comparatively low acute toxicity than acetaminophen, with life-threatening side effects. Both toxic and bioavailability issues can be resolved with controlled drug-release properties [18] , drug-delivery nanocarriers [19] and drug-dose calibration [20] .

To the best of our knowledge, our study predicted the antiviral mechanism of acetaminophen and naproxen for the first time. These drugs are highly effective against essential SARS-CoV-2 protein targets and possess moderate toxicity and pharmacokinetic properties. However, more studies are essential to validate the findings of this research in animals and clinical trials. The findings of this study have opened a new and safer gateway for using these drugs as a prophylactic and antiviral medication that can be of great significance to clinicians and biomedical scientists struggling to curb the ongoing pandemic.

• SARS-CoV-2, also known as novel coronavirus-19 (COVID-19) is a highly contagious and life-threatening disease and employs spike glycoprotein, main proteases and RNA-dependent RNA polymerase enzyme to interact with angiotensin-converting enzyme-2 receptors on host cells, process viral polyproteins and replicate to produce viral copies of themselves inside the infected cell. • Most patients infected with this virus recover and develop immunity against this disease by taking fever relieving medicines and supportive care. • Structure-based drug repositioning studies confirmed that acetaminophen and naproxen act as antiviral, anti-inflammatory and immunostimulatory agents against this disease. • These drugs interact with active ligand binding site residues of spike glycoprotein, main protease and RdRp enzyme. • Further biological activities predicted by PASS indicate that these drugs are active against coronavirus and various other pathogenic viruses and can prevent the onset of cytokine storm. • Toxicity and pharmacokinetic properties of these drugs are cautionary but can be ameliorated with drug delivery, controlled drug release and drug dose calibration methods. • To the best of our knowledge, this study is the first of its kind to describe the antiviral activity of acetaminophen and naproxen against essential SARS-CoV-2 receptors.

To view the supplementary data that accompany this paper please visit the journal website at: www.futuremedicine.com/doi/sup pl/10.2217/fvl-2020-0398

All authors equally contributed for this study.

This work was supported by the National Research Foundation of Korea (NRF) funded by the Korean Government (MEST) (2020R1I1A3069699). The authors have no other relevant affiliations or financial involvement with any organization or entity with a financial interest in or financial conflict with the subject matter or materials discussed in the manuscript apart from those disclosed.

No writing assistance was utilized in the production of this manuscript.

Virtual screening of immunomodulatory medicinal compounds as promising anti-SARS-CoV-2 inhibitors

SARS-CoV-2/COVID-19: viral genomics, epidemiology, vaccines, and therapeutic interventions

Inhibition of translation and immune responses by the virulence factor Nsp1 of SARS-CoV-2

COVID-19: consider cytokine storm syndromes and immunosuppression

Discusses the mechanism of SARS CoV-2 in the development of cytokine storm and immunosuppression

Supporting recovery from COVID-19

Pharmacology IS of NSAIDs in patients with viral infections, including Covid-19: victims or perpetrators?

• Describes the protective effects of nonsteroidal anti-inflammatory drugs in viral infections

Improving the physical realism and structural accuracy of protein models by a two-step atomic-level energy minimization

PRODRG: a tool for high-throughput crystallography of protein-ligand complexes

Proteins Plus: a web portal for structure analysis of macromolecules

PrankWeb: a web server for ligand binding site prediction and visualization

• A bioinformatic database used to predict potential ligand binding sites in a protein structure

PASS: prediction of activity spectra for biologically active substances

• An online webserver used to predict different biological activities of a compound

admetSAR 2.0: web-service for prediction and optimization of chemical ADMET properties

• A bioinformatic database used to predict pharmacokinetic and drug likeness properties of a compound

pkCSM: predicting small-molecule pharmacokinetic and toxicity properties using graph-based signatures

• A bioinformatic database used to predict pharmacokinetic and toxicity properties of a compound

Potential therapeutic effect of turmeric contents against SARS-CoV-2 compared with experimental COVID-19 therapies: in silico study

• An in silico study describing the role of Turmeric Compounds against SARS-CoV-2 targets

Coronavirus (COVID-19) Update: FDA Authorizes Drug Combination for Treatment of COVID-19

Broad-spectrum antiviral activity of naproxen: from influenza A to SARS-CoV-2 coronavirus

Describe the antiviral activity of Naproxen against SARS-CoV-2 and Influenza A

Inhibitory effects of anti-inflammatory drugs on interleukin-6 bioactivity

Controlled drug delivery systems: past forward and future back

Nano based drug delivery systems: recent developments and future prospects

Predicting the need for a reduced drug dose, at first prescription