Tomographic dual tilt series alignment and reconstruction

Tomo dual flow chart

The aim in aligning a dual tilt series is similar to that of two single tilt series. The difference is that the numbering of the fiducial markers in two series should correspond. The program bdual provides a way to transfer the seed markers from the first to the second series after fitting the zero-tilt micrographs to each other. After tracking the markers and aligning the micrographs, the 3D marker models are fit to each other and the transformations for the second series of micrographs adjusted to correspond to the first, again with the program bdual. A single reconstruction is then calculated from both series.

1. Setup and seed markers

The two tilt series should be set up separtely with the normalized micrographs:

btomo -verb 7 -axis 90 -tilt -60,2 -gold 7 -out ser1.star ser1.pif

btomo -verb 7 -axis 90 -tilt -60,2 -gold 7 -out ser2.star ser2.pif

A fiducial marker seed is then selected for the zero-tilt micrograph of the first series in bshow, and saved to a new parameter file (ser1_seed.star). The marker seed is transfered to the second series:

bdual -v 1 -seed 1,-100,-80 -refine -out ser2_seed.star ser1_seed.star ser2.star

The first argument of the -seed option sets the angular increment (in degrees) to search for the rotation angle that will superimpose the first zero-tilt micrograph onto the second zero-tilt micrograph. The other two arguments set the range of angles to be tested (the default is 360 degrees). The -refine option refines the positions of the markers after the transfer. The markers for the second series are then checked in bshow to correspond to the fiducials.


2. Tracking the markers

The two series are tracked independently:

btrack -verb 1 -reset -track 5 -refine markers -resol 30,300 -out ser1_trk.star ser1_seed.star

btrack -verb 1 -reset -track 5 -refine markers -resol 30,300 -out ser2_trk.star ser2_seed.star

The -resetmodel options ensures that any adjustments in the seed marker positions are transferred to the 3D marker model.


3. Refining the alignment

Each alignment can be refined in bshow separately for the two tilt series (ser1_ref.star and ser2_ref.star). The parameter files are then combined with proper adjustments for the second series derived from a fit of the 3D marker models:

bdual -v 1 -transform -zcomp -bild dual_tx.bld -Post dual_tx.ps -out sets.star set1_ref.star set2_ref.star

The -bild output generates a BILD file with the two marker models that can be displayed in UCSF Chimera to assess the quality of the fit. The -Postscript option plots the the x, y and z coordinates of the fitted models to further assess the quality of the fit.

The combined alignment for both series can then be further refined in bshow. While it is dangerous to use the scale refinement too much, it may be necessary to use it one final time in the refinement to adjust for small scaling distortions between the two tilt series.


4. Dual tilt reconstruction

The reconstruction is done as for a single-tilt series, but using the parameter file referencing both tilt series:

btomrec -v 1 -rescale 0,1 -trans full -resol 30 -size 1024,1024,260 -rec dual_rec.pif -out dual_rec.star sets.star