Regulatory potential, phyletic distribution and evolution of ancient, intracellular small-molecule-binding domains

CPred: a web server for predicting viable circular permutations in proteins

Circular permutation (CP) is a protein structural rearrangement phenomenon, through which nature allows structural homologs to have different locations of termini and thus varied activities, stabilities and functional properties. It can be applied in many fields of protein research and bioengineering. The limitation of applying CP lies in its technical complexity, high cost and uncertainty of the viability of the resulting protein variants. Not every position in a protein can be used to create a viable circular permutant, but there is still a lack of practical computational tools for evaluating the positional feasibility of CP before costly experiments are carried out. We have previously designed a comprehensive method for predicting viable CP cleavage sites in proteins. In this work, we implement that method into an efficient and user-friendly web server named CPred (CP site predictor), which is supposed to be helpful to promote fundamental researches and biotechnological applications of CP. The CPred is accessible at http://sarst.life.nthu.edu.tw/CPred.

Crystal structure of tandem ACT domain-containing protein ACTP from Galdieria sulphuraria

The ACT domain is a structurally conserved small molecule binding domain which is mostly involved in amino acid and purine metabolism. Here, we report the crystal structure of a tandem ACT domain-containing protein (ACTP) from Galdieria sulphuraria. The two ACTP monomers in the asymmetric unit form a dimer with a non-crystallographic twofold axis in a domain-swapped manner, showing a horseshoe-like structure with a central crevice. This structure contributes to expand our knowledge on the structural diversity of ACT domain-containing proteins.

Study and reengineering of the binding sites and allosteric regulation of biosynthetic threonine deaminase by isoleucine and valine in Escherichia coli

Biosynthetic threonine deaminase (TD) is a key enzyme for the synthesis of isoleucine which is allosterically inhibited and activated by Ile and Val, respectively. The binding sites of Ile and Val and the mechanism of their regulations in TD are not clear, but essential for a rational design of efficient productive strain(s) for Ile and related amino acids. In this study, structure-based computational approach and site-directed mutagenesis were combined to identify the potential binding sites of Ile and Val in Escherichia coli TD. Our results demonstrated that each regulatory domain of the TD monomer possesses two nonequivalent effector-binding sites. The residues R362, E442, G445, A446, Y369, I460, and S461 only interact with Ile while E347, G350, and F352 are involved not only in the Ile binding but also in the Val binding. By further considering enzyme kinetic data, we propose a concentration-dependent mechanism of the allosteric regulation of TD by Ile and Val. For the construction of Ile overproducing strain, a novel TD mutant with double mutation of F352A/R362F was also created, which showed both higher activity and much stronger resistance to Ile inhibition comparing to those of wild-type enzyme. Overexpression of this mutant TD in E. coli JW3591 significantly increased the production of ketobutyrate and Ile in comparison to the reference strains overexpressing wild-type TD or the catabolic threonine deaminase (TdcB). This work builds a solid basis for the reengineering of TD and related microorganisms for Ile production.

Novel bacterial gas sensor proteins with transition metal-containing prosthetic groups as active sites

SIGNIFICANCE

Gas molecules function as signaling molecules in many biological regulatory systems responsible for transcription, chemotaxis, and other complex physiological processes. Gas sensor proteins play a crucial role in regulating such biological systems in response to gas molecules.

RECENT ADVANCES

New sensor proteins that sense oxygen or nitric oxide have recently been found, and they have been characterized by X-ray crystallographic and/or spectroscopic analysis. It has become clear that the interaction between a prosthetic group and gas molecules triggers dynamic structural changes in the protein backbone when a gas sensor protein senses gas molecules. Gas sensor proteins employ novel mechanisms to trigger conformational changes in the presence of a gas.

CRITICAL ISSUES

In gas sensor proteins that have iron-sulfur clusters as active sites, the iron-sulfur clusters undergo structural changes, which trigger a conformational change. Heme-based gas sensor proteins reconstruct hydrogen-bonding networks around the heme and heme-bound ligand.

FUTURE DIRECTION

Gas sensor proteins have two functional states, on and off, which are active and inactive, respectively, for subsequent signal transduction in response to their physiological effector molecules. To fully understand the structure-function relationships of gas sensor proteins, it is vital to perform X-ray crystal structure analyses of full-length proteins in both the on and off states.

Identification of multicomponent histidine-aspartate phosphorelay system controlling flagellar and motility gene expression in Geobacter species

Geobacter species play an important role in the natural biogeochemical cycles of aquatic sediments and subsurface environments as well as in subsurface bioremediation by oxidizing organic compounds with the reduction of insoluble Fe(III) oxides. Flagellum-based motility is considered to be critical for Geobacter species to locate fresh sources of Fe(III) oxides. Functional and comparative genomic approaches, coupled with genetic and biochemical methods, identified key regulators for flagellar gene expression in Geobacter species. A master transcriptional regulator, designated FgrM, is a member of the enhancer-binding protein family. The fgrM gene in the most studied strain of Geobacter species, Geobacter sulfurreducens strain DL-1, is truncated by a transposase gene, preventing flagellar biosynthesis. Integrating a functional FgrM homolog restored flagellar biosynthesis and motility in G. sulfurreducens DL-1 and enhanced the ability to reduce insoluble Fe(III) oxide. Interrupting the fgrM gene in G. sulfurreducens strain KN400, which is motile, removed the capacity for flagellar production and inhibited Fe(III) oxide reduction. FgrM, which is also a response regulator of the two-component His-Asp phosphorelay system, was phosphorylated by histidine kinase GHK4, which was essential for flagellar production and motility. GHK4, which is a hybrid kinase with a receiver domain at the N terminus, was phosphorylated by another histidine kinase, GHK3. Therefore, the multicomponent His-Asp phosphorelay system appears to control flagellar gene expression in Geobacter species.

Deciphering the preference and predicting the viability of circular permutations in proteins

Circular permutation (CP) refers to situations in which the termini of a protein are relocated to other positions in the structure. CP occurs naturally and has been artificially created to study protein function, stability and folding. Recently CP is increasingly applied to engineer enzyme structure and function, and to create bifunctional fusion proteins unachievable by tandem fusion. CP is a complicated and expensive technique. An intrinsic difficulty in its application lies in the fact that not every position in a protein is amenable for creating a viable permutant. To examine the preferences of CP and develop CP viability prediction methods, we carried out comprehensive analyses of the sequence, structural, and dynamical properties of known CP sites using a variety of statistics and simulation methods, such as the bootstrap aggregating, permutation test and molecular dynamics simulations. CP particularly favors Gly, Pro, Asp and Asn. Positions preferred by CP lie within coils, loops, turns, and at residues that are exposed to solvent, weakly hydrogen-bonded, environmentally unpacked, or flexible. Disfavored positions include Cys, bulky hydrophobic residues, and residues located within helices or near the protein's core. These results fostered the development of an effective viable CP site prediction system, which combined four machine learning methods, e.g., artificial neural networks, the support vector machine, a random forest, and a hierarchical feature integration procedure developed in this work. As assessed by using the hydrofolate reductase dataset as the independent evaluation dataset, this prediction system achieved an AUC of 0.9. Large-scale predictions have been performed for nine thousand representative protein structures; several new potential applications of CP were thus identified. Many unreported preferences of CP are revealed in this study. The developed system is the best CP viability prediction method currently available. This work will facilitate the application of CP in research and biotechnology.

Ryanodine receptor calcium release channels: an evolutionary perspective

Ryanodine receptors (RyRs), along with the related inositol 1,4,5-trisphosphate receptors (IP(3)Rs), mediate the release of Ca(2+) from intracellular organelles of eukaryotes. As discussed in other chapters, such increases in intracellular Ca(2+) levels act a fundamental second messenger, regulating a diverse array of cellular processes. For over two decades, it has been reported that vertebrates express multiple RYR genes, whereas non-vertebrate multicellular organisms possess a single homologue within their genomes. Recently, the existence of RyR-like channels in unicellular organisms has also been reported. This chapter exploits recent expansions in available genome data to generate an overview of the expression of RyR-like genes in organisms representing a broad range of viral, archaeal, bacterial and eukaryotic taxa. Analyses of the multidomain structures and phylogenetic relationships of these proteins has lead to a model in which, early during eukaryotic evolution, IP(3)R-like ancestral Ca(2+) release channels were converted to RyR proteins via the addition of promiscuous protein domains, possibly via horizontal gene transfer mechanisms.

Structure and evolution of ubiquitin and ubiquitin-related domains

Since its discovery over three decades ago, it has become abundantly clear that the ubiquitin (Ub) system is a quintessential feature of all aspects of eukaryotic biology. At the heart of the system lies the conjugation and deconjugation of Ub and Ub-like (Ubls) proteins to proteins or lipids drastically altering the biochemistry of the targeted molecules. In particular, it represents the primary mechanism by which protein stability is regulated in eukaryotes. Ub/Ubls are typified by the β-grasp fold (β-GF) that has additionally been recruited for a strikingly diverse range of biochemical functions. These include catalytic roles (e.g., NUDIX phosphohydrolases), scaffolding of iron-sulfur clusters, binding of RNA and other biomolecules such as co-factors, sulfur transfer in biosynthesis of diverse metabolites, and as mediators of key protein-protein interactions in practically every conceivable cellular context. In this chapter, we present a synthetic overview of the structure, evolution, and natural classification of Ub, Ubls, and other members of the β-GF. The β-GF appears to have differentiated into at least seven clades by the time of the last universal common ancestor of all extant organisms, encompassing much of the structural diversity observed in extant versions. The β-GF appears to have first emerged in the context of translation-related RNA-interactions and subsequently exploded to occupy various functional niches. Most biochemical diversification of the fold occurred in prokaryotes, with the eukaryotic phase of its evolution mainly marked by the expansion of the Ubl clade of the β-GF. Consequently, at least 70 distinct Ubl families are distributed across eukaryotes, of which nearly 20 families were already present in the eukaryotic common ancestor. These included multiple protein and one lipid conjugated forms and versions that functions as adapter domains in multimodule polypeptides. The early diversification of the Ubl families in eukaryotes played a major role in the emergence of characteristic eukaryotic cellular substructures and systems pertaining to nucleo-cytoplasmic compartmentalization, vesicular trafficking, lysosomal targeting, protein processing in the endoplasmic reticulum, and chromatin dynamics. Recent results from comparative genomics indicate that precursors of the eukaryotic Ub-system were already present in prokaryotes. The most basic versions are those combining an Ubl and an E1-like enzyme involved in metabolic pathways related to metallopterin, thiamine, cysteine, siderophore and perhaps modified base biosynthesis. Some of these versions also appear to have given rise to simple protein-tagging systems such as Sampylation in archaea and Urmylation in eukaryotes. However, other prokaryotic systems with Ubls of the YukD and other families, including one very close to Ub itself, developed additional elements that more closely resemble the eukaryotic state in possessing an E2, a RING-type E3, or both of these components. Additionally, prokaryotes have evolved conjugation systems that are independent of Ub ligases, such as the Pup system.

The natural history of ubiquitin and ubiquitin-related domains

The ubiquitin (Ub) system is centered on conjugation and deconjugation of Ub and Ub-like (Ubls) proteins by a system of ligases and peptidases, respectively. Ub/Ubls contain the beta-grasp fold, also found in numerous proteins with biochemically distinct roles unrelated to the conventional Ub-system. The beta-GF underwent an early radiation spawning at least seven clades prior to the divergence of extant organisms from their last universal common ancestor, first emerging in the context of translation-related RNA-interactions and subsequently exploding to occupy various functional niches. Most beta-GF diversification occurred in prokaryotes, with the Ubl clade showing dramatic expansion in the eukaryotes. Diversification of Ubl families in eukaryotes played a major role in emergence of characteristic eukaryotic cellular sub-structures and systems. Recent comparative genomics studies indicate precursors of the eukaryotic Ub-system emerged in prokaryotes. The simplest of these combine an Ubl and an E1-like enzyme in metabolic pathways. Sampylation in archaea and Urmylation in eukaryotes appear to represent recruitment of such systems as simple protein-tagging apparatuses. However, other prokaryotic systems incorporated further components and mirror the eukaryotic condition in possessing an E2, a RING-type E3 or both of these components. Additionally, prokaryotes have evolved conjugation systems independent of Ub ligases, such as the Pup system.

Insights from the architecture of the bacterial transcription apparatus

We provide a portrait of the bacterial transcription apparatus in light of the data emerging from structural studies, sequence analysis and comparative genomics to bring out important but underappreciated features. We first describe the key structural highlights and evolutionary implications emerging from comparison of the cellular RNA polymerase subunits with the RNA-dependent RNA polymerase involved in RNAi in eukaryotes and their homologs from newly identified bacterial selfish elements. We describe some previously unnoticed domains and the possible evolutionary stages leading to the RNA polymerases of extant life forms. We then present the case for the ancient orthology of the basal transcription factors, the sigma factor and TFIIB, in the bacterial and the archaeo-eukaryotic lineages. We also present a synopsis of the structural and architectural taxonomy of specific transcription factors and their genome-scale demography. In this context, we present certain notable deviations from the otherwise invariant proteome-wide trends in transcription factor distribution and use it to predict the presence of an unusual lineage-specifically expanded signaling system in certain firmicutes like Paenibacillus. We then discuss the intersection between functional properties of transcription factors and the organization of transcriptional networks. Finally, we present some of the interesting evolutionary conundrums posed by our newly gained understanding of the bacterial transcription apparatus and potential areas for future explorations.

Joint scaling laws in functional and evolutionary categories in prokaryotic genomes

We propose and study a class-expansion/innovation/loss model of genome evolution taking into account biological roles of genes and their constituent domains. In our model, numbers of genes in different functional categories are coupled to each other. For example, an increase in the number of metabolic enzymes in a genome is usually accompanied by addition of new transcription factors regulating these enzymes. Such coupling can be thought of as a proportional ‘recipe’ for genome composition of the type ‘a spoonful of sugar for each egg yolk’. The model jointly reproduces two known empirical laws: the distribution of family sizes and the non-linear scaling of the number of genes in certain functional categories (e.g. transcription factors) with genome size. In addition, it allows us to derive a novel relation between the exponents characterizing these two scaling laws, establishing a direct quantitative connection between evolutionary and functional categories. It predicts that functional categories that grow faster-than-linearly with genome size to be characterized by flatter-than-average family size distributions. This relation is confirmed by our bioinformatics analysis of prokaryotic genomes. This proves that the joint quantitative trends of functional and evolutionary classes can be understood in terms of evolutionary growth with proportional recipes.

The role of the GAF and central domains of the transcriptional activator VnfA in Azotobacter vinelandii

