Dear Justin,<br />
<br />
Thanks a lot for your reply.<br />
Here I am providing the data and explaining the problem in detail:<br />
<br />
I hope that it will help you in understanding in a better way.
Considering the tetramer as 'abcd'. 'ab' as one dimer and 'cd' as
another, I want to separate 'ab' and 'cd'. Length of each dimer is 8 nm,
so the tetramer (initial configuration) is approximately 16 nm and I
have used a cubic box of 30 nm. Two groups (reference group and pull
group) are decided on the basis of residue numbers forming 'ab' as
reference grp and 'bc' as pull group by making an index file. <br />
<br />
I have verified that the initial configuration does not split across
periodic boundaries, it is always inside the box throughout the
simulation. Moreover, the initial distance is not equal to box vector. <br />
<br />
By <span style="font-weight: bold;">g_dist</span> , I have tried to calculate the distance between the two groups and it came out to be approx. 15 nm (Please have a look at the <span style="font-weight: bold;">dist.xvg</span>
file at the end). Do you think that in the first ps time frame, should
it get separated by 15 nm?? Please correct me if I have done anything
wrong.<br />
<br />
By mistake, in my last mail I wrote that I need to pull along Y axis. In
actual, it is Z-axis and I have used pull_dim N N Y. So it is just the
mistake in writing the mail.<br />
<br />
Please let me know do you need any other data or parameters that I have used during this simulation.<br />
An early reply would me much appreciated.<br />
<br />
<br />
-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------<br />
-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------<br />
dist.xvg<br />
------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------<br />
<br />
# This file was created Tue Sep 13 11:09:23 2011<br />
# by the following command:<br />
# g_dist -f tub_pull.xtc -s tub_pull.tpr -n index.ndx <br />
#<br />
# g_dist is part of G R O M A C S:<br />
#<br />
# Giving Russians Opium May Alter Current Situation<br />
#<br />
@ title "Distance"<br />
@ xaxis label "Time (ps)"<br />
@ yaxis label "Distance (nm)"<br />
@TYPE xy<br />
@ view 0.15, 0.15, 0.75, 0.85<br />
@ legend on<br />
@ legend box on<br />
@ legend loctype view<br />
@ legend 0.78, 0.8<br />
@ legend length 2<br />
@ s0 legend "|d|"<br />
@ s1 legend "d\sx\N"<br />
@ s2 legend "d\sy\N"<br />
@ s3 legend "d\sz\N"<br />
0.0000000 15.0342455 -0.7714367 -0.9718266 14.9829559<br />
0.5000000 15.0345469 -0.7726345 -0.9766102 -14.9828863<br />
1.0000000 15.0345631 -0.7665710 -0.9862967 -14.9825792<br />
1.5000000 15.0343399 -0.7575226 -0.9916229 14.9824638<br />
2.0000000 15.0335970 -0.7490749 -0.9957008 -14.9818726<br />
2.5000000 15.0322971 -0.7446661 -0.9996347 -14.9805260<br />
3.0000000 15.0341187 -0.7364883 -1.0029345 -14.9825373<br />
3.5000000 15.0331717 -0.7328205 -1.0078001 14.9814405<br />
4.0000000 15.0311899 -0.7334194 -1.0118599 14.9791489<br />
4.5000000 15.0334759 -0.7351160 -1.0128736 14.9812908<br />
5.0000000 15.0336533 -0.7327585 -1.0102243 14.9817638<br />
5.5000000 15.0327044 -0.7342167 -1.0076447 14.9809132<br />
6.0000000 15.0335188 -0.7365904 -1.0090990 -14.9815168<br />
6.5000000 15.0296631 -0.7375288 -1.0111809 -14.9774609<br />
7.0000000 15.0294476 -0.7318525 -1.0155020 -14.9772310<br />
7.5000000 15.0322533 -0.7323456 -1.0173597 -14.9798965<br />
8.0000000 15.0353861 -0.7333517 -1.0172930 14.9829950<br />
8.5000000 15.0323915 -0.7285290 -1.0122719 14.9805660<br />
9.0000000 15.0284719 -0.7241497 -1.0072327 14.9771843<br />
9.5000000 15.0285101 -0.7232113 -1.0051012 14.9774113<br />
10.0000000 15.0306396 -0.7257242 -1.0020084 14.9796333<br />
10.5000000 15.0324621 -0.7242069 -0.9959688 14.9819384<br />
11.0000000 15.0324860 -0.7183657 -0.9881392 -14.9827623<br />
11.5000000 15.0308609 -0.7121477 -0.9848318 -14.9816465<br />
12.0000000 15.0318832 -0.7091122 -0.9842033 14.9828577<br />
12.5000000 15.0316772 -0.7063818 -0.9834919 -14.9828262<br />
13.0000000 15.0291939 -0.7024584 -0.9846754 -14.9804420<br />
13.5000000 15.0307999 -0.7052488 -0.9864941 14.9818020<br />
14.0000000 15.0286407 -0.7049980 -0.9902172 14.9794016<br />
14.5000000 15.0301456 -0.7010813 -0.9916439 -14.9810019<br />
15.0000000 15.0305882 -0.6967783 -0.9922857 -14.9816036<br />
15.5000000 15.0319242 -0.6941795 -0.9950562 14.9828815<br />
16.0000000 15.0315027 -0.6892128 -0.9928131 14.9828358<br />
16.5000000 15.0314274 -0.6889334 -0.9883585 -14.9830685<br />
17.0000000 15.0296593 -0.6906672 -0.9852924 14.9814167<br />
17.5000000 15.0277662 -0.6964712 -0.9844646 -14.9793034<br />
18.0000000 15.0301647 -0.7006731 -0.9863977 -14.9813862<br />
18.5000000 15.0285692 -0.7027130 -0.9882107 14.9795704<br />
19.0000000 15.0304384 -0.7040758 -0.9874868 -14.9814291<br />
19.5000000 15.0310335 -0.7008915 -0.9824295 -14.9825077<br />
20.0000000 15.0304193 -0.6977320 -0.9794283 14.9822359<br />
20.5000000 15.0306473 -0.6955061 -0.9762669 14.9827747<br />
21.0000000 15.0293531 -0.6917162 -0.9746027 14.9817610<br />
21.5000000 15.0295181 -0.6870613 -0.9746943 -14.9821339<br />
22.0000000 15.0302134 -0.6855507 -0.9771519 14.9827404<br />
22.5000000 15.0306492 -0.6869965 -0.9785366 14.9830208<br />
23.0000000 15.0285063 -0.6840801 -0.9742985 -14.9812803<br />
23.5000000 15.0265427 -0.6824503 -0.9738894 14.9794121<br />
24.0000000 15.0262260 -0.6803436 -0.9716139 14.9793377<br />
24.5000000 15.0266447 -0.6778135 -0.9723339 -14.9798260<br />
25.0000000 15.0236931 -0.6710443 -0.9711208 -14.9772491<br />
<br />
<br />
<br />
<br />
Best regards,<br />
Shilpi Chaurasia<br />
<br />
<br />
<br />
<br />
<br />
<br /><br /><span>On 09/09/11, <b class="name">"Justin A. Lemkul" </b> <jalemkul@vt.edu> wrote:</span><blockquote cite="mid:4E6A4D6E.2090505@vt.edu" class="iwcQuote" style="border-left: 1px solid rgb(0, 0, 255); padding-left: 13px; margin-left: 0pt;" type="cite"><div class="mimepart text plain"><br /><br />Shilpi Chaurasia wrote:<br />> I am trying to pull apart two dimers of tubulin protein joined together, to form two separate dimers by steered molecular dynamics using Pull code. The objective is to separate the dimers by pulling along one axis only (Y axis). I ran the simulation for 25 ps and it accomplished well but while visualizing the trajectory in VMD, in the very beginning (first frame), dimers are at distance which is much greater than expected. I expected that in first frame they should not have separated too much (i.e. in just 1 ps, they are at greater distance) I guess this is unusual. Please let me know where I am doing wrong? I have used the following pull code:<br />><br /><br />Actual numbers would help. Is your initial configuration simply split across periodic boundaries? If the initial distance is roughly equivalent to one of the box vectors, then that's the case. Otherwise, please provide real data and check to make sure your pull groups are what you think they are.<br /><br />>pull = umbrella<br />>pull_geometry = distance<br />>pull_dim = N N Y<br /><br />Note that you're not actually pulling along Y as you state above. Here, you're pulling along Z. Perhaps that's one issue.<br /><br />-Justin<br /><br />>pull_start = yes pull_init1 = 0.0<br />>pull_ngroups = 1<br />>pull_group0 = dimer_ref ; reference group (based on index file)<br />>pull_group1 = dimer_pull ; pull group<br />>pull_rate1 = 0.1 pull_k1 = 1000 pull_nstfout = 10 pull_nstxout = 10<br />><br />>best regards,<br />>Shilpi Chaurasia<br />><br /><br />-- <br />========================================<br /><br />Justin A. Lemkul<br />Ph.D. Candidate<br />ICTAS Doctoral Scholar<br />MILES-IGERT Trainee<br />Department of Biochemistry<br />Virginia Tech<br />Blacksburg, VA<br />jalemkul[at]vt.edu | (540) 231-9080<br /><a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="l">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br /><br />========================================<br />-- <br />gmx-users mailing list gmx-users@gromacs.org<br /><a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="l">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br />Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="l">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br />Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-request@gromacs.org.<br />Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="l">http://www.gromacs.org/Support/Mailing_Lists</a><br /></div></blockquote>