Equilibrium and kinetic selectivity profiling on the human adenosine receptors

The translational value of ligand-receptor binding kinetics in drug discovery

The translation of in vitro potency of a candidate drug, as determined by traditional pharmacology metrics (such as EC50 /IC50 and KD /Ki values), to in vivo efficacy and safety is challenging. Residence time, which represents the duration of drug-target interaction, can be part of a more comprehensive understanding of the dynamic nature of drug-target interactions in vivo, thereby enabling better prediction of drug efficacy and safety. As a consequence, a prolonged residence time may help in achieving sustained pharmacological activity, while transient interactions with shorter residence times may be favourable for targets associated with side effects. Therefore, integration of residence time into the early stages of drug discovery and development has yielded a number of clinical candidates with promising in vivo efficacy and safety profiles. Insights from residence time research thus contribute to the translation of in vitro potency to in vivo efficacy and safety. Further research and advances in measuring and optimizing residence time will bring a much-needed addition to the drug discovery process and the development of safer and more effective drugs. In this review, we summarize recent research progress on residence time, highlighting its importance from a translational perspective.

Unbinding Kinetics of Muscarinic M3 Receptor Antagonists Explained by Metadynamics Simulations

The residence time (RT), the time for which a drug remains bound to its biological target, is a critical parameter for drug design. The prediction of this key kinetic property has been proven to be challenging and computationally demanding in the framework of atomistic simulations. In the present work, we setup and applied two distinct metadynamics protocols to estimate the RTs of muscarinic M3 receptor antagonists. In the first method, derived from the conformational flooding approach, the kinetics of unbinding is retrieved from a physics-based parameter known as the acceleration factor α (i.e., the running average over time of the potential deposited in the bound state). Such an approach is expected to recover the absolute RT value for a compound of interest. In the second method, known as the t_META-D approach, a qualitative estimation of the RT is given by the time of simulation required to drive the ligand from the binding site to the solvent bulk. This approach has been developed to reproduce the change of experimental RTs for compounds targeting the same target. Our analysis shows that both computational protocols are able to rank compounds in agreement with their experimental RTs. Quantitative structure-kinetics relationship (SKR) models can be identified and employed to predict the impact of a chemical modification on the experimental RT once a calibration study has been performed.

A guide to enzyme kinetics in early drug discovery

Drugs interact with their target of interest to bring about the desired phenotypic outcome that results in disease alleviation. Traditionally, most lead optimization exercises were driven by affinity measures (like IC₅₀ ) to inform structure-activity relationship (SAR)-guided medicinal chemistry. However, an IC₅₀ value is a thermodynamic estimate measured under equilibrium conditions that can vary as a function of substrate concentration and/or time (the latter especially for nonequilibrium modalities). Further, like other thermodynamic estimates, it is a state-function that is indifferent to the path traversed from the initial state to the final state. This can be a cause for concern in drug discovery given the predominance of nonequilibrium interactions and the open thermodynamic nature of the human system. Under such situations, employing rates along with equilibrium constants (or IC₅₀ values) would be far more relevant to capture the time evolution of the small molecule's interaction with the target of interest. These rates are generally typified by the rate of association, rate of dissociation and the residence time of the small molecule on the target (target occupancy). These parameters, when combined with the concept of target vulnerability, therapeutic window, pharmacokinetic profile of the small molecule, estimates of endogenous ligand and target turnover, will shed critical insights into the kinetics and dynamics of a small molecule's interaction with the protein, and allow realistic modelling of the system to enable optimizations and dosing decisions. With that aim, this guide will attempt to introduce the traditional role of mechanistic enzymology within drug discovery and emphasize the importance of kinetics in guiding SAR-based optimizations. It will also present initial ideas on how kinetic investigation should be positioned relative to the temporal span of a drug-discovery pipeline to leverage maximal utility from the investment in time and effort.

Predicting Residence Time of GPCR Ligands with Machine Learning

Drug-target residence time, the duration of binding at a given protein target, has been shown in some protein families to be more significant for conferring efficacy than binding affinity. To carry out efficient optimization of residence time in drug discovery, machine learning models that can predict that value need to be developed. One of the main challenges with predicting residence time is the paucity of data. This chapter outlines all of the currently available ligand kinetic data, providing a repository that contains the largest publicly available source of GPCR-ligand kinetic data to date. To help decipher the features of kinetic data that might be beneficial to include in computational models for the prediction of residence time, the experimental evidence for properties that influence residence time are summarized. Finally, two different workflows for predicting residence time with machine learning are outlined. The first is a single-target model trained on ligand features; the second is a multi-target model trained on features generated from molecular dynamics simulations.

