Chimera Hands-On Tutorial: Fitting a ParM Filament Map

Presenter:Tom Goddard, Chimera developer
Date:Tuesday, November 30, 2011
Time:10:15 - 11:45 AM
Room:MRGN 129, Purdue University

Tutorial Setup

Bring your laptop with the following Chimera version and data files.

Older Chimera versions won't work. They do not have many new features we will use.


Open Question and Answer Session

30 minutes will be reserved for answering any Chimera questions with live demonstration. Questions need not be related to the tutorial.

What does ParM do?

Modeling a ParM Filament

Analysis Steps

Show Molecule and Map Fit Molecule in Map
  1. Open PDB 1mwk, EMDB 5129.
  2. Show Command Line (Favorites menu).
  3. Deactivate map (model 0) below command-line to move 1mwk away from map.
  4. Command "rainbow chain" to color two ParM monomers.
  5. Delete chain B. Use menu Select / Chains / B, then menu Actions / Atoms / Delete.
  1. Fit 1mwk in map.
  2. Move 1mwk into map.
  3. Press Fit button in Fit in Map dialog (volume dialog Tools menu).
  4. Fit using correlation: Fit dialog Options button, enable "Use map simulated from atoms..." resolution 19. Press Fit.
  5. Show simulated map (Volume dialog eye icon) as mesh.
  6. Correlation depends on domain of calculation. Change simulated map threshold and press Update in Fit dialog.
  7. Spend a few minutes trying alternative fit positions.
Global Fit Search Calculate Map Symmetry Fit Asymmetric Unit
  1. Search for best fit using 30 tries with command "fit #1 #0 search 30"
  2. Fits appear all along filament, many are equivalent due to symmetry of the filament.
  1. Determine map symmetry to eliminate equivalent fits.
  2. Command "measure symmetry #0 helix 20,180,opt minimumCorrelation 0.95".
  3. The "helix" option gives Chimera a hint about helical parameters.
  4. The "minimumCorrelation" option accounts for this unusual map where the helix does not extend to the edges of the volume box.
  5. View symmetry copies of molecule with command "sym #1 group #0 surf true".
  6. Remove symmetry copies with "~sym #1" (note leading tilde character which means "undo" in Chimera commands).
  1. Clear fit list.
  2. Rerun previous fitting command "fit #1 #0 search 50".
  3. Clear fit list.
  4. Use correlation optimization "fit #1 #0 search 50 res 19".
Make Filament Model Reduce Molecular Clashes Calculate Predicted Map
  1. Show symmetric molecule copies for best fit. Command "sym #1 group #0 update true".
  2. Color molecules distinctly. Command "rainbow model".
  1. Inspect clashes between adjacent ParM molecules.
  2. Fit asymmetric unit including all overlapped symmetric molecules. Command "fit #1 #0 sym true res 19"
  3. Note increased space between molecules.
  1. Compute difference map between experimental map and predicted map for molecular model.
  2. First delete extra ParM molecules outside experimental map. Ctrl-drag mouse to select outside molecules. Press up-arrow key to extend selection to full molecules. Menu Actions / Atoms / Delete to delete.
  3. Calculate predicted map. Command "molmap #1,2, 19"
Difference Map Compare "Open" and "Closed" Filaments.
  1. Subtract two maps, scaling the second to minimize difference. Command "vop subtract #0 #3 minRMS true"
  2. Adjust difference map contour level. Add negative contour with ctrl-click on histogram. Adjust contour colors.
  3. Hide molecular models with Model Panel (menu Favorites).
  1. Load "closed" filament (GDP bound), EMDB 5128, menu File / Open...
  2. Flip closed map 180 degrees and fit to open map.
  3. Compute closed map on same grid as open map. Command "vop resample #5 onGrid #0"
  4. Use Morph Map (Tools menu of volume dialog) to morph between maps. Movie.