Back to resolve table of contents

Running resolve


Summary

To run resolve, you need:

setenv SYMOP /usr/local/lib/ccp4/symop.lib
setenv CCP4_OPEN UNKNOWN

The most basic resolve script  (see the sample scripts for more cases):
 

#!/bin/csh
#
# Here is a minimal script to run resolve on MAD/MIR/SAD etc data:
#
# Set CCP4 variables for symmetry information and
# for file handling:
#
setenv SYMOP /usr/local/lib/ccp4/symop.lib
setenv CCP4_OPEN UNKNOWN
#
# Now run resolve:
#
resolve<<EOD
hklin solve.mtz
LABIN FP=FP PHIB=PHIB FOM=FOM HLA=HLA HLB=HLB HLC=HLC HLD=HLD
hklout resolve.mtz
solvent_content 0.4             ! your solvent content goes here.
EOD
#
# Now "resolve.mtz" has the output amplitudes, phases,
# and figure of merit in columns labelled: FP PHIM FOMM.  
#

Notes on input and output mtz files (see the sample scripts )

Keywords for resolve

KEYWORD         DEFAULT                 WHAT IT IS

access_file     solve2.access            Name of solve2.access file.  If it is not in the
                                              /usr/local/lib/solve/ directory or in the
                                              current directory or in the directory $SOLVEDIR
                                              then you will want to tell resolve where it is 

hklin           solve.mtz               mtz file with input phases and phase probabilities

hklout          resolve.mtz             mtz file with output phases

LABIN FP=FP PHIB=PHIB FOM=FOM HLA=HLA HLB=HLB HLC=HLC HLD=HLD
                                        LABIN statement identifying the columns of 
                                        data in the hklin mtz file
LABIN FP=FP FC=FC PHIC=PHIC FOM=WCMB  FWT=FWT    LABIN statement suitable for prime-and-switch. 
                                        FWT is optional; if present initial map calculated with FWT exp(i PHIC).


mask_cycles     5                       Number of cycles in which masks are redone and
                                        images are compared to map

minor_cycles    10                      Number of minor cycles per mask_cycle


prime_and_switch                        Use prime-and-switch phasing. Use input phases only to calculate initial map. 
                                        Input phase probabilities not used at all.

no_erase_protein                        Do not let P(protein) be less than the starting value in prime-and-switch phasing.

solvent_content 0.3                     Fraction of unit cell in solvent region
resolution                              Resolution limits (default=whatever is in input mtz file)

wang_radius                             (default varies with resolution and FOM)
                                        Radius (A) for averaging rms density in Wang method for getting solvent mask
hklstart                                mtz file with a starting set of phases. Resolve will start with
                                        these phases (but use probabilities from hklin).  Useful for
                                        running a few cycles, getting an output resolve_1.mtz, then
                                        continuing on from there.  Goes with labstart. 

labstart  FP=FP PHIB=PHIB FOM=FOM       LABIN statement for hklstart.  only PHIB and FOM used.
cc_ratio      0.1                       Cut off prime-and-switch phasing at resolution where CC 
                                        of FP with FC < cc_ratio
verbose                                 print out a lot of data every cycle, not just at first and last.
nohl                                    don't calculate HL coefficients at end of resolve
rota_matrix                             rotation matrix for NCS symmetry. Three rota_matrix lines define a matrix:
                                        rEnter the identity as first molecule 1.
                                        rota_matrix 1 0 0
                                        rota_matrix 0 1 0
                                        rota_matrix 0 0 1
                                        The rotation matrix applies to orthogonal Angstrom coordinates.
                                        The matrix and translation maps molecule j on to molecule 1.
                                        This is what you get from the ccp4 program lsqkab if molecule 1 is
                                        the reference (xyzin1) and molecule j is the working molecule (xyzin2)
tran_orth                               Translation vector for NCS symmetry element, in orthogonal Angstroms
                                        Enter after rota_matrix
center_orth                             Approximate center of mass of NCS symmetry element, in orthogonal Angstroms
                                        Enter after rota_matrix. Required only for element 1
fraction_ncs  0.15                      fraction of asymmetric unit occupied by 1 molecule (default = fraction protein / N )
invert                                  Invert NCS matrices (they correspond to mapping molecule 1 on to molecule j)
ncs_begin  mask_cycle small_cycle       Begin using NCS on mask cycle mask_cycle, small cycle small_cycle
                                        Resolve works best if the solvent-flattening is mostly done before starting NCS
                                        Default is to start halfway through first mask_cycle if solvent content <50%, halfway
                                          through second mask_cycle if solvent content is 50% or greater.
fix_ncs                                 Do not refine NCS operators
force_ncs                               Use input or heavy-atom-site NCS symmetry even if there is very low correlation
ha_file    ha.pdb                       Use the entries in this file (PDB format) to look for non-crystallographic symmetry
no_find_ncs                             Don't try to find NCS from heavy-atom sites even if the file "ha.pdb" exists