Rename HIS residues based on protonation

Question

[davidsaez]

Hello everyone,

I have an already protonated .pdb file. There are different types of histidine, (delta and epsilon-protonated) that can be identified based on the presence of HD1 or HE2, but all the residues are named HIS. I want to, based on the presence or absence of these hydrogen atoms, rename these residues to HSD or HSE in order to write a new pdb file.

Is there a way to rename residues in BioPython?

Thanks for your help.