<div>My peptide is 26 residue alpha helix obtained from crystal structure  .pdb file. I am posting energy minimization, position restarint and full MD simulation .mdp file</div>
<div> </div>
<div>Energy minimization</div><font size="2">
<p>cpp = /usr/bin/cpp</p>
<p>define = -DFLEX_SPC</p>
<p>constraints = none</p>
<p>integrator = steep </p>
<p>nsteps = 3000</p>
<p>;</p>
<p>; Energy minimizing stuff</p>
<p>;</p>
<p>emtol = 1000</p>
<p>emstep = 0.01</p>
<p>nstcomm = 1</p>
<p>ns_type = grid</p>
<p>rlist = 0.9</p>
<p>coulombtype = PME</p>
<p>rcoulomb = 0.9</p>
<p>rvdw = 0.9</p>
<p>fourierspacing = 0.14</p>
<p>fourier_nx = 0</p>
<p>fourier_ny = 0</p>
<p>fourier_nz = 0</p>
<p>pme_order = 4</p>
<p>ewald_rtol = 1e-5</p>
<p>optimize_fft = yes</p>
<p>Tcoupl = no</p>
<p>Pcoupl = no</p>
<p>gen_vel = no</p>
<div></div></font>Position_restraint.mdp
<div> </div><font size="2">
<p>cpp = /usr/bin/cpp</p>
<p>define = -DPOSRES</p>
<p>constraints = all-bonds</p>
<p>constraintalgorithm = LINCS</p>
<p>integrator = md</p>
<p>dt = 0.002 ; ps !</p>
<p>nsteps = 25000 ; total 50 ps.</p>
<p>nstcomm = 1</p>
<p>nstxout = 500</p>
<p>nstvout = 1000</p>
<p>nstfout = 0</p>
<p>nstlog = 10</p>
<p>nstenergy = 10</p>
<p>nstlist = 10</p>
<p>ns_type = grid</p>
<p>rlist = 0.9</p>
<p>coulombtype = PME</p>
<p>rcoulomb = 0.9</p>
<p>rvdw = 0.9</p>
<p>fourierspacing = 0.14</p>
<p>fourier_nx = 0</p>
<p>fourier_ny = 0</p>
<p>fourier_nz = 0</p>
<p>pme_order = 4</p>
<p>ewald_rtol = 1e-5</p>
<p>optimize_fft = yes</p>
<p>; Berendsen temperature coupling is on in two groups</p>
<p>Tcoupl = berendsen</p>
<p>tc-grps = Protein Non-protein </p>
<p>tau_t = 0.1 0.1 </p>
<p>ref_t = 300 300 </p>
<p>; Energy monitoring</p>
<p>energygrps = Protein Non-protein</p>
<p>; Pressure coupling is not on</p>
<p>Pcoupl = no</p>
<p>tau_p = 0.5</p>
<p>compressibility = 4.5e-5</p>
<p>ref_p = 1.0</p>
<p>; Generate velocites is on at 300 K.</p>
<p>gen_vel = yes</p>
<p>gen_temp = 300.0</p>
<p>gen_seed = 173529</p>
<div></div></font>Full_MD.mdp
<div><font size="2">
<p>cpp = /usr/bin/cpp</p>
<p>constraints = all-bonds</p>
<p>integrator = md</p>
<p>dt = 0.002 ; ps !</p>
<p>nsteps = 5000000 ; total 10000 ps.</p>
<p>nstcomm = 1</p>
<p>nstxout = 5000</p>
<p>nstvout = 40000</p>
<p>nstfout = 0</p>
<p>nstlog = 500</p>
<p>nstenergy = 500</p>
<p>nstlist = 10</p>
<p>ns_type = grid</p>
<p>rlist = 0.9</p>
<p>coulombtype = PME</p>
<p>rcoulomb = 0.9</p>
<p>rvdw = 0.9</p>
<p>fourierspacing = 0.14</p>
<p>fourier_nx = 0</p>
<p>fourier_ny = 0</p>
<p>fourier_nz = 0</p>
<p>pme_order = 4</p>
<p>ewald_rtol = 1e-5</p>
<p>optimize_fft = yes</p>
<p>; Berendsen temperature coupling is on in two groups</p>
<p>Tcoupl = berendsen</p>
<p>tc-grps = Protein Non-protein</p>
<p>tau_t = 0.1 0.1 </p>
<p>ref_t = 500 500 </p>
<p>; Energy monitoring</p>
<p>energygrps = Protein Non-protein</p>
<p>; Isotropic pressure coupling is now on</p>
<p>Pcoupl = berendsen</p>
<p>Pcoupltype = isotropic</p>
<p>tau_p = 0.5</p>
<p>compressibility = 4.5e-5</p>
<p>ref_p = 1.0</p>
<p>; Generate velocites is off at 500 K.</p>
<p>gen_vel = no</p>
<p>gen_temp = 500.0</p>
<p>gen_seed = 173529</p>
<p> shahid</p></font></div>
<div><br><br> </div>
<div><span class="gmail_quote">On 4/27/10, <b class="gmail_sendername">Justin A. Lemkul</b> &lt;<a href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a>&gt; wrote:</span>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid"><span class="q"><br><br>shahid nayeem wrote:<br>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">Dear Mark<br>Following your advice I started using three peptide in one simulation box. Iwas able to add these with genconf as previously in ordered manner, generated .gro with genconf, solvated it and after energy minimization I did MD run for 10ns. Everything ran well. In the end when I see the trajectory I find unfolding of the original chain but the two additional peptide introduced through genconf show appearance of new secondary structures. Even in these two the secondary structure do not develop at the same point. Why the three equivalent peptide behave differently in similar environment. How can I explain this observation. why the first peptide does not show any new secondary structure. Sholud I go with higher number of molecule. Will it make any difference if peptides are added in disordered manner and then simulated.<br>
</blockquote><br></span>Initial orientation should likely have nothing to do with it.  Perhaps this is even the proper behavior for whatever your peptide is.  Is its structure dynamic?  Is the size of your peptides large enough to even believe that they would be intrinsically stable?  Many model peptides, in isolation, have very transient structures.<br>
<br>It could also be that your simulation parameters are poorly chosen, so the force field is breaking down.  If you want comments on your .mdp file, please post it.<br><br>-Justin<br><br>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">Shahid <br>
<div><span class="e" id="q_1283eed72006d203_2"><br> On 4/23/10, *Mark Abraham* &lt;<a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:Mark.Abraham@anu.edu.au" target="_blank">Mark.Abraham@anu.edu.au</a> &lt;mailto:<a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:Mark.Abraham@anu.edu.au" target="_blank">Mark.Abraham@anu.edu.au</a>&gt;&gt; wrote:<br>
<br>   On 23/04/10 13:16, shahid nayeem wrote:<br><br>       Dear All<br>       I am trying to study inter peptide interaction fpr which I need<br>       to put<br>       more than one peptide in one simulation box. I did it with genconf<br>
       command but this inserts peptide in a regular ordered manner I want<br>       these to be in irregular disordered insertion. Even after using<br>       genconf<br><br><br>   Well that&#39;s a difficult and atypical scenario. genconf -shuffle will<br>
   allow you to stack the same peptide in a regular array with random<br>   rotations of the whole box. Then you can solvate, equilibrate and<br>   run MD at a high temperature to give yourself a quasi-disordered<br>   starting state.<br>
