<P>
&nbsp; <BR>
Thanks for the response<BR>
Just diverting this topic to about specific number of popc molecules.<BR>
<BR>
I created the bilayer by using genconf command <BR>
genconf -f popc128a.pdb -o out.gro -dist 0 0 0 -nbox 2 1 1 (as I posted&nbsp;  in my previous mail) generated output file contain 128 popc in each leaflet of bilayer. <BR>
<BR>
If you see original popc box dimensions 6.1x6.2x6.9 (means in all dimensions popc number almost same)but with genconf command above mentioned options created box values 12x6.1x6.9. I dont want that many popc molecules because in X-dimension too many popc molecules are present.&nbsp;  <BR>
<BR>
1.is there anyway to reduce those popc molecules from 128 to 80/90 popc molecules? or <BR>
2.I wanted to create popc molecules 80 or 90 in eachleaflet is it possible to generate?<BR>
<BR>
These are may be trivial queries<BR>
Could you give suggest me please<BR>
Thanks in advance.<BR>
<BR>
<BR>
&gt;&gt;&gt;Thanks Justin and Nicolas for gave suggestions.<BR>
&gt;&gt;&gt;I have tried with Nicolas suggested command in VMD Tkconsole, its showing whole popc molecules but not leaflet. I typed the command in Tkconsole like this<BR>
&gt;&gt;&gt;<BR>
&gt;&gt;&gt;[atomselect top &quot;name P8 and z&gt;0&quot;] num<BR>
&gt;&gt;&gt;it has showed 201, means the total number popc molecues in the .gro file.<BR>
&gt;&gt;&gt;Could you tell any suggestion<BR>
&gt;&gt;&gt;<BR>
&gt;&gt;It's probably because your system is not center on 0.0. By default, Gromacs center boxes on lenght/2<BR>
&gt;It's probably because your system is not centered on 0.0. By default, Gromacs centers boxes on lenght/2<BR>
&gt;.. Sorry for the grammar, I'm not well awake this morning.<BR>
&gt;&gt;&gt;<BR>
&gt;&gt;&gt;<BR>
&gt;&gt;&gt;Thanks in advance.<BR>
&gt;&gt;&gt;e<BR>
&gt;&gt;&gt; &gt; &gt;&gt;<BR>
&gt;&gt;&gt; &gt; &gt;&gt; Hi Jochen thanks for your reply<BR>
&gt;&gt;&gt; &gt; &gt;&gt; I have gone through this recent mail<BR>
&gt;&gt;&gt; &gt; &gt;&gt; http://www.gromacs.org/pipermail/gmx-users/2008-September/036508.html<BR>
&gt;&gt;&gt; &gt; &gt;&gt; more over if I use genconf command like this<BR>
&gt;&gt;&gt; &gt; &gt;&gt; genconf -f .gro -o out -nbox 2 1 1 -dist 0 0 0 its adding 128 in<BR>
&gt;&gt;&gt; &gt; &gt;&gt; eachleaflet I dont wany that many popc molecules.<BR>
&gt;&gt;&gt; &gt; &gt;&gt;<BR>
&gt;&gt;&gt; &gt; &gt;&gt; 1.Is it wrong if I increase the popc molecules by using genbox?<BR>
&gt;&gt;&gt; &gt; &gt;<BR>
&gt;&gt;&gt; &gt; &gt; It is best to use genconf, because then the periodic images of the<BR>
&gt;&gt;&gt; &gt; &gt; unit cell remain intact, that is, since you're using a<BR>
&gt;&gt;&gt; &gt; &gt; pre-equilibrated bilayer, it's better to not snip chunks out of it.<BR>
&gt;&gt;&gt; &gt; &gt; You can deal with that by sufficient equilibration, however.<BR>
&gt;&gt;&gt; &gt; &gt;<BR>
&gt;&gt;&gt; &gt; &gt; It is also easier to use genconf, because you then know exactly how<BR>
&gt;&gt;&gt; &gt; &gt; many lipids you are dealing with (in regards to your previous<BR>
&gt;&gt;&gt; &gt; &gt; message).&nbsp; You could probably write some script to tell you which<BR>
&gt;&gt;&gt; &gt; &gt; lipid is in a given leaflet based on whether a certain atom (i.e., P8<BR>
&gt;&gt;&gt; &gt; &gt; or something else) is above or below the center of the bilayer.<BR>
&gt;&gt;&gt; &gt;In case you use VMD, you can get the number of phospholipid per leaflet<BR>
&gt;&gt;&gt; &gt;with the following command:<BR>
&gt;&gt;&gt; &gt;<BR>
&gt;&gt;&gt; &gt;&nbsp; &nbsp; [atomselect top &quot;name P8 and z&gt;0&quot;] num<BR>
&gt;&gt;&gt; &gt;<BR>
&gt;&gt;&gt; &gt;This will give you the number of PC in the upper leaflet, assuming 1)<BR>
&gt;&gt;&gt; &gt;the phosphorus atom is named P8 and 2) the bilayer is center on 0.0<BR>
&gt;&gt;&gt; &gt;along the z axis.<BR>
&gt;&gt;&gt; &gt;<BR>
&gt;&gt;&gt; &gt;Nicolas<BR>
&gt;&gt;&gt;<BR>
&gt;&gt;&gt;<BR>
&gt;&gt;&gt; &gt;<BR>
&gt;&gt;&gt; &gt; &gt;<BR>
&gt;&gt;&gt; &gt; &gt;&gt; 2.Is there anyway to increase popc and water numbers by mentioning<BR>
&gt;&gt;&gt; &gt; &gt;&gt; specific molecules number?<BR>
&gt;&gt;&gt; &gt; &gt;<BR>
&gt;&gt;&gt; &gt; &gt; Not that I'm aware of.&nbsp; There is a -maxsol option in genbox, but that<BR>
&gt;&gt;&gt; &gt; &gt; is for capping the amount of water molecules added to a box.<BR>
&gt;&gt;&gt; &gt; &gt;<BR>
&gt;&gt;&gt; &gt; &gt; -Justin<BR>
&gt;&gt;&gt; &gt; &gt;<BR>
&gt;&gt;&gt; &gt; &gt;&gt; Could you suggest me<BR>
&gt;&gt;&gt; &gt; &gt;&gt; Thanks in advance.<BR>
&gt;&gt;&gt; &gt; &gt;&gt;<BR>
&gt;&gt;&gt; &gt; &gt;&gt; On Fri, 19 Sep 2008 Jochen Hub wrote :<BR>
&gt;&gt;&gt; &gt; &gt;&gt;&nbsp; &gt;minnale wrote:<BR>
&gt;&gt;&gt; &gt; &gt;&gt;&nbsp; &gt; &gt;<BR>
&gt;&gt;&gt; &gt; &gt;&gt;&nbsp; &gt; &gt; Hi all,<BR>
&gt;&gt;&gt; &gt; &gt;&gt;&nbsp; &gt; &gt;&nbsp; &nbsp; I have extended popc bilayer(intial popc.pdb from Dr.Tielmen<BR>
&gt;&gt;&gt; &gt; &gt;&gt; site) by using genbox command, I issued<BR>
&gt;&gt;&gt; &gt; &gt;&gt;&nbsp; &gt; &gt; genbox -cs popc128a.gro -o out.gro -box 9.2 9.2 6.9 it ran<BR>
&gt;&gt;&gt; &gt; &gt;&gt; successfully with increase of popc and water molecules.<BR>
&gt;&gt;&gt; &gt; &gt;&gt;&nbsp; &gt; &gt;&nbsp; Now I want to visualise this out file in VMD in a way that in<BR>
&gt;&gt;&gt; &gt; &gt;&gt; eachleaflet how many popc molecules and water residues are there, May<BR>
&gt;&gt;&gt; &gt; &gt;&gt; be this is trivial query.<BR>
&gt;&gt;&gt; &gt; &gt;&gt;&nbsp; &gt; &gt; Could you give me suggestion.<BR>
&gt;&gt;&gt; &gt; &gt;&gt;&nbsp; &gt;<BR>
&gt;&gt;&gt; &gt; &gt;&gt;&nbsp; &gt;If you want to enlarge a membrane patch, use genconf. Not genbox!<BR>
&gt;&gt;&gt; &gt; &gt;&gt;&nbsp; &gt;<BR>
&gt;&gt;&gt; &gt; &gt;&gt;&nbsp; &gt;jochen<BR>
&gt;&gt;&gt;<BR>
&gt;&gt;&gt;<BR>
&gt;&gt;&gt;<BR>
&gt;&gt;&gt;Ebay &lt;http://adworks.rediff.com/cgi-bin/AdWorks/click.cgi/www.rediff.com/signature-home.htm/1050715198@Middle5/2401775_2394076/2397136/1?PARTNER=3&amp;OAS_QUERY=null&gt; <BR>
&gt;&gt;&gt;<BR>
&gt;&gt;&gt;------------------------------------------------------------------------<BR>
&gt;&gt;&gt;<BR>
&gt;&gt;&gt;_______________________________________________<BR>
&gt;&gt;&gt;gmx-users mailing list&nbsp; &nbsp; gmx-users@gromacs.org<BR>
&gt;&gt;&gt;http://www.gromacs.org/mailman/listinfo/gmx-users<BR>
&gt;&gt;&gt;Please search the archive at http://www.gromacs.org/search before posting!<BR>
&gt;&gt;&gt;Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-request@gromacs.org.<BR>
&gt;&gt;&gt;Can't post? Read http://www.gromacs.org/mailing_lists/users.php<BR>
&gt;&gt;<BR>
&gt;&gt;_______________________________________________<BR>
&gt;&gt;gmx-users mailing list&nbsp; &nbsp; gmx-users@gromacs.org<BR>
&gt;&gt;http://www.gromacs.org/mailman/listinfo/gmx-users<BR>
&gt;&gt;Please search the archive at http://www.gromacs.org/search before posting!<BR>
&gt;&gt;Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-request@gromacs.org.<BR>
&gt;&gt;Can't post? Read http://www.gromacs.org/mailing_lists/users.php<BR>
&gt;<BR>

</P>
<br><br>
<Table border=0 Width=644 Height=57 cellspacing=0 cellpadding=0 style='font-family:Verdana;font-size:11px;line-height:15px;'><TR><td><a href='http://adworks.rediff.com/cgi-bin/AdWorks/click.cgi/www.rediff.com/signature-home.htm/1050715198@Middle5/2401775_2394076/2397136/1?PARTNER=3&OAS_QUERY=null' target=new ><img src ='http://imadworks.rediff.com/cgi-bin/AdWorks/adimage.cgi/2401775_2394076/creative_2397136.gif'  alt='Ebay'  border=0></a></td></TR></Table>