The N-terminal octapeptide acts as a dimerization inhibitor of SARS coronavirus 3C-like proteinase

Competitive Interaction of the SGFRKMAF Peptide with 3CLpro Dimerization Intermediates: A Brownian Dynamics Investigation

The SGFRKMAF peptide is known to inhibit the dimerization of 3CLpro monomers, which is essential for SARS-CoV-2 replication. The mechanism behind this, however, is largely unknown. In this work, we used Brownian dynamics simulations to compare and contrast 3CLpro monomer-monomer interactions and 3CLpro monomer-SGFRKMAF peptide interactions. We found that formation of the 3CLpro wild-type dimer could potentially involve formation of three intermediates that are primarily stabilized by G11-G124, S1-S301, and T118-G278 interactions. Analysis of 3CLpro monomer interaction with the SGFRKMAF peptide, however, revealed the presence of eight basins of interactions where the peptide assumes the highest local densities at the 3CLPro monomer surface. The second highest-density basin was found to coincide with the interface region of the wild-type 3CLpro dimer, thereby directly blocking the 3CLpro dimer-dimer interactions. The other basins, however, were found to lie far from the interface region. Notably, we found that only 6% of the BD trajectories end up directly into the basin at the interface region and ∼39% of the trajectories end up into those basins lying away from the interface region, indicating a greater role for peptide binding at sites away from the dimer interface region. Importantly, the locations of the basins lying away from the interface were found to coincide with the 3CLpro residues involved in stabilization of the 3CLpro monomer-monomer intermediates. Given that the rate constant of the peptide reaching the monomer surface was found to be almost an order of magnitude higher than the rate constant of monomer-monomer association, the SGFRKMAF peptide has the potential to inhibit dimerization of 3CLpro monomers not only through blocking the interface region but also through blocking the formation of the intermediates involved in the dimerization process. This could potentially open new avenues for 3CLpro dimerization inhibitors that transcend traditional X-ray-based discovery approaches.

SARS-CoV-2 Mpro oligomerization as a potential target for therapy

The main protease (Mpro) of SARS-CoV-2 is critical in the virus's replication cycle, facilitating the maturation of polyproteins into functional units. Due to its conservation across taxa, Mpro is a promising target for broad-spectrum antiviral drugs. Targeting Mpro with small molecule inhibitors, such as nirmatrelvir combined with ritonavir (Paxlovid™), which the FDA has approved for post-exposure treatment and prophylaxis, can effectively interrupt the replication process of the virus. A key aspect of Mpro's function is its ability to form a functional dimer. However, the mechanics of dimerization and its influence on proteolytic activity remain less understood. In this study, we utilized biochemical, structural, and molecular modelling approaches to explore Mpro dimerization. We evaluated critical residues, specifically Arg4 and Arg298, that are essential for dimerization. Our results show that changes in the oligomerization state of Mpro directly affect its enzymatic activity and dimerization propensity. We discovered a synergistic relationship influencing dimer formation, involving both intra- and intermolecular interactions. These findings highlight the potential for developing allosteric inhibitors targeting Mpro, offering promising new directions for therapeutic strategies.

On Inactivation of the Coronavirus Main Protease

A deeper understanding of the inactive conformations of the coronavirus main protease (MPro) could inform the design of allosteric drugs. Based on extensive molecular dynamics simulations, we built a Markov State Model to investigate structural changes that can inactivate the SARS-CoV-2 MPro. In a subset of structures, one subunit of the homodimer assumes an inactive conformation that resembles an inactive crystal structure. However, contradicting the widely held half-of-sites activity hypothesis, the most populated enzyme structures have two active subunits. We then used transition path theory (TPT) and the Jensen-Shannon Divergence (JSD) to pinpoint residues involved in the inactivation process. A π stack between Phe140 and His163 is a key feature that can distinguish active and inactive conformations of MPro. Each subunit has unique inactive conformations stabilized by π stacking interactions involving residues Phe140, Tyr118, His163, and His172, a hydrogen bonding network centered around His163 and His172, and a modified network of interactions in the dimer interface. The importance of these residues in maintaining an active structure explains the sensitivity of enzymatic activity to site-directed mutagenesis.

AI-driven covalent drug design strategies targeting main protease (mpro) against SARS-CoV-2: structural insights and molecular mechanisms

The emergence of new SARS-CoV-2 variants has raised concerns about the effectiveness of COVID-19 vaccines. To address this challenge, small-molecule antivirals have been proposed as a crucial therapeutic option. Among potential targets for anti-COVID-19 therapy, the main protease (Mpro) of SARS-CoV-2 is important due to its essential role in the virus's life cycle and high conservation. The substrate-binding region of the core proteases of various coronaviruses, including SARS-CoV-2, SARS-CoV, and Middle East respiratory syndrome coronavirus (MERS-CoV), could be used for the generation of new protease inhibitors. Various drug discovery methods have employed a diverse range of strategies, targeting both monomeric and dimeric forms, including drug repurposing, integrating virtual screening with high-throughput screening (HTS), and structure-based drug design, each demonstrating varying levels of efficiency. Covalent inhibitors, such as Nirmatrelvir and MG-101, showcase robust and high-affinity binding to Mpro, exhibiting stable interactions confirmed by molecular docking studies. Development of effective antiviral drugs is imperative to address potential pandemic situations. This review explores recent advances in the search for Mpro inhibitors and the application of artificial intelligence (AI) in drug design. AI leverages vast datasets and advanced algorithms to streamline the design and identification of promising Mpro inhibitors. AI-driven drug discovery methods, including molecular docking, predictive modeling, and structure-based drug repurposing, are at the forefront of identifying potential candidates for effective antiviral therapy. In a time when COVID-19 potentially threat global health, the quest for potent antiviral solutions targeting Mpro could be critical for inhibiting the virus.Communicated by Ramaswamy H. Sarma.

A naturally occurring G11S mutation in the 3C-like protease from the SARS-CoV-2 virus dramatically weakens the dimer interface

The 3C-like protease (3CLpro ) is crucial to the replication of SARS-CoV-2, the causative agent of COVID-19, and is the target of several successful drugs including Paxlovid and Xocova. Nevertheless, the emergence of viral resistance underlines the need for alternative drug strategies. 3CLpro only functions as a homodimer, making the protein-protein interface an attractive drug target. Dimerization is partly mediated by a conserved glycine at position 11. However, some naturally occurring SARS-CoV-2 sequences contain a serine at this position, potentially disrupting the dimer. We have used concentration-dependent activity assays and mass spectrometry to show that indeed the G11S mutation reduces the stability of the dimer by 600-fold. This helps to set a quantitative benchmark for the minimum potency required of any future protein-protein interaction inhibitors targeting 3CLpro and raises interesting questions regarding how coronaviruses bearing such weakly dimerizing 3CLpro enzymes are capable of replication.

Structure and function of SARS-CoV and SARS-CoV-2 main proteases and their inhibition: A comprehensive review

Severe acute respiratory syndrome-associated coronavirus (SARS-CoV) identified in 2003 infected ∼8000 people in 26 countries with 800 deaths, which was soon contained and eradicated by syndromic surveillance and enhanced quarantine. A closely related coronavirus SARS-CoV-2, the causative agent of COVID-19 identified in 2019, has been dramatically more contagious and catastrophic. It has infected and caused various flu-like symptoms of billions of people in >200 countries, including >6 million people died of or with the virus. Despite the availability of several vaccines and antiviral drugs against SARS-CoV-2, finding new therapeutics is needed because of viral evolution and a possible emerging coronavirus in the future. The main protease (Mpro) of these coronaviruses plays important roles in their life cycle and is essential for the viral replication. This article represents a comprehensive review of the function, structure and inhibition of SARS-CoV and -CoV-2 Mpro, including structure-activity relationships, protein-inhibitor interactions and clinical trial status.

