Difference between revisions of "Task 5: Homology Modeling"

From Bioinformatikpedia
Line 44: Line 44:
 
<xr id="Modeller single"/> lists the selected templates and the Modeller results for the different template structures and alignment methods. In addition to the standard pairwise sequence alignment based on dynamic programming, we also used Modeller's alig2dn() method to improve the alignment by including secondary structure information of the template. The method tries to place gaps outside secondary structure segments. This is especially useful for alignments between more distant related sequences, because those alignments usually contain more gaps than alignments between more close related sequences.
 
<xr id="Modeller single"/> lists the selected templates and the Modeller results for the different template structures and alignment methods. In addition to the standard pairwise sequence alignment based on dynamic programming, we also used Modeller's alig2dn() method to improve the alignment by including secondary structure information of the template. The method tries to place gaps outside secondary structure segments. This is especially useful for alignments between more distant related sequences, because those alignments usually contain more gaps than alignments between more close related sequences.
 
We also tried to improve the alignments manually.
 
We also tried to improve the alignments manually.
  +
 
As Modeller quality score, we chose the [http://salilab.org/modeller/9.11/manual/node253.html DOPE score], which is a statistical potential that was optimized for the assessment of model quality. The DOPE score has an arbitrary scale, but scores for structures of the same protein are comparable and can be used to select the best model from a collection of structures. The lower the score, the better the model.
 
As Modeller quality score, we chose the [http://salilab.org/modeller/9.11/manual/node253.html DOPE score], which is a statistical potential that was optimized for the assessment of model quality. The DOPE score has an arbitrary scale, but scores for structures of the same protein are comparable and can be used to select the best model from a collection of structures. The lower the score, the better the model.
In addition to the DOPE score, we also computed the RMSD and GDT score. The RMSD is a a good measure of the average distance between all pairs of corresponding atoms in two structures. Therefore, the lower the RMSD the better.
+
In addition to the DOPE score, we also computed the RMSD and GDT scores. The RMSD is a a good measure of the average distance between all pairs of corresponding atoms in two structures. Therefore, the lower the RMSD the better.
 
For the GDT score, the average coverage of the target sequence under four defined distance cutoffs is computed. Normally, 1, 2, 4 and 8 Å are used as distance thresholds. The GDT score ranges between 0 and 1, with random superpositions of unrelated structures having a score of 0.1 to 0.2.
 
For the GDT score, the average coverage of the target sequence under four defined distance cutoffs is computed. Normally, 1, 2, 4 and 8 Å are used as distance thresholds. The GDT score ranges between 0 and 1, with random superpositions of unrelated structures having a score of 0.1 to 0.2.
   
Line 56: Line 57:
 
</figtable>
 
</figtable>
   
<xr id="pymol modeller 1QVO"/> shows a visualisation of the two models (purple) created from the template 1QVO_A (red) with the closest homology to the target (green). The first model a) is obviously much better than the second, because its secondary structure features match the target quite good and the position of the alpha helices of the second model differs more from the target.
+
<xr id="pymol modeller 1QVO"/> shows a visualisation of the two models (purple) created from the template 1QVO_A (red) with the closest homology to the target (green). The first model a) is obviously very good ant much better than the second model b), because its secondary structure features match the target quite good and the position of the alpha helices of the second model differs a lot from the target.
   
 
<figtable id="pymol modeller 1S7X">
 
<figtable id="pymol modeller 1S7X">
Line 66: Line 67:
 
</figtable>
 
</figtable>
   
The models created from the template 1S7X_A are worse than those from the more closely related 1QVO_A. The 3D representation in <xr id="pymol modeller 1S7X"/> shows several regions in both models, where the secondary structure elements could not be superimposed correctly to the reference.
+
The models created from the template 1S7X_A are worse than those from the more closely related 1QVO_A. The 3D representation in <xr id="pymol modeller 1S7X"/> shows several regions in both models, where the secondary structure elements could not be superimposed correctly to the reference. The beta sheets in the second model b) match those of the immunoglobuline domain (lower part of the protein) in the target, but there are no beta sheets in that region in the first model a). Although the RMSD is better for model a) , we would agree with the DOPE score that the model b) is better.
 
   
 
<figtable id="pymol modeller 1CD1">
 
<figtable id="pymol modeller 1CD1">
Line 76: Line 76:
 
|}
 
|}
 
</figtable>
 
</figtable>
  +
  +
<xr id="pymol modeller 1CD1"/> shows 3D representations for the models created form the most distant related homolog 1CD1. Both model a) and b) are very different from the target and also the template. The 21% sequence identity between the two proteins is obvious too low to create good alignments. Nevertheless, the secondary structure guided alignment method from Modeller was especially designed for this task: the alignment of sequences with low identity. But
  +
Including the secondary structure information in the alignments did only slightly improve the alignment for the template 1S7X.
  +
