Generating conformational transition paths with low potential-energy barriers for proteins

International audienceThe knowledge of conformational transition paths in proteins can be useful for understanding protein mechanisms. Recently, we have introduced the As-Rigid-As-Possible (ARAP) interpolation method, for generating interpolation paths between two protein conformations. The method was shown to preserve well the rigidity of the initial conformation along the path. However, because the method is totally geometry-based, the generated paths may be inconsistent because the atom interactions are ignored. Therefore, in this article, we would like to introduce a new method to generate conformational transition paths with low potential-energy barriers for proteins. The method is composed of three processing stages. First, ARAP interpolation is used for generating an initial path. Then, the path conformations are enhanced by a clash remover. Finally, Nudged Elastic Band, a path-optimization method, is used to produce a low-energy path. Large energy reductions are found in the paths obtained from the method than in those obtained from the ARAP interpolation method alone. The results also show that ARAP interpolation is a good candidate for generating an initial path because it leads to lower potential-energy paths than two other common methods for path interpolation

A graph theoretical analysis of the energy landscape of model polymers

In systems characterized by a rough potential energy landscape, local energetic minima and saddles define a network of metastable states whose topology strongly influences the dynamics. Changes in temperature, causing the merging and splitting of metastable states, have non trivial effects on such networks and must be taken into account. We do this by means of a recently proposed renormalization procedure. This method is applied to analyze the topology of the network of metastable states for different polypeptidic sequences in a minimalistic polymer model. A smaller spectral dimension emerges as a hallmark of stability of the global energy minimum and highlights a non-obvious link between dynamic and thermodynamic properties.Comment: 15 pages, 15 figure

Exploration of Reaction Pathways and Chemical Transformation Networks

For the investigation of chemical reaction networks, the identification of all relevant intermediates and elementary reactions is mandatory. Many algorithmic approaches exist that perform explorations efficiently and automatedly. These approaches differ in their application range, the level of completeness of the exploration, as well as the amount of heuristics and human intervention required. Here, we describe and compare the different approaches based on these criteria. Future directions leveraging the strengths of chemical heuristics, human interaction, and physical rigor are discussed.Comment: 48 pages, 4 figure

Constraint methods for determining pathways and free energy of activated processes

Activated processes from chemical reactions up to conformational transitions of large biomolecules are hampered by barriers which are overcome only by the input of some free energy of activation. Hence, the characteristic and rate-determining barrier regions are not sufficiently sampled by usual simulation techniques. Constraints on a reaction coordinate r have turned out to be a suitable means to explore difficult pathways without changing potential function, energy or temperature. For a dense sequence of values of r, the corresponding sequence of simulations provides a pathway for the process. As only one coordinate among thousands is fixed during each simulation, the pathway essentially reflects the system's internal dynamics. From mean forces the free energy profile can be calculated to obtain reaction rates and insight in the reaction mechanism. In the last decade, theoretical tools and computing capacity have been developed to a degree where simulations give impressive qualitative insight in the processes at quantitative agreement with experiments. Here, we give an introduction to reaction pathways and coordinates, and develop the theory of free energy as the potential of mean force. We clarify the connection between mean force and constraint force which is the central quantity evaluated, and discuss the mass metric tensor correction. Well-behaved coordinates without tensor correction are considered. We discuss the theoretical background and practical implementation on the example of the reaction coordinate of targeted molecular dynamics simulation. Finally, we compare applications of constraint methods and other techniques developed for the same purpose, and discuss the limits of the approach

Theory and Simulation of Rare Events in Stochastic Systems

Activated processes driven by rare fluctuations are discussed in this thesis. Understanding the dynamics of these activated processes is important for understanding chemical and biological reactions, drug design and many other important applications. First, theoretical tools including the Langevin equation, the Fokker-Planck equation and the path integral technique are reviewed. Based on these theories, simulation methods have been developed to sample the activated processes by a number of investigators. Several of the most important path sampling and path generating approaches are introduced. A combination of analytic and numerical techniques are applied to study the distribution of the durations of transition events over a barrier in a one-dimensional system undergoing over-damped Langevin dynamics. Then we employ the ``weighted ensemble' path sampling method to generate an unbiased ensemble of paths for a conformational transition in a 210-dimensional model of the protein calmodulin, and also find the reaction rate. The results show that the weighted ensemble approach is a remarkably straightforward and successful method. At last, systems with multiple channels are studied by the weighted ensemble approach and the more common transition path sampling approach. The weighted ensemble method is distinguished by its ability to perform complete path sampling for systems with multiplechannels at reasonable cost