ViewDock ViewDock icon and HearDock HearDock icon

BACKGROUND

DOCK calculates possible binding modes, given the structures of ligand and receptor molecules. In a DOCK search, each compound (typically from a database of many thousands) is treated as a possible ligand, and the structure of a target macromolecule is treated as the receptor. Various scoring methods are used to identify the most favorable binding modes of a given molecule and then to rank the molecules. The output consists of many candidate ligands in their best-scoring poses. It is then up to human users to look through the results and decide which compounds should be tested in the real world. For further details, please consult the DOCK web site.

ViewDock facilitates the interactive selection of compounds from the output of DOCK. With a single mouseclick, a user can view an individual docked molecule and simultaneously display its name, scores, and other information. The list of molecules can be sorted in various ways, edited to remove compounds deemed uninteresting, and saved for future reference. For a step-by-step example, see the ViewDock tutorial.

HearDock (only available for SGI systems) includes the functionality of ViewDock, but also allows sounds to be associated with properties such as score.

The state of ViewDock is included in saved sessions. DOCK (versions 4, 5, 6) scoring grids can be visualized with Volume Viewer.

STARTUP AND INPUT

There are several ways to start ViewDock (or HearDock), a tool in the Surface/Binding Analysis category. Starting the tool brings up a dialog for opening the file of docked molecules. Format options include:

Another way to start ViewDock is with the viewdock command, which loads a specified file of docking results.

THE LISTBOX

When the file of docked molecules has been read, the ViewDock ListBox will appear and the first molecule will be displayed in the main Chimera window. In most cases, a user will also display the receptor structure, which can be opened either before or after the file of docked molecules. ViewDock is controlled via the ListBox and its menu.

ViewDock ListBox

The upper panel lists the compounds along with their status settings (S) and the values of various descriptors. Available descriptors can be shown or hidden using the Column menu. Descriptors are simply read from the file of docked ligands; what descriptors are present and available depends on the DOCK version and options used. In addition, the HBonds... Add Count... options can be used to identify hydrogen bonds and add their counts as new descriptors. The listing can be sorted by the values in any column by clicking the column header. Clicking the header once sorts the entries in order of increasing value and places an up arrow in the header. Clicking again sorts the entries in decreasing order and places a down arrow in the header.

A compound can be chosen by clicking on its line in the upper panel. Multiple compounds may be chosen at once. Chosen lines are highlighted and only the chosen compounds are displayed. Ctrl-click adds to an existing choice rather than replacing it. To highlight a block of compounds without having to hold down the mouse button, click on the first (or last) and then Shift-click on the last (or first) in the desired block. For results from DOCK versions 4 and later, the lower panel shows detailed information when a single compound is chosen.

Any numerical descriptors from the input file are assigned as attributes of models (each docked molecule is treated as a separate model), with the attribute name derived from the descriptor name by substituting underscores for spaces, changing any capital letters to lowercase, and appending dock_ at the beginning. For example, the descriptor Energy score is assigned as an attribute of models named dock_energy_score. Such assignments allow docked compounds to be rendered (with colors or radii) or selected by their descriptor values.

Hide dismisses the ListBox without exiting from ViewDock. Quit exits from ViewDock, leaving the structures open (like File... Exit in the ListBox menu). Help brings up this manual page in a browser window. When the ListBox is hidden or obscured by other windows, it can be resurrected using the Raise option for the ViewDock instance in the Tools menu. The instance's Hide and Quit options have the same effect as the corresponding ListBox buttons.

COMPOUND STATUS

Three mutually exclusive states are possible for each compound: Viable, Deleted, or Purged. The initial state is viable, but if a molecule is deemed uninteresting, it can be changed to deleted or purged by clicking on the checkboxes near the bottom of the ListBox. Clicking a status checkbox changes the status of all chosen compounds.

In general, Viable compounds are interesting (or have not been looked at yet), Deleted compounds are less interesting but may deserve another look, and Purged compounds are meant to be discarded. The Compounds menu controls (based on status) which compounds are listed in the upper panel.

When structures are written out using File... Rewrite, viable and deleted but not purged structures are included in the output file. In contrast, structures in all three states are included in a file written out using File... Save or File... Save As. The states of the molecules are also recorded in these output files, which can be input to later sessions of ViewDock (or HearDock). No matter how the molecules are sorted in the ListBox, they remain in their original order in a saved PDB file.

There are other ways to change status besides clicking a status checkbox:

HYDROGEN BONDS

The HBonds menu allows incorporation of hydrogen bonding information calculated with FindHBond. If FindHBond has not already been run, its interface will appear. If FindHBond has been run, but one wishes to run it again with different parameters, it can be started from the main Chimera Tools menu or by using another of the standard ways to start tools. When FindHBond is used to help screen compounds in ViewDock, the inter-model setting is most efficient.

The HBonds... Add Count... options add hydrogen-bonding statistics as descriptor columns:

The hydrogen bond descriptors can then be used to choose compounds by value or guide manual changes in compound status.

HBonds... Change Compound State allows changes in compound status based on the number of hydrogen bonds to the receptor or to a selected portion of the the receptor. To use different hydrogen-bonding criteria, remember to change the compounds back to their previous states before re-running the hydrogen bond calculation. (If many status changes have been performed manually, it may be helpful to use File... Save or File... Save As before making further changes.)

CHOOSING BY VALUE

Compounds can be chosen based on the values of any (or all) of the available descriptors, including hydrogen-bond statistics calculated with FindHBond (see the HBonds menu).

Compounds... Choose by Value brings up a dialog for defining value ranges. Whichever descriptors are shown in the ListBox are available for choosing compounds. To add a descriptor to the Choose by Value dialog, simply show the descriptor in the ListBox using the Column menu. Any descriptor shown in the ListBox but not of interest for choosing compounds can be disregarded by unchecking the corresponding checkbox.

Checkboxes under Choose from designate which compounds should be considered, on the basis of their status. The values displayed in the histograms and/or lists only come from this designated set of compounds.

Numerical descriptors are depicted in a histogram. Within a histogram, two markers are shown as vertical bars. Clicking on a marker shows its Value in black. Clicking elsewhere within the histogram shows the Value (X-coordinate) of the mouseclick in gray. A marker can be moved by changing its Value and then pressing Enter (return) or by dragging it horizontally with the left mouse button. Holding the Shift key down reduces the speed (mouse sensitivity) of marker dragging tenfold, allowing finer control. Only compounds with values between the markers (inclusive) will be chosen.

String-valued descriptors are displayed as a list. Clicking a line highlights just that line; Ctrl-click adds to an existing choice rather than replacing it. To highlight a block of lines without having to hold down the mouse button, click on the first (or last) and then Shift-click on the last (or first) in the desired block. Only compounds with highlighted values will be chosen.

A message of the form N of M compounds applicable reports how many of the compounds (N) of those designated to choose from (M) meet all current criteria. These compounds will be chosen when Apply or OK (which also dismisses the dialog) is pressed.

COMPOUND DISPLAY

Display can controlled in many ways. Because of the potential for confusion, however, it is recommended that the ViewDock ListBox, not commands, be used to control the display of the docked compounds.

The remainder of this section explains the interplay among display controls, but can be skipped if only the ViewDock ListBox will be used to control compound display.

Individual compounds are submodels of a single model (model 0, for example) and are specified #0.1, #0.2, etc. (see atom specification syntax). Even though they are termed submodels, the display of each compound can be independently disabled or enabled at the model level (see display hierarchy). Within ViewDock, there are two ways to control model-level display:

Compounds whose display has been disabled at the model level are unresponsive to the command display. Ways of enabling display at the model level (and restoring responsiveness to the command display) include choosing the compound in the ListBox upper panel, using one of the Chimera... Show options, and using the command modeldisplay.

For example, if only the first five compounds are display-enabled,

Command: ~disp #0.2
and
Command: disp #0.2
but not
Command: disp #0.8
will work. In contrast, other commands like
Command: label #0.7
and
Command: color pink #0.6
are executed, although their effects may not be visible until the corresponding molecules are displayed. The command
Command: disp #0
will display only the compounds whose display is enabled at the model level.


LISTBOX MENU

File

Compounds

The Compounds menu controls which compounds are listed in the top panel of the ListBox. Any combination of the first three entries may be checked to indicate which status types should be listed.

Column

Selection

The following checkboxes are mutually exclusive and control how picking from the graphics window is handled by ViewDock.

Chimera

The Hide and Show options toggle the display status of compounds in the main Chimera window.

HBonds

The HBonds menu allows incorporation of H-bond information calculated with FindHBond.

Movie

Sonify (in HearDock only)

HearDock is only available for SGI systems. Any listed descriptor can be selected for expression as sound. Upon selection of a descriptor, the Sonification Parameters panel (shown below) appears and allows specification of instrument, minimum and maximum pitch, and optionally, what range of the descriptor value to sonify. Lower pitch is associated with lower values of the descriptor (smaller positive or larger negative values).

Sonify can be invoked several times to set sonification parameters for multiple descriptors, or to change parameters for a single descriptor. For example, a trombone could be used for the electrostatic interaction energy and a tuba for the van der Waals interaction energy. To subsequently hear only the electrostatic term, set the instrument representing the van der Waals term to None. To hear the descriptors, it is necessary to use Movie... Play. As each molecule is displayed, the tone(s) associated with it is (are) played. The length of time per molecule can be adjusted, or the movie stopped, as described above.


EXAMPLE FILES

Note that DOCK 4 interprets information about each molecule as Name or Description based on where this information occurs in the input database file. Here are examples of a single molecule from (A) a Mol2 database file input to DOCK 4, (B) a Mol2 output file from DOCK 4, and (C) a PDB output file from DOCK 4. Multiple-molecule files are simply concatenations of the data for single molecules. Number is sequential number in order of occurrence in the file from DOCK, and Source num is the sequential number in order of occurrence in the database file input to DOCK.


UCSF Computer Graphics Laboratory / March 2007