The platelet integrin ␣ IIb ␤ 3 is essential for hemostasis and thrombosis through its binding of adhesive plasma proteins. We have determined crystal structures of the

G protein-coupled receptors play a pivotal role in many physiological signaling pathways. Mounting evidence suggests that G protein-coupled receptors, including opioid receptors, form dimers, and dimerization is necessary for receptor maturation, signaling, and trafficking. However, the physiological role of dimerization in vivo has not been well-explored because of the lack of tools to study these dimers in endogenous systems. To address this problem, we previously generated antibodies to μ-δ opioid receptor (μOR-δOR) dimers and used them to study the pharmacology and signaling by this heteromer. We also showed that the heteromer exhibits restricted distribution in the brain and that its abundance is increased in response to chronic morphine administration. Thus, the μOR-δOR heteromer represents a potentially unique target for the development of therapeutics to treat pain. Here, we report the identification of compounds targeting μOR-δOR heteromers through high-throughput screening of a small-molecule library. These compounds exhibit activity in μOR-δOR cells but not μOR or δOR cells alone. Among them, CYM51010 was found to be a μOR-δOR–biased ligand, because its activity is blocked by the μOR-δOR heteromer antibody. Notably, systemic administration of CYM51010 induced antinociceptive activity similar to morphine, and chronic administration of CYM51010 resulted in lesser antinociceptive tolerance compared with morphine. Taken together, these results suggest that CYM51010, a μOR-δOR–biased ligand, could serve as a scaffold for the development of a unique type (heteromer-biased) of drug that is more potent and without the severe side effects associated with conventional clinical opioids.

The rapid growth of structural information for G-protein-coupled receptors (GPCRs) has led to a greater understanding of their structure, function, selectivity, and ligand binding. Although novel ligands have been identified using methods such as virtual screening, computationally driven lead optimization has been possible only in isolated cases because of challenges associated with predicting binding free energies for related compounds. Here, we provide a systematic characterization of the performance of free-energy perturbation (FEP) calculations to predict relative binding free energies of congeneric ligands binding to GPCR targets using a consistent protocol and no adjustable parameters. Using the FEP+ package, first we validated the protocol, which includes a full lipid bilayer and explicit solvent, by predicting the binding affinity for a total of 45 different ligands across four different GPCRs (adenosine A2AAR, β1 adrenergic, CXCR4 chemokine, and δ opioid receptors). Comparison with experimental binding affinity measurements revealed a highly predictive ranking correlation (average spearman ρ = 0.55) and low root-mean-square error (0.80 kcal/mol). Next, we applied FEP+ in a prospective project, where we predicted the affinity of novel, potent adenosine A2A receptor (A2AR) antagonists. Four novel compounds were synthesized and tested in a radioligand displacement assay, yielding affinity values in the nanomolar range. The affinity of two out of the four novel ligands (plus three previously reported compounds) was correctly predicted (within 1 kcal/mol), including one compound with approximately a tenfold increase in affinity compared to the starting compound. Detailed analyses of the simulations underlying the predictions provided insights into the structural basis for the two cases where the affinity was overpredicted. Taken together, these results establish a protocol for systematically applying FEP+ to GPCRs and provide guidelines for identifying potent molecules in drug discovery lead optimization projects.

BACKGROUND AND PURPOSEPM01183 is a new synthetic tetrahydroisoquinoline alkaloid that is currently in phase I clinical development for the treatment of solid tumours. In this study we have characterized the interactions of PM01183 with selected DNA molecules of defined sequence and its in vitro and in vivo cytotoxicity. EXPERIMENTAL APPROACHDNA binding characteristics of PM01183 were studied using electrophoretic mobility shift assays, fluorescence-based melting kinetic experiments and computational modelling methods. Its mechanism of action was investigated using flow cytometry, Western blot analysis and fluorescent microscopy. In vitro anti-tumour activity was determined by 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide assay and the in vivo activity utilized several human cancer models. KEY RESULTSElectrophoretic mobility shift assays demonstrated that PM01183 bound to DNA. Fluorescence-based thermal denaturation experiments showed that the most favourable DNA triplets providing a central guanine for covalent adduct formation are AGC, CGG, AGG and TGG. These binding preferences could be rationalized using molecular modelling. PM01183-DNA adducts in living cells give rise to double-strand breaks, triggering S-phase accumulation and apoptosis. The potent cytotoxic activity of PM01183 was ascertained in a 23-cell line panel with a mean GI50 value of 2.7 nM. In four murine xenograft models of human cancer, PM01183 inhibited tumour growth significantly with no weight loss of treated animals. CONCLUSIONS AND IMPLICATIONSPM01183 is shown to bind to selected DNA sequences and promoted apoptosis by inducing double-strand breaks at nanomolar concentrations. The potent anti-tumour activity of PM01183 in several murine models of human cancer supports its development as a novel anti-neoplastic agent. Abbreviations

Extensive experimental information supports the formation of ligand-specific conformations of G protein-coupled receptors (GPCRs) as a possible molecular basis for their functional selectivity for signaling pathways. Taking advantage of the recently published inactive and active crystal structures of GPCRs, we have implemented an all-atom computational strategy that combines different adaptive biasing techniques to identify ligand-specific conformations along pre-determined activation pathways. Using the prototypic GPCR β2-adrenergic receptor as a suitable test case for validation, we show that ligands with different efficacies (either inverse agonists, neutral antagonists, or agonists) modulate the free-energy landscape of the receptor by shifting the conformational equilibrium towards active or inactive conformations depending on their elicited physiological response. Notably, we provide for the first time a quantitative description of the thermodynamics of the receptor in an explicit atomistic environment, which accounts for the receptor basal activity and the stabilization of different active-like states by differently potent agonists. Structural inspection of these metastable states reveals unique conformations of the receptor that may have been difficult to retrieve experimentally.