A method to locate protein coding sequences in DNA of prokaryotic systems

PAAR proteins act as the 'sorting hat' of the type VI secretion system

Bacteria exist in polymicrobial environments and compete to prevail in a niche. The type VI secretion system (T6SS) is a nanomachine employed by Gram-negative bacteria to deliver effector proteins into target cells. Consequently, T6SS-positive bacteria produce a wealth of antibacterial effector proteins to promote their survival among a prokaryotic community. These toxins are loaded onto the VgrG–PAAR spike and Hcp tube of the T6SS apparatus and recent work has started to document the specificity of effectors for certain spike components. Pseudomonas aeruginosa encodes several PAAR proteins, whose roles have been poorly investigated. Here we describe a phospholipase family antibacterial effector immunity pair from Pseudomonas aeruginosa and demonstrate that a specific PAAR protein is necessary for the delivery of the effector and its cognate VgrG. Furthermore, the PAAR protein appears to restrict the delivery of other phospholipase effectors that utilise distinct VgrG proteins. We provide further evidence for competition for PAAR protein recruitment to the T6SS apparatus, which determines the identities of the delivered effectors.

A Type VI Secretion System Trans-Kingdom Effector Is Required for the Delivery of a Novel Antibacterial Toxin in Pseudomonas aeruginosa

Pseudomonas aeruginosa has evolved multiple strategies to disarm and take advantage of its host. For this purpose, this opportunist pathogen has particularly developed protein secretion in the surrounding medium or injection into host cells. Among this, the type VI secretion system (T6SS) is utilized to deliver effectors into eukaryotic host as well as target bacteria. It assembles into a contractile bacteriophage tail-like structure that functions like a crossbow, injecting an arrow loaded with effectors into the target cell. The repertoire of T6SS antibacterial effectors of P. aeruginosa is remarkably broad to promote environmental adaptation and survival in various bacterial communities, and presumably in the eukaryotic host too. Here, we report the discovery of a novel pair of antibacterial effector and immunity of P. aeruginosa, Tle3 and Tli3. Tli3 neutralizes the toxicity of Tle3 in the periplasm to protect from fratricide intoxication. The characterization of the secretion mechanism of Tle3 indicates that it requires a cytoplasmic adaptor, Tla3, to be targeted and loaded onto the VgrG2b spike and thus delivered by the H2-T6SS machinery. Tla3 is different from the other adaptors discovered so far and defines a novel family among T6SS with a DUF2875. Interestingly, this led us to discover that VgrG2b that we previously characterized as an anti-eukaryotic effector possesses an antibacterial activity as well, as it is toxic towards Escherichia coli. Excitingly Tli3 can counteract VgrG2b toxicity. VgrG2b is thus a novel trans-kingdom effector targeting both bacteria and eukaryotes. VgrG2b represents an interesting target for fighting against P. aeruginosa in the environment and in the context of host infection.

Genome wide identification and experimental validation of Pseudomonas aeruginosa Tat substrates

In bacteria, the twin-arginine translocation (Tat) pathway allows the export of folded proteins through the inner membrane. Proteins targeted to this system are synthesized with N-terminal signal peptides bearing a conserved twin-arginine motif. The Tat pathway is critical for many bacterial processes including pathogenesis and virulence. However, the full set of Tat substrates is unknown in many bacteria, and the reliability of in silico prediction methods largely uncertain. In this work, we performed a combination of in silico analysis and experimental validation to identify a core set of Tat substrates in the opportunistic pathogen Pseudomonas aeruginosa. In silico analysis predicted 44 putative Tat signal peptides in the P. aeruginosa PA14 proteome. We developed an improved amidase-based Tat reporter assay to show that 33 of these are real Tat signal peptides. In addition, in silico analysis of the full translated genome revealed a Tat candidate with a missassigned start codon. We showed that it is a new periplasmic protein in P. aeruginosa. Altogether we discovered and validated 34 Tat substrates. These show little overlap with Escherichia coli Tat substrates, and functional analysis points to a general role for the P. aeruginosa Tat system in the colonization of environmental niches and pathogenicity.

Contribution of the Twin Arginine Translocation system to the exoproteome of Pseudomonas aeruginosa

The opportunistic pathogen Pseudomonas aeruginosa uses secretion systems to deliver exoproteins into the environment. These exoproteins contribute to bacterial survival, adaptation, and virulence. The Twin arginine translocation (Tat) export system enables the export of folded proteins into the periplasm, some of which can then be further secreted outside the cell. However, the full range of proteins that are conveyed by Tat is unknown, despite the importance of Tat for the adaptability and full virulence of P. aeruginosa. In this work, we explored the P. aeruginosa Tat-dependent exoproteome under phosphate starvation by two-dimensional gel analysis. We identified the major secreted proteins and new Tat-dependent exoproteins. These exoproteins were further analyzed by a combination of in silico analysis, regulation studies, and protein localization. Altogether we reveal that the absence of the Tat system significantly affects the composition of the exoproteome by impairing protein export and affecting gene expression. Notably we discovered three new Tat exoproteins and one novel type II secretion substrate. Our data also allowed the identification of two new start codons highlighting the importance of protein annotation for subcellular predictions. The new exoproteins that we identify may play a significant role in P. aeruginosa pathogenesis, host interaction and niche adaptation.

BtpB, a novel Brucella TIR-containing effector protein with immune modulatory functions

Several bacterial pathogens have TIR domain-containing proteins that contribute to their pathogenesis. We identified a second TIR-containing protein in Brucella spp. that we have designated BtpB. We show it is a potent inhibitor of TLR signaling, probably via MyD88. BtpB is a novel Brucella effector that is translocated into host cells and interferes with activation of dendritic cells. In vivo mouse studies revealed that BtpB is contributing to virulence and control of local inflammatory responses with relevance in the establishment of chronic brucellosis. Together, our results show that BtpB is a novel Brucella effector that plays a major role in the modulation of host innate immune response during infection.

Role of thymine in protein coding frames of mRNA sequences

Distribution of thymine in protein coding mRNA sequences has been studied here. Our study suggest that thymine in protein coding sequences are not randomly distributed but with probability. Frame1 prefers to have definite amount of thymine. It is observed that the thymine content of frame 4 is also involved in protein coding. Frame 3 prefers to have least amount of thymine. However, frame 2 and frame 6 shows a variable degree of thymine content. The mRNA sequences of heterosexual animals, particularly, the human show a different distribution profile (less thymine in frame 1) compared to that of yeast and plants.

A novel secondary acyl chain in the lipopolysaccharide of Bordetella pertussis required for efficient infection of human macrophages

Lipopolysaccharide is one of the major constituents of the Gram-negative bacterial outer membrane and is a potent stimulator of the host innate immune response. The biosynthesis of the lipid A moiety of lipopolysaccharide is a complex process in which multiple gene products are involved. Two late lipid A acyl transferases, LpxL and LpxM, were first identified in Escherichia coli and shown to be responsible for the addition of secondary acyl chains to the 2' and 3' positions of lipid A, respectively. Here, we describe the identification of two lpxL homologues in the genome of Bordetella pertussis. We show that one of them, LpxL2, is responsible for the addition of the secondary myristate group that is normally present at the 2' position of B. pertussis lipid A, whereas the other one, LpxL1, mediates the addition of a previously unrecognized secondary 2-hydroxy laurate at the 2 position. Increased expression of lpxL1 results in the appearance of a hexa-acylated lipopolysaccharide form with strongly increased endotoxic activity. In addition, we show that an lpxL1-deficient mutant of B. pertussis displays a defect in the infection of human macrophages.

Function of neisserial outer membrane phospholipase a in autolysis and assessment of its vaccine potential

Outer membrane phospholipase A (OMPLA) is an outer membrane-localized enzyme, present in many gram-negative bacterial species. It is implicated in the virulence of several pathogens. Here, we investigated the presence, function, and vaccine potential of OMPLA in the human pathogen Neisseria meningitidis. Immunoblot analysis showed the presence of OMPLA in 28 out of 33 meningococcal strains investigated. The OMPLA-negative strains all contained a pldA gene, but these alleles contained premature stop codons. All six Neisseria gonorrhoeae strains tested, but only two out of seven commensal neisserial strains investigated, expressed OMPLA, showing that OMPLA is expressed by, but not limited to, many pathogenic neisserial strains. The function of OMPLA was investigated by assessing the phenotypes of isogenic strains, expressing no OMPLA, expressing wild-type levels of OMPLA, or overexpressing OMPLA. OMPLA exhibited phospholipase activity against endogenous phospholipids. Furthermore, OMPLA was characterized as an autolysin that acted under specific conditions, such as prolonged growth of the bacteria. The vaccine potential of the protein was investigated by immunizing mice with in vitro refolded, recombinant OMPLA. High levels of antibody titers were obtained, but the murine sera were neither bactericidal nor protective. Also, convalescent patients and vaccinee sera did not contain detectable levels of anti-OMPLA antibodies, indicating that OMPLA may not be sufficiently immunogenic to be included in a meningococcal vaccine.

Complete sequence of the gene encoding a chorionic gonadotropin-like protein from Xanthomonas maltophilia