VnfA is a transcriptional activator that is required for the expression of the structural genes encoding nitrogenase-2 in Azotobacter vinelandii. VnfA consists of three domains: an N-terminal regulatory domain termed GAF, including a Cys-rich motif; a central domain from the AAA+ family; and a C-terminal domain for DNA binding. Previously, we reported that transcriptionally active VnfA harboring an Fe-S cluster (presumably of the 3Fe-4S type) as a prosthetic group and the Cys-rich motif were possibly associated with coordination of the Fe-S cluster. In the present study, we have investigated the roles of the GAF and central domains in the regulatory function of VnfA using truncated variants: ΔN15(VnfA) and ΔGAF(VnfA) that lack the N-terminal 15 residues and whole GAF domain, respectively, and GAF(VnfA) consisting of only the GAF domain. ΔN15(VnfA) and ΔGAF(VnfA) lost the ability to bind the Fe-S cluster, whereas GAF(VnfA) was still able to bind to the cluster, consistent with the hypothesis that the Cys-rich motif is essential for Fe-S cluster binding. The GAF domain showed an inhibitory effect on the transcriptional activity of VnfA, which was reversed in the presence of the Fe-S cluster, and reactivated upon disassembly of the cluster. The inhibitory activity of the GAF domain acts on the NTPase activity of the central domain, whereas the binding ability of VnfA to DNA was not significantly affected, when VnfA retains its tetrameric conformation. The results imply that a major pathway, by which VnfA function is regulated, operates via the control of NTPase activity by the GAF domain.

Comparative analyses imply that the enigmatic Sigma factor 54 is a central controller of the bacterial exterior

Background

Sigma-54 is a central regulator in many pathogenic bacteria and has been linked to a multitude of cellular processes like nitrogen assimilation and important functional traits such as motility, virulence, and biofilm formation. Until now it has remained obscure whether these phenomena and the control by Sigma-54 share an underlying theme.

Results

We have uncovered the commonality by performing a range of comparative genome analyses. A) The presence of Sigma-54 and its associated activators was determined for all sequenced prokaryotes. We observed a phylum-dependent distribution that is suggestive of an evolutionary relationship between Sigma-54 and lipopolysaccharide and flagellar biosynthesis. B) All Sigma-54 activators were identified and annotated. The relation with phosphotransfer-mediated signaling (TCS and PTS) and the transport and assimilation of carboxylates and nitrogen containing metabolites was substantiated. C) The function annotations, that were represented within the genomic context of all genes encoding Sigma-54, its activators and its promoters, were analyzed for intra-phylum representation and inter-phylum conservation. Promoters were localized using a straightforward scoring strategy that was formulated to identify similar motifs. We found clear highly-represented and conserved genetic associations with genes that concern the transport and biosynthesis of the metabolic intermediates of exopolysaccharides, flagella, lipids, lipopolysaccharides, lipoproteins and peptidoglycan.

Conclusion

Our analyses directly implicate Sigma-54 as a central player in the control over the processes that involve the physical interaction of an organism with its environment like in the colonization of a host (virulence) or the formation of biofilm.

Transcriptome analysis shows differential gene expression in the saprotrophic to parasitic transition of Pochonia chlamydosporia

Expression profiles were identified in the fungus Pochonia chlamydosporia, a biological control agent of plant parasitic nematodes, through a cDNA-amplified fragment length polymorphism approach. Two isolates with different host ranges, IMI 380407 and IMI 331547, were assayed in conditions of saprotrophic-to-parasitic transition, through in vitro assays. Gene expression profiles from three different nutritional conditions and four sampling times were generated, with eggs of host nematodes Globodera pallida and Meloidogyne incognita. Expression of transcripts changed in RNA fingerprints obtained under different nutritional stresses (starvation in presence/absence of eggs, or rich growth media). Transcript derived fragments (TDFs) obtained from the expression profiles corresponded to 6,800 products. A subset was sequenced and their expression profile confirmed through RT PCR. A total of 57 TDFs were selected for further analysis, based on similarities to translated or annotated sequences. Genes expressed during egg parasitism for both IMI 380407 and IMI 331547 were involved in metabolic functions, cellular signal regulation, cellular transport, regulation of gene expression, DNA repair, and other unknown functions. Multivariate analysis of TDF expression showed three groups for IMI 380407 and one for IMI 331547, each characterized by expression of genes related to eggs parasitism. Common amplification profiles among TDF clusters from both isolates also reflected a pool of constitutive genes, not affected by the nutritional conditions and nematode associations, related to general metabolic functions. The differential expression of parasitism related genes suggest a network of induced/repressed products, playing a role in fungal signaling and infection, with partial overlaps in host infection and parasitism traits.

Light-induced transcriptional responses associated with proteorhodopsin-enhanced growth in a marine flavobacterium

Proteorhodopsin (PR) is a photoprotein that functions as a light-driven proton pump in diverse marine Bacteria and Archaea. Recent studies have suggested that PR may enhance both growth rate and yield in some flavobacteria when grown under nutrient-limiting conditions in the light. The direct involvement of PR, and the metabolic details enabling light-stimulated growth, however, remain uncertain. Here, we surveyed transcriptional and growth responses of a PR-containing marine flavobacterium during carbon-limited growth in the light and the dark. As previously reported (Gómez-Consarnau et al., 2007), Dokdonia strain MED134 exhibited light-enhanced growth rates and cell yields under low carbon growth conditions. Inhibition of retinal biosynthesis abolished the light-stimulated growth response, supporting a direct role for retinal-bound PR in light-enhanced growth. Among protein-coding transcripts, both PR and retinal biosynthetic enzymes showed significant upregulation in the light. Other light-associated proteins, including bacterial cryptochrome and DNA photolyase, were also expressed at significantly higher levels in the light. Membrane transporters for Na(+)/phosphate and Na(+)/alanine symporters, and the Na(+)-translocating NADH-quinone oxidoreductase (NQR) linked electron transport chain, were also significantly upregulated in the light. Culture experiments using a specific inhibitor of Na(+)-translocating NQR indicated that sodium pumping via NQR is a critical metabolic process in the light-stimulated growth of MED134. In total, the results suggested the importance of both the PR-enabled, light-driven proton gradient, as well as the generation of a Na(+) ion gradient, as essential components for light-enhanced growth in these flavobacteria.

Mammalian Cyclic Nucleotide Phosphodiesterases: Molecular Mechanisms and Physiological Functions

The superfamily of cyclic nucleotide (cN) phosphodiesterases (PDEs) is comprised of 11 families of enzymes. PDEs break down cAMP and/or cGMP and are major determinants of cellular cN levels and, consequently, the actions of cN-signaling pathways. PDEs exhibit a range of catalytic efficiencies for breakdown of cAMP and/or cGMP and are regulated by myriad processes including phosphorylation, cN binding to allosteric GAF domains, changes in expression levels, interaction with regulatory or anchoring proteins, and reversible translocation among subcellular compartments. Selective PDE inhibitors are currently in clinical use for treatment of erectile dysfunction, pulmonary hypertension, intermittent claudication, and chronic pulmonary obstructive disease; many new inhibitors are being developed for treatment of these and other maladies. Recently reported x-ray crystallographic structures have defined features that provide for specificity for cAMP or cGMP in PDE catalytic sites or their GAF domains, as well as mechanisms involved in catalysis, oligomerization, autoinhibition, and interactions with inhibitors. In addition, major advances have been made in understanding the physiological impact and the biochemical basis for selective localization and/or recruitment of specific PDE isoenzymes to particular subcellular compartments. The many recent advances in understanding PDE structures, functions, and physiological actions are discussed in this review.

Identification of succinate exporter in Corynebacterium glutamicum and its physiological roles under anaerobic conditions

Corynebacterium glutamicum produces succinate from glucose via the reductive tricarboxylic acid cycle under microaerobic and anaerobic conditions. We identified a NCgl2130 gene of C. glutamicum as a novel succinate exporter that functions in succinate production, and designated sucE1. sucE1 expression levels were higher under microaerobic conditions than aerobic conditions, and overexpression or disruption of sucE1 respectively increased or decreased succinate productivity during fermentation. Under microaerobic conditions, the sucE1 disruptant sucE1Δ showed 30% less succinate productivity and a lower sugar-consumption rate than the parental strain. Under anaerobic conditions, succinate production by sucE1Δ ceased. The intracellular succinate and fructose-1,6-bisphosphate levels of sucE1Δ under microaerobic conditions were respectively 1.7-fold and 1.6-fold higher than those of the parental strain, suggesting that loss of SucE1 function caused a failure of succinate removal from the cells, leading to intracellular accumulation that inhibited upstream sugar metabolism. Homology and transmembrane helix searches identified SucE1 as a membrane protein belonging to the aspartate:alanine exchanger (AAE) family. Partially purified 6x-histidine-tagged SucE1 (SucE1-[His](6)) reconstituted in succinate-loaded liposomes clearly demonstrated counterflow and self-exchange activities for succinate. Together, these findings suggest that sucE1 encodes a novel succinate exporter that is induced under microaerobic conditions, and is important for succinate production under both microaerobic and anaerobic conditions.

Natural history of the eukaryotic chromatin protein methylation system

In eukaryotes, methylation of nucleosomal histones and other nuclear proteins is a central aspect of chromatin structure and dynamics. The past 15 years have seen an enormous advance in our understanding of the biochemistry of these modifications, and of their role in establishing the epigenetic code. We provide a synthetic overview, from an evolutionary perspective, of the main players in the eukaryotic chromatin protein methylation system, with an emphasis on catalytic domains. Several components of the eukaryotic protein methylation system had their origins in bacteria. In particular, the Rossmann fold protein methylases (PRMTs and DOT1), and the LSD1 and jumonji-related demethylases and oxidases, appear to have emerged in the context of bacterial peptide methylation and hydroxylation systems. These systems were originally involved in synthesis of peptide secondary metabolites, such as antibiotics, toxins, and siderophores. The peptidylarginine deiminases appear to have been acquired by animals from bacterial enzymes that modify cell-surface proteins. SET domain methylases, which display the β-clip fold, apparently first emerged in prokaryotes from the SAF superfamily of carbohydrate-binding domains. However, even in bacteria, a subset of the SET domains might have evolved a chromatin-related role in conjunction with a BAF60a/b-like SWIB domain protein and topoisomerases. By the time of the last eukaryotic common ancestor, multiple SET and PRMT methylases were already in place and are likely to have mediated methylation at the H3K4, H3K9, H3K36, and H4K20 positions, and carried out both asymmetric and symmetric arginine dimethylation. Inference of H3K27 methylation in the ancestral eukaryote appears uncertain, though it was certainly in place a little later in eukaryotic evolution. Current data suggest that unlike SET methylases, which are universally present in eukaryotes, demethylases are not. They appear to be absent in the earliest-branching eukaryotic lineages, and emerged later along with several other chromatin proteins, such as the Dot1-methylase, prior to divergence of the kinetoplastid-heterolobosean lineage from the remaining eukaryotes. This period also corresponds to the point of origin of DNA cytosine methylation by DNMT1. Origin of major lineages of SET domains such as the Trithorax, Su(var)3-9, Ash1, SMYD, and TTLL12 and E(Z) might have played the initial role in the establishment of multiple distinct heterochromatic and euchromatic states that are likely to have been present, in some form, through much of eukaryotic evolution. Elaboration of these chromatin states might have gone hand-in-hand with acquisition of multiple jumonji-related and LSD1-like demethylases, and functional linkages with the DNA methylation and RNAi systems. Throughout eukaryotic evolution, there were several lineage-specific expansions of SET domain proteins, which might be related to a special transcription regulation process in trypanosomes, acquisition of new meiotic recombination hotspots in animals, and methylation and associated modifications of the diatom silaffin proteins involved in silica biomineralization. The use of specific domains to "read" the methylation marks appears to have been present in the ancestral eukaryote itself. Of these the chromo-like domains appear to have been acquired from bacterial secreted proteins that might have a role in binding cell-surface peptides or peptidoglycan. Domain architectures of the primary enzymes involved in the eukaryotic protein methylation system indicate key features relating to interactions with each other and other modifications in chromatin, such as acetylation. They also emphasize the profound functional distinction between the role of demethylation and deacetylation in regulation of chromatin dynamics.

An overview of prokaryotic transcription factors : a summary of function and occurrence in bacterial genomes

Transcriptional initiation is arguably the most important control point for gene expression. It is regulated by a combination of factors, including DNA sequence and its three-dimensional topology, proteins and small molecules. In this chapter, we focus on the trans-acting factors of bacterial regulation. Initiation begins with the recruitment of the RNA polymerase holoenzyme to a specific locus upstream of the gene known as its promoter. The sigma factor, which is a component of the holoenzyme, provides the most fundamental mechanisms for orchestrating broad changes in gene expression state. It is responsible for promoter recognition as well as recruiting the holoenzyme to the promoter. Distinct sigma factors compete with for binding to a common pool of RNA polymerases, thus achieving condition-dependent differential expression. Another important class of bacterial regulators is transcription factors, which activate or repress transcription of target genes typically in response to an environmental or cellular trigger. These factors may be global or local depending on the number of genes and range of cellular functions that they target. The activities of both global and local transcription factors may be regulated either at a post-transcriptional level via signal-sensing protein domains or at the level of their own expression. In addition to modulating polymerase recruitment to promoters, several global factors are considered as "nucleoid-associated proteins" that impose structural constraints on the chromosome by altering the conformation of the bound DNA, thus influencing other processes involving DNA such as replication and recombination. This chapter concludes with a discussion of how regulatory interactions between transcription factors and their target genes can be represented as a network.

Ligand-binding PAS domains in a genomic, cellular, and structural context

Per-Arnt-Sim (PAS) domains occur in proteins from all kingdoms of life. In the bacterial kingdom, PAS domains are commonly positioned at the amino terminus of signaling proteins such as sensor histidine kinases, cyclic-di-GMP synthases/hydrolases, and methyl-accepting chemotaxis proteins. Although these domains are highly divergent at the primary sequence level, the structures of dozens of PAS domains across a broad section of sequence space have been solved, revealing a conserved three-dimensional architecture. An all-versus-all alignment of 63 PAS structures demonstrates that the PAS domain family forms structural clades on the basis of two principal variables: (a) topological location inside or outside the plasma membrane and (b) the class of small molecule that they bind. The binding of a chemically diverse range of small-molecule metabolites is a hallmark of the PAS domain family. PAS ligand binding either functions as a primary cue to initiate a cellular signaling response or provides the domain with the capacity to respond to secondary physical or chemical signals such as gas molecules, redox potential, or photons. This review synthesizes the current state of knowledge of the structural foundations and evolution of ligand recognition and binding by PAS domains.

Phylogenetic analysis and evolution of aromatic amino acid hydroxylase

