[gmx-users] position restraints

[Steven Neumann]

In your topology:

; Include Protein topology
#include "protein.itp"

; Include Position restraint file of ligand
#ifdef POSRES
#include "posre_prot.itp"
#endif


where posre_prot.itp

; In this topology include file, you will find position restraint
; entries for all the heavy atoms in your original pdb file.
; This means that all the protons which were added by pdb2gmx are
; not restrained.

[ position_restraints ]
; atom  type      fx      fy      fz
1    1    1000    1000    1000
2...
...
X

here you specify atoms you wish to restraints from topology of your
molecule in given direction X, Y, Z. The higher the fx, fy, fz the more
stable it will be.
Then in mdp:

define = -DPOSRES

Steven



On Wed, Feb 5, 2014 at 1:38 PM, Gurunath Katagi
<gurunath.katagi at gmail.com>wrote:

> Hi everyone,
> I am doing a simulation of a protein using gromacs with charmm forcefields.
> From the crystallographic analysis, it is found that certain residues are
> rigid and few are flexible. I want to apply restraints on these rigid
> residues.
> Can anyone let me know how can I apply restrains on only these rigid
> residues, so that do dont move much during simulation and leave other
> residues without any restraints..
>
> Thank you
> Gurunath
> --
> Gromacs Users mailing list
>
> * Please search the archive at
> http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
> posting!
>
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>
> * For (un)subscribe requests visit
> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
> send a mail to gmx-users-request at gromacs.org.
>