Our laboratory has previously reported that: (i) Xanthomonas maltophilia (Xm) produces a protein which has immunological resemblance to the beta-subunit of human chorionic gonadotropin (hCG) and (ii) possesses a high-affinity receptor which binds holo hCG, and the endogenous ligand, Xm chorionic gonadotropin (xCG), but does not bind human luteinizing hormone (hLH). We have also previously published a 492-bp partial nucleotide sequence of the gene (xcg) coding for xCG. We report herein the entire xcg sequence of 1362 bp, which codes for a 48-kDa protein. This sequence confirmed the 492-bp sequence, as well as two partial amino acid (aa) sequences which we have previously reported. The sequence has a region which is homologous to aa 56-139 of the beta-subunit of hCG, and a second region homologous to the C-terminal tail of hCG. This is the first report of a prokaryotic gene homologous to the hCG beta-subunit-encoding gene.

A statistical analytical approach to decipher information from biological sequences: application to murine splice-site analysis and prediction

A simple statistical approach for the analysis of biological sequences, such as splice-sites, promoter regions, helices and extended structure forming regions or any other sequence dependent functional entities in proteins, is presented. The approach has been proved useful to develop a method for prediction of such entities in newly available sequences. We first search for invariant sequence features of each functional entity from the experimentally available sequences and identify a set of 'like' sequences with similar sequence features. In the next step, concrete features of sequence entities in terms of occurrences of smaller subsequences are identified at various positions which are used as a knowledge base to select potential functional entities from the identified 'like' sequences. The third step consists of refinement of this pattern learning, statistical improvements of the knowledge base weight matrices, and finally its application to predict functional entities in newly available sequences. Such an analysis is operationally described for murine splice-site predictions. Regions comprising -30 to +30 nucleotides from the splice-junction at the murine splice-sites (donors and acceptors), reported earlier, were analyzed. Invariant sequence-specific features in terms of monomer frequency average were used to identify splice-site-like sequences in the EMBL murine DNA sequence data base. The frequencies of occurrence of mono-, di-, tri- and tetranucleotides in the known splice-sites were studied in comparison with the splice-site-like sequences; the significant differences in their occurrences were extracted as statistical knowledge coded in weight matrices for computer to identify potential splice-sites. The algorithm was refined and a method was developed to predict potential splice-sites in a given murine DNA; the analysis was also extended to human DNA. The success rate of the method to predict correct splice-sites in these species is found to be 80% and 85%, respectively. The major strength of this method lies in reducing significantly the number of false positives which are normally picked up in such analysis.

Genetic rearrangement associated with in vivo mucoid conversion of Pseudomonas aeruginosa PAO is due to insertion elements

The conversion of Pseudomonas aeruginosa PAO to the mucoid phenotype has been reported for a chronic pulmonary infection model in rats (D. E. Woods, P. A. Sokol, L. E. Bryan, D. G. Storey, S. J. Mattingly, H. J. Vogel, and H. Ceri, J. Infect. Dis. 163:143-149, 1991). This conversion was associated with a genetic rearrangement upstream of the exotoxin A gene. To characterize the genetic rearrangement, the region upstream of the toxA gene was cloned from PAO, PAO-muc (a mucoid strain), and PAO-rev (a nonmucoid revertant strain). The nucleotide sequence of a 4.8-kb fragment from PAO-muc was determined. A+T-rich regions of approximately 2 kb (IS-PA-4) and 0.4 kb (IS-PA-5) were identified in this fragment. DNA probes constructed internal to these regions hybridized to PAO-muc but not to PAO or PAO-rev, suggesting that PAO-muc contains an insertion element. Sequence analysis of the nonmucoid clones indicated that a 2,561-bp fragment corresponding to IS-PA-4 and a 992-bp fragment corresponding to IS-PA-5 were not present in PAO or PAO-rev. Both nonmucoid clones, however, contained in the same location as IS-PA-4, a 1,313-bp region which was not present in PAO-muc. DNA probes complementary to this sequence, designated IS-PA-6, did not hybridize with PAO-muc, indicating that this sequence had been replaced upon conversion to the mucoid phenotype. Between IS-PA-4 and IS-PA-5 there was a 500-bp sequence which was 94% identical to the 500-bp sequence downstream of IS-PA-6. These insertion elements had some DNA sequence similarity to plasmid and transposon sequences, suggesting that they may be of plasmid origin. IS-PA-4 and IS-PA-5 were shown also to be present in two mucoid isolates from cystic fibrosis patients. The insertions occurred in the same location upstream of the toxA gene, suggesting that this type of genetic recombination may also be associated with mucoid conversion in some P. aeruginosa clinical isolates.

Combined use of sequence similarity and codon bias for coding region identification

A computer program called BLASTX was previously shown to be effective in identifying and assigning putative function to likely protein coding regions by detecting significant similarity between a conceptually translated nucleotide query sequence and members of a protein sequence database. We present and assess the sensitivity of a new option to this software tool, herein called BLASTC, which employs information obtained from biases in codon utilization, along with the information obtained from sequence similarity. A rationale for combining these diverse information sources was derived, and analyses of the information available from codon utilization in several species were performed, with wide variation seen. Codon bias information was found on average to improve the sensitivity of detection of short coding regions of human origin by about a factor of 5. The implications of combining information sources on the interpretation of positive findings are discussed.

Transfer functions of the conjugative integrating element pSAM2 from Streptomyces ambofaciens: characterization of a kil-kor system associated with transfer

pSAM2 is an 11-kb integrating element from Streptomyces ambofaciens. During matings, pSAM2 can be transferred at high frequency, forming pocks, which are zones of growth inhibition of the recipient strain. The nucleotide sequences of the regions involved in pSAM2 transfer, pock formation, and maintenance have been determined. Seven putative open reading frames with the codon usage typical of Streptomyces genes have been identified: traSA (306 amino acids [aa]), orf84 (84 aa), spdA (224 aa), spdB (58 aa), spdC (51 aa), spdD (104 aa), and korSA (259 aa). traSA is essential for pSAM2 intermycelial transfer and pock formation. It could encode a protein with similarities to the major transfer protein, Tra, of pIJ101. TraSA protein contains a possible nucleotide-binding sequence and a transmembrane segment. spdA, spdB, spdC, and spdD influence pock size and transfer efficiency and may be required for intramycelial transfer. A kil-kor system similar to that of pIJ101 is associated with pSAM2 transfer: the korSA (kil-override) gene product could control the expression of the traSA gene, which has lethal effects when unregulated (Kil phenotype). The KorSA protein resembles KorA of pIJ101 and repressor proteins belonging to the GntR family. Thus, the integrating element pSAM2 possesses for transfer general features of nonintegrating Streptomyces plasmids: different genes are involved in the different steps of the intermycelial and intramycelial transfer, and a kil-kor system is associated with transfer. However, some differences in the functional properties, organization, and sizes of the transfer genes compared with those of other Streptomyces plasmids have been found.

Identification and characterization of the pupB gene encoding an inducible ferric-pseudobactin receptor of Pseudomonas putida WCS358

Pseudomonas putida WCS358 can transport iron complexed to a wide variety of pseudobactins produced by other Pseudomonas strains. The pupB gene encoding an outer membrane ferric-pseudobactin receptor was isolated from a genomic library of P. putida WCS358. The PupB receptor facilitated iron transport via two distinct heterologous siderophores, i.e. pseudobactin BN8 and pseudobactin BN7. The amino acid sequence deduced from the nucleotide sequence consisted of 804 amino acids (molecular weight 88,369) of which the N-terminal part was very similar to a prokaryotic leader peptide. The mature protein shared significant homology with the receptor for ferric-pseudobactin 358 (PupA) and contained three regions common to TonB-dependent receptor proteins of Escherichia coli. Interestingly, PupB expression was only observed in cells cultured in iron-deficient medium containing pseudobactin BN8 or pseudobactin BN7. This expression required a transcriptional unit, pupR, identified upstream of the structural pupB gene. Transposon Tn5 insertion mutants defective in PupB production still exhibited uptake of iron via pseudobactin BN8, although with reduced efficiency. Apparently, an additional transport system for this ferric-siderophore complex operates in this strain. In addition to pseudobactin BN8 also other heterologous siderophores were capable of inducing synthesis of specific high-molecular-weight outer membrane proteins in strain WCS358, which suggests the existence of multiple siderophore-inducible iron transport systems in this strain.

Identification and characterization of the exbB, exbD and tonB genes of Pseudomonas putida WCS358: their involvement in ferric-pseudobactin transport

Catechol-cephalosporins are siderophore-like antibiotics which are taken up by cells of Pseudomonas putida WCS358 via the ferric-siderophore transport pathway. Mutants of strain WCS358 were isolated that are resistant to high concentrations of these antibiotics. These mutants failed to grow under iron-limiting conditions, and could not utilize different ferric-siderophores. The mutants fall in three complementation groups. The nucleotide sequence determination identified three contiguous open reading frames, which were homologous to the exbB, exbD and tonB genes of Escherichia coli respectively. The deduced amino acid sequence of P. putida ExbB showed 58.6% homology with its E. coli homologue, but, unlike the E. coli protein, it has a N-terminal extension of 91 amino acids. The ExbD proteins are 64.8% homologous, whereas the TonB proteins only show 27.7% homology. The P. putida exbB gene could complement an E. coli exbB mutation, but the TonB proteins were not interchangeable between the species. It is concluded that P. putida WCS358 contains an energy-coupling system between the membranes, for active transport across the outer membrane, which is comprised of a TonB-like energy-transducing protein and two accessory proteins. This system is similar to, but not completely compatible with, the E. coli system.

Assessment of protein coding measures

A number of methods for recognizing protein coding genes in DNA sequence have been published over the last 13 years, and new, more comprehensive algorithms, drawing on the repertoire of existing techniques, continue to be developed. To optimize continued development, it is valuable to systematically review and evaluate published techniques. At the core of most gene recognition algorithms is one or more coding measures--functions which produce, given any sample window of sequence, a number or vector intended to measure the degree to which a sample sequence resembles a window of 'typical' exonic DNA. In this paper we review and synthesize the underlying coding measures from published algorithms. A standardized benchmark is described, and each of the measures is evaluated according to this benchmark. Our main conclusion is that a very simple and obvious measure--counting oligomers--is more effective than any of the more sophisticated measures. Different measures contain different information. However there is a great deal of redundancy in the current suite of measures. We show that in future development of gene recognition algorithms, attention can probably be limited to six of the twenty or so measures proposed to date.

Cloning and sequencing of IS1086, an Alcaligenes eutrophus insertion element related to IS30 and IS4351

A new insertion sequence (IS), designated IS1086, was isolated from Alcaligenes eutrophus CH34 by being trapped in plasmid pJV240, which contains the Bacillus subtilis sacB and sacR genes. The 1,106-bp IS1086 element contains partially matched (22 of 28 bp) terminal-inverted repeats and a long open reading frame. Hybridization data suggest the presence of one copy of IS1086 in the strain CH34 heavy-metal resistance plasmid pMOL28 and at least two copies in its chromosome. Analysis of the IS1086 nucleotide sequence revealed striking homology with two other IS elements, IS30 and IS4351, suggesting that they are three close members in a family of phylogenetically related insertion sequences. One open reading frame of the Spiroplasma citri phage SpV1-R8A2 B was also found to be related to this IS family but to a lesser extent. Comparison of the G+C contents of IS30 and IS1086 revealed that they conform to their respective hosts (46 versus 50% for IS30 and Escherichia coli and 64.5% for IS1086 and A. eutrophus). The pressure on the AT/GC ratio led to a very different codon usage in these two closely related IS elements. Results suggesting that IS1086 transposition might be activated by some forms of stress are discussed.

Rhodobacter sphaeroides rdxA, a homolog of Rhizobium meliloti fixG, encodes a membrane protein which may bind cytoplasmic [4Fe-4S] clusters

In the photosynthetic bacterium Rhodobacter sphaeroides, a chromosomal gene, rdxA, which encodes a 52-kDa protein, was found to be homologous to fixG, the first gene of a Rhizobium meliloti nitrogen fixation operon on the pSym plasmid (D. Kahn, M. David, O. Domergue, M.-L. Daveran, J. Ghai, P. R. Hirsch, and J. Batut, J. Bacteriol. 171:929-939, 1989). The deduced amino acid sequences of RdxA and FixG are 53% identical and 73% similar; sequence analyses suggested that each has five transmembrane helices and a central region resembling bacterial-type ferredoxins. Translational fusion proteins with an alkaline phosphatase reporter group were expressed in both R. sphaeroides and Escherichia coli and were used to assess the membrane topology of RdxA. Its ferredoxinlike sequence, which may bind two [4Fe-4S] centers, was found to be cytoplasmically located. Genetic disruptions showed that rdxA is not essential for nitrogen fixation in R. sphaeroides. Immediately downstream of rdxA, an open reading frame (ORFT2) that encoded a 48-kDa protein was found. This DNA sequence was not homologous to any region of the R. meliloti fixG operon. The N-terminal sequence of the ORFT2 gene product resembled amino acid sequences found in members of the GntR family of regulatory proteins (D. J. Haydon and J. R. Guest, FEMS Microbiol. Lett. 79:291-296, 1991). The rdxA gene was localized to the smaller of two R. sphaeroides chromosomes, upstream of and divergently transcribed from hemT, which encodes one of two 5-aminolevulinate synthase isozymes. The rdxA and hemT genes may share a transcriptional regulatory region. Southern hybridization analysis demonstrated the presence of an rdxA homolog on the R. sphaeroides large chromosome. The functions of this homolog, like those of rdxA, remain to be determined, but roles in oxidation-reduction processes are likely.

Determination of eukaryotic protein coding regions using neural networks and information theory

Our previous work applied neural network techniques to the problem of discriminating open reading frame (ORF) sequences taken from introns versus exons. The method counted the codon frequencies in an ORF of a specified length, and then used this codon frequency representation of DNA fragments to train a neural net (essentially a Perceptron with a sigmoidal, or "soft step function", output) to perform this discrimination. After training, the network was then applied to a disjoint "predict" set of data to assess accuracy. The resulting accuracy in our previous work was 98.4%, exceeding accuracies reported in the literature at that time for other algorithms. Here, we report even higher accuracies stemming from calculations of mutual information (a correlation measure) of spatially separated codons in exons, and in introns. Significant mutual information exists in exons, but not in introns, between adjacent codons. This suggests that dicodon frequencies of adjacent codons are important for intron/exon discrimination. We report that accuracies obtained using a neural net trained on the frequency of dicodons is significantly higher at smaller fragment lengths than even our original results using codon frequencies, which were already higher than simple statistical methods that also used codon frequencies. We also report accuracies obtained from including codon and dicodon statistics in all six reading frames, i.e. the three frames on the original and complement strand. Inclusion of six-frame statistics increases the accuracy still further. We also compare these neural net results to a Bayesian statistical prediction method that assumes independent codon frequencies in each position. The performance of the Bayesian scheme is poorer than any of the neural based schemes, however many methods reported in the literature either explicitly, or implicitly, use this method. Specifically, Bayesian prediction schemes based on codon frequencies achieve 90.9% accuracy on 90 codon ORFs, while our best neural net scheme reaches 99.4% accuracy on 60 codon ORFs. "Accuracy" is defined as the average of the exon and intron sensitivities. Achievement of sufficiently high accuracies on short fragment lengths can be useful in providing a computational means of finding coding regions in unannotated DNA sequences such as those arising from the mega-base sequencing efforts of the Human Genome Project. We caution that the high accuracies reported here do not represent a complete solution to the problem of identifying exons in "raw" base sequences. The accuracies are considerably lower from exons of small length, although still higher than accuracies reported in the literature for other methods. Short exon lengths are not uncommon.(ABSTRACT TRUNCATED AT 400 WORDS)

Construction and expression of diphtheria toxin-encoding gene derivatives in Escherichia coli

We reported earlier that in the periplasmic space of Escherichia coli, truncated derivatives of diphtheria toxin undergo limited proteolysis [Zdanovsky et al., Mol. Biol. 22 (1988) 1037-1293]. Here, we present data indicating that this proteolysis is reduced in cells bearing a mutation in the degP gene. We have also constructed hybrid genes whose products are not secreted into the periplasm. These hybrid genes were expressed in E. coli from both the pR promoter, controlled by the heat-inducible CI857 repressor, and from the P(lac) promoter, controlled by the IPTG-inducible LacI repressor. The latter system proved to be more productive.

Identification of replication and stability functions in the complete nucleotide sequence of plasmid pUH24 from the cyanobacterium Synechococcus sp. PCC 7942

The complete nucleotide sequence is presented for pUH24, the small plasmid of Synechococcus sp. PCC 7942. pUH24 consists of 7835bp and has a G + C content of 59%. The distribution of translation start and stop codons in the sequence allows 36 open reading frames that potentially encode polypeptides of 50 or more amino acids. We postulate that eight of these open reading frames are actual coding sequences. A region has been identified, by experiment, that contains two functions, designated pmaA and pmaB, involved in the segregational stability of the plasmid. The minimal region of pUH24 fully capable of supporting autonomous replication consists of a 3.6kb DNA fragment, which is almost entirely occupied by two overlapping genes most likely coding for essential replication proteins (repA and repB).

Nucleotide sequences of the Acinetobacter calcoaceticus benABC genes for benzoate 1,2-dioxygenase reveal evolutionary relationships among multicomponent oxygenases

The nucleotide sequences of the Acinetobacter calcoaceticus benABC genes encoding a multicomponent oxygenase for the conversion of benzoate to a nonaromatic cis-diol were determined. The enzyme, benzoate 1,2-dioxygenase, is composed of a hydroxylase component, encoded by benAB, and an electron transfer component, encoded by benC. Comparison of the deduced amino acid sequences of BenABC with related sequences, including those for the multicomponent toluate, toluene, benzene, and naphthalene 1,2-dioxygenases, indicated that the similarly sized subunits of the hydroxylase components were derived from a common ancestor. Conserved cysteine and histidine residues may bind a [2Fe-2S] Rieske-type cluster to the alpha-subunits of all the hydroxylases. Conserved histidines and tyrosines may coordinate a mononuclear Fe(II) ion. The less conserved beta-subunits of the hydroxylases may be responsible for determining substrate specificity. Each dioxygenase had either one or two electron transfer proteins. The electron transfer component of benzoate dioxygenase, encoded by benC, and the corresponding protein of the toluate 1,2-dioxygenase, encoded by xylZ, were each found to have an N-terminal region which resembled chloroplast-type ferredoxins and a C-terminal region which resembled several oxidoreductases. These BenC and XylZ proteins had regions similar to certain monooxygenase components but did not appear to be evolutionarily related to the two-protein electron transfer systems of the benzene, toluene, and naphthalene 1,2-dioxygenases. Regions of possible NAD and flavin adenine dinucleotide binding were identified.

