Genomic analysis of protein kinases,protein phosphatases and two-component regulatory systems of the cyanobacterium Anabaena sp. strain PCC 7120

Anabaena sp. PCC 7120 is a cyanobacterium capable of performing several important biological functions: photosynthesis, nitrogen fixation, cell differentiation, cell-cell communication, etc. These activities require an extensive signaling capability in order to respond to the changing environment. Based on the genomic data, we have retrieved several gene families encoding signaling components. It is estimated that 211 genes encode two-component signaling elements, and 66 genes encode Ser/Thr kinases and phosphatases. These genes together represent 4.2% of the coding capacity of the whole genome, making Anabaena PCC 7120 a leading member among prokaryotes in terms of its signaling potential. It is known that two-component systems are composed of a few basic modules that can arrange into different structures best adapted for each signaling system. Many proteins in Anabaena PCC 7120 have incorporated both modules of two-component systems and catalytic domains of either Ser/Thr kinases or phosphatases. A family of 13 genes encode proteins with both a Ser/Thr kinase domain and a His kinase domain, and another four genes were also found whose products have both a response regulator domain and a Ser/Thr phosphatase domain. Of all the signaling proteins in Anabaena PCC 7120, about one third (35%) are conserved in the genome of the unicellular cyanobacterium strain Synechocystis sp. PCC 6803. Interestingly, one subfamily of His kinases and two subfamilies of response regulators are found in Anabaena PCC 7120 but are absent in Synechocystis PCC 6803. This study constitutes a basis for analyses of signal transduction in Anabaena PCC 7120 using functional genomic approaches.     

A genomic analysis of two-component signal transduction in Streptococcus pneumoniae

A genomics-based approach was used to identify the entire gene complement of putative two-component signal transduction systems (TCSTSs) in Streptococcus pneumoniae. A total of 14 open reading frames (ORFs) were identified as putative response regulators, 13 of which were adjacent to genes encoding probable histidine kinases. Both the histidine kinase and response regulator proteins were categorized into subfamilies on the basis of phylogeny. Through a systematic programme of mutagenesis, the importance of each novel TCSTS was determined with respect to viability and pathogenicity. One TCSTS was identified that was essential for the growth of S. pneumoniaeThis locus was highly homologous to the yycFG gene pair encoding the essential response regulator/histidine kinase proteins identified in Bacillus subtilis and Staphylococcus aureus. Separate deletions of eight other loci led in each case to a dramatic attenuation of growth in a mouse respiratory tract infection model, suggesting that these signal transduction systems are important for the in vivo adaptation and pathogenesis of S. pneumoniae. The identification of conserved TCSTSs important for both pathogenicity and viability in a Gram-positive pathogen highlights the potential of two-component signal transduction as a multicomponent target for antibacterial drug discovery.     

Mycobacterium tuberculosis Rv3651 is a triple sensor‐domain protein

The genome of the human pathogen Mycobacterium tuberculosis (Mtb) encodes ～4,400 proteins, but one third of them have unknown functions. We solved the crystal structure of Rv3651, a hypothetical protein with no discernible similarity to proteins with known function. Rv3651 has a three‐domain architecture that combines one cG MP‐specific phosphodiesterases, a denylyl cyclases and F hlA (GAF) domain and two P er‐A RNT‐S im (PAS) domains. GAF and PAS domains are sensor domains that are typically linked to signaling effector molecules. Unlike these sensor‐effector proteins, Rv3651 is an unusual sensor domain‐only protein with highly divergent sequence. The structure suggests that Rv3651 integrates multiple different signals and serves as a scaffold to facilitate signal transfer.     

Novel domains of the prokaryotic two-component signal transduction systems

The archetypal two-component signal transduction systems include a sensor histidine kinase and a response regulator, which consists of a receiver CheY-like domain and a DNA-binding domain. Sequence analysis of the sensor kinases and response regulators encoded in complete bacterial and archaeal genomes revealed complex domain architectures for many of them and allowed the identification of several novel conserved domains, such as PAS, GAF, HAMP, GGDEF, EAL, and HD-GYP. All of these domains are widely represented in bacteria, including 19 copies of the GGDEF domain and 17 copies of the EAL domain encoded in the Escherichia coli genome. In contrast, these novel signaling domains are much less abundant in bacterial parasites and in archaea, with none at all found in some archaeal species. This skewed phyletic distribution suggests that the newly discovered complexity of signal transduction systems emerged early in the evolution of bacteria, with subsequent massive loss in parasites and some horizontal dissemination among archaea. Only a few proteins containing these domains have been studied experimentally, and their exact biochemical functions remain obscure; they may include transformations of novel signal molecules, such as the recently identified cyclic diguanylate. Recent experimental data provide the first direct evidence of the participation of these domains in signal transduction pathways, including regulation of virulence genes and extracellular enzyme production in the human pathogens Bordetella pertussis and Borrelia burgdorferi and the plant pathogen Xanthomonas campestris. Gene-neighborhood analysis of these new domains suggests their participation in a variety of processes, from mercury and phage resistance to maintenance of virulence plasmids. It appears that the real picture of the complexity of phosphorelay signal transduction in prokaryotes is only beginning to unfold.     

Proteomic identification of M. tuberculosis protein kinase substrates: PknB recruits GarA, a FHA domain-containing protein, through activation loop-mediated interactions

Genes for functional Ser/Thr protein kinases (STPKs) are ubiquitous in prokaryotic genomes, but little is known about their physiological substrates and their actual involvement in bacterial signal transduction pathways. We report here the identification of GarA (Rv1827), a Forkhead-associated (FHA) domain-containing protein, as a putative physiological substrate of PknB, an essential Ser/Thr protein kinase from Mycobacterium tuberculosis. Using a global proteomic approach, GarA was found to be the best detectable substrate of the PknB catalytic domain in non-denatured whole-cell protein extracts from M. tuberculosis and the saprophyte Mycobacterium smegmatis. Enzymological and binding studies of the recombinant proteins demonstrate that docking interactions between the activation loop of PknB and the C-terminal FHA domain of GarA are required to enable efficient phosphorylation at a single N-terminal threonine residue, Thr22, of the substrate. The predicted amino acid sequence of the garA gene, including both the N-terminal phosphorylation motif and the FHA domain, is strongly conserved in mycobacteria and other related actinomycetes, suggesting a functional role of GarA in putative STPK-mediated signal transduction pathways. The ensuing model of PknB-GarA interactions suggests a substrate recruitment mechanism that might apply to other mycobacterial kinases bearing multiple phosphorylation sites in their activation loops.     

Molecular evolution of PAS domain-containing proteins of filamentous cyanobacteria through domain shuffling and domain duplication.

When the entire genome of a filamentous heterocyst-forming N2-fixing cyanobacterium, Anabaena sp. PCC 7120 (Anabaena) was determined in 2001, a large number of PAS domains were detected in signal-transducing proteins. The draft genome sequence is also available for the cyanobacterium, Nostoc punctiforme strain ATCC 29133 (Nostoc), that is closely related to Anabaena. In this study, we extracted all PAS domains from the Nostoc genome sequence and analyzed them together with those of Anabaena. Clustering analysis of all the PAS domains gave many specific pairings, indicative of evolutionary conservations. Ortholog analysis of PAS-containing proteins showed composite multidomain architecture in some cases of conserved domains and domains of disagreement between the two species. Further inspection of the domains of disagreement allowed us to trace them back in evolution. Thus, multidomain proteins could have been generated by duplication or shuffling in these cyanobacteria. The conserved PAS domains in the orthologous proteins were analyzed by structural fitting to the known PAS domains. We detected several subclasses with unique sequence features, which will be the target of experimental analysis.     

Evolution of prokaryotic two-component system signaling pathways: gene fusions and fissions

Two-component systems (TCSs) are common signal transduction systems, typically comprising paired histidine protein kinase (HK) and response regulator (RR) proteins. In many examples, it appears RR and HK genes have fused, producing a "hybrid kinase " We have characterized a set of prokaryotic genes encoding RRs, HKs, and hybrid kinases, enabling characterization of gene fusion and fission. Primary factors correlating with fusion rates are the presence of transmembrane helices in HKs and the presence of DNA-binding domains in RRs, features that require correct (and separate) spatial location. In the absence of such features, there is a relative abundance of fused genes. The order of paired HK and RR genes and the nucleotide distance between encoded domains also correlate with apparent gene fusion rates. We propose that localization requirements and relative positioning of encoded domains within TCS genes affect the function (and therefore retention) of hybrid kinases resulting from gene fusion.     

Survey, analysis and genetic organization of genes encoding eukaryotic-like signaling proteins on a cyanobacterial genome.

Bacteria usually use two-component systems for signal transduction, while eukaryotic organisms employ Ser/Thr and Tyr kinases and phosphatases for the same purpose. Many prokaryotes turn out to harbor Ser/Thr and Tyr kinases, Ser/Thr and Tyr phosphatases, and their accessory components as well. The sequence determination of the genome of the cyanobacterium Synechocystis sp. strain PCC 6803 offers the possibility to survey the extent of such molecules in a prokaryotic organism. This cyanobacterium possesses seven Ser/Thr kinases, seven Ser/Thr and Tyr phosphatases, one protein kinase interacting protein, one protein kinase regulatory subunit and several WD40-repeat-containing proteins. The majority of the protein phosphatases presented in this study were previously reported as hypothetical proteins. We analyze here the structure and genetic organization of these ORFs in the hope of providing a guidance for their functional analysis. Unlike their eukaryotic counterparts, many of these genes are clustered on the chromosome, and this genetic organization offers the opportunity to study their possible interaction. In several cases, genes of two-component transducers are found within the same cluster as those encoding a Ser/Thr kinase or a Ser/Thr phosphatase; the implication for signal transduction mechanism will be discussed.     

