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

Computational analysis of propeptide-containing proteins and prediction of their post-cleavage conformation changes

A propeptide is removed from a precursor protein to generate its active or mature form. Propeptides play essential roles in protein folding, transportation, and activation and are present in about 2.3% of reviewed proteins in the UniProt database. They are often found in secreted or membrane-bound proteins including proteolytic enzymes, hormones, and toxins. We identified a variety of globular and nonglobular Pfam domains in protein sequences designated as propeptides, some of which form intramolecular interactions with other domains in the mature proteins. Propeptide-containing enzymes mostly function as proteases, as they are depleted in other enzyme classes such as hydrolases acting on DNA and RNA, isomerases, and lyases. We applied AlphaFold to generate structural models for over 7000 proteins with propeptides having no less than 20 residues. Analysis of residue contacts in these models revealed conformational changes for over 300 proteins before and after the cleavage of the propeptide. Examples of conformation change occur in several classes of proteolytic enzymes in the families of subtilisins, trypsins, aspartyl proteases, and thermolysin-like metalloproteases. In most of the observed cases, cleavage of the propeptide releases the constraints imposed by the covalent bond between the propeptide and the mature protein, and cleavage enables stronger interactions between the propeptide and the mature protein. These findings suggest that post-cleavage propeptides could play critical roles in regulating the activity of mature proteins.

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.

Biochemical and structural characterization of SplD protease from Staphylococcus aureus

Staphylococcus aureus is a dangerous human pathogen. A number of the proteins secreted by this bacterium are implicated in its virulence, but many of the components of its secretome are poorly characterized. Strains of S. aureus can produce up to six homologous extracellular serine proteases grouped in a single spl operon. Although the SplA, SplB, and SplC proteases have been thoroughly characterized, the properties of the other three enzymes have not yet been investigated. Here, we describe the biochemical and structural characteristics of the SplD protease. The active enzyme was produced in an Escherichia coli recombinant system and purified to homogeneity. P1 substrate specificity was determined using a combinatorial library of synthetic peptide substrates showing exclusive preference for threonine, serine, leucine, isoleucine, alanine, and valine. To further determine the specificity of SplD, we used high-throughput synthetic peptide and cell surface protein display methods. The results not only confirmed SplD preference for a P1 residue, but also provided insight into the specificity of individual primed- and non-primed substrate-binding subsites. The analyses revealed a surprisingly narrow specificity of the protease, which recognized five consecutive residues (P4-P3-P2-P1-P1’) with a consensus motif of R-(Y/W)-(P/L)-(T/L/I/V)↓S. To understand the molecular basis of the strict substrate specificity, we crystallized the enzyme in two different conditions, and refined the structures at resolutions of 1.56 Å and 2.1 Å. Molecular modeling and mutagenesis studies allowed us to define a consensus model of substrate binding, and illustrated the molecular mechanism of protease specificity.

SKEMPI: a Structural Kinetic and Energetic database of Mutant Protein Interactions and its use in empirical models

MOTIVATION

Empirical models for the prediction of how changes in sequence alter protein-protein binding kinetics and thermodynamics can garner insights into many aspects of molecular biology. However, such models require empirical training data and proper validation before they can be widely applied. Previous databases contained few stabilizing mutations and no discussion of their inherent biases or how this impacts model construction or validation.

RESULTS

We present SKEMPI, a database of 3047 binding free energy changes upon mutation assembled from the scientific literature, for protein-protein heterodimeric complexes with experimentally determined structures. This represents over four times more data than previously collected. Changes in 713 association and dissociation rates and 127 enthalpies and entropies were also recorded. The existence of biases towards specific mutations, residues, interfaces, proteins and protein families is discussed in the context of how the data can be used to construct predictive models. Finally, a cross-validation scheme is presented which is capable of estimating the efficacy of derived models on future data in which these biases are not present.

AVAILABILITY

The database is available online at http://life.bsc.es/pid/mutation_database/.

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.

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.

Structural Insights into the Non-additivity Effects in the Sequence-to-Reactivity Algorithm for Serine Peptidases and their Inhibitors

Sequence-to-reactivity algorithms (SRAs) for proteins have the potential of being broadly applied in molecular design. Recently, Laskowski et al. have reported an additivity-based SRA that accurately predicts most of the standard free energy changes of association for variants of turkey ovomucoid third domain (OMTKY3) with six serine peptidases, one of which is streptogrisin B (commonly known as Streptomyces griseus peptidase B, SGPB). Non-additivity effects for residues 18I and 32I, and for residues 20I and 32I of OMTKY3 occurred when the associations with SGPB were predicted using the SRA. To elucidate precisely the mechanics of these non-additivity effects in structural terms, we have determined the crystal structures of the unbound OMTKY3 (with Gly32I as in the wild-type amino acid sequence) at a resolution of 1.16 A, the unbound Ala32I variant of OMTKY3 at a resolution of 1.23 A, and the SGPB:OMTKY3-Ala32I complex (equilibrium association constant K(a)=7.1x10(9) M(-1) at 21(+/-2) C degrees, pH 8.3) at a resolution of 1.70 A. Extensive comparisons with the crystal structure of the unbound OMTKY3 confirm our understanding of some previously addressed non-additivity effects. Unexpectedly, SGPB and OMTKY3-Ala32I form a 1:2 complex in the crystal. Comparison with the SGPB:OMTKY3 complex shows a conformational change in the SGPB:OMTKY3-Ala32I complex, resulting from a hinged rigid-body rotation of the inhibitor caused by the steric hindrance between the methyl group of Ala32IA of the inhibitor and Pro192BE of the peptidase. This perturbs the interactions among residues 18I, 20I, 32I and 36I of the inhibitor, probably resulting in the above non-additivity effects. This conformational change also introduces residue 10I as an additional hyper-variable contact residue to the SRA.

Defining residues involved in human rhinovirus 2A proteinase substrate recognition

The 2A proteinase (2A(pro)) of human rhinoviruses (HRVs) initiates proteolytic processing by cleaving between the C-terminus of VP1 and its own N-terminus. It subsequently cleaves the host protein eIF4GI. HRV2 and HRV14 2A(pro) cleave at IITTA *GPSD and DIKSY *GLGP on their respective polyproteins. The HRV2 2A(pro) cleavage site on eIF4GI is TLSTR *GPPR. We show that HRV2 2A(pro) can self-process at the eIF4GI cleavage sequence whereas HRV14 2A(pro) cannot, due to the presence of the arginine residue at P1. The mutations A104C or A104S in HRV14 2A(pro) restored cleavage when arginine was present at P1, although not to wild-type levels. These experiments define residues which determine substrate recognition in rhinoviral 2A(pro).

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.

Structural consequences of accommodation of four non-cognate amino acid residues in the S1 pocket of bovine trypsin and chymotrypsin

Crystal structures of P1 Gly, Val, Leu and Phe bovine pancreatic trypsin inhibitor (BPTI) variants in complex with two serine proteinases, bovine trypsin and chymotrypsin, have been determined. The association constants for the four mutants with the two enzymes show that the enlargement of the volume of the P1 residue is accompanied by an increase of the binding energy, which is more pronounced for bovine chymotrypsin. Since the conformation of the P1 side-chains in the two S1 pockets is very similar, we suggest that the difference in DeltaG values between the enzymes must arise from the more polar environment of the S1 site of trypsin. This results mainly from the substitutions of Met192 and Ser189 observed in chymotrypsin with Gln192 and Asp189 present in trypsin. The more polar interior of the S1 site of trypsin is reflected by a much higher order of the solvent network in the empty pocket of the enzyme, as is observed in the complexes of the two enzymes with the P1 Gly BPTI variant. The more optimal binding of the large hydrophobic P1 residues by chymotrypsin is also reflected by shrinkage of the S1 pocket upon the accommodation of the cognate residues of this enzyme. Conversely, the S1 pocket of trypsin expands upon binding of such side-chains, possibly to avoid interaction with the polar residues of the walls. Further differentiation between the two enzymes is achieved by small differences in the shape of the S1 sites, resulting in an unequal steric hindrance of some of the side-chains, as observed for the gamma-branched P1 Leu variant of BPTI, which is much more favored by bovine chymotrypsin than trypsin. Analysis of the discrimination of beta-branched residues by trypsin and chymotrypsin is based on the complexes with the P1 Val BPTI variant. Steric repulsion of the P1 Val residue by the walls of the S1 pocket of both enzymes prevents the P1 Val side-chain from adopting the most optimal chi1 value.

Predicting the reactivity of proteins from their sequence alone: Kazal family of protein inhibitors of serine proteinases

An additivity-based sequence to reactivity algorithm for the interaction of members of the Kazal family of protein inhibitors with six selected serine proteinases is described. Ten consensus variable contact positions in the inhibitor were identified, and the 19 possible variants at each of these positions were expressed. The free energies of interaction of these variants and the wild type were measured. For an additive system, this data set allows for the calculation of all possible sequences, subject to some restrictions. The algorithm was extensively tested. It is exceptionally fast so that all possible sequences can be predicted. The strongest, the most specific possible, and the least specific inhibitors were designed, and an evolutionary problem was solved.

Predicting the reactivity of proteins from their sequence alone: Kazal family of protein inhibitors of serine proteinases

An additivity-based sequence to reactivity algorithm for the interaction of members of the Kazal family of protein inhibitors with six selected serine proteinases is described. Ten consensus variable contact positions in the inhibitor were identified, and the 19 possible variants at each of these positions were expressed. The free energies of interaction of these variants and the wild type were measured. For an additive system, this data set allows for the calculation of all possible sequences, subject to some restrictions. The algorithm was extensively tested. It is exceptionally fast so that all possible sequences can be predicted. The strongest, the most specific possible, and the least specific inhibitors were designed, and an evolutionary problem was solved.