A tree-based algorithm for determining the effects of solvation on the structure of salivary gland tripeptide NH3+-D-PHE-D-GLU-GLY-COO-

A D-enantiomeric analog of the submandibular gland rat-1 tripeptide FEG (Seq: NH(3)(+)-Phe-Glu-Gly-COO(-)) called feG (Seq: NH(3)(+)-D-Phe-D-Glu-Gly-COO(-)) was examined by molecular dynamics simulations in water. Previous in vacuo simulations suggested a conformation consisting predominantly of interactions between the Phe side chain and glutamyl-carboxyl group and a carboxyl/amino termini interaction. The solvated peptide was simulated using two approaches which were compared-a single 400-ns simulation and a "simulation tree." The "tree" approach utilized 45 10-ns simulations with different conformations used as initial structures for given trajectories. We demonstrate that multiple short duration simulations are able to describe the same conformational space as that described by longer simulations. Furthermore, previously described in vacuo interactions were confirmed with amendments: the previously described head-to-tail arrangement of the amino and carboxyl termini, was not observed; the interaction between the glutamyl carboxyl and Phe side chain describes only one of a continuum of conformations present wherein the aromatic residue remains in close proximity to the glutamyl carbonyl group, and also interacts with either of the two available carboxyl groups. Finally, utilizing only two separate 10-ns trajectories, we were able to better describe the conformational space than a single 60-ns trajectory, realizing a threefold decrease in the computational complexity of the problem.     

Features of the pp60v-src carboxyl terminus that are required for transformation.

Analysis of the biological and biochemical activities of pp60recombinant-src proteins encoded by 12 carboxyl-terminal mutants showed that a wide family of alternate src carboxyl termini permit complete transforming and kinase activities. src proteins having carboxyl termini which are up to 10 amino acids longer than that of pp60c-src (17 amino acids longer than that of pp60v-src) still permit transformation. Transformation-positive mutations preserve leucine-516, a residue which is highly conserved in protein-tyrosine kinase sequences; removal causes in vivo protein instability. Successive deletion mutants show that this residue is at the boundary of a region required for kinase activity. pp60src which is truncated just outside this point still transforms cells and binds both pp50 and pp90 cellular proteins.     

Accurate mass comparison coupled with two endopeptidases enables identification of protein termini

Protein termini play important roles in biological processes, but there have been few methods for comprehensive terminal proteomics. We have developed a new method that can identify both the amino and the carboxyl termini of proteins. The method independently uses two proteases, (lysyl endopeptidase) Lys-C and peptidyl-Lys metalloendopeptidase (Lys-N), to digest proteins, followed by LC-MS/MS analysis of the two digests. Terminal peptides can be identified by comparing the peptide masses in the two digests as follows: (i) the amino terminal peptide of a protein in Lys-C digest is one lysine residue mass heavier than that in Lys-N digest; (ii) the carboxyl terminal peptide in Lys-N digest is one lysine residue mass heavier than that in Lys-C digest; and (iii) all internal peptides give exactly the same molecular masses in both the Lys-C and the Lys-N digest, although amino acid sequences of Lys-C and Lys-N peptides are different (Lys-C peptides end with lysine, whereas Lys-N peptides begin with lysine). The identification of terminal peptides was further verified by examining their MS/MS spectra to avoid misidentifying pairs as termini. In this study, we investigated the usefulness of this method using several protein and peptide mixtures. Known protein termini were successfully identified. Acetylation on N-terminus and protein isoforms, which have different termini, was also determined. These results demonstrate that our new method can confidently identify terminal peptides in protein mixtures.     

Interaction between the noncatalytic region of Sid1p kinase and Cdc14p is required for full catalytic activity and localization of Sid1p

Sid1p is a group II p21-activated kinase/germinal center kinase family member that is part of a signaling network required for cytokinesis in fission yeast. Germinal center kinases are characterized by well conserved amino-terminal catalytic domains followed by less conserved carboxyl termini. The carboxyl termini among group I germinal center kinases are moderately conserved and thought to be regulatory regions. Little is known about the carboxyl termini of group II family members. Sid1p has been shown to bind the novel protein Cdc14p; however, the functional significance of this interaction is unknown. Here we report that the carboxyl terminus of Sid1p is an essential regulatory region. Our results indicate that this region contains the binding domain for Cdc14p, and this association is required for full Sid1p catalytic activity as well as intracellular localization. Furthermore, overexpression of the carboxyl terminus of Sid1p alone compromises the signaling of cytokinesis. We conclude that Cdc14p positively regulates the Sid1p kinase by binding the noncatalytic carboxyl-terminal region of the protein.     

