[BioRuby] How to renumber a PDB file using BioRuby ?

K. Shameer shameer at ncbs.res.in
Thu Apr 10 08:17:49 UTC 2008


Dear Alex,

Thanks for your help

> You have to alter the number by atom rather than by residue. Altering
> the residue object directly used to work, but has been broken. Anyway,
> this seems to work for your purpose:
>
> require 'bio'
>
> pdb = Bio::PDB.new(File.read(ARGV[0]))
>
> pdb.each_atom do |atom|
>    atom.resSeq += 12
> end
> puts pdb

Sorry, this is exactly I want to do, I want to apply it in a whole atom
level, not residue level - this works well for me !

Once again Thank you Alex and Chen for quick and great help,
K. Shameer

>
> It would be good (IMHO) to have:
>
> require 'bio'
>
> pdb = Bio::PDB.new(File.read(ARGV[0]))
>
> pdb.each_residue do |residue|
>    residue.resSeq += 10
> end
>
> puts pdb
>
> Do the same thing, but currently changes to residue objects don't
> propagate back down to the atoms so the above code doesn't do anything
> useful.
>
> Alex Gutteridge
>
> Department of Biochemistry
> University of Cambridge





More information about the BioRuby mailing list