Stabilization of apoflavodoxin by replacing hydrogen-bonded charged Asp or Glu residues by the neutral isosteric Asn or Gln

Knowledge of protein stability principles provides a means to increase protein stability in a rational way. Here we explore the feasibility of stabilizing proteins by replacing solvent-exposed hydrogen-bonded charged Asp or Glu residues by the neutral isosteric Asn or GLN: The rationale behind this is a previous observation that, in some cases, neutral hydrogen bonds may be more stable that charged ones. We identified, in the apoflavodoxin from Anabaena PCC 7119, three surface-exposed aspartate or glutamate residues involved in hydrogen bonding with a single partner and we mutated them to asparagine or glutamine, respectively. The effect of the mutations on apoflavodoxin stability was measured by both urea and temperature denaturation. We observed that the three mutant proteins are more stable than wild-type (on average 0.43 kcal/mol from urea denaturation and 2.8 degrees C from a two-state analysis of fluorescence thermal unfolding data). At high ionic strength, where potential electrostatic repulsions in the acidic apoflavodoxin should be masked, the three mutants are similarly more stable (on average 0.46 kcal/mol). To rule out further that the stabilization observed is due to removal of electrostatic repulsions in apoflavodoxin upon mutation, we analysed three control mutants and showed that, when the charged residue mutated to a neutral one is not hydrogen bonded, there is no general stabilizing effect. Replacing hydrogen-bonded charged Asp or Glu residues by Asn or Gln, respectively, could be a straightforward strategy to increase protein stability.     

Histidine residues at the N- and C-termini of alpha-helices: perturbed pKas and protein stability.

A single histidine residue has been placed at either the N-terminus or the C-terminus of each of the two alpha-helices of barnase. The pKa of that histidine residue in each of the four mutants has been determined by 1H NMR. The pKas of the two residues at the C-terminus are, on average, 0.5 unit higher, and those of the residues at the N-terminus are 0.8 unit lower, than the pKa of histidines in unfolded barnase at low ionic strength. The conformational stability of the mutant proteins at different values of pH has been measured by urea denaturation. C-Terminal histidine mutants are approximately 0.6 kcal mol-1 more stable when the introduced histidine is protonated, both at low and high ionic strength. N-Terminal mutants with a protonated histidine residue are approximately 1.1 kcal mol-1 less stable at low ionic strength and 0.5 kcal mol-1 less stable at high ionic strength (1 M NaCl). The low-field 1H NMR spectra of the mutant proteins at low pH suggest that the C-terminal histidines form hydrogen bonds with the protein while the N-terminal histidines do not form the same. The perturbations of pKa and stability result from a combination of different electrostatic environments and hydrogen-bonding patterns at either ends of helices. The value of 0.6 kcal mol-1 represents a lower limit to the favorable electrostatic interaction between the alpha-helix dipole and a protonated histidine residue at the C-terminal end of the helix.(ABSTRACT TRUNCATED AT 250 WORDS)     

The role of a conserved tyrosine residue in high-potential iron sulfur proteins.

