I am interested in understanding the allosteric effect caused by mutating a particular residue on ligand binding by tracking force propagation through the protein upon ligand binding. Could gromacs or any other program calculate force per residue from MD trajectories? Please advise.