I have simulated a protein-ligand complex by using gromacs software and gromos 54A7 force field was deployed for this purpose. Ligand topology (.itp) file was prepared by usimg ATB server. Now my questin is that can I use single trajectory method of gmx_MMPBSA tool to calculate protein-ligand binding free energy?