Evaluation of methods for the prediction of membrane spanning regions

Mouse proteome analysis

A general overview of the protein sequence set for the mouse transcriptome produced during the FANTOM2 sequencing project is presented here. We applied different algorithms to characterize protein sequences derived from a nonredundant representative protein set (RPS) and a variant protein set (VPS) of the mouse transcriptome. The functional characterization and assignment of Gene Ontology terms was done by analysis of the proteome using InterPro. The Superfamily database analyses gave a detailed structural classification according to SCOP and provide additional evidence for the functional characterization of the proteome data. The MDS database analysis revealed new domains which are not presented in existing protein domain databases. Thus the transcriptome gives us a unique source of data for the detection of new functional groups. The data obtained for the RPS and VPS sets facilitated the comparison of different patterns of protein expression. A comparison of other existing mouse and human protein sequence sets (e.g., the International Protein Index) demonstrates the common patterns in mammalian proteomes. The analysis of the membrane organization within the transcriptome of multiple eukaryotes provides valuable statistics about the distribution of secretory and transmembrane proteins

The transmembrane domain and PXXP motifs of ApoE receptor 2 exclude it from carrying out clathrin-mediated endocytosis

The low density lipoprotein (LDL) receptor family comprises several proteins with similar structures including the LDL receptor and apoE receptor 2 (apoER2). The human brain expresses two major splice variants of apoER2 mRNA, one of which includes an additional exon that encodes 59 residues in the cytoplasmic domain. This exon is absent from the LDL receptor and contains three proline-rich (PXXP) motifs that may allow apoER2 to function as a signal transducer. To investigate the role of this insert, we took advantage of the well characterized low density lipoprotein receptor pathway. Chimeras comprising the ectodomain and transmembrane domain of the LDL receptor fused to the cytoplasmic domain of apoER2 lacking the PXXP-containing residues are able to mediate clathrin-dependent endocytosis of LDL as effectively as cells expressing the LDL receptor but not if the PXXP insert is present in the protein. Although expressed on the cell surface, the PXXP-containing chimeric receptor is excluded from clathrin vesicles as judged by its failure to co-localize with adaptor protein-2 possibly due to interaction with intracellular adaptors or scaffolding proteins. Chimeras with the transmembrane domain of apoER2, predicted to be longer than that of the LDL receptor by several residues, fail to mediate endocytosis of LDL or to co-localize with adaptor protein-2 regardless of the presence or absence of the PXXP insert. Thus features of apoER2 that distinguish it as a signaling receptor, rather than as an endocytosis receptor like the LDL receptor, reside in or near the transmembrane domain and in the proline-rich motifs.

Calcineurin/nuclear factors of activated T cells (NFAT)-activating and immunoreceptor tyrosine-based activation motif (ITAM)-containing protein (CNAIP), a novel ITAM-containing protein that activates the calcineurin/NFAT-signaling pathway

We report in this study the identification and characterization of a novel protein that we designated as calcineurin/NFAT-activating and immunoreceptor tyrosine-based activation motif (ITAM)-containing protein (CNAIP). The predicted 270-amino acid sequence contains an N-terminal signal peptide, an immunoglobin domain in the extracellular region, a transmembrane domain and an ITAM in the cytoplasmic tail. Quantitative reverse transcription-PCR showed that CNAIP was preferentially expressed in neutrophils, monocytes, mast cells, and other immune-related cells. Co-transfection of CNAIP expression constructs with luciferase reporter plasmids in HMC-1 cells resulted in activation of interleukin-13 and tumor necrosis factor-alpha promoters, which was mediated through the calcineurin/NFAT-signaling pathway. Mutation of either or both tyrosines in the ITAM abolished transcriptional activation induced by CNAIP, indicating that the ITAM is indispensable for CNAIP function in activating cytokine gene promoters. Thus, it is concluded that CNAIP is a novel ITAM-containing protein that activates the calcineurin/NFAT-signaling pathway and the downstream cytokine gene promoters.

An analysis of the Candida albicans genome database for soluble secreted proteins using computer-based prediction algorithms

We sought to identify all genes in the Candida albicans genome database whose deduced proteins would likely be soluble secreted proteins (the secretome). While certain C. albicans secretory proteins have been studied in detail, more data on the entire secretome is needed. One approach to rapidly predict the functions of an entire proteome is to utilize genomic database information and prediction algorithms. Thus, we used a set of prediction algorithms to computationally define a potential C. albicans secretome. We first assembled a validation set of 47 C. albicans proteins that are known to be secreted and 47 that are known not to be secreted. The presence or absence of an N-terminal signal peptide was correctly predicted by SignalP version 2.0 in 47 of 47 known secreted proteins and in 47 of 47 known non-secreted proteins. When all 6165 C. albicans ORFs from CandidaDB were analysed with SignalP, 495 ORFs were predicted to encode proteins with N-terminal signal peptides. In the set of 495 deduced proteins with N-terminal signal peptides, 350 were predicted to have no transmembrane domains (or a single transmembrane domain at the extreme N-terminus) and 300 of these were predicted not to be GPI-anchored. TargetP was used to eliminate proteins with mitochondrial targeting signals, and the final computationally-predicted C. albicans secretome was estimated to consist of up to 283 ORFs. The C. albicans secretome database is available at http://info.med.yale.edu/intmed/infdis/candida/

A method for the comprehensive proteomic analysis of membrane proteins

We describe a method that allows for the concurrent proteomic analysis of both membrane and soluble proteins from complex membrane-containing samples. When coupled with multidimensional protein identification technology (MudPIT), this method results in (i) the identification of soluble and membrane proteins, (ii) the identification of post-translational modification sites on soluble and membrane proteins, and (iii) the characterization of membrane protein topology and relative localization of soluble proteins. Overlapping peptides produced from digestion with the robust nonspecific protease proteinase K facilitates the identification of covalent modifications (phosphorylation and methylation). High-pH treatment disrupts sealed membrane compartments without solubilizing or denaturing the lipid bilayer to allow mapping of the soluble domains of integral membrane proteins. Furthermore, coupling protease protection strategies to this method permits characterization of the relative sidedness of the hydrophilic domains of membrane proteins.

Expression profiling and bioinformatic analyses of a novel stress-regulated multispanning transmembrane protein family from cereals and Arabidopsis

Cold acclimation is a multigenic trait that allows hardy plants to develop efficient tolerance mechanisms needed for winter survival. To determine the genetic nature of these mechanisms, several cold-responsive genes of unknown function were identified from cold-acclimated wheat (Triticum aestivum). To identify the putative functions and structural features of these new genes, integrated genomic approaches of data mining, expression profiling, and bioinformatic predictions were used. The analyses revealed that one of these genes is a member of a small family that encodes two distinct groups of multispanning transmembrane proteins. The cold-regulated (COR)413-plasma membrane and COR413-thylakoid membrane groups are potentially targeted to the plasma membrane and thylakoid membrane, respectively. Further sequence analysis of the two groups from different plant species revealed the presence of a highly conserved phosphorylation site and a glycosylphosphatidylinositol-anchoring site at the C-terminal end. No homologous sequences were found in other organisms suggesting that this family is specific to the plant kingdom. Intraspecies and interspecies comparative gene expression profiling shows that the expression of this gene family is correlated with the development of freezing tolerance in cereals and Arabidopsis. In addition, several members of the family are regulated by water stress, light, and abscisic acid. Structure predictions and comparative genome analyses allow us to propose that the cor413 genes encode putative G-protein-coupled receptors.

Membrane topology of PssT, the transmembrane protein component of the type I exopolysaccharide transport system in Rhizobium leguminosarum bv. trifolii strain TA1

The pssT gene was identified as the fourth gene located upstream of the pssNOP gene cluster possibly involved in the biosynthesis, polymerization, and transport of exopolysaccharide (EPS) in Rhizobium leguminosarum bv. trifolii strain TA1. The hydropathy profile and homology searches indicated that PssT belongs to the polysaccharide-specific transport family of proteins, a component of the type I system of the polysaccharide transport. The predicted membrane topology of the PssT protein was examined with a series of PssT-PhoA fusion proteins and a complementary set of PssT-LacZ fusions. The results generally support a predicted topological model for PssT consisting of 12 transmembrane segments, with amino and carboxyl termini located in the cytoplasm. A mutant lacking the C-terminal part of PssT produced increased amounts of total EPS with an altered distribution of high- and low-molecular-weight forms in comparison to the wild-type RtTA1 strain. The PssT mutant produced an increased number of nitrogen fixing nodules on clover.

etramps, a new Plasmodium falciparum gene family coding for developmentally regulated and highly charged membrane proteins located at the parasite-host cell interface

After invasion of erythrocytes, the human malaria parasite Plasmodium falciparum resides within a parasitophorous vacuole and develops from morphologically and metabolically distinct ring to trophozoite stages. During these developmental phases, major structural changes occur within the erythrocyte, but neither the molecular events governing this development nor the molecular composition of the parasitophorous vacuole membrane (PVM) is well known. Herein, we describe a new family of highly cationic proteins from P. falciparum termed early transcribed membrane proteins (ETRAMPs). Thirteen members were identified sharing a conserved structure, of which six were found only during ring stages as judged from Northern and Western analysis. Other members showed different stage-specific expression patterns. Furthermore, ETRAMPs were associated with the membrane fractions in Western blots, and colocalization and selective permeabilization studies demonstrated that ETRAMPs were located in the PVM. This was confirmed by immunoelectron microscopy where the PVM and tubovesicular extensions of the PVM were labeled. Early expressed ETRAMPs clearly defined separate PVM domains compared with the negatively charged integral PVM protein EXP-1, suggesting functionally different domains in the PVM with an oppositely charged surface coat. We also show that the dynamic change of ETRAMP composition in the PVM coincides with the morphological changes during development. The P. falciparum PVM is an important structure for parasite survival, and its analysis might provide better understanding of the requirements of intracellular parasites.

Reliability measures for membrane protein topology prediction algorithms