Phytonutrient Inhibitors of SARS-CoV-2/NSP5-Encoded Main Protease (Mpro) Autocleavage Enzyme Critical for COVID-19 Pathogenesis

The genomic reshuffling, mutagenicity, and high transmission rate of the SARS-CoV-2 pathogen highlights an urgent need for effective antiviral interventions for COVID-19 control. Targeting the highly conserved viral genes and/or gene-encoded viral proteins such as main proteinase (M^pro), RNA-dependent RNA polymerase (RdRp) and helicases are plausible antiviral approaches to prevent replication and propagation of the SARS-CoV-2 infection. Coronaviruses (CoVs) are prone to extensive mutagenesis; however, any genetic alteration to its highly conserved M^pro enzyme is often detrimental to the viral pathogen. Therefore, inhibitors that target the M^pro enzyme could reduce the risk of mutation-mediated drug resistance and provide effective antiviral protection. Several existing antiviral drugs and dietary bioactives are currently repurposed to treat COVID-19. Dietary bioactives from three ayurvedic medicinal herbs, 18 β-glycyrrhetinic acid (ΔG = 8.86 kcal/mol), Solanocapsine (ΔG = 8.59 kcal/mol), and Vasicoline (ΔG = 7.34 kcal/mol), showed high-affinity binding to M^pro enzyme than the native N3 inhibitor (ΔG = 5.41 kcal/mol). Flavonoids strongly inhibited SARS-CoV-2 M^pro with comparable or higher potency than the antiviral drug, remdesivir. Several tannin hydrolysates avidly bound to the receptor-binding domain and catalytic dyad (His₄₁ and Cys₁₄₅) of SARS-CoV-2 M^pro through H-bonding forces. Quercetin binding to M^pro altered the thermostability of the viral protein through redox-based mechanism and inhibited the viral enzymatic activity. Interaction of quercetin-derivatives with the M^pro seem to be influenced by the 7-OH group and the acetoxylation of sugar moiety on the ligand molecule. Based on pharmacokinetic and ADMET profiles, several phytonutrients could serve as a promising redox nutraceutical for COVID-19 management.

Venetoclax: a promising repurposed drug against SARS-CoV-2 main protease

Global health care emergency caused by a new coronavirus (severe acute respiratory syndrome coronavirus 2 or SARS-CoV-2) demands urgent need to repurpose the approved pharmaceutical drugs. Main protease, M^pro of SARS-CoV-2 draws significant attention as a drug target. Herein, we have screened FDA approved organosulfur drugs (till 2016) and our laboratory synthesized organosulfur and organoselenium compounds (L1-L306) against M^pro-apo using docking followed by classical MD simulations. Additionally, a series of compounds (L307-L364) were chosen from previous experimental studies, which were reported to exhibit inhibitory potentials towards M^pro. We found several organosulfur drugs, particularly Venetoclax (FDA approved organosulfur drug for Leukemia) to be a high-affinity binders to the M^pro of SARS-CoV-2. The results reveal that organosulfur compounds including Venetoclax preferentially bind (non-covalently) to the non-catalytic pocket of the protein located in the dimer interface. We found that the ligand binding is primarily favoured by ligand-protein van der Waals interaction and penalized by desolvation effect. Interestingly, Venetoclax binding alters the local flexibility of M^pro and exerts pronounced effect in the C-terminal as well as two loop regions (Loop-A and Loop-B) that play important roles in catalysis. These findings highlighted the importance of drug repurposing and explored the non-catalytic pockets of M^pro in combating COVID-19 infection in addition to the importance of catalytic binding pocket of the protein.Communicated by Ramaswamy H. Sarma.

Structural differences in 3C-like protease (Mpro) from SARS-CoV and SARS-CoV-2: molecular insights revealed by Molecular Dynamics Simulations

Novel coronavirus SARS-CoV-2 has infected millions of people with thousands of mortalities globally. The main protease (Mpro) is vital in processing replicase polyproteins. Both the CoV's Mpro shares 97% identity, with 12 mutations, but none are present in the active site. Although many therapeutics and vaccines are available to combat SARS-CoV-2, these treatments may not be practical due to their high mutational rate. On the other hand, Mpro has a high degree of conservation throughout variants, making Mpro a stout drug target. Here, we report a detailed comparison of both the monomeric Mpro and the biologically active dimeric Mpro using MD simulation to understand the impact of the 12 divergent residues (T35V, A46S, S65N, L86V, R88K, S94A, H134F, K180N, L202V, A267S, T285A and I286L) on the molecular microenvironment and the interaction between crucial residues. The present study concluded that the change in the microenvironment of residues at the entrance (T25, T26, M49 and Q189), near the catalytic site (F140, H163, H164, M165 and H172) and in the substrate-binding site (V35, N65, K88 and N180) is due to 12 mutations in the SARS-CoV-2 Mpro. Furthermore, the involvement of F140, E166 and H172 residues in dimerization stabilizes the Mpro dimer, which should be considered. We anticipate that networks and microenvironment changes identified here might guide repurposing attempts and optimization of new Mpro inhibitors.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11224-022-02089-6.

The SARS-CoV-2 main protease (Mpro): Structure, function, and emerging therapies for COVID-19

The main proteases (Mpro), also termed 3-chymotrypsin-like proteases (3CLpro), are a class of highly conserved cysteine hydrolases in β-coronaviruses. Increasing evidence has demonstrated that 3CLpros play an indispensable role in viral replication and have been recognized as key targets for preventing and treating coronavirus-caused infectious diseases, including COVID-19. This review is focused on the structural features and biological function of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main protease Mpro (also known as 3CLpro), as well as recent advances in discovering and developing SARS-CoV-2 3CLpro inhibitors. To better understand the characteristics of SARS-CoV-2 3CLpro inhibitors, the inhibition activities, inhibitory mechanisms, and key structural features of various 3CLpro inhibitors (including marketed drugs, peptidomimetic, and non-peptidomimetic synthetic compounds, as well as natural compounds and their derivatives) are summarized comprehensively. Meanwhile, the challenges in this field are highlighted, while future directions for designing and developing efficacious 3CLpro inhibitors as novel anti-coronavirus therapies are also proposed. Collectively, all information and knowledge presented here are very helpful for understanding the structural features and inhibitory mechanisms of SARS-CoV-2 3CLpro inhibitors, which offers new insights or inspiration to medicinal chemists for designing and developing more efficacious 3CLpro inhibitors as novel anti-coronavirus agents.

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

The current pandemic caused by the COVID-19 disease has reached everywhere in the world and has affected every aspect of our lives. As of the current data, the World Health Organization (WHO) has reported more than 300 million confirmed COVID-19 cases worldwide and more than 5 million deaths. M^pro is an enzyme that plays a key role in the life cycle of the SARS-CoV-2 virus, and it is vital for the disease progression. The M^pro enzyme seems to have several allosteric sites that can hinder the enzyme catalytic activity. Furthermore, some of these allosteric sites are located at or nearby the dimerization interface which is essential for the overall M^pro activity. In this review paper, we investigate the potential of the M^pro allosteric site to act as a drug target, especially since they interestingly appear to be resistant to mutation. The work is illustrated through three subsequent sections: First, the two main categories of M^pro allosteric sites have been explained and discussed. Second, a total of six pockets have been studied and evaluated for their druggability and cavity characteristics. Third, the experimental and computational attempts for the discovery of new allosteric inhibitors have been illustrated and discussed. To sum up, this review paper gives a detailed insight into the feasibility of developing new M^pro inhibitors to act as a potential treatment for the COVID-19 disease.

A ricin-based peptide BRIP from Hordeum vulgare inhibits Mpro of SARS-CoV-2

