Effects of bowel cleansing on the composition of the gut microbiota in inflammatory bowel disease patients and healthy controls.

Background: In patients with inflammatory bowel disease (IBD), Crohn's disease (CD), and ulcerative colitis (UC), numerous cases of exacerbations could be observed after colonoscopy, raising the possible pathogenetic effect of colonic microbiota alterations in IBD flare. Objectives: We aimed to investigate the changes in the fecal microbiota composition in IBD patients influenced by the bowel preparation with sodium picosulfate. Design: We enrolled patients with IBD undergoing bowel preparation for colonoscopy in the prospective cohort study. The control group (Con) comprised non-IBD patients who underwent colonoscopy. Clinical data, blood, and stool samples were collected before colonoscopy (timepoint A), 3 days later (timepoint B), and 4 weeks later (timepoint C). Methods: Disease activity and gut microbiota changes were assessed at each timepoint. Fecal microbiota structure - at family level - was determined by sequencing the V4 region of the 16S rRNA gene. Statistical analysis included differential abundance analysis and Mann-Whitney tests. Results: Forty-one patients (9 CD, 13 UC, and 19 Con) were included. After bowel preparation, alpha diversity was lower in the CD group than in the UC (p = 0.01) and Con (p = 0.02) groups at timepoint B. Alpha diversity was significantly higher in the UC group than in the CD and Con (p = 0.03) groups at timepoint C. Beta diversity difference differed between the IBD and Con (p = 0.001) groups. Based on the differential abundance analysis, the Clostridiales family was increased, whereas the Bifidobacteriaceae family was decreased in CD patients compared to the Con at timepoint B. Conclusions: Bowel preparation may change the fecal microbial composition in IBD patients, which may have a potential role in disease exacerbation after bowel cleansing.

Accessing the In Vivo Efficiency of Clinically Isolated Phages against Uropathogenic and Invasive Biofilm-Forming Escherichia coli Strains for Phage Therapy.

Escherichia coli is one of the most common members of the intestinal microbiota. Many of its strains are associated with various inflammatory infections, including urinary or gut infections, especially when displaying antibiotic resistance or in patients with suppressed immune systems. According to recent reports, the biofilm-forming potential of E. coli is a crucial factor for its increased resistance against antibiotics. To overcome the limitations of using antibiotics against resistant E. coli strains, the world is turning once more towards bacteriophage therapy, which is becoming a promising candidate amongst the current personalized approaches to target different bacterial infections. Although matured and persistent biofilms pose a serious challenge to phage therapy, they can still become an effective alternative to antibiotic treatment. Here, we assess the efficiency of clinically isolated phages in phage therapy against representative clinical uropathogenic and invasive biofilm-forming E. coli strains. Our results demonstrate that irrespective of host specificity, bacteriophages producing clear plaques with a high burst size, and exhibiting depolymerizing activity, are good candidates against biofilm-producing E. coli pathogens as verified from our in vitro and in vivo experiments using Galleria mellonella where survival was significantly increased for phage-therapy-treated larvae.

Bacterial RNA virus MS2 exposure increases the expression of cancer progression genes in the LNCaP prostate cancer cell line.

Bacteriophages effectively counteract diverse bacterial infections, and their ability to treat most types of cancer has been explored using phage engineering or phage-virus hybrid platforms. In the present study, it was demonstrated that the bacteriophage MS2 can affect the expression of genes associated with the proliferation and survival of LNCaP prostate epithelial cells. LNCaP cells were exposed to bacteriophage MS2 at a concentration of 1×107 plaque forming units/ml for 24-48 h. After exposure, various cellular parameters, including cell viability, morphology, and changes in gene expression, were examined. MS2 affected cell viability adversely, reducing viability by 25% in the first 4 h of treatment; however, cell viability recovered within 24-48 h. Similarly, the AKT, androgen receptor, integrin α5, integrin ß1, MAPK1, MAPK3, STAT3, and peroxisome proliferator-activated receptor-γ coactivator 1α genes, which are involved in various normal cellular processes and tumor progression, were significantly upregulated, whereas the expression levels of HSP90, ITGB5, ITGB3, HSP27, ITGAV, and PI3K genes were unchanged. Therefore, based on viability and gene expression changes, bacteriophage MS2 severely impaired LNCaP cells by reducing anchorage-dependent survival and androgen signaling. A caveolin-mediated endocytosis mechanism for MS2-mediated signaling in prostate cancer cells was proposed based on reports involving bacteriophages T4, M13, and MS2, and their interactions with LNCaP and PC3 cell lines.

SpyPhage: A Cell-Free TXTL Platform for Rapid Engineering of Targeted Phage Therapies.

The past decade has seen the emergence of multidrug resistant pathogens as a leading cause of death worldwide, reigniting interest in the field of phage therapy. Modern advances in the genetic engineering of bacteriophages have enabled several useful results including host range alterations, constitutive lytic growth, and control over phage replication. However, the slow licensing process of genetically modified organisms clearly inhibits the rapid therapeutic application of novel engineered variants necessary to fight mutant pathogens that emerge throughout the course of a pandemic. As a solution to this problem, we propose the SpyPhage system where a "scaffold" bacteriophage is engineered to incorporate a SpyTag moiety on its capsid head to enable rapid postsynthetic modification of their surfaces with SpyCatcher-fused therapeutic proteins. As a proof of concept, through CRISPR/Cas-facilitated phage engineering and whole genome assembly, we targeted a SpyTag capsid fusion to K1F, a phage targeting the pathogenic strain Escherichia coli K1. We demonstrate for the first time the cell-free assembly and decoration of the phage surface with two alternative fusion proteins, SpyCatcher-mCherry-EGF and SpyCatcher-mCherry-Rck, both of which facilitate the endocytotic uptake of the phages by a urinary bladder epithelial cell line. Overall, our work presents a cell-free phage production pipeline for the generation of multiple phenotypically distinct phages with a single underlying "scaffold" genotype. These phages could become the basis of next-generation phage therapies where the knowledge-based engineering of numerous phage variants would be quickly achievable without the use of live bacteria or the need to repeatedly license novel genetic alterations.

A novel Bacteroides metallo-ß-lactamase (MBL) and its gene (crxA) in Bacteroides xylanisolvens revealed by genomic sequencing and functional analysis.

OBJECTIVES: We sought to characterize the carbapenem resistance mechanism of Bacteroides xylanisolvens 14880, an imipenem-resistant strain from Germany, and assess its prevalence. METHODS: Antimicrobial susceptibilities were determined using agar dilution or Etest methodology and specific imipenemase activity was detected. The genomic sequence of B. xylanisolvens 14880 was determined and analysed for antibiotic resistance genes and genomic islands. We also used gene transfer to a carbapenem susceptible host, along with 5'-RACE, conventional PCR with capillary sequencing and RT-PCR-based screening. RESULTS: B. xylanisolvens 14880 displayed resistance to carbapenems and produced high specific imipenemase activity. Its genomic sequence was 6.1 Mbp and a class B1 ß-lactamase gene (termed crxA) was detected in it. crxA was carried on a putative genomic island with insertion sequence (IS) elements and a putative GNAT (Gcn5-like acetyltransferase) toxin gene. Promoter localization by 5'-RACE and gene targeting to an imipenem-susceptible Bacteroides host indicated that it is activated by an IS1380-like IS element and it can confer carbapenem resistance. The PCR screening of Bacteroides strains showed that crxA was specific to B. xylanisolvens with a carriage rate of 16.7%. CONCLUSIONS: B. xylanisolvens strains can harbour a carbapenem resistance gene, which has many similarities to the 'cfiA system': metallo-ß-lactamase (MBL), IS element activation, carriage of a GNAT toxin gene, specific for a unique Bacteroides species with a significant prevalence.

inPOSE: A Flexible Toolbox for Chromosomal Cloning and Amplification of Bacterial Transgenes.

