Dear RG experts,
I'm playing with Meta dynamics for the first time and I'm struggling to come up with a way of rotating my ligand in the catalytic site of an enzyme. The ligand is a long narrow peptide, an alpha helix embedded within the catalytic site. I want to investigate the possibility of non-specific binding using the metadynamics capabilities of PLUMED patched onto Gromacs. I am currently trying the idea of defining four points, two on the ligand and two on the enzyme, then using them to define a torsion angle.
I've trawled the literature and I only seem to find 'in house' or mathematics to define specific actions or rotation/spin of one protein against the other. Unfortunately, none of these have been implemented in PLUMED as far as I can see.
I would really appreciate and thoughts on this problem.
Thanks
Anthony