Titration Curves of Amino Acids

Drop a spoonful of glycine into water. Add acid until the solution is syrupy with protons. Then titrate slowly with sodium hydroxide. A pH electrode records every hundredth of a milliliter. The resulting curve is not a smooth ramp. It is a staircase punctuated by two flat plateaus. Two inflection points thread through it. Every kink in that staircase is a direct readout. It shows a specific proton leaving a specific atom. Learning to read that staircase is a key skill. It is one of the most information-dense skills in the acid-base chemistry of biomolecules. A single experiment, done with nothing more than a buret and a pH meter, reveals every pKapK_a an amino acid owns. This matters far beyond a teaching lab. The shape of a titration curve tells you how many ionizable groups a molecule carries. It also shows what those groups are chemically. Critically, it reveals the net charge of the molecule at any pH the cell might present it with. Enzyme active sites, ion-exchange chromatography, protein solubility, isoelectric focusing, and the folding of proteins into charge-balanced structures all trace back to exactly this bookkeeping. A biochemist can look at a titration curve. They can immediately name the groups being deprotonated — carboxyl, ammonium, imidazole, thiol, phenol. Doing so means reading the molecule's electrostatic personality directly off the graph. The big picture is a marriage of two ideas already in hand. The Acid-Base Chemistry of the Alpha-Amino and Alpha-Carboxyl Groups gives every amino acid at least two ionizable sites. The Acid-Base Chemistry of Ionizable Side Chains adds a third (or more). This applies to the seven amino acids with ionizable R groups. A titration curve is simply the Henderson-Hasselbalch equation applied group by group. It is stitched together along a single pH axis. For a monoprotic system (a simple weak acid) the curve has one buffering region. It also has one equivalence point. For glycine it has two of each. For a trifunctional amino acid like glutamate or histidine, it has three. This article builds the machinery to interpret any such curve. It shows how to identify the number of equivalence points and half-equivalence points. It shows how to assign each inflection to a specific functional group. This assignment uses each group's characteristic pKapK_a. It also shows how to compute the buffering capacity and titration volumes quantitatively. Finally, it shows how the whole curve sets up a key concept. That concept is the isoelectric point, developed fully elsewhere.

Key takeaways
  • A titration curve plots pH versus equivalents of strong base (or acid) added, and it is built from one Henderson-Hasselbalch segment per ionizable group.
  • Each buffering region is centered on a pKapK_a; each equivalence point falls exactly between two buffering regions, at the pH of the pure intermediate species.
  • Monoprotic amino acids in the strict sense do not exist among the 20 standard residues — even glycine has two groups — but any single weak-acid/base pair titrates with one plateau and one equivalence point.
  • Diprotic amino acids (no ionizable side chain) show two plateaus, two equivalence points, and the isoelectric point sits at the volume/pH exactly halfway between the two pKapK_a values flanking the neutral zwitterion.
  • Triprotic amino acids (ionizable side chain) show three plateaus and three equivalence points; the isoelectric point calculation depends on which two pKapK_a values bracket the zero-net-charge species.
  • Buffering capacity is maximal at each pKapK_a (pH = pKapK_a) and vanishes at each equivalence point, where the slope dpH/dVd\text{pH}/dV is steepest.
  • Key idea: a titration curve is a chart of protons leaving specific atoms in order of increasing pKapK_a, and every flat region and every steep jump has a one-to-one chemical meaning.

Where this lives today

Ion-exchange chromatography, a workhorse of protein and peptide purification, is built entirely on titration-curve logic: a cation-exchange resin retains proteins whose net charge (computed exactly as in the derivations above, summed over all ionizable side chains) is positive at the running buffer's pH, and releases them as the pH or salt gradient shifts that net charge through zero and negative. Isoelectric focusing, used routinely in 2D gel electrophoresis for proteomics, physically separates proteins by migrating them through a pH gradient until each stops at its own isoelectric point — the direct multi-group extension of the pI concept introduced here. In pharmaceutical formulation, the ionization state of an amino acid side chain or a drug's amine/carboxyl groups (governed by the same Henderson-Hasselbalch bookkeeping) determines aqueous solubility, membrane permeability, and even the counter-ion selected for a salt form — a huge fraction of oral drug formulation chemistry is precisely this titration-curve reasoning applied to a small molecule instead of a free amino acid. In enzymology, the pH-dependence of catalytic rate for enzymes like chymotrypsin or lysozyme is interpreted using exactly this framework: a bell-shaped rate-versus-pH profile reflects two ionizable active-site residues (often a general acid and a general base) each behaving as an independent Henderson-Hasselbalch titration, and the inflection points of that rate profile are read off just like inflection points on an amino acid titration curve. Finally, buffer selection in any biochemistry laboratory — choosing Tris, phosphate, HEPES, or acetate to hold a reaction at a target pH — is a direct application of the buffering-capacity result β=CTln⁡10 f(1−f)\beta = C_T\ln10\,f(1-f), maximized by choosing a buffer whose pKapK_a sits within one unit of the desired pH.

:::quiz[At the exact half-equivalence volume of a titration step, what is true of the ratio [A−]/[HA][\text{A}^-]/[\text{HA}] for that group?]

  • It equals 1, so pH = pKapK_a ✓
  • It equals 0
  • It approaches infinity
  • It equals pKapK_a numerically

Ask the tutor: How do I combine every side-chain and terminal pKapK_a in a full titration curve to calculate the isoelectric point of a trifunctional amino acid like histidine or lysine?

The physical picture — reading the staircase

Consider glycine, H3−CH2−COOH\text{H}_3\overset{+}{\text{N}}-\text{CH}_2-\text{COOH}, dissolved as the fully protonated cation by adding excess strong acid before the titration begins. As NaOH is added, the strongest acid in solution reacts first: the α\alpha-carboxyl proton, with pKa1≈2.34pK_{a1} \approx 2.34, leaves before the ammonium proton even begins to ionize, because the two pKapK_a values are separated by more than four units and the equilibria do not overlap significantly. The curve therefore shows a first buffering region centered near pH 2.34, a first equivalence point at pH ≈5.97\approx 5.97 (the average of the two pKapK_a values, where the zwitterion dominates), a second buffering region centered near pH 9.60 as the ammonium group ionizes, and a second equivalence point at high pH where the fully deprotonated anion dominates. The chemical logic behind the ordering is straightforward Brønsted acidity: a proton on an α\alpha-carboxyl group sitting adjacent to a protonated ammonium cation is strongly acidified by the inductive electron withdrawal of that positive charge, giving pKa1pK_{a1} values around 1.8–2.4 for all 20 amino acids — far lower than acetic acid's 4.76, because acetic acid has no adjacent cationic group. The α\alpha-ammonium group, by contrast, is a comparatively weak acid, with pKa2pK_{a2} typically 9.0–10.8, similar to other primary ammonium ions. Side-chain groups fall in between or overlap these ranges depending on their chemistry: the imidazolium of histidine at pKR≈6.0pK_R \approx 6.0, the thiol of cysteine at pKR≈8.3pK_R \approx 8.3, the phenol of tyrosine at pKR≈10.5pK_R \approx 10.5, the carboxyls of aspartate and glutamate at pKR≈3.65pK_R \approx 3.65 and 4.254.25, and the guanidinium of arginine at pKR≈12.5pK_R \approx 12.5.

Figure
Titration curve of glycine: two buffering plateaus, two equivalence points
OH⁻ equivalents added pH 0 0.5 1.0 1.5 2.0 2 6 9.6 12 pKa1 = 2.34 (–COOH) 1st equiv. pt, pH 5.97 (zwitterion) pKa2 = 9.60 (–NH3+) 2nd equiv. pt, pH ≈ 11.3 (anion)

The half-equivalence points are the most chemically informative landmarks on the curve. At exactly half of one equivalent of base added — 0.5 equivalents for the first group, 1.5 equivalents for the second — exactly half of that group's population has been converted from acid form to base form. By the Henderson-Hasselbalch relation this means [A−]=[HA][\text{A}^-] = [\text{HA}] for that group, so log⁡([A−]/[HA])=0\log([\text{A}^-]/[\text{HA}]) = 0 and pH =pKa= pK_a exactly. This is why a titration curve is the standard experimental method for measuring pKapK_a values in the first place: locate the half-equivalence volume, read the pH.

Figure
Triprotic curve for glutamate: three plateaus, three equivalence points, pI between pKa1 and pKR
OH⁻ equivalents added pH 0 1.0 2.0 3.0 2.19 4.25 9.67 pKa1 2.19 (α-COOH) pI = 3.22 = (2.19+4.25)/2 pKR 4.25 (γ-COOH) 2nd equiv. pt (net –1 species) pKa2 9.67 (α-NH3+) 3rd equiv. pt (net –2)

The core ideas — the governing equations

The entire curve rests on one relationship, applied independently to each ionizable group ii with its own pKa,ipK_{a,i}:

pH=pKa,i+log⁡10 ⁣([Ai−][HAi])\text{pH} = pK_{a,i} + \log_{10}\!\left(\frac{[\text{A}_i^{-}]}{[\text{HA}_i]}\right)

where:

  • pKa,ipK_{a,i} is the acid dissociation constant of group ii (e.g. α\alpha-COOH, α\alpha-NH₃⁺, or an R-group),
  • [Ai−][\text{A}_i^-] and [HAi][\text{HA}_i] are the deprotonated and protonated concentrations of that specific group,
  • the equation holds group-by-group even when multiple groups are present, provided their pKapK_a values are separated by roughly two units or more so that the equilibria do not appreciably overlap.

Reading: whenever the ratio of base form to acid form for a given group equals 1, the log term vanishes and the pH equals that group's pKapK_a — this is the flattest, most buffered point in that group's titration segment, and it occurs at the half-equivalence volume for that group. The fraction of base-form species αi\alpha_i (degree of deprotonation) for group ii follows directly:

