Angel Gonzalez

Molecular Basis of Ligand Dissociation in β-Adrenergic Receptors

The important and diverse biological functions of β-adrenergic receptors (βARs) have promoted the search for compounds to stimulate or inhibit their activity. In this regard, unraveling the molecular basis of ligand binding/unbinding events is essential to understand the pharmacological properties of these G protein-coupled receptors. In this study, we use the steered molecular dynamics simulation method to describe, in atomic detail, the unbinding process of two inverse agonists, which have been recently co-crystallized with β1 and β2ARs subtypes, along four different channels. Our results indicate that this type of compounds likely accesses the orthosteric binding site of βARs from the extracellular water environment. Importantly, reconstruction of forces and energies from the simulations of the dissociation process suggests, for the first time, the presence of secondary binding sites located in the extracellular loops 2 and 3 and transmembrane helix 7, where ligands are transiently retained by electrostatic and Van der Waals interactions. Comparison of the residues that form these new transient allosteric binding sites in both βARs subtypes reveals the importance of non-conserved electrostatic interactions as well as conserved aromatic contacts in the early steps of the binding process.

A biased ligand for OXE-R uncouples Gα and Gβγ signaling within a heterotrimer

Differential targeting of heterotrimeric G protein versus β-arrestin signaling are emerging concepts in G protein-coupled receptor (GPCR) research and drug discovery, and biased engagement by GPCR ligands of either β-arrestin or G protein pathways has been disclosed. Herein we report on a new mechanism of ligand bias to titrate the signaling specificity of a cell-surface GPCR. Using a combination of biomolecular and virtual screening, we identified the small-molecule modulator Gue1654, which inhibits Gβγ but not Gα signaling triggered upon activation of Gα(i)-βγ by the chemoattractant receptor OXE-R in both recombinant and human primary cells. Gue1654 does not interfere nonspecifically with signaling directly at or downstream of Gβγ. This hitherto unappreciated mechanism of ligand bias at a GPCR highlights both a new paradigm for functional selectivity and a potentially new strategy to develop pathway-specific therapeutics.

Impact of Helix Irregularities on Sequence Alignment and Homology Modeling of G Protein‐Coupled Receptors

Comparison of the crystal structures of G protein‐coupled receptors (GPCRs) revealed backbone irregularities in the majority of the transmembrane (TM) helices. Among these, wide (π bulge) and tight (310) helical turns on TM2 and TM5 deserve special attention because of their proximity to the ligand binding site. These irregularities are related to residue insertion or deletion (reflected by inclusion of gaps in sequence alignments) accumulated during the evolution of these two helices. These findings have direct implications for the sequence alignments, phylogeny reconstruction, and homology modeling of class A GPCRs.

Genetic variability of porcine circovirus 2 in vaccinating and non-vaccinating commercial farms.

Vaccines against porcine circovirus 2 (PCV2) are now widely used to control the diseases caused by the virus. Although the vaccines protect pigs against the disease, they do not lead to sterilizing immunity and therefore infections with PCV2 continue in farms. It is expected that, due to its high evolutionary rate, PCV2 can adapt quickly to environmental pressures such as vaccination. The goal of this study was to elucidate the molecular variation of PCV2 in relation to vaccination. PCV2 variability was investigated from samples of infected pigs from five farms where vaccination had never been applied and two farms where pigs had been vaccinated for at least 2 years. For the genetic analysis, full PCV2 genomes were amplified and subsequently pooled by vaccination status from serum of eight vaccinated, infected pigs and 16 non-vaccinated, infected pigs. Variability of viral populations was quantified using next-generation sequencing and subsequent bioinformatics analysis. The number of segregating sites was similar in the non-vaccinated (n=109) and vaccinated pools (n=96), but the distribution of these sites in the genome differed. Most notably, in the capsid gene, the number of segregating sites was observed only in the non-vaccinated population. Based on the structural analysis, it is expected that some low-frequency amino acids result in biologically low-fit viruses. On the contrary, D294 in replicase represents a novel amino acid which was dominant and unique in the vaccinated pool. This work showed that variable PCV2 populations were circulating in commercial farms, and that this variability was different in samples obtained from vaccinating and non-vaccinating farms.

The importance of solvation in the design of ligands targeting membrane proteins

A crucial contribution to the ligand-receptor binding affinity is, in addition to their electrostatic and van der Waals interactions, the desolvation of the ligand. This is of special relevance in membrane proteins because the ligand has to be transferred from the aqueous environment to the transmembrane binding site crevice. Herein we report the synthesis of new serotonin 5-HT4receptor antagonists that replace a key carbonyl group by the thiocarbonyl bioisoster. This modification enhances experimental 5-HT4 receptor binding affinities by as much as 91 times. Free energy perturbation calculations have shown that the significant decrease of the penalty of desolvation, facilitating the entrance of the ligands into the binding site crevice, compensates for the weaker ligand-receptor interaction.

The role of hydrophobic amino acids in the structure and function of the rhodopsin family of G protein-coupled receptors.

