@*~RSC_HELP_TUTORIAL_MUTATIONS
Introduction
Visualization of residue positions in 3D structures is an
important step towards the interpretation of
mutations or polymorphisms in terms of protein function,
interaction and thermodynamic stability.
It helps to identify phenotype - genotype correlations.
For an increasing number of proteins, that may be analyzed and
screened for mutations, three-dimensional protein structures can be
found in the PDB-database.
Selecting and highlighting large numbers of residue positions in a
protein structure can be time-consuming and tedious and is error-prone when performed by
hand.
Strap facilitates mapping of mutations onto
three-dimensional protein structures.
This tutorial explains how all known point mutations of the human cardiac beta WIKI:myosin
are visualized on the 3D-structure of an homologous protein.
Mutations in this gene usually cause WIKI:Hypertrophic_cardiomyopathy.
Only a very few mutations produce a completely different heart disease, WIKI:Dilated_cardiomyopathy.
Also see publication PUBMED:16322575.
With the techniques described here it is possible to compare the 3D
location of both groups of mutations.
Overview of the program modules used in this tutorial
Strap offers several program modules which are involved.
-
If the designations of the mutations refer to sites in genomic DNA, Strap
translates the sequence file from DNA to amino acids:
HTMLDOC_DIALOG:BUT_C1(DialogGenbank). The reading orientation and the intron/exon boundaries can be altered manually:
HTMLDOC_DIALOG:BUT_C1(DialogEditDna)
-
If the structure of the protein of interest is not yet available, a
related protein can frequently be found in the structure databases
using
HTMLDOC_DIALOG:BUT_C1(DialogBlast).
In this case the alignment of both proteins becomes the
crucial part of the analysis:
HTMLDOC_DIALOG:BUT_C1(DialogAlign)
-
The mutations found in public mutation databases or mutations identified in the own lab
are imported as a text list of mutation designations.
HTMLDOC_DIALOG:BUT_C1(DialogResidueAnnotationList)
-
The structure is viewed in Astex or another 3D-viewer:
HTMLDOC_DIALOG:BUT_C1(Dialog3DViewer)
Tutorial Start
Download human beta cardiac myosin M57965 from the Genbank ( GB:M57965 ) or from EMBL ( EMBL:M57965 ) either by clicking these links or using
HTMLDOC_BUTTON:BUT_C1(DialogFetchSRS)!.
Exercise:
- View the protein file and find the CDS-text string describing the positions of exons.
Use HTMLDOC_BUTTON:M_SEL_EDIT from the context menu.
Translate the Genbank or Embl file
The file M57965 contains a nucleotide sequence.
The button HTMLDOC_BUTTON:SBUT_NUCL! above the alignment becomes active.
The nucleotide sequence is translated into amino acids
using the CDS-expression join(5721..5921,6225..6368,6657..6813,6945.. ) contained in this file.
In this case there is only one gene record and Strap initially translates the nucleotides into amino acids.
Import the mutations or SNPs
Several online resources contain lists of mutations for various genes and diseases.
Open the dialog HTMLDOC_BUTTON:BUT_C1(DialogResidueAnnotationList)!
Select the protein M57965 in the dialog and copy the following list of mutations into the multi-line text field of the dialog:
R054X V059I T124I R143Q Y162C N187K R190T Q222K N232S
F244L K246Q R249Q G256E I263T A326P K383N L390V
R403L R403Q R403W V406M R453C R453H R453L E483K Q499K
F513C G584R D587V L601V N602S V606M K615N
R663C R663H R694C N696S R712L G716R R719Q R719W R723C R723G P731L
I736M G741A G741R G741W D778G S782N A797T E846K E846Q R869C
R870C R870H M877K L908V E924K E930K E930del E935K E949K L961R
F764L S532P
Then press HTMLDOC_BUTTON:BUTS_GO to generate the residue selection in M57965 for each mutation.
These mutations had been obtained from http://genetics.med.harvard.edu/~seidman/cg3/
and http://www.angis.org.au/pbin/Databases/Heart/fhcquery.cgi.
Other resources are for example.:
- http://www.hgvs.org/
- https://www.ncbi.nlm.nih.gov/SNP/
- WIKI:OMIA Online Mendelian Inheritance in Animals other than human and mouse
- WIKI:OMIM Online Mendelian Inheritance in Man
- WIKI:SNP Single nucleotide polymorphism.
- http://www.umd.be/
These residue selection have a context menu (right-click) and a Balloon message.
When HTMLDOC_ITEM:M_SEL_EDIT in the context menu is clicked, a table is opened which shows the name, sequence positions
and additional annotations.
Exercise:
- Choose one mutation and change the color.
- Edit one mutation and add comments.
- Simultaneously select two mutations and change their color. Two mutations can be selected in the alignment with the Ctrl-key.
- Select several mutations by dragging or with the rubber band.
Referring to nucleotide sequence positions
The dialog HTMLDOC_BUTTON:BUT_C1(DialogResidueAnnotationList)! can also be used when positions refer to nucleotides instead of amino acids.
In this case "nucleotides" must be selected in the choice menu.
Please add the following two mutations:
C12164T G7799A
You might use a different color for those.
Since these mutations are defined by nucleotide positions they can be
best seen in the nucleotide view.
Select M57965 and click the menu item HTMLDOC_BUTTON:BUT_C1(DialogEditDna)! or the tool-button HTMLDOC_BUTTON:SBUT_NUCL
.
The entire nucleotide sequence in the file is very large (26689 bp) and navigation with the scroll-bar is relatively
difficult. This is because the coding exons are relatively short compared to the non-coding introns.
The arrow keys allows to jump from one exon/intron start to the next.
View DNA sequence in Alignment Pane
Nucleotides can be shown in the alignment pane simultaneously with amino acids.
Select the last item "nucleotide" in
HTMLDOC_JCOMPONENT:HTMLDOC_JC_Strap_comboColorScheme!.
In M57965 each amino acid corresponds to three nucleotides.
The nucleotides are shown as colored bars, guanine=GRAY, adenine=green, thymine=red and cytosine=blue.
Exercise:
- Identify the sequence position of the mutations C12164T and G7799A in the alignment pane.
Loading the 3D-structure file
Because a three-dimensional structure for the cardiac protein is not yet available, we use
the structure of myosin of skeletal muscle instead.
PDB:1ALM is a theoretical model of the relative position of the backbones of myosin and five actin monomers.
The myosin heavy chain is chain A in 1ALM and the actins are chains V,
W, X, Y and Z.
The X-ray structure 2MYS, however does not contain actin but contains the amino acid side chains.
Therefore both are combined to one file for this tutorial and can
be loaded by clicking
HTMLDOC_BUTTON:SBUT_TUTORIAL_LOAD_EXAMPLE_MUT!.
Aligning myosin from skeletal muscle and heart muscle
Before copying the mutations from M57965 to 1alm_2mys, both must be aligned.
Select both (Ctrl+click) and push the tool-button "Align 2 proteins".
3D-visualization
In order to visualize the protein 3D trace, mark all alignment rows with
a pdb structure file and open the 3D-wire representation
HTMLDOC_BUTTON:SBUT_3D!.
Full featured 3D view
Load the protein into the Astex 3D viewer by clicking
HTMLDOC_BUTTON:M_PV_DEFAULT_VIEWER above the
protein backbone view.
You can drag residue annotations such as mutations directly into Astex, even if they are placed on a different protein.
See below for how to transfer mutations.
Transfer the mutations to 3D
Currently the mutations are shown in M57965 but we want to see
them in the protein structure file "1alm_2mys_A.ent".
First select the residue annotations of M57965.
Selected annotations are indicated by WIKI:Marching_ants.
For large numbers, clicking each selection would be tedious and there are various
approaches which are more efficient.
-
Open a red/white rubber band by dragging the mouse.
Initially, the rectangle must comprise more than one row but can be narrowed afterwards.
All residue selections in the rectangle are selected.
- The menu "Select" offers a string pattern search.
-
The graphical tree lists all proteins and their associated objects. The tree is initially not shown, because the vertical
divider is at the left edge of the application frame.
-
Ctrl+click the first residue annotation and Shift+click the last residue annotation.
The selected residue annotations are indicated by marching ants and can be dragged to the 3D-structure.
#if CPP_WITH_PYMOL
Using Pymol instead if Astex
If you are used to WIKI:Pymol rather then Astex you may want to use Pymol.
It can be opened in HTMLDOC_DIALOG:BUT_C1(Dialog3DViewer).
Pymol will initially not display the actins since pdb1alm is a
model containing only C-alpha atoms but no amino acid side chains.
To make the C-alpha atoms visible, select "non-bonded" in
Pymol. This is done in the popup menus "[S]" of the Pymol side pane.
Since Pymol is not a Java application, you cannot drag the mutations onto the Pymol panel.
Therefore there is a Pymol button which brings up the 3D generic menu bar in Strap and which acts as a drop target.
#endif //CPP_WITH_PYMOL
*@