On 10/17/12 5:44 AM, Netaly Khazanov wrote:
Of course I read Justin's tutorial.
The reason I am asking is that I am not sure how to do it.


As far as I am aware, there is no way to conduct umbrella sampling in this way using Gromacs. The PMF calculations that Gromacs does are based on restrained distances or angles. I do not believe one can set an arbitrary variable like RMSD to be a reaction coordinate in the same sense, and certainly the pull code is not applicable in this instance.

-Justin

On Wed, Oct 17, 2012 at 11:09 AM, Arman Mahboubi Soufiani <
armansoufi...@gmail.com> wrote:

Hi,

I recommend you Justin's tutorial on Umbrella Sampling

http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin/gmx-tutorials/umbrella/index.html

Regards

Arman

On Wed, Oct 17, 2012 at 10:47 AM, Netaly Khazanov <neta...@gmail.com>
wrote:

Dear All,

I've performed TMD simulation using NAMDprogram. I'd like to calculate
PMF
plot based on the frames that were taken from TMD simulations by using
GROMACS.
The  coordinate  of reaction is RMSD of backbone.
The first step is to do Umbrella sampling. The question is how fix the
position of backbone during the Umbrella sampling?
Is this looks reasonable?

pull            = umbrella
pull_geometry   = position
;pull_dim        = Y Y Y
pull_start      = yes
pull_ngroups    = 1
;pull_group0     =0
pull_group1     = backbone
pull_init1      = 0
pull_rate1      = 0.0
pull_k1         = 1000      ; kJ mol^-1 nm^-2
pull_nstxout    = 1000      ; every 2 ps
pull_nstfout    = 1000      ; every 2 ps

Any ideas how to convert the reaction coordinate from position of
backbone
to RMSD of backbone in PMF plot?

Thank you in advance.

Netaly Khazanov
--
Netaly
--
gmx-users mailing list    gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

--
gmx-users mailing list    gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists





--
========================================

Justin A. Lemkul, Ph.D.
Research Scientist
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin

========================================
--
gmx-users mailing list    gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Reply via email to