Characteristic features of thermal stability map of DNA in Escherichia coli and eukaryotic genes

Distribution of double-helix thermal stability of Escherichia coli and eukaryotic DNAs was analyzed. The results confirmed the previous propositions based on the study of the stability distribution in phage DNAs: (1) stability fluctuation appears near the boundaries of protein coding regions (PCRs) and non protein coding regions (NPCRs); (2) PCRs have less fluctuation than NPCRs. The present analysis also revealed that the local G + C content is lower in the beginning of PCRs of E. coli than the average G + C content of PCR and that deviations in the amino acid composition and the third letter usage PCRs are involved in the low G + C content; the biological meaning of this is discussed in relation to mRNA structure.     

Evidence for effect of random genetic drift on G+C content after lateral transfer of fucose pathway genes to Escherichia coli K-12

The cps cluster of Escherichia coli K-12 comprises genes involved in synthesis of capsular polysaccharide colanic acid. Part of the E. coli K-12 cps region has been cloned and sequenced and compared to its Salmonella enterica LT2 counterpart. The cps genes from the two organisms are homologous; in the case of the LT2 genes, with G+C content of 0.61 and codons characteristic of high G+C species, it seems clear that they have been acquired relatively recently by lateral transfer from a high G+C species. The K-12 form of these cps genes is closely related to those of LT2 so must derive from the same high G+C species, but it appears to have transferred much earlier such that random genetic drift has brought P3 (the corrected G+C content of codon base 3) down from 0.77 to 0.64, more than halfway to the E. coli average of 0.57. We estimate, using an equation developed by Sueoka, that the lateral transfer to E. coli took place approximately 45 million years ago. This is the first report we are aware of demonstrating the expected adjustment of P3 after lateral transfer between species with different G+C content DNA.      

Escherichia coli genome is composed of two distinct types of nucleotide sequences

We calculated correlations of the nucleotide distributions along the E. coli genome. Subsequent cluster analysis of the correlation distributions showed that the genome was composed of two qualitatively different types of nucleotide sequences. The first type exhibited strong correlations of the genomic distributions of A with T and G with C, and high anticorrelations of A with C and G with T. In contrast, the second type was characterized by weak or negligible correlations typical of randomized sequences. Both types of sequences were almost equally abundant in the E. coli genome and their length varied from several hundred nucleotides to about 70 kilobases. They were not disjunct with respect to their (G + C) content but the high correlations and anticorrelations were rather characteristic for (A + T)-rich genomic segments. We offer possible explanations of the mosaic structure of the E. coli genome.     

Identification and characterization of E.coli ribosomal binding sites by free energy computation.

Sequences upstream from translational initiation sites of different E.coli genes show various degrees of complementarity to the Shine-Dalgarno (SD) sequence at the 3' end of the 16S rRNA. We propose a quantitative measure for the SD region on the mRNA, that reflects its degree of complementarity to the rRNA. This measure is based on the stability of the rRNA-mRNA duplex as established by free energy computations. The free energy calculations are based on the same principles that are used for folding a single RNA molecule, and are executed by similar algorithms. Bulges and internal loops in the rRNA and mRNA are allowed. The mRNA string with maximum free energy gain upon binding to the rRNA is selected as the most favorable SD sequence of a gene. The free energy value that represents the SD region provides a quantitative measure that can be used for comparing SD sequences of different genes. The distribution of this measure in more than 1000 E.coli genes is presented and discussed.     

mRNAs have greater negative folding free energies than shuffled or codon choice randomized sequences

An examination of 51 mRNA sequences in GenBank has revealed that calculated mRNA folding is more stable than expected by chance. Free energy minimization calculations of native mRNA sequences are more negative than randomized mRNA sequences with the same base composition and length. Randomization of the coding region of genes yields folding free energies of less negative magnitude than the original native mRNA sequence. Randomization of codon choice, while still preserving original base composition, also results in less stable mRNAs. This suggests that a bias in the selection of codons favors the potential formation of mRNA structures which contribute to folding stability.     

Organization and codon usage of the streptomycin operon in Micrococcus luteus, a bacterium with a high genomic G + C content.

The DNA sequence of the Micrococcus luteus str operon, which includes genes for ribosomal proteins S12 (str or rpsL) and S7 (rpsG) and elongation factors (EF) G (fus) and Tu (tuf), has been determined and compared with the corresponding sequence of Escherichia coli to estimate the effect of high genomic G + C content (74%) of M. luteus on the codon usage pattern. The gene organization in this operon and the deduced amino acid sequence of each corresponding protein are well conserved between the two species. The mean G + C content of the M. luteus str operon is 67%, which is much higher than that of E. coli (51%). The codon usage pattern of M. luteus is very different from that of E. coli and extremely biased to the use of G and C in silent positions. About 95% (1,309 of 1,382) of codons have G or C at the third position. Codon GUG is used for initiation of S12, EF-G, and EF-Tu, and AUG is used only in S7, whereas GUG initiates only one of the EF-Tu's in E. coli. UGA is the predominant termination codon in M. luteus, in contrast to UAA in E. coli.     

Cloning and expression in Escherichia coli of two additional amylase genes of a strictly anaerobic thermophile, Dictyoglomus thermophilum, and their nucleotide sequences with extremely low guanine-plus-cytosine contents

An obligately anaerobic and extremely thermophilic bacterium, Dictyoglomus thermophilum, produces multiple extracellular amylases. In addition to one of the amylase genes, amyA, which we previously cloned and characterized, we have cloned two additional genes, amyB and amyC, coding for amylases of this thermophile, into Escherichia coli and determined their nucleotide sequences. The two amylase genes were expressed under the control of E. coli promoters. Almost all activity was detected in the intracellular fraction in the E. coli cells. The molecular mass and NH2-terminal amino acid sequence of the AmyB enzyme, which was purified from an E. coli transformant containing the amyB gene, confirmed that the reading frame of amyB consisted of 562 amino acids (Mr 67,000). The molecular mass of the AmyC enzyme, estimated by activity staining of a crude extract of E. coli containing amyC, confirmed that AmyC consisted of 498 amino acids (Mr 59,000). The optimal temperatures for AmyB and AmyC activities on soluble starch were 80 degrees C and 70 degrees C, respectively. Both AmyB and AmyC showed a pH optimum of 5.5. AmyB and AmyC showed a different pattern of starch hydrolysis when examined by thin-layer chromatography. Some homology in the amino acid sequences with the functional regions of Taka-amylase A was found in both AmyB and AmyC. The codon usage in the amyA, amyB and amyC genes was highly biased, which reflects the fact that the guanine-plus-cytosine (G + C) content of DNA of D. thermophilum is 29 mol%. The distribution of G and C at each position of the codons was non-random; the G + C content of the first position of codons is significantly high, whereas that of the third position is somewhat low. In addition, codons consisting only of A and T were preferentially used in this thermophile.     

The relation between mRNA folding and protein structure

About 200 mRNA sequences of Escherichia coli and human with matching protein secondary structure data were studied. The mRNA folding for each native sequence and for corresponding randomized sequences was calculated through free energy minimization. We have found that the folding energy of mRNA segments in different protein secondary structures is significantly different. The average Z score is more negative for regular secondary structure (alpha-helix and beta-strand) than that for coil. This suggests that the codon choice in native mRNA sequence coding for protein regular structure contributes more to the mRNA folding stability.