[Chimera-users] B-Factor Difference Plot

[Elaine Meng]

Hi Amanda,
You could try various types of normalization based on the mean and/or spread of values, but I don’t know what is most appropriate biophysically for comparing B-factors between different crystal structures.  You could try asking on some other mailing lists frequented by crystallographers (maybe PDB-L, see https://lists.sdsc.edu/mailman/listinfo/pdb-l or CCP4 bulletin board, see http://www.ccp4.ac.uk/ccp4bb.php ) if you don’t get advice on this list.  My only other idea is to find papers that make similar comparisons and read what those authors did, but I don’t have any specific knowledge of such papers.
Best,
Elaine

> On Aug 28, 2017, at 3:34 PM, Amanda Evans Constantinides <aconstantinides1 at student.gsu.edu> wrote:
> 
> Hi Elaine,
> Thanks for your detailed reply. I successfully generated Chimera output files containing lists of average B factor values, aligned the residues, plotted residue vs. B factor difference, and saw the graphs existed mainly on one side or the other of the x axis, but were not distributed across the x axis as expected based on what I saw using B putty. I tried also using Baverage in CCP4 and made Main Chain B factor difference plots and saw the same thing. My question, hopefully written more clearly, is how to normalize between the structures in comparison. I am looking at the equation for b-factor to see how I can accomplish this. If you have any ideas or know anyone who does, please send them my way!
> 
> Best regards,
> 
> Amanda Constantinides
> Graduate Research Assistant
> Georgia State University 
> From: Elaine Meng <meng at cgl.ucsf.edu>
> Sent: Monday, August 28, 2017 5:37:36 PM
> To: Amanda Evans Constantinides
> Cc: chimera-users at cgl.ucsf.edu
> Subject: Re: [Chimera-users] B-Factor Difference Plot
>  
> Hi Amanda,
> Visually comparing the putty models (worms in Chimera) between the two structures is not meaningful unless you are sure you mapped the same value to the same radius in both structures.  Even so, it is hard to see in this way.  
> 
> The more direct way is to just look at the PDB file in a text editor and see what the B-factor values are in the corresponding residues.  I don’t know how you are calculating the differences, but you would have to make sure your script or whatever is mathematically correct and that it is also correctly pairing the residues between the two structures.  I.e. the residue numbering might be different in the two structures.
> 
> Now, comparing all the atomic values is probably too hard, so you could use Chimera to output the average values per residue so that you have fewer numbers to compare.  In the Render by Attribute dialog (the same one used to make the worms display) you could choose File… Save Attributes and then choose to save the attribute of “residues” named “average->bfactor” to a file.  You would need to do this for each of the two structures.  So that it doesn’t include other stuff besides the protein,  you could use command “select protein” and then choose the option to “Restrict save to current selection” (or alternatively, delete all the stuff you don’t care about for this comparison like solvent before saving the values).
> <https://na01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.rbvi.ucsf.edu%2Fchimera%2Fdocs%2FContributedSoftware%2Frender%2Frender.html&data=02%7C01%7Caconstantinides1%40student.gsu.edu%7Ce0956536835143ee7b8d08d4ee5d0094%7C704d822c358a47849a1649e20b75f941%7C0%7C0%7C636395530590758186&sdata=Q1j%2BbcqhfE3KREkJ8xy9rchlkFYvZav1ikqIy1XeRtg%3D&reserved=0>
> 
> Now, if you are sure your script is mathematically correct and also that the correct residues are paired between the two structures, then maybe your question is really how to compare B-factors between two structures when they seem to be on a different scale or systematically offset.  The question of whether or how best to normalize the values between different crystal structures is not a Chimera question per se, but a more general biophysical one to which I do not know the answer.
> 
> Best,
> Elaine
> ----------
> Elaine C. Meng, Ph.D. 
> UCSF Chimera(X) team
> Department of Pharmaceutical Chemistry
> University of California, San Francisco
> 
> 
> 
> > On Aug 28, 2017, at 9:00 AM, Amanda Evans Constantinides <aconstantinides1 at student.gsu.edu> wrote:
> > 
> > Hello, 
> > I am plotting the B-factor differences between two similar structures. I do not understand why the plot is almost entirely on the positive side of the X-axis. If I open B factor putty models in PYMOL and in Chimera, there are some visible areas that should translate to negative B-factor differences on a quantitative plot. There is no way one structure is entirely more flexible than the other. Can you help explain this to me? Is there away to adjust the zero point to more accurately reflect these differences on the plot? 
> > Thank you, 
> > Amanda Constantinides
> > Graduate Research Assistant
> > Georgia State University
> 
> _______________________________________________
> Chimera-users mailing list: Chimera-users at cgl.ucsf.edu
> Manage subscription: http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users