The alignments of 1QVO_A and the target show, that including the 2D information did increase the number of gaps and thus decreased the quality of the model.
  +
  +
  +
<figtable id="pymol modeller 1CD1">
  +
{| align="center"
  +
! std alignment || 2D alignment
  +
|-
  +
| >P1;1qvo
  +
GSHSMRYFYTSVSRPGRGEPRFIAVGYVDDTQFVRFDSDAASQRMEPRAPWIEQEGPEYWDQETRNVKAQSQTDR
  +
VDLGTLRGYYNQSEDGSHTIQIMYGCDVGPDGRFLRGYRQDAYDGKDYIALNEDLRSWTAADMAAQITKRKWEAA
  +
HAAEQQRAYLEGRCVEWLRRYLENGKETLQRTDPPKTHMTHHPISDHEATLRCWALGFYPAEITLTWQRDGEDQT
  +
QDTELVETRPAGDGTFQKWAAVVVPSGEEQRYTCHVQHEGLPKPLTLRWE/MIQRTPKIQVYSRHPAENGKSNFL
  +
NCYVSGFHPSDIEVDLLKNGERIEKVEHSDLSFSKDWSFYLLYYTEFTPTEKDEYACRVNHVTLSQPKIVKWDRD
  +
M/QVPLRPMTYK/GSHSMRYFYTSVSRPGRGEPRFIAVGYVDDTQFVRFDSDAASQRMEPRAPWIEQE-GPEYWD
  +
QETRNVKAQSQTDRVDLGTLRGYYNQSEDGSHTIQIMYGCDVGPDGRFLRGYRQDAYDGKDYIALNEDLRSWTAA
  +
DMAAQITKRKWEAAHA-AEQQRAYLEGRCVEWLRRYLENGKETLQRTDPPKTHMTHHPISDHEATLRCWALGFYP
  +
AEITLTWQRDGED-QTQDTELVETRPAGDGTFQKWAAVVVPSGEEQRYTCHVQHEGLPKPLTLRWE/MIQRTPKI
  +
QVYSRHPAENGKSNFLNCYVSGFHPSDIEVDLLKNGERIEKVEHSDLSFSKDWSFYLLYYTEFTPTEKDEYACRV
  +
NHVTLSQPKIVKWDRDM/QVPLRPMTYK*
  +
  +
>P1;1a6z
  +
sequence:1a6z: : : : ::: 0.00: 0.00
  +
---------------------------------------------------------------------------
  +
--------------------------------RLLR---------------------------------------
  +
---------------------------------------------------------------------------
  +
---------------------------------------------------------------------------
  +
---------------------------------------------------------------------------
  +
--------------SHSLHYLFMGASEQDLGLSLFEALGYVDDQLFVFYDHE--SRRVEPRTPWVSSRISSQMWL
  +
QLSQSLKGWDHMFTVDFWTIMENHNHSKE-SHTLQVILGCEMQEDNS-TEGYWKYGYDGQDHLEFCPDTLDWRAA
  +
EPRAWPTKLEWERHKIRARQNRAYLERDCPAQLQQLLELGRGVLDQQVPPLVKVTHH-VTSSVTTLRCRALNYYP
  +
QNITMKWLKDKQPMDAKEFEPKDVLPNGDGTYQGWITLAVPPGEEQRYTCQVEHPGLDQPLIVIW----------
  +
---------------------------------------------------------------------------
  +
----------------------------*
  +
  +
||
  +
  +
>P1;1qvo
  +
---GSHSMRYFYTSVSRPGRGEPRFIAVGYVDDTQFVRFDSDAASQRMEPRAPWIEQE-GPEYWDQETRNVKAQS
  +
QTDRVDLGTLRGYYNQSEDGSHTIQIMYGCDVGPDGRFLRGYRQDAYDGKDYIALNEDLRSWTAADMAAQITKRK
  +
WEAAHAAEQQRAYLEGRCVEWLRRYLENGKETLQRTDPPKTHMTHHPISDHEATLRCWALGFYPAEITLTWQRDG
  +
EDQTQDTELVETRPAGDGTFQKWAAVVVPSGEEQRYTCHVQHEGLPKPLTLRWE/MIQRTPKIQVYSRHPAENGK
  +
SNFLNCYVSGFHPSDIEVDLLKNGERIEKVEHSDLSFSKDWSFYLLYYTEFTPTEKDEYACRVNHVTLSQPKIVK
  +
WDRDM/QVPLRPMTYK/GSHSMRYFYTSVSRPGRGEPRFIAVGYVDDTQFVRFDSDAASQRMEPRAPWIEQEGPE
  +
YWDQETRNVKAQSQTDRVDLGTLRGYYNQSEDGSHTIQIMYGCDVGPDGRFLRGYRQDAYDGKDYIALNEDLRSW
  +
TAADMAAQITKRKWEA-AHAAEQQRAYLEGRCVEWLRRYLENGKETLQRTDPPKTHMTHHPISDHEATLRCWALG
  +
FYPAEITLTWQRDGEDQT-QDTELVETRPAGDGTFQKWAAVVVPSGEEQRYTCHVQHEGLPKPLTLRWE/MIQRT
  +
PKIQVYSRHPAENGKSNFLNCYVSGFHPSDIEVDLLKNGERIEKVEHSDLSFSKDWSFYLLYYTEFTPTEKDEYA
  +
CRVNHVTLSQPKIVKWDRDM/QVPLRPMTYK*
  +
  +
>P1;1a6z
  +
sequence:1a6z: : : : ::: 0.00: 0.00
  +
RLLRSHSLHYLFMGASEQDLGLSLFEALGYVDDQLFVFYD--HESRRVEPRTPWVSSRISSQMWLQLSQSLKGWD
  +
HMFTVDFWTIMENHNHSKE-SHTLQVILG----------------------------------------------
  +
-----------------C---------------------------------------------------------
  +
---------------------------------------------------------------------------
  +
---------------------------------------------------------------------------
  +
------EMQ--E-D-N-ST--------------------------------------------------------
  +
-------------------------------------------E---------GYWKYGYDGQDHLEFCPDTLDW
  +
RAAEPRAWPTKLEWERHKIRARQNRAYLERDCPAQLQQLLELGRGVLDQQVPPLVKVTHHVTSS-VTTLRCRALN
  +
YYPQNITMKWLKDKQPMDAKEFEPKDVLPNGDGTYQGWITLAVPPGEEQRYTCQVEHPGLDQPLIVIW-------
  +
---------------------------------------------------------------------------
  +
-------------------------------*
  +
|+ style="caption-side: bottom; text-align: left" |<font size=1.5>'''Figure 4:''' Standard and 2d alignment of 1QVO_A and 1A6Z_A.
  +
|}
  +
</figtable>
  +
   
   
Including the secondary structure information in the alignment did only improve the model of the most distant homolog 1CD1_A. <xr id="pymol modeller 1CD1"/> shows, that the second model b) is comparable to the models created from 1QVO_A, but one end of model a) in the upper left stands out from the protein.
+
improve the model of the most distant homolog 1CD1_A. <xr id="pymol modeller 1CD1"/> shows, that the second model b) is comparable to the models created from 1QVO_A, but one end of model a) in the upper left stands out from the protein.
 
The alignments between the target sequence and the two more close related template sequences 1QVO and 1S7X are probably already quite good so that including secondary structure information could not improve those alignments.
 
The alignments between the target sequence and the two more close related template sequences 1QVO and 1S7X are probably already quite good so that including secondary structure information could not improve those alignments.
   

Revision as of 14:04, 28 August 2013

lab journal task 5

1A6Z chain A was used as modeling target for all three methods.

Modeller

We used Modeller to create models based on a single template and also multiple templates.

Single template

<css> table.colBasic2 { margin-left: auto; margin-right: auto; border: 2px solid black; border-collapse:collapse; width: 70%; } .colBasic2 th,td { padding: 3px; border: 2px solid black; } .colBasic2 td { text-align:left; } .colBasic2 tr th { background-color:#efefef; color: black;} .colBasic2 tr:first-child th { background-color:#adceff; color:black;} </css>

<figtable id="Modeller single">

Template Seq. identity std Alignment 2d alignment curated Alignment
DOPE score RMSD GDT score DOPE score RMSD GDT score DOPE score RMSD GDT score
1QVO_A 39% -27772 3.647 0.6241 -20467 18.789 0.3410
1S7X_A 29% -19941 15.806 0.3355 -22681 18.101 0.3649
1CD1_A 21% -19034 18.066 0.3640 -18938 21.670 0.2822
Table 1: Template structures and their sequence identity to the target, as computed by Blast. The DOPE score, RMSD and GDT score are given as a quality measure. The different models were created based on pairwise sequence alignments with dynamic programming (std Alignment), pairwise sequence alignment with additional secondary structure information (2d Alignment) and manually curated alignments (curated Alignment).

</figtable>

<xr id="Modeller single"/> lists the selected templates and the Modeller results for the different template structures and alignment methods. In addition to the standard pairwise sequence alignment based on dynamic programming, we also used Modeller's alig2dn() method to improve the alignment by including secondary structure information of the template. The method tries to place gaps outside secondary structure segments. This is especially useful for alignments between more distant related sequences, because those alignments usually contain more gaps than alignments between more close related sequences. We also tried to improve the alignments manually.

As Modeller quality score, we chose the DOPE score, which is a statistical potential that was optimized for the assessment of model quality. The DOPE score has an arbitrary scale, but scores for structures of the same protein are comparable and can be used to select the best model from a collection of structures. The lower the score, the better the model. In addition to the DOPE score, we also computed the RMSD and GDT scores. The RMSD is a a good measure of the average distance between all pairs of corresponding atoms in two structures. Therefore, the lower the RMSD the better. For the GDT score, the average coverage of the target sequence under four defined distance cutoffs is computed. Normally, 1, 2, 4 and 8 Å are used as distance thresholds. The GDT score ranges between 0 and 1, with random superpositions of unrelated structures having a score of 0.1 to 0.2.

<figtable id="pymol modeller 1QVO">

a) classical pairwise sequence alignment
b) inclusion of secondary structure information in the alignment
Figure 1: Superposition of the target 1A6Z_A (green), the template 1QVO_A (red) and the model (purple). Two different alignment methods were used to create the input alignment for Modeller.

