1
|
Erill I, O'Neill MC. A reexamination of information theory-based methods for DNA-binding site identification. BMC Bioinformatics 2009; 10:57. [PMID: 19210776 PMCID: PMC2680408 DOI: 10.1186/1471-2105-10-57] [Citation(s) in RCA: 29] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2008] [Accepted: 02/11/2009] [Indexed: 11/10/2022] Open
Abstract
Background Searching for transcription factor binding sites in genome sequences is still an open problem in bioinformatics. Despite substantial progress, search methods based on information theory remain a standard in the field, even though the full validity of their underlying assumptions has only been tested in artificial settings. Here we use newly available data on transcription factors from different bacterial genomes to make a more thorough assessment of information theory-based search methods. Results Our results reveal that conventional benchmarking against artificial sequence data leads frequently to overestimation of search efficiency. In addition, we find that sequence information by itself is often inadequate and therefore must be complemented by other cues, such as curvature, in real genomes. Furthermore, results on skewed genomes show that methods integrating skew information, such as Relative Entropy, are not effective because their assumptions may not hold in real genomes. The evidence suggests that binding sites tend to evolve towards genomic skew, rather than against it, and to maintain their information content through increased conservation. Based on these results, we identify several misconceptions on information theory as applied to binding sites, such as negative entropy, and we propose a revised paradigm to explain the observed results. Conclusion We conclude that, among information theory-based methods, the most unassuming search methods perform, on average, better than any other alternatives, since heuristic corrections to these methods are prone to fail when working on real data. A reexamination of information content in binding sites reveals that information content is a compound measure of search and binding affinity requirements, a fact that has important repercussions for our understanding of binding site evolution.
Collapse
Affiliation(s)
- Ivan Erill
- Department of Biological Sciences, University of Maryland-Baltimore County, Baltimore, MD, USA.
| | | |
Collapse
|
2
|
Perera A, Vallverdu M, Claria F, Soria JM, Caminal P. DNA binding site characterization by means of Rényi entropy measures on nucleotide transitions. IEEE Trans Nanobioscience 2008; 7:133-41. [PMID: 18556261 DOI: 10.1109/tnb.2008.2000744] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
In this work, parametric information-theory measures for the characterization of binding sites in DNA are extended with the use of transitional probabilities on the sequence. We propose the use of parametric uncertainty measures such as Rényi entropies obtained from the transition probabilities for the study of the binding sites, in addition to nucleotide frequency-based Rényi measures. Results are reported in this work comparing transition frequencies (i.e., dinucleotides) and base frequencies for Shannon and parametric Rényi entropies for a number of binding sites found in E. Coli, lambda and T7 organisms. We observe that the information provided by both approaches is not redundant. Furthermore, under the presence of noise in the binding site matrix we observe overall improved robustness of nucleotide transition-based algorithms when compared with nucleotide frequency-based method.
Collapse
Affiliation(s)
- A Perera
- Centre for Biomedical Engineering Research, Technical University of Catalonia, Barcelona, Spain.
| | | | | | | | | |
Collapse
|
3
|
Abstract
In contrast to the vast majority of the members of the domain Bacteria, several Pseudomonas and Xanthomonas species have two lexA genes, whose products have been shown to recognize different LexA binding motifs, making them an interesting target for studying the interplay between cohabiting LexA regulons in a single species. Here we report an analysis of the genetic composition of the two LexA regulons of Pseudomonas putida KT2440 performed with a genomic microarray. The data obtained indicate that one of the two LexA proteins (LexA1) seems to be in control of the conventional Escherichia coli-like SOS response, while the other LexA protein (LexA2) regulates only its own transcriptional unit, which includes the imuA, imuB, and dnaE2 genes, and a gene (PP_3901) from a resident P. putida prophage. Furthermore, PP_3901 is also regulated by LexA1 and is required for DNA damage-mediated induction of several P. putida resident prophage genes. In silico searches suggested that this marked asymmetry in regulon contents also occurs in other Pseudomonas species with two lexA genes, and the implications of this asymmetry in the evolution of the SOS network are discussed.
Collapse
|
4
|
Abstract
Information theory was used to build a promoter model that accounts for the -10, the -35 and the uncertainty of the gap between them on a common scale. Helical face assignment indicated that base -7, rather than -11, of the -10 may be flipping to initiate transcription. We found that the sequence conservation of sigma70 binding sites is 6.5 +/- 0.1 bits. Some promoters lack a -35 region, but have a 6.7 +/- 0.2 bit extended -10, almost the same information as the bipartite promoter. These results and similarities between the contacts in the extended -10 binding and the -35 suggest that the flexible bipartite sigma factor evolved from a simpler polymerase. Binding predicted by the bipartite model is enriched around 35 bases upstream of the translational start. This distance is the smallest 5' mRNA leader necessary for ribosome binding, suggesting that selective pressure minimizes transcript length. The promoter model was combined with models of the transcription factors Fur and Lrp to locate new promoters, to quantify promoter strengths, and to predict activation and repression. Finally, the DNA-bending proteins Fis, H-NS and IHF frequently have sites within one DNA persistence length from the -35, so bending allows distal activators to reach the polymerase.
Collapse
Affiliation(s)
| | | | | | - Thomas D. Schneider
- To whom correspondence should be addressed. Tel: +1 301 846 5581; Fax: +1 301 846 5598;
| |
Collapse
|
5
|
Elloumi M, Maddouri M. New voting strategies designed for the classification of nucleic sequences. Knowl Inf Syst 2005. [DOI: 10.1007/s10115-004-0151-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
6
|
Kiryu H, Oshima T, Asai K. Extracting relations between promoter sequences and their strengths from microarray data. Bioinformatics 2004; 21:1062-8. [PMID: 15513998 DOI: 10.1093/bioinformatics/bti094] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
MOTIVATION The relations between the promoter sequences and their strengths were extensively studied in the 1980s. Although these studies uncovered strong sequence-strength correlations, the cost of their elaborate experimental methods have been too high to be applied to a large number of promoters. On the contrary, a recent increase in the microarray data allows us to compare thousands of gene expressions with their DNA sequences. RESULTS We studied the relations between the promoter sequences and their strengths using the Escherichia coli microarray data. We modeled those relations using a simple weight matrix, which was optimized with a novel support vector regression method. It was observed that several non-consensus bases in the '-35' and '-10' regions of promoter sequences act positively on the promoter strength and that certain consensus bases have a minor effect on the strength. We analyzed outliers for which the observed gene expressions deviate from the promoter strength predictions, and identified several genes with enhanced expressions due to multiple promoters and genes under strong regulation by transcription factors. Our method is applicable to other procaryotes for which both the promoter sequences and the microarray data are available.
Collapse
Affiliation(s)
- Hisanori Kiryu
- Graduate School of Information Sciences, Nara Institute of Science and Technology 8916-5 Takayama-cho, Ikoma, Nara 630-0192, Japan.
| | | | | |
Collapse
|
7
|
O'Neill MC. A Procedure for Identifying Loosely Conserved Protein-Binding DNA Sequences. Methods Enzymol 2003; 370:229-37. [PMID: 14712648 DOI: 10.1016/s0076-6879(03)70020-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Michael C O'Neill
- Department of Biological Sciences, University of Maryland, Baltimore County (UMBC), 1000 Hilltop Circle, Baltimore, Maryland 21250, USA
| |
Collapse
|
8
|
Munson GP, Holcomb LG, Scott JR. Novel group of virulence activators within the AraC family that are not restricted to upstream binding sites. Infect Immun 2001; 69:186-93. [PMID: 11119505 PMCID: PMC97871 DOI: 10.1128/iai.69.1.186-193.2001] [Citation(s) in RCA: 51] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022] Open
Abstract
Several regulators within the AraC family control the expression of genes known or thought to be required for virulence of bacterial pathogens. One of these, Rns, activates transcription from an unprecedented variety of binding-site locations. Although nearly all prokaryotic activators bind within a small region upstream and adjacent to the promoter that they regulate, Rns does not bind within this region to activate its own promoter, Prns. Instead, to activate Prns, Rns requires one binding site 224.5 bp upstream and one downstream of the transcription start site. We show in this study that several other AraC family activators recognize the same binding sites as Rns and share with it the ability to utilize a downstream binding site. Like Rns, other members of this group of activators positively regulate the expression of virulence factors in pathogenic bacteria. These regulators are also able to activate transcription from promoter-proximal upstream binding sites since they are able to substitute for Rns at Pcoo, a promoter with only upstream binding sites. Thus, Rns is the prototype for a group of regulators, which include CfaR, VirF, AggR, and CsvR and which activate transcription from locations that are more diverse than those of any other known activator.
Collapse
Affiliation(s)
- G P Munson
- Department of Microbiology and Immunology, Emory University School of Medicine, Atlanta, Georgia 30322, USA
| | | | | |
Collapse
|
9
|
Abstract
This paper presents a survey of currently available mathematical models and algorithmical methods for trying to identify promoter sequences. The methods concern both searching in a genome for a previously defined consensus and extracting a consensus from a set of sequences. Such methods were often tailored for either eukaryotes or prokaryotes although this does not preclude use of the same method for both types of organisms. The survey therefore covers all methods; however, emphasis is placed on prokaryotic promoter sequence identification. Illustrative applications of the main extracting algorithms are given for three bacteria.
Collapse
Affiliation(s)
- A Vanet
- Institut de biologie physico-chimique, Paris, France
| | | | | |
Collapse
|
10
|
Abstract
Availability of complete bacterial genomes opens the way to the comparative approach to the recognition of transcription regulatory sites. Assumption of regulon conservation in conjunction with profile analysis provides two lines of independent evidence making it possible to make highly specific predictions. Recently this approach was used to analyze several regulons in eubacteria and archaebacteria. The present review covers recent advances in the comparative analysis of transcriptional regulation in prokaryotes and phylogenetic fingerprinting techniques in eukaryotes, and describes the emerging patterns of the evolution of regulatory systems.
Collapse
Affiliation(s)
- M S Gelfand
- State Scientific Center for Biotechnology 'NIIGenetika', Moscow, Russia.
| |
Collapse
|
11
|
Wang JT, Rozen S, Shapiro BA, Shasha D, Wang Z, Yin M. New techniques for DNA sequence classification. J Comput Biol 1999; 6:209-18. [PMID: 10421523 DOI: 10.1089/cmb.1999.6.209] [Citation(s) in RCA: 22] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
DNA sequence classification is the activity of determining whether or not an unlabeled sequence S belongs to an existing class C. This paper proposes two new techniques for DNA sequence classification. The first technique works by comparing the unlabeled sequence S with a group of active motifs discovered from the elements of C and by distinction with elements outside of C. The second technique generates and matches gapped fingerprints of S with elements of C. Experimental results obtained by running these algorithms on long and well conserved Alu sequences demonstrate the good performance of the presented methods compared with FASTA. When applied to less conserved and relatively short functional sites such as splice-junctions, a variation of the second technique combining fingerprinting with consensus sequence analysis gives better results than the current classifiers employing text compression and machine learning algorithms.
Collapse
Affiliation(s)
- J T Wang
- Department of Computer and Information Science, New Jersey Institute of Technology, University Heights, Newark 07102, USA.
| | | | | | | | | | | |
Collapse
|
12
|
Robison K, McGuire AM, Church GM. A comprehensive library of DNA-binding site matrices for 55 proteins applied to the complete Escherichia coli K-12 genome. J Mol Biol 1998; 284:241-54. [PMID: 9813115 DOI: 10.1006/jmbi.1998.2160] [Citation(s) in RCA: 263] [Impact Index Per Article: 10.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
A major mode of gene regulation occurs via the binding of specific proteins to specific DNA sequences. The availability of complete bacterial genome sequences offers an unprecedented opportunity to describe networks of such interactions by correlating existing experimental data with computational predictions. Of the 240 candidate Escherichia coli DNA-binding proteins, about 55 have DNA-binding sites identified by DNA footprinting. We used these sites to construct recognition matrices, which we used to search for additional binding sites in the E. coli genomic sequence. Many of these matrices show a strong preference for non-coding DNA. Discrepancies are identified between matrices derived from natural sites and those derived from SELEX (Systematic Evolution of Ligands by Exponential enrichment) experiments. We have constructed a database of these proteins and binding sites, called DPInteract (available at http://arep.med.harvard.edu/dpinteract).
Collapse
Affiliation(s)
- K Robison
- Department of Genetics, Harvard Medical School, Boston, MA 02115, USA
| | | | | |
Collapse
|
13
|
O'Neill MC. A general procedure for locating and analyzing protein-binding sequence motifs in nucleic acids. Proc Natl Acad Sci U S A 1998; 95:10710-5. [PMID: 9724769 PMCID: PMC27960 DOI: 10.1073/pnas.95.18.10710] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
In the last decade, two tools, one drawn from information theory and the other from artificial neural networks, have proven particularly useful in many different areas of sequence analysis. The work presented herein indicates that these two approaches can be joined in a general fashion to produce a very powerful search engine that is capable of locating members of a given nucleic acid sequence family in either local or global sequence searches. This program can, in turn, be queried for its definition of the motif under investigation, ranking each base in context for its contribution to membership in the motif family. In principle, the method used can be applied to any binding motif, including both DNA and RNA sequence families, given sufficient family size.
Collapse
Affiliation(s)
- M C O'Neill
- Department of Biological Sciences, University of Maryland, Baltimore County, 1000 Hilltop Circle, Baltimore, MD 21250, USA
| |
Collapse
|
14
|
Marolda CL, Valvano MA. Promoter region of the Escherichia coli O7-specific lipopolysaccharide gene cluster: structural and functional characterization of an upstream untranslated mRNA sequence. J Bacteriol 1998; 180:3070-9. [PMID: 9620955 PMCID: PMC107806 DOI: 10.1128/jb.180.12.3070-3079.1998] [Citation(s) in RCA: 44] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
We report the identification of the promoter region of the Escherichia coli O7-specific lipopolysaccharide (LPS) gene cluster (wbEcO7). Typical -10 and -35 sequences were found to be located in the intervening region between galF and rlmB, the first gene of the wbEcO7 cluster. Data from RNase protection experiments revealed the existence of an untranslated leader mRNA segment of 173 bp, including the JUMPStart and two ops sequences. We characterized the structure of this leader mRNA by using the program Mfold and a combination of nested and internal deletions transcriptionally fused to a promoterless lac operon. Our results indicated that the leader mRNA may fold into a series of complex stem-loop structures, one of which includes the JUMPStart element. We have also found that one of the ops sequences resides on the predicted stem and the other resides on the loop region, and we confirmed that these sequences are essential for the RfaH-mediated regulation of the O polysaccharide cluster. A very similar stem-loop structure could be predicted in the promoter region of the LPS core operon encoding the waaQGPSBIJYZK genes. We observed another predicted stem-loop, located immediately downstream from the wbEcO7 transcription initiation site, which appeared to be involved in premature termination of transcription. This putative stem-loop is common to many other O polysaccharide gene clusters but is not present in core oligosaccharide genes. wbEcO7-lac transcriptional fusions in single copy numbers were also used to determine the effects of various environmental cues in the transcriptional regulation of O polysaccharide synthesis. No effects were detected with temperature, osmolarity, Mg2+ concentration, and drugs inducing changes in DNA supercoiling. We therefore conclude that the wbEcO7 promoter activity may be constitutive and that regulation takes place at the level of elongation of the mRNA in a RfaH-mediated manner.
Collapse
Affiliation(s)
- C L Marolda
- Department of Microbiology and Immunology, University of Western Ontario, London, Canada
| | | |
Collapse
|
15
|
Zhu Q, Zhao S, Somerville RL. Expression, purification, and functional analysis of the TyrR protein of Haemophilus influenzae. Protein Expr Purif 1997; 10:237-46. [PMID: 9226720 DOI: 10.1006/prep.1997.0757] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
The gene that was inferred to encode the TyrR protein of Haemophilus influenzae Rd was synthesized by polymerase chain reaction and inserted into a T7-based expression vector. Methods were developed to overexpress the TyrR protein of H. influenzae in Escherichia coli and to purify the protein on a large scale. Both in vitro and in vivo functional comparisons of the H. influenzae and E. coli TyrR proteins were carried out. The TyrR protein of H. influenzae was able to bind in vitro to an operator target upstream of the aroF-tyrA gene of E. coli. In the presence of [gamma-S]ATP, the DNA binding ability of the H. influenzae TyrR protein was drastically reduced. Despite the much shorter peptide chain length (318 amino acid residues vs 513), the TyrR protein of H. influenzae was as active in repressing the aroF promoter as the TyrR protein of E. coli. Repression by both proteins was enhanced in the presence of tyrosine; however, the transcriptional activation function associated with the TyrR protein of E. coli could not be detected when the H. influenzae TyrR protein was expressed in E. coli. By computer analysis, at least five operator targets for TyrR were identified within the genomic DNA of H. influenzae. These observations show that the assignment of function to the tyrR gene of H. influenzae was correctly made. Further studies of the H. influenzae TyrR protein in comparison to its E. coli counterpart should provide valuable mechanistic information on transcriptional regulation in this system.
Collapse
Affiliation(s)
- Q Zhu
- Department of Biochemistry, Purdue University, West Lafayette, Indiana 47907, USA
| | | | | |
Collapse
|
16
|
Marolda CL, Valvano MA. Genetic analysis of the dTDP-rhamnose biosynthesis region of the Escherichia coli VW187 (O7:K1) rfb gene cluster: identification of functional homologs of rfbB and rfbA in the rff cluster and correct location of the rffE gene. J Bacteriol 1995; 177:5539-46. [PMID: 7559340 PMCID: PMC177362 DOI: 10.1128/jb.177.19.5539-5546.1995] [Citation(s) in RCA: 82] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023] Open
Abstract
The O-repeating unit of the Escherichia coli O7-specific lipopolysaccharide is made of galactose, mannose, rhamnose, 4-acetamido-4,6-dideoxyglucose, and N-acetyglucosamine. We have recently characterized the genes involved in the biosynthesis of the sugar precursor GDP-mannose occurring in the E. coli O7:K1 strain VW187 (C. L. Marolda and M. A. Valvano, J. Bacteriol. 175:148-158, 1993). In the present study, we identified and sequenced the rfbBDAC genes encoding the enzymes for the biosynthesis of another precursor, dTDP-rhamnose. These genes are localized on the upstream end of the rfbEcO7 region, and they are strongly conserved compared with similar genes found in various enteric and nonenteric bacteria. Upstream of rfbB we identified a DNA segment containing the rfb promoter and a highly conserved untranslated leader sequence also present in the promoter regions of other surface polysaccharide gene clusters. Also, we have determined that rfbB and rfbA have homologs, rffG (o355) and rffH (o292), respectively, located on the rff cluster, which is involved in the synthesis of enterobacterial common antigen. We provide biochemical evidence that rffG and rffH encode dTDP-glucose dehydratase and glucose-1-phosphate thymidylyltransferase activities, respectively, and we also show that rffG complemented the rfbB defect in the O7+ cosmid pJHCV32. We also demonstrate that rffG is distinct from rffE and map the rffE gene to the second gene of the rff cluster.
Collapse
Affiliation(s)
- C L Marolda
- Department of Microbiology and Immunology, University of Western Ontario, London, Canada
| | | |
Collapse
|
17
|
Abstract
Lrp (leucine-responsive regulatory protein) is a major regulatory protein involved in the expression of numerous operons in Escherichia coli. For ilvIH, one of the operons positively regulated by Lrp, Lrp binds to multiple sites upstream of the transcriptional start site and activates transcription. An alignment of 12 Lrp binding sites within ilvIH DNA from two different organisms revealed a tentative consensus sequence AGAAT TTTATTCT (Q. Wang, M. Sacco, E. Ricca, C.T. Lago, M. DeFelice, and J.M. Calvo, Mol. Microbiol. 7:883-891, 1993). To further characterize the binding specificity of Lrp, we used a variation of the Selex procedure of C. Tuerk and L. Gold (Science 249:505-510, 1990) to identify sequences that bound Lrp out of a pool of 10(12) different DNA molecules. We identified 63 related DNA sequences that bound Lrp and estimated their relative binding affinities for Lrp. A consensus sequence derived from analysis of these sequences, YAGHAWATTWT DCTR, where Y = C or T, H = not G, W = A or T, D = not C, and R = A or G, contains clear dyad symmetry and is very similar to the one defined earlier. To test the idea that Lrp in the presence of leucine might bind to a different subset of DNA sequences, we carried out a second selection experiment with leucine present during the binding reactions. DNA sequences selected in the presence or absence of leucine were similar, and leucine did not stimulate binding to any of the sequences that were selected in the presence of leucine. Therefore, it is unlikely that leucine changes the specificity of Lrp binding.
Collapse
Affiliation(s)
- Y Cui
- Section of Biochemistry, Molecular and Cell Biology, Cornell University, Ithaca, New York 14853, USA
| | | | | | | |
Collapse
|
18
|
Holberton DV, Marshall J. Analysis of consensus sequence patterns in Giardia cytoskeleton gene promoters. Nucleic Acids Res 1995; 23:2945-53. [PMID: 7659516 PMCID: PMC307134 DOI: 10.1093/nar/23.15.2945] [Citation(s) in RCA: 50] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
Abstract
Protein-coding genes in the ancient eukaryote Giardia lamblia lack typical promoter consensus elements. We have analysed the immediate 5' flanking sequences of seven genes of related function (structural cytoskeleton proteins) to identify shared DNA motifs that might have a role in transcription initiation. Transcription start sites for five genes have been determined previously. Genomic mapping and mRNA primer extension experiments demonstrate additionally that the genes for beta-giardin and median body protein are (i) present as single copies in the genome, (ii) transcribed with very short 5' leader sequences. Two search algorithms designed to extract conserved motifs from either aligned or non-aligned sequences independently discovered three sites constituting a common pattern in all seven promoters. Sites were optimally aligned using weight matrix building trials to achieve the maximum 'information content'. Profiling the information content of best alignments defines the extent of the homologies as: a 9 bp box (initiator) at the start site and upstream 18 and 6 bp boxes. The initiator is the most highly conserved element and contains a universal Py-A-Pu motif at which transcription starts. We show that the best matrices can be combined in a search pattern that correctly locates transcription start sites in genomic DNA sequences.
Collapse
Affiliation(s)
- D V Holberton
- Department of Life Science, Nottingham University, UK
| | | |
Collapse
|
19
|
Nalin R, Domenach AM, Normand P. Molecular structure of the Frankia spp. nifD-K intergenic spacer and design of Frankia genus compatible primer. Mol Ecol 1995; 4:483-91. [PMID: 8574444 DOI: 10.1111/j.1365-294x.1995.tb00242.x] [Citation(s) in RCA: 23] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
The nifD-K intergenic spacer (IGS) of ArI3 and ACoN24d were found to have a length 265 and 199 nucleotides, respectively. They are markedly less conserved than the two neighbouring genes and have, in some instances, a repeated structure reminiscent of an insertion event. The repeated sequence and the IGSs have no detectable homology with sequences in DNA databanks. The IGS has a stem-loop structure with a low folding energy, lower than that between nifH and nifD. No convincing alignment of IGS sequences could be obtained among Frankia strains. Only between ACoN24d and ArI3, which belong to the same genomic species, was the alignment good enough to permit detection of a doubly repeated structure. No promoter could be detected in the IGSs. The putative nifK open reading frame (ORF) in Frankia strain ArI3 has a length of 1587 nucleotides, starting with a GTG codon, preceded by a ribosome binding site of a structure similar to that of nifH (GGAGGN7). The codon usage was similar to that of previously sequenced Frankia genes with a strong bias toward G- and C-ending codons except in the case of glycine where GGT is frequent. Alignment of the three Frankia nifK sequences (EUN1f; ArI3 and ACoN24d) with those of other nitrogen-fixing bacteria permitted detection of a sequence conserved among the three Frankia strains but absent in the other sequences. A primer targeted to that region in combination with FGPD807-85 amplified the nifD-KIGS sequences of all Frankia strains (except the non-nitrogen-fixing Frankia strains CN3 and AgB1-9) and yet failed to amplify DNA of all other nitrogen-fixing bacteria.(ABSTRACT TRUNCATED AT 250 WORDS)
Collapse
Affiliation(s)
- R Nalin
- Laboratoire d'Ecologie Microbienne du Sol, URA CNRS 1977, Université Claude Bernard Lyon I, Villeurbanne, France
| | | | | |
Collapse
|
20
|
Bisant D, Maizel J. Identification of ribosome binding sites in Escherichia coli using neural network models. Nucleic Acids Res 1995; 23:1632-9. [PMID: 7784221 PMCID: PMC306908 DOI: 10.1093/nar/23.9.1632] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023] Open
Abstract
This study investigated the use of neural networks in the identification of Escherichia coli ribosome binding sites. The recognition of these sites based on primary sequence data is difficult due to the multiple determinants that define them. Additionally, secondary structure plays a significant role in the determination of the site and this information is difficult to include in the models. Efforts to solve this problem have so far yielded poor results. A new compilation of E. coli ribosome binding sites was generated for this study. Feedforward backpropagation networks were applied to their identification. Perceptrons were also applied, since they have been the previous best method since 1982. Evaluation of performance for all the neural networks and perceptrons was determined by ROC analysis. The neural network provided significant improvement in the recognition of these sites when compared with the previous best method, finding less than half the number of false positives when both models were adjusted to find an equal number of actual sites. The best neural network used an input window of 101 nucleotides and a single hidden layer of 9 units. Both the neural network and the perceptron trained on the new compilation performed better than the original perceptron published by Stormo et al. in 1982.
Collapse
Affiliation(s)
- D Bisant
- Neuroscience Program (151 B), Stanford University, CA 94305, USA
| | | |
Collapse
|
21
|
Abstract
Recognition of function of newly sequenced DNA fragments is an important area of computational molecular biology. Here we present an extensive review of methods for prediction of functional sites, tRNA, and protein-coding genes and discuss possible further directions of research in this area.
Collapse
Affiliation(s)
- M S Gelfand
- Institute of Protein Research, Russian Academy of Sciences, Pushchino, Moscow region, Russia
| |
Collapse
|
22
|
Titgemeyer F, Reizer J, Reizer A, Tang J, Parr TR, Saier MH. Nucleotide sequence of the region between crr and cysM in Salmonella typhimurium: five novel ORFs including one encoding a putative transcriptional regulator of the phosphotransferase system. DNA SEQUENCE : THE JOURNAL OF DNA SEQUENCING AND MAPPING 1995; 5:145-52. [PMID: 7612925 DOI: 10.3109/10425179509029354] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
A 4471 bp region between crr and cysM on the Salmonella typhimurium chromosome (49.5 min) has been sequenced. Five ORFs were found within this region, one of which is likely to be the putative regulatory gene, ptsJ, that corresponds in map position to a gene which when mutated allows expression of a cryptic Enzyme I of the phosphotransferase system. The deduced amino acid sequence of the encoded protein is similar to those of several open reading frames (ORFs) including ORFT2 of Rhodobacter spheroides with which it is 28% identical throughout most of its length (comparison score of 21 S.D.). PtsJ exhibits a putative, N-terminal, helix-turn-helix, DNA binding domain that is similar in sequence to those in members of the GntR family of transcriptional regulators. Analyses of the sequences of the ORFs encoded within this region are presented.
Collapse
Affiliation(s)
- F Titgemeyer
- Department of Biology, University of California at San Diego, La Jolla 92093-0116, USA
| | | | | | | | | | | |
Collapse
|
23
|
Grahn AM, Bamford JK, O'Neill MC, Bamford DH. Functional organization of the bacteriophage PRD1 genome. J Bacteriol 1994; 176:3062-8. [PMID: 8188607 PMCID: PMC205464 DOI: 10.1128/jb.176.10.3062-3068.1994] [Citation(s) in RCA: 23] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023] Open
Abstract
PRD1 is a broad-host-range virus that infects Escherichia coli cells. It has a linear double-stranded DNA genome that replicates by a protein-primed mechanism. The virus particle is composed of a protein coat enclosing a lipid membrane. On the basis of this structure, PRD1 is being used as a membrane biosynthesis and structure model. In this investigation, we constructed the transcription map of the 15-kb-long phage genome. This was achieved by a computer search of putative promoters, which were then tested for activity by primer extension and for the capability to promote the synthesis of chloramphenicol acetyltransferase.
Collapse
Affiliation(s)
- A M Grahn
- Department of Genetics, University of Helsinki, Finland
| | | | | | | |
Collapse
|
24
|
Silk GW, Wu M. Posttranscriptional accumulation of chloroplast tufA (elongation factor gene) mRNA during chloroplast development in Chlamydomonas reinhardtii. PLANT MOLECULAR BIOLOGY 1993; 23:87-96. [PMID: 8219059 DOI: 10.1007/bf00021422] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/22/2023]
Abstract
Light induces chloroplast (Cp) differentiation in dark-grown y-1 strains of Chlamydomonas reinhardtii. Slot blot analysis was used to quantitate tufA, psbA, psbK, rbcL, and 16S rRNA transcript accumulation and transcription during Cp differentiation. When etiolated cc-125 y-1 cells were illuminated for 5 h, a 1710 bp tufA mRNA accumulated up to 5-fold while the psbA, rbcL, and 16S rRNA transcripts accumulated less than 1.5-fold. The tufA gene encodes translational elongation factor EF-Tu. The light-induced accumulation of tufA mRNA did not occur in cc-1931, a strain that does not become etiolated in darkness. Pulse labelling was used to measure the transcription of Cp transcripts during tufA mRNA accumulation, and no detectable change in tufA transcription was observed. These results imply that the half life of the tufA transcript increases during the greening process.
Collapse
Affiliation(s)
- G W Silk
- Department of Biological Sciences, University of Maryland Baltimore County 21228
| | | |
Collapse
|
25
|
Abstract
An updated compilation of 300 E. coli mRNA promoter sequences is presented. For each sequence the most recent relevant paper was checked, to verify the location of the transcriptional start position as identified experimentally. We comment on the reliability of the sequence databanks and analyze the conservation of known promoter features in the current compilation. This database is available by E-mail.
Collapse
Affiliation(s)
- S Lisser
- Department of Molecular Genetics, Hebrew University-Hadassah Medical School, Jerusalem, Israel
| | | |
Collapse
|
26
|
Kessler AC, Haase A, Reeves PR. Molecular analysis of the 3,6-dideoxyhexose pathway genes of Yersinia pseudotuberculosis serogroup IIA. J Bacteriol 1993; 175:1412-22. [PMID: 8444803 PMCID: PMC193228 DOI: 10.1128/jb.175.5.1412-1422.1993] [Citation(s) in RCA: 52] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023] Open
Abstract
Salmonella enterica and Yersinia pseudotuberculosis are the only examples in nature known to use a variety of 3,6-dideoxyhexose derivatives as O antigen constituents. To allow a comparison of the responsible biosynthetic genes of the two organisms, we have sequenced a section of the Y. pseudotuberculosis serogroup IIA rfb region that contained the genes for the abequose biosynthetic pathway. Comparison of the identified genes with the rfb region of S. enterica LT2 showed that the two dideoxyhexose pathway gene clusters are related. The arrangement of the genes was largely conserved, and the G + C compositions of the two DNA regions were strikingly similar; however, the degree of conservation of nucleotide and protein sequences suggested that the two gene clusters have been evolving independently for considerable time. Hybridization experiments showed that the dideoxyhexose pathway genes are widespread throughout the various serogroups of Y. pseudotuberculosis.
Collapse
Affiliation(s)
- A C Kessler
- Department of Microbiology, University of Sydney, New South Wales, Australia
| | | | | |
Collapse
|
27
|
Horton PB, Kanehisa M. An assessment of neural network and statistical approaches for prediction of E. coli promoter sites. Nucleic Acids Res 1992; 20:4331-8. [PMID: 1508724 PMCID: PMC334144 DOI: 10.1093/nar/20.16.4331] [Citation(s) in RCA: 54] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
We have constructed a perceptron type neural network for E. coli promoter prediction and improved its ability to generalize with a new technique for selecting the sequence features shown during training. We have also reconstructed five previous prediction methods and compared the effectiveness of those methods and our neural network. Surprisingly, the simple statistical method of Mulligan et al. performed the best amongst the previous methods. Our neural network was comparable to Mulligan's method when false positives were kept low and better than Mulligan's method when false negatives were kept low. We also showed the correlation between the prediction rates of neural networks achieved by previous researchers and the information content of their data sets.
Collapse
Affiliation(s)
- P B Horton
- Institute for Chemical Research, Kyoto University, Japan
| | | |
Collapse
|
28
|
O'Neill MC. Escherichia coli promoters: neural networks develop distinct descriptions in learning to search for promoters of different spacing classes. Nucleic Acids Res 1992; 20:3471-7. [PMID: 1630917 PMCID: PMC312504 DOI: 10.1093/nar/20.13.3471] [Citation(s) in RCA: 40] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
Back-propagation neural networks were trained to recognize promoter sequences of each of the three major spacing classes found in E. coli. These networks were trained with the object of maximizing their ability to generalize while maintaining the level of false positive identifications at a fraction of 1 percent. These objectives were generally met. Networks for the 16 base spacing class captured between 78 and 100% of previously unseen promoters in different tests; networks for the 17 base class identified 97% of the test promoters; networks for the 18 base class identified 79% of the test promoters. A tandem poll of networks for all three spacing classes produced a cumulative false positive level of less than 0.5%. In each case, the weight matrices used by the networks in their classification were analyzed to determine the relative weight assigned to the occurrence of a given base at a given position within the promoter. In this fashion, an approximate description of the network's definition of the promoter can be obtained.
Collapse
Affiliation(s)
- M C O'Neill
- Department of Biological Sciences, University of Maryland-Baltimore County 21228
| |
Collapse
|
29
|
Lewis LK, Jenkins ME, Mount DW. Isolation of DNA damage-inducible promoters in Escherichia coli: regulation of polB (dinA), dinG, and dinH by LexA repressor. J Bacteriol 1992; 174:3377-85. [PMID: 1577702 PMCID: PMC206008 DOI: 10.1128/jb.174.10.3377-3385.1992] [Citation(s) in RCA: 72] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
A new genetic screening method has been developed to isolate Escherichia coli promoters which are components of the SOS regulon. Plasmids containing the regulatory regions of polB (dinA) and two new loci, dinG and dinH, were characterized. Galactokinase gene fusion experiments indicated that transcription of these genes is inducible by treatment with mitomycin and conforms to a classical model of SOS regulation involving simple LexA repression. Mapping studies using the E. coli DNA library of Kohara et al. (Y. Kohara, K. Akiyama, and K. Isono, Cell 50:495-508, 1987) revealed that dinG and dinH are located at 17.8 and 19.8 min on the chromosome, respectively. The nucleotide sequence of the dinH regulatory region contains a segment which is very similar to previously characterized binding sites for LexA protein. An asymmetric, noncanonical 20-bp LexA operator in the cloned dinG promoter region was identified. Additional experiments have revealed that the nucleotide sequence of the gene immediately downstream of the DNA damage-inducible polB locus encodes a polypeptide which has extensive sequence homology to several known and putative DNA and RNA helicase proteins. This gene, which is not regulated by the LexA repressor, has been designated hepA. The predicted amino acid sequence of the product of hepA contains several highly conserved sequence motifs that are also found in enzymes such as the RecQ and UvrB proteins of E. coli and the Rad3 protein of Saccharomyces cerevisiae.
Collapse
Affiliation(s)
- L K Lewis
- Molecular and Cellular Biology Department, University of Arizona, Tucson 85721
| | | | | |
Collapse
|
30
|
Cardon LR, Stormo GD. Expectation maximization algorithm for identifying protein-binding sites with variable lengths from unaligned DNA fragments. J Mol Biol 1992; 223:159-70. [PMID: 1731067 DOI: 10.1016/0022-2836(92)90723-w] [Citation(s) in RCA: 97] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
An Expectation Maximization algorithm for identification of DNA binding sites is presented. The approach predicts the location of binding regions while allowing variable length spacers within the sites. In addition to predicting the most likely spacer length for a set of DNA fragments, the method identifies individual sites that differ in spacer size. No alignment of DNA sequences is necessary. The method is illustrated by application to 231 Escherichia coli DNA fragments known to contain promoters with variable spacings between their consensus regions. Maximum-likelihood tests of the differences between the spacing classes indicate that the consensus regions of the spacing classes are not distinct. Further tests suggest that several positions within the spacing region may contribute to promoter specificity.
Collapse
Affiliation(s)
- L R Cardon
- Institute for Behavior Genetics, University of Colorado Boulder 80309-0447
| | | |
Collapse
|
31
|
Rowley DL, Fawcett WP, Wolf RE. Molecular characterization of mutations affecting expression level and growth rate-dependent regulation of the Escherichia coli zwf gene. J Bacteriol 1992; 174:623-6. [PMID: 1729252 PMCID: PMC205758 DOI: 10.1128/jb.174.2.623-626.1992] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
We characterized three cis dominant mutations which elevate glucose 6-phosphate dehydrogenase level. Growth rate-dependent regulation and oxidative stress control of enzyme level were altered by the mutations. DNA sequencing and transcript mapping showed that the "up" mutations created new promoters whose hyperactive expression overrides the normal regulation of the native promoter.
Collapse
Affiliation(s)
- D L Rowley
- Department of Biological Sciences, University of Maryland Baltimore County, Catonsville 21228-5398
| | | | | |
Collapse
|
32
|
Gartmann CJ, Grob U. SQUIRREL: Sequence QUery, Information Retrieval and REporting Library. A program package for analyzing signals in nucleic acid sequences for the VAX. Nucleic Acids Res 1991; 19:6033-40. [PMID: 1945887 PMCID: PMC329063 DOI: 10.1093/nar/19.21.6033] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
A computer tool is described for comparison, analysis and search of genetic signals. The method is based on sequence consensus matrices. It assumes that a genetic signal (such as a promoter, enhancer or whatever) is composed of several signal blocks separated from each other by variable distances. A set of programs is presented to perform the analysis. The result of such an analysis is a description of the investigated signal including matrices for each signal block, distances between each block and distribution of the values. Programs are provided to search for a signal using results from previous analysis. The method is able to align large sets of sequences within a few minutes and to check the quality of the alignment. An analysis of E.coli promoters is provided as an example.
Collapse
Affiliation(s)
- C J Gartmann
- Max-Planck-Institut für Immunbiologie, Freiburg, FRG
| | | |
Collapse
|
33
|
Vartak NB, Reizer J, Reizer A, Gripp JT, Groisman EA, Wu LF, Tomich JM, Saier MH. Sequence and evolution of the FruR protein of Salmonella typhimurium: a pleiotropic transcriptional regulatory protein possessing both activator and repressor functions which is homologous to the periplasmic ribose-binding protein. Res Microbiol 1991; 142:951-63. [PMID: 1805309 DOI: 10.1016/0923-2508(91)90005-u] [Citation(s) in RCA: 34] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
The repressor of the fructose (fru) operon of Salmonella typhimurium (FruR) has been implicated in the transcriptional regulation of dozens of genes concerned with central metabolic pathways of carbon utilization. We here report the nucleotide sequence of the gene encoding FruR and analyse both its operator-promoter region and its deduced amino acyl sequence. The FruR protein was overexpressed and was shown to have a molecular weight of about 36 kDa in agreement with the molecular weight deduced from the gene sequence. Sequence analyses revealed that FruR is homologous to 9 distinct bacterial DNA-binding proteins, most of which recognize sugar inducers and all of which possess helix-turn-helix motifs within their N-terminal regions and exhibit sequence identity throughout most of their lengths. FruR is also homologous to the periplasmic ribose-binding protein which serves as a constituent of the ribose transport/chemoreception system. The ribose-binding protein is in turn homologous to binding proteins specific for arabinose and galactose. The periplasmic binding proteins, the structures of some of which have been elucidated in three dimensions, lack the N-terminal helix-turn-helix region, but instead possess N-terminal hydrophobic signal sequences which target them to the periplasm. A phylogenetic tree for the more closely related proteins of this superfamily was constructed, and a signature motif was identified which should facilitate future detection of additional transcriptional regulatory proteins belonging to this family.
Collapse
Affiliation(s)
- N B Vartak
- Department of Biology, University of California at San Diego, La Jolla 92093-0116
| | | | | | | | | | | | | | | |
Collapse
|
34
|
Wu LF, Reizer A, Reizer J, Cai B, Tomich JM, Saier MH. Nucleotide sequence of the Rhodobacter capsulatus fruK gene, which encodes fructose-1-phosphate kinase: evidence for a kinase superfamily including both phosphofructokinases of Escherichia coli. J Bacteriol 1991; 173:3117-27. [PMID: 1850730 PMCID: PMC207905 DOI: 10.1128/jb.173.10.3117-3127.1991] [Citation(s) in RCA: 70] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
The fruK gene encoding fructose-1-phosphate kinase (FruK), located within the fructose (fru)-catabolic operon of Rhodobacter capsulatus, was sequenced. FruK of R. capsulatus (316 amino acids; molecular weight = 31,232) is the same size as and is homologous to FruK of Escherichia coli, phosphofructokinase B (PfkB) of E. coli, phosphotagatokinase of Staphylococcus aureus, and ribokinase of E. coli. These proteins therefore make up a family of homologous proteins, termed the PfkB family. A phylogenetic tree for this new family was constructed. Sequence comparisons plus chemical inactivation studies suggested the lack of involvement of specific residues in catalysis. Although the Rhodobacter FruK differed markedly from the other enzymes within the PfkB family with respect to amino acid composition, these enzymes exhibited similar predicted secondary structural features. A large internal segment of the Rhodobacter FruK was found to be similar in sequence to the domain bearing the sugar bisphosphate-binding region of the large subunit of ribulose 1,5-bisphosphate carboxylase/oxygenase of plants and bacteria. Proteins of the PfkB family did not exhibit statistically significant sequence identity with PfkA of E. coli. PfkA, however, is homologous to other prokaryotic and eukaryotic ATP- and PPi-dependent Pfks (the PfkA family). These eukaryotic, ATP-dependent enzymes each consist of a homotetramer (mammalian) or a heterooctamer (yeasts), with each subunit containing an internal duplication of the size of the entire PfkA protein of E. coli. In some of these enzymes, additional domains are present. A phylogenetic tree was constructed for the PfkA family and revealed that the bacterial enzymes closely resemble the N-terminal domains of the eukaryotic enzyme subunits whereas the C-terminal domains have diverged more extensively. The PPi-dependent Pfk of potato is only distantly related to the ATP-dependent enzymes. On the basis of their similar functions, sizes, predicted secondary structures, and sequences, we suggest that the PfkA and PfkB families share a common evolutionary origin.
Collapse
Affiliation(s)
- L F Wu
- Department of Biology, University of California at San Diego, La Jolla 92093-0116
| | | | | | | | | | | |
Collapse
|
35
|
The sequence and reactive site of ecotin. A general inhibitor of pancreatic serine proteases from Escherichia coli. J Biol Chem 1991. [DOI: 10.1016/s0021-9258(18)38162-6] [Citation(s) in RCA: 46] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022] Open
|
36
|
Jiang XM, Neal B, Santiago F, Lee SJ, Romana LK, Reeves PR. Structure and sequence of the rfb (O antigen) gene cluster of Salmonella serovar typhimurium (strain LT2). Mol Microbiol 1991; 5:695-713. [PMID: 1710759 DOI: 10.1111/j.1365-2958.1991.tb00741.x] [Citation(s) in RCA: 279] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
The rfb gene cluster of Salmonella LT2 has been cloned and sequenced. The genes rfbA, rfbB, rfbD, rfbF, rfbG, rfbK, rfbM and rfbP were located individually and the gene rfbL was located outside the cluster. Approximately 16 open reading frames were found in the region which is essential for the expression of O antigen. The gene products of rfbB and rfbG were found to have homology with the group of dehydrogenase and related enzymes described previously. Analysis of the G + C ratio of the rfb cluster extended the area of low-G + C composition previously found in the sequence of rfbJ to the whole rfb gene cluster. Three to five segments with discrete G + C contents and codon adaptation indices are present in the rfb region, indicating a heterogeneous origin of these segments. Potential promoters were found near the start of the rfb region, supporting the possibility that the rfb gene cluster is an operon.
Collapse
Affiliation(s)
- X M Jiang
- Department of Microbiology, University of Sydney, New South Wales, Australia
| | | | | | | | | | | |
Collapse
|
37
|
Rowley DL, Wolf RE. Molecular characterization of the Escherichia coli K-12 zwf gene encoding glucose 6-phosphate dehydrogenase. J Bacteriol 1991; 173:968-77. [PMID: 1704005 PMCID: PMC207213 DOI: 10.1128/jb.173.3.968-977.1991] [Citation(s) in RCA: 75] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
In Escherichia coli K-12, expression of zwf, the gene for glucose 6-phosphate dehydrogenase, is coordinated with the cellular growth rate and induced by superoxide-generating agents. To initiate the study of the molecular mechanisms regulating its expression, the gene was cloned and its DNA sequence was determined. The 5' ends of zwf mRNA isolated from cells growing in glucose and acetate minimal media were mapped. The map was complex in that transcripts mapped to -45, -52, and -62, with respect to the beginning of the coding sequence. Three analytical methods were used to search the DNA sequence for putative promoters. Only one sequence for a promoter recognized by the sigma 70 form of RNA polymerase was found by all three search routines that could be aligned with a mapped transcript, indicating that the other transcripts arise by processing of the mRNA. A computer-assisted search did not reveal a thermodynamically stable long-range mRNA secondary structure that is capable of sequestering the translation initiation region, which suggests that growth-rate-dependent regulation of glucose 6-phosphate dehydrogenase level may not be carried out by a mechanism similar to the one for the gene (gnd) for 6-phosphogluconate dehydrogenase. The DNA segment between the -10 hexamer and the start point of transcription resembles the discriminator sequence of stable RNA genes, which has been implicated in stringent control and growth-rate-dependent regulation.
Collapse
Affiliation(s)
- D L Rowley
- Department of Biological Sciences, University of Maryland Baltimore County, Catonsville 21228
| | | |
Collapse
|
38
|
O'Neill MC. Training back-propagation neural networks to define and detect DNA-binding sites. Nucleic Acids Res 1991; 19:313-8. [PMID: 2014171 PMCID: PMC333596 DOI: 10.1093/nar/19.2.313] [Citation(s) in RCA: 63] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
A three layered back-propagation neural network was trained to recognize E. coli promoters of the 17 base spacing class. To this end, the network was presented with 39 promoter sequences and derivatives of those sequences as positive inputs; 60% A + T random sequences and sequences containing 2 promoter-down point mutations were used as negative inputs. The entire promoter sequence of 58 bases, approximately -50 to +8, was entered as input. The network was asked to associate an output of 1.0 with promoter sequence input and 0.0 with non-promoter input. Generally, after 100,000 input cycles, the network was virtually perfect in classifying the training set. A trained network was about 80% effective in recognizing 'new' promoters which were not in the training set, with a false positive rate below 0.1%. Network searches on pBR322 and on the lambda genome were also performed. Overall the results were somewhat better than the best rule-based procedures. The trained network can be analyzed both for its choice of base and relative weighting, positive and negative, at each position of the sequence. This method, which requires only appropriate input/output training pairs, can be used to define and search for any DNA regulatory sequence for which there are sufficient exemplars.
Collapse
Affiliation(s)
- M C O'Neill
- Department of Biological Sciences, University of Maryland, Baltimore County 21228
| |
Collapse
|
39
|
Nègre D, Cortay JC, Old IG, Galinier A, Richaud C, Saint Girons I, Cozzone AJ. Overproduction and characterization of the iclR gene product of Escherichia coli K-12 and comparison with that of Salmonella typhimurium LT2. Gene 1991; 97:29-37. [PMID: 1995431 DOI: 10.1016/0378-1119(91)90006-w] [Citation(s) in RCA: 28] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
The iclR gene of Escherichia coli K-12, which encodes a regulatory protein (repressor) for the aceBAK operon, is located between that operon and metH in the 91-min region of the chromosome. The iclR gene was cloned and expressed in a coupled T7 RNA polymerase/promoter system and the gene product was identified by specific binding to a fragment containing the aceBAK operator region. The iclR gene product is a polypeptide of 274 amino acids (aa) with a calculated Mr of 29,741. Comparison of the deduced IclR aa sequence to that of Salmonella typhimurium revealed that the two IclR repressors exhibit 89% identity. A possible helix-turn-helix motif characteristic of DNA-binding proteins was found within the IclR sequence. A search in protein data banks revealed that IclR has a score of similarity of 43.7% with GylR, a transcriptional regulator of the glycerol operon of Streptomyces coelicolor.
Collapse
Affiliation(s)
- D Nègre
- Institut de Biologie et Chimie des Protéines, C.N.R.S., Université de Lyon, France
| | | | | | | | | | | | | |
Collapse
|
40
|
Lawrence CE, Reilly AA. An expectation maximization (EM) algorithm for the identification and characterization of common sites in unaligned biopolymer sequences. Proteins 1990; 7:41-51. [PMID: 2184437 DOI: 10.1002/prot.340070105] [Citation(s) in RCA: 319] [Impact Index Per Article: 9.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Statistical methodology for the identification and characterization of protein binding sites in a set of unaligned DNA fragments is presented. Each sequence must contain at least one common site. No alignment of the sites is required. Instead, the uncertainty in the location of the sites is handled by employing the missing information principle to develop an "expectation maximization" (EM) algorithm. This approach allows for the simultaneous identification of the sites and characterization of the binding motifs. The reliability of the algorithm increases with the number of fragments, but the computations increase only linearly. The method is illustrated with an example, using known cyclic adenosine monophosphate receptor protein (CRP) binding sites. The final motif is utilized in a search for undiscovered CRP binding sites.
Collapse
Affiliation(s)
- C E Lawrence
- Biometrics Laboratory, Wadsworth Center for Laboratories and Research, New York State Department of Health, Albany 12201
| | | |
Collapse
|