<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;">Hi Rama,<br><br> I was looking at the tdb files for the amber03 force field ( I didn't change anything, i just read with the cat command ), anyways I didn't see anything written in the tdb -n or -c terminal files. After reading the gromacs manual on topology files, it said that residue connection types and such are specified through the tdb file. I'm only guessing that since the chromophore isn't an amino acid that it may have a different connection deviating from the traditional -n or -c terminals. So I was wondering if the amber force field was using the .tdb files that comes standard with gromacs ,... which I'm pretty sure lacks the type of connection that pyp makes with the protein. So to fix this problem, do we need to make or modify the .tdb file for the chromophore, or is the .tdb file of no concern? <br><br>Thanks,<br><br>Taylor<br><br>---
On <b>Wed, 4/20/11, Ramachandran G <i><gtrama@gmail.com></i></b> wrote:<br><blockquote style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; padding-left: 5px;"><br>From: Ramachandran G <gtrama@gmail.com><br>Subject: Re: [gmx-users] PYP chromophore force field<br>To: "Peter C. Lai" <pcl@uab.edu><br>Cc: "gmx-users@gromacs.org" <gmx-users@gromacs.org><br>Date: Wednesday, April 20, 2011, 5:12 PM<br><br><div id="yiv1713269628"><br>Yes, you are right, the chromophore(p-coumaric acid) is covalently bonded to the Sulphur(S) atom of the Cys-69 residue.<br>But I specified the covalent bond in my topology.<br><br>[ bonds]<br> C +SG<br><br>[ angles ]<br>
O C +SG <br> CA1 C +SG <br> C +SG +CB <br><br>[ dihedrals ]<br> O C +SG +CB<br> CA1 C +SG +CB<br> HA2 CA1 C +SG<br> CA2 CA1 C +SG<br> C +SG +CB +CA<br>
C +SG +CB +HB1<br> C +SG +CB +HB2<br><br>After doing pdb2gmx, i got the right structure. But the chromophore flies away after minimising. <br><br><br>Thanks<br>Rama<br><br><div class="yiv1713269628gmail_quote">On Wed, Apr 20, 2011 at 4:57 PM, Peter C. Lai <span dir="ltr"><<a rel="nofollow" ymailto="mailto:pcl@uab.edu" target="_blank" href="/mc/compose?to=pcl@uab.edu">pcl@uab.edu</a>></span> wrote:<br>
<blockquote class="yiv1713269628gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Isn't the chromophore supposed to be covalently bonded to the protein?<br>
My cursory search shows Vengris et al 2004 in Biophys. J. citing Baca et al<br>
1994 and van Beeumen et al 1993: "This small 125-residue protein contains a<br>
p-coumaric acid chromophore that is covalently bound to the protein backbone<br>
via a thiol-ester cysteine linkage Cys-69"...<br>
<br>
Maybe you forgot to specify a covalent bond in your topology...<br>
<div><div></div><div class="yiv1713269628h5"><br>
On 2011-04-20 06:46:54PM -0500, Ramachandran G wrote:<br>
> Hi gromacs users,<br>
> I am working on Photo active yellow protein.<br>
><br>
> Although i successfully build the force field and patched the chromophore to the protein.<br>
> After energy minimizing the protein, the chromophore flies away separately. I don't know whether i am missing anything?<br>
> Please help.<br>
> Rama<br>
><br>
><br>
<br>
</div></div>> --<br>
> gmx-users mailing list <a rel="nofollow" ymailto="mailto:gmx-users@gromacs.org" target="_blank" href="/mc/compose?to=gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
> <a rel="nofollow" target="_blank" href="http://lists.gromacs.org/mailman/listinfo/gmx-users">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
> Please search the archive at <a rel="nofollow" target="_blank" href="http://www.gromacs.org/Support/Mailing_Lists/Search">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
> Please don't post (un)subscribe requests to the list. Use the<br>
> www interface or send it to <a rel="nofollow" ymailto="mailto:gmx-users-request@gromacs.org" target="_blank" href="/mc/compose?to=gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>
> Can't post? Read <a rel="nofollow" target="_blank" href="http://www.gromacs.org/Support/Mailing_Lists">http://www.gromacs.org/Support/Mailing_Lists</a><br>
<font color="#888888"><br>
<br>
--<br>
===============================================================<br>
Peter C. Lai | University of Alabama-Birmingham<br>
Programmer/Analyst | BEC 257<br>
Genetics, Div. of Research | 1150 10th Avenue South<br>
<a rel="nofollow" ymailto="mailto:pcl@uab.edu" target="_blank" href="/mc/compose?to=pcl@uab.edu">pcl@uab.edu</a> | Birmingham AL 35294-4461<br>
<a rel="nofollow">(205) 690-0808</a> |<br>
===============================================================<br>
<br>
</font></blockquote></div><br><br clear="all"><br><br>
</div><br>-----Inline Attachment Follows-----<br><br><div class="plainMail">-- <br>gmx-users mailing list <a ymailto="mailto:gmx-users@gromacs.org" href="/mc/compose?to=gmx-users@gromacs.org">gmx-users@gromacs.org</a><br><a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>Please don't post (un)subscribe requests to the list. Use the <br>www interface or send it to <a ymailto="mailto:gmx-users-request@gromacs.org" href="/mc/compose?to=gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a></div></blockquote></td></tr></table>