</figtable>

<xr id="pymol modeller 1QVO"/> shows a visualisation of the two models (purple) created from the template 1QVO_A (red) with the closest homology to the target (green). The first model a) is obviously very good ant much better than the second model b), because its secondary structure features match the target quite good and the position of the alpha helices of the second model differs a lot from the target.

<figtable id="pymol modeller 1S7X">

a) classical pairwise sequence alignment
b) inclusion of secondary structure information in the alignment
Figure 2: Superposition of the target 1A6Z_A (green), the template 1S7X_A (red) and the model (purple). Two different alignment methods were used to create the input alignment for Modeller.

</figtable>

The models created from the template 1S7X_A are worse than those from the more closely related 1QVO_A. The 3D representation in <xr id="pymol modeller 1S7X"/> shows several regions in both models, where the secondary structure elements could not be superimposed correctly to the reference. The beta sheets in the second model b) match those of the immunoglobuline domain (lower part of the protein) in the target, but there are no beta sheets in that region in the first model a). Although the RMSD is better for model a) , we would agree with the DOPE score that the model b) is better.

<figtable id="pymol modeller 1CD1">

a) classical pairwise sequence alignment
b) inclusion of secondary structure information in the alignment
Figure 3: Superposition of the target 1A6Z_A (green), the template 1CD1_A (red) and the model (purple). Two different alignment methods were used to create the input alignment for Modeller.

</figtable>

<xr id="pymol modeller 1CD1"/> shows 3D representations for the models created form the most distant related homolog 1CD1. Both model a) and b) are very different from the target and also the template. The 21% sequence identity between the two proteins is obvious too low to create good alignments. Nevertheless, the secondary structure guided alignment method from Modeller was especially designed for this task: the alignment of sequences with low identity. But Including the secondary structure information in the alignments did only slightly improve the alignment for the template 1S7X. The alignments of 1QVO_A and the target show, that including the 2D information did increase the number of gaps and thus decreased the quality of the model.


<figtable id="pymol modeller 1CD1">

std alignment 2D alignment
>P1;1qvo

GSHSMRYFYTSVSRPGRGEPRFIAVGYVDDTQFVRFDSDAASQRMEPRAPWIEQEGPEYWDQETRNVKAQSQTDR VDLGTLRGYYNQSEDGSHTIQIMYGCDVGPDGRFLRGYRQDAYDGKDYIALNEDLRSWTAADMAAQITKRKWEAA HAAEQQRAYLEGRCVEWLRRYLENGKETLQRTDPPKTHMTHHPISDHEATLRCWALGFYPAEITLTWQRDGEDQT QDTELVETRPAGDGTFQKWAAVVVPSGEEQRYTCHVQHEGLPKPLTLRWE/MIQRTPKIQVYSRHPAENGKSNFL NCYVSGFHPSDIEVDLLKNGERIEKVEHSDLSFSKDWSFYLLYYTEFTPTEKDEYACRVNHVTLSQPKIVKWDRD M/QVPLRPMTYK/GSHSMRYFYTSVSRPGRGEPRFIAVGYVDDTQFVRFDSDAASQRMEPRAPWIEQE-GPEYWD QETRNVKAQSQTDRVDLGTLRGYYNQSEDGSHTIQIMYGCDVGPDGRFLRGYRQDAYDGKDYIALNEDLRSWTAA DMAAQITKRKWEAAHA-AEQQRAYLEGRCVEWLRRYLENGKETLQRTDPPKTHMTHHPISDHEATLRCWALGFYP AEITLTWQRDGED-QTQDTELVETRPAGDGTFQKWAAVVVPSGEEQRYTCHVQHEGLPKPLTLRWE/MIQRTPKI QVYSRHPAENGKSNFLNCYVSGFHPSDIEVDLLKNGERIEKVEHSDLSFSKDWSFYLLYYTEFTPTEKDEYACRV NHVTLSQPKIVKWDRDM/QVPLRPMTYK*

>P1;1a6z sequence:1a6z:  : :  : ::: 0.00: 0.00



RLLR---------------------------------------



