Re: [AMBER] internal geometry restraint

From: case <>
Date: Tue, 26 Nov 2013 17:56:11 -0500

On Tue, Nov 26, 2013, Jio M wrote:
> I have three linear polymers and I have placed them parallel to each
> other. I am not sure they have been placed 'very close' because my
> script just considered the atoms as point coordinates while placing
> them. So when I run a simulation they are repelled to each other in
> just 50 ps and the internal geometry of each polymer chain is destroyed
> badly.

> What I am thinking is to run first some simulation considering three
> polymer chains as 'hard rods', that is keeping internal geometry
> coordinates of a polymer constrained (all atoms in a polymer chain being
> fixed with respect to each other; keeping internal geometry of polymer
> chain fixed) so that they just act as hard rods and can adjust their
> positions and later I release the restraints.

You can carry out a version of this sort of calculation using "EMAP"
constraints, as explained in section 6.13 of the manual. This is not exactly
what you describe, but it might suffice. You'll probably have to experiment a
bit (or a lot!) with force constants and the resolution of the maps you are
using. Some other dynamics programs do a better job of allowing you to define
rigid groups that can move relative to each other.


AMBER mailing list
Received on Tue Nov 26 2013 - 15:00:03 PST
Custom Search