AMBER: Torsion terms in xleap

From: Lili Peng <lpeng.ucsd.edu>
Date: Fri, 14 Sep 2007 15:35:23 -0700

Hi AMBER users,

I have a quick question about preparing a mol2 file in antechamber and
uploading it into xleap. I currently have a mol2 file (see attachment) and
here are the commands I'm running:

1. Generate frcmod file using antechamber: /parmchk -i glu4_linear_PTX.mol2
-f mol2 -o glu4_linear_PTX.frcmod
2. Upload the mol2 file into xleap:
    source leaprc.gaff
    loadamberparams glu4_linear_PTX.frcmod
    loadamberparams gaff.dat
    loadamberparams parm99.dat
    x = loadmol2 glu4_linear_PTX.mol2
    solvatebox x TIP3PBOX 10
    saveamberparm x glu4_linear_PTX.prmtop glu4_linear_PTX.inpcrd

....But then I get the error messages:
1-4: angle 30 35 duplicates bond ('triangular' bond) or angle ('square'
bond)
1-4: angle 31 34 duplicates bond ('triangular' bond) or angle ('square'
bond)
** No torsion terms for H-C.3-C.3-O.3
** No torsion terms for H-C.3-C.3-H
** No torsion terms for H-C.3-C.3-H
** No torsion terms for O.3-C.3-O.3-H

What is the problem? A colleague of mine suggested that I manually add in
the torsion terms into gaff.dat. How should I go about doing this?

If you have any insight on this problem, I would appreciate it.

Thanks,
Lili



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

Received on Sun Sep 16 2007 - 06:07:43 PDT
Custom Search