SHSLHYLFMGASEQDLGLSLFEALGYVDDQLFVFYDHE--SRRVEPRTPWVSSRISSQMWL

QLSQSLKGWDHMFTVDFWTIMENHNHSKE-SHTLQVILGCEMQEDNS-TEGYWKYGYDGQDHLEFCPDTLDWRAA EPRAWPTKLEWERHKIRARQNRAYLERDCPAQLQQLLELGRGVLDQQVPPLVKVTHH-VTSSVTTLRCRALNYYP QNITMKWLKDKQPMDAKEFEPKDVLPNGDGTYQGWITLAVPPGEEQRYTCQVEHPGLDQPLIVIW----------



*

>P1;1qvo ---GSHSMRYFYTSVSRPGRGEPRFIAVGYVDDTQFVRFDSDAASQRMEPRAPWIEQE-GPEYWDQETRNVKAQS QTDRVDLGTLRGYYNQSEDGSHTIQIMYGCDVGPDGRFLRGYRQDAYDGKDYIALNEDLRSWTAADMAAQITKRK WEAAHAAEQQRAYLEGRCVEWLRRYLENGKETLQRTDPPKTHMTHHPISDHEATLRCWALGFYPAEITLTWQRDG EDQTQDTELVETRPAGDGTFQKWAAVVVPSGEEQRYTCHVQHEGLPKPLTLRWE/MIQRTPKIQVYSRHPAENGK SNFLNCYVSGFHPSDIEVDLLKNGERIEKVEHSDLSFSKDWSFYLLYYTEFTPTEKDEYACRVNHVTLSQPKIVK WDRDM/QVPLRPMTYK/GSHSMRYFYTSVSRPGRGEPRFIAVGYVDDTQFVRFDSDAASQRMEPRAPWIEQEGPE YWDQETRNVKAQSQTDRVDLGTLRGYYNQSEDGSHTIQIMYGCDVGPDGRFLRGYRQDAYDGKDYIALNEDLRSW TAADMAAQITKRKWEA-AHAAEQQRAYLEGRCVEWLRRYLENGKETLQRTDPPKTHMTHHPISDHEATLRCWALG FYPAEITLTWQRDGEDQT-QDTELVETRPAGDGTFQKWAAVVVPSGEEQRYTCHVQHEGLPKPLTLRWE/MIQRT PKIQVYSRHPAENGKSNFLNCYVSGFHPSDIEVDLLKNGERIEKVEHSDLSFSKDWSFYLLYYTEFTPTEKDEYA CRVNHVTLSQPKIVKWDRDM/QVPLRPMTYK*

>P1;1a6z sequence:1a6z:  : :  : ::: 0.00: 0.00 RLLRSHSLHYLFMGASEQDLGLSLFEALGYVDDQLFVFYD--HESRRVEPRTPWVSSRISSQMWLQLSQSLKGWD HMFTVDFWTIMENHNHSKE-SHTLQVILG----------------------------------------------


C---------------------------------------------------------


EMQ--E-D-N-ST--------------------------------------------------------
E---------GYWKYGYDGQDHLEFCPDTLDW

RAAEPRAWPTKLEWERHKIRARQNRAYLERDCPAQLQQLLELGRGVLDQQVPPLVKVTHHVTSS-VTTLRCRALN YYPQNITMKWLKDKQPMDAKEFEPKDVLPNGDGTYQGWITLAVPPGEEQRYTCQVEHPGLDQPLIVIW-------



*
Figure 4: Standard and 2d alignment of 1QVO_A and 1A6Z_A.

</figtable>


improve the model of the most distant homolog 1CD1_A. <xr id="pymol modeller 1CD1"/> shows, that the second model b) is comparable to the models created from 1QVO_A, but one end of model a) in the upper left stands out from the protein. The alignments between the target sequence and the two more close related template sequences 1QVO and 1S7X are probably already quite good so that including secondary structure information could not improve those alignments.

The DOPE score has a positive correlation with the the GDT score, as well as with the RMSD. The scores differ a bit in some cases, but all in all, they agree that the model created from 1QVO_A and the standard alignment is the best. This is not surprising, since the 3D visualisations show that 1QVO_A is already very similar to the reference, whereas the other two templates differ much more.

Multiple templates

We also created templates using more than one template in a single modeling step. Therefore, we created three sets of structures, one with close homologs, one with distant homologs and one mixed set.

<figtable id="multiple sets">

close homology distant homology mixed
Template Seq. identity Template Seq. identity Template Seq. identity
1QVO_A 39% 3HUJ_C 23% 1QVO_A 39%
1ZAG_A 36% 1CD1_A 21% 1CD1_A 21%
1RJZ_D 34% 1VZY_A 14%
Table 2: The three different sets used as templates for Modeller: two sets of close and distant homologs and a mixed set.

</figtable>

<xr id="multiple sets"/> specifies the three sets.

