Justin,<br><br>thanks again for help.<br><br>Finally is there any generall solution to parametrise hetero-groups covalently bonded with the protein ? Many proteins consist of such groups e.g chromophore in GFP, retinall in rhodopsin as well as some prostetic groups in the enzymes.<br>
<br>I&#39;ve tried to make something like you&#39;ve told me via inclusion of pre-parametrised residues in the existing gromacs ff but forced with some problems due to the atom order in new ITP and gro files provided by ATb or PRODRG are different from initial pdb file so pdb2gmx on the whole protein where het-group in the old order would not work properly :(<br>
<br>James<br><br><div class="gmail_quote">2012/6/10 Justin A. Lemkul <span dir="ltr">&lt;<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>&gt;</span><br><blockquote class="gmail_quote" style="margin:0pt 0pt 0pt 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div class="im"><br>
<br>
On 6/10/12 3:06 AM, James Starlight wrote:<br>
<blockquote class="gmail_quote" style="margin:0pt 0pt 0pt 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
Justin<br>
<br>
I have one extra question about parametrisation of the bond type<br>
<br>
Initially I had -c-c- bond so both atoms were in the sp3 hybridization. I&#39;ve<br>
changed bond type between that atoms to the gb_16 ( from phe ring) to define<br>
them as the &gt;c=c&lt; but during the simulation I&#39;ve noticed that the bond between<br>
that atoms is not planar ( sp2 form) and still is in the sp3. What addition<br>
changees should I do in my topology to change the planarity of new bond ?<br>
<br>
<br>
</blockquote>
<br></div>
Defining appropriate atom types will be necessary here, as well.  For a single bond, the atom types would be -CH2-CH2- while a double bond will necessarily be -CH=CH-<br>
<br>
Also note that with a united atom representation, it&#39;s very hard to make such assertions about hybridization state.  There are no hydrogens, so the geometry of these central carbons is not easily defined.<div class="im">
<br>
<br>
<blockquote class="gmail_quote" style="margin:0pt 0pt 0pt 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
Should I make some changes ia angles term as well?<br>
<br>
E.g I&#39;ve found<br>
<br>
   612   613   614     2    ga_13  ;<br>
<br>
where 612 and 613 are the atoms wich I would like to present as the rigid. I&#39;ve<br>
change to the ga_27 ( 120) but during simulation the bond was still in sp3. What<br>
I&#39;ve done wrong ?<br>
<br>
</blockquote>
<br></div>
As there are no parameters for linear C=C bonds in the Gromos force fields, you may have to do some parameterization yourself, particularly for appropriate dihedrals.  One option is to look at how these bonds are treated in lipid force fields like the Berger parameter set, distributed by Peter Tieleman&#39;s group. There are double bonds defined for lipids like POPC.<span class="HOEnZb"><font color="#888888"><br>

<br>
-Justin</font></span><div class="HOEnZb"><div class="h5"><br>
<br>
-- <br>
==============================<u></u>==========<br>
<br>
Justin A. Lemkul, Ph.D.<br>
Research Scientist<br>
Department of Biochemistry<br>
Virginia Tech<br>
Blacksburg, VA<br>
jalemkul[at]<a href="http://vt.edu" target="_blank">vt.edu</a> | (540) 231-9080<br>
<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.<u></u>vt.edu/Pages/Personal/justin</a><br>
<br>
==============================<u></u>==========<br>
-- <br>
gmx-users mailing list    <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/<u></u>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/<u></u>Support/Mailing_Lists/Search</a> before posting!<br>
Please don&#39;t post (un)subscribe requests to the list. Use the www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can&#39;t post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists</a><br>
</div></div></blockquote></div><br>