Phylogenetic relationships among bacteriorhodopsins

Structural and Functional Studies of a Newly Grouped Haloquadratum walsbyi Bacteriorhodopsin Reveal the Acid-resistant Light-driven Proton Pumping Activity

Retinal bound light-driven proton pumps are widespread in eukaryotic and prokaryotic organisms. Among these pumps, bacteriorhodopsin (BR) proteins cooperate with ATP synthase to convert captured solar energy into a biologically consumable form, ATP. In an acidic environment or when pumped-out protons accumulate in the extracellular region, the maximum absorbance of BR proteins shifts markedly to the longer wavelengths. These conditions affect the light-driven proton pumping functional exertion as well. In this study, wild-type crystal structure of a BR with optical stability under wide pH range from a square halophilic archaeon, Haloquadratum walsbyi (HwBR), was solved in two crystal forms. One crystal form, refined to 1.85 Å resolution, contains a trimer in the asymmetric unit, whereas another contains an antiparallel dimer was refined at 2.58 Å. HwBR could not be classified into any existing subgroup of archaeal BR proteins based on the protein sequence phylogenetic tree, and it showed unique absorption spectral stability when exposed to low pH values. All structures showed a unique hydrogen-bonding network between Arg⁸² and Thr²⁰¹, linking the BC and FG loops to shield the retinal-binding pocket in the interior from the extracellular environment. This result was supported by R82E mutation that attenuated the optical stability. The negatively charged cytoplasmic side and the Arg⁸²–Thr²⁰¹ hydrogen bond may play an important role in the proton translocation trend in HwBR under acidic conditions. Our findings have unveiled a strategy adopted by BR proteins to solidify their defenses against unfavorable environments and maintain their optical properties associated with proton pumping.

Bioinformatic characterization of the trimeric intracellular cation-specific channel protein family

Trimeric intracellular cation-specific (TRIC) channels are integral to muscle excitation-contraction coupling. TRIC channels provide counter-ionic flux when calcium is rapidly transported from intracellular stores to the cell cytoplasm. Until recently, knowledge of the presence of these proteins was limited to animals. We analyzed the TRIC family and identified a profusion of prokaryotic family members with topologies and motifs similar to those of their eukaryotic counterparts. Prokaryotic members far outnumber eukaryotic members, and although none has been functionally characterized, the evidence suggests that they function as secondary carriers. The presence of fused N- or C-terminal domains of known biochemical functions as well as genomic context analyses provide clues about the functions of these prokaryotic homologs. They are proposed to function in metabolite (e.g., amino acid/nucleotide) efflux. Phylogenetic analysis revealed that TRIC channel homologs diverged relatively early during evolutionary history and that horizontal gene transfer was frequent in prokaryotes but not in eukaryotes. Topological analyses of TRIC channels revealed that these proteins possess seven putative transmembrane segments (TMSs), which arose by intragenic duplication of a three-TMS polypeptide-encoding genetic element followed by addition of a seventh TMS at the C terminus to give the precursor of all current TRIC family homologs. We propose that this family arose in prokaryotes.

Homologues of archaeal rhodopsins in plants, animals and fungi: structural and functional predications for a putative fungal chaperone protein

The microbial rhodopsins (MR) are homologous to putative chaperone and retinal-binding proteins of fungi. These proteins comprise a coherent family that we have termed the MR family. We have used modeling techniques to predict the structure of one of the putative yeast chaperone proteins, YRO2, based on homology with bacteriorhodopsins (BR). Availability of the structure allowed depiction of conserved residues that are likely to be of functional significance. The results lead us to predict an extracellular protein folding function and a transmembrane proton transport pathway. We suggest that protein folding is energized by a novel mechanism involving the proton motive force. We further show that MR family proteins are distantly related to a family of fungal, animal and plant proteins that include the human lysosomal cystine transporter (LCT) of man (cystinosin), mutations in which cause cystinosis. Sequence and phylogenetic analyses of both the MR family and the LCT family are reported. Proteins in both families are of the same approximate size, exhibit seven putative transmembrane alpha-helical spanners (TMSs) and show limited sequence similarity. We show that the LCT family arose by an internal gene duplication event and that TMSs 1-3 are homologous to TMSs 5-7. Although the same could not be demonstrated statistically for MR family members, homology with the LCT family suggests (but does not prove) a common evolutionary pathway. Thus, TMSs 1-3 and 5-7 in both LCT and MR family members may share a common origin, accounting for their shared structural features.

Phylogenetic approaches to the identification and characterization of protein families and superfamilies

