Intracellular allosteric antagonism of the CCR9 receptor.

Chemokines and their G-protein-coupled receptors play a diverse role in immune defence by controlling the migration, activation and survival of immune cells. They are also involved in viral entry, tumour growth and metastasis and hence are important drug targets in a wide range of diseases. Despite very significant efforts by the pharmaceutical industry to develop drugs, with over 50 small-molecule drugs directed at the family entering clinical development, only two compounds have reached the market: maraviroc (CCR5) for HIV infection and plerixafor (CXCR4) for stem-cell mobilization. The high failure rate may in part be due to limited understanding of the mechanism of action of chemokine antagonists and an inability to optimize compounds in the absence of structural information. CC chemokine receptor type 9 (CCR9) activation by CCL25 plays a key role in leukocyte recruitment to the gut and represents a therapeutic target in inflammatory bowel disease. The selective CCR9 antagonist vercirnon progressed to phase 3 clinical trials in Crohn's disease but efficacy was limited, with the need for very high doses to block receptor activation. Here we report the crystal structure of the CCR9 receptor in complex with vercirnon at 2.8 Å resolution. Remarkably, vercirnon binds to the intracellular side of the receptor, exerting allosteric antagonism and preventing G-protein coupling. This binding site explains the need for relatively lipophilic ligands and describes another example of an allosteric site on G-protein-coupled receptors that can be targeted for drug design, not only at CCR9, but potentially extending to other chemokine receptors.

Accurate Prediction of GPCR Ligand Binding Affinity with Free Energy Perturbation.

The computational prediction of relative binding free energies is a crucial goal for drug discovery, and G protein-coupled receptors (GPCRs) are arguably the most important drug target class. However, they present increased complexity to model compared to soluble globular proteins. Despite breakthroughs, experimental X-ray crystal and cryo-EM structures are challenging to attain, meaning computational models of the receptor and ligand binding mode are sometimes necessary. This leads to uncertainty in understanding ligand-protein binding induced changes such as, water positioning and displacement, side chain positioning, hydrogen bond networks, and the overall structure of the hydration shell around the ligand and protein. In other words, the very elements that define structure activity relationships (SARs) and are crucial for accurate binding free energy calculations are typically more uncertain for GPCRs. In this work we use free energy perturbation (FEP) to predict the relative binding free energies for ligands of two different GPCRs. We pinpoint the key aspects for success such as the important role of key water molecules, amino acid ionization states, and the benefit of equilibration with specific ligands. Initial calculations following typical FEP setup and execution protocols delivered no correlation with experiment, but we show how results are improved in a logical and systematic way. This approach gave, in the best cases, a coefficient of determination (R2) compared with experiment in the range of 0.6-0.9 and mean unsigned errors compared to experiment of 0.6-0.7 kcal/mol. We anticipate that our findings will be applicable to other difficult-to-model protein ligand data sets and be of wide interest for the community to continue improving FE binding energy predictions.

X-Ray Crystallography and Free Energy Calculations Reveal the Binding Mechanism of A2A Adenosine Receptor Antagonists.

We present a robust protocol based on iterations of free energy perturbation (FEP) calculations, chemical synthesis, biophysical mapping and X-ray crystallography to reveal the binding mode of an antagonist series to the A2A adenosine receptor (AR). Eight A2A AR binding site mutations from biophysical mapping experiments were initially analyzed with sidechain FEP simulations, performed on alternate binding modes. The results distinctively supported one binding mode, which was subsequently used to design new chromone derivatives. Their affinities for the A2A AR were experimentally determined and investigated through a cycle of ligand-FEP calculations, validating the binding orientation of the different chemical substituents proposed. Subsequent X-ray crystallography of the A2A AR with a low and a high affinity chromone derivative confirmed the predicted binding orientation. The new molecules and structures here reported were driven by free energy calculations, and provide new insights on antagonist binding to the A2A AR, an emerging target in immuno-oncology.

Understanding Ligand Binding Selectivity in a Prototypical GPCR Family.

