From: Vermaas, Josh (
Date: Wed Nov 17 2021 - 20:34:02 CST

Hi Georcki,

Serial number and indices aren't mutable, but are instead inferred from the order the atoms appear in the pdb. The way I'd do this is to use topotools to reorder the atoms so that atoms that match segid PROB end up first. Borrowed from;!!DZ3fjg!vFQmpulJhkl0pGH0Eq5OMI2ePWX4esX-m58Mm2ALa80HbLk-1S3p3GYlB7_V3UmX6g$ :

package require topotools
set prob [atomselect top "segid PROB"]
set notprob [atomselect top "not segid PROB"]
set mol [::TopoTools::selections2mol [list $prob $notprob]
animate write pdb reordered.pdb $mol


On 11/17/21 8:04 PM, Georcki Ropon wrote:
Dear user,

I want change serial atom number , my number start with 3000 but I want all number (para chain B) start 1 2 3 4 ….

Bellow my code but not fix it.

mol new prot_orient_add_ter.pdb

set prob [atomselect top "segid PROB"]

set serials [$prob get serial]

set count 0
foreach atom $serials {
$prob set serial $count
incr $count


set all [atomselect top all]
$all writepdb renum.pdb

Josh Vermaas
Assistant Professor, MSU-DOE Plant Research Lab and Department of Biochemisty and Molecular Biology<>;!!DZ3fjg!vFQmpulJhkl0pGH0Eq5OMI2ePWX4esX-m58Mm2ALa80HbLk-1S3p3GYlB7-qvz2IIg$