Molecular Basis for the Long Duration of Action and Kinetic Selectivity of Tiotropium for the Muscarinic M3 Receptor

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.

Structural basis for the ligand recognition and G protein subtype selectivity of kisspeptin receptor

Kisspeptin receptor (KISS1R), belonging to the class A peptide-GPCR family, plays a key role in the regulation of reproductive physiology after stimulation by kisspeptin and is regarded as an attractive drug target for reproductive diseases. Here, we demonstrated that KISS1R can couple to the Gi/o pathway besides the well-known Gq/11 pathway. We further resolved the cryo-electron microscopy (cryo-EM) structure of KISS1R-Gq and KISS1R-Gi complexes bound to the synthetic agonist TAK448 and structure of KISS1R-Gq complex bound to the endogenous agonist KP54. The high-resolution structures provided clear insights into mechanism of KISS1R recognition by its ligand and can facilitate the design of targeted drugs with high affinity to improve treatment effects. Moreover, the structural and functional analyses indicated that conformational differences in the extracellular loops (ECLs), intracellular loops (ICLs) of the receptor, and the "wavy hook" of the Gα subunit may account for the specificity of G protein coupling for KISS1R signaling.

Discovery of a brain penetrant small molecule antagonist targeting LPA1 receptors to reduce neuroinflammation and promote remyelination in multiple sclerosis

Multiple sclerosis (MS) is a chronic neurological disease characterized by inflammatory demyelination that disrupts neuronal transmission resulting in neurodegeneration progressive disability. While current treatments focus on immunosuppression to limit inflammation and further myelin loss, no approved therapies effectively promote remyelination to mitigate the progressive disability associated with chronic demyelination. Lysophosphatidic acid (LPA) is a pro-inflammatory lipid that is upregulated in MS patient plasma and cerebrospinal fluid (CSF). LPA activates the LPA1 receptor, resulting in elevated CNS cytokine and chemokine levels, infiltration of immune cells, and microglial/astrocyte activation. This results in a neuroinflammatory response leading to demyelination and suppressed remyelination. A medicinal chemistry effort identified PIPE-791, an oral, brain-penetrant, LPA1 antagonist. PIPE-791 was characterized in vitro and in vivo and was found to be a potent, selective LPA1 antagonist with slow receptor off-rate kinetics. In vitro, PIPE-791 induced OPC differentiation and promoted remyelination following a demyelinating insult. PIPE-791 further mitigated the macrophage-mediated inhibition of OPC differentiation and inhibited microglial and fibroblast activation. In vivo, the compound readily crossed the blood-brain barrier and blocked LPA1 in the CNS after oral dosing. Direct dosing of PIPE-791 in vivo increased oligodendrocyte number, and in the mouse experimental autoimmune encephalomyelitis (EAE) model of MS, we observed that PIPE-791 promoted myelination, reduced neuroinflammation, and restored visual evoked potential latencies (VEP). These findings support targeting LPA1 for remyelination and encourage development of PIPE-791 for treating MS patients with advantages not seen with current immunosuppressive disease modifying therapies.

Analysing the effect caused by increasing the molecular volume in M1-AChR receptor agonists and antagonists: a structural and computational study

M1 muscarinic acetylcholine receptor (M1-AChR), a member of the G protein-coupled receptors (GPCR) family, plays a crucial role in learning and memory, making it an important drug target for Alzheimer's disease (AD) and schizophrenia. M1-AChR activation and deactivation have shown modifying effects in AD and PD preclinical models, respectively. However, understanding the pharmacology associated with M1-AChR activation or deactivation is complex, because of the low selectivity among muscarinic subtypes, hampering their therapeutic applications. In this regard, we constructed two quantitative structure-activity relationship (QSAR) models, one for M1-AChR agonists (total and partial), and the other for the antagonists. The binding mode of 59 structurally different compounds, including agonists and antagonists with experimental binding affinity values (pKi), were analyzed employing computational molecular docking over different structures of M1-AChR. Furthermore, we considered the interaction energy (Einter), the number of rotatable bonds (NRB), and lipophilicity (ilogP) for the construction of the QSAR model for agonists (R2 = 89.64, QLMO2 = 78, and Qext2 = 79.1). For the QSAR model of antagonists (R2 = 88.44, QLMO2 = 82, and Qext2 = 78.1) we considered the Einter, the fraction of sp3 carbons fCsp3, and lipophilicity (MlogP). Our results suggest that the ligand volume is a determinant to establish its biological activity (agonist or antagonist), causing changes in binding energy, and determining the affinity for M1-AChR.

Xanomeline displays concomitant orthosteric and allosteric binding modes at the M4 mAChR

The M4 muscarinic acetylcholine receptor (M4 mAChR) has emerged as a drug target of high therapeutic interest due to its expression in regions of the brain involved in the regulation of psychosis, cognition, and addiction. The mAChR agonist, xanomeline, has provided significant improvement in the Positive and Negative Symptom Scale (PANSS) scores in a Phase II clinical trial for the treatment of patients suffering from schizophrenia. Here we report the active state cryo-EM structure of xanomeline bound to the human M4 mAChR in complex with the heterotrimeric Gi1 transducer protein. Unexpectedly, two molecules of xanomeline were found to concomitantly bind to the monomeric M4 mAChR, with one molecule bound in the orthosteric (acetylcholine-binding) site and a second molecule in an extracellular vestibular allosteric site. Molecular dynamic simulations supports the structural findings, and pharmacological validation confirmed that xanomeline acts as a dual orthosteric and allosteric ligand at the human M4 mAChR. These findings provide a basis for further understanding xanomeline's complex pharmacology and highlight the myriad of ways through which clinically relevant ligands can bind to and regulate GPCRs.

Role of Conserved Tyrosine Lid Residues in the Activation of the M2 Muscarinic Acetylcholine Receptor

The development of subtype selective small molecule drugs for the muscarinic acetylcholine receptor (mAChR) family has been challenging. The design of more selective ligands can be improved by understanding the structure and function of key amino acid residues that line ligand binding sites. Here we study the role of three conserved key tyrosine residues [Y1043.33, Y4036.51, and Y4267.39 (Ballesteros and Weinstein numbers in superscript)] at the human M2 mAChR, located at the interface between the orthosteric and allosteric binding sites of the receptor. We specifically focused on the role of the three tyrosine hydroxyl groups in the transition between the inactive and active conformations of the receptor by making phenylalanine point mutants. Single-point mutation at either of the three positions was sufficient to reduce the affinity of agonists by ∼100-fold for the M2 mAChR, whereas the affinity of antagonists remained largely unaffected. In contrast, neither of the mutations affected the efficacy of orthosteric agonists. When mutations were combined into double and triple M2 mAChR mutants, the affinity of antagonists was reduced by more than 100-fold compared with the wild-type M2 receptor. In contrast, the affinity of allosteric modulators, either negative or positive, was retained at all single and multiple mutations, but the degree of allosteric effect exerted on the endogenous ligand acetylcholine was affected at all mutants containing Y4267.39F. These findings will provide insights to consider when designing future mAChR ligands. SIGNIFICANCE STATEMENT: Structural studies demonstrated that three tyrosine residues between the orthosteric and allosteric sites of the M2 muscarinic acetylcholine receptor (mAChR) had different hydrogen bonding networks in the inactive and active conformations. The role of hydroxyl groups of the tyrosine residues on orthosteric and allosteric ligand pharmacology was unknown. We found that hydroxyl groups of the tyrosine residues differentially affected the molecular pharmacology of orthosteric and allosteric ligands. These results provide insights to consider when designing future mAChR ligands.

Investigating the Unbinding of Muscarinic Antagonists from the Muscarinic 3 Receptor

Patient symptom relief is often heavily influenced by the residence time of the inhibitor-target complex. For the human muscarinic receptor 3 (hMR3), tiotropium is a long-acting bronchodilator used in conditions such as asthma or chronic obstructive pulmonary disease (COPD). The mechanistic insights into this inhibitor remain unclear; specifically, the elucidation of the main factors determining the unbinding rates could help develop the next generation of antimuscarinic agents. Using our novel unbinding algorithm, we were able to investigate ligand dissociation from hMR3. The unbinding paths of tiotropium and two of its analogues, N-methylscopolamin and homatropine methylbromide, show a consistent qualitative mechanism and allow us to identify the structural bottleneck of the process. Furthermore, our machine learning-based analysis identified key roles of the ECL2/TM5 junction involved in the transition state. Additionally, our results point to relevant changes at the intracellular end of the TM6 helix leading to the ICL3 kinase domain, highlighting the closest residue L482. This residue is located right between two main protein binding sites involved in signal transduction for hMR3's activation and regulation. We also highlight key pharmacophores of tiotropium that play determining roles in the unbinding kinetics and could aid toward drug design and lead optimization.

Structural basis of efficacy-driven ligand selectivity at GPCRs

A drug's selectivity for target receptors is essential to its therapeutic utility, but achieving selectivity between similar receptors is challenging. The serendipitous discovery of ligands that stimulate target receptors more strongly than closely related receptors, despite binding with similar affinities, suggests a solution. The molecular mechanism of such 'efficacy-driven selectivity' has remained unclear, however, hindering design of such ligands. Here, using atomic-level simulations, we reveal the structural basis for the efficacy-driven selectivity of a long-studied clinical drug candidate, xanomeline, between closely related muscarinic acetylcholine receptors (mAChRs). Xanomeline's binding mode is similar across mAChRs in their inactive states but differs between mAChRs in their active states, with divergent effects on active-state stability. We validate this mechanism experimentally and use it to design ligands with altered efficacy-driven selectivity. Our results suggest strategies for the rational design of ligands that achieve efficacy-driven selectivity for many pharmaceutically important G-protein-coupled receptors.

Unbinding Kinetics of Muscarinic M3 Receptor Antagonists Explained by Metadynamics Simulations

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

Exploring the kinetic selectivity of drugs targeting the β1 -adrenoceptor

In this study, we report the β1 -adrenoceptor binding kinetics of several clinically relevant β1/2 -adrenoceptor (β1/2 AR) agonists and antagonists. [3 H]-DHA was used to label CHO-β1 AR for binding studies. The kinetics of ligand binding was assessed using a competition association binding method. Ligand physicochemical properties, including logD7.4 and the immobilized artificial membrane partition coefficient (KIAM ), were assessed using column-based methods. Protein Data Bank (PDB) structures and hydrophobic and electrostatic surface maps were constructed in PyMOL. We demonstrate that the hydrophobic properties of a molecule directly affect its kinetic association rate (kon ) and affinity for the β1 AR. In contrast to our findings at the β2 -adrenoceptor, KIAM , reflecting both hydrophobic and electrostatic interactions of the drug with the charged surface of biological membranes, was no better predictor than simple hydrophobicity measurements such as clogP or logD7.4 , at predicting association rate. Bisoprolol proved kinetically selective for the β1 AR subtype, dissociating 50 times slower and partly explaining its higher measured affinity for the β1 AR. We speculate that the association of positively charged ligands at the β1 AR is curtailed somewhat by its predominantly neutral/positive charged extracellular surface. Consequently, hydrophobic interactions in the ligand-binding pocket dominate the kinetics of ligand binding. In comparison at the β2 AR, a combination of hydrophobicity and negative charge attracts basic, positively charged ligands to the receptor's surface promoting the kinetics of ligand binding. Additionally, we reveal the potential role kinetics plays in the on-target and off-target pharmacology of clinically used β-blockers.

The Impact of the Secondary Binding Pocket on the Pharmacology of Class A GPCRs

G-protein coupled receptors (GPCRs) are considered important therapeutic targets due to their pathophysiological significance and pharmacological relevance. Class A receptors represent the largest group of GPCRs that gives the highest number of validated drug targets. Endogenous ligands bind to the orthosteric binding pocket (OBP) embedded in the intrahelical space of the receptor. During the last 10 years, however, it has been turned out that in many receptors there is secondary binding pocket (SBP) located in the extracellular vestibule that is much less conserved. In some cases, it serves as a stable allosteric site harbouring allosteric ligands that modulate the pharmacology of orthosteric binders. In other cases it is used by bitopic compounds occupying both the OBP and SBP. In these terms, SBP binding moieties might influence the pharmacology of the bitopic ligands. Together with others, our research group showed that SBP binders contribute significantly to the affinity, selectivity, functional activity, functional selectivity and binding kinetics of bitopic ligands. Based on these observations we developed a structure-based protocol for designing bitopic compounds with desired pharmacological profile.

Developing inhaled drugs for respiratory diseases: A medicinal chemistry perspective

Despite the devastating impact of many lung diseases on human health, there is still a significant unmet medical need in respiratory diseases, for which inhaled delivery represents a crucial strategy. More guidance on how to design and carry out multidisciplinary inhaled projects is needed. When designing inhaled drugs, the medicinal chemist must carefully balance the physicochemical properties of the molecule to achieve optimal target engagement in the lung. Although the medicinal chemistry strategy is unique for each project, and will change depending on multiple factors, such as the disease, target, systemic risk, delivery device, and formulation, general guidelines aiding inhaled drug design can be applied and are summarised in this review.

Guanidine Derivatives: How Simple Structural Modification of Histamine H3R Antagonists Has Led to the Discovery of Potent Muscarinic M2R/M4R Antagonists

This article describes the discovery of novel potent muscarinic receptor antagonists identified during a search for more active histamine H₃ receptor (H₃R) ligands. The idea was to replace the flexible seven methylene linker with a semirigid 1,4-cyclohexylene or p-phenylene substituted group of the previously described histamine H₃R antagonists ADS1017 and ADS1020. These simple structural modifications of the histamine H₃R antagonist led to the emergence of additional pharmacological effects, some of which unexpectedly showed strong antagonist potency at muscarinic receptors. This paper reports the routes of synthesis and pharmacological characterization of guanidine derivatives, a novel chemotype of muscarinic receptor antagonists binding to the human muscarinic M₂ and M₄ receptors (hM₂R and hM₄R, respectively) in nanomolar concentration ranges. The affinities of the newly synthesized ADS10227 (1-{4-{4-{[4-(phenoxymethyl)cyclohexyl]methyl}piperazin-1-yl}but-1-yl}-1-(benzyl)guanidine) at hM₂R and hM₄R were 2.8 nM and 5.1 nM, respectively.

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.

Analyzing GPCR-Ligand Interactions with the Fragment Molecular Orbital (FMO) Method

G-protein-coupled receptors (GPCRs) have enormous physiological and biomedical importance, and therefore it is not surprising that they are the targets of many prescribed drugs. Further progress in GPCR drug discovery is highly dependent on the availability of protein structural information. However, the ability of X-ray crystallography to guide the drug discovery process for GPCR targets is limited by the availability of accurate tools to explore receptor-ligand interactions. Visual inspection and molecular mechanics approaches cannot explain the full complexity of molecular interactions. Quantum mechanics (QM) approaches are often too computationally expensive to be of practical use in time-sensitive situations, but the fragment molecular orbital (FMO) method offers an excellent solution that combines accuracy, speed, and the ability to reveal key interactions that would otherwise be hard to detect. Integration of GPCR crystallography or homology modelling with FMO reveals atomistic details of the individual contributions of each residue and water molecule toward ligand binding, including an analysis of their chemical nature. Such information is essential for an efficient structure-based drug design (SBDD) process. In this chapter, we describe how to use FMO in the characterization of GPCR-ligand interactions.

