How robust is the ligand binding transition state?

14 August 2023, Version 1
This content is a preprint and has not undergone peer review at the time of posting.

Abstract

For many drug targets, it has been shown that the kinetics of drug binding (e.g. on rate and off rate) is more predictive of drug efficacy than thermodynamic quantities alone. This motivates the development of predictive computational models that can be used to optimize compounds on the basis of their kinetics. The structural details underpinning these computational models are not only found in the bound state, but also in the short-lived ligand binding transition state: the highest free energy point along the (un)binding pathway. Although this transition state cannot be directly observed experimentally, due to its extremely short lifetime, recent successes have demonstrated that modeling of the ligand binding transition state is possible with the help of enhanced sampling methods for molecular dynamics. In our previous work we analyzed the transition state ensemble for an inhibitor of soluble epoxide hydrolase (sEH) with a residence time of 11 minutes. Here we computationally modeled unbinding events for five additional inhibitors of sEH with residence times ranging from 14.25 to 31.75 minutes, with our results recapitulating these experimental residence times to within an order of magnitude. The unbinding ensembles are analyzed in detail, focusing on features of the ligand binding transition state. We find that ligands with similar structures and similar bound poses can show significant differences in their ligand binding transition states, in terms of their spatial location and their interactions with specific protein residues. However we also find similarities across the transition state ensembles when examining more general features such as ligand degrees of freedom. Together these findings show significant challenges for rational, kinetics-based drug design.

Keywords

binding kinetics
molecular dynamics
protein-ligand binding
enhanced sampling
weighted ensemble
drug design

Supplementary materials

Title
Description
Actions
Title
Supplemental Information
Description
A set of Figures (S1 - S12) and Tables (S1 - S3).
Actions
Title
Supplemental Methods
Description
Additional information describing the experimental procedures for determining transition rates.
Actions

Comments

Comments are not moderated before they are posted, but they can be removed by the site moderators if they are found to be in contravention of our Commenting Policy [opens in a new tab] - please read this policy before you post. Comments should be used for scholarly discussion of the content in question. You can find more information about how to use the commenting feature here [opens in a new tab] .
This site is protected by reCAPTCHA and the Google Privacy Policy [opens in a new tab] and Terms of Service [opens in a new tab] apply.