On the pH dependence of protein stability

This paper treats the free energy contribution of ionizable groups to protein stability. A method is presented for the calculation of the pH dependence of the denaturation free energy of a protein, which yields results that can be compared directly to experiment. The first step in the treatment is the determination of the average charges of all the ionizable groups in both the folded and unfolded protein. An expression due to Tanford then relates the pH dependence of the unfolding free energy to the difference in net charge between the two states. In order to determine absolute rather than relative unfolding free energies, it is necessary to calculate the total contribution of ionizable groups to protein stability at some reference pH. This is accomplished through a statistical mechanical treatment similar to the one used previously in the calculation of pKas. The treatment itself is rigorous but it suffers from uncertainties in the pKa calculations. Nevertheless, the overall shape of experimentally observed plots of denaturation free energy as a function of pH are reasonably well reproduced by the calculations. A number of general conclusions that arise from the analysis are: (1) knowledge of titration curves and/or effective pKa values of ionizable groups in proteins is sufficient to calculate the pH dependence of the denaturation free energy with respect to some reference pH value. However, in order to calculate the absolute contribution of ionizable groups to protein stability, it is necessary to also know the intrinsic pKa of each group. This is defined as the pKa of a group in a hypothetical state of the protein where all other groups are neutral. (2) Due to desolvation effects, ionizable groups destabilize proteins, although the effect is strongly dependent on pH. There are however, strongly stabilizing pairwise Coulombic interactions on the surface of proteins. (3) Plots of stability versus pH should not be interpreted in terms of a group whose pKa corresponds to the titration midpoint, but rather to a group with different pKas (that correspond approximately to the titration end points) in each state. (4) Any residual structure in the GuHCl-denatured state of proteins appears to have little effect on the pH dependence of stability. (5) pH-dependent unfolding, for example to the "molten globule" state, appears due to individual groups with anomalous pKas whose locations on the protein surface may determine the nature of the unfolded state.     

Effect of pH on formation of a nativelike intermediate on the unfolding pathway of a Lys 73 --> His variant of yeast iso-1-cytochrome c

Previous work on a Lys 73 --> His (H73) variant of iso-1-cytochrome c at pH 7.5 [Godbole et al. (1997) Biochemistry 36, 119-126] showed that this variant unfolds through a nativelike intermediate that has properties consistent with replacement of the Met 80 heme ligand by His 73. Here, the pH dependence of the equilibrium unfolding of the wild type (WT) and H73 proteins have been investigated, since a characteristic pH dependence is expected for the stability of an intermediate stabilized by histidine-heme ligation. Stability has been evaluated using guanidine hydrochloride and pH denaturation methods. Above pH 5, the m-values from guanidine hydrochloride denaturation of the WT and H73 variants remain significantly different, consistent with continued population of this intermediate. At pH 4.5 the m-values for the two proteins are within error the same. To assess stability at lower pH, acid denaturation was carried out. The midpoint is about 3.3 for both proteins but the transition is broader for the H73 protein, suggestive of intermediates again being populated during the unfolding of the H73 protein at this lower pH. Heme ligation by Met 80 was monitored (695 nm absorbance) during gdnHCl (pH 4.5 and 5.0) and acid denaturation, confirming, respectively, the absence and presence of intermediates. A thermodynamic analysis demonstrates that this complex pH dependence for the presence of histidine ligation induced intermediates is expected and implicates a titratable group with a pKa of approximately 6.6. The analysis also demonstrates when the pH dependences of global stability and stability of an intermediate differ significantly, population of folding intermediates as a function of pH will show novel behavior.     

Chemical denaturation: potential impact of undetected intermediates in the free energy of unfolding and m-values obtained from a two-state assumption

