The Pentameric Ligand-Gated Ion Channel Family: A New Member of the Voltage Gated Ion Channel Superfamily?

In this report we present seven lines of bioinformatic evidence supporting the conclusion that the Pentameric Ligand-gated Ion Channel (pLIC) Family is a member of the Voltage-gated Ion Channel (VIC) Superfamily. In our approach, we used the Transporter Classification Database (TCDB) as a reference and applied a series of bioinformatic methods to search for similarities between the pLIC family and members of the VIC superfamily. These include: (1) sequence similarity, (2) compatibility of topology and hydropathy profiles, (3) shared domains, (4) conserved motifs, (5) similarity of Hidden Markov Model profiles between families, (6) common 3D structural folds, and (7) clustering analysis of all families. Furthermore, sequence and structural comparisons as well as the identification of a 3-TMS repeat unit in the VIC superfamily suggests that the sixth transmembrane segment evolved into a re-entrant loop. This evidence suggests that the voltage-sensor domain and the channel domain have a common origin. The classification of the pLIC family within the VIC superfamily sheds light onto the topological origins of this family and its evolution, which will facilitate experimental verification and further research into this superfamily by the scientific community.

Sequence Similarity among Structural Repeats in the Piezo Family of Mechanosensitive Ion Channels.

Members of the Piezo family of mechanically activated cation channels are involved in multiple physiological processes in higher eukaryotes, including vascular development, cell differentiation, touch perception, hearing, and more, but they are also common in single-celled eukaryotic microorganisms. Mutations in these proteins in humans are associated with a variety of diseases, such as colorectal adenomatous polyposis, dehydrated hereditary stomatocytosis, and hereditary xerocytosis. Available 3D structures for Piezo proteins show nine regions of four transmembrane segments each that have the same fold. Despite the remarkable similarity among the nine characteristic structural repeats in the family, no significant sequence similarity among them has been reported. Using bioinformatics approaches and the Transporter Classification Database (TCDB) as reference, we reliably identified sequence similarity among repeats based on four lines of evidence: (1) hidden Markov model-profile similarities across repeats at the family level, (2) pairwise sequence similarities between different repeats across Piezo homologs, (3) Piezo-specific conserved sequence signatures that consistently identify the same regions across repeats, and (4) conserved residues that maintain the same orientation and location in 3D space.

Discovery and Characterization of the Phospholemman/SIMP/Viroporin Superfamily.

Using bioinformatic approaches, we present evidence of distant relatedness among the Ephemerovirus Viroporin family, the Rhabdoviridae Putative Viroporin U5 family, the Phospholemman family, and the Small Integral Membrane Protein family. Our approach is based on the transitivity property of homology complemented with five validation criteria: (1) significant sequence similarity and alignment coverage, (2) compatibility of topology of transmembrane segments, (3) overlap of hydropathy profiles, (4) conservation of protein domains, and (5) conservation of sequence motifs. Our results indicate that Pfam protein domains PF02038 and PF15831 can be found in or projected onto members of all four families. In addition, we identified a 26-residue motif conserved across the superfamily. This motif is characterized by hydrophobic residues that help anchor the protein to the membrane and charged residues that constitute phosphorylation sites. In addition, all members of the four families with annotated function are either responsible for or affect the transport of ions into and/or out of the cell. Taken together, these results justify the creation of the novel Phospholemman/SIMP/Viroporin superfamily. Given that transport proteins can be found not just in cells, but also in viruses, the ability to relate viroporin protein families with their eukaryotic and bacterial counterparts is an important development in this superfamily.

Comparative population genomic analyses of transporters within the Asgard archaeal superphylum.

Upon discovery of the first archaeal species in the 1970s, life has been subdivided into three domains: Eukarya, Archaea, and Bacteria. However, the organization of the three-domain tree of life has been challenged following the discovery of archaeal lineages such as the TACK and Asgard superphyla. The Asgard Superphylum has emerged as the closest archaeal ancestor to eukaryotes, potentially improving our understanding of the evolution of life forms. We characterized the transportomes and their substrates within four metagenome-assembled genomes (MAGs), that is, Odin-, Thor-, Heimdall- and Loki-archaeota as well as the fully sequenced genome of Candidatus Prometheoarchaeum syntrophicum strain MK-D1 that belongs to the Loki phylum. Using the Transporter Classification Database (TCDB) as reference, candidate transporters encoded within the proteomes were identified based on sequence similarity, alignment coverage, compatibility of hydropathy profiles, TMS topologies and shared domains. Identified transport systems were compared within the Asgard superphylum as well as within dissimilar eukaryotic, archaeal and bacterial organisms. From these analyses, we infer that Asgard organisms rely mostly on the transport of substrates driven by the proton motive force (pmf), the proton electrochemical gradient which then can be used for ATP production and to drive the activities of secondary carriers. The results indicate that Asgard archaea depend heavily on the uptake of organic molecules such as lipid precursors, amino acids and their derivatives, and sugars and their derivatives. Overall, the majority of the transporters identified are more similar to prokaryotic transporters than eukaryotic systems although several instances of the reverse were documented. Taken together, the results support the previous suggestions that the Asgard superphylum includes organisms that are largely mixotrophic and anaerobic but more clearly define their metabolic potential while providing evidence regarding their relatedness to eukaryotes.

The Transporter Classification Database (TCDB): 2021 update.

The Transporter Classification Database (TCDB; tcdb.org) is a freely accessible reference resource, which provides functional, structural, mechanistic, medical and biotechnological information about transporters from organisms of all types. TCDB is the only transport protein classification database adopted by the International Union of Biochemistry and Molecular Biology (IUBMB) and now (October 1, 2020) consists of 20 653 proteins classified in 15 528 non-redundant transport systems with 1567 tabulated 3D structures, 18 336 reference citations describing 1536 transporter families, of which 26% are members of 82 recognized superfamilies. Overall, this is an increase of over 50% since the last published update of the database in 2016. This comprehensive update of the database contents and features include (i) adoption of a chemical ontology for substrates of transporters, (ii) inclusion of new superfamilies, (iii) a domain-based characterization of transporter families for the identification of new members as well as functional and evolutionary relationships between families, (iv) development of novel software to facilitate curation and use of the database, (v) addition of new subclasses of transport systems including 11 novel types of channels and 3 types of group translocators and (vi) the inclusion of many man-made (artificial) transmembrane pores/channels and carriers.

Expansion of the Transporter-Opsin-G protein-coupled receptor superfamily with five new protein families.

Here we provide bioinformatic evidence that the Organo-Arsenical Exporter (ArsP), Endoplasmic Reticulum Retention Receptor (KDELR), Mitochondrial Pyruvate Carrier (MPC), L-Alanine Exporter (AlaE), and the Lipid-linked Sugar Translocase (LST) protein families are members of the Transporter-Opsin-G Protein-coupled Receptor (TOG) Superfamily. These families share domains homologous to well-established TOG superfamily members, and their topologies of transmembranal segments (TMSs) are compatible with the basic 4-TMS repeat unit characteristic of this Superfamily. These repeat units tend to occur twice in proteins as a result of intragenic duplication events, often with subsequent gain/loss of TMSs in many superfamily members. Transporters within the ArsP family allow microbial pathogens to expel toxic arsenic compounds from the cell. Members of the KDELR family are involved in the selective retrieval of proteins that reside in the endoplasmic reticulum. Proteins of the MPC family are involved in the transport of pyruvate into mitochondria, providing the organelle with a major oxidative fuel. Members of family AlaE excrete L-alanine from the cell. Members of the LST family are involved in the translocation of lipid-linked glucose across the membrane. These five families substantially expand the range of substrates of transport carriers in the superfamily, although KDEL receptors have no known transport function. Clustering of protein sequences reveals the relationships among families, and the resulting tree correlates well with the degrees of sequence similarity documented between families. The analyses and programs developed to detect distant relatedness, provide insights into the structural, functional, and evolutionary relationships that exist between families of the TOG superfamily, and should be of value to many other investigators.

Expansion of the Major Facilitator Superfamily (MFS) to include novel transporters as well as transmembrane-acting enzymes.

The Major Facilitator Superfamily (MFS) is currently the largest characterized superfamily of transmembrane secondary transport proteins. Its diverse members are found in essentially all organisms in the biosphere and function by uniport, symport, and/or antiport mechanisms. In 1993 we first named and described the MFS which then consisted of 5 previously known families that had not been known to be related, and by 2012 we had identified a total of 74 families, classified phylogenetically within the MFS, all of which included only transport proteins. This superfamily has since expanded to 89 families, all included under TC# 2.A.1, and a few transporter families outside of TC# 2.A.1 were identified as members of the MFS. In this study, we assign nine previously unclassified protein families in the Transporter Classification Database (TCDB; http://www.tcdb.org) to the MFS based on multiple criteria and bioinformatic methodologies. In addition, we find integral membrane domains distantly related to partial or full-length MFS permeases in Lysyl tRNA Synthases (TC# 9.B.111), Lysylphosphatidyl Glycerol Synthases (TC# 4.H.1), and cytochrome b561 transmembrane electron carriers (TC# 5.B.2). Sequence alignments, overlap of hydropathy plots, compatibility of repeat units, similarity of complexity profiles of transmembrane segments, shared protein domains and 3D structural similarities between transport proteins were analyzed to assist in inferring homology. The MFS now includes 105 families.

Cryo-EM structure of OSCA1.2 from Oryza sativa elucidates the mechanical basis of potential membrane hyperosmolality gating.

Sensing and responding to environmental water deficiency and osmotic stresses are essential for the growth, development, and survival of plants. Recently, an osmolality-sensing ion channel called OSCA1 was discovered that functions in sensing hyperosmolality in Arabidopsis Here, we report the cryo-electron microscopy (cryo-EM) structure and function of an OSCA1 homolog from rice (Oryza sativa; OsOSCA1.2), leading to a model of how it could mediate hyperosmolality sensing and transport pathway gating. The structure reveals a dimer; the molecular architecture of each subunit consists of 11 transmembrane (TM) helices and a cytosolic soluble domain that has homology to RNA recognition proteins. The TM domain is structurally related to the TMEM16 family of calcium-dependent ion channels and lipid scramblases. The cytosolic soluble domain possesses a distinct structural feature in the form of extended intracellular helical arms that are parallel to the plasma membrane. These helical arms are well positioned to potentially sense lateral tension on the inner leaflet of the lipid bilayer caused by changes in turgor pressure. Computational dynamic analysis suggests how this domain couples to the TM portion of the molecule to open a transport pathway. Hydrogen/deuterium exchange mass spectrometry (HDXMS) experimentally confirms the conformational dynamics of these coupled domains. These studies provide a framework to understand the structural basis of proposed hyperosmolality sensing in a staple crop plant, extend our knowledge of the anoctamin superfamily important for plants and fungi, and provide a structural mechanism for potentially translating membrane stress to transport regulation.

Bioinformatic characterization of the Anoctamin Superfamily of Ca2+-activated ion channels and lipid scramblases.

Our laboratory has developed bioinformatic strategies for identifying distant phylogenetic relationships and characterizing families and superfamilies of transport proteins. Results using these tools suggest that the Anoctamin Superfamily of cation and anion channels, as well as lipid scramblases, includes three functionally characterized families: the Anoctamin (ANO), Transmembrane Channel (TMC) and Ca2+-permeable Stress-gated Cation Channel (CSC) families; as well as four families of functionally uncharacterized proteins, which we refer to as the Anoctamin-like (ANO-L), Transmembrane Channel-like (TMC-L), and CSC-like (CSC-L1 and CSC-L2) families. We have constructed protein clusters and trees showing the relative relationships among the seven families. Topological analyses suggest that the members of these families have essentially the same topologies. Comparative examination of these homologous families provides insight into possible mechanisms of action, indicates the currently recognized organismal distributions of these proteins, and suggests drug design potential for the disease-related channel proteins.

Comparative Analyses of Transport Proteins Encoded within the Genomes of Bdellovibrio bacteriovorus HD100 and Bdellovibrio exovorus JSS.

Bdellovibrio, δ-proteobacteria, including B. bacteriovorus (Bba) and B. exovorus (Bex), are obligate predators of other Gram-negative bacteria. While Bba grows in the periplasm of the prey cell, Bex grows externally. We have analyzed and compared the transport proteins of these 2 organisms based on the current contents of the Transporter Classification Database (TCDB; www.tcdb.org). Bba has 103 transporters more than Bex, 50% more secondary carriers, and 3 times as many MFS carriers. Bba has far more metabolite transporters than Bex as expected from its larger genome, but there are 2 times more carbohydrate uptake and drug efflux systems, and 3 times more lipid transporters. Bba also has polyamine and carboxylate transporters lacking in Bex. Bba has more than twice as many members of the Mot-Exb family of energizers, but both may have energizers for gliding motility. They use entirely different types of systems for iron acquisition. Both contain unexpectedly large numbers of pseudogenes and incomplete systems, suggesting that they are undergoing genome size reduction. Interestingly, all 5 outer-membrane receptors in Bba are lacking in Bex. The 2 organisms have similar numbers and types of peptide and amino acid uptake systems as well as protein and carbohydrate secretion systems. The differences observed correlate with and may account, in part, for the different lifestyles of these 2 bacterial predators.

The Membrane Attack Complex/Perforin Superfamily.

The membrane attack complex/perforin (MACPF) superfamily consists of a diverse group of proteins involved in bacterial pathogenesis and sporulation as well as eukaryotic immunity, embryonic development, neural migration and fruiting body formation. The present work shows that the evolutionary relationships between the members of the superfamily, previously suggested by comparison of their tertiary structures, can also be supported by analyses of their primary structures. The superfamily includes the MACPF family (TC 1.C.39), the cholesterol-dependent cytolysin (CDC) family (TC 1.C.12.1 and 1.C.12.2) and the pleurotolysin pore-forming (pleurotolysin B) family (TC 1.C.97.1), as revealed by expansion of each family by comparison against a large protein database, and by the comparisons of their hidden Markov models. Clustering analyses demonstrated grouping of the CDC homologues separately from the 12 MACPF subfamilies, which also grouped separately from the pleurotolysin B family. Members of the MACPF superfamily revealed a remarkably diverse range of proteins spanning eukaryotic, bacterial, and archaeal taxonomic domains, with notable variations in protein domain architectures. Our strategy should also be helpful in putting together other highly divergent protein families.

Characterization of the Tetraspan Junctional Complex (4JC) superfamily.

Connexins or innexins form gap junctions, while claudins and occludins form tight junctions. In this study, statistical data, derived using novel software, indicate that these four junctional protein families and eleven other families of channel and channel auxiliary proteins are related by common descent and comprise the Tetraspan (4 TMS) Junctional Complex (4JC) Superfamily. These proteins all share similar 4 transmembrane α-helical (TMS) topologies. Evidence is presented that they arose via an intragenic duplication event, whereby a 2 TMS-encoding genetic element duplicated tandemly to give 4 TMS proteins. In cases where high resolution structural data were available, the conclusion of homology was supported by conducting structural comparisons. Phylogenetic trees reveal the probable relationships of these 15 families to each other. Long homologues containing fusions to other recognizable domains as well as internally duplicated or fused domains are reported. Large "fusion" proteins containing 4JC domains proved to fall predominantly into family-specific patterns as follows: (1) the 4JC domain was N-terminal; (2) the 4JC domain was C-terminal; (3) the 4JC domain was duplicated or occasionally triplicated and (4) mixed fusion types were present. Our observations provide insight into the evolutionary origins and subfunctions of these proteins as well as guides concerning their structural and functional relationships.

Inferring molecular function: contributions from functional linkages.

In the current era of high-throughput sequencing and structure determination, functional annotation has become a bottleneck in biomedical science. Here, we show that automated inference of molecular function using functional linkages among genes increases the accuracy of functional assignments by > or =8% and enriches functional descriptions in > or =34% of top assignments. Furthermore, biochemical literature supports >80% of automated inferences for previously unannotated proteins. These results emphasize the benefit of incorporating functional linkages in protein annotation.

The partitioned Rhizobium etli genome: genetic and metabolic redundancy in seven interacting replicons.

We report the complete 6,530,228-bp genome sequence of the symbiotic nitrogen fixing bacterium Rhizobium etli. Six large plasmids comprise one-third of the total genome size. The chromosome encodes most functions necessary for cell growth, whereas few essential genes or complete metabolic pathways are located in plasmids. Chromosomal synteny is disrupted by genes related to insertion sequences, phages, plasmids, and cell-surface components. Plasmids do not show synteny, and their orthologs are mostly shared by accessory replicons of species with multipartite genomes. Some nodulation genes are predicted to be functionally related with chromosomal loci encoding for the external envelope of the bacterium. Several pieces of evidence suggest an exogenous origin for the symbiotic plasmid (p42d) and p42a. Additional putative horizontal gene transfer events might have contributed to expand the adaptive repertoire of R. etli, because they include genes involved in small molecule metabolism, transport, and transcriptional regulation. Twenty-three putative sigma factors, numerous isozymes, and paralogous families attest to the metabolic redundancy and the genomic plasticity necessary to sustain the lifestyle of R. etli in symbiosis and in the soil.

Evolutionary, structural and functional relationships revealed by comparative analysis of syntenic genes in Rhizobiales.

BACKGROUND: Comparative genomics has provided valuable insights into the nature of gene sequence variation and chromosomal organization of closely related bacterial species. However, questions about the biological significance of gene order conservation, or synteny, remain open. Moreover, few comprehensive studies have been reported for rhizobial genomes. RESULTS: We analyzed the genomic sequences of four fast growing Rhizobiales (Sinorhizobium meliloti, Agrobacterium tumefaciens, Mesorhizobium loti and Brucella melitensis). We made a comprehensive gene classification to define chromosomal orthologs, genes with homologs in other replicons such as plasmids, and those which were species-specific. About two thousand genes were predicted to be orthologs in each chromosome and about 80% of these were syntenic. A striking gene colinearity was found in pairs of organisms and a large fraction of the microsyntenic regions and operons were similar. Syntenic products showed higher identity levels than non-syntenic ones, suggesting a resistance to sequence variation due to functional constraints; also, an unusually high fraction of syntenic products contained membranal segments. Syntenic genes encode a high proportion of essential cell functions, presented a high level of functional relationships and a very low horizontal gene transfer rate. The sequence variability of the proteins can be considered the species signature in response to specific niche adaptation. Comparatively, an analysis with genomes of Enterobacteriales showed a different gene organization but gave similar results in the synteny conservation, essential role of syntenic genes and higher functional linkage among the genes of the microsyntenic regions. CONCLUSION: Syntenic bacterial genes represent a commonly evolved group. They not only reveal the core chromosomal segments present in the last common ancestor and determine the metabolic characteristics shared by these microorganisms, but also show resistance to sequence variation and rearrangement, possibly due to their essential character. In Rhizobiales and Enterobacteriales, syntenic genes encode a high proportion of essential cell functions and presented a high level of functional relationships.

Sequencing and analysis of common bean ESTs. Building a foundation for functional genomics.

Although common bean (Phaseolus vulgaris) is the most important grain legume in the developing world for human consumption, few genomic resources exist for this species. The objectives of this research were to develop expressed sequence tag (EST) resources for common bean and assess nodule gene expression through high-density macroarrays. We sequenced a total of 21,026 ESTs derived from 5 different cDNA libraries, including nitrogen-fixing root nodules, phosphorus-deficient roots, developing pods, and leaves of the Mesoamerican genotype, Negro Jamapa 81. The fifth source of ESTs was a leaf cDNA library derived from the Andean genotype, G19833. Of the total high-quality sequences, 5,703 ESTs were classified as singletons, while 10,078 were assembled into 2,226 contigs producing a nonredundant set of 7,969 different transcripts. Sequences were grouped according to 4 main categories, metabolism (34%), cell cycle and plant development (11%), interaction with the environment (19%), and unknown function (36%), and further subdivided into 15 subcategories. Comparisons to other legume EST projects suggest that an entirely different repertoire of genes is expressed in common bean nodules. Phaseolus-specific contigs, gene families, and single nucleotide polymorphisms were also identified from the EST collection. Functional aspects of individual bean organs were reflected by the 20 contigs from each library composed of the most redundant ESTs. The abundance of transcripts corresponding to selected contigs was evaluated by RNA blots to determine whether gene expression determined by laboratory methods correlated with in silico expression. Evaluation of root nodule gene expression by macroarrays and RNA blots showed that genes related to nitrogen and carbon metabolism are integrated for ureide production. Resources developed in this project provide genetic and genomic tools for an international consortium devoted to bean improvement.

