key: cord-0923086-sl26ab4h
authors: Federico, Leonardo Bruno; Silva, Guilherme Martins; da Silva Hage-Melim, Lorane Izabel; Gomes, Suzane Quintana; Barcelos, Mariana Pegrucci; Galindo Francischini, Isaque Antônio; Tomich de Paula da Silva, Carlos Henrique
title: Identification of known drugs as potential SARS-CoV-2 Mpro inhibitors using ligand- and structure-based virtual screening
date: 2021-06-25
journal: Future medicinal chemistry
DOI: 10.4155/fmc-2021-0025
sha: 185c223f6c917823e9cc01607c238c87fbca0c06
doc_id: 923086
cord_uid: sl26ab4h

Background: The new coronavirus pandemic has had a significant impact worldwide, and therapeutic treatment for this viral infection is being strongly pursued. Efforts have been undertaken by medicinal chemists to discover molecules or known drugs that may be effective in COVID-19 treatment – in particular, targeting the main protease (Mpro) of the virus. Materials & methods: We have employed an innovative strategy – application of ligand- and structure-based virtual screening – using a special compilation of an approved and diverse set of SARS-CoV-2 crystallographic complexes that was recently published. Results and conclusion: We identified seven drugs with different original indications that might act as potential Mpro inhibitors and may be preferable to other drugs that have been repurposed. These drugs will be experimentally tested to confirm their potential Mpro inhibition and thus their effectiveness against COVID-19.

The new coronavirus pandemic has had a significant impact worldwide, and therapeutic treatment for this viral infection is being strongly pursued. Efforts have been undertaken by medicinal chemists to discover molecules or known drugs that may be effective in COVID-19 treatment -in particular, targeting the main protease (Mpro) of the virus. Materials & methods : We have employed an innovative strategy -application of ligand-and structure-based virtual screening -using a special compilation of an approved and diverse set of SARS-CoV-2 crystallographic complexes that was recently published. Results and conclusion: We identified seven drugs with different original indications that might act as potential Mpro inhibitors and may be preferable to other drugs that have been repurposed. These drugs will be experimentally tested to confirm their potential Mpro inhibition and thus their effectiveness against COVID-19. 19) and in March of the same year a pandemic was declared. Sadly, to date, there is no specific treatment against the disease and research in the field is ongoing [1] [2] [3] . Molecularly, the virus infection mechanism in host cells starts after the virus surface protein, known as Spike protein, binds to the receptor on the surface of human cells, ACE-2. This connection generates a conformational change in the viral structure, allowing the fusion between the virus and the host membrane and, thus, permitting the pathogen to release its genetic material in the infected organism [4] . Two proteases play a fundamental role in the infection process of SARS-CoV-2, with the sequencing of its genome, these proteases were identified as being, respectively, PLpro and 3CLpro -the latter also known as main protease, Mpro [5] .

SARS-CoV 3C-like protease (SARS-CoV 3CLpro), as a part of the replicase polyproteins, cleaves a functional polypeptide and consequently leads to the maturation of SARS-CoV. Because of its functional importance in the SARS-CoV replication cycle, SARS-CoV 3CLpro is considered a potential target to develop novel anti-SARS drugs [6] .

Because COVID-19 has a fast transmission speed, several research groups have been working together and disseminating their results quickly. The rapid dissemination of proteomics results, for example, has allowed researchers in the field of drug development to use these data to propose structures capable of stopping the progression of COVID-19 [7, 8] .

Several structures of the Mpro protease are already available in the protein structure repository Protein Data Bank (PDB) [9] . From these data, it is possible to use virtual screening (VS) techniques, a fast and efficient methodology for development of new drugs that can be linked to a database of drugs that have already been approved by regulatory agencies. In this way, it is possible to shorten the development time for new drugs, allowing these compounds that have known activity against a certain biological target to be studied against another biological target and go through the process known as repositioning of drugs [10] [11] [12] [13] [14] [15] .

In the present work, we sought to identify known drugs (approved by regulatory agencies around the world) with the potential to interact and inhibit the SARS-CoV-2 Mpro (i.e., by applying the so-called drug repositioning or drug repurposing strategy). This could be done by means of employing computational chemistry techniques [15] to identify potential drug candidates for COVID-19 treatment. We used an innovative strategy to scan an inhouse compilation of known approved drugs, specifically by application of hybrid ligand-and structure-based VS approaches. The idea is to speed up the workflow toward in vitro tests, making it possible to foresee data concerning the biological activity of known drugs [5, 11] and, in this way, possibly helping combat the pandemic. Therefore, we determined seven drugs, acting by inhibitory mechanism of enzyme Mpro, as potential candidates for COVID-19 treatment. Furthermore, the binding mode toward Mpro of these seven drugs was compared with known antivirals previously repurposed for COVID-19 (lopinavir, ritonavir and darunavir), which indeed inhibit the enzyme and showed promising in vitro results but have not demonstrated improvement of disease treatment in subsequent clinical studies. Next, these seven drugs will be tested in vitro to confirm their inhibition potential toward Mpro and thus their effectiveness against the new coronavirus.

