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
Received on Fri Oct 21 2016 - 07:30:03 PDT
