The GPCR crystallography boom: providing an invaluable source of structural information and expanding the scope of homology modeling

Computational study of the structural ensemble of CC chemokine receptor type 5 (CCR5) and its interactions with different ligands

CC Chemokine receptor 5 (CCR5), a member of the Superfamily of G Protein-Coupled Receptors (GPCRs), is an important effector in multiple physiopathological processes such as inflammatory and infectious entities, including central nervous system neuroinflammatory diseases such as Alzheimer's disease, recovery from nervous injuries, and in the HIV-AIDS infective processes. Thus, CCR5 is an attractive target for pharmacological modulation. Since maraviroc was described as a CCR5 ligand that modifies the HIV-AIDS progression, multiple efforts have been developed to describe the functionality of the receptor. In this work, we characterized key structural features of the CCR5 receptor employing extensive atomistic molecular dynamics (MD) in its apo form and in complex with an endogenous agonist, the chemokine CCL5/RANTES, an HIV entry inhibitor, the partial inverse agonist maraviroc, and the experimental antagonists Compound 21 and 34, aiming to elucidate the structural features and mechanistic processes that constitute its functional states, contributing with structural details and a general understanding of this relevant system.

Thyroid-Stimulating Hormone Receptor: the Role in the Development of Thyroid Pathology and Its Correction

Abstract

One of the key elements responsible for the thyroid response to thyroid-stimulating hormone (TSH) is the TSH receptor (TSHR), which belongs to the G protein-coupled receptor superfamily. Binding of TSH or stimulatory autoantibodies to the TSHR extracellular domain triggers multiple signaling pathways in target cells that are mediated through various types of G proteins and β-arrestins. Inhibitory autoantibodies, in contrast, suppress TSHR activity, inducing hypothyroid states. Activating mutations lead to constitutively active TSHR forms and can trigger cancer. Therefore, the TSHR is one of the key targets for the regulation of thyroid function and thyroid status, as well as correction of diseases caused by changes in TSHR activity (autoimmune hyper- and hypothyroidism, Graves’ ophthalmopathy, thyroid cancer). TSH preparations are extremely rarely used in medicine due to their immunogenicity and severe side effects. Most promising is the development of low-molecular allosteric TSHR regulators with an activity of full and inverse agonists and neutral antagonists, which are able to penetrate into the allosteric site located in the TSHR transmembrane domain and specifically bind to it, thus controlling the ability of the receptor to interact with G proteins and β-arrestins. Allosteric regulators do not affect the binding of TSH and autoantibodies to the receptor, which enables mild and selective regulation of thyroid function, while avoiding critical changes in TSH and thyroid hormone levels. The present review addresses the current state of the problem of regulating TSHR activity, including the possibility of using ligands of its allosteric sites.

Molecular Modeling of Histamine Receptors-Recent Advances in Drug Discovery

The recent developments of fast reliable docking, virtual screening and other algorithms gave rise to discovery of many novel ligands of histamine receptors that could be used for treatment of allergic inflammatory disorders, central nervous system pathologies, pain, cancer and obesity. Furthermore, the pharmacological profiles of ligands clearly indicate that these receptors may be considered as targets not only for selective but also for multi-target drugs that could be used for treatment of complex disorders such as Alzheimer's disease. Therefore, analysis of protein-ligand recognition in the binding site of histamine receptors and also other molecular targets has become a valuable tool in drug design toolkit. This review covers the period 2014-2020 in the field of theoretical investigations of histamine receptors mostly based on molecular modeling as well as the experimental characterization of novel ligands of these receptors.

Capturing Peptide-GPCR Interactions and Their Dynamics

Many biological functions of peptides are mediated through G protein-coupled receptors (GPCRs). Upon ligand binding, GPCRs undergo conformational changes that facilitate the binding and activation of multiple effectors. GPCRs regulate nearly all physiological processes and are a favorite pharmacological target. In particular, drugs are sought after that elicit the recruitment of selected effectors only (biased ligands). Understanding how ligands bind to GPCRs and which conformational changes they induce is a fundamental step toward the development of more efficient and specific drugs. Moreover, it is emerging that the dynamic of the ligand–receptor interaction contributes to the specificity of both ligand recognition and effector recruitment, an aspect that is missing in structural snapshots from crystallography. We describe here biochemical and biophysical techniques to address ligand–receptor interactions in their structural and dynamic aspects, which include mutagenesis, crosslinking, spectroscopic techniques, and mass-spectrometry profiling. With a main focus on peptide receptors, we present methods to unveil the ligand–receptor contact interface and methods that address conformational changes both in the ligand and the GPCR. The presented studies highlight a wide structural heterogeneity among peptide receptors, reveal distinct structural changes occurring during ligand binding and a surprisingly high dynamics of the ligand–GPCR complexes.

Chemistry and Structure-Activity Relationships of Psychedelics

This chapter will summarize structure-activity relationships (SAR) that are known for the classic serotonergic hallucinogens (aka psychedelics), focusing on the three chemical types: tryptamines, ergolines, and phenethylamines. In the brain, the serotonin 5-HT_2A receptor plays a key role in regulation of cortical function and cognition, and also appears to be the principal target for hallucinogenic/psychedelic drugs such as LSD. It is one of the most extensively studied of the 14 known types of serotonin receptors. Important structural features will be identified for activity and, where possible, those that the psychedelics have in common will be discussed. Because activation of the 5-HT_2A receptor is the principal mechanism of action for psychedelics, compounds with 5-HT_2A agonist activity generally are quickly discarded by the pharmaceutical industry. Thus, most of the research on psychedelics can be related to activation of 5-HT_2A receptors. Therefore, much of the discussion will include not only clinical or anecdotal studies, but also will consider data from animal models as well as a certain amount of molecular pharmacology where it is known.

Assessing GPCR homology models constructed from templates of various transmembrane sequence identities: Binding mode prediction and docking enrichment

GPCR crystal structures have become more readily accessible in recent years. However, homology models of GPCRs continue to play an important role as many GPCR structures remain unsolved. The new crystal structures now available provide not only additional templates for homology modelling but also the opportunity to assess the performance of homology models against their respective crystal structures and gain insight into the performance of such models. In this study we have constructed homology models from templates of various transmembrane sequence identities for eight GPCR targets to better understand the relationship between transmembrane sequence identity and model quality. Model quality was assessed relative to the crystal structure in terms of structural accuracy as well as performance in two typical structure-based drug design applications: ligand binding pose prediction and docking enrichment in virtual screening. Crystal structures significantly outperformed homology models in both assessments. Accurate ligand binding pose prediction was possible but difficult to achieve using homology models, even with the use of induced fit docking. In virtual screening using homology models still conferred significant enrichment compared to random selection, with a clear benefit also observed in using models optimized through induced fit docking. Our results indicate that while homology models that are reasonably accurate structurally can be constructed, without significant refinement homology models will be outperformed by crystal structures in ligand binding pose prediction and docking enrichment regardless of the template used, primarily due to the extremely high level of structural accuracy needed for such applications.

The Science Behind G Protein-Coupled Receptors (GPCRs) and Their Accurate Visual Representation in Scientific Research

G Protein-Coupled Receptors (GPCRs) are transmembrane (TM) proteins that span the cell membrane seven times, and contain intracellular and extracellular domains, comprised of connecting loops, as well as terminal extension sequences. GPCRs bind ligands within their transmembrane and/or extracellular domains. Ligand binding elicits conformational changes that initiate downstream intracellular signaling events through arrestins and G proteins. GPCRs play central roles in many physiological processes, from sensory to neurological, cardiovascular, endocrine, and reproductive functions. This paper strives to provide an entry point to current GPCR science, and to identify visual approaches to communicate select aspects of GPCR structure and function with clarity and accuracy. The overall GPCR structure, primary sequence and the implications of sequence for membrane topology, ligand binding and helical rearrangements accompanying activation are considered and discussed in the context of visualization strategies, including two-dimensional topological diagrams, three-dimensional representations, and common errors that arise from these representations.

Structural-Functional Features of the Thyrotropin Receptor: A Class A G-Protein-Coupled Receptor at Work

The thyroid-stimulating hormone receptor (TSHR) is a member of the glycoprotein hormone receptors, a sub-group of class A G-protein-coupled receptors (GPCRs). TSHR and its endogenous ligand thyrotropin (TSH) are of essential importance for growth and function of the thyroid gland and proper function of the TSH/TSHR system is pivotal for production and release of thyroid hormones. This receptor is also important with respect to pathophysiology, such as autoimmune (including ophthalmopathy) or non-autoimmune thyroid dysfunctions and cancer development. Pharmacological interventions directly targeting the TSHR should provide benefits to disease treatment compared to currently available therapies of dysfunctions associated with the TSHR or the thyroid gland. Upon TSHR activation, the molecular events conveying conformational changes from the extra- to the intracellular side of the cell across the membrane comprise reception, conversion, and amplification of the signal. These steps are highly dependent on structural features of this receptor and its intermolecular interaction partners, e.g., TSH, antibodies, small molecules, G-proteins, or arrestin. For better understanding of signal transduction, pathogenic mechanisms such as autoantibody action and mutational modifications or for developing new pharmacological strategies, it is essential to combine available structural data with functional information to generate homology models of the entire receptor. Although so far these insights are fragmental, in the past few decades essential contributions have been made to investigate in-depth the involved determinants, such as by structure determination via X-ray crystallography. This review summarizes available knowledge (as of December 2016) concerning the TSHR protein structure, associated functional aspects, and based on these insights we suggest several receptor complex models. Moreover, distinct TSHR properties will be highlighted in comparison to other class A GPCRs to understand the molecular activation mechanisms of this receptor comprehensively. Finally, limitations of current knowledge and lack of information are discussed highlighting the need for intensified efforts toward TSHR structure elucidation.

Urolinin: The First Linear Peptidic Urotensin-II Receptor Agonist

This study investigated the role of individual U-II amino acid positions and side chain characteristics important for U-IIR activation. A complete permutation library of 209 U-II variants was studied in an activity screen that contained single substitution variants of each position with one of the other 19 proteinogenic amino acids. Receptor activation was measured using a cell-based high-throughput fluorescence calcium mobilization assay. We generated the first complete U-II substitution map for U-II receptor activation, resulting in a detailed view into the structural features required for receptor activation, accompanied by complementary information from receptor modeling and ligand docking studies. On the basis of the systematic SAR study of U-II, we created 33 further short and linear U-II variants from eight to three amino acids in length, including d- and other non-natural amino acids. We identified the first high-potency linear U-II analogues. Urolinin, a linear U-II agonist (nWWK-Tyr(3-NO₂)-Abu), shows low nanomolar potency as well as improved metabolic stability.

Homology modeling of a Class A GPCR in the inactive conformation: A quantitative analysis of the correlation between model/template sequence identity and model accuracy

With the present work we quantitatively studied the modellability of the inactive state of Class A G protein-coupled receptors (GPCRs). Specifically, we constructed models of one of the Class A GPCRs for which structures solved in the inactive state are available, namely the β₂ AR, using as templates each of the other class members for which structures solved in the inactive state are also available. Our results showed a detectable linear correlation between model accuracy and model/template sequence identity. This suggests that the likely accuracy of the homology models that can be built for a given receptor can be generally forecasted on the basis of the available templates. We also probed whether sequence alignments that allow for the presence of gaps within the transmembrane domains to account for structural irregularities afford better models than the classical alignment procedures that do not allow for the presence of gaps within such domains. As our results indicated, although the overall differences are very subtle, the inclusion of internal gaps within the transmembrane domains has a noticeable a beneficial effect on the local structural accuracy of the domain in question.

The Dynamic Process of Drug-GPCR Binding at Either Orthosteric or Allosteric Sites Evaluated by Metadynamics

Major advances in G Protein-Coupled Receptor (GPCR) structural biology over the past few years have yielded a significant number of high-resolution crystal structures for several different receptor subtypes. This dramatic increase in GPCR structural information has underscored the use of automated docking algorithms for the discovery of novel ligands that can eventually be developed into improved therapeutics. However, these algorithms are often unable to discriminate between different, yet energetically similar, poses because of their relatively simple scoring functions. Here, we describe a metadynamics-based approach to study the dynamic process of ligand binding to/unbinding from GPCRs with a higher level of accuracy and yet satisfying efficiency.

Homology modeling, docking, and molecular dynamics simulation of the receptor GALR2 and its interactions with galanin and a positive allosteric modulator

Galanin receptor type 2 (GALR2) is a class A G-protein-coupled receptor (GPCR), and it has been reported that orthosteric ligands and positive allosteric modulators (PAMs) of GALR2 could potentially be used to treat epilepsy. So far, the X-ray structure of this receptor has not been resolved, and knowledge of the 3D structure of GALR2 may prove informative in attempts to design novel ligands and to explore the mechanism for the allosteric modulation of this receptor. In this study, homology modeling was used to obtain several GALR2 models using known templates. ProSA-web Z-scores and Ramachandran plots as well as pre-screening against a test dataset of known compounds were all utilized to select the best model of GALR2. Molecular dockings of galanin (a peptide) and a nonpeptide ligand were carried out to choose the (GALR2 model)-galanin complex that showed the closest agreement with the corresponding experimental data. Finally, a 50-ns MD simulation was performed to study the interactions between the GALR2 model and the synthetic and endogenous ligands. The results from docking and MD simulation showed that, besides the reported residues, Tyr160(4.60), Ile105(3.32), Ala274(7.35), and Tyr163(ECL2) also appear to play important roles in the binding of galanin. The potential allosteric binding pockets in the GALR2 model were then investigated via MD simulation. The results indicated that the mechanism for the allosteric modulation caused by PAMs is the binding of the PAM at pocket III, which is formed by galanin, ECL2, TM2, TM3, and ECL1; this results in the disruption of the Na(+)-binding site and/or the Na(+) ion pathway, leading to GALR2 agonism.

