Long-range cooperative interactions modulate dimerization in SARS 3CLpro

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.

An orally bioavailable SARS-CoV-2 main protease inhibitor exhibits improved affinity and reduced sensitivity to mutations

Inhibitors of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main protease (Mpro) such as nirmatrelvir (NTV) and ensitrelvir (ETV) have proven effective in reducing the severity of COVID-19, but the presence of resistance-conferring mutations in sequenced viral genomes raises concerns about future drug resistance. Second-generation oral drugs that retain function against these mutants are thus urgently needed. We hypothesized that the covalent hepatitis C virus protease inhibitor boceprevir (BPV) could serve as the basis for orally bioavailable drugs that inhibit SARS-CoV-2 Mpro more efficiently than existing drugs. Performing structure-guided modifications of BPV, we developed a picomolar-affinity inhibitor, ML2006a4, with antiviral activity, oral pharmacokinetics, and therapeutic efficacy similar or superior to those of NTV. A crucial feature of ML2006a4 is a derivatization of the ketoamide reactive group that improves cell permeability and oral bioavailability. Last, ML2006a4 was found to be less sensitive to several mutations that cause resistance to NTV or ETV and occur in the natural SARS-CoV-2 population. Thus, anticipatory design can preemptively address potential resistance mechanisms to expand future treatment options against coronavirus variants.

The SARS-CoV-2 Mpro Dimer-Based Screening System: A Synthetic Biology Tool for Identifying Compounds with Dimerization Inhibitory Potential

The COVID-19 endemic remains a global concern. The search for effective antiviral candidates is still needed to reduce disease risk. However, the availability of high biosafety level laboratory facilities for drug screening is limited in number. To address this issue, a screening system that could be utilized at lower biosafety levels remains essential. This study aimed to develop a novel SARS-CoV-2 main protease (Mpro) dimer-based screening system (DBSS) utilizing synthetic biology in Escherichia coli BL21(DE3). We linked the SARS-CoV-2 Mpro with the DNA-binding domain of AraC regulatory protein, which regulates the reporter gene expression. Protein modeling and molecular docking showed that saquinavir could bind to AraC-Mpro both in its monomer and dimer forms. The constructed DBSS assay indicated the screening system could detect saquinavir inhibitory activity at a concentration range of 4-10 μg/mL compared to the untreated control (P ≤ 0.05). The Vero E6 cell assay validated the DBSS result that saquinavir at 4-10 μg/mL exhibited antiviral activity against SARS-CoV-2. Our DBSS could be used for preliminary screening of numerous drug candidates that possess a dimerization inhibitor activity of SARS-CoV-2 Mpro and also minimize the use of a high biosafety level laboratory.

New Insights into the Cooperativity and Dynamics of Dimeric Enzymes

A survey of protein databases indicates that the majority of enzymes exist in oligomeric forms, with about half of those found in the UniProt database being homodimeric. Understanding why many enzymes are in their dimeric form is imperative. Recent developments in experimental and computational techniques have allowed for a deeper comprehension of the cooperative interactions between the subunits of dimeric enzymes. This review aims to succinctly summarize these recent advancements by providing an overview of experimental and theoretical methods, as well as an understanding of cooperativity in substrate binding and the molecular mechanisms of cooperative catalysis within homodimeric enzymes. Focus is set upon the beneficial effects of dimerization and cooperative catalysis. These advancements not only provide essential case studies and theoretical support for comprehending dimeric enzyme catalysis but also serve as a foundation for designing highly efficient catalysts, such as dimeric organic catalysts. Moreover, these developments have significant implications for drug design, as exemplified by Paxlovid, which was designed for the homodimeric main protease of SARS-CoV-2.

Systematic Analyses of the Resistance Potential of Drugs Targeting SARS-CoV-2 Main Protease

Drugs that target the main protease (Mpro) of SARS-CoV-2 are effective therapeutics that have entered clinical use. Wide-scale use of these drugs will apply selection pressure for the evolution of resistance mutations. To understand resistance potential in Mpro, we performed comprehensive surveys of amino acid changes that can cause resistance to nirmatrelvir (Pfizer), and ensitrelvir (Xocova) in a yeast screen. We identified 142 resistance mutations for nirmatrelvir and 177 for ensitrelvir, many of which have not been previously reported. Ninety-nine mutations caused apparent resistance to both inhibitors, suggesting likelihood for the evolution of cross-resistance. The mutation with the strongest drug resistance score against nirmatrelvir in our study (E166V) was the most impactful resistance mutation recently reported in multiple viral passaging studies. Many mutations that exhibited inhibitor-specific resistance were consistent with the distinct interactions of each inhibitor in the substrate binding site. In addition, mutants with strong drug resistance scores tended to have reduced function. Our results indicate that strong pressure from nirmatrelvir or ensitrelvir will select for multiple distinct-resistant lineages that will include both primary resistance mutations that weaken interactions with drug while decreasing enzyme function and compensatory mutations that increase enzyme activity. The comprehensive identification of resistance mutations enables the design of inhibitors with reduced potential of developing resistance and aids in the surveillance of drug resistance in circulating viral populations.

Transmissible SARS-CoV-2 variants with resistance to clinical protease inhibitors

Vaccines and drugs have helped reduce disease severity and blunt the spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). However, ongoing virus transmission, continuous evolution, and increasing selective pressures have the potential to yield viral variants capable of resisting these interventions. Here, we investigate the susceptibility of natural variants of the main protease [Mpro; 3C-like protease (3CLpro)] of SARS-CoV-2 to protease inhibitors. Multiple single amino acid changes in Mpro confer resistance to nirmatrelvir (the active component of Paxlovid). An additional clinical-stage inhibitor, ensitrelvir (Xocova), shows a different resistance mutation profile. Importantly, phylogenetic analyses indicate that several of these resistant variants have pre-existed the introduction of these drugs into the human population and are capable of spreading. These results encourage the monitoring of resistance variants and the development of additional protease inhibitors and other antiviral drugs with different mechanisms of action and resistance profiles for combinatorial therapy.

Mutation Effects on Structure and Dynamics: Adaptive Evolution of the SARS-CoV-2 Main Protease

The main protease of SARS-CoV-2 (Mpro) plays a critical role in viral replication; although it is relatively conserved, Mpro has nevertheless evolved over the course of the COVID-19 pandemic. Here, we examine phenotypic changes in clinically observed variants of Mpro, relative to the originally reported wild-type enzyme. Using atomistic molecular dynamics simulations, we examine effects of mutation on protein structure and dynamics. In addition to basic structural properties such as variation in surface area and torsion angles, we use protein structure networks and active site networks to evaluate functionally relevant characters related to global cohesion and active site constraint. Substitution analysis shows a continuing trend toward more hydrophobic residues that are dependent on the location of the residue in primary, secondary, tertiary, and quaternary structures. Phylogenetic analysis provides additional evidence for the impact of selective pressure on mutation of Mpro. Overall, these analyses suggest evolutionary adaptation of Mpro toward more hydrophobicity and a less-constrained active site in response to the selective pressures of a novel host environment.

Transmissible SARS-CoV-2 variants with resistance to clinical protease inhibitors

Vaccines and drugs have helped reduce disease severity and blunt the spread of SARS-CoV-2. However, ongoing virus transmission, continuous evolution, and increasing selective pressures have the potential to yield viral variants capable of resisting these interventions. Here, we investigate the susceptibility of natural variants of the main protease (Mpro/3CLpro) of SARS-CoV-2 to protease inhibitors. Multiple single amino acid changes in Mpro confer resistance to nirmatrelvir (the active component of Paxlovid). An additional clinical-stage inhibitor, ensitrelvir (Xocova), shows a different resistance mutation profile. Importantly, phylogenetic analyses indicate that several of these resistant variants have pre-existed the introduction of these drugs into the human population and are capable of spreading. These results encourage the monitoring of resistance variants and the development of additional protease inhibitors and other antiviral drugs with different mechanisms of action and resistance profiles for combinatorial therapy.

pH profiles of 3-chymotrypsin-like protease (3CLpro) from SARS-CoV-2 elucidate its catalytic mechanism and a histidine residue critical for activity

3-Chymotrypsin-like protease (3CLpro) is a promising drug target for coronavirus disease 2019 and related coronavirus diseases because of the essential role of this protease in processing viral polyproteins after infection. Understanding the detailed catalytic mechanism of 3CLpro is essential for designing effective inhibitors of infection by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Molecular dynamics studies have suggested pH-dependent conformational changes of 3CLpro, but experimental pH profiles of SARS-CoV-2 3CLpro and analyses of the conserved active-site histidine residues have not been reported. In this work, pH-dependence studies of the kinetic parameters of SARS-CoV-2 3CLpro revealed a bell-shaped pH profile with 2 pK_a values (6.9 ± 0.1 and 9.4 ± 0.1) attributable to ionization of the catalytic dyad His41 and Cys145, respectively. Our investigation of the roles of conserved active-site histidines showed that different amino acid substitutions of His163 produced inactive enzymes, indicating a key role of His163 in maintaining catalytically active SARS-CoV-2 3CLpro. By contrast, the H164A and H172A mutants retained 75% and 26% of the activity of WT, respectively. The alternative amino acid substitutions H172K and H172R did not recover the enzymatic activity, whereas H172Y restored activity to a level similar to that of the WT enzyme. The pH profiles of H164A, H172A, and H172Y were similar to those of the WT enzyme, with comparable pK_a values for the catalytic dyad. Taken together, the experimental data support a general base mechanism of SARS-CoV-2 3CLpro and indicate that the neutral states of the catalytic dyad and active-site histidine residues are required for maximum enzyme activity.

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.

Functional map of SARS-CoV-2 3CL protease reveals tolerant and immutable sites

The SARS-CoV-2 3CL protease (3CLpro) is an attractive therapeutic target, as it is essential to the virus and highly conserved among coronaviruses. However, our current understanding of its tolerance to mutations is limited. Here, we develop a yeast-based deep mutational scanning approach to systematically profile the activity of all possible single mutants of the 3CLpro and validate a subset of our results within authentic viruses. We reveal that the 3CLpro is highly malleable and is capable of tolerating mutations throughout the protein. Yet, we also identify specific residues that appear immutable, suggesting that these may be targets for future 3CLpro inhibitors. Finally, we utilize our screening as a basis to identify E166V as a resistance-conferring mutation against the clinically used 3CLpro inhibitor, nirmatrelvir. Collectively, the functional map presented herein may serve as a guide to better understand the biological properties of the 3CLpro and for drug development against coronaviruses.

A genetically encoded BRET-based SARS-CoV-2 Mpro protease activity sensor

The main protease, Mpro, is critical for SARS-CoV-2 replication and an appealing target for designing anti-SARS-CoV-2 agents. Therefore, there is a demand for the development of improved sensors to monitor its activity. Here, we report a pair of genetically encoded, bioluminescence resonance energy transfer (BRET)-based sensors for detecting Mpro proteolytic activity in live cells as well as in vitro. The sensors were generated by sandwiching peptides containing the Mpro N-terminal autocleavage sites, either AVLQSGFR (short) or KTSAVLQSGFRKME (long), in between the mNeonGreen and NanoLuc proteins. Co-expression of the sensors with Mpro in live cells resulted in their cleavage while mutation of the critical C145 residue (C145A) in Mpro completely abrogated their cleavage. Additionally, the sensors recapitulated the inhibition of Mpro by the well-characterized pharmacological agent GC376. Further, in vitro assays with the BRET-based Mpro sensors revealed a molecular crowding-mediated increase in the rate of Mpro activity and a decrease in the inhibitory potential of GC376. The sensors developed here will find direct utility in studies related to drug discovery targeting the SARS-CoV-2 Mpro and functional genomics application to determine the effect of sequence variation in Mpro.

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.

Comprehensive fitness landscape of SARS-CoV-2 Mpro reveals insights into viral resistance mechanisms

With the continual evolution of new strains of SARS-CoV-2 that are more virulent, transmissible, and able to evade current vaccines, there is an urgent need for effective anti-viral drugs SARS-CoV-2 main protease (M^pro) is a leading target for drug design due to its conserved and indispensable role in the viral life cycle. Drugs targeting M^pro appear promising but will elicit selection pressure for resistance. To understand resistance potential in M^pro, we performed a comprehensive mutational scan of the protease that analyzed the function of all possible single amino acid changes. We developed three separate high-throughput assays of M^pro function in yeast, based on either the ability of M^pro variants to cleave at a defined cut-site or on the toxicity of their expression to yeast. We used deep sequencing to quantify the functional effects of each variant in each screen. The protein fitness landscapes from all three screens were strongly correlated, indicating that they captured the biophysical properties critical to M^pro function. The fitness landscapes revealed a non-active site location on the surface that is extremely sensitive to mutation making it a favorable location to target with inhibitors. In addition, we found a network of critical amino acids that physically bridge the two active sites of the M^pro dimer. The clinical variants of M^pro were predominantly functional in our screens, indicating that M^pro is under strong selection pressure in the human population. Our results provide predictions of mutations that will be readily accessible to M^pro evolution and that are likely to contribute to drug resistance. This complete mutational guide of M^pro can be used in the design of inhibitors with reduced potential of evolving viral resistance.

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.

Human Superantibodies to 3CLpro Inhibit Replication of SARS-CoV-2 across Variants

Broadly effective and safe anti-coronavirus agent is existentially needed. Major protease (3CL^pro) is a highly conserved enzyme of betacoronaviruses. The enzyme plays pivotal role in the virus replication cycle. Thus, it is a good target of a broadly effective anti-Betacoronavirus agent. In this study, human single-chain antibodies (HuscFvs) of the SARS-CoV-2 3CL^pro were generated using phage display technology. The 3CL^pro-bound phages were used to infect Escherichia coli host for the production the 3CL^pro-bound HuscFvs. Computerized simulation was used to guide the selection of the phage infected-E. coli clones that produced HuscFvs with the 3CL^pro inhibitory potential. HuscFvs of three phage infected-E. coli clones were predicted to form contact interface with residues for 3CL^pro catalytic activity, substrate binding, and homodimerization. These HuscFvs were linked to a cell-penetrating peptide to make them cell-penetrable, i.e., became superantibodies. The superantibodies blocked the 3CL^pro activity in vitro, were not toxic to human cells, traversed across membrane of 3CL^pro-expressing cells to co-localize with the intracellular 3CL^pro and most of all, they inhibited replication of authentic SARS-CoV-2 Wuhan wild type and α, β, δ, and Omicron variants that were tested. The superantibodies should be investigated further towards clinical application as a safe and broadly effective anti-Betacoronavirus agent.

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.

High-throughput design of symmetrical dimeric SARS-CoV-2 main protease: structural and physical insights into hotspots for adaptation and therapeutics

Dimerization of SARS-CoV-2 main protease (M^pro) is a prerequisite for its processing activity. With >2000 mutations already reported in M^pro, SARS-CoV-2 may accumulate mutations in the M^pro dimeric interface to stabilize it further. We employed high-throughput protein design strategies to design the symmetrical dimeric interface of M^pro (300 000 designs) to identify mutational hotspots that render the M^pro more stable. We found that ∼22% of designed mutations that yield stable M^pro dimers already exist in SARS-CoV-2 genomes and are currently circulating. Our multi-parametric analyses highlight potential M^pro mutations that SARS-CoV-2 may develop, providing a foundation for assessing viral adaptation and mutational surveillance.

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.

Structure-Function Characteristics of SARS-CoV-2 Proteases and Their Potential Inhibitors from Microbial Sources

The COVID-19 pandemic, caused by Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), is considered the greatest challenge to the global health community of the century as it continues to expand. This has prompted immediate urgency to discover promising drug targets for the treatment of COVID-19. The SARS-CoV-2 viral proteases, 3-chymotrypsin-like protease (3CLpro) and papain-like cysteine protease (PLpro), have become the promising target to study due to their essential functions in spreading the virus by RNA transcription, translation, protein synthesis, processing and modification, virus replication, and infection of the host. As such, understanding of the structure and function of these two proteases is unavoidable as platforms for the development of inhibitors targeting this protein which further arrest the infection and spread of the virus. While the abundance of reports on the screening of natural compounds such as SARS-CoV-2 proteases inhibitors are available, the microorganisms-based compounds (peptides and non-peptides) remain less studied. Indeed, microorganisms-based compounds are also one of the potent antiviral candidates against COVID-19. Microbes, especially bacteria and fungi, are other resources to produce new drugs as well as nucleosides, nucleotides, and nucleic acids. Thus, we have compiled various reported literature in detail on the structures, functions of the SARS-CoV-2 proteases, and potential inhibitors from microbial sources as assistance to other researchers working with COVID-19. The compounds are also compared to HIV protease inhibitors which suggested the microorganisms-based compounds are advantageous as SARS-CoV2 proteases inhibitors. The information should serve as a platform for further development of COVID-19 drug design strategies.

Selection of SARS-CoV-2 main protease inhibitor using structure-based virtual screening

Background: Conserved domains within SARS coronavirus 2 nonstructural proteins represent key targets for the design of novel inhibitors. Methods: The authors aimed to identify potential SARS coronavirus 2 NSP5 inhibitors using the ZINC database along with structure-based virtual screening and molecular dynamics simulation. Results: Of 13,840 compounds, 353 with robust docking scores were initially chosen, of which ten hit compounds were selected as candidates for detailed analyses. Three compounds were selected as coronavirus NSP5 inhibitors after passing absorption, distribution, metabolism, excretion and toxicity study; root and mean square deviation; and radius of gyration calculations. Conclusion: ZINC000049899562, ZINC000169336666 and ZINC000095542577 are potential NSP5 protease inhibitors that warrant further experimental studies.

Michaelis-like complex of SARS-CoV-2 main protease visualized by room-temperature X-ray crystallography

SARS-CoV-2 emerged at the end of 2019 to cause an unprecedented pandemic of the deadly respiratory disease COVID-19 that continues to date. The viral main protease (Mpro) is essential for SARS-CoV-2 replication and is therefore an important drug target. Understanding the catalytic mechanism of Mpro, a cysteine protease with a catalytic site comprising the noncanonical Cys145-His41 dyad, can help in guiding drug design. Here, a 2.0 Å resolution room-temperature X-ray crystal structure is reported of a Michaelis-like complex of Mpro harboring a single inactivating mutation C145A bound to the octapeptide Ac-SAVLQSGF-CONH2 corresponding to the nsp4/nsp5 autocleavage site. The peptide substrate is unambiguously defined in subsites S5 to S3' by strong electron density. Superposition of the Michaelis-like complex with the neutron structure of substrate-free Mpro demonstrates that the catalytic site is inherently pre-organized for catalysis prior to substrate binding. Induced fit to the substrate is driven by P1 Gln binding in the predetermined subsite S1 and rearrangement of subsite S2 to accommodate P2 Leu. The Michaelis-like complex structure is ideal for in silico modeling of the SARS-CoV-2 Mpro catalytic mechanism.

Precursors of Viral Proteases as Distinct Drug Targets

Viral proteases are indispensable for successful virion maturation, thus making them a prominent drug target. Their enzyme activity is tightly spatiotemporally regulated by expression in the precursor form with little or no activity, followed by activation via autoprocessing. These cleavage events are frequently triggered upon transportation to a specific compartment inside the host cell. Typically, precursor oligomerization or the presence of a co-factor is needed for activation. A detailed understanding of these mechanisms will allow ligands with non-canonical mechanisms of action to be designed, which would specifically modulate the initial irreversible steps of viral protease autoactivation. Binding sites exclusive to the precursor, including binding sites beyond the protease domain, can be exploited. Both inhibition and up-regulation of the proteolytic activity of viral proteases can be detrimental for the virus. All these possibilities are discussed using examples of medically relevant viruses including herpesviruses, adenoviruses, retroviruses, picornaviruses, caliciviruses, togaviruses, flaviviruses, and coronaviruses.

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.

An Integrated Computational and Experimental Approach to Identifying Inhibitors for SARS-CoV-2 3CL Protease

The newly evolved SARS-CoV-2 has caused the COVID-19 pandemic, and the SARS-CoV-2 main protease 3CLpro is essential for the rapid replication of the virus. Inhibiting this protease may open an alternative avenue toward therapeutic intervention. In this work, a computational docking approach was developed to identify potential small-molecule inhibitors for SARS-CoV-2 3CLpro. Totally 288 potential hits were identified from a half-million bioactive chemicals via a protein-ligand docking protocol. To further evaluate the docking results, a quantitative structure activity relationship (QSAR) model of 3CLpro inhibitors was developed based on existing small molecule inhibitors of the 3CLproSARS- CoV- 1 and their corresponding IC50 data. The QSAR model assesses the physicochemical properties of identified compounds and estimates their inhibitory effects on 3CLproSARS- CoV- 2. Seventy-one potential inhibitors of 3CLpro were selected through these computational approaches and further evaluated via an enzyme activity assay. The results show that two chemicals, i.e., 5-((1-([1,1'-biphenyl]-4-yl)-2,5-dimethyl-1H-pyrrol-3-yl)methylene)pyrimidine-2,4,6(1H,3H,5H)-trione and N-(4-((3-(4-chlorophenylsulfonamido)quinoxalin-2-yl)amino)phenyl)acetamide, effectively inhibited 3CLpro SARS-CoV-2 with IC50's of 19 ± 3 μM and 38 ± 3 μM, respectively. The compounds contain two basic structures, pyrimidinetrione and quinoxaline, which were newly found in 3CLpro inhibitor structures and are of high interest for lead optimization. The findings from this work, such as 3CLpro inhibitor candidates and the QSAR model, will be helpful to accelerate the discovery of inhibitors for related coronaviruses that may carry proteases with similar structures to SARS-CoV-2 3CLpro.

Potential 3-chymotrypsin-like cysteine protease cleavage sites in the coronavirus polyproteins pp1a and pp1ab and their possible relevance to COVID-19 vaccine and drug development

Coronavirus (CoV) 3-chymotrypsin (C)-like cysteine protease (3CLpro ) is a target for anti-CoV drug development and drug repurposing because along with papain-like protease, it cleaves CoV-encoded polyproteins (pp1a and pp1ab) into nonstructural proteins (nsps) for viral replication. However, the cleavage sites of 3CLpro and their relevant nsps remain unclear, which is the subject of this perspective. Here, we address the subject from three standpoints. First, we explore the inconsistency in the cleavage sites and relevant nsps across CoVs, and investigate the function of nsp11. Second, we consider the nsp16 mRNA overlapping of the spike protein mRNA, and analyze the effect of this overlapping on mRNA vaccines. Finally, we study nsp12, whose existence depends on ribosomal frameshifting, and investigate whether 3CLpro requires a large number of inhibitors to achieve full inhibition. This perspective helps us to clarify viral replication and is useful for developing anti-CoV drugs with 3CLpro as a target in the current coronavirus disease 2019 (COVID-19) pandemic.

Crystallographic models of SARS-CoV-2 3CLpro: in-depth assessment of structure quality and validation

The appearance at the end of 2019 of the new SARS-CoV-2 coronavirus led to an unprecedented response by the structural biology community, resulting in the rapid determination of many hundreds of structures of proteins encoded by the virus. As part of an effort to analyze and, if necessary, remediate these structures as deposited in the Protein Data Bank (PDB), this work presents a detailed analysis of 81 crystal structures of the main protease 3CLpro, an important target for the design of drugs against COVID-19. The structures of the unliganded enzyme and its complexes with a number of inhibitors were determined by multiple research groups using different experimental approaches and conditions; the resulting structures span 13 different polymorphs representing seven space groups. The structures of the enzyme itself, all determined by molecular replacement, are highly similar, with the exception of one polymorph with a different inter-domain orientation. However, a number of complexes with bound inhibitors were found to pose significant problems. Some of these could be traced to faulty definitions of geometrical restraints for ligands and to the general problem of a lack of such information in the PDB depositions. Several problems with ligand definition in the PDB itself were also noted. In several cases extensive corrections to the models were necessary to adhere to the evidence of the electron-density maps. Taken together, this analysis of a large number of structures of a single, medically important protein, all determined within less than a year using modern experimental tools, should be useful in future studies of other systems of high interest to the biomedical community.

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.

Development of a Cell-Based Luciferase Complementation Assay for Identification of SARS-CoV-2 3CLpro Inhibitors

The 3C-like protease (3CLpro) of SARS-CoV-2 is considered an excellent target for COVID-19 antiviral drug development because it is essential for viral replication and has a cleavage specificity distinct from human proteases. However, drug development for 3CLpro has been hindered by a lack of cell-based reporter assays that can be performed in a BSL-2 setting. Current efforts to identify 3CLpro inhibitors largely rely upon in vitro screening, which fails to account for cell permeability and cytotoxicity of compounds, or assays involving replication-competent virus, which must be performed in a BSL-3 facility. To address these limitations, we have developed a novel cell-based luciferase complementation reporter assay to identify inhibitors of SARS-CoV-2 3CLpro in a BSL-2 setting. The assay is based on a lentiviral vector that co-expresses 3CLpro and two luciferase fragments linked together by a 3CLpro cleavage site. 3CLpro-mediated cleavage results in a loss of complementation and low luciferase activity, whereas inhibition of 3CLpro results in 10-fold higher levels of luciferase activity. The luciferase reporter assay can easily distinguish true 3CLpro inhibition from cytotoxicity, a powerful feature that should reduce false positives during screening. Using the assay, we screened 32 small molecules for activity against SARS-CoV-2 3CLpro, including HIV protease inhibitors, HCV protease inhibitors, and various other compounds that have been reported to inhibit SARS-CoV-2 3CLpro. Of these, only five exhibited significant inhibition of 3CLpro in cells: GC376, boceprevir, Z-FA-FMK, calpain inhibitor XII, and GRL-0496. This assay should greatly facilitate efforts to identify more potent inhibitors of SARS-CoV-2 3CLpro.

Structure-Based Virtual Screening to Discover Potential Lead Molecules for the SARS-CoV-2 Main Protease

The COVID-19 disease is caused by a new strain of the coronavirus family (SARS-CoV-2), and it has affected at present millions of people all over the world. The indispensable role of the main protease (Mpro) in viral replication and gene expression makes this enzyme an attractive drug target. Therefore, inhibition of SARS-CoV-2 Mpro as a proposition to halt virus ingression is being pursued by scientists globally. Here we carried out a study with two objectives: the first being to perform comparative protein sequence and 3D structural analysis to understand the effect of 12 point mutations on the active site. Among these, two mutations, viz., Ser46 and Phe134, were found to cause a significant change at the active sites of SARS-CoV-2. The Ser46 mutation present at the entrance of the S5 subpocket of SARS-CoV-2 increases the contribution of other two hydrophilic residues, while the Phe134 mutation, present in the catalytic cysteine loop, can cause an increase in catalytic efficiency of Mpro by facilitating fast proton transfer from the Cys145 to His41 residue. It was observed that active site remained conserved among Mpro of both SARS-CoVs, except at the entrance of the S5 subpocket, suggesting sustenance of substrate specificity. The second objective was to screen the inhibitory effects of three different data sets (natural products, coronaviruses main protease inhibitors, and FDA-approved drugs) using a structure-based virtual screening approach. A total of 73 hits had a combo score >2.0. Eight different structural scaffold classes were identified, such as one/two tetrahydropyran ring(s), dipeptide/tripeptide/oligopeptide, large (approximately 20 atoms) cyclic peptide, and miscellaneous. The screened hits showed key interactions with subpockets of the active site. Further, molecular dynamics studies of selected screened compounds confirmed their perfect fitting into the subpockets of the active site. This study suggests promising structures that can fit into the SARS-CoV-2 Mpro active site and also offers direction for further lead optimization and rational drug design.

Biochemical and biophysical characterization of the main protease, 3-chymotrypsin-like protease (3CLpro) from the novel coronavirus SARS-CoV 2

Severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2) is responsible for the novel coronavirus disease 2019 (COVID-19). An appealing antiviral drug target is the coronavirus 3C-like protease (3CLpro) that is responsible for the processing of the viral polyproteins and liberation of functional proteins essential for the maturation and infectivity of the virus. In this study, multiple thermal analytical techniques have been implemented to acquire the thermodynamic parameters of 3CLpro at different buffer conditions. 3CLpro exhibited relatively high thermodynamic stabilities over a wide pH range; however, the protease was found to be less stable in the presence of salts. Divalent metal cations reduced the thermodynamic stability of 3CLpro more than monovalent cations; however, altering the ionic strength of the buffer solution did not alter the stability of 3CLpro. Furthermore, the most stable thermal kinetic stability of 3CLpro was recorded at pH 7.5, with the highest enthalpy of activation calculated from the slope of Eyring plot. The biochemical and biophysical properties of 3CLpro explored here may improve the solubility and stability of 3CLpro for optimum conditions for the setup of an enzymatic assay for the screening of inhibitors to be used as lead candidates in the discovery of drugs and design of antiviral therapeutics against COVID-19.

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.

Crystallographic and electrophilic fragment screening of the SARS-CoV-2 main protease

COVID-19, caused by SARS-CoV-2, lacks effective therapeutics. Additionally, no antiviral drugs or vaccines were developed against the closely related coronavirus, SARS-CoV-1 or MERS-CoV, despite previous zoonotic outbreaks. To identify starting points for such therapeutics, we performed a large-scale screen of electrophile and non-covalent fragments through a combined mass spectrometry and X-ray approach against the SARS-CoV-2 main protease, one of two cysteine viral proteases essential for viral replication. Our crystallographic screen identified 71 hits that span the entire active site, as well as 3 hits at the dimer interface. These structures reveal routes to rapidly develop more potent inhibitors through merging of covalent and non-covalent fragment hits; one series of low-reactivity, tractable covalent fragments were progressed to discover improved binders. These combined hits offer unprecedented structural and reactivity information for on-going structure-based drug design against SARS-CoV-2 main protease.

Unraveling the SARS-CoV-2 Main Protease Mechanism Using Multiscale Methods

We present a detailed theoretical analysis of the reaction mechanism of proteolysis catalyzed by the main protease of SARS-CoV-2. Using multiscale simulation methods, we have characterized the interactions established by a peptidic substrate in the active site, and then we have explored the free energy landscape associated with the acylation and deacylation steps of the proteolysis reaction, characterizing the transition states of the process. Our mechanistic proposals can explain most of the experimental observations made on the highly similar ortholog protease of SARS-CoV. We point to some key interactions that may facilitate the acylation process and thus can be crucial in the design of more specific and efficient inhibitors of the main protease activity. In particular, from our results, the P1' residue can be a key factor to improve the thermodynamics and kinetics of the inhibition process.

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.

Treatments for COVID-19: emerging drugs against the coronavirus

The Coronavirus disease 19 (COVID-19) outbreak has been recognized as a global threat to public health. It is caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and no effective therapies currently exist against this novel viral agent. Along with extensive public health measures, an unprecedented global effort in identifying effective drugs for the treatment is being implemented. Potential drug targets are emerging as the result of a fast-evolving understanding of SARS-CoV-2 virology, host response to the infection, and clinical course of the disease. This brief review focuses on the latest and most promising pharmacological treatments against COVID-19 currently under investigation and discuss their potential use based on either documented efficacy in similar viral infections, or their activity against inflammatory syndromes. Ongoing clinical trials are also emphasized.

NBCZone: Universal three-dimensional construction of eleven amino acids near the catalytic nucleophile and base in the superfamily of (chymo)trypsin-like serine fold proteases

(Chymo)trypsin-like serine fold proteases belong to the serine/cysteine proteases found in eukaryotes, prokaryotes, and viruses. Their catalytic activity is carried out using a triad of amino acids, a nucleophile, a base, and an acid. For this superfamily of proteases, we propose the existence of a universal 3D structure comprising 11 amino acids near the catalytic nucleophile and base - Nucleophile-Base Catalytic Zone (NBCZone). The comparison of NBCZones among 169 eukaryotic, prokaryotic, and viral (chymo)trypsin-like proteases suggested the existence of 15 distinct groups determined by the combination of amino acids located at two "key" structure-functional positions 54_T and 55_T near the catalytic base His57_T. Most eukaryotic and prokaryotic proteases fell into two major groups, [ST]A and TN. Usually, proteases of [ST]A group contain a disulfide bond between cysteines Cys42_T and Cys58_T of the NBCZone. In contrast, viral proteases were distributed among seven groups, and lack this disulfide bond. Furthermore, only the [ST]A group of eukaryotic proteases contains glycine at position 43_T, which is instrumental for activation of these enzymes. In contrast, due to the side chains of residues at position 43_T prokaryotic and viral proteases do not have the ability to carry out the structural transition of the eukaryotic zymogen-zyme type.

Prediction of the SARS-CoV-2 (2019-nCoV) 3C-like protease (3CL pro) structure: virtual screening reveals velpatasvir, ledipasvir, and other drug repurposing candidates

We prepared the three-dimensional model of the SARS-CoV-2 (aka 2019-nCoV) 3C-like protease (3CL pro) using the crystal structure of the highly similar (96% identity) ortholog from the SARS-CoV. All residues involved in the catalysis, substrate binding and dimerisation are 100% conserved. Comparison of the polyprotein PP1AB sequences showed 86% identity. The 3C-like cleavage sites on the coronaviral polyproteins are highly conserved. Based on the near-identical substrate specificities and high sequence identities, we are of the opinion that some of the previous progress of specific inhibitors development for the SARS-CoV enzyme can be conferred on its SARS-CoV-2 counterpart.  With the 3CL pro molecular model, we performed virtual screening for purchasable drugs and proposed 16 candidates for consideration. Among these, the antivirals ledipasvir or velpatasvir are particularly attractive as therapeutics to combat the new coronavirus with minimal side effects, commonly fatigue and headache.  The drugs Epclusa (velpatasvir/sofosbuvir) and Harvoni (ledipasvir/sofosbuvir) could be very effective owing to their dual inhibitory actions on two viral enzymes.

Prediction of the SARS-CoV-2 (2019-nCoV) 3C-like protease (3CL pro) structure: virtual screening reveals velpatasvir, ledipasvir, and other drug repurposing candidates

We prepared the three-dimensional model of the SARS-CoV-2 (aka 2019-nCoV) 3C-like protease (3CL pro) using the crystal structure of the highly similar (96% identity) ortholog from the SARS-CoV. All residues involved in the catalysis, substrate binding and dimerisation are 100% conserved. Comparison of the polyprotein PP1AB sequences showed 86% identity. The 3C-like cleavage sites on the coronaviral polyproteins are highly conserved. Based on the near-identical substrate specificities and high sequence identities, we are of the opinion that some of the previous progress of specific inhibitors development for the SARS-CoV enzyme can be conferred on its SARS-CoV-2 counterpart.  With the 3CL pro molecular model, we performed virtual screening for purchasable drugs and proposed 16 candidates for consideration. Among these, the antivirals ledipasvir or velpatasvir are particularly attractive as therapeutics to combat the new coronavirus with minimal side effects, commonly fatigue and headache.  The drugs Epclusa (velpatasvir/sofosbuvir) and Harvoni (ledipasvir/sofosbuvir) could be very effective owing to their dual inhibitory actions on two viral enzymes.

Structurally- and dynamically-driven allostery of the chymotrypsin-like proteases of SARS, Dengue and Zika viruses

Coronavirus 3C-like and Flavivirus NS2B-NS3 proteases utilize the chymotrypsin fold to harbor their catalytic machineries but also contain additional domains/co-factors. Over the past decade, we aimed to decipher how the extra domains/co-factors mediate the catalytic machineries of SARS 3C-like, Dengue and Zika NS2B-NS3 proteases by characterizing their folding, structures, dynamics and inhibition with NMR, X-ray crystallography and MD simulations, and the results revealed: 1) the chymotrypsin fold of the SARS 3C-like protease can independently fold, while, by contrast, those of Dengue and Zika proteases lack the intrinsic capacity to fold without co-factors. 2) Mutations on the extra domain of SARS 3C-like protease can transform the active catalytic machinery into the inactive collapsed state by structurally-driven allostery. 3) Amazingly, even without detectable structural changes, mutations on the extra domain are sufficient to either inactivate or enhance the catalytic machinery of SARS 3C-like protease by dynamically-driven allostery. 4) Global networks of correlated motions have been identified: for SARS 3C-like protease, N214A inactivates the catalytic machinery by decoupling the network, while STI/A and STIF/A enhance by altering the patterns of the network. The global networks of Dengue and Zika proteases are coordinated by their NS2B-cofactors. 5) Natural products were identified to allosterically inhibit Zika and Dengue proteases through binding a pocket on the back of the active site. Therefore, by introducing extra domains/cofactors, nature develops diverse strategies to regulate the catalytic machinery embedded on the chymotrypsin fold through folding, structurally- and dynamically-driven allostery, all of which might be exploited to develop antiviral drugs.

X-Ray Structure and Inhibition of 3C-like Protease from Porcine Epidemic Diarrhea Virus

Porcine epidemic diarrhea virus (PEDV) is a coronavirus that infects pigs and can have mortality rates approaching 100% in piglets, causing serious economic impact. The 3C-like protease (3CL(pro)) is essential for the coronaviral life cycle and is an appealing target for the development of therapeutics. We report the expression, purification, crystallization and 2.10 Å X-ray structure of 3CL(pro) from PEDV. Analysis of the PEDV 3CL(pro) structure and comparison to other coronaviral 3CL(pro)'s from the same alpha-coronavirus phylogeny shows that the overall structures and active site architectures across 3CL(pro)'s are conserved, with the exception of a loop that comprises the protease S2 pocket. We found a known inhibitor of severe acute respiratory syndrome coronavirus (SARS-CoV) 3CL(pro), (R)-16, to have inhibitor activity against PEDV 3CL(pro), despite that SARS-3CL(pro) and PEDV 3CL(pro) share only 45.4% sequence identity. Structural comparison reveals that the majority of residues involved in (R)-16 binding to SARS-3CL(pro) are conserved in PEDV-3CL(pro); however, the sequence variation and positional difference in the loop forming the S2 pocket may account for large observed difference in IC50 values. This work advances our understanding of the subtle, but important, differences in coronaviral 3CL(pro) architecture and contributes to the broader structural knowledge of coronaviral 3CL(pro)'s.

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.

