Correct protonation states of HIS

Find HIS ND1 within 4 A of any metal, if ND1 is closer than NE2, write
Protonate/Deprotonate commands to inConstr, which is then used by DMD
in the creation of new.pdb to give correct protonation states.