Classification of 29 families of secondary transport proteins into a single structural class using hydropathy profile analysis

The ups and downs of elevator-type di-/tricarboxylate membrane transporters

The divalent anion sodium symporter (DASS) family contains both sodium-driven anion cotransporters and anion/anion exchangers. The family belongs to a broader ion transporter superfamily (ITS), which comprises 24 families of transporters, including those of AbgT antibiotic efflux transporters. The human proteins in the DASS family play major physiological roles and are drug targets. We recently determined multiple structures of the human sodium-dependent citrate transporter (NaCT) and the succinate/dicarboxylate transporter from Lactobacillus acidophilus (LaINDY). Structures of both proteins show high degrees of structural similarity to the previously determined VcINDY fold. Conservation between these DASS protein structures and those from the AbgT family indicates that the VcINDY fold represents the overall protein structure for the entire ITS. The new structures of NaCT and LaINDY are captured in the inward- or outward-facing conformations, respectively. The domain arrangements in these structures agree with a rigid body elevator-type transport mechanism for substrate translocation across the membrane. Two separate NaCT structures in complex with a substrate or an inhibitor allowed us to explain the inhibition mechanism and propose a detailed classification scheme for grouping disease-causing mutations in the human protein. Structural understanding of multiple kinetic states of DASS proteins is a first step toward the detailed characterization of their entire transport cycle.

Structural insights into the elevator-like mechanism of the sodium/citrate symporter CitS

The sodium-dependent citrate transporter of Klebsiella pneumoniae (KpCitS) belongs to the 2-hydroxycarboxylate transporter (2-HCT) family and allows the cell to use citrate as sole carbon and energy source in anaerobic conditions. Here we present crystal structures of KpCitS in citrate-bound outward-facing, citrate-bound asymmetric, and citrate-free inward-facing state. The structures reveal that the KpCitS dimerization domain remains stationary throughout the transport cycle due to a hydrogen bond network as well as extensive hydrophobic interactions. In contrast, its transport domain undergoes a ~35° rigid-body rotation and a ~17 Å translocation perpendicular to the membrane to expose the substrate-binding site alternately to either side of the membrane. Furthermore, homology models of two other 2-HCT proteins based on the KpCitS structure offer structural insights into their differences in substrate specificity at a molecular level. On the basis of our results and previous biochemical data, we propose that the activity of the 2-HCT CitS involves an elevator-like movement in which the transport domain itself traverses the lipid bilayer, carrying the substrate into the cell in a sodium-dependent manner.

Convergent evolution of the arginine deiminase pathway: the ArcD and ArcE arginine/ornithine exchangers

The arginine deiminase (ADI) pathway converts L-arginine into L-ornithine and yields 1 mol of ATP per mol of L-arginine consumed. The L-arginine/L-ornithine exchanger in the pathway takes up L-arginine and excretes L-ornithine from the cytoplasm. Analysis of the genomes of 1281 bacterial species revealed the presence of 124 arc gene clusters encoding the pathway. About half of the clusters contained the gene encoding the well-studied L-arginine/L-ornithine exchanger ArcD, while the other half contained a gene, termed here arcE, encoding a membrane protein that is not a homolog of ArcD. The arcE gene product of Streptococcus pneumoniae was shown to take up L-arginine and L-ornithine with affinities of 0.6 and 1 μmol/L, respectively, and to catalyze metabolic energy-independent, electroneutral exchange. ArcE of S. pneumoniae could replace ArcD in the ADI pathway of Lactococcus lactis and provided the cells with a growth advantage. In contrast to ArcD, ArcE catalyzed translocation of the pathway intermediate L-citrulline with high efficiency. A short version of the ADI pathway is proposed for L-citrulline catabolism and the presence of the evolutionary unrelated arcD and arcE genes in different organisms is discussed in the context of the evolution of the ADI pathway.

Structure and elevator mechanism of the Na+-citrate transporter CitS

The recently determined crystal structure of the bacterial Na⁺-citrate symporter CitS provides unexpected structural and mechanistic insights. The protein has a fold that has not been seen in other proteins, but the oligomeric state, domain organization and proposed transport mechanism strongly resemble those of the sodium-dicarboxylate symporter vcINDY, and the putative exporters YdaH and MtrF, thus hinting at convergence in structure and function. CitS and the related proteins are predicted to translocate their substrates by an elevator-like mechanism, in which a compact transport domain slides up and down through the membrane while the dimerization domain is stably anchored. Here we review the large body of available biochemical data on CitS in the light of the new crystal structure. We show that the biochemical data are fully consistent with the proposed elevator mechanism, but also demonstrate that the current structural data cannot explain how strict coupling of citrate and Na⁺ transport is achieved. We propose a testable model for the coupling mechanism.

Structural fold and binding sites of the human Na⁺-phosphate cotransporter NaPi-II

