Dear Sir, <br><br> May be this protein is a disulphide bonded dimer where both the monomers <br>are identical and only 1 disulphide bond is there between them. But other proteins <br>were not like that. May be that is the MAJOR reason for this. but strangely temp, pres, <br>
vol, density, kinetic, potential and total energy is stable. STRANGE!!<br><br>Thanks for Suggestions<br><br>With regards<br>M. Kavyashree<br><br><br><div class="gmail_quote">On Mon, Jun 13, 2011 at 6:09 PM, Justin A. Lemkul <span dir="ltr"><<a href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div><div></div><div class="h5"><br>
<br>
Kavyashree M wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Dear Sir,<br>
<br>
What does gmxcheck tell you about the .edr file that is giving weird<br>
results? Do the plots look normal? Perhaps a frame got corrupted<br>
somewhere along the way. The screen output should print how many<br>
frames were considered in the analysis; if it does not match your<br>
expectations based on nstenergy and the length of the simulation<br>
then something went wrong.<br>
<br>
<br>
gmxcheck gave proper number of frames corresponding to 100ns, other plots<br>
temp, pressure, volume, density look fine<br>
<br>
Most analyses do not need re-imaging. Keeping the protein within<br>
the confines of one unit cell is typically just a convenience for<br>
visualization.<br>
<br>
<br>
Ok.<br>
<br>
I had done simulations of four similar protein using the same mdp file. in one of them<br>
the minimum distance between the periodic images went near 0.9nm, I had used 1.0nm<br>
as distance between protein atoms and box wall. and cut offs were 1.0nm for vdw and 1.4nm<br>
for columb. Till some 17ns the minimum distance was above 2nm the gradually there was a<br>
dip after around 20-25ns. Now I ran 100ns simulation and I have to discard this trajectory<br>
because of this error. I thought distance of 2nm between protein atoms was enough as 1.4nm<br>
was the max cutoff.<br>
<br>
How can we know prior to starting the simulation that we may get some such errors for using<br>
such a parameter. I could have used larger box size but it will increase the time. Is it trial and<br>
error basis to find out the optimum box size?<br>
<br>
</blockquote>
<br></div></div>
Generally, no, you don't have to waste lots of time fiddling with the box size before you find the right one. The choice is based on the cutoffs and the nature of the system. For a well-folded, stable protein, I see no reason why what you've done isn't appropriate (unless you've set up the box incorrectly and only think you've set certain dimensions). For disordered proteins or those capable of large conformational changes, then the box needs to be large to accommodate these possible motions.<br>
<br>
The only other possibility I can think of is that the starting configuration compressed a lot over time, shrinking the box. I don't know why this would happen for a protein in water, but I suppose anything is possible. Most condensed phase systems should not be very compressible, but any such change would be obvious from plotting the volume over time. If it is stable, then the protein must be doing something unexpected.<br>
<font color="#888888">
<br>
-Justin</font><div><div></div><div class="h5"><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]<a href="http://vt.edu" target="_blank">vt.edu</a> | (540) 231-9080<br>
<a 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 <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a 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 href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
Please don't post (un)subscribe requests to the list. Use the www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
</div></div></blockquote></div><br>