<div class="aju"><span class="gD">Tsjerk</span>,<br></div><br>Thank you for suggestions!<br><br><br>Indeed the hight temperature ( I'm using 700K) which I use for enhansing sampling rate resulted in destabilisation of the secondary structure.<br>
<br>To prevent this I've used two slightly different aproaches based on the restraint. But in both cases I've used slightly soft restrains with fs=10-200.<br><br>The first aproach is the ussage of the posres applied on each backbone atom with fc=10-50. I've tested case with fc=50 and found that such restrains were very hight. I've not noticed any conformation sampling of my protein ( rmsd of backbone < 0.3 nm) during 10ns of such simulation. So I've decided to test a case with fc=10 ( under calculation)<br>
<br>The second approach is the application of the harmonic distance restrainse wich I've applied on each backbone atom pair in the CUTOFF radius of 1.0 nm. The Rc value was chosen because of my protein is the 7 buddle of alpha helices so this aproach could be usefull but exactly value for R have been chosen empirically. I've selected deviation value wich are equal to 1\2 of cutoff radius = 1.5 nm. I have not realise whaat I could obtain yet from this because I'm not quite sure about coccect values of such disres.<br>
<br><br>James<br><br><div class="gmail_quote">11 апреля 2012 г. 13:35 пользователь Tsjerk Wassenaar <span dir="ltr"><<a href="mailto:tsjerkw@gmail.com">tsjerkw@gmail.com</a>></span> написал:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hey :)<br>
<br>
I'd say a protein should be loosing structure at 700K. Can't even say<br>
the force field is wrong there, even though it hasn't been<br>
parameterized for such temperatures for certain. You'd have to check<br>
with experiments.<br>
Trying to do high-temperature simulations is fine. But trying to do<br>
high temperature simulations to get results that match a room/body<br>
temperature ensemble is completely bogus. If you aim to use this<br>
approach for enhanced sampling, you're in for some serious<br>
reparameterization. Part of the deal may indeed be adding<br>
long(er)-range distance restraints. However, the force constants to<br>
use will increase with temperature, and with very high temperatures<br>
the force constants may end up so high that they give rise to fast<br>
oscillations, which will require using a smaller time step.<br>
<br>
Just my 2 cents...<br>
<span class="HOEnZb"><font color="#888888"><br>
Tsjerk<br>
<br>
--<br>
Tsjerk A. Wassenaar, Ph.D.<br>
<br>
post-doctoral researcher<br>
Molecular Dynamics Group<br>
* Groningen Institute for Biomolecular Research and Biotechnology<br>
* Zernike Institute for Advanced Materials<br>
University of Groningen<br>
The Netherlands<br>
</font></span><div class="HOEnZb"><div class="h5">--<br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org">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<br>
www interface or send it to <a href="mailto:gmx-users-request@gromacs.org">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>