Computer-aided design and implementation of efficient biosynthetic pathways to produce high added-value products derived from tyrosine in Escherichia coli.

Developing efficient bioprocesses requires selecting the best biosynthetic pathways, which can be challenging and time-consuming due to the vast amount of data available in databases and literature. The extension of the shikimate pathway for the biosynthesis of commercially attractive molecules often involves promiscuous enzymes or lacks well-established routes. To address these challenges, we developed a computational workflow integrating enumeration/retrosynthesis algorithms, a toolbox for pathway analysis, enzyme selection tools, and a gene discovery pipeline, supported by manual curation and literature review. Our focus has been on implementing biosynthetic pathways for tyrosine-derived compounds, specifically L-3,4-dihydroxyphenylalanine (L-DOPA) and dopamine, with significant applications in health and nutrition. We selected one pathway to produce L-DOPA and two different pathways for dopamine-one already described in the literature and a novel pathway. Our goal was either to identify the most suitable gene candidates for expression in Escherichia coli for the known pathways or to discover innovative pathways. Although not all implemented pathways resulted in the accumulation of target compounds, in our shake-flask experiments we achieved a maximum L-DOPA titer of 0.71 g/L and dopamine titers of 0.29 and 0.21 g/L for known and novel pathways, respectively. In the case of L-DOPA, we utilized, for the first time, a mutant version of tyrosinase from Ralstonia solanacearum. Production of dopamine via the known biosynthesis route was accomplished by coupling the L-DOPA pathway with the expression of DOPA decarboxylase from Pseudomonas putida, resulting in a unique biosynthetic pathway never reported in literature before. In the context of the novel pathway, dopamine was produced using tyramine as the intermediate compound. To achieve this, tyrosine was initially converted into tyramine by expressing TDC from Levilactobacillus brevis, which, in turn, was converted into dopamine through the action of the enzyme encoded by ppoMP from Mucuna pruriens. This marks the first time that an alternative biosynthetic pathway for dopamine has been validated in microbes. These findings underscore the effectiveness of our computational workflow in facilitating pathway enumeration and selection, offering the potential to uncover novel biosynthetic routes, thus paving the way for other target compounds of biotechnological interest.

Bioproduction of methylated phenylpropenes and isoeugenol in Escherichia coli.

Phenylpropenes are a class of natural products that are synthesised by a vast range of plant species and hold considerable promise in the flavour and fragrance industries. Many in vitro studies have been carried out to elucidate and characterise the enzymes responsible for the production of these volatile compounds. However, there is a scarcity of studies demonstrating the in vivo production of phenylpropenes in microbial cell factories. In this study, we engineered Escherichia coli to produce methylchavicol, methyleugenol and isoeugenol from their respective phenylacrylic acid precursors. We achieved this by extending and modifying a previously optimised heterologous pathway for the biosynthesis of chavicol and eugenol. We explored the potential of six S-adenosyl l-methionine (SAM)-dependent O-methyltransferases to produce methylchavicol and methyleugenol from chavicol and eugenol, respectively. Additionally, we examined two isoeugenol synthases for the production of isoeugenol from coniferyl acetate. The best-performing strains in this study were able to achieve titres of 13 mg L-1 methylchavicol, 59 mg L-1 methyleugenol and 361 mg L-1 isoeugenol after feeding with their appropriate phenylacrylic acid substrates. We were able to further increase the methyleugenol titre to 117 mg L-1 by supplementation with methionine to facilitate SAM recycling. Moreover, we report the biosynthesis of methylchavicol and methyleugenol from l-tyrosine through pathways involving six and eight enzymatic steps, respectively.

SelenzymeRF: updated enzyme suggestion software for unbalanced biochemical reactions.

Selenzyme is a retrobiosynthesis tool that suggests candidate enzymes for user query reactions. Enzyme suggestions are based on identical reactions, as well as similar reactions, since enzymes are often capable of promiscuous substrate binding. Selenzyme is a user-friendly, widely used web-tool for ranking enzymes based on reaction similarity and additional features, including the phylogenetic distance between the source species of the enzyme and the intended host. While Selenzyme has proved invaluable in assisting with enzyme selection for known reactions, as well as many novel or orphan reactions, weaknesses have been exposed in its ability to rank functionally related enzymes. Within this update, we introduce a new reaction similarity scoring algorithm, which is used in conjunction with the previous similarity calculation, to improve the accuracy of enzyme suggestions based on non-identical similar reactions, across a range of EC reaction classes. This allows enzymes to be suggested for reactions not found within the database, even if the reaction is unbalanced. A database update was also carried out, to ensure that reaction and enzyme knowledge remains current. This update can be accessed at http://selenzymeRF.synbiochem.co.uk/.

