Estimation of the dissociation rate of unlabelled ligand-receptor complexes by a 'two-step' competition binding approach

Ligand with Two Modes of Interaction with the Dopamine D2 Receptor-An Induced-Fit Mechanism of Insurmountable Antagonism

A solid understanding of the mechanisms governing ligand binding is crucial for rational design of therapeutics targeting the dopamine D₂ receptor (D₂R). Here, we use G protein-coupled inward rectifier potassium (GIRK) channel activation in Xenopus oocytes to measure the kinetics of D₂R antagonism by a series of aripiprazole analogues, as well as the recovery of dopamine (DA) responsivity upon washout. The aripiprazole analogues comprise an orthosteric and a secondary pharmacophore and differ by the length of the saturated carbon linker joining these two pharmacophores. Two compounds containing 3- and 5-carbon linkers allowed for a similar extent of recovery from antagonism in the presence of 1 or 100 μM DA (>25 and >90% of control, respectively), whereas recovery was less prominent (∼20%) upon washout of the 4-carbon linker compound, SV-III-130, both with 1 and 100 μM DA. Prolonging the coincubation time with SV-III-130 further diminished recovery. Curve-shift experiments were consistent with competition between SV-III-130 and DA. Two mutations in the secondary binding pocket (V91A and E95A) of D₂R decreased antagonistic potency and increased recovery from SV-III-130 antagonism, whereas a third mutation (L94A) only increased recovery. Our results suggest that the secondary binding pocket influences recovery from inhibition by the studied aripiprazole analogues. We propose a mechanism, supported by in silico modeling, whereby SV-III-130 initially binds reversibly to the D₂R, after which the drug-receptor complex undergoes a slow transition to a second ligand-bound state, which is dependent on secondary binding pocket integrity and irreversible during the time frame of our experiments.

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.

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.

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.

Development of a population pharmacokinetic model to predict brain distribution and dopamine D2 receptor occupancy of raclopride in non-anesthetized rat

BACKGROUND

Raclopride is a selective antagonist of the dopamine D2 receptor. It is one of the most frequently used in vivo D2 tracers (at low doses) for assessing drug-induced receptor occupancy (RO) in animals and humans. It is also commonly used as a pharmacological blocker (at high doses) to occupy the available D2 receptors and antagonize the action of dopamine or drugs on D2 in preclinical studies. The aims of this study were to comprehensively evaluate its pharmacokinetic (PK) profiles in different brain compartments and to establish a PK-RO model that could predict the brain distribution and RO of raclopride in the freely moving rat using a LC-MS based approach.

METHODS

Rats (n=24) received a 10-min IV infusion of non-radiolabeled raclopride (1.61μmol/kg, i.e. 0.56mg/kg). Plasma and the brain tissues of striatum (with high density of D2 receptors) and cerebellum (with negligible amount of D2 receptors) were collected. Additional microdialysis experiments were performed in some rats (n=7) to measure the free drug concentration in the extracellular fluid of the striatum and cerebellum. Raclopride concentrations in all samples were analyzed by LC-MS. A population PK-RO model was constructed in NONMEM to describe the concentration-time profiles in the unbound plasma, brain extracellular fluid and brain tissue compartments and to estimate the RO based on raclopride-D2 receptor binding kinetics.

RESULTS

In plasma raclopride showed a rapid distribution phase followed by a slower elimination phase. The striatum tissue concentrations were consistently higher than that of cerebellum tissue throughout the whole experimental period (10-h) due to higher non-specific tissue binding and D2 receptor binding in the striatum. Model-based simulations accurately predicted the literature data on rat plasma PK, brain tissue PK and D2 RO at different time points after intravenous or subcutaneous administration of raclopride at tracer dose (RO <10%), sub-pharmacological dose (RO 10%-30%) and pharmacological dose (RO >30%).

CONCLUSION

For the first time a predictive model that could describe the quantitative in vivo relationship between dose, PK and D2 RO of raclopride in non-anesthetized rat was established. The PK-RO model could facilitate the selection of optimal dose and dosing time when raclopride is used as tracer or as pharmacological blocker in various rat studies. The LC-MS based approach, which doses and quantifies a non-radiolabeled tracer, could be useful in evaluating the systemic disposition and brain kinetics of tracers.

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.

Cell membranes… and how long drugs may exert beneficial pharmacological activity in vivo

