Ionizable P1 residues in serine proteinase inhibitors undergo large pK shifts on complex formation

Dissecting the Energetics of Intrinsically Disordered Proteins via a Hybrid Experimental and Computational Approach

Intrinsically disordered proteins (IDPs) play important roles in biology, but little is known about the energetics of their inter-residue interactions. Methods that have been successfully applied to analyze the energetics of globular proteins are not applicable to the fluctuating partially ordered ensembles populated by IDPs. A combined computational experimental strategy is introduced for analyzing the energetic role of individual residues in the free state of IDPs. The approach combines experimental measurements of the binding of wild-type and mutant IDPs to their partners with alchemical free energy calculations of the structured complexes. These data allow quantitative information to be deduced about the free state via a thermodynamic cycle. The approach is validated by the analysis of the effects of mutations upon the binding free energy of the ovomucoid inhibitor third binding domain to its partners and is applied to the C-terminal domain of the measles virus nucleoprotein, a 125-residue IDP involved in the RNA transcription and replication of measles virus. The analysis reveals significant inter-residue interactions in the unbound IDP and suggests a biological role for them. The work demonstrates that advances in force fields and computational hardware have now led to the point where it is possible to develop methods, which integrate experimental and computational techniques to reveal insights that cannot be studied using either technique alone.

Biomolecular NMR: Past and future

The editors of this special volume suggested this topic, presumably because of the perspective lent by our combined >90-year association with biomolecular NMR. What follows is our personal experience with the evolution of the field, which we hope will illustrate the trajectory of change over the years. As for the future, one can confidently predict that it will involve unexpected advances. Our narrative is colored by our experience in using the NMR Facility for Biomedical Studies at Carnegie-Mellon University (Pittsburgh) and in developing similar facilities at Purdue (1977-1984) and the University of Wisconsin-Madison (1984-). We have enjoyed developing NMR technology and making it available to collaborators and users of these facilities. Our group's association with the Biological Magnetic Resonance data Bank (BMRB) and with the Worldwide Protein Data Bank (wwPDB) has also been rewarding. Of course, many groups contributed to the early growth and development of biomolecular NMR, and our brief personal account certainly omits many important milestones.

Additivity-based design of the strongest possible turkey ovomucoid third domain inhibitors for porcine pancreatic elastase (PPE) and Streptomyces griseus protease B (SGPB)

We describe here successful designs of strong inhibitors for porcine pancreatic elastase (PPE) and Streptomyces griseus protease B (SGPB). For each enzyme two inhibitor variants were designed. In one, the reactive site residue (position 18) was retained and the best residues were substituted at contact positions 13, 14, and 15. In the other variant the best residues were substituted at all contact positions except the reactive site where a Gly was substituted. The four designed variants were: for PPE, T(13)E(14)Y(15)-OMTKY3 and T(13)E(14)Y(15)G(18)M(21)P(32)V(36)-OMTKY3, and for SGPB, S(13)D(14)Y(15)-OMTKY3 and S(13)D(14)Y(15)G(18)I(19)K(21)-OMTKY3. The free energies of association (ΔG(0)) of expressed variants have been measured with the proteases for which they were designed as well as with five other serine proteases and the results are discussed.

Protonation and pK changes in protein-ligand binding

Formation of protein-ligand complexes causes various changes in both the receptor and the ligand. This review focuses on changes in pK and protonation states of ionizable groups that accompany protein-ligand binding. Physical origins of these effects are outlined, followed by a brief overview of the computational methods to predict them and the associated corrections to receptor-ligand binding affinities. Statistical prevalence, magnitude and spatial distribution of the pK and protonation state changes in protein-ligand binding are discussed in detail, based on both experimental and theoretical studies. While there is no doubt that these changes occur, they do not occur all the time; the estimated prevalence varies, both between individual complexes and by method. The changes occur not only in the immediate vicinity of the interface but also sometimes far away. When receptor-ligand binding is associated with protonation state change at particular pH, the binding becomes pH dependent: we review the interplay between sub-cellular characteristic pH and optimum pH of receptor-ligand binding. It is pointed out that there is a tendency for protonation state changes upon binding to be minimal at physiologically relevant pH for each complex (no net proton uptake/release), suggesting that native receptor-ligand interactions have evolved to reduce the energy cost associated with ionization changes. As a result, previously reported statistical prevalence of these changes - typically computed at the same pH for all complexes - may be higher than what may be expected at optimum pH specific to each complex. We also discuss whether proper account of protonation state changes appears to improve practical docking and scoring outcomes relevant to structure-based drug design. An overview of some of the existing challenges in the field is provided in conclusion.

pH-selective mutagenesis of protein-protein interfaces: in silico design of therapeutic antibodies with prolonged half-life

Understanding the effects of mutation on pH-dependent protein binding affinity is important in protein design, especially in the area of protein therapeutics. We propose a novel method for fast in silico mutagenesis of protein–protein complexes to calculate the effect of mutation as a function of pH. The free energy differences between the wild type and mutants are evaluated from a molecular mechanics model, combined with calculations of the equilibria of proton binding. The predicted pH-dependent energy profiles demonstrate excellent agreement with experimentally measured pH-dependency of the effect of mutations on the dissociation constants for the complex of turkey ovomucoid third domain (OMTKY3) and proteinase B. The virtual scanning mutagenesis identifies all hotspots responsible for pH-dependent binding of immunoglobulin G (IgG) to neonatal Fc receptor (FcRn) and the results support the current understanding of the salvage mechanism of the antibody by FcRn based on pH-selective binding. The method can be used to select mutations that change the pH-dependent binding profiles of proteins and guide the time consuming and expensive protein engineering experiments. As an application of this method, we propose a computational strategy to search for mutations that can alter the pH-dependent binding behavior of IgG to FcRn with the aim of improving the half-life of therapeutic antibodies in the target organism.

Cleavage of peptide bonds bearing ionizable amino acids at P(1) by serine proteases with hydrophobic S(1) pocket

Enzymatic hydrolysis of the synthetic substrate succinyl-Ala-Ala-Pro-Xxx-pNA (where Xxx=Leu, Asp or Lys) catalyzed by bovine chymotrypsin (CHYM) or Streptomyces griseus protease B (SGPB) has been studied at different pH values in the pH range 3-11. The pH optima for substrates having Leu, Asp, and Lys have been found to be 7.5-8.0, 5.5-6.0, and ∼10, respectively. At the normally reported pH optimum (pH 7-8) of CHYM and SGPB, the substrate with Leu at the reactive site is more than 25,000-fold more reactive than that with Asp. However, when fully protonated, Asp is nearly as good a substrate as Leu. The pK values of the side chains of Asp and Lys in the hydrophobic S(1) pocket of CHYM and SGPB have been calculated from pH-dependent hydrolysis data and have been found to be about 9 for Asp and 7.4 and 9.7 for Lys for CHYM and SGPB, respectively. The results presented in this communication suggest a possible application of CHYM like enzymes in cleaving peptide bonds contributed by acidic amino acids between pH 5 and 6.

Protein-protein binding is often associated with changes in protonation state

pK(a) values of ionizable residues have been calculated using the PROPKA method and structures of 75 protein-protein complexes and their corresponding free forms. These pK(a) values were used to compute changes in protonation state of individual residues, net changes in protonation state of the complex relative to the uncomplexed proteins, and the correction to a binding energy calculated assuming standard protonation states at pH 7. For each complex, two different structures for the uncomplexed form of the proteins were used: the X-ray structures determined for the proteins in the absence of the other protein and the individual protein structures taken from the structure of the complex (referred to as unbound and bound structures, respectively). In 28 and 77% of the cases considered here, protein-protein binding is accompanied by a complete (>95%) or significant (>50%) change in protonation state of at least one residue using unbound structures. Furthermore, in 36 and 61% of the cases, protein-protein binding is accompanied by a complete or significant net change in protonation state of the complex relative to the separated monomers. Using bound structures, the corresponding values are 12, 51, 20, and 48%. Comparison to experimental data suggest that using unbound and bound structures lead to over- and underestimation of binding-induced protonation state changes, respectively. Thus, we conclude that protein-protein binding is often associated with changes in protonation state of amino acid residues and with changes in the net protonation state of the proteins. The pH-dependent correction to the binding energy contributes at least one order of magnitude to the binding constant in 45 and 23%, using unbound and bound structures, respectively.

1.2A-resolution crystal structures reveal the second tetrahedral intermediates of streptogrisin B (SGPB)

Streptogrisin B (SGPB) has served as one of the models for studying the catalytic activities of serine peptidases. Here we report its native crystal structures at pH 4.2 at a resolution of 1.18A, and at pH 7.3 at a resolution of 1.23A. Unexpectedly, outstanding electron density peaks occurred in the active site and the substrate-binding region of SGPB in the computed maps at both pHs. The densities at pH 4.2 were assigned as a tetrapeptide, Asp-Ala-Ile-Tyr, whereas those at pH 7.3 were assigned as a tyrosine molecule and a leucine molecule existing at equal occupancies in both of the SGPB molecules in the asymmetric unit. Refinement with relaxed geometric restraints resulted in molecular structures representing mixtures of the second tetrahedral intermediates and the enzyme-product complexes of SGPB existing in a pH-dependent equilibrium. Structural comparisons with the complexes of SGPB with turkey ovomucoid third domain (OMTKY3) and its variants have shown that, upon the formation of the tetrahedral intermediate, residues Glu192A to Gly193 of SGPB move towards the alpha-carboxylate O of residue P1 of the bound species, and adjustments in the side-chain conformational angles of His57 and Ser195 of SGPB favor the progression of the catalytic mechanism of SGPB.

Free energy calculations show that acidic P1 variants undergo large pKa shifts upon binding to trypsin

Serine proteinases and their protein inhibitors belong to one of the most comprehensively studied models of protein-protein interactions. It is well established that the narrow trypsin specificity is caused by the presence of a negatively charged aspartate at the specificity pocket. X-ray crystallography as well as association measurements revealed, surprisingly, that BPTI with glutamatic acid as the primary binding (P1) residue was able to bind to trypsin. Previous free energy calculations showed that there was a substantially unfavorable binding free energy associated with accommodation of ionized P1 Glu at the S1-site of trypsin. In this study, the binding of P1 Glu to trypsin has been systematically investigated in terms of the protonation states of P1 Glu and Asp189, the orientation of Gln192, as well as the possible presence of counterions using the linear interaction energy (LIE) approach and the free energy perturbation (FEP) method. Twenty-four conceivable binding arrangements were evaluated and quantitative agreement with experiments is obtained when the P1 Glu binds in its protonated from. The results suggest that P1 Glu is one of the variants of BPTI that inhibit trypsin strongest at low pH, contrary to the specificity profile of trypsin, suggesting a new regulation mechanism of trypsin-like enzymes.

Predicting proteinase specificities from free energy calculations

The role of the primary binding residue (P1) in complexes between three different subtilases (subtilisin Carlsberg, thermitase and proteinase K) and their canonical protein inhibitor eglin c have been studied by free energy calculations. Based on the crystal structures of eglin c in complex with subtilisin Carlsberg and thermitase, and a homology model of the eglin c-proteinase K complex, a total of 57 mutants have been constructed and docked into their host proteins. The binding free energy was then calculated using molecular dynamics (MD) simulations combined with the linear interaction energy (LIE) method for all complexes differing only in the nature of the amino acid at the P1 position. LIE calculations for 19 different complexes for each subtilase were thus carried out excluding proline. The effects of substitutions at the P1 position on the binding free energies are found to be very large, and positively charged residues (Arg, Lys and His) are particularly deleterious for all three enzymes. The charged variants of the acidic side chains are found to bind more favorably as compared to their protonated states in all three subtilases. Furthermore, hydrophobic amino acids are accommodated most favorably at the S1-site in all three enzymes. Comparison of the three series of binding free energies shows only minor differences in the 19 computed relative binding free energies among these subtilases. This is further reflected in the correlation coefficient between the 23 relative binding free energies obtained, including the possible protonation states of ionizable side chains, but excluding the P1 Pro, for subtilisin Carlsberg versus thermitase (0.95), subtilisin versus proteinase K (0.94) and thermitase versus proteinase K (0.96).

Probing the effect of point mutations at protein-protein interfaces with free energy calculations

We have studied the effect of point mutations of the primary binding residue (P1) at the protein-protein interface in complexes of chymotrypsin and elastase with the third domain of the turkey ovomucoid inhibitor and in trypsin with the bovine pancreatic trypsin inhibitor, using molecular dynamics simulations combined with the linear interaction energy (LIE) approach. A total of 56 mutants have been constructed and docked into their host proteins. The free energy of binding could be reliably calculated for 52 of these mutants that could unambiguously be fitted into the binding sites. We find that the predicted binding free energies are in very good agreement with experimental data with mean unsigned errors between 0.50 and 1.03 kcal/mol. It is also evident that the standard LIE model used to study small drug-like ligand binding to proteins is not suitable for protein-protein interactions. Three different LIE models were therefore tested for each of the series of protein-protein complexes included, and the best models for each system turn out to be very similar. The difference in parameterization between small drug-like compounds and protein point mutations is attributed to the preorganization of the binding surface. Our results clearly demonstrate the potential of free energy calculations for probing the effect of point mutations at protein-protein interfaces and for exploring the principles of specificity of hot spots at the interface.

Crystal structures of five bovine chymotrypsin complexes with P1 BPTI variants

The bovine chymotrypsin-bovine pancreatic trypsin inhibitor (BPTI) interaction belongs to extensively studied models of protein-protein recognition. The accommodation of the inhibitor P1 residue in the S1 binding site of the enzyme forms the hot spot of this interaction. Mutations introduced at the P1 position of BPTI result in a more than five orders of magnitude difference of the association constant values with the protease. To elucidate the structural aspects of the discrimination between different P1 residues, crystal structures of five bovine chymotrypsin-P1 BPTI variant complexes have been determined at pH 7.8 to a resolution below 2 A. The set includes polar (Thr), ionizable (Glu, His), medium-sized aliphatic (Met) and large aromatic (Trp) P1 residues and complements our earlier studies of the interaction of different P1 side-chains with the S1 pocket of chymotrypsin. The structures have been compared to the complexes of proteases with similar and dissimilar P1 preferences, including Streptomyces griseus proteases B and E, human neutrophil elastase, crab collagenase, bovine trypsin and human thrombin. The S1 sites of these enzymes share a common general shape of significant rigidity. Large and branched P1 residues adapt in their complexes similar conformations regardless of the polarity and size differences between their S1 pockets. Conversely, long and flexible residues such as P1 Met are present in the disordered form and display a conformational diversity despite similar inhibitory properties with respect to most enzymes studied. Thus, the S1 specificity profiles of the serine proteases appear to result from the precise complementarity of the P1-S1 interface and minor conformational adjustments occurring upon the inhibitor binding.

NMR determination of pKa values for Asp, Glu, His, and Lys mutants at each variable contiguous enzyme-inhibitor contact position of the turkey ovomucoid third domain

From the larger set of 191 variants at all the variable contact positions in the turkey ovomucoid third domain, we selected a subset that consists of Asp, Glu, His, and Lys residues at eight of the nine contiguous P6-P3' positions (residues 13-21), the exception being P3-Cys16 which is involved in a conserved disulfide bridge. Two-dimensional [1H,1H]-TOCSY data were collected for each variant as a function of sample pH. This allowed for the evaluation of 31 of the 32 pK(a) values for these residues, the exception being that of P5-Lys14, whose signals at high pH could not be resolved from those of other Lys residues in the molecule. Only two of the titrating residues are present in the wild-type protein (P6-Lys13 and P1'-Glu19); hence, these measurements complement earlier measurements by A. D. Robertson and co-workers. This data set was supplemented with results from the pH dependence of NMR spectra of four additional single mutants, P1-Leu18Gly, P1-Leu18Ala, P2-Thr17Val, and P3'-Arg21Ala, and two double mutants, P2-Thr17Val/P3'-Arg21Ala and P8-Tyr11Phe/P6-Lys13Asp. Probably the most striking result was observation of a P2-Thr17...P1'-Glu19 hydrogen bond and a P1'-Glu19-P3'-Arg21 electrostatic interaction within the triad of P2, P1', and P3' (residues 17, 19, and 21, respectively). In several cases, the pK(a) of a particular residue was sensed by resonances not only in that residue but also in residue(s) with which it interacts. Remarkably, in several interacting systems, resonances from different protons within the same residue yielded different pHmid values.

Additivity-based prediction of equilibrium constants for some protein-protein associations

For many protein families, such as serine proteinases or serine proteinase inhibitors, the family assignment predicts reactivity only in general terms. Both detailed specificity and quantitative reactivity are lacking. We believe that, for many such protein families, algorithms can be devised by defining the subset of n functionally important sequence positions, making the 19n possible single mutants and measuring their reactivity. Given the assumption that the contributions of the n positions are additive, the reactivities of the 20(n) variants can be predicted. This is illustrated by an almost complete algorithm for the Kazal family of protein inhibitors of serine proteinases.

Computational analysis of binding of P1 variants to trypsin

The binding of P1 variants of bovine pancreatic trypsin inhibitor (BPTI) to trypsin has been investigated by means of molecular dynamics simulations. The specific interaction formed between the amino acid at the primary binding (P1) position of the binding loop of BPTI and the specificity pocket of trypsin was estimated by use of the linear interaction energy (LIE) method. Calculations for 13 of the naturally occurring amino acids at the P1 position were carried out, and the results obtained were found to correlate well with the experimental binding free energies. The LIE calculations rank the majority of the 13 variants correctly according to the experimental association energies and the mean error between calculated and experimental binding free energies is only 0.38 kcal/mole, excluding the Glu and Asp variants, which are associated with some uncertainties regarding protonation and the possible presence of counter-ions. The three-dimensional structures of the complex with three of the P1 variants (Asn, Tyr, and Ser) included in this study have not at present been solved by any experimental techniques and, therefore, were modeled on the basis of experimental data from P1 variants of similar size. Average structures were calculated from the MD simulations, from which specific interactions explaining the broad variation in association energies were identified. The present study also shows that explicit treatment of the complex water-mediated hydrogen bonding network at the protein-protein interface is of crucial importance for obtaining reliable binding free energies. The successful reproduction of relative binding energies shows that this type of methodology can be very useful as an aid in rational design and redesign of biologically active macromolecules.

NMR chemical shift mapping of the binding site of a protein proteinase inhibitor: changes in the (1)H, (13)C and (15)N NMR chemical shifts of turkey ovomucoid third domain upon binding to bovine chymotrypsin A(alpha)

The substrate-like inhibition of serine proteinases by avian ovomucoid domains has provided an excellent model for protein inhibitor-proteinase interactions of the standard type. 1H,15N and 13C NMR studies have been undertaken on complexes formed between turkey ovomucoid third domain (OMTKY3)2 and chymotrypsin A(alpha) (Ctr) in order to characterize structural changes occurring in the Ctr binding site of OMTKY3. 15N and 13C were incorporated uniformly into OMTKY3, allowing backbone resonances to be assigned for OMTKY3 in both its free and complex states. Chemical shift perturbation mapping indicates that the two regions, K13-P22 and N33-A40, are the primary sites in OMTKY3 involved in Ctr binding, in full agreement with the 12 consensus proteinase-contact residues of OMTKY3 defined previously on the basis of X-ray crystallographic and mutational analysis. Smaller chemical shift perturbations in selected other regions may result from minor structural changes on binding. Through-bond 15N-13C correlations between P1-13C' and P1'-15N in two-dimensional H(N)CO and HN(CO) NMR spectra of selectively labeled OMTKY3 complexed with Ctr indicate that the scissile peptide bond between L18 and E19 of the inhibitor is intact in the complex. The chemical shifts of the reactive site peptide bond indicate that it is predominantly trigonal, although the data are not inconsistent with a slight perturbation of the hybridization of the peptide bond toward the first tetrahedral state along the reaction coordinate.

The ionization of a buried glutamic acid is thermodynamically linked to the stability of Leishmania mexicana triose phosphate isomerase

The amino acid sequence of Leishmania mexicana triose phosphate isomerase is unique in having at position 65 a glutamic acid instead of a glutamine. The stability properties of LmTIM and the E65Q mutant were investigated by pH and guanidinium chloride-induced unfolding. The crystal structure of E65Q was determined. Three important observations were made: (a) there are no structural rearrangements as the result of the substitution; (b) the mutant is more stable than the wild-type; and (c) the stability of the wild-type enzyme shows strong pH dependence, which can be attributed to the ionization of Glu65. Burying of the Glu65 side chain in the uncharged environment of the dimer interface results in a shift in pKa of more than 3 units. The pH-dependent decrease in overall stability is due to weakening of the monomer-monomer interactions (in the dimer). The E65Q substitution causes an increase in stability as the result of the formation of an additional hydrogen bond in each subunit (DeltaDeltaG degrees of 2 kcal.mol-1 per monomer) and the elimination of a charged group in the dimer interface (DeltaDeltaG degrees of at least 9 kcal.mol-1 per dimer). The computated shift in pKa and the stability of the dimer calculated from the charge distribution in the protein structure agree closely with the experimental results. The guanidinium chloride dependence of the unfolding constant was smaller than expected from studies involving monomeric model proteins. No intermediates could be identified in the unfolding equilibrium by combining fluorescence and CD measurements. Study of a stable monomeric triose phosphate isomerase variant confirmed that the phenomenon persists in the monomer.

Deleterious effects of beta-branched residues in the S1 specificity pocket of Streptomyces griseus proteinase B (SGPB): crystal structures of the turkey ovomucoid third domain variants Ile18I, Val18I, Thr18I, and Ser18I in complex with SGPB

Turkey ovomucoid third domain (OMTKY3) is a canonical inhibitor of serine proteinases. Upon complex formation, the inhibitors fully exposed P1 residue becomes fully buried in the preformed cavity of the enzyme. All 20 P1 variants of OMTKY3 have been obtained by recombinant DNA technology and their equilibrium association constants have been measured with six serine proteinases. To rationalize the trends observed in this data set, high resolution crystal structures have been determined for OMTKY3 P1 variants in complex with the bacterial serine proteinase, Streptomyces griseus proteinase B (SGPB). Four high resolution complex structures are being reported in this paper; the three beta-branched variants, Ile18I, Val18I, and Thr18I, determined to 2.1, 1.6, and 1.7 A resolution, respectively, and the structure of the Ser18I variant complex, determined to 1.9 A resolution. Models of the Cys18I, Hse18I, and Ape18I variant complexes are also discussed. The beta-branched side chains are not complementary to the shape of the S1 binding pocket in SGPB, in contrast to that of the wild-type gamma-branched P1 residue for OMTKY3, Leu18I. Chi1 angles of approximately 40 degrees are imposed on the side chains of Ile18I, Val18I, and Thr18I within the S1 pocket. Dihedral angles of +60 degrees, -60 degrees, or 180 degrees are more commonly observed but 40 degrees is not unfavorable for the beta-branched side chains. Thr18I Ogamma1 also forms a hydrogen bond with Ser195 Ogamma in this orientation. The Ser18I side chain adopts two alternate conformations within the S1 pocket of SGPB, suggesting that the side chain is not stable in either conformation.

Computational analysis of the binding of P1 variants of domain 3 of turkey ovomucoid inhibitor to Streptomyces griseus protease B

Binding constants for complexes of variants of the ovomucoid inhibitor domain 3 from turkey (OMTKY3) and Streptomyces griseus protease B (SGPB) have been computed. On the basis of the crystallographically determined structures of the complexes, continuum electrostatic calculations have been carried out to evaluate the electrostatic contribution to the binding energy. The hydrophobic component was computed based on the change in the solvent accessible surface area on complex formation. These two terms were combined linearly and the parameters for the protein dielectric, atomic solvation parameter and a constant term were derived using a multivariate fit to the observed binding energies. The resulting fit shows a high correlation with a multiple coefficient of determination of 0.79. This indicates that 79% of the variation in the observed binding energies is explained by the electrostatic and hydrophobic terms. The analysis results in a protein dielectric of 8.2 and an atomic solvation parameter of 30 cal/mol A2. As a test, these parameters were used to calculate the binding energies of complexes of chymotrypsin and of leukocyte elastase OMTKY3, as well as three other variants of OMTKY3 bound to SGPB. As these structures were not used for the multivariate fit, they serve as an independent check on the derived parameters. The calculated energies for the three new variants of OMTKY3 are in good agreement with the observed values. However, the binding energies of the other complexes are poorly predicted. This implies that the parameters that were obtained are not transferable. The possible causes for this lack of transferability are discussed.

Purification and biochemical analysis of WprA, a 52-kDa serine protease secreted by B. subtilis as an active complex with its 23-kDa propeptide

