Poly(L-histidyl-L-alanyl-α-L-glutamic acid). II. Catalysis of p-nitrophenyl acetate hydrolysis

The hydrolysis of p-nitrophenyl acetate is catalyzed by imidazole, free in solution or as the side chain in poly(His-Ala-Glu). This is based on the observations that the reaction is first order in ester and first order in nonprotonated imidazole. Catalysis of p-nitrophenyl acetate hydrolysis is dependent on solvent conditions. The effect of low concentrations of ethanol, dioxane, and trifluoroethanol were investigated. As the concentration of organic solvent is increased, the second-order rate constant for imidazole catalysis decreases. The decrease, however, is greater for imidazole than for poly(His-Ala-Glu). In 2% trifluoroethanol/water solution, free imidazole has twice the catalytic activity of polymeric imidazole, while in 40% trifluoroethanol/water they have equal activity. Since under the latter solvent conditions poly(His-Ala-Glu) is partially α-helical, the relative improvement in polymeric–imidazole catalysis may be attributed to imidazole hydrogen-bonded to a carboxylate ion. With this assumption the carboxylate–imidazole hydrogen-bonded system has been calculated to have three times the base catalytic activity of imidazole.     

Identification of Catalytically Active Groups of Wheat Germ Lipoxygenase

Dependences of lipoxygenase activity on pH, ionization enthalpies of various chemical groups, photoinactivation of the enzyme, and effects of specific reagents (p-CMB, DPF, and PMSF) were studied to identify catalytically active groups of the enzyme. The data suggest that the catalytic site of lipoxygenase includes imidazole and hydroxyl groups of histidine and serine residues, respectively.     

Mutation of non-conserved amino acids surrounding catalytic site to shift pH optimum of Bacillus circulans xylanase

Improvement of enzyme function by engineering pH dependence of enzymatic activity is of importance for industrial application of Bacillus circulans xylanases. Target mutation sites were selected by structural alignment between B. circulans xylanase and other xylanases having different pH optima. We selected non-conserved mutant sites within 8 Å from the catalytic residues, to see whether these residues have some role in modulating pK_as of the catalytic residues. We hypothesized that the non-conserved residues which may not have any role in enzyme catalysis might perturb pK_as of the catalytic residues. Change in pK_a of a titratable group due to change in electrostatic potential of a mutation was calculated and the change in pH optimum was predicted from the change in pK_a of the catalytic residues. Our strategy is proved to be useful in selection of promising mutants to shift the pH optimum of the xylanases towards desired side.     

A study of the histidine residues of human carbonic anhydrase B using 270 MHz proton magnetic resonance

Nine resonances in the 270 MHz proton magnetic resonance spectrum of human carbonic anhydrase B have been identified with imidazole C₍₂₎ protons of histidine residues, six of which are observed to titrate with pK_a values in the range 4.7 to 7.4. The behaviour of the nine resonances has been studied in the presence of the inhibitors, iodide, cyanide, acetate, hexacyanochromate, and imidazole. Measurements have also been made of the enzyme in its apo, cobalt, and mono-alkylated forms. Used in conjunction with the crystal structure, these results have enabled the tentative assignment of all nine resonances to particular histidine residues in the amino-acid sequence. Three of the active-site histidines at positions 64, 67, and 200 have low pK_a values and cannot be directly linked to the activity of the enzyme. However, the resonances assigned to the three metal-liganding histidines do exhibit changes on anion binding and with pH, which parallel changes in the esterase activity. These results are consistent with the model of an ionizable water molecule bound to the zinc ion.Linewidth measurements of the resonances of the histidine residues on the enzyme surface are used to estimate pseudo-first-order rate constants of the order of 4 × 10³ s^?1 for D⁺ exchange between imidazole N and solvent in the absence of buffer. These rates are observed to increase in the presence of small amounts of the buffers Tris and imidazole.     

1H nuclear magnetic resonance studies of histidine-containing di- and tripeptides. Estimation of the effects of charged groups on the pKa value of the imidazole ring.

The nmr titration curves of chemical shifts versus pH were observed for the protons of various histidine-containing di- and tripeptides. With these results, the macroscopic pK_a values and the chemical shifts intrinsic to each ionic species were determined by a computer curve-fitting based on a simple acid dissociation sequence. The pK_a value of the imidazole ring in N-acetyl-L -histidine methylamide was assumed to represent the intrinsic (or unperturbed) pK_a of the imidazole rings of histidine having peptide linkages at both the CO and NH sides. The pK_a values of the imidazole rings observed for most di- and tripeptides were reasonably reproduced by simple calculations using the intrinsic value and the perturbations due to the CO₂^? and NH₃⁺ groups located at various positions. Some other factors affecting the pK_a value of the imidazole ring are also discussed.     

Identification of Catalytically Active Groups of Penicillium canescens F-436 β-Galactosidase

The functional groups of Penicillium canescens F-436 b-galactosidase have been identified. The pK values and heats of ionization of these groups and photoinactivation of the enzyme with methylene blue indicate that the active site contains carboxyl and imidazole groups. A mechanism for the participation of these groups in the cleavage of the glycoside bond in lactose is proposed.     

Function of four tryptophan residues on Cel7A catalytic efficiency: Insight into cellulase screening strategy based on natural cellulose or cellulose analogs

There is a high level of conservation of tryptophans within the active site architecture of the cellulase family, whereas the function of the four tryptophans in the catalytic domain of Cel7A is unclear. By mutating four tryptophan residues in the catalytic domain of Cel7A from Penicillium piceum (PpCel7A), the binding affinity between PpCel7A and p-nitrophenol-d -cellobioside (pNPC) was reduced as determined by Michaelis–Menten constants, molecular dynamics simulations, and fluorescence spectroscopy. Furthermore, PpCel7A variants showed a reduced level of cellobiohydrolase (CBH) activity against cellulose analogs or natural cellulose. Therefore, it could be concluded four tryptophan residues in Cel7A played a critical role in substrate binding. Mutagenesis results indicated that the W390 stacking interactions at the −2 site played an essential role in facilitating substrate distortion to the −1 site. As soon as the function was altered, the mutation would inevitably affect the catalytic activity against the natural substrate. Interestingly, no clear relationship was found between the CBH activity of PpCel7A variants against pNPC and Avicel. p-Nitrophenol contains many electrophilic groups that may result in overestimation of the binding constant between tryptophan residues and pNPC in comparison with the natural substrate. Consequently, screening improved cellulase using cellulose analogs would divert attention from the target direction for lignocellulose biorefinery. Clarifying mechanism of catalytic diversity on the natural cellulose or cellulose analogs may give better insight into cellulase screening and selecting strategy.     

Mutation of non-conserved amino acids surrounding catalytic site to shift pH optimum of Bacillus circulans xylanase

Improvement of enzyme function by engineering pH dependence of enzymatic activity is of importance for industrial application of Bacillus circulans xylanases. Target mutation sites were selected by structural alignment between B. circulans xylanase and other xylanases having different pH optima. We selected non-conserved mutant sites within 8 Å from the catalytic residues, to see whether these residues have some role in modulating pK_as of the catalytic residues. We hypothesized that the non-conserved residues which may not have any role in enzyme catalysis might perturb pK_as of the catalytic residues. Change in pK_a of a titratable group due to change in electrostatic potential of a mutation was calculated and the change in pH optimum was predicted from the change in pK_a of the catalytic residues. Our strategy is proved to be useful in selection of promising mutants to shift the pH optimum of the xylanases towards desired side.     

Functional Expression of <Emphasis Type="Italic">Spirulina</Emphasis>-Δ<Superscript>6</Superscript> Desaturase Gene in Yeast, <Emphasis Type="Italic">Saccharomyces cerevisiae</Emphasis>

Spirulina-acyl-lipid desaturases are membrane-bound enzymes found in thylakoid and plasma membranes. These enzymes carry out the fatty acid desaturation process of Spirulina to yield γ-linolenic acid (GLA) as the final desaturation product. In this study, Spirulina-Δ⁶ desaturase encoded by the desD gene was heterologously expressed and characterized in Saccharomyces cerevisiae. We then conducted site-directed mutagenesis of the histidine residues in the three histidine boxes to determine the role of these amino acid residues in the enzyme function. Our results showed that while four mutants showed complete loss of Δ⁶-desaturase activity and two mutants showed only trace of the activity, the enzyme activity could be partially restored by chemical rescue using exogenously provided imidazole. This study reveals that the histidine residues (which have imidazole as their functional group) in the conserved clusters play a critical role in Δ⁶-desaturase activity, possibly by providing a di-iron catalytic center. In our previous study, this enzyme was expressed in Escherichia coli. The results reveal that the enzyme can function only in the presence of an exogenous cofactor, ferredoxin, provided in vitro. This evidence suggests that baker’s yeast has a cofactor that can complement ferredoxin, thought to act as an electron donor for the Δ⁶ desaturation in cyanobacteria, including Spirulina. The electron donor of the Spirulina-Δ⁶ desaturation in yeast is more likely to be cytochrome b₅, which is absent in E. coli. This means that the enzyme expressed in S. cerevisiae can catalyze the biosynthesis of the product, GLA, in vivo.     

Application of high-performance liquid chromatography to competitive labeling studies: The chemical properties of functional groups of glucaton

A competitive-labeling study of glucagon was carried out using [³H]- and [¹⁴C]-1-fluoro-2,4-dinitrobenzene to determine simultaneously the chemical properties of the α-amino and imidazole groups of the N-terminal histidine residue, and the lysine and tyrosine residues, under conditions where glucagon is in its physiologically active monomer form. The dinitrophenyl derivatives of these groups were purified by high-performance liquid chromatography which greatly simplified the separation steps of the procedure. The results showed the α-amino and tyrosine groups to have relatively normal behavior, with pK values of 7.98 and 10.22, respectively, while the lysine had a low pK of 8.46. The imidazole function had an apparent pK of 7.84, substantially higher than previous estimates. This difference may be accounted for by the effect of the charged form of the adjacent α-amino group on the nucleophilicity of the imidazole group.     

Structural characterization of 2,6‐dichloro‐p‐hydroquinone 1,2‐dioxygenase (PcpA) from Sphingobium chlorophenolicum,a new type of aromatic ring‐cleavage enzyme

PcpA (2,6‐dichloro‐p‐hydroquinone 1,2‐dioxygenase) from Sphingobium chlorophenolicum, a non‐haem Fe(II) dioxygenase capable of cleaving the aromatic ring of p‐hydroquinone and its substituted variants, is a member of the recently discovered p‐hydroquinone 1,2‐dioxygenases. Here we report the 2.6 Å structure of PcpA, which consists of four βαβββ motifs, a hallmark of the vicinal oxygen chelate superfamily. The secondary co‐ordination sphere of the Fe(II) centre forms an extensive hydrogen‐bonding network with three solvent exposed residues, linking the catalytic Fe(II) to solvent. A tight hydrophobic pocket provides p‐hydroquinones access to the Fe(II) centre. The p‐hydroxyl group is essential for the substrate‐binding, thus phenols and catechols, lacking a p‐hydroxyl group, do not bind to PcpA. Site‐directed mutagenesis and kinetic analysis confirm the critical catalytic role played by the highly conserved His10, Thr13, His226 and Arg259. Based on these results, we propose a general reaction mechanism for p‐hydroquinone 1,2‐dioxygenases.     

