AMBER Archive (2009)

Subject: Re: [AMBER] Re: COOH troubling

From: Dean Cuebas (deancuebas_at_missouristate.edu)
Date: Wed Dec 09 2009 - 18:00:40 CST


Dear Jiomm

Your structure is definitely incorrect. You cannot have a carbonyl (C=O) as
part of an aromatic ring system.

Either get rid of the double bond and add a negative charge to make a
phenolate ion, or make a C-O-H by adding a hydrogen to the oxygen to make a
phenolic group, like you have on the opposite side of your anthracene ring
system.

-- 
Dr. Dean Cuebas, Associate Prof of Chemistry
deancuebas_at_missouristate.edu, Ph 417-836-8567 FAX 417-836-5507
Dept. of Chemistry, Missouri State University
Springfield, Missouri 65897

> From: Junmei Wang <junmwang_at_gmail.com> > Reply-To: AMBER Mailing List <amber_at_ambermd.org> > Date: Wed, 9 Dec 2009 17:17:12 -0600 > To: AMBER Mailing List <amber_at_ambermd.org> > Subject: Re: [AMBER] Re: COOH troubling > > Hi, > It seems to me that your ligand (lig.mol2) has some problem. If the net > charge of this molecule is not -1, then it should be a radical according to > chemists' intuition. What is the spin multiplicity when you did QM > calculation? > > Best > > Junmei > > On Wed, Dec 9, 2009 at 6:16 AM, Jio M <jiomm_at_yahoo.com> wrote: > >> Dear AMBER users; >> >> >> >> I am attaching files of ligand mol2 (lig.mol2) and its frcmod (lig.frcmod). >> and solvent mol2 (dcm.mol2) and frcmod as dcm.frcmod >> >> >> >> Now I tried with dynamics of single ligand with dcm solvent. I >> found that same problem of COOH atoms is occuring. I have used GAFF. >> >> I solvated like this (with single solvent molecule) >> >> solvateoct ligand dcm 10 >> >> approx. 180 dcm were added. >> >> >> >> 1) Using prmtop and inpcrd file I minimised first with restraints on >> ligand only. It minimised without error as restraints were applied on >> solute so atoms of COOH did not trouble. >> >> 2) But after removing constraints, error came indicating atoms of COOH. >> >> >> >> 3) Then I tried directly NVP dynamics using rst file of restraint >> minimisation (1 step) and this input: >> >> &cntrl >> >> imin = 0, >> >> ntb = 2, pres0 = 1.0, ntp = 1, >> >> taup = 2.0, >> >> cut = 10.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 >> >> / >> >> >> >> here I can see in trajectory that >> >> a) the hydrogen of COOH is approaching towards carbonyl oxygen of COOH. >> >> >> >> b) ALSO all the rings are distorted and I think rings are behaving as >> Cyclohexane. >> >> >> >> I think some parameters are missing? Please help and suggest which >> parameters I am missing here. >> >> >> >> thanks and regards; >> >> JIomm >> >> >> >> >> >> >> >> _______________________________________________ >> AMBER mailing list >> AMBER_at_ambermd.org >> http://lists.ambermd.org/mailman/listinfo/amber >> >> > _______________________________________________ > AMBER mailing list > AMBER_at_ambermd.org > http://lists.ambermd.org/mailman/listinfo/amber

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