Just a note regarding OP1,OP2,O3P naming: this happens upon writing a pdb
with ambpdb. I presume the name is O1P, O2P, O3P in the prmtop. When you
write a pdb with ptraj, the names from the prmtop will be used.