Successful lateral transfer requires codon usage compatibility between foreign genes and recipient genomes.

We present evidence supporting the notion that codon usage (CU) compatibility between foreign genes and recipient genomes is an important prerequisite to assess the selective advantage of imported functions, and therefore to increase the fixation probability of horizontal gene transfer (HGT) events. This contrasts with the current tendency in research to predict recent HGTs in prokaryotes by assuming that acquired genes generally display poor CU. By looking at the CU level (poor, typical, or rich) exhibited by putative xenologs still resembling their original CU, we found that most alien genes predominantly present typical CU immediately upon introgression, thereby suggesting that the role of CU amelioration in HGT has been overemphasized. In our strategy, we first scanned a representative set of 103 complete prokaryotic genomes for all pairs of candidate xenologs (exported/imported genes) displaying similar CU. We applied additional filtering criteria, including phylogenetic validations, to enhance the reliability of our predictions. Our approach makes no assumptions about the CU of foreign genes being typical or atypical within the recipient genome, thus providing a novel unbiased framework to study the evolutionary dynamics of HGT.

Evolutionary and structural diversity of fungal laccases.

Fungal laccases have been extensively exploited for industrial purposes and there is a wealth of information available regarding their reaction mechanism, biological role and several molecular aspects, including cloning, heterologous expression and transcriptional analyses. Here we present the reconstruction of the fungal laccase loci evolution inferred from the comparative analysis of 48 different sequences. The topology of the phylogenetic trees indicate that a single monophyletic branch exists for fungal laccases and that laccase isozyme genes may have evolved independently, possibly through duplication-divergence events. Laccases are copper-containing enzymes generally identified by the utilization of substituted p-diphenol substrates. Interestingly, our approach permitted the assignment of two copper-containing oxidases, preliminarily catalogued as laccases, to a different evolutionary group, distantly related to the main branch of bona fide laccases.

The mosaic structure of the symbiotic plasmid of Rhizobium etli CFN42 and its relation to other symbiotic genome compartments.

BACKGROUND: Symbiotic bacteria known as rhizobia interact with the roots of legumes and induce the formation of nitrogen-fixing nodules. In rhizobia, essential genes for symbiosis are compartmentalized either in symbiotic plasmids or in chromosomal symbiotic islands. To understand the structure and evolution of the symbiotic genome compartments (SGCs), it is necessary to analyze their common genetic content and organization as well as to study their differences. To date, five SGCs belonging to distinct species of rhizobia have been entirely sequenced. We report the complete sequence of the symbiotic plasmid of Rhizobium etli CFN42, a microsymbiont of beans, and a comparison with other SGC sequences available. RESULTS: The symbiotic plasmid is a circular molecule of 371,255 base-pairs containing 359 coding sequences. Nodulation and nitrogen-fixation genes common to other rhizobia are clustered in a region of 125 kilobases. Numerous sequences related to mobile elements are scattered throughout. In some cases the mobile elements flank blocks of functionally related sequences, thereby suggesting a role in transposition. The plasmid contains 12 reiterated DNA families that are likely to participate in genomic rearrangements. Comparisons between this plasmid and complete rhizobial genomes and symbiotic compartments already sequenced show a general lack of synteny and colinearity, with the exception of some transcriptional units. There are only 20 symbiotic genes that are shared by all SGCs. CONCLUSIONS: Our data support the notion that the symbiotic compartments of rhizobia genomes are mosaic structures that have been frequently tailored by recombination, horizontal transfer and transposition.