Cloning the genes and operons encoding heterologous functions in bacterial hosts is now almost exclusively carried out using plasmid vectors. This has multiple drawbacks, including the need for constant selection and variation in copy numbers. The chromosomal integration of transgenes has always offered a viable alternative; however, to date, it has been of limited use due to its tedious nature and often being limited to a single copy. We introduce here a strategy that uses bacterial insertion sequences, which are the simplest autonomous transposable elements to insert and amplify genetic cargo into a bacterial chromosome. Transgene insertion can take place either as transposition or homologous recombination, and copy number amplification is achieved using controlled copy-paste transposition. We display the successful use of IS1 and IS3 for this purpose in Escherichia coli cells using various selection markers. We demonstrate the insertion of selectable genes, an unselectable gene and a five-gene operon in up to two copies in a single step. We continue with the amplification of the inserted cassette to double-digit copy numbers within two rounds of transposase induction and selection. Finally, we analyze the stability of the cloned genetic constructs in the lack of selection and find it to be superior to all investigated plasmid-based systems. Due to the ubiquitous nature of transposable elements, we believe that with proper design, this strategy can be adapted to numerous other bacterial species.

ProtTrans: Toward Understanding the Language of Life Through Self-Supervised Learning.

Computational biology and bioinformatics provide vast data gold-mines from protein sequences, ideal for Language Models (LMs) taken from Natural Language Processing (NLP). These LMs reach for new prediction frontiers at low inference costs. Here, we trained two auto-regressive models (Transformer-XL, XLNet) and four auto-encoder models (BERT, Albert, Electra, T5) on data from UniRef and BFD containing up to 393 billion amino acids. The protein LMs (pLMs) were trained on the Summit supercomputer using 5616 GPUs and TPU Pod up-to 1024 cores. Dimensionality reduction revealed that the raw pLM-embeddings from unlabeled data captured some biophysical features of protein sequences. We validated the advantage of using the embeddings as exclusive input for several subsequent tasks: (1) a per-residue (per-token) prediction of protein secondary structure (3-state accuracy Q3=81%-87%); (2) per-protein (pooling) predictions of protein sub-cellular location (ten-state accuracy: Q10=81%) and membrane versus water-soluble (2-state accuracy Q2=91%). For secondary structure, the most informative embeddings (ProtT5) for the first time outperformed the state-of-the-art without multiple sequence alignments (MSAs) or evolutionary information thereby bypassing expensive database searches. Taken together, the results implied that pLMs learned some of the grammar of the language of life. All our models are available through https://github.com/agemagician/ProtTrans.

Frequency, composition and mobility of Escherichia coli-derived transposable elements in holdings of plasmid repositories.

By providing the scientific community with uniform and standardized resources of consistent quality, plasmid repositories play an important role in enabling scientific reproducibility. Plasmids containing insertion sequence elements (IS elements) represent a challenge from this perspective, as they can change the plasmid structure and function. In this study, we conducted a systematic analysis of a subset of plasmid stocks distributed by plasmid repositories (The Arabidopsis Biological Resource Center and Addgene) which carry unintended integrations of bacterial mobile genetic elements. The integration of insertion sequences was most often found in, but not limited to, pBR322-derived vectors, and did not affect the function of the specific plasmids. In certain cases, the entire stock was affected, but the majority of the stocks tested contained a mixture of the wild-type and the mutated plasmids, suggesting that the acquisition of IS elements likely occurred after the plasmids were acquired by the repositories. However, comparison of the sequencing results of the original samples revealed that some plasmids already carried insertion mutations at the time of donation. While an extensive BLAST analysis of 47 877 plasmids sequenced from the Addgene repository uncovered IS elements in only 1.12%, suggesting that IS contamination is not widespread, further tests showed that plasmid integration of IS elements can propagate in conventional Escherichia coli hosts over a few tens of generations. Use of IS-free E. coli hosts prevented the emergence of IS insertions as well as that of small indels, suggesting that the use of IS-free hosts by donors and repositories could help limit unexpected and unwanted IS integrations into plasmids.

Analysing Parallel Strategies to Alter the Host Specificity of Bacteriophage T7.

The recognition and binding of host bacteria by bacteriophages is most often enabled by a highly specific receptor-ligand type of interaction, with the receptor-binding proteins (RBPs) of phages being the primary determinants of host specificity. Specifically modifying the RBPs could alter or extend the host range of phages otherwise exhibiting desired phenotypic properties. This study employed two different strategies to reprogram T7 phages ordinarily infecting commensal K12 Escherichia coli strains to infect pathogen-associated K1-capsule-expressing strains. The strategies were based on either plasmid-based homologous recombination or bacteriophage recombineering using electroporated DNA (BRED). Our work pursued the construction of two genetic designs: one replacing the gp17 gene of T7, the other replacing gp11, gp12, and gp17 of T7 with their K1F counterparts. Both strategies displayed successful integration of the K1F sequences into the T7 genome, detected by PCR screening. Multiple methods were utilised to select or enrich for chimeric phages incorporating the K1F gp17 alone, including trxA, host-specificity, and CRISPR-Cas-based selection. Irrespective of the selection method, the above strategy yielded poorly reproducible phage propagation on the new host, indicating that the chimeric phage was less fit than the wild type and could not promote continual autonomous reproduction. Chimeric phages obtained from BRED incorporating gp11-12 and gp17, however, all displayed infection in a 2-stage pattern, indicating the presence of both K1F and T7 phenotypes. This study shows that BRED can be used as a tool to quickly access the potential of new RBP constructs without the need to engineer sustainably replicating phages. Additionally, we show that solely repurposing the primary RBP is, in some cases, insufficient to produce a viable chimeric phage.

Deuterium-depletion has no significant impact on the mutation rate of Escherichia coli, deuterium abundance therefore has a probabilistic, not deterministic effect on spontaneous mutagenesis.

Deuterium (D), the second most abundant isotope of hydrogen is present in natural waters at an approximate concentration of 145-155 ppm (ca. 1.5E-4 atom/atom). D is known to influence various biological processes due to its physical and chemical properties, which significantly differ from those of hydrogen. For example, increasing D-concentration to >1000-fold above its natural abundance has been shown to increase the frequency of genetic mutations in several species. An interesting deterministic hypothesis, formulated with the intent of explaining the mechanism of D-mutagenicity is based on the calculation that the theoretical probability of base pairs to comprise two adjacent D-bridges instead of H-bridges is 2.3E-8, which is equal to the mutation rate of certain species. To experimentally challenge this hypothesis, and to infer the mutagenicity of D present at natural concentrations, we investigated the effect of a nearly 100-fold reduction of D concentration on the bacterial mutation rate. Using fluctuation tests, we measured the mutation rate of three Escherichia coli genes (cycA, ackA and galK) in media containing D at either <2 ppm or 150 ppm concentrations. Out of 15 pair-wise fluctuation analyses, nine indicated a significant decrease, while three marked the significant increase of the mutation/culture value upon D-depletion. Overall, growth in D-depleted minimal medium led to a geometric mean of 0.663-fold (95% confidence interval: 0.483-0.911) change in the mutation rate. This falls nowhere near the expected 10,000-fold reduction, indicating that in our bacterial systems, the effect of D abundance on the formation of point mutations is not deterministic. In addition, the combined results did not display a statistically significant change in the mutation/culture value, the mutation rate or the mutant frequency upon D-depletion. The potential mutagenic effect of D present at natural concentrations on E. coli is therefore below the limit of detection using the indicated methods.

Myopia-26, the female-limited form of early-onset high myopia, occurring in a European family.

BACKGROUND: Female-limited early-onset high myopia, also called Myopia-26 is a rare monogenic disorder characterized by severe short sightedness starting in early childhood and progressing to blindness potentially by the middle ages. Despite the X-linked locus of the mutated ARR3 gene, the disease paradoxically affects females only, with males being asymptomatic carriers. Previously, this disease has only been observed in Asian families and has not gone through detailed investigation concerning collateral symptoms or pathogenesis. RESULTS: We found a large Hungarian family displaying female-limited early-onset high myopia. Whole exome sequencing of two individuals identified a novel nonsense mutation (c.214C>T, p.Arg72*) in the ARR3 gene. We carried out basic ophthalmological testing for 18 family members, as well as detailed ophthalmological examination (intraocular pressure, axial length, fundus appearance, optical coherence tomography, visual field- testing) as well as colour vision- and electrophysiology tests (standard and multifocal electroretinography, pattern electroretinography and visual evoked potentials) for eight individuals. Ophthalmological examinations did not reveal any signs of cone dystrophy as opposed to animal models. Electrophysiology and colour vision tests similarly did not evidence a general cone system alteration, rather a central macular dysfunction affecting both the inner and outer (postreceptoral and receptoral) retinal structures in all patients with ARR3 mutation. CONCLUSIONS: This is the first description of a Caucasian family displaying Myopia-26. We present two hypotheses that could potentially explain the pathomechanism of this disease.

A single plasmid based CRISPR interference in Synechocystis 6803 - A proof of concept.

We developed a simple method to apply CRISPR interference by modifying an existing plasmid pCRISPathBrick containing the native S. pyogenes CRISPR assembly for Synechocystis PCC6803 and named it pCRPB1010. The technique presented here using deadCas9 is easier to implement for gene silencing in Synechocystis PCC6803 than other existing techniques as it circumvents the genome integration and segregation steps thereby significantly shortens the construction of the mutant strains. We executed CRISPR interference against well characterized photosynthetic genes to get a clear phenotype to validate the potential of pCRPB1010 and presented the work as a "proof of concept". Targeting the non-template strand of psbO gene resulted in decreased amount of PsbO and 50% decrease in oxygen evolution rate. Targeting the template strand of psbA2 and psbA3 genes encoding the D1 subunit of photosystem II (PSII) using a single spacer against the common sequence span of the two genes, resulted in full inhibition of both genes, complete abolition of D1 protein synthesis, complete loss of oxygen evolution as well as photoautotrophic growth arrest. This is the first report of a single plasmid based, completely lesion free and episomal expression and execution of CRISPR interference in Synechocystis PCC6803.

CRISPR-interference-based modulation of mobile genetic elements in bacteria.

Spontaneous mutagenesis of synthetic genetic constructs by mobile genetic elements frequently results in the rapid loss of engineered functions. Previous efforts to minimize such mutations required the exceedingly time-consuming manipulation of bacterial chromosomes and the complete removal of insertional sequences (ISes). To this aim, we developed a single plasmid-based system (pCRIS) that applies CRISPR-interference to inhibit the transposition of bacterial ISes. pCRIS expresses multiple guide RNAs to direct inactivated Cas9 (dCas9) to simultaneously silence IS1, IS3, IS5 and IS150 at up to 38 chromosomal loci in Escherichia coli, in vivo. As a result, the transposition rate of all four targeted ISes dropped to negligible levels at both chromosomal and episomal targets. Most notably, pCRIS, while requiring only a single plasmid delivery performed within a single day, provided a reduction of IS-mobility comparable to that seen in genome-scale chromosome engineering projects. The fitness cost of multiple IS-knockdown, detectable in flask-and-shaker systems was readily outweighed by the less frequent inactivation of the transgene, as observed in green fluorescent protein (GFP)-overexpression experiments. In addition, global transcriptomics analysis revealed only minute alterations in the expression of untargeted genes. Finally, the transposition-silencing effect of pCRIS was easily transferable across multiple E. coli strains. The plasticity and robustness of our IS-silencing system make it a promising tool to stabilize bacterial genomes for synthetic biology and industrial biotechnology applications.

Engineered K1F bacteriophages kill intracellular Escherichia coli K1 in human epithelial cells.

Bacterial infections can be treated with bacteriophages that show great specificity towards their bacterial host and can be genetically modified for different applications. However, whether and how bacteriophages can kill intracellular bacteria in human cells remains elusive. Here, using CRISPR/Cas selection, we have engineered a fluorescent bacteriophage specific for E. coli K1, a nosocomial pathogen responsible for urinary tract infections, neonatal meningitis and sepsis. By confocal and live microscopy, we show that engineered bacteriophages K1F-GFP and E. coli EV36-RFP bacteria displaying the K1 capsule, enter human cells via phagocytosis. Importantly, we show that bacteriophage K1F-GFP efficiently kills intracellular E. coli EV36-RFP in T24 human urinary bladder epithelial cells. Finally, we provide evidence that bacteria and bacteriophages are degraded by LC3-associated phagocytosis and xenophagy.

Genome-Wide Abolishment of Mobile Genetic Elements Using Genome Shuffling and CRISPR/Cas-Assisted MAGE Allows the Efficient Stabilization of a Bacterial Chassis.

The ideal bacterial chassis provides a simplified, stable and predictable host environment for synthetic biological circuits. Mutability and evolution can, however, compromise stability, leading to deterioration of artificial genetic constructs. By eliminating certain sources of instability, these undesired genetic changes can be mitigated. Specifically, deletion of prophages and insertion sequences, nonessential constituents of bacterial genomes, has been shown to be beneficial in cellular and genetic stabilization. Here, we sought to establish a rapid methodology to improve the stability of microbial hosts. The novel workflow involves genome shuffling between a mobile genetic element-free strain and the target cell, and subsequent rounds of CRISPR/Cas-assisted MAGE on multiplex targets. The power and speed of the procedure was demonstrated on E. coli BL21(DE3), a host routinely used for plasmid-based heterologous protein expression. All 9 prophages and 50 insertion elements were efficiently deleted or inactivated. Together with additional targeted manipulations (e.g., inactivation of error-prone DNA-polymerases), the changes resulted in an improved bacterial host with a hybrid (harboring segments of K-12 DNA), 9%-downsized and clean genome. The combined capacity of phage-mediated generalized transduction and CRISPR/Cas-selected MAGE offers a way for rapid, large scale editing of bacterial genomes.

System-level genome editing in microbes.

The release of the first complete microbial genome sequences at the end of the past century opened the way for functional genomics and systems-biology to uncover the genetic basis of various phenotypes. The surge of available sequence data facilitated the development of novel genome editing techniques for system-level analytical studies. Recombineering allowed unprecedented throughput and efficiency in microbial genome editing and the recent discovery and widespread use of RNA-guided endonucleases offered several further perspectives: (i) previously recalcitrant species became editable, (ii) the efficiency of recombineering could be elevated, and as a result (iii) diverse genomic libraries could be generated more effectively. Supporting recombineering by RNA-guided endonucleases has led to success stories in metabolic engineering, but their use for system-level analysis is mostly unexplored. For the full exploitation of opportunities that are offered by the genome editing proficiency, future development of large scale analytical procedures is also vitally needed.

Indispensability of Horizontally Transferred Genes and Its Impact on Bacterial Genome Streamlining.

Why are certain bacterial genomes so small and compact? The adaptive genome streamlining hypothesis posits that selection acts to reduce genome size because of the metabolic burden of replicating DNA. To reveal the impact of genome streamlining on cellular traits, we reduced the Escherichia coli genome by up to 20% by deleting regions which have been repeatedly subjects of horizontal transfer in nature. Unexpectedly, horizontally transferred genes not only confer utilization of specific nutrients and elevate tolerance to stresses, but also allow efficient usage of resources to build new cells, and hence influence fitness in routine and stressful environments alike. Genome reduction affected fitness not only by gene loss, but also by induction of a general stress response. Finally, we failed to find evidence that the advantage of smaller genomes would be due to a reduced metabolic burden of replicating DNA or a link with smaller cell size. We conclude that as the potential energetic benefit gained by deletion of short genomic segments is vanishingly small compared with the deleterious side effects of these deletions, selection for reduced DNA synthesis costs is unlikely to shape the evolution of small genomes.

A Sense of Balance: Experimental Investigation and Modeling of a Malonyl-CoA Sensor in Escherichia coli.

