Revisiting catalysis by chymotrypsin family serine proteases using peptide substrates and inhibitors with unnatural main chains

Identifying signatures of proteolytic stability and monomeric propensity in O-glycosylated insulin using molecular simulation

Insulin has been commonly adopted as a peptide drug to treat diabetes as it facilitates the uptake of glucose from the blood. The development of oral insulin remains elusive over decades owing to its susceptibility to the enzymes in the gastrointestinal tract and poor permeability through the intestinal epithelium upon dimerization. Recent experimental studies have revealed that certain O-linked glycosylation patterns could enhance insulin's proteolytic stability and reduce its dimerization propensity, but understanding such phenomena at the molecular level is still difficult. To address this challenge, we proposed and tested several structural determinants that could potentially influence insulin's proteolytic stability and dimerization propensity. We used these metrics to assess the properties of interest from [Formula: see text] aggregate molecular dynamics of each of 12 targeted insulin glyco-variants from multiple wild-type crystal structures. We found that glycan-involved hydrogen bonds and glycan-dimer occlusion were useful metrics predicting the proteolytic stability and dimerization propensity of insulin, respectively, as was in part the solvent-accessible surface area of proteolytic sites. However, other plausible metrics were not generally predictive. This work helps better explain how O-linked glycosylation influences the proteolytic stability and monomeric propensity of insulin, illuminating a path towards rational molecular design of insulin glycoforms.

The role of N-terminal heterocycles in hydrogen bonding to α-chymotrypsin

A series of dipeptide aldehydes containing different N-terminal heterocycles was prepared and assayed in vitro against α-chymotrypsin to ascertain the importance of the heterocycle in maintaining a β-strand geometry while also providing a hydrogen bond donor equivalent to the backbone amide nitrogen of the surrogate amino acid. The dipeptide containing a pyrrole constraint (10) was the most potent inhibitor, with >30-fold improved activity over dipeptides which lacked a nitrogen hydrogen bond donor (namely thiophene 11, furan 12 and pyridine 13). Molecular docking studies of 10 bound to α-chymotrypsin demonstrates a hydrogen bond between the pyrrole nitrogen donor and the backbone carbonyl of Gly₂₁₆ located in the S₃ pocket which is proposed to be critical for overall binding.

The rhomboid protease GlpG has weak interaction energies in its active site hydrogen bond network

Rhomboid proteases are membrane-integrated enzymes that hydrolyze peptide bonds in the transmembrane domains of protein substrates. Gaffney and Hong experimentally determine interaction energies between active site residues to reveal weak coupling, which may explain the slow proteolysis mediated by GlpG.

Intramembrane rhomboid proteases are of particular interest because of their function to hydrolyze a peptide bond of a substrate buried in the membrane. Crystal structures of the bacterial rhomboid protease GlpG have revealed a catalytic dyad (Ser201-His254) and oxyanion hole (His150/Asn154/the backbone amide of Ser201) surrounded by the protein matrix and contacting a narrow water channel. Although multiple crystal structures have been solved, the catalytic mechanism of GlpG is not completely understood. Because it is a serine protease, hydrogen bonding interactions between the active site residues are thought to play a critical role in the catalytic cycle. Here, we dissect the interaction energies among the active site residues His254, Ser201, and Asn154 of Escherichia coli GlpG, which form a hydrogen bonding network. We combine double mutant cycle analysis with stability measurements using steric trapping. In mild detergent, the active site residues are weakly coupled with interaction energies (ΔΔG_Inter) of ‒1.4 kcal/mol between His254 and Ser201 and ‒0.2 kcal/mol between Ser201 and Asn154. Further, by analyzing the propagation of single mutations of the active site residues, we find that these residues are important not only for function but also for the folding cooperativity of GlpG. The weak interaction between Ser and His in the catalytic dyad may partly explain the unusually slow proteolysis by GlpG compared with other canonical serine proteases. Our result suggests that the weak hydrogen bonds in the active site are sufficient to carry out the proteolytic function of rhomboid proteases.

KLK14 interactions with HAI-1 and HAI-2 serine protease inhibitors: A molecular dynamics and relative free-energy calculations study

Kallikrein 14 (KLK14) is a serine protease linked to several pathologies including prostate cancer and positively correlates with Gleason score. Though KLK14 functioning in cancer is poorly understood, it has been implicated in HGF/Met signaling, given that KLK14 proteolytically inhibits HGF activator-inhibitor 1 (HAI-1), which strongly inhibits pro-HGF activators, thereby contributing to tumor progression. In this work, KLK14 binding to either hepatocyte growth factor activator inhibitor type-1 (HAI-1) or type-2 (HAI-2) was essayed using homology modeling, molecular dynamic simulations and free-energy calculations through MM/PBSA and MM/GBSA. KLK14 was successfully modeled. Calculated free energies suggested higher binding affinity for the KLK14/HAI-1 interaction than for KLK14/HAI-2. This difference in binding affinity is largely explained by the higher stability of the hydrogen-bond networks in KLK14/HAI-1 along the simulation trajectory. A key arginine residue in both HAI-1 and HAI-2 is responsible for their interaction with the S1 pocket in KLK14. Additionally, MM/GBSA free-energy decomposition postulates that KLK14 Asp174 and Trp196 are hotspots for binding HAI-1 and HAI-2.

Amide Rotation Hindrance Predicts Proteolytic Resistance of Cystine-Knot Peptides

Cystine-knot peptides have remarkable stability against protease degradation and are attractive scaffolds for peptide-based therapeutic and diagnostic agents. In this work, by studying the hydrolysis reaction of a cystine-knot inhibitor MCTI-A and its variants with ab initio QM/MM molecular dynamics simulations, we have elucidated an amide rotation hindrance mechanism for proteolysis resistance: The proteolysis of MCTI-A is retarded due to the higher free energy cost during the rotation of NH group around scissile peptide bond at the tetrahedral intermediate of acylation, and covalent constraint provided by disulfide bonds is the key factor to hinder this rotation. A nearly linear correlation has been revealed between free energy barriers of the peptide hydrolysis reaction and the amide rotation free energy changes at the protease-peptide Michaelis complex state. This suggests that amide rotation hindrance could be one useful feature to estimate peptide proteolysis stability.

Insights into the serine protease mechanism based on structural observations of the conversion of a peptidyl serine protease inhibitor to a substrate

BACKGROUND

Serine proteases are one of the most studied group of enzymes. Despite the extensive mechanistic studies, some crucial details remain controversial, for example, how the cleaved product is released in the catalysis reaction. A cyclic peptidyl inhibitor (CSWRGLENHRMC, upain-1) of a serine protease, urokinase-type plasminogen activator (uPA), was found to become a slow substrate and cleaved slowly upon the replacement of single residue (W3A).

METHODS

By taking advantage of the unique property of this peptide, we report the high-resolution structures of uPA in complex with upain-1-W3A peptide at four different pH values by X-ray crystallography.

RESULTS

In the structures obtained at low pH (pH4.6 and 5.5), the cyclic peptide upain-1-W3A was found to be intact and remained in the active site of uPA. At 7.4, the scissile bond of the peptide was found cleaved, showing that the peptide became a uPA substrate. At pH9.0, the C-terminal part of the substrate was no longer visible, and only the P1 residue occupying the S1 pocket was identified.

CONCLUSIONS

The analysis of these structures provides explanations why the upain-1-W3A is a slow substrate. In addition, we clearly identified the cleaved fragments of the peptide at both sides of the scissile bond in the active site of the enzyme, showing a slow release of the cleaved peptide.

GENERAL SIGNIFICANCE

This work indicates that the quick release of the cleaved P' fragment after the first step of hydrolysis may not always be needed for the second hydrolysis.

Structure of a proteolytically resistant analogue of (NLys)5SFTI-1 in complex with trypsin: evidence for the direct participation of the Ser214 carbonyl group in serine protease-mediated proteolysis

Peptide-peptoid hybrids are found to be potent inhibitors of serine proteases. These engineered peptidomimetics benefit from both types of units of the biopolymeric structure: the natural inhibitor part serves as a good binding template, while the P1-positioned peptoid component provides complete resistance towards proteolysis. In this report, the mechanism of proteolytic resistance of a P1 peptoid-containing analogue is postulated based on the crystal structure of the (NLys)(5)-modified sunflower trypsin inhibitor SFTI-1 in complex with bovine trypsin solved at 1.29 Å resolution. The structural differences between the (NLys)(5)SFTI-1-trypsin complex and the native SFTI-1-trypsin complex are surprisingly small and reveal the key role of the carbonyl group of the Ser214 residue of the enzyme, which is crucial for binding of the inhibitor and plays a crucial role in proteolysis mediated by serine proteases. The incorporated NLys5 peptoid residue prevents Ser214 from forming a hydrogen bond to the P1 residue, and in turn Gln192 does not form a hydrogen bond to the carbonyl group of the P2 residue. It also increases the distance between the Ser214 carbonyl group and the Ser195 residue, thus preventing proteolysis. The hybrid inhibitor structure reported here provides insight into protein-protein interaction, which can be efficiently and selectively probed with the use of peptoids incorporated within endogenous peptide ligands.

Time-resolved botulinum neurotoxin A activity monitored using peptide-functionalized Au nanoparticle energy transfer sensors

A peptide mimicking SNAP-25 was immobilized on Au nanoparticles for the monitoring of botulinum neurotoxin light chain A activity.

Substrate phage display for protease substrate sequence characterization: bovine factor Xa as a model system

Regulatory proteases modulate proteomic dynamics with a spectrum of specificities against substrate proteins. Substrate phage display is one of the key methodologies in producing substrate sequence information in vitro. Factor Xa, a key regulatory protease in the blood coagulation system, is used as a model system to demonstrate a high-throughput procedure to quantitatively characterize substrate sequences and their susceptibilities for enzymatic cleavage. This methodology can be generalized to proteases for which the active forms (not necessarily purified forms) are available for the in vitro experiments.

Preparation and evaluation at the delta opioid receptor of a series of linear leu-enkephalin analogues obtained by systematic replacement of the amides

Leu-enkephalin analogues, in which the amide bonds were sequentially and systematically replaced either by ester or N-methyl amide bonds, were prepared using classical organic chemistry as well as solid phase peptide synthesis (SPPS). The peptidomimetics were characterized using competition binding, ERK1/2 phosphorylation, receptor internalization, and contractility assays to evaluate their pharmacological profile over the delta opioid receptor (DOPr). The lipophilicity (LogD7.4) and plasma stability of the active analogues were also measured. Our results revealed that the last amide bond can be successfully replaced by either an ester or an N-methyl amide bond without significantly decreasing the biological activity of the corresponding analogues when compared to Leu-enkephalin. The peptidomimetics with an N-methyl amide function between residues Phe and Leu were found to be more lipophilic and more stable than Leu-enkephalin. Findings from the present study further revealed that the hydrogen-bond donor properties of the fourth amide of Leu-enkephalin are not important for its biological activity on DOPr. Our results show that the systematic replacement of amide bonds by isosteric functions represents an efficient way to design and synthesize novel peptide analogues with enhanced stability. Our findings further suggest that such a strategy can also be useful to study the biological roles of amide bonds.

Probing the role of backbone hydrogen bonds in protein-peptide interactions by amide-to-ester mutations

One of the most frequent protein-protein interaction modules in mammalian cells is the postsynaptic density 95/discs large/zonula occludens 1 (PDZ) domain, involved in scaffolding and signaling and emerging as an important drug target for several diseases. Like many other protein-protein interactions, those of the PDZ domain family involve formation of intermolecular hydrogen bonds: C-termini or internal linear motifs of proteins bind as β-strands to form an extended antiparallel β-sheet with the PDZ domain. Whereas extensive work has focused on the importance of the amino acid side chains of the protein ligand, the role of the backbone hydrogen bonds in the binding reaction is not known. Using amide-to-ester substitutions to perturb the backbone hydrogen-bonding pattern, we have systematically probed putative backbone hydrogen bonds between four different PDZ domains and peptides corresponding to natural protein ligands. Amide-to-ester mutations of the three C-terminal amides of the peptide ligand severely affected the affinity with the PDZ domain, demonstrating that hydrogen bonds contribute significantly to ligand binding (apparent changes in binding energy, ΔΔG = 1.3 to >3.8 kcal mol(-1)). This decrease in affinity was mainly due to an increase in the dissociation rate constant, but a significant decrease in the association rate constant was found for some amide-to-ester mutations suggesting that native hydrogen bonds have begun to form in the transition state of the binding reaction. This study provides a general framework for studying the role of backbone hydrogen bonds in protein-peptide interactions and for the first time specifically addresses these for PDZ domain-peptide interactions.

Assembly of a concentric Förster resonance energy transfer relay on a quantum dot scaffold: characterization and application to multiplexed protease sensing

Semiconductor nanocrystals, or quantum dots (QDs), are one of the most widely utilized nanomaterials for biological applications. Their cumulative physicochemical and optical properties are both unique among nanomaterials and highly advantageous. In particular, Förster resonance energy transfer (FRET) has been widely utilized as a spectroscopic tool with QDs, whether for characterizing QD bioconjugates as a "molecular ruler" or for modulating QD luminescence "on" and "off" in biosensing configurations. Here, we investigate the assembly and utility of a new "concentric" FRET relay that comprises a central QD conjugated with multiple copies of two different peptides, each labeled with one of two fluorescent dyes, Alexa Fluor 555 (A555) or Alexa Fluor 647 (A647). Energy transfer occurs from the QD to the A555 (FRET(1)) then to the A647 (FRET(2)) and, to a lesser extent, directly from the QD to the A647 (FRET(3)). We show that such an arrangement can provide insight into the interfacial distribution of peptides assembled to the QD and can further be utilized for sensing proteolytic activity. In the latter, progress curves for digestion of the assembled peptides by two prototypical proteases, trypsin and chymotrypsin, were measured from the relative QD, A555 and A647 PL contributions, and used to extract Michaelis-Menten kinetic parameters. We further show that the concentric FRET relay, as a single nanoparticle vector, can track the tryptic activation of a proenzyme, chymotrypsinogen, to active chymotrypsin. The concentric FRET relay is thus a potentially powerful tool for the characterization of QD bioconjugates and multiplexed sensing of coupled biological activity.

Proteolytic activity at quantum dot-conjugates: kinetic analysis reveals enhanced enzyme activity and localized interfacial "hopping"

Recent studies show that polyvalent, ligand-modified nanoparticles provide significantly enhanced binding characteristics compared to isolated ligands. Here, we assess the ability of substrate-modified nanoparticles to provide enhanced enzymatic activity. Energy transfer assays allowed quantitative, real-time measurement of proteolytic digestion at polyvalent quantum dot-peptide conjugates. Enzymatic progress curves were analyzed using an integrated Michaelis-Menten (MM) formalism, revealing mechanistic details, including deviations from classic MM-behavior. A "hopping" mode of proteolysis at the nanoparticle was identified, confirming enhanced activity.

Characterization and inhibition of norovirus proteases of genogroups I and II using a fluorescence resonance energy transfer assay

Noroviruses are the major cause of food- or water-borne gastroenteritis outbreaks in humans. The norovirus protease that cleaves a large viral polyprotein to nonstructural proteins is essential for virus replication and an attractive target for antiviral drug development. Noroviruses show high genetic diversity with at least five genogroups, GI-GV, of which GI and GII are responsible for the majority of norovirus infections in humans. We cloned and expressed proteases of Norwalk virus (GI) and MD145 virus (GII) and characterized the enzymatic activities with fluorescence resonance energy transfer substrates. We demonstrated that the GI and GII proteases cleaved the substrates derived from the naturally occurring cleavage site in the open reading frame (ORF) 1 of G1 norovirus with similar efficiency, and that enzymatic activity of both proteases was inhibited by commercial protease inhibitors including chymostatin. The interaction of chymostatin to Norwalk virus protease was validated by nuclear magnetic resonance (NMR) spectroscopy.

Mechanism of heparin acceleration of tissue inhibitor of metalloproteases-1 (TIMP-1) degradation by the human neutrophil elastase

Heparin has been shown to regulate human neutrophil elastase (HNE) activity. We have assessed the regulatory effect of heparin on Tissue Inhibitor of Metalloproteases-1 [TIMP-1] hydrolysis by HNE employing the recombinant form of TIMP-1 and correlated FRET-peptides comprising the TIMP-1 cleavage site. Heparin accelerates 2.5-fold TIMP-1 hydrolysis by HNE. The kinetic parameters of this reaction were monitored with the aid of a FRET-peptide substrate that mimics the TIMP-1 cleavage site in pre-steady-state conditionsby using a stopped-flow fluorescence system. The hydrolysis of the FRET-peptide substrate by HNE exhibits a pre-steady-state burst phase followed by a linear, steady-state pseudo-first-order reaction. The HNE acylation step (k₂ = 21±1 s⁻¹) was much higher than the HNE deacylation step (k₃ = 0.57±0.05 s⁻¹). The presence of heparin induces a dramatic effect in the pre-steady-state behavior of HNE. Heparin induces transient lag phase kinetics in HNE cleavage of the FRET-peptide substrate. The pre-steady-state analysis revealed that heparin affects all steps of the reaction through enhancing the ES complex concentration, increasing k₁ 2.4-fold and reducing k₋₁ 3.1-fold. Heparin also promotes a 7.8-fold decrease in the k₂ value, whereas the k₃ value in the presence of heparin was increased 58-fold. These results clearly show that heparin binding accelerates deacylation and slows down acylation. Heparin shifts the HNE pH activity profile to the right, allowing HNE to be active at alkaline pH. Molecular docking and kinetic analysis suggest that heparin induces conformational changes in HNE structure. Here, we are showing for the first time that heparin is able to accelerate the hydrolysis of TIMP-1 by HNE. The degradation of TIMP-1is associated to important physiopathological states involving excessive activation of MMPs.

Aminopeptidase fingerprints, an integrated approach for identification of good substrates and optimal inhibitors

Aminopeptidases process the N-terminal amino acids of target substrates by sequential cleavage of one residue at a time. They are found in all cell compartments of prokaryotes and eukaryotes, being implicated in the major proteolytic events of cell survival, defense, growth, and development. We present a new approach for the fast and reliable evaluation of the substrate specificity of individual aminopeptidases. Using solid phase chemistry with the 7-amino-4-carbamoylmethylcoumarin fluorophore, we have synthesized a library of 61 individual natural and unnatural amino acids substrates, chosen to cover a broad spectrum of the possible interactions in the S1 pocket of this type of protease. As proof of concept, we determined the substrate specificity of human, pig, and rat orthologs of aminopeptidase N (CD13), a highly conserved cell surface protease that inactivates enkephalins and other bioactive peptides. Our data reveal a large and hydrophobic character for the S1 pocket of aminopeptidase N that is conserved with aminopeptidase Ns. Our approach, which can be applied in principle to all aminopeptidases, yields useful information for the design of specific inhibitors, and more importantly, reveals a relationship between the kinetics of substrate hydrolysis and the kinetics of enzyme inhibition.

Depsipeptide synthesis

Naturally occurring cyclic depsipeptides, peptides that contain one or more ester bonds in addition to the amide bonds, have emerged as an important source of pharmacologically active compounds or promising lead structures for the development of novel synthetically derived drugs. This class of natural products has been found in many organisms, such as fungi, bacteria, and marine organisms. It is very well known that cyclic depsipeptides and their derivatives exhibit a diverse spectrum of biological activities, including insecticidal, antiviral, antimicrobial, antitumor, tumor-promotive, anti-inflammatory, and immunosuppressive actions. However, they have shown the greatest therapeutic potential as anticancer and particularly antimicrobial agents. Difficulties associated with isolation and purification of larger quantities of this class of natural products and, particularly, unlimited access to their synthetic analogs significantly hampered cyclic depsipeptides exploitation as lead compounds for development of new drugs. As an alternative, total solution or solid-phase peptide synthesis of these important natural products and combinatorial chemistry approaches can be employed to elucidate structure-activity relationships and to find new potent compounds of this class. In this chapter, methods for formation of depsipeptide ester bonds, hydroxyl group protection, and solid-phase reaction monitoring are described.

Functional and structural roles of the Cys14-Cys38 disulfide of bovine pancreatic trypsin inhibitor

The disulfide bond between Cys14 and Cys38 of bovine pancreatic trypsin inhibitor lies on the surface of the inhibitor and forms part of the protease-binding region. The functional properties of three variants lacking this disulfide, with one or both of the Cys residues replaced with Ser, were examined, and X-ray crystal structures of the complexes with bovine trypsin were determined and refined to the 1.58-A resolution limit. The crystal structure of the complex formed with the mutant with both Cys residues replaced was nearly identical with that of the complex containing the wild-type protein, with the Ser oxygen atoms positioned to replace the disulfide bond with a hydrogen bond. The two structures of the complexes with single replacements displayed small local perturbations with alternate conformations of the Ser side chains. Despite the absence of the disulfide bond, the crystallographic temperature factors show no evidence of increased flexibility in the complexes with the mutant inhibitors. All three of the variants were cleaved by trypsin more rapidly than the wild-type inhibitor, by as much as 10,000-fold, indicating that the covalent constraint normally imposed by the disulfide contributes to the remarkable resistance to hydrolysis displayed by the wild-type protein. The rates of hydrolysis display an unusual dependence on pH over the range of 3.5-8.0, decreasing at the more alkaline values, as compared with the increased hydrolysis rates for normal substrates under these conditions. These observations can be accounted for by a model for inhibition in which an acyl-enzyme intermediate forms at a significant rate but is rapidly converted back to the enzyme-inhibitor complex by nucleophilic attack by the newly created amino group. The model suggests that a lack of flexibility in the acyl-enzyme intermediate, rather than the enzyme-inhibitor complex, may be a key factor in the ability of bovine pancreatic trypsin inhibitor and similar inhibitors to resist hydrolysis.

Factor Xa active site substrate specificity with substrate phage display and computational molecular modeling

Structural origin of substrate-enzyme recognition remains incompletely understood. In the model enzyme system of serine protease, canonical anti-parallel beta-structure substrate-enzyme complex is the predominant hypothesis for the substrate-enzyme interaction at the atomic level. We used factor Xa (fXa), a key serine protease of the coagulation system, as a model enzyme to test the canonical conformation hypothesis. More than 160 fXa-cleavable substrate phage variants were experimentally selected from three designed substrate phage display libraries. These substrate phage variants were sequenced and their specificities to the model enzyme were quantified with quantitative enzyme-linked immunosorbent assay for substrate phage-enzyme reaction kinetics. At least three substrate-enzyme recognition modes emerged from the experimental data as necessary to account for the sequence-dependent specificity of the model enzyme. Computational molecular models were constructed, with both energetics and pharmacophore criteria, for the substrate-enzyme complexes of several of the representative substrate peptide sequences. In contrast to the canonical conformation hypothesis, the binding modes of the substrates to the model enzyme varied according to the substrate peptide sequence, indicating that an ensemble of binding modes underlay the observed specificity of the model serine protease.

Peptidomimetics: Fmoc solid-phase pseudopeptide synthesis

Peptidomimetic modifications or cyclization of linear peptides are frequently used as attractive methods to provide more conformationally constrained and thus more stable and bioactive peptides. Among numerous peptidomimetic approaches described recently in the literature, particularly attractive are pseudopeptides or peptide bond surrogates in which peptide bonds have been replaced with other chemical groups. In these peptidomimetics the amide bond surrogates possess three-dimensional structures similar to those of natural peptides, yet with significant differences in polarity, hydrogen bonding capability, and acid-base character. The introduction of such modifications to the peptide sequence is expected to completely prevent protease cleavage of amide bond and significantly improve peptides' metabolic stability. In this chapter we consider Fmoc solid-phase synthesis of peptide analogs containing the amide surrogate that tend to be isosteric with the natural amide. This includes synthesis of peptidosulfonamides, phosphonopeptides, oligoureas, depsides, depsipeptides, and peptoids.

Probing the substrate specificity of the dengue virus type 2 NS3 serine protease by using internally quenched fluorescent peptides

The NS3 (dengue virus non-structural protein 3) serine protease of dengue virus is an essential component for virus maturation, thus representing an attractive target for the development of antiviral drugs directed at the inhibition of polyprotein processing. In the present study, we have investigated determinants of substrate specificity of the dengue virus NS3 protease by using internally quenched fluorogenic peptides containing Abz (o-aminobenzoic acid; synonymous to anthranilic acid) and 3-nitrotyrosine (nY) representing both native and chimaeric polyprotein cleavage site sequences. By using this combinatorial approach, we were able to describe the substrate preferences and determinants of specificity for the dengue virus NS2B(H)-NS3pro protease. Kinetic parameters (kcat/K(m)) for the hydrolysis of peptide substrates with systematic truncations at the prime and non-prime side revealed a length preference for peptides spanning the P4-P3' residues, and the peptide Abz-RRRRSAGnY-amide based on the dengue virus capsid protein processing site was discovered as a novel and efficient substrate of the NS3 protease (kcat/K(m)=11087 M(-1) x s(-1)). Thus, while having confirmed the exclusive preference of the NS3 protease for basic residues at the P1 and P2 positions, we have also shown that the presence of basic amino acids at the P3 and P4 positions is a major specificity-determining feature of the dengue virus NS3 protease. Investigation of the substrate peptide Abz-KKQRAGVLnY-amide based on the NS2B/NS3 polyprotein cleavage site demonstrated an unexpected high degree of cleavage efficiency. Chimaeric peptides with combinations of prime and non-prime sequences spanning the P4-P4' positions of all five native polyprotein cleavage sites revealed a preponderant effect of non-prime side residues on the K(m) values, whereas variations at the prime side sequences had higher impact on kcat.

Enzyme:substrate hydrogen bond shortening during the acylation phase of serine protease catalysis

Atomic resolution (<or=1.2 A) serine protease intermediate structures revealed that the strength of the hydrogen bonds between the enzyme and the substrate changed during catalysis. The well-conserved hydrogen bonds of antiparallel beta-sheet between the enzyme and the substrate become significantly shorter in the transition from a Michaelis complex analogue (Pontastacus leptodactylus (narrow-fingered crayfish) trypsin (CFT) in complex with Schistocerca gregaria (desert locust) trypsin inhibitor (SGTI) at 1.2 A resolution) to an acyl-enzyme intermediate (N-acetyl-Asn-Pro-Ile acyl-enzyme intermediate of porcine pancreatic elastase at 0.95 A resolution) presumably synchronously with the nucleophilic attack on the carbonyl carbon atom of the scissile peptide bond. This is interpreted as an active mechanism that utilizes the energy released from the stronger hydrogen bonds to overcome the energetic barrier of the nucleophilic attack by the hydroxyl group of the catalytic serine. In the CFT:SGTI complex this hydrogen bond shortening may be hindered by the 27I-32I disulfide bridge and Asn-15I of SGTI. The position of the catalytic histidine changes slightly as it adapts to the different nucleophilic attacker during the transition from the Michaelis complex to the acyl-enzyme state, and simultaneously its interaction with Asp-102 and Ser-214 becomes stronger. The oxyanion hole hydrogen bonds provide additional stabilization for acyl-ester bond in the acyl-enzyme than for scissile peptide bond of the Michaelis complex. Significant deviation from planarity is not observed in the reactive bonds of either the Michaelis complex or the acyl-enzyme. In the Michaelis complex the electron distribution of the carbonyl bond is distorted toward the oxygen atom compared to other peptide bonds in the structure, which indicates the polarization effect of the oxyanion hole.