The Gram-positive bacterium Bacillus subtilis produces numerous proteases that are secreted to the extracellular milieu, and as strains are generated which lack the more prominent proteases, minor ones become detectable. We have isolated a 52-kDa secreted protease from the protease-deficient strain WB600. It is encoded by the wprA gene which encompasses a signal sequence, a 46-kDa propeptide further processed to 23 kDa, and the 52-kDa mature protease. The 52-kDa and 23-kDa polypeptides were previously detected in cell-wall preparations of a wild-type strain. We have co-purified these proteins from culture supernatant, and confirmed the same N-termini and molecular weights as the membrane-bound species. The WprA protease domain has 28.5% identity to subtilisin A, and like other subtilisins, it displays a broad substrate specificity. WprA and subtilisin A have similar pH profiles, showing optimal activity near pH 7.5 for substrates with Met, Gln, or Lys residues at P1. Using a substrate with Asp at P1, another peak of activity was observed for WprA at pH 5 and at pH 6 for subtilisin A. The pH dependence of some bacterial proteases in their interaction with substrates and inhibitors may be biologically relevant.

Dissecting the energetics of a protein-protein interaction: the binding of ovomucoid third domain to elastase

An understanding of the structural basis for protein-protein interactions, and molecular recognition in general, requires complete characterization of binding energetics. Not only does this include quantification of the changes that occur in all of the thermodynamic parameters upon binding, including the enthalpy, entropy and heat capacity, but a description of how these changes are modulated by environmental conditions, most notably pH. Here, we have investigated the binding of turkey ovomucoid third domain (OMTKY3), a potent serine protease inhibitor, to the serine protease porcine pancreatic elastase (PPE) using isothermal titration calorimetry and structure-based thermodynamic calculations. We find that near neutral pH the binding energetics are influenced by a shift in the pKa of an ionizable group, most likely histidine 57 in the protease active site. Consequently, the observed binding energetics are strongly dependent upon solution conditions. Through a global analysis, the intrinsic energetics of binding have been determined, as have those associated with the pKa shift. The protonation energetics suggest that the drop in pKa is largely due to desolvation of the histidine residue. The resulting deprotonation is necessary for the enzymatic function of elastase. Intrinsically, at 25 degrees C the binding of OMTKY3 to PPE is characterized by an almost negligible enthalpy change, a large positive entropy change, and a large negative heat capacity change. These parameters are consistent with a model of the OMTKY3-PPE complex, which shows a large and significantly apolar protein-protein interface. Thermodynamic calculations based upon changes that occur in polar and apolar solvent-accessible surface area are in very good agreement with the measured intrinsic binding energetics.

The refined crystal structure of the 3C gene product from hepatitis A virus: specific proteinase activity and RNA recognition

The virally encoded 3C proteinases of picornaviruses process the polyprotein produced by the translation of polycistronic viral mRNA. The X-ray crystallographic structure of a catalytically active mutant of the hepatitis A virus (HAV) 3C proteinase (C24S) has been determined. Crystals of this mutant of HAV 3C are triclinic with unit cell dimensions a = 53.6 A, b = 53.5 A, c = 53.2 A, alpha = 99.1 degrees, beta = 129.0 degrees, and gamma = 103.3 degrees. There are two molecules of HAV 3C in the unit cell of this crystal form. The structure has been refined to an R factor of 0.211 (Rfree = 0.265) at 2.0-A resolution. Both molecules fold into the characteristic two-domain structure of the chymotrypsin-like serine proteinases. The active-site and substrate-binding regions are located in a surface groove between the two beta-barrel domains. The catalytic Cys 172 S(gamma) and His 44 N(epsilon2) are separated by 3.9 A; the oxyanion hole adopts the same conformation as that seen in the serine proteinases. The side chain of Asp 84, the residue expected to form the third member of the catalytic triad, is pointed away from the side chain of His 44 and is locked in an ion pair interaction with the epsilon-amino group of Lys 202. A water molecule is hydrogen bonded to His 44 N(delta1). The side-chain phenolic hydroxyl group of Tyr 143 is close to this water and to His 44 N(delta1) and may be negatively charged. The glutamine specificity for P1 residues of substrate cleavage sites is attributed to the presence of a highly conserved His 191 in the S1 pocket. A very unusual environment of two water molecules and a buried glutamate contribute to the imidazole tautomer believed to be important in the P1 specificity. HAV 3C proteinase has the conserved RNA recognition sequence KFRDI located in the interdomain connection loop on the side of the molecule diametrically opposite the proteolytic site. This segment of polypeptide is located between the N- and C-terminal helices, and its conformation results in the formation of a well-defined surface with a strongly charged electrostatic potential. Presumably, this surface of HAV 3C participates in the recognition of the 5' and 3' nontranslated regions of the RNA genome during viral replication.

Binding of amino acid side-chains to S1 cavities of serine proteinases

The P1 or primary specificity residue of standard mechanism canonical protein inhibitors of serine proteinases, inserts into the S1 primary specificity cavity of the cognate enzyme upon enzyme-inhibitor complex formation. Both natural evolution and protein engineering often change the P1 residue to greatly alter the specificity and the binding strength. To systematize such results we have obtained all 20 coded P1 variants of one such inhibitor, turkey ovomucoid third domain, by recombinant DNA technology. The variants were extensively characterized. The association equilibrium constants were measured at pH 8.30, 21 (+/-2) degrees C, for interaction of these variants with six well characterized serine proteinases with hydrophobic S1, cavities. The enzyme names are followed by the best, worst and most specific coded residue for each. Bovine chymotrypsin A alpha (Tyr, Pro, Trp), porcine pancreatic elastase (Leu/Ala, Arg, Ala), subtilisin Carlsberg (Cys, Pro, Glu), Streptomyces griseus proteinase A (Cys, Pro, Leu) and B (Cys, Pro, Lys) and human leukocyte elastase (Ile, Asp, Ile). The data set was merged with Ka values for five non-coded variants at P1 of turkey ovomucoid third domain obtained in our laboratory by enzymatic semisynthesis. The ratios of the highest to the lowest Ka for each of the six enzymes range from 10(6) to 10(8). The dominant force for binding to these pockets is the hydrophobic interaction. Excess steric bulk (too large for the pocket), awkward shape (Pro, Val and Ile), polarity (Ser) oppose interaction. Ionic charges, especially negative charges on Glu- and Asp- are strongly unfavorable. The Pearson pro duct moment correlations for all the 15 enzyme pairs were calculated. We suggest that these may serve as a quantitative description of the specificity of the enzymes at P1. The sets of Streptomyces griseus proteinases A and B and of the two elastases are strongly positively correlated. Strikingly, chymotrypsin and pancreatic elastase are negatively correlated (-0.10). Such correlations can be usefully extended to many other enzymes and to many other binding pockets to provide a general measure of pocket binding specificity.

Interscaffolding additivity. Association of P1 variants of eglin c and of turkey ovomucoid third domain with serine proteinases

Standard mechanism protein inhibitors of serine proteinases share a common mechanism of interaction with their cognate enzymes. The P1 residue of the inhibitor interacts with the enzyme in a substrate-like manner. Its side chain becomes imbedded in the S1 cavity of the enzyme. The nature of P1, the primary specificity residue, greatly affects the strength and specificity of the enzyme inhibitor association. In canonical inhibitors, residues P4-P2'(P3'), where P1-P1' is the reactive site, share a common main chain conformation that does not change on complex formation. The remainder of the inhibitor's structure, the scaffolding, is not always common. Instead, there are at least 20 inhibitor families, each with a different scaffolding. In this paper, we ask whether the differences in standard free energy of association of enzyme-inhibitor complexes upon P1 mutations are independent of the nature of the scaffolding. We have already reported on 25 P1 variants of turkey ovomucoid third domain, a member of the Kazal inhibitor family, interacting with six different serine proteinases. Here, we report on seven different P1 variants of eglin c, a potato I family member, interacting with the same six serine proteinases under the same conditions. The differences in standard free energy on P1 mutations in the eglin c system agree very well, when P1 Pro is omitted. Complete agreement indicates that these P1 residues are interscaffolding additive. This is consistent with the superimposition of the high-resolution structures of eglin c and of turkey ovomucoid third domain with chymotrypsin. In both cases, the P1 Leu side chain is similarly oriented in almost indistinguishable specificity pockets of the enzyme.