An optimized acetylcholine sensor for monitoring in vivo cholinergic activity

The ability to directly measure acetylcholine (ACh) release is an essential step toward understanding its physiological function. Here we optimized the GRAB_ACh (GPCR-activation-based ACh) sensor to achieve substantially improved sensitivity in ACh detection, as well as reduced downstream coupling to intracellular pathways. The improved version of the ACh sensor retains the subsecond response kinetics, physiologically relevant affinity and precise molecular specificity for ACh of its predecessor. Using this sensor, we revealed compartmental ACh signals in the olfactory center of transgenic flies in response to external stimuli including odor and body shock. Using fiber photometry recording and two-photon imaging, our ACh sensor also enabled sensitive detection of single-trial ACh dynamics in multiple brain regions in mice performing a variety of behaviors.

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.

How Do Molecular Dynamics Data Complement Static Structural Data of GPCRs

G protein-coupled receptors (GPCRs) are implicated in nearly every physiological process in the human body and therefore represent an important drug targeting class. Advances in X-ray crystallography and cryo-electron microscopy (cryo-EM) have provided multiple static structures of GPCRs in complex with various signaling partners. However, GPCR functionality is largely determined by their flexibility and ability to transition between distinct structural conformations. Due to this dynamic nature, a static snapshot does not fully explain the complexity of GPCR signal transduction. Molecular dynamics (MD) simulations offer the opportunity to simulate the structural motions of biological processes at atomic resolution. Thus, this technique can incorporate the missing information on protein flexibility into experimentally solved structures. Here, we review the contribution of MD simulations to complement static structural data and to improve our understanding of GPCR physiology and pharmacology, as well as the challenges that still need to be overcome to reach the full potential of this technique.

Regiospecific Introduction of Halogens on the 2-Aminobiphenyl Subunit Leading to Highly Potent and Selective M3 Muscarinic Acetylcholine Receptor Antagonists and Weak Inverse Agonists

Muscarinic M₃ receptor antagonists and inverse agonists displaying high affinity and subtype selectivity over the antitarget M₂ are valuable pharmacological tools and may enable improved treatment of chronic obstructive pulmonary disease (COPD), asthma, or urinary incontinence. On the basis of known M₃ antagonists comprising a piperidine or quinuclidine unit attached to a biphenyl carbamate, 5-fluoro substitution was responsible for M₃ subtype selectivity over M₂, while 3'-chloro substitution substantially increased affinity through a σ-hole interaction. Resultantly, two piperidinyl- and two quinuclidinium-substituted biphenyl carbamates OFH243 (13n), OFH244 (13m), OFH3911 (14n), and OFH3912 (14m) were discovered, which display two-digit picomolar affinities with K_i values from 0.069 to 0.084 nM, as well as high selectivity over the M₂ subtype (46- to 68-fold). While weak inverse agonistic properties were determined for the biphenyl carbamates 13m and 13n, neutral antagonism was observed for 14m and 14n and tiotropium under identical assay conditions.

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.

G Protein-Coupled Receptors in Asthma Therapy: Pharmacology and Drug Action

Asthma is a heterogeneous inflammatory disease of the airways that is associated with airway hyperresponsiveness and airflow limitation. Although asthma was once simply categorized as atopic or nonatopic, emerging analyses over the last few decades have revealed a variety of asthma endotypes that are attributed to numerous pathophysiological mechanisms. The classification of asthma by endotype is primarily routed in different profiles of airway inflammation that contribute to bronchoconstriction. Many asthma therapeutics target G protein-coupled receptors (GPCRs), which either enhance bronchodilation or prevent bronchoconstriction. Short-acting and long-acting β 2-agonists are widely used bronchodilators that signal through the activation of the β 2-adrenergic receptor. Short-acting and long-acting antagonists of muscarinic acetylcholine receptors are used to reduce bronchoconstriction by blocking the action of acetylcholine. Leukotriene antagonists that block the signaling of cysteinyl leukotriene receptor 1 are used as an add-on therapy to reduce bronchoconstriction and inflammation induced by cysteinyl leukotrienes. A number of GPCR-targeting asthma drug candidates are also in different stages of development. Among them, antagonists of prostaglandin D2 receptor 2 have advanced into phase III clinical trials. Others, including antagonists of the adenosine A2B receptor and the histamine H4 receptor, are in early stages of clinical investigation. In the past decade, significant research advancements in pharmacology, cell biology, structural biology, and molecular physiology have greatly deepened our understanding of the therapeutic roles of GPCRs in asthma and drug action on these GPCRs. This review summarizes our current understanding of GPCR signaling and pharmacology in the context of asthma treatment. SIGNIFICANCE STATEMENT: Although current treatment methods for asthma are effective for a majority of asthma patients, there are still a large number of patients with poorly controlled asthma who may experience asthma exacerbations. This review summarizes current asthma treatment methods and our understanding of signaling and pharmacology of G protein-coupled receptors (GPCRs) in asthma therapy, and discusses controversies regarding the use of GPCR drugs and new opportunities in developing GPCR-targeting therapeutics for the treatment of asthma.

