Molecular dynamic simulations analysis of ritonavir and lopinavir as SARS-CoV 3CL(pro) inhibitors

Since the emergence of the severe acute respiratory syndrome (SARS) to date, neither an effective antiviral drug nor a vaccine against SARS is available. However, it was found that a mixture of two HIV-1 proteinase inhibitors, lopinavir and ritonavir, exhibited some signs of effectiveness against the SARS virus. To understand the fine details of the molecular interactions between these proteinase inhibitors and the SARS virus via complexation, molecular dynamics simulations were carried out for the SARS-CoV 3CL^pro free enzyme (free SARS) and its complexes with lopinavir (SARS-LPV) and ritonavir (SARS-RTV). The results show that flap closing was clearly observed when the inhibitors bind to the active site of SARS-CoV 3CL^pro. The binding affinities of LPV and RTV to SARS-CoV 3CL^pro do not show any significant difference. In addition, six hydrogen bonds were detected in the SARS-LPV system, while seven hydrogen bonds were found in SARS-RTV complex.     

Dieckol,a SARS-CoV 3CLpro inhibitor,isolated from the edible brown algae Ecklonia cava

SARS-CoV 3CL^pro plays an important role in viral replication. In this study, we performed a biological evaluation on nine phlorotannins isolated from the edible brown algae Ecklonia cava. The nine isolated phlorotannins (1–9), except phloroglucinol (1), possessed SARS-CoV 3CL^pro inhibitory activities in a dose-dependently and competitive manner. Of these phlorotannins (1–9), two eckol groups with a diphenyl ether linked dieckol (8) showed the most potent SARS-CoV 3CL^pro trans/cis-cleavage inhibitory effects (IC₅₀s = 2.7 and 68.1 μM, respectively). This is the first report of a (8) phlorotannin chemotype significantly blocking the cleavage of SARS-CoV 3CL^pro in a cell-based assay with no toxicity. Furthermore, dieckol (8) exhibited a high association rate in the SPR sensorgram and formed extremely strong hydrogen bonds to the catalytic dyad (Cys145 and His41) of the SARS-CoV 3CL^pro.     

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.     

Virtual screening identification of novel severe acute respiratory syndrome 3C-like protease inhibitors and in vitro confirmation

The 3C-like protease (3CL^pro) of severe acute respiratory syndrome associated coronavirus (SARS-CoV) is vital for SARS-CoV replication and is a promising drug target. Structure based virtual screening of 308 307 chemical compounds was performed using the computation tool Autodock 3.0.5 on a WISDOM Production Environment. The top 1468 ranked compounds with free binding energy ranging from −14.0 to −17.09 kcal mol⁻¹ were selected to check the hydrogen bond interaction with amino acid residues in the active site of 3CL^pro. Fifty-three compounds from 35 main groups were tested in an in vitro assay for inhibition of 3CL^pro expressed by Escherichia coli. Seven of the 53 compounds were selected; their IC₅₀ ranged from 38.57 ± 2.41 to 101.38 ± 3.27 μM. Two strong 3CL^pro inhibitors were further identified as competitive inhibitors of 3CL^pro with K_i values of 9.11 ± 1.6 and 9.93 ± 0.44 μM. Hydrophobic and hydrogen bond interactions of compound with amino acid residues in the active site of 3CL^pro were also identified.     

Synthesis,modification and docking studies of 5-sulfonyl isatin derivatives as SARS-CoV 3C-like protease inhibitors

The Severe Acute Respiratory Syndrome (SARS) is a serious life-threatening and strikingly mortal respiratory illness caused by SARS-CoV. SARS-CoV which contains a chymotrypsin-like main protease analogous to that of the main picornavirus protease, 3CL^pro. 3CL^pro plays a pivotal role in the viral replication cycle and is a potential target for SARS inhibitor development. A series of isatin derivatives as possible SARS-CoV 3CL^pro inhibitors was designed, synthesized, and evaluated by in vitro protease assay using fluorogenic substrate peptide, in which several showed potent inhibition against the 3CL^pro. Structure–activity relationship was analyzed, and possible binding interaction modes were proposed by molecular docking studies. Among all compounds, 8k₁ showed most potent inhibitory activity against 3CL^pro (IC₅₀ = 1.04 μM). These results indicated that these inhibitors could be potentially developed into anti-SARS drugs.     

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.     

Scutellaria baicalensis extract and baicalein inhibit replication of SARS-CoV-2 and its 3C-like protease in vitro

