This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Getting sequence information by using atoms

I have a cif format file and its contents are as follows.

My purpose is to get sequence information (amino acid alphabet) from atoms information. Is there a function written in Python to achieve this?

enter image description here

atoms python cif protein
• 1,888 views • link
updated by Joe Joe • written by drbioengineer • drbioengineer 
awk '($1=="ATOM") {print $6}' < in.cif | paste -s -d '-'

??

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Please do not post images of the data. Always post data it self and expected output for better understanding the issue.

• 0 views • link
written by cpad0112 cpad0112

1 answer

You should be able to use BioPython for this, and it has functions already for this type of thing. Alternatively you could always use something like UCSF Chimera or pymol.

If you want to use generic commandline tools, as Pierre pointed out, the information you need can be obtained from column 6.

Note however, it may not be this simple to extract a sequence from this, as crystal files often have discontinuities in the sequences when compared to their genomic annotations.

• 0 views • link
written by Joe Joe

Log in to answer this question.

More posts like this