Structures and solution properties of two novel periplasmic sensor domains with c-type heme from chemotaxis proteins of Geobacter sulfurreducens: implications for signal transduction

Periplasmic sensor domains from two methyl-accepting chemotaxis proteins from Geobacter sulfurreducens (encoded by genes GSU0935 and GSU0582) were expressed in Escherichia coli. The sensor domains were isolated, purified, characterized in solution, and their crystal structures were determined. In the crystal, both sensor domains form swapped dimers and show a PAS-type fold. The swapped segment consists of two helices of about 45 residues at the N terminus with the hemes located between the two monomers. In the case of the GSU0582 sensor, the dimer contains a crystallographic 2-fold symmetry and the heme is coordinated by an axial His and a water molecule. In the case of the GSU0935 sensor, the crystals contain a non-crystallographic dimer, and surprisingly, the coordination of the heme in each monomer is different; monomer A heme has His-Met ligation and monomer B heme has His-water ligation as found in the GSU0582 sensor. The structures of these sensor domains are the first structures of PAS domains containing covalently bound heme. Optical absorption, electron paramagnetic resonance and NMR spectroscopy have revealed that the heme groups of both sensor domains are high-spin and low-spin in the oxidized and reduced forms, respectively, and that the spin-state interconversion involves a heme axial ligand replacement. Both sensor domains bind NO in their ferric and ferrous forms but bind CO only in the reduced form. The binding of both NO and CO occurs via an axial ligand exchange process, and is fully reversible. The reduction potentials of the sensor domains differ by 95 mV (− 156 mV and − 251 mV for sensors GSU0582 and GSU0935, respectively). The swapped dimerization of these sensor domains and redox-linked ligand switch might be related to the mechanism of signal transduction by these chemotaxis proteins.     

ABC-type neutral amino acid permease N-I is required for optimal diazotrophic growth and is repressed in the heterocysts of Anabaena sp. strain PCC 7120

Anabaena sp. strain PCC 7120 is a filamentous cyanobacterium that can fix N2 in differentiated cells called heterocysts. The products of Anabaena open reading frames (ORFs) all1046, all1047, all1284, alr1834 and all2912 were identified as putative elements of a neutral amino acid permease. Anabaena mutants of these ORFs were strongly affected (1-12% of the wild-type activity) in the transport of Pro, Phe, Leu and Gly and also impaired (17-30% of the wild-type activity) in the transport of Ala and Ser. These results identified those ORFs as the nat genes encoding the N-I neutral amino acid permease. According to amino acid sequence homologies, natA (all1046) and natE (all2912) encode ATPases, natC (all1047) and natD (all1284) encode transmembrane proteins, and natB (alr1834) encodes a periplasmic substrate-binding protein of an ABC-type uptake transporter. The natA, natC, natD and natE mutants showed defects in Gln and His uptake that were not observed in the natB mutant suggesting that NatB is not a binding protein for Gln or His. The nat mutants released hydrophobic amino acids to the medium, and amino acid release took place at higher levels in cultures incubated in the absence of combined N than in the presence of nitrate. Alanine was the amino acid released at highest levels, and its release was impaired in a mutant unable to develop heterocysts. The nat mutants were also impaired in diazotrophic growth, with natA, natC, natD and natE mutants showing more severe defects than the natB mutant. Expression of natA and natC, which constitute an operon, natCA, as well as of natB was studied and found to take place in vegetative cells but not in the heterocysts. These results indicate that the N-I permease is necessary for normal growth of Anabaena sp. strain PCC 7120 on N2, and that this permease has a role in the diazotrophic filament specifically in the vegetative cells.     

The hybrid histidine kinase Slr1759 of the cyanobacterium Synechocystis sp. PCC 6803 contains FAD at its PAS domain

The cyanobacterium Synechocystis sp. PCC 6803 harbours 47 histidine kinases (Hiks). Among these are hybrid histidine kinases with one or two response regulator domains as well as numerous Hiks with several sensory domains. One example is the hybrid histidine kinase Slr1759 (Hik14) that has two PAS domains arranged in tandem linked to a predicted GAF domain. Here, we show that a Slr1759 derivative recombinantly expressed in Escherichia coli has a flavin cofactor. Using truncated Slr1759 variants, it is shown that the flavin associates with the first PAS domain. The cofactor reconstitutes the activity of d-amino acid oxidase apoprotein from pig kidney, indicating that the flavin derivative is FAD. Furthermore, the Slr1759 histidine kinase domain indeed undergoes autophosphorylation in vitro. The phosphorylated product of a recombinant Slr1759 derivative is sensitive to acids, pointing to a histidine residue as the phosphate-accepting group.     

Programmed translational frameshifting is likely required for expressions of genes encoding putative nuclear protein kinases of the ciliate Euplotes octocarinatus

Three macronuclear genes encoding putative nuclear protein kinases of the ciliate Euplotes octocarinatus syngen 1 were isolated and sequenced. All three deduced gene products share significant properties with a group of recently identified nuclear serine/threonine protein kinases named Ndr. The three predicted proteins contain the twelve conserved catalytic subdomains of protein kinases and 22 near universally-conserved amino acids residues that are characteristic of serine/threonine protein kinases. In addition, there is an approximately 30 amino acid-peptide insertion between subdomains VII and VIII that contains a potential nuclear localization signal. Sequence analysis suggests that expression of the Eondr2 gene requires a + 1 programmed translational frameshift for its translation. Comparison of the deduced EoNdr2 with other known Ndr protein kinases implies that a + 1 ribosomal frameshift occurs at the motif AAATAA.     

Intrinsic and extrinsic approaches for detecting genes in a bacterial genome.

The unannotated regions of the Escherichia coli genome DNA sequence from the EcoSeq6 database, totaling 1,278 'intergenic' sequences of the combined length of 359,279 basepairs, were analyzed using computer-assisted methods with the aim of identifying putative unknown genes. The proposed strategy for finding new genes includes two key elements: i) prediction of expressed open reading frames (ORFs) using the GeneMark method based on Markov chain models for coding and non-coding regions of Escherichia coli DNA, and ii) search for protein sequence similarities using programs based on the BLAST algorithm and programs for motif identification. A total of 354 putative expressed ORFs were predicted by GeneMark. Using the BLASTX and TBLASTN programs, it was shown that 208 ORFs located in the unannotated regions of the E. coli chromosome are significantly similar to other protein sequences. Identification of 182 ORFs as probable genes was supported by GeneMark and BLAST, comprising 51.4% of the GeneMark 'hits' and 87.5% of the BLAST 'hits'. 73 putative new genes, comprising 20.6% of the GeneMark predictions, belong to ancient conserved protein families that include both eubacterial and eukaryotic members. This value is close to the overall proportion of highly conserved sequences among eubacterial proteins, indicating that the majority of the putative expressed ORFs that are predicted by GeneMark, but have no significant BLAST hits, nevertheless are likely to be real genes. The majority of the putative genes identified by BLAST search have been described since the release of the EcoSeq6 database, but about 70 genes have not been detected so far. Among these new identifications are genes encoding proteins with a variety of predicted functions including dehydrogenases, kinases, several other metabolic enzymes, ATPases, rRNA methyltransferases, membrane proteins, and different types of regulatory proteins.     

Screening for glycosylphosphatidylinositol (GPI)-dependent cell wall proteins in Saccharomyces cerevisiae

Open reading frames in the genome of Saccharomyces cerevisiae were screened for potential glycosylphosphatidylinositol (GPI)-attached proteins. The identification of putative GPI-attached proteins was based on three criteria: the presence of a GPI-attachment signal sequence, a signal sequence for secretion and a serine- or threonine-rich sequence. In all, 53 ORFs met these three criteria and 38 were further analyzed as follows. The sequence encoding the 40 C-terminal amino acids of each was fused with the structural gene for a reporter protein consisting of a secretion signal, α-galactosidase and a hemagglutinin (HA) epitope, and examined for the ability to become incorporated into the cell wall. On this basis, 14 of fusion proteins were classified as GPI-dependent cell wall proteins because cells expressing these fusion proteins: (i) had high levels of α-galactosidase activity on their surface; (ii) released significant amounts of the fusion proteins from the membrane on treatment with phosphatidylinositol-specific phospholipase C (PI-PLC); and (iii) released fusion proteins from the cell wall following treatment with laminarinase. Of the 14 identified putative GPI-dependent cell wall proteins, 12 had novel ORFs adjacent to their GPI-attachment signal sequence. Amino acid sequence alignment of the C-terminal sequences of the 12 ORFs, together with those of known cell wall proteins, reveals some sequence similarities among them. Received: 1 September 1997 / Accepted: 20 November 1997