Hello all,
While there is definitely scientific value in reaching perfect convergence by "cooking" refinement jobs "à point", in this day and age I think we should all reflect on whether this value is worth it, given that computing is also contributing to cooking the planet. I am not advocating for poorly refined models, but for finding a reasonable baseline with no more cooking than is necessary to achieve satisfactory models. I know, define "necessary" and "satisfactory"... they are not the same if you want a reasonable model to answer a biological question or if you want to benchmark refinement strategies... Here is an interesting resource if you worry about the sustainability of your computing: https://www.green-algorithms.org/ The calculator is nice too: https://calculator.green-algorithms.org/ Cheers, Guillaume --- Guillaume Gaullier, PhD Researcher, Blikstad group Molecular Biomimetics / Microbial Chemistry Department of Chemistry - Ångström Uppsala University Lägerhyddsvägen 1 752 37 Uppsala Sweden ________________________________ From: CCP4 bulletin board <CCP4BB@JISCMAIL.AC.UK> on behalf of Nigel Moriarty <nwmoria...@lbl.gov> Sent: Friday, January 19, 2024 2:10:42 AM To: CCP4BB@JISCMAIL.AC.UK Subject: Re: [ccp4bb] Automated refinement convergence I, too, love the smell of cooking (or cooked?) jobs in the morning. Cheers Nigel --- Nigel W. Moriarty Building 33R0349, Molecular Biophysics and Integrated Bioimaging Lawrence Berkeley National Laboratory Berkeley, CA 94720-8235 Email : nwmoria...@lbl.gov Web : CCI.LBL.gov<http://CCI.LBL.gov> ORCID : orcid.org/0000-0001-8857-9464<https://orcid.org/0000-0001-8857-9464> On Thu, Jan 18, 2024 at 1:16 PM James Holton <jmhol...@lbl.gov<mailto:jmhol...@lbl.gov>> wrote: Hey there Robert, Refmac has a keyword called "kill" that I think is what you are looking for. It is documented here: https://www2.mrc-lmb.cam.ac.uk/groups/murshudov/content/refmac/refmac_keywords.html You can specify a conditional exit based on R factor, etc. Or you can just create a specified file containing "stop Y" from an external process. I use it when running refmac on a cluster that has run time limits but difficult-to-predict CPU speeds. Phenix, I don't think has a checkpointing feature. Not that I know of. Amber does support checkpointing and now counts as a refinement program since support for structure factor restraints was added in 22. Personally, when I do refinements I do dozens to hundreds of macro-macro-cycles. As in, take the pdb file output by one run and feed it into another run. There is an instantiation overhead to doing this, as you note, but I like my models to be super converged. I define convergence as the x,y,z,B and occ values in the pdb file are not changed by the refinement program. This does not happen quickly, but it does eventually happen. Yes, you can get oscillations, but one way to deal with those is to add a bit more damping, or to adjust the x-ray weight down and then up and then back to auto again. This "weight snap" tends to take things that were dangling from a cliff in the energy landscape and knock them to the ground. After that, the oscillations are less common. And like an equilibrated chromatography column, an xyz-converged model is the best way to know that when you edit and re-refine, everything you see is due to the edit, and not some other process that just wasn't finished yet. That's what I do. Maybe I just want to feel like I've got something cooking while I sleep... Cheers, -James Holton MAD Scientist On 1/18/2024 3:04 AM, Robert Oeffner wrote: > Hi, > > I am wondering if authors of refinement programs would like to consider > putting on their users wish list the ability of refinement programs to > automatically terminate once the refinement has reached convergence. Various > refinement metrics such as R factors, CC or RMS values typically will reach a > plateau once the refinement of a macromolecular structure with X-ray or > EM-data has converged and further macro-cycles of refinement will no longer > improve the structure. The default number of macro-cycles in programs such as > Phenix-refine and Refmac are probably sensible for most cases but in some > cases it would be nice if the programs automatically extended the number of > macro-cycles as needed (or decreased the number). > > The user can of course examine log files from refinement themselves and > decide whether to continue refinement. But since starting a new session of > refinement appears to always create an initial fluctuation in the refinement > metrics before they align with the values of the last macro-cycles in the > previous refinement session, the user is compelled to do at least, say 3 or > more macrocycles in addition to whatever may be needed for reaching > convergence. I guess it would therefore be more efficient if this was > implemented directly in the refinement programs and presented as an option > for the user to choose. > > There could be cases where alternate conformations of a structure will > repeatedly be oscillating in and out of density thus causing the refinement > metrics also to oscillate. Hopefully such cases could be covered by gauging > the level of fluctuations of the refinement metrics and terminate the > refinement accordingly. > > Many thanks, > > Robert > > ######################################################################## > > 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<http://www.jiscmail.ac.uk/CCP4BB>, a mailing list > hosted by www.jiscmail.ac.uk<http://www.jiscmail.ac.uk>, terms & conditions > are available at https://www.jiscmail.ac.uk/policyandsecurity/ ######################################################################## 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<http://www.jiscmail.ac.uk/CCP4BB>, a mailing list hosted by www.jiscmail.ac.uk<http://www.jiscmail.ac.uk>, terms & conditions are available at https://www.jiscmail.ac.uk/policyandsecurity/ ________________________________ To unsubscribe from the CCP4BB list, click the following link: https://www.jiscmail.ac.uk/cgi-bin/WA-JISC.exe?SUBED1=CCP4BB&A=1 VARNING: Klicka inte på länkar och öppna inte bilagor om du inte känner igen avsändaren och vet att innehållet är säkert. CAUTION: Do not click on links or open attachments unless you recognise the sender and know the content is safe. När du har kontakt med oss på Uppsala universitet med e-post så innebär det att vi behandlar dina personuppgifter. För att läsa mer om hur vi gör det kan du läsa här: http://www.uu.se/om-uu/dataskydd-personuppgifter/ E-mailing Uppsala University means that we will process your personal data. For more information on how this is performed, please read here: http://www.uu.se/en/about-uu/data-protection-policy ######################################################################## 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/
