Switching to h-bond constraints is likely incorrect for the pmx based ligand hybrid topologies. In case, a hydrogen atom is mapped to a dummy in stateA, the constraint will not be generated for the bonds involving this dummy.
I would suggest trying an SD integrator (note that you need to change tau-t in this case).