[gmx-users] restraints definition in mdp

Sana Saeed bioinformatic.lady at yahoo.com
Wed May 4 10:13:59 CEST 2016


hi gromacs experts
i have applied distance, angle and dihedral restraints to my ligand in a protein's pocket so that it doesnt roam around when i turn off all interactions during simulation. do i need to mention about these restraints in the mdp file? i am using gromacs 5.0 and 5.1. i have used intermolecular_interactions directive for defining these restraints. 
MY SECOND QUESTION is that the force constant i used is 41.840 kj/mol/K for angle and dihedral restr. does that mean that the Force constant for distance restraint is also 41.840? as i am calculating free energy of restraints (part of absolute binding F.E) so i need force constant for distance also. 

[ intermolecular-interactions ]
[ distance_restraints ]
;   i     j type label      funct         lo        up1        up2     weight
  768  2131 1     0          1   0.640512   0.840512    1.84051          1

[ angle_restraints ]
;   ai    aj    ak    al  type    thA      fcA    multA  thB      fcB    multB
   770   768    2131  768  1      43.420    0.0   1      43.420   41.840    1
   768   2131   2125  2131 1     114.89    0.0   1      114.89   41.840    1

[ dihedral_restraints ]
;   ai    aj    ak    al  type    phiA     dphiA  fcA    phiB      dphiB  fcB
  782  770   768   2131  1       -152.840  0.0  0.0    -152.840    0.0  41.840
  770  768   2131  2125  1       -46.638  0.0  0.0     -46.638     0.0  41.840
  768  2131  2125  2132  1       -42.465   0.0  0.0    -42.465     0.0  41.840

 
Sana Saeed Khan,
Research Assistant
Chemoinformatics Lab
Graduate Student, MS bioinfo
Department of Bioinformatics
Soongsil University, Seoul, South Korea.


More information about the gromacs.org_gmx-users mailing list