Mining Enzyme Diversity of Transcriptome Libraries through DNA

Jul 21, 2016 - Phylogenetic trees were generated from multiple sources of sequence data to identify a library of candidate enzymes that were purchased...
0 downloads 3 Views 2MB Size
Subscriber access provided by La Trobe University Library

Article

Mining enzyme diversity of transcriptome libraries through DNA synthesis for benzylisoquinoline alkaloid pathway optimization in yeast Lauren Narcross, Leanne Bourgeois, Elena Fossati, Euan Burton, and Vincent J. J. Martin ACS Synth. Biol., Just Accepted Manuscript • DOI: 10.1021/acssynbio.6b00119 • Publication Date (Web): 21 Jul 2016 Downloaded from http://pubs.acs.org on July 22, 2016

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

ACS Synthetic Biology is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Synthetic Biology

1

Mining enzyme diversity of transcriptome libraries through DNA synthesis for

2

benzylisoquinoline alkaloid pathway optimization in yeast

3 4

Lauren Narcross1,2, Leanne Bourgeois1,2, Elena Fossati, Euan Burton1,2, Vincent J.J. Martin1,2*

5 6

1

Department of Biology, Concordia University, Montréal, Québec, Canada, H4B 1R6

7 2Centre for Structural and Functional Genomics, Concordia University, Montréal, Québec, Canada, H4B 1R6

8 9 10 11 12 13 14 15 16 17 18 19 20 21 22

* Corresponding author:

23

E-mail: [email protected]

ACS Paragon Plus Environment

ACS Synthetic Biology

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

1

Page 2 of 38

Abstract

2

The ever-increasing quantity of data deposited to GenBank is a valuable resource for

3

mining new enzyme activities. Falling costs of DNA synthesis enables metabolic engineers to

4

take advantage of this resource for identifying superior or novel enzymes for pathway

5

optimization.

6

dihydrosanguinarine in yeast from norlaudanosoline at a molar conversion of 1.5%. Molar

7

conversion could be improved by reduction of the side-product N-methylcheilanthifoline, a key

8

bottleneck in dihydrosanguinarine biosynthesis. Two pathway enzymes, an N-methyltransferase

9

and a cytochrome P450 of the CYP719A subfamily, were implicated in the synthesis of the side-

10

product. Here, we conducted an extensive screen to identify enzyme homologs whose co-

11

expression reduces side-product synthesis. Phylogenetic trees were generated from multiple

12

sources of sequence data to identify a library of candidate enzymes that were purchased codon-

13

optimized and pre-cloned into expression vectors designed to facilitate high-throughput analysis

14

of gene expression as well as activity assay. Simple in vivo assays were sufficient to guide the

15

selection of superior enzyme homologs that ablated the synthesis of the side-product, and

16

improved molar conversion of norlaudanosoline to dihydrosanguinarine to 10%.

Previously,

we

reported

synthesis

of

the

benzylisoquinoline

alkaloid

17 18

Keywords

19

Synthetic DNA; transcriptome mining; benzylisoquinoline alkaloids; dihydrosanguinarine;

20

Saccharomyces cerevisiae; pathway optimization

2

ACS Paragon Plus Environment

Page 3 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

1

ACS Synthetic Biology

Introduction

2

A vast amount of genome and transcriptome data is deposited in publicly available

3

resources such as GenBank, which reached a milestone of one trillion base pairs of sequence data

4

in January 2015.1 More targeted databases like the Thousand Plant2 and PhytoMetaSyn Projects3-

5

4

6

evolutionary analysis but have traditionally presented few opportunities for metabolic engineers

7

due to the lack of physical DNA available to them.5 Until recently, the RNA used to generate

8

transcriptome sequence databases was also the source of cDNA used for the targeted

9

amplification of putative ORFs and gene discovery.6-13 With the cost of DNA synthesis falling

10

from $1/bp in 2006 to $0.12/bp in 2014,14 digital sources of DNA sequences are becoming

11

broadly-accessible primary resources of unique enzymes for the purposes of pathway

12

optimization and the identification of novel activities. This information represents an attractively

13

simple alternative to more traditional methods of pathway optimization through protein

14

engineering approaches such as directed evolution15-16 or rational modification.17-18 For example,

15

heterologous synthesis of methyl halides was enabled through the screening of 89 putative

16

methyl halide transferases from metagenomics data deposited to NCBI.

17

enzyme bottleneck in the heterologous synthesis of coumarate was alleviated through the

18

screening of a library of both putative and published enzymes purchased entirely from

19

GenBank.20

provide further sources of sequence information. Such in silico resources are valuable for

19

More recently, an

20

Previously, we reported the reconstitution in Saccharomyces cerevisiae of a 10-gene

21

