Allosteric Binding Sites of the SARS-CoV-2 Main Protease: Potential Targets for Broad-Spectrum Anti-Coronavirus Agents

Introduction

Since January 2020, the world has been afflicted by an unprecedented and emergent pandemic of coronavirus disease 2019 (COVID-19), which is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2).¹ The newly discovered coronavirus has spread to every country, with cases reported in Asia, Europe, North and South America, Australia, and Africa.² As of February 01, 2022, The World Health Organization (WHO) has reported 373,229,380 confirmed COVID-19 cases worldwide and 5,658,702 deaths (https://covid19.who.int/table). Besides SARS-CoV-2, coronaviruses from the betacoronavirus genus such as the severe acute respiratory syndrome coronavirus (SARS-CoV) and the Middle East respiratory syndrome coronavirus (MERS-CoV), have several cases reported previously in China and the Middle East, respectively.^3,4 However, COVID19 consequences on global health and the economy have made it an unprecedented public health crisis, prompting pharmaceutical companies and research institutes to devote significant time and resources to developing antiviral medications and vaccines as quickly as possible.⁵

The SARS-CoV-2 virus is a single-stranded positive-sense RNA genome enclosed within a membrane envelope belonging to the Coronaviridae family of viruses.^6,7 Similar to other viral infections, the entry of the genetic material into the host cell for replication and release of virions is a critical stage in the virus life cycle (Figure 1).⁸ The viral RNA is translated once the SARS-CoV-2 is internalized in the host cell resulting in virus-encoded proteins of various open reading frames (ORFs). Subsequently, the ORFs are translated into two viral replicase polyproteins (pp1a and pp1ab).⁹ The cleavage of these viral polyproteins into individual nonstructural proteins is essential for viral genome replication. With papain-like protease (PL^pro) cleaving two non-structural proteins at the N-terminus and main protease or 3C-Like protease (M^pro or 3CL^pro) recognizing 11 separate cleavage sites at the C-terminus, this cleaving process generates 16 mature non-structural proteins.^6,8–10 The remaining ORFs encode four structural proteins (spike, membrane, envelope, and nucleocapsid) and several accessory proteins.^11,12

Figure 1 The life cycle of the SARS-CoV-2 coronavirus inside the host cell. Created with BioRender.com.

Various important enzymes and proteins are involved in the viral replication and infectious capacity of the SARS-CoV-2 virus.¹³ Among them, two key proteases, such as main protease (M^pro) and papain-like protease (PL^pro), are required for viral replication.^14,15 PL^pro is known to cleave nsp1, nsp2, and nsp3, whereas M^pro cleaves the remaining 13 non-structural proteins (Figure 1).¹⁶ In addition, PL^pro is involved in antagonizing the host’s immune response upon viral infection through its deubiquitinating and deISGylating activities.^17,18 However, because both PL^pro and human deubiquitinases (DUBs) bind ubiquitin at the extended C-terminus with the consensus sequence Leu-X-Gly-Gly, PL^pro inhibitors may exhibit off-target effects on human DUBs.¹⁹ Hence, the issue of target selection must be addressed early on in the developmental process. The M^pro enzyme, on the other hand, is unique in that it exclusively cleaves polypeptides after a glutamine (Gln) residue, and no known human protease possesses the same cleavage specificity, making it a viable target for the development of highly selective antiviral agents.^20–23 With new COVID-19 variants and the susceptibility of RNA viruses for genetic reshuffling, mutations, and interspecies transmission, it has become more challenging to rely only on vaccines to effectively protect individuals and limit down the community spread of the infection.^24,25 Thus, small molecules that target conserved viral proteases like the main protease could lock key steps of the SARS-CoV-2 life cycle providing a broad-spectrum antiviral effect with minimal side effects.^20,26 In December 2021, the FDA issued emergency use authorization (EUA) for SARS-CoV-2 main protease (M^pro) inhibitor PAXLOVID (nirmatrelvir and ritonavir).²⁷ This is the first antiviral drug to be approved as an oral treatment of mild to moderate COVID-19 in patients who are at high risk of developing severe illness.²⁷

The M^pro enzyme encompasses several pockets on its surface that were proven to be important for its catalytic activity; some of them exist in distal areas from the main catalytic pocket, and others have been identified on the dimerization interface. The latter type, in particular, has shown promising results and seems to hinder the overall catalytic activity of the M^pro enzyme if a suitable organic molecule can successfully bind to it. Hence, in this review, we showcase all allosteric sites identified so far and illustrate all drug discovery efforts to target them.

M^pro Structure and Dimerization

The SARS-CoV-2 M^pro is a homodimer made up of two monomers, arranged almost perpendicular to one another.²³ Each monomer consists of three domains – domains I, II, and III. Domains I (residues 10–99) and II (residues 100–182) assume a chymotrypsin-like structure with a six-stranded antiparallel b-barrel fold.²³ Domains II and III are linked by a loop (residues 183–198) that takes on two different conformations: a long helix shape or a short helix shape.⁶

The active site is located in a cleft between the two N-terminal domains (Domain I and II) of the three-domain structure of the monomer.²³ The C-terminal helical domain III (residues 198–303) is involved in the regulation and dimerization of the enzyme via an intermolecular salt-bridge interaction between Glu290 and Arg4.^28,29 Upon dimerization and activation (Figure 2), both monomers orient themselves in the appropriate conformation to perform the catalytic function, and the N finger of each monomer interacts with Glu166 of the other monomer.²³ This step is vital for the enzyme’s catalytic activity as it shapes the S1 pocket of the substrate-binding site, which contains the catalytic dyad, His41, and Cys145.²⁵ Given the importance of dimerization, targeting the dimer site might be used to attenuate the catalytic activity.³⁰ Earlier MD simulations and mutational analyses studies revealed key residues involved in stabilizing the catalytically active dimeric structure of the SARS-CoV-2 M^pro enzyme, including Arg4, Ser10, Gly11, Glu14, Asn28, Ser139, Phe140, Ser147, Glu290, and Arg298.²⁵ These residues operate cooperatively to control dimerization and maintain the integrity of the dimer interface which is important for dimer stability.³¹ Several inhibitors have been identified to bind to the dimer site, inducing weak to moderate antiviral activity.^32,33 Thus, interfering with the residues on the dimer site is likely to deter the enzymes` catalytic activity and tendency to dimerize.³⁴

Figure 2 The Monomer (PDB: 6LU723) and Dimer (PDB:7CAM79) forms of the SARS-CoV-2 Mpro.

Allosteric Binding Sites on the M^pro Enzyme

Over the past two years, there have been extensive research efforts worldwide to develop new M^pro inhibitors, most of which targeted the enzyme catalytic pocket.²³ However, several inhibitors appeared to bind at sites distant from the catalytic pocket, which opens the door for an alternative approach to target the M^pro enzyme via an allosteric mechanism.³² Nevertheless, there are limited data on this area so far, which hinders the search for allosteric inhibitors.

Currently, two pockets of varying size and druggability properties have been identified as possible therapeutic targets: distal site and dimerization sites.^25,35,36 Computational studies have predicted that targeting either sites will result in allosteric modulation by inducing either conformational changes or thermal fluctuations around a fixed, mean conformation.^37–42 Given the difficulty associated with targeting the active site, investigating distal and dimer, allosteric sites will provide a better understanding of their characteristics and potential as molecular targets, ultimately assisting in the development of novel broad-spectrum inhibitors of SARS-CoV-2 M^pro.

Allosteric Sites on the Distal Region of the M^pro Enzyme

Analysis of the key contacts in the SARS-CoV-2 M^pro, which include Arg131-Thr199, Arg131-Asp289, Pro132-Thr196, Asp197, and Thr198-Asn238, Tyr239-Leu287, shows the presence of a pocket at the interface between domains II and III.⁴³ This pocket, which is located far from the active site, is known as the distal site. The distal site contains a number of residues surrounding the beta-sheet and alpha-helices. Due to the position of these residues, allosteric inhibitors targeting this site operate via reversible non-competitive inhibition.⁴⁴

Up to this point, no inhibitors of the distal site have been reported, suggesting that a druggability analysis would be necessary to determine the viability of the cavities distant from the catalytic pocket. When evaluated using the PockDrug⁴⁵ web server, these pockets druggability scores ranged from 0.37 to 1.00, 60% of which possess scores greater than 0.90.³⁶ As a result of its ability to bind drug-like molecules, the distal site provides a promising molecular target for potential allosteric inhibitors of the M^pro enzyme.

Further studies have looked into the coupled movement of residues from aggregate 1 μs MD simulations of the apo and dimer forms of M^pro to assess the effect of small molecules binding to the distal site on catalytic activity.³⁶ Regardless of whether the correlation is positive or negative, the distal and active sites appear to be dynamically dependent on one another in both the apo and dimer forms, demonstrating long-range communication between the distal and active sites.^36,46

Allosteric Sites on the M^pro Dimerization Interface

Located between domains II and III, the dimer site encompasses residues from the N finger and the C-terminal helix. These residues are essential for the dimerization and subsequent activation and stabilization of the protease.^47,48 The SARS-CoV-2 M^pro, similar to other viral proteases, is made up of two monomers that dimerize to produce an active homodimer. The dimerization process is vital for enzyme activity, as the monomer form of M^pro cannot fulfill the enzymatic role in catalysis.⁴⁹

Due to the complexities of dimer dynamics, it is frequently assumed that one monomer is active in the dimer while the other is inactive.⁵⁰ Protease activity can be disrupted by mutating key residues in the dimer interface.^25,51 These residues operate cooperatively to control dimerization and maintain dimer stability.³¹ Mutations of Arg298 in the dimer interface, in particular, are detrimental to viral activity because they prevent dimerization and inactivate the enzyme irreversibly.⁴⁷ Hence, compounds that block dimerization may operate as allosteric inhibitors of protease activity by impairing its catalytic activity, inducing weak to moderate antiviral activity.^52,53

Given the allosteric linkage between the dimerization site and the catalytic region, targeting protease dimerization could potentially impact the substrate pocket and thus inhibit M^pro activity.^54,55 x1187 and x1086 are two fragments that disrupt dimer formation by binding to hydrophobic pockets at the dimer site, as illustrated in Figure 3.^32,44,51 By disrupting dimerization, these fragments can be used to develop more potent drugs that allosterically modulate M^pro activity. Such fragments can induce antiviral activity by disrupting the dimer interface and conserving M^pro in its inactive monomer state.^44,56 Jiménez-Avalos et al⁵³ have found that the binding of CHEMBL2171598 at the dimer site prevents subsequent dimerization and catalytic activity, mirroring the effect of the detrimental mutation Arg298Ala.

Figure 3 (A) The Mpro enzyme in the dimer form with co-crystallized ligands, (B) x1086 (PDB:5RGQ51) and (C) x1187 (PDB:5RFA51) at the dimer interface indicates opportunities for allosteric modulation. For clarity, superimposition on the dimer form of Mpro (PDB:7CAM79) was performed to demonstrate binding onto the dimerization site.

Due to the fact that dimerization stabilizes both the protein and its active site, apo-structures tend to lack any distinctive pockets at the dimer interface.⁵³ However, when preserved in the monomer state, the protease features a well-defined pocket with desirable size and polarity that may bind ligand-like molecules. As a result, ligand-bound structures that maintain the monomer state are more likely to produce a druggable pocket at the dimer interface than apo-structures. Investigating these pockets represents a promising strategy for inhibiting the M^pro via allosteric means.

Druggability Assessment of Allosteric Sites

As multiple allosteric pockets have been identified experimentally,^32,51 it was necessary to perform a proper evaluation of these sites and assess them in terms of druggability or the ability of the target site to recruit and bind drug-like ligands. At the time of writing this review, six allosteric sites were experimentally published along with eight co-crystallized ligands or fragments, as shown in Figure 4.^32,51 Accordingly, we run a well-established module in the Schrodinger⁵⁷ software that is called SiteMap.⁵⁸ In this tool, we can identify the co-crystallized ligand as the center of the pocket that has to be assessed. Druggability scores (Dscore) will be then generated for each of these pockets; the higher the scoring, the better the druggability (typically pockets with a Dscore ≥1.0 are considered very druggable, whereas sites with a Dscore less than 0.8 are classified as difficult non-drug binding sites).⁵⁹ Such scoring is calculated according to three major factors: pocket size, hydrophobicity, and enclosure.⁵⁹

Figure 4 Representation of the Mpro enzyme in the monomer form labeled with its six experimentally proven allosteric sites, along with their co-crystallized ligands and Dscores.

Interestingly, we found that only two sites have shown very druggable pockets (ie, site #2 and site #5), while one showed poor druggability (ie, site #3), and the other three sites could not be detected at all by SiteMap (ie, site #1, site #4 and site #6), obviously because of their very shallow nature.

From Table 1, it can be seen that the highest-ranking site, site #2, possesses a large pocket (n = 102 spheres) which is well-defined (e = 0.71) with high hydrophilicity (p = 1.083), as illustrated in Figure 5. Next to site 2, also at the dimer interface, is another “druggable” site, site #5, which has comparable features as it is a large (n = 100 spheres), less defined cavity (e = 0.71), and lower hydrophilicity (p = 0.92). Similarly, site #3 is also located at the dimer interface; however, it is considered a difficult target due to its very small (n = 14 spheres), and shallow cavity (e = 0.46) that is significantly less hydrophilic (p = 0.433). Finally, no well-defined pocket was detected by SiteMap for sites #1, #4, and #6, as they showed a relatively flat surface. Visual inspection of these sites confirms their superficial character as they appear on the interface (Figure 5). With protease movement, these sites may become more exposed; nevertheless, further dynamic analysis is required. Hence, out of the experimentally found allosteric sites, only sites #2 and #5 seem to be promising and able to bind druglike inhibitors.

Table 1 Druggability Assessment of the Six Experimentally Proven Allosteric Sites, Along with Their Co-Crystallized Ligands and Dscores. (ND: Not Detected)

Figure 5 The overview shows the surface of the Mpro and ligand-receptor interaction diagrams across all reported allosteric sites. (A) Site 1 (x0390; PDB:5REC51); (B) Site 2 (AT7519; PDB:7AGA32); (C) Site 3 (Pelitinib; PDB:7AXM32); (D) Site 4 (x0425; PDB:5RGJ51); (E) Site 5 (x1086; PDB:5RGQ51); (F) Site 6 (x0887; PDB:5RF051).

Recent Advances in the Development of Allosteric SARS-CoV-2 M^pro Inhibitors

Several strategies were employed to find therapeutic agents targeting the M^pro enzyme, which involve screening approaches, such as high throughput virtual screening, mass spectrometry, X-ray crystallographic, electrophilic, and NMR fragment screening.^32,44,51,60 Most discovered inhibitors block the active site and bind covalently to its nucleophilic cysteine residue. Interestingly, few appeared to bind at sites distant to the catalytic pocket, opening the door for an alternative approach that possibly targets the M^pro via an allosteric mechanism. It is important to note that only one attempt has currently been made to target the allosteric site through both computational and experimental means, which is discussed in the section.

Experimental Studies

Limited work has been done to evaluate reported M^pro allosteric inhibitors in vitro as shown in Figure 6. Gunther et al³² used X-Ray screening to identify drugs that bind at the M^pro active site. Surprisingly, they have identified five compounds that bind at two distinct allosteric sites. The first allosteric site (site #3 in our paper) binds five compounds: pelitinib, RS-102895, Ifenprodil, PD-168568, and tofogliflozin. Pelitinib, an anti-cancer drug, has shown considerable anti-viral activity with EC₅₀ of 1.25 μM, outperforming other inhibitors in this study as the top-ranking allosteric inhibitor. The remaining four drugs had little to no anti-viral activity. The second allosteric site is found near the dimerization site (site #5 in our analysis), binds CDK inhibitor AT-7519, and has exhibited weak anti-viral activity with EC₅₀ of 25.2 μM. Although this validates their anti-viral activity and affinity for allosteric sites, they have not been experimentally tested for their enzyme inhibition property. When evaluating the binding mode of the co-crystallized ligands AT-7519 (Figure 7), this ligand shows a reasonable fitting; however, many of its functional groups seem to have no role in binding (in line with its experimentally proven weak activity). Similarly, Pelitinib (Figure 7) exhibits a similar binding pattern and weak interactions. Moreover, both described sites are in the dimer site, which is essential for protease dimerization and activation.

Figure 6 Chemical structure of experimentally evaluated allosteric Mpro inhibitor.

Figure 7 (A) The Mpro enzyme in the dimer form with co-crystallized ligands, (B) AT-7519 (PDB:7AGA32) and (C) Pelitinib (PDB:7AXM32) at the dimer interface indicates opportunities for allosteric modulation. For clarity, superimposition on the dimer form of Mpro (PDB:7CAM79) was performed to demonstrate binding onto the dimerization site.

Broad-spectrum antiviral agents, chebulagic acid and punicalagin are another set of compounds that were previously identified as non-competitive M^pro inhibitors. Such compounds demonstrated a non-druglike character as they exhibit large sizes with EC₅₀ and IC₅₀ values in the micromolar range.³³ There is no information about their binding site, but they are suspected to bind near the dimer site. Nonetheless, these results demonstrate that chebulagic acid and pelitinib are reversible inhibitors that bind non-competitively at distinct allosteric sites on the M^pro.^32,33

Another notable drug discovery strategy is fragment screening, which uses small molecular fragments to probe new cavities on the SARS-CoV-2 M^pro interface. Depending on their structural features, these cavities serve as potential targets for therapeutic development. So far, electrophilic and noncovalent fragment screening appears to be the most effective and informative strategy for defining these sites.^32,51,60

Douangamath et al⁵¹ used a combination of mass spectrometry and crystallographic techniques, resulting in six M^pro complexes with co-crystallized allosteric fragments (Figure 8), two of which bind at the dimerization site (x1086 and x1087), and four of which bind in distinct hydrophilic allosteric cavities (x0390, x0464, x0425, and x0887). Recently, NMR Fragment-based screening found fragment F15 to bind at the M^pro dimer interface, where it is predicted to mimic the binding of fragments x1086 and x1187 in the hydrophobic pocket.⁶⁰

Figure 8 The so far known co-crystallized fragments bound to different allosteric sites after superposing them onto a single crystal structure of the SARS-CoV-2 Mpro.

Since there is no testing performed, it is highly predictable that these small molecular fragments have little to no antiviral activity; yet, they form multiple interactions with their respective binding site, shown in Figure 7, which boosts the binding energy. Researchers can benefit from these interactions by using the fragment as the lead compound, which could be later optimized to derive a series of small molecules with higher affinity and in-vitro and in-vivo antiviral activity.

Computational Studies

Given the similarities between the SARS-CoV-2 M^pro and other viral cysteine proteases, repurposing FDA-approved drugs as anti-COVID-19 therapeutics provides a time and cost-effective alternative to novel drug design.²⁰ The initial concept was to re-propose current antiviral medications used in the treatment of HIV, hepatitis B and C, influenza, and the common cold, as they have shown to be effective in treating past coronavirus infections.⁶¹ A few of these drugs, listed in Table 2, appear to bind at high binding energies to sites distant from the active site, indicating a preference for allosteric sites.

Table 2 List of Repurposed Drugs That Bind to Allosteric Sites on the SARS-CoV-2 Mpro

Ebselen, a synthetic organoselenium molecule, exhibits higher binding energy for the potential allosteric site (−8.87 kcal/mol) compared to the active site (−5.55 kcal/mol). When ebselen binds to the site between the II and III domains, it exerts a significant allosteric effect that limits catalytic site access via surface-loop interactions.⁶²

The antiviral drug elbasvir has a higher docking score at both potential allosteric sites (−10.8 kcal/mol and −11.1 kcal/mol) compared to the active site (−8.8 kcal/mol), as well as a higher binding affinity towards both allosteric sites (−51.08 kcal/mol and −42.37 kcal/mol), making it an ideal candidate for clinical evaluation.⁶³ However, it is important to remember that binding at higher energy does not necessarily mean a higher binding affinity for the allosteric sites. Some drugs, such as dihydroergotamine, bromocriptine, and diosmin, have higher binding energies for potential allosteric sites, but a higher binding affinity for the active site.⁶³

Although drug repurposing has frequently shown to be a successful approach, the majority of these studies were carried out computationally through virtual screening of existing drugs, meaning that they require additional in-vitro and in-vivo testing for viral activity. The lack of in-vitro and in-vivo testing data makes it difficult to draw firm conclusions about whether the presupposing drug approach acts on the M^pro or other biological pathways involved with disease progression and symptoms.

Advantages of Targeting Allosteric Sites on the SARS-CoV-2 M^pro

Allosteric modulation has proven to be a viable alternative to orthosteric ligands, contributing significantly to the development of GPCR and protein kinase modulators.⁶⁴ As targeting the SARS-CoV-2 active site appears to be a challenging endeavor, allosteric sites have emerged as attractive sites for the design of small-molecule inhibitors, offering several advantages over their orthosteric counterpart.^25,35,36 All allosteric inhibitors identified thus far act via reversible means, making them potentially safer than many of the previously discovered M^pro inhibitors, which contain reactive groups and form a covalent bond with the enzyme, and hence possibly with other enzymes.^65,66 Therefore, the M^pro allosteric sites appear to be an interesting target for the discovery of broad-spectrum anti-corona viral agents that can act alone or in combination with the other M^pro competitive inhibitors.

Dimerization-Targeted Allosteric Inhibition

Acting on the dimerization site seems to be rather promising as this strategy has previously been reported to modulate several viral enzymes, such as HIV reverse transcriptase, integrase, and herpes simplex virus.⁶⁷ The dimeric nature of SARS-CoV-2 M^pro allows targeting dimerization at certain hot spots by small molecule inhibitors, which seems to disrupt and inactivate the enzyme.⁶⁸ The dimerization interface, located between domains II and III, is densely packed with high-energy hot spots.⁶⁹ A few of these residues are essential for M^pro dimerization (Ser1, Arg4, His172, and Glu290), as they form salt bridges between the two monomers.^28,69,70

Beyond promoting dimer dissociation, the mechanism by which M^pro dimerization inhibitors interfere with the catalytic activity is not entirely clear.³⁰ Some studies have linked enzyme inhibition to the compound’s tendency to induce dimer dissociation.²⁵ Others have proposed the presence of an allosteric connection between the sizes of the catalytic and allosteric dimerization sites.⁷¹ AMOEBA trajectory analysis verifies this, revealing that the studied structures had an open catalytic site and a small allosteric dimerization site, and vice versa.⁷¹ It also implies that although inhibiting M^pro dimerization disrupts its activation, it does not necessarily block its interaction with the substrate. Nonetheless, current findings provide encouraging data on dimerization-targeted M^pro allosteric inhibition.

Maximizing Antiviral Activity Through Synergy

As previously stated, inhibiting M^pro dimerization does not necessarily hamper its interaction with the substrate,³⁰ as there will always be other dimeric molecules that will exert their enzymatic activity upon substrate binding. This, however, hints at the possibility of synergism, in which a combination of active and allosteric site inhibitors can be employed to induce maximum inhibition of M^pro activity.^30,72

In fact, a previous study investigated the potential synergistic effect of combining allosteric and active site inhibitors.⁷² When tested alone, Quinacrine demonstrated competitive inhibition with an IC₅₀ of 6.3 μM, while suramin exhibited non-competitive inhibition with an IC₅₀ of 7.8 μM. Interestingly, combining them together has successfully boosted the M^pro inhibition activity by more than ten folds (new IC₅₀ = 0.46 μM).⁷² Thus, M^pro allosteric inhibitors seem to have the potential to act alone or in conjunction with their competitive peers, to allow lower doses and hence greater safety for the COVID19 patients.^30,72

Conservation of Allosteric Sites

Targeting conserved regions on the protease can vastly enhance treatment efficiency.^73,74 The M^pro is distinct from other viral proteins that it is highly conserved across the coronavirus family, as any mutation can be detrimental to the virus and its activity.^75,76 However, the active site is of high plasticity and is regulated by a flexible C44-P52 loop, raising concerns that any mutations to its amino acid sequence may impede small-molecule entry.⁷⁷ Accordingly, resistance mutations to these active site inhibitors should be expected. Structural changes that confer resistance may occur in any site throughout the protease and disrupt the complex interactions required for enzymatic function.⁴⁵ Hence, targeting multiple regions on the protease, including allosteric sites, might be a better strategy than using a single compound that targets the active site.

Most of the residues implicated in allosteric sites were reported to be conserved in most SARS-CoV-2 structures studies. This could be interpreted as a reassuring sign of the viability of the newly identified allosteric targets, as potential drugs targeting these sites are assumed to be robust against mutations.^46,78 Moreover, hydrophobic residues at the allosteric sites are substantially conserved throughout the M^pro of all human coronaviruses, indicating that mutations are unlikely.⁴³ As a result, the likelihood of mutation-mediated drug resistance is quite low, and inhibitors will have broad antiviral activity.³⁰

Conclusions and Future Prospective

This review paper summarizes the work done so far on the M^pro allosteric site and it illustrates the potential of such pockets to act as a drug target (Figure 5). There have been several attempts to target the M^pro allosteric sites; however, many of these are computational and do not provide us with very conclusive data. Having said that, a few experimental works have confirmed the presence of multiple pockets, at least two of them appear to have a reasonable druggability potential (as per an in-house computational assessment). Both pockets are located at the dimerization site or in its vicinity, offering more potential to interfere with the M^pro catalytic site through a dimerization-based allosteric mechanism.

As a future perspective, computer-aided drug design seems to be in a good position for the discovery of new inhibitors, especially since several crystal structures have been already revealed, some of them co-crystallized with small ligands in their allosteric sites. For instance, compounds x1086 and x1187 both together exhibited a complementary binding mode in the allosteric site #5, having their sulfonamide and the hetero aromatic ring sitting in two distinct subpockets (Figure 9). Based on the aforementioned observations, Site #5 can be targeted by docking-guided derivatization through synthesizing a series of compounds designed based on both x1086 and x1187, where distinguishingly bound moieties will be joined together. To sum up, druglike allosteric inhibitors seem to be achievable and can provide us with an extra weapon in our war against the COVID-19 disease, and possibly against other coronaviruses-caused pandemics.

Figure 9 (A) The binding mode of compounds x1086 (pink) and x1187 (blue) aligned on top of each other inside the allosteric site #5, (B) 2D representation of compounds x1086 (pink) and x1187 (blue) aligned on top of each other.