Difference between revisions of "Fabry:Structure-based mutation analysis"

From Bioinformatikpedia
(Create mutation)
Line 219: Line 219:
   
 
== Create mutation ==
 
== Create mutation ==
  +
=== Pymol ===
SCWRL
 
   
  +
<div style="float:left; border:thin solid lightgrey; margin-right: 20px;">
  +
<figtable id="tab:Pymol">
  +
<caption>Mutagenesis of the 10 selected SNPs performed with pymol. This was done on the basis of a backbone independant library. Usually the rotamer with the least atomic clashes was chosen. Clashes are depicted as red and green disks. The wildtype amino acid is shown in green, the mutated one in red. Hydrogen bonds of the mutant to the surrounding are depicted in blue. If shown, the active site (residues 170 and 231) is colored pink, substrate binding site (position 203-207) cyan and the existing five disulfide bonds (52 ↔ 94, 56 ↔ 63, 142 ↔ 172, 202 ↔ 223, 378 ↔ 382) are highlighted in yellow. Ligand atoms are represented as spheres in dark cyan.</caption>
  +
{| style="border-style: none"
  +
| [[File:FABRY_pymol_A143T.png|right|280px|thumb| A143T]]
  +
| [[File:FABRY_pymol_I289V.png|right|280px|thumb| I289V]]
  +
|-
  +
| [[File:FABRY_pymol_N215S.png|right|280px|thumb| N215S]]
  +
| [[File:FABRY_pymol_P40S.png|right|280px|thumb| P40S]]
  +
|-
  +
| [[File:FABRY_pymol_P323T.png|right|280px|thumb| P323T]]
  +
| [[File:FABRY_pymol_Q279E.png|right|280px|thumb| Q279E]]
  +
|-
  +
| [[File:FABRY_pymol_R118H.png|right|280px|thumb| R118H]]
  +
| [[File:FABRY_pymol_R356W.png|right|280px|thumb| R356W]]
  +
|-
  +
| [[File:FABRY_pymol_S65T.png|right|280px|thumb| S65T]]
  +
| [[File:FABRY_pymol_V316I.png|right|280px|thumb| V316I]]
  +
|-
  +
|}
  +
</figtable>
  +
</div>
  +
  +
Taking a closer look at each of the mutations, which are shown in the pictures in <xr id="tab:Pymol"/> we can underline most of our assumptions posed in section [[Vizualisation |Vizualisation]]. The SNPs shown there were created with Pymol, a backbone independant library was used and we tried to pick that rotamer, that seemed to produce the least atom clashes. Inspite of that, A143T interferes with the ligand and a cysteine, that forms a disulfide bond and thus can be considered as an important site. <br>
  +
In this view, the mutation S65T does not seem to be that malign, since most of the important hydrogen bonds are kept and only little clashes are caused. The same applies again for V316I, Q279E and I289V. The most atom clashes are produced by the introduced Glu at position 279, but in all three point mutations the important h-bonds that stabilize the alpha-helix are still present.<br>
  +
The assumption that Asn 215 has an important role in binding of the ligand can be BESTÄTIGT, but considering the newly introduced hydrogen bonds when mutating it to a Serine, that form a contact to N-Acetyl-D-Glucosamine the mutation might though be not disease causing. Again, both mutations of the prolines (P40S and P323T) seem to be not malign, since there appear hardly any atom clashes and all hydrogen bonds are still present. <br>
  +
The closer look at the SNP R118H reveals that this mutation has to be disease causing, since ZUM EINEN the size and shape of the amino acid is changed completely and also the formed non covalent connections to the adjoint alpha-helices are changed completely. This SCHLUSSFOLGERUNG can be GEZOGEN WERDEN for the mutation of Arginine 356 to Tryptophane, because the crucial bonds to the alpha helix on the right are destroyed by the introduced aromatic rings.
  +
<br style="clear:both">
  +
=== SCWRL ===
  +
  +
<div style="float:left; border:thin solid lightgrey; margin-right: 20px;">
  +
<figtable id="tab:MutSCWRL">
  +
<caption>Mutagenesis of the 10 selected SNPs performed with ''SCWRL''. This was done on the basis of a backbone dependant library. The wildtype amino acid is shown in green, the mutated one in red. The whole structure of the wildtype is colored light blue, the mutant light red. Hydrogen bonds of the mutant to the surrounding are depicted in red, those of the not mutated site in green. </caption>
  +
{| style="border-style: none"
  +
| [[File:FABRY_mutant_A143T_SCWRL.png|right|280px|thumb| A143T]]
  +
| [[File:FABRY_mutant_I289V_SCWRL.png|right|280px|thumb| I289V]]
  +
|-
  +
| [[File:FABRY_mutant_N215S_SCWRL.png|right|280px|thumb| N215S]]
  +
| [[File:FABRY_mutant_P40S_SCWRL.png|right|280px|thumb| P40S]]
  +
|-
  +
| [[File:FABRY_mutant_P323T_SCWRL.png|right|280px|thumb| P323T]]
  +
| [[File:FABRY_mutant_Q279E_SCWRL.png|right|280px|thumb| Q279E]]
  +
|-
  +
| [[File:FABRY_mutant_R118H_SCWRL.png|right|280px|thumb| R118H]]
  +
| [[File:FABRY_mutant_R356W_SCWRL.png|right|280px|thumb| R356W]]
  +
|-
  +
| [[File:FABRY_mutant_S65T_SCWRL.png|right|280px|thumb| S65T]]
  +
| [[File:FABRY_mutant_V316I_SCWRL.png|right|280px|thumb| V316I]]
  +
|-
  +
|}
  +
</figtable>
  +
</div>
  +
<div style="float:left; border:thin solid lightgrey; margin: 0px 20px 0px 0px;">
  +
<figtable id="tab:EnergySCWRL">
  +
<caption>
  +
Minimal energy introduced by the mutated site.<br>
  +
Shown is the energy if only chain A is mutated, as well<br>
  +
as the energy when both sites are mutated.</caption>
  +
{| class="wikitable sortable" style="border-collapse: collapse; border-spacing: 0; border-width: 1px; border-style: solid; padding-left:5px; padding-right:5px; border-color: #000; padding: 0"
  +
! style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 2px 0;"| SNP
  +
! style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 2px 0;"| Minimal energy<br>Chain A
  +
! style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 2px 0;"| Minimal energy<br>Both chains
  +
|-
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| A143T
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| 142.857
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| 290.65
  +
|-
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| I289V
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| 138.101
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| 277.376
  +
|-
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| N215S
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| 143.36
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| 287.975
  +
|-
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| P323T
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| 143.375
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| 287.394
  +
|-
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| P40S
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| 144.203
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| 289.775
  +
|-
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| Q279E
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| 143.082
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| 286.738
  +
|-
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| R118H
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| 142.636
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| 285.207
  +
|-
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| R356W
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| 149.729
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| 305.003
  +
|-
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| S65T
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| 145.551
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| 293.644
  +
|-
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| V316I
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| 148.135
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| 294.203
  +
|-
  +
|}
  +
</figtable>
  +
</div>
  +
  +
The mutations created by ''SCWRL'' are shown in <xr id="tab:MutSCWRL"/>(only chain A was mutated in the structures shown in these pictures), where they are also compared to the wildtype of this residue. Here the focus IS GERICHTET AUF the newly introduced and the deleted hydrogen bonds. The amount of new connections correlates with the size of the minimal energy shown in <xr id="tab:EnergySCWRL"/>, the higher the number, the more new non-covalent bondings are formed and the more severe the changes are. The smallest value is a free energy of 138.101, which WIRD ERZIELT VON the SNP I289V. Looking at the corresponding picture, we see that none of the h-bonds formed by the wildtype is changed by the mutant. The biggest value can be observed when examining the mutation R356W, where all hydrogen bonds on the right hand side are deleted due to the aromatic rings of the Tryptophane. The smallest values for single mutation sites are I289V, R118H and A143T.<br>
  +
Comparing the values of the structures where only one chain was mutated to these where both chains were mutated, a doubling of the values can be observed, but in most cases, the latter is roughly twice as big as the first term, but usually a little bigger.
  +
  +
<br style="clear:both">
   
 
== Comparison energies ==
 
== Comparison energies ==

Revision as of 18:12, 23 June 2012

Fabry Disease » Structure-based mutation analysis


The following analyses were performed on the basis of the α-Galactosidase A sequence. Please consult the journal for the commands used to generate the results.

Preparation

<figtable id="tab:Prep"> All available PDB structures assigned to the Uniprot entry P06280 along with the according Resolution,
Coverage and R-factor. The R-factor was obtained from the PDBsum page. The chosen structure 3S5Y is highlighted.

Entry Method Resolution (Å) Chain Positions (up to 429) R-factor<ref>R value http://www.proteopedia.org/wiki/index.php?title=R_value&oldid=569063, last checked on June 21, 2012</ref> R-free<ref>Free R http://www.proteopedia.org/wiki/index.php?title=Free_R&oldid=1390871, last checked on June 21, 2012</ref> pH PDBsum PDB
1R46 X-ray 3.25 A/B 32-422 0.262 0.301 8.0 ] ]
1R47 X-ray 3.45 A/B 32-422 0.285 0.321 8.0 ] ]
3GXN X-ray 3.01 A/B 32-421 0.239 0.301 4.5 ] ]
3GXP X-ray 2.20 A/B 32-422 0.204 0.265 4.5 ] ]
3GXT X-ray 2.70 A/B 32-422 0.245 0.306 4.5 ] ]
3HG2 X-ray 2.30 A/B 32-422 0.178 0.202 4.6 ] ]
3HG3 X-ray 1.90 A/B 32-426 0.167 0.197 6.5 ] ]
3HG4 X-ray 2.30 A/B 32-423 0.166 0.221 4.6 ] ]
3HG5 X-ray 2.30 A/B 32-422 0.192 0.227 4.6 ] ]
3LX9 X-ray 2.04 A/B 32-422 0.178 0.218 6.5 ] ]
3LXA X-ray 3.04 A/B 32-426 0.216 0.244 6.5 ] ]
3LXB X-ray 2.85 A/B 32-427 0.227 0.264 6.5 ] ]
3LXC X-ray 2.35 A/B 32-422 0.186 0.237 6.5 ] ]
3S5Y X-ray 2.10 A/B 32-422 0.195 0.230 5.1 ] ]
3S5Z X-ray 2.00 A/B 32-421 0.211 0.234 5.1 ] ]
3TV8 X-ray 2.64 A/B 32-422 0.203 0.239 4.6 ] ]

</figtable>

We did not choose the structure 3HG3, although it has the best resolution (1.90 Å) and the second best R-factor (see <xr id="tab:Prep"/>), which is a measure of the agreement between the crystallographic model and the experimental X-ray diffraction data <ref>R-factor (crystallography) (May 17, 2012‎) http://en.wikipedia.org/wiki/R-factor_%28crystallography%29, June 20, 2012</ref>, since it has an Alanin at position 170 (part of the active site) instead of an Aspartic acid. After excluding those structures, that had deviations in the sequence, we had to choose between ten sequences (1R46, 1R47, 3GXN, 3GXP, 3GXT, 3HG2, 3HG4, 3HG5, 3S5Y, 3S5Z) and decided to use 3S5Y. This structure has the advantage of a good pH, very good coverage and still reasonable resolution and R-factor.

Vizualisation

<figure id="fig:allSNPs">

All SNPs mapped onto the structure 3S5Y, chain A. Mutated sites are shown in red. The active site (residues 170 and 231) is shown in pink, substrate binding site (position 203-207) in cyan and the existing five disulfide bonds (52 ↔ 94, 56 ↔ 63, 142 ↔ 172, 202 ↔ 223, 378 ↔ 382) are highlighted in yellow. Ligands are depicted in lines representation in dark cyan.

</figure>

From the first glance, we would consider A143T as disease causing, because it is located right next to a Cysteine, that forms a disulfide bond (yellow) and might interfere with it. The same applies for S65T, which is in the structural neighborhood of this bond and also might cause atom clashes.
There is no mutation that can interfere with the active site(purple) or the substrate binding site (cyan).
The mutation N215S seems to play an important role in the binding of the ligand N-Acetyl-D-Glucosamine.
Both prolines Pro 323 and Pro 40 do not give the impression as if they would be very important. Also Arg 356 is on the surface of the protein and is not involved in any binding, although it might be influencing the nearby helix. This also is applicable for the SNP R118H.
It seems possible, that the mutation V316I does not have a big impact on the helix it is located in.
Neither Q279E nor I289V have an obvious reason to be considered as crucial in this structure.

Create mutation

Pymol

<figtable id="tab:Pymol"> Mutagenesis of the 10 selected SNPs performed with pymol. This was done on the basis of a backbone independant library. Usually the rotamer with the least atomic clashes was chosen. Clashes are depicted as red and green disks. The wildtype amino acid is shown in green, the mutated one in red. Hydrogen bonds of the mutant to the surrounding are depicted in blue. If shown, the active site (residues 170 and 231) is colored pink, substrate binding site (position 203-207) cyan and the existing five disulfide bonds (52 ↔ 94, 56 ↔ 63, 142 ↔ 172, 202 ↔ 223, 378 ↔ 382) are highlighted in yellow. Ligand atoms are represented as spheres in dark cyan.