The chemical unfolding transition of a protein was simulated, including the presence of an intermediate (I) in equilibrium with the native (N) and unfolded (U) states. The calculations included free energies of unfolding, DeltaGuw, in the range of 1.4 kcal/mol to 10 kcal/mol and three different global m-values. The simulations included a broad range of equilibrium constants for the N left arrow over right arrow I process. The dependence of the N <--> I equilibrium on the concentration of denaturant was also included in the simulations. Apparent DeltaGuw and m-values were obtained from the simulated unfolding transitions by fitting the data to a two-state unfolding process. The potential errors were calculated for two typical experimental situations: 1) the unfolding is monitored by a physical property that does not distinguish between native and intermediate states (case I), and 2) the physical property does not distinguish between intermediate and unfolded states (case II). The results obtained indicated that in the presence of an intermediate, and in both experimental situations, the free energy of unfolding and the m-values could be largely underestimated. The errors in DeltaGuw and m-values do not depend on the m-values that characterize the global N <--> U transition. They are dependent on the equilibrium constant for the N <--> I transition and its characteristic m1-value. The extent of the underestimation increases for higher energies of unfolding. Including no random error in the simulations, it was estimated that the underestimation in DeltaGuw could range between 25% and 35% for unfolding transitions of 3-10 kcal/mol (case I). In case II, the underestimation in DeltaGuw could be even larger than in case I. In the same energy range, a 50% error in the m-value could also take place. The fact that most of the mutant proteins are characterized by both a lower m-value and a lower stability than the wild-type protein suggests that in some cases the results could have been underestimated due to the application of the two-state assumption.     

Redox properties of human medium-chain acyl-CoA dehydrogenase, modulation by charged active-site amino acid residues

The modulation of the electron-transfer properties of human medium-chain acyl-CoA dehydrogenase (hwtMCADH) has been studied using wild-type and site-directed mutants by determining their midpoint potentials at various pH values and estimating the involved pKs. The mutants used were E376D, in which the negative charge is retained; E376Q, in which one negative charge (pKa approximately 6. 0) is removed from the active center; E99G, in which a different negative charge (pKa approximately 7.3) also is affected; and E376H (pKa approximately 9.3) in which a positive charge is present. Em for hwtMCADH at pH 7.6 is -0.114 V. Results for the site-directed mutants indicate that loss of a negative charge in the active site causes a +0.033 V potential shift. This is consistent with the assumption that electrostatic interactions (as in the case of flavodoxins) and specific charges are important in the modulation of the electron-transfer properties of this class of dehydrogenases. Specifically, these charge interactions appear to correlate with the positive Em shift observed upon binding of substrate/product couple to MCADH [Lenn, N. D., Stankovich, M. T., and Liu, H. (1990) Biochemistry 29, 3709-3715], which coincides with a pK increase of Glu376-COOH from approximately 6 to 8-9 [Rudik, I., Ghisla, S., and Thorpe, C. (1998) Biochemistry 37, 8437-8445]. From the pH dependence of the midpoint potentials of hwtMCADH two mechanistically important ionizations are estimated. The pKa value of approximately 6.0 is assigned to the catalytic base, Glu376-COOH, in the oxidized enzyme based on comparison with the pH behavior of the E376H mutant, it thus coincides with the pK value recently estimated [Vock, P., Engst, S., Eder, M., and Ghisla, S. (1998) Biochemistry 37, 1848-1860]. The pKa of approximately 7.1 is assigned to Glu376-COOH in reduced hwtMCADH. Comparable values for these pKas for Glu376-COOH in pig kidney MCADH are pKox = 6.5 and pKred = 7.9. The Em measured for K304E-MCADH (a major mutant resulting in a deficiency syndrome) is essentially identical to that of hwtMCADH, indicating that the disordered enzyme has an intact active site.     

Effects of pH on rhodopsin photointermediates from lumirhodopsin to metarhodopsin II

Time-resolved absorption difference spectra of membrane suspensions of bovine rhodopsin at pH 5, 6, 7, 8, 9, and 10 were collected in the time range from 1 micro s to 200 ms after laser photolysis with 7-ns pulses of 477-nm light. The data were analyzed using singular value decomposition (SVD) and global exponential fitting. At pH 7 the data agree well with previously obtained data (Thorgeirsson et al. (1993) Biochemistry 32, 13861-13872) with fits improved at all pH's by inclusion of a small component due to an absorbance change caused by rotational diffusion which is detectable even at magic angle polarization. A "square scheme" suggested to best explain the previous data, which involves two branches following decay of the lumi intermediate with pathways (1) lumi --> MI480 right harpoon over left harpoon MII and (2) lumi right harpoon over left harpoon MI380 --> MII, could be confirmed throughout the entire pH range. However, to account for the increased rate of the MII --> MI480 reaction in path 1 for rising pH values, we propose that the MII in the square scheme consists of deprotonated MII and protonated MIIH+ forms in rapid equilibrium with each other, resulting in an extended square scheme and increasing the number of 380-nm products from two to three. In addition to the kinetic processes described by the extended square scheme, above pH 8 fast ( approximately 10 micro s) and slow ( approximately 50 ms) components were found. The fast component was assigned to the decay of a blue-shifted lumi intermediate, and the slow component, resolvable only at pH 10, was assigned to formation of a 450 nm absorbing photoproduct.     

Clathrin self-assembly is regulated by three light-chain residues controlling the formation of critical salt bridges

Clathrin self-assembly into a polyhedral lattice mediates membrane protein sorting during endocytosis and organelle biogenesis. Lattice formation occurs spontaneously in vitro at low pH and, intracellularly, is triggered by adaptors at physiological pH. To begin to understand the cellular regulation of clathrin polymerization, we analyzed molecular interactions during the spontaneous assembly of recombinant hub fragments of the clathrin heavy chain, which bind clathrin light-chain subunits and mimic the self-assembly of intact clathrin. Reconstitution of hubs using deletion and substitution mutants of the light-chain subunits revealed that the pH dependence of clathrin self-assembly is controlled by only three acidic residues in the clathrin light-chain subunits. Salt inhibition of hub assembly identified two classes of salt bridges which are involved and deletion analysis mapped the clathrin heavy-chain regions participating in their formation. These combined observations indicated that the negatively charged regulatory residues, identified in the light-chain subunits, inhibit the formation of high-affinity salt bridges which would otherwise induce clathrin heavy chains to assemble at physiological pH. In the presence of light chains, clathrin self-assembly depends on salt bridges that form only at low pH, but is exquisitely sensitive to regulation. We propose that cellular clathrin assembly is controlled via the simple biochemical mechanism of reversing the inhibitory effect of the light-chain regulatory sequence, thereby promoting high-affinity salt bridge formation.     

Determination of the pKa values of titratable groups of an antigen-antibody complex, HyHEL-5-hen egg lysozyme

The titration behavior of the ionizable residues of the HyHEL-5-hen egg lysozyme complex and its individual components has been studied using continuum electrostatic calculations. Several residues of HyHEL-5 had pKa values shifted away from model values for isolated residues by more than three pH units. Shifts away from the model values were smaller for the residues of hen egg lysozyme. A moderate variation in the pKa values of the titratable groups was observed upon increase of the ionic strength from 0 to 100 mM, amounting to 1-2 pH units in most cases. Under physiological conditions, the net charge of HyHEL-5 was opposite that for hen egg lysozyme. Several residues, including those involved in the Arg-Glu salt bridges that have been proposed to be important in antibody-antigen binding, had pKa values that were changed significantly upon binding. The main titration event upon antibody-antigen binding appears to be loss of a proton from residue GluH50 of the Fv molecule. The limitations of our calculation methods and the role they might play in the design of antibodies for use in assays, sensors and separations are discussed.     

Kinetic role of electrostatic interactions in the unfolding of hyperthermophilic and mesophilic rubredoxins

The temperature dependence of the unfolding kinetics of rubredoxins from the hyperthermophile Pyrococcus furiosus (RdPf) and the mesophile Clostridium pasteurianum (RdCp) has been studied. Results show that RdPf unfolds much more slowly, under all experimentally accessible temperature regimes, than RdCp and other typical mesophilic proteins. Rates of RdCp and RdPf unfolding decrease upon increasing the pH above 2 and diverge dramatically at pH 7. As shown by detailed electrostatic energy calculations, this is the result of a differential degree of protonation of the negatively charged amino acids, which causes distinct electrostatic configurations as a function of pH. We propose that ion pairs, particularly those that are placed in key surface positions, may play a kinetic role by mildly clamping the protein and thereby influencing the nature and the number of the vibrational normal modes that are thermally accessible upon unfolding. More generally, these modes are also likely to be affected by the favorable electrostatic configurations, which we have shown to be directly linked to the extremely slow unfolding rates of RdPf at neutral pH. Even at pH 2, in the absence of any salt bridges, the unfolding rates of RdPf are much smaller than those of RdCp. This is ascribed to presently unidentified structural elements of nonelectrostatic nature. Since electrostatic effects influence the unfolding kinetics of both mesophilic and thermophilic rubredoxins, these findings may be of general significance for proteins.     

Sulfhydryl oxidation of mutants with cysteine in place of acidic residues in the lactose permease

To examine further the role of charge-pair interactions in the structure and function of lactose permease, Asp237 (helix VII), Asp240 (helix VII), Glu126 (cytoplasmic loop IV/V), Glu269 (helix VIII), and Glu325 (helix X) were replaced individually with Cys in a functional mutant devoid of Cys residues. Each mutant was then oxidized with H2O2 in order to generate a sulfinic and/or sulfonic acid at these positions. Due to the isosteric relationship between aspartate and sulfinate, in particular, and the lower pKa of the sulfinic and sulfonic acid side chains, oxidized derivatives of Cys are useful probes for examining the role of carboxylates. Asp237-->Cys or Asp240-->Cys permease is inactive, as shown previously, but H2O2 oxidation restores activity to an extent similar to that observed when a negative charge is reintroduced by other means. Glu126-->Cys, Glu269-->Cys, or Glu325-->Cys permease is inactive, but oxidation does not restore active lactose transport. The data are consistent with previous observations indicating that Asp237 and Asp240 are not critical for active lactose transport, while Glu126, Glu269, and Glu325 are irreplaceable. Although Glu269-->Cys permease does not transport lactose, the oxidized mutant exhibits significant transport of beta,D-galactosylpyranosyl 1-thio-beta,D-galactopyranoside, a property observed with Glu269-->Asp permease. The observation supports the idea that an acidic residue at position 269 is important for substrate recognition. Finally, oxidized Glu325-->Cys permease catalyzes equilibrium exchange with an apparent pKa of about 6.5, more than a pH unit lower than that observed with Glu325-->Asp permease, thereby providing strong confirmatory evidence that a negative charge at position 325 determines the rate of translocation of the ternary complex between the permease, substrate, and H+.     

Intramolecular charge transfer in the bacteriorhodopsin mutants Asp85-->Asn and Asp212-->Asn: effects of pH and anions

