DNA Cauldron (complete documentation here) is a Python library to simulate restriction-based assembly operations. Provided the sequences of genetic parts and receptor vectors, DNA Cauldron will compute the assembli(es) that could result from the mix.
DNA Cauldron was written for Synthetic Biology applications - typically, to predict and validate batches of parts-based assemblies. It is simple to use, plays well with BioPython, can import and export Genbank (it conserves all features), and provides advanced methods such as connector part auto-selection, backbone selection for linear parts, methods to select constructs subsets when dealing with large combinatorial assemblies.
Try it online ! Use this web service to predict the outcome of a batch of (possibly combinatorial) Type 2S assemblies.
You can install DnaCauldron through PIP
sudo pip install dnacauldron
Alternatively, you can unzip the sources in a folder and type
sudo python setup.py install
It works better with the Networkx development version, that you install with
sudo pip3 install git+https://github.com/networkx/networkx.git
To assemble several parts and a receptor plasmid into a single construct, use single_assembly. The parts can be provided either as paths to genbank files or as Biopython records. Dna Cauldron returns a Biopython record of the final assembly, and (optionally) writes it to a Genbank file.
from dnacauldron.utils import single_assembly final_construct = single_assembly( parts=["partA.gb", "partB.gb", "partC.gb", "partD.gb", "receptor.gb"] outfile="final_construct.gb", # Name of the output enzyme="BsmBI" # enzyme used for the assembly )
The following example imports parts from Genbank files and outputs all possible outcomes of BsmBI-based Golden-Gate assembly as new genbank files 001.gb, 002.gb, etc. We ignore the final assemblies containing a BsmBI site as these are unstable.
from dnacauldron import (RestrictionLigationMix, NoRestrictionSiteFilter, load_record, write_record) # Load all the parts (including the receptor) parts_files = ["partA.gb", "partA2.gb", "partB.gb", "partB2.gb", "partC.gb", "receptor.gb"] parts = [load_record(filename, linear=False, id=filename) for filename in parts_files] # Create the "reaction mix" mix = RestrictionLigationMix(parts, enzyme='BsmBI') # Find all final assemblies (containing no sites from the restriction enzyme) assemblies = mix.compute_circular_assemblies() # Iter through all possible constructs and write them on disk as Genbanks. for i, assembly in enumerate(assemblies): out_path = os.path.join("..", "%03d.gb" % i) write_record(assembly, out_path, "genbank")
Full Assembly report
DNA Cauldron also implements routine to generate reports on the assemblies, featuring the resulting constructs (in genbank and PDF format) as well as figures for verifying that the parts assembled as expected and help troubleshoot if necessary.
The following code produces a structured directory with various reports:
from dnacauldron import load_genbank, full_assembly_report parts = [ load_genbank("partA.gb", linear=False, name="PartA"), load_genbank("partB.gb", linear=False, name="PartB"), load_genbank("partC.gb", linear=False, name="PartC"), load_genbank("receptor.gb", linear=False, name="Receptor") ] dc.full_assembly_report(parts, target="./my_report", enzyme="BsmBI", max_assemblies=40, fragments_filters='auto', assemblies_prefix='asm')
How it works
Dna Cauldron simulates enzyme digestions and computes sticky ends, then generates a graph of the fragments that bind together, and explores circular paths in this graph (which correspond to circular constructs), an idea also used in PyDNA and first described in Pereira et al. Bioinf. 2015 . DNA Cauldron adds methods to deal with combinatorial assemblies, selecting constructs based on a marker, routines for report generation, etc.
Dna Cauldron is an open-source software originally written at the Edinburgh Genome Foundry by Zulko and released on Github under the MIT licence (Copyright Edinburgh Genome Foundry). Everyone is welcome to contribute !
More biology software
DNA Cauldron is part of the EGF Codons synthetic biology software suite for DNA design, manufacturing and validation.