Structures of two coronavirus main proteases: implications for substrate binding and antiviral drug design

Coronaviruses (CoVs) can infect humans and multiple species of animals, causing a wide spectrum of diseases. The coronavirus main protease (M^pro), which plays a pivotal role in viral gene expression and replication through the proteolytic processing of replicase polyproteins, is an attractive target for anti-CoV drug design. In this study, the crystal structures of infectious bronchitis virus (IBV) M^pro and a severe acute respiratory syndrome CoV (SARS-CoV) M^pro mutant (H41A), in complex with an N-terminal autocleavage substrate, were individually determined to elucidate the structural flexibility and substrate binding of M^pro. A monomeric form of IBV M^pro was identified for the first time in CoV M^pro structures. A comparison of these two structures to other available M^pro structures provides new insights for the design of substrate-based inhibitors targeting CoV M^pros. Furthermore, a Michael acceptor inhibitor (named N3) was cocrystallized with IBV M^pro and was found to demonstrate in vitro inactivation of IBV M^pro and potent antiviral activity against IBV in chicken embryos. This provides a feasible animal model for designing wide-spectrum inhibitors against CoV-associated diseases. The structure-based optimization of N3 has yielded two more efficacious lead compounds, N27 and H16, with potent inhibition against SARS-CoV M^pro.     

Profiling of substrate specificities of 3C-like proteases from group 1, 2a, 2b, and 3 coronaviruses

Background

Coronaviruses (CoVs) can be classified into alphacoronavirus (group 1), betacoronavirus (group 2), and gammacoronavirus (group 3) based on diversity of the protein sequences. Their 3C-like protease (3CL^pro), which catalyzes the proteolytic processing of the polyproteins for viral replication, is a potential target for anti-coronaviral infection.

Methodology/Principal Findings

Here, we profiled the substrate specificities of 3CL^pro from human CoV NL63 (group 1), human CoV OC43 (group 2a), severe acute respiratory syndrome coronavirus (SARS-CoV) (group 2b) and infectious bronchitis virus (IBV) (group 3), by measuring their activity against a substrate library of 19×8 of variants with single substitutions at P5 to P3'' positions. The results were correlated with structural properties like side chain volume, hydrophobicity, and secondary structure propensities of substituting residues. All 3CL^pro prefer Gln at P1 position, Leu at P2 position, basic residues at P3 position, small hydrophobic residues at P4 position, and small residues at P1'' and P2'' positions. Despite 3CL^pro from different groups of CoVs share many similarities in substrate specificities, differences in substrate specificities were observed at P4 positions, with IBV 3CL^pro prefers P4-Pro and SARS-CoV 3CL^pro prefers P4-Val. By combining the most favorable residues at P3 to P5 positions, we identified super-active substrate sequences ‘VARLQ↓SGF’ that can be cleaved efficiently by all 3CL^pro with relative activity of 1.7 to 3.2, and ‘VPRLQ↓SGF’ that can be cleaved specifically by IBV 3CL^pro with relative activity of 4.3.

Conclusions/Significance

The comprehensive substrate specificities of 3CL^pro from each of the group 1, 2a, 2b, and 3 CoVs have been profiled in this study, which may provide insights into a rational design of broad-spectrum peptidomimetic inhibitors targeting the proteases.     

Design of wide-spectrum inhibitors targeting coronavirus main proteases

The genus Coronavirus contains about 25 species of coronaviruses (CoVs), which are important pathogens causing highly prevalent diseases and often severe or fatal in humans and animals. No licensed specific drugs are available to prevent their infection. Different host receptors for cellular entry, poorly conserved structural proteins (antigens), and the high mutation and recombination rates of CoVs pose a significant problem in the development of wide-spectrum anti-CoV drugs and vaccines. CoV main proteases (M^pros), which are key enzymes in viral gene expression and replication, were revealed to share a highly conservative substrate-recognition pocket by comparison of four crystal structures and a homology model representing all three genetic clusters of the genus Coronavirus. This conclusion was further supported by enzyme activity assays. Mechanism-based irreversible inhibitors were designed, based on this conserved structural region, and a uniform inhibition mechanism was elucidated from the structures of M^pro-inhibitor complexes from severe acute respiratory syndrome-CoV and porcine transmissible gastroenteritis virus. A structure-assisted optimization program has yielded compounds with fast in vitro inactivation of multiple CoV M^pros, potent antiviral activity, and extremely low cellular toxicity in cell-based assays. Further modification could rapidly lead to the discovery of a single agent with clinical potential against existing and possible future emerging CoV-related diseases.     

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.     

Critical Assessment of the Important Residues Involved in the Dimerization and Catalysis of MERS Coronavirus Main Protease

Background

A highly pathogenic human coronavirus (CoV), Middle East respiratory syndrome coronavirus (MERS-CoV), has emerged in Jeddah and other places in Saudi Arabia, and has quickly spread to European and Asian countries since September 2012. Up to the 1^st October 2015 it has infected at least 1593 people with a global fatality rate of about 35%. Studies to understand the virus are necessary and urgent. In the present study, MERS-CoV main protease (M^pro) is expressed; the dimerization of the protein and its relationship to catalysis are investigated.

Methods and Results

The crystal structure of MERS-CoV M^pro indicates that it shares a similar scaffold to that of other coronaviral M^pro and consists of chymotrypsin-like domains I and II and a helical domain III of five helices. Analytical ultracentrifugation analysis demonstrated that MERS-CoV M^pro undergoes a monomer to dimer conversion in the presence of a peptide substrate. Glu169 is a key residue and plays a dual role in both dimerization and catalysis. The mutagenesis of other residues found on the dimerization interface indicate that dimerization of MERS-CoV M^pro is required for its catalytic activity. One mutation, M298R, resulted in a stable dimer with a higher level of proteolytic activity than the wild-type enzyme.

Conclusions

MERS-CoV M^pro shows substrate-induced dimerization and potent proteolytic activity. A critical assessment of the residues important to these processes provides insights into the correlation between dimerization and catalysis within the coronaviral M^pro family.     

C‐terminal domain of SARS‐CoV main protease can form a 3D domain‐swapped dimer

SARS coronavirus main protease (M^pro) plays an essential role in the extensive proteolytic processing of the viral polyproteins (pp1a and pp1ab), and it is an important target for anti-SARS drug development. We have reported that both the M^pro C-terminal domain alone (M^pro-C) and the N-finger deletion mutant of M^pro (M^pro-Δ7) exist as a stable dimer and a stable monomer (Zhong et al., J Virol 2008; 82:4227-4234). Here, we report structures of both M^pro-C monomer and dimer. The structure of the M^pro-C monomer is almost identical to that of the C-terminal domain in the crystal structure of M^pro. Interestingly, the M^pro-C dimer structure is characterized by 3D domain-swapping, in which the first helices of the two protomers are interchanged and each is enwrapped by four other helices from the other protomer. Each folding subunit of the M^pro-C domain-swapped dimer still has the same general fold as that of the M^pro-C monomer. This special dimerization elucidates the structural basis for the observation that there is no exchange between monomeric and dimeric forms of M^pro-C and M^pro-Δ7.     

Recent progress in the discovery of inhibitors targeting coronavirus proteases

Coronaviruses (CoVs) can cause highly prevalent diseases in humans and animals. The fatal outbreak of severe acute respiratory syndrome (SARS) and Middle East respiratory syndrome (MERS) highlights the threat posed by this unique virus subfamily. However, no specific drugs have been approved to treat CoV-associated diseases to date. The CoV proteases, which play pivotal roles in viral gene expression and replication through a highly complex cascade involving the proteolytic processing of replicase polyproteins, are attractive targets for drug design. This review summarizes the recent advances in biological and structural studies, together with the development of inhibitors targeting CoV proteases, particularly main proteases (M^pros), which could help develop effective treatments to prevent CoV infection.

Open image in new window

     

Improved SARS-CoV-2 main protease high-throughput screening assay using a 5-carboxyfluorescein substrate

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) as a global threat to human health has highlighted the need for the development of novel therapies targeting current and emerging coronaviruses with pandemic potential. The coronavirus main protease (M^pro, also called 3CL^pro) is a validated drug target against coronaviruses and has been heavily studied since the emergence of SARS-CoV-2 in late 2019. Here, we report the biophysical and enzymatic characterization of native M^pro, then characterize the steady-state kinetics of several commonly used FRET substrates, fluorogenic substrates, and six of the 11 reported SARS-CoV-2 polyprotein cleavage sequences. We then assessed the suitability of these substrates for high-throughput screening. Guided by our assessment of these substrates, we developed an improved 5-carboxyfluorescein-based FRET substrate, which is better suited for high-throughput screening and is less susceptible to interference and false positives than existing substrates. This study provides a useful framework for the design of coronavirus M^pro enzyme assays to facilitate the discovery and development of therapies targeting M^pro.     

N-Terminal Finger Stabilizes the S1 Pocket for the Reversible Feline Drug GC376 in the SARS-CoV-2 Mpro Dimer

The main protease (M^pro, also known as 3CL protease) of SARS-CoV-2 is a high priority drug target in the development of antivirals to combat COVID-19 infections. A feline coronavirus antiviral drug, GC376, has been shown to be effective in inhibiting the SARS-CoV-2 main protease and live virus growth. As this drug moves into clinical trials, further characterization of GC376 with the main protease of coronaviruses is required to gain insight into the drug’s properties, such as reversibility and broad specificity. Reversibility is an important factor for therapeutic proteolytic inhibitors to prevent toxicity due to off-target effects. Here we demonstrate that GC376 has nanomolar K_i values with the M^pro from both SARS-CoV-2 and SARS-CoV strains. Restoring enzymatic activity after inhibition by GC376 demonstrates reversible binding with both proteases. In addition, the stability and thermodynamic parameters of both proteases were studied to shed light on physical chemical properties of these viral enzymes, revealing higher stability for SARS-CoV-2 M^pro. The comparison of a new X-ray crystal structure of M^pro from SARS-CoV complexed with GC376 reveals similar molecular mechanism of inhibition compared to SARS-CoV-2 M^pro, and gives insight into the broad specificity properties of this drug. In both structures, we observe domain swapping of the N-termini in the dimer of the M^pro, which facilitates coordination of the drug’s P1 position. These results validate that GC376 is a drug with an off-rate suitable for clinical trials.     

1H, 13C and 15N resonance assignments of SARS-CoV main protease N-terminal domain

The main protease (M^pro) of severe acute respiratory syndrome coronavirus (SARS-CoV) plays an essential role in the extensive proteolytic processing of the viral polyproteins (pp1a and pp1ab), and it is an important target for anti-SARS drug development. SARS-CoV M^pro is composed of a catalytic N-terminal domain and an α-helical C-terminal domain linked by a long loop. Even though the N-terminal domain of SARS-CoV M^pro adopts a similar chymotrypsin-like fold as that of piconavirus 3C protease, the extra C-terminal domain is required for SARS-CoV M^pro to be enzymatically active. Here, we reported the NMR assignments of the SARS-CoV M^pro N-terminal domain alone, which are essential for its solution structure determination.     

Identification of natural compounds as potent inhibitors of SARS-CoV-2 main protease using combined docking and molecular dynamics simulations

Coronavirus disease 2019 (COVID-19) has emerged from China and globally affected the entire population through the human-to-human transmission of a newly emerged virus called severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The genome of SARS-CoV-2 encodes several proteins that are essential for multiplication and pathogenesis. The main protease (M^pro or 3CL^pro) of SARS-CoV-2 plays a central role in its pathogenesis and thus is considered as an attractive drug target for the drug design and development of small-molecule inhibitors. We have employed an extensive structure-based high-throughput virtual screening to discover potential natural compounds from the ZINC database which could inhibit the M^pro of SARS-CoV-2. Initially, the hits were selected on the basis of their physicochemical and drug-like properties. Subsequently, the PAINS filter, estimation of binding affinities using molecular docking, and interaction analyses were performed to find safe and potential inhibitors of SARS-CoV-2 M^pro. We have identified ZINC02123811 (1-(3-(2,5,9-trimethyl-7-oxo-3-phenyl-7H-furo[3,2-g]chromen-6-yl)propanoyl)piperidine-4-carboxamide), a natural compound bearing appreciable affinity, efficiency, and specificity towards the binding pocket of SARS-CoV-2 M^pro. The identified compound showed a set of drug-like properties and preferentially binds to the active site of SARS-CoV-2 M^pro. All-atom molecular dynamics (MD) simulations were performed to evaluate the conformational dynamics, stability and interaction mechanism of M^pro with ZINC02123811. MD simulation results indicated that M^pro with ZINC02123811 forms a stable complex throughout the trajectory of 100 ns. These findings suggest that ZINC02123811 may be further exploited as a promising scaffold for the development of potential inhibitors of SARS-CoV-2 M^pro to address COVID-19.     

The Molecular Mechanism of Domain Swapping of the C-Terminal Domain of the SARS-Coronavirus Main Protease

In three-dimensional domain swapping, two protein monomers exchange a part of their structures to form an intertwined homodimer, whose subunits resemble the monomer. Several viral proteins domain swap to increase their structural complexity or functional avidity. The main protease (M^pro) of the severe acute respiratory syndrome (SARS) coronavirus proteolyzes viral polyproteins and has been a target for anti-SARS drug design. Domain swapping in the α-helical C-terminal domain of M^pro (M^proC) locks M^pro into a hyperactive octameric form that is hypothesized to promote the early stages of viral replication. However, in the absence of a complete molecular understanding of the mechanism of domain swapping, investigations into the biological relevance of this octameric M^pro have stalled. Isolated M^proC can exist as a monomer or a domain-swapped dimer. Here, we investigate the mechanism of domain swapping of M^proC using coarse-grained structure-based models and molecular dynamics simulations. Our simulations recapitulate several experimental features of M^proC folding. Further, we find that a contact between a tryptophan in the M^proC domain-swapping hinge and an arginine elsewhere forms early during folding, modulates the folding route, and promotes domain swapping to the native structure. An examination of the sequence and the structure of the tryptophan containing hinge loop shows that it has a propensity to form multiple secondary structures and contacts, indicating that it could be stabilized into either the monomer- or dimer-promoting conformations by mutations or ligand binding. Finally, because all residues in the tryptophan loop are identical in SARS-CoV and SARS-CoV-2, mutations that modulate domain swapping may provide insights into the role of octameric M^pro in the early-stage viral replication of both viruses.     

A COVID moonshot: assessment of ligand binding to the SARS-CoV-2 main protease by saturation transfer difference NMR spectroscopy

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the etiological cause of the coronavirus disease 2019, for which no effective antiviral therapeutics are available. The SARS-CoV-2 main protease (M^pro) is essential for viral replication and constitutes a promising therapeutic target. Many efforts aimed at deriving effective M^pro inhibitors are currently underway, including an international open-science discovery project, codenamed COVID Moonshot. As part of COVID Moonshot, we used saturation transfer difference nuclear magnetic resonance (STD-NMR) spectroscopy to assess the binding of putative M^pro ligands to the viral protease, including molecules identified by crystallographic fragment screening and novel compounds designed as M^pro inhibitors. In this manner, we aimed to complement enzymatic activity assays of M^pro performed by other groups with information on ligand affinity. We have made the M^pro STD-NMR data publicly available. Here, we provide detailed information on the NMR protocols used and challenges faced, thereby placing these data into context. Our goal is to assist the interpretation of M^pro STD-NMR data, thereby accelerating ongoing drug design efforts.

     

