Difference between revisions of "Canavan Task 9 - Normal Mode Analysis"

From Bioinformatikpedia
(Webnm@)
(Comparison 2I3C and 2O53)
Line 195: Line 195:
   
 
<figure id="comp_def_ener">[[File:comp_def_energy_2o53_2i3c.png|thumb|250px|<b><xr nolink id="comp_def_ener"/></b><br>]]</figure>
 
<figure id="comp_def_ener">[[File:comp_def_energy_2o53_2i3c.png|thumb|250px|<b><xr nolink id="comp_def_ener"/></b><br>]]</figure>
  +
  +
====Fluctuation====
  +
  +
As you can see in <xr id="comparison_2o53_2i3c_fluctuations"/>, the fluctuations for all modes for both proteins are identical.
  +
  +
<figure id="comparison_2o53_2i3c_fluctuations">[[File:comparison_2o53_2i3c_fluctuations.png|thumb|250px|<b><xr nolink id="comparison_2o53_2i3c_fluctuations"/></b><br>]]</figure>
   
 
==ElNemo==
 
==ElNemo==

Revision as of 16:21, 10 July 2012

Protocol

Further information can be found in the protocol.


Choosing structures

For human Aspartoacylase, there are four structures available in the PDB. For a closer description of the structures have a look at task 4.

One important difference between the structures is the conformation of a loop formed by residues 158 - 164. This loop represents the gate for the binding site, which can either be closed or open. 2O53 and 2O4H both represent the closed conformation, whereas 2I3C and 2Q51 represent the open conformation.

Therefore we thought it would be interesting to investigate the normal modes of a structure with the closed and the open conformation.

We also looked at 2O53 and 2O4H in order to examine if there are differences between the ligand bound and apo-structure.


Webnm@

  • Deformation Energies and Eigenvalue Plot:

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.

  • Atomic Displacement Analysis:

Plots the displacement of each Calpha atom, i.e. highlights which parts of the protein are the most displaced for each mode.

  • Correlation Matrix Analysis:

Plots the correlation of motions between all the Calphas in the protein structure.

  • Mode Visualisation:

Webnm@ provides the possibility to visualize the modes as well as the download option of dcd and vmd files.



2O53

We used the PDB identifier 2O53. Therefore Webnm@ used the dimer from the PDB. 2O53 is a Holodimer without bound substrate and is in the closed conformation.

Deformation Energies

Mode IndexDeformation Energy Mode IndexDeformation Energy
7870.03 143740.78
8976.54 154168.69
91562.41 165263.58
102800.99 175461.10
112798.65 185906.92
122954.75 196126.73
133448.49 206356.04

Eigenvalues Plot

The eigenvalues are inversely related to the motion of the modes. Therefore mode 7 describes the highest amplitude of motion in the protein and the other modes describe a decreasing amplitude of motion.

The eigenvalue plot can be seen in <xr id="2o53_eigenvaluesplot"/>


Atomic Displacements

In the lowest mode (7) there are quite random displacements. They are scattered all over the protein. The same can be said for mode 9. For mode 8, 10 and 12 there are peaks around residue 400, which is a long solvent exposed loop. For modes 8,10,11 there is a peak around residue 100 in chain A which corresponds to the peak for the residues around position 400 in chain B. For mode 11 and 12 there are two sharp peaks around residues 250 and 265. This regions are composed of many solvent exposed flexible loops of the C-terminal domain. [AD for modes 7 to 12]


In the fluctuations plot averaged over all modes, it can be seen that the atomic movements are different for the two monomers. The peaks in this plot are mapped to the dimer in figure .... As one can see, these highly motile regions correlate with exposed loop regions.

Correlation Matrix Analysis

From the correlation matrix it can be found, that the C-terminal domain (residue 211-311) shows correlated movements. In chain B this domain can be found starting at residue 520.

There are some other smaller correalted motions, that are not very emphasized.


Plot

<figtable id="2o53_plots">

<figure id="2o53_eigenvaluesplot">
<xr nolink id="2o53_eigenvaluesplot"/>
</figure>
<figure id="2o53_fluctuationsplot">
<xr nolink id="2o53_fluctuationsplot"/>
</figure>
<figure id="2o53_correlation_matrix">
<xr nolink id="2o53_correlation_matrix"/>
</figure>

</figtable>

2O4H

In 2O4H an intermediate substrate has been cocrystallized and the enzyme is in the closed conformation. We uploaded a modified pdb file, where we changed the HETATM entry into ATOM, so that the ligand will be considered during normal mode analysis.

Deformation Energies

Mode IndexDeformation Energy Mode IndexDeformation Energy
7890.47 143813.56
81001.20 154209.21
91576.14 165299.13
102832.63 175374.76
112837.90 185750.94
122995.56 196278.71
133417.11 206396.79

Eigenvalues Plot

The eigenvalues are inversely related to the motion of the modes. Therefore mode 7 describes the highest amplitude of motion in the protein and the other modes describe a decreasing amplitude of motion.

The eigenvalue plot can be seen in <xr id="2o4h_eigenvaluesplot"/>


Atomic Displacements

The plots for mode 7 to 12 are identical with the plots for 2O53. Therefore the ligand in the structure does not influence the movement of the protein at all. The atomic displacements for all modes can be seen her: [AD for modes 7 to 12]


The fluctuations plot averaged over all modes is also identical with the plot for 2o53. Even very small fluctuations are exactly the same.


Correlation Matrix Analysis

And again, the correlation matrix is also totally identical with the correlation matrix for 2O53.

Plot

<figtable id="2o4h_plots">

<figure id="2o4h_eigenvaluesplot">
<xr nolink id="2o4h_eigenvaluesplot"/>
</figure>
<figure id="2o4h_fluctuationsplot">
<xr nolink id="2o4h_fluctuationsplot"/>
</figure>
<figure id="2o4h_correlation_matrix">
<xr nolink id="2o4h_correlation_matrix"/>
</figure>

</figtable>


2I3C

This structure shows the binding site in an open conformation.

Some amino acids have been replaced by artificial amino acids and are not considered by webnm@ in the normal mode analysis.

Deformation Energies

Mode IndexDeformation Energy Mode IndexDeformation Energy
7819.58 143831.12
8897.71 154275.48
91423.78 165370.37
102661.79 175549.08
112485.43 185226.29
122576.98 196459.76
133362.10 205641.86

Eigenvalues Plot

The eigenvalue plot looks as for the proteins analysed before: almost linear increase in eigenvalue for the modes.

The eigenvalue plot can be seen in <xr id="2i3c_eigenvalues"/>


Plot

<figtable id="2i3c_plots">

<figure id="2i3c_eigenvalues">
<xr nolink id="2i3c_eigenvalues"/>
</figure>
<figure id="2I3C_fluctuationsplot">
<xr nolink id="2i3c_fluctuationsplot"/>
</figure>
<figure id="2o4h_correlation_matrix">
<xr nolink id="2o4h_correlation_matrix"/>
</figure>

</figtable>


Comparison 2I3C and 2O53

Deformation Energy

For mode 7 - 13, 2I3C has the lower deformation energies. Then it is almost vice.versa and 2O53 has lower deformation energies.


<figure id="comp_def_ener">

<xr nolink id="comp_def_ener"/>

</figure>

Fluctuation

As you can see in <xr id="comparison_2o53_2i3c_fluctuations"/>, the fluctuations for all modes for both proteins are identical.

<figure id="comparison_2o53_2i3c_fluctuations">

<xr nolink id="comparison_2o53_2i3c_fluctuations"/>

</figure>

ElNemo

In order for ElNemo to accept our proteins, we had to modify the pdf files. We eliminated all records except for the ATOM records.

2O53

Normal Mode Analysis for ID 12070820013211761

Correlation= 0.322 for 604 C-alpha atoms.


<R2> 	        frequency 	collectivity
mode 7 	1.00 		0.6719 	
mode 8 	1.26 		0.6285 	
mode 9 	1.76 		0.7056 	
mode 10 	2.16 		0.7016 	
mode 11 	2.31 		0.5394 	
mode 12 	2.39 		0.6059 	
mode 13 	2.46 		0.4075 	
mode 14 	2.66 		0.6475 	
mode 15 	2.75 		0.4065 	
mode 16 	3.06 		0.3493 	
mode 17 	3.08 		0.4132 	
mode 18 	3.19 		0.4097 	
mode 19 	3.26 		0.3425 	
mode 20 	3.37 		0.2945