Structure-affinity and structure-residence time relationships of macrocyclic Gαq protein inhibitors

The macrocyclic depsipeptides YM-254890 (YM) and FR900359 (FR) are potent inhibitors of Gα_q/11 proteins. They are important pharmacological tools and have potential as therapeutic drugs. The hydrogenated, tritium-labeled YM and FR derivatives display largely different residence times despite similar structures. In the present study we established a competition-association binding assay to determine the dissociation kinetics of unlabeled G_q protein inhibitors. Structure-affinity and structure-residence time relationships were analyzed. Small structural modifications had a large impact on residence time. YM and FR exhibited 4- to 10-fold higher residence times than their hydrogenated derivatives. While FR showed pseudo-irreversible binding, YM displayed much faster dissociation from its target. The isopropyl anchor present in FR and some derivatives was essential for slow dissociation. These data provide a basis for future drug design toward modulating residence times of macrocyclic G_q protein inhibitors, which has been recognized as a crucial determinant for therapeutic outcome.

Competition in drug binding and … the race to equilibrium

Binding kinetics has become a popular topic in pharmacology due to its potential contribution to the selectivity and duration of drug action. Yet, the overall kinetic aspects of complex binding mechanisms are still merely described in terms of elaborate algebraic equations. Interestingly, it has been recommended some 10 years ago to examine such mechanisms in terms of binding fluxes instead of the conventional rate constants. Alike the velocity of product formation in enzymology, those fluxes refer to the velocity by which one target species converts into another one. Novel binding flux-based approaches are utilized to get a better visual insight into the "competition" between two drugs/ligands for a single target as well as between induced fit- and conformational selection pathways for a single ligand within a thermodynamic cycle. The present data were obtained by differential equation-based simulations. Early on, the ligand-binding steps "race" to equilibrium (i.e., when their forward and reverse fluxes are equal) at their individual pace. The overall/global equilibrium is only reached later on. For the competition association assays, this parting might produce a transient "overshoot" of one of the bound target species. A similar overshoot may also show up within a thermodynamic cycle and, at first glance, suggest that the induced fit pathway dominates. Yet, present findings show that under certain circumstances, it could rather be the other way round. Novel binding flux-based approaches offer visually attractive insights into crucial aspects of "complex" binding mechanisms under non-equilibrium conditions.

Optical-Controlled Kinetic Switch: Fine-Tuning of the Residence Time of an Antagonist Binding to the Vasopressin V2 Receptor in In Vitro, Ex Vivo, and In Vivo Models of ADPKD

The pharmacological activity of a small-molecule ligand is linked to its receptor residence time. Therefore, precise control of the duration for which a ligand binds to its receptor is highly desirable. Herein, we designed photoswitchable ligands targeting the vasopressin V₂ receptor (V₂R), a validated target for autosomal dominant polycystic kidney disease (ADPKD). We adapted the photoswitching trait of azobenzene to the parent V₂R antagonist lixivaptan (LP) to generate azobenzene lixivaptan derivatives (aLPs). Among them, aLPs-5g was a potential optical-controlled kinetic switch. Upon irradiation, cis-aLPs-5g displayed a 4.3-fold prolonged V₂R residence time compared to its thermally stable trans configuration. The optical-controlled kinetic variations led to distinct inhibitory effects on cellular functional readout. Furthermore, conversion of the cis/trans isomer of aLPs-5g resulted in different efficacies of inhibiting renal cystogenesis ex vivo and in vivo. Overall, aLPs-5g represents a photoswitch for precise control of ligand-receptor residence time and, consequently, the pharmacological activity.

Dual A1/A3 Adenosine Receptor Antagonists: Binding Kinetics and Structure-Activity Relationship Studies Using Mutagenesis and Alchemical Binding Free Energy Calculations

Drugs targeting adenosine receptors (AR) can provide treatment for diseases. We report the identification of 7-(phenylamino)-pyrazolo[3,4-c]pyridines L2-L10, A15, and A17 as low-micromolar to low-nanomolar A₁R/A₃R dual antagonists, with 3-phenyl-5-cyano-7-(trimethoxyphenylamino)-pyrazolo[3,4-c]pyridine (A17) displaying the highest affinity at both receptors with a long residence time of binding, as determined using a NanoBRET-based assay. Two binding orientations of A17 produce stable complexes inside the orthosteric binding area of A₁R in molecular dynamics (MD) simulations, and we selected the most plausible orientation based on the agreement with alanine mutagenesis supported by affinity experiments. Interestingly, for drug design purposes, the mutation of L250^6.51 to alanine increased the binding affinity of A17 at A₁R. We explored the structure-activity relationships against A₁R using alchemical binding free energy calculations with the thermodynamic integration coupled with the MD simulation (TI/MD) method, applied on the whole G-protein-coupled receptor-membrane system, which showed a good agreement (r = 0.73) between calculated and experimental relative binding free energies.

International Union of Basic and Clinical Pharmacology. CXII: Adenosine Receptors: A Further Update

Our previous International Union of Basic and Clinical Pharmacology report on the nomenclature and classification of adenosine receptors (2011) contained a number of emerging developments with respect to this G protein-coupled receptor subfamily, including protein structure, protein oligomerization, protein diversity, and allosteric modulation by small molecules. Since then, a wealth of new data and results has been added, allowing us to explore novel concepts such as target binding kinetics and biased signaling of adenosine receptors, to examine a multitude of receptor structures and novel ligands, to gauge new pharmacology, and to evaluate clinical trials with adenosine receptor ligands. This review should therefore be considered a further update of our previous reports from 2001 and 2011. SIGNIFICANCE STATEMENT: Adenosine receptors (ARs) are of continuing interest for future treatment of chronic and acute disease conditions, including inflammatory diseases, neurodegenerative afflictions, and cancer. The design of AR agonists ("biased" or not) and antagonists is largely structure based now, thanks to the tremendous progress in AR structural biology. The A2A- and A2BAR appear to modulate the immune response in tumor biology. Many clinical trials for this indication are ongoing, whereas an A2AAR antagonist (istradefylline) has been approved as an anti-Parkinson agent.

