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.     

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.     

In silico structure-based discovery of a SARS-CoV-2 main protease inhibitor

The Coronavirus Disease 2019 (COVID-19) pandemic caused by the novel lineage B betacoroanvirus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has resulted in significant mortality, morbidity, and socioeconomic disruptions worldwide. Effective antivirals are urgently needed for COVID-19. The main protease (M^pro) of SARS-CoV-2 is an attractive antiviral target because of its essential role in the cleavage of the viral polypeptide. In this study, we performed an in silico structure-based screening of a large chemical library to identify potential SARS-CoV-2 M^pro inhibitors. Among 8,820 compounds in the library, our screening identified trichostatin A, a histone deacetylase inhibitor and an antifungal compound, as an inhibitor of SARS-CoV-2 M^pro activity and replication. The half maximal effective concentration of trichostatin A against SARS-CoV-2 replication was 1.5 to 2.7µM, which was markedly below its 50% effective cytotoxic concentration (75.7µM) and peak serum concentration (132µM). Further drug compound optimization to develop more stable analogues with longer half-lives should be performed. This structure-based drug discovery platform should facilitate the identification of additional enzyme inhibitors of SARS-CoV-2.     

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.     

Chalcones isolated from Angelica keiskei inhibit cysteine proteases of SARS-CoV

Two viral proteases of severe acute respiratory syndrome coronavirus (SARS-CoV), a chymotrypsin-like protease (3CL^pro) and a papain-like protease (PL^pro) are attractive targets for the development of anti-SARS drugs. In this study, nine alkylated chalcones (1–9) and four coumarins (10–13) were isolated from Angelica keiskei, and the inhibitory activities of these constituents against SARS-CoV proteases (3CL^pro and PL^pro) were determined (cell-free/based). Of the isolated alkylated chalcones, chalcone 6, containing the perhydroxyl group, exhibited the most potent 3CL^pro and PL^pro inhibitory activity with IC₅₀ values of 11.4 and 1.2?µM. Our detailed protein-inhibitor mechanistic analysis of these species indicated that the chalcones exhibited competitive inhibition characteristics to the SARS-CoV 3CL^pro, whereas noncompetitive inhibition was observed with the SARS-CoV PL^pro.     

Tanshinones as selective and slow-binding inhibitors for SARS-CoV cysteine proteases

In the search for anti-SARS-CoV, tanshinones derived from Salvia miltiorrhiza were found to be specific and selective inhibitors for the SARS-CoV 3CL^pro and PL^pro, viral cysteine proteases. A literature search for studies involving the seven isolated tanshinone hits showed that at present, none have been identified as coronaviral protease inhibitors. We have identified that all of the isolated tanshinones are good inhibitors of both cysteine proteases. However, their activity was slightly affected by subtle changes in structure and targeting enzymes. All isolated compounds (1–7) act as time dependent inhibitors of PL^pro, but no improved inhibition was observed following preincubation with the 3CL^pro. In a detail kinetic mechanism study, all of the tanshinones except rosmariquinone (7) were identified as noncompetitive enzyme isomerization inhibitors. However, rosmariquinone (7) showed a different kinetic mechanism through mixed-type simple reversible slow-binding inhibition. Furthermore, tanshinone I (5) exhibited the most potent nanomolar level inhibitory activity toward deubiquitinating (IC₅₀ = 0.7 μM). Additionally, the inhibition is selective because these compounds do not exert significant inhibitory effects against other proteases including chymotrysin, papain, and HIV protease. These findings provide potential inhibitors for SARS-CoV viral infection and replication.     

Crystal Structure of the Middle East Respiratory Syndrome Coronavirus (MERS-CoV) Papain-like Protease Bound to Ubiquitin Facilitates Targeted Disruption of Deubiquitinating Activity to Demonstrate Its Role in Innate Immune Suppression

Middle East respiratory syndrome coronavirus (MERS-CoV) is a newly emerging human pathogen that was first isolated in 2012. MERS-CoV replication depends in part on a virus-encoded papain-like protease (PL^pro) that cleaves the viral replicase polyproteins at three sites releasing non-structural protein 1 (nsp1), nsp2, and nsp3. In addition to this replicative function, MERS-CoV PL^pro was recently shown to be a deubiquitinating enzyme (DUB) and to possess deISGylating activity, as previously reported for other coronaviral PL^pro domains, including that of severe acute respiratory syndrome coronavirus. These activities have been suggested to suppress host antiviral responses during infection. To understand the molecular basis for ubiquitin (Ub) recognition and deconjugation by MERS-CoV PL^pro, we determined its crystal structure in complex with Ub. Guided by this structure, mutations were introduced into PL^pro to specifically disrupt Ub binding without affecting viral polyprotein cleavage, as determined using an in trans nsp3↓4 cleavage assay. Having developed a strategy to selectively disable PL^pro DUB activity, we were able to specifically examine the effects of this activity on the innate immune response. Whereas the wild-type PL^pro domain was found to suppress IFN-β promoter activation, PL^pro variants specifically lacking DUB activity were no longer able to do so. These findings directly implicate the DUB function of PL^pro, and not its proteolytic activity per se, in the inhibition of IFN-β promoter activity. The ability to decouple the DUB activity of PL^pro from its role in viral polyprotein processing now provides an approach to further dissect the role(s) of PL^pro as a viral DUB during MERS-CoV infection.     

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.

     

A Crystallographic Snapshot of SARS-CoV-2 Main Protease Maturation Process

SARS-CoV-2 is the causative agent of COVID-19. The dimeric form of the viral M^pro is responsible for the cleavage of the viral polyprotein in 11 sites, including its own N and C-terminus. The lack of structural information for intermediary forms of M^pro is a setback for the understanding its self-maturation process. Herein, we used X-ray crystallography combined with biochemical data to characterize multiple forms of SARS-CoV-2 M^pro. For the immature form, we show that extra N-terminal residues caused conformational changes in the positioning of domain-three over the active site, hampering the dimerization and diminishing its activity. We propose that this form preludes the cis and trans-cleavage of N-terminal residues. Using fragment screening, we probe new cavities in this form which can be used to guide therapeutic development. Furthermore, we characterized a serine site-directed mutant of the M^pro bound to its endogenous N and C-terminal residues during dimeric association stage of the maturation process. We suggest this form is a transitional state during the C-terminal trans-cleavage. This data sheds light in the structural modifications of the SARS-CoV-2 main protease during its self-maturation process.     

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.     

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.     

Evaluation of polyphenols from Broussonetia papyrifera as coronavirus protease inhibitors

The current study was designed to assess the inhibitory activity of Broussonetia papyrifera-derived polyphenols against 3-chymotrypsin-like and papain-like coronavirus cysteine proteases. The isolated compounds were broussochalcone B (1), broussochalcone A (2), 4-hydroxyisolonchocarpin (3), papyriflavonol A (4), 3′-(3-methylbut-2-enyl)-3′,4,7-trihydroxyflavane (5), kazinol A (6), kazinol B (7), broussoflavan A (8), kazinol F (9), and kazinol J (10). All polyphenols were more potent against papain-like protease (PL^pro) than against 3-chymotripsin-like protease (3CL^pro); therefore, we investigated their structural features that were responsible for this selectivity. Compound 4 was the most potent inhibitor of PL^pro with an IC₅₀ value of 3.7?μM. The active compounds displayed kinetic behaviors, and the binding constants of their interaction with PL^pro were determined from surface plasmon resonance analysis. Our results suggest B. papyrifera constituents as promising candidates for development into potential anti-coronaviral agents.     

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

The main protease (3CL M^pro) from SARS–CoV-2, the etiological agent of COVID-19, is an essential enzyme for viral replication. 3CL M^pro possesses an unusual catalytic dyad composed of Cys¹⁴⁵ and His⁴¹ 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 M^pro, 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 Cys¹⁴⁵ is in the negatively charged thiolate state and His⁴¹ 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.     

