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.

A statistical model for predicting thermal chemical reaction rate: Application to bimolecule reactions

A statistical model for predicting thermal chemical reaction rate

Multi-Scale Coupling Model Involving Bed and Catalyst Particle for MTP Reactor Simulation

Utility computable modeling of a Boltzmann model equation for bimolecular chemical reactions and numerical application

Molecular dynamics simulation of a chemical reaction in solution

Bimolecular Nucleophilic Substitution Reactions: Predictive Models for Rate Constants and Molecular Reaction Pairs Analysis

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

Simple statistical model for predicting thermal atom diffusion on crystal surfaces

Correction of Reaction Models Using Collision Limit Violation Analyses: Application to a Silane Reaction Model

Simulating Bistable Biochemical Systems by Means of Reactive Multiparticle Collision Dynamics

On Kinetic Phase Transition of the Dimer-Monomer Reaction Model

Thermodynamic model for electron emission and negative- and positive-ion formation in keV molecular collisions

A Comparison of Bimolecular Reaction Models for Stochastic Reaction Diffusion Systems

Accuracy of reaction coordinate based rate theories for modelling chemical reactions: insights from the thermal isomerization in retinal

A generalized interface reaction kinetic model for describing heterogeneous processes driven by contracting mechanisms

Theoretical Study of the Thermal Rate Coefficients of the H3+ + C2H4 Reaction: Dynamics Study on a Full-Dimensional Potential Energy Surface

Model-Free Kinetic Determination Of Pre-Exponential Factor And Reaction Mechanism In Accelerating Rate Calorimetry

A Discrete Stochastic Formulation for Reversible Bimolecular Reactions via Diffusion Encounter

A Generalized Kinetic Model for Heterogeneous Gas-Solid Reactions

Exactly Solvable Model of Monomer-Monomer Reactions on a Two-Dimensional Random Catalytic Substrate

Brownian dynamics simulation of a chemical reaction in solution