key: cord-0864698-nsaz7ksw authors: Baby, Krishnaprasad; Maity, Swastika; Mehta, Chetan H.; Suresh, Akhil; Nayak, Usha Y.; Nayak, Yogendra title: SARS-CoV-2 entry inhibitors by dual targeting TMPRSS2 and ACE2: An in silico drug repurposing study date: 2021-04-05 journal: Eur J Pharmacol DOI: 10.1016/j.ejphar.2021.173922 sha: 11c4c013851b7f221ef28e7960b12cebc3912414 doc_id: 864698 cord_uid: nsaz7ksw The coronavirus disease (COVID-19) is spreading between human populations mainly through nasal droplets. Currently, the vaccines have great hope, but it takes years for testing its efficacy in human. As there is no specific drug treatment available for COVID-19 pandemic, we explored in silico repurposing of drugs with dual inhibition properties by targeting transmembrane serine protease 2 (TMPRSS2) and human angiotensin-converting enzyme 2 (ACE2) from FDA-approved drugs. The TMPRSS2 and ACE2 dual inhibitors in COVID-19 would be a novel antiviral class of drugs called “entry inhibitors.” For this purpose, approximately 2800 US-FDA approved drugs were docked using a virtual docking tool with the targets TMPRSS2 and ACE2. The best-fit drugs were selected as per docking scores and visual outcomes. Later on, drugs were selected on the basis of molecular dynamics simulations. The drugs alvimopan, arbekacin, dequalinum, fleroxacin, lopinavir, and valrubicin were shortlisted by visual analysis and molecular dynamics simulations. Among these, lopinavir and valrubicin were found to be superior in terms of dual inhibition. Thus, lopinavir and valrubicin have the potential of dual-target inhibition whereby preventing SARS-CoV-2 entry to the host. For repurposing of these drugs, further screening in vitro and in vivo would help in exploring clinically. Coronaviruses, once regarded as relatively innocuous respiratory pathogens in human beings. The previously reported two instances of respiratory tract infections, SARS, and MERS are both zoonotic diseases. SARS-coronavirus (SARS-CoV) was initially found in civets in markets and majorly transmitted through bats. MERS coronavirus typically spreads from camels to humans. SARS-CoV-2 erupted in Wuhan, China, and became a worldwide pandemic COVID-19 (Ye et al., 2020) . The frequent indications of COVID-19 are varying degrees of fever, dry cough, headache, anosmia and dysgeusia, sore throat, and diarrhoea. Usually, death is due to comorbidity such as diabetes, hypertension, renal and liver disorders, asthma etc. (Lovato and de Filippis, 2020) . The inflammatory response in alveoli of the lung makes it difficult for gas exchange. It is now evident that the death in COVID-19 is because of the acute release of inflammatory mediators and "cytokine storm" (Coperchini et al., 2020) . Hence, immunomodulators such as dexamethasone and prednisolone were approved by regulatory agencies across the world. The SARS-CoV-2 generally spreads by nasal droplets. Hence social distancing, isolation, and quarantine are the best-suggested methods for preventing the disease . Currently, prevention of infection, and supportive treatments such as oxygen therapy and mechanical ventilation, are the strategies for managing COVID-19 till the specific treatment is established (Sanders et al., 2020) . Vaccines are currently under trial, and the efficacy of convalescent plasma is limited. The approved specific antiviral agents such as remdesivir and favipiravir do not produce a complete cure . Hence, the fast-tracking of drug discovery and development is inevitable. Many scientists have proposed drug repurposing. In drug-repurposing, the in silico screening offers an excellent and rapid drug screening (Deshpande et al., 2020) . We have used computational tools to study the drug-targeting "main protease" (M pro ) and "RNA-dependent RNA polymerase" (RdRp) for the repurposing of US-FDA-approved drugs (Baby et al., 2020a (Baby et al., , 2020b (Baby et al., , 2021 . The crucial step in computational drug repurposing is the selection of validated drug-target. We envisaged that inhibiting the host cell entry would prevent viral load. The access of coronavirus into the host cell is enabled by two main drug targets, "human angiotensin-converting enzyme 2" (ACE2) and transmembrane serine protease 2 (TMPRSS2) (Hoffmann et al., 2020) . Hence we selected these two targets for our studies. The SARS-CoV-2 spike protein (S) attaches to ACE2 with 10-20 folds higher affinity than the SARS-CoV spike glycoprotein. It might contribute to the spread of COVID-19 from humans to humans at high-speed (Wrapp et al., 2020) . The TMPRSS2, a host-cell serine-protease enzyme, and ACE2 facilitate endocytosis and endosome-containing coronavirus (Hou et al., 2020) . The dual-targeting, as a viral entry inhibitor, is a novel approach for anti-SARS-CoV-2 treatment. Hence, in this study, we examined FDA approved drugs for inhibiting serine protease TMPRSS2 and human ACE2. We prepared the protein backbone of the TMPRSS2 enzyme by homology modelling and used ACE2 for docking and molecular dynamics simulations. The in silico modelling, screening, and simulation experiments were executed on Ubuntu desktop workstation, on the Maestro graphical user interface of Schrödinger Suite (www.schrodinger.com) in an Intel® Xenon® Gold 6130 CPU @ 2.10 GHz x 64 processors, Quadro P620/ PCle/SSE2 graphics card, and 8 GB RAM. There was no available and curated crystal structure for TMPRSS2 in the PDB database. Hence, the amino acid sequence of the catalytic domain of TMPRSS2 protein was imported from www.uniprot.org with the code of UniProtKB-O15393 (Boutet et al., 2016) . BLAST Homology search tool was employed to look for an acceptable prototype to create a homology TMPRSS2 model. The homology model was created using the human plasma kallikrein (PDB ID: 5TJX) model as a reference utilising Schrodinger Prime Module (Jacobson et al., 2004; Vyas et al., 2012) . The modelled structure was later modified by Refine Loops software. Thus the optimised homology model was validated by the "Ramachandran plot" (Ramachandran et al., 1963) . Based on the resolution, the best TMPRSS2 homology model was later used in molecular docking experiments. The protein TMPRSS2 was created by homology-modelling and was further optimised (Madhavi Sastry et al., 2013) . A proper ionisation status was defined at pH 7.4 through PROPKA3 Rostkowski et al., 2011) . The protein structure was then minimised using the restricted minimisation technique. The X-ray diffraction crystallographic human ACE2 associated inhibitor carboxypeptidase protein structure (PDB ID: 1R4L) was acquired from PDB (www.rcsb.com) with a resolution of 3.00 Å (Madhavi Sastry et al., 2013) . The Protein Prep Wizard (PPW) was used to refine and optimise the protein structure. The Prime tool in Schrodinger was applied for adding the missing hydrogen bonds, side chains, and missing residues. The bulky water molecules were removed Rostkowski et al., 2011) . The active site and essential catalytic residues were preserved in the protein structure. The chemical structure of compounds licensed by USFDA (~2800 Nos) was saved from www.drugbank.ca. Proper ligand alignments were produced on the molecules, wherein they formed 3D coordinates. (Chen and Foloppe, 2010) . The accurate degree of ionisation was assessed at pH 8.0. Through LigPrep, 2-dimensional chemical structures were converted to 3-dimensional structures with multiple configurations of each input chemical structure with diverse ionisation states, tautomers, stereochemistries, ring conformations. The resulting structures were then saved in Maestro format (Shelley et al., 2007) . Finally, the moiety was converted into a lower energy level by lowering energy using the OPLS3e force field (Roos et al., 2019) . The Ligand structures so prepared by LigPrep were then used for docking studies. The docking was executed utilising the Glide Module in Schrodinger suit. The linking position of molecular ligands was analysed based on the bound ligand and its centroid identified by the Glide grid tool . The co-crystalline ligand was separated during ligand docking. Initially, ~2800 FDA-acquiesced molecules were docked independently into the receptor grids developed on TMPRSS2 and ACE2 using Schrodinger's High Throughput Virtual Screening (HTVS) (Damm-Ganamet et al., 2019) . Later, the Standard Precision (SP) docking technique was used on top molecules, high docking scores. Also, their protein-ligand interaction contacts were analysed independently against both targets. Afterwards, based on the SP docking score, ligand-protein contacts and visual analysis of specific molecules were examined, and the molecules were chosen for further analyses. This study was performed using the Schrodinger Desmond platform, and it helps assess patterns of interaction among proteins and ligands and time-dependent protein-ligand interaction associations and conformational dynamics of complex systems. The predicted structures of human ACE2 and TMPRSS2 combined with ligands observed during Glide SP docking were used in MD simulations as explained in our previous work (Baby et al., 2020b) . The TMPRSS2 homology model was found suitable for the docking studies when it was validated by the Ramachandran plot ( Fig. 1) . Further, human ACE2 protein preparation was found to be suitable for docking studies. Original ligand docking of ~2800 FDA approved drugs were performed separately to recognise how they bind to host proteins TMPRSS2 and ACE2. The ligand docking offered ample knowledge of the affinity and orientation of ligand-protein complexes which impairs the protein activity. Initially, ~2800 molecules were subjected to HTVS separately on both the targets. Then top drugs obtained on HTVS (top-500) based on the docking scores were further docked under Standard Precision (SP) docking. All the top-500 molecules were observed by visual examination for intermolecular interactions. The specific drugs, those interacted strongly with the TMPRSS2 homology model, and ACE2 were identified and shortlisted. Tables 1 and 2 represents intramolecular interactions in SP-docking with TMPRSS2 and ACE2, respectively. The frequent hits against target TMPRSS2 and ACE2 were identified upon visual examination. The drugs such as valrubicin, lopinavir, fleroxacin, alvimopan, arbekacin, and dequalinium exhibited a compelling binding mode against both the targets. Fig. 2 represents the structures of the shortlisted molecules. These molecules interacted with the various residues of protein through H-bond, Salt-bridge, π-π interaction, and π-cation interaction with TMPRSS2 and ACE2. Figs. 3 and 4 depict valrubicin and lopinavir bound against TMPRSS2 and ACE2 as superior to other ligands. The molecular dynamics (MD) simulation investigation established structural insights into the protein-ligand complex's molecular interaction based on observed variations in root mean square deviation (RMSD) and their stability under simulated physiological conditions using the Maestro Desmond module. MD simulation was performed on the ligands that demonstrated stronger ligand-protein interactions, namely valrubicin, lopinavir, fleroxacin, alvimopan, arbekacin, and dequalinium following SP-docking. Examination of RMSD variations demonstrated showing catalytic sites, because the TMPRSS2 crystal structure was not available in PDB, the sequence was obtained from UniProtKB-O15393 and the homology model was created using the human plasma kallikrein PDB ID-5TJX; B) Ramachandran plot for optimised homology model of TMPRSS2, here the number of outliers in homology generated molecule of TMPRSS2 is less, which proves that the developed homology model is fit for further studies; C) ACE2 with catalytic site (PDB ID: 1R4L); th figure are obtained from the Maestro graphical user interface of Schrödinger Suite. The TMPRSS2 bound with valrubicin complex displayed hydrophobic and hydrophilic interactions during MD simulations. The narrative of the protein structure backbone specified during the MD simulations was compared with the original structure for the RMSD analysis. After initial variations due to equilibration for 50 nsec, the RMSD for protein structures was maintained from 0.40 to 4.2 Å until the simulation. (Fig. 5) . It was observed that the fluctuations for protein have remained within 0.20 Å, which indicated a highly stable protein structure. Similarly, the fluctuations remained within 1.5 Å, which indicated a highly stable ligand in the protein complex. Fig. 7 illustrates the ligand-protein interactions recorded during the MD simulation. The amino acid residue TRP461 made maximum contact with the ligand throughout simulations. In the simulation, bridged hydrogen bonding, π-cation interactions, and direct hydrogen bonding with protein's residues were observed. The residue TRP461 interacted with the ligand via a π-π stacking, which accounted for 30 and 40% of the interactions. The valrubicin-bound ACE2 illustrated hydrophilic interactions during the MD simulations. A plot of protein backbone structure described during the MD simulations were compared to the original structures for the RMSD. After initial variations due to balancing, for a period of 50 nsec, until the end of the simulation, the RMSD for protein structures was held between 0.30 and 1.8 Å (Fig. 6) . The protein fluctuations that persisted within 0.30 Å proved a stable protein structure. Fluctuations for the ligand was persisted within 2.1 Å, suggesting a stable ligand in the protein complex. Fig. 8 displays the ligand-protein contacts identified during the MD simulation. The residue ASN149 made full interaction with valrubicin during the simulations. Bridged hydrogen bonding with different amino acid residues, namely GLU145, ASN149, ASP367, THR371, were observed during the simulation. TMPRSS2 bound to the lopinavir complex demonstrated mainly hydrophilic associations during molecular simulations. The plot of the protein structure of the backbone listed during the MD simulations matched with the initial structure of the RMSD. After preliminary variations due to equilibration, for initial 50 nsec, the RMSD for protein structure stayed between 1.4 and 2.8 Å till the completion of simulations (Fig. 5) . The fluctuations that persisted within 0.8 Å showed a stable protein structure. The variations persisted within 7.6 Å, suggesting a relatively unstable ligand within the protein complex. Fig. 9 depicts the ligand-protein contacts and the protein-ligand contacts recorded during the MD simulations. The TRP461 and LYS300 residues made full interaction with the lopinavir. A direct hydrogen bonding was demonstrated with amino acid residue TRP461 and the unsaturated double-bonded ligand oxygen, which accounted for 70% of the interactions. Further, the bridged hydrogen bonding with HIS296, LYS300, GLY439, and GLY462 formed by 30, 77, 34, and 39% of the interactions were observed. Similarly, a direct hydrogen bond was formed, which forms 70% of the interactions. Many of the associations found amongst protein and ligand throughout docking was reliably maintained during MD simulation. The protein ACE2 linked to the lopinavir complex demonstrated mostly hydrophilic interactions during MD simulation. The plot of the target protein backbone described during MD simulations were compared with the original structure for the RMSD analysis. After preliminary variations due to equilibration, for a duration of 50 nsec, the RMSD between 0.25 and 1.75 Å until the end of the simulation Fig. 6 . The fluctuations for the protein that persisted within 0.25 Å suggested a highly stable protein structure. The fluctuations remained within 1.8 Å, indicating a highly unstable ligand in the protein complex. Fig. 10 depicts the ligand-protein contacts and the protein-ligand contacts recorded during the MD simulation. The residue ARG273 made maximum contact with the ligand throughout simulations. Bridged H-bonding, π-cation interactions, and direct hydrogen bonding with protein residues were observed in the simulation. The residue ARG273 interacted with the ligand via a π-cation and two direct hydrogen bonds, which accounted for 52, 68, and 79%, respectively. Similarly, bridged hydrogen bonding was with ASP367 and GLU406, which accounted for 62 and 46%, respectively. These results suggest that drugs like valrubicin and lopinavir have superior binding among shortlisted drugs. Further, the other drugs such as fleroxacin, arbekacin, alvimopan, and dequalinium were also identified as useful dual inhibitors by SP-docking, binding orientation, and interaction patterns. These ligands, whose molecular interaction pattern is represented in the supplementary files (Figs. S1-S8). Evolutionarily, SARS-CoV-2 comprises four major structural proteins, namely spike (S), envelope (E), membrane (M), and nucleocapsid A: The protein-ligand interaction diagram. B: Amino acid residues that interact with the ligand in each trajectory frame. The residues making strong contacts are shown in darker colour shade. C: Ligand interaction with the amino acid residues of protein during MD simulation. Interactions that occur more than 30% of the simulation time are shown. Fig. 9 . Interaction diagram of lopinavir with TMPRSS2 during MD simulations. A: The protein-ligand interaction diagram. B: Amino acid residues that interact with the ligand in each trajectory frame. The residues making strong contacts are shown in darker colour shade. C: Ligand interaction with the amino acid residues of protein during MD simulation. Interactions that occur more than 30% of the simulation time are shown. (N) protein, together with many accessory proteins (Jiang et al., 2020) . The S-protein is a transmembrane protein in the outer part of the virus, with a molecular weight of around 150 kDa (Walls et al., 2020) . Infection with SARS-CoV-2 relies on human ACE2 for entry into the cell and TMPRSS2 (Hoffmann et al., 2020) . SARS-CoV-2 primarily infects the pneumocytes and macrophages of the lungs. The S-protein is cleaved and activated by TMPRSS2, before membrane fusion. This makes the transfer of viral contents into the cytoplasm of the pneumocytes (Hoffmann et al., 2020) . ACE2 has been shown to transform angiotensin-2 into angiotensin 1-7, which attaches to the Mas-receptor to counteract the angiotensin 2 receptor leading to vasoconstriction and pro-inflammatory effects (Simões E Silva et al., 2013) . Angiotensin 2 moderates its effects by activating MAPK at p38 (Park et al., 2007) . Angiotensin 1-7 Mas-receptor stimulation reduces p38 MAPK activation to ameliorate inflammation (Yu et al., 2018) . Therefore, blocking ACE2 after coronavirus entry might cause Angiotensin 2 induced p38 activation to dominate in heart and lungs because of Angiotensin 1-7 downregulation. This in turn, stimulates the inflammatory process and creates a positive feedback loop due to p38 mediated upregulation of ADAM17, a protease which severs the ACE2 ectodomain to further the defensive function of local ACE2 (Zipeto et al., 2020) . In the present computational investigation, we observed six molecules, namely alvimopan, arbekacin, dequalinium, fleroxacin, lopinavir, and valrubicin, as a potent inhibitor of both TMPRSS2 and ACE2 by docking and MD-simulations. These shortlisted drugs will be helpful in mitigating the inflammatory processes associated with TMPRSS2 and ACE2. The usage of these established drugs to SARS-CoV-2 might render them superior over other currently used drugs. The MD simulations indicated the association of valrubicin and lopinavir against TMPRSS2 and human ACE2 demonstrated stronger binding stability. The RMSDs of valrubicin and lopinavir complexes with both the proteins TMPRSS2 and ACE2 were stable than other shortlisted drugs. Moreover, both the ligands exhibited an excellent binding interaction with amino acid residues of TMPRSS2 throughout the simulations. The same was observed with valrubicin and lopinavir interactions with the catalytic triad of human ACE2. Hence, we have ranked valrubicin, and lopinavir as best hits against the targets TMPRSS2 and ACE2. Alvimopan is a peripheral opioid mu-receptor antagonist for improving postoperative gastrointestinal recovery in patients with ileus and is devoid of severe adverse effects (Xu et al., 2016) . The drug arbekacin is efficacious towards gram-positive bacteria, and methicillin-resistant Staphylococcus aureus, without any chronic, dangerous, or life-threatening adverse events. It is efficacious when administered as inhalation in patients with pneumonia caused by P. aeruginosa and MRSA (Lee and Lee, 2016) . Dequalinium chloride is an antimicrobial agent generally used for skin and vagina (Mendling et al., 2016) . Advantages of dequalinium is that it also has modulations on inflammation, especially pharyngitis (Cooper et al., 2008) . Fleroxacin is a tri-fluorinated quinolone with a broad-spectrum activity, often used in managing intra-abdominal abscesses and of traveller's diarrhoea as well as in pharmacokinetic studies. Fluoroquinolone displayed antimicrobial activity against several infections caused by mycobacteria, mycoplasma, chlamydia, legionellae, etc. (Blouin et al., 1992; Steffen et al., 1993) . We observed a good in silico binding of valrubicin and lopinavir among shortlisted drugs could be beneficial in preventing the entry of SARS-CoV-2 and in reducing molecular virulence. Lopinavir was discovered to be a type-1 aspartate protease antagonist of the HIV, which has shown in vitro inhibition effect against SARS-CoV-2 (Choy et al., 2020; Kang et al., 2020) . The combination of lopinavir-ritonavir was tried for the management of COVID-19 in India, and the reports are not discouraging (clinical trial registry CTRI/2020/06/026196 (http://ctri.nic.in/) (Bhatnagar et al., 2020) . Likewise, ritonavir co-administration impair lopinavir metabolism and improve ritonavir plasma levels. The same study reports the possibilities of combining lopinavir-ritonavir with ribavirin in COVID-19 (Zeitlinger et al., 2020) . The triple antiviral therapy, lopinavir-ritonavir, and ribavirin with interferon beta-1b was found safe and superior over the lopinavir-ritonavir multicentric, prospective, open-labelled, randomised, phase-2 trial in adults with COVID-19 (Hung et al., 2020) . However, Cao et al., 2020 did not see much advantage from treatment with lopinavir/ritonavir relative to normal care in the hospitalised adult with severe Covid-19 symptoms (Cao et al., 2020) . Though A: The protein-ligand interaction diagram. B: Amino acid residues that interact with the ligand in each trajectory frame. The residues making strong contacts are shown in darker colour shade. C: Ligand interaction with the amino acid residues of protein during MD simulation. Interactions that occur more than 30% of the simulation time are shown. lopinavir-ritonavir combination was not successful with standard medical treatment in a randomised, controlled, open-label study in patients with SARS-CoV-2 infection and positive symptoms, Cao et al., 2020 have not entirely rejected the use in the treatment of COVID-19. The subgroup analysis suggests that the lopinavir-ritonavir combination still has efficacy (Cao et al., 2020) . Further, the emergency use of lopinavir-ritonavir combination is accepted in the USA, Japan, Italy, China, and India (Frediansyah et al., 2020) . Additionally, valrubicin influences many intracellular physiological functions affecting nucleic acid metabolism, which was earlier found to be a semi-synthetic, anthracycline doxorubicin counterpart could be an alternative drug based on its binding stabilities on dual targets such as ACE2 and TMPRSS2. Wang J., 2020 has reported the bindings to major protease (M pro ) by computational studies for in silico repurposing in COVID-19 and shortlisted compounds, including lopinavir (Wang, 2020) . Our current study found it also has actions as dual inhibition as an entry inhibitor. Thus valrubicin has added advantages if it is tried further for in vitro and in vivo studies despite its limitations as cardiotoxicity at higher doses (Lotrionte et al., 2013) . Hence lopinavir and valrubicin, which might act as potential dual inhibitors of TMPRSS2 and ACE2, could be further studied in vitro and in vivo models of TMPRSS2 and ACE2 inhibition to confirm the in silico predictions. In this study, we provided possible SARS-CoV-2 entry inhibitors through inhibiting TMPRSS2 and ACE2. We conducted virtual screening of commercially available drugs licensed by the FDA, and analysis through MD simulations shortlisted six drugs. Valrubicin and lopinavir were found superior among the drugs shortlisted. Since these drugs have the least degree of toxicity, they may be explicitly tested for their beneficial role in COVID-19 infection followed in vitro, in vivo, and further clinical investigations. This study has not received any special funding, but authors would like to acknowledge other important funding supports during the study. Thanks to Manipal Academy of Higher Education, Manipal, Karnataka India, for supporting with "TMA Pai fellowship" to Krishnaprasad B and Akhil Suresh. The authors are also thankful to ICMR, New Delhi, for Senior Research Fellow ICMR-SRF to Swastika Maity (45/33/2019/ PHA/BMS). The authors would like to thank the Department of Pharmaceutics, Manipal College of Pharmaceutical Sciences Manipal Academy of Higher Education, Manipal, Karnataka India, for Schrodinger Systems, and Computer simulations which was procured by Usha YN under a grant from DST-SERB, New Delhi, India (EMR/2016/007006). No ethical approval was necessary for these experiments. There are figures supplied as supportive or supplementary files. Other than this, the manuscript has all the data generated by the experiments. But, if any data other than these, authors are ready to supply. Please see the MD simulation data of alvimopan, arbekacin, dequalinium, and fleroxacin as supplementary figures (Figs. S1-S8). "Authors declare that there is no interest in publishing the data generated by this research work". In silico drug repurposing of penicillins to target main protease Mpro of SARS-CoV-2 Targeting SARS-CoV-2 RNA-dependent RNA polymerase: An in silico drug repurposing for COVID-19, 9 Targeting SARS-CoV-2 main protease: a computational drug repurposing study Lopinavir/ritonavir combination therapy amongst symptomatic coronavirus disease 2019 patients in India: protocol for restricted public health emergency use Fleroxacin pharmacokinetics in patients with liver cirrhosis Uniprotkb/swiss-prot, the manually annotated section of the uniprot knowledgebase: how to use the entry view A trial of lopinavir-ritonavir in adults hospitalized with severe covid-19 Drug-like bioactive structures and conformational coverage with the ligprep/confgen suite: comparison to programs MOE and catalyst Remdesivir, lopinavir, emetine, and homoharringtonine inhibit SARS-CoV-2 replication in vitro Physical distancing, face masks, and eye protection to prevent person-to-person transmission of SARS-CoV-2 and COVID-19: a systematic review and meta-analysis Diagnosis and treatment of infectious pharyngeal inflammation The cytokine storm in COVID-19: an overview of the involvement of the chemokine/chemokine-receptor system Accelerating lead identification by high Throughput virtual screening: prospective case studies from the pharmaceutical industry In silico molecular docking analysis for repurposing therapeutics against multiple proteins from SARS-CoV-2 Antivirals for COVID-19: a critical review Glide: a New approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy SARS-CoV-2 cell entry depends on ACE2 and TMPRSS2 and is blocked by a clinically proven protease inhibitor New insights into genetic susceptibility of COVID-19: an ACE2 and TMPRSS2 polymorphism analysis Triple combination of interferon beta-1b, lopinavir-ritonavir A hierarchical approach to all-atom protein loop prediction Neutralizing antibodies against SARS-CoV-2 and other human coronaviruses In vitro activity of lopinavir/ritonavir and hydroxychloroquine against severe acute respiratory syndrome coronavirus 2 at concentrations achievable by usual doses Clinical usefulness of arbekacin Review and meta-analysis of incidence and clinical predictors of anthracycline cardiotoxicity Clinical presentation of COVID-19: a systematic review focusing on upper airway symptoms Potential therapeutic agents against COVID-19: what we know so far Protein and ligand preparation: parameters, protocols, and influence on virtual screening enrichments Use of locally delivered dequalinium chloride in the treatment of vaginal infections: a review PROPKA3: consistent treatment of internal and surface residues in empirical p K a predictions p38 mitogen-activated protein kinase inhibition ameliorates angiotensin II-induced target organ damage Stereochemistry of polypeptide chain configurations OPLS3e: extending Force Field Coverage for Drug-Like Small Molecules Graphical analysis of pH-dependent properties of proteins predicted using PROPKA Pharmacologic treatments for coronavirus disease 2019 (COVID-19): a review Epik: a software program for pKa prediction and protonation state generation for drug-like molecules ACE2, angiotensin-(1-7) and Mas receptor axis in inflammation and fibrosis Treatment of travellers' diarrhoea with fleroxacin: a case study Homology modeling a fast tool for drug discovery: current perspectives Structure, function, and antigenicity of the SARS-CoV-2 spike glycoprotein Fast identification of possible drug treatment of coronavirus disease-19 (COVID-19) through computational drug repurposing study Cryo-EM structure of the 2019-nCoV spike in the prefusion conformation Alvimopan combined with enhanced recovery strategy for managing postoperative ileus after open abdominal surgery: a systematic review and meta-analysis Zoonotic origins of human coronaviruses Angiotensin-converting enzyme 2-angiotensin (1-7)-Mas axis prevents pancreatic acinar cell inflammatory response via inhibition of the p38 mitogen-activated protein kinase/nuclear factor-?B pathway PK/PD of Anti-Infectives Study Group (EPASG) of the European Society of Clinical Microbiology, Infectious Diseases (ESCMID), 2020. Pharmacokinetics/pharmacodynamics of antiviral agents used to treat SARS-CoV-2 and their potential interaction with drugs and other supportive measures: a comprehensive review by the PK/PD of antiinfectives study group of the European society of antimicrobial agents ACE2/ADAM17/ TMPRSS2 interplay may Be the main risk factor for COVID-19 Supplementary data to this article can be found online at https://doi. org/10.1016/j.ejphar.2021.173922.