Hello,
After doing 100 ns MD simulations to a peptide, I want to determine the trajectory between the side chains of some residue pairs.
For example:
RSFFSFLGEAGFDG - I want to plot distance vs time for the guanidino side chain of Arg-1 and carboxyl side chain of Asp-12.
The best thing I can do right now is to determine distance between alpha-carbons using the following commands:
gmx make_ndx -f md_protein.tpr -o 1_12_distance.ndx
3 & ri 1 12
gmx distance -f md_protein.xtc -s md_protein.tpr -oall 1_12_dist.xvg -oh 1_12_disthist.xvg -n 1_12_distance.ndx -tu ns -len 2 -binw 0.01
I chose 3 for alpha-C. I also tried choosing "Side chain" but it used all atoms in the side chains. Thank you for those who will give any suggestions.