Dear Gromacs users,<br><br>I am doing SMD of a ligand pathway, and then want ot do the PMF analysis. Initially I pulled the ligand from extra to intracellular side of the protein. The pull code used are given below.<br><br>
pull                     = umbrella<br>pull_geometry            = distance<br>pull_dim                 =  N N Y<br>pull_start               = yes<br>pull_nstxout             =  10<br>pull_nstfout             =  10<br>pull_ngroups             =  1<br>
pull_group0              =  r_57<br>pull_group1              =  r_C1<br>pull_rate1               =  0.005<br>pull_k1                  =  1000<br><br>here group0 is one residue at the extracellular region and r_C1 is the side chain Carbon atom of the ligand.<br>
<br>My question is , the total length of the channel is almost 40 Angstrom and as per my knowledge, when we did the pulling the pullx file will give the coordinate of the ligand through the channel. <br>Even though ligand reaches the opposite end of the channel, in pullx file I am getting around 19 Angstrom in total.<br>
<br>Again I calculate the g_dist of the ligand and the center of the channel, it shows around 40 Ang in total . Please find the links provided<br><br>Am I misunderstanding anything about the pullx file. <br><br>Could you please suggest me where thing get wrong?<br>
<br><a href="http://docs.google.com/fileview?id=0B1PyTWWGrqt6N2IyNjNhM2UtMzRiZC00MTJmLTg5ODItODJlMzQ5YTM3OGMy&amp;hl=en">http://docs.google.com/fileview?id=0B1PyTWWGrqt6N2IyNjNhM2UtMzRiZC00MTJmLTg5ODItODJlMzQ5YTM3OGMy&amp;hl=en</a><br>
<br><a href="http://docs.google.com/fileview?id=0B1PyTWWGrqt6ZWRlZGJkOTgtMGEwYi00MzE0LThiZTEtM2JjNmJkMWU5NDYy&amp;hl=en">http://docs.google.com/fileview?id=0B1PyTWWGrqt6ZWRlZGJkOTgtMGEwYi00MzE0LThiZTEtM2JjNmJkMWU5NDYy&amp;hl=en</a><br>
<br><br clear="all">Thank you,<br>-- <br>-Aswathy<br>