Structural Alignments (Phenylketonuria)
From Bioinformatikpedia
Contents
Summary
...
Explore structural alignments
Dataset generation
Our protein has the CATH Code: 1.10.800.10 (Phenylalanine Hydroxylase)
- reference structure of PAH: 2PAH (96,41% identity)
- identical sequence with filled binding site: 1LRM (--> pdb entry: looked at 3D structure and saw a filled binding site, two ligands: FE and HBI)
- identical sequence with unfilled binding site: not found anyone
- low sequence identity: 3LUY (32,2% - no pdb ID under 30%)
- high sequence identity: pdb ID: 2PHM (89,7%)
- CAT: 1J8U (CATH Code: 1.10.800.10) - there is no other category than 1.10.800.10 for CAT
- CA: 2B5U (CATH Code: 1.10.287.620)
- C: 3BQO (CATH Code: 1.25.40.210)
- other CATH category: 1V8H (CATH Code: 2.60.40.10)
Pymol
...
LGA
LGA ...
SSAP / CATHEDRAL (used by CATH)
SSAP
...
uses Cβ
TopMatch
TopMatch
...
uses Cα
SAP or CE
SAP? ->Error!
->
CE
CE-PDB
...
Modelling scores
To compare the different models, the RMSDs (root-mean-square deviation) are compared. In TopMatch the same formular is taken but called Er (root-mean-square error). The RMSD gives the squared distance between corresponding positions of two superimposed proteins in Ångström. The results are shown in <xr id="rmsd"/>. <figtable id="rmsd">
| RMSD results | |||||||
|---|---|---|---|---|---|---|---|
| Method | 1lrm | 3luy | 2phm | 1j8u | 2b5u | 3bqo | 1vh8 |
| LGA-RMSD | 0.81 | 3.30 | 0.88 | 0.73 | 3.07 | 3.59 | 3.42 |
| SSAP-RMSD | 0.99 | 18.77 | 1.24 | 1.02 | 39.16 | 22.39 | 7.27 |
| TopMatch-Er | 0.60 | 1.98 | 0.81 | 0.63 | 1.21 | 1.12 | 3.25 |
| CE-RMSD | 0.65 | 5.13 | 0.95 | 0.68 | 4.06 | 4.68 | 5.92 |
</figtable>
- lowest RMSDs: TopMatch
- LGA and CE sometimes the one sometimes the other is better. For very similar structures CE better, otherwise LGA???
- worst/highest RMSDs: SSAP, maybe the Cβ are more distant???
- careful about the sidechains: here always 2pah.A as query is taken and xx.A as target
...
Evaluate sequence alignments
Lab journal
...
References
<references/>
