Hi Alexandra,
I am having the same issue with a lone pair on a chlorine atom on my ligand generated by CHARMM-GUI. Did you ever solve this? I am currently trying to increase the rcoulomb to bigger than the distance excluded distance for all the minimization/equilibration steps as was suggested at the end of this thread (link). But I’m not sure if this is a good idea.