AMBER: problem with restraints

From: Barbault Florent <florent.barbault.paris7.jussieu.fr>
Date: Thu, 10 Jul 2008 20:59:15 +0200

Hello,

I am facing a problem with AMBER 9.0. I am doing molecular dynamics
with explicit solvent of a ligand/protein structure. I use harmonic
restraints to fix a part of the protein. These restraints are
expressed with this mask: '(:111 >.8.0 & :1-100!.H=)'

The mask contains 718 atoms and should correspond to the heavy atoms
which are far from 8 A of the residue 111.

Generally, I launch 20 runs of 100 ps in order to get 2 ns trajectory.
However, after each restart (irest=1, ntx=7) the restraint energy is
increasing and then the system is instable. This doesn't happen if I
launch 2ns simulation in one time.

Since the reference structure is identical in every running steps, I
don't understand why I get this problem. Any help will be highly
appreciated.

Regards
Florent Barbault




-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
      to majordomo.scripps.edu
Received on Sun Jul 13 2008 - 06:07:33 PDT
Custom Search