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.

More Marvin Bilog's questions See All
Similar questions and discussions