Author(s)

• dock_in_map: Tom Terwilliger

Purpose

The routine dock_in_map will automatically dock a model or models into a map.

Usage

How dock_in_map works:

Dock-in-map uses a convolution-based shape search to find a part of a map that is similar to the model. The key elements of the search are:

Initial search at lower resolution that focuses on overall shape of the
molecule.  This allows a quick search but is less effective if the
map contains more than one molecule

Initial search without rotation.  This allows a fast search that can place
a molecule that is just shifted by a translation

Optional initial search based on matching moments of inertia of model and map.
This is also a fast search but requires an accurate map that looks a lot
like the model and does not have extra density.

Full search at the full resolution of the map.  This search can be run on
multiple processors to speed it up.

Examples

Standard run of dock_in_map:

Running dock_in_map is easy. From the command-line you can type:

phenix.dock_in_map 1ss8_A.pdb emd_8750.map resolution=4 nproc=4 \
   pdb_out=placed_model_from_emd_8750.pdb

where 1ss8_A.pdb is the model you would like to place, emd_8750.map is a CCP4, mrc or other related map format, and you specify the nominal resolution of the map and the number of processors to use.

Possible Problems

Specific limitations and problems:

Literature

Additional information

List of all available keywords

• job_title = None Job title in PHENIX GUI, not used on command line
• input_files
  • map_file = None File with CCP4-style map.
  • seq_file = None Optional sequence file
  • search_model = None Input PDB file(s) to be placed in the map.
  • search_model_copies = None You can specify how many copies for each search model by supplying one value of search_model_copies for each of your search models.
  • fixed_model = None Fixed model marking region not to be considered. Optionally to be added to placed search model to create new final model if append_to_fixed_model=True
  • ncs_file = None Symmetry file (.ncs_spec format or MATRX records) with reconstruction symmetry. Used in identification of unique part of map.
• output_files
  • pdb_out = placed_model.pdb Output PDB file matching the map.
• crystal_info
  • resolution = None High-resolution limit for main search. This can be lower resolution than the data. The search is quicker at lower resolution. .short_caption = High-resolution limit
  • scattering_table = n_gaussian wk1995 it1992 *electron neutron Choice of scattering table for structure factor calculations. Standard for X-ray is n_gaussian, for cryoEM is electron.
  • sequence = None Sequence as text string. Normally supply a sequence file instead
• search
  • low_res_search = True Try to fit by searching at low resolution first
  • resolution_radius_scale = 0.5 Resolution for low-res search will be resolution_radius_scale times the radius of gyration of the search model.
  • align_moments = False Try to fit by aligning moments of inertia if size of density region and molecule are similar
  • max_radius_ratio = 2. Radius of gyration of molecule and density must be within max_radius_ratio of each other
  • radius_scale = 1.5 Mask for density and atoms will be radius_scale times the radius of gyration of model
  • rigid_body_refinement = True Run rigid-body refinement on final model
  • rigid_body_refinement_single_unit = None Run rigid-body refinement with just one unit (do not break up into chains)
  • rigid_body_refinement_resolution = None Run rigid-body refinement at this resolution if specified
  • append_to_fixed_model = True Append input model to fixed model (if any) after search
  • min_cc = 0.4 If quick run, stop if minimum CC is achieved in local search.
  • run_in_boxes = True Run on sub-boxes and combine at the end
  • target_boxes = None Try to get this many boxes. Default is same as nproc.
  • box_to_run = None Run only this box
  • box_overlap_scale = 1 box overlap (overlap of boxes) will be box_overlap_scale times the density radius
  • edge_ratio = 10 box edge box_overlap times edge_ratio
  • density_radius = None Radius for density to be cut out and compared. Default is 6 times the resolution.
  • model_radius = 3 Radius for removing density near fixed_model
  • zero_value = 0 Value to set map in regions overlapping fixed model
  • density_peaks = 20 Number of NCS-related peaks of density to check
  • delta_phi = 20 Angular spacing of search
  • max_rot = None Maximum rotations to try
  • single_positions_to_try = 10 Number of offset positions to try in optimizing orientation. Positions along the chain are selected as centers and a local fit near each position is carried out. The resulting offsets relative to the original placement are used to optimize the overall orientation and position.
  • max_position_shift_frac = 0.05 Maximum fractional positional shift in single_positions run
  • min_relative_cc = 0.67 Minimum local CC relative to original CC to keep a local search. This is a way to reject local searches that are completely wrong.
  • sieve_fit = None Use sieve_fit fraction of single positions in fitting. If None, use all
  • ncs_copies_max = None Maximum number of matching models to write. If more than one they will be written as MODEL records in the output PDB file. You can get them individually afterwards with phenix.pdbtools placed_model.pdb keep="model 1" etc.
  • start_rot = None Three numbers rotx, rotz, rotx defining the starting rotation of the search model. Normally used along with delta_phi=1000 to generate exactly one defined rotation.
  • search_center = None Optional coordinates in search model for centering search. Note this is different from search_near_xyz which is the location xyz in the map to look.
  • search_center_selection = None Optional selection defining coordinates in search model for centering search.
  • rot_id_n = None Number of rotation groups. Along with rot_id_group, allows defining groups of rotations to be carried out in one run. .short_caption = Number of rotation groups
  • rot_id_group = None rotation group to include. See rot_id_n.
  • map_box = True Run map_box to extract useful part of map before search
  • fix_search_position = False You can choose to not move your model center of mass to the origin by fixing the search position
  • search_box_size = None You can choose the size of the search box for local searches. Normally set by default corresponding to 3 times density_radius
  • lower_bounds = None You can select a part of your map for analysis with lower_bounds and upper_bounds.
  • upper_bounds = None You can select a part of your map for analysis with lower_bounds and upper_bounds.
  • keep_search_order = None Keep search order as input
  • remove_water = True Remove waters from input files
• control
  • nproc = 1 Number of processors to use
  • random_seed = 171731 Random seed
  • verbose = False Verbose output
  • quick = True Try translation search first and stop if CC is at least min_cc
• guiGUI-specific parameter required for output directory
  • output_dir = None