Difference between revisions of "Normal mode analysis Gaucher Disease"
(→Introduction) |
(→Deformation Energies) |
||
(47 intermediate revisions by 2 users not shown) | |||
Line 1: | Line 1: | ||
== Introduction == |
== Introduction == |
||
+ | The protein function often depends on its conformation and dynamical properties. To get a complete picture of the dynamic properties of proteins, the traditional de novo method is known to be '''Molecular dynamics''' (MD) simulation. These simulations consider both harmonic and anharmonic motions, and can provide insights into the dynamic of proteins. By taking into account all atoms of a protein, MD simulations shed also light on small motions. However, fine grained MD simulations are computationally very costly which restricts their use to rather short time frames. |
||
+ | An alternative method is '''Normal mode analysis''' (NMA) which has become a popular and often used theoretical tool in the study of functional motions in enzymes, viruses, and large protein assemblies<ref name="NMA1">Eric C Dykeman and Otto F Sankey.(2010). [http://iopscience.iop.org/0953-8984/22/42/423202/ Normal mode analysis and applications in biological physics]. JOURNAL OF PHYSICS.</ref>. NMA is based on a physical theory about the motion of an oscillation system where all parts within the system move sinusoidally with the same frequency and with a fixed phase relation. By using it to study the protein dynamical motion, the atoms are considered as point masses connected by springs to simulate the inter-atomic forces. NMA uses only harmonic approximations and anharmonic motions are neglected. <ref name="NMA2">Adam D. Schuyler, Gregory S. Chirikjian.(2003). [http://custer.lcsr.jhu.edu/wiki/images/7/76/Schuyler03.pdf Normal mode analysis of proteins: a comparison of rigid |
||
− | The protein function often depends on its motion or dynamical properties. To get a complete picture of the dynamic properties of proteins, the traditional de novo mehthod is known to be '''Molecular dynamics''' (MD) simulation. These simulations consider both harmonic and anharmonic motions, and produce great insights into protein dynamics, and when the timescales of the dynamics is short enough, MD is a great tool. However this method is still very computationally consuming. |
||
+ | cluster modes with Cα coarse graining]. Journal of Molecular Graphics and Modelling.</ref>. Thus, the simulations are not as detailed as MD simulations. Instead, a normal mode refers to the harmonic motion of larger parts of a protein such as domains. Since the normal models can be computed efficiently by matrix decomposition, motions can be simulated over a larger time frame compared to MD simulations. The shape of each normal mode is determined by an eigenvector whose eigenvalue corresponds to the frequency of the normal mode. Normal modes are numbered by their frequency and the lowest frequency mode 7 is often closets to the biological motion of an protein. |
||
+ | In this task, we used different normal mode analysis servers to study the normal modes of chain A of glucocerebrosidase ([http://www.rcsb.org/pdb/explore/explore.do?structureId=2nt0 2nt0]) and the mutation structures of W209R and L470P. 2nt0_A is has N=496 residues such that there are altogether 3N-6=1482 normal modes. Each mode is defined by an eigenvector and its corresponding eigenfrequency. The eigenvector contains the amplitude and direction of motion for each atom. In mode i, all N atoms oscillate at the same frequency. |
||
− | An alternative method is '''Normal mode analysis''' which has become a popular and often used theoretical tool in the study of functional motions in enzymes, viruses, and large protein assemblies<ref name="NMA">Eric C Dykeman and Otto F Sankey.(2010). [http://iopscience.iop.org/0953-8984/22/42/423202/ Normal mode analysis and applications in biological physics]. JOURNAL OF PHYSICS.</ref>. |
||
Technical details are reported in our [[Gaucher_Task09_Protocol|protocol]]. |
Technical details are reported in our [[Gaucher_Task09_Protocol|protocol]]. |
||
== WEBnm@ == |
== WEBnm@ == |
||
+ | |||
+ | === Introduction === |
||
+ | |||
+ | [http://apps.cbu.uib.no/webnma/home WEBnm@]<ref name="webnma">Siv Hollup, Gisle Salensminde, and Nathalie Reuter. WEBnm@: a web application for normal mode analyses of proteins. BMC Bioinformatics, 6(1):52+, 2005.</ref> is a web server for automatically computing and analysing low-frequency normal modes of proteins. WEBnm@ is based on the Elastic Network Model (ENM) and only takes into account the C-alpha atoms. For computing the 200 lowest frequency mode, WEBnm@ employs the MMTK package[http://dirac.cnrs-orleans.fr/MMTK/]. The server offers various methods for investigating the resulting normal modes: |
||
+ | |||
+ | *Visualization of the modes |
||
+ | *Deformation Energies |
||
+ | *Atomic displacements and normalized squared fluctuations |
||
+ | *Correlation matrix for investigating the which residue movements are correlated |
||
+ | |||
+ | <br/> |
||
+ | |||
=== Deformation Energies === |
=== Deformation Energies === |
||
+ | The deformation energy is the potential energy of the motion described by a normal mode. <xr id="tab:defor_enery"/> to <xr id="tab:defor_enery_L470P"/> list the deformation energy of mode 7 to 20 for the normal structure and the two mutated ones and <xr id="fig:Eigenvalues_plot"/> to <xr id="fig:Eigenvalues_plot_L470P"/> show the respective eigenvalue which corresponds to the frequency of the mode. The deformation energy is positive correlated with the eigenvalue, i.e. modes with a higher frequency also exhibit a higher deformation energy. This can be motivated by the fact that low modes usually characterize movements of large parts of the protein whereas several smaller parts are moving in opposite directions in case of higher modes which results in a higher deformation energy. |
||
− | "Deformation energies and eigenvalues reflect the energy associated with each mode and are inversely related to the amplitude of the motion described by a the corresponding modes. The protein structure moves along all the normal modes at once. A mode should only be interpreted in isolation if it is energetically well separated from other modes. " - from Webnma@ |
||
+ | <table> |
||
− | Below are the values of the deformation energy for modes 7 to 20. |
||
+ | <tr> |
||
+ | <td> |
||
<figtable id="tab:defor_enery"> |
<figtable id="tab:defor_enery"> |
||
{| style="border-collapse: separate; border-style: solid; border-spacing: 0; border-width: 2px 0 2px 0; text-align:center" width="400px" |
{| style="border-collapse: separate; border-style: solid; border-spacing: 0; border-width: 2px 0 2px 0; text-align:center" width="400px" |
||
|- style="background-color: lightgrey" |
|- style="background-color: lightgrey" |
||
− | ! Mode Index !! Deformation Energy |
+ | ! Mode Index !! Deformation Energy |
|- |
|- |
||
| colspan=6 style="border-collapse: separate; border-style: solid; border-spacing: 0; border-width: 1px 0 0 0"| |
| colspan=6 style="border-collapse: separate; border-style: solid; border-spacing: 0; border-width: 1px 0 0 0"| |
||
|- |
|- |
||
− | |7 || 1663.91 |
+ | |7 || 1663.91 |
|- |
|- |
||
− | | 8 || 2377.48 |
+ | | 8 || 2377.48 |
|- |
|- |
||
− | | 9 || 2720.71 |
+ | | 9 || 2720.71 |
|- |
|- |
||
− | | 10 || 5191.86 |
+ | | 10 || 5191.86 |
|- |
|- |
||
− | | 11 || 5033.67 |
+ | | 11 || 5033.67 |
|- |
|- |
||
− | | 12 || 6174.70 |
+ | | 12 || 6174.70 |
|- |
|- |
||
− | | 13 || 6360.72 |
+ | | 13 || 6360.72 |
+ | |- |
||
+ | | 14 || 6698.31 |
||
+ | |- |
||
+ | | 15 || 9791.69 |
||
+ | |- |
||
+ | | 16 || 9534.10 |
||
+ | |- |
||
+ | | 17 || 10022.74 |
||
+ | |- |
||
+ | | 18 || 11137.44 |
||
+ | |- |
||
+ | | 19 || 11592.80 |
||
+ | |- |
||
+ | | 20 || 12045.03 |
||
|} |
|} |
||
− | <caption>The deformation energy |
+ | <caption>The deformation energy of mode 7 to 20 <br/>of protein structure 2NT0 chain A computed by Webnm@.</caption> |
</figtable> |
</figtable> |
||
+ | </td><td> |
||
− | <br style="clear:both;"> |
||
− | < |
+ | <figtable id="tab:defor_enery_W209R"> |
+ | {| style="border-collapse: separate; border-style: solid; border-spacing: 0; border-width: 2px 0 2px 0; text-align:center" width="400px" |
||
− | [[File:Eingenvalues_plot_2NT0.png|thumb|400px|left|<caption> The Eigenvalues Plot for for protein structure 2NT0 chain A from Webnma@.</caption>]] |
||
+ | |- style="background-color: lightgrey" |
||
+ | ! Mode Index !! Deformation Energy |
||
+ | |- |
||
+ | | colspan=6 style="border-collapse: separate; border-style: solid; border-spacing: 0; border-width: 1px 0 0 0"| |
||
+ | |- |
||
+ | |7 || 1663.94 |
||
+ | |- |
||
+ | | 8 || 2377.58 |
||
+ | |- |
||
+ | | 9 || 2720.58 |
||
+ | |- |
||
+ | | 10 || 5188.56 |
||
+ | |- |
||
+ | | 11 || 5035.65 |
||
+ | |- |
||
+ | | 12 || 6174.92 |
||
+ | |- |
||
+ | | 13 || 6353.38 |
||
+ | |- |
||
+ | | 14 || 6706.38 |
||
+ | |- |
||
+ | | 15 || 9791.22 |
||
+ | |- |
||
+ | | 16 || 9533.93 |
||
+ | |- |
||
+ | | 17 || 10019.91 |
||
+ | |- |
||
+ | | 18 || 11138.34 |
||
+ | |- |
||
+ | | 19 || 11592.61 |
||
+ | |- |
||
+ | | 20 || 12043.20 |
||
+ | |} |
||
+ | <caption>The deformation energy of mode 7 to 20 <br/>of 2NT0 chain A with mutation W209R computed by Webnm@.</caption> |
||
+ | </figtable> |
||
+ | |||
+ | </td><td> |
||
+ | |||
+ | <figtable id="tab:defor_enery_L470P"> |
||
+ | {| style="border-collapse: separate; border-style: solid; border-spacing: 0; border-width: 2px 0 2px 0; text-align:center" width="400px" |
||
+ | |- style="background-color: lightgrey" |
||
+ | ! Mode Index !! Deformation Energy |
||
+ | |- |
||
+ | | colspan=6 style="border-collapse: separate; border-style: solid; border-spacing: 0; border-width: 1px 0 0 0"| |
||
+ | |- |
||
+ | | 7 || 1663.66 |
||
+ | |- |
||
+ | | 8 || 2376.36 |
||
+ | |- |
||
+ | | 9 || 2721.02 |
||
+ | |- |
||
+ | | 10 || 5185.20 |
||
+ | |- |
||
+ | | 11 || 5039.42 |
||
+ | |- |
||
+ | | 12 || 6175.61 |
||
+ | |- |
||
+ | | 13 || 6357.13 |
||
+ | |- |
||
+ | | 14 || 6702.52 |
||
+ | |- |
||
+ | | 15 || 9791.92 |
||
+ | |- |
||
+ | | 16 || 9532.10 |
||
+ | |- |
||
+ | | 17 || 10022.80 |
||
+ | |- |
||
+ | | 18 || 11124.30 |
||
+ | |- |
||
+ | | 19 || 11589.96 |
||
+ | |- |
||
+ | | 20 || 12043.17 |
||
+ | |} |
||
+ | <caption>The deformation energy of mode 7 to 20 <br/>of 2NT0 chain A with mutation L470P computed by Webnm@.</caption> |
||
+ | </figtable> |
||
+ | |||
+ | </td><td> |
||
+ | |||
+ | <tr> |
||
+ | <td> |
||
+ | |||
+ | <figure id="fig:Eigenvalues_plot"> |
||
+ | [[File:Eingenvalues_plot_2NT0.png|thumb|350px|right|<caption>Normal mode eigenvalues of 2NT0 chain A computed by Webnma@.</caption>]] |
||
</figure> |
</figure> |
||
+ | |||
+ | <td> |
||
+ | |||
+ | <figure id="fig:Eigenvalues_plot_W209R"> |
||
+ | [[File:Eingenvalues_plot_2NT0_W209R.png|thumb|350px|right|<caption>Normal mode eigenvalues of 2NT0 chain A with mutation W209R computed by Webnma@.</caption>]] |
||
+ | </figure> |
||
+ | |||
+ | <td> |
||
+ | |||
+ | <figure id="fig:Eigenvalues_plot_L470P"> |
||
+ | [[File:Eingenvalues_plot_2NT0_L470P.png|thumb|350px|right|<caption>Normal mode eigenvalues of 2NT0 chain A with mutation L470P computed by Webnma@.</caption>]] |
||
+ | </figure> |
||
+ | |||
+ | </td></tr></table> |
||
+ | |||
<br style="clear:both;"> |
<br style="clear:both;"> |
||
− | === |
+ | === Correlation Matrix === |
+ | <xr id="fig:corr_matr"/> depicts correlated motions of C-alpha. Each cell shows the isotropic correlation of two residues in the protein from -1 (anti-correlated) over 0 (uncorrelated) to 1 (correlated). |
||
− | " Atomic displacements: |
||
+ | The red areas along the diagonal suggest smaller local correlations. In addition, the N terminal region (1-100) is correlated with the C terminal region (350-500), that is, the beginning and the end of the protein exhibit similar motions. |
||
+ | <br style="clear:both;"> |
||
− | The square of the displacement of each Calpha atom (for modes 7 to 12), normalized so that the sum over all residues is equal to 100. |
||
+ | <figure id="fig:corr_matr"> |
||
− | Highest values correspond to the most displaced regions. On the plots, one should look for cluster of peaks, those identify significantly big regions. Isolated peaks reflect local flexibility and are not relevant. |
||
+ | [[File:2NT0_correlation_matrix.png|thumb|400px|left|<caption> The correlation matrix for protein structure 2NT0 chain A from Webnma@. Each cell in the plot shows the isotropic correlation of two residues in the protein on a range from -1 (anti-correlated) via 0 (uncorrelated) to 1 (correlated).</caption>]] |
||
− | Fluctuations |
||
+ | </figure> |
||
+ | <br style="clear:both;"> |
||
− | The square of the fluctuation of each Calpha atom (for all non-trivial modes), normalized so that the sum over all residues is equal to 100. |
||
+ | |||
+ | === Atomic Displacement === |
||
+ | <xr id="fig:atom_disp_7to12"/> shows the normalized square of the displacement of each C-alpha atom for mode 7 to 12. Narrow peaks denote local movements whereas broad peaks indicate larger flexible protein regions. <xr id="fig:fluc_plot"/> depicts the normalized fluctuation of each C-alpha atom calculated by averaging over the atomic displacements of all modes which are weighted by their eigenvalue. |
||
− | The fluctuations are the sum of the atomic displacements in each mode weighted by the inverse of their corresponding eigenvalues. These are equivalent to the normalized temperature factors. |
||
+ | Except for mode 9 and mode 10, the most flexible regions are at the beginning and at the end of the protein. All modes have clear peaks around range 55-65 and 439-442, which are also pronounces in <xr id="fig:fluc_plot"/>. Furthermore, these two regions seem to be correlated (<xr id="fig:corr_matr"/>). We therefore visualized residue 55-65 and 439-442 (<xr id="fig:webnma_displacement"/>) and found that they are part of two loop regions on the surface. These two loop regions are significantly more flexible than other parts of the protein such as the active site residues E235 and E340 which both exhibit small displacement values for most modes. |
||
− | Results can be retrieved as plots at the PNG format or raw data (first column: resid, second column: normalized squared atomic displacement/fluctuation) " |
||
− | - from Webnma@ |
||
<br style="clear:both;"> |
<br style="clear:both;"> |
||
Line 71: | Line 202: | ||
<br style="clear:both;"> |
<br style="clear:both;"> |
||
+ | {|align="left" |
||
+ | |<figure id="fig:fluc_plot">[[File:2NT0_fluctuationsplot.png|thumb|left|200px|<caption> The normalized squared fluctuations for all modes for protein structure 2NT0 chain A from Webnma@.</caption>]]</figure> |
||
+ | |<figure id="fig:webnma_displacement">[[File:webnma_displacemnent.png|thumb|left|200px|<caption>Visualization of the two peaked regions of Figure 4. Red: residue 55-65; Orange: residue 439-443.</caption>]]</figure> |
||
+ | |} |
||
+ | <br style="clear:both;"> |
||
+ | === Mode Visualization === |
||
− | <figure id="fig:fluc_plot"> |
||
+ | |||
− | [[File:2NT0_fluctuationsplot.png|thumb|400px|left|<caption> The normalized squared fluctuations for all modes for protein structure 2NT0 chain A from Webnma@.</caption>]] |
||
+ | |||
+ | <figure id="fig:webnm_model7"> |
||
+ | <gallery widths=200px heights=200px perrow=3> |
||
+ | File:mode7_animation.gif|Animation of normal mode 7 |
||
+ | File:mode7_vector.png|Vectors of normal mode 7 |
||
+ | Mode 7 exhibits a bending movement. |
||
+ | </gallery> |
||
+ | <caption>The visualization of mode 7 and the vectors showing the direction of the movements.</caption> |
||
</figure> |
</figure> |
||
+ | <br/> |
||
− | <br style="clear:both;"> |
||
− | === Correlation Matrix === |
||
+ | <figure id="fig:webnm_model8"> |
||
− | "The correlation matrix shows the correlated movement of the Calphas in the protein. Both axis denote the Calphas of the protein in sequential order, so that each cell in the plot shows the isotropic correlation of two residues in the protein on a range from -1 (anti-correlated) via 0 (uncorrelated) to 1 (correlated). " - from Webnma@ |
||
+ | <gallery widths=200px heights=200px perrow=3> |
||
+ | File:mode8_animation.gif|Animation of normal mode 8 |
||
+ | File:mode8_vector.png|Vectors of normal mode 8 |
||
+ | In Mode 8, loop regions are moving most. The amplitude is smaller compared to mode 7. |
||
+ | </gallery> |
||
+ | <caption>The visualization of mode 8 and the vectors showing the direction of the movements.</caption> |
||
+ | </figure> |
||
+ | <br/> |
||
− | <br style="clear:both;"> |
||
− | <figure id="fig: |
+ | <figure id="fig:webnm_model9"> |
+ | <gallery widths=200px heights=200px perrow=3> |
||
− | [[File:2NT0_correlation_matrix.png|thumb|400px|left|<caption> The correlation matrix for protein structure 2NT0 chain A from Webnma@.</caption>]] |
||
+ | File:mode9_animation.gif|Animation of normal mode 9 |
||
+ | File:mode9_vector.png|Vectors of normal mode 9 |
||
+ | Mode 9 exhibits a shearing moving of the two protein domains. |
||
+ | </gallery> |
||
+ | <caption>The visualization of mode 9 and the vectors showing the direction of the movements.</caption> |
||
</figure> |
</figure> |
||
+ | <br/> |
||
− | <br style="clear:both;"> |
||
+ | <figure id="fig:webnm_model10"> |
||
− | === Mode Visualization === |
||
+ | <gallery widths=200px heights=200px perrow=3> |
||
+ | File:mode10_animation.gif|Animation of normal mode 10 |
||
+ | File:mode10_vector.png|Vectors of normal mode 10 |
||
+ | For mode 10, also alpha helical regions in the middle of the protein are moving with a small amplitude. |
||
+ | </gallery> |
||
+ | <caption>The visualization of mode 10 and the vectors showing the direction of the movements.</caption> |
||
+ | </figure> |
||
+ | |||
+ | <br/> |
||
+ | |||
+ | <figure id="fig:webnm_model11"> |
||
+ | <gallery widths=200px heights=200px perrow=3> |
||
+ | File:mode11_animation.gif|Animation of normal mode 11 |
||
+ | File:mode11_vector.png|Vectors of normal mode 11 |
||
+ | In mode 11, all parts of the protein a slightly moving similar to mode 10. |
||
+ | </gallery> |
||
+ | <caption>The visualization of mode 11 and the vectors showing the direction of the movements.</caption> |
||
+ | </figure> |
||
+ | |||
+ | <br/> |
||
+ | |||
+ | <figure id="fig:webnm_model12"> |
||
+ | <gallery widths=200px heights=200px perrow=3> |
||
+ | File:mode12_animation.gif|Animation of normal mode 12 |
||
+ | File:mode12_vector.png|Vectors of normal mode 12 |
||
+ | Mode 12 exhibits a hinge like movements which displaces in particular loop region 55-65. |
||
+ | </gallery> |
||
+ | <caption>The visualization of mode 11 and the vectors showing the direction of the movements.</caption> |
||
+ | </figure> |
||
+ | |||
+ | <br/> |
||
+ | |||
+ | === Comparison of different conformations === |
||
+ | |||
+ | To do: compare different conformations, for example: [http://www.rcsb.org/pdb/explore/explore.do?structureId=2nt0 2nt0] and [http://www.rcsb.org/pdb/explore/explore.do?structureId=3gxi 3gxi]. |
||
== ElNemo == |
== ElNemo == |
||
+ | |||
+ | |||
+ | === Introduction === |
||
+ | |||
+ | [http://www.igs.cnrs-mrs.fr/elnemo/start.html elNémo]<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 Research, 32(suppl 2):W610–W614, July 2004.</ref> is a normal mode web server for protein movement analysis. Just as Webnm@, ElNeno employs the Elastic Network Model (ENM) for computing the lowest frequency modes. Interaction forces are computes between C-alpha atoms only and a cut-off distance (default 8 A) specifies the number of neighbouring atoms which are taken into account in the computations. The cut-off distance can be arbitrarily adjusted by the user and is a trade-off between runtime and approximation accuracy. ElNema further employs the Rotational Translation Block (RTB) algorithm for speeding up the calculations by grouping a several residues to super residues. The ElNemo web-interface takes a PDB structure as input and computes the 100 lowest frequency modes which can be investigated in different ways: |
||
+ | *viewing 3-D animations of the protein movement for each mode |
||
+ | *comparing the collectivity of the modes |
||
+ | *identifying residues that have the largest distance fluctuations in a given mode. |
||
+ | *comparing B-factors derived from the normal mode decomposition to experimentally determined B-factors for analysing to which extend their flexibility differs |
||
+ | |||
+ | If two conformations of the same protein are uploaded, the user will see the contribution of each mode to the conformational changes (overlap between a protein motion and a normal mode). If two homologous proteins are uploaded, the root mean square distance (RMSD) between all residues as well as the number of residues that are closer than 3A to the function of mode and perturbation can be computed. |
||
+ | |||
+ | === Mode Visualization and distance fluctuation === |
||
+ | |||
+ | In the following figures, the visualization of model 7 to 11 and their distance fluctuation are shown. The structure is presented from three different orthologuous viewpoints where the secondary structures are determined from the C-alpha atom positions of the uploaded protein (N-terminal blue, C-terminal red). In the distance fluctuation map, each cell measures the relative movement between a pair of residues. One can find the rigid and the flexible regions of the protein, as well as their relative movements. The regions where amino acid residues move as rigid bodies are shown in black whereas flexible segments are colored blue or red. The blue parts indicate where the distance between two C-alpha atoms increases significantly, and the red parts where the distance decreases. |
||
<figure id="fig:elne_model7"> |
<figure id="fig:elne_model7"> |
||
− | <gallery |
+ | <gallery widths=200px heights=200px perrow=5> |
− | File:2nt0.7.x.0.gif| |
+ | File:2nt0.7.x.0.gif|Normal Mode 7: view through Z axis |
− | File:2nt0.7.x.90.gif| |
+ | File:2nt0.7.x.90.gif|Normal Mode 7: view through X axis |
− | File:2nt0.7.y.90.gif| |
+ | File:2nt0.7.y.90.gif|Normal Mode 7: view through Y axis |
+ | File:Example.distance_fluctuations_7.png|The distance fluctuations between all C-alpha atoms. |
||
+ | Mode 7 presents a bending and stretching movement. |
||
</gallery> |
</gallery> |
||
− | <caption>The |
+ | <caption>The visualization of mode 7 and its distance fluctuation.</caption> |
</figure> |
</figure> |
||
+ | |||
+ | <br/> |
||
<figure id="fig:elne_model8"> |
<figure id="fig:elne_model8"> |
||
− | <gallery |
+ | <gallery widths=200px heights=200px perrow=5> |
− | File:2nt0.8.x.0.gif| |
+ | File:2nt0.8.x.0.gif|Normal Mode 8: view through Z axis |
− | File:2nt0.8.x.90.gif| |
+ | File:2nt0.8.x.90.gif|Normal Mode 8: view through X axis |
− | File:2nt0.8.y.90.gif| |
+ | File:2nt0.8.y.90.gif|Normal Mode 8: view through Y axis |
+ | File:Example.distance_fluctuations_8.png|The distance fluctuations between all C-alpha atoms. |
||
+ | Mode 8 presents a bending and stretching movement as similar to mode 7 but towards different direction. |
||
</gallery> |
</gallery> |
||
− | <caption>The |
+ | <caption>The visualization of mode 8 and its distance fluctuation.</caption> |
</figure> |
</figure> |
||
+ | |||
+ | <br/> |
||
<figure id="fig:elne_model9"> |
<figure id="fig:elne_model9"> |
||
− | <gallery |
+ | <gallery widths=200px heights=200px perrow=5> |
− | File:2nt0.9.x.0.gif| |
+ | File:2nt0.9.x.0.gif|Normal Mode 9: view through Z axis |
− | File:2nt0.9.x.90.gif| |
+ | File:2nt0.9.x.90.gif|Normal Mode 9: view through X axis |
− | File:2nt0.9.y.90.gif| |
+ | File:2nt0.9.y.90.gif|Normal Mode 9: view through Y axis |
+ | File:Example.distance_fluctuations_9.png|The distance fluctuations between all C-alpha atoms. |
||
+ | Mode 9 presents a twisting movement. |
||
</gallery> |
</gallery> |
||
− | <caption>The |
+ | <caption>The visualization of mode 9 and its distance fluctuation.</caption> |
</figure> |
</figure> |
||
+ | |||
+ | <br/> |
||
<figure id="fig:elne_model10"> |
<figure id="fig:elne_model10"> |
||
− | <gallery |
+ | <gallery widths=200px heights=200px perrow=5> |
− | File:2nt0.10.x.0.gif| |
+ | File:2nt0.10.x.0.gif|Normal Mode 10: view through Z axis |
− | File:2nt0.10.x.90.gif| |
+ | File:2nt0.10.x.90.gif|Normal Mode 10: view through X axis |
− | File:2nt0.10.y.90.gif| |
+ | File:2nt0.10.y.90.gif|Normal Mode 10: view through Y axis |
+ | File:Example.distance_fluctuations_10.png|The distance fluctuations between all C-alpha atoms. |
||
+ | Mode 10 presents a movement of a mixture of bending,stretching and twisting. |
||
</gallery> |
</gallery> |
||
− | <caption>The |
+ | <caption>The visualization of mode 10 and its distance fluctuation.</caption> |
</figure> |
</figure> |
||
+ | |||
+ | <br/> |
||
<figure id="fig:elne_model11"> |
<figure id="fig:elne_model11"> |
||
− | <gallery |
+ | <gallery widths=200px heights=200px perrow=5> |
− | File:2nt0.11.x.0.gif| |
+ | File:2nt0.11.x.0.gif|Normal Mode 11 |
− | File:2nt0.11.x.90.gif| |
+ | File:2nt0.11.x.90.gif|Normal Mode 11: view through X axis |
− | File:2nt0.11.y.90.gif| |
+ | File:2nt0.11.y.90.gif|Normal Mode 11: view through Y axis |
+ | File:Example.distance_fluctuations_11.png|The distance fluctuations between all C-alpha atoms. |
||
+ | Mode 11 presents a twisting movement as similar to mode 9. |
||
</gallery> |
</gallery> |
||
− | <caption>The |
+ | <caption>The visualization of mode 11 and its distance fluctuation.</caption> |
</figure> |
</figure> |
||
+ | === Comparison of different conformations === |
||
− | == Anisotropic Network Model == |
||
+ | To do: compare different conformations, for example: [http://www.rcsb.org/pdb/explore/explore.do?structureId=2nt0 2nt0] and [http://www.rcsb.org/pdb/explore/explore.do?structureId=3gxi 3gxi]. |
||
− | == Discussion == |
||
== References == |
== References == |
Latest revision as of 23:02, 2 September 2012
Contents
Introduction
The protein function often depends on its conformation and dynamical properties. To get a complete picture of the dynamic properties of proteins, the traditional de novo method is known to be Molecular dynamics (MD) simulation. These simulations consider both harmonic and anharmonic motions, and can provide insights into the dynamic of proteins. By taking into account all atoms of a protein, MD simulations shed also light on small motions. However, fine grained MD simulations are computationally very costly which restricts their use to rather short time frames.
An alternative method is Normal mode analysis (NMA) which has become a popular and often used theoretical tool in the study of functional motions in enzymes, viruses, and large protein assemblies<ref name="NMA1">Eric C Dykeman and Otto F Sankey.(2010). Normal mode analysis and applications in biological physics. JOURNAL OF PHYSICS.</ref>. NMA is based on a physical theory about the motion of an oscillation system where all parts within the system move sinusoidally with the same frequency and with a fixed phase relation. By using it to study the protein dynamical motion, the atoms are considered as point masses connected by springs to simulate the inter-atomic forces. NMA uses only harmonic approximations and anharmonic motions are neglected. <ref name="NMA2">Adam D. Schuyler, Gregory S. Chirikjian.(2003). [http://custer.lcsr.jhu.edu/wiki/images/7/76/Schuyler03.pdf Normal mode analysis of proteins: a comparison of rigid cluster modes with Cα coarse graining]. Journal of Molecular Graphics and Modelling.</ref>. Thus, the simulations are not as detailed as MD simulations. Instead, a normal mode refers to the harmonic motion of larger parts of a protein such as domains. Since the normal models can be computed efficiently by matrix decomposition, motions can be simulated over a larger time frame compared to MD simulations. The shape of each normal mode is determined by an eigenvector whose eigenvalue corresponds to the frequency of the normal mode. Normal modes are numbered by their frequency and the lowest frequency mode 7 is often closets to the biological motion of an protein.
In this task, we used different normal mode analysis servers to study the normal modes of chain A of glucocerebrosidase (2nt0) and the mutation structures of W209R and L470P. 2nt0_A is has N=496 residues such that there are altogether 3N-6=1482 normal modes. Each mode is defined by an eigenvector and its corresponding eigenfrequency. The eigenvector contains the amplitude and direction of motion for each atom. In mode i, all N atoms oscillate at the same frequency.
Technical details are reported in our protocol.
WEBnm@
Introduction
WEBnm@<ref name="webnma">Siv Hollup, Gisle Salensminde, and Nathalie Reuter. WEBnm@: a web application for normal mode analyses of proteins. BMC Bioinformatics, 6(1):52+, 2005.</ref> is a web server for automatically computing and analysing low-frequency normal modes of proteins. WEBnm@ is based on the Elastic Network Model (ENM) and only takes into account the C-alpha atoms. For computing the 200 lowest frequency mode, WEBnm@ employs the MMTK package[1]. The server offers various methods for investigating the resulting normal modes:
- Visualization of the modes
- Deformation Energies
- Atomic displacements and normalized squared fluctuations
- Correlation matrix for investigating the which residue movements are correlated
Deformation Energies
The deformation energy is the potential energy of the motion described by a normal mode. <xr id="tab:defor_enery"/> to <xr id="tab:defor_enery_L470P"/> list the deformation energy of mode 7 to 20 for the normal structure and the two mutated ones and <xr id="fig:Eigenvalues_plot"/> to <xr id="fig:Eigenvalues_plot_L470P"/> show the respective eigenvalue which corresponds to the frequency of the mode. The deformation energy is positive correlated with the eigenvalue, i.e. modes with a higher frequency also exhibit a higher deformation energy. This can be motivated by the fact that low modes usually characterize movements of large parts of the protein whereas several smaller parts are moving in opposite directions in case of higher modes which results in a higher deformation energy.
</figtable> </figtable> </figtable></figure>
</figure>
</figure>
<figtable id="tab:defor_enery">
| ||||||||||||||||||||||||||||||||||||
<figtable id="tab:defor_enery_W209R">
| ||||||||||||||||||||||||||||||||||||
<figtable id="tab:defor_enery_L470P">
| ||||||||||||||||||||||||||||||||||||
<figure id="fig:Eigenvalues_plot"> | ||||||||||||||||||||||||||||||||||||
<figure id="fig:Eigenvalues_plot_W209R"> | ||||||||||||||||||||||||||||||||||||
<figure id="fig:Eigenvalues_plot_L470P"> |
Correlation Matrix
<xr id="fig:corr_matr"/> depicts correlated motions of C-alpha. Each cell shows the isotropic correlation of two residues in the protein from -1 (anti-correlated) over 0 (uncorrelated) to 1 (correlated). The red areas along the diagonal suggest smaller local correlations. In addition, the N terminal region (1-100) is correlated with the C terminal region (350-500), that is, the beginning and the end of the protein exhibit similar motions.
<figure id="fig:corr_matr">
</figure>
Atomic Displacement
<xr id="fig:atom_disp_7to12"/> shows the normalized square of the displacement of each C-alpha atom for mode 7 to 12. Narrow peaks denote local movements whereas broad peaks indicate larger flexible protein regions. <xr id="fig:fluc_plot"/> depicts the normalized fluctuation of each C-alpha atom calculated by averaging over the atomic displacements of all modes which are weighted by their eigenvalue.
Except for mode 9 and mode 10, the most flexible regions are at the beginning and at the end of the protein. All modes have clear peaks around range 55-65 and 439-442, which are also pronounces in <xr id="fig:fluc_plot"/>. Furthermore, these two regions seem to be correlated (<xr id="fig:corr_matr"/>). We therefore visualized residue 55-65 and 439-442 (<xr id="fig:webnma_displacement"/>) and found that they are part of two loop regions on the surface. These two loop regions are significantly more flexible than other parts of the protein such as the active site residues E235 and E340 which both exhibit small displacement values for most modes.
<figure id="fig:atom_disp_7to12">
</figure>
<figure id="fig:fluc_plot"> |
<figure id="fig:webnma_displacement"> |
Mode Visualization
<figure id="fig:webnm_model7">
- Mode 7 exhibits a bending movement.
The visualization of mode 7 and the vectors showing the direction of the movements. </figure>
<figure id="fig:webnm_model8">
- In Mode 8, loop regions are moving most. The amplitude is smaller compared to mode 7.
The visualization of mode 8 and the vectors showing the direction of the movements. </figure>
<figure id="fig:webnm_model9">
- Mode 9 exhibits a shearing moving of the two protein domains.
The visualization of mode 9 and the vectors showing the direction of the movements. </figure>
<figure id="fig:webnm_model10">
- For mode 10, also alpha helical regions in the middle of the protein are moving with a small amplitude.
The visualization of mode 10 and the vectors showing the direction of the movements. </figure>
<figure id="fig:webnm_model11">
- In mode 11, all parts of the protein a slightly moving similar to mode 10.
The visualization of mode 11 and the vectors showing the direction of the movements. </figure>
<figure id="fig:webnm_model12">
- Mode 12 exhibits a hinge like movements which displaces in particular loop region 55-65.
The visualization of mode 11 and the vectors showing the direction of the movements. </figure>
Comparison of different conformations
To do: compare different conformations, for example: 2nt0 and 3gxi.
ElNemo
Introduction
elNémo<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 Research, 32(suppl 2):W610–W614, July 2004.</ref> is a normal mode web server for protein movement analysis. Just as Webnm@, ElNeno employs the Elastic Network Model (ENM) for computing the lowest frequency modes. Interaction forces are computes between C-alpha atoms only and a cut-off distance (default 8 A) specifies the number of neighbouring atoms which are taken into account in the computations. The cut-off distance can be arbitrarily adjusted by the user and is a trade-off between runtime and approximation accuracy. ElNema further employs the Rotational Translation Block (RTB) algorithm for speeding up the calculations by grouping a several residues to super residues. The ElNemo web-interface takes a PDB structure as input and computes the 100 lowest frequency modes which can be investigated in different ways:
- viewing 3-D animations of the protein movement for each mode
- comparing the collectivity of the modes
- identifying residues that have the largest distance fluctuations in a given mode.
- comparing B-factors derived from the normal mode decomposition to experimentally determined B-factors for analysing to which extend their flexibility differs
If two conformations of the same protein are uploaded, the user will see the contribution of each mode to the conformational changes (overlap between a protein motion and a normal mode). If two homologous proteins are uploaded, the root mean square distance (RMSD) between all residues as well as the number of residues that are closer than 3A to the function of mode and perturbation can be computed.
Mode Visualization and distance fluctuation
In the following figures, the visualization of model 7 to 11 and their distance fluctuation are shown. The structure is presented from three different orthologuous viewpoints where the secondary structures are determined from the C-alpha atom positions of the uploaded protein (N-terminal blue, C-terminal red). In the distance fluctuation map, each cell measures the relative movement between a pair of residues. One can find the rigid and the flexible regions of the protein, as well as their relative movements. The regions where amino acid residues move as rigid bodies are shown in black whereas flexible segments are colored blue or red. The blue parts indicate where the distance between two C-alpha atoms increases significantly, and the red parts where the distance decreases.
<figure id="fig:elne_model7">
- Mode 7 presents a bending and stretching movement.
The visualization of mode 7 and its distance fluctuation. </figure>
<figure id="fig:elne_model8">
- Mode 8 presents a bending and stretching movement as similar to mode 7 but towards different direction.
The visualization of mode 8 and its distance fluctuation. </figure>
<figure id="fig:elne_model9">
- Mode 9 presents a twisting movement.
The visualization of mode 9 and its distance fluctuation. </figure>
<figure id="fig:elne_model10">
- Mode 10 presents a movement of a mixture of bending,stretching and twisting.
The visualization of mode 10 and its distance fluctuation. </figure>
<figure id="fig:elne_model11">
- Mode 11 presents a twisting movement as similar to mode 9.
The visualization of mode 11 and its distance fluctuation. </figure>
Comparison of different conformations
To do: compare different conformations, for example: 2nt0 and 3gxi.
References
<references/>