AMBER: RESP of negatively charged molecule

From: Rhoad, Jonathan S. <rhoadj.missouri.edu>
Date: Tue, 1 Mar 2005 15:39:05 -0600

I was attempting to obtain RESP charges for glucose-6-phosphate. I have
included the method and optimized geometry below. The structure was
optimized at mp2/6-31+G* with PCM (without PCM a proton migrated to the
phosphate). When I use antechamber or the resp(gen) programs
individually, I get the same nonsense charges. The prep file is after
the Gaussian input file.

 

antechamber -i output.log -fi gout -o g6p.prep -fo prepi -rn G6P -at
amber -c resp

 

where output.log is the output of the Gaussian file included.

 

Any suggestions? Would the ESP charges output in the Gaussian file be
good enough? I am trying to simulate the dynamic motion of sugar
phosphates in the active site of enzymes.

 

Thanks,

Jonathan Rhoad

University of Missouri

 

 

# hf/6-31+G(d) sp pop=mk scf=tight IOp(6/33=2)

 

-2 1

C -2.218496 -1.514421 0.250925

C -2.841025 -0.236926 0.810923

C -2.456779 0.955423 -0.051291

C -0.948398 0.992892 -0.228695

C -0.429008 -0.326598 -0.783213

C 1.074543 -0.354423 -0.931119

O -2.860871 -1.763714 -0.980495

O -4.247791 -0.385538 0.970549

O -2.922062 2.136351 0.601125

O -0.658237 2.074725 -1.122705

O -0.809237 -1.390773 0.116253

O 1.692802 0.083376 0.267180

P 3.404618 -0.066814 0.273953

O 3.873914 0.726611 -0.958070

O 3.668239 -1.580765 0.194876

O 3.740833 0.577286 1.628032

H -2.363305 -2.350348 0.948609

H -2.438877 -0.073120 1.819327

H -2.927300 0.855711 -1.044632

H -0.489311 1.167197 0.753295

H -0.879028 -0.504575 -1.771736

H 1.386325 -1.379133 -1.177377

H 1.356781 0.294608 -1.771237

H -2.617198 -2.682830 -1.264005

H -4.617197 -0.572358 0.075266

H -2.609697 2.891499 0.049849

H 0.223790 2.439407 -0.865811

 

    0 0 2

 

This is a remark line

molecule.res

G6P XYZ 0

CORRECT OMIT DU BEG

  0.0000

   1 DUMM DU M 0 -1 -2 0.000 .0 .0 .00000

   2 DUMM DU M 1 0 -1 1.449 .0 .0 .00000

   3 DUMM DU M 2 1 0 1.522 111.1 .0 .00000

   4 O4 OH M 3 2 1 1.540 111.208 180.000 1.224

   5 H11 HO E 4 3 2 0.993 152.287 -148.100 0.430

   6 C4 CT M 4 3 2 1.436 78.469 -39.722 -1.121

   7 H4 H1 E 6 4 3 1.099 110.233 -122.639 2.167

   8 C3 CT M 6 4 3 1.521 106.979 -4.616 -13.618

   9 O3 OH S 8 6 4 1.432 111.076 -65.289 2.625

  10 H10 HO E 9 8 6 0.991 107.756 58.034 1.309

  11 H3 H1 E 8 6 4 1.103 109.191 55.801 5.178

  12 C2 CT M 8 6 4 1.523 109.444 176.132 -2.363

  13 O2 OH S 12 8 6 1.428 113.121 -177.859 0.955

  14 H9 HO E 13 12 8 0.992 108.656 63.148 2.084

  15 H2 H1 E 12 8 6 1.099 108.379 65.602 3.678

  16 C1 CT M 12 8 6 1.529 110.650 -52.563 -12.046

  17 O1 OH S 16 12 8 1.413 106.989 -68.254 0.310

  18 H8 HO E 17 16 12 0.996 108.491 -167.564 2.299

  19 H1 H2 E 16 12 8 1.098 110.642 170.717 4.529

  20 O5 OS M 16 12 8 1.424 111.048 54.580 3.759

  21 C5 CT M 20 16 12 1.446 113.172 -59.830 -0.233

  22 H5 H1 E 21 20 16 1.101 109.345 -58.210 1.798

  23 C6 CT M 21 20 16 1.512 107.769 -175.474 -5.644

  24 H6 H1 E 23 21 20 1.099 108.599 49.872 1.605

  25 H7 H1 E 23 21 20 1.099 108.881 167.584 1.605

  26 O6 OS M 23 21 20 1.426 109.774 -71.486 0.182

  27 P1 P M 26 23 21 1.710 116.411 173.721 -5.788

  28 O8 O E 27 26 23 1.550 105.802 -64.142 1.691

  29 O9 O E 27 26 23 1.548 101.227 176.266 1.691

  30 O7 O M 27 26 23 1.551 105.534 56.542 1.691

 

 

LOOP

   C5 C4

 

IMPROPER

 

DONE

STOP


-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Tue Mar 01 2005 - 21:53:01 PST
Custom Search