Exploring Biginelli-based scaffolds as A2B adenosine receptor antagonists: Unveiling novel structure-activity relationship trends, lead compounds, and potent colorectal anticancer agents

Antagonists of the A2B adenosine receptor have recently emerged as targeted anticancer agents and immune checkpoint inhibitors within the realm of cancer immunotherapy. This study presents a comprehensive evaluation of novel Biginelli-assembled pyrimidine chemotypes, including mono-, bi-, and tricyclic derivatives, as A2BAR antagonists. We conducted a comprehensive examination of the adenosinergic profile (both binding and functional) of a large compound library consisting of 168 compounds. This approach unveiled original lead compounds and enabled the identification of novel structure-activity relationship (SAR) trends, which were supported by extensive computational studies, including quantum mechanical calculations and free energy perturbation (FEP) analysis. In total, 25 molecules showed attractive affinity (Ki < 100 nM) and outstanding selectivity for A2BAR. From these, five molecules corresponding to the new benzothiazole scaffold were below the Ki < 10 nM threshold, in addition to a novel dual A2A/A2B antagonist. The most potent compounds, and the dual antagonist, showed enantiospecific recognition in the A2BAR. Two A2BAR selective antagonists and the dual A2AAR/A2BAR antagonist reported in this study were assessed for their impact on colorectal cancer cell lines. The results revealed a significant and dose-dependent reduction in cell proliferation. Notably, the A2BAR antagonists exhibited remarkable specificity, as they did not impede the proliferation of non-tumoral cell lines. These findings support the efficacy and potential that A2BAR antagonists as valuable candidates for cancer therapy, but also that they can effectively complement strategies involving A2AAR antagonism in the context of immune checkpoint inhibition.

N-(Heteroaryl)thiophene sulfonamides as angiotensin AT2 receptor ligands

Two series of N-(heteroaryl)thiophene sulfonamides, encompassing either a methylene imidazole group or a tert-butylimidazolylacetyl group in the meta position of the benzene ring, have been synthesized. An AT2R selective ligand with a Ki of 42 nM was identified in the first series and in the second series, six AT2R selective ligands with significantly improved binding affinities and Ki values of <5 nM were discovered. The binding modes to AT2R were explored by docking calculations combined with molecular dynamics simulations. Although some of the high affinity ligands exhibited fair stability in human liver microsomes, comparable to that observed with C21 undergoing clinical trials, most ligands displayed a very low metabolic stability with t½ of less than 10 min in human liver microsomes. The most promising ligand, with an AT2R Ki value of 4.9 nM and with intermediate stability in human hepatocytes (t½ = 77 min) caused a concentration-dependent vasorelaxation of pre-contracted mouse aorta.

Exploring the Effect of Halogenation in a Series of Potent and Selective A2B Adenosine Receptor Antagonists

The modulation of the A_2B adenosine receptor is a promising strategy in cancer (immuno) therapy, with A_2BAR antagonists emerging as immune checkpoint inhibitors. Herein, we report a systematic assessment of the impact of (di- and mono-)halogenation at positions 7 and/or 8 on both A_2BAR affinity and pharmacokinetic properties of a collection of A_2BAR antagonists and its study with structure-based free energy perturbation simulations. Monohalogenation at position 8 produced potent A_2BAR ligands irrespective of the nature of the halogen. In contrast, halogenation at position 7 and dihalogenation produced a halogen-size-dependent decay in affinity. Eight novel A_2BAR ligands exhibited remarkable affinity (K_i < 10 nM), exquisite subtype selectivity, and enantioselective recognition, with some eutomers eliciting sub-nanomolar affinity. The pharmacokinetic profile of representative derivatives showed enhanced solubility and microsomal stability. Finally, two compounds showed the capacity of reversing the antiproliferative effect of adenosine in activated primary human peripheral blood mononuclear cells.

Slow conformational dynamics of the human A2A adenosine receptor are temporally ordered

A more complete depiction of protein energy landscapes includes the identification of different function-related conformational states and the determination of the pathways connecting them. We used total internal reflection fluorescence (TIRF) imaging to investigate the conformational dynamics of the human A2A adenosine receptor (A2AAR), a class A G protein-coupled receptor (GPCR), at the single-molecule level. Slow, reversible conformational exchange was observed among three different fluorescence emission states populated for agonist-bound A2AAR. Transitions among these states predominantly occurred in a specific order, and exchange between inactive and active-like conformations proceeded through an intermediate state. Models derived from molecular dynamics simulations with available A2AAR structures rationalized the relative fluorescence emission intensities for the highest and lowest emission states but not the transition state. This suggests that the functionally critical intermediate state required to achieve activation is not currently visualized among available A2AAR structures.

Exploring Non-orthosteric Interactions with a Series of Potent and Selective A3 Antagonists

A library of potent and highly A₃AR selective pyrimidine-based compounds was designed to explore non-orthosteric interactions within this receptor. Starting from a prototypical orthosteric A₃AR antagonist (ISVY130), the structure-based design explored functionalized residues at the exocyclic amide L1 region and aimed to provide additional interactions outside the A₃AR orthosteric site. The novel ligands were assembled through an efficient and succinct synthetic approach, resulting in compounds that retain the A₃AR potent and selective profile while improving the solubility of the original scaffold. The experimentally demonstrated tolerability of the L1 region to structural functionalization was further assessed by molecular dynamics simulations, giving hints of the non-orthosteric interactions explored by these series. The results pave the way to explore newly functionalized A₃AR ligands, including covalent drugs and molecular probes for diagnostic and delivery purposes.

Optimization of 2-Amino-4,6-diarylpyrimidine-5-carbonitriles as Potent and Selective A1 Antagonists

We herein document a large collection of 108 2-amino-4,6-disubstituted-pyrimidine derivatives as potent, structurally simple, and highly selective A₁AR ligands. The most attractive ligands were confirmed as antagonists of the canonical cyclic adenosine monophosphate pathway, and some pharmacokinetic parameters were preliminarilly evaluated. The library, built through a reliable and efficient three-component reaction, comprehensively explored the chemical space allowing the identification of the most prominent features of the structure–activity and structure–selectivity relationships around this scaffold. These included the influence on the selectivity profile of the aromatic residues at positions R⁴ and R⁶ of the pyrimidine core but most importantly the prominent role to the unprecedented A₁AR selectivity profile exerted by the methyl group introduced at the exocyclic amino group. The structure–activity relationship trends on both A₁ and A_2AARs were conveniently interpreted with rigorous free energy perturbation simulations, which started from the receptor-driven docking model that guided the design of these series.

