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

Kinetic profiling of novel spirobenzo-oxazinepiperidinone derivatives as equilibrative nucleoside transporter 1 inhibitors

Evaluation of kinetic parameters of drug-target binding, kon, koff, and residence time (RT), in addition to the traditional in vitro parameter of affinity is receiving increasing attention in the early stages of drug discovery. Target binding kinetics emerges as a meaningful concept for the evaluation of a ligand's duration of action and more generally drug efficacy and safety. We report the biological evaluation of a novel series of spirobenzo-oxazinepiperidinone derivatives as inhibitors of the human equilibrative nucleoside transporter 1 (hENT1, SLC29A1). The compounds were evaluated in radioligand binding experiments, i.e., displacement, competition association, and washout assays, to evaluate their affinity and binding kinetic parameters. We also linked these pharmacological parameters to the compounds' chemical characteristics, and learned that separate moieties of the molecules governed target affinity and binding kinetics. Among the 29 compounds tested, 28 stood out with high affinity and a long residence time of 87 min. These findings reveal the importance of supplementing affinity data with binding kinetics at transport proteins such as hENT1.

Computational Workflow for Refining AlphaFold Models in Drug Design Using Kinetic and Thermodynamic Binding Calculations: A Case Study for the Unresolved Inactive Human Adenosine A3 Receptor

A structure-based drug design pipeline that considers both thermodynamic and kinetic binding data of ligands against a receptor will enable the computational design of improved drug molecules. For unresolved GPCR-ligand complexes, a workflow that can apply both thermodynamic and kinetic binding data in combination with alpha-fold (AF)-derived or other homology models and experimentally resolved binding modes of relevant ligands in GPCR-homologs needs to be tested. Here, as test case, we studied a congeneric set of ligands that bind to a structurally unresolved G protein-coupled receptor (GPCR), the inactive human adenosine A3 receptor (hA3R). We tested three available homology models from which two have been generated from experimental structures of hA1R or hA2AR and one model was a multistate alphafold 2 (AF2)-derived model. We applied alchemical calculations with thermodynamic integration coupled with molecular dynamics (TI/MD) simulations to calculate the experimental relative binding free energies and residence time (τ)-random accelerated MD (τ-RAMD) simulations to calculate the relative residence times (RTs) for antagonists. While the TI/MD calculations produced, for the three homology models, good Pearson correlation coefficients, correspondingly, r = 0.74, 0.62, and 0.67 and mean unsigned error (mue) values of 0.94, 1.31, and 0.81 kcal mol-1, the τ-RAMD method showed r = 0.92 and 0.52 for the first two models but failed to produce accurate results for the multistate AF2-derived model. With subsequent optimization of the AF2-derived model by reorientation of the side chain of R1735.34 located in the extracellular loop 2 (EL2) that blocked ligand's unbinding, the computational model showed r = 0.84 for kinetic data and improved performance for thermodynamic data (r = 0.81, mue = 0.56 kcal mol-1). Overall, after refining the multistate AF2 model with physics-based tools, we were able to show a strong correlation between predicted and experimental ligand relative residence times and affinities, achieving a level of accuracy comparable to an experimental structure. The computational workflow used can be applied to other receptors, helping to rank candidate drugs in a congeneric series and enabling the prioritization of leads with stronger binding affinities and longer residence times.

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.

Structure-Kinetic Relationship Studies for the Development of Long Residence Time LpxC Inhibitors

UDP-3-O-(R-3-hydroxymyristoyl)-N-acetylglucosamine deacetylase (LpxC) is a promising drug target in Gram-negative bacteria. Previously, we described a correlation between the residence time of inhibitors on Pseudomonas aeruginosa LpxC (paLpxC) and the post-antibiotic effect (PAE) caused by the inhibitors on the growth of P. aeruginosa. Given that drugs with prolonged activity following compound removal may have advantages in dosing regimens, we have explored the structure-kinetic relationship for paLpxC inhibition by analogues of the pyridone methylsulfone PF5081090 (1) originally developed by Pfizer. Several analogues have longer residence times on paLpxC than 1 (41 min) including PT913, which has a residence time of 124 min. PT913 also has a PAE of 4 h, extending the original correlation observed between residence time and PAE. Collectively, the studies provide a platform for the rational modulation of paLpxC inhibitor residence time and the potential development of antibacterial agents that cause prolonged suppression of bacterial growth.

