{"id":296,"date":"2019-07-12T17:39:43","date_gmt":"2019-07-12T17:39:43","guid":{"rendered":"http:\/\/statisticalbiophysicsblog.org\/?p=296"},"modified":"2019-07-12T17:39:45","modified_gmt":"2019-07-12T17:39:45","slug":"learn-your-lesson-every-method-has-its-limits","status":"publish","type":"post","link":"https:\/\/statisticalbiophysicsblog.org\/?p=296","title":{"rendered":"Learn your lesson: Every method has its limits"},"content":{"rendered":"<p>It was at breakfast during a recent conference that Prof. X leaned toward me and quietly said, \u201cWe tried using weighted ensemble but it didn\u2019t work.\u201d  I got the sense he was trying not to broadcast this to other conference attendees, as a courtesy.<\/p>\n<p><!--more--><\/p>\n<p>Needless to say, I pressed him on what he meant.  After all, the weighted ensemble (WE) method has been shown to work for quite a number of challenging systems \u2013 including ligand binding\/unbinding and protein folding.  On the other hand, in my group, we\u2019ve seen WE fall short in notable instances: protein folding from certain initial conformations and for both unbinding and conformational change in the estrogen receptor\u2019s ligand binding domain (LBD).  I\u2019ve been interested in the estrogen receptor for a number of years, but it\u2019s a tricky system.<\/p>\n<p>By coincidence, it turned out Prof. X had attempted to study the famous conformational change in the estrogen receptor LBD, when one of its helices undergoes a huge conformational swing.  But no luck.<\/p>\n<p>Are such failures something to whisper about, only?  Absolutely not.<\/p>\n<p>Every method has its limitations, and these fall into two categories.  There are general limitations common to a broad class of approaches.  For PMF calculations, for instance, there is the generic issue of sampling orthogonal coordinates, as <a href=\"https:\/\/statisticalbiophysicsblog.org\/?p=160\">discussed previously<\/a>.  For path sampling approaches like WE, there is the intrinsic cost of the path ensemble itself, which includes the path length and the number of instances obtained, as was also <a href=\"https:\/\/statisticalbiophysicsblog.org\/?p=115\">noted early on in this blog<\/a>.  Then there are limitations specific to a particular method \u2013 e.g., based on the binning scheme in WE.  We have tried to be open about the challenges of using WE, for instance in a <a href=\"https:\/\/www.annualreviews.org\/doi\/full\/10.1146\/annurev-biophys-070816-033834\">review article<\/a> and in our <a href=\"https:\/\/westpa.github.io\/westpa\/overview.html\">online WESTPA documentation<\/a>.<\/p>\n<p>And yet, generic or abstract warnings may not be enough.  They\u2019re not so scary when you\u2019re feeling ambitious and have a large compute allocation in your hands.<\/p>\n<p>I think we need concrete examples to learn from, and that is the motivation for the \u201c<a href=\"https:\/\/www.livecomsjournal.org\/post\/130\">Lessons Learned<\/a>\u201d category of the new <a href=\"https:\/\/www.livecomsjournal.org\/\">LiveCoMS<\/a> journal, which I helped Michael Shirts and David Mobley in founding.  Let\u2019s face it: our community would greatly benefit from airing out the <em>numerous<\/em> shortcomings of prior efforts \u2013 most of which, I believe, are buried in <em>didntwork\/<\/em> directories.  We are doomed to burn ever more computing time attempting impossible tasks.  Ten times more computing won\u2019t help if pertinent sampling times have been under-estimated by multiple orders of magnitude.<\/p>\n<p>I strongly encourage you to consider the value of your failed efforts in educating the community.  You might even get a highly cited paper out of it!  Help us learn our lessons!<\/p>\n","protected":false},"excerpt":{"rendered":"<p>It was at breakfast during a recent conference that Prof. X leaned toward me and quietly said, \u201cWe tried using weighted ensemble but it didn\u2019t work.\u201d I got the sense he was trying not to broadcast this to other conference attendees, as a courtesy.<\/p>\n","protected":false},"author":6,"featured_media":298,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[14,18],"tags":[],"class_list":["post-296","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-path-sampling","category-statistical-uncertainty"],"_links":{"self":[{"href":"https:\/\/statisticalbiophysicsblog.org\/index.php?rest_route=\/wp\/v2\/posts\/296","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/statisticalbiophysicsblog.org\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/statisticalbiophysicsblog.org\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/statisticalbiophysicsblog.org\/index.php?rest_route=\/wp\/v2\/users\/6"}],"replies":[{"embeddable":true,"href":"https:\/\/statisticalbiophysicsblog.org\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=296"}],"version-history":[{"count":3,"href":"https:\/\/statisticalbiophysicsblog.org\/index.php?rest_route=\/wp\/v2\/posts\/296\/revisions"}],"predecessor-version":[{"id":300,"href":"https:\/\/statisticalbiophysicsblog.org\/index.php?rest_route=\/wp\/v2\/posts\/296\/revisions\/300"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/statisticalbiophysicsblog.org\/index.php?rest_route=\/wp\/v2\/media\/298"}],"wp:attachment":[{"href":"https:\/\/statisticalbiophysicsblog.org\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=296"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/statisticalbiophysicsblog.org\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=296"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/statisticalbiophysicsblog.org\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=296"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}