COVID-19 has become a global pandemic and there is an urgent call for developing drugs against the virus (SARS-CoV-2). The 3C-like protease (3CL^pro) of SARS-CoV-2 is a preferred target for broad spectrum anti-coronavirus drug discovery. We studied the anti-SARS-CoV-2 activity of S. baicalensis and its ingredients. We found that the ethanol extract of S. baicalensis and its major component, baicalein, inhibit SARS-CoV-2 3CL^pro activity in vitro with IC₅₀’s of 8.52 µg/ml and 0.39 µM, respectively. Both of them inhibit the replication of SARS-CoV-2 in Vero cells with EC₅₀’s of 0.74 µg/ml and 2.9 µM, respectively. While baicalein is mainly active at the viral post-entry stage, the ethanol extract also inhibits viral entry. We further identified four baicalein analogues from other herbs that inhibit SARS-CoV-2 3CL^pro activity at µM concentration. All the active compounds and the S. baicalensis extract also inhibit the SARS-CoV 3CL^pro, demonstrating their potential as broad-spectrum anti-coronavirus drugs.     

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.     

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.     

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₅₀ 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.     

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.     

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.     

Maturation Mechanism of Severe Acute Respiratory Syndrome (SARS) Coronavirus 3C-like Proteinase

The 3C-like proteinase (3CL^pro) of the severe acute respiratory syndrome (SARS) coronavirus plays a vital role in virus maturation and is proposed to be a key target for drug design against SARS. Various in vitro studies revealed that only the dimer of the matured 3CL^pro is active. However, as the internally encoded 3CL^pro gets matured from the replicase polyprotein by autolytic cleavage at both the N-terminal and the C-terminal flanking sites, it is unclear whether the polyprotein also needs to dimerize first for its autocleavage reaction. We constructed a large protein containing the cyan fluorescent protein (C), the N-terminal flanking substrate peptide of SARS 3CL^pro (XX), SARS 3CL^pro (3CLP), and the yellow fluorescent protein (Y) to study the autoprocessing of 3CL^pro using fluorescence resonance energy transfer. In contrast to the matured 3CL^pro, the polyprotein, as well as the one-step digested product, 3CLP-Y-His, were shown to be monomeric in gel filtration and analytic ultracentrifuge analysis. However, dimers can still be induced and detected when incubating these large proteins with a substrate analog compound in both chemical cross-linking experiments and analytic ultracentrifuge analysis. We also measured enzyme activity under different enzyme concentrations and found a clear tendency of substrate-induced dimer formation. Based on these discoveries, we conclude that substrate-induced dimerization is essential for the activity of SARS-3CL^pro in the polyprotein, and a modified model for the 3CL^pro maturation process was proposed. As many viral proteases undergo a similar maturation process, this model might be generally applicable.     

新型冠状病毒3C样蛋白酶结构和功能特征分析

采用生物信息学方法分析新型冠状病毒(Severe acute respiratory syndrome coronavirus 2, SARS-CoV-2)3C样蛋白酶(3-chymotrypsin-like protease, 3CL^pro)的理化性质、结构与功能,为抗SARS-CoV-2药物研发提供参考。通过ProtParam、ProtScale、Bioedit服务器对3CL^pro进行一级结构如氨基酸理化性质、疏水性的预测分析;COILS Server、SignalP、TMPred、TargetP Server、NetPhos Server、NetNGlyc Server服务器对3CL^pro结构进行如卷曲螺旋区、信号肽、跨膜结构域、亚细胞定位、磷酸化位点、糖基化位点的预测分析;SOPMA、SWISS MODEL服务器对3CL^pro进行二级结构、三级结构的预测分析;IEBD对3CL^pro进行B细胞表位的预测分析。3CL^pro由306个氨基酸组成,其中亮氨酸占比最高,分子质量为33 796.64,理论等电点值为5.95,半衰期为1.9 h,脂肪系数为82.12;亲水性较高,不具有卷曲螺旋区与信号肽特点,含一个跨膜区;具有4个磷酸化位点,2个糖基化修饰点;二级结构中无规则卷曲占据主导地位,三级结构能与已知的6y2g.1（SMTL ID）模型同源建模;存在4个潜在的B细胞表位,位于92~101位的氨基酸区域应答频率最高。利用生物信息学技术分析3CL^pro的结构和功能特征,可为新型冠状肺炎药物的研发提供参考。     

Identification,synthesis and evaluation of SARS-CoV and MERS-CoV 3C-like protease inhibitors

Severe acute respiratory syndrome (SARS) led to a life-threatening form of atypical pneumonia in late 2002. Following that, Middle East Respiratory Syndrome (MERS-CoV) has recently emerged, killing about 36% of patients infected globally, mainly in Saudi Arabia and South Korea. Based on a scaffold we reported for inhibiting neuraminidase (NA), we synthesized the analogues and identified compounds with low micromolar inhibitory activity against 3CL^pro of SARS-CoV and MERS-CoV. Docking studies show that a carboxylate present at either R¹ or R⁴ destabilizes the oxyanion hole in the 3CL^pro. Interestingly, 3f, 3g and 3m could inhibit both NA and 3CL^pro and serve as a starting point to develop broad-spectrum antiviral agents.     

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.     

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.     

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.     

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.