Project

General

Profile

Helical Image Processing (PHOELIX) » History » Version 32

Lauren Fisher, 09/01/2011 11:10 AM

1 1 Lauren Fisher
h1. Helical Image Processing (PHOELIX)
2
3 2 Lauren Fisher
h2. General Workflow:
4
5 9 Lauren Fisher
# Make sure that the appropriate run name and directory tree are specified.  Appion increments names automatically, but users are free to specify proprietary names and directories.
6 27 Lauren Fisher
# Enter the directory tree where the mandatory files are located.  These files include llbo.sa, strong.sa, range.sa, cutfit1.dek, cutfit2.dek, cutfit3.dek, chop1.dek, chop2.dek, and template. If you do not have these files already, complete the steps for [[Helical Image Processing (PHOELIX)#Running-PreHIP| Running preHIP]] and follow the guided protocol to generate them. 
7 2 Lauren Fisher
# Enter a description of your run into the description box.
8
# Select the stack to process from the drop down menu.  Note that stacks can be identified in this menu by stack name, stack ID, and that the number of particles, pixel and box sizes are listed for each.
9 16 Lauren Fisher
# Make sure that "Commit to Database" box is checked. (For test runs in which you do not wish to store results in the database this box can be unchecked.  This box is automatically deactivated when running preHIP).
10 2 Lauren Fisher
# Enter the Filament Parameters or verify the default values. 
11 17 Lauren Fisher
#* Number of Filament Segments - The number of particles to use in the run.
12
#* Repeat Length - (Helical Repeat) The 360 degree repeat distance along the filament, in Ångströms. Calculated in Fourier space by the distance to the first layer line or in real space by the distance it takes one subunit to travel 360°.
13
#* Subunits per Repeat - Number of subunits in one helical repeat. Determined by calculating (Helical Repeat/Rise) or finding the Layer Line with a J0 intercept.
14
#* Diameter - The approximate outer diameter of the filament, in Ångströms
15
#* Inner Diameter - The approximate inner diameter of the filament, in Ångströms.
16 2 Lauren Fisher
# Enter the Processing Parameters or verify the default values.
17 17 Lauren Fisher
#* Filament Segment Length - Desired length of the filament segments, usually a certain number of helical repeats, in pixels.  Should be the same as your stack box size.
18
#* Box Height - Height of the filament box after alignments, usually the power of 2 value greater than the filament diameter, in pixels.
19
#* Pad Value - Value to pad the filament box to during alignments, usually the power of 2 value greater than the filament length/box size.
20 18 Lauren Fisher
#* Phase Residual Cutoff - Maximum phase residual for post averaging.  A phase residual is calculated between the layer lines from the template and the layer line from each filament segment.  Any files with phase residual above this cutoff will not be included in the final average, map, overplot, etc.
21 17 Lauren Fisher
#* Filter Value - Lowpass filter value for filament segments used during alignments to find filament edges. 150-200 seems to work well, but should be tested using preHIP before processing the full stack
22
#* Binning - Binning of the particle images. This quickly shrinks the image to help make the processing faster.
23
#* Contrast Change - If your filament stack has black on white density, the final reconstruction will need a contrast change.
24 2 Lauren Fisher
# Click on "Run HIP" to submit your job to the cluster. Alternatively, click on "Just Show Command" to obtain a command that can be pasted into a UNIX shell.
25 6 Lauren Fisher
26 8 Lauren Fisher
!http://emg.nysbc.org/attachments/959/hipwiki_fig1.png!
27 3 Lauren Fisher
28 2 Lauren Fisher
h2. Running PreHIP:
29
30 24 Lauren Fisher
After completing steps 1-7 in the [[Helical Image Processing (PHOELIX)#General-Workflow|General Workflow]] complete these additional steps:
31 2 Lauren Fisher
# Click the "Run PreHIP" checkbox which will activate the indexing parameters below it and display a warning message.
32 17 Lauren Fisher
# Fill in the Indexing Parameters.
33
#* N-fold - Order of the symmetry axis. Default is 1-fold for no symmetry.
34
#* Maximum Layer Line - Maximum layer line to include in the llbo.sa file.
35
#* Maximum Bessel Order - Maximum Bessel Order to include in the computation for llbo.sa.
36 20 Lauren Fisher
# Select either the "Use rise and twist" or "Use layer line/bessel order" radio button which will activate the parameters below it. You can only choose one of these options, therefore the input parameter boxes under the non-selected option become deactivated.
37 17 Lauren Fisher
#* Rise and Twist - Use the rise distance and twist angle between protein subunits to determine selection rule and generate llbo.sa file. Rise is the vertical distance between two subunits in the one-start helix and twist is the angle between two adjacent subunits.
38
#* (1,0),(0,1) LLBO - Use the layer line and Bessel Order for the (1,0) and (0,1) layer lines to determine selection rule and generate llbo.sa file.
39 21 Lauren Fisher
# Click "Just Show Command" to generate the UNIX shell command.  Do NOT click "Run HIP", it will not run because the interactive GUI must be executed from a command line. 
40 14 Lauren Fisher
# Note that the "Commit to Database" box is deactivated. preHIP is used only for generating the mandatory files and therefore the preliminary results are not stored to the database.
41
!http://emg.nysbc.org/attachments/962/hipwiki_fig2.png!
42 1 Lauren Fisher
# Copy and paste the command into a UNIX shell.
43 29 Lauren Fisher
!http://emg.nysbc.org/attachments/967/hipwiki_fig3.png!
44 30 Lauren Fisher
# Review the llbo.sa file when it is displayed and either approve or reject the LLBO assignments. If you reject it, you have the option of either restarting preHIP with different rise/twist/helical repeat values or supplying your own llbo.sa (if you choose this option make sure it is in the proper format! You can use the failed llbo.sa as a reference.)
45 31 Lauren Fisher
!http://emg.nysbc.org/attachments/970/hipwiki_fig4.png!
46 32 Lauren Fisher
# Scroll through the files using the Tkir viewer to see if the filaments were centered properly (files with extension .s are raw files and .sbut are filtered).  When finished, close both Tkir windows.  Enter either "yes" or "no" after the prompt in the terminal.  If you enter "no", you will be prompted to enter a new filter value and the filaments will be recentered with this value.  Continue this process until you find a filter value that works. 
47
# All files are diffracted and summed together to find the strongest layer lines. Review the strong layer lines when they are displayed and either approve or reject them.  If they are not correct you have the option to change, remove, or add lines.  Follow the instructions as prompted until you are satisfied. 
48
!http://emg.nysbc.org/attachments/1003/hipwiki_fig5.png!
49
# The strong layer lines from step 9 are used to select ranges for out of plane tilt and shift correction. Review/adjust the ranges when they are displayed in the Tkll viewer. Right click to remove points, left click to add a point.  Each layer line can only have 2 points and they must be ordered 0-1. If you modify the file you must save it, do not change the filename. (File > Save Data > Ok)
50
!http://emg.nysbc.org/attachments/1004/hipwiki_fig6.png!