Evaluation of peptide-aldehyde inhibitors using R188I mutant of SARS 3CL protease as a proteolysis-resistant mutant

The 3C-like (3CL) protease of the severe acute respiratory syndrome (SARS) coronavirus is a key enzyme for the virus maturation. We found for the first time that the mature SARS 3CL protease is subject to degradation at 188Arg/189Gln. Replacing Arg with Ile at position 188 rendered the protease resistant to proteolysis. The R188I mutant digested a conserved undecapeptide substrate with a K(m) of 33.8 microM and k(cat) of 4753 s(-1). Compared with the value reported for the mature protease containing a C-terminal His-tag, the relative activity of the mutant was nearly 10(6). Novel peptide-aldehyde derivatives containing a side-chain-protected C-terminal Gln efficiently inhibited the catalytic activity of the R188I mutant. The results indicated for the first time that the tetrapeptide sequence is enough for inhibitory activities of peptide-aldehyde derivatives.     

Mechanism of Dissociative Inhibition of HIV Protease and Its Autoprocessing from a Precursor

Dimerization is indispensible for release of the human immunodeficiency virus protease (PR) from its precursor (Gag-Pol) and ensuing mature-like catalytic activity that is crucial for virus maturation. We show that a single-chain Fv fragment (scFv) of a previously reported monoclonal antibody (mAb1696), which recognizes the N-terminus of PR, dissociates a dimeric mature D25N PR mutant with an enhanced dimer dissociation constant (K(d)) in the sub-micromolar range to form predominantly a monomer-scFv complex at a 1:1 ratio, along with small (5-10%) amounts of a dimer-scFv complex. Enzyme kinetics indicate a mixed mechanism of inhibition of the wild-type PR, which exhibits a K(d)<10nM, with effects both on K(m) and k(cat) at an scFv-to-PR ratio of 10:1. ScFv binds to the N-terminal peptide P(1)QITLW(6) of PR and to PR monomers with dissociation constants of ≤30nM and ~100nM, respectively. Consistent with an ~400-fold increase in the dissociation of the antibody (K(Ab)) on even addition of an acetyl group to P(1) of the peptide, the antibody fails to inhibit N-terminal autoprocessing of the PR from a model precursor (at ~5μM). However, subsequent to this cleavage, it sequesters the PR, thus blocking autoprocessing at its C-terminus. A second monoclonal antibody [PRM1 (human monoclonal antibody to PR)], which recognizes part of the flap region (residues 41-47) of the mature PR and its precursor, does not inhibit autoprocessing and ensuing catalytic activity. However, its failure to recognize drug-resistant clinical mutants of PR may be beneficial to monitor the selection of mutations in this region under drug pressure.     

Insight into the activity of SARS main protease: Molecular dynamics study of dimeric and monomeric form of enzyme

The phenomenon that SARS coronavirus main protease (SARS M(pro)) dimer is the main functional form has been confirmed by experiment. However, because of the absence of structural information of the monomer, the reasons for this remain unknown. To investigate it, two molecular dynamics (MD) simulations in water for dimer and monomer models have been carried out, using the crystal structure of protomer A of the dimer as the starting structure for the monomer. During the MD simulation of dimer, three interest phenomena of protomer A have been observed: (i) the distance between NE2 of His41 and SG of Cys145 averages 3.72 A, which agrees well with the experimental observations made by X-ray crystallography; (ii) His163 and Glu166 form the "tooth" conformational properties, resulting in the specificity for glutamine at substrate P1 site; and (iii) the substrate-binding pocket formed by loop 140-146 and loop 184-197 is large enough to accommodate the substrate analog. However, during the MD simulation of the monomer complex, the three structural characteristics are all absent, which results directly in the inactivation of the monomer. Throughout the MD simulation of the dimer, the N-terminus of protomer B forms stable hydrogen bonds with Phe140 and Glu166, through which His163, Glu166, and loop 140-146 are kept active form. Furthermore, a water-bridge has been found between the N-terminus of protomer B and Gly170, which stabilizes His172 and avoids it moving toward Tyr161 to disrupt the H-bond between Tyr161 and His163, stabilizing the conformation of His163. The interactions between the N-terminus and another monomer maintain the activity of dimer.     

Molecular cloning, expression, purification, and mass spectrometric characterization of 3C-like protease of SARS coronavirus

