Four-body potentials reveal protein-specific correlations to stability changes caused by hydrophobic core mutations

Mutational experiments show how changes in the hydrophobic cores of proteins affect their stabilities. Here, we estimate these effects computationally, using four-body likelihood potentials obtained by simplicial neighborhood analysis of protein packing (SNAPP). In this procedure, the volume of a known protein structure is tiled with tetrahedra having the center of mass of one amino acid side-chain at each vertex. Log-likelihoods are computed for the 8855 possible tetrahedra with equivalent compositions from structural databases and amino acid frequencies. The sum of these four-body potentials for tetrahedra present in a given protein yields the SNAPP score. Mutations change this sum by changing the compositions of tetrahedra containing the mutated residue and their related potentials. Linear correlation coefficients between experimental mutational stability changes, Delta(DeltaG(unfold)), and those based on SNAPP scoring range from 0.70 to 0.94 for hydrophobic core mutations in five different proteins. Accurate predictions for the effects of hydrophobic core mutations can therefore be obtained by virtual mutagenesis, based on changes to the total SNAPP likelihood potential. Significantly, slopes of the relation between Delta(DeltaG(unfold)) and DeltaSNAPP for different proteins are statistically distinct, and we show that these protein-specific effects can be estimated using the average SNAPP score per residue, which is readily derived from the analysis itself. This result enhances the predictive value of statistical potentials and supports previous suggestions that "comparable" mutations in different proteins may lead to different Delta(DeltaG(unfold)) values because of differences in their flexibility and/or conformational entropy.     

Urea-induced conformational changes in cold- and heat-denatured states of a protein, Streptomyces subtilisin inhibitor.

Streptomyces subtilisin inhibitor (SSI) is known to exist in at least two distinct denatured states, cold-denatured (D') and heat-denatured (D) under acidic conditions. In the present work, we investigated the manner how increasing urea concentration from 0 to 8 M changes the polypeptide chain conformation of SSI that exists initially in the D' and D states as well as in the native state (N), in terms of the secondary structure, the tertiary structure, and the chain form, based on the results of the experiments using circular dichroism (CD), small-angle X-ray scattering (SAXS) and 1H-NMR spectroscopy. Our results indicate that the urea-induced conformational transitions of SSI under typical conditions of D' (pH 1.8, 3 degrees C) occur at least in two steps. In the urea concentration range of 0-2 M (step 1), a cooperative destruction of the tertiary structure occurs, resulting in a mildly denatured state (DU), which may still contain a little amount of secondary structures. In the concentration range of 2-4 M urea (step 2), the DU state gradually loses its residual secondary structure, and increases the radius of gyration nearly to a maximum value. At 4 M urea, the polypeptide chain is highly disordered with highly mobile side chains. Increasing the urea concentration up to 8 M probably results in the more highly denatured or alternatively the stiffer chain conformations. The conformational transition starting from the N state proceeds essentially the same way as in the above scheme in which D' is replaced with N. The conformational transition starting from the D state lacks step 1 because the D state contains no tertiary structures and is similar to the DU state. The fact that similar conformations are reached at urea concentrations above 2 M from different conformations of D', D, and N indicates that the effect of urea dominates in determining the polypeptide conformation of SSI in the denatured states rather than the pH and temperature.     

Hydrogen exchange kinetics of core peptide protons in Streptomyces subtilisin inhibitor

The hydrogen isotope exchange kinetics of the 10 slowest exchanging resonances in the 1H NMR spectrum of Streptomyces subtilisin inhibitor (SSI) have been determined at pH 7-11 and 30-60 degrees C. These resonances are assigned to peptide amide protons in the beta-sheet core that comprises the extensive protein-protein interface of the tightly bound SSI dimer. The core protons are atypical in that their exchange rates are orders of magnitude slower than those for all other SSI protons. When they do exchange at temperatures greater than 50 degrees C, they do so as a set and with a very high temperature coefficient. The pH dependence of the exchange rate constants is also atypical. Exchange rates are approximately first order in hydroxyl ion dependence at pH less than 8.5 and greater than 9.5 and pH independent between pH 8.5 and 9.5. The pH dependence and temperature dependence of the SSI proton exchange rates are interpreted by the two-process model [Woodward, C. K., & Hilton, B. D. (1980) Biophys. J. 32, 561-575]. The results suggest that in the average solution structure of SSI, an unusual mobility of secondary structural elements at the protein surface is, in a sense, compensated by an unusual rigidity and inaccessibility of the beta-sheet core at the dimer interface.     

Inhibition of subtilisin BPN' by reaction site P1 mutants of Streptomyces subtilisin inhibitor

It has been shown that the P1 site (the center of the reactive site) of protease inhibitors corresponds to the specificity of the cognate protease, and consequently specificity of Streptomyces subtilisin inhibitor (SSI) can be altered by substitution of a single amino acid at the P1 site. In this paper, to investigate whether similar correlation between inhibitory activity of mutated SSI and substrate preference of protease is observed for subtilisin BPN', which has broad substrate specificity, a complete set of mutants of SSI at the reaction site P1 (position 73) was constructed by cassette and site-directed mutagenesis and their inhibitory activities toward subtilisin BPN' were measured. Mutated SSIs which have a polar (Ser, Thr, Gln, Asn), basic (Lys, Arg), or aromatic amino acid (Tyr, Phe, Trp, His), or Ala or Leu, at the P1 site showed almost the same strong inhibitory activity toward subtilisin as the wild type (Met) SSI. However, the inhibitory activity of SSI variants with an acidic (Glu, Asp), or a beta-branched aliphatic amino acid (Val, Ile), or Gly or Pro, at P1 was decreased. The values of the inhibitor constant (Ki) of mutated SSIs toward subtilisin BPN' were consistent with the substrate preference of subtilisin BPN'. A linear correlation was observed between log(1/Ki) of mutated SSIs and log(1/Km) of synthetic substrates. These results demonstrate that the inhibitory activities of P1 site mutants of SSI are linearly related to the substrate preference of subtilisin BPN', and indicate that the binding mode of the inhibitors with the protease may be similar to that of substrates, as in the case of trypsin and chymotrypsin.(ABSTRACT TRUNCATED AT 250 WORDS)     

Kinetic studies on the binding of Streptomyces subtilisin inhibitor with subtilisin BPN′

The binding mechanism of Streptomyces subtilisin inhibitor and subtilisin BPN′ was studied kinetically with the stopped-flow method by monitoring the protein fluorescence increase due to complex formation. In the lower concentration range of proteins, the reaction followed the second-order kinetics. The pH dependence of the apparent second-order rate constant, k_on, suggested the involvement of the two ionizable groups of pK_a of 7.8 and 10 in the binding. The activation parameters were calculated from the temperature dependence of the apparent second-order rate constants. The value of the apparent activation energy (E_A = 39.7 kJ · mol^?1, 9.50 kcal · mol^?1) and insensitivity of k_on to the viscosity of the medium suggest that the binding is not a simple diffusion-controlled bimolecular association. Further studies with a much broader range of protein concentrations have revealed that the reaction tends to approach first-order kinetics as the inhibitor concentration increases. The binding reaction is, therefore, reconcilable with a two-step mechanism, in which a fast bimolecular association is followed by a slow unimolecular isomerization step; the dissociation constant of the first step, K_L, is estimated to be 1.2 × 10^?4m and the rate constant of the second step, k₊₂, to be 770 s^?1. It was also found that the increase of tryptophan fluorescence due to the complex formation occurs solely in the rate-determining unimolecular process.     

A protease inhibitor produced by Streptomyces lividans 66 exhibits inhibitory activities toward both subtilisin BPN' and trypsin.

A proteinaceous protease inhibitor was isolated from the culture broth of Streptomyces lividans 66 by a series of purification steps (salting out by ammonium sulfate, ion-exchange chromatography on DEAE-cellulose, hydrophobic chromatography on Phenyl-Sepharose, and gel-filtration on Sephacryl S-200), and was named S. lividans protease inhibitor (SLPI). The purified SLPI existed in a dimeric form consisting of two identical subunits, each of which was composed of 107 amino acids. SLPI exhibited strong inhibitory activity toward subtilisin BPN'. These features were similar to those of protein protease inhibitors produced by other Streptomyces (SSI family inhibitor). In addition, SLPI was capable of inhibiting trypsin with an inhibitor constant (Ki) of about 10(-9) M. The primary structure of SLPI and location of two disulfide bridges were homologous to those of the other serine protease inhibitors of Streptomyces. The reactive site of SLPI was found to be Arg67-Glu68 from the sequence analysis of cleaved SLPI which was produced by acidification of subtilisin-SLPI complex. An Arg residue at the P1 site was consistent with the trypsin-inhibitory property of SLPI. Sequence comparison with other members of the SSI family revealed that amino acid replacements in SLPI were mainly localized on the surface of the SLPI molecule, and many of the amino acid residues in beta-sheets and hydrophobic core were well conserved.     

Stabilizing mutations and calcium-dependent stability of subtilisin

Stability is a property of subtilisin which has proven particularly amenable to enhancement via random mutagenesis and screening, yet the effects of most stabilizing mutations are not understood in structural and energetic detail. This paper seeks to explain the longstanding observation that stabilizing mutations are usually calcium-dependent in their stabilizing effect, irrespective of their proximity to the calcium binding sites. Stabilizing mutations in subtilisin fall into one of three classes. The largest class of mutations stabilize only in the presence of excess calcium. A smaller number of mutations stabilize independently of [calcium], and a few mutations stabilize only in the presence of chelating agents, such as EDTA. This study compares the effects of mutations from each class when introduced into subtilisin BPN' and two calcium-free versions of subtilisin. The calcium-dependent effects of mutations can be explained by considering subtilisin to be in conformational equilibrium between two structurally similar but energetically distinct states: N and N*. The equilibrium from the N* to the N state can be altered either by calcium binding to site A or by mutation. Mutations which stabilize only in the presence of calcium stabilize the N state relative to N*. Mutations which stabilize only in the presence of chelants stabilize the N* state relative to N. As a byproduct of this analysis, we have developed a hyperstable variant of subtilisin whose inactivation at high temperature in the presence of EDTA is 10(5) times slower than wild-type subtilisin.     

Structure and fluctuation of a Streptomyces subtilisin inhibitor.

The kinetics of the hydrogen-deuterium exchange reaction in a subtilisin inhibitor from Streptomyces albogriseolus has been examined by infrared absorption measurement in aqueous solutions at various pH values and temperatures. In the analysis of each piece of kinetic data, it was assumed that the total 104 peptide hydrogen atoms are classified into three kinetic classes A, B1, and B2, and that the sizes of these classes are 72, 15, and 17, respectively at every pH and at every temperature examined. On the basis of the peak position determined for the amide II band in each stage of the exchange reaction, an approximate assignment was suggested of the A, B1 and B2 respectively to an unordered structure, a beta-structure,and an alpha-helical structure in the molecule. This assignment was supported by infrared absorption measurement of a film of this protein and by circular dichroic study of the solutions. On the basis of the temperature effect on the hydrogen-exchange rate constants and on the basis of ultraviolet absorption study in the higher temperature region (40 to 90 degrees C), a discussion has been made on the nature of the fluctuation of the molecular structure of this protein.     

Thermodynamics of the binding of Streptomyces subtilisin inhibitor to alpha-chymotrypsin

The binding of Streptomyces subtilisin inhibitor (SSI) to alpha-chymotrypsin (CT) (EC 3.4.21.1) was studied by isothermal and differential scanning calorimetry at pH 7.0. Thermodynamic quantities for the binding of SSI to the enzyme were derived as functions of temperature from binding constants (S. Matsumori, B. Tonomura, and K. Hiromi, private communication) and isothermal calorimetric experiments at 5-30 degrees C. At 25 degrees C, the values are delta G degrees b = -29.9 kJ mol-1, delta Hb = +18.7 (+/- 1.3) kJ mol-1, delta S degrees b = +0.16 kJ K-1 mol-1, and delta C p,b = -1.08 (+/- 0.11) kJ mol-1. The binding of SSI to CT is weak compared with its binding to subtilisin [Uehara, Y., Tonomura, B., & Hiromi, K. (1978) J. Biochem. (Tokyo) 84, 1195-1202; Takahashi, K., & Fukada, H. (1985) Biochemistry 24, 297-300]. This difference is due primarily to a less favorable enthalpy change in the formation of the complex with CT. The hydrophobic effect is presumably the major source of the entropy and heat capacity changes which accompany the binding process. The unfolding temperature of the complex is about 7 degrees C higher than that of the free enzyme. The enthalpy and the heat capacity changes for the unfolding of CT were found to be 814 kJ mol-1 and 17.3 kJ K-1 mol-1 at 49 degrees C. The same quantities for the unfolding of the SSI-CT complex are 1183 kJ mol-1 and 39.2 kJ K-1 mol-1 at 57 degrees C.     

Refined crystal structure of the complex of subtilisin BPN' and Streptomyces subtilisin inhibitor at 1.8 A resolution

The crystal structure of subtilisin BPN' complexed with a proteinaceous inhibitor SSI (Streptomyces subtilisin inhibitor) was refined at 1.8 A resolution to an R-factor of 0.177 with a root-mean-square deviation from ideal bond lengths of 0.014 A. The work finally established that the SSI-subtilisin complex is a Michaelis complex with a distance between the O gamma of active Ser221 and the carbonyl carbon of the scissile peptide bond being an intermediate value between a covalent bond and a van der Waals' contact, 2.7 A. This feature, as well as the geometry of the catalytic triad and the oxyanion hole, is coincident with that found in other highly refined crystal structures of the complex of subtilisin Novo, subtilisin Carlsberg, bovine trypsin or Streptomyces griseus protease B with their proteinaceous inhibitors. The enzyme-inhibitor beta-sheet interaction is composed of two separate parts: that between the P1-P3 residues of SSI and the 125-127 chain segment (the "S1-3 site") of subtilisin and that between the P4-P6 residues of SSI and th 102-104 chain segment (the "S4-6 site") of subtilisin. The latter beta-interaction is unique to subtilisin. In contrast, the beta-sheet interaction previously found in the complex of subtilisin Novo and chymotrypsin inhibitor 2 or in the complex of subtilisin Carlsberg and Eglin C is distinct from the present complex in that the two types of beta-interactions are not separate. As for the flexibility of the molecules comprising the present complex, the following observations were made by comparing the B-factors for free and complexed SSI and comparing those for free and complexed subtilisin BPN'. The rigidification of the component molecules upon complex formation occurs in a very localized region: in SSI, the "primary" and "secondary" contact regions and the flanking region; in subtilisin BPN', the S1-3 and S4-6 sites and the flanking region.     

Calorimetric studies of the binding of Streptomyces subtilisin inhibitor to subtilisin of Bacillus subtilis strain N'

The binding of Streptomyces subtilisin inhibitor (SSI) to subtilisin of Bacillus subtilis strain N' (subtilisin BPN', EC 3.4.21.14) was studied by isothermal calorimetry at pH 7.0 and at various temperatures ranging from 5 to 30 degrees C. Thermodynamic quantities for the binding reaction were derived as a function of temperature by combining the data reported for the dissociation constant with the present calorimetric results. At 25 degrees C, the values are delta G degrees = -57.9 kJ mol-1, delta H = -19.8 kJ mol-1, delta S degree = 0.13 kJ K-1 mol-1, and delta Cp = -1.02 kJ K-1 mol-1. The entropy and the heat capacity changes are discussed in terms of the contributions from the changes in vibrational modes and in hydrophobic interactions.     

Thermal stability and conformational changes of transglutaminase from a newly isolated Streptomyces hygroscopicus

Thermal stability and conformational changes of transglutaminase (TGase) from a newly isolated Streptomyces hygroscopicus were investigated in this study. The inactivation kinetics of the microbial transglutaminase (MTGase) was fitted using one-step inactivation model. It was much more stable under 40 degrees C. The half-lives for the MTGase at 50 degrees C and 60 degrees C were only 20 min and 8 min, respectively. Spectroscopic studies of the enzyme suggested conformational transition from ordered secondary structural elements (alpha/beta-protein) to unordered structure during thermal denaturation. Some polyols could improve the thermal stability of the enzyme. Among the polyols examined, the prolonged half-lives of 40 min at 50 degrees C and 20 min at 60 degrees C were gained by adding 10% glycerol. The results of differential scanning calorimetric (DSC) analysis showed a distinct transition peak with a significant greater Tm and DeltaH for the MTGase mixed with polyols in comparison with the control, which indicated that the polyols could maintain the natural structure of the enzyme to some extent. The SDS-PAGE electrophoresis of cross-linked casein confirmed that the stabilizers could protect the MTGase from thermal denaturation.     

Complex of subtilisin BPN' with Streptomyces subtilisin inhibitor. Complex formation concomitant with change in reducibility of disulfide bonds in the inhibitor

Structure of the complex of Streptomyces subtilisin inhibitor (SSI) with subtilisin BPN' was studied by examining the thermal denaturation and reducibility of disulfide bonds. The denaturation temperature of the complex was significantly higher than that of the enzyme. Two disulfide bonds localized in the inhibitor side were completely reduced in the complex, whereas only one of them was reduced in the free SSI. Gel filtration of the reduced complex solution showed clearly that the main products of reduction of the complex were two peptide fragments of SSI divided at the active site. The resistive disulfide bond in the complexed inhibitor became accessible as a result of a large conformational change due to splitting of the half-reduced inhibitor.     

High-level expression in Streptomyces lividans 66 of a gene encoding Streptomyces subtilisin inhibitor from Streptomyces albogriseolus S-3253

A secretory expression system for Streptomyces subtilisin inhibitor (SSI) was established in a heterologous host, Streptomyces lividans 66, by introducing the 1.8-kbp BglII/SalI fragment containing SSI gene into the Streptomyces multicopy vector, pIJ 702. The expression of SSI did not depend on the orientation of the 1.8-kbp BglII/SalI fragment or on the promoter for tyrosinase gene (mel) in pIJ 702, which suggested that this fragment also carries the SSI promoter. The expressed SSI in S.lividans 66 was secreted into the culture medium in a large amount, as observed with the original strain, S. albogriseolus S-3253. Amino acid sequence analysis showed that the SSI secreted from S. lividans 66 contained three additional amino acid residues in the NH2-terminal region. The inhibitory activity toward subtilisin BPN' and the antigenic activity of the SSI secreted from S. lividans 66 were found to be identical with those of authentic SSI.     

The states of tyrosyl and tryptophyl residues in a protein proteinase inhibitor (Streptomyces subtilisin inhibitor

The states of tyrosyl and tryptophyl residues of a dimeric protein proteinase inhibitor, Streptomyces subtilisin inhibitor (Sato, S & Murao, S. (1973), Agric. Biol. Chem. 37, 1067) were studies by solvent perturbation difference spectroscopy with methanol, ethylene glycol, polyethylene glycol, and deuterium oxide as perturbants, and by spectrophotometric titration at alkaline pH. It appeared that all three tyrosyl residues per monomer of the inhibitor were exposed on the surface of the molecule, and their apparent pK values were estimated separately to be 9.58, 11.10, and 12.42. The single tryptophyl residue per monomer of the inhibitor appeared to be partially buried in the protein molecule.     

Crystal structure of the complex of subtilisin BPN' with its protein inhibitor Streptomyces subtilisin inhibitor. The structure at 4.3 Angstroms resolution

The crystal structure of the complex of subtilisin BPN′ (EC 3.4.21.14) with its protein inhibitor (Streptomyces subtilisin inhibitor) was solved at 4.3 Å resolution, thus establishing the following. (1) Two subtilisin BPN′ molecules (2E) associate with one dimeric inhibitor molecule (I₂) to form the complex molecule E₂I₂. (2) The conformation of neither the inhibitor nor subtilisin BPN′ undergoes any detectable change at this resolution upon complex formation. (3) The inhibitor binds to subtilisin to form an antiparallel β-sheet, as in the case of trypsin/ trypsin inhibitor complexes. (4) The scissible bond of the inhibitor is between Met73′ and Val74′, as proposed earlier (Ikenaka et al., 1974). (5) The protein inhibitor and the substrates bind to subtilisin BPN′ in essentially the same way.     

Catalytic and conformational changes induced by limited subtilisin cleavage of bovine carboxypeptidase A.

Limited proteolysis of carboxypeptidase A from bovine pancreas with subtilisin Carlsberg generates a stable intermediate, carboxypeptidase S, whose esterase and peptidase activities are increased and decreased, respectively, under standard assay conditions. Carboxypeptidase S was isolated by affinity chromatography. Sequence analysis shows that it is cleaved solely at the Ala154-Gly155 bond. Its enzymatic properties were determined under stopped-flow conditions with Dns-Gly-Ala-Phe and its ester analogue Dns-Gly-Ala-OPhe. For both substrates, the Km values are increased 30-40-fold. The kcat value for peptide hydrolysis is virtually unaffected whereas that for ester hydrolysis is increased 10-fold. The magnitude of the Km effect is equivalent to a loss of 9 kJ/mol of binding energy and likely reflects a disruption of the network of hydrogen bonds that links Tyr-248 and Arg-145 to the backbone carbonyls of Ala-154 and Gly-155. The difference in kcat effects for the two substrate classes is related to differences in the chemical nature of the rate-determining step. Product release is rate determining for catalytic hydrolysis of ester substrates, and hence, the increase in kcat indicates that dissociation of products is facilitated as a result of the Ala154-Gly155 bond scission. The changes in enzymatic activity accompanying limited proteolysis are due to conformational alterations in the vicinity of the active center of the molecule. The affinity of a monoclonal antibody, mAb 100, directed toward the antigenic determinant located between residues 209 and 218 in carboxypeptidase A is diminished considerably for carboxypeptidase S.(ABSTRACT TRUNCATED AT 250 WORDS)     

Secondary structures and structural fluctuation in a dimeric protein, Streptomyces subtilisin inhibitor.

Based on the nuclear magnetic resonance assignments of a dimeric protein, Streptomyces subtilisin inhibitor (SSI), microscopic details of secondary structures in solution have been elucidated. The chemical shift index of C(alpha) signals, together with information on the hydrogen exchange rates of the backbone amide protons, were used to identify secondary structures. The locations of these secondary structures were found to be different in some critical points from those determined earlier by X-ray crystallography of the crystal. Notably, the beta3 strand is completely missing and the alpha2 helix is extended toward the C-terminus. Furthermore, hydrogen exchange experiments of individual peptide NH protons under strongly folding conditions revealed mechanisms of global and local structural fluctuation within the dimeric structure. It has been suggested that the global fluctuation of the monomeric unit occurs without affecting the accompanying monomer, in contrast to the equilibrium thermal unfolding, which is cooperative. Higher protection against hydrogen exchange for residues in part of the beta4 strand implies that this region might serve as a folding core.     

pH-induced conformational change in an alpha-helical coiled-coil is controlled by His residues in the hydrophobic core

An alpha-helical coiled-coil structure is one of the basic structural units in proteins. Hydrophilic residues at the hydrophobic positions in the coiled-coil structure play important roles in structures and functions of natural proteins. We reported here a peptide that formed a triple stranded alpha-helical coiled-coil showing the pH-dependent structural change. The peptide was designed to have two His residues at the hydrophobic positions of the center of the coiled-coil structure. The peptide folded into a triple stranded coiled-coil at neutral pH, while it unfolded at acidic pH. This construct is useful to create a protein that the structure or function is controlled by pH.