AMBER Archive (2005)

Subject: RE: AMBER: minimize a ptraj average pdb structure

From: Bill Ross (ross_at_cgl.ucsf.edu)
Date: Wed Feb 09 2005 - 13:39:27 CST


> >Hi, I have obtained an average structure for a protein
> >in water. Now the question is how to minimize the
> >average structure. I suppose that I should limit the
> >number of steepest descent (NCYC) while keeping the
> >average water molecule. Can people give some
> >suggestions based on their experience?

If your trajectory ran any appreciable amount of time, the waters
will need to be discarded.

> Since you don't want the structure to change that much, in my opinion you
> could simply control it by switching rmsd of gradients....another opinion is
> to proceed with a gradual progressive minimization, like minimizing first
> hydrogens, then side chains and then all the solute with backbone fixed by a
> force constant. In this last step you could choose an RMSD value say
> 0.1......and not 0.0001.....you will then have just a clean structure,
> without any bumps, that is perhaps you need..

The same 500 steps of steepest descent that one would normally
use to prepare for initial dynamics will probably suffice.

Unless there were major motions that result in average positions
that involve, say, bonds threading through rings or the like. Which
is the sort of thing that may be found with the waters.

Bill
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu