Rebinding: or why drugs may act longer in vivo than expected from their in vitro target residence time

Application of Generative Artificial Intelligence in Predicting Membrane Partitioning of Drugs: Combining Denoising Diffusion Probabilistic Models and MD Simulations Reduces the Computational Cost to One-Third

The optimal interaction of drugs with plasma membranes and membranes of subcellular organelles is a prerequisite for desirable pharmacology. Importantly, for drugs targeting the transmembrane lipid-facing sites of integral membrane proteins, the relative affinity of a drug to the bilayer lipids compared to the surrounding aqueous phase affects the partitioning, access, and binding of the drug to the target site. Molecular dynamics (MD) simulations, including enhanced sampling techniques such as steered MD, umbrella sampling (US), and metadynamics, offer valuable insights into the interactions of drugs with the membrane lipids and water in atomistic detail. However, these methods are computationally prohibitive for the high-throughput screening of drug candidates. This study shows that applying denoising diffusion probabilistic models (DDPMs), a generative AI method, to US simulation data reduces the computational cost significantly. Specifically, the models used only partial (one-third) data from the US simulations and reproduced the complete potential of mean force (PMF) profiles for three FDA-approved drugs (β2-adrenergic agonists) and ∼20 biologically relevant chemicals with known experimentally characterized bilayer locations. Intriguingly, the model can predict the solvation-free energies for partitioning and crossing the bilayer, preferred bilayer locations (low-energy well), and orientations of the ligands with high accuracy. The results indicate that DDPMs can be used to characterize the complete membrane partitioning profile of drug molecules using fewer umbrella sampling simulations at select positions along the bilayer normal (z-axis), irrespective of their amphiphilic-lipophilic-cephalophilic characteristics.

From Molecular Dynamics to Taste Sensory Perception: A Comprehensive Study on the Interaction of Umami Peptides with the T1R1/T1R3-VFT Receptor

The research on the umami receptor-ligand interaction is crucial for understanding umami perception. This study integrated molecular simulations, sensory evaluation, and biosensor technology to analyze the interaction between umami peptides and the umami receptor T1R1/T1R3-VFT. Molecular dynamics simulations were used to investigate the dissociation process of seven umami peptides with the umami receptor T1R1/T1R3-VFT, and by calculating the potential mean force curve using the Jarzynski equation, it was found that the binding free energy of umami peptide is between -58.80 and -12.17 kcal/mol, which had a strong correlation with the umami intensity obtained by time intensity sensory evaluation. Through correlation analysis, the dissociation rate constants (0.0126-0.394 1/s) of umami peptides were found to have a great impact on umami perception. The faster the dissociation rate of umami peptides from receptors, the stronger the perceived intensity of the umami taste. This research aims to elucidate the relationship between the umami peptide-receptor interaction and umami perception, providing theoretical support for the exploration of umami perception mechanisms.

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

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

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.

Synthesis and single-molecule imaging reveal stereospecific enhancement of binding kinetics by the antitumour eEF1A antagonist SR-A3

Ternatin-family cyclic peptides inhibit protein synthesis by targeting the eukaryotic elongation factor-1α. A potentially related cytotoxic natural product ('A3') was isolated from Aspergillus, but only 4 of its 11 stereocentres could be assigned. Here, we synthesized SR-A3 and SS-A3-two out of 128 possible A3 epimers-and discovered that synthetic SR-A3 is indistinguishable from naturally derived A3. Relative to SS-A3, SR-A3 exhibits an enhanced residence time and rebinding kinetics, as revealed by single-molecule fluorescence imaging of elongation reactions catalysed by eukaryotic elongation factor-1α in vitro. An increased residence time-stereospecifically conferred by the unique β-hydroxyl in SR-A3-was also observed in cells. Consistent with its prolonged duration of action, thrice-weekly dosing with SR-A3 led to a reduced tumour burden and increased survival in an aggressive Myc-driven mouse lymphoma model. Our results demonstrate the potential of SR-A3 as a cancer therapeutic and exemplify an evolutionary mechanism for enhancing cyclic peptide binding kinetics via stereospecific side-chain hydroxylation.

A Scintillation Proximity Assay for Real-Time Kinetic Analysis of Chemokine–Chemokine Receptor Interactions

Chemokine receptors are extensively involved in a broad range of physiological and pathological processes, making them attractive drug targets. However, despite considerable efforts, there are very few approved drugs targeting this class of seven transmembrane domain receptors to date. In recent years, the importance of including binding kinetics in drug discovery campaigns was emphasized. Therefore, kinetic insight into chemokine–chemokine receptor interactions could help to address this issue. Moreover, it could additionally deepen our understanding of the selectivity and promiscuity of the chemokine–chemokine receptor network. Here, we describe the application, optimization and validation of a homogenous Scintillation Proximity Assay (SPA) for real-time kinetic profiling of chemokine–chemokine receptor interactions on the example of ACKR3 and CXCL12. The principle of the SPA is the detection of radioligand binding to receptors reconstituted into nanodiscs by scintillation light. No receptor modifications are required. The nanodiscs provide a native-like environment for receptors and allow for full control over bilayer composition and size. The continuous assay format enables the monitoring of binding reactions in real-time, and directly accounts for non-specific binding and potential artefacts. Minor adaptations additionally facilitate the determination of equilibrium binding metrics, making the assay a versatile tool for the study of receptor–ligand interactions.

Long Residence Time at the Vasopressin V2 Receptor Translates into Superior Inhibitory Effects in Ex Vivo and In Vivo Models of Autosomal Dominant Polycystic Kidney Disease

Prevailing strategies directing early-phase drug discovery heavily rely on equilibrium-based metrics such as affinity, which overlooks the kinetic process of a drug molecule interacting with its target. Herein, we developed a number of vasopressin V₂ receptor (V₂R) antagonists with divergent binding affinities and kinetics for autosomal dominant polycystic kidney disease (ADPKD). Surprisingly, the residence time of the V₂R antagonists, but not their affinity, was correlated with the efficacy in both ex vivo and in vivo models of ADPKD. We envision that the kinetics-directed drug candidate selection and development may have general applicability for ADPKD and other therapeutic areas as well.

Development of an intracellular quantitative assay to measure compound binding kinetics

Contemporary drug discovery typically quantifies the effect of a molecule on a biological target using the equilibrium-derived measurements of IC50, EC50, or KD. Kinetic descriptors of drug binding are frequently linked with the effectiveness of a molecule in modulating a disease phenotype; however, these parameters are yet to be fully adopted in early drug discovery. Nanoluciferase bioluminescence resonance energy transfer (NanoBRET) can be used to measure interactions between fluorophore-conjugated probes and luciferase fused target proteins. Here, we describe an intracellular NanoBRET competition assay that can be used to quantify cellular kinetic rates of compound binding to nanoluciferase-fused bromodomain and extra-terminal (BET) proteins. Comparative rates are generated using a cell-free NanoBRET assay and by utilizing orthogonal recombinant protein-based methodologies. A screen of known pan-BET inhibitors is used to demonstrate the value of this approach in the investigation of kinetic selectivity between closely related proteins.

Relative Binding Free-Energy Calculations at Lipid-Exposed Sites: Deciphering Hot Spots

Relative binding free-energy (RBFE) calculations are experiencing resurgence in the computer-aided drug design of novel small molecules due to performance gains allowed by cutting-edge molecular mechanic force fields and computer hardware. Application of RBFE to soluble proteins is becoming a routine, while recent studies outline necessary steps to successfully apply RBFE at the orthosteric site of membrane-embedded G-protein-coupled receptors (GPCRs). In this work, we apply RBFE to a congeneric series of antagonists that bind to a lipid-exposed, extra-helical site of the P2Y1 receptor. We find promising performance of RBFE, such that it may be applied in a predictive manner on drug discovery programs targeting lipid-exposed sites. Further, by the application of the microkinetic model, binding at a lipid-exposed site can be split into (1) membrane partitioning of the drug molecule followed by (2) binding at the extra-helical site. We find that RBFE can be applied to calculate the free energy of each step, allowing the uncoupling of observed binding free energy from the influence of membrane affinity. This protocol may be used to identify binding hot spots at extra-helical sites and guide drug discovery programs toward optimizing intrinsic activity at the target.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

Mapping Arginase Expression with 18F-Fluorinated Late-Generation Arginase Inhibitors Derived from Quaternary α-Amino Acids

Arginase hydrolyzes L-arginine and influences levels of polyamines and nitric oxide. Arginase overexpression is associated with inflammation and tumorigenesis. Thus, radiolabeled arginase inhibitors may be suitable PET tracers for staging arginase-related pathophysiologies. We report the synthesis and evaluation of 2 radiolabeled arginase inhibitors, 18F-FMARS and 18F-FBMARS, developed from α-substituted-2-amino-6-boronohexanoic acid derivatives. Methods: Arylboronic ester-derived precursors were radiolabeled via copper-mediated fluorodeboronation. Binding assays using arginase-expressing PC3 and LNCaP cells were performed. Autoradiography of lung sections from a guinea pig model of asthma overexpressing arginase and dynamic small-animal PET imaging with PC3-xenografted mice evaluated the radiotracers' specific binding and pharmacokinetics. Results:18F-fluorinated compounds were obtained with radiochemical yields of up to 5% (decay-corrected) and an average molar activity of 53 GBq⋅μmol-1 Cell and lung section experiments indicated specific binding that was blocked up to 75% after pretreatment with arginase inhibitors. Small-animal PET studies indicated fast clearance of the radiotracers (7.3 ± 0.6 min), arginase-mediated uptake, and a selective tumor accumulation (SUV, 3.0 ± 0.7). Conclusion: The new 18F-fluorinated arginase inhibitors have the potential to map increased arginase expression related to inflammatory and tumorigenic processes. 18F-FBMARS showed the highest arginase-mediated uptake in PET imaging and a significant difference between uptake in control and arginase-inhibited PC3 xenografted mice. These results encourage further research to examine the suitability of 18F-FBMARS for selecting patients for treatments with arginase inhibitors.