Severe acute respiratory syndrome (SARS) is an acute respiratory illness, which has broken out in China. It has been known that SARS coronavirus (SARS_CoV) is a novel human coronavirus and is responsible for SARS infection. Belonging to one of the major proteins associated with SARS_CoV, SARS 3C-like protease (SARS_3CL(pro)) functions as a cysteine protease engaging in the proteolytic cleavage of the viral precursor polyprotein to a series of functional proteins required for coronavirus replication and is considered as an appealing target for designing anti-SARS agents. To facilitate the studies regarding the functions and structures of SARS_3CL(pro), in this report the synthetic genes encoding 3CL(pro) of SARS_CoV were assembled, and the plasmid was constructed using pQE30 as vector and expressed in Escherichia coli M15 cells. The highly yielded ( approximately 15mg/L) expressed protease was purified by use of NTA-Ni(2+) affinity chromatography and FPLC system, and its sequence was determined by LC/MS with the residue coverage of 46.4%.     

Mutation of Gly-11 on the dimer interface results in the complete crystallographic dimer dissociation of severe acute respiratory syndrome coronavirus 3C-like protease: crystal structure with molecular dynamics simulations

SARS-CoV 3C-like protease (3CL(pro)) is an attractive target for anti-severe acute respiratory syndrome (SARS) drug discovery, and its dimerization has been extensively proved to be indispensable for enzymatic activity. However, the reason why the dissociated monomer is inactive still remains unclear due to the absence of the monomer structure. In this study, we showed that mutation of the dimer-interface residue Gly-11 to alanine entirely abolished the activity of SARS-CoV 3CL(pro). Subsequently, we determined the crystal structure of this mutant and discovered a complete crystallographic dimer dissociation of SARS-CoV 3CL(pro). The mutation might shorten the alpha-helix A' of domain I and cause a mis-oriented N-terminal finger that could not correctly squeeze into the pocket of another monomer during dimerization, thus destabilizing the dimer structure. Several structural features essential for catalysis and substrate recognition are severely impaired in the G11A monomer. Moreover, domain III rotates dramatically against the chymotrypsin fold compared with the dimer, from which we proposed a putative dimerization model for SARS-CoV 3CL(pro). As the first reported monomer structure for SARS-CoV 3CL(pro), the crystal structure of G11A mutant might provide insight into the dimerization mechanism of the protease and supply direct structural evidence for the incompetence of the dissociated monomer.     

Folding pathway mediated by an intramolecular chaperone: dissecting conformational changes coincident with autoprocessing and the role of Ca(2+) in subtilisin maturation.

Subtilisin is produced as a precursor that requires its N-terminal propeptide to chaperone the folding of its protease domain. Once folded, subtilisin adopts a remarkably stable conformation, which has been attributed to a high affinity Ca(2+) binding site. We investigated the role of the metal ligand in the maturation of pro-subtilisin, a process that involves folding, autoprocessing and partial degradation. Our results establish that although Ca(2+) ions can stabilize the protease domain, the folding and autoprocessing of pro-subtilisin take place independent of Ca(2+) ion. We demonstrate that the stabilizing effect of calcium is observed only after the completion of autoprocessing and that the metal ion appears to be responsible for shifting the folding equilibrium towards the native conformation in both mature subtilisin and the autoprocessed propeptide:subtilisin complex. Furthermore, the addition of active subtilisin to unautoprocessed pro-subtilisin in trans does not facilitate precursor maturation, but rather promotes rapid autodegradation. The primary cleavage site that initiates this autodegradation is at Gln19 in the N-terminus of mature subtilisin. This corresponds to the loop that links alpha-helix-2 and beta-strand-1 in mature subtilisin and has indirect effects on the formation of the Ca(2+) binding site. Our results show that the N-terminus of mature subtilisin undergoes rearrangement subsequent to propeptide autoprocessing. Since this structural change enhances the proteolytic stability of the precursor, our results suggest that the autoprocessing reaction must be completed before the release of active subtilisin in order to maximize folding efficiency.     

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.     

Identification of novel inhibitors of the SARS coronavirus main protease 3CLpro

SARS (severe acute respiratory syndrome) is caused by a newly discovered coronavirus. A key enzyme for the maturation of this virus and, therefore, a target for drug development is the main protease 3CL(pro) (also termed SARS-CoV 3CL(pro)). We have cloned and expressed in Escherichia coli the full-length SARS-CoV 3CL(pro) as well as a truncated form containing only the catalytic domains. The recombinant proteins have been characterized enzymatically using a fluorescently labeled substrate; their structural stability in solution has been determined by differential scanning calorimetry, and novel inhibitors have been discovered. Expression of the catalytic region alone yields a protein with a reduced catalytic efficiency consistent with the proposed regulatory role of the alpha-helical domain. Differential scanning calorimetry indicates that the alpha-helical domain does not contribute to the structural stability of the catalytic domains. Analysis of the active site cavity reveals the presence of subsites that can be targeted with specific chemical functionalities. In particular, a cluster of serine residues (Ser139, Ser144, and Ser147) was identified near the active site cavity and was susceptible to being targeted by compounds containing boronic acid. This cluster is highly conserved in similar proteases from other coronaviruses, defining an attractive target for drug development. It was found that bifunctional aryl boronic acid compounds were particularly effective at inhibiting the protease, with inhibition constants as strong as 40 nM. Isothermal titration microcalorimetric experiments indicate that these inhibitors bind reversibly to 3CL(pro) in an enthalpically favorable fashion, implying that they establish strong interactions with the protease molecule, thus defining attractive molecular scaffolds for further optimization.     

In silico analysis of ORF1ab in coronavirus HKU1 genome reveals a unique putative cleavage site of coronavirus HKU1 3C-like protease

Recently we have described the discovery and complete genome sequence of a novel coronavirus associated with pneumonia, coronavirus HKU1 (CoV-HKU1). In this study, a detailed in silico analysis of the ORF1ab, encoding the 7,182-amino acid replicase polyprotein in the CoV-HKU1 genome showed that the replicase polyprotein of CoV-HKU1 is cleaved by its papain-like proteases and 3C-like protease (3CL(pro)) into 16 polypeptides homologous to the corresponding polypeptides in other coronaviruses. Surprisingly, analysis of the putative cleavage sites of the 3CL(pro) revealed a unique putative cleavage site. In all known coronaviruses, the P1 positions at the cleavage sites of the 3CL(pro) are occupied by glutamine. This is also observed in CoV-HKU1, except for one site at the junction between nsp10 (helicase) and nsp11 (member of exonuclease family), where the P1 position is occupied by histidine. This amino acid substitution is due to a single nucleotide mutation in the CoV-HKU1 genome, CAG/A to CAT. This probably represents a novel cleavage site because the same mutation was consistently observed in CoV-HKU1 sequences from multiple specimens of different patients; the P2 and P1'-P12' positions of this cleavage site are consistent between CoV-HKU1 and other coronaviruses; and as the helicase is one of the most conserved proteins in coronaviruses, cleavage between nsp10 and nsp11 should be an essential step for the generation of the mature functional helicase. Experiments, including purification and C-terminal amino acid sequencing of the CoV-HKU1 helicase and trans-cleavage assays of the CoV-HKU1 3CL(pro) will confirm the presence of this novel cleavage site.     

3C-like proteinase from SARS coronavirus catalyzes substrate hydrolysis by a general base mechanism

SARS 3C-like proteinase has been proposed to be a key enzyme for drug design against SARS. Lack of a suitable assay has been a major hindrance for enzyme kinetic studies and a large-scale inhibitor screen for SARS 3CL proteinase. Since SARS 3CL proteinase belongs to the cysteine protease family (family C3 in clan CB) with a chymotrypsin fold, it is important to understand the catalytic mechanism of SARS 3CL proteinase to determine whether the proteolysis proceeds through a general base catalysis mechanism like chymotrypsin or an ion pair mechanism like papain. We have established a continuous colorimetric assay for SARS 3CL proteinase and applied it to study the enzyme catalytic mechanism. The proposed catalytic residues His41 and Cys145 were confirmed to be critical for catalysis by mutating to Ala, while the Cys145 to Ser mutation resulted in an active enzyme with a 40-fold lower activity. From the pH dependency of catalytic activity, the pK(a)'s for His41 and Cys145 in the wild-type enzyme were estimated to be 6.38 and 8.34, while the pK(a)'s for His41 and Ser145 in the C145S mutant were estimated to be 6.15 and 9.09, respectively. The C145S mutant has a normal isotope effect in D(2)O for general base catalysis, that is, reacts slower in D(2)O, while the wild-type enzyme shows an inverse isotope effect which may come from the lower activation enthalpy. The pK(a) values measured for the active site residues and the activity of the C145S mutant are consistent with a general base catalysis mechanism and cannot be explained by a thiolate-imidazolium ion pair model.     

Mutation of Glu-166 Blocks the Substrate-Induced Dimerization of SARS Coronavirus Main Protease