Production of value-added chemicals in microorganisms is regarded as a viable alternative to chemical synthesis. In the past decade, several engineered pathways producing such chemicals, including plant secondary metabolites in microorganisms have been reported; upscaling their production yields, however, was often challenging. Here, we analyze a modular device designed for sensing malonyl-CoA, a common precursor for both fatty acid and flavonoid biosynthesis. The sensor can be used either for high-throughput pathway screening in synthetic biology applications or for introducing a feedback circuit to regulate production of the desired chemical. Here, we used the sensor to compare the performance of several predicted malonyl-CoA-producing pathways, and validated the utility of malonyl-CoA reductase and malonate-CoA transferase for malonyl-CoA biosynthesis. We generated a second-order dynamic linear model describing the relation of the fluorescence generated by the sensor to the biomass of the host cell representing a filter/amplifier with a gain that correlates with the level of induction. We found the time constants describing filter dynamics to be independent of the level of induction but distinctively clustered for each of the production pathways, indicating the robustness of the sensor. Moreover, by monitoring the effect of the copy-number of the production plasmid on the dose-response curve of the sensor, we managed to coarse-tune the level of pathway expression to maximize malonyl-CoA synthesis. In addition, we provide an example of the sensor's use in analyzing the effect of inducer or substrate concentrations on production levels. The rational development of models describing sensors, supplemented with the power of high-throughput optimization provide a promising potential for engineering feedback loops regulating enzyme levels to maximize productivity yields of synthetic metabolic pathways.

Validation of RetroPath, a computer-aided design tool for metabolic pathway engineering.

Metabolic engineering has succeeded in biosynthesis of numerous commodity or high value compounds. However, the choice of pathways and enzymes used for production was many times made ad hoc, or required expert knowledge of the specific biochemical reactions. In order to rationalize the process of engineering producer strains, we developed the computer-aided design (CAD) tool RetroPath that explores and enumerates metabolic pathways connecting the endogenous metabolites of a chassis cell to the target compound. To experimentally validate our tool, we constructed 12 top-ranked enzyme combinations producing the flavonoid pinocembrin, four of which displayed significant yields. Namely, our tool queried the enzymes found in metabolic databases based on their annotated and predicted activities. Next, it ranked pathways based on the predicted efficiency of the available enzymes, the toxicity of the intermediate metabolites and the calculated maximum product flux. To implement the top-ranking pathway, our procedure narrowed down a list of nine million possible enzyme combinations to 12, a number easily assembled and tested. One round of metabolic network optimization based on RetroPath output further increased pinocembrin titers 17-fold. In total, 12 out of the 13 enzymes tested in this work displayed a relative performance that was in accordance with its predicted score. These results validate the ranking function of our CAD tool, and open the way to its utilization in the biosynthesis of novel compounds.

Computer-aided design for metabolic engineering.

The development and application of biotechnology-based strategies has had a great socio-economical impact and is likely to play a crucial role in the foundation of more sustainable and efficient industrial processes. Within biotechnology, metabolic engineering aims at the directed improvement of cellular properties, often with the goal of synthesizing a target chemical compound. The use of computer-aided design (CAD) tools, along with the continuously emerging advanced genetic engineering techniques have allowed metabolic engineering to broaden and streamline the process of heterologous compound-production. In this work, we review the CAD tools available for metabolic engineering with an emphasis, on retrosynthesis methodologies. Recent advances in genetic engineering strategies for pathway implementation and optimization are also reviewed as well as a range of bionalytical tools to validate in silico predictions. A case study applying retrosynthesis is presented as an experimental verification of the output from Retropath, the first complete automated computational pipeline applicable to metabolic engineering. Applying this CAD pipeline, together with genetic reassembly and optimization of culture conditions led to improved production of the plant flavonoid pinocembrin. Coupling CAD tools with advanced genetic engineering strategies and bioprocess optimization is crucial for enhanced product yields and will be of great value for the development of non-natural products through sustainable biotechnological processes.