<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;"><DIV>Dear sir:</DIV>
<DIV>Thank you for your help! <BR><BR><BR><A class=plink href="http://cn.webmessenger.yahoo.com/index.php?t=1&to=eWlkPXpoYW94aWl0YzIwMDI-&sig=703fa929658518b2720b087c59cd85f2dabf8844" target=_blank rel=nofollow><IMG class=pimg alt=4 src="http://opi.yahoo.com/online?u=zhaoxiitc2002&t=4&l=cn" border=0></A><BR><BR>--- <B>10年4月11日，周日, shayamra@post.tau.ac.il <I><shayamra@post.tau.ac.il></I></B> 写道：<BR></DIV>
<BLOCKQUOTE style="PADDING-LEFT: 5px; MARGIN-LEFT: 5px; BORDER-LEFT: rgb(16,16,255) 2px solid"><BR>发件人: shayamra@post.tau.ac.il <shayamra@post.tau.ac.il><BR>主题: Re: [gmx-users] InflateGRO and pentamer protein simulation<BR>收件人: jalemkul@vt.edu, "Discussion list for GROMACS users" <gmx-users@gromacs.org><BR>日期: 2010年4月11日,周日,下午10:01<BR><BR>
<DIV class=plainMail>Dear Xi Zhao,<BR><BR>For inserting protein into membrane you might like to try use<BR>g_membed<BR>(link <A href="http://wwwuser.gwdg.de/~ggroenh/membed.html" target=_blank>http://wwwuser.gwdg.de/~ggroenh/membed.html</A>)<BR>This is not strictly (yet) a gromacs package, but it can be very easily be complied and used per instructions of the authors.<BR><BR>Technically, the package 'compress' the protein to be a linear line, then inserts it into the membrane and starts 'pumping' it back to the original size.<BR>Hope that help,<BR><BR>-Shay<BR><BR>Quoting "Justin A. Lemkul" <<A href="http://cn.mc151.mail.yahoo.com/mc/compose?to=jalemkul@vt.edu" ymailto="mailto:jalemkul@vt.edu">jalemkul@vt.edu</A>>:<BR><BR>> <BR>> <BR>> xi zhao wrote:<BR>>> I would like to run a simulation of a pentamers in a POPC membrane, and using new inflategro with doughnut mode, but I have met similar result as the links <A
href="http://lists.gromacs.org/pipermail/gmx-users/2010-January/047977.html" target=_blank>http://lists.gromacs.org/pipermail/gmx-users/2010-January/047977.html</A> ,I also <<A href="http://lists.gromacs.org/pipermail/gmx-users/2010-January/047977.html" target=_blank>http://lists.gromacs.org/pipermail/gmx-users/2010-January/047977.html</A> ,I also> got a lot of errors about uninitialized values. When the script finishes, the membrane is rescaled, but the protein is untouched in the corner of the new box. Please give me some suggestions! Thank you in advance!<BR>>> > best regards!<BR>> <BR>> I think the problem comes from pattern matching in the script. With large bilayers, the atom names and numbers will not be split correctly; this could be causing a problem for you. I'd again suggest that you contact the developer of the program and discuss this issue, as it is not really a Gromacs problem and only a few
people on this list really use such programs. I think you'll find a resolution faster that way.<BR>> <BR>> -Justin<BR>> <BR>>> <BR>>> <BR>>> <BR>>> 4 <<A href="http://cn.webmessenger.yahoo.com/index.php?t=1&to=eWlkPXpoYW94aWl0YzIwMDI-&sig=703fa929658518b2720b087c59cd85f2dabf8844" target=_blank>http://cn.webmessenger.yahoo.com/index.php?t=1&to=eWlkPXpoYW94aWl0YzIwMDI-&sig=703fa929658518b2720b087c59cd85f2dabf8844</A>><BR>>> <BR>>> --- *10年4月9日，周五, Justin A. Lemkul /<<A href="http://cn.mc151.mail.yahoo.com/mc/compose?to=jalemkul@vt.edu" ymailto="mailto:jalemkul@vt.edu">jalemkul@vt.edu</A>>/* 写道：<BR>>> <BR>>> <BR>>> 发件人: Justin A. Lemkul <<A href="http://cn.mc151.mail.yahoo.com/mc/compose?to=jalemkul@vt.edu" ymailto="mailto:jalemkul@vt.edu">jalemkul@vt.edu</A>><BR>>> 主题: Re: [gmx-users] InflateGRO
and pentamer protein simulation<BR>>> 收件人: "Discussion list for GROMACS users" <<A href="http://cn.mc151.mail.yahoo.com/mc/compose?to=gmx-users@gromacs.org" ymailto="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</A>><BR>>> 日期: 2010年4月9日,周五,上午2:14<BR>>> <BR>>> <BR>>> <BR>>> xi zhao wrote:<BR>>> ><BR>>> ><BR>>> > Dear sir :<BR>>> > I want to a pentamer membrane protein, when I used InflateGRO.pl<BR>>> with /DOUGHNUT Mode, the results were wrong, please help me! /<BR>>> <BR>>> If you want any useful help, you'll have to do a whole lot better<BR>>> than simply saying "the results were wrong." No one on this list<BR>>> will have any
idea what you mean. If you believe there is some<BR>>> error in the script itself, you're better off contacting its author.<BR>>> <BR>>> -Justin<BR>>> <BR>>> > 4<BR>>> <<A href="http://cn.webmessenger.yahoo.com/index.php?t=1&to=eWlkPXpoYW94aWl0YzIwMDI-&sig=703fa929658518b2720b087c59cd85f2dabf8844" target=_blank>http://cn.webmessenger.yahoo.com/index.php?t=1&to=eWlkPXpoYW94aWl0YzIwMDI-&sig=703fa929658518b2720b087c59cd85f2dabf8844</A><BR>>> <<A href="http://cn.webmessenger.yahoo.com/index.php?t=1&to=eWlkPXpoYW94aWl0YzIwMDI-&sig=703fa929658518b2720b087c59cd85f2dabf8844" target=_blank>http://cn.webmessenger.yahoo.com/index.php?t=1&to=eWlkPXpoYW94aWl0YzIwMDI-&sig=703fa929658518b2720b087c59cd85f2dabf8844</A>>>show<BR>>> " [user@localhost protein-tutorial]$ perl
inflategro.txt kkr.gro 4<BR>>> POPC 14 kk_inflated.gro 5 area.dat doughnut protein_subunits<BR>>> > Doughnut mode activated. Protein coordinates will be translated<BR>>> by subunit ....<BR>>> > Reading.....<BR>>> > Reading chain identifiers....<BR>>> > Subunit 1: atom 1 to atom 3238<BR>>> > Subunit 2: atom 3239 to atom 6476<BR>>> > Subunit 3: atom 6477 to atom 9714<BR>>> > Subunit 4: atom 9715 to atom 12952<BR>>> > Subunit 5: atom 12953 to atom 16190<BR>>> > There are 5 protein subunits....<BR>>> > Scaling
lipids....<BR>>> > There are 512 lipids...<BR>>> > with 65 atoms per lipid..<BR>>> > Determining upper and lower leaflet...<BR>>> > 256 lipids in the upper...<BR>>> > 256 lipids in the lower leaflet<BR>>> > Checking for overlap....<BR>>> > ...this might actually take a while...<BR>>> > .......<BR>>> > Argument "N4" isn't numeric in printf at inflategro.txt line 708,<BR>>> <CHAINS> line 5.<BR>>> > Argument "C5" isn't numeric in printf at inflategro.txt line 708,<BR>>> <CHAINS> line 5.<BR>>> > Argument "C6" isn't numeric in printf at inflategro.txt
line 708,<BR>>> <CHAINS> line 5.<BR>>> > Argument "O7" isn't numeric in printf at inflategro.txt line 708,<BR>>> <CHAINS> line 5.<BR>>> > Argument "P8" isn't numeric in printf at inflategro.txt line 708,<BR>>> <CHAINS> line 5.<BR>>> > Argument "O9" isn't numeric in printf at inflategro.txt line 708,<BR>>> <CHAINS> line 5.<BR>>> > Argument "A1" isn't numeric in printf at inflategro.txt line 708,<BR>>> <CHAINS> line 5.<BR>>> > Argument "A2" isn't numeric in printf at inflategro.txt line 708,<BR>>> <CHAINS> line 5.<BR>>> > Argument "C1" isn't numeric in printf at inflategro.txt line 708,<BR>>>
<CHAINS> line 5.<BR>>> > Argument "C2" isn't numeric in printf at inflategro.txt line 708,<BR>>> <CHAINS> line 5.<BR>>> > Argument "C3" isn't numeric in printf at inflategro.txt line 708,<BR>>> <CHAINS> line 5.<BR>>> > Argument "N4" isn't numeric in printf at inflategro.txt line 708,<BR>>> <CHAINS> line 5.<BR>>> > Argument "C5" isn't numeric in printf at inflategro.txt line 708,<BR>>> <CHAINS> line 5.<BR>>> > Argument "C6" isn't numeric in printf at inflategro.txt line 708,<BR>>> <CHAINS> line 5.<BR>>> > Argument "O7" isn't numeric in printf at inflategro.txt line 708,<BR>>> <CHAINS> line
5.<BR>>> > Argument "P8" isn't numeric in printf at inflategro.txt line 708,<BR>>> <CHAINS> line 5.<BR>>> > Argument "O9" isn't numeric in printf at inflategro.txt line 708,<BR>>> <CHAINS> line 5.<BR>>> > Argument "A1" isn't numeric in printf at inflategro.txt line 708,<BR>>> <CHAINS> line 5.<BR>>> > Argument "A2" isn't numeric in printf at inflategro.txt line 708,<BR>>> <CHAINS> line 5.<BR>>> > Calculating Area per lipid...<BR>>> > Protein X-min/max: 23 105<BR>>> > Protein Y-min/max: 18 100<BR>>> > X-range: 82 A Y-range: 82 A<BR>>>
> Building 82 X 82 2D grid on protein coordinates...<BR>>> > Calculating area occupied by protein..<BR>>> > full TMD..<BR>>> > upper TMD....<BR>>> > lower TMD....<BR>>> > Area per protein: 49 nm^2<BR>>> > Area per lipid: 9.83158146861789 nm^2<BR>>> > Area per protein, upper half: 38.75 nm^2<BR>>> > Area per lipid, upper leaflet : 9.79362516645161 nm^2<BR>>> > Area per protein, lower half: 45.75 nm^2<BR>>> > Area per lipid, lower leaflet : 9.92548787409836 nm^2<BR>>> > Writing Area per lipid...<BR>>> > Done!<BR>>>
><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 <A href="http://cn.mc151.mail.yahoo.com/mc/compose?to=gmx-users@gromacs.org"
ymailto="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</A><BR>>> <<A href="http://cn.mc151.mail.yahoo.com/mc/compose?to=gmx-users@gromacs.org" target=_blank>http://cn.mc151.mail.yahoo.com/mc/compose?to=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<BR>>> posting!<BR>>> Please don't post (un)subscribe requests to the list. Use the www<BR>>> interface or send it to <A href="http://cn.mc151.mail.yahoo.com/mc/compose?to=gmx-users-request@gromacs.org" ymailto="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</A><BR>>> <<A
href="http://cn.mc151.mail.yahoo.com/mc/compose?to=gmx-users-request@gromacs.org" target=_blank>http://cn.mc151.mail.yahoo.com/mc/compose?to=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>>> <BR>>> <BR>>> <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 <A
href="http://cn.mc151.mail.yahoo.com/mc/compose?to=gmx-users@gromacs.org" ymailto="mailto: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 href="http://cn.mc151.mail.yahoo.com/mc/compose?to=gmx-users-request@gromacs.org" ymailto="mailto: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>> <BR><BR><BR>--<BR>gmx-users mailing list <A
href="http://cn.mc151.mail.yahoo.com/mc/compose?to=gmx-users@gromacs.org" ymailto="mailto: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<BR>www interface or send it to <A href="http://cn.mc151.mail.yahoo.com/mc/compose?to=gmx-users-request@gromacs.org" ymailto="mailto: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>