Related contribution of specific helix 2 and 7 residues to conformational activation of the serotonin 5-HT2A receptor

A career in numbers: A citation network analysis of the work of RP Millar and his contribution to GnRH research

Here, we reflect on the long career in neuroendocrinology of a single, highly productive scientist ('Bob' Millar), by analysing his oeuvre of published papers through the lens of citation metrics. We use citation network analysis in a novel manner to identify the specific topics to which his papers have made a particular contribution, allowing us to compare the citations of his papers with those of contemporary papers on the same topic, rather than on the same broad field as generally used to normalise citations. It appears that citation rates are highest for topics on which Bob has published a relatively large number of papers that have become core to a tightly-knit community of authors that cite each other. This analysis shows that an author's impact depends on the existence of a receptive community that is alert to the potential utility of papers from that author, and which uses, amplifies, extends and qualifies the contents of their papers-activities that entail reciprocal citation between authors. The obvious conclusion is that a scientist's impact depends on the use that his or her contemporaries make of his or her contributions, rather than on the contributions in themselves.

Therapeutic potentials of nonpeptidic V2R agonists for partial cNDI-causing V2R mutants

Loss-of-function mutations in the type 2 vasopressin receptor (V2R) are a major cause of congenital nephrogenic diabetes insipidus (cNDI). In the context of partial cNDI, the response to desmopressin (dDAVP) is partially, but not entirely, diminished. For those with the partial cNDI, restoration of V2R function would offer a prospective therapeutic approach. In this study, we revealed that OPC-51803 (OPC5) and its structurally related V2R agonists could functionally restore V2R mutants causing partial cNDI by inducing prolonged signal activation. The OPC5-related agonists exhibited functional selectivity by inducing signaling through the Gs-cAMP pathway while not recruiting β-arrestin1/2. We found that six cNDI-related V2R partial mutants (V882.53M, Y1283.41S, L1614.47P, T2736.37M, S3298.47R and S3338.51del) displayed varying degrees of plasma membrane expression levels and exhibited moderately impaired signaling function. Several OPC5-related agonists induced higher cAMP responses than AVP at V2R mutants after prolonged agonist stimulation, suggesting their potential effectiveness in compensating impaired V2R-mediated function. Furthermore, docking analysis revealed that the differential interaction of agonists with L3127.40 caused altered coordination of TM7, potentially contributing to the functional selectivity of signaling. These findings suggest that nonpeptide V2R agonists could hold promise as potential drug candidates for addressing partial cNDI.

Mechanisms and molecular targets surrounding the potential therapeutic effects of psychedelics

Psychedelics, also known as classical hallucinogens, have been investigated for decades due to their potential therapeutic effects in the treatment of neuropsychiatric and substance use disorders. The results from clinical trials have shown promise for the use of psychedelics to alleviate symptoms of depression and anxiety, as well as to promote substantial decreases in the use of nicotine and alcohol. While these studies provide compelling evidence for the powerful subjective experience and prolonged therapeutic adaptations, the underlying molecular reasons for these robust and clinically meaningful improvements are still poorly understood. Preclinical studies assessing the targets and circuitry of the post-acute effects of classical psychedelics are ongoing. Current literature is split between a serotonin 5-HT2A receptor (5-HT2AR)-dependent or -independent signaling pathway, as researchers are attempting to harness the mechanisms behind the sustained post-acute therapeutically relevant effects. A combination of molecular, behavioral, and genetic techniques in neuropharmacology has begun to show promise for elucidating these mechanisms. As the field progresses, increasing evidence points towards the importance of the subjective experience induced by psychedelic-assisted therapy, but without further cross validation between clinical and preclinical research, the why behind the experience and its translational validity may be lost.

All-Atom Molecular Dynamics Simulations Indicated the Involvement of a Conserved Polar Signaling Channel in the Activation Mechanism of the Type I Cannabinoid Receptor

The type I cannabinoid G protein-coupled receptor (CB1, GPCR) is an intensely investigated pharmacological target, owing to its involvement in numerous physiological functions as well as pathological processes such as cancers, neurodegenerative diseases, metabolic disorders and neuropathic pain. In order to develop modern medications that exert their effects through binding to the CB1 receptor, it is essential to understand the structural mechanism of activation of this protein. The pool of atomic resolution experimental structures of GPCRs has been expanding rapidly in the past decade, providing invaluable information about the function of these receptors. According to the current state of the art, the activity of GPCRs involves structurally distinct, dynamically interconverting functional states and the activation is controlled by a cascade of interconnecting conformational switches in the transmembrane domain. A current challenge is to uncover how different functional states are activated and what specific ligand properties are responsible for the selectivity towards those different functional states. Our recent studies of the μ-opioid and β₂-adrenergic receptors (MOP and β₂AR, respectively) revealed that the orthosteric binding pockets and the intracellular surfaces of these receptors are connected through a channel of highly conserved polar amino acids whose dynamic motions are in high correlation in the agonist- and G protein-bound active states. This and independent literature data led us to hypothesize that, in addition to consecutive conformational transitions, a shift of macroscopic polarization takes place in the transmembrane domain, which is furnished by the rearrangement of polar species through their concerted movements. Here, we examined the CB1 receptor signaling complexes utilizing microsecond scale, all-atom molecular dynamics (MD) simulations in order to see if our previous assumptions could be applied to the CB1 receptor too. Apart from the identification of the previously proposed general features of the activation mechanism, several specific properties of the CB1 have been indicated that could possibly be associated with the signaling profile of this receptor.

Evolutionary information helps understand distinctive features of the angiotensin II receptors AT1 and AT2 in amniota

In vertebrates, the octopeptide angiotensin II (AngII) is an important in vivo regulator of the cardiovascular system. It acts mainly through two G protein-coupled receptors, AT1 and AT2. To better understand distinctive features of these receptors, we carried out a phylogenetic analysis that revealed a mirror evolution of AT1 and AT2, each one split into two clades, separating fish from terrestrial receptors. It also revealed that hallmark mutations occurred at, or near, the sodium binding site in both AT1 and AT2. Electrostatics computations and molecular dynamics simulations support maintained sodium binding to human AT1 with slow ingress from the extracellular side and an electrostatic component of the binding free energy around -3kT, to be compared to around -2kT for human AT2 and the δ opioid receptor. Comparison of the sodium binding modes in wild type and mutated AT1 and AT2 from humans and eels indicates that the allosteric control by sodium in both AT1 and AT2 evolved during the transition from fish to amniota. The unusual S7.46N mutation in AT1 is mirrored by a L3.36M mutation in AT2. In the presence of sodium, the N7.46 pattern in amniota AT1 stabilizes the inward orientation of N3.35 in the apo receptor, which should contribute to efficient N3.35 driven biased signaling. The M3.36 pattern in amniota AT2 favours the outward orientation of N3.35 and the receptor promiscuity. Both mutations have physiological consequences for the regulation of the renin-angiotensin system.

The analysis of protein sequences from different species can reveal interesting trends in the structural and functional evolution of a protein family. Here, we analyze the evolution of two G protein-coupled receptors, AT1 and AT2, which bind the angiotensin II peptide and are important regulators of the cardiovascular system. We show that these receptors underwent a mirror evolution. Specific mutations at, or near, the sodium binding pocket occurred in both AT1 and AT2 during the transition to terrestrial life. We carried out electrostatics computations and molecular dynamics simulations to decipher the details of the sodium binding mode in eel and human receptors, as prototypes of fish and amniota receptors. Our results indicate that sodium binding is kinetically slow but thermodynamically stable. Comparison of the sodium binding modes in eel and human receptors reveals that an unusual mutation in the sodium binding pocket of AT1 is critical for biased signaling of amniota AT1 whereas a mutation in AT2 promotes promiscuity of amniota AT2. In turn, these data indicate that a few mutations at a strategic position (here the sodium binding pocket) are an efficient way to gain functional evolution.

Correlated Motions of Conserved Polar Motifs Lay out a Plausible Mechanism of G Protein-Coupled Receptor Activation

Recent advancements in the field of experimental structural biology have provided high-resolution structures of active and inactive state G protein-coupled receptors (GPCRs), a highly important pharmaceutical target family, but the process of transition between these states is poorly understood. According to the current theory, GPCRs exist in structurally distinct, dynamically interconverting functional states of which populations are shifted upon binding of ligands and intracellular signaling proteins. However, explanation of the activation mechanism, on an entirely structural basis, gets complicated when multiple activation pathways and active receptor states are considered. Our unbiased, atomistic molecular dynamics simulations of the μ opioid receptor (MOP) revealed that transmission of external stimulus to the intracellular surface of the receptor is accompanied by subtle, concerted movements of highly conserved polar amino acid side chains along the 7th transmembrane helix. This may entail the rearrangement of polar species and the shift of macroscopic polarization in the transmembrane domain, triggered by agonist binding. Based on our observations and numerous independent indications, we suggest amending the widely accepted theory that the initiation event of GPCR activation is the shift of macroscopic polarization between the ortho- and allosteric binding pockets and the intracellular G protein-binding interface.

Dissecting Role of Charged Residue from Transmembrane Domain 5 of Latent Membrane Protein 1 via In Silico Simulations and Wet-Lab Experiments

Charged residues are frequently found in the transmembrane segments of membrane proteins, which reside in the hydrophobic bilayer environment. Charged residues are critical for the function of membrane protein. However, studies of their role in protein oligomerization are limited. By taking the fifth transmembrane domain (TMD5) of latent membrane protein 1 from the Epstein-Barr virus as a prototype model, in silico simulations and wet-lab experiments were performed to investigate how the charged states affect transmembrane domain oligomerization. Molecular dynamics (MD) simulations showed that the D150-protonated TMD5 trimer was stable, whereas unprotonated D150 created bends in the helices which distort the trimeric structure. D150 was mutated to asparagine to mimic the protonated D150 in TMD5, and the MD simulations of different D150N TMD5 trimers supported that the protonation state of D150 was critical for the trimerization of TMD5. In silico mutations found that D150N TMD5 preferred to interact with TMD5 to form the heterotrimer (1 D150N TMD5:2 protonated TMD5s) rather than the heterotrimer (2 D150N TMD5s:1 protonated TMD5). D150R TMD5 interacted with TMD5 to form the heterotrimer (1 D150R TMD5:2 protonated TMD5). These in silico results imply that D150N TMD5 and D150R TMD5 peptides may be probes for disrupting TMD5 trimerization, which was supported by the dominant-negative ToxR assay in bacterial membranes. In all, this study elucidates the role of charged residues at the membrane milieu in membrane protein oligomerization and provides insight into the development of oligomerization-regulating peptides for modulating transmembrane domain lateral interactions.

Overview of the major classes of new psychoactive substances, psychoactive effects, analytical determination and conformational analysis of selected illegal drugs

The misuse of psychoactive substances is attracting a great deal of attention from the general public. An increase use of psychoactive substances is observed among young people who do not have enough awareness of the harmful effects of these substances. Easy access to illicit drugs at low cost and lack of effective means of routine screening for new psychoactive substances (NPS) have contributed to the rapid increase in their use. New research and evidence suggest that drug use can cause a variety of adverse psychological and physiological effects on human health (anxiety, panic, paranoia, psychosis, and seizures). We describe different classes of these NPS drugs with emphasis on the methods used to identify them and the identification of their metabolites in biological specimens. This is the first review that thoroughly gives the literature on both natural and synthetic illegal drugs with old known data and very hot new topics and investigations, which enables the researcher to use it as a starting point in the literature exploration and planning of the own research. For the first time, the conformational analysis was done for selected illegal drugs, giving rise to the search of the biologically active conformations both theoretically and using lab experiments.

International Union of Basic and Clinical Pharmacology. CX. Classification of Receptors for 5-hydroxytryptamine; Pharmacology and Function

5-HT receptors expressed throughout the human body are targets for established therapeutics and various drugs in development. Their diversity of structure and function reflects the important role 5-HT receptors play in physiologic and pathophysiological processes. The present review offers a framework for the official receptor nomenclature and a detailed understanding of each of the 14 5-HT receptor subtypes, their roles in the systems of the body, and, where appropriate, the (potential) utility of therapeutics targeting these receptors. SIGNIFICANCE STATEMENT: This review provides a comprehensive account of the classification and function of 5-hydroxytryptamine receptors, including how they are targeted for therapeutic benefit.

Homology Modeling of Class A G-Protein-Coupled Receptors in the Age of the Structure Boom

With 700 members, G protein-coupled receptors (GPCRs) of the rhodopsin family (class A) form the largest membrane receptor family in humans and are the target of about 30% of presently available pharmaceutical drugs. The recent boom in GPCR structures led to the structural resolution of 57 unique receptors in different states (39 receptors in inactive state only, 2 receptors in active state only and 16 receptors in different activation states). In spite of these tremendous advances, most computational studies on GPCRs, including molecular dynamics simulations, virtual screening and drug design, rely on GPCR models obtained by homology modeling. In this protocol, we detail the different steps of homology modeling with the MODELLER software, from template selection to model evaluation. The present structure boom provides closely related templates for most receptors. If, in these templates, some of the loops are not resolved, in most cases, the numerous available structures enable to find loop templates with similar length for equivalent loops. However, simultaneously, the large number of putative templates leads to model ambiguities that may require additional information based on multiple sequence alignments or molecular dynamics simulations to be resolved. Using the modeling of the human bradykinin receptor B1 as a case study, we show how several templates are managed by MODELLER, and how the choice of template(s) and of template fragments can improve the quality of the models. We also give examples of how additional information and tools help the user to resolve ambiguities in GPCR modeling.

Deciphering collaborative sidechain motions in proteins during molecular dynamics simulations

The dynamic structure of proteins is essential for their functions and may include large conformational transitions which can be studied by molecular dynamics (MD) simulations. However, details of these transitions are difficult to automatically track. To facilitate their analysis, we developed two scores of correlation between sidechain dihedral angles. The CIRCULAR and OMES scores are computed from, respectively, dihedral angle values and rotamer distributions. As a case study, we applied our methods to an activation-like transition of the chemokine receptor CXCR4, observed during accelerated MD simulations. The principal component analysis of the correlation matrices was consistent with the networking structure of the top ranking pairs. Both scores identify a set of residues whose “collaborative” sidechain rotamerization immediately preceded or accompanied the conformational transition of CXCR4. Detailed analysis of the sequential order of these rotamerizations suggests that an allosteric mechanism, involving the outward motion of an asparagine residue in transmembrane helix 3, might be a prerequisite to the large scale conformational transition of CXCR4. This case study provides the proof-of-concept that the correlation methods developed here are valuable exploratory techniques to help decipher complex reactional pathways.

Two highly related odorant receptors specifically detect α-bile acid pheromones in sea lamprey (Petromyzon marinus)

Pheromones play critical roles in habitat identification and reproductive behavior synchronization in the sea lamprey (Petromyzon marinus). The bile acid 3-keto petromyzonol sulfate (3kPZS) is a major component of the sex pheromone mixture from male sea lamprey that induces specific olfactory and behavioral responses in conspecific individuals. Olfactory receptors interact directly with pheromones, which is the first step in their detection, but identifying the cognate receptors of specific pheromones is often challenging. Here, we deorphanized two highly related odorant receptors (ORs), OR320a and OR320b, of P. marinus that respond to 3kPZS. In a heterologous expression system coupled to a cAMP-responsive CRE-luciferase, OR320a and OR320b specifically responded to C24 5α-bile acids, and both receptors were activated by the same set of 3kPZS analogs. OR320a displayed larger responses to all 3kPZS analogs than did OR320b. This difference appeared to be largely determined by a single amino acid residue, Cys-792.56, the C-terminal sixth residue relative to the most conserved residue in the second transmembrane domain (2.56) of OR320a. This region of TM2 residues 2.56-2.60 apparently is critical for the detection of steroid compounds by odorant receptors in lamprey, zebrafish, and humans. Finally, we identified OR320 orthologs in Japanese lamprey (Lethenteron camtschaticum), suggesting that the OR320 family may be widely present in lamprey species and that OR320 may be under purifying selection. Our results provide a system to examine the origin of olfactory steroid detection in vertebrates and to define a highly conserved molecular mechanism for steroid-ligand detection by G protein-coupled receptors.

5-HT2A receptor-dependent phosphorylation of mGlu2 receptor at Serine 843 promotes mGlu2 receptor-operated Gi/o signaling

The serotonin 5-HT_2A and glutamate mGlu₂ receptors continue to attract particular attention, given their implication in psychosis associated with schizophrenia and the mechanism of action of atypical antipsychotics and a new class of antipsychotics, respectively. A large body of evidence indicates a functional crosstalk between both receptors in the brain, but the underlying mechanisms are not entirely elucidated. Here, we have explored the influence of 5-HT_2A receptor upon the phosphorylation pattern of mGlu₂ receptor in light of the importance of specific phosphorylation events in regulating G protein-coupled receptor signaling and physiological outcomes. Among the five mGlu₂ receptor-phosphorylated residues identified in HEK-293 cells, the phosphorylation of Ser⁸⁴³ was enhanced upon mGlu₂ receptor stimulation by the orthosteric agonist LY379268 only in cells co-expressing the 5-HT_2A receptor. Likewise, administration of LY379268 increased mGlu₂ receptor phosphorylation at Ser⁸⁴³ in prefrontal cortex of wild-type mice but not 5-HT_2A^-/- mice. Exposure of HEK-293 cells co-expressing mGlu₂ and 5-HT_2A receptors to 5-HT also increased Ser⁸⁴³ phosphorylation state to a magnitude similar to that measured in LY379268-treated cells. In both HEK-293 cells and prefrontal cortex, Ser⁸⁴³ phosphorylation elicited by 5-HT_2A receptor stimulation was prevented by the mGlu₂ receptor antagonist LY341495, while the LY379268-induced effect was abolished by the 5-HT_2A receptor antagonist M100907. Mutation of Ser⁸⁴³ into alanine strongly reduced G_i/o signaling elicited by mGlu₂ or 5-HT_2A receptor stimulation in cells co-expressing both receptors. Collectively, these findings identify mGlu₂ receptor phosphorylation at Ser⁸⁴³ as a key molecular event that underlies the functional crosstalk between both receptors.

Harnessing Ion-Binding Sites for GPCR Pharmacology

Endogenous ions play important roles in the function and pharmacology of G-protein coupled receptors (GPCRs). Historically the evidence for ionic modulation of GPCR function dates to 1973 with studies of opioid receptors, where it was demonstrated that physiologic concentrations of sodium allosterically attenuated agonist binding. This Na+-selective effect was distinct from effects of other monovalent and divalent cations, with the latter usually counteracting sodium's negative allosteric modulation of binding. Since then, numerous studies documenting the effects of mono- and divalent ions on GPCR function have been published. While ions can act selectively and nonselectively at many sites in different receptors, the discovery of the conserved sodium ion site in class A GPCR structures in 2012 revealed the unique nature of Na+ site, which has emerged as a near-universal site for allosteric modulation of class A GPCR structure and function. In this review, we synthesize and highlight recent advances in the functional, biophysical, and structural characterization of ions bound to GPCRs. Taken together, these findings provide a molecular understanding of the unique roles of Na+ and other ions as GPCR allosteric modulators. We will also discuss how this knowledge can be applied to the redesign of receptors and ligand probes for desired functional and pharmacological profiles. SIGNIFICANCE STATEMENT: The function and pharmacology of GPCRs strongly depend on the presence of mono and divalent ions in experimental assays and in living organisms. Recent insights into the molecular mechanism of this ion-dependent allosterism from structural, biophysical, biochemical, and computational studies provide quantitative understandings of the pharmacological effects of drugs in vitro and in vivo and open new avenues for the rational design of chemical probes and drug candidates with improved properties.

A Machine Learning Approach for the Discovery of Ligand-Specific Functional Mechanisms of GPCRs

G protein-coupled receptors (GPCRs) play a key role in many cellular signaling mechanisms, and must select among multiple coupling possibilities in a ligand-specific manner in order to carry out a myriad of functions in diverse cellular contexts. Much has been learned about the molecular mechanisms of ligand-GPCR complexes from Molecular Dynamics (MD) simulations. However, to explore ligand-specific differences in the response of a GPCR to diverse ligands, as is required to understand ligand bias and functional selectivity, necessitates creating very large amounts of data from the needed large-scale simulations. This becomes a Big Data problem for the high dimensionality analysis of the accumulated trajectories. Here we describe a new machine learning (ML) approach to the problem that is based on transforming the analysis of GPCR function-related, ligand-specific differences encoded in the MD simulation trajectories into a representation recognizable by state-of-the-art deep learning object recognition technology. We illustrate this method by applying it to recognize the pharmacological classification of ligands bound to the 5-HT_2A and D2 subtypes of class-A GPCRs from the serotonin and dopamine families. The ML-based approach is shown to perform the classification task with high accuracy, and we identify the molecular determinants of the classifications in the context of GPCR structure and function. This study builds a framework for the efficient computational analysis of MD Big Data collected for the purpose of understanding ligand-specific GPCR activity.

Evolution of chemokine receptors is driven by mutations in the sodium binding site

Chemokines and their receptors (members of the GPCR super-family) are involved in a wide variety of physiological processes and diseases; thus, understanding the specificity of the chemokine receptor family could help develop new receptor specific drugs. Here, we explore the evolutionary mechanisms that led to the emergence of the chemokine receptors. Based on GPCR hierarchical classification, we analyzed nested GPCR sets with an eigen decomposition approach of the sequence covariation matrix and determined three key residues whose mutation was crucial for the emergence of the chemokine receptors and their subsequent divergence into homeostatic and inflammatory receptors. These residues are part of the allosteric sodium binding site. Their structural and functional roles were investigated by molecular dynamics simulations of CXCR4 and CCR5 as prototypes of homeostatic and inflammatory chemokine receptors, respectively. This study indicates that the three mutations crucial for the evolution of the chemokine receptors dramatically altered the sodium binding mode. In CXCR4, the sodium ion is tightly bound by four protein atoms and one water molecule. In CCR5, the sodium ion is mobile within the binding pocket and moves between different sites involving from one to three protein atoms and two to five water molecules. Analysis of chemokine receptor evolution reveals that a highly constrained sodium binding site characterized most ancient receptors, and that the constraints were subsequently loosened during the divergence of this receptor family. We discuss the implications of these findings for the evolution of the chemokine receptor functions and mechanisms of action.

Sodium functions as a negative allosteric modulator of the oxytocin receptor

The oxytocin receptor, a class A G protein coupled receptor (GPCR), is essentially involved in the physiology of reproduction. Two parameters are crucially important to support high-affinity agonist binding of the receptor: Mg²⁺ and cholesterol, both acting as positive modulators. Using displacement assays with a high-affinity fluorescent antagonist (OTAN-A647), we now show that sodium functions as a negative allosteric modulator of the oxytocin receptor. In membranes from HEK293 cells stably expressing the oxytocin receptor, oxytocin binding occurred with about 15-fold lower affinity when sodium chloride was increased from 0 to 300 mM, whereas antagonist binding remained largely unchanged. The effect was concentration-dependent, sodium-specific, and it was also observed for oxytocin receptors endogenously expressed in Hs578T breast cancer cells. A conserved Asp (Asp 85) is known to stabilize the sodium binding site in other GCPRs. Mutations of this residue into Ala or Asn are known to yield non-functional oxytocin receptors. When Asp 85 was exchanged for Glu, most of the oxytocin receptors were localized in intracellular structures, but a faint plasma membrane labeling with OTAN-A647 and the appearance of oxytocin-induced calcium responses indicated that these receptors were functional. However, a sodium effect was not detectable for the mutant D85E oxytocin receptors. Thus, the oxytocin receptor is allosterically controlled by sodium similar to other GPCRs, but it behaves differently concerning the involvement of the conserved Asp 85. In case of the oxytocin receptor, Asp 85 is obviously essential for proper localization in the plasma membrane.

GPCRs: What Can We Learn from Molecular Dynamics Simulations?

Advances in the structural biology of G-protein Coupled Receptors have resulted in a significant step forward in our understanding of how this important class of drug targets function at the molecular level. However, it has also become apparent that they are very dynamic molecules, and moreover, that the underlying dynamics is crucial in shaping the response to different ligands. Molecular dynamics simulations can provide unique insight into the dynamic properties of GPCRs in a way that is complementary to many experimental approaches. In this chapter, we describe progress in three distinct areas that are particularly difficult to study with other techniques: atomic level investigation of the conformational changes that occur when moving between the various states that GPCRs can exist in, the pathways that ligands adopt during binding/unbinding events and finally, the influence of lipids on the conformational dynamics of GPCRs.

In Vitro Effects of the Endocrine Disruptor p,p'-DDT on Human Follitropin Receptor

BACKGROUND

1-chloro-4-[2,2,2-trichloro-1-(4-chlorophenyl)ethyl]benzene (p,p'-DDT) is a persistent environmental endocrine disruptor (ED). Several studies have shown an association between p,p'-DDT exposure and reproductive abnormalities.

OBJECTIVES

To investigate the putative effects of p,p'-DDT on the human follitropin receptor (FSHR) function.

METHODS AND RESULTS

We used Chinese hamster ovary (CHO) cells stably expressing human FSHR to investigate the impact of p,p'-DDT on FSHR activity and its interaction with the receptor. At a concentration of 5 μM p,p'-DDT increased the maximum response of the FSHR to follitropin by 32 ± 7.45%. However, 5 μM p,p'-DDT decreased the basal activity and did not influence the maximal response of the closely related LH/hCG receptor to human chorionic gonadotropin (hCG). The potentiating effect of p,p'-DDT was specific for the FSHR. Moreover, in cells that did not express FSHR, p,p'-DDT had no effect on cAMP response. Thus, the potentiating effect of p,p'-DDT was dependent on the FSHR. In addition, p,p'-DDT increased the sensitivity of FSHR to hCG and to a low molecular weight agonist of the FSHR, 3-((5methyl)-2-(4-benzyloxy-phenyl)-5-{[2-[3-ethoxy-4-methoxy-phenyl)-ethylcarbamoyl]-methyl}-4-oxo-thiazolidin-3-yl)-benzamide (16a). Basal activity in response to p,p'-DDT and potentiation of the FSHR response to FSH by p,p'-DDT varied among FSHR mutants with altered transmembrane domains (TMDs), consistent with an effect of p,p'-DDT via TMD binding. This finding was corroborated by the results of simultaneously docking p,p'-DDT and 16a into the FSHR transmembrane bundle.

CONCLUSION

p,p'-DDT acted as a positive allosteric modulator of the FSHR in our experimental model. These findings suggest that G protein-coupled receptors are additional targets of endocrine disruptors.

CITATION

Munier M, Grouleff J, Gourdin L, Fauchard M, Chantreau V, Henrion D, Coutant R, Schiøtt B, Chabbert M, Rodien P. 2016. In vitro effects of the endocrine disruptor p,p'-DDT on human follitropin receptor. Environ Health Perspect 124:991-999; http://dx.doi.org/10.1289/ehp.1510006.

Activation is Hallucinogenic and Antagonism is Therapeutic: Role of 5-HT2A Receptors in Atypical Antipsychotic Drug Actions

This review summarizes recent studies with 5-hydroxytryptamine2A (5-HT2A) receptors, which represent the major site of action of hallucinogens and a likely site for atypical antipsychotic drug actions. We present evidence demonstrating that atypical antipsychotic drugs, as a group, have a preferentially high affinity for 5-HT2A receptors, compared with their affinities for other neurotransmitter receptors. The 5-HT2A receptor blockade seen with atypical antipsychotic drugs is probably an essential factor in explaining many of the unique features of atypical antipsychotic drugs. Atypical antipsychotic drugs have high affinities for several other 5-HT receptors (5-HT2C, 5-HT6, and 5-HT7), and the potential role of these novel 5-HT receptors in atypical antipsychotic drug action is also summarized.

Molecular Insights into the Transmembrane Domain of the Thyrotropin Receptor

The thyrotropin receptor (TSHR) is a G protein-coupled receptor (GPCR) that is member of the leucine-rich repeat subfamily (LGR). In the absence of crystal structure, the success of rational design of ligands targeting the receptor internal cavity depends on the quality of the TSHR models built. In this subfamily, transmembrane helices (TM) 2 and 5 are characterized by the absence of proline compared to most receptors, raising the question of the structural conformation of these helices. To gain insight into the structural properties of these helices, we carried out bioinformatics and experimental studies. Evolutionary analysis of the LGR family revealed a deletion in TM5 but provided no information on TM2. Wild type residues at positions 2.58, 2.59 or 2.60 in TM2 and/or at position 5.50 in TM5 were substituted to proline. Depending on the position of the proline substitution, different effects were observed on membrane expression, glycosylation, constitutive cAMP activity and responses to thyrotropin. Only proline substitution at position 2.59 maintained complex glycosylation and high membrane expression, supporting occurrence of a bulged TM2. The TSHR transmembrane domain was modeled by homology with the orexin 2 receptor, using a protocol that forced the deletion of one residue in the TM5 bulge of the template. The stability of the model was assessed by molecular dynamics simulations. TM5 straightened during the equilibration phase and was stable for the remainder of the simulations. Our data support a structural model of the TSHR transmembrane domain with a bulged TM2 and a straight TM5 that is specific of glycoprotein hormone receptors.

Biased Gs versus Gq proteins and β-arrestin signaling in the NK1 receptor determined by interactions in the water hydrogen bond network

X-ray structures, molecular dynamics simulations, and mutational analysis have previously indicated that an extended water hydrogen bond network between trans-membranes I-III, VI, and VII constitutes an allosteric interface essential for stabilizing different active and inactive helical constellations during the seven-trans-membrane receptor activation. The neurokinin-1 receptor signals efficiently through Gq, Gs, and β-arrestin when stimulated by substance P, but it lacks any sign of constitutive activity. In the water hydrogen bond network the neurokinin-1 has a unique Glu residue instead of the highly conserved AspII:10 (2.50). Here, we find that this GluII:10 occupies the space of a putative allosteric modulating Na(+) ion and makes direct inter-helical interactions in particular with SerIII:15 (3.39) and AsnVII:16 (7.49) of the NPXXY motif. Mutational changes in the interface between GluII:10 and AsnVII:16 created receptors that selectively signaled through the following: 1) Gq only; 2) β-arrestin only; and 3) Gq and β-arrestin but not through Gs. Interestingly, increased constitutive Gs but not Gq signaling was observed by Ala substitution of four out of the six core polar residues of the network, in particular SerIII:15. Three residues were essential for all three signaling pathways, i.e. the water-gating micro-switch residues TrpVI:13 (6.48) of the CWXP motif and TyrVII:20 (7.53) of the NPXXY motif plus the totally conserved AsnI:18 (1.50) stabilizing the kink in trans-membrane VII. It is concluded that the interface between position II:10 (2.50), III:15 (3.39), and VII:16 (7.49) in the center of the water hydrogen bond network constitutes a focal point for fine-tuning seven trans-membrane receptor conformations activating different signal transduction pathways.

Scanning mutagenesis in a yeast system delineates the role of the NPxxY(x)(5,6)F motif and helix 8 of the adenosine A(2B) receptor in G protein coupling

The adenosine receptor subfamily includes four subtypes: the A1, A2A, A2B and A3 receptors, which all belong to the superfamily of G protein-coupled receptors (GPCRs). The adenosine A2B receptor is the least investigated of the adenosine receptors, and the molecular mechanisms of its activation have hardly been explored. We used a single-GPCR-one-G protein yeast screening method in combination with mutagenesis studies, molecular modeling and bio-informatics to investigate the importance of the different amino acid residues of the NPxxY(x)6F motif and helix 8 in the human adenosine A2B receptor (hA2BR) activation. A scanning mutagenesis protocol was employed, yielding 11 single mutations and one double mutation of the NPxxY(x)6F motif and 16 single mutations of helix 8. The amino acid residues P287(7.50), Y290(7.53), R293(7.56) and I304(8.57) were found to be essential, since mutation of these amino acid residues to alanine led to a complete loss of function. Western blot analysis showed that mutant receptor R293(7.56)A was not expressed, whereas the other proteins were. Amino acid residues that are also important in receptor activation are: N286(7.49), V289(7.52), Y292(7.55), N294(8.47), F297(8.50), R298(8.51), H302(8.55) and R307(8.60). The mutation Y290(7.53)F lost 50% of efficacy, while F297(8.50)A behaved similar to wild type receptor. The double mutation, Y290(7.53)F/F297(8.50)Y, lost around 70% of efficacy and displayed a lower potency for the reference agonist 5'-(N-ethylcarboxamido)adenosine (NECA). This study provides new insight into the molecular interplay and impact of TM7 and helix 8 for hA2B receptor activation, which may be extrapolated to other adenosine receptors and possibly to other GPCRs.

Allosteric sodium in class A GPCR signaling

Despite their functional and structural diversity, G-protein-coupled receptors (GPCRs) share a common mechanism of signal transduction via conformational changes in the seven-transmembrane (7TM) helical domain. New major insights into this mechanism come from the recent crystallographic discoveries of a partially hydrated sodium ion that is specifically bound in the middle of the 7TM bundle of multiple class A GPCRs. This review discusses the remarkable structural conservation and distinct features of the Na(+) pocket in this most populous GPCR class, as well as the conformational collapse of the pocket upon receptor activation. New insights help to explain allosteric effects of sodium on GPCR agonist binding and activation, and sodium's role as a potential co-factor in class A GPCR function.

Novel insights into CB1 cannabinoid receptor signaling: a key interaction identified between the extracellular-3 loop and transmembrane helix 2

Activation of the cannabinoid CB1 receptor (CB1) is modulated by aspartate residue D2.63(176) in transmembrane helix (TMH) 2. Interestingly, D2.63 does not affect the affinity for ligand binding at the CB1 receptor. Studies in class A G protein-coupled receptors have suggested an ionic interaction between residues of TMH2 and 7. In this report, modeling studies identified residue K373 in the extracellular-3 (EC-3) loop in charged interactions with D2.63. We investigated this possibility by performing reciprocal mutations and biochemical studies. D2.63(176)A, K373A, D2.63(176)A-K373A, and the reciprocal mutant with the interacting residues juxtaposed D2.63(176)K-K373D were characterized using radioligand binding and guanosine 5'-3-O-(thio)triphosphate functional assays. None of the mutations resulted in a significant change in the binding affinity of N-(piperidiny-1-yl)-5-(4-chlorophenyl)-1-(2,4-dichloro-phenyl)-4-methyl-1H-pyrazole-3-carboxamide hydrochloride (SR141716A) or (-)-3cis -[2-hydroxyl-4-(1,1-dimethyl-heptyl)phenyl]-trans-4-[3-hydroxyl-propyl] cyclohexan-1-ol (CP55,940). Modeling studies indicated that binding-site interactions and energies of interaction for CP55,940 were similar between wild-type and mutant receptors. However, the signaling of CP55,940, and (R)-(+)-[2,3-dihydro-5-methyl-3-[(4-morpholinyl)methyl]-pyrrolo[1,2,3-de]-1,4-benzoxazin-6-yl](1-naphthalenyl)-methanone mesylate (WIN55,212-2) was impaired at the D2.63(176)A-K373A and the single-alanine mutants. In contrast, the reciprocal D2.63(176)K-K373D mutant regained function for both CP55,940 and WIN55,212-2. Computational results indicate that the D2.63(176)-K373 ionic interaction strongly influences the conformation(s) of the EC-3 loop, providing a structure-based rationale for the importance of the EC-3 loop to signal transduction in CB1. The putative ionic interaction results in the EC-3 loop pulling over the top (extracellular side) of the receptor; this EC-3 loop conformation may serve protective and mechanistic roles. These results suggest that the ionic interaction between D2.63(176) and K373 is important for CB1 signal transduction.

Structural evolution of G-protein-coupled receptors: a sequence space approach

Class A G-protein-coupled receptors (GPCRs) provide a fascinating example of evolutionary success. In this review, we discuss how metric multidimensional scaling (MDS), a multivariate analysis method, complements traditional tree-based phylogenetic methods and helps decipher the mechanisms that drove the evolution of class A GPCRs. MDS provides low-dimensional representations of a distance matrix. Applied to a multiple sequence alignment, MDS represents the sequences in a Euclidean space as points whose interdistances are as close as possible to the distances in the alignment (the so-called sequence space). We detail how to perform the MDS analysis of a multiple sequence alignment and how to analyze and interpret the resulting sequence space. We also show that the projection of supplementary data (a property of the MDS method) can be used to straightforwardly monitor the evolutionary drift of specific subfamilies. The sequence space of class A GPCRs reveals the key role of mutations at the level of the TM2 and TM5 proline residues in the evolution of class A GPCRs.

Action of molecular switches in GPCRs--theoretical and experimental studies

G protein coupled receptors (GPCRs), also called 7TM receptors, form a huge superfamily of membrane proteins that, upon activation by extracellular agonists, pass the signal to the cell interior. Ligands can bind either to extracellular N-terminus and loops (e.g. glutamate receptors) or to the binding site within transmembrane helices (Rhodopsin-like family). They are all activated by agonists although a spontaneous auto-activation of an empty receptor can also be observed. Biochemical and crystallographic methods together with molecular dynamics simulations and other theoretical techniques provided models of the receptor activation based on the action of so-called "molecular switches" buried in the receptor structure. They are changed by agonists but also by inverse agonists evoking an ensemble of activation states leading toward different activation pathways. Switches discovered so far include the ionic lock switch, the 3-7 lock switch, the tyrosine toggle switch linked with the nPxxy motif in TM7, and the transmission switch. The latter one was proposed instead of the tryptophan rotamer toggle switch because no change of the rotamer was observed in structures of activated receptors. The global toggle switch suggested earlier consisting of a vertical rigid motion of TM6, seems also to be implausible based on the recent crystal structures of GPCRs with agonists. Theoretical and experimental methods (crystallography, NMR, specific spectroscopic methods like FRET/BRET but also single-molecule-force-spectroscopy) are currently used to study the effect of ligands on the receptor structure, location of stable structural segments/domains of GPCRs, and to answer the still open question on how ligands are binding: either via ensemble of conformational receptor states or rather via induced fit mechanisms. On the other hand the structural investigations of homoand heterodimers and higher oligomers revealed the mechanism of allosteric signal transmission and receptor activation that could lead to design highly effective and selective allosteric or ago-allosteric drugs.

Understanding functional residues of the cannabinoid CB1

The brain cannabinoid (CB(1)) receptor that mediates numerous physiological processes in response to marijuana and other psychoactive compounds is a G protein coupled receptor (GPCR) and shares common structural features with many rhodopsin class GPCRs. For the rational development of therapeutic agents targeting the CB(1) receptor, understanding of the ligand-specific CB(1) receptor interactions responsible for unique G protein signals is crucial. For a more than a decade, a combination of mutagenesis and computational modeling approaches has been successfully employed to study the ligand-specific CB(1) receptor interactions. In this review, after a brief discussion about recent advances in understanding of some structural and functional features of GPCRs commonly applicable to the CB(1) receptor, the CB(1) receptor functional residues reported from mutational studies are divided into three different types, ligand binding (B), receptor stabilization (S) and receptor activation (A) residues, to delineate the nature of the binding pockets of anandamide, CP55940, WIN55212-2 and SR141716A and to describe the molecular events of the ligand-specific CB(1) receptor activation from ligand binding to G protein signaling. Taken these CB(1) receptor functional residues, some of which are unique to the CB(1) receptor, together with the biophysical knowledge accumulated for the GPCR active state, it is possible to propose the early stages of the CB(1) receptor activation process that not only provide some insights into understanding molecular mechanisms of receptor activation but also are applicable for identifying new therapeutic agents by applying the validated structure-based approaches, such as virtual high throughput screening (HTS) and fragment-based approach (FBA).

Multidimensional scaling reveals the main evolutionary pathways of class A G-protein-coupled receptors

Class A G-protein-coupled receptors (GPCRs) constitute the largest family of transmembrane receptors in the human genome. Understanding the mechanisms which drove the evolution of such a large family would help understand the specificity of each GPCR sub-family with applications to drug design. To gain evolutionary information on class A GPCRs, we explored their sequence space by metric multidimensional scaling analysis (MDS). Three-dimensional mapping of human sequences shows a non-uniform distribution of GPCRs, organized in clusters that lay along four privileged directions. To interpret these directions, we projected supplementary sequences from different species onto the human space used as a reference. With this technique, we can easily monitor the evolutionary drift of several GPCR sub-families from cnidarians to humans. Results support a model of radiative evolution of class A GPCRs from a central node formed by peptide receptors. The privileged directions obtained from the MDS analysis are interpretable in terms of three main evolutionary pathways related to specific sequence determinants. The first pathway was initiated by a deletion in transmembrane helix 2 (TM2) and led to three sub-families by divergent evolution. The second pathway corresponds to the differentiation of the amine receptors. The third pathway corresponds to parallel evolution of several sub-families in relation with a covarion process involving proline residues in TM2 and TM5. As exemplified with GPCRs, the MDS projection technique is an important tool to compare orthologous sequence sets and to help decipher the mutational events that drove the evolution of protein families.

Biased signaling of the angiotensin II type 1 receptor can be mediated through distinct mechanisms

BACKGROUND

Seven transmembrane receptors (7TMRs) can adopt different active conformations facilitating a selective activation of either G protein or β-arrestin-dependent signaling pathways. This represents an opportunity for development of novel therapeutics targeting selective biological effects of a given receptor. Several studies on pathway separation have been performed, many of these on the Angiotensin II type 1 receptor (AT1R). It has been shown that certain ligands or mutations facilitate internalization and/or recruitment of β-arrestins without activation of G proteins. However, the underlying molecular mechanisms remain largely unresolved. For instance, it is unclear whether such selective G protein-uncoupling is caused by a lack of ability to interact with G proteins or rather by an increased ability of the receptor to recruit β-arrestins. Since uncoupling of G proteins by increased ability to recruit β-arrestins could lead to different cellular or in vivo outcomes than lack of ability to interact with G proteins, it is essential to distinguish between these two mechanisms.

METHODOLOGY/PRINCIPAL FINDINGS

We studied five AT1R mutants previously published to display pathway separation: D74N, DRY/AAY, Y292F, N298A, and Y302F (Ballesteros-Weinstein numbering: 2.50, 3.49-3.51, 7.43, 7.49, and 7.53). We find that D74N, DRY/AAY, and N298A mutants are more prone to β-arrestin recruitment than WT. In contrast, receptor mutants Y292F and Y302F showed impaired ability to recruit β-arrestin in response to Sar1-Ile4-Ile8 (SII) Ang II, a ligand solely activating the β-arrestin pathway.

CONCLUSIONS/SIGNIFICANCE

Our analysis reveals that the underlying conformations induced by these AT1R mutants most likely represent principally different mechanisms of uncoupling the G protein, which for some mutants may be due to their increased ability to recruit β-arrestin2. Hereby, these findings have important implications for drug discovery and 7TMR biology and illustrate the necessity of uncovering the exact molecular determinants for G protein-coupling and β-arrestin recruitment, respectively.

Theoretical studies on the interaction of partial agonists with the 5-HT2A receptor

A series of 51 5-HT(2A) partial agonistic arylethylamines (primary or benzylamines) from different structural classes (indoles, methoxybenzenes, quinazolinediones) was investigated by fragment regression analysis (FRA), docking and 3D-QSAR approaches. The data, pEC(50) values and intrinsic activities (E(max)) on rat arteries, show high variability of pEC(50) from 4 to 10 and of E(max) from 15 to 70%. FRA indicates which substructures affect potency or intrinsic activity. The high contribution of halogens in para position of phenethylamines to pEC(50) points to a specific hydrophobic pocket. Other results suggest the significance of hydrogen bonds of the aryl moiety for activation and the contrary effect of benzyl groups on affinity (increasing) and intrinsic activity (decreasing). Results from fragment regression and data on all available mutants were considered to derive a common binding site at the rat 5-HT(2A) receptor. After generation and MD simulations of a receptor model based on the β(2)-adrenoceptor structure, typical derivatives were docked, leading to the suggestion of common interactions, e.g., with serines in TM3 and TM5 and with a cluster of aromatic amino acids in TM5 and TM6. The whole series was aligned by docking and minimization of the complexes. The pEC(50) values correlate well with Sybyl docking energies and hydrophobicity of the aryl moieties. With this alignment, CoMFA and CoMSIA approaches based on a training set of 36 and a test set of 15 compounds were performed. The correlation of pEC(50) with steric, electrostatic, hydrophobic and H-bond acceptor fields resulted in sufficient fit (q (2): 0.75-0.8, r (2): 0.92-0.95) and predictive power (r (pred) (2) : 0.85-0.88). The important interaction regions largely reflect the patterns provided by the putative binding site. In particular, the fit of the aryl moieties and benzyl substituents to two hydrophobic pockets is evident.

Conserved water-mediated hydrogen bond network between TM-I, -II, -VI, and -VII in 7TM receptor activation

Five highly conserved polar residues connected by a number of structural water molecules together with two rotamer micro-switches, TrpVI:13 and TyrVII:20, constitute an extended hydrogen bond network between the intracellular segments of TM-I, -II, -VI, and -VII of 7TM receptors. Molecular dynamics simulations showed that, although the fewer water molecules in rhodopsin were relatively movable, the hydrogen bond network of the beta2-adrenergic receptor was fully loaded with water molecules that were surprisingly immobilized between the two rotamer switches, both apparently being in their closed conformation. Manipulations of the rotamer state of TyrVII:20 and TrpVI:13 demonstrated that these residues served as gates for the water molecules at the intracellular and extracellular ends of the hydrogen bond network, respectively. TrpVI:13 at the bottom of the main ligand-binding pocket was shown to apparently function as a catching trap for water molecules. Mutational analysis of the beta2-adrenergic receptor demonstrated that the highly conserved polar residues of the hydrogen bond network were all important for receptor signaling but served different functions, some dampening constitutive activity (AsnI:18, AspII:10, and AsnVII:13), whereas others (AsnVII:12 and AsnVII:16) located one helical turn apart and sharing a water molecule were shown to be essential for agonist-induced signaling. It is concluded that the conserved water hydrogen bond network of 7TM receptors constitutes an extended allosteric interface between the transmembrane segments being of crucial importance for receptor signaling and that part of the function of the rotamer micro-switches, TyrVII:20 and TrpVI:13, is to gate or trap the water molecules.

Molecular modeling and docking studies of human 5-hydroxytryptamine 2A (5-HT2A) receptor for the identification of hotspots for ligand binding

The serotonergic system has been implicated in emotional and cognitive function. In particular, 5-HT(2A) (5-hydroxytrytamine receptor 2A) is attributed to a number of disorders like schizophrenia, depression, eating disorders and anxiety. 5-HT(2A), being a GPCR (G-protein coupled receptor), is important in the pharmaceutical industry as a proven target for these disorders. Despite their extensive clinical importance, the structural studies of this protein is lacking due to difficulties in determining its crystal structure. We have performed sequence analysis and molecular modeling of 5-HT(2A) that has revealed a set of conserved residues and motifs considered to play an important role in maintaining structural integrity and function of the receptor. The analysis also revealed a set of residues specific to the receptor which distinguishes them from other members of the subclass and their orthologs. Further, starting from the model structure of human 5-HT(2A) receptor, docking studies were attempted to envisage how it might interact with eight of its ligands (such as serotonin, dopamine, DOI, LSD, haloperidol, ketanserin, risperidone and clozapine). The binding studies of dopamine to 5-HT(2A) receptor can bring up better understanding in the etiology of a number of neurological disorders involving both these two receptors. Our sequence analysis and study of interactions of this receptor with other ligands reveal additional residue hotspots such as Asn 363 and Tyr 370. The function of these residues can be further analyzed by rational design of site-directed mutagenesis. Two distinct binding sites are identified which could play important roles in ligand binding and signaling.

An indel in transmembrane helix 2 helps to trace the molecular evolution of class A G-protein-coupled receptors

Class A G-protein-coupled receptors (GPCRs) constitute a large family of transmembrane receptors. Helical distortions play a major role in the overall fold of these receptors. Most are related to conserved proline residues. However, in transmembrane helix 2, the proline pattern is not conserved, and when present, proline may be located at position 2.58, 2.59, or 2.60. Sequence analysis, three-dimensional data mining, and molecular modeling were undertaken to investigate the origin of this unusual pattern. Taken together, the data strongly support the assumption that an indel led to two structural motifs for helix 2: a bulged structure in P2.59 and P2.60 receptors and a "typical" proline kink in P2.58 receptors. The proline pattern of helix 2 can be used as an evolutionary marker and helps to trace the molecular evolution of class A GPCRs. Two indel events yielding functional receptors occurred independently. One indel arose very early in GPCR evolution, in a bilaterian ancestor, before the protostome-deuterostome divergence. This indel led to the split between the P2.58 somatostatin/opioid receptors and other peptide receptors with the P2.59 pattern. A second indel also occurred in insect opsins and corresponds to a deletion. Subfamilies with proline at position 2.59 or no proline expanded earlier, whereas P2.60 receptors remained marginal throughout evolution. P2.58 receptors underwent rapid expansion in vertebrates with the development of the chemokine and purinergic receptor subfamilies from somatostatin/opioid-related ancestors.

Fatty acid binding receptors and their physiological role in type 2 diabetes

G-protein-coupled receptors (GPCRs) respond to various physiological ligands such as photons, ions, and small molecules that include amines, fatty acids, and amino acids to peptides, proteins and steroids. Therefore, this family of proteins represents an attractive target for biopharmaceutical research. The physiological role of fatty acids and other lipid molecules as important signal mediators is well studied in various metabolic pathways. Acute administration of free fatty acids (FFAs) stimulates insulin release. Conversely, chronic exposure to high levels of free fatty acids leads to impairment of beta cell function and lipotoxicity. However, the receptors through which these fatty acids and lipids act were unknown, until the identification of fatty acid binding receptors: GPR40, GPR41, GPR43, and GPR119. Based on their tissue-expression profile, and pharmacologic analysis, the fatty acid binding receptors along with lipid binding receptor GPR119 are linked to diabetes and obesity. They play a critical role in the metabolic regulation of insulin release and glucose homeostasis. In this review, the mechanism of receptor activation, pharmacology, and the physiological functions of the fatty acid binding receptors will be discussed.

Mutational analysis of the conserved Asp2.50 and ERY motif reveals signaling bias of the urotensin II receptor

Class A (rhodopsin-like) G protein-coupled receptors possess conserved residues and motifs that are important for their specific activity. In the present study, we examined the role of residue Asp97(2.50) as well as residues Glu147(3.49), Arg148(3.50), and Tyr149(3.51) of the ERY motif on the functionality of the urotensin II receptor (UT). Mutations D97(2.50)A, R148(3.50)A, and R148(3.50)H abolished the ability of UT to activate phospholipase C, whereas mutations E147(3.49)A and Y149(3.51)A reduced the ability to activate PLC by 50%. None of the mutants exhibited constitutive activity. However, R148(3.50)A and R148(3.50)H promoted ERK1/2 activation, which was abolished by 4-(3-chloroanilino)-6,7-dimethoxyquinazoline (AG1478), an inhibitor of epidermal growth factor receptor (EGFR) tyrosine kinase activity. Both these mutants were capable of directly activating EGFR, which confirmed that they activated the mitogen-activated protein kinase (MAPK) pathway by a Galpha(q/11)-independent transactivation of EGFR. The D97(2.50)A, R148(3.50)A, and R148(3.50)H mutants did not readily internalize and did not promote translocation or colocalize with beta-arrestin2-GFP. Finally, the agonist-induced internalization of the E147(3.49)A mutant receptor was significantly increased compared with wild-type receptor. This study highlights the major contribution of the conserved Asp(2.50) residue to the functionality of the UT receptor. The Arg residue in the ERY motif of UT is an important structural element in signaling crossroads that determine whether Galpha(q/11)-dependent and -independent events can occur.

An intracellular interaction network regulates conformational transitions in the dopamine transporter

Neurotransmitter:sodium symporters (NSS)(1) mediate sodium-dependent reuptake of neurotransmitters from the synaptic cleft and are targets for many psychoactive drugs. The crystal structure of the prokaryotic NSS protein, LeuT, was recently solved at high resolution; however, the mechanistic details of regulation of the permeation pathway in this class of proteins remain unknown. Here we combine computational modeling and experimental probing in the dopamine transporter (DAT) to demonstrate the functional importance of a conserved intracellular interaction network. Our data suggest that a salt bridge between Arg-60 in the N terminus close to the cytoplasmic end of transmembrane segment (TM) 1 and Asp-436 at the cytoplasmic end of TM8 is stabilized by a cation-pi interaction between Arg-60 and Tyr-335 at the cytoplasmic end of TM6. Computational probing illustrates how the interactions may determine the flexibility of the permeation pathway, and mutagenesis within the network and results from assays of transport, as well as the state-dependent accessibility of a substituted cysteine in TM3, support the role of this network in regulating access between the substrate binding site and the intracellular milieu. The mechanism that emerges from these findings may be unique to the NSS family, where the local disruption of ionic interactions modulates the transition of the transporter between the outward- and inward-facing conformations.

Computational modeling study of functional microdomains in cannabinoid receptor type 1

The seven transmembrane helices (TMH) G-protein-coupled receptors (GPCRs) constitute one of the largest superfamily of signaling proteins found in mammals. Some of its members, in which the cannabinoid (CB) receptors are included, stand out because their functional states can be modulated by a broad spectrum of effector molecules. The relative ligand promiscuity exhibited by these receptors could be related with particular attributes conferred by their molecular architecture and represents a motivating issue to be explored. In this regard, this study represents an effort to investigate the cannabinoid receptor type 1 (CB1) ligand recognition plasticity, using comparative modeling, molecular dynamics (MD) simulations and docking. Our results suggest that a cooperative set of subtle structural rearrangements within the TMHs provide to the CB1 protein the plasticity to reach alternate configurations. These changes include the relaxation of intramolecular constraints, the rotations, translations and kinks of the majority of TMHs and the reorganization of the ligand binding cavities.

Mapping the structural requirements in the CB1 cannabinoid receptor transmembrane helix II for signal transduction

Amino acid residues in the transmembrane domains of the CB(1) receptor are important for ligand recognition and signal transduction. We used site-directed mutagenesis to identify the role of two novel and adjacent residues in the transmembrane helix II domain, Ile2.62 and Asp2.63. We investigated the role of the conserved, negatively charged aspartate at position 2.63 in cannabinoid receptor (CB(1)) function by substituting it with asparagine (D2.63N) and glutamate (D2.63E). In addition, the effect of the mutant I2.62T alone and in combination with D2.63N (double mutant) on the affinity and potency of structurally diverse ligands was investigated. Recombinant human CB(1) receptors, stably expressed in human embryonic kidney 293 cells, were assayed for ligand affinity and agonist-stimulated guanosine 5'-3-O-(thio)triphosphate (GTPgammaS) binding. The charge-conserved mutant D2.63E behaved similar to wild type. The charge-neutralization mutation D2.63N attenuated the potency of (-)-3-[2-hydroxyl-4-(1,1-dimethylheptyl)phenyl]-4-[3-hydroxylpropyl] cyclohexan-1-ol (CP,55940), (R)-(-)-[2,3-dihydro-5-methyl-3-[(4-morpholinyl)methyl]pyrrolo[1,2,3-de]-1,4-benzoxazin-6-yl](1-naphthalenyl)methanone (WIN55212-2), (-)-11beta-hydroxy-3-(1',1'-dimethylheptyl) hexahydrocannabinol (AM4056), and (-)-11-hydroxyldimethylheptyl-Delta(8)-tetrahydrocannabinol (HU210) for the stimulation of GTPgammaS binding, without affecting their binding affinities. Likewise, the I2.62T mutant selectively altered agonist potency without altering agonist affinity. It was surprising to note that the double mutant (I2.62T-D2.63N) displayed a drastic and synergistic increase (by approximately 50-fold) in the EC(50) for agonist-mediated activation. The profound loss of function in the I2.62T-D2.63N double mutant suggests that, although these residues are not obligatory for agonist recognition, they play a synergistic and crucial role in modulating signal transduction.

Genes encoding putative biogenic amine receptors in the parasitic nematode Brugia malayi

Filarial nematodes, such as Brugia malayi, cause major health problems worldwide. The lack of a vaccine against B. malayi, combined with ineffective chemotherapy against the adult has prompted the examination of biogenic amine receptors (BARs) as possible targets for drug discovery. We employed bioinformatics to identify genes encoding putative B. malayi BARs. Surprisingly, the B. malayi genome contains half of the genes predicted to encode BARs in the genomes of free-living nematodes such as Caenorhabditis elegans or C. briggsae; however, all of the predicted B. malayi receptors have clear orthologues in C. elegans. The B. malayi genes encode each of the major BAR subclasses, including three serotonin, two dopamine and two tyramine/octopamine receptors and the structure of orthologous BAR genes is conserved. We find that potential G-protein coupling and ligand-specificity of individual BARs may be predicted by phylogenetic comparisons. Our results provide a framework for how G-protein coupled receptors may be targeted for drug development in medically important parasitic nematodes.

Functional analysis of transmembrane domain 2 of the M1 muscarinic acetylcholine receptor

Ala substitution scanning mutagenesis has been used to probe the functional role of amino acids in transmembrane (TM) domain 2 of the M1 muscarinic acetylcholine receptor, and of the highly conserved Asn43 in TM1. The mutation of Asn43, Asn61, and Leu64 caused an enhanced ACh affinity phenotype. Interpreted using a rhodopsin-based homology model, these results suggest the presence of a network of specific contacts between this group of residues and Pro415 and Tyr418 in the highly conserved NPXXY motif in TM7 that exhibit a similar mutagenic phenotype. These contacts may be rearranged or broken when ACh binds. D71A, like N414A, was devoid of signaling activity. We suggest that formation of a direct hydrogen bond between the highly conserved side chains of Asp71 and Asn414 may be a critical feature stabilizing the activated state of the M1 receptor. Mutation of Leu67, Ala70, and Ile74 also reduced the signaling efficacy of the ACh-receptor complex. The side chains of these residues are modeled as an extended surface that may help to orient and insulate the proposed hydrogen bond between Asp71 and Asn414. Mutation of Leu72, Gly75, and Met79 in the outer half of TM2 primarily reduced the expression of functional receptor binding sites. These residues may mediate contacts with TM1 and TM7 that are preserved throughout the receptor activation cycle. Thermal inactivation measurements confirmed that a reduction in structural stability followed the mutation of Met79 as well as Asp71.

The Angiotensin II AT1 Receptor Structure-Activity Correlations in the Light of Rhodopsin Structure

The most prevalent physiological effects of ANG II, the main product of the renin-angiotensin system, are mediated by the AT1 receptor, a rhodopsin-like AGPCR. Numerous studies of the cardiovascular effects of synthetic peptide analogs allowed a detailed mapping of ANG II's structural requirements for receptor binding and activation, which were complemented by site-directed mutagenesis studies on the AT1 receptor to investigate the role of its structure in ligand binding, signal transduction, phosphorylation, binding to arrestins, internalization, desensitization, tachyphylaxis, and other properties. The knowledge of the high-resolution structure of rhodopsin allowed homology modeling of the AT1 receptor. The models thus built and mutagenesis data indicate that physiological (agonist binding) or constitutive (mutated receptor) activation may involve different degrees of expansion of the receptor's central cavity. Residues in ANG II structure seem to control these conformational changes and to dictate the type of cytosolic event elicited during the activation. 1) Agonist aromatic residues (Phe8 and Tyr4) favor the coupling to G protein, and 2) absence of these residues can favor a mechanism leading directly to receptor internalization via phosphorylation by specific kinases of the receptor's COOH-terminal Ser and Thr residues, arrestin binding, and clathrin-dependent coated-pit vesicles. On the other hand, the NH2-terminal residues of the agonists ANG II and [Sar1]-ANG II were found to bind by two distinct modes to the AT1 receptor extracellular site flanked by the COOH-terminal segments of the EC-3 loop and the NH2-terminal domain. Since the [Sar1]-ligand is the most potent molecule to trigger tachyphylaxis in AT1 receptors, it was suggested that its corresponding binding mode might be associated with this special condition of receptors.

Critical role for polar residues in coupling leukotriene B4 binding to signal transduction in BLT1

Leukotriene B(4) (LTB(4)) mediates a variety of inflammatory diseases such as asthma, arthritis, atherosclerosis, and cancer through activation of the G-protein-coupled receptor, BLT1. Using in silico molecular dynamics simulations combined with site-directed mutagenesis we characterized the ligand binding site and activation mechanism for BLT1. Mutation of residues predicted as potential ligand contact points in transmembrane domains (TMs) III (H94A and Y102A), V (E185A), and VI (N241A) resulted in reduced binding affinity. Analysis of arginines in extracellular loop 2 revealed that mutating arginine 156 but not arginine 171 or 178 to alanine resulted in complete loss of LTB(4) binding to BLT1. Structural models for the ligand-free and ligand-bound states of BLT1 revealed an activation core formed around Asp-64, displaying multiple dynamic interactions with Asn-36, Ser-100, and Asn-281 and a triad of serines, Ser-276, Ser-277, and Ser-278. Mutagenesis of many of these residues in BLT1 resulted in loss of signaling capacity while retaining normal LTB(4) binding function. Thus, polar residues within TMs III, V, and VI and extracellular loop 2 are critical for ligand binding, whereas polar residues in TMs II, III, and VII play a central role in transducing the ligand-induced conformational change to activation. The delineation of a validated binding site and activation mechanism should facilitate structure-based design of inhibitors targeting BLT1.

Receptor-based pharmacophores for serotonin 5-HT7R antagonists-implications to selectivity

A set of 31 diversified 5-HT7 receptor antagonists was automatically docked to a conformational ensemble of rhodopsin-based 5-HT7R models (flexible docking). It was found that ergolines, aporphines, and tricyclic psychotropic agents were always docked in a pocket formed by TMHs 4-6, and besides the main ionic bond with Asp3.32, they had specific interactions with Phe6.52, Phe6.51, Trp6.48, and Ser5.42. The arylpiperidine, arylpiperazine, or beta-carboline fragment of the complex ligands occupied the same pocket, whereas the terminal amide/imide part of those compounds reached the second cavity formed by TMHs 7-3 and interacted with Phe3.28, Arg7.36, and Tyr7.43. A similar orientation of selective antagonists of the group of arylsulfonamidoalkylamines was observed, that is, the sulfonamide part was located in the second pocket. Coherent docking results allowed the generation of two receptor-based pharmacophores: first containing features necessary for high 5-HT7R affinity and the other defining selectivity for this receptor subtype. The latter model indicated the importance of specific interactions with the residues of the TMHs 7-3 pocket (especially nonconserved Arg7.36) for selectivity over other monoamine GPCRs.