A study was performed to investigate the phylogenetic relationship among AAAH members and to statistically evaluate sequence conservation and functional divergence. In total, 161 genes were identified from 103 species. Phylogenetic analysis showed that well-conserved subfamilies exist. Exon-intron structure analysis showed that the gene structures of AAAH were highly conserved across some different lineage species, while some species-specific introns were also found. The dynamic distribution of ACT domain suggested one gene fusion event has occurred in eukaryota. Significant functional divergence was found between some subgroups. Analysis of the site-specific profiles revealed critical amino acid residues for functional divergence. This study highlights the molecular evolution of this family and may provide a starting point for further experimental verifications.

A Rhodobacter capsulatus member of a universal permease family imports molybdate and other oxyanions

Molybdenum (Mo) is an important trace element that is toxic at high concentrations. To resolve the mechanisms underlying Mo toxicity, Rhodobacter capsulatus mutants tolerant to high Mo concentrations were isolated by random transposon Tn5 mutagenesis. The insertion sites of six independent isolates mapped within the same gene predicted to code for a permease of unknown function located in the cytoplasmic membrane. During growth under Mo-replete conditions, the wild-type strain accumulated considerably more Mo than the permease mutant. For mutants defective for the permease, the high-affinity molybdate importer ModABC, or both transporters, in vivo Mo-dependent nitrogenase (Mo-nitrogenase) activities at different Mo concentrations suggested that ModABC and the permease import molybdate in nanomolar and micromolar ranges, respectively. Like the permease mutants, a mutant defective for ATP sulfurylase tolerated high Mo concentrations, suggesting that ATP sulfurylase is the main target of Mo inhibition in R. capsulatus. Sulfate-dependent growth of a double mutant defective for the permease and the high-affinity sulfate importer CysTWA was reduced compared to those of the single mutants, implying that the permease plays an important role in sulfate uptake. In addition, permease mutants tolerated higher tungstate and vanadate concentrations than the wild type, suggesting that the permease acts as a general oxyanion importer. We propose to call this permease PerO (for oxyanion permease). It is the first reported bacterial molybdate transporter outside the ABC transporter family.

The structure of SSO2064, the first representative of Pfam family PF01796, reveals a novel two-domain zinc-ribbon OB-fold architecture with a potential acyl-CoA-binding role

SSO2064 is the first structural representative of PF01796 (DUF35), a large prokaryotic family with a wide phylogenetic distribution. The structure reveals a novel two-domain architecture comprising an N-terminal, rubredoxin-like, zinc ribbon and a C-terminal, oligonucleotide/oligosaccharide-binding (OB) fold domain. Additional N-terminal helical segments may be involved in protein-protein interactions. Domain architectures, genomic context analysis and functional evidence from certain bacterial representatives of this family suggest that these proteins form a novel fatty-acid-binding component that is involved in the biosynthesis of lipids and polyketide antibiotics and that they possibly function as acyl-CoA-binding proteins. This structure has led to a re-evaluation of the DUF35 family, which has now been split into two entries in the latest Pfam release (v.24.0).

The structure of Jann_2411 (DUF1470) from Jannaschia sp. at 1.45 Å resolution reveals a new fold (the ABATE domain) and suggests its possible role as a transcription regulator

The crystal structure of Jann_2411 from Jannaschia sp. strain CCS1, a member of the Pfam PF07336 family classified as a domain of unknown function (DUF1470), was solved to a resolution of 1.45 Å by multiple-wavelength anomalous dispersion (MAD). This protein is the first structural representative of the DUF1470 Pfam family. Structural analysis revealed a two-domain organization, with the N-terminal domain presenting a new fold called the ABATE domain that may bind an as yet unknown ligand. The C-terminal domain forms a treble-clef zinc finger that is likely to be involved in DNA binding. Analysis of the Jann_2411 protein and the broader ABATE-domain family suggests a role as stress-induced transcriptional regulators.

Identification of the membrane protein SucE and its role in succinate transport in Corynebacterium glutamicum

Succinic acid is excreted during anaerobiosis by many bacteria, and manifold applications are known making the biotechnological production of succinate attractive. Although the pathways for succinate formation are known, succinate export is not understood in most of the succinate producing bacteria. Here, we present a bioinformatic approach for identification of a putative succinate export system in Corynebacterium glutamicum. The subsequent screening revealed that a mutant in the gene cg2425 is impaired in succinate production or transport under anaerobic conditions. A function of the Cg2425 protein as import system was excluded. In contrast, a role of the Cg2425 protein as succinate export system was indicated by accumulation of increased amounts of internal succinate under anaerobic conditions in a Cg2425-dependent manner and a concomitant impairment of external succinate accumulation. In conclusion, we propose that Cg2425 participates in succinate export in C. glutamicum and suggest the name SucE for the protein.

Cooperative and uncooperative cyclic-nucleotide-gated ion channels

Ion channels gated by cyclic nucleotides serve multiple functions in sensory signaling in diverse cell types ranging from neurons to sperm. Newly discovered members from bacteria and marine invertebrates provide a wealth of structural and functional information on this channel family. A hallmark of classical tetrameric cyclic-nucleotide-gated channels is their cooperative activation by binding of several ligands. By contrast, the new members seem to be uncooperative, and binding of a single ligand molecule suffices to open these channels. These new findings provide a fresh look at the mechanism of allosteric activation of ion channels.

The nature of protein domain evolution: shaping the interaction network

The proteomes that make up the collection of proteins in contemporary organisms evolved through recombination and duplication of a limited set of domains. These protein domains are essentially the main components of globular proteins and are the most principal level at which protein function and protein interactions can be understood. An important aspect of domain evolution is their atomic structure and biochemical function, which are both specified by the information in the amino acid sequence. Changes in this information may bring about new folds, functions and protein architectures. With the present and still increasing wealth of sequences and annotation data brought about by genomics, new evolutionary relationships are constantly being revealed, unknown structures modeled and phylogenies inferred. Such investigations not only help predict the function of newly discovered proteins, but also assist in mapping unforeseen pathways of evolution and reveal crucial, co-evolving inter- and intra-molecular interactions. In turn this will help us describe how protein domains shaped cellular interaction networks and the dynamics with which they are regulated in the cell. Additionally, these studies can be used for the design of new and optimized protein domains for therapy. In this review, we aim to describe the basic concepts of protein domain evolution and illustrate recent developments in molecular evolution that have provided valuable new insights in the field of comparative genomics and protein interaction networks.

Novel regulatory cascades controlling expression of nitrogen-fixation genes in Geobacter sulfurreducens

Geobacter species often play an important role in bioremediation of environments contaminated with metals or organics and show promise for harvesting electricity from waste organic matter in microbial fuel cells. The ability of Geobacter species to fix atmospheric nitrogen is an important metabolic feature for these applications. We identified novel regulatory cascades controlling nitrogen-fixation gene expression in Geobacter sulfurreducens. Unlike the regulatory mechanisms known in other nitrogen-fixing microorganisms, nitrogen-fixation gene regulation in G. sulfurreducens is controlled by two two-component His–Asp phosphorelay systems. One of these systems appears to be the master regulatory system that activates transcription of the majority of nitrogen-fixation genes and represses a gene encoding glutamate dehydrogenase during nitrogen fixation. The other system whose expression is directly activated by the master regulatory system appears to control by antitermination the expression of a subset of the nitrogen-fixation genes whose transcription is activated by the master regulatory system and whose promoter contains transcription termination signals. This study provides a new paradigm for nitrogen-fixation gene regulation.

