I find the things you are saying to be confusing. You sent a pdb file
(2.pdb) that has only normal amino acids and should work fine, but then in
your email you don't show that this is the pdb file that you used for Leap.
Also this is just a peptide, so why would you save the files as "complex"?
What error do you get if you try to use only the pdb file that you sent? If
no error, then what is the pdb file that gives you the error? I still don't
see anything in what you sent that explains why you are trying to use gaff.
On Thu, Oct 20, 2016 at 10:46 AM, Thakur, Abhishek <axt651.miami.edu> wrote:
> yes I have tried that
>
>
> What I did was
>
> > > source leaprc.gaff
> > > loadamberparams PLX4720.frcmod
> > > solvateBox complex TIP3PBOX 10
> > > charge complex
> > > addIons complex Cl- 7
> > > saveamberparm complex V600E_4720.prmtop V600E_4720.inpcrd
>
>
> and it gave me following error
>
> For your reference I have attached my pdb file
>
>
> That is why I was trying to make frcmod file as it is asking me those
> parameters
>
>
> ________________________________
> From: Carlos Simmerling <carlos.simmerling.gmail.com>
> Sent: Thursday, October 20, 2016 3:40:39 AM
> To: AMBER Mailing List
> Subject: Re: [AMBER] frcmod file error
>
> if it is standard amino acids, I'm not sure I understand why you are using
> a frcmod file.. Did you try without it?
>
> On Thu, Oct 20, 2016 at 10:37 AM, Thakur, Abhishek <axt651.miami.edu>
> wrote:
>
> > Hey Dr.David,
> >
> >
> > Thank you for your suggestion, I think I have a cyclic peptide of around
> > 75 atoms, so I think that is why antechamber is not working.
> >
> >
> > So I tried making it as protein is prepared. As I just want to simulate
> > this peptide by solvating it, so I have not frcmod, lib files of
> heteroatom
> > to load.
> >
> > So as these are peptide amber should build it as a protein? No hetero
> atom
> > only small chain of 9 amino acid in cyclic form.
> >
> >
> > > source leaprc.gaff
> > > loadamberparams PLX4720.frcmod
> > > solvateBox complex TIP3PBOX 10
> > > charge complex
> > > addIons complex Cl- 7
> > > saveamberparm complex V600E_4720.prmtop V600E_4720.inpcrd
> >
> >
> > The error that I am getting is
> >
> >
> > Could not find bond parameter for: N3 - C
> > Building angle parameters.
> > Could not find angle parameter: CX - N3 - C
> > Could not find angle parameter: H - N3 - C
> > Could not find angle parameter: H - N3 - C
> > Could not find angle parameter: H - N3 - C
> > Could not find angle parameter: N3 - C - CX
> > Could not find angle parameter: N3 - C - O2
> > Could not find angle parameter: N3 - C - O2
> > Building proper torsion parameters.
> > ** No torsion terms for CX-N3-C-CX
> > ** No torsion terms for CX-N3-C-O2
> > ** No torsion terms for CX-N3-C-O2
> > ** No torsion terms for H-N3-C-CX
> > ** No torsion terms for H-N3-C-O2
> > ** No torsion terms for H-N3-C-O2
> > ** No torsion terms for H-N3-C-CX
> > ** No torsion terms for H-N3-C-O2
> > ** No torsion terms for H-N3-C-O2
> > ** No torsion terms for H-N3-C-CX
> > ** No torsion terms for H-N3-C-O2
> > ** No torsion terms for H-N3-C-O2
> >
> >
> > As antechamber doe not work so I am unable to generate frcmod file
> > containing this files.
> >
> > But as these are just amino acid do I need to give frcmod file or lib for
> > these things?
> >
> > If yes then how do I generate it as antechamber is not creating a mol2
> > file for my peptide and amber gives error if I use mol2 file generated by
> > some other software and make frcmod file.
> >
> > Unknown keyword: C.3-C.3-S. in parameter file. Perhaps a format error?
> >
> >
> >
> > Thanking you,
> >
> > Abhsihek
> >
> >
> > ________________________________
> > From: David A Case <david.case.rutgers.edu>
> > Sent: Wednesday, October 19, 2016 1:06:36 PM
> > To: AMBER Mailing List
> > Subject: Re: [AMBER] frcmod file error
> >
> > On Wed, Oct 19, 2016, Thakur, Abhishek wrote:
> > >
> > > > antechamber -i PLX4720.pdb -fi pdb -o PLX4720.mol2 -fo mol2 -c bcc
> -s 2
> > > But It was continuously giving me error
> >
> > > ---Judge bond type for Residue 1 with ID of 1 and Name of GLU ---
> > >
> > > ---Judge bond type for Residue 2 with ID of 2 and Name of ASN ---
> > >
> > > Warning: the assigned bond types may be wrong, please :
> > > (1) double check the structure (the connectivity) and/or
> > > (2) adjust atom valence penalty parameters in APS.DAT, and/or
> > > (3) increase PSCUTOFF in define.h and recompile bondtype.c
> > > Be cautious, use a large value of PSCUTOFF (>100) will
> significantly
> > increase the computation time
> >
> > Note: antechamber is intended for *small molecules*, not for polypeptides
> > or
> > proteins. Amber builds proteins in a completely different way: please
> see
> > the
> > manual or tutorials.
> >
> > [Scott: does acdoctor properly complain if a user tries to use a protein
> > pdb
> > file?]
> >
> > ...good luck...dac
> >
> >
> > _______________________________________________
> > AMBER mailing list
> > AMBER.ambermd.org
> > https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.
> > ambermd.org_mailman_listinfo_amber&d=DQICAg&c=y2w-uYmhgFWijp_IQN0DhA&r=
> > 0Hah93XWYYBgmROOVcaccg&m=tzk1nJYUywLBe_8SN5Q89lVTcSAFz0940JNz3RzCCPs&
> > s=lDeec35G74H2q_6oa6VVDM9N2Za-76ygiTtgT0UnC5c&e=
> > _______________________________________________
> > AMBER mailing list
> > AMBER.ambermd.org
> > https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.
> ambermd.org_mailman_listinfo_amber&d=DQICAg&c=y2w-uYmhgFWijp_IQN0DhA&r=
> 0Hah93XWYYBgmROOVcaccg&m=jbzG8RFnibcn6VYCd16zuRwF8ZQKF_
> brz7oG1zbD7Jc&s=Pec7P_U0dS67_QGlaK_7WrvmaUwplPhrtsDmcCBCnVc&e=
> >
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.
> ambermd.org_mailman_listinfo_amber&d=DQICAg&c=y2w-uYmhgFWijp_IQN0DhA&r=
> 0Hah93XWYYBgmROOVcaccg&m=jbzG8RFnibcn6VYCd16zuRwF8ZQKF_
> brz7oG1zbD7Jc&s=Pec7P_U0dS67_QGlaK_7WrvmaUwplPhrtsDmcCBCnVc&e=
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>
>
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Oct 20 2016 - 08:00:05 PDT