COVID-19 pandemic caused by SARS-CoV-2 led to the research aiming to find the inhibitors of this virus. Towards this world problem, an attempt was made to identify SARS-CoV-2 main protease (Mpro) inhibitory peptides from ricin domains. The ricin-based peptide from barley (BRIP) was able to inhibit Mpro in vitro with an IC50 of 0.52 nM. Its low and no cytotoxicity upto 50 µM suggested its therapeutic potential against SARS-CoV-2. The most favorable binding site on Mpro was identified by molecular docking and steered molecular dynamics (MD) simulations. The Mpro-BRIP interactions were further investigated by evaluating the trajectories for microsecond timescale MD simulations. The structural parameters of Mpro-BRIP complex were stable, and the presence of oppositely charged surfaces on the binding interface of BRIP and Mpro complex further contributed to the overall stability of the protein-peptide complex. Among the components of thermodynamic binding free energy, Van der Waals and electrostatic contributions were most favorable for complex formation. Our findings provide novel insight into the area of inhibitor development against COVID-19.

The Functional Landscape of SARS-CoV-2 3CL Protease

SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) as the etiologic agent of COVID-19 (coronavirus disease 2019) has drastically altered life globally. Numerous efforts have been placed on the development of therapeutics to treat SARS-CoV-2 infection. One particular target is the 3CL protease (3CL ^pro ), which holds promise as it is essential to the virus and highly conserved among coronaviruses, suggesting that it may be possible to find broad inhibitors that treat not just SARS-CoV-2 but other coronavirus infections as well. While the 3CL protease has been studied by many groups for SARS-CoV-2 and other coronaviruses, our understanding of its tolerance to mutations is limited, knowledge which is particularly important as 3CL protease inhibitors become utilized clinically. Here, we develop a yeast-based deep mutational scanning approach to systematically profile the activity of all possible single mutants of the SARS-CoV-2 3CL ^pro , and validate our results both in yeast and in authentic viruses. We reveal that the 3CL ^pro is highly malleable and is capable of tolerating mutations throughout the protein, including within the substrate binding pocket. Yet, we also identify specific residues that appear immutable for function of the protease, suggesting that these interactions may be novel targets for the design of future 3CL ^pro inhibitors. Finally, we utilize our screening results as a basis to identify E166V as a resistance-conferring mutation against the therapeutic 3CL ^pro inhibitor, nirmatrelvir, in clinical use. Collectively, the functional map presented herein may serve as a guide for further understanding of the biological properties of the 3CL protease and for drug development for current and future coronavirus pandemics.

Allosteric Modulation of the Main Protease (MPro) of SARS-CoV-2 by Casticin—Insights from Molecular Dynamics Simulations

Inhibition of the main protease (M^pro) of SARS-CoV-2 has been suggested to be vital in shutting down viral replication in a host. Most efforts aimed at inhibiting M^Pro activity have been channeled into competitive inhibition at the active site, but this strategy will require a high inhibitor concentration and impressive inhibitor-M^Pro binding affinity. Allosteric inhibition can potentially serve as an effective strategy for alleviating these limitations. In this study, the ability of antiviral natural products to inhibit M^Pro in an allosteric fashion was explored with in silico techniques. Molecular docking revealed a strong interaction between casticin, an antiviral flavonoid, and M^pro at a site distant from the active site. This site, characterized as a distal site, has been shown to have an interdependent dynamic effect with the active site region. M^pro only, M^pro-peptide (binary) and M^pro-peptide-casticin (ternary) complexes were subjected to molecular dynamics simulations for 50 ns to investigate the modulatory activity of casticin binding on M^pro. Molecular dynamic simulations revealed that binding of casticin at the distal site interferes with the proper orientation of the peptide substrate in the oxyanion hole of the active site, and this could lead to a halt or decrease in catalytic activity. This study therefore highlights casticin as a potential allosteric modulator of the SARS-CoV-2 main protease, which could be optimized and developed into a potential lead compound for anti-SARS-CoV-2 chemotherapy.

Maturation of the SARS-CoV-2 virus is regulated by dimerization of its main protease

SARS-CoV-2 main protease (M^pro) involved in COVID-19 is required for maturation of the virus and infection of host cells. The key question is how to block the activity of M^pro. By combining atomistic simulations with machine learning, we found that the enzyme regulates its own activity by a collective allosteric mechanism that involves dimerization and binding of a single substrate. At the core of the collective mechanism is the coupling between the catalytic site residues, H41 and C145, which direct the activity of M^pro dimer, and two salt bridges formed between R4 and E290 at the dimer interface. If these salt bridges are mutated, the activity of M^pro is blocked. The results suggest that dimerization of main proteases is a general mechanism to foster coronavirus proliferation, and propose a robust drug-based strategy that does not depend on the frequently mutating spike proteins at the viral envelope used to develop vaccines.

Key dimer interface residues impact the catalytic activity of 3CLpro, the main protease of SARS-CoV-2

3C-like protease (3CLpro) is one of two proteases that process and liberate functional viral proteins essential for the maturation and infectivity of severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2), the virus responsible for COVID-19. It has been suggested that 3CLpro is catalytically active as a dimer, making the dimerization interface a target for antiviral development. Guided by structural analysis, here we introduced single amino acid substitutions at nine residues at three key sites of the dimer interface to assess their impact on dimerization and activity. We show that at site 1, alanine substitution of S1 or E166 increased by 2-fold or reduced relative activity, respectively. At site 2, alanine substitution of S10 or E14 eliminated activity, whereas K12A exhibited ∼60% relative activity. At site 3, alanine substitution of R4, E290, or Q299 eliminated activity, whereas S139A exhibited 46% relative activity. We further found the oligomerization states of the dimer interface mutants varied; the inactive mutants R4A, R4Q, S10A/C, E14A/D/Q/S, E290A, and Q299A/E were present as dimers, demonstrating that dimerization is not an indication of catalytically active 3CLpro. In addition, present mostly as monomers, K12A displayed residual activity, which could be attributed to the conspicuous amount of dimer present. Finally, differential scanning calorimetry did not reveal a direct relationship between the thermodynamic stability of mutants with oligomerization or catalytic activity. These results provide insights on two allosteric sites, R4/E290 and S10/E14, that may promote the design of antiviral compounds that target the dimer interface rather than the active site of SARS-CoV-2 3CLpro.

Targeting two potential sites of SARS-CoV-2 main protease through computational drug repurposing

Before the rise of SARS-CoV-2, emergence of different coronaviruses such as SARS-CoV and MERS-CoV has been reported that indicates possibility of the future novel pathogen from the coronavirus family at a pandemic level. In this context, explicit studies on identifying inhibitors focused on the coronavirus life cycle, are immensely important. The main protease is critical for the life cycle of coronaviruses. Majority of the work done on the inhibitor studies on the catalytically active dimeric SARS-CoV-2 main protease (M^pro), primarily focussed on the catalytic site of a single protomer, with a few targeting the dimeric site. In this study, we have exploited the FDA-approved drugs, for a computational drug repurposing study against the M^pro. A virtual screening approach was employed with docking and molecular dynamics (MD) methods. Out of 1576, FDA-approved compounds, our study suggests three compounds: netupitant, paliperidone and vilazodone as possible inhibitors with a potential to inhibit both sites (monomeric and dimeric) of the M^pro. These compounds were found to be stable during the MD simulations and their post simulation binding energies were also correlated for both the targeted sites, suggesting equal binding capacity. This unique efficiency of the reported compounds might support further experimental studies on developing inhibitors against SARS-CoV-2 main protease.Communicated by Ramaswamy H. Sarma.

Modulation of the monomer-dimer equilibrium and catalytic activity of SARS-CoV-2 main protease by a transition-state analog inhibitor

The role of dimer formation for the onset of catalytic activity of SARS-CoV-2 main protease (MPro^WT) was assessed using a predominantly monomeric mutant (MPro^M). Rates of MPro^WT and MPro^M catalyzed hydrolyses display substrate saturation kinetics and second-order dependency on the protein concentration. The addition of the prodrug GC376, an inhibitor of MPro^WT, to MPro^M leads to an increase in the dimer population and catalytic activity with increasing inhibitor concentration. The activity reaches a maximum corresponding to a dimer population in which one active site is occupied by the inhibitor and the other is available for catalytic activity. This phase is followed by a decrease in catalytic activity due to the inhibitor competing with the substrate. Detailed kinetics and equilibrium analyses are presented and a modified Michaelis-Menten equation accounts for the results. These observations provide conclusive evidence that dimer formation is coupled to catalytic activity represented by two equivalent active sites.

A new inactive conformation of SARS-CoV-2 main protease

The SARS‐CoV‐2 main protease (M^pro) has a pivotal role in mediating viral genome replication and transcription of the coronavirus, making it a promising target for drugs against the COVID‐19 pandemic. Here, a crystal structure is presented in which M^pro adopts an inactive state that has never been observed before, called new‐inactive. It is shown that the oxyanion loop, which is involved in substrate recognition and enzymatic activity, adopts a new catalytically incompetent conformation and that many of the key interactions of the active conformation of the enzyme around the active site are lost. Solvation/desolvation energetic contributions play an important role in the transition from the inactive to the active state, with Phe140 moving from an exposed to a buried environment and Asn142 moving from a buried environment to an exposed environment. In new‐inactive M^pro a new cavity is present near the S2′ subsite, and the N‐terminal and C‐terminal tails, as well as the dimeric interface, are perturbed, with partial destabilization of the dimeric assembly. This novel conformation is relevant both for comprehension of the mechanism of action of M^pro within the catalytic cycle and for the successful structure‐based drug design of antiviral drugs.

Allosteric Regulation of 3CL Protease of SARS-CoV-2 and SARS-CoV Observed in the Crystal Structure Ensemble

The 3C-like protease (3CL^pro) of SARS-CoV-2 is a potential therapeutic target for COVID-19. Importantly, it has an abundance of structural information solved as a complex with various drug candidate compounds. Collecting these crystal structures (83 Protein Data Bank (PDB) entries) together with those of the highly homologous 3CL^pro of SARS-CoV (101 PDB entries), we constructed the crystal structure ensemble of 3CL^pro to analyze the dynamic regulation of its catalytic function. The structural dynamics of the 3CL^pro dimer observed in the ensemble were characterized by the motions of four separate loops (the C-loop, E-loop, H-loop, and Linker) and the C-terminal domain III on the rigid core of the chymotrypsin fold. Among the four moving loops, the C-loop (also known as the oxyanion binding loop) causes the order (active)-disorder (collapsed) transition, which is regulated cooperatively by five hydrogen bonds made with the surrounding residues. The C-loop, E-loop, and Linker constitute the major ligand binding sites, which consist of a limited variety of binding residues including the substrate binding subsites. Ligand binding causes a ligand size dependent conformational change to the E-loop and Linker, which further stabilize the C-loop via the hydrogen bond between the C-loop and E-loop. The T285A mutation from SARS-CoV 3CL^pro to SARS-CoV-2 3CL^pro significantly closes the interface of the domain III dimer and allosterically stabilizes the active conformation of the C-loop via hydrogen bonds with Ser1 and Gly2; thus, SARS-CoV-2 3CL^pro seems to have increased activity relative to that of SARS-CoV 3CL^pro.

Pharmacophore Model for SARS-CoV-2 3CLpro Small-Molecule Inhibitors and in Vitro Experimental Validation of Computationally Screened Inhibitors

Among the biomedical efforts in response to the current coronavirus (COVID-19) pandemic, pharmacological strategies to reduce viral load in patients with severe forms of the disease are being studied intensively. One of the main drug target proteins proposed so far is the SARS-CoV-2 viral protease 3CLpro (also called Mpro), an essential component for viral replication. Ongoing ligand- and receptor-based computational screening efforts would be facilitated by an improved understanding of the electrostatic, hydrophobic, and steric features that characterize small-molecule inhibitors binding stably to 3CLpro and by an extended collection of known binders. Here, we present combined virtual screening, molecular dynamics (MD) simulation, machine learning, and in vitro experimental validation analyses, which have led to the identification of small-molecule inhibitors of 3CLpro with micromolar activity and to a pharmacophore model that describes functional chemical groups associated with the molecular recognition of ligands by the 3CLpro binding pocket. Experimentally validated inhibitors using a ligand activity assay include natural compounds with the available prior knowledge on safety and bioavailability properties, such as the natural compound rottlerin (IC50 = 37 μM) and synthetic compounds previously not characterized (e.g., compound CID 46897844, IC50 = 31 μM). In combination with the developed pharmacophore model, these and other confirmed 3CLpro inhibitors may provide a basis for further similarity-based screening in independent compound databases and structural design optimization efforts to identify 3CLpro ligands with improved potency and selectivity. Overall, this study suggests that the integration of virtual screening, MD simulations, and machine learning can facilitate 3CLpro-targeted small-molecule screening investigations. Different receptor-, ligand-, and machine learning-based screening strategies provided complementary information, helping to increase the number and diversity of the identified active compounds. Finally, the resulting pharmacophore model and experimentally validated small-molecule inhibitors for 3CLpro provide resources to support follow-up computational screening efforts for this drug target.

Regulation of the Dimerization and Activity of SARS-CoV-2 Main Protease through Reversible Glutathionylation of Cysteine 300

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent for coronavirus disease 2019 (COVID-19), encodes two proteases required for replication. The main protease (M^pro), encoded as part of two polyproteins, pp1a and pp1ab, is responsible for 11 different cleavages of these viral polyproteins to produce mature proteins required for viral replication. M^pro is therefore an attractive target for therapeutic interventions. Certain proteins in cells under oxidative stress undergo modification of reactive cysteines. We show M^pro is susceptible to glutathionylation, leading to inhibition of dimerization and activity. Activity of glutathionylated M^pro could be restored with reducing agents or glutaredoxin. Analytical studies demonstrated that glutathionylated M^pro primarily exists as a monomer and that modification of a single cysteine with glutathione is sufficient to block dimerization and inhibit its activity. Gel filtration studies as well as analytical ultracentrifugation confirmed that glutathionylated M^pro exists as a monomer. Tryptic and chymotryptic digestions of M^pro as well as experiments using a C300S M^pro mutant revealed that Cys300, which is located at the dimer interface, is a primary target of glutathionylation. Moreover, Cys300 is required for inhibition of activity upon M^pro glutathionylation. These findings indicate that M^pro dimerization and activity can be regulated through reversible glutathionylation of a non-active site cysteine, Cys300, which itself is not required for M^pro activity, and provides a novel target for the development of agents to block M^pro dimerization and activity. This feature of M^pro may have relevance to the pathophysiology of SARS-CoV-2 and related bat coronaviruses.

Interfacial Water Many-Body Effects Drive Structural Dynamics and Allosteric Interactions in SARS-CoV-2 Main Protease Dimerization Interface

Following our previous work ( Chem. Sci. 2021, 12, 4889-4907), we study the structural dynamics of the SARS-CoV-2 Main Protease dimerization interface (apo dimer) by means of microsecond adaptive sampling molecular dynamics simulations (50 μs) using the AMOEBA polarizable force field (PFF). This interface is structured by a complex H-bond network that is stable only at physiological pH. Structural correlations analysis between its residues and the catalytic site confirms the presence of a buried allosteric site. However, noticeable differences in allosteric connectivity are observed between PFFs and non-PFFs. Interfacial polarizable water molecules are shown to appear at the heart of this discrepancy because they are connected to the global interface H-bond network and able to adapt their dipole moment (and dynamics) to their diverse local physicochemical microenvironments. The water-interface many-body interactions appear to drive the interface volume fluctuations and to therefore mediate the allosteric interactions with the catalytic cavity.

