Optimize DNA sequences under constraints.
Project description
![[logo]](https://pypi-camo.freetls.fastly.net/663f4c97dca85c430dba69fefcfd72fb0989da44/68747470733a2f2f7261772e67697468756275736572636f6e74656e742e636f6d2f4564696e62757267682d47656e6f6d652d466f756e6472792f446e6143686973656c2f6d61737465722f646f63732f5f7374617469632f696d616765732f7469746c652e706e67)
DnaChisel
(Documentation in progress, come back later !)
DnaChisel (complete documentation here) is a Python library to optimize the nucleotides of DNA sequences with respect to a set of constraints and optimization objectives.
It can be used for codon-optimizing the genes of a sequence for a particular micro-organism, modifying a sequence to meet the constraints of a DNA provider while preserving genes, and other sensible patterns.
DnaChisel also provides much freedom to define optimization problems and model new kinds of specifications, making it suitable for either automated sequence design, or for complex custom design projects.
Example of use
In this basic example we generate a random sequence and optimize it so that
It will be rid of BsaI sites.
GC content will be between 30% and 70% on every 50bp window.
The reading frame at position 500-1300 will be codon-optimized for E. coli.
Here is the code to achieve that:
from dnachisel import *
# DEFINE THE OPTIMIZATION PROBLEM
problem = DnaOptimizationProblem(
sequence=random_dna_sequence(10000),
constraints=[AvoidPattern(enzyme_="BsaI"),
EnforceGCContent(mini=0.3, maxi=0.7, window=50)],
objectives=[CodonOptimize(species='e_coli', location=(500, 1300))]
)
# SOLVE THE CONSTRAINTS, OPTIMIZE WITH RESPECT TO THE OBJECTIVE
problem.resolve_constraints()
problem.optimize()
# PRINT SUMMARIES TO CHECK THAT CONSTRAINTS PASS
print(problem.constraints_text_summary())
print(problem.objectives_text_summary())DnaChisel implements advanced constraints such as the preservation of coding sequences, or the inclusion or exclusion of advanced patterns, as well as some common biological objectives (such as codon optimization, GC content), but it is also very easy to implement new constraints and objectives.
Installation
You can install DnaChisel through PIP
sudo pip install dnachiselAlternatively, you can unzip the sources in a folder and type
sudo python setup.py installTo be able to generate plots and reports, run
sudo pip install dna_features_viewer weasyprintLicense = MIT
DnaChisel is an open-source software originally written at the Edinburgh Genome Foundry by Zulko and released on Github under the MIT licence (¢ Edinburg Genome Foundry).
Everyone is welcome to contribute !
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