Nuclear magnetic resonance titration curves of histidine ring protons. Conformational transition affecting three of the histidine residues of ribonuclease.

NMR titration curves are reported for the 4 histidine residues of ribonuclease A in sodium acetate and for ribonuclease S in sodium acetate, phosphate, and sulfate solutions. Evidence is presented that the imidazole side chain of histidine residue 48 undergoes a conformational change, probably also involving the carboxyl side chain of aspartic acid residue 14. This group is considered to be responsible for the low pH inflection with pKa 4.2 present in the NMR titration curve of the C-2 proton resonance of histidine 48. The NMR titration curves of the active site histidine residues 12 and 119 also exhibit inflections at low pH values, although there is no carboxyl group within 9 A of the imidazole side chain of histidine residue 12 in the structure of ribonuclease S determined by x-ray crystallography (Wyckoff, H. W., Tsernoglou, D., Hanson, A. W. Knox, J. R., Lee, B., and Richards, F. M. (1970) J. Biol. Chem. 245, 305-328). Curve fitting was carried out on 11 sets of NMR titration data using a model in which the 3 histidine residues 12, 119, and 48 are assumed to be affected by a common carboxyl group. The results obtained indicate that such a model with fewer parameters gives as good a representation of the data as the model in which each histidine residue is assumed to interact separately with a different carboxyl group. Therefore, it is concluded that the ionization of aspartic acid residue 14 is indirectly experienced by the active site histidine residues through the conformational change at histidine 48. A model assuming mutual interaction of the active site histidine residues does not account for the low pH inflections in these curves.     

Identification of regulatory domains in ADP-ribosyltransferase-1 that determine transferase and NAD glycohydrolase activities

Mono-ADP-ribosyltransferases (ART1-7) transfer ADP-ribose from NAD+ to proteins (transferase activity) or water (NAD glycohydrolase activity). The mature proteins contain two domains, an alpha-helical amino terminus and a beta-sheet-rich carboxyl terminus. A basic region in the carboxyl termini is encoded in a separate exon in ART1 and ART5. Structural motifs are conserved among ART molecules. Successive amino- or carboxyl-terminal truncations of ART1, an arginine-specific transferase, identified regions that regulated transferase and NAD glycohydrolase activities. In mouse ART1, amino acids 24-38 (ART-specific extension) were needed to inhibit both activities; amino acids 39-45 (common ART coil) were required for both. Successive truncations of the alpha-helical region reduced transferase and NAD glycohydrolase activities; however, truncation to residue 106 enhanced both. Removal of the carboxyl-terminal basic domain decreased transferase, but enhanced NAD glycohydrolase, activity. Thus, amino- and carboxyl-terminal regions of ART1 are required for transferase activity. The enhanced glycohydrolase activity of the shorter mutants indicates that sequences, which are not part of the NAD binding, core catalytic site, exert structural constraints, modulating substrate specificity and catalytic activity. These functional domains, defined by discrete exons or structural motifs, are found in ART1 and other ARTs, consistent with conservation of structure and function across the ART family.     

Specificity and enzyme kinetics of the quorum-quenching N-Acyl homoserine lactone lactonase (AHL-lactonase)

N-Acyl homoserine lactone (AHL) quorum-sensing signals are the vital elements of bacterial quorum-sensing systems, which regulate diverse biological functions, including virulence. The AHL-lactonase, a quorumquenching enzyme encoded by aiiA from Bacillus sp., inactivates AHLs by hydrolyzing the lactone bond to produce corresponding N-acyl homoserines. To characterize the enzyme, the recombinant AHL-lactonase and its four variants were purified. Kinetic and substrate specificity analysis showed that AHL-lactonase had no or little residue activity to non-acyl lactones and noncyclic esters, but displayed strong enzyme activity toward all tested AHLs, varying in length and nature of the substitution at the C3 position of the acyl chain. The data also indicate that the amide group and the ketone at the C1 position of the acyl chain of AHLs could be important structural features in enzyme-substrate interaction. Surprisingly, although carrying a (104)HX- HXDH(109) short sequence identical to the zinc-binding motif of several groups of metallohydrolytic enzymes, AHL-lactonase does not contain or require zinc or other metal ions for enzyme activity. Except for the amino acid residue His-104, which was shown previously to not be required for catalysis, kinetic study and conformational analysis using circular dichroism spectrometry showed that substitution of the other key residues in the motif (His-106, Asp-108, and His-109), as well as His-169 with serine, respectively, caused conformational changes and significant loss of enzyme activity. We conclude that AHL-lactonase is a highly specific enzyme and that the (106)HXDH(109) approximately H(169) of AHL-lactonase represents a novel catalytic motif, which does not rely on zinc or other metal ions for activity.     

Chemical properties of the histidine residue of secretin: evidence for a specific intramolecular interaction

Secretin has a single histidine residue located at the amino terminus which plays a crucial role in its biological activity. The chemical properties, viz. pK and reactivity, of the alpha-amino and imidazole groups of this residue were determined at a secretin concentration of 10(-6) M in 0.1 M KCl at 37 degrees C. Competitive labelling using tritiated 1-fluoro-2,4-dinitrobenzene (DNP-F) as the labelling reagent was the experimental approach employed. The alpha-amino group was found to have a pK value of 8.83 and a reactivity 5-times that of the alpha-amino group in the model compound, histidylglycine. For the imidazole function a pK value of 8.24 and a reactivity 26-times that of the imidazole function in histidylglycine was found. Both these groups in secretin had pK values which were shifted one pK unit higher than in histidylglycine, but like the model compound the reactivity of the imidazole function was still linked to the state of ionization of the alpha-amino group. These observations are interpreted as evidence for the existence of a major conformational state in dilute aqueous solution in which the amino-terminal histidine of secretion is interacting with a negatively charged carboxyl group.     

Studies on the activity and activation of rat liver microsomal glutathione transferase with a series of glutathione analogues

The substrate specificity of rat liver microsomal glutathione transferase toward glutathione has been examined in a systematic manner. Out of a glycyl-modified and eight gamma-glutamyl-modified glutathione analogues, it was found that four (glutaryl-L-Cys-Gly, alpha-L-Glu-L-Cys-Gly, alpha-D-Glu-L-Cys-Gly, and gamma-L-Glu-L-Cys-beta-Ala) function as substrates. The kinetic parameters for three of these substrates (the alpha-D-Glu-L-Cys-Gly analogue gave very low activity) were compared with those of GSH with both unactivated and the N-ethylmaleimide-activated microsomal glutathione transferase. The alpha-L-Glu-L-Cys-Gly analogue is similar to GSH in that it has a higher kcat (6.9 versus 0.6 s-1) value with the activated enzyme compared with the unactivated enzyme but displays a high Km (6 versus 11 mM) with both forms. Glutaryl-L-Cys-Gly, in contrast, exhibited a similar kcat (8.9 versus 6.7 s-1) with the N-ethylmaleimide-treated enzyme but retains a higher Km value (50 versus 15 mM). Thus, the alpha-amino group of the glutamyl residue in GSH is important for the activity of the activated microsomal glutathione transferase. These observations were quantitated by analyzing the changes in the Gibbs free energy of binding calculated from the changes in kcat/Km values, comparing the analogues to GSH and each other. It is estimated that the binding energy of the alpha-amino group of the glutamyl residue in GSH contributes 9.7 kJ/mol to catalysis by the activated enzyme, whereas the corresponding value for the unactivated enzyme is 3.2 kJ/mol. The importance of the acidic functions in glutathione is also evident as shown by the lack of activity with 4-aminobutyric acid-L-Cys-Gly and the low kcat/Km values with gamma-L-Glu-L-Cys-beta-Ala (0.03 and 0.01 mM-1s-1 for unactivated and activated enzyme, respectively). Utilization of binding energy from a correctly positioned carboxyl group in the glycine residue (10 and 17 kJ/mol for unactivated and activated enzyme, respectively) therefore also appears to be required for optimal activity and activation. A conformational change in the microsomal glutathione transferase upon treatment with N-ethylmaleimide or trypsin, which allows utilization of binding energy from the alpha-amino group of GSH as well as the glycine carboxyl in catalysis, is suggested to account for at least part of the activation of the enzyme.     

Crystal structure of Escherichia coli CheY refined at 1.7-A resolution

