BUSTER frequently-asked questions


I get an error message 'OMP abort...

Answer from Claus Flensburg

OMP abort: Unable to set worker thread stack size to 4194304 bytes

A MakeTNT application does not run.

Answer from Peter Keller.


grade_PDB_ligand gives a message 'Tool babel not found on path'

Answer from Tom Womack


I get an error message 'rms gradient is NaN'

This rather unhelpful message will be improved in the next release. For the August 2009 academic release, it normally means that your PDB file has two atoms in exactly the same place, generally because you went through a list of peaks in coot and inserted a water atom in one of them twice.

To find these, you could do

egrep -e "^(ATOM  |HETATM)" model.pdb | awk '{print substr($_,30,24)}' | sort | uniq -d

which will tell you the coordinates, and then grep model.pdb for the coordinates and delete one of the coincident atoms.


I get a 'error code 11' from grade when using mol2 input generated by CORINA

This is a known issue in grade in the February release; adding

-o t=mol2,fcharges

to your CORINA command line should produce mol2 output that grade can parse. This will be fixed in the next release.


I get a 'server: bind: Operation not supported' message from refine

BUSTER requires that the file system it's running on supports the creation of UNIX sockets - the Linux and Darwin local filesystems and NFS do this, but AFS does not. If you see this message, try running BUSTER in a directory on a local filesystem.


How can I check my BUSTER is correctly installed? How fast is it supposed to run?

Answer from Tom Womack

We distribute a test script refine_test.sh which runs a small refinement, using the options that we would recommend for regularising a structure and generating a new map after you've done a bit of rebuilding in coot; information on how to run it and on how long it takes to run on several machines we have tested it on is at the BusterShortRefineTest page.


What BUSTER output should I look at?

Please see dedicated page on interpreting BUSTER output


When should I start using TLS?

Answer from Tom Womack:

For relatively high-resolution X-ray data, I've found that TLS is useful for making difference density clearer at the end of a refinement; once modelling to take account of difference density is starting to become frustrating, running BUSTER with -M TLSbasic sometimes makes the density enough more clear that you can do another couple of model-and-run-BUSTER cycles. Once you've started using TLS, it's probably best to keep using it for the rest of the refinement.

For low-resolution data and models where there is a well-understood domain structure, TLS may well be useful from the very beginning, and we'd be very interested in hearing any success stories, particularly ones non-proprietary enough that they could at some point be converted into tutorial material.


Why does BUSTER use one B-factor per atom even at 3Å resolution?

Answer by Tom Womack

The BUSTER geometry function incorporates a very strong BCORREL term which couples the B-factors of bonded atoms; so neighbouring atoms are restrained to having similar B-factors. In tests, this 'soft grouping' works well from high resolution at least as far as 3Å.

At very high resolutions (say 1.2Å and below) you may want to use GeometryWeight_bcorrel=0 to allow the B-factor to be refined separately and in an unrestrained fashion for every atom.


In visualise-geometry-coot, what is an "Unhappy atom"?

Answer by Oliver Smart

The "unhappy atom" category uses the information from the geometry screen output (you can find this by follow link from LIST.html):

Atoms that make the largest contribution to overall function value.
    Report the (function contribution)/(median contribution). Median= 8.093
        5.5     A|175:CA(LEU)
        4.8     D|175:CA(LEU)
 ...

So it is a measure of the atom's total geometry function contribution including BOND, ANGLE, STORS, PLANE, BCORREL, IDEAL, BAD terms. LSSR is excluded from the measure. So if an atom is "unhappy" this can be because a combination of things add to make it so. Generally if the number is large the area will need rebuilding.


Why does BUSTER refinement result in a molprobity score with a large number of C-beta deviations? Is there anything that I can do about this?

Answer by Oliver Smart

-Gelly protgeo_option_chiralrestraint_from_equilib.dat

Why do the Rfree and Rfact values change abruptly at the start of the final cycle of a BUSTER refinement?

Answer by Andrew Sharff

If the parameter AnalyseVoids is set to 'yes', as it is by default, then the solvent-correction process in the final cycle attempts to account for hydrophobic voids. We have found that this is often useful in removing patches of negative difference density from, for example, the interior of helix bundles; if you find that the change is for the worse on one of your structures, try running with refine AnalyseVoids=no.


How should I run a BUSTER refinement after some model rebuilding with coot?

Answer by Andrew Sharff

The February 2011 release includes macros specially for BUSTER refinement after some model rebuilding with coot: see AutoBusterShortRunMacros page.


What is an easy way to get started with occupancy refinement?

We provide a simple tool to generate appropriate commands (in Gelly-syntax) describing the most common classes of occupancy refinement:

Analysis is based on some basic assumptions, like

Running

% pdb2occ -p your.pdb -o occref.gelly
% refine -p your.pdb -Gelly occref.gelly ...

will first analyse the existing PDB file (creating the text file occref.gelly), which is subsequently used in refinement. The text file can also be used as a starting template for defining more complicated situations. For further information and a full example please see AutoBusterLigandAndOccRef.


How do I specify the value used to select the set of reflections used for computing R-free?

Answer by Claus Flensburg

The October 2010 release introduced the parameter BusterFreeFlagValue:

% refine ... BusterFreeFlagValue=0

It defaults to '0' to follow the usual CCP4 selection. For releases prior to that one has to use the indirect specification:

% refine ... BusterExtraKwds="FREFLG=0"

You can also see AutoBusterFreeRFlag


What exactly do the 24 columns in the output MTZ file mean?

Answer by Tom Womack

There are a few circumstances (PDB deposition in particular) in which you have to pick columns from the MTZ file; the MTZcolumns page aims to explain what each column means and under what circumstances you should pick it.


Handling LINKs and covalent ligands in BUSTER

Answer by Oliver Smart


How do I set up a disulfide bond to a symmetry-related CYS?

Answer by Oliver Smart

This can be done by adding a SSBOND card before the CRYST1 line in the input PDB file. The card must have correct column spacing.

SSBOND   1 CYS A  429    CYS A  476                          1555   4445  2.04
key:           c nnnn        d mmmm                               oooooo
change 'c' and 'nnnn' to the chain id residue number of the identity copy CYS
       'd' and 'mmmm' to the chain id residue number of the symmetry related CYS
leave  'oooooo' alone, provide it is not 1555 BUSTER will work out the correct symmetry operator.

PDB entry http://www.rcsb.org/pdb/files/4J0U.pdb provides an example with correct column spacing and a pretty disulfide across a symmetry contact:

Currently BUSTER refine will suggest adding a SSBOND lacking a symmetry operator that will not work! This bug will be fixed in the next release this note dated April 2013


Is there anything I can do with grade if all I have as input is a PDB file with no hydrogens?

If that's all you have, then the general answer is no; it's entirely possible for a PDB file arising from a sufficiently badly driven refinement process to lack critical information about the chemistry of the ligand. But if you also have a good idea as to what the chemistry is, have a look at HydrogenatePDBwithOpenBabel


What is the order and the priorities for setting specific options?

The order in which settings can be changed from the defaults when running BUSTER:

1. installation-wide

This means a .autoBUSTER file in the same directory as the command (e.g. 'refine'). Or in more practical terms, these files should be located in

    % dir=`which refine`
    % ls -l `dirname $dir`/.autoBUSTER
    % set dir=`which refine`
    % ls -l `dirname $dir`/.autoBUSTER

2. user-wide

A file in the users home directory, i.e.

$HOME/.autoBUSTER

3. directory-wide

A file in the current directory (where the command, eg. 'refine' is being issued):

./.autoBUSTER

4. environment

Setting environemental variables in the usual shell-dependent way, eg.

    % MyParameter="XYZ"
    % export MyParameter
    % refine ...
    % setenv MyParameter "XYZ"
    % refine ...              

5. command-line

% refine MyParameter="XYZ" ...

or

% refine -M MyMacro ...

where a file MyMacro.macro contains something like

# my fancy stuff
MyParameter="XYZ"

and the directory containing that MyMacro.macro file being specified in the autoBUSTER_MacroDirs environmental variable.

Remember, the command-line is read in order, so

% refine MyParameter="ABC" -M MyMacro ...

will have a different effect to

% refine -M MyMacro MyParameter="ABC" ...

There is a big advantage of just using command-line arguments of the form par="value": the complete command-line is written into the header of the final refine.pdb file - so you always know exactly how it was run.


What is the optimal number of processors to use with BUSTER?

Please see the dedicated page Optimal number of processors for BUSTER.


BusterFAQ page. Please address problems, corrections and clarifications to buster-develop@globalphasing.com