Cloning and characterization of the Zymobacter palmae pyruvate decarboxylase gene (pdc) and comparison to bacterial homologues

Pyruvate decarboxylase (PDC) is the key enzyme in all homo-ethanol fermentations. Although widely distributed among plants, yeasts, and fungi, PDC is absent in animals and rare in bacteria (established for only three organisms). Genes encoding the three known bacterial pdc genes have been previously described and expressed as active recombinant proteins. The pdc gene from Zymomonas mobilis has been used to engineer ethanol-producing biocatalysts for use in industry. In this paper, we describe a new bacterial pdc gene from Zymobacter palmae. The pattern of codon usage for this gene appears quite similar to that for Escherichia coli genes. In E. coli recombinants, the Z. palmae PDC represented approximately 1/3 of the soluble protein. Biochemical and kinetic properties of the Z. palmae enzyme were compared to purified PDCs from three other bacteria. Of the four bacterial PDCs, the Z. palmae enzyme exhibited the highest specific activity (130 U mg of protein(-1)) and the lowest Km for pyruvate (0.24 mM). Differences in biochemical properties, thermal stability, and codon usage may offer unique advantages for the development of new biocatalysts for fuel ethanol production.     

Expression and surface display of Cellulomonas endoglucanase in the ethanologenic bacterium Zymobacter palmae

In order to reduce the cost of bioethanol production from lignocellulosic biomass, we developed a tool for cell surface display of cellulolytic enzymes on the ethanologenic bacterium Zymobacter palmae. Z. palmae is a novel ethanol-fermenting bacterium capable of utilizing a broad range of sugar substrates, but not cellulose. Therefore, to express and display heterologous cellulolytic enzymes on the Z. palmae cell surface, we utilized the cell-surface display motif of the Pseudomonas ice nucleation protein Ina. The gene encoding Ina from Pseudomonas syringae IFO3310 was cloned, and its product was comprised of three functional domains: an N-terminal domain, a central domain with repeated amino acid residues, and a C-terminal domain. The N-terminal domain of Ina was shown to function as the anchoring motif for a green fluorescence protein fusion protein in Escherichia coli. To express a heterologous cellulolytic enzyme extracellularly in Z. palmae, we fused the N-terminal coding sequence of Ina to the coding sequence of an N-terminal-truncated Cellulomonas endoglucanase. Z. palmae cells carrying the fusion endoglucanase gene were shown to degrade carboxymethyl cellulose. Although a portion of the expressed fusion endoglucanase was released from Z. palmae cells into the culture broth, we confirmed the display of the protein on the cell surface by immunofluorescence microscopy. The results indicate that the N-terminal anchoring motif of Ina from P. syringae enabled the translocation and display of the heterologous cellulase on the cell surface of Z. palmae.     

Transfer RNA Gene Numbers may not be Completely Responsible for the Codon Usage Bias in Asparagine, Isoleucine, Phenylalanine, and Tyrosine in the High Expression Genes in Bacteria

It is generally believed that the effect of translational selection on codon usage bias is related to the number of transfer RNA genes in bacteria, which is more with respect to the high expression genes than the whole genome. Keeping this in the background, we analyzed codon usage bias with respect to asparagine, isoleucine, phenylalanine, and tyrosine amino acids. Analysis was done in seventeen bacteria with the available gene expression data and information about the tRNA gene number. In most of the bacteria, it was observed that codon usage bias and tRNA gene number were not in agreement, which was unexpected. We extended the study further to 199 bacteria, limiting to the codon usage bias in the two highly expressed genes rpoB and rpoC which encode the RNA polymerase subunits β and β′, respectively. In concordance with the result in the high expression genes, codon usage bias in rpoB and rpoC genes was also found to not be in agreement with tRNA gene number in many of these bacteria. Our study indicates that tRNA gene numbers may not be the sole determining factor for translational selection of codon usage bias in bacterial genomes.     

Genetic Engineering of Zymobacter palmae for Production of Ethanol from Xylose