With the advent of megabase genome sequencing, the need for computational analyses increases exponentially. Sequencing errors must be corrected, encoded proteins must be identified, functions must be assigned to these proteins, and distant phylogenetic relationships must be recognized in order to maximize the yield of information obtainable from genome sequencing projects. Both the computer and the human brain have their limitations, but using them in combination, the biologist can vastly extend his or her analytic capabilities. Computer techniques can be used to estimate protein structure, function, biogenesis, and evolution. In this review, the application of available computer programs to several protein families, particularly transport, receptor, and transcriptional regulatory protein families, illustrate our current capabilities and limitations. Although some multidomain protein families are evolutionarily homogeneous, others have mosaic origins. Evidence concerning the nature and frequency of occurrence of domain shuffling, splicing, fusion, deletion, and duplication during evolution of specific protein families is evaluated. It is shown that specific families of enzymes, receptors, transport proteins, and transcriptional regulatory proteins share a common evolutionary origin, frequently diverging in function because of domain splicing and ligation. Some large families arose gradually over evolutionary time, whereas others developed suddenly, due to bursts of intragenic or intergenic (or both) duplication events occurring over relatively short periods of time. It is argued that energy coupling to transport was a late occurrence, superimposed on preexisting mechanisms of solute facilitation. It is also shown that several transport protein families have evolved independently of each other, employing different routes, at different times in evolutionary history, to give topologically similar transmembrane protein complexes.

A functional-phylogenetic classification system for transmembrane solute transporters

A comprehensive classification system for transmembrane molecular transporters has been developed and recently approved by the transport panel of the nomenclature committee of the International Union of Biochemistry and Molecular Biology. This system is based on (i) transporter class and subclass (mode of transport and energy coupling mechanism), (ii) protein phylogenetic family and subfamily, and (iii) substrate specificity. Almost all of the more than 250 identified families of transporters include members that function exclusively in transport. Channels (115 families), secondary active transporters (uniporters, symporters, and antiporters) (78 families), primary active transporters (23 families), group translocators (6 families), and transport proteins of ill-defined function or of unknown mechanism (51 families) constitute distinct categories. Transport mode and energy coupling prove to be relatively immutable characteristics and therefore provide primary bases for classification. Phylogenetic grouping reflects structure, function, mechanism, and often substrate specificity and therefore provides a reliable secondary basis for classification. Substrate specificity and polarity of transport prove to be more readily altered during evolutionary history and therefore provide a tertiary basis for classification. With very few exceptions, a phylogenetic family of transporters includes members that function by a single transport mode and energy coupling mechanism, although a variety of substrates may be transported, sometimes with either inwardly or outwardly directed polarity. In this review, I provide cross-referencing of well-characterized constituent transporters according to (i) transport mode, (ii) energy coupling mechanism, (iii) phylogenetic grouping, and (iv) substrates transported. The structural features and distribution of recognized family members throughout the living world are also evaluated. The tabulations should facilitate familial and functional assignments of newly sequenced transport proteins that will result from future genome sequencing projects.

Evolution of the archaeal rhodopsins: evolution rate changes by gene duplication and functional differentiation

The amino acid sequences of 25 archaeal retinal proteins from 13 different strains of extreme halophiles were analyzed to establish their molecular phylogenetic relationship. On the basis of amino acid sequence similarity, these proteins apparently formed a distinct family designated as the archaeal rhodopsin family (ARF), which was not related to other known proteins, including G protein-coupled receptors. The archaeal rhodopsin family was further divided into four clusters with different functions; H+ pump (bacteriorhodopsin), Cl- pump (halorhodopsin), and two kinds of sensor (sensory rhodopsin and phoborhodopsin). These four rhodopsin clusters seemed to have occurred by gene duplication(s) before the generic speciation of halophilic archaea, based on phylogenetic analysis. Therefore, the degrees of differences in amino acid sequences within each cluster simply reflected the divergent evolution of halophilic archaea. By comparing the branch lengths after speciation points of the reconstituted tree, we calculated the relative evolution rates of the four archaeal rhodopsins bacteriorhodopsin:halorhodopsin:sensory rhodopsin: phoborhodopsin to be 5:4:3:10. From these values, the degrees of functional and structural restriction of each protein can be inferred. The branching topology of four clusters grouped bacteriorhodopsin and halorhodopsin versus sensory rhodopsin and phoborhodopsin by likelihood mapping. Using bacteriorhodopsin (and halorhodopsin) as an outgroup, the gene duplication point of sensory rhodopsin/phoborhodopsin was determined. By calculating the branch lengths between the gene duplication point and each halophilic archaea speciation point, we could speculate upon the relative evolution rate of pre-sensory rhodopsin and pre-phoborhodopsin. The evolution rate of pre-sensory rhodopsin was fivefold faster than that of pre-phoborhodopsin, which suggests that the original function of the ancestral sensor was similar to that of phoborhodopsin, and that sensory rhodopsin evolved from pre-sensory rhodopsin by the accumulation of mutations. The changes in evolution rate by gene duplication and functional differentiation were demonstrated in the archaeal rhodopsin family using the gene duplication date and halobacterial speciation date as common time stamps.

An analysis of the conserved residues between halobacterial retinal proteins and G-protein coupled receptors: implications for GPCR modeling

An alignment of the transmembrane domains of halobacterial retinal proteins (including bacteriorhodopsin) and G-protein coupled receptors (GPCRs) is presented based on the commonality of conserved residues between families. Due to the limited sequence homology displayed by these proteins, an alternative strategy is proposed for sequence alignment that correlates residues within secondary structure elements. The nonsequential alignment developed identifies three proline and two aspartates residues that share common positions and, in the former case, similar functions in the transmembrane domain. The alignment is further applied to model the packing of transmembrane helices 5 and 6 of the beta-adrenergic receptor based on the backbone coordinates of bacteriorhodopsin helices 3 and 2, respectively. Unlike models derived from standard sequential alignments, the approach developed here allows the key structural features conferred by the proline residues to be captured during model building. The structure described is also compared with available site directed mutagenesis results as well as existing GPCR models. In addition to the implications to model building, the commonality observed suggests a potential relationship among the GPCRs and retinal proteins.

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

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

A simple flexible program for the computational analysis of amino acyl residue distribution in proteins: application to the distribution of aromatic versus aliphatic hydrophobic amino acids in transmembrane alpha-helical spanners of integral membrane transport proteins

We describe a simple, flexible program (AAD) with a primary function of depicting the distribution of aliphatic and aromatic amino acid residues along the linear aligned sequence of a family of homologous proteins and a secondary function of depicting the distribution of all amino acids along the same linear sequence. The program is used to examine the distribution of aromatic versus aliphatic residues in representative well-characterized families of polytopic membrane proteins. Many but not all such protein families are shown to exhibit a predominance of aliphatic residues in the central regions of their transmembrane spanners but a predominance of aromatic residues at the peripheries of their spanners. We propose that this distribution stabilizes the hydrophobic-hydrophilic interface and renders the centers of these integral membrane proteins more fluid than their peripheries.

Topological and segmental phylogenetic analyses of the anion exchanger (band 3) family of transporters

Eleven sequenced anion exchanger (AE; band 3) proteins, including five AE1, four AE2 and two AE3 proteins, comprise the anion exchanger family (AEF) of homologous proteins. Eliminating the rat and rabbit proteins that are nearly identical to the corresponding mouse proteins, seven dissimilar members of this family were selected for study, divided into N-terminal, central and C-terminal segments (designated segments 0, 1 and 2, respectively) and analysed separately for sequence similarity and phylogenetic relatedness. Segments 0 are variable in length and sequence, are essentially lacking in some of the members of the AEF, and are not demonstrably homologous in other members of the family. All segments 1 and 2 are homologous, but they exhibit widely differing degrees of sequence divergence. Segments 2 are highly conserved in all AEF proteins. Segments 1 of the AE2 and AE3 proteins are as conserved as are segments 2, but segments 1 of the AE1 proteins have diverged from each other and from the AE2 and AE3 segments 1 much more than have segments 2 of these same proteins. The distributions of various types of amino acid residues in the putative transmembrane helical spanners of the seven dissimilar members of the AEF, based on a modification of the 14-spanner model of Wang et al. (1994) was determined, and this distribution was compared with those of other transmembrane transport proteins of known structure (bacterial rhodopsins, outer membrane porins of Gram-negative bacteria and bacterial photosynthetic reaction centres.(ABSTRACT TRUNCATED AT 250 WORDS)

Differential codon usage: a safeguard against inappropriate expression of specialized genes?

Recent work has suggested that rare codons are sometimes used for the regulation of specialized gene expression in bacteria. Moreover, the cellular levels of certain tRNAs may fluctuate with growth conditions. Evidence implicating such mechanisms in the control of photosynthesis in Rhodobacter, solventogenesis in Clostridium, sporulation in Streptomyces, and fimbrial phase variation in E. coli is summarized. It is suggested that such mechanisms will prove applicable to the control of numerous additional specialized functions, and that the empirical tools for testing this possibility are currently available.