Psilocybin

In vivo production of psilocybin in E. coli

This synthesis article (2019) explains the gram-scale production of psilocybin in E. coli (a prokaryotic host), which is a 32-fold improvement over earlier techniques. However, it's still much more expensive than growing it in mushrooms/truffles.

Authors

  • Adams, A. M.
  • Brinton, J. D.
  • Enacopol, A. L.

Published

Metabolic Engineering
individual Study

Abstract

Psilocybin, the prodrug of the psychoactive molecule psilocin, has demonstrated promising results in clinical trials for the treatment of addiction, depression, and post-traumatic stress disorder. The development of a psilocybin production platform in a highly engineerable microbe could lead to rapid advances towards the bioproduction of psilocybin for use in ongoing clinical trials. Here, we present the development of a modular biosynthetic production platform in the model microbe, Escherichia coli. Efforts to optimize and improve pathway performance using multiple genetic optimization techniques were evaluated, resulting in a 32-fold improvement in psilocybin titer. Further enhancements to this genetically superior strain were achieved through fermentation optimization, ultimately resulting in a fed-batch fermentation study, with a production titer of 1.16 g/L of psilocybin. This is the highest psilocybin titer achieved to date from a recombinant organism and a significant step towards demonstrating the feasibility of industrial production of biologically-derived psilocybin.

Unlocked with Blossom Pro

Research Summary of 'In vivo production of psilocybin in E. coli'

Introduction

Adams and colleagues situate their work against renewed clinical interest in psilocybin, a prodrug of psilocin that has shown therapeutic promise in trials for conditions such as depression, anxiety in terminal cancer and post‑traumatic stress disorder. Chemical synthesis of psilocybin is possible but remains multi‑step and costly; recent work by other groups identified the fungal biosynthetic genes and demonstrated in vitro and eukaryotic in vivo production, but reported titres from recombinant fungal hosts remained modest. The authors frame a gap around the need for a more easily engineered, cost‑effective microbial production platform that could support larger scale and potentially industrial manufacture for clinical use. This study aimed to heterologously express the key psilocybin biosynthetic genes from Psilocybe cubensis in the prokaryotic host Escherichia coli and to optimise production through parallel genetic library strategies and fermentation parameter tuning. The objective was to demonstrate proof‑of‑principle in E. coli, identify a genetically superior strain, and scale production in a fed‑batch bioreactor to assess whether a bacterial host could reach competitive titres compared with previously reported fungal systems.

Methods

The investigators used Escherichia coli BL21 star (DE3) as the production host and E. coli DH5α for plasmid propagation. Heterologous genes encoding L‑tryptophan decarboxylase (PsiD), a kinase (PsiK), and an S‑adenosyl‑L‑methionine (SAM)‑dependent N‑methyltransferase (PsiM) were obtained from Psilocybe cubensis and cloned into a set of modified ePathBrick/ePathOptimize vectors. Three parallel genetic optimisation strategies were applied: (1) a defined three‑level copy‑number library in which each gene was placed on low, medium or high copy plasmids (27 combinations), (2) a random basic operon library in which all three genes were expressed from a single high‑copy plasmid with randomized T7 promoter variants, and (3) a pseudooperon library with different mutant promoters upstream of individual genes. Libraries were screened in medium throughput using 2 mL cultures in 48‑well plates under standard screening conditions. Standard screening conditions were AMM (Andrew’s Magic Media) supplemented with serine (1 g/L) and 4‑hydroxyindole (350 mg/L), incubation at 37 °C, induction with IPTG (1 mM) four hours after inoculation, and sampling at 24 hours for HPLC analysis. The defined copy‑number library intentionally maintained the same antibiotic burden across strains to control selection pressure. The basic operon library used five mutant T7 promoters (G6, H9, H10, C4, Consensus) generated by standard ePathOptimize methods. For fermentation optimisation, the study varied induction timing (1–6 hours post‑inoculation), IPTG concentration (0.1–1.0 mM), growth/production temperatures (evaluated 30–42 °C with a protocol that initially grew at 37 °C then shifted), media and carbon source identity, and targeted supplements (4‑hydroxyindole, serine, methionine) across several concentrations. Scale‑up was performed as a fed‑batch bioreactor study in which 4‑hydroxyindole was supplied initially at a low concentration and then continuously fed via a syringe pump; feed rates were adjusted using frequent HPLC measurements of pathway intermediates to inform a feedback strategy. The extraction does not fully report detailed bioreactor parameters (agitation, dissolved oxygen set points, pH control and exact feed rates are said to be in supplementary figures). Analytical methods comprised HPLC with diode array detection and refractive index detection for quantification, and high‑resolution LC‑MS/LC‑MS‑MS on an Orbitrap instrument for compound identification. Authentic standards were used where available; for some intermediates expensive standards were substituted with calibration to structurally related compounds. Sample preparation involved simple extraction and direct injection of supernatant for analysis.

Results

Initial heterologous expression of psiD, psiK and psiM under a T7 system in BL21 star (DE3) produced 2.19 ± 0.02 mg/L psilocybin upon IPTG induction. High‑accuracy LC‑MS and MS/MS confirmed the presence of psilocybin and pathway intermediates (4‑hydroxytryptophan, 4‑hydroxytryptamine, norbaeocystin and baeocystin) with mass accuracy better than 5 ppm and expected retention times. Screening the defined copy‑number library yielded only modest improvements over the original ‘‘all‑high’’ construct; the best combination produced 4.0 ± 0.2 mg/L. The basic operon promoter library produced substantially better results: after screening and recloning to exclude plasmid artifacts, operon clone #16 (designated pPsilo16), containing the H10 mutant promoter (a medium‑strength T7 variant), was selected for further study. The top operon mutants showed about a 17‑fold improvement in titre relative to the best clones from the copy‑number library. Further optimisation of fermentation conditions for pPsilo16 raised small‑scale production to 139 ± 2.7 mg/L psilocybin under the identified optimal screening conditions. Key process findings were: maximal production with induction around 3–4 hours post‑inoculation and low sensitivity to induction timing; improved performance with higher IPTG (0.5–1.0 mM versus 0.1 mM); a clear preference for an isothermal 37 °C process rather than lower or higher temperatures; strong sensitivity to media composition and carbon source (rich undefined media such as LB produced low psilocybin and an insoluble coloured product); toxicity and growth inhibition at high concentrations of 4‑hydroxyindole; minimal effect of serine supplementation alone; and a significant enhancement of psilocybin titre when methionine was added in the presence of >350 mg/L 4‑hydroxyindole (p < 0.05). The pseudooperon library mostly underperformed: roughly 95% of mutants showed little or no psilocybin production. The authors report some false positives during screening, which were addressed by recloning and rescreening. Overexpression of native TrpAB was attempted but native levels were judged sufficient, as supported by consistent buildup of 4‑hydroxytryptophan in many fermentations. No significant intracellular accumulation of target metabolites was detected and transmembrane transport was assumed passive but not experimentally investigated. Applying the fermentation knowledge in a fed‑batch bioreactor with an HPLC‑informed feed strategy for 4‑hydroxyindole produced a final psilocybin titre of 1.16 g/L (reported as 1.16 g/L elsewhere in the text). Maximum and final molar yields from the supplied 4‑hydroxyindole substrate were reported as 0.60 and 0.38 mol/mol, respectively. The authors state this is the highest psilocybin titre reported from a recombinant host to date.

Discussion

Adams and colleagues interpret their findings as proof that a prokaryotic host can produce psilocybin at titres that rival and, upon scale‑up, exceed those reported for recombinant fungal systems. They attribute the large gains to the combined strategy of parallel genetic library screening to find a transcriptionally favourable construct (pPsilo16) and iterative fermentation optimisation, culminating in a fed‑batch process that used HPLC feedback to manage feeding of the toxic substrate 4‑hydroxyindole. The authors place their results relative to prior work by noting that small batch titres were similar to those previously obtained in Aspergillus nidulans, but that the fed‑batch scale‑up achieved roughly a 10‑fold enhancement over the published fungal host titres. They highlight practical insights from the genetic screens: the defined copy‑number approach was laborious and yielded limited benefit (likely due to metabolic burden of multiple plasmids), the pseudooperon library largely failed for unknown reasons, and the simple basic‑operon design with a medium‑strength promoter gave the best outcome, emphasising the current limitations of predictive pathway design and the need for screening. Limitations acknowledged in the paper include the reliance on chemical supplementation (4‑hydroxyindole, serine and methionine) rather than de novo biosynthesis from simple carbon sources; the assumption of passive transport across the membrane without direct verification; and several procedural or numerical details relegated to supplementary figures (for example precise bioreactor feed rates and some fermentation parameters are not fully reported in the extracted text). The authors also note the occurrence of false positives during library screening and the unexplained poor performance of most pseudooperon variants. In terms of implications, the study suggests that further metabolic engineering to enable de novo synthesis (removing the need for substrate supplementation) and to improve SAM availability could make biologically produced psilocybin more competitive with chemical synthesis. The authors indicate that such engineering steps are feasible but beyond the scope of the present study, and propose that their platform and the fed‑batch, HPLC‑informed feeding approach provide a foundation for future strain and process development aimed at industrial production.

View full paper sections

INTRODUCTION/BACKGROUND

Psilocybin (4-phosphoryloxy-N,N-dimethyltryptamine) has gained attention in pharmaceutical markets as a result of recent clinical studies. The efficacy of psilocybin has been demonstrated for the treatment of anxiety in terminal cancer patientsand alleviating the symptoms of post-traumatic stress disorder. Most recently, the United States Food and Drug Administration (U.S. FDA) has approved the first Phase IIb clinical trialfor the use of psilocybin as a treatment for depression that is not well controlled with current available interventions such as antidepressants and cognitive behavioral therapies. Psilocybin was first purified from the Psilocybe mexicana mushroom by the Swiss chemist,. The first reports of the complete chemical synthesis of psilocybin were published in 1959; however, large-scale synthesis methods were not developed until the early 2000's by Shirota and colleagues at the National Institute of Sciences in Tokyo. Despite significant improvements over early synthetic routes, current methods remain tedious and costly, involving numerous intermediate separation and purification steps resulting in an overall yield of 49% from 4-hydroxyindole, incurring an estimated cost of $2 USD per milligram for pharmaceutical-grade psilocybin. Many scientists are interested in psilocybin because of its biosynthetic precursorsnorbaeocystin and baeocystin). These compounds have structural similarity to the neurotransmitter serotoninand sparked the interest of researchers who were curious to understand the mechanism behind their hallucinogenic properties. After being named a Schedule I compound in the US with implementation of the Controlled Substance Act of 1970, research efforts involving psilocybin were abandoned for other less regulated bioactive molecules; however, experts in the field have suggested a reclassification to schedule IV would be appropriate if a psilocybin-containing medicine were to be approved in the future. Carhart-Harris and colleagues are conducting ongoing clinical trials with psilocybin as a medication for individuals struggling with treatment-resistant depression. Treatment-resistant forms of depression represent about 20% of the 350 million cases of depression worldwide. During the clinical trials, the medication is administered in a controlled environment, with a team of psychiatrists and psychologists guiding the patients through the psilocybin experience. After a single treatment, 42% of the patients were considered to be in remission, which was classified as having reported improved mood and reduced depressive symptoms for up to three months. These promising results have motivated researchers to investigate biosynthetic routes for the production of psilocybin. The Hoffmeister group recently identified the enzymatic pathway from Psilocybe cubensis (a native fungal psilocybin producer) and demonstrated in vitro synthesis of psilocybin from 4-hydroxytryptophan. In a subsequent work, they utilized the promiscuity of the tryptophan synthase from P. cubensis to make psilocybin from 4hydroxyindole and serine. Additionally, the Hoffmeister group produced psilocybin in vivo using an eukaryotic fungal host, Aspergillus nidulans, at titers reported near 100 mg/L. These breakthroughs in the biosynthesis of psilocybin have resulted in a newfound interest in understanding how to efficiently synthesize the promising and pharmaceutically-relevant molecule for the long-term treatment of multiple chronic disease that affect millions. These new advancements fueled our interest in developing a more cost-effective and easily manipulated host for the biosynthetic production of psilocybin. Utilizing the gene sequences recently identified by the Hoffmeister group from P. cubensis encoding an L-tryptophan decarboxylase (PsiD), a kinase (PsiK), and an S-adensoyl-L-methionine (SAM)-dependent Nmethyltransferase (PsiM), together with the promiscuity of the native Escherichia coli tryptophan synthase (TrpAB), the biosynthesis pathway capable of psilocybin production from 4-hydroxyindole, was expressed in the prokaryotic, model organism E. coli BL21 star TM (DE3). For this proof-of-principle study, E. coli was chosen for the vast amount of metabolic knowledge and genetic tools available to assist in the engineering effort, unlike the limited set for Aspergillus nidulans. Psilocybin production was optimized through the application of a series of 3 parallel genetic optimization methods including: (1) a defined three-level copy number library, (2) a random 5-member operon library, and (3) a random 125-member pseudooperon library. After transcriptional optimization methods were employed, the best strain, pPsilo16, underwent a thorough optimization of fermentation conditions, resulting in the production of ~139 ± 2.7 mg/L of psilocybin from 4-hydroxyindole. We conclude our work with a fed-batch bioreactor scale-up study, which resulted in the production of ~1160 mg/L of psilocybin, the highest titer reported to date from a recombinant host.

PSILOCYBIN PRODUCTION IN E. COLI

The previously confirmed psilocybin production genes (psiD, psiK, and psiM) from P. cubensis were heterologously expressed in E. coli using the strong T7 promoter system. Induction with IPTG allowed for the production of 2.19 ± 0.02 mg/L psilocybin. To confirm compound identities, culture media from our psilocybin production host was subjected to liquid chromatography-mass spectroscopy analysis on a Thermo Orbitrap XL LC-MS system. Psilocybin, as well as all precursor and intermediate compounds in the biosynthetic pathway, were identified with better than 5 ppm mass accuracy as shown in Fig S3 . The sample was then subjected to additional MS/MS fragmentation analysis to further support structural identification of all indole derived intermediates and final products. In each case, fragmentation products for the deamination, dephosphorylation (if applicable), and loss of both functional groups were observed, confirming the identification of psilocybin, and its intermediates: 4hydroxytryptophan, 4-hydroxytryptamine, norbaeocystin, and baeocystin, with better than 5 ppm mass accuracy (Fig). Additionally, expected retention times and order of elution were consistent with previously published efforts. The overexpression of the native tryptophan synthase (TrpAB) was also performed in an attempt to push flux through the heterologous production pathway (data not shown). The native expression level was determined to be sufficient to maintain the necessary pathway flux, as supported by the buildup of 4-hydroxytryptophan in nearly all fermentation studies performed (Fig.).

DEFINED COPY NUMBER LIBRARY

A defined 27-member copy number library consisting of the 3 heterologous biosynthesis genes (psiD, psiK, and psiM) each expressed on 3 different copy number plasmids was constructed and screened in 48-well plates as shown in Fig 2a . Each member of the library contained each of the three genes spread across a low (pACM4-SDM2x), medium (pCDM4-SDM2x), or high (pETM6-SDM2x) copy number plasmid (Fig). The functional expression levelsfor this set of plasmid expression vectors has been previously reported in the literature. This library screen realized minor improvements over the original All-High construct (2.19 ± 0.02 mg/L), where final titers of 4.0 ± 0.2 mg/L were achieved with the combination of psiK expressed from the pETM6-SDM2x vector, psiD expressed from the pCDM4-SDM2x vector, and psiM expressed from the pACM4-SDM2x vector in the BL21 star TM (DE3) expression host. Operon library clones #13 and #15 (Fig.) demonstrated a large reduction in product titer and were identified as false positives in the original screen. Operon clone #16 (pPsilo16, purple) was selected for further study. All combinations were screened in 48-well plates under standard screening conditions and quantified using HPLC analysis. Error bars represent ±1 standard deviation from the mean of replicate samples. *Psilocybin not detected.

PSEUDOOPERON LIBRARY

The pseudooperon library constructs have a different mutant promoter in front of each of the

BASIC OPERON LIBRARY

In the operon configuration, the three-gene pathway was expressed from a single high-copy plasmid under the control of a single promoter and terminator where each gene has an identical ribosome binding site (RBS) (Fig). The promoter sequence was randomized to one of five mutant T7 promoters (G6, H9, H10, C4, Consensus) using the ePathOptimize approach (J. Andrew, resulting in a library that contains 5 potential promoter combinations. After screening nearly 50x the library size, the top 10 variants were selected for further screening. These top variants were re-cloned into an empty plasmid backbone and transformed to eliminate the possibility of spurious plasmid or strain mutations (Fig). Mutant #16 (pPsilo16) was selected for further investigation due to its top production and high reproducibility across multiple fermentations. The sequencing results revealed that pPsilo16 contains the H10 mutant promoter which has been previously characterized as a medium strength promoter, with between 40% and 70% of the effective expression strength of the consensus T7 sequence (J. Andrew. The top mutants from the basic operon screen show a 17-fold improvement in titer over the best performing mutants from the defined copy number library study.

FERMENTATION CONDITIONS OPTIMIZATION

After identifying pPsilo16 as the best strain with respect to highest psilocybin production, low buildup of intermediate products, and consistent reproducibility, the strain underwent a series of optimization experiments to determine the best fermentation conditions for the production of psilocybin. All genetic optimization experiments were conducted under standard conditions (as described in the methods) determined from initial screening in our earlier proofof-principle studies (data not shown). Many studies in the metabolic engineering literature have demonstrated high sensitivity to variations in induction point for pathways controlled by the T7lac inducible promoter. Additionally, induction timing can have a large impact on overall cell growth and can lead to difficulties achieving reproducible production upon scale-up (J Andrew. Upon evaluation of induction sensitivity for pPsilo16, we found that the cells demonstrate low sensitivity to induction point, with the maximum production achieved with induction 3 to 4 hours post inoculation (Fig). No psilocybin production was observed in the non-induced controls (data not shown). Next, optimization of base media selection, carbon source identity, and inducer concentration was evaluated. Since these variables can affect cellular growth rate and corresponding optimal induction points, each of these variables were evaluated across a range of induction points from 1 to 6 hours (Fig S8). As demonstrated in Fig 3bproduction was very sensitive to both media and carbon source selection (p < 0.05). When production was attempted in a rich undefined media such as LB, a dark colored insoluble product was observed along with low psilocybin production. Similarly, low production was also observed when grown on glycerol, however no colored products were observed. pPsilo16 demonstrated moderate sensitivity to IPTG concentration, with higher final concentrations of 0.5 and 1.0 mM outperforming 0.1 mM over a range of induction time conditions (p < 0.05) (Fig). This trend is likely influenced by the initial library screening, which was performed at 1.0 mM IPTG. Production temperatures of 30, 37, 40, and 42 ˚C were also evaluated for their effect on psilocybin production. In an attempt to minimize the effect on changing optimal induction points, all fermentations were started at 37 ˚C through the growth phase of the fermentation before being shifted to the production temperature 1 hour prior to induction. A significant preference (p < 0.05) was seen for maintaining an isothermal fermentation temperature of 37 ˚C throughout both growth and production phases (Fig). The fermentation optimization was completed by evaluating the effects of the targeted media supplements: 4-hydroxyindole, serine, and methionine. Each media supplement was provided at high, medium, and low levels: 4-hydroxyindole (150, 350, and 500 mg/L), serine and methionine (0, 1, and 5 g/L). At high concentrations of 4-hydroxyindole, the cells demonstrated noticeable growth decline due to presumed cellular toxicity leading to reduced productivity. Serine addition showed minimal effects on psilocybin production, however, the addition of methionine in the presence of greater than 350 mg/L of 4-hydroxyindole resulted in a significant enhancement of psilocybin titer (p < 0.05). Under the identified optimal screening conditions, psilocybin was produced at 139 ± 2.7 mg/L, which represents a 63-fold improvement through the synergistic efforts of genetic and fermentation optimization.

SCALE-UP STUDY

After identification of optimized production conditions for pPsilo16, a fed-batch scale up study was completed as described in the methods. This study resulted in the production of 1.16 In the initial proof-of-principle study, the production of psilocybin and all pathway intermediates were confirmed through the use of high mass accuracy LC-MS (Fig). HPLC analysis of fermentation broth from strains containing incomplete pathways (i.e. psiDM and psiDK) was consistent with the conclusions of previous studiesaimed at identifying the order of specific biosynthetic steps in the synthesis pathway (data not shown). Multiple genetic optimization methods were utilized in parallel to identify a genetically superior mutant. Starting with the copy number-based approach, we constructed a 27-member library of 3 pathway genes, each at 3 discrete copy numbers (Fig). Of the three genetic optimization screens presented, this method was the most tedious to construct, requiring each plasmid to be independently cloned and verified prior to screening. This defined library approach also yielded the lowest product titers with the best mutants demonstrating small but statistically significant (p<0.05) improvements over the All-High initial construct. Although a similar modular optimization approach has previously proven successful for a range of products including taxadiene, free fatty acids, and others. We suspect the limited titer improvement from this approach is likely due to the increased metabolic burden associated with selection for and propagation of three independent plasmids. Subsequent screening of two independent single-plasmid transcriptionally-varied promoter libraries with pathway genes in basic operon (Fig) and pseudooperon (Fig) configuration yielded considerably improved results over the initial copy number library. In each case, the library was screened using a medium-throughput HPLC-based screen. Each of these transcriptionally varied libraries were constructed using the high copy pETM6 plasmid vector. This enabled a wide range of expression levels to be screened, resulting in greater coverage of the psilocybin transcriptional landscape. The pseudooperon library screen demonstrated that a large majority of mutants (~95%) showed low or no psilocybin production. The reason for this widespread underperformance is unknown; however, it does motivate the use of random libraries coupled with variant screening for the identification of genetically superior mutants as the current predictive power of a priori pathway design is still lacking for most applications. Surprisingly, the simplistic basic-operon pathway design yielded the highest titer psilocybin production in this study. This coupled with the smallest library size of only 5 mutants, enabled rapid screening of several times the theoretical library size, resulting in high confidence of complete coverage of the full transcriptional landscape. Upon recloning and rescreening the top mutants from the operon library screen, several false positives were identified as shown in Fig. 2d . The source of error for these false positive mutants was not investigated as the false positive rate was at an acceptable level for the study design. Additional increases in titer and yield were achieved through careful optimization of fermentation conditions (Fig). The genetically superior strain, pPsilo16, demonstrated low sensitivity to induction timing as compared to that of other amino acid derived high-value products (Ahmadi and Pfeifer, 2016; J. Andrew; however, this could also be due to the supplementation of both 4-hydroxyindole and serine to the fermentation media, reducing the requirement for high flux through amino acid metabolism. Therefore, all additional fermentation optimization experiments were performed under a range of induction times. Little variation from the induction optimum of 4 hours post inoculation was observed, strengthening the observation of reduced sensitivity to induction timing. The psilocybin production host demonstrated high sensitivity to media composition, carbon source identity, fermentation temperature, and inducer concentration (Fig). In each case, this preferred level was similar to that of the standard screening conditions. This is likely not a coincidence, as some basic initial screening was performed to identify conditions under which our proof-of-principle strain best performed. Furthermore, the initial genetic screening studies were performed under standard screening conditions, which also self-selects for mutants with top performance under the test conditions. The largest gains in the fermentation optimization aspect of this study were achieved through the media supplementation studies. In this study, the concentrations of 4hydroxyindole, serine, and methionine were varied. These supplements were selected specifically for their direct effect on the psilocybin production pathway (Fig.). 4hydroxyindole and serine are condensed by TrpAB in the first dedicated step of the pathway to form the intermediate 4-hydroxytryptophan. Although E. coli can produce serine and indole naturally, it lacks the ability to express the P450 hydroxylase that oxidizes indole into 4hydroxyindole. Additionally, with the high fluxes through our engineered pathway, we hypothesized that the cellular supply of serine would be quickly depleted, requiring additional supplementation to not limit pathway flux. Finally, methionine was supplemented to enhance intercellular pools of the activated methyl donor, SAM. The final two biosynthetic steps are both catalyzed by the SAM-dependent methyltransferase, PsiM. Previous studies with SAMdependent methylations in E. coli have documented SAM-limited flux to final products. In each of these cases, the need for the media supplementation can be addressed using metabolic engineering approaches that are beyond the scope of this study, but certainly applicable in future studies to improve the efficiency of psilocybin production in vivo. The information gained from the genetic and fermentation optimization studies was applied in a scale-up study for the production of psilocybin in a fed-batch bioreactor. In this study, many of the optimization parameters such as temperature, inducer concentration, and induction timing were applied as previously optimized. Information from the supplement addition studies was used but applied with modification from the 2 mL batch studies. In the fedbatch studies, both serine and methionine were supplemented at the high level of 5 g/L to account for higher cellular demand due to enhanced cell growth. Furthermore, in the small-scale studies a growth deficit was observed at higher concentrations of 4-hydroxyindole and 4hydroxytryptophan. To counter this, a low amount of 4-hydroxyindole (150 mg/L) was added initially to the media, while a low-flow syringe pump, containing a 40 mg/mL 4-hydroxyindole solution, was connected for slow external supplementation. To determine the optimal feed rate, the pathway flux through the bottleneck point, PsiD, was estimated through frequent HPLC analysis of the fermentation broth. As 4-hydroxytryptophan titers fell, the flux of 4hydroxyindole was increased to meet the high flux demand, and vise-versa. This strategy resulted in an oscillatory concentration profile for 4-hydroxytryptophan and maintained all intermediates at low levels, enabling robust and extended growth and psilocybin production (Fig). Additional details on substrate feed rate, cumulative yield, and bioreactor parameters (agitation, DO, pH, and temperature) can be found in Supplementary Figure.

ANALYSIS

In small batch fermentation studies, the work presented above resulted in a similar titer of psilocybin to that presented previously in the A. nidulans host.This indicates that both bacterial and fungal hosts show potential as production platforms for this important chemical. However, upon scale-up to a fed batch reactor our bacterial host demonstrated greatly enhanced psilocybin production resulting in a 10-fold enhancement over previously published results. In conclusion, we have presented the first example of psilocybin production in a prokaryotic organism and the highest psilocybin titer to date from a recombinant host from any kingdom. This was accomplished through the combination of genetic and fermentation optimization in small scale, coupled with a scaled-up fed-batch study utilizing a unique HPLC informed substrate feeding strategy. The fed-batch study resulted in a psilocybin titer of 1.16 g/L with maximum and final molar yields from the 4-hydroxyindole substrate of 0.60 and 0.38 mol/mol, respectively. Future efforts to enhance the strain background to enable de novo production of psilocybin, without the need for chemical supplementation, could enable a biologically synthesized psilocybin pharmaceutical product to compete with chemical synthesis strategies for applications in medicine, psychology, and neurobiology.

METHODS

Bacterial strains, vectors, and media E. coli DH5α was used to propagate all plasmids, while BL21 star TM (DE3) was used as the host for all chemical production experiments. Plasmid transformations were completed using standard electro and chemical competency protocols as specified. Unless noted otherwise, Andrew's Magic Media (AMM)was used for both overnight growth and production media, while Luria Broth (LB) was used for plasmid propagation during cloning. The antibiotics ampicillin (80 µg/mL), chloramphenicol (25 µg/mL), and streptomycin (50 µg/mL) were added at their respective concentrations to the culture media when using pETM6, pACM4, and pCDM4-derived vectors, respectively. A description of all plasmids and strains used in this study can be found in Table. The exogenous pathway genes encoding the enzymes PsiD, PsiK, and PsiM contained on plasmids pJF24, pJF23, and pFB13, respectively, were a generous gift from the Hoffmeister group.

PLASMID CONSTRUCTION

The original ePathBrick expression vectors, #4, #5, and #6 (Table) were modified through two rounds of site directed mutagenesis with primers 1 through 4 (Table) to result in the corresponding 'SDM2x' series of vectors: #7, #8, and #9 (Table). This mutagenesis was performed to swap the positions of the isocaudomer restriction enzyme pair XmaJI/XbaI in the vector. This change allows for the monocistronic and pseudooperon pathway configurations to be constructed more cost efficiently by avoiding the use of the costly XmaJI restriction enzyme. This series of vectors was then used to construct the vectors used in the defined copy number library study #10 -#27 (Table). Plasmids #1 -#3 containing psiD, psiK, and psiM, respectively, were restriction enzyme digested with NdeI and HindIII, gel extracted, and ligated into the pETM6-SDM2x (#7, Table) plasmid backbone, resulting in plasmids #10, #11, and #12 (Table). All multigene expression plasmids were constructed in pseudooperon configuration using a modified version of the previously published ePathBrick methods as described above, while all transcriptional libraries were constructed using standard ePathOptimize methods (J. Andrew.

STANDARD SCREENING CONDITIONS

Standard screening was performed in 2 mL working volume cultures in 48-well plates at 37 ˚C. AMM supplemented with serine (1 g/L), 4-hydroxyindole (350 mg/L), and appropriate antibiotics were used unless otherwise noted. Overnight cultures were grown from either an agar plate or freezer stock culture in AMM with appropriate antibiotics and supplements for 14-16 hours in a shaking 37 ˚C incubator. Induction with 1 mM isopropyl β-D-1thiogalactopyranoside (IPTG) occurred four hours after inoculation, unless otherwise noted. Cultures were then sampled 24 hours post inoculation and subjected to HPLC analysis as described in analytical methods below.

LIBRARY CONSTRUCTION

The defined copy number library was constructed using plasmid #7 (High), #8 (Medium), and #9 (Low). The pathway genes were modulated in either the high, medium, and low copy number vectors, as shown in Figure. The BL21 star TM (DE3) production host was transformed with the appropriate plasmids such that each strain had all three vectors, even if some were empty, to enable the same antibiotic resistance burden to be present in all defined library members (Fig). In the cases where multiple genes were present at a single expression level, the plasmids were constructed in pseudooperon configuration as described above. Random promoter libraries were assembled using standard ePathOptimize methods with the five original mutant T7 promoters: G6, H9, H10, C4, and consensus. Random libraries were built in pseudooperon (Fig) and basic operon (Fig) forms, maintaining a sufficient number of colonies at each cloning step as to not limit library size. was defined under the conditions of 37 ˚C, pH 7.0, 250 rpm agitation, and 3 lpm of standard air. The zero-oxygen set point was achieved by a nitrogen gas flush. Samples were collected periodically for measurement of OD 600 and metabolite analysis. The bioreactor was induced with analyzed via HPLC on an approximate 45-minute delay and were used as feedback into the feeding strategy described above.

ANALYTICAL METHODS

Samples were prepared by adding an equal volume of 100% ethanol or 100% deionized water and fermentation broth, vortexed briefly, and then centrifuged at 12000 x g for 10 minutes. 2 µL of the resulting supernatant was then injected for HPLC or LC-MS analysis. Analysis was performed on a Thermo Scientific Ultimate 3000 High-Performance Liquid Chromatography (HPLC) system equipped with Diode Array Detector (DAD) and Refractive Index Detector (RID). Authentic standards were purchased for glucose (Sigma), psilocybin (Cerilliant), and 4hydroxyindole (BioSynth). Standards for baeocystin, norbaeocystin, 4-hydroxytryptamine, and 4-hydroxytryptophan were quantified using a standard for a similar analog due to limited commercial availability and extremely high cost, approx. $2000 USD for 1 mg of the authentic standard. Baeocystin and norbaeocystin were quantified on the psilocybin standard curve, while 4-hydroxytryptamine and 4-hydroxytryptophan were quantified on the standard curves of 5hydroxytryptamine (Alfa Aesar) and 5-hydroxytryptophan (Alfa Aesar), respectively. No significant intracellular accumulation of target metabolites was observed upon analysis with and without cell lysis. Transport across the cell membrane was assumed to be passive, however, specific investigation into this phenomenon was not undertaken for this work. Glucose analysis was performed using an Aminex HPX-87H column maintained at 30 ˚C followed by a refractive index detector (RID) held at 35 ˚C. The mobile phase was 5 mM H 2 SO 4 in water at a flow rate of 0.6 mL/min. Glucose was quantified using a standard curve with a retention time of 8.8 min. UV absorbance at 280 nm was used to quantify all aromatic compounds. Analysis was performed using an Agilent ZORBAX Eclipse XDB-C18 analytical column (3.0 mm x 250 mm, 5 µm) with mobile phases of acetonitrile (A) and water (B) both containing 0.1% formic acid at a flow rate of 1 mL/min: 0 min, 5% A; 0.43 min, 5% A; 5.15 min, 19% A; 6.44 min, 100 % A; 7.73 min 100% A; 7.73 min, 5% A; 9.87 min, 5% A. This method resulted in the following observed retention times: psilocybin (2.2 min), baeocystin (1.7 min), norbaeocystin (1.9 min), 4hydroxytryptamine (3.4 min), 4-hydroxytryptophan (3.6 min), and 4-hydroxyindole (6.6 min). High Resolution Liquid Chromatography Mass Spectrometry (LC-MS) and Mass Spectrometry-Mass Spectrometry (LC-MS/MS) data were measured on a Thermo Scientific LTQ Orbitrap XL mass spectrometer equipped with an Ion Max ESI source using the same mobile phases and column described above. The flow rate was adjusted to 0.250 mL/min resulting in a method with the following gradient: 0 min, 5% A; 1 min, 5% A; 24 min, 19% A; 30 min, 100 % A; 36 min 100% A; 36 min, 5% A; 46 min, 5% A. This method resulted in the following observed retention times: psilocybin (8.7 min), baeocystin (7.6 min), norbaeocystin (6.4 min), 4-hydroxytryptamine

Study Details

Your Library