Difference between revisions of "Normal mode analysis GLA"

From Bioinformatikpedia
(Discussion)
(BPTI)
Line 157: Line 157:
   
 
==BPTI==
 
==BPTI==
  +
{| class="wikitable" border="1" style="text-align:center; border-spacing:0;"
  +
|-
  +
! Mode || NOMAD All Atom (600K) || NOMAD All Atom (2000K) || NOMAD Ref
  +
|-
  +
| Mode 7 || [[File:BPTI_nomad_atom_600K_7.gif|thumb|center|185px|]] || [[File:BPTI_nomad_atom_2000K_7.gif|thumb|center|185px|]] || [[File:BPTI_nomad_ref_7.gif|thumb|center|185px|]]
  +
|-
  +
| Mode 7 || [[File:BPTI_nomad_atom_600K_8.gif|thumb|center|185px|]] || [[File:BPTI_nomad_atom_2000K_8.gif|thumb|center|185px|]] || [[File:BPTI_nomad_ref_8.gif|thumb|center|185px|]]
  +
|-
  +
| Mode 7 || [[File:BPTI_nomad_atom_600K_9.gif|thumb|center|185px|]] || [[File:BPTI_nomad_atom_2000K_9.gif|thumb|center|185px|]] || [[File:BPTI_nomad_ref_9.gif|thumb|center|185px|]]
  +
|-
  +
| Mode 7 || [[File:BPTI_nomad_atom_600K_10.gif|thumb|center|185px|]] || [[File:BPTI_nomad_atom_2000K_10.gif|thumb|center|185px|]] || [[File:BPTI_nomad_ref_10.gif|thumb|center|185px|]]
  +
|-
  +
| Mode 7 || [[File:BPTI_nomad_atom_600K_11.gif|thumb|center|185px|]] || [[File:BPTI_nomad_atom_2000K_11.gif|thumb|center|185px|]] || [[File:BPTI_nomad_ref_11.gif|thumb|center|185px|]]
  +
|-
  +
| Mode 7 || [[File:BPTI_nomad_atom_600K_12.gif|thumb|center|185px|]] || [[File:BPTI_nomad_atom_2000K_12.gif|thumb|center|185px|]] || [[File:BPTI_nomad_ref_12.gif|thumb|center|185px|]]
  +
|}
   
 
=Discussion=
 
=Discussion=

Revision as of 01:49, 2 September 2011

by Benjamin Drexler and Fabian Grandke

Introduction

Methods

The following sections give a brief explanation about the methods and information about their usage.

WEBnm@

The normal mode analysis webserver WEBnm@ was published by Siv Midtun Hollup, Gisle Salensminde and Nathalie Reuter in 2005<ref name=webnma>Siv Midtun Hollup, Gisle Salensminde and Nathalie Reuter. "WEBnm@: a web application for normal mode analyses of proteins". BMC Bioinformatics 2005, 6:52. PubMed</ref>. It allows the calculation of the normal modes and offers serveral types of analyses, i.e. deformation energy, animation of the vibration, atomic squared displacements and vector field analysis.


Usage

  • Webserver: http://apps.cbu.uib.no/webnma/home
  • Input
    • PDB structure (ID or file)
    • Specification of chains
  • Output
    • Animation of vibration
    • Deformation energy
    • Squared atomic displacements

ElNemo

Karsten Suhre and Yves-Henri Sanejouand published the webserver ElNemo in 2004<ref name=elnemo>Karsten Suhre and Yves-Henri Sanejouand. "ElNémo: a normal mode web server for protein movement analysis and the generation of templates for molecular replacement". Nucleic Acids Res. 2004 Jul 1;32(Web Server issue):W610-4. PubMed</ref>. It is used for the calculation of normal modes and is able to process very large proteins due to a building block approximation. This is, several residues are grouped into a single super residue.


Usage

  • Webserver: http://www.igs.cnrs-mrs.fr/elnemo/start.html
  • Input
    • PDB structure (file or pasted sequence)
    • Number of modes to calculate (the trivial modes 1 to 6 are excluded in this number)
    • Range of perturbation (DQMIN, DQMAX and DQSTEP)
    • Cutoff used to dentify elastic interactions
  • Output
    • Animation of the vibration
    • Distance fluctations between the C-alpha atoms
    • Mean square displacement of all C-alpha atoms

Anisotropic Network Model

