<br><br><div class="gmail_quote">On Mon, May 17, 2010 at 12:14 PM, Aswathy <span dir="ltr"><<a href="mailto:ammasachu@gmail.com">ammasachu@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hi,<br><br>I tried to calculate PMF for the ligand transport using pull_geometry=position.<br><br>Let me explain you what I did so far,<br><br>I picked a small collection of backbone atoms nearly from the centre of the channel.My SMD was displacing the ligand from extracellular to intracellular. <br>
I checked the pullx file for getting the window spacing (pull_init ). <br>Initial structure is at =7.61<br>Centre of the channel(reference) =8.20<br>Last displacement is at =8.91<br>
<br>Therfore I selected pul_init as 0.6, 0.5, 0.4 ,....0....-0.4,-0.5,-0.6<br><br>For example, the first frame(at 7.61) was sampled using the following parameters,<div class="im"><br>pull = umbrella<br>
pull_geometry = position<br>
pull_dim = N N Y<br>pull_start = no<br></div>pull_nstxout = 10<br>pull_nstfout = 10<br>pull_ngroups = 1<br>pull_group0 = U_ref<br>pull_pbcatom0 = 0<br>
pull_group1 = r_C1<br>pull_pbcatom1 = 0<br>pull_init1 = 0 0 0.6<br>pull_k1 = 1000<br>pull_rate1 = 0<div class="im"><br>pull_vec1 = 0 0 0<br>
<br></div>The problem is that the ligand undergoes a displacement to wrong direction during sampling. My idea was pull_k1 = 1000 is the force to keep ligand in one place where it get sampled?<br></blockquote><div><br>I want to add some more points to this. I agree that as I am using the "Position" , it should be displaced to the position of reference. But my worry is that, when i visualize this trajectory, within the first two or three frames, ligand get displaced. As per my u understanding, we need to sample well at the same position to get the PMF at that point.(Is it so?) in that case, the sudden drop from the initial position will cause any problem in the PMF?<br>
<br><br></div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><br>Is there anything wrong that I am doing?<br>
any suggestions will be appreciated.<br><br>-Aswathy<div><div></div><div class="h5"><br><br><br><div class="gmail_quote">On Tue, May 11, 2010 at 11:15 AM, Aswathy <span dir="ltr"><<a href="mailto:ammasachu@gmail.com" target="_blank">ammasachu@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Hi Chris, <br><br>Thank you very much for your detailed mail.<br><br>Now I have a doubt on this pull_init parameter. i read your previous posts regarding this, but still have a confusion. <br>
<br>My query is that for each configuration when I run umbrella sampling, will this pull_init value needs to change?(I suppose so, if its true how?)<br>
<br>When it should be negative and positive?<br><br><br>Could you please explain this. Thanks for your valuable time<br><br>Thanks & Regards,<br>Aswathy<div><div></div><div><br><br><br><div class="gmail_quote">
On Mon, May 10, 2010 at 9:55 PM, Chris Neale <span dir="ltr"><<a href="mailto:chris.neale@utoronto.ca" target="_blank">chris.neale@utoronto.ca</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Pick a small collection of backbone atoms near the center of your channel and use them as your reference group. Overcome the sign problem by optimal selection of pull options (see below). pull_pbcatom values should not be important<br>
if you select your groups as I suggest -- otherwise be sure to understand how they work.<br>
<br>
<br>
; COM PULLING<br>
pull = umbrella<br>
pull_geometry = position<br>
pull_dim = N N Y<br>
pull_start = no<br>
pull_nstxout = 250<br>
pull_nstfout = 250<br>
pull_ngroups = 1<br>
pull_group0 = MY_SELECTION_OF_BACKBONE_ATOMS<br>
pull_pbcatom0 = 0<br>
pull_group1 = LIGAND<br>
pull_pbcatom1 = 0<br>
pull_init1 = 0 0 THISDIST<br>
pull_rate1 = 0<br>
pull_k1 = 500.0<br>
pull_vec1 = 0 0 0<br>
<br>
Chris.<br>
<br>
-- original message --<br>
<br>
I think Justin meant that you have various positions of the ligand in the channel (from the SMD), so you don't need to make a new run to determine new positions in the channel. You need only new umbrella sampling simulations.<br>
<br>
Yep, the movement of the particle will also matter, because if the particle moves much on the z-axis, the distance between the particle and the ligand will change. So you would want the particle fixed relative to the channel.<br>
Two ideas:<br>
Place the particle above the entrance of the channel. Pick three atoms from the entrance of the channel and determine the distance between the atoms and the particle. Then use distance_restraints or constraints with a 'bondlength' equal to the measured distance. If everything goes right the particle would stay where you placed it, since it does not interact with the enviroment it should not really influence your simulation (only through the constraints or distance_restraints).<br>
I don't now how big your system is, but it would probably be a good idea to make a short test simulation, to look if the particle changes the system behavior.<br>
But it's only an idea, i hope other people comment it.<br>
<br>
Greetings<br>
Thomas<br><font color="#888888">
<br>
<br>
-- <br></font><div><div></div><div>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before posting!<br>
Please don't post (un)subscribe requests to the list. Use the www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can't post? Read <a href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>
</div></div></blockquote></div><br><br clear="all"><br></div></div>-- <br><font color="#888888">Aswathy<br>
</font></blockquote></div><br><br clear="all"><br></div></div>-- <br><font color="#888888">Aswathy<br>
</font></blockquote></div><br><br clear="all"><br>-- <br>Aswathy<br>