key: cord-0897124-n6q409km authors: Ismail, Muhammad I.; Ragab, Hanan M.; Bekhit, Adnan A.; Ibrahim, Tamer M. title: Targeting Multiple Conformations of SARS-CoV2 Papain-Like Protease for Drug Repositioning: An In-silico Study date: 2021-02-24 journal: Comput Biol Med DOI: 10.1016/j.compbiomed.2021.104295 sha: aad579c2bc3589c13e2cdb1c890109d516007571 doc_id: 897124 cord_uid: n6q409km Papain-Like Protease (PLpro) is a key protein for SARS-CoV-2 viral replication which is the cause of the emerging COVID-19 pandemic. Targeting PLpro can suppress viral replication and provide treatment options for COVID-19. Due to the dynamic nature of its binding site loop, PLpro multiple conformations were generated through a long-range 1 micro-second molecular dynamics (MD) simulation. Clustering the MD trajectory enabled us to extract representative structures for the conformational space generated. Adding to the MD representative structures, X-ray structures were involved in an ensemble docking approach to screen the FDA approved drugs for drug repositioning endeavor. Guided by our recent benchmarking study of SARS-CoV-2 PLpro, FRED docking software was selected for such virtual screening task. The results highlighted potential consensus binders to many of the MD clusters as well as the newly introduced X-ray structure of PLpro complexed with a small molecule. For instance, three drugs Benserazide, Dobutamine and Masoprocol showed a superior consensus enrichment against the PLpro conformations. Further MD simulations for these drugs complexed with PLpro suggested the superior stability and binding of dobutamine and masoprocol inside the binding site compared to Benserazide. Generally, this approach can facilitate identifying drugs for repositioning via targeting multiple conformations of a crucial target for the rapidly emerging COVID-19 pandemic. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the relevant virus for the Coronavirus disease (COVID- 19) and was originally discovered in Wuhan, China [1] [2] [3] [4] . Based on the situation report of the World Health Organization (WHO) -November 2020, COVID-19 is extremely spreading worldwide over 230 countries and accountable so far for > 96 million cases and > 2 million fatalities (https://covid19.who.int/, accessed on 23 rd Jan 2021). Coronaviruses is part of a large family of enveloped single stranded RNA genome (ssRNA) that belong to the Coronaviridae family. They are categorized into four genera: alpha, beta, gamma, and delta coronaviruses [5] . Some of coronaviruses prompted several respiratory diseases, such as SARS-CoV [6] , middle east respiratory syndrome coronavirus (MERS-CoV) [7] and the pandemic COVID-19 [3] . SARS-CoV-2 are beta coronaviruses [3, 8] with symptoms commonly like other respiratory viruses infection, such as influenza and rhinovirus. [9] Papain-like protease (PLpro) and 3C-like / main protease (3CLpro / Mpro) [10, 11] are critical for the release of 16 non-structural proteins (nsps1-16) after processing the two large polyproteins, pp1a and pp1ab. [12] These polyproteins are created following the virion entry to the host cell where their production is initiated via translation of 5'-terminal open reading frames (ORF1a and ORF1ab). [13] It is known that the establishment of the replicase complex necessary for viral genome replication is conditional on nsps [14] . Nonetheless, PLpro performs a key role for the release of nsp1-3 from the viral polyprotein which is critical for viral replication. Furthermore, PLpro has been reported to negatively regulate the host innate immune response towards the viral infection by its deubiquinating and deISGylating effect. [15, 16] Therefore, PLpro has been identified as an imperative target for viral replication suppression attempts in SARS-CoV and SARS-CoV-2. [15, 17] Structure-based virtual screening (SBVS) remains an essential technique in modern drug discovery. [18] [19] [20] [21] Molecular docking is widely employed in SBVS campaigns and computational drug repositioning for COVID-19, [22] [23] [24] [25] [26] [27] [28] [29] [30] utilizing the structural information of the protein targets to gauge molecular databases and predict the favored binding of compounds. However, the docking tool should be evaluated utilizing benchmarking molecular sets. [31, 32] Also, ensemble docking has emerged as a popular approach which incorporates the protein J o u r n a l P r e -p r o o f flexibility and tackles the problem of false positive results of rigid docking [33] . This can be achieved via Molecular dynamics (MD) simulations that has shown to be a powerful approach for sampling different protein conformations. Usually, screening this ensemble of conformations gives better results than single crystallographic structure [34] . Accordingly, it is the target of this study to investigate samples of the conformational space of SARS-CoV2 PLpro. One reason for this is the limited number of X-ray structures deposited in the Protein Data Bank (PDB) for PLpro (around 24 structures) compared to the main protease (Mpro, around 200 structures). Additionally, compared to Mpro [35] [36] [37] [38] [39] [40] [41] [42] [43] [44] [45] , fewer reports about targeting PLpro is presented in literature. [46] [47] [48] [49] [50] [51] [52] In the context of drug repurposing to treat COVID-19, Balasubramaniam et al. [46] reported potential targeting of three SARS-CoV-2 proteins, RNA-dependent RNA polymerase, papain-like proteinase and helicase, by the antiviral drug elbasvir through virtual screening of 54 FDA-approved antivirals and 3300 investigational drugs. Kandeel et al. proposed some FDA-approved drugs and supplements to target PLpro. Their methodology was based on molecular dynamics (MD) simulations followed by molecular mechanics/generalized Born surface area (MM/GBSA) binding energy calculations. [51] Also, Cavasotto and Di Filippo presented a docking-based screening using a quantum mechanical scoring of a library built from approved drugs and compounds undergoing clinical trials, against three SARS-CoV-2 target proteins, including PLpro. [52] Kouznetsova et al. [48] repurposed FDA-approved drugs along with inhibitors of Hepatitis C Virus (HCV) and Cytomegalovirus (CMV). De Vita et al. [49] adopted an FDA-approved drugs screening approach combining MM-GBSA calculations for ranking drugs along three different SARS-CoV-2 targets, including PLpro, and MD simulations for the best ligand in each target member. Delre et al. [50] suggested known drugs as both covalent and non-covalent inhibitors for PLpro depending on covalent docking, non-covalent docking and MM-GBSA calculations. The mentioned approaches lacked to objectively consider different conformations of the highly flexible loop (T265-H272). This encouraged us to explore the conformational space of PLpro for screening purposes. For this we conducted a long-range molecular dynamics (MD) simulation of PLpro for 1 micro-second. Then the simulation trajectory was clustered to extract representative structures for ensemble docking purposes. Furthermore, guided by our recent J o u r n a l P r e -p r o o f benchmarking study of SARS-CoV2 PLpro [13] , we employed FRED [53, 54] for ensemble screening of DrugBank database against the cluster representatives. The results highlight potential consensus binders to many of the MD clusters as well as the newly introduced X-ray structure complexed with small molecule (PDB ID: 7JRN). Generally, this approach can help in identifying drugs for repositioning via targeting multiple conformations of a crucial target for the rapidly evolving COVID-19 pandemic. The 24 structures of PLpro comprises the following PDB IDs: 6W9C, 6WRH, 6YVA, 6XA9, 6XAA, 6WZU, 6XG3, 7CJD, 7D6H and 7D47 for the apo structures, while 6WUU and 6WX4 for complexes with a peptide-like inhibitors. Also, PDB IDs of 7JRN, 7JN2, 7JIR, 7JIT, 7JIV, 7JIW, 7CJM, 7KOK, 7KOL, 7KOJ, 7D7T and 7CMD have been recently introduced for protein -small molecule complexes (Resolution and R-Value Free are shown in Table T1 in the Based on previous reports and the recently deposited X-ray PDB structures of PLpro, Tyr268 and Gln269 were proposed to be the key residues for small molecule recognition. [13, 17] They encompass the main part of a flexible loop that can adapt various backbone and side chain conformations, as seen in Fig. 1 . However, their unbound X-ray structures appeared to show close conformations (pairwise RMSD < 1A) and their side chains mostly appear to point outwards to the solvent exposed area. [13] Unlike the unbound conformation, the structures complexed with small molecules appeared to show more distinct and ordered conformations for the key Tyr268 and Gln269 where their side chains pointed inwards to offer a hydrophobic wall for interactions with the small molecules. [13] Therefore, based on the data available from the PDB, two main X-ray structures can be extracted representing distinct conformations of the binding site loop (see Hence, we aim in this section to investigate samples of the conformational space of SARS-CoV2 PLpro with a focus on the conformational changes of the key binding site residues. Molecular dynamics simulations were carried out to explore the protein conformational space of the apo/native form (PDB ID: 6W9C) and give insights into the dynamic nature of the binding site. Long-range molecular dynamics simulations for 1 micro-second ware performed to give robust sampling of the protein conformations. The simulation system was subjected to a preliminary minimization step followed by an equilibration step of 6 ns. The system was then subjected to production run of 1 micro-second. Fig. 2 shows Radius of Gyration (Rg) and Root Mean Square Deviation (RMSD) for the protein backbone atoms through the 1 micro-second simulation, where radius of gyration is a measure of the change in the compactness of the protein during the simulation. Its low range of 1Å (48-49Å) indicates the stability of the protein during the simulation and the absence of major conformational change in protein structure. While RMSD is a measure of the deviation of the backbone atoms from the initial protein structure. It is observed that RMSD has a range of 3Å with lowest and highest values of 1Å and 4Å, respectively, which gives an indication of the stability of the protein throughout the simulation. It is also observed that the trajectories are divided into three main ensembles. The first ensemble has RMSD of 1-3 Å (0 -300 ns), the second and third parts are interchangeable from 300-1000 ns showing fluctuation in their RMSD from 2 to 4 Å. Clustering for the simulation trajectory was performed to extract representative structures of the protein different conformations. The J o u r n a l P r e -p r o o f clustering algorithm used calculates the average RMSD between trajectory frames and clusters close frames together. In a first run, we assigned the RMSD measurements to be calculated on the backbone atoms. Fig. 2 shows bar plot of timeline frame distribution of the three resulting clusters with a frame distribution between clusters similar to the distribution of backbone RMSD through the simulation run shown in Fig. 2 , which supports the aforementioned observation. Top: Protein radius of gyration (R g ) through the 1 micro-second simulation time showing low fluctuation of 48-49 Å with low range of 1 Å indicating the protein compactness during the simulation . Middle: protein RMSD of backbone atoms during the simulation time showing three clusters for the protein with the first cluster from 0 to 300 ns, and the second and the third clusters from 300 ns to 1 micro-second. Bottom: Timeline bar plot of protein cluster distribution during the 1 micro-second simulation assigning the RMSD matrix to be measured on the backbone atoms. It shows three clusters corresponding to the three clusters concluded from protein RMSD measurement, with the first cluster in red and second and third clusters in blue and green, respectively. Fig. 1 and reported earlier. [13] The rest of the binding site amino acids (L162-E167, P247-P248, Y264-Y273 and T301) showed lower RMSF values < 1 Å indicating their stability during the simulation. The rest of the protein amino acids displayed low RMSF values highlighting their stability as well (Fig. 3) . Bottom: per residue RMSF with color notation corresponding to the protein (top) color notation. It shows the binding site residues RMSF (red and blue rectangles) and the highly flexible residues (yellow rectangles). The MD simulation resulted in a trajectory file of 14285 protein frames. Considering all the frames in an ensemble docking approach is certainly time-consuming, and redundant outcome can be produced due to high frame similarities in some regions. Therefore, partitioning MD frames into clusters of similar frames addresses the problem of their large number. TTClust program [55] was used for this clustering task, where it calculates the average RMSD between frames and cluster similar frames together. In the first run, the RMSD was assigned to be calculated on the backbone atoms ( Twelve protein structures comprising the ten MD clusters representatives, the starting X-ray structure (PDB ID: 6W9C) and the X-ray structure of PLpro-small molecule complex (PDB ID: 7JRN) were utilized for ensemble docking for screening the DrugBank database [56] . Such set up is to represent various states of the conformation space of PLpro. This approach would enable us to extract a consensus list of drugs that is able to target PLpro in different conformations simultaneously. Based on the results of our previous benchmarking study [13] , we applied FRED docking tool for this ensemble screening approach since it showed promising screening performance in recognizing bioactive molecules in a pool of challenging decoys for both the bound and unbound states of PLpro. [13, 31] Also, we biased our consensus selection for drugs enriched in the newly introduced X-ray structure complexed with a small molecule (PDB ID: 7JRN), because it showed the best benchmarking results among other PLpro X-ray structures. [13] Top enriched 1% of the screened DrugBank database [56] (25 compounds) were retrieved for each docking run and concatenated together to extract commonly enriched drugs in a consensus list. Drugs with 1% enrichment and in consensus rank for 5, 6, 7 and 8 different conformations are listed in Table 1 . As a consensus, ten drugs were shown to commonly target different conformations of PLpro. and -10.38, respectively. Interestingly, their docking score is observed to be the best in the X-ray structure (7JRN) compared to other conformations, therefore, we inspected their postulated binding poses in such X-ray structure, as shown in Fig. 6 and Fig. 7 . Benserazide appeared to be enriched for 5 conformations involving both X-ray structures (PDB IDs: 7JRN and 6W9C). It is an inhibitor of levodopa (L-DOPA) decarboxylation that acts peripherally and does not cross the blood-brain barrier. It is used in combination with the dopamine precursor L-DOPA in Parkinson's disease management to prevent its decarboxylation and facilitate its blood-brain barrier crossing. Upon investigation of Benserazide postulated binding mode and comparing it with the co-crystallized ligand in (PDB ID: 7JRN), two hydrogen bonds with Asp164 side chain and Gln269 backbone atoms are conserved in both poses. These two hydrogen bonds are formed by the hydrazide group of Benserazide and the amide group of the co-crystallized ligand. Two additional hydrogen bonds with Tyr268 and Tyr273 side chain atoms are formed with Benserazide. Additionally, the T-shaped pi-pi stack between the cocrystallized ligand naphthyl ring and Tyr 268, and its pi-alkyl interaction with Pro248, are also conserved with the trihydroxy phenyl ring of Benserazide (Fig. 6A) . can be valuable in prospective drug design projects targeting PLpro (Fig. 6B) . Dobutamine is observed to be enriched for 5 conformations including the X-ray structures (PDB IDs: 7JRN and 6W9C). It is a β 1 agonist used in the treatment of congestive heart failure. As shown in Fig. 7 , its dihydroxy phenyl group maintains a pi-pi stacking with Tyr264 phenyl group, and a hydrogen bond with Tyr273. Additionally, the phenol substituent forms a pi-anion interaction with Asp164 amino acid (Fig. 7A) . Fig. S4 and Fig. S5 . Masoprocol is enriched for 6 conformations including the X-ray structure (PDB ID: 7JRN). It is a lipoxygenase inhibitor used as antineoplastic drug for precancerous skin growths. It possesses a structure similarity with Dobutamine, hence, its docking pose reproduced similar binding pi-anion interactions with Asp164 by both phenyl rings. Two additional hydrogen bonds were also predicted between both Leu162 backbone and Gln269 side chain and the dihydroxy phenyl ring, which also forms a T-shaped pi-pi stack with Tyr268. The other dihydroxyphenyl ring also forms a hydrogen bond with Arg166 side chain. The proposed docking poses of the three drugs (benserazide, dobutamine and masoprocol) were subjected to molecular dynamics simulations for 50 ns to test their stability inside the binding site. The values of R g , RMSD, and RMSF indicate stable simulation for these ligand-PLpro complexed forms (Fig. S6 in the Supplementary Material) . It is noteworthy mentioning that seven other drugs were predicted to target different conformations of PLpro (Table 1) , namely: Amifostine, Carglumic, Pemirolast, Penciclovir, 5-O-phosphono-alpha-D-ribofuranosyl diphosphate, Pyridoxal and Adenosine phosphate. Amifostine is a prodrug which is dephosphorylated by alkaline dephosphatases to yield the free thiol active drug. It is used as cytoprotective in cancer chemotherapy and radiotherapy. Carglumic acid is N-acetylglutamate (NAG) synthetic analog which is used in treatment of hyperammonaemia in NAG-deficient patients through activation of carbamoyl phosphate synthetase 1 (CPS1) liver enzyme to start the urea cycle and convert ammonia into urea. Table T2 in the Supplementary Material. The PLpro apo crystal structure (PDB ID: 6W9C) was subjected to 1 microsecond molecular dynamics simulation at NPT ensemble. The simulation system was prepared using OpenMM Setup utility [57] . AMBER14 [58] force field was used for the protein parameters and TIP3P-FB water model [59] was used for water molecules. A preliminary clean up step was done for the protein PDB file by adding missing OXT atom for LYS315 amino acid. Hydrogen atoms were added at pH value of 7.0. The protein was immersed in a cubic water box with a padding Amber Force Field (GAFF) [60] parameters were used for the three drugs. Both GAFF parameters and ELF10 charges were generated on the fly by the OpenEye toolkits 2019.4.1. [61] Particle Mesh Ewald (PME) method [62] was used for the calculation of long-range electrostatic interactions with a 0.0005 error tolerance for truncating the Ewald summation. A 1.2 nm cut off was used for both PME direct space interactions and Lennard-Jones interactions. The length of all bonds that involve a hydrogen atom, and the water molecules, involving their bond length and angles, were constrained. Langevin integrator was utilized for integrating force equations with step size of 2 ns, temperature of 310 K and friction coefficient of 1 ps -1 . Mote Carlo barostat [63] was used for pressure coupling with assigning 310 K as constant temperature value, 1 atm as constant pressure value and 25 steps as pressure update frequency. The dynamics simulations were carried out using OpenMM python application layer [57] . We started the simulation with a minimization step with 10kJ/mole as energy convergence criteria, then, 3,000,000 equilibration steps (6 ns). Finally, the simulation was carried out for 500 million steps (1 microsecond). A trajectory snapshot was reported every 35,000 steps (70 ps), resulting in a final trajectory file of 14,285 total frames. For the three ligand-PLpro MD runs with benserazide, dobutamine and masoprocol, equilibration was run for 500,000 steps (1 ns), the production was run for 25,000,000 steps (50 ns), and trajectory reporting was set to one snapshot per 1000 steps resulting in a trajectory file of 25,000 snapshots. RMSD, RMSF and radius of gyration was calculated out using ProDy python library [64, 65] . As a first step of trajectory clustering, the trajectory frames were aligned on the protein backbone atoms of the starting structure. Then, pairwise RMSD between frames was calculated and stored in a matrix, for backbone in the first clustering run, and heavy atoms of residues T265-H272 in the second run. This matrix was then used in hierarchical clustering of frames implementing the ward method in calculating the distance between the newly formed clusters. A third clustering run was carried out assigning 10 clusters as output cluster numbers. Clustering was carried out using TTClust program [55] . J o u r n a l P r e -p r o o f DrugBank database [66] was downloaded from www.drugbank.ca (DrugBank -release March 2020) with 2454 ligands. It was prepared by MOE [67] with retrieving one conformer of one protonation state at pH 7.0 implementing Amber: 10EHT forcefield. Multi-conformer file generation for the database was carried out using OMEGA conformer generation module of OpenEye software [68, 69] with a default number of 200 conformers generated for each ligand. The 10 representative protein structures, along with the downloaded protein structures (6W9C and 7JRN), were prepared using MakeReceptor GUI module of OpenEye software [70] . With respect to the 10 representative structures, possible cavities on the protein surface was detected by cavity detection utility of MakeReceptor implementing a molecular probe. Then, the search box dimensions were assigned according to the detected cavity volume of the binding site surrounding the binding site amino acids: L162-E167, P247-P248, Y264-Y273 and T301. Virtual screening was carried out using FRED module of OpenEye software [53, 54] , by docking the multi-conformer DrugBank database file into the prepared protein structures. Protein and protein-ligand 3D interaction figures were generated using UCSF Chimera [71] . Drug repurposing of DrugBank database was carried out by adopting ensemble docking approach for multiple conformations of PLpro SARS-CoV-2 molecular target. Ten PLpro different conformations sampled from 1 µs MD simulation, in addition to two reported PDB PLpro structures, were used for ensemble virtual screening. Based on consensus ranking of the top 1% of the screening hits, three drugs showed to target different conformations with high docking scores and interaction with Tyr268 and Gln269 key binding residues, namely: Benserazide, Dobutamine and Masoprocol. Further MD simulations for these drugs complexed with PLpro indicated higher stability and binding of dobutamine and masoprocol inside the binding site compared to benserazide. Moreover, Metoclopramide, Isocarboxazid, Procainamide and Remoxipride showed the best docking scores over all docking runs combined despite their lower consensus over different conformations. These proposed drugs can be biologically evaluated as potential SARS-CoV-2 treatment candidates. Globally, this protocol can help in J o u r n a l P r e -p r o o f identifying drugs for repositioning via targeting multiple conformations of a crucial target for the rapidly evolving COVID-19 pandemic. The 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. TMI conceived the idea and supervised the project. MII caried out all experiments. MII and TMI wrote the manuscript. All authors have given approval to the final version of the manuscript. All authors agree to be accountable for the content of the work. Additional files providing data about: a) The superposition of three main groups of the PLpro from the PDB (Fig. S1) , a) TTClust results applied to the frames of 1micro-second MD for the backbone heavy atoms for the calculation of average RMSD between frames (Fig. S2) , b) 2D representation for the docking poses of drug-PLpro interactions (Fig. S3 -S5) , c) Radius of gyration, RMSD and RMSF of the protein backbone atoms of the MD runs with benserazide, dobutamine and masoprocol (Fig. S6) , d) Resolution and R-Value Free of the 24 PLpro PDB structures (Table T1) , and e) indications of the proposed drugs (Table T2) . Transmission dynamics and evolutionary history of 2019-nCoV Fighting against the common enemy of COVID-19: a practice of building a community with a shared future for mankind SARS-CoV-2 and coronavirus disease 2019: What we know so far, Pathogens Comparative computational analysis of SARS-CoV-2 nucleocapsid protein epitopes in taxonomically related coronaviruses The structure and functions of coronavirus genomic 3′ and 5′ ends Identification of a novel coronavirus in patients with severe acute respiratory syndrome Isolation of a novel coronavirus from a man with pneumonia in Saudi Arabia Molecular immune pathogenesis and diagnosis of COVID-19 Featuring COVID-19 cases via screening symptomatic patients with epidemiologic link during flu season in a medical center of central Taiwan The Papain-Like Protease of Severe Acute Respiratory Syndrome Coronavirus Has Deubiquitinating Activity Murine Coronavirus Ubiquitin-Like Domain Is Important for Papain-Like Protease Stability and Viral Pathogenesis Corona virus versus existence of human on the earth: A computational and biophysical approach Supporting SARS-CoV-2 Papain-Like Protease Drug Discovery: In silico Methods and Benchmarking Coronaviruses: An Overview of Their Replication and Pathogenesis The SARS-coronavirus papain-like protease: structure, function and inhibition by designed antiviral compounds Deubiquitinating Enzymes in Coronaviruses and Possible Therapeutic Opportunities for COVID-19 Characterization and Noncovalent Inhibition of the Deubiquitinase and deISGylase Activity of SARS-CoV-2 Papain-Like Protease Virtual screening: an endless staircase? Recognizing pitfalls in virtual screening: a critical review Nuclear hormone receptor targeted virtual screening Virtual target screening: validation using kinase inhibitors Virtual screening and repurposing of FDA approved drugs against COVID-19 main protease Anti-HCV, nucleotide inhibitors, repurposing against COVID-19 Molecular Docking and Virtual Screening based prediction of drugs for COVID-19 Déjà vu: Stimulating open drug discovery for SARS Fast Identification of Possible Drug Treatment of Coronavirus Disease-19 (COVID-19) through Computational Drug Repurposing Study Virtual screening, ADME/Tox predictions and the drug repurposing concept for future use of old drugs against the COVID-19 Computational search for potential COVID-19 drugs from FDAapproved drugs and small molecules of natural origin identifies several anti-virals and plant products M-pro) of SARS-CoV-2 through Consensus Docking and Drug Reposition In Silico Structure-Based Repositioning of Approved Drugs for Spike Glycoprotein S2 Domain Fusion Peptide of SARS-CoV-2: Rationale from Molecular Dynamics and Binding Free Energy Calculations, mSystems Evaluation and optimization of virtual screening workflows with DEKOIS 2.0-a public library of challenging docking benchmark sets Applying DEKOIS 2.0 in structure-based virtual screening to probe the impact of preparation procedures and score normalization Potential and Limitations of Ensemble Docking Ensemble Docking in Drug Discovery Quinolines-Based SARS-CoV-2 3CLpro and RdRp Inhibitors and Spike-RBD-ACE2 Inhibitor for Drug-Repurposing Against COVID-19: An in silico Analysis Crystal structure of SARS-CoV-2 main protease provides a basis for design of improved α-ketoamide inhibitors Ligand and structure-based virtual screening applied to the SARS-CoV-2 main protease: an in silico repurposing study Putative Inhibitors of SARS-CoV-2 Main Protease from A Library of Marine Natural Products: A Virtual Screening and Molecular Modeling Study High Throughput Virtual Screening to Discover Inhibitors of the Main Protease of the Coronavirus SARS-CoV-2 Fast Identification of Possible Drug Treatment of Coronavirus Disease-19 (COVID-19) through Computational Drug Repurposing Study Potential SARS-CoV-2 main protease inhibitors Recent Progress in the Drug Development Targeting SARS-CoV-2 Main Protease as Treatment for COVID-19 Overview of antiviral drug candidates targeting coronaviral 3C-like main proteases A comprehensive review on promising anti-viral therapeutic candidates identified against main protease from SARS-CoV-2 through various computational methods Progress in Developing Inhibitors of SARS-CoV-2 3C-Like Protease, Microorganisms Computational target-based drug repurposing of elbasvir, an antiviral drug predicted to bind multiple SARS-CoV-2 proteins Chemical-informatics approach to COVID-19 drug discovery: Monte Carlo based QSAR, virtual screening and molecular docking study of some inhouse molecules as papain-like protease (PLpro) inhibitors Potential COVID-19 papain-like protease PL(pro) inhibitors: repurposing FDA-approved drugs Accelerating the repurposing of FDAapproved drugs against coronavirus disease-19 (COVID-19) Repurposing Known Drugs as Covalent and Non-covalent Inhibitors of the SARS-CoV-2 Papain-Like Protease Repurposing of FDA-approved antivirals, antibiotics, anthelmintics, antioxidants, and cell protectives against SARS-CoV-2 papain-like protease silico Drug Repurposing for COVID-19: Targeting SARS-CoV-2 Proteins through Docking and Consensus Ranking FRED Pose Prediction and Virtual Screening Accuracy FRED and HYBRID docking performance on standardized datasets TTClust: A Versatile Molecular Simulation Trajectory Clustering Program with Graphical Summaries OpenMM 7: Rapid development of high performance algorithms for molecular dynamics Simmerling, ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB Building Force Fields: An Automatic, Systematic, and Reproducible Approach Development and testing of a general amber force field OpenEye Toolkits 2019.4.1 OpenEye Scientific Software A smooth particle mesh Ewald method Isothermal-isobaric molecular dynamics simulations with Monte Carlo volume sampling Evol and ProDy for bridging protein sequence evolution and structural dynamics ProDy: Protein Dynamics Inferred from Theory and Experiments OMEGA 3.1.2.2: OpenEye Scientific Software Conformer Generation with OMEGA: Algorithm and Validation Using High Quality Structures from the Protein Databank and Cambridge Structural Database OpenEye Scientific Software A visualization system for exploratory research and analysis We would like to gratefully acknowledge the support of NVIDIA Corporation with the donation of the Titan Xp GPU used for molecular dynamics calculations. We gratefully acknowledge the support of OpenEye Scientific Software Inc. for offering an academic license.J o u r n a l P r e -p r o o f