Entropy in Biomolecular Systems
DACAM, Max F. Perutz Laboratories, University of Vienna, Dr. Bohrgasse 9, A-1030, Vienna, Austria
Venue: Max F. Perutz Laboratories, University of Vienna, Dr. Bohrgasse 9, Vienna, A-1030, Austria.
Registration Fee: 350 EUR
Workshop registration 21 April 2014
Registration payment deadline 30 April 2014
Summary: This workshop brings together the world's experts to address the challenges of determining the entropy of biomolecular systems, either by experiment or computer simulation. Entropy is one the main driving forces for any biological process such as binding, folding, partitioning and reacting. Our deficient understandng of entropy, however, means that such important processes remain controversial and only partially understood. Contributions of water, ions, cofactors, and biomolecular flexibility are actively examined but yet to be resolved. The state-of-the-art of each entropy method will be presented and explained, highlighting its capabilities and deficiencies. This will be followed by intensive discussion on the main areas that need improving, leading suitable actions and collaborations to address the main biological and industrial questions.
Introduction and Motivation
Entropy and enthalpy comprise the two main driving forces for any physical, chemical or biological process. Taken together, these two quantities constitute the free energy, which determines the extent of any process and, as a free energy barrier, the rate at which it occurs. Historically, a significant part of the effort in molecular simulation has been dedicated to the calculation of free energy differences including a long-standing series of CECAM workshops related to this topic. However, separating entropic and enthalpic contributions is key to understanding the nature of phenomena such as the hydrophobic effect, protein folding or ligand binding. Energy quantifies a system’s molecular interactions and entropy quantifies its structural variation. While energy can be readily obtained from a simulation, the determination of entropy, and, in particular, its conformational and solvent components, is much more challenging than calculating a free energy difference. This is to a large extent due to the fact that for the calculation of absolute and relative entropies one in principle needs full coverage of the phase space. This being essentially impossible because of the astronomical size of phase space, ingenious ways are required to make this feasible. There are numerous techniques that qualitatively probe a system’s flexibility and structural heterogeneity as measured in a simulation or experimentally. However, determining the quantitative and complete link between structure and entropy remains a major unsolved problem. This is especially so for biomolecules, which, together with the surrounding aqueous medium, display considerable structural heterogeneity. Solving this problem will substantially improve our ability to understand the intricacy and subtlety of biological molecular function as well as rationally design our own molecular systems.
The recent years have seen the emergence of a number of promising computational and experimental techniques to address entropy-related issues in various biologically relevant processes. However, major challenges of a fundamental nature exist on both fronts, as detailed below, before entropy calculations can become routinely used in biomolecular systems. Surprisingly, despite the fundamental nature of entropy and the inadequacy of current approaches, very few scientific meetings have been entirely dedicated to its determination and its role in biomolecular systems. By bringing together theorists, computational scientists and experimentalists, the proposed workshop “Entropy in Biomolecular Systems” aims to fill this gap and play a catalytic role in delineating the main challenges in the field, enabling cross-talk across the theory/experiment divide and stimulating active collaborations. We believe that such a meeting would play a transformative role in establishing biomolecular entropy, as well as entropy calculations more generally, as a rich, propulsive, and multidisciplinary field of research.
The important and extremely diverse roles played by entropic effects in different biologically relevant processes are being recognized with an ever-increasing frequency. Here, we mention three examples. Firstly, in the area of rational drug design, it has long been recognized that solvent entropy and conformational entropy changes upon binding of a ligand to a biomolecule represent an important and yet still largely uncharted space of design possibilities (1-4). Understanding the basic quantitative principles that define this space may propel rational drug design to completely new levels of success. Secondly, in protein folding, chaperones can be understood as agents, which decrease the entropy of the unfolded state by reducing the number of conformations available to the protein. Any quantitative assessment of the chaperone’s role must accurately take this entropic contribution into account (5). Finally, dynamically driven allostery in which ligand binding affects protein behavior at distant sites is now seen as the main mechanism underlying signaling in the absence of clear structural changes (6, 7). Again, quantifying conformational entropy changes is key to understanding allostery, as well as many other biomolecular processes underlying biological function.
Despite such clear and compelling motivation, theoretical treatment of conformational and solvent entropies in particular is still largely underdeveloped, and the same goes for experimental measurement of these quantities. Part of the problem is a lack of recognition of the importance of entropy and the inadequacy of the current methods for its determination. One only has to look at the advanced and sophisticated theories that have been developed for electronic structure calculations of energy to realize how inferior is our theory for entropy. When it comes to the prediction of entropies, a combination of theory and computer simulation is required in all but the simplest cases of ideal gases, polymers or crystals. Depending on the problem at hand, the key quantity in question may be either the relative entropy between initial and final states of a given process or the absolute entropy of a given state. Changes in entropy can be derived using computer simulation from the free energy and enthalpy or from the temperature derivative of the free energy (8). While these methods are widely studied and widely used, they require sampling along a pathway between the two states, something that becomes computationally unfeasible for complex systems such as biomolecules. Moreover, they do not give the separation of entropy into its structural components, and most importantly its conformational and solvent parts.
The two traditional methods to calculate the absolute conformational entropy of a molecule are normal mode analysis (9-11) and quasi-harmonic analysis (12-16). Both of these approaches assume that the conformational state is restricted to a single basin of the potential energy surface, which is locally approximated by a harmonic oscillator. Normal mode analysis derives the associated force constants from the Hessian matrix. This is the preferred method in the electronic structure community because it only requires a single optimized structure (11). On the other hand, quasi-harmonic analysis derives force constants from the variance-covariance matrix of the displacements measured in constant temperature molecular dynamics simulations (12-16). Both normal mode and quasi-harmonic methods fail if the system is inherently anharmonic or if non-linear correlations exists between atomic motions, both conditions being frequently met in large biomolecules. The contributions of the surrounding solvent and ions have commonly been treated using continuum and empirical terms, which do not adequately capture local solvation effects (17, 18). Finally, while changes in total entropy are experimentally accessible in many cases using calorimetric techniques such as isothermal titration calorimetry, they remain silent when it comes to entropy components (19).
Over the past twenty years, a number of new methods have been proposed both computationally and experimentally to treat these inadequacies (16, 20). This is particularly true for biomolecular systems because of their central scientific importance and because of the conspicuous shortcomings of the traditional methods for treating entropy in such systems. On the simulation front, recent methodological developments account for the multiple conformational states, the anharmonic nature of these states, the presence of correlated motions and the entropy of the solvent (15, 21-32). In particular, different corrections have been developed to go beyond the quasi-harmonic approximation including, for example, those in third moments of coordinates (21) or in pairwise, linear correlations (15, 27). Most importantly, mutual information expansion techniques have been developed to allow for a treatment of supra-linear correlations, albeit with notorious convergence challenges (15, 22, 26, 30). For example, the minimally coupled subspace approach combines full correlation analysis, adaptive kernel density estimation and mutual information expansion to give accurate conformational entropy estimates at least for small peptides (31). Finally, important advances have also been made when it comes to the treatment of solvent entropy. From cell theory-based approaches (23, 25) to morphometric and integral equation theory (24, 29) to permutation reduction techniques (28), different methods have been developed in this context. While all of these approaches represent welcome advances in conformational and solvent entropy treatment, numerous computational and conceptual bottlenecks still prevent their wider usage. One of the main goals of the meeting proposed herein would be to define, categorize and prioritize these bottlenecks and delineate strategies that might lead to their resolution.
On the other hand, advances in the nuclear magnetic resonance (NMR) techniques, most notably relaxation techniques of multinuclear and multidimensional NMR, have recently led to a revolutionary new view of the role of dynamics in general and conformational entropy, in particular, in biomolecular processes (7, 33-39). Specifically, most progress has been made with NMR by deriving conformational entropy from generalized order parameters of isolated bond vectors or methyl symmetry axes (33-36, 39). Importantly, however, NMR analyses must rely on dynamical models in order to give a parametric relationship between different NMR observables, physical rearrangements and conformational entropy (33-36, 39). The absolute entropies obtained in this way very much depend on the particulars of the potential energy function used. What is more, the extremely important correlated motions are typically ignored in such approaches, primarily because they are inaccessible to most experiments. On the other hand, molecular dynamics simulations not only help interpret NMR relaxation experiments in terms of entropic contributions, but also give a detailed, atomistic picture of such contributions even on their own (40). As both simulation and experiment converge on calculating the same quantities, the proposed meeting will provide an exciting and valuable platform for comparing the approximations and assumptions of the two approaches, exploring the potential synergies, and making new methodologies and capabilities available to the wider community.
University of Vienna, Dr. Bohrgasse 9, Vienna, A-1030, Austria
Chris Oostenbrink (University of Natural Resources and Life Sciences) - Organiser & speaker
Bojan Zagrovic (University of Vienna) - Organiser
Michel Cuendet (Swiss Institute of Bioinformatics / Division of Analytical Personalized Oncology, Lausanne University Hospital / Department of Physiology and BIophysics, Weill Cornell Medicine) - Organiser & speaker
Richard Henchman (University of Manchester) - Organiser & speaker