Insight into Thermodynamic and Kinetic Profiles in Small-Molecule Optimization

Structure-activity relationships (SARs) and structure-property relationships (SPRs) have been considered the most important factors during the drug optimization process. For medicinal chemists, improvements in the potencies and druglike properties of small molecules are regarded as their major goals. Among them, the binding affinity and selectivity of small molecules on their targets are the most important indicators. In recent years, there has been growing interest in using thermodynamic and kinetic profiles to analyze ligand-receptor interactions, which could provide not only binding affinities but also detailed binding parameters for small-molecule optimization. In this perspective, we are trying to provide an insight into thermodynamic and kinetic profiles in small-molecule optimization. Through a highlight of strategies on the small-molecule optimization with specific cases, we aim to put forward the importance of structure-thermodynamic relationships (STRs) and structure-kinetic relationships (SKRs), which could provide more guidance to find safe and effective small-molecule drugs.

Selective A3 Adenosine Receptor Antagonist Radioligand for Human and Rodent Species

The A₃ adenosine receptor (A₃AR) is a target for pain, ischemia, and inflammatory disease therapy. Among the ligand tools available are selective agonists and antagonists, including radioligands, but most high-affinity non-nucleoside antagonists are limited in selectivity to primate species. We have explored the structure-activity relationship of a previously reported A₃AR antagonist DPTN 9 (N-[4-(3,5-dimethylphenyl)-5-(4-pyridyl)-1,3-thiazol-2-yl]nicotinamide) for radiolabeling, including 3-halo derivatives (3-iodo, MRS7907), and characterized 9 as a high -affinity radioligand [³H]MRS7799. A₃AR K _d values were (nM): 0.55 (human), 3.74 (mouse), and 2.80 (rat). An extended methyl acrylate (MRS8074, 19) maintained higher affinity (18.9 nM) than a 3-((5-chlorothiophen-2-yl)ethynyl) derivative 20. Compound 9 had an excellent brain distribution in rats (brain/plasma ratio ∼1). Receptor docking predicted its orthosteric site binding by engaging residues that were previously found to be essential for AR binding. Thus the new radioligand promises to be a useful species-general antagonist tracer for receptor characterization and drug discovery.

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.

In Vitro and In Silico Analysis of the Residence Time of Serotonin 5-HT7 Receptor Ligands with Arylpiperazine Structure: A Structure-Kinetics Relationship Study

During the last decade, the kinetics of drug-target interaction has received increasing attention as an important pharmacological parameter in the drug development process. Several studies have suggested that the lipophilicity of a molecule can play an important role. To date, this aspect has been studied for several G protein-coupled receptors (GPCRs) ligands but not for the 5-HT₇ receptor (5-HT₇R), a GPCR proposed as a valid therapeutic target in neurodevelopmental and neuropsychiatric disorders associated with abnormal neuronal connectivity. In this study, we report on structure-kinetics relationships of a set of arylpiperazine-based 5-HT₇R ligands. We found that it is not the overall lipophilicity of the molecule that influences drug-target interaction kinetics but rather the position of polar groups within the molecule. Next, we performed a combination of molecular docking studies and molecular dynamics simulations to gain insights into structure-kinetics relationships. These studies did not suggest specific contact patterns between the ligands and the receptor-binding site as determinants for compounds kinetics. Finally, we compared the abilities of two 5-HT₇R agonists with similar receptor-binding affinities and different residence times to stimulate the 5-HT₇R-mediated neurite outgrowth in mouse neuronal primary cultures and found that the compounds induced the effect with different timing. This study provides the first insights into the binding kinetics of arylpiperazine-based 5-HT₇R ligands that can be helpful to design new 5-HT₇R ligands with fine-tuning of the kinetic profile.

Unraveling binding mechanism and kinetics of macrocyclic Gαq protein inhibitors

