Difference between revisions of "Normal mode analysis"
m (→All-atom NMA using Gromacs on the NOMAD-Ref server) |
|||
Line 39: | Line 39: | ||
Result: |
Result: |
||
{| class="centered" |
{| class="centered" |
||
− | |[[Image:Wnma_1a6z_7plot.png|thumb|Figure 1.1: Mode 7 by WEBnm@: |
+ | |[[Image:Wnma_1a6z_7plot.png|thumb|Figure 1.1: Mode 7 by WEBnm@: plot]] |
|[[Image:Wnma_1a6z_7.gif|thumb|Figure 1.2: Mode 7 by WEBnm@: vibrations]] |
|[[Image:Wnma_1a6z_7.gif|thumb|Figure 1.2: Mode 7 by WEBnm@: vibrations]] |
||
− | |[[Image:Wnma_1a6z_8plot.png|thumb|Figure 2.1: Mode 8 by WEBnm@: |
+ | |[[Image:Wnma_1a6z_8plot.png|thumb|Figure 2.1: Mode 8 by WEBnm@: plot]] |
|[[Image:Wnma_1a6z_8.gif|thumb|Figure 2.2: Mode 8 by WEBnm@: vibrations]] |
|[[Image:Wnma_1a6z_8.gif|thumb|Figure 2.2: Mode 8 by WEBnm@: vibrations]] |
||
+ | |[[Image:Wnma_1a6z_9plot.png|thumb|Figure 3.1: Mode 9 by WEBnm@: plot]] |
||
+ | |[[Image:Wnma_1a6z_9.gif|thumb|Figure 3.2: Mode 9 by WEBnm@: vibrations]] |
||
|} |
|} |
||
{| class="centered" |
{| class="centered" |
||
− | |[[Image: |
+ | |[[Image:Wnma_1a6z_10plot.png|thumb|Figure 4.1: Mode 10 by WEBnm@: plot]] |
− | |[[Image:Wnma_1a6z_9.gif|thumb|Figure 3.2: Mode 9 by WEBnm@: vibrations]] |
||
− | |[[Image:Wnma_1a6z_10plot.png|thumb|Figure 4.1: Mode 10 by WEBnm@: nsad]] |
||
|[[Image:Wnma_1a6z_10.gif|thumb|Figure 4.2: Mode 10 by WEBnm@: vibrations]] |
|[[Image:Wnma_1a6z_10.gif|thumb|Figure 4.2: Mode 10 by WEBnm@: vibrations]] |
||
+ | |[[Image:Wnma_1a6z_11plot.png|thumb|Figure 5.1: Mode 11 by WEBnm@: plot]] |
||
− | |} |
||
− | |||
− | {| class="centered" |
||
− | |[[Image:Wnma_1a6z_11plot.png|thumb|Figure 5.1: Mode 11 by WEBnm@: nsad]] |
||
|[[Image:Wnma_1a6z_11.gif|thumb|Figure 5.2: Mode 11 by WEBnm@: vibrations]] |
|[[Image:Wnma_1a6z_11.gif|thumb|Figure 5.2: Mode 11 by WEBnm@: vibrations]] |
||
− | |[[Image:Wnma_1a6z_12plot.png|thumb|Figure 6.1: Mode 12 by WEBnm@: |
+ | |[[Image:Wnma_1a6z_12plot.png|thumb|Figure 6.1: Mode 12 by WEBnm@: plot]] |
|[[Image:Wnma_1a6z_12.gif|thumb|Figure 6.2: Mode 12 by WEBnm@: vibrations]] |
|[[Image:Wnma_1a6z_12.gif|thumb|Figure 6.2: Mode 12 by WEBnm@: vibrations]] |
||
|} |
|} |
||
Line 123: | Line 120: | ||
|[[Image:Anm_Mode_3.gif|thumb|Figure 14: Mode 3 by ANM]] |
|[[Image:Anm_Mode_3.gif|thumb|Figure 14: Mode 3 by ANM]] |
||
|[[Image:Anm_Mode_4.gif|thumb|Figure 15: Mode 4 by ANM]] |
|[[Image:Anm_Mode_4.gif|thumb|Figure 15: Mode 4 by ANM]] |
||
− | |} |
||
− | |||
− | {| class="centered" |
||
|[[Image:Anm_Mode_5.gif|thumb|Figure 16: Mode 5 by ANM]] |
|[[Image:Anm_Mode_5.gif|thumb|Figure 16: Mode 5 by ANM]] |
||
|[[Image:Anm_Mode_6.gif|thumb|Figure 17: Mode 6 by ANM]] |
|[[Image:Anm_Mode_6.gif|thumb|Figure 17: Mode 6 by ANM]] |
||
Line 149: | Line 143: | ||
|[[Image:Ogn_mode_2_plot.png|thumb|Figure 19.1: Mode 2 by oGNM: plot]] |
|[[Image:Ogn_mode_2_plot.png|thumb|Figure 19.1: Mode 2 by oGNM: plot]] |
||
|[[Image:Ogn_mode_2.png|thumb|Figure 19.2: Mode 2 by oGNM: vibrations]] |
|[[Image:Ogn_mode_2.png|thumb|Figure 19.2: Mode 2 by oGNM: vibrations]] |
||
+ | |[[Image:Ogn_mode_3_plot.png|thumb|Figure 20.1: Mode 3 by oGNM: plot]] |
||
+ | |[[Image:Ogn_mode_3.png|thumb|Figure 20.2: Mode 3 by oGNM: vibrations]] |
||
|} |
|} |
||
{| class="centered" |
{| class="centered" |
||
− | |[[Image:Ogn_mode_3_plot.png|thumb|Figure 20.1: Mode 3 by oGNM: plot]] |
||
− | |[[Image:Ogn_mode_3.png|thumb|Figure 20.2: Mode 3 by oGNM: vibrations]] |
||
|[[Image:Ogn_mode_4_plot.png|thumb|Figure 21.1: Mode 4 by oGNM: plot]] |
|[[Image:Ogn_mode_4_plot.png|thumb|Figure 21.1: Mode 4 by oGNM: plot]] |
||
|[[Image:Ogn_mode_4.png|thumb|Figure 21.2: Mode 4 by oGNM: vibrations]] |
|[[Image:Ogn_mode_4.png|thumb|Figure 21.2: Mode 4 by oGNM: vibrations]] |
||
− | |} |
||
− | |||
− | {| class="centered" |
||
|[[Image:Ogn_mode_5_plot.png|thumb|Figure 22.1: Mode 5 by oGNM: plot]] |
|[[Image:Ogn_mode_5_plot.png|thumb|Figure 22.1: Mode 5 by oGNM: plot]] |
||
|[[Image:Ogn_mode_5.png|thumb|Figure 22.2: Mode 5 by oGNM: vibrations]] |
|[[Image:Ogn_mode_5.png|thumb|Figure 22.2: Mode 5 by oGNM: vibrations]] |
||
Line 177: | Line 168: | ||
|[[Image:nomadref_mode_7_0.gif|thumb|Figure 24: Mode 7 by ANM]] |
|[[Image:nomadref_mode_7_0.gif|thumb|Figure 24: Mode 7 by ANM]] |
||
|[[Image:nomadref_mode_8_0.gif|thumb|Figure 25: Mode 8 by ANM]] |
|[[Image:nomadref_mode_8_0.gif|thumb|Figure 25: Mode 8 by ANM]] |
||
− | |} |
||
− | |||
− | {| class="centered" |
||
|[[Image:nomadref_mode_9_0.gif|thumb|Figure 26: Mode 9 by ANM]] |
|[[Image:nomadref_mode_9_0.gif|thumb|Figure 26: Mode 9 by ANM]] |
||
|[[Image:nomadref_mode_10_0.gif|thumb|Figure 27: Mode 10 by ANM]] |
|[[Image:nomadref_mode_10_0.gif|thumb|Figure 27: Mode 10 by ANM]] |
||
− | |} |
||
− | |||
− | {| class="centered" |
||
|[[Image:nomadref_mode_11_0.gif|thumb|Figure 28: Mode 11 by ANM]] |
|[[Image:nomadref_mode_11_0.gif|thumb|Figure 28: Mode 11 by ANM]] |
||
|[[Image:nomadref_mode_12_0.gif|thumb|Figure 29: Mode 12 by ANM]] |
|[[Image:nomadref_mode_12_0.gif|thumb|Figure 29: Mode 12 by ANM]] |
Revision as of 16:52, 24 August 2011
Contents
TODO
WORKY!
re-check all pictures
re-reference all pictures
add all references/quotes
Introduction
NMA (normal mode analysis) is a time-independent apprach to simulate low-frequency motions and vibrations of protein. These simulation are all based on the harmonic approximation and therefore ignore the influence of the solvent. The proteins are seen as models made out of springs and point masses, which are connected and represent the interatomic forces. Simulation done this way are very easy to do, but are no more than a slight insight into the protein flexibility.
Task?
For each server, analyze at least the lowest five normal modes.
- What information do the different servers provide?
- Which regions of your protein are most flexible, most stable?
- When you visualize the modes (provided by server or using for example PyMol or VMD), try to describe what movements you observe? Hinge-movement, “breathing”…
- Can you observe notable differences between the normal modes calculated by different servers?
- Out of the servers, chose one or two favorites and discuss the results of these in more detail. Why do you like these?
- When your MD simulations are finished, compare the lowest-frequency normal modes with your MD simulation using visualization software, e.g. PyMol or VMD. Can you observe different movements or similar dynamics? If possible, compare an overlay of the lowest-frequency modes to your MD simulation. You can superimpose the normal modes for example in VMD.
- What are the advantages and disadvantages of NMA compared to MD?
Parameters
- If possible, use a cutoff for Cα atom pairs of 15 Å.
- Calculate the 10 lowest-frequency normal modes (the six zero modes have to be considered for a few applications).
- In most cases you can upload the original .pdb file from the Protein Data Bank. In some cases, however, you can upload only the structure itself (ATOM lines of the .pdb file).
WEBnm@
WEBnm@ is a webserver based application that allows computation and low-frequency analysis of normal nodes of proteins. This computation is fully automated and only different types of results are presented to the user.
Webserver:
Input:
- 1a6z - all chains
Result:
Discussion:
All animated gifs have to be created the hard way, frame after frame, because WEBnm@ does not allow the concurrent saving of more than one frame.
The Normalized Squared Atomic Displacements (nsad) plots show the vibrations according to the amino acid position.
Except for mode 11 there is no special movement inside the alpha helix of chains A and C. The movement is almost everytime between the chains or inside/around the beta strands of chain B and D. This behaviour is also visible by analyzing the plots; the regions of low movement are always around the chains A and C with their corresponding alpha helices and the high movement regions lies within the beta strands of chain B and D.
The movement/vibrations can be described mostly as repulsive or flattening, stretching and twisting.
There seems to be some strange behaviour at figure 4.2 mode 10; it is slighty twitching and we do not know why. Maybe it is because of a wrong frame or some other aspect of visual glitches, we will check that again, if there is time.
ElNemo
ElNémo is a webserver based to work with the Elastic Network Model. It calculates and analyses low-frequency normal modes of proteins.
Webserver:
Input:
- 1a6z
Result:
Discussion:
For all generated models the vibrations are shown in three different perspectives.
The Movement/Vibrations are very similar to these obtained by WEBnm@. There is almost no movement inside the alpha helices of chain A and C and much movement inside and outside the the beta strands of chain B and D. Vibrations between chains can also be observed but these are mostly between A+C and B+D because they form a subunit.
Anisotropic Network Model web server
The Anisotropic Network Model web server uses the fast approach anisotropic network model (elastic network) to calculate the global modes.
Webserver:
Params:
- distance weight: 3
Result:
Discussion:
oGNM – Gaussian network model
Webserver:
Input:
- 1a6z
Params:
- cutoff: 15 Å
Result:
NOMAD-Ref
Webserver:
Params:
- distance weight: 3.0
- cutoff: 15 Å
Result:
Discussion:
All-atom NMA using Gromacs on the NOMAD-Ref server
Webserver:
Params:
- temperature: 600K and 2000K
- pdb ID: 1BPT
Information: We used the given protein 1BPT because our HFE protein (1A6Z) has around 6080 ATOM lines, which are 4080 too much.
Result:
- at 600K
- at 2000K
- with Elastic Network
Discussion: