The role of water and steric constraints in the kinetics of cavity-ligand unbinding

Abstract

A key factor influencing a drug's efficacy is its residence time in the binding pocket of the host protein. Using atomistic computer simulation to predict this residence time and the associated dissociation process is a desirable but extremely difficult task due to the long timescales involved. This gets further complicated by the presence of biophysical factors such as steric and solvation effects. In this work, we perform molecular dynamics (MD) simulations of the unbinding of a popular prototypical hydrophobic cavity-ligand system using a metadynamics based approach that allows direct assessment of kinetic pathways and parameters. When constrained to move in an axial manner, we find the unbinding time to be on the order of 4000 sec. In accordance with previous studies, we find that the ligand must pass through a region of sharp dewetting transition manifested by sudden and high fluctuations in solvent density in the cavity. When we remove the steric constraints on ligand, the unbinding happens predominantly by an alternate pathway, where the unbinding becomes 20 times faster, and the sharp dewetting transition instead becomes continuous. We validate the unbinding timescales from metadynamics through a Poisson analysis, and by comparison through detailed balance to binding timescale estimates from unbiased MD. This work demonstrates that enhanced sampling can be used to perform explicit solvent molecular dynamics studies at timescales previously unattainable, obtaining direct and reliable pictures of the underlying physio-chemical factors including free energies and rate constants.

7 pages, 4 figures, supplementary PDF file, submitted

References in corpus (5)

• Canonical sampling through velocity-rescaling
• Well-Tempered Metadynamics: A Smoothly Converging and Tunable Free-Energy Method
• PLUMED: a portable plugin for free-energy calculations with molecular dynamics
• Pathways to dewetting in hydrophobic confinement
• Are hydrodynamic interactions important in the kinetics of hydrophobic collapse?

Cited by in corpus (16)

• Caliber based spectral gap optimization of order parameters (SGOOP) for sampling complex molecular systems
• Understanding Hydrophobic Effects: Insights from Water Density Fluctuations
• Frequency adaptive metadynamics for the calculation of rare-event kinetics
• Overcoming timescale and finite-size limitations to compute nucleation rates from small scale Well Tempered Metadynamics simulations
• Predicting reaction coordinates in energy landscapes with diffusion anisotropy
• How wet should be the reaction coordinate for ligand unbinding?
• Assessing models of force-dependent unbinding rates via infrequent metadynamics
• Variational implicit-solvent predictions of the dry-wet transition pathways for ligand-receptor binding and unbinding kinetics
• Kramers turnover: from energy diffusion to spatial diffusion using metadynamics
• Solvent-mediated interactions between nanostructures: from water to Lennard-Jones liquid
• On the Role of Solvent in Hydrophobic Cavity-ligand Recognition Kinetics
• Solute particle near a nanopore: influence of size and surface properties on the solvent-mediated forces
• Thermodynamically Optimized Machine-learned Reaction Coordinates for Hydrophobic Ligand Dissociation
• Affinity, Kinetics, and Pathways of Anisotropic Ligands Binding to Hydrophobic Model Pockets
• First-Passage Approach to Optimizing Perturbations for Improved Training of Machine Learning Models
• Computing rates and understanding unbinding mechanism in host-guest systems