The human mitochondrial transport protein family: Identification and protein regions significant for transport function and substrate specificity

Protein sequence similarities and predicted structures identified 75 mitochondrial transport proteins (37 subfamilies) from among the 28,994 human RefSeq (NCBI) protein sequences. All, except two, have an E-value of less than 4e−05 with respect to the structure of the single subunit bovine ADP/ATP carrier/carboxyatractyloside complex (bAAC/CAT) (mGenThreader program). The two 30-kDa exceptions have E-values of 0.003 and 0.005. 21 have been functionally identified and belong to 14 subfamilies. A subset of subfamilies with sequence similarities for each of 12 different protein regions was identified. Many of the 12 protein regions for each tested protein yielded different size subsets. The sum of subfamilies in the 12 subsets was lowest for the phosphate transport protein (PTP) and highest for aralar 1. Transmembrane sequences are most unique. Sequence similarities are highest near the membrane center and matrix. They are highest for the region of transmembrane helices H1, H2 and connecting matrix loop 12 and smallest for transmembrane helices H3, H4 and loop 34. These sequence similarities and the predicted high similarities to the bAAC/CAT structure point to common structural/functional elements that could include subunit/subunit contact sites as they have been identified for PTP and AAC. The four residues protein segment (SerLysGlnIle) of loop 12 is the only segment projecting into the center of the funnel-like structure of the bAAC/CAT. It is present in its entirety only in the AACs and with some replacements in the large Ca2+-modulated aspartate/glutamate transporters. Other transporters have deletions and replacements in this region of loop 12. This protein segment with its central location and variation in size and composition likely contributes to the substrate specificity of the transporters.