Analysis and manipulation of aspartate pathway genes for L-lysine overproduction from methanol by Bacillus methanolicus

Giving a signal: how protein phosphorylation helps Bacillus navigate through different life stages

Protein phosphorylation is a universal mechanism regulating a wide range of cellular responses across all domains of life. The antagonistic activities of kinases and phosphatases can orchestrate the life cycle of an organism. The availability of bacterial genome sequences, particularly Bacillus species, followed by proteomics and functional studies have aided in the identification of putative protein kinases and protein phosphatases, and their downstream substrates. Several studies have established the role of phosphorylation in different physiological states of Bacillus species as they pass through various life stages such as sporulation, germination, and biofilm formation. The most common phosphorylation sites in Bacillus proteins are histidine, aspartate, tyrosine, serine, threonine, and arginine residues. Protein phosphorylation can alter protein activity, structural conformation, and protein-protein interactions, ultimately affecting the downstream pathways. In this review, we summarize the knowledge available in the field of Bacillus signaling, with a focus on the role of protein phosphorylation in its physiological processes.

Efficient cell factories for the production of N-methylated amino acids and for methanol-based amino acid production

The growing world needs commodity amino acids such as L-glutamate and L-lysine for use as food and feed, and specialty amino acids for dedicated applications. To meet the supply a paradigm shift regarding their production is required. On the one hand, the use of sustainable and cheap raw materials is necessary to sustain low production cost and decrease detrimental effects of sugar-based feedstock on soil health and food security caused by competing uses of crops in the feed and food industries. On the other hand, the biotechnological methods to produce functionalized amino acids need to be developed further, and titres enhanced to become competitive with chemical synthesis methods. In the current review, we present successful strain mutagenesis and rational metabolic engineering examples leading to the construction of recombinant bacterial strains for the production of amino acids such as L-glutamate, L-lysine, L-threonine and their derivatives from methanol as sole carbon source. In addition, the fermentative routes for bioproduction of N-methylated amino acids are highlighted, with focus on three strategies: partial transfer of methylamine catabolism, S-adenosyl-L-methionine dependent alkylation and reductive methylamination of 2-oxoacids.

Discovery of Early-Branching Wolbachia Reveals Functional Enrichment on Horizontally Transferred Genes

Wolbachia is a widespread endosymbiont of insects and filarial nematodes that profoundly influences host biology. Wolbachia has also been reported in rhizosphere hosts, where its diversity and function remain poorly characterized. The discovery that plant-parasitic nematodes (PPNs) host Wolbachia strains with unknown roles is of interest evolutionarily, ecologically, and for agriculture as a potential target for developing new biological controls. The goal of this study was to screen communities for PPN endosymbionts and analyze genes and genomic patterns that might indicate their role. Genome assemblies revealed 1 out of 16 sampled sites had nematode communities hosting a Wolbachia strain, designated wTex, that has highly diverged as one of the early supergroup L strains. Genome features, gene repertoires, and absence of known genes for cytoplasmic incompatibility, riboflavin, biotin, and other biosynthetic functions placed wTex between mutualist C + D strains and reproductive parasite A + B strains. Functional terms enriched in group L included protoporphyrinogen IX, thiamine, lysine, fatty acid, and cellular amino acid biosynthesis, while dN/dS analysis suggested the strongest purifying selection on arginine and lysine metabolism, and vitamin B6, heme, and zinc ion binding, suggesting these as candidate roles in PPN Wolbachia. Higher dN/dS pathways between group L, wPni from aphids, wFol from springtails, and wCfeT from cat fleas suggested distinct functional changes characterizing these early Wolbachia host transitions. PPN Wolbachia had several putative horizontally transferred genes, including a lysine biosynthesis operon like that of the mitochondrial symbiont Midichloria, a spirochete-like thiamine synthesis operon shared only with wCfeT, an ATP/ADP carrier important in Rickettsia, and a eukaryote-like gene that may mediate plant systemic acquired resistance through the lysine-to-pipecolic acid system. The Discovery of group L-like variants from global rhizosphere databases suggests diverse PPN Wolbachia strains remain to be discovered. These findings support the hypothesis of plant-specialization as key to shaping early Wolbachia evolution and present new functional hypotheses, demonstrating promise for future genomics-based rhizosphere screens.

Aerobic Utilization of Methanol for Microbial Growth and Production

Methanol is a reduced one-carbon (C1) compound. It supports growth of aerobic methylotrophs that gain ATP from reduced redox equivalents by respiratory phosphorylation in their electron transport chains. Notably, linear oxidation of methanol to carbon dioxide may yield three reduced redox equivalents if methanol oxidation is NAD-dependent as, e.g., in Bacillus methanolicus. Methanol has a higher degree of reduction per carbon than glucose (6 vs. 4), and thus, lends itself as an ideal carbon source for microbial production of reduced target compounds. However, C-C bond formation in the RuMP or serine cycle, a prerequisite for production of larger molecules, requires ATP and/or reduced redox equivalents. Moreover, heat dissipation and a high demand for oxygen during catabolic oxidation of methanol may pose challenges for fermentation processes. In this chapter, we summarize metabolic pathways for aerobic methanol utilization, aerobic methylotrophs as industrial production hosts, strain engineering, and methanol bioreactor processes. In addition, we provide technological and market outlooks.

Evaluation of Heterologous Biosynthetic Pathways for Methanol-Based 5-Aminovalerate Production by Thermophilic Bacillus methanolicus

The use of methanol as carbon source for biotechnological processes has recently attracted great interest due to its relatively low price, high abundance, high purity, and the fact that it is a non-food raw material. In this study, methanol-based production of 5-aminovalerate (5AVA) was established using recombinant Bacillus methanolicus strains. 5AVA is a building block of polyamides and a candidate to become the C5 platform chemical for the production of, among others, δ-valerolactam, 5-hydroxy-valerate, glutarate, and 1,5-pentanediol. In this study, we test five different 5AVA biosynthesis pathways, whereof two directly convert L-lysine to 5AVA and three use cadaverine as an intermediate. The conversion of L-lysine to 5AVA employs lysine 2-monooxygenase (DavB) and 5-aminovaleramidase (DavA), encoded by the well-known Pseudomonas putida cluster davBA, among others, or lysine α-oxidase (RaiP) in the presence of hydrogen peroxide. Cadaverine is converted either to γ-glutamine-cadaverine by glutamine synthetase (SpuI) or to 5-aminopentanal through activity of putrescine oxidase (Puo) or putrescine transaminase (PatA). Our efforts resulted in proof-of-concept 5AVA production from methanol at 50°C, enabled by two pathways out of the five tested with the highest titer of 0.02 g l^-1. To our knowledge, this is the first report of 5AVA production from methanol in methylotrophic bacteria, and the recombinant strains and knowledge generated should represent a valuable basis for further improved 5AVA production from methanol.

Developing a Riboswitch-Mediated Regulatory System for Metabolic Flux Control in Thermophilic Bacillus methanolicus

Genome-wide transcriptomic data obtained in RNA-seq experiments can serve as a reliable source for identification of novel regulatory elements such as riboswitches and promoters. Riboswitches are parts of the 5′ untranslated region of mRNA molecules that can specifically bind various metabolites and control gene expression. For that reason, they have become an attractive tool for engineering biological systems, especially for the regulation of metabolic fluxes in industrial microorganisms. Promoters in the genomes of prokaryotes are located upstream of transcription start sites and their sequences are easily identifiable based on the primary transcriptome data. Bacillus methanolicus MGA3 is a candidate for use as an industrial workhorse in methanol-based bioprocesses and its metabolism has been studied in systems biology approaches in recent years, including transcriptome characterization through RNA-seq. Here, we identify a putative lysine riboswitch in B. methanolicus, and test and characterize it. We also select and experimentally verify 10 putative B. methanolicus-derived promoters differing in their predicted strength and present their functionality in combination with the lysine riboswitch. We further explore the potential of a B. subtilis-derived purine riboswitch for regulation of gene expression in the thermophilic B. methanolicus, establishing a novel tool for inducible gene expression in this bacterium.

Identification of a cytosine methyltransferase that improves transformation efficiency in Methylomonas sp. DH-1

BACKGROUND

Industrial biofuels and other value-added products can be produced from metabolically engineered microorganisms. Methylomonas sp. DH-1 is a candidate platform for bioconversion that uses methane as a carbon source. Although several genetic engineering techniques have been developed to work with Methylomonas sp. DH-1, the genetic manipulation of plasmids remains difficult because of the restriction-modification (RM) system present in the bacteria. Therefore, the RM system in Methylomonas sp. DH-1 must be identified to improve the genetic engineering prospects of this microorganism.

RESULTS