G proteins represent intracellular switches that transduce signals relayed from G protein-coupled receptors. The structurally related macrocyclic depsipeptides FR900359 (FR) and YM-254890 (YM) are potent, selective inhibitors of the Gα_q protein family. We recently discovered that radiolabeled FR and YM display strongly divergent residence times, which translates into significantly longer antiasthmatic effects of FR. The present study is aimed at investigating the molecular basis for this observed disparity. Based on docking studies, we mutated amino acid residues of the Gα_q protein predicted to interact with FR or YM, and recombinantly expressed the mutated Gα_q proteins in cells in which the native Gα_q proteins had been knocked out by CRISPR-Cas9. Both radioligands showed similar association kinetics, and their binding followed a conformational selection mechanism, which was rationalized by molecular dynamics simulation studies. Several mutations of amino acid residues near the putative binding site of the "lipophilic anchors" of FR, especially those predicted to interact with the isopropyl group present in FR but not in YM, led to dramatically accelerated dissociation kinetics. Our data indicate that the long residence time of FR depends on lipophilic interactions within its binding site. The observed structure-kinetic relationships point to a complex binding mechanism of FR, which likely involves snap-lock- or dowel-like conformational changes of either ligand or protein, or both. These experimental data will be useful for the design of compounds with a desired residence time, a parameter that has now been recognized to be of utmost importance in drug development.

Evolution of the drug-target residence time model

IntroductionThe pharmacological action of a drug is linked to its affinity for a specific molecular target as quantified by in vitro equilibrium measurements. However, it is clear that for many highly effective drugs, interactions with their molecular targets do not conform to simple, equilibrium conditions in vivo and this results in a temporal discordance between pharmacokinetics and pharmacodynamics. The drug-target residence time model was developed to provide a theoretical framework with which to understand cases in which very slow dissociation of the drug-target complex in vivo results in durable PD effects even after systemic concentrations of drug have waned.Area coveredIn this article, the author provides a brief description of the drug-target residence time model and focuses on the refinements that have been made to the original model to incorporate the influences of compound rebinding in cells and pharmacokinetic properties of drug molecules.Expert opinionThere is now overwhelming evidence for the utility of the drug-target residence time model as a framework for understanding in vivo drug action. The in vitro measured residence time (τ_R) must be used in concert with equilibrium measures of drug-target affinity (e.g. IC₅₀) and with in vivo measures of pharmacokinetic half-life, to afford the researcher a powerful approach to compound optimization for clinical effect. Despite the significant use and refinement of this model, continued studies are required to better understand the dynamic interplay between residence time, target pathobiology, drug distribution and drug pharmacokinetics.

The Medicinal Chemistry of Caffeine

The purine alkaloid caffeine is the most widely consumed psychostimulant drug in the world and has multiple beneficial pharmacological activities, for example, in neurodegenerative diseases. However, despite being an extensively studied bioactive natural product, the mechanistic understanding of caffeine's pharmacological effects is incomplete. While several molecular targets of caffeine such as adenosine receptors and phosphodiesterases have been known for decades and inspired numerous medicinal chemistry programs, new protein interactions of the xanthine are continuously discovered providing potentially improved pharmacological understanding and a molecular basis for future medicinal chemistry. In this Perspective, we gather knowledge on the confirmed protein interactions, structure activity relationship, and chemical biology of caffeine on well-known and upcoming targets. The diversity of caffeine's molecular activities on receptors and enzymes, many of which are abundant in the CNS, indicates a complex interplay of several mechanisms contributing to neuroprotective effects and highlights new targets as attractive subjects for drug discovery.

Exploring the Effect of Cyclization of Histamine H1 Receptor Antagonists on Ligand Binding Kinetics

There is an increasing interest in guiding hit optimization by considering the target binding kinetics of ligands. However, compared to conventional structure-activity relationships, structure-kinetics relationships have not been as thoroughly explored, even for well-studied archetypical drug targets such as the histamine H1 receptor (H1R), a member of the family A G-protein coupled receptor. In this study, we show that the binding kinetics of H1R antagonists at the H1R is dependent on the cyclicity of both the aromatic head group and the amine moiety of H1R ligands, the chemotypes that are characteristic for the first-generation H1R antagonists. Fusing the two aromatic rings of H1R ligands into one tricyclic aromatic head group prolongs the H1R residence time for benchmark H1R ligands as well as for tailored synthetic analogues. The effect of constraining the aromatic rings and the basic amines is systematically explored, leading to a coherent series and detailed discussions of structure-kinetics relationships. This study shows that cyclicity has a pronounced effect on the binding kinetics.

Chemoenzymatic synthesis of enantiomerically enriched diprophylline and xanthinol nicotinate

A concise chemoenzymatic route toward enantiomerically enriched active pharmaceutical ingredients (API) - diprophylline and xanthinol nicotinate - is reported for the first time. The decisive step is an enantioselective lipase-mediated methanolysis of racemic chlorohydrin-synthon acetate, namely 1-chloro-3-(1,3-dimethyl-2,6-dioxo-1,2,3,6-tetrahydro-7H-purin-7-yl)propan-2-yl acetate, performed under kinetically-controlled conditions on a preparative 500 mg-scale. The best results in terms of reaction enantioselectivity (E = 14) were obtained for the enantiomers resolution performed with lipase type B from Candida antarctica immobilized on acrylic resin (CAL-B, Novozym 435) suspended in homophasic acetonitrile-methanol mixture. The elaborated biocatalytic system furnished the key chlorohydrin intermediate (in 71% ee and 38% yield), which was then smoothly converted into enantioenriched active agents: (R)-(-)-diprophylline (57% ee) and (S)-(+)-xanthinol nicotinate (65% ee). To support the assignment of absolute configurations of EKR-products as well as to confirm the stereochemical outcome of the remaining reaction steps, docking studies toward the prediction of enantiomers binding selectivity in CAL-B active site as well as the respective chemical correlations with enantiomerically enriched analytical standards obtained from commercially available (R)-(-)-epichlorohydrin, were applied. In addition, single-crystal X-ray diffraction (XRD) analyses were performed for the synthesized optically active APIs furnishing by this manner a first crystal structures of nicotinic acid salt of xanthinol.

Rhodium-catalyzed iminoiodane-mediated oxyamidation studies of 5-vinyluracil derivatives using aryl and alkyl sulfamates

The dirhodium tetraacetate-catalyzed iminoiodane-mediated reaction of 1,3-dimethyl-5-vinyluracil with phenyl sulfamate provided a high yield of 5-(1-acetyl-2-phenylsulfamoyl)ethyluracil via regioselective nucleophilic ring opening by acetate anion of the transiently formed 5-(1,2)-N-phenylsulfonylaziridine intermediate. This 1,2-oxyamidation reaction was found to be general for a variety of aryl- and alkylsulfamates as well as for various 1,3-dialkyl-5-vinyluracil derivatives. Addition of an alcohol to the reaction mixture allowed formation of the corresponding 1-alkoxy products. A selection of the substituted uracil derivatives prepared by this procedure was evaluated for cytotoxic activities in HCT-116 and HepG2 cancer cell lines and showed either no or modest activities. Further evaluation for α-glucosidase inhibition revealed that compounds 15ca and 15da were more active than acarbose, the reference inhibitor. This methodology thus allows efficient preparation of highly functionalized uracil derivatives thereby providing a synthetic route to novel compounds with potentially useful biological activities.

Structure-Kinetic Profiling of Haloperidol Analogues at the Human Dopamine D2 Receptor

Haloperidol is a typical antipsychotic drug (APD) associated with an increased risk of extrapyramidal side effects (EPSs) and hyperprolactinemia relative to atypical APDs such as clozapine. Both drugs are dopamine D₂ receptor (D₂R) antagonists, with contrasting kinetic profiles. Haloperidol displays fast association/slow dissociation at the D₂R, whereas clozapine exhibits relatively slow association/fast dissociation. Recently, we have provided evidence that slow dissociation from the D₂R predicts hyperprolactinemia, whereas fast association predicts EPS. Unfortunately, clozapine can cause severe side effects independent of its D₂R action. Our results suggest an optimal kinetic profile for D₂R antagonist APDs that avoids EPS. To begin exploring this hypothesis, we conducted a structure-kinetic relationship study of haloperidol and revealed that subtle structural modifications dramatically change binding kinetic rate constants, affording compounds with a clozapine-like kinetic profile. Thus, optimization of these kinetic parameters may allow development of novel APDs based on the haloperidol scaffold with improved side-effect profiles.