Simplified models of biological networks

The function of living cells is controlled by complex regulatory networks that are built of a wide diversity of interacting molecular components. The sheer size and intricacy of molecular networks of even the simplest organisms are obstacles toward understanding network functionality. This review discusses the achievements and promise of a bottom-up approach that uses well-characterized subnetworks as model systems for understanding larger networks. It highlights the interplay between the structure, logic, and function of various types of small regulatory circuits. The bottom-up approach advocates understanding regulatory networks as a collection of entangled motifs. We therefore emphasize the potential of negative and positive feedback, as well as their combinations, to generate robust homeostasis, epigenetics, and oscillations.

Addition at the molecular level: signal integration in designed Per-ARNT-Sim receptor proteins

Survival of organisms in dynamic environments requires accurate perception and integration of signals. At the molecular level, signal detection is mediated by signal receptor proteins that largely are of modular composition. Sensor modules, such as the widespread Per-ARNT-Sim (PAS) domains, detect signals and, in response, regulate the biological activity of effector modules. Here, we exploit the modularity of signal receptors to design and engineer synthetic receptors that comprise two PAS sensor domains responsive to different signals, and we use these signals to control the activity of a histidine kinase effector. Designed two-input PAS receptors detected oxygen and blue light in a positive cooperative manner. The extent of the response to the signals was dictated by domain topology: the dominant regulatory effect was exerted by the PAS domain proximal to the effector domain. The presence of one sensor domain modulated the signal response function of the other. Sequence and structural data on natural receptors with tandem PAS domains show that these are predominantly linked by short amphipathic alpha-helices. Signals from multiple sensor domains could be integrated and propagated to the effector domain as torques. Our results inform the rational design of receptors that integrate multiple signals to modulate cellular behavior.

Origin and evolution of peptide-modifying dioxygenases and identification of the wybutosine hydroxylase/hydroperoxidase

Unlike classical 2-oxoglutarate and iron-dependent dioxygenases, which include several nucleic acid modifiers, the structurally similar jumonji-related dioxygenase superfamily was only known to catalyze peptide modifications. Using comparative genomics methods, we predict that a family of jumonji-related enzymes catalyzes wybutosine hydroxylation/peroxidation at position 37 of eukaryotic tRNAPhe. Identification of this enzyme raised questions regarding the emergence of protein- and nucleic acid-modifying activities among jumonji-related domains. We addressed these with a natural classification of DSBH domains and reconstructed the precursor of the dioxygenases as a sugar-binding domain. This precursor gave rise to sugar epimerases and metal-binding sugar isomerases. The sugar isomerase active site was exapted for catalysis of oxygenation, with a radiation of these enzymes in bacteria, probably due to impetus from the primary oxygenation event in Earth’s history. 2-Oxoglutarate-dependent versions appear to have further expanded with rise of the tricarboxylic acid cycle. We identify previously under-appreciated aspects of their active site and multiple independent innovations of 2-oxoacid-binding basic residues among these superfamilies. We show that double-stranded β-helix dioxygenases diversified extensively in biosynthesis and modification of halogenated siderophores, antibiotics, peptide secondary metabolites and glycine-rich collagen-like proteins in bacteria. Jumonji-related domains diversified into three distinct lineages in bacterial secondary metabolism systems and these were precursors of the three major clades of eukaryotic enzymes. The specificity of wybutosine hydroxylase/peroxidase probably relates to the structural similarity of the modified moiety to the ancestral amino acid substrate of this superfamily.

Origin and diversification of basic-helix-loop-helix proteins in plants

Basic helix-loop-helix (bHLH) proteins are a class of transcription factors found throughout eukaryotic organisms. Classification of the complete sets of bHLH proteins in the sequenced genomes of Arabidopsis thaliana and Oryza sativa (rice) has defined the diversity of these proteins among flowering plants. However, the evolutionary relationships of different plant bHLH groups and the diversity of bHLH proteins in more ancestral groups of plants are currently unknown. In this study, we use whole-genome sequences from nine species of land plants and algae to define the relationships between these proteins in plants. We show that few (less than 5) bHLH proteins are encoded in the genomes of chlorophytes and red algae. In contrast, many bHLH proteins (100-170) are encoded in the genomes of land plants (embryophytes). Phylogenetic analyses suggest that plant bHLH proteins are monophyletic and constitute 26 subfamilies. Twenty of these subfamilies existed in the common ancestors of extant mosses and vascular plants, whereas six further subfamilies evolved among the vascular plants. In addition to the conserved bHLH domains, most subfamilies are characterized by the presence of highly conserved short amino acid motifs. We conclude that much of the diversity of plant bHLH proteins was established in early land plants, over 440 million years ago.

Cyclic nucleotide binding GAF domains from phosphodiesterases: structural and mechanistic insights

GAF domains regulate the catalytic activity of certain vertebrate cyclic nucleotide phosphodiesterases (PDEs) by allosteric, noncatalytic binding of cyclic nucleotides. GAF domains arranged in tandem are found in PDE2, -5, -6, -10, and -11, all of which regulate the cellular concentrations of the second messengers cAMP and/or cGMP. Nucleotide binding to GAF domains affects the overall conformation and the catalytic activity of full-length PDEs. The cyclic nucleotide-bound GAF domains from PDE2, -5, -6, and -10 all adopt a conserved fold but show subtle differences within the binding pocket architecture that account for a large range of nucleotide affinities and selectivity. NMR data and details from the structure of full-length nucleotide-free PDE2A reveal the dynamic nature and magnitude of the conformational change that accompanies nucleotide binding. The discussed GAF domain structures further reveal differences in dimerization properties and highlight the structural diversity within GAF domain-containing PDEs.

Exploring the molecular basis for selective binding of Mycobacterium tuberculosis Asp kinase toward its natural substrates and feedback inhibitors: a docking and molecular dynamics study

Tuberculosis (TB) is still a major public health problem, compounded by the human immunodeficiency virus (HIV)-TB co-infection and recent emergence of multidrug-resistant (MDR) and extensively drug resistant (XDR)-TB. In this context, aspartokinase of mycobacterium tuberculosis has drawn attention for designing novel anti-TB drugs. Asp kinase is an enzyme responsible for the synthesis of 4-phospho-L-aspartate from L-aspartate and involved in the branched biosynthetic pathway leading to the synthesis of amino acids lysine, threonine, methionine and isoleucine. An intermediate of lysine biosynthetic branch, mesodiaminopimelate is also a component of the peptidoglycan which is a component of bacterial cell wall. To interfere with the production of all these amino acids and cell wall, it is possible to inhibit Asp kinase activity. This can be achieved using Asp kinase inhibitors. In order to design novel Asp kinase inhibitors as effective anti-TB drugs, it is necessary to have an understanding of the binding sites of Asp kinase. As no crystal structure of the enzyme has yet been published, we built a homology model of Asp kinase using the crystallized Asp kinase from M. Jannaschii, as template structures (2HMF and 3C1M). After the molecular dynamics refinement, the optimized homology model was assessed as a reliable structure by PROCHECK, ERRAT, WHAT-IF, PROSA2003 and VERIFY-3D. The results of molecular docking studies with natural substrates, products and feedback inhibitors are in agreement with the published data and showed that ACT domain plays an important role in binding to ligands. Based on the docking conformations, pharmacophore model can be developed by probing the common features of ligands. By analyzing the results, ACT domain architecture, certain key residues that are responsible for binding to feedback inhibitors and natural substrates were identified. This would be very helpful in understanding the blockade mechanism of Asp kinase and providing insights into rational design of novel Asp kinase inhibitors for M.tuberculosis.

Structure and signaling mechanism of Per-ARNT-Sim domains

Per-ARNT-Sim (PAS) domains serve as versatile sensor and interaction modules in signal transduction proteins. PAS sensors detect chemical and physical stimuli and regulate the activity of functionally diverse effector domains. In contrast to this chemical, physical, and functional diversity, the structure of the core of PAS domains is broadly conserved and comprises a five-stranded antiparallel beta sheet and several alpha helices. Signals originate within the conserved core and generate structural and dynamic changes predominantly within the beta sheet, from which they propagate via amphipathic alpha-helical and coiled-coil linkers at the N or C termini of the core to the covalently attached effector domain. Effector domains are typically dimeric; their activity appears to be largely regulated by signal-dependent changes in quaternary structure and dynamics. The signaling mechanisms of PAS and other signaling domains share common features, and these commonalities can be exploited to enable structure-based design of artificial photosensors and chemosensors.

Structure and function of plant photoreceptors

Signaling photoreceptors use the information contained in the absorption of a photon to modulate biological activity in plants and a wide range of organisms. The fundamental-and as yet imperfectly answered-question is, how is this achieved at the molecular level? We adopt the perspective of biophysicists interested in light-dependent signal transduction in nature and the three-dimensional structures that underpin signaling. Six classes of photoreceptors are known: light-oxygen-voltage (LOV) sensors, xanthopsins, phytochromes, blue-light sensors using flavin adenine dinucleotide (BLUF), cryptochromes, and rhodopsins. All are water-soluble proteins except rhodopsins, which are integral membrane proteins; all are based on a modular architecture except cryptochromes and rhodopsins; and each displays a distinct, light-dependent chemical process based on the photochemistry of their nonprotein chromophore, such as isomerization about a double bond (xanthopsins, phytochromes, and rhodopsins), formation or rupture of a covalent bond (LOV sensors), or electron transfer (BLUF sensors and cryptochromes).

A novel hybrid kinase is essential for regulating the sigma(B)-mediated stress response of Bacillus cereus

A common bacterial strategy for monitoring environmental challenges is to use two-component systems, which consist of a sensor histidine kinase (HK) and a response regulator (RR). In the food-borne pathogen Bacillus cereus, the alternative sigma factor sigma(B) is activated by the RR RsbY. Here we present strong indications that the PP2C-type phosphatase RsbY receives its input from the multi-sensor hybrid kinase BC1008 (renamed RsbK). Genome analyses revealed that, across bacilli, rsbY and rsbK are located in a conserved gene cluster. A B. cereus rsbK deletion strain was shown to be incapable of inducing sigma(B) upon stress conditions and was impaired in its heat adaptive response. Comparison of the wild-type and rsbK mutant transcriptomes upon heat shock revealed that RsbK was primarily involved in the activation of the sigma(B)-mediated stress response. Truncation of the RsbK RR receiver domain demonstrated the importance of this domain for sigma(B) induction upon stress. The domain architecture of RsbK suggests that in the B. cereus group and in other bacilli, environmental and intracellular stress signalling routes are combined into one single protein. This strategy is markedly different from the sigma(B) activation pathway in other low-GC Gram-positives.

Molecular evolution of multisubunit RNA polymerases: structural analysis

Comprehensive multiple sequence alignments of the multisubunit DNA-dependent RNA polymerase (RNAP) large subunits, including the bacterial beta and beta' subunits and their homologs from archaebacterial RNAPs, eukaryotic RNAPs I-III, nuclear-cytoplasmic large double-stranded DNA virus RNAPs, and plant plastid RNAPs, were created [Lane, W. J. and Darst, S. A. (2009). Molecular evolution of multisubunit RNA polymerases: sequence analysis. In press]. The alignments were used to delineate sequence regions shared among all classes of multisubunit RNAPs, defining common, fundamental RNAP features as well as identifying highly conserved positions. Here, we present a systematic, detailed structural analysis of these shared regions and highly conserved positions in terms of the RNAP structure, as well as the RNAP structure/function relationship, when known.

Analysis of two gene clusters involved in the degradation of 4-fluorophenol by Arthrobacter sp. strain IF1

Arthrobacter sp. strain IF1 is able to grow on 4-fluorophenol (4-FP) as a sole source of carbon and energy. To clone the 4-FP degradation genes, DNA libraries were constructed and screened with a probe obtained by PCR using primers designed on the basis of conserved regions of aromatic two-component monooxygenases. Sequencing of positive clones yielded two gene clusters, each harboring a gene encoding a monooxygenase with high sequence similarity to the oxygenase component of 4-nitrophenol and 4-chlorophenol monooxygenase systems. Both these monooxygenase genes were differentially expressed during growth on 4-FP, as revealed by Northern blotting and reverse transcription-PCR. One cluster also contained a gene for a flavin reductase. The monooxygenase and reductase were purified from Escherichia coli cells expressing the corresponding genes, and together they catalyzed NADH-dependent hydroxylation and dehalogenation of 4-halophenols. The results indicate that strain IF1 transforms 4-FP to hydroquinone by a two-component monooxygenase system of which one component provides reduced flavin adenine dinucleotide at the expense of NADH and the other catalyzes para-hydroxylation of 4-FP and other 4-substituted phenols.

Receiver domains control the active-state stoichiometry of Aquifex aeolicus sigma54 activator NtrC4, as revealed by electrospray ionization mass spectrometry

A common challenge with studies of proteins in vitro is determining which constructs and conditions are most physiologically relevant. sigma(54) activators are proteins that undergo regulated assembly to form an active ATPase ring that enables transcription by sigma(54)-polymerase. Previous studies of AAA(+) ATPase domains from sigma(54) activators have shown that some are heptamers, while others are hexamers. Because active oligomers assemble from off-state dimers, it was thought that even-numbered oligomers should dominate, and that heptamer formation would occur when individual domains of the activators, rather than the intact proteins, were studied. Here we present results from electrospray ionization mass spectrometry experiments characterizing the assembly states of intact NtrC4 (a sigma(54) activator from Aquifex aeolicus, an extreme thermophile), as well as its ATPase domain alone, and regulatory-ATPase and ATPase-DNA binding domain combinations. We show that the full-length and activated regulatory-ATPase proteins form hexamers, whereas the isolated ATPase domain, unactivated regulatory-ATPase, and ATPase-DNA binding domain form heptamers. Activation of the N-terminal regulatory domain is the key factor stabilizing the hexamer form of the ATPase, relative to the heptamer.

Toolbox model of evolution of prokaryotic metabolic networks and their regulation

It has been reported that the number of transcription factors encoded in prokaryotic genomes scales approximately quadratically with their total number of genes. We propose a conceptual explanation of this finding and illustrate it using a simple model in which metabolic and regulatory networks of prokaryotes are shaped by horizontal gene transfer of coregulated metabolic pathways. Adapting to a new environmental condition monitored by a new transcription factor (e.g., learning to use another nutrient) involves both acquiring new enzymes and reusing some of the enzymes already encoded in the genome. As the repertoire of enzymes of an organism (its toolbox) grows larger, it can reuse its enzyme tools more often and thus needs to get fewer new ones to master each new task. From this observation, it logically follows that the number of functional tasks and their regulators increases faster than linearly with the total number of genes encoding enzymes. Genomes can also shrink, e.g., because of a loss of a nutrient from the environment, followed by deletion of its regulator and all enzymes that become redundant. We propose several simple models of network evolution elaborating on this toolbox argument and reproducing the empirically observed quadratic scaling. The distribution of lengths of pathway branches in our model agrees with that of the real-life metabolic network of Escherichia coli. Thus, our model provides a qualitative explanation for broad distributions of regulon sizes in prokaryotes.