Conserved tyrosine-12 of Ectothiorhodospira halophila high-potential iron sulphur protein (HiPIP) iso-I was substituted with phenylalanine (Y12F), histidine (Y12H), tryptophan (Y12W), isoleucine (Y12I), and alanine (Y12A). Variants Y12A and Y12I were expressed to reasonable levels in cells grown at lower temperatures, but decomposed during purification. Variants Y12F, Y12H, and Y12W were substantially destabilized with respect to the recombinant wild-type HiPIP (rcWT) as determined by differential scanning calorimetry over a pH range of 7.0-11.0. Characterization of the Y12F variant by NMR indicates that the principal structural differences between this variant and the rcWT HiPIP result from the loss of the two hydrogen bonds of the Tyr-12 hydroxyl group with Asn-14 O delta 1 and Lys-59 NH, respectively. The effect of the loss of the latter interaction is propagated through the Lys-59/Val-58 peptide bond, thereby perturbing Gly-46. The delta delta GDapp of Y12F of 2.3 kcal/mol with respect to rcWT HiPIP (25 degrees C, pH 7.0) is entirely consistent with the contribution of these two hydrogen bonds to the stability of the latter. CD measurements show that Tyr-12 influences several electronic transitions within the cluster. The midpoint reduction potentials of variants Y12F, Y12H, and Y12W were 17, 19, and 22 mV (20 mM MOPS, 0.2 M sodium chloride, pH 6.98, 25 degrees C), respectively, higher than that of rcWT HiPIP. The current results indicate that, although conserved Tyr-12 modulates the properties of the cluster, its principle function is to stabilize the HiPIP through hydrogen bonds involving its hydroxyl group and electrostatic interactions involving its aromatic ring.     

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.     

Engineering of the pH optimum of Bacillus cereus beta-amylase: conversion of the pH optimum from a bacterial type to a higher-plant type

The optimum pH of Bacillus cereus beta-amylase (BCB, pH 6.7) differs from that of soybean beta-amylase (SBA, pH 5.4) due to the substitution of a few amino acid residues near the catalytic base residue (Glu 380 in SBA and Glu 367 in BCB). To explore the mechanism for controlling the optimum pH of beta-amylase, five mutants of BCB (Y164E, Y164F, Y164H, Y164Q, and Y164Q/T47M/Y164E/T328N) were constructed and characterized with respect to enzymatic properties and X-ray structural crystal analysis. The optimum pH of the four single mutants shifted to 4.2-4.8, approximately 2 pH units and approximately 1 pH unit lower than those of BCB and SBA, respectively, and their k(cat) values decreased to 41-3% of that of the wild-type enzyme. The X-ray crystal analysis of the enzyme-maltose complexes showed that Glu 367 of the wild type is surrounded by two water molecules (W1 and W2) that are not found in SBA. W1 is hydrogen-bonded to both side chains of Glu 367 and Tyr 164. The mutation of Tyr 164 to Glu and Phe resulted in the disruption of the hydrogen bond between Tyr 164 Oeta and W1 and the introduction of two additional water molecules near position 164. In contrast, the triple mutant of BCB with a slightly decreased pH optimum at pH 6.0 has no water molecules (W1 and W2) around Glu 367. These results suggested that a water-mediated hydrogen bond network (Glu 367...W1...Tyr 164...Thr 328) is the primary requisite for the increased pH optimum of wild-type BCB. This strategy is completely different from that of SBA, in which a hydrogen bond network (Glu 380...Thr 340...Glu 178) reduces the optimum pH in a hydrophobic environment.     

Crystal structure of a pH-stabilized mutant of villin headpiece

Villin-type headpiece domains are compact F-actin-binding motifs that have been used extensively as a model system to investigate protein folding by both experimental and computational methods. Villin headpiece (HP67) harbors a highly helical, thermostable, and autonomously folding subdomain in the C terminus (HP35), and because of this feature, HP67 is usually considered to be composed of a N- and C-terminal subdomain. Unlike the C-terminal subdomain, the N-terminal subdomain consists mainly of loops and turns, and the folding is dependent upon the presence of the C-terminal subdomain. The pH sensitivity of this subdomain is thought to arise from, at least partially, protonation of H41 buried in the hydrophobic core. Substitution of this histidine with tyrosine, another permissive residue at this position for naturally occurring sequences, increases not only the pH stability of HP67 but also the thermal stability and the cooperativity of thermal unfolding over a wide pH range (0.9-7.5). The crystal structures of wild-type HP67 and the H41Y mutant, determined under the same conditions, indicate that the H41Y substitution causes only localized rearrangement around the mutated residue. The F-actin-binding motif remains essentially the same after the mutation, accounting for the negligible effect of the mutation on F-actin affinity. The hydrogen bond formed between the imidazole ring of H41 and the backbone carbonyl of E14 of HP67 is eliminated by the H41Y mutation, which renders the extreme N terminus of H41Y more mobile; the hydrogen bond formed between the imidazole ring of H41 and the backbone nitrogen of D34 is replaced with that between the hydroxyl group of Y41 and the backbone nitrogen of D34 after the H41Y substitution. The increased hydrophobicity of tyrosine compensates for the loss of hydrogen bonds in the extreme N terminus and accounts for the increased stability and cooperativity of the H41Y mutant.     

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.     