Backbone-Backbone H-Bonds Make Context-Dependent Contributions to Protein Folding Kinetics and Thermodynamics: Lessons from Amide-to-Ester Mutations

The contribution of backbone-backbone hydrogen bonds (H-bonds) to protein folding energetics has been controversial. This is due, at least in part, to the inability to perturb backbone-backbone H-bonds by traditional methods of protein mutagenesis. Recently, however, protein backbone mutagenesis has become possible with the development of chemical and biological methods to replace individual amides in the protein backbone with esters. Here, we review the use of amide-to-ester mutation as a tool to evaluate the contribution of backbone-backbone H-bonds to protein folding kinetics and thermodynamics.

Toward Assessing the Position-Dependent Contributions of Backbone Hydrogen Bonding to β-Sheet Folding Thermodynamics Employing Amide-to-Ester Perturbations

An amide-to-ester backbone substitution in a protein is accomplished by replacing an alpha-amino acid residue with the corresponding alpha-hydroxy acid, preserving stereochemistry, and conformation of the backbone and the structure of the side chain. This substitution replaces the amide NH (a hydrogen bond donor) with an ester O (which is not a hydrogen bond donor) and the amide carbonyl (a strong hydrogen bond acceptor) with an ester carbonyl (a weaker hydrogen bond acceptor), thus perturbing folding energetics. Amide-to-ester perturbations were used to evaluate the thermodynamic contribution of each hydrogen bond in the PIN WW domain, a three-stranded beta-sheet protein. Our results reveal that removing a hydrogen bond donor destabilizes the native state more than weakening a hydrogen bond acceptor and that the degree of destabilization is strongly dependent on the location of the amide bond replaced. Hydrogen bonds near turns or at the ends of beta-strands are less influential than hydrogen bonds that are protected within a hydrophobic core. Beta-sheet destabilization caused by an amide-to-ester substitution cannot be directly related to hydrogen bond strength because of differences in the solvation and electrostatic interactions of amides and esters. We propose corrections for these differences to obtain approximate hydrogen bond strengths from destabilization energies. These corrections, however, do not alter the trends noted above, indicating that the destabilization energy of an amide-to-ester mutation is a good first-order approximation of the free energy of formation of a backbone amide hydrogen bond.

Binding of nonphysiological protein and peptide substrates to proteases: differences between urokinase-type plasminogen activator and trypsin and contributions to the evolution of regulated proteolysis

Understanding the regulation of physiological processes requires detailed knowledge of the recognition of substrates by enzymes. One of the most productive model systems for the study of enzyme-substrate interactions is the serine protease family; however, most studies of protease action have used small substrates that contain an activated, non-natural scissile bond. Because few kinetic or structural studies have used protein substrates, the physiologically relevant target of most proteases, it seems likely that important mechanisms of substrate recognition and processing by proteases have not yet been fully elucidated. Consistent with this hypothesis, we have observed that K(m) values for protein substrates are reduced as much as 200-15000-fold relative to those of analogous peptide substrates. Here we examine the thermodynamic consequences of interactions between proteases and their substrates using staphylococcal nuclease (SNase) and SNase variants as model protein substrates. We have obtained values for enthalpy, entropy, and K(d) for binding of proteins and peptides by the nonspecific protease trypsin and the highly specific protease urokinase-type plasminogen activator (u-PA). To avoid cleavage of substrates during these measurements, we used inactive variants of trypsin and u-PA whose catalytic serine S195 had been replaced by alanine. Differences in the K(d) values for binding of protein and peptide substrates closely approximate the large differences observed in the corresponding K(m) values. Improved binding of protein substrates is due to decreased enthalpy, and this effect is pronounced for the selective protease u-PA. Fundamental differences in recognition of analogous protein and peptide substrates may have influenced the evolution of protease specificity.

Generation of ligand conformations in continuum solvent consistent with protein active site topology: application to thrombin

Using the crystal structure of an inhibitor complexed with the serine protease thrombin (PDB code ) and the functional group definitions contained within the Catalyst software, a representation of the enzyme's active site was produced (structure-based pharmacophore model). A training set of 16 homologous non-peptide inhibitors whose conformations had been generated in continuum solvent (MacroModel) and clustered into conformational families (XCluster) was regressed against this pharmacophore so as to obtain a 3D-QSAR model. To test the robustness of the resulting QSAR model, the synthesis of a series of non-peptide thrombin inhibitors based on arylsuphonyl derivatives of an aminophenol ring linked to a pyridyl-based S1 binding group was undertaken. These compounds served as a test set (20-24). The crystal structure for the novel symmetrical disulfonyl compound 24, in complex with thrombin, has been solved. Its calculated binding mode is in general agreement with the crystallographically observed one, and the predicted K(i) value is in close accord with the experimental value.

Probing the role of backbone hydrogen bonding in beta-amyloid fibrils with inhibitor peptides containing ester bonds at alternate positions

Protein-protein interactions are frequently mediated by stable, intermolecular beta-sheets. A number of cytokines and the HIV Protease, for example, dimerize through beta-sheet motifs. Evidence also suggests that the macromolecular assemblies of peptides and proteins in amyloid fibrils are stabilized by intermolecular beta-sheets. In this paper, we report that interfering with the backbone hydrogen bonding of an amyloidgenic peptide (Abeta16-20) by replacing amide bonds with ester bonds prevents the aggregation of the peptide. The ester bonds were incorporated in an alternating fashion so that the peptide presents two unique hydrogen bonding faces when arrayed in an extended, beta-strand conformation; one face of the peptide has normal hydrogen bonding capabilities, but the other face is missing amide protons and its ability to hydrogen bond is severely limited. Analytical ultracentrifugation experiments demonstrate that this ester peptide, Abeta16-20e, is predominantly monomeric under solution conditions, unlike the fibril-forming Abeta16-20 peptide. Abeta16-20e also inhibits the aggregation of the Abeta1-40 peptide and disassembles preformed Abeta1-40 fibrils. These results suggest that backbone hydrogen bonding is critical for the assembly of amyloid fibrils.

A clogged gutter mechanism for protease inhibitors

A classical peptide inhibitor of serine proteases that is hydrolyzed approximately 10(7) times more slowly than a good substrate is shown to form an acyl-enzyme intermediate rapidly. Despite this quick first step, further reaction is slowed dramatically because of tight and oriented binding of the cleaved peptide, preventing acyl-enzyme hydrolysis and favoring the reverse reaction. Moreover, this mechanism appears to be common to a large class of tight-binding serine protease inhibitors that mimic good substrates. The arrest of enzymatic reaction at the intermediate stage allowed us to determine that the consensus nucleophilic attack angle is close to 90 degrees in the reactive Michaelis complexes.

The methyl group of N(alpha)(Me)Arg-containing peptides disturbs the active-site geometry of thrombin, impairing efficient cleavage

Bivalent peptidic thrombin inhibitors consisting of an N-terminal d-cyclohexylalanine-Pro-N(alpha)(Me)Arg active-site fragment, a flexible polyglycine linker, and a C-terminal hirugen-like segment directed towards the fibrinogen recognition exosite inhibit thrombin with K(i) values in the picomolar range, remaining stable in buffered solution at pH 7.8 for at least 15 hours. In order to investigate the structural basis of this increased stability, the most potent of these inhibitors, I-11 (K(i)=37pM), containing an N(alpha)(Me)Arg-Thr bond, was crystallized in complex with human alpha-thrombin. X-ray data were collected to 1.8A resolution and the crystal structure of this complex was determined. The Fourier map displays clear electron density for the N-terminal fragment and for the exosite binding segment. It indicates, however, that in agreement with Edman sequencing, the peptide had been cleaved in the crystal, presumably due to the long incubation time of 14 days needed for crystallization and data collection. The N(alpha)(Me) group is directed toward the carbonyl oxygen atom of Ser214, pushing the Ser195 O(gamma) atom out of its normal site. This structure suggests that upon thrombin binding, the scissile peptide bond of the intact peptide and the Ser195 O(gamma) are separated from each other, impairing the nucleophilic attack of the Ser195 O(gamma) toward the N(alpha)(Me)Arg carbonyl group. In the time-scale of two weeks, however, cleavage geometries favoured by the crystal allow catalysis at a slow rate.

In vitro selection of RNA molecules that inhibit the activity of ricin A-chain

The cytotoxin ricin disables translation by depurinating a conserved site in eukaryotic rRNA. In vitro selection has been used to generate RNA ligands (aptamers) specific for the catalytic ricin A-chain (RTA). The anti-RTA aptamers bear no resemblance to the normal RTA substrate, the sarcin-ricin loop (SRL), and were not depurinated by RTA. An initial 80-nucleotide RNA ligand was minimized to a 31-nucleotide aptamer that contained all sequences and structures necessary for interacting with RTA. This minimal RNA formed high affinity complexes with RTA (K(d) = 7.3 nM) which could compete directly with the SRL for binding to RTA. The aptamer inhibited RTA depurination of the SRL and could partially protect translation from RTA inhibition. The IC(50) of the aptamer for RTA in an in vitro translation assay is 100 nM, roughly 3 orders of magnitude lower than a small molecule inhibitor of ricin, pteroic acid, and 2 orders of magnitude lower than the best known RNA inhibitor. The novel anti-RTA aptamers may find application as diagnostic reagents for a potential biological warfare agent and hold promise as scaffolds for the development of strong ricin inhibitors.