Computational-guided discovery of UDP-glycosyltransferases for lauryl glucoside production using engineered E. coli

Abstract

Defining suitable enzymes for reaction steps in novel synthetic pathways is crucial for developing microbial cell factories for non-natural products. Here, we developed a computational workflow to identify C12 alcohol-active UDP-glycosyltransferases. The workflow involved three steps: (1) assembling initial candidates of putative UDP-glycosyltransferases, (2) refining selection by examining conserved regions, and (3) 3D structure prediction and molecular docking. Genomic sequences from Candida, Pichia, Rhizopus, and Thermotoga, known for lauryl glucoside synthesis via whole-cell biocatalysis, were screened. Out of 240 predicted glycosyltransferases, 8 candidates annotated as glycosyltransferases were selected after filtering out those with signal peptides and identifying conserved UDP-glycosyltransferase regions. These proteins underwent 3D structure prediction and molecular docking with 1-dodecanol. RO3G, a candidate from Rhizopus delemar RA 99–880 with a relatively high ChemPLP fitness score, was selected and expressed in Escherichia coli BL21 (DE3). It was further characterized using a feeding experiment with 1-dodecanol. Results confirmed that the RO3G-expressing strain could convert 1-dodecanol to lauryl glucoside, as quantified by HPLC and identified by targeted LC-MS. Monitoring the growth and fermentation profiles of the engineered strain revealed that RO3G expression did not affect cell growth. Interestingly, acetate, a major fermentation product, was reduced in the RO3G-expressing strain compared to the GFP-expressing strain, suggesting a redirection of flux from acetate to other pathways. Overall, this work presents a successful workflow for discovering UDP-glycosyltransferase enzymes with confirmed activity toward 1-dodecanol for lauryl glucoside production.

Graphical abstract

Keywords

UDP-glycosyltransferase / Lauryl glucoside / Computational method / 1-Dodecanol / Engineered E. coli / Novel biosynthetic pathway / Genome mining

References

[1] Bhatwa A, Wang W, Hassan YI et al (2021) Challenges associated with the formation of recombinant protein inclusion bodies in Escherichia coli and strategies to address them for industrial applications. Front Bioeng Biotechnol 9:1–18. https://doi.org/10.3389/fbioe.2021.630551

[2] BIOVIA DS (2017) ‘Discovery studio’, Discovery studio modeling environment, Release 2017

[3] Blum M, Chang HY, Chuguransky S et al (2021) The InterPro protein families and domains database: 20 years on. Nucleic Acids Res 49:D344–D354. https://doi.org/10.1093/nar/gkaa977

[4] Brockhausen I, Hu B, Liu B et al (2008) Characterization of two β-1,3-glucosyltransferases from Escherichia coli serotypes O56 and O152. J Bacteriol 190:4922–4932. https://doi.org/10.1128/JB.00160-08

[5] Chen G, Seukep AJ, Guo M (2020) Recent advances in molecular docking for the research and discovery of potential marine drugs. Mar Drugs 18(11). https://doi.org/10.3390/md18110545

[6] Cho JS, Kim GB, Eun H et al (2022) Designing microbial cell factories for the production of chemicals. JACS Au 2(8):1781–1799. https://doi.org/10.1021/jacsau.2c00344

[7] El-Sharkawy SH (1996) Production of useful chemicals from sunflower oil by microbial biotransformation with Rhizopus stolonifer NRRL 1478. Australas. Biotechnol 6:8–12.

[8] Fiume MM, Heldreth B, Bergfeld WF et al (2013) Safety assessment of decyl glucoside and other alkyl glucosides as used in cosmetics. Int J Toxicol 32:22S–48S. https://doi.org/10.1177/1091581813497764

[9] Gharabli H, Della Gala V, Welner DH (2023) The function of UDP-glycosyltransferases in plants and their possible use in crop protection. Biotechnol Adv 67:108182. https://doi.org/10.1016/j.biotechadv.2023.108182

[10] Gueguen Y, Chemardin P, Pommares P. Enzymatic synthesis of dodecyl.β-d-glucopyranoside catalyzed by Candida Molischiana 35M5N β-glucosidase. Synthesis, 1995, 8524(95): 263-267

[11] He B, Bai X, Tan Y et al (2022) Glycosyltransferases: mining, engineering and applications in biosynthesis of glycosylated plant natural products. Synth Syst Biotechnol 7:602–620. https://doi.org/10.1016/j.synbio.2022.01.001

[12] Hsieh S, Wang J, Lai Y et al (2018) Production of 1-Dodecanol, 1-Tetradecanol, and 1,12- dodecanediol through whole-cell biotransformation in Escherichia coli. Appl Environ Microbiol 84(4)

[13] Kaewkod T, Bovonsombut S, Tragoolpua Y (2019) Efficacy of kombucha obtained from green, oolongand black teas on inhibition of pathogenic bacteria, antioxidation, and toxicity on colorectal cancer cell line. Microorganisms 7(12):1–18. https://doi.org/10.3390/microorganisms7120700

[14] Lidholt K, Fjelstad M, Jann K, Lindahl U (1994) Substrate specificities of glycosyltransferases involved in formation of heparin precursor and E. coli K5 capsular polysaccharides. Carbohydr Res 255:87–101. https://doi.org/10.1016/S0008-6215(00)90972-8

[15] Ma LJ, Ibrahim AS, Skory C et al (2009) Genomic analysis of the basal lineage fungus Rhizopus oryzae reveals a whole-genome duplication. PLoS Geneti 5(7). https://doi.org/10.1371/journal.pgen.1000549

[16] Meech R, et al. The UDP-Glycosyltransferase (UGT) superfamily: new members, new functions, and novel paradigms. Physiol Rev, 2019, 99(2): 1153-1222 Available at

[17] Mirdita M, Schütze K, Moriwaki Y et al (2022) ColabFold: making protein folding accessible to all. Nat Methods 19:679–682. https://doi.org/10.1038/s41592-022-01488-1

[18] Mohammadi T, Karczmarek A, Crouvoisier M et al (2007) The essential peptidoglycan glycosyltransferase MurG forms a complex with proteins involved in lateral envelope growth as well as with proteins involved in cell division in Escherichia coli. Mol Microbiol 65:1106–1121. https://doi.org/10.1111/j.1365-2958.2007.05851.x

[19] Pandurangan AP, Stahlhacke J, Oates ME et al (2019) The SUPERFAMILY 2.0 database: a significant proteome update and a new webserver. Nucleic Acids Res 47:D490–D494. https://doi.org/10.1093/nar/gky1130

[20] Papanikolaou A, Chatzikonstantinou AV, Zarafeta D et al (2023) Substrate specificity of the highly thermostable esterase EstDZ3. ChemBioChem 24(5). https://doi.org/10.1002/cbic.202200642

[21] Pei T, Yan M, Li T et al (2022) Characterization of UDP-glycosyltransferase family members reveals how major flavonoid glycoside accumulates in the roots of Scutellaria baicalensis. BMC Genomics 23:1–14. https://doi.org/10.1186/s12864-022-08391-1

[22] Pozzo T et al (2010) Structural and functional analyses of β-glucosidase 3B from Thermotoga neapolitana: a thermostable three-domain representative of glycoside hydrolase 3. J Mol Biol 397(3):724–739. https://doi.org/10.1016/j.jmb.2010.01.072

[23] Pulungan I, Sugijanto K, Prawita A. Validation Gas Chromatography-Fid Method for Analysis of Ethanol Content in Vinegar Metode Validasi Gas Chromatography-Fid Untuk Analisis Konten Ethanol Pada Cuka. J Halal Prod Res, 2018, 01(02): 22-32

[24] Ross J, et al. Higher plant glycosyltransferases. Genome Biol, 2001, 2(2): 1-6

[25] Sattayawat P, Sofian Yunus I, Jones PR (2020) Bioderivatization as a concept for renewable production of chemicals that are toxic or poorly soluble in the liquid phase. Proc Natl Acad Sci U S A 117:1404–1413. https://doi.org/10.1073/pnas.1914069117

[26] Shahbaaz M, Hassan MI, Ahmad F (2013) Functional annotation of conserved hypothetical proteins from Haemophilus influenzae Rd KW20. PLoS ONE 8(12). https://doi.org/10.1371/journal.pone.0084263

[27] Shahbaaz M, Ahmad F, Imtaiyaz Hassan M (2015) Structure-based functional annotation of putative conserved proteins having lyase activity from Haemophilus influenzae. 3 Biotech 5(3):317–336. https://doi.org/10.1007/s13205-014-0231-z

[28] Sorrell S, Speirs J, Bentley R et al (2010) Global oil depletion: a review of the evidence. Energy Policy 38:5290–5295. https://doi.org/10.1016/j.enpol.2010.04.046

[29] Taujale R, Venkat A, Huang LC et al (2020) Deep evolutionary analysis reveals the design principles of fold a glycosyltransferases Elife 9:1–24. https://doi.org/10.7554/eLife.54532

[30] Teufel F, Almagro Armenteros JJ, Johansen AR et al (2022) SignalP 6.0 predicts all five types of signal peptides using protein language models. Nat Biotechnol 40:1023–1025. https://doi.org/10.1038/s41587-021-01156-3

[31] Tripathi D, Kapoor A, Bulbul et al (2024) Functional annotation of Candida albicans hypothetical proteins: a bioinformatics approach. Arch Microbiol 206:1–11. https://doi.org/10.1007/s00203-024-03840-9

[32] Uttarotai T, Mukjang N, Chaisoung N et al (2022) Putative protein discovery from microalgal genomes as a synthetic biology protein library for heavy metal bio-removal. Biology (Basel) 11:1226. https://doi.org/10.3390/biology11081226

[33] Verdonk ML, Cole JC, Hartshorn MJ et al (2003) Improved protein-ligand docking using GOLD. Proteins Struct Funct Genet 52:609–623. https://doi.org/10.1002/prot.10465

[34] Vivijs B, Haberbeck LU, Mfortaw Mbong VB et al (2015) Formate hydrogen lyase mediates stationary-phase deacidification and increases survival during sugar fermentation in acetoin-producing enterobacteria. Front Microbiol 6:1–11. https://doi.org/10.3389/fmicb.2015.00150

[35] Wang D, Wang J, Shi Y et al (2020) Elucidation of the complete biosynthetic pathway of the main triterpene glycosylation products of Panax notoginseng using a synthetic biology platform. Metab Eng 61:131–140. https://doi.org/10.1016/j.ymben.2020.05.007

[36] Wang M, Ji Q, Lai B et al (2023a) Structure-function and engineering of plant UDP-glycosyltransferase. Comput Struct Biotechnol J 21:5358–5371. https://doi.org/10.1016/j.csbj.2023.10.046

[37] Wang HT, Wang ZL, Chen K et al (2023b) Insights into the missing apiosylation step in flavonoid apiosides biosynthesis of Leguminosae plants. Nat Commun 14. https://doi.org/10.1038/s41467-023-42393-1

[38] Xiao Z, Grosse S, Bergeron H, Lau PCK (2014) Cloning and characterization of the first GH10 and GH11 xylanases from Rhizopus oryzae.Appl Microbiol Biotechnol 98:8211–8222. https://doi.org/10.1007/s00253-014-5741-4

[39] Xie K, Dou X, Chen R et al (2017) Two Novel Fungal Phenolic UDP glycosyltransferases from Absidia coerulea and Rhizopus japonicus. Appl Environ Microbiol 83(8):1–14

[40] Xue Y et al (2021) Engineering Thermotoga maritima β-glucosidase for improved alkyl glycosides synthesis by site-directed mutagenesis. J Ind Microbiol Biotechnol 48(5–6):1–11. https://doi.org/10.1093/jimb/kuab031

[41] Yin Q, Shen G, Chang Z et al (2017) Involvement of three putative glucosyltransferases from the UGT72 family in flavonol glucoside/rhamnoside biosynthesis in Lotus japonicus seeds. J Exp Bot 68:597–612. https://doi.org/10.1093/jxb/erw420

[42] Younis Rather M, Mishra S, Verma V, Chand S (2012) Biotransformation of methyl-β-d-glucopyranoside to higher chain alkyl glucosides by cell bound β-glucosidase of Pichia etchellsii. Bioresour Technol 107:287–294. https://doi.org/10.1016/j.biortech.2011.11.061

[43] Yunus IS, Jones PR (2018) Photosynthesis-dependent biosynthesis of medium chain-length fatty acids and alcohols. Metab Eng 49:59–68. https://doi.org/10.1016/j.ymben.2018.07.015

[44] Zou L, Geng X, Li Z, Li T (2023) Design of highly active substrates using molecular docking for microbial transglutaminase detection. RSC Adv 13:5259–5265. https://doi.org/10.1039/d2ra06467g

Funding

Office of the Permanent Secretary, Ministry of Higher Education, Science, Research and Innovation, Thailand(RGNS 64 – 069)