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