A143T
I289V
N215S
P40S
P323T
Q279E
R118H
R356W
S65T
V316I

</figtable>

Taking a closer look at each of the mutations, which are shown in the pictures in <xr id="tab:Pymol"/> we can underline most of our assumptions posed in section Vizualisation. The SNPs shown there were created with Pymol, a backbone independant library was used and we tried to pick that rotamer, that seemed to produce the least atom clashes. Inspite of that, A143T interferes with the ligand and a cysteine, that forms a disulfide bond and thus can be considered as an important site.
In this view, the mutation S65T does not seem to be that malign, since most of the important hydrogen bonds are kept and only little clashes are caused. The same applies again for V316I, Q279E and I289V. The most atom clashes are produced by the introduced Glu at position 279, but in all three point mutations the important h-bonds that stabilize the alpha-helix are still present.
The assumption that Asn 215 has an important role in binding of the ligand can be BESTÄTIGT, but considering the newly introduced hydrogen bonds when mutating it to a Serine, that form a contact to N-Acetyl-D-Glucosamine the mutation might though be not disease causing. Again, both mutations of the prolines (P40S and P323T) seem to be not malign, since there appear hardly any atom clashes and all hydrogen bonds are still present.
The closer look at the SNP R118H reveals that this mutation has to be disease causing, since ZUM EINEN the size and shape of the amino acid is changed completely and also the formed non covalent connections to the adjoint alpha-helices are changed completely. This SCHLUSSFOLGERUNG can be GEZOGEN WERDEN for the mutation of Arginine 356 to Tryptophane, because the crucial bonds to the alpha helix on the right are destroyed by the introduced aromatic rings.

SCWRL

<figtable id="tab:MutSCWRL"> Mutagenesis of the 10 selected SNPs performed with SCWRL. This was done on the basis of a backbone dependant library. The wildtype amino acid is shown in green, the mutated one in red. The whole structure of the wildtype is colored light blue, the mutant light red. Hydrogen bonds of the mutant to the surrounding are depicted in red, those of the not mutated site in green.

</figtable>

<figtable id="tab:EnergySCWRL"> Minimal energy introduced by the mutated site.
Shown is the energy if only chain A is mutated, as well
as the energy when both sites are mutated.

SNP Minimal energy
Chain A
Minimal energy
Both chains
A143T 142.857 290.65
I289V 138.101 277.376
N215S 143.36 287.975
P323T 143.375 287.394
P40S 144.203 289.775
Q279E 143.082 286.738
R118H 142.636 285.207
R356W 149.729 305.003
S65T 145.551 293.644
V316I 148.135 294.203

</figtable>

The mutations created by SCWRL are shown in <xr id="tab:MutSCWRL"/>(only chain A was mutated in the structures shown in these pictures), where they are also compared to the wildtype of this residue. Here the focus IS GERICHTET AUF the newly introduced and the deleted hydrogen bonds. The amount of new connections correlates with the size of the minimal energy shown in <xr id="tab:EnergySCWRL"/>, the higher the number, the more new non-covalent bondings are formed and the more severe the changes are. The smallest value is a free energy of 138.101, which WIRD ERZIELT VON the SNP I289V. Looking at the corresponding picture, we see that none of the h-bonds formed by the wildtype is changed by the mutant. The biggest value can be observed when examining the mutation R356W, where all hydrogen bonds on the right hand side are deleted due to the aromatic rings of the Tryptophane. The smallest values for single mutation sites are I289V, R118H and A143T.
Comparing the values of the structures where only one chain was mutated to these where both chains were mutated, a doubling of the values can be observed, but in most cases, the latter is roughly twice as big as the first term, but usually a little bigger.


Comparison energies

foldX

<figtable id="tab:energyStart"> ADD CAPTION HERE

Type Energy
BackHbond -528.41
SideHbond -149.96
Energy_VdW -1013.85
Electro -25.51
Energy_SolvP 1310.89
Energy_SolvH -1318.53
Energy_vdwclash 99.21
energy_torsion 26.73
backbone_vdwclash 478.11
Entropy_sidec 482.20
Entropy_mainc 1219.31
water bonds -29.53
helix dipole -12.51
loop_entropy 0.00
cis_bond 4.50
disulfide -29.93
kn electrostatic -0.98
partial covalent interactions 0.00
Energy_Ionisation 1.20
Entropy Complex 0.00
Total 34.82

</figtable>


Minimise

Gromacs

References

<references />