Re: [gmx-users] more than one peptide in one simulation box

2010-04-28 Thread Justin A. Lemkul
shahid nayeem wrote: Hi Justin Should I try to do position restraint at 500k and then full MD simulation. Always equilibrate under the conditions you wish to use prior to collecting data under those conditions. As for whether or not your force field is even valid at such high temperature

Re: [gmx-users] more than one peptide in one simulation box

2010-04-27 Thread shahid nayeem
Dear Mark 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

Re: [gmx-users] more than one peptide in one simulation box

2010-04-27 Thread Justin A. Lemkul
shahid nayeem wrote: Dear Mark 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.

Re: [gmx-users] more than one peptide in one simulation box

2010-04-27 Thread shahid nayeem
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 Energy minimization cpp = /usr/bin/cpp define = -DFLEX_SPC constraints = none integrator = steep nsteps = 3000 ; ; Energy

Re: [gmx-users] more than one peptide in one simulation box

2010-04-27 Thread Justin A. Lemkul
shahid nayeem wrote: 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 snip ref_t = 300 300 Here, you're equilibrating at 300 K... snip ref_t = 500 500 and here,

Re: [gmx-users] more than one peptide in one simulation box

2010-04-27 Thread Mark Abraham
On 27/04/2010 8:58 PM, Justin A. Lemkul wrote: shahid nayeem wrote: Dear Mark 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

Re: [gmx-users] more than one peptide in one simulation box

2010-04-27 Thread shahid nayeem
Hi Justin Should I try to do position restraint at 500k and then full MD simulation. shahid On 4/27/10, Justin A. Lemkul jalem...@vt.edu wrote: shahid nayeem wrote: My peptide is 26 residue alpha helix obtained from crystal structure .pdb file. I am posting energy minimization, position

Re: [gmx-users] more than one peptide in one simulation box

2010-04-27 Thread shahid nayeem
Hi Mark How one should be certain that this much trajectory is long enough to get coverged ensemble. Shahid On 4/27/10, Mark Abraham mark.abra...@anu.edu.au wrote: On 27/04/2010 8:58 PM, Justin A. Lemkul wrote: shahid nayeem wrote: Dear Mark Following your advice I started using three

Re: [gmx-users] more than one peptide in one simulation box

2010-04-27 Thread Mark Abraham
On 28/04/10 15:13, shahid nayeem wrote: Hi Mark How one should be certain that this much trajectory is long enough to get coverged ensemble. When the observables of interest aren't changing... This is a how long is a piece of string?-type question. Read some literature about simulations of

[gmx-users] more than one peptide in one simulation box

2010-04-22 Thread shahid nayeem
Dear All I am trying to study inter peptide interaction fpr which I need to put more than one peptide in one simulation box. I did it with genconf command but this inserts peptide in a regular ordered manner I want these to be in irregular disordered insertion. Even after using genconf , I tried

Re: [gmx-users] more than one peptide in one simulation box

2010-04-22 Thread Mark Abraham
On 23/04/10 13:16, shahid nayeem wrote: Dear All I am trying to study inter peptide interaction fpr which I need to put more than one peptide in one simulation box. I did it with genconf command but this inserts peptide in a regular ordered manner I want these to be in irregular disordered