[ccp4bb] B-factor distributions
Posted on behalf of Gert Vriend: This article didn't make it to Nature Methods (...), but might be of interest to theoreticians interested in B-factor distributions: https://journals.calstate.edu/pump/article/view/2409/2168 Gert Vriend To unsubscribe from the CCP4BB list, click the following link: https://www.jiscmail.ac.uk/cgi-bin/WA-JISC.exe?SUBED1=CCP4BB&A=1 This message was issued to members of www.jiscmail.ac.uk/CCP4BB, a mailing list hosted by www.jiscmail.ac.uk, terms & conditions are available at https://www.jiscmail.ac.uk/policyandsecurity/
Re: [ccp4bb] B-factor standardization
> If I am not wrong, I remember that someone proposed to standardize > B-factors of protein atoms as “BS = B - Bave”, where Bave is the average > B-factor of the protein. This will make some of BS negative (if B
Re: [ccp4bb] B-factor standardization
Dear Oliviero, On one aspect of your query, this analysis dissected the different sources of disorder contributions to B factors:- Diamond, R Acta Cryst (1990). *A46*, 425-435 All best wishes, John On Thu, Apr 5, 2018 at 2:49 PM, Oliviero Carugo < oliviero.car...@univie.ac.at> wrote: > Dears, > > everybody knows that B-factors may change amongst different crystal > structures and that they need to be standardized when different protein > crystal structures are compared. > > If I am not wrong, I remember that someone proposed to standardize > B-factors of protein atoms as “BS = B - Bave”, where Bave is the average > B-factor of the protein. Such standardization is based on the hypothesis > that independent sources of disorder add in determining the final B-factor. > BS should represent atomic B-factors depurated by all factors different > from atom oscillation, since Bave differences are neutralized. > > Does anyone can help me in finding a publication (80s or 90s) where these > BS values are used? > > Thanks! > > Oliviero > -- Professor John R Helliwell DSc
Re: [ccp4bb] B-factor standardization
On Thursday, 05 April 2018 15:49:44 Oliviero Carugo wrote: > Dears, > > everybody knows that B-factors may change amongst different crystal > structures and that they need to be standardized when different protein > crystal structures are compared. > > If I am not wrong, I remember that someone proposed to standardize > B-factors of protein atoms as “BS = B - Bave”, where Bave is the average > B-factor of the protein. Such standardization is based on the hypothesis > that independent sources of disorder add in determining the final > B-factor. BS should represent atomic B-factors depurated by all factors > different from atom oscillation, since Bave differences are neutralized. That sounds like flawed logic to me. If you were going to do this at all, I think you would have to start by comparing the residual B factors after removing TLS contributions. > Does anyone can help me in finding a publication (80s or 90s) where > these BS values are used? I think there are better tools available now than there were 25 years ago. I wouldn't go back that far to choose one without first considering the work and thought that has been put into structure analysis in the meantime. If nothing else, you might consider that "Bave" as calculated from the "Biso" (really Beq) column in typical PDB files is a less than perfect approximation [Acta Cryst. 2012 A67:512-516] Ethan > Thanks! > > Oliviero -- Ethan A Merritt, Dept of Biochemistry Biomolecular Structure Center, K-428 Health Sciences Bldg MS 357742, University of Washington, Seattle 98195-7742
Re: [ccp4bb] B-factor standardization
Oliviero, We published a paper in 2003 in which we normalized B-factors to do a comparison of relative mobility or flexibility. The reference is: Gourinath et al. Structure 11:1621-1627 (2003). The jargon we used for Bave of the protein is . In our case, to be conservative in our interpretation of B-factor comparisons, we excluded random turns from the calculation and only included alpha helices and beta sheets, because we were trying to get a handle on the degree of excess flexibility of the less ordered regions of a large structure. -Daniel On Thu, Apr 5, 2018 at 9:49 AM, Oliviero Carugo < oliviero.car...@univie.ac.at> wrote: > Dears, > > everybody knows that B-factors may change amongst different crystal > structures and that they need to be standardized when different protein > crystal structures are compared. > > If I am not wrong, I remember that someone proposed to standardize > B-factors of protein atoms as “BS = B - Bave”, where Bave is the average > B-factor of the protein. Such standardization is based on the hypothesis > that independent sources of disorder add in determining the final B-factor. > BS should represent atomic B-factors depurated by all factors different > from atom oscillation, since Bave differences are neutralized. > > Does anyone can help me in finding a publication (80s or 90s) where these > BS values are used? > > Thanks! > > Oliviero >
[ccp4bb] B-factor standardization
Dears, everybody knows that B-factors may change amongst different crystal structures and that they need to be standardized when different protein crystal structures are compared. If I am not wrong, I remember that someone proposed to standardize B-factors of protein atoms as “BS = B - Bave”, where Bave is the average B-factor of the protein. Such standardization is based on the hypothesis that independent sources of disorder add in determining the final B-factor. BS should represent atomic B-factors depurated by all factors different from atom oscillation, since Bave differences are neutralized. Does anyone can help me in finding a publication (80s or 90s) where these BS values are used? Thanks! Oliviero
Re: [ccp4bb] Fwd: [ccp4bb] B-factor blurring
Dear Mike, I am not sure what you mean by 'reference'. When you multiply your data with exp(-Bs^2) with B > 0 you weight your reflections depending on resolution so that high resolution data (probably including noise) are supressed while low resolution data a underlined. Using the term 'B' with its relationship to the ADPs is probably only a way to make it sound more fancy. You can call this textbook knowledge, maybe that's why there is not a real reference? Best, Tim On 11/19/2014 01:34 AM, Mike Lawrence wrote: > My sincere thanks to all who are responding to my request below. > > To be explicit, my question relates to B-factor blurring (+B correction), not > to B-factor sharpening (-B correction). > > thanks > > Mike > > > Begin forwarded message: > >> From: Mike Lawrence >> Subject: [ccp4bb] B-factor blurring >> Date: 18 November 2014 12:01:07 pm AEDT >> To: CCP4BB@JISCMAIL.AC.UK >> Reply-To: Mike Lawrence >> >> Dear all >> >> can anyone help me with literatures references to B-factor blurring as a >> technique to reveal low resolution features in an electron density map? I >> have seen the poster from Andrea Thorn at >> >> http://shelx.uni-ac.gwdg.de/~athorn/pdf/thorn_iucr2014_poster.pdf >> >> but was wondering if there were any alternative references? >> >> with many thanks! >> >> Mike >> >> > > > > __ > The information in this email is confidential and intended solely for the > addressee. > You must not disclose, forward, print or use it without the permission of the > sender. > __ > -- Dr Tim Gruene Institut fuer anorganische Chemie Tammannstr. 4 D-37077 Goettingen GPG Key ID = A46BEE1A signature.asc Description: OpenPGP digital signature
[ccp4bb] Fwd: [ccp4bb] B-factor blurring
My sincere thanks to all who are responding to my request below. To be explicit, my question relates to B-factor blurring (+B correction), not to B-factor sharpening (-B correction). thanks Mike Begin forwarded message: > From: Mike Lawrence > Subject: [ccp4bb] B-factor blurring > Date: 18 November 2014 12:01:07 pm AEDT > To: CCP4BB@JISCMAIL.AC.UK > Reply-To: Mike Lawrence > > Dear all > > can anyone help me with literatures references to B-factor blurring as a > technique to reveal low resolution features in an electron density map? I > have seen the poster from Andrea Thorn at > > http://shelx.uni-ac.gwdg.de/~athorn/pdf/thorn_iucr2014_poster.pdf > > but was wondering if there were any alternative references? > > with many thanks! > > Mike > > __ The information in this email is confidential and intended solely for the addressee. You must not disclose, forward, print or use it without the permission of the sender. __
Re: [ccp4bb] B-factor blurring
Another useful reference: Liu C & Xiong Y (2014). Electron Density Sharpening as a General Crystallographic Technique. *J. Mol. Biol.* 426, 980-993. http://www.ncbi.nlm.nih.gov/pubmed/24269527 Ho Leung Ng University of Hawaii at Manoa Assistant Professor, Department of Chemistry h...@hawaii.edu
Re: [ccp4bb] B-factor blurring
Mike, you can have a look here (from page 52 and on): http://phenix-online.org/presentations/fem_06MAY2014.pdf There it is shown that map kurtosis can be used as an optimization criterion for choosing the best sharpening B. Pavel On Mon, Nov 17, 2014 at 5:01 PM, Mike Lawrence wrote: > Dear all > > can anyone help me with literatures references to B-factor blurring as a > technique to reveal low resolution features in an electron density map? I > have seen the poster from Andrea Thorn at > > http://shelx.uni-ac.gwdg.de/~athorn/pdf/thorn_iucr2014_poster.pdf > > but was wondering if there were any alternative references? > > with many thanks! > > Mike > > > > > Mike Lawrence, PhD > > Associate Professor and WEHI Fellow > Structural Biology Division > Walter and Eliza Hall Institute of Medical Research > 1G Royal Parade, Parkville > Victoria 3052, AUSTRALIA > > Tel. 61-3-9345-2693 > Fax 61-3-9345-2686 > Email: lawre...@wehi.edu.au > > > > > __ > The information in this email is confidential and intended solely for the > addressee. > You must not disclose, forward, print or use it without the permission of > the sender. > __ >
Re: [ccp4bb] B-factor blurring
Hi Mike, my favourite summary on B factor sharpening/blurring is here: Acta Crystallogr D Biol Crystallogr. 2006 Aug;62(Pt 8):923-32. Epub 2006 Jul 18. Considerations for the refinement of low-resolution crystal structures. DeLaBarre B, Brunger AT. Hope that helps, Tomas On Mon, Nov 17, 2014 at 8:01 PM, Mike Lawrence wrote: > Dear all > > can anyone help me with literatures references to B-factor blurring as a > technique to reveal low resolution features in an electron density map? I > have seen the poster from Andrea Thorn at > > http://shelx.uni-ac.gwdg.de/~athorn/pdf/thorn_iucr2014_poster.pdf > > but was wondering if there were any alternative references? > > with many thanks! > > Mike > > > > > Mike Lawrence, PhD > > Associate Professor and WEHI Fellow > Structural Biology Division > Walter and Eliza Hall Institute of Medical Research > 1G Royal Parade, Parkville > Victoria 3052, AUSTRALIA > > Tel. 61-3-9345-2693 > Fax 61-3-9345-2686 > Email: lawre...@wehi.edu.au > > > > > __ > The information in this email is confidential and intended solely for the > addressee. > You must not disclose, forward, print or use it without the permission of > the sender. > __ >
[ccp4bb] B-factor blurring
Dear all can anyone help me with literatures references to B-factor blurring as a technique to reveal low resolution features in an electron density map? I have seen the poster from Andrea Thorn at http://shelx.uni-ac.gwdg.de/~athorn/pdf/thorn_iucr2014_poster.pdf but was wondering if there were any alternative references? with many thanks! Mike Mike Lawrence, PhD Associate Professor and WEHI Fellow Structural Biology Division Walter and Eliza Hall Institute of Medical Research 1G Royal Parade, Parkville Victoria 3052, AUSTRALIA Tel. 61-3-9345-2693 Fax 61-3-9345-2686 Email: lawre...@wehi.edu.au __ The information in this email is confidential and intended solely for the addressee. You must not disclose, forward, print or use it without the permission of the sender. __
[ccp4bb] B-factor statistics around an ion - Summary
Dear all, I received a few replies to my request, all but Pavel's privately. Those included suggestions of using moleman2, pymol, baverage (from ccp4). Several required manual intervention like visual selection of the responding atoms, which made them less appealing to me. I used Robbie Joosten's suggestion who pointed out that the PDB updated all PDB entries to contain link records and that these link records match the atom description on the ATOM/HETATM cards. As my perl script had a bug I could not resolve within an acceptable period of time I wrote a c++ program which makes use of the STL map to create per-Atom statistics of the atoms within the link record. I attach the program in case anyone is interested. It is not debugged and utterly slow (21s for 50PDB files), but it served my purpose. Thanks to everyone who responded. Best, Tim On 01/23/2014 01:36 PM, Tim Gruene wrote: > Dear all, > > could anyone suggest a way to get the average B-factor from a PDB-file > of those atoms a specific ion binds to (e.g. as judged by header LINK > records or a distance interval)? > > I would like to get this number for all K-ions from a set of > PDB-files, and I hope there is a quicker way than using coot to click > on the surrounding atoms and transfer the numbers into a chart. > > Best, > Tim > > -- Dr Tim Gruene Institut fuer anorganische Chemie Tammannstr. 4 D-37077 Goettingen GPG Key ID = A46BEE1A #include #include #include #include #include class KData { public: double BK_; double enviB_; int numenvi_; KData() : BK_(0.0), enviB_(0.0), numenvi_(0) { } }; void printKs(const std::map > & Ks) { std::map >::const_iterator it; for (it = Ks.begin(); it != Ks.end(); ++it) { std::cout << "K: " << it->first << ":\n"; std::vector::const_iterator it2; for (it2 = it->second.begin(); it2 != it->second.end(); ++it2) { std::cout << " ---> " << *it2 << "\n"; } } } void printData (const std::map& data) { std::map::const_iterator it; for (it = data.begin(); it != data.end(); ++it) { std::cout << "# ---> \"" << it->first << "\": " << " number of LINKed atoms: " << it->second.numenvi_ << ";" << " B/ average B:\n" << it->second.BK_ << " " << it->second.enviB_ / it->second.numenvi_<< "\n"; } } void checkline(const std::string& line, const std::map > & Ks, std::map& data) { std::map >::const_iterator it; for (it = Ks.begin(); it != Ks.end(); ++it) { if (line.find(it->first) != std::string::npos) { std::string B(line.substr(60, 6)); double b(std::stod(B)); data[it->first].BK_ = b; } std::vector::const_iterator it2; for (it2 = it->second.begin(); it2 != it->second.end(); ++it2) { if (line.find(*it2) != std::string::npos) { std::string B(line.substr(60, 6)); double b(std::stod(B)); data[it->first].enviB_ += b; ++(data[it->first].numenvi_); } } } } int main(int argc, char* argv[]) { std::map > Kenvi; std::map data; if (argc < 2) { std::cout << "*** Error: please provide one of more PDB filenames.\n"; } char* line; int linewidth(100); line = new char[linewidth + 1]; for (int i = 1; i < argc; ++i) { std::ifstream inp(argv[i]); if (!inp.is_open()) { std::cout << "---> Cannot open file " << argv[i] << "\n"; } // extract link records while (true) { inp.getline(line, linewidth); if (!inp.good()) break; std::string myline(line); // break on CRYST1 if (myline.substr(0, 6) == "CRYST1") { break; } if (myline.substr(0, 4) == "LINK") { std::string atom1 = myline.substr(12, 14); std::string atom2 = myline.substr(42, 14); if (atom1.find(" K") != std::string::npos) { Kenvi[atom1].push_back(atom2); } else if (atom2.find(" K") != std::string::npos) { Kenvi[atom2].push_back(atom1); } } } // map is set, find atoms while (true) { inp.getline(line, linewidth); if (!inp.good()) break; std::string myline(line); if (myline.substr(0, 6) == "ATOM " || myline.substr (0, 6) == "HETATM") { // check each line against map content to extract K-values checkline (myline, Kenvi, data); } } } delete [] line; #ifdef DEBUG printKs(Kenvi); #endif printData(data); return 0; } signature.asc Description: OpenPGP digital signature
Re: [ccp4bb] B-factor statistics around an ion
Hi Tim, a non-ccp4 solution (since you haven't gotten any suggestion yet).. 1) Get atom selection of atoms involved into ion coordination: phenix.metal_coordination model.pdb 2) Using atom selection from above extract portion of PDB that contains atoms in question: phenix.pdb_atom_selection model.pdb "atom_selection_string_here" --write-pdb-file="partial.pdb" 3) Get average B: phenix.pdbtools partial.pdb model_statistics=true Remark: steps 2)-3) can be done in one go: phenix.b_factor_statistics model.pdb selection="atom_selection_string_here" If this is of general interest I can wrap 1)-3) into one command. Good luck, Pavel On Thu, Jan 23, 2014 at 4:36 AM, Tim Gruene wrote: > -BEGIN PGP SIGNED MESSAGE- > Hash: SHA1 > > Dear all, > > could anyone suggest a way to get the average B-factor from a PDB-file > of those atoms a specific ion binds to (e.g. as judged by header LINK > records or a distance interval)? > > I would like to get this number for all K-ions from a set of > PDB-files, and I hope there is a quicker way than using coot to click > on the surrounding atoms and transfer the numbers into a chart. > > Best, > Tim > > - -- > - -- > Dr Tim Gruene > Institut fuer anorganische Chemie > Tammannstr. 4 > D-37077 Goettingen > > GPG Key ID = A46BEE1A > > -BEGIN PGP SIGNATURE- > Version: GnuPG v1.4.12 (GNU/Linux) > Comment: Using GnuPG with Icedove - http://www.enigmail.net/ > > iD8DBQFS4QziUxlJ7aRr7hoRAtlAAKCei/Ehfg71Xir3pkvWvTQg93YwEQCePFM6 > 2kPFPJw8aZySVYnLaQCjKR4= > =q9pT > -END PGP SIGNATURE- >
[ccp4bb] B-factor statistics around an ion
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Dear all, could anyone suggest a way to get the average B-factor from a PDB-file of those atoms a specific ion binds to (e.g. as judged by header LINK records or a distance interval)? I would like to get this number for all K-ions from a set of PDB-files, and I hope there is a quicker way than using coot to click on the surrounding atoms and transfer the numbers into a chart. Best, Tim - -- - -- Dr Tim Gruene Institut fuer anorganische Chemie Tammannstr. 4 D-37077 Goettingen GPG Key ID = A46BEE1A -BEGIN PGP SIGNATURE- Version: GnuPG v1.4.12 (GNU/Linux) Comment: Using GnuPG with Icedove - http://www.enigmail.net/ iD8DBQFS4QziUxlJ7aRr7hoRAtlAAKCei/Ehfg71Xir3pkvWvTQg93YwEQCePFM6 2kPFPJw8aZySVYnLaQCjKR4= =q9pT -END PGP SIGNATURE-
Re: [ccp4bb] B factor of the loop
Indeed! If the B factors are rather large compared to the globular protein core (assuming there is a globular core being that the protein crystallized), one can make the assumption, especially within a loop region, that this is an indirect measurement of flexibility. However, as Jürgen pointed out, it IS imperative to take a close look at the crystal packing in the unit cell. For instance, if the loop region were to make hydrogen bond or electrostatic interactions with a symmetry mate, you must be careful in your conclusion. Might I recommend a paper that uses B factors as a direct correlation with heteronuclear NOEs to compare two almost identical structures (both of which contain disordered regions? Wang, Y., Fisher, J.C., Assem, M., Matthew, R., Sublet, J, Xiao, L., Roussel, M.F., and Kriwacki, R.W. (2011). Structural basis for the diverse cell cycle regulatory functions of the intrinsically disordered protein, p21 Cip1. Nature Chemical Biology, 7, 214-221. You can always confirm disordered or flexible loop segments using limited proteolysis. Best, John John Fisher, M.D./PhD St. Jude Children's Research Hospital Department of Structural Biology Department of Oncology On Sun, Mar 3, 2013 at 3:52 PM, Bosch, Juergen wrote: > yes - but keep in mind your protein is in the context of the crystal > lattice, so flexible regions in solution are likely to be stabilized in the > crystal lattice. So if you color by B also look at the symmetry mates. > And you should also submit both structures to the TLSMD server and look at > those results. > http://skuld.bmsc.washington.edu/~tlsmd/ > > Jürgen > > On Mar 3, 2013, at 4:35 PM, Faisal Tarique wrote: > > > Dear all > > > > Can B factor in the crystal structure be the criteria to look into the > flexibility of a region or domain.? Also if two structures are at > different resolutions. > > > > Faisal > > -- > > > > .. > Jürgen Bosch > Johns Hopkins University > Bloomberg School of Public Health > Department of Biochemistry & Molecular Biology > Johns Hopkins Malaria Research Institute > 615 North Wolfe Street, W8708 > Baltimore, MD 21205 > Office: +1-410-614-4742 > Lab: +1-410-614-4894 > Fax: +1-410-955-2926 > http://lupo.jhsph.edu >
Re: [ccp4bb] B factor of the loop
yes - but keep in mind your protein is in the context of the crystal lattice, so flexible regions in solution are likely to be stabilized in the crystal lattice. So if you color by B also look at the symmetry mates. And you should also submit both structures to the TLSMD server and look at those results. http://skuld.bmsc.washington.edu/~tlsmd/ Jürgen On Mar 3, 2013, at 4:35 PM, Faisal Tarique wrote: > Dear all > > Can B factor in the crystal structure be the criteria to look into the > flexibility of a region or domain.? Also if two structures are at different > resolutions. > > Faisal > -- > .. Jürgen Bosch Johns Hopkins University Bloomberg School of Public Health Department of Biochemistry & Molecular Biology Johns Hopkins Malaria Research Institute 615 North Wolfe Street, W8708 Baltimore, MD 21205 Office: +1-410-614-4742 Lab: +1-410-614-4894 Fax: +1-410-955-2926 http://lupo.jhsph.edu
[ccp4bb] B factor of the loop
Dear all Can B factor in the crystal structure be the criteria to look into the flexibility of a region or domain.? Also if two structures are at different resolutions. Faisal --
Re: [ccp4bb] B factor
Hi Amit, you can manually define the range of a spectrum in PyMol: spectrum [parameter], [spectrum type], minimum=[min], maximum=[max] --> if you want to color according to B-factor this will translate to e.g. : spectrum b, blue_white_red, minimum=0, maximum=50 hth, Arjen On May 12, 2011, at 8:31 PM, Luthra,Amit wrote: Hey all, I would like to make one ribbon diagram of the structure showing the B factor in different colors (contour levels). I am using pymol but it is not representing the perfect contour level of B factor. Is any other program where I can change the color according to B factor? Thanks Amit Luthra, Ph.D. Post-Doctoral Fellow The Radolf Laboratory Department of Medicine University of Connecticut Health Center [w] http://spirocheteresearch.uchc.edu
Re: [ccp4bb] B factor
Hi Amit, you can manually define the range of a spectrum in PyMol: spectrum [parameter], [spectrum type], minimum=[min], maximum=[max] --> if you want to color according to B-factor this will translate to e.g. : spectrum b, blue_white_red, minimum=0, maximum=50 hth, Arjen On May 12, 2011, at 8:31 PM, Luthra,Amit wrote: Hey all, I would like to make one ribbon diagram of the structure showing the B factor in different colors (contour levels). I am using pymol but it is not representing the perfect contour level of B factor. Is any other program where I can change the color according to B factor? Thanks Amit Luthra, Ph.D. Post-Doctoral Fellow The Radolf Laboratory Department of Medicine University of Connecticut Health Center [w] http://spirocheteresearch.uchc.edu
Re: [ccp4bb] B factor
hello sir , in chimera its available in module Render by Attribute. go through this link http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/render/render.html . regards On Fri, May 13, 2011 at 12:01 AM, Luthra,Amit wrote: > Hey all, > > I would like to make one ribbon diagram of the structure showing the B > factor in different colors (contour levels). > I am using pymol but it is not representing the perfect contour level of B > factor. Is any other program where I can change the color according to B > factor? > Thanks > > > Amit Luthra, Ph.D. > Post-Doctoral Fellow > The Radolf Laboratory > Department of Medicine > University of Connecticut Health Center > [w] http://spirocheteresearch.uchc.edu > -- Vandana Kukshal Postdoc Fellow structure biology group international Center For Genetic Engineering and Biotechnolgy India
[ccp4bb] B factor
Hey all, I would like to make one ribbon diagram of the structure showing the B factor in different colors (contour levels). I am using pymol but it is not representing the perfect contour level of B factor. Is any other program where I can change the color according to B factor? Thanks Amit Luthra, Ph.D. Post-Doctoral Fellow The Radolf Laboratory Department of Medicine University of Connecticut Health Center [w] http://spirocheteresearch.uchc.edu
[ccp4bb] B-factor-weighted Superposition
Dear Crystallographers, do any of the commonly-used superposing algorithms weight superpositions by b-factors? Seems like it would be a good idea, so I assume it is implemented somewhere? Jacob Keller -- *** Jacob Pearson Keller Northwestern University Medical Scientist Training Program cel: 773.608.9185 email: j-kell...@northwestern.edu ***
Re: [ccp4bb] B factor Coloring in Pymol
Hi All Seems Roberts suggestion has worked. Thank you Robert Best Gina -Original Message- From: CCP4 bulletin board [mailto:ccp...@jiscmail.ac.uk] On Behalf Of Clayton, Gina Martyn Sent: Wednesday, March 31, 2010 4:34 PM To: CCP4BB@JISCMAIL.AC.UK Subject: Re: [ccp4bb] B factor Coloring in Pymol Hi Robert That seems like a great suggestion! I will give it a shot and let CCP4BB know if it works Thanks Gina -Original Message- From: CCP4 bulletin board [mailto:ccp...@jiscmail.ac.uk] On Behalf Of Robert Campbell Sent: Wednesday, March 31, 2010 4:29 PM To: CCP4BB@JISCMAIL.AC.UK Subject: Re: [ccp4bb] B factor Coloring in Pymol Hi Gina, On Wed, 31 Mar 2010 15:00:34 -0400 "Clayton, Gina Martyn" wrote: > I am trying to colour a structure by B factor in Pymol. More > specifically I am colouring conserved residues (value in b factor). I > want to use 4 colours - say white, yellow, orange, red. However it > seems that Pymol now only allows 3 colours to be used - I tried the > script from the Wiki site, but Pymol will not accept the 4 rgb colours. > Does anyone have an idea how I can colour with 4 colours (but not > rainbow too confusing). Say, your b-factor values range from 0-10, and your object is called "protein", then you could do: color white, protein color yellow, protein & b > 2.5 color orange, protein & b > 5.0 color red, protein & b > 7.5 If you want to automate that process, then you need to write a script that determines, the b-factor limits for the 4 colours based on the b-factors present in your structure. Cheers, Rob -- Robert L. Campbell, Ph.D. Senior Research Associate/Adjunct Assistant Professor Botterell Hall Rm 644 Department of Biochemistry, Queen's University, Kingston, ON K7L 3N6 Canada Tel: 613-533-6821Fax: 613-533-2497 http://pldserver1.biochem.queensu.ca/~rlc Notice: This e-mail message, together with any attachments, contains information of Merck & Co., Inc. (One Merck Drive, Whitehouse Station, New Jersey, USA 08889), and/or its affiliates Direct contact information for affiliates is available at http://www.merck.com/contact/contacts.html) that may be confidential, proprietary copyrighted and/or legally privileged. It is intended solely for the use of the individual or entity named on this message. If you are not the intended recipient, and have received this message in error, please notify us immediately by reply e-mail and then delete it from your system. Notice: This e-mail message, together with any attachments, contains information of Merck & Co., Inc. (One Merck Drive, Whitehouse Station, New Jersey, USA 08889), and/or its affiliates Direct contact information for affiliates is available at http://www.merck.com/contact/contacts.html) that may be confidential, proprietary copyrighted and/or legally privileged. It is intended solely for the use of the individual or entity named on this message. If you are not the intended recipient, and have received this message in error, please notify us immediately by reply e-mail and then delete it from your system.
Re: [ccp4bb] B factor Coloring in Pymol
Hi Robert That seems like a great suggestion! I will give it a shot and let CCP4BB know if it works Thanks Gina -Original Message- From: CCP4 bulletin board [mailto:ccp...@jiscmail.ac.uk] On Behalf Of Robert Campbell Sent: Wednesday, March 31, 2010 4:29 PM To: CCP4BB@JISCMAIL.AC.UK Subject: Re: [ccp4bb] B factor Coloring in Pymol Hi Gina, On Wed, 31 Mar 2010 15:00:34 -0400 "Clayton, Gina Martyn" wrote: > I am trying to colour a structure by B factor in Pymol. More > specifically I am colouring conserved residues (value in b factor). I > want to use 4 colours - say white, yellow, orange, red. However it > seems that Pymol now only allows 3 colours to be used - I tried the > script from the Wiki site, but Pymol will not accept the 4 rgb colours. > Does anyone have an idea how I can colour with 4 colours (but not > rainbow too confusing). Say, your b-factor values range from 0-10, and your object is called "protein", then you could do: color white, protein color yellow, protein & b > 2.5 color orange, protein & b > 5.0 color red, protein & b > 7.5 If you want to automate that process, then you need to write a script that determines, the b-factor limits for the 4 colours based on the b-factors present in your structure. Cheers, Rob -- Robert L. Campbell, Ph.D. Senior Research Associate/Adjunct Assistant Professor Botterell Hall Rm 644 Department of Biochemistry, Queen's University, Kingston, ON K7L 3N6 Canada Tel: 613-533-6821Fax: 613-533-2497 http://pldserver1.biochem.queensu.ca/~rlc Notice: This e-mail message, together with any attachments, contains information of Merck & Co., Inc. (One Merck Drive, Whitehouse Station, New Jersey, USA 08889), and/or its affiliates Direct contact information for affiliates is available at http://www.merck.com/contact/contacts.html) that may be confidential, proprietary copyrighted and/or legally privileged. It is intended solely for the use of the individual or entity named on this message. If you are not the intended recipient, and have received this message in error, please notify us immediately by reply e-mail and then delete it from your system.
Re: [ccp4bb] B factor Coloring in Pymol
Hi Ron Mmmm... not really but thanks for the info:) I have my Bfactor column assigned in percent conservation so the B value ranges from say 30% to 100%. I want to colour my molecule accordingly. Pymol allows B factor colouring but seems only to allow 3 colour choices. It used to allow 4 colour choices say "WYOR" I was wondering if any one had an idea how I could trick Pymol into accepting say 4 colours. Sorry if I was not clearer in the previous. Thanks again advance. Gina -Original Message- From: CCP4 bulletin board [mailto:ccp...@jiscmail.ac.uk] On Behalf Of Rob Nicholls Sent: Wednesday, March 31, 2010 3:14 PM To: CCP4BB@JISCMAIL.AC.UK Subject: Re: [ccp4bb] B factor Coloring in Pymol Hi, I'm aware of two ways to do this: by writing the colour as a word, e.g: color white, (mypdbfile//A/*/*) which will colour a whole chain. Or by using rgb colours, e.g: set_color newcolor0 = [1,1,1]; color newcolor0, (mypdbfile//A/2/*) which will colour individual residues (residue 2 in this case) Pymol should accept a script containing such lines. Hope that helps, Rob On 31 Mar 2010, at 15:00, Clayton, Gina Martyn wrote: > > Dear Everyone > > Slightly off topic... > > I am trying to colour a structure by B factor in Pymol. More > specifically I am colouring conserved residues (value in b factor). I > want to use 4 colours - say white, yellow, orange, red. However it > seems that Pymol now only allows 3 colours to be used - I tried the > script from the Wiki site, but Pymol will not accept the 4 rgb > colours. > Does anyone have an idea how I can colour with 4 colours (but not > rainbow too confusing). > > > Thanks in advance for any help > > Best > Gina > > Notice: This e-mail message, together with any attachments, > contains information of Merck & Co., Inc. (One Merck Drive, > Whitehouse Station, New Jersey, USA 08889), and/or its affiliates > Direct contact information for affiliates is available at http://www.merck.com/contact/contacts.html) > that may be confidential, proprietary copyrighted and/or legally > privileged. It is intended solely for the use of the individual or > entity named on this message. If you are not the intended recipient, > and have received this message in error, please notify us > immediately by reply e-mail and then delete it from your system. Notice: This e-mail message, together with any attachments, contains information of Merck & Co., Inc. (One Merck Drive, Whitehouse Station, New Jersey, USA 08889), and/or its affiliates Direct contact information for affiliates is available at http://www.merck.com/contact/contacts.html) that may be confidential, proprietary copyrighted and/or legally privileged. It is intended solely for the use of the individual or entity named on this message. If you are not the intended recipient, and have received this message in error, please notify us immediately by reply e-mail and then delete it from your system.
Re: [ccp4bb] B factor Coloring in Pymol
Hi Gina, On Wed, 31 Mar 2010 15:00:34 -0400 "Clayton, Gina Martyn" wrote: > I am trying to colour a structure by B factor in Pymol. More > specifically I am colouring conserved residues (value in b factor). I > want to use 4 colours - say white, yellow, orange, red. However it > seems that Pymol now only allows 3 colours to be used - I tried the > script from the Wiki site, but Pymol will not accept the 4 rgb colours. > Does anyone have an idea how I can colour with 4 colours (but not > rainbow too confusing). Say, your b-factor values range from 0-10, and your object is called "protein", then you could do: color white, protein color yellow, protein & b > 2.5 color orange, protein & b > 5.0 color red, protein & b > 7.5 If you want to automate that process, then you need to write a script that determines, the b-factor limits for the 4 colours based on the b-factors present in your structure. Cheers, Rob -- Robert L. Campbell, Ph.D. Senior Research Associate/Adjunct Assistant Professor Botterell Hall Rm 644 Department of Biochemistry, Queen's University, Kingston, ON K7L 3N6 Canada Tel: 613-533-6821Fax: 613-533-2497 http://pldserver1.biochem.queensu.ca/~rlc
Re: [ccp4bb] B factor Coloring in Pymol
Hi, I'm aware of two ways to do this: by writing the colour as a word, e.g: color white, (mypdbfile//A/*/*) which will colour a whole chain. Or by using rgb colours, e.g: set_color newcolor0 = [1,1,1]; color newcolor0, (mypdbfile//A/2/*) which will colour individual residues (residue 2 in this case) Pymol should accept a script containing such lines. Hope that helps, Rob On 31 Mar 2010, at 15:00, Clayton, Gina Martyn wrote: Dear Everyone Slightly off topic... I am trying to colour a structure by B factor in Pymol. More specifically I am colouring conserved residues (value in b factor). I want to use 4 colours - say white, yellow, orange, red. However it seems that Pymol now only allows 3 colours to be used - I tried the script from the Wiki site, but Pymol will not accept the 4 rgb colours. Does anyone have an idea how I can colour with 4 colours (but not rainbow too confusing). Thanks in advance for any help Best Gina Notice: This e-mail message, together with any attachments, contains information of Merck & Co., Inc. (One Merck Drive, Whitehouse Station, New Jersey, USA 08889), and/or its affiliates Direct contact information for affiliates is available at http://www.merck.com/contact/contacts.html) that may be confidential, proprietary copyrighted and/or legally privileged. It is intended solely for the use of the individual or entity named on this message. If you are not the intended recipient, and have received this message in error, please notify us immediately by reply e-mail and then delete it from your system.
[ccp4bb] B factor Coloring in Pymol
Dear Everyone Slightly off topic... I am trying to colour a structure by B factor in Pymol. More specifically I am colouring conserved residues (value in b factor). I want to use 4 colours - say white, yellow, orange, red. However it seems that Pymol now only allows 3 colours to be used - I tried the script from the Wiki site, but Pymol will not accept the 4 rgb colours. Does anyone have an idea how I can colour with 4 colours (but not rainbow too confusing). Thanks in advance for any help Best Gina Notice: This e-mail message, together with any attachments, contains information of Merck & Co., Inc. (One Merck Drive, Whitehouse Station, New Jersey, USA 08889), and/or its affiliates Direct contact information for affiliates is available at http://www.merck.com/contact/contacts.html) that may be confidential, proprietary copyrighted and/or legally privileged. It is intended solely for the use of the individual or entity named on this message. If you are not the intended recipient, and have received this message in error, please notify us immediately by reply e-mail and then delete it from your system.
Re: [ccp4bb] B-factor problem
Sampath Natarajan wrote: Dear All, I am refining a structure with 2.5A resolution by refmac5. I could find the solution by MR using molrep. After fitting the model, I refined the structure again with 0.3 weighting term, but the output PDB file shows many splits in the residues. So I used 'auto' as a weighting term, then I didn't find any splits. At the same time I noticed that figure of merits is not increasing after model fitting. But the correlation coefficient is more than 80%. Also while refining with auto weighting term, the figure of merit goes down about 25% than previous. Finally I find that some of the residues show high B factor. Could anyone give suggestions to reduce the B-factor? Thanks in advance! Yours sincerely, Sampath There are several Qs here but you must remember that initial refinement protocols straight from an MR solution are somewhat different to those you would use later. High B factors may be useful - they often indicate where something is wrong, and ditto poor geometry. These are the regions you will probably need to rebuild. Eleanor , to the
Re: [ccp4bb] B-factor problem
Hi, What is the lower limit of the weighting term one can use. If the data is around 2A can one use 0.04 or less and which type of refinement is more useful when one is doing the initial refinement (isotropic, aisotropic, overall or mixed). sincerely Debajyoti On Tue, 17 Jun 2008 Roger Rowlett wrote : >Sampath Natarajan wrote: >>Dear All, >>I am refining a structure with 2.5A resolution by refmac5. I >> could find the solution by MR using molrep. After fitting the model, I >> refined the structure again with 0.3 weighting term, but the output PDB file >> shows many splits in the residues. So I used 'auto' as a weighting term, >> then I didn't find any splits. At the same time I noticed that figure of >> merits is not increasing after model fitting. But the correlation >> coefficient is more than 80%. Also while refining with auto weighting term, >> the figure of merit goes down about 25% than previous. >> >> Finally I find that some of the residues show high B factor. Could >> anyone give suggestions to reduce the B-factor? >>Thanks in advance! >> Yours sincerely, >>Sampath >> >For 2.5A data, the refmac default X-ray weighting factor of 0.3 is typically >much too high, and does not appropriately apply geometric constraints. (That >is why your proteins chains are breaking up.) An alternative to the auto >weighting function is to set your own weighting factor and monitor rmsBOND and >rmsANGLE deviations in the final model. If the final rmsBOND is on the order >of 0.015 A and rmsBOND is on the order of 1.5 degrees then your choice of >weighting factor is reasonable. For medium resolution data in the 2.2-2.8A >range, I would start with an X-ray weighting factor of about 0.05 or so. > >Cheers, > > >-- >Roger S. Rowlett >Professor >Colgate University Presidential Scholar >Department of Chemistry >Colgate University >13 Oak Drive >Hamilton, NY 13346 > >tel: (315)-228-7245 >ofc: (315)-228-7395 >fax: (315)-228-7935 >email: [EMAIL PROTECTED]
Re: [ccp4bb] B-factor problem
On Tuesday 17 June 2008 09:08:24 am Sampath Natarajan wrote: > Dear All, >I am refining a structure with 2.5A resolution by refmac5. I > could find the solution by MR using molrep. After fitting the model, I > refined the structure again with 0.3 weighting term, but the output PDB file > shows many splits in the residues. So I used 'auto' as a weighting term, > then I didn't find any splits. At the same time I noticed that figure of > merits is not increasing after model fitting. But the correlation > coefficient is more than 80%. Also while refining with auto weighting term, > the figure of merit goes down about 25% than previous. > > Finally I find that some of the residues show high B factor. Could > anyone give suggestions to reduce the B-factor? It is worth remembering that B is a measure of the uncertainty in the atomic position. So one generic answer is that to reduce B you need to be more certain about where your atoms are :-) If/when your model correctly describes the mean position for each atom, then B can be interpreting as representing 'motion' or more properly speaking displacement about that mean position. But if your atoms are not in the right place, then their high B factors are simply telling you that the calculated X-ray density doesn't support your current model. So following automated MR, high B factors are best viewed as a hint as to where your model needs to be rebuilt.
Re: [ccp4bb] B-factor problem
On 17 Jun 2008, at 19:56, Roger Rowlett wrote: Sampath Natarajan wrote: Dear All, I am refining a structure with 2.5A resolution by refmac5. I could find the solution by MR using molrep. After fitting the model, I refined the structure again with 0.3 weighting term, but the output PDB file shows many splits in the residues. So I used 'auto' as a weighting term, then I didn't find any splits. At the same time I noticed that figure of merits is not increasing after model fitting. But the correlation coefficient is more than 80%. Also while refining with auto weighting term, the figure of merit goes down about 25% than previous. Finally I find that some of the residues show high B factor. Could anyone give suggestions to reduce the B-factor? Thanks in advance! Yours sincerely, Sampath For 2.5A data, the refmac default X-ray weighting factor of 0.3 is typically much too high, and does not appropriately apply geometric constraints. (That is why your proteins chains are breaking up.) ***Restraints*** really. Constraints are something really different (for one things they imply that they cannot be violated but must be exactly reproduced, so a weighting term for a constraint does not exists). I think we must avoid confusion, since implementation with bond/angle constraints (torsional refinement) are available. An alternative to the auto weighting function is to set your own weighting factor and monitor rmsBOND and rmsANGLE deviations in the final model. If the final rmsBOND is on the order of 0.015 A and rmsBOND is on the order of 1.5 degrees then your choice of weighting factor is reasonable. For medium resolution data in the 2.2-2.8A range, I would start with an X-ray weighting factor of about 0.05 or so. Although that the general advice above is of course valid and it reflects common practice, it would be fare to iterate once more that this remained for a while a somewhat controversial issue. In my opinion it is getting settled lately, given the RMSZ in more recent versions of refmac which might be a better quantity to follow - as suggested in the papers by Ian Tickle, eg: Acta Crystallogr D Biol Crystallogr. 2007 Dec;63(Pt 12):1274-81; author reply 1282-3. Epub 2007 Nov 16. Experimental determination of optimal root-mean-square deviations of macromolecular bond lengths and angles from their restrained ideal values. Tickle IJ. I claim no expertise and the math is above my head mostly, but I got pretty much convinced by that work for what I want to be doing and its good reading even for the mathematically challenged, like myself. Tassos Cheers, -- -- -- Roger S. Rowlett Professor Colgate University Presidential Scholar Department of Chemistry Colgate University 13 Oak Drive Hamilton, NY 13346 tel: (315)-228-7245 ofc: (315)-228-7395 fax: (315)-228-7935 email: [EMAIL PROTECTED]
Re: [ccp4bb] B-factor problem
Sampath Natarajan wrote: Dear All, I am refining a structure with 2.5A resolution by refmac5. I could find the solution by MR using molrep. After fitting the model, I refined the structure again with 0.3 weighting term, but the output PDB file shows many splits in the residues. So I used 'auto' as a weighting term, then I didn't find any splits. At the same time I noticed that figure of merits is not increasing after model fitting. But the correlation coefficient is more than 80%. Also while refining with auto weighting term, the figure of merit goes down about 25% than previous. Finally I find that some of the residues show high B factor. Could anyone give suggestions to reduce the B-factor? Thanks in advance! Yours sincerely, Sampath For 2.5A data, the refmac default X-ray weighting factor of 0.3 is typically much too high, and does not appropriately apply geometric constraints. (That is why your proteins chains are breaking up.) An alternative to the auto weighting function is to set your own weighting factor and monitor rmsBOND and rmsANGLE deviations in the final model. If the final rmsBOND is on the order of 0.015 A and rmsBOND is on the order of 1.5 degrees then your choice of weighting factor is reasonable. For medium resolution data in the 2.2-2.8A range, I would start with an X-ray weighting factor of about 0.05 or so. Cheers, -- Roger S. Rowlett Professor Colgate University Presidential Scholar Department of Chemistry Colgate University 13 Oak Drive Hamilton, NY 13346 tel: (315)-228-7245 ofc: (315)-228-7395 fax: (315)-228-7935 email: [EMAIL PROTECTED]
[ccp4bb] B-factor problem
Dear All, I am refining a structure with 2.5A resolution by refmac5. I could find the solution by MR using molrep. After fitting the model, I refined the structure again with 0.3 weighting term, but the output PDB file shows many splits in the residues. So I used 'auto' as a weighting term, then I didn't find any splits. At the same time I noticed that figure of merits is not increasing after model fitting. But the correlation coefficient is more than 80%. Also while refining with auto weighting term, the figure of merit goes down about 25% than previous. Finally I find that some of the residues show high B factor. Could anyone give suggestions to reduce the B-factor? Thanks in advance! Yours sincerely, Sampath
Re: [ccp4bb] B factor calculation
in ccp4, "baverage" does the job. Mark Mark J. van Raaij Dpto de Bioquímica, Facultad de Farmacia Universidad de Santiago 15782 Santiago de Compostela Spain http://web.usc.es/~vanraaij/ On 4 Jun 2008, at 21:40, xu zhen wrote: Hi, everyone, I am preparing the table of data collection and refinement statics. In this table, I need to list the average B factor of protein, ligand and water seperately, can you tell me how to calculate thoes numbers? Zhen _ News, entertainment and everything you care about at Live.com. Get it now! http://www.live.com/getstarted.aspx
Re: [ccp4bb] B factor calculation
xu zhen wrote: > Hi, everyone, > > I am preparing the table of data collection and refinement statics. In this > table, I need to list the average B factor of protein, ligand and water > seperately, can you tell me how to calculate thoes numbers? > > Zhen > _ > News, entertainment and everything you care about at Live.com. Get it now! > http://www.live.com/getstarted.aspx > MOLEMAN2 can do this with the "BF ST TY" command. Cheers, -- Roger S. Rowlett Professor Colgate University Presidential Scholar Department of Chemistry Colgate University 13 Oak Drive Hamilton, NY 13346 tel: (315)-228-7245 ofc: (315)-228-7395 fax: (315)-228-7935 email: [EMAIL PROTECTED]
Re: [ccp4bb] B factor calculation
For example, this will show you complete statistics about B-factors: phenix.pbdtools your_model.pdb --show-adp-statistics this will show you complete statistics about stereochemistry: phenix.pbdtools your_model.pdb --show-geometry-statistics For more information: http://www.phenix-online.org/ or ask me directly. Cheers, Pavel. On 6/4/2008 12:40 PM, xu zhen wrote: > Hi, everyone, > > I am preparing the table of data collection and refinement statics. In this > table, I need to list the average B factor of protein, ligand and water > seperately, can you tell me how to calculate thoes numbers? > > Zhen > _ > News, entertainment and everything you care about at Live.com. Get it now! > http://www.live.com/getstarted.aspx >
[ccp4bb] B factor calculation
Hi, everyone, I am preparing the table of data collection and refinement statics. In this table, I need to list the average B factor of protein, ligand and water seperately, can you tell me how to calculate thoes numbers? Zhen _ News, entertainment and everything you care about at Live.com. Get it now! http://www.live.com/getstarted.aspx
[ccp4bb] b-factor sharpening and FFT/CAD
Hi all, I know this was most addressed by Eleanor back on June 14, 2007... and I've had no problems making such "sharpened" maps. But I'm confused as to why FFT decides to put "SCALE 2.0 -100.0" in the command file by default, and not "SCALE 1.0 -100.0"; I've been correcting this manually. Any thoughts? Thanks, Roni
Re: [ccp4bb] B-factor sharpening in CCP4 or Coot
There is a CAD option to apply a scale and B factor to any set of data, and an FFT option to apply a scale and B factor to any column in the map. Check the documentation for whether you need SCALE 1 -10.0 or SCALE 1 +10.0 ! Eleanor Jeff Lee wrote: Hi, I have a question for everyone. I wanted to up-weight my high resolution terms in my electron density map by applying a B-factor sharpening term. I have mtz files produced from DM and I usually use Coot to generate my electron density maps. I was wondering is there an easy way to apply a B-factor sharpening term in Coot? If not, is there a program in ccp4 where by I can set the B-sharp factor? Thanks Jeff
Re: [ccp4bb] B-factor sharpening in CCP4 or Coot
Senility is, in fact, setting in. It was a modified version of CAD! Sorry for the confusion. To further compound the confusion, I am not sure if my solution is what you are really looking for, although I suspect you should try it. What we did was use an overall negative B-factor to achieve a sharpening effect (Cell 119(3):393-405). Others had done this before (David Borhani comes to mind)anyway, it can be very, very useful. Steve -Original Message- From: CCP4 bulletin board [mailto:[EMAIL PROTECTED] On Behalf Of Soisson, Stephen Michael Sent: Wednesday, June 13, 2007 1:23 PM To: CCP4BB@JISCMAIL.AC.UK Subject: Re: [ccp4bb] B-factor sharpening in CCP4 or Coot Eleanor was kind enough to modify truncate years ago for me to do this - I would guess the feature is still there. Thanks again Eleanor! Steve -Original Message- From: CCP4 bulletin board [mailto:[EMAIL PROTECTED] On Behalf Of Jeff Lee Sent: Wednesday, June 13, 2007 1:19 PM To: CCP4BB@JISCMAIL.AC.UK Subject: [ccp4bb] B-factor sharpening in CCP4 or Coot Hi, I have a question for everyone. I wanted to up-weight my high resolution terms in my electron density map by applying a B-factor sharpening term. I have mtz files produced from DM and I usually use Coot to generate my electron density maps. I was wondering is there an easy way to apply a B-factor sharpening term in Coot? If not, is there a program in ccp4 where by I can set the B-sharp factor? Thanks Jeff -- Notice: This e-mail message, together with any attachments, contains information of Merck & Co., Inc. (One Merck Drive, Whitehouse Station, New Jersey, USA 08889), and/or its affiliates (which may be known outside the United States as Merck Frosst, Merck Sharp & Dohme or MSD and in Japan, as Banyu - direct contact information for affiliates is available at http://www.merck.com/contact/contacts.html) that may be confidential, proprietary copyrighted and/or legally privileged. It is intended solely for the use of the individual or entity named on this message. If you are not the intended recipient, and have received this message in error, please notify us immediately by reply e-mail and then delete it from your system. -- -- Notice: This e-mail message, together with any attachments, contains information of Merck & Co., Inc. (One Merck Drive, Whitehouse Station, New Jersey, USA 08889), and/or its affiliates (which may be known outside the United States as Merck Frosst, Merck Sharp & Dohme or MSD and in Japan, as Banyu - direct contact information for affiliates is available at http://www.merck.com/contact/contacts.html) that may be confidential, proprietary copyrighted and/or legally privileged. It is intended solely for the use of the individual or entity named on this message. If you are not the intended recipient, and have received this message in error, please notify us immediately by reply e-mail and then delete it from your system. -- -- Notice: This e-mail message, together with any attachments, contains information of Merck & Co., Inc. (One Merck Drive, Whitehouse Station, New Jersey, USA 08889), and/or its affiliates (which may be known outside the United States as Merck Frosst, Merck Sharp & Dohme or MSD and in Japan, as Banyu - direct contact information for affiliates is available at http://www.merck.com/contact/contacts.html) that may be confidential, proprietary copyrighted and/or legally privileged. It is intended solely for the use of the individual or entity named on this message. If you are not the intended recipient, and have received this message in error, please notify us immediately by reply e-mail and then delete it from your system. --
Re: [ccp4bb] B-factor sharpening in CCP4 or Coot
Eleanor was kind enough to modify truncate years ago for me to do this - I would guess the feature is still there. Thanks again Eleanor! Steve -Original Message- From: CCP4 bulletin board [mailto:[EMAIL PROTECTED] On Behalf Of Jeff Lee Sent: Wednesday, June 13, 2007 1:19 PM To: CCP4BB@JISCMAIL.AC.UK Subject: [ccp4bb] B-factor sharpening in CCP4 or Coot Hi, I have a question for everyone. I wanted to up-weight my high resolution terms in my electron density map by applying a B-factor sharpening term. I have mtz files produced from DM and I usually use Coot to generate my electron density maps. I was wondering is there an easy way to apply a B-factor sharpening term in Coot? If not, is there a program in ccp4 where by I can set the B-sharp factor? Thanks Jeff -- Notice: This e-mail message, together with any attachments, contains information of Merck & Co., Inc. (One Merck Drive, Whitehouse Station, New Jersey, USA 08889), and/or its affiliates (which may be known outside the United States as Merck Frosst, Merck Sharp & Dohme or MSD and in Japan, as Banyu - direct contact information for affiliates is available at http://www.merck.com/contact/contacts.html) that may be confidential, proprietary copyrighted and/or legally privileged. It is intended solely for the use of the individual or entity named on this message. If you are not the intended recipient, and have received this message in error, please notify us immediately by reply e-mail and then delete it from your system. --
[ccp4bb] B-factor sharpening in CCP4 or Coot
Hi, I have a question for everyone. I wanted to up-weight my high resolution terms in my electron density map by applying a B-factor sharpening term. I have mtz files produced from DM and I usually use Coot to generate my electron density maps. I was wondering is there an easy way to apply a B-factor sharpening term in Coot? If not, is there a program in ccp4 where by I can set the B-sharp factor? Thanks Jeff
Re: [ccp4bb] B-factor & Space gr questions!
On Tuesday 05 June 2007 12:19, Edward A Berry wrote: > You have a good point there and I would be interested in hearing > some other opinions, so I take the liberty of reposting- > > My instinctive preference is that each structure should be > supported solely by the data that is deposited with it - > (one dataset one structure) but in terms of good science > we want to produce the best model we can, and that might be > the rigid-body-located structure from another dataset. I don't think that is quite the right way to look at it. In general we refine our model so that it both - agrees with the data - agrees with a priori knowledge In maximum likelihood terms: we want to find the model that is the most likely explanation for our observed data. An inherently unlikely model is also an inherently unlikely explanation. Therefore we focus on likely models. We impose geometric restraints because we believe that we have a better a priori expectation for bond lengths and angles than can be determined de novo from the data in this one experiment. Similarly we impose the known sequence of our protein on the model, even if the maps are not sufficiently good to identify each amino acid directly from the electron density. If we have an a priori expectation for the conformation of the whole protein, or large pieces of it, then we should account for this in the model, even if the data is not sufficiently good to reproduce this expectation de novo. Therefore if you have a high-resolution structure available, the best treatment of low-resolution data may well be to place the known structure as a rigid body. If you suspect hinge motions or other large scale inter-domain shifts, you might want to refine the hinge angle explicitly, but unfortunately our usual refinement programs are not really set up for this. These are important issues, and are close to the heart of the Maximum Likelihood approach to model refinement. Ethan > cdekker wrote: > > Hi, > > > > Your reply to the ccp4bb has confused me a bit. I am currently refining > > a low res structure and realise that I don't know what to expect for > > final R and Rfree - it is definitely not what most people would publish. > > So the absolute values of R and Rfree are not telling me much, the only > > gauge I have is that as long as both R and Rfree are decreasing I am > > improving the model (and yes, at the moment that is only rigid body > > refinement). > > In your email reply you suggest that even though a refinement to > > convergence that will lead to an increased Rfree (and lower R? - a > > classic case of overfitting!) would be a better model than the > > rigid-body-refined only model. This is what confuses me. > > I can see your reasoning that starting with an atomic model to solve > > low-res data can lead to this behaviour, but then should the solution > > not be a modification of the starting model (maybe high B-factors?) to > > compensate for the difference in resolution of model and data? > > > > Carien > > > > On 4 Jun 2007, at 19:38, Edward A Berry wrote: > > > >> Ibrahim M. Moustafa wrote: > >>> The last question: In the same paper, for the complex structure R and > >>> Rfree are equal (30%) is that an indication for improper refinement > >>> in these published structure? I'd love to hear your comments on that > >>> too. > >> Several times I solved low resolution structures using high resolution > >> models, and noticed that R-free increased during atomic positional > >> refinement. This could be expected from the assertion that after > >> refinement to convergence, the final values should not depend on > >> the starting point: If I had started with a crude model and refined > >> against low resolution data, Rfree would not have gone as low as the > >> high-resolution model, so if I start with the high resolution model > >> and refine, Rfree should worsen to the same value as the structure > >> converges to the same point. > >> > >> Thinking about the main purpose of the Rfree statistic, in a very > >> real way this tells me that the model was better before this step > >> of refinement, and it would be better to omit the minimization step. > >> Perhaps this is what the authors did. > >> > >>On the other hand it does not seem quite right submit a model that > >> has simply been rigid-body-refined against the data- I would prefer to > >> refine to convergence and submit the best model that can be supported > >> by the data alone, rather than a better model which is really the model > >> from a better dataset repositioned in the new crystal. > >> > >> Ed > > > > > > The Institute of Cancer Research: Royal Cancer Hospital, a charitable > > Company Limited by Guarantee, Registered in England under Company No. > > 534147 with its Registered Office at 123 Old Brompton Road, London SW7 3RP. > > > > This e-mail message is confidential and for use by the addressee only. > > If the message is received by any
Re: [ccp4bb] B-factor & Space gr questions!
I think the relevant point in this discussion is that the original paper discussed the apo and substrate complexes of the protein. For the structure with lower resolution data you may indeed get a better model by taking the high resolution model and just applying rigid body refinement to it. After that step you would like to find and model the differences between the two structures. This includes the bound substrate (or the lack thereof) and any significant structural changes that accompany substrate binding. Significant meaning those changes that can be reliably determined at the lower resolution. For most of the structure that may mean you are best off by simply taking the rigid-body-refined coordinates of the higher resolution structure without further refinement. I see no problem in doing so and as long as interesting differences between the structures can be clearly defined and the procedure is explicity described in publications this should be perfectly reasonable. Bart Edward A Berry wrote: You have a good point there and I would be interested in hearing some other opinions, so I take the liberty of reposting- My instinctive preference is that each structure should be supported solely by the data that is deposited with it - (one dataset one structure) but in terms of good science we want to produce the best model we can, and that might be the rigid-body-located structure from another dataset. In particular the density for the ligand might be clearer before overfitting with the low resolution data. Even if the free-R set is not preserved for the new crystal, R and R-free tend to diverge rapidly once any kind of fitting with a low data/param is performed, so I think the new structure must not have been refined much beyond rigid body (and over-all B which is included in any kind of refinement). And that choice may be well justified. Ed cdekker wrote: Hi, Your reply to the ccp4bb has confused me a bit. I am currently refining a low res structure and realise that I don't know what to expect for final R and Rfree - it is definitely not what most people would publish. So the absolute values of R and Rfree are not telling me much, the only gauge I have is that as long as both R and Rfree are decreasing I am improving the model (and yes, at the moment that is only rigid body refinement). In your email reply you suggest that even though a refinement to convergence that will lead to an increased Rfree (and lower R? - a classic case of overfitting!) would be a better model than the rigid-body-refined only model. This is what confuses me. I can see your reasoning that starting with an atomic model to solve low-res data can lead to this behaviour, but then should the solution not be a modification of the starting model (maybe high B-factors?) to compensate for the difference in resolution of model and data? Carien On 4 Jun 2007, at 19:38, Edward A Berry wrote: Ibrahim M. Moustafa wrote: The last question: In the same paper, for the complex structure R and Rfree are equal (30%) is that an indication for improper refinement in these published structure? I'd love to hear your comments on that too. Several times I solved low resolution structures using high resolution models, and noticed that R-free increased during atomic positional refinement. This could be expected from the assertion that after refinement to convergence, the final values should not depend on the starting point: If I had started with a crude model and refined against low resolution data, Rfree would not have gone as low as the high-resolution model, so if I start with the high resolution model and refine, Rfree should worsen to the same value as the structure converges to the same point. Thinking about the main purpose of the Rfree statistic, in a very real way this tells me that the model was better before this step of refinement, and it would be better to omit the minimization step. Perhaps this is what the authors did. On the other hand it does not seem quite right submit a model that has simply been rigid-body-refined against the data- I would prefer to refine to convergence and submit the best model that can be supported by the data alone, rather than a better model which is really the model from a better dataset repositioned in the new crystal. Ed The Institute of Cancer Research: Royal Cancer Hospital, a charitable Company Limited by Guarantee, Registered in England under Company No. 534147 with its Registered Office at 123 Old Brompton Road, London SW7 3RP. This e-mail message is confidential and for use by the addressee only. If the message is received by anyone other than the addressee, please return the message to the sender by replying to it and then delete the message from your computer and network. -- == Bart Hazes (Assistant Professor) Dept. of Medical Microbiology & Immunology U
Re: [ccp4bb] B-factor & Space gr questions!
Wouldn't the desirability of this depend on the extent to which the molecule has moved between the high-resolution and low-resolution datasets ? I would have thought that there was an effective information transfer between R-work and R-free once the rigid body movements became too large, which might provide one with an over-optimistic idea of what the R-free would be with the high-resolution model with the low-resolution data. Phil Princeton NJ Edward A Berry wrote: Even if the free-R set is not preserved for the new crystal, R and R-free tend to diverge rapidly once any kind of fitting with a low data/param is performed, so I think the new structure must not have been refined much beyond rigid body (and over-all B which is included in any kind of refinement). And that choice may be well justified. Ed cdekker wrote: Hi, Your reply to the ccp4bb has confused me a bit. I am currently refining a low res structure and realise that I don't know what to expect for final R and Rfree - it is definitely not what most people would publish. So the absolute values of R and Rfree are not telling me much, the only gauge I have is that as long as both R and Rfree are decreasing I am improving the model (and yes, at the moment that is only rigid body refinement). In your email reply you suggest that even though a refinement to convergence that will lead to an increased Rfree (and lower R? - a classic case of overfitting!) would be a better model than the rigid-body-refined only model. This is what confuses me. I can see your reasoning that starting with an atomic model to solve low-res data can lead to this behaviour, but then should the solution not be a modification of the starting model (maybe high B-factors?) to compensate for the difference in resolution of model and data? Carien On 4 Jun 2007, at 19:38, Edward A Berry wrote: Ibrahim M. Moustafa wrote: The last question: In the same paper, for the complex structure R and Rfree are equal (30%) is that an indication for improper refinement in these published structure? I'd love to hear your comments on that too. Several times I solved low resolution structures using high resolution models, and noticed that R-free increased during atomic positional refinement. This could be expected from the assertion that after refinement to convergence, the final values should not depend on the starting point: If I had started with a crude model and refined against low resolution data, Rfree would not have gone as low as the high-resolution model, so if I start with the high resolution model and refine, Rfree should worsen to the same value as the structure converges to the same point. Thinking about the main purpose of the Rfree statistic, in a very real way this tells me that the model was better before this step of refinement, and it would be better to omit the minimization step. Perhaps this is what the authors did. On the other hand it does not seem quite right submit a model that has simply been rigid-body-refined against the data- I would prefer to refine to convergence and submit the best model that can be supported by the data alone, rather than a better model which is really the model from a better dataset repositioned in the new crystal. Ed
Re: [ccp4bb] B-factor & Space gr questions!
You have a good point there and I would be interested in hearing some other opinions, so I take the liberty of reposting- My instinctive preference is that each structure should be supported solely by the data that is deposited with it - (one dataset one structure) but in terms of good science we want to produce the best model we can, and that might be the rigid-body-located structure from another dataset. In particular the density for the ligand might be clearer before overfitting with the low resolution data. Even if the free-R set is not preserved for the new crystal, R and R-free tend to diverge rapidly once any kind of fitting with a low data/param is performed, so I think the new structure must not have been refined much beyond rigid body (and over-all B which is included in any kind of refinement). And that choice may be well justified. Ed cdekker wrote: Hi, Your reply to the ccp4bb has confused me a bit. I am currently refining a low res structure and realise that I don't know what to expect for final R and Rfree - it is definitely not what most people would publish. So the absolute values of R and Rfree are not telling me much, the only gauge I have is that as long as both R and Rfree are decreasing I am improving the model (and yes, at the moment that is only rigid body refinement). In your email reply you suggest that even though a refinement to convergence that will lead to an increased Rfree (and lower R? - a classic case of overfitting!) would be a better model than the rigid-body-refined only model. This is what confuses me. I can see your reasoning that starting with an atomic model to solve low-res data can lead to this behaviour, but then should the solution not be a modification of the starting model (maybe high B-factors?) to compensate for the difference in resolution of model and data? Carien On 4 Jun 2007, at 19:38, Edward A Berry wrote: Ibrahim M. Moustafa wrote: The last question: In the same paper, for the complex structure R and Rfree are equal (30%) is that an indication for improper refinement in these published structure? I'd love to hear your comments on that too. Several times I solved low resolution structures using high resolution models, and noticed that R-free increased during atomic positional refinement. This could be expected from the assertion that after refinement to convergence, the final values should not depend on the starting point: If I had started with a crude model and refined against low resolution data, Rfree would not have gone as low as the high-resolution model, so if I start with the high resolution model and refine, Rfree should worsen to the same value as the structure converges to the same point. Thinking about the main purpose of the Rfree statistic, in a very real way this tells me that the model was better before this step of refinement, and it would be better to omit the minimization step. Perhaps this is what the authors did. On the other hand it does not seem quite right submit a model that has simply been rigid-body-refined against the data- I would prefer to refine to convergence and submit the best model that can be supported by the data alone, rather than a better model which is really the model from a better dataset repositioned in the new crystal. Ed The Institute of Cancer Research: Royal Cancer Hospital, a charitable Company Limited by Guarantee, Registered in England under Company No. 534147 with its Registered Office at 123 Old Brompton Road, London SW7 3RP. This e-mail message is confidential and for use by the addressee only. If the message is received by anyone other than the addressee, please return the message to the sender by replying to it and then delete the message from your computer and network.
Re: [ccp4bb] B-factor & Space gr questions!
Hi All, Thanks a lot for all reply with valuable inputs. In my original post: I meant a = b "does not equal" c. I used # for "does not equal". Many asked where is that paper published! Actually the paper is under revision. When reading, I assumed the unit cell dimensions (or the space group) is a typo as others thought. The low B value for the low resolution structure makes me suspicious that something is wrong. In my little experience, and as others mentioned, B-factor is expected to be around 70-80 for 2.8 A structure and very likely higher for 3.0 A structure. David Briggs suggested that they reported the Wilson B-factor; however, clearly, it is reported as the B-factor of the refined structure. Also, the Rwork = Rfree indicated that something is not right with the refinement protocol but I was not sure what that could be! The suspect that they did not transfer the FreeR sounds reasonable explanation. thanks a lot, Ibrahim At 03:50 AM 6/5/2007, Eleanor Dodson wrote: Yes; a==b for P6i - prob. a typo.. B factors at 3.2A are hard to fix - it will depend on scaling convention to some extent.. Can you download the data and re-run refinement for your own satisfaction. If R ==Rfree for the complex then I suspect they did not transfer the FreeR flags from the apo-protein data to the complex. Again if the data is available you may be able to check this. Eleanor Ibrahim M. Moustafa wrote: Hi all, While reading a crystallographic paper describing the structure of an apo-protein and its complex I noticed that the authors described the space goup as P6122 for the unit cell: a=141.9, b=143.9, c=380.4 ! Could this be considered as a typo or I'm missing something here! the requirement for the hexagonal is a = b # Cright? Another observation in that paper too: the B-factors for the 2.4 A and 3.2 A structures are 39 and 40?? Does this make sense to anyone?? The last question: In the same paper, for the complex structure R and Rfree are equal (30%) is that an indication for improper refinement in these published structure? I'd love to hear your comments on that too. thanks, Ibrahim -- Ibrahim M. Moustafa, Ph.D. Biochemistry and Molecular Biology Dept. 201 Althouse Lab., Uinversity Park Pennsylvania State University, PA16802 Tel. (814)863-8703 Fax. (814)865-7927 -- -- Ibrahim M. Moustafa, Ph.D. Biochemistry and Molecular Biology Dept. 201 Althouse Lab., University Park Pennsylvania State University, PA16802 Tel. (814)863-8703 Fax. (814)865-7927 --
Re: [ccp4bb] B-factor & Space gr questions!
Yes; a==b for P6i - prob. a typo.. B factors at 3.2A are hard to fix - it will depend on scaling convention to some extent.. Can you download the data and re-run refinement for your own satisfaction. If R ==Rfree for the complex then I suspect they did not transfer the FreeR flags from the apo-protein data to the complex. Again if the data is available you may be able to check this. Eleanor Ibrahim M. Moustafa wrote: Hi all, While reading a crystallographic paper describing the structure of an apo-protein and its complex I noticed that the authors described the space goup as P6122 for the unit cell: a=141.9, b=143.9, c=380.4 ! Could this be considered as a typo or I'm missing something here! the requirement for the hexagonal is a = b # Cright? Another observation in that paper too: the B-factors for the 2.4 A and 3.2 A structures are 39 and 40?? Does this make sense to anyone?? The last question: In the same paper, for the complex structure R and Rfree are equal (30%) is that an indication for improper refinement in these published structure? I'd love to hear your comments on that too. thanks, Ibrahim -- Ibrahim M. Moustafa, Ph.D. Biochemistry and Molecular Biology Dept. 201 Althouse Lab., Uinversity Park Pennsylvania State University, PA16802 Tel. (814)863-8703 Fax. (814)865-7927 --
Re: [ccp4bb] B-factor & Space gr questions!
Hi Ibrahim, On 04/06/07, Ibrahim M. Moustafa <[EMAIL PROTECTED]> wrote: Hi all, While reading a crystallographic paper describing the structure of an apo-protein and its complex I noticed that the authors described the space goup as P6122 for the unit cell: a=141.9, b=143.9, c=380.4 Could this be considered as a typo or I'm missing something here! the requirement for the hexagonal is a = b # Cright? You are correct, for Hexagonal, a=b - so It's got to be a typo - data most processing software wouldn't let you do this. Another observation in that paper too: the B-factors for the 2.4 A and 3.2 A structures are 39 and 40?? Does this make sense to anyone? They're quoting Wilson B-factors, I imagine. A small but rather important difference - where was this published? The last question: In the same paper, for the complex structure R and Rfree are equal (30%) is that an indication for improper refinement in these published structure? I'd love to hear your comments on that too. Well, it certainly is a little suspicious looking - I've had similar experiences to Ed, regarding similar R & Rfrees from rigid rigid body refinement prior to positional refinement. Have the authors deposited the Structure factors? I would use EDS to check the maps out: eds.bmc.uu.se/eds/ thanks, Ibrahim HTH, Dave -- --- David Briggs, PhD. Father & Crystallographer www.dbriggs.talktalk.net iChat AIM ID: DBassophile --- Anyone who is capable of getting themselves made President should on no account be allowed to do the job. - Douglas Adams
Re: [ccp4bb] B-factor & Space gr questions!
Ibrahim M. Moustafa wrote: The last question: In the same paper, for the complex structure R and Rfree are equal (30%) is that an indication for improper refinement in these published structure? I'd love to hear your comments on that too. Several times I solved low resolution structures using high resolution models, and noticed that R-free increased during atomic positional refinement. This could be expected from the assertion that after refinement to convergence, the final values should not depend on the starting point: If I had started with a crude model and refined against low resolution data, Rfree would not have gone as low as the high-resolution model, so if I start with the high resolution model and refine, Rfree should worsen to the same value as the structure converges to the same point. Thinking about the main purpose of the Rfree statistic, in a very real way this tells me that the model was better before this step of refinement, and it would be better to omit the minimization step. Perhaps this is what the authors did. On the other hand it does not seem quite right submit a model that has simply been rigid-body-refined against the data- I would prefer to refine to convergence and submit the best model that can be supported by the data alone, rather than a better model which is really the model from a better dataset repositioned in the new crystal. Ed
[ccp4bb] B-factor & Space gr questions!
Hi all, While reading a crystallographic paper describing the structure of an apo-protein and its complex I noticed that the authors described the space goup as P6122 for the unit cell: a=141.9, b=143.9, c=380.4 ! Could this be considered as a typo or I'm missing something here! the requirement for the hexagonal is a = b # Cright? Another observation in that paper too: the B-factors for the 2.4 A and 3.2 A structures are 39 and 40?? Does this make sense to anyone?? The last question: In the same paper, for the complex structure R and Rfree are equal (30%) is that an indication for improper refinement in these published structure? I'd love to hear your comments on that too. thanks, Ibrahim -- Ibrahim M. Moustafa, Ph.D. Biochemistry and Molecular Biology Dept. 201 Althouse Lab., Uinversity Park Pennsylvania State University, PA16802 Tel. (814)863-8703 Fax. (814)865-7927 --