Chance Favors the Perplexed Mind: The Critical Role of Mechanistic Biochemistry in Drug Discovery

Scientific discoveries often start with an observation that does not quite make sense, within the framework of a well-established hypothesis. It is when researchers delve deeply to understand such perplexing data that established hypotheses are modified or replaced, and new and expanded knowledge of the system can be gained. This is often the case in the field of drug discovery. In this Perspective, case studies demonstrate how an understanding of perplexing data can lead to novel discoveries regarding the biological function of drug targets, or the mechanisms of compound-target interactions, that can ultimately result in new drugs entering the clinic. These case studies reinforce two interdependent themes: (1) that understanding the pathophysiological context in which drug targets function and the mechanistic details of drug-target interactions are critical to efficient and effective drug discovery and (2) that investing time and energy into following up on perplexing data can lead to novel discoveries that can drive the development of new and improved medicines.

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.

Accurate quantitative determination of affinity and binding kinetics for tight binding inhibition of xanthine oxidase

The accurate quantitative determination of affinity and binding kinetics (BK) for tight binding inhibition is extraordinary important from both the continuous optimization of compounds, particularly in developing structure-activity relationships (SAR), and the prediction of in vivo target occupancy (TO). Due to the unique properties for tight binding inhibition that the inhibitors are characterized by the ultrahigh-affinity, relatively fast association to the target enzyme combined with extremely slow dissociation of the inhibitor-enzyme binary complex, the classical steady state equilibrium methods are no longer valid. Here, we made several recommendations of how to design the optimal experiments and apply special mathematical calculation approaches to quantitatively evaluate the accurate affinity and BK as the examples of two tight binding inhibitors against the xanthine oxidase (XO), as well as compared the differences in the results calculated from the different data analytical methods and analyzed the influence of these differences on the XO engagement in human. Analysis of the results displayed that the accurate apparent dissociation constant (K_i^*,app) was 0.2 ± 0.06 nM for topiroxotstat and was 0.45 ± 0.2 nM for febuxostat; that on-rate (k_on) was (4.3 ± 1.1) × 10⁶ M^-1s^-1 for topiroxotstat and was(133.3 ± 3.5) × 10⁶ M^-1s^-1 for febuxostat, and off-rate (k_off) was (1.0±0.2) × 10^-5 s^-1 for topiroxotstat and was ≤ 0.16 × 10^-5 s^-1for febuxostat. Moreover, there were significant differences in the K_i^*,app and k_off values estimated using the appropriate specialized methods for tight binding inhibition versus classical steady state equilibrium methods, with the substantial differences of 14-fold and 32-fold reduction for topiroxostat, respectively, and of 9.6-fold and ≥ 213-fold reduction for febuxostat, while the k_on values remain the moderate differences for the two inhibitors. The obvious greater AUC of XO engagement time courses and longer durations of above 70% engagement by the appropriate specialized methods for tight binding inhibition were observed that the results display the differences of 70.1% and 88%, respectively for topiroxostat and of 38.1% and 35.0%, respectively for febuxostat in human liver cell than by classical steady state equilibrium methods. Again, our studies provide several valuable recommendations of the optimal experiment protocols and appropriate analytical approaches for accurately quantitatively assessing the affinity and BK parameters as well as demonstrate the ability of our recommended methods to generate reliable data for tight binding inhibitors against XO.

Perspectives on High-Throughput Ligand/Protein Docking With Martini MD Simulations

Molecular docking is central to rational drug design. Current docking techniques suffer, however, from limitations in protein flexibility and solvation models and by the use of simplified scoring functions. All-atom molecular dynamics simulations, on the other hand, feature a realistic representation of protein flexibility and solvent, but require knowledge of the binding site. Recently we showed that coarse-grained molecular dynamics simulations, based on the most recent version of the Martini force field, can be used to predict protein/ligand binding sites and pathways, without requiring any a priori information, and offer a level of accuracy approaching all-atom simulations. Given the excellent computational efficiency of Martini, this opens the way to high-throughput drug screening based on dynamic docking pipelines. In this opinion article, we sketch the roadmap to achieve this goal.

Adapting ADME and Pharmacokinetic Analysis to the Next Generation of Therapeutic Modalities

The development of multiple drug modalities over the past 20 years has dramatically expanded the therapeutic space for intervention in disease processes. Rather than being alternative therapeutic approaches, these modalities tend to be complimentary both in the scope of target space and the biological mechanisms harnessed for disease control. Realization of these therapeutic opportunities requires an understanding of the physiological, biochemical and biological barriers that control exposure to the drug target and resulting biological response. Consequently, successful application of ADME and PK/PD to characterization of novel therapeutics needs to consider the unique attributes conferred by the therapeutic modality and the desired and potential off-target biological responses. The discussion that follows provides examples of how barriers to exposure, and translation of exposure to efficacy can change across different modalities. Additionally, recommendations are made for ADME analysis in which biological barriers and mechanistic properties unique to specific modalities are used to focus ADME PK optimization and characterization.

Slow-binding reversible inhibitor of acetylcholinesterase with long-lasting action for prophylaxis of organophosphate poisoning

Organophosphorus (OP) compounds represent a serious health hazard worldwide. The dominant mechanism of their action results from covalent inhibition of acetylcholinesterase (AChE). Standard therapy of acute OP poisoning is partially effective. However, prophylactic administration of reversible or pseudo-irreversible AChE inhibitors before OP exposure increases the efficiency of standard therapy. The purpose of the study was to test the duration of the protective effect of a slow-binding reversible AChE inhibitor (C547) in a mouse model against acute exposure to paraoxon (POX). It was shown that the rate of inhibition of AChE by POX in vitro after pre-inhibition with C547 was several times lower than without C547. Ex vivo pre-incubation of mouse diaphragm with C547 significantly prevented the POX-induced muscle weakness. Then it was shown that pre-treatment of mice with C547 at the dose of 0.01 mg/kg significantly increased survival after poisoning by 2xLD50 POX. The duration of the pre-treatment was effective up to 96 h, whereas currently used drug for pre-exposure treatment, pyridostigmine at a dose of 0.15 mg/kg was effective less than 24 h. Thus, long-lasting slow-binding reversible AChE inhibitors can be considered as new potential drugs to increase the duration of pre-exposure treatment of OP poisoning.

Application of the Nested Enzyme-Within-Enterocyte (NEWE) Turnover Model for Predicting the Time Course of Pharmacodynamic Effects

The gut wall consists of many biological elements, including enterocytes. Rapid turnover, a prominent feature of the enterocytes, has generally been ignored in the development of enterocyte-targeting drugs, although it has a comparable rate to other kinetic rates. Here, we investigated the impact of enterocyte turnover on the pharmacodynamics of enterocyte-targeting drugs by applying a model accounting for turnover of enterocytes and target proteins. Simulations showed that the pharmacodynamics depend on enterocyte lifespan when drug-target affinity is strong and half-life of target protein is long. Interindividual variability of enterocyte lifespan, which can be amplified by disease conditions, has a substantial impact on the variability of response. However, our comprehensive literature search showed that the enterocyte turnover causes a marginal impact on currently approved enterocyte-targeting drugs due to their relatively weak target affinities. This study proposes a model-informed drug development approach for selecting enterocyte-targeting drugs and their optimal dosage regimens.

A structure-kinetic relationship study using matched molecular pair analysis

The lifetime of a binary drug–target complex is increasingly acknowledged as an important parameter for drug efficacy and safety. With a better understanding of binding kinetics and better knowledge about kinetic parameter optimization, intentionally induced prolongation of the drug–target residence time through structural changes of the ligand could become feasible. In this study we assembled datasets from 21 publications and the K4DD (Kinetic for Drug Discovery) database to conduct large scale data analysis. This resulted in 3812 small molecules annotated to 78 different targets from five protein classes (GPCRs: 273, kinases: 3238, other enzymes: 240, HSPs: 160, ion channels: 45). Performing matched molecular pair (MMP) analysis to further investigate the structure–kinetic relationship (SKR) in this data collection allowed us to identify a fundamental contribution of a ligand's polarity to its association rate, and in selected cases, also to its dissociation rate. However, we furthermore observed that the destabilization of the transition state introduced by increased polarity is often accompanied by simultaneous destabilization of the ground state resulting in an unaffected or even worsened residence time. Supported by a set of case studies, we provide concepts on how to alter ligands in ways to trigger on-rates, off-rates, or both.

A large-scale study employing matched molecular pair (MMP) analysis to uncover the contribution of a compound's polarity to its association and dissociation rates.

Quantifying CDK inhibitor selectivity in live cells

Concerted multidisciplinary efforts have led to the development of Cyclin-Dependent Kinase inhibitors (CDKi’s) as small molecule drugs and chemical probes of intracellular CDK function. However, conflicting data has been reported on the inhibitory potency of CDKi’s and a systematic characterization of affinity and selectivity against intracellular CDKs is lacking. We have developed a panel of cell-permeable energy transfer probes to quantify target occupancy for all 21 human CDKs in live cells, and present a comprehensive evaluation of intracellular isozyme potency and selectivity for a collection of 46 clinically-advanced CDKi’s and tool molecules. We observed unexpected intracellular activity profiles for a number of CDKi’s, offering avenues for repurposing of highly potent molecules as probes for previously unreported targets. Overall, we provide a broadly applicable method for evaluating the selectivity of CDK inhibitors in living cells, and present a refined set of tool molecules to study CDK function.

Cyclin-dependent kinase (CDK) inhibitors are widely used both in the clinic and for basic research aimed at dissecting the specific cellular functions of specific CDKs. Here, the authors report the development of a panel of fluorescent reporter probes and provide a comprehensive profile of the inhibitory activity of several CDK inhibitors towards all 21 CDKs in living cells.

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.

[Ligand-receptor binding kinetics in drug design]

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

Drug-Target Residence Time Affects in Vivo Target Occupancy through Multiple Pathways

The drug discovery and development process is greatly hampered by difficulties in translating in vitro potency to in vivo efficacy. Recent studies suggest that the long-neglected drug-target residence time parameter complements classical drug affinity parameters (K _I, K _d, IC₅₀, or EC₅₀) and is a better predictor of in vivo efficacy. Compounds with a long drug-target residence time are often more efficacious in vivo. The impact, however, of the drug-target residence time on in vivo efficacy remains controversial due to difficulties in experimentally determining the in vivo target occupancy during drug treatment. To tackle this problem, an in vivo displacement assay was developed using soluble epoxide hydrolase as a biological model. In this report, we experimentally demonstrated that drug-target residence time affects the duration of in vivo drug-target binding. In addition, the drug-target residence time plays an important role in modulating the rate of drug metabolism which also affects the efficacy of the drug.

Drug-Target Association Kinetics in Drug Discovery

The important role of ligand-receptor binding kinetics in drug design and discovery is increasingly recognized by the drug research community. Over the past decade, accumulating evidence has shown that optimizing the ligand's dissociation rate constant can lead to desirable duration of in vivo target occupancy and, hence, improved pharmacodynamic properties. However, the association rate constant as a pharmacological principle remains less investigated, whereas it can play an equally important role in the selection of drug candidates. This review provides a compilation and discussion of otherwise scarce and dispersed information on this topic, bringing to light the importance of drug-target association in kinetics-directed drug design and discovery.

Does the Lipid Bilayer Orchestrate Access and Binding of Ligands to Transmembrane Orthosteric/Allosteric Sites of G Protein-Coupled Receptors?

The ligand-binding sites of many G protein-coupled receptors (GPCRs) are situated around and deeply embedded within the central pocket formed by their seven transmembrane-spanning α-helical domains. Generally, these binding sites are assumed accessible to endogenous ligands from the aqueous phase. Recent advances in the structural biology of GPCRs, along with biophysical and computational studies, suggest that amphiphilic and lipophilic molecules may gain access to these receptors by first partitioning into the membrane and then reaching the binding site via lateral diffusion through the lipid bilayer. In addition, several crystal structures of class A and class B GPCRs bound to their ligands offer unprecedented details on the existence of lipid-facing allosteric binding sites outside the transmembrane helices that can only be reached via lipid pathways. The highly organized structure of the lipid bilayer may direct lipophilic or amphiphilic drugs to a specific depth within the bilayer, changing local concentration of the drug near the binding site and affecting its binding kinetics. Additionally, the constraints of the lipid bilayer, including its composition and biophysical properties, may play a critical role in "pre-organizing" ligand molecules in an optimal orientation and conformation to facilitate receptor binding. Despite its clear involvement in molecular recognition processes, the critical role of the membrane in binding ligands to lipid-exposed transmembrane binding sites remains poorly understood and warrants comprehensive investigation. Understanding the mechanistic basis of the structure-membrane interaction relationship of drugs will not only provide useful insights about receptor binding kinetics but will also enhance our ability to take advantage of the apparent membrane contributions when designing drugs that target transmembrane proteins with improved efficacy and safety. In this minireview, we summarize recent structural and computational studies on membrane contributions to binding processes, elucidating both lipid pathways of ligand access and binding mechanisms for several orthosteric and allosteric ligands of class A and class B GPCRs.

Binding kinetics of ligands acting at GPCRs

The influence of drug-receptor binding kinetics has often been overlooked during the development of new therapeutics that target G protein-coupled receptors (GPCRs). Over the last decade there has been a growing understanding that an in-depth knowledge of binding kinetics at GPCRs is required to successfully target this class of proteins. Ligand binding to a GPCR is often not a simple single step process with ligand freely diffusing in solution. This review will discuss the experiments and equations that are commonly used to measure binding kinetics and how factors such as allosteric regulation, rebinding and ligand interaction with the plasma membrane may influence these measurements. We will then consider the molecular characteristics of a ligand and if these can be linked to association and dissociation rates.

