Hugo Gutiérrez-de-Terán

The Role of a Sodium Ion Binding Site in the Allosteric Modulation of the A2A Adenosine G Protein-Coupled Receptor

The function of G protein-coupled receptors (GPCRs) can be modulated by a number of endogenous allosteric molecules. In this study, we used molecular dynamics, radioligand binding, and thermostability experiments to elucidate the role of the recently discovered sodium ion binding site in the allosteric modulation of the human A(2A) adenosine receptor, conserved among class A GPCRs. While the binding of antagonists and sodium ions to the receptor was noncompetitive in nature, the binding of agonists and sodium ions appears to require mutually exclusive conformational states of the receptor. Amiloride analogs can also bind to the sodium binding pocket, showing distinct patterns of agonist and antagonist modulation. These findings suggest that physiological concentrations of sodium ions affect functionally relevant conformational states of GPCRs and can help to design novel synthetic allosteric modulators or bitopic ligands exploiting the sodium ion binding pocket.

Community-wide assessment of GPCR structure modelling and ligand docking

Recent breakthroughs in the determination of the crystal structures of G protein-coupled receptors (GPCRs) have provided new opportunities for structure-based drug design strategies targeting this protein family. With the aim of evaluating the current status of GPCR structure prediction and ligand docking, a community-wide, blind prediction assessment — GPCR Dock 2008 — was conducted in coordination with the publication of the crystal structure of the human adenosine A2A receptor bound to the ligand ZM241385. Twenty-nine groups submitted 206 structural models before the release of the experimental structure, which were evaluated for the accuracy of the ligand binding mode and the overall receptor model compared with the crystal structure. This analysis highlights important aspects for success and future development, such as accurate modelling of structurally divergent regions and use of additional biochemical insight such as disulphide bridges in the extracellular loops.

Crystal structure of thioflavin-T and its binding to amyloid fibrils: insights at the molecular level

Combining X-ray data on thioflavin-T and theoretical calculations on its binding to a peptide model for Abeta(1-42) fibrils gives evidence of main stabilizing interactions, which influence the dihedral angle between the two moieties of thioflavin-T and thereby its fluorescence properties; these results shed new light on possible strategies for the design of dyes to bind amyloid fibrils more specifically.

Pyrimidine Derivatives as Potent and Selective A3 Adenosine Receptor Antagonists

Two regioisomeric series of diaryl 2- or 4-amidopyrimidines have been synthesized and their adenosine receptor affinities were determined in radioligand binding assays at the four human adenosine receptors (hARs). Some of the ligands prepared herein exhibit remarkable affinities (K(i) < 10 nm) and, most noticeably, the absence of activity at the A(1), A(2A), and A(2B) receptors. The structural determinants that support the affinity and selectivity profiles of the series were highlighted through an integrated computational approach, combining a 3D-QSAR model built on the second generation of GRid INdependent Descriptors (GRIND2) with a novel homology model of the hA(3) receptor. The robustness of the computational model was subsequently evaluated by the design of new derivatives exploring the alkyl substituent of the exocyclic amide group. The synthesis and evaluation of the novel compounds validated the predictive power of the model, exhibiting excellent agreement between predicted and experimental activities.

Computational inhibitor design against malaria plasmepsins

Abstract.Plasmepsins are aspartic proteases involved in the degradation of the host cell hemoglobin that is used as a food source by the malaria parasite. Plasmepsins are highly promising as drug targets, especially when combined with the inhibition of falcipains that are also involved in hemoglobin catabolism. In this review, we discuss the mechanism of plasmepsins I–IV in view of the interest in transition state mimetics as potential compounds for lead development. Inhibitor development against plasmepsin II as well as relevant crystal structures are summarized in order to give an overview of the field. Application of computational techniques, especially binding affinity prediction by the linear interaction energy method, in the development of malarial plasmepsin inhibitors has been highly successful and is discussed in detail. Homology modeling and molecular docking have been useful in the current inhibitor design project, and the combination of such methods with binding free energy calculations is analyzed.

Molecular Modelling of G Protein-Coupled Receptors Through the Web

With the recent crystallization of several G Protein‐Coupled receptors (GPCRs), homology modelling and all atom molecular dynamics (MD) simulations have proven their usefulness for exploring the structure and function of this superfamily of membrane receptors. Subsequently, automated computational protocols have been implemented as web‐based servers in the recent years to produce reliable models of GPCRs, providing partial or global solutions for the structural characterization and molecular simulation of GPCRs. These dedicated modelling services represent an attractive tool for the broader community of public researchers and pharmaceutical companies, in order to assist in the structure‐based drug design of GPCRs. We here collect and analyze the existing web servers, among which a previously unreported service, GPCR‐ModSim, offers for the first time full atom MD simulations in the pipeline for GPCR molecular modelling.

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.

Characterization of the dynamic events of GPCRs by automated computational simulations

The recent advances in membrane protein crystallography have provided extremely valuable structural information of the superfamily of GPCRs (G-protein-coupled receptors). This has been particularly true for a few receptors whose structure was solved several times under different biochemical conditions. It follows that the mechanisms of receptor conformational equilibrium and related dynamic events can be explored by computational simulations. In the present article, we summarize our recent understanding of several dynamic features of GPCRs, accomplished through the use of MD (molecular dynamics) simulations. Our pipeline for the MD simulations of GPCRs, implemented in the web service http://gpcr.usc.es, is updated in the present paper and illustrated by recent applications. Special emphasis is put on the A2A adenosine receptor, one of the selected cases where crystal structures in several conformations and conditions exist, and on the dimerization process of the CXCR4 (CXC chemokine receptor 4).

In silico directed chemical probing of the adenosine receptor family

One of the grand challenges in chemical biology is identifying a small-molecule modulator for each individual function of all human proteins. Instead of targeting one protein at a time, an efficient approach to address this challenge is to target entire protein families by taking advantage of the relatively high levels of chemical promiscuity observed within certain boundaries of sequence phylogeny. We recently developed a computational approach to identifying the potential protein targets of compounds based on their similarity to known bioactive molecules for almost 700 targets. Here, we describe the direct identification of novel antagonists for all four adenosine receptor subtypes by applying our virtual profiling approach to a unique synthesis-driven chemical collection composed of 482 biologically-orphan molecules. These results illustrate the potential role of in silico target profiling to guide efficiently screening campaigns directed to discover new chemical probes for all members of a protein family.

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 (D52A2.50 and N284A7.49) or reduce (S91A3.39, W246A6.48, and N280A7.45) the negative allosteric effect of sodium ions on agonist binding. Mutations D52A2.50 and N284A7.49 completely abolished receptor signaling, whereas mutations S91A3.39 and N280A7.45 elevated basal activity and mutations S91A3.39, W246A6.48, and N280A7.45 decreased agonist-stimulated receptor signaling. In molecular dynamics simulations D52A2.50 directly affected the mobility of sodium ions, which readily migrated to another pocket formed by Glu131.39 and His2787.43. The D52A2.50 mutation also decreased the potency of amiloride with respect to ligand displacement but did not change orthosteric ligand affinity. In contrast, W246A6.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.