key: cord-0693400-ism9brmx
authors: Verma, Shalja; Pandey, Anand Kumar
title: Factual insights of the allosteric inhibition mechanism of SARS-CoV-2 main protease by quercetin: an in silico analysis
date: 2021-01-12
journal: 3 Biotech
DOI: 10.1007/s13205-020-02630-6
sha: caceec76a99e120bad327a8a9424d2869560697c
doc_id: 693400
cord_uid: ism9brmx

SARS-CoV-2 main protease (M(pro)) cleaves the viral polypeptide 1a and 1ab in a site-specific ((L-Q|(S, A, G)) manner and produce functional enzymes for mediating viral replication. Numerous studies have reported synthetic competitive inhibitors against this target enzyme but increase in substrate concentration often reduces the effectiveness of such inhibitors. Allosteric inhibition by natural compound can provide safe and effective treatment by alleviating this limitation. Present study deals with in silico allosteric inhibition analysis of quercetin, against SARS-CoV-2-M(pro). Molecular docking of quercetin with M(pro) revealed consistent binding of quercetin at a site other than active site in multiple runs, with the highest binding energy of − 8.31 kcal/mol, forming 6 H-bonds with residues Gln127, Cys128, Lys137, Asp289 and Glu290. Molecular dynamic simulation of 50 ns revealed high stability of M(pro)-quercetin complex with RMSD values ranging from 0.1 to 0.25 nm. Moreover, native-M(pro) and M(pro)-quercetin complex conformations extracted at different time points from simulation trajectories were subjected to active site-specific docking with modelled substrate peptide (AVLQSGFR) by ZDOCK server. Results displayed site-specific cleavage of peptide when docked with native-M(pro). While substrate peptide remained intact when docked with M(pro)-quercetin complex, also the binding energy of peptide reduced from 785 to 86 from 1 to 50 ns as quercetin induced alterations in the active site cavity reducing its affinity for the substrate. Further, no interactions were noticed between peptide and active site residues of M(pro)-quercetin complex conformations at 40 and 50 ns. Hence, quercetin displayed effective allosteric inhibition potential against SARS-CoV-2 M(pro), and can be developed into an efficient treatment for COVID-19.

From December 2019, when SARS-CoV-2 first introduced its devastating effect to human race in China, to till date pandemic situation of COVID-19, around 1.6 million deaths and more than 72 million confirmed cases have been reported by World Health Organisation (WHO) at a global level pandemic World Health Organisation accessed on 18 December 2020). Though the current stats and highly contagious nature of the disease convey urgent need for an incredibly potent treatment, vaccines or synthetic antiviral agent developed till date have been reported to have many side effects.

Severe acute respiratory syndrome corona virus-2 (SARS-CoV-2), the novel RNA virus and the etiological agent of COVID-19 possesses about 82% identity with the sequence of SARS-CoV (Severe acute respiratory syndrome coronavirus), which victimised thousands of lives in 2003, and 50% sequence identity with Middle-east respiratory syndrome corona virus (MERS-CoV) (Ma et al. 2020 ). This single-stranded, positive (+) sense, enveloped RNA virus belongs to genus beta-coronavirus and family coronaviridae. The genome of SARS-CoV-2 is the largest among RNA viruses and comprises around 26-32 kilo bases (Chen et al. 2020) . The large genome codes for many non-structural and Shalja Verma and Anand Kumar Pandey have contributed equally.

Page 2 of 10 structural proteins but the mutation prone nature of viral genome makes it difficult to decide the target for the development of therapies. Till date, around 116 mutations have been reported in SARS-CoV-2 (Khailany et al. 2020) . The structural proteins such as the spike, nucleocapsid, membrane and glycoproteins are known to be more susceptible to mutations as compared to non-structural proteins (Wu and Yan 2005; McBride et al. 2014; Ortega et al. 2020 ). Thus, non-structural protein-like main protease or 3-chymotrypsin-like protease, RNA-dependent RNA polymerase, helicase and accessary proteins can be of immense importance for therapy development. Targeting conserved regions in these non-structural counterparts can in turn increase the effectiveness of the treatment (Wu et al. 2020) . Main protease (M pro ) of SARS-CoV-2, a 33.8 kDa enzyme, plays a vital role in the cleavage of viral polyproteins (pp1a and pp1ab) and results in the release of functional replicase enzyme essential for transcription and replication of virus (Ma et al. 2020; Chen et al. 2020; Zhang et al. 2020) . By autolytic cleavage, M pro first cleaves itself from the polyproteins and subsequently by catalytic cleavage at about 11 highly conserved sites, it cleaves off the other essential proteins. M pro has an exclusive substrate preference for amino acid Q (glutamine) at site P1 (L-Q|(S, A, G)) of polyproteins, and such substrate specificity has not yet been reported for any human protease. (Ma et al. 2020; Chen et al. 2020; Zhang et al. 2020; Zhu et al. 2011) . Thus, the extremely essential contribution of M pro in viral replication along with lack of homologous counterparts in the human host, makes it evident that M pro can be a highly effective target for the development of enzyme inhibitor based therapeutics.

