<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;"><br>Thanx Justin,<br><br>I am somewhat faraway from the topic so can't remember clearly how to define. may be you are right. But this the case where i stuck before.<br><br><br><br>regards,<br><div><br><span style="color: rgb(56, 83, 118); font-family: comic sans ms;">Ravi Datta Sharma</span></div>  <div><span style="color: rgb(56, 83, 118); font-family: comic sans ms;">Lecturer,</span></div>  <div><span style="color: rgb(56, 83, 118); font-family: comic sans ms;">Bioinformatics,</span><br style="color: rgb(56, 83, 118); font-family: comic sans ms;"><span style="color: rgb(56, 83, 118); font-family: comic sans ms;">Department of Microbiology,<br>CCS Unversity,<br>Meerut</span><span style="color: rgb(56, 83, 118); font-family: comic sans
 ms;"><br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br></span></div><br><br>--- On <b>Sat, 27/3/10, Justin A. Lemkul <i>&lt;jalemkul@vt.edu&gt;</i></b> wrote:<br><blockquote style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; padding-left: 5px;"><br>From: Justin A. Lemkul &lt;jalemkul@vt.edu&gt;<br>Subject: Re: [gmx-users] POPC Membrane + protein dynamics<br>To: "Gromacs Users' List" &lt;gmx-users@gromacs.org&gt;<br>Date: Saturday, 27 March, 2010, 4:59 PM<br><br><div class="plainMail"><br><br>ravi sharma wrote:<br>&gt; <br>&gt; Dear,<br>&gt; This is the case when you added some ions in your pdb but not define it in your topology, Go in the end
 of you topology and subtract the number of molecule added. Also define the number of ions added and the name of the ion.<br>&gt; like :<br>&gt; NA&nbsp; &nbsp;&nbsp;&nbsp;14<br>&gt; <br><br>Addition of ions is one possibility, but is not definitively the answer here, especially since the difference in atoms is on the order of 10,000 (roughly 20% of the atoms in the system).&nbsp; In the case of the original post, it is far more likely that there is something more fundamentally wrong with the topology.<br><br>Addition of molecules to the topology requires the molecule name, not the atom name, be added to the [molecules] directive.&nbsp; I don't know any force field in Gromacs that uses "NA" as a molecule name; they are all either Na, Na+, or NA+.&nbsp; Any number of molecules subtracted must come from the group in which they are embedded, not just some generic "number of molecules" (stated for the sake of clarity).<br><br>-Justin<br><br>&gt; <br>&gt;
 <br>&gt; regards,<br>&gt; Ravi<br>&gt; <br>&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; <br>&gt; <br>&gt; --- On *Sat, 27/3/10, Justin A. Lemkul /&lt;<a ymailto="mailto:jalemkul@vt.edu" href="/mc/compose?to=jalemkul@vt.edu">jalemkul@vt.edu</a>&gt;/* wrote:<br>&gt; <br>&gt; <br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;From: Justin A. Lemkul &lt;<a ymailto="mailto:jalemkul@vt.edu" href="/mc/compose?to=jalemkul@vt.edu">jalemkul@vt.edu</a>&gt;<br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;Subject: Re: [gmx-users] POPC Membrane + protein dynamics<br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;To: "Gromacs Users' List" &lt;<a ymailto="mailto:gmx-users@gromacs.org" href="/mc/compose?to=gmx-users@gromacs.org">gmx-users@gromacs.org</a>&gt;<br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;Date: Saturday, 27 March, 2010, 1:50 AM<br>&gt; <br>&gt; <br>&gt;
 <br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;<a ymailto="mailto:padhu@ncbs.res.in" href="/mc/compose?to=padhu@ncbs.res.in">padhu@ncbs.res.in</a> &lt;/mc/compose?to=<a ymailto="mailto:padhu@ncbs.res.in" href="/mc/compose?to=padhu@ncbs.res.in">padhu@ncbs.res.in</a>&gt; wrote:<br>&gt;&nbsp; &nbsp; &nbsp; &gt; hi<br>&gt;&nbsp; &nbsp; &nbsp; &gt; i have tried as per your justin lemkul's tutorial.<br>&gt;&nbsp; &nbsp; &nbsp; &gt;<br>&gt;&nbsp; &nbsp; &nbsp; &gt; after aligning the protein with membrane, subsequently we tried<br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;to use<br>&gt;&nbsp; &nbsp; &nbsp; &gt; grompp to generate a tpr file i got the following error.<br>&gt;&nbsp; &nbsp; &nbsp; &gt;<br>&gt;&nbsp; &nbsp; &nbsp; &gt; -------------------------------------------------------<br>&gt;&nbsp; &nbsp; &nbsp; &gt; Program grompp_mpi, VERSION 3.3.1<br>&gt;&nbsp; &nbsp; &nbsp; &gt; Source code file: grompp.c, line: 448<br>&gt;&nbsp; &nbsp; &nbsp; &gt;<br>&gt;&nbsp; &nbsp; &nbsp; &gt;
 Fatal error:<br>&gt;&nbsp; &nbsp; &nbsp; &gt; number of coordinates in coordinate file (pope.pdb, 37867)<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; does not match topology (topol_pope.top, 46779)<br>&gt;&nbsp; &nbsp; &nbsp; &gt; -------------------------------------------------------<br>&gt;&nbsp; &nbsp; &nbsp; &gt;<br>&gt;&nbsp; &nbsp; &nbsp; &gt; can any one have idea about how to solve this error.<br>&gt;&nbsp; &nbsp; &nbsp; &gt;<br>&gt; <br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;Please see the general advice:<br>&gt; <br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;<a href="http://www.gromacs.org/Documentation/Errors#Number_of_coordinates_in_coordinate_file_does_not_match_topology" target="_blank">http://www.gromacs.org/Documentation/Errors#Number_of_coordinates_in_coordinate_file_does_not_match_topology</a><br>&gt; <br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;Beyond that, you'll have to post greater details of what step is<br>&gt;&nbsp;
 &nbsp;&nbsp;&nbsp;giving rise to this error, what you have done as far as inserting<br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;the protein, adding solvent and/or ions, etc.<br>&gt; <br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;I would also recommend using the latest version of Gromacs (4.0.7),<br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;unless for some reason you have to use an ancient version (3.3.1).&nbsp; &nbsp;&nbsp;&nbsp;The speed upgrade is incredible :)<br>&gt; <br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;-Justin<br>&gt; <br>&gt;&nbsp; &nbsp; &nbsp; &gt; thanks in advance,<br>&gt;&nbsp; &nbsp; &nbsp; &gt; Padhu.<br>&gt;&nbsp; &nbsp; &nbsp; &gt;<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt;<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt; <a ymailto="mailto:padhu@ncbs.res.in" href="/mc/compose?to=padhu@ncbs.res.in">padhu@ncbs.res.in</a> &lt;/mc/compose?to=<a ymailto="mailto:padhu@ncbs.res.in" href="/mc/compose?to=padhu@ncbs.res.in">padhu@ncbs.res.in</a>&gt; wrote:<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt;&gt; Dear
 All,<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt;&gt; I'm a new comer to gromacs. I need to perform molecular dynamics<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt;&gt; simulation of my protein within the POPC membrane. I have<br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;downloaded the<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt;&gt; 128a popc lipid from Prof.Tieleman's group along with the required<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt;&gt; popc.itp. My protein of interest is 458 residues. Since the<br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;128a popc is<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt;&gt; already Simulated for 1600 ps, i haven't done membrane only<br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;simulation.<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt;&gt; So<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt;&gt; i directly aligned my protein with the membrane using VMD and<br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;have saved<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt;&gt; the coordinates. However after i aligned my protein with the<br>&gt;&nbsp;
 &nbsp;&nbsp;&nbsp;membrane, i<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt;&gt; felt the membrane was not big enough for my protein of<br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;interest, also,<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt;&gt; the<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt;&gt; N-terminal region and the loop regions in the C-terminal region and<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt;&gt; stretching outside water molecule above and below<br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;respectively.ow, i<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt;&gt; have<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt;&gt; two questions here...<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt; Yep, if the protein is sticking "out" of the box on both sides,<br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;the box is<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt; too<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt; small.<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt;<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt;&gt; 1. How would i increase the number of lipids from the existing one?<br>&gt;&nbsp;
 &nbsp; &nbsp; &gt;&gt; Two options:<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt;<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt; 1. genconf -nbox will give integer multiples of the lipid system<br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;(i.e.,<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt; -nbox 2<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt; 2 1 will expand the bilayer in the x and y directions and give<br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;512 lipids)<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt;<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt; 2. genbox -cs popc128a.pdb -box (whatever dimensions) allows you<br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;a bit<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt; finer<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt; control over how large the box is and how many lipids will result.<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt;<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt;&gt; 2. Also, i should increase increase the number of water<br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;molecules above<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt;&gt; and below the
 membrane so that it covers the whole region of<br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;N-terminal<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt;&gt; and also the loop regions of the C-terminal.<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt; Expand the box using editconf and add solvent with genbox.<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt;<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt;&gt; Now, How to perform these tasks??If i do so will it be<br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;compatible with<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt;&gt; the<br>&gt;&nbsp; &nbsp; &nbsp; &gt;&gt;&gt; p<br>&gt;&nbsp; &nbsp; &nbsp; &gt;<br>&gt;&nbsp; &nbsp; &nbsp; &gt;<br>&gt; <br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;-- ========================================<br>&gt; <br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;Justin A. Lemkul<br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;Ph.D. Candidate<br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;ICTAS Doctoral Scholar<br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;MILES-IGERT Trainee<br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;Department of
 Biochemistry<br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;Virginia Tech<br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;Blacksburg, VA<br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;jalemkul[at]vt.edu | (540) 231-9080<br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;<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;&nbsp; &nbsp;&nbsp;&nbsp;========================================<br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;-- gmx-users mailing list&nbsp; &nbsp; <a ymailto="mailto:gmx-users@gromacs.org" href="/mc/compose?to=gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;&lt;/mc/compose?to=<a ymailto="mailto:gmx-users@gromacs.org" href="/mc/compose?to=gmx-users@gromacs.org">gmx-users@gromacs.org</a>&gt;<br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>&gt;&nbsp;
 &nbsp;&nbsp;&nbsp;Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before<br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;posting!<br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;Please don't post (un)subscribe requests to the list. Use the www<br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;interface or send it to <a ymailto="mailto:gmx-users-request@gromacs.org" href="/mc/compose?to=gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a><br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;&lt;/mc/compose?to=<a ymailto="mailto:gmx-users-request@gromacs.org" href="/mc/compose?to=gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>&gt;.<br>&gt;&nbsp; &nbsp;&nbsp;&nbsp;Can'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>&gt; <br>&gt; <br>&gt; ------------------------------------------------------------------------<br>&gt; The INTERNET now has
 a personality. YOURS! See your Yahoo! Homepage &lt;<a href="http://in.rd.yahoo.com/tagline_yyi_1/*http://in.yahoo.com/" target="_blank">http://in.rd.yahoo.com/tagline_yyi_1/*http://in.yahoo.com/</a>&gt;.<br><br>-- ========================================<br><br>Justin A. Lemkul<br>Ph.D. Candidate<br>ICTAS Doctoral Scholar<br>MILES-IGERT Trainee<br>Department of Biochemistry<br>Virginia Tech<br>Blacksburg, VA<br>jalemkul[at]vt.edu | (540) 231-9080<br><a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br><br>========================================<br>-- gmx-users mailing list&nbsp; &nbsp; <a ymailto="mailto:gmx-users@gromacs.org" href="/mc/compose?to=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't post (un)subscribe requests to the list. Use the www interface or send it to <a ymailto="mailto:gmx-users-request@gromacs.org" href="/mc/compose?to=gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>Can'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></blockquote></td></tr></table><br>



      <!--1--><hr size=1></hr> 
The INTERNET now has a personality. YOURS! <a href="http://in.rd.yahoo.com/tagline_yyi_1/*http://in.yahoo.com/" target="_blank">See your Yahoo! Homepage</a>.