Carboxylic acid reductase-dependent biosynthesis of eugenol and related allylphenols.

BACKGROUND: (Hydroxy)cinnamyl alcohols and allylphenols, including coniferyl alcohol and eugenol, are naturally occurring aromatic compounds widely utilised in pharmaceuticals, flavours, and fragrances. Traditionally, the heterologous biosynthesis of (hydroxy)cinnamyl alcohols from (hydroxy)cinnamic acids involved CoA-dependent activation of the substrate. However, a recently explored alternative pathway involving carboxylic acid reductase (CAR) has proven efficient in generating the (hydroxy)cinnamyl aldehyde intermediate without the need for CoA activation. In this study, we investigated the application of the CAR pathway for whole-cell bioconversion of a range of (hydroxy)cinnamic acids into their corresponding (hydroxy)cinnamyl alcohols. Furthermore, we sought to extend the pathway to enable the production of a variety of allylphenols and allylbenzene. RESULTS: By screening the activity of several heterologously expressed enzymes in crude cell lysates, we identified the combination of Segniliparus rugosus CAR (SrCAR) and Medicago sativa cinnamyl alcohol dehydrogenase (MsCAD2) as the most efficient enzymatic cascade for the two-step reduction of ferulic acid to coniferyl alcohol. To optimise the whole-cell bioconversion in Escherichia coli, we implemented a combinatorial approach to balance the gene expression levels of SrCAR and MsCAD2. This optimisation resulted in a coniferyl alcohol yield of almost 100%. Furthermore, we extended the pathway by incorporating coniferyl alcohol acyltransferase and eugenol synthase, which allowed for the production of eugenol with a titre of up to 1.61 mM (264 mg/L) from 3 mM ferulic acid. This improvement in titre surpasses previous achievements in the field employing a CoA-dependent coniferyl alcohol biosynthesis pathway. Our study not only demonstrated the successful utilisation of the CAR pathway for the biosynthesis of diverse (hydroxy)cinnamyl alcohols, such as p-coumaryl alcohol, caffeyl alcohol, cinnamyl alcohol, and sinapyl alcohol, from their corresponding (hydroxy)cinnamic acid precursors but also extended the pathway to produce allylphenols, including chavicol, hydroxychavicol, and methoxyeugenol. Notably, the microbial production of methoxyeugenol from sinapic acid represents a novel achievement. CONCLUSION: The combination of SrCAR and MsCAD2 enzymes offers an efficient enzymatic cascade for the production of a wide array of (hydroxy)cinnamyl alcohols and, ultimately, allylphenols from their respective (hydroxy)cinnamic acids. This expands the range of value-added molecules that can be generated using microbial cell factories and creates new possibilities for applications in industries such as pharmaceuticals, flavours, and fragrances. These findings underscore the versatility of the CAR pathway, emphasising its potential in various biotechnological applications.

Expanding flavone and flavonol production capabilities in Escherichia coli.

Flavones and flavonols are important classes of flavonoids with nutraceutical and pharmacological value, and their production by fermentation with recombinant microorganisms promises to be a scalable and economically favorable alternative to extraction from plant sources. Flavones and flavonols have been produced recombinantly in a number of microorganisms, with Saccharomyces cerevisiae typically being a preferred production host for these compounds due to higher yields and titers of precursor compounds, as well as generally improved ability to functionally express cytochrome P450 enzymes without requiring modification to improve their solubility. Recently, a rapid prototyping platform has been developed for high-value compounds in E. coli, and a number of gatekeeper (2S)-flavanones, from which flavones and flavonols can be derived, have been produced to high titers in E. coli using this platform. In this study, we extended these metabolic pathways using the previously reported platform to produce apigenin, chrysin, luteolin and kaempferol from the gatekeeper flavonoids naringenin, pinocembrin and eriodictyol by the expression of either type-I flavone synthases (FNS-I) or type-II flavone synthases (FNS-II) for flavone biosynthesis, and by the expression of flavanone 3-dioxygenases (F3H) and flavonol synthases (FLS) for the production of the flavonol kaempferol. In our best-performing strains, titers of apigenin and kaempferol reached 128 mg L-1 and 151 mg L-1 in 96-DeepWell plates in cultures supplemented with an additional 3 mM tyrosine, though titers for chrysin (6.8 mg L-1) from phenylalanine, and luteolin (5.0 mg L-1) from caffeic acid were considerably lower. In strains with upregulated tyrosine production, apigenin and kaempferol titers reached 80.2 mg L-1 and 42.4 mg L-1 respectively, without the further supplementation of tyrosine beyond the amount present in the rich medium. Notably, the highest apigenin, chrysin and luteolin titers were achieved with FNS-II enzymes, suggesting that cytochrome P450s can show competitive performance compared with non-cytochrome P450 enzymes in prokaryotes for the production of flavones.

Microbial production of the plant flavanone hesperetin from caffeic acid.

OBJECTIVE: Hesperetin is an important O-methylated flavonoid produced by citrus fruits and of potential pharmaceutical relevance. The microbial biosynthesis of hesperetin could be a viable alternative to plant extraction, as plant extracts often yield complex mixtures of different flavonoids making it challenging to isolate pure compounds. In this study, hesperetin was produced from caffeic acid in the microbial host Escherichia coli. We combined a previously optimised pathway for the biosynthesis of the intermediate flavanone eriodictyol with a combinatorial library of plasmids expressing three candidate flavonoid O-methyltransferases. Moreover, we endeavoured to improve the position specificity of CCoAOMT7, a flavonoid O-methyltransferase from Arabidopsis thaliana that has been demonstrated to O-methylate eriodictyol in both the para- and meta-position, thus leading to a mixture of hesperetin and homoeriodictyol. RESULTS: The best performing flavonoid O-methyltransferase in our screen was found to be CCoAOMT7, which could produce up to 14.6 mg/L hesperetin and 3.8 mg/L homoeriodictyol from 3 mM caffeic acid in E. coli 5-alpha. Using a platform for enzyme engineering that scans the mutational space of selected key positions, predicting their structures using homology modelling and inferring their potential catalytic improvement using docking simulations, we were able to identify a CCoAOMT7 mutant with a two-fold higher position specificity for hesperetin. The mutant's catalytic activity, however, was considerably diminished. Our findings suggest that hesperetin can be created from central carbon metabolism in E. coli following the introduction of a caffeic acid biosynthesis pathway.

TFBMiner: A User-Friendly Command Line Tool for the Rapid Mining of Transcription Factor-Based Biosensors.

Transcription factors responsive to small molecules are essential elements in synthetic biology designs. They are often used as genetically encoded biosensors with applications ranging from the detection of environmental contaminants and biomarkers to microbial strain engineering. Despite our efforts to expand the space of compounds that can be detected using biosensors, the identification and characterization of transcription factors and their corresponding inducer molecules remain labor- and time-intensive tasks. Here, we introduce TFBMiner, a new data mining and analysis pipeline that enables the automated and rapid identification of putative metabolite-responsive transcription factor-based biosensors (TFBs). This user-friendly command line tool harnesses a heuristic rule-based model of gene organization to identify both gene clusters involved in the catabolism of user-defined molecules and their associated transcriptional regulators. Ultimately, biosensors are scored based on how well they fit the model, providing wet-lab scientists with a ranked list of candidates that can be experimentally tested. We validated the pipeline using a set of molecules for which TFBs have been reported previously, including sensors responding to sugars, amino acids, and aromatic compounds, among others. We further demonstrated the utility of TFBMiner by identifying a biosensor for S-mandelic acid, an aromatic compound for which a responsive transcription factor had not been found previously. Using a combinatorial library of mandelate-producing microbial strains, the newly identified biosensor was able to distinguish between low- and high-producing strain candidates. This work will aid in the unraveling of metabolite-responsive microbial gene regulatory networks and expand the synthetic biology toolbox to allow for the construction of more sophisticated self-regulating biosynthetic pathways.

A biobank of small cell lung cancer CDX models elucidates inter- and intratumoral phenotypic heterogeneity.

Although small cell lung cancer (SCLC) is treated as a homogeneous disease, biopsies and preclinical models reveal heterogeneity in transcriptomes and morphology. SCLC subtypes were recently defined by neuroendocrine transcription factor (NETF) expression. Circulating-tumor-cell-derived explant models (CDX) recapitulate donor patients' tumor morphology, diagnostic NE marker expression and chemotherapy responses. We describe a biobank of 38 CDX models, including six CDX pairs generated pretreatment and at disease progression revealing complex intra- and intertumoral heterogeneity. Transcriptomic analysis confirmed three of four previously described subtypes based on ASCL1, NEUROD1 and POU2F3 expression and identified a previously unreported subtype based on another NETF, ATOH1. We document evolution during disease progression exemplified by altered MYC and NOTCH gene expression, increased 'variant' cell morphology, and metastasis without strong evidence of epithelial to mesenchymal transition. This CDX biobank provides a research resource to facilitate SCLC personalized medicine.

Using set theory to reduce redundancy in pathway sets.