We have developed reliability scores for five widely used membrane protein topology prediction methods, and have applied them both on a test set of 92 bacterial plasma membrane proteins with experimentally determined topologies and on all predicted helix bundle membrane proteins in three fully sequenced genomes: Escherichia coli, Saccharomyces cerevisiae and Caenorhabditis elegans. We show that the reliability scores work well for the TMHMM and MEMSAT methods, and that they allow the probability that the predicted topology is correct to be estimated for any protein. We further show that the available test set is biased towards high-scoring proteins when compared to the genome-wide data sets, and provide estimates for the expected prediction accuracy of TMHMM across the three genomes. Finally, we show that the performance of TMHMM is considerably better when limited experimental information (such as the in/out location of a protein's C terminus) is available, and estimate that at least ten percentage points in overall accuracy in whole-genome predictions can be gained in this way.

Two overlapping divergent transcription units in the human genome: the FEN1/C11orf10 locus

Flap endonuclease 1 (FEN-1) is a nuclear enzyme involved in DNA metabolism, such as replication, repair, and recombination. Here, we report the comparative genomic organization of the chicken, mouse, and human FEN1 genes as well as the comparative organization of a small gene (C11orf10) located immediately upstream of the FEN1 gene in reverse orientation. Immunostaining revealed that the C11orf10 protein, unlike FEN-1, is located in the cytoplasm, suggesting that these two proteins do not form a physical complex. Importantly, in the human genome, the two mRNAs are overlapping (14 bp) in their 5' ends. Thus, the FEN1/C11orf10 locus is a new example of two overlapping, divergent transcription units in the human genome.

Distribution analysis of nonsynonymous polymorphisms within the G-protein-coupled receptor gene family

The G-protein-coupled receptor (GPCR) superfamily is one of the largest classes of proteins in mammalian genomes. GPCRs mediate diverse physiological functions and are the targets of >50% of all clinical drugs. The sequencing of the human genome and large-scale polymorphism discovery efforts have established an abundant source of single nucleotide polymorphisms (SNPs), particularly those that result in a change in the encoded amino acids (cSNPs), many are of which in GPCRs. Although the majority of these cSNPs are assumed not to be disease-causing (nDCs), experimental data on their functional impact are lacking. Here, we have computationally analyzed the distribution of 454 cSNPs within the GPCR gene family and have found that disease-causing cSNPs (DCs) are overrepresented, whereas nDCs are underrepresented or neutral in transmembrane and extracellular loop domains, respectively. This finding reflects the relative importance of these domains to GPCR function and implies different biological characteristics for the two sets of human polymorphisms.

Prediction of the prototype of the human Toll-like receptor gene family from the pufferfish, Fugu rubripes, genome

The insect Toll family of proteins and their mammalian counterparts seemingly shared one common ancestor and evolved independently. Here we demonstrated that the prototype of the mammalian-type (M-type) Toll family is shared by the fish and humans. According to the draft of the pufferfish Fugu genome project, the signature Toll-IL-1 receptor homology domain (TIR domain) has been conserved during evolution. FuguTLR2, 3, 5, 7, 8 and 9 members correspond structurally to respective mammalian TLRs. One Fugu TLR showed equally high amino acid identity to human TLR1, 6 and 10, and we named it FuguTLR1. Fugu rubripes has genes for TLR21 and 22, which are unique to fish. One possible interpretation of these findings is that TLR1, 2, 3, 4, 5, 7, 8, 9, 21 and 22 existed in the ancestral genome common to fish and mammals, and that TLR4 was lost in the fish lineage, while TLR21 and 22 were lost in the mammalian lineage. Strikingly, a solitary ascidian, Halocynthia roretzi, has only a few Toll-like proteins, which, like Caenorhabditis elegans Toll, represent primitive ones before the expansion of the Toll family. Therefore, the expansion of TLR genes should have occurred earlier than fish, but not C. intestinalis, separated evolutionarily from mammals. These results infer that the appearance of the M-type innate system was completed before or concomitant with the appearance of acquired immunity. We interpret the present data to mean that the differences of TLRs identified in this study between fishes and humans may be rather peripheral, partially due to selection pressure exerted by pathogens in distinct environments.

Sequence analysis reveals new membrane anchor of reaction centre-bound cytochromes possibly related to PufX

Most of the bacterial photosynthetic reaction centres known to date contain a cytochrome subunit with four covalently bound haem groups. In the case of Blastochloris viridis, this reaction centre subunit is anchored in the membrane by a lipid molecule covalently attached to the cysteine which forms the N-terminus of the mature protein after processing by a signal peptidase. We show that posttranslational N-terminal cleavage of the cytochrome subunit does not occur in the aerobic photosynthetic bacterium Roseobacter denitrificans. From sequence analysis of the resulting elongated N-terminus it follows that a transmembrane helix is anchoring the reaction centre-bound cytochrome in the membrane. Comparative sequence analysis strongly suggests that all cytochrome subunits lacking the lipid coupling cysteine share this structural feature. Comparison of the N-terminal segment of the cytochrome subunit of Roseobacter denitrificans with the sequences of the PufX proteins from Rhodobacter sphaeroides and Rhodobacter capsulatus suggests a phylogenetic relation.

Comprehensive analysis of transmembrane topologies in prokaryotic genomes

We analysed comprehensively transmembrane (TM) topologies of TM proteins of 50 selected prokaryotic genomes, by discriminating between TM and soluble proteins by using SOSUI, then detecting and removing signal peptides by applying 'DetecSig', and finally predicting TM topologies by employing 'ConPred'. Estimated fraction of TM proteins in proteome averaged over the 50 genomes is approximately 22%. About 13% of TM proteins were predicted to have a signal peptide, and the fraction of soluble proteins with signal peptide (secretory proteins) ranges from 8 to 18% for most majority of the genomes. The N(in)-type TM proteins with 2-, 4-, 6- and 12-tms (number of transmembrane segments) are predominant among multi-spanning TM proteins, and correspondingly, significantly higher fractions of N(out)-type TM proteins with 1-, 3-, 5- and 11-tms have a signal peptide. It is also found that the TM proteins with signal peptide tend to have a long N-tail loop. The averaged sequence length of TM proteins increases linearly with the increase of the number of TM segments, with the increasing rate of about 35 residues, suggesting a possibility that TM topologies might have been evolved by the 'internal gene duplication' mechanism. Datasets of TM topologies predicted in this study are available at http://bioinfo.si.hirosaki-u.ac.jp/ approximately TMPinGS/.

SeqVISTA: a graphical tool for sequence feature visualization and comparison

BACKGROUND

Many readers will sympathize with the following story. You are viewing a gene sequence in Entrez, and you want to find whether it contains a particular sequence motif. You reach for the browser's "find in page" button, but those darn spaces every 10 bp get in the way. And what if the motif is on the opposite strand? Subsequently, your favorite sequence analysis software informs you that there is an interesting feature at position 13982-14013. By painstakingly counting the 10 bp blocks, you are able to examine the sequence at this location. But now you want to see what other features have been annotated close by, and this information is buried several screenfuls higher up the web page.

RESULTS

SeqVISTA presents a holistic, graphical view of features annotated on nucleotide or protein sequences. This interactive tool highlights the residues in the sequence that correspond to features chosen by the user, and allows easy searching for sequence motifs or extraction of particular subsequences. SeqVISTA is able to display results from diverse sequence analysis tools in an integrated fashion, and aims to provide much-needed unity to the bioinformatics resources scattered around the Internet. Our viewer may be launched on a GenBank record by a single click of a button installed in the web browser.

CONCLUSION

SeqVISTA allows insights to be gained by viewing the totality of sequence annotations and predictions, which may be more revealing than the sum of their parts. SeqVISTA runs on any operating system with a Java 1.4 virtual machine. It is freely available to academic users at http://zlab.bu.edu/SeqVISTA.

G Protein-Coupled Receptors in Invertebrates: A State of the Art

G protein-coupled receptors (GPCRs) constitute one of the largest and most ancient superfamilies of membrane-spanning proteins. We focus on neuropeptide GPCRs, in particular on those of invertebrates. In general, such receptors mediate the responses of signaling molecules that constitute the highest hierarchical position in the regulation of physiological processes. Until recently, only a few of these receptors were identified in invertebrates. However, the availability of a plethora of genomic information has boosted the discovery of novel members in several invertebrate species, such as Drosophila, in which 18 neuropeptide GPCRs have been characterized. The finalization of genomic projects in other invertebrates will lead to a similar expansion of GPCR understanding. Many new insights regarding neuropeptide regulation have followed from the discovery of their cognate receptors. Furthermore, information on GPCR signaling is still fragmentary and the elucidation of these pathways in model insects such as Drosophila will lead to further insights in other species, including mammals. In this review we present the current status of what is known about invertebrate GPCRs, discuss some novel perceptions that follow from the identified members, and, finally, present some future prospects.

TMPDB: a database of experimentally-characterized transmembrane topologies

TMPDB is a database of experimentally-characterized transmembrane (TM) topologies. TMPDB release 6.2 contains a total of 302 TM protein sequences, in which 276 are alpha-helical sequences, 17 beta-stranded, and 9 alpha-helical sequences with short pore-forming helices buried in the membrane. The TM topologies in TMPDB were determined experimentally by means of X-ray crystallography, NMR, gene fusion technique, substituted cysteine accessibility method, N-linked glycosylation experiment and other biochemical methods. TMPDB would be useful as a test and/or training dataset in improving the proposed TM topology prediction methods or developing novel methods with higher performance, and as a guide for both the bioinformaticians and biologists to better understand TM proteins. TMPDB and its subsets are freely available at the following web site: http://bioinfo.si.hirosaki-u.ac.jp/~TMPDB/.

ARAMEMNON, a novel database for Arabidopsis integral membrane proteins

A specialized database (DB) for Arabidopsis membrane proteins, ARAMEMNON, was designed that facilitates the interpretation of gene and protein sequence data by integrating features that are presently only available from individual sources. Using several publicly available prediction programs, putative integral membrane proteins were identified among the approximately 25,500 proteins in the Arabidopsis genome DBs. By averaging the predictions from seven programs, approximately 6,500 proteins were classified as transmembrane (TM) candidate proteins. Some 1,800 of these contain at least four TM spans and are possibly linked to transport functions. The ARAMEMNON DB enables direct comparison of the predictions of seven different TM span computation programs and the predictions of subcellular localization by eight signal peptide recognition programs. A special function displays the proteins related to the query and dynamically generates a protein family structure. As a first set of proteins from other organisms, all of the approximately 700 putative membrane proteins were extracted from the genome of the cyanobacterium Synechocystis sp. and incorporated in the ARAMEMNON DB. The ARAMEMNON DB is accessible at the URL http://aramemnon.botanik.uni-koeln.de.

SadA, a novel adhesion receptor in Dictyostelium

Little is known about cell-substrate adhesion and how motile and adhesive forces work together in moving cells. The ability to rapidly screen a large number of insertional mutants prompted us to perform a genetic screen in Dictyostelium to isolate adhesion-deficient mutants. The resulting substrate adhesion-deficient (sad) mutants grew in plastic dishes without attaching to the substrate. The cells were often larger than their wild-type parents and displayed a rough surface with many apparent blebs. One of these mutants, sadA-, completely lacked substrate adhesion in growth medium. The sadA- mutant also showed slightly impaired cytokinesis, an aberrant F-actin organization, and a phagocytosis defect. Deletion of the sadA gene by homologous recombination recreated the original mutant phenotype. Expression of sadA-GFP in sadA-null cells restored the wild-type phenotype. In sadA-GFP-rescued mutant cells, sadA-GFP localized to the cell surface, appropriate for an adhesion molecule. SadA contains nine putative transmembrane domains and three conserved EGF-like repeats in a predicted extracellular domain. The EGF repeats are similar to corresponding regions in proteins known to be involved in adhesion, such as tenascins and integrins. Our data combined suggest that sadA is the first substrate adhesion receptor to be identified in Dictyostelium.

Reliability of transmembrane predictions in whole-genome data

Transmembrane prediction methods are generally benchmarked on a set of proteins with experimentally verified topology. We have investigated if the accuracy measured on such datasets can be expected in an unbiased genomic analysis, or if there is a bias towards 'easily predictable' proteins in the benchmark datasets. As a measurement of accuracy, the concordance of the results from five different prediction methods was used (TMHMM, PHD, HMMTOP, MEMSAT, and TOPPRED). The benchmark dataset showed significantly higher levels (up to five times) of agreement between different methods than in 10 tested genomes. We have also analyzed which programs are most prone to make mispredictions by measuring the frequency of one-out-of-five disagreeing predictions.

GapA and CrmA coexpression is essential for Mycoplasma gallisepticum cytadherence and virulence

It was previously demonstrated that avirulent Mycoplasma gallisepticum strain R(high) (passage 164) is lacking three proteins that are expressed in its virulent progenitor, strain R(low) (passage 15). These proteins were identified as the cytadhesin molecule GapA, the putative cytadhesin-related molecule CrmA, and a component of a high-affinity transporter system, HatA. Complementation of R(high) with wild-type gapA restored expression in the transformant (GT5) but did not restore the cytadherence phenotype and maintained avirulence in chickens. These results suggested that CrmA might play an essential role in the M. gallisepticum cytadherence process. CrmA is encoded by the second gene in the gapA operon and shares significant sequence homology to the ORF6 gene of Mycoplasma pneumoniae, which has been shown to play an accessory role in the cytadherence process. Complementation of R(high) with wild-type crmA resulted in the transformant (SDCA) that lacked the cytadherence and virulence phenotype comparable to that found in R(high) and GT5. In contrast, complementation of R(high) with the entire wild-type gapA operon resulted in the transformant (GCA1) that restored cytadherence to the level found in wild-type R(low). In vivo pathogenesis trials revealed that GCA1 had regained virulence, causing airsacculitis in chickens. These results demonstrate that both GapA and CrmA are required for M. gallisepticum cytadherence and pathogenesis.

Transmembrane helix predictions revisited

Methods that predict membrane helices have become increasingly useful in the context of analyzing entire proteomes, as well as in everyday sequence analysis. Here, we analyzed 27 advanced and simple methods in detail. To resolve contradictions in previous works and to reevaluate transmembrane helix prediction algorithms, we introduced an analysis that distinguished between performance on redundancy-reduced high- and low-resolution data sets, established thresholds for significant differences in performance, and implemented both per-segment and per-residue analysis of membrane helix predictions. Although some of the advanced methods performed better than others, we showed in a thorough bootstrapping experiment based on various measures of accuracy that no method performed consistently best. In contrast, most simple hydrophobicity scale-based methods were significantly less accurate than any advanced method as they overpredicted membrane helices and confused membrane helices with hydrophobic regions outside of membranes. In contrast, the advanced methods usually distinguished correctly between membrane-helical and other proteins. Nonetheless, few methods reliably distinguished between signal peptides and membrane helices. We could not verify a significant difference in performance between eukaryotic and prokaryotic proteins. Surprisingly, we found that proteins with more than five helices were predicted at a significantly lower accuracy than proteins with five or fewer. The important implication is that structurally unsolved multispanning membrane proteins, which are often important drug targets, will remain problematic for transmembrane helix prediction algorithms. Overall, by establishing a standardized methodology for transmembrane helix prediction evaluation, we have resolved differences among previous works and presented novel trends that may impact the analysis of entire proteomes.

Proteomics of Synechocystis sp. strain PCC 6803: identification of plasma membrane proteins

Cyanobacteria are unique prokaryotes since they in addition to outer and plasma membranes contain the photosynthetic membranes (thylakoids). The plasma membranes of Synechocystis 6803, which can be completely purified by density centrifugation and polymer two-phase partitioning, have been found to be more complex than previously anticipated, i.e. they appear to be essential for assembly of the two photosystems. A proteomic approach for the characterization of cyanobacterial plasma membranes using two-dimensional gel electrophoresis and mass spectrometry analysis revealed a total of 57 different membrane proteins of which 17 are integral membrane spanning proteins. Among the 40 peripheral proteins 20 are located on the periplasmic side of the membrane, while 20 are on the cytoplasmic side. Among the proteins identified are subunits of the two photosystems as well as Vipp1, which has been suggested to be involved in vesicular transport between plasma and thylakoid membranes and is thus relevant to the possibility that plasma membranes are the initial site for photosystem biogenesis. Four subunits of the Pilus complex responsible for cell motility were also identified as well as several subunits of the TolC and TonB transport systems. Several periplasmic and ATP-binding proteins of ATP-binding cassette transporters were also identified as were two subunits of the F(0) membrane part of the ATP synthase.

The involvement of a multicopper oxidase in iron uptake by the green algae Chlamydomonas reinhardtii

In the unicellular green algae Chlamydomonas reinhardtii, high-affinity uptake of iron (Fe) requires an Fe(3+)-chelate reductase and an Fe transporter. Neither of these proteins nor their corresponding genes have been isolated. We previously identified, by analysis of differentially expressed plasma membrane proteins, an approximately 150-kD protein whose synthesis was induced under conditions of Fe-deficient growth. Based on homology of internal peptide sequences to the multicopper oxidase hephaestin, this protein was proposed to be a ferroxidase. A nucleotide sequence to the full-length cDNA clone for this ferroxidase-like protein has been obtained. Analysis of the primary amino acid sequence revealed a putative transmembrane domain near the amino terminus of the protein and signature sequences for two multicopper oxidase I motifs and one multicopper oxidase II motif. The ferroxidase-like gene was transcribed under conditions of Fe deficiency. Consistent with the role of a copper (Cu)-containing protein in Fe homeostasis, growth of cells in Cu-depleted media eliminated high-affinity Fe uptake, and Cu-deficient cells that were grown in optimal Fe showed greatly reduced Fe accumulation compared with control, Cu-sufficient cells. Reapplication of Cu resulted in the recovery of Fe transport activity. Together, these results were consistent with the participation of a ferroxidase in high-affinity Fe uptake in C. reinhardtii.

Prediction of partial membrane protein topologies using a consensus approach

We have developed a method to reliably identify partial membrane protein topologies using the consensus of five topology prediction methods. When evaluated on a test set of experimentally characterized proteins, we find that approximately 90% of the partial consensus topologies are correctly predicted in membrane proteins from prokaryotic as well as eukaryotic organisms. Whole-genome analysis reveals that a reliable partial consensus topology can be predicted for approximately 70% of all membrane proteins in a typical bacterial genome and for approximately 55% of all membrane proteins in a typical eukaryotic genome. The average fraction of sequence length covered by a partial consensus topology is 44% for the prokaryotic proteins and 17% for the eukaryotic proteins in our test set, and similar numbers are found when the algorithm is applied to whole genomes. Reliably predicted partial topologies may simplify experimental determinations of membrane protein topology.

Long membrane helices and short loops predicted less accurately

Low-resolution experiments suggest that most membrane helices span over 17-25 residues and that most loops between two helices are longer than 15 residues. Both constraints have been used explicitly in the development of prediction methods. Here, we compared the largest possible sequence-unique data sets from high- and low-resolution experiments. For the high-resolution data, we found that only half of the helices fall into the expected length interval and that half of the loops were shorter than 10 residues. We compared the accuracy of detecting short loops and long helices for 28 advanced and simple prediction methods: All methods predicted short loops less accurately than longer ones. In particular, loops shorter than 7 residues appeared to be very difficult to detect by current methods. Similarly, all methods tended to be more accurate for longer than for shorter helices. However, helices with more than 32 residues were predicted less accurately than all other helices. Our findings may suggest particular strategies for improving predictions of membrane helices.

Identification and characterization of a conserved baculoviral structural protein ODVP-6E/ODV-E56 from Choristoneura fumiferana granulovirus

A gene that encodes a homologue to baculoviral ODVP-6E/ODV-E56, a baculoviral envelope-associated viral structural protein, has been identified and sequenced on the genome of Choristoneura fumiferana granulovirus (ChfuGV). The ChfuGV odvp-6e/odv-e56 gene was located on an 11-kb BamHI subgenomic fragment using different sets of degenerated primers, which were designed using the results of the protein sequencing of a major 39 kDa structural protein that is associated with the occlusion-derived virus (ODV). The gene has a 1062 nucleotide (nt) open-reading frame (ORF) that encodes a protein with 353 amino acids with a predicted molecular mass of 38.5 kDa. The amino acid sequence data that was derived from the nucleotide sequence in ChfuGV was compared to those of other baculoviruses. ChfuGV ODVP-6E/ODV-E56, along with other baculoviral ODVP-6E/ODV-E56 proteins, all contained two putative transmembrane domains at their C-terminus. Several putative N- and O-glycosylation, N-myristoylation, and phosphorylation sites were detected in the ChfuGV ODVP-6E/ODV-E56 protein. A similar pattern was detected when a hydrophobicity-plots comparison was performed on ChfuGV ODVP-6E/ODV-E56 with other baculoviral homologue proteins. At the nucleotide level, a late promoter motif (GTAAG) was located at -14 nt upstream to the start codon of the ChfuGV odvp-6e/odv-e56 gene. A slight variant of the polyadenylation signal, AATAAT, was detected at the position +10 nt that is downstream from the termination signal. A phylogenetic tree for baculoviral ODVP-6E/ODV-E56 was constructed using a maximum parsimony analysis. The phylogenetic estimation demonstrated that ChfuGV ODVP-6E/ODV-E56 is most closely related to those of Cydia pomonella granulovirus (CpGV) and Plutella xylostella granulovirus (PxGV).

The predicted candidates of Arabidopsis plastid inner envelope membrane proteins and their expression profiles

Plastid envelope proteins from the Arabidopsis nuclear genome were predicted using computational methods. Selection criteria were: first, to find proteins with NH(2)-terminal plastid-targeting peptides from all annotated open reading frames from Arabidopsis; second, to search for proteins with membrane-spanning domains among the predicted plastidial-targeted proteins; and third, to subtract known thylakoid membrane proteins. Five hundred forty-one proteins were selected as potential candidates of the Arabidopsis plastid inner envelope membrane proteins (AtPEM candidates). Only 34% (183) of the AtPEM candidates could be assigned to putative functions based on sequence similarity to proteins of known function (compared with the 69% function assignment of the total predicted proteins in the genome). Of the 183 candidates with assigned functions, 40% were classified in the category of "transport facilitation," indicating that this collection is highly enriched in membrane transporters. Information on the predicted proteins, tissue expression data from expressed sequence tags and microarrays, and publicly available T-DNA insertion lines were collected. The data set complements proteomic-based efforts in the increased detection of integral membrane proteins, low-abundance proteins, or those not expressed in tissues selected for proteomic analysis. Digital northern analysis of expressed sequence tags suggested that the transcript levels of most AtPEM candidates were relatively constant among different tissues in contrast to stroma and the thylakoid proteins. However, both digital northern and microarray analyses identified a number of AtPEM candidates with tissue-specific expression patterns.

Complete genome sequence and analyses of the subgenomic RNAs of sweet potato chlorotic stunt virus reveal several new features for the genus Crinivirus

The complete nucleotide sequences of genomic RNA1 (9,407 nucleotides [nt]) and RNA2 (8,223 nt) of Sweet potato chlorotic stunt virus (SPCSV; genus Crinivirus, family Closteroviridae) were determined, revealing that SPCSV possesses the second largest identified positive-strand single-stranded RNA genome among plant viruses after Citrus tristeza virus. RNA1 contains two overlapping open reading frames (ORFs) that encode the replication module, consisting of the putative papain-like cysteine proteinase, methyltransferase, helicase, and polymerase domains. RNA2 contains the Closteroviridae hallmark gene array represented by a heat shock protein homologue (Hsp70h), a protein of 50 to 60 kDa depending on the virus, the major coat protein, and a divergent copy of the coat protein. This grouping resembles the genome organization of Lettuce infectious yellows virus (LIYV), the only other crinivirus for which the whole genomic sequence is available. However, in striking contrast to LIYV, the two genomic RNAs of SPCSV contained nearly identical 208-nt-long 3' terminal sequences, and the ORF for a putative small hydrophobic protein present in LIYV RNA2 was found at a novel position in SPCSV RNA1. Furthermore, unlike any other plant or animal virus, SPCSV carried an ORF for a putative RNase III-like protein (ORF2 on RNA1). Several subgenomic RNAs (sgRNAs) were detected in SPCSV-infected plants, indicating that the sgRNAs formed from RNA1 accumulated earlier in infection than those of RNA2. The 5' ends of seven sgRNAs were cloned and sequenced by an approach that provided compelling evidence that the sgRNAs are capped in infected plants, a novel finding for members of the Closteroviridae.

The novel EPTP repeat defines a superfamily of proteins implicated in epileptic disorders

Recent studies suggest that mutations in the LGI1/Epitempin gene cause autosomal dominant lateral temporal epilepsy. This gene encodes a protein of unknown function, which we postulate is secreted. The LGI1 protein has leucine-rich repeats in the N-terminal sequence and a tandem repeat (which we named EPTP) in its C-terminal region. A redefinition of the C-terminal repeat and the application of sensitive sequence analysis methods enabled us to define a new superfamily of proteins carrying varying numbers of the novel EPTP repeats in combination with various extracellular domains. Genes encoding proteins of this family are located in genomic regions associated with epilepsy and other neurological disorders.

Identification of Xenorhabdus nematophila genes required for mutualistic colonization of Steinernema carpocapsae nematodes

One stage in the symbiotic interaction between the bacterium Xenorhabdus nematophila and its nematode host, Steinernema carpocapsae, involves the species-specific colonization of the nematode intestinal vesicle by the bacterium. To characterize the bacterial molecular determinants that are essential for vesicle colonization, we adapted and applied a signature-tagged mutagenesis (STM) screen to this system. We identified 15 out of 3000 transposon mutants of X. nematophila with at least a 15-fold reduction in average vesicle colonization. These 15 mutants harbour disruptions in nine separate loci. Three of these loci have predicted open reading frames (ORFs) with similarity to genes (rpoS, rpoE, lrp) encoding regulatory proteins; two have predicted ORFs with similarity to genes (aroA, serC) encoding amino acid biosynthetic enzymes; one, designated nilB (nematode intestine localization), has an ORF with similarity to a gene encoding a putative outer membrane protein (OmpU) in Neisseria; and three, nilA, nilC and nilD, have no apparent homologues in the public database. nilA, nilB and nilC are linked on a single 4 kb locus. nilB and nilC are > 104-fold reduced in their ability to colonize the nematode vesicle and are predicted to encode membrane-localized proteins. The nilD locus contains an extensive repeat region and several small putative ORFs. Other than reduced colonization, the nilB, nilC and nilD mutants did not display alterations in any other phenotype tested, suggesting a specific role for these genes in allowing X. nematophila to associate with the nematode host.

On filtering false positive transmembrane protein predictions

While helical transmembrane (TM) region prediction tools achieve high (>90%) success rates for real integral membrane proteins, they produce a considerable number of false positive hits in sequences of known nontransmembrane queries. We propose a modification of the dense alignment surface (DAS) method that achieves a substantial decrease in the false positive error rate. Essentially, a sequence that includes possible transmembrane regions is compared in a second step with TM segments in a sequence library of documented transmembrane proteins. If the performance of the query sequence against the library of documented TM segment-containing sequences in this test is lower than an empirical threshold, it is classified as a non-transmembrane protein. The probability of false positive prediction for trusted TM region hits is expressed in terms of E-values. The modified DAS method, the DAS-TMfilter algorithm, has an unchanged high sensitivity for TM segments ( approximately 95% detected in a learning set of 128 documented transmembrane proteins). At the same time, the selectivity measured over a non-redundant set of 526 soluble proteins with known 3D structure is approximately 99%, mainly because a large number of falsely predicted single membrane-pass proteins are eliminated by the DAS-TMfilter algorithm.

Predicting protein cellular localization using a domain projection method

We investigate the co-occurrence of domain families in eukaryotic proteins to predict protein cellular localization. Approximately half (300) of SMART domains form a "small-world network", linked by no more than seven degrees of separation. Projection of the domains onto two-dimensional space reveals three clusters that correspond to cellular compartments containing secreted, cytoplasmic, and nuclear proteins. The projection method takes into account the existence of "bridging" domains, that is, instances where two domains might not occur with each other but frequently co-occur with a third domain; in such circumstances the domains are neighbors in the projection. While the majority of domains are specific to a compartment ("locale"), and hence may be used to localize any protein that contains such a domain, a small subset of domains either are present in multiple locales or occur in transmembrane proteins. Comparison with previously annotated proteins shows that SMART domain data used with this approach can predict, with 92% accuracy, the localizations of 23% of eukaryotic proteins. The coverage and accuracy will increase with improvements in domain database coverage. This method is complementary to approaches that use amino-acid composition or identify sorting sequences; these methods may be combined to further enhance prediction accuracy.

TRPV3 is a calcium-permeable temperature-sensitive cation channel

Transient receptor potential (TRP) proteins are cation-selective channels that function in processes as diverse as sensation and vasoregulation. Mammalian TRP channels that are gated by heat and capsaicin (>43 degrees C; TRPV1 (ref. 1)), noxious heat (>52 degrees C; TRPV2 (ref. 2)), and cooling (< 22 degrees C; TRPM8 (refs 3, 4)) have been cloned; however, little is known about the molecular determinants of temperature sensing in the range between approximately 22 degrees C and 40 degrees C. Here we have identified a member of the vanilloid channel family, human TRPV3 (hTRPV3) that is expressed in skin, tongue, dorsal root ganglion, trigeminal ganglion, spinal cord and brain. Increasing temperature from 22 degrees C to 40 degrees C in mammalian cells transfected with hTRPV3 elevated intracellular calcium by activating a nonselective cationic conductance. As in published recordings from sensory neurons, the current was steeply dependent on temperature, sensitized with repeated heating, and displayed a marked hysteresis on heating and cooling. On the basis of these properties, we propose that hTRPV3 is thermosensitive in the physiological range of temperatures between TRPM8 and TRPV1.

Molecular cloning of rabbit organic cation transporter rbOCT2 and functional comparisons with rbOCT1

Multiple organic cation transporters (OCTs) are present in rabbit kidney and may play different functional roles. We cloned rabbit OCT2 (rbOCT2) and compared its function with that of rabbit OCT1 (rbOCT1). In transiently transfected COS-7 cells, rbOCT1 and rbOCT2 mediated uptake of [3H]tetraethylammonium (TEA) with K(t) values of 188 and 125 microM, respectively. n-Tetraalkylammonium compounds showed similar affinities for the two homologs, with IC50 values for inhibition of OCT1- and OCT2-mediated [3H]TEA transport, respectively, of 4,538 and 1,395 microM for tetramethylammonium, 88.5 and 3.9 microM for tetrapropylammonium, 13.9 and 5.3 microM for tetrabutylammonium, and 8.8 and 7.6 microM for tetrapentylammonium. However, the transporters had very different affinities for cimetidine (CIM): IC50 of 916 and 5.7 microM for rbOCT1 and rbOCT2, respectively. CIM inhibition of TEA uptake into single S2 segments of rabbit proximal tubule was used to estimate the contributions of OCT1 and OCT2 to basolateral organic cation uptake. The median IC50 for CIM inhibition of TEA uptake was 12.3 microM, suggesting that OCT2 is the major contributor to basolateral organic cation transport in the S2 segment of proximal tubule in rabbit kidney.

Identification of signal peptide peptidase, a presenilin-type aspartic protease

Signal peptide peptidase (SPP) catalyzes intramembrane proteolysis of some signal peptides after they have been cleaved from a preprotein. In humans, SPP activity is required to generate signal sequence-derived human lymphocyte antigen-E epitopes that are recognized by the immune system, and to process hepatitis C virus core protein. We have identified human SPP as a polytopic membrane protein with sequence motifs characteristic of the presenilin-type aspartic proteases. SPP and potential eukaryotic homologs may represent another family of aspartic proteases that promote intramembrane proteolysis to release biologically important peptides.

Predicting the transmembrane secondary structure of ligand-gated ion channels

Recent mutational analyses of ligand-gated ion channels (LGICs) have demonstrated a plausible site of anesthetic action within their transmembrane domains. Although there is a consensus that the transmembrane domain is formed from four membrane-spanning segments, the secondary structure of these segments is not known. We utilized 10 state-of-the-art bioinformatics techniques to predict the transmembrane topology of the tetrameric regions within six members of the LGIC family that are relevant to anesthetic action. They are the human forms of the GABA alpha 1 receptor, the glycine alpha 1 receptor, the 5HT3 serotonin receptor, the nicotinic AChR alpha 4 and alpha 7 receptors and the Torpedo nAChR alpha 1 receptor. The algorithms utilized were HMMTOP, TMHMM, TMPred, PHDhtm, DAS, TMFinder, SOSUI, TMAP, MEMSAT and TOPPred2. The resulting predictions were superimposed on to a multiple sequence alignment of the six amino acid sequences created using the CLUSTAL W algorithm. There was a clear statistical consensus for the presence of four alpha helices in those regions experimentally thought to span the membrane. The consensus of 10 topology prediction techniques supports the hypothesis that the transmembrane subunits of the LGICs are tetrameric bundles of alpha helices.

Thiol-disulfide oxidoreductases are essential for the production of the lantibiotic sublancin 168

Thiol-disulfide oxidoreductases are required for disulfide bond formation in proteins that are exported from the cytoplasm. Four enzymes of this type, termed BdbA, BdbB, BdbC, and BdbD, have been identified in the Gram-positive eubacterium Bacillus subtilis. BdbC and BdbD have been shown to be critical for the folding of a protein required for DNA uptake during natural competence. In contrast, no function has been assigned so far to the BdbA and BdbB proteins. The bdbA and bdbB genes are located in one operon that also contains the genes specifying the lantibiotic sublancin 168 and the ATP-binding cassette transporter SunT. Interestingly sublancin 168 contains two disulfide bonds. The present studies demonstrate that SunT and BdbB, but not BdbA, are required for the production of active sublancin 168. In addition, the BdbB paralogue BdbC is at least partly able to replace BdbB in sublancin 168 production. These observations show the unprecedented involvement of thiol-disulfide oxidoreductases in the synthesis of a peptide antibiotic. Notably BdbB cannot complement BdbC in competence development, showing that these two closely related thiol-disulfide oxidoreductases have different, but partly overlapping, substrate specificities.

Basic charge clusters and predictions of membrane protein topology

The topology predictor SPLIT 4.0 (http://pref.etfos.hr) predicts the sequence location of transmembrane helices by performing an automatic selection of optimal amino acid attribute and corresponding preference functions. The best topological model is selected by choosing the highest absolute bias parameter that combines the bias in basic charge motifs and the bias in positive residues (the "positive inside rule") with the charge difference across the first transmembrane segment. Basic charge motifs, such as the BBB, BXXBB, and BBXXB motifs in alpha-helical integral membrane proteins, are significantly more frequent near cytoplasmic membrane surface than expected from the Arg/Lys (B) frequency. The predictor's accuracy is 99% for predicting 178 transmembrane helices in all membrane proteins or subunits of known 3D structure.

Identification of novel membrane proteins by searching for patterns in hydropathy profiles

A technique has been developed to search a proteome database for new members of a functional class of membrane protein. It takes advantage of the highly conserved secondary structure of functionally related membrane proteins. Such proteins typically have the same number of transmembrane domains located at similar relative positions in their polypeptide sequence. This gives rise to a characteristic pattern of peaks in their hydropathy profiles. To conduct a search, each member of a polypeptide database is converted to a hydropathy profile, peaks are automatically detected, and the pattern of peaks is compared with a template. A template was designed for the acetylcholine (ACh) and glycine receptors of the cys-loop receptor superfamily. The key feature was a closely spaced triplet of hydropathy peaks bracketed by deep valleys. When applied to the human proteome the search procedure retrieved 153 profiles with a receptor-like triplet of peaks. The approach was highly selective with 70% of the retrieved profiles annotated as known or putative receptors. These included ACh, glycine, gamma-amino butyric acid and serotonin receptors, which are all related by sequence. However, ionotropic glutamate receptors, which have almost no sequence homology with ACh receptors, were also retrieved. Thus, the strategy can find members of a functional class that cannot be identified by sequence alignment. To demonstrate that the strategy can easily be extended to other membrane protein families, a template was developed for the neurotransmitter/Na+ symporter family, and similar results were obtained. This approach should prove a useful adjunct to sequence-based retrieval tools when searching for novel membrane proteins.

Rapid topology mapping of Escherichia coli inner-membrane proteins by prediction and PhoA/GFP fusion analysis

We present an approach that allows rapid determination of the topology of Escherichia coli inner-membrane proteins by a combination of topology prediction and limited fusion-protein analysis. We derive new topology models for 12 inner-membrane proteins: MarC, PstA, TatC, YaeL, YcbM, YddQ, YdgE, YedZ, YgjV, YiaB, YigG, and YnfA. We estimate that our approach should make it possible to arrive at highly reliable topology models for roughly 10% of the approximately 800 inner-membrane proteins thought to exist in E. coli.

Dominant and recessive deafness caused by mutations of a novel gene, TMC1, required for cochlear hair-cell function

Positional cloning of hereditary deafness genes is a direct approach to identify molecules and mechanisms underlying auditory function. Here we report a locus for dominant deafness, DFNA36, which maps to human chromosome 9q13-21 in a region overlapping the DFNB7/B11 locus for recessive deafness. We identified eight mutations in a new gene, transmembrane cochlear-expressed gene 1 (TMC1), in a DFNA36 family and eleven DFNB7/B11 families. We detected a 1.6-kb genomic deletion encompassing exon 14 of Tmc1 in the recessive deafness (dn) mouse mutant, which lacks auditory responses and has hair-cell degeneration. TMC1 and TMC2 on chromosome 20p13 are members of a gene family predicted to encode transmembrane proteins. Tmc1 mRNA is expressed in hair cells of the postnatal mouse cochlea and vestibular end organs and is required for normal function of cochlear hair cells.

Conservation of amino acid transporters in fungi, plants and animals

When comparing the transporters of three completely sequenced eukaryotic genomes--Saccharomyces cerevisiae, Arabidopsis thaliana and Homo sapiens--transporter types can be distinguished according to phylogeny, substrate spectrum, transport mechanism and cell specificity. The known amino acid transporters belong to five different superfamilies. Two preferentially Na(+)-coupled transporter superfamilies are not represented in the yeast and Arabidopsis genomes, whereas the other three groups, which often function as H(+)-coupled systems, have members in all investigated genomes. Additional superfamilies exist for organellar transport, including mitochondrial and plastidic carriers. When used in combination with phylogenetic analyses, functional comparison might aid our prediction of physiological functions for related but uncharacterized open reading frames.

Beethoven, a mouse model for dominant, progressive hearing loss DFNA36

Despite recent progress in identifying genes underlying deafness, there are still relatively few mouse models of specific forms of human deafness. Here we describe the phenotype of the Beethoven (Bth) mouse mutant and a missense mutation in Tmc1 (transmembrane cochlear-expressed gene 1). Progressive hearing loss (DFNA36) and profound congenital deafness (DFNB7/B11) are caused by dominant and recessive mutations of the human ortholog, TMC1 (ref. 1), for which Bth and deafness (dn) are mouse models, respectively.

Prediction of organellar targeting signals

The subcellular location of a protein is an important characteristic with functional implications, and hence the problem of predicting subcellular localization from the amino acid sequence has received a fair amount of attention from the bioinformatics community. This review attempts to summarize the present state of the art in the field.

Formation of helical hairpins during membrane protein integration into the endoplasmic reticulum membrane. Role of the N and C-terminal flanking regions

The helical hairpin, two closely spaced transmembrane helices separated by a short turn, is a common structural element in integral membrane proteins. Previous studies on the sequence determinants of helical hairpin formation have focussed on the role of polar and charged residues placed centrally in a long stretch of hydrophobic residues, and have yielded a "propensity scale" for the relative efficiency with which different residues promote the formation of helical hairpins. In this study, we shift our attention to the role of charged residues flanking the hydrophobic stretch. Clusters of charged residues are known to hinder membrane translocation, and thus flanking charged residues may conceivably force a long hydrophobic segment to form a helical hairpin even if there are no or only weakly turn-promoting residues in the hydrophobic stretch. We indeed find that Lys and, more surprisingly, Asp residues strongly affect helical hairpin formation when placed next to a poly-Leu-based transmembrane segment. We also find that a cluster of four consecutive Lys residues can affect the efficiency of helical hairpin formation even when placed approximately 30 residues downstream of the hydrophobic stretch. These observations have interesting implications for the way we picture membrane protein topogenesis within the context of the endoplasmic reticulum (ER) translocon.