Prediction of Drug-Target Binding Kinetics for Flexible Proteins by Comparative Binding Energy Analysis

04 June 2021, Version 2
This content is a preprint and has not undergone peer review at the time of posting.


There is growing consensus that the optimization of the kinetic parameters for drug-protein binding leads to improved drug efficacy. Therefore, computational methods have been developed to predict kinetic rates and to derive quantitative structure-kinetic relationships (QSKRs). Many of these methods are based on crystal structures of ligand-protein complexes. However, a drawback is that each protein-ligand complex is usually treated as having a single structure. Here, we present a modification of COMparative BINding Energy (COMBINE) analysis, which uses the structures of protein-
ligand complexes to predict binding parameters. We introduce the option to use multiple structures to describe each ligand-protein complex into COMBINE analysis and
apply this to study the effects of protein flexibility on the derivation of dissociation rate constants (koff) for inhibitors of p38 mitogen-activated protein (MAP) kinase, which has a flexible binding site. Multiple structures were obtained for each ligand-protein complex by performing docking to an ensemble of protein configurations obtained from molecular dynamics simulations. Coefficients to scale ligand-protein interaction energies determined from energy-minimized structures of ligand-protein complexes were obtained by partial least squares regression and allowed the computation of koff values. The QSKR model obtained using single, energy minimized crystal structures for each ligand-protein complex had a higher predictive power than the QSKR model obtained with multiple structures from ensemble docking. However, the incorporation of protein-ligand flexibility helped to highlight additional ligand-protein interactions that lead to longer residence times, like interactions with residues Arg67 and Asp168, which are close to the ligand in many crystal structures. These results show that COMBINE analysis is a promising method to guide the design of compounds that bind to flexible proteins with improved binding kinetics.


ligand binding kinetics
drug design


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.