<html><head><style type="text/css"><!-- DIV {margin:0px;} --></style></head><body><div style="font-family:times new roman, new york, times, serif;font-size:12pt"><DIV style="FONT-SIZE: 12pt; FONT-FAMILY: times new roman, new york, times, serif">
<DIV style="FONT-SIZE: 12pt; FONT-FAMILY: times new roman, new york, times, serif">---- Original Message ----<BR>From: Mohamed Osman <osmangmx@gmail.com><BR>To: gmx-users@gromacs.org<BR>Sent: Thursday, January 11, 2007 10:49:29 PM<BR>Subject: [gmx-users] box size for simulation of membrane protein<BR><BR>My questions: <BR>1. Protein is closer to one edge than the other. Does this make any difference? <BR> The solvants are part of equilibritaed lipid.<BR></DIV>
<DIV style="FONT-SIZE: 12pt; FONT-FAMILY: times new roman, new york, times, serif">-Shouldn't do. It's still the same distance from its periodic image on each side, which is what matters. As long as it's in the correct place in the bilayer, I don't believe there are any 'edge effects' in the implementation of PBC in gromacs you need to worry about.</DIV>
<DIV style="FONT-SIZE: 12pt; FONT-FAMILY: times new roman, new york, times, serif"><BR>Along protein axis:<BR> a. The periodic image is 1.43 nm away. (along protein axis) <BR> b. The distance between the protein and its mirror image is: 0.8 nm on one side and 2.0 nm on the other side? (These affect PME - coulmb interactions). Will these cause trouble?</DIV>
<DIV style="FONT-SIZE: 12pt; FONT-FAMILY: times new roman, new york, times, serif"> </DIV>
<DIV style="FONT-SIZE: 12pt; FONT-FAMILY: times new roman, new york, times, serif">-The minimum protein-image distance is all that matters. It doesn't matter how close the protein is to one side of the box - if it's closer on one side, it's further on the other. Unless you mean different distances on the x and y axis?</DIV>
<DIV style="FONT-SIZE: 12pt; FONT-FAMILY: times new roman, new york, times, serif">-As far as I know, PME sums out to infinity, so you can't really change that. BUT:<BR><BR>2. Is it necessary to keep the protein at least 1.5 nm away from the edges normal to z-axis? This more than doubles the number of solvant molecules. ---- <BR><BR>I have simulations with less than 50,000 atoms (protein, solvants, and KcsA). But no one realy says much about the box sizes. <BR><BR>What is the common practice? Please give examples?<BR>and <BR>How much compromise is acceptable ? </DIV>
<DIV style="FONT-SIZE: 12pt; FONT-FAMILY: times new roman, new york, times, serif">-The official requirement is for the minimim distance from protein to its periodic image to be double the cutoff, so that no atom feels forces from more than one image of the protein. However, you can probably get away with reducing that minimum distance to slightly more than the cutoff (1.5nm in your case) which just stops the protein from interacting with its own periodic images. Most people work on this basis, on the grounds it probably doesn't matter if the water molecules are 'feeling' more than one protein. It's a judgement call really, depending on how badly you want it to run faster vs how much you lose by making it do that.</DIV>
<DIV style="FONT-SIZE: 12pt; FONT-FAMILY: times new roman, new york, times, serif"> </DIV>
<DIV style="FONT-SIZE: 12pt; FONT-FAMILY: times new roman, new york, times, serif"> </DIV>
<DIV style="FONT-SIZE: 12pt; FONT-FAMILY: times new roman, new york, times, serif">"<SPAN style="FONT-STYLE: italic">It is common to compormize in this respect !!!"</SPAN></DIV><SPAN style="FONT-STYLE: italic"></SPAN></DIV>
<P><SPAN style="FONT-STYLE: italic"></SPAN> </P>
<P><SPAN style="FONT-STYLE: italic"></SPAN> </P>
<P><SPAN style="FONT-STYLE: italic"></SPAN> </P>
<P><SPAN style="FONT-STYLE: italic"> </P>
<DIV style="FONT-SIZE: 12pt; FONT-FAMILY: times new roman, new york, times, serif"><BR><BR>Thanks<BR>Mohamed Osman <BR><SPAN style="FONT-STYLE: italic"><BR>=======================<BR></SPAN>Professor of Electrical Engineering<BR>Washington State University<BR>Pullman, WA 99164-2752<BR></SPAN><BR></DIV>
<DIV style="FONT-SIZE: 12pt; FONT-FAMILY: times new roman, new york, times, serif">_______________________________________________<BR>gmx-users mailing list gmx-users@gromacs.org<BR><A href="http://www.gromacs.org/mailman/listinfo/gmx-users" target=_blank>http://www.gromacs.org/mailman/listinfo/gmx-users</A><BR>Please don't post (un)subscribe requests to the list. Use the <BR>www interface or send it to gmx-users-request@gromacs.org.<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></DIV>
<DIV style="FONT-SIZE: 12pt; FONT-FAMILY: times new roman, new york, times, serif"><BR></DIV></div><br>
<hr size=1>Looking for earth-friendly autos? <br> <a href="http://autos.yahoo.com/green_center/;_ylc=X3oDMTE4MGw4Z2hlBF9TAzk3MTA3MDc2BHNlYwNtYWlsdGFncwRzbGsDZ3JlZW5jZW50ZXI-">Browse Top Cars by "Green Rating"</a> at Yahoo! Autos' Green Center. </body></html>