Cytochromes and Heme Electron Transfer

A single amino acid substitution eighteen ångströms from an iron atom can shift a cytochrome's reduction potential by 100 mV. That's enough to redirect an electron down one branch of a respiratory chain instead of another. This is the central puzzle of heme-based electron transfer. The iron center itself, cycling quietly between Fe(II) and Fe(III), is almost bystander-simple. Meanwhile, the protein wrapped around it does all the interesting physical chemistry. Understanding how that wrapping works — which axial ligands, which pocket polarities, which solvent exposures — is the difference between a rusting nail and a mitochondrion. Cytochromes are the workhorses of biological respiration and photosynthesis. They ferry single electrons over distances too great for direct collision chemistry. They use the iron–porphyrin unit introduced in Porphyrins and Metalloporphyrin Chemistry as an outer-sphere electron reservoir. Unlike the oxygen-carrying hemoglobins, cytochromes almost never bind exogenous small molecules at iron. The axial coordination sites are permanently capped by protein side chains. The entire redox event is a clean one-electron Fe(III)/Fe(II) couple with no bond breaking at the metal. That structural fact — a coordinatively saturated, kinetically inert metal center — is precisely what makes the iron a good electron-transfer switch rather than a catalytic center. It is why cytochromes are the paradigm system for the physical framework built in Electron-Transfer Metalloproteins. The a/b/c classification was introduced by David Keilin from optical absorption spectra in the 1920s–30s. It turns out to encode real chemical distinctions. Different heme substituents exist (b = protoporphyrin IX, c = covalently thioether-linked heme, a = formyl- and farnesyl-substituted heme). Different axial ligand sets also exist (His/His, His/Met, His/imidazolate). Together these span a potential range from about −400 mV to +400 mV. This range emerges using chemically the same Fe–porphine core. This article works through the electronic structure of low-spin heme iron. It covers the four principal levers a protein pulls on to tune $E°'$ . It also covers the machinery — cytochrome $c$ , the $bc_1$ complex, cytochrome $c$ oxidase — that assembles these tuned pieces into a working respiratory chain. The rate theory governing how fast electrons actually hop between these tuned centers is deferred to Marcus Theory and Biological Electron Transfer. But every structural feature discussed here is precisely a knob on that rate.

Key takeaways

Cytochromes use Fe-protoporphyrin (heme) held in low-spin, six-coordinate geometry with no bond making/breaking — pure outer-sphere Fe(III)/Fe(II) electron transfer.
Class letters (a, b, c) denote heme substituents/attachment (formyl+farnesyl; unmodified protoporphyrin IX; thioether-linked to Cys), not just optical bands, though Soret/Q-band shifts historically defined them.
Axial ligand set is the dominant tuning knob: bis-His hemes sit near 0 to +100 mV, His/Met (cytochrome $c$ ) near +250 mV, bis-Met or His/Lys can go higher or lower depending on partner ligand basicity.
Heme solvent exposure and local electrostatics (charged residues, dipoles, heme propionate ionization) shift $E°'$ by up to several hundred mV by differentially stabilizing Fe(III) vs. Fe(II).
Hydrogen bonding to axial ligands (e.g., to the His Nδ or Met lone pair) and heme ruffling/doming perturb metal–ligand overlap and add tens of mV of fine control.
Respiratory and photosynthetic chains stack cytochromes of increasing $E°'$ so that electron flow is thermodynamically downhill at every step, converting free energy into a proton-motive force.
Key idea: the iron/porphyrin chemistry is nearly invariant across cytochromes — essentially all redox tuning is done by the protein matrix acting on a fixed low-spin Fe(III)/Fe(II) couple.

Where this lives today

