Re: [AMBER] RESP charges

From: Sadaf Rani <sadafrani6.gmail.com>
Date: Thu, 11 Jul 2019 16:49:25 +0100

Dear Kellon
thank you
it makes sence for mol2 file but not for prepi
as I did following:-

antechamber -fi gout -fo prepi -i .gout -o bg6.prepi -c resp -at gaff -rn
bg6
Here it is different. If you please could have a look at these 2 files.
If I want to use respgen program for the output of gaussian is it the same
file(gout) to convert fisrt in to .ac file and then run in to respgen for 2
stage fitting?

thank you


On Thu, Jul 11, 2019 at 3:19 PM Kellon Belfon <kellonbelfon.gmail.com>
wrote:

> Hi Sadaf,
>
> The second column in the resp.out file represents the atomic number of the
> element. The index of the first column in the resp.out file correspond to
> the first column in your mol2 file. For example, first index is 1 which has
> an atomic number of 15 (resp.out file) and hence it is a P (mol2 file). You
> will replace the charge 1.352900 in the last column of your mol2 file for
> index 1 with the optimized charge of 1.304538 from your resp fitting.
> 1 P 3.8390 123.5530 13.4480 p5 1 LIG 1.352900
> will become:
> 1 P 3.8390 123.5530 13.4480 p5 1 LIG 1.30453
>
> Hope this is clear. Let me know if you have any other questions.
>
> Regards,
>
> Kellon
>
> On Thu, Jul 11, 2019 at 8:37 AM Sadaf Rani <sadafrani6.gmail.com> wrote:
>
> > Hi David and AMBER
> > I am still confused as how to recognise atoms in resp.out and add charges
> > accordingly?
> > as format of resp.out is:-
> > Resp charges for organic molecule
> >
> >
> > Point Charges Before & After Optimization
> >
> > no. At.no. q(init) q(opt) ivary d(rstr)/dq
> > 1 15 0.000000 1.304538 0 0.000382
> > 2 8 0.000000 -0.937716 0 0.000530
> > 3 8 0.000000 -0.937716 2 0.000530
> > 4 8 0.000000 -0.937716 2 0.000530
> > 5 6 0.000000 0.518773 0 0.000946
> > 6 8 0.000000 -0.705468 0 0.000702
> > 7 6 0.000000 0.238018 0 0.001937
> > 8 8 0.000000 -0.710723 0 0.000697
> > 9 6 0.000000 0.245035 0 0.001889
> > 10 8 0.000000 -0.734824 0 0.000674
> > 11 6 0.000000 0.104696 0 0.003454
> >
> >
> > and for mol2:-
> >
> > .<TRIPOS>ATOM
> > 1 P 3.8390 123.5530 13.4480 p5 1 LIG
> > 1.352900
> > 2 O1P 3.1070 124.2090 14.5470 o 1 LIG
> > -0.747167
> > 3 O2P 5.1500 124.1140 13.0540 o 1 LIG
> > -0.747167
> > 4 O3P 3.7560 122.0800 13.4850 o 1 LIG
> > -0.747167
> > 5 C1 0.9050 125.5120 9.0780 c3 1 LIG
> > 0.355900
> > 6 O1 1.5480 125.2890 7.8690 oh 1 LIG
> > -0.627800
> > 7 C2 -0.5830 125.2840 8.9230 c3 1 LIG
> > 0.090100
> > 8 O2 -1.0930 126.1730 7.9430 oh 1 LIG
> > -0.631800
> > 9 C3 -1.2720 125.5090 10.2640 c3 1 LIG
> > 0.015100
> > 10 O3 -2.6840 125.2920 10.1360 oh 1 LIG
> > -0.665800
> > 11 C4 -0.6670 124.5400 11.2730 c3 1 LIG
> > 0.283100
> >
> >
> > Also I am look on this tutorial
> > http://ambermd.org/tutorials/advanced/tutorial1/section1.htm but from
> > where
> > resp.in file is comming in step iv)?
> > Any one could you please suggest me.
> >
> > thank you
> >
> > Sadaf
> >
> > On Tue, Jul 9, 2019 at 2:28 PM David Case <david.case.rutgers.edu>
> wrote:
> >
> > > On Sun, Jul 07, 2019, Sadaf Rani wrote:
> > >
> > > >I calculated charges for my ligand
> > > >As per para 10 they have asked to edit atomic charges in prepi file
> from
> > > >resp.out file. However I have confusion in recognising atom numbers in
> > > resp
> > > >out and in my prepi file.
> > >
> > > It may be easier to ask antechamber to make a mol2 file (rather than
> > > prepi), and edit the charges there. But you'll still have to be sure
> > > that the atom ordering in the resp and mol2 files are the same.
> > >
> > > ...dac
> > >
> > >
> > > _______________________________________________
> > > AMBER mailing list
> > > AMBER.ambermd.org
> > > http://lists.ambermd.org/mailman/listinfo/amber
> > >
> > _______________________________________________
> > AMBER mailing list
> > AMBER.ambermd.org
> > http://lists.ambermd.org/mailman/listinfo/amber
> >
>
>
> --
> Kellon A. A. Belfon, Graduate Student
> Carlos Simmerling Laboratory
> The Laufer Center for Physical and Quantitative Biology
> The Department of Chemistry, Stony Brook University
> Stony Brook, New York 11794
> Phone: (347) 546-4237 <(347)+546+4237> Email: kellon.belfon.stonybrook.
> <kellon.belfon.stonybrook.edu>edu
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Jul 11 2019 - 09:00:02 PDT
Custom Search