<P>
&nbsp; <BR>
<BR>
<BR>
Thanks Justin and Nicolas for gave suggestions.<BR>
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>
<BR>
[atomselect top &quot;name P8 and z&gt;0&quot;] num <BR>
it has showed 201, means the total number popc molecues in the .gro file.<BR>
Could you tell any suggestion<BR>
<BR>
Thanks in advance. <BR>
<BR>
&gt; &gt;&gt;<BR>
&gt; &gt;&gt; Hi Jochen thanks for your reply<BR>
&gt; &gt;&gt; I have gone through this recent mail<BR>
&gt; &gt;&gt; http://www.gromacs.org/pipermail/gmx-users/2008-September/036508.html<BR>
&gt; &gt;&gt; more over if I use genconf command like this<BR>
&gt; &gt;&gt; genconf -f .gro -o out -nbox 2 1 1 -dist 0 0 0 its adding 128 in<BR>
&gt; &gt;&gt; eachleaflet I dont wany that many popc molecules.<BR>
&gt; &gt;&gt;<BR>
&gt; &gt;&gt; 1.Is it wrong if I increase the popc molecules by using genbox?<BR>
&gt; &gt;<BR>
&gt; &gt; It is best to use genconf, because then the periodic images of the<BR>
&gt; &gt; unit cell remain intact, that is, since you're using a<BR>
&gt; &gt; pre-equilibrated bilayer, it's better to not snip chunks out of it.<BR>
&gt; &gt; You can deal with that by sufficient equilibration, however.<BR>
&gt; &gt;<BR>
&gt; &gt; It is also easier to use genconf, because you then know exactly how<BR>
&gt; &gt; many lipids you are dealing with (in regards to your previous<BR>
&gt; &gt; message).&nbsp; You could probably write some script to tell you which<BR>
&gt; &gt; lipid is in a given leaflet based on whether a certain atom (i.e., P8<BR>
&gt; &gt; or something else) is above or below the center of the bilayer.<BR>
&gt;In case you use VMD, you can get the number of phospholipid per leaflet<BR>
&gt;with the following command:<BR>
&gt;<BR>
&gt;&nbsp; &nbsp;  [atomselect top &quot;name P8 and z&gt;0&quot;] num<BR>
&gt;<BR>
&gt;This will give you the number of PC in the upper leaflet, assuming 1)<BR>
&gt;the phosphorus atom is named P8 and 2) the bilayer is center on 0.0<BR>
&gt;along the z axis.<BR>
&gt;<BR>
&gt;Nicolas<BR>
<BR>
<BR>
&gt;<BR>
&gt; &gt;<BR>
&gt; &gt;&gt; 2.Is there anyway to increase popc and water numbers by mentioning<BR>
&gt; &gt;&gt; specific molecules number?<BR>
&gt; &gt;<BR>
&gt; &gt; Not that I'm aware of.&nbsp; There is a -maxsol option in genbox, but that<BR>
&gt; &gt; is for capping the amount of water molecules added to a box.<BR>
&gt; &gt;<BR>
&gt; &gt; -Justin<BR>
&gt; &gt;<BR>
&gt; &gt;&gt; Could you suggest me<BR>
&gt; &gt;&gt; Thanks in advance.<BR>
&gt; &gt;&gt;<BR>
&gt; &gt;&gt; On Fri, 19 Sep 2008 Jochen Hub wrote :<BR>
&gt; &gt;&gt;&nbsp; &gt;minnale wrote:<BR>
&gt; &gt;&gt;&nbsp; &gt; &gt;<BR>
&gt; &gt;&gt;&nbsp; &gt; &gt; Hi all,<BR>
&gt; &gt;&gt;&nbsp; &gt; &gt;&nbsp; &nbsp; I have extended popc bilayer(intial popc.pdb from Dr.Tielmen<BR>
&gt; &gt;&gt; site) by using genbox command, I issued<BR>
&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; successfully with increase of popc and water molecules.<BR>
&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; eachleaflet how many popc molecules and water residues are there, May<BR>
&gt; &gt;&gt; be this is trivial query.<BR>
&gt; &gt;&gt;&nbsp; &gt; &gt; Could you give me suggestion.<BR>
&gt; &gt;&gt;&nbsp; &gt;<BR>
&gt; &gt;&gt;&nbsp; &gt;If you want to enlarge a membrane patch, use genconf. Not genbox!<BR>
&gt; &gt;&gt;&nbsp; &gt;<BR>
&gt; &gt;&gt;&nbsp; &gt;jochen<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>