Email updates

Keep up to date with the latest news and content from BMC Bioinformatics and BioMed Central.

Open Access Highly Accessed Research article

Lipid exposure prediction enhances the inference of rotational angles of transmembrane helices

Jhih-Siang Lai1, Cheng-Wei Cheng1, Allan Lo2*, Ting-Yi Sung1* and Wen-Lian Hsu1

Author affiliations

1 Institute of Information Science, Academia Sinica, Taipei, Taiwan

2 Verinata Health, an Illumina Company, Redwood City, CA 94063, USA

For all author emails, please log on.

Citation and License

BMC Bioinformatics 2013, 14:304  doi:10.1186/1471-2105-14-304


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2105/14/304


Received:24 May 2013
Accepted:1 October 2013
Published:11 October 2013

© 2013 Lai et al.; licensee BioMed Central Ltd.

This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Since membrane protein structures are challenging to crystallize, computational approaches are essential for elucidating the sequence-to-structure relationships. Structural modeling of membrane proteins requires a multidimensional approach, and one critical geometric parameter is the rotational angle of transmembrane helices. Rotational angles of transmembrane helices are characterized by their folded structures and could be inferred by the hydrophobic moment; however, the folding mechanism of membrane proteins is not yet fully understood. The rotational angle of a transmembrane helix is related to the exposed surface of a transmembrane helix, since lipid exposure gives the degree of accessibility of each residue in lipid environment. To the best of our knowledge, there have been few advances in investigating whether an environment descriptor of lipid exposure could infer a geometric parameter of rotational angle.

Results

Here, we present an analysis of the relationship between rotational angles and lipid exposure and a support-vector-machine method, called TMexpo, for predicting both structural features from sequences. First, we observed from the development set of 89 protein chains that the lipid exposure, i.e., the relative accessible surface area (rASA) of residues in the lipid environment, generated from high-resolution protein structures could infer the rotational angles with a mean absolute angular error (MAAE) of 46.32˚. More importantly, the predicted rASA from TMexpo achieved an MAAE of 51.05˚, which is better than 71.47˚ obtained by the best of the compared hydrophobicity scales. Lastly, TMexpo outperformed the compared methods in rASA prediction on the independent test set of 21 protein chains and achieved an overall Matthew’s correlation coefficient, accuracy, sensitivity, specificity, and precision of 0.51, 75.26%, 81.30%, 69.15%, and 72.73%, respectively. TMexpo is publicly available at http://bio-cluster.iis.sinica.edu.tw/TMexpo webcite.

Conclusions

TMexpo can better predict rASA and rotational angles than the compared methods. When rotational angles can be accurately predicted, free modeling of transmembrane protein structures in turn may benefit from a reduced complexity in ensembles with a significantly less number of packing arrangements. Furthermore, sequence-based prediction of both rotational angle and lipid exposure can provide essential information when high-resolution structures are unavailable and contribute to experimental design to elucidate transmembrane protein functions.

Background

Integral membrane proteins participate in diverse cellular functions such as signal transductions, bioenergetics, ion transport, cell adhesion, and cell-cell recognition. It has also been estimated that about 20-30% of a typical genome encode for proteins with a transmembrane (TM) domain [1,2]. Despite their biological importance and abundance, the mechanism by which TM proteins fold into native structures remains poorly understood due to a limited number of solved structures, accounting for less than 1% of all deposited structures in the Protein Data Bank (PDB) [3]. Therefore, computational methods play an important role in deciphering the sequence-to-structure relationships and advancing our knowledge in this particular class of proteins.

Though recently solved structures of several amino-acid transporters, e.g., eukaryotic CLC Transporter (coded as 3ORG [4] in PDB) and potassium ion transporter (coded as 3PJZ [5] in PDB) revealed the existence of short helices in the reentrant region [6], the canonical topologies of TM proteins can be viewed as pairs of interacting transmembrane helices (TMHs), connecting loops and extramembraneous domains. In particular, the interaction between TMHs is an important determinant of folding and stability by the proposed two-stage model [7,8]. Such an interaction is mediated by structural contacts at the helical interfaces with the protein itself, the ligands, as well as the lipid environment. From the perspective of structural modeling, the rotational angle of a TMH is a strong determinant of its interacting faces with the rest of the protein structure and the lipids. At the stage of conformation space sampling, we could filter out decoys that severely deviate from the predicted rotational angles. To elucidate rotational angles, Eisenberg et al. [9,10] showed that hydrophobic scales can be used to estimate the hydrophobic moment direction to approximate the lipid-facing direction and proposed equations to calculate the rotational angles of TMHs based on such property. Later, several hydrophobicity scales or propensities have been proposed [11-15] to predict exposed residues or faces. To the best of our knowledge, there have been few advances to use lipid exposure, specifically the relative accessible surface area (rASA) in the lipid environment, to predict rotational angles. Henceforth, we use rASA for convenience to represent rASA in the lipid environment since in this paper we focus on the residues in such environment.

To determine the rotational angle of each TMH in the tertiary structure of a TM protein requires the information of lipid-facing direction, which has been defined differently in the literature. Pilpel et al. [12] described the lipid-facing direction as the vector opposite to the bisector of the acute angle formed by the two lines from the geometric center of the target TMH pointing to the geometric centers of the two nearest TMHs in the whole molecule. Stevens and Arkin [16] defined lipid-facing direction of a TMH as the vector connecting the geometric centers of the target helix and the whole molecule. The molecule could be a single chain or the complete protein. Dastmalchi et al. [15] accepted both the above definitions for lipid-facing direction.

Lipid exposure of each TMH has been shown useful in distinguishing between surfaces and interior interfaces, identifying potential functional residues buried in the protein core and exposed residues for protein-protein interactions [14,17], and therefore facilitates the prediction of helix-packing conformations [17,18]. In this work, we propose to use predicted rASA to estimate lipid-facing direction and then determine rotational angles of TMHs.

In order to train machine learning models by observed rASA from solved structures, we calculated the accessible surface area (ASA) by rolling a spherical probe along the van der Waal’s (VDW) surface of a protein molecule [19]. The observed rASA for each residue is defined as dividing ASA by its reference value in an extended Gly-X-Gly tripeptide conformation. Several methods have been proposed to predict ASA instead of rASA in the lipid environment. For example, ASAP [20] uses evolutionary profiles to predict solvent accessibility by the support vector regression (SVR) and reports a Pearson correlation coefficient (PCC) of 0.62 for ASA prediction. The MPRAP [21] uses the support vector machine (SVM) to predict ASA of complete TM proteins by evolutionary profiles and residue distance from membrane center [22]. A number of methods have been proposed to predict the burial or exposed status of each TM residue, where the status is determined by rASA whether below a predefined threshold. Several methods predict the status of TM residues directly from sequences without predicting their ASA; most of them rely on sequence conservation or knowledge-based propensities, including kPROT [12], ProperTM [11], TMLIP [13], MO [14] and TMX [23]. RHYTHM [24] predicts the burial status of TM residues by matrix-based helix-helix contact prediction and sequence conservation, but this method requires prior knowledge of TMHs such as membrane coils or transporter/channels.

In this paper, we present TMexpo, a method to predict rotational angles of TMHs. For each TM residue, TMexpo first predicts rASA by SVR and predicts the burial or exposed status by SVM; both models use evolutionary profiles, sequence conservation, helix insertion energy and biochemical properties as features. Next, TMexpo determines rotational angles of TMHs based on the predicted rASA. In rotational angle prediction, TMexpo outperformed predictors using hydrophobicity scales and propensities by at least 19.2˚ in terms of mean absolute angular error (MAAE) on the independent test set of 21 protein chains. Notably, the prediction results showed that rotational angles of TMHs could be better inferred by predicted rASA than by existing scales. We expect the rotational angle prediction could benefit the structure prediction, especially for free modeling of transmembrane protein structures due to its difficulty and the necessities of reducing the number of packing arrangements.

Results and discussion

Observed rASA can better infer the lipid-facing direction than hydrophobicity scales and lipid-facing propensities

The packing mechanism in which TM protein assembles in the lipids is not fully understood. One significant advance in this area is seen in the mechanism of the Sec translocon which demonstrates how TM proteins enter the membrane [25-27]. A commonly accepted model of TM protein folding is the two-stage model [7]. Later, White and colleagues extended the two-stage model by including a four-step thermodynamic cycle of folding energy description [28]. Several scales and propensities were developed to understand the lipid-facing direction of TMHs based on sequence analysis or knowledge-based information; however, using these scales to predict the rotational angle of TM proteins often results in low prediction accuracy. Therefore, we investigated whether an environment descriptor such as lipid exposure could better infer a geometric parameter such as the rotational angles of TMHs.

To gain insights into the relationship between lipid exposure and rotational angles, we analyzed a dataset of 110 multi-spanning protein chains with high-resolution structures, which were divided into a development set of 89 protein chains and an independent test set of 21 protein chains (described in Methods). We first used the NACCESS program [19,29] to determine the observed rASA of each TM residue. We then calculated the rotational angle of each TMH by the rASA moment (Equations 6, 7, 8, 9 described in Methods) derived from the observed rASA and compared it with the rotational angle calculated from existing hydrophobicity scales or propensities. The comparison results on the dataset in terms of MAAE are shown in Table 1. It can be observed from Table 1 that the rotational angle calculated from the observed rASA achieves the best MAAE, outperforming the propensity-based methods by at least 25.15˚ on the development set and 21.96˚ on the independent test set. It demonstrates that the observed rASA could better infer the rotational angles of TMHs.

Table 1. The MAAE of rotational angles determined by various approaches

The success of determining the rotational angle via the observed rASA calculated by NACCESS may be due to the following two reasons. First, the description of observed rASA is derived from known protein structures, but the descriptions of hydrophobicity and lipid-facing propensities are derived from the sequence. Therefore, hydrophobicity and lipid-facing propensities alone are insufficient for accurate inference of helical packing, thereby rendering worse rotational angle estimation. Second, researchers had a simplified view of membrane proteins being “inside-out” proteins, which have interior polar core and exterior apolar surface [30,31]. However, this paradigm was challenged since biased distribution of hydrophobic residues could not be detected in every membrane protein [32]. As more solved structures become available, statistical analyses on these structures also support the above finding [33-36]. The canonical view of the “inside-out” property of membrane proteins based on hydrophobicity is challenged. On the contrary, observed rASA is a structural environment descriptor of helical packing and better infer the rotational angle.

Relative accessible surface area predicted by TMexpo can also better infer the lipid-facing direction than hydrophobicity scales and lipid-facing propensities

To evaluate the capability of our proposed method for TM proteins with unknown structures, we used rASA predicted by TMexpo to determine the rotational angles of TMHs. Then we evaluated how predicted rASA could infer the rotational angle in terms of MAAE by the following two experiments. First, we used the development set of 89 chains to develop the TMexpo prediction model and tested on the independent test set of 21 chains. Second, we performed leave-one-out cross validation (LOOCV) on the development set. The results of both experiments are shown in Table 1. Both rASA prediction results of the development set and the independent test set are provided in the Additional file 1: Dataset S1. Particularly, detailed prediction results of 188 TMHs in the independent test set are reported in the Additional file 2: Table S1. Notably, 44 of 188 TMHs had the angular error less than 15˚ and their rotational angles were predicted precisely by TMexpo.

Additional file 1: Dataset S1. The rASA prediction results of the development set and the independent test set.

Format: ZIP Size: 1.5MB Download fileOpen Data

Additional file 2: Table S1. Sequence, observed angle, predicted angle, moment lengths and MAAE of the 188 TMHs in the independent test set of 21 proteins.

Format: DOCX Size: 52KB Download fileOpen Data

Observed from the first experiment, TMexpo-predicted rASA is shown to be comparable to the observed rASA derived from NACESS for inferring the rational angles; the predicted rASA achieved an MAAE of 48.31˚, slightly worse than the MAAE of 45.55˚ achieved by the observed rASA. Nevertheless, TMexpo-predicted rASA achieved much better MAAE than the other predictors using different hydrophobicity scales and propensities, including Eisenberg et al.’s consensus hydrophobicity scale (ES) [10], kPROT [12], ProperTM [11], TMLIP [13], and MO [14] by at least 19.2˚. The second experiment reported consistent results with the first experiment. Specifically, performing LOOCV on the development set resulted in the MAAE of 51.05˚, slightly worse than that inferred by the observed rASA, but better than the compared predictors by at least 20.42˚. Our results show that without known structures, TMexpo can effectively infer the rotational angle within a close margin to that inferred by the observed rASA and improve the prediction compared to hydrophobicity scales or lipid-facing propensities.

Comparison of relative accessible surface area prediction methods

The proposed method TMexpo is capable of predicting not only the rotational angles of TMHs but also rASA of TM residues. We evaluated real-number rASA and binary classification of burial and exposed status predicted by TMexpo on the independent test set of 21 protein chains. At first, we demonstrated the performances of TMexpo in which TM residues participate in interchain contacts were excluded from the independent test set. The results are shown in Table 2, where exposed residues defined by observed rASA ≥5% are regarded as positive cases for evaluating the classification performance. TMexpo achieved an overall PCC of 0.66, mean absolute error (MAE) of 0.12 and root mean squared error (RMSE) of 0.16 in predicting rASA. With respect to the binary classification, TMexpo achieved an overall MCC, accuracy, sensitivity, specificity, and precision of 0.51, 75.26%, 81.30%, 69.15%, and 72.73%, respectively.

Table 2. Comparison of different methods for classifying exposed/buried residues on the independent test set without interface TM residues

Next, we compared TMexpo’s performance of predicting rASA with existing methods, including TMX [23], RHYTHM [24] and MPRAP [21]. One distinction among these methods is that RHYTHM requires prior knowledge of protein types as membrane-coil or channel for prediction. For comparison, we retrieved the prediction results for TMX, RHYTHM and MPRAP from their web servers by using their default parameters. Table 2 shows that TMexpo outperforms the compared methods across most of the measures except a slightly lower specificity compared to MPRAP by 1.16%. The specificity of TMexpo is lower than its sensitivity by 12.15%, and most of the predictors except MPRAP have the same trend as TMexpo’s results. This observation implies that the detection of buried residues may be more difficult than that of exposed residues in the TM domains. To further gain insights into this issue, we extracted buried residues in our dataset of 110 chains from the helix-packing database TMPad [37], and found that over 77% residues have at least one interhelical contact. This suggests that prediction of buried residues could very likely be improved by detecting contacts of interhelical interactions; however, TMexpo and most of the compared predictors retrieve features by local information of subsequences only. Interestingly, interhelical contacts may be conserved in sequences and discovered from evolutionary information such as PSSM profiles [38]. We consider that evolutionary information is an effective feature for capturing interhelical interactions that contributes to rASA/burial status prediction; however, interhelical interaction prediction is still a challenging problem.