<br>       , I tried to proceed furthe after solvation with spc water. The<br>       energy<br>       minimization (steepest descent) failed to converge even after<br>       5000 steps<br>       and theirafter position restraint dynamics failed giving<br>
       segmentation<br>       fault. Introducing more peptide after generating .gro with -ci -nmol<br>       gives error showing more than one residue in insert molecule.<br>       Please help me and write  commands which I should follow.<br>
<br><br>   No, because that&#39;s an impossible task. We can&#39;t begin to guess the<br>   reasons for things failing without seeing the actual output (was the<br>   EM energy large and negative? what was the actual error message<br>
    from -ci -nmol?).<br><br>   You should be careful to start with a small test case so that you<br>   can learn the workflow with a manageable problem. Can you get a<br>   single peptide to equilibrate? Two stacked peptides? It is best to<br>
   learn to walk before trying to run :-)<br><br>   Mark<br>   --     gmx-users mailing list    <a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
</span></div>   &lt;mailto:<a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>&gt;<span class="q"><br>   <a onclick="return top.js.OpenExtLink(window,event,this)" 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 onclick="return top.js.OpenExtLink(window,event,this)" href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before<br>   posting!<br>   Please don&#39;t post (un)subscribe requests to the list. Use the www<br>
   interface or send it to <a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br></span>   &lt;mailto:<a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>&gt;.<span class="q"><br>
   Can&#39;t post? Read <a onclick="return top.js.OpenExtLink(window,event,this)" href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br><br><br></span></blockquote>
<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]<a onclick="return top.js.OpenExtLink(window,event,this)" href="http://vt.edu/" target="_blank">vt.edu</a> | (540) 231-9080<br>
<a onclick="return top.js.OpenExtLink(window,event,this)" href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br><br>======================================== 
<div><span class="e" id="q_1283eed72006d203_8"><br>-- <br>gmx-users mailing list    <a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a onclick="return top.js.OpenExtLink(window,event,this)" 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 onclick="return top.js.OpenExtLink(window,event,this)" href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before posting!<br>
Please don&#39;t post (un)subscribe requests to the list. Use the www interface or send it to <a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can&#39;t post? Read <a onclick="return top.js.OpenExtLink(window,event,this)" href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br></span></div></blockquote>
</div><br>