Hi Ed, Thanks for this. I will create the report asap. Also, as mentioned in the earlier email, in my protein structure, I see a high backbone RMSD for first 5 residues and last 2 residues. When I use just the low energy structure I get one set of diffusion tensor data, which is quite different from what I get then i use the ensemble. I am abit puzzled about why this happens. Is it possible that the change from oblate to ellipsoid, when a 15 structure ensemble is used, might be due to over estimation of the tensor caused by having 15 different structures having highly flexible N and C termini? Regards, Vineet ________________________________________ Von: edward.dauvergne@xxxxxxxxx [edward.dauvergne@xxxxxxxxx]" im Auftrag von "Edward d'Auvergne [edward@xxxxxxxxxxxxx] Gesendet: Mittwoch, 14. Mai 2014 17:31 An: Panwalkar, Vineet Cc: relax-users@xxxxxxx Betreff: Re: Dear Vineet, The problem you are seeing is a simple bug - Dpar is simply being displayed twice. Could you please create a bug report for this? Simply fill out the details at https://gna.org/bugs/?func=additem&group=relax. I will then fix the problem and the solution will be available with the next relax release. The bug report will be useful for other users Cheers, Edward P. S. I can already see the problem in the relax source code. See if you can see it too in the display() function at http://www.nmr-relax.com/api/3.1/pipe_control.diffusion_tensor-pysrc.html#display. On 14 May 2014 17:00, Panwalkar, Vineet <v.panwalkar@xxxxxxxxxxxxx> wrote:
Dear users, I have recently carried out two seperate model-free analysis on a ~5 kDa protein with relaxation data measured at 600, 800 and 900 MHz. I used the fully automated analysis scripts for the analysis. In the first run, I used the lowest energy structure from my 15 structure ensemble. I get an oblate diffusion tensor. When I read the results and ask for diffusion_tensor.display() in relax, I get the table with all the data which is as follows... Diffusion type spheroid tm (s) 4.63577e-09 Diso (rad/s) 3.59523e+07 Da (rad/s) -3.36136e+07 Dpar (rad/s) 1.35432e+07 Dper (rad/s) 1.35432e+07 Dratio 0.287196 theta (rad) 0.633624 phi (rad) 2.61738 Fixed flag True What is confusing me is that, even though the Dratio is 0.287196, the individual Dpar and Dper values are exactly the same! I have seen the same when I run relax on a complex of the 5 kDa protein and a 1.1 kDa peptide, where relaxation data was recorded only for the protein. I ran another run, this time using the pdb file containing all the 15 lowest energy structures. Analysis of the same relaxation data with this structure gives me an ellipsoid diffusion model with following parameters.. Diffusion type ellipsoid tm (s) 4.63094e-09 Diso (rad/s) 3.59898e+07 Da (rad/s) 2.4722e+07 Dr 0.475806 Dx (rad/s) 1.59863e+07 Dy (rad/s) 3.9512e+07 Dz (rad/s) 5.24711e+07 alpha (rad) 2.59643 beta (rad) 2.087 gamma (rad) 0.438009 Fixed flag True I remember reading an old post where Ed mentions that relax is not set up to carry out analysis for residues with multiple NH bond vectors since it lies in the model-free theory itself. In the same post it is also mentioned that if that is the case then relax should give an error. Now, about 5 residues at the N terminus and 2 at the C terminus of my 43 residue protein are unstructured and show high backbone RMSD. It would be great if i get some help regarding the understanding of why in such a case, there was no error given from relax. And also, what measures are to be taken to run relax analysis on such a protein using the 15 structure ensemble as an input structure. Is this change from oblate to ellipsoid a case of over-estimation of the diffusion properties occuring due to highly flexible termini in the protein? Regards, Vineet ------------------------------------------------------------------------------------------------ ------------------------------------------------------------------------------------------------ Forschungszentrum Juelich GmbH 52425 Juelich Sitz der Gesellschaft: Juelich Eingetragen im Handelsregister des Amtsgerichts Dueren Nr. HR B 3498 Vorsitzender des Aufsichtsrats: MinDir Dr. Karl Eugen Huthmacher Geschaeftsfuehrung: Prof. Dr. Achim Bachem (Vorsitzender), Karsten Beneke (stellv. Vorsitzender), Prof. Dr.-Ing. Harald Bolt, Prof. Dr. Sebastian M. Schmidt ------------------------------------------------------------------------------------------------ ------------------------------------------------------------------------------------------------ _______________________________________________ relax (http://www.nmr-relax.com) This is the relax-users mailing list relax-users@xxxxxxx To unsubscribe from this list, get a password reminder, or change your subscription options, visit the list information page at https://mail.gna.org/listinfo/relax-users