Difference density map

In X-ray crystallography, a difference density map or Fo–Fc map shows the spatial distribution of the difference between the measured electron density of the crystal and the electron density explained by the current model.^[1]

A way to compute this map has been formulated for cyro-EM.^[2]

Display

Conventionally, they are displayed as isosurfaces with positive density—electron density where there's nothing in the model, usually corresponding to some constituent of the crystal that hasn't been modelled, for example a ligand or a crystallisation adjutant -- in green, and negative density—parts of the model not backed up by electron density, indicating either that an atom has been disordered by radiation damage or that it is modelled in the wrong place—in red. The typical contouring (display threshold) is set at 3σ.^[3]

Calculation

Difference density maps are usually calculated using Fourier coefficients which are the differences between the observed structure factor amplitudes from the X-ray diffraction experiment and the calculated structure factor amplitudes from the current model, using the phase from the model for both terms (since no phases are available for the observed data). The two sets of structure factors must be on the same scale.

${\displaystyle C_{diffmap}=(|F_{obs}|-|F_{calc}|)exp(2\pi i\phi _{calc})}$

It is now normal to also include maximum-likelihood weighting terms which take into account the estimated errors in the current model:

${\displaystyle C_{diffmap}=(m|F_{obs}|-D|F_{calc}|)exp(2\pi i\phi _{calc})}$

where m is a figure of merit which is an estimate of the cosine of the error in the phase, and D is a "σ_A" scale factor. These coefficients are derived from the gradient of the likelihood function of the observed structure factors on the basis of the current model. A difference map built with m and D is known as a mFo – DFc map.^[3]

The use of ML weighting reduces model bias (due to using the model's phase) in the 2 Fo–Fc map, which is the main estimate of the true density. However, it does not fully eliminate such bias.^[4]

References

Further reading

• Electron density maps on Proteopedia