I-sites are short sequence-structure motifs that are mined from the Protein Data Bank (PDB) that correlate strongly with three-dimensional structural elements. These sequence-structure motifs are used for the local structure prediction of proteins. Local structure can be expressed as fragments or as backbone angles. Locations in the protein sequence that have high confidence I-sites predictions may be the initiation sites of folding. I-sites have also been identified as discrete models for folding pathways. I-sites consist of about 250 motifs. Each motif has an amino acid profile, a fragment structure (represented by a "paradigm" fragment chosen from a protein in the PDB) and optionally, a 4-dimensional tensor of pairwise sequence covariance.
Construction of I-site Library The sequence and structure database The database initially consisted of 471 protein sequence families from the HSSP database, with an average of 47 aligned sequences per family. Each family contained a single known structure (parent) from the Brookhaven protein Data Bank. These were a subset of the PDBSelect-25 list, having no more than 25% sequence identity between any two alignments. Disordered loops were omitted. Gaps and insertions in the sequence were ignored. Clustering of sequence segments Each position in the database is described by a weighted amino acid frequency. A similarity measure in sequence space between a segment (p) and a cluster of segments (q) is defined as:
D p q = ∑ i j l o g [ P i j ( p ) + α F i ( 1 + α ) F i ] l o g [ ∑ k ϵ q P i j ( k ) + α ′ F i ( N q + α ′ ) F i ] {\displaystyle D_{pq}=\sum _{ij}log\left[{\dfrac {P_{ij}(p)+\alpha F_{i}}{(1+\alpha )F_{i}}}\right]log\left[{\dfrac {\sum _{k\epsilon q}P_{ij}(k)+\alpha 'F_{i}}{(N_{q}+\alpha ')F_{i}}}\right]} where Pij(p) is the frequency of amino acid i in position j within the segment p. Nq is the number of sequence segments k in the cluster q. Fi is the frequency of amino acid type i in the database overall. The optimal values of a and a0 were determined empirically to be 0.5 and 15, respectively. Using this similarity measure, segments of a given length (3 to 15) were clustered via the k-means algorithm. Assessing structure within a cluster; choice of paradigm The structural similarity between any two peptide segments was evaluated using a combination of the RMS distance matrix error (dme):
d m e = ∑ i = 1 L ∑ j = i − 5 i + 5 ( α i → j s 1 − α i → j s 2 ) 2 N {\displaystyle dme={\sqrt {\dfrac {\sum \limits _{i=1}^{L}\sum \limits _{j=i-5}^{i+5}(\alpha _{i\rightarrow j}^{s1}-\alpha _{i\rightarrow j}^{s2})^{2}}{N}}}} where ai->j is the distance between a-carbon atoms i and j in the segment s1 of length L, and the maximum deviation in backbone torsion angles (mda) over the length of the segment is given by:
… excerpt ends here. Continue reading the full article.
