[AMBER] parameterization with crystal water

From: Hira Jabeen <hira_bioinfo.yahoo.com>
Date: Wed, 25 Mar 2020 17:53:52 +0000 (UTC)

Hi,
I am trying to parameterize protein and ligand complex along with crystal water.
Can you please direct me how to generate? There are 560 water in crystal structure of protein.
I have already run it through this command:
antechamber -fi pdb -fo mol2 -i WAT_H.pdb -o WAT.mol2-at amber -c bcc -pf y  -dr no
Here is the input for sqm:Run semi-empirical minimization
 &qmmm
    qm_theory='AM1', grms_tol=0.0005,
 scfconv=1.d-10, ndiis_attempts=700,   qmcharge=0,
 /
   8     O      -14.6250       -1.5500       21.3030
   1    H1      -13.6680       -1.5500       21.3030
   1    H2      -14.8650       -0.6230       21.3030
   8     O      -26.6640      -27.0540       16.5150
   1    H1      -25.7070      -27.0540       16.5150
   1    H2      -26.9040      -26.1270       16.5150
   8     O      -19.0830      -17.0780       21.6480
   1    H1      -18.1260      -17.0780       21.6480
   1    H2      -19.3230      -16.1510       21.6480
   8     O       -7.9350       -6.3420        1.2950....

And the output file (sqm.out) stopped at this:
 …….
  QMMM:   989      989      H      -17.4370   -0.9490   36.9620
  QMMM:   990      990      H      -18.6340   -0.0220   36.9620
  QMMM:   991      991      O       -9.7060   10.7040   17.9070
  QMMM:   992      992      H       -8.7490   10.7040   17.9070
  QMMM:   993      993      H       -9.9460   11.6310   17.9070
  QMMM:   994      994      O      -26.1090    6.2790   19.7550
  QMMM:   995      995      H      -25.1520    6.2790   19.7550
  QMMM:   996      996      H      -26.3490    7.2060   19.7550
  QMMM:   997      997      O      -23.4160    9.1110   18.6120
  QMMM:   998      998      H      -22.4590    9.1110   18.6120
  QMMM:   999      999      H      -23.6560   10.0380   18.6120--------------------------------------------------------------------------------
  RESULTS
--------------------------------------------------------------------------------      iter         sqm energy              rms gradient
      ----    -------------------    -----------------------
xmin    10   -20278.8965 kcal/mol        2.8979 kcal/(mol*A)
xmin    20   -20820.3242 kcal/mol        2.7131 kcal/(mol*A)
xmin    30   -21129.5884 kcal/mol        1.3588 kcal/(mol*A)
xmin    40   -21327.1301 kcal/mol        1.8634 kcal/(mol*A)
Its like this from last 24 hours, stopped at this last line (no other log or output file for error). I don't know, is this the right way to parameterize along with all water or should I just take 5 angstrom crystal water around substrate? 
Looking forward to see your response.
Thanks,Hira Jabeen


_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed Mar 25 2020 - 11:00:02 PDT
Custom Search