Ligand-induced Dimerization of Middle East Respiratory Syndrome (MERS) Coronavirus nsp5 Protease (3CLpro): IMPLICATIONS FOR nsp5 REGULATION AND THE DEVELOPMENT OF ANTIVIRALS

All coronaviruses, including the recently emerged Middle East respiratory syndrome coronavirus (MERS-CoV) from the β-CoV subgroup, require the proteolytic activity of the nsp5 protease (also known as 3C-like protease, 3CL^pro) during virus replication, making it a high value target for the development of anti-coronavirus therapeutics. Kinetic studies indicate that in contrast to 3CL^pro from other β-CoV 2c members, including HKU4 and HKU5, MERS-CoV 3CL^pro is less efficient at processing a peptide substrate due to MERS-CoV 3CL^pro being a weakly associated dimer. Conversely, HKU4, HKU5, and SARS-CoV 3CL^pro enzymes are tightly associated dimers. Analytical ultracentrifugation studies support that MERS-CoV 3CL^pro is a weakly associated dimer (K_d ∼52 μm) with a slow off-rate. Peptidomimetic inhibitors of MERS-CoV 3CL^pro were synthesized and utilized in analytical ultracentrifugation experiments and demonstrate that MERS-CoV 3CL^pro undergoes significant ligand-induced dimerization. Kinetic studies also revealed that designed reversible inhibitors act as activators at a low compound concentration as a result of induced dimerization. Primary sequence comparisons and x-ray structural analyses of two MERS-CoV 3CLpro and inhibitor complexes, determined to 1.6 Å, reveal remarkable structural similarity of the dimer interface with 3CL^pro from HKU4-CoV and HKU5-CoV. Despite this structural similarity, substantial differences in the dimerization ability suggest that long range interactions by the nonconserved amino acids distant from the dimer interface may control MERS-CoV 3CL^pro dimerization. Activation of MERS-CoV 3CL^pro through ligand-induced dimerization appears to be unique within the genogroup 2c and may potentially increase the complexity in the development of MERS-CoV 3CL^pro inhibitors as antiviral agents.

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.

Development of potent dipeptide-type SARS-CoV 3CL protease inhibitors with novel P3 scaffolds: design, synthesis, biological evaluation, and docking studies

We report the design and synthesis of a series of dipeptide-type inhibitors with novel P3 scaffolds that display potent inhibitory activity against SARS-CoV 3CLpro. A docking study involving binding between the dipeptidic lead compound 4 and 3CLpro suggested the modification of a structurally flexible P3 N-(3-methoxyphenyl)glycine with various rigid P3 moieties in 4. The modifications led to the identification of several potent derivatives, including 5c-k and 5n with the inhibitory activities (Ki or IC50) in the submicromolar to nanomolar range. Compound 5h, in particular, displayed the most potent inhibitory activity, with a Ki value of 0.006 μM. This potency was 65-fold higher than the potency of the lead compound 4 (Ki=0.39 μM). In addition, the Ki value of 5h was in very good agreement with the binding affinity (16 nM) observed in isothermal titration calorimetry (ITC). A SAR study around the P3 group in the lead 4 led to the identification of a rigid indole-2-carbonyl unit as one of the best P3 moieties (5c). Further optimization showed that a methoxy substitution at the 4-position on the indole unit was highly favorable for enhancing the inhibitory potency.

Chimeric exchange of coronavirus nsp5 proteases (3CLpro) identifies common and divergent regulatory determinants of protease activity

Human coronaviruses (CoVs) such as severe acute respiratory syndrome CoV (SARS-CoV) and Middle East respiratory syndrome CoV (MERS-CoV) cause epidemics of severe human respiratory disease. A conserved step of CoV replication is the translation and processing of replicase polyproteins containing 16 nonstructural protein domains (nsp's 1 to 16). The CoV nsp5 protease (3CLpro; Mpro) processes nsp's at 11 cleavage sites and is essential for virus replication. CoV nsp5 has a conserved 3-domain structure and catalytic residues. However, the intra- and intermolecular determinants of nsp5 activity and their conservation across divergent CoVs are unknown, in part due to challenges in cultivating many human and zoonotic CoVs. To test for conservation of nsp5 structure-function determinants, we engineered chimeric betacoronavirus murine hepatitis virus (MHV) genomes encoding nsp5 proteases of human and bat alphacoronaviruses and betacoronaviruses. Exchange of nsp5 proteases from HCoV-HKU1 and HCoV-OC43, which share the same genogroup, genogroup 2a, with MHV, allowed for immediate viral recovery with efficient replication albeit with impaired fitness in direct competition with wild-type MHV. Introduction of MHV nsp5 temperature-sensitive mutations into chimeric HKU1 and OC43 nsp5 proteases resulted in clear differences in viability and temperature-sensitive phenotypes compared with MHV nsp5. These data indicate tight genetic linkage and coevolution between nsp5 protease and the genomic background and identify differences in intramolecular networks regulating nsp5 function. Our results also provide evidence that chimeric viruses within coronavirus genogroups can be used to test nsp5 determinants of function and inhibition in common isogenic backgrounds and cell types.

Design, synthesis, and biological evaluation of novel dipeptide-type SARS-CoV 3CL protease inhibitors: structure-activity relationship study

This work describes the design, synthesis, and evaluation of low-molecular weight peptidic SARS-CoV 3CL protease inhibitors. The inhibitors were designed based on the potent tripeptidic Z-Val-Leu-Ala(pyrrolidone-3-yl)-2-benzothiazole (8; K_i = 4.1 nM), in which the P3 valine unit was substituted with a variety of distinct moieties. The resulting series of dipeptide-type inhibitors displayed moderate to good inhibitory activities against 3CL^pro. In particular, compounds 26m and 26n exhibited good inhibitory activities with K_i values of 0.39 and 0.33 μM, respectively. These low-molecular weight compounds are attractive leads for the further development of potent peptidomimetic inhibitors with pharmaceutical profiles. Docking studies were performed to model the binding interaction of the compound 26m with the SARS-CoV 3CL protease. The preliminary SAR study of the peptidomimetic compounds with potent inhibitory activities revealed several structural features that boosted the inhibitory activity: (i) a benzothiazole warhead at the S1′ position, (ii) a γ-lactam unit at the S1-position, (iii) an appropriately hydrophobic leucine moiety at the S2-position, and (iv) a hydrogen bond between the N-arylglycine unit and a backbone hydrogen bond donor at the S3-position.

Design and synthesis of new tripeptide-type SARS-CoV 3CL protease inhibitors containing an electrophilic arylketone moiety

We describe here the design, synthesis and biological evaluation of a series of molecules toward the development of novel peptidomimetic inhibitors of SARS-CoV 3CL(pro). A docking study involving binding between the initial lead compound 1 and the SARS-CoV 3CL(pro) motivated the replacement of a thiazole with a benzothiazole unit as a warhead moiety at the P1' site. This modification led to the identification of more potent derivatives, including 2i, 2k, 2m, 2o, and 2p, with IC(50) or K(i) values in the submicromolar to nanomolar range. In particular, compounds 2i and 2p exhibited the most potent inhibitory activities, with K(i) values of 4.1 and 3.1 nM, respectively. The peptidomimetic compounds identified through this process are attractive leads for the development of potential therapeutic agents against SARS. The structural requirements of the peptidomimetics with potent inhibitory activities against SARS-CoV 3CL(pro) may be summarized as follows: (i) the presence of a benzothiazole warhead at the S1'-position; (ii) hydrogen bonding capabilities at the cyclic lactam of the S1-site; (iii) appropriate stereochemistry and hydrophobic moiety size at the S2-site and (iv) a unique folding conformation assumed by the phenoxyacetyl moiety at the S4-site.