Binding kinetics of ligands acting at GPCRs

Kinetic analysis of fluorescent ligand binding to cell surface receptors: Insights into conformational changes and allosterism in living cells

Equilibrium binding assays are one of the mainstays of current drug discovery efforts to evaluate the interaction of drugs with receptors in membranes and intact cells. However, in recent years, there has been increased focus on the kinetics of the drug-receptor interaction to gain insight into the lifetime of drug-receptor complexes and the rate of association of a ligand with its receptor. Furthermore, drugs that act on topically distinct sites (allosteric) from those occupied by the endogenous ligand (orthosteric site) can induce conformational changes in the orthosteric binding site leading to changes in the association and/or dissociation rate constants of orthosteric ligands. Conformational changes in the orthosteric ligand binding site can also be induced through interaction with neighbouring accessory proteins and receptor homodimerisation and heterodimerisation. In this review, we provide an overview of the use of fluorescent ligand technologies to interrogate ligand-receptor kinetics in living cells and the novel insights that they can provide into the conformational changes induced by drugs acting on a variety of cell surface receptors including G protein-coupled receptors (GPCRs), receptor tyrosine kinases (RTKs) and cytokine receptors.

Structural basis of psychedelic LSD recognition at dopamine D1 receptor

Understanding the kinetics of LSD in receptors and subsequent induced signaling is crucial for comprehending both the psychoactive and therapeutic effects of LSD. Despite extensive research on LSD's interactions with serotonin 2A and 2B receptors, its behavior on other targets, including dopamine receptors, has remained elusive. Here, we present cryo-EM structures of LSD/PF6142-bound dopamine D1 receptor (DRD1)-legobody complexes, accompanied by a β-arrestin-mimicking nanobody, NBA3, shedding light on the determinants of G protein coupling versus β-arrestin coupling. Structural analysis unveils a distinctive binding mode of LSD in DRD1, particularly with the ergoline moiety oriented toward TM4. Kinetic investigations uncover an exceptionally rapid dissociation rate of LSD in DRD1, attributed to the flexibility of extracellular loop 2 (ECL2). Moreover, G protein can stabilize ECL2 conformation, leading to a significant slowdown in ligand's dissociation rate. These findings establish a solid foundation for further exploration of G protein-coupled receptor (GPCR) dynamics and their relevance to signal transduction.

IS4-FAM, a fluorescent tool to study CXCR4 affinity and competitive antagonism in native cancer cells

The ability to accurately measure drug-target interaction is critical for the discovery of new therapeutics. Classical pharmacological bioassays such as radioligand or fluorescent ligand binding assays can define the affinity or Kd of a ligand for a receptor with the lower the Kd, the stronger the binding and the higher the affinity. However, in many drug discovery laboratories today, the target of interest if often artificially upregulated by means of transfection to modify the host cell's genetic makeup. This then potentially invalidates the assumptions of classical pharmacology affinity calculations as the receptor of interest is no longer at normal physiological densities. The CXCR4 receptor is expressed on many different cancer cell types and is associated with metastasis and poor prognosis. Therefore, the CXCR4 receptor is a desirable target for novel therapeutics. In this study, we explore the applicability of the newly developed fluorescently tagged CXCR4 antagonists, IS4-FAM as an investigative tool to study CXCR4 affinity and competitive antagonism in native, non-transfected cancer cells using confocal microscopy and flow cytometry. IS4-FAM directly labels CXCR4 in several cell lines including high CXCR4 expressing SK-MEL-28 (malignant melanoma) and PC3 (metastatic prostate cancer) and lower CXCR4 expressing THP-1 (acute monocytic leukemia) and was competitive with the established CXCR4 antagonist, AMD3100. This highlights the potential of IS4-FAM as a pharmacological tool for drug discovery in native cells lines and tissues.

Collagen-targeted protein nanomicelles for the imaging of non-alcoholic steatohepatitis

In vivo molecular imaging tools hold immense potential to drive transformative breakthroughs by enabling researchers to visualize cellular and molecular interactions in real-time and/or at high resolution. These advancements will facilitate a deeper understanding of fundamental biological processes and their dysregulation in disease states. Here, we develop and characterize a self-assembling protein nanomicelle called collagen type I binding - thermoresponsive assembled protein (Col1-TRAP) that binds tightly to type I collagen in vitro with nanomolar affinity. For ex vivo visualization, Col1-TRAP is labeled with a near-infrared fluorescent dye (NIR-Col1-TRAP). Both Col1-TRAP and NIR-Col1-TRAP display approximately a 3.8-fold greater binding to type I collagen compared to TRAP when measured by surface plasmon resonance (SPR). We present a proof-of-concept study using NIR-Col1-TRAP to detect fibrotic type I collagen deposition ex vivo in the livers of mice with non-alcoholic steatohepatitis (NASH). We show that NIR-Col1-TRAP demonstrates significantly decreased plasma recirculation time as well as increased liver accumulation in the NASH mice compared to mice without disease over 4 hours. As a result, NIR-Col1-TRAP shows potential as an imaging probe for NASH with in vivo targeting performance after injection in mice. STATEMENT OF SIGNIFICANCE: Direct molecular imaging of fibrosis in NASH patients enables the diagnosis and monitoring of disease progression with greater specificity and resolution than do elastography-based methods or blood tests. In addition, protein-based imaging probes are more advantageous than alternatives due to their biodegradability and scalable biosynthesis. With the aid of computational modeling, we have designed a self-assembled protein micelle that binds to fibrillar and monomeric collagen in vitro. After the protein was labeled with near-infrared fluorescent dye, we injected the compound into mice fed on a NASH diet. NIR-Col1-TRAP clears from the serum faster in these mice compared to control mice, and accumulates significantly more in fibrotic livers.This work advances the development of targeted protein probes for in vivo fibrosis imaging.

Application of FRET- and BRET-based live-cell biosensors in deorphanization and ligand discovery studies on orphan G protein-coupled receptors

Bioluminescence- and fluorescence-based resonance energy transfer assays have gained considerable attention in pharmacological research as high-throughput scalable tools applicable to drug discovery. To this end, G protein-coupled receptors represent the biggest target class for marketed drugs, and among them, orphan G protein-coupled receptors have the biggest untapped therapeutic potential. In this review, the cases where biophysical methods, BRET and FRET, were employed for deorphanization and ligand discovery studies on orphan G protein-coupled receptors are listed and discussed.

Highly Selective Drug-Derived Fluorescent Probes for the Cannabinoid Receptor Type 1 (CB1R)

The cannabinoid receptor type 1 (CB1R) is pivotal within the endocannabinoid system regulating various signaling cascades with effects in appetite regulation, pain perception, memory formation, and thermoregulation. Still, understanding of CB1R's cellular signaling, distribution, and expression dynamics is very fragmentary. Real-time visualization of CB1R is crucial for addressing these questions. Selective drug-like CB1R ligands with a defined pharmacological profile were investigated for the construction of CB1R fluorescent probes using a reverse design-approach. A modular design concept with a diethyl glycine-based building block as the centerpiece allowed for the straightforward synthesis of novel probe candidates. Validated by computational docking studies, radioligand binding, and cAMP assay, this systematic approach allowed for the identification of novel pyrrole-based CB1R fluorescent probes. Application in fluorescence-based target-engagement studies and live cell imaging exemplify the great versatility of the tailored CB1R probes for investigating CB1R localization, trafficking, pharmacology, and its pathological implications.

Time-dependent ligand-receptor binding kinetics and functionality in a heterodimeric receptor model

GPCRs heteromerize both in CNS and non-CNS regions. The cell uses receptor heteromerization to modulate receptor functionality and to provide fine tuning of receptor signaling. In order for pharmacologists to explore these mechanisms for therapeutic purposes, quantitative receptor models are needed. We have developed a time-dependent model of the binding kinetics and functionality of a preformed heterodimeric receptor involving two drugs. Two cases were considered: both or only one of the drugs are in excess with respect to the total concentration of the receptor. The latter case can be applied to those situations in which a drug causes unwanted side effects that need to be reduced by decreasing its concentration. The required efficacy can be maintained by the allosteric effects mutually exerted by the two drugs in the two-drug combination system. We discuss this concept assuming that the drug causing unwanted side effects is an opioid and that analgesia is the therapeutic effect. As additional points, allosteric modulation by endogenous compounds and synthetic bivalent ligands was included in the study. Receptor heteromerization offers a mechanistic understanding and quantification of the pharmacological effects elicited by combinations of two drugs at different doses and with different efficacies and cooperativity effects, thus providing a conceptual framework for drug combination therapy.

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.

Drugging the entire human proteome: Are we there yet?

Each of the ∼20,000 proteins in the human proteome is a potential target for compounds that bind to it and modify its function. The 3D structures of most of these proteins are now available. Here, we discuss the prospects for using these structures to perform proteome-wide virtual HTS (VHTS). We compare physics-based (docking) and AI VHTS approaches, some of which are now being applied with large databases of compounds to thousands of targets. Although preliminary proteome-wide screens are now within our grasp, further methodological developments are expected to improve the accuracy of the results.

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

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

Ligand selectivity hotspots in serotonin GPCRs

Serotonin is a neurotransmitter regulating numerous physiological processes also modulated by drugs, for example, schizophrenia, depression, migraine, and obesity. However, these drugs typically have adverse effects caused by promiscuous binding across 12 serotonin and more than 20 homologous receptors. Recently, structures of the entire serotonin receptor family uncovered molecular ligand recognition. Here, we present a map of 19 'selectivity hotspots', that is, nonconserved binding site residues governing selectivity via favorable target interactions or repulsive 'off-target' contacts. Furthermore, we review functional rationale from observed ligand-binding affinities and mutagenesis effects. Unifying knowledge underlying specific probes and drugs is critical toward the functional characterization of different receptors and alleviation of adverse effects.

Subtle Structural Modification of a Synthetic Cannabinoid Receptor Agonist Drastically Increases its Efficacy at the CB1 Receptor

The emergence of synthetic cannabinoid receptor agonists (SCRAs) as illicit psychoactive substances has posed considerable public health risks, including fatalities. Many SCRAs exhibit much higher efficacy and potency compared with the phytocannabinoid Δ9-tetrahydrocannabinol (THC) at the cannabinoid receptor 1 (CB1R), leading to dramatic differences in signaling levels that can be toxic. In this study, we investigated the structure-activity relationships of aminoalkylindole SCRAs at CB1Rs, focusing on 5F-pentylindoles containing an amide linker attached to different head moieties. Using in vitro bioluminescence resonance energy transfer assays, we identified a few SCRAs exhibiting significantly higher efficacy in engaging the Gi protein and recruiting β-arrestin than the reference CB1R full agonist CP55940. Importantly, the extra methyl group on the head moiety of 5F-MDMB-PICA, as compared to that of 5F-MMB-PICA, led to a large increase in efficacy and potency at the CB1R. This pharmacological observation was supported by the functional effects of these SCRAs on glutamate field potentials recorded in hippocampal slices. Molecular modeling and simulations of the CB1R models bound with both of the SCRAs revealed critical structural determinants contributing to the higher efficacy of 5F-MDMB-PICA and how these subtle differences propagated to the receptor-G protein interface. Thus, we find that apparently minor structural changes in the head moiety of SCRAs can cause major changes in efficacy. Our results highlight the need for close monitoring of the structural modifications of newly emerging SCRAs and their potential for toxic drug responses in humans.

Encounter Complexes Between the N-terminal of Neurotensin with the Extracellular Loop 2 of the Neurotensin Receptor 1 Steer Neurotensin to the Orthosteric Binding Pocket

