Re: AMBER: RESP charges

From: Youyi Peng <pengyo.UMDNJ.EDU>
Date: Thu, 18 Dec 2003 10:20:43 -0500

Hi Junmei,
Following your instructions, I went through all the procedure to calculate
RESP charges. However, the resp charges from the first and second stage
fitting didn't change at all and have very high value. Something must be
wrong.
My system includes six single molecules (a Mg coordination system). But all
g98, antechamber and resp treat it as one single molecule. All tutorial and
information seem suitable for only single molecule. I don't know what to deal
with the situation. I list my mol2, gcrt,ac and respin files here, and also
the resp chareges I got for this system. Would you please take a while to look
at how to modify these files.
Thanks a lot.

Youyi Peng

mol2 file
.<TRIPOS>MOLECULE
MGSYS
   45 39 6 0 0
SMALL
USER_CHARGES
.<TRIPOS>ATOM
      1 CE 38.0430 60.6410 7.1880 C.3 1 LYS+33 -0.1030
****
      2 HE3 38.1579 61.6791 6.8761 H 1 LYS+33 0.0660
****
      3 HE1 37.2002 60.1703 6.6818 H 1 LYS+33 0.0660
****
      4 HE2 37.9358 60.5554 8.2693 H 1 LYS+33 0.0660
****
      5 NZ 39.2610 59.9180 6.7990 N.3 1 LYS+33 -0.0730
****
      6 HZ1 39.1779 58.9298 7.0773 H 1 LYS+33 0.3260
****
      7 HZ2 39.3849 59.9775 5.7782 H 1 LYS+33 0.3260
****
      8 HZ3 40.0755 60.3390 7.2684 H 1 LYS+33 0.3260
****
      9 PG 37.5270 54.4619 3.8299 P 2 ATP381 1.4930
****
     10 O1G 36.0790 54.0039 3.8239 O.2 2 ATP381 -0.9130
****
     11 O2G 38.5130 53.6969 2.9709 O.2 2 ATP381 -0.9130
****
     12 O3G 37.9270 54.1519 5.3139 O.2 2 ATP381 -0.9130
****
     13 PB 36.9230 57.0059 4.7339 P 2 ATP381 1.3340
****
     14 O1B 35.5110 57.0769 4.2769 O.2 2 ATP381 -0.9130
****
     15 O2B 37.0600 56.6369 6.1029 O.2 2 ATP381 -0.9130
****
     16 O3B 37.5990 56.0519 3.6339 O.3 2 ATP381 -0.5080
****
     17 PA 38.6220 58.8639 3.3589 P 2 ATP381 1.3340
****
     18 O1A 39.6740 59.7669 3.8329 O.2 2 ATP381 -0.9130
****
     19 O2A 39.3220 57.7569 2.5729 O.2 2 ATP381 -0.9130
****
     20 O3A 37.6620 58.3999 4.4999 O.3 2 ATP381 -0.5080
****
     21 O5* 37.6940 59.4409 2.2009 O.3 2 ATP381 -0.4910
****
     22 C5* 36.5940 58.9739 1.4229 C.3 2 ATP381 0.0390
****
     23 H5*3 36.7413 59.2754 0.3858 H 2 ATP381 0.0660
****
     24 H5*1 35.6748 59.4119 1.8117 H 2 ATP381 0.0660
****
     25 H5*2 36.5479 57.8871 1.4924 H 2 ATP381 0.0660
****
     26 MG 39.5010 55.7059 2.7639 Mg 3 MG383 2.0000
****
     27 CB 43.3161 55.2559 0.8887 C.3 4 ASN132 -0.2930
****
     28 HB3 42.9768 55.9158 0.0903 H 4 ASN132 0.0660
****
     29 HB1 43.9554 55.7978 1.5857 H 4 ASN132 0.0660
****
     30 HB2 43.8686 54.3951 0.5121 H 4 ASN132 0.0660
****
     31 CG 42.0841 54.7459 1.6547 C.2 4 ASN132 0.7260
****
     32 OD1 41.0571 55.4209 1.5367 O.2 4 ASN132 -0.4920
****
     33 ND2 42.1431 53.6889 2.5147 N.3 4 ASN132 -0.6870
****
     34 HD21 43.0105 53.1396 2.5969 H 4 ASN132 0.2740
****
     35 HD22 41.3202 53.4400 3.0820 H 4 ASN132 0.2740
****
     36 O 38.2370 55.5959 1.0426 O.3 5 HOH111 -0.8200
****
     37 HO1 39.2670 55.5959 1.0426 H 5 HOH111 0.4100
****
     38 HO2 37.8932 56.5668 1.0426 H 5 HOH111 0.4100
****
     39 CB 42.1651 57.1282 3.9506 C.3 6 ASP-145 -0.2930
****
     40 HB3 42.7444 57.8504 4.5259 H 6 ASP-145 0.0660
****
     41 HB1 42.8388 56.4040 3.4926 H 6 ASP-145 0.0660
****
     42 HB2 41.6064 57.6476 3.1720 H 6 ASP-145 0.0660
****
     43 CG 41.2021 56.4102 4.8676 C.2 6 ASP-145 0.6790
****
     44 OD1 41.2281 56.6912 6.0816 O.2 6 ASP-145 -0.7920
****
     45 OD2 40.4021 55.5492 4.4356 O.2 6 ASP-145 -0.7920
****
.<TRIPOS>BOND
    1 1 5 1
    2 1 3 1
    3 1 4 1
    4 1 2 1
    5 5 6 1
    6 5 7 1
    7 5 8 1
    8 9 10 1
    9 9 11 1
   10 9 12 1
   11 9 16 1
   12 13 14 1
   13 13 15 1
   14 13 16 1
   15 13 20 1
   16 17 18 1
   17 17 19 1
   18 17 20 1
   19 17 21 1
   20 21 22 1
   21 22 24 1
   22 22 25 1
   23 22 23 1
   24 27 31 1
   25 27 29 1
   26 27 30 1
   27 27 28 1
   28 31 32 2
   29 31 33 ar
   30 33 34 1
   31 33 35 1
   32 36 37 1
   33 36 38 1
   34 39 43 1
   35 39 41 1
   36 39 42 1
   37 39 40 1
   38 43 44 ar
   39 43 45 ar
.<TRIPOS>SUBSTRUCTURE
      1 LYS+33 1 RESIDUE 0 A LYS+
      2 ATP381 9 **** 0 A ATP
      3 MG383 26 **** 0 A MG
      4 ASN132 27 RESIDUE 0 A ASN
      5 HOH111 36 **** 0 Z HOH
      6 ASP-145 39 RESIDUE 0 A ASP-

gcrt file
#HF/6-31G* SCF=tight Test Pop=MK iop(6/33=2) iop(6/42=6) opt

remark line goes here

-2 1
    C 38.0430 60.6410 7.1880
    H 38.1579 61.6791 6.8761
    H 37.2002 60.1703 6.6818
    H 37.9358 60.5554 8.2693
    N 39.2610 59.9180 6.7990
    H 39.1779 58.9298 7.0773
    H 39.3849 59.9775 5.7782
    H 40.0755 60.3390 7.2684
    P 37.5270 54.4619 3.8299
    O 36.0790 54.0039 3.8239
    O 38.5130 53.6969 2.9709
    O 37.9270 54.1519 5.3139
    P 36.9230 57.0059 4.7339
    O 35.5110 57.0769 4.2769
    O 37.0600 56.6369 6.1029
    O 37.5990 56.0519 3.6339
    P 38.6220 58.8639 3.3589
    O 39.6740 59.7669 3.8329
    O 39.3220 57.7569 2.5729
    O 37.6620 58.3999 4.4999
    O 37.6940 59.4409 2.2009
    C 36.5940 58.9739 1.4229
    H 36.7413 59.2754 0.3858
    H 35.6748 59.4119 1.8117
    H 36.5479 57.8871 1.4924
   Mg 39.5010 55.7059 2.7639
    C 43.3161 55.2559 0.8887
    H 42.9768 55.9158 0.0903
    H 43.9554 55.7978 1.5857
    H 43.8686 54.3951 0.5121
    C 42.0841 54.7459 1.6547
    O 41.0571 55.4209 1.5367
    N 42.1431 53.6889 2.5147
    H 43.0105 53.1396 2.5969
    H 41.3202 53.4400 3.0820
    O 38.2370 55.5959 1.0426
    H 39.2670 55.5959 1.0426
    H 37.8932 56.5668 1.0426
    C 42.1651 57.1282 3.9506
    H 42.7444 57.8504 4.5259
    H 42.8388 56.4040 3.4926
    H 41.6064 57.6476 3.1720
    C 41.2021 56.4102 4.8676
    O 41.2281 56.6912 6.0816
    O 40.4021 55.5492 4.4356

ac file
CHARGE 0.00
Formula: H19C6N2O14Mg1P3
ATOM 1 C1 MOL 1 -4.703 -2.167 -1.268 0.000000
C C
ATOM 2 H1 MOL 1 -5.178 -2.874 -0.588 0.000000
H H
ATOM 3 H2 MOL 1 -4.670 -1.168 -0.834 0.000000
H H
ATOM 4 H3 MOL 1 -5.200 -2.147 -2.238 0.000000
H H
ATOM 5 N1 MOL 1 -3.319 -2.611 -1.481 0.000000
N N
ATOM 6 H4 MOL 1 -2.843 -1.956 -2.118 0.000000
H H
ATOM 7 H5 MOL 1 -2.823 -2.634 -0.579 0.000000
H H
ATOM 8 H6 MOL 1 -3.320 -3.553 -1.896 0.000000
H H
ATOM 9 P1 MOL 1 0.581 2.471 -0.734 0.000000
P P
ATOM 10 O1 MOL 1 0.017 3.875 -0.609 0.000000
O O
ATOM 11 O2 MOL 1 2.066 2.283 -0.499 0.000000
O O
ATOM 12 O3 MOL 1 0.355 2.198 -2.261 0.000000
O O
ATOM 13 P2 MOL 1 -1.868 1.242 -0.356 0.000000
P P
ATOM 14 O4 MOL 1 -2.557 2.335 0.379 0.000000
O O
ATOM 15 O5 MOL 1 -2.161 1.221 -1.750 0.000000
O O
ATOM 16 O6 MOL 1 -0.336 1.438 0.081 0.000000
O O
ATOM 17 P3 MOL 1 -1.453 -1.087 1.267 0.000000
P P
ATOM 18 O7 MOL 1 -1.631 -2.516 0.996 0.000000
O O
ATOM 19 O8 MOL 1 0.053 -0.847 1.357 0.000000
O O
ATOM 20 O9 MOL 1 -2.239 -0.178 0.269 0.000000
O O
ATOM 21 O10 MOL 1 -1.868 -0.610 2.728 0.000000
O O
ATOM 22 C2 MOL 1 -1.869 0.615 3.459 0.000000
C C
ATOM 23 H7 MOL 1 -1.510 0.423 4.470 0.000000
H H
ATOM 24 H8 MOL 1 -2.885 1.007 3.491 0.000000
H H
ATOM 25 H9 MOL 1 -1.211 1.325 2.958 0.000000
H H
ATOM 26 Mg1 MOL 1 1.427 0.282 0.303 0.000000
Mg Mg
ATOM 27 C3 MOL 1 4.862 -2.175 0.966 0.000000
C C
ATOM 28 H10 MOL 1 4.584 -2.242 2.018 0.000000
H H
ATOM 29 H11 MOL 1 4.573 -3.083 0.436 0.000000
H H
ATOM 30 H12 MOL 1 5.932 -2.020 0.827 0.000000
H H
ATOM 31 C4 MOL 1 4.113 -0.987 0.339 0.000000
C C
ATOM 32 O11 MOL 1 3.104 -0.607 0.941 0.000000
O O
ATOM 33 N2 MOL 1 4.457 -0.461 -0.871 0.000000
N N
ATOM 34 H13 MOL 1 5.303 -0.794 -1.355 0.000000
H H
ATOM 35 H14 MOL 1 3.871 0.269 -1.301 0.000000
H H
ATOM 36 O12 MOL 1 1.519 1.510 2.052 0.000000
O O
ATOM 37 H15 MOL 1 2.137 0.716 1.832 0.000000
H H
ATOM 38 H16 MOL 1 0.672 1.166 2.525 0.000000
H H
ATOM 39 C5 MOL 1 1.552 -2.792 -0.728 0.000000
C C
ATOM 40 H17 MOL 1 1.163 -3.754 -1.063 0.000000
H H
ATOM 41 H18 MOL 1 2.641 -2.807 -0.765 0.000000
H H
ATOM 42 H19 MOL 1 1.226 -2.604 0.295 0.000000
H H
ATOM 43 C6 MOL 1 1.033 -1.697 -1.631 0.000000
C C
ATOM 44 O13 MOL 1 0.315 -2.024 -2.596 0.000000
O O
ATOM 45 O14 MOL 1 1.316 -0.493 -1.434 0.000000
O O
BOND 1 1 2 1 C1 H1
BOND 2 1 3 1 C1 H2
BOND 3 1 4 1 C1 H3
BOND 4 1 5 1 C1 N1
BOND 5 5 6 1 N1 H4
BOND 6 5 7 1 N1 H5
BOND 7 5 8 1 N1 H6
BOND 8 9 10 2 P1 O1
BOND 9 9 11 9 P1 O2
BOND 10 9 12 9 P1 O3
BOND 11 9 16 1 P1 O6
BOND 12 13 14 2 P2 O4
BOND 13 13 15 9 P2 O5
BOND 14 13 16 1 P2 O6
BOND 15 13 20 1 P2 O9
BOND 16 17 18 2 P3 O7
BOND 17 17 19 9 P3 O8
BOND 18 17 20 1 P3 O9
BOND 19 17 21 1 P3 O10
BOND 20 21 22 1 O10 C2
BOND 21 22 23 1 C2 H7
BOND 22 22 24 1 C2 H8
BOND 23 22 25 1 C2 H9
BOND 24 27 28 1 C3 H10
BOND 25 27 29 1 C3 H11
BOND 26 27 30 1 C3 H12
BOND 27 27 31 1 C3 C4
BOND 28 31 32 2 C4 O11
BOND 29 31 33 1 C4 N2
BOND 30 33 34 1 N2 H13
BOND 31 33 35 1 N2 H14
BOND 32 36 37 1 O12 H15
BOND 33 36 38 1 O12 H16
BOND 34 39 40 1 C5 H17
BOND 35 39 41 1 C5 H18
BOND 36 39 42 1 C5 H19
BOND 37 39 43 1 C5 C6
BOND 38 43 44 9 C6 O13
BOND 39 43 45 9 C6 O14

respin1
Resp charges for organic molecule

 &cntrl

 nmol = 1,
 ihfree = 1,
 ioutopt = 1,

 &end
    1.0
Resp charges for organic molecule
   -2 45
    6 0
    1 0
    1 0
    1 0
    7 0
    1 0
    1 6
    1 6
   15 0
    8 0
    8 10
    8 10
   15 0
    8 0
    8 14
    8 0
   15 0
    8 0
    8 18
    8 0
    8 0
    6 0
    1 0
    1 0
    1 0
   12 0
    6 0
    1 0
    1 0
    1 0
    6 0
    8 0
    7 0
    1 0
    1 34
    8 0
    1 0
    1 37
    6 0
    1 0
    1 0
    1 0
    6 0
    8 0
    8 44

respin2
Resp charges for organic molecule

 &cntrl

 nmol = 1,
 ihfree = 1,
 ioutopt = 1,
 iqopt = 2,
 qwt = 0.001,

 &end
    1.0
Resp charges for organic molecule
   -2 45
    6 0
    1 0
    1 2
    1 2
    7 -99
    1 -99
    1 -99
    1 -99
   15 -99
    8 -99
    8 -99
    8 -99
   15 -99
    8 -99
    8 -99
    8 -99
   15 -99
    8 -99
    8 -99
    8 -99
    8 -99
    6 0
    1 0
    1 23
    1 23
   12 -99
    6 0
    1 0
    1 28
    1 28
    6 -99
    8 -99
    7 -99
    1 -99
    1 -99
    8 -99
    1 -99
    1 -99
    6 0
    1 0
    1 40
    1 40
    6 -99
    8 -99
    8 -99

Resp charges for organic molecule

          Point Charges Before & After Optimization

    no. At.no. q(init) q(opt) ivary d(rstr)/dq
    1 6 0.000000 -0.002792 0 0.004998
    2 1 0.000000 0.031850 0 0.000000
    3 1 0.000000 0.115049 0 0.000000
    4 1 0.000000 0.055049 0 0.000000
    5 7 0.000000 -0.210419 0 0.002146
    6 1 0.000000 0.295469 0 0.000000
    7 1 0.000000 0.295469 6 0.000000
    8 1 0.000000 0.295469 6 0.000000
    9 15 0.000000 1.242491 0 0.000401
   10 8 0.000000 -0.874121 0 0.000568
   11 8 0.000000 -0.874121 10 0.000568
   12 8 0.000000 -0.874121 10 0.000568
   13 15 0.000000 1.295660 0 0.000385
   14 8 0.000000 -0.808153 0 0.000614
   15 8 0.000000 -0.808153 14 0.000614
   16 8 0.000000 -0.504444 0 0.000972
   17 15 0.000000 1.367824 0 0.000365
   18 8 0.000000 -0.838262 0 0.000592
   19 8 0.000000 -0.838262 18 0.000592
   20 8 0.000000 -0.547496 0 0.000898
   21 8 0.000000 -0.529450 0 0.000928
   22 6 0.000000 0.008818 0 0.004981
   23 1 0.000000 0.073118 0 0.000000
   24 1 0.000000 0.056336 0 0.000000
   25 1 0.000000 0.093358 0 0.000000
   26 12 0.000000 1.298719 0 0.000384
   27 6 0.000000 -0.509647 0 0.000963
   28 1 0.000000 0.177629 0 0.000000
   29 1 0.000000 0.153961 0 0.000000
   30 1 0.000000 0.091579 0 0.000000
   31 6 0.000000 0.983691 0 0.000506
   32 8 0.000000 -0.776470 0 0.000639
   33 7 0.000000 -0.950083 0 0.000523
   34 1 0.000000 0.429671 0 0.000000
   35 1 0.000000 0.429671 34 0.000000
   36 8 0.000000 -0.744648 0 0.000665
   37 1 0.000000 0.365327 0 0.000000
   38 1 0.000000 0.365327 37 0.000000
   39 6 0.000000 -0.469241 0 0.001042
   40 1 0.000000 0.057543 0 0.000000
   41 1 0.000000 0.052121 0 0.000000
   42 1 0.000000 0.153304 0 0.000000
   43 6 0.000000 0.941899 0 0.000528
   44 8 0.000000 -0.783262 0 0.000633
   45 8 0.000000 -0.783262 44 0.000633

 Sum over the calculated charges: -2.000




Junmei Wang wrote:

> In this situation, prepgen can't work properly. You may solve this problem
> as the following: generate an ac file from gaussian output file; generate
> resp input file with respgen and do modifications (you want to the metal to
> bear an integer charge); run resp to get resp charges and read in to the ac
> file with antechamber; For each residue, extract the corresponding part
> from the ac file (you may ignore the bond connection information) to form a
> mpdb file; for each residue, read in the mpdb file and generate the prep
> input file with antechamber. This procedure should work for you.
>
> Best
>
> Junmei
>
> ===============================================================
> Dr. Junmei Wang
> Chemistry & Biophysics
> Encysive Pharmaceuticals
> 7000 Fannin, Houston TX 77030
> Tel: 713-5786649
> Email: jwang.tbc.com
> Homepage: Http://sigyn.compchem.ucsf.edu/members/jmwang/index.html
> ===============================================================
>
>
>


-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Wed Jan 14 2004 - 15:53:11 PST
Custom Search