- for i in range(0,len(ADN_5),3):
- Codons.append (ARNm[i]+ARNm[i+1]+ARNm[i+2])
Well this is now the most fastidious part of the programming typing ;)
This is the result in the Pytho shell:
We have now our Codons list and ready to translate each of them to a START,END,or an aminoacid working should follow this rule table according to the combination of the the triplet:
Ala GCU GCC GCA GCG
Arg CGU CGC CGA CGG AGA AGG
Asn AAU AAC
Asp GAU GAC
Cys UGU UGC
Gln CAA CAG
Glu GAA GAG
Gly GGU GGC GGA GGG
His CAU CAC
Ile AUU AUC AUA
Leu CUU CUC CUA CUG UUA UUG
Lys AAA AAG
Met/START AUG
Phe UUU UUC
Pro CCU CCC CCA CCG
Ser UCU UCC UCA UCG AGU AGC
Thr ACU ACC ACA ACG
Trp UGG
Typ UAU UAC
Val GUU GUC GUA GUG
Met AUG
Stopcodons UAA UAG UGA
- for i in range(len(Codons)):
- if Codons[i]=="AUG":
- Met_nb=Met_nb + 1
- for i in range(i,len(Codons)):
- if Codons[i]=="UUU" or Codons[i]=="UUC":
- display( "----------->Phe")
- Protide.append("Phe")
- Phe_nb=Phe_nb + 1
- if Codons[i]=="AUG":
- display( "----------->Met")
- Protide.append("Met")
- if Codons[i]=="UGA" or Codons[i]=="UAG" or Codons[i]=="UAA":
- display( "Traduction Stopped")
- Protide.append("STOP")
- STOP_nb=STOP_nb + 1
We created Counter's variables like Met_nb for the Met numbers,Phe_nb for the Phe numbers wich they increment each time we found a codon that much,like that we know the number of each aminoacid generated by this DNA sequence.
In the next Code fragment we will look for functionnel protides wich are those who begin with (Met/START AUG) and end with (STOP UGA UAG UAA) each time we found the START followed by the STOP codon we add to All_Protide list the whole fragment betwin STARTand STOP:
- for i in range(len(Protide)):
- if Protide[i]=="Met":
- j=i+1
- for j in range(len(Protide)):
- if Protide[j]=="STOP":
- All_Protide.append(Protide[i:j])
- break
- else:
- pass
Python Shell result:
Until this point we got our All_Protide list filled by the functionnel protides,the next step is to do some basic statistc graph to have a quantitative represention of this DNA fragment in term of nucleotid and aminoacid frequency.
Aucun commentaire:
Enregistrer un commentaire