Re: [AMBER] Problem about Sander with large system

From: Adrian Roitberg <roitberg.qtp.ufl.edu>
Date: Mon, 24 May 2010 23:03:24 -0400

What are seeing is a 'bubble'.

Try imagining the system with its neighboring boxes and you will see a
'hole' instead of a funny shape.

This simply means that you had too much volume for the amount of water
you had originally, so water decided to have the correct density itself,
regardless of the box.

There are numerous threads in the list over the years, but what you need
to is to be VERY carefully minimizing, heating and changing from to NPT
very slowly, restraining the protein, the water, and releasing them.

Adrian


On 5/24/10 10:58 PM, Tom Williams wrote:
> Dear All,
>
> I ran a MD of a case soaking a protein molecule with Sander in a rectangle
> PBC box,
> but finally the result on file *.rst, are totally with not a rectangle
> box. The are 4 missing parts
> at the 4 corners of the projection of all the atoms in the the X,Y
> directions.
> In the direction of X, Z, the projection is still with rectangle shape. Pls
> see the attached
> projected images. (Water atoms are red; protein atoms are blue).
>
> I don't think it is because of the PBC images because all the atoms look
> still to be within the area of this rectangle box and my protein is at the
> center of the box.
> Can anyone help explain this strange solvation box shape?
>
> The following is the details:
> ----------------------------------------
> I ran a simple case of soaking a protein molecule in a rectangle box:
> solvatebox LYZ TIP3PBOX 27
>
> My system is with big dimensions. In the file
> of *.inpcrd, the bottom line states:
> 107.0781340 101.4056200 99.1742170 90.0000000 90.0000000 90.0000000
>
> I ran MD with Sander
>
> 0-300K constant temp MD
> &cntrl
> imin=0,
> ntb=1,cut=18.0,
> ntc=2, ntf=2,
> tempi=0.0, temp0=300.0,
> ntt=1,iwrap=1,
> nstlim=10000, dt=0.002,
> ntpr=500, ntwx=1000
> /
>
> Thanks for your help!
>
> Tom
>
>
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber

-- 
                            Dr. Adrian E. Roitberg
                              Associate Professor
                             Quantum Theory Project
                            Department of Chemistry
                  Senior Editor. Journal of Physical Chemistry
                           American Chemical Society
University of Florida                         PHONE 352 392-6972
P.O. Box 118435                               FAX   352 392-8722
Gainesville, FL 32611-8435                    Email adrian.qtp.ufl.edu
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon May 24 2010 - 20:30:03 PDT
Custom Search