[gmx-users] fail to pull
Albert
mailmd2011 at gmail.com
Sat Apr 6 20:47:56 CEST 2013
Dear:
I am trying to pull my ligand outside of the binding pocket with
following configurations:
title = Umbrella pulling simulation
define = -DPOSRES
; Pull code
pull = umbrella
pull_geometry = distance ; simple distance increase
pull_dim = Y N N
pull_start = yes ; define initial COM distance > 0
pull_ngroups = 1
pull_group0 = Protein
pull_group1 = LIG
pull_rate1 = 0.001 ; 0.01 nm per ps = 10 nm per ns
pull_k1 = 1000 ; kJ mol^-1 nm^-2
Tcoupl = v-rescale
tc_grps = Protein_LIG Water_and_ions
tau_t = 0.5 0.5
ref_t = 310 310
; Pressure coupling is on
Pcoupl = Parrinello-Rahman
pcoupltype = isotropic
tau_p = 1.0 1.0
compressibility = 4.5e-5
ref_p = 1.0 1.0
refcoord_scaling = com
; Generate velocities is off
gen_vel = no
; Periodic boundary conditions are on in all directions
pbc = xyz
; Long-range dispersion correction
DispCorr = EnerPres
It is quite strange, the ligand is still in place and not outside the
pocket at the end of simulations. I am just wondering where is the problem?
thank you very much
best
Albert
More information about the gromacs.org_gmx-users
mailing list