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.

Tutorial on Umbrella Sampling Simulation with a Combined QM/MM Potential: The Potential of Mean Force for an S N 2 Reaction in Water

Tutorial on Umbrella Sampling Simulation with a Combined QM/MM Potential: The Potential of Mean Force for an SN2 Reaction in Water

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

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

A combined quantum mechanical and molecular mechanical potential for molecular dynamics simulations

Combine Umbrella Sampling with Integrated Tempering Method for Efficient and Accurate Calculation of Free Energy Changes of Complex Energy Surface

Learning the Quantum Centroid Force Correction in Molecular Systems: A Localized Approach

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

A Quantum-Mechanics Molecular-Mechanics scheme for extended systems

A quantum mechanical computational method for modeling electrostatic and solvation effects of protein

Free energy calculations of a proton transfer reaction by simulated tempering umbrella sampling first principles molecular dynamics simulations

Study on the Application of the Combination of TMD Simulation and Umbrella Sampling in PMF Calculation for Molecular Conformational Transitions

Analytic gradients for the electrostatic embedding QM/MM in periodic boundary conditions using particle-mesh Ewald sums and electrostatic potential fitted charge operators

Assessment of free energies from electrostatic embedding density functional tight binding-based/molecular mechanics in periodic boundary conditions

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

Ab initio statistical mechanics of surface adsorption and desorption: II. Nuclear quantum effects

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.

Accelerated Quantum Mechanics/Molecular Mechanics Simulations Via Neural Networks Incorporated with Mechanical Embedding Scheme.

The adaptive buffered force QM/MM method in the CP2K and AMBER software packages

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