key: cord-0933546-71wj0ajq
authors: Zhu, Lili; George, Shyla; Schmidt, Marco F.; Al-Gharabli, Samer I.; Rademann, Jörg; Hilgenfeld, Rolf
title: Peptide aldehyde inhibitors challenge the substrate specificity of the SARS-coronavirus main protease
date: 2011-08-11
journal: Antiviral Res
DOI: 10.1016/j.antiviral.2011.08.001
sha: ecfb3c3040888f2ce3fffad5624f86bb261c7096
doc_id: 933546
cord_uid: 71wj0ajq

SARS coronavirus main protease (SARS-CoV M(pro)) is essential for the replication of the virus and regarded as a major antiviral drug target. The enzyme is a cysteine protease, with a catalytic dyad (Cys-145/His-41) in the active site. Aldehyde inhibitors can bind reversibly to the active-site sulfhydryl of SARS-CoV M(pro). Previous studies using peptidic substrates and inhibitors showed that the substrate specificity of SARS-CoV M(pro) requires glutamine in the P1 position and a large hydrophobic residue in the P2 position. We determined four crystal structures of SARS-CoV M(pro) in complex with pentapeptide aldehydes (Ac-ESTLQ-H, Ac-NSFSQ-H, Ac-DSFDQ-H, and Ac-NSTSQ-H). Kinetic data showed that all of these aldehydes exhibit inhibitory activity towards SARS-CoV M(pro), with K(i) values in the μM range. Surprisingly, the X-ray structures revealed that the hydrophobic S2 pocket of the enzyme can accommodate serine and even aspartic-acid side-chains in the P2 positions of the inhibitors. Consequently, we reassessed the substrate specificity of the enzyme by testing the cleavage of 20 different tetradecapeptide substrates with varying amino-acid residues in the P2 position. The cleavage efficiency for the substrate with serine in the P2 position was 160-times lower than that for the original substrate (P2 = Leu); furthermore, the substrate with aspartic acid in the P2 position was not cleaved at all. We also determined a crystal structure of SARS-CoV M(pro) in complex with aldehyde Cm-FF-H, which has its P1-phenylalanine residue bound to the relatively hydrophilic S1 pocket of the enzyme and yet exhibits a high inhibitory activity against SARS-CoV M(pro), with K(i) = 2.24 ± 0.58 μM. These results show that the stringent substrate specificity of the SARS-CoV M(pro) with respect to the P1 and P2 positions can be overruled by the highly electrophilic character of the aldehyde warhead, thereby constituting a deviation from the dogma that peptidic inhibitors need to correspond to the observed cleavage specificity of the target protease.

SARS coronavirus main protease (SARS-CoV M pro ) is essential for the replication of the virus and regarded as a major antiviral drug target Yang et al., 2003; Tan et al., 2005; Steuber and Hilgenfeld, 2010) . Peptide aldehyde inhibitors are widely used as research tools to characterize the substrate specificity of serine and cysteine proteases. In complex with the crystalline target enzyme, they provide a wealth of information on the specificity-defining subsites of the substrate-binding site. For this purpose, we have previously described peptide aldehydes as inhibitors of the SARS-coronavirus main protease (SARS-CoV M pro ) (Al-Gharabli et al., 2006) . A major advantage of aldehydes over other peptide electrophiles is their reversible binding to the 0166-3542/$ -see front matter Ó 2011 Elsevier B.V. All rights reserved. doi:10.1016/j.antiviral.2011.08.001

Abbreviations: SARS-CoV M pro , SARS coronavirus main protease; Cm-FF-H, Cinnamoyl-Phe-Phe-H; Boc, tert-butyl oxycarbonyl; IBX, 2-iodoxybenzoic acid; FRET, fluorescence resonance energy transfer; PEG, polyethylene glycol; MPD, 2methyl-2,4-pentanediol; MES, 2-(N-morpholino)ethanesulfonic acid; DMSO, dimethylsulfoxide. active-site sulfhydryl of cysteine proteases. We have previously used this property to screen a small library of non-peptide nucleophiles (mostly amines) for molecules that would efficiently enhance the inhibition of the aldehyde by formation of a stronger binding aldehyde-nucleophile ligation product. In a second step of this process that we named ''Dynamic Ligation Screening'' (Schmidt et al., 2008) , the most efficient compound in this process had its amino group replaced by an aldehyde moiety and was reacted with the target enzyme, and the resulting covalent adduct was subsequently used to screen the same library of nucleophiles for those that would react with the aldehyde. This way, starting from a substrate-like peptide aldehyde, we obtained a low-lM non-peptidic, reversible inhibitor from two fragments which by themselves had no or only very low inhibitory activity (Schmidt et al., 2008) .

Interestingly, it turned out that some peptide aldehydes with an amino-acid sequence deviating from the consensus sequence embracing the cleavage sites of polyprotein substrates were surprisingly efficient as inhibitors. In particular, peptide aldehydes with P2 = Asp or Ser inhibited the main protease with a relatively low IC 50 (Al-Gharabli et al., 2006) , although it is generally agreed that the S2 specificity subsite of the enzyme has a strong preference for large hydrophobic side-chains such as Leu, Phe, or Met (Fan et al., 2004 (Fan et al., , 2005 Lai et al., 2006) . We therefore assumed that the hydrophilic P2 side-chain of these inhibitors would be oriented towards the solvent, rather than occupy the S2 pocket, and that the hydrophobic P3 residue would be binding to that pocket (Al-Gharabli et al., 2006; Schmidt et al., 2008) . This model was based on a crystal structure of a peptidyl chloromethyl ketone bound to the SARS-CoV main protease, in which this arrangement had been observed (Yang et al., 2003) . Here we show by X-ray crystallography of a number of peptide aldehyde complexes with the protein that this is not the case. Instead, the hydrophilic Ser or Asp residues in the P2 position bind to the hydrophobic S2 subsite. In order to understand these binding modes, the atomic interactions were analyzed in detail. Also, we report the inhibition kinetics of these compounds and re-evaluate the cleavage specificity of the enzyme as far as the P2 position is concerned. In addition, we determined a crystal structure of SARS-CoV M pro in complex with the aldehyde Cm-FF-H (cinnamoyl-Phe-Phe-H), which has a phenylalanine residue in the P1 position and exhibits high inhibitory activity against SARS-CoV M pro , with K i = 2.24 ± 0.58 lM. These results show that the stringent substrate specificity of the SARS-CoV M pro with respect to the P1 and P2 positions can be overruled by the highly electrophilic character of the aldehyde warhead.

Chemical synthesis of peptide aldehydes Ac-ESTLQ-H, Ac-NSFSQ-H, Ac-DSFDQ-H, and Ac-NSTSQ-H was performed employing a solid-state method (Al-Gharabli et al., 2006) . Briefly, protected glutamine aldehyde obtained by racemization-free oxidation of the corresponding amino alcohol with Dess-Martin periodinane was immobilized on a threonyl resin as oxazolidine. Following N-tert-butyl oxycarbonyl (Boc)-protection of the ring nitrogen to yield the N-protected oxazolidine linker, peptide synthesis was performed on the resin.

Synthesis of Cm-FF-H was carried out by amidation of Boc-L-phenylalanine with L-phenylalanine methyl ester followed by deprotection of the Boc group and acylation of the corresponding product with cinnamoyl chloride to provide N-cinnamoyl-L-phenylalanyl-L-phenylalanine methyl ester. Cm-FF-H was then obtained by reduction of the ester with NaBH 4 /CaCl 2 and alcohol oxidation with IBX.

The recombinant production and purification of SARS-CoV M pro with authentic N and C termini were performed as described previously (Xue et al., 2007; Verschueren et al., 2008) . The substrate Dabcyl-KTSAVLQ;SGFRKME-(Edans)-amide (95% purity; Biosyntan GmbH, Berlin, Germany), which contains a main-protease cleavage site (indicated by the arrow), was used as the substrate in the fluorescence resonance energy transfer (FRET)-based cleavage assay. The substrate stock was prepared by dissolving 1 mM of the peptide in DMSO. The dequenching of the Dabcyl fluorescence due to the cleavage of the substrate as catalyzed by the SARS-CoV M pro was monitored at 490 nm with excitation at 340 nm, using a Cary Eclipse fluorescence spectrophotometer. The experiments were performed in the buffer consisting of 20 mM Tris-HCl (pH 7.3), 100 mM NaCl, and 1 mM EDTA. 

. In order to estimate the population of catalytically active M pro dimers, the respective monomer and dimer concentrations were calculated according to Graziano et al. (2006) . The dimerization of the SARS-CoV M pro follows the scheme M þ MD The equilibrium dissociation constant, K D , is defined by

where [M] and [D] are the molar concentrations of the monomer and dimer, respectively. The total protein concentration [M T ] expressed in terms of molar monomer equivalents is

substituting the above equation into the expression for the K D yields

Solving the above quadratic equation for [M] gives

It then follows that

Aldehyde stock solutions were prepared by dissolving the compounds in DMSO at 1 mM. All aldehydes form covalent bonds between the aldehyde (-CHO) group of the inhibitor and the sulfhydryl (-SH) group of Cys145 of SARS-CoV M pro . The binding of these inhibitors is reversible (Schmidt et al., 2008) , although it is strong. Therefore, they were treated as reversible tight-binding inhibitors. For the measurement of the inhibition constant K i , SARS-CoV M pro was incubated with the aldehyde inhibitor in reaction buffer at room temperature for 30 min. Then the FRET peptide substrate, Dabcyl-KTSAVLQ;SGFRKME-Edans, was added and the reaction velocity was calculated according to substrate cleavage. (Ascenzi et al., 1987; Copeland, 2000) . Values of the intrinsic inhibition constant (K i ) for aldehyde binding to SARS-CoV M pro were calculated according to the equation K app (Ascenzi et al., 1987; Copeland, 2000) .

Twenty peptide substrates harboring alternative amino-acid residues in P2 (=X; SWTSAVXQ;SGFRKWA) were purchased from GL Biochemistry Ltd. (Shanghai, China). These peptides correspond to the N-terminal autocleavage site of the SARS-CoV Mpro, with the exception of the P7 Ile which had been replaced by Trp, and the P6 0 Met which had also been replaced by Trp. All peptide substrate stocks were dissolved in DMSO at 10 mM. To determine the k cat /K m for the substrate, 0.2 mM substrate peptide was incubated with 10 lM SARS-CoV Mpro in 40 mM Tris-HCl buffer, pH 7.3. Aliquots of reactions were removed at different times, stopped by the addition of 1% trichloroacetic acid. Separation of products and substrate was carried out using a reverse-phase (RP) HPLC and a linear gradient (1-90%) of acetonitrile in 0.1% trifluoroacetic acid. The absorbance was determined at 280 nm, and peak areas were calculated by integration. The (k cat /K m ) app ratio was determined by plotting the substrate peak area using the equation lnPA = C À (k cat / K m ) app  C E  t, where PA is the peak area of the substrate peptide, C E is the total concentration of SARS-CoV M pro , and C is an experimental constant (Fan et al., 2004 (Fan et al., , 2005 .

SARS-CoV M pro with authentic chain termini was concentrated to 10 mg/ml and crystallized by vapor diffusion using sitting drops (Xue et al., 2007) . The crystals grew overnight at 20°C by equilibration against a reservoir containing 6-8% polyethylene glycol (PEG) 6000, 0.1 M MES (pH 6.0), 3% 2-methyl-2,4-pentanediol (MPD), and 3% DMSO. All aldehydes were dissolved in 8% PEG 6000, 0.1 M MES (pH 6.0), 3% MPD, and 10% DMSO to a concentration of 10 mM. Crystals of the aldehyde complexes of SARS-CoV M pro were obtained either by adding a 4-ll aliquot of aldehyde solution to the drop and soaking of the crystals for 12 h, or by incubating the enzyme for 2 h at 20°C with a 7-fold excess of the aldehyde solution and subsequent cocrystallizing at 20°C against a reservoir containing 8% PEG 6000, 0.1 M MES (pH 6.0), 3% MPD, and 3% DMSO. In the latter case, nucleation was initiated by microseeding using crushed monoclinic (space group C2) crystals of the SARS-CoV M pro . Prior to diffraction data collection, the aldehyde complex crystals obtained from soaking and cocrystallization were transferred for a few seconds to a cryoprotectant solution containing the crystallization ingredients and 20% MPD.

All diffraction data were collected at 100 K at the Joint EMBL/ University of Hamburg/University of Lübeck synchrotron beamline X13 at DESY (Hamburg, Germany), using a 165-mm MAR CCD detector (Mar Research, Hamburg, Germany), or at synchrotron beamline BL14.1 at BESSY (Berlin, Germany), using an MX225 CCD detector (Rayonics, Evanston, IL). Statistics of data collection, processing, and refinement are summarized in Table 1 . Data were processed with MOSFLM (Leslie, 1992) , and scaled using the SCALA program from the CCP4 suite (Collaborative Computational Project, 1994; Evans, 2006) .

Structure elucidation and refinement were carried out using the CCP4 suite of programs (Collaborative Computational Project, 1994; Potterton et al., 2003) . Crystal structures were determined by molecular replacement, using the original X-ray structure of the SARS-CoV M pro with authentic chain termini (PDB ID: 2H2Z) (Xue et al., 2007) as the initial model. REFMAC (Murshudov et al., 1997) was employed for structure refinement. The computer graphics program Coot (Emsley and Cowtan, 2004) was used for interpretation of electron density maps and model building. Gap volumes between inhibitor and protein in the subsites of the enzyme were calculated using UCSF-Chimera (Pettersen et al., 2004) . The molecular graphics package PyMOL (DeLano, 2002) was used to generate the figures.

The kinetic parameters of SARS-CoV M pro with authentic chain termini were determined using the FRET substrate Dabcyl-KTSAVLQ;SGFRKME-Edans amide. The K m value for this substrate was determined as 24.5 lM and the apparent k cat /K m value was 2359 M À1 Á19 s À1 . However, the dissociation of the catalytically active SARS-CoV M pro dimer into inactive monomers (Fan et al., 2004) has to be taken into account. A wide range of K D values has been reported for M pro dimer dissociation (see Grum-Tokars et al. (2008) for an overview). For the enzyme with authentic chain termini, the latter authors reported a K D of 0.25 to 1.0 lM. As the M pro dimer tends to be stabilized by the presence of substrate (Cheng et al., 2010) , we used the lower limit of this range for estimation of the necessary corrections and obtained a k cat /K m value of 7863 M À1 s À1 . The K m value determined in our study is of the same magnitude as other values derived from data not corrected for dissociation, which are typically in the range from 10 to 50 lM (Grum-Tokars et al., 2008) . This suggests that much higher values reported previously (e.g., Verschueren et al., 2008) should probably be reconsidered.

Initially, four pentapeptide substrate-analogous aldehydes containing the canonical P1 residue, glutamine, namely Ac-ESTLQ-H, Ac-NSTSQ-H, Ac-DSFDQ-H, and Ac-NSFSQ-H, were tested for inhibition of SARS-CoV M pro . Overall, the four analogues were found to be reversible tight-binding inhibitors. Harboring the canonical P2 Leu, aldehyde Ac-ESTLQ-H exhibited inhibition with a relatively low K i of 8.27 ± 1.52 lM. Aldehydes Ac-NSTSQ-H, Ac-DSFDQ-H, and Ac-NSFSQ-H, all with a non-canonical P2 residue, exhibited moderate inhibition with K i values of 40.98 ± 2.63, 41.24 ± 2.25, and 72.73 ± 3.60 lM. Surprisingly, aldehyde CmFF-H, carrying a cinnamoyl group in the P3 and a Phe residue in the P1 position, had an even higher inhibitory activity against SARS-CoV M pro than the four pentapeptide aldehydes, with a K i of 2.24 ± 0.58 lM.

The aldehydes Ac-ESTLQ-H, Ac-NSTSQ-H, Ac-DSFDQ-H, Ac-NSFSQ-H, and Cm-FF-H were separately soaked into crystals of SARS-CoV Mpro. The crystals were all of space group C2, which is often observed for SARS-CoV M pro (Lee et al., 2005; Xue et al., 2007; Verschueren et al., 2008) . These crystals contain one SARS-CoV M pro monomer per asymmetric unit and the dimer (which is the enzymatically active species) is formed through the symmetry of the crystal. The four pentapeptide aldehydes Ac-ESTLQ-H, Ac-NSTSQ-H, Ac-DSFDQ-H, and Ac-NSFSQ-H are bound in extended conformations in the S6-S1 specificity subsites of SARS-CoV M pro . Cm-FF-H occupies sites S3-S1. Remarkably, the P1 phenylalanine side chain of this inhibitor is bound deeply in the S1 pocket, which is generally considered to be specific for glutamine. 2F o À F c electron density maps of these aldehyde inhibitors are shown in Fig. 1 . In all complexes, continuous electron density between the aldehyde carbonyl C-atom of the inhibitor and Cys145-Sc of SARS-CoV M pro (Fig. 1) indicates the formation of a thiohemiacetal, as a result of the nucleophilic attack of the catalytic cysteine onto the C-terminal aldehyde of the inhibitor. The main-chain conformations of SARS-CoV M pro in the five complexes are basically identical, with overall root mean-square deviations (RMSD) of 0.16 À 0.36 Å for Ca atoms. In addition to the crystal soaking experiments, we also tried to cocrystallize SARS-CoV M pro with the aldehydes. The crystals obtained were predominantly of space group P2 1 , with the exception of the complex with Ac-NSFSQ-H, which still displayed space group C2. However, in most of the P2 1 crystals, no electron density for the aldehyde could be observed; thus, while the presence of the inhibitors induced a change of space group in most cases, the compounds themselves were not detected in the crystals. An exception was Ac-ESTLQ-H. The crystals of its complex with the M pro diffracted to 1.89 Å, but in both enzyme monomers in the asymmetric unit, electron density could only be seen for residues P1 (Gln) and P2 (Leu) of the inhibitor (Fig. 1E and F) . Crystallographic data and refinement statistics for all six crystal structures are summarized in Table 1 .

In the SARS-CoV M pro complex with aldehyde inhibitor Ac-NSFSQ-H, the thiohemiacetal adopts alternative configurations ( Fig. 2A and B) . This is likely the result of nucleophilic attack by Cys145 onto the planar carbonyl from either side. In the (R)-isomer, the thiohemiacetal oxygen forms a hydrogen bond with the imidazole of His41 (3.30 Å) of the catalytic dyad (Fig. 2C) , whereas in the (S)-isomer, the oxygen atom points away from His41 and into the oxyanion hole, forming H-bonds with the main-chain amides of Gly143 (2.71 Å) and Cys 145 (3.15 Å) (Fig. 2D) . Previous X-ray and NMR studies of aldehyde binding to various cysteine and serine proteases demonstrated that one or both configurations were present in the structures of the complexes (Delbaere and Brayer, 1985; Webber et al., 1998; Robin et al., 2009) .

The S2 subsite is a large pocket lined by residues His41, Met49, Cys145, His164, Met165, Asp187, Arg188mc (mc: contribution from main-chain atoms only), and Gln 189 (Fig. 1) . At the cleavage sites of the SARS-CoV M pro in the viral polyproteins, the P2 position is mostly found to be occupied by Leu or Phe. Accordingly, almost all inhibitors designed for the enzyme carry a large hydrophobic group in the P2 position. We were therefore very surprised to find the Ser or Asp side chains in the P2 position of the aldehydes bound to this hydrophobic site (Fig. 1) . In the SARS-CoV M pro complexes with Ac-NSTSQ-H or Ac-NSFSQ-H, there is no hydrogen-bonding interaction between the side chain of P2-Ser and residues of the S2 subsite. As the Ser side chain is small, there is no steric barrier for Ser binding in the large S2 pocket. The serine does not quite reach the ''bottom'' of the S2 pocket, but the distance between its Oc atom and the carbonyl oxygen of residue Gln192 is 4.0 Å, leaving no space for a water molecule to locate in between. Accordingly, we did not find any electron density suggesting the presence of a well-ordered water in the S2 pocket, although disordered water is probably present in the pocket on both sides of the serine side-chain, where free spaces with a total volume of >70 Å 3 are present between the substrate and the protein. It should be noted that the side-chain of Gln189, which contributes to one wall of the S2 pocket, is quite flexible and adopts different conformations in the various inhibitor complexes, allowing or preventing the access of water to the pocket.

In case of the SARS-CoV M pro complex with Ac-DSFDQ-H, the P2 Asp side-chain is situated between the side chains of Met49 and Met165, which form opposite walls of the S2 pocket. Most interestingly, its carboxylate oxygens undergo close interactions with the sulfur atoms of the two methionines. The distances (r S. . .O ) between the Sd atoms of Met49 and Met165 on the one hand and the Od2 and Od1 atoms, respectively, of P2-Asp on the other are 3.36 and 3.45 Å (Fig. 3) . Then the relative distances (d S. between the methionine sulfur atoms and main-chain carbonyl oxygens or carboxylate side-chains have been detected previously in the hydrophobic cores of proteins and were proposed to stabilize the protein fold (Pal and Chakrabarti, 2001) . It has also been suggested that S. . .O interactions should be taken into account in protein engineering studies (Iwaoka et al., 2002; Pal and Chakrabarti, 2001 ), but to the best of our knowledge, we provide here the first description of a methionine-carboxylate interaction in a protein-ligand complex. The unexpected finding of Ser and Asp binding in the S2 subsite constitutes a deviation from the dogma that peptide inhibitors of proteases should contain amino-acid residues corresponding to the sequence specificity of the target enzyme.

The unexpected observation of the P2-Asp residue of aldehyde Ac-DSFDQ-H binding to the hydrophobic S2 pocket prompted us to re-determine the cleavage specificity of the SARS-CoV M pro with respect to the P2 position. Within the context of the peptide substrate SWTSAVXQ;SGFRKWA, all 20 proteinogenic amino acids were tested in position P2 (=X). In the HPLC-based assay, the protease concentration was kept constant at 0.5 lM. The canonical substrate with Leu in the P2 position was completely cleaved within two minutes and was therefore used as the standard to compare the proteolytic activities with other substrates. The results are listed in Fig. 4 . Our study confirmed that the most favored residues in the P2 position of the substrate (after Leu) are the hydrophobic amino acids Phe, Met, Val, and Ile. Substrates with polar amino acids in P2 are poorly cleaved. The cleavage efficiency for the substrate with P2 = Ser was 160-times lower than for the best substrate (P2 = Leu); furthermore, the peptide with P2 = Asp was not cleaved at all after incubation with SARS-CoV M pro for 16 h.

The S1 subsite is lined by residues Phe140, Asn142, Gly143, Ser144, Cys145, His163, Glu166, and His172. From the substrate specificity results (Fan et al., 2004 (Fan et al., , 2005 Lai et al., 2006) , P1 has to be Gln. Moreover, in all SARS-CoV M pro complex structures described so far (Lee et al., 2005; Xue et al., 2007; Yang et al., 2005 Yang et al., , 2003 Yin et al., 2007) , the S1 subsite is occupied by a polar group, i.e. the side chain of Gln or a five-membered lactam (used as a glutamine surrogate). Shie et al. synthesized a series of a,b-unsaturated esters containing both P1 and P2 phenylalanine residues which showed modest inhibitory activity against the SARS-CoV M pro (IC 50 = 11-39 lM) (Shie et al., 2005) . The a,b-unsaturated ethyl ester of 4-(dimethylamino)cinnamoyl-Phe-Phe was reported to be a potent inhibitor, with an inhibition constant of 0.52 lM (Shie et al., 2005) . The computer model for the complex of SARS-CoV M pro with the compound proposed that the P1 and P2 phenyl groups occupy the S2 and S3 pockets, respectively (Shie et al., 2005) . However, in the crystal structure of the complex with the aldehyde Cm-FF-H that we describe in this communication, the side chain of P1-Phe is inserted into the S1 subsite (Fig. 1G) . In order to understand this unexpected deviation from the commonly observed specificity, we compared the interaction with the S1 subsite observed in our structures of the complexes with Cm-FF-H and Ac-ESTLQ-H, the latter of which corresponds exactly to the described cleavage specificity of the M pro . The predominant S1 specificity of the enzyme for Gln is determined primarily by His163. In the complex formed with Ac-ESTLQ-H, Ne2 of His163 donates a 2.83-Å hydrogen bond to the side-chain carbonyl oxygen (Oe1) of P1-Gln (Fig. 1) . In the complex with Cm-FF-H, the side chain of Asn142 undergoes an 83.8°rotation about its v1 torsion angle compared to the conformation in the complex with Ac-ESTLQ-H (Fig. 5) , leading to an opening of the S1 subsite towards the solvent. A similar movement of Asn142 has been observed in the crystal structure of SARS-CoV M pro in complex with an inhibitor called ''N1'' (PDB ID: 1WOF), where both the regular conformation of this residue and the one observed in our structure were found . The P1-benzyl group of Cm-FF-H binds to the S1 subsite by making hydrophobic interactions with Phe140, Leu141, Asn142, and the P3 cinnamoyl group of Cm-FF-H. As the S1 subsite specificity does not seem to be as stringent as previously thought. It may well be possible that other chemical groups can be accommodated here, thereby expanding the opportunities for designing and synthesizing efficient inhibitors of the SARS-CoV M pro .

In this study, we report six crystal structures of SARS-Cov M pro in complex with peptide aldehydes and the inhibition kinetics of these compounds. The crystal structures reveal the mechanism of aldehyde binding to the active site of SARS-Cov M pro and provide detailed information on the atomic interactions. Since Asp or Ser were found to bind to the hydrophobic S2 subsite, and Phe was located in the hydrophilic S1 subsite of SARS-Cov M pro , we conclude that the stringent substrate specificity of the SARS-CoV M pro with respect to the P1 and P2 positions can be overruled by the highly electrophilic character of the aldehyde warhead. This constitutes a deviation from the dogma that peptidic protease inhibitors should comprise an amino-acid sequence corresponding to the cleavage specificity of the target enzyme. The observed nonbonded interaction of the carboxylate oxygen atoms of an aspartate residue of an inhibitor with the thioether moieties of two methionines forming part of an hydrophobic pocket of the target protein is remarkable and suggests that such S. . .O interactions should be added to the repertoire of computer-aided drug design.

Atomic coordinates and structure factors have been deposited in the RCSB Protein Data Bank under ID codes 3SN8, 3SNA, 3SNB, 3SNC, 3SND, and 3SNE (see Table 1 ).

An efficient method for the synthesis of peptide aldehyde libraries employed in the discovery of reversible SARS coronavirus main protease (SARS-CoV M pro ) inhibitors

Coronavirus main proteinase (3CL pro ) structure: basis for design of anti-SARS drugs

Enzyme competitive-inhibitiongraphical determination of K i and presentation of data in comparative studies

Mutation of Glu-166 blocks the substrateinduced dimerization of SARS coronavirus main protease

The CCP4 suite: programs for protein crystallography

Enzymes: A Practical Introduction to Structure, Mechanism, and Data Analysis

The PyMOL Molecular Graphics System Delano Scientific

The 1.8 A structure of the complex between chymostatin and Streptomyces griseus protease A. A model for serine protease catalytic tetrahedral intermediates

Coot: model-building tools for molecular graphics

Scaling and assessment of data quality

Biosynthesis, purification, and substrate specificity of severe acute respiratory syndrome coronavirus 3C-like proteinase

The substrate specificity of SARS coronavirus 3C-like proteinase

SARS CoV main proteinase: the monomer-dimer equilibrium dissociation constant

Evaluating the 3C-like protease activity of SARS-Coronavirus: recommendations for standardized assays for drug discovery

Statistical and theoretical investigations on the directionality of nonbonded S. . .O interactions. Implications for molecular design and protein engineering

Quaternary structure, substrate selectivity and inhibitor design for SARS 3C-like proteinase

Crystal structures of the main peptidase from the SARS coronavirus inhibited by a substrate-like aza-peptide epoxide

Recent changes to the MOSFLM package for processing film and image plate data

Refinement of macromolecular structures by the maximum-likelihood method

Non-hydrogen bond interactions involving the methionine sulfur atom

UCSF Chimera -a visualization system for exploratory research and analysis

A graphical user interface to the CCP4 program suite

Structure of West Nile virus NS3 protease: ligand stabilization of the catalytic conformation

Sensitized detection of inhibitory fragments and iterative development of non-peptidic protease inhibitors by dynamic ligation screening

Inhibition of the severe acute respiratory syndrome 3CL protease by peptidomimetic alpha, beta-unsaturated esters

Recent advances in targeting viral proteases for the discovery of novel antivirals

PHdependent conformational flexibility of the SARS-CoV main proteinase (M pro ) dimer: molecular dynamics simulations and multiple X-ray structure analyses

A structural view of the inactivation of the SARS coronavirus main proteinase by benzotriazole esters

Tripeptide aldehyde inhibitors of human rhinovirus 3C protease: design, synthesis, biological evaluation, and cocrystal structure solution of P1 glutamine isosteric replacements

Production of authentic SARS-CoV M pro with enhanced activity: application as a novel tag-cleavage endopeptidase for protein overproduction

Design of wide-spectrum inhibitors targeting coronavirus main proteases

The crystal structures of severe acute respiratory syndrome virus main protease and its complex with an inhibitor

A mechanistic view of enzyme inhibition and peptide hydrolysis in the active site of the SARS-CoV 3C-like peptidase

We thank Sebastian Schwarz for technical assistance during the re-determination of the cleavage rates of SARS-CoV M pro with respect to the P2 position of the substrate, and Stefan Anemüller, Qingjun Ma, Jeroen R. Mesters as well as Jiajie Zhang for discussion. This work was initially supported by the Deutsche Forschungsgemeinschaft (Grants Hi611/4, Ra895/2), the Sino-German Center for the Promotion of Research, Beijing, and the Sino-European Project on SARS Diagnostics and Antivirals (SEPSDA) of the European Commission (contract number SP22-CT-2004-003831), and subsequently, by the SILVER project of the European Commission (contract number HEALTH-F3-2010-260644). RH and JR acknowledge continuous support by the Fonds der Chemischen Industrie. RH is also supported by a Chinese Academy of Sciences Visiting Professorship for Senior International Scientists, Grant No. 2010T1S6.