The maturation of SARS coronavirus involves the autocleavage of polyproteins 1a and 1ab by the main protease (Mpro) and a papain-like protease; these represent attractive targets for the development of anti-SARS drugs. The functional unit of Mpro is a homodimer, and each subunit has a His-41?Cys-145 catalytic dyad. Current thinking in this area is that Mpro dimerization is essential for catalysis, although the influence of the substrate binding on the dimer formation has never been explored. Here, we delineate the contributions of the peptide substrate to Mpro dimerization. Enzyme kinetic assays indicate that the monomeric mutant R298A/L exhibits lower activity but in a cooperative manner. Analytical ultracentrifugation analyses indicate that in the presence of substrates, the major species of R298A/L shows a significant size shift toward the dimeric form and the monomer-dimer dissociation constant of R298A/L decreases by 12- to 17-fold, approaching that for wild-type. Furthermore, this substrate-induced dimerization was found to be reversible after substrates were removed. Based on the crystal structures, a key residue, Glu-166, which is responsible for recognizing the Gln-P1 of the substrate and binding to Ser-1 of another protomer, will interact with Asn-142 and block the S1 subsite entrance in the monomer. Our studies indicate that mutation of Glu-166 in the R298A mutant indeed blocks the substrate-induced dimerization. This demonstrates that Glu-166 plays a pivotal role in connecting the substrate binding site with the dimer interface. We conclude that protein-ligand and protein-protein interactions are closely correlated in Mpro.     

Sequential autoprocessing of Marek's disease herpesvirus protease differs from that of other herpesviruses

Herpesviruses encode a unique serine protease essential for viral capsid maturation. This protease undergoes autoprocessing at two sites, R and M, at the consensus sequence (V, L, I)(P3)-X(P2)-A(P1)/S(P1') (where X is a polar amino acid). We observed complete autoprocessing at the R and M sites of Marek's disease virus (MDV) protease following production of the polyprotein in Escherichia coli. Site-directed mutagenesis confirmed the predicted sequence of the R and M sites, with the M site sequence being nonconsensual: M(P3)-N(P2)-A(P1)/S(P1'). Mutagenesis and expression kinetics studies suggested that the atypical MDV M site was cleaved exclusively by the processed short protease, a feature making MDV unique among herpesviruses.     

Folded monomer of HIV-1 protease

The mature human immunodeficiency virus type 1 protease rapidly folds into an enzymatically active stable dimer, exhibiting an intricate interplay between structure formation and dimerization. We now show by NMR and sedimentation equilibrium studies that a mutant protease containing the R87K substitution (PR(R87K)) within the highly conserved Gly(86)-Arg(87)-Asn(88) sequence forms a monomer with a fold similar to a single subunit of the dimer. However, binding of the inhibitor DMP323 to PR(R87K) produces a stable dimer complex. Based on the crystal structure and our NMR results, we postulate that loss of specific interactions involving the side chain of Arg(87) destabilizes PR(R87K) by perturbing the inner C-terminal beta-sheet (residues 96-99 from each monomer), a region that is sandwiched between the two beta-strands formed by the N-terminal residues (residues 1-4) in the mature protease. We systematically examined the folding, dimerization, and catalytic activities of mutant proteases comprising deletions of either one of the terminal regions (residues 1-4 or 96-99) or both. Although both N- and C-terminal beta-strands were found to contribute to dimer stability, our results indicate that the inner C-terminal strands are absolutely essential for dimer formation. Knowledge of the monomer fold and regions critical for dimerization may aid in the rational design of novel inhibitors of the protease to overcome the problem of drug resistance.     

The 3C-like proteinase of an invertebrate nidovirus links coronavirus and potyvirus homologs

Gill-associated virus (GAV), a positive-stranded RNA virus of prawns, is the prototype of newly recognized taxa (genus Okavirus, family Roniviridae) within the order NIDOVIRALES: In this study, a putative GAV cysteine proteinase (3C-like proteinase [3CL(pro)]), which is predicted to be the key enzyme involved in processing of the GAV replicase polyprotein precursors, pp1a and pp1ab, was characterized. Comparative sequence analysis indicated that, like its coronavirus homologs, 3CL(pro) has a three-domain organization and is flanked by hydrophobic domains. The putative 3CL(pro) domain including flanking regions (pp1a residues 2793 to 3143) was fused to the Escherichia coli maltose-binding protein (MBP) and, when expressed in E. coli, was found to possess N-terminal autoprocessing activity that was not dependent on the presence of the 3CL(pro) C-terminal domain. N-terminal sequence analysis of the processed protein revealed that cleavage occurred at the location (2827)LVTHE downward arrow VRTGN(2836). The trans-processing activity of the purified recombinant 3CL(pro) (pp1a residues 2832 to 3126) was used to identify another cleavage site, (6441)KVNHE downward arrow LYHVA(6450), in the C-terminal pp1ab region. Taken together, the data tentatively identify VxHE downward arrow (L,V) as the substrate consensus sequence for the GAV 3CL(pro). The study revealed that the GAV and potyvirus 3CL(pro)s possess similar substrate specificities which correlate with structural similarities in their respective substrate-binding sites, identified in sequence comparisons. Analysis of the proteolytic activities of MBP-3CL(pro) fusion proteins carrying replacements of putative active-site residues provided evidence that, in contrast to most other 3C/3CL(pro)s but in common with coronavirus 3CL(pro)s, the GAV 3CL(pro) employs a Cys(2968)-His(2879) catalytic dyad. The properties of the GAV 3CL(pro) define a novel RNA virus proteinase variant that bridges the gap between the distantly related chymotrypsin-like cysteine proteinases of coronaviruses and potyviruses.     

Only one protomer is active in the dimer of SARS 3C-like proteinase

The severe acute respiratory syndrome coronavirus 3C-like protease has been proposed to be a key target for structurally based drug design against SARS. The enzyme exists as a mixture of dimer and monomer, and only the dimer was considered to be active. In this report, we have investigated, using molecular dynamics simulation and mutational studies, the problems as to why only the dimer is active and whether both of the two protomers in the dimer are active. The molecular dynamics simulations show that the monomers are always inactive, that the two protomers in the dimer are asymmetric, and that only one protomer is active at a time. The enzyme activity of the hybrid severe acute respiratory syndrome coronavirus 3C-like protease of the wild-type protein and the inactive mutant proves that the dimerization is important for enzyme activity and only one active protomer in the dimer is enough for the catalysis. Our simulations also show that the right conformation for catalysis in one protomer can be induced upon dimer formation. These results suggest that the enzyme may follow the association, activation, catalysis, and dissociation mechanism for activity control.     

Polyprotein cleavage mechanism of SARS CoV Mpro and chemical modification of the octapeptide

The cleavage mechanism of severe acute respiratory syndrome (SARS) coronavirus main proteinase (M(pro) or 3CL(pro)) for the octapeptide AVLQSGFR is studied using molecular mechanics (MM) and quantum mechanics (QM). The catalytic dyad His-41 and Cys-145 in the active pocket between domain I and II seem to polarize the pi-electron density of the peptide bond between Gln and Ser in the octapeptide, leading to an increase of positive charge on C(CO) of Gln and negative charge on N(NH) of Ser. The possibility of enhancing the chemical bond between Gln and Ser based on the "distorted key" theory [Anal. Biochem. 233 (1996) 1] is examined. The scissile peptide bond between Gln and Ser is found to be solidified through "hybrid peptide bond" by changing the carbonyl group CO of Gln to CH(2) or CF(2). This leads to a break of the pi-bond system for the peptide bond, making the octapeptide (AVLQSGFR) a "distorted key" and a potential starting system for the design of anti SARS drugs.     

Without its N-finger, the main protease of severe acute respiratory syndrome coronavirus can form a novel dimer through its C-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. It was found that SARS-CoV M(pro) exists in solution as an equilibrium of both monomeric and dimeric forms, and the dimeric form is the enzymatically active form. However, the mechanism of SARS-CoV M(pro) dimerization, especially the roles of its N-terminal seven residues (N-finger) and its unique C-terminal domain in the dimerization, remain unclear. Here we report that the SARS-CoV M(pro) C-terminal domain alone (residues 187 to 306; M(pro)-C) is produced in Escherichia coli in both monomeric and dimeric forms, and no exchange could be observed between them at room temperature. The M(pro)-C dimer has a novel dimerization interface. Meanwhile, the N-finger deletion mutant of SARS-CoV M(pro) also exists as both a stable monomer and a stable dimer, and the dimer is formed through the same C-terminal-domain interaction as that in the M(pro)-C dimer. However, no C-terminal domain-mediated dimerization form can be detected for wild-type SARS-CoV M(pro). Our study results help to clarify previously published controversial claims about the role of the N-finger in SARS-CoV M(pro) dimerization. Apparently, without the N-finger, SARS-CoV M(pro) can no longer retain the active dimer structure; instead, it can form a new type of dimer which is inactive. Therefore, the N-finger of SARS-CoV M(pro) is not only critical for its dimerization but also essential for the enzyme to form the enzymatically active dimer.