The three-dimensional structure of wild-type CheY from Escherichia coli has been refined by stereochemically restrained least squares minimization to a crystallographic R-factor of 15.1% at 1.7-A resolution. The structure contains 1165 atoms, including all atoms of the protein, 147 water molecules, and three sulfate ions. The final model has root mean square deviations of 0.018 and 0.049 A from idealized bond lengths and angle distances, respectively. Seven amino acid side chains have been modeled in dual conformations. CheY folds as a compact (beta/alpha)5 globular protein, with the phosphorylation region contained in a cavity on one face of the molecule. This active site area is bordered by the carboxyl termini of the three central beta-strands, by alpha 1, and by the loop connecting beta 5 to alpha 5. The Lys-109 side chain of this loop extends into the active site by virtue of its cis peptide bond conformation preceding Pro-110. The epsilon-amino group of Lys-109 is in close bonding contact with the carboxyl group of Asp-57, the residue that is phosphorylated in the activation process of CheY. The details of the hydrogen bonding network in the phosphorylation region indicate that structural rearrangements must accompany the phosphorylation of Asp-57.     

Critical Arginine Residue for Maintaining the Bacteriophage Tail Structure

The addition of 0.2 m l-arginine to various T-even bacteriophage preparations inactivated the virus preparations irreversibly. The virus particles were even more sensitive to added d-arginine and l-homoarginine than to l-arginine but were unaffected by arginine analogues with either an altered carboxyl group or guanidyl group. Treatment of phage T2H with 2,3-butanedione, a reagent which specifically reacts with the guanidyl portion of arginine residues, resulted in the apparent in-activation of most of the virus particles. However, after incubation of the treated particles at pH 7.5 at 37 C for 1 hr in the absence of butanedione, the original virus titer almost completely returned. The reactivation was completely inhibited by the presence of 0.2 m d-arginine. It appeared that the virus protein coat was sufficiently plastic so that the initial conformational change resulting from the alteration of an arginine residue (to possibly an ornithine residue) was at least partially reversible and that the virus tail proteins then refolded to produce a stable and active virus particle. These reactivated virus particles were not sensitive to inactivation by d-arginine but could now be rapidly inactivated by l-ornithine. Virus particles inactivated by arginine have altered tail structures. They have contracted tail sheaths still attached to tail plates and still contain tail cores. These properties of virus particles indicate that there is a free carboxyl group and a guanidyl group spatially equivalent to an arginine residue on one component of the virus tail which bind reversibly by means of polar linkages to another tail component. These bonds maintain the integrity of the virus tail. Added arginine appears to compete with this endogenous viral arginine for the binding sites and then to favor an irreversible conformational change.     

Structure-activity relationship of endothelin: importance of charged groups

Endothelin (ET)-related peptides including ET-1 (1-39) were synthesized, and their constricting activity in rat pulmonary artery rings and pressor activity in unanesthetized rat were measured to elucidate their structure-activity relationship. The vasoconstrictor activities of ET-2, ET-3 and sarafotoxin S6b were one-half, one-60th and one-third that of ET-1, respectively. Such differences in biological activities should mainly arise from sequence heterogeneity at the N-terminal portion, especially at positions 4 to 7. All of the blocked ETs at the amino or carboxyl termini showed greatly decreased activities. A monocyclic analog, in which Cys3 and Cys11 were replaced by Ala, showed one-third the activity of ET-1; however, its deamino dicarba analog was almost completely inactive. Significant activities were retained even with replacement of amino acids at positions Ser4, Ser5, Leu6, Met7, Lys9, Tyr13, and Trp21 by Ala, Ala, Gly, Met(0), Leu, Phe, and Tyr or Phe, respectively. On the other hand, replacement of Asp8, Glu10 and Phe14 by Asn, Gln and Ala, respectively, resulted in complete loss of the biological activity. These results indicated that two disulfide bonds in ET molecule were not essential for the expression of vasoconstricting activity. Both terminal amino and carboxyl groups, carboxyl groups of Asp8 and Glu10, and the aromatic group of Phe14 seemed to be contributing, more or less, to the expression of the biological activities.     

Comparison of the backbone dynamics of the apo- and holo-carboxy-terminal domain of the biotin carboxyl carrier subunit of Escherichia coli acetyl-CoA carboxylase

