spoel at xray.bmc.uu.se
Fri Jul 11 21:17:01 CEST 2003
On Fri, 2003-07-11 at 20:39, Nagy, Peter I. wrote:
> Dear David!
> I assume that when you write "make it 10 times as big" you mean the
> volume. I increased the edges by a factor of 2.2-2.5, and resolvated
> sytem. More than 10,000 water were added that I think to be enormously
Forget about the water! To test whether you have a PBC problem, you
don't need any water.
Now, if your molecule is elongated (as it probably is) you can use
editconf -princ to align it along the principal components, i.e. the
long axis. You could the use a rectangular box and save on the amount of
water, however there is a possibility that in a simulation the molecule
will rotate by 90 degrees and interact with itself, which you don't want
> many. With a smaller box I still got 5758 waters. I think it is still
> too many. In Sybyl cca 1500 were added. I try to find a box including
> about 3000
> water molecules, I assume they would suffice? Do you agree?
> Peter Nagy
> -----Original Message-----
> From: David [mailto:spoel at xray.bmc.uu.se]
> Sent: Fri 7/11/2003 2:15 PM
> To: gmx-users at gromacs.org
> Subject: RE: [gmx-users] helices
> On Fri, 2003-07-11 at 20:06, Nagy, Peter I. wrote:
> > Dear David!
> > Of course, it may be a periodic boundary effect. But then what
> > would I do?
> > Should I remarkably increase the box-size accomodating much more
> > molecules and repeating the procedure?
> Just the box size will do for a test, just edit the last line in the
> file, make it 10 times as big. mdrun will always keep the proteins in
> the box, but the algorithm makes that proteins can easily jump. If it
> a pbc effect you box is definitely too small, unless this is a
> Is it maybe the ice-binding protein?
> > Peter Nagy
> > -----Original Message-----
> > From: David [mailto:spoel at xray.bmc.uu.se]
> > Sent: Fri 7/11/2003 2:03 PM
> > To: gmx-users at gromacs.org
> > Cc:
> > Subject: Re: [gmx-users] helices
> > On Fri, 2003-07-11 at 19:44, Nagy, Peter I. wrote:
> > > Dear Gromacs Users!
> > >
> > > I performed a test run for two close, nearly parallel alpha
> > helices
> > > with
> > > 26+26+NAC+ACE residues. The net charge of the system was just zero
> > > because of having only one ASP and one LYS residues. The helices
> > were
> > > slightly shifted with respect to each other, and the header of the
> > > second helix was at a height of 1/3 total length of the first one
> > (so
> > > 2/3 of both helices were close and parallel). The system was
> > solvated
> > > by 620 SPC waters.
> > > In performing the energy minimization, I used the steepes
> > > method, setting rlist, rcoulomb, and rvdw to 0.9 (nm).
> > > Minimization stopped at 780 out of 1000 steps but did not reach
> > the
> > > emtol=
> > > 100 limit. Reading the recent mails on this list, I was not
> > I
> > > started
> > > a 200 ps MD run, and checked the .xtc trajectory in ngmx.
> > > Ngmx starts displaying structures at t=0 ps, thus with the
> > > structure obtained as the output of the energy minimization. This
> > is,
> > > at least, that I
> > > presumed. Then I was very much surprised that the structure at
> > > point was
> > > close to a collinear pair of two alpha helices. I understood it
> > > this structure was created by the energy minimization.
> > > I performed an energy minimization of the original dimer
> > structure
> > > using
> > > the Sybyl package, as well. The system in this case consisted of
> > about
> > > 1500
> > > water molecules + the dimer. I used periodic boundary, TIP3P water
> > > molecules
> > > and allowd 1000 steps with the conjugate gradient method. Although
> > the
> > > two
> > > energy minizmizations differ in some technical elements, but it
> > > not explain for me the result: the dimer structure was pactically
> > > unchanged in
> > > Sybyl after 1000 steps.
> > > Then why did the GROMACS energy minimizer change (if it really
> > > changed, as I suspect) the relative positions of the two helices?
> > Simplest suggestion first: isn't it a periodic boundary effect?
> > minimization will not move atoms more than a fraction of an
> > >
> > >
> > > Peter Nagy
> > > Dept. Medicinal and Biological Chemistry
> > > The University of Toledo
> > > Toledo, OH 43606 USA
> > pnagy at utnet.utoledo.edu
> > --
> > Groeten, David.
> > Dr. David van der Spoel, Dept. of Cell and Molecular Biology
> > Husargatan 3, Box 596, 75124 Uppsala, Sweden
> > phone: 46 18 471 4205 fax: 46 18 511 755
> > spoel at xray.bmc.uu.se spoel at gromacs.org
> > http://xray.bmc.uu.se/~spoel
> > _______________________________________________
> > gmx-users mailing list
> > gmx-users at gromacs.org
> > http://www.gromacs.org/mailman/listinfo/gmx-users
> > Please don't post (un)subscribe requests to the list. Use the
> > www interface or send it to gmx-users-request at gromacs.org.
> Groeten, David.
> Dr. David van der Spoel, Dept. of Cell and Molecular Biology
> Husargatan 3, Box 596, 75124 Uppsala, Sweden
> phone: 46 18 471 4205 fax: 46 18 511 755
> spoel at xray.bmc.uu.se spoel at gromacs.org
> gmx-users mailing list
> gmx-users at gromacs.org
> Please don't post (un)subscribe requests to the list. Use the
> www interface or send it to gmx-users-request at gromacs.org.
Dr. David van der Spoel, Dept. of Cell and Molecular Biology
Husargatan 3, Box 596, 75124 Uppsala, Sweden
phone: 46 18 471 4205 fax: 46 18 511 755
spoel at xray.bmc.uu.se spoel at gromacs.org http://xray.bmc.uu.se/~spoel
More information about the gromacs.org_gmx-users