Difference between revisions of "ASPA Normal Mode Analysis"

From Bioinformatikpedia
(NOMAD-ref)
(NOMAD-ref)
Line 137: Line 137:
 
|[[File:canavan_nomadref_mode7_anim.gif |thumb|center|300px| Nomad-ref Normal Mode 7]]
 
|[[File:canavan_nomadref_mode7_anim.gif |thumb|center|300px| Nomad-ref Normal Mode 7]]
 
|[[File:Mode_7.gif |thumb|center|300px| cRMS plot for mode 7]]
 
|[[File:Mode_7.gif |thumb|center|300px| cRMS plot for mode 7]]
  +
|}
  +
  +
{| cellpadding="2" style="border: 0px;" align="center"
 
|[[File:canavan_nomadref_mode8_anim.gif |thumb|center|300px| Nomad-ref Normal Mode 8]]
 
|[[File:canavan_nomadref_mode8_anim.gif |thumb|center|300px| Nomad-ref Normal Mode 8]]
 
|[[File:Mode_8.gif |thumb|center|300px| cRMS plot for mode 8]]
 
|[[File:Mode_8.gif |thumb|center|300px| cRMS plot for mode 8]]
Line 146: Line 149:
 
|[[File:canavan_nomadref_mode9_anim.gif |thumb|center|300px| Nomad-ref Normal Mode 9]]
 
|[[File:canavan_nomadref_mode9_anim.gif |thumb|center|300px| Nomad-ref Normal Mode 9]]
 
|[[File:Mode_9.gif |thumb|center|300px| cRMS plot for mode 9]]
 
|[[File:Mode_9.gif |thumb|center|300px| cRMS plot for mode 9]]
  +
|}
  +
  +
  +
  +
{| cellpadding="2" style="border: 0px;" align="center"
 
|[[File:canavan_nomadref_mode10_anim.gif |thumb|center|300px| Nomad-ref Normal Mode 10]]
 
|[[File:canavan_nomadref_mode10_anim.gif |thumb|center|300px| Nomad-ref Normal Mode 10]]
 
|[[File:Mode_10.gif |thumb|center|300px| cRMS plot for mode 10]]
 
|[[File:Mode_10.gif |thumb|center|300px| cRMS plot for mode 10]]
 
|}
 
