AMBER Archive (2007)

Subject: Re: AMBER: Lone Pairs

From: David A. Case (case_at_scripps.edu)
Date: Mon Feb 12 2007 - 15:24:48 CST


On Mon, Feb 12, 2007, Beale, John wrote:

> How are lone pairs handled in AMBER 8.0 (specifically in FF03)?

There are no lone pairs in FF03.

> I am
> trying to input a peptide into LEaP. The structure contains one
> methionine. The pdb file that I am loading was saved through SYBYL, and
> two lone pairs are present in the pdb file, designated as "LPD". . When
> I do a loadPdb in LEaP the program tells me that two lone pairs were
> added to the structure.

This makes sense: your pdb file has lone pairs in it, but the FF03 library
files don't know anything about lone pairs. Hence, LEaP thinks it needs
to "add" these to the structure. When the time comes to save a prmtop
file, LEaP should complain that it doesn't have enough information about
the lone-pairs to create a force field for them. (I'm guessing here:
you just said the "LEaP doesn't save the prmtop file", without any
indication of what the actual message is.)

There are some Amber force fields that include lone pairs (look for
leaprc files with "EP" in their filenames.) However, they are all quite
old, and are not really recommeded for any routine use. [Yong or Piotr
might want to chime in here: leaprc.ff02polEP.r1 is in fact "recent",
but I don't believe that it has been extensively tested.]

...good luck...dac
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu