UCSF-Chimera
Short Playlist:
So you need to visualize protein structures, analyze binding sites, or understand why a mutation causes disease? Welcome to Chimera — your molecular visualization workhorse.
What is Chimera?
UCSF Chimera — a free molecular visualization program from UC San Francisco. It lets you:
- Visualize 3D protein/DNA/RNA structures
- Analyze protein-ligand interactions
- Measure distances and angles
- Compare structures (superposition)
- Color by various properties (charge, hydrophobicity, conservation, flexibility)
- Generate publication-quality images
Getting Started
Opening a Structure
From PDB (online):
File → Fetch by ID → Enter PDB code (e.g., 1a6m) → Fetch
From file:
File → Open → Select your .pdb file
Representation Styles
The Main Styles
| Style | What it shows | Use for |
|---|---|---|
| Ribbon/Cartoon | Secondary structure (helices, sheets) | Overall fold |
| Sticks | All bonds as sticks | Detailed view of residues |
| Ball and Stick | Atoms as balls, bonds as sticks | Ligands, active sites |
| Sphere/Spacefill | Atoms as van der Waals spheres | Space-filling, surfaces |
| Wire | Thin lines for bonds | Large structures |
How to Change Representation
Actions → Atoms/Bonds → [stick/ball & stick/sphere/wire]
Actions → Ribbon → [show/hide]
Ribbon for protein backbone + Sticks for ligand/active site residues = best of both worlds
Selection: The Most Important Skill
Everything in Chimera starts with selection. Select what you want, then do something to it.
Selection Methods
| Method | How | Example |
|---|---|---|
| Click | Ctrl + Click on atom | Select one atom |
| Menu | Select → ... | Various options |
| Chain | Select → Chain → A | Select chain A |
| Residue type | Select → Residue → HIS | All histidines |
| Command line | select :153 | Residue 153 |
Useful Selection Menu Options
Select → Chain → [A, B, C...] # Select by chain
Select → Residue → [ALA, HIS, IHP...] # Select by residue type
Select → Structure → Protein # All protein
Select → Structure → Ligand # All ligands
Select → Chemistry → Side chain # Just sidechains
Select → Clear Selection # Deselect everything
Select → Invert (all models) # Select everything NOT selected
Zone Selection (Within Distance)
Select everything within X Å of current selection:
Select → Zone...
→ Set distance (e.g., 6 Å)
→ OK
This is super useful for finding binding site residues!
Coloring
Color by Element (Default)
Actions → Color → by element
| Element | Color |
|---|---|
| Carbon | Gray |
| Oxygen | Red |
| Nitrogen | Blue |
| Sulfur | Yellow |
| Hydrogen | White |
| Iron | Orange-brown |
| Phosphorus | Orange |
Color by Hydrophobicity
Tools → Depiction → Render by Attribute
→ Attribute: kdHydrophobicity
→ OK
| Color | Meaning |
|---|---|
| Blue/Cyan | Hydrophilic (polar) |
| White | Intermediate |
| Orange/Red | Hydrophobic (nonpolar) |
Why use this? To see the hydrophobic core of proteins — nonpolar residues hide inside, polar residues face the water.
Color by Electrostatic Potential (Coulombic)
This is the red-white-blue coloring from your exercise!
Step 1: Generate surface first
Actions → Surface → Show
Step 2: Color by charge
Tools → Surface/Binding Analysis → Coulombic Surface Coloring → OK
| Color | Charge | Attracts... |
|---|---|---|
| Blue | Positive (+) | Negative molecules |
| Red | Negative (−) | Positive molecules |
| White | Neutral | Hydrophobic stuff |
The OK button is disabled if no surface exists. Always do Actions → Surface → Show first!
What to look for:
- Binding pockets often have complementary charge to ligand
- DNA-binding proteins have positive (blue) surfaces to attract negative DNA
- Negatively charged ligands (like phosphates) bind in positive (blue) pockets
Color by B-factor (Flexibility)
B-factor = temperature factor = how much an atom "wiggles" in the crystal.
Tools → Depiction → Render by Attribute
→ Attribute: bfactor
→ OK
| Color | B-factor | Meaning |
|---|---|---|
| Blue | Low | Rigid, well-ordered |
| Red | High | Flexible, mobile |
What to expect:
- Protein core: Blue (rigid)
- Loops and termini: Red (floppy)
- Active sites: Often intermediate
Color by Conservation
When you have multiple aligned structures:
Tools → Sequence → Multialign Viewer
→ (structures get aligned)
Structure → Render by Conservation
| Color | Conservation |
|---|---|
| Blue/Purple | Highly conserved |
| Red | Variable |
Conserved residues = functionally important (active sites, structural core)
Molecular Surfaces
Show/Hide Surface
Actions → Surface → Show
Actions → Surface → Hide
Transparency
Actions → Surface → Transparency → [0-100%]
Use ~50-70% transparency to see ligands through the surface.
Cross-Section (Clipping)
To see inside the protein:
Tools → Depiction → Per-Model Clipping
→ Enable clipping
→ Adjust plane position
Or use the Side View panel:
Tools → Viewing Controls → Side View
Measuring Distances
Method 1: Distance Tool
Tools → Structure Analysis → Distances
Then Ctrl+Shift+Click on first atom, Ctrl+Shift+Click on second atom.
Distance appears as a yellow dashed line with measurement.
Method 2: Command Line
distance :169@OG :301@O34
What Distances Mean
| Distance | Interaction Type |
|---|---|
| ~1.0–1.5 Å | Covalent bond |
| ~1.8–2.1 Å | Coordination bond (metal) |
| ~2.5–3.5 Å | Hydrogen bond |
| ~2.8–4.0 Å | Salt bridge |
| > 4 Å | No direct interaction |
Hydrogen Bonds
What is a Hydrogen Bond?
Donor—H · · · · Acceptor
↑
H-bond (~2.5-3.5 Å)
- Donor: Has hydrogen to give (—OH, —NH)
- Acceptor: Has lone pair to receive (O=, N)
Find H-Bonds Automatically
Tools → Structure Analysis → FindHBond
Options:
- ✓ Include intra-molecule (within protein)
- ✓ Include inter-molecule (protein-ligand)
H-bonds appear as blue/green lines.
Common H-Bond Donors in Proteins
| Amino Acid | Donor Atom | Group |
|---|---|---|
| Serine | OG | —OH |
| Threonine | OG1 | —OH |
| Tyrosine | OH | —OH |
| Histidine | NE2, ND1 | Ring —NH |
| Lysine | NZ | —NH₃⁺ |
| Arginine | NH1, NH2, NE | Guanidinium |
| Backbone | N | Amide —NH |
Common H-Bond Acceptors
| Group | Atoms |
|---|---|
| Phosphate | O atoms |
| Carboxylate | OD1, OD2 (Asp), OE1, OE2 (Glu) |
| Carbonyl | O (backbone) |
| Hydroxyl | O (can be both donor AND acceptor) |
Salt Bridges (Ionic Interactions)
A salt bridge = electrostatic attraction between opposite charges.
| Positive (basic) | Negative (acidic) |
|---|---|
| Lysine (NZ) | Aspartate (OD1, OD2) |
| Arginine (NH1, NH2) | Glutamate (OE1, OE2) |
| Histidine (when protonated) | C-terminus |
| N-terminus | Phosphate groups |
Typical distance: ~2.8–4.0 Å between charged atoms
Coordination Bonds (Metals)
Metals like Fe, Zn, Mg are coordinated by specific atoms:
| Metal | Common Ligands | Distance |
|---|---|---|
| Fe (heme) | His NE2, O₂ | ~2.0–2.2 Å |
| Zn | Cys S, His N | ~2.0–2.3 Å |
| Mg | Asp/Glu O, water | ~2.0–2.2 Å |
Example: In myoglobin (1a6m), the proximal histidine coordinates Fe at ~2.1 Å.
Ramachandran Plot
Shows allowed backbone angles (φ/ψ) for amino acids.
Tools → Structure Analysis → Ramachandran Plot
Regions of the Plot
| Region | Location | Structure |
|---|---|---|
| Lower left | φ ≈ -60°, ψ ≈ -45° | α-helix |
| Upper left | φ ≈ -120°, ψ ≈ +130° | β-sheet |
| Upper right | Positive φ | Left-handed helix (rare) |
Why Glycine is Special
Glycine has no sidechain → no steric clashes → can be in "forbidden" regions (positive φ).
Select → Residue → GLY
Glycines often appear in the right half of the Ramachandran plot where other residues can't go.
Structural Superposition
Compare two similar structures by overlaying them.
Method 1: MatchMaker (Sequence-based)
Tools → Structure Comparison → MatchMaker
→ Reference: structure 1
→ Match: structure 2
→ OK
Output tells you:
- RMSD (Root Mean Square Deviation): How well they align
- < 1 Å = very similar
- 1–2 Å = similar fold
-
3 Å = significant differences
- Sequence identity %: How similar the sequences are
Method 2: Match (Command)
match #1 #0
Restricting Alignment to a Region
To align just the active site (e.g., within 4 Å of ligand):
sel #1:hem #0:hem zr < 4
match sel
Working with Chains
Delete Unwanted Chains
Select → Chain → B
Actions → Atoms/Bonds → Delete
Select Specific Chain
Select → Chain → A
Or command:
select #0:.A
AlphaFold Structures and pLDDT
What is pLDDT?
AlphaFold stores its confidence score (pLDDT) in the B-factor column.
| pLDDT | Confidence | Typical regions |
|---|---|---|
| > 90 | Very high | Structured core |
| 70–90 | Confident | Most of protein |
| 50–70 | Low | Loops, uncertain |
| < 50 | Very low | Disordered regions |
Color by pLDDT
Since pLDDT is in B-factor column, use:
Tools → Depiction → Render by Attribute → bfactor
Or select low-confidence regions:
select @@bfactor<70
Low pLDDT regions in AlphaFold often correspond to regions that are ALSO missing in experimental structures — they're genuinely disordered/flexible, not just bad predictions.
The Hydrophobic Core
Soluble proteins organize with:
- Hydrophobic residues (Leu, Ile, Val, Phe, Met) → inside (core)
- Polar/charged residues (Lys, Glu, Ser, Asp) → outside (surface)
Visualizing the Core
- Color by hydrophobicity
- Use cross-section/clipping to see inside
- Orange/tan inside, blue/cyan outside = correct fold
Protein-Ligand Interaction Analysis
General Workflow
-
Isolate the binding site:
Select → Residue → [ligand name] Select → Zone → 5-6 Å -
Delete or hide everything else:
Select → Invert Actions → Atoms/Bonds → Delete (or Hide) -
Show interactions:
Tools → Structure Analysis → FindHBond -
Measure specific distances:
Tools → Structure Analysis → Distances -
Look at electrostatics:
Actions → Surface → Show Tools → Surface/Binding Analysis → Coulombic Surface Coloring
What to Report
For protein-ligand interactions, describe:
| Interaction Type | How to Identify |
|---|---|
| Hydrogen bonds | Distance 2.5–3.5 Å, involves N-H or O-H |
| Salt bridges | Opposite charges, distance ~2.8–4 Å |
| Hydrophobic | Nonpolar residues surrounding nonpolar parts of ligand |
| Coordination | Metal ion with specific geometry |
| Electrostatic complementarity | Blue pocket for negative ligand (or vice versa) |
Example: Analyzing a Binding Site (3eeb)
This is the exercise you did!
The Setup
1. Fetch 3eeb
2. Delete chain B (Select → Chain → B, then Delete)
3. Show surface, color by electrostatics
Result: Blue (positive) binding pocket for the negative IHP (6 phosphates).
The Details
1. Hide surface
2. Select IHP, then Zone 6 Å
3. Invert selection, Delete
4. Show sidechains, keep ribbon
5. Measure distances
Result:
- Ser 169 OG ↔ IHP O34: ~2.8 Å = hydrogen bond (Ser donates H)
- His 55 NE2 ↔ IHP O22: ~2.9 Å = hydrogen bond (His donates H)
The Interpretation
"IHP binding is driven by electrostatic attraction (positive pocket, negative ligand) and stabilized by specific hydrogen bonds from Ser 169 and His 55 to phosphate oxygens."
Cancer Mutations in p53 (1tup)
Example from your lectures showing how to analyze mutation hotspots:
The Hotspot Residues
| Residue | Type | Role |
|---|---|---|
| R248 | Contact | Directly touches DNA |
| R273 | Contact | Directly touches DNA |
| R175 | Structural | Stabilizes DNA-binding loop |
| H179 | Structural | Stabilizes DNA-binding loop |
Analysis Approach
1. Open 1tup, keep chain B
2. Show R175, R248, R273, H179 in spacefill
3. Color surface by electrostatics
Result:
- R248 and R273 are right at the DNA interface (positive surface touching negative DNA)
- R175 and H179 are buried, maintaining the fold
- Mutations here → lose DNA binding → lose tumor suppression → cancer
Common Chimera Workflows
Quick Look at a Structure
1. File → Fetch by ID
2. Actions → Ribbon → Show
3. Presets → Interactive 1 (ribbons)
4. Rotate, zoom, explore
Analyze Active Site
1. Select ligand
2. Select → Zone → 5 Å
3. Actions → Atoms/Bonds → Show (for selection)
4. Tools → Structure Analysis → FindHBond
Compare Two Structures
1. Open both structures
2. Tools → Structure Comparison → MatchMaker
3. Check RMSD and sequence identity
Make a Figure
1. Set up your view
2. Presets → Publication 1
3. File → Save Image
Command Line Quick Reference
The command line is at the bottom of the Chimera window. Faster than menus once you know commands.
| Command | What it does |
|---|---|
open 1a6m | Fetch and open PDB |
select :153 | Select residue 153 |
select :HIS | Select all histidines |
select #0:.A | Select chain A of model 0 |
select :hem zr<5 | Select within 5 Å of heme |
display sel | Show selected atoms |
~display ~sel | Hide unselected atoms |
color red sel | Color selection red |
represent sphere | Spacefill for selection |
distance :169@OG :301@O34 | Measure distance |
match #1 #0 | Superpose model 1 onto 0 |
surface | Show surface |
~surface | Hide surface |
del sel | Delete selection |
Keyboard Shortcuts
| Key | Action |
|---|---|
| Ctrl + Click | Select atom |
| Ctrl + Shift + Click | Add to selection / measure distance |
| Scroll wheel | Zoom |
| Right-drag | Translate |
| Left-drag | Rotate |
| Middle-drag | Zoom (alternative) |
Troubleshooting Common Issues
"Nothing selected"
You tried to do something but nothing happened:
- Check: Is anything actually selected? (Green highlighting)
- Fix: Select → [what you want] first
Surface coloring disabled
- Check: Does a surface exist?
- Fix:
Actions → Surface → Showfirst
Can't see ligand
- Check: Is it hidden?
- Fix:
Select → Residue → [ligand], thenActions → Atoms/Bonds → Show
Structure looks weird after operations
- Fix:
Presets → Interactive 1to reset to default view
Atoms showing when you want ribbon only
Actions → Atoms/Bonds → Hide
Actions → Ribbon → Show
External Resources for Structure Analysis
| Resource | URL | Use for |
|---|---|---|
| RCSB PDB | rcsb.org | US PDB, structure info |
| PDBe | ebi.ac.uk/pdbe | European PDB, ligand interactions |
| PLIP | plip-tool.biotec.tu-dresden.de | Automated interaction analysis |
| AlphaFold DB | alphafold.ebi.ac.uk | Predicted structures |
| COSMIC | cancer.sanger.ac.uk/cosmic | Cancer mutations |
TL;DR
| Task | How |
|---|---|
| Open structure | File → Fetch by ID |
| Select | Select → [Chain/Residue/Zone] |
| Delete | Select, then Actions → Atoms/Bonds → Delete |
| Show surface | Actions → Surface → Show |
| Color by charge | Surface first, then Tools → Surface/Binding Analysis → Coulombic |
| Color by flexibility | Tools → Depiction → Render by Attribute → bfactor |
| Measure distance | Tools → Structure Analysis → Distances, then Ctrl+Shift+Click |
| Find H-bonds | Tools → Structure Analysis → FindHBond |
| Compare structures | Tools → Structure Comparison → MatchMaker |
Key distances:
- ~2.0 Å = coordination bond
- ~2.5–3.5 Å = hydrogen bond
- ~2.8–4.0 Å = salt bridge
Electrostatic colors:
- Blue = positive
- Red = negative
- White = neutral
Now go visualize some proteins! 🧬