pathway for the synthesis of the benzylisoquinoline alkaloid (BIA) dihydrosanguinarine – a

22

reduced form of the antimicrobial sanguinarine – from the precursor norlaudanosoline (Figure

23

1A).21 De novo synthesis of benzylisoquinolines in yeast is currently at the microgram/liter level

3

ACS Paragon Plus Environment

ACS Synthetic Biology

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 38

1

due to low precursor titers and poor performance of the enzyme catalyzing the committed step,

2

22-24

3

When supplemented to culture medium, the highest-reported molar conversion of the substrate

4

norlaudanosoline to the key branch point reticuline in yeast is 20%.21,

5

expression of the 7 enzymes necessary for dihydrosanguinarine synthesis from reticuline drops

6

molar conversion to 1.5%21, with the accumulation of the dead-end intermediate N-

7

methylcheilanthifoline contributing to the drop in yield (Figure 1A).21

although the conversion of norlaudanosoline to downstream products can also be inefficient.

24

The additional co-

8

Conversion of the dihydrosanguinarine pathway intermediates scoulerine to stylopine is

9

catalyzed by two cytochrome P450 enzymes of the CYP719A subfamily: cheilanthifoline

10

synthase (CFS) converts scoulerine to cheilanthifoline, and stylopine synthase (SPS) converts

11

cheilanthifoline to stylopine (Figure 1A). However, in in vitro and heterologous in vivo systems,

12

stylopine synthase activity is insufficient and leads to cheilanthifoline accumulation.21, 25 Since

13

cheilanthifoline

14

tetrahydroprotoberberine N-methyltransferase (TNMT), accumulation of cheilanthifoline results

15

in the synthesis of the non-productive intermediate N-methylcheilanthifoline.

is

also

a

substrate

for

the

promiscuous

N-methylating

enzyme

16

Enzyme promiscuity leading to non-productive intermediates is a common problem in

17

heterologous pathway reconstitution.24, 26-30 Many strategies for reducing side-reactions, such as

18

compartmentalization of competing reactions or enzyme engineering for improved specificity,

19

work within the constraints of enzymes currently in use.24,

20

dihydrosanguinarine pathway, we demonstrate the power of mining transcriptome libraries

21

combined with gene synthesis as an effective strategy for pathway engineering. We postulated

22

that either an SPS that is able to outcompete TNMT for cheilanthifoline, or a TNMT with

4

ACS Paragon Plus Environment

31-34

Here, using the

Page 5 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Synthetic Biology

1

narrower substrate selectivity, or both, would prevent N-methylcheilanthifoline synthesis and

2

greatly improve current dihydrosanguinarine yields.

3

In this work, two enzyme libraries, one of TNMTs and one of CYP719s were purchased

4

as codon-optimized synthetic genes and screened individually and in combinations. In assaying

5

these 73 enzymes, a new activity was discovered, which inspired a new route to stylopine

6

synthesis (Figure 1B). Consequently, synthesis of N-methylcheilanthifoline was ablated. The

7

newly engineered dihydrosanguinarine pathway now reaches 10% conversion in yeast cultures

8

supplemented with the precursor norlaudanosoline. The strategy described here is a simple

9

alternative to more rational methods that can be applied to any pathway that requires

10

optimization.

11 12

Results

13

Generation of CYP719 and NMT enzyme libraries

14

Synthesis of the BIA dihydrosanguinarine from norlaudanosoline requires nine enzymatic

15

reactions (Figure 1A). Two of these are catalyzed by N-methyltransferases: conversion of 6-O-

16

methylnorlaudanosoline to 3’-hydroxy-N-methylcoclaurine by coclaurine N-methyltransferase

17

(CNMT), and conversion of stylopine to N-methylstylopine by tetrahydroprotoberberine N-

18

methyltransferase (TNMT). NMTs from BIA-producing plants can accept a variety of BIAs as

19

substrates (Table 1). Nevertheless, NMTs have also been demonstrated to differentiate between

20

BIAs that differ by a single methyl group or methylenedioxy bridge.7 Thus, one approach to

21

reducing N-methylcheilanthifoline synthesis was to identify a TNMT that accepted stylopine but

22

not cheilanthifoline. CNMTs were also considered, as some can N-methylate downstream

23

dihydrosanguinarine pathway intermediates (Table 1). The reverse-translated PhytoMetaSyn

5

ACS Paragon Plus Environment

ACS Synthetic Biology

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

1

transcriptome database was queried using a conserved TNMT/CNMT motif. Putative ORFs

2

selected from the transcriptome database were aligned with published TNMTs/CNMTs using

3

MUSCLE and a phylogenetic tree was generated with the program MEGA6 (Figure 2A).35 The

4

phylogenetic tree served as a guide for the choice of enzyme candidates to be screened. A total of

5

15 published and putative NMTs were purchased.

6

Conversion of the dihydrosanguinarine pathway intermediate scoulerine to stylopine

7

requires the formation of two methylenedioxy bridges, indicated by “A” and “B” in Figure 1B.

8

While theoretically the reactions could occur in either order, it has been experimentally

9

determined in planta that Ring B closure (catalyzed by CFS) occurs before Ring A closure

10

(catalyzed by SPS).36 Both CFS and SPS are cytochrome P450s in the CYP719A subfamily. This

11

subfamily also includes other members that catalyze methylenedioxy bridge formations on other

12

BIAs and other alkaloids (Table 2), and still other methylenedioxy bridge-containing alkaloids

13

have been identified for which the appropriate methylenedioxy bridge-forming enzymes are still

14

unknown. Diversity amongst methylenedioxy-bridge containing alkaloids and CYP719 substrate

15

acceptance profiles suggests that the CYP719 enzyme family is extensive and may include SPS

16

enzyme homologs that are more appropriate for heterologous reconstitution of the

17

dihydrosanguinarine pathway. Traditionally, the naming scheme for CYP719s is based on an

18

identified product (i.e. stylopine synthase and cheilanthifoline synthase). However, this naming

19

scheme becomes untenable when the same enzyme can synthesize multiple products. Here, we

20

refer to CYP719s by the location of methylenedioxy bridge formation: Ring A-closing CYP719s

21

and Ring B-closing CYP719s.

22

Reverse-translated transcriptome data from the PhytoMetaSyn database was queried for a

23

conserved heme-binding cytochrome P450 motif and an N-terminal motif conserved amongst

6

ACS Paragon Plus Environment

Page 6 of 38

Page 7 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Synthetic Biology

1

published CYP719s. Queries were narrowed down using BLASTclust. BLASTclust sorts

2

sequences into groups using the criterion of percent sequence identity, which is convenient for

3

the study of CYPs because CYP families and subfamilies are defined based on this criterion

4

(45% amino acid identity defines a family, 55% amino acid identity defines a subfamily).37

5

Stringency was set to 50% in order to include the CYP719B subfamily, which also has activity

6

on BIAs.38 Putative CYP719s that clustered with published CYP719s were aligned with

7

MUSCLE, and a phylogenetic tree was generated using MEGA6 (Figure 2B).35 Three clades

8

were observed, which were assigned predicted activities based on the co-alignment with

9

characterized CYP719s: Ring A-closing CYP719s, further segregated into two subclades of

10

CYP719s predicted to act on cheilanthifoline (stylopine synthases) or on the BIA

11

tetrahydrocolumbamine (canadine synthases); Ring B-closing CYP719s, further segregated into

12

one subclade of cheilanthifoline synthases and one subclade of CYP719Bs; and CYP719s with

13

unknown activities. A total of 54 characterized and putative CYP719s were purchased for

14

screening.

15 16

Selection of replacement Ring A-closing CYP719s

17

The library of CYP719s included enzymes with characterized activity on relevant BIAs

18

(scoulerine and cheilanthifoline), other BIAs or other alkaloids, as well as enzymes with

19

predicted activity or no predicted activity (Table 2, Figure 2B). Following a qualitative

20

assessment of CYP719 expression through comparison of the fluorescence of CYP719-GFP

21

fusion proteins (see Supporting Results & Discussion), an initial activity screen was performed

22

to validate predicted Ring A- and Ring B-closing activities within the CYP719 library. CYP719s

23

were expressed in yeast and supplemented with the dihydrosanguinarine pathway intermediate

7

ACS Paragon Plus Environment

ACS Synthetic Biology

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 38

1

scoulerine. Scoulerine is a substrate for both Ring A- and Ring B-closing CYP719s, forming

2

nandinine

3

cheilanthifoline have the same mass and similar structures, they are distinguishable by HPLC-

4

MS both through elution time and by MS/MS profile (Figure 3B). Analysis of nandinine and

5

cheilanthifoline synthesis confirmed that predicted Ring A- and Ring B-closing activities were

6

generally accurate with no cases of Ring B-closure where Ring A-closure was predicted, or vice

7

versa (Figures 3C, 3D). In addition, no conversion of scoulerine was detected for predicted

8

CYP719Bs, CYP719s with characterized activity on non-BIAs, or putative CYP719s with no

9

predicted activity. Scoulerine was widely accepted amongst Ring B-closing CYP719s with 10 of

10

12 candidates converting >95% of the scoulerine to cheilanthifoline (Figures 3C, 3D).

11

Conversely, a greater range of nandinine synthesis was observed amongst predicted Ring A-

12

closing CYP719s. Of these, 7 of 18 predicted stylopine synthases and 3 of 16 predicted canadine

13

synthases converted >95% of the scoulerine to nandinine (Figures 3C, 3D). These Ring A-

14

closing CYP719s were considered for further characterization.

and

cheilanthifoline,

respectively

(Figure

3A).

Although

nandinine

and

15

Next, selected Ring A-closing CYP719s were co-expressed with the Ring B-closing

16

CYP719 PsCFS, previously used in the heterologous reconstitution of dihydrosanguinarine

17

synthesis.21, 39 When supplemented with scoulerine, the expected product is stylopine (Figure

18

4A). In each yeast strain, scoulerine was entirely consumed but different proportions of

19

cheilanthifoline, nandinine, and stylopine were observed depending on the co-expressed Ring A-

20

closing CYP719 (Figure 4A). Expression of Ring A-closing CYP719s predicted to be canadine

21

synthases resulted in residual cheilanthifoline. These candidates were not considered for further

22

screening. With one exception, cheilanthifoline was not detected when the Ring A-closing

23

CYP719 was predicted to be a stylopine synthase.

8

ACS Paragon Plus Environment

Page 9 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Synthetic Biology

1

Unlike previously reported combinations of cheilanthifoline and stylopine synthases,

2

nandinine was a product observed in each combination of Ring A-closing CYP719s and

3

PsCFS.21,

4

PsCFS for scoulerine. Accumulation of nandinine also indicated that it is not a preferred

5

substrate of PsCFS. Improved activity of Ring A-closing CYP719s relative to PsCFS is a desired

6

quality, but the generation of a new side product is not. Nandinine accumulation could be

7

avoided by limiting the pool of potential SPS’s to those with activity on cheilanthifoline but not

8

scoulerine. Because many Ring A-closing CYP719s in the library can accept scoulerine, this is

9

not an ideal limitation. Alternatively, if a Ring B-closing CYP719 could be identified that also

10

accepts nandinine as substrate, then nandinine could be re-captured into the main pathway,

11

shifting from a side-product to a pathway intermediate (Figure 1B).

25, 39

Nandinine synthesis resulted from Ring A-closing CYP719s out-competing

12 13

Engineering of a non-natural stylopine synthesis pathway

14

To simplify further assessment of Ring A- and Ring B-closing CYP719 activity, the

15

appropriate substrates were supplied directly to yeast strains. As cheilanthifoline and nandinine

16

were not available commercially, they were generated from scoulerine by incubation with yeast

17

expressing an appropriate CYP719. Supernatant containing cheilanthifoline was used as

18

substrate to test the activity of Ring A-closing CYP719s (Figure 4B), while supernatant

19

containing nandinine was used to test Ring B-closing CYP719s (Figure 4C). As was suggested in

20

co-expression analysis, all selected Ring A-closing CYP719s were able to convert >95% of

21

cheilanthifoline to stylopine. Many Ring B-closing CYP719s had some activity on nandinine,

22

with 2 of 10 candidates (Sdi-1 and Cma-2) converting >95% of supplemented nandinine to

23

stylopine. Acceptance of nandinine by the two Ring B-closing CYP719s enables the non-natural

9

ACS Paragon Plus Environment

ACS Synthetic Biology

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

1

“Ring A first” pathway for stylopine synthesis (Figure 1B). Hence, the six Ring A-closing

2

CYP719s and two Ring B-closing CYP719s that converted >95% of their supplemented BIA to

3

stylopine were selected for combinatorial testing in the presence of TNMT.

Page 10 of 38

4 5

Combinatorial testing of Ring A- and Ring B-closing CYP719s

6

Combinations of Ring A- and Ring B-closing CYP719s were next co-expressed in the

7

presence and absence of TNMT to measure production of downstream N-methylated BIAs

8

(Figure 5). When supplemented with scoulerine, stylopine should be the product in the absence

9

of TNMT and N-methylstylopine should be the product in the presence of TNMT. Any

10

accumulated cheilanthifoline, nandinine, or their N-methylated derivatives, would indicate an

11

undesired combination of CYP719s. Nandinine was observed during co-expression of PsCFS

12

with Ring A-closing CYP719s, but not with either of the two selected Ring B-closing CYP719s.

13

Residual cheilanthifoline (and N-methylcheilanthifoline in the presence of TNMT) was observed

14

in samples expressing 2 of the 6 Ring A-closing CYP719s, which was not expected because

15

these enzymes previously converted >95% of cheilanthifoline to stylopine (Figure 4B). Between

16

experiments, Ring A-closing CYP719s had been placed under the control of a new promoter /

17

terminator pair in order to allow homology-mediated cloning of a double CYP719 gene cassette.

18

The other 4 Ring A-closing CYP719s, when co-expressed with either of the 2 Ring B-closing

19

CYP719s, resulted in >95% conversion of scoulerine to stylopine in the absence of TNMT, and

20

>95% conversion of scoulerine to N-methylstylopine in the presence of TNMT. These

21

combinations were selected for integration into the dihydrosanguinarine pathway.

22 23

Selection of replacement TNMT

10

ACS Paragon Plus Environment

Page 11 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Synthetic Biology

1

Before retrofitting a final dihydrosanguinarine production strain with the newly

2

discovered CYP719s, the NMT library was screened for NMTs with greater substrate preference

3

for stylopine. Yeast strains harboring the NMT library were supplemented with scoulerine

4

(Figure 6A) and stylopine (Figure 6B) in order to assess relevant substrate acceptance profiles. In

5

addition, the activity of NMTs on cheilanthifoline was also assayed through the supplementation

6

of scoulerine to yeast strains co-expressing the NMT library and the Ring B-closing CYP719

7

PsCFS (Figure 6C). Included, as negative controls, were two O-methyltransferases (OMTs) in

8

the dihydrosanguinarine pathway that were not expected to methylate the supplemented BIAs.40

9

In general, if an NMT N-methylated stylopine, it also N-methylated cheilanthifoline and

10

scoulerine (Figure 6). These NMTs all aligned with published TNMTs, but not every enzyme

11

aligning with TNMTs had activity on the BIAs tested here. Amongst CNMTs, PsCNMT is

12

uniquely able to N-methylate scoulerine and cheilanthifoline, but no CNMT was able to N-

13

methylate stylopine. Furthermore, when co-expressed with PsCFS, NMTs with no measurable

14

activity on the BIAs tested appeared to interfere with cheilanthifoline synthesis, as conversion of

15

scoulerine to cheilanthifoline was lower than either empty vector or OMT control strains. Since

16

we did not identify an NMT that methylated stylopine without also methylating scoulerine

17

and/or cheilanthifoline, the strategy of identifying an NMT with a different substrate acceptance

18

profile was not pursued further.

19 20

Retrofitting and testing an optimized dihydrosanguinarine producing strain

21

Through CRISPR-directed homologous recombination, most of the norlaudanosoline-to-

22

dihydrosanguinarine pathway was chromosomally integrated into a single yeast strain. Missing

23

were the CYP719 Ring A- and Ring B-closing enzymes, which were combinatorially co-

11

ACS Paragon Plus Environment

ACS Synthetic Biology

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 38

1

expressed from plasmids. Cultures of the resulting strains were then supplemented with

2

norlaudanosoline at increasing concentrations, and levels of N-methylated intermediates and

3

dihydrosanguinarine were measured. In the absence of CYP719s, N-methylscoulerine was

4

produced as a result of the accumulation of scoulerine (Figure 7A). Normalized N-

5

methylscoulerine levels remained constant across increasing substrate concentrations, which

6

indicated that there were no measurable bottlenecks in the pathway from norlaudanosoline to

7

scoulerine. N-methylscoulerine was not observed when CYP719s were expressed. As expected,

8

BIAs extracted from yeast cultures expressing the original CYP719 combination (PsSPS and

9

PsCFS) showed N-methylcheilanthifoline accumulation at every norlaudanosoline concentration.

10

Further, as norlaudanosoline concentration increased, N-methylcheilanthifoline levels rose

11

relative to normalized values. In contrast, N-methylcheilanthifoline did not accumulate with any

12

combination of the newly selected CYP719s, at any concentration of norlaudanosoline. When the

13

strain co-expressing PsSPS and PsCFS was supplemented with 10 µM norlaudanosoline, 5%

14

was converted to dihydrosanguinarine and sanguinarine (Figure 7B). Multiple combinations of

15

Ring A- and Ring B-closing CYP719s resulted in improved levels of dihydrosanguinarine and

16

sanguinarine, reaching ~10% conversion of norlaudanosoline to dihydrosanguinarine with no

17

intermediate or side product observed. At higher concentrations of norlaudanosoline, the

18

pathway intermediate N-methylstylopine accumulated, increasing relative to normalized

19

norlaudanosoline conversion values.

20 21 22 23

Discussion With the cost of next-generation DNA sequencing less than $1/million base pairs,14 transcriptome databases can be cheaply generated and used for enzyme discovery.6,

12

ACS Paragon Plus Environment

8, 41-43

Page 13 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Synthetic Biology

1

Historically, harnessing the power of transcriptome databases was difficult without access to the

2

physical RNA used to generate the libraries. Now, advances in DNA synthesis technologies44

3

and molecular biology techniques for efficient heterologous gene expression45-48 can accelerate

4

and improve the enzyme discovery process. Here, we demonstrate the power of combined

5

accessibility to transcriptome data and affordable gene synthesis for pathway engineering and

6

optimization. Not only did our synthetic gene library contain multiple enzymes capable of

7

improving the pathway, it also included a novel activity enabling pathway redesign to prevent

8

the synthesis of side products. As the cost of DNA synthesis continues to drop, we foresee this

9

strategy becoming increasingly common for pathway engineering and optimization.

10

The work presented here is an example of the power and limits of predictive search for

11

enzymatic activities from sequence data. Predicted activities for putative ORFs were assigned

12

based on co-alignment of characterized enzymes within clades. For both NMTs and CYP719s,

13

broad activities were accurate (i.e. N-methylase activity or Ring A- vs. Ring B-closure), but

14

enzymes co-aligned within a sub-clade did not necessarily have the same substrate acceptance

15

profiles. This is not an uncommon phenomenon,6, 49-50 which highlights the value in using large

16

libraries of orthologous enzymes to increase the chances of finding an activity of interest.

17

For both published NMTs and CYP719s, our analysis of substrate preferences included

18

more positive hits than previous characterization. For instance, TNMTs have been demonstrated

19

to have much lower activity on scoulerine than stylopine (0-10% of relative activity depending

20

on the homolog),7, 51-53 whereas no difference was observed here (Figures 6A, 6C). Similarly,

21

some CYP719s such as CYP719A1 have been shown to have little to no activity on scoulerine54

22

but here were able to convert >95% of supplemented scoulerine to nandinine. Unlike traditional

23

biochemical assays, supplementation and bioconversion assays often occur over a longer

13

ACS Paragon Plus Environment

ACS Synthetic Biology

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

1

duration (1 hr vs. >16 hrs) and are a measure of reaction progression, not reaction speed.55-56 For

2

the identification of enzymes with novel activities for heterologous pathway engineering, longer

3

incubation times in in vivo conditions may flag candidates that may otherwise be discarded.

Page 14 of 38

4

In this work, all CYP719s were co-expressed with a single cytochrome P450 reductase

5

(CPR) from P. somniferum (PsCPR). This setup was sufficient to identify both Ring A- and Ring

6

B-closing CYP719s from multiple plant species with activity on supplemented BIAs. However,

7

PsCPR may not be an ideal partner for all the CYPs in the library. It has long been recognized

8

that any one CPR cannot support the activity of all CYPs.57 While a CPR from the same species

9

is often used if available,58-60 the existence of multiple CPR in plants complicates the selection

10

process.61 Further, the actual relationship between CYP and CPR is unpredictable, with CPRs

11

from other plants like Arabidopsis thaliana38-39, 62 and even the native yeast CPR22 supporting

12

heterologous activity of some CYPs. Therefore we presume that combinatorial co-expression of

13

the CYP719 library with a putative CPR library from other organisms may improve the activity

14

of some of the CYPs in our library.

15

Enzyme promiscuity presents serious challenges for the reconstitution of heterologous

16

pathways,24, 34, 63 which can compound as a pathway increases in size. For example, scoulerine is

17

a substrate for 4 of 9 enzymatic steps in the dihydrosanguinarine pathway (Tables 1 and 2).

18

Nevertheless, in vivo combinatorial screens identified multiple enzyme combinations in which no

19

scoulerine side-products were found to accumulate. Initially, nandinine synthesis was an

20

undesired activity. Although nandinine acceptance by Ring B-closing CYP719s had not been

21

characterized (Table 2), and although published Ring B-closing CYP719s had little activity on

22

nandinine when assayed here, our enzyme library contained CYP719s that could synthesize

23

stylopine from nandinine. Thus, a side-product generated through enzyme promiscuity was re-

14

ACS Paragon Plus Environment

Page 15 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Synthetic Biology

1

introduced back into the pathway through the activity of a second promiscuous enzyme. Using

2

complementary substrate acceptance profiles to achieve a single product of interest is still new,34

3

although we expect this approach will become more prominent as large enzyme libraries become

4

increasingly accessible through gene synthesis.

5

The conversion of norlaudanosoline to dihydrosanguinarine in this work compares

6

favorably to other reports of heterologous dihydrosanguinarine/sanguinarine synthesis.

7

Norlaudanosoline supplementation in this work is lower than the system published by Trenchard

8

et al. (10 µM vs. 2 mM), but yield is higher (10% conversion to dihydrosanguinarine and

9

sanguinarine vs. 0.012% conversion to sanguinarine).39 Following our identification of a

10

bottleneck at SPS,21 Trenchard et al. performed a 2x2 combinatorial search for a new CFS and

11

SPS, ultimately selecting CYP719A5 and CYP719A2, respectively.39 While both enzymes were

12

included in our screen, neither enzyme was ultimately selected; CYP719A5 did not have

13

sufficient activity on nandinine, and CYP719A2 did not have sufficient activity on scoulerine.

14

