This vector is identical to pMEC1126 (pYPK0_TEF1_PsXYL1_TDH3tp) but for a N272D mutation in the XYL1 reading frame.
In [1]:
from pydna.all import *
In [2]:
pYPK0_TEF1_PsXYL1_TDH3 =read("pYPK0_TEF1_PsXYL1_TDH3.gb")
In [3]:
len(pYPK0_TEF1_PsXYL1_TDH3)
Out[3]:
In [4]:
# NBVAL_IGNORE_OUTPUT
from IPython.display import IFrame
IFrame('http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2988607', width='100%', height=250)
Out[4]:
In [5]:
pYPK0_TEF1_PsXYL1_TDH3.list_features()
Out[5]:
In [6]:
xyl1 = pYPK0_TEF1_PsXYL1_TDH3.extract_feature(5)
In [7]:
wt_protein_seq = str(xyl1.seq.translate())
In [8]:
print(wt_protein_seq)
In [9]:
print(wt_protein_seq[271])
In [10]:
pYPK0_TEF1_PsXYL1_TDH3.features[6].location.start
Out[10]:
In [11]:
codon = pYPK0_TEF1_PsXYL1_TDH3[1163+271*3:1163+271*3+3]
In [12]:
print(codon.seq)
print(codon.seq.translate())
In [13]:
candidate = (pYPK0_TEF1_PsXYL1_TDH3[:1163+271*3]+
Dseqrecord("GAC")+
pYPK0_TEF1_PsXYL1_TDH3[1163+271*3+3:]).looped()
In [14]:
new_feature = pYPK0_TEF1_PsXYL1_TDH3.features[8]
new_feature.qualifiers = {"label":"PsXYL1_N272D"}
candidate.features.append(new_feature)
In [15]:
candidate
Out[15]:
In [16]:
candidate.cseguid()
Out[16]:
In [17]:
candidate.id = "pMEC1127"
#candidate.description = "pYPK0_TEF1_PsXYL1_N272D_TDH3tp"
candidate.stamp()
Out[17]:
In [18]:
candidate.write("pYPK0_TEF1_PsXYL1_N272D_TDH3.gb")
In [19]:
reloaded = read("pYPK0_TEF1_PsXYL1_N272D_TDH3.gb")
In [20]:
reloaded.cseguid()
Out[20]: