NP-Sticky: Model-Guided Ligation Calculator


Sticky End Free Energy Calculator

There are two methods to obtain the free energy of ligation:
(1) Estimate using DNA thermodynamics (recommended for routine cloning).
(2) Calculate from two-piece ligation experiments by measuring output product distribution (recommended for library assembly applications).
The second method is more accurate for specific ligations conditions and is especially useful when the same sticky ends are used in multiple occasions (e.g., a standard display vector used for various applications).

By DNA Thermodynamics

In both examples above, the sticky end overhang is TCCA with a C before and an A after. (A sticky end overhang of TGGA with a T before and a G after would give the same result.) Mismatch scores are most accurate between sticky ends with the same number of overhang bases.

Please input 1 or 2 sticky ends (5'→3' upper strand) into the form below.

Sticky end overhang 1:
1 base before (5' side) sticky end 1:       1 base after (3' side) sticky end 1:
5' overhang 3' overhang

Sticky end overhang 2: (optional)
1 base before (5' side) sticky end 2:       1 base after (3' side) sticky end 2:
5' overhang 3' overhang

Custom stacking energy
(Stacking energy is the stabilizing interaction between adjacent base pairs that forms when the nicks between DNA fragments are repaired. The default is from Forián et al. J Phys Chem B, 103, 884-892 (1999); alternately, you can choose stacking energy from Protozanova et al. J Mol Biol, 342, 775-785 (2004).)

Optional: Custom free energy for the repair of phosphodiester bonds (default is -6.3 kcal/mol)
 kcal/mol



By Experiments

In the form below, please enter the experimental conditions and the results of 1-10 two-piece linear ligations experiments. The output mass ratio can be easily quantified by imaging the ligation products on an agarose gel. This algorithm calculates the sum of the squared residuals between model and experiment for a range of free energies. The free energy of the sticky end is the value that minimizes the sum of the squared residuals. A graph of the resulting residuals is provided to show if a true minimum is obtained.

Please input experimental conditions (up to 10 different trials) into the form below.

Number of base pairs of DNA 1:       Number of base pairs of DNA 2:
Number of overhang bases of sticky end:
Optional: Custom range of free energy: from kcal/mol to  kcal/mol

Experiment Input concentration (nM) Output mass ratio
DNA 1 DNA 2 DNA 1 DNA 2 DNA 1+2
1
2
3
4
5
6
7
8
9
10


Creative Commons License This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.
The Python source codes are available here: se_eval.cgi; se_fitting.cgi

Last modified by Daphne Ng on Jan 14th, 2014