From: Chenglong Li <>
Date: Wed, 19 Nov 2003 10:44:19 -0800

Dear AMBER users,

In MM_PB/GB/SA analysis, I have problems with both un-reduced and
reduced complex approaches:

1) With the whole system, the igb=4 minimizations are too slow because
the large size of the protein. If I change to cut=16 from without
it converges much faster so I can use much fewer maxcyc.
But I don't know if it's a good choice because a) the energies are
b) MM may be off quite a bit. Even if deltaG may not got affacted much
due to
cancelling-out, still may be larger than it should be.

2) If I cut a reduced system out centered around the ligand, the new
contains a lot of short pieces of individual peptide strands. The
re-creaction of
the receptor and complex prmtops introduces some new H in the
N-terminals and
new OXT in the C-terminals. These new atoms of course are not in the
snapshots taken out of MD trajectories. It looks that you cannot get a
with leap without these new atoms because leap would complain the large
distances between the terminal atoms. I wonder how to get a workable
complex model.

Ant solution is welcome.

Chenglong Li

The AMBER Mail Reflector
To post, send mail to
To unsubscribe, send "unsubscribe amber" to
Received on Wed Nov 19 2003 - 18:53:00 PST
Custom Search