Inverse electrostatic effect: electrostatic repulsion in the unfolded state stabilizes a leucine zipper

The pH-dependent stability of a protein is strongly affected by electrostatic interactions between ionizable residues in the folded as well as unfolded state. Here we characterize the individual contributions of charged Glu and His residues to stability and determine the NMR structure of the designed, heterodimeric leucine zipper AB consisting of an acidic A chain and a basic B chain. Thermodynamic parameters are compared with those of the homologous leucine zipper AB(SS) in which the A and B chains are disulfide-linked. NMR structures of AB based on (1)H NMR data collected at 600 MHz converge, and formation of the same six interchain salt bridges found previously in disulfide-linked AB(SS) [Marti, D. N., and Bosshard, H. R. (2003) J. Mol. Biol. 330, 621-637] is indicated. While the structures of AB and AB(SS) are very similar, their pH-dependent relative stabilities are strikingly different. The stability of AB peaks at pH approximately 4.5 and is higher at pH 8 than at pH 2. In contrast, AB(SS) is most stable at acidic pH where no interhelical salt bridges are formed. The different energetic contributions of charged Glu and His residues to stability of the two coiled coil structures were evaluated from pK(a) shifts induced by folding. The six charged Glu residues involved in salt bridges stabilize leucine zipper AB by 4.5 kJ/mol yet destabilize disulfide-linked AB(SS) by -1.1 kJ/mol. Two non-ion-paired Glu charges destabilize AB by only -1.8 kJ/mol but AB(SS) by -5.6 kJ/mol. The higher relative stability of AB at neutral pH is not caused by more favorable electrostatic interactions in the folded leucine zipper. It is due mainly to unfavorable electrostatic interactions in the unfolded A and B chains and may therefore be called an inverse electrostatic effect. This study illustrates the importance of residual interactions in the unfolded state and how the energetics of the unfolded state affect the stability of the folded protein.     

Salt bridges destabilize a leucine zipper designed for maximized ion pairing between helices

Interhelical salt bridges are common in leucine zippers and are thought to stabilize the coiled coil conformation. Here we present a detailed thermodynamic investigation of the designed, disulfide-linked leucine zipper AB(SS) whose high-resolution NMR structure shows six interhelical ion pairs between heptad positions g of one helix and e' of the other helix but no ion pairing within single helices. The average pK(a) value of the Glu side chain carboxyl groups of AB(SS) is slightly higher than the pK(a) of a freely accessible Glu in an unfolded peptide [Marti, D. N., Jelesarov, I., and Bosshard, H. R. (2000) Biochemistry 39, 12804-12818]. This indicates that the salt bridges are destabilizing, a prediction we now have confirmed by determining the pH +/- stability profile of AB(SS). Circular dichroism-monitored unfolding by urea and by heating and differential scanning calorimetry show that the coiled coil conformation is approximately 5 kJ/mol more stable when salt bridges are broken by protonation of the carboxyl side chains. Using guanidinium chloride as the denaturant, the increase in the free energy of unfolding on protonation of the carboxyl side chains is larger, approximately 17 kJ/mol. The discrepancy between urea and guanidinium chloride unfolding can be ascribed to the ionic nature of guanidinium chloride, which screens charge-charge interactions. This work demonstrates the difficulty of predicting the energetic contribution of salt bridges from structural data alone even in a case where the ion pairs are seen in high-resolution NMR structures. The reason is that the contribution to stability results from a fine balance between energetically favorable Coulombic attractions and unfavorable desolvation of charges and conformational constraints of the residues involved in ion pairing. The apparent discrepancy between the results presented here and mutational studies indicating stabilization by salt bridges is discussed and resolved. An explanation is proposed for why interhelical salt bridges are frequently found in natural coiled coils despite evidence that they do not directly contribute to stability.     

Interhelical ion pairing in coiled coils: solution structure of a heterodimeric leucine zipper and determination of pKa values of Glu side chains

Residues of opposite charge often populate heptad positions g (heptad i on chain 1) and e' (heptad i + 1 on chain 2) in dimeric coiled coils and may stabilize the dimer by formation of interchain ion pairs. To investigate the contribution to stability of such electrostatic interactions we have designed a disulfide-linked heterodimeric zipper (AB zipper) consisting of the acidic chain Ac-E-VAQLEKE-VAQAEAE-NYQLEQE-VAQLEHE-CG-NH(2) and the basic chain Ac-E-VQALKKR-VQALKAR-NYAAKQK-VQALRHK-CG-NH(2) in which all e and g positions are occupied by either E or K/R to form a maximum of seven interhelical salt bridges. Temperature-induced denaturation experiments monitored by circular dichroism reveal a stable coiled coil conformation below 50 degrees C and in the pH range 1.2-8.0. Stability is highest at pH approximately 4.0 [DeltaG(U) (37 degrees C) = 5.18 +/- 0.51 kcal mol(-)(1)]. The solution structure of the AB zipper at pH 5.65 has been elucidated on the basis of homonuclear (1)H NMR data collected at 800 MHz [heavy atom rmsd's for the ensemble of 50 calculated structures are 0.47 +/- 0.13 A (backbone) and 0.95 +/- 0.16 A (all)]. Both chains of the AB zipper are almost entirely in alpha-helical conformation and form a superhelix with a left-handed twist. Overhauser connectivities reveal close contacts between g position residues (heptad i on chain 1) and residues d/f (heptad i on chain 1), residues a/d (heptad i + 1 on chain 1), and residue a' (heptad i + 1 on chain 2). Residues in position e (heptad i on chain 1) are in contact with residues a/b/d/f (heptad i on chain 1) and residue d' (heptad i on chain 2). These connectivities hint at a relatively defined alignment of the side chains across the helix interface. Partial H-bond formation between the functional groups of residues g and e'(+1) is observed in the calculated structures. NMR pH titration experiments disclose pK(a) values for Glu delta-carboxylate groups: 4.14 +/- 0.02 (E(1)), 4.82 +/- 0.07 (E(6)), 4.52 +/- 0.01 (E(8)), 4.37 +/- 0.03 (E(13)), 4.11 +/- 0.02 (E(15)), 4.41 +/- 0.07 (E(20)), 4.82 +/- 0.03 (E(22)), 4.65 +/- 0.04 (E(27)), 4.63 +/- 0.03 (E(29)), 4.22 +/- 0.02 (E(1)(')). By comparison with pK(a) of Glu in unfolded peptides ( approximately 4. 3 +/- 0.1), our pK(a) data suggest marginal or even unfavorable contribution of charged Glu to the stability of the AB zipper. The electrostatic energy gained from interhelical ion pairs is likely to be surpassed by hydrophobic energy terms upon protonation of Glu, due to increased hydrophobicity of uncharged Glu and, thus, better packing against apolar residues at the chain interface.     

Optimized variants of the cold shock protein from in vitro selection: structural basis of their high thermostability

The bacterial cold shock proteins (Csp) are widely used as models for the experimental and computational analysis of protein stability. In a previous study, in vitro evolution was employed to identify strongly stabilizing mutations in Bs-CspB from Bacillus subtilis. The best variant found by this approach contained the mutations M1R, E3K and K65I, which raised the midpoint of thermal unfolding of Bs-CspB from 53.8 degrees C to 83.7 degrees C, and increased the Gibbs free energy of stabilization by 20.9 kJ mol(-1). Another selected variant with the two mutations A46K and S48R was stabilized by 11.1 kJ mol(-1). To elucidate the molecular basis of these stabilizations, we determined the crystal structures of these two Bs-CspB variants. The mutated residues are generally well ordered and provide additional stabilizing interactions, such as charge interactions, additional hydrogen bonds and improved side-chain packing. Several mutations improve the electrostatic interactions, either by the removal of unfavorable charges (E3K) or by compensating their destabilizing interactions (A46K, S48R). The stabilizing mutations are clustered at a contiguous surface area of Bs-CspB, which apparently is critically important for the stability of the beta-barrel structure but not well optimized in the wild-type protein.     

Key role of coulombic interactions for the folding transition state of the cold shock protein

The cold shock protein CspB shows a five-stranded beta-sheet structure, and it folds rapidly via a native-like transition state. A previous Phi value analysis showed that most of the residues with Phi values close to one reside in strand beta1, and two of them, Lys5 and Lys7 are partially exposed charged residues. To elucidate how coulombic interactions of these two residues contribute to the energetic organisation of the folding transition state we performed comparative folding experiments in the presence of an ionic denaturant (guanidinium chloride) and a non-ionic denaturant (urea) and a double-mutant analysis. Lys5 contributes 6.6 kJ mol(-1) to the stability of the transition state, and half of it originates from screenable coulombic interactions. Lys7 contributes 5.3 kJ mol(-1), and 3.4 kJ mol(-1) of it are screened by salt. In the folded protein Lys7 interacts with Asp25, and the screenable coulombic interaction between these two residues is fully formed in the transition state. This suggests that long-range coulombic interactions such as those originating from Lys5 and Lys7 of CspB can be important for organizing and stabilizing native-like structure early in protein folding.     

Are acidic and basic groups in buried proteins predicted to be ionized?

Ionizable residues play essential roles in proteins, modulating protein stability, fold and function. Asp, Glu, Arg, and Lys make up about a quarter of the residues in an average protein. Multi-conformation continuum electrostatic (MCCE) calculations were used to predict the ionization states of all acidic and basic residues in 490 proteins. Of all 36,192 ionizable residues, 93.5% were predicted to be ionized. Thirty-five percent have lost 4.08 kcal/mol solvation energy (DeltaDeltaG(rxn)) sufficient to shift a pK(a) by three pH units in the absence of other interactions and 17% have DeltaDeltaG(rxn) sufficient to shift pK(a) by five pH units. Overall 85% of these buried residues (DeltaDeltaG(rxn)>5DeltapK units) are ionized, including 92% of the Arg, 86% of the Asp, 77% of the Glu, and 75% of the Lys. Ion-pair interactions stabilize the ionization of both acids and bases. The backbone dipoles stabilize anions more than cations. The interactions with polar side-chains are also different for acids and bases. Asn and Gln stabilize all charges, Ser and Thr stabilize only acids while Tyr rarely stabilize Lys. Thus, hydroxyls are better hydrogen bond donors than acceptors. Buried ionized residues are more likely to be conserved than those on the surface. There are 3.95 residues buried per 100 residues in an average protein.     

Carboxyl pK(a) values, ion pairs, hydrogen bonding, and the pH-dependence of folding the hyperthermophile proteins Sac7d and Sso7d

Sac7d and Sso7d are homologous, hyperthermophile proteins with a high density of charged surface residues and potential ion pairs. To determine the relative importance of specific amino acid side-chains in defining the stability and function of these Archaeal chromatin proteins, pK(a) values were measured for the acidic residues in both proteins using (13)C NMR chemical shifts. The stability of Sso7d enabled titrations to pH 1 under low-salt conditions. Two aspartate residues in Sso7d (D16 and D35) and a single glutamate residue (G54) showed significantly perturbed pK(a) values in low salt, indicating that the observed pH-dependence of stability was primarily due to these three residues. The pH-dependence of backbone amide NMR resonances demonstrated that perturbation of all three pK(a) values was primarily the result of side-chain to backbone amide hydrogen bonds. Few of the significantly perturbed acidic pK(a) values in Sac7d and Sso7d could be attributed to primarily ion pair or electrostatic interactions. A smaller perturbation of E48 (E47 in Sac7d) was ascribed to an ion pair interaction that may be important in defining the DNA binding surface. The small number (three) of significantly altered pK(a) values was in good agreement with a linkage analysis of the temperature, pH, and salt-dependence of folding. The linkage of the ionization of two or more side-chains to protein folding led to apparent cooperativity in the pH-dependence of folding, although each group titrated independently with a Hill coefficient near unity. These results demonstrate that the acid pH-dependence of protein stability in these hyperthermophile proteins is due to independent titration of acidic residues with pK(a) values perturbed primarily by hydrogen bonding of the side-chain to the backbone. This work demonstrates the need for caution in using structural data alone to argue the importance of ion pairs in stabilizing hyperthermophile proteins.     

Stabilization of the cold shock protein CspB from Bacillus subtilis by evolutionary optimization of Coulombic interactions

The bacterial cold shock proteins (Csp) are used by both experimentalists and theoreticians as model systems for analyzing the Coulombic contributions to protein stability. We employ Proside, a method of directed evolution, to identify stabilized variants of Bs-CspB from Bacillus subtilis. Proside links the increased protease resistance of stabilized protein variants to the infectivity of a filamentous phage. Here, three cspB libraries were used for in vitro selections to explore the stabilizing potential of charged amino acids in Bs-CspB. In the first library codons for nine selected surface residues were partially randomized, in the second one random mutations were introduced non-specifically by error-prone PCR, and in the third one the spontaneous mutation rate of the phage in Escherichia coli was used. Stabilizing mutations were found at the surface positions 1, 3, 46, 48, 65, and 66. The contributions of these mutations to stability were characterized by analyzing them individually and in combination. The best combination (M1R, E3K, K65I, and E66L) increased the midpoint of thermal unfolding of Bs-CspB from 53.8 to 85.0 degrees C. The effects of most mutations are strongly context dependent. A good example is provided by the E3R mutation. It is strongly stabilizing (DeltaDeltaGD=11.1kJ mol(-1)) in the wild-type protein, but destabilizing (DeltaDeltaGD=-4.0kJ mol(-1)) in the A46K/S48R/E66L variant. The stabilizations by charge mutations did not correlate well with the corresponding changes in the protein net charge, and they could not be ascribed to the formation of ion pairs. Previous theoretical analyses did not identify the stabilization caused by the mutations at positions 1, 46, and 48. Also, electrostatics calculations based on protein net charge or charge asymmetry did not predict well the stability changes that occur when charged residues in Bs-CspB are mutated. It remains a challenge to model the Coulombic interactions of charged residues in a protein and to determine their contributions to the Gibbs free energy of protein folding.     

Asp79 makes a large, unfavorable contribution to the stability of RNase Sa

The two most buried carboxyl groups in ribonuclease Sa (RNase Sa) are Asp33 (99% buried; pK 2.4) and Asp79 (85% buried; pK 7.4). Above these pK values, the stability of the D33A variant is 6kcal/mol less than wild-type RNase Sa, and the stability of the D79A variant is 3.3kcal/mol greater than wild-type RNase Sa. The key structural difference between the carboxyl groups is that Asp33 forms three intramolecular hydrogen bonds, and Asp79 forms no intramolecular hydrogen bond. Here, we focus on Asp79 and describe studies of 11 Asp79 variants. Most of the variants were at least 2kcal/mol more stable than wild-type RNase Sa, and the most interesting was D79F. At pH 3, below the pK of Asp79, RNase Sa is 0.3kcal/mol more stable than the D79F variant. At pH 8.5, above the pK of Asp79, RNase Sa is 3.7kcal/mol less stable than the D79F variant. The unfavorable contribution of Asp79 to the stability appears to result from the Born self-energy of burying the charge and, more importantly, from unfavorable charge-charge interactions. To counteract the effect of the negative charge on Asp79, we prepared the Q94K variant and the crystal structure showed that the amino group of the Lys formed a hydrogen-bonded ion pair (distance, 2.71A; angle, 100 degrees ) with the carboxyl group of Asp79. The stability of the Q94K variant was about the same as the wild-type at pH 3, where Asp79 is uncharged, but 1kcal/mol greater than that of wild-type RNase Sa at pH 8.5, where Asp79 is charged. Differences in hydrophobicity, steric strain, Born self-energy, and electrostatic interactions all appear to contribute to the range of stabilities observed in the variants. When it is possible, replacing buried, non-hydrogen bonded, ionizable side-chains with non-polar side-chains is an excellent means of increasing protein stability.     

Stabilizing and destabilizing clusters in the hydrophobic core of long two-stranded alpha-helical coiled-coils

Detailed sequence analyses of the hydrophobic core residues of two long two-stranded alpha-helical coiled-coils that differ dramatically in sequence, function, and length were performed (tropomyosin of 284 residues and the coiled-coil domain of the myosin rod of 1086 residues). Three types of regions were present in the hydrophobic core of both proteins: stabilizing clusters and destabilizing clusters, defined as three or more consecutive core residues of either stabilizing (Leu, Ile, Val, Met, Phe, and Tyr) or destabilizing (Gly, Ala, Cys, Ser, Thr, Asn, Gln, Asp, Glu, His, Arg, Lys, and Trp) residues, and intervening regions that consist of both stabilizing and destabilizing residues in the hydrophobic core but no clusters. Subsequently, we designed a series of two-stranded coiled-coils to determine what defines a destabilizing cluster and varied the length of the destabilizing cluster from 3 to 7 residues to determine the length effect of the destabilizing cluster on protein stability. The results showed a dramatic destabilization, caused by a single Leu to Ala substitution, on formation of a 3-residue destabilizing cluster (DeltaT(m) of 17-21 degrees C) regardless of the stability of the coiled-coil. Any further substitution of Leu to Ala that increased the size of the destabilizing cluster to 5 or 7 hydrophobic core residues in length had little effect on stability (DeltaT(m) of 1.4-2.8 degrees C). These results suggested that the contribution of Leu to protein stability is context-dependent on whether the hydrophobe is in a stabilizing cluster or its proximity to neighboring destabilizing and stabilizing clusters.     