Molecular modelling studies on 2-substituted octahydropyrazinopyridoindoles for histamine H2 receptor antagonism

The human histamine H2 receptor (hH2HR) is a G-protein coupled receptor protein with seven transmembrane (TM)-spanning helices primarily involved in regulation of gastric acid secretion. Antagonists targeting hH2HR are useful in the treatment of hyperacidic conditions such as peptic ulcers, gastresophageal reflux disease and gastrointestinal bleeding. We have previously reported the antagonism of 2-substituted pyrazinopyridoindoles at the human histamine H1 receptor and mode of binding of these compounds at the hH1HR using in silico methods. Interestingly, some of the compounds in the series also showed promising activity towards hH2HR that prompted us to investigate the mode of binding of these compounds at hH2HR. In the absence of the crystal structure of hH2HR a homology model has been constructed using multiple sequence alignment, using the X-ray crystal structures of Turkey β1-adrenergic receptor (tβ1AR), Human histamine H1 receptor (hH1HR), Human β2-adrenergic receptor (hβ2AR) and Human D3 dopamine receptor (hD3R). The important residues for binding were depicted in TMIII, TMV, TMVI and TMVII by the homology modelled hH2HR for 2-substituted pyrazinopyridoindoles. A comparative study for deducing the selectivity regarding the binding towards hH1HR and hH2HR has been carried out, which may be useful in designing of selective hH1HR/hH2HR antagonists in these classes of compounds.

Expanding the horizons of G protein-coupled receptor structure-based ligand discovery and optimization using homology models

We show the key role of structural homology models in GPCR structure-based lead discovery and optimization, highlighting methodological aspects, recent progress and future directions.

On the modularity of the intrinsic flexibility of the µ opioid receptor: a computational study

The µ opioid receptor (µOR), the principal target to control pain, belongs to the G protein-coupled receptors (GPCRs) family, one of the most highlighted protein families due to their importance as therapeutic targets. The conformational flexibility of GPCRs is one of their essential characteristics as they take part in ligand recognition and subsequent activation or inactivation mechanisms. It is assessed that the intrinsic mechanical properties of the µOR, more specifically its particular flexibility behavior, would facilitate the accomplishment of specific biological functions, at least in their first steps, even in the absence of a ligand or any chemical species usually present in its biological environment. The study of the mechanical properties of the µOR would thus bring some indications regarding the highly efficient ability of the µOR to transduce cellular message. We therefore investigate the intrinsic flexibility of the µOR in its apo-form using all-atom Molecular Dynamics simulations at the sub-microsecond time scale. We particularly consider the µOR embedded in a simplified membrane model without specific ions, particular lipids, such as cholesterol moieties, or any other chemical species that could affect the flexibility of the µOR. Our analyses highlighted an important local effect due to the various bendability of the helices resulting in a diversity of shape and volume sizes adopted by the µOR binding site. Such property explains why the µOR can interact with ligands presenting highly diverse structural geometry. By investigating the topology of the µOR binding site, a conformational global effect is depicted: the correlation between the motional modes of the extra- and intracellular parts of µOR on one hand, along with a clear rigidity of the central µOR domain on the other hand. Our results show how the modularity of the µOR flexibility is related to its pre-ability to activate and to present a basal activity.

Computational studies to predict or explain G protein coupled receptor polypharmacology

Since G protein-coupled receptors (GPCRs) belong to a very large superfamily of evolutionarily related receptors (>800 members in humans), and due to the rapid progress on their structural biology, they are ideal candidates for polypharmacology studies. Broad screening and bioinformatics/chemoinformatics have been applied to understanding off-target effects of GPCR ligands. It is now feasible to approach the question of GPCR polypharmacology using molecular modeling and the available X-ray GPCR structures. As an example, large and sterically constrained adenosine derivatives (potent adenosine receptor ligands with low conformational freedom and multiple extended substituents) were screened for binding at diverse receptors. Unanticipated off-target interactions, including at biogenic amine receptors, were then modeled using a structure-based approach to provide a consistent understanding of recognition. A conserved Asp in TM3 changed its role from counterion for biogenic amines to characteristic H-bonding to adenosine. The same systematic approach could potentially be applied to many GPCRs or other receptors using other sets of congeneric ligands.

