AMBER Archive (2009)

Subject: Re: [AMBER] Re: COOH troubling

From: Jio M (jiomm_at_yahoo.com)
Date: Sat Dec 12 2009 - 00:27:25 CST


>I am not sure how you did QM calculations and how

>you set the radii for Iodine. I calculated HF/CEP-121G (the radii for I was

>set to 1.98) optimization for your ligand molecule and derived the resp

>charges

I used this for lig.mol2 that I provided as attached file: 
hf/lanl2dz SCF=Tight Pop=MK IOp(6/33=2,6/41=10,6/42=17)
geom=connectivity pop=(mk,readrad)

Iodine radius also set to 1.98 Angstrom and derived RESP with antechamber.

The attached lig2.mol2  was having charges provided by your method
HF/CEP-121G (pls correct me if I am wrong). I was able to run dynamics
not only of single ligand but also in docked state by loading your
lig2.mol2 file. But I did parmchk on lig2.mol2 to get frcmod file. I
used .in file as given in min6.out and md7.out.

But when I use my conditions, same thing happened , COOH trouble and
ring distortion.Using my conditions even single molecule dynamics is
crashing.

this is input file:

&cntrl

  imin = 0,

  ntb = 2, pres0 = 1.0, ntp = 1,

  taup = 2.0,

  cut = 8.0, ntr = 0,

  tempi = 0.0, temp0 = 0.0,

  ntt = 3, gamma_ln = 1.0,

  nstlim = 100000, dt = 0.0001,

  ntpr = 100, ntwx = 100, ntwr = 1000

 /

But by USING your input file and charges I was succesfull always.

It would add to my knowledge if you can please share how to decide and
become confident that HF/CEP-121G or other basis will fit. Does it
depends on relative RMS SQRT while doing RESP calculations.

One problem still hanging, in starting few trajectories I can see that
COOH group is cis but after some ps it becomes trans. Does the force
constant of CO--OH bond rotation is weaker.

>You also need to question why the article used the protonated carboxyl and

>protonated phenolic form.

I checked cross ref. of the experimental procedure they are following
to simulate. The release pattern was studied by dialysis with solvents
which I think donot have pH of 2. there is no clear cut pH conditions
mentioned. May be I am still missing something to understand why
protonated structure was used.

>You need to seriously consider the pH at which you want to
simulate, because that will >completely determine which of the three
protonation forms you need to parameterize

Yes, honestly, I was not so much serious before, regarding such
structures. Now I will consider the things. Thanks for giving such a
detailed information for my ligand.

>Good luck with your work,

I am highly thankfull to Dr. Dean, Dr. Case and Dr. Junmei, for the explanation

thanks and regards;

JIomm

      
_______________________________________________
AMBER mailing list
AMBER_at_ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber