Designing plasmids encoding predicted pathways by using the BASIC assembly method

Overview
Creative Commons License: CC-BY Questions:
  • How to design plasmids encoding predicted pathways using BASIC assembly method ?

  • How to generate scripts enabling the automatic build of plasmids ?

Objectives:
  • Annotate the SBML pathway using Selenzy tool.

  • Design plasmids encoding predicted pathways using rpbasicdesign tool.

  • Generate dnabot scripts enabling the automatic build of plasmids.

Requirements:
Time estimation: 15 minutes
Supporting Materials:
Published: May 12, 2022
Last modification: Nov 3, 2023
License: Tutorial Content is licensed under Creative Commons Attribution 4.0 International License. The GTN Framework is licensed under MIT
purl PURL: https://gxy.io/GTN:T00273
version Revision: 18

Synthetic biology is a novel engineering discipline which requires computational tools for the design of metabolic pathways for the production of chemicals such as SynBioCAD portal which is the first Galaxy set of tools for synthetic biology and metabolic engineering (Hérisson et al. 2022).

In this tutorial, we will use a set of tools from the Genetic Design - BASIC Assembly Workflow (https://galaxy-synbiocad.org) which will enable you to design plasmids implementing metabolic pathways for the bioproduction of lycopene in E.coli (one of the preferred host cell for microbial biochemicals production).

Lycopene is a potent antioxidant and has been widely used in the fields of pharmaceuticals, nutraceuticals, and cosmetics. It’s widely found in fruits including tomato, watermelon, guava, and papaya but the extraction method of lycopene from these natural sources is expensive, complicated and cannot match the large market demand.

To address this demand, synthetic biology and metabolic engineering have been employed to develop microbial cell factories (e.g. E.coli strains) for lycopene production.

To design plasmids encoding lycopene bioproducing pathways, we will use the BASIC assembly method (as described in Storch et al. 2015) which relies on orthogonal linkers and type IIs restriction enzyme cleavage to provide a robust and accurate assembly of DNA parts into plasmid constructs. From these construct definition, the workflow will generate scripts enabling the automatic build of the plasmids as well as the transformation of strains using an Opentrons liquid handler robot. After downloading these scripts onto a computer connected to an Opentrons, one can perform the automated construction of the plasmids at the bench.

The workflow scheme we will use is shown below. First, we will run the steps of this workflow individually so as not to neglect the understanding of the intermediate steps as well. Then, we will run the workflow automatically so that it itself retrieves the outputs from the previous step and gives them as input to the next tool.

Genetic Design - BASIC Assembly Workflow.

Agenda

In this tutorial, we will cover:

  1. Data Preparation
    1. Get data
  2. Find enzymes using Selenzyme
  3. Generate potential constructs from SBML file using BasicDesign
  4. Generating python scripts using DNA-Bot
  5. Run the Genetic Design - BASIC Assembly Workflow

Data Preparation

First we need to upload and prepare the following inputs to analyze:

  • One SBML (Systems Biology Markup Language) file modeling a heterologous pathway producing lycopene such as those produced by the Pathway Analysis Workflow (https://galaxy-synbiocad.org).

  • The parts_for_lycopene.csv file listing the parts to be used (linkers, backbone and promoters) in the constructions.

  • Two YAML files providing two examples of settings, i.e. providing the identifiers of the laboratory equipment and the parameters to be used in the Opentrons scripts.

Get data

Hands-on: Data upload
  1. Create a new history for this tutorial named Genetic Design - BASIC Assembly Analysis.
  2. Import the input files from Zenodo:

    https://zenodo.org/record/6123887/files/dnabot_london_settings.yaml
    https://zenodo.org/record/6123887/files/dnabot_paris_settings.yaml
    https://zenodo.org/record/6123887/files/parts_for_lycopene.csv
    https://zenodo.org/record/6123887/files/rp_002_0011.xml
    
    • Copy the link location
    • Click galaxy-upload Upload Data at the top of the tool panel

    • Select galaxy-wf-edit Paste/Fetch Data
    • Paste the link(s) into the text field

    • Press Start

    • Close the window

  3. Rename the datasets if you want to shorten the file names. Rename rp_002_0011.xml to pathway.xml.

    • Click on the galaxy-pencil pencil icon for the dataset to edit its attributes
    • In the central panel, change the Name field
    • Click the Save button

Find enzymes using Selenzyme

At first, a pathway generated by the Pathway Analysis workflow (https://galaxy-synbiocad.org) is provided as input to the Selenzyme tool (Carbonell et al. 2018). Selenzyme searches for enzymes corresponding to each reaction of the pathway. It performs a reaction similarity search in the reference reaction database MetaNetX and outputs an updated SBML file annotated with the enzyme UniProt IDs. The tool provides several scores that can be combined in order to define an overall score. Scores are given for reaction similarity, conservation based on a multiple sequence alignment, phylogenetic distance between source organism and host, and additional scores calculated from sequence properties as shown in the example bellow.

Selenzyme Results.

Hands-on: Annotate enzymes with Uniprot IDs from a SBML pathway
  1. Run Selenzyme ( Galaxy version 0.2.0) with the following parameters:
    • param-file “Pathway (SBML)”: Select pathway.xml from your current history
    • In “Advanced Options”:
      • “Host taxon ID”: Leave the default value 83333. This stands for using E. coli as the chassis host.
      • “Comma separated taxon IDs of output enzyme sequences”: enter 553 which is the taxon ID of the Pantoea ananatis strain from which we want to extact enzymes
    Comment

    This tool will produce annotated pathway with UniProt IDs in XML and CSV format.

Question
  1. How are identified the enzyme in the SBML file?
  1. Enzyme are identified by their UniProt ID

Generate potential constructs from SBML file using BasicDesign

Second, the BasicDesign tool extracts enzyme IDs contained in the SBML (previously produced by Selenzyme) to generate genetic constructs compliant with the BASIC assembly approach. It uses as input an SBML file annotated with enzyme IDs for each reaction, and optionally one or several CSV files listing by their IDs the linkers, the promoters and the backbone used (parts_for_lycopene.csv). Example, below:

id,type,sequence,comment
L1,neutral linker,,
L2,neutral linker,,
L3,neutral linker,,

For linkers, the type annotation should be one of neutral linker, methylated linker, peptide fusion linker or RBS linker. For user parts, type should be one of backbone or constitutive promoter. Other type will raise a warning and will be omited.

BasicDesign converts the SBML file into CSV files describing the DNA-parts to be included into each construct (in an operon format, i.e. with only one promoter) and enumerate possible combinations of promoters, RBSs and enzymes into constructs. Depending on the numbers of enzymes per reaction, of RBSs and promoters available, and whether or not to perform CDS permutation within the operon, the number of constructs may vary.

Hands-on: Generate genetic constructs from the previously annotated SBML file
  1. Run BasicDesign ( Galaxy version 0.3.4) with the following parameters:
    • param-file “rpSBML file”: uniprot_ids (output of Selenzyme tool in xml format)
    • In “Advanced Options”:
      • param-file “Linkers and user parts”: Select parts_for_lycopene.csv from your current history in csv format
    Comment

    This tool will output a CSV constructs file (listing the constructs to be built), where each row corresponds to one construct and consists of a sequence of BASIC linker and DNA part IDs. The two other outputted CSV files provide the plate coordinates of (i) the BASIC linkers and (ii) the DNA-parts that the user will need to provide. Additionally, one SBOL (Synthetic Biology Open Language) file is produced for each construct generated. These files can be downloaded and visualized using online tools such as VisBOL.

Question
  1. How many constructs did you get ?
  1. 88 construct designs were generated in CSV and SBOL format.

Generating python scripts using DNA-Bot

In the last step, the DNA-Bot tool (Storch et al. 2020) reads the list of constructs (previously produced by BasicDesign) and the DNA-parts position on the source plates and generates a set of python scripts to drive an Opentrons liquid handling robot for building the the plasmids. Optional parameters can be set by the user to define the plastic labwares to be used, and set protocol parameters such as washing or incubation times for purification step (dnabot_paris_settings.yaml).

Hands-on: Generate DNA Bot python scripts
  1. Run DNA-Bot ( Galaxy version 3.1.0) with the following parameters:
    • param-file “Source Construct”: Constructs (output of BasicDesign tool)
    • param-files “Plate files”: Select these two files: User parts plate (output of BasicDesign tool), Biolegio plate (output of BasicDesign tool) with Multiple datasets.
    • In “Advanced Options”:
      • param-file “Yaml file providing labware IDs and parameter to be used (Optionally)”: Select either dnabot_paris_settings.yaml or dnabot_london_settings.yaml from your history)
  1. Click on param-files Multiple datasets
  2. Select several files by keeping the Ctrl (or COMMAND) key pressed and clicking on the files of interest

Comment

This tool will output DNA-Bot scripts in tar format. You need to download it and decompress the archive. After downloading these scripts onto a computer connected to an Opentrons, one can perform the automated construction of the plasmids at the bench. Additional metadata meaningful to keep track of parameters are also outputted by the tool. Dnabot Outputs. BASIC assembly workflow.

Question
  1. Looking at the scripts names, can you figure out the main steps of DNA-Bot ?
  1. DNA-Bot outputs python scripts that implement the 4 assembly steps: clip reactions, purification, assembly and strain transformation. In short, the Clip reactions step prepares the mixes for the ligation of the individual DNA parts with the linkers; the Purification step purifies the linker-ligated DNA parts using magnetic beads and the Opentrons magnetic module; the Assembly step mixes the DNA purified parts to build the final constructs; while the Transformation step transforms the chassis micro-organism with the plasmid and inoculates onto agar.

Run the Genetic Design - BASIC Assembly Workflow

In this section, you can run the Genetic Design - BASIC Assembly Workflow more easily and fastly following these instructions:

Hands-on: Execute the entire workflow in one go.
  1. Import the workflow into Galaxy

    Hands-on: Importing and Launching a workflow
    Click here to run Genetic Design - Basic Assembly Workflow workflow.

    Or copy the URL (e.g. via right-click) of this workflow, and import the workflow into Galaxy

    • Click on Workflow on the top menu bar of Galaxy. You will see a list of all your workflows.
    • Click on galaxy-upload Import at the top-right of the screen
    • Provide your workflow
      • Option 1: Paste the URL of the workflow into the box labelled “Archived Workflow URL”
      • Option 2: Upload the workflow file in the box labelled “Archived Workflow File”
    • Click the Import workflow button

    Below is a short video demonstrating how to import a workflow from GitHub using this procedure:

  2. Click on Workflow on the top menu bar of Galaxy. You will see Genetic Design - Basic Assembly Workflow
  3. Click on the workflow-run (Run workflow) button next to your workflow
  4. Provide the workflow with the following parameters:
    • param-file “Pathway (SBML)”: Select pathway.xml from your current history
    • “Host taxon ID”: Leave the default value 83333. This stands for using E. coli as the chassis host.
    • “Comma separated taxon IDs of output enzyme sequences”: enter 553 which is the taxon ID of the Pantoea ananatis strain from which we want to extact enzymes
    • param-file “Linkers and user parts”: Select parts_for_lycopene.csv from your current history in csv format
    • param-file “DNA-Bot settings (Optionally)”: Select either dnabot_paris_settings.yaml or dnabot_london_settings.yaml from your history)
    Comment

    All the outputs will be automatically generated and identical to the previous ones.