> Good energies or not, two of my subunits do still not fit into the
> rest of the complex properly. Total rmsd of 70Å when I do not reimage
> the wrapped restart file, and 1.6Å when I do reimage the wrapped
> restart file (with ptraj).
>
> You mentioned that my subunit complex could be spanning the primary
> unit cell, and that two subunits are therefore individually
> wrapped(?). This makes sense, but should it not be possible to reimage
> this with ptraj and not get an rmsd of 1,6Å? I can clearly see that
> two subunits are shifted.
Lars sent me his coordinates and I can verify that the coordinates before
and after imaging are essentially equivalent with the exception that one
of the subunits is in an adjacent periodic image cell after PMEMD imaging
and not ptraj imaging. The difference likely results due to the use of
the "center" command in ptraj to center the protein prior to imaging. I
can also verify that this is an incredibly large system!
Lars, if you do not center before imaging in ptraj likely you will get the
same thing as you saw with PMEMD. What is happening is that the whole
protein complex is moving in the box (as expected). One of the subunits
has its center of mass in an adjacent unit cell so when imaged moves to
the other side of the unit cell. There is no way to easily prevent this
except by (a) restraining the whole complex to stay at the center of the
box or by (b) altering the prmtop to treat all the protein as a single
molecule. You can however, most likely, through clever imaging in ptraj
"fix" this after the fact likely using the center/image commands.
So, I think you can continue with your PMEMD wrapping (i.e. restart from
the b4 restrt and continue running turning on IWRAP=1) and then later use
ptraj to put the complex back together again. [e-mail me offline if you
want to have the prmtop modified so that the protein will never
"dissociate" to periodic images...]
--tom
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
to majordomo.scripps.edu
Received on Sun Aug 31 2008 - 06:07:40 PDT
