SHARP User Manual previous next
Appendix 1

What can I do with results from SHARP?

Copyright    © 2001-2006 by Global Phasing Limited
 
  All rights reserved.
 
  This software is proprietary to and embodies the confidential technology of Global Phasing Limited (GPhL). Possession, use, duplication or dissemination of the software is authorised only pursuant to a valid written licence from GPhL.
 
Documentation    (2001-2006)  Clemens Vonrhein
 
Contact sharp-develop@GlobalPhasing.com


This describes the various files produced by SHARP and autoSHARP and how best to use them. It is especially important if you want to use these outside the graphical user interface.

Contents


resid.mtz

This output file from SHARP contains the columns to calculate log-likelihood gradient residual maps. The coefficients are stored as pairs of real/imaginary part of a "structure factor". They follow the same hierarchical syntax as the Sharp INput file: CC_XX_WW_BB_SRC_TYP_RI, where

   CC :  compound number
   XX :  crystal number
   WW :  wavelength number
   BB :  batch number
   SRC :  source of gradient or model information
   
  • ISO :  isomorphous/dispersive differences between this dataset and the reference dataset
  • ANO :  anomalous differences within this dataset
   TYP :  type of map
   
  • HAT :  heavy atom model map for this dataset
  • RES :  log-likelihood gradient residual map
   RI :  real or imaginary part
   
  • RL :  real part
  • IM :  imaginary part

ColumnTypeDescription
H, K, LMiller indices 
CC_XX_WW_BB_ISO_HAT__Frealamplitude of calculated heavy atom density for isomorphous/dispersive difference
CC_XX_WW_BB_ISO_HAT__Prealphase of calculated heavy atom density for isomorphous/dispersive difference
CC_XX_WW_BB_ISO_RES__Frealamplitude of residual map for isomorphous/dispersive difference
CC_XX_WW_BB_ISO_RES__Prealphase of residual map for isomorphous/dispersive difference
CC_XX_WW_BB_ANO_HAT__Frealamplitude of calculated heavy atom density for anomalous difference
CC_XX_WW_BB_ANO_HAT__Prealphase of calculated heavy atom density for anomalous difference
CC_XX_WW_BB_ANO_RES__Frealamplitude of residual map for anomalous difference
CC_XX_WW_BB_ANO_RES__Prealphase of residual map for anomalous difference

These columns can be used in a variety of ways:

  1. calculation of heavy atom maps: this is done with columns CC_XX_WW_BB_ISO_HAT__F/CC_XX_WW_BB_ISO_HAT__P or CC_XX_WW_BB_ANO_HAT__F/CC_XX_WW_BB_ANO_HAT__P. To do this in the CCP4 program FFT:
          fft hklin resid.mtz mapout 01_01_02_01_ISO_HAT.map <<end_ip
          LABI F1=01_01_02_01_ISO_HAT__F PHI=01_01_02_01_ISO_HAT__P
          end_ip
          
    This would calculate a heavy atom model map for wavelength 2 in the first crystal of the first compound (e.g. the second wavelength in a MAD experiment), based on dispersive differences.

  2. calculation of residual maps: this is done with columns CC_XX_WW_BB_ISO_RES__F/CC_XX_WW_BB_ISO_RES__P or CC_XX_WW_BB_ANO_RES__F/CC_XX_WW_BB_ANO_RES__P. To do this in the CCP4 program FFT:
          fft hklin resid.mtz mapout 01_01_01_01_ANO_RES.map <<end_ip
          LABI F1=01_01_01_01_ANO_RES__F PHI=01_01_01_01_ANO_RES__P
          end_ip
          
    This would calculate an anomalous log-likelihood gradient residual map for the first dataset/batch in the first crystal of the first compound (e.g. the anomalous data in a SAD experiment).


eden.mtz

The main output file from SHARP contains the following columns:

ColumnTypeDescription
H, K, LMiller indices 
FPamplitude"native", unperturbed structure factor amplitude
SIGFPstandard deviationstandard deviation of FP
HLA, HLB, HLC, HLDHendrickson-Lattman coefficientsdescription of phase probability
FBamplitudecentroid electron density structure factor amplitude
PHIBphasecentroid electron density structure factor phase
FOMfigure-of-merit 
NUMORB multiplicity description

These columns can be used in a variety of ways:

  1. calculation of electron density maps: this is done with amplitude FB and phase PHIB (no need to include a weight like FOM). To do this in the CCP4 program FFT:
          fft hklin eden.mtz mapout eden.map <<end_ip
          LABI F1=FB PHI=PHIB
          end_ip
          

  2. density modification: the starting map should always be calculated using FB/PHIB. The native structure factor amplitudes (and its standard deviation) are FP/SIGFP. For phase combination, the Hendrickson-Lattman coefficients HLA, HLB, HLC and HLD are used. A (typical) DM run could look like this:
          dm hklin eden.mtz hklout dm.mtz <<end_ip
          SOLC 0.50
          LABI FP=FP SIGFP=SIGFP PHIO=PHIB FOMO=FOM -
               HLA=HLA HLB=HLB HLC=HLC HLD=HLD -
               FDM=FB PHIDM=PHIB
          LABO FDM=FDM PHIDM=PHIDM
          end_ip
          
    (please check your documentation for complete syntax)

    Note : the phase improvement and interpretation tools provided through the interface (e.g. using SOLOMON for solvent flattening) are highly optimised for use with SHARP. We strongly recommend the use of these for phase improvement after phasing with SHARP. However, if you find a program/strategy that performs better in your particular case, we would be glad to hear about it.


eden-unique.mtz

This is the same as eden.mtz but an additional column (FreeR_flag) has been added to distinguish test and working sets (Brünger, 1992). The various phase improvement and interpretation tools available through the interface are using this file (so that the FreeR_flag is always carried through).

eden_flat_<solvent fraction>pc.mtz

Any solvent flattening run will produce a final MTZ file eden_flat_<solvent fraction>pc.mtz (where <solvent fraction> is the solvent fraction used for this run - which also distinguishes different solvent flattening results). This file contains the following columns:

ColumnTypeDescription
H, K, LMiller indices 
FPshaamplitude"native", same as FP in eden.mtz
SIGFPshastandard deviationstandard deviation of FPsha (same as SIGFP in eden.mtz)
HLA, HLB, HLC, HLDHendrickson-Lattman coefficientsdescription of phase probability (same as in eden.mtz)
FBshaamplitudecentroid electron density structure factor (including average heavy atom contribution) amplitude (same as FB in eden.mtz)
PHIBshaphasecentroid electron density structure factor (including average heavy atom contribution) phase (same as PHIB in eden.mtz)
FOMshafigure-of-merit(same as FOM in eden.mtz)
FreeR_flagtest flag 
FBshasolamplitudesolvent flattened structure factor (including average heavy atom contribution) amplitude
PHIBshasolphasesolvent flattened structure factor (including average heavy atom contribution) phase
Fcentshasolamplitudesolvent flattened structure factor (light atoms only) amplitude
PHIcentshasolphasesolvent flattened structure factor (light atoms only) phase
PHIshasolphasecombined phase
FOMshasolweightcombined figure-of-merit
HLAshasol, HLBshasol, HLCshasol, HLDshasolHendrickson-Lattman coefficientsdescription of phase probability (after density modification with SOLOMON)

These columns can be used to calculate solvent flattened electron density maps, e.g. with the CCP4 program FFT:

fft hklin eden_flat_50.0pc.mtz mapout eden_flat_50.0pc.map <<end_ip
LABI F1=FBshasol PHI=PHIBshasol
end_ip

wARP_<solvent fraction>pc.mtz

If some automatic building was done using ARP/wARP, a final MTZ file is available as wARP_<solvent fraction>pc.mtz (where <solvent fraction> is the solvent fraction of the solvent flattening run that this automatic interpretation was started from). The resulting file contains the following columns:

ColumnTypeDescription
H, K, LMiller indices 
FPshaamplitude"native", unperturbed structure factor amplitude (same as FP in eden.mtz)
SIGFPshastandard deviationstandard deviation of FPsha (same as SIGFP in eden.mtz)
FreeR_flagtest flag 
FCamplitudecalculated structure factor amplitude of (auto-built) model
PHICphasecalculated structure factor phase of (auto-built) model
2FOFCWTamplitudeFourier amplitude for (2mFo-DFCalc) map
PH2FOFCWTphaseFourier phase for (2mFo-DFCalc) map
FOFCWTamplitudeFourier amplitude for (mFo-DFCalc) map
PHFOFCWTphaseFourier phase for (mFo-DFCalc) map
FOMfigure-of-merit 
PHCOMBphasecombined (with experimental phases from SHARP) phase
(see also your local documentation for REFMAC)

The best electron density map can be calculated with FFT:

fft hklin wARP_50.0pc.mtz mapout wARP_50.0pc.map <<end_ip
LABI F1=2FOFCWT PHI=PH2FOFCWT
end_ip

model_<project>.<id>_wARP_<solvent content>pc.pdb

After automatic building with ARP/wARP, the final model PDB file is (linked) in the sharpfiles/datafiles directory. The naming convention is model_<project>.<id>_wARP_<solvent content>pc.pdb, where
Last modification: 28.04.2015