With its ability to catabolize a wide variety of carbon sources and a growing engineering toolkit, Pseudomonas putida KT2440 is emerging as an important chassis organism for metabolic engineering. Despite advances in our understanding of the organism, many gaps remain in our knowledge of the genetic basis of its metabolic capabilities. The gaps are particularly noticeable in our understanding of both fatty acid and alcohol catabolism, where many paralogs putatively coding for similar enzymes coexist, making biochemical assignment via sequence homology difficult. To rapidly assign function to the enzymes responsible for these metabolisms, we leveraged random barcode transposon sequencing (RB-Tn-Seq). Global fitness analyses of transposon libraries grown on 13 fatty acids and 10 alcohols produced strong phenotypes for hundreds of genes. Fitness data from mutant pools grown on fatty acids of varying chain lengths indicated specific enzyme substrate preferences and enabled us to hypothesize that DUF1302/DUF1329 family proteins potentially function as esterases. From the data, we also postulate catabolic routes for the two biogasoline molecules isoprenol and isopentanol, which are catabolized via leucine metabolism after initial oxidation and activation with coenzyme A (CoA). Because fatty acids and alcohols may serve as both feedstocks and final products of metabolic-engineering efforts, the fitness data presented here will help guide future genomic modifications toward higher titers, rates, and yields. IMPORTANCE To engineer novel metabolic pathways into P. putida, a comprehensive understanding of the genetic basis of its versatile metabolism is essential. Here, we provide functional evidence for the putative roles of hundreds of genes involved in the fatty acid and alcohol metabolism of the bacterium. These data provide a framework facilitating precise genetic changes to prevent product degradation and to channel the flux of specific pathway intermediates as desired.

Mitchell G. Thompson,a,b,c Matthew R. Incha,a,b,d Allison N. Pearson,a,b,d Matthias Schmidt,a,b,e William A. Sharpless,a,b Christopher B. Eiben,a,b,f Pablo Cruz-Morales,a,b,g Jacquelyn M. Blake-Hedges,a,b,h Yuzhong Liu,a,b Catharine A. Adams,b Robert W. Haushalter,a,b Rohith N. Krishna,a,b Patrick Lichtner,a,b Lars M. Blank,e Aindrila Mukhopadhyay,a,b AdamM. Deutschbauer,d,i Patrick M. Shih,a,b,c,j Jay D. Keaslinga,b,f,k,l,m,n

Author(s): Thompson, Mitchell G; Incha, Matthew R; Pearson, Allison N; Schmidt, Matthias; Sharpless, William A; Eiben, Christopher B; Cruz-Morales, Pablo; Blake-Hedges, Jacquelyn M; Liu, Yuzhong; Adams, Catharine A; Haushalter, Robert W; Krishna, Rohith N; Lichtner, Patrick; Blank, Lars M; Mukhopadhyay, Aindrila; Deutschbauer, Adam M; Shih, Patrick M; Keasling, Jay D

Mycobacterium sp. strain JC1 DSM3803 is one of the few known bacteria predicted to possess the xylulose monophosphate (XuMP) pathway of C1 carbon assimilation. The draft genome is 7,921,603 bp with a GC content of 66.88% and will allow more in depth investigation of the bacterium’s unique metabolisms.

Terpenes constitute the largest class of natural products with over 55,000 compounds with versatile applications including drugs and biofuels. Introducing structural modifications to terpenes through metabolic engineering is an efficient and sustainable way to improve their properties. Here, we report the optimization of the lepidopteran mevalonate (LMVA) pathway towards the efficient production of isopentenyl pyrophosphate (IPP) analogs as terpene precursors. First, we linked the LMVA pathway to NudB, a promiscuous phosphatase, resulting in the production of the six-carbon analog of 3-methyl-3-buten-1-ol (isoprenol), 3-ethyl-3-buten-1-ol (C6-isoprenol). Using C6-isoprenol as the final product, we then engineered the LMVA pathway by redirecting its upstream portion from a thiolase-dependent pathway to a beta-oxidation pathway. The beta-oxidation LMVA pathway transforms valeric acid, a platform chemical that can be produced from biomass, into C6-isoprenol at a titer of 110.3 mg/L, improved from 5.5 mg/L by the thiolase LMVA pathway, which used propionic acid as a feedstock. Knockout of the E. coli endogenous thiolase genes further improved the C6-isoprenol titer to 390 mg/L, implying efficient production of homo isopentenyl pyrophosphate (HIPP). The beta-oxidation LMVA-NudB pathway also converts butanoic acid and hexanoic acid into isoprenol and isoprenol's seven-carbon analog, 3-propyl-3-buten1-ol (C7-isoprenol), respectively, suggesting the beta-oxidation LMVA pathway produces IPP and C7-IPP from the corresponding fatty acids. Fuel property tests revealed the longer chain isoprenol analogs have lower water solubilities, similar or higher energy densities, and comparable research octane number (RON) boosting effects to isopentenols. This work not only optimizes the LMVA pathway, setting the basis for homoterpene biosynthesis to expand terpene chemical space, but provides an efficient pathway to produce isoprenol analogs as nextgeneration biofuels from sustainable feedstocks.

Despite intensive study, plant lysine catabolism beyond the 2-oxoadipate (2OA) intermediate remains unvalidated. Recently we described a missing step in the D-lysine catabolism of Pseudomonas putida in which 2OA is converted to D-2-hydroxyglutarate (2HG) via hydroxyglutarate synthase (HglS), a DUF1338 family protein. Here we solve the structure of HglS to 1.1 Å resolution in substrate-free form and in complex with 2OA. We propose a successive decarboxylation and intramolecular hydroxylation mechanism forming 2HG in a Fe(II)- and O 2 -dependent manner. Specificity is mediated by a single arginine, highly conserved across most DUF1338 proteins. An Arabidopsis thaliana HglS homolog coexpresses with known lysine catabolism enzymes, and mutants show phenotypes consistent with disrupted lysine catabolism. Structural and biochemical analysis of Oryza sativa homolog FLO7 reveals identical activity to HglS despite low sequence identity. Our results suggest DUF1338-containing enzymes catalyze the same biochemical reaction, exerting the same physiological function across bacteria and eukaryotes.

With its ability to catabolize a wide variety of carbon sources and a growing engineering toolkit, Pseudomonas putida KT2440 is emerging as an important chassis organism for metabolic engineering. Despite advances in our understanding of this organism, many gaps remain in our knowledge of the genetic basis of its metabolic capabilities. These gaps are particularly noticeable in our understanding of both fatty acid and alcohol catabolism, where many paralogs putatively coding for similar enzymes co-exist making biochemical assignment via sequence homology difficult. To rapidly assign function to the enzymes responsible for these metabolisms, we leveraged Random Barcode Transposon Sequencing (RB-TnSeq). Global fitness analyses of transposon libraries grown on 13 fatty acids and 10 alcohols produced strong phenotypes for hundreds of genes. Fitness data from mutant pools grown on varying chain length fatty acids indicated specific enzyme substrate preferences, and enabled us to hypothesize that DUF1302/DUF1329 family proteins potentially function as esterases. From the data we also postulate catabolic routes for the two biogasoline molecules isoprenol and isopentanol, which are catabolized via leucine metabolism after initial oxidation and activation with CoA. Because fatty acids and alcohols may serve as both feedstocks or final products of metabolic engineering efforts, the fitness data presented here will help guide future genomic modifications towards higher titers, rates, and yields. IMPORTANCE To engineer novel metabolic pathways into P. putida , a comprehensive understanding of the genetic basis of its versatile metabolism is essential. Here we provide functional evidence for the putative roles of hundreds of genes involved in the fatty acid and alcohol metabolism of this bacterium. These data provide a framework facilitating precise genetic changes to prevent product degradation and channel the flux of specific pathway intermediates as desired. ### Competing Interest Statement J.D.K. has financial interests in Amyris, Lygos, Demetrix, Napigen, Maple Bio, and Apertor Labs. C.B.E has a financial interest in Perlumi Chemicals.

Climate change necessitates the development of CO2 neutral or negative routes to chemicals currently produced from fossil carbon. In this paper we demonstrate a pathway from the renewable resource glucose to next generation biofuel isopentanol by pairing the isovaleryl-CoA biosynthesis pathway from Myxococcus xanthus and a butyryl-CoA reductase from Clostridium acetobutylicum. The best plasmid and Escherichia coli strain combination makes 80.50 ± 8.08 (SD) mg/L of isopentanol after 36 h under microaerobic conditions with an oleyl alcohol overlay. In addition, the system also shows a strong preference for isopentanol production over prenol in microaerobic conditions. Finally, the pathway requires zero adenosine triphosphate and can be paired theoretically with nonoxidative glycolysis, the combination being redox balanced from glucose thus avoiding unnecessary carbon loss as CO2. These pathway properties make the isovaleryl-CoA pathway an attractive isopentanol production route for further optimization.

Lepidoptera (butterflies and moths) make the six-carbon compounds homoisopentenyl pyrophosphate (HIPP) and homodimethylallyl pyrophosphate (HDMAPP) that are incorporated into 16, 17, and 18 carbon farnesyl pyrophosphate (FPP) analogues. In this work we heterologously expressed the lepidopteran modified mevalonate pathway, a propionyl-CoA ligase, and terpene cyclases in E. coli to produce several novel terpenes containing 16 carbons. Changing the terpene cyclase generated different novel terpene product profiles. To further validate the new compounds we confirmed 13C propionate was incorporated, and that the masses and fragmentation patterns were consistent with novel 16 carbon terpenes by GC-QTOF. On the basis of the available farnesyl pyrophosphate analogues lepidoptera produce, this approach should greatly expand the reachable biochemical space with applications in areas where terpenes have traditionally found uses.

P. putida lysine metabolism can produce multiple commodity chemicals, conferring great biotechnological value. Despite much research, the connection of lysine catabolism to central metabolism in P. putida remained undefined. Here, we used random barcode transposon sequencing to fill the gaps of lysine metabolism in P. putida . We describe a route of 2-oxoadipate (2OA) catabolism, which utilizes DUF1338-containing protein P. putida 5260 (PP_5260) in bacteria. Despite its prevalence in many domains of life, DUF1338-containing proteins have had no known biochemical function. We demonstrate that PP_5260 is a metalloenzyme which catalyzes an unusual route of decarboxylation of 2OA to d -2-hydroxyglutarate ( d -2HG). Our screen also identified a recently described novel glutarate metabolic pathway. We validate previous results and expand the understanding of glutarate hydroxylase CsiD by showing that can it use either 2OA or 2KG as a cosubstrate. Our work demonstrated that biological novelty can be rapidly identified using unbiased experimental genetics and that RB-TnSeq can be used to rapidly validate previous results.

Advances in retooling microorganisms have enabled bioproduction of 'drop-in' biofuels, fuels that are compatible with existing spark-ignition, compression-ignition, and gas-turbine engines. As the majority of petroleum consumption in the United States consists of gasoline (47%), diesel fuel and heating oil (21%), and jet fuel (8%), 'drop-in' biofuels that replace these petrochemical sources are particularly attractive. In this review, we discuss the application of aldehyde decarbonylases to produce gasoline substitutes from fatty acid products, a recently crystallized reductase that could hydrogenate jet fuel precursors from terpene synthases, and the exquisite control of polyketide synthases to produce biofuels with desired physical properties (e.g., lower freezing points). With our increased understanding of biosynthetic logic of metabolic pathways, we discuss the unique advantages of fatty acid, terpene, and polyketide synthases for the production of bio-based gasoline, diesel and jet fuel.

•Fatty acid, isoprenoid and polyketide pathways for production of drop-in biofuels.•FAS pathways produce alkanes via newly discovered fatty aldehyde decarbonylase.•Crystallized geranylgeranyl reductase hydrogenates isoprenoid-based double bonds.•Polyketide synthases for generation of designer biofuels.

We describe a computationally designed enzyme, formolase (FLS), which catalyzes the carboligation of three one-carbon formaldehyde molecules into one three-carbon dihydroxyacetone molecule. The existence of FLS enables the design of a new carbon fixation pathway, the formolase pathway, consisting of a small number of thermodynamically favorable chemical transformations that convert formate into a three-carbon sugar in central metabolism. The formolase pathway is predicted to use carbon more efficiently and with less backward flux than any naturally occurring one-carbon assimilation pathway. When supplemented with enzymes carrying out the other steps in the pathway, FLS converts formate into dihydroxyacetone phosphate and other central metabolites in vitro. These results demonstrate how modern protein engineering and design tools can facilitate the construction of a completely new biosynthetic pathway.

Baker and colleagues have already shown that protein folding can be turned into an online game played by nonscientists. Now, Foldit players tackle the problem of increasing an enzyme's catalytic activity—with promising results. Computational enzyme design holds promise for the production of renewable fuels, drugs and chemicals. De novo enzyme design has generated catalysts for several reactions, but with lower catalytic efficiencies than naturally occurring enzymes1,2,3,4. Here we report the use of game-driven crowdsourcing to enhance the activity of a computationally designed enzyme through the functional remodeling of its structure. Players of the online game Foldit5,6 were challenged to remodel the backbone of a computationally designed bimolecular Diels-Alderase3 to enable additional interactions with substrates. Several iterations of design and characterization generated a 24-residue helix-turn-helix motif, including a 13-residue insertion, that increased enzyme activity >18-fold. X-ray crystallography showed that the large insertion adopts a helix-turn-helix structure positioned as in the Foldit model. These results demonstrate that human creativity can extend beyond the macroscopic challenges encountered in everyday life to molecular-scale design problems.

Past anthrax attacks in the United States have highlighted the need for improved measures against bioweapons. The virulence of anthrax stems from the shielding properties of the Bacillus anthracis poly-γ-d-glutamic acid capsule. In the presence of excess CapD, a B. anthracis γ-glutamyl transpeptidase, the protective capsule is degraded, and the immune system can successfully combat infection. Although CapD shows promise as a next generation protein therapeutic against anthrax, improvements in production, stability, and therapeutic formulation are needed. In this study, we addressed several of these problems through computational protein engineering techniques. We show that circular permutation of CapD improved production properties and dramatically increased kinetic thermostability. At 45 °C, CapD was completely inactive after 5 min, but circularly permuted CapD remained almost entirely active after 30 min. In addition, we identify an amino acid substitution that dramatically decreased transpeptidation activity but not hydrolysis. Subsequently, we show that this mutant had a diminished capsule degradation activity, suggesting that CapD catalyzes capsule degradation through a transpeptidation reaction with endogenous amino acids and peptides in serum rather than hydrolysis.