The biotin carboxyl carrier protein (BCCP) is a subunit of acetyl-CoA carboxylase, a biotin-dependent enzyme that catalyzes the first committed step of fatty acid biosynthesis. In its functional cycle, this protein engages in heterologous protein-protein interactions with three distinct partners, depending on its state of post-translational modification. Apo-BCCP interacts specifically with the biotin holoenzyme synthetase, BirA, which results in the post-translational attachment of biotin to a single lysine residue on BCCP. Holo-BCCP then interacts with the biotin carboxylase subunit of acetyl-CoA carboxylase, which leads to the addition of the carboxylate group of bicarbonate to biotin. Finally, the carboxy-biotinylated form of BCCP interacts with transcarboxylase in the transfer of the carboxylate to acetyl-CoA to form malonyl-CoA. The determinants of protein-protein interaction specificity in this system are unknown. The NMR solution structure of the unbiotinylated form of an 87 residue C-terminal domain fragment (residue 70-156) of BCCP (holoBCCP87) and the crystal structure of the biotinylated form of a C-terminal fragment (residue 77-156) of BCCP from Escherichia coli acetyl-CoA carboxylase have previously been determined. Comparative analysis of these structures provided evidence for small, localized conformational changes in the biotin-binding region upon biotinylation of the protein. These structural changes may be important for regulating specific protein-protein interactions. Since the dynamic properties of proteins are correlated with local structural environments, we have determined the relaxation parameters of the backbone 15N nuclear spins of holoBCCP87, and compared these with the data obtained for the apo protein. The results indicate that upon biotinylation, the inherent mobility of the biotin-binding region and the protruding thumb, with which the biotin group interacts in the holo protein, are significantly reduced.     

Use of engineered proteins with internal tryptophan reporter groups and pertubation techniques to probe the mechanism of ligand-protein interactions: investigation of the mechanism of calcium binding to calmodulin.

Stopped-flow kinetic and fluorescence spectroscopic analyses, including solvent and temperature perturbations, of five isofunctional structural mutants of calmodulin indicate that calcium binding to calmodulin follows the order site III, site IV, site I, site II, with dissociation occurring in the reverse order. Each of the isofunctional structural mutants contains a single tryptophan residue, introduced by site-specific mutagenesis, as an internal spectroscopic reporter group that was used as a probe of local conformational change. Calcium binding was studied by using flow dialysis or by using fluorescence spectroscopy and monitoring the change in the single tryptophan residue in each calcium-binding site. Calcium removal was examined by using EDTA and monitoring tryptophan fluorescence or by using Quin 2 and monitoring the change in the chromophoric chelator. Computational analysis of the data suggests a rate-limiting step for dissociation between calcium removal from sites I/II and sites III/IV. Unexpected results with the site IV isofunctional mutant (Q135W-CaM) indicated cross-talk between the amino and carboxyl terminal halves of CaM during the calcium-binding mechanism. Studies with ethylene glycol provided empirical data that suggest the functional importance of the electrostatic potential of CaM, or the molarity of water, in the calcium-binding process. Altogether, the data allowed a kinetic extension of the sequential, cooperative model for calcium binding to calmodulin and provided values for additional parameters in the model of calcium binding to CaM, a prototypical member of the family of proteins required for calcium signal transduction in eukaryotic cells.(ABSTRACT TRUNCATED AT 250 WORDS)     

CPSARST: an efficient circular permutation search tool applied to the detection of novel protein structural relationships

Circular permutation of a protein can be visualized as if the original amino- and carboxyl termini were linked and new ones created elsewhere. It has been well-documented that circular permutants usually retain native structures and biological functions. Here we report CPSARST (Circular Permutation Search Aided by Ramachandran Sequential Transformation) to be an efficient database search tool. In this post-genomics era, when the amount of protein structural data is increasing exponentially, it provides a new way to rapidly detect novel relationships among proteins.     

A single activity carboxyl methylates both farnesyl and geranylgeranyl cysteine residues.

Members of the Ras superfamily of small GTP-binding proteins, gamma-subunits of heterotrimeric G proteins and nuclear lamin B are subject to a series of post-translational modifications that produce prenylcysteine methylester groups at their carboxyl termini. The thioether-linked polyisoprenoid substituent can be either farnesyl (C15) or geranylgeranyl (C20). Small molecule prenylcysteine derivatives with either the C15 or C20 modification, such as N-acetyl-S-trans,trans-farnesyl-L-cysteine (AFC), S-trans,trans-farnesylthiopropionate (FTP), as well as the corresponding geranylgeranyl derivatives (AGGC and GGTP) are substrates for the carboxyl methyltransferase. Saccharomyces cerevisiae ste14 mutants that lack RAS and a-factor carboxyl methyltransferase activity are also unable to methylate farnesyl and geranylgeranylcysteine derivatives. Moreover, C20-substituted cysteine analogs directly compete for carboxyl methylation with the C15-substituted cysteine analogs and vice versa. Finally, AGGC is even more effective than AFC as an inhibitor of Ras carboxyl methylation, despite the fact that Ras is methylated at a farnesylcysteine rather than a geranylgeranylcysteine residue.     