Catalytic activity of linear-, cyclic- and polypeptides having a sequence of -Asp-β-Ala-Gly-ser-β-Ala-Gly-His-β-Ala-Gly in the hydrolysis of p-nitrophenyl acetate

Catalytic activities of Boc-Asp-β-Ala-Gly-Ser-β-Ala-Gly-His-β-Ala-Gly-OEt(Boc-9-Oet), Boc-Asp-β-Ala-Gly-Ser-β-Ala-Gly-His-β-Ala-Gly-OH(Boc-9-OH), cyclo(Asp-β-Ala-Gly-Ser-β-Ala-Gly-His-β-Ala-Gly) (Cyclic 9) and poly(Asp-β-Ala-Gly-Ser-β-Ala-Gly-His-β-Ala-Gly) (Poly 9) in the Hydrolysis of p-nitrophenyl acetate were investigated in detail and compared with each other and with poly(His-β-Ala-Gly) (Poly 3) which has no Ser and Asp residues. Generally, Poly 3 was less active than the others, which contain Ser and Asp residues together with the His residue. The reaction rate-substrate concentration for Boc-9-OEt, Boc-9-OH. Cyclic 9 and Poly 3 gave straight lines, while that for Poly 9 showed slightly the tendency of saturation at high substrate concentration. The reaction rates were all proportional to the concentration of the peptides. All peptides gave similar, sigmoid-type pH-k_cat profiles. The pK values obtained from these pH-k_cat profiles agreed fairly well with those of histidine residues obtained from ¹³C n.m.r. chemical shifts, which suggests that the predominant participating functional group in the catalytic reaction is the imidazole group in the histidine residue. The pK values of the His residue in peptides with the -Asp-β-Ala-Gly-Ser-β-Ala-Gly-His-β-Ala-Gly- sequence were shifted to higher pH region compared with Poly 3, suggesting that the effect of the carboxyl group in the Asp residue and the extents of pK-shift for linear peptides were larger than for Cyclic 9 or Poly 9. The catalytic reaction rates by Boc-9-OEt or Cyclic 9 increased steadily with increase in temperature, while the reaction rate-temperature profiles for Poly 9 and Poly 3 gave the optimum temperatures at around 40–50°C.     

Chemical modification of a catalytic antibody that accelerates the hydrolysis of carbonate esters

Catalytic antibody, 4A1, catalyzes the hydrolysis of p-nitrophenyl alkyl carbonate. To determine the amino acid residues related to the catalytic activity of the antibody, we studied the effect of Tyr-, Trp-, and Lys-selective reagents on the catalytic activity and determined the amino acid sequences around the modified amino acid residues. We found that the Tyr-selective reagent is the most effective one and the modification of one Tyr residue results in the complete loss of the catalytic activity. The modified Tyr residue is identified to be Tyr-32 in the CDR-1 of the L chain.     

Functional tuning of the catalytic residue pKa in a de novo designed esterase

AlleyCatE is a de novo designed esterase that can be allosterically regulated by calcium ions. This artificial enzyme has been shown to hydrolyze p‐nitrophenyl acetate (pNPA) and 4‐nitrophenyl‐(2‐phenyl)‐propanoate (pNPP) with high catalytic efficiency. AlleyCatE was created by introducing a single‐histidine residue (His¹⁴⁴) into a hydrophobic pocket of calmodulin. In this work, we explore the determinants of catalytic properties of AlleyCatE. We obtained the pK_a value of the catalytic histidine using experimental measurements by NMR and pH rate profile and compared these values to those predicted from electrostatics pK_a calculations (from both empirical and continuum electrostatics calculations). Surprisingly, the pK_a value of the catalytic histidine inside the hydrophobic pocket of calmodulin is elevated as compared to the model compound pK_a value of this residue in water. We determined that a short‐range favorable interaction with Glu¹²⁷ contributes to the elevated pK_a of His¹⁴⁴. We have rationally modulated local electrostatic potential in AlleyCatE to decrease the pK_a of its active nucleophile, His¹⁴⁴, by 0.7 units. As a direct result of the decrease in the His¹⁴⁴ pK_a value, catalytic efficiency of the enzyme increased by 45% at pH 6. This work shows that a series of simple NMR experiments that can be performed using low field spectrometers, combined with straightforward computational analysis, provide rapid and accurate guidance to rationally improve catalytic efficiency of histidine‐promoted catalysis. Proteins 2017; 85:1656–1665. © 2017 Wiley Periodicals, Inc.     

Catalytic amidolysis of amino acid p-nitroanilides using transition state analogue imprinted artificial enzymes: Cooperative effect of pyridine moiety

Enzyme-like polymer catalysts with the imprints of phosphonate transition state analogue (TSA) lined along with imidazole and pyridine moieties were synthesized using methacryloyl-l-histidine and 4-vinylpyridine as the functional monomers and phenyl-1-(N-benzyloxycarbonylamino)-2-(phenyl)ethyl phosphonate – the TSA of hydrolytic reaction as the template for the amidolysis of N-benzyloxycarbonyl-l-phenylalanine p-nitroanilide (Z-l-Phe-PNA). Polymers containing different functional groups can act together to provide catalytic activity and selectivity superior to what can be obtained from monofunctional analogues. The higher rate acceleration exhibited by the bifunctional polymer over the monofunctional polymers indicates cooperative catalysis of imidazole and pyridine moieties. The optimum catalytic competence is shown by the bifunctional polymer containing imidazole and pyridine moieties in 2:1 M ratio which may be due to alignment of the functional groups in proper H-bond distance. In addition to the non-covalent interactions like hydrogen bonding or π-stacking interactions between the functional groups of the polymer and the template, 3D-microcavities complementary to the geometry of the template are necessary for effective shape selective binding. Michaelis-Menten kinetics implies that only the catalysts with imidazole moieties act as enzyme-like catalysts and imidazole is the key catalytic function of the enzyme mimics.     

A comprehensive analysis of the computed tautomer fractions of the imidazole ring of histidines in Loligo vulgaris

A recently introduced electrostatic-based method to determine the pKa values of ionizable residues and fractions of ionized and tautomeric forms of histidine (His) and acid residues in proteins, at a given fixed pH, is applied here to the analysis of a His-rich protein, namely Loligo vulgaris (pdb id 1E1A), a 314-residue all-β protein. The average tautomeric fractions for the imidazole ring of each of the six histidines in the sequence were computed using an approach that includes, but is not limited to, molecular dynamic simulations coupled with calculations of the ionization states for all 94 ionizable residues of protein 1E1A in water at pH 6.5 and 300 K. The electrostatic-calculated tautomeric fractions of the imidazole ring of His were compared with predictions obtained from an existent NMR-based methodology. Our results indicate that: (i) the averaged electrostatic-based tautomeric predictions for the imidazole ring of all histidines of Loligo vulgaris are dominated by the N^ε2-H rather than the N^δ1-H form, although such preferences from the NMR-based methodology are not so well defined; (ii) the computed average absolute difference between the electrostatic- and the NMR-based tautomeric predictions among all six histidines vary among 0% to 17%; (iii) for the His showing the largest fraction of the neutral form (81%), the absolute difference between the NMR- and electrostatic-based computed tautomeric predictions is only 3%; and (iv) the tautomeric predictions for the imidazole ring of His computed with the NMR-based methodology are stable within a certain, well-defined, range of variations of a tautomer-related parameter.     

Characterization of midgut trypsin-like enzymes and three trypsinogen cDNAs from the lesser grain borer, Rhyzopertha dominica (Coleoptera: Bostrichidae)

Protein digestion in the lesser grain borer, Rhyzopertha dominica (F.) (Coleoptera: Bostrichidae), results from the action of a complex of serine proteinases present in the midgut. In this study we partially characterized trypsin-like enzyme activity against N-α-benzoyl- -arginine p-nitroanilide (BApNA) in midgut preparations and cloned and sequenced three cDNAs for trypsinogen-like proteins. BApNAase activity in R. dominica midgut was significantly reduced by serine proteinase inhibitors and specific inhibitors of trypsin, whereas BApNAase activity was not sensitive to specific inhibitors of chymotrypsin or aspartic proteinases. However, trans-epoxysuccinyl- -leucylamido-(4-guanidino) butane (E-64) inhibited BApNAase activity by about 30%. BApNAase was most active in a broad pH range from about pH 7 to 9.5. The gut of R. dominica is a tubular tract approximately 2.5 mm in length. BApNAase activity was primarily located in the midgut region with about 1.5-fold more BApNAase activity in the anterior region compared to that in the posterior region. Proteinases with apparent molecular masses of 23–24 kDa that were visualized on casein zymograms following electrophoresis were inhibited by TLCK.Three cDNAs for trypsinogen-like proteins were cloned and sequenced from mRNA of R. dominica midgut. The full cDNA sequences consisted of open reading frames encoding 249, 293, and 255 amino acid residues for RdoT1, RdoT2, and RdoT3, respectively. cDNAs RdoT1, RdoT2, and RdoT3 shared 77–81% sequence identity. The three encoded trypsinogens shared 54–62% identity in their amino acid sequences and had 16–18 residues of signal peptides and 12–15 residues of activation peptides. The three predicted mature trypsin-like enzymes had molecular masses of 23.1, 28, and 23.8 kDa for RdoT1, RdoT2, and RdoT3, respectively. Typical features of these trypsin-like enzymes included the conserved N-terminal residues IVGG^62–65, the catalytic amino acid triad of serine proteinase active sites (His¹⁰⁹, Asp¹⁵⁶, Ser²⁵⁷), three pairs of conserved cysteine residues for disulfide bridges, and the three residues (Asp²⁵¹, Gly²⁷⁴, Gly²⁸⁴) that determine specificity in trypsin-like enzymes. In addition, RdoT2 has both a PEST-like sequence at the C-terminus and a free Cys¹⁵⁸ near the active site, suggesting instability of this enzyme and/or sensitivity to thiol reagents. The sequences have been deposited in GenBank database (accession numbers AF130840 for RdoT1, AF130841 for RdoT2, and AF130842 for RdoT3).     

Protein selectivity in immobilized metal affinity chromatography based on the surface accessibility of aspartic and glutamic acid residues

The interaction of different species variants of cytochrome c and myoglobin, as well as hen egg white lysozyme, with the hard Lewis metal ions Al³⁺, Ca²⁺, Fe³⁺, and Yb³⁺ and the borderline metal ion Cu²⁺, immobilized to iminodiacetic acid (IDA)-Sepharose CL-4B, has been investigated over the rangepH 5.5–8.0. With appropriately chosen buffer and metal ion conditions, these proteins can be bound to the immobilized M ⁿ +-IDA adsorbents via negatively charged amino acid residues accessible on the protein surface. For example, tuna heart cytochrome c, which lacks surface-accessible histidine residues, readily bound to the Fe³⁺-IDA adsorbent, while the other proteins also showed affinity toward immobilized Fe³⁺-IDA adsorbents when buffers containing 30 mM of imidazole were used. These studies document that protein selectivity can be achieved with hard-metalion immobilized metal ion affinity chromatography (IMAC) systems through the interaction of surfaceexposed aspartic and glutamic acid residues on the protein with the immobilized M ⁿ +-IDA complex. These investigations have also documented that the so-called soft or borderline immobilized metal ions such as the Cu²⁺-IDA adsorbent can also interact with surface-accessible aspartic and glutamic acid residues in a protein-dependent manner. A relationship is evident between the number and extent of clustering of the surfaceaccessible aspartic and glutamic acid residues and protein selectivity with these IMAC systems. The use of elution buffers which contain organic compound modifiers which replicate the carboxyl group moieties of these amino acids on the surface of proteins is also described.Abbreviations IDA iminodiacetic acid - IDA-Mⁿ⁺ iminodiacetic acid chelated to metal ion - IMAC immobilized metal affinity chromatography - DHCC dog heart cytochrome c - HHCC horse heart cytochrome c, THCC, tuna heart cytochrome c - HMYO horse skeletal muscle myoglobin - SMYO sheep skeletal muscle myoglobin - HEWL hen egg white lysozyme     

