Modifying molecule coordinates to evaluate energetics

GROMACS version: 2024.1-Homebrew
GROMACS modification: No
I essentially want to vary the distance between two molecules in my simulation box and evaluate the energetics at each point. I plan on writing a python script to manually update the coordinates of the molecules of interest, but wondered if there was a better way of doing this.

that is what pullling simulation does, and umbrella sampling.