|}
  +
  +
   
 
{| cellpadding="2" style="border: 0px;" align="center"
 
{| cellpadding="2" style="border: 0px;" align="center"
 
|[[File:canavan_nomadref_mode11_anim.gif |thumb|center|300px| Nomad-ref Normal Mode 11]]
 
|[[File:canavan_nomadref_mode11_anim.gif |thumb|center|300px| Nomad-ref Normal Mode 11]]
 
|[[File:Mode_11.gif |thumb|center|300px| cRMS plot for mode 11]]
 
|[[File:Mode_11.gif |thumb|center|300px| cRMS plot for mode 11]]
  +
|}
  +
  +
  +
  +
{| cellpadding="2" style="border: 0px;" align="center"
 
|[[File:canavan_nomadref_mode12_anim.gif |thumb|center|300px| Nomad-ref Normal Mode 12]]
 
|[[File:canavan_nomadref_mode12_anim.gif |thumb|center|300px| Nomad-ref Normal Mode 12]]
 
|[[File:Mode_12.gif |thumb|center|300px| cRMS plot for mode 12]]
 
|[[File:Mode_12.gif |thumb|center|300px| cRMS plot for mode 12]]

Revision as of 21:16, 31 August 2011

WebNM@

WebNM@ is a Normal Mode analysis tool developed by Hollup, Sælensminde and Reuter. It is one of the methods uses that only look at C_Alpha atoms. We ran it on our protein (PDB Code 2O53) using the default settings and using all chains.

These are the deformation energies computed by WebNM@:

Mode Index 	Deformation Energy 	Mode Index 	Deformation Energy
7 	         870.03 	        14              3740.78
8 	         976.54 	        15              4168.69
9 	        1562.41 	        16              5263.58
10 	        2800.99 	        17      	5461.10
11 	        2798.65 	        18 	        5906.92
12 	        2954.75 	        19 	        6126.73
13 	        3448.49 	        20 	        6356.04

Since the actual energy of the protein is a composite of all the normal modes (after all, all of the motions happen simultaneously), comparing the energies is of little significance. What does meet the eye, however, is that there is a fairly steady increase of energy with higher normal mode numbers.

WebNM@ also generates atomic displacement charts. These indicate how much movement individual residues undergo in the various normal modes.

ASPA webnma displacement.png

These are quite interesting in that they differ so strongly, yet there are quite a few places (around residue 100, 200) where almost all of the graphs show peaks to some degree. These seem to be regions that are very flexible. Mode 7 shows a rather strong cluster of peaks in between 350-400. This cluster is recognizable in Mode 8 and Mode 10 too, and to some degree in Mode 12. This seems to indicate a highly mobile region of the protein. We expected some kind of faint symmetry to be visible since we analyzed both domains, but there is none to be seen in the plots, which leads us to assume that this hypothesis is faulty. We also expected to see low mobility around the active site at 117 resp. 427; we found that this supposition holds.

Also returned was a graph of the fluctuations in position of the individual C Alpha atoms:

ASPA webnma fluctuations.png

There seem to be two mobile regions, around ca. 80 and ca. 380, and a couple of more isolated peaks that we assume to denote local mobility e.g. of coils. It is interesting to see that the active site around 117 resp. 427 is in a region of very low mobility, as we would have expected to find.

Finally, animations of the individual movements for the Normal Modes:

WebNM@ Normal Mode 7
WebNM@ Normal Mode 8
WebNM@ Normal Mode 9
WebNM@ Normal Mode 10
WebNM@ Normal Mode 11
WebNM@ Normal Mode 12

ElNemo

The ElNemo web server was published in 2004 by Suhre and Sanejouand. It uses an elastic network model to compute a large number of normal modes and produce animations for the first five non-trivial normal modes for a given protein structure. We fed our protein, 2O53, into the ElNemo webserver and ran a calculation with default parameters.


ElNemo Normal Mode 7
ElNemo Fluctuations for Normal Mode 7
ElNemo Normal Mode 8
ElNemo Fluctuations for Normal Mode 8

In normal mode 7 we observe a shearing motion of the two domains against each other, roughly around an axis that transverses both domains and the connection between them. The Fluctuation Chart gives two big red areas; these denote amino acid pairs that undergo significant decreases in distance. The numerous blue spots denote increases in distance.

In normal mode 8, again a shearing motion is visible, albeit around a different axis this time, seemingly rotated somewhat towards the viewer. The fluctuation chart shows no such distinct patterns as the normal mode 7 one, instead a larger number of small splotches can be seen. There are more distance increases and less decreases, too.

ElNemo Normal Mode 9
ElNemo Fluctuations for Normal Mode 9
ElNemo Normal Mode 10
ElNemo Fluctuations for Normal Mode 10

Normal mode 9 exhibits a jaw-like or pliers-like open-and-close movement. The fluctuation graph contains no large areas of either distinct overall increase or overall decrease of distances; instead, a number of smaller patches of either can be seen.

Normal mode 10 displays an expanding motion. The fluctuation chart has odd cross-like red segments and a few clustered blue patches.

ElNemo Normal Mode 11
ElNemo Fluctuations for Normal Mode 11

Normal mode 11 appears to undergo two different simultaneous shearing motions. The fluctuation graph is largely composed of red patches of varying size.

Anisotropic Network Model web server

In Normal Mode 1, there are two axes around which rotation happens. Both the parallel axes go through one of the two domains, and the domains undergo aboud 30 degrees of rotation around them, being constrained by the inter-domain binding site. The motion gives the impression of closing and opening jaws or like the head of a pair of pliers.

The distance matrix shows strong correlation within two distinct parts of the structure, visible in the form of two big square red regions. These are the domains that move independently from each other.

The deformation energy graph shows one very mobile cluster at around 60 and little in the way of mobility around the active site at 117.


In the case of the second Normal Mode 2 there is but one axis of rotation which transverses both domains and their connecting site. The domains undergo around 30 degrees of rotation in opposing directions, again constrained by the binding site. The motion makes a sort-of dancing impression.

The distance matrix shows weaker correlation, though the two domains are still discernible as two mostly-red square shapes. There is much more correlation in the rest of the graph than in Normal Mode 1, though.

The deformation energy graph again shows little deformation near the active site and again a mobile cluster at around 60. It is less pronounced than in Normal Mode 1, but still quite prominent.

In Normal Mode 3 the motion that occurs is strikingly similar to that in Normal Mode 1, except that the direction of motion slightly differs.

The distance matrix for Normal Mode 3 contains two smaller but well-defined rectangles that correspond to regions within the domains that have strongly correlated motion.

In the deformation energy graph, a number of clustered peaks is visible on the right side of the graph; the 60 peak cluster is largely gone, and the active site still shows little in the way of deformation.

In Normal Mode 4, both domains undergo a clamshell-like open/close motion.

The distance matrix is quite uniformly correlated. There are no readily discernible areas of circumscribed correlation, but a few small areas of anticorrelation are visible. These likely correspond to the edges of the clamshell-like cavity.

The deformation energy graph, a strong cluster of deformation is visible around 60. Another smaller cluster can be seen around 280, and at the active site little to no deformation happens. One isolated peak at 160 probably corresponds to a small isolated portion of a coil section that is highly mobile.

In Normal Mode 5 we observe a shearing motion with clamshell-like characteristics.

Again, the distance matrix shows larges areas of strong correlation with interspersed small areas of strong anticorrelation. We again surmise that the latter parts of the plot correspond to the clamshell cavity which undergoes strongly anticorrelated motion.

The deformation graph shows little deformation around the active site at 117, and no really strong clusters of motion. There are some weaker clusters at around 70, 160 and 280.

ANM web server Normal Mode 1
ANM web server Normal Mode 1 Distance Matrix
ANM web server Normal Mode 1 Deformation Energy
ANM web server Normal Mode 2
ANM web server Normal Mode 2 Distance Matrix
ANM web server Normal Mode 2 Deformation Energy
ANM web server Normal Mode 3
ANM web server Normal Mode 3 Distance Matrix
ANM web server Normal Mode 3 Deformation Energy
ANM web server Normal Mode 4
ANM web server Normal Mode 4 Distance Matrix
ANM web server Normal Mode 4 Deformation Energy
ANM web server Normal Mode 5
ANM web server Normal Mode 5 Distance Matrix
ANM web server Normal Mode 5 Deformation Energy

NOMAD-ref

Nomad-ref Normal Mode 7
cRMS plot for mode 7
Nomad-ref Normal Mode 8
cRMS plot for mode 8

In Normal Mode 7 we observe a jaw-like open-close motion; the individual domains stay fairly rigid, and the motion occurs at a number of connecting coils which act as a sort-of hinge.

Nomad-ref Normal Mode 9
cRMS plot for mode 9


Nomad-ref Normal Mode 10
cRMS plot for mode 10


Nomad-ref Normal Mode 11
cRMS plot for mode 11


Nomad-ref Normal Mode 12
cRMS plot for mode 12

All-Atom NMA

Temperature 600K

Nomad-ref Normal Mode 7
Nomad-ref Normal Mode 8
Nomad-ref Normal Mode 9
Nomad-ref Normal Mode 10
Nomad-ref Normal Mode 11
Nomad-ref Normal Mode 12

Temperature 2000K

Nomad-ref Normal Mode 7
Nomad-ref Normal Mode 8
Nomad-ref Normal Mode 9
Nomad-ref Normal Mode 10
Nomad-ref Normal Mode 11
Nomad-ref Normal Mode 12