There is not a single official database of drugs approved by the US FDA; moreover, those available with this classification have differing numbers of compounds. Hence, we compiled a pool of compounds originating from drug databases approved by different agencies or organizations, comprising approximately 8,000 molecules. Our drug library was compiled into three main datasets, Drugbank, the FDA and ZincWorld; these were gathered from the following databases: DrugBank (Release Version 5.1.5) [16] , Zinc15 (subset: world not fda) [17] , Selleckchem subset FDA-approved Drug Library (https://www.selleckchem.com/), Enamine subset approved Drugs Collection (https://www.enaminestore.com/platedsets/fda-approved) e-Drug3D [18] and the BindingDB subset of FDA approved drugs [19] .

Preparation of the drug molecules contained in these databases was carried out by removing duplicate/redundant molecules using OpenBabel [20] , considering InChI (IUPAC InchI identifier), adding hydrogens, setting partial charges and adjusting bond orders with Discovery Studio [21] and saving the corresponding output files in mol2 extension format. Compounds were also processed with FILTER [22] using no properties filtering, but its default functions of canonicalization, pKa normalization and salt were removed. In sequence, drug databases were processed using OMEGA [23, 24] for conformer generation, with the flipper flag turned on (for specification of unassigned stereocenters), adjustment of the strain energy (above the energy of global minimum conformer) be-low 9.0 kcal/mol and a root mean square deviation (RMSD) of 0.6Å as a cutoff for conformer identity. For ROCS [25, 26] , as the ligand-based VS (LBVS) approach, up to 300 conformers generated per molecule; for docking, as structure-based VS (SBVS) approach, only one 3D conformer was generated per molecule.

Selection of crystallographic structures to be used as targets in these structure-based studies, as well as corresponding complex ligands in ligand-based studies, was carried out by inspecting more than 90 available structures in the PDB, considering those of Mpro of the 'old coronavirus' (SARS-CoV) and those available for the new coronavirus (SARS-CoV-2). Initially, we avoided apo structures, MPro structures in complex with ligands bound to different sites or those with covalent bond within the active site of the enzyme. Finally, we selected the following protein-ligand complexes: PDB ID 6W63 ('Structure of COVID-19 main protease bound to potent broad-spectrum non-covalent inhibitor X77,' at 2.10-Å resolution [27] ); PDB ID 5R83 ('PanDDA analysis group deposition -crystal structure of SARS-CoV-2 Mpro in complex with Z44592329,' at 1.58-Å resolution [28] ); and PDB ID 4MDS ('Discovery of N-(benzo [1, 2, 3] triazol-1-yl)-N-(benzyl)acetamido)phenyl) carboxamides as severe acute respiratory syndrome coronavirus SARS-CoV 3CLpro inhibitors -bound to ML300,' at 1.60-Å resolution [29] ).

In addition, we selected one of the most potent noncovalent inhibitors of SARS-CoV MPro, as reported by Ghosh et al. [30] , with IC 50 value of 30 nM and here named C10 (5-chloropyridin-3-yl 1H-indole-4-carboxylate). This ligand has been also used as query/reference in our ligand-based approaches.

MPro protein structures (retrieved from complexes of PDB IDs 6W63, 5R83 and 4MDS) were prepared before docking runs by keeping only the corresponding primary chain; removing water, ligands, metals and cofactors; and adding hydrogens using the Protein Preparation Wizard tool implemented in Maestro software [31, 32] .

It is worth noting that four of the ligands were here used in ligand-based methodologies. As previously cited, three were retrieved from corresponding complexes PDB IDs 6W63, 5R83 and 4MDS, respectively X77 (N-

]acetamide) and the fourth named here C10, as mentioned earlier. Crystallographic poses (conformation and orientation) from the three complexes ligands were kept unchanged. For C10, we performed a docking run using GOLD [33, 34] in which the chosen pose was obtained with best score and in consensus within 10 best scoring poses; thus, this pose was used as query/template/reference in subsequent ligand-based strategies.

Validation of docking procedures on SBVS & queries performance on LBVS SBVS (with docking) strategies have been validated by two procedures: redocking and receiver operating characteristic (ROC) curve analysis. Redocking was performed using GOLD default settings and the previously prepared protein from each of the three selected PDB mentioned earlier. Previous to redocking simulations, simple preparation of compounds X77, Z44592329 and ML300 was carried out by same 'standard' procedure adopted in this work (i.e., that mentioned for preparation of the drug databases).

ROC curve analysis has been applied to the docking strategies used here as a second validation procedure. This docking validation was performed using GOLD software for each one of the three selected and prepared Mpro structures retrieved from PDB. The purpose of this validation procedure is to assess the ability of such docking methodology to recover known active compounds from a dataset of actives contaminated with 'decoys.' One dataset containing approximately 40 actives was used for each of the three inhibitor-enzyme complexes selected in this work. To do this, we randomly selected forty actives for each target from a small library that was built with ∼140 known Mpro inhibitors of SARS-CoV and SARS-CoV-2, compiled from literature [27] [28] [29] [30] [35] [36] [37] [38] [39] [40] [41] [42] [43] [44] , considering either reported data concerning their IC 50 values or existence of respective PDB files (see Supplementary Table 1) . Decoy generation was executed from the webserver Database of Useful Decoys Enhanced (DUD-E) [45] , which furnished 50 decoys per active compound. Each dataset (40 actives and 2000 decoys) was then submitted as input in the GOLD docking software, using the exact same settings as those further used in GOLD runs of VS studies. ROC curves were built using the Screening Explorer webserver (http://stats.drugdesign.f r/) [46] and subsequently evaluated by analysis of the corresponding area under the curve (AUC) values.

Validation of LBVS was also performed by ROC curve analysis, that is, by analyzing the performance of each selected 'query' in the shape similarity comparison screened by ROCS. Selected queries for this purpose were the three ligands X77, Z44592329 and ML300 in their crystallographic poses (conformation), as well as the best Docking + visual inspection scored docking pose of C10, as previously mentioned. For each ligand, we used the same corresponding datasets (as stated in the previous paragraph) compiled from literature, consisting of 40 actives, plus 50 decoys per active compound (generated by DUD-E), and these were used as input database to run vROCS [47] . Results were ranked and evaluated according to the Tanimoto Combo (shape and color) similarity indices, obtained for each one of the corresponding four queries, and thus used to plot ROC curves (using vROCS and the Screening Explorer webserver as well), which were then evaluated according the AUC values obtained.

We performed seven independent VS: three only using structure-based methodology (docking) and four using mixed ligand-and structure-based methodologies (ROCS + docking) (see Table 1 ). Structure-based strategies (i.e., mainly conducted by docking) were done using GOLD software [33, 34] . For this, we used individual drug databases previously prepared as well as the three selected Mpro enzymes. Hence, we selected the best 100 compounds scored by ChemPLP fitness function, that were subsequently inputted into DEREK to assess their toxicity profiles [48, 49] . From these analyses, we filtered compounds that afforded alerts as 'certain' and 'plausible,' according to Lhasa expert system predictions. Survivors compounds were then inspected by docking about probable interactions within active sites of corresponding enzymes so that we only kept molecules that showed consensus poses as well as interactions especially regarding those with His41, and Cys145 [50] and Glu 166 observed in the crystallographic complexes [28, 29, 51] .

Furthermore, we applied mixed VS workflows by means of ligand-based shape and electrostatic potential similarities, using ROCS [25] and EON [52], respectively. For this, we used all the drug databases noted earlier individually, as well as the four ligand queries. We then retrieved the best 100 scored compounds regarding their Tanimoto similarity indices (ROCS TanimotoCombo and EON ET Combo), obtained for each query. These compounds were then assessed with regard their toxicity profiles using DEREK [48, 49] , filtering out compounds that afforded alerts as 'certain' and 'plausible.' Next, survivor compounds were submitted to docking (SBVS) using GOLD to inspect their interactions within the active sites of the corresponding enzymes, especially regarding those with His41 and Cys145, as well as Glu166 observed in the crystallographic complexes.

Finally, the predicted binding modes for the selected compounds were compared with those of the drugs lopinavir, ritonavir, darunavir and boceprevir, which were approved by the FDA for the SARS-CoV-2 treatment.

The decision to use the PDB ID 4MDS complex [29] (of SARS-CoV Mpro), in this work, is due to the existence of few SARS-Cov-2 MPro structures deposited in the PDB in complex with ligands showing high potency or binding affinity. Reliability of using this structure was confirmed by the high sequential similarity and identity shared between the enzymes of the two viruses (sequence identity of 96.1 and 99.02% sequence similarity, with only 12 mutations, all of which are outside the region of interest -the catalytic site) [50, 53] .

The high similarity shared by Mpro of SARS-CoV and Mpro of SARS-CoV-2, together with the fact of that they present practically identical catalytic sites, allows the use of both available co-crystals to elucidate interactions between the enzyme active site and SARS-CoV-2 inhibitors [53] [54] [55] .

Surrounding both Mpro catalytic sites, there is only one difference/mutation, represented by a serine in SARS-CoV-2 and by an alanine in SARS-CoV; nevertheless, their sidechains point to the outer region of the binding site [53] . Moreover, recent studies have shown that both binding sites possess similar substract specificities, suggesting that drugs developed for SARS-CoV might also be efficient in inhibiting the SARS-CoV-2 enzyme [56] .

Furthermore, superposition of the PDB ID 4MDS enzyme structure (MPro SARS-CoV) with the PDB ID 5R83 (Mpro SARS-CoV-2) revealed that the region of the site is composed mostly of the same amino acid residues (and rotamers), with the exception of different rotamers observed for the Asn142 and Gln189 residues of the two structures ( Figure 1) .

The PDB ID 4MDS (of SARS-CoV) enzyme structure has a cavity slightly larger than the corresponding SARS-CoV-2 cavity, where a main chain movement due to the larger and more branched inhibitor is observed compared with that observed in the PDB ID 5R83 structure. The difference in cavity sizes between the MPro structures of the two viruses is, in fact, only due to the size of the corresponding complexed inhibitors. This hypothesis is confirmed by the PDB ID 6W63 (SARS-CoV-2 Mpro) structure, which presents an inhibitor equally large and branched compared with the PDB ID 4MDS inhibitor structure of SARS-CoV Mpro.

Another important aspect that motivated us to use the PDB ID 4MDS structure was the presence of an inhibitor chemical and that it is structurally different from the inhibitor present in the PDB ID 6W63 structure, even though the cavities are very similar and correspondent in size and shape. The virus-1 enzyme inhibitor is more branched than the SARS-CoV-2 enzyme and similar to some of the main protease inhibitors on the market. Furthermore, our comparative analyzes here performed with both the Mpro structures of SARS-CoV and SARS-CoV-2 contribute to the clarification of their structural similarity, corroborating with studies that were previously carried out with this objective [24, 50, 51] and thus allowing the use of previous results from the literature (since 2003) on relevant structural characteristics of compounds to act as Mpro SARS-CoV-2 inhibitors as well. Validation of the docking & shape-queries performances Regarding to the docking validation procedures, a standard redocking protocol was carried out. Redocking was performed using the GOLD software to verify the ability of the docking methodology used here to reproduce the crystallographic poses (in terms of conformation and orientation -or 'pose') of the original ligands located at the three proteins binding sites. Results are shown in Figure 2A -C, where a suitable convergence of the crystallographic poses for X77, Z44592329, and ML300 (respectively of the PDB IDs 6W63, 5R83 and 4MDS) can be seen. In fact, the top-ranked poses of each ligand, obtained using GOLD, have RMSD values of 0.6, 0.9 and 1.7Å, respectively, in agreement with that expected for a redocking validation procedure [57, 58] .

A second docking validation procedure was also employed, as described earlier, to assess the ability of this docking methodology and recover known active compounds from a dataset containing a large number of DUD-E-generated decoys.

It is worth noting that we randomly compiled one individual dataset of 40 actives for each validation with each PDB structure. In fact, the three crystallographic Mpro inhibitors (two of SARS-CoV-2 and one of SARS-CoV) have different molecular weights, volumes and shapes, inducing and requiring an apparent flexibility of the active site of such enzyme to occupy and interact with residues of interest. Moreover, the small library compiled from the literature (see Supplementary Table 1 ), from where we retrieved 40 actives to rationalize the validation datasets, includes inhibitors of both the SARS-CoV-2 and the SARS-CoV Mpro structures and reasonable structure diversity. Thus, due to the structural similarity of the enzymes and inhibitors of the two SARS-CoV Mpro crystallographic structures used here, the random selection of actives in the validation datasets proved to be efficient.

For the GOLD docking results obtained for each enzyme and each selected dataset, compounds were ranked using the ChemPLP score function and the corresponding binaries codes (1 to active, 0 to decoy) used to plot ROC curves, correlating sensitivity (y-axis) against 1-specificity (x-axis). In general, sensitivity means the 'true positive,' whereas 1-specificity signifies a 'false positive' [59] . In this way, compounds were classified regarding true-false positive-negative ratios and represented in terms of sensitivity and specificity to build the ROC curves. Figure 3 demonstrates the resulting curves and respective AUC values, which were 0.645, 0.732 and 0.655 for the PDB IDs 6W63, 5R83 and 4MDS, respectively.

We also validated our ROCS queries before applying them to large databases in VS prospections by means of analyses of ROC curves generated by vROCS. When using the ligand X77 we obtained an AUC value of 0.801, for Z44592329 an AUC value of 0.839, for ML300 an AUC value of 0.895, and for C10 an AUC value of 0.824, as depicted in Figure 3 .

Docking AUC results that were not close to the unit could be attributed to trivial differences observed between the active compounds and the decoys, to active compounds belonging to limited chemical classes or to their possessing similar physicochemical properties. However, because these limitations should rely principally on the selection of the compounds dataset, which was indeed constituted of variable chemical classes (see Supplementary  Table 1 ), we believe it represents a satisfactory result because it is relatively far from a 0.5 AUC random value.

Finally, both the structure-based (by docking, with GOLD) and ligand-based (by shape similarity, with ROCS) methodologies employed were successful in each validation procedure according to their ROC curve analyses. They showed sufficiently good performance to ensure reliable use of the corresponding protocols and queries in prospective VS studies.

Considering the previously described methodological strategies (see Materials & methods), including database preparation and validation of methodologies, we carried out the VS workflows. The databases were used separately so that results were obtained for each database and for each of three proteins (three SBVS workflows) and four ligands (four LBVS + SBVS workflows), for a total of seven VS workflow results for each database. Thus, we could analyze consensus between results in each VS workflow as a way to select final potential hits. Supplementary Figures 1-3) .

