Flexible programs for the prediction of average amphipathicity of multiply aligned homologous proteins: application to integral membrane transport proteins

Simple flexible programs (TREEMOMENT and PILEUPMOMENT) are described for depicting the average amphipathicity (hydrophobic moment) along multiply aligned sequences of a family of evolutionarily related proteins. The programs are applicable to any number of aligned sequences and can be set for any desired angle corresponding to a residue repeat unit in a protein secondary structural element such as 100 per residue for an alpha- helix or 180 per residue for a beta-strand. These programs can be used to identify amphipathic regions common to the members of a protein family. The use of these programs is exemplified by showing that some families of integral membrane transport proteins (i.e. permeases of the bacterial phosphotransferase system (PTS) and the anion exchangers of animals) exhibit strikingly amphipathic alpha-helical structures immediately preceding the first hydrophobic transmembrane segment of their membrane-embedded domain(s). Other families, such as the major facilitator superfamily of uniporters, symporters and antiporters, do not exhibit this structural feature. The amphipathic structures in PTS permeases have been implicated in membrane insertion during biogenesis.     

Evolutionary relationships among the permease proteins of the bacterial phosphoenolpyruvate: Sugar phosphotransferase system. Construction of phylogenetic trees and possible relatedness to proteins of eukaryotic mitochondria

Summary The amino acid sequences of 15 sugar permeases of the bacterial phosphoenolpyruvatedependent phosphotransferase system (PTS) were divided into four homologous segments, and these segments were analyzed to give phylogenetic trees. The permease segments fell into four clusters: the lactose-cellobiose cluster, the fructose-mannitol cluster, the glucose-N-acetylglucosamine cluster, and the sucrose--glucoside cluster. Sequences of the glucitol and mannose permeases (clusters 5 and 6, respectively) were too dissimilar to establish homology with the other permeases, but short regions of statistically significant sequence similarities were noted. The functional and structural relationships of these permease segments are discussed.Some of the homologous PTS permeases were found to exhibit sufficient sequence similarity to subunits 4 and 5 of the eukaryotic mitochondrial NADH dehydrogenase complex to suggest homology. Moreover, subunits 4 and 5 of this complex appeared to be homologous to each other, suggesting that these PTS and mitochondrial proteins comprise a superfamily. The integral membrane subunits of the evolutionarily divergent mannose PTS permease, the P and M subunits, exhibited limited sequence similarity to subunit 6 of the mitochondrial F₁F₀-ATPase and subunit 5b of cytochrome oxidase, respectively. These results suggest that PTS sugar permeases and mitochondrial proton-translocating proteins may be related, although the possibility of convergent evolution cannot be ruled out.     

Bacterial proteins with N-terminal leader sequences resembling mitochondrial targeting sequences of eukaryotes

Amphipathic, alpha-helical, leader sequences, analogous to those that direct nuclear-encoded eukaryotic proteins into mitochondria, have been found in one and only one class of bacterial integral membrane proteins. These bacterial proteins are the sugar permeases of the phosphoenolpyruvate-dependent phosphotransferase system. The amphipathic leader sequence in each of these proteins is terminated by a helix breaker, either a prolyl residue or 2 adjacent glycyl residues. Preliminary evidence suggests that these leader sequences function to target the proteins to the envelope fraction of the prokaryotic cell during their biosynthesis.     

Sugar permeases of the bacterial phosphoenolpyruvate-dependent phosphotransferase system: sequence comparisons

The amino acyl sequences of eight permeases (enzymes II and enzyme II-III pairs) of the bacterial phosphoenolpyruvate:sugar phosphotransferase system (PTS) have been analyzed. All systems show similar sizes, and six of these systems exhibit the same molecular weight +/- 2%. Several exhibit sequence homology. Characteristic NH2-terminal and COOH-terminal sequences were found. The NH2-terminal leader sequences are believed to function in targeting of the permeases to the membrane, whereas the characteristic COOH-terminal sequences are postulated to mediate interaction with the energy-coupling protein phospho HPr. One of the systems, the one specific for mannose, exhibits distinctive characteristics. A pair of probable phosphorylation sites was detected in each of the five most similar systems, those specific for beta-glucosides, sucrose, glucose, N-acetylglucosamine, and mannitol. One of the two equivalent phosphorylation sites (proposed phosphorylation site 1) was located approximately 80 residues from the COOH terminus of each system. The other site (proposed phosphorylation site 2) was located approximately 440 residues from the COOH termini of the glucose and N-acetylglucosamine systems, approximately 320 residues from the COOH termini of the beta-glucoside and sucrose systems, and 381 residues from the COOH terminus of the mannitol system. Intragenic rearrangement during evolutionary history may account for the different positions of phosphorylation sites 2 in the different PTS permeases. More extensive intragenic rearrangements may have given rise to entirely different positions of phosphorylation in the glucitol, mannose, and lactose systems. A single, internal amphipathic alpha-helix with characteristic features was found in each of seven of the eight enzymes II. The lactose-specific enzyme III of Staphylococcus aureus was unique in possessing a COOH-terminal amphipathic alpha-helix rich in basic amino acyl residues. Possible functions for these amphipathic segments are discussed.     

The SMR family: a novel family of multidrug efflux proteins involved with the efflux of lipophilic drugs

The sequenced members of a novel family of small, hydrophobic, bacterial multidrug-resistance efflux proteins, which we have designated the small multidrug resistance (SMR) protein family, are identified and analysed. Two distinct clusters of proteins were identified within this family: (i) small multidrug efflux systems; and (ii) Sug proteins, potentially involved in the suppression of groEL mutations. Hydropathy and residue distribution analyses of this family suggest a structural model in which the polypeptide chain spans the membrane four times as mildly amphipathic α-helices. The roles of specific residues, a possible mechanistic model of drug efflux, and the primary physiological role(s) of the SMR proteins are discussed.     

Sequence relationships between integral inner membrane proteins of binding protein-dependent transport systems: evolution by recurrent gene duplications.

Periplasmic binding protein-dependent transport systems are composed of a periplasmic substrate-binding protein, a set of 2 (sometimes 1) very hydrophobic integral membrane proteins, and 1 (sometimes 2) hydrophilic peripheral membrane protein that binds and hydrolyzes ATP. These systems are members of the superfamily of ABC transporters. We performed a molecular phylogenetic analysis of the sequences of 70 hydrophobic membrane proteins of these transport systems in order to investigate their evolutionary history. Proteins were grouped into 8 clusters. Within each cluster, protein sequences displayed significant similarities, suggesting that they derive from a common ancestor. Most clusters contained proteins from systems transporting analogous substrates such as monosaccharides, oligopeptides, or hydrophobic amino acids, but this was not a general rule. Proteins from diverse bacteria are found within each cluster, suggesting that the ancestors of current clusters were present before the divergence of bacterial groups. The phylogenetic trees computed for hydrophobic membrane proteins of these permeases are similar to those described for the periplasmic substrate-binding proteins. This result suggests that the genetic regions encoding binding protein-dependent permeases evolved as whole units. Based on the results of the classification of the proteins and on the reconstructed phylogenetic trees, we propose an evolutionary scheme for periplasmic permeases. According to this model, it is probable that these transport systems derive from an ancestral system having only 1 hydrophobic membrane protein.(ABSTRACT TRUNCATED AT 250 WORDS)     

Identification of novel families of membrane proteins from the model plant Arabidopsis thaliana

MOTIVATION: The completion of the Arabidopsis genome offers the first opportunity to analyze all of the membrane protein sequences of a plant. The majority of integral membrane proteins including transporters, channels, and pumps contain hydrophobic alpha-helices and can be selected based on TransMembrane Spanning (TMS) domain prediction. By clustering the predicted membrane proteins based on sequence, it is possible to sort the membrane proteins into families of known function, based on experimental evidence or homology, or unknown function. This provides a way to identify target sequences for future functional analysis. RESULTS: An automated approach was used to select potential membrane protein sequences from the set of all predicted proteins and cluster the sequences into related families. The recently completed sequence of Arabidopsis thaliana, a model plant, was analyzed. Of the 25,470 predicted protein sequences 4589 (18%) were identified as containing two or more membrane spanning domains. The membrane protein sequences clustered into 628 distinct families containing 3208 sequences. Of these, 211 families (1764 sequences) either contained proteins of known function or showed homology to proteins of known function in other species. However, 417 families (1444 sequences) contained only sequences with no known function and no homology to proteins of known function. In addition, 1381 sequences did not cluster with any family and no function could be assigned to 1337 of these.     

Topological predictions for integral membrane permeases of the phosphoenolpyruvate:sugar phosphotransferase system

We report bioinformatic analyses of the largest superfamily of integral membrane permeases of the bacterial phosphotransferase system (PTS), the Enzyme IIC constituents of the Glc superfamily. Phylogenetic analyses reveal that this superfamily consists of five equally distant families, the Glucose (Glc), beta-Glucoside (Bgl), Fructose (Fru), Mannitol (Mtl) and Lactose (Lac) families. Average hydropathy, amphipathicity and similarity plots were generated for these five families as well as for the entire superfamily. Charged residue distribution was analyzed, and the most conserved sequence motif, common to all five families, was identified. The results show that the members of all five families exhibit similar average hydropathy plots with regions of average amphipathicity and relative conservation also being similar. Evidence is presented suggesting that the Glucitol (Gut) family of Enzyme IIC constituents is a distant member of the Glc superfamily. Based on our analyses we offer a topological model that resembles, but differs in detail from the two previously proposed models.     

Structural similarity between legumin and vicilin storage proteins from legumes.

The primary structures for several members of both the vicilin and legumin families of storage proteins were examined using a computer routine based on amino acid physical characteristics. The comparison algorithm revealed that sequences from the two families could be aligned and share a number of predicted secondary structural features. The COOH-terminal half of the subunits in both families displayed a highly conserved core region that was largely hydrophobic and in which a high proportion of the residues were predicted to be in beta-sheet conformations. The central region of the molecules which contained mixed areas of predicted helical and sheet conformations showed more variability in residue selection than the COOH-terminal regions. The NH2-terminal segments of subunits from the two different families could not be aligned though they characteristically had a high proportion of residues predicted to be in helical conformations. The feature which most clearly distinguished subunits between the two families was an inserted span in the legumin group with a high proportion of acidic amino acids located between the central and COOH-terminal domains. Residues in this insertion were predicted to exist mainly in helical conformation. Since considerable size variation occurs in this area amongst the legumin subunits, alterations in this region may have a minimal detrimental effect on the structure of the proteins.     

A conserved glutamate residue, Glu-257, is important for substrate binding and transport by the Escherichia coli mannitol permease.

The mannitol permease, or D-mannitol-specific enzyme II of the phosphoenolpyruvate-dependent carbohydrate phosphotransferase system (PTS) of Escherichia coli, both transports and phosphorylates its substrate. Previous analyses of the amino acid sequences of PTS permeases specific for various carbohydrates in different species of bacteria revealed several regions of similarity. The most highly conserved region includes a GIXE motif, in which the glutamate residue is completely conserved among the permeases that contain this motif. The corresponding residue in the E. coli mannitol permease is Glu-257, which is located in a large putative cytoplasmic loop of the transmembrane domain of the protein. We used site-directed mutagenesis to investigate the role of Glu-257. The properties of proteins with mutations at position 257 suggest that a carboxylate side chain at this position is essential for mannitol binding. E257A and E257Q mutant proteins did not bind mannitol detectably, while the E257D mutant could still bind this substrate. Kinetic studies with the E257D mutant protein also showed that a glutamate residue at position 257 of this permease is specifically required for efficient mannitol transport. While the E257D permease phosphorylated mannitol with kinetic parameters similar to those of the wild-type protein, the Vmax for mannitol uptake by this mutant protein is less than 5% that of the wild type. These results suggest that Glu-257 of the mannitol permease and the corresponding glutamate residues of other PTS permeases play important roles both in binding the substrate and in transporting it through the membrane.     

Statistical and functional analyses of viral and cellular proteins with N-terminal amphipathic alpha-helices with large hydrophobic moments: importance to macromolecular recognition and organelle targeting.

A total of 1,911 proteins with N-terminal methionyl residues were computer screened for potential N-terminal alpha-helices with strong amphipathic character. By the criteria of D. Eisenberg (Annu. Rev. Biochem. 53:595-623, 1984), only 3.5% of nonplastid, nonviral proteins exhibited potential N-terminal alpha-helices, 18 residues in length, with hydrophobic moment values per amino acyl residue ([muH]) in excess of 0.4. By contrast, 10% of viral proteins exhibited corresponding [muH] values in excess of 0.4. Of these viral proteins with known functions, 55% were found to interact functionally with nucleic acids, 30% were membrane-interacting proteins or their precursors, and 15% were structural proteins, primarily concerned with host cell interactions. These observations suggest that N-terminal amphipathic alpha-helices of viral proteins may (i) function in nucleic acid binding, (ii) facilitate membrane insertion, and (iii) promote host cell interactions. Analyses of potential amphipathic N-terminal alpha-helices of cellular proteins are also reported, and their significance to organellar or envelope targeting is discussed.     