Identification of quinone analogues as potential inhibitors of picornavirus 3C protease in vitro

Picornaviruses are non-enveloped viruses that represent a large family of positive-sense single-stranded RNA viruses including a number of causative agents of many human and animal diseases such as coxsackievirus B3 (CVB3) and rhinoviruses (HRV). In this study, we performed a high-throughput screening of a compound library composed of ～6000 small molecules in search of potential picornavirus 3C protease (3C^pro) inhibitors. As results, we identified quinone analogues that effectively inhibited both CVB3 3C^pro and HRV 3C^pro with IC₅₀ values in low micromolar range. Together with predicted binding modes of these compounds to the active site of the viral protease, it is implied that structural features of these non-peptidic inhibitors may act as useful scaffold for further anti-picornavirus drug design and development.     

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.     

Identification and design of novel small molecule inhibitors against MERS-CoV papain-like protease via high-throughput screening and molecular modeling

The development of new therapeutic agents against the coronavirus causing Middle East Respiratory Syndrome (MERS) is a continuing imperative. The initial MERS-CoV epidemic was contained entirely through public health measures, but episodic cases continue, as there are currently no therapeutic agents effective in the treatment of MERS-CoV, although multiple strategies have been proposed. In this study, we screened 30,000 compounds from three different compound libraries against one of the essential proteases, the papain-like protease (PL^pro), using a fluorescence-based enzymatic assay followed by surface plasmon resonance (SPR) direct binding analysis for hit confirmation. Mode of inhibition assays and competition SPR studies revealed two compounds to be competitive inhibitors. To improve upon the inhibitory activity of the best hit compounds, a small fragment library consisting of 352 fragments was screened in the presence of each hit compound, resulting in one fragment that enhanced the IC₅₀ value of the best hit compound by 3-fold. Molecular docking and MM/PBSA binding energy calculations were used to predict potential binding sites, providing insight for design and synthesis of next-generation compounds.     

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.     

Structural basis for the in vitro efficacy of nirmatrelvir against SARS-CoV-2 variants

The COVID-19 pandemic continues to be a public health threat with emerging variants of SARS-CoV-2. Nirmatrelvir (PF-07321332) is a reversible, covalent inhibitor targeting the main protease (M^pro) of SARS-CoV-2 and the active protease inhibitor in PAXLOVID (nirmatrelvir tablets and ritonavir tablets). However, the efficacy of nirmatrelvir is underdetermined against evolving SARS-CoV-2 variants. Here, we evaluated the in vitro catalytic activity and potency of nirmatrelvir against the M^pro of prevalent variants of concern (VOCs) or variants of interest (VOIs): Alpha (α, B.1.1.7), Beta (β, B.1.351), Delta (δ, B1.617.2), Gamma (γ, P.1), Lambda (λ, B.1.1.1.37/C37), Omicron (ο, B.1.1.529), as well as the original Washington or wildtype strain. These VOCs/VOIs carry prevalent mutations at varying frequencies in the M^pro specifically for α, β, γ (K90R), λ (G15S), and ο (P132H). In vitro biochemical enzymatic assay characterization of the enzyme kinetics of the mutant M^pros demonstrates that they are catalytically comparable to wildtype. We found that nirmatrelvir has similar potency against each mutant M^pro including P132H that is observed in the Omicron variant with a Ki of 0.635 nM as compared to a Ki of 0.933 nM for wildtype. The molecular basis for these observations were provided by solution-phase structural dynamics and structural determination of nirmatrelvir bound to the ο, λ, and β M^pro at 1.63 to 2.09 Å resolution. These in vitro data suggest that PAXLOVID has the potential to maintain plasma concentrations of nirmatrelvir many-fold times higher than the amount required to stop the SARS-CoV-2 VOC/VOI, including Omicron, from replicating in cells.