Would you recommend position restrained of the protein backbone atoms or e.g. residues from the lower part so that the protien will not roatate?<br><br><div class="gmail_quote">On Fri, Feb 17, 2012 at 10:41 AM, Steven Neumann <span dir="ltr"><<a href="mailto:s.neumann08@gmail.com">s.neumann08@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello Justin,<br><br>As you recommended I run longer pulling of my ligand. I pull the ligand which is on the top of my protein so the top of the protein is pulled so that the whole protein rotated app. 30 degrees - the low part of the protein came out of the box due to the rotation. After app 600 ps the ligand does not move any more and the force applied increase linearly. As I saw the trajectory, ligand does not collide with a periodic image but it does not move. Please, see attacheg plot force vs time. Can you explain why the force increase linearly and ligand is not pulled any more? This is my mdp file for pulling:<br>
<br>title = Umbrella pulling simulation <br>; Run parameters<br>integrator = md<br>dt = 0.002<br>tinit = 0<br>nsteps = 500000 ; 500 ps<br>nstcomm = 10<br>; Output parameters<br>nstxout = 5000 ; every 10 ps<br>
nstvout = 5000 <br>nstfout = 500<br>nstxtcout = 500 ; every 1 ps<br>nstenergy = 500<br>; Bond parameters<br>constraint_algorithm = lincs<br>constraints = all-bonds<br>continuation = yes ; continuing from NPT <br>
; Single-range cutoff scheme<br>nstlist = 5<br>ns_type = grid <br>rlist = 0.9<br>rcoulomb = 0.9<br>rvdw = 0.9<br>; PME electrostatics parameters<br>coulombtype = PME<br>fourierspacing = 0.12<br>
fourier_nx = 0<br>fourier_ny = 0<br>fourier_nz = 0<br>pme_order = 4<br>ewald_rtol = 1e-5<br>optimize_fft = yes<br>; Temperature coupling is on<br>tcoupl = V-rescale ; modified Berendsen thermostat<br>
tc_grps = Protein_LIG Water_and_ions ; two coupling groups - more accurate<br>tau_t = 0.1 0.1 ; time constant, in ps<br>ref_t = 298 298 ; reference temperature, one for each group, in K<br>
; Pressure coupling is on<br>Pcoupl = Parrinello-Rahman <br>pcoupltype = isotropic<br>tau_p = 1.0 <br>compressibility = 4.5e-5<br>ref_p = 1.0<br>; Generate velocities is off<br>gen_vel = no <br>
; Periodic boundary conditions are on in all directions<br>pbc = xyz<br>; Long-range dispersion correction<br>DispCorr = EnerPres<br>; Pull code<br>pull = umbrella<br>pull_geometry = distance ; simple distance increase <br>
pull_dim = N N Y<br>pull_start = yes ; define initial COM distance > 0<br>pull_ngroups = 1<br>pull_group0 = Protein <br>pull_group1 = LIG182 <br>pull_rate1 = 0.01 ; 0.01 nm per ps = 10 nm per ns<br>
pull_k1 = 200 ; kJ mol^-1 nm^-2<br><br>Thank you,<br><br>Steven<div class="HOEnZb"><div class="h5"><br><br><br><div class="gmail_quote">On Mon, Feb 13, 2012 at 2:53 PM, Steven Neumann <span dir="ltr"><<a href="mailto:s.neumann08@gmail.com" target="_blank">s.neumann08@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Thank you Justin!<div><div><br><br>
<div class="gmail_quote">On Mon, Feb 13, 2012 at 2:44 PM, Justin A. Lemkul <span dir="ltr"><<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>></span> wrote:<br>
<blockquote style="BORDER-LEFT:#ccc 1px solid;MARGIN:0px 0px 0px 0.8ex;PADDING-LEFT:1ex" class="gmail_quote">
<div><br><br>Steven Neumann wrote:<br>
<blockquote style="BORDER-LEFT:#ccc 1px solid;MARGIN:0px 0px 0px 0.8ex;PADDING-LEFT:1ex" class="gmail_quote">Thank you Justin. I run my pulling using two force constant for pulling. K1=100 and K1=200<br>Please, see attached plots of force vs time. Is there any criteria to adjust pulling constant? Would you suggest running it for a longer time?<br>
</blockquote><br></div>I know of no systematic study for choosing a force constant. Guessing wildly at what's going on, I'd say you need longer simulations as it appears you have only just caused dissociation towards the end of the 500 ps.<br>
<br>-Justin<br><br>
<blockquote style="BORDER-LEFT:#ccc 1px solid;MARGIN:0px 0px 0px 0.8ex;PADDING-LEFT:1ex" class="gmail_quote">Steven
<div><br>On Mon, Feb 13, 2012 at 1:11 PM, Justin A. Lemkul <<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a> <mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>>> wrote:<br>
<br><br><br> Steven Neumann wrote:<br><br> Dear Gmx Users,<br> Is it always required to restrained positions of the protein<br> while pulling your ligand? My system is made of 10 ligands<br> attached to my protein surface. I am pulling one of them.<br>
<br><br> No, it is not required. I assume you've gotten this idea from my<br> tutorial - the restraints there were used for a very specific<br> purpose (detailed in the paper linked from the tutorial).<br><br>
<br>
I have just seen trajectory of pulling my ligand without<br> restraining positions of protein and 9 remaining ligands. My<br> ligand while pulling also pulled the protein with itself (for 1<br> nm distance) and then splited. Is is this approach more reliable?<br>
<br><br> If your goal is umbrella sampling, you need only generate a series<br> of reasonable starting configurations along a defined reaction<br> coordinate. The absolute positions are irrelevant; it is the<br>
relative distance that matters.<br><br> -Justin<br><br></div> -- ==============================<u></u>__==========
<div><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></div> jalemkul[at]<a href="http://vt.edu/" target="_blank">vt.edu</a> <<a href="http://vt.edu/" target="_blank">http://vt.edu/</a>> | <a href="tel:%28540%29%20231-9080" value="+15402319080" target="_blank">(540) 231-9080</a><br>
<tel:%28540%29%20231-9080><br> <a href="http://www.bevanlab.biochem./" target="_blank">http://www.bevanlab.biochem.</a>__<a href="http://vt.edu/Pages/Personal/justin" target="_blank"><u></u>vt.edu/Pages/Personal/justin</a><br>
<<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.<u></u>vt.edu/Pages/Personal/justin</a>><br><br> ==============================<u></u>__==========
<div><br> -- gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br></div> <mailto:<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/__<u></u>mailman/listinfo/gmx-users</a>
<div><br> <<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/<u></u>mailman/listinfo/gmx-users</a>><br> Please search the archive at<br></div> <a href="http://www.gromacs.org/__Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/__<u></u>Support/Mailing_Lists/Search</a>
<div><br> <<a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists/Search</a>> before posting!<br> Please don't post (un)subscribe requests to the list. Use the www<br>
interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br></div> <mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@<u></u>gromacs.org</a>>.<br>
Can't post? Read <a href="http://www.gromacs.org/__Support/Mailing_Lists" target="_blank">http://www.gromacs.org/__<u></u>Support/Mailing_Lists</a><br> <<a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists</a>><br>
<br><br><br>------------------------------<u></u>------------------------------<u></u>------------<br><br><br>------------------------------<u></u>------------------------------<u></u>------------<br><br></blockquote>
<div>
<div><br>-- <br>==============================<u></u>==========<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]<a href="http://vt.edu/" target="_blank">vt.edu</a> | <a href="tel:%28540%29%20231-9080" value="+15402319080" target="_blank">(540) 231-9080</a><br><a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.<u></u>vt.edu/Pages/Personal/justin</a><br>
<br>==============================<u></u>==========<br>-- <br>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/<u></u>mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/<u></u>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 <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/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists</a><br></div></div></blockquote></div><br>
</div></div></blockquote></div><br>
</div></div></blockquote></div><br>