We gave you advice before but it seems you keep repeating the same error
and not trying to find out why it is happening. You'll need to look into
your leap logs and see what it says. Sometimes doing research simulations
takes some troubleshooting and testing on your part.
To summarize - making a coordinate file in another program that had a bond
is only part of the process. You also need to explain to leap that it is
cyclic, which prof case already described. If you follow these directions
and still have problems, look into the logs and see if you can find out
why, or tell us why leap is not doing what is expected from the commands
you gave.
On Oct 21, 2016 10:04 AM, "Thakur, Abhishek" <axt651.miami.edu> wrote:
> Hi I have managed to make it cyclic by making bond in gaussian.
>
> But as soon as I am solvating it and making prmtop, inpcrd and pdb it
> seems that N terminal is breaking bond and making it -NH3
>
>
> I don't know why it is creating so much of trouble.
>
> ________________________________
> From: Bill Ross <ross.cgl.ucsf.edu>
> Sent: Friday, October 21, 2016 1:43:38 AM
> To: david.case.rutgers.edu; AMBER Mailing List
> Subject: Re: [AMBER] frcmod file error
>
>
> >> 3. generate a reasonable starting structure for the cyclized peptide.
> > Amber is not a very good tool (generally speaking) for carrying out step
> 3.
> >
> Though not the best at it, xleap will let you construct a cyclic
> peptide. In the cmd line:
>
> > mymol = sequence { ALA ALA .. ALA } # not a chemist btw
> > edit mymol
>
> And there get used to the controls so you can select the end atoms of
> the linear molecule one by one and and crudely move each out of line, so
> that you can draw a long bond between those atoms. Then select all the
> atoms, relax selection, all the weird bonds will be fixed as the
> molecule reforms, and with luck you are done, though iterative
> selection/relaxation might be needed if kinks develop. The energy
> minimization function is based on atom hybridizations, and lacks vdw.
>
> In the end you might want to have a more powerful and used molecular
> editor program at your fingertips, but xleap is there and once upon a
> time we put some time into making it work.
>
> Bill
>
> _______________________________________________
> 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=HLgnadgQYVqBGkQSIdb5YMS4Szk5sHPJcjLLme_cFVs&s=
> JwaD4ttp4K2xCyXKfEBi7g9Z0uuv0qcVQ5ouK1Baqgk&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 Fri Oct 21 2016 - 07:30:04 PDT
