Pseudo amino acid composition

To predict the subcellular localization of proteins and other attributes based on their sequence, two kinds of models are generally used to represent protein samples: (1) the sequential model, and (2) the non-sequential model or discrete model.

The most typical sequential representation for a protein sample is its entire amino acid (AA) sequence, which can contain its most complete information. This is an obvious advantage of the sequential model. To get the desired results, the sequence-similarity-search-based tools are usually utilized to conduct the prediction.

Given a protein sequence P with ${\displaystyle L}$ amino acid residues, i.e.,

• ${\displaystyle \mathbf {P} ={\begin{bmatrix}\mathrm {R} _{1}\mathrm {R} _{2}\mathrm {R} _{3}\mathrm {R} _{4}\mathrm {R} _{5}\mathrm {R} _{6}\mathrm {R} _{7}\cdots \mathrm {R} _{L}\end{bmatrix}}\qquad {\text{(1)}}}$

where R₁ represents the 1st residue of the protein P, R₂ the 2nd residue, and so forth. This is the representation of the protein under the sequential model.

However, this kind of approach fails when a query protein does not have significant homology to the known protein(s). Thus, various discrete models were proposed that do not rely on sequence-order. The simplest discrete model is using the amino acid composition (AAC) to represent protein samples. Under the AAC model, the protein P of Eq.1 can also be expressed by

• ${\displaystyle \mathbf {P} ={\begin{bmatrix}f_{1}&f_{2}&\cdots &f_{20}\end{bmatrix}}^{\mathbf {T} }\qquad {\text{(2)}}}$

where ${\displaystyle \,f_{u}\,(u=1,2,\cdots ,20)}$ are the normalized occurrence frequencies of the 20 native amino acids in P, and T the transposing operator. The AAC of a protein is trivially derived with the protein primary structure known like given in Eq.1; it is also possible by hydrolysis without knowing the exact sequence, and such a step in fact is often a prerequisite for protein sequencing.[3]

Owing to its simplicity, the amino acid composition (AAC) model was widely used in many earlier statistical methods for predicting protein attributes. However, all the sequence-order information is lost. This is its main shortcoming.

To avoid completely losing the sequence-order information, the concept of PseAAC (pseudo amino acid composition) was proposed. [1] In contrast with the conventional amino acid composition (AAC) that contains 20 components with each reflecting the occurrence frequency for one of the 20 native amino acids in a protein, the PseAAC contains a set of greater than 20 discrete factors, where the first 20 represent the components of its conventional amino acid composition while the additional factors incorporate some sequence-order information via various pseudo components.

The additional factors are a series of rank-different correlation factors along a protein chain, but they can also be any combinations of other factors so long as they can reflect some sorts of sequence-order effects one way or the other. Therefore, the essence of PseAAC is that on one hand it covers the AA composition, but on the other hand it contains the information beyond the AA composition and hence can better reflect the feature of a protein sequence through a discrete model.

Meanwhile, various modes to formulate the PseAAC vector have also been developed, as summarized in a 2009 review article.[2]

Figure 1. A schematic drawing to show (a) the 1st-tier, (b) the 2nd-tier, and (c) the 3rd-tier sequence-order-correlation mode along a protein sequence, where R₁ represents the amino acid residue at the sequence position 1, R₂ at position 2, and so forth (cf. Eq.1), and the coupling factors ${\displaystyle J_{i,j}}$ are given by Eq.6. Panel (a) reflects the correlation mode between all the most contiguous residues, panel (b) that between all the 2nd most contiguous residues, and panel (c) that between all the 3rd most contiguous residues.

According to the PseAAC model, the protein P of Eq.1 can be formulated as

• ${\displaystyle \mathbf {P} ={\begin{bmatrix}p_{1},\,p_{2},\,\ldots ,\,p_{20},\,p_{20+1},\,\ldots ,\,p_{20+\lambda }\end{bmatrix}}^{\mathbf {T} },\,\,\,(\lambda <L)\qquad {\text{(3)}}}$

where the (${\displaystyle 20+\lambda }$) components are given by

• ${\displaystyle p_{u}={\begin{cases}{\dfrac {f_{u}}{\sum _{i=1}^{20}f_{i}\,+\,w\sum _{k=1}^{\lambda }\tau _{k}}},&(1\leq u\leq 20)\\[10pt]{\dfrac {w\tau _{u-20}}{\sum _{i=1}^{20}f_{i}\,+\,w\sum _{k=1}^{\lambda }\tau _{k}}},&(20+1\leq u\leq 20+\lambda )\end{cases}}\qquad {\text{(4)}}}$

where ${\displaystyle w}$ is the weight factor, and ${\displaystyle \tau _{k}}$ the ${\displaystyle k}$-th tier correlation factor that reflects the sequence order correlation between all the ${\displaystyle k}$-th most contiguous residues as formulated by

• ${\displaystyle \tau _{k}={\frac {1}{L-k}}\sum _{i=1}^{L-k}\,\mathrm {J} _{i,i+k},\,\,\,(k<L)\qquad {\text{(5)}}}$

with

• ${\displaystyle \mathrm {J} _{i,i+k}={\frac {1}{\Gamma }}\sum _{q=1}^{\Gamma }\left[\Phi _{q}\left(\mathrm {R} _{i+k}\right)-\Phi _{q}\left(\mathrm {R} _{i}\right)\right]^{2}\qquad {\text{(6)}}}$

where ${\displaystyle \Phi _{q}\left(\mathrm {R} _{i}\right)}$ is the ${\displaystyle {q}}$-th function of the amino acid ${\displaystyle \mathrm {R} _{i}\,}$, and ${\displaystyle \Gamma \,}$ the total number of the functions considered. For example, in the original paper by Chou,[1] ${\displaystyle \Phi _{1}\left(\mathrm {R} _{i}\right)}$, ${\displaystyle \Phi _{2}\left(\mathrm {R} _{i}\right)}$ and ${\displaystyle \Phi _{3}\left(\mathrm {R} _{i}\right)}$ are respectively the hydrophobicity value, hydrophilicity value, and side chain mass of amino acid ${\displaystyle \mathrm {R} _{i}\,}$; while ${\displaystyle \Phi _{1}\left(\mathrm {R} _{i+1}\right)}$, ${\displaystyle \Phi _{2}\left(\mathrm {R} _{i+1}\right)}$ and ${\displaystyle \Phi _{3}\left(\mathrm {R} _{i+1}\right)}$ the corresponding values for the amino acid ${\displaystyle \mathrm {R} _{i+1}\,}$. Therefore, the total number of functions considered there is ${\displaystyle \Gamma =3\,}$. It can be seen from Eq.3 that the first 20 components, i.e. ${\displaystyle p_{1},\,p_{2},\,\cdots ,\,p_{20}}$ are associated with the conventional AA composition of protein, while the remaining components ${\displaystyle p_{20+1},\,\cdots ,\,p_{20+\lambda }}$ are the correlation factors that reflect the 1st tier, 2nd tier, …, and the ${\displaystyle \lambda \,}$-th tier sequence order correlation patterns (Figure 1). It is through these additional ${\displaystyle \lambda \,}$ factors that some important sequence-order effects are incorporated.

${\displaystyle \lambda \,}$ in Eq.3 is a parameter of integer and that choosing a different integer for ${\displaystyle \lambda \,}$ will lead to a dimension-different PseAA composition.[4]

Using Eq.6 is just one of the many modes for deriving the correlation factors in PseAAC or its components. The others, such as the physicochemical distance mode[5] and amphiphilic pattern mode,[6] can also be used to derive different types of PseAAC, as summarized in a 2009 review article.[2] In 2011, the formulation of PseAAC (Eq.3) was extended to a form of the general PseAAC as given by:[7]

• ${\displaystyle \mathbf {P} ={\begin{bmatrix}\psi _{1},\,\psi _{2},\,\ldots ,\,\psi _{u},\,\ldots ,\psi _{\Omega }\end{bmatrix}}^{\mathbf {T} }\,\,\,\qquad {\text{(7)}}}$

where the subscript ${\displaystyle \Omega }$ is an integer, and its value and the components ${\displaystyle \psi _{1},\,\psi _{2},\,\ldots }$ will depend on how to extract the desired information from the amino acid sequence of P in Eq.1.

The general PseAAC can be used to reflect any desired features according to the targets of research, including those core features such as functional domain, sequential evolution, and gene ontology to improve the prediction quality for the subcellular localization of proteins.[8][9] as well as their many other important attributes.

• PseAAC web server