Interactions involved in the realignment of membrane-associated helices. An investigation using oriented solid-state NMR and attenuated total reflection Fourier transform infrared spectroscopies

A series of histidine-containing peptides (LAH4X6) was designed to investigate the membrane interactions of selected side chains. To this purpose, their pH-dependent transitions from in-plane to transmembrane orientations were investigated by attenuated total reflection Fourier transform infrared and oriented solid-state NMR spectroscopies. Peptides of the same family have previously been shown to exhibit antibiotic and DNA transfection activities. Solution NMR spectroscopy indicates that these peptides form amphipathic helical structures in membrane environments, and the technique was also used to characterize the pK values of all histidines in the presence of detergent micelles. Whereas one face of the amphipathic helix is clearly hydrophobic, the opposite side is flanked by four histidines surrounding six leucine, alanine, glycine, tryptophan, or tyrosine residues, respectively. This diversity in peptide composition causes pronounced shifts in the midpoint pH of the in-plane to transmembrane helical transition, which is completely abolished for the peptides carrying the most hydrophilic amino acid residues. These properties open up a conceptually new approach to study in a quantitative manner the hydrophobic as well as specific interactions of amino acids in membranes. Notably, the resulting scale for whole residue transitions from the bilayer interface to the hydrophobic membrane interior is obtained from extended helical sequences in lipid bilayers.     

Evolution of permease diversity and energy-coupling mechanisms with special reference to the bacterial phosphotransferase system

Different classes of apparently unrelated permeases couple different forms of energy to solute transport. While the energy coupling mechanisms utilized by the different permease classes are clearly distinct, it is proposed, based on structural comparisons, that many of these permeases possess transmembrane, hydrophobic domains which are evolutionarily related. Carriers may have arisen from transmembrane pore-forming proteins, and the protein constituents or domains which are specifically responsible for energy coupling may have had distinct origins. Thus, complex permeases may possess mosaic structures. This suggestion is substantiated by recent findings regarding the evolutionary origins of the bacterial phosphoenolpyruvate-dependent phosphotransferase system (PTS). Mechanistic implications of this proposal are presented.     

Whole genome analyses of transporters in spirochetes: Borrelia burgdorferi and Treponema pallidum

The completely sequenced genomes of two spirochetes, Borrelia burgdorferi(Bbu) and Treponema pallidum (Tpa) were analyzed for the distribution of transporter types. Both organisms exhibited fewer proteins with >7 alpha-helical transmembrane spanners (TMSs), and fewer identified transport systems per megabase pair of DNA than most other prokaryotes analyzed. Each organism exhibits one recognizable ion channel protein of the MscS family. Tpa has twice as many primary carriers as Bbu but lacks PTS permeases that are plentiful in Bbu. Tpa is the only prokaryote so far sequenced which has two F-type ATPases. Large families of secondary nutrient uptake carriers (MFS and APC) that are prevalent in other organisms are essentially lacking in Spirochetes. The largest Spirochete secondary carrier families consist of efflux systems. While both Bbu and Tpa exhibit an unusual degree of transporter diversity, major differences in specificity exist between these two organisms.     

Topology of the hydrophobic membrane-bound components of the histidine periplasmic permease. Comparison with other members of the family.

The membrane-bound complex of periplasmic permeases comprises two hydrophobic proteins which have been hypothesized to be integral membrane-spaninning proteins. We have investigated the topological organization of the hydrophobic components of the Salmonella typhimurium histidine permease, HisQ and HisM. Both proteins are digested by trypsin and proteinase K when either inside-out or right-side-out membrane vesicles are used. Therefore, these proteins are exposed to both surfaces of the membrane. Digestion with carboxypeptidase and binding studies with antibodies directed against the carboxyl terminus of HisQ and HisM have localized their carboxyl termini to the inside surface of the cytoplasmic membrane. Aminopeptidase digestion suggests periplasmic localization of their amino termini. Alkaline phosphatase fusions to HisQ and HisM indicate the existence of five spanners in both proteins. The periodicity and orientation of spanners and loops in HisQ and HisM match those of the five carboxyl-terminal spanners of MalF, the only other hydrophobic component of the periplasmic permeases for which topological information is available. An alignment of the sequences of all known hydrophobic components of periplasmic permeases is presented which indicates clear conservation of secondary structure and some conservation of primary sequence. The structural conservation of the components is discussed, and a role for a hydrophilic loop containing a conserved sequence (the EAA loop) is proposed.     