Link between a high k on for drug binding and a fast clinical action: to be or not to be?

Review articles on binding kinetics essentially focus on drugs that dissociate slowly from their target since this is required for the successful treatment of many pathophysiological conditions. Recently, the therapeutic benefit of a high k on (i.e. the second order association rate constant) has also been linked to fast association and to a fast clinical action. Other studies, however, called this assertion into question since additional factors, like the dosing paradigm and the binding mechanism, are important as well. The still ongoing reticence about integrating binding kinetics in lead optimization programs motivated us to critically review the link between the drug's kinetic rate constants and their in vitro and in vivo target occupancy profile, with special focus on k on. The presented simulations tally with a positive link between a drug's effective/observed association rate (which is quite easy to determine in vitro) and the swiftness of its clinical action. On the other hand, the simulations show that the k on-concept should not be confounded with the effective association process since increasing this parameter only enhances the drug's in vitro and in vivo association under certain conditions: the binding mechanism should be suitable, rebinding (and thus the factors within the target's micro-environment that favour this mechanism) should not be too prominent and the dosage should not be kept in par with the drug's affinity. Otherwise, increasing k on could be ineffective or even be counter-productive.

KinITC-One Method Supports both Thermodynamic and Kinetic SARs as Exemplified on FimH Antagonists

Affinity data, such as dissociation constants (K_D ) or inhibitory concentrations (IC₅₀ ), are widely used in drug discovery. However, these parameters describe an equilibrium state, which is often not established in vivo due to pharmacokinetic effects and they are therefore not necessarily sufficient for evaluating drug efficacy. More accurate indicators for pharmacological activity are the kinetics of binding processes, as they shed light on the rate of formation of protein-ligand complexes and their half-life. Nonetheless, although highly desirable for medicinal chemistry programs, studies on structure-kinetic relationships (SKR) are still rare. With the recently introduced analytical tool kinITC this situation may change, since not only thermodynamic but also kinetic information of the binding process can be deduced from isothermal titration calorimetry (ITC) experiments. Using kinITC, ITC data of 29 mannosides binding to the bacterial adhesin FimH were re-analyzed to make their binding kinetics accessible. To validate these kinetic data, surface plasmon resonance (SPR) experiments were conducted. The kinetic analysis by kinITC revealed that the nanomolar affinities of the FimH antagonists arise from both (i) an optimized interaction between protein and ligand in the bound state (reduced off-rate constant k_off ) and (ii) a stabilization of the transition state or a destabilization of the unbound state (increased on-rate constant k_on ). Based on congeneric ligand modifications and structural input from co-crystal structures, a strong relationship between the formed hydrogen-bond network and k_off could be concluded, whereas electrostatic interactions and conformational restrictions upon binding were found to have mainly an impact on k_on .

Autoregulation of Acetylcholine Release and Micro-Pharmacodynamic Mechanisms at Neuromuscular Junction: Selective Acetylcholinesterase Inhibitors for Therapy of Myasthenic Syndromes

Neuromuscular junctions (NMJs) are directly involved into such indispensable to life processes as respiration and locomotion. However, motor nerve forms only one synaptic contact at each muscle fiber. This unique configuration requires specific properties and constrains to be effective. The very high density of acetylcholine receptors (AChRs) of muscle type in synaptic cleft and an excess of acetylcholine (ACh) released under physiological conditions make this synapse extremely reliable. Nevertheless, under pathological conditions such as myasthenia gravis and congenital myasthenic syndromes, the safety factor can be markedly reduced. Drugs used for short-term symptomatic therapy of these pathological states, cause partial inhibition of cholinesterases (ChEs). These enzymes catalyze the hydrolysis of ACh, thus terminate its action on AChRs. Extension of the lifetime of ACh molecules compensates muscular AChRs abnormalities and, consequently, rescues muscle contractions. In this mini review, we will first outline the functional organization of the NMJ, and then, consider the concept of the safety factor and how it may be changed. This will be followed by a look at autoregulation of ACh release that influences the safety factor of NMJs. Finally, we will consider the morphological features of NMJs as a putative reserve to increase effectiveness of pathological muscle weakness therapy by ChEs inhibitors due to opportunity to use micro-pharmacodynamic mechanisms.

Distinct in vivo target occupancy by bivalent- and induced-fit-like binding drugs

BACKGROUND AND PURPOSE

Optimal drug therapy often requires long-lasting target occupancy While this attribute was usually linked to the drug's pharmacokinetic properties, the dissociation rate is now increasingly recognized to contribute as well. Nearly all the earlier pharmacokinetic-pharmacodynamic (PK-PD) simulations encompassed single-step binding drugs and focused on koff . However, 'micro'-PK mechanisms and more complex binding mechanisms like bivalent- and induced-fit binding may contribute as well. Corresponding binding models are presently explored.

EXPERIMENTAL APPROACH

We compared the 24 h in vivo occupancy over time profiles of prototype bivalent- and induced-fit-like binding drugs (A and B) after one or repeated daily dosings, both without and with rebinding. Special attention was focused on the effect of each of the microscopic rate constants on the occupancy profiles and on the metrics to represent those profiles.

KEY RESULTS

Although both models can be represented by the same mathematical formulation, drugs A and B display quite different occupancy profiles, even though they have the same potency. These differences can be attributed to the different effects of their microscopic rate constants on their composite koff and also on their susceptibility to experience rebinding. This also affects how the occupancy profiles of bivalent- and induced-fit-like binders progress when repeating the dosings and by changing the dosage.

CONCLUSIONS AND IMPLICATIONS

Closer attention should be paid to more complex binding models in PK-PD simulations. This may help pharmacologists and medicinal chemists to improve the translation of in vitro kinetic measurements from preclinical screening programmes into clinical efficiency.

Interactions between β2-Adrenoceptor Ligands and Membrane: Atomic-Level Insights from Magic-Angle Spinning NMR

To understand the relationship between structural properties of the β2-adrenoceptor ligands and their interactions with membranes, we have investigated the location and distribution of five β2 agonists with distinct clinical durations and onsets of action (indacaterol, two indacaterol analogues, salmeterol and formoterol) in monounsaturated model membranes using magic angle spinning NMR to measure these interactions through both ¹H nuclear Overhauser enhancement (NOE) and paramagnetic relaxation enhancement (PRE) techniques. The hydrophilic aromatic groups of all five β2 agonists show maximum distribution in the lipid/water interface, but distinct location and dynamic behavior were observed for the lipophilic aromatic rings. Our study elucidates at atomic level that the hydrophobicity and substitution geometry of lipophilic groups play important roles in compound-lipid interactions.

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.

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.

Use of a new proximity assay (NanoBRET) to investigate the ligand-binding characteristics of three fluorescent ligands to the human β1-adrenoceptor expressed in HEK-293 cells

Previous research has indicated that allosteric interactions across the dimer interface of β₁‐adrenoceptors may be responsible for a secondary low affinity binding conformation. Here we have investigated the potential for probe dependence, in the determination of antagonist pK_i values at the human β₁‐adenoceptor, which may result from such allosterism interactions. Three fluorescent β₁‐adrenoceptor ligands were used to investigate this using bioluminescence energy transfer (BRET) between the receptor‐bound fluorescent ligand and the N‐terminal NanoLuc tag of a human β₁‐adrenoceptor expressed in HEK 293 cells (NanoBRET). This proximity assay showed high‐affinity‐specific binding to the NanoLuc‐ β₁‐adrenoceptor with each of the three fluorescent ligands yielding K_D values of 87.1 ± 10 nmol/L (n = 8), 38.1 ± 12 nmol/L (n = 7), 13.4 ± 2 nmol/L (n = 14) for propranolol‐Peg8‐BY630, propranolol‐ β(Ala‐Ala)‐BY630 and CGP‐12177‐TMR, respectively. Parallel radioligand‐binding studies with ³H‐CGP12177 and TIRF microscopy, to monitor NanoLuc bioluminescence, confirmed a high cell surface expression of the NanoLuc‐ β₁‐adrenoceptor in HEK 293 cells (circa 1500 fmol.mg protein⁻¹). Following a 1 h incubation with fluorescent ligands and β₁‐adrenoceptor competing antagonists, there were significant differences (P < 0.001) in the pK_i values obtained for CGP20712a and CGP 12177 with the different fluorescent ligands and ³H‐CGP 12177. However, increasing the incubation time to 2 h removed these significant differences. The data obtained show that the NanoBRET assay can be applied successfully to study ligand‐receptor interactions at the human β₁‐adrenoceptor. However, the study also emphasizes the importance of ensuring that both the fluorescent and competing ligands are in true equilibrium before interpretations regarding probe dependence can be made.

Effects of target binding kinetics on in vivo drug efficacy: koff , kon and rebinding

BACKGROUND AND PURPOSE

Optimal drug therapy often requires continuing high levels of target occupancy. Besides the traditional pharmacokinetic contribution, target binding kinetics is increasingly considered to play an important role as well. While most attention has been focused on the dissociation rate of the complex, recent reports expressed doubt about the unreserved translatability of this pharmacodynamic property into clinical efficacy. 'Micro'-pharmacokinetic mechanisms like drug rebinding and partitioning into the cell membrane may constitute a potential fix.

EXPERIMENTAL APPROACH

Simulations were based on solving differential equations.

KEY RESULTS

Based on a selected range of association and dissociation rate constants, kon and koff , and rebinding potencies of the drugs as variables, their effects on the temporal in vivo occupancy profile of their targets, after one or multiple repetitive dosings, have here been simulated.

CONCLUSIONS AND IMPLICATIONS

Most strikingly, the simulations show that, when rebinding is also taken into account, increasing kon may produce closely the same outcome as decreasing koff when dosing is performed in accordance with the therapeutically most relevant constant [Lmax ]/KD ratio paradigm. Also, under certain conditions, rebinding may produce closely the same outcome as invoking slow diffusion of the drug between the plasma compartment and a target-containing 'effect' compartment. Although the present simulations should only be regarded as a 'proof of principle', these findings may help pharmacologists and medicinal chemists to devise ex vivo and in vitro binding kinetic assays that are more relevant and translatable to in vivo settings.

Uncoupling the Structure-Activity Relationships of β2 Adrenergic Receptor Ligands from Membrane Binding

Ligand binding to membrane proteins may be significantly influenced by the interaction of ligands with the membrane. In particular, the microscopic ligand concentration within the membrane surface solvation layer may exceed that in bulk solvent, resulting in overestimation of the intrinsic protein-ligand binding contribution to the apparent/measured affinity. Using published binding data for a set of small molecules with the β2 adrenergic receptor, we demonstrate that deconvolution of membrane and protein binding contributions allows for improved structure-activity relationship analysis and structure-based drug design. Molecular dynamics simulations of ligand bound membrane protein complexes were used to validate binding poses, allowing analysis of key interactions and binding site solvation to develop structure-activity relationships of β2 ligand binding. The resulting relationships are consistent with intrinsic binding affinity (corrected for membrane interaction). The successful structure-based design of ligands targeting membrane proteins may require an assessment of membrane affinity to uncouple protein binding from membrane interactions.

Mechanism of Action of Thymosinα1: Does It Interact with Membrane by Recognition of Exposed Phosphatidylserine on Cell Surface? A Structural Approach

Thymosinα1 is a peptidic hormone with pleiotropic activity, which is used in the therapy of several diseases. It is unstructured in water solution and interacts with negative regions of micelles and vesicles assuming two tracts of helical conformation with a structural flexible break in between. The studies of the interaction of Thymosinα1 with micelles of mixed dipalmitoylphosphatidylcholine and sodium dodecylsulfate and vesicles with mixed dipalmitoylphosphatidylcholine/dipalmitoylphosphatidylserine, the latter the negative component of the membranes, by (1)H and natural abundance (15)N NMR are herewith reported, reviewed, and discussed. The results indicate that the preferred interactions are those where the surface is negatively charged due to sodium dodecylsulfate or due to the presence of dipalmitoylphosphatidylserine exposed on the surface. In fact the unbalance of dipalmitoylphosphatidylserine on the cellular surface is an important phenomenon present in pathological conditions of cells. Moreover, the direct interaction of Thymosinα1 with K562 cells presenting an overexposure of phosphatidylserine as a consequence of resveratrol-induced apoptosis was carried out.

The drug-target residence time model: a 10-year retrospective

The drug-target residence time model was first introduced in 2006 and has been broadly adopted across the chemical biology, biotechnology and pharmaceutical communities. While traditional in vitro methods view drug-target interactions exclusively in terms of equilibrium affinity, the residence time model takes into account the conformational dynamics of target macromolecules that affect drug binding and dissociation. The key tenet of this model is that the lifetime (or residence time) of the binary drug-target complex, and not the binding affinity per se, dictates much of the in vivo pharmacological activity. Here, this model is revisited and key applications of it over the past 10 years are highlighted.

Target engagement and drug residence time can be observed in living cells with BRET

The therapeutic action of drugs is predicated on their physical engagement with cellular targets. Here we describe a broadly applicable method using bioluminescence resonance energy transfer (BRET) to reveal the binding characteristics of a drug with selected targets within intact cells. Cell-permeable fluorescent tracers are used in a competitive binding format to quantify drug engagement with the target proteins fused to Nanoluc luciferase. The approach enabled us to profile isozyme-specific engagement and binding kinetics for a panel of histone deacetylase (HDAC) inhibitors. Our analysis was directed particularly to the clinically approved prodrug FK228 (Istodax/Romidepsin) because of its unique and largely unexplained mechanism of sustained intracellular action. Analysis of the binding kinetics by BRET revealed remarkably long intracellular residence times for FK228 at HDAC1, explaining the protracted intracellular behaviour of this prodrug. Our results demonstrate a novel application of BRET for assessing target engagement within the complex milieu of the intracellular environment.

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.