Moreover, since we excluded TM residues that participate in interchain contacts, i.e., interface TM residues, in the training of TMexpo (as described in Model Development of the subsection entitled “An SVM-based predictor for lipid exposure of TM helices” in Methods), we also compared the methods for classifying burial status of each interface TM residue in the independent test set. After removing TM residues with missing atoms, there are 392 interface TM residues out of the remained entire 3,553 TM residues in the independent test set of 21 protein chains. Since observed rASA derived from subunit structure and complete structure, respectively, can be different for interface TM residues, we used both the observed rASA to determine the “true” status of each residue. We reported the performance of each method based on the “true” observed rASA for both the 392 interface TM residues and the 3,553 entire TM residues in Table 3 and Table 4, respectively. The tables demonstrate that most rASA predictors, except MPRAP, have lower specificity than sensitivity. For the 392 interface TM residues, most predictors have better MCC calculated from complete structure than that calculated from subunit structure. Furthermore, Table 4 shows that TMexpo achieved slightly better MCC, sensitivity, and accuracy on the 3,553 entire TM residues.

Table 3. Comparison of different methods for classifying exposed/buried residues on 392 interface TM residues of the independent test set by rASA derived from both subunit structure and complete structure (in parentheses)

Table 4. Comparison of different methods for classifying exposed/buried residues on 3,553 entire TM residues of the independent test set by rASA derived from both subunit structure and complete structure (in parentheses)

Rotational angle can help in determining helical packing in transmembrane proteins

Harrington and Ben-Tal [39] characterized five structural features of interhelical interactions, namely, aromatic interactions, hydrogen bonds, salt bridges, and two interactions from packing motifs, that are useful for helical packing. They proposed an algorithm to pack the TMHs of TM proteins, as follows: First, the algorithm ordered the TMHs by the sequence from the N-terminus, and then iteratively grouped sequential TMHs by a scoring function based on the five types of interactions. They demonstrated helical packing on 15 diverse proteins, and the average RMSD of Cα in the native structure of the 15 reconstructed TM proteins ranged from 0.51 Å to 1.35 Å. In this subsection, we reexamined these proteins to study the rotational angle of TMHs and its relationship to helical packing. Since the protein 1AFO discussed in their work contains only one TMH, we excluded this protein in our analysis.

Table 5 shows angular error of the 14 proteins. The MAAE of all 73 TMHs of the 14 proteins is 41.04˚, and 60% of all TMHs were predicted with MAAE ≤43˚. The worst MAAE (113.71˚) was found in protein 2OAR, and the best case is 2UUH (10.47˚). Out of all 73 TMHs, 13 TMHs were predicted with angular errors less than 10˚. With further investigation, we observed some common features from TMPad. First, most of TMHs are linear or curved, and three TMHs are slightly kinked in their structures. Second, we observed strong interhelical interactions constrained in a single helical interface in most of these cases. For example, in the acid-sensing ion channel 2QTS:A [40], its TMH1 (curved) shares 7 contact pairs with TMH2. The second example is an avian mitochondrial complex II (2H88:C) [41], where its TMH3 (curved) shares 13 contact pairs with TMH2. Lastly, in the structure of aquaporin-0 (2B6O:A) [42], its TMH6 (linear) associates with TMH4 via a dense cluster of 16 contacts. Based on the above findings, regular TMHs with strong and constrained interhelical interactions in a single helical interface reveal periodicity in rASA of their residues, thus making rotational angle prediction highly accurate.

Table 5. Rotational angle prediction on the protein chains from Harrington and Ben-Tal’s work

On the other hand, since rotational angle prediction is strongly correlated with the periodicity of a helix, predicted rotational angles may not work well for packing of helices that deviate from regular periodicities of rASA, such as those severely kinked, disrupted, highly tilted, or associated with a reentrant loop. Six out of all 73 TMHs were poorly predicted with angular errors over 100.29˚. We observed these six TMHs being classified as kinked or containing partial non-helical structure in the TM domain, and therefore the moment-based prediction performed poorly.

Rotational angles prediction based on both predicted topology and predicted relative accessible surface area

To illustrate the capability of the proposed method for rotational angle prediction based on predicted topology information, we submitted sequences of independent test set to three web servers including SVMSignal [43], TOPCONS [44] and MemBrain [45]. As a pre-processing step, we removed predicted N-terminal signal peptide sequences by SVMSignal. For all predicted TMHs, a correctly predicted TMH is defined as a one-to-one overlap with the observed TMH of PDBTM, and the minimum residue number of overlaps between predicted TMHs and observed TMHs of PDBTM is eight. Within 188 TMHs of independent test set, the recall and precision of TOPCONS is 93.09% (175/188) and 100% (188/188), respectively. The recall and precision of MemBrain is 97.87% (184/188) and 95.34% (184/193), respectively. For fair comparison in rotational angles, we only discuss 175 TMHs that have been predicted well by both of the two topology predictors. Prediction results from SVMSignal, TOPCONS and MemBrain are available in the Additional file 3: Dataset S2.

Additional file 3: Dataset S2. Signal peptide and topology prediction results of the independent test set from SVMSignal, TOPCONS and MemBrain.

Format: ZIP Size: 5.8MB Download fileOpen Data

To obtain observed rotational angles corresponding to each predicted TMH, we removed atoms which were predicted outside the membrane by the topology prediction for each protein chain, and then we followed the definition in this work to calculate the rotational angles. Specifically, we did not directly assign rotational angles calculated by TMH of PDBTM to predicted TMHs, but we recalculated rotational angles based on atoms of protein structure selected by predicted topology within the TM region. There are two reasons to do that. First, since the sequence of a predicted TMH is not identical to that defined by PDBTM and the definition of rotational angle depends on the helical principal axis and the Cα vector of the first residue to its lipid-facing direction vector, we cannot simply assume their structural property is similar. Second, for a predicted 3D protein structure, the TM region information comes from the topology predictor, and the rotational angle of each TMH is established on the atoms within predicted TM region, not from the PDBTM. Therefore, we have to recalculate the rotational angles of predicted TMH for comparison. Finally, we ignored any predicted TMH which includes residues that do not have structural data within PDB entity, and 155 TMHs were left for comparison.

The Additional file 4: Table S2 and the Additional file 5: Table S3 provide sequences of predicted TMHs corresponding to TMHs annotated in PDBTM, observed angle defined by residues of predicted topology, predicted angle, and moment lengths. For all 155 TMHs, the MAAE of TOPCONS is 43.04˚ and MAAE of MemBrain is 56.59˚. These two tables demonstrate the ability of TMexpo to predict rotational angles based on predicted topology. Interestingly, the MAAE of 155 TMHs based on topology predicted by TOPCONS has better results than topology annotated in PDBTM. There are two possible explanations for this observation. First, TMHs predicted by TOPCONS is longer than annotated in PDBTM, and this may help calculating the helical principle axis. Second, we excluded TMHs that have partially incomplete structural data, and the performance of 155 TMHs would differ from the dataset of 188TMHs. We conclude that rotational angles calculated by both predicted TMHs and predicted rASA are still consistent with the observed rotational angles defined by the predicted TMH. Therefore, while TMH boundary is not perfectly predicted, the predicted rotational angle can still provide useful information to the interior side of a TM protein and constrain decoys of predicted 3D structure.

