[gmx-users] Position Restraints MD

ISHRAT JAHAN jishrat17 at gmail.com
Tue Sep 24 06:05:15 CEST 2019


Thanks Justin Sir for the response.

On Mon, Sep 23, 2019 at 4:18 PM Justin Lemkul <jalemkul at vt.edu> wrote:

>
>
> On 9/20/19 2:10 AM, ISHRAT JAHAN wrote:
> >   Dear all,
> > I am doing MD simulation of protein and metal atom. I had applied
> position
> > restraints on metal atom. I had applied it as follows- ( kindly correct
> me
> > if i am wrong)
> > first made an index of metal atom and generated the posre_ion.itp file
> > using gmx genrestr and then added the force and atom index into the top
> > file of protein in [position restraint] part. In full.mdp file , added
> the
> > line define= -DPOSRES.
> > Then performed full md. But when i calculated the rmsd of backbone its
> > value is found to be around 0.02 nm(very small). I do not know whether
> > restraints have been applied on metal or on the whole system.  Am I doing
> > it correct or not?
>
> Using -DPOSRES is the default keyword for enabling position restraints
> on the entire protein. Given your tiny RMSD, that should confirm it.
> Visualize your trajectory and you should see that the protein isn't moving.
>
> Applying an absolute restraint to a metal ion is a curious approach.
> What may end up happening is your protein will diffuse but the metal
> will stay behind. If you're trying to enforce some type of coordination,
> use distance restraints instead.
>
> -Justin
>
> --
> ==================================================
>
> Justin A. Lemkul, Ph.D.
> Assistant Professor
> Office: 301 Fralin Hall
> Lab: 303 Engel Hall
>
> Virginia Tech Department of Biochemistry
> 340 West Campus Dr.
> Blacksburg, VA 24061
>
> jalemkul at vt.edu | (540) 231-3129
> http://www.thelemkullab.com
>
> ==================================================
>
> --
> 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.
>


-- 
Ishrat Jahan
Research Scholar
Department Of Chemistry
A.M.U Aligarh


More information about the gromacs.org_gmx-users mailing list