De novo ssRNA Aptamers against the SARS-CoV-2 Main Protease: In Silico Design and Molecular Dynamics Simulation

Herein, we have generated ssRNA aptamers to inhibit SARS-CoV-2 M^pro, a protease necessary for the SARS-CoV-2 coronavirus replication. Because there is no aptamer 3D structure currently available in the databanks for this protein, first, we modeled an ssRNA aptamer using an entropic fragment-based strategy. We refined the initial sequence and 3D structure by using two sequential approaches, consisting of an elitist genetic algorithm and an RNA inverse process. We identified three specific aptamers against SARS-CoV-2 M^pro, called MApta^pro, MApta^pro-IR1, and MApta^pro-IR2, with similar 3D conformations and that fall in the dimerization region of the SARS-CoV-2 M^pro necessary for the enzymatic activity. Through the molecular dynamic simulation and binding free energy calculation, the interaction between the MApta^pro-IR1 aptamer and the SARS-CoV-2 M^pro enzyme resulted in the strongest and the highest stable complex; therefore, the ssRNA MApta^pro-IR1 aptamer was selected as the best potential candidate for the inhibition of SARS-CoV-2 M^pro and a perspective therapeutic drug for the COVID-19 disease.

Catalytic Dyad Residues His41 and Cys145 Impact the Catalytic Activity and Overall Conformational Fold of the Main SARS-CoV-2 Protease 3-Chymotrypsin-Like Protease

Coronaviruses are responsible for multiple pandemics and millions of deaths globally, including the current pandemic of coronavirus disease 2019 (COVID-19). Development of antivirals against coronaviruses, including the severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2) responsible for COVID-19, is essential for containing the current and future coronavirus outbreaks. SARS-CoV-2 proteases represent important targets for the development of antivirals because of their role in the processing of viral polyproteins. 3-Chymotrypsin-like protease (3CLpro) is one such protease. The cleavage of SARS-CoV-2 polyproteins by 3CLpro is facilitated by a Cys145–His41 catalytic dyad. We here characterized the catalytic roles of the cysteine–histidine pair for improved understanding of the 3CLpro reaction mechanism, to inform the development of more effective antivirals against Sars-CoV-2. The catalytic dyad residues were substituted by site-directed mutagenesis. All substitutions tested (H41A, H41D, H41E, C145A, and C145S) resulted in a complete inactivation of 3CLpro, even when amino acids with a similar catalytic function to that of the original residues were used. The integrity of the structural fold of enzyme variants was investigated by circular dichroism spectroscopy to test if the catalytic inactivation of 3CLpro was caused by gross changes in the enzyme secondary structure. C145A, but not the other substitutions, shifted the oligomeric state of the enzyme from dimeric to a higher oligomeric state. Finally, the thermodynamic stability of 3CLpro H41A, H41D, and C145S variants was reduced relative the wild-type enzyme, with a similar stability of the H41E and C145A variants. Collectively, the above observations confirm the roles of His41 and Cys145 in the catalytic activity and the overall conformational fold of 3CLpro SARS-CoV-2. We conclude that the cysteine–histidine pair should be targeted for inhibition of 3CLpro and development of antiviral against COVID-19 and coronaviruses.

Proposition of a new allosteric binding site for potential SARS-CoV-2 3CL protease inhibitors by utilizing molecular dynamics simulations and ensemble docking

The SARS-CoV-2 3CL protease (3CLpro) shows a high similarity with 3CL proteases of other beta-coronaviruses, such as SARS and MERS. It is the main enzyme involved in generating various non-structural proteins that are important for viral replication and is one of the most important proteins responsible for SARS-CoV-2 virulence. In this study, we have conducted an ensemble docking of molecules from the DrugBank database using both the crystallographic structure of the SARS-CoV-2 3CLpro, as well as five conformations obtained after performing a cluster analysis of a 300 ns molecular dynamics (MD) simulation. This procedure elucidated the inappropriateness of the active site for non-covalent inhibitors, but it has also shown that there exists an additional, more favorable, allosteric binding site, which could be a better target for non-covalent inhibitors, as it could prevent dimerization and activation of SARS-CoV-2 3CLpro. Two such examples are radotinib and nilotinib, tyrosine kinase inhibitors already in use for treatment of leukemia and which binding to the newly found allosteric binding site was also confirmed using MD simulations. Communicated by Ramaswamy H. Sarma.

Targeting SARS-CoV-2 viral proteases as a therapeutic strategy to treat COVID-19

The 21st century has witnessed three outbreaks of coronavirus (CoVs) infections caused by severe acute respiratory syndrome (SARS)-CoV, Middle East respiratory syndrome (MERS)-CoV, and SARS-CoV-2. Coronavirus disease 2019 (COVID-19), caused by SARS-CoV-2, spreads rapidly and since the discovery of the first COVID-19 infection in December 2019, has caused 1.2 million deaths worldwide and 226,777 deaths in the United States alone. The high amino acid similarity between SARS-CoV and SARS-CoV-2 viral proteins supports testing therapeutic molecules that were designed to treat SARS infections during the 2003 epidemic. In this review, we provide information on possible COVID-19 treatment strategies that act via inhibition of the two essential proteins of the virus, 3C-like protease (3CLpro ) or papain-like protease (PLpro ).

Targeting protein-protein interaction interfaces in COVID-19 drug discovery

To date, the COVID-19 pandemic has claimed over 1 million human lives, infected another 50 million individuals and wreaked havoc on the global economy. The crisis has spurred the ongoing development of drugs targeting its etiological agent, the SARS-CoV-2. Targeting relevant protein-protein interaction interfaces (PPIIs) is a viable paradigm for the design of antiviral drugs and enriches the targetable chemical space by providing alternative targets for drug discovery. In this review, we will provide a comprehensive overview of the theory, methods and applications of PPII-targeted drug development towards COVID-19 based on recent literature. We will also highlight novel developments, such as the successful use of non-native protein-protein interactions as targets for antiviral drug screening. We hope that this review may serve as an entry point for those interested in applying PPIIs towards COVID-19 drug discovery and speed up drug development against the pandemic.

A review of the latest research on Mpro targeting SARS-COV inhibitors

Since the outbreak of COVID-19, the pandemic caused by SARS-CoV-2 infection is still spreading at an alarming rate and has caused huge loss of life and economic damage worldwide. Although more than one year has passed, effective treatments for COVID-19 and other pathogenic coronaviruses have not yet been developed. Therefore, the development of SARS-CoV-2 inhibitors is an urgent priority. Given that the M^pro sequences of SARS-CoV-2 and SARS-CoV-1 are 100% identical in the catalytic domain for protein cleavage, the viral main protease (M^pro) is one of the most extensive drug targets in all the drug targets being investigated for SARS-CoV-2. To provide scientific researchers with timely anti-SARS-CoV drug development information for M^pro, we focus on the past and current drug design and development strategies for M^Pro in this review. We believe that this review will provide meaningful guidance for the design and development of innovative drugs against COVID-19 and other pathogenic coronaviruses in the future.

Regulation of the Dimerization and Activity of SARS-CoV-2 Main Protease through Reversible Glutathionylation of Cysteine 300

SARS-CoV-2 encodes main protease (M^pro), an attractive target for therapeutic interventions. We show M^pro is susceptible to glutathionylation leading to inhibition of dimerization and activity. Activity of glutathionylated M^pro could be restored with reducing agents or glutaredoxin. Analytical studies demonstrated that glutathionylated M^pro primarily exists as a monomer and that a single modification with glutathione is sufficient to block dimerization and loss of activity. Proteolytic digestions of M^pro revealed Cys300 as a primary target of glutathionylation, and experiments using a C300S M^pro mutant revealed that Cys300 is required for inhibition of activity upon M^pro glutathionylation. These findings indicate that M^pro dimerization and activity can be regulated through reversible glutathionylation of Cys300 and provides a novel target for the development of agents to block Mpro dimerization and activity. This feature of M^pro may have relevance to human disease and the pathophysiology of SARS-CoV-2 in bats, which develop oxidative stress during flight.

Targeting novel structural and functional features of coronavirus protease nsp5 (3CLpro, Mpro) in the age of COVID-19

Coronavirus protease nsp5 (Mpro, 3CLpro) remains a primary target for coronavirus therapeutics due to its indispensable and conserved role in the proteolytic processing of the viral replicase polyproteins. In this review, we discuss the diversity of known coronaviruses, the role of nsp5 in coronavirus biology, and the structure and function of this protease across the diversity of known coronaviruses, and evaluate past and present efforts to develop inhibitors to the nsp5 protease with a particular emphasis on new and mostly unexplored potential targets of inhibition. With the recent emergence of pandemic SARS-CoV-2, this review provides novel and potentially innovative strategies and directions to develop effective therapeutics against the coronavirus protease nsp5.

An update review of emerging small-molecule therapeutic options for COVID-19

The SARS-CoV-2 outbreak and pandemic that began near the end of 2019 has posed a challenge to global health. At present, many candidate small-molecule therapeutics have been developed that can inhibit both the infection and replication of SARS-CoV-2 and even potentially relieve cytokine storms and other related complications. Meanwhile, host-targeted drugs that inhibit cellular transmembrane serine protease (TMPRSS2) can prevent SARS-CoV-2 from entering cells, and its combination with chloroquine and dihydroorotate dehydrogenase (DHODH) inhibitors can limit the spread of SARS-CoV-2 and reduce the morbidity and mortality of patients with COVID-19. The present article provides an overview of these small-molecule therapeutics based on insights from medicinal chemistry research and focuses on RNA-dependent RNA polymerase (RdRp) inhibitors, such as the nucleoside analogues remdesivir, favipiravir and ribavirin. This review also covers inhibitors of 3C-like protease (3CL^pro), papain-like protease (PL^pro) and other potentially innovative active ingredient molecules, describing their potential targets, activities, clinical status and side effects.

Impact of dimerization and N3 binding on molecular dynamics of SARS-CoV and SARS-CoV-2 main proteases

SARS-CoV-2 main protease is one of the major targets in drug development efforts against Covid-19. Even though several structures were reported to date, its dynamics is not understood well. In particular, impact of dimerization and ligand binding on the dynamics is an important issue to investigate. In this study, we performed molecular dynamics simulations of SARS-CoV and SARS-CoV-2 main proteases to investigate influence of dimerization on the dynamics by modeling monomeric and dimeric apo and holo forms. The dimerization causes an organization of the interdomain dynamics as well as some local structural changes. Moreover, we investigated impact of a peptide mimetic (N3) on the dynamics of SARS-CoV and SARS-CoV-2 Mpro. The ligand binding to the dimeric forms causes some key local changes at the dimer interface and it causes an allosteric interaction between the active sites of two protomers. Our results support the idea that only one protomer is active on SARS-CoV-2 due to this allosteric interaction. Additionally, we analyzed the molecular dynamics trajectories from graph theoretical perspective and found that the most influential residues – as measured by eigenvector centrality – are a group of residues in active site and dimeric interface of the protease. This study may form a bridge between what we know about the dynamics of SARS-CoV and SARS-CoV-2 Mpro. We think that enlightening allosteric communication of the active sites and the role of dimerization in SARS-CoV-2 Mpro can contribute to development of novel drugs against this global health problem as well as other similar proteases.

Communicated by Ramaswamy H. Sarma

Overview of antiviral drug candidates targeting coronaviral 3C-like main proteases

Coronaviruses (CoVs) are positive single-stranded RNA viruses that cause severe respiratory syndromes in humans, including severe acute respiratory syndrome (SARS) and Middle East respiratory syndrome (MERS). Coronavirus disease 2019 (COVID-19) caused by a novel severe acute respiratory syndrome CoV (SARS-CoV-2) at the end of 2019 became a global pandemic. The 3C-like cysteine protease (3CLpro) processes viral polyproteins to yield mature non-structural proteins, thus playing an important role in the CoV life cycle, and therefore is considered as a prominent target for antiviral drugs. To date, many 3CLpro inhibitors have been reported, and their molecular mechanisms have been illustrated. Here, we briefly introduce the structural features of 3CLpro of the human-related SARS-CoV, MERS-CoV and SARS-CoV-2, and explore the potency and mechanism of their cognate inhibitors. This information will shed light on the development and optimization of CoV 3CLpro inhibitors, which may benefit the further designation of therapeutic strategies for treating CoV diseases.

What coronavirus 3C-like protease tells us: From structure, substrate selectivity, to inhibitor design

The emergence of a variety of coronaviruses (CoVs) in the last decades has posed huge threats to human health. Especially, the ongoing pandemic of coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has led to more than 70 million infections and over 1.6 million of deaths worldwide in the past few months. None of the efficacious antiviral agents against human CoVs have been approved yet. 3C-like protease (3CL^pro ) is an attractive target for antiviral intervention due to its essential role in processing polyproteins translated from viral RNA, and its conserved structural feature and substrate specificity among CoVs in spite of the sequence variation. This review focuses on all available crystal structures of 12 CoV 3CL^pro s and their inhibitors, and intends to provide a comprehensive understanding of this protease from multiple aspects including its structural features, substrate specificity, inhibitor binding modes, and more importantly, to recapitulate the similarity and diversity among different CoV 3CL^pro s and the structure-activity relationship of various types of inhibitors. Such an attempt could gain a deep insight into the inhibition mechanisms and drive future structure-based drug discovery targeting 3CL^pro s.

Unusual zwitterionic catalytic site of SARS-CoV-2 main protease revealed by neutron crystallography

The main protease (3CL Mpro) from SARS-CoV-2, the etiological agent of COVID-19, is an essential enzyme for viral replication. 3CL Mpro possesses an unusual catalytic dyad composed of Cys145 and His41 residues. A critical question in the field has been what the protonation states of the ionizable residues in the substrate-binding active-site cavity are; resolving this point would help understand the catalytic details of the enzyme and inform rational drug development against this pernicious virus. Here, we present the room-temperature neutron structure of 3CL Mpro, which allowed direct determination of hydrogen atom positions and, hence, protonation states in the protease. We observe that the catalytic site natively adopts a zwitterionic reactive form in which Cys145 is in the negatively charged thiolate state and His41 is doubly protonated and positively charged, instead of the neutral unreactive state usually envisaged. The neutron structure also identified the protonation states, and thus electrical charges, of all other amino acid residues and revealed intricate hydrogen-bonding networks in the active-site cavity and at the dimer interface. The fine atomic details present in this structure were made possible by the unique scattering properties of the neutron, which is an ideal probe for locating hydrogen positions and experimentally determining protonation states at near-physiological temperature. Our observations provide critical information for structure-assisted and computational drug design, allowing precise tailoring of inhibitors to the enzyme's electrostatic environment.

The development of Coronavirus 3C-Like protease (3CLpro) inhibitors from 2010 to 2020