Several researchers have investigated and developed competitive inhibitors against SARS-CoV-2 M pro which mechanistically inhibit the enzyme, when present in a considerable amount, by competing with the substrate polypeptide to occupy the enzyme's active site. Comparatively high concentration or exposure of the substrate to enzyme may limit the efficiency of competitive inhibitors (Ma et al. 2020; Zhang et al. 2020; Zhu et al. 2011; Jin et al. 2020a, b; Mengist et al. 2020) . In contrary, allosteric inhibitors do not compete with the substrate; thus, the high concentration of substrate does not affect their efficiency. They bind at a different site other than the active site and pose inhibition by passing perturbations at the active site which result in conformational changes in active site. These modulations reduce the affinity of the active site for the substrate (Liu et al. 2008; Strelow et al. 2012; Wenthur et al. 2014 ). Recent researches have reported very high sensitivity (~ 1900%) of SARS-CoV-2 M pro in transmittance of structural alterations across the enzyme via long-range connections, compared to SARS-CoV M pro . Such a high tendency of the SARS-CoV-2 M pro enzyme to convey any perturbation occurring at one site to other seems to be promising and encourage investigation of potential allosteric inhibitors (Estradaa 2020; Doshi et al. 2016; Negre et al. 2018) .

Quercetin, a natural flavonol, has been reported for its antiviral effect against a wide variety of influenza virus strains and provide inhibition at an early infection stage (Wu et al. 2015) . It also possesses antioxidant, anti-inflammatory, anti-cancer, anti-bacterial and anti-proliferative effects and thus contributes to the treatment of numerous diseases (Anand David et al. 2016; Batiha et al. 2020) . Fruits (such as apple, raspberries, cherries, and red grapes), leafy green vegetables and some herbs are major natural sources of quercetin and can be easily exploited for its extraction. Being a natural compound, it shows least toxicity in vivo compared to the till date known synthetic antiviral compounds (Panche et al. 2016; Salehi et al. 2020 ).

In the present study, we analysed the effective allosteric inhibition of SARS-CoV-2 M pro by a natural well-established nutraceutical compound, quercetin. Molecular docking and molecular dynamic (MD) simulations were performed to evaluate the structural binding affinity of quercetin for M pro . Allosteric inhibitory effect of quercetin was revealed by analysing the binding affinity of M pro and M pro -quercetin complexes (conformations extracted at different time points from MD simulation trajectory), for the short essential peptide sequence (AVLQSGFR containing the cleavage site of M pro ) of polypeptide 1a and 1ab at the enzyme's active site.

The 3D coordinate structure file of SARS-CoV-2 M pro (PDB ID: 6Y84) was procured form RCSB Protein Data Bank. The 3D coordinates of quercetin were generated using online smile translator tool (https ://cactu s.nci.nih.gov/trans late/). AutoDockTools-1.5.6 (http://autod ock.scrip ps.edu/ resou rces/adt) developed by "The Scripps Research Institute" was employed for preparing M pro protein structure file and performing blind docking analysis. The water molecules were deleted from the protein structure, all atoms were assigned to AD4 type, Kollman charge of − 135.203 was added, polar hydrogens were added and non-polar hydrogen were merged. Molecular docking of M pro with quercetin was performed using Lamarckian Genetic Algorithm parameters: number of GA runs 10, population size 150, max number of evaluations 2,500,000, max number of generations 27,000, gene mutation rate 0.02 and crossover rate of 0.8. The lowest binding energy conformation was considered best for further analysis. The docking complex conformation and 2D interactions were studied using system-based freeware Pymol 

The best conformation of M pro -quercetin docked complex having the lowest binding energy was considered for MD simulation study using GROMACS 2020 (http://www.groma cs.org/) (Pronk et al. 2013; Abraham et al. 2015) . MD simulations of M pro -quercetin complex, Mpro and quercetin were performed individually for 50 ns each at a temperature of 300 K. The topology of quercetin was prepared using online server CgenFF (https ://cgenff .umary land.edu/) and the topology of protein was made using GROMACS (Vanommeslaeghe et al. 2010 ). The CHARMM36 force field was used for all three simulations. Solvation was performed in a rhombic dodecahedron box using TIP3P three-point water model followed by ions addition to neutralise the system. Energy minimization was performed using steepest descent minimization algorithm for 50,000 steps. NVT and NPT equilibration simulations were done for 100 ps with position restrain at 300 K. Finally, equilibrated system was subjected to MD production run of 50000 ps at pressure and temperature of 1 bar and 300 K, respectively. During MD production run, trajectory coordinates were recorded at every 0.02 ps. The root mean square deviation (RMSD), root mean square fluctuation (RMSF), radius of gyration (Rg), and number of hydrogen bonds (h-bond) formed with respect to time were evaluated and plots were prepared using Origin software [version 2020 (https ://www.origi nlab.com)]. The MD simulation trajectories were analysed by Visual Molecular Dynamic (VMD) software (https ://www.ks.uiuc.edu/Resea rch/vmd/) developed by University of Illinois at Urbana-Champaign (Humphrey et al. 1996) . The 3D coordinates at every 100 ps were extracted from the trajectory files into a single PDB file for further evaluation by Pymol freeware (https ://pymol .org/2/). The 3D images of M pro -quercetin complex were generated at different time points using Pymol software. Eventually, individual 3D coordinates depicting M pro -quercetin complex conformations and M pro conformations at different time points (1 ns, 10 ns, 20 ns, 30 ns, 40 ns and 50 ns) during MD run were extracted from the full 3D coordinates PDB files made from simulation trajectories for further docking analysis to study allosteric inhibition effect of quercetin.

The substrate peptide (AVLQSGFR), which is a part of substrate polypeptide 1a and 1ab and contains the cleavage site of SARS-CoV-2 Mpro, was modelled using PEPFOLD3 online tool (https ://biose rv.rpbs.univ-paris -dider ot.fr/servi ces/PEP-FOLD3 /) that works on Hidden Markov Model and was energy minimised using GROMACS minimizer (https ://mobyl e.rpbs.univ-paris -dider ot.fr/cgi-bin/porta l.py#forms ::Groma cs_Minim izer) associated with the server (Lamiable et al. 2016; Lindhal et al. 2001) .

To study the allosteric inhibition effect of quercetin against SARS-CoV-2 M pro , the binding affinities of SARS-CoV-2 M pro conformations for its substrate peptide at its active site were analysed when quercetin was bound to it at an independent site (in M pro ) other than the active site. The binding affinity of M pro conformations in absence of quercetin were also analysed for reference. The 3D coordinate structures of SARS-CoV-2 M pro and SARS-CoV-2 M pro -quercetin complex conformations were extracted form MD simulation trajectories at 1 ns, 10 ns, 20 ns, 30 ns, 40 ns and 50 ns and the reference 3D structure of SARS-CoV-2 M pro (PDB ID: 6Y84) was procured from RCSB Protein Data Bank. These structures were docked with the modelled peptide (AVLQSGFR) structure using ZDOCK online server (Pierce et al. 2011 (Pierce et al. , 2014 Vreven et al. 2011 ). The active site residues (His41, Cys145, Asp187, Arg188, Gln189 and Thr190) were selected to assure specific binding of the substrate peptide at the enzyme's active site. The ZDOCK score of the best docking cluster was considered in each docking performed for the evaluation of binding affinity of M pro (either in complexed conformation with quercetin or in isolation) for its substrate peptide. The binding interactions were analysed using Pymol freeware (https ://pymol .org/2/). The ZDOCK scores and the docking interactions of M pro and M pro -quercetin complex structures for the substrate peptide were recorded and interaction images were generated using Pymol software. The changes in structure, interactions and ZDOCK scores were compared between M pro -peptide complex and M pro -quercetin-peptide complex conformations.

The present study deals with the in-depth analysis of allosteric inhibition effect of quercetin against the SARS-CoV-2 M pro enzyme for the development of a natural, non-toxic, high potential inhibitor. Molecular docking of SARS-CoV-2 M pro with quercetin resulted in binding of quercetin at an independent site other than the active site in majority of docking runs proving high binding affinity of quercetin for the site. The best docking conformation resulted in binding energy of − 8.31 kcal/mol and Ki value of 810.43 nM. The 3D binding cavity of quercetin in M pro and the 2D interactions revealed effective interaction of quercetin resulting in 67 Page 4 of 10 formation of 6 hydrogen bonds with the conserved residues Gln127, Cys128, Lys137, Asp289 and Glu290 of M pro protein (Fig. 1) . Sencanski et al., in their recent in silico study on SARS-CoV-2 M pro inhibitors predicted allosteric site comprising residues Arg4, Lys5, Arg131, Lys137, Trp207, Asp289, Leu287, Glu288, Gln127, Asp197, Gly138, Glu290, Tyr126, Phe291, and Leu286 which is same as depicted for M pro -quercetin complex in Fig. 1 . They reported Raltegravir, Rolitetracycline, Tolvaptan, Ciclesonide and Rescinnamine as high potential allosteric inhibitors with binding energies ranging from − 7.9 to − 7.2 kcal/mol which proves high binding potential of quercetin to the allosteric site of M pro (Sencanski et al. 2020 ). El-Baba et al., in their mass spectrometry based study for the investigation of allosteric inhibitor of M pro revealed binding of potential allosteric inhibitor x1187 at the binding site comprising residues Phe8, Met6, Gln127, Ser139, Asp295, Arg298, and Gln299. Gln127 being common to the mentioned studies and to our analysis may be considered as effective residue for allosteric inhibition thus depicting the accuracy of allosteric site selection in the present study (El-Baba et al. 2020) . Iftikhar et al. in their in-silico study on inhibitors against three potential enzyme targets of SARS-CoV-2 reported Bagrosin as an allosteric inhibitor, binding with an energy of − 7.9 kcal/ mol at an allosteric site lying close to the site reported in the present study, of SARS-CoV-2 proteinase with Lys5 as common interacting residue (Iftikhar et al. 2020) . The Ramachandran plot of the docked M pro -quercetin complex in turn displayed high stability of the M pro -quercetin complex as most of the residues lie in favourable regions Page 5 of 10 67 (Fig. 1) . MD simulations each of 50 ns were performed for the best docked SARS-CoV-2 M pro -quercetin complex (protein-ligand complex), SARS-CoV-2 M pro (protein) and quercetin (ligand) using GROMACS software (Pronk et al. 2013; Abraham et al. 2015) . The MD simulation trajectory of SARS-CoV-2 M pro -quercetin complex was analysed and interaction images at different time points are presented in Fig. 2 . The RMSD, RMSF, Rg and hydrogen bonds plots were generated. The RMSD values of both SARS-CoV-2 M pro and SARS-CoV-2 M pro -quercetin complex fluctuated within a range of 0.15 nm with respect to time thus displayed high stability of protein-ligand complex during the simulation (Fig. 3a ) (Kato et al. 2017) . The variations in Rg of the M pro -quercetin complex and M pro alone, during simulation period showed similar fashion (Fig. 3b) . Very low variation in Rg values with respect to time revealed consistent compactness of the protein during simulation even when present in the complexed form with quercetin, hence proved high stability of the complex (Lobanov et al. 2008) . The root mean square fluctuations calculated for the protein backbone of M pro -quercetin complex were found to be similar to that of the reference M pro protein during simulation thus provide evidence that the integrity of protein backbone was maintained and least damaging displacements in backbone were encountered by the protein when complexed with quercetin (Fig. 3c) . Along with this, the high amplitude of fluctuations at the terminal residues were due to the high degree of freedom of these residues compared to that of the intermediate region of the protein chain (Martínez 2015) . Further, the fluctuations in residues of loop regions were higher compared to other forms of secondary structures. The number of hydrogen bonds formed between quercetin and M pro protein in complex form during the simulation was plotted which proved effective interaction of quercetin with the M pro protein during the time period of simulation (Fig. 3d) . The higher number of hydrogen bond display significant interaction between the quercetin and M pro during the simulation. All these results proved highly stable binding of quercetin with SARS-CoV-2 M pro , thus providing great scope for its potential to modulate the structural conformation of M pro to a greater extent.

To study the allosteric inhibitory effect of quercetin against M pro , the substrate peptide (AVLQSGFR) sequence containing the cleavage site of M pro was modelled and Ramachandran plot was generated to check the stability of structure (Fig. 4) . The structures of SARS-CoV-2 M pro and SARS-CoV-2 M pro -quercetin complex achieved different conformations at different time points during MD simulation. The structural conformations were extracted from MD simulation trajectory at 1 ns, 10 ns, 20 ns, 30 ns, 40 ns and 50 ns. Molecular docking of these conformations with the modelled substrate peptide at the active site of M pro was performed. The docking analysis of M pro conformations with substrate peptide at enzyme's active site showed insignificant conformational change in active site at different time points (Fig. 5) . In turn, cleavage of peptide was noticed when docked with 1 ns, 10 ns, 20 ns, 30 ns, 40 ns and 50 ns conformations proving significant activity of the enzyme in absence of quercetin. The ZDock scores of cleaved peptide with M pro conformations were in the range of 233-137 due to reduce surface area of contact of cleaved peptide fragments with the enzyme. In case of docking analysis of Page 6 of 10 M pro -quercetin complex conformations with substrate peptide at M pro active site, significant changes in the active site conformation were observed at 1 ns, 10 ns, 20 ns, 30 ns, 40 ns and 50 ns (Fig. 6) . The leading active site residues (His41, Cys145, Asp187, Arg188, Gln189 and Thr190) were found to be displaced resulting in change in area of substrate binding cavity thus no peptide cleavage occurred in the presence of quercetin. No interaction of the substrate peptide was noticed with the leading active site residue Cys145 and His41 in the presence of quercetin, as was observed when M pro was docked with substrate peptide in the absence of quercetin. This was due to significant dislocation of residues from their native positions and opening of loops containing Cys45 and His41. The peptide was found intact at the active site in 1 ns, 10 ns, 20 ns and 30 ns conformations with significant reduction in ZDock scores from 785 to 158. In 40 ns and 50 ns conformations, the severe changes in the position of the loop containing Asp187, Arg188, Gln189 and Thr190 and its penetration into the cavity along with the change in relative conformation of beta sheets caused reduction in the affinity of the substrate peptide for the active site and no bonding of substrate was found at the catalytic site resulting in low ZDock scores of 114 and 86, respectively. The significant changes in the positioning of active site residues and reduction in the exposed area due to perturbations induced by quercetin interactions at an allosteric site are apparent by the M pro -quercetin 3D conformations, which were extracted at a different time point during the simulation period of 50 ns, docked with substrate peptide (Fig. 6 ). Many researches have reported active site inhibitors for SARS-CoV-2 M pro which compete with the substrate polypeptide to occupy the active site (Ma et al. 2020; Zhang et al. 2020; Zhu et al. 2011; Jin et al. 2020a, b; Mengist et al. 2020 ). An adequate concentration of such inhibitors, greater than the substrate concentration, is often required so that the inhibitor can compete with substrate and can effectively bind to the enzyme's active site thereby blocking the binding of the substrate (Strelow et al. 2012) . Ma et al. in their in vitro study, reported GC376, boceprevir and calpain inhibitors as potent active site inhibitors of SARS-CoV-2 M pro . These (Jin et al. 2020a ). In another study, an antineoplastic agent carmofur was reported to inhibit viral replication by targeting the M pro at its active site with an EC50 value of 24.30 μM and IC50 value of 1.82 μM (Jin et al. 2020b ). The fact that effectiveness of competitive inhibitors or active site-binding inhibitors depends on substrate concentration and decreases with the increase in the concentration of substrate, put forward a demand for a concentration-independent inhibitor (Ramsay and Tipton 2017) . In our analysis, we revealed quercetin as an effective allosteric inhibitor which can bind to the M pro even when the active site is occupied by the substrate and will not be affected by any change in substrate concentration (Ramsay and Tipton 2017; Robinson 2015) . Although very few studies have recently reported effective allosteric inhibitors against SARS-CoV-2 M pro , a topological analysis of the M pro revealed its very high sensitivity against any perturbation which it encounters and displayed excessively high tendency of the enzyme in transmitting it across the structure (Estradaa 2020; Doshi et al. 2016; Negre et al. 2018 ). These evidences can provide immense support to our present analysis of allosteric inhibition of M pro by quercetin where binding of the compound in consideration, at a different site modulates the active site to reduce its affinity for the substrate. Hence, all the above results and facts prove that quercetin by binding at an allosteric site in SARS-CoV-2 M pro can effectively inhibit the binding and thus the catalytic cleavage of substrate polypeptides 1a and 1ab. Such exclusive allosteric inhibition by a natural flavonol can massively hamper the replication of SARS-CoV-2 virus and can be developed into an exceedingly efficient and safe treatment against the current pandemic of COVID-19.

In the present in silico study, we analysed a natural flavonol quercetin for its effective allosteric inhibition against the SARS-CoV-2 M pro . This flavonol allosteric inhibitor can be effective in both free and substrate conjugated form of M pro . Binding energy of -8.31 kcal/mol for M pro -quercetin complex in molecular docking analysis provided evidence for effective binding of the compound with M pro . MD simulation of M pro -quercetin complex proved stable interactions of quercetin with M pro for a duration of 50 ns. In further molecular docking analysis of substrate peptide with M pro (in both quercetin conjugated and free form) proved a significant reduction in the binding affinity of M pro , when complexed with quercetin, for its substrate polypeptide. The docking ZDock scores of substrate peptide to M pro binding reduced from 1043 for native state (0 ns) M pro in the absence of quercetin to 86 for 50 ns conformation in the presence of quercetin. Henceforth, quercetin being a natural compound can be an efficient allosteric inhibitor against the effective target SARS-CoV-2 M pro . Furthermore, exhaustive research in this concern is highly demanded to develop a natural, safe and effective treatment for the extremely dreaded respiratory syndrome of COVID-19.

GROMACS: high performance molecular simulations through multi-level parallelism from laptops to supercomputers

Overviews of biological importance of quercetin: a bioactive flavonoid

The pharmacological activity, biochemical properties, and pharmacokinetics of the major natural polyphenolic flavonoid: quercetin

Dassault Systèmes, discovery studio modeling environment, release 2019

COVID-19) pandemic World Health Organisation

Dynamical network of residue-residue contacts reveals coupled allosteric effects in recognition, catalysis, and mutation

Allosteric inhibition of the SARS-CoV-2 main protease: insights from mass spectrometry based assays

Topological analysis of SARS CoV-2 main protease. Chaos: an Interdisciplinary

VMD: visual molecular dynamics

Identification of potential inhibitors of three key enzymes of SARS-CoV2 using computational approach

Structure of Mpro from SARS-CoV-2 and discovery of its inhibitors

Structural basis for the inhibition of SARS-CoV-2 main protease by antineoplastic drug carmofur

Validation of molecular dynamics simulations for prediction of three-dimensional structures of small proteins

Genomic characterization of a novel SARS-CoV-2

PEP-FOLD3: faster de novo structure prediction for linear peptides in solution and in complex

GROMACS 3.0: a package for molecular simulation and trajectory analysis

The structural basis for allosteric inhibition of a threonine-sensitive aspartokinase

Radius of gyration as an indicator of protein structure compactness

Boceprevir, GC-376, and calpain inhibitors II, XII inhibit SARS-CoV-2 viral replication by targeting the viral main protease

Automatic identification of mobile and rigid substructures in molecular dynamics simulations and fractional structural fluctuation analysis

The coronavirus nucleocapsid is a multifunctional protein

Designing of improved drugs for COVID-19: Crystal structure of SARS-CoV-2 main protease M pro

Eigenvector centrality for characterization of protein allosteric pathways

Role of changes in SARS-CoV-2 spike protein in the interaction with the human ACE2 receptor: an in silico analysis

Flavonoids: an overview

Accelerating protein docking in ZDOCK using an advanced 3D convolution library

ZDOCK server: interactive docking prediction of protein-protein complexes and symmetric multimers

GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit

Assessment of enzyme inhibition: a review with examples from the development of monoamine oxidase and cholinesterase inhibitory drugs

Enzymes: principles and biotechnological applications

Therapeutic potential of quercetin: new insights and perspectives for human health

Drug repurposing for candidate SARS-CoV-2 main protease inhibitors by a novel in silico method

CHARMM general force field: a force field for drug-like molecules compatible with the CHARMM all-atom additive biological force fields

Integrating atom-based and residue-based scoring functions for protein-protein docking

Drugs for allosteric sites on receptors

Reasoning of spike glycoproteins being more vulnerable to mutations among 158 coronavirus proteins from different species

Quercetin as an antiviral agent inhibits influenza A virus (IAV) entry

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

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

Peptide aldehyde inhibitors challenge the substrate specificity of the SARS-coronavirus main protease

Acknowledgements This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.