Revision as of 18:10, 17 October 2018

CAT-Seq Esterase (Ready for Expression with RBS34)

CAT-Seq Esterase is a hydrolase that was used in Catalytic Activity Sequencing system for its promiscuous capability to hydrolyse N4-acyl-2'-deoxycytidine triphosphates (Substrate Nucleotides). It has been found that the enzyme accepts various N4-acyl-2'-deoxycytidine triphosphates as substrates for hydrolysis, leaving a 2'-deoxycytidine triphosphate (Product Nucleotide).

Sequence-wise it is the most similar to tannases (EC 3.1.1.20) and feruloyl esterases (EC 3.1.1.73), however to date there has been no enzyme characterized hydrolysing N4-amidic bond in modified deoxycytidine triphosphates.

It is the main component of Catalytic Activity Sequencing (CAT-Seq) method. CAT-Seq is a method for high-throughput catalytic biomolecule and genetic regulatory part activity-sequence relationship assessment toolkit.

See how this part is used in the CAT-Seq by pressing here!

Sequence and Features

Introduction

Biology

Description of the CAT-Seq esterase

CAT-Seq Esterase is a hydrolase that has a promiscuous capability to hydrolyse N4-acyl-2'-deoxycytidine triphosphates (Substrate Nucleotide). Sequence-wise it is the most similar to tannases (EC 3.1.1.20) and feruloyl esterases (EC 3.1.1.73), however to date there has been no enzyme characterized hydrolysing N4-amidic bond in modified deoxycytidine triphosphates.

To further elaborate on the fact, it has been found that the enzyme accepts various N4-acyl-2'-deoxycytidine triphosphates as substrates for hydrolysis, leaving a 2'-deoxycytidine triphosphate (Product Nucleotide).

To better understand the nature of this biological part, we have generated a structural homology model based on solved structures of homologous feruloyl esterases. Currently known feruloyl esterases contain a catalytic triad in their active sites consisting of Serine, Aspartate and Histidine. Aligning CAT-Seq Esterase to the structural 3WMT template, we could also identify the catalytic triad consisting of Ser195, Asp429, His467. The sequence similarities and the positions Catalytic Triad amino acids are highlighted in the illustration below (Fig. 1).

Figure 1. Superposition of the generated homology model (light blue) and the template used (yellow, PDB ID: 6G21). Catalytic triad residues are marked in red boxes with numbering. Calcium ion shown as grey sphere.

In-Silico design of the CAT-Seq esterase mutants

The reason of making esterase mutants was to benchmark the accuracy of CAT-Seq system. We have aimed to make Esterase mutants that would have slightly different activities. After deriving the mutations, we have measured their activities using standard methods, and then compared them with CAT-Seq measurements. Therefore, the esterase and its mutants can now be used for accuracy calibration when setting up the CAT-Seq for the first time.

The residues chosen for mutation were scattered around the catalytic triad of the esterase and involved mainly polar and aromatic side chains. In literature, some of these side chains are assumed to form bonds stabilizing the binding of the substrate (Fig. 2).

Figure 2. <b>View of the CAT-Seq Esterase Active Center. Residue coloring: catalytic triad - red, mutation positions - orange, calcium-coordinating residues - cyan. Catalytic residues as well as successful mutations are numbered in boxes. Hydrogen bonds in the catalytic triad represented as black dashed lines.

After carefully selecting the mutations in-silico, we have synthesized and measured the Activity of those esterase mutants in the laboratory. The mutations have successfully altered the activity of the original Esterase. Also, there were no cases of complete inactivation of enzyme catalytic activity from the mutants we managed test. Please head to the results section for further information.

Usage with CAT-Seq (Catalytic Activity Sequencing)

About CAT-Seq

CAT-Seq stands for Catalytic Activity Sequencing - a system designed and built for high-speed activity and interaction characterization of Catalytic and Regulatory biological parts. You can learn more about CAT-Seq [http://2018.igem.org/Team:Vilnius-Lithuania-OG by clicking this link]

Catalytic Activity Sequencing Overview

Figure 1. Abstract scheme of the Catalytic Activity Sequencing method

1. Library preparation - A library of catalytic biomolecules is prepared.
2. Library encapsulation into droplets - Every library fragment is physically separated by encapsulating them into picoliter water droplets. Also, substrate nucleotides, the targets for catalytic biomolecules, are encapsulated.
3. Catalytic biomolecule production - In each droplet catalytic biomolecules are produced.
4. Catalysis of the substrate conversion - Catalytic biomolecules may recognise the Substrate Nucleotides as a target for chemical reaction catalysis. Depending on biomolecule activity, a specific number of nucleotides with removed substrates (product nucleotides) is established in each droplet.
5. Activity Recording
   1. Droplet Merging - each of prior droplet is merged with new droplet that contains DNA amplification mix and reference nucleotides. The reference nucleotides are helping to tracking the Product Nucleotide number.
   2. DNA amplification - DNA is amplified using the different unique catalytic biomolecule DNA in each droplet. During the amplification, the Product Nucleotides and the Reference Nucleotides are incorporated into the DNA sequence.
6. Activity Reading by Nanopore Sequencing - All of the droplets are broken and the amplified DNA is sequenced. During the sequencing, biomolecule’s activity is retrieved by calculating reference and Product Nucleotides (substrate removed), together with the sequence of particular biomolecule variant.

Determining the accuracy of CAT-Seq

Trying to build a CAT-Seq pipeline in your own laboratory will require the CAT-Seq esterase in order to troubleshoot the system and assess the measurement accuracy and precision. In other words, the esterase can be used to calibrate the CAT-Seq.

Together with the Esterase, its substrate attached to a nucleotide is required (substrate nucleotide). In the standard case, the Substrate Nucleotide is N4-benzoyl-2'-deoxycytidine triphosphate. If the Esterase catalyzes the removal of the substrate from the nucleotide, it becomes the Product Nucleotide - 2'-deoxycytidine triphosphate.

Genetic Regulatory Part activity and cross-interaction assessment

While Catalytic Activity Sequencing began as a method for catalytic biomolecule activity recording, we have also create a way to adjust CAT-Seq to record activities of regulatory part. In addition to the activities, cross-interactions of different regulatory parts can also be measured.

When assessing the activities and sequences of libraries of catalytic biomolecules in CAT-Seq , the activity is measured and recorded as a function of Product Nucleotide that was produced in each droplet.

Yet, the activity of the catalytic biomolecule is not the only aspect that can influence the amount of Product Nucleotides that are produced. If all of the droplets would contain the same catalytic biomolecule , but each droplet would have a different concentration of that biomolecule, we would in result get different amounts of Product Nucleotides. For example, droplets with large amount of biomolecules may produce a large number of Product Nucleotides and vice-versa. The default and well-characterized Catalytic Biomolecule in CAT-Seq for regulatory part charectation would be the CAT-Seq Esterase.

Figure 1. A scheme which illustrates the Product Nucleotide concentration on the Catalytic Biomolecule amount in droplets.

Yet, the activity of the catalytic biomolecule is not the only aspect that can influence the amount of Product Nucleotides that are produced. If all of the droplets would contain the same catalytic biomolecule , but each droplet would have a different concentration of that biomolecule, we would in result get different amounts of Product Nucleotides. For example, droplets with large amount of biomolecules may produce a large number of Product Nucleotides and vice-versa. The default and well-characterized Catalytic Biomolecule in CAT-Seq for regulatory part charectation would be the CAT-Seq Esterase.

Example for RBS activity strength determination

We can now explore an example in which the CAT-Seq esterase can be used to assess the strengths of different ribosome binding site (RBS) sequences.

While keeping all the workflow from the main Catalytic Activity Sequencing design, all we need to do is to change the way we prepare the library. Firstly, instead of using a library of catalytic biomolecules, we can use a single CAT-Seq Esterase enzyme. Then, specific or random RBS sequence library can be prepared and attached to the esterase enzyme sequence. After that, every step of CAT-Seq method is the same.

Figure 1. Comparison between Standard CAT-Seq library constructs and Regulatory part adjusted CAT-Seq library constructs.

After the library is encapsulated into droplets, depending on the RBS sequence strength some droplets may have different amounts of the esterase. For example, the droplets with more esterase will produce more Product Nucleotides and vice-versa. After the amplification and sequencing steps we can learn two main things about the measured library - the RBS sequence and its relative RBS strength in the form of Product Nucleotide concentration for every library variant screened.

Example for Toehold cross-interaction determination

Next, we want to give an example on how to record regulatory part interactions using CAT-Seq. In this case, we will be using Toehold Switches.

The Toehold Switch systems are composed of two RNA strands referred to as the Switch and Trigger. The Switch RNA contains the coding sequence of the gene being regulated. The Switch RNA forms a hairpin structure that includes the RBS site. While the hairpin structure is formed, the translation is inhibited. The Trigger RNA is a molecule that can selectively bind to the Switch RNA region and expose the gene RBS site for ribosomes. Once that happens, the protein translation can be initiated.

Once again, only the first part of general CAT-Seq design needs to be changed - the library preparation. Instead of using a catalytic biomolecule library, a single enzyme is used. Then, libraries need to be prepared - one for Toehold Switches and another for Trigger RNA. The Trigger RNA libraries also require a separate T7 promoter for RNA expression. Then, both of those libraries must be ligated to the enzyme DNA fragment. Resulting library contains fragments which have the same catalytic biomolecule (In the general case - the CAT-Seq Esterase), yet each of them have a random combination of a specific Toehold Switch and RNA Trigger.

Figure 1. Comparison between Standard CAT-Seq library constructs and Regulatory part adjusted CAT-Seq library constructs (For cross-interaction recording).

After the encapsulation, the droplets which have the correct Trigger RNA and Toehold Switch combination (in which Trigger RNA binds to that specific Toehold Switch), can produce catalytic biomolecules . In turn, those catalytic biomolecules can produce the Product Nucleotides.

In other droplets, where the Trigger RNA binds to the Toehold Switch with weaker affinity, there are also less catalytic biomolecules and Product Nucleotides produced. Finally, the droplets in which the Trigger RNA does not bind to the Toehold Switch, there are no biomolecules or product nucleotides produced.

Part Characterization (Vilnius-Lithuania Overgraduate 2018)

Kinetic characterization of the CAT-Seq Esterase

Spectrofotometic kinetic data based on decay of absorbance at 310 nm due to substrate nucleotide catalytic conversion was gathered using a range of starting substrate nucleotide concentrations. A Michelis Menten curve for the CAT-Seq esterase was plotted using the aquired data. Data shown in 1 Fig., show perfect fit (R2 = 0.9449) to a standard Michaelis Menten curve.

Figure 1. Michaelis Menten Curve generated for CAT-Seq Esterase enzyme. The initial velocity of the enzyme was determined spectrophotometrically as a decay of absorbance at 310 nm at different starting substrate nucleotide concentrations. The graph shows the average values of three independent experiments.

Michaelis Menten plot transformations were generated for CAT-Seq esterase enzyme (BBa_K2621000). 2 Fig., display Lineweaver-Burk and Hanes-Woolf transformation plot. Based on the equation koeficients in the Hanes-Woolf transformation, the expermentally determined Vmax value is 17.2 µM/min, Km value is 86 µM.

Figure 1. Two Michaelis Menten Curve transformations for CAT-Seq Esterase enzyme. Lineweaver-Burk and Hanes-Woolf transformations for Michaelis Menten curve generated earlier. Both of the transformations were performed on data acquired from 3 independent experiments. Graphs show the fitted linear function and its correlation coefficient.

Esterase and its mutants catalytic activity determination

The 10 esterase mutants housing mutations at bioinformatically predicted sites were created. To see why and how mutations were designed in-silico, [http://2018.igem.org/Team:Vilnius-Lithuania-OG/Model please click this link]. Each of the mutant was constructed utilizing PCR and synthesized using In vitro transcription and translation kit and their catalytic activity towards N4-benzoyl-2'-deoxycytidine triphosphate were tested.

Figure 3. Relative activity of In silico generated CAT-Seq esterase mutants CAT-Seq esterase mutant sequences were generated In silico and synthesised using in vitro transcription and translation kit. The graph shows the relative catalytic activity of each generated mutant measured spectrophotometrically and the corresponding mutation site.The activity was normalized to wild type (WT) enzyme.

The reaction kinetics were measured using the spectrophotometer as a decrease of absorbance. Figure 3 displays the relative hydrolysis speed of each mutant generated. As seen from these results, a variety of mutants, showing different catalysis speeds were produced. Some of the amino acids changes affected the activity drastically, for example Trp224 to Tyr, Lys227 to Arg or Glu509 to Lys. Other mutations only modulated the activity Asn107 to Asp or Glu194 to Ala. Additionally, large 8 amino acid deletion at position Pro348-Hy356 caused only a moderate decrease in enzymes activity.

Esterase and mutants activity assessment using CAT-Seq

The constructed in silico designed mutant library was subjected to catalytic activity sequencing. By applying the data preparation and analysis pipeline, the mean methylation scores arising from different ratios of catalytically converted and reference nucleotides for each barcoded mutant DNA template were filtered and extracted from the DNA embedded with catalytic activity information (in a form of incorporated reference to catalytically converted nucleotide ratio). The collected data was normalized over Wild Type Esterase and K227R mutant (lowest activity).

Figure 1. Comparison of In bulk and CAT-Seq measured esterase mutant relative activity. In silico generated Esterase mutant library was subjected to catalytic activity sequencing. The mean methylation scores for each barcoded mutant DNA template were filtered and extracted. The collected data was normalized over Wild Type CAT seq Esterase and K227R mutant (lowest activity). The relative activity, extracted from the mean methylation score of each mutant read is compared to measurement data gathered in standard sized reactions (in bulk) .

The relative methylation score (reference nucleotide count) of each mutant read corresponds to the activity of the enzyme it encodes. The higher the activity of the expressed enzyme, the lower methylation score are assigned, due to catalytic conversion of substrate nucleotides. The comparison of the results, gathered with CAT-Seq catalytic activity sequencing method and in standard sized reactions (in bulk) spectrophotometric data (Fig 1.) conclude the viability of CAT-Seq approach. The activity reading, extracted from the DNA sequence correlates with the in bulk measurement data perfectly. The activity of the each Esterase mutant is measured accurately and is assigned to the corresponding DNA sequence.

Ribosome binding site activity measurements

The three ribose binding sites (BBa_B0030, BBa_B0032, BBa_B0034) already found in biobrick registry and a new CAT-Seq RBS BBa_K2621038 were place upstream the CAT-Seq esterase gene. The strength of each ribosome binding site was investigated by expressing the esterase enzyme using the In vitro transcription and translation kit and spectrophotometrically measuring the as a decrease of absorbance due to the hydrolysis of substrate nucleotides.

Figure 4. Relative activity esterase enzyme expressed using different RBS sites. CAT-Seq esterase mutant sequences were placed under 4 different RBS control. The enzyme was synthesised using cell free expression system and catalytic activity of each generated variant was measured spectrophotometrically as a decrease of substrate absorbance at 310 nm. The shown data is normalized to BBa_B0034 RBS activity.

Figure 4 displays the relative hydrolysis speed (normalized to BBa_B0034) of CAT-Seq esterase expressed using 4 different ribosome binding sites. Based on the fact that the same enzyme and substrate nucleotide concentration are used, the difference in activity corresponds to difference in synthesised biomolecule concentrations. The expression of an enzyme depends on the particular strength of used ribosome binding sites.

Taking this into consideration, the measured relative activity corresponds to strength of ribosome binding sites.

Ribosome binding site activity measurements using CAT-Seq

The constructed Ribosome Binding site library (BBa_B0030, BBa_B0032, BBa_B0034, BBa_K2621038 with a downstream CAT-Seq esterase gene BBa_K2621000) was subjected to catalytic activity sequencing method. The mean methylation scores (reference nucleotide count) for each barcoded DNA template, housing different RBS were filtered and extracted from the DNA embedded with catalytic activity information (in a form of incorporated reference to catalytically converted nucleotide ratio). The collected activity data was normalized to BBa_B0034 data and is shown in Fig. 2.

Figure 2. Comparison of in bulk and CAT-Seq measured ribosome binding site relative strength.The catalytic activity of esterase gene, regulated by a library of ribosome binding sites was measured using cell free expression system in standard sized reactions (in bulk) or CAT-Seq approach and compared side by side. The mean methylation scores for each barcoded mutant DNA template were filtered and extracted. The collected data was normalized BBa_B0034 corresponding to mean strength of 1.

Stronger ribosome binding sites increase the yield of translated proteins and in turn increase the number of catalytically converted substrate nucleotides. This increase is inversely proportional to the assigned mean methylation score. Based on this fact, the activity results can be extracted from mean methylation scores (reference nucleotide count) and correspond to ribosome binding site strength. The catalytic activity sequencing results were compared to earlier measured in bulk RBS strength results.

The comparison once again concludes the viability of CAT-Seq approach. The ribosome binding site strength, extracted from the DNA sequence reference nucleotide count correlates with the in bulk measurement data. These results display the validity of CAT-Seq as a method for screening the strength of regulatory sequences and its ability to assign accurate phenotype to genotype linkage.

Cross-interaction measurements of Toehold Switches

The 9 library members, composed of 3 unique Toehold Switch sequences and 3 unique activating RNA sequences, termed Trigger RNA, were designed to test the orthogonality and regulatory characteristics of each part.

Each of the regulatory part, consisting of one toehold sequence upstream of the CAT-Seq esterase gene (BBa_K2621000) and one trigger sequence were constructed. First of all, the orthogonality of each toehold:activating RNA pair and they regulatory characteristic have been tested in bulk.

The constructed library members were synthesized using In vitro transcription and translation kit and their catalytic activity towards N4-benzoyl-2'-deoxycytidine triphosphate were tested. The reaction kinetics were measured using the spectrophotometer as a decrease of absorbance due tue hydrolyzed substrate nucleotide.

Figure 5. The catalytic activity of each Toehold:Trigger RNA construct matrix.CAT-Seq esterase mutant sequences were placed downstream the riboregulatory sequences with corresponding trigger part and they catalytic activity was measured. The decrease of absorbance corresponds to catalytically active enzyme. The graph shown as a matrix concludes that only Toehold sequences expressed with their corresponding Trigger RNA produce an active enzyme molecule.

Figure 5 displays the relative hydrolysis speed of each regulatory part variant in a form of matrix. The decrease of absorbance shown in Y axis corresponds to the catalytic conversion of substrate nucleotides. As seen from control experiment, in which standard esterase was expressed, the decline of absorbance over time is seen.

Taking these results into consideration, the same decrease of absorbance is only seen in the diagonal of the matrix. This means, that active catalytic molecules is expressed only when both regulatory molecules of the same group are present - Toehold1 (BBa_K2621011) with Trigger1 (BBa_K2621014), Toehold2 (BBa_K2621012) with Trigger2 (BBa_K2621015) and Toehold3 (BBa_K2621013) with Trigger3 (BBa_K2621016). None of the regulatory sequences show any cross talk with the other group. These results conclude the generation of working toehold riboswitches control library for regulatory sequences parameter and orthogonality screening.

Cross-interaction measurements of Toehold Switches using CAT-Seq

In addition to ribosome binding sites, we have constructed Toehold regulatory sequence library constituted of different toehold and triggers pairs was constructed subjected to catalytic activity sequencing method:

The mean methylation scores (reference nucleotide count) for each barcoded DNA template, housing different regulatory sequence were filtered and extracted from the DNA embedded with catalytic activity information (in a form of incorporated reference to catalytically converted nucleotide ratio).

Figure 23. The evaluation of Toehold-Trigger riboregulatory sequence orthogonality using CAT-Seq.The catalytic activity of esterase genes, regulated by different Toehold switches were measured using CAT-Seq. The mean methylation scores for each barcoded regulatory construct DNA template was filtered and essigned. Low methylation scores correspond to actively expressed protein and are only assigned when both Toehold and trigger sequences from the same group are present verifying the already measured orthogonality of regulatory parts.

The graph displays the mean methylation (reference nucleotide) scores assigned to each barcoded toehold-trigger construct read. Based on the results, low methylation score are only assigned when both Toehold and trigger sequences from the same group are present, due tue esterase being expressed. These results correlate perfectly to the standard (not in droplet) measurement results, carried out earlier.

Based on this fact it can be concluded that CAT-Seq activity sequencing method can be utilized as a precise and accurate way to screen and assign the activity and orthogonality of regulatory sequences.

References