This review fully describes the coronavirus 3CL^pro peptidomimetic inhibitors and nonpeptidic small molecule inhibitors developed from 2010 to 2020. Specifically, the structural characteristics, binding modes and SARs of these 3CL^pro inhibitors are expounded in detail by division into two categories: peptidomimetic inhibitors mainly utilize electrophilic warhead groups to covalently bind the 3CL^pro Cys145 residue and thereby achieve irreversible inhibition effects, whereas nonpeptidic small molecule inhibitors mainly interact with residues in the S1', S1, S2 and S4 pockets via hydrogen bonds, hydrophobic bonds and van der Waals forces. Based on the emerging PROTAC technology and the existing 3CL^pro inhibitors, 3CL^pro PROTAC degraders are hypothesised to be next-generation anti-coronavirus drugs.

In silico pharmacokinetic and molecular docking studies of natural flavonoids and synthetic indole chalcones against essential proteins of SARS-CoV-2

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is distinctly infective and there is an ongoing effort to find a cure for this pandemic. Flavonoids exist in many diets as well as in traditional medicine, and their modern subset, indole-chalcones, are effective in fighting various diseases. Hence, these flavonoids and structurally similar indole chalcones derivatives were studied in silico for their pharmacokinetic properties including absorption, distribution, metabolism, excretion, toxicity (ADMET) and anti-SARS-CoV-2 properties against their proteins, namely, RNA dependent RNA polymerase (rdrp), main protease (M^pro) and Spike (S) protein via homology modelling and docking. Interactions were studied with respect to biology and function of SARS-CoV-2 proteins for activity. Functional/structural roles of amino acid residues of SARS-CoV-2 proteins and, the effect of flavonoid and indole chalcone interactions which may cause disease suppression are discussed. The results reveal that out of 23 natural flavonoids and 25 synthetic indole chalcones, 30 compounds are capable of M^pro deactivation as well as potentially lowering the efficiency of M^pro function. Cyanidin may inhibit RNA polymerase function and, Quercetin is found to block interaction sites on the viral spike. These results suggest flavonoids and their modern pharmaceutical cousins, indole chalcones are capable of fighting SARS-CoV-2. The in vitro anti-SARS-CoV-2 activity of these 30 compounds needs to be studied further for complete understanding and confirmation of their inhibitory potential.

Impact of Early Pandemic Stage Mutations on Molecular Dynamics of SARS-CoV-2 Mpro

A new coronavirus (SARS-CoV-2) is a global threat to world health and economy. Its dimeric main protease (Mpro), which is required for the proteolytic cleavage of viral precursor proteins, is a good candidate for drug development owing to its conservation and the absence of a human homolog. Improving our understanding of Mpro behavior can accelerate the discovery of effective therapies to reduce mortality. All-atom molecular dynamics (MD) simulations (100 ns) of 50 mutant Mpro dimers obtained from filtered sequences from the GISAID database were analyzed using root-mean-square deviation, root-mean-square fluctuation, Rg, averaged betweenness centrality, and geometry calculations. The results showed that SARS-CoV-2 Mpro essentially behaves in a similar manner to its SAR-CoV homolog. However, we report the following new findings from the variants: (1) Residues GLY15, VAL157, and PRO184 have mutated more than once in SARS CoV-2; (2) the D48E variant has lead to a novel "TSEEMLN"" loop at the binding pocket; (3) inactive apo Mpro does not show signs of dissociation in 100 ns MD; (4) a non-canonical pose for PHE140 widens the substrate binding surface; (5) dual allosteric pockets coinciding with various stabilizing and functional components of the substrate binding pocket were found to display correlated compaction dynamics; (6) high betweenness centrality values for residues 17 and 128 in all Mpro samples suggest their high importance in dimer stability-one such consequence has been observed for the M17I mutation whereby one of the N-fingers was highly unstable. (7) Independent coarse-grained Monte Carlo simulations suggest a relationship between the rigidity/mutability and enzymatic function. Our entire approach combining database preparation, variant retrieval, homology modeling, dynamic residue network (DRN), relevant conformation retrieval from 1-D kernel density estimates from reaction coordinates to other existing approaches of structural analysis, and data visualization within the coronaviral Mpro is also novel and is applicable to other coronaviral proteins.

The SARS-CoV-2 main protease as drug target

The unprecedented pandemic of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is threatening global health. The virus emerged in late 2019 and can cause a severe disease associated with significant mortality. Several vaccine development and drug discovery campaigns are underway. The SARS-CoV-2 main protease is considered a promising drug target, as it is dissimilar to human proteases. Sequence and structure of the main protease are closely related to those from other betacoronaviruses, facilitating drug discovery attempts based on previous lead compounds. Covalently binding peptidomimetics and small molecules are investigated. Various compounds show antiviral activity in infected human cells.

Progress in Developing Inhibitors of SARS-CoV-2 3C-Like Protease

Coronavirus disease 2019 (COVID-19) is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The viral outbreak started in late 2019 and rapidly became a serious health threat to the global population. COVID-19 was declared a pandemic by the World Health Organization in March 2020. Several therapeutic options have been adopted to prevent the spread of the virus. Although vaccines have been developed, antivirals are still needed to combat the infection of this virus. SARS-CoV-2 is an enveloped virus, and its genome encodes polyproteins that can be processed into structural and nonstructural proteins. Maturation of viral proteins requires cleavages by proteases. Therefore, the main protease (3 chymotrypsin-like protease (3CLpro) or Mpro) encoded by the viral genome is an attractive drug target because it plays an important role in cleaving viral polyproteins into functional proteins. Inhibiting this enzyme is an efficient strategy to block viral replication. Structural studies provide valuable insight into the function of this protease and structural basis for rational inhibitor design. In this review, we describe structural studies on the main protease of SARS-CoV-2. The strategies applied in developing inhibitors of the main protease of SARS-CoV-2 and currently available protein inhibitors are summarized. Due to the availability of high-resolution structures, structure-guided drug design will play an important role in developing antivirals. The availability of high-resolution structures, potent peptidic inhibitors, and diverse compound scaffolds indicate the feasibility of developing potent protease inhibitors as antivirals for COVID-19.

Targeting severe acute respiratory syndrome-coronavirus (SARS-CoV-1) with structurally diverse inhibitors: a comprehensive review

Coronaviruses, which were discovered in 1968, can lead to some human viral disorders, like severe acute respiratory syndrome (SARS), Middle East respiratory syndrome-related (MERS), and, recently, coronavirus disease 2019 (COVID-19). The coronavirus that leads to COVID-19 is rapidly spreading all over the world and is the reason for the deaths of thousands of people. Recent research has revealed that there is about 80% sequence homology between the coronaviruses that cause SARS and COVID-19. Considering this fact, we decided to collect the maximum available information on targets, structures, and inhibitors reported so far for SARS-CoV-1 that could be useful for researchers who work on closely related COVID-19. There are vital proteases, like papain-like protease 2 (PL2pro) and 3C-like protease (3CLpro), or main protease (Mpro), that are involved in and are essential for the replication of SARS coronavirus and so are valuable targets for the treatment of patients affected by this type of virus. SARS-CoV-1 NTPase/helicase plays an important role in the release of several non-structural proteins (nsps), so it is another essential target relating to the viral life cycle. In this paper, we provide extensive information about diverse molecules with anti-SARS activity. In addition to traditional medicinal chemistry outcomes, HTS, virtual screening efforts, and structural insights for better understanding inhibitors and SARS-CoV-1 target complexes are also discussed. This study covers a wide range of anti-SARS agents, particularly SARS-CoV-1 inhibitors, and provides new insights into drug design for the deadly SARS-CoV-2 virus.

Targeting the Dimerization of the Main Protease of Coronaviruses: A Potential Broad-Spectrum Therapeutic Strategy

