physiology

VMD plugins at WMC Physiology & Biophysics

RMSDTT: RMSD Trajectory Tool v2.0

Description | Options | Download | Installation | Author

Description

The RMSDTT plugin adds support for rmsd calculations and alignment of trajectories to the RMSD Tool plugin. Other feature added include:

  • Plotting data via the Multiplot VMD plugin (in all systems), Xmgrace (on unix systems) or MS Excel (on windows).
  • Saving data to a file.
  • Swapping equivalent atoms (i.e. Glu OE1 and OE2) to reduce RMSD.
  • Statistics (average, standard deviation, minima and maxima).
  • Multi-line atoms selection with support for comments.

The RMSDTT plugin is now distributed as part of VMD since version 1.8.4 (to be released shortly). Development snapshots of the RMSDTT plugin can be obtained here.

As an example, in the following figure the RMSD of the trace of minimized structure min and trajectories equilibration and md are calculated with the xray structure as reference (top molecule). Statistics for each molecule and totals are shown in the table. In addition a plot was requested (plot not shown). Note how lines can be commented in the atom selection box.

Options

  • Menu:
    • File: the RMSD values can be plotted or saved to file after the calculation (options Save data and Plot data) or during the calculation (through plot and save to file options in the main window). Do not change the trajectory options in the main window for this to work (except the Time options, which can be changed). Save summary will save the summary of the data as it appears in the table to file.
    • The Options menu allows you to select different programs for plotting, to change the definition of the backbone, switch on/off the output of statistics, and color the table by molecule number (the table will only be colored if this switch in on and the plot option is on during the calculation).
  • Atoms selection:
    • Multi-line atoms selection box. The selection can be split in different lines. The character '#' can be used to comment out a portion of the selection, namely from the '#' to the end of that line (see figure above).
    • trace will add 'and name CA' to the atoms selection.
    • backbone will add 'and name C CA N' (default) or 'and name C CA N O' to the atoms selection. The backbone definition to use can be choose from menu (Options | Backbone def...).
    • noh will exclude hydrogens from the atoms selection (adds 'and noh' to the atoms selection).
    • history is a drop down menu of previously used selections. It is updated when pressing RMSD (not anymore by pressing Return).
  • Swapping atoms:
    • Swap_atoms will try to swap equivalent atoms to reduce the RMSD value. Sets of equivalent atoms can be selected from the drop down menu (defaults to All defined pairs). How to define or delete new pairs is explained in the Swap plugin webpage. You need to have the Swap plugin installed to see this options. It does not work with Align.
    • The verbose option will print the out the status of the swapped atoms during the RMSD calculations.
    • The List button will print out the current definitions.
  • Calculation type:
    • RMSD will calculate the RMSD using the reference molecule defined in Ref and the structures in the table as targets. Note that the RMSD button has to be pressed again to obtain the new RMSD values after the alignment.
    • Align will align the structures in the table using the molecule defined in Ref as reference. Align is not available if the All in the trajectory options is active or Average is used as reference.
  • Reference molecule:
    • Top will use the actual top molecule as reference (it doesn't need to be one included in the table).
    • Average will use the average structure calculated from all the molecules in the table (and all frames if trajectory is activated) as reference. It does not work for Align calculations.
    • Selected will use the molecule selected in the table as reference. If more than one molecule is selected, only the first one will be used as reference.
  • Trajectory options:
    • Trajectory activates options to include the frames of a molecule in the calculations.
    • Frame ref sets the frame of reference for the reference molecule. if the reference molecule is set to Average, all frames are used.
    • All activates the calculation of the RMSD of all frames in the selected structures against all frames in the reference molecule.
    • Skip allows to skip frames in a regular fashion, starting from Start and skipping every Steps.
    • Time activates the use of time units. Ini is the initial time, and Step is the time between frames. This changes will be reflected in the plot and the output file.
    • Plot will use the Multiplot VMD plugin (available in all systems), Xmgrace (only on unix) or MS Excel (only windows) to plot the data RMSD vs frame/time. This options is not available if the All option is on. The same plot can be obtained from the File menu. Xmgrace or MS Excel have to be installed on the system by the user. MS Excel needs the tcom tcl library (see installation instructions). VMD gives a warning message when exiting on windows systems if MS Excel is open with data from VMD, but you still can continue using MS Excel.
    • Save to file will save the data to the specified file. The format of the data depends on the status of the All option. If it is off, the format is a table in which each row represents a frame and each column a molecule or trajectory. If All is on, the data is a table in which each row represents the RMSD between the corresponding ref_mol/ref_frame and target_mol/target_frame pairs. If the trajectories have differing numbers of frames, NA (not available) is used. The data can also be saved from the File menu.
      All off All on
      frame   mol 1   mol 2   mol 3
          0   0.652   2.828   0.652
          1      NA      NA   0.659
          2      NA      NA   0.674
          3      NA      NA   0.692
          4      NA      NA   0.722
          5      NA      NA   0.743
        ...     ...     ...     ...
      ref_mol ref_frame       mol     frame      rmsd
            0         0         1         0     0.652
            0         0         2         0     2.828
            0         0         3         0     0.652
            0         0         3         1     0.659
            0         0         3         2     0.674
            0         0         3         3     0.692
          ...       ...       ...       ...       ...
  • Table of target molecules: This table shows the results by molecule. It presents the average, standard deviation, minimum, maximum and number of frames used (this number does not need to be the same as the number of frames in the molecule, since the reference molecule/frame is not included in the calculation). The Overall row shows the overall values as if all frames were independent of each other, i. e., it is not an average of the averages.

    The following buttons can be used to add/remove molecules from the table:

    • Erase all deletes all molecules from the list.
    • Erase selected deletes the selected molecule from the list.
    • Add all adds all molecules to the list.
    • Add active adds only the active molecules to the list.

Download

Send a mail to the author if you would like to be informed about new updates of this plugin.
version 2.01 (24 Apr '06) [download]
Change log:
  • Using Average as the reference molecule was not working. Syntax of VMD command measure avpos updated (thanks to Christopher Lanci for reporting this bug).
version 2.0 (3 Apr '06) [download]
Change log:
  • Included in VMD 1.8.4.
  • Complete rewrite to speed up the calculations, add Multiplot support and redesign the gui.

Previous versions can be found here

Installation

The RMSDTT plugin v2.0 has been included in the VMD distribution since VMD version 1.8.4. If you have a previous VMD version I recommend you to upgrade. If you still need the RMSDTT plugin in a previous VMD version or if you want to update to a most recent version of the RMSDTT plugin, follow this instructions:

  1. Download the plugin from here to the directory where you store your VMD plugins.
  2. Add the following to your .vmdrc startup file (or create one) (for unix the path should be $HOME/.vmdrc, and for windows %USERPROFILE%\vmd.rc):
    1. For VMD 1.8.4 and up (will be used instead of the RMSDTT plugin version distributed with VMD 1.8.4):
      set auto_path [linsert $auto_path 0 {/path/to/your/plugins/directory}]
      
    2. For VMD 1.8.3
      lappend auto_path {/path/to/your/plugins/directory}
      vmd_install_extension rmsdtt rmsdtt_tk_cb "WMC PhysBio/RMSDTT"
      
    3. For VMD 1.8.2 and down
      lappend auto_path {/path/to/your/plugins/directory}
      if { [catch "package require rmsdtt" msg] } {
        puts "VMD RMSDTT package could not be loaded:\n$msg"
      } elseif { [catch {menu tk register "rmsdtt" rmsdtt} msg] } {
        puts "VMD RMSDTT could not be started:\n$msg"
      }
      

    Remember to use the correct path to your plugins directory; i.e, if you saved the plugin in /home/user/myplugins, write that same path in the auto_path line.

    Note: If you created the .vmdrc file, remember to add menu main on to get the main menu back.

  3. To be able to use MS Excel on Windows you will need the tcom tcl library. Download the library and install it (i.e. in your plugins directories) (in general, only the tcom/lib/tcom directory is needed). Then add the following to your vmd.rc file:
    lappend auto_path {/path/to/tcom/lib/tcom}
    if { [catch "package require tcom" msg] } {
      puts "VMD tcom package could not be loaded:\n$msg"
    }
    
  4. Start VMD. The RMSDTT plugin should be accessible from the Extensions menu.

Author

Luis Gracia, PhD
Weill Medical College, Cornell University
1300 York Ave, New York, 10021 NY

Send comments and feedback to Luis Gracia.

Luis Gracia | last updated: Tuesday, 4 Apr 2006