Kinase and GPCR polypharmacological approach for the identification of efficient anticancer medicines

Discovery of an anticancer medicine using a single target protein has often been unsuccessful due to the complexity of pathogenic mechanisms as well as the presence of redundant signaling pathways. In this work, we attempted to find promising anticancer drug candidates by simultaneously targeting casein kinase 1 delta (CK1δ) and muscarinic acetylcholine receptor M3 (M3R). Through the structure-based virtual screening and de novo design with the modified potential function for protein-ligand binding, a series of benzo[4,5]imidazo[1,2-a][1,3,5]triazine-2-amine (BITA) derivatives were identified as CK1δ inhibitors and also as M3R antagonists. The biochemical potencies of these bifunctional molecules reached the nanomolar and low-micromolar levels with respect to CK1δ and M3R, respectively. A common interaction feature in the calculated CK1δ-inhibitor and M3R-antagonist complexes is that the BITA moiety is well-stabilized in the orthosteric site of M3R and the hinge region of CK1δ through the establishment of the three hydrogen bonds and the hydrophobic contacts in the vicinity. The computational and experimental results found in this work exemplify the efficiency of kinase and GPCR polypharmacology in developing anticancer medicines.

[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.

Crystal structure of the M5 muscarinic acetylcholine receptor

The 5 subtypes of the muscarinic acetylcholine receptors (mAChRs) are expressed throughout the central and peripheral nervous system where they play a vital role in physiology and pathologies. Recently, the M₅ mAChR subtype has emerged as an exciting drug target for the treatment of drug addiction. We have determined the atomic structure of the M₅ mAChR bound to the clinically used inverse agonist tiotropium. The M₅ mAChR structure now allows for a full comparison of all 5 mAChR subtypes and reveals that small differences in the extracellular loop regions can mediate orthosteric and allosteric ligand selectivity. Together, these findings open the door for future structure-based design of selective drugs that target this therapeutically important class of receptors.

The human M₅ muscarinic acetylcholine receptor (mAChR) has recently emerged as an exciting therapeutic target for treating a range of disorders, including drug addiction. However, a lack of structural information for this receptor subtype has limited further drug development and validation. Here we report a high-resolution crystal structure of the human M₅ mAChR bound to the clinically used inverse agonist, tiotropium. This structure allowed for a comparison across all 5 mAChR family members that revealed important differences in both orthosteric and allosteric sites that could inform the rational design of selective ligands. These structural studies, together with chimeric swaps between the extracellular regions of the M₂ and M₅ mAChRs, provided structural insight into kinetic selectivity, where ligands show differential residency times between related family members. Collectively, our study provides important insights into the nature of orthosteric and allosteric ligand interaction across the mAChR family that could be exploited for the design of selective drugs.

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

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

Second M3 muscarinic receptor binding site contributes to bronchoprotection by tiotropium

BACKGROUND AND PURPOSE

The bronchodilator tiotropium binds not only to its main binding site on the M₃ muscarinic receptor but also to an allosteric site. Here, we have investigated the functional relevance of this allosteric binding and the potential contribution of this behaviour to interactions with long-acting β-adrenoceptor agonists, as combination therapy with anticholinergic agents and β-adrenoceptor agonists improves lung function in chronic obstructive pulmonary disease.

EXPERIMENTAL APPROACH

ACh, tiotropium, and atropine binding to M₃ receptors were modelled using molecular dynamics simulations. Contractions of bovine and human tracheal smooth muscle strips were studied.

KEY RESULTS

Molecular dynamics simulation revealed extracellular vestibule binding of tiotropium, and not atropine, to M₃ receptors as a secondary low affinity binding site, preventing ACh entry into the orthosteric binding pocket. This resulted in a low (allosteric binding) and high (orthosteric binding) functional affinity of tiotropium in protecting against methacholine-induced contractions of airway smooth muscle, which was not observed for atropine and glycopyrrolate. Moreover, antagonism by tiotropium was insurmountable in nature. This behaviour facilitated functional interactions of tiotropium with the β-agonist olodaterol, which synergistically enhanced bronchoprotective effects of tiotropium. This was not seen for glycopyrrolate and olodaterol or indacaterol but was mimicked by the interaction of tiotropium and forskolin, indicating no direct β-adrenoceptor-M₃ receptor crosstalk in this effect.

CONCLUSIONS AND IMPLICATIONS

We propose that tiotropium has two binding sites at the M₃ receptor that prevent ACh action, which, together with slow dissociation kinetics, may contribute to insurmountable antagonism and enhanced functional interactions with β-adrenoceptor agonists.

Metadynamics simulations of ligand binding to GPCRs

Recent developments in metadynamics simulation techniques for ligand binding to Class A GPCRs are described and the results obtained elucidated. The computational protocol makes good use of modern massively parallel hardware, making simulations of the binding/unbinding process routine. The simulations reveal unprecedented details of the ligand-binding pathways, including multiple binding sites in many cases. Free energies of binding are reproduced very well and the simulations allow prediction of the efficacy (agonist, antagonist etc.) of ligands.

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.

Enhancing Drug Residence Time by Shielding of Intra-Protein Hydrogen Bonds: A Case Study on CCR2 Antagonists

The target residence time (RT) for a given ligand is one of the important parameters that have to be optimized during drug design. It is well established that shielding the receptor-ligand hydrogen bond (H-bond) interactions from water has been one of the factors in increasing ligand RT. Building on this foundation, here we report that shielding an intra-protein H-bond, which confers rigidity to the binding pocket and which is not directly involved in drug-receptor interactions, can strongly influence RT for CCR2 antagonists. Based on our recently solved CCR2 structure with MK-0812 and molecular dynamics (MD) simulations, we show that the RT for this and structurally related ligands is directly dependent on the shielding of the Tyr120-Glu291 H-bond from the water. If solvated this H-bond is often broken, making the binding pocket flexible and leading to shorter RT.

Peeking at G-protein-coupled receptors through the molecular dynamics keyhole

Molecular dynamics is a state of the art computational tool for the investigation of biophysics phenomenon at a molecular scale, as it enables the modeling of dynamic processes, such as conformational motions, molecular solvation and ligand binding. The recent advances in structural biology have led to a bloom in published G-protein-coupled receptor structures, representing a solid and valuable resource for molecular dynamics studies. During the last decade, indeed, a plethora of physiological and pharmacological facets of this membrane protein superfamily have been addressed by means of molecular dynamics simulations, including the activation mechanism, allosterism and, very recently, biased signaling. Here, we try to recapitulate some of the main contributions that molecular dynamics has recently produced in the field.

M3-subtype muscarinic receptor activation stimulates intracellular calcium oscillations and aldosterone production in human adrenocortical HAC15 cells

A previous body of work in bovine and rodent models shows that cholinergic agonists modulate the secretion of steroid hormones from the adrenal cortex. In this study we used live-cell Ca²⁺ imaging to investigate cholinergic activity in the HAC15 human adrenocortical carcinoma cell line. The cholinergic agonists carbachol and acetylcholine triggered heterogeneous Ca²⁺ oscillations that were strongly inhibited by antagonists with high affinity for the M₃ muscarinic receptor subtype, while preferential block of M₁ or M₂ receptors was less effective. Acute exposure to carbachol and acetylcholine modestly elevated aldosterone secretion in HAC15 cells, and this effect was also diminished by M₃ inhibition. HAC15 cells expressed relatively high levels of mRNA for M₃ and M₂ receptors, while M₁ and M₅ mRNA were much lower. In conclusion, our data extend previous findings in non-human systems to implicate the M₃ receptor as the dominant muscarinic receptor in the human adrenal cortex.

Aminergic GPCR-Ligand Interactions: A Chemical and Structural Map of Receptor Mutation Data

The aminergic family of G protein-coupled receptors (GPCRs) plays an important role in various diseases and represents a major drug discovery target class. Structure determination of all major aminergic subfamilies has enabled structure-based ligand design for these receptors. Site-directed mutagenesis data provides an invaluable complementary source of information for elucidating the structural determinants of binding of different ligand chemotypes. The current study provides a comparative analysis of 6692 mutation data points on 34 aminergic GPCR subtypes, covering the chemical space of 540 unique ligands from mutagenesis experiments and information from experimentally determined structures of 52 distinct aminergic receptor-ligand complexes. The integrated analysis enables detailed investigation of structural receptor-ligand interactions and assessment of the transferability of combined binding mode and mutation data across ligand chemotypes and receptor subtypes. An overview is provided of the possibilities and limitations of using mutation data to guide the design of novel aminergic receptor ligands.

Evolution of the Muscarinic Acetylcholine Receptors in Vertebrates

The family of muscarinic acetylcholine receptors (mAChRs) consists of five members in mammals, encoded by the CHRM1-5 genes. The mAChRs are G-protein-coupled receptors, which can be divided into the following two subfamilies: M2 and M4 receptors coupling to G_i/o; and M1, M3, and M5 receptors coupling to G_q/11. However, despite the fundamental roles played by these receptors, their evolution in vertebrates has not yet been fully described. We have combined sequence-based phylogenetic analyses with comparisons of exon–intron organization and conserved synteny in order to deduce the evolution of the mAChR receptors. Our analyses verify the existence of two ancestral genes prior to the two vertebrate tetraploidizations (1R and 2R). After these events, one gene had duplicated, resulting in CHRM2 and CHRM4; and the other had triplicated, forming the CHRM1, CHRM3, and CHRM5 subfamily. All five genes are still present in all vertebrate groups investigated except the CHRM1 gene, which has not been identified in some of the teleosts or in chicken or any other birds. Interestingly, the third tetraploidization (3R) that took place in the teleost predecessor resulted in duplicates of all five mAChR genes of which all 10 are present in zebrafish. One of the copies of the CHRM2 and CHRM3 genes and both CHRM4 copies have gained introns in teleosts. Not a single separate (nontetraploidization) duplicate has been identified in any vertebrate species. These results clarify the evolution of the vertebrate mAChR family and reveal a doubled repertoire in zebrafish, inviting studies of gene neofunctionalization and subfunctionalization.

Utility of an "Allosteric Site-Impaired" M2 Muscarinic Acetylcholine Receptor as a Novel Construct for Validating Mechanisms of Action of Synthetic and Putative Endogenous Allosteric Modulators

Muscarinic acetylcholine receptors (mAChRs) are exemplar models for understanding G protein-coupled receptor (GPCR) allostery, possessing a "common" allosteric site in an extracellular vestibule (ECV) for synthetic modulators including gallamine, strychnine, and brucine. In addition, there is intriguing evidence of endogenous peptides/proteins that may target this region at the M₂ mAChR. A common feature of synthetic and endogenous M₂ mAChR negative allosteric modulators (NAMs) is their cationic nature. Using a structure-based approach, we previously designed a mutant M₂ mAChR (N410K+T423K) to specifically abrogate binding of ECV cationic modulators (Dror et al., 2013). Herein, we used this "allosteric site-impaired" receptor to investigate allosteric interactions of synthetic modulators as well as basic peptides (poly-l-arginine, endogenously produced protamine, and major basic protein). Using [³H]N-methylscopolamine equilibrium and kinetic binding and functional assays of guanosine 5'-O-[γ-thio]triphosphate [³⁵S] binding and extracellular signal-regulated kinases 1 and 2 phosphorylation, we found modest effects of the mutations on potencies of orthosteric antagonists and an increase in the affinity of the cognate agonist, acetylcholine, likely reflecting the effect of the mutations on the access/egress of these ligands into the orthosteric pocket. More importantly, we noted a significant abrogation in affinity for all synthetic or peptidic modulators at the mutant mAChR, validating their allosteric nature. Collectively, these findings provide evidence for a hitherto-unappreciated role of endogenous cationic peptides interacting allosterically at the M₂ mAChR and identify the allosteric site-impaired GPCR as a tool for validating NAM activity as well as a potential candidate for future chemogenetic strategies to understand the physiology of endogenous allosteric substances.

Structure-kinetic relationships that control the residence time of drug-target complexes: insights from molecular structure and dynamics

Time-dependent target occupancy is a function of both the thermodynamics and kinetics of drug-target interactions. However, while the optimization of thermodynamic affinity through approaches such as structure-based drug design is now relatively straight forward, less is understood about the molecular interactions that control the kinetics of drug complex formation and breakdown since this depends on both the ground and transition state energies on the binding reaction coordinate. In this opinion we highlight several recent examples that shed light on current approaches that are elucidating the factors that control the life-time of the drug-target complex.

Pharmacological properties of revefenacin (TD-4208), a novel, nebulized long-acting, and lung selective muscarinic antagonist, at human recombinant muscarinic receptors and in rat, guinea pig, and human isolated airway tissues

Revefenacin (TD‐4208) is a novel, long‐acting, and lung‐selective muscarinic cholinergic receptor (mAChR) antagonist in development as a nebulized inhalation solution for the treatment of chronic obstructive pulmonary disease (COPD) patients. This study evaluated the pharmacology of revefenacin at human recombinant mAChRs and in airway tissues from rats, guinea pigs, and humans. At human recombinant mAChRs, revefenacin displayed high affinity (pK_I = 8.2‐9.8) and behaved as a competitive antagonist (pKI, apparent = 9.4‐10.9) at the five human recombinant mAChRs. Kinetic studies demonstrated that revefenacin dissociated significantly slower from the hM₃ (t_1/2 = 82 minutes) compared to the hM₂ (t_1/2 = 6.9 minutes) mAChR at 37°C, thereby making it kinetically selective for the former subtype. Similarly, in functional studies, revefenacin‐mediated antagonism of acetylcholine (ACh)‐evoked calcium mobilization responses were reversed less rapidly at hM₃ compared to the hM₂ mAChR. In isolated tracheal tissues from rat and guinea pig and isolated bronchial tissues from humans, revefenacin potently antagonized mAChR‐mediated contractile responses. Furthermore, the antagonistic effects of revefenacin in rat, guinea pig, and human airway tissues were slowly reversible (t_1/2 of 13.3, >16, and >10 hours, respectively). These data demonstrate that revefenacin is a potent, high affinity, and selective functional mAChR antagonist with kinetic selectivity for the hM₃ receptor and produces potent and long‐lasting antagonism of mAChR‐mediated contractile responses in rat, guinea pig, and human airway tissue. These data suggest that revefenacin has the potential to be a potent once‐daily dosed inhaled bronchodilator in COPD patients.

Drug-Target Kinetics in Drug Discovery

The development of therapies for the treatment of neurological cancer faces a number of major challenges including the synthesis of small molecule agents that can penetrate the blood-brain barrier (BBB). Given the likelihood that in many cases drug exposure will be lower in the CNS than in systemic circulation, it follows that strategies should be employed that can sustain target engagement at low drug concentration. Time dependent target occupancy is a function of both the drug and target concentration as well as the thermodynamic and kinetic parameters that describe the binding reaction coordinate, and sustained target occupancy can be achieved through structural modifications that increase target (re)binding and/or that decrease the rate of drug dissociation. The discovery and deployment of compounds with optimized kinetic effects requires information on the structure-kinetic relationships that modulate the kinetics of binding, and the molecular factors that control the translation of drug-target kinetics to time-dependent drug activity in the disease state. This Review first introduces the potential benefits of drug-target kinetics, such as the ability to delineate both thermodynamic and kinetic selectivity, and then describes factors, such as target vulnerability, that impact the utility of kinetic selectivity. The Review concludes with a description of a mechanistic PK/PD model that integrates drug-target kinetics into predictions of drug activity.

Molecular Basis of Ligand Dissociation from G Protein-Coupled Receptors and Predicting Residence Time

G protein-coupled receptors (GPCRs) are integral membrane proteins and represent the largest class of drug targets. During the past decades progress in structural biology has enabled the crystallographic elucidation of the architecture of these important macromolecules. It also provided atomic-level visualization of ligand-receptor interactions, dramatically boosting the impact of structure-based approaches in drug discovery. However, knowledge obtained through crystallography is limited to static structural information. Less information is available showing how a ligand associates with or dissociates from a given receptor, whose importance is in fact increasingly recognized by the drug research community. Owing to recent advances in computer power and algorithms, molecular dynamics stimulations have become feasible that help in analyzing the kinetics of the ligand binding process. Here, we review what is currently known about the dynamics of GPCRs in the context of ligand association and dissociation, as determined through both crystallography and computer simulations. We particularly focus on the molecular basis of ligand dissociation from GPCRs and provide case studies that predict ligand dissociation pathways and residence time.

Recent Advances in Structure-Based Drug Design Targeting Class A G Protein-Coupled Receptors Utilizing Crystal Structures and Computational Simulations

G protein-coupled receptors (GPCRs) represent the largest and most physiologically important integral membrane protein family, and these receptors respond to a wide variety of physiological and environmental stimuli. GPCRs are among the most critical therapeutic targets for numerous human diseases, and approximately one-third of the currently marketed drugs target this receptor family. The recent breakthroughs in GPCR structural biology have significantly contributed to our understanding of GPCR function, ligand binding, and pharmacological action as well as to the design of new drugs. This perspective highlights the latest advances in GPCR structures with a focus on the receptor-ligand interactions of each receptor family in class A nonrhodopsin GPCRs as well as the structural features for their activation, biased signaling, and allosteric mechanisms. The current state-of-the-art methodologies of structure-based drug design (SBDD) approaches in the GPCR research field are also discussed.

Binding Kinetics and Pathways of Ligands to GPCRs

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

Impact, determination and prediction of drug-receptor residence times for GPCRs

The residence time of a ligand on a GPCR of interest has become an optimization parameter in many examples in drug design. Long residence times can counterbalance unfavorable pharmacokinetic parameters, contributing to compound safety, and short residence times can be a tool to avoid target related side effects. Unlike the prediction and interpretation of the structure-activity relationship (SAR) of a ligand class on a receptor, the understanding and prediction of the structure-kinetics relationship (SKR) is much more demanding. Experimental and computational approaches are described, which serve to either rationalize SKR or to predict the kinetic parameters such as on-rates and off-rates.

G protein-coupled receptors as promising cancer targets

G protein-coupled receptors (GPCRs) regulate an array of fundamental biological processes, such as growth, metabolism and homeostasis. Specifically, GPCRs are involved in cancer initiation and progression. However, compared with the involvement of the epidermal growth factor receptor in cancer, that of GPCRs have been largely ignored. Recent findings have implicated many GPCRs in tumorigenesis, tumor progression, invasion and metastasis. Moreover, GPCRs contribute to the establishment and maintenance of a microenvironment which is permissive for tumor formation and growth, including effects upon surrounding blood vessels, signaling molecules and the extracellular matrix. Thus, GPCRs are considered to be among the most useful drug targets against many solid cancers. Development of selective ligands targeting GPCRs may provide novel and effective treatment strategies against cancer and some anticancer compounds are now in clinical trials. Here, we focus on tumor related GPCRs, such as G protein-coupled receptor 30, the lysophosphatidic acid receptor, angiotensin receptors 1 and 2, the sphingosine 1-phosphate receptors and gastrin releasing peptide receptor. We also summarize their tissue distributions, activation and roles in tumorigenesis and discuss the potential use of GPCR agonists and antagonists in cancer therapy.

Crystal structures of the M1 and M4 muscarinic acetylcholine receptors

Muscarinic M1-M5 acetylcholine receptors are G-protein-coupled receptors that regulate many vital functions of the central and peripheral nervous systems. In particular, the M1 and M4 receptor subtypes have emerged as attractive drug targets for treatments of neurological disorders, such as Alzheimer's disease and schizophrenia, but the high conservation of the acetylcholine-binding pocket has spurred current research into targeting allosteric sites on these receptors. Here we report the crystal structures of the M1 and M4 muscarinic receptors bound to the inverse agonist, tiotropium. Comparison of these structures with each other, as well as with the previously reported M2 and M3 receptor structures, reveals differences in the orthosteric and allosteric binding sites that contribute to a role in drug selectivity at this important receptor family. We also report identification of a cluster of residues that form a network linking the orthosteric and allosteric sites of the M4 receptor, which provides new insight into how allosteric modulation may be transmitted between the two spatially distinct domains.