<div style="font-size: 16px; font-family: 'times new roman'"><div><br /></div><div>On 01/01/10, <b class="name">&quot;Justin A. Lemkul&quot; </b>&lt;jalemkul@vt.edu&gt; wrote:</div><blockquote cite="mid:4B3D1E21.8090601@vt.edu" class="iwcQuote" style="padding-left: 13px; margin-left: 0px; border-left: #00f 1px solid" type="cite"><br /><br />Arik Cohen wrote:<br />&gt;No, sorry for the confusion. The images are only from a simulation of one protein(tmRBP_Unliganded, PDB: 2FN9). The problem seen with BPTI was a bit different in a way that only 1 C-alpha was &quot;detached&quot; (cell size ?).<br />&gt;<br /><br />If there is supposed to be only a single protein in the images you've shown (and after having a look at the 2FN9 structure from the PDB), it seems pretty clear to me that part of your protein has simply crossed a periodic boundary and trjconv -pbc mol (or some such command) should fix it.<br /></blockquote><div>Indeed. Further, note that the reappearance of a subset of the seemingly detached atoms back with the main group is totally normal. With domain decomposition in GROMACS 4, such &quot;broken&quot; molecules can be written to the trajectory. Judicious use of trjconv fixes the appearance of this non-problem.</div><div ><br/></div><div>Mark</div><div ><br/></div><blockquote cite="mid:4B3D1E21.8090601@vt.edu" class="iwcQuote" style="padding-left: 13px; margin-left: 0px; border-left: #00f 1px solid" type="cite"><div class="mimepart text plain"><br />&gt;Arik<br />&gt;<br />&gt;Justin A. Lemkul wrote:<br />&gt;&gt;<br />&gt;&gt;<br />&gt;&gt;Arik Cohen wrote:<br />&gt;&gt;&gt;Which two proteins ? I have at least in beginning only one protein which some how is divided into two along the calculation.<br />&gt;&gt;&gt;Any way I'll try both increasing the cell and fix it with trjconv.<br />&gt;&gt;&gt;<br />&gt;&gt;<br />&gt;&gt;Quoting your original message:<br />&gt;&gt;<br />&gt;&gt;&quot;While running a simple MD simulation with both a small protein such as BPTI and a larger one such as tmRBP_Unliganded_2FN9.pdb...&quot;<br />&gt;&gt;<br />&gt;&gt;I assumed that what I was seeing in the images was a set of two proteins.  My concern was that you defined a box relative to the larger protein, then inserted the smaller one (BPTI?), leaving insufficient space in the box to satisfy the minimum image convention.  If that's not what we're looking at, then that'd be useful to know :)<br />&gt;&gt;<br />&gt;&gt;If you have a single protein, &quot;divided into two&quot; then the problem is almost certainly a simple periodicity artifact.  Bonds do not break in a normal MD calculation (in fact they can't using the standard MM approximations).<br />&gt;&gt;<br />&gt;&gt;-Justin<br />&gt;&gt;<br />&gt;&gt;&gt;Thanks a lot<br />&gt;&gt;&gt;<br />&gt;&gt;&gt;Arik<br />&gt;&gt;&gt;<br />&gt;&gt;&gt;Justin A. Lemkul wrote:<br />&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;Arik Cohen wrote:<br />&gt;&gt;&gt;&gt;&gt;I'm using dodecahedron -d 0.7<br />&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;Was that distance specified with respect to both of the protein molecules in the unit cell?  You can check for spurious PBC interactions with g_mindist -pi. Anyway, I'd be curious to see how you do with trjconv.<br />&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;-Justin<br />&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;Justin A. Lemkul wrote:<br />&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;Arik Cohen wrote:<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;Hi,<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;I have not tried yet to fix it with trjconv which I will . Attached is a picture with 4 snapshots taken from the simulation. The C-alphas in question are emphasized with red color.<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;Is your unit cell sufficiently large?  It looks like the C-alphas indicated are simply crossing the periodic boundary on the &quot;left&quot; of the frame and interacting with the protein molecule in the &quot;right&quot; of the frame, which would indicate to me that the unit cell is too small and you're seeing spurious PBC interactions (i.e., violation of the minimum image convention).<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;&gt;Thanks<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;Arik<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;Justin A. Lemkul wrote:<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;Arik Cohen wrote:<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;Hi,<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;Sorry to bother you again ,but its not only a periodic effect since only *some of the atoms* in the  &quot;Detached&quot; group are vanishing from this group and reappearing in the main protein group. The rest of the atoms are either always in the detached or the main group.<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;In addition, the &quot;detached&quot; group includes three segments of the protein(8 residues(126-131), 8 residues(157-164) and 4 residues186-189).<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;From your description, this sounds exactly like a periodicity problem - some of the atoms are crossing the periodic boundary and are appearing in strange locations.  Have you even tried trjconv to fix it?  That would be useful information, as I see that Mark long ago also suggested the same sort of fix.<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;It is hard for me to envision what you are seeing.  It would be enormously helpful if you could post images (screenshots, etc) of the problematic structures to get a more expedient resolution.<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;-Justin<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;Thanks a lot<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;Arik<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;Justin A. Lemkul wrote:<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;Arik Cohen wrote:<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;Hi,<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;With regards to your question I do see some periodicity in which for a section of time in the trajectory some of the Calphas in the &quot;detached group&quot; are vanishing from it and reappear in the main protein.<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;In addition,<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;I would appreciate as before any suggestion you might have in the matter.<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;If this is just a periodicity artifact, fix it with trjconv.<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;-Justin<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;Thanks<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;Arik<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;Mark Abraham wrote:<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;Arik Cohen wrote:<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;Hi,<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;Thanks for answering so quickly !. Apparently whole residues have detached from the protein.<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;So... like I asked last time, are you seeing a periodicity artefact? &quot;Detached&quot; covers a whole gamut of possibilities.<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;Another strange thing that happens in pyMol and VMD is that when I select an atom or a residue in the detached group the selection appears twice: one in the detached group and one in the main part.<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;If you've got atoms duplicated, then it sounds like something's going wrong with how they're interpreting the structure file, or how you're manipulating it afterwards. Either way, it's not a problem for the GROMACS mailing list unless you can demonstrate the atoms are duplicated in the structure file (which they aren't!).<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;Mark<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;Arik<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;Mark Abraham wrote:<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;Arik Cohen wrote:<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;Dear GROMACS users,<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;While running a simple MD simulation with both a small protein such as BPTI and a larger one such as tmRBP_Unliganded_2FN9.pdb, I'm encountering an odd situation in which one (in the case of BPTI) or several Calphas (in the later case) are &quot;detaching them selfs&quot; from the main group.<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&quot;main group&quot; of what? Do the atoms bound to them move also? Are you seeing a periodicity artefact?<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;Mark<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;The problem appeared only after adding salt to the simulation(at least in the case of BPTI).<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;I would appreciate any suggestions and comments on the matter.<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;Thanks<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;Arik<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;The run files are:<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;*em.mdp:*<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; title               =  tmRBP_Unliganded_2FN9 Minimization<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;integrator          =  steep      ; (steep)using steepest descent<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;nsteps              =  50000<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;nstlist             =  1<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;rlist               =  1.0<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;coulombtype         =  PME<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;rcoulomb            =  1.0<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;vdw-type            =  cut-off<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;rvdw                =  1.0<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;nstenergy           =  10<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;emtol               =  5.0 ; tolerance kJ/(Mol -1 nm-1) instead of 10.0<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;*pr.mdp<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;*<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;title               =  tmRBP_Unliganded_2FN9 PR<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;integrator          =  md<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;nsteps              =  50000<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;dt                  =  0.002 ;(in ps) doing a 100ps traj.<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;constraints         =  all-bonds<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;nstlist             =  10 ; neighbour list updates every number of steps<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;rlist               =  1.0<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;coulombtype         =  PME<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;rcoulomb            =  1.0<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;vdw-type            =  cut-off<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;rvdw                =  1.0<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;tcoupl              =  Berendsen<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;tc-grps             =  Protein non-protein<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;tau-t               =  0.1 0.1<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;ref-t               =  298 298<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;Pcoupl              =  Berendsen<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;tau-p               =  1.0<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;compressibility     =  5e-5 5e-5 5e-5 0 0 0<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;ref-p               =  1.0<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;nstenergy           =  100<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;define              =  -DPOSRES ; include posre.itp(position restraint) file<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;*run.md<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;*title               =  tmRBP_Unliganded_2FN9<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;integrator          =  md<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;nsteps              =  300000<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;dt                  =  0.001<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;constraints         =  all-bonds<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;nstlist             =  10<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;rlist               =  1.0<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;coulombtype         =  PME<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;rcoulomb            =  1.0<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;vdw-type            =  cut-off<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;rvdw                =  1.0<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;tcoupl              =  V-rescale  ;V-rescale<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;tc-grps             =  Protein non-protein<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;tau-t               =  0.8 0.8<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;ref-t               =  298 298<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;nstxout             =  1000<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;nstvout             =  1000<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;nstxtcout           =  1000<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;nstenergy           =  1000<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;The runs commands are(integrated inside a C++ code):<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;SysCommand1 = &quot;echo 6 | pdb2gmx -f &quot; + FileName + &quot; -water tip3p&quot;;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; system(&quot;editconf -f conf.gro -bt dodecahedron -d 0.7 -o box.gro&quot;);<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;system(&quot;genbox -cp box.gro -cs spc216.gro -p topol.top -o solvated.gro&quot;);<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;minimization:<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;--------<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; if(Mode == &quot;NoSalt&quot;)<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;    {<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;     system(&quot;grompp -f MDP/em.mdp -p topol.top -c solvated.gro -o em.tpr&quot;);<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;      //system(&quot;mpirun -np 4 mdrun -v -deffnm em&quot;);<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;    }<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;  if(Mode == &quot;WithSalt&quot;)<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;    {<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;      system(&quot;grompp -f MDP/em.mdp -p topol.top -c solvated.gro -o em.tpr&quot;);<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;           system(&quot;mpirun -np 4 mdrun -v -deffnm em&quot;);<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;    }<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;Salting:<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;--------<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; system(&quot;echo 12 | genion -s em.tpr -conc 0.1 -neutral -o solvated.gro&quot;);<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; pr:<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;----<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;system(&quot;grompp -f MDP/prmd.mdp -p topol.top -c em.gro -o pr.tpr&quot;);<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;  /* The actual run*/<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;  system(&quot;mpirun -np 4 mdrun -v -deffnm pr&quot;); <br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;------------------------------------------------------------------------ <br />&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;&gt;<br />&gt;&gt;&gt;<br />&gt;&gt;<br />&gt;<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 _zipIdx="1001" 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 />-- <br />gmx-users mailing list    gmx-users@gromacs.org<br /><a _zipIdx="1001" 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 _zipIdx="1001" 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 gmx-users-request@gromacs.org.<br />Can't post? Read <a _zipIdx="1001" href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br /></div></blockquote></div>