Pseudo amino acid composition, or PseAAC, in molecular biology, was originally introduced by Kuo-Chen Chou in 2001 to represent protein samples for improving protein subcellular localization prediction and membrane protein type prediction.[1] Like the vanilla amino acid composition (AAC) method, it characterizes the protein mainly using a matrix of amino-acid frequencies, which helps with dealing with proteins without significant sequential homology to other proteins. Compared to AAC, additional information are also included in the matrix to represent some local features, such as correlation between residues of a certain distance.[2] When dealing the cases of PseAAC, the Chou's invariance theorem has been often used.
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
L
P={\begin{bmatrix}R1R2R3R4R5R6R7 … RL\end{bmatrix}} (1)
where R1 represents the 1st residue of the protein P, R2 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
P={\begin{bmatrix}f1&f2& … &f20\end{bmatrix}}T (2)
where
fu(u=1,2, … ,20)
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. 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]
According to the PseAAC model, the protein P of Eq.1 can be formulated as
P={\begin{bmatrix}p1,p2,\ldots,p20,p20+1,\ldots,p20+λ\end{bmatrix}}T,(λ<L) (3)
where the (
20+λ
pu=\begin{cases} \dfrac{fu}
| 20 | |
| {\sum | |
| i=1 |
fi+
| λ | |
| w\sum | |
| k=1 |
\tauk},&(1\leu\le20) \\[10pt] \dfrac{w\tauu-20
where
w
\tauk
k
k
\tauk=
| 1 | |
| L-k |
| L-k | |
| \sum | |
| i=1 |
Ji,,(k<L) (5)
with
Ji,=
| 1 | |
| \Gamma |
| \Gamma | |
| \sum | |
| q=1 |
\left[\Phiq\left(Ri+k\right)-\Phiq\left(Ri\right)\right]2 (6)
where
\Phiq\left(Ri\right)
{q}
Ri
\Gamma
\Phi1\left(Ri\right)
\Phi2\left(Ri\right)
\Phi3\left(Ri\right)
Ri
\Phi1\left(Ri+1\right)
\Phi2\left(Ri+1\right)
\Phi3\left(Ri+1\right)
Ri+1
\Gamma=3
p1,p2, … ,p20
p20+1, … ,p20+λ
λ
λ
λ
λ
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]
P={\begin{bmatrix}\psi1,\psi2,\ldots,\psiu,\ldots,\psi\Omega\end{bmatrix}}T (7)
where the subscript
\Omega
\psi1,\psi2,\ldots
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.