Photomodulation of G protein-coupled adenosine receptors by a novel light-switchable ligand

The adenosinergic system operates through G protein-coupled adenosine receptors, which have become promising therapeutic targets for a wide range of pathological conditions. However, the ubiquity of adenosine receptors and the eventual lack of selectivity of adenosine-based drugs have frequently diminished their therapeutic potential. Accordingly, here we aimed to develop a new generation of light-switchable adenosine receptor ligands that change their intrinsic activity upon irradiation, thus allowing the spatiotemporal control of receptor functioning (i.e., receptor activation/inactivation dependent on location and timing). Therefore, we synthesized an orthosteric, photoisomerizable, and nonselective adenosine receptor agonist, nucleoside derivative MRS5543 containing an aryl diazo linkage on the N⁶ substituent, which in the dark (relaxed isomer) behaved as a full adenosine A₃ receptor (A₃R) and partial adenosine A_2A receptor (A_2AR) agonist. Conversely, upon photoisomerization with blue light (460 nm), it remained a full A₃R agonist but became an A_2AR antagonist. Interestingly, molecular modeling suggested that structural differences encountered within the third extracellular loop of each receptor could modulate the intrinsic, receptor subtype-dependent, activity. Overall, the development of adenosine receptor ligands with photoswitchable activity expands the pharmacological toolbox in support of research and possibly opens new pharmacotherapeutic opportunities.

G protein-coupled receptors: computer-aided ligand discovery and computational structural analyses in the 2010s

G protein-coupled receptors, or GPCRs, are a large superfamily of proteins found on the plasma membrane of cells. They are involved in most physiological and pathophysiological functions and constitute the target of the majority of marketed drugs. Although these receptors have been historically elusive to attempts of structural determination, GPCR crystallography is now in full blossom, opening the way to structure-based drug discovery and enabling homology modeling. This thematic issue of the journal In Silico Pharmacology, which illustrates how the expanding body of structural knowledge is fostering complex computational analyses of the structure-function relationships of the receptors and their interactions with their ligands, stems from the 31st Camerino-Cyprus-Noordwijkerhout Symposium held in Italy, in May 2013, at the University of Camerino. Specifically, it originates from a session of the symposium entitled “Structure-Based Discovery of Ligands of G Protein-Coupled Receptors: Finally a Reality”, and features a mix of research articles and reviews on the application of computational modeling to the analysis of the structure of GPCRs and the interactions of the receptors with their ligands.

Graph analysis of β2 adrenergic receptor structures: a "social network" of GPCR residues

Purpose

G protein-coupled receptors (GPCRs) are a superfamily of membrane proteins of vast pharmaceutical interest. Here, we describe a graph theory-based analysis of the structure of the β₂ adrenergic receptor (β₂ AR), a prototypical GPCR. In particular, we illustrate the network of direct and indirect interactions that link each amino acid residue to any other residue of the receptor.

Methods

Networks of interconnected amino acid residues in proteins are analogous to social networks of interconnected people. Hence, they can be studied through the same analysis tools typically employed to analyze social networks – or networks in general – to reveal patterns of connectivity, influential members, and dynamicity. We focused on the analysis of closeness-centrality, which is a measure of the overall connectivity distance of the member of a network to all other members.

Results

The residues endowed with the highest closeness-centrality are located in the middle of the seven transmembrane domains (TMs). In particular, they are mostly located in the middle of TM2, TM3, TM6 or TM7, while fewer of them are located in the middle of TM1, TM4 or TM5. At the cytosolic end of TM6, the centrality detected for the active structure is markedly lower than that detected for the corresponding residues in the inactive structures. Moreover, several residues acquire centrality when the structures are analyzed in the presence of ligands. Strikingly, there is little overlap between the residues that acquire centrality in the presence of the ligand in the blocker-bound structures and the agonist-bound structures.

Conclusions

Our results reflect the fact that the receptor resembles a bow tie, with a rather tight knot of closely interconnected residues and two ends that fan out in two opposite directions: one toward the extracellular space, which hosts the ligand binding cavity, and one toward the cytosol, which hosts the G protein binding cavity. Moreover, they underscore how interaction network is by the conformational rearrangements concomitant with the activation of the receptor and by the presence of agonists or blockers.

Electronic supplementary material

The online version of this article (doi:10.1186/2193-9616-1-16) contains supplementary material, which is available to authorized users.