Anisotropic network model (ANM) is an elastic network and was introduced by Doruker et al. and Atilgan et al. in 2000<ref name=anm_doruker>Doruker P, Atilgan AR, Bahar I. "Dynamics of proteins predicted by molecular dynamics simulations and analytical approaches: application to alpha-amylase inhibitor.". Proteins. 2000 Aug 15;40(3):512-24. PubMed</ref><ref name=anm_atilgan>Atilgan AR, Durell SR, Jernigan RL, Demirel MC, Keskin O, Bahar I. "Anisotropy of fluctuation dynamics of proteins with an elastic network model.". Biophys J. 2001 Jan;80(1):505-15. PubMed</ref>. The ANM webserver is used to calculate the global modes and was published by Eyal et al. in 2006<ref name=anm_eyal>Eyal E, Yang LW, Bahar I. "Anisotropic network model: systematic evaluation and a new web interface.". Bioinformatics. 2006 Nov 1;22(21):2619-27. Epub 2006 Aug 23. PubMed</ref>. The nodes are represented by the C-alpha atoms and the calculation is based on the spring force constant γ.

Usage

  • Webserver: http://ignmtest.ccbb.pitt.edu/cgi-bin/anm/anm1.cgi
  • Input
    • PDB structure (ID or file)
    • Cutoff for interactions
    • Distance weight for interactions
  • Output
    • Animation of the vibration
    • Fluctuation profiles
    • Inter-residue distance fluctuations
    • Deformation energies

oGNM

NOMAD-Ref

All-atom NMA using Gromacs on the NOMAD-Ref server

Results & Discussion

α-galactosidase A

We applied five webserver (WEBnm@, ElNemo, ANM, oGNM and NOMAD Ref) with the protein α-galactosidase A. We used the PDB structure 1R47 as input and restricted the computation to chain A. The results for the modes 7 to 12 are discussed in the following sections.

Mode 7

WEBnm@ ElNemo ANM oGNM NOMAD Ref
GLA webnma nma mode 7.gif
GLA webnma elnemo mode 7 1.gif
GLA webnma anm mode 7.gif
GLA webnma ognm mode 7.png
GLA webnma nomad ref mode 7.gif


The results of the webservers display a high degree of consensus. Especially the beta strands on the right side show a strong movement in the down direction. This movement can also be expressed as a twist movement which is the strongest in the results of ANM, ElNemo and NOMAD Ref. It seems like that there is also a fluctuation of a large part on the left side which could be driven by the twist movement on the right side.

These obversations in the animations can also be seen in the presentation of oGNM that highlights the beta strands on the right side and also indicates some sort of flexibility on the left side. It is very likely that these flexibility results in a motion which is similar to the movement of the results of the other webservers.

Mode 8

WEBnm@ ElNemo ANM oGNM NOMAD Ref
GLA webnma nma mode 8.gif
GLA webnma elnemo mode 8 1.gif
GLA webnma anm mode 8.gif
GLA webnma ognm mode 8.png
GLA webnma nomad ref mode 8.gif


Similar to mode 7, there is a strong type of twist movement in the outer parts of the molecule for the most part of the results. This time the fluctuation goes upwards on the right side in the results of ANM, ElNemo and WEBnm@. These three results differ only in the fluctuations of the left side. The range of motion decreases from high (ANMN) to almost none (WEBnm@).

The mode 8 of NOMAD Ref is very similar to the mode 7. Mode 8 shows an increase of flexibility which is especially the case for the left part. The result of NOMAD Ref is similar to the three results of ANM, ElNemo and WEBnm@, but the right side moves slightly in a different direction.

The presentation of oGNM differs totally from the other results. oGNM predicts that the most flexible regions are alpha-helices on the upper and lower part. These areas barely move in the results of the other webservers.


Mode 9

WEBnm@ ElNemo ANM oGNM NOMAD Ref
GLA webnma nma mode 9.gif
GLA webnma elnemo mode 9 1.gif
GLA webnma anm mode 9.gif
GLA webnma ognm mode 9.png
GLA webnma nomad ref mode 9.gif


Again, we have the most kind of movement in the outer parts of the molecule and the five results exhibit a high degree of consensus. This time, the fluctuation resembles more the movement of a hinge than a twist movement. This is especially observable for the results of ANM and NOMAD Ref, because the range of motions seems to be bigger than in the results of WEBnm@ and ElNemo.

The presentation of oGNM also indicates a flexible region of an alpha-helix in the middle, which is also noticeable in the results of ElNemo, ANM and NOMAD Ref.


