26. Phenix Search Fit protocol
Phenix Search and Fit a sequence in the map densityProtocol designed to fit in Scipion a small sequence for which we were able to trace the skeleton of \alpha carbons in the experimental map. The process finishes with the refinement in the real space of the structure candidates. These potential fitted structures will be shown ranked by its map-model fitting score. For more information about the real space refinement have a look to Appendix Phenix Real space refine [Afonine et al., 2018].
- Requirements to run this protocol and visualize results:
- Scipion plugin: scipion-em
- Scipion plugin: scipion-em-phenix
- PHENIX software suite (tested for versions 1.18.2-3874, 1.19.2-4158 and 1.20.1-4487)
- Scipion plugin: scipion-em-ccp4
- CCP4 software suite (from version 7.0.056 to 7.1)
- Scipion plugin: scipion-em-chimera
- Scipion menu: Model building -> Flexible fitting (Fig. 26.1 (A))
- Protocol form parameters (Fig. 26.1 (B)):
NOTE:
Since the process of fitting several structures in the density map takes quite a long time we suggest to increase the number of MPI to accelerate the execution of the protocol.- Input volume: Electron density map previously downloaded or generated in Scipion to fit the atomic structure.
- Resolution (Å): Electron density map resolution.
- Input atom structure: Atomic structure previously downloaded or generated in Scipion to be fitted to an electron density map. Structures from PDB do not make sense here due to the purpose of the protocol of fitting partially traced models. Some residues of this model could be traced as ALA.
- Extra Params: Advance param. First line of extra params that the user would like to add to the PHENIX real space refiment process that will be performed in the last step before retrieving the structure candidates of the atomic structure fitted to the density map.
- Test sequence: Sequence of the atomic structure that we want to fit in the density map. This sequence has previously imported in Scipion.
- Residues: Fill in this box with the residues of the *Test sequence that you are iterested in. The wizard on the right allows to delimit exactly the part of the sequence that has to be considered by selecting the first and last residues of this sequence (position and type of residue). Press the wizard and every residue of the sequence will be shown. Select the first one clicking it, then press
CTRL
, and select the second one by clicking it. The aminoacid sequence contained between those two selected first and last residues will be included in the box. - Extra Params: Advanced param. Second line of extra params that can be added to the Coot script. This script is generated in the mutation process that allows to retrieve all possible atomic structures of a certain length starting in the first residue selected and ending in the last one selected.
- Additional Advanced Params for real space refinement: Go to Appendix Real space refine.
Protocol execution:
Adding specific protocol label is recommended in Run name section, at the form top. To add the label, open the protocol form, press the pencil symbol at the right side of Run name box, complete the label in the new opened window, press OK and, finally, close the protocol. This label will be shown in the output summary content (see below). If you want to run again this protocol, do not forget to set to Restart the Run mode.Press the Execute red button at the form bottom.Visualization of protocol results:
After executing the protocol, press Analyze Results and the pop up viewer window will open. This window allows to observe the retrieved models as well as the map-model fitting scores for all structures (by default) or for a selection of them (Fig. 26.2).
- Show volume and atomic structures in ChimeraX: ChimeraX graphics window will be opened by default. Map and atomic structures are referred to the origin of coordinates in ChimeraX. To show the relative position of the density volume and the atomic structures, the three coordinate axes are represented; X axis (red), Y axis (yellow), and Z axis (blue) (Fig. 5.3). Coordinate axes, map, initial untraced and final fitted atomic structures are models numbers #1, #2, #3 and from #4 to #n, respectively, in ChimeraX Models panel. The different retrieved fitted models are ordered according to the map-model fitting score value.
- Max. Number Atom Structs.: To observe the total number of retrieved fitted models, maintain the default value (1000). Otherwise, write in the box the desired maximum number of models to see. We presume that the atomic structures selected will be the best ones since they have the highest map-model fitting scores. The selected number of structures will be open in ChimeraX instead the total number of them.
- Show Area around input atomic structure (Å): In order to make easier the observation of the atomic structures, a small fraction of the map density can be selected (3Å by default).
- Summary Plot: Plot of map-model fitting scores of the protocol retrieved models. The average and standard deviation values are also shown. The values have been computed considering the total number of models. When only a certain number of the models has been selected in the item Max. Number Atom Structs., the fitting score values of those selected models will be remarked in red.
- Summary content:
- Protocol output (below Scipion framework):phenix - search fit -> ouputAtomStruct 0;AtomStruct (pseudoatoms=False, volume=False).phenix - search fit -> ouputAtomStruct 1;AtomStruct (pseudoatoms=False, volume=False).phenix - search fit -> ouputAtomStruct 2;AtomStruct (pseudoatoms=False, volume=False).phenix - search fit -> ouputAtomStruct 3;AtomStruct (pseudoatoms=False, volume=False).phenix - search fit -> ouputAtomStruct 4;AtomStruct (pseudoatoms=False, volume=False).Pseudoatoms is set to True when the structure is made of pseudoatoms instead of atoms. Volume is set to True when an electron density map is associated to the atomic structure.
- SUMMARY box:No summary information