Heavy-atom labeling

Heavy-atom labeling introduces atoms of relatively high atomic number at ordered crystal sites, giving a strong X-ray scattering contribution so that their signal can supply information for solving its phase problem. In X-ray crystallography, measured intensities constrain the magnitudes of structure factors, while their phases are also needed to reconstruct electron density. A heavy-atom derivative provides an additional measurement linked to a model of the label sites.

Under an ideal single isomorphous replacement model, the derivative keeps the parent structure unchanged apart from the added label contribution. For one reflection, write the complex amplitudes asThe heavy-atom model supplies $F_H$, while measurements supply $|F_P|$ and $|F_{PH}|$. These are vector constraints in the complex plane; subtracting intensity magnitudes does not directly reveal the protein phase.^[ref]

For an original example, take $|F_P|=5$, $|F_H|=3$ and $|F_{PH}|=7$ in consistently scaled amplitude units. If $\Delta \varphi ={\varphi }_P-{\varphi }_H$, thenso $\cos \Delta \varphi =1/2$ and $\Delta \varphi =\pm \pi /3$. The one-derivative measurement leaves two possible phases. In particular, $7-5=2$ is not the heavy-atom amplitude3: amplitudes must be combined with their phases.

Multiple isomorphous replacement adds measurements from further derivatives with consistently referenced label positions. These can distinguish the candidate phases in favorable circumstances.^[ref] Large changes in the crystal or protein would violate the stated parent-plus-label model, so derivatization must be assessed rather than assumed to be structurally harmless. Ordered occupancy and useful diffraction matter alongside the chemical act of adding a heavy atom.

Multiple isomorphous replacement combines a native X-ray diffraction data set with data from two or more heavy-atom derivatives to estimate the native structure factor phases. In the ideal crystal isomorphism model, each derivative adds a constraint $|F_P+F_{H,j}|=|F_{PH,j}|$. Consistently referenced heavy-atom models can distinguish the two alternatives left by single isomorphous replacement.^[ref]

Continue the original example in single isomorphous replacement, where $F_P=4\pm 3i$ both fit the native amplitude5 and first derivative with $F_{H,1}=3$. Choose a second reference $F_{H,2}=2i$. Its predicted squared amplitudes differ:A measured second-derivative amplitude $\sqrt{41}$ therefore selects $4+3i$ in this exact model. If the second reference had instead been real, both candidates would still have equal magnitudes after adding it. An additional derivative must contribute useful geometry, rather than merely repeat a constraint.

The heavy-atom models must share a consistent coordinate origin and handedness. Resolving a circle-crossing ambiguity does not by itself determine absolute handedness: conjugating all the amplitudes preserves their magnitudes. CCP4's diagrams explicitly retain this distinction.^[ref]

The method has a concrete experimental history. The opening of the1960 myoglobin structure report describes several mercury- or gold-containing derivatives whose X-ray diffraction patterns were compared with the unsubstituted protein. These measurements had supplied phases for the earlier map at spacings greater than6Å; the new report advanced the analysis to2Å.^[ref] Different labeled crystals supplied information unavailable from native intensities alone.

Single isomorphous replacement estimates phases in X-ray crystallography by comparing a native crystal with one heavy-atom derivative. Under the ideal crystal isomorphism assumption, a reflection's derivative structure factor is $F_{PH}=F_P+F_H$, where the heavy-atom model supplies the complex reference $F_H$. The measured magnitudes constrain $F_P$ to two circles in the complex plane: one centered at zero with radius $|F_P|$, the other centered at $-F_H$ with radius $|F_{PH}|$. A general intersection gives two phase choices.^[ref]

An original exact example uses $|F_P|=5$, $F_H=3$ real and positive, and $|F_{PH}|=\sqrt{58}$. Writing $F_P=x+iy$ givesSubtracting the equations gives $6x+9=33$, hence $x=4$ and $y=\pm 3$. The alternatives areBoth alternatives predict exactly the same native and derivative intensities. Greater precision in these two magnitudes alone cannot select between them. The phase is the missing direction of a complex amplitude, not an unmeasured additive intensity.

Multiple isomorphous replacement can supply an additional circle from another derivative. The numerical example there resolves these particular candidates. In actual measurements, uncertainty, imperfect crystal isomorphism and an imperfect heavy-atom model broaden the constraints, so exact circle crossings are an explanatory idealization.

A heavy-atom derivative is a modified crystal containing ordered heavy-atom sites that can provide a reference scattering contribution for X-ray crystallography. The crystal before modification is called native. For single isomorphous replacement or multiple isomorphous replacement, the desired crystal isomorphism means that the parent crystal structure is sufficiently preserved for its contribution to be shared between data sets. Merely detecting the added chemical element does not establish that model.

The early myoglobin work provides an experimental example. Several crystals with mercury or gold at single sites on the protein were compared with unsubstituted crystals to infer phases. These are chemically modified reference specimens, not computational substitute structures.^[ref]

Crystallographic site occupancy describes the fraction present at a modeled site. The IUCr heavy-atom dictionary defines an occupancy separately for each site and derivative.^[ref] In a simplified single-site model with fixed position and displacement behavior, let $A_H$ be the full-occupancy complex amplitude and let $q$ be occupancy. ThenFor an original illustrative reflection, choose $F_P=4+3i$ and $A_H=6$. With $q=1/2$, the derivative amplitude is $7+3i$ and its squared magnitude is58. At full occupancy it is $10+3i$ and its squared magnitude is109. The native squared magnitude is25. Thus doubling occupancy doubles this modeled heavy-atom amplitude, but does not generally double the measured intensity change: interference with the parent contribution matters.

Useful derivative assessment therefore concerns ordered scattering and comparability as well as chemical uptake. If modification changes the parent structure substantially, the assumed decomposition into unchanged parent plus added reference is inadequate. Heavy-atom labeling describes the motivation for introducing the reference; a heavy-atom derivative is the resulting experimental specimen.