iSARST: an integrated SARST web server for rapid protein structural similarity searches

iSARST is a web server for efficient protein structural similarity searches. It is a multi-processor, batch-processing and integrated implementation of several structural comparison tools and two database searching methods: SARST for common structural homologs and CPSARST for homologs with circular permutations. iSARST allows users submitting multiple PDB/SCOP entry IDs or an archive file containing many structures. After scanning the target database using SARST/CPSARST, the ordering of hits are refined with conventional structure alignment tools such as FAST, TM-align and SAMO, which are run in a PC cluster. In this way, iSARST achieves a high running speed while preserving the high precision of refinement engines. The final outputs include tables listing co-linear or circularly permuted homologs of the query proteins and a functional summary of the best hits. Superimposed structures can be examined through an interactive and informative visualization tool. iSARST provides the first batch mode structural comparison web service for both co-linear homologs and circular permutants. It can serve as a rapid annotation system for functionally unknown or hypothetical proteins, which are increasing rapidly in this post-genomics era. The server can be accessed at http://sarst.life.nthu.edu.tw/iSARST/.

Genomic and structural aspects of protein evolution

It has been known for more than 35 years that, during evolution, new proteins are formed by gene duplications, sequence and structural divergence and, in many cases, gene combinations. The genome projects have produced complete, or almost complete, descriptions of the protein repertoires of over 600 distinct organisms. Analyses of these data have dramatically increased our understanding of the formation of new proteins. At the present time, we can accurately trace the evolutionary relationships of about half the proteins found in most genomes, and it is these proteins that we discuss in the present review. Usually, the units of evolution are protein domains that are duplicated, diverge and form combinations. Small proteins contain one domain, and large proteins contain combinations of two or more domains. Domains descended from a common ancestor are clustered into superfamilies. In most genomes, the net growth of superfamily members means that more than 90% of domains are duplicates. In a section on domain duplications, we discuss the number of currently known superfamilies, their size and distribution, and superfamily expansions related to biological complexity and to specific lineages. In a section on divergence, we describe how sequences and structures diverge, the changes in stability produced by acceptable mutations, and the nature of functional divergence and selection. In a section on domain combinations, we discuss their general nature, the sequential order of domains, how combinations modify function, and the extraordinary variety of the domain combinations found in different genomes. We conclude with a brief note on other forms of protein evolution and speculations of the origins of the duplication, divergence and combination processes.

Capturing the VirA/VirG TCS of Agrobacterium tumefaciens

Two-component systems (TCS) regulate pathogenic commitment in many interactions and provide an opportunity for unique therapeutic intervention. The VirA/VirG TCS of Agrobacterium tumefaciens mediates inter-kingdom gene transfer in the development of host tumors and sets in motion the events that underlie the great success of this multi-host plant pathogen. Significant proof for the feasibility of interventions has now emerged with the discovery of a natural product that effectively "blinds" the pathogen to the host via inhibition of VirA/VirG signal transduction. Moreover, the emerging studies on the mechanism of signal perception have revealed general sites suitable for intervention of TCS signaling. Given the extensive functional homology, it should now be possible to transfer the models discovered for VirA/VirG broadly to other pathogenic interactions.

A genomic analysis of the archaeal system Ignicoccus hospitalis-Nanoarchaeum equitans

Sequencing of the complete genome of Ignicoccus hospitalis gives insight into its association with another species of Archaea, Nanoarchaeum equitans.

Background

The relationship between the hyperthermophiles Ignicoccus hospitalis and Nanoarchaeum equitans is the only known example of a specific association between two species of Archaea. Little is known about the mechanisms that enable this relationship.

Results

We sequenced the complete genome of I. hospitalis and found it to be the smallest among independent, free-living organisms. A comparative genomic reconstruction suggests that the I. hospitalis lineage has lost most of the genes associated with a heterotrophic metabolism that is characteristic of most of the Crenarchaeota. A streamlined genome is also suggested by a low frequency of paralogs and fragmentation of many operons. However, this process appears to be partially balanced by lateral gene transfer from archaeal and bacterial sources.

Conclusions

A combination of genomic and cellular features suggests highly efficient adaptation to the low energy yield of sulfur-hydrogen respiration and efficient inorganic carbon and nitrogen assimilation. Evidence of lateral gene exchange between N. equitans and I. hospitalis indicates that the relationship has impacted both genomes. This association is the simplest symbiotic system known to date and a unique model for studying mechanisms of interspecific relationships at the genomic and metabolic levels.

Role of cAMP in mediating AHR signaling

Regulation of the nuclear import of many transcription factors represents a step in gene regulation which is crucial for a number of cellular processes. The aryl hydrocarbon receptor (AHR), a basic helix-loop-helix protein of the PAS (PER-ARNT-SIM) family of transcriptional regulators is a cytosol-associated and ligand-activated receptor. The environmental toxin dioxin binds with high affinity to AHR rendering it nuclear and leading to the activation of AHR sensitive genes. However, the fact, that the AHR mediates a large variety of physiological events without the involvement of any known exogenous ligand, including liver and vascular system development, maturation of the immune system, regulation of genes involved in cellular growth, cell differentiation and circadian rhythm, speaks for an important role of AHR in cell biology independent of the presence of an exogenous ligand. Different approaches were applied to study mechanism(s) which render AHR nuclear and design its function in absence of exogenous ligands. We found that AHR is sensitive to cAMP signaling mediated by cAMP-dependent protein kinase (PKA) which fundamentally differs from AHR signaling mediated by the exogenous ligand dioxin. It has been shown that PKA mediated signaling can be confined by compartmentalization of signaling components in microdomains conferring specificity to signaling by the ubiquitous second messenger cAMP. Moreover, A-kinase-anchoring proteins (AKAPs) and newly discovered cAMP receptors, Epac (exchange protein directly activated by cAMP), may give us a further chance to enter into new dimensions of cAMP signal transmissions that potentially may bring us closer to AHR physiology.

CPDB: a database of circular permutation in proteins

Circular permutation (CP) in a protein can be considered as if its sequence were circularized followed by a creation of termini at a new location. Since the first observation of CP in 1979, a substantial number of studies have concluded that circular permutants (CPs) usually retain native structures and functions, sometimes with increased stability or functional diversity. Although this interesting property has made CP useful in many protein engineering and folding researches, large-scale collections of CP-related information were not available until this study. Here we describe CPDB, the first CP DataBase. The organizational principle of CPDB is a hierarchical categorization in which pairs of circular permutants are grouped into CP clusters, which are further grouped into folds and in turn classes. Additions to CPDB include a useful set of tools and resources for the identification, characterization, comparison and visualization of CP. Besides, several viable CP site prediction methods are implemented and assessed in CPDB. This database can be useful in protein folding and evolution studies, the discovery of novel protein structural and functional relationships, and facilitating the production of new CPs with unique biotechnical or industrial interests. The CPDB database can be accessed at http://sarst.life.nthu.edu.tw/cpdb