Identify restriction sites of a cDNA sequence using a dictionary
2
0
Entering edit mode
7.5 years ago

Hi.

I need to identify restriction sites of a cDNA sequence in a fasta file. First, I transform the sequence from FASTA to a string and I have created a dictionary which contains recognition sequences and names of all restriction enzymes.

I want to identify these restriction sites using the dictionary... how can I do that?

Thank you.

restriction enzymes python biopython dictionary • 4.3k views
ADD COMMENT
0
Entering edit mode

Do your restriction sites have different length?

ADD REPLY
0
Entering edit mode

Yes, the sequence contains different restriction sites of different lengths

ADD REPLY
0
Entering edit mode

Then this is really hard to do because you don't know how much of your cDNA string you would want to use as a key for the dictionary (because of the different length). I would vote for trying something similar to what Torst is suggesting.

ADD REPLY
0
Entering edit mode

Thank you for your answer.

ADD REPLY
0
Entering edit mode
7.5 years ago
Torst ▴ 960

You need to search for each recognition sequence string within the cDNA string.

In python you would use the find() method in the String object.

ADD COMMENT
0
Entering edit mode

How can I do that if my string is a whole cDNA sequence and the dictionary contains ALL retsriction enzymes? I need to identify all the enzymes that cut the sequence...

ADD REPLY
0
Entering edit mode
7.5 years ago
Biojl ★ 1.7k

Is this a real task for your lab or just some class practice?

If you are in the first case you don't have to reinvent the wheel. This is already implemented in biopython...

http://biopython.org/DIST/docs/cookbook/Restriction.html

ADD COMMENT
0
Entering edit mode

It's a class practice. I've tried to do it with the Restriction module, but I couldn't, so I decided to create a dictionary with the enzymes...

Could you tell me how I should do this practice with the Retsriction module?

ADD REPLY

Login before adding your answer.

Traffic: 1755 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6