Reaction coordinates and transition states in enzymatic catalysis

Abstract

Enzymatic reactions are complex chemical processes taking place in complex dynamic environments. Theoretical characterization of these reactions requires the determination of the reaction coordinate and the transition state ensemble. This is not an easy task because many degrees of freedom may be involved in principle. We present recent efforts to find good enzymatic reaction coordinates and the implications of these findings in the interpretation of enzymatic efficiency. In particular, we analyze different strategies based on the use of minimum free energy paths and direct localization of the dividing surface on multidimensional free energy surfaces. Another strategy is based on the generation of reactive trajectories, using the transition path sampling method, from which transition state configurations can be harvested. Most of the applications carried out until now coincide to stress the change in the nature of the reaction coordinate, in terms of the participation of the chemical and environmental degrees of freedom, as the reaction advances. The degrees of freedom of the chemical system are dominant at the transition state while environmental participation can be more important at early or late stages of the process.

Impact of Warhead Modulations on the Covalent Inhibition of SARS-CoV-2 Mpro Explored by QM/MM Simulations

“…The potential of mean force for each chemical step was obtained using a combination of the umbrella sampling (US) approach 43 with the weighted histogram analysis method (WHAM). 44 A series of MD simulations were performed on adding a restraint along the collective reaction coordinate s , 45 with an umbrella force constant of 3000 kJ mol –1 Å –2 . In every window, QM/MM MD-NVT simulations were performed with a total of 3.25 ps at 310 K with a time step of 0.5 fs (a total of 6500 steps).…”

How a Second Mg2+ Ion Affects the Phosphoryl-Transfer Mechanism in a Protein Kinase: A Computational Study

“…These points, or nodes, move to the regions of lower free energy, converging to the MFEP. 74 In this way, the calculation of the complete free energy surface is not needed assuming that reactive trajectories rarely visit regions far from the MFEP. As in previous studies, the string method has been combined with the use of a path CV 75,76 , s coordinate in this work, that measures the progress along the reaction.…”

Optimal Paths Between Parameter Estimates in Non-linear ODE Systems Using the Nudged Elastic Band Method

“…In order to be able to further extend the usage of the waterfall plot by checking the connectivity of parameter estimates, efficient methods are necessary to identify also non-trivial connections between points in the parameter space. Identification of minimal energy paths and free energy barriers in high-dimensional landscapes is a well-known and actively discussed task, especially in the field of computational chemistry and physics of molecular dynamics for finding smooth transitions states of chemical components [27][28][29][30]. There, in analogy to the described ODE models, evaluation of the objective function is rather expensive in terms of computational resources, such that sampling the highdimensional parameter space is not practical and path finding remains a non-trivial task.…”