The time course of the beneficial pharmacological effect of a drug has long been considered to depend merely on the temporal fluctuation of its free concentration. Only in the last decade has it become widely accepted that target-binding kinetics can also affect in vivo pharmacological activity. Although current reviews still essentially focus on genuine dissociation rates, evidence is accumulating that additional micro-pharmacokinetic (PK) and -pharmacodynamic (PD) mechanisms, in which the cell membrane plays a central role, may also increase the residence time of a drug on its target. The present review provides a compilation of otherwise widely dispersed information on this topic. The cell membrane can intervene in drug binding via the following three major mechanisms: (i) by acting as a sink/repository for the drug; (ii) by modulating the conformation of the drug and even by participating in the binding process; and (iii) by facilitating the approach (and rebinding) of the drug to the target. To highlight these mechanisms, we focus on drugs that are currently used in clinical therapy, such as the antihypertensive angiotensin II type 1 receptor antagonist candesartan, the atypical antipsychotic agent clozapine and the bronchodilator salmeterol. Although the role of cell membranes in PK-PD modelling is gaining increasing interest, many issues remain unresolved. It is likely that novel biophysical and computational approaches will provide improved insights in the near future.

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.

The fast-off hypothesis revisited: A functional kinetic study of antipsychotic antagonism of the dopamine D2 receptor

Newer, "atypical" antipsychotics carry a lower risk of motor side-effects than older, "typical" compounds. It has been proposed that a ~100-fold faster dissociation from the dopamine D2 receptor (D2R) distinguishes atypical from typical antipsychotics. Furthermore, differing antipsychotic D2R affinities have been suggested to reflect differences in dissociation rate constants (koff), while association rate constants (kon) were assumed to be similar. However, it was recently demonstrated that lipophilic accumulation of ligand in the cell interior and/or membrane can cause underestimation of koff, and as high-affinity D2R antagonists are frequently lipophilic, this may have been a confounding factor in previous studies. In the present work, a functional electrophysiology assay was used to measure the recovery of dopamine-mediated D2R responsivity from antipsychotic antagonism, using elevated concentrations of dopamine to prevent the potential bias of re-binding of lipophilic ligands. The variability of antipsychotic kon was also reexamined, capitalizing on the temporal resolution of the assay. kon was estimated from the experimental recordings using a simple mathematical model assumed to describe the binding process. The time course of recovery from haloperidol (typical antipsychotic) was only 6.4- to 2.5-fold slower than that of the atypical antipsychotics, amisulpride, clozapine, and quetiapine, while antipsychotic kons were found to vary more widely than previously suggested. Finally, affinities calculated using our kon and koff estimates correlated well with functional potency and with affinities reported from radioligand binding studies. In light of these findings, it appears unlikely that typical and atypical antipsychotics are primarily distinguished by their D2R binding kinetics.

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.

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.

On the 'micro'-pharmacodynamic and pharmacokinetic mechanisms that contribute to long-lasting drug action

INTRODUCTION

Optimal drug therapy often requires continuing high levels of target occupancy. Besides the traditional pharmacokinetic (PK) contribution thereto, drug-target interactions that comprise successive 'microscopic' steps as well as the intervention of the cell membrane and other 'micro'-anatomical structures nearby may help attaining this objective.

AREAS COVERED

This article reviews the 'micro'-pharmacodynamic (PD) and PK mechanisms that may increase a drug's residence time. Special focus is on induced-fit- and bivalent ligand binding models as well as on the ability of the plasma membrane surrounding the target to act as a repository for the drug (e.g., microkinetic model), to actively participate in the binding process (e.g., exosite model) and, along with microanatomical elements like synapses and interstitial spaces, to act on the drug's diffusion properties (reduction in dimensionality and drug-rebinding models).

EXPERT OPINION

The PK profile, as well as the target dissociation kinetics of a drug, may fail to account for its long-lasting efficiency in intact tissues and in vivo. This lacuna could potentially be alleviated by incorporating some of the enumerated 'microscopic' mechanisms and, to unveil them, dedicated experiments on sufficiently physiologically relevant biological material like cell monolayers can already be implemented early on in the lead optimization process.

Typical and atypical antipsychotics do not differ markedly in their reversibility of antagonism of the dopamine D2 receptor

It has been suggested that the favorable side-effect profiles of atypical antipsychotics (e.g. clozapine and amisulpride) are related to their ∼100-fold faster dissociation from dopamine D2 receptors (D2R) compared with typical antipsychotics (e.g. haloperidol and chlorpromazine). Fast dissociation would entail rapidly reversible antagonism; however, this has not been thoroughly studied using functional assays. We compared the reversibilities of D2R antagonism by 17 compounds using an electrophysiological method to measure dopamine-evoked potassium channel activation via D2R. Varying rates and amplitudes of D2R response recovery were observed following antagonist removal. Whereas recovery rates differed 15-fold between atypical drugs, recovery from clozapine and amisulpride antagonism was, unexpectedly, less than twofold faster than from chlorpromazine. The recovery amplitude correlated with calculated water solubility and lipid/water distribution coefficients, suggesting variable drug partitioning into cell membranes. Our data do not support the notion that the rate of reversibility of D2R antagonism is what distinguishes atypical from typical antipsychotics.

Structure-kinetic relationships--an overlooked parameter in hit-to-lead optimization: a case of cyclopentylamines as chemokine receptor 2 antagonists