Histidine hydrogen-deuterium exchange mass spectrometry for probing the microenvironment of histidine residues in dihydrofolate reductase

Background

Histidine Hydrogen-Deuterium Exchange Mass Spectrometry (His-HDX-MS) determines the HDX rates at the imidazole C₂-hydrogen of histidine residues. This method provides not only the HDX rates but also the pK _a values of histidine imidazole rings. His-HDX-MS was used to probe the microenvironment of histidine residues of E. coli dihydrofolate reductase (DHFR), an enzyme proposed to undergo multiple conformational changes during catalysis.

Methodology/Principal Findings

Using His-HDX-MS, the pK _a values and the half-lives (t _1/2) of HDX reactions of five histidine residues of apo-DHFR, DHFR in complex with methotrexate (DHFR-MTX), DHFR in complex with MTX and NADPH (DHFR-MTX-NADPH), and DHFR in complex with folate and NADP⁺ (DHFR-folate-NADP⁺) were determined. The results showed that the two parameters (pK _a and t _1/2) are sensitive to the changes of the microenvironment around the histidine residues. Although four of the five histidine residues are located far from the active site, ligand binding affected their pK _a, t _1/2 or both. This is consistent with previous observations of ligand binding-induced distal conformational changes on DHFR. Most of the observed pK _a and t _1/2 changes could be rationalized using the X-ray structures of apo-DHFR, DHFR-MTX-NADPH, and DHFR-folate-NADP⁺. The availability of the neutron diffraction structure of DHFR-MTX enabled us to compare the protonation states of histidine imidazole rings.

Conclusions/Significance

Our results demonstrate the usefulness of His-HDX-MS in probing the microenvironments of histidine residues within proteins.     

Potentiometric titration curves of oxidized and reduced horse heart cytochrome c

Potentiometric titration curves of oxidized and reduced horse heart cytochrome c in 0.15M KCl at 20°C have been obtained by timed titration (0.125–0.500 μmol/sec) from the isoionic points (pH 10.2–10.4) to pH 3 and back to the isoionic point. Computer-assisted (PROPHET) data acquisition and blank corrections give curves with good precision with a maximum standard deviation of 0.3 groups for an average error of 1%. The potentiometric titration curve of reduced cytochrome c is reversible within the precision of the method and for the pH range studied. The potentiometric curves for oxidized cytochrome c titrated upscale (pH 3–10) and downscale (pH 10–3) are not reversible. However, they show the same ionization behavior after the initial downscale titration. This is probably the result of a conformational change. Comparison of the data herein reported with the titration curves of oxidized cytochrome c already published by others indicates good agreement on the basis of a normalization of the concentration of protein or on the basis of 25 titrable groups between the acid end point and the isoionic pH. Titration of the 2 μmol imidazole in the upscale or downscale direction gives the correct analytical concentration and pK′ after correction for the solvent titration. Titration of reduced cytochrome c in the presence and absence of an additional equivalent of imidazole gave a difference titration curve, which indicates that a group on the protein shifts from pK′ 5.8 to pK′ 5.3 in the presence of imidazole. The pK′ of imidazole, in the presence of the protein, remains at a nearly normal value of 7.34.