Thanks David,
You're right, it's from Leap, but I mixed things because now in Sleap,
'antechamber' is called automatically (set default fastbld on).
Anyway, I maybe still need something more robust. It's because I have a DB
with almost 9000 ligands and acpypi/antechamber/mopac worked for ~75%, but I
have several inputs there that doesn't have a good initial structure or
Corina didn't do a good job.
Alan.
On Thu, Jun 4, 2009 at 16:06, David A. Case <case.biomaps.rutgers.edu>wrote:
> On Thu, Jun 04, 2009, Alan wrote:
>
> > But I still want to understand better this "Antechamber actually computes
> a
> > minimized structure" assertion because I would like to know why I got
> this
> > msg sometimes when running antechamber:
> >
> > Warning: Close contact of 1.932973 angstroms between .R<MOL 1>.A<C7 8>
> and
> > .R<MOL 1>.A<H1 21>
> >
> > I got this warning after 'minimisation'? I definitely would like to find
> a
> > way of calculating myself QM/semi-empirical minimisation. As you
> mentioned,
> > I will try to find how to do it with mopac explicitly.
>
> As I said earlier, antechamber uses the AM1-minimized structure to get
> the am1-bcc charges, but then discards that, and uses the original input
> coordinates in the output mol2 or prepin files. The message you quote
> above looks like it comes from LEaP, not from antechamber.
>
> I've often thought that we should have an option to output
> the minimized structure -- I'll cc this to junmei for comment.
>
> ...dac
>
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>
--
Alan Wilter S. da Silva, D.Sc. - CCPN Research Associate
Department of Biochemistry, University of Cambridge.
80 Tennis Court Road, Cambridge CB2 1GA, UK.
>>http://www.bio.cam.ac.uk/~awd28<<
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Fri Jun 05 2009 - 01:11:47 PDT