Preclinical models of inflammatory diseases (e.g., neuropathic pain, rheumatoid arthritis, and multiple sclerosis) have pointed to a critical role of the chemokine receptor 2 (CCR2) and chemokine ligand 2 (CCL2). However, one of the biggest problems of high-affinity inhibitors of CCR2 is their lack of efficacy in clinical trials. We report a new approach for the design of high-affinity and long-residence-time CCR2 antagonists. We developed a new competition association assay for CCR2, which allows us to investigate the relation of the structure of the ligand and its receptor residence time [i.e., structure-kinetic relationship (SKR)] next to a traditional structure-affinity relationship (SAR). By applying combined knowledge of SAR and SKR, we were able to re-evaluate the hit-to-lead process of cyclopentylamines as CCR2 antagonists. Affinity-based optimization yielded compound 1 with good binding (Ki = 6.8 nM) but very short residence time (2.4 min). However, when the optimization was also based on residence time, the hit-to-lead process yielded compound 22a, a new high-affinity CCR2 antagonist (3.6 nM), with a residence time of 135 min.

Dual-Point Competition Association Assay

The concept of ligand-receptor binding kinetics is emerging as an important parameter in the early phase of drug discovery. Since the currently used kinetic assays are laborious and low throughput, we developed a method that enables fast and large format screening. It is a so-called dual-point competition association assay, which measures radioligand binding at two different time points in the absence or presence of unlabeled competitors. Specifically, this assay yields the kinetic rate index (KRI), which is a measure for the binding kinetics of the unlabeled ligands screened. As a prototypical drug target, the adenosine A1 receptor (A1R) was chosen for assay validation and optimization. A screen with 35 high-affinity A1R antagonists yielded seven compounds with a KRI value above 1.0, which indicated a relatively slow dissociation from the target. All other compounds had a KRI value below or equal to 1.0, predicting a relatively fast dissociation rate. Several compounds were selected for follow-up kinetic quantifications in classical kinetic assays and were shown to have kinetic rates that corresponded to their KRI values. The dual-point assay and KRI value may have general applicability at other G-protein-coupled receptors, as well as at drug targets from other protein families.

Clozapine, atypical antipsychotics, and the benefits of fast-off D2 dopamine receptor antagonism

Drug-receptor interactions are traditionally quantified in terms of affinity and efficacy, but there is increasing awareness that the drug-on-receptor residence time also affects clinical performance. While most interest has hitherto been focused on slow-dissociating drugs, D(2) dopamine receptor antagonists show less extrapyramidal side effects but still have excellent antipsychotic activity when they dissociate swiftly. Fast dissociation of clozapine, the prototype of the "atypical antipsychotics", has been evidenced by distinct radioligand binding approaches both on cell membranes and intact cells. The surmountable nature of clozapine in functional assays with fast-emerging responses like calcium transients is confirmatory. Potential advantages and pitfalls of the hitherto used techniques are discussed, and recommendations are given to obtain more precise dissociation rates for such drugs. Surmountable antagonism is necessary to allow sufficient D(2) receptor stimulation by endogenous dopamine in the striatum. Simulations are presented to find out whether this can be achieved during sub-second bursts in dopamine concentration or rather during much slower, activity-related increases thereof. While the antagonist's dissociation rate is important to distinguish between both mechanisms, this becomes much less so when contemplating time intervals between successive drug intakes, i.e., when pharmacokinetic considerations prevail. Attention is also drawn to the divergent residence times of hydrophobic antagonists like haloperidol when comparing radioligand binding data on cell membranes with those on intact cells and clinical data.

Binding of "AT4 receptor" ligands to insulin regulated aminopeptidase (IRAP) in intact Chinese hamster ovary cells

Insulin regulated aminopeptidase (IRAP) recognises "AT(4)-receptor" ligands like angiotensin IV (Ang IV) and peptidomimetics like AL-11. The metabolic stability and high affinity of [(3)H]AL-11 for catalytically active IRAP allowed its detection in Chinese hamster ovary (CHO-K1) cell membranes in the absence of chelators (Demaegdt et al., 2009). Here, we show that, contrary to [(3)H]Ang IV, [(3)H]AL-11 displays high affinity and specificity for IRAP in intact CHO-K1 cells as well. After binding to IRAP at the surface, [(3)H]AL-11 is effectively internalized by an endocytotic process. Unexpectedly, surface binding and internalization of [(3)H]AL-11 was not affected by pretreating the cells with Ang IV but declined with AL-11. In the latter case surface expression of IRAP even increased. After elimination of simpler explanations, it is proposed that metabolically stable "AT(4)-receptor" ligands undergo semi-continuous cycling between the cell surface and endosomal compartments. The in vivo efficacy of stable and unstable "AT(4)-receptor" ligands could therefore differ.