The Problems of Applying Classical Pharmacology Analysis to Modern In Vitro Drug Discovery Assays: Slow Binding Kinetics and High Target Concentration

The analysis framework used to quantify drug potency in vitro (e.g., K_d or K_i) was initially developed for classical pharmacology bioassays, for example, organ bath experiments testing moderate-affinity natural products. Modern drug discovery can infringe the assumptions of the classical pharmacology analysis equations, owing to the reduction of assay volume in miniaturization, target overexpression, and the increase of compound-target affinity in medicinal chemistry. These assumptions are that (1) the compound concentration greatly exceeds the target concentration (i.e., minimal ligand depletion), and (2) the compound is at equilibrium with the receptor (i.e., rapid ligand binding kinetics). Unappreciated infringement of these assumptions can lead to substantial underestimation of compound affinity, which negatively impacts the drug discovery process, from early-stage lead optimization to prediction of human dosing. This study evaluates the real-world impact of these factors on the target interaction assays used in drug discovery using literature examples, database searches, and simulations. The ranges of compound affinity and the assay types that are prone to depletion and equilibration artifacts are identified. Importantly, the highest-affinity compounds, usually the highest value chemical matter in drug discovery, are the most affected. Methods and simulation tools are provided to enable investigators to evaluate, manage, and minimize depletion or equilibration artifacts. This study enables the correct application of pharmacological data analysis to accurately quantify affinity using modern drug discovery assay technology.

Binding pathway determines norepinephrine selectivity for the human β1AR over β2AR

Beta adrenergic receptors (βARs) mediate physiologic responses to the catecholamines epinephrine and norepinephrine released by the sympathetic nervous system. While the hormone epinephrine binds β₁AR and β₂AR with similar affinity, the smaller neurotransmitter norepinephrine is approximately tenfold selective for the β₁AR. To understand the structural basis for this physiologically important selectivity, we solved the crystal structures of the human β₁AR bound to an antagonist carazolol and different agonists including norepinephrine, epinephrine and BI-167107. Structural comparison revealed that the catecholamine-binding pockets are identical between β₁AR and β₂AR, but the extracellular vestibules have different shapes and electrostatic properties. Metadynamics simulations and mutagenesis studies revealed that these differences influence the path norepinephrine takes to the orthosteric pocket and contribute to the different association rates and thus different affinities.

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.

Design and pharmacological profile of a novel covalent partial agonist for the adenosine A1 receptor

Partial agonists for G protein-coupled receptors (GPCRs) provide opportunities for novel pharmacotherapies with enhanced on-target safety compared to full agonists. For the human adenosine A₁ receptor (hA₁AR) this has led to the discovery of capadenoson, which has been in phase IIa clinical trials for heart failure. Accordingly, the design and profiling of novel hA₁AR partial agonists has become an important research focus. In this study, we report on LUF7746, a capadenoson derivative bearing an electrophilic fluorosulfonyl moiety, as an irreversibly binding hA₁AR modulator. Meanwhile, a nonreactive ligand bearing a methylsulfonyl moiety, LUF7747, was designed as a control probe in our study. In a radioligand binding assay, LUF7746's apparent affinity increased to nanomolar range with longer pre-incubation time, suggesting an increasing level of covalent binding over time. Moreover, compared to the reference full agonist CPA, LUF7746 was a partial agonist in a hA₁AR-mediated G protein activation assay and resistant to blockade with an antagonist/inverse agonist. An in silico structure-based docking study combined with site-directed mutagenesis of the hA₁AR demonstrated that amino acid Y271^7.36 was the primary anchor point for the covalent interaction. Additionally, a label-free whole-cell assay was set up to identify LUF7746's irreversible activation of an A₁ receptor-mediated cell morphological response. These results led us to conclude that LUF7746 is a novel covalent hA₁AR partial agonist and a valuable chemical probe for further mapping the receptor activation process. It may also serve as a prototype for a therapeutic approach in which a covalent partial agonist may cause less on-target side effects, conferring enhanced safety compared to a full agonist.

In Silico Prediction of the Dissociation Rate Constants of Small Chemical Ligands by 3D-Grid-Based VolSurf Method

Accumulated evidence suggests that binding kinetic properties—especially dissociation rate constant or drug-target residence time—are crucial factors affecting drug potency. However, quantitative prediction of kinetic properties has always been a challenging task in drug discovery. In this study, the VolSurf method was successfully applied to quantitatively predict the k_off values of the small ligands of heat shock protein 90α (HSP90α), adenosine receptor (AR) and p38 mitogen-activated protein kinase (p38 MAPK). The results showed that few VolSurf descriptors can efficiently capture the key ligand surface properties related to dissociation rate; the resulting models demonstrated to be extremely simple, robust and predictive in comparison with available prediction methods. Therefore, it can be concluded that the VolSurf-based prediction method can be widely applied in the ligand-receptor binding kinetics and de novo drug design researches.

Perspective: Implications of Ligand-Receptor Binding Kinetics for Therapeutic Targeting of G Protein-Coupled Receptors

