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.

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.

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.

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.

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.

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.

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.

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.

Controlling the Dissociation of Ligands from the Adenosine A2A Receptor through Modulation of Salt Bridge Strength.

The association and dissociation kinetics of ligands binding to proteins vary considerably, but the mechanisms behind this variability are poorly understood, limiting their utilization for drug discovery. This is particularly so for G protein-coupled receptors (GPCRs) where high resolution structural information is only beginning to emerge. Engineering the human A2A adenosine receptor has allowed structures to be solved in complex with the reference compound ZM241385 and four related ligands at high resolution. Differences between the structures are limited, with the most pronounced being the interaction of each ligand with a salt bridge on the extracellular side of the receptor. Mutagenesis experiments confirm the role of this salt bridge in controlling the dissociation kinetics of the ligands from the receptor, while molecular dynamics simulations demonstrate the ability of ligands to modulate salt bridge stability. These results shed light on a structural determinant of ligand dissociation kinetics and identify a means by which this property may be optimized.

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.

Probing GPCR structure: adenosine and P2Y nucleotide receptors.

The adenosine receptors (ARs) provide an example of how to accurately predict ligand recognition, even prior to the availability of a crystallographic structure. Homology modeling has been used to gain structural insight, in conjunction with site-directed mutagenesis, and structure-activity relationships of small molecular ligands. Recent X-ray structures greatly improved the accuracy of knowledge of AR ligand recognition and furthermore characterized conformational changes induced by receptor activation. Now, homology modeling extends these structural insights to related GPCRs and suggests new ligand structures. This strategy is also being applied to the eight subtypes of P2Y receptors for extracellular nucleotides, which lack X-ray structures and are best modeled by homology to the CXCR4 (peptide) receptor. Neoceptors, as studied for three of the four AR subtypes, create a molecular complementarity between a mutant receptor and a chemically tailored agonist ligand to selectively enhance affinity, implying direct physical contact and thus validating docking hypotheses.

Structural sweet spot for A1 adenosine receptor activation by truncated (N)-methanocarba nucleosides: receptor docking and potent anticonvulsant activity.

A(1) adenosine receptor (AR) agonists display antiischemic and antiepileptic neuroprotective activity, but peripheral cardiovascular side effects impeded their development. SAR study of N(6)-cycloalkylmethyl 4'-truncated (N)-methanocarba-adenosines identified 10 (MRS5474, N(6)-dicyclopropylmethyl, K(i) = 47.9 nM) as a moderately A(1)AR-selective full agonist. Two stereochemically defined N(6)-methynyl group substituents displayed narrow SAR; groups larger than cyclobutyl greatly reduced AR affinity, and those larger or smaller than cyclopropyl reduced A(1)AR selectivity. Nucleoside docking to A(1)AR homology model characterized distinct hydrophobic cyclopropyl subpockets, the larger "A" forming contacts with Thr270 (7.35), Tyr271 (7.36), Ile274 (7.39), and carbon chains of glutamates (EL2) and the smaller subpocket "B" forming contacts between TM6 and TM7. 10 suppressed minimal clonic seizures (6 Hz mouse model) without typical rotarod impairment of A(1)AR agonists. Truncated nucleosides, an appealing preclinical approach, have more druglike physicochemical properties than other A(1)AR agonists. Thus, we identified highly restricted regions for substitution around N(6) suitable for an A(1)AR agonist with anticonvulsant activity.

Farnesyl pyrophosphate is an endogenous antagonist to ADP-stimulated P2Y12 receptor-mediated platelet aggregation.

Farnesyl pyrophosphate (FPP) is an intermediate in cholesterol biosynthesis, and it has also been reported to activate platelet LPA (lysophosphatidic acid) receptors. The aim of this study was to investigate the role of extracellular FPP in platelet aggregation. Human platelets were studied with light transmission aggregometry, flow cytometry and [³5S]GTPγS binding assays. As shown previously, FPP could potentiate LPA-stimulated shape change. Surprisingly, FPP also acted as a selective insurmountable antagonist to ADP-induced platelet aggregation. FPP inhibited ADP-induced expression of P-selectin and the activated glycoprotein (Gp)IIb/IIIa receptor. FPP blocked ADP-induced inhibition of cAMP accumulation and [³5S]GTPγS binding in platelets. In Chinese hamster ovary cells expressing the P2Y12 receptor, FPP caused a rightward shift of the [³5S]GTPγS binding curve. In Sf9 insect cells expressing the human P2Y12 receptor, FPP showed a concentration-dependent, although incomplete inhibition of [³H]PSB-0413 binding. Docking of FPP in a P2Y12 receptor model revealed molecular similarities with ADP and a good fit into the binding pocket for ADP. In conclusion, FPP is an insurmountable antagonist of ADP-induced platelet aggregation mediated by the P2Y12 receptor. It could be an endogenous antithrombotic factor modulating the strong platelet aggregatory effects of ADP in a manner similar to the use of clopidogrel, prasugrel or ticagrelor in the treatment of ischaemic heart disease.

Structure-guided design of A(3) adenosine receptor-selective nucleosides: combination of 2-arylethynyl and bicyclo[3.1.0]hexane substitutions.

(N)-Methanocarba adenosine 5'-methyluronamides containing known A(3) AR (adenosine receptor)-enhancing modifications, i.e., 2-(arylethynyl)adenine and N(6)-methyl or N(6)-(3-substituted-benzyl), were nanomolar full agonists of human (h) A(3)AR and highly selective (K(i) ∼0.6 nM, N(6)-methyl 2-(halophenylethynyl) analogues 13 and 14). Combined 2-arylethynyl-N(6)-3-chlorobenzyl substitutions preserved A(3)AR affinity/selectivity in the (N)-methanocarba series (e.g., 3,4-difluoro full agonist MRS5698 31, K(i) 3 nM, human and mouse A(3)) better than that for ribosides. Polyaromatic 2-ethynyl N(6)-3-chlorobenzyl analogues, such as potent linearly extended 2-p-biphenylethynyl MRS5679 34 (K(i) hA(3) 3.1 nM; A(1), A(2A), inactive) and fluorescent 1-pyrene adduct MRS5704 35 (K(i) hA(3) 68.3 nM), were conformationally rigid; receptor docking identified a large, mainly hydrophobic binding region. The vicinity of receptor-bound C2 groups was probed by homology modeling based on recent X-ray structure of an agonist-bound A(2A)AR, with a predicted helical rearrangement requiring an agonist-specific outward displacement of TM2 resembling opsin. Thus, the X-ray structure of related A(2A)AR is useful in guiding the design of new A(3)AR agonists.

Optimization of adenosine 5'-carboxamide derivatives as adenosine receptor agonists using structure-based ligand design and fragment screening.

Structures of G protein-coupled receptors (GPCRs) have a proven utility in the discovery of new antagonists and inverse agonists modulating signaling of this important family of clinical targets. Applicability of active-state GPCR structures to virtual screening and rational optimization of agonists, however, remains to be assessed. In this study of adenosine 5' derivatives, we evaluated the performance of an agonist-bound A(2A) adenosine receptor (AR) structure in retrieval of known agonists and then employed the structure to screen for new fragments optimally fitting the corresponding subpocket. Biochemical and functional assays demonstrate high affinity of new derivatives that include polar heterocycles. The binding models also explain modest selectivity gain for some substituents toward the closely related A(1)AR subtype and the modified agonist efficacy of some of these ligands. The study suggests further applicability of in silico fragment screening to rational lead optimization in GPCRs.

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.

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.

Evaluation of molecular modeling of agonist binding in light of the crystallographic structure of an agonist-bound A2A adenosine receptor.

Molecular modeling of agonist binding to the human A(2A) adenosine receptor (AR) was assessed and extended in light of crystallographic structures. Heterocyclic adenine nitrogens of cocrystallized agonist overlaid corresponding positions of the heterocyclic base of a bound triazolotriazine antagonist, and ribose moiety was coordinated in a hydrophilic region, as previously predicted based on modeling using the inactive receptor. Automatic agonist docking of 20 known potent nucleoside agonists to agonist-bound A(2A)AR crystallographic structures predicted new stabilizing protein interactions to provide a structural basis for previous empirical structure activity relationships consistent with previous mutagenesis results. We predicted binding of novel C2 terminal amino acid conjugates of A(2A)AR agonist CGS21680 and used these models to interpret effects on binding affinity of newly synthesized agonists. d-Amino acid conjugates were generally more potent than l-stereoisomers and free terminal carboxylates more potent than corresponding methyl esters. Amino acid moieties were coordinated close to extracellular loops 2 and 3. Thus, molecular modeling is useful in probing ligand recognition and rational design of GPCR-targeting compounds with specific pharmacological profiles.