Additional file 4: Table S2. Sequences, observed angles defined by residues of predicted topology (by TOPCONS), predicted angles, moment lengths and MAAE on the 155 TMHs of the independent test set.

Format: DOCX Size: 47KB Download fileOpen Data

Additional file 5: Table S3. Sequences, observed angles defined by residues of predicted topology (by MemBrain), predicted angles, moment lengths and MAAE on the 155 TMHs of the independent test set.

Format: DOCX Size: 47KB Download fileOpen Data

An application to an amino acid antiporter, AdiC

We selected from our independent test set an amino acid antiporter, called AdiC, of E. coli strain O157:H7 to demonstrate how predicted rASA and rotational angle in the helical wheel presentation of TMHs facilitate the analysis of TM proteins. The E. coli strain O157:H7 is a pathogen and causes hemorrhagic diarrhea, and AdiC is a multi-spanning TM protein that enables E. coli to resist the acidic environment via exchanging extracellular arginine and intracellular agmatine [46,47]. An arginine-bounded structure of AdiC was solved and is coded as 3L1L [46] in PDB.

The protein AdiC has 445 amino acids and contains 12 TMHs, including an inner layer containing 5 TMHs, i.e., transmembrane segments 1, 3, 6, 8 and 10 from the N-terminus, and an outer layer containing 7 TMHs, i.e., transmembrane segments 2, 4, 5, 7, 9, 11 and 12. Figure 1 illustrates the top view of 3L1L, in which each residue in a TMH is shown by a color-coded gradient scale of its observed rASA, with dark red and dark blue representing extremely buried and extremely exposed, respectively. The red arrow inside each helical wheel indicates the lipid-facing direction, and the green clock-wise arrow starts at the first amino acid. A gray line connects two TMHs if they have interhelical interactions with more than five VDWs as annotated in TMPad. It can be observed that most of residues in the five TMHs of the inner layer are colored as dark red, and the TMHs of the outer layer have residues in both exposed and buried sides colored by blue and red, respectively.

thumbnailFigure 1. Top view of the structure 3L1L shown by helical wheels, interhelical interactions and observed rASA of the TM residues. Each TM residue is shown by a color-coded gradient scale of its observed rASA, with dark red and dark blue representing rASA value 1.0 for extremely buried and 0 for extremely exposed, respectively. The red arrow inside each helical wheel indicates the lipid-facing direction, and the green clock-wise arrow starts at the first amino acid.

Interestingly, several binding sites and functional mutagenesis were discovered in TMHs of the inner layer. According to annotations in the UniProt/Swiss-Prot, the amino acids 22, 26, 93, 208 and 365 are substrate binding sites located in the TMHs 1, 3, 6 and 10 of the inner layer. Furthermore, four reported mutations, Y87A, Y93A, Y93K and Y365A, related to transporter activity also occur in TMHs 3 and 10. Therefore, TMHs of the inner layer are important to transporter function study. All of the above six residues are classified as buried in the crystallized structure 3L1L, and their observed rASA range from 0% to 2.13%. TMexpo predicted five of the above six residues as buried. Only one residue 87Y was predicted as marginally exposed and had 9.03% rASA. Figure 2 shows that the observed and predicted lipid-facing directions of the 12 TMHs are quite close. The MAAE of 12 TMHs is 36.72˚, and 7 out of 12 TMHs were accurately predicted with angular errors <30˚. The best angular error is 0.49˚ obtained by the linear TMH 7, and the worst angular error is 153.48˚ obtained by the kinked TMH 6. The PCC between the observed and predicted rASA of the entire 266 TM residues is 0.75; the MCC of the TMexpo’s classification of all the TM residues is 0.60. Our analysis shows that the predicted rASA of each residue could facilitate the design of site-directed mutagenesis experiments. Therefore, sequence-based prediction of both rotational angle and rASA can provide useful information in transporter function study, and thus contribute to experimental design to elucidate TM protein functions when high-resolution structure is unavailable.

thumbnailFigure 2. Observed and predicted rASA and lipid-facing directions of the 12 TMHs in 3L1L represented by helical wheels. The angular error of 12 TMHs ranges from 0.49˚ to 153.48˚; the MAAE of 12 TMHs is 36.72˚.

Interhelical contacts play an important role in relative accessible surface area prediction

To understand how intrachain interhelical interactions affect the burial status of residues is important for helical packing and stability of TM proteins [48-51]. At first we observed residues carrying more intrachain interhelical contacts with other TMHs tend to be buried. We investigated on the TMHs of the development set of 89 protein chains whether buried residues are hard to predict their rASA values. Using the annotations in TMPad, we separated residues into two sets according to different thresholds c of VDW contact numbers; one is called contact-enriched set of residues with at least c contacts, and the other set is called reference set consisting of the remaining residues. For each set, we calculated the PCC between the observed rASA and the predicted rASA obtained from LOOCV. The results are shown in Table 6, where different thresholds c = 1, c = 2, and c = 3 are considered. By comparing correlations of predicted rASA against observed rASA for different thresholds in LOOCV, we found that the set containing contact-enriched residues has lower correlation than the reference set by 40.66% when c = 3 and by 9.53% even when c = 1. This observation implies rASA of residues that have a large number of interhelical contacts may be more difficult to be predicted, and then the information of interhelical contacts should be integrated for improving rASA prediction.

Table 6. Comparing Pearson correlation coefficients between contact-enriched set and reference set defined by different thresholds ( c )

Conclusions

Sequence-based prediction of both rotational angle and rASA can provide indispensable information for structure prediction when high-resolution structures are unavailable and contribute to experimental design to elucidate TM protein functions. In this paper, we present a novel concept of using lipid exposure to infer rotational angles and have developed a machine learning approach to predict rotational angles of TMHs. Significantly, using predicted rASA from our sequence-based model achieved an MAAE of 48.31˚ on the independent test set, which is better than that obtained by the best of the compared knowledge-based propensities (67.51˚). Furthermore, we demonstrate an application for structural analysis via an amino acid antiporter. We believe improving prediction of rotational angle can benefit the structure prediction because free modeling of TM protein structures is a tough task and reducing the number of packing arrangements is necessary.

Methods

Evaluation measures

The metric used for evaluating rotational angle prediction in this work is mean absolute angular error (MAAE). To evaluate the classification model, i.e., classifying burial and exposed status, we used the following performance measures, including Matthew’s correlation coefficient (MCC), accuracy, sensitivity, specificity, and precision. With respect to the regression model, i.e., rASA prediction, we used mean absolute error (MAE), root mean squared error (RMSE), and Pearson correlation coefficient (PCC).

We evaluate the MAAE, which is the average angular error for all TMHs. Specifically, we defined observed rotational angle as x and predicted one as y for each helix, and both of them range from 0˚ to 360˚. The absolute angular error θerror, ranging from 0˚ to 180˚, is defined as the difference between two rotational angles of TMHs as Equation 1. The extreme error between observed and predicted angles occurs when they indicate opposite direction i.e. 180˚. The MAAE is defined as the average of all absolute angular error θerror within n helices as Equation 2.

<a onClick="popup('http://www.biomedcentral.com/1471-2105/14/304/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/14/304/mathml/M4">View MathML</a>

(1)

<a onClick="popup('http://www.biomedcentral.com/1471-2105/14/304/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/14/304/mathml/M5">View MathML</a>

(2)

Let predictedi and observedi denote the predicted and observed rASA, respectively, of the ith TMH in the sample dataset. MAE and RMSE are defined as Equations 3 and 4.

