Abstract:Understanding the interaction of ligands with biomolecules is an integral component of drug discovery and development. Challenges for computing thermodynamic and kinetic quantities for pharmaceutically relevant receptor-ligand complexes include the size and flexibility of the ligands, large-scale conformational rearrangements of the receptor, accurate force field parameters, simulation efficiency, and sufficient sampling associated with rare events. Our recently developed multiscale milestoning simulation approach, SEEKR2 (Simulation Enabled Estimation of Kinetic Rates v.2), has demonstrated success in predicting unbinding (koff) kinetics by employing molecular dynamics (MD) simulations in regions closer to the binding site. The MD region is further subdivided into smaller Voronoi tessellations to improve the simulation efficiency and parallelization. To date, all MD simulations are run using general molecular mechanics (MM) force fields. The accuracy of calculations can be further improved by incorporating quantum mechanical (QM) methods into generating system-specific force fields through reparameterizing ligand partial charges in the bound state. The force field reparametrization process modifies the potential energy landscape of the bimolecular complex, enabling a more accurate representation of the intermolecular interactions and polarization effects at the bound state. We present QMrebind (Quantum Mechanical force field reparameterization at the receptor-ligand binding site), an ORCA-based software that facilitates reparameterizing the potential energy function within the phase space representing the bound state in a receptor-ligand complex. With SEEKR2 koff estimates and experimentally determined kinetic rates, we compare and interpret the receptor-ligand unbinding kinetics obtained using the newly reparameterized force fields for model host-guest systems and HSP90-inhibitor complexes. This method provides an opportunity to achieve higher accuracy in predicting receptor-ligand koff rate constants.

Microsecond timescale MD simulations at the transition state of Pm HMGR predict remote allosteric residues

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.

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

Ten-Microsecond Molecular Dynamics Simulation of a Fast-Folding WW Domain

pH-dependent reaction triggering in PmHMGR crystals for time-resolved crystallography

Experimental Confirmation Of An Ntrc Transition Pathway Predicted By Targeted Md

Exploration of Multi-State Conformational Dynamics and Underlying Global Functional Landscape of Maltose Binding Protein

Protein-peptide association kinetics beyond the seconds timescale from atomistic simulations

QMrebind: Incorporating quantum mechanical force field reparameterization at the ligand binding site for improved drug-target kinetics through milestoning simulations

An application of QM/MM simulation: the second protonation of cytochrome P450.

Quantum molecular dynamics simulation of proton transfer in cytochrome c oxidase

Multisecond ligand dissociation dynamics from atomistic simulations

Millisecond Allosteric Dynamics Of Activated Ras Reproduced With A Slowly Hydrolyzable Gtp Analogue

Studying functional dynamics in bio-molecules using accelerated molecular dynamics

Simulation of Spontaneous G Protein Activation Reveals a New Intermediate Driving GDP Unbinding

Hidden GPCR structural transitions addressed by multiple walker supervised molecular dynamics (mwSuMD)

Accurate Machine Learned Quantum-Mechanical Force Fields for Biomolecular Simulations

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

Enhanced Sampling in Molecular Dynamics Simulations: How Many MD Snapshots can be Needed to Reproduce the Biological Behavior?

Ligand Binding Rate Constants in Heme Proteins Using Markov State Models and Molecular Dynamics Simulations

An Evaluation of Force Field Accuracy for the Mini-Protein Chignolin using Markov State Models