Abstract:The behaviour of biomolecular systems is governed by their thermodynamic and kinetic properties. It is thus important to be able to calculate, for example, both the affinity and rate of binding and dissociation of a protein-ligand complex, or the populations and exchange rates between distinct conformational states. Because these are typically rare events, calculating these properties from long molecular dynamics simulations remains extremely difficult. Instead, one often adopts a divide-and-conquer strategy in which equilibrium free-energy differences and the fastest state-to-state transition (e.g. ligand association or minor-to-major state conversion) are combined to estimate the slow rate (e.g. ligand dissociation) using a two-state assumption. Here we instead address these problems by using a previously developed method to calculate both the forward and backward rates directly from simulations. We then estimate the thermodynamics from the rates, and validate these values by independent means. We applied the approach to three systems of increasing complexity, including the association and dissociation of benzene to a fully buried cavity inside the L99A mutant variant of T4 lysozyme. In particular, we were able to determine both millisecond association and dissociation rates, and the affinity, of the protein-ligand system by directly observing dozens of rare events in atomic detail. Our approach both sheds light on the precision of methods for calculating kinetics and further provides a generally useful test for the internal consistency of kinetics and thermodynamics. We also expect our route to be useful for obtaining both the kinetics and thermodynamics at the same time in more challenging cases.

Multisecond ligand dissociation dynamics from atomistic simulations

Biomolecular conformational changes and ligand binding: from kinetics to thermodynamics.

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

Monitoring receptor-ligand interactions between surfaces by thermal fluctuations.

Mapping transiently formed and sparsely populated conformations on a complex energy landscape

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

Effects of electrostatic interactions on ligand dissociation kinetics

Ligand Binding and Release Investigated by Contact-Guided Iterative Multiple Independent Molecular Dynamics Simulations

Investigation on Structure and Dynamics of Heterodimeric SAM:SAM Complexes Using Micro-Second Molecular Dynamics Simulations

Dynamic coarse-graining fills the gap between atomistic simulations and experimental investigations of mechanical unfolding

A Workflow for Exploring Ligand Dissociation from a Macromolecule: Efficient Random Acceleration Molecular Dynamics Simulation and Interaction Fingerprints Analysis of Ligand Trajectories

Binding free energy of protein/ligand complexes calculated using dissociation Parallel Cascade Selection Molecular Dynamics and Markov state model

Coarse-grained dynamics of transiently-bound fast linkers

Enhanced-Sampling Simulations for the Estimation of Ligand Binding Kinetics: Current Status and Perspective

Prediction of Threonine-Tyrosine Kinase Receptor-Ligand Unbinding Kinetics with Multiscale Milestoning and Metadynamics

Identifying ligand binding sites and poses using GPU-accelerated Hamiltonian replica exchange molecular dynamics

Kinetics and Free Energy of Ligand Dissociation Using Weighted Ensemble Milestoning

Heterogeneous and rate-dependent streptavidin-biotin unbinding revealed by high-speed force spectroscopy and atomistic simulations

Ligand Gaussian accelerated Molecular Dynamics 3 (LiGaMD3): Improved Calculations of Binding Thermodynamics and Kinetics of Both Small Molecules and Flexible Peptides

MSM/RD: Coupling Markov state models of molecular kinetics with reaction-diffusion simulations

Anomalous relaxation kinetics of biological lattice-ligand binding models