[gmx-users] few doubts regarding 1-4 interaction

[Ms. Aswathy S]

Thanks Justin.

I am trying to do the simulation of protein-ligand complex.

1. Minimized the protein(but converged before  the Fmax reached the threshold)
2. Did a 50 steps of NVT equilibration
3. Then an NPT equilibration reached at 650 steps (checked the RMSD with the initial structure) 
4. Then tried to do a 3 ns production run. But failed  ~ after 170 ps.

I am attching the NPT file I have used for equilibration. Also the structure after equilibration does not have any clashes or bad contact in the rough analysis. How can I do a more close checking of this?

Once again Thank you  for your response.
Aswathy

----- Original Message -----
From: "Justin A. Lemkul" <jalemkul at vt.edu>
To: "Discussion list for GROMACS users" <gmx-users at gromacs.org>
Sent: Tuesday, June 30, 2009 8:33:19 PM GMT +05:30 Chennai, Kolkata, Mumbai, New Delhi
Subject: Re: [gmx-users] few doubts regarding 1-4 interaction



Ms. Aswathy S wrote:
> Hi gromacs users,
> 
> I am facing the same error during Production run that faced by so many
> gromacs users and I know this query has been answered for several times. But
> Could you please tell me some specific doubt regarding the 1-4 interaction
> error that usually occurs?
> 
> I have equilibrated my protein ligand system of around 600 ps and shown a
> satisfied level of equilibration. Then I moved to the Production run of 3ns.
> But shows the same 1-4 interaction error. My doubts are,
> 
> 1. It was advised by gromacs experts that better I should go for another
> round of minimization. Is it on the equilibrated structure?

It would be good if you could remind us of what you are doing, as well as cite
the actual post where you were told to do this.

> 2. If so Do i need to equilibrate again after minimization?

Yes.

> 3. I have did 2500 steps of minimization on the structure which is
> equilibrated and then again run for 3ns production run But occured the same
> error at 3 lakh step? if I reduce my time step how that will affect this
> error(I have seen in the error section of gromacs that another possibility is
> to reduce the time step?)

Reducing the timestep can in some cases finesse the problem into behaving.  If
there is something more deeply wrong with your model, then this probably won't
help out a whole lot.

If you can post a more complete description of what it is you are doing, what
parameters you are using, etc. you might get more detailed suggestions.  The
bottom line is, if your system is blowing up, something about the model physics
is unreasonable.

-Justin

> Can you please tell me your suggestion on these??
> 
> Thank you. Aswathy
> 
> Dept. Biotechnology Ext. 3108 _______________________________________________
>  gmx-users mailing list    gmx-users at gromacs.org 
> http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive
> at http://www.gromacs.org/search before posting! Please don't post
> (un)subscribe requests to the list. Use the www interface or send it to
> gmx-users-request at gromacs.org. Can't post? Read
> http://www.gromacs.org/mailing_lists/users.php
> 

-- 
========================================

Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin

========================================
_______________________________________________
gmx-users mailing list    gmx-users at gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-request at gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
-------------- next part --------------
A non-text attachment was scrubbed...
Name: md_NPT.mdp
Type: application/octet-stream
Size: 1073 bytes
Desc: not available
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20090701/c2d03aeb/attachment.obj>