Role of hydrogen bond networks and dynamics in positive and negative cooperative stabilization of a protein

Cooperativity mediated through hydrogen bond networks in yeast iso-1-cytochrome c was studied using a thermodynamic triple mutant cycle. Three known stabilizing mutations, Asn 26 to His, Asn 52 to Ile, and Tyr 67 to Phe, were used to construct the triple mutant cycle. The side chain of His 26, a wild-type residue, forms two hydrogen bonds that bridge two substructures of the wild-type protein, and Tyr 67 and Asn 52 are part of an extensive buried hydrogen bond network. The stabilities of all variants in the triple mutant cycle were determined by guanidine hydrochloride denaturation methods and used to determine the pairwise, Delta(2)G(int), and triple interaction energies. His 26 and Ile 52 interact cooperatively (Delta(2)G(int) is 1-2 kcal/mol), whereas the two other pairs of mutations interact anticooperatively (Delta(2)G(int) is -0.5 to -1.5 kcal/mol). Previously reported structural data for iso-1-cytochrome c variants containing these mutations show that changes in the strength of the His 26 to Glu 44 hydrogen bond, apparently caused by changes in main chain dynamics, provide a mechanism for the long distance (His 26 to Phe 67 and His 26 to Ile 52) propagation of pairwise interaction energies. Opposing changes in the strength of the His 26 to Glu 44 hydrogen bond caused by the N52I and Y67F mutations generate a negative triple interaction energy (-0.9 +/-0.7 kcal/mol) that combined with cancellation of cooperative and anticooperative pairwise interactions produce apparent additivity for the stabilizing effects of the single mutations in the triple mutant variant.     

Removing a hydrogen bond in the dimer interface of Escherichia coli manganese superoxide dismutase alters structure and reactivity

Among manganese superoxide dismutases, residues His30 and Tyr174 are highly conserved, forming part of the substrate access funnel in the active site. These two residues are structurally linked by a strong hydrogen bond between His30 NE2 from one subunit and Tyr174 OH from the other subunit of the dimer, forming an important element that bridges the dimer interface. Mutation of either His30 or Tyr174 in Escherichia coli MnSOD reduces the superoxide dismutase activity to 30--40% of that of the wt enzyme, which is surprising, since Y174 is quite remote from the active site metal center. The 2.2 A resolution X-ray structure of H30A-MnSOD shows that removing the Tyr174-->His30 hydrogen bond from the acceptor side results in a significant displacement of the main-chain segment containing the Y174 residue, with local rearrangement of the protein. The 1.35 A resolution structure of Y174F-MnSOD shows that disruption of the same hydrogen bond from the donor side has much greater consequences, with reorientation of F174 having a domino effect on the neighboring residues, resulting in a major rearrangement of the dimer interface and flipping of the His30 ring. Spectroscopic studies on H30A, H30N, and Y174F mutants show that (like the previously characterized Y34F mutant of E. coli MnSOD) all lack the high pH transition of the wt enzyme. This observation supports assignment of the pH sensitivity of MnSOD to coordination of hydroxide ion at high pH rather than to ionization of the phenolic group of Y34. Thus, mutations near the active site, as in the Y34F mutant, as well as at remote positions, as in Y174F, similarly affect the metal reactivity and alter the effective pK(a) for hydroxide ion binding. These results imply that hydrogen bonding of the H30 imidazole N--H group plays a key role in substrate binding and catalysis.     

Probing electrostatic interactions along the reaction pathway of a glycoside hydrolase: histidine characterization by NMR spectroscopy

We have characterized by NMR spectroscopy the three active site (His80, His85, and His205) and two non-active site (His107 and His114) histidines in the 34 kDa catalytic domain of Cellulomonas fimi xylanase Cex in its apo, noncovalently aza-sugar-inhibited, and trapped glycosyl-enzyme intermediate states. Due to protection from hydrogen exchange, the level of which increased upon inhibition, the labile 1Hdelta1 and 1H epsilon1 atoms of four histidines (t1/2 approximately 0.1-300 s at 30 degrees C and pH approximately 7), as well as the nitrogen-bonded protons in the xylobio-imidazole and -isofagomine inhibitors, could be observed with chemical shifts between 10.2 and 17.6 ppm. The histidine pKa values and neutral tautomeric forms were determined from their pH-dependent 13C epsilon1-1H epsilon1 chemical shifts, combined with multiple-bond 1H delta2/epsilon1-15N delta1/epsilon2 scalar coupling patterns. Remarkably, these pKa values span more than 8 log units such that at the pH optimum of approximately 6 for Cex activity, His107 and His205 are positively charged (pKa > 10.4), His85 is neutral (pKa < 2.8), and both His80 (pKa = 7.9) and His114 (pKa = 8.1) are titrating between charged and neutral states. Furthermore, upon formation of the glycosyl-enzyme intermediate, the pKa value of His80 drops from 7.9 to <2.8, becoming neutral and accepting a hydrogen bond from an exocyclic oxygen of the bound sugar moiety. Changes in the pH-dependent activity of Cex due to mutation of His80 to an alanine confirm the importance of this interaction. The diverse ionization behaviors of the histidine residues are discussed in terms of their structural and functional roles in this model glycoside hydrolase.     

Increasing Protein Stability by Polar Surface Residues: Domain-Wide Consequences of Interactions Within a Loop

We have examined the influence of surface hydrogen bonds on the stability of proteins by studying the effects of mutations of human immunoglobulin light chain variable domain (V_L). In addition to the variants Y27dD, N28F, and T94H of protein κIV Len that were previously described, we characterized mutants M4L, L27cN, L27cQ, and K39T, double mutant M4L/Y27dD, and triple mutant M4L/Y27dD/T94H. The triple mutant had an enhanced thermodynamic stability of 4.2 kcal/mol. We determined the structure of the triple mutant by x-ray diffraction and correlated the changes in stability due to the mutations with changes in the three-dimensional structure. Y27dD mutant had increased stability of Len by 2.7 kcal/mol, a large value for a single mutation. Asp27d present in CDR1 formed hydrogen bonds with the side-chain and main-chain atoms within the loop. In the case of the K39T mutant, which reduces stability by 2 kcal/mol, Lys39 in addition to forming a hydrogen bond with a carbonyl oxygen of a neighboring loop may also favorably influence the surface electrostatics of the molecule. We showed that hydrogen bonds between residues in surface loops can add to the overall stability of the V_L domains. The contribution to stability is further increased if the surface residue makes more than one hydrogen bond or if it forms a hydrogen bond between neighboring turns or loops separated from each other in the amino acid sequence. Based on our experiments we suggest that stabilization of proteins might be systematically accomplished by introducing additional hydrogen bonds on the surface. These substitutions are more straightforward to predict than core-packing interactions and can be selected to avoid affecting the protein’s function.     

Solution conformation of the His-47 to Ala-47 mutant of Pseudomonas stutzeri ZoBell ferrocytochrome c-551

