Sentra: a database of signal transduction proteins for comparative genome analysis

PA0575 (RmcA) interacts with other c-di-GMP metabolizing proteins in Pseudomonas aeruginosa PAO1

As a central signaling molecule, c-di-GMP (bis-(3,5)-cyclic diguanosine monophosphate) is becoming the focus for research in bacteria physiology. Pseudomonas aeruginosa PAO1 genome contains highly complicated c-di-GMP metabolizing genes and a number of these proteins have been identified and investigated. Especially, a sophisticated network of these proteins is emerging. In current study, mainly through Bacteria-2-Hybrid assay, we found PA0575 (RmcA), a GGDEF-EAL dual protein, to interact with two other dual proteins of PA4601 (MorA) and PA4959 (FimX). These observations imply the intricacy of c-di-GMP metabolizing protein interactions. Our work thus provides one piece of data to increase the understandings to c-di-GMP signaling.

Enhancing multiplex genome editing by natural transformation (MuGENT) via inactivation of ssDNA exonucleases

Recently, we described a method for multiplex genome editing by natural transformation (MuGENT). Mutant constructs for MuGENT require large arms of homology (>2000 bp) surrounding each genome edit, which necessitates laborious in vitro DNA splicing. In Vibrio cholerae, we uncover that this requirement is due to cytoplasmic ssDNA exonucleases, which inhibit natural transformation. In ssDNA exonuclease mutants, one arm of homology can be reduced to as little as 40 bp while still promoting integration of genome edits at rates of ∼50% without selection in cis. Consequently, editing constructs are generated in a single polymerase chain reaction where one homology arm is oligonucleotide encoded. To further enhance editing efficiencies, we also developed a strain for transient inactivation of the mismatch repair system. As a proof-of-concept, we used these advances to rapidly mutate 10 high-affinity binding sites for the nucleoid occlusion protein SlmA and generated a duodecuple mutant of 12 diguanylate cyclases in V. cholerae. Whole genome sequencing revealed little to no off-target mutations in these strains. Finally, we show that ssDNA exonucleases inhibit natural transformation in Acinetobacter baylyi. Thus, rational removal of ssDNA exonucleases may be broadly applicable for enhancing the efficacy and ease of MuGENT in diverse naturally transformable species.

Different evolutionary modifications as a guide to rewire two-component systems

Two-component systems (TCS) are short signalling pathways generally occurring in prokaryotes. They frequently regulate prokaryotic stimulus responses and thus are also of interest for engineering in biotechnology and synthetic biology. The aim of this study is to better understand and describe rewiring of TCS while investigating different evolutionary scenarios. Based on large-scale screens of TCS in different organisms, this study gives detailed data, concrete alignments, and structure analysis on three general modification scenarios, where TCS were rewired for new responses and functions: (i) exchanges in the sequence within single TCS domains, (ii) exchange of whole TCS domains; (iii) addition of new components modulating TCS function. As a result, the replacement of stimulus and promotor cassettes to rewire TCS is well defined exploiting the alignments given here. The diverged TCS examples are non-trivial and the design is challenging. Designed connector proteins may also be useful to modify TCS in selected cases.

Extending signaling pathways with protein-interaction networks. Application to apoptosis

Cells exploit signaling pathways during responses to environmental changes, and these processes are often modulated during disease. Particularly, relevant human pathologies such as cancer or viral infections require downregulating apoptosis signaling pathways to progress. As a result, the identification of proteins responsible for these changes is essential for the diagnostics and development of therapeutics. Transferring functional annotation within protein interaction networks has proven useful to identify such proteins, although this is not a trivial task. Here, we used different scoring methods to transfer annotation from 53 well-studied members of the human apoptosis pathways (as known by 2005) to their protein interactors. All scoring methods produced significant predictions (compared to a random negative model), but its number was too large to be useful. Thus, we made a final prediction using specific combinations of scoring methods and compared it to the proteins related to apoptosis signaling pathways during the last 5 years. We propose 273 candidate proteins that may be relevant in apoptosis signaling pathways. Although some of them have known functions consistent with their proposed apoptotsis involvement, the majority have not been annotated yet, leaving room for further experimental studies. We provide our predictions at http://sbi.imim.es/web/Apoptosis.php.

Pathogenic and mutualistic plant-bacteria interactions: ever increasing similarities

Plant-interacting bacteria can establish either mutualistic or pathogenic interactions that cause beneficial or detrimental effects respectively, to their hosts. In spite of the completely different outcomes, accumulating evidence indicates that similar molecular bases underlie the establishment of these two contrasting plant-bacteria associations. Recent findings observed in the mutualistic nitrogen-fixing Rhizobium-legume symbiosis add new elements to the increasing list of similarities. Amongst these, in this review we describe the role of plant resistance proteins in determining host specificity in the Rhizobium-legume symbiosis that resemble the gene-for-gene resistance of plant-pathogen interactions, and the production of antimicrobial peptides by certain legumes to control rhizobial proliferation within nodules. Amongst common bacterial strategies, cyclic diguanylate (c-di-GMP) appears to be a second messenger used by both pathogenic and mutualistic bacteria to regulate key features for interaction with their plant hosts.

A multi-repeat adhesin of the phytopathogen, Pectobacterium atrosepticum, is secreted by a Type I pathway and is subject to complex regulation involving a non-canonical diguanylate cyclase

Cyclic diguanylate (c-di-GMP) is a second messenger controlling many important bacterial processes. The phytopathogen Pectobacterium atrosepticum SCRI1043 (Pba1043) possesses a Type I secretion system (T1SS) essential for the secretion of a proteinaceous multi-repeat adhesin (MRP) required for binding to the host plant. The genes encoding the MRP and the T1SS are tightly linked to genes encoding several putative c-di-GMP regulatory components. We show that c-di-GMP regulates secreted MRP levels in Pba1043 through the action of two genes encoding predicted diguanylate cyclase (DGC) and phosphodiesterase proteins (ECA3270 and ECA3271). Phenotypic analyses and quantification of c-di-GMP levels demonstrated that ECA3270 and ECA3271 regulate secreted MRP levels by increasing and decreasing, respectively, the intracellular levels of c-di-GMP. Moreover, ECA3270 represents the first active DGC reported to have an alternative active-site motif from the 'canonical' GG[D/E]EF. ECA3270 has an A-site motif of SGDEF and analysis of single amino acid replacements demonstrated that the first position of this motif can tolerate functional substitution. Serine in position one of the A-site is also observed in many other DGCs. Finally, another T1SS-linked regulator (ECA3265) also plays an important role in regulating secreted MRP, with an altered localization of MRP observed in an ECA3265 mutant background. Mutants defective in these three T1SS-linked regulators exhibit a reduction in root binding and virulence, confirming that this complex, finely tuned regulation system is crucial in the interaction with host plants.

Evolution and applications of plant pathway resources and databases

Plants are important sources of food and plant products are essential for modern human life. Plants are increasingly gaining importance as drug and fuel resources, bioremediation tools and as tools for recombinant technology. Considering these applications, database infrastructure for plant model systems deserves much more attention. Study of plant biological pathways, the interconnection between these pathways and plant systems biology on the whole has in general lagged behind human systems biology. In this article we review plant pathway databases and the resources that are currently available. We lay out trends and challenges in the ongoing efforts to integrate plant pathway databases and the applications of database integration. We also discuss how progress in non-plant communities can serve as an example for the improvement of the plant pathway database landscape and thereby allow quantitative modeling of plant biosystems. We propose Good Database Practice as a possible model for collaboration and to ease future integration efforts.

P2CS: a database of prokaryotic two-component systems

P2CS (http://www.p2cs.org) is a specialized database for prokaryotic two-component systems (TCSs), virtually ubiquitous signalling proteins which regulate a wide range of physiological processes. The primary aim of the database is to annotate and classify TCS proteins from completely sequenced prokaryotic genomes and metagenomes. Information within P2CS can be accessed through a variety of routes-TCS complements can be browsed by metagenome, replicon or sequence cluster (and these genesets are available for download by users). Alternatively a variety of database-wide or taxon-specific searches are supported. Each TCS protein is fully annotated with sequence-feature information including replicon context, while properties of the predicted proteins can be queried against several external prediction servers to suggest homologues, interaction networks, sub-cellular localization and domain complements. Another unique feature of P2CS is the analysis of ORFeomes to identify TCS genes missed during genome annotation. Recent innovations for P2CS include a CGView representation of the distribution of TCS genes around a replicon, categorization of TCS genes based on gene organization, an expanded domain-based classification scheme, a P2CS 'gene cart' and categorization on the basis of sequence clusters.

A framework for classification of prokaryotic protein kinases

BACKGROUND

Overwhelming majority of the Serine/Threonine protein kinases identified by gleaning archaeal and eubacterial genomes could not be classified into any of the well known Hanks and Hunter subfamilies of protein kinases. This is owing to the development of Hanks and Hunter classification scheme based on eukaryotic protein kinases which are highly divergent from their prokaryotic homologues. A large dataset of prokaryotic Serine/Threonine protein kinases recognized from genomes of prokaryotes have been used to develop a classification framework for prokaryotic Ser/Thr protein kinases.

METHODOLOGY/PRINCIPAL FINDINGS

We have used traditional sequence alignment and phylogenetic approaches and clustered the prokaryotic kinases which represent 72 subfamilies with at least 4 members in each. Such a clustering enables classification of prokaryotic Ser/Thr kinases and it can be used as a framework to classify newly identified prokaryotic Ser/Thr kinases. After series of searches in a comprehensive sequence database we recognized that 38 subfamilies of prokaryotic protein kinases are associated to a specific taxonomic level. For example 4, 6 and 3 subfamilies have been identified that are currently specific to phylum proteobacteria, cyanobacteria and actinobacteria respectively. Similarly subfamilies which are specific to an order, sub-order, class, family and genus have also been identified. In addition to these, we also identify organism-diverse subfamilies. Members of these clusters are from organisms of different taxonomic levels, such as archaea, bacteria, eukaryotes and viruses.

CONCLUSION/SIGNIFICANCE

Interestingly, occurrence of several taxonomic level specific subfamilies of prokaryotic kinases contrasts with classification of eukaryotic protein kinases in which most of the popular subfamilies of eukaryotic protein kinases occur diversely in several eukaryotes. Many prokaryotic Ser/Thr kinases exhibit a wide variety of modular organization which indicates a degree of complexity and protein-protein interactions in the signaling pathways in these microbes.

Comparative genomic and phylogenetic analyses reveal the evolution of the core two-component signal transduction systems in enterobacteria

The two-component signal transduction system (TCST) consists of a histidine kinase (HK) and a response regulator (RR). TCSTs play important roles in sensing and reacting to environmental changes, and in bacterial pathogenesis. Previously, we have identified and characterized TCSTs in Erwinia amylovora, a severe plant enterobacterial pathogen, at genome-wide level. Here we conducted a comparative genomic analysis of TCSTs in 53 genomes of 16 enterobacterial species. These species include important plant, animal, human, and insect pathogenic, saprophytic or symbiotic microorganisms. Comparative genomic analysis revealed that enterobacteria contain eight pairs of core TCSTs. Phylogenetic trees reconstructed from a concatenation of the core set of TCSTs from enterobacteria and for individual TCST proteins from species in Proteobacteria showed that most TCST protein trees in the Enterobacteriaceae or in species of the γ-Proteobacteria agreed well with that of the corresponding 16S rRNA gene. It also showed that co-evolutionary relationships existed between cognate partners of the HKs and RRs. Several core TCSTs were quite ancient and universal based on phylogenomic analysis of protein structures. These results indicate that the core TCSTs are relatively conserved, and suggest that these enterobacteria may have maintained their ancient core TCSTs and might acquire specific new TCSTs for their survival in different environments or hosts, or may have evolved new functionalities of the core TCSTs for adaptation to different ecological niches.

P2CS: a two-component system resource for prokaryotic signal transduction research

BACKGROUND

With the escalation of high throughput prokaryotic genome sequencing, there is an ever-increasing need for databases that characterise, catalogue and present data relating to particular gene sets and genomes/metagenomes. Two-component system (TCS) signal transduction pathways are the dominant mechanisms by which micro-organisms sense and respond to external as well as internal environmental changes. These systems respond to a wide range of stimuli by triggering diverse physiological adjustments, including alterations in gene expression, enzymatic reactions, or protein-protein interactions.

DESCRIPTION

We present P2CS (Prokaryotic 2-Component Systems), an integrated and comprehensive database of TCS signal transduction proteins, which contains a compilation of the TCS genes within 755 completely sequenced prokaryotic genomes and 39 metagenomes. P2CS provides detailed annotation of each TCS gene including family classification, sequence features, functional domains, as well as genomic context visualization. To bypass the generic problem of gene underestimation during genome annotation, we also constituted and searched an ORFeome, which improves the recovery of TCS proteins compared to searches on the equivalent proteomes.

CONCLUSION

P2CS has been developed for computational analysis of the modular TCSs of prokaryotic genomes and metagenomes. It provides a complete overview of information on TCSs, including predicted candidate proteins and probable proteins, which need further curation/validation. The database can be browsed and queried with a user-friendly web interface at http://www.p2cs.org/.

Systems level analysis of two-component signal transduction systems in Erwinia amylovora: role in virulence, regulation of amylovoran biosynthesis and swarming motility

BACKGROUND

Two-component signal transduction systems (TCSTs), consisting of a histidine kinase (HK) and a response regulator (RR), represent a major paradigm for signal transduction in prokaryotes. TCSTs play critical roles in sensing and responding to environmental conditions, and in bacterial pathogenesis. Most TCSTs in Erwinia amylovora have either not been identified or have not yet been studied.

RESULTS

We used a systems approach to identify TCST and related signal transduction genes in the genome of E. amylovora. Comparative genomic analysis of TCSTs indicated that E. amylovora TCSTs were closely related to those of Erwinia tasmaniensis, a saprophytic enterobacterium isolated from apple flowers, and to other enterobacteria. Forty-six TCST genes in E. amylovora including 17 sensor kinases, three hybrid kinases, 20 DNA- or ligand-binding RRs, four RRs with enzymatic output domain (EAL-GGDEF proteins), and two kinases were characterized in this study. A systematic TCST gene-knockout experiment was conducted, generating a total of 59 single-, double-, and triple-mutants. Virulence assays revealed that five of these mutants were non-pathogenic on immature pear fruits. Results from phenotypic characterization and gene expression experiments indicated that several groups of TCST systems in E. amylovora control amylovoran biosynthesis, one of two major virulence factors in E. amylovora. Both negative and positive regulators of amylovoran biosynthesis were identified, indicating a complex network may control this important feature of pathogenesis. Positive (non-motile, EnvZ/OmpR), negative (hypermotile, GrrS/GrrA), and intermediate regulators for swarming motility in E. amylovora were also identified.

CONCLUSION

Our results demonstrated that TCSTs in E. amylovora played major roles in virulence on immature pear fruit and in regulating amylovoran biosynthesis and swarming motility. This suggested presence of regulatory networks governing expression of critical virulence genes in E. amylovora.

Relating gene expression data on two-component systems to functional annotations in Escherichia coli

Background

Obtaining physiological insights from microarray experiments requires computational techniques that relate gene expression data to functional information. Traditionally, this has been done in two consecutive steps. The first step identifies important genes through clustering or statistical techniques, while the second step assigns biological functions to the identified groups. Recently, techniques have been developed that identify such relationships in a single step.

Results

We have developed an algorithm that relates patterns of gene expression in a set of microarray experiments to functional groups in one step. Our only assumption is that patterns co-occur frequently. The effectiveness of the algorithm is demonstrated as part of a study of regulation by two-component systems in Escherichia coli. The significance of the relationships between expression data and functional annotations is evaluated based on density histograms that are constructed using product similarity among expression vectors. We present a biological analysis of three of the resulting functional groups of proteins, develop hypotheses for further biological studies, and test one of these hypotheses experimentally. A comparison with other algorithms and a different data set is presented.

Conclusion

Our new algorithm is able to find interesting and biologically meaningful relationships, not found by other algorithms, in previously analyzed data sets. Scaling of the algorithm to large data sets can be achieved based on a theoretical model.

Hierarchical binding of the TodT response regulator to its multiple recognition sites at the tod pathway operon promoter

The TodS and TodT proteins form a highly specific two-component regulatory system that controls the expression of genes involved in the degradation of toluene, benzene, and ethylbenzene via the toluene dioxygenase pathway. The catabolic genes of the toluene dioxygenase pathway are transcribed from a single promoter called P(todX) once the response regulator TodT is phosphorylated by the TodS sensor kinase in response to pathway substrates. We show here that TodT is a monomer in solution and that it binds to three specific sites in the P(todX) promoter, centered at -57, -85, and -106 with respect to the transcription start site. The -85 and -106 sites are pseudopalindromic, whereas the -57 site is half a palindrome. TodT binding to its target sites is sequential, as shown by electrophoresis mobility gel shift assays and footprinting. The binding affinity values of TodT, as determined by isothermal titration calorimetry, are 1.8+/-0.2, 5+/-0.4, and 6.3+/-0.8 microM for the -106, -85, and -57 sites, respectively, and the binding stoichiometry is one monomer per half-palindromic element. Mutational analysis revealed that all three sites contribute to P(todX) strength, although the most relevant site is the distal one with respect to the -10 extended element of the downstream promoter element. The C-TodT [C-terminal TodT fragment (amino acids 154-206)], a truncated variant of TodT that contains the C-terminal half of the protein bearing the DNA binding domain, binds in vitro to all three sites with affinity similar to that of the full-length protein. However, C-TodT, in contrast to the full-length regulator, does not activate in vitro transcription from P(todX). We discuss the consequences of the organization of the binding sites on transcriptional control and propose that the N-terminal domain of TodT is necessary for appropriate interactions with other transcriptional elements.

Structure and function predictions of the Msa protein in Staphylococcus aureus

BACKGROUND

Staphylococcus aureus is a human pathogen that causes a wide variety of life-threatening infections using a large number of virulence factors. One of the major global regulators used by S. aureus is the staphylococcal accessory regulator (sarA). We have identified and characterized a new gene (modulator of sarA: msa) that modulates the expression of sarA. Genetic and functional analysis shows that msa has a global effect on gene expression in S. aureus. However, the mechanism of Msa function is still unknown. Function predictions of Msa are complicated by the fact that it does not have a homologous partner in any other organism. This work aims at predicting the structure and function of the Msa protein.

RESULTS

Preliminary sequence analysis showed that Msa is a putative membrane protein. It would therefore be very difficult to purify and crystallize Msa in order to acquire structure information about this protein. We have used several computational tools to predict the physico-chemical properties, secondary structural features, topology, 3D tertiary structure, binding sites, motifs/patterns/domains and cellular location. We have built a consensus that is derived from analysis using different algorithms to predict several structural features. We confirm that Msa is a putative membrane protein with three transmembrane regions. We also predict that Msa has phosphorylation sites and binding sites suggesting functions in signal transduction.

CONCLUSION

Based on our predictions we hypothesise that Msa is a novel signal transducer that might be involved in the interaction of the S. aureus with its environment.

Identification of sensory and signal-transducing domains in two-component signaling systems

The availability of complete genome sequences of diverse bacteria and archaea makes comparative sequence analysis a powerful tool for analyzing signal transduction systems encoded in these genomes. However, most signal transduction proteins consist of two or more individual protein domains, which significantly complicates their functional annotation and makes automated annotation of these proteins in the course of large-scale genome sequencing projects particularly unreliable. This chapter describes certain common-sense protocols for sequence analysis of two-component histidine kinases and response regulators, as well as other components of the prokaryotic signal transduction machinery: Ser/Thr/Tyr protein kinases and protein phosphatases, adenylate and diguanylate cyclases, and c-di-GMP phosphodiesterases. These protocols rely on publicly available computational tools and databases and can be utilized by anyone with Internet access.

A census of membrane-bound and intracellular signal transduction proteins in bacteria: bacterial IQ, extroverts and introverts

BACKGROUND

Analysis of complete microbial genomes showed that intracellular parasites and other microorganisms that inhabit stable ecological niches encode relatively primitive signaling systems, whereas environmental microorganisms typically have sophisticated systems of environmental sensing and signal transduction.

RESULTS

This paper presents results of a comprehensive census of signal transduction proteins--histidine kinases, methyl-accepting chemotaxis receptors, Ser/Thr/Tyr protein kinases, adenylate and diguanylate cyclases and c-di-GMP phosphodiesterases--encoded in 167 bacterial and archaeal genomes, sequenced by the end of 2004. The data have been manually checked to avoid false-negative and false-positive hits that commonly arise during large-scale automated analyses and compared against other available resources. The census data show uneven distribution of most signaling proteins among bacterial and archaeal phyla. The total number of signal transduction proteins grows approximately as a square of genome size. While histidine kinases are found in representatives of all phyla and are distributed according to the power law, other signal transducers are abundant in certain phylogenetic groups but virtually absent in others.

CONCLUSION

The complexity of signaling systems differs even among closely related organisms. Still, it usually can be correlated with the phylogenetic position of the organism, its lifestyle, and typical environmental challenges it encounters. The number of encoded signal transducers (or their fraction in the total protein set) can be used as a measure of the organism's ability to adapt to diverse conditions, the 'bacterial IQ', while the ratio of transmembrane receptors to intracellular sensors can be used to define whether the organism is an 'extrovert', actively sensing the environmental parameters, or an 'introvert', more concerned about its internal homeostasis. Some of the microorganisms with the highest IQ, including the current leader Wolinella succinogenes, are found among the poorly studied beta-, delta- and epsilon-proteobacteria. Among all bacterial phyla, only cyanobacteria appear to be true introverts, probably due to their capacity to conduct oxygenic photosynthesis, using a complex system of intracellular membranes. The census data, available at http://www.ncbi.nlm.nih.gov/Complete_Genomes/SignalCensus.html, can be used to get an insight into metabolic and behavioral propensities of each given organism and improve prediction of the organism's properties based solely on its genome sequence.