Development of Covalent Ligands for G Protein-Coupled Receptors: A Case for the Human Adenosine A3 Receptor

The development of covalent ligands for G protein-coupled receptors (GPCRs) is not a trivial process. Here, we report a streamlined workflow thereto from synthesis to validation, exemplified by the discovery of a covalent antagonist for the human adenosine A₃ receptor (hA₃AR). Based on the 1 H,3 H-pyrido[2,1- f]purine-2,4-dione scaffold, a series of ligands bearing a fluorosulfonyl warhead and a varying linker was synthesized. This series was subjected to an affinity screen, revealing compound 17b as the most potent antagonist. In addition, a nonreactive methylsulfonyl derivative 19 was developed as a reversible control compound. A series of assays, comprising time-dependent affinity determination, washout experiments, and [³⁵S]GTPγS binding assays, then validated 17b as the covalent antagonist. A combined in silico hA₃AR-homology model and site-directed mutagenesis study was performed to demonstrate that amino acid residue Y265^7.36 was the unique anchor point of the covalent interaction. This workflow might be applied to other GPCRs to guide the discovery of covalent ligands.

A live cell NanoBRET binding assay allows the study of ligand-binding kinetics to the adenosine A3 receptor

There is a growing interest in understanding the binding kinetics of compounds that bind to G protein-coupled receptors prior to progressing a lead compound into clinical trials. The widely expressed adenosine A₃ receptor (A₃AR) has been implicated in a range of diseases including immune conditions, and compounds that aim to selectively target this receptor are currently under development for arthritis. Kinetic studies at the A₃AR have been performed using a radiolabelled antagonist, but due to the kinetics of this probe, they have been carried out at 10 °C in membrane preparations. In this study, we have developed a live cell NanoBRET ligand binding assay using fluorescent A₃AR antagonists to measure kinetic parameters of labelled and unlabelled compounds at the A₃AR at physiological temperatures. The kinetic profiles of four fluorescent antagonists were determined in kinetic association assays, and it was found that XAC-ser-tyr-X-BY630 had the longest residence time (RT = 288 ± 62 min) at the A₃AR. The association and dissociation rate constants of three antagonists PSB-11, compound 5, and LUF7565 were also determined using two fluorescent ligands (XAC-ser-tyr-X-BY630 or AV039, RT = 6.8 ± 0.8 min) as the labelled probe and compared to those obtained using a radiolabelled antagonist ([³H]PSB-11, RT = 44.6 ± 3.9 min). There was close agreement in the kinetic parameters measured with AV039 and [³H]PSB-11 but significant differences to those obtained using XAC-S-ser-S-tyr-X-BY630. These data indicate that selecting a probe with the appropriate kinetics is important to accurately determine the kinetics of unlabelled ligands with markedly different kinetic profiles.

Computational modeling approaches to quantitative structure-binding kinetics relationships in drug discovery

Simple comparative correlation analyses and quantitative structure-kinetics relationship (QSKR) models highlight the interplay of kinetic rates and binding affinity as an essential feature in drug design and discovery. The choice of the molecular series, and their structural variations, used in QSKR modeling is fundamental to understanding the mechanistic implications of ligand and/or drug-target binding and/or unbinding processes. Here, we discuss the implications of linear correlations between kinetic rates and binding affinity constants and the relevance of the computational approaches to QSKR modeling.

Determination of Affinity and Residence Time of Potent Drug-Target Complexes by Label-free Biosensing

Prolonged drug-target occupancy has become increasingly important in lead optimization, and biophysical assays that measure residence time are in high demand. Here we report a practical label-free assay methodology that provides kinetic and affinity measurements suitable for most target classes without long preincubations and over comparatively short sample contact times. The method, referred to as a "chaser" assay, has been applied to three sets of unrelated kinase/inhibitor panels in order to measure the residence times, where correlation with observed efficacy was suspected. A lower throughput chaser assay measured a residence time of 3.6 days ±3.4% (95% CI) and provided single digit pM sensitivity. A higher throughput chaser methodology enabled a maximum capacity of 108 compounds in duplicate/day with an upper residence time limit of 9 h given an assay dissociation time of 34 min.

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.