Hemochromatosis Normal modes

From Bioinformatikpedia
Revision as of 21:27, 1 September 2013 by Betza (talk | contribs) (WebNM@)

Labjournal_NMA-Hemochromatosis

Introduction

<figure id="timescale">

Figure 1: Timescale of protein motions <ref name="timescale_motions"> Normal Mode Analysis of Biomolecular Structures: Functional Mechanisms of Membrane Proteins Ivet Bahar, Timothy R. Lezon, Ahmet Bakan, and Indira H. Shrivastava Chemical Reviews 2010 110 (3), 1463-1497</ref>

</figure>

TODO for what do we use NMA? -> reference to <xr id="timescale"/>

WebNMA and elNemo are two different web servers for normal mode analysis. Both provide different information:

  • WebNMA
    • normal modes + animations
    • average residue fluctuations
    • correlated residue movement map
  • elNemo
    • normal modes + animations
    • average residue fluctuations
    • pairwise correlated residue movement as a list

The modes are only calculated from the C-alpha atoms of the proteins residues. In principal, 3N-6 normal modes could be computed, where N is the number of atoms in the molecule. The reasoning behind this number is, that each atom in a molecule has three coordinates and the molecule thus has 3N degrees of freedom. However, three degrees of freedom have to be omitted to account for translation and another three have to be omitted to account for rotational symmetry. For the HFE protein with 272 residues, 810 normal modes can be calculated, but only the low frequency normal modes are of interest, because they resemble slow motions.


elNemo

<figure id="elnemo">

Mode 7 Mode 8 Mode 11 Mode 12 Mode 13
Mode7 hfe.gif Mode8 hfe.gif Mode11 hfe.gif Mode12 hfe.gif Mode13.gif
This mode shows a hinge motion about the connection between the two domains. The immunoglobulin domain and the MHC-like domain move towards each other. This mode also shows a hinge motion about the domain connection, but in this case, the Ig domain twists away from the upward moving MHC-like domain. A hinge motion of the Ig domain and a breathing motion of the MHC-like domain can be observed in mode 11. Mode 12 and 13 both show a slight hinge motion of the Ig domain and a strong breathing motion of the MHC-like domain. This observation can be expected, because the MHC domain is part of the binding interface to the ferritin receptor, thus this interface needs to be flexible.
Figure 2: Visualisation of the five most interesting and distinct normal modes of the HFE protein computed with elNemo. The immunoglobulin (Ig) domain is colored in yellow-red and the MHC-like domain in blue-green.

</figure>

<xr id="elnemo"/> shows the five most most interesting and distinct normal modes of the HFE protein. The first six are omitted because they are of low complexity.

WebNM@

<css> table.basic2 { margin-left: auto; margin-right: auto; border: 0px solid black; border-collapse:collapse; width: 20%; }

.basic2 th,td { padding: 3px; border: 1px solid black; }

.basic2 th { border-bottom: 2px solid black; background-color: #fff; }

.basic2 td { text-align:left; }

.basic2 tr:first-child th {

 border-top: 0;

} .basic2 tr:last-child td {

 border-bottom: 0;

} .basic2 tr td:first-child, .basic2 tr th:first-child {

 border-left: 0;

} .basic2 tr td:last-child, .basic2 tr th:last-child {

 border-right: 0;

}

</css>

<figtable id="energy">

Mode Deformation Energy
7 325.91
8 817.13
9 849.13
10 1452.95
11 1876.97
12 2972.27
Table 1: Deformation energy of the six first normal modes computed with webNMA.

</figtable>

Apart from elNemo, we also used webNMA to compute the normal modes of the HFE protein. The deformation energy of the six first normal modes is given in <xr id="energy"/>.

TODO: explanation of deformation energy, what does it tell us?

<figtable id="webnma_modes7to12">

Mode 7
Mode 8
Mode 9
Mode 10
Mode 11
Mode 12
Table 2: Visualization of the normal modes 7 to 12 for 1A6Z chain A computed with WebNMA. The original protein structure is shown in green and the movement of the C alpha atoms is shown as a red ribbon.

</figtable>

The movements of the normal modes 7 to 12 are visualised in <xr id="webnma_modes7to12"/>.

Domain definition

<figure id="corr_motions">

Figure 3: Correlated residue motions for the HFE protein

</figure>

<xr id="corr_motions"/> shows two distinct zones with correlated movement. The first zone starts at the beginning of the sequence and ends at residue 175. The second zone starts at residue 176 and extends until the end of the sequence. Based on this observation, two domains can be assigned to the two zones. SCOP, CATH and Pfam all assign the same two domains to the HFE protein. <xr id="1a6z_domains" /> shows a 3D visualisation of the HFE protein and the MHC I domain (green) and the immunoglobulin domain (blue). This is in complete accordance with the domain assignment based on the correlated motion matrix.

References:

CATH MHC domain CATH Ig domain

SCOP MHC domain SCOP Ig domain

Pfam domains

<figure id="1a6z_domains">

Side view
Front view
Top view
Figure 4: Pymol visualisation of the two domains of the HFE protein from different perspectives. The MHC I domain is colored in green and the immunoglobulin domain in blue.

</figure>

Flexibility

<figure id="fluctuations">

Figure 5: The square of the fluctuation of each C-alpha atom (for all non-trivial modes), normalized so that the sum over all residues is equal to 100. The fluctuations are the sum of the atomic displacements in each mode weighted by the inverse of their corresponding eigenvalues.

</figure>

<figure id="3d_fluc">

Side
Front
Figure 6: The flexibility of the C-alpha atoms of the HFE protein is shown color coded, where red indicates no flexibility, white indicates medium flexibility and blue indicates high flexibility.

</figure>

Figure <xr id="fluctuations"/> shows some highly flexible regions connected by rigid stretches. Another way to look at the C-alpha fluctuation data is shown in figure <xr id="3d_fluc"/>. The more flexible regions are almost exclusively located at the ends of the protein, that are farthest away from the hinge, that is the stretch, that connects the two domains.

Comparison of Molecular Dynamics and Normal Mode Analysis

In comparison to Molecular Dynamics, normal mode analysis can grasp the slow and globular motions of a protein. Although it can only grasp the harmonic part of these motions, it gives a good insight into them at comparably low computational cost. MD takes the influence of the side chains into account by modeling them explicitly, whereas NMA only takes them into account implicitly through the position of the C-alphas. Thus, MD can model the effect of side chain movements and SNPs on the proteins dynamics, whereas NMA cannot do this.

References

<references />