Hi everyone,<br><br>I am doing MD simulation with a protein-ligand system, and I want to pick out the water molecules (their residue numbers or coordinates in any frame) that simultaneously contact (within 0.4 nm range for heavy atoms) the ligand and the protein, so that I could plot the lifetime of these resident water molecules within a trajectory.<br>

<br>Does anyone know how to achieve this within GROMACS?<br><br>Thanks for any advice.<br>Yun<br>