Question: After Extracting The Primary Structure Of A Multi-Chain Protein From Its Pdb File, How To Display In Fasta Format??
0
gravatar for Coli
9.2 years ago by
Coli0
Coli0 wrote:

After extracting the primary structure of a multi-chain protein from its PDB file cannot find a way to display it in FASTA format?? If the pdb text file has only one chain is pretty easy to put all the info in a here document and print staff on the screen as well as storing the same output in file.

However if you have more than one chain in the pdb file I can get the desired chains with each chain letter and length of each chain as well and print the info on the screen while iterating with a loop to get each chain. How do I send this output in a file that I wanna create?? with one chain u simply do:

my $fasta =<<FASTA;
>$id|$title|$chain_length aa
$sequence
FASTA
print $fasta;

How you do that when you have multiple chain and you need to display one chain at a time in a fasta format??????

ADD COMMENTlink modified 6.9 years ago by bluewoodtree60 • written 9.2 years ago by Coli0

Why did you (only) tag this question 'bioperl'? Do you want to use Bioperl to solve the problem? There is no mention of it in your question.

ADD REPLYlink written 9.2 years ago by Neilfws48k

Also, I reformatted your code, but it does not make much sense. Perhaps you could edit it to show what you are trying to achieve. Lines with code need 4 leading spaces and blank lines top and bottom.

ADD REPLYlink written 9.2 years ago by Neilfws48k
3
gravatar for Neilfws
9.2 years ago by
Neilfws48k
Sydney, Australia
Neilfws48k wrote:

Some (mostly non-coding) solutions to this problem are discussed in this blog post at the P212121 blog.

The simplest is probably to use the direct FASTA download from the PDB. For example, to get PDB code 4HHB, chain A:

curl -o 4HHBA.fa "http://www.rcsb.org/pdb/download/downloadFile.do?fileFormat=fastachain&compression=NO&structureId=4HHB&chainId=A"
ADD COMMENTlink modified 4 months ago by RamRS25k • written 9.2 years ago by Neilfws48k
0
gravatar for bluewoodtree
6.9 years ago by
bluewoodtree60
bluewoodtree60 wrote:

It might be a little bit late, but if someone is still looking for a solution: I wrote a web app for that:

http://bluewoodtree.zxq.net/pdb_to_fasta.html

ADD COMMENTlink written 6.9 years ago by bluewoodtree60
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 764 users visited in the last hour