Recent advances in crystallization methods have permitted to resolve the molecular structure of several members of the rhodopsin family of G protein-coupled receptors (GPCRs). Comparison among these structures revealed a number of conserved polar and charged residues implicated in the receptor transduction pathways. These residues function as micro-switches in the process of receptor activation and has been the object of study of many research groups. However, hydrophobic forces, usually underappreciated, also play a major role in GPCR function. Conserved hydrophobic residues contribute significantly to receptor activation, G protein coupling, and oligomerization processes. This review focuses on the impact of the hydrophobic amino acids observed in the structure of class A GPCRs necessary for their function. This information represents a fundamental piece to complete a holistic view of the GPCR signal transduction machinery.

Modeling of G Protein-Coupled Receptors Using Crystal Structures: From Monomers to Signaling Complexes

G protein-coupled receptors constitute a large and functionally diverse family of transmembrane proteins. They are fundamental in the transfer of extracellular stimuli to intracellular signaling pathways and are among the most targeted proteins in drug discovery. Recent advances in crystallization methods have permitted to resolve the molecular structure of several members of the family. This chapter focuses on the impact of these structures in the use of homology modeling techniques for building three-dimensional models of homologous G protein-coupled receptors, higher order oligomers, and their complexes with ligands and signaling proteins.

GPCRtm: An amino acid substitution matrix for the transmembrane region of class A G Protein-Coupled Receptors

BackgroundProtein sequence alignments and database search methods use standard scoring matrices calculated from amino acid substitution frequencies in general sets of proteins. These general-purpose matrices are not optimal to align accurately sequences with marked compositional biases, such as hydrophobic transmembrane regions found in membrane proteins. In this work, an amino acid substitution matrix (GPCRtm) is calculated for the membrane spanning segments of the G protein-coupled receptor (GPCR) rhodopsin family; one of the largest transmembrane protein family in humans with great importance in health and disease.ResultsThe GPCRtm matrix reveals the amino acid compositional bias distinctive of the GPCR rhodopsin family and differs from other standard substitution matrices. These membrane receptors, as expected, are characterized by a high content of hydrophobic residues with regard to globular proteins. On the other hand, the presence of polar and charged residues is higher than in average membrane proteins, displaying high frequencies of replacement within themselves.ConclusionsAnalysis of amino acid frequencies and values obtained from the GPCRtm matrix reveals patterns of residue replacements different from other standard substitution matrices. GPCRs prioritize the reactivity properties of the amino acids over their bulkiness in the transmembrane regions. A distinctive role is that charged and polar residues seem to evolve at different rates than other amino acids. This observation is related to the role of the transmembrane bundle in the binding of ligands, that in many cases involve electrostatic and hydrogen bond interactions. This new matrix can be useful in database search and for the construction of more accurate sequence alignments of GPCRs.

Ligand chain length drives activation of lipid G protein-coupled receptors

Sphingosine-1-phosphate (S1P) is a lipid mediator that can activate five cell membrane G protein-coupled receptors (GPCRs) which carry a variety of essential functions and are promising drug targets. S1P is composed of a polar zwitterionic head-group and a hydrophobic alkyl chain. This implies an activation mechanism of its cognate receptor that must be significantly different from what is known for prototypical GPCRs (ie receptor to small hydrophilic ligands). Here we aim to identify the structural features responsible for S1P agonism by combining molecular dynamics simulations and functional assays using S1P analogs of different alkyl chain lengths. We propose that high affinity binding involves polar interactions between the lipid head-group and receptor side chains while activation is due to hydrophobic interactions between the lipid tail and residues in a distinct binding site. We observe that ligand efficacy is directly related to alkyl chain length but also varies with receptor subtypes in correlation with the size of this binding pocket. Integrating experimental and computational data, we propose an activation mechanism for the S1P receptors involving agonist-induced conformational events that are conserved throughout class A GPCRs.

Development of a nucleotide exchange inhibitor that impairs Ras oncogenic signaling

Despite more than three decades of intense effort, no anti-Ras therapies have reached clinical application. Contributing to this failure has been an underestimation of Ras complexity and a dearth of structural information. In this regard, recent studies have revealed the highly dynamic character of the Ras surface and the existence of transient pockets suitable for small-molecule binding, opening up new possibilities for the development of Ras modulators. Herein, a novel Ras inhibitor (compound 12) is described that selectively impairs mutated Ras activity in a reversible manner without significantly affecting wild-type Ras, reduces the Ras-guanosine triphosphate (GTP) levels, inhibits the activation of the mitogen-activated protein kinase (MAPK) pathway, and exhibits remarkable cytotoxic activity in Ras-driven cellular models. The use of molecular dynamics simulations and NMR spectroscopy experiments has enabled the molecular bases responsible for the interactions between compound 12 and Ras protein to be explored. The new Ras inhibitor binds partially to the GTP-binding region and extends into the adjacent hydrophobic pocket delimited by switch II. Hence, Ras inhibitor 12 could represent a new compound for the development of more efficacious drugs to target Ras-driven cancers; a currently unmet clinical need.