Core side-chain packing and backbone conformation in Lpp-56 coiled-coil mutants

Native proteins exhibit precise geometric packing of atoms in their hydrophobic interiors. Nonetheless, controversy remains about the role of core side-chain packing in specifying and stabilizing the folded structures of proteins. Here we investigate the role of core packing in determining the conformation and stability of the Lpp-56 trimerization domain. The X-ray crystal structures of Lpp-56 mutants with alanine substitutions at two and four interior core positions reveal trimeric coiled coils in which the twist of individual helices and the helix-helix spacing vary significantly to achieve the most favored superhelical packing arrangement. Introduction of each alanine "layer" into the hydrophobic core destabilizes the superhelix by 1.4 kcal mol(-1). Although the methyl groups of the alanine residues pack at their optimum van der Waals contacts in the coiled-coil trimer, they provide a smaller component of hydrophobic interactions than bulky hydrophobic side-chains to the thermodynamic stability. Thus, specific side-chain packing in the hydrophobic core of coiled coils are important determinants of protein main-chain conformation and stability.     

Rearrangement of charge-charge interactions in variant ubiquitins as detected by double-mutant cycles and NMR

Previous studies of ubiquitin disclosed numerous charge-charge interactions on the protein's surface. To investigate how neighboring residues influence the strength of these interactions, double-mutant cycles are combined with pK(a) determinations by 2D NMR. More specifically, the environment around the Asp21-Lys29 ion pair has been altered through mutations at position 25, which is an asparagine in mammalian ubiquitin and a positively-charged residue in many other ubiquitin-like proteins. The pK(a) value of Asp21 decreases by 0.4 to 0.7 pH unit when Asn25 is substituted with a positively charged residue, suggesting a new and favorable ion pair interaction between positions 21 and 25. However, analysis of double mutants reveals that the favorable interaction between Asp21 and Lys29 is weakened when position 25 is a positively charged residue. Interestingly, while the pK(a) value of His25 in the N25H variant agrees with model compound values, additional mutants reveal that this agreement is fortuitous, resulting from a balance of favorable and unfavorable interactions; similar results were observed previously for Glu34 in ubiquitin and His8 in staphylococcal nuclease. Ionizable groups may thus have pK(a) values similar to model compound values and yet still be involved in significant interactions with other protein groups. One surprising result of introducing positively charged residues at position 25 is a new interaction between Lys29 and Glu18, an interaction not present in wild-type ubiquitin. This unanticipated result illustrates a key advantage of using NMR to determine pK(a) values for many residues simultaneously in the variant proteins. Overall, the strength of an interaction between two residues at the surface of ubiquitin is sensitive to the identity of neighboring residues. The results also demonstrate that relatively conservative and common point mutations such as substitutions of polar with charged residues and vice versa can have effects on interactions beyond the site of mutation per se.     

Surface electrostatic interactions contribute little of stability of barnase

Electrostatic interactions are believed to play an important role in stabilizing the native structure of proteins. We have quantified the contribution to stability of an interaction between two oppositely charged side-chains on the surface of barnase. Using site-directed mutagenesis, glutamate 28 and lysine 32 were introduced onto the solvent-accessible side of the second alpha-helix in barnase. These two residues are separated by one turn of the helix, and so are ideally situated for their opposite charges to interact. Double mutant cycle analysis reveals that the interaction between Glu28 and Lys32 contributes only approximately 0.2 kcal/mol to stability of the protein. All other interactions between exposed charged side-chains in barnase examined so far also contribute little to stability. We explain this low value by their location on the surface, rather than in the interior, of the protein.     

Hydrogen bonding markedly reduces the pK of buried carboxyl groups in proteins

The ionizable groups in proteins with the lowest pKs are the carboxyl groups of aspartic acid side-chains. One of the lowest, pK=0.6, is observed for Asp76 in ribonuclease T1. This low pK appeared to result from hydrogen bonds to a water molecule and to the side-chains of Asn9, Tyr11, and Thr91. The results here confirm this by showing that the pK of Asp76 increases to 1.7 in N9A, to 4.0 in Y11F, to 4.2 in T91V, to 4.4 in N9A+Y11F, to 4.9 in N9A+T91V, to 5.9 in Y11F+T91V, and to 6.4 in the triple mutant: N9A+Y11F+T91V. In ribonuclease Sa, the lowest pK=2.4 for Asp33. This pK increases to 3.9 in T56A, which removes the hydrogen bond to Asp33, and to 4.4 in T56V, which removes the hydrogen bond and replaces the -OH group with a -CH(3) group. It is clear that hydrogen bonds are able to markedly lower the pK values of carboxyl groups in proteins. These same hydrogen bonds make large contributions to the conformational stability of the proteins. At pH 7, the stability of D76A ribonuclease T1 is 3.8 kcal mol(-1) less than wild-type, and the stability of D33A ribonuclease Sa is 4.1 kcal mol(-1) less than wild-type. There is a good correlation between the changes in the pK values and the changes in stability. The results suggest that the pK values for these buried carboxyl groups would be greater than 8 in the absence of hydrogen bonds, and that the hydrogen bonds and other interactions of the carboxyl groups contribute over 8 kcal mol(-1) to the stability.     

Origins of the high stability of an in vitro-selected cold-shock protein

In previous work, we had identified stabilized forms of the cold-shock protein Bs-CspB from Bacillus subtilis in a combinatorial library by an in vitro selection procedure. In this library, the sequence positions 2, 3, 46, 64, 66, and 67 had been randomized, because Bs-CspB differs from the naturally thermostable homolog Bc-Csp from Bacillus caldolyticus, among others, at these six positions. For the most stable selected variant, the midpoint of thermal unfolding (tM) increased by 28.2 deg. C and the Gibbs free energy of unfolding (deltaG(D)) by 19 kJ/mol. Here, we analyzed by site-directed mutagenesis how the selected residues contribute individually to this strong stabilization. Val3 and Val66, which replace Glu3 and Glu66 of wild-type Bs-CspB, each contribute about 7 kJ/mol to stability, the Thr64Arg substitution contributes 4.5 kJ/mol, and 3.2 kJ/mol originate from the Ala46Leu replacement. Gly67 at the carboxy terminus is unimportant for stability, the Arg selected at position 2 is overall slightly destabilizing but improves the coulombic interactions. The best variant differs from Bc-Csp at all six positions; nevertheless, natural and in vitro selection followed similar principles. In both cases, negatively charged residues at the adjacent positions 3 and 66 are avoided, and a positively charged residue is introduced into this area of the protein surface. Its exact location is unimportant. It can be at position 3, as in the thermophilic Bc-Csp, or at positions 2 or 64, as in the most stable selected variant. These positively charged residues contribute to stability not by engaging in pairwise coulombic interactions with a specific carboxyl group, but by generally improving the charge distribution in this particular region of the protein surface. These coulombic effects contribute significantly to the thermostability of the cold-shock proteins. They are only weakly interdependent and best explained by the presence of a flexible ion network at the protein surface. Our results emphasize that surface positions are very good candidates for optimizing protein stability.     

Metal cofactors play a dual role in Mycobacterium tuberculosis inorganic pyrophosphatase

Inorganic pyrophosphatase from Mycobacterium tuberculosis (Mt-PPase) is one of the possible targets for the rational design of anti-tuberculosis agents. In this paper, functional properties of this enzyme are characterized in the presence of the most effective activators--Mg2+ and Mn2+. Dissociation constants of Mt-PPase complexed with Mg2+ or Mn2+ are essentially similar to those of Escherichia coli PPase. Stability of a hexameric form of Mt-PPase has been characterized as a function of pH both for the metal-free enzyme and for Mg2+- or Mn2+-enzyme. Hexameric metal-free Mt-PPase has been shown to dissociate, forming monomers at pH below 4 or trimers at pH from 8 to 10. Mg2+ or Mn2+ shift the hexamer-trimer equilibrium found for the apo-Mt-PPase at pH 8-10 toward the hexameric form by stabilizing intertrimeric contacts. The pK(a) values have been determined for groups that control the observed hexamer-monomer (pK(a) 5.4), hexamer-trimer (pK(a) 7.5), and trimer-monomer (pK(a) 9.8) transitions. Our results demonstrate that due to the non-conservative amino acid residues His21 and His86 in the active site of Mt-PPase, substrate specificity of this enzyme, in contrast to other typical PPases, does not depend on the nature of the metal cofactor.     

Crystal structure of the Geobacillus stearothermophilus carboxylesterase Est55 and its activation of prodrug CPT-11

