pYPKa_A_ScRKI1

Freezer ID

Strain#..: ?
Box......: ?
Position.: ?

Importing the pydna package. Pydna is open source, documentated here and has a support forum as well as a publication:

Pereira F, Azevedo F, Carvalho Â, Ribeiro GF, Budde MW, Johansson B: Pydna: a simulation and documentation tool for DNA assembly strategies using python. BMC Bioinformatics 2015, 16:142.


In [1]:
from pydna.all import *

Read the cloning vector from a local file.


In [2]:
pYPKa =read("pYPKa.gb")

This vector should be a circular 3128 bp DNA molecule.


In [3]:
pYPKa


Out[3]:

The circular seguid checksum of pYPKa should be

aV1eIrzOiCjvw01yvKkxDXHKLMk


In [4]:
pYPKa.cseguid()


Out[4]:
aV1eIrzOiCjvw01yvKkxDXHKLMk

Importing the restriction enzyme to be used for cloning from Biopython.


In [5]:
from Bio.Restriction import AjiI

With the recognition site


In [6]:
print(AjiI.site)


CACGTC

This enzyme should cut only once in pYPKa.


In [7]:
pYPKa.cut(AjiI)


Out[7]:
(Dseqrecord(-3128),)

pYPKa is digested with AjiI to form a linear DNA fragment.


In [8]:
pYPKa_AjiI = pYPKa.linearize(AjiI)

The insert comes from a Genbank record. Access to Genbank is needed in order to download the template. If you execute this script, change the email address below to your own. Always tell Genbank who you are, when using the service.


In [9]:
gb =Genbank("myemail@home.se")

The template is downloaded from Genbank below.


In [10]:
template = gb.nucleotide("NM_001183514.1")

The template length and a representation:


In [11]:
print(len(template))
print(str(template.seq))


777
ATGGCTGCCGGTGTCCCAAAAATTGATGCGTTAGAATCTTTGGGCAATCCTTTGGAGGATGCCAAGAGAGCTGCAGCATACAGAGCAGTTGATGAAAATTTAAAATTTGATGATCACAAAATTATTGGAATTGGTAGTGGTAGCACAGTGGTTTATGTTGCCGAAAGAATTGGACAATATTTGCATGACCCTAAATTTTATGAAGTAGCGTCTAAATTCATTTGCATTCCAACAGGATTCCAATCAAGAAACTTGATTTTGGATAACAAGTTGCAATTAGGCTCCATTGAACAGTATCCTCGCATTGATATAGCGTTTGACGGTGCTGATGAAGTGGATGAGAATTTACAATTAATTAAAGGTGGTGGTGCTTGTCTATTTCAAGAAAAATTGGTTAGTACTAGTGCTAAAACCTTCATTGTCGTTGCTGATTCAAGAAAAAAGTCACCAAAACATTTAGGTAAGAACTGGAGGCAAGGTGTTCCCATTGAAATTGTACCTTCCTCATACGTGAGGGTCAAGAATGATCTATTAGAACAATTGCATGCTGAAAAAGTTGACATCAGACAAGGAGGTTCTGCTAAAGCAGGTCCTGTTGTAACTGACAATAATAACTTCATTATCGATGCGGATTTCGGTGAAATTTCCGATCCAAGAAAATTGCATAGAGAAATCAAACTGTTAGTGGGCGTGGTGGAAACAGGTTTATTCATCGACAACGCTTCAAAAGCCTACTTCGGTAATTCTGACGGTAGTGTTGAAGTTACCGAAAAGTGA

The seguid checksum of the template should be

insertseguid

Two primers are used to amplify the insert:


In [12]:
f,r =parse(""">555_RKI1_YPK_fwd
                      CGCGCCAAAATGGCTGCCGGTGTCCCAAA
                     >554_RKI1_YPK_rev
                      GGTCACTTTTCGGTAACTTCAACACT""", ds=False)

insert =pcr(f, r, template)

The primers anneal on the template like this.


In [13]:
insert.figure()


Out[13]:
         5ATGGCTGCCGGTGTCCCAAA...AGTGTTGAAGTTACCGAAAAGTGA3
                                 |||||||||||||||||||||||| tm 58.0 (dbd) 64.0
                                3TCACAACTTCAATGGCTTTTCACTGG5
5CGCGCCAAAATGGCTGCCGGTGTCCCAAA3
          |||||||||||||||||||| tm 65.2 (dbd) 75.2
         3TACCGACGGCCACAGGGTTT...TCACAACTTCAATGGCTTTTCACT5

A recombinant plasmid is formed by ligating the insert PCR product to the linear vector.


In [14]:
plasmid = (pYPKa_AjiI + insert).looped()

The plasmid sequence is rotated so that the origin is in the same position as for the cloning vector sequence.


In [15]:
pYPKa_A_ScRKI1 = plasmid.synced(pYPKa)

Stamp sequence with cSEGUID checksum


In [16]:
pYPKa_A_ScRKI1.stamp()


Out[16]:
cSEGUID_lWV6QjgzeLo047BeX0mXZF_PJ-8

Write sequence to a local file.


In [17]:
pYPKa_A_ScRKI1.write("pYPKa_A_ScRKI1.gb")




Calculate cseguid:


In [18]:
pYPKa_A_ScRKI1.cseguid()


Out[18]:
lWV6QjgzeLo047BeX0mXZF_PJ-8

Download

pYPKa_A_ScRKI1


In [19]:
from pydna.all import *
pYPKa_A_ScRKI1 =read("pYPKa_A_ScRKI1.gb")
assert pYPKa_A_ScRKI1.cseguid() in pYPKa_A_ScRKI1.definition