We identified a DNA methylation site, TGGCCA, and its corresponding cytosine methyltransferase for the first time in Methylomonas sp. DH-1 through whole-genome bisulfite sequencing. The methyltransferase was confirmed to methylate the fourth nucleotide of TGGCCA. In general, methylated plasmids exhibited better transformation efficiency under the protection of the RM system than non-methylated plasmids did. As expected, when we transformed Methylomonas sp. DH-1 with plasmid DNA harboring the psy gene, the metabolic flux towards carotenoid increased. The methyltransferase-treated plasmid exhibited an increase in transformation efficiency of 2.5 × 103 CFU/μg (124%). The introduced gene increased the production of carotenoid by 26%. In addition, the methyltransferase-treated plasmid harboring anti-psy sRNA gene exhibited an increase in transformation efficiency by 70% as well. The production of carotenoid was decreased by 40% when the psy gene was translationally repressed by anti-psy sRNA.

CONCLUSIONS

Plasmid DNA methylated by the discovered cytosine methyltransferase from Methylomonas sp. DH-1 had a higher transformation efficiency than non-treated plasmid DNA. The RM system identified in this study may facilitate the plasmid-based genetic manipulation of methanotrophs.

Omics-driven identification and elimination of valerolactam catabolism in Pseudomonas putida KT2440 for increased product titer

Pseudomonas putida is a promising bacterial chassis for metabolic engineering given its ability to metabolize a wide array of carbon sources, especially aromatic compounds derived from lignin. However, this omnivorous metabolism can also be a hindrance when it can naturally metabolize products produced from engineered pathways. Herein we show that P. putida is able to use valerolactam as a sole carbon source, as well as degrade caprolactam. Lactams represent important nylon precursors, and are produced in quantities exceeding one million tons per year (Zhang et al., 2017). To better understand this metabolism we use a combination of Random Barcode Transposon Sequencing (RB-TnSeq) and shotgun proteomics to identify the oplBA locus as the likely responsible amide hydrolase that initiates valerolactam catabolism. Deletion of the oplBA genes prevented P. putida from growing on valerolactam, prevented the degradation of valerolactam in rich media, and dramatically reduced caprolactam degradation under the same conditions. Deletion of oplBA, as well as pathways that compete for precursors L-lysine or 5-aminovalerate, increased the titer of valerolactam from undetectable after 48 h of production to ~90 mg/L. This work may serve as a template to rapidly eliminate undesirable metabolism in non-model hosts in future metabolic engineering efforts.

Enzymatic characterization and molecular mechanism of a novel aspartokinase mutant M372I/T379W from Corynebacterium pekinense

A novel aspartokinase mutant M372I/T379W from Corynebacterium pekinense was constructed by using site-directed mutagenesis. The enzyme was then purified, characterized, and its molecular mechanism was comprehensively analyzed. Compared with wild-type AK, the catalytic activity of M372I/T379W AK was 16.51 fold higher and the optimum temperature increased from 28 to 35 °C. The thermostability of M372I/T379W AK was significantly improved. Microscale thermophoresis analysis indicated that M372I/T379W AK not only weakened the inhibitory effect of Lys, but also had stronger binding force with Asp. Molecular dynamics simulation showed that mutations M372I and T379W could regulate the activity of CpAK through affecting the flexibility of Asp and ATP binding pocket residues and the hydrogen bond between CpAK and Asp. In addition, mutations could affect the relative position of protein domains. The width of the Asp binding pocket entrance gate Arg169-Ala60 of M372I/T379W AK was greater than that in wild-type AK and the CpAK switched from T-state to R-state, which promoted the binding of the enzyme to Asp and improving the catalytic efficiency of this enzyme. These results explain the molecular mechanism of M372I/T379W AK, which will greatly facilitate the rational design of more aspartokinase mutants, with have potential applications in aspartic acid metabolism.

l-lysine production by Bacillus methanolicus: Genome-based mutational analysis and l-lysine secretion engineering

