Thank you Justin for your suggestion. I searched into the archive to
find something about GFP chromophore before writing to the list, but I
did not find anything, I' d try to search better.
Anna
______________________________________________
Anna Marabotti, Ph.D.
Assistant Professor
Department of Chemistry and Biology
University of Salerno
Via Ponte don Melillo
84084 Fisciano (SA)
Italy
Phone: +39 089 969583
Fax: +39 089 969603
E-mail: [email protected]
Skype: annam1972
"When a man with a gun meets a man with a pen, the man with the gun is a dead
man"
(Roberto Benigni, about Roberto Saviano)
Message: 3
Date: Wed, 06 Mar 2013 07:09:58 -0500
From: Justin Lemkul<[email protected]>
Subject: Re: [gmx-users] problem with a complex cofactor
To: Discussion list for GROMACS users<[email protected]>
Message-ID:<[email protected]>
Content-Type: text/plain; charset=ISO-8859-15; format=flowed
On 3/6/13 6:29 AM, Anna Marabotti wrote:
Dear gmx-users,
I'm working on a protein of the family of green fluorescent proteins, which has
a cofactor formed by the rearrangement of three residues that cyclized forming a
complex heterocyclic fluorescent probe. This probe is covalently bound to the
protein and is part of the backbone; in the PDB file it is marked as HETATM in
the middle of the protein sequence which is: ...S61-H62-V63-CFY66-H68...
I really do not know how to deal with the topology of this residue and with the
fact that the protein sequence is interrupted. I don't know if I have to manage
this residue as a new residue to be added in the force field; please take into
account that this residue is bound to V63 and H68 with covalent bonds similar,
but not strictly identical, to peptide bonds. I found in the gmx-users archive a
suggestion about how to manage the topology of a Lysine bound to pyridoxal
phosphate, but that situation is different from mine, since Lysine is a
"regular" residue "regularly" forming peptide bonds with its adjacent residues,
and which is attached to PLP with its terminal side chain NZ atom, whereas in my
case a heterocyclic molecules is forming peptide-like bonds with adjacent
residues. Moreover, mine is not the case of a N-ter or a C-ter residue somehow
modified: this heterocyclic compound is in the middle of protein sequence.
Can somebody suggest me a strategy on how to proceed?
Weird bonds can be assigned with specbond.dat, but you may be able to do
everything within the .rtp file, anyway, since you're still dealing with a
linear polypeptide. Usually specbond.dat is only necessary for branch points or
other complex structures. Please note that we have discussed the GFP
chromophore at length several times, even recently, so there is bound to be
relevant information in the list archive.
-Justin
--
gmx-users mailing list [email protected]
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the
www interface or send it to [email protected].
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
