Team:KAIST-Korea/Project/Modeling

From 2010.igem.org

(Difference between revisions)
(→Result : Figures)
 
(26 intermediate revisions not shown)
Line 8: Line 8:
----
----
<br>
<br>
-
== Single chain antibody structural alignment ==
+
&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.
-
=== Protocol ===
+
-
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).
+
<br>
<br>
<br>
<br>
-
==== Data source ====
+
# Antigen binds to the single chain antibody of modified FGFR and form FGFR-complex(FGFRC).
-
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.
+
# Two FGFR-complexes bind each other and form FGFR-complex-dimer(FGFRCD).
 +
# FGFRCD phosphorylate each other and form FGFR-complex-dimer-phosphorylate (FGFRCDnP)
 +
# FGFRCDnP phosphorylate STAT1 to for STAT1P.
 +
# Two STAT1Ps bind each other and form STAT1P-dimer.
 +
# STAT1P-dimer near cell surface diffuse to inside of nucleus of yeast.
 +
# STAT1P-dimer in nucleus binds to GAS element and activate GFP gene.
<br>
<br>
<br>
<br>
-
==== Single chain antibody synthesis ====
+
&nbsp;&nbsp;And these steps may be formulated as these reactions.
-
We combine the antibody variable region sequences in order of V<sub>H</sub>-linker-V<sub>L</sub> to make single chain antibody sequence. The sequence of linker is GGGGSGGGGS
+
<br>
<br>
 +
<html>
 +
<img src = "https://static.igem.org/mediawiki/2010/b/bd/Modeling1.png" width = 400>
 +
</html>
<br>
<br>
-
 
-
==== Structure Prediction ====
 
-
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.
 
<br>
<br>
 +
&nbsp;&nbsp;And we can expect that the concentration of FGFRCDnP is constant at equilibrium state. The concentration of FGFRCDnP is<br>
 +
<html>
 +
<img src = "https://static.igem.org/mediawiki/2010/a/a6/KAISTModeling2.png" width = 200>
 +
</html>
<br>
<br>
-
==== Structure Alignment ====
 
-
In this step, we check the structural similarity between single chain antibodies 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.
 
<br>
<br>
 +
&nbsp;&nbsp;And the concentration of FGFRCD is
<br>
<br>
-
 
+
<html>
-
=== Result : Table===
+
<img src = "https://static.igem.org/mediawiki/2010/5/5e/KAISTModeling3.png" width = 150>
-
{| class="sortable wikitable" border="0" align="center" style="border: 1px solid #999; background-color:#FFFFFF"
+
</html>
-
|-align="left" bgcolor="#CCCCCC"
+
-
! Name
+
-
!Core Residues||Core RMSD||Raw Score||P-value
+
-
|-
+
-
| 2VXT
+
-
| 95
+
-
| 2.676
+
-
| 84.306
+
-
| 0.000222
+
-
|-
+
-
| 2VXU
+
-
| 98
+
-
| 2.563
+
-
| 92.118
+
-
| 0.0000823
+
-
|-
+
-
| 2VXV
+
-
| 77
+
-
| 2.64
+
-
| 52.269
+
-
| 0.0644
+
-
|-
+
-
| 2ZKH
+
-
| 96
+
-
| 2.855
+
-
| 93.411
+
-
| 0.001
+
-
|-
+
-
| 3AAZ
+
-
| 84
+
-
| 4.197
+
-
| 68.078
+
-
| 0.0642
+
-
|-
+
-
| 3D69
+
-
| 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
+
-
|-
+
-
|}
+
-
* '''Name''' -- Except 16A1(Name of antibody itself), all name of single chain antibodies are came from PDB ID of its source.;
+
<br>
<br>
<br>
<br>
-
 
+
&nbsp;&nbsp;And the concentration of FGFRC is
-
=== Result : Figures===
+
<br>
-
Using arranged pdb files and protein 3D structure drawing program - PyMOL, again, we predict structural similarlity between the structure of FGF binding domain and single chain antibodies associated with the table. In order to show the structures clearly, we control the shape setting 'cartoon' and 'chain'. <font color=#00c5c5><b> Cyan Color </b> </font>is <b>the structure of FGF binding domain</b> and <font color=#25b925><b> Green Color </b> </font> is <b>the comparative antibodies</b>. In the result, we can easily see all antibodies' structure are very similar completely or symmetrically.
+
<html>
 +
<img src = "https://static.igem.org/mediawiki/2010/c/c0/KAISTModeling4.png" width = 200>
 +
</html>
<br>
<br>
<br>
<br>
 +
&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>
<html>
<html>
-
<center><img src="https://static.igem.org/mediawiki/2010/d/d8/Antibody_alignment.jpg"></center>
+
<img src = "https://static.igem.org/mediawiki/2010/0/02/KAISTModeling5.png" width = 500>
</html>
</html>
<br>
<br>
<br>
<br>
-
 
-
===Result Analysis===
 
-
The table above shows the result of alignment of 65 single chain antibody for control group and 1 single chain antibody for our real experiment (16A1). Higher core residue and Raw score and lower core RMSD and p-value means more similar. 16A1 antibody shows 26th highest in core residue, 2nd lowest in core RMSD, 13th highest Raw score and 13th lowest p-value. These result means that our 16A1 antibody has higher similarity than average antibodies.
 

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"