Mode 10

WEBnm@ ElNemo ANM oGNM NOMAD Ref
GLA webnma nma mode 10.gif
GLA webnma elnemo mode 10 1.gif
GLA webnma anm mode 10.gif
GLA webnma ognm mode 10.png
GLA webnma nomad ref mode 10.gif


The results of ElNemo, ANM and NOMAD Ref display some kind of twist movement in the outer parts of the molecule. This twist movement is achieved by the flexibility of three small areas in the result of ANM and probably also in the result of ElNemo, whereas the result of NOMAD Ref seems to be more global. That is, the range of motion seems to be the biggest for NOMAD Ref and the smallest for ElNemo.

A totally different kind of motion is displayed by the result of WEBnm@. The result shows some sort of a stretch motion in the horizontal axis with the right part being the most flexible one.

oGNM indicate flexible regions in alpha-helices in the upper and lower part of the molecule. This is very similar to mode 8, but this time the flexibility is lower which can be seen by the less intense red. Just like at mode 7, the result shows distinction to the result of the other webservers.


Mode 11

WEBnm@ ElNemo ANM oGNM NOMAD Ref
GLA webnma nma mode 11.gif
GLA webnma elnemo mode 11 1.gif
GLA webnma anm mode 11.gif
GLA webnma ognm mode 11.png
GLA webnma nomad ref mode 11.gif


It is very tough to recognize a clear trend in the movements for this mode. The result of ANM shows the same stretching motion like mode 10 of WEBnm@. It could be the case that the molecule undergoes the same fluctuations in the result of ElNemo, but it is not that noticeable.

The result of WEBnm@ is very similar to the twist motion of mode 10 of ElNemo, ANM and NOMAD Ref.

The presentation of oGNM is tough to interpret. There are two very flexible regions, i.e. a coil on the left side and a coil and a part of a beta strand on the right side. We do not assume that these flexible regions result in a motion witch is similar to a mode of the other webservers.


Mode 12

WEBnm@ ElNemo ANM oGNM NOMAD Ref
GLA webnma nma mode 12.gif
GLA webnma elnemo mode 12 1.gif
GLA webnma anm mode 12.gif
GLA webnma ognm mode 12.png
GLA webnma nomad ref mode 12.gif

A distinct coil region on the left side shows a very high amount of motion in the result of ANM. This movement can also be observed in the result of ElNemo. The presentation of oGNM is almost entirely blue which suggest an rigid molecule. Only the coil on the left side and a coil region on the right side have a red color and this fits very well with the observed motions in the result of ANM and ElNemo.

The result of NOMAD Ref displays the twist movement which was also observed in ealier modes and it covers a very wide range of motion.

The molecule in the result of WEBnm@ undergoes a stretching motion which is very similar to the movement in the mode 10 of WEBnm@.

BPTI

Mode NOMAD All Atom (600K) NOMAD All Atom (2000K) NOMAD Ref
Mode 7
BPTI nomad atom 600K 7.gif
BPTI nomad atom 2000K 7.gif
BPTI nomad ref 7.gif
Mode 7
BPTI nomad atom 600K 8.gif
BPTI nomad atom 2000K 8.gif
BPTI nomad ref 8.gif
Mode 7
BPTI nomad atom 600K 9.gif
BPTI nomad atom 2000K 9.gif
BPTI nomad ref 9.gif
Mode 7
BPTI nomad atom 600K 10.gif
BPTI nomad atom 2000K 10.gif
BPTI nomad ref 10.gif
Mode 7
BPTI nomad atom 600K 11.gif
BPTI nomad atom 2000K 11.gif
BPTI nomad ref 11.gif
Mode 7
BPTI nomad atom 600K 12.gif
BPTI nomad atom 2000K 12.gif
BPTI nomad ref 12.gif

Discussion

First of all, we discuss a general observation about the results. It was noticeable that the higer the examined mode is the more divergent the results get. But it was still possible to see similarities between the results of the different webservers. Not only in the same mode, but also in nearby modes. This is not surprising, since the higher the mode gets the more degrees of freedom are available and hence it gets more complex.

The most abundant type of motion was a twist movement of the molecule. This is, the left and the right parts move in opposite directions. It is tough to connect this sort of movement with the biological meaning of the protein. It is somehow noteworthy that the region where the ligand of the proteins binds (referring to the animations, this region is nearby some alpha helices on the left side, see figure X) was rarely involved in the movements.

References

<references/>