Phosphate plays essential biological roles and its plasma level in humans requires tight control to avoid bone loss (insufficiency) or vascular calcification (excess). Intestinal absorption and renal reabsorption of phosphate are mediated by members of the SLC34 family of sodium-coupled transporters (NaPi-IIa,b,c) whose membrane expression is regulated by various hormones, circulating proteins, and phosphate itself. Consequently, NaPi-II proteins are also potentially important pharmaceutical targets for controlling phosphate levels. Their crucial role in Pi homeostasis is underscored by pathologies resulting from naturally occurring SLC34 mutations and SLC34 knockout animals. SLC34 isoforms have been extensively studied with respect to transport mechanism and structure-function relationships; however, the three-dimensional structure is unknown. All SLC34 transporters share a duplicated motif comprising a glutamine followed by a stretch of threonine or serine residues, suggesting the presence of structural repeats as found in other transporter families. Nevertheless, standard bioinformatic approaches fail to clearly identify a suitable template for molecular modeling. Here, we used hydrophobicity profiles and hidden Markov models to define a structural repeat common to all SLC34 isoforms. Similar approaches identify a relationship with the core regions in a crystal structure of Vibrio cholerae Na(+)-dicarboxylate transporter VcINDY, from which we generated a homology model of human NaPi-IIa. The aforementioned SLC34 motifs in each repeat localize to the center of the model, and were predicted to form Na(+) and Pi coordination sites. Functional relevance of key amino acids was confirmed by biochemical and electrophysiological analysis of expressed, mutated transporters. Moreover, the validity of the predicted architecture is corroborated by extensive published structure-function studies. The model provides key information for elucidating the transport mechanism and predicts candidate substrate binding sites.

AlignMe--a membrane protein sequence alignment web server

We present a web server for pair-wise alignment of membrane protein sequences, using the program AlignMe. The server makes available two operational modes of AlignMe: (i) sequence to sequence alignment, taking two sequences in fasta format as input, combining information about each sequence from multiple sources and producing a pair-wise alignment (PW mode); and (ii) alignment of two multiple sequence alignments to create family-averaged hydropathy profile alignments (HP mode). For the PW sequence alignment mode, four different optimized parameter sets are provided, each suited to pairs of sequences with a specific similarity level. These settings utilize different types of inputs: (position-specific) substitution matrices, secondary structure predictions and transmembrane propensities from transmembrane predictions or hydrophobicity scales. In the second (HP) mode, each input multiple sequence alignment is converted into a hydrophobicity profile averaged over the provided set of sequence homologs; the two profiles are then aligned. The HP mode enables qualitative comparison of transmembrane topologies (and therefore potentially of 3D folds) of two membrane proteins, which can be useful if the proteins have low sequence similarity. In summary, the AlignMe web server provides user-friendly access to a set of tools for analysis and comparison of membrane protein sequences. Access is available at http://www.bioinfo.mpg.de/AlignMe

Asymmetric perturbations of signalling oligomers

This review focuses on rapid and reversible noncovalent interactions for symmetric oligomers of signalling proteins. Symmetry mismatch, transient symmetry breaking and asymmetric perturbations via chemical (ligand binding) and physical (electric or mechanic) effects can initiate the signalling events. Advanced biophysical methods can reveal not only structural symmetries of stable membrane-bound signalling proteins but also asymmetric functional transition states. Relevant techniques amenable to distinguish between symmetric and asymmetric architectures are discussed including those with the capability of capturing low-populated transient conformational states. Typical examples of signalling proteins are overviewed for symmetry breaking in dimers (GPCRs, growth factor receptors, transcription factors); trimers (acid-sensing ion channels); tetramers (voltage-gated cation channels, ionotropic glutamate receptor, CNG and CHN channels); pentameric ligand-gated and mechanosensitive channels; higher order oligomers (gap junction channel, chaperonins, proteasome, virus capsid); as well as primary and secondary transporters. In conclusion, asymmetric perturbations seem to play important functional roles in a broad range of communicating networks.

Structure and substrate-induced conformational changes of the secondary citrate/sodium symporter CitS revealed by electron crystallography

The secondary Na+/citrate symporter CitS of Klebsiella pneumoniae is the best-characterized member of the 2-hydroxycarboxylate transporter family. The recent projection structure gave insight into its overall structural organization. Here, we present the three-dimensional map of dimeric CitS obtained with electron crystallography. Each monomer has 13 a-helical transmembrane segments; six are organized in a distal helix cluster and seven in the central dimer interface domain. Based on structural analyses and comparison to VcINDY, we propose a molecular model for CitS, assign the helices, and demonstrate the internal structural symmetry. We also present projections of CitS in several conformational states induced by the presence and absence of sodium and citrate as substrates. Citrate binding induces a defined movement of a helices within the distal helical cluster. Based on this, we propose a substrate translocation site and conformational changes that are in agreement with the transport model of ‘‘alternating access’’.

Classification of α-helical membrane proteins using predicted helix architectures

Despite significant methodological advances in protein structure determination high-resolution structures of membrane proteins are still rare, leaving sequence-based predictions as the only option for exploring the structural variability of membrane proteins at large scale. Here, a new structural classification approach for α-helical membrane proteins is introduced based on the similarity of predicted helix interaction patterns. Its application to proteins with known 3D structure showed that it is able to reliably detect structurally similar proteins even in the absence of any sequence similarity, reproducing the SCOP and CATH classifications with a sensitivity of 65% at a specificity of 90%. We applied the new approach to enhance our comprehensive structural classification of α-helical membrane proteins (CAMPS), which is primarily based on sequence and topology similarity, in order to find protein clusters that describe the same fold in the absence of sequence similarity. The total of 151 helix architectures were delineated for proteins with more than four transmembrane segments. Interestingly, we observed that proteins with 8 and more transmembrane helices correspond to fewer different architectures than proteins with up to 7 helices, suggesting that in large membrane proteins the evolutionary tendency to re-use already available folds is more pronounced.

Molecular signatures identify a candidate target of balancing selection in an arcD-like gene of Staphylococcus epidermidis

A comparative population genetics study revealed high levels of nucleotide polymorphism and intermediate-frequency alleles in an arcC gene of Staphylococcus epidermidis, but not in a homologous gene of the more aggressive human pathogen, Staphylococcus aureus. Further investigation showed that the arcC genes used in the multilocus sequence typing schemes of these two species were paralogs. Phylogenetic analyses of arcC-containing loci, including the arginine catabolic mobile element, from both species, suggested that these loci had an eventful history involving gene duplications, rearrangements, deletions, and horizontal transfers. The peak signatures in the polymorphic S. epidermidis locus were traced to an arcD-like gene adjacent to arcC; these signatures consisted of unusually elevated Tajima's D and π/K ratios, which were robust to assumptions about recombination and species divergence time and among the most elevated in the S. epidermidis genome. Amino acid polymorphisms, including one that differed in polarity and hydropathy, were located in the peak signatures and defined two allelic lineages. Recombination events were detected between these allelic lineages and potential donors and recipients of S. epidermidis were identified in each case. By comparison, the orthologous gene of S. aureus showed no unusual signatures. The ArcD-like protein belonged to the unknown ion transporter 3 family and appeared to be unrelated to ArcD from the arginine deiminase pathway. These studies report the first comparative population genetics results for staphylococci and the first statistical evidence for a candidate target of balancing selection in S. epidermidis.

Projection structure of the secondary citrate/sodium symporter CitS at 6 Å resolution by electron crystallography

CitS from Klebsiella pneumoniae acts as a secondary symporter of citrate and sodium ions across the inner membrane of the host. The protein is the best characterized member of the 2-hydroxycarboxylate transporter family, while no experimental structural information at sub-nanometer resolution is available on this class of membrane proteins. Here, we applied electron crystallography to two-dimensional crystals of CitS. Carbon-film-adsorbed tubular two-dimensional crystals were studied by cryo-electron microscopy, producing the 6-Å-resolution projection structure of the membrane-embedded protein. In the p22(1)2(1)-symmetrized projection map, the predicted dimeric structure is clearly visible. Each monomeric unit can tentatively be interpreted as being composed of 11 transmembrane α-helices. In projection, CitS shows a high degree of structural similarity to NhaP1, the Na(+)/H(+) antiporter of Methanococcus jannaschii. We discuss possible locations for the dimer interface and models for the helical arrangements and domain organizations of the symporter based on existing models.

Turnover and accessibility of a reentrant loop of the Na+-glutamate transporter GltS are modulated by the central cytoplasmic loop

GltS of Escherichia coli is a secondary transporter that catalyzes Na+-glutamate symport. The structural model of GltS shows two homologous domains with inverted membrane topology that are connected by a central loop that resides in the cytoplasm. Each domain contains a reentrant loop structure. Accessibility of the Cys residues in two GltS mutants in which Pro351 and Asn356 in the reentrant loop in the C-terminal domain were replaced by Cys is demonstrated to be sensitive to the catalytic state supporting a role for the reentrant loops in catalysis. Saturating concentrations of the substrate L-glutamate protected both mutants against inactivation by thiol reagents, while the presence of the co-ion Na+ stimulated the inactivation of both mutants. Insertion of the 10 kDa biotin acceptor domain (BAD) of oxaloacetate decarboxylase of Klebsiella pneumoniae in the central cytoplasmic loop blocked the access pathway from the periplasmic side of the membrane to the cysteine residues in mutants P351C and N356C in the reentrant loop. Kinetically, insertion of BAD increased the maximal rate of uptake 2.7-fold while leaving the apparent affinity constants for L-glutamate and Na+ unaltered. The data suggests that insertion of BAD in the central loop results in conformational changes at the translocation site that lower the activation energy of the translocation step without affecting the access pathway from the periplasmic side for substrate and co-ions. It is concluded that changes in the central loop that connects the two domains may have a regulatory function on the activity of the transporter.

Membrane topology screen of secondary transport proteins in structural class ST[3] of the MemGen classification. Confirmation and structural diversity

The MemGen structural classification of membrane proteins groups families of proteins by hydropathy profile alignment. Class ST[3] of the MemGen classification contains 32 families of transporter proteins including the IT superfamily. Transporters from 19 different families in class ST[3] were evaluated by the TopScreen experimental topology screening method to verify the structural classification by MemGen. TopScreen involves the determination of the cellular disposition of three sites in the polypeptide chain of the proteins which allows for discrimination between different topology models. For nearly all transporters at least one of the predicted localizations is different in the models produced by MemGen and predictor TMHMM. Comparison to the experimental data showed that in all cases the prediction by MemGen was correct. It is concluded that the structural model available for transporters of the [st324]ESS and [st326]2HCT families is also valid for the other families in class ST[3]. The core structure of the model consists of two homologous domains, each containing 5 transmembrane segments, which have an opposite orientation in the membrane. A reentrant loop is present in between the 4th and 5th segments in each domain. Nearly all of the identified and experimentally confirmed structural variations involve additions of transmembrane segments at the boundaries of the core model, at the N- and C-termini or in between the two domains. Most remarkable is a domain swap in two subfamilies of the [st312]NHAC family that results in an inverted orientation of the proteins in the membrane.

Mechanism of citrate metabolism by an oxaloacetate decarboxylase-deficient mutant of Lactococcus lactis IL1403

Citrate metabolism in resting cells of Lactococcus lactis IL1403(pFL3) results in the formation of two end products from the intermediate pyruvate, acetoin and acetate (A. M. Pudlik and J. S. Lolkema, J. Bacteriol. 193:706-714, 2011). Pyruvate is formed from citrate following uptake by the transporter CitP through the subsequent actions of citrate lyase and oxaloacetate decarboxylase. The present study describes the metabolic response of L. lactis when oxaloacetate accumulates in the cytoplasm. The oxaloacetate decarboxylase-deficient mutant ILCitM(pFL3) showed nearly identical rates of citrate consumption, but the end product profile in the presence of glucose shifted from mainly acetoin to only acetate. In addition, in contrast to the parental strain, the mutant strain did not generate proton motive force. Citrate consumption by the mutant strain was coupled to the excretion of oxaloacetate, with a yield of 80 to 85%. Following citrate consumption, oxaloacetate was slowly taken up by the cells and converted to pyruvate by a cryptic decarboxylase and, subsequently, to acetate. The transport of oxaloacetate is catalyzed by CitP. The parental strain IL1403(pFL3) containing CitP consumed oxaloacetate, while the original strain, IL1403, not containing CitP, did not. Moreover, oxaloacetate consumption was enhanced in the presence of L-lactate, indicating exchange between oxaloacetate and L-lactate catalyzed by CitP. Hence, when oxaloacetate inadvertently accumulates in the cytoplasm, the physiological response of L. lactis is to excrete oxaloacetate in exchange with citrate by an electroneutral mechanism catalyzed by CitP. Subsequently, in a second step, oxaloacetate is taken up by CitP and metabolized to pyruvate and acetate.

Tripartite ATP-independent periplasmic (TRAP) transporters in bacteria and archaea

The tripartite ATP-independent periplasmic (TRAP) transporters are the best-studied family of substrate-binding protein (SBP)-dependent secondary transporters and are ubiquitous in prokaryotes, but absent from eukaryotes. They are comprised of an SBP of the DctP or TAXI families and two integral membrane proteins of unequal sizes that form the DctQ and DctM protein families, respectively. The SBP component has a structure comprised of two domains connected by a hinge that closes upon substrate binding. In DctP-TRAP transporters, substrate binding is mediated through a conserved and specific arginine/carboxylate interaction in the SBP. While the SBP component has now been relatively well characterized, the membrane components of TRAP transporters are still poorly understood both in terms of their structure and function. We review the expanding repertoire of substrates and physiological roles for experimentally characterized TRAP transporters in bacteria and discuss mechanistic aspects of these transporters using data primarily from the sialic acid-specific TRAP transporter SiaPQM from Haemophilus influenzae, which suggest that TRAP transporters are high-affinity, Na(+)-dependent unidirectional secondary transporters.

Role of isoleucine-554 in lithium binding by the Na+/dicarboxylate cotransporter NaDC1

Sodium-coupled transport of citric acid cycle intermediates, such as succinate and citrate, is mediated by the NaDC1 transporter located on the apical membrane of kidney proximal tubule and small intestine cells. Our previous study showed that transmembrane helix (TM) 11 of NaDC1 is important for sodium and lithium binding, as well as for determining citrate affinity [Kahn and Pajor (1999) Biochemistry 38, 6151]. In the present study, 21 amino acids in TM11 and the extracellular loop of NaDC1 were mutated one at a time to cysteine. All of the mutants were well expressed on the plasma membrane, but many of them had decreased transport activity. The G550C, W561C, and L568C mutants were inactive, suggesting that these residues may be critical for function. None of the cysteine mutants was sensitive to inhibition by the membrane-impermeant cysteine reagents, MTSET or MTSES, suggesting that the helix is inaccessible to the extracellular solvent. Although NaDC1 is inhibited by low concentrations of lithium in the presence of sodium, the I554C mutant was stimulated by lithium with a K(0.5) of 4.8 mM. The I554C mutant also had decreased affinity for sodium. We conclude that TM11 is likely to be an outer helix in NaDC1 that contains several residues critical for transport. Ile-554 in the middle of the helix may be an important determinant of cation affinity and selectivity, in particular the high affinity cation binding site that recognizes lithium.

Structural perspectives on secondary active transporters

Secondary active transporters catalyze the concentrative transport of substrates across lipid membranes by harnessing the energy of electrochemical ion gradients. These transporters bind their ligands on one side of the membrane, and undergo a global conformational change to release them on the other side of the membrane. Over the last few years, crystal structures have captured several bacterial secondary transporters in different states along their transport cycle, providing insight into possible molecular mechanisms. In this review, we summarize recent findings focusing on the emerging structural and mechanistic similarities between evolutionary diverse transporters. We also discuss the structural basis of substrate binding, ion coupling and inhibition viewed from the perspective of these similarities.

Cross-linking of trans reentrant loops in the Na(+)-citrate transporter CitS of Klebsiella pneumoniae

The membrane topology model of the Na(+)-citrate transporter CitS of Klebsiella pneumoniae shows a core of two homologous domains with opposite orientation in the membrane and each containing a so-called reentrant loop. A split version of CitS was constructed to study domain interactions and proximity relationships of the putative reentrant loops. Split CitS retained 50% transport activity of the wild-type version in membrane vesicles. Unspecific cross-linking of the purified complex with glutaraldehyde revealed a tetrameric complex with two N and two C domains corresponding to dimeric CitS. The separately expressed domains were not detected in the membrane. Strong interaction between the two domains followed from successful purification of the whole complex by Ni(2+)-NTA chromatography when only one domain was His-tagged. The presence of citrate and/or the co-ion Na(+) during purification did not seem to affect the interaction significantly. Successful disulfide cross-linking was obtained between single cysteine residues introduced in the highly conserved GGNG sequence motif at the vertex of the reentrant loop in the N domain and either of two endogenous cysteine residues at the base of the reentrant loop in the C domain. The disulfide bond was formed within one subunit in the dimer. A model is proposed in which the reentrant loops in the N and C domains are overlapping at the domain interface in the 3D structure where they form (part of) the translocation pathway for substrate and co-ions.

Structures of membrane proteins

In reviewing the structures of membrane proteins determined up to the end of 2009, we present in words and pictures the most informative examples from each family. We group the structures together according to their function and architecture to provide an overview of the major principles and variations on the most common themes. The first structures, determined 20 years ago, were those of naturally abundant proteins with limited conformational variability, and each membrane protein structure determined was a major landmark. With the advent of complete genome sequences and efficient expression systems, there has been an explosion in the rate of membrane protein structure determination, with many classes represented. New structures are published every month and more than 150 unique membrane protein structures have been determined. This review analyses the reasons for this success, discusses the challenges that still lie ahead, and presents a concise summary of the key achievements with illustrated examples selected from each class.

Rapid screening of membrane topology of secondary transport proteins

Limited experimental data may be very useful to discriminate between membrane topology models of membrane proteins derived from different methods. A membrane topology screening method is proposed by which the cellular disposition of three positions in a membrane protein are determined, the N- and the C-termini and a position in the middle of the protein. The method involves amplification of the encoding genes or gene fragments by PCR, rapid cloning in dedicated vectors by ligation independent cloning, and determination of the cellular disposition of the three sites using conventional techniques. The N-terminus was determined by labeling with a fluorescent probe, the central position and the C-terminus by the reporter fusion technique using alkaline phosphatase (PhoA) and green fluorescence protein (GFP) as reporters. The method was evaluated using 16 transporter proteins of known function from four different structural classes. For 13 proteins a complete set of three localizations was obtained. The experimental data was used to discriminate between membrane topology models predicted by TMHMM, a widely used predictor using the amino acid sequence as input and by MemGen that uses hydropathy profile alignment and known 3D structures or existing models. It follows that in those cases where the models from the two methods were similar, the models were consistent with the experimental data. In those cases where the models differed, the MemGen model agreed with the experimental data. Three more recent predictors, MEMSAT3, OCTOPUS and TOPCONS showed a significantly higher consistency with the experimental data than observed with TMHMM.

Projection structure by single-particle electron microscopy of secondary transport proteins GltT, CitS, and GltS

The structure of three secondary transporter proteins, GltT of Bacillus stearothermophilus, CitS of Klebsiella pneumoniae, and GltS of Escherichia coli, was studied. The proteins were purified to homogeneity in detergent solution by Ni(2+)-NTA affinity chromatography, and the complexes were determined by BN-PAGE to be trimeric, dimeric, and dimeric for GltT, CitS, and GltS, respectively. The subunit stoichiometry correlated with the binding affinity of the Ni(2+)-NTA resin for the protein complexes. Projection maps of negatively stained transporter particles were obtained by single-particle electron microscopy. Processing of the GltT particles revealed a projection map possessing 3-fold rotational symmetry, in good agreement with the trimer observed in the crystal structure of a homologous protein, Glt(Ph) of Pyrococcus horikoshii. The CitS protein showed up in two main views: as a kidney-shaped particle and a biscuit-shaped particle, both with a long axis of 160 A. The latter has a width of 84 A, the former of 92 A. Symmetry considerations identify the biscuit shape as a top view and the kidney shape as a side view from within the membrane. Combining the two images shows that the CitS dimer is a protein with a strong curvature at one side of the membrane and, at the opposite side, an indentation in the middle at the subunit interface. The GltS protein was shaped like CitS with dimensions of 145 A x 84 A. The shapes and dimensions of the CitS and GltS particles are consistent with a similar structure of these two unrelated proteins.

The major amino acid transporter superfamily has a similar core structure as Na+-galactose and Na+-leucine transporters

The sodium solute symporters (SSS) and neurotransmitter sodium symporters (NSS) are two families of secondary transporters that are not related in amino acid sequence. Nonetheless, recent crystal structures showed that the Na(+)/galactose (SSS) and Na(+)/leucine (NSS) transporters have similar core structures. The structural relatedness highlights the need for classification methods for membrane protein structures based on other criteria than amino acid similarity. Here, we demonstrate that a method based on hydropathy profile alignments convincingly identifies structural similarity between the NSS and SSS families. Most importantly, the method shows that one of the largest transporter families for which a crystal structure is elusive (the amino acid/polyamine/organocation or APC superfamily), also shares the similar core structure observed for the Na(+)/galactose and Na(+)/leucine transporters. The APC superfamily contains the major amino acid transporter families that are found throughout life. Insight into their structure will significantly facilitate the studies of this important group of transporters.

Evolution of antiparallel two-domain membrane proteins: tracing multiple gene duplication events in the DUF606 family

X-ray crystallography has revealed that many integral membrane proteins consist of two domains with a similar fold but opposite (antiparallel) orientation in the membrane. The proteins are believed to have evolved by gene duplication and gene fusion events from a dual topology ancestral membrane protein, that adapted both orientations in the membrane and formed antiparallel homodimers. Here, we present a detailed analysis of the DUF606 family of bacterial membrane proteins that contains the entire collection of intermediate states of such an evolutionary pathway: single genes that would code for dual topology homodimeric proteins, paired genes coding for homologous proteins with a fixed but opposite orientation in the membrane that would form heterodimers, and fused genes that encode antiparallel two-domain fusion proteins. Two types of paired genes can be discriminated corresponding to the order in which the genes coding for the two oppositely oriented proteins occur in the operon. On the protein level, the heterodimers resulting from the two types of gene pairs are indistinguishable. In contrast, two types of fused genes corresponding to the two possible orders in which the oppositely oriented domains are present in the encoded proteins, do result in discernible types of proteins. The large number of genetic and protein states in the DUF606 family allowed for a detailed phylogenic analysis that revealed a total of nine independent duplication events in the DUF606 family, five of which resulted in paired genes, and four resulted in fused genes. Noticeably, there was no evidence for a sequential mechanism in which fusions evolve from a pair of genes. Rather, an evolutionary mechanism is proposed by which antiparallel two-domain proteins are the direct result of a gene duplication event. Combining the phylogeny of proteins and hosting microorganisms allowed for a reconstruction of the evolutionary pathway.

Secondary transport of amino acids in prokaryotes

Amino acid transport is a ubiquitous phenomenon and serves a variety of functions in prokaryotes, including supply of carbon and nitrogen for catabolic and anabolic processes, pH homeostasis, osmoprotection, virulence, detoxification, signal transduction and generation of electrochemical ion gradients. Many of the participating proteins have eukaryotic relatives and are successfully used as model systems for exploration of transporter structure and function. Distribution, physiological roles, functional properties, and structure-function relationships of prokaryotic alpha-amino acid transporters are discussed.

Membrane protein prediction methods

We survey computational approaches that tackle membrane protein structure and function prediction. While describing the main ideas that have led to the development of the most relevant and novel methods, we also discuss pitfalls, provide practical hints and highlight the challenges that remain. The methods covered include: sequence alignment, motif search, functional residue identification, transmembrane segment and protein topology predictions, homology and ab initio modeling. In general, predictions of functional and structural features of membrane proteins are improving, although progress is hampered by the limited amount of high-resolution experimental information available. While predictions of transmembrane segments and protein topology rank among the most accurate methods in computational biology, more attention and effort will be required in the future to ameliorate database search, homology and ab initio modeling.

Membrane topology prediction by hydropathy profile alignment: membrane topology of the Na(+)-glutamate transporter GltS

Structural classification of families of membrane proteins by bioinformatics techniques has become a critical aspect of membrane protein research. We have proposed hydropathy profile alignment to identify structural homology between families of membrane proteins. Here, we demonstrate experimentally that two families of secondary transporters, the ESS and 2HCT families, indeed share similar folds. Members of the two families show highly similar hydropathy profiles but cannot be shown to be homologous by sequence similarity. A structural model was predicted for the ESS family transporters based upon an existing model of the 2HCT family transporters. In the model, the transporters fold into two domains containing five transmembrane segments and a reentrant or pore-loop each. The two pore-loops enter the membrane embedded part of the proteins from opposite sides of the membrane. The model was verified by accessibility studies of cysteine residues in single-Cys mutants of the Na+-glutamate transporter GltS of Escherichia coli, a member of the ESS family. Cysteine residues positioned in predicted periplasmic loops were accessible from the periplasm by a bulky, membrane-impermeable thiol reagent, while cysteine residues in cytoplasmic loops were not. Furthermore, two cysteine residues in the predicted pore-loop entering the membrane from the cytoplasmic side were shown to be accessible for small, membrane-impermeable thiol reagents from the periplasm, as was demonstrated before for the Na+-citrate transporter CitS of Klebsiella pneumoniae, a member of the 2HCT family. The data strongly suggests that GltS of the ESS family and CitS of the 2HCT family share the same fold as was predicted by comparing the averaged hydropathy profiles of the two families.

Functional characterization and Me2+ion specificity of a Ca2+?citrate transporter from Enterococcus faecalis

Secondary transporters of the bacterial CitMHS family transport citrate in complex with a metal ion. Different members of the family are specific for the metal ion in the complex and have been shown to transport Mg(2+)-citrate, Ca(2+)-citrate or Fe(3+)-citrate. The Fe(3+)-citrate transporter of Streptococcus mutans clusters on the phylogenetic tree on a separate branch with a group of transporters found in the phylum Firmicutes which are believed to be involved in anaerobic citrate degradation. We have cloned and characterized the transporter from Enterococcus faecalis EfCitH in this cluster. The gene was functionally expressed in Escherichia coli and studied using right-side-out membrane vesicles. The transporter catalyzes proton-motive-force-driven uptake of the Ca(2+)-citrate complex with an affinity constant of 3.5 microm. Homologous exchange is catalyzed with a higher efficiency than efflux down a concentration gradient. Analysis of the metal ion specificity of EfCitH activity in right-side-out membrane vesicles revealed a specificity that was highly similar to that of the Bacillus subtilis Ca(2+)-citrate transporter in the same family. In spite of the high sequence identity with the S. mutans Fe(3+)-citrate transporter, no transport activity with Fe(3+) (or Fe(2+)) could be detected. The transporter of E. faecalis catalyzes translocation of citrate in complex with Ca(2+), Sr(2+), Mn(2+), Cd(2+) and Pb(2+) and not with Mg(2+), Zn(2+), Ni(2+) and Co(2+). The specificity appears to correlate with the size of the metal ion in the complex.

