<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 TRANSITIONAL//EN">
<HTML>
<HEAD>
  <META HTTP-EQUIV="Content-Type" CONTENT="text/html; CHARSET=UTF-8">
  <META NAME="GENERATOR" CONTENT="GtkHTML/3.28.1">
</HEAD>
<BODY>
Thanks again Justin<BR>
<BR>
-----Original Message-----<BR>
<B>From</B>: Justin A. Lemkul &lt;<A HREF="mailto:%22Justin%20A.%20Lemkul%22%20%3cjalemkul@vt.edu%3e">jalemkul@vt.edu</A>&gt;<BR>
<B>Reply-to</B>: &quot;jalemkul@vt.edu&quot; &lt;jalemkul@vt.edu&gt;, Discussion list for GROMACS users &lt;gmx-users@gromacs.org&gt;<BR>
<B>To</B>: Gromacs Users' List &lt;<A HREF="mailto:Gromacs%20Users'%20List%20%3cgmx-users@gromacs.org%3e">gmx-users@gromacs.org</A>&gt;<BR>
<B>Subject</B>: Re: [gmx-users] g_density CHARMM<BR>
<B>Date</B>: Tue, 11 Jan 2011 22:08:22 +0100<BR>
<BR>
<PRE>

Ramon Guix&#224; wrote:
&gt; Justin,
&gt; Thanks for your response.
&gt; The membranes I've already run are made of POPC and cholesterol, so the 
&gt; sphingomyelin membrane is my next simulation but I haven't managed to 
&gt; find a proper forcefield for it. Would you know where to find this?

Not for CHARMM, no, but there are several available for different Gromos force 
fields.

<A HREF="http://www.gromacs.org/Downloads/User_contributions/Molecule_topologies">http://www.gromacs.org/Downloads/User_contributions/Molecule_topologies</A>
<A HREF="http://www.apmaths.uwo.ca/~mkarttu/downloads.shtml">http://www.apmaths.uwo.ca/~mkarttu/downloads.shtml</A> (contact Prof. Karttunen for 
the topology)

-Justin

&gt; Ramon
&gt; 
&gt; -----Original Message-----
&gt; *From*: Justin A. Lemkul &lt;<A HREF="mailto:jalemkul@vt.edu">jalemkul@vt.edu</A> 
&gt; &lt;<A HREF="mailto:%22Justin%20A.%20Lemkul%22%20%3cjalemkul@vt.edu">mailto:%22Justin%20A.%20Lemkul%22%20%3cjalemkul@vt.edu</A>%3e&gt;&gt;
&gt; *Reply-to*: &quot;<A HREF="mailto:jalemkul@vt.edu">jalemkul@vt.edu</A>&quot; &lt;<A HREF="mailto:jalemkul@vt.edu">jalemkul@vt.edu</A>&gt;, Discussion list for 
&gt; GROMACS users &lt;<A HREF="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</A>&gt;
&gt; *To*: Discussion list for GROMACS users &lt;<A HREF="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</A> 
&gt; &lt;<A HREF="mailto:Discussion%20list%20for%20GROMACS%20users%20%3cgmx-users@gromacs.org">mailto:Discussion%20list%20for%20GROMACS%20users%20%3cgmx-users@gromacs.org</A>%3e&gt;&gt;
&gt; *Subject*: Re: [gmx-users] g_density CHARMM
&gt; *Date*: Tue, 11 Jan 2011 19:45:58 +0100
&gt; 
&gt; 
&gt; Ramon Guix&#224; wrote:
&gt;&gt; 
&gt;&gt; Dear all,
&gt;&gt; 
&gt;&gt; I am just a newbie to Gromacs so sorry for being naive with my questions...
&gt;&gt; 
&gt;&gt; I would like to yield electron density profiles out of different 
&gt;&gt; trajectories of membranes simulated with NAMD and last CHARMM forcefield 
&gt;&gt; (all_36_lipid_cholesterol). I first tried to use pdb2gmx in order to 
&gt;&gt; generate a GROMACS topology file that I could use for generating the 
&gt;&gt; .tpr file needed for g_density. Well, I did not manage but is this the 
&gt;&gt; proper way though?
&gt; 
&gt; Not likely.  pdb2gmx is useful for creating topologies of macromolecules built 
&gt; of repeating units (i.e. proteins and polymers).  For lipid bilayers, it is far 
&gt; more common to construct the topology as a set of #include statements, i.e.:
&gt; 
&gt; #include (forcefield)
&gt; 
&gt; #include &quot;first_lipid_type.itp&quot;
&gt; #include &quot;next_lipid_type.itp&quot;
&gt; 
&gt; etc.
&gt; 
&gt;&gt; Also, so far all my trajectories are .dcd files; does the g_density 
&gt;&gt; function needs a .pdb file per frame? Otherwise, is it possible to 
&gt;&gt; directly convert .dcd files into the gromacs .xtc ones?
&gt;&gt; 
&gt; 
&gt; Newer versions of Gromacs (&gt;4.5) can read any file format that VMD can read, so 
&gt; conversion is not necessary, provided you have the VMD libraries installed.
&gt; 
&gt;&gt; On another note, would anyone know where could I get proper topology and 
&gt;&gt; parameter files of sphingomyelin?
&gt;&gt; 
&gt; 
&gt; If you already ran the simulations, do you not already have this?  With the 
&gt; information in Chapter 5 of the manual, you should be able to convert your 
&gt; topologies from whatever format you have it in into the one required by Gromacs.
&gt; 
&gt; -Justin
&gt; 
&gt;&gt; Many thanks in advance!
&gt;&gt; 
&gt;&gt; Ramon
&gt;&gt; 
&gt; 
&gt; -- 
&gt; ========================================
&gt; 
&gt; Justin A. Lemkul
&gt; Ph.D. Candidate
&gt; ICTAS Doctoral Scholar
&gt; MILES-IGERT Trainee
&gt; Department of Biochemistry
&gt; Virginia Tech
&gt; Blacksburg, VA
&gt; jalemkul[at]vt.edu | (540) 231-9080
&gt; <A HREF="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</A>
&gt; 
&gt; ========================================
&gt; 

-- 
========================================

Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
MILES-IGERT Trainee
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
<A HREF="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</A>

========================================
</PRE>
</BODY>
</HTML>