Re: [AMBER] restrained minimization

From: David A. Case <case.biomaps.rutgers.edu>
Date: Tue, 14 Apr 2009 03:03:11 +0100

On Mon, Apr 13, 2009, quantum_mania.yahoo.com wrote:
>
> I wanted to minimize a system with all heavy atoms frozen, so I had the input file:
>
> &cntrl
> ntxo = 1,
> ntf = 1, ntb = 0,
> cut = $13,
> imin = 1, maxcyc = 500,
> ncyc = 200, drms = 0.0001,
> igb = 2, gbsa = 1,
> saltcon= 0.2,
> ntr = 1, restraint_wt=100,
> restraintmask="!.H="
> &end

You can use the ambmask program to check out masks to make sure they are
doing what you want.

>
> The output, though, showed most atoms minimized were the heavy atoms:

?? I don't see how your output shows what you claim. Have you examined
the final coordinates themselves to see what moved? Note that you have
very bad contacts in the initial structure (very large VDW terms), so
you should expect some motion from heavy atoms, even with a very large
restraint weight.

...dac


_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed Apr 15 2009 - 01:12:30 PDT
Custom Search