Neurotensin (NT) is a linear disordered peptide that activates two different class A GPCRs, neurotensin receptor 1 (NTS1) and NTS2. Resolved structures of the complex of the C-terminal fragment of NT, NT8-13, with NTS1 shows the peptide takes a well-defined structure in the bound state. However, the mechanisms underlying NT recognition of NTS1, and the conformational transition of NT upon binding NTS1 is an open question that if answered may aid discovery of highly selective drugs and reveal potential secondary binding sites on the surface of the receptor. Herein we investigated the interactions guiding NT to the orthosteric binding pocket of NTS1 by combining NMR experiments with kinetic analysis of the binding pathway using stopped-flow fluorescence and mutagenesis on both NT and NTS1. We show the presence of transient structures in the middle part of NT that kinetically regulate the binding of NT to NTS1. Moreover, our results indicate that the binding pathway of NT onto NTS1 is mediated via electrostatic interactions between the N-terminal region of NT with the extracellular loop 2 of NTS1. These interactions induce backbone conformational changes in neurotensin similar to the bound-state neurotensin, suggesting that the N-terminal region of NT and these interactions should be considered for development of selective drugs against NTS1.

Comparative Saturation Binding Analysis of 64Cu-Labeled Somatostatin Analogues Using Cell Homogenates and Intact Cells

The development of novel ligands for G-protein-coupled receptors (GPCRs) typically entails the characterization of their binding affinity, which is often performed with radioligands in a competition or saturation binding assay format. Since GPCRs are transmembrane proteins, receptor samples for binding assays are prepared from tissue sections, cell membranes, cell homogenates, or intact cells. As part of our investigations on modulating the pharmacokinetics of radiolabeled peptides for improved theranostic targeting of neuroendocrine tumors with a high abundance of the somatostatin receptor sub-type 2 (SST₂), we characterized a series of ⁶⁴Cu-labeled [Tyr³]octreotate (TATE) derivatives in vitro in saturation binding assays. Herein, we report on the SST₂ binding parameters measured toward intact mouse pheochromocytoma cells and corresponding cell homogenates and discuss the observed differences taking the physiology of SST₂ and GPCRs in general into account. Furthermore, we point out method-specific advantages and limitations.

A subtle structural modification of a synthetic cannabinoid receptor agonist drastically increases its efficacy at the CB1 receptor

The emergence of synthetic cannabinoid receptor agonists (SCRAs) as illicit psychoactive substances has posed considerable public health risks that include fatalities. Many SCRAs exhibit much higher efficacy and potency, compared with the phytocannabinoid Δ9-tetrahydrocannabinol (THC), at the cannabinoid receptor 1 (CB1R), a G protein-coupled receptor involved in modulating neurotransmitter release. In this study, we investigated structure activity relationships (SAR) of aminoalkylindole SCRAs at CB1Rs, focusing on 5F-pentylindoles containing an amide linker attached to different head moieties. Using in vitro bioluminescence resonance energy transfer (BRET) assays, we identified a few of SCRAs exhibiting significantly higher efficacy in engaging the Gi protein and recruiting β-arrestin than the reference CB1R full agonist CP55940. Importantly, adding a methyl group at the head moiety of 5F-MMB-PICA yielded 5F-MDMB-PICA, an agonist exhibiting a large increase in efficacy and potency at the CB1R. This pharmacological observation was supported by a functional assay of the effects of these SCRAs on glutamate field potentials recorded in hippocampal slices. Molecular modeling and simulations of the CB1R bound with either of the SCRAs revealed critical structural determinants contributing to the higher efficacy of 5F-MDMB-PICA, and how these subtle differences propagated to the receptor-G protein interface. Thus, we find that apparently minor structural changes in the head moiety of SCRAs can cause major changes in efficacy. Our results highlight the need for close monitoring of structural modifications of newly emerging SCRAs and their potential for toxic drug responses in humans.

An increase in the cytosolic concentration of free calcium ions activates the neutrophil NADPH-oxidase provided that the free fatty acid receptor 2 has been allosterically modulated

Signals generated by free fatty acid receptor 2 (FFA2R) can activate the neutrophil NADPH-oxidase without involvement of any orthosteric FFA2R agonist. The initiating signals may be generated by P2Y₂R, the receptor for ATP. An FFA2R specific allosteric modulator (PAM; Cmp58) was required for this response and used to investigate the mechanism by which signals generated by ATP/P2Y₂R activate an FFA2R dependent process. The P2Y₂R induced signal that together with the modulated FFA2R activates neutrophils, was generated downstream of the Gα_q containing G protein coupled to P2Y₂R. A rise in the cytosolic concentration of ionized calcium ([Ca²⁺]_i) was hypothesized to be the important signal. The hypothesis gained support from the finding that the modulator transferred the neutrophils to a Ca²⁺sensitive state. The rise in [Ca²⁺]_i induced by the Ca²⁺ specific ionophore ionomycin, activated the neutrophils provided that an allosteric modulator was bound to FFA2R. The activity of the superoxide generating NADPH-oxidase induced by ionomycin was rapidly terminated and the FFA2Rs could then no longer be activated by the FFA2R agonist propionate or by the signal generated by ATP/P2Y₂R. The non-responding state of FFA2R was, however, revoked by a cross-activating allosteric FFA2R modulator. The [Ca²⁺]_i mediated activation of neutrophils with their FFA2Rs allosterically modulated, represent a unique regulatory receptor crosstalk mechanism by which the activation potency of a G protein coupled receptor is controlled by a receptor-crosstalk signaling system operating from the cytosolic side of the plasma membrane.

Novel Cannabinoid Receptor 2 (CB2) Low Lipophilicity Agonists Produce Distinct cAMP and Arrestin Signalling Kinetics without Bias

Cannabinoid Receptor 2 (CB2) is a promising target for treating inflammatory diseases. We designed derivatives of 3-carbamoyl-2-pyridone and 1,8-naphthyridin-2(1H)-one-3-carboxamide CB2-selective agonists with reduced lipophilicity. The new compounds were measured for their affinity (radioligand binding) and ability to elicit cyclic adenosine monophosphate (cAMP) signalling and β-arrestin-2 translocation with temporal resolution (BRET-based biosensors). For the 3-carbamoyl-2-pyridone derivatives, we found that modifying the previously reported compound UOSS77 (also known as S-777469) by appending a PEG2-alcohol via a 3-carbomylcyclohexyl carboxamide (UOSS75) lowered lipophilicity, and preserved binding affinity and signalling profile. The 1,8-naphthyridin-2(1H)-one-3-carboxamide UOMM18, containing a cis configuration at the 3-carboxamide cyclohexyl and with an alcohol on the 4-position of the cyclohexyl, had lower lipophilicity but similar CB2 affinity and biological activity to previously reported compounds of this class. Relative to CP55,940, the new compounds acted as partial agonists and did not exhibit signalling bias. Interestingly, while all compounds shared similar temporal trajectories for maximal efficacy, differing temporal trajectories for potency were observed. Consequently, when applied at sub-maximal concentrations, CP55,940 tended to elicit sustained (cAMP) or increasing (arrestin) responses, whereas responses to the new compounds tended to be transient (cAMP) or sustained (arrestin). In future studies, the compounds characterised here may be useful in elucidating the consequences of differential temporal signalling profiles on CB2-mediated physiological responses.

Biophysical investigations of class A GPCRs

G protein-coupled receptors (GPCRs) play a central role in cellular communication, converting external stimuli into intracellular responses. GPCRs bind a very broad panel of ligands, such as hormones, neurotransmitters, peptides and lipids. Ligand binding triggers a series of receptor conformational rearrangements, enabling the coupling to intracellular partners and the activation of signaling cascades. The major breakthrough in GPCRs structural biology of the past decade has considerably advanced our understanding of GPCR activation. However, structural information cannot fully explain the molecular details of GPCRs pharmacology. Biophysical investigations reveal that GPCRs are very dynamic proteins, capable of exploring a wide range of conformational states. Binding to ligands of various pharmacological classes, as well as intracellular effectors and allosteric modulators, can shift the equilibrium between these states and the kinetic of interconversions among the different conformers. Investigation of GPCR dynamic interplay is therefore important to better understand the complex pharmacology and signaling profile of these receptors.

Allosteric binding cooperativity in a kinetic context

Allosteric modulators are of prime interest in drug discovery. These drugs regulate the binding and function of endogenous ligands, with some advantages over orthosteric ligands. A typical pharmacological parameter in allosteric modulation is binding cooperativity. This property can yield unexpected but illuminating results when decomposed into its kinetic parameters. Using two reference models (the allosteric ternary complex receptor model and a heterodimer receptor model), a relationship has been derived for the cooperativity rate constant parameters. This relationship allows many combinations of the cooperativity kinetic parameters for a single binding cooperativity value obtained under equilibrium conditions. This assessment may help understand striking experimental results involving allosteric modulation and suggest further investigations in the field.

A Single Point Mutation Blocks the Entrance of Ligands to the Cannabinoid CB2 Receptor via the Lipid Bilayer

Molecular dynamic (MD) simulations have become a common tool to study the pathway of ligand entry to the orthosteric binding site of G protein-coupled receptors. Here, we have combined MD simulations and site-directed mutagenesis to study the binding process of the potent JWH-133 agonist to the cannabinoid CB₂ receptor (CB₂R). In CB₂R, the N-terminus and extracellular loop 2 fold over the ligand binding pocket, blocking access to the binding cavity from the extracellular environment. We, thus, hypothesized that the binding pathway is a multistage process consisting of the hydrophobic ligand diffusing in the lipid bilayer to contact a lipid-facing vestibule, from which the ligand enters an allosteric site inside the transmembrane bundle through a tunnel formed between TMs 1 and 7 and finally moving from the allosteric to the orthosteric binding cavity. This pathway was experimentally validated by the Ala282^7.36Phe mutation that blocks the entrance of the ligand, as JWH-133 was not able to decrease the forskolin-induced cAMP levels in cells expressing the mutant receptor. This proposed ligand entry pathway defines transient binding sites that are potential cavities for the design of synthetic modulators.

Discovery and Structure-Activity Relationship Studies of Novel Adenosine A1 Receptor-Selective Agonists

A series of benzyloxy and phenoxy derivatives of the adenosine receptor agonists N6-cyclopentyl adenosine (CPA) and N6-cyclopentyl 5'-N-ethylcarboxamidoadenosine (CP-NECA) were synthesized, and their potency and selectivity were assessed. We observed that the most potent were the compounds with a halogen in the meta position on the aromatic ring of the benzyloxy- or phenoxycyclopentyl substituent. In general, the NECA-based compounds displayed greater A1R selectivity than the adenosine-based compounds, with N6-2-(3-bromobenzyloxy)cyclopentyl-NECA and N6-2-(3-methoxyphenoxy)cyclopentyl-NECA showing ∼1500-fold improved A1R selectivity compared to NECA. In addition, we quantified the compounds' affinity and kinetics of binding at both human and rat A1R using a NanoBRET binding assay and found that the halogen substituent in the benzyloxy- or phenoxycyclopentyl moiety seems to confer high affinity for the A1R. Molecular modeling studies suggested a hydrophobic subpocket as contributing to the A1R selectivity displayed. We believe that the identified selective potent A1R agonists are valuable tool compounds for adenosine receptor research.

Development of fluorescent peptide G protein-coupled receptor activation biosensors for NanoBRET characterization of intracellular allosteric modulators

G protein-coupled receptors (GPCRs) are widely therapeutically targeted, and recent advances in allosteric modulator development at these receptors offer further potential for exploitation. Intracellular allosteric modulators (IAM) represent a class of ligands that bind to the receptor-effector interface (e.g., G protein) and inhibit agonist responses noncompetitively. This potentially offers greater selectivity between receptor subtypes compared to classical orthosteric ligands. However, while examples of IAM ligands are well described, a more general methodology for assessing compound interactions at the IAM site is lacking. Here, fluorescent labeled peptides based on the Gα peptide C terminus are developed as novel binding and activation biosensors for the GPCR-IAM site. In TR-FRET binding studies, unlabeled peptides derived from the Gαs subunit were first characterized for their ability to positively modulate agonist affinity at the β₂ -adrenoceptor. On this basis, a tetramethylrhodamine (TMR) labeled tracer was synthesized based on the 19 amino acid Gαs peptide (TMR-Gαs19cha18, where cha = cyclohexylalanine). Using NanoBRET technology to detect binding, TMR-Gαs19cha18 was recruited to Gs coupled β₂ -adrenoceptor and EP₂ receptors in an agonist-dependent manner, but not the Gi-coupled CXCR2 receptor. Moreover, NanoBRET competition binding assays using TMR-Gαs19cha18 enabled direct assessment of the affinity of unlabeled ligands for β₂ -adrenoceptor IAM site. Thus, the NanoBRET platform using fluorescent-labeled G protein peptide mimetics offers novel potential for medium-throughput screens to identify IAMs, applicable across GPCRs coupled to a G protein class. Using the same platform, Gs peptide biosensors also represent useful tools to probe orthosteric agonist efficacy and the dynamics of receptor activation.

Measurement of Protein Kinetics Using a Liquid Phase-Based Biosensing Platform

Macromolecular association is crucial to many fields in biomedical sciences, including drug development, gene editing, and diagnostics. In particular, protein-protein association and dissociation rate constants are typically determined using surface plasmon resonance systems, which require costly instrumentation and cumbersome procedures (e.g., blocking, washing, and separation). Herein, we demonstrate that protein-binding constants can be readily determined using a real-time biosensing platform facilitated by graphene oxide-modified microwell plates and fluorophore-labeled proteins, where the fluorescent probes remain highly fluorescent during protein association, whereas fluorescent bioprobes that are not associated with their counterparts are quenched by graphene oxide. Binding data of three pairs of proteins were systematically determined employing this single-step platform and compared with those data reported by the suppliers or the literature, suggesting that this approach is comparable and consistent with the existing ones. Such pairs include (i) human immunoglobulin G (H-IgG)-fluorophore-labeled anti-H-IgG, (ii) prostate-specific antigen (PSA)-quantum dot-labeled anti-PSA, and (iii) anti-RBD-fluorophore-labeled SARS-CoV-2 spike receptor-binding domain recombinant protein. We also offer an open-source software that automatically determines the binding kinetics constants of proteins. This Technical Note introduces a simple, yet effective, platform to determine relevant information on protein kinetics, which can be performed using a microwell plate reader and economical materials like graphene oxide. We foresee a new generation of diagnostics based on our affordable protein kinetics analysis.

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

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

Membrane Lipids Are an Integral Part of Transmembrane Allosteric Sites in GPCRs: A Case Study of Cannabinoid CB1 Receptor Bound to a Negative Allosteric Modulator, ORG27569, and Analogs

A growing number of G-protein-coupled receptor (GPCR) structures reveal novel transmembrane lipid-exposed allosteric sites. Ligands must first partition into the surrounding membrane and take lipid paths to these sites. Remarkably, a significant part of the bound ligands appears exposed to the membrane lipids. The experimental structures do not usually account for the surrounding lipids, and their apparent contribution to ligand access and binding is often overlooked and poorly understood. Using classical and enhanced molecular dynamics simulations, we show that membrane lipids are critical in the access and binding of ORG27569 and its analogs at the transmembrane site of cannabinoid CB1 receptor. The observed differences in the binding affinity and cooperativity arise from the functional groups that interact primarily with lipids. Our results demonstrate the significance of incorporating membrane lipids as an integral component of transmembrane sites for accurate characterization, binding-affinity calculations, and lead optimization in drug discovery.

Activation and Speciation Mechanisms in Class A GPCRs

Accurate development of allosteric modulators of GPCRs require a thorough assessment of their sequence, structure, and dynamics, toward gaining insights into their mechanisms of actions shared by family members, as well as dynamic features that distinguish subfamilies. Building on recent progress in the characterization of the signature dynamics of proteins, we analyzed here a dataset of 160 Class A GPCRs to determine their sequence similarities, structural landscape, and dynamic features across different species (human, bovine, mouse, squid, and rat), different activation states (active/inactive), and different subfamilies. The two dominant directions of variability across experimentally resolved structures, identified by principal component analysis of the dataset, shed light to cooperative mechanisms of activation, subfamily differentiation, and speciation of Class A GPCRs. The analysis reveals the functional significance of the conformational flexibilities of specific structural elements, including: the dominant role of the intracellular loop 3 (ICL3) together with the cytoplasmic ends of the adjoining helices TM5 and TM6 in enabling allosteric activation; the role of particular structural motifs at the extracellular loop 2 (ECL2) connecting TM4 and TM5 in binding ligands specific to different subfamilies; or even the differentiation of the N-terminal conformation across different species. Detailed analyses of the modes of motions accessible to the members of the dataset and their variations across members demonstrate how the active and inactive states of GPCRs obey distinct conformational dynamics. The collective fluctuations of the GPCRs are robustly defined in the active state, while the inactive conformers exhibit broad variance among members.

Aβ peptides stabilize GPCRs in inactive form and trigger inverse agonism in Alzheimer's disease

Several G-protein coupled receptors (GPCR) are upregulated in Alzheimer's Disease (AD), which ought to facilitate neurotransmission, and improve cognition. Yet, despite this upregulation, associated physiological effects are not observed in AD patients. This paradox solicits urgent attention to find a suitable justification for disturbed neurotransmission in AD. This article focuses on the role of Aβ granules and their possible interaction with GPCRs that modulate neurotransmission and AD progression.

Micellar Composition Affects Lipid Accretion Kinetics in Molecular Dynamics Simulations: Support for Lipid Network Reproduction

Mixed lipid micelles were proposed to facilitate life through their documented growth dynamics and catalytic properties. Our previous research predicted that micellar self-reproduction involves catalyzed accretion of lipid molecules by the residing lipids, leading to compositional homeostasis. Here, we employ atomistic Molecular Dynamics simulations, beginning with 54 lipid monomers, tracking an entire course of micellar accretion. This was done to examine the self-assembly of variegated lipid clusters, allowing us to measure entry and exit rates of monomeric lipids into pre-micelles with different compositions and sizes. We observe considerable rate-modifications that depend on the assembly composition and scrutinize the underlying mechanisms as well as the energy contributions. Lastly, we describe the measured potential for compositional homeostasis in our simulated mixed micelles. This affirms the basis for micellar self-reproduction, with implications for the study of the origin of life.

Biochemical characterization of cell-free synthesized human β1 adrenergic receptor cotranslationally inserted into nanodiscs

Cell-free expression enables direct cotranslational insertion of G protein coupled receptors (GPCRs) and other membrane proteins into the defined membrane environments of nanodiscs. This technique avoids GPCR contacts with detergents and allows rapid identification of lipid effects on GPCR function as well as fast screening of receptor derivatives. Critical steps of conventional GPCR preparation from cellular membranes followed by detergent-based reconstitution into nanodisc membranes are thus eliminated. We report the efficient cotranslational insertion of full-length human β₁-adrenergic receptor and of a truncated derivative into preformed nanodisc membranes. Their biochemical characterization revealed significant differences in lipid requirements, dimer formation and ligand binding activity. The truncated receptor showed a higher affinity to most tested ligands, in particular in presence of choline-containing lipids. However, introducing the naturally occurring G389R polymorphism in the full-length receptor resulted into an increased affinity to the antagonists alprenolol and carvedilol. Receptor quality was generally improved by coexpression with the agonist isoproterenol and the percentage of the ligand binding active fraction was twofold increased. Specific coupling of full-length and truncated human receptors in nanodisc membranes to Mini-Gα_s protein as well as to purified G_s heterotrimer could be demonstrated and homogeneity of purified GPCR/G_s protein complexes in nanodiscs was demonstrated by negative stain single particle analysis.

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.

Kinetic profiling and functional characterization of 8-phenylxanthine derivatives as A2B adenosine receptor antagonists

A_2B adenosine receptor (A_2BAR) antagonists have therapeutic potential in inflammation-related diseases such as asthma, chronic obstructive pulmonary disease and cancer. However, no drug is currently clinically approved, creating a demand for research on novel antagonists. Over the last decade, the study of target binding kinetics, along with affinity and potency, has been proven valuable in early drug discovery stages, as it is associated with improved in vivo drug efficacy and safety. In this study, we report the synthesis and biological evaluation of a series of xanthine derivatives as A_2BAR antagonists, including an isothiocyanate derivative designed to bind covalently to the receptor. All 28 final compounds were assessed in radioligand binding experiments, to evaluate their affinity and for those qualifying, kinetic binding parameters. Both structure-affinity and structure-kinetic relationships were derived, providing a clear relationship between affinity and dissociation rate constants. Two structurally similar compounds, 17 and 18, were further evaluated in a label-free assay due to their divergent kinetic profiles. An extended cellular response was associated with long A_2BAR residence times. This link between a ligand's A_2BAR residence time and its functional effect highlights the importance of binding kinetics as a selection parameter in the early stages of drug discovery.

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.

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

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

G Protein-Coupled Receptor-Ligand Dissociation Rates and Mechanisms from τRAMD Simulations

There is a growing appreciation of the importance of drug-target binding kinetics for lead optimization. For G protein-coupled receptors (GPCRs), which mediate signaling over a wide range of time scales, the drug dissociation rate is often a better predictor of in vivo efficacy than binding affinity, although it is more challenging to compute. Here, we assess the ability of the τ-Random Acceleration Molecular Dynamics (τRAMD) approach to reproduce relative residence times and reveal dissociation mechanisms and the effects of allosteric modulation for two important membrane-embedded drug targets: the β2-adrenergic receptor and the muscarinic acetylcholine receptor M2. The dissociation mechanisms observed in the relatively short RAMD simulations (in which molecular dynamics (MD) simulations are performed using an additional force with an adaptively assigned random orientation applied to the ligand) are in general agreement with much more computationally intensive conventional MD and metadynamics simulations. Remarkably, although decreasing the magnitude of the random force generally reduces the number of egress routes observed, the ranking of ligands by dissociation rate is hardly affected and agrees well with experiment. The simulations also reproduce changes in residence time due to allosteric modulation and reveal associated changes in ligand dissociation pathways.

Exploring Ligand Binding to Calcitonin Gene-Related Peptide Receptors

Class B1 G protein-coupled receptors (GPCRs) are important targets for many diseases, including cancer, diabetes, and heart disease. All the approved drugs for this receptor family are peptides that mimic the endogenous activating hormones. An understanding of how agonists bind and activate class B1 GPCRs is fundamental for the development of therapeutic small molecules. We combined supervised molecular dynamics (SuMD) and classic molecular dynamics (cMD) simulations to study the binding of the calcitonin gene-related peptide (CGRP) to the CGRP receptor (CGRPR). We also evaluated the association and dissociation of the antagonist telcagepant from the extracellular domain (ECD) of CGRPR and the water network perturbation upon binding. This study, which represents the first example of dynamic docking of a class B1 GPCR peptide, delivers insights on several aspects of ligand binding to CGRPR, expanding understanding of the role of the ECD and the receptor-activity modifying protein 1 (RAMP1) on agonist selectivity.

Measuring the rapid kinetics of receptor-ligand interactions in live cells using NanoBRET

The importance of receptor-ligand binding kinetics has often been overlooked during drug development, however, over the past decade it has become increasingly clear that a better understanding of the kinetic parameters is crucial for fully evaluating pharmacological effects of a drug. One technique enabling us to measure the real-time kinetics of receptor-ligand interactions in live cells is NanoBRET, which is a bioluminescence resonance energy transfer (BRET)-based assay that uses Nano luciferase. The assay described here allows the measurement of kinetic parameters of a fluorescent ligand and an unlabeled ligand binding to the same place at the receptor, as well as monitoring the effects of another compound like an allosteric modulator on the ligand binding.

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

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

Subtype selective fluorescent ligands based on ICI 118,551 to study the human β2-adrenoceptor in CRISPR/Cas9 genome-edited HEK293T cells at low expression levels

Fluorescent ligand technologies have proved to be powerful tools to improve our understanding of ligand‐receptor interactions. Here we have characterized a small focused library of nine fluorescent ligands based on the highly selective β₂‐adrenoceptor (β₂AR) antagonist ICI 118,551. The majority of fluorescent ICI 118,551 analogs had good affinity for the β₂AR (pK_D >7.0) with good selectivity over the β₁AR (pK_D <6.0). The most potent and selective ligands being 8c (ICI 118,551‐Gly‐Ala‐BODIPY‐FL‐X; β₂AR pK_D 7.48), 9c (ICI 118,551‐βAla‐βAla‐BODIPY‐FL‐X; β₂AR pK_D 7.48), 12a (ICI 118,551‐PEG‐BODIPY‐X‐630/650; β₂AR pK_D 7.56), and 12b (ICI 118,551‐PEG‐BODIPY‐FL; β₂AR pK_D 7.42). 9a (ICI 118,551‐βAla‐βAla‐BODIPY‐X‐630/650) had the highest affinity at recombinant β₂ARs (pK_D 7.57), but also exhibited significant binding affinity to the β₁AR (pK_D 6.69). Nevertheless, among the red fluorescent ligands, 9a had the best imaging characteristics in recombinant HEK293 T cells and labeling was mostly confined to the cell surface. In contrast, 12a showed the highest propensity to label intracellular β₂ARs in HEK293 T cell expressing exogenous β₂ARs. This suggests that a combination of the polyethylene glycol (PEG) linker and the BODIPY‐X‐630/650 makes this ICI 118,551 derivative particularly susceptible to crossing the cell membrane to access the intracellular β₂ARs. We have also used these ligands in combination with CRISPR/Cas9 genome‐edited HEK293 T cells to undertake for the first time real‐time ligand binding to native HEK293 T β₂ARs at low native receptor expression levels. These studies provided quantitative data on ligand‐binding characteristics but also allowed real‐time visualization of the ligand‐binding interactions in genome‐edited cells using NanoBRET luminescence imaging.

An Integrated Approach toward NanoBRET Tracers for Analysis of GPCR Ligand Engagement

Gaining insight into the pharmacology of ligand engagement with G-protein coupled receptors (GPCRs) under biologically relevant conditions is vital to both drug discovery and basic research. NanoLuc-based bioluminescence resonance energy transfer (NanoBRET) monitoring competitive binding between fluorescent tracers and unmodified test compounds has emerged as a robust and sensitive method to quantify ligand engagement with specific GPCRs genetically fused to NanoLuc luciferase or the luminogenic HiBiT peptide. However, development of fluorescent tracers is often challenging and remains the principal bottleneck for this approach. One way to alleviate the burden of developing a specific tracer for each receptor is using promiscuous tracers, which is made possible by the intrinsic specificity of BRET. Here, we devised an integrated tracer discovery workflow that couples machine learning-guided in silico screening for scaffolds displaying promiscuous binding to GPCRs with a blend of synthetic strategies to rapidly generate multiple tracer candidates. Subsequently, these candidates were evaluated for binding in a NanoBRET ligand-engagement screen across a library of HiBiT-tagged GPCRs. Employing this workflow, we generated several promiscuous fluorescent tracers that can effectively engage multiple GPCRs, demonstrating the efficiency of this approach. We believe that this workflow has the potential to accelerate discovery of NanoBRET fluorescent tracers for GPCRs and other target classes.

Detection of genome-edited and endogenously expressed G protein-coupled receptors

G protein-coupled receptors (GPCRs) are the largest family of membrane receptors and major targets for FDA-approved drugs. The ability to quantify GPCR expression and ligand binding characteristics in different cell types and tissues is therefore important for drug discovery. The advent of genome editing along with developments in fluorescent ligand design offers exciting new possibilities to probe GPCRs in their native environment. This review provides an overview of the recent technical advances employed to study the localisation and ligand binding characteristics of genome-edited and endogenously expressed GPCRs.

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.

Development and Application of Subtype-Selective Fluorescent Antagonists for the Study of the Human Adenosine A1 Receptor in Living Cells

The adenosine A₁ receptor (A₁AR) is a G-protein-coupled receptor (GPCR) that provides important therapeutic opportunities for a number of conditions including congestive heart failure, tachycardia, and neuropathic pain. The development of A₁AR-selective fluorescent ligands will enhance our understanding of the subcellular mechanisms underlying A₁AR pharmacology facilitating the development of more efficacious and selective therapies. Herein, we report the design, synthesis, and application of a novel series of A₁AR-selective fluorescent probes based on 8-functionalized bicyclo[2.2.2]octylxanthine and 3-functionalized 8-(adamant-1-yl) xanthine scaffolds. These fluorescent conjugates allowed quantification of kinetic and equilibrium ligand binding parameters using NanoBRET and visualization of specific receptor distribution patterns in living cells by confocal imaging and total internal reflection fluorescence (TIRF) microscopy. As such, the novel A₁AR-selective fluorescent antagonists described herein can be applied in conjunction with a series of fluorescence-based techniques to foster understanding of A₁AR molecular pharmacology and signaling in living cells.

Deciphering the Agonist Binding Mechanism to the Adenosine A1 Receptor

Despite being among the most characterized G protein-coupled receptors (GPCRs), adenosine receptors (ARs) have always been a difficult target in drug design. To date, no agonist other than the natural effector and the diagnostic regadenoson has been approved for human use. Recently, the structure of the adenosine A1 receptor (A1R) was determined in the active, Gi protein complexed state; this has important repercussions for structure-based drug design. Here, we employed supervised molecular dynamics simulations and mutagenesis experiments to extend the structural knowledge of the binding of selective agonists to A1R. Our results identify new residues involved in the association and dissociation pathway, they suggest the binding mode of N6-cyclopentyladenosine (CPA) related ligands, and they highlight the dramatic effect that chemical modifications can have on the overall binding mechanism, paving the way for the rational development of a structure-kinetics relationship of A1R agonists.

Development of a Testing Funnel for Identification of Small-Molecule Modulators Targeting Secretin Receptors

The secretin receptor (SCTR), a prototypical class B G protein-coupled receptor (GPCR), exerts its effects mainly by activating Gαs proteins upon binding of its endogenous peptide ligand secretin. SCTRs can be found in a variety of tissues and organs across species, including the pancreas, stomach, liver, heart, lung, colon, kidney, and brain. Beyond that, modulation of SCTR-mediated signaling has therapeutic potential for the treatment of multiple diseases, such as heart failure, obesity, and diabetes. However, no ligands other than secretin and its peptide analogs have been described to regulate SCTRs, probably due to inherent challenges in family B GPCR drug discovery. Here we report creation of a testing funnel that allowed targeted detection of SCTR small-molecule activators. Pursuing the strategy to identify positive allosteric modulators (PAMs), we established a unique primary screening assay employing a mixture of three orthosteric stimulators that was compared in a screening campaign testing 12,000 small-molecule compounds. Beyond that, we developed a comprehensive set of secondary assays, such as a radiolabel-free target engagement assay and a NanoBiT (NanoLuc Binary Technology)-based approach to detect β-arrestin-2 recruitment, all feasible in a high-throughput environment as well as capable of profiling ligands and hits regarding their effect on binding and receptor function. This combination of methods enabled the discovery of five promising scaffolds, four of which have been validated and further characterized with respect to their allosteric activities. We propose that our results may serve as starting points for developing the first in vivo active small molecules targeting SCTRs.

Molecular Dynamics Simulations in Drug Discovery and Pharmaceutical Development

Molecular dynamics (MD) simulations have become increasingly useful in the modern drug development process. In this review, we give a broad overview of the current application possibilities of MD in drug discovery and pharmaceutical development. Starting from the target validation step of the drug development process, we give several examples of how MD studies can give important insights into the dynamics and function of identified drug targets such as sirtuins, RAS proteins, or intrinsically disordered proteins. The role of MD in antibody design is also reviewed. In the lead discovery and lead optimization phases, MD facilitates the evaluation of the binding energetics and kinetics of the ligand-receptor interactions, therefore guiding the choice of the best candidate molecules for further development. The importance of considering the biological lipid bilayer environment in the MD simulations of membrane proteins is also discussed, using G-protein coupled receptors and ion channels as well as the drug-metabolizing cytochrome P450 enzymes as relevant examples. Lastly, we discuss the emerging role of MD simulations in facilitating the pharmaceutical formulation development of drugs and candidate drugs. Specifically, we look at how MD can be used in studying the crystalline and amorphous solids, the stability of amorphous drug or drug-polymer formulations, and drug solubility. Moreover, since nanoparticle drug formulations are of great interest in the field of drug delivery research, different applications of nano-particle simulations are also briefly summarized using multiple recent studies as examples. In the future, the role of MD simulations in facilitating the drug development process is likely to grow substantially with the increasing computer power and advancements in the development of force fields and enhanced MD methodologies.

Equilibrium and Kinetic Measurements of Ligand Binding to HiBiT-tagged GPCRs on the Surface of Living Cells

G-protein coupled receptors (GPCRs) remain at the forefront of drug discovery efforts. Detailed assessment of features contributing to GPCR ligand engagement in a physiologically relevant environment is imperative to the development of new therapeutics with improved efficacy. Traditionally, binding properties such as affinity and kinetics were obtained using biochemical radioligand binding assays. More recently, the high specificity of resonance energy transfer has been leveraged toward the development of homogeneous cell-based proximity assays with capacity for real-time kinetic measurements. This suite of ligand binding protocols couples the specificity of bioluminescent resonance energy transfer (BRET) with the sensitivity afforded by the luminescent HiBiT peptide. The BRET format is used to quantify dynamic interactions between ligands and their cognate HiBiT-tagged GPCRs through competitive binding with fluorescent Tracers. At the same time, high affinity complementation of HiBiT with the cell impermeable LgBiT limits the bright bioluminescence donor signal to the cell surface and eliminates luminescence background from unoccupied receptors present in intracellular compartments.

Biased Ligands Differentially Shape the Conformation of the Extracellular Loop Region in 5-HT2B Receptors

G protein-coupled receptors are linked to various intracellular transducers, each pathway associated with different physiological effects. Biased ligands, capable of activating one pathway over another, are gaining attention for their therapeutic potential, as they could selectively activate beneficial pathways whilst avoiding those responsible for adverse effects. We performed molecular dynamics simulations with known β-arrestin-biased ligands like lysergic acid diethylamide and ergotamine in complex with the 5-HT_2B receptor and discovered that the extent of ligand bias is directly connected with the degree of closure of the extracellular loop region. Given a loose allosteric coupling of extracellular and intracellular receptor regions, we delineate a concept for biased signaling at serotonin receptors, by which conformational interference with binding pocket closure restricts the signaling repertoire of the receptor. Molecular docking studies of biased ligands gathered from the BiasDB demonstrate that larger ligands only show plausible docking poses in the ergotamine-bound structure, highlighting the conformational constraints associated with bias. This emphasizes the importance of selecting the appropriate receptor conformation on which to base virtual screening workflows in structure-based drug design of biased ligands. As this mechanism of ligand bias has also been observed for muscarinic receptors, our studies provide a general mechanism of signaling bias transferable between aminergic receptors.