From trajectories to reaction coordinates: making sense of molecular simulation data
CECAM-AT, Erwin Schrödinger Institute, University of Vienna, Austria
The goal of this workshop is to define and advance the state of the art in constructing reliable and meaningful models from complex data, and in using such models to gain mechanistic insight as a guide for applications and further study.
Having appropriate reaction coordinates is key for the efficient sampling both in space and in time. Knowledge of the relevant degrees of freedom for transition between metastable states helps improve the efficiency in free energy calculations. Exploring the system along these most slowly relaxing degrees of freedom minimizes the hysteresis effects associated with unresolved slow degrees of freedom [Rosta 2009]. One can also use the reaction coordinates to increase the efficiency of the dynamic sampling of rare events in methods such as forward-flux sampling [Allen 2009] or string methods [E 2005], where the appropriate choice of coordinates is critical. Even for methods that do not rely on prior knowledge of reaction coordinates such as transition path sampling [Dellago 2002] and transition-interface sampling [van Erp 2003], improved coordinates are beneficial.
The identification of reliable order parameters and reaction coordinates also constitutes a key step in the construction of mechanistic models. They help define both the state space in a coarse-grained manner more accessible to humans (e.g., reducing the protein folding problem to transitions between just two states), and the dynamics in state space. By contrast, if one has only poor order parameters and reaction coordinates, the resulting models will not accurately capture the populations and their dynamics, since the metastable states will be mixed to together and the transitions between them poorly resolved.
A key element of the workshop will also be the interpretation and use of reaction coordinate information in areas such as control. An emerging practical application is in the design of novel pharmaceuticals that no longer aim to inhibit enzyme binding and function directly, in part to avoid cross-reactivity. Instead, often-critical allosteric couplings are targeted that relay the effects of distal binding and conformational dynamics [D’Abramo 2012]. Understanding the relevant reaction coordinates of associated conformational changes should aid in rational design approaches for allosteric inhibitors.
We expect much progress from adaptively optimized reaction coordinates in simulation applications such as targeted molecular dynamics [Ovchinnikov 2012] and ratchet-and-pawl simulations [Marchi 1999], where large-scale conformational changes are induced through time-dependent biases. Similarly, the efficient use of information about slow degrees of freedom in free energy calculations will be a major theme, with relevance in all areas of molecular simulation, from simple physical systems to biomolecules and complex materials. Importantly, the exploration of free energy surfaces and the sampling of rare events should benefit greatly from having good coordinates. How the expected gains can be realized will be an important theme of the workshop.
As much work has been done the development and use of Markov state models to describe complex processes, we expect a lively discussion on whether and when one should try identifying a single (possibly complicated) reaction coordinate, or try to construct an MSM using many metastable states. The development and application of transition path theory [Vanden-Eijnden 2006] plays a key role in this discussion.
Another important topic that will be discussed during the workshop is the use of reaction coordinates for dynamical coarse-graining. Since reaction coordinates (often) constitute the slow variables of a process, a natural coarse-graining of the dynamics would be to replace the many-particle system by a diffusive (overdamped Langevin) process on the reaction coordinate. When and how this procedure is allowed is hotly debated. Several novel developments in this area address these issues [Peters2013, Berezhkovskii 2013, Lu2014], and hence we expect a lively discussion about this topic.
With reaction coordinates being at the center of a wide range of problems, it will be important to pool the cumulative knowledge of diverse but related fields. Applied mathematicians interested in areas such as transition path theory, computer scientists in areas such as machine learning, physicists interested in dynamical systems and phase transitions, chemists working on reaction mechanisms, and biologists elucidating the function of molecules and cellular networks are representative of the many fields having contributed to and benefiting from an improved understanding of reaction coordinates. While the focus is mainly on theory and simulation, the discussion will also benefit from exchange with experimentalists. As the time and space resolution of experiments reaches the single-particle level (atoms, molecules, or colloids), experimentalist face the same problems as computer simulators in gaining insight from their data. A main goal of the meeting is to enhance communication in the community, with the hope of developing new ideas and computational approaches.
Christoph Dellago (University of Vienna) - Organiser
Gerhard Hummer (Max Planck Institute of Biophysics) - Organiser
Peter Bolhuis (University of Amsterdam) - Organiser