<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div><br></div>Well I do not know how g_wham is doing but it is likely that it does align <div>different pmf obtained through the bootstraap the their minimum.</div><div><br></div><div>Justin just answered. The other alternative is to inverse your distances </div><div>with the data files.</div><div><br></div><div><div><div><div>On Jul 21, 2011, at 2:57 PM, Rebeca García Fandiño wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; font-size: medium; "><div class="hmmessage" style="font-size: 10pt; font-family: Tahoma; "><div dir="ltr">Thanks for your answer.<br>My peptide has 6 residues and it is cyclic.<br>I calculated the errors through bootstrap, using this:<br><br>g_wham -it tpr-files.dat -if pullf-files.dat -o -hist -unit kCal -b 500 -n Bootstrap 200<br><br>Are the pmfs aligned on the minimum distance value? Is that the value by default? How could I calculate the errors based on the PMFs aligned on the long distance value?<br><br>Thanks a lot for all the help.<br><br>Best wishes,<br><br>Rebeca.<br><br><br><div>> From:<span class="Apple-converted-space"> </span><a href="mailto:x.periole@rug.nl">x.periole@rug.nl</a><br>> To:<span class="Apple-converted-space"> </span><a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>> Subject: Re: [gmx-users] large error bars in PMF<br>> Date: Thu, 21 Jul 2011 14:33:28 -0600<br>><span class="Apple-converted-space"> </span><br>><span class="Apple-converted-space"> </span><br>> One potential problem you have is that as Justin mentioned your minimum<br>> is not well defined and certain much less well sampled than the long<span class="Apple-converted-space"> </span><br>> distances<br>> windows. Small peptides (depends the size) may sample relevant phase<br>> space to get reasonable convergence within 8 ns when free in solution;<br>> in contact certainly not ...<br>><span class="Apple-converted-space"> </span><br>> How long is you peptide?<br>><span class="Apple-converted-space"> </span><br>> You might get a better estimate of your error by plotting all the pmfs<span class="Apple-converted-space"> </span><br>> you've<br>> got through bootstrap (if this is the case) ... and get the errors<span class="Apple-converted-space"> </span><br>> based on the<br>> pmfs aligned on the long distance value and not on the "minimum" which<span class="Apple-converted-space"> </span><br>> is<br>> probably what you did ... if the minimum is not well defined then it<span class="Apple-converted-space"> </span><br>> does not<br>> make sense.<br>><span class="Apple-converted-space"> </span><br>> On Jul 21, 2011, at 1:36 PM, Justin A. Lemkul wrote:<br>><span class="Apple-converted-space"> </span><br>> ><br>> ><br>> > Rebeca García Fandiño wrote:<br>> >> OK,<br>> >> I will try a dodecahedral box and also to extend my actual<span class="Apple-converted-space"> </span><br>> >> simulations.<br>> >> Could you give me some advice about starting to learn about 3D PMF?<span class="Apple-converted-space"> </span><br>> >> I have not seen this in the manual, and I have never used it<span class="Apple-converted-space"> </span><br>> >> before. I have only found your tutorial about how to calculate PMF<span class="Apple-converted-space"> </span><br>> >> in Gromacs 4...<br>> ><br>> > In practice, there's basically nothing different between 3D and 1D.<span class="Apple-converted-space"> </span><br>> > You specify the dimensions to which the biasing potential is applied<span class="Apple-converted-space"> </span><br>> > with pull_dim or pull_vec (depending on the pull_geometry used).<span class="Apple-converted-space"> </span><br>> > Right now you're calculating the PMF along the z-dimension only<span class="Apple-converted-space"> </span><br>> > ("pull_dim = N N Y"), which may be appropriate for one-dimensional<span class="Apple-converted-space"> </span><br>> > processes or those in which the reference group does not rotate much<span class="Apple-converted-space"> </span><br>> > in the timeframe of the sampling. Setting "pull_dim = Y Y Y" will<span class="Apple-converted-space"> </span><br>> > apply the restraint in all dimensions.<br>> ><br>> > -Justin<br>> ><br>> >> Thanks a lot again for your help.<br>> >> Best wishes,<br>> >> Rebeca.<br>> >> > Date: Thu, 21 Jul 2011 15:16:52 -0400<br>> >> > From:<span class="Apple-converted-space"> </span><a href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a><br>> >> > To:<span class="Apple-converted-space"> </span><a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>> >> > Subject: Re: FW: [gmx-users] large error bars in PMF<br>> >> ><br>> >> ><br>> >> ><br>> >> > Rebeca García Fandiño wrote:<br>> >> > ><br>> >> > ><br>> >> > > I am trying to achieve the binding energy of the dimer composed<span class="Apple-converted-space"> </span><br>> >> by the<br>> >> > > two small cyclic peptides, to compare it with experimental. What<br>> >> > > advantages would I have using 3D PMF instead only 1D for this<span class="Apple-converted-space"> </span><br>> >> calculation?<br>> >> ><br>> >> > Intuitively, two molecules diffuse through solution until they<span class="Apple-converted-space"> </span><br>> >> find one another,<br>> >> > which to me sounds a lot like a 3D path. Further, using a<span class="Apple-converted-space"> </span><br>> >> dodecahedral box for<br>> >> > your umbrella sampling removes the problems you're having with<span class="Apple-converted-space"> </span><br>> >> the peptides<br>> >> > rotating. It sounds like you're trying to pull in one direction<span class="Apple-converted-space"> </span><br>> >> along a<br>> >> > rectangular box, but the peptides are not playing nice. I feel<span class="Apple-converted-space"> </span><br>> >> like this<br>> >> > discussion has come up at least once or twice before, though...<br>> >> ><br>> >> > -Justin<br>> >> ><br>> >> > > Thanks a lot!<br>> >> > > Rebeca.<br>> >> > ><br>> >> > > > Date: Thu, 21 Jul 2011 14:14:44 -0400<br>> >> > > > From:<span class="Apple-converted-space"> </span><a href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a><br>> >> > > > To:<span class="Apple-converted-space"> </span><a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>> >> > > > Subject: Re: [gmx-users] large error bars in PMF<br>> >> > > ><br>> >> > > ><br>> >> > > ><br>> >> > > > Rebeca García Fandiño wrote:<br>> >> > > > > Hi,<br>> >> > > > > thanks a lot for your quick answer.<br>> >> > > > > What I am trying to pull are two small peptides one from<span class="Apple-converted-space"> </span><br>> >> another (r_1<br>> >> > > > > and r_2).<br>> >> > > > > I did not understand very well your last suggestion: "...if<span class="Apple-converted-space"> </span><br>> >> you want<br>> >> > > > > reasonable error bars you will not lots of well-converged<span class="Apple-converted-space"> </span><br>> >> data".<br>> >> > > ><br>> >> > > > Oops, that should have read "you will *need* lots of well-<span class="Apple-converted-space"> </span><br>> >> converged<br>> >> > > data."<br>> >> > > ><br>> >> > > > > Do you mean I will need also more windows besides extending<span class="Apple-converted-space"> </span><br>> >> the<br>> >> > > simulations?<br>> >> > > ><br>> >> > > > I doubt you need more windows. Likely you just need more time<span class="Apple-converted-space"> </span><br>> >> in each.<br>> >> > > ><br>> >> > > > > I think the problem could be also that the peptides I am<span class="Apple-converted-space"> </span><br>> >> using<br>> >> > > rotate in<br>> >> > > > > the box and they do not remain flat one respect to the<span class="Apple-converted-space"> </span><br>> >> other. They<br>> >> > > > > gyrate freely and some parts of their structure interact<span class="Apple-converted-space"> </span><br>> >> along the<br>> >> > > > > pulling...<br>> >> > > ><br>> >> > > > Interactions are part of the dissociation process and are not<br>> >> > > problematic per<br>> >> > > > se. But if you're trying to obtain only a one-dimensional PMF<span class="Apple-converted-space"> </span><br>> >> then your<br>> >> > > > rotation could be a problem. Is there some reason you need a<br>> >> > > one-dimensional<br>> >> > > > PMF and not a three-dimensional PMF? What are you trying to<span class="Apple-converted-space"> </span><br>> >> achieve?<br>> >> > > ><br>> >> > > > -Justin<br>> >> > > ><br>> >> > > > > Thanks a lot again for your help.<br>> >> > > > > Best wishes,<br>> >> > > > > Rebeca.<br>> >> > > > ><br>> >> > > > ><br>> >> > > > ><br>> >> > > > ><br>> >> > ><span class="Apple-converted-space"> </span><br>> >> ------------------------------------------------------------------------<br>> >> > > > > From:<span class="Apple-converted-space"> </span><a href="mailto:regafan@hotmail.com">regafan@hotmail.com</a><br>> >> > > > > To:<span class="Apple-converted-space"> </span><a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>> >> > > > > Date: Thu, 21 Jul 2011 16:36:59 +0000<br>> >> > > > > Subject: [gmx-users] large error bars in PMF<br>> >> > > > ><br>> >> > > > ><br>> >> > > > > Hi,<br>> >> > > > > I am trying to calculate the binding energy of two<span class="Apple-converted-space"> </span><br>> >> molecules using the<br>> >> > > > > PMF (Umbrella Sampling method) and Gromacs 4.0.<br>> >> > > > > Some weeks ago I have written to the list because changing<span class="Apple-converted-space"> </span><br>> >> the<br>> >> > > number of<br>> >> > > > > windows used in the Umbrella Sampling calculations<span class="Apple-converted-space"> </span><br>> >> different results<br>> >> > > > > were obtained, and I was suggested to extend my simulations<span class="Apple-converted-space"> </span><br>> >> since the<br>> >> > > > > error bars associated to each windows were too high.<br>> >> > > > > I have now extended my simulations from 1 ns to 8 ns,<span class="Apple-converted-space"> </span><br>> >> however, I<br>> >> > > cannot<br>> >> > > > > see much different from the shorter calculations. I send<span class="Apple-converted-space"> </span><br>> >> you the<br>> >> > > > > comparison of the two PMF including the error bars<span class="Apple-converted-space"> </span><br>> >> (attached).<br>> >> > > > > Now I am using 50 windows, but the shorter simulations were<span class="Apple-converted-space"> </span><br>> >> done using<br>> >> > > > > 100 windows, so I don't think increasing the number of<span class="Apple-converted-space"> </span><br>> >> windows<br>> >> > > could help.<br>> >> > > > > My system has about 29200 atoms (where 29000 are chloroform<span class="Apple-converted-space"> </span><br>> >> atoms).<br>> >> > > The<br>> >> > > > > *mdp file I am using is copied below.<br>> >> > > > > Would you have any suggestion to improve the results and<span class="Apple-converted-space"> </span><br>> >> decrease the<br>> >> > > > > error bars in the calculations?<br>> >> > > > ><br>> >> > > > > ----------------------------MDP<span class="Apple-converted-space"> </span><br>> >> file---------------------------<br>> >> > > > > title = Umbrella pulling simulation<br>> >> > > > > define =<br>> >> > > > > define =<br>> >> > > > > ; Run parameters<br>> >> > > > > integrator = md<br>> >> > > > > dt = 0.002<br>> >> > > > > tinit = 0<br>> >> > > > > nsteps = 500000 ; 1 ns<br>> >> > > > > nstcomm = 10<br>> >> > > > > ; Output parameters<br>> >> > > > > nstxout = 5000 ; every 10 ps<br>> >> > > > > nstvout = 5000<br>> >> > > > > nstfout = 5000<br>> >> > > > > nstxtcout = 5000 ; every 10 ps<br>> >> > > > > nstenergy = 5000<br>> >> > > > > ; Bond parameters<br>> >> > > > > constraint_algorithm = lincs<br>> >> > > > > constraints = all-bonds<br>> >> > > > > continuation = yes<br>> >> > > > > ; Single-range cutoff scheme<br>> >> > > > > nstlist = 5<br>> >> > > > > ns_type = grid<br>> >> > > > > rlist = 1.4<br>> >> > > > > rcoulomb = 1.4<br>> >> > > > > rvdw = 1.4<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>> >> > > > > ; Berendsen temperature coupling is on in two groups<br>> >> > > > > Tcoupl = Nose-Hoover<br>> >> > > > > tc_grps = ACH CL3<br>> >> > > > > tau_t = 0.5 0.5<br>> >> > > > > ref_t = 300 300<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<br>> >> > > > > pull_dim = N N Y<br>> >> > > > > pull_start = yes<br>> >> > > > > pull_ngroups = 1<br>> >> > > > > pull_group0 = r_1<br>> >> > > > > pull_group1 = r_2<br>> >> > > > > pull_init1 = 0<br>> >> > > > > pull_rate1 = 0.0<br>> >> > > > > pull_k1 = 1000 ; kJ mol^-1 nm^-2<br>> >> > > > > pull_nstxout = 1000 ; every 2 ps<br>> >> > > > > pull_nstfout = 1000 ; every 2 ps<br>> >> > > > ><br>> >> > > > > -----------------------------------------------<br>> >> > > > ><br>> >> > > > ><br>> >> > > > > Thanks a lot in advance.<br>> >> > > > ><br>> >> > > > > Best wishes,<br>> >> > > > ><br>> >> > > > > Dr. Rebeca Garcia<br>> >> > > > > Santiago de Compostela University<br>> >> > > > > Spain<br>> >> > > > ><br>> >> > > > ><br>> >> > > > ><br>> >> > > > > -- gmx-users mailing list<span class="Apple-converted-space"> </span><a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>> >> > > > ><span class="Apple-converted-space"> </span><a href="http://lists.gromacs.org/mailman/listinfo/gmx-users">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><span class="Apple-converted-space"> </span>Please<span class="Apple-converted-space"> </span><br>> >> search the<br>> >> > > > > archive at<span class="Apple-converted-space"> </span><a href="http://www.gromacs.org/Support/Mailing_Lists/Search">http://www.gromacs.org/Support/Mailing_Lists/Search</a><span class="Apple-converted-space"> </span><br>> >> before<br>> >> > > > > posting! Please don't post (un)subscribe requests to the<span class="Apple-converted-space"> </span><br>> >> list. Use the<br>> >> > > > > www interface or send it to<span class="Apple-converted-space"> </span><a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<span class="Apple-converted-space"> </span><br>> >> Can't post?<br>> >> > > > > Read<span class="Apple-converted-space"> </span><a href="http://www.gromacs.org/Support/Mailing_Lists">http://www.gromacs.org/Support/Mailing_Lists</a><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>> >> > > ><span class="Apple-converted-space"> </span><a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br>> >> > > ><br>> >> > > > ========================================<br>> >> > > > --<br>> >> > > > gmx-users mailing list<span class="Apple-converted-space"> </span><a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>> >> > > ><span class="Apple-converted-space"> </span><a href="http://lists.gromacs.org/mailman/listinfo/gmx-users">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>> >> > > > Please search the archive at<br>> >> > ><span class="Apple-converted-space"> </span><a href="http://www.gromacs.org/Support/Mailing_Lists/Search">http://www.gromacs.org/Support/Mailing_Lists/Search</a><span class="Apple-converted-space"> </span>before<span class="Apple-converted-space"> </span><br>> >> posting!<br>> >> > > > Please don't post (un)subscribe requests to the list. Use the<br>> >> > > > www interface or send it to<span class="Apple-converted-space"> </span><a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>> >> > > > Can't post? Read<span class="Apple-converted-space"> </span><a href="http://www.gromacs.org/Support/Mailing_Lists">http://www.gromacs.org/Support/Mailing_Lists</a><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>> >> ><span class="Apple-converted-space"> </span><a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br>> >> ><br>> >> > ========================================<br>> >> > --<br>> >> > gmx-users mailing list<span class="Apple-converted-space"> </span><a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>> >> ><span class="Apple-converted-space"> </span><a href="http://lists.gromacs.org/mailman/listinfo/gmx-users">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>> >> > Please search the archive at<span class="Apple-converted-space"> </span><a href="http://www.gromacs.org/Support/Mailing_Lists/Search">http://www.gromacs.org/Support/Mailing_Lists/Search</a><span class="Apple-converted-space"> </span><br>> >> before posting!<br>> >> > Please don't post (un)subscribe requests to the list. Use the<br>> >> > www interface or send it to<span class="Apple-converted-space"> </span><a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>> >> > Can't post? Read<span class="Apple-converted-space"> </span><a href="http://www.gromacs.org/Support/Mailing_Lists">http://www.gromacs.org/Support/Mailing_Lists</a><br>> ><br>> > --<span class="Apple-converted-space"> </span><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>> ><span class="Apple-converted-space"> </span><a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br>> ><br>> > ========================================<br>> > --<span class="Apple-converted-space"> </span><br>> > gmx-users mailing list<span class="Apple-converted-space"> </span><a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>> ><span class="Apple-converted-space"> </span><a href="http://lists.gromacs.org/mailman/listinfo/gmx-users">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>> > Please search the archive at<span class="Apple-converted-space"> </span><a href="http://www.gromacs.org/Support/Mailing_Lists/Search">http://www.gromacs.org/Support/Mailing_Lists/Search</a><span class="Apple-converted-space"> </span><br>> > before posting!<br>> > Please don't post (un)subscribe requests to the list. Use the www<span class="Apple-converted-space"> </span><br>> > interface or send it to<span class="Apple-converted-space"> </span><a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>> > Can't post? Read<span class="Apple-converted-space"> </span><a href="http://www.gromacs.org/Support/Mailing_Lists">http://www.gromacs.org/Support/Mailing_Lists</a><br>><span class="Apple-converted-space"> </span><br>> --<span class="Apple-converted-space"> </span><br>> gmx-users mailing list<span class="Apple-converted-space"> </span><a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>><span class="Apple-converted-space"> </span><a href="http://lists.gromacs.org/mailman/listinfo/gmx-users">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>> Please search the archive at<span class="Apple-converted-space"> </span><a href="http://www.gromacs.org/Support/Mailing_Lists/Search">http://www.gromacs.org/Support/Mailing_Lists/Search</a><span class="Apple-converted-space"> </span>before posting!<br>> Please don't post (un)subscribe requests to the list. Use the<span class="Apple-converted-space"> </span><br>> www interface or send it to<span class="Apple-converted-space"> </span><a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>> Can't post? Read<span class="Apple-converted-space"> </span><a href="http://www.gromacs.org/Support/Mailing_Lists">http://www.gromacs.org/Support/Mailing_Lists</a><br></div></div>--<span class="Apple-converted-space"> </span><br>gmx-users mailing list <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br><a href="http://lists.gromacs.org/mailman/listinfo/gmx-users">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>Please search the archive at<span class="Apple-converted-space"> </span><a href="http://www.gromacs.org/Support/Mailing_Lists/Search">http://www.gromacs.org/Support/Mailing_Lists/Search</a><span class="Apple-converted-space"> </span>before posting!<br>Please don't post (un)subscribe requests to the list. Use the<span class="Apple-converted-space"> </span><br>www interface or send it to<span class="Apple-converted-space"> </span><a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>Can't post? Read<span class="Apple-converted-space"> </span><a href="http://www.gromacs.org/Support/Mailing_Lists">http://www.gromacs.org/Support/Mailing_Lists</a></div></span></blockquote></div><br></div></div></body></html>