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.

Faster Sampling in Molecular Dynamics Simulations with TIP3P-F Water

Smoothed Molecular Dynamics for Large Step Time Integration

Improving the Efficiency of Non-equilibrium Sampling in the Aqueous Environment via Implicit-Solvent Simulations.

A local resampling trick for focused molecular dynamics

Editorial: Advanced Sampling and Modeling in Molecular Simulations for Slow and Large-Scale Biomolecular Dynamics

Path Integral Molecular Dynamics within the Grand Canonical-like Adaptive Resolution Technique: Simulation of Liquid Water

Modeling diffusive dynamics in adaptive resolution simulation of liquid water

Comparative Study of the Water Response to External Force at Nanoscale and Mesoscale

Exploring A Multi-Scale Method for Molecular Simulations in Continuum Solvent Model: Explicit Simulation of Continuum Solvent As An Incompressible Fluid

A Review of Enhanced Sampling Approaches for Accelerated Molecular Dynamics

Enhanced solvation force extrapolation for speeding up molecular dynamics simulations of complex biochemical liquids

Enhanced Sampling Method in Molecular Simulations

Accelerating the 3D-RISM theory of molecular solvation with treecode summation and cut-offs

Developing Semi-Empirical Water Model for Efficiently Simulating Temperature-Dependent Chemisorption of CO2 in Amine Solvents

Coupling Accelerated Molecular Dynamics Methods with Thermodynamic Integration Simulations

Accelerating ab initio QM/MM Molecular Dynamics Simulations with Multiple Time Step Integration and a Recalibrated Semi-empirical QM/MM Hamiltonian

Kinetics-Optimized Enhanced Sampling Using Mean First Passage Times

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

Speeding-up Hybrid Functional based Ab Initio Molecular Dynamics using Multiple Time-stepping and Resonance Free Thermostat

Accelerated Atomistic Simulations of a Supramolecular Polymer in Water

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