<span class="gD">Hi Tsjerk!<br><br>First, I'd like also thanks you for yourš help.<br><br><br>Today I tried to make PCA of my X-ray data as well as comparison between results of such PCA and EDA ( PCA wich is based on the MD trajectory of the same protein).<br>
<br>In generaly I have no any problems with the X-ray PCA but I've forced with some during comparison of results of both PCAs due to the different atom numbers in both datasets ( X-ray structures consist of missing atoms ). As I understood this problem could be solved by selection atoms for the MD trajectory ( before PCA calculation of that data) wich are present in all X-ray structures. <br>
<br>Because of the X-ray dataset consist of some missing residues in the loops of structures I've had to reduce atom numbers in initial tpr topology.<br>So by means of tpbconvš I've made new .tpr file for this X-ray structures. This new TPR was based on theš .tpr file wich I've obtained from the MD of full-atomic model of protein. In new TPR I include only mainchain atoms as well as residues presented in the X-ray structures ( I've defined it by means of index.ndx file ). Does this aproach correct ? Are there any extra ways to obtaint TPR file for my X-ray datasetš without redusing topology of the MD structure ?<br>
<br>Also I'd like to ask some more about PCA results.<br><br>1) Firstly, what exactly is the new compresed trajectory made by<br>-filt filtered.xtc<br><br>In case of x-ray PCA this trajectory correspond to the initial numbers of pdb's files but in case of MD_based PCA this trajectory consist of redused number of atoms.<br>
<br>2) Also I'd like to know what actyally is the <br>-extr extremePCA.pdb ?<br><br>In case of MD_based PCA I've obtaned 20 different extreme.pdb trajectories where 20 was the number of calculated Principal components.<br>
<br>But in case of X-ray PCA I've obtained only one such file wich was similar to visualisation of the motion along softest PC althrough I've calculated 10 eigenvectors. Why in case of PCA I've obtain only one such file and what exactly this trajectory is ? Is there any extra methods to visualise motions along specified components ( not for ensembles of components )?<br>
<br>3) How I could specify exactly number of PC in the projection graphs ? AS I understood 2d and 3d projections are made along the -first and -last components. How I can make such projections based on two/ three another specified components (e.g along 2 and 7 modes from 20 calculated) ?<br>
<br><br>Thanks for help,<br><br><br>James<br><br><br><br></span><br><div class="gmail_quote">28 ÍÁÒÔÁ 2012šÇ. 18:02 ÐÏÌØÚÏ×ÁÔÅÌØ 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">Hi Thomas,<br>
<div class="im"><br>
> Thanks for all the clarifications about PCA you make on the mailing list!<br>
<br>
</div>Thank you for the appreciation :)<br>
<div class="im"><br>
> I have a question about the commandlines you wrote. Why do you use the .tpr<br>
> file with the "-s" flag? Is it because you want to compare the<br>
> mass-wheighted covariance matrices? I use to calculate the covariance<br>
> matrices by giving to g_covar a .pdb file with the "-s" flag and then<br>
> calculate the RMSIP without giving any structure file. I guess no masses are<br>
> used in that covariance analysis, right? Do you recommend using atom masses<br>
> for PCA in general?<br>
<br>
</div>Well, I admit that in most cases I don't use mass-weighting myself.<br>
Unless you include hydrogens, it also doesn't matter much, as the<br>
masses are not very different. Only if you want to calculate<br>
frequencies, e.g. to connect to NMA and/or IR spectroscopy you would<br>
really need masses.<br>
If you use a .pdb or .gro file, you don't get mass-weighting. And<br>
you're right that for calculating the RMSIP, and the subspace overlap,<br>
and the martix of inner products, you don't need a structure filel,<br>
but only the eigenvectors, and possibly the eigenvalues.<br>
<br>
Cheers,<br>
<br>
Tsjerk<br>
<div class="HOEnZb"><div class="h5"><br>
<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>
--<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>