Bacillus methanolicus is a methylotrophic bacterium with an increasing interest in academic research and for biotechnological applications. This bacterium was previously applied for methanol-based production of l-glutamate, l-lysine and the five-carbon diamine cadaverine by wild type, classical mutant and recombinant strains. The genomes of two different l-lysine secreting B. methanolicus classical mutant strains, NOA2#13A52-8A66 and M168-20, were sequenced. We focused on mutational mapping in genes present in l-lysine and other relevant amino acid biosynthetic pathways, as well as in the primary cell metabolism important for precursor supply. In addition to mutations in the aspartate pathway genes dapG, lysA and hom-1, new mutational target genes like alr, proA, proB1, leuC, odhA and pdhD were identified. Surprisingly, no mutations were found in the putative l-lysine transporter gene lysE^MGA3. Inspection of the wild type B. methanolicus strain PB1 genome sequence identified two homologous putative l-lysine transporter genes, lysE^PB1 and lysE2^PB1. The biological role of these putative l-lysine transporter genes, together with the heterologous l-lysine exporter gene lysE^Cg from Corynebacterium glutamicum, were therefore investigated. Our results demonstrated that the titer of secreted l-lysine in B. methanolicus was significantly increased by overexpression of lysE^Cg while overexpression of lysE^MGA3, lysE^PB1 and lysE2^PB1 had no measurable effect.

Finding approximate gene clusters with Gecko 3

Gene-order-based comparison of multiple genomes provides signals for functional analysis of genes and the evolutionary process of genome organization. Gene clusters are regions of co-localized genes on genomes of different species. The rapid increase in sequenced genomes necessitates bioinformatics tools for finding gene clusters in hundreds of genomes. Existing tools are often restricted to few (in many cases, only two) genomes, and often make restrictive assumptions such as short perfect conservation, conserved gene order or monophyletic gene clusters. We present Gecko 3, an open-source software for finding gene clusters in hundreds of bacterial genomes, that comes with an easy-to-use graphical user interface. The underlying gene cluster model is intuitive, can cope with low degrees of conservation as well as misannotations and is complemented by a sound statistical evaluation. To evaluate the biological benefit of Gecko 3 and to exemplify our method, we search for gene clusters in a dataset of 678 bacterial genomes using Synechocystis sp. PCC 6803 as a reference. We confirm detected gene clusters reviewing the literature and comparing them to a database of operons; we detect two novel clusters, which were confirmed by publicly available experimental RNA-Seq data. The computational analysis is carried out on a laptop computer in <40 min.

Mutagenesis of Key Residues in the Binding Center of l-Aspartate-b-Semialdehyde Dehydrogenase from Escherichia coli Enhances Utilization of the Cofactor NAD(H)

L-Aspartate-β-semialdehyde dehydrogenase (ASADH) is a key enzyme in the aspartate pathway. In bacteria, ASADH is highly specific for the cofactor NADP(+) rather than NAD(+). Limited information on cofactor utilization is available, and neither the wild-type protein nor the available mutants could utilize NAD(+) efficiently. In this study, we identified several residues crucial for cofactor utilization by Escherichia coli ASADH (ecASADH) by mutating residues within the cofactor binding center. Among the investigated mutants, ecASADH-Q350N and ecASADH-Q350N/H171A, which exhibited markedly improved NAD(+) utilization, were further investigated by various biochemical approaches and molecular modeling. Relative to the wild type, the two mutants showed approximately 44-fold and 66-fold increases, respectively, in the constant kcat /Km of NAD(+). As desired, they could also utilize NADH efficiently to synthesize l-homoserine in cascade reactions in vitro.

Methanol-based cadaverine production by genetically engineered Bacillus methanolicus strains

Methanol is regarded as an attractive substrate for biotechnological production of value-added bulk products, such as amino acids and polyamines. In the present study, the methylotrophic and thermophilic bacterium Bacillus methanolicus was engineered into a microbial cell factory for the production of the platform chemical 1,5-diaminopentane (cadaverine) from methanol. This was achieved by the heterologous expression of the Escherichia coli genes cadA and ldcC encoding two different lysine decarboxylase enzymes, and by increasing the overall L-lysine production levels in this host. Both CadA and LdcC were functional in B. methanolicus cultivated at 50°C and expression of cadA resulted in cadaverine production levels up to 500 mg l⁻¹ during shake flask conditions. A volume-corrected concentration of 11.3 g l⁻¹ of cadaverine was obtained by high-cell density fed-batch methanol fermentation. Our results demonstrated that efficient conversion of L-lysine into cadaverine presumably has severe effects on feedback regulation of the L-lysine biosynthetic pathway in B. methanolicus. By also investigating the cadaverine tolerance level, B. methanolicus proved to be an exciting alternative host and comparable to the well-known bacterial hosts E. coli and Corynebacterium glutamicum. This study represents the first demonstration of microbial production of cadaverine from methanol.

Methylotrophy in the thermophilic Bacillus methanolicus, basic insights and application for commodity production from methanol

Using methanol as an alternative non-food feedstock for biotechnological production offers several advantages in line with a methanol-based bioeconomy. The Gram-positive, facultative methylotrophic and thermophilic bacterium Bacillus methanolicus is one of the few described microbial candidates with a potential for the conversion of methanol to value-added products. Its capabilities of producing and secreting the commercially important amino acids L-glutamate and L-lysine to high concentrations at 50 °C have been demonstrated and make B. methanolicus a promising target to develop cell factories for industrial-scale production processes. B. methanolicus uses the ribulose monophosphate cycle for methanol assimilation and represents the first example of plasmid-dependent methylotrophy. Recent genome sequencing of two physiologically different wild-type B. methanolicus strains, MGA3 and PB1, accompanied with transcriptome and proteome analyses has generated fundamental new insight into the metabolism of the species. In addition, multiple key enzymes representing methylotrophic and biosynthetic pathways have been biochemically characterized. All this, together with establishment of improved tools for gene expression, has opened opportunities for systems-level metabolic engineering of B. methanolicus. Here, we summarize the current status of its metabolism and biochemistry, available genetic tools, and its potential use in respect to overproduction of amino acids.

Functional characterization of key enzymes involved in L-glutamate synthesis and degradation in the thermotolerant and methylotrophic bacterium Bacillus methanolicus

Bacillus methanolicus wild-type strain MGA3 secretes 59 g/liter(-1) of l-glutamate in fed-batch methanol cultivations at 50°C. We recently sequenced the MGA3 genome, and we here characterize key enzymes involved in l-glutamate synthesis and degradation. One glutamate dehydrogenase (GDH) that is encoded by yweB and two glutamate synthases (GOGATs) that are encoded by the gltAB operon and by gltA2 were found, in contrast to Bacillus subtilis, which has two different GDHs and only one GOGAT. B. methanolicus has a glutamine synthetase (GS) that is encoded by glnA and a 2-oxoglutarate dehydrogenase (OGDH) that is encoded by the odhAB operon. The yweB, gltA, gltB, and gltA2 gene products were purified and characterized biochemically in vitro. YweB has a low Km value for ammonium (10 mM) and a high Km value for l-glutamate (250 mM), and the Vmax value is 7-fold higher for l-glutamate synthesis than for the degradation reaction. GltA and GltA2 displayed similar Km values (1 to 1.4 mM) and Vmax values (4 U/mg) for both l-glutamate and 2-oxoglutarate as the substrates, and GltB had no effect on the catalytic activities of these enzymes in vitro. Complementation assays indicated that GltA and not GltA2 is dependent on GltB for GOGAT activity in vivo. To our knowledge, this is the first report describing the presence of two active GOGATs in a bacterium. In vivo experiments indicated that OGDH activity and, to some degree, GOGAT activity play important roles in regulating l-glutamate production in this organism.

Methylotrophic Bacillus methanolicus encodes two chromosomal and one plasmid born NAD+ dependent methanol dehydrogenase paralogs with different catalytic and biochemical properties

Bacillus methanolicus can utilize methanol as the sole carbon source for growth and it encodes an NAD⁺-dependent methanol dehydrogenase (Mdh), catalyzing the oxidation of methanol to formaldehyde. Recently, the genomes of the B. methanolicus strains MGA3 (ATCC53907) and PB1 (NCIMB13113) were sequenced and found to harbor three different putative Mdh encoding genes, each belonging to the type III Fe-NAD⁺-dependent alcohol dehydrogenases. In each strain, two of these genes are encoded on the chromosome and one on a plasmid; only one chromosomal act gene encoding the previously described activator protein ACT was found. The six Mdhs and the ACT proteins were produced recombinantly in Escherichia coli, purified, and characterized. All Mdhs required NAD⁺ as cosubstrate, were catalytically stimulated by ACT, exhibited a broad and different substrate specificity range and displayed both dehydrogenase and reductase activities. All Mdhs catalyzed the oxidation of methanol; however the catalytic activity for methanol was considerably lower than for most other alcohols tested, suggesting that these enzymes represent a novel class of alcohol dehydrogenases. The kinetic constants for the Mdhs were comparable when acting as pure enzymes, but together with ACT the differences were more pronounced. Quantitative PCR experiments revealed major differences with respect to transcriptional regulation of the paralogous genes. Taken together our data indicate that the repertoire of methanol oxidizing enzymes in thermotolerant bacilli is larger than expected with complex mechanisms involved in their regulation.

Methanol oxidation by temperate soils and environmental determinants of associated methylotrophs

The role of soil methylotrophs in methanol exchange with the atmosphere has been widely overlooked. Methanol can be derived from plant polymers and be consumed by soil microbial communities. In the current study, methanol-utilizing methylotrophs of 14 aerated soils were examined to resolve their comparative diversities and capacities to utilize ambient concentrations of methanol. Abundances of cultivable methylotrophs ranged from 10(6)-10(8) gsoilDW(-1). Methanol dissimilation was measured based on conversion of supplemented (14)C-methanol, and occurred at concentrations down to 0.002 μmol methanol gsoilDW(-1). Tested soils exhibited specific affinities to methanol (a(0)s=0.01 d(-1)) that were similar to those of other environments suggesting that methylotrophs with similar affinities were present. Two deep-branching alphaproteobacterial genotypes of mch responded to the addition of ambient concentrations of methanol (0.6 μmol methanol gsoilDW(-1)) in one of these soils. Methylotroph community structures were assessed by amplicon pyrosequencing of genes of mono carbon metabolism (mxaF, mch and fae). Alphaproteobacteria-affiliated genotypes were predominant in all investigated soils, and the occurrence of novel genotypes indicated a hitherto unveiled diversity of methylotrophs. Correlations between vegetation type, soil pH and methylotroph community structure suggested that plant-methylotroph interactions were determinative for soil methylotrophs.

Genome sequence of thermotolerant Bacillus methanolicus: features and regulation related to methylotrophy and production of L-lysine and L-glutamate from methanol

Bacillus methanolicus can utilize methanol as its sole carbon and energy source, and the scientific interest in this thermotolerant bacterium has focused largely on exploring its potential as a biocatalyst for the conversion of methanol into L-lysine and L-glutamate. We present here the genome sequences of the important B. methanolicus model strain MGA3 (ATCC 53907) and the alternative wild-type strain PB1 (NCIMB13113). The physiological diversity of these two strains was demonstrated by a comparative fed-batch methanol cultivation displaying highly different methanol consumption and respiration profiles, as well as major differences in their L-glutamate production levels (406 mmol liter(-1) and 11 mmol liter(-1), respectively). Both genomes are small (ca 3.4 Mbp) compared to those of other related bacilli, and MGA3 has two plasmids (pBM19 and pBM69), while PB1 has only one (pBM20). In particular, we focus here on genes representing biochemical pathways for methanol oxidation and concomitant formaldehyde assimilation and dissimilation, the important phosphoenol pyruvate/pyruvate anaplerotic node, the tricarboxylic acid cycle including the glyoxylate pathway, and the biosynthetic pathways for L-lysine and L-glutamate. Several unique findings were made, including the discovery of three different methanol dehydrogenase genes in each of the two B. methanolicus strains, and the genomic analyses were accompanied by gene expression studies. Our results provide new insight into a number of peculiar physiological and metabolic traits of B. methanolicus and open up possibilities for system-level metabolic engineering of this bacterium for the production of amino acids and other useful compounds from methanol.

Expression of recombinant green fluorescent protein in Bacillus methanolicus

Microbial biocatalysts are used in a wide range of industries to produce large scale quantities of proteins, amino acids, and commodity chemicals. While the majority of these processes use glucose or other low-cost sugars as the substrate, Bacillus methanolicus is one example of a biocatalyst that has shown sustained growth on methanol as a carbon source at elevated temperature (50-53°C optimum) resulting in reduced feed and utility costs. Specifically, the complete chemical process enabled by this approach takes methane from natural gas, and following a low-cost conversion to methanol, can be used for the production of high value products. In this study, production of recombinant green fluorescent protein (GFPuv) by B. methanolicus is explored. A plasmid was constructed that incorporates the methanol dehydrogenase (mdh) promoter of B. methanolicus MGA3 together with the GFPuv gene. The plasmid, pNW33N, was shown to be effective for expression in other Bacillus strains, although not previously in B. methanolicus. A published electroporation protocol for transformation of B. methanolicus was modified to result in expression of GFP using plasmid pNW33N-mdh-GFPuv (pNmG). Transformation was confirmed by both agarose gel electrophoresis and by observation of green fluorescence under UV light exposure. The mass yield of cells and protein were measured in shake flask experiments. The optimum concentration of methanol for protein production was found to be at 200 mM. Higher concentrations than 200 mM resulted in slightly higher biomass production but lower amounts of recombinant protein.