The concept of ligand-receptor binding kinetics has been broadly applied in drug development pipelines focusing on G protein-coupled receptors (GPCRs). The ligand residence time (RT) for a receptor describes how long a ligand-receptor complex exists, and is defined as the reciprocal of the dissociation rate constant (k _off). RT has turned out to be a valuable parameter for GPCR researchers focusing on drug development as a good predictor of in vivo efficacy. The positive correlation between RT and in vivo efficacy has been established for several drugs targeting class A GPCRs (e.g., the neurokinin-1 receptor (NK₁R), the β₂ adrenergic receptor (β₂AR), and the muscarinic 3 receptor (M₃R)) and for drugs targeting class B1 (e.g., the glucagon-like peptide 1 receptor (GLP-1R)). Recently, the association rate constant (k _on) has gained similar attention as another parameter affecting in vivo efficacy. In the current perspective, we address the importance of studying ligand-receptor binding kinetics for therapeutic targeting of GPCRs, with an emphasis on how binding kinetics can be altered by subtle molecular changes in the ligands and/or the receptors and how such changes affect treatment outcome. Moreover, we speculate on the impact of binding kinetic parameters for functional selectivity and sustained receptor signaling from endosomal compartments; phenomena that have gained increasing interest in attempts to improve therapeutic targeting of GPCRs.

Considerations for improved performance of competition association assays analysed with the Motulsky-Mahan's "kinetics of competitive binding" model

Background and Purpose

Target engagement dynamics can influence drugs' pharmacological effects. Kinetic parameters for drug:target interactions are often quantified by evaluating competition association experiments—measuring simultaneous protein binding of labelled tracers and unlabelled test compounds over time—with Motulsky–Mahan's “kinetics of competitive binding” model. Despite recent technical improvements, the current assay formats impose practical limitations to this approach. This study aims at the characterisation, understanding and prevention of these experimental constraints, and associated analytical challenges.

Experimental Approach

Monte Carlo simulations were used to run virtual kinetic and equilibrium tracer binding and competition experiments in both normal and perturbed assay conditions. Data were fitted to standard equations derived from the mass action law (including Motulsky–Mahan's) and to extended versions aiming to cope with frequently observed deviations of the canonical traces. Results were compared to assess the precision and accuracy of these models and identify experimental factors influencing their performance.

Key Results

Key factors influencing the precision and accuracy of the Motulsky–Mahan model are the interplay between compound dissociation rates, measurement time and interval frequency, tracer concentration and binding kinetics and the relative abundance of equilibrium complexes in vehicle controls. Experimental results produced recommendations for better design of tracer characterisation experiments and new strategies to deal with systematic signal decay.

Conclusions and Implications

Our data advances our comprehension of the Motulsky–Mahan kinetics of competitive binding models and provides experimental design recommendations, data analysis tools, and general guidelines for its practical application to in vitro pharmacology and drug screening.

Route to Prolonged Residence Time at the Histamine H1 Receptor: Growing from Desloratadine to Rupatadine

Drug–target binding kinetics are an important predictor of in vivo drug efficacy, yet the relationship between ligand structures and their binding kinetics is often poorly understood. We show that both rupatadine (1) and desloratadine (2) have a long residence time at the histamine H₁ receptor (H₁R). Through development of a [³H]levocetirizine radiolabel, we find that the residence time of 1 exceeds that of 2 more than 10-fold. This was further explored with 22 synthesized rupatadine and desloratadine analogues. Methylene-linked cycloaliphatic or β-branched substitutions of desloratadine increase the residence time at the H₁R, conveying a longer duration of receptor antagonism. However, cycloaliphatic substituents directly attached to the piperidine amine (i.e., lacking the spacer) have decreased binding affinity and residence time compared to their methylene-linked structural analogues. Guided by docking studies, steric constraints within the binding pocket are hypothesized to explain the observed differences in affinity and binding kinetics between analogues.

Kinetics of Drug Binding and Residence Time

The kinetics of drug binding and unbinding is assuming an increasingly crucial role in the long, costly process of bringing a new medicine to patients. For example, the time a drug spends in contact with its biological target is known as residence time (the inverse of the kinetic constant of the drug-target unbinding, 1/ koff). Recent reports suggest that residence time could predict drug efficacy in vivo, perhaps even more effectively than conventional thermodynamic parameters (free energy, enthalpy, entropy). There are many experimental and computational methods for predicting drug-target residence time at an early stage of drug discovery programs. Here, we review and discuss the methodological approaches to estimating drug binding kinetics and residence time. We first introduce the theoretical background of drug binding kinetics from a physicochemical standpoint. We then analyze the recent literature in the field, starting from the experimental methodologies and applications thereof and moving to theoretical and computational approaches to the kinetics of drug binding and unbinding. We acknowledge the central role of molecular dynamics and related methods, which comprise a great number of the computational methods and applications reviewed here. However, we also consider kinetic Monte Carlo. We conclude with the outlook that drug (un)binding kinetics may soon become a go/no go step in the discovery and development of new medicines.

Probe dependency in the determination of ligand binding kinetics at a prototypical G protein-coupled receptor

Drug-target binding kinetics are suggested to be important parameters for the prediction of in vivo drug-efficacy. For G protein-coupled receptors (GPCRs), the binding kinetics of ligands are typically determined using association binding experiments in competition with radiolabelled probes, followed by analysis with the widely used competitive binding kinetics theory developed by Motulsky and Mahan. Despite this, the influence of the radioligand binding kinetics on the kinetic parameters derived for the ligands tested is often overlooked. To address this, binding rate constants for a series of histamine H₁ receptor (H₁R) antagonists were determined using radioligands with either slow (low k_off) or fast (high k_off) dissociation characteristics. A correlation was observed between the probe-specific datasets for the kinetic binding affinities, association rate constants and dissociation rate constants. However, the magnitude and accuracy of the binding rate constant-values was highly dependent on the used radioligand probe. Further analysis using recently developed fluorescent binding methods corroborates the finding that the Motulsky-Mahan methodology is limited by the employed assay conditions. The presented data suggest that kinetic parameters of GPCR ligands depend largely on the characteristics of the probe used and results should therefore be viewed within the experimental context and limitations of the applied methodology.

SNAPshots of the MCHR1: a Comparison Between the PET-Tracers [18F]FE@SNAP and [11C]SNAP-7941

PURPOSE

The melanin-concentrating hormone receptor 1 (MCHR1) has become an important pharmacological target, since it may be involved in various diseases, such as diabetes, insulin resistance, and obesity. Hence, a suitable positron emission tomography radiotracer for the in vivo assessment of the MCHR1 pharmacology is imperative. The current paper contrasts the extensive in vitro, in vivo, and ex vivo assessments of the radiotracers [18F]FE@SNAP and [11C]SNAP-7941 and provides comprehensive information about their biological and physicochemical properties. Furthermore, it examines their suitability for first-in-man imaging studies.

PROCEDURES

Kinetic real-time cell-binding studies with [18F]FE@SNAP and [11C]SNAP-7941 were conducted on adherent Chines hamster ovary (CHO-K1) cells stably expressing the human MCHR1 and MCHR2. Small animal imaging studies on mice and rats were performed under displacement and baseline conditions, as well as after pretreatment with the P-glycoprotein/breast cancer resistant protein inhibitor tariquidar. After the imaging studies, detailed analyses of the ex vivo biodistribution were performed. Ex vivo metabolism was determined in rat blood and brain and analyzed at various time points using a quantitative radio-HPLC assay.

RESULTS

[11C]SNAP-7941 demonstrates high uptake on CHO-K1-hMCHR1 cells, whereas no uptake was detected for the CHO-K1-hMCHR2 cells. In contrast, [18F]FE@SNAP evinced binding to CHO-K1-hMCHR1 and CHO-K1-hMCHR2 cells. Imaging studies with [18F]FE@SNAP and [11C]SNAP-7941 showed an increased brain uptake after tariquidar pretreatment in mice, as well as in rats, and exhibited a significant difference between the time-activity curves of the baseline and blocking groups. Biodistribution of both tracers demonstrated a decreased uptake after displacement. [11C]SNAP-7941 revealed a high metabolic stability in rats, whereas [18F]FE@SNAP was rapidly metabolized.

CONCLUSIONS

Both radiotracers demonstrate appropriate imaging properties for the MCHR1. However, the pronounced metabolic stability as well as superior selectivity and affinity of [11C]SNAP-7941 underlines the decisive superiority over [18F]FE@SNAP.

A two-state model for the kinetics of competitive radioligand binding

Background and Purpose

Ligand–receptor binding kinetics is receiving increasing attention in the drug research community. The Motulsky and Mahan model, a one‐state model, offers a method for measuring the binding kinetics of an unlabelled ligand, with the assumption that the labelled ligand has no preference while binding to distinct states or conformations of a drug target. As such, the one‐state model is not applicable if the radioligand displays biphasic binding kinetics to the receptor.

Experimental Approach

We extended the Motulsky and Mahan model to a two‐state model, in which the kinetics of the unlabelled competitor binding to different receptor states (R₁ and R₂) can be measured. With this extended model, we determined the binding kinetics of unlabelled N‐5′‐ethylcarboxamidoadenosine (NECA), a representative agonist for the adenosine A₁ receptor. Subsequently, an application of the model was exemplified by measuring the binding kinetics of other A₁ receptor ligands. In addition, limitations of the model were investigated as well.

Key Results

The kinetic rate constants of unlabelled NECA were comparable with the results of kinetic radioligand binding assays in which [³H]‐NECA was used. The model was further validated by good correlation between simulated results and the experimental data.

Conclusion

The two‐state model is sufficient to analyse the binding kinetics of an unlabelled ligand, when a radioligand shows biphasic association characteristics. We expect this two‐state model to have general applicability for other targets as well.

Covalent Allosteric Probe for the Metabotropic Glutamate Receptor 2: Design, Synthesis, and Pharmacological Characterization

Covalent labeling of G protein-coupled receptors (GPCRs) by small molecules is a powerful approach to understand binding modes, mechanism of action, pharmacology, and even facilitate structure elucidation. We report the first covalent positive allosteric modulator (PAM) for a class C GPCR, the mGlu₂ receptor. Three putatively covalent mGlu₂ PAMs were designed and synthesized. Pharmacological characterization identified 2 to bind the receptor covalently. Computational modeling combined with receptor mutagenesis revealed T791^7.29×30 as the likely position of covalent interaction. We show how this covalent ligand can be used to characterize the PAM binding mode and that it is a valuable tool compound in studying receptor function and binding kinetics. Our findings advance the understanding of the mGlu₂ PAM interaction and suggest that 2 is a valuable probe for further structural and chemical biology approaches.

A binding kinetics study of human adenosine A3 receptor agonists

The human adenosine A₃ (hA₃) receptor has been suggested as a viable drug target in inflammatory diseases and in cancer. So far, a number of selective hA₃ receptor agonists (e.g. IB-MECA and 2-Cl-IB-MECA) inducing anti-inflammatory or anticancer effects are under clinical investigation. Drug-target binding kinetics is increasingly recognized as another pharmacological parameter, next to affinity, for compound triage in the early phases of drug discovery. However, such a kinetics-driven analysis has not yet been performed for the hA₃ receptor. In this study, we first validated a competition association assay for adenosine A₃ receptor agonists to determine the target interaction kinetics. Affinities and Kinetic Rate Index (KRI) values of 11 ribofurano and 10 methanocarba nucleosides were determined in radioligand binding assays. Afterwards, 15 analogues were further selected (KRI <0.70 or KRI >1.35) for full kinetics characterization. The structure-kinetics relationships (SKR) were derived and longer residence times were associated with methanocarba and enlarged adenine N⁶ and C2 substitutions. In addition, from a k_on-k_off-K_D kinetic map we divided the agonists into three subgroups. A residence time "cliff" was observed, which might be relevant to (N)-methanocarba derivatives' rigid C2-arylalkynyl substitutions. Our findings provide substantial evidence that, next to affinity, additional knowledge of binding kinetics is useful for developing and selecting new hA₃R agonists in the early phase of the drug discovery process.

Structure-Affinity Relationships and Structure-Kinetic Relationships of 1,2-Diarylimidazol-4-carboxamide Derivatives as Human Cannabinoid 1 Receptor Antagonists

We report on the synthesis and biological evaluation of a series of 1,2-diarylimidazol-4-carboxamide derivatives developed as CB₁ receptor antagonists. These were evaluated in a radioligand displacement binding assay, a [³⁵S]GTPγS binding assay, and in a competition association assay that enables the relatively fast kinetic screening of multiple compounds. The compounds show high affinities and a diverse range of kinetic profiles at the CB₁ receptor and their structure–kinetic relationships (SKRs) were established. Using the recently resolved hCB₁ receptor crystal structures, we also performed a modeling study that sheds light on the crucial interactions for both the affinity and dissociation kinetics of this family of ligands. We provide evidence that, next to affinity, additional knowledge of binding kinetics is useful for selecting new hCB₁ receptor antagonists in the early phases of drug discovery.

Kinetics of human cannabinoid 1 (CB1) receptor antagonists: Structure-kinetics relationships (SKR) and implications for insurmountable antagonism

While equilibrium binding affinities and in vitro functional antagonism of CB1 receptor antagonists have been studied in detail, little is known on the kinetics of their receptor interaction. In this study, we therefore conducted kinetic assays for nine 1-(4,5-diarylthiophene-2-carbonyl)-4-phenylpiperidine-4-carboxamide derivatives and included the CB1 antagonist rimonabant as a comparison. For this we newly developed a dual-point competition association assay with [³H]CP55940 as the radioligand. This assay yielded Kinetic Rate Index (KRI) values from which structure-kinetics relationships (SKR) of hCB1 receptor antagonists could be established. The fast dissociating antagonist 6 had a similar receptor residence time (RT) as rimonabant, i.e. 19 and 14 min, respectively, while the slowest dissociating antagonist (9) had a very long RT of 2222 min, i.e. pseudo-irreversible dissociation kinetics. In functional assays, 9 displayed insurmountable antagonism, while the effects of the shortest RT antagonist 6 and rimonabant were surmountable. Taken together, this study shows that hCB1 receptor antagonists can have very divergent RTs, which are not correlated to their equilibrium affinities. Furthermore, their RTs appear to define their mode of functional antagonism, i.e. surmountable vs. insurmountable. Finally, based on the recently resolved hCB1 receptor crystal structure, we propose that the differences in RT can be explained by a different binding mode of antagonist 9 from short RT antagonists that is able to displace unfavorable water molecules. Taken together, these findings are of importance for future design and evaluation of potent and safe hCB1 receptor antagonists.

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.

Structure-Affinity Relationships and Structure-Kinetics Relationships of Pyrido[2,1-f]purine-2,4-dione Derivatives as Human Adenosine A3 Receptor Antagonists

We expanded on a series of pyrido[2,1-f]purine-2,4-dione derivatives as human adenosine A₃ receptor (hA₃R) antagonists to determine their kinetic profiles and affinities. Many compounds showed high affinities and a diverse range of kinetic profiles. We found hA₃R antagonists with very short residence time (RT) at the receptor (2.2 min for 5) and much longer RTs (e.g., 376 min for 27 or 391 min for 31). Two representative antagonists (5 and 27) were tested in [³⁵S]GTPγS binding assays, and their RTs appeared correlated to their (in)surmountable antagonism. From a k_on–k_off–K_D kinetic map, we divided the antagonists into three subgroups, providing a possible direction for the further development of hA₃R antagonists. Additionally, we performed a computational modeling study that sheds light on the crucial receptor interactions, dictating the compounds’ binding kinetics. Knowledge of target binding kinetics appears useful for developing and triaging new hA₃R antagonists in the early phase of drug discovery.

Binding Kinetics and Pathways of Ligands to GPCRs

Previously, drugs were developed focusing on target affinity and selectivity. However, it is becoming evident that the drug-target residence time, related to the off-rate, is an important parameter for successful drug development. The residence time influences both the on-rate and overall effectiveness of drugs. Furthermore, ligand binding is now appreciated to be a multistep process because metastable and/or intermediate binding sites in the extracellular region have been identified. In this review, we summarize experimental ligand-binding data for G-protein-coupled receptors (GPCRs), and their binding pathways, analyzed by molecular dynamics (MD). The kinetics of drug binding to GPCRs are complex and depend on several factors, including charge distribution on the receptor surface, ligand-receptor interactions in the binding channel and the binding site, or solvation.

Kinetics for Drug Discovery: an industry-driven effort to target drug residence time

A considerable number of approved drugs show non-equilibrium binding characteristics, emphasizing the potential role of drug residence times for in vivo efficacy. Therefore, a detailed understanding of the kinetics of association and dissociation of a target-ligand complex might provide crucial insight into the molecular mechanism-of-action of a compound. This deeper understanding will help to improve decision making in drug discovery, thus leading to a better selection of interesting compounds to be profiled further. In this review, we highlight the contributions of the Kinetics for Drug Discovery (K4DD) Consortium, which targets major open questions related to binding kinetics in an industry-driven public-private partnership.

New approaches for the reliable in vitro assessment of binding affinity based on high-resolution real-time data acquisition of radioligand-receptor binding kinetics

Background

Resolving the kinetic mechanisms of biomolecular interactions have become increasingly important in early-phase drug development. Since traditional in vitro methods belong to dose-dependent assessments, binding kinetics is usually overlooked. The present study aimed at the establishment of two novel experimental approaches for the assessment of binding affinity of both, radiolabelled and non-labelled compounds targeting the A₃R, based on high-resolution real-time data acquisition of radioligand-receptor binding kinetics. A novel time-resolved competition assay was developed and applied to determine the K_i of eight different A₃R antagonists, using CHO-K1 cells stably expressing the hA₃R. In addition, a new kinetic real-time cell-binding approach was established to quantify the rate constants k_on and k_off, as well as the dedicated K_d of the A₃R agonist [¹²⁵I]-AB-MECA. Furthermore, lipophilicity measurements were conducted to control influences due to physicochemical properties of the used compounds.

Results

Two novel real-time cell-binding approaches were successfully developed and established. Both experimental procedures were found to visualize the kinetic binding characteristics with high spatial and temporal resolution, resulting in reliable affinity values, which are in good agreement with values previously reported with traditional methods. Taking into account the lipophilicity of the A₃R antagonists, no influences on the experimental performance and the resulting affinity were investigated.

Conclusions

Both kinetic binding approaches comprise tracer administration and subsequent binding to living cells, expressing the dedicated target protein. Therefore, the experiments resemble better the true in vivo physiological conditions and provide important markers of cellular feedback and biological response.

A novel CCR2 antagonist inhibits atherogenesis in apoE deficient mice by achieving high receptor occupancy

CC Chemokine Receptor 2 (CCR2) and its endogenous ligand CCL2 are involved in a number of diseases, including atherosclerosis. Several CCR2 antagonists have been developed as potential therapeutic agents, however their in vivo clinical efficacy was limited. In this report, we aimed to determine whether 15a, an antagonist with a long residence time on the human CCR2, is effective in inhibiting the development of atherosclerosis in a mouse disease model. First, radioligand binding assays were performed to determine affinity and binding kinetics of 15a on murine CCR2. To assess the in vivo efficacy, western-type diet fed apoE-/- mice were treated daily with 15a or vehicle as control. Treatment with 15a reduced the amount of circulating CCR2+ monocytes and the size of the atherosclerotic plaques in both the carotid artery and the aortic root. We then showed that the long pharmacokinetic half-life of 15a combined with the high drug concentrations ensured prolonged CCR2 occupancy. These data render 15a a promising compound for drug development and confirms high receptor occupancy as a key parameter when targeting chemokine receptors.

Structure-Activity Relationships of the Sustained Effects of Adenosine A2A Receptor Agonists Driven by Slow Dissociation Kinetics

The duration of action of adenosine A_2A receptor (A2A) agonists is critical for their clinical efficacy, and we sought to better understand how this can be optimized. The in vitro temporal response profiles of a panel of A2A agonists were studied using cAMP assays in recombinantly (CHO) and endogenously (SH-SY5Y) expressing cells. Some agonists (e.g., 3cd; UK-432,097) but not others (e.g., 3ac; CGS-21680) demonstrated sustained wash-resistant agonism, where residual receptor activation continued after washout. The ability of an antagonist to reverse pre-established agonist responses was used as a surrogate read-out for agonist dissociation kinetics, and together with radioligand binding studies suggested a role for slow off-rate in driving sustained effects. One compound, 3ch, showed particularly marked sustained effects, with a reversal t_1/2 > 6 hours and close to maximal effects that remained for at least 5 hours after washing. Based on the structure-activity relationship of these compounds, we suggest that lipophilic N6 and bulky C2 substituents can promote stable and long-lived binding events leading to sustained agonist responses, although a high compound logD is not necessary. This provides new insight into the binding interactions of these ligands and we anticipate that this information could facilitate the rational design of novel long-acting A2A agonists with improved clinical efficacy.

Long Residence Time at the Neurosteroidogenic 18 kDa Translocator Protein Characterizes the Anxiolytic Ligand XBD173

Recent data have demonstrated a positive correlation between the residence time (RT) and neurosteroidogenic efficacy of a ligand at the translocator protein (TSPO), an attractive anxyolitic target. To explore the potential impact of RT on TSPO ligand anxiolytic activity, the RT and the steroidogenic activity of XBD173, a ligand exerting anxiolytic activity in humans, were retrospectively evaluated. To this aim, XBD173 association and dissociation rate constants were measured (1.23 × 10(7) M(-1) min(-1) and 0.0079 min(-1), respectively). XBD173 resulted to have a long RT (127 min) and to stimulate efficaciously neurosteroidogenesis, in terms of pregnenolone production. The present findings corroborate the importance of TSPO ligand RT to predict their effective neurosteroidogenic activity and promising anxiolytic action. These positive results prompted us to set up a fast and high-throughput kinetic method to improve the efficiency of RT-based TSPO drug-discovery process.

Characterization of Kinetic Binding Properties of Unlabeled Ligands via a Preincubation Endpoint Binding Approach

The dissociation rates of unlabeled drugs have been well studied by kinetic binding analyses. Since kinetic assays are laborious, we developed a simple method to determine the kinetic binding parameters of unlabeled competitors by a preincubation endpoint assay. The probe binding after preincubation of a competitor can be described by a single equation as a function of time. Simulations using the equation revealed the degree of IC50 change induced by preincubation of a competitor depended on the dissociation rate koff of the competitor but not on the association rate kon To validate the model, an in vitro binding assay was performed using a smoothened receptor (SMO) and [(3)H]TAK-441, a SMO antagonist. The equilibrium dissociation constants (KI) and koff of SMO antagonists determined by globally fitting the model to the concentration-response curves obtained with and without 24 h preincubation correlated well with those determined by other methods. This approach could be useful for early-stage optimization of drug candidates by enabling determination of binding kinetics in a high-throughput manner because it does not require kinetic measurements, an intermediate washout step during the reaction, or prior determination of competitors' KI values.

Binding Kinetics in Drug Discovery

Over the last years, researchers have increasingly become interested in measuring and understanding drugs' binding kinetics, namely the time in which drug and its target associate and dissociate. Historically, drug discovery programs focused on the optimization of target affinity as a proxy of in-vivo efficacy. However, often the efficacy of a ligand is not appropriately described by the in-vitro measured drug-receptor affinity, but rather depends on the lifetime of the in-vivo drug-receptor interaction. In this review we review recent works that highlight the importance of binding kinetics, molecular determinants for rational optimization and the recent emergence of computational methods as powerful tools in measuring and understanding binding kinetics.

A Four-Point Screening Method for Assessing Molecular Mechanism of Action (MMOA) Identifies Tideglusib as a Time-Dependent Inhibitor of Trypanosoma brucei GSK3β

Background

New therapeutics are needed for neglected tropical diseases including Human African trypanosomiasis (HAT), a progressive and fatal disease caused by the protozoan parasites Trypanosoma brucei gambiense and T. b. rhodesiense. There is a need for simple, efficient, cost effective methods to identify new molecules with unique molecular mechanisms of action (MMOAs). The mechanistic features of a binding mode, such as competition with endogenous substrates and time-dependence can affect the observed inhibitory IC₅₀, and differentiate molecules and their therapeutic usefulness. Simple screening methods to determine time-dependence and competition can be used to differentiate compounds with different MMOAs in order to identify new therapeutic opportunities.

Methodology/Principal Findings

In this work we report a four point screening methodology to evaluate the time-dependence and competition for inhibition of GSK3β protein kinase isolated from T. brucei. Using this method, we identified tideglusib as a time-dependent inhibitor whose mechanism of action is time-dependent, ATP competitive upon initial binding, which transitions to ATP non-competitive with time. The enzyme activity was not recovered following 100-fold dilution of the buffer consistent with an irreversible mechanism of action. This is in contrast to the T. brucei GSK3β inhibitor GW8510, whose inhibition was competitive with ATP, not time-dependent at all measured time points and reversible in dilution experiments. The activity of tideglusib against T. brucei parasites was confirmed by inhibition of parasite proliferation (GI₅₀ of 2.3 μM).

Conclusions/Significance

Altogether this work demonstrates a straightforward method for determining molecular mechanisms of action and its application for mechanistic differentiation of two potent TbGSK3β inhibitors. The four point MMOA method identified tideglusib as a mechanistically differentiated TbGSK3β inhibitor. Tideglusib was shown to inhibit parasite growth in this work, and has been reported to be well tolerated in one year of dosing in human clinical studies. Consequently, further supportive studies on the potential therapeutic usefulness of tideglusib for HAT are justified.

Drug discovery for neglected tropical diseases must use efficient methods due to limited resources. One preferred drug discovery strategy is target-based drug discovery. In this strategy it is assumed that drug action begins with binding of a drug to its target. However, while binding is required, it is not sufficient to describe all the molecular interactions that translate binding to a therapeutically useful response. The contribution of aspects of the molecular mechanism of action (MMOA) such as time-dependence and substrate competition can influence concentration response relationships. To address this, a four point MMOA methodology was developed to evaluate time-dependence and substrate competition. We used this method to evaluate the MMOA for T.brucei GSK3β inhibitors, and observed tideglusib to have a time-dependent, ATP-competitive mechanism that differentiated it from rapidly reversible inhibitors, such as GW8510. Adjusting the enzyme assays to account for these mechanisms showed that GW8510 and tideglusib had similar activities for TbGSK3β. However, this similarity did not translate to cellular activity, where GW-8510 was more active than tideglusib (0.12 μM to 2.3 μM, respectively). These data suggest that factors other than TbGSK3β MMOA differentiate the effect of these molecules against T. brucei.

Scintillation proximity assay (SPA) as a new approach to determine a ligand's kinetic profile. A case in point for the adenosine A1 receptor

Scintillation proximity assay (SPA) is a radio-isotopic technology format used to measure a wide range of biological interactions, including drug-target binding affinity studies. The assay is homogeneous in nature, as it relies on a “mix and measure” format. It does not involve a filtration step to separate bound from free ligand as is the case in a traditional receptor-binding assay. For G protein-coupled receptors (GPCRs), it has been shown that optimal binding kinetics, next to a high affinity of a ligand, can result in more desirable pharmacological profiles. However, traditional techniques to assess kinetic parameters tend to be cumbersome and laborious. We thus aimed to evaluate whether SPA can be an alternative platform for real-time receptor-binding kinetic measurements on GPCRs. To do so, we first validated the SPA technology for equilibrium binding studies on a prototypic class A GPCR, the human adenosine A₁ receptor (hA₁R). Differently to classic kinetic studies, the SPA technology allowed us to study binding kinetic processes almost real time, which is impossible in the filtration assay. To demonstrate the reliability of this technology for kinetic purposes, we performed the so-called competition association experiments. The association and dissociation rate constants (k_on and k_off) of unlabeled hA₁R ligands were reliably and quickly determined and agreed very well with the same parameters from a traditional filtration assay performed simultaneously. In conclusion, SPA is a very promising technique to determine the kinetic profile of the drug-target interaction. Its robustness and potential for high-throughput may render this technology a preferred choice for further kinetic studies.

Mechanistic models enable the rational use of in vitro drug-target binding kinetics for better drug effects in patients

INTRODUCTION

Drug-target binding kinetics are major determinants of the time course of drug action for several drugs, as clearly described for the irreversible binders omeprazole and aspirin. This supports the increasing interest to incorporate newly developed high-throughput assays for drug-target binding kinetics in drug discovery. A meaningful application of in vitro drug-target binding kinetics in drug discovery requires insight into the relation between in vivo drug effect and in vitro measured drug-target binding kinetics.

AREAS COVERED

In this review, the authors discuss both the relation between in vitro and in vivo measured binding kinetics and the relation between in vivo binding kinetics, target occupancy and effect profiles.

EXPERT OPINION

More scientific evidence is required for the rational selection and development of drug-candidates on the basis of in vitro estimates of drug-target binding kinetics. To elucidate the value of in vitro binding kinetics measurements, it is necessary to obtain information on system-specific properties which influence the kinetics of target occupancy and drug effect. Mathematical integration of this information enables the identification of drug-specific properties which lead to optimal target occupancy and drug effect in patients.

Mass spectrometry-based ligand binding assays on adenosine A1 and A2A receptors

Conventional methods to measure ligand-receptor binding parameters typically require radiolabeled ligands as probes. Despite the robustness of radioligand binding assays, they carry inherent disadvantages in terms of safety precautions, expensive synthesis, special lab requirements, and waste disposal. Mass spectrometry (MS) is a method that can selectively detect ligands without the need of a label. The sensitivity of MS equipment increases progressively, and currently, it is possible to detect low ligand quantities that are usually found in ligand binding assays. We developed a label-free MS ligand binding (MS binding) assay on the adenosine A₁ and A_2A receptors (A₁AR and A_2AAR), which are well-characterized members of the class A G protein-coupled receptor (GPCR) family. Radioligand binding assays for both receptors are well established, and ample data is available to compare and evaluate the performance of an MS binding assay. 1,3-Dipropyl-8-cyclopentyl-xanthine (DPCPX) and 4-(2-((7-amino-2-(furan-2-yl)-[1,2,4]triazolo[1,5-a]-[1,3,5]triazin-5-yl)amino)ethyl)phenol (ZM-241,385) are high-affinity ligands selective for the A₁AR and A_2AAR, respectively. To proof the feasibility of MS binding on the A₁AR and A_2AAR, we first developed an MS detection method for unlabeled DPCPX and ZM-241,385. To serve as internal standards, both compounds were also deuterium-labeled. Subsequently, we investigated whether the two unlabeled compounds could substitute for their radiolabeled counterparts as marker ligands in binding experiments, including saturation, displacement, dissociation, and competition association assays. Furthermore, we investigated the accuracy of these assays if the use of internal standards was excluded. The results demonstrate the feasibility of the MS binding assay, even in the absence of a deuterium-labeled internal standard, and provide great promise for the further development of label-free assays based on MS for other GPCRs.

Radioligand binding to intact cells as a tool for extended drug screening in a representative physiological context

Radioligand binding assays on intact cells offer distinct advantages to those on membrane suspensions. Major pharmacological properties like drug affinity and binding kinetics are more physiologically relevant. Complex mechanisms can be studied with a wider choice of experimental approaches and so provide insights into induced-fit type binding, receptor internalisation and even into pharmacomicrokinetic phenomena like drug rebinding and partitioning into the membrane. Hence, intact cell binding constitutes a valuable addition to the pharmacologist's toolbox.

Partial agonism and fast dissociation of LASSBio-579 at dopamine D2 receptor

In an attempt to better understand the molecular mechanism of action of the antipsychotic lead LASSBio-579 and of its main metabolite LQFM 037, the aim of this work was to evaluate their intrinsic activity and binding kinetics at the dopamine D2 receptor. In transfected HEK cells expressing the D2L receptor under an inducible promoter, LASSBio-579 and LQFM 037, but not clozapine, behaved as weak partial agonists in [(35)S]-GTPγS binding assays performed in optimized conditions previously shown to evidence the partial agonist profile of aripiprazole. Besides, data obtained in radioligand competition assays on rat striatal membranes suggested a rapid association to and dissociation from the D2-like receptors. Using the kinetic rate index based on the strategy of the dual-point competition association assay, we showed that our compounds share a similar kinetic profile with clozapine, distinct from the typical antipsychotic haloperidol. These two characteristics could contribute to the atypical-like profile observed after administration of LASSBio-579 to rodents, in models of positive and negative symptoms of schizophrenia.

Affinity and kinetics study of anthranilic acids as HCA2 receptor agonists

Structure-affinity relationship (SAR) and structure-kinetics relationship (SKR) studies were combined to investigate a series of biphenyl anthranilic acid agonists for the HCA2 receptor. In total, 27 compounds were synthesized and twelve of them showed higher affinity than nicotinic acid. Two compounds, 6g (IC50=75nM) and 6z (IC50=108nM) showed a longer residence time profile compared to nicotinic acid, exemplified by their kinetic rate index (KRI) values of 1.31 and 1.23, respectively. The SAR study resulted in the novel 2-F, 4-OH derivative (6x) with an IC50 value of 23nM as the highest affinity HCA2 agonist of the biphenyl series, although it showed a similar residence time as nicotinic acid. The SAR and SKR data suggest that an early compound selection based on binding kinetics is a promising addition to the lead optimization process.