αi=[Ai−][Ai−]+[HAi]=10 pH−pKa,i1+10 pH−pKa,i\alpha_i = \frac{[\text{A}_i^-]}{[\text{A}_i^-]+[\text{HA}_i]} = \frac{10^{\,\text{pH}-pK_{a,i}}}{1+10^{\,\text{pH}-pK_{a,i}}}

where:

  • αi\alpha_i ranges from 0 (fully protonated) to 1 (fully deprotonated),
  • this is a sigmoidal (logistic) function of pH centered at pKa,ipK_{a,i}.

Reading: the sum of all group fractions, weighted by charge, gives the net charge of the molecule at any pH — this sum is exactly what generates the Isoelectric Point when set to zero. For a molecule with nn ionizable groups, the equivalence-point pH between the kk-th and (k+1)(k+1)-th group (when their pKapK_a values are well separated) is the simple average:

pHeq,k=pKa,k+pKa,k+12\text{pH}_{\text{eq},k} = \frac{pK_{a,k} + pK_{a,k+1}}{2}

where:

  • this is only exact when the flanking pKapK_a values differ by roughly ≥ 3 units, so that at that pH each group is either fully deprotonated or fully protonated to a good approximation,
  • for closely spaced pKapK_a values (e.g., aspartate's two carboxyls at 1.88 and 3.65 are separated by less than 2 units), corrections apply, though the amino-acid side chains discussed here are separated widely enough for the simple average to be an excellent approximation.
Derivation ALocating the half-equivalence pH from the Henderson-Hasselbalch equation

Goal: show algebraically why the pH at the half-equivalence volume of any titration step equals that group's pKapK_a, and derive the buffering-capacity behavior around that point.

Step 1. Take a single ionizable group HA with total concentration CT=[HA]+[A−]C_T = [\text{HA}] + [\text{A}^-], fixed because no group is created or destroyed during titration — only converted from HA to A⁻ by added strong base. Let ff be the fraction of the equivalents required to fully deprotonate this group that have been added so far (0≤f≤10 \le f \le 1 for one titration step). Since each mole of added OH⁻ converts one mole of HA to A⁻ (until exhausted),

[A−]=f CT,[HA]=(1−f) CT[\text{A}^-] = f\,C_T, \qquad [\text{HA}] = (1-f)\,C_T

Step 2. Substitute into the Henderson-Hasselbalch equation:

pH=pKa+log⁡10 ⁣(f CT(1−f) CT)=pKa+log⁡10 ⁣(f1−f)\text{pH} = pK_a + \log_{10}\!\left(\frac{f\,C_T}{(1-f)\,C_T}\right) = pK_a + \log_{10}\!\left(\frac{f}{1-f}\right)

The total concentration CTC_T cancels exactly, which is the algebraic reason the titration curve's shape is independent of the amino acid's concentration (dilute solutions give the same curve shape, only affecting total buffer capacity in absolute terms).

Step 3. Set f=1/2f = 1/2, the half-equivalence point:

pH=pKa+log⁡10 ⁣(1/21/2)=pKa+log⁡10(1)=pKa+0=pKa\text{pH} = pK_a + \log_{10}\!\left(\frac{1/2}{1/2}\right) = pK_a + \log_{10}(1) = pK_a + 0 = pK_a

This confirms that at exactly half the equivalents needed to titrate one group, the measured pH is that group's pKapK_a — the operational basis for reading pKapK_a values directly off the curve.

Step 4. Buffering capacity β\beta is defined as the moles of strong base absorbed per unit pH change, β=dCb/d(pH)\beta = dC_b/d(\text{pH}), where Cb=f CTC_b = f\,C_T is the concentration of base added. Differentiate pH with respect to ff from Step 2:

d(pH)df=ddf[log⁡10 ⁣(f1−f)]=1ln⁡10(1f+11−f)=1ln⁡10⋅1f(1−f)\frac{d(\text{pH})}{df} = \frac{d}{df}\left[\log_{10}\!\left(\frac{f}{1-f}\right)\right] = \frac{1}{\ln 10}\left(\frac{1}{f}+\frac{1}{1-f}\right) = \frac{1}{\ln 10}\cdot\frac{1}{f(1-f)}

Inverting and multiplying by CTC_T (since Cb=fCTC_b = fC_T, dCb=CT dfdC_b = C_T\,df):

β=CT dfd(pH)=CT ln⁡10⋅f(1−f)\beta = C_T\,\frac{df}{d(\text{pH})} = C_T\,\ln 10 \cdot f(1-f)

Step 5. This is a parabola in ff, maximized at f=1/2f = 1/2 (giving βmax⁡=CTln⁡10/4≈0.576 CT\beta_{\max} = C_T \ln10/4 \approx 0.576\,C_T) and equal to zero at f=0f=0 and f=1f=1 — the endpoints of each titration step. This is the mechanistic reason plateaus (high buffering, small pH change per drop of titrant) occur at pKapK_a values and steep jumps (near-zero buffering) occur at equivalence points, where ff for one group is 1 and for the next group is still 0.

Go deeper · Derivation B · 2 minNet charge and the equivalence-point spacing for a diprotic amino acid
Derivation BNet charge and the equivalence-point spacing for a diprotic amino acid

Goal: derive the net charge as a function of pH for a simple diprotic amino acid (glycine-type, no ionizable side chain) and show why the equivalence point separating the cation and zwitterion sits exactly at (pKa1+pKa2)/2(pK_{a1}+pK_{a2})/2 when the two pKapK_a values are well separated.

Step 1. Glycine has three possible protonation states: the fully protonated cation H3N+-CH2-COOH\text{H}_3\text{N}^+\text{-CH}_2\text{-COOH} (net charge +1+1), the zwitterion H3N+-CH2-COO−\text{H}_3\text{N}^+\text{-CH}_2\text{-COO}^- (net charge 00), and the anion H2N-CH2-COO−\text{H}_2\text{N-CH}_2\text{-COO}^- (net charge −1-1). Using the degree-of-deprotonation formula for each group, let α1\alpha_1 be the fraction of the carboxyl group deprotonated (COOH → COO⁻) with pKa1=2.34pK_{a1}=2.34, and α2\alpha_2 the fraction of the ammonium group deprotonated (NH₃⁺ → NH₂) with pKa2=9.60pK_{a2}=9.60.

Step 2. The net average charge ZZ is the charge of the fully protonated species plus the charge changes contributed by each deprotonation:

Z(pH)=(+1)−α1(pH)−α2(pH)Z(\text{pH}) = (+1) - \alpha_1(\text{pH}) - \alpha_2(\text{pH})

because losing the carboxyl proton changes charge by −1-1 (from the reference +1+1 cation) and losing the ammonium proton changes charge by another −1-1.

Step 3. At low pH (pH ≪ 2.34), α1→0\alpha_1 \to 0 and α2→0\alpha_2 \to 0, giving Z→+1Z \to +1: the fully protonated cation, consistent with the picture. At high pH (pH ≫ 9.60), both α1,α2→1\alpha_1, \alpha_2 \to 1, giving Z→−1Z \to -1: the fully deprotonated anion.

Step 4. At the first equivalence point, exactly one equivalent of base per mole of amino acid has been added, converting the carboxyl group essentially fully (since pKa1pK_{a1} and pKa2pK_{a2} differ by more than 7 units, at any pH near 5–7 the ammonium group is still >99.9% protonated). So α1≈1\alpha_1 \approx 1 and α2≈0\alpha_2 \approx 0 at this point, and

Z≈(+1)−1−0=0Z \approx (+1) - 1 - 0 = 0

confirming the solution is dominated by the neutral zwitterion. Because α1\alpha_1 transitions from 0 to 1 symmetrically around pKa1pK_{a1} and α2\alpha_2 is still negligible there, the pH at which Z=0Z=0 is essentially pKa1pK_{a1} itself is not quite right — rather, the isoelectric species (zwitterion) is the dominant form over a range, and the pH at which the population-averaged charge is exactly zero, which is the true isoelectric point, requires balancing the small residual protonation of the carboxyl (favoring ++, since pKa1pK_{a1} pulls pH up slightly below full dissociation) against the small residual deprotonation of the ammonium (favoring −-). Setting Z=0Z=0 explicitly:

1−10pH−pKa11+10pH−pKa1−10pH−pKa21+10pH−pKa2=01 - \frac{10^{\text{pH}-pK_{a1}}}{1+10^{\text{pH}-pK_{a1}}} - \frac{10^{\text{pH}-pK_{a2}}}{1+10^{\text{pH}-pK_{a2}}} = 0

Step 5. For pH between pKa1pK_{a1} and pKa2pK_{a2} and well separated from both (as is the case here, since pKa2−pKa1=7.26pK_{a2}-pK_{a1} = 7.26), the first fraction is close to 1 and the second is close to 0; write α1=1−ϵ1\alpha_1 = 1-\epsilon_1 and α2=ϵ2\alpha_2 = \epsilon_2 where ϵ1,ϵ2≪1\epsilon_1, \epsilon_2 \ll 1. Then Z≈ϵ1−ϵ2Z \approx \epsilon_1 - \epsilon_2. Using ϵ1=1/(1+10pH−pKa1)=10pKa1−pH\epsilon_1 = 1/(1+10^{\text{pH}-pK_{a1}}) = 10^{pK_{a1}-\text{pH}} for pH well above pKa1pK_{a1}, and ϵ2=10pH−pKa2\epsilon_2 = 10^{\text{pH}-pK_{a2}} for pH well below pKa2pK_{a2}, setting ϵ1=ϵ2\epsilon_1=\epsilon_2 gives

10pKa1−pH=10pH−pKa2  ⟹  pKa1−pH=pH−pKa2  ⟹  pH=pKa1+pKa2210^{pK_{a1}-\text{pH}} = 10^{\text{pH}-pK_{a2}} \implies pK_{a1}-\text{pH} = \text{pH}-pK_{a2} \implies \text{pH} = \frac{pK_{a1}+pK_{a2}}{2}

This confirms the isoelectric/equivalence-point pH is the simple arithmetic mean of the two flanking pKapK_a values, and shows explicitly why — it is the pH at which the exponentially small leakage of charge from each side balances exactly.

Worked examples

Example 1 — Foundational: locating every landmark on glycine's titration curve. Glycine has pKa1=2.34pK_{a1}=2.34 (α-COOH) and pKa2=9.60pK_{a2}=9.60 (α-NH₃⁺). Suppose we titrate 10.00 mL of a 0.100 M solution of the fully protonated cationic form (glycine hydrochloride) with 0.100 M NaOH. Total moles of glycine =0.010 L×0.100 M=1.00×10−3= 0.010\text{ L} \times 0.100\text{ M} = 1.00\times10^{-3} mol. One full equivalent of NaOH per ionizable group therefore requires 1.00×10−31.00\times10^{-3} mol NaOH, i.e., 10.0010.00 mL of 0.100 M NaOH per group.

  • At V=0V=0 mL (0 equivalents): pure cation, pH set by pKa1pK_{a1} and initial concentration (compute via pH=12(pKa1−log⁡C0)\text{pH} = \tfrac12(pK_{a1}-\log C_0) for a weak acid approximation using C0=0.100C_0=0.100 M): pH≈12(2.34+1.00)=1.67\text{pH} \approx \tfrac12(2.34+1.00)=1.67. (Reported experimental starting pH for 0.1 M glycine·HCl is close to 1.7–1.8, consistent.)
  • At V=5.00V=5.00 mL (0.5 equivalents, half of the first group): pH=pKa1=2.34\text{pH} = pK_{a1} = 2.34 exactly, by Step 3 of the first derivation.
  • At V=10.00V=10.00 mL (1.0 equivalent): first equivalence point. pH=(pKa1+pKa2)/2=(2.34+9.60)/2=5.97\text{pH} = (pK_{a1}+pK_{a2})/2 = (2.34+9.60)/2 = 5.97 — the isoelectric point, dominated by the neutral zwitterion.
  • At V=15.00V=15.00 mL (1.5 equivalents, half of the second group): pH=pKa2=9.60\text{pH} = pK_{a2} = 9.60 exactly.
  • At V=20.00V=20.00 mL (2.0 equivalents): second equivalence point, pH determined by hydrolysis of the fully deprotonated anion (a weak base); estimate via pH=14−12(pKb−log⁡C)\text{pH} = 14 - \tfrac12(pK_b-\log C) where pKb=14−pKa2=4.40pK_b = 14-pK_{a2}=4.40 and C≈0.0333C \approx 0.0333 M (diluted, since total volume is now 30 mL): pOH=12(4.40+1.48)=2.94\text{pOH}=\tfrac12(4.40+1.48)=2.94, so pH≈11.06\text{pH}\approx 11.06.

Every one of these five landmark points is directly measurable with a pH meter and buret, and together they fully specify both pKapK_a values and the isoelectric point without needing any other technique.

Example 2 — Applied: titrating aspartate and identifying an unknown amino acid from its curve. Suppose a biochemistry student is handed an unlabeled amino acid and titrates 10.0 mL of a 0.0500 M solution of its fully protonated form with 0.0500 M NaOH, recording three inflection points at cumulative volumes of 10.0 mL, 20.0 mL, and 30.0 mL, with half-equivalence pH readings of 1.88, 3.65, and 9.60. Three buffering regions and three equivalence points immediately signal a trifunctional (triprotic) amino acid — one with an ionizable side chain in addition to the α-carboxyl and α-amino groups. The half-equivalence pH values read directly as the three pKapK_a values: pKa1=1.88pK_{a1}=1.88 (α-COOH, unusually low, since standard α\alpha-COOH values are ~2.0–2.4, but suppressed further here by proximity to a second nearby carboxylate), pKR=3.65pK_{R}=3.65 (side-chain carboxyl), and pKa2=9.60pK_{a2}=9.60 (α-NH₃⁺). Matching against standard reference values (pKa1=1.88pK_{a1}=1.88, pKR=3.65pK_R=3.65, pKa2=9.60pK_{a2}=9.60 for aspartate) identifies the unknown unambiguously as aspartate. The equivalence points in this three-group system occur at:

pHeq,1=pKa1+pKR2=1.88+3.652=2.77 (net charge 0→cation to neutral zwitterion boundary)\text{pH}_{\text{eq},1} = \frac{pK_{a1}+pK_R}{2} = \frac{1.88+3.65}{2} = 2.77 \text{ (net charge } 0 \to \text{cation to neutral zwitterion boundary)}
pHeq,2=pKR+pKa22=3.65+9.602=6.63 — this is the isoelectric point, pI\text{pH}_{\text{eq},2} = \frac{pK_R+pK_{a2}}{2} = \frac{3.65+9.60}{2} = 6.63 \text{ — this is the isoelectric point, pI}
pHeq,3: beyond the third equivalence point, dominated by fully deprotonated trianion\text{pH}_{\text{eq},3}: \text{ beyond the third equivalence point, dominated by fully deprotonated trianion}

Note carefully which two pKapK_a values flank the neutral (net charge zero) species: for aspartate, the zero-charge zwitterion has the α-ammonium protonated (+1) and both carboxyls deprotonated (−1, −1), giving net charge +1−1−1=−1+1-1-1=-1 at pH between pKa1pK_{a1} and pKRpK_R — that is not neutral. Working through the charge states explicitly: below pKa1pK_{a1}, charge is +1+1 (both COOH protonated, NH₃⁺); between pKa1pK_{a1} and pKRpK_R, charge is 00 (α-COOH deprotonated, side-chain COOH still protonated, NH₃⁺ intact); between pKRpK_R and pKa2pK_{a2}, charge is −1-1 (both COOH deprotonated, NH₃⁺ intact); above pKa2pK_{a2}, charge is −2-2. The neutral (net charge zero) species dominates between pKa1pK_{a1} and pKRpK_R, so the true isoelectric point is

pI=pKa1+pKR2=1.88+3.652=2.77\text{pI} = \frac{pK_{a1}+pK_R}{2} = \frac{1.88+3.65}{2} = 2.77

This example is worked fully because it demonstrates a common trap: the isoelectric point is not automatically the average of the middle two pKapK_a values — it is the average of whichever two pKapK_a values flank the neutral species, found by explicitly tracking charge state group by group (this logic is developed fully as the general method in the treatment of the Isoelectric Point).

Example 3 — Exam/competition level: back-calculating pKRpK_R from a partial titration curve of histidine. A 20.0 mL sample of 0.0800 M histidine, prepared in its fully protonated triprotonated cationic form (α-COOH protonated, imidazolium protonated, α-NH₃⁺ protonated; net charge +2+2), is titrated with 0.1600 M NaOH. The student records the pH at several volumes: at V=5.00V=5.00 mL, pH =1.82=1.82; the first equivalence point occurs at V=10.0V=10.0 mL; at V=17.5V=17.5 mL, pH =6.00= 6.00; the second equivalence point occurs at V=20.0V=20.0 mL; a further equivalence point occurs at V=30.0V = 30.0 mL. First, find total moles of histidine: 0.0200 L×0.0800 M=1.60×10−30.0200\text{ L}\times0.0800\text{ M} = 1.60\times10^{-3} mol. Moles of NaOH needed for one full equivalent: 1.60×10−31.60\times10^{-3} mol, requiring 1.60×10−3 mol/0.1600 M=10.01.60\times10^{-3}\text{ mol}/0.1600\text{ M} = 10.0 mL — consistent with the first equivalence point at V=10.0V=10.0 mL, confirming the stoichiometry and group count (three groups titrated total, since equivalence points appear at 10.0, 20.0, and 30.0 mL, each spaced by 10.0 mL). The half-equivalence volume for the first group is V=5.00V=5.00 mL (half of 10.0 mL), where pH =1.82=1.82, so

pKa1=1.82(α-COOH)pK_{a1} = 1.82 \quad (\alpha\text{-COOH})

This is the group that titrates first because the α-carboxyl is the most acidic (lowest pKapK_a) of the three histidine groups. The half-equivalence volume for the second group falls at the midpoint between the first equivalence point (10.0 mL) and the second (20.0 mL), i.e., at V=10.0+5.00=15.0V = 10.0+5.00=15.0 mL — but the problem gives a reading at V=17.5V=17.5 mL instead, which is not the half-equivalence point for the second group; it is three-quarters of the way through the second titration step. To extract pKRpK_R correctly, use the general Henderson-Hasselbalch form for that step. Within the second titration step (from 10.0 to 20.0 mL), define ff as the fraction of that step's equivalents delivered:

f=V−10.020.0−10.0=17.5−10.010.0=0.75f = \frac{V-10.0}{20.0-10.0} = \frac{17.5-10.0}{10.0} = 0.75

Then from Step 2 of the first derivation,

pH=pKR+log⁡10 ⁣(f1−f)=pKR+log⁡10 ⁣(0.750.25)=pKR+log⁡10(3)\text{pH} = pK_R + \log_{10}\!\left(\frac{f}{1-f}\right) = pK_R + \log_{10}\!\left(\frac{0.75}{0.25}\right) = pK_R + \log_{10}(3)
6.00=pKR+0.477  ⟹  pKR=6.00−0.477=5.526.00 = pK_R + 0.477 \implies pK_R = 6.00-0.477 = 5.52

So the imidazole side chain of this particular histidine sample has pKR=5.52pK_R = 5.52 (a reasonable value; free histidine's imidazolium pKRpK_R is commonly cited near 6.0, and the modest downward shift here is chemically plausible depending on solution conditions such as ionic strength). By elimination and standard chemical assignment (the remaining group must be the α-ammonium, the last and least acidic group to titrate), pKa2pK_{a2} for the α-NH₃⁺ can be located at the third half-equivalence volume, V=25.0V=25.0 mL, and would need an additional pH reading to pin down numerically — illustrating that a complete triprotic curve requires three independent half-equivalence readings, one per group, and that equivalence-point volumes alone fix the group count and stoichiometry while half-equivalence pH values fix each pKapK_a individually.

Limiting cases, scope & common traps

CaseBehaviorWhy
pKapK_a values separated by ≥ 3 unitsEquivalence point = simple average of flanking pKapK_a'sEach group is essentially fully reacted before the next begins; negligible overlap in Step 5 of the second derivation
pKapK_a values separated by < 2 units (e.g. some diprotic dicarboxylic acids)Equivalence point deviates from the simple average; curve plateaus merge/blurOverlapping equilibria mean neither group is "done" before the next starts; requires solving simultaneous equilibria
Very dilute solution (CT→0C_T \to 0)Buffering capacity β=CTln⁡10 f(1−f)→0\beta = C_T\ln10\, f(1-f) \to 0 everywhereβ\beta scales linearly with total concentration; curve shape (pH vs. equivalents) is unchanged, but absolute buffering strength vanishes
f→0f \to 0 or f→1f \to 1 (near start/end of any step)β→0\beta \to 0, steepest pH change per dropParabolic β(f)\beta(f) derived above vanishes at the endpoints of each titration step
f=1/2f = 1/2 (half-equivalence)β\beta maximal, pH flattest, pH =pKa= pK_a exactlyDerived exactly in Step 3–4 of the first derivation
Side chain absent (Gly, Ala, Val, Leu, Ile, Ser, Thr, etc.)Two-step curve, two equivalence pointsOnly α-COOH and α-NH₃⁺ are ionizable
Side chain ionizable and acidic (Asp, Glu) or basic (Lys, Arg, His)Three-step curve, three equivalence pointsThird group adds a third buffering region and shifts pI away from the simple 2-group average
Cysteine and TyrosineAlso three-step, but the third pKRpK_R (8.3, 10.5) is close to pKa2pK_{a2} (~9–10.3)Groups can partially overlap, softening the distinction between the second and third plateaus

The chemist behind it

Lawrence Joseph Henderson and Karl Albert Hasselbalch, whose names are permanently welded to the equation underlying every titration curve, approached the buffer problem from opposite ends of physiology and chemistry. Henderson, working at Harvard in the early 1900s, was consumed by the question of how blood maintains a nearly constant pH despite constant metabolic acid production — he treated the blood bicarbonate system as a chemical equilibrium and, in 1908, wrote down the equilibrium relationship connecting [H+][\text{H}^+], [HCO3−][\text{HCO}_3^-], and dissolved CO2\text{CO}_2 essentially in the form we would now recognize, though without logarithms. Hasselbalch, a Danish physician and physiologist working on respiratory physiology in Copenhagen, took Henderson's relationship and in 1916 recast it using the newly popularized logarithmic pH scale (introduced by Sørensen in 1909), producing the compact logarithmic form used universally today. Neither man was working on amino acids specifically — both were physiologists chasing the chemistry of blood and respiration — but their equation turned out to be exactly the right tool for interpreting any weak acid/base equilibrium, amino acid titration curves included. The broader amino acid titration methodology was developed and systematized over subsequent decades by protein chemists measuring individual pKapK_a values for free amino acids and, later, for ionizable groups buried within folded proteins — work that underlies modern techniques from isoelectric focusing to NMR-based pKapK_a determination in enzyme active sites.

Origins · 1 minLawrence Joseph Henderson and how this was discovered
LJ
Lawrence Joseph Henderson
1878–1942 · Harvard physiological chemist who derived the equilibrium relationship for blood buffering, later extended to all weak acid/base systems

Henderson trained as a physician but spent his career as a physiological chemist, insisting that living systems obey ordinary physical chemistry — his 1913 book "The Fitness of the Environment" argued that the chemical properties of water, carbon, and carbonic acid make life on Earth possible. His 1908 equilibrium analysis of blood as a carbonic acid/bicarbonate buffer system, though originally expressed without logarithms, is the direct ancestor of the equation used to construct and interpret every amino acid titration curve.

History
- 1908 — Henderson derives the equilibrium expression for the blood bicarbonate buffer system.
- 1909 — Sørensen introduces the pH scale, providing the logarithmic framework buffer chemistry needed.
- 1916 — Hasselbalch recasts Henderson's relationship in logarithmic form, giving the modern Henderson-Hasselbalch equation.
- 1920s–1930s — Systematic potentiometric titration of free amino acids establishes reference pKapK_a tables still used today.
- Mid-20th century onward — Titration-curve methodology extended to individual ionizable residues within folded proteins, underlying modern enzyme mechanism studies.
Why does the α-carboxyl group of an amino acid have such a low pKapK_a (~2) compared to acetic acid's pKapK_a of 4.76?

The adjacent, positively charged α-ammonium group (-NH3+\text{-NH}_3^+) inductively withdraws electron density through the intervening α\alpha-carbon, stabilizing the developing negative charge on the carboxylate oxygen as the proton leaves. This electrostatic stabilization lowers the free energy of the deprotonated form relative to an isolated carboxylic acid like acetic acid, which has no such adjacent cation, shifting the equilibrium toward dissociation and lowering the pKapK_a by roughly 2.4 units. Once the ammonium group itself is deprotonated at higher pH, this inductive effect disappears, which is part of why the second carboxyl-like ionization (if present, as in Asp/Glu side chains) sits at a more "normal" pKapK_a of 3.65–4.25, much closer to acetic acid's value.

Can you determine the number of ionizable groups on an unknown amino acid just from equivalence-point volumes, without knowing any pKapK_a values in advance?

Yes. Each ionizable group consumes exactly one equivalent of titrant, so the number of equally spaced equivalence-point volumes directly counts the number of ionizable groups. If equivalence points appear at 10.0, 20.0, and 30.0 mL (equal 10.0 mL spacings), the molecule has three ionizable groups. The half-equivalence pH readings between consecutive equivalence points then supply the individual pKapK_a values, and comparing those values against standard reference tables (or side-chain chemistry) identifies which functional groups they belong to and, often, the specific amino acid.

Why is the buffering capacity exactly zero at an equivalence point rather than merely small?

At an equivalence point between two well-separated pKapK_a values, the group that finished titrating is essentially 100% in its base form (f=1f=1 for that step) and the next group has not yet begun accepting the base's deprotonating effect (f=0f=0 for the next step). The derived formula β=CTln⁡10 f(1−f)\beta = C_T\ln10\,f(1-f) gives exactly zero at f=0f=0 or f=1f=1. In practice it is never exactly zero because ff never reaches precisely 0 or 1 (there is always a vanishingly small residual population), but the buffering capacity becomes so small that a single drop of titrant produces a large, easily observed pH jump — which is precisely how equivalence points are detected experimentally.

Does dilution change the shape of the titration curve (pH vs. equivalents), or only its buffering strength?

Dilution changes only the buffering strength (the absolute amount of acid or base needed per unit pH change), not the shape of the curve when plotted as pH versus fraction of equivalents or versus ff. This falls directly out of the derivation: CTC_T cancels completely in the Henderson-Hasselbalch expression once concentrations are written as fCTfC_T and (1−f)CT(1-f)C_T, so pH depends only on ff and pKapK_a, never on CTC_T itself. What does change with dilution is the plateau's flatness in absolute pH-per-milliliter terms, because β=CTln⁡10 f(1−f)\beta = C_T\ln10\,f(1-f) scales linearly with CTC_T — more dilute solutions are more easily pushed off their buffered plateau by the same volume of titrant.

How does this titration-curve framework connect to the isoelectric point of a full protein, not just a free amino acid?

A folded protein carries dozens of ionizable groups — every Asp, Glu, His, Cys, Tyr, Lys, and Arg side chain, plus the terminal α-amino and α-carboxyl groups at the two ends of the chain. The net charge as a function of pH is computed by summing the same Henderson-Hasselbalch degree-of-deprotonation function αi(pH)\alpha_i(\text{pH}) over every group, weighted by its charge sign, exactly as done here for single amino acids. Setting that sum to zero and solving numerically (since dozens of overlapping pKapK_a values make an algebraic average impossible) gives the protein's isoelectric point, the working principle behind isoelectric focusing gels used throughout proteomics.

Why do textbook pKapK_a tables sometimes list slightly different values for the same amino acid group?

Measured pKapK_a values depend on temperature, ionic strength, and the specific method used (potentiometric titration, NMR chemical-shift titration, spectrophotometric titration for groups with a UV/visible chromophore like tyrosine's phenol or histidine's imidazole). Values reported in different reference tables can differ by a few tenths of a pH unit for this reason. What is robust across all conditions is the relative ordering of pKapK_a values (α-COOH lowest, then acidic side chains, then histidine's imidazolium, then α-NH₃⁺, then cysteine's thiol, tyrosine's phenol, and finally arginine's guanidinium highest), and the qualitative shape of the titration curve this produces.

What happens to the titration curve of histidine specifically, given that its side-chain pKR≈6.0pK_R \approx 6.0 is close to physiological pH?

Histidine's imidazole side chain is unusual among the 20 standard amino acids because its pKRpK_R (about 6.0) sits close to physiological pH (7.4), meaning the side chain exists as a significant mixture of protonated (imidazolium, charged) and deprotonated (imidazole, neutral) forms under cellular conditions — most other side-chain pKapK_a values sit far from 7.4 and are therefore almost entirely in one form. On histidine's titration curve this appears as the second buffering plateau (after the α-COOH near pH 1.8) occurring at a notably lower pH than in other basic amino acids, making histidine's curve visually distinct: two closely spaced low-pH events (α-COOH near 1.8, imidazolium near 6.0) followed by a widely separated third event (α-NH₃⁺ near 9.2). This proximity to physiological pH is exactly why histidine is so frequently found at enzyme active sites acting as a general acid/base catalyst.

Before this

Where this leads