Re: [AMBER] problem with tleap

From: Rejwan <rejwan_88.yahoo.com>
Date: Wed, 19 Jun 2013 08:06:05 -0700 (PDT)

Hi

Thanks for your input. The G09 input file I used to extract RESP via Antechamber reads as following. Should I make changes and re-run Gaussian in the following line?

#HF/6-31G* SCF=tight Test Pop=MK iop(6/33=2) iop(6/42=6) opt iop(6/50=1)

Thanks for your suggestions.

//Rejwan


--Link1--
%chk=molecule
#HF/6-31G* SCF=tight Test Pop=MK iop(6/33=2) iop(6/42=6) opt iop(6/50=1)

remark line goes here

0 1
    C 6.6215000000 5.1461000000 0.1163000000
    C 5.7259000000 4.2389000000 -0.3176000000
    C 6.2373000000 6.3368000000 0.9704000000
--------------------------------
 H 11.0017000000 6.8296000000 2.6926000000
    H 10.5984000000 6.6920000000 -1.6004000000
    H 11.7998000000 9.1604000000 2.5323000000
    H 11.3920000000 9.0243000000 -1.7406000000
    H 11.9987000000 10.2632000000 0.3195000000

g09.gesp

g09.gesp






--- On Tue, 6/18/13, Rejwan wrote:

&gt; From: Rejwan
&gt; Subject: Re: [AMBER] problem with tleap
&gt; To: &quot;AMBER Mailing List&quot;
&gt; Cc: jason.swails.gmail.com
&gt; Date: Tuesday, June 18, 2013, 9:46 PM
&gt; --- On Tue, 6/18/13, Jason
&gt; Swails  wrote:
&gt;
&gt; &gt; From: Jason Swails
&gt; &gt; Subject: Re: [AMBER] problem with tleap
&gt; &gt; To: &quot;AMBER Mailing List&quot;
&gt; &gt; Date: Tuesday, June 18, 2013, 5:09 PM
&gt; &gt; On Tue, Jun 18, 2013 at 4:06 PM,
&gt; &gt; Rejwan
&gt; &gt; wrote:
&gt; &gt;
&gt; &gt; &gt; I&#39;m trying to make the topology file
&gt; of a 
&gt; &gt; ligand-protein system.
&gt; &gt; &gt;
&gt; &gt; &gt; In tLeap suit of Amber I fed both .prepin
&gt; and .frcmod
&gt; &gt; of the ligand,  I
&gt; &gt; &gt; used RESP charge obtained via Gaussian 09
&gt; Geometry
&gt; &gt; Optimization for PREPIN
&gt; &gt; &gt; file. The .frcmod file was created by
&gt; parmchk routine.
&gt; &gt; &gt;
&gt; &gt;
&gt; &gt; Did you look at the resulting frcmod file?  Did
&gt; you see
&gt; &gt; any statements like
&gt; &gt; [ATTN: Needs Revision]?
&gt;
&gt; I LOOKED AT FRCMOD FILE. IT DOES NOT HAVE ANY WARNING AT
&gt; ALL.
&gt;
&gt; &gt;
&gt; &gt;
&gt; &gt; &gt; The tLeap error showed only 2 out of 10
&gt; torsion terms
&gt; &gt; could be
&gt; &gt; &gt; implemented. Could you please help how can
&gt; I extract
&gt; &gt; the missing Angle and
&gt; &gt; &gt; Torsion terms ? Also where should I plug
&gt; the terms (Is
&gt; &gt; it in  the .frmod
&gt; &gt; &gt; file ?) to rerun Leap ?  I&#39;m copying
&gt; the error
&gt; &gt; screen.
&gt; &gt; &gt;
&gt; &gt;
&gt; &gt; You should put any missing parameters into the
&gt; frcmod
&gt; &gt; file.  However,
&gt; &gt; parmchk should have done this for you (unless there
&gt; were
&gt; &gt; bonds between your
&gt; &gt; ligand and protein, for instance, but this seems
&gt; unlikely
&gt; &gt; given that all
&gt; &gt; missing parameters involve gaff atom types).  Even
&gt; if
&gt; &gt; parmchk couldn&#39;t find
&gt; &gt; any parameters, it will put zeroes in for
&gt; everything and
&gt; &gt; indicate that it
&gt; &gt; needs revision. (You should always check the frcmod
&gt; file
&gt; &gt; after parmchk
&gt; &gt; creates it).  Therefore, the fact that parameters
&gt; are
&gt; &gt; still missing
&gt; &gt; indicates that you&#39;ve probably made some kind
&gt; of error (did
&gt; &gt; you remember
&gt; &gt; to load the frcmod into your tleap session via the
&gt; [
&gt; &gt; loadAmberParams ]
&gt; &gt; command?)
&gt;
&gt; PARMCHK DID NOT HAVE ANY COMPLAIN AND OF COURSE I DID IMPORT
&gt; FRCMOD INTO TLEAP VIA loadAmberParams
&gt; &gt;
&gt; __________________
&gt;
&gt;
&gt; IF YOU WANT TO HAVE A LOOK AT THOSE FILES I CAN SEND THEM.
&gt; HOWEVER, WHAT ELSE COULD BE THE PITFALL ? INTERESTINGLY I
&gt; RUN LIGEAND-PROTEIN  MD  BEFORE BY SAME
&gt; ANTECHAMBER PROTOCOL WITHOUT PROBLEM. THANKS.
&gt;
&gt; -rejwan
&gt;
&gt;
&gt; &gt;
&gt; &gt; --
&gt; &gt; Jason M. Swails
&gt; &gt; Quantum Theory Project,
&gt; &gt; University of Florida
&gt; &gt; Ph.D. Candidate
&gt; &gt; 352-392-4032
&gt; &gt; _______________________________________________
&gt; &gt; AMBER mailing list
&gt; &gt; AMBER.ambermd.org
&gt; &gt; http://lists.ambermd.org/mailman/listinfo/amber
&gt; &gt;
&gt;

_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed Jun 19 2013 - 08:30:02 PDT
Custom Search