In the cytochrome c-551 family, the heme 17-propionate caboxylate group is always hydrogen bonded to an invariant Trp-56 and conserved residues (His and Arg mainly, Lys occasionally) at position 47. The mutation of His-47 to Ala-47 for Pseudomas stutzeri ZoBell cytochrome c-551 removes this otherwise invariant hydrogen bond. The solution structure of ferrous H47A has been solved based on NMR-derived constraints. Results indicate that the mutant has very similar main chain folding compared to wild-type. However, less efficient packing of residues in the mutant surrounding the heme propionates leads to more solvent exposure for both propionate groups, which may account for decreased stability of the mutant. The mutant has a reduction potential different from wild-type, and furthermore, the pH dependence of this potential is not the same as for wild-type. The structure of the mutant suggests that these changes are related to the loss of the residue-47 propionate hydrogen bond and the loss of charge on the side chain of residue 47.     

Trans-substitution of the proximal hydrogen bond in myoglobin: II. Energetics, functional consequences, and implications for hemoglobin allostery

The trans-substituted histidine to glycine mutant of sperm whale myoglobin (H93G Mb) is used to study energetics of proximal hydrogen bonding, proximal ligand-heme interactions, and coupling to distal ligand binding. Comparison of mono- and dimethylimidazole structural isomers shows that the hydrogen bond between the proximal ligand and the neighboring Ser92 hydroxyl (position F7) is stabilizing. The range of hydrogen bond stabilities measured here for different distal ligand complexes ranges from -0.7 kcal/mol (monomethylimidazole isomers to MbCO) to -4.1 kcal/mol (dimethylimidazole isomers to MbCN). This range of hydrogen bond stabilities, which is similar to that seen in protein mutagenesis unfolding studies, demonstrates the high sensitivity of the hydrogen bond to modest structural perturbations. The degree to which the 2-methyl group destabilizes proximal ligand binding is found to depend inversely on the total electronic spin. For monomethylimidazole proximal ligands, distal ligand binding weakens the proximal hydrogen bond compared to deoxyMb. Surprisingly, this trend is largely reversed for the dimethylimidazole proximal ligands. These results demonstrate strong coupling between the proximal protein matrix and distal ligand binding. These results provide an explanation for the strong avoidance of hydrogen bonding residues at position F7 in hemoglobin sequences.     

Removal of distal protein-water hydrogen bonds in a plant epoxide hydrolase increases catalytic turnover but decreases thermostability

A putative proton wire in potato soluble epoxide hydrolase 1, StEH1, was identified and investigated by means of site-directed mutagenesis, steady-state kinetic measurements, temperature inactivation studies, and X-ray crystallography. The chain of hydrogen bonds includes five water molecules coordinated through backbone carbonyl oxygens of Pro(186), Leu(266), His(269), and the His(153) imidazole. The hydroxyl of Tyr(149) is also an integrated component of the chain, which leads to the hydroxyl of Tyr(154). Available data suggest that Tyr(154) functions as a final proton donor to the anionic alkylenzyme intermediate formed during catalysis. To investigate the role of the putative proton wire, mutants Y149F, H153F, and Y149F/H153F were constructed and purified. The structure of the Y149F mutant was solved by molecular replacement and refined to 2.0 A resolution. Comparison with the structure of wild-type StEH1 revealed only subtle structural differences. The hydroxyl group lost as a result of the mutation was replaced by a water molecule, thus maintaining a functioning hydrogen bond network in the proton wire. All mutants showed decreased catalytic efficiencies with the R,R-enantiomer of trans-stilbene oxide, whereas with the S,S-enantiomer, k (cat)/K (M) was similar or slightly increased compared with the wild-type reactions. k (cat) for the Y149F mutant with either TSO enantiomer was increased; thus the lowered enzyme efficiencies were due to increases in K (M). Thermal inactivation studies revealed that the mutated enzymes were more sensitive to elevated temperatures than the wild-type enzyme. Hence, structural alterations affecting the hydrogen bond chain caused increases in k (cat) but lowered thermostability.     

A general approach for detecting folding intermediates from steady-state and time-resolved fluorescence of single-tryptophan-containing proteins

During denaturant-induced equilibrium (un)folding of wild-type apoflavodoxin from Azotobacter vinelandii, a molten globule-like folding intermediate is formed. This wild-type protein contains three tryptophans. In this study, we use a general approach to analyze time-resolved fluorescence and steady-state fluorescence data that are obtained upon denaturant-induced unfolding of a single-tryptophan-containing variant of apoflavodoxin [i.e., W74/F128/F167 (WFF) apoflavodoxin]. The experimental data are assembled in matrices, and subsequent singular-value decomposition of these matrices (i.e., based on either steady-state or time-resolved fluorescence data) shows the presence of three significant, and independent, components. Consequently, to further analyze the denaturation trajectories, we use a three-state protein folding model in which a folding intermediate and native and unfolded protein molecules take part. Using a global analysis procedure, we determine the relative concentrations of the species involved and show that the stability of WFF apoflavodoxin against global unfolding is ～4.1 kcal/mol. Analysis of time-resolved anisotropy data of WFF apoflavodoxin unfolding reveals the remarkable observation that W74 is equally well fixed within both the native protein and the molten globule-like folding intermediate. Slight differences between the direct environments of W74 in the folding intermediate and native protein cause different rotameric populations of the indole in both folding species as fluorescence lifetime analysis reveals. Importantly, thermodynamic analyses of the spectral denaturation trajectories of the double-tryptophan-containing protein variants WWF apoflavodoxin and WFW apoflavodoxin show that these variants are significantly more stable (5.9 kcal/mol and 6.8 kcal/mol, respectively) than WFF apoflavodoxin (4.1 kcal/mol) Hence, tryptophan residues contribute considerably to the 10.5 kcal/mol thermodynamic stability of native wild-type apoflavodoxin.     

The role of Tyr41 and His155 in the functional properties of superoxide dismutase from the archaeon Sulfolobus solfataricus

We have examined and compared the effects of mutating Y41 and H155 in the iron superoxide dismutase (SOD) from the archaeon Sulfolobus solfataricus (Ss). These two neighboring residues in the active site are known to have crucial functions in structurally related SODs from different sources. The metal analysis indicates a slightly lower iron content after either Y41F or H155Q replacement, without any significant substitution of iron for manganese. The specific activity of SsSOD referred to the iron content is 17-fold reduced in the Y41F mutant, whereas it is less than 2-fold reduced by the H155Q mutation. The noticeable pH dependence of the activity of SsSOD and H155Q-SsSOD, due to the ionization of Y41 (pK 8.4), is lost in Y41F-SsSOD. After H155Q and even more after the Y41F substitution, the archaeal enzyme acquires a moderate sensitivity to sodium azide inhibition. The hydrogen peroxide inactivation of SsSOD is significantly increased after H155Q replacement; however, both mutants are insensitive to the modification of residue 41 by phenylmethanesulfonyl fluoride. Heat inactivation studies showed that the high stability of SsSOD is reduced by the H155Q mutation; however, upon the addition of SDS, a much faster inactivation kinetics is observed both with wild-type and mutant SsSOD forms. The detergent is also required to follow thermal denaturation of the archaeal enzyme by Fourier transform infrared spectroscopy; these studies gave information about the effect of mutations and modification on flexibility and compactness of the protein structure. The crystal structure of Y41F mutant revealed an uninterrupted hydrogen bond network including three solvent molecules connecting the iron-ligating hydroxide ion via H155 with F41 and H37, which is not present in structures of the corresponding mutant SODs from other sources. These data suggest that Y41 and H155 are important for the structural and functional properties of SsSOD; in particular, Y41 seems to be a powerful regulator of the activity of SsSOD, whereas H155 is apparently involved in the organization of the active site of the enzyme.     

Mutations of surface residues in Anabaena vegetative and heterocyst ferredoxin that affect thermodynamic stability as determined by guanidine hydrochloride denaturation.

The stability properties of oxidized wild-type (wt) and site-directed mutants in surface residues of vegetative (Vfd) and heterocyst (Hfd) ferredoxins from Anabaena 7120 have been characterized by guanidine hydrochloride (Gdn-HCl) denaturation. For Vfd it was found that mutants E95K, E94Q, F65Y, F65W, and T48A are quite similar to wt in stability. E94K is somewhat less stable, whereas E94D, F65A, F65I, R42A, and R42H are substantially less stable than wt. R42H is a substitution found in all Hfds, and NMR comparison of the Anabaena 7120 Vfd and Hfd showed the latter to be much less stable on the basis of hydrogen exchange rates (Chae YK, Abildgaard F, Mooberry ES, Markley JL, 1994, Biochemistry 33:3287-3295); we also find this to be true with respect to Gdn-HCl denaturation. Strikingly, the Hfd mutant H42R is more stable than the wt Hfd by precisely the amount of stability lost in Vfd upon mutating R42 to H (2.0 kcal/mol). On the basis of comparison of the X-ray crystal structures of wt Anabaena Vfd and Hfd, the decreased stabilities of F65A and F65I can be ascribed to increased solvent exposure of interior hydrophobic groups. In the case of Vfd mutants E94K and E94D, the decreased stabilities may result from disruption of a hydrogen bond between the E94 and S47 side chains. The instability of the R42 mutants is also most probably due to decreased hydrogen bonding capabilities.(ABSTRACT TRUNCATED AT 250 WORDS)     

Contribution of a single hydrogen bond between betaHis81 of MHC class II I-E(k) and the bound peptide to the pH-dependent thermal stability

To determine the energetic contribution of the hydrogen bond between betaHis81 of the major histocompatibility complex class II (MHC II) molecule, I-E(k), and the bound hemoglobin peptide (Hb), we analyzed the thermal stability of the hydrogen bond-disrupted mutant, I-E(k)-Hb betaH81Y, in which the betaHis81 residue was replaced with Tyr, by differential scanning calorimetry. The thermal stability of the I-E(k)-Hb betaH81Y mutant was lower than that of the I-E(k)-Hb wild-type, mainly due to the decreased enthalpy change. The difference in the denaturation temperature of the I-E(k)-Hb betaH81Y mutant as compared with that of the I-E(k)-Hb wild-type at pH 5.5 was only slightly smaller than that at pH 7.4, in agreement with the increased stability at an acidic pH, a unique characteristic of MHC II. Thus, the hydrogen bond contributed by betaHis81 is critical for peptide binding, and is independent of pH, which can alter the hydrophilicity of the His residue.     

Increasing protein stability by altering long-range coulombic interactions.

It is difficult to increase protein stability by adding hydrogen bonds or burying nonpolar surface. The results described here show that reversing the charge on a side chain on the surface of a protein is a useful way of increasing stability. Ribonuclease T1 is an acidic protein with a pI approximately 3.5 and a net charge of approximately -6 at pH 7. The side chain of Asp49 is hyperexposed, not hydrogen bonded, and 8 A from the nearest charged group. The stability of Asp49Ala is 0.5 kcal/mol greater than wild-type at pH 7 and 0.4 kcal/mol less at pH 2.5. The stability of Asp49His is 1.1 kcal/mol greater than wild-type at pH 6, where the histidine 49 side chain (pKa = 7.2) is positively charged. Similar results were obtained with ribonuclease Sa where Asp25Lys is 0.9 kcal/mol and Glu74Lys is 1.1 kcal/mol more stable than the wild-type enzyme. These results suggest that protein stability can be increased by improving the coulombic interactions among charged groups on the protein surface. In addition, the stability of RNase T1 decreases as more hydrophobic aromatic residues are substituted for Ala49, indicating a reverse hydrophobic effect.     

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.