<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=utf-8"><meta name=Generator content="Microsoft Word 12 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
{font-family:Helvetica;
panose-1:2 11 6 4 2 2 2 2 2 4;}
@font-face
{font-family:Helvetica;
panose-1:2 11 6 4 2 2 2 2 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
@font-face
{font-family:Consolas;
panose-1:2 11 6 9 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman","serif";
color:black;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
pre
{mso-style-priority:99;
mso-style-link:"HTML Preformatted Char";
margin:0cm;
margin-bottom:.0001pt;
font-size:10.0pt;
font-family:"Courier New";
color:black;}
span.HTMLPreformattedChar
{mso-style-name:"HTML Preformatted Char";
mso-style-priority:99;
mso-style-link:"HTML Preformatted";
font-family:Consolas;
color:black;}
span.EmailStyle19
{mso-style-type:personal-reply;
font-family:"Calibri","sans-serif";
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body bgcolor=white lang=EN-US link=blue vlink=purple><div class=WordSection1><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>It seems to me that it will be for the reasons mentioned in the paper <a href="http://dx.doi.org/10.1063/1.2393240">http://dx.doi.org/<span style='font-size:12.0pt;font-family:"Times New Roman","serif"'>10.1063/1.2393240</span></a></span><o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>By focusing on the phosphate atom, you are measuring the actual movement of the entire lipid molecule within the membrane, which is the diffusion that is measured typically experimentally. If you use the entire molecule, then you get the movement of the alkane chains as well, which increase the diffusion coefficient by more than a factor of ten.<span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><div><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>Catch ya,<br><br>Dr. Dallas Warren<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>Medicinal Chemistry and Drug Action<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>Monash Institute of Pharmaceutical Sciences</span><span style='font-size:10.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>, Monash University<br>381 Royal Parade, Parkville VIC 3010<br>dallas.warren@monash.edu<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>+61 3 9903 9304<br>---------------------------------<br>When the only tool you own is a hammer, every problem begins to resemble a nail.</span><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'> <o:p></o:p></span></p></div><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><div style='border:none;border-left:solid blue 1.5pt;padding:0cm 0cm 0cm 4.0pt'><div><div style='border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm'><p class=MsoNormal><b><span style='font-size:10.0pt;font-family:"Tahoma","sans-serif";color:windowtext'>From:</span></b><span style='font-size:10.0pt;font-family:"Tahoma","sans-serif";color:windowtext'> gmx-users-bounces@gromacs.org [mailto:gmx-users-bounces@gromacs.org] <b>On Behalf Of </b>Javier Cerezo<br><b>Sent:</b> Friday, 3 December 2010 2:46 AM<br><b>To:</b> gmx-users@gromacs.org<br><b>Subject:</b> Re: [gmx-users] Re: individual lateral diffusion coefficients<o:p></o:p></span></p></div></div><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Sorry for the misspell...<br><br><o:p></o:p></p><p class=MsoNormal>Thanks Justin. <br><br>Do you <b>know</b> the reason behind? <br><br>I am trying following that protocol and my P curve in not as linear as All-lipid-atoms one (DPPC, FF: G53a6, t=50ns). Furthermore, regarding the linear region, the slopes are not the same. So which one do you think is more accurate? <br><br>Thanks again <br><br>Javier <br><br>El 02/12/10 16:09, Justin A. Lemkul escribió: <br><br><o:p></o:p></p><p class=MsoNormal><br><br>Ángel Piñeiro wrote: <br><br><o:p></o:p></p><p class=MsoNormal style='margin-bottom:12.0pt'>Hi Javier <br>I think I saw this in several mails of this list and it is also implicit in the Justin tutorial for analysis of bilayers. I am not sure whether or not this has also been published... I do not remember any paper. I think this is reasonable for lipids in contact with membrane proteins because only a part of the lipid could be "tied" to the protein... then the diffusion for different parts of the lipid could be different. <br><br>I think I will calculate the diffusion for each lipid individually... <br><br>If you are interested in comparing results you could contact me off the list. <o:p></o:p></p><p class=MsoNormal><br>I haven't followed this thread at all, but I saw my name come up :) This is what I have always based my g_msd calculations on: <br><br><a href="http://lists.gromacs.org/pipermail/gmx-users/2008-January/031804.html">http://lists.gromacs.org/pipermail/gmx-users/2008-January/031804.html</a> <br><br>-Justin <br><br><br><o:p></o:p></p><p class=MsoNormal>Saludos, <br><br>Ángel. <br><br><br><br><br>On Thu, 2010-12-02 at 14:22 +0100, Javier Cerezo wrote: <br><br><o:p></o:p></p><p class=MsoNormal>Hi Ángel <br><br>Can you provide a citation about the use of only PO4 atoms to calculate the diffusion constant? Is it always recommended or just with CG simulations? I'm also working on diffusion calculation and that will be interesting. <br><br>By the way, regarding the index files I mentioned, it might be better to have a group of lipids that are at a certain distance from the protein in the same index, again to improve sampling (maybe this is not a way to improve sampling, I don't know). <br><br>Thanks <br><br>Javier <br><br><br>El 02/12/10 13:56, Ángel Piñeiro escribió: <br><br><o:p></o:p></p><p class=MsoNormal>Hi Javier <br>1.- you are right! the diff_mol.xvg file I reported was from a previous attempt in which I used the whole lipid molecules with the -mol option on, instead of the PO4 beads with -mol off. Sorry for this confusion <br><br>2.- As I said above, I did attempts using both the whole molecule and the PO4 beads. Yes I saw the figure 6 in the Wolhert &Edholm's paper but I read in several other references that the calculation is more accurate by using only the P atom... what makes sense to me mainly for the lipids which are in contact with proteins <br><br>3.- I agree that removing jumps does not change anything. I decided to give this information in my message to avoid a reply saying "try to remove jumps" ;) <br><br>4.- Yes I agree that I could do the calculation by creating an index for each lipid... I guess that is the safest way to proceed... <br><br>Thanks for your reply! <br><br>Ángel. <br><br><br><br>On Thu, 2010-12-02 at 13:30 +0100, Javier Cerezo wrote: <br><br><o:p></o:p></p><p class=MsoNormal>Hello Ángel. <br><br>Well, there are a some things that I don't understand about your calculation, but might be just a problem of mine. Here you have my comments: <br><br>1. How do you get the diff_mol.xvg file if you are not using -mol in your command line input (and you index file has broken molecules). <br><br>2. Why do you select just an atom to calculate the diffusion? According to Wolhert and Edholm (JCP, 125, 204703) the MSD for all lipids atoms reach the same slope, so I guess using them all could improve sampling (I'm not sure). <br><br>3. I think that reprocessing of your trajectory to remove jumps is no longer needed (I got the same results in a recent test using version 4.5.1). <br><br>4. What I would do to calculate D as funtion to the distance to the membrane protein is generate different index files containing lipids according to this distance (and hoping they don't move a lot during the simulation) and run different msd calculations. I think I have read in the mailing list about a script to make such a selection regarding distances to construct the index file or just make your own one. <br><br>Good luck <br><br>Javier <br>El 02/12/10 12:50, Ángel Piñeiro escribió: <br><br><o:p></o:p></p><p class=MsoNormal>I want to add that the MSD as a function of time (msd.xvg file) looks completely linear <br><br>Greetings, <br><br>Ángel Piñeiro. <br><br>On Thu, 2010-12-02 at 12:45 +0100, Ángel Piñeiro wrote: <br><br><o:p></o:p></p><p class=MsoNormal style='margin-bottom:12.0pt'>Dear all, <br>I aim to calculate the lateral diffusion coefficients of lipids as a function of the distance to a membrane protein using the Martini force field. For this I guess I could use the diff_mol.xvg output file of the g_msd command which provides the list of diffusion coefficients for each lipid (I guess the lipids are ordered as in the trajectory file). Then I would calculate the protein-lipid distance for each lipid and I would generate the diffusion vs distance file. Before starting the calculations on the membrane protein system I tested the g_msd command on a DPPC bilayer. In my bilayer simulation I removed the COM of lipids and water separately. Before analyzing it I removed jumps over the box boundaries using trjconv -pbc nojump and I created a index file with the PO4 atoms as a new group. Then I executed the following command: <br><br>g_msd -s topol.tpr -f trajnojump.xtc -n p.ndx -lateral z -rmcomm <br><br>from which I get the following output: <br>D[ PO4] 0.0958 (+/- 0.0135) 1e-5 cm^2/s <br><br>I think the value is not crazy for DPPC at 323 K using Martini... but I noticed that the D values for the independent lipids reported in the diff_mol.xvg file range from 0.0021959 to 0.482909 cm^2/s. If the differences are so high for a single lipid bilayer I suspect that I will not observe significant differences as a function of the distance to the protein in my simulations of the whole system... probably I am doing something wrong¿? <br><br>Thanks for any advice <br><br>Ángel Piñeiro. <o:p></o:p></p><p class=MsoNormal><br>-- <br>Javier CEREZO BASTIDA <br>Estudiante de Doctorado <br>--------------------- <br>Dpto. Química-Física <br>Universidad de Murcia <br>30100 MURCIA (España) <br>Tlf.(+34)868887434 <br>-- <br>gmx-users mailing list <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a> <a href="mailto:gmx-users@gromacs.org"><mailto:gmx-users@gromacs.org></a> <br><a 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 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 href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a> <a href="mailto:gmx-users-request@gromacs.org"><mailto:gmx-users-request@gromacs.org></a>. <br>Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists">http://www.gromacs.org/Support/Mailing_Lists</a> <o:p></o:p></p><p class=MsoNormal><br>-- <br>Javier CEREZO BASTIDA <br>Estudiante de Doctorado <br>--------------------- <br>Dpto. Química-Física <br>Universidad de Murcia <br>30100 MURCIA (España) <br>Tlf.(+34)868887434 <br>-- <br>gmx-users mailing list <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a> <a href="mailto:gmx-users@gromacs.org"><mailto:gmx-users@gromacs.org></a> <br><a 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 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 www interface or send it to <a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a> <a href="mailto:gmx-users-request@gromacs.org"><mailto:gmx-users-request@gromacs.org></a>. <br>Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists">http://www.gromacs.org/Support/Mailing_Lists</a> <o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><br><br><o:p></o:p></p><pre>-- <o:p></o:p></pre><pre>Javier CEREZO BASTIDA<o:p></o:p></pre><pre>Estudiante de Doctorado<o:p></o:p></pre><pre>---------------------<o:p></o:p></pre><pre>Dpto. Química-Física<o:p></o:p></pre><pre>Universidad de Murcia<o:p></o:p></pre><pre>30100 MURCIA (España)<o:p></o:p></pre><pre>Tlf.(+34)868887434<o:p></o:p></pre></div></div></body></html>