Adenosine receptors are involved in many pathological conditions and are thus promising drug targets. However, developing drugs that target this GPCR subfamily is a challenging task. A number of drug candidates fail due to lack of selectivity which results in unwanted side effects. The extensive structural similarity of adenosine receptors complicates the design of selective ligands. The problem of selective targeting is a general concern in GPCRs, and in this respect adenosine receptors are a prototypical example. Here we use enhanced sampling simulations to decipher the determinants of selectivity of ligands in A2a and A1 adenosine receptors. Our model shows how small differences in the binding pocket and in the water network around the ligand can be leveraged to achieve selectivity.

Decoding the Role of Water Dynamics in Ligand-Protein Unbinding: CRF1R as a Test Case.

The residence time of a ligand-protein complex is a crucial aspect in determining biological effect in vivo. Despite its importance, the prediction of ligand koff still remains challenging for modern computational chemistry. We have developed aMetaD, a fast and generally applicable computational protocol to predict ligand-protein unbinding events using a molecular dynamics (MD) method based on adiabatic-bias MD and metadynamics. This physics-based, fully flexible, and pose-dependent ligand scoring function evaluates the maximum energy (RTscore) required to move the ligand from the bound-state energy basin to the next. Unbinding trajectories are automatically analyzed and translated into atomic solvation factor (SF) values representing the water dynamics during the unbinding event. This novel computational protocol was initially tested on two M3 muscarinic receptor and two adenosine A2A receptor antagonists and then evaluated on a test set of 12 CRF1R ligands. The resulting RTscores were used successfully to classify ligands with different residence times. Additionally, the SF analysis was used to detect key differences in the degree of accessibility to water molecules during the predicted ligand unbinding events. The protocol provides actionable working hypotheses that are applicable in a drug discovery program for the rational optimization of ligand binding kinetics.

Alchemical Free Energy Calculations on Membrane-Associated Proteins.

Membrane proteins have diverse functions within cells and are well-established drug targets. The advances in membrane protein structural biology have revealed drug and lipid binding sites on membrane proteins, while computational methods such as molecular simulations can resolve the thermodynamic basis of these interactions. Particularly, alchemical free energy calculations have shown promise in the calculation of reliable and reproducible binding free energies of protein-ligand and protein-lipid complexes in membrane-associated systems. In this review, we present an overview of representative alchemical free energy studies on G-protein-coupled receptors, ion channels, transporters as well as protein-lipid interactions, with emphasis on best practices and critical aspects of running these simulations. Additionally, we analyze challenges and successes when running alchemical free energy calculations on membrane-associated proteins. Finally, we highlight the value of alchemical free energy calculations calculations in drug discovery and their applicability in the pharmaceutical industry.

Click modification in the N6 region of A3 adenosine receptor-selective carbocyclic nucleosides for dendrimeric tethering that preserves pharmacophore recognition.

Adenosine derivatives were modified with alkynyl groups on N(6) substituents for linkage to carriers using Cu(I)-catalyzed click chemistry. Two parallel series, both containing a rigid North-methanocarba (bicyclo[3.1.0]hexane) ring system in place of ribose, behaved as A(3) adenosine receptor (AR) agonists: (5'-methyluronamides) or partial agonists (4'-truncated). Terminal alkynyl groups on a chain at the 3 position of a N(6)-benzyl group or simply through a N(6)-propargyl group were coupled to azido derivatives, which included both small molecules and G4 (fourth-generation) multivalent poly(amidoamine) (PAMAM) dendrimers, to form 1,2,3-triazolyl linkers. The small molecular triazoles probed the tolerance in A(3)AR binding of distal, sterically bulky groups such as 1-adamantyl. Terminal 4-fluoro-3-nitrophenyl groups anticipated nucleophilic substitution for chain extension and (18)F radiolabeling. N(6)-(4-Fluoro-3-nitrophenyl)-triazolylmethyl derivative 32 displayed a K(i) of 9.1 nM at A(3)AR with ∼1000-fold subtype selectivity. Multivalent conjugates additionally containing click-linked water-solubilizing polyethylene glycol groups potently activated A(3)AR in the 5'-methyluronamide, but not 4' truncated series. N(6)-Benzyl nucleoside conjugate 43 (apparent K(i) 24 nM) maintained binding affinity of the monomer better than a N(6)-triazolylmethyl derivative. Thus, the N(6) region of 5'-methyluronamide derivatives, as modeled in receptor docking, is suitable for functionalization and tethering by click chemistry to achieve high A(3)AR agonist affinity and enhanced selectivity.

G protein-coupled adenosine (P1) and P2Y receptors: ligand design and receptor interactions.

The medicinal chemistry and pharmacology of the four subtypes of adenosine receptors (ARs) and the eight subtypes of P2Y receptors (P2YRs, activated by a range of purine and pyrimidine mono- and dinucleotides) has recently advanced significantly leading to selective ligands. X-ray crystallographic structures of both agonist- and antagonist-bound forms of the A(2A)AR have provided unprecedented three-dimensional detail concerning molecular recognition in the binding site and the conformational changes in receptor activation. It is apparent that this ubiquitous cell signaling system has implications for understanding and treating many diseases. ATP and other nucleotides are readily released from intracellular sources under conditions of injury and organ stress, such as hypoxia, ischemia, or mechanical stress, and through channels and vesicular release. Adenosine may be generated extracellularly or by cellular release. Therefore, depending on pathophysiological factors, in a given tissue, there is often a tonic activation of one or more of the ARs or P2YRs that can be modulated by exogenous agents for a beneficial effect. Thus, this field has provided fertile ground for pharmaceutical development, leading to clinical trials of selective receptor ligands as imaging agents or for conditions including cardiac arrhythmias, ischemia/reperfusion injury, diabetes, pain, thrombosis, Parkinson's disease, rheumatoid arthritis, psoriasis, dry eye disease, pulmonary diseases such as cystic fibrosis, glaucoma, cancer, chronic hepatitis C, and other diseases.

Novel Macrocyclic Antagonists of the CGRP Receptor Part 2: Stereochemical Inversion Induces an Unprecedented Binding Mode.

The diastereomeric macrocyclic calcitonin gene-related peptide (CGRP) antagonists HTL0029881 (3) and HTL0029882 (4), in which the stereochemistry of a spiro center is reversed, surprisingly demonstrate comparable potency. X-ray crystallographic characterization demonstrates that 3 binds to the CGRP receptor in a precedented manner but that 4 binds in an unprecedented, unexpected, and radically different manner. The observation of this phenomenon is noteworthy and may open novel avenues for CGRP receptor antagonist design.

Novel Macrocyclic Antagonists of the Calcitonin Gene-Related Peptide Receptor: Design, Realization, and Structural Characterization of Protein-Ligand Complexes.

A series of macrocyclic calcitonin gene-related peptide (CGRP) receptor antagonists identified using structure-based design principles, exemplified by HTL0028016 (1) and HTL0028125 (2), is described. Structural characterization by X-ray crystallography of the interaction of two of the macrocycle antagonists with the CGRP receptor ectodomain is described, along with structure-activity relationships associated with point changes to the macrocyclic antagonists. The identification of non-peptidic/natural product-derived, macrocyclic ligands for a G protein coupled receptor (GPCR) is noteworthy.

Pharmacochemistry of the platelet purinergic receptors.

Platelets contain at least five purinergic G protein-coupled receptors, e.g., the pro-aggregatory P2Y(1) and P2Y(12) receptors, a P2Y(14) receptor (GPR105) of unknown function, and anti-aggregatory A(2A) and A(2B) adenosine receptor (ARs), in addition to the ligand-gated P2X1 ion channel. Probing the structure-activity relationships (SARs) of the P2X and P2Y receptors for extracellular nucleotides has resulted in numerous new agonist and antagonist ligands. Selective agents derived from known ligands and novel chemotypes can be used to help define the subtypes pharmacologically. Some of these agents have entered into clinical trials in spite of the challenges of drug development for these classes of receptors. The functional architecture of P2 receptors was extensively explored using mutagenesis and molecular modeling, which are useful tools in drug discovery. In general, novel drug delivery methods, prodrug approaches, allosteric modulation, and biased agonism would be desirable to overcome side effects that tend to occur even with receptor subtype-selective ligands. Detailed SAR analyses have been constructed for nucleotide and non-nucleotide ligands at the P2Y(1), P2Y(12), and P2Y(14) receptors. The thienopyridine antithrombotic drugs Clopidogrel and Prasugrel require enzymatic pre-activation in vivo and react irreversibly with the P2Y(12) receptor. There is much pharmaceutical development activity aimed at identifying reversible P2Y(12) receptor antagonists. The screening of chemically diverse compound libraries has identified novel chemotypes that act as competitive, non-nucleotide antagonists of the P2Y(1) receptor or the P2Y(12) receptor, and antithrombotic properties of the structurally optimized analogues were demonstrated. In silico screening at the A(2A) AR has identified antagonist molecules having novel chemotypes. Fluorescent and other reporter groups incorporated into ligands can enable new technology for receptor assays and imaging. The A(2A) agonist CGS21680 and the P2Y(1) receptor antagonist MRS2500 were derivatized for covalent attachment to polyamidoamine dendrimeric carriers of MW 20,000, and the resulting multivalent conjugates inhibited ADP-promoted platelet aggregation. In conclusion, a wide range of new pharmacological tools is available to control platelet function by interacting with cell surface purine receptors.

Molecular probes for the A2A adenosine receptor based on a pyrazolo[4,3-e][1,2,4]triazolo[1,5-c]pyrimidin-5-amine scaffold.

Pyrazolo[4,3-e][1,2,4]triazolo[1,5-c]pyrimidin-5-amine derivatives such as SCH 442416 display high affinity and selectivity as antagonists for the human A(2A) adenosine receptor (AR). We extended ether-linked chain substituents at the p-position of the phenyl group using optimized O-alkylation. The conjugates included an ester, carboxylic acid and amines (for amide condensation), an alkyne (for click chemistry), a fluoropropyl group (for (18)F incorporation), and fluorophore reporter groups (e.g., BODIPY conjugate 14, K(i) 15 nM). The potent and A(2A)AR-selective N-aminoethylacetamide 7 and N-[2-(2-aminoethyl)-aminoethyl]acetamide 8 congeners were coupled to polyamidoamine (PAMAM) G3.5 dendrimers, and the multivalent conjugates displayed high A(2A)AR affinity. Theoretical docking of an AlexaFluor conjugate to the receptor X-ray structure highlighted the key interactions between the heterocyclic core and the binding pocket of the A(2A)AR as well as the distal anchoring of the fluorophore. In conclusion, we have synthesized a family of high affinity functionalized congeners as pharmacological probes for studying the A(2A)AR.

Comparison of three GPCR structural templates for modeling of the P2Y12 nucleotide receptor.

The P2Y(12) receptor (P2Y(12)R) is an ADP-activated G protein-coupled receptor (GPCR) that is an important target for antithrombotic drugs. Three homology models of P2Y(12)R were compared, based on different GPCR structural templates: bovine rhodopsin (bRHO), human A(2A) adenosine receptor (A(2A)AR), and human C-X-C chemokine receptor type 4 (CXCR4). By criteria of sequence analysis (25.6% identity in transmembrane region), deviation from helicity in the second transmembrane helix (TM2), docked poses of ligands highlighting the role of key residues, accessibility of a conserved disulfide bridge that is reactive toward irreversibly-binding antagonists, and the presence of a shared disulfide bridge between the third extracellular loop (EL3) and the N-terminus, the CXCR4-based model appeared to be the most consistent with known characteristics of P2Y(12)R. The docked poses of agonist 2MeSADP and charged anthraquinone antagonist PSB-0739 in the binding pocket of P2Y(12)R-CXC agree with previously published site-directed mutagenesis studies of Arg256 and Lys280. A sulfonate at position 2 of the anthraquinone core created a strong interaction with the Lys174(EL2) side chain. The docking poses of the irreversibly-binding, active metabolite (existing as two diastereoisomers in vivo) of the clinically utilized antagonist Clopidogrel were compared. The free thiol group of the 4S diastereoisomer, but not the 4R isomer, was found in close proximity (~4.7 Å) to the sulfur atom of a disulfide bridge involving Cys175, suggesting greater activity in covalent binding. Therefore, ligand docking to the CXCR4-based model of the P2Y(12)R predicted poses of both reversibly and irreversibly-binding small molecules, consistent with observed pharmacology and mutagenesis studies.

Synthesis and evaluation of 1,2,4-triazolo[1,5-c]pyrimidine derivatives as A2A receptor-selective antagonists.

Movement disorders such as Parkinson's disease and Huntington's disease are serious life-limiting and debilitating movement disorders. Their onset typically occurs from mid-life to late in life, and effective diagnostic techniques for detecting and following the disease course are lacking. Our goal is to develop receptor imaging agents for positron emission tomography (PET) that selectively target the most relevant subtype of adenosine receptors (AR) that are highly expressed in the striatum, that is, the A(2A) AR. To further this goal, we have synthesized and characterized pharmacologically a family of high affinity A(2A) AR ligands, based on the known antagonist, SCH 442416 (R=-Me), which have structural variability on the terminus (R=-Et, -i-Pr, -allyl, and others). A O-fluoroethyl analogue suitable for use as a PET tracer had a K(i) value of 12.4 nM and was highly selective for the A(2A) AR in comparison to the A(1) and A(3) ARs.

Structure-Based Drug Discovery of N-((R)-3-(7-Methyl-1H-indazol-5-yl)-1-oxo-1-(((S)-1-oxo-3-(piperidin-4-yl)-1-(4-(pyridin-4-yl)piperazin-1-yl)propan-2-yl)amino)propan-2-yl)-2'-oxo-1',2'-dihydrospiro[piperidine-4,4'-pyrido[2,3-d][1,3]oxazine]-1-carboxamide (HTL22562): A Calcitonin Gene-Related Peptide Receptor Antagonist for Acute Treatment of Migraine.

Structure-based drug design enabled the discovery of 8, HTL22562, a calcitonin gene-related peptide (CGRP) receptor antagonist. The structure of 8 complexed with the CGRP receptor was determined at a 1.6 Å resolution. Compound 8 is a highly potent, selective, metabolically stable, and soluble compound suitable for a range of administration routes that have the potential to provide rapid systemic exposures with resultant high levels of receptor coverage (e.g., subcutaneous). The low lipophilicity coupled with a low anticipated clinically efficacious plasma exposure for migraine also suggests a reduced potential for hepatotoxicity. These properties have led to 8 being selected as a clinical candidate for acute treatment of migraine.

Structural basis of the selectivity of the beta(2)-adrenergic receptor for fluorinated catecholamines.

The important and diverse biological functions of adrenergic receptors, a subclass of G protein-coupled receptors (GPCRs), have made the search for compounds that selectively stimulate or inhibit the activity of different adrenergic receptor subtypes an important area of medicinal chemistry. We previously synthesized 2-, 5-, and 6-fluoronorepinehprine (FNE) and 2-, 5-, and 6-fluoroepinephrine (FEPI) and found that 2FNE and 2FEPI were selective beta-adrenergic agonists and that 6FNE and 6FEPI were selective alpha-adrenergic agonists, while 5FNE and 5FEPI were unselective. Agonist potencies correlated well with receptor binding affinities. Here, through a combination of molecular modeling and site-directed mutagenesis, we have identified N293 in the beta(2)-adrenergic receptor as a crucial residue for the selectivity of the receptor for catecholamines fluorinated at different positions.

Understanding the structural and functional differences between mouse thyrotropin-releasing hormone receptors 1 and 2.

Multiple computational methods have been employed in a comparative study of thyrotropin-releasing hormone receptors 1 and 2 (TRH-R1 and TRH-R2) to explore the structural bases for the different functional properties of these G protein-coupled receptors. Three-dimensional models of both murine TRH receptors have been built and optimized by means of homology modeling based on the crystal structure of bovine rhodopsin, molecular dynamics simulations, and energy minimizations in a membrane-aqueous environment. The comparison between the two models showed a correlation between the higher flexibility and higher basal activity of TRH-R2 versus the lesser flexibility and lower basal activity of TRH-R1 and supported the involvement of the highly conserved W6.48 in the signaling process. A correlation between the level of basal activity and conformational changes of TM5 was detected also. Comparison between models of the wild type receptors and their W6.48A mutants, which have reversed basal activities compared with their respective wild types, further supported these correlations. A flexible molecular docking procedure revealed that TRH establishes a direct interaction with W6.48 in TRH-R2 but not in TRH-R1. We designed and performed new mutagenesis experiments that strongly supported these observations.

Structures of Human A1 and A2A Adenosine Receptors with Xanthines Reveal Determinants of Selectivity.

The adenosine A1 and A2A receptors belong to the purinergic family of G protein-coupled receptors, and regulate diverse functions of the cardiovascular, respiratory, renal, inflammation, and CNS. Xanthines such as caffeine and theophylline are weak, non-selective antagonists of adenosine receptors. Here we report the structure of a thermostabilized human A1 receptor at 3.3 Å resolution with PSB36, an A1-selective xanthine-based antagonist. This is compared with structures of the A2A receptor with PSB36 (2.8 Å resolution), caffeine (2.1 Å), and theophylline (2.0 Å) to highlight features of ligand recognition which are common across xanthines. The structures of A1R and A2AR were analyzed to identify the differences that are important selectivity determinants for xanthine ligands, and the role of T2707.35 in A1R (M2707.35 in A2AR) in conferring selectivity was confirmed by mutagenesis. The structural differences confirmed to lead to selectivity can be utilized in the design of new subtype-selective A1R or A2AR antagonists.

Impact of protein-ligand solvation and desolvation on transition state thermodynamic properties of adenosine A2A ligand binding kinetics.

Ligand-protein binding kinetic rates are growing in importance as parameters to consider in drug discovery and lead optimization. In this study we analysed using surface plasmon resonance (SPR) the transition state (TS) properties of a set of six adenosine A2A receptor inhibitors, belonging to both the xanthine and the triazolo-triazine scaffolds. SPR highlighted interesting differences among the ligands in the enthalpic and entropic components of the TS energy barriers for the binding and unbinding events. To better understand at a molecular level these differences, we developed suMetaD, a novel molecular dynamics (MD)-based approach combining supervised MD and metadynamics. This method allows simulation of the ligand unbinding and binding events. It also provides the system conformation corresponding to the highest energy barrier the ligand is required to overcome to reach the final state. For the six ligands evaluated in this study their TS thermodynamic properties were linked in particular to the role of water molecules in solvating/desolvating the pocket and the small molecules. suMetaD identified kinetic bottleneck conformations near the bound state position or in the vestibule area. In the first case the barrier is mainly enthalpic, requiring the breaking of strong interactions with the protein. In the vestibule TS location the kinetic bottleneck is instead mainly of entropic nature, linked to the solvent behaviour.

Ligand-based homology modeling as attractive tool to inspect GPCR structural plasticity.

G protein-coupled receptors (GPCRs) represent the largest family known of signal-transducing molecules. They convey signals for light and many extracellular regulatory molecules. GPCRs have been found to be dysfunctional/dysregulated in a growing number of human diseases and they have been estimated to be the targets of more than 40% of the drugs used in clinical medicine today. The crystal structure of rhodopsin provides the first three-dimensional GPCR information, which now supports homology modeling studies and structure-based drug design approaches. Here, we review our recent work on adenosine receptors, a family of GPCRs and, in particular, on A(3) adenosine receptor subtype antagonists. We will focus on an alternative approach to computationally explore the multi-conformational space of the antagonist-like state of the human A(3) receptor. We define ligand-based homology modeling as new approach to simulate the reorganization of the receptor induced by the ligand binding. The success of this approach is due to the synergic interaction between theory and experiment.