Metal Binding Affinity and Selectivity

Introduction

A protein site does not simply "choose" its metal. Every metal ion in the cell competes with protons, with other metals and with thousands of other ligands for the same donor atoms. Which ion ends up in the site depends on binding constants, on how much of each metal is actually available, and on how fast ions arrive and leave. This page turns the language of stability constants from coordination chemistry into a tool for predicting metal occupancy in living systems, and shows why pure thermodynamics would often put the wrong metal in the wrong protein.

Core explanation

Affinity as an equilibrium. For a protein P binding a metal ion M, P + M ⇌ PM, the association constant is Ka = [PM] / ([P][M]) and the dissociation constant is Kd = 1/Ka. Kd has units of concentration and has a useful meaning: it is the free metal concentration at which half the sites are occupied. The fraction of occupied sites is θ = [M] / (Kd + [M]). A zinc site with Kd = 10⁻¹¹ M is half-filled when free Zn²⁺ is about 10 pM, and nearly full when free Zn²⁺ is ten times higher.

Conditional constants. Binding constants measured in the laboratory are often "conditional": they include the effects of pH and buffer. Protons compete with metals for histidine and cysteine donors, because a protonated imidazole or thiol has no free lone pair. Lowering the pH therefore weakens apparent metal binding. Buffers such as phosphate or Tris can themselves bind metals, lowering the free concentration. A Kd quoted without its pH and competing ligands is incomplete, and comparing values from different conditions can mislead.

Competition between metals. When two metals compete for one site, the ratio of occupancies depends on both affinity and free concentration: [PM₁]/[PM₂] = (Ka,1[M₁]) / (Ka,2[M₂]). A metal with weaker binding can still win if it is far more available. This is the central idea behind biological selectivity.

The Irving-Williams series. For high-spin divalent first-row ions binding a given ligand set, stability constants follow the order Mn²⁺ < Fe²⁺ < Co²⁺ < Ni²⁺ < Cu²⁺ > Zn²⁺. Two factors explain it. Ionic radius decreases across the row, increasing electrostatic attraction, and ligand field stabilisation energy rises from zero at d⁵ Mn²⁺ to a maximum near d⁸ Ni²⁺. Cu²⁺ (d⁹) gains extra stabilisation from Jahn-Teller distortion, which strengthens four equatorial bonds. Zn²⁺ (d¹⁰) has no ligand field stabilisation, so it falls back below copper, though it still binds tightly.

The selectivity problem. The series holds for most N, O and S donor sets, so a protein cannot easily build a site that prefers Mn²⁺ over Cu²⁺ by donor choice alone. Geometry, donor type and the rigidity of the protein fold can shift preferences somewhat, but not by the many orders of magnitude needed. If all metals were equally available, copper and zinc would fill most sites. Cells solve this by controlling availability: free Cu⁺ and Cu²⁺ are held at extremely low concentrations, free Zn²⁺ in the picomolar range, and Mn²⁺ and Fe²⁺ at higher levels. The result is an inverted availability that roughly cancels the Irving-Williams affinity order.

Kinetics and delivery. Some sites are loaded before the protein finishes folding, or in a separate cellular compartment where a competitor is scarce. Metallochaperones hand a specific metal directly to a target by ligand exchange between protein partners, avoiding a free-ion step. Once a metal is buried inside a folded protein, its exchange may be very slow, so the site can remain in a non-equilibrium state that reflects where and when it was filled.

Formulae

Kd = [P][M] / [PM]; fractional occupancy θ = [M] / (Kd + [M]). For two competing metals, [PM₁]/[PM₂] = (Kd,2 [M₁]) / (Kd,1 [M₂]). ΔG° = −RT ln Ka = RT ln Kd.

Step-by-step reasoning

1. Write the binding equilibrium and identify Kd for each candidate metal at the relevant pH. 2. Estimate the free (not total) concentration of each metal in the compartment where the site is loaded. 3. Calculate θ for each metal alone, then the occupancy ratio when they compete. 4. Ask whether the site equilibrates, or whether delivery by a chaperone or slow exchange controls the outcome. 5. Compare the prediction with measured metal content of the purified or in-cell protein.

Visual explanation

Draw a horizontal axis of free metal concentration on a log scale. Plot sigmoidal occupancy curves for a site binding Cu²⁺, Zn²⁺ and Mn²⁺, with the copper curve furthest left and manganese furthest right. Then mark vertical lines showing each metal's typical free concentration in the cytosol. The lines fall so that each curve is only partly or appropriately filled, showing how availability offsets affinity.

Real-world analogy

Imagine a club where the most eager guests are also the rarest. If the doorman only let in whoever arrived most keen, the rare enthusiasts would crowd out everyone else — unless the city keeps those enthusiasts almost entirely at home. Controlling who is on the street matters as much as who wants to come in. Unlike guests, ions have no preferences; the "eagerness" is simply a binding constant.

Real-world example

In cyanobacteria, the manganese-binding protein MncA and the copper protein CucA both end up in the periplasm. Studies showed that MncA binds several metals more tightly than Mn²⁺ in the test tube, yet it acquires manganese because it folds in the cytoplasm, where Mn²⁺ is relatively available and competing ions such as Cu⁺ and Zn²⁺ are tightly buffered, before export. The folded site then traps the metal, illustrating control by location and kinetics rather than affinity alone.

Why?

Why does Zn²⁺ often bind more tightly than Fe²⁺ even though both are divalent and similar in size? Fe²⁺ high-spin gains modest ligand field stabilisation, but Zn²⁺ is slightly smaller and forms more covalent bonds with thiolate and imidazole donors, and its flexible d¹⁰ configuration tolerates tetrahedral protein sites well. The net result places Zn²⁺ above Fe²⁺ in typical binding constants.

Common misconception

"A metalloprotein contains the metal it binds most tightly." Many enzymes bind a non-native metal more strongly than their functional one in isolation. Occupancy in the cell reflects free metal concentrations, compartment and delivery pathways as well as Kd. Another error is to use total cellular metal in θ calculations instead of the free, exchangeable concentration.

Worked example

Question: A site has Kd(Zn²⁺) = 1 × 10⁻¹² M and Kd(Mn²⁺) = 1 × 10⁻⁶ M. The free concentrations are [Zn²⁺] = 1 × 10⁻¹¹ M and [Mn²⁺] = 1 × 10⁻⁵ M. What is the ratio of zinc-bound to manganese-bound sites?

Reasoning: Ratio = (Kd,Mn × [Zn]) / (Kd,Zn × [Mn]) = (10⁻⁶ × 10⁻¹¹) / (10⁻¹² × 10⁻⁵) = 10⁻¹⁷ / 10⁻¹⁷ = 1.

Answer: About 1 : 1. A million-fold stronger affinity for zinc is exactly cancelled by a million-fold higher manganese availability.

Quick check

1. What does a Kd of 10⁻⁹ M tell you about a metal site at a free metal concentration of 10⁻⁹ M? Answer: The site is half occupied, because Kd equals the free concentration giving 50% occupancy.

Exam focus

State the Irving-Williams order correctly, including the peak at Cu²⁺, and explain it with ionic radius and ligand field stabilisation. Use free, not total, metal concentrations in occupancy calculations. Be ready to explain why pH changes apparent affinity and to calculate competitive occupancy ratios from Kd values.

Advanced insight

Metal-sensing transcription factors have tuned Kd values that match the free metal concentration the cell wants to maintain; each sensor switches gene expression when its metal rises above or falls below a set point. Measuring these sensor affinities has allowed researchers to estimate "availability" as a free energy for each metal in a cell, and to design engineered proteins that load the intended metal by matching their affinity to that buffered pool.

Summary

Metal occupancy of a protein site depends on Kd, pH and the free concentration of every competing ion. The Irving-Williams series, Mn²⁺ < Fe²⁺ < Co²⁺ < Ni²⁺ < Cu²⁺ > Zn²⁺, means donor sets alone cannot select weakly binding metals. Cells overcome this by keeping strongly binding metals scarce, loading sites in suitable compartments and using chaperones and slow exchange to trap the correct metal.

Practice questions

1. Explain why stability constants rise from Mn²⁺ to Cu²⁺ in the Irving-Williams series. Answer: Ionic radius decreases, increasing electrostatic attraction, and ligand field stabilisation increases towards d⁸, with Jahn-Teller distortion adding extra stabilisation for d⁹ Cu²⁺. 2. Why does lowering the pH usually weaken metal binding to a histidine or cysteine site? Answer: Protons compete for the donor lone pairs; protonated imidazole or thiol groups cannot coordinate the metal as effectively. 3. A site has Kd = 10⁻⁸ M for a metal whose free concentration is 10⁻⁷ M. Calculate the fractional occupancy. Answer: θ = 10⁻⁷ / (10⁻⁸ + 10⁻⁷) = 0.91, so about 91% of sites are occupied. 4. Give two ways a cell can ensure a manganese enzyme receives Mn²⁺ rather than Zn²⁺. Answer: Keep free Zn²⁺ far lower than free Mn²⁺, and fold or load the enzyme in a compartment or via a delivery protein where manganese is available and zinc is scarce.