Hit-to-lead and lead optimization binding free energy calculations for G protein-coupled receptors

We apply the hit-to-lead ESMACS (enhanced sampling of molecular dynamics with approximation of continuum solvent) and lead-optimization TIES (thermodynamic integration with enhanced sampling) methods to compute the binding free energies of a series of ligands at the A1 and A2A adenosine receptors, members of a subclass of the GPCR (G protein-coupled receptor) superfamily. Our predicted binding free energies, calculated using ESMACS, show a good correlation with previously reported experimental values of the ligands studied. Relative binding free energies, calculated using TIES, accurately predict experimentally determined values within a mean absolute error of approximately 1 kcal mol-1. Our methodology may be applied widely within the GPCR superfamily and to other small molecule-receptor protein systems.

[Ligand-receptor binding kinetics in drug design]

Traditionally, the thermodynamic values of affinity are considered as the main criterion for the development of new drugs. Usually, these values for drugs are measured <i>in vitro</i> at steady concentrations of the receptor and ligand, which are differed from <i>in vivo</i> environment. Recent studies have shown that the kinetics of the process of drug binding to its receptor make significant contribution in the drug effectiveness. This has increased attention in characterizing and predicting the rate constants of association and dissociation of the receptor ligand at the stage of preclinical studies of drug candidates. A drug with a long residence time can determine ligand-receptor selectivity (kinetic selectivity), maintain pharmacological activity of the drug at its low concentration in vivo. The paper discusses the theoretical basis of protein-ligand binding, molecular determinants that control the kinetics of the drug-receptor binding. Understanding the molecular features underlying the kinetics of receptor-ligand binding will contribute to the rational design of drugs with desired properties.

Ensemble-Based Steered Molecular Dynamics Predicts Relative Residence Time of A2A Receptor Binders

Drug-target residence time, the length of time for which a small molecule stays bound to its receptor target, has increasingly become a key property for optimization in drug discovery programs. However, its in silico prediction has proven difficult. Here we describe a method, using atomistic ensemble-based steered molecular dynamics (SMD), to observe the dissociation of ligands from their target G protein-coupled receptor in a time scale suitable for drug discovery. These dissociation simulations accurately, precisely, and reproducibly identify ligand-residue interactions and quantify the change in ligand energy values for both protein and water. The method has been applied to 17 ligands of the A_2A adenosine receptor, all with published experimental kinetic binding data. The residues that interact with the ligand as it dissociates are known experimentally to have an effect on binding affinities and residence times. There is a good correlation ( R² = 0.79) between the computationally calculated change in water-ligand interaction energy and experimentally determined residence time. Our results indicate that ensemble-based SMD is a rapid, novel, and accurate semi-empirical method for the determination of drug-target relative residence time.

Homogeneous, Real-Time NanoBRET Binding Assays for the Histamine H3 and H4 Receptors on Living Cells

Receptor-binding affinity and ligand-receptor residence time are key parameters for the selection of drug candidates and are routinely determined using radioligand competition-binding assays. Recently, a novel bioluminescence resonance energy transfer (BRET) method utilizing a NanoLuc-fused receptor was introduced to detect fluorescent ligand binding. Moreover, this NanoBRET method gives the opportunity to follow fluorescent ligand binding on intact cells in real time, and therefore, results might better reflect in vivo conditions as compared with the routinely used cell homogenates or purified membrane fractions. In this study, a real-time NanoBRET-based binding assay was established and validated to detect binding of unlabeled ligands to the histamine H₃ receptor (H₃R) and histamine H₄ receptor on intact cells. Obtained residence times of clinically tested H₃R antagonists were reflected by their duration of H₃R antagonism in a functional receptor recovery assay.

Target and Tissue Selectivity Prediction by Integrated Mechanistic Pharmacokinetic-Target Binding and Quantitative Structure Activity Modeling

Selectivity is an important attribute of effective and safe drugs, and prediction of in vivo target and tissue selectivity would likely improve drug development success rates. However, a lack of understanding of the underlying (pharmacological) mechanisms and availability of directly applicable predictive methods complicates the prediction of selectivity. We explore the value of combining physiologically based pharmacokinetic (PBPK) modeling with quantitative structure-activity relationship (QSAR) modeling to predict the influence of the target dissociation constant (K _D) and the target dissociation rate constant on target and tissue selectivity. The K _D values of CB1 ligands in the ChEMBL database are predicted by QSAR random forest (RF) modeling for the CB1 receptor and known off-targets (TRPV1, mGlu5, 5-HT1a). Of these CB1 ligands, rimonabant, CP-55940, and Δ⁸-tetrahydrocanabinol, one of the active ingredients of cannabis, were selected for simulations of target occupancy for CB1, TRPV1, mGlu5, and 5-HT1a in three brain regions, to illustrate the principles of the combined PBPK-QSAR modeling. Our combined PBPK and target binding modeling demonstrated that the optimal values of the K _D and k _off for target and tissue selectivity were dependent on target concentration and tissue distribution kinetics. Interestingly, if the target concentration is high and the perfusion of the target site is low, the optimal K _D value is often not the lowest K _D value, suggesting that optimization towards high drug-target affinity can decrease the benefit-risk ratio. The presented integrative structure-pharmacokinetic-pharmacodynamic modeling provides an improved understanding of tissue and target selectivity.

From receptor binding kinetics to signal transduction; a missing link in predicting in vivo drug-action

An important question in drug discovery is how to overcome the significant challenge of high drug attrition rates due to lack of efficacy and safety. A missing link in the understanding of determinants for drug efficacy is the relation between drug-target binding kinetics and signal transduction, particularly in the physiological context of (multiple) endogenous ligands. We hypothesized that the kinetic binding parameters of both drug and endogenous ligand play a crucial role in determining cellular responses, using the NK1 receptor as a model system. We demonstrated that the binding kinetics of both antagonists (DFA and aprepitant) and endogenous agonists (NKA and SP) have significantly different effects on signal transduction profiles, i.e. potency values, in vitro efficacy values and onset rate of signal transduction. The antagonistic effects were most efficacious with slowly dissociating aprepitant and slowly associating NKA while the combination of rapidly dissociating DFA and rapidly associating SP had less significant effects on the signal transduction profiles. These results were consistent throughout different kinetic assays and cellular backgrounds. We conclude that knowledge of the relationship between in vitro drug-target binding kinetics and cellular responses is important to ultimately improve the understanding of drug efficacy in vivo.

A Pluridimensional View of Biased Agonism

When studying G protein-coupled receptor (GPCR) signaling and ligand-biased agonism, at least three dimensional spaces must be considered, as follows: 1) the distinct conformations that can be stabilized by different ligands promoting the engagement of different signaling effectors and accessory regulators; 2) the distinct subcellular trafficking that can be conferred by different ligands, which results in spatially distinct signals; and 3) the differential binding kinetics that maintain the receptor in specific conformation and/or subcellular localization for different periods of time, allowing for the engagement of distinct signaling effector subsets. These three pluridimensional aspects of signaling contribute to different faces of functional selectivity and provide a complex, interconnected way to define the signaling profile of each individual ligand acting at GPCRs. In this review, we discuss how each of these aspects may contribute to the diversity of signaling, but also how they shed light on the complexity of data analyses and interpretation. The impact of phenotype variability as a source of signaling diversity, and the influence of novel and more sensitive assays in the detection and analysis of signaling pluridimensionality, is also discussed. Finally, we discuss perspectives for the use of the concept of pluridimensional signaling in drug discovery, in which we highlight future predictive tools that may facilitate the identification of compounds with optimal therapeutic and safety properties based on the signaling signatures of drug candidates.

The Added Value of Assessing Ligand-Receptor Binding Kinetics in Drug Discovery

In the past decade drug research community has started to appreciate the indispensable role of ligand-receptor binding kinetics (BK) in drug discovery. Next to the classical equilibrium-based drug evaluation process with affinity and potency values as outcomes, kinetic investigation of the ligand-receptor interaction can aid compound triage in the hit-to-lead campaign and provide additional information to understand the molecular mechanism of drug action. Translational models incorporating BK are emerging as well, which represent powerful tools for the prediction of in vivo effects. In this viewpoint we will summarize some recent findings and discuss and emphasize the added value of ligand-receptor binding kinetics in drug research.