Multalign Viewer Multalign Viewer icon

Multalign Viewer allows viewing and manipulation of sequence alignments. If any corresponding structures are open in Chimera, highlighting regions in the sequence alignment will highlight portions of the structures and vice versa. In addition, structures can be superimposed based on the sequence alignment. The state of Multalign Viewer (along with sequence alignment data) is included in saved sessions. For an informal introduction, see the Sequences and Structures tutorial.

STARTUP AND INPUT

There are several ways to start Multalign Viewer, a tool in the Structure Comparison category. Explicitly starting Multalign Viewer brings up a dialog for opening files. File formats that can be read and written are:

Format Prefix Suffix
Aligned FASTA afasta: .afasta
.afa
.fasta
.fa
Aligned NBRF/PIR pir: .pir
Clustal ALN aln: .aln
.clustal
.clustalw
.clustalx
GCG MSF msf: .msf
GCG RSF rsf: .rsf
Stockholm pfam:
hmmer:
.pfam
.sto

One of these types can be chosen from the File type: pulldown menu near the bottom of the dialog, or the type can be inferred from the suffix(es) of the chosen file(s). Once a file has been read, a sequence window containing the alignment appears; if multiple files are read, each is shown in a separate sequence window.

Multalign Viewer menu items are explained within the following sections, and a full listing with short descriptions is included at the end of this page.

THE SEQUENCE WINDOW

Multalign Viewer sequence window

The figure shows part of the sequence window contents for the input file apoex.fa. Font size and sequence wrapping behavior can be controlled with the Layout preferences. Alignment regions (colored boxes enclosing one or more residues) are described below. The coloring of residue one-letter codes can be controlled with the Analysis preferences. One can search for occurrences of a particular string or pattern of residues in one or all of the sequences using Tools... Find Subsequence or Tools... Find PROSITE Pattern in the Multalign Viewer menu.

When the cursor is in the sequence window, the Page Down key (or space) moves the view down to start with the block below the topmost block whose beginning is currently visible; Page Up (or Shift-space) moves the view up to start with the block above the topmost block whose beginning is currently visible.

The Hide button closes the sequence window without changing the state of Multalign Viewer. The sequence window can be reinvoked using the Raise option for the Multalign Viewer instance in the Tools menu (abbreviated MAV). This is also useful when the window has become obscured by other windows. Help opens this manual page in a browser window, and Quit exits from Multalign Viewer.

ALIGNMENT HEADERS

By default, header lines named Consensus and Conservation are shown above the alignment in the sequence window. The consensus sequence shows the residue type most prevalent (not necessarily >50%) at each position in the alignment; highly conserved residues (80% or greater) are capitalized and shown in purple, except that the completely conserved residues are shown in red. The Analysis preferences control how consensus and conservation are determined and represented.

Arbitrary, user-defined header lines may also be shown. New header line contents can be loaded from an external header definition file (Headers... Load in the Multalign Viewer menu), while existing header line contents can be saved to the same type of file (Headers... Save in the Multalign Viewer menu). Loading values for a header with the same name as an existing header will create an additional header line, not replace the values in the existing header line.

GC (Generic per-Column) annotations in alignments in Stockholm format are also treated as headers.

The display of headers in the sequence window can be turned on and off with their respective checkboxes in the Multalign Viewer Headers menu. Clicking in a header line selects any structure residues associated with that column of the alignment.

Values of headers are available as residue attributes named mavHeader_name in any associated structures. Numerical attributes will appear in the attribute lists of Render/Select by Attribute; character attributes will be listed only in the Select by Attribute portion. It may be necessary to Refresh the attribute menus or histogram/list in these tools to update them with any changes in Multalign Viewer header information.

ALIGNMENT REGIONS

Alignment regions can be defined manually or by any of several operations within Multalign Viewer. A region is shown using a colored rectangle or outline in the sequence window; a single region can contain any number of disjoint and/or abutting rectangular blocks.

A region can be created manually by dragging with the left mouse button within the sequence window. Dragging downward into the following block highlights to the end of the preceding block. The initial region colors can be controlled in the Regions preferences. Shift-dragging with the left mouse button adds to the active region. Ctrl-dragging creates a new region and makes it the active region.

The active region is the region most recently created manually, clicked on in the sequence window, or designated as Active in the Region Browser. The corresponding parts of any associated structures are highlighted in the main Chimera window by becoming selected. Only one region can be active at a time. In the sequence window, the active region is indicated with a dashed outline; clicking the active region deactivates it, and clicking a different region deactivates the former active region and makes the new region active. A region with no interior color is only responsive to clicks on its borders. Where regions overlap, only the highest is responsive to clicks.

The Region Browser (Tools... Region Browser in the Multalign Viewer menu) can be used to manage the existing regions. Checkboxes reflect and control which region is Active (if any) and which are Shown in the sequence window. Clicking an entry in the Region Browser highlights the line and chooses the region. More than one region can be chosen at a time. The chosen region(s) are not necessarily active.

The Border color and Interior color of the chosen region(s) can be changed using the adjacent color wells.

The Delete button on the Region Browser deletes the chosen region(s). The region named Chimera selection cannot be deleted, but it can include zero residues. Since regions may overlap, a region can be considered higher or lower than another. Raise puts the chosen region(s) in front of any other regions in the sequence window and moves the corresponding entries in the Region Browser to the top. Lower puts the chosen region(s) behind any overlapping regions in the sequence window and moves the corresponding entries in the Region Browser to to the bottom. Copy copies the chosen region(s) and requests a new name for each copy, while Rename requests a new name for each chosen region.

Close dismisses the Region Browser and Help opens this manual page in a browser window.

Alternatively, when the cursor focus is in the sequence window, the up arrow and down arrow keys can be used to raise and lower the active region, respectively; pressing the Delete key will delete the active region.

Besides manual creation, there are several ways to generate or change the contents of regions:

EDITING AND SAVING THE ALIGNMENT

The active region can be moved using Ctrl-left arrow and Ctrl-right arrow. Residues cannot be added or removed, but gaps can be created, extended, reduced, or removed. To accommodate the shift, sequence flanking the active region will be pushed after intervening gaps are closed. The alignment generally cannot be restored by reversing the shift because the positions of gaps are forgotten as the gaps are closed.

Choosing File... Save As from the Multalign Viewer menu brings up a dialog for saving the sequence alignment to a file. It is possible to save just the active region to a file; the region may consist of disjoint sections, but each sequence (row) included must contain the same set of columns. The option to omit all-gap columns from the output is especially useful for saving regions with fewer sequences than the starting alignment. The formats available for saving are the same as those that can be read. An alignment saved in Stockholm format will automatically include annotations describing the secondary structure of any associated structures.

The sequence window contents can also be saved as an EPS file (File... Save EPS in the Multalign Viewer menu).

PERCENT IDENTITY

The percent identity between any pair of sequences in the alignment can be calculated. Tools... Percent Identity in the Multalign Viewer menu opens a dialog for specifying the sequences (by name from the pulldown menus labeled Compare and with) and what to divide by:

The order in which the sequences are specified is not important.

Apply calculates the percent identity without dismissing the dialog, while OK calculates the percent identity and dismisses the dialog. The percent identity is reported in the status area near the bottom of the sequence window and the Chimera Reply Log. Close simply dismisses the dialog, and Help opens this manual page in a browser window.

ANNOTATIONS

Sequence alignment files may contain various annotations.

Alignment annotations (those describing the alignment as a whole) can be viewed and altered by choosing Tools... Alignment Annotations from the Multalign Viewer menu. The dialog shows any annotations from the input file, each of which can be changed by typing in the adjacent field. Changes are not assessed for correctness.

Each annotation in the top section of the panel has two parts, a name and a value (generally text). For example, MSF length is a name and 366 is a possible associated value. A new annotation can be added by clicking New and specifying a name. The name will appear in the top section of the panel, and a value can be entered by typing in the adjacent field. Values can have multiple lines even though only one is shown; hitting return when typing into a value field retains the previous information and starts a new line. Clicking Delete and then choosing an annotation name removes the corresponding annotation. Multiple lines of free-form text can be placed in the Comments area in the lower section of the panel.

Not all formats can accommodate annotations. Thus, saved files may not include this information or reflect any changes that have been made. Only the Stockholm format accommodates arbitrary annotations, whereas both Stockholm and RSF formats can include comments. When an alignment in Stockholm format is read into Multalign Viewer, any GF (Generic per-File) markups are interpreted as alignment annotations and GC (Generic per-Column) markups are treated as headers: they can be shown above the alignment and are available as residue attributes in any associated structures.

When an alignment is saved in Stockholm format, the saved file will automatically include GR (Generic per-sequence and per-column) markup lines describing the secondary structure of any associated structures. The GR markup line for a sequence and its associated structure starts with

#=GR seq-name Chimera_actual_SS_struct-name
The structure name is included because there can be more than one structure associated with a given sequence. Chimera-generated GR secondary structure markups may include the symbols:

symbol meaning
. gap in sequence
H helix
E strand
C other structure
X sequence residue not associated with structure residue

If a sequence did not have a GR SS markup in the original input file and is associated with only one structure, another markup line with the same contents but named SS instead of Chimera_actual_SS_struct-name will also be included. GR markups are not displayed by Multalign Viewer.

SEQUENCE-STRUCTURE ASSOCIATION

A sequence in the alignment will be associated automatically with a structure chain if their sequences can be aligned without too many mismatches. The allowable number of mismatches is user-specified as a proportion of the total number of residues in the structure chain (1/10 by default; see the Structure preferences). For automatic association, gaps in the structure sequence relative to the sequence in the alignment file can only occur where residues are missing from the structure (for example, a flexible loop with insufficient density for coordinates to be determined). The order in which the sequence and structure files are opened does not matter. Associations are reported in the status area near the bottom of the sequence window and the Chimera Reply Log. A structure (even if it has multiple chains) cannot be associated with more than one sequence, but a single sequence can be associated with more than one structure. If more than one sequence matches a given structure chain, the single best-matching sequence is associated.

The names of structure-associated sequences are shown in bold over a box indicating the model-level color of the structure. For example, APE_HUMAN in the sequence window figure is associated with a structure whose model-level color is white. If multiple models are associated with the same sequence, the box is shown as a dark green dashed outline. When the cursor is placed over the name of a sequence, information on the sequence and any associated structure(s) is shown in the status area near the bottom of the sequence window. When the cursor is placed over a sequence residue that is associated with a residue in one or more structures, information on the structure residue(s) is given. Association information can be saved to a file.

Sequence-structure association may create new regions. If all residues of a sequence are matched, no region is created. Otherwise, by default

The region names also specify the structure and chain that were sequence-associated and the number of gaps or mismatches (indicated by "..." above).

Association enables several types of sequence-structure crosstalk:

Adding or Changing Associations

Associations can be added or changed using the Structure/Sequence Associations panel (Structure... Associations in the Multalign Viewer menu). The name of each open structure model is shown next to pulldown menus for specifying chain (if applicable) and sequence name. An association is indicated by choosing the name of the sequence to be associated with the model and (if applicable) with which chain the association should be made. The choices for sequence to associate include those in the alignment displayed by Multalign Viewer and none. The specified associations will be made no matter how inappropriate they are.

Apply performs the indicated association(s) without dismissing the panel. OK performs the indicated association(s) and dismisses the panel, Close dismisses the panel without changing the associations, and Help opens this manual page in a browser window.

Automatic Structure Loading

Structure... Load Structures in the Multalign Viewer menu can be used to open structure files corresponding to sequences that are not already structure-associated. The corresponding structure files are determined from the sequence names using the rules given in the Structure preferences, then retrieved and opened (as described for Fetch by ID). If Automatically load Structures is turned on (also in the Structure preferences), this will occur as soon as an alignment is read.

Regions Defined by Secondary Structure

Regions can be defined within structure-associated sequences based on helices and sheets in the corresponding structures. Structure... Secondary Structure... show actual in the Multalign Viewer menu creates regions named structure helices and structure strands colored goldenrod and lime green (see the named colors), respectively. Helices and strands are defined by HELIX and SHEET records in the input file, or if these are not present, using ksdssp. (Regardless of whether the regions exist, however, an alignment saved in Stockholm format will automatically include annotations describing the secondary structure of any associated structures.)

In sequences not associated with any structures, Structure... Secondary Structure... show predicted in the Multalign Viewer menu creates regions named predicted helices and predicted strands colored gold and light green (see the named colors) that correspond to helices and strands, respectively, predicted using GOR.

RESIDUE ATTRIBUTES

Structure residues associated with residues in a sequence alignment shown by Multalign Viewer are assigned attributes:

Superposition assessment generates an additional residue attribute. Structure residues not associated with residues in the sequence alignment are not assigned attribute values.

Structure... Select (or Render) by Conservation in the Multalign Viewer menu opens the corresponding portion of Render/Select by Attribute, set to the residue attribute mavPercentConserved.

Numerical attributes will appear in the attribute lists of Render/Select by Attribute; character attributes will be listed only Select by Attribute portion. It may be necessary to Refresh the attribute menus or histogram/list in these tools to update them with any changes in Multalign Viewer header information.

STRUCTURAL SUPERPOSITION

Structure... Match in the Multalign Viewer menu allows structures to be superimposed based on the alignment of their associated sequences. The dialog contains two subpanels, each listing the structures associated with sequences in the alignment. One reference structure should be chosen in the left side, but any number of structures to be matched (superimposed) with it can be chosen in the right side. Residues in each Structure to match are paired with the aligned (in the sequence alignment) residues of the Reference structure. If both structures are associated with the same sequence, the correspondence is even more obvious. One point per residue is used for the least-squares fitting: CA in amino acid residues and C4' in nucleic acid residues. The number of atom pairs fitted and the resulting RMSD are reported in the Chimera Reply Log and the status area near the bottom of the sequence window.

Match highly conserved residues only causes only the well-conserved (at least 80%) positions in the alignment to be used for the least-squares fit. These are the positions shown as capital letters in the consensus sequence.

Match active region only causes only the positions in the current active region of the alignment to be used for matching.

Use pseudobonds to show matched atoms indicates that lines (pseudobonds) should be drawn between the matched atoms. For each matched pair of structures, a pseudobond group is created and colored uniquely (in order, the named colors dark green, dodger blue, sienna, yellow, spring green, purple, gray, and coral are used). Each group is named matches of..., where the rest of the name indicates the structures and chains that were matched. The PseudoBond Panel can be used to change the appearance of or delete the pseudobonds.

Create region showing matched residues indicates that the residues paired to generate the fit should be shown as a region named matched residues.

Iterate by pruning long atom pairs until no pair exceeds [x] angstroms refers to an iterative fitting procedure: in each cycle, atom pairs are removed from the match list and the remaining pairs are fitted, until no matched pair is more than x angstroms apart (x=2.0 by default). The atom pairs removed are either the 10% farthest apart of all pairs or the 50% farthest apart of all pairs exceeding the cutoff, whichever is the lesser number of pairs. The result is that the best-matching "core" regions are maximally superimposed; conformationally dissimilar regions such as flexible loops are not included in the final fit, even though they may be aligned in the sequence alignment.

Apply performs the matching (superposition) without dismissing the dialog. OK performs the matching and dismisses the dialog, Cancel dismisses the dialog without performing a match, and Help opens this manual page in a browser window.

SUPERPOSITION ASSESSMENT

Structure... Assess Match in the Multalign Viewer menu opens a dialog for evaluating how well the sequence-aligned residues of structures are superimposed in space. A Reference structure should be chosen from the pulldown menu of structures associated with sequences in the alignment. One or more Structures to evaluate should be chosen from the list of remaining structures. The structures should already be superimposed, but the superposition can have been generated in any way, including manually (not necessarily using Structure... Match).

One point per residue is used for the comparison: CA in amino acid residues and C4' in nucleic acid residues. The distance between each reference-evaluation pair of residues aligned in the sequence alignment (or corresponding to the same residue, if the reference and evaluation structures are associated with the same sequence) is measured. The distances are assigned as a residue attribute of the evaluation structure(s), named matchDist by default. After attribute values have been assigned, Select by Attribute will appear, set to the new residue attribute.

OK performs the comparison and dismisses the dialog, while Close just dismisses the dialog. Help opens this manual page in a browser window.

PREFERENCES

The Preferences section of the Multalign Viewer menu controls preferences specific to Multalign Viewer; settings are saved with other preferences in the Chimera preferences file. Multalign Viewer preferences are grouped together in sections shown as index cards:

Close dismisses the Multalign Viewer preferences tool, and Help opens this manual page in a browser window.

Note that a typed-in value will not be applied until Apply, OK (which also dismisses the preferences tool), or the Enter (return) key has been pressed. Other types of settings take effect right away.

The Layout section of the Multalign Viewer preferences controls the font and arrangement of text in the sequence window.

The Structure section of the Multalign Viewer preferences allows adjustment of criteria for automatic sequence-structure association and automatic loading of structures corresponding to the sequences in an alignment.

The Analysis section of the Multalign Viewer preferences controls how the Consensus and Conservation headers are determined and represented, and how residue one-letter codes are colored.

The Regions section of the Multalign Viewer preferences controls whether a selection will be shown as a region and sets initial region colors.

MENU LISTING

The Multalign Viewer menu includes the following sections:

File

Structure

Headers

The available headers are listed in the bottom part of the menu; checkboxes control which are shown in the sequence window:

Tools

Preferences


UCSF Computer Graphics Laboratory / January 2006