The photovoltage kinetics of the bacteriorhodopsin mutants Asp212-->Asn and Asp85-->Asn after excitation at 580 nm have been investigated in the pH range from 0 to 11. With the mutant Asp85-->Asn (D85N) at pH 7 no net charge translocation is observed and the signal is the same, both in the presence of Cl- (150 mM) and in its absence (75 mM SO4(2-)). Under both conditions the color of the pigment is blue (lambda max = 615 nm). The time course of the photovoltage kinetics is similar to that of the acid-blue form of wild-type, except that an additional transient charge motion occurs with time constants of 60 microseconds and 1.3 ms, indicating the transient deprotonation and reprotonation of an unknown group to and from the extracellular side of the membrane. It is suggested that this is the group XH, which is responsible for proton release in wild-type. At pH 1, the photovoltage signal of D85N changes upon the addition of Cl- from that characteristic for the acid-blue state of wild-type to that characteristic for the acid-purple state. Therefore, the protonation of the group at position at 85 is necessary, but not sufficient for the chloride-binding. At pH 11, well above the pKa of the Schiff base, there is a mixture of "M-like" and "N-like" states. Net proton transport in the same direction as in wild-type is restored in D85N from this N-like state. With the mutant Asp212-->Asn (D212N), time-resolved photovoltage measurements show that in the absence of halide ions the signal is similar to that of the acid-blue form of wild-type and that no net charge translocation occurs in the entire pH range from 0 to 11. Upon addition of Cl- in the pH range from 3.8 to 7.2 the color of the pigment returns to purple and the photovoltage experiments indicate that net proton pumping is restored. However, this Cl(-)-induced activation of net charge-transport in D212N is only partial. Outside this pH range, no net charge transport is observed even in the presence of chloride, and the photovoltage shows the same chloride-dependent features as those accompanying the acid-blue to acid-purple transition of the wild-type.     

Hydrophobic core substitutions in calbindin D9k: effects on stability and structure

The effects of hydrophobic core mutations on the stability and structure of the four-helix calcium-binding protein, calbindin D9k, have been investigated. Eleven mutations involving eight residues distributed within the hydrophobic core of calbindin D9k were examined. Stabilities were measured by denaturant and thermal induced unfolding monitored by circular dichroism spectroscopy. The mutations were found to exert large effects on the stability with midpoints in the urea induced unfolding varying from 1.8 M for Leu23 --> Gly up to 6.6 M for Val70 --> Leu and free energies of unfolding in the absence of denaturant ranging from 6.6 to 27.4 kJ/mol for the Phe66 --> Ala mutant and the wild-type, respectively. A significant correlation was found between the difference in free energy of unfolding (Delta Delta GNU) and the change in the surface area of the side chain caused by the mutation, in agreement with other studies. Notably, both increases and decreases in side-chain surface area caused quantitatively equivalent effects on the stability. In other words, a correlation between the absolute value of the change in the surface of the side chain and Delta DeltaGNU was observed with a value of approximately 0.14 kJ M-1 A-2. The generality of this observation is discussed. Significant effects on the cooperativity of the unfolding reaction were also observed. However, a correlation between the cooperativity and Delta Delta GNU, which has been reported in other systems as an indication of effects of mutations on the unfolded state, was not observed for calbindin D9k. Despite the large effects on Delta Delta GNU and cooperativity, the structures of the mutants in the native form remained intact as indicated by circular dichroism, NMR, and fluorescence measurements. The structural response to calcium-binding was also conserved. The following paper in this issue [Kragelund, B. B., et al. (1998) Biochemistry 37, 8926-8937] examines the effects of these mutations on the calcium binding properties of calbindin D9k.     

Neutralization of conservative charged transmembrane residues in the Na+/glucose cotransporter SGLT1

Our goal was to identify pairs of charged residues in the membrane domains of the Na+/glucose cotransporter (SGLT1) that form salt bridges, to obtain information about packing of the transmembrane helices. The strategy was to neutralize Glu225, Asp273, Asp294, and Lys321 in helices 6-8, express the mutants in oocytes, measure [14C]-alphaMDG uptake, and then attempt to find second-site mutations of opposite charge that restored function. alphaMDG uptake by E225A was identical to that by SGLT1, whereas transport was reduced by over 90% for D273A, D294A, and K321A and was not restored in the double mutants D273A/K321A or D294A/K321A. This suggested that K321 did not form salt bridges with D273 or D294 and that E225 was not involved in salt-bridging. Neutralization of K321 dramatically changed the Na+ uniport and Na+/glucose cotransport kinetics. The maximum rate of uniport in K321A increased 3-5-fold with a decrease in the apparent affinity for Na+ (70 vs 3 mM) and no change in apparent H+ affinity (0.5 microM). The change in Na+ affinity caused a +50 mV shift in the charge/voltage (Q/V) and relaxation time constant (tau)/voltage curves in the presteady-state kinetics. The presteady-state kinetics in H+ remained unchanged. The lower Na+ affinity resulted also in a 200-fold increase in the apparent K0.5 for alphaMDG and phlorizin. Replacements of K321 with alanine, valine, glutamine, arginine, or glutamic acid residues changed the steady-state kinetics in a similar way. Therefore, we suggest that K321 determines, directly or indirectly, (i) the rate and selectivity of SGLT1 uniport activity and (ii) the apparent affinities of SGLT1 for Na+, and indirectly sugar in the cotransport mode.     

Cooperativity of folding of the apomyoglobin pH 4 intermediate studied by glycine and proline mutations

The apomyoglobin pH 4 folding intermediate contains the A, G, and H helices of myoglobin. Helix destabilizing mutations in the A and G helices are used to test whether the pH 4 folding intermediate of apomyoglobin folds cooperatively. Single glycine or proline mutations destabilize the intermediate substantially, showing that intrinsic helix propensities are important for stability of the intermediate. The A and G helices interact to stabilize each other, as shown by the effect of mutations in the G helix on the unfolding of the A helix, which can be monitored by tryptophan fluorescence. Wild type and the most stable mutant unfold in a two-state reaction, as shown by superposition of the unfolding curves measured by two probes (far-UV circular dichroism and Trp fluorescence), while unfolding of the less stable mutants is more complex. Cooperativity and stability of folding are linked also when stabilizing anions (sulphate, perchlorate) are used to adjust stability.     

Proton transfer from Asp-96 to the bacteriorhodopsin Schiff base is caused by a decrease of the pKa of Asp-96 which follows a protein backbone conformational change

