Apart from the standard procedures for detecting ligands (difference Fourier maps), BUSTER has one particular feature that needs a bit further explanation - to explain what it can do, what it can't do and what potential bias it might introduce.

This feature is triggered by the –L and –Lpdb flags to the refine command:

% refine -p some.pdb -m some.mtz -L ...
  - or -
% refine -p some.pdb -m some.mtz -Lpdb bindingsite.pdb ...

It will treat a certain region of the model differently during the last big cycle of refinement: that region will neither contain an atomic model nor a contribution from the bulk solvent. However, if there is some electron density present in this region, it should show up in difference (Fo-Fc) maps as strong positive density. The interpretation, what this electron density might represent (atomic model, bulk solvent or a mixture) is up to the user.

Here we're going to explain the typical usage of this BUSTER feature, their assumptions and caveats.

Unknown location

The least biased assumption is that there might be some ligand bound, but its location is unknown. In that case one would use the

% refine -p some.pdb -m some.mtz -L ...

option. It will also trigger the automatic addition of 'dummy' atoms (currently described in terms of 'water' residues). Just before the last big cycle, an analysis of regions badly explained by those dummy atoms is done. This could result in a report (on standard output) about a certain number of waters being removed in a so-called 'interesting' region. This region is then treated as described above.

Known location

If the location of ligand binding is known, a PDB file describing this region can be used explicitly with

% refine -p some.pdb -m some.mtz -Lpdb bindingsite.pdb ...

The same logic as before will be applied just before the last big cycle - with one difference: instead of searching for regions with unexplained difference density, the region described by the user-given PDB file is being used.


Some points need to be considered here:

1. the placed 'waters'

Since the goal of this procedure is to find, highlight and enhance one or more (potential) ligand-binding sites, the placed 'waters' need to be taken not too serious. They are placed to explain density, which could be non-water density (missing loops, disordered side-chains, unmodelled compounds) but also just noise in the data. It is not recommended to use the resulting solvent structure as such.

2. false negatives

It is possible to miss genuine bound compounds - especially if they are weakly bound (partial occupancy and/or partially disordered) or rather small (fragment searches). There are options available for changing the minimum required volume.

3. false positives

It is possible to get false positives with that procedure - especially the –Lpdb flag. Some typical examples are:


While being aware of the limitations and potential biasing issue, the procedures described above can still be very useful. It is important to know what question one is actually asking:

Once those questions have been answered, the ligand can be modelled into the binding site and refined normally.