Extra domains in secondary transport carriers and channel proteins

"Extra" domains in members of the families of secondary transport carrier and channel proteins provide secondary functions that expand, amplify or restrict the functional nature of these proteins. Domains in secondary carriers include TrkA and SPX domains in DASS family members, DedA domains in TRAP-T family members (both of the IT superfamily), Kazal-2 and PDZ domains in OAT family members (of the MF superfamily), USP, IIA(Fru) and TrkA domains in ABT family members (of the APC superfamily), ricin domains in OST family members, and TrkA domains in AAE family members. Some transporters contain highly hydrophilic domains consisting of multiple repeat units that can also be found in proteins of dissimilar function. Similarly, transmembrane alpha-helical channel-forming proteins contain unique, conserved, hydrophilic domains, most of which are not found in carriers. In some cases the functions of these domains are known. They may be ligand binding domains, phosphorylation domains, signal transduction domains, protein/protein interaction domains or complex carbohydrate-binding domains. These domains mediate regulation, subunit interactions, or subcellular targeting. Phylogenetic analyses show that while some of these domains are restricted to closely related proteins derived from specific organismal types, others are nearly ubiquitous within a particular family of transporters and occur in a tremendous diversity of organisms. The former probably became associated with the transporters late in the evolutionary process; the latter probably became associated with the carriers much earlier. These domains can be located at either end of the transporter or in a central region, depending on the domain and transporter family. These studies provide useful information about the evolution of extra domains in channels and secondary carriers and provide novel clues concerning function.

The 2-hydroxycarboxylate transporter family: physiology, structure, and mechanism

The 2-hydroxycarboxylate transporter family is a family of secondary transporters found exclusively in the bacterial kingdom. They function in the metabolism of the di- and tricarboxylates malate and citrate, mostly in fermentative pathways involving decarboxylation of malate or oxaloacetate. These pathways are found in the class Bacillales of the low-CG gram-positive bacteria and in the gamma subdivision of the Proteobacteria. The pathways have evolved into a remarkable diversity in terms of the combinations of enzymes and transporters that built the pathways and of energy conservation mechanisms. The transporter family includes H+ and Na+ symporters and precursor/product exchangers. The proteins consist of a bundle of 11 transmembrane helices formed from two homologous domains containing five transmembrane segments each, plus one additional segment at the N terminus. The two domains have opposite orientations in the membrane and contain a pore-loop or reentrant loop structure between the fourth and fifth transmembrane segments. The two pore-loops enter the membrane from opposite sides and are believed to be part of the translocation site. The binding site is located asymmetrically in the membrane, close to the interface of membrane and cytoplasm. The binding site in the translocation pore is believed to be alternatively exposed to the internal and external media. The proposed structure of the 2HCT transporters is different from any known structure of a membrane protein and represents a new structural class of secondary transporters.

Sequence and hydropathy profile analysis of two classes of secondary transporters

A structural class in the MemGen classification of membrane proteins is a set of evolutionary related proteins sharing a similar global fold. A structural class contains both closely related pairs of proteins for which homology is clear from sequence comparison and very distantly related pairs, for which it is not possible to establish homology based on sequence similarity alone. In the latter case the evolutionary link is based on hydropathy profile analysis. Here, we use these evolutionary related sets of proteins to analyze the relationship between E-values in BLAST searches, sequence similarities in multiple sequence alignments and structural similarities in hydropathy profile analyses. Two structural classes of secondary transporters termed ST[3], which includes the Ion Transporter (IT) superfamily and ST[4], which includes the DAACS family (TC# 2.A.23) were extracted from the NCBI protein database. ST[3] contains 2051 unique sequences distributed over 32 families and 59 subfamilies. ST[4] is a smaller class containing 399 unique sequences distributed over 2 families and 7 subfamilies. One subfamily in ST[4] contains a new class of binding protein dependent secondary transporters. Comparison of the averaged hydropathy profiles of the subfamilies in ST[3] and ST[4] revealed that the two classes represent different folds. Divergence of the sequences in ST[4] is much smaller than observed in ST[3], suggesting different constraints on the proteins during evolution. Analysis of the correlation between the evolutionary relationship of pairs of proteins in a class and the BLAST E-value revealed that: (i) the BLAST algorithm is unable to pick up the majority of the links between proteins in structural class ST[3], (ii) "low complexity filtering" and "composition based statistics" improve the specificity, but strongly reduce the sensitivity of BLAST searches for distantly related proteins, indicating that these filters are too stringent for the proteins analyzed, and (iii) the E-value cut-off, which may be used to evaluate evolutionary significance of a hit in a BLAST search is very different for the two structural classes of membrane proteins.

Secondary transporters of the 2HCT family contain two homologous domains with inverted membrane topology and trans re-entrant loops

The 2-hydroxycarboxylate transporter (2HCT) family of secondary transporters belongs to a much larger structural class of secondary transporters termed ST3 which contains about 2000 transporters in 32 families. The transporters of the 2HCT family are among the best studied in the class. Here we detect weak sequence similarity between the N- and C-terminal halves of the proteins using a sensitive method which uses a database containing the N- and C-terminal halves of all the sequences in ST3 and involves blast searches of each sequence in the database against the whole database. Unrelated families of secondary transporters of the same length and composition were used as controls. The sequence similarity involved major parts of the N- and C-terminal halves and not just a small stretch. The membrane topology of the homologous N- and C-terminal domains was deduced from the experimentally determined topology of the members of the 2HCT family. The domains consist of five transmembrane segments each and have opposite orientations in the membrane. The N terminus of the N-terminal domain is extracellular, while the N terminus of the C-terminal domain is cytoplasmic. The loops between the fourth and fifth transmembrane segment in each domain are well conserved throughout the class and contain a high fraction of residues with small side chains, Gly, Ala and Ser. Experimental work on the citrate transporter CitS in the 2HCT family indicates that the loops are re-entrant or pore loops. The re-entrant loops in the N- and C-terminal domains enter the membrane from opposite sides (trans-re-entrant loops). The combination of inverted membrane topology and trans-re-entrant loops represents a new fold for secondary transporters and resembles the structure of aquaporins and models proposed for Na+/Ca2+ exchangers.

Structural and mechanistic diversity of secondary transporters

Recent reports on the three-dimensional structure of secondary transporters have dramatically increased our knowledge of the translocation mechanism of ions and solutes. The structures of five transporters at atomic resolution have yielded four different folds and as many different translocation mechanisms. The structure of the glutamate transporter homologue GltPh confirmed the role of pore-loop structures as essential parts of the translocation mechanism in one family of secondary transporters. Biochemical evidence for pore-loop structures in several other families suggest that they might be common in secondary transporters, adding to the structural and mechanistic diversity of secondary transporters.

Loop VIII/IX of the Na+-Citrate Transporter CitS of Klebsiella pneumoniae Folds into an Amphipathic Surface Helix

The sodium ion-dependent citrate transporter CitS of Klebsiella pneumoniae is a member of the 2-hydroxycarboxylate transporter (2HCT) family whose members transport divalent citrate in symport with two sodium ions. Profiles of the hydrophobic moment suggested the presence of an amphipathic helical structure in the cytoplasmic loop between transmembrane segments (TMSs) VIII and IX (the AH loop) in all members of the family. Cysteine-scanning mutagenesis was used to study the secondary structure of the AH loop. We have mutated 20 successive residues into cysteine residues, characterized each of the mutants for its transport activity, and determined the accessibility of the residues. Three of the mutants, G324C, F331C, and F332C, had very low citrate transport activity, and two others, I321C and S333C, exhibited significantly decreased activity after treatment of right-side-out membranes with membrane permeable thiol reagent N-ethylmaleimide (NEM), but not with membrane impermeable 4-acetamido-4'-maleimidylstilbene-2,2'-disulfonic acid (AmdiS) and [2-(trimethylammonium)ethyl]methanethiosulfonate (MTSET). No protection against NEM was observed with citrate or sodium ions. Labeling of the cysteine residues in the 20 mutants with the fluorescent probe fluorescein 5-maleimide, in membrane vesicles with an inverted orientation, resulted in a clear periodicity in the accessibility of the residues. Residues expected to be at the hydrophobic face of the putative alpha-helix were not accessible for the label, whereas those at the hydrophilic face were easily accessed and labeled. Pretreatment of whole cells and inside-out membranes expressing the mutants with the membrane impermeable reagent AmdiS confirmed the cytoplasmic localization of the AH region. It is concluded that the loop between TMSs VIII and IX folds into an amphipathic surface helix.

The tyrosine decarboxylase operon of Lactobacillus brevis IOEB 9809: characterization and conservation in tyramine-producing bacteria

Bacterial genes of tyrosine decarboxylases were recently identified. Here we continued the sequencing of the tyrosine decarboxylase locus of Lactobacillus brevis IOEB 9809 and determined a total of 7979 bp. The sequence contained four complete genes encoding a tyrosyl-tRNA synthetase, the tyrosine decarboxylase, a probable tyrosine permease and a Na+/H+ antiporter. Rapid amplification of cDNA ends (RACE) was employed to determine the 5'-end of mRNAs containing the tyrosine decarboxylase gene. It was located only 34-35 nucleotides upstream of the start codon, suggesting that the preceding tyrosyl-tRNA synthetase gene was transcribed separately. In contrast, reverse transcription-polymerase chain reactions (RT-PCRs) carried out with primers designed to amplify regions spanning gene junctions showed that some mRNAs contained the four genes. Homology searches revealed similar clusters of four genes in the genome sequences of Enterococcus faecalis and Enterococcus faecium. Phylogenetic analyses supported the hypothesis that these genes evolved all together. These data suggest that bacterial tyrosine decarboxylases are encoded in an operon containing four genes.

The ion transporter superfamily

We define a novel superfamily of secondary carriers specific for cationic and anionic compounds, which we have termed the ion transporter (IT) superfamily. Twelve recognized and functionally defined families constitute this superfamily. We provide statistical sequence analyses demonstrating that these families were in fact derived from a common ancestor. Further, we characterize the 12 families in terms of (1) the known substrates transported, (2) the modes of transport and energy coupling mechanisms used, (3) the family sizes (in numbers of sequenced protein members in the current NCBI database), (4) the organismal distributions of the members of each family, (5) the size ranges of the constituent proteins, (6) the predicted topologies of these proteins, and (7) the occurrence of non-homologous auxiliary proteins that may either facilitate or be required for transport. No member of the superfamily is known to function in a capacity other than transport. Proteins in several of the constituent families are shown to have arisen by tandem intragenic duplication events, but topological variation has resulted from a variety of dissimilar genetic fusion, splicing and insertional events. The evolutionary relationships between the members of each family are defined, leading to predictions of functionally relevant orthologous relationships. Some but not all of the families include functionally dissimilar paralogues that arose by early extragenic duplication events.