[AMBER] amber tutorial13 :constructing and simulating crystals with amber12

From: Yongxiu Li <lyx.itcs.ecnu.edu.cn>
Date: Mon, 22 Apr 2013 20:16:14 +0800

hello everyone,
          I want to construct the unit cell as the tutorial13 and do the
same work as the D.S.Cerutti's paper "Simulations of a protein crystal
with a high resolution X-ray structure: Evaluation of force fields and
water models." /J. Phys. Chem. B/ *114*:12811-12824.
  As the paper, before MD simulation, we need to construct the
simulation Cell from 12 individual unit cells in a 2*2*3 arrangement and
then add ions and water. I do next using amber12:
1) ${AMBERHOME}/bin/PropPDB -p x1AHO.pdb -o x12_1AHO.pdb -X 45.9 -Y 40.7
-Z 30.1 -ix 2 -iy 2 -iz 3
2) ${AMBERHOME}/bin/AddToBox -c x12_1AHO.pdb -a Acetate.pdb -na 84 -o
xa1AHO.pdb -P 49727 -RP 4.0 -RW 6.0 -G 0.1 -V 1 -X 91.8 -Y 81.4 -Z 90.3
-al 90.0 -bt 90.0 -gm 90.0
3) ${AMBERHOME}/bin/AddToBox -c xa1AHO.pdb -a Ammonium.pdb -na 36 -o
xm1AHO.pdb -P 50315 -RP 4.0 -RW 6.0 -G 0.1 -V 1 -X 91.8 -Y 81.4 -Z 90.3
-al 90.0 -bt 90.0 -gm 90.0
4) ${AMBERHOME}/bin/AddToBox -c xm1AHO.pdb -a spce.pdb -na 9000 -o
solv1AHO.pdb -P 50315 -RP 4.0 -RW 6.0 -G 0.1 -V 1 -X 91.8 -Y 81.4 -Z
90.3 -al 90.0 -bt 90.0 -gm 90.0
5) tleap -s -f solvate.tLEaP
6) ${AMBERHOME}/bin/ChBox -c solv1AHO.crd -o solv1AHO.crd -X 91.8 -Y
81.4 -Z 90.3 -al 90.0 -bt 90.0 -gm 90.0

    Question 1: I don't know whether these steps are right expecially
the (-X -Y -Z) value when proppdb ,addtobox and chbox------unit cell
box edge lengths (-X 45.9 -Y 40.7 -Z 30.1 )or supercell box edge
lengths (-X 91.8 -Y 81.4 -Z 90.3)?
    Question2: I do "ambpdb -p solv1AHO.top <solv1AHO.crd
>solv1AHO1.pdb" and then use vmd to open solv1AHO1.pdb .But I see many
warning news such as "Warning) Unusual bond between residues: 3192
(none) and 1342 (waters) " and many water displaying gather together.Why?
     After I finished the six step, I will do minimization using
solv1AHO.top and solv1AHO.crd .We do three minimization as the paper. In
the fist minimization,the final result energy is ****** .In the second
minimization, the min2.out result is also that and we got ************in
min2.rst. Then we can't do the third minimization .
Question3: why these energy are ***** and the bond energy is 0.000? Are
the top-file and crd-file I got wrong?
=====================
The tLEap file:
source leaprc.ff99SB
loadAmberPrep Acetate.prp
loadAmberPrep Ammonium.prp
ammparms = loadAmberParams Ammonium.frcmod
x = loadPdb "solv1AHO.pdb"
bond x.11.SG x.62.SG
bond x.15.SG x.35.SG
bond x.21.SG x.45.SG
bond x.25.SG x.47.SG
bond x.204.SG x.255.SG
bond x.208.SG x.228.SG
....................(there are many lines I can't display)
bond x.8710.SG x.8732.SG
bond x.8889.SG x.8940.SG
bond x.8893.SG x.8913.SG
bond x.8899.SG x.8923.SG
bond x.8903.SG x.8925.SG
bond x.9082.SG x.9133.SG
bond x.9086.SG x.9106.SG
bond x.9092.SG x.9116.SG
bond x.9096.SG x.9118.SG
setBox x vdw 1.0
saveAmberParm x solv1AHO.top solv1AHO.crd
quit
================================
  The min1.in file :
  initial minimisation restraint solvent particles ,added hydrogens
  &cntrl
   imin = 1,
   maxcyc = 40000,
   ncyc = 2000,
   ntb = 1,
   cut = 9.0
   ntr = 1,restraintmask=':1-3072!.H=',
   restraint_wt=1000.0,
  /
  and the min1.out file:
    NSTEP ENERGY RMS GMAX NAME NUMBER
       1 -3.4207E+09 1.5420E+06 8.3886E+06 O 87904

  BOND = 188466.4833 ANGLE = 48742.9387 DIHED = 27311.6719
  VDWAALS = ************* EEL = 2893801.6968 HBOND = 0.0000
  1-4 VDW = 19497.2575 1-4 EEL = 141875.8611 RESTRAINT =
0.0000


    NSTEP ENERGY RMS GMAX NAME NUMBER
      50 -2.8425E+08 1.5407E+06 8.3886E+06 O 71233

  BOND = 188465.4924 ANGLE = 48742.9560 DIHED = 27311.6719
  VDWAALS = ************* EEL = 2893818.9021 HBOND = 0.0000
  1-4 VDW = 19497.2543 1-4 EEL = 141875.8649 RESTRAINT =
0.5613
  EAMBER = *************
.............................................(there are many lines I
can't display)
    NSTEP ENERGY RMS GMAX NAME NUMBER
    2200 -8.4000E+09 2.8326E+05 8.3886E+06 O 32863

  BOND = 0.0000 ANGLE = 138665.0194 DIHED = 28167.4538
  VDWAALS = ************* EEL = 89742545.8658 HBOND = 0.0000
  1-4 VDW = 85150.6573 1-4 EEL = 141892.0475 RESTRAINT =
0.0000
  EAMBER = *************

      .... RESTARTED DUE TO LINMIN FAILURE ...


                     FINAL RESULTS



    NSTEP ENERGY RMS GMAX NAME NUMBER
    2221 -8.3999E+09 2.8326E+05 8.3886E+06 O 32863

  BOND = 0.0000 ANGLE = 138665.0203 DIHED = 28167.4538
  VDWAALS = ************* EEL = 89742545.4937 HBOND = 0.0000
  1-4 VDW = 85150.5481 1-4 EEL = 141892.0475 RESTRAINT =
0.0000
  EAMBER = *************

      ***** REPEATED LINMIN FAILURE *****

===========================
The min2.in file:
  initial minimisation restraint solvent particle
  &cntrl
   imin = 1,
   maxcyc = 40000,
   ncyc = 2000,
   ntb = 1,
   cut = 9.0
   ntr = 1,restraintmask=':3073-18384',
   restraint_wt=1000.0,
  /
AND the min2.out file:
    NSTEP ENERGY RMS GMAX NAME NUMBER
    2050 -1.0153E+10 8.0112E+05 8.3886E+06 HB3 39147

  BOND = 0.0000 ANGLE = 4108575.1870 DIHED = 116768.8191
  VDWAALS = ************* EEL = -154872.1287 HBOND = 0.0000
  1-4 VDW = ************* 1-4 EEL = 68200.0632 RESTRAINT =
0.0000
  EAMBER = *************

      .... RESTARTED DUE TO LINMIN FAILURE ...

      .... RESTARTED DUE TO LINMIN FAILURE ...

      .... RESTARTED DUE TO LINMIN FAILURE ...


                     FINAL RESULTS



    NSTEP ENERGY RMS GMAX NAME NUMBER
    2085 -1.0154E+10 8.0112E+05 8.3886E+06 HB3 39147

  BOND = 0.0000 ANGLE = 4108575.1969 DIHED = 116768.8190
  VDWAALS = ************* EEL = -154872.1019 HBOND = 0.0000
  1-4 VDW = ************* 1-4 EEL = 68200.0630 RESTRAINT =
0.0000
  EAMBER = *************

      ***** REPEATED LINMIN FAILURE *****

==========================
The vmd warning when loading solv1AHO1.pdb :
Warning) Unusual bond between residues: 3183 (none) and 2330 (waters)
Warning) Unusual bond between residues: 3183 (none) and 3109 (waters)
Warning) Unusual bond between residues: 3183 (none) and 5109 (waters)
Warning) Unusual bond between residues: 3183 (none) and 9587 (waters)
Warning) Unusual bond between residues: 3183 (none) and 2048 (waters)
Warning) Unusual bond between residues: 3183 (none) and 2078 (waters)
Warning) Unusual bond between residues: 3183 (none) and 2721 (waters)
Warning) Unusual bond between residues: 3183 (none) and 9587 (waters)
Warning) Unusual bond between residues: 3183 (none) and 213 (waters)
Warning) Unusual bond between residues: 3183 (none) and 982 (waters)
Warning) Unusual bond between residues: 3183 (none) and 1310 (waters)
Warning) Unusual bond between residues: 3183 (none) and 2048 (waters)
Warning) Unusual bond between residues: 3183 (none) and 2078 (waters)
Warning) Unusual bond between residues: 3183 (none) and 2721 (waters)
Warning) Unusual bond between residues: 3183 (none) and 3085 (waters)
Warning) Unusual bond between residues: 3183 (none) and 3364 (waters)
Warning) Unusual bond between residues: 3183 (none) and 5154 (waters)
Warning) Unusual bond between residues: 3183 (none) and 5562 (waters)
Warning) Unusual bond between residues: 3184 (none) and 9548 (waters)
Warning) Unusual bond between residues: 3184 (none) and 9548 (waters)
================

          Thank you very much!
          Yongxiu Li

_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon Apr 22 2013 - 05:30:03 PDT
Custom Search