<a onClick="popup('http://www.biomedcentral.com/1471-2105/14/304/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/14/304/mathml/M6">View MathML</a>

(3)

<a onClick="popup('http://www.biomedcentral.com/1471-2105/14/304/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/14/304/mathml/M7">View MathML</a>

(4)

We defined exposed residues as positive data and buried residues as negative data, and consequently MCC are defined as Equation 5.

<a onClick="popup('http://www.biomedcentral.com/1471-2105/14/304/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/14/304/mathml/M8">View MathML</a>

(5)

See the Additional file 6: Table S4 for definitions of Pearson correlation coefficient, accuracy, sensitivity, specificity and precision.

Additional file 6: Table S4. Evaluation measures used in this work.

Format: DOCX Size: 25KB Download fileOpen Data

Data preparation

We retrieved the α-helical TM proteins from the PDBTM database [52,53], a collection of automatically identified TM proteins from the Protein Data Bank (PDB), and obtained 4,202 chains from 1,174 PDB entities. Among them, we only kept multi-spanning α-helical TM protein structures solved by X-rays with a resolution less than 4 Å, and then 2,305 protein chains remained. Next, the 2,305 chains were reduced at mutual sequence identity of less than 30% using CD-HIT [54]. Finally, the representative 110 multi-spanning protein chains listed in Table 7 were divided into two datasets for model development and independent test. We took 89 (554 TMHs) from 110 TM protein chains solved before October 20, 2008 as the development set and optimized the parameters in the models by leave-one-out cross-validation (LOOCV). The other 21 (188 TM helices) protein chains solved after October 20, 2008 were used for independent test. The two datasets are provided on the TMexpo server website (http://bio-cluster.iis.sinica.edu.tw/TMexpo webcite) and also in the Additional file 7: Dataset S3.

Table 7. The list of all protein chains (PDB:Chain) included in the development set and the independent test set

Additional file 7: Dataset S3. The representative 110 multi-spanning protein chains used in this work.

Format: ZIP Size: 1004KB Download fileOpen Data

Calculation of relative accessible surface area from structures

To calculate lipid exposure or exposed area of a structure, we used NACCESS program [19,29] with the probe radius set to 2.0 Å. The size of probe radius 2.0 Å was selected to mimic the -CH2 of hydrocarbon chains, and it is identical to that used in Yuan et al. [20], Illergård et al. [21] and Lo et al. [37]. The ASA for a residue was the sum of ASA from all atoms belonging to that residue. To extract the helical boundaries from the protein chains, we used the annotations of PDBTM. From the 89 protein chains in the development set, we obtained a total of 10,441 residues in TM domain. For the independent test set of newly solved proteins, we obtained 3,581 residues in TM domain. To annotate missing residues and missing atoms, we used PDB Validation Suite [55]. In order to obtain rASA as a normalized measure for a TM residue, we divided the ASA values by their reference values in a Gly-X-Gly tripeptide in an extended conformation. The reference values were derived from Samantha et al. [56]. To classify burial status of each residue for model training and testing, we followed the rASA threshold defined in Miller et al.’s work [57], i.e., rASA <5% to characterize buried residues and otherwise exposed, though different thresholds have been used in the literature.

An SVM-based predictor for lipid exposure of TM helices

Model development

We proposed residue-wise predictors based on support vector machines (SVMs), i.e., an SVM classifier to predict the burial/exposed status and a support vector regression (SVR) model to predict rASA values of each residue in TM domain. Specifically, C-SVC and epsilon-SVR implemented in LIBSVM [58] were used to develop the models, and both of them used the RBF kernel function. The parameters of the models were optimized by chain-wise LOOCV procedure on the development set. In LOOCV procedure, the best set of parameters to train the burial/exposed status classification model is of cost c = 21 and gamma g = 2-4; and the best set of parameters to train the real-number rASA regression model is of cost c = 2-1, gamma g = 2-5, loss function p = 10-3 and tolerance of termination criterion e = 10-2. Details of LOOCV performances can be obtained in the Additional file 8: Text S1.

Additional file 8: Text S1. Details of LOOCV performances during model development.

Format: PDF Size: 29KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Given a TM domain of a protein chain, each residue to be predicted was located at the center of a sliding window of length 17 and features were generated according to the 17-mer sequence. To train the classification model, exposed residues with label “E” were considered as positive data, and buried residues with label “B” as negative data. To train the regression model, the input was taken from the real-number rASA. We searched parameters by LOOCV procedure for the classification model and the regression model based on optimizing the MCC and the PCC, respectively. We did not directly predict ASA values because they are not normalized in a zero to one interval and this could produce bias in the presence of an outlier.

In training and testing, we excluded residues that participate in interchain contacts and the rationale is as follows: A sequence-based rASA predictor, which accepts the sequence of a structural subunit as input, can only describe structural properties of one subunit, not of the complete structure. Thus predicted rASA of those residues may be drastically different depending on their locations in the interacting interfaces. In the case of residues residing on the interchain surface, we observed rASA of these residues in a single chain may be significantly different from those seen in the complete structure with multiple subunits. Out of 110 representative protein chains used in our work, 86 protein chains are multimeric. Among all of their 9,800 TM residues without any missing atom, 2,167 (22.11%) residues have two different rASA values calculated from the single subunit and the complete structure, respectively; and the former-derived rASA is always larger than the latter-derived rASA. Notably, the maximum and average differences of the two rASA values of these 2,167 residues are 82.28% and 23.43%, respectively. Furthermore, 831 out of the 2,167 residues would be assigned inconsistent burial/exposed status according to their two different rASA values. In other words, 8.48% of the overall 9,800 transmembrane residues were considered as exposed from the perspective of a single subunit but turned out to be buried in their complete structures. For example, in 2OAR:A, 41 residues of 52 TM residues have different rASA values, and 13 residues are calculated as being exposed in the single chain but as being buried in the complete structure. It is noteworthy that the 60S (i.e., 37S by PDB indexing) and the 45 V (i.e., 22 V by PDB indexing) have drastic differences in their rASA values, i.e., 60.77% vs. 2.42% and 60.02% vs. 1.68%, as calculated by single chain and by complete structure, respectively. Since we did not know the native state of amino acids lying on the interchain surface, we excluded these residues from our training and testing data. For each protein chain, we calculated rASA for both single subunit structure and complete structure. Later, we excluded residues which were not identical in rASA by comparing the above two calculations, and also excluded residues that were missing partial or entire atoms.

In the testing stage, we performed a simple post-processing by rounding off their upper and lower bound to 1 and 0 because rASA values are contained in this interval. To derive the ASA values for each residue, we multiplied the predicted rASA values by the reference values [56].

Input features for predictors

In the design of TMexpo, we did not use a specific feature selection technique, and all the features used in TMexpo belong to one of the three feature groups. The first group is about interhelical contacts, specifically volume, polarity, charge and residue interhelical contact propensity. Since we have observed buried residues tend to have more interhelical contacts, and therefore we examined features related to interhelical contacts. For example, the well-known GxxxG motif can be regarded as small-xxx-small motifs [39], and we use volume profiles to incorporate such feature in the machine learning model. The polarity and charge can also be seen as features related to hydrogen bonding [50] and cation-pi interaction [59], respectively. To encode features into TMexpo, the volume [60] of each residue was divided by their maximum value 237.2 of tyrosine. The polarity was also encoded by the sigmoidal functions 1–1/(1 + e-po), where po denotes the mean residue polarity calculated by Radzicka and Wolfenden’s method [61]. We defined positively charged residues as 1, neutral residues as 0.5, and negatively charged residues as 0 based on the index used by Klein et al. [62]. The residue interhelical contact propensity were developed by Lo et al. [38], and we used in TMexpo the normalized propensity by division of the maximum value 1.43 of cysteine.

The second group provides evolutionary information as position-specific scoring matrix (PSSM) profiles and conservation score to machine learning model. Evolutionary information is an important feature and has been incorporated in interhelical interaction predictors [38,63]. To encode PSSM as features, the matrix was generated by performing PSI-BLAST against NCBI’s non-redundant database. This feature of a 17-mer peptide was encoded by a vector of size 17 × 20, where each entry was normalized by 1- 1/(1 + e-PSSM). The conservation score was calculated by an algorithm developed by Capra and Singh [64] on the multiple sequence alignment generated by MAFFT [65,66] based on the 17-mer peptide. We used the raw scores without using the local Z-score transformation described in their method.

The third group includes the TMH insertion energy, amphiphilicity of residues and turn propensities, which relate to structural information. The first two features can reveal residue position toward hydrophobic membrane or water interface, which is akin to Zpred features used in MPRAP that directly predict relative position from the center of membrane for each residue. The position-specific free energy of TMH insertion, termed as “free energy” to describe the hydrophobic core, was encoded by a sigmoidal function as 1-1/(1 + e-energy), where energy denotes the free energy of TMH insertion estimated by Hessa et al.’s method [67]. The amphiphilicity was encoded by the sigmoidal functions 1–1/(1 + eam), where am denotes the amphiphilicity derived by Mitaku et al.’s method [68]. We also considered the helix turn propensities in order to capture sequence information related to tight turns in naturally occurring TM helices from Monné et al. [69]. This feature was normalized propensities to [0, 1] by dividing the maximum value 2.7 of proline.

All of the above features were normalized to a closed [0, 1] interval. A feature value close to 1 means the corresponding residue is more hydrophobic, more amphiphilic, higher polarity, positively charged, larger volume, more conserved, tends to have turns and interhelical contacts. We filled 0.5 as features for nonexistent residues in windows, except charge, interhelical contacts, and volume, we filled zeroes.

Predicting rotational angle based on relative accessible surface area

Determination of rotational angle of a transmembrane helix

The rotational angle of a TMH was calculated as follows: First, we removed atoms which were annotated outside the membrane. Second, we computed the helical principal axis of the TMH of interest and aligned it with the z-axis with the N-terminal facing the screen, creating a top-view of the protein with respect to the target helix. Third, we identified the geometric centers of the molecule and each of the individual helices in the two-dimensional plane from the average x and y coordinates of Cα in the constituent TMH residues. We defined the lipid-facing direction of a TMH as in the opposite direction circumscribed by the geometric center of the target TMH connects to the molecular geometric center of the protein chain unit. The rotational angle of the target helix was measured as the angle rotated from the Cα vector of the first residue to its lipid-facing direction vector by clockwise motion viewed from the helix N-terminal to C-terminal. The angle ranges from 0˚ to 360˚.

Calculation of relative accessible surface area moment direction

The moment M was computed for several propensity scales with the moment length |M| and the moment direction θ relative to the angular direction of the Cα vector of the first residue of the TMH [9]. In this work, we used the rASA moment direction to describe lipid-facing direction and used the definition and post-processing techniques for both outliers and ramps described in Donnelly et al. [70] to resolve the moment direction θ. They were defined as follows:

<a onClick="popup('http://www.biomedcentral.com/1471-2105/14/304/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/14/304/mathml/M9">View MathML</a>

(6)

<a onClick="popup('http://www.biomedcentral.com/1471-2105/14/304/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/14/304/mathml/M10">View MathML</a>

(7)

<a onClick="popup('http://www.biomedcentral.com/1471-2105/14/304/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/14/304/mathml/M11">View MathML</a>

(8)

<a onClick="popup('http://www.biomedcentral.com/1471-2105/14/304/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/14/304/mathml/M12">View MathML</a>

(9)

For each residue, the rASA values can be seen as degree of directional lipid-facing. Therefore, for one TMH, the summation of all TMH residues’ lipid-facing tendency can characterize its rotational angle. In Equation 6, the x and the y terms are the vector summation of over n residues in a TMH. The moment length |M| is defined as Equation 7. The angle γ was solved first by inverse cosine function as Equation 8, and we determined moment direction θ by taking (360 - γ) as a result as Equation 9 if the sign of y term is negative.

Predicting rotational angles of transmembrane helices

In a previous study, propensity values, such as hydrophobicity or lipid-facing propensities, were used to describe the moment direction [9,11,12,14]. The procedure of calculating rotational angles of TMHs from predicted rASA values is described as follows: First, we predict rASA values from sequence by TMexpo. Second, we select each TMH according to TM domain annotated by PDBTM. To predict rotational angle, we further assign predicted rASA values of each TMH using Equation 6 and Equation 7 to obtain predicted rASA moment length (<a onClick="popup('http://www.biomedcentral.com/1471-2105/14/304/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/14/304/mathml/M13">View MathML</a>). Finally, we calculate the angles of moments for each TMH by Equation 8 and Equation 9, then take these angles as predicted rotational angles as lipid-facing directions of TMHs. Figure 3 illustrates the workflow of TMexpo for predicting rotational angles from sequences. We followed the above procedure for several propensity scales or rASA to obtain and compare the predicted rotational angles.

thumbnailFigure 3. The workflow of TMexpo for predicting rotational angles from sequences.

Abbreviations

ASA: Accessible surface area; LOOCV: Leave-one-out cross validation; MCC: Matthew’s correlation coefficient; MAAE: Mean absolute angular error; MAE: Mean absolute error; PCC: Pearson correlation coefficient; PDB: Protein data bank; rASA: Relative accessible surface area; RMSE: Root mean squared error; SVM: Support vector machine; SVR: Support vector regression; TM: Transmembrane; TMH: Transmembrane helix; VDW: Van der Waal’s.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

JSL, CCW, AL and TYS conceived and designed the experiments. JSL and CCW performed the experiments. JSL analyzed the data. JSL and CCW developed analysis tools and web services. JSL, CCW, AL, TYS and WLH wrote the manuscript. All authors have read and approved the manuscript.

Acknowledgements

This work was supported in part by the National Science Council under grants NSC100-2319-B-010-002 and NSC101-2221-E-001-022. No additional external funding received for this study. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

References

  1. Wallin E, Von Heijne G: Genome-wide analysis of integral membrane proteins from eubacterial, archaean, and eukaryotic organisms.

    Protein Sci 1998, 7:1029-1038. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  2. Stevens TJ, Arkin IT: Do more complex organisms have a greater proportion of membrane proteins in their genomes?

    Proteins: Struct, Funct, Bioinform 2000, 39:417-420. Publisher Full Text OpenURL

  3. Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE: The protein data bank.

    Nucleic Acids Res 2000, 28:235-242. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  4. Feng L, Campbell EB, Hsiung Y, MacKinnon R: Structure of a eukaryotic CLC transporter defines an intermediate state in the transport cycle.

    Science 2010, 330:635-641. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  5. Cao Y, Jin X, Huang H, Derebe MG, Levin EJ, Kabaleeswaran V, Pan Y, Punta M, Love J, Weng J: Crystal structure of a potassium ion transporter, TrkH.

    Nature 2011, 471:336-340. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. Von Heijne G: Membrane-protein topology.

    Nat Rev Mol Cell Biol 2006, 7:909-918. PubMed Abstract | Publisher Full Text OpenURL

  7. Popot JL, Engelman DM: Membrane protein folding and oligomerization: the two-stage model.

    Biochemistry 1990, 29:4031-4037. PubMed Abstract | Publisher Full Text OpenURL

  8. Popot J-L, Engelman DM: Helical membrane protein folding, stability, and evolution.

    Annu Rev Biochem 2000, 69:881-922. PubMed Abstract | Publisher Full Text OpenURL

  9. Eisenberg D, Weiss RM, Terwilliger TC: The helical hydrophobic moment: a measure of the amphiphilicity of a helix.

    Nature 1982, 299:371-374. PubMed Abstract | Publisher Full Text OpenURL

  10. Eisenberg D, Weiss RM, Terwilliger TC, Wilcox W: Hydrophobic moments and protein structure.

    Faraday Symp Chem Soc 1982, 17:109-120. OpenURL

  11. Beuming T, Weinstein H: A knowledge-based scale for the analysis and prediction of buried and exposed faces of transmembrane domain proteins.

    Bioinformatics 2004, 20:1822-1835. PubMed Abstract | Publisher Full Text OpenURL

  12. Pilpel Y, Ben-Tal N, Lancet D: kPROT: a knowledge-based scale for the propensity of residue orientation in transmembrane segments, application to membrane protein structure prediction.

    J Mol Biol 1999, 294:921-935. PubMed Abstract | Publisher Full Text OpenURL

  13. Adamian L, Nanda V, DeGrado WF, Liang J: Empirical lipid propensities of amino acid residues in multispan alpha helical membrane proteins.

    Proteins: Struct, Funct, Bioinformatics 2005, 59:496-509. Publisher Full Text OpenURL

  14. Park Y, Helms V: On the derivation of propensity scales for predicting exposed transmembrane residues of helical membrane proteins.

    Bioinformatics 2007, 23:701-708. PubMed Abstract | Publisher Full Text OpenURL

  15. Dastmalchi S, Beheshti S, Morris MB, Bret Church W: Prediction of rotational orientation of transmembrane helical segments of integral membrane proteins using new environment-based propensities for amino acids derived from structural analyses.

    FEBS J 2007, 274:2653-2660. PubMed Abstract | Publisher Full Text OpenURL

  16. Stevens TJ, Arkin IT: Substitution rates in α-helical transmembrane proteins.

    Protein Sci 2001, 10:2507-2517. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Adamian L, Liang J: Prediction of transmembrane helix orientation in polytopic membrane proteins.

    BMC Struct Biol 2006, 6:13. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  18. Hildebrand PW, Lorenzen S, Goede A, Preissner R: Analysis and prediction of helix-helix interactions in membrane channels and transporters.

    Proteins: Struct, Funct, Bioinformatics 2006, 64:253-262. Publisher Full Text OpenURL

  19. Lee B, Richards FM: The interpretation of protein structures: estimation of static accessibility.

    J Mol Biol 1971, 55:379-400. PubMed Abstract | Publisher Full Text OpenURL

  20. Yuan Z, Zhang F, Davis MJ, Bodén M, Teasdale RD: Predicting the solvent accessibility of transmembrane residues from protein sequence.

    J Proteome Res 2006, 5:1063-1070. PubMed Abstract | Publisher Full Text OpenURL

  21. Illergård K, Callegari S, Elofsson A: MPRAP: an accessibility predictor for a-helical transmem-brane proteins that performs well inside and outside the membrane.

    BMC Bioinforma 2010, 11:333. BioMed Central Full Text OpenURL

  22. Granseth E, Viklund H, Elofsson A: ZPRED: predicting the distance to the membrane center for residues in α-helical membrane proteins.

    Bioinformatics 2006, 22:e191-e196. PubMed Abstract | Publisher Full Text OpenURL

  23. Park Y, Hayat S, Helms V: Prediction of the burial status of transmembrane residues of helical membrane proteins.

    BMC Bioinforma 2007, 8:302. BioMed Central Full Text OpenURL

  24. Rose A, Lorenzen S, Goede A, Gruening B, Hildebrand PW: RHYTHM—a server to predict the orientation of transmembrane helices in channels and membrane-coils.

    Nucleic Acids Res 2009, 37:W575-W580. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. White SH, Von Heijne G: How translocons select transmembrane helices.

    Annu Rev Biophys 2008, 37:23-42. PubMed Abstract | Publisher Full Text OpenURL

  26. Driessen AJ, Manting EH, van der Does C: The structural basis of protein targeting and translocation in bacteria.

    Nat Struct Mol Biol 2001, 8:492-498. Publisher Full Text OpenURL

  27. Bibi E: The role of the ribosome-translocon complex in translation and assembly of polytopic membrane proteins.

    Trends Biochem Sci 1998, 23:51-55. PubMed Abstract | Publisher Full Text OpenURL

  28. White SH, Wimley WC: Membrane protein folding and stability: physical principles.

    Annu Rev Biophys Biomol Struct 1999, 28:319-365. PubMed Abstract | Publisher Full Text OpenURL

  29. Hubbard SJ, Thornton JM: ‘NACCESS’ Computer program, Department of Biochemistry and Molecular Biology. University College London; 1993. OpenURL

  30. Eisenberg D, Schwarz E, Komaromy M, Wall R: Analysis of membrane and surface protein sequences with the hydrophobic moment plot.

    J Mol Biol 1984, 179:125-142. PubMed Abstract | Publisher Full Text OpenURL

  31. Rees DC, DeAntonio L, Eisenberg D: Hydrophobic organization of membrane proteins.

    Science 1989, 245:510-513. PubMed Abstract | Publisher Full Text OpenURL

  32. Stevens TJ, Arkin IT: Are membrane proteins “inside-out” proteins?

    Proteins: Struct, Funct, Bioinformatics 1999, 36:135-143. Publisher Full Text OpenURL

  33. Mokrab Y, Stevens TJ, Mizuguchi K: Lipophobicity and the residue environments of the transmembrane α-helical bundle.

    Proteins: Struct, Funct, Bioinformatics 2009, 74:32-49. Publisher Full Text OpenURL

  34. Park Y, Helms V: How strongly do sequence conservation patterns and empirical scales correlate with exposure patterns of transmembrane helices of membrane proteins?

    Biopolymers 2006, 83:389-399. PubMed Abstract | Publisher Full Text OpenURL

  35. Hill JR, Kelm S, Shi J, Deane CM: Environment specific substitution tables improve membrane protein alignment.

    Bioinformatics 2011, 27:i15-i23. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  36. Ulmschneider MB, Sansom MSP: Amino acid distributions in integral membrane protein structures,Biochimica et Biophysica Acta (BBA).

    Biomembranes 2001, 1512:1-14. Publisher Full Text OpenURL

  37. Lo A, Cheng CW, Chiu YY, Sung TY, Hsu WL: TMPad: an integrated structural database for helix-packing folds in transmembrane proteins.

    Nucleic Acids Res 2011, 39:D347-D355. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  38. Lo A, Chiu YY, Rødland EA, Lyu PC, Sung TY, Hsu WL: Predicting helix-helix interactions from residue contacts in membrane proteins.

    Bioinformatics 2009, 25:996-1003. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  39. Harrington SE, Ben-Tal N: Structural determinants of transmembrane helical proteins.

    Structure 2009, 17:1092-1103. PubMed Abstract | Publisher Full Text OpenURL

  40. Jasti J, Furukawa H, Gonzales EB, Gouaux E: Structure of acid-sensing ion channel 1 at 1.9 Å resolution and low pH.

    Nature 2007, 449:316-323. PubMed Abstract | Publisher Full Text OpenURL

  41. Huang L-S, Shen JT, Wang AC, Berry EA: Crystallographic studies of the binding of ligands to the dicarboxylate site of complex II, and the identity of the ligand in the “oxaloacetate-inhibited” state, Biochimica et Biophysica Acta (BBA).

    Bioenergetics 2006, 1757:1073-1083. Publisher Full Text OpenURL

  42. Gonen T, Cheng Y, Sliz P, Hiroaki Y, Fujiyoshi Y, Harrison SC, Walz T: Lipid-protein interactions in double-layered two-dimensional AQP0 crystals.

    Nature 2005, 438:633-638. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  43. Lai JS, Cheng CW, Sung TY, Hsu WL: Computational comparative study of tuberculosis proteomes using a model learned from signal peptide structures.

    Plos One 2012, 7:e35018. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  44. Bernsel A, Viklund H, Hennerdal A, Elofsson A: TOPCONS: consensus prediction of membrane protein topology.

    Nucleic Acids Res 2009, 37:W465-W468. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  45. Shen H, Chou JJ: MemBrain: improving the accuracy of predicting transmembrane helices.

    Plos One 2008, 3:e2399. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  46. Gao X, Zhou L, Jiao X, Lu F, Yan C, Zeng X, Wang J, Shi Y: Mechanism of substrate recognition and transport by an amino acid antiporter.

    Nature 2010, 463:828-832. PubMed Abstract | Publisher Full Text OpenURL

  47. Gao X, Lu F, Zhou L, Dang S, Sun L, Li X, Wang J, Shi Y: Structure and mechanism of an amino acid antiporter.

    Science 2009, 324:1565-1568. PubMed Abstract | Publisher Full Text OpenURL

  48. Faham S, Yang D, Bare E, Yohannan S, Whitelegge JP, Bowie JU: Side-chain contributions to membrane protein structure and stability.

    J Mol Biol 2004, 335:297-305. PubMed Abstract | Publisher Full Text OpenURL

  49. Fleming KG, Engelman DM: Specificity in transmembrane helix-helix interactions can define a hierarchy of stability for sequence variants.

    Proc Natl Acad Sci U S A 2001, 98:14340-14344. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  50. Bowie JU: Membrane protein folding: how important are hydrogen bonds?

    Curr Opin Struct Biol 2011, 21:42-49. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  51. Joh NH, Min A, Faham S, Whitelegge JP, Yang D, Woods VL, Bowie JU: Modest stabilization by most hydrogen-bonded side-chain interactions in membrane proteins.

    Nature 2008, 453:1266-1270. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  52. Tusnády GE, Dosztányi Z, Simon I: PDB_TM: selection and membrane localization of transmembrane proteins in the protein data bank.

    Nucleic Acids Res 2005, 33:D275-D278. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  53. Tusnády GE, Dosztányi Z, Simon I: Transmembrane proteins in the protein data bank: identification and classification.

    Bioinformatics 2004, 20:2964-2972. PubMed Abstract | Publisher Full Text OpenURL

  54. Li W, Godzik A: Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences.

    Bioinformatics 2006, 22:1658-1659. PubMed Abstract | Publisher Full Text OpenURL

  55. Westbrook J, Feng Z, Burkhardt K, Berman HM: Validation of protein structures for protein data bank.

    Methods Enzymol 2003, 374:370-385. PubMed Abstract | Publisher Full Text OpenURL

  56. Samanta U, Bahadur RP, Chakrabarti P: Quantifying the accessible surface area of protein residues in their local environment.

    Protein Eng 2002, 15:659-667. PubMed Abstract | Publisher Full Text OpenURL

  57. Miller S, Janin J, Lesk AM, Chothia C: Interior and surface of monomeric proteins.

    J Mol Biol 1987, 196:641-656. PubMed Abstract | Publisher Full Text OpenURL

  58. Chang C-C, Lin C-J: LIBSVM: a library for support vector machines.

    ACM Trans Intell Syst Technol 2011, 2:1-27. OpenURL

  59. Johnson RM, Hecht K, Deber CM: Aromatic and cation-π interactions enhance helix-helix association in a membrane environment.

    Biochemistry 2007, 46:9208-9214. PubMed Abstract | Publisher Full Text OpenURL

  60. Pontius J, Richelle J, Wodak SJ: Deviations from standard atomic volumes as a quality measure for protein crystal structures.

    J Mol Biol 1996, 264:121-136. PubMed Abstract | Publisher Full Text OpenURL

  61. Radzicka A, Wolfenden R: Comparing the polarities of the amino acids: side-chain distribution coefficients between the vapor phase, cyclohexane, 1-octanol, and neutral aqueous solution.

    Biochemistry 1988, 27:1664-1670. Publisher Full Text OpenURL

  62. Klein P, Kanehisa M, DeLisi C: Prediction of protein function from sequence properties: discriminant analysis of a data base.

    Biochim Biophys Acta Protein Struct Mol Enzymol 1984, 787:221-226. Publisher Full Text OpenURL

  63. Fuchs A, Martin-Galiano AJ, Kalman M, Fleishman S, Ben-Tal N, Frishman D: Co-evolving residues in membrane proteins.

    Bioinformatics 2007, 23:3312-3319. PubMed Abstract | Publisher Full Text OpenURL

  64. Capra JA, Singh M: Predicting functionally important residues from sequence conservation.

    Bioinformatics 2007, 23:1875-1882. PubMed Abstract | Publisher Full Text OpenURL

  65. Katoh K, Misawa K, Kuma KI, Miyata T: MAFFT: a novel method for rapid multiple sequence alignment based on fast fourier transform.

    Nucleic Acids Res 2002, 30:3059-3066. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  66. Katoh K, Kuma KI, Toh H, Miyata T: MAFFT version 5: improvement in accuracy of multiple sequence alignment.

    Nucleic Acids Res 2005, 33:511-518. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  67. Hessa T, Meindl-Beinker NM, Bernsel A, Kim H, Sato Y, Lerch-Bader M, Nilsson I, White SH, Von Heijne G: Molecular code for transmembrane-helix recognition by the Sec61 translocon.

    Nature 2007, 450:1026-1030. PubMed Abstract | Publisher Full Text OpenURL

  68. Mitaku S, Hirokawa T, Tsuji T: Amphiphilicity index of polar amino acids as an aid in the characterization of amino acid preference at membrane-water interfaces.

    Bioinformatics 2002, 18:608-616. PubMed Abstract | Publisher Full Text OpenURL

  69. Monné M, Nilsson I, Elofsson A, Von Heijne G: Turns in transmembrane helices: determination of the minimal length of a “helical hairpin” and derivation of a fine-grained turn propensity scale.

    J Mol Biol 1999, 293:807-814. PubMed Abstract | Publisher Full Text OpenURL

  70. Donnelly D, Overington JP, Ruffle SV, Nugent JHA, Blundell TL: Modeling α-helical transmembrane domains: the calculation and use of substitution tables for lipid-facing residues.

    Protein Sci 1993, 2:55-70. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL