Hi,<br><br>I already have both your publication and supplementary material from your website. The files were also submitted to Gromacs user contribution section:<br><a href="http://www.gromacs.org/@api/deki/files/23/=LipidsForGro96_53a6.zip">http://www.gromacs.org/@api/deki/files/23/=LipidsForGro96_53a6.zip</a><br>
Have you noticed anything suspicious in alpha helix behaviour in G53a6 force field?<br><br>Christopher<br><br><div class="gmail_quote">2010/1/21 Kukol, Andreas <span dir="ltr">&lt;<a href="mailto:a.kukol@herts.ac.uk">a.kukol@herts.ac.uk</a>&gt;</span><br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Hi,<br>
<br>
Lipid parameters for Gromos96 53a6 are available based on the publication<br>
<br>
Kukol, A., 2009. Lipid models for united-atom molecular dynamics simulations of proteins. J. Chem. Theor. Comput.,5(3), 615-626.<br>
<br>
Unfortunately, I did not manage yet to upload it on the new Gromacs website. Please contact me off the list for the itp-files for DPPC, DMPC, POPC and POPG for 53a6 forcefield.<br>
<font color="#888888"><br>
Andreas<br>
</font><div><div></div><div class="h5"><br>
<br>
&gt; -----Original Message-----<br>
&gt; From: <a href="mailto:gmx-users-bounces@gromacs.org">gmx-users-bounces@gromacs.org</a> [mailto:<a href="mailto:gmx-users-bounces@gromacs.org">gmx-users-bounces@gromacs.org</a>]<br>
&gt; On Behalf Of Justin A. Lemkul<br>
&gt; Sent: 21 January 2010 01:36<br>
&gt; To: Discussion list for GROMACS users<br>
&gt; Subject: Re: [gmx-users] Lipid parameters for GROMOS96 force fields<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; Justin A. Lemkul wrote:<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt; Krzysztof Mlynarczyk wrote:<br>
&gt; &gt;&gt; Hello,<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; On Peter Tieleman&#39;s website we can find topologies for various popular<br>
&gt; &gt;&gt; lipids. However those seem to be valid for the deprecated GROMOS87<br>
&gt; &gt;&gt; force field. The g96_lipids.itp file is still unavailable. I&#39;d like to<br>
&gt; &gt;&gt; run a simulation of 7TM protein embedded into POPC membrane. G43a2x<br>
&gt; &gt;&gt; seems to be a reasonable choice for this purpose. I still need the<br>
&gt; &gt;&gt; parameters for my lipids. I am aware that combining parameters from<br>
&gt; &gt;&gt; different force fields might lead to a disaster, that&#39;s why I want to<br>
&gt; &gt;&gt; ask you:<br>
&gt; &gt;&gt; 1. Do you know any papers where POPC membrane and G43a2 force field<br>
&gt; &gt;&gt; was used?<br>
&gt; &gt;<br>
&gt; &gt; I&#39;ve not seen one, because it has been known for some time that early<br>
&gt; &gt; Gromos96 parameter sets did not reproduce long alkyl chain properties<br>
&gt; &gt; particularly well.<br>
&gt; &gt;<br>
&gt;<br>
&gt; Quick amendment - modifications to 43a2 and 45a3 have been made to better<br>
&gt; reproduce lipid properties.  Both are available for download on the Gromacs<br>
&gt; site<br>
&gt; and you can probably track down the references, as well.<br>
&gt;<br>
&gt; Sorry, read that one wrong the first time :)<br>
&gt;<br>
&gt; -Justin<br>
&gt;<br>
&gt; &gt;&gt; 2. If not, is there any way to derive the proper parameters for the<br>
&gt; &gt;&gt; force field of my choice using the lipid parameters from Peter<br>
&gt; &gt;&gt; Tieleman&#39;s website or e.g. the parameters published by Andreas Kukol<br>
&gt; &gt;&gt; for G53a6?<br>
&gt; &gt;&gt;<br>
&gt; &gt;<br>
&gt; &gt; I don&#39;t see why you need to do such reverse engineering.  The Kukol<br>
&gt; &gt; parameters for lipids under 53a6 can be directly combined with a G53a6<br>
&gt; &gt; protein without any issues; I believe that was the purpose of the whole<br>
&gt; &gt; new derivation :)<br>
&gt; &gt;<br>
&gt; &gt; As an aside, you are quite right that multiple force fields within the<br>
&gt; &gt; same simulation is incorrect.  However, the Berger lipid parameters may<br>
&gt; &gt; be an exception to this rule, since they are really a hybridized version<br>
&gt; &gt; of OPLS-UA and Gromos87 parameters (some of which were modified anyway),<br>
&gt; &gt; so they really don&#39;t belong to any one particular force field.  The<br>
&gt; &gt; Berger/G87 combination is widely used, but essentially amounts to the<br>
&gt; &gt; following: lipid interactions are Berger-Berger or OPLS-OPLS<br>
&gt; &gt; interactions, while protein-lipid interations are Berger-G87, and<br>
&gt; &gt; protein-protein interactions are G87-G87.  You can see quite quickly why<br>
&gt; &gt; things become complicated!<br>
&gt; &gt;<br>
&gt; &gt; Based on a discussion I had with Dr. Tieleman, it seems to be reasonable<br>
&gt; &gt; to use the G96 parameter set of your choice in conjunction with<br>
&gt; &gt; lipid.itp (Berger lipids), although other approaches may be more<br>
&gt; &gt; rigorously correct (pure G96 parameters such as those by Kukol, pure<br>
&gt; &gt; OPLS recently derived by Ulmschneider, or the modifications to the<br>
&gt; &gt; Berger parameters from the Tieleman group, to name a few).  If you want<br>
&gt; &gt; to use a G96-lipid.itp combination, I created a tutorial that teaches<br>
&gt; &gt; you how to build the system and properly prepare the topology.  It is<br>
&gt; &gt; linked from the Tutorials page of the Gromacs site.<br>
&gt; &gt;<br>
&gt; &gt; -Justin<br>
&gt; &gt;<br>
&gt; &gt;&gt; Thank you in advance,<br>
&gt; &gt;&gt; Christopher<br>
&gt; &gt;&gt;<br>
&gt; &gt;<br>
&gt;<br>
&gt; --<br>
&gt; ========================================<br>
&gt;<br>
&gt; Justin A. Lemkul<br>
&gt; Ph.D. Candidate<br>
&gt; ICTAS Doctoral Scholar<br>
&gt; MILES-IGERT Trainee<br>
&gt; Department of Biochemistry<br>
&gt; Virginia Tech<br>
&gt; Blacksburg, VA<br>
&gt; jalemkul[at]<a href="http://vt.edu" target="_blank">vt.edu</a> | (540) 231-9080<br>
&gt; <a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br>
&gt;<br>
&gt; ========================================<br>
&gt; --<br>
&gt; gmx-users mailing list    <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
&gt; <a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
&gt; Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before posting!<br>
&gt; Please don&#39;t post (un)subscribe requests to the list. Use the<br>
&gt; www interface or send it to <a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>
&gt; Can&#39;t post? Read <a href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>
--<br>
gmx-users mailing list    <a href="mailto: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/search" target="_blank">http://www.gromacs.org/search</a> before posting!<br>
Please don&#39;t post (un)subscribe requests to the list. Use the<br>
www interface or send it to <a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>
Can&#39;t post? Read <a href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>
</div></div></blockquote></div><br>