<figtable id="multiple sets">

close homology distant homology mixed homology
Template 1QVO_A, 1ZAG_A 1QVO_A, 1ZAG_A, 1RJZ_D 3HUJ_C, 1CD1_A 3HUJ_C, 1CD1_A, 1VZY_A 1QVO_A, 1CD1_A
DOPE score -28073 -27460 -25967 -20588 -25894
RMSD 3.432 2.431 4.130 7.741 3.974
GDT score 0.6553 0.7638 0.5607 0.3814 0.5846
Pymol visualisation
Visualisation of the target (green) and the model created from 1QVO_A and 1ZAG_A (purple).
Visualisation of the target (green) and model created from 1QVO_A, 1ZAG_A and 1RJZ_D (purple).
Visualisation of the target (green) and model model created from 3HUJ_C and 1CD1_A (purple).
Visualisation of the target (green) and model created from 3HUJ_C, 1CD1_A and 1VZY_A (purple).
Visualisation of the target (green) and the model created from 1QVO_A and 1CD1_A (purple).
Table 3: Results of the modeling with multiple templates. We computed models using 2 structures as templates and also using three structures.

</figtable>

Including more templates improves the quality of the models. However, we got the best results with two template sequences, because three templates led to a bit worse model than two templates. Surprisingly, two templates with low sequence identity to the target led to a good model with an RMSD of 4.12 which is nearly as good as the model created with the mixed set.

Swiss-Model

We used Swiss-Model to create models using 1QVO_A, 1S7X_A and 1CD1_A as template.

Swiss-Model outputs a raw score and also a Z-score that represents an absolute measure of the model quality. It relates the model's raw score to the scores that high-resolution X-ray structures get and thus gives an estimate of how likely the model has a quality comparable to an experimental structure. A low quality model is indicated by a strong negative Z-score, which means that the raw score is several standard deviations lower as the scores of experimental structures with similar size (see Swiss-Model help).

Swiss-Model also provides plots that help to analyse the local energy of the model. For this, the atomic empirical mean force potential (Anolea) and the Gromos simulation package are used. Both are used calculate the energy of each amino acid in the sequence. The two plots show the protein sequence on the x-axis and the calculated energy of each residue on the y-axis. A low energy corresponds to a favorable energy environment for an amino acid and a positive energy represents an unfavorable energy environment.

<figtable id="swiss-model">

1QVO_A 1S7X_A 1CD1_A
Seq. identity 39% 29% 21%
Z-score -1.977 -2.005 -2.707
RMSD 2.847 2.757 3.604
GDT score 0.6774 0.7086 0.6121
Pymol visualisation
Visualisation of the target (green), the template 1QVO_A and the model (purple).
Visualisation of the target (green), the template 1S7X-A and the model (purple).
Visualisation of the target (green), the template 1CD1_A and the model (purple).
Anolea and Gromos energy
1qvo anolea.png
1s7x anolea.png
1cd1 anolea.png
Table 4: Overview of the Swiss-Model results for the three different templates.

</figtable>

<xr id="swiss-model"/> contains the Swiss-Model results. The models from 1QVO_A and 1S7X_A are both very good having a RMSD of 2.847 and 2.757. The GDT score agrees with the RMSD that the model from 1S7X_A is slightly better than the one from 1QVO_A. However, the Z-score decreases with the sequence identity. Looking at the two models in 3D does not clearly reveal the best model. Some secondary structure segments are more correct in one model and some are better modeled in the oher one. We therefore would rank both models as equally good. The Anolea and Gromos plots for each model are also given in the table. For the first model, they show that the residues at the N and C termini are in a favorable energy environments, but there are some less favorable segments in the middle of the model. Since the first two models are better than the third, we expected to see this also trend represented in the energy landscape. But the energy plot of the third model is not very different from the other two.

I-TASSER

I-Tasser was used to create models from two different templates. Due to I-Tassers very long runtime of over 60h for one protein and because we were only allowed to run one job at a time, we only created 2 models.

<figtable id="i-tasser">

1QVO_A 1CD1_A
Seq. identity 39% 21%
C-score 1.73
RMSD 3.062
GDT score 0.6719
Pymol visualisation
Visualisation of the target (green), the template 1QVO_A and the model (purple).
Visualisation of the target (green), the template 1CD1_A and the model (purple).

</figtable>

I-Tasser uses threading in the first step to search for several template structures with high secondary structure similarity to the target in addition to the user specified template(s). Fragments of those structures are then reassembled to create several models for the target. The models are clustered and the lowest energy models are reported. We only selected the first and best model, because we wanted to make it comparable to the results from Modeller and Swiss-Model, both methods only report a single model.

I-Tasser computes a confidence score (C-score) as quality measure for the created models. It ranges between -5 and 2, with a high score indicating a high confidence (see cscore.txt).

Discussion