Deciphering conformational selectivity in the A2A adenosine G protein-coupled receptor by free energy simulations

Transmembranal G Protein-Coupled Receptors (GPCRs) transduce extracellular chemical signals to the cell, via conformational change from a resting (inactive) to an active (canonically bound to a G-protein) conformation. Receptor activation is normally modulated by extracellular ligand binding, but mutations in the receptor can also shift this equilibrium by stabilizing different conformational states. In this work, we built structure-energetic relationships of receptor activation based on original thermodynamic cycles that represent the conformational equilibrium of the prototypical A_2A adenosine receptor (AR). These cycles were solved with efficient free energy perturbation (FEP) protocols, allowing to distinguish the pharmacological profile of different series of A_2AAR agonists with different efficacies. The modulatory effects of point mutations on the basal activity of the receptor or on ligand efficacies could also be detected. This methodology can guide GPCR ligand design with tailored pharmacological properties, or allow the identification of mutations that modulate receptor activation with potential clinical implications.

The design of new ligands as chemical modulators of G protein-coupled receptors (GPCRs) has benefited considerably during the last years of advances in both the structural and computational biology disciplines. Within the last area, the use of free energy calculation methods has arisen as a computational tool to predict ligand affinities to explain structure-affinity relationships and guide lead optimization campaigns. However, our comprehension of the structural determinants of ligands with different pharmacological profile is scarce, and knowledge of the chemical modifications associated with an agonistic or antagonistic profile would be extremely valuable. We herein report an original implementation of the thermodynamic cycles associated with free energy perturbation (FEP) simulations, to mimic the conformational equilibrium between active and inactive GPCRs, and establish a framework to describe pharmacological profiles as a function of the ligands selectivity for a given receptor conformation. The advantage of this method resides into its simplicity of use, and the only consideration of active and inactive conformations of the receptor, with no simulation of the transitions between them. This model can accurately predict the pharmacological profile of series of full and partial agonists as opposed to antagonists of the A_2A adenosine receptor, and moreover, how certain mutations associated with modulation of basal activity can influence this pharmacological profiles, which enables our understanding of such clinically relevant mutations.

Potent Inhibition of Nicotinamide N-Methyltransferase by Alkene-Linked Bisubstrate Mimics Bearing Electron Deficient Aromatics

Nicotinamide N-methyltransferase (NNMT) methylates nicotinamide (vitamin B3) to generate 1-methylnicotinamide (MNA). NNMT overexpression has been linked to a variety of diseases, most prominently human cancers, indicating its potential as a therapeutic target. The development of small-molecule NNMT inhibitors has gained interest in recent years, with the most potent inhibitors sharing structural features based on elements of the nicotinamide substrate and the S-adenosyl-l-methionine (SAM) cofactor. We here report the development of new bisubstrate inhibitors that include electron-deficient aromatic groups to mimic the nicotinamide moiety. In addition, a trans-alkene linker was found to be optimal for connecting the substrate and cofactor mimics in these inhibitors. The most potent NNMT inhibitor identified exhibits an IC₅₀ value of 3.7 nM, placing it among the most active NNMT inhibitors reported to date. Complementary analytical techniques, modeling studies, and cell-based assays provide insights into the binding mode, affinity, and selectivity of these inhibitors.

Design, synthesis, HER2 inhibition and anticancer evaluation of new substituted 1,5-dihydro-4,1-benzoxazepines

A series of 11 new substituted 1,5-dihydro-4,1-benzoxazepine derivatives was synthesised to study the influence of the methyl group in the 1-(benzenesulphonyl) moiety, the replacement of the purine by the benzotriazole bioisosteric analogue, and the introduction of a bulky substituent at position 6 of the purine, on the biological effects. Their inhibition against isolated HER2 was studied and the structure–activity relationships have been confirmed by molecular modelling studies. The most potent compound against isolated HER2 is 9a with an IC₅₀ of 7.31 µM. We have investigated the effects of the target compounds on cell proliferation. The most active compound (7c) against all the tumour cell lines studied (IC₅₀ 0.42–0.86 µM) does not produce any modification in the expression of pro-caspase 3, but increases the caspase 1 expression, and promotes pyroptosis.

Potent and Subtype-Selective Dopamine D2 Receptor Biased Partial Agonists Discovered via an Ugi-Based Approach

Using a previously unexplored, efficient, and versatile multicomponent method, we herein report the rapid generation of novel potent and subtype-selective DRD₂ biased partial agonists. This strategy exemplifies the search for diverse and previously unexplored moieties for the secondary/allosteric pharmacophore of the common phenyl-piperazine scaffold. The pharmacological characterization of the new compound series led to the identification of several ligands with excellent DRD₂ affinity and subtype selectivity and remarkable functional selectivity for either the cAMP (22a and 24d) or the β-arrestin (27a and 29c) signaling pathways. These results were further interpreted on the basis of molecular models of these ligands in complex with the recent DRD₂ crystal structures, highlighting the critical role of the secondary/allosteric pharmacophore in modulating the functional selectivity profile.

Structural Basis of Inhibition of Human Insulin-Regulated Aminopeptidase (IRAP) by Benzopyran-Based Inhibitors