A new coronavirus (CoV) caused a pandemic named COVID-19, which has become a global health care emergency in the present time. The virus is referred to as SARS-CoV-2 (severe acute respiratory syndrome-coronavirus-2) and has a genome similar (∼82%) to that of the previously known SARS-CoV (SARS coronavirus). An attractive therapeutic target for CoVs is the main protease (Mpro) or 3-chymotrypsin-like cysteine protease (3CLpro), as this enzyme plays a key role in polyprotein processing and is active in a dimeric form. Further, Mpro is highly conserved among various CoVs, and a mutation in Mpro is often lethal to the virus. Thus, drugs targeting the Mpro enzyme significantly reduce the risk of mutation-mediated drug resistance and display broad-spectrum antiviral activity. The combinatorial design of peptide-based inhibitors targeting the dimerization of SARS-CoV Mpro represents a potential therapeutic strategy. In this regard, we have compiled the literature reports highlighting the effect of mutations and N-terminal deletion of residues of SARS-CoV Mpro on its dimerization and, thus, catalytic activity. We believe that the present review will stimulate research in this less explored yet quite significant area. The effect of the COVID-19 epidemic and the possibility of future CoV outbreaks strongly emphasize the urgent need for the design and development of potent antiviral agents against CoV infections.

An Overview of Severe Acute Respiratory Syndrome-Coronavirus (SARS-CoV) 3CL Protease Inhibitors: Peptidomimetics and Small Molecule Chemotherapy

Severe acute respiratory syndrome (SARS) is caused by a newly emerged coronavirus that infected more than 8000 individuals and resulted in more than 800 (10-15%) fatalities in 2003. The causative agent of SARS has been identified as a novel human coronavirus (SARS-CoV), and its viral protease, SARS-CoV 3CL(pro), has been shown to be essential for replication and has hence been recognized as a potent drug target for SARS infection. Currently, there is no effective treatment for this epidemic despite the intensive research that has been undertaken since 2003 (over 3500 publications). This perspective focuses on the status of various efficacious anti-SARS-CoV 3CL(pro) chemotherapies discovered during the last 12 years (2003-2015) from all sources, including laboratory synthetic methods, natural products, and virtual screening. We describe here mainly peptidomimetic and small molecule inhibitors of SARS-CoV 3CL(pro). Attempts have been made to provide a complete description of the structural features and binding modes of these inhibitors under many conditions.

Structural basis for the dimerization and substrate recognition specificity of porcine epidemic diarrhea virus 3C-like protease

Porcine epidemic diarrhea virus (PEDV), a member of the genus Alphacoronavirus, has caused significant damage to the Asian and American pork industries. Coronavirus 3C-like protease (3CL^pro), which is involved in the processing of viral polyproteins for viral replication, is an appealing antiviral drug target. Here, we present the crystal structures of PEDV 3CL^pro and a molecular complex between an inactive PEDV 3CL^pro variant C144A bound to a peptide substrate. Structural characterization, mutagenesis and biochemical analysis reveal the substrate-binding pockets and the residues that comprise the active site of PEDV 3CL^pro. The dimerization of PEDV 3CL^pro is similar to that of other Alphacoronavirus 3CL^pros but has several differences from that of SARS-CoV 3CL^pro from the genus Betacoronavirus. Furthermore, the non-conserved motifs in the pockets cause different cleavage of substrate between PEDV and SARS-CoV 3CL^pros, which may provide new insights into the recognition of substrates by 3CL^pros in various coronavirus genera.

•

The substrate binding mechanism of PEDV 3CL^pro has been characterized.

•

The large buried surface area is responsible for dimerization of PEDV 3CL^pro.

•

Non-conserved motifs cause different cleavage between PEDV and SARS-CoV 3CL^pros.

Conformational Flexibility of a Short Loop near the Active Site of the SARS-3CLpro is Essential to Maintain Catalytic Activity

The SARS 3C-like proteinase (SARS-3CLpro), which is the main proteinase of the SARS coronavirus, is essential to the virus life cycle. This enzyme has been shown to be active as a dimer in which only one protomer is active. However, it remains unknown how the dimer structure maintains an active monomer conformation. It has been observed that the Ser139-Leu141 loop forms a short 3(10)-helix that disrupts the catalytic machinery in the inactive monomer structure. We have tried to disrupt this helical conformation by mutating L141 to T in the stable inactive monomer G11A/R298A/Q299A. The resulting tetra-mutant G11A/L141T/R298A/Q299A is indeed enzymatically active as a monomer. Molecular dynamics simulations revealed that the L141T mutation disrupts the 3(10)-helix and helps to stabilize the active conformation. The coil-3(10)-helix conformational transition of the Ser139-Leu141 loop serves as an enzyme activity switch. Our study therefore indicates that the dimer structure can stabilize the active conformation but is not a required structure in the evolution of the active enzyme, which can also arise through simple mutations.

Potential Broad Spectrum Inhibitors of the Coronavirus 3CLpro: A Virtual Screening and Structure-Based Drug Design Study

Human coronaviruses represent a significant disease burden; however, there is currently no antiviral strategy to combat infection. The outbreak of severe acute respiratory syndrome (SARS) in 2003 and Middle East respiratory syndrome (MERS) less than 10 years later demonstrates the potential of coronaviruses to cross species boundaries and further highlights the importance of identifying novel lead compounds with broad spectrum activity. The coronavirus 3CL(pro) provides a highly validated drug target and as there is a high degree of sequence homology and conservation in main chain architecture the design of broad spectrum inhibitors is viable. The ZINC drugs-now library was screened in a consensus high-throughput pharmacophore modeling and molecular docking approach by Vina, Glide, GOLD and MM-GBSA. Molecular dynamics further confirmed results obtained from structure-based techniques. A highly defined hit-list of 19 compounds was identified by the structure-based drug design methodologies. As these compounds were extensively validated by a consensus approach and by molecular dynamics, the likelihood that at least one of these compounds is bioactive is excellent. Additionally, the compounds segregate into 15 significantly dissimilar (p < 0.05) clusters based on shape and features, which represent valuable scaffolds that can be used as a basis for future anti-coronaviral inhibitor discovery experiments. Importantly though, the enriched subset of 19 compounds identified from the larger library has to be validated experimentally.

Dynamically-driven enhancement of the catalytic machinery of the SARS 3C-like protease by the S284-T285-I286/A mutations on the extra domain

Previously we revealed that the extra domain of SARS 3CLpro mediated the catalysis via different mechanisms. While the R298A mutation completely abolished the dimerization, thus resulting in the inactive catalytic machinery, N214A inactivated the enzyme by altering its dynamics without significantly perturbing its structure. Here we studied another mutant with S284-T285-I286 replaced by Ala (STI/A) with a 3.6-fold activity increase and slightly enhanced dimerization. We determined its crystal structure, which still adopts the dimeric structure almost identical to that of the wild-type (WT), except for slightly tighter packing between two extra-domains. We then conducted 100-ns molecular dynamics (MD) simulations for both STI/A and WT, the longest reported so far for 3CLpro. In the simulations, two STI/A extra domains become further tightly packed, leading to a significant volume reduction of the nano-channel formed by residues from both catalytic and extra domains. The enhanced packing appears to slightly increase the dynamic stability of the N-finger and the first helix residues, which subsequently triggers the redistribution of dynamics over residues directly contacting them. This ultimately enhances the dynamical stability of the residues constituting the catalytic dyad and substrate-binding pockets. Further correlation analysis reveals that a global network of the correlated motions exists in the protease, whose components include all residues identified so far to be critical for the dimerization and catalysis. Most strikingly, the N214A mutation globally decouples this network while the STI/A mutation alters the correlation pattern. Together with previous results, the present study establishes that besides the classic structural allostery, the dynamic allostery also operates in the SARS 3CLpro, which is surprisingly able to relay the perturbations on the extra domain onto the catalytic machinery to manifest opposite catalytic effects. Our results thus imply a promising avenue to design specific inhibitors for 3CL proteases by disrupting their dynamic correlation network.