Re: [AMBER] Reorienting the protein during simulation

From: Daniel Roe <>
Date: Tue, 16 Jan 2018 09:54:33 -0500


The important thing to remember is that with periodic boundary
conditions your protein is *never* really outside the box. It may
appear to be outside the boundaries of your unit cell, but it is
always interacting with its images. If your protein is rotating so
that it is no longer aligned with your unit cell axes it is because it
is favorable to do so. Your box may be too small. If you absolutely
need the protein to remain aligned with e.g. the Z axis you'll
probably have to apply some external restraint and make sure you
correct for it in any later analysis if necessary.

On Tue, Jan 16, 2018 at 5:52 AM, Ahmed Mashaly <> wrote:
> This is what I use for recentering>
> trajin my_rst_file.r>
> autoimage
> center: 1-1500 mass origin>
> image origin center

The 'center' and 'image' commands in this case will override whatever
was done with 'autoimage', so it's not useful to combine the two.
Either use 'autoimage' (recommended) *or* use 'image'/'center'.

> I also tried to add>
> prinicipal :1-1500 dorotation mass

Note that the rotation done by 'principal' means that further imaging
after this will not be possible.


Daniel R. Roe
Laboratory of Computational Biology
National Institutes of Health, NHLBI
5635 Fishers Ln, Rm T900
Rockville MD, 20852
AMBER mailing list
Received on Tue Jan 16 2018 - 07:00:04 PST
Custom Search