Data Processing after WS-MR

1 Verify Results for Search Models

The first follow-up step for a search model resulting from WSMR is to repeat MR locally.

To download one of the search models; first click the "Load Top 100 Hits" button, click one of the yellow circles in the resulting plot, and then use the download link for the model.




A typical input deck (MR01.in for this example) for phaser would look something like this (with the paths to your input MTZ and PDB files adjusted appropriately):

ROOT            MR01
MODE            MR_AUTO
HKLIn           reflections.mtz
LABIn           F=FP SIGF=SIGFP
ENSEmble        SearchModel PDB model.pdb IDENtity 0.3
COMPosition     PERCENTAGE 0.50
RESOlution      HIGH 3.1
SEARch          ENSEmble SearchModel NUM 1
SGALternative   SELECT ALL

This can be run with versions of phaser distributed with CCP4 or PHENIX (with phaser < MR01.in > MR01.out or phenix.phaser < MR01.in > MR01.out respectively). The statistics (LLG and Z-scores) from local runs may differ due to differences in phaser versions.

The SGALternative card is only necessary if you're checking alternative space-groups.

2 Expand Search

WSMR searches for a single copy of the search model - however there may be multiple copies of the model in the ASU. If so, the next step is the increase the NUM value for the SEARch card. When in doubt, it's better to set this too hight rather than too low - phaser will stop when unable to locate additional copies.

3 Autobuild

Once MR has been exhausted (no more copies of the search model can be placed), one possible next step is to use an autobuilder for an initial attempt to correct the sequence of the model.

Using the PHENIX autobuilder:

phenix.autobuild reflections.mtz MR_Results.pdb project.fasta nproc=4 > AutoBuild01.out

The input model (MR_Results.pdb) for the autobuilder is the final result modeles from the initial or expanded MR search.

The input MTZ file should be the original amplitudes (or intensities) used for searching, not the results from phaser. If phaser was run with the SGALternative card, make sure the spacegroup of this MTZ matches that of the search solutions.

The nproc=4 option should only be used if there are four (or more) available cores - when using the rebuild_in_place option, phenix.autobuild can't take advantage of more than 4 processors.