Instead of trying to fix the molecule in place during step6_equilibration, could I let those steps run unconstrained, do one step in production and then edit the coordinates/velocities from my gas phase dynamics into the restart file? I realize that this will require some efforts with positioning in the membrane.

My other idea is read the initial coord/velocities into comparison, do the equilibration and then read the comp set back into the system. Would this work?