In the bacteriorhodopsin photocycle the transported proton crosses the major part of the hydrophobic barrier during the M to N reaction; in this step the Schiff base near the middle of the protein is reprotonated from D96 located near the cytoplasmic surface. In the recombinant D212N protein at pH > 6, the Schiff base remains protonated throughout the photocycle [Needleman, Chang, Ni, Váró, Fornés, White, & Lanyi (1991) J. Biol. Chem. 266, 11478-11484]. Time-resolved difference spectra in the visible and infrared are described by the kinetic scheme BR-->K<==>L<==>N (-->N')-->BR. As evidenced by the large negative 1742-cm-1 band of the COOH group of the carboxylic acid, deprotonation of D96 in the N state takes place in spite of the absence of the unprotonated Schiff base acceptor group of the M intermediate. Instead of internal proton transfer to the Schiff base, the proton is released to the bulk, and can be detected with the indicator dye pyranine during the accumulation of N'. The D212N/D96N protein has a similar photocycle, but no proton is released. As in wild-type, deprotonation of D96 in the N state is accompanied by a protein backbone conformational change indicated by characteristic amide I and II bands. In D212N the residue D96 can thus deprotonate independent of the Schiff base, but perhaps dependent on the detected protein conformational change. This could occur through increased charge interaction between D96 and R227 and/or increased hydration near D96. We suggest that the proton transfer from D96 to the Schiff base in the wild-type photocycle is driven also by such a decrease in the pKa of D96.     

Electrostatic contributions to the binding free energy of the lambdacI repressor to DNA

A model based on the nonlinear Poisson-Boltzmann (NLPB) equation is used to study the electrostatic contribution to the binding free energy of the lambdacI repressor to its operator DNA. In particular, we use the Poisson-Boltzmann model to calculate the pKa shift of individual ionizable amino acids upon binding. We find that three residues on each monomer, Glu34, Glu83, and the amino terminus, have significant changes in their pKa and titrate between pH 4 and 9. This information is then used to calculate the pH dependence of the binding free energy. We find that the calculated pH dependence of binding accurately reproduces the available experimental data over a range of physiological pH values. The NLPB equation is then used to develop an overall picture of the electrostatics of the lambdacI repressor-operator interaction. We find that long-range Coulombic forces associated with the highly charged nucleic acid provide a strong driving force for the interaction of the protein with the DNA. These favorable electrostatic interactions are opposed, however, by unfavorable changes in the solvation of both the protein and the DNA upon binding. Specifically, the formation of a protein-DNA complex removes both charged and polar groups at the binding interface from solvent while it displaces salt from around the nucleic acid. As a result, the electrostatic contribution to the lambdacI repressor-operator interaction opposes binding by approximately 73 kcal/mol at physiological salt concentrations and neutral pH. A variety of entropic terms also oppose binding. The major force driving the binding process appears to be release of interfacial water from the protein and DNA surfaces upon complexation and, possibly, enhanced packing interactions between the protein and DNA in the interface. When the various nonelectrostatic terms are described with simple models that have been applied previously to other binding processes, a general picture of protein/DNA association emerges in which binding is driven by the nonpolar interactions, whereas specificity results from electrostatic interactions that weaken binding but are necessary components of any protein/DNA complex.     

Site-directed mutagenesis of the active site glutamate in human matrilysin: investigation of its role in catalysis

Glu-198 of human matrilysin is a conserved residue in the matrix metalloproteinases and is considered to play an important role in catalysis by acting as a general base catalyst toward the zinc-bound water molecule, on the basis of mechanistic proposals for other zinc proteases. In the present study, Glu-198 is mutated into Asp, Cys, Gln, and Ala, and the zinc binding properties, kinetic parameters, and pH dependence of each mutant are determined in order to examine the role of Glu-198 in catalysis. The mutations chosen either modify (C and D) or eliminate (A and Q) the general base properties of residue-198. All the mutants bind 2 mol of zinc per mol of enzyme, indicating that Glu-198 is not crucial to the binding of the catalytic zinc to the enzyme. The value of kcat/Km for the E198D mutant is only 4-fold lower than that of wild-type enzyme at the pH optimum of 7.5, while that for the E198C mutant is decreased by 160-fold. The E198Q and E198A enzymes containing the mutations that have eliminated the nucleophilic and acid/base properties of the residue are still active, having lower kcat/Km values of 590- and 1900-fold, respectively. The decrease in activity of all the mutants is essentially due to a decrease in kcat. The kcat/Km values of the mutants as a function of pH display broad bell-shaped curves that are similar to the wild-type enzyme. The acidic pKa value is not greatly affected by the change in the chemical properties of residue-198. The similarity in the pH profiles for the mutant and wild-type enzymes indicates that the ionization of Glu-198 is not responsible for the acidic pKa. Ionization of the zinc-bound water may be responsible for this pKa since the three His ligands and the scaffolding of the matrilysin catalytic zinc site are different from that observed in carboxypeptidase A and would predict a lower pKa for the metal-bound water. If the zinc-bound water is the nucleophile in the reaction, the role of Glu-198 in catalysis may be to stabilize the transition state or act as a general acid catalyst after the rate-determining step.     

Surface salt bridges stabilize the GCN4 leucine zipper

We present a study of the role of salt bridges in stabilizing a simplified tertiary structural motif, the coiled-coil. Changes in GCN4 sequence have been engineered that introduce trial patterns of single and multiple salt bridges at solvent exposed sites. At the same sites, a set of alanine mutants was generated to provide a reference for thermodynamic analysis of the salt bridges. Introduction of three alanines stabilizes the dimer by 1.1 kcal/mol relative to the wild-type. An arrangement corresponding to a complex type of salt bridge involving three groups stabilizes the dimer by 1.7 kcal/ mol, an apparent elevation of the melting temperature relative to wild type of about 22 degrees C. While identifying local from nonlocal contributions to protein stability is difficult, stabilizing interactions can be identified by use of cycles. Introduction of alanines for side chains of lower helix propensity and complex salt bridges both stabilize the coiled-coil, so that combining the two should yield melting temperatures substantially higher than the starting species, approaching those of thermophilic sequences.     

Green fluorescent protein as a noninvasive intracellular pH indicator

It was found that the absorbance and fluorescence of green fluorescent protein (GFP) mutants are strongly pH dependent in aqueous solutions and intracellular compartments in living cells. pH titrations of purified recombinant GFP mutants indicated >10-fold reversible changes in absorbance and fluorescence with pKa values of 6.0 (GFP-F64L/S65T), 5.9 (S65T), 6.1 (Y66H), and 4.8 (T203I) with apparent Hill coefficients of 0.7 for Y66H and approximately 1 for the other proteins. For GFP-S65T in aqueous solution in the pH range 5-8, the fluorescence spectral shape, lifetime (2.8 ns), and circular dichroic spectra were pH independent, and fluorescence responded reversibly to a pH change in <1 ms. At lower pH, the fluorescence response was slowed and not completely reversed. These findings suggest that GFP pH sensitivity involves simple protonation events at a pH of >5, but both protonation and conformational changes at lower pH. To evaluate GFP as an intracellular pH indicator, CHO and LLC-PK1 cells were transfected with cDNAs that targeted GFP-F64L/S65T to cytoplasm, mitochondria, Golgi, and endoplasmic reticulum. Calibration procedures were developed to determine the pH dependence of intracellular GFP fluorescence utilizing ionophore combinations (nigericin and CCCP) or digitonin. The pH sensitivity of GFP-F64L/S65T in cytoplasm and organelles was similar to that of purified GFP-F64L/S65T in saline. NH4Cl pulse experiments indicated that intracellular GFP fluorescence responds very rapidly to a pH change. Applications of intracellular GFP were demonstrated, including cytoplasmic and organellar pH measurement, pH regulation, and response of mitochondrial pH to protonophores. The results establish the application of GFP as a targetable, noninvasive indicator of intracellular pH.     

Pathway of proton transfer in bacterial reaction centers: role of aspartate-L213 in proton transfers associated with reduction of quinoneto dihydroquinone

The role of Asp-L213 in proton transfer to reduced quinone QB in the reaction center (RC) from Rhodobacter sphaeroides was studied by site-directed replacement of Asp with residues having different proton donor properties. Reaction centers (RCs) with Asn, Leu, Thr, and Ser at L213 had greatly reduced (approximately 6000-fold) proton-coupled electron transfer [kAB(2)] and proton uptake rates associated with the second electron reduction of QB (QA- QB- + 2H(+)-->QAQBH2) compared to native RCs. RCs containing Glu at L213 showed faster (approximately 90-fold) electron and proton transfer rates than the other mutant RCs but were still reduced (approximately 70-fold) compared with native RCs. These results show that kAB(2) is larger when a carboxylic acid occupies the L213 site, consistent with the proposal that Asp-L213 is a component of a proton transfer chain. The reduced kAB(2) observed with Glu versus Asp at L213 suggests that Asp at L213 is important for proton transfer for some other reason in addition to its proton transfer capabilities. Glu-L213 is estimated to have a higher apparent pKa (pKa > or = 7) than Asp-L213 (pKa < or = 4), as indicated by the slower rate of charge recombination (D+QAQB(-)-->DQAQB) in the mutant RCs. The importance of the pKa and charge of the residue at L213 for proton transfer are discussed. Based on these studies, a model for proton transfer is proposed in which Asp-L213 contributes to proton transfer in native RCs in two ways: (1) it is a component of a proton transfer chain connecting the buried QB molecule with the solvent and/or (2) it provides a negative charge that stabilizes a proton on or near QB.