Characterization of two genes, glpQ and ugpQ, encoding glycerophosphoryl diester phosphodiesterases of Escherichia coli

The nucleotide sequences of the glpQ and ugpQ genes of Escherichia coli, which both encode glycerophosphoryl diester phosphodiesterases, were determined. The glpQ gene encodes a periplasmic enzyme of 333 amino acids, produced initially with a 25 residue long signal sequence, while ugpQ codes for a cytoplasmic protein of 247 amino acids. Despite differences in size and cellular location, significant similarity in the primary structures of the two enzymes was found suggesting a common evolutionary origin. The 3' end of the ugpQ gene overlaps an open reading frame that is transcribed in the opposite direction. This open reading frame encodes a polypeptide with an unusual composition, i.e., 46 of the 146 amino acids are Gln or Asn. This polypeptide and the UgpQ protein were identified in an in vitro transcription/translation system as proteins with apparent molecular weights of 19.5 and 27 kDa, respectively.

IS801, an insertion sequence element isolated from Pseudomonas syringae pathovar phaseolicola

A transposable element, designated IS801, was isolated from strain LR781 of Pseudomonas syringae pathovar phaseolicola in two independent events using the entrapment plasmid, pUCD800. IS801 is 1517 base pairs in length and contains open reading frames that potentially encode proteins of 311 and 172 amino acids, as well as smaller proteins. Unlike most other prokaryotic transposable elements, IS801 lacks terminal repeats. Sequence analysis revealed two target pentamers for IS801 insertion that differ by one base pair. One copy of IS801 generated a perfect duplication of its target, TGAAC. The second copy of IS801 was flanked by the target, TGGAC, at one end, and TGAAC at the other end. A third copy of IS801 was cloned from pMMC7105, an indigenous plasmid of strain LR781, and it was flanked by copies of the pentamer TGAAC.

The ferric-pseudobactin receptor PupA of Pseudomonas putida WCS358: homology to TonB-dependent Escherichia coli receptors and specificity of the protein

The initial step in the uptake of iron via ferric pseudobactin by the plant-growth-promoting Pseudomonas putida strain WCS358 is binding to a specific outer-membrane protein. The nucleotide sequence of the pupA structural gene, which codes for a ferric pseudobactin receptor, was determined. It contains a single open reading frame which potentially encodes a polypeptide of 819 amino acids, including a putative N-terminal signal sequence of 47 amino acids. Significant homology, concentrated in four boxes, was found with the TonB-dependent receptor proteins of Escherichia coli. The pupA mutant MH100 showed a residual efficiency of 30% in the uptake of 55Fe3+ complexed to pseudobactin 358, whereas the iron uptake of four other pseudobactins was not reduced at all. Cells of strain WCS374 supplemented with the pupA gene of strain WCS358 could transport ferric pseudobactin 358 but showed no affinity for three other pseudobactins. It is concluded that PupA is a specific receptor for ferric pseudobactin 358, and that strain WCS358 produces at least one other receptor for other pseudobactins.

Nitrous oxide reductase from denitrifying Pseudomonas stutzeri. Genes for copper-processing and properties of the deduced products, including a new member of the family of ATP/GTP-binding proteins

Nitrous oxide (N2O) respiration by the denitrifying bacterium Pseudomonas stutzeri requires the synthesis of the multi-copper enzyme N2O reductase, nosZ, the structural gene for this enzyme, is part of a DNA region of approximately 8 kbp that carries several essential genes. Insertional mutagenesis of the region downstream of nosZ generates apoenzyme-synthesizing strains, which argues for the existence of functions for copper acquisition or copper processing, or both, for N2O reductase. The relevant DNA region of approximately 3.2 kbp was sequenced and found to consist of three genes, nosDFY, presumably within a single transcriptional unit. Cellular concentration, copper content, and location of the N2O reductase protein were studied with mutants which were affected in the three genes. Interactions of the deduced gene products among each other and with the cytoplasmic membrane appear to be analogous to those of the components of osmotic-shock-sensitive bacterial transport systems. The 33.8-kDa product of the nosF gene belongs to the family of nucleotide-binding proteins [C. Higgins et al. (1986) Nature 323, 448-450]. Its amino acid sequence shows two canonical nucleotide-binding motifs, and the positional identity of amino acids to members of this family is around 30%. The 29.4-kDa product of the nosY gene is a hydrophobic protein with six predicted transmembrane helices and an export signal. The 48.2-kDa product of the nosD gene is a periplasmic component; it carries an export signal and is a hydrophilic protein. N2O reductase itself is a periplasmic enzyme. Our results provide evidence for an auxiliary system of Cu processing and suggest its involvement in the periplasmic biosynthesis of the Cu centers of N2O reductase.

Nucleotide sequence of the ugp genes of Escherichia coli K-12: homology to the maltose system

The nucleotide sequence of the ugp genes of Escherichia coli K-12, which encode a phosphate-limitation inducible uptake system for sn-glycerol-3-phosphate and glycerophosphoryl diesters, was determined. The genetic organization of the operon differed from previously published results. A single promoter, containing a putative pho box, was detected by S1-nuclease mapping. The promoter is followed by four open reading frames, designated ugpB, A, E and C, which encode a periplasmic binding protein, two hydrophobic membrane proteins and a protein that is likely to couple energy to the transport system, respectively. The sequences of the proteins contain the characteristics of several other binding protein-dependent transport systems, but they seem to be particularly closely related to the maltose system.

More than 150 nucleotides flanking the initiation codon contribute to the efficiency of the ribosomal binding site from bacteriophage T7 gene 1

The ribosomal binding site (RBS) from gene 1 of bacteriophage T7 was isolated on fragments of differing length and cloned upstream of the mouse dihydrofolate reductase gene to control the translation of its sequence. A 29 base pair sequence containing all elements generally believed to be essential for the RBS's showed extremely low activity. Additional upstream and downstream sequences were required to obtain a several orders of magnitude higher efficiency. By contrast, areas further downstream than +112 nucleotides from the initiator proved to be inhibitory, whereas the presence of an upstream RNaseIII cleavage site showed a strong stimulatory effect. This suggests that tertiary structures are involved in the function of the RBS studied. The efficient RBS's were complexed by ribosomes at much lower concentrations of the mRNA than the weak ones.

At least two nodD genes are necessary for efficient nodulation of alfalfa by Rhizobium meliloti

A Rhizobium meliloti DNA region (nodD1) involved in the regulation of other early nodulation genes has been delimited by directed Tn5 mutagenesis and its nucleotide sequence has been determined. The sequence data indicate a large open reading frame with opposite polarity to nodA, -B and -C, coding for a protein of 308 (or 311) amino acid residues. Tn5 insertion within the gene caused a delay in nodulation of Medicago sativa from four to seven days. Hybridization of nodD1 to total DNA of Rhizobium meliloti revealed two additional nodD sequences (nodD2 and nodD3) and both were localized on the megaplasmid pRme41b in the vicinity of the other nod genes. Genetic and DNA hybridization data, combined with nucleotide sequencing showed that nodD2 is a functional gene, while requirement of nodD3 for efficient nodulation of M. sativa could not be detected under our experimental conditions. The nodD2 gene product consists of 310 amino acid residues and shares 86.4% homology with the nodD1 protein. Single nodD2 mutants had the same nodulation phenotype as the nodD1 mutants, while a double nodD1-nodD2 mutant exhibited a more severe delay in nodulation. These results indicate that at least two functional copies of the regulatory gene nodD are necessary for the optimal expression of nodulation genes in R. meliloti.

Contextual constraints on codon pair usage: structural and biological implications

Complementary DNA sequence data of 278 protein coding genes from prokaryotic systems have been analysed at the level of near neighbour codon pairs. Our analysis points out that constraints exist even at the level of near neighbour codon pairs. These constraints are in addition to those which arise due to relative levels of tRNA. Codon pairs, which in the data base have different occurrence values from their expected values, neither have common secondary structure nor do have better stabilization due to high base stacking. Our study points out that there are strong interaction between constituent codons in these codon pairs. These strongly interacting codon pairs, we suggest, are involved in the formation of three dimensional structural elements of cDNA/mRNA and interact with ribosome and thus modulate translation.

Delineation of coding areas in DNA sequences through assignment of codon probabilities

Codon usage tables have been produced for E. coli, yeast, human, and mouse. The nonrandom employment of codons allows assignment of probability values to trinucleotides in any DNA sequence. These values represent the probability that a given trinucleotide is used as a codon in the organism from which the table is derived. For the graphical delineation of coding areas in DNA sequences, a probability is assigned to each trinucleotide equal to its frequency in the codon table. Averaging and smoothing procedures then greatly enhance the detectability of areas of high average codon probability and better represent the mean codon probability. These manipulations increase graphical clarity without altering the overall magnitude of probabilities. Averaging introduces an error of less than 0.5% between "raw" and smoothed data. This graphical delineation of coding sequences does not depend on the presence of punctuation, ribosomal binding sites, etc: moreover the delineation of introns and exons is also possible.