The GDT score correates with the RMSD for most models, but there are some exceptions where the GDT score is much worse than the RMSD.

The DOPE score from Modeller is also suited to evaluate the models. lab journal task 5

1A6Z chain A was used as modeling target for all three methods.

Modeller

We used Modeller to create models based on a single template and also multiple templates.

Single template

<css> table.colBasic2 { margin-left: auto; margin-right: auto; border: 2px solid black; border-collapse:collapse; width: 70%; } .colBasic2 th,td { padding: 3px; border: 2px solid black; } .colBasic2 td { text-align:left; } .colBasic2 tr th { background-color:#efefef; color: black;} .colBasic2 tr:first-child th { background-color:#adceff; color:black;} </css>

<figtable id="Modeller single">

Template Seq. identity std Alignment 2d alignment curated Alignment
DOPE score RMSD GDT score DOPE score RMSD GDT score DOPE score RMSD GDT score
1QVO_A 39% -27772 3.647 0.6241 -27169 4.994 0.5653
1S7X_A 29% -19941 15.806 0.3355 -18667 18.099 0.2509
1CD1_A 21% -19034 18.066 0.3640 -24213 5.640 0.4697
Table 1: Template structures and their sequence identity to the target, as computed by Blast. The DOPE score, RMSD and GDT score are given as a quality measure. The different models were created based on pairwise sequence alignments with dynamic programming (std Alignment), pairwise sequence alignment with additional secondary structure information (2d Alignment) and manually curated alignments (curated Alignment).

</figtable>

<xr id="Modeller single"/> lists the selected templates and the Modeller results for the different template structures and alignment methods. In addition to the standard pairwise sequence alignment based on dynamic programming, we also used Modeller's alig2dn() method to improve the alignment by including secondary structure information of the template. The method tries to place gaps outside secondary structure segments. This is especially useful for alignments between more distant related sequences, because those alignments usually contain more gaps than alignments between more close related sequences. We also tried to improve the alignments manually. As Modeller quality score, we chose the DOPE score, which is a statistical potential that was optimized for the assessment of model quality. The DOPE score has an arbitrary scale, but scores for structures of the same protein are comparable and can be used to select the best model from a collection of structures. The lower the score, the better the model. In addition to the DOPE score, we also computed the RMSD and GDT score. The RMSD is a a good measure of the average distance between all pairs of corresponding atoms in two structures. Therefore, the lower the RMSD the better. For the GDT score, the average coverage of the target sequence under four defined distance cutoffs is computed. Normally, 1, 2, 4 and 8 Å are used as distance thresholds. The GDT score ranges between 0 and 1, with random superpositions of unrelated structures having a score of 0.1 to 0.2.

<figtable id="pymol modeller 1QVO">

a) classical pairwise sequence alignment
b) inclusion of secondary structure information in the alignment
Figure 1: Superposition of the target 1A6Z_A (green), the template 1QVO_A (red) and the model (purple). Two different alignment methods were used to create the input alignment for Modeller.

</figtable>

<xr id="pymol modeller 1QVO"/> shows a visualisation of the two models (purple) created from the template 1QVO_A (red) with the closest homology to the target (green). The first model a) is obviously much better than the second, because its secondary structure features match the target quite good and the position of the alpha helices of the second model differs more from the target.

<figtable id="pymol modeller 1S7X">

a) classical pairwise sequence alignment
b) inclusion of secondary structure information in the alignment
Figure 2: Superposition of the target 1A6Z_A (green), the template 1S7X_A (red) and the model (purple). Two different alignment methods were used to create the input alignment for Modeller.

</figtable>

The models created from the template 1S7X_A are worse than those from the more closely related 1QVO_A. The 3D representation in <xr id="pymol modeller 1S7X"/> shows several regions in both models, where the secondary structure elements could not be superimposed correctly to the reference.


<figtable id="pymol modeller 1CD1">

a) classical pairwise sequence alignment
b) inclusion of secondary structure information in the alignment
Figure 3: Superposition of the target 1A6Z_A (green), the template 1CD1_A (red) and the model (purple). Two different alignment methods were used to create the input alignment for Modeller.

</figtable>


Including the secondary structure information in the alignment did only improve the model of the most distant homolog 1CD1_A. <xr id="pymol modeller 1CD1"/> shows, that the second model b) is comparable to the models created from 1QVO_A, but one end of model a) in the upper left stands out from the protein. The alignments between the target sequence and the two more close related template sequences 1QVO and 1S7X are probably already quite good, so that including secondary structure information could not improve those alignments.

The DOPE score has a positive correlation with the the GDT score, as well as with the RMSD. The scores differ a bit in some cases, but all in all, they agree that the model created from 1QVO_A and the standard alignment is the best. This is not surprising, since the 3D visualisations show that 1QVO_A is already very similar to the reference, whereas the other two templates differ much more.


Multiple templates

We also created templates using more than one template in a single modeling step. Therefore, we created three sets of structures, one with close homologs, one with distant homologs and one mixed set.

<figtable id="multiple sets">

close homology distant homology mixed
Template Seq. identity Template Seq. identity Template Seq. identity
1QVO_A 39% 3HUJ_C 23% 1QVO_A 39%
1ZAG_A 36% 1CD1_A 21% 1CD1_A 21%
1RJZ_D 34% 1VZY_A 14%
Table 2: The three different sets used as templates for Modeller: two sets of close and distant homologs and a mixed set.

</figtable>

<xr id="multiple sets"/> specifies the three sets.

<figtable id="multiple sets">

close homology distant homology mixed homology
Template 1QVO_A, 1ZAG_A 1QVO_A, 1ZAG_A, 1RJZ_D 3HUJ_C, 1CD1_A 3HUJ_C, 1CD1_A, 1VZY_A 1QVO_A, 1CD1_A
DOPE score -28073 -27460 -25967 -20588 -25894
RMSD 3.432 2.431 4.130 7.741 3.974
GDT score 0.6553 0.7638 0.5607 0.3814 0.5846
Pymol visualisation
Visualisation of the target (green) and the model created from 1QVO_A and 1ZAG_A (purple).
Visualisation of the target (green) and model created from 1QVO_A, 1ZAG_A and 1RJZ_D (purple).
Visualisation of the target (green) and model model created from 3HUJ_C and 1CD1_A (purple).
Visualisation of the target (green) and model created from 3HUJ_C, 1CD1_A and 1VZY_A (purple).
Visualisation of the target (green) and the model created from 1QVO_A and 1CD1_A (purple).
Table 3: Results of the modeling with multiple templates. We computed models using 2 structures as templates and also using three structures.

</figtable>

Including more templates improves the quality of the model, too many templates

Swiss-Model

We used Swiss-Model to create models from the same templates as the ones used for the were used for Modeller, namely 1QVO_A, 1S7X_A and 1CD1_A.

Swiss-Model outputs a raw score and also a Z-score that represents an absolute measure of the model quality. It relates the model's raw score to the scores that high-resolution X-ray structures get and thus gives an estimate of how likely the model has a quality comparable to an experimental structure. A low quality model is indicated by a strong negative Z-score, which means that the raw score is several standard deviations lower as the scores of experimental structures with similar size (see Swiss-Model help).

Swiss-Model also provides plots that help to analyse the local energy of the model. For this, the atomic empirical mean force potential (Anolea) and the Gromos simulation package are used. Both are used calculate the energy of each amino acid in the sequence. The two plots show the protein sequence on the x-axis and the calculated energy of each residue on the y-axis. A low energy corresponds to a favorable energy environment for an amino acid and a positive energy represents an unfavorable energy environment.

<figtable id="swiss-model">

1QVO_A 1S7X_A 1CD1_A
Seq. identity 39% 29% 21%
Z-score -1.977 -2.005 -2.707
RMSD 2.847 2.757 3.604
GDT score 0.6774 0.7086 0.6121
Pymol visualisation
Visualisation of the target (green), the template 1QVO_A and the model (purple).
Visualisation of the target (green), the template 1S7X-A and the model (purple).
Visualisation of the target (green), the template 1CD1_A and the model (purple).
Anolea and Gromos energy
1qvo anolea.png
1s7x anolea.png
1cd1 anolea.png
Table 4: Overview of the Swiss-Model results for the three different templates.

</figtable>

The Swiss-Model results are summarised in <xr id="swiss-model"/>. The best model is the one from


I-TASSER

I-Tasser was used to create models from two different templates. Due to I-Tassers very long runtime of over 60h for one protein and because we were only allowed to run one job at a time, we only created 2 models.

<figtable id="i-tasser">

1QVO_A 1CD1_A
Seq. identity 39% 21%
C-score 1.73
RMSD 3.062
GDT score 0.6719
Pymol visualisation
Visualisation of the target (green), the template 1QVO_A and the model (purple).
Visualisation of the target (green), the template 1CD1_A and the model (purple).

</figtable>

I-Tasser uses threading in the first step to search for several template structures with high secondary structure similarity to the target in addition to the user specified template(s). Fragments of those structures are then reassembled to create several models for the target. The models are clustered and the lowest energy models are reported. We only selected the first and best model, because we wanted to make it comparable to the results from Modeller and Swiss-Model, both methods only report a single model.

I-Tasser computes a confidence score (C-score) as quality measure for the created models. It ranges between -5 and 2, with a high score indicating a high confidence (see cscore.txt).

Discussion

The GDT score correates with the RMSD for most models, but there are some exceptions where the GDT score is much worse than the RMSD.

All methods were ableo to create good models, although ??? created the best. The models from Modeller have the highest RMSD and lowes GDT score and are thus the worst.