Normal Mode Analysis of Glucocerebrosidase
- 1 Introduction
- 2 Elastic and Gaussian Network Models
- 3 Comparison of All-atom NMA to Elastic Network NMA
- 4 Comparison NMA / MD
- 5 References
Nomal mode analysis (NMA) is a time-independent method to calculate vibrational modes and protein flexibility and therefore helps to describe slow large-amplitude/ low-frequency motions. These motions describe conformational changes which are essential for the protein function. Each NMA technique uses harmonic approximations and neglect solvent damping. For the analysis atoms are modeled as point masses connected by springs, which represent the interatomic forces. The normal mode describes a pattern of motion in which all parts of the system (all atoms or only Cα atoms) move with the same frequency and in phase. These modes are no actual protein motions, but can give insight in the possible movements and regions of flexibility. <ref> Alexandrov et al., Normal modes for predicting protein motions: A comprehensive database assessment and associated Web tool. 2004. Protein Science, Vol. 14 </ref> Normal modes can either be calculated with all atoms or can be restrained to Cα atoms only.
In this section different normal mode analyses were carried out for glucocerebrosidase (2NT0). Furthermore, the difference between all-atom and elastic network normal mode analysis gets illustrated on a small protein (1BPT).
Elastic and Gaussian Network Models
The Elastic Network Model is a particular type of NMA in which the springs connecting each node to the neighbouring nodes are of equal strength and only the atom pairs within a cutoff distance get considered. In the Gaussian Network Model the proteins are represented by point masses of the Cα atoms and a harmonic approximation is used to model the interactions between point masses.
WEBnm@ is a web-server which provides automated computation and analysis of low-frequency normal modes of proteins. For an input structure file in PDB format, the normal modes are calculated and a series of automated calculations (normalized squared atomic displacements, vector field representation and an animation of the first six vibrational modes). Additional to the single analysis (calculating lowest frequency normal modes for one protein), the server offers a comperative analysis. In this analysis, the normal modes of a set of aligned sequences are calculated and compared. This feature is still under development. <ref>Hollup SM, Sælensminde G, Reuter N. WEBnm@: a web application for normal mode analysis of proteins BMC Bioinformatics. 2005 Mar 11;6(1):52 </ref>
- Webserver: http://apps.cbu.uib.no/webnma/home
- Single Analysis: structure file in PDB format
- Comparative Analysis: structure files in PDB format and sequence files in fasta format
Figure 1 to 5 show the first five vibrational modes (modes 7 to 11) calculated with WEBnm@ for the structure of glucocerebrosidase (2NTO). On the left side of each figure the vectors are shown, whereas the movement is illustrated to the right.
The atomic displacement of the different models is shown in the following pdf file: File:Webnma atomic displacement glucocerebrosidase.pdf. The plots show the normalized square of the displacement of each Calpha atom, whereof the highest values correspond to the most displaced regions.
Having a look at the atomic displacements, one can see that modes 7, 9, 10 and 11 have peaks and high values all over the protein. This indicates, that most parts of glucocerebrosidase are flexible. Mode 8 describes instead local flexibilities at the N- and C- terminus of the protein. The movement consists movements and displacements of the different domains. These movements consist of streching, bending, twisting, attraction or repulsion of the different glucocerebrosidase domains.
ElNémo is a web-server based on the Elastic Network Model. The tool computes, visualizes and analyses low-frequency normal modes of large macro-molecules. Any size of proteins can be treated as ElNémo uses a 'rotation-translation-block' approximation. <ref>K. Suhre & Y.H. Sanejouand, ElNemo: a normal mode web-server for protein movement analysis and the generation of templates for molecular replacement. Nucleic Acids Research, 32, W610-W614, 2004. </ref>
- Webserver: http://www.igs.cnrs-mrs.fr/elnemo/start.html
- structure file in PDB format
- supplementary options for NMA calculation (number of normal modes to be calculated, minimum and maximum perturbation and stepsize between DQMIN and DQMAX)
- options for computing the eigenmodes (NRBL and cutoff to identify elastic interactions)
Figure 6 to 9 show for each of the five first vibrational normal modes calculated with ElNémo three different perspectives.
The normal modes obtained with ElNémo seem to describe the same movements as the resulting normal modes of Webnm@. Most of the movement seems to be between the different domains. As observed in the results of Webnm@ the movements consists of stretching, bending, attraction, repulsion and twisting.
Anisotropic Network Model
The ANM Webserver provides NMA Analysis with the anisotropic network model (ANM) which is an elastic network (EN) introduced in 2000. It is a very fast method which predicts the global modes. The ANM adopts a uniform force constant γ to all springs and nodes are refered as the Cα atoms.<ref>Anisotropic network model: systematic evaluation and a new web interface, Eyal E, Yang LW, Bahar I. Bioinformatics. 22, 2619-2627, (2006)</ref>
- Webserver: http://ignmtest.ccbb.pitt.edu/cgi-bin/anm/anm1.cgi
- structure file in PDB format or PDB id and chain
- cutoff for interaction between Cα atoms in Å
- distance weight for interaction between Cα atoms
After submitting the job you get a results page where you can analyze the different normal modes more precisely. You can change the amplitude and the frequency of the motions, you can visualize the vectors and change the appearance of the protein. Furthermore you have several possibilities to choose:
- download files
- create PDB (motion)
- create PyMol script
- get anisotrpic temp. factors
- B-factors/mode fluctuations
- Distance fluctuations and deformation energy
The movement, distance fluctuations, the deformation energy per position and mode fluctuations are shown for five different modes in the following figures. Additional modes can be found here.
At the end you can see the Cα anisotropic temperature factors as ellipsoids calculated in different ways.
|Mode 1||Mode 2||Mode 3||Mode 4||Mode 5|
|Inter Distance Analysis and Deformation Analysis|
In each of the five best normal modes obtained, the whole protein seems to be in motion. This is confirmed by the various peaks in the mode fluctuations and the deformation energy plots. Various movements can be observed: twisting, turning, repulsion, attraction and bending.
- Webserver: http://ignm.ccbb.pitt.edu/Online_GNM.htm
- structure file in PDB format or PDB id
- number of nodes to represent a nucleotide
- cutoff distances for amino acids and nucleotide pairs
- preferred visualization engine
also available here: http://ignm.ccbb.pitt.edu/ognm/795639148/temp/index.htm
|Mode 1||Mode 2||Mode 3||Mode 4||Mode 5|
|Mobility-Color-Coded Ribbon Diagram for modes 1 to 5. Red means high mobility, blue low mobility.|
|Fluctuations per residue for modes 1 to 5|
|Cross-correlation plot for mode 1 to 5|
NOMAD-Ref is a webserver for normal mode analysis. By simplifying the potential energy between atoms all atoms can be used. This is possible by using the Tirion calculation of normal modes.<ref>Tirion, Monique M. Large Amplitude Elastic Motions in Proteins from a Single-Parameter, Atomic Analysis. Phys Rev Lett, 77, 1905--1908 (1996).</ref> NOMAD-Ref provides pdb-files for pymol and the elastic network of different normal modes, a graph with the RMSD values per residue and a modes.dat file with gives information about the frequency of the modes and the eigenvectors.
- Webserver: http://lorentz.dynstr.pasteur.fr/nma/submission.php
- PDB file
- number of modes to calculate
- distance weight parameter in Å
- cutoff to use for mode calculation in Å
- average RMSD (Angstrom) in output trajectories
Modes 1-6 are rigid body motions, so the interesting ones are mode 7-11. For each mode we have two figures:
- a) shows the motion of the protein - for modes 1-6 there are no real motions because these are rigid
- c) shows the rmsd per residue
|Mode 7||Mode 8||Mode 9||Mode 10||Mode 11|
Mode 7 to 9 show a movement all over the protein (see plots for rmsd per residue), whereas Mote 10 and 11 have a high flexibility at the N-terminus of the protein.
Comparison of All-atom NMA to Elastic Network NMA
All-atom NMA considers each atom for the calculation of the modes, not only the Cα atoms and is therefore quite time conusming. In this section, all-atom normal mode analyses are compared to the normal modes based on elastic networks. Instead of glucocerebrosidase (our main protein of interest) 1BPT was chosen, as it is quite small and therefore enables a quite fast all-atom normal mode analysis.
The all-atom normal mode analysis was carried out two times with two different temperatures (600 and 2000 K). The temperature was modified as low frequencies dynamics of proteins are temperature dependent and it would therefore be interesting to see the differences.
1BPT at 600 K
1BPT at 2000 K
Elastic Network NMA
The following images show the normal modes 7 to 10 for 1BPT created with Nomad-Ref.
Comparison NMA / MD