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.

High-Performance Molecular Dynamics Simulations for Native Mass Spectrometry of Large Protein Complexes with the Fast Multipole Method

Quantum Fragment Based ab Initio Molecular Dynamics for Proteins.

A Force Balanced Fragmentation Method for Ab Initio Molecular Dynamic Simulation of Protein

Proton Dynamics in Protein Mass Spectrometry

Studying functional dynamics in bio-molecules using accelerated molecular dynamics

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

Constant-pH Molecular Dynamics Simulations for Large Biomolecular Systems.

Molecular Dynamics Simulations of Peptides and Proteins with Amplified Collective Motions.

Combining Fourier Transform Ion Mobility with Charge Detection Mass Spectrometry for the Analysis of Multimeric Protein Complexes

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.

NNP/MM: Accelerating Molecular Dynamics Simulations with Machine Learning Potentials and Molecular Mechanics

Scalable Indirect Free Energy Method Applied to Divalent Cation-Metalloprotein Binding

Constant pH Simulation with FMM Electrostatics in GROMACS. (A) Design and Applications

Native mass spectrometry for understanding dynamic protein complex

Molecular dynamics-based approaches for enhanced sampling of long-time, large-scale conformational changes in biomolecules

Running and analyzing massively parallel molecular simulations

NNP/MM: Accelerating molecular dynamics simulations with machine learning potentials and molecular mechanic

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

Protein Simulation Using Coarse-Grained Two-Bead Multipole Force Field with Polarizable Water Models.

Breaking the Million-Electron and 1 EFLOP/s Barriers: Biomolecular-Scale Ab Initio Molecular Dynamics Using MP2 Potentials

An Integrated Native Mass Spectrometry and Top-Down Proteomics Method That Connects Sequence to Structure and Function of Macromolecular Complexes