Its metabolic characteristics suggest that Zymobacter palmae gen. nov., sp. nov. could serve as a useful new ethanol-fermenting bacterium, but its biotechnological exploitation will require certain genetic modifications. We therefore engineered Z. palmae so as to broaden the range of its fermentable sugar substrates to include the pentose sugar xylose. The Escherichia coli genes encoding the xylose catabolic enzymes xylose isomerase, xylulokinase, transaldolase, and transketolase were introduced into Z. palmae, where their expression was driven by the Zymomonas mobilis glyceraldehyde-3-phosphate dehydrogenase promoter. When cultured with 40 g/liter xylose, the recombinant Z. palmae strain was able to ferment 16.4 g/liter xylose within 5 days, producing 91% of the theoretical yield of ethanol with no accumulation of organic acids as metabolic by-products. Notably, xylose acclimation enhanced both the expression of xylose catabolic enzymes and the rate of xylose uptake into recombinant Z. palmae, which enabled the acclimated organism to completely and simultaneously ferment a mixture of 40 g/liter glucose and 40 g/liter xylose within 8 h, producing 95% of the theoretical yield of ethanol. Thus, efficient fermentation of a mixture of glucose and xylose to ethanol can be accomplished by using Z. palmae expressing E. coli xylose catabolic enzymes.     

Multiscale Modeling of Metabolism and Macromolecular Synthesis in E. coli and Its Application to the Evolution of Codon Usage

Biological systems are inherently hierarchal and multiscale in time and space. A major challenge of systems biology is to describe biological systems as a computational model, which can be used to derive novel hypothesis and drive experiments leading to new knowledge. The constraint-based reconstruction and analysis approach has been successfully applied to metabolism and to the macromolecular synthesis machinery assembly. Here, we present the first integrated stoichiometric multiscale model of metabolism and macromolecular synthesis for Escherichia coli K12 MG1655, which describes the sequence-specific synthesis and function of almost 2000 gene products at molecular detail. We added linear constraints, which couple enzyme synthesis and catalysis reactions. Comparison with experimental data showed improvement of growth phenotype prediction with the multiscale model over E. coli’s metabolic model alone. Many of the genes covered by this integrated model are well conserved across enterobacters and other, less related bacteria. We addressed the question of whether the bias in synonymous codon usage could affect the growth phenotype and environmental niches that an organism can occupy. We created two classes of in silico strains, one with more biased codon usage and one with more equilibrated codon usage than the wildtype. The reduced growth phenotype in biased strains was caused by tRNA supply shortage, indicating that expansion of tRNA gene content or tRNA codon recognition allow E. coli to respond to changes in codon usage bias. Our analysis suggests that in order to maximize growth and to adapt to new environmental niches, codon usage and tRNA content must co-evolve. These results provide further evidence for the mutation-selection-drift balance theory of codon usage bias. This integrated multiscale reconstruction successfully demonstrates that the constraint-based modeling approach is well suited to whole-cell modeling endeavors.     

Codon usage in highly expressed genes of Haemophillus influenzae and Mycobacterium tuberculosis: translational selection versus mutational bias

Biases in the codon usage and base compositions at three codon sites in different genes of A+T-rich Gram-negative bacterium Haemophillus influenzae and G+C-rich Gram-positive bacterium Mycobacterium tuberculosis have been examined to address the following questions: (1) whether the synonymous codon usage in organisms having highly skewed base compositions is totally dictated by the mutational bias as reported previously (Sharp, P.M., Devine, K.M., 1989. Codon usage and gene expression level in Dictyostelium discoideum: highly expressed genes do `prefer' optimal codons. Nucleic Acids Res. 17, 5029–5039), or is also controlled by translational selection; (2) whether preference of G in the first codon positions by highly expressed genes, as reported in Escherichia coli (Gutierrez, G., Marquez, L., Marin, A., 1996. Preference for guanosine at first codon position in highly expressed Escherichia coli genes. A relationship with translational efficiency. Nucleic Acids Res. 24, 2525–2527), is true in other bacteria; and (3) whether the usage of bases in three codon positions is species-specific. Result presented here show that even in organisms with high mutational bias, translational selection plays an important role in dictating the synonymous codon usage, though the set of optimal codons is chosen in accordance with the mutational pressure. The frequencies of G-starting codons are positively correlated to the level of expression of genes, as estimated by their Codon Adaptation Index (CAI) values, in M. tuberculosis as well as in H. influenzae in spite of having an A+T-rich genome. The present study on the codon preferences of two organisms with oppositely skewed base compositions thus suggests that the preference of G-starting codons by highly expressed genes might be a general feature of bacteria, irrespective of their overall G+C contents. The ranges of variations in the frequencies of individual bases at the first and second codon positions of genes of both H. influenzae and M. tuberculosis are similar to those of E. coli, implying that though the composition of all three codon positions is governed by a selection-mutation balance, the mutational pressure has little influence in the choice of bases at the first two codon positions, even in organisms with highly biased base compositions.     

Cloning and expression of the structural gene for pyruvate decarboxylase of Zymomonas mobilis in Escherichia coli

A genomic library of Zymomonas mobilis DNA was constructed in Escherichia coli using cosmid vector pHC79. Immunological screening of 483 individual E. coli strains revealed two clones expressing pyruvate decarboxylase, the key enzyme for efficient ethanol production of Z. mobilis. The two plasmids, pZM1 and pZM2, isolated from both E. coli strains were found to be related and to exhibit a common 4.6 kb SphI fragment on which the gene coding for pyruvate decarboxylase, pdc, was located.The pdc gene was similarily well expressed in both aerobically and anaerobically grown E. coli cells, and exerted a considerable effect on the amount of fermentation products formed. During fermentative growth on 25 mM glucose, plasmid-free E. coli lacking a pdc gene produced 6.5 mM ethanol, 8.2 mM acetate, 6.5 mM lactate, 0.5 mM succinate, and about 1 mM formate leaving 10.4 mM residual glucose. In contrast, recombinant E. coli harbouring a cloned pdc gene from Z. mobilis completely converted 25 mM glucose to up to 41.5 mM ethanol while almost no acids were formed.     

Evolutionary implications of potential frameshift hot-spots, methylation sites, and codon usage differences between genes and between organisms in Escherichia coli and Salmonella typhimurium

The published nucleotide sequences of the E. coli and S. typhimurium trp A and trp B genes show a high degree of similarity between homologous genes of the two organisms, and an even greater degree of similarity between the amino acid sequences of the gene products. In spite of this, analysis of the nucleotide sequences reveals that there are marked differences between E. coli and S. typhimurium genes with respect to potential frameshift mutation hot-spots and dam and mec, mutationally important, methylation sites. Such existing differences may well lead to divergent evolution of these two, presently closely related, bacteria. Codon usage patterns in the trp A and trp B genes of E. coli and S. typhimurium, and the lac I gene of E. coli, have been re-analysed in terms of AT-rich, GC-rich, neutral, or unique codons and marked preferences found. In some cases particular amino acids are most often specified by AT-rich, in others by the GC-rich, alternative codons. In still other cases the codon preference depends on the gene studied. These patterns can be interpreted in terms of enteric bacterial evolution, via hybridizations, from ancestral bacteria with AT- or GC-rich DNA.     

The Pyruvate Dehydrogenase Complexes: Structure-based Function and Regulation

The pyruvate dehydrogenase complexes (PDCs) from all known living organisms comprise three principal catalytic components for their mission: E1 and E2 generate acetyl-coenzyme A, whereas the FAD/NAD⁺-dependent E3 performs redox recycling. Here we compare bacterial (Escherichia coli) and human PDCs, as they represent the two major classes of the superfamily of 2-oxo acid dehydrogenase complexes with different assembly of, and interactions among components. The human PDC is subject to inactivation at E1 by serine phosphorylation by four kinases, an inactivation reversed by the action of two phosphatases. Progress in our understanding of these complexes important in metabolism is reviewed.     

Characterization of synonymous codon usage bias in the pseudorabies virus US1 gene

In the present study, we examined the codon usage bias between pseudorabies virus (PRV) US1 gene and the US1-like genes of 20 reference alphaherpesviruses. Comparative analysis showed noticeable disparities of the synonymous codon usage bias in the 21 alphaherpesviruses, indicated by codon adaptation index, effective number of codons (ENc) and GC3s value. The codon usage pattern of PRV US1 gene was phylogenetically conserved and similar to that of the US1-like genes of the genus Varicellovirus of alphaherpesvirus, with a strong bias towards the codons with C and G at the third codon position. Cluster analysis of codon usage pattern of PRV US1 gene with its reference alphaherpesviruses demonstrated that the codon usage bias of US1-like genes of 21 alphaherpesviruses had a very close relation with their gene functions. ENc-plot revealed that the genetic heterogeneity in PRV US1 gene and the 20 reference alphaherpesviruses was constrained by G+C content, as well as the gene length. In addition, comparison of codon preferences in the US1 gene of PRV with those of E. coli, yeast and human revealed that there were 50 codons showing distinct usage differences between PRV and yeast, 49 between PRV and human, but 48 between PRV and E. coli. Although there were slightly fewer differences in codon usages between E.coli and PRV, the difference is unlikely to be statistically significant, and experimental studies are necessary to establish the most suitable expression system for PRV US1. In conclusion, these results may improve our understanding of the evolution, pathogenesis and functional studies of PRV, as well as contributing to the area of herpesvirus research or even studies with other viruses.     

Knockout of the p-Coumarate Decarboxylase Gene from Lactobacillus plantarum Reveals the Existence of Two Other Inducible Enzymatic Activities Involved in Phenolic Acid Metabolism

Lactobacillus plantarum NC8 contains a pdc gene coding for p-coumaric acid decarboxylase activity (PDC). A food grade mutant, designated LPD1, in which the chromosomal pdc gene was replaced with the deleted pdc gene copy, was obtained by a two-step homologous recombination process using an unstable replicative vector. The LPD1 mutant strain remained able to weakly metabolize p-coumaric and ferulic acids into vinyl derivatives or into substituted phenyl propionic acids. We have shown that L. plantarum has a second acid phenol decarboxylase enzyme, better induced with ferulic acid than with p-coumaric acid, which also displays inducible acid phenol reductase activity that is mostly active when glucose is added. Those two enzymatic activities are in competition for p-coumaric and ferulic acid degradation, and the ratio of the corresponding derivatives depends on induction conditions. Moreover, PDC appeared to decarboxylate ferulic acid in vitro with a specific activity of about 10 nmol · min⁻¹ · mg⁻¹ in the presence of ammonium sulfate. Finally, PDC activity was shown to confer a selective advantage on LPNC8 grown in acidic media supplemented with p-coumaric acid, compared to the LPD1 mutant devoid of PDC activity.     

Screening of yeasts for cell-free production of (R)-phenylacetylcarbinol

105 yeast strains from 10 genera and 40 species were evaluated for cell-free production of (R)-phenylacetylcarbinol (PAC), the chiral precursor in the manufacture of the pharmaceuticals ephedrine and pseudoephedrine. Carboligase activity of pyruvate decarboxylase (PDC), forming PAC from benzaldehyde and pyruvate, was found in extracts of 98 strains. PAC was not formed from benzaldehyde and acetaldehyde, an activity of bacterial PDCs from Zymomonas mobilis and Zymobacter palmae. Two interesting groups of candidates were identified in the yeast screening: carboligase activities of Schizosaccharomyces pombe PDCs were very low but showed best resistance to pre-incubation with acetaldehyde and benzaldehyde; and highest carboligase activities combined with medium resistance were found in strains of Candida utilis, C. tropicalis and C. albicans.     

Evolutionary Lability of Context-Dependent Codon Bias in Bacteria

In bacteria, synonymous codon usage can be considerably affected by base composition at neighboring sites. Such context-dependent biases may be caused by either selection against specific nucleotide motifs or context-dependent mutation biases. Here we consider the evolutionary conservation of context-dependent codon bias across 11 completely sequenced bacterial genomes. In particular, we focus on two contextual biases previously identified in Escherichia coli; the avoidance of out-of-frame stop codons and AGG motifs. By identifying homologues of E. coli genes, we also investigate the effect of gene expression level in Haemophilus influenzae and Mycoplasma genitalium. We find that while context-dependent codon biases are widespread in bacteria, few are conserved across all species considered. Avoidance of out-of-frame stop codons does not apply to all stop codons or amino acids in E. coli, does not hold for different species, does not increase with gene expression level, and is not relaxed in Mycoplasma spp., in which the canonical stop codon, TGA, is recognized as tryptophan. Avoidance of AGG motifs shows some evolutionary conservation and increases with gene expression level in E. coli, suggestive of the action of selection, but the cause of the bias differs between species. These results demonstrate that strong context-dependent forces, both selective and mutational, operate on synonymous codon usage but that these differ considerably between genomes. Received: 6 May 1999 / Accepted: 29 October 1999     

Enhancing the recombinant protein expression of halohydrin dehalogenase HheA in Escherichia coli by applying a codon optimization strategy

Halohydrin dehalogenases are attractive biocatalysts in producing a series of important chiral building blocks. Recombinant expression of halohydrin dehalogenase from Arthrobacter sp. AD2 (HheA) in Escherichia coli using T7 promoter-based pGEF(+) system revealed much lower expression level than that of the well-studied halohydrin dehalogenase from Agrobacterium radiobacter AD1 (HheC). In this study, we changed the codon usage in the 5′-end of hheA gene to improve the expression yield of HheA. Our results showed that the expression of HheA could be largely improved by the replacement of G-rich +2 codon (adjacent to the start codon) with less G-containing codons. The expression of one of the resulting mutants HheA-D1 (replaced +2 codon GTG with CCA) was about 4-fold higher and purified yields about 8-fold greater than that of the wild-type HheA. Moreover, the expression level of the resulting HheA variants correlated well with the minimal folding free energy (ΔG) of the mRNA secondary structure surrounding the 5′-end region of the genes. These findings suggested that the G-rich +2 codon of hheA gene might be the main suppressive factor for limiting the recombinant expression of HheA and that +2 codon optimization strategy could be used as a general tool in modulating recombinant protein production in E. coli.     

Impact of translational selection on codon usage bias in the archaeon Methanococcus maripaludis

Patterns of codon usage have been extensively studied among Bacteria and Eukaryotes, but there has been little investigation of species from the third domain of life, the Archaea. Here, we examine the nature of codon usage bias in a methanogenic archaeon, Methanococcus maripaludis. Genome-wide patterns of codon usage are dominated by a strong A + T bias, presumably largely reflecting mutation patterns. Nevertheless, there is variation among genes in the use of a subset of putatively translationally optimal codons, which is strongly correlated with gene expression level. In comparison with Bacteria such as Escherichia coli, the strength of selected codon usage bias in highly expressed genes in M. maripaludis seems surprisingly high given its moderate growth rate. However, the pattern of selected codon usage differs between M. maripaludis and E. coli: in the archaeon, strongly selected codon usage bias is largely restricted to twofold degenerate amino acids (AAs). Weaker bias among the codons for fourfold degenerate AAs is consistent with the small number of tRNA genes in the M. maripaludis genome.     

Construction of expression vectors for the gram-negative bacterium Zymomonas mobilis

Summary A set of vectors was constructed for the cloning and expression of heterologous genes in the Gramnegative bacterium Zymomonas mobilis under the control of the pdc promoter of Z. mobilis. The vectors pPTZ1, pPTZ3, and pPTZ4 are based on the cryptic Z. mobilis plasmid pZM02 and on parts of the Escherichia coli plasmids pKK223-3 and pBR322 together with the multiple cloning site of phage Ml3mp18. DNA fragments can be readily inserted immediately downstream from the pdc promoter at unique restriction sites for KpnI, XbaI and PstI in pPTZl and additionally for SmaI and BamHI in pPTZ3. In pPTZ4, the 5 terminal codons of pdc were deleted allowing the formation of gene fusions. Expression of a promoterless chloramphenicol acetyltransferase gene (cat) controlled by the pdc gene promoter resulted in enzyme activities of up to 5.5 U/mg total cell protein in Z. mobilis cells.     

Statistical analysis of pentose phosphate pathway genes from eubacteria and eukarya reveals translational selection as a major force in shaping codon usage pattern

Comparative analysis of metabolic pathways among widely diverse species provides an excellent opportunity to extract information about the functional relation of organisms and pentose phosphate pathway exemplifies one such pathway. A comparative codon usage analysis of the pentose phosphate pathway genes of a diverse group of organisms representing different niches and the related factors affecting codon usage with special reference to the major forces influencing codon usage patterns was carried out. It was observed that organism specific codon usage bias percolates into vital metabolic pathway genes irrespective of their near universality. A clear distinction in the codon usage pattern of gram positive and gram negative bacteria, which is a major classification criterion for bacteria, in terms of pentose phosphate pathway was an important observation of this study. The codon utilization scheme in all the organisms indicates the presence of translation selection as a major force in shaping codon usage. Another key observation was the segregation of the H. sapiens genes as a separate cluster by correspondence analysis, which is primarily attributed to the different codon usage pattern in this genus along with its longer gene lengths. We have also analyzed the amino acid distribution comparison of transketolase protein primary structures among all the organisms and found that there is a certain degree of predictability in the composition profile except in A. fumigatus and H. sapiens, where few exceptions are prominent. In A. fumigatus, a human pathogen responsible for invasive aspergillosis, a significantly different codon usage pattern, which finally translated into its amino acid composition model portraying a unique profile in a key pentose phosphate pathway enzyme transketolase was observed.     

Study of codon bias perspective of fungal xylanase gene by multivariate analysis

Fungal xylanases has important applications in food, baking, pulp and paper industries in addition to various other industries. Xylanases are produced extensively by both bacterial and fungal sources and has tremendous potential of being active at extremes of temperature and pH. In the present study an effort has been made to explore the codon bias perspective of this potential enzyme using bioinformatics tools. Multivariate analysis has been used as a tool to study codon bias perspectives of xylanases. It was further observed that the codon usage of xylanases genes from different fungal sources is not similar and to reveal this phenomenon the relative synonymous codon usage (RSCU) and base composition variation in fungal xylanase genes were also studied. The codon biasing data like GC content at third position (GC_3S), effective codon number (N_C), codon adaptive index (CAI) were further analyzed with statistical softwares like Sigma1plot 9.0 and Systat 11.0. Furthermore, study of translation selection was also performed to verify the influences of codon usage variation among the 94 xylanase genes. In the present study xylanase gene from 12 organisms were analyzed and codon usages of all xylanases from each organism were compared separately. Analysis indicates biased codon among all 12 fungi taken for study with Aspergillus nidulans, Chaetomium globosum, Aspergillus terreus and Aspergillus clavatus showing maximum biasing. N_C plot and correspondence analysis on relative synonymous codon usage indicate that mutation bias and translation selection influences codon usage variation in fungal xylanase gene. To reveal the relative synonymous codon usage and base composition variation in xylanase, 94 genes from 12 fungi were used as model system.     

Direct ethanol production from cellulosic materials by Zymobacter palmae carrying Cellulomonas endoglucanase and Ruminococcus β-glucosidase genes

In order to reduce the cost of bioethanol production from lignocellulosic biomass, we conferred the ability to ferment cellulosic materials directly on Zymobacter palmae by co-expressing foreign endoglucanase and β-glucosidase genes. Z. palmae is a novel ethanol-fermenting bacterium capable of utilizing a broad range of sugar substrates, but not cellulose. Therefore, the six genes encoding the cellulolytic enzymes (CenA, CenB, CenD, CbhA, CbhB, and Cex) from Cellulomonas fimi were introduced and expressed in Z. palmae. Of these cellulolytic enzyme genes cloned, CenA degraded carboxymethylcellulose and phosphoric acid-swollen cellulose (PASC) efficiently. The extracellular CenA catalyzed the hydrolysis of barley β-glucan and PASC to liberate soluble cello-oligosaccharides, indicating that CenA is the most suitable enzyme for cellulose degradation among those cellulolytic enzymes expressed in Z. palmae. Furthermore, the cenA gene and β-glucosidase gene (bgl) from Ruminococcus albus were co-expressed in Z. palmae. Of the total endoglucanase and β-glucosidase activities, 57.1 and 18.1 % were localized in the culture medium of the strain. The genetically engineered strain completely saccharified and fermented 20 g/l barley β-glucan to ethanol within 84 h, producing 79.5 % of the theoretical yield. Thus, the production and secretion of CenA and BGL enabled Z. palmae to efficiently ferment a water-soluble cellulosic polysaccharide to ethanol.