By comparing results obtained with workflows consisting of only SBVS and of both LBVS + SBVS, the LBVS + SBVS strategies were shown to be more consistent because results obtained using LBVS workflows showed drugs as potential hits with a higher number of consensus docking poses (higher homogeneity between them) and also the indication that these compacts can interact with residues Glu166 and His41 in their corresponding final visual inspection. Use of these two methodologies together has shown significant improvement in results obtained in research from our group [60, 61] , allowing elimination of extensive or bulky ligands in the initial stages of the screening protocols so that only molecules with similar characteristics to the queries of interest proceed to later stages of each workflow.

Importantly, the docking algorithm used to perform conformer generation was the same as that used in all the docking simulations, and 10 docking runs (poses) were conducted for each ligand. Hence, drugs with potential to be repurposed for Mpro inhibition that presented a homogeneity among the corresponding sets of docked poses were selected (some ligands presented unanimity for 10 poses in this analysis; for examples, see metyrapone and zonisamide in Table 2 ).

Moreover, these seven selected drugs maintained key interactions within the active site of enzymes in a similar manner as observed in the crystallographic complex structures, particularly with regard to the hydrogen bond with Glu166. Furthermore, in the majority of poses, the key interaction with His41 was established, which was also observed in the crystallographic complex PDB ID 4MDS; when this is not present, as in the case of the drug felodipine, it was substituted in one interaction with Gly143. Table 2 shows the potential drugs to be repurposed as Mpro inhibitors, retrieved from corresponding databases and obtained from respective VS workflows. These same compounds are also illustrated in Figure 4 , showing their 3D and 2D representations and their potential as Mpro inhibitor candidates.

The first three drugs shown in Table 2 -metyrapone, rufinamide and zonisamide -had the most consistent results. Metyrapone, which has been used for the diagnosis of adrenal insufficiency and treatment of Cushing syndrome, was recovered in three simulations, each with a different query, demonstrating the strong tendency of this drug to interact effectively with the Mpro enzyme, even with small conformational changes in it. The anticonvulsant rufinamide was also recovered in three simulations. The adjuvant zonisamide used in the treatment of epilepsy also showed interesting results, achieving consensus in all the docking poses in two simulation routes in two drug databases.

The other four drugs selected, despite having been recovered by only one of the simulations performed, presented a high number of docking poses in consensus. Lacosamide, which is used as an adjuvant in epilepsy treatment, presented a consensus of nine poses in the simulation using the FDA database and the reference structure PDB ID 6W63. Apatinib and Pramiracetam, used to treat gastric cancer and as a CNS stimulant, respectively, demosntrated consensus of eight poses among their corresponding results. Finally, the calcium channel blocker Felodipine, originally to control hypertension, remains a potential drug to treat of COVID-19; this drug showed a consensus of seven poses using the FDA database and the structure of the PDB ID 4MDS inhibitor as a reference.

Results of docking simulations allowed us to infer that proximity between Glu166 and Hys41 residues at the catalytic site of Mpro favored the selection of compounds bearing groups that are able to act as hydrogen bond acceptors next to other groups capable of transferring electrons by pi-pi or cation-pi interactions. In this way, interactions presented by the majority of drugs in this study show oxygen with electron pairs near an aromatic ring.

In addition, we evaluated the likelihood that these seven drug hits indeed interact with the binding pocket from each SARS-CoV-2 Mpro used in docking studies. For this, we employed the software DoGSiteScorer [62, 63] , which is a tool for pocket detection and protein druggability assessment. Preliminarily, we detected the most significant pockets in each Mpro protein (PDB IDs 4MDS, 5R83 and 6W63) and found that the respective active sites (where corresponding original crystallographic ligands were bound) presented favorable values of clusters' number, chemical descriptors and druggability scores (see Supplementary Table 2) .

Furthermore, as shown in the last column of Table 2 , when retrieving the best docking poses for each of the seven drug hits and evaluating which pocket should preferably host these ligands, all of them presented unique and highest scores toward potential binding in the active sites (see also Supplementary Figures 1-3) . Worth noting that the main active site (with highest number of clusters) predicted using DoGSiteScorer are represented by P1 for PDB ID 4MDS, P2 for 5R83, and P2 + P6 for 6W63. Thus, this evaluation should strongly indicate that top poses prefer the active site over other sites present in the surface of SARS-CoV-2 Mpro.

The presence of a hydrogen bond acceptor group is represented by a carbonyl in metyrapone ( Figure 4A ); an amide oxygen in rufinamide, lacosamide and pramiracetam ( Figure 4B , D & F, respectively); a sulfonamide oxygen in zonisamide ( Figure 4C ); and an ester in felodipine ( Figure 4G ). In the case of apatinibe (Figure 4e) , the drug presents a secondary amine in place of carbonyl that is able to interact with amide oxygen of Glu166 by means of donating a hydrogen bond. All these groups present distances ranging from 2.0 to 2.5Å from the amide nitrogen of Glu166 -or from the corresponding oxygen for apatinibe -which therefore should have strong hydrogen bonding (represented by yellow dotted lines in 3D images and by purple arrows in 2D images).

Moreover, selected drugs present an aromatic ring or a protonated amine next to the aforementioned oxygen that forms a hydrogen bond with Glu166. Both groups are capable of interacting by electron transfer with imidazole from Hys41 (represented by blue dotted lines in 3D and green lines in 2D images), with the exception of felodipine, which does not interact with this residue. There is a pyrimidine ring in metyrapone, a triazole in rufinamide, benzioxazole in zonisamide and benzene in lacosamide and apatinibe; all of these are distanced at a maximum of 4.4Å when face-to-face or 5.5Å when edge-to-face with Hys41 imidazole. These groups and distances allow the occurrence of pi-pi interactions with residue when drugs bear aromatic rings or cation-pi interactions when they bear a protonated amine.

Other interactions are yet to be observed, such as lacosamide, which other than previously mentioned interactions, is still able to establish a hydrogen bond between its ether acceptor group and the side chain nitrogen from Gly143 and also donate a hydrogen bond from its amide nitrogen to carbonyl of Asn142.

To enhance the quality of selection of these seven compounds, we compared their binding interactions with the corresponding ones from known drugs that have been reported as promising repurposing candidates for COVID-19. We used lopinavir, ritonavir and darunavir, three antivirals known to act as protease inhibitors that prevent viral replication of HIV and also have been shown to inhibit SARS-CoV-2 Mpro in vitro [64] . Nonetheless, the use of these drugs in clinical studies to treat patients with COVID-19 have not improved the clinical picture of hospitalized patients or patients with severe cases of the disease [65] [66] [67] [68] [69] .

In addition, these three drugs did not show a high affinity for Mpro inhibition when evaluating the corresponding docking results, obtained by applying the same protocol and conditions performed in the present work ( Figure 5 ).

Lopinavir shows in vitro data concerning the inhibition of viral replication, with an EC 50 value of 5.73 μM, CC 50 of 74.44 μM and SI of 12.99 [70] . However, in our docking studies, it had no consensus results within 10 poses generated using GOLD ( Figure 5A ), and none of the poses had simultaneous key interactions with more than one of the main residues (Glu166, His41 and Gly143) present at the enzyme catalytic site. The same behavior was observed for ritonavir ( Figure 5B ), which is usually used in clinical practice with lopinavir-ritonavir and had an EC 50 value of 8.63 μM, CC 50 of 74.11 μM and SI of 8.59 [70] .

Darunavir, which had EC 50 value of 46.41 μM, CC 50 >81 μM and SI >1.75 [70] , presented consensus in three of 10 poses generated by docking simulation, with variation of the positions only the sulfonamide and isobutyl groups (oriented toward the entrance of the active site), but it did not show key interactions with the main residues of the catalytic site ( Figure 5C ). This drug also did not have satisfactory clinical data and has not been approved as a possible COVID-19 treatment [68, 69] .

These drugs have presented significant and promising results concerning the inhibition of COVID-19 viral replication; however, this was not noted for in vivo or clinical studies. This could be due to the low bioavailability of free drug (not bound to plasma proteins) -that is, in their corresponding pharmacological active forms [69] . Hence, we expect that the seven drugs repurposed here will have an even greater affinity toward inhibition of SARS-CoV-2 Mpro and overcome the drawbacks related to clinical and in vivo studies, given the results we hope to obtain with appropriate assays that will be performed in the future.

