We will use UCSF Chimera (Pettersen et al. 2004) for molecular graphics visualization, but the software is capable of many complex operations (see also Acknowledgments below.)
The material presented here is inspired by the RCSB Protein Data Bank YouTube1 tutorial series provided as a list2 created by Dr Shuchismita Dutta, Ph.D.
Chimera is free for academic use and is available for download at https://www.cgl.ucsf.edu/chimera/download.html but is already installed on the iMacs under the Macintosh OS.
If you are using the iMac for the first time you may need to review and update the mouse button assignment by going to the “Apple” menu at the top left:
Apple
> System Preferences...
> Mouse
“Standard” setting should be set as (if you are right-handed) :
Note: If you are left-handed inverting Left and Right assignment might be best for you.
This button will invite you to act on suggested actions.
Example file: 1mbo
. Structure and refinement of oxymyoglobin at 1.6 angstroms resolution (Phillips 1980).
On a Macintosh Chimera
is installed in /Applications
.
You can also find Chimera
by using Spotlight Search
: the “magnifying glass” icon on the top right of the Mac.
On a Windows computer you would type Chimera
in the search area next to the Start button in Windows 7, or within a “magnifying glass” text entry, also bottom left, with “Cortna” helper in Windows 10.
Double click on the Chimera icon.
For the very first launch, Chimera will open a black window with a menu bar at the top. (Note that on a Macintosh these menus are also found on the Finder
menu bar at the very top of the computer screen.)
After this first time, new launches will show a blue screen with a history list of previously opened structures shown on the right hand side. You can click on any of these to reopen the structures in a new session.
The File
menu offers 2 ways to load a structure:
File
> Open...
that is already saved on your computerFetch by ID...
to fetch by an identifier. In that case the computer must be connected to the Internet.In both cases you would need to know the PDB ID code of four characters gleaned e.g. from the PDB site or a publication.
When the structure is loaded, it is presented in a default preset mode showing ribbons and select atoms as stick representation.
Fetch PDB id 1mbo
:
Note that the 1mbo
protein structure is shown as a ribbon, while the ligands and some specific amino acids are shown as stick models in this default preset.
Later we will review this default presentation mode with the “Presets
” options.
Note: When using
Fetch
Chimera will save the downloaded file within the user’sDownloads/Chimera
directory. In this case the file will be saved as1MBO.pdb
and can be explored with a simple word processor as it is a “plain text” file. On a Mac one would use the built-in text editorTextEdit
and on WindowsWordPad
or perhapsNotePad
.
Hovering the mouse over a depiction will reveal the residue name and the name of the chain pointed out by the mouse.
Hover the mouse over at N-t to reveal info.
For example, if you place the mouse over the N-terminus and stop moving a temporary label within a box will reveal: VAL 1.A
which means: chain A, residue 1 which is a valine. (See above illustration.)
Later we will learn how to select and label an amino acid.
An external, 3-button mouse is best to interact with the displayed molecule.
x
and y
axes.The same movement when performed far from the molecule, e.g. far on the left side, will rotate the molecule around the z
axis.
Right button: Zoom. Press the right button to zoom in (left-to-right or top-to-bottom) and out (right-to-left or bottom-to-top)
Note: if these button do not appear to work as planned, review the mouse preference settings as described above.
Once you have found a suitable view, you can try the “Presets” menus that alter the display of the molecule with preset options.
The first one ** Presets
> Interactive1(ribbon)
** looks very similar to the current display when a molecule is opened.
The difference is that the ribbon is rainbow-colored from a blue-colored N-termial blue to a red-colored C-terminus.
Engage the Presets
menus to explore the preset representations that are offered.
Presets menus |
check if seen | Notes |
---|---|---|
Interactive 1 (ribbons) | [ ] | ribbon is rainbow colored |
Interactive 2 (all atoms) | [ ] | protein is wireframe, ligands are spheres |
Interactive 3 (hydrophobicity surface) | [ ] | blue: hydrophilic; orange: hydrophobic |
Publication 1 (silhouette, rounded ribbon) | [ ] | |
Publication 2 (silhouette, licorice) | [ ] | |
Publication 3 (depth-cued, rounded ribbon) | [ ] | |
Publication 4 (depth-cued, licorice) | [ ] |
Note: the Interactive menu options create a view. The Publication menus alter the view prepared by the Interactive menu. Background is white.
To save an image, use the menu cascade: File
> Save Image...
File name:
Provide a name for the file; if needed navigate to where the file will be saved.File type:
default format is PNG
, alter format if desired from this pull-down menu.Image width:
and Image height:
image size is same as current Chimera display by default; alter if needed with these menues.Save
Note other options such as Transparent background
button etc.
Save a test image on the Desktop
A session file contains most of the information to recreate the current state of the Chimera display. It retains selections and molecules currently uploaded within Chimera.
To save a session use the menu cascade: File
> Save Session As...
A session file can be reopened with the menu option: File
> Restore Session...
Save the current state of Chimera in a session file on the Desktop
When you are done with a session you can use the menu File
> Close Session
to close the session an keep Chimera running and open and start a new session.
Close this session
If you are finished, you can can use File
> Quit
to exit Chimera gracefully and return to your computer.
Example file seen previously: myoglobin structure 1mbo
.
If you are continuing from the previous section you can simply close the session and re-open the previous session with 1mbo
by clicking on the right hand panel on 1mbo [PDB]
.
The current view of the protein should be similar to the first Presets
menu version with the protein shown as ribbon (perhaps rainbow colored) and the heme and amino acids close to the heme shown as stick models.
In the following exercise we will learn to select the amino acids that are shown as stick by distance criterion and alter their color.
In this exercise we will select amino acid residues found within 5 Å of the heme group.
Select
> Residue
> HEM
to select the heme group.Select
> Zone...
to select atoms that are within 5 Å of the heme groupNote that if the bottom button is selected, if any one atom is selected with the distance criterion then all of the atoms of the amino acid residue that it is part of will also be selected, therefore selecting whole amino acids in the process.
Action
> Color
> by element
to alter the color of the stick representations
In this exercise we will measure the distance between the iron in the heme group and the coordinating histidine. We will select the specific atoms with the mouse.
Measure the distance between FE
in heme group and his 93.A
.
Action
> Focus
, or rotate with left button and zoom with right button drag.)Ctrl
+Shift
and then left-button click to select and add to the selection. The heme iron should appear as an orange sphere and the closest atom of histidine 93 as stick.Struture Measurements
panel: Tools
> Structure Analysis
> Distances
and click the Create
button.In this exercise we will explore how to show hydrogen bonds and using selection methods to show them only for specific residues, in our example on the hydroxyl of Serine 92 (SER 92.A OC
.)
The first step is to find all hydrogen bonds with the defaults options.
Select hydrogen bonds.
Tools
> Structure Analysis
> FindHBond
to find all hydrogen bonds in the structure using default parameters already selected.Apply
button at the bottom to reveal all H-bonds that can be displayed: within the helices, displayed side-chains and water molecules.We will now display only specific H bonds for a residue of interest, and as an example it will be the hydroxyl of Serine 92. In order to do that, we will use a menu within the H-Bond Parameters
panel:
Ctrl
-Click on the hydroxyl group of Serine 92 (it should be displayed as stick)Only find H-bonds
with at least one end selected
located in the middle of the H-Bond Parameters
panel.Apply
button: now only the relevant H-bonds are shown.
Show molecular surface.
Actions
> Surface
> Show
will display a molecular surface.The color may be a default color (e.g. beige) or might be a rainbow color depending on the state of your molecule when you activate the menu.
Note: An alternate way to bring about a surface, colored by hydrophobicity by default is the preset setting obtained with the menu cascade:
Presets
> Interactive 3 (hydrophobicity surface)
Exerpt from the Chimera web page on Coulombic Surface Coloring3
Coulombic Surface Coloring calculates electrostatic potential according to Coulomb’s law: \[\Phi = \Sigma[q_i/\epsilon d_i]\] \(\Phi\) is the potential (which varies in space), \(q\) are the atomic partial charges, \(d\) are the distances from the atoms, and \(\epsilon\) is the dielectric, representing screening by the medium or solvent. A distance-dependent dielectric (\(\epsilon = Cd\) where \(C\) is some constant) is sometimes used to approximate screening by implicit solvent.
Note: Poisson-Boltzmann calculations are more complex and, if done correctly, more accurate than simple Coulomb’s law approaches. However, a Coulombic potential may suffice for visualization. Software that calculate a Poisson-Boltzmann grid are not included in Chimera.However, Chimera does include interfaces to such programs: DelPhiController requires a local (user-installed) copy of DelPhi, and the APBS tool can use either a web service or a locally installed copy of APBS (Adaptive Poisson-Boltzmann Solver).
Coulombic Surface Coloring.
Tools
> Surface/Binding Analysis
> Coulombic Surface Coloring
will open a window meant to color surfaces. The she top part of the window will list the surface(s) currently available within Chimera.Apply
button o create an approximate electrostatic Colombic surface coloring. Keep the default values as well as the default suggested 3-colors of red, white and blue.Coulombic Surface Coloring | Notes |
---|---|
Note 1: to draw a scale click the button at the very bottom Create corresponding color key and use the mouse to draw at the desired size. Use the button Reverse ordering of above if necessary. Note 2: You can optionally use the Viewing panel (side view) to alter the image as we have seen in a previous section to add a silhouette contour or change the lighting option e.g. to “ambient.” |
.
Colored surface “2-lights” | Colored surface “ambient” |
---|---|
Example files: myoglobin 1mbo
and hemoglogobin 2hhb
both used in previous sections.
Open hemoglobin structure 2hhb
.
2hhb
structure.Tools
> Depiction
> Rainbow
chain
in the Rainbow
control panelApply
buttonClose
buttonChains will be colored one solid color per chain as shown within the Color range
options on the left side of the Rainbow
control panel (see image.)
Fetch myoglobin structure 1mbo
and superimpose to a hemoglobin chain.
File
> Fetch
menu to add the structure 1mbo to the display. It will appear partly overlapping the hemoglobin structure.We will now superimpose the myoblogin structure to one of the hemoglobin chain. For now we use default parameters
Tools
> Structure Comparison
> MatchMaker
to open the MatchMaker
control panel.2hhb (#0)
and1mbo (#1)
OK
The 1mbo
myoglobin structure will be be moved to overlap chain D
of the hemoglobin structure.
Note: it is possible to select the exact chain to match by changing the option to:
Specific chain(s) in reference structure with specific chain(s) in match structure
which is the 3rd button within theChain pairing
box within the panel.
Hide unmatched chains.
Here we will select unmatched chains and hide them:
Select
> Selection mode
> append
to select multiple chains easilyActions
> Ribbon
> hide
Actions
> Atoms/Bonds
> hide
Select
> Selection mode
> replace
to return selection to replace
Note that for chain A, there is now a choice between the 2hhb
an 1mbo
structures.
Now we can focus on the superimposed chains:
Select
> Chain
> A
> 1mbo #1
to select myoglobin chainAction
> Focus
to zoom on the superimposed chainsCtrl
-Click on the black background to unselect.You can now rotate the structure and note where there are similarities and differences.
For example you may look for the conserved histidine coordinating the heme group, His 92 in 2hhb
and His 93 in 1mbo
.
It may be useful to view each chain separately rather than as a superimposed composite. For this we will use the Model Panel
to turn on or off the display of molecules. The Model Panel
can be activated from the menu cascade Tools
> General Controls
> Model Panel
or simply from Favorites
> Model Panel
Display or hide molecules.
Model Panel
(see menu cascades above)Shown
columnIn this exercise we will compare the sequences of the superimposed structures.
Compare aligned sequences
Tools
> Structure Comparison
> Match -> Align
will open the Create Alignment from Superimposition
control panel.A
to all structures to the alignment that was performed above:2hhb (#0) chain D
1mbo (#1) chain A
OK
The alignment will be shown:
Can you spot the conserved H
histidines? [_] Yes [_] No
Can you spot other conserved, identical residues ?
You can write their name/number ________________________
Compute % identity
We can compute a “percent identity” with the menu contained within the matched sequence panels:
Info
> Percent identity...
will open the Compute Percent Identities
control panelall sequences
to the desired chain. For example:2hhb, chain D
1mbo
, chain A`OK
- the % identity will appear at the bootom of the sequence alignment window as 24.66%.What value did you get? _________________________
We have explored the Favorites
options execept Command Line
The command line options will be explored in a separate workshop.
However, you can learn here a very simple and useful command to alter the background color. The default color is black, in this exercise we will change it to white.
Change background color by line command
Favorites
> Command Line
will open the command line at the bottom of the 3D view panel.set bg_color white
and the background will change accordingly.set bg_color black
There are many settings that can be changed on the fly, but also some that can be set as preferences.
The Favorites
> Preferences
control panel allows to change many settings that may become permanent in the Chimera installation.
This is mentioned here as there is one change that we alluded before for the label font size. This is the place where many aspects of labels can be changed.
Change label size
Favorites
> Preferences
Category
pull-down menu to Labels
16
to e.g. 36
From the UCSF Chimera web site:
UCSF Chimera is a highly extensible program for interactive visualization and analysis of molecular structures and related data, including density maps, supramolecular assemblies, sequence alignments, docking results, trajectories, and conformational ensembles. High-quality images and animations can be generated. Chimera includes complete documentation and several tutorials, and can be downloaded free of charge for academic, government, nonprofit, and personal use. Chimera is developed by the Resource for Biocomputing, Visualization, and Informatics (RBVI), funded by the National Institutes of Health (NIGMS P41-GM103311).
UCSF ChimeraX (or simply ChimeraX) is the next-generation molecular visualization program from the RBVI, following UCSF Chimera.
Fermi, G., M. F. Perutz, B. Shaanan, and R. Fourme. 1984. “The crystal structure of human deoxyhaemoglobin at 1.74 A resolution.” J. Mol. Biol. 175 (2): 159–74.
Pettersen, E. F., T. D. Goddard, C. C. Huang, G. S. Couch, D. M. Greenblatt, E. C. Meng, and T. E. Ferrin. 2004. “UCSF Chimera–a visualization system for exploratory research and analysis.” J Comput Chem 25 (13): 1605–12. http://www.cgl.ucsf.edu/chimera.
Phillips, S. E. 1980. “Structure and refinement of oxymyoglobin at 1.6 A resolution.” J. Mol. Biol. 142 (4): 531–54.