Re: AMBER: Problem with leaprc

From: yen li <chem_me2000.yahoo.com>
Date: Wed, 4 Jul 2007 05:23:49 -0700 (PDT)

Dear Dr. Case,
As suggested I have done 'rdparm->dihedral' analysis for a tripeptide. Below I am giving the first two non hydrogen dihedral details. (Using leaprc.ff03)
 
      113: 2.500 3.14 2.0 :1.O :1.C :2.N :2.CA (19,18,20,22)
      114: 0.778 3.14 1.0 :1.C :2.N :2.CA :2.CB (18,20,22,24)
E 115: 0.066 3.14 2.0 :1.C :2.N :2.CA :2.CB (18,20,22,24)
E 116: 0.056 0.00 3.0 :1.C :2.N :2.CA :2.CB (18,20,22,24)
E 117: 0.000 3.14 4.0 :1.C :2.N :2.CA :2.CB (18,20,22,24)


(a) first dihedral is O -C -N -CT or CT-N -C -O, I am unable to find such dihed parameters in parm99.dat or frcmod.ff03
(b) 2nd, 3rd, 4th dihedrals are C -N -CT-CT but the above parameters are of CT-CT-C -N of frcmod.ff03.

Is this a bug?

I am attaching the corresponding pdb file.

Regards,
Yen


"David A. Case" <case.scripps.edu> wrote: On Sun, Jul 01, 2007, yen li wrote:

> (2) I do "desc ff03"

I've never done this, so I'm not sure what to expect. To really know what is
being computed, set up a small molecule, and "dihed" command in rdparm to
print out the dihedrals from the prmtop file. That is what sander and pmemd
really see.

...dac

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu


       
---------------------------------
Looking for a deal? Find great prices on flights and hotels with Yahoo! FareChase.

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu

Received on Sun Jul 08 2007 - 06:07:07 PDT
Custom Search