The additional methionine residue at the N-terminus of bacterially expressed human interleukin-2 affects the interaction between the N- and C-termini

To gain insight into the origin of the difference in isoelectric point (pI) values for wild-type human interleukin-2 (IL-2) and IL-2 with an additional methionine residue at the N-terminus (Met-IL-2), conformational properties of the two molecular forms of IL-2 were compared by utilizing 1H NMR spectroscopy. Although overall conformations were conserved in the two forms, the presence of the additional methionine residue at the N-terminus induced chemical shift changes for residues Ala1 to Lys8 as well as for Thr133, which is located at the C-terminus. These observations indicate that the effect of the additional methionine residue is confined to the N- and C-terminal regions and unveil the existence of an interaction between the N- and C-terminal regions. The chemical shift change observed for Thr133 can be interpreted in terms of a change in pKa of the C-terminal carboxyl group, which interacts differently with the N-terminal amino group in the two forms of IL-2. It seems to be reasonable to conclude that the difference in pI values for the two forms of IL-2 is the consequence of the different interactions between the C- and N-terminal residues.     

Structural and energetic consequences of disruptive mutations in a protein core.

We have characterized the properties of a set of variants of the N-terminal domain of lambda repressor bearing disruptive mutations in the hydrophobic core. These mutations include some that dramatically alter the total core residue volume (by up to six methylene groups) and some that place a single polar residue into the otherwise hydrophobic core. The structural properties of the purified proteins have been studied by CD spectroscopy, biological activity, recognition by conformation-specific monoclonal antibodies, and 1H NMR spectroscopy. The stabilities of the proteins have been measured by thermal and guanidine hydrochloride denaturation. Proteins with disruptive core mutations are found to display a continuum of increasingly nonnative properties. Large internal volume changes cause both significant conformational rearrangements and destabilization by up to 5 kcal/mol. Variants with polar substitutions at core positions no longer behave like well-folded proteins but rather display characteristics of molten globules. However, even proteins bearing some of the most disruptive mutations retain many of the crude secondary and tertiary structural features of the wild-type protein. These results indicate that primitive elements of native structure can form in the absence of normal core packing.     

Conformational studies of chemotactic HCO-Met-Leu-Phe-OMe analogues

Summary In order to investigate the proper peptide backbone conformation able to elicit a biological activity, HCO-Met-Pro-Phe-OMe, HCO-Met-[COO]Leu-Phe-OMe, and HCO-Met-OLeu-Phe-OMe, analogues of the prototypical chemotactic peptide HCO-Met-Leu-Phe-OMe, were studied by CD and IR techniques. The results obtained comparing biological and conformational data evidence the critical presence of (i) the NH group at position 2, (ii) a rather flexible backbone, (iii) the chemical structure of the central residue which can affect the stability of a possible active conformer.     

Interaction between the carboxyl groups of Asp127 and Asp129 in the active site of Escherichia coli phosphofructokinase.

The pH dependence of the enzymic properties of the phosphofructokinase from Escherichia coli was compared to those of two mutants in which one carboxyl group of the active site has been removed from either Asp127 or Asp129. All measurements of activity were made in the presence of allosteric activator ADP or GDP to eliminate any cooperative process. Asp129 is a crucial residue for the activity of phosphofructokinase since its conversion to Ser decreases the catalytic activity by 2-3 orders of magnitude in both the forward and reverse reactions, but the ionization of Asp129 is not directly related the pH dependence of phosphofructokinase activity. This pH dependence is however modified by the Asp129----Ser mutation, which decreases the pK of another residue, Asp127, by as much as pH of 1.5. The side chain of Asp127 has the catalytic role proposed earlier: its deprotonated form acts as a base in the forward reaction, and its protonated form acts as an acid in the reverse reaction. The protonated form of Asp127 is also required for the binding of fructose 1,6-bisphosphate. The electrostatic interaction between the carboxyl groups of Asp127 and Asp129 seems different in free phosphofructokinase to that in enzyme/substrate complexes, suggesting that a conformational change occurs upon substrate binding. The pH dependence of phosphofructokinase activity involves one other ionizable group with a pK of approximately 6 which does not belong to the side chains of Asp127 or Asp129.