Cytochrome $c$ is the textbook electron shuttle between Complex III (cytochrome $bc_1$ ) and Complex IV (cytochrome $c$ oxidase) in the mitochondrial inner membrane, and its release from mitochondria into the cytosol is also the trigger for apoptosis — cytosolic cytochrome $c$ binds Apaf-1 to nucleate the apoptosome, meaning the same tuned Fe(III)/Fe(II) heme that ferries electrons in healthy cells doubles as a cell-death signal once mislocalized. Cytochrome $c$ oxidase's binuclear heme $a_3$ /Cu $_B$ center, tuned to a high potential exactly so it can be the terminal electron sink that reduces O $_2$ to water, is the target of cyanide and azide poisoning (both bind the ferric heme $a_3$ site, blocking O $_2$ reduction and halting aerobic respiration). Photosynthetic electron transport uses an analogous cytochrome, cytochrome $f$ (a c-type heme despite the name, tuned to roughly +365 mV), inside the cytochrome $b_6f$ complex to link photosystem II to photosystem I. Engineered and de novo heme proteins — "maquettes" built from simple four-helix bundles with bis-His or His/Met ligation designed in — are now used as a testbed to verify the tuning rules described here quantitatively, since they let chemists vary axial ligand and pocket polarity independently in a minimal scaffold. Cytochrome $c$ 's extreme thermodynamic and kinetic stability alongside its well-defined, tunable potential has also made it a standard component in bioelectrochemistry and biosensor design, where surface-immobilized cytochrome $c$ films are used as a well-behaved redox couple for calibrating protein-film voltammetry — a technique whose rate interpretation depends directly on the reorganization-energy ideas that this heme chemistry sets up for Marcus Theory and Biological Electron Transfer.

Ask the tutor: How does the small reorganization energy of cytochrome hemes quantitatively translate into their fast Marcus-theory electron-transfer rates compared to a typical small-molecule inorganic redox couple?

The physical picture

Every cytochrome heme is a low-spin, six-coordinate iron center: four equatorial nitrogens from the porphyrin macrocycle, and two axial ligands supplied by the protein. In the resting oxidized state, Fe(III) is $d^5$ ; in the reduced state, Fe(II) is $d^6$ . Both states are essentially always low-spin in functional cytochromes, meaning the $t_{2g}$ -like set ( $d_{xy}, d_{xz}, d_{yz}$ in the local octahedral-ish field) holds all five or six electrons, and the $e_g$ -like $\sigma^*$ orbitals ( $d_{z^2}, d_{x^2-y^2}$ ) stay empty. This is the electronic-structure fact that makes cytochromes fast, reversible one-electron switches: because no strongly antibonding orbital is populated in either oxidation state, the Fe–N bond lengths and the overall geometry change only slightly on reduction (typically a contraction of a few hundredths of an ångström at Fe–N $_{porph}$ , and a small stiffening of the axial bonds). Small geometric reorganization means a small inner-sphere reorganization energy $\lambda_i$ , which is exactly the ingredient that makes electron self-exchange and cross reactions fast in the Marcus framework used downstream. The axial ligand pair is the single most important structural variable across the three classes. Cytochrome $b$ hemes (and many others) use bis-histidine axial ligation — two neutral imidazole nitrogens donating into the $d_{z^2}$ -type orbital. Cytochrome $c$ uses a His/Met pair: one imidazole nitrogen and one thioether sulfur, an asymmetric, weaker-field, more polarizable ligand set. Cytochrome $a$ (as in cytochrome $c$ oxidase) pairs one histidine with a second histidine that is often deprotonated to an anionic imidazolate character in some contexts, or otherwise sits in a highly asymmetric heme (heme $a$ , bearing a formyl group conjugated into the porphyrin π system and a long farnesyl hydrophobic tail). Each substitution changes the σ-donor/π-acceptor balance at iron and therefore the relative stabilization of $d^5$ vs. $d^6$ — which is exactly the reduction potential. Superimposed on the ligand-field picture is electrostatics: the heme propionates, nearby charged residues (Lys, Arg, Asp, Glu), backbone dipoles, and the degree to which the heme edge is exposed to bulk solvent versus buried in a hydrophobic pocket. Because Fe(III) is more positively charged than Fe(II), anything that stabilizes positive charge (nearby negative charges, high dielectric/solvent exposure) will stabilize Fe(III) relative to Fe(II) and lower $E°'$ (make oxidation easier, i.e., destabilize the reduced state relatively less — precisely, it raises the energy of removing an electron less, so $E°'$ falls). Conversely, burying the heme in a low-dielectric, positively-polarized pocket stabilizes Fe(II) relative to Fe(III) and raises $E°'$ . This electrostatic tuning, layered on top of the ligand-field tuning, is how nature spans nearly a full volt of chemically similar Fe/heme couples using one invariant metal-ligand core.

Figure

Low-spin Fe(III) and Fe(II) d-orbital occupancy in a heme with axial His/Met ligation

Figure

Reduction potential ladder: axial ligation and pocket electrostatics stack to span ~800 mV

The core ideas

The Nernst relation anchors everything quantitatively:

E = E°' - \frac{RT}{nF}\ln\frac{[\text{Fe(II)}]}{[\text{Fe(III)}]}

where: $E°'$ is the standard (biochemical, pH 7) reduction potential of the Fe(III)/Fe(II) couple, $R$ is the gas constant, $T$ the temperature, $n=1$ for a one-electron heme couple, $F$ the Faraday constant, and the log ratio is reduced-to-oxidized heme concentration (or activity).

Reading: at $[\text{Fe(II)}]=[\text{Fe(III)}]$ the electrode/heme potential equals $E°'$ exactly; this is the number tabulated for each cytochrome and the number the protein matrix is tuning. The free energy of electron transfer between two hemes of differing potential sets the thermodynamic driving force for a respiratory chain step:

\Delta G°' = -nF\,\Delta E°' = -nF\left(E°'_{\text{acceptor}} - E°'_{\text{donor}}\right)

where: $\Delta E°'$ is the potential difference between the electron acceptor and donor couples, and $\Delta G°'$ is the standard free energy released per electron transferred; a negative $\Delta G°'$ (i.e., acceptor potential higher than donor potential) means the transfer is thermodynamically favorable.

Reading: a chain of cytochromes with monotonically increasing $E°'$ (cyt $b \to$ cyt $c_1 \to$ cyt $c \to$ cyt $a \to$ cyt $a_3 \to$ O $_2$ ) guarantees each hop is downhill, converting the overall $\Delta G°'$ of NADH/O $_2$ chemistry into a series of small, kinetically manageable steps rather than one large uncontrolled drop. The ligand-field/electrostatic tuning itself is best organized as a thermodynamic cycle (Born–Haber-like) separating intrinsic (gas-phase/ligand-field) contributions from protein-environment contributions:

E°'_{\text{protein}} = E°'_{\text{intrinsic}} + \Delta E°'_{\text{axial ligand}} + \Delta E°'_{\text{electrostatic}} + \Delta E°'_{\text{H-bond/solvent}}

where: $E°'_{\text{intrinsic}}$ is a reference value for Fe-porphine with a fixed, idealized ligand set, and each $\Delta E°'$ term is the perturbation from axial ligand identity, local charge/dipole environment, and hydrogen bonding or solvent accessibility, respectively.

Reading: this decomposition is not exact quantum chemistry but a bookkeeping device chemists use constantly — it lets you predict the direction and rough magnitude of a mutation's effect on $E°'$ from structural reasoning alone, which is exactly what site-directed mutagenesis studies on cytochrome $c$ have repeatedly validated.

Derivation AWhy axial ligand basicity controls E°′: a σ-donation argument

Goal: show, from ligand-field/donor-strength reasoning, why replacing an axial His (neutral imidazole) with a more σ-donating ligand (e.g., anionic imidazolate, or a stronger-field ligand) lowers the Fe(III)/Fe(II) reduction potential, while a weaker-donor/more polarizable ligand (e.g., thioether Met) raises it relative to bis-His.

Step 1 — identify what changes upon reduction. Fe(III) ( $d^5$ , low spin) has a net positive charge and a half-filled $t_{2g}$ -type set; Fe(II) ( $d^6$ , low spin) has one more electron, filling $t_{2g}$ completely, with the $\sigma^*$ ( $d_{z^2}$ / $d_{x^2-y^2}$ ) orbitals empty in both states. The reduction is effectively adding an electron into a metal-centered, largely non-bonding $t_{2g}$ -type orbital that has some $\pi$ -backbonding character with the axial ligand and the porphyrin.

Step 2 — consider a strong σ-donor axial ligand (e.g., imidazolate, or a strongly basic amine). A strong σ-donor raises the energy of metal orbitals with which it overlaps most, and — more importantly for the redox thermodynamics — it stabilizes higher formal positive charge on iron by donating electron density into empty/partially-empty metal orbitals. Fe(III), being more electron-poor, benefits disproportionately from a strong donor: the ligand "pushes" electron density onto iron, which is exactly what a $d^5$ , cationic, electron-hungry center wants. This preferentially stabilizes the Fe(III) state relative to Fe(II).

Step 3 — translate stabilization into a potential shift. Reduction potential is (by definition, from Step Goal's cycle) related to the free-energy difference $G(\text{Fe(III)}) - G(\text{Fe(II)})$ : a more negative value of this quantity (Fe(III) state more stable) corresponds to harder reduction, i.e., a lower $E°'$ . Because a strong σ-donor axial ligand disproportionately stabilizes Fe(III) (Step 2), it lowers $G(\text{Fe(III)})$ more than $G(\text{Fe(II)})$ , and therefore lowers $E°'$ .

\Delta E°' \propto -\left[\Delta G_{\text{stab}}(\text{Fe(III)}) - \Delta G_{\text{stab}}(\text{Fe(II)})\right]

where: $\Delta G_{\text{stab}}$ is the extra stabilization from the stronger donor in each oxidation state; since Fe(III) is stabilized more, the bracket is positive, and $\Delta E°'$ (the shift relative to a weaker-donor reference) is negative.

Step 4 — apply to His vs. Met vs. imidazolate. Histidine imidazole is a moderate σ-donor; thioether Met is a weaker σ-donor but more polarizable π-acceptor-tolerant ligand that does not stabilize Fe(III) as strongly. Replacing one His with Met (as in cytochrome $c$ 's His/Met pair versus a hypothetical bis-His analog) therefore removes some of the extra Fe(III) stabilization that a second His would have provided, raising $E°'$ — consistent with cytochrome $c$ 's characteristic +250 mV versus roughly 0 to −100 mV for many bis-His cytochromes $b$ . Conversely, deprotonating an axial His to imidazolate (a strong anionic σ-donor, seen transiently in some peroxidase-type or engineered systems) stabilizes Fe(III) strongly and drives $E°'$ sharply negative — by several hundred mV in extreme cases.

Step 5 — state the qualitative rule. Stronger axial σ-donation → lower $E°'$ (Fe(III) preferentially stabilized, reduction harder). Weaker/softer axial donation (thioether, weak-field amines) → higher $E°'$ (Fe(II) relatively favored, reduction easier). This single rule organizes the entire bis-His/His-Met/His-Lys spread of biological heme potentials.

Go deeper · Derivation B · 2 minElectrostatic tuning by heme burial and charged residues

Derivation BElectrostatic tuning by heme burial and charged residues

Goal: quantify, using a simple Born-type solvation/Coulomb argument, how burying the heme edge in a low-dielectric pocket (versus solvent exposure) and placing a charged residue near the propionates shifts $E°'$ , and show the direction of the shift is consistent with observed mutagenesis data on cytochrome $c$ and cytochrome $b_5$ .

Step 1 — model the heme charge states. Treat Fe(III)-heme (net charge $+1$ relative to Fe(II)-heme for the iron alone, though the full macrocycle/propionate charges are more complex) and Fe(II)-heme as two point-charge-like states differing by one unit of positive charge, $\Delta q = +1$ , at a fixed position (the iron). The reduction event moves the "extra charge" from iron.

Step 2 — apply the Born solvation formula for placing charge in a dielectric. The electrostatic free energy of a charge $q$ in a spherical cavity of radius $a$ embedded in a medium of dielectric constant $\varepsilon$ is

G_{\text{Born}} = -\frac{q^2}{8\pi\varepsilon_0 a}\left(1 - \frac{1}{\varepsilon}\right)

where: $q$ is the ionic charge, $a$ the effective cavity radius, $\varepsilon_0$ the vacuum permittivity, and $\varepsilon$ the medium's dielectric constant (large, ~80, for water; small, ~4, for protein interior).

Step 3 — compare the two oxidation states. Fe(III) (larger effective positive charge) has a more negative $G_{\text{Born}}$ (i.e., is more stabilized) in a high-dielectric environment than Fe(II), because $G_{\text{Born}}$ scales as $q^2$ and $q_{\text{Fe(III)}}^2 > q_{\text{Fe(II)}}^2$ relative to a common reference — put simply, solvent (or any high- $\varepsilon$ environment) preferentially stabilizes the more charged, more polar Fe(III) state. This lowers $E°'$ for solvent-exposed hemes.

Step 4 — invert for buried hemes. A heme buried in a hydrophobic, low-dielectric pocket loses this preferential stabilization of Fe(III); the relative free energy of Fe(II) improves (destabilization of Fe(III) is larger than destabilization of Fe(II), because Fe(III)'s stabilization in a high-ε medium was larger to begin with). The reduction becomes thermodynamically easier, and $E°'$ rises. This is precisely why deeply buried, hydrophobically packed hemes (heme $a_3$ /cytochrome oxidase region, many cytochromes $c$ ) sit at higher potentials than similarly-ligated but solvent-exposed hemes.

Step 5 — add a discrete charge (mutagenesis case). Now place a single point charge $q_r$ (e.g., a Lys $^+$ or Asp $^-$ side chain) at distance $r$ from the iron, treated by Coulomb's law with an effective dielectric $\varepsilon_{\text{eff}}$ (typically 10–40 for a protein interior/surface region):

\Delta E°' \approx \frac{1}{nF}\cdot\frac{q_r\,\Delta q}{4\pi\varepsilon_0\varepsilon_{\text{eff}}\, r}

where: $\Delta q = +1$ is the charge change on iron upon oxidation (Fe(II) $\to$ Fe(III)), $q_r$ is the sign and magnitude (in elementary charge units) of the introduced residue charge, $r$ is the distance from that charge to the iron, and the sign convention is such that a nearby negative charge ( $q_r<0$ ) lowers $E°'$ (stabilizes Fe(III), the more positive state, making oxidation easier/reduction harder) while a nearby positive charge raises $E°'$ .

Step 6 — sanity-check magnitude. For $q_r = -1$ (a single anionic Asp), $r = 8\times10^{-10}\ \text{m}$ , and $\varepsilon_{\text{eff}}=20$ :

\Delta E°' \approx \frac{(1.60\times10^{-19})^2}{4\pi(8.854\times10^{-12})(20)(8\times10^{-10})(96485)}\ \text{V}

Evaluating the constants gives $\Delta E°' \approx -90\ \text{mV}$ , of the same order as experimentally observed shifts (tens to ~100 mV) for single charged-residue mutations near cytochrome hemes — confirming the sign and rough scale of the electrostatic lever without needing full electronic structure theory.

Worked examples

1. Foundational — assigning spin state and counting the redox-active electron. Cytochrome $b_5$ carries a bis-His-ligated heme. In the Fe(III) resting state, is the iron more likely low-spin or high-spin, and how many unpaired electrons does that imply? Two strong-to-moderate field imidazole nitrogens plus the four porphyrin nitrogens constitute a strong, symmetric octahedral-like ligand field; bis-His hemes in cytochromes are essentially always low-spin. For $d^5$ low-spin, the configuration is $t_{2g}^5 e_g^0$ , giving one unpaired electron, $S=1/2$ — consistent with the characteristic EPR signal ( $g$ values around 2.9, 2.3, 1.5, the classic low-spin ferriheme rhombic pattern) used diagnostically to confirm bis-His ligation in unknown heme proteins. Upon reduction to Fe(II), the configuration becomes $t_{2g}^6 e_g^0$ , $S=0$ , EPR silent — the basis for using EPR to titrate cytochrome oxidation state in mitochondrial extracts.

2. Applied — predicting the direction of a potential shift from a real mutation. In yeast iso-1-cytochrome $c$ , replacing the wild-type Asn52 (near the heme edge, hydrogen-bonding to a heme propionate) with an Ile removes a hydrogen bond that normally helps neutralize/stabilize the propionate's negative charge, indirectly affecting the electron density at iron. Using the electrostatic argument (Step 5 above): removing a stabilizing hydrogen bond to a negatively charged propionate is like removing a nearby effective negative charge's stabilizing interaction with the metal. Since a nearby negative charge normally lowers $E°'$ (stabilizes Fe(III)), removing/weakening that interaction should raise $E°'$ — the reduced state is now relatively favored because the destabilizing effect on Fe(III) diminishes less than it does on Fe(II) once the compensating H-bond is gone... more directly: losing stabilization of the anionic propionate makes the whole heme environment slightly more electropositive-tolerant, favoring Fe(II). Experimentally, Asn52Ile-type mutants of iso-1-cytochrome $c$ shift $E°'$ upward by several tens of millivolts (reported values cluster around +20 to +50 mV depending on exact construct), in the direction predicted by the electrostatic-lever argument, and this family of mutations is a standard teaching example for how propionate hydrogen-bonding fine-tunes potential without touching the axial ligands at all.

3. Exam/competition-level — building a full three-heme thermodynamic ladder and computing overall driving force. A designed electron-transfer chain has three sequential hemes: heme A ( $E°'_A = -0.15\ \text{V}$ , bis-His, solvent-exposed), heme B ( $E°'_B = +0.10\ \text{V}$ , bis-His, partially buried), and heme C ( $E°'_C = +0.26\ \text{V}$ , His/Met, buried, cytochrome- $c$ -like), wired A $\to$ B $\to$ C $\to$ O $_2$ (using $E°'_{\text{O}_2/\text{H}_2\text{O}} = +0.815\ \text{V}$ at pH 7 as the terminal acceptor). Compute the standard free energy change for each hop and confirm the chain is everywhere downhill, then compute the total free energy released per electron over the whole chain. Hop A→B:

\Delta E°'_{AB} = E°'_B - E°'_A = 0.10 - (-0.15) = +0.25\ \text{V}

\Delta G°'_{AB} = -nF\Delta E°'_{AB} = -(1)(96485\ \text{C/mol})(0.25\ \text{V}) = -24.1\ \text{kJ/mol}

Hop B→C:

\Delta E°'_{BC} = 0.26 - 0.10 = +0.16\ \text{V}, \qquad \Delta G°'_{BC} = -(96485)(0.16) = -15.4\ \text{kJ/mol}

Hop C→O $_2$ :

\Delta E°'_{C,\text{O}_2} = 0.815 - 0.26 = +0.555\ \text{V}, \qquad \Delta G°'_{C,\text{O}_2} = -(96485)(0.555) = -53.5\ \text{kJ/mol}

Every hop has $\Delta E°' > 0$ , hence $\Delta G°' < 0$ throughout — the chain is thermodynamically downhill at each transfer, exactly the design principle. Summing (equivalent to computing directly from the endpoints since $F$ and $n=1$ are common):

\Delta G°'_{\text{total}} = -nF\left(E°'_{\text{O}_2} - E°'_A\right) = -(96485)(0.815-(-0.15)) = -93.1\ \text{kJ/mol}

which indeed equals $-24.1 - 15.4 - 53.5 = -93.0\ \text{kJ/mol}$ (rounding), confirming that the sum of stepwise free energies reproduces the overall driving force — the free energy is a state function independent of how many intermediate hemes the electron passes through, which is exactly why a respiratory chain can be built from many small, kinetically fast, reversible-looking hops without losing any of the total available free energy to the reaction as a whole (it is instead partitioned, hop by hop, partly stored as proton-motive force at the coupling sites and partly dissipated as heat at each electron-transfer step depending on how far $\Delta G°'$ exceeds the reorganization energy at that step).

Limiting cases, scope & common traps

Regime / limit	What happens	Why
Fully solvent-exposed bis-His heme model (e.g., microperoxidase in water)	$E°'$ drops to strongly negative values (often below −200 mV)	High dielectric preferentially stabilizes Fe(III) (Born argument); no protein shielding
Heme deeply buried, hydrophobic pocket, His/Met	$E°'$ rises toward +250 to +300 mV (typical soluble cytochrome $c$ )	Low dielectric removes preferential Fe(III) stabilization; weaker donor (Met) also favors Fe(II)
Axial ligand deprotonated to imidazolate (His $^-$ )	$E°'$ drops sharply, sometimes >300 mV negative vs. neutral His	Strong anionic σ-donor stabilizes Fe(III) strongly (ligand-field argument)
$n \ne 1$ assumed for a heme couple	Formula and interpretation break	Cytochrome hemes are strictly one-electron couples; never treat as 2-electron unless describing a coupled catalytic cycle (e.g., whole $\text{O}_2 \to \text{H}_2\text{O}$ at a binuclear center)
High-spin heme (5- or 4-coordinate, as in some peroxidases/globins)	Larger structural reorganization on redox change	Occupation of $\sigma^$ ( $e_g$ -like) orbitals lengthens Fe–N bonds; this is not* the typical cytochrome situation
Treating $E°'$ as an intrinsic "iron" property	Wrong — ignores protein context entirely	The bare Fe-porphine couple is not what's measured; always decompose into intrinsic + ligand + electrostatic + H-bond terms

Common trapCommon misconception: "cytochromes bind and release electrons the way hemoglobin binds and releases O₂"

Hemoglobin/myoglobin chemistry involves a genuine coordination change at iron (an open axial site binding O $_2$ , CO, or nothing) and often a spin-state change (high-spin deoxy-Fe(II) to low-spin oxy-Fe(II)). Cytochrome electron transfer, by contrast, involves no ligand binding or dissociation at iron at all in the normal catalytic cycle — both axial sites are permanently occupied by protein side chains, and only the metal's formal oxidation state and the accompanying spin-state-conserving electron count change. Students often assume the heme "does chemistry" the same way in both protein families; in cytochromes, the heme's job is purely as an electron reservoir with fixed coordination number, and that coordinative rigidity is the entire reason the inner-sphere reorganization energy is small and electron transfer is fast.

The chemist behind it

Origins · 1 minDavid Keilin and how this was discovered

David Keilin

1887–1963 · Cambridge, Molteno Institute, 1920s–1940s

Keilin, a Polish-born British parasitologist-turned-biochemist, rediscovered and systematically characterized the cytochromes using simple absorption spectroscopy — literally observing the four-banded spectrum of insect flight muscle through a hand spectroscope and recognizing that the bands shifted reversibly with the physiological state of the tissue (respiring versus at rest). He established the a, b, c nomenclature from the position of their characteristic α-absorption bands and demonstrated that these pigments formed a linked chain capable of passing electrons stepwise from substrate to oxygen — the conceptual skeleton of the respiratory chain that biochemistry still uses. His patient, largely single-handed spectroscopic work, done with minimal instrumentation compared to what followed, is a model of how far rigorous observation and clear physical reasoning can go before any crystal structure existed to confirm the chemistry.

History

- 1925 — Keilin rediscovers and names cytochromes a, b, c from tissue absorption spectra, building on earlier 19th-century observations by MacMunn.

- 1930s — Keilin and Hartree establish the respiratory chain concept: sequential cytochrome oxidation/reduction linking substrate dehydrogenation to O

_2

reduction.

- 1958 — Margoliash and coworkers determine the amino acid sequence of horse heart cytochrome

c

, launching comparative sequence/evolution studies of electron-transfer proteins.

- 1971 — Dickerson and coworkers solve the first high-resolution cytochrome

c

crystal structure, revealing the His/Met axial ligation directly.

- 1990s–2000s — Systematic mutagenesis (Bill Durham, Judith Hazzard, and others) and computational electrostatics (Warshel-style continuum models) quantitatively map axial ligand and charge effects on

E°'

across cytochrome families.

Why do cytochromes use iron specifically, rather than another transition metal, for electron transfer?

Iron's Fe(III)/Fe(II) couple sits in a convenient potential window (roughly −400 mV to +400 mV depending on ligands) that matches the biologically relevant range between NADH/NAD $^+$ (~−320 mV) and O $_2$ /H $_2$ O (+815 mV), and iron's abundance made it evolutionarily accessible early on. Additionally, both Fe(III) $d^5$ and Fe(II) $d^6$ are readily stabilized as low-spin states by a porphyrin plus two axial ligands, giving a small structural change on redox cycling — the key to fast, low-barrier electron transfer. Manganese or cobalt porphyrins can be made synthetically and have interesting chemistry, but they do not naturally give this combination of tunable potential window and minimal reorganization energy as cheaply as iron does.

What exactly distinguishes cytochrome a, b, and c chemically, beyond historical optical bands?

Cytochrome $b$ hemes use unmodified protoporphyrin IX (heme $b$ ), non-covalently held in the protein pocket by ligand coordination alone. Cytochrome $c$ hemes are heme $c$ , covalently attached to the protein through thioether bonds from two cysteine residues to the vinyl groups of the porphyrin, making the heme essentially inseparable from the polypeptide without proteolysis. Cytochrome $a$ hemes use heme $a$ , which carries a formyl substituent conjugated into the π system (redshifting the optical spectrum and altering electronics at iron) and a long hydrophobic farnesyl chain that anchors it deep in the membrane. All three still use the same Fe-porphine core electronic structure; the differences are in covalent attachment, substituents, and consequently local electronics/packing.

How does hydrogen bonding to axial ligands fine-tune the potential without changing the ligand itself?

A hydrogen bond donated to the lone pair on an axial ligand (for example, a backbone NH or a side chain hydrogen-bonding to the Nδ of an axial histidine, or to a lone pair on the coordinating methionine sulfur) effectively withdraws electron density from that lone pair before it reaches iron, weakening the ligand's net donor strength as felt by the metal. By the same logic used for changing ligand identity outright, a weaker effective donor stabilizes Fe(III) less strongly, which raises $E°'$ modestly — typically tens of millivolts per hydrogen bond, rather than the hundreds of millivolts seen for full ligand substitution. This is one of the finest-grained tuning mechanisms available to a protein and is frequently used for small potential adjustments between otherwise very similar homologous cytochromes.

Is the Fe(III)/Fe(II) heme couple ever coupled to proton transfer (a "redox-Bohr" effect)?

Yes, in some cytochromes (notably certain multi-heme cytochromes $c_3$ and some cytochrome $c$ oxidase intermediates) the reduction potential is measurably pH-dependent because a nearby ionizable group (often a heme propionate or a nearby His) changes protonation state depending on the iron's oxidation state — the "redox-Bohr effect," named by analogy to hemoglobin's Bohr effect. When this coupling exists, the simple Nernst equation used here needs a pH-dependent correction term, and the electron transfer is formally proton-coupled electron transfer (PCET) rather than pure outer-sphere electron transfer. Most simple monoheme cytochromes like cytochrome $c$ show only weak redox-Bohr behavior, so the pure Nernst treatment is an excellent approximation for them.

Why does cytochrome c leaking into the cytosol trigger apoptosis rather than just causing a metabolic problem?

Cytosolic cytochrome $c$ is recognized specifically by the protein Apaf-1, which upon binding cytochrome $c$ and ATP/dATP undergoes a conformational change and oligomerizes into the apoptosome, a heptameric platform that recruits and activates procaspase-9. This is a repurposing of the same tuned, stable, well-folded heme protein used in respiration — its stability and its precisely tuned surface (which normally docks onto Complex III and Complex IV) turn out to also be a good recognition motif for Apaf-1. The redox chemistry of the heme itself is not directly what triggers apoptosis; it is the protein's release from the mitochondrial intermembrane space (normally prevented by the outer membrane) and its subsequent protein–protein interaction that matters.

Can you predict the exact E°′ of a new cytochrome from sequence alone?

Not exactly, but you can predict the right order of magnitude and the correct qualitative trend using the rules developed here: identify the axial ligand pair first (this sets the largest single contribution, easily 200–400 mV of spread), then assess heme burial/solvent exposure (tens to a couple hundred mV), then look for nearby charged residues and hydrogen bonds to axial ligands or propionates (tens of mV each). Full quantitative prediction (to within ~20 mV) generally requires either continuum electrostatics calculations (Poisson–Boltzmann or similar, parametrized on the specific fold) or empirical correlation with homologous, experimentally characterized proteins — sequence alone, without structural/electrostatic modeling, is not sufficient for high precision.

Why is the inner-sphere reorganization energy important for more than just heme electron transfer?

A small $\lambda_i$ is exactly the structural feature that Marcus theory identifies as necessary (along with modest outer-sphere/solvent reorganization $\lambda_o$ ) for very fast electron self-exchange and cross-reaction rates; cytochromes are in fact one of the primary experimental systems used historically to validate Marcus-theory rate predictions in biology, because their $\lambda$ values are unusually small and well-characterized compared to many small-molecule inorganic couples. This heme chemistry is therefore not just descriptive background — it is the structural reason biological electron transfer chains can move electrons over nanometer distances on microsecond-or-faster timescales without needing large activation barriers, a link made quantitative in the rate theory that follows this material.