CYP719A2 also displays low fluorescence as a GFP fusion protein, whereas enzymes selected

15

for integration into the dihydrosanguinarine pathway tended to display the highest levels of

16

fluorescence amongst enzymes with any particular desired activity (see Supporting Results and

17

Discussion).

18

While the system here represents a 10-fold improvement over previous work,21 the

19

pathway can be improved further. Conversion of norlaudanosoline to dihydrosanguinarine was

20

10% at 10 µM norlaudanosoline. However, as supplemented norlaudanosoline concentrations

21

increased, so did buildup of the intermediate N-methylstylopine. This points to the downstream

22

enzyme N-methylstylopine hydroxylase (MSH), a member of the CYP82 family, as the next

23

target for improvement. A single MSH has been identified thus far,64 but other members of the

15

ACS Paragon Plus Environment

ACS Synthetic Biology

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

1

same family have been demonstrated to have activity on BIAs.65-66 The strategy of enzyme

2

selection and screening employed here may be able to improve this next step as well.

3

Commonly-cited strategies for heterologous pathway optimization include control of

4

enzyme transcription and translation, as well as spatial control of enzymes through the use of

5

scaffolds or targeting signals67. The work presented here adds to the growing number of studies

6

that highlight the potential of transcriptome libraries to provide new solutions to these problems.

7

As the cost of DNA synthesis continues to drop, the screening of enzyme homolog libraries

8

should be considered an integral part of pathway engineering.

Page 16 of 38

9 10

Materials and Methods

11

Transcriptome data analysis

12

The PhytoMetaSyn database (www.phytometasyn.ca) of assembled transcriptome data

13

from BIA-producing plants3 was reverse-translated into putative ORFs using OrfPredictor.68

14

Translated ORFs were scanned for motifs of interest generated from sequence alignments of

15

published proteins and candidate sequences kindly provided by Dr. Peter Facchini (University of

16

Calgary) (NMTs: ERAQI(K/Q)DG; CYP719s: FxxGxxxCxG, PxIGN). Putative NMTs

17

identified from the library were aligned with published sequences and candidate sequences and

18

phylogenetic trees were generated and a subset was manually selected for testing. Putative

19

CYP719s, published CYP719s, and CYP719s deposited online on GenBank and the Cytochrome

20

P450 Homepage69 were grouped by BLASTclust into groups with 55% sequence identity at the

21

amino acid level of over 95% of the sequence. These groups were then aligned and phylogenetic

22

trees were generated in the same manner as for NMTs, and a subset was manually chosen for

23

testing. All alignments were performed with MUSCLE and phylogenetic trees were built by the

16

ACS Paragon Plus Environment

Page 17 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Synthetic Biology

1

neighbor-joining method using the program MEGA6.035 with a bootstrap value of 1000. When

2

indicated, branches with less than 80% confidence values were condensed to build condensed

3

phylogenetic trees. Enzymes selected for study were codon-optimized for expression in yeast and

4

synthesized by Gen9 (Cambridge, MA).

5 6

Construction of plasmids

7

All cloning was performed via yeast homologous recombination46 using regions of

8

homology added to DNA during PCR (see primer list, Supporting Table S2). PCR of DNA to be

9

cloned was performed with Phusion polymerase. When appropriate, E. coli was cultivated in LB

10

medium at 37°C with shaking at 200 rpm with supplementation of 100 µg/mL ampicillin as

11

necessary. A series of vectors designated as pBOT (Supporting Figure S1) were designed to

12

facilitate gene expression and enzyme activity assays. Each pBOT vector has a unique

13

combination of yeast selection marker, promoter, and terminator (Supporting Table S3). Details

14

on the construction of the pBOT vectors is outlined in Supporting Materials and Methods. To

15

switch selectable markers of Ring B-closing CYP719s, the promoter-gene-terminator cassette

16

from pBOT-Trp was introduced to pBOT-Leu by digesting both vectors with NotI/AscI, gel

17

purifying the pBOT-Trp insert and the pBOT-Leu vector, and ligating both fragments.

18

Dihydrosanguinarine pathway genes other than those purchased in this study were cloned

19

into either pGREG or pYES vectors (Supporting Table S3), where they could be used for activity

20

assays and/or genomic integration. Promoters, genes, and terminators introduced into pGREG or

21

pYES vectors were amplified with overlapping homology regions (Supporting Table S2) and

22

cloned by homologous recombination. Heterologous DNA (linkers C1, C6, H1, and H2) were

23

added as previously described21 (indicated in bold in Supporting Tables S2, S3 and S4).

17

ACS Paragon Plus Environment

ACS Synthetic Biology

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 38

1 2

Chromosomal integration of genes and multi-gene pathways

3

To facilitate activity assays and to build a dihydrosanguinarine synthesis strain, some

4

genes and multi-gene pathways were integrated into the genome of S. cerevisiae into sites

5

previously determined to allow high levels of gene expression70 (sites and genes are indicated in

6

Supporting Table S4). Strains built to facilitate activity assays (strains GCY1333, GCY1270, and

7

GCY1317) were built using homologous recombination and selected with antibiotic resistance to

8

200 µg/L geneticin and/or hygromycin using the antibiotic markers kanMX and hphNT1,

9

respectively.71-72 The dihydrosanguinarine synthesis strain (strain GCY1440) was built using

10

homologous recombination and CRISPR-Cas9 (Supporting Table S4). Regions of DNA (~500

11

bp) upstream and downstream of the integration site (UP and DOWN regions, respectively) were

12

amplified with homology to heterologous DNA to guide gene integration. Dihydrosanguinarine

13

synthesis genes used to build strains GCY133, GCY1270 and GCY1317 were excised from

14

pGREG vectors using AscI/NotI and gel purified, while dihydrosanguinarine synthesis genes

15

used to build strain GCY1440 were amplified from plasmids and heterologous linkers were

16

added (LV3, LTP1, LTP2, LV5). When applicable, Cas9 was directed to the 5´ and 3´ ends of

17

integration sites using two guide RNAs (gRNAs). The 20-bp targeting regions of gRNAs were

18

introduced through splice overlap extension using the primers indicated in Supporting Table S2.

19

Both gRNAs, along with linearized vector containing Cas9 (pCAS-Tyr)73 were co-transformed

20

into S. cerevisiae. Successful gene integration by both methods was confirmed by PCR. All

21

primers used are listed in Supporting Table S2.

22 23

BIA culture supplementation assays

18

ACS Paragon Plus Environment

Page 19 of 38

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Synthetic Biology

1

Activity of CYP719s and NMTs on BIAs was analyzed using culture substrate

2

supplementation assays. Yeast cultures were grown in yeast nitrogen base with 2% glucose and

3

amino acid dropout media as appropriate (YNB) at 30°C and 200 rpm. Yeast cells harboring

4

enzymes of interest were inoculated in triplicate into 100 µL of media in 96-well 2 mL deep-well

5

plates and incubated overnight. The following day, 900 µL of fresh media was added (1:10

6

dilution) and cultures were incubated for an additional 6 hrs. Cells were pelleted by

7

centrifugation for 5 min at 3,200 g and supernatants were aspirated. Cell pellets were suspended

8

in 300 µL TE (10 mM Tris, 1 mM EDTA, pH 8) containing BIAs as appropriate, at a

9

concentration of 5 µM unless otherwise specified, and incubated overnight at 30°C with shaking.

10

The following day, cells were pelleted at 3,200 g. Supernatants were transferred to 96-well

11

microtiter plates, diluted 1:1 in 100% methanol and clarified at 3,200 g prior to analysis by LC-

12

MS. To extract BIAs from cells, pellets were suspended in 300 µL methanol, vortexed at 1,000

13

rpm and 4°C for 30 min and centrifuged at 3,200 g. The resulting extracts were analyzed by LC-

14

MS as described below. To stay in the linear range of the LC-MS, samples supplemented with

15

>10 µM of a BIA were diluted before analysis. (S)-Scoulerine and (S)-stylopine were purchased

16

from ChromaDex (Irvine, CA, USA) and (R,S)-norlaudanosoline was purchased from Enamine

17

Ltd. (Kiev, Ukraine). Dihydrosanguinarine was prepared from sanguinarine by NaBH4

18

reduction74.

19 20

Liquid chromatography-mass spectrometry

21

Using a Perkin Elmer SERIES 200 Micropump, 5 µL of samples were injected onto an

22

Agilent Zorbax Rapid Resolution HT C18 2.1*30mm, 1.8 micron column and analytes were

23

separated using reverse-phase HLPC using the following gradient: Solvent A, 0.1% formic acid;

19

ACS Paragon Plus Environment

ACS Synthetic Biology

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

1

Solvent B, 100% acetonitrile, 0.1% formic acid; 0-1 min 95% A, 1-8 min 5 to 100% B (linear

2

gradient), 8-9 min 100% B; 9-9.1 min 95% A, followed by a 2 min equilibration at 95% A.

3

HPLC-grade methanol and acetonitrile were purchased from Fischer Scientific, and HPLC-grade

4

water and formic acid were purchased from Fluka. Sample elution was followed by injection into

5

the 7T-LTQ FT ICR mass spectrometer (Thermo Scientific) under the following conditions:

6

resolution, 50000 at 400 m/z; scanning range, 150-500 AMU; source voltage, 4.9 kV; source

7

temperature, 380°C; AGC target for full mass spectrum was set to 1 x 106 ions. Retention time,

8

exact mass (