Sequence and structure conservation in a protein core

In order to study structural aspects of sequence conservation in families of homologous proteins, we have analyzed structurally aligned sequences of 585 proteins grouped into 128 homologous families. The conservation of a residue in a family is defined as the average residue similarity in a given position of aligned sequences. The residue similarities were expressed in the form of log-odd substitution tables that take into account the environments of amino acids in three-dimensional structures. The protein core is defined as those residues that have less then 7% solvent accessibility. The density of a protein core is described in terms of atom packing, which is investigated as a criterion for residue substitution and conservation. Although there is no significant correlation between sequence conservation and average atom packing around nonpolar residues such as leucine, valine and isoleucine, a significant correlation is observed for polar residues in the protein core. This may be explained by the hydrogen bonds in which polar residues are involved; the better their protection from water access the more stable should be the structure in that position. Proteins 33:358–366, 1998. © 1998 Wiley-Liss, Inc.     

Structural basis for NHERF recognition by ERM proteins

The Na+/H+ exchanger regulatory factor (NHERF) is a key adaptor protein involved in the anchoring of ion channels and receptors to the actin cytoskeleton through binding to ERM (ezrin/radixin/moesin) proteins. NHERF binds the FERM domain of ERM proteins, although NHERF has no signature Motif-1 sequence for FERM binding found in adhesion molecules. The crystal structures of the radixin FERM domain complexed with the NHERF-1 and NHERF-2 C-terminal peptides revealed a peptide binding site of the FERM domain specific for the 13 residue motif MDWxxxxx(L/I)Fxx(L/F) (Motif-2), which is distinct from Motif-1. This Motif-2 forms an amphipathic alpha helix for hydrophobic docking to subdomain C of the FERM domain. This docking causes induced-fit conformational changes in subdomain C and affects binding to adhesion molecule peptides, while the two binding sites are not overlapped. Our studies provide structural paradigms for versatile ERM linkages between membrane proteins and the cytoskeleton.     

Amino acid sequence analysis of the annexin super-gene family of proteins.

The annexins are a widespread family of calcium-dependent membrane-binding proteins. No common function has been identified for the family and, until recently, no crystallographic data existed for an annexin. In this paper we draw together 22 available annexin sequences consisting of 88 similar repeat units, and apply the techniques of multiple sequence alignment, pattern matching, secondary structure prediction and conservation analysis to the characterisation of the molecules. The analysis clearly shows that the repeats cluster into four distinct families and that greatest variation occurs within the repeat 3 units. Multiple alignment of the 88 repeats shows amino acids with conserved physicochemical properties at 22 positions, with only Gly at position 23 being absolutely conserved in all repeats. Secondary structure prediction techniques identify five conserved helices in each repeat unit and patterns of conserved hydrophobic amino acids are consistent with one face of a helix packing against the protein core in predicted helices a, c, d, e. Helix b is generally hydrophobic in all repeats, but contains a striking pattern of repeat-specific residue conservation at position 31, with Arg in repeats 4 and Glu in repeats 2, but unconserved amino acids in repeats 1 and 3. This suggests repeats 2 and 4 may interact via a buried saltbridge. The loop between predicted helices a and b of repeat 3 shows features distinct from the equivalent loop in repeats 1, 2 and 4, suggesting an important structural and/or functional role for this region. No compelling evidence emerges from this study for uteroglobin and the annexins sharing similar tertiary structures, or for uteroglobin representing a derivative of a primordial one-repeat structure that underwent duplication to give the present day annexins. The analyses performed in this paper are re-evaluated in the Appendix, in the light of the recently published X-ray structure for human annexin V. The structure confirms most of the predictions and shows the power of techniques for the determination of tertiary structural information from the amino acid sequences of an aligned protein family.