Team:KAIST-Korea/Project/Modeling

From 2010.igem.org

(Difference between revisions)
(→Structure Alignment)
 
(46 intermediate revisions not shown)
Line 1: Line 1:
__NOTOC__
__NOTOC__
{{:Team:KAIST/header}}
{{:Team:KAIST/header}}
-
== Single chain antibody structural alignment ==
+
<table width="100%" border="0" cellpadding="20px">
-
=== Protocol ===
+
<tr>
-
There are four steps to compare structure of single chain antibody and FGF binding domain of FGFR. First step is taking variable region sequences of antibodies. Next step is combining these variable region sequences with linker sequence to make single chain antibody sequence. Third step is predicting the structure of single chain antibody with structure prediction program like modeler. Final step is to structural align these structures of antibodies with structure of FGF binding domain of FGFR(PDB ID : 1EVT).
+
<td valign="top" width="75%">
-
==== Data source ====
+
<span style=font-size:30px><font face="Maiandra GD"><b> Modeling </b></font></span>
-
Single chain antibody is the combination of variable regions of known antibodies with linker sequence which can bind to the antigens. We need to know the V<sub>L</sub> and V<sub>H</sub> sequences to make single chain antibody. The source of these antibody sequences are NCBI, Uniprot and RCSB PDB. NCBI and Uniprot provide the single chain sequence of variable regions (V<sub>L</sub> and V<sub>H</sub>) and antigen binding fragments (F<sub>ab</sub>). RCSB provide the structure of the complexes of antigenbinding fragment which binds to its antigens. But we only need the sequence of variable region. So we get the last 120~150 reside and assume them as the variable region. And data from RCSB contain not only sequence of antibody, but also antigens. S we filter them based on label of files.to get heavy chains and light chains of antibody.
+
<br>
-
 
+
----
-
==== Single chain antibody synthesis ====
+
<br>
-
We combine the antibody variable region sequences with order of V<sub>H</sub>-linker-V<sub>L</sub> to make single chain antibody sequence. The sequence of linker is GGGGSGGGGS
+
&nbsp;&nbsp;Leaving cell transmit signal through chemical modifications of signal-related molecules like phosphorylation, methylation, acetylation or ubiquitination. And these chemical modifications are basically chemical reactions. Therefore, we can assume signal transduction pathway as sequential chemical reactions. We applied this assumption to signal transduction initiated by modified FGFR-antigen binding. Signal transduction follows these steps.
-
 
+
<br>
-
==== Structure Prediction ====
+
<br>
-
We used the Modeller program to predict the structure of single chain antibody from its sequence. Modeller predict 3D structure of protein with structure of know similar proteins based on homology model. Input file is the sequence of single chain antibody with fasta format and output file is the structure of single chain antibody with pdb format.
+
# Antigen binds to the single chain antibody of modified FGFR and form FGFR-complex(FGFRC).
-
 
+
# Two FGFR-complexes bind each other and form FGFR-complex-dimer(FGFRCD).
-
==== Structure Alignment ====
+
# FGFRCD phosphorylate each other and form FGFR-complex-dimer-phosphorylate (FGFRCDnP)
-
In this step, we check the structural similarity between single chain antibody and FGF binding domain to align the structure of single chain antibody with FGF binding domain of FGFR. The structure of FGF binding domain of FGFR is provided by RCSB PDB(PDB ID : 1EVT) We used Matt structural alignment program to do this job. Matt do the structural alignment which minimize the distance between α-carbon chain of two proteins based on the common structure (α helix). Input file is the structure of single chain antibody with pdb format and output file is the text file which contain the number of amino acids which compose shared structure (Core residue), average distance between alpha carbon chains of two proteins (Core RMSD), the score of similarity which is calculated by Matt(Raw score) and the probability that this similarity is just the product of random(p-value) and pdb files which contain the alignment result of single chain antibody with FGF binding domain of FGFR.
+
# FGFRCDnP phosphorylate STAT1 to for STAT1P.
-
 
