Python-based Hierarchical ENvironment for Integrated Xtallography |
Documentation Home |
Replace side chains in a model to match the sequence with replace_side_chains
Author(s)
Purposereplace_side_chains is a tool to take a main-chain model and to add side-chains from a rotamer library based on an electron density map. NOTE: Normally this method is called from phenix.phase_and_build but you can use it interactively if you want. UsageHow replace_side_chains works:replace_side_chains uses resolve to replace the side chains in a model. You can choose to use the default resolve rotamer library or the Lovell et al penultimate rotamer library. Output files from replace_side_chainsreplace_side_chains.pdb: A PDB file with side chains replaced. ExamplesStandard run of replace_side_chains:Running replace_side_chains is easy. From the command-line you can type: phenix.replace_side_chains model.pdb map_coeffs.mtz sequence.datYou should supply a sequence file that matches the sequence of the model. Possible ProblemsSpecific limitations and problems:LiteratureAdditional informationList of all replace_side_chains keywords------------------------------------------------------------------------------- Legend: black bold - scope names black - parameter names red - parameter values blue - parameter help blue bold - scope help Parameter values: * means selected parameter (where multiple choices are available) False is No True is Yes None means not provided, not predefined, or left up to the program "%3d" is a Python style formatting descriptor ------------------------------------------------------------------------------- replace_side_chains input_files seq_file= None File with 1-letter code sequence of molecule. Chains separated by blank line or greater-than sign pdb_in= None PDB file with coordinates to evaluate mtz_in= None MTZ file with coefficients for a map map_coeff_labels= None If map coefficients cannot be identified automatically from your MTZ file, you can specify the label or labels for them. (Please separate labels with blank space, MTZ columns grouped together separated by commas with no blanks.) You can specify: map_coeff_labels (e.g., FWT,PHIFWT) amplitudes and phases (e.g., FP,SIGFP PHIB) or amplitudes, phases, weights (e.g., FP,SIGFP PHIB FOM) labin= "" For backward compatibility only. Labin line for MTZ file with map coefficients. This is optional if replace_side_chains can guess the correct coefficients for FP PHI and FOM. Otherwise specify: LABIN FP=myFP PHIB=myPHI FOM=myFOM where myFP is your column label for FP output_files pdb_out= replace_side_chains.pdb Output PDB file log= None Output log file params_out= replace_side_chains_params.eff Parameters file to rerun replace_side_chains crystal_info resolution= 0. high-resolution limit for map calculation chain_type= *PROTEIN DNA RNA Chain type (for identifying main-chain and side-chain atoms) NOTE: residues that are not this type will be ignored and not written out. side_chains richardson_rotamers= False Use rotamers from Lovell et al., 2000 library force_replace= True Replace all side chains, whether better or not if False, then only replace if it improves fit directories temp_dir= "temp_dir" Optional temporary work directory control verbose= False Verbose output raise_sorry= False Raise sorry if problems debug= False Debugging output dry_run= False Just read in and check parameter names |