[gmx-users] rotation while pulling
gromacs query
gromacsquery at gmail.com
Tue Mar 28 13:37:16 CEST 2017
Hi All,
{apologies if posted multiple times; seems like some trouble with my
account}
I want to pull linear molecule across the membrane in two ways allowing:
- translation along Z and rotation around Z
- translation along Z and rotation in all axes
My current mdp-setting pulls molecular along Z but it translates and
rotates in all axes. Do I need to use some advanced options?
; Pull Code
pull = umbrella
pull_geometry = direction
pull-vec1 = 0 0 -1
pull_dim = N N Y ;
pull_start = yes ; define initial COM distance > 0
pull_ngroups = 1
pull_group0 = MEM ; restrained group
pull_group1 = LIG; group to pull
pull_init1 = 0
pull_rate1 = 0.00001; nm per ps
pull_k1 = 1000 ; kJ mol^-1 nm^-2
Also in this case would it be wise to use/say absolute Z coordinate as
reaction coordinate? rather center of mass distances between membrane and
ligand.
Thanks
Jiom
More information about the gromacs.org_gmx-users
mailing list