Abstract:Water-mediated proton transfer reactions are central for catalytic processes in a wide range of biochemical systems, ranging from biological energy conversion to chemical transformations in the metabolism. Yet, the accurate computational treatment of such complex bio-chemical reactions is highly challenging and requires the application of multiscale methods, in particular hybrid quantum/classical (QM/MM) approaches combined with free energy simulations. Here we combine the unique exploration power of new advanced sampling methods with density functional theory (DFT)-based QM/MM free energy methods for multiscale simulations of long-range protonation dynamics in biological systems. In this regard, we show that combining multiple walkers/well-tempered metadynamics with an extended-system adaptive biasing force method (MWE), provides a powerful approach for exploration of water-mediated proton transfer reactions in complex biochemical systems. We compare and combine the MWE method also with QM/MM-umbrella sampling and explore the sampling of the free energy landscape with both geometric (linear combination of proton transfer distances) and physical (center of excess charge) reaction coordinates, and show how these affect the convergence of the potential of mean force (PMF) and the activation free energy. We find that the QM/MM-MWE method can efficiently explore both direct and water-mediated proton transfer pathways together with forward and reverse hole transfer mechanisms in the highly complex proton channel of respiratory Complex I, while the QM/MM-US approach shows a systematic convergence of selected long-range proton transfer pathways. In this regard, we show that the PMF along multiple proton transfer pathways is recovered by combining the strengths of both approaches in a QM/MM-MWE/focused US (FUS) scheme, and revealing new mechanistic insight into the proton transfer principles of Complex I. Our findings provide a promising basis for the quantitative multi-scale simulations of long-range proton transfer reactions in biological systems.

Improved Sampling of Adaptive Path Collective Variables by Stabilized Extended-System Dynamics

Overcoming Sampling Issues and Improving Computational Efficiency in Collective-Variable-Based Enhanced-Sampling Simulations: A Tutorial

Frequency Adaptive Metadynamics for the Calculation of Rare-Event Kinetics.

Integrating Path Sampling with Enhanced Sampling for Rare-event Kinetics

Analyzing Complex Reaction Mechanisms Using Path Sampling

Metadynamics to Enhance Sampling in Biomolecular Simulations

Kinetics-Optimized Enhanced Sampling Using Mean First Passage Times

Enhanced Sampling with Sub-optimal Collective Variables: Reconciling Accuracy and Convergence Speed

Enhanced sampling of chemical and biochemical reactions with metadynamics

Enhanced sampling techniques in molecular dynamics simulations of biological systems

Effective Data-Driven Collective Variables for Free Energy Calculations from Metadynamics of Paths

Bayesian Parameter Estimation for Biochemical Reaction Networks Using Region-Based Adaptive Parallel Tempering

From Enhanced Sampling to Reaction Profiles

Enhanced Sampling for Free Energy Profiles with Post-Transition-State Bifurcations

Ab initio molecular dynamics with enhanced sampling in heterogeneous catalysis

Discover, Sample, and Refine: Exploring Chemistry with Enhanced Sampling Techniques

Efficient Sampling over Rough Energy Landscapes with High Barriers: A Combination of Metadynamics with Integrated Tempering Sampling.

Sampling free energy surfaces as slices by combining umbrella sampling and metadynamics

QM/MM Free Energy Calculations of Long-Range Biological Protonation Dynamics by Adaptive and Focused Sampling

Enhanced Sampling in Molecular Dynamics Using Metadynamics, Replica-Exchange, and Temperature-Acceleration

Combining transition path sampling with data-driven collective variables through a reactivity-biased shooting algorithm