To my knowledge, David Chandler and coworkers were the first to highlight dangers of the PMF. They described energy landscapes \u2013 potential energy functions \u2013 with mechanisms that could not be described by obvious coordinate choices for PMF calculations (thus motivating path-sampling techniques<\/a>). Expanding on those ideas, consider the two-dimensional landscapes below.<\/p>\n Although some of the landscapes perhaps could be described by some (tortuous) one-dimensional coordinate, in other instances that just wouldn\u2019t be possible. Think how much worse the situation would have to be in a system like a protein with thousands of degrees of freedom.<\/p>\n Now take one more step on the road to skepticism and consider how PMFs tend to be constructed in computational studies. Because the calculations are so expensive, we tend to choose one or two coordinates in advance for study. These coordinate choices represent our pre-conceived ideas about how a system might function.<\/p>\n Based on the features of the PMF landscape, we build a story about mechanism \u2026 the dwells are here and the transition states there. And yet we know that all basins are not alike \u2013 some may be energy stabilized and other favored by entropy. The heights of barriers, their widths and likewise the dimensions of basins all are influenced by the particular coordinate(s) we chose in advance, as illustrated in the two-dimensional examples above.<\/p>\n So I would say the quantitative analysis of a PMF is inherently misleading, but the \u201cstories\u201d we build are almost more dangerous. Because our minds are naturally drawn to stories, such narratives are easy to retain. This wouldn\u2019t be a problem if the stories were not biased by the subjective coordinate choices and perhaps also by inadequate sampling.<\/p>\n Is there a good way to think about mechanism? I believe the starting point has to be a trajectory ensemble of continuous unbiased transition events<\/a>. Trajectories tell us the full sequence of events from an initial to a final state. And the ensemble reports on the diversity of mechanisms, something a PMF could never do. It\u2019s true that trajectory ensembles are difficult to obtain<\/a>, but at least they offer a potentially unbiased way to describe mechanism.<\/p>\n As a final perspective on mechanism, think about the goal of discovery. We use large amounts of computing resources and we would like to be able to discover something we did not know before. But if we only look at coordinates we presumed from the start to be important, we are severely impairing our ability to discover novel phenomena.<\/p>\n Perhaps it\u2019s worth a moment to consider that our subjective coordinate choices typically are based on a further assumption \u2013 that we have good or complete knowledge of our system\u2019s biological function. Is this really true?<\/p>\n Ideally, although it\u2019s a challenge, we would perform unbiased analyses of path\/trajectory ensembles to discover important coordinates. Significant work has already been done on automated discovery of coordinates, and I think such methodologies should play an increasingly important role in the future.<\/p>\n Does a PMF predict kinetics?<\/em><\/p>\n We\u2019ve all done it \u2013 looked at a PMF (free energy landscape), estimated a barrier height, and tried to guess a rate constant for a process. If we haven\u2019t tried to guess an absolute rate, we\u2019ve at least said to ourselves something like, \u201cWell, that barrier is higher than the other, so one process is slower than the other.\u201d<\/p>\n But in principle, the PMF may not yield any information about kinetics at all! Look again at the two-dimensional landscapes above. Would you trust a one-dimensional PMF of any of these to give reliable rates? So why trust a projection from 1,000+ dimensions to one or two, as is usually done for biomolecular systems?<\/p>\n Again, I think part of the reason we over-interpret PMFs is because free energy landscapes speak to us like stories. We can\u2019t resist the narrative that arises in the stat-mech part of our brains, even when we know better.<\/p>\n To drive home the point further, note that it\u2019s even possible to construct a PMF that is exactly constant and yet does not exhibit diffusive dynamics. One way would be with a potential consisting of `anti-parallel\u2019 valleys as in landscape (c) above, where the well depths and widths were tuned at each x value to exactly have equal probability when integrated over y.<\/p>\n Another constant-in-x PMF based on a non-trivial potential energy in x and y results from the following potential energy, which is a double-well potential in x modulated by a harmonic y component with x-varying width:<\/p>\n You can check this yields a constant PMF by integrating over y and multiplying the result by the Boltzmann factor of the energy minimum at any x value. Below is a sample trajectory of x values (y not shown) simulated with overdamped Langevin (a.k.a. Brownian) dynamics.<\/p>\n Clearly the behavior is not simple diffusion in x, even though the PMF(x) is constant. (Histogramming a lot of these trajectories into x bins numerically confirms the PMF is flat.) And that\u2019s no surprise, from a theory standpoint. The PMF generally does not govern dynamics, except in rare cases where an ideal reaction coordinate has been used \u2026 and suitable effective dynamics have been specifically derived for the selected coordinate.<\/p>\n Are biomolecular PMFs well sampled?<\/em><\/p>\n Very sophisticated simulation and analysis techniques are used to calculate PMFs. But should we trust the results? To be clear, I\u2019m not questioning the rigor behind the methods, but rather the likelihood that the PMFs are well-sampled.<\/p>\n Consider the popular analysis approach WHAM (weighted histogram analysis method). WHAM seeks to provide the best possible PMF given the data available from the different simulation windows. This is a very different goal from assessing whether sufficient sampling has been performed. For better and for worse, WHAM almost always provides a smooth estimate of the free energy profile. Our minds tend to confuse smoothness with reliable, well-sampled data. But the two things are completely independent.<\/p>\n Let\u2019s be more concrete. How many nsec are typically used in a single window of a WHAM protein calculation? Often the answer seems to be 10 nsec or less. Is that really enough time to sample even a constrained protein process? When careful studies are done on small peptides, they suggest 10s of nsec are needed for good sampling of these tiny systems!<\/p>\n Now consider a thought experiment. If indeed every window of a WHAM-like PMF calculation is well sampled, then we have an equilibrium ensemble for that window. And if the PMF is accurate, we also know the relative free energy of that coordinate value. We therefore can generate an overall equilibrium ensemble by combining all the window ensembles, weighting each by the Boltzmann factor of the window-specific PMF. Would you trust this equilibrium ensemble? Further, an equilibrium ensemble can be projected onto any coordinate to generate a new PMF. Would you trust that new PMF?<\/p>\n Pessimism<\/em><\/p>\n To sum up: Be careful! I am dubious that protein PMFs are well-sampled. But even when a PMF is exactly accurate, it reveals a landscape based on a subjectively chosen coordinate set. Physical scientists are adept at building a story from a landscape, but any given landscape is likely to be deceptive both in terms of mechanism and kinetics. Perhaps it\u2019s time for us to move on from over-reliance on the PMF.<\/p>\n Acknowledgement<\/em>.<\/p>\n I very much appreciate comments on a draft given by Alan Grossfield.<\/p>\n Further reading<\/strong><\/p>\n Dellago, C.; Bolhuis, P. G.; Csajka, F. S. & Chandler, D., Transition path sampling and the calculation of rate constants, J. Chem. Phys., 1998, 108, 1964-1977<\/a><\/p>\n Grossfield, A. & Zuckerman, D. M., Quantifying uncertainty and sampling quality in biomolecular simulations, Annu Rep Comput Chem, 2009, 5, 23-48. Lyman, E. & Zuckerman, D. M., On the Structural Convergence of Biomolecular Simulations by Determination of the Effective Sample Size, J. Phys. Chem. B, 2007, 111, 12876-12882<\/a><\/p>\n McGibbon, Robert T., Brooke E. Husic, and Vijay S. Pande, Identification of simple reaction coordinates from complex dynamics, The Journal of Chemical Physics 146, 044109 (2017)<\/a><\/p>\n![\"\"](\"https:\/\/statisticalbiophysicsblog.org\/wp-content\/uploads\/2017\/05\/saddlepoint-300x210.png\")
<\/a>![\"\"](\"https:\/\/statisticalbiophysicsblog.org\/wp-content\/uploads\/2017\/05\/saddlepoint2-300x197.png\")
<\/a>![\"\"](\"https:\/\/statisticalbiophysicsblog.org\/wp-content\/uploads\/2017\/05\/saddlepoint3-300x202.png\")
<\/a><\/p>\n![\"\"](\"https:\/\/statisticalbiophysicsblog.org\/wp-content\/uploads\/2017\/05\/bias_equation-300x52.png\")
<\/a><\/p>\n![\"\"](\"https:\/\/statisticalbiophysicsblog.org\/wp-content\/uploads\/2017\/05\/bias_equation2-300x42.png\")
<\/a><\/p>\n![\"\"](\"https:\/\/statisticalbiophysicsblog.org\/wp-content\/uploads\/2017\/05\/bias_equation3-300x26.png\")
<\/a><\/p>\n![\"\"](\"https:\/\/statisticalbiophysicsblog.org\/wp-content\/uploads\/2017\/05\/bias_Langevin.png\")
<\/a><\/p>\n
\nKumar, S.; Rosenberg, J. M.; Bouzida, D.; Swendsen, R. H. & Kollman, P. A., Multidimensional free-energy calculations using the weighted histogram analysis method, J. Comput. Chem., 1995, 16, 1339-135<\/a><\/p>\n