+
# Two STAT1Ps bind each other and form STAT1P-dimer.
-
=== Result===
+
# STAT1P-dimer near cell surface diffuse to inside of nucleus of yeast.
-
{| class="sortable wikitable" border="0" align="center" style="border: 1px solid #999; background-color:#FFFFFF"
+
# STAT1P-dimer in nucleus binds to GAS element and activate GFP gene.
-
|-align="left" bgcolor="#CCCCCC"
+
<br>
-
! Name
+
<br>
-
!Core Residues||Core RMSD||Raw Score||P-value
+
&nbsp;&nbsp;And these steps may be formulated as these reactions.
-
|-
+
<br>
-
| 2VXT
+
<html>
-
| 95
+
<img src = "https://static.igem.org/mediawiki/2010/b/bd/Modeling1.png" width = 400>
-
| 2.676
+
</html>
-
| 84.306
+
<br>
-
| 0.000222
+
<br>
-
|-
+
&nbsp;&nbsp;And we can expect that the concentration of FGFRCDnP is constant at equilibrium state. The concentration of FGFRCDnP is<br>
-
| 2VXU
+
<html>
-
| 98
+
<img src = "https://static.igem.org/mediawiki/2010/a/a6/KAISTModeling2.png" width = 200>
-
| 2.563
+
</html>
-
| 92.118
+
<br>
-
| 0.0000823
+
<br>
-
|-
+
&nbsp;&nbsp;And the concentration of FGFRCD is
-
| 2VXV
+
<br>
-
| 77
+
<html>
-
| 2.64
+
<img src = "https://static.igem.org/mediawiki/2010/5/5e/KAISTModeling3.png" width = 150>
-
| 52.269
+
</html>
-
| 0.0644
+
<br>
-
|-
+
<br>
-
| 2ZKH
+
&nbsp;&nbsp;And the concentration of FGFRC is
-
| 96
+
<br>
-
| 2.855
+
<html>
-
| 93.411
+
<img src = "https://static.igem.org/mediawiki/2010/c/c0/KAISTModeling4.png" width = 200>
-
| 0.001
+
</html>
-
|-
+
<br>
-
| 3AAZ
+
<br>
-
| 84
+
&nbsp;&nbsp;Therefore, the concentration of FGFRCDnP as the function of concentration of FGFR [FGFR], that of antigen [antigen] and phosphorylation number n is<br>
-
| 4.197
+
<html>
-
| 68.078
+
<img src = "https://static.igem.org/mediawiki/2010/0/02/KAISTModeling5.png" width = 500>
-
| 0.0642
+
</html>
-
|-
+
<br>
-
| 3D69
+
<br>
-
| 137
+
-
| 12.224
+
-
| 98.134
+
-
| 0.00541
+
-
|-
+
-
| 3EO9
+
-
| 80
+
-
| 2.282
+
-
| 84.665
+
-
| 0.00295
+
-
|-
+
-
| 3EOA
+
-
| 116
+
-
| 10.402
+
-
| 98.544
+
-
| 0.0547
+
-
|-
+
-
| 3EOB
+
-
| 68
+
-
| 4.685
+
-
| 56.357
+
-
| 0.3753
+
-
|-
+
-
| 3EYV
+
-
| 85
+
-
| 4.704
+
-
| 66.669
+
-
| 0.1653
+
-
|-
+
-
| 3FMG
+
-
| 81
+
-
| 2.629
+
-
| 73.298
+
-
| 0.0129
+
-
|-
+
-
| 3FOG
+
-
| 90
+
-
| 7.243
+
-
| 74.031
+
-
| 0.0785
+
-
|-
+
-
| 3G6D
+
-
| 72
+
-
| 4.403
+
-
| 69.774
+
-
| 0.1281
+
-
|-
+
-
| 3GBM
+
-
| 119
+
-
| 10.331
+
-
| 115.929
+
-
| 0.0314
+
-
|-
+
-
| 3GBN
+
-
| 77
+
-
| 2.788
+
-
| 57.204
+
-
| 0.0592
+
-
|-
+
-
| 3GHB
+
-
| 101
+
-
| 12.861
+
-
| 81.189
+
-
| 0.3984
+
-
|-
+
-
| 3GHE
+
-
| 96
+
-
| 3.302
+
-
| 109.208
+
-
| 0.00117
+
-
|-
+
-
| 3GI8
+
-
| 127
+
-
| 6.837
+
-
| 113.467
+
-
| 0.00038
+
-
|-
+
-
| 3GI9
+
-
| 52
+
-
| 3.534
+
-
| 43.56
+
-
| 0.166
+
-
|-
+
-
| 3GIZ
+
-
| 89
+
-
| 8.217
+
-
| 68.702
+
-
| 0.3172
+
-
|-
+
-
| 3GK8
+
-
| 91
+
-
| 6.84
+
-
| 84.048
+
-
| 0.00811
+
-
|-
+
-
| 3GKW
+
-
| 74
+
-
| 3.205
+
-
| 47.096
+
-
| 0.1215
+
-
|-
+
-
| 3GNM
+
-
| 135
+
-
| 4.739
+
-
| 103.812
+
-
| 0.00000297
+
-
|-
+
-
| 3GO1
+
-
| 103
+
-
| 8.842
+
-
| 77.631
+
-
| 0.2587
+
-
|-
+
-
| 3GRW
+
-
| 107
+
-
| 9.819
+
-
| 108.405
+
-
| 0.0449
+
-
|-
+
-
| 3H42
+
-
| 69
+
-
| 4.389
+
-
| 59.165
+
-
| 0.3421
+
-
|-
+
-
| 3HC0
+
-
| 67
+
-
| 3.065
+
-
| 40.982
+
-
| 0.1028
+
-
|-
+
-
| 3HC3
+
-
| 67
+
-
| 3.173
+
-
| 49.109
+
-
| 0.0944
+
-
|-
+
-
| 3HC4
+
-
| 36
+
-
| 1.998
+
-
| 42.538
+
-
| 0.206
+
-
|-
+
-
| 3HI5
+
-
| 77
+
-
| 3.222
+
-
| 58.096
+
-
| 0.1832
+
-
|-
+
-
| 3HI6
+
-
| 81
+
-
| 3.646
+
-
| 63.526
+
-
| 0.0394
+
-
|-
+
-
| 3HMW
+
-
| 83
+
-
| 4.159
+
-
| 62.598
+
-
| 0.1976
+
-
|-
+
-
| 3HMX
+
-
| 70
+
-
| 5.215
+
-
| 62.814
+
-
| 0.2228
+
-
|-
+
-
| 3HNT
+
-
| 84
+
-
| 2.62
+
-
| 74.666
+
-
| 0.00118
+
-
|-
+
-
| 3HNV
+
-
| 84
+
-
| 2.998
+
-
| 77.17
+
-
| 0.00374
+
-
|-
+
-
| 3HR5
+
-
| 78
+
-
| 4.158
+
-
| 58.026
+
-
| 0.2978
+
-
|-
+
-
| 3I50
+
-
| 65
+
-
| 2.206
+
-
| 52.252
+
-
| 0.1156
+
-
|-
+
-
| 3I9G
+
-
| 84
+
-
| 5.026
+
-
| 64.155
+
-
| 0.1506
+
-
|-
+
-
| 3IU3
+
-
| 75
+
-
| 3.192
+
-
| 56.006
+
-
| 0.2552
+
-
|-
+
-
| 3IXT
+
-
| 78
+
-
| 2.951
+
-
| 68.627
+
-
| 0.0378
+
-
|-
+
-
| 3KDM
+
-
| 87
+
-
| 4.273
+
-
| 65.722
+
-
| 0.1683
+
-
|-
+
-
| 3KS0
+
-
| 81
+
-
| 3.039
+
-
| 55.562
+
-
| 0.0192
+
-
|-
+
-
| 3KYK
+
-
| 88
+
-
| 4.728
+
-
| 68.681
+
-
| 0.1016
+
-
|-
+
-
| 3KYM
+
-
| 85
+
-
| 6.52
+
-
| 66.102
+
-
| 0.5135
+
-
|-
+
-
| 3L1O
+
-
| 93
+
-
| 2.545
+
-
| 85.909
+
-
| 0.0014
+
-
|-
+
-
| 3L5W
+
-
| 89
+
-
| 4.925
+
-
| 72.178
+
-
| 0.0637
+
-
|-
+
-
| 3L5X
+
-
| 88
+
-
| 4.796
+
-
| 65.262
+
-
| 0.1111
+
-
|-
+
-
| 3L5Y
+
-
| 64
+
-
| 4.689
+
-
| 65.227
+
-
| 0.2007
+
-
|-
+
-
| 3L95
+
-
| 73
+
-
| 5.105
+
-
| 70.312
+
-
| 0.0515
+
-
|-
+
-
| 3LMJ
+
-
| 76
+
-
| 2.699
+
-
| 56.627
+
-
| 0.0237
+
-
|-
+
-
| 3LQA
+
-
| 75
+
-
| 8.371
+
-
| 58.374
+
-
| 0.8315
+
-
|-
+
-
| 3LS4
+
-
| 79
+
-
| 2.458
+
-
| 72.68
+
-
| 0.00694
+
-
|-
+
-
| 3LS5
+
-
| 81
+
-
| 2.443
+
-
| 73.043
+
-
| 0.0014
+
-
|-
+
-
| 3LZF
+
-
| 96
+
-
| 2.468
+
-
| 93.292
+
-
| 0.00155
+
-
|-
+
-
| 3MLR
+
-
| 71
+
-
| 2.574
+
-
| 64.979
+
-
| 0.0165
+
-
|-
+
-
| 3MLS
+
-
| 75
+
-
| 2.534
+
-
| 57.386
+
-
| 0.0202
+
-
|-
+
-
| 3MLU
+
-
| 127
+
-
| 10.025
+
-
| 73.766
+
-
| 0.0305
+
-
|-
+
-
| 3MLV
+
-
| 75
+
-
| 2.454
+
-
| 62.321
+
-
| 0.0218
+
-
|-
+
-
| 3MLW
+
-
| 83
+
-
| 8.94
+
-
| 59.029
+
-
| 0.4846
+
-
|-
+
-
| 3MLX
+
-
| 118
+
-
| 11.485
+
-
| 105.391
+
-
| 0.0961
+
-
|-
+
-
| 3MLY
+
-
| 84
+
-
| 4.704
+
-
| 63.962
+
-
| 0.1023
+
-
|-
+
-
| 3MLZ
+
-
| 82
+
-
| 3.995
+
-
| 64.563
+
-
| 0.123
+
-
|-
+
-
| 3MUG
+
-
| 79
+
-
| 2.827
+
-
| 68.72
+
-
| 0.0223
+
-
|-
+
-
| 3MXV
+
-
| 125
+
-
| 10.857
+
-
| 121.911
+
-
| 0.01
+
-
|-
+
-
| 3MXW
+
-
| 102
+
-
| 11.262
+
-
| 78.348
+
-
| 0.4213
+
-
|--bgcolor="#C3C4D8"
+
-
| 16A1
+
-
| 86
+
-
| 2.139
+
-
| 89.893
+
-
| 0.00443
+
-
|-
+
-
|}
+
-
 
+
-
===Result Analysis===
+

Latest revision as of 05:22, 16 August 2010

 

Modeling


  Leaving cell transmit signal through chemical modifications of signal-related molecules like phosphorylation, methylation, acetylation or ubiquitination. And these chemical modifications are basically chemical reactions. Therefore, we can assume signal transduction pathway as sequential chemical reactions. We applied this assumption to signal transduction initiated by modified FGFR-antigen binding. Signal transduction follows these steps.

  1. Antigen binds to the single chain antibody of modified FGFR and form FGFR-complex(FGFRC).
  2. Two FGFR-complexes bind each other and form FGFR-complex-dimer(FGFRCD).
  3. FGFRCD phosphorylate each other and form FGFR-complex-dimer-phosphorylate (FGFRCDnP)
  4. FGFRCDnP phosphorylate STAT1 to for STAT1P.
  5. Two STAT1Ps bind each other and form STAT1P-dimer.
  6. STAT1P-dimer near cell surface diffuse to inside of nucleus of yeast.
  7. STAT1P-dimer in nucleus binds to GAS element and activate GFP gene.



  And these steps may be formulated as these reactions.


  And we can expect that the concentration of FGFRCDnP is constant at equilibrium state. The concentration of FGFRCDnP is


  And the concentration of FGFRCD is


  And the concentration of FGFRC is


  Therefore, the concentration of FGFRCDnP as the function of concentration of FGFR [FGFR], that of antigen [antigen] and phosphorylation number n is
Retrieved from "http://2010.igem.org/Team:KAIST-Korea/Project/Modeling"