Individual and common inhibitors of coronavirus and picornavirus main proteases

Picornaviruses (PV) and coronaviruses (CoV) are positive-stranded RNA viruses which infect millions of people worldwide each year, resulting in a wide range of clinical outcomes. As reported in this study, using high throughput screening against ∼6800 small molecules, we have identified several novel inhibitors of SARS-CoV 3CL^pro with IC₅₀ of low μM. Interestingly, one of them equally inhibited both 3C^pro and 3CL^pro from PV and CoV, respectively. Using computer modeling, the structural features of these compounds as individual and common protease inhibitors were elucidated to enhance our knowledge for developing anti-viral agents against PV and CoV.     

Structural Basis of Main Proteases of Coronavirus Bound to Drug Candidate PF-07304814

New variants of the severe acute respiratory syndrome Coronavirus 2 (SARS-CoV-2) emerged and spread rapidly all over the world, which strongly supports the need for pharmacological options to complement vaccine strategies. Main protease (M^pro or 3CL^pro) is a critical enzyme in the life cycle of SARS-CoV-2 and appears to be highly conserved among different genera of coronaviruses, making it an ideal target for the development of drugs with broad-spectrum property. PF-07304814 developed by Pfizer is an intravenously administered inhibitor targeting SARS-CoV-2 M^pro. Here we showed that PF-07304814 displays broad-spectrum inhibitory activity against M^pros from multiple coronaviruses. Crystal structures of M^pros of SARS-CoV-2, SARS-CoV, MERS-CoV, and HCoV-NL63 bound to the inhibitor PF-07304814 revealed a conserved ligand-binding site, providing new insights into the mechanism of inhibition of viral replication. A detailed analysis of these crystal structures complemented by comprehensive comparison defined the key structural determinants essential for inhibition and illustrated the binding mode of action of M^pros from different coronaviruses. In view of the importance of M^pro for the medications of SARS-CoV-2 infection, insights derived from the present study should accelerate the design of pan-coronaviral main protease inhibitors that are safer and more effective.     

Two-way antigenic cross-reactivity between severe acute respiratory syndrome coronavirus (SARS-CoV) and group 1 animal CoVs is mediated through an antigenic site in the N-terminal region of the SARS-CoV nucleoprotein

In 2002, severe acute respiratory syndrome-associated coronavirus (SARS-CoV) emerged in humans, causing a global epidemic. By phylogenetic analysis, SARS-CoV is distinct from known CoVs and most closely related to group 2 CoVs. However, no antigenic cross-reactivity between SARS-CoV and known CoVs was conclusively and consistently demonstrated except for group 1 animal CoVs. We analyzed this cross-reactivity by an enzyme-linked immunosorbent assay (ELISA) and Western blot analysis using specific antisera to animal CoVs and SARS-CoV and SARS patient convalescent-phase or negative sera. Moderate two-way cross-reactivity between SARS-CoV and porcine CoVs (transmissible gastroenteritis CoV [TGEV] and porcine respiratory CoV [PRCV]) was mediated through the N but not the spike protein, whereas weaker cross-reactivity occurred with feline (feline infectious peritonitis virus) and canine CoVs. Using Escherichia coli-expressed recombinant SARS-CoV N protein and fragments, the cross-reactive region was localized between amino acids (aa) 120 to 208. The N-protein fragments comprising aa 360 to 412 and aa 1 to 213 reacted specifically with SARS convalescent-phase sera but not with negative human sera in ELISA; the fragment comprising aa 1 to 213 cross-reacted with antisera to animal CoVs, whereas the fragment comprising aa 360 to 412 did not cross-react and could be a potential candidate for SARS diagnosis. Particularly noteworthy, a single substitution at aa 120 of PRCV N protein diminished the cross-reactivity. We also demonstrated that the cross-reactivity is not universal for all group 1 CoVs, because HCoV-NL63 did not cross-react with SARS-CoV. One-way cross-reactivity of HCoV-NL63 with group 1 CoVs was localized to aa 1 to 39 and at least one other antigenic site in the N-protein C terminus, differing from the cross-reactive region identified in SARS-CoV N protein. The observed cross-reactivity is not a consequence of a higher level of amino acid identity between SARS-CoV and porcine CoV nucleoproteins, because sequence comparisons indicated that SARS-CoV N protein has amino acid identity similar to that of infectious bronchitis virus N protein and shares a higher level of identity with bovine CoV N protein within the cross-reactive region. The TGEV and SARS-CoV N proteins are RNA chaperons with long disordered regions. We speculate that during natural infection, antibodies target similar short antigenic sites within the N proteins of SARS-CoV and porcine group 1 CoVs that are exposed to an immune response. Identification of the cross-reactive and non-cross-reactive N-protein regions allows development of SARS-CoV-specific antibody assays for screening animal and human sera.     

Structural and functional characterization of MERS coronavirus papain-like protease

Backgrounds

A new highly pathogenic human coronavirus (CoV), Middle East respiratory syndrome coronavirus (MERS-CoV), has emerged in Jeddah and Saudi Arabia and quickly spread to some European countries since September 2012. Until 15 May 2014, it has infected at least 572 people with a fatality rate of about 30% globally. Studies to understand the virus and to develop antiviral drugs or therapy are necessary and urgent. In the present study, MERS-CoV papain-like protease (PL^pro) is expressed, and its structural and functional consequences are elucidated.

Results

Circular dichroism and Tyr/Trp fluorescence analyses indicated that the secondary and tertiary structure of MERS-CoV PL^pro is well organized and folded. Analytical ultracentrifugation analyses demonstrated that MERS-CoV PL^pro is a monomer in solution. The steady-state kinetic and deubiquitination activity assays indicated that MERS-CoV PL^pro exhibits potent deubiquitination activity but lower proteolytic activity, compared with SARS-CoV PL^pro. A natural mutation, Leu105, is the major reason for this difference.

Conclusions

Overall, MERS-CoV PL^pro bound by an endogenous metal ion shows a folded structure and potent proteolytic and deubiquitination activity. These findings provide important insights into the structural and functional properties of coronaviral PL^pro family, which is applicable to develop strategies inhibiting PL^pro against highly pathogenic coronaviruses.     

Structural Bases of Coronavirus Attachment to Host Aminopeptidase N and Its Inhibition by Neutralizing Antibodies

The coronaviruses (CoVs) are enveloped viruses of animals and humans associated mostly with enteric and respiratory diseases, such as the severe acute respiratory syndrome and 10–20% of all common colds. A subset of CoVs uses the cell surface aminopeptidase N (APN), a membrane-bound metalloprotease, as a cell entry receptor. In these viruses, the envelope spike glycoprotein (S) mediates the attachment of the virus particles to APN and subsequent cell entry, which can be blocked by neutralizing antibodies. Here we describe the crystal structures of the receptor-binding domains (RBDs) of two closely related CoV strains, transmissible gastroenteritis virus (TGEV) and porcine respiratory CoV (PRCV), in complex with their receptor, porcine APN (pAPN), or with a neutralizing antibody. The data provide detailed information on the architecture of the dimeric pAPN ectodomain and its interaction with the CoV S. We show that a protruding receptor-binding edge in the S determines virus-binding specificity for recessed glycan-containing surfaces in the membrane-distal region of the pAPN ectodomain. Comparison of the RBDs of TGEV and PRCV to those of other related CoVs, suggests that the conformation of the S receptor-binding region determines cell entry receptor specificity. Moreover, the receptor-binding edge is a major antigenic determinant in the TGEV envelope S that is targeted by neutralizing antibodies. Our results provide a compelling view on CoV cell entry and immune neutralization, and may aid the design of antivirals or CoV vaccines. APN is also considered a target for cancer therapy and its structure, reported here, could facilitate the development of anti-cancer drugs.