AMBER Archive (2007)Subject: AMBER: Water cap MD error
From: Rima Chaudhuri (rima.chaudhuri_at_gmail.com)
Date: Mon Oct 15 2007 - 17:17:07 CDT
Hello!
Does amber9 support solvating only a small portion of the protein (say
active site with ligand) with a TIP3PBOX water cap (25A) around it, using
igb=10? or does the entire protein have to be within the cap?
I am working on a dimer with a ligand, a system of ~6000 atoms. If I try to
use explicit solvent (using a water box), a mere 8A box hikes the number to
almost 90k atoms, which is too big for our computer system. I tried to run
the water cap option with the following min and equil steps which ran fine,
on trying the production run, the md output gives errors of "vlimit
exceeded". Earlier I tried using the implicit solvent option of igb=2,
gbsa=1, but does not work for my system, one of the subunits of the dimer
seem to move away from the other quite rapidly. So, am I missing any options
that I could try and also what am I doing wrong here?..
min.in
&cntrl
imin=1, maxcyc=5000, ncyc=1000, ntb=0, ntpr=20, scee=1.2, nsnb=25,
ntmin=2, cut=12,
&end
FINAL RESULTS
NSTEP ENERGY RMS GMAX NAME NUMBER
5000 -1.2772E+04 2.7340E-01 4.6052E+00 O 10911
BOND = 564.0388 ANGLE = 1303.0728 DIHED =
5820.9651
VDWAALS = -4605.3428 EEL = -47067.8395 HBOND =
0.0000
1-4 VDW = 2088.6990 1-4 EEL = 29124.6525 RESTRAINT =
0.2275
EAMBER = -12771.7542
equil.in
&cntrl
imin=0, irest=0, ntx=1, ntb=0, cut=12, ntc=2, scee=1.2,
ntf=2, tempi=0.0, temp0=300.0, ntt=1, taup=0.1, dtemp=2.0, ntp=0, nsnb=25,
nstlim=20000, dt=0.002, ntpr=250, ntwx=250, ntwr=1000,
&end
A V E R A G E S O V E R 5000 S T E P S
NSTEP = 5000 TIME(PS) = 10.000 TEMP(K) = 237.81 PRESS =
0.0
Etot = -5968.9244 EKtot = 5828.3659 EPtot = -
11797.2903
BOND = 1192.5247 ANGLE = 3215.4731 DIHED =
6228.2326
1-4 NB = 2200.3827 1-4 EEL = 29157.2263 VDWAALS = -
5406.3059
EELEC = -41625.3623 EGB = -9040.6615 RESTRAINT =
2149.6171
ESURF= 131.5828
EAMBER (non-restraint) = -13946.9074
------------------------------------------------------------------------------
R M S F L U C T U A T I O N S
NSTEP = 5000 TIME(PS) = 10.000 TEMP(K) = 74.91 PRESS =
0.0
Etot = 3736.7361 EKtot = 1835.9481 EPtot =
1902.7262
BOND = 356.7281 ANGLE = 756.1285 DIHED =
177.5090
1-4 NB = 57.5145 1-4 EEL = 39.3693 VDWAALS =
61.2788
EELEC = 339.9597 EGB = 328.2173 RESTRAINT =
526.3146
ESURF= 0.3270
EAMBER (non-restraint) = 1376.4116
Now, the troublesome MD run:
&cntrl
ntx=5, irest=1, imin=0,
ntpr=500, ntwx=500, ntwr=5000,
nstlim=100000, dt=0.001,
ntt=1, tempi=300.0, temp0=300.0, tautp=0.1,
igb=10, ntb=0, cut=0, fcap=10.0, ivcap=0,
ntc=2, ntf=2, tol=0.000001,
&end
RESULTS
vlimit exceeded for step 29; vmax = 27.5043
vlimit exceeded for step 30; vmax = 93.7706
vlimit exceeded for step 31; vmax = 38.9037
vlimit exceeded for step 32; vmax = 80.5579
vlimit exceeded for step 33; vmax = 137.4696
Thanks, any help will be much appreciated!
--
-Rima
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu
|