Several mammalian carboxylesterases were shown to activate the prodrug irinotecan (CPT-11) to produce 7-ethyl-10-hydroxycamptothecin (SN-38), a topoisomerase inhibitor used in cancer therapy. However, the potential use of bacterial carboxylesterases, which have the advantage of high stability, has not been explored. We present the crystal structure of the carboxyesterase Est55 from Geobacillus stearothermophilus and evaluation of its enzyme activity on CPT-11. Crystal structures were determined at pH 6.2 and pH 6.8 and resolution of 2.0 A and 1.58 A, respectively. Est55 folds into three domains, a catalytic domain, an alpha/beta domain and a regulatory domain. The structure is in an inactive form; the side-chain of His409, one of the catalytic triad residues, is directed away from the other catalytic residues Ser194 and Glu310. Moreover, the adjacent Cys408 is triply oxidized and lies in the oxyanion hole, which would block the binding of substrate, suggesting a regulatory role. However, Cys408 is not essential for enzyme activity. Mutation of Cys408 showed that hydrophobic side-chains were favorable, while polar serine was unfavorable for enzyme activity. Est55 was shown to hydrolyze CPT-11 into the active form SN-38. The mutant C408V provided a more stable enzyme for activation of CPT-11. Therefore, engineered thermostable Est55 is a candidate for use with irinotecan in enzyme-prodrug cancer therapy.     

Tautomerism,acid-base equilibria,and H-bonding of the six histidines in subtilisin BPN' by NMR

We have determined by (15)N, (1)H, and (13)C NMR, the chemical behavior of the six histidines in subtilisin BPN' and their PMSF and peptide boronic acid complexes in aqueous solution as a function of pH in the range of from 5 to 11, and have assigned every (15)N, (1)H, C(epsilon 1), and C(delta2) resonance of all His side chains in resting enzyme. Four of the six histidine residues (17, 39, 67, and 226) are neutrally charged and do not titrate. One histidine (238), located on the protein surface, titrates with pK(a) = 7.30 +/- 0.03 at 25 degrees C, having rapid proton exchange, but restricted mobility. The active site histidine (64) in mutant N155A titrates with a pK(a) value of 7.9 +/- 0.3 and sluggish proton exchange behavior, as shown by two-site exchange computer lineshape simulation. His 64 in resting enzyme contains an extremely high C(epsilon 1)-H proton chemical shift of 9.30 parts per million (ppm) owing to a conserved C(epsilon 1)-H(.)O=C H-bond from the active site imidazole to a backbone carbonyl group, which is found in all known serine proteases representing all four superfamilies. Only His 226, and His 64 at high pH, exist as the rare N(delta1)-H tautomer, exhibiting (13)C(delta1) chemical shifts approximately 9 ppm higher than those for N(epsilon 2)-H tautomers. His 64 in the PMSF complex, unlike that in the resting enzyme, is highly mobile in its low pH form, as shown by (15)N-(1)H NOE effects, and titrates with rapid proton exchange kinetics linked to a pK(a) value of 7.47 +/- 0.02.     

Contribution of histidine residues to the conformational stability of ribonuclease T1 and mutant Glu-58----Ala

The pK values of the histidine residues in ribonuclease T1 (RNase T1) are unusually high: 7.8 (His-92), 7.9 (His-40), and 7.3 (His-27) [Inagaki et al. (1981) J. Biochem. 89, 1185-1195]. In the RNase T1 mutant Glu-58----Ala, the first two pK values are reduced to 7.4 (His-92) and 7.1 (His-40). These lower pKs were expected since His-92 (5.5 A) and His-40 (3.7 A) are in close proximity to Glu-58 at the active site. The conformational stability of RNase T1 increases by over 4 kcal/mol between pH 9 and 5, and this can be entirely accounted for by the greater affinity for protons by the His residues in the folded protein (average pK = 7.6) than in the unfolded protein (pk approximately 6.6). Thus, almost half of the net conformational stability of RNase T1 results from a difference between the pK values of the histidine residues in the folded and unfolded conformations. In the Glu-58----Ala mutant, the increase in stability between pH 9 and 5 is halved (approximately 2 kcal/mol), as expected on the basis of the lower pK values for the His residues in the folded protein (average pK = 7.1). As a consequence, RNase T1 is more stable than the mutant below pH 7.5, and less stable above pH 7.5. These results emphasize the importance of measuring the conformational stability as a function of pH when comparing proteins differing in structure.