Sequence Alignments TSD

Sequence searches

There are several alignment methods provided by various initiatives, who tackle the problem of sequence searches. Here some of them are applied for the Hex A protein and analyzed. For the searches non redundant protein databases are used. The outputs are adapted to each other and put in comparison in order to determine the best results. A protocol containing the basic steps taken is available.

Blast

The first sequence similarity search with the Hex A protein was run with Blast. Here the default settings which provide an output of 250 alignments cover a just a small fraction of similar proteins as the e-value of the last hit receives a significantly low e-value of 3e-48. This shows that the sequence search can be continued and more sequences added safely. This is especially important because there are sequences with a comparably low sequence identity of 20% needed for the multiple sequence alignment. The sequence identity correlates with the hit rank of blast, meaning that with a worse sequence identity the e-value is overall expected to increase. To manage between quality deterioration with a worse e-value and on the other hand the need for low sequence identity a limitation of the output sequences was chosen of 1200. Here the e-value does not go beyond 1e-4 and thus the quality of the alignment is still sufficient but there are also sequences aligned with the required low sequence identity.

The results from the BIG80 database contain only Uniprot sequences, which can be explained by the clustering used for big_80 (based on CD_hit), where long sequences are preferred. All hits are unique.

PSI-Blast

The PSI-Blast alignment was assessed in different constellations of e-value cutoff and iteration number for the profile compilation. An appropriate output cutoff was also chosen to avoid irrelevant hits. This threshold was set to 1200 to limit the results to relevant hits just like in the simple Blast search.
After the profile creation from the BIG80 database the search space could be extended.Thus the profiles were used to run Psi-Blast against the larger BIG database. Here the results were restricted to 3800 sequences.

All results are unique hits from the BIG80 as well as the BIG database. In the result set from BIG80 there are only Uniprot sequences whereas in the BIG hits a range between 117-122 proteins are from PDB.
Within the iterations for the profile compilation the proteins are chosen from a pool of ca. 1350 unique proteins within the 10 iteration runs. The 2 iterations consist each of around 1240 unique hits. This shows that in process of each iteration new sequences can be added to the alignments.

The performance for the different combinations of e-value and iterations for the search in the BIG80 database as well as the BIG database are shown in <xr id="tab:psiblast"/>. This comparison depicts the general trend of a higher runtime with a higher number of iterations. As there are other circumstances which also affect the performance such as other concurrent processes on the computing ressources (and the ensuing bottlenecks in I/O and CPU ressources), there is no further reliable inference from the run time results possible.

<figtable id="tab:psiblast">

Iterations	2	2	10	10
E-value	0.002	10E-10	0.002	10E-10
BIG80	3m53	4m3	18m57	21m9
BIG	17m19	11m13	16m39	11m13

Table 1: Different performances of PSI-Blast. </figtable>

... HISTOGRAMM? ...

The sequence overlap from all profile combinations applied to the search in BIG is displayed in <xr id="fig:psi-blast combinations big"/>. To adjuts the different sets for comparability the pdb identifiers were mapped to Uniprot whenever possible. The adjustments of the iteration number and the e-value cut-off clearly accounts for differences in the result set. Each result set expresses an unique intersection with the other sets.

<figure id="fig:psi-blast combinations big">

Overlap of all PSI-Blast runs against the BIG database: The result sets are comprised of 3800 unique proteins. With the mapping to Uniprot there are some additions due to possible ambiguity. 2 Iterations, Evalue: 0.002 = purple; 2 Iterations, Evalue: 10e-10 = yellow; 10 Iterations, Evalue: 0.002 = green; 10 Iterations, Evalue: 10e-10 = blue.

</figure>

HHBlits

The next step was the performance of the sequence search with HHBlits against Uniprot. The results were limited to 460 representatives in accordance with the previous alignments.

Evaluation

<figure id="fig:overlap all">

Overlap of all Uniprot based searchesPSI-Blast search with 10 Iterations .

</figure> <figure id="fig:overlap all">

Overlap of all Uniprot based searchesPSI-Blast search with 2 Iterations .

</figure>