In view of the urgent need for a safe and reliable COVID-19 treatment, drug repositioning studies prove to be a quick and important alternative for discovery of drugs that are already approved or regulated and could also act against the new coronavirus. The use of this approach eliminates a long and important step during drug developmentpharmacokinetics evaluation -which is responsible for a large number of failures in drug development.

Thus, we have developed and described two systematic VS workflows that present excellent validation results, with the aim of repositioning drugs to act as potential inhibitors of the SARS-CoV-2 Mpro. More specifically, seven workflows were conducted using LBVS approaches, and we ultimately selected seven promising drugs with strong indications of acting against the coronavirus.

Following the same docking protocols adopted in the selection of the seven drugs with repurposing potential, we compared the methodology application of known repurposed drugs -lopinavir, ritonavir and darunavir. These three drugs are described in the literature as potential and expressive inhibitors of the enzyme under study; however ,they did not provide clinical improvement in patients hospitalized with COVID-19, which may be due to their low bioavailability in the human body. Furthermore, analysis of the main docking interactions for the seven compounds that we selected have shown more consensual results, pointing to a greater affinity for Mpro than lopinavir, ritonavir and darunavir, which have not shown an equal consensus of poses in docking simulations. These drugs will be followed by in vitro inhibition tests with the SARS-CoV-2 Mpro and further in vitro assays to reduce viral load in cells infected with the new coronavirus and prove their potential.

The worldwide impact of the COVID-19 pandemic has pushed the development of novel strategies to contain this viral spread. In addition to known vaccination plans that countries have been adopting, it would be worthwhile to determine efficient therapeutic strategies to circumvent viral mortality in serious cases. Thus, drug repurposing is a strategy that could be successfully applied in our work, along with the use of hybrid virtual screening. Hence, we identified 7 known drugs that provide interesting insights for SARS-CoV-2 Mpro inhibition and binding interactions, that shall, after corresponding assays are developed, help this field of drug discovery to evolve.

• Hybrid ligand-and structure-based virtual screening • Drug repurposing using a compilation of approximately 8,000 approved drugs • A variety of recently published SARS-CoV-2 crystallographic complexes • Robust validation of both 3D shape-similarity and docking methodologies • Consensus analysis of docking poses from ligand based virtual screening • Seven known drugs studied as potential Mpro inhibitors • Comparison with known drugs repurposed for SARS-CoV-2, such as lopinavir and ritonavir • Further experimental assays will confirm whether these seven drugs have the potential to treat COVID-19 future science group 10 .4155/fmc-2021-0025

To view the supplementary data that accompany this paper please visit the journal website at: www.futurescience.com/doi/suppl/10.4155/fmc-2021-0025

Financial & competing interests disclosure

The authors have no 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. This includes employment, consultancies, honoraria, stock ownership or options, expert testimony, grants or patents received or pending, or royalties.

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

Epidemiologic features and clinical course of patients infected with SARS-CoV-2 in Singapore

An interactive web-based dashboard to track COVID-19 in real time

Discovery of small molecule PLpro inhibitor against COVID-19 using structure-based virtual screening, molecular dynamics simulation, and molecular mechanics/Generalized Born surface area (MM/GBSA) calculation

Immune responses in COVID-19 and potential vaccines: Lessons learned from SARS and MERS epidemic

Potential inhibitors against papain-like protease of novel coronavirus (SARS-CoV-2) from FDA approved drugs

Discovery of a novel family of SARS-CoV protease inhibitors by virtual screening and 3D-QSAR studies

Structural genomics of SARS-CoV-2 indicates evolutionary conserved functional regions of viral proteins

Structural basis of SARS-CoV-2 3CLpro and anti-COVID-19 drug discovery from medicinal plants

The Protein Data Bank

Repositioning FDA drugs as potential cruzain inhibitors from Trypanosoma cruzi: virtual screening, in vitro and in vivo studies

Virtual screening in drug discovery -a computational perspective

Network-based drug repurposing for novel coronavirus 2019-nCoV/SARS-CoV-2

Repurposing of known anti-virals as potential inhibitors for SARS-CoV-2 main protease using molecular docking analysis

QSAR Modeling of SARS-CoV M pro inhibitors identifies sufugolix, cenicriviroc, proglumetacin, and other drugs as candidates for repurposing against SARS-CoV-2

• Good example of Computer-Aided Drug Design (CADD) and Quantitative Structure-Activity Relationship (QSAR) studies performed in the same context

Analysis of therapeutic targets for SARS-CoV-2 and discovery of potential drugs by computational methods

• Good example of structure-based virtual screening in same context

DrugBank: a comprehensive resource for in silico drug discovery and exploration

ZINC 15 -ligand discovery for everyone

Data sets representative of the structures and experimental properties of FDA-approved drugs

BindingDB: a web-accessible database of experimentally determined protein-ligand binding affinities

Open Babel: an open chemical toolbox