Inhibition of the insulin-regulated aminopeptidase (IRAP) improves memory and cognition in animal models. The enzyme has recently been crystallized and several series of inhibitors reported. We herein focused on one series of benzopyran-based inhibitors of IRAP known as the HFI series, with unresolved binding mode to IRAP, and developed a robust computational model to explain the structure-activity relationship (SAR) and potentially guide their further optimization. The binding model here proposed places the benzopyran ring in the catalytic binding site, coordinating the Zn²⁺ ion through the oxygen in position 3, in contrast to previous hypothesis. The whole series of HFI compounds was then systematically simulated, starting from this binding mode, using molecular dynamics and binding affinity estimated with the linear interaction energy (LIE) method. The agreement with experimental affinities supports the binding mode proposed, which was further challenged by rigorous free energy perturbation (FEP) calculations. Here, we found excellent correlation between experimental and calculated binding affinity differences, both between selected compound pairs and also for recently reported experimental data concerning the site directed mutagenesis of residue Phe544. The computationally derived structure-activity relationship of the HFI series and the understanding of the involvement of Phe544 in the binding of this scaffold provide valuable information for further lead optimization of novel IRAP inhibitors.

Free Energy Calculations for Protein-Ligand Binding Prediction

Computational prediction of protein-ligand binding involves initial determination of the binding mode and subsequent evaluation of the strength of the protein-ligand interactions, which directly correlates with ligand binding affinities. As a consequence of increasing computer power, rigorous approaches to calculate protein-ligand binding affinities, such as free energy perturbation (FEP) methods, are becoming an essential part of the toolbox of computer-aided drug design. In this chapter, we provide a general overview of these methods and introduce the QFEP modules, which are open-source API workflows based on our molecular dynamics (MD) package Q. The module QligFEP allows estimation of relative binding affinities along ligand series, while QresFEP is a module to estimate binding affinity shifts caused by single-point mutations of the protein. We herein provide guidelines for the use of each of these modules based on data extracted from ligand-design projects. While these modules are stand-alone, the combined use of the two workflows in a drug-design project yields complementary perspectives of the ligand binding problem, providing two sides of the same coin. The selected case studies illustrate how to use QFEP to approach the two key questions associated with ligand binding prediction: identifying the most favorable binding mode from different alternatives and establishing structure-affinity relationships that allow the rational optimization of hit compounds.

3,4-Dihydropyrimidin-2(1H)-ones as Antagonists of the Human A2B Adenosine Receptor: Optimization, Structure-Activity Relationship Studies, and Enantiospecific Recognition

We present and thoroughly characterize a large collection of 3,4-dihydropyrimidin-2(1H)-ones as A_2BAR antagonists, an emerging strategy in cancer (immuno) therapy. Most compounds selectively bind A_2BAR, with a number of potent and selective antagonists further confirmed by functional cyclic adenosine monophosphate experiments. The series was analyzed with one of the most exhaustive free energy perturbation studies on a GPCR, obtaining an accurate model of the structure-activity relationship of this chemotype. The stereospecific binding modeled for this scaffold was confirmed by resolving the two most potent ligands [(±)-47, and (±)-38 K_i = 10.20 and 23.6 nM, respectively] into their two enantiomers, isolating the affinity on the corresponding (S)-eutomers (K_i = 6.30 and 11.10 nM, respectively). The assessment of the effect in representative cytochromes (CYP3A4 and CYP2D6) demonstrated insignificant inhibitory activity, while in vitro experiments in three prostate cancer cells demonstrated that this pair of compounds exhibits a pronounced antimetastatic effect.

N-(Methyloxycarbonyl)thiophene sulfonamides as high affinity AT2 receptor ligands

A series of meta-substituted acetophenone derivatives, encompassing N-(alkyloxycarbonyl)thiophene sulfonamide fragments have been synthesized. Several selective AT2 receptor ligands were identified, among those a tert-butylimidazole derivative (20) with a K_i of 9.3 nM, that demonstrates a high stability in human liver microsomes (t_½ = 62 min) and in human hepatocytes (t_½ = 194 min). This methyloxycarbonylthiophene sulfonamide is a 20-fold more potent binder to the AT2 receptor and is considerably more stable in human liver microsomes, than a previously reported and broadly studied structurally related AT₂R prototype antagonist 3 (C38). Ligand 20 acts as an AT₂R agonist and caused an AT₂R mediated concentration-dependent vasorelaxation of pre-contracted mouse aorta. Furthermore, in contrast to imidazole derivative C38, the tert-butylimidazole derivative 20 is a poor inhibitor of CYP3A4, CYP2D6 and CYP2C9. It is demonstrated herein that smaller alkyloxycarbonyl groups make the ligands in this series of AT₂R selective compounds less prone to degradation and that a high AT2 receptor affinity can be retained after truncation of the alkyloxycarbonyl group. Binding modes of the most potent AT₂R ligands were explored by docking calculations combined with molecular dynamics simulations.

Nitrogen-Walk Approach to Explore Bioisosteric Replacements in a Series of Potent A2B Adenosine Receptor Antagonists

A systematic exploration of bioisosteric replacements for furan and thiophene cores in a series of potent A_2BAR antagonists has been carried out using the nitrogen-walk approach. A collection of 42 novel alkyl 4-substituted-2-methyl-1,4-dihydrobenzo[4,5]imidazo[1,2-a]pyrimidine-3-carboxylates, which contain 18 different pentagonal heterocyclic frameworks at position 4, was synthesized and evaluated. This study enabled the identification of new ligands that combine remarkable affinity (K_i < 30 nM) and exquisite selectivity. The structure-activity relationship (SAR) trends identified were substantiated by a molecular modeling study, based on a receptor-driven docking model and including a systematic free energy perturbation (FEP) study. Preliminary evaluation of the CYP3A4 and CYP2D6 inhibitory activity in optimized ligands evidenced weak and negligible activity, respectively. The stereospecific interaction between hA_2BAR and the eutomer of the most attractive novel antagonist (S)-18g (K_i = 3.66 nM) was validated.

Evolution of Angiotensin Peptides and Peptidomimetics as Angiotensin II Receptor Type 2 (AT2) Receptor Agonists

Angiotensin II receptor type 1 and 2 (AT1R and AT2R) are two G-protein coupled receptors that mediate most biological functions of the octapeptide Angiotensin II (Ang II). AT2R is upregulated upon tissue damage and its activation by selective AT2R agonists has become a promising approach in the search for new classes of pharmaceutical agents. We herein analyzed the chemical evolution of AT2R agonists starting from octapeptides, through shorter peptides and peptidomimetics to the first drug-like AT2R-selective agonist, C21, which is in Phase II clinical trials and aimed for idiopathic pulmonary fibrosis. Based on the recent crystal structures of AT1R and AT2R in complex with sarile, we identified a common binding model for a series of 11 selected AT2R agonists, consisting of peptides and peptidomimetics of different length, affinity towards AT2R and selectivity versus AT1R. Subsequent molecular dynamics simulations and free energy perturbation (FEP) calculations of binding affinities allowed the identification of the bioactive conformation and common pharmacophoric points, responsible for the key interactions with the receptor, which are maintained by the drug-like agonists. The results of this study should be helpful and facilitate the search for improved and even more potent AT2R-selective drug-like agonists.

Synthesis, Evaluation and Proposed Binding Pose of Substituted Spiro-Oxindole Dihydroquinazolinones as IRAP Inhibitors

Insulin-regulated aminopeptidase (IRAP) is a new potential macromolecular target for drugs aimed for treatment of cognitive disorders. Inhibition of IRAP by angiotensin IV (Ang IV) improves the memory and learning in rats. The majority of the known IRAP inhibitors are peptidic in character and suffer from poor pharmacokinetic properties. Herein, we present a series of small non-peptide IRAP inhibitors derived from a spiro-oxindole dihydroquinazolinone screening hit (pIC50 5.8). The compounds were synthesized either by a simple microwave (MW)-promoted three-component reaction, or by a two-step one-pot procedure. For decoration of the oxindole ring system, rapid MW-assisted Suzuki-Miyaura cross-couplings (1 min) were performed. A small improvement of potency (pIC50 6.6 for the most potent compound) and an increased solubility could be achieved. As deduced from computational modelling and MD simulations it is proposed that the S-configuration of the spiro-oxindole dihydroquinazolinones accounts for the inhibition of IRAP.

Macrocyclic peptidomimetics as inhibitors of insulin-regulated aminopeptidase (IRAP)

Macrocyclic analogues of the linear hexapeptide, angiotensin IV (AngIV) have proved to be potent inhibitors of insulin-regulated aminopeptidase (IRAP, oxytocinase, EC 3.4.11.3). Along with higher affinity, macrocycles may also offer better metabolic stability, membrane permeability and selectivity, however predicting the outcome of particular cycle modifications is challenging. Here we describe the development of a series of macrocyclic IRAP inhibitors with either disulphide, olefin metathesis or lactam bridges and variations of ring size and other functionality. The binding mode of these compounds is proposed based on molecular dynamics analysis. Estimation of binding affinities (ΔG) and relative binding free energies (ΔΔG) with the linear interaction energy (LIE) method and free energy perturbation (FEP) method showed good general agreement with the observed inhibitory potency. Experimental and calculated data highlight the cumulative importance of an intact N-terminal peptide, the specific nature of the macrocycle, the phenolic oxygen and the C-terminal functionality.

Inhibition of translation termination by small molecules targeting ribosomal release factors

The bacterial ribosome is an important drug target for antibiotics that can inhibit different stages of protein synthesis. Among the various classes of compounds that impair translation there are, however, no known small-molecule inhibitors that specifically target ribosomal release factors (RFs). The class I RFs are essential for correct termination of translation and they differ considerably between bacteria and eukaryotes, making them potential targets for inhibiting bacterial protein synthesis. We carried out virtual screening of a large compound library against 3D structures of free and ribosome-bound RFs in order to search for small molecules that could potentially inhibit termination by binding to the RFs. Here, we report identification of two such compounds which are found both to bind free RFs in solution and to inhibit peptide release on the ribosome, without affecting peptide bond formation.

Functional characterization in vitro of twelve naturally occurring variants of the human pancreatic polypeptide receptor NPY4R

Obesity has become a global health problem and therefore understanding of the mechanisms regulating hunger and satiety is of utmost importance for the development of new treatment strategies. The Y4 receptor, encoded by the NPY4R gene, and its ligand pancreatic polypeptide (PP) have been reported to mediate a satiety signal. Multiple genetic studies have reported an association between NPY4R copy number and body weight. The gene also displays several SNP variants, many of which lead to amino acid differences, making it interesting to study. We have investigated the functional properties of 12 naturally occurring amino acid sequence variants of the Y4 and interpret the results in relation to sequence conservation and our structural model of the human Y4 receptor protein. Three receptor variants, Cys201^ECL2Tyr, Val271^6.41Leu and Asn318^7.49Asp, were found to completely lose functional response, measured as inositol phosphate turnover, while retaining membrane expression. They display high sequence conservation and have important roles in the receptor structure. For two receptor variants the potency of PP was significantly decreased, Cys34^NTSer (EC50 = 2.9 nM, p < .001) and Val135^3.46Met (EC50 = 3.0 nM, p < .01), compared to wild-type Y4 (EC50 = 0.68 nM). Cys34 forms a disulphide bond with Cys298, linking the N-terminal part to ECL3. The Val135^3.46Met variant has an amino acid replacement located in the TM3 helix, one helix turn above the highly conserved ERH motif. This position has influence on the network of residues involved in receptor activation and subsequent inactivation. Sequence conservation and the structural model are consistent with these results. The remaining seven positions had no significant effect on the receptor's functional response compared to wild-type Y4. These positions display more variation during evolution. Understanding of the interactions between the Y4 receptor and its native PP agonist and the effects of amino acid variation on its functional response will hopefully lead to future therapeutic possibilities.

QligFEP: an automated workflow for small molecule free energy calculations in Q

The process of ligand binding to a biological target can be represented as the equilibrium between the relevant solvated and bound states of the ligand. This which is the basis of structure-based, rigorous methods such as the estimation of relative binding affinities by free energy perturbation (FEP). Despite the growing capacity of computing power and the development of more accurate force fields, a high throughput application of FEP is currently hampered due to the need, in the current schemes, of an expert user definition of the "alchemical" transformations between molecules in the series explored. Here, we present QligFEP, a solution to this problem using an automated workflow for FEP calculations based on a dual topology approach. In this scheme, the starting poses of each of the two ligands, for which the relative affinity is to be calculated, are explicitly present in the MD simulations associated with the (dual topology) FEP transformation, making the perturbation pathway between the two ligands univocal. We show that this generalized method can be applied to accurately estimate solvation free energies for amino acid sidechain mimics, as well as the binding affinity shifts due to the chemical changes typical of lead optimization processes. This is illustrated in a number of protein systems extracted from other FEP studies in the literature: inhibitors of CDK2 kinase and a series of A2A adenosine G protein-coupled receptor antagonists, where the results obtained with QligFEP are in excellent agreement with experimental data. In addition, our protocol allows for scaffold hopping perturbations to identify the binding affinities between different core scaffolds, which we illustrate with a series of Chk1 kinase inhibitors. QligFEP is implemented in the open-source MD package Q, and works with the most common family of force fields: OPLS, CHARMM and AMBER.

Structural mechanism of AadA, a dual-specificity aminoglycoside adenylyltransferase from Salmonella enterica

Streptomycin and spectinomycin are antibiotics that bind to the bacterial ribosome and perturb protein synthesis. The clinically most prevalent bacterial resistance mechanism is their chemical modification by aminoglycoside-modifying enzymes such as aminoglycoside nucleotidyltransferases (ANTs). AadA from Salmonella enterica is an aminoglycoside (3″)(9) adenylyltransferase that O-adenylates position 3″ of streptomycin and position 9 of spectinomycin. We previously reported the apo-AadA structure with a closed active site. To clarify how AadA binds ATP and its two chemically distinct drug substrates, we here report crystal structures of WT AadA complexed with ATP, magnesium, and streptomycin and of an active-site mutant, E87Q, complexed with ATP and streptomycin or the closely related dihydrostreptomycin. These structures revealed that ATP binding induces a conformational change that positions the two domains for drug binding at the interdomain cleft and disclosed the interactions between both domains and the three rings of streptomycin. Spectinomycin docking followed by molecular dynamics simulations suggested that, despite the limited structural similarities with streptomycin, spectinomycin makes similar interactions around the modification site and, in agreement with mutational data, forms critical interactions with fewer residues. Using structure-guided sequence analyses of ANT(3″)(9) enzymes acting on both substrates and ANT(9) enzymes active only on spectinomycin, we identified sequence determinants for activity on each substrate. We experimentally confirmed that Trp-173 and Asp-178 are essential only for streptomycin resistance. Activity assays indicated that Glu-87 is the catalytic base in AadA and that the nonadenylating E87Q mutant can hydrolyze ATP in the presence of streptomycin.

Structural Basis of Inhibition of Human Insulin-Regulated Aminopeptidase (IRAP) by Aryl Sulfonamides

The insulin-regulated aminopeptidase (IRAP) is a membrane-bound zinc metallopeptidase with many important regulatory functions. It has been demonstrated that inhibition of IRAP by angiotensin IV (Ang IV) and other peptides, as well as more druglike inhibitors, improves cognition in several rodent models. We recently reported a series of aryl sulfonamides as small-molecule IRAP inhibitors and a promising scaffold for pharmacological intervention. We have now expanded with a number of derivatives, report their stability in liver microsomes, and characterize the activity of the whole series in a new assay performed on recombinant human IRAP. Several compounds, such as the new fluorinated derivative 29, present submicromolar affinity and high metabolic stability. Starting from the two binding modes previously proposed for the sulfonamide scaffold, we systematically performed molecular dynamics simulations and binding affinity estimation with the linear interaction energy method for the full compound series. The significant agreement with experimental affinities suggests one of the binding modes, which was further confirmed by the excellent correlation for binding affinity differences between the selected pair of compounds obtained by rigorous free energy perturbation calculations. The new experimental data and the computationally derived structure-activity relationship of the sulfonamide series provide valuable information for further lead optimization of novel IRAP inhibitors.

Molecular Mechanisms in the Selectivity of Nonsteroidal Anti-Inflammatory Drugs

Nonsteroidal anti-inflammatory drugs (NSAIDs) inhibit cyclooxygenase (COX) 1 and 2 with varying degrees of selectivity. A group of COX-2 selective inhibitors-coxibs-binds in a time-dependent manner through a three-step mechanism, utilizing a side pocket in the binding site. Coxibs have been extensively probed to identify the structural features regulating the slow tight-binding mechanism responsible for COX-2 selectivity. In this study, we further probe a structurally and kinetically diverse data set of COX inhibitors in COX-2 by molecular dynamics and free energy simulations. We find that the features regulating the high affinities associated with time-dependency in COX depend on the inhibitor kinetics. In particular, most time-dependent inhibitors share a common structural binding mechanism, involving an induced-fit rotation of the side-chain of Leu531 in the main binding pocket. The high affinities of two-step slow tight-binding inhibitors and some slow reversible inhibitors can thus be explained by the increased space in the main binding pocket after this rotation. Coxibs that belong to a separate class of slow tight-binding inhibitors benefit more from the displacement of the neighboring side-chain of Arg513, exclusive to the COX-2 side-pocket. This displacement further stabilizes the aforementioned rotation of Leu531 and can explain the selectivity of coxibs for COX-2.

Elucidation of the Binding Mode of the Carboxyterminal Region of Peptide YY to the Human Y2 Receptor

Understanding the agonist-receptor interactions in the neuropeptide Y (NPY)/peptide YY (PYY) signaling system is fundamental for the design of novel modulators of appetite regulation. We report here the results of a multidisciplinary approach to elucidate the binding mode of the native peptide agonist PYY to the human Y₂ receptor, based on computational modeling, peptide chemistry and in vitro pharmacological analyses. The preserved binding orientation proposed for full-length PYY and five analogs, truncated at the amino terminus, explains our pharmacological results where truncations of the N-terminal proline helix showed little effect on peptide affinity. This was followed by receptor mutagenesis to investigate the roles of several receptor positions suggested by the modeling. As a complement, PYY-(3-36) analogs were synthesized with modifications at different positions in the common PYY/NPY C-terminal fragment (³²TRQRY³⁶-amide). The results were assessed and interpreted by molecular dynamics and Free Energy Perturbation (FEP) simulations of selected mutants, providing a detailed map of the interactions of the PYY/NPY C-terminal fragment with the transmembrane cavity of the Y₂ receptor. The amidated C-terminus would be stabilized by polar interactions with Gln288^6.55 and Tyr219^5.39, while Gln130^3.32 contributes to interactions with Q³⁴ in the peptide and T³² is close to the tip of TM7 in the receptor. This leaves the core, α-helix of the peptide exposed to make potential interactions with the extracellular loops. This model agrees with most experimental data available for the Y₂ system and can be used as a basis for optimization of Y₂ receptor agonists.

Characterization of Ligand Binding to GPCRs Through Computational Methods

The recent increase in available G protein-coupled receptor structures now contributes decisively to the structure-based ligand design. In this context, computational approaches in combination with medicinal chemistry and pharmacology are extremely helpful. Here, we provide an update on our structure-based computational protocols, used to answer key questions related to GPCR-ligand binding. All combined, these techniques can shed light on ligand binding modes, determine the molecular basis of conformational selection, for agonists and antagonists, as well as of subtype selectivity. To illustrate each of these questions, we will consider examples from existing projects on three families of class A (rhodopsin-like) GPCRs: one small-molecule (nucleotide-like) family, i.e., the adenosine receptors, and two peptide-binding receptors: neuropeptide-Y and angiotensin II receptors. The successful application of the same computational protocols to investigate this diverse group of receptor families gives an idea of the general applicability of our methodology in the characterization of GPCR-ligand binding.

The GPR139 reference agonists 1a and 7c, and tryptophan and phenylalanine share a common binding site

GPR139 is an orphan G protein-coupled receptor expressed in the brain, in particular in the habenula, hypothalamus and striatum. It has therefore been suggested that GPR139 is a possible target for metabolic disorders and Parkinson's disease. Several surrogate agonist series have been published for GPR139. Two series published by Shi et al. and Dvorak et al. included agonists 1a and 7c respectively, with potencies in the ten-nanomolar range. Furthermore, Isberg et al. and Liu et al. have previously shown that tryptophan (Trp) and phenylalanine (Phe) can activate GPR139 in the hundred-micromolar range. In this study, we produced a mutagenesis-guided model of the GPR139 binding site to form a foundation for future structure-based ligand optimization. Receptor mutants studied in a Ca2+ assay demonstrated that residues F1093×33, H1875×43, W2416×48 and N2717×38, but not E1083×32, are highly important for the activation of GPR139 as predicted by the receptor model. The initial ligand-receptor complex was optimized through free energy perturbation simulations, generating a refined GPR139 model in agreement with experimental data. In summary, the GPR139 reference surrogate agonists 1a and 7c, and the endogenous amino acids L-Trp and L-Phe share a common binding site, as demonstrated by mutagenesis, ligand docking and free energy calculations.

Probing the Time Dependency of Cyclooxygenase-1 Inhibitors by Computer Simulations

Time-dependent inhibition of the cyclooxygenases (COX) by a range of nonsteroidal anti-inflammatory drugs has been described since the first experimental assays of COX were performed. Slow tight-binding inhibitors of COX-1 bind in a two-step mechanism in which the EI → EI* transition is slow and practically irreversible. Since then, various properties of the inhibitors have been proposed to cause or affect the time dependency. Conformational changes in the enzyme have also been proposed to cause the time dependency, but no particular structural feature has been identified. Here, we investigated a series of inhibitors of COX-1 that are either time-independent or time-dependent using a combination of molecular dynamics simulations, binding free energy calculations, and potential of mean force calculations. We find that the time-dependent inhibitors stabilize a conformational change in the enzyme mainly identified by the rotation of a leucine side chain adjacent to the binding pocket. The induced conformation has been previously shown to be essential for the high binding affinities of tight-binding inhibitors in COX-1. The results of this work show that the structural features of the enzyme involved in both time-dependent and tight-binding inhibition are identical and further identify a structural mechanism responsible for the transition between the two enzyme-inhibitor complexes characteristic of slow tight-binding COX-1 inhibitors.

Structure-Based Rational Design of Adenosine Receptor Ligands

The family of adenosine receptors (ARs) is focus of several medicinal chemistry programs aimed to find new potent and selective drugs. Each receptor subtype has been proposed as a relevant drug target in the treatment of, e.g., cardiovascular or inflammatory diseases, asthma or Parkinson's disease. Until recently, most of these efforts have been dominated by ligand-based or empirical approaches. However, the latest advances in G protein-coupled receptor (GPCR) crystallography allowed for a thorough structural characterization of the A2AAR subtype, which has been crystalized with a number of agonists and antagonists. Consequently, the ligand discovery of AR ligands has been enriched with a number of structure-based approaches. These include the generation of higher-confident homology models for the remaining AR subtypes, virtual screening identification of novel chemotypes, structure-based lead-optimization programs, rationalization of selectivity profiles, or the structural characterization of novel binding sites that enable the design of novel allosteric modulators. Computational methodologies have importantly contributed to the success of these structure-based approaches, and the recent advances in the field are also analyzed in this review. We conclude that the design of adenosine receptor ligands has improved dramatically with the consideration of structure- based approaches, which is paving the way to a better understanding of the biology and pharmacological modulation of this relevant family of receptors.

Aryl Sulfonamide Inhibitors of Insulin-Regulated Aminopeptidase Enhance Spine Density in Primary Hippocampal Neuron Cultures

The zinc metallopeptidase insulin regulated aminopeptidase (IRAP), which is highly expressed in the hippocampus and other brain regions associated with cognitive function, has been identified as a high-affinity binding site of the hexapeptide angiotensin IV (Ang IV). This hexapeptide is thought to facilitate learning and memory by binding to the catalytic site of IRAP to inhibit its enzymatic activity. In support of this hypothesis, low molecular weight, nonpeptide specific inhibitors of IRAP have been shown to enhance memory in rodent models. Recently, it was demonstrated that linear and macrocyclic Ang IV-derived peptides can alter the shape and increase the number of dendritic spines in hippocampal cultures, properties associated with enhanced cognitive performance. After screening a library of 10 500 drug-like substances for their ability to inhibit IRAP, we identified a series of low molecular weight aryl sulfonamides, which exhibit no structural similarity to Ang IV, as moderately potent IRAP inhibitors. A structural and biological characterization of three of these aryl sulfonamides was performed. Their binding modes to human IRAP were explored by docking calculations combined with molecular dynamics simulations and binding affinity estimations using the linear interaction energy method. Two alternative binding modes emerged from this analysis, both of which correctly rank the ligands according to their experimental binding affinities for this series of compounds. Finally, we show that two of these drug-like IRAP inhibitors can alter dendritic spine morphology and increase spine density in primary cultures of hippocampal neurons.

Computational exploration of the binding mode of heme-dependent stimulators into the active catalytic domain of soluble guanylate cyclase

Soluble guanylate cyclase (sGC), the main target of nitric oxide (NO), has been proven to have a significant role in coronary artery disease, pulmonary hypertension, erectile dysfunction, and myocardial infarction. One of its agonists, BAY 41-2272 (Riociguat), has been recently approved for treatment of pulmonary arterial hypertension (PHA), while some others are in clinical phases of development. However, the location of the binding sites for the two known types of agonists, heme-dependent stimulators and heme-independent activators, is a matter of debate, particularly for the first group where both a location on the regulatory (H-NOX) and on the catalytic domain have been suggested by different authors. Here, we address its potential location on the catalytic domain, the unique well characterized at the structural level, by an "in silico" approach. Homology models of the catalytic domain of sGC in "inactive" or "active" conformations were constructed using the structure of previously described crystals of the catalytic domains of "inactive" sGCs (2WZ1, 3ET6) and of "active" adenylate cyclase (1CJU). Each model was submitted to six independent molecular dynamics simulations of about 1 μs. Docking of YC-1, a classic heme-dependent stimulator, to all frames of representative trajectories of "inactive" and "active" conformations, followed by calculation of absolute binding free energies with the linear interaction energy (LIE) method, revealed a potential high-affinity binding site on the "active" structure. The site, located between the pseudo-symmetric and the catalytic site just over the loop β2 -β3 , does not overlap with the forskolin binding site on adenylate cyclases. Proteins 2016; 84:1534-1548. © 2016 Wiley Periodicals, Inc.

GPCR-ModSim: A comprehensive web based solution for modeling G-protein coupled receptors

