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.

Using A Classical Potential As an Efficient Importance Function for Sampling from an Ab Initio Potential

Ab Initio Validation on the Connection Between Atomistic and Hydrodynamic Description to Unravel the Ion Dynamics of Warm Dense Matter

Probing the self-ionization of liquid water with ab initio deep potential molecular dynamics

An \emph{ab initio} method for locating characteristic potential energy minima of liquids

Sampling of Organic Solutes in Aqueous and Heterogeneous Environments Using Oscillating Excess Chemical Potentials in Grand Canonical-like Monte Carlo-Molecular Dynamics Simulations

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

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

Quasi-Chemical Theory with Cluster Sampling from Ab Initio Molecular Dynamics: Fluoride (F$^-$) Anion Hydration

Importance Sampling for Equilibrium Properties in Non-Adiabatic Systems: Application to Anomalous Interfacial Charge Transfer

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

An Efficient Umbrella Potential for the Accurate Calculation of Free Energies by Molecular Simulation.

A Quantum-Mechanics Molecular-Mechanics scheme for extended systems

On the use of QM/MM Surface Hopping simulations to understand thermally-activated rare event nonadiabatic transitions in the condensed phase

Redox potentials and acidity constants from density functional theory based molecular dynamics.

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

Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based $Ab$ $Initio$ Molecular Dynamics I: Theory, Algorithm, and Performance

Iterative derivation of effective potentials to sample the conformational space of proteins at atomistic scale

The solvation and hydrophobic interaction of non-polar molecules in water in the approximation of interatomic potentials: The Monte Carlo method

Achieving an Order of Magnitude Speed-up in Hybrid Functional and Plane Wave based Ab Initio Molecular Dynamics: Applications to Proton Transfer Reactions in Enzymes and in Solution

Ab Initio Potential Energy Surface for the Reactions Between H2O and H

Accelerating Ab Initio Quantum Mechanical and Molecular Mechanical (QM/MM) Molecular Dynamics Simulations with Multiple Time Step Integration and a Recalibrated Semiempirical QM/MM Hamiltonian.