Systèmes -Discovery Studio Visualizer

Conformer generation with OMEGA: algorithm and validation using high quality structures from the Protein Databank and Cambridge Structural Database

Comparison of shape-matching and docking as virtual screening tools

Structure of COVID-19 main protease bound to potent broad-spectrum non-covalent inhibitor X77 -6W63

PDB ID:5R83 PanDDA analysis group deposition -crystal structure of SARS-CoV-2 main protease in complex with Z44592329

Discovery of N-(benzo[1,2,3]triazol-1-yl)-N-(benzyl)acetamido)phenyl) carboxamides as severe acute respiratory syndrome coronavirus (SARS-CoV) 3CLpro inhibitors: identification of ML300 and noncovalent nanomolar inhibitors with an induced-fit binding

Design, synthesis and antiviral efficacy of a series of potent chloropyridyl ester-derived SARS-CoV 3CLpro inhibitors

Protein and ligand preparation: parameters, protocols, and influence on virtual screening enrichments

Improved protein-ligand docking using GOLD

Development and validation of a genetic algorithm for flexible docking

Discovery of potent anilide inhibitors against the severe acute respiratory syndrome 3CL protease

Structure of M pro fromSARS-CoV-2 and discovery of its inhibitors

α-Ketoamides as broad-spectrum inhibitors of coronavirus and enterovirus replication: structure-based design, synthesis, and activity assessment

Identification of antiviral drug candidates against SARS-CoV-2 from FDA-approved drugs

• From a previous screen of 3000 approved drugs against SARS-CoV, 50 were selected to perform antiviral activity assays with SARS-CoV-2, finding potential in 24

Crystal structure of SARS-CoV-2 main protease provides a basis for design of improved α-ketoamide inhibitors

Structure-guided design of potent and permeable inhibitors of MERS coronavirus 3CL protease that utilize a piperidine moiety as a novel design element

High-throughput screening identifies inhibitors of the SARS coronavirus main proteinase

Discovery, synthesis, and structure-based optimization of a series of N-(tert-Butyl)-2-(N-arylamido)-2-(pyridin-3-yl) acetamides (ML188) as potent noncovalent small molecule inhibitors of the severe acute respiratory syndrome coronavirus (SARS-CoV) 3CL

Synthesis and evaluation of pyrazolone compounds as SARS-coronavirus 3C-like protease inhibitors

Design, synthesis, and evaluation of inhibitors for severe acute respiratory syndrome 3C-like protease based on phthalhydrazide ketones or heteroaromatic esters

Directory of Useful Decoys, Enhanced (DUD-E): better ligands and decoys for better benchmarking

Screening explorer -an interactive tool for the analysis of screening results

In silico tools for sharing data and knowledge on toxicity and metabolism: Derek for Windows, Meteor, and Vitic

Distinguishing between expert and statistical systems for application under ICH M7

Rapid identification of potential inhibitors of SARS-CoV-2 main protease by deep docking of 1.3 billion compounds

6W63 structure of COVID-19 main protease bound to potent broad-spectrum non-covalent inhibitor X77

Conserved interactions required for inhibition of the main protease of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)

• Describes the similarity between the SARS-CoV and SARS-CoV-2 structures

A search for medications to treat COVID-19 via in silico molecular docking models of the SARS-CoV-2 spike glycoprotein and 3CL protease

The continuing 2019-nCoV epidemic threat of novel coronaviruses to global health -the latest 2019 novel coronavirus outbreak in Wuhan

Substrate specificity profiling of SARS-CoV-2 main protease enables design of activity-based probes for patient-sample imaging

Comparing protein-ligand docking programs is difficult

Validation of molecular docking programs for virtual screening against dihydropteroate synthase

Virtual screening workflow development guided by the "receiver operating characteristic" curve approach

Discovery of novel dual acetylcholinesterase inhibitors with antifibrillogenic activity related to Alzheimer's disease

Identification of novel αβ-tubulin modulators with antiproliferative activity directed to cancer therapy using ligand and structure-based virtual screening

Combining global and local measures for structure-based druggability predictions

Analyzing the topology of active sites: on the prediction of pockets and subpockets

Therapeutic drugs for SARS-CoV-2 treatment: current state and perspective

No benefit for lopinavir-ritonavir in severe COVID-19

A trial of lopinavir-ritonavir in adults hospitalized with severe Covid-19

Early administration of lopinavir/ritonavir plus hydroxychloroquine does not alter the clinical course of SARS-CoV-2 infection: a retrospective cohort study

Lack of antiviral activity of darunavir against SARS-CoV-2

SARS-CoV-2 and HIV protease inhibitors: why lopinavir/ritonavir will not work for COVID-19 infection

Nelfinavir inhibits replication of severe acute respiratory syndrome coronavirus 2 in vitro