GPCR-ModSim (http://open.gpcr-modsim.org) is a centralized and easy to use service dedicated to the structural modeling of G-protein Coupled Receptors (GPCRs). 3D molecular models can be generated from amino acid sequence by homology-modeling techniques, considering different receptor conformations. GPCR-ModSim includes a membrane insertion and molecular dynamics (MD) equilibration protocol, which can be used to refine the generated model or any GPCR structure uploaded to the server, including if desired non-protein elements such as orthosteric or allosteric ligands, structural waters or ions. We herein revise the main characteristics of GPCR-ModSim and present new functionalities. The templates used for homology modeling have been updated considering the latest structural data, with separate profile structural alignments built for inactive, partially-active and active groups of templates. We have also added the possibility to perform multiple-template homology modeling in a unique and flexible way. Finally, our new MD protocol considers a series of distance restraints derived from a recently identified conserved network of helical contacts, allowing for a smoother refinement of the generated models which is particularly advised when there is low homology to the available templates. GPCR- ModSim has been tested on the GPCR Dock 2013 competition with satisfactory results.

Binding to and Inhibition of Insulin-Regulated Aminopeptidase by Macrocyclic Disulfides Enhances Spine Density

Angiotensin IV (Ang IV) and related peptide analogs, as well as nonpeptide inhibitors of insulin-regulated aminopeptidase (IRAP), have previously been shown to enhance memory and cognition in animal models. Furthermore, the endogenous IRAP substrates oxytocin and vasopressin are known to facilitate learning and memory. In this study, the two recently synthesized 13-membered macrocyclic competitive IRAP inhibitors HA08 and HA09, which were designed to mimic the N terminus of oxytocin and vasopressin, were assessed and compared based on their ability to bind to the IRAP active site, and alter dendritic spine density in rat hippocampal primary cultures. The binding modes of the IRAP inhibitors HA08, HA09, and of Ang IV in either the extended or γ-turn conformation at the C terminus to human IRAP were predicted by docking and molecular dynamics simulations. The binding free energies calculated with the linear interaction energy method, which are in excellent agreement with experimental data and simulations, have been used to explain the differences in activities of the IRAP inhibitors, both of which are structurally very similar, but differ only with regard to one stereogenic center. In addition, we show that HA08, which is 100-fold more potent than the epimer HA09, can enhance dendritic spine number and alter morphology, a process associated with memory facilitation. Therefore, HA08, one of the most potent IRAP inhibitors known today, may serve as a suitable starting point for medicinal chemistry programs aided by MD simulations aimed at discovering more drug-like cognitive enhancers acting via augmenting synaptic plasticity.

Origin of the Enigmatic Stepwise Tight-Binding Inhibition of Cyclooxygenase-1

Nonsteroidal anti-inflammatory drugs (NSAIDs) are widely used for the treatment of pain, fever, inflammation, and some types of cancers. Their mechanism of action is the inhibition of isoforms 1 and 2 of the enzyme cyclooxygenase (COX-1 and COX-2, respectively). However, both nonselective and selective NSAIDs may have side effects that include gastric intestinal bleeding, peptic ulcer formation, kidney problems, and occurrences of myocardial infarction. The search for selective high-affinity COX inhibitors resulted in a number of compounds characterized by a slow, tight-binding inhibition that occurs in a two-step manner. It has been suggested that the final, only very slowly reversible, tight-binding event is the result of conformational changes in the enzyme. However, the nature of these conformational changes has remained elusive. Here we explore the structural determinants of the tight-binding phenomenon in COX-1 with molecular dynamics and free energy simulations. The calculations reveal how different classes of inhibitors affect the equilibrium between two conformational substates of the enzyme in distinctly different ways. The class of tight-binding inhibitors is found to exclusively stabilize an otherwise unfavorable enzyme conformation and bind significantly stronger to this state than to that normally observed in crystal structures. By also computing free energies of binding to the two enzyme conformations for 16 different NSAIDs, we identify an induced-fit mechanism and the key structural features associated with high-affinity tight binding. These results may facilitate the rational development of new COX inhibitors with improved selectivity profiles.

Structural determinants of subtype selectivity and functional activity of angiotensin II receptors

Agonists of the angiotensin II receptor type 2 (AT2), a G-protein coupled receptor, promote tissue protective effects in cardiovascular and renal diseases, while antagonists reduce neuropathic pain. We here report detailed molecular models that explain the AT2 receptor selectivity of our recent series of non-peptide ligands. In addition, minor structural changes of these ligands that provoke different functional activity are rationalized at a molecular level, and related to the selectivity for the different receptor conformations. These findings should pave the way to structure based drug discovery of AT2 receptor ligands.

Free energy calculations of A2Aadenosine receptor mutation effects on agonist binding

A general computational scheme to evaluate the effects of single point mutations on ligand binding is reported.

Sodium ion binding pocket mutations and adenosine A2A receptor function

Recently we identified a sodium ion binding pocket in a high-resolution structure of the human adenosine A2A receptor. In the present study we explored this binding site through site-directed mutagenesis and molecular dynamics simulations. Amino acids in the pocket were mutated to alanine, and their influence on agonist and antagonist affinity, allosterism by sodium ions and amilorides, and receptor functionality was explored. Mutation of the polar residues in the Na(+) pocket were shown to either abrogate (D52A(2.50) and N284A(7.49)) or reduce (S91A(3.39), W246A(6.48), and N280A(7.45)) the negative allosteric effect of sodium ions on agonist binding. Mutations D52A(2.50) and N284A(7.49) completely abolished receptor signaling, whereas mutations S91A(3.39) and N280A(7.45) elevated basal activity and mutations S91A(3.39), W246A(6.48), and N280A(7.45) decreased agonist-stimulated receptor signaling. In molecular dynamics simulations D52A(2.50) directly affected the mobility of sodium ions, which readily migrated to another pocket formed by Glu13(1.39) and His278(7.43). The D52A(2.50) mutation also decreased the potency of amiloride with respect to ligand displacement but did not change orthosteric ligand affinity. In contrast, W246A(6.48) increased some of the allosteric effects of sodium ions and amiloride, whereas orthosteric ligand binding was decreased. These new findings suggest that the sodium ion in the allosteric binding pocket not only impacts ligand affinity but also plays a vital role in receptor signaling. Because the sodium ion binding pocket is highly conserved in other class A G protein-coupled receptors, our findings may have a general relevance for these receptors and may guide the design of novel synthetic allosteric modulators or bitopic ligands.

Structural and energetic effects of A2A adenosine receptor mutations on agonist and antagonist binding

To predict structural and energetic effects of point mutations on ligand binding is of considerable interest in biochemistry and pharmacology. This is not only useful in connection with site-directed mutagenesis experiments, but could also allow interpretation and prediction of individual responses to drug treatment. For G-protein coupled receptors systematic mutagenesis has provided the major part of functional data as structural information until recently has been very limited. For the pharmacologically important A(2A) adenosine receptor, extensive site-directed mutagenesis data on agonist and antagonist binding is available and crystal structures of both types of complexes have been determined. Here, we employ a computational strategy, based on molecular dynamics free energy simulations, to rationalize and interpret available alanine-scanning experiments for both agonist and antagonist binding to this receptor. These computer simulations show excellent agreement with the experimental data and, most importantly, reveal the molecular details behind the observed effects which are often not immediately evident from the crystal structures. The work further provides a distinct validation of the computational strategy used to assess effects of point-mutations on ligand binding. It also highlights the importance of considering not only protein-ligand interactions but also those mediated by solvent water molecules, in ligand design projects.

Toward an optimal docking and free energy calculation scheme in ligand design with application to COX-1 inhibitors

Cyclooxygenase-1 (COX-1) is one of the main targets of most pain-relieving pharmaceuticals. Although the enzyme is well characterized, it is known to be a difficult target for automated molecular docking and scoring. We collected from the literature a structurally diverse set of 45 nonsteroidal anti-inflammatory drugs (NSAIDs) and COX-2-selective inhibitors (coxibs) with a wide range of binding affinities for COX-1. The binding of this data set to a homology model of human COX-1 was analyzed with different combinations of molecular docking algorithms, scoring functions, and the linear interaction energy (LIE) method for estimating binding affinities. It is found that the computational protocols for estimation of binding affinities are extremely sensitive to the initial orientations of the ligands in the binding pocket. To overcome this limitation, we propose a systematic exploration of docking poses using the LIE calculations as a postscoring function. This scheme yields predictions in excellent agreement with experiment, with a mean unsigned error of 0.9 kcal/mol for binding free energies and structures of high quality. A significant improvement of the results is also seen when averaging over experimental data from several independent measurements.

Computational prediction of alanine scanning and ligand binding energetics in G-protein coupled receptors

Site-directed mutagenesis combined with binding affinity measurements is widely used to probe the nature of ligand interactions with GPCRs. Such experiments, as well as structure-activity relationships for series of ligands, are usually interpreted with computationally derived models of ligand binding modes. However, systematic approaches for accurate calculations of the corresponding binding free energies are still lacking. Here, we report a computational strategy to quantitatively predict the effects of alanine scanning and ligand modifications based on molecular dynamics free energy simulations. A smooth stepwise scheme for free energy perturbation calculations is derived and applied to a series of thirteen alanine mutations of the human neuropeptide Y1 receptor and series of eight analogous antagonists. The robustness and accuracy of the method enables univocal interpretation of existing mutagenesis and binding data. We show how these calculations can be used to validate structural models and demonstrate their ability to discriminate against suboptimal ones.

Discovery of 3,4-Dihydropyrimidin-2(1H)-ones As a Novel Class of Potent and Selective A2B Adenosine Receptor Antagonists

We describe the discovery and optimization of 3,4-dihydropyrimidin-2(1H)-ones as a novel family of (nonxanthine) A2B receptor antagonists that exhibit an unusually high selectivity profile. The Biginelli-based hit optimization process enabled a thoughtful exploration of the structure-activity and structure-selectivity relationships for this chemotype, enabling the identification of ligands that combine structural simplicity with excellent hA2B AdoR affinity and remarkable selectivity profiles.