Sequence and hydropathy profile analysis of two classes of secondary transporters

Genome-Wide Transcriptional Excavation of Dipsacus asperoides Unmasked both Cryptic Asperosaponin Biosynthetic Genes and SSR Markers

BACKGROUND

Dipsacus asperoides is a traditional Chinese medicinal crop. The root is generally used as a medicine and is frequently prescribed by Chinese doctors for the treatment of back pain, limb paralysis, flutter trauma, tendon injuries, and fractures. With the rapid development of bioinformatics, research has been focused on this species at the gene or molecular level. For purpose of fleshing out genome information about D. asperoides, in this paper we conducted transcriptome analysis of this species.

PRINCIPAL FINDINGS

To date, many genes encoding enzymes involved in the biosynthesis of triterpenoid saponins in D.asperoides have not been elucidated. Illumina paired-end sequencing was employed to probe D. asperoides's various enzymes associated with the relevant mesostate. A total of 30, 832,805 clean reads and de novo spliced 43,243 unigenes were obtained. Of all unigenes, only 8.27% (3578) were successfully annotated in total of seven public databases: Nr, Nt, Swiss-Prot, GO, KOG, KEGG, and Pfam, which might be attributed to the poor studies on D. asperoides. The candidate genes encoding enzymes involved in triterpenoid saponin biosynthesis were identified and experimentally verified by reverse transcription qPCR, encompassing nine cytochrome P450s and 17 UDP-glucosyltransferases. Specifically, unearthly putative genes involved in the glycosylation of hederagenin were acquired. Simultaneously, 4490 SSRs from 43,243 examined sequences were determined via bioinformatics analysis.

CONCLUSION

This study represents the first report on the use of the Illumina sequence platform on this crop at the transcriptome level. Our findings of candidate genes encoding enzymes involved in Dipsacus saponin VI biosynthes is provide novel information in efforts to further understand the triterpenoid metabolic pathway on this species. The initial genetics resources in this study will contribute significantly to the genetic breeding program of D. asperoides, and are beneficial for clinical diagnosis and treatment.

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

Alignment of helical membrane protein sequences using AlignMe

Few sequence alignment methods have been designed specifically for integral membrane proteins, even though these important proteins have distinct evolutionary and structural properties that might affect their alignments. Existing approaches typically consider membrane-related information either by using membrane-specific substitution matrices or by assigning distinct penalties for gap creation in transmembrane and non-transmembrane regions. Here, we ask whether favoring matching of predicted transmembrane segments within a standard dynamic programming algorithm can improve the accuracy of pairwise membrane protein sequence alignments. We tested various strategies using a specifically designed program called AlignMe. An updated set of homologous membrane protein structures, called HOMEP2, was used as a reference for optimizing the gap penalties. The best of the membrane-protein optimized approaches were then tested on an independent reference set of membrane protein sequence alignments from the BAliBASE collection. When secondary structure (S) matching was combined with evolutionary information (using a position-specific substitution matrix (P)), in an approach we called AlignMePS, the resultant pairwise alignments were typically among the most accurate over a broad range of sequence similarities when compared to available methods. Matching transmembrane predictions (T), in addition to evolutionary information, and secondary-structure predictions, in an approach called AlignMePST, generally reduces the accuracy of the alignments of closely-related proteins in the BAliBASE set relative to AlignMePS, but may be useful in cases of extremely distantly related proteins for which sequence information is less informative. The open source AlignMe code is available at https://sourceforge.net/projects/alignme/, and at http://www.forrestlab.org, along with an online server and the HOMEP2 data set.

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.

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.

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.

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.