Re: [AMBER] Error in setting irest=1

From: Soumendranath Bhakat <bhakatsoumendranath.gmail.com>
Date: Wed, 21 May 2014 09:19:58 +0200

Huh!!

Your script is wrong during minimisation. Use these scripts for
minimisation. Once you into MD step after some initial velocity if you put
irest=0, and ig=-1 and tempi and temp0 equal to 300K then it will assign
initial velocities based on boltzman weighed.

Attach herewith all the inputs that might help you

Partial_Mini.in
Initial minimization of MMP3 (MMMM): solvent molecules and added ions
 &cntrl
  imin = 1,
  maxcyc = 2500,
  ncyc = 750,
  ntb = 1,
  ntr = 1,
  cut = 12.0,
 /
Hold the Protein fixed
10.0
RES 1 557
END
END

Full_Mini.in
full minimization of MMP3 (MMMM): protein, ligand, solvent molecules and
added ions
 &cntrl
  imin = 1,
  maxcyc = 200,
  ncyc = 50,
  ntb = 1,
  ntr = 0,
  cut = 12.0,
  Drms = 0.0001,
/
END

Heating.in
Heating Step of MMP3 (MMMM): stage-5
 &cntrl
  imin= 0,
  irest=1,
  NTX=1,
  ntb= 1,
  NTPR=500,
  NTWX=500,
  NTWR=500,
  ntr=1,
  Tempi=0.0,
  Temp0=300.0,
  NTT=3,
  gamma_ln=1.0,
  NTC=2,
  NTF=2,
  cut= 12.0,
  nstlim=2500,
  dt=0.002,
/
Keep Protein and inhibitor fixed with weak restraints
10.0
RES 1 557
END
END

equil.in
Equilibration Step of MMP3 (MMMM): stage-1
 &cntrl
  imin= 0,
  irest=1,
  NTX=7,
  ntb=2,
  ntp=1,
  PRES0=1.0,
  TAUP=2.0,
  NTPR=500,
  NTWX=500,
  ntr=0,
  Tempi=300.0,
  Temp0=300.0,
  NTT=3,
  gamma_ln=1.0,
  NTC=2,
  NTF=2,
  cut=12.0,
  nstlim=250000,
  dt=0.002
 /

md.in
Equilibration Step of MMP3 (MMMM): stage-1
 &cntrl
  imin= 0,
  irest=0,
  ig=-1
  NTX=7,
  ntb=2,
  ntp=1,
  PRES0=1.0,
  TAUP=2.0,
  NTPR=500,
  NTWX=500,
  ntr=0,
  Tempi=300.0,
  Temp0=300.0,
  NTT=3,
  gamma_ln=1.0,
  NTC=2,
  NTF=2,
  cut=12.0,
  nstlim=2500000,
  dt=0.002
 /



On Tue, May 20, 2014 at 6:39 PM, Daniel Roe <daniel.r.roe.gmail.com> wrote:

> Hi,
>
> Minimization does not use velocity information, only MD (which is why
> you don't have velocity info in your restart to begin with), so you
> don't need to set 'irest' if you are continuing a minimization.
>
> -Dan
>
> On Tue, May 20, 2014 at 10:14 AM, Valentina Romano
> <valentina.romano.unibas.ch> wrote:
> > Dear users
> >
> > I ran a min. step for 2500 steps.
> > Since the energy was not stable and the RMS still high i decided to run
> a second min. (2500 additional steps) using as starting point the files i
> got from the previous step.
> >
> > In the first step:
> > -r PknGAde-min.rst and -ref PknGAde-solv-min.rst
> >
> > The 2nd step input file was the following:
> >
> > Initial XMIN minimisation of PknG-Adenine complex: complex, water and
> ions
> > &cntrl
> > imin=1,
> > maxcyc=2500,
> > ntmin=3,
> > irest=1,
> > ntx=5,
> > ntb=1,
> > igb=0,
> > ntr=0,
> > cut=10,
> > ntpr=10, ntwx=10, ntwr=1000
> >
> > And i ran sander as follow:
> > sander -O -i PknGAde-min.in -o PknGAde-min.out -p
> ../PknGAde_params/PknGHAdeH_ion_wt.prmtop -c PknGAde-min.rst -r
> PknGAde-min.rst -ref PknGAde-solv-min.rst -x PknGAde-min.mdcrd &
> >
> > Thus I used as coord. the file I got from the previuos min.
> >
> > I got the following error:
> > ERROR: I could not find enough velocities in
> > PknGAde-min.rst
> >
> > I do not understand why i got that error.
> > Anyone can help me in understanding what happened?
> >
> > Thank you
> > Valentina
> >
> >
> > ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
> > Valentina Romano | PhD Student | Biozentrum, University of Basel & SIB
> Swiss Institute of Bioinformatics
> > Klingelbergstrasse 61 | CH-4056 Basel |
> >
> > Phone: +41 61 267 15 80
> >
> >
> >
> > _______________________________________________
> > AMBER mailing list
> > AMBER.ambermd.org
> > http://lists.ambermd.org/mailman/listinfo/amber
>
>
>
> --
> -------------------------
> Daniel R. Roe, PhD
> Department of Medicinal Chemistry
> University of Utah
> 30 South 2000 East, Room 201
> Salt Lake City, UT 84112-5820
> http://home.chpc.utah.edu/~cheatham/
> (801) 587-9652
> (801) 585-6208 (Fax)
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>



-- 
Thanks & Regards;
Soumendranath Bhakat
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed May 21 2014 - 00:30:02 PDT
Custom Search