Dear Tserk and the rest of GROMACS users,<br><br>Last time I measured the cosine content of different time intervals from the PCs of the whole trajectory. This time I did PCA for each time interval and measured the cc which is the right way I suppose. Maisuradze et al., 2009 claim that a CC value below 0.5 indicates that the protein-water system stops behaving as one with Brownian motions. Yet in none of my time intervals the CC dropped below 0.6. Can this happen for an equilibrated system?<br>
<br>It is known that RMSIP can be used to measure the overlap between essential subspaces. My question is can RMSIP<br>help me decide at which point to consider a system equilibrated from a thermodynamical point of view, namely the point at which the system stops doing Brownian motion? Does this point necessarily coincide with that where the RMSD is stabilized?<br>
<br>Our protein is co-complexed in another one in the crystal structure, which is believed to hold it in an "inactive state". However, we simulate the monomer state, hence we expect a large domain motion to occur in the first few nanosecond. Provided that the first 20ns are considered the time for equilibration, does it make sense to do PCA in this time interval?<br>
<br><br>Thanks in advance for your help.<br><br><br>Thomas<br><br><br><br><div class="gmail_quote">On 12 March 2012 15:44, Tsjerk Wassenaar <span dir="ltr"><<a href="mailto:tsjerkw@gmail.com">tsjerkw@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Thomas,<br>
<br>
Whether or not it makes sense to do PCA on the domain only depends on<br>
the question you ask. It may well make sense if you aim at<br>
characterizing the intra-domain motions. But be aware that you will<br>
view those motions within the context of the rest of the protein. It<br>
is quite likely that the internal motions depend on the interaction<br>
with the surroundings.<br>
<br>
Concerning the cosine content, there indeed seems to ba a transition.<br>
But that is already evident from your RMSD plot. A lot of stuff is<br>
still happening in the second half of the simulation. To see if the<br>
last part of your simulation is in a (local) equilibrium, you can<br>
calculate the average structure from that part, say the alst 10 or 20<br>
ns, and determine the RMSD against the average. The RMSD should then<br>
decrease towards the average and level off or oscillate around it to<br>
assure you've reached some sort of equilibrium. On the basis of that<br>
you could then decide whether it's worthwhile performing PCA on that<br>
part only.<br>
<br>
Cheers,<br>
<br>
Tsjerk<br>
<div><div></div><br clear="all"></div></blockquote></div><br>-- <br>
<p style="margin-bottom:0cm" align="LEFT">======================================================================</p>
<p style="margin-bottom:0cm" align="LEFT">Thomas Evangelidis</p>
<p style="margin-bottom:0cm" align="LEFT">PhD student</p>
<p style="margin-bottom:0cm" align="LEFT">Biomedical Research
Foundation, Academy of Athens</p>
<p style="margin-bottom:0cm" align="LEFT">4 Soranou Ephessiou , 115 27
Athens, Greece<br><br>email: <a href="mailto:tevang@bioacademy.gr" target="_blank">tevang@bioacademy.gr</a></p>
<p style="margin-bottom:0cm" align="LEFT"> <a href="mailto:tevang3@gmail.com" target="_blank">tevang3@gmail.com</a></p>
<p style="margin-bottom:0cm" align="LEFT"><br>website:
<a href="https://sites.google.com/site/thomasevangelidishomepage/" target="_blank">https://sites.google.com/site/thomasevangelidishomepage/</a></p>
<p style="margin-bottom:0cm" align="LEFT"><br>
</p>
<br>