AMBER Archive (2008)

Subject: AMBER: problem with restraints

From: Barbault Florent (florent.barbault_at_paris7.jussieu.fr)
Date: Thu Jul 10 2008 - 13:59:15 CDT

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_at_scripps.edu
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
      to majordomo_at_scripps.edu