Inputs needed:

1. Enzyme name(s)

2. How many cuts does the target DNA have for each selected enzyme? (optional, default 1)

3. How long is the target DNA (in base pairs)? (optional, default 5000)

4. How long are you going to incubate the reaction (in hours)? (optional, default 1 hour)

5. How much DNA do you want to cut (in µg)? (optional, default 1 µg)

Output:

1. Possible buffers in the order of suitability (or the result “simultaneous digest not recommended”)

2. Amount of each enzyme needed for all possible buffers

Requirements:

biopython

Since even the latest Biopython distribution doesn’t contain all enzymes sold by NEB, you need to manually update the Restriction_Dictionary.py manually with the file Restriction_Dictionary.py by copying it into the Bio/Restriction folder of the folder, where your python3 stores the python packages. On Ubuntu 16.04, this would be /usr/lib/python3/dist-packages/Bio/Restriction.

Files:

reoptimize.py The script that does the calculations. Usage examples:

This is a double digest with AflIII and HindIII, where the target DNA has two AflIII sites and one HindIII site:

reoptimize -e ‘AflIII 2’ ‘HindIII 1’

This gives all necessary parameters via the command line:

reoptimize -e ‘EcoRI 2’ ‘HindIII 3’ -l 3000 -t 4 -m 2

-l (length of target dna, in base pairs) -t (incubation time, in hours) -m (amount of DNA, in µg)

make_sqlite_database.py This script fetches all the data for NEB enzymes from the NEB web pages and assembles the database that is needed for the script to run. Running it results in the database file “REsqlite3.db”

assay_DNAs.fasta This files contains the full DNA sequences of all assay DNAs used by NEB. We include it here to avoid querying the “Frequency of restriction sites” table (which is anyway incomplete).