BACKGROUND: The consolidation of pathway databases, such as KEGG, Reactome and ConsensusPathDB, has generated widespread biological interest, however the issue of pathway redundancy impedes the use of these consolidated datasets. Attempts to reduce this redundancy have focused on visualizing pathway overlap or merging pathways, but the resulting pathways may be of heterogeneous sizes and cover multiple biological functions. Efforts have also been made to deal with redundancy in pathway data by consolidating enriched pathways into a number of clusters or concepts. We present an alternative approach, which generates pathway subsets capable of covering all of genes presented within either pathway databases or enrichment results, generating substantial reductions in redundancy. RESULTS: We propose a method that uses set cover to reduce pathway redundancy, without merging pathways. The proposed approach considers three objectives: removal of pathway redundancy, controlling pathway size and coverage of the gene set. By applying set cover to the ConsensusPathDB dataset we were able to produce a reduced set of pathways, representing 100% of the genes in the original data set with 74% less redundancy, or 95% of the genes with 88% less redundancy. We also developed an algorithm to simplify enrichment data and applied it to a set of enriched osteoarthritis pathways, revealing that within the top ten pathways, five were redundant subsets of more enriched pathways. Applying set cover to the enrichment results removed these redundant pathways allowing more informative pathways to take their place. CONCLUSION: Our method provides an alternative approach for handling pathway redundancy, while ensuring that the pathways are of homogeneous size and gene coverage is maximised. Pathways are not altered from their original form, allowing biological knowledge regarding the data set to be directly applicable. We demonstrate the ability of the algorithms to prioritise redundancy reduction, pathway size control or gene set coverage. The application of set cover to pathway enrichment results produces an optimised summary of the pathways that best represent the differentially regulated gene set.

Mapping biological process relationships and disease perturbations within a pathway network.

Molecular interaction networks are routinely used to map the organization of cellular function. Edges represent interactions between genes, proteins, or metabolites. However, in living cells, molecular interactions are dynamic, necessitating context-dependent models. Contextual information can be integrated into molecular interaction networks through the inclusion of additional molecular data, but there are concerns about completeness and relevance of this data. We developed an approach for representing the organization of human cellular processes using pathways as the nodes in a network. Pathways represent spatial and temporal sets of context-dependent interactions, generating a high-level network when linked together, which incorporates contextual information without the need for molecular interaction data. Analysis of the pathway network revealed linked communities representing functional relationships, comparable to those found in molecular networks, including metabolism, signaling, immunity, and the cell cycle. We mapped a range of diseases onto this network and find that pathways associated with diseases tend to be functionally connected, highlighting the perturbed functions that result in disease phenotypes. We demonstrated that disease pathways cluster within the network. We then examined the distribution of cancer pathways and showed that cancer pathways tend to localize within the signaling, DNA processes and immune modules, although some cancer-associated nodes are found in other network regions. Altogether, we generated a high-confidence functional network, which avoids some of the shortcomings faced by conventional molecular models. Our representation provides an intuitive functional interpretation of cellular organization, which relies only on high-quality pathway and Gene Ontology data. The network is available at https://data.mendeley.com/datasets/3pbwkxjxg9/1.

Disentangling the multigenic and pleiotropic nature of molecular function.

BACKGROUND: Biological processes at the molecular level are usually represented by molecular interaction networks. Function is organised and modularity identified based on network topology, however, this approach often fails to account for the dynamic and multifunctional nature of molecular components. For example, a molecule engaging in spatially or temporally independent functions may be inappropriately clustered into a single functional module. To capture biologically meaningful sets of interacting molecules, we use experimentally defined pathways as spatial/temporal units of molecular activity. RESULTS: We defined functional profiles of Saccharomyces cerevisiae based on a minimal set of Gene Ontology terms sufficient to represent each pathway's genes. The Gene Ontology terms were used to annotate 271 pathways, accounting for pathway multi-functionality and gene pleiotropy. Pathways were then arranged into a network, linked by shared functionality. Of the genes in our data set, 44% appeared in multiple pathways performing a diverse set of functions. Linking pathways by overlapping functionality revealed a modular network with energy metabolism forming a sparse centre, surrounded by several denser clusters comprised of regulatory and metabolic pathways. Signalling pathways formed a relatively discrete cluster connected to the centre of the network. Genetic interactions were enriched within the clusters of pathways by a factor of 5.5, confirming the organisation of our pathway network is biologically significant. CONCLUSIONS: Our representation of molecular function according to pathway relationships enables analysis of gene/protein activity in the context of specific functional roles, as an alternative to typical molecule-centric graph-based methods. The pathway network demonstrates the cooperation of multiple pathways to perform biological processes and organises pathways into functionally related clusters with interdependent outcomes.