<br><br>----- Original Message -----<br>From: Ryogo Sugitani &lt;rsugitani@ucdavis.edu&gt;<br>Date: Wednesday, September 3, 2008 5:46 pm<br>Subject: Re: [gmx-users] Re: Re: Bonds break while Minimising using distance restraints<br>To: plmallip@mail.uh.edu<br><br>&gt; Latha,<br>&gt; <br>&gt; I'm not sure why there is 1-4 interaction between 6151 and 6160<br>&gt; by just making distance restraint between 1703 and 1712<br>&gt; (because they shouldn't make a physical bond...)<br>&gt; <br>&gt; but maybe adding space before and after distace_restraints <br>&gt; (between brackets) might help?<br>&gt; Also, check the trajectory on VMD. <br>&gt; The warning is given, but md (or minimization) itself might be <br>&gt; running fine...<br>&gt; <br>&gt; Other than this, I probably won't be able to help you any further.<br>&gt; <br>&gt; Also, please post it back to gmx-ml for me...<br>&gt; (just reply to this message and change the address to gmx's list)<br>&gt; <br>&gt; Thanks,<br>&gt; <br>&gt; Ryogo<br>&gt; <br>&gt; <br>&gt; <br>&gt; On Wed, 03 Sep 2008 16:43:45 -0500, plmallip@mail.uh.edu wrote:<br>&gt; &gt; Hi Ryogo,<br>&gt; &gt; <br>&gt; &gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Thanks <br>&gt; for your suggestion. You are right. But, I did do the <br>&gt; &gt; energy minimization assuming the units are in nm But, I am <br>&gt; still <br>&gt; &gt; facing the same problem of the residues to which I imposed <br>&gt; distance <br>&gt; &gt; restraint breaking into fragments. I get this warning -&nbsp; <br>&gt; 1-4 <br>&gt; &gt; interaction between 6151 and 6160 at distance 1.387 which is <br>&gt; larger <br>&gt; &gt; than the 1-4 table size 1.000 &gt; nm. These are ignored for the <br>&gt; rest of <br>&gt; &gt; the simulation. This <br>&gt; &gt; usually means your system is exploding. <br>&gt; &gt; <br>&gt; &gt;&nbsp;&nbsp;&nbsp;&nbsp; Do you have any other suggestions?<br>&gt; &gt; <br>&gt; &gt; Thanks &amp; regards,<br>&gt; &gt; Latha.<br>&gt; &gt; <br>&gt; &gt; ----- Original Message -----<br>&gt; &gt; From: Ryogo Sugitani &lt;rsugitani@ucdavis.edu&gt;<br>&gt; &gt; Date: Wednesday, September 3, 2008 2:58 pm<br>&gt; &gt; Subject: Re: [gmx-users] Re: Re: Bonds break while Minimising <br>&gt; using <br>&gt; &gt; distance restraints<br>&gt; &gt; To: plmallip@mail.uh.edu<br>&gt; &gt; <br>&gt; &gt;&gt; Hi Latha,<br>&gt; &gt;&gt; <br>&gt; &gt;&gt; Since I cannot post it to the gmx ML directly for some <br>&gt; reason, <br>&gt; &gt;&gt; I'll just give this&nbsp; e-mail to you.<br>&gt; &gt;&gt; <br>&gt; &gt;&gt; I believe the unit for the distance used in <br>&gt; distance_restraints <br>&gt; &gt;&gt; is nm.<br>&gt; &gt;&gt; So, you should change it to something like 1.15 1.20 1.20 for <br>&gt; &gt;&gt; your setting.<br>&gt; &gt;&gt; Currently you have 115A, 120A, 120A...<br>&gt; &gt;&gt; <br>&gt; &gt;&gt; Hope that helps,<br>&gt; &gt;&gt; <br>&gt; &gt;&gt; Best,<br>&gt; &gt;&gt; <br>&gt; &gt;&gt; Ryogo<br>&gt; &gt;&gt; <br>&gt; &gt;&gt; <br>&gt; &gt;&gt; On Wed, 03 Sep 2008 13:27:06 -0500, plmallip@mail.uh.edu wrote:<br>&gt; &gt;&gt;&gt; Hi,<br>&gt; &gt;&gt;&gt; <br>&gt; &gt;&gt;&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; There is no chance of <br>&gt; steric <br>&gt; &gt;&gt; clashes between 6151 &amp; 6160. They <br>&gt; &gt;&gt;&gt; are seperated by 12 A. ("Warning: 1-4 interaction between <br>&gt; 6151 <br>&gt; &gt;&gt; and <br>&gt; &gt;&gt;&gt; 6160 at distance 1.387 which is larger than the 1-4 table <br>&gt; size <br>&gt; &gt;&gt; 1.000 <br>&gt; &gt;&gt;&gt; nm. These are ignored for the rest of the simulation. This <br>&gt; &gt;&gt; usually <br>&gt; &gt;&gt;&gt; means your system is exploding") I cannot build the missing <br>&gt; &gt;&gt; residues <br>&gt; &gt;&gt;&gt; without knowing the secondary structure. I am already <br>&gt; running <br>&gt; &gt;&gt; one <br>&gt; &gt;&gt;&gt; simulation with built residues. But, my ultimate goal is to <br>&gt; &gt;&gt; run <br>&gt; &gt;&gt;&gt; dynamics with distance restraints.<br>&gt; &gt;&gt;&gt;&nbsp;&nbsp; <br>&gt; &gt;&gt;&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; My sytem is <br>&gt; an <br>&gt; &gt;&gt; pentamer. Here is the toplogy file<br>&gt; &gt;&gt;&gt; <br>&gt; &gt;&gt;&gt; ; Include forcefield parameters<br>&gt; &gt;&gt;&gt; #include "ffG43a1.itp"<br>&gt; &gt;&gt;&gt; <br>&gt; &gt;&gt;&gt; ; Include chain topologies<br>&gt; &gt;&gt;&gt; #include "chnrc_A.itp"<br>&gt; &gt;&gt;&gt; #include "chnrc_B.itp"<br>&gt; &gt;&gt;&gt; #include "chnrc_C.itp"<br>&gt; &gt;&gt;&gt; #include "chnrc_D.itp"<br>&gt; &gt;&gt;&gt; #include "chnrc_E.itp"<br>&gt; &gt;&gt;&gt; <br>&gt; &gt;&gt;&gt; ; Include water topology<br>&gt; &gt;&gt;&gt; #include "spce.itp"<br>&gt; &gt;&gt;&gt; <br>&gt; &gt;&gt;&gt; #ifdef POSRES_WATER<br>&gt; &gt;&gt;&gt; ; Position restraint for each water oxygen<br>&gt; &gt;&gt;&gt; [ position_restraints ]<br>&gt; &gt;&gt;&gt; ;&nbsp; i funct&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>&gt; &gt;&gt; fcx&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>&gt; fcy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; fcz<br>&gt; &gt;&gt;&gt;&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp; <br>&gt; &gt;&gt; 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>&gt; 1000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1000<br>&gt; &gt;&gt;&gt; #endif<br>&gt; &gt;&gt;&gt; <br>&gt; &gt;&gt;&gt; ; Include generic topology for ions<br>&gt; &gt;&gt;&gt; #include "ions.itp"<br>&gt; &gt;&gt;&gt; <br>&gt; &gt;&gt;&gt; [ system ]<br>&gt; &gt;&gt;&gt; ; Name<br>&gt; &gt;&gt;&gt; Protein in water<br>&gt; &gt;&gt;&gt; <br>&gt; &gt;&gt;&gt; [ molecules ]<br>&gt; &gt;&gt;&gt; ; Compound&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; #mols<br>&gt; &gt;&gt;&gt; <br>&gt; &gt;&gt; <br>&gt; Protein_A&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1<br>&gt; &gt;&gt;&gt; <br>&gt; &gt;&gt; <br>&gt; Protein_B&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1<br>&gt; &gt;&gt;&gt; <br>&gt; &gt;&gt; <br>&gt; Protein_C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1<br>&gt; &gt;&gt;&gt; <br>&gt; &gt;&gt; <br>&gt; Protein_D&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1<br>&gt; &gt;&gt;&gt; <br>&gt; &gt;&gt; <br>&gt; Protein_E&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1<br>&gt; &gt;&gt;&gt; <br>&gt; &gt;&gt; <br>&gt; SOL&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 55419<br>&gt; &gt;&gt;&gt; <br>&gt; &gt;&gt; <br>&gt; NA+&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 54<br>&gt; &gt;&gt;&gt; <br>&gt; &gt;&gt;&gt; <br>&gt; &gt;&gt;&gt; Thanks &amp; regards,<br>&gt; &gt;&gt;&gt; Latha.<br>&gt; &gt;&gt;&gt; <br>&gt; &gt;&gt;&gt;&gt;&gt; plmallip@mail.uh.edu wrote:<br>&gt; &gt;&gt;&gt;&gt;&gt; <br>&gt; &gt;&gt;&gt;&gt;&gt; Hi Justin, thanks for your response. You are right. I get <br>&gt; &gt;&gt; the <br>&gt; &gt;&gt;&gt;&gt; following <br>&gt; &gt;&gt;&gt;&gt;&gt; message in the .job file<br>&gt; &gt;&gt;&gt;&gt;&gt; <br>&gt; &gt;&gt;&gt;&gt;&gt; "Warning: 1-4 interaction between 6151 and 6160 at <br>&gt; distance <br>&gt; &gt;&gt;&gt;&gt; 1.387 which <br>&gt; &gt;&gt;&gt;&gt;&gt; is larger than the 1-4 table size 1.000 nm<br>&gt; &gt;&gt;&gt;&gt;&gt; These are ignored for the rest of the simulation<br>&gt; &gt;&gt;&gt;&gt;&gt; This usually means your system is exploding"<br>&gt; &gt;&gt;&gt;&gt;&gt; <br>&gt; &gt;&gt;&gt;&gt;&gt; The atoms I restrained are 1703 &amp; 1712 (6151 &amp; 6160 atoms <br>&gt; &gt;&gt; are <br>&gt; &gt;&gt;&gt;&gt; also from <br>&gt; &gt;&gt;&gt;&gt;&gt; the same residues). The distance between them is 12 A <br>&gt; (there <br>&gt; &gt;&gt;&gt;&gt; are 14 <br>&gt; &gt;&gt;&gt;&gt;&gt; missing residues in between). My aim is to use distance <br>&gt; &gt;&gt;&gt;&gt; restraints, <br>&gt; &gt;&gt;&gt;&gt;&gt; without building missing residues in between. Is there any <br>&gt; &gt;&gt; way <br>&gt; &gt;&gt;&gt;&gt; I can <br>&gt; &gt;&gt;&gt;&gt;&gt; overcome this warning why minimizing the system?<br>&gt; &gt;&gt;&gt;&gt; <br>&gt; &gt;&gt;&gt;&gt; It's very hard to say, because none of us knows what's in <br>&gt; &gt;&gt; your <br>&gt; &gt;&gt;&gt;&gt; system, topology, <br>&gt; &gt;&gt;&gt;&gt; or how you built things :)<br>&gt; &gt;&gt;&gt;&gt; <br>&gt; &gt;&gt;&gt;&gt; You have some type of nasty steric clash that's driving <br>&gt; atoms <br>&gt; &gt;&gt;&gt;&gt; 6151 and 6160 <br>&gt; &gt;&gt;&gt;&gt; apart.&nbsp; Visualize the trajectory (.trr) to see if you <br>&gt; &gt;&gt; can <br>&gt; &gt;&gt;&gt;&gt; identify where things <br>&gt; &gt;&gt;&gt;&gt; start to break down.&nbsp; I don't know if the distance <br>&gt; &gt;&gt;&gt;&gt; restraint has anything to do <br>&gt; &gt;&gt;&gt;&gt; with the problem or not.<br>&gt; &gt;&gt;&gt;&gt; <br>&gt; &gt;&gt;&gt;&gt; Is there a problem with building missing residues?&nbsp; <br>&gt; That <br>&gt; &gt;&gt;&gt;&gt; might make life quite a <br>&gt; &gt;&gt;&gt;&gt; bit easier in the long run.<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; <br>&gt; &gt;&gt;&gt;&gt;&gt; Thanks &amp; regards,<br>&gt; &gt;&gt;&gt;&gt;&gt; Latha.<br>&gt; &gt;&gt;&gt;&gt;&gt; <br>&gt; &gt;&gt;&gt;&gt;&gt; <br>&gt; &gt;&gt;&gt;&gt;&gt;&nbsp;&nbsp;&nbsp; &gt; Dear colleagues,<br>&gt; &gt;&gt;&gt;&gt;&gt;&nbsp;&nbsp;&nbsp; &gt;<br>&gt; &gt;&gt;&gt;&gt;&gt;&nbsp;&nbsp;&nbsp; &gt; I need to use simple distance <br>&gt; restraints of <br>&gt; &gt;&gt; 12.5 A <br>&gt; &gt;&gt;&gt;&gt; between two CA atoms<br>&gt; &gt;&gt;&gt;&gt;&gt;&nbsp;&nbsp;&nbsp; &gt; of two residues. I am using the <br>&gt; following <br>&gt; &gt;&gt; lines in the <br>&gt; &gt;&gt;&gt;&gt; .itp file<br>&gt; &gt;&gt;&gt;&gt;&gt;&nbsp;&nbsp;&nbsp; &gt;<br>&gt; &gt;&gt;&gt;&gt;&gt;&nbsp;&nbsp;&nbsp; &gt; #ifdef DDISRES<br>&gt; &gt;&gt;&gt;&gt;&gt;&nbsp;&nbsp;&nbsp; &gt; [distance_restraints]<br>&gt; &gt;&gt;&gt;&gt;&gt;&nbsp;&nbsp;&nbsp; &gt; ;ai&nbsp;&nbsp; aj type index type' <br>&gt; &gt;&gt; low&nbsp;&nbsp; <br>&gt; &gt;&gt;&gt;&gt; up1&nbsp; up2&nbsp; fac<br>&gt; &gt;&gt;&gt;&gt;&gt;&nbsp;&nbsp;&nbsp; &gt; 1703 1712 1&nbsp;&nbsp;&nbsp; <br>&gt; &gt;&gt;&gt;&gt; 0&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp; 11.5&nbsp; <br>&gt; 12.0 <br>&gt; &gt;&gt;&gt;&gt; 12.5 1.0<br>&gt; &gt;&gt;&gt;&gt;&gt;&nbsp;&nbsp;&nbsp; &gt;<br>&gt; &gt;&gt;&gt;&gt;&gt;&nbsp;&nbsp;&nbsp; &gt; The first feww line of .mdp file for <br>&gt; &gt;&gt; minimisation of <br>&gt; &gt;&gt;&gt;&gt; protein alone is<br>&gt; &gt;&gt;&gt;&gt;&gt;&nbsp;&nbsp;&nbsp; &gt;<br>&gt; &gt;&gt;&gt;&gt;&gt;&nbsp;&nbsp;&nbsp; &gt; ; Preprocessing<br>&gt; &gt;&gt;&gt;&gt;&gt;&nbsp;&nbsp;&nbsp; &gt; ;<br>&gt; &gt;&gt;&gt;&gt;&gt;&nbsp;&nbsp;&nbsp; &gt; <br>&gt; &gt;&gt;&gt;&gt; <br>&gt; &gt;&gt; <br>&gt; title&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; ${MOL}<br>&gt; &gt;&gt;&gt;&gt;&gt;&nbsp;&nbsp;&nbsp; &gt; <br>&gt; &gt;&gt;&gt;&gt; <br>&gt; &gt;&gt; <br>&gt; cpp&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; /lib/cpp&nbsp;&nbsp;&nbsp; ;Preprocessor<br>&gt; &gt;&gt;&gt;&gt;&gt;&nbsp;&nbsp;&nbsp; &gt; <br>&gt; &gt;&gt;&gt;&gt; <br>&gt; &gt;&gt; <br>&gt; define&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = -DDISRES&nbsp; ;For cg, and also steep<br>&gt; &gt;&gt;&gt;&gt;&gt;&nbsp;&nbsp;&nbsp; &gt;<br>&gt; &gt;&gt;&gt;&gt;&gt; <br>&gt; &gt;&gt;&gt;&gt;&gt; You are not actually applying your distance <br>&gt; restraint.&nbsp; <br>&gt; &gt;&gt;&gt;&gt; If you have "#ifdef<br>&gt; &gt;&gt;&gt;&gt;&gt; DDISRES," then you would have to "define = -DDDISRES" in <br>&gt; the <br>&gt; &gt;&gt;&gt;&gt; .mdp file.&nbsp; <br>&gt; &gt;&gt;&gt;&gt;&gt; What<br>&gt; &gt;&gt;&gt;&gt;&gt; you probably meant to define was "#ifdef DISRES" in the <br>&gt; topology.&gt;&gt;&gt;&gt;&gt; <br>&gt; &gt;&gt;&gt;&gt;&gt;&nbsp;&nbsp;&nbsp; <br>&gt; &gt;&gt;&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>&gt; &gt;&gt;&gt;&gt; After minimisation, the restarined residues &amp; the adjacent<br>&gt; &gt;&gt;&gt;&gt;&gt;&nbsp;&nbsp;&nbsp; &gt; bonds break. This results in fragmnets <br>&gt; - <br>&gt; &gt;&gt; residues <br>&gt; &gt;&gt;&gt;&gt; alone, peptide bond<br>&gt; &gt;&gt;&gt;&gt;&gt;&nbsp;&nbsp;&nbsp; &gt; alone and the rest of the protein. The <br>&gt; &gt;&gt; distance <br>&gt; &gt;&gt;&gt;&gt; restrained residues<br>&gt; &gt;&gt;&gt;&gt;&gt;&nbsp;&nbsp;&nbsp; &gt; seems to try to move towards each <br>&gt; other (6.04 <br>&gt; &gt;&gt; A after <br>&gt; &gt;&gt;&gt;&gt; minimisation) and<br>&gt; &gt;&gt;&gt;&gt;&gt;&nbsp;&nbsp;&nbsp; &gt; this might have caused fragmentation. <br>&gt; I tried <br>&gt; &gt;&gt; to use <br>&gt; &gt;&gt;&gt;&gt; various upper and<br>&gt; &gt;&gt;&gt;&gt;&gt;&nbsp;&nbsp;&nbsp; &gt; lower values for bond length so as to <br>&gt; increase <br>&gt; &gt;&gt;&gt;&gt; flexibility. But, still I<br>&gt; &gt;&gt;&gt;&gt;&gt;&nbsp;&nbsp;&nbsp; &gt; end up in the fragments.<br>&gt; &gt;&gt;&gt;&gt;&gt;&nbsp;&nbsp;&nbsp; &gt;<br>&gt; &gt;&gt;&gt;&gt;&gt; <br>&gt; &gt;&gt;&gt;&gt;&gt; Bonds don't break in classical MD, this is just an <br>&gt; artifact <br>&gt; &gt;&gt; of <br>&gt; &gt;&gt;&gt;&gt;&gt; visualization,<br>&gt; &gt;&gt;&gt;&gt;&gt; probably from nasty steric clashes within your structure.<br>&gt; &gt;&gt;&gt;&gt;&gt; <br>&gt; &gt;&gt;&gt;&gt;&gt; <br>&gt; &gt;&gt;&gt;&gt;&gt;&nbsp;&nbsp;&nbsp; <br>&gt; &gt;&gt;&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Some <br>&gt; &gt;&gt;&gt;&gt; of the suggestions in the archive says VMD doesn't show<br>&gt; &gt;&gt;&gt;&gt;&gt;&nbsp;&nbsp;&nbsp; &gt; bonds if they r above threshold value. <br>&gt; When I <br>&gt; &gt;&gt; checked <br>&gt; &gt;&gt;&gt;&gt; the distances<br>&gt; &gt;&gt;&gt;&gt;&gt;&nbsp;&nbsp;&nbsp; &gt; between the atoms, one of them is <br>&gt; really long <br>&gt; &gt;&gt; CO-CA <br>&gt; &gt;&gt;&gt;&gt; bond 3.23 A<br>&gt; &gt;&gt;&gt;&gt;&gt;&nbsp;&nbsp;&nbsp; &gt; (normally its 1.59A). This means the <br>&gt; bond is <br>&gt; &gt;&gt; no longer <br>&gt; &gt;&gt;&gt;&gt; there.&gt;&nbsp; &gt;<br>&gt; &gt;&gt;&gt;&gt;&gt; <br>&gt; &gt;&gt;&gt;&gt;&gt; No, the bond is there, VMD just isn't smart enough to see <br>&gt; it <br>&gt; &gt;&gt;&gt;&gt; :)&nbsp; You are<br>&gt; &gt;&gt;&gt;&gt;&gt; probably well on your way to an explosion if you try to <br>&gt; &gt;&gt;&gt;&gt; constrain bond <br>&gt; &gt;&gt;&gt;&gt;&gt; lengths<br>&gt; &gt;&gt;&gt;&gt;&gt; with LINCS, however.<br>&gt; &gt;&gt;&gt;&gt;&gt; <br>&gt; &gt;&gt;&gt;&gt;&gt; -Justin<br>&gt; &gt;&gt;&gt;&gt;&gt; <br>&gt; &gt;&gt;&gt;&gt;&gt; <br>&gt; &gt;&gt;&gt;&gt;&gt; -----------------------------------------------------------<br>&gt; --<br>&gt; &gt;&gt; --<br>&gt; &gt;&gt;&gt;&gt; ---------<br>&gt; &gt;&gt;&gt;&gt;&gt; <br>&gt; &gt;&gt;&gt;&gt;&gt; _______________________________________________<br>&gt; &gt;&gt;&gt;&gt;&gt; gmx-users mailing list&nbsp;&nbsp;&nbsp; gmx-users@gromacs.org<br>&gt; &gt;&gt;&gt;&gt;&gt; http://www.gromacs.org/mailman/listinfo/gmx-users<br>&gt; &gt;&gt;&gt;&gt;&gt; Please search the archive at http://www.gromacs.org/search <br>&gt; &gt;&gt;&gt;&gt; before posting!<br>&gt; &gt;&gt;&gt;&gt;&gt; Please don't post (un)subscribe requests to the list. Use <br>&gt; &gt;&gt; the <br>&gt; &gt;&gt;&gt;&gt;&gt; www interface or send it to gmx-users-request@gromacs.org.<br>&gt; &gt;&gt;&gt;&gt;&gt; Can't post? Read http://www.gromacs.org/mailing_lists/users.php<br>&gt; &gt;&gt;&gt;&gt; <br>&gt; &gt;&gt;&gt;&gt; -- <br>&gt; &gt;&gt;&gt;&gt; ========================================<br>&gt; &gt;&gt;&gt;&gt; <br>&gt; &gt;&gt;&gt;&gt; Justin A. Lemkul<br>&gt; &gt;&gt;&gt;&gt; Graduate Research Assistant<br>&gt; &gt;&gt;&gt;&gt; Department of Biochemistry<br>&gt; &gt;&gt;&gt;&gt; Virginia Tech<br>&gt; &gt;&gt;&gt;&gt; Blacksburg, VA<br>&gt; &gt;&gt;&gt;&gt; jalemkul[at]vt.edu | (540) 231-9080<br>&gt; &gt;&gt;&gt;&gt; http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin<br>&gt; &gt;&gt;&gt;&gt; <br>&gt; &gt;&gt;&gt;&gt; ========================================<br>&gt; &gt;&gt;&gt;&nbsp;&nbsp; <br>&gt; &gt;&gt;&gt; _______________________________________________<br>&gt; &gt;&gt;&gt; gmx-users mailing list&nbsp;&nbsp;&nbsp; gmx-users@gromacs.org<br>&gt; &gt;&gt;&gt; http://www.gromacs.org/mailman/listinfo/gmx-users<br>&gt; &gt;&gt;&gt; Please search the archive at http://www.gromacs.org/search <br>&gt; &gt;&gt; before posting!<br>&gt; &gt;&gt;&gt; Please don't post (un)subscribe requests to the list. Use <br>&gt; the <br>&gt; &gt;&gt;&gt; www interface or send it to gmx-users-request@gromacs.org.<br>&gt; &gt;&gt;&gt; Can't post? Read http://www.gromacs.org/mailing_lists/users.php<br>&gt; &gt;&gt; <br>&gt;