Crystal structure of squid rhodopsin with intracellularly extended cytoplasmic region

G Protein-Coupled Receptor-Ligand Pose and Functional Class Prediction

G protein-coupled receptor (GPCR) transmembrane protein family members play essential roles in physiology. Numerous pharmaceuticals target GPCRs, and many drug discovery programs utilize virtual screening (VS) against GPCR targets. Improvements in the accuracy of predicting new molecules that bind to and either activate or inhibit GPCR function would accelerate such drug discovery programs. This work addresses two significant research questions. First, do ligand interaction fingerprints provide a substantial advantage over automated methods of binding site selection for classical docking? Second, can the functional status of prospective screening candidates be predicted from ligand interaction fingerprints using a random forest classifier? Ligand interaction fingerprints were found to offer modest advantages in sampling accurate poses, but no substantial advantage in the final set of top-ranked poses after scoring, and, thus, were not used in the generation of the ligand-receptor complexes used to train and test the random forest classifier. A binary classifier which treated agonists, antagonists, and inverse agonists as active and all other ligands as inactive proved highly effective in ligand function prediction in an external test set of GPR31 and TAAR2 candidate ligands with a hit rate of 82.6% actual actives within the set of predicted actives.

Transfer of knowledge from model organisms to evolutionarily distant non-model organisms: The coral Pocillopora damicornis membrane signaling receptome

With the ease of gene sequencing and the technology available to study and manipulate non-model organisms, the extension of the methodological toolbox required to translate our understanding of model organisms to non-model organisms has become an urgent problem. For example, mining of large coral and their symbiont sequence data is a challenge, but also provides an opportunity for understanding functionality and evolution of these and other non-model organisms. Much more information than for any other eukaryotic species is available for humans, especially related to signal transduction and diseases. However, the coral cnidarian host and human have diverged over 700 million years ago and homologies between proteins in the two species are therefore often in the gray zone, or at least often undetectable with traditional BLAST searches. We introduce a two-stage approach to identifying putative coral homologues of human proteins. First, through remote homology detection using Hidden Markov Models, we identify candidate human homologues in the cnidarian genome. However, for many proteins, the human genome alone contains multiple family members with similar or even more divergence in sequence. In the second stage, therefore, we filter the remote homology results based on the functional and structural plausibility of each coral candidate, shortlisting the coral proteins likely to have conserved some of the functions of the human proteins. We demonstrate our approach with a pipeline for mapping membrane receptors in humans to membrane receptors in corals, with specific focus on the stony coral, P. damicornis. More than 1000 human membrane receptors mapped to 335 coral receptors, including 151 G protein coupled receptors (GPCRs). To validate specific sub-families, we chose opsin proteins, representative GPCRs that confer light sensitivity, and Toll-like receptors, representative non-GPCRs, which function in the immune response, and their ability to communicate with microorganisms. Through detailed structure-function analysis of their ligand-binding pockets and downstream signaling cascades, we selected those candidate remote homologues likely to carry out related functions in the corals. This pipeline may prove generally useful for other non-model organisms, such as to support the growing field of synthetic biology.

Rhodopsins: An Excitingly Versatile Protein Species for Research, Development and Creative Engineering

The first member and eponym of the rhodopsin family was identified in the 1930s as the visual pigment of the rod photoreceptor cell in the animal retina. It was found to be a membrane protein, owing its photosensitivity to the presence of a covalently bound chromophoric group. This group, derived from vitamin A, was appropriately dubbed retinal. In the 1970s a microbial counterpart of this species was discovered in an archaeon, being a membrane protein also harbouring retinal as a chromophore, and named bacteriorhodopsin. Since their discovery a photogenic panorama unfolded, where up to date new members and subspecies with a variety of light-driven functionality have been added to this family. The animal branch, meanwhile categorized as type-2 rhodopsins, turned out to form a large subclass in the superfamily of G protein-coupled receptors and are essential to multiple elements of light-dependent animal sensory physiology. The microbial branch, the type-1 rhodopsins, largely function as light-driven ion pumps or channels, but also contain sensory-active and enzyme-sustaining subspecies. In this review we will follow the development of this exciting membrane protein panorama in a representative number of highlights and will present a prospect of their extraordinary future potential.

Additive and epistatic effects influence spectral tuning in molluscan retinochrome opsin

The relationship between genotype and phenotype is nontrivial due to often complex molecular pathways that make it difficult to unambiguously relate phenotypes to specific genotypes. Photopigments, an opsin apoprotein bound to a light-absorbing chromophore, present an opportunity to directly relate the amino acid sequence to an absorbance peak phenotype (λmax). We examined this relationship by conducting a series of site-directed mutagenesis experiments of retinochrome, a non-visual opsin, from two closely related species: the common bay scallop, Argopecten irradians, and the king scallop, Pecten maximus. Using protein folding models, we identified three amino acid sites of likely functional importance and expressed mutated retinochrome proteins in vitro. Our results show that the mutation of amino acids lining the opsin binding pocket are responsible for fine spectral tuning, or small changes in the λmax of these light sensitive proteins Mutations resulted in a blue or red shift as predicted, but with dissimilar magnitudes. Shifts ranged from a 16 nm blue shift to a 12 nm red shift from the wild-type λmax. These mutations do not show an additive effect, but rather suggests the presence of epistatic interactions. This work highlights the importance of binding pocket shape in the evolution of spectral tuning and builds on our ability to relate genotypic changes to phenotypes in an emerging model for opsin functional analysis.

The C-Terminus and Third Cytoplasmic Loop Cooperatively Activate Mouse Melanopsin Phototransduction

Melanopsin, an atypical vertebrate visual pigment, mediates non-image-forming light responses including circadian photoentrainment and pupillary light reflexes and contrast detection for image formation. Melanopsin-expressing intrinsically photosensitive retinal ganglion cells are characterized by sluggish activation and deactivation of their light responses. The molecular determinants of mouse melanopsin's deactivation have been characterized (i.e., C-terminal phosphorylation and β-arrestin binding), but a detailed analysis of melanopsin's activation is lacking. We propose that an extended third cytoplasmic loop is adjacent to the proximal C-terminal region of mouse melanopsin in the inactive conformation, which is stabilized by the ionic interaction of these two regions. This model is supported by site-directed spin labeling and electron paramagnetic resonance spectroscopy of melanopsin, the results of which suggests a high degree of steric freedom at the third cytoplasmic loop, which is increased upon C-terminus truncation, supporting the idea that these two regions are close in three-dimensional space in wild-type melanopsin. To test for a functionally critical C-terminal conformation, calcium imaging of melanopsin mutants including a proximal C-terminus truncation (at residue 365) and proline mutation of this proximal region (H377P, L380P, Y382P) delayed melanopsin's activation rate. Mutation of all potential phosphorylation sites, including a highly conserved tyrosine residue (Y382), into alanines also delayed the activation rate. A comparison of mouse melanopsin with armadillo melanopsin-which has substitutions of various potential phosphorylation sites and a substitution of the conserved tyrosine-indicates that substitution of these potential phosphorylation sites and the tyrosine residue result in dramatically slower activation kinetics, a finding that also supports the role of phosphorylation in signaling activation. We therefore propose that melanopsin's C-terminus is proximal to intracellular loop 3, and C-terminal phosphorylation permits the ionic interaction between these two regions, thus forming a stable structural conformation that is critical for initiating G-protein signaling.

Protein Interactors and Trafficking Pathways That Regulate the Cannabinoid Type 1 Receptor (CB1R)

The endocannabinoid system (ECS) acts as a negative feedback mechanism to suppress synaptic transmission and plays a major role in a diverse range of brain functions including, for example, the regulation of mood, energy balance, and learning and memory. The function and dysfunction of the ECS are strongly implicated in multiple psychiatric, neurological, and neurodegenerative diseases. Cannabinoid type 1 receptor (CB1R) is the most abundant G protein-coupled receptor (GPCR) expressed in the brain and, as for any synaptic receptor, CB1R needs to be in the right place at the right time to respond appropriately to changing synaptic circumstances. While CB1R is found intracellularly throughout neurons, its surface expression is highly polarized to the axonal membrane, consistent with its functional expression at presynaptic sites. Surprisingly, despite the importance of CB1R, the interacting proteins and molecular mechanisms that regulate the highly polarized distribution and function of CB1R remain relatively poorly understood. Here we set out what is currently known about the trafficking pathways and protein interactions that underpin the surface expression and axonal polarity of CB1R, and highlight key questions that still need to be addressed.

Melanopsin Carboxy-terminus phosphorylation plasticity and bulk negative charge, not strict site specificity, achieves phototransduction deactivation

Melanopsin is a visual pigment expressed in a small subset of ganglion cells in the mammalian retina known as intrinsically photosensitive retinal ganglion cells (ipRGCs) and is implicated in regulating non-image forming functions such as circadian photoentrainment and pupil constriction and contrast sensitivity in image formation. Mouse melanopsin's Carboxy-terminus (C-terminus) possesses 38 serine and threonine residues, which can potentially serve as phosphorylation sites for a G-protein Receptor Kinase (GRK) and be involved in the deactivation of signal transduction. Previous studies suggest that S388, T389, S391, S392, S394, S395 on the proximal region of the C-terminus of mouse melanopsin are necessary for melanopsin deactivation. We expressed a series of mouse melanopsin C-terminal mutants in HEK293 cells and using calcium imaging, and we found that the necessary cluster of six serine and threonine residues, while being critical, are insufficient for proper melanopsin deactivation. Interestingly, the additional six serine and threonine residues adjacent to the required six sites, in either proximal or distal direction, are capable of restoring wild-type deactivation of melanopsin. These findings suggest an element of plasticity in the molecular basis of melanopsin phosphorylation and deactivation. In addition, C-terminal chimeric mutants and molecular modeling studies support the idea that the initial steps of deactivation and β-arrestin binding are centered around these critical phosphorylation sites (S388-S395). The degree of functional versatility described in this study, along with ipRGC biophysical heterogeneity and the possible use of multiple signal transduction cascades, might contribute to the diverse ipRGC light responses for use in non-image and image forming behaviors, even though all six sub types of ipRGCs express the same melanopsin gene OPN4.

Update of the GRIP web service

G protein-coupled receptors (GPCRs) can form homodimers, heterodimers, or higher-order molecular complexes (oligomers). The reports on the change of functions through the oligomerization have been accumulated. Inhibition of GPCR oligomerization without affecting the protomer's overall structure would clarify the oligomer-specific functions although inhibition experiments are costly and require accurate information about the interface location. Unfortunately, the number of experimentally determined interfaces is limited. The precise prediction of the oligomerization interfaces is, therefore, useful for inhibition experiments to examine the oligomer-specific functions, which would accelerate investigations of the GPCR signaling. However, interface prediction for GPCR oligomerization is difficult because different GPCR subtypes belonging to the same subfamily often use different structural regions as their interfaces. We previously developed a high-performance method to predict the interfaces for GPCR oligomerization, by identifying the conserved surfaces with the sequence and structure information. Then, the structural characteristic of a GPCR structure is regarded to be a thick-tube like conformation that is approximately perpendicular to the membrane plane. Our method had successfully predicted all of the interfaces available on that day. We had launched a web server for our interface prediction of GPCRs (GRIP). We have improved the previous version of GRIP server and enhanced its usability. First, we discarded the approximation of the GPCR structure as the thick-tube-like conformation. This improvement increased the number of structures for the prediction. Second, the FUGUE-based template recommendation service was introduced to facilitate the choice of an appropriate structure for the prediction. The new prediction server is available at http://grip.b.dendai.ac.jp/∼grip/.

Improved homology modeling of the human & rat EP4 prostanoid receptors

Background

The EP₄ prostanoid receptor is one of four GPCRs that mediate the diverse actions of prostaglandin E₂ (PGE₂). Novel selective EP₄ receptor agonists would assist to further elucidate receptor sub-type function and promote development of therapeutics for bone healing, heart failure, and other receptor associated conditions. The rat EP₄ (rEP₄) receptor has been used as a surrogate for the human EP₄ (hEP₄) receptor in multiple SAR studies. To better understand the validity of this traditional approach, homology models were generated by threading for both receptors using the RaptorX server. These models were fit to an implicit membrane using the PPM server and OPM database with refinement of intra and extracellular loops by Prime (Schrödinger). To understand the interaction between the receptors and known agonists, induced-fit docking experiments were performed using Glide and Prime (Schrödinger), with both endogenous agonists and receptor sub-type selective, small-molecule agonists. The docking scores and observed interactions were compared with radioligand displacement experiments and receptor (rat & human) activation assays monitoring cAMP.

Results

Rank-ordering of in silico compound docking scores aligned well with in vitro activity assay EC₅₀ and radioligand binding K_i. We observed variations between rat and human EP₄ binding pockets that have implications in future small-molecule receptor-modulator design and SAR, specifically a S103G mutation within the rEP4 receptor. Additionally, these models helped identify key interactions between the EP₄ receptor and ligands including PGE₂ and several known sub-type selective agonists while serving as a marked improvement over the previously reported models.

Conclusions

This work has generated a set of novel homology models of the rEP₄ and hEP₄ receptors. The homology models provide an improvement upon the previously reported model, largely due to improved solvation. The hEP₄ docking scores correlates best with the cAMP activation data, where both data sets rank order Rivenprost>CAY10684 > PGE₁ ≈ PGE₂ > 11-deoxy-PGE₁ ≈ 11-dexoy-PGE₂ > 8-aza-11-deoxy-PGE₁. This rank-ordering matches closely with the rEP₄ receptor as well. Species-specific differences were noted for the weak agonists Sulprostone and Misoprostol, which appear to dock more readily within human receptor versus rat receptor.

Electronic supplementary material

The online version of this article (10.1186/s12860-019-0212-5) contains supplementary material, which is available to authorized users.

Crystal structure of jumping spider rhodopsin-1 as a light sensitive GPCR

Light-sensitive G protein-coupled receptors (GPCRs)-rhodopsins-absorb photons to isomerize their covalently bound retinal, triggering conformational changes that result in downstream signaling cascades. Monostable rhodopsins release retinal upon isomerization as opposed to the retinal in bistable rhodopsins that "reisomerize" upon absorption of a second photon. Understanding the mechanistic differences between these light-sensitive GPCRs has been hindered by the scarcity of recombinant models of the latter. Here, we reveal the high-resolution crystal structure of a recombinant bistable rhodopsin, jumping spider rhodopsin-1, bound to the inverse agonist 9-cis retinal. We observe a water-mediated network around the ligand hinting toward the basis of their bistable nature. In contrast to bovine rhodopsin (monostable), the transmembrane bundle of jumping spider rhodopsin-1 as well that of the bistable squid rhodopsin adopts a more "activation-ready" conformation often observed in other nonphotosensitive class A GPCRs. These similarities suggest the role of jumping spider rhodopsin-1 as a potential model system in the study of the structure-function relationship of both photosensitive and nonphotosensitive class A GPCRs.

The counterion-retinylidene Schiff base interaction of an invertebrate rhodopsin rearranges upon light activation

Animals sense light using photosensitive proteins-rhodopsins-containing a chromophore-retinal-that intrinsically absorbs in the ultraviolet. Visible light-sensitivity depends primarily on protonation of the retinylidene Schiff base (SB), which requires a negatively-charged amino acid residue-counterion-for stabilization. Little is known about how the most common counterion among varied rhodopsins, Glu181, functions. Here, we demonstrate that in a spider visual rhodopsin, orthologue of mammal melanopsins relevant to circadian rhythms, the Glu181 counterion functions likely by forming a hydrogen-bonding network, where Ser186 is a key mediator of the Glu181-SB interaction. We also suggest that upon light activation, the Glu181-SB interaction rearranges while Ser186 changes its contribution. This is in contrast to how the counterion of vertebrate visual rhodopsins, Glu113, functions, which forms a salt bridge with the SB. Our results shed light on the molecular mechanisms of visible light-sensitivity relevant to invertebrate vision and vertebrate non-visual photoreception.

Sequence, Structure, and Expression of Opsins in the Monochromatic Stomatopod Squilla empusa

Most stomatopod crustaceans have complex retinas in their compound eyes, with up to 16 spectral types of photoreceptors, but members of the superfamily Squilloidea have much simpler retinas, thought to contain a single photoreceptor spectral class. In the Atlantic stomatopod Squilla empusa, microspectrophotometry shows that all photoreceptors absorb light maximally at 517 nm, indicating that a single visual pigment is present in all photoreceptors in the retina. However, six distinct, but partial, long wavelength sensitive (LWS) opsin transcripts, which encode the protein component of the visual pigment, have been previously isolated through RT-PCR. In order to investigate the spectral and functional differences among S. empusa's opsins, we used RT-PCR to complete the 3' end of sequences for five of the six expressed opsins. The extended sequences spanned from the first transmembrane (TM1) helix to the 3' end of the coding region. Using homology-based modeling, we predicted the three-dimensional structure of the amino acid translation of the S. empusa opsins. Based on these analyses, S. empusa LWS opsins share a high sequence identity in TM regions and in amino acids within 15 Å of the chromophore-binding lysine on TM helix 7 (TM7), suggesting that these opsins produce spectrally similar visual pigments in agreement with previous results. However, we propose that these spectrally similar opsins differ functionally, as there are non-conservative amino acid substitutions found in intracellular loop 2 (ICL2) and TM5/ICL3, which are critical regions for G-protein binding, and substitutions in extracellular regions suggest different chromophore attachment affinities. In situ hybridization of two of the opsins (Se5 and Se6) revealed strong co-expression in all photoreceptors in both midband and peripheral regions of the retina as well as in selected ocular and cerebral ganglion neuropils. These data suggest the expression of multiple opsins-likely spectrally identical, but functionally different-in multiple types of neuronal cells in S. empusa. This suggests that the multiple opsins characteristic of other stomatopod species may have similar functional specialization.

The Two-Photon Reversible Reaction of the Bistable Jumping Spider Rhodopsin-1

Bistable opsins are photopigments expressed in both invertebrates and vertebrates. These light-sensitive G-protein-coupled receptors undergo a reversible reaction upon illumination. A first photon initiates the cis to trans isomerization of the retinal chromophore—attached to the protein through a protonated Schiff base—and a series of transition states that eventually results in the formation of the thermally stable and active Meta state. Excitation by a second photon reverts this process to recover the original ground state. On the other hand, monostable opsins (e.g., bovine rhodopsin) lose their chromophore during the decay of the Meta II state (i.e., they bleach). Spectroscopic studies on the molecular details of the two-photon cycle in bistable opsins are limited. Here, we describe the successful expression and purification of recombinant rhodopsin-1 from the jumping spider Hasarius adansoni (JSR1). In its natural configuration, spectroscopic characterization of JSR1 is hampered by the similar absorption spectra in the visible spectrum of the inactive and active states. We solved this issue by separating their absorption spectra by replacing the endogenous 11-cis retinal chromophore with the blue-shifted 9-cis JSiR1. With this system, we used time-resolved ultraviolet-visible spectroscopy after pulsed laser excitation to obtain kinetic details of the rise and decay of the photocycle intermediates. We also used resonance Raman spectroscopy to elucidate structural changes of the retinal chromophore upon illumination. Our data clearly indicate that the protonated Schiff base is stable throughout the entire photoreaction. We additionally show that the accompanying conformational changes in the protein are different from those of monostable rhodopsin, as recorded by light-induced FTIR difference spectroscopy. Thus, we envisage JSR1 as becoming a model system for future studies on the reaction mechanisms of bistable opsins, e.g., by time-resolved x-ray crystallography.

A Novel Polar Core and Weakly Fixed C-Tail in Squid Arrestin Provide New Insight into Interaction with Rhodopsin

Photoreceptors of the squid Loligo pealei contain a G-protein-coupled receptor (GPCR) signaling system that activates phospholipase C in response to light. Analogous to the mammalian visual system, signaling of the photoactivated GPCR rhodopsin is terminated by binding of squid arrestin (sArr). sArr forms a light-dependent, high-affinity complex with squid rhodopsin, which does not require prior receptor phosphorylation for interaction. This is at odds with classical mammalian GPCR desensitization where an agonist-bound phosphorylated receptor is needed to break stabilizing constraints within arrestins, the so-called "three-element interaction" and "polar core" network, before a stable receptor-arrestin complex can be established. Biophysical and mass spectrometric analysis of the squid rhodopsin-arrestin complex indicates that in contrast to mammalian arrestins, the sArr C-tail is not involved in a stable three-element interaction. We determined the crystal structure of C-terminally truncated sArr that adopts a basal conformation common to arrestins and is stabilized by a series of weak but novel polar core interactions. Unlike mammalian arrestin-1, deletion of the sArr C-tail does not influence kinetic properties of complex formation of sArr with the receptor. Hydrogen-deuterium exchange studies revealed the footprint of the light-activated rhodopsin on sArr. Furthermore, double electron-electron resonance spectroscopy experiments provide evidence that receptor-bound sArr adopts a conformation different from the one known for arrestin-1 and molecular dynamics simulations reveal the residues that account for the weak three-element interaction. Insights gleaned from studying this system add to our general understanding of GPCR-arrestin interaction.

Identifying G protein-coupled receptor dimers from crystal packings

Dimers of G protein-coupled receptors (GPCRs) are believed to be important for signaling with their associated G proteins. Low-resolution electron microscopy has shown rhodopsin dimers in native retinal membranes, and CXCR4 dimers have been found in several different crystal structures. Evidence for dimers of other GPCRs is more indirect. An alternative to computational modeling studies is to search for parallel dimers in the packing environments of the reported crystal structures of GPCRs. Two major structural types of GPCR dimers exist (as predicted by others), but there is considerable structural variation within each cluster. The different structural variants described here might reflect different functional properties and should provide a range of model structures for computational and experimental examination.

In Silico Studies Targeting G-protein Coupled Receptors for Drug Research Against Parkinson's Disease

Parkinson's Disease (PD) is a long-term neurodegenerative brain disorder that mainly affects the motor system. The causes are still unknown, and even though currently there is no cure, several therapeutic options are available to manage its symptoms. The development of novel antiparkinsonian agents and an understanding of their proper and optimal use are, indeed, highly demanding. For the last decades, L-3,4-DihydrOxyPhenylAlanine or levodopa (L-DOPA) has been the gold-standard therapy for the symptomatic treatment of motor dysfunctions associated to PD. However, the development of dyskinesias and motor fluctuations (wearing-off and on-off phenomena) associated with long-term L-DOPA replacement therapy have limited its antiparkinsonian efficacy. The investigation for non-dopaminergic therapies has been largely explored as an attempt to counteract the motor side effects associated with dopamine replacement therapy. Being one of the largest cell membrane protein families, G-Protein-Coupled Receptors (GPCRs) have become a relevant target for drug discovery focused on a wide range of therapeutic areas, including Central Nervous System (CNS) diseases. The modulation of specific GPCRs potentially implicated in PD, excluding dopamine receptors, may provide promising non-dopaminergic therapeutic alternatives for symptomatic treatment of PD. In this review, we focused on the impact of specific GPCR subclasses, including dopamine receptors, adenosine receptors, muscarinic acetylcholine receptors, metabotropic glutamate receptors, and 5-hydroxytryptamine receptors, on the pathophysiology of PD and the importance of structure- and ligand-based in silico approaches for the development of small molecules to target these receptors.

Optogenetic interrogation reveals separable G-protein-dependent and -independent signalling linking G-protein-coupled receptors to the circadian oscillator

BACKGROUND

Endogenous circadian oscillators distributed across the mammalian body are synchronised among themselves and with external time via a variety of signalling molecules, some of which interact with G-protein-coupled receptors (GPCRs). GPCRs can regulate cell physiology via pathways originating with heterotrimeric G-proteins or β-arrestins. We applied an optogenetic approach to determine the contribution of these two signalling modes on circadian phase.

RESULTS

We employed a photopigment (JellyOp) that activates Gαs signalling with better selectivity and higher sensitivity than available alternatives, and a point mutant of this pigment (F112A) biased towards β-arrestin signalling. When expressed in fibroblasts, both native JellyOp and the F112A arrestin-biased mutant drove light-dependent phase resetting in the circadian clock. Shifts induced by the two opsins differed in their circadian phase dependence and the degree to which they were associated with clock gene induction.

CONCLUSIONS

Our data imply separable G-protein and arrestin inputs to the mammalian circadian clock and establish a pair of optogenetic tools suitable for manipulating Gαs- and β-arrestin-biased signalling in live cells.

Structural differences and differential expression among rhabdomeric opsins reveal functional change after gene duplication in the bay scallop, Argopecten irradians (Pectinidae)

Background

Opsins are the only class of proteins used for light perception in image-forming eyes. Gene duplication and subsequent functional divergence of opsins have played an important role in expanding photoreceptive capabilities of organisms by altering what wavelengths of light are absorbed by photoreceptors (spectral tuning). However, new opsin copies may also acquire novel function or subdivide ancestral functions through changes to temporal, spatial or the level of gene expression. Here, we test how opsin gene copies diversify in function and evolutionary fate by characterizing four rhabdomeric (G_q-protein coupled) opsins in the scallop, Argopecten irradians, identified from tissue-specific transcriptomes.

Results

Under a phylogenetic analysis, we recovered a pattern consistent with two rounds of duplication that generated the genetic diversity of scallop G_q-opsins. We found strong support for differential expression of paralogous G_q-opsins across ocular and extra-ocular photosensitive tissues, suggesting that scallop G_q-opsins are used in different biological contexts due to molecular alternations outside and within the protein-coding regions. Finally, we used available protein models to predict which amino acid residues interact with the light-absorbing chromophore. Variation in these residues suggests that the four G_q-opsin paralogs absorb different wavelengths of light.

Conclusions

Our results uncover novel genetic and functional diversity in the light-sensing structures of the scallop, demonstrating the complicated nature of G_q-opsin diversification after gene duplication. Our results highlight a change in the nearly ubiquitous shadow response in molluscs to a narrowed functional specificity for visual processes in the eyed scallop. Our findings provide a starting point to study how gene duplication may coincide with eye evolution, and more specifically, different ways neofunctionalization of G_q-opsins may occur.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0823-9) contains supplementary material, which is available to authorized users.

Targeting Liver Fibrosis with a Cell-penetrating Protease-activated Receptor-2 (PAR2) Pepducin

Chronic liver inflammation and fibrosis in nonalcoholic steatohepatitis can lead to cirrhosis and liver failure for which there are currently no approved treatments. Protease-activated receptor-2 (PAR2) is an emerging new target expressed on liver stellate cells and hepatocytes that regulates the response to liver injury and inflammation. Here, we identified a pepducin to block the deleterious actions of PAR2 in promoting liver fibrosis. Non-alcoholic fatty liver disease and early fibrosis were induced by the methionine-choline-deficient diet in mice. Fibrotic liver disease was induced by administering carbon tetrachloride for 8 weeks. Mice were treated with the pepducin PZ-235 either from onset of the experiment or after fibrosis was established. Hepatic fibrosis, collagen content, inflammatory cytokines, steatosis, triglycerides, and NAFLD activity score were assessed as primary outcome parameters depending on the model. The activity of the PAR2 pepducin on cultured stellate cell activation and hepatocyte reactive oxygen species production was evaluated. PZ-235 significantly suppressed liver fibrosis, collagen deposition, inflammatory cytokines, NAFLD activity score, steatosis, triglycerides, aspartate transaminase, alanine transaminase, and stellate cell proliferation by up to 50-100%. The PAR2 inhibitor afforded significant protective effects against hepatocellular necrosis and attenuated PAR2-mediated reactive oxygen species production in hepatocytes. PZ-235 was distributed to liver and other mouse tissues and was found to form a well structured α-helix that closely resembles the juxtamembrane helical region of the analogous TM6 and third intracellular region of the intact receptor that is critical for coupling to internal G proteins. The ability of PZ-235 to effectively suppress fibrosis, hepatocellular necrosis, reactive oxygen species production, steatosis, and inflammation indicates the potential for PAR2 pepducin inhibitors to be broadly efficacious in the treatment of liver fibrosis.

The leukotriene B4 receptor BLT1 is stabilized by transmembrane helix capping mutations

In this study, we introduced structure-based rational mutations in the guinea pig leukotriene B₄ receptor (gpBLT1) in order to enhance the stabilization of the protein. Elements thought to be unfavorable for the stability of gpBLT1 were extracted based on the stabilization elements established in soluble proteins, determined crystal structures of G-protein-coupled receptors (GPCRs), and multiple sequence alignment. The two unfavorable residues His83^2.67 and Lys88^3.21, located at helix capping sites, were replaced with Gly (His83Gly^2.67 and Lys88Gly^3.21). The modified protein containing His83Gly^2.67/Lys88Gly^3.21 was highly expressed, solubilized, and purified and exhibited improved thermal stability by 4 °C in comparison with that of the original gpBLT1 construct. Owing to the double mutation, the expression level increased by 6-fold (B_max=311 pmol/mg) in the membrane fraction of Pichia pastoris. The ligand binding affinity was similar to that of the original gpBLT1 without the mutations. Similar unfavorable residues have been observed at helix capping sites in many other GPCRs; therefore, the replacement of such residues with more favorable residues will improve stabilization of the GPCR structure for the crystallization.

•

Point mutations were rationally designed to stabilize LTB₄ receptor (BLT1).

•

The stability of mutant His83Gly^2.67/Lys88Gly^3.21 improved by 5 °C.

•

BLT1 expression by P. pastoris was increased 6-fold.

•

Mutations were designed to replace unfavorable residues at the helix capping site.

•

This method would be useful for the stabilization of the other membrane proteins.

Large-scale production and protein engineering of G protein-coupled receptors for structural studies

Structural studies of G protein-coupled receptors (GPCRs) gave insights into molecular mechanisms of their action and contributed significantly to molecular pharmacology. This is primarily due to technical advances in protein engineering, production and crystallization of these important receptor targets. On the other hand, NMR spectroscopy of GPCRs, which can provide information about their dynamics, still remains challenging due to difficulties in preparation of isotopically labeled receptors and their low long-term stabilities. In this review, we discuss methods used for expression and purification of GPCRs for crystallographic and NMR studies. We also summarize protein engineering methods that played a crucial role in obtaining GPCR crystal structures.

Predicted structures for kappa opioid G-protein coupled receptor bound to selective agonists

Human kappa opioid receptor (κ-OR), a G protein-coupled receptor (GPCR), has been identified as a drug target for treatment of such human disorders as pain perception, neuroendocrine physiology, affective behavior, and cognition. In order to find more selective and active agonists, one would like to do structure based drug design. Indeed, there is an X-ray structure for an antagonist bound to κ-OR, but structures for activated GPCRs are quite different from those for the inactive GPCRs. Here we predict the ensemble of 24 low-energy structures of human kappa opioid receptor (κ-OR), obtained by application of the GEnSeMBLE (GPCR Ensemble of Structures in Membrane Bilayer Environment) complete sampling method, which evaluates 13 trillion combinations of tilt and rotation angles for κ-OR to select the best 24. To validate these structures, we used the DarwinDock complete sampling method to predict the binding sites for five known agonists (ethylketocyclazocine, bremazocine, pentazocine, nalorphine, and morphine) bound to all 24 κ-OR conformations. We find that some agonists bind selectively to receptor conformations that lack the salt bridge between transmembrane domains 3 and 6 as expected for active conformations. These 3D structures for κ-OR provide a structural basis for understanding ligand binding and activation of κ-OR, which should be useful for guiding subtype specific drug design.

Molecular basis for dramatic changes in cannabinoid CB1 G protein-coupled receptor activation upon single and double point mutations

There is considerable interest in determining the activation mechanism of G protein-coupled receptors (GPCRs), one of the most important types of proteins for intercellular signaling. Recently, it was demonstrated for the cannabinoid CB1 GPCR, that a single mutation T210A could make CB1 completely inactive whereas T210I makes it essentially constitutively active. To obtain an understanding of this dramatic dependence of activity on mutation, we used first-principles-based methods to predict the ensemble of low-energy seven-helix conformations for the wild-type (WT) and mutants (T210A and T210I). We find that the transmembrane (TM) helix packings depend markedly on these mutations, leading for T210A to both TM3+TM6 and TM2+TM6 salt-bridge couplings in the cytoplasmic face that explains the inactivity of this mutant. In contrast T210I has no such couplings across the receptor explaining the ease in activating this mutant. WT has just the TM3+TM6 coupling, known to be broken upon GPCR activation. To test this hypothesis on activity, we predicted double mutants that would convert the inactive mutant to normal activity and then confirmed this experimentally. This CB1 activation mechanism, or one similar to it, is expected to play a role in other constitutively active GPCRs as well.

Direct imaging electron microscopy (EM) methods in modern structural biology: overview and comparison with X-ray crystallography and single-particle cryo-EM reconstruction in the studies of large macromolecules

Determining the structure of macromolecules is important for understanding their function. The fine structure of large macromolecules is currently studied primarily by X-ray crystallography and single-particle cryo-electron microscopy (EM) reconstruction. Before the development of these techniques, macromolecular structure was often examined by negative-staining, rotary-shadowing and freeze-etching EM, which are categorised here as 'direct imaging EM methods'. In this review, the results are summarised by each of the above techniques and compared with respect to four macromolecules: the ryanodine receptor, cadherin, rhodopsin and the ribosome-translocon complex (RTC). The results of structural analysis of the ryanodine receptor and cadherin are consistent between each technique. The results obtained for rhodopsin vary to some extent within each technique and between the different techniques. Finally, the results for RTC are inconsistent between direct imaging EM and other analytical techniques, especially with respect to the space within RTC, the reasons for which are discussed. Then, the role of direct imaging EM methods in modern structural biology is discussed. Direct imaging methods should support and verify the results obtained by other analytical methods capable of solving three-dimensional molecular architecture, and they should still be used as a primary tool for studying macromolecule structure in vivo.

Computational analysis of the 3-D structure of human GPR87 protein: implications for structure-based drug design

The G-protein coupled receptor 87 (GPR87) is a recently discovered orphan GPCR which means that the search of their endogenous ligands has been a novel challenge. GPR87 has been shown to be overexpressed in squamous cell carcinomas (SCCs) or adenocarcinomas in lungs and bladder. The 3D structure of GPR87 was here modeled using two templates (2VT4 and 2ZIY) by a threading method. Functional assignment of GPR87 by SVM revealed that along with transporter activity, various novel functions were predicted. The 3D structure was further validated by comparison with structural features of the templates through Verify-3D, ProSA and ERRAT for determining correct stereochemical parameters. The resulting model was evaluated by Ramachandran plot and good 3D structure compatibility was evidenced by DOPE score. Molecular dynamics simulation and solvation of protein were studied through explicit spherical boundaries with a harmonic restraint membrane water system. A DRY-motif (Asp-Arg-Tyr sequence) was found at the end of transmembrane helix3, where GPCR binds and thus activation of signals is transduced. In a search for better inhibitors of GPR87, in silico modification of some substrate ligands was carried out to form polar interactions with Arg115 and Lys296. Thus, this study provides early insights into the structure of a major drug target for SCCs.

The heterodimeric glycoprotein hormone, GPA2/GPB5, regulates ion transport across the hindgut of the adult mosquito, Aedes aegypti

A family of evolutionarily old hormones is the glycoprotein cysteine knot-forming heterodimers consisting of alpha- (GPA) and beta-subunits (GPB), which assemble by noncovalent bonds. In mammals, a common glycoprotein hormone alpha-subunit (GPA1) pairs with unique beta-subunits that establish receptor specificity, forming thyroid stimulating hormone (GPA1/TSHβ) and the gonadotropins luteinizing hormone (GPA1/LHβ), follicle stimulating hormone (GPA1/FSHβ), choriogonadotropin (GPA1/CGβ). A novel glycoprotein heterodimer was identified in vertebrates by genome analysis, called thyrostimulin, composed of two novel subunits, GPA2 and GPB5, and homologs occur in arthropods, nematodes and cnidarians, implying that this neurohormone system existed prior to the emergence of bilateral metazoans. In order to discern possible physiological roles of this hormonal signaling system in mosquitoes, we have isolated the glycoprotein hormone genes producing the alpha- and beta-subunits (AedaeGPA2 and AedaeGPB5) and assessed their temporal expression profiles in the yellow and dengue-fever vector, Aedes aegypti. We have also isolated a putative receptor for this novel mosquito hormone, AedaeLGR1, which contains features conserved with other glycoprotein leucine-rich repeating containing G protein-coupled receptors. AedaeLGR1 is expressed in tissues of the alimentary canal such as the midgut, Malpighian tubules and hindgut, suggesting that this novel mosquito glycoprotein hormone may regulate ionic and osmotic balance. Focusing on the hindgut in adult stage A. aegypti, where AedaeLGR1 was highly enriched, we utilized the Scanning Ion-selective Electrode Technique (SIET) to determine if AedaeGPA2/GPB5 modulated cation transport across this epithelial tissue. Our results suggest that AedaeGPA2/GPB5 does indeed participate in ionic and osmotic balance, since it appears to inhibit natriuresis and promote kaliuresis. Taken together, our findings imply this hormone may play an important role in ionic balance when levels of Na⁺ are limited and levels of K⁺ are in excess – such as during the digestion and assimilation of erythrocytes following vertebrate blood-feeding by females.

Structural approaches to understanding retinal proteins needed for vision

The past decade has witnessed an impressive expansion of our knowledge of retinal photoreceptor signal transduction and the regulation of the visual cycle required for normal eyesight. Progress in human genetics and next generation sequencing technologies have revealed the complexity behind many inherited retinal diseases. Structural studies have markedly increased our understanding of the visual process. Moreover, technical innovations and improved methodologies in proteomics, macromolecular crystallization and high resolution imaging at different levels set the scene for even greater advances. Pharmacology combined with structural biology of membrane proteins holds great promise for developing innovative accessible therapies for millions robbed of their sight or progressing toward blindness.

GPCR activation: protonation and membrane potential

GPCR proteins represent the largest family of signaling membrane proteins in eukaryotic cells. Their importance to basic cell biology, human diseases, and pharmaceutical interventions is well established. Many crystal structures of GPCR proteins have been reported in both active and inactive conformations. These data indicate that agonist binding alone is not sufficient to trigger the conformational change of GPCRs necessary for binding of downstream G-proteins, yet other essential factors remain elusive. Based on analysis of available GPCR crystal structures, we identified a potential conformational switch around the conserved Asp2.50, which consistently shows distinct conformations between inactive and active states. Combining the structural information with the current literature, we propose an energy-coupling mechanism, in which the interaction between a charge change of the GPCR protein and the membrane potential of the living cell plays a key role for GPCR activation.

Structural features of the G-protein/GPCR interactions

BACKGROUND

The details of the functional interaction between G proteins and the G protein coupled receptors (GPCRs) have long been subjected to extensive investigations with structural and functional assays and a large number of computational studies.

SCOPE OF REVIEW

The nature and sites of interaction in the G-protein/GPCR complexes, and the specificities of these interactions selecting coupling partners among the large number of families of GPCRs and G protein forms, are still poorly defined.

MAJOR CONCLUSIONS

Many of the contact sites between the two proteins in specific complexes have been identified, but the three dimensional molecular architecture of a receptor-Gα interface is only known for one pair. Consequently, many fundamental questions regarding this macromolecular assembly and its mechanism remain unanswered.

GENERAL SIGNIFICANCE

In the context of current structural data we review the structural details of the interfaces and recognition sites in complexes of sub-family A GPCRs with cognate G-proteins, with special emphasis on the consequences of activation on GPCR structure, the prevalence of preassembled GPCR/G-protein complexes, the key structural determinants for selective coupling and the possible involvement of GPCR oligomerization in this process.

From atomic structures to neuronal functions of g protein-coupled receptors

G protein-coupled receptors (GPCRs) are essential mediators of signal transduction, neurotransmission, ion channel regulation, and other cellular events. GPCRs are activated by diverse stimuli, including light, enzymatic processing of their N-termini, and binding of proteins, peptides, or small molecules such as neurotransmitters. GPCR dysfunction caused by receptor mutations and environmental challenges contributes to many neurological diseases. Moreover, modern genetic technology has helped identify a rich array of mono- and multigenic defects in humans and animal models that connect such receptor dysfunction with disease affecting neuronal function. The visual system is especially suited to investigate GPCR structure and function because advanced imaging techniques permit structural studies of photoreceptor neurons at both macro and molecular levels that, together with biochemical and physiological assessment in animal models, provide a more complete understanding of GPCR signaling.

X-ray structural information of GPCRs in drug design: what are the limitations and where do we go?

INTRODUCTION

In 2007, the X-ray structural determination of non-rhodopsin G-Protein coupled receptors (GPCRs), considered the most extensively targeted protein class for marketed drugs, commenced. With the relatively rapid availability of additional structures, an assessment of the progression made is needed in addition to the assessment of the understandings gleaned, deployment successes and forthcoming prospects.

AREAS COVERED

The author reviews the approaches and tools that have made it possible to determine the three dimensional structures of GPCRs using X-ray crystallography. Furthermore, the author describes the methods suited for crystallization of membrane bound GPCR proteins including the lipidic cubic phase and various protein modification approaches. The author also provides highlights, from the literature, of the structures determined to date including targets solved, the nature of the content provided (such as selectivity, activating vs. inactivating determinants) and how these structural features relate to drug design strategies.

EXPERT OPINION

The GPCR X-ray structures that have been so far determined have yielded significant information. This has presented dramatic evidence concerning their ability to impact the discovery of compounds through their action as traditional, orthosteric modulators. It is, however, noted that more challenging design strategies, such as identifying biased agonists and the use of sites remote from the orthosteric site for allosteric modulation, are still in their infancy.

A Strategy Combining Differential Low-Throughput Screening and Virtual Screening (DLS-VS) Accelerating the Discovery of new Modulators for the Orphan GPR34 Receptor

The DLS-VS strategy was developed as an integrated method for identifying chemical modulators for orphan GPCRs. It combines differential low-throughput screening (DLS) and virtual screening (VS). The two cascaded techniques offer complementary advantages and allow the experimental testing of a minimal number of compounds. First, DLS identifies modulators specific for the considered receptor among a set of receptors, through the screening of a small library with diverse chemical compounds. Then, an active molecular model of the receptor is built by homology to a validated template, and it is progressively refined by rotamers modification for key side-chains, by VS of the already screened library, and by iterative selection of the model generating the best enrichment. The refined active model is finally used for the VS of a large chemical library and the selection of a small set of compounds for experimental testing. Applied to the orphan receptor GPR34, the DLS-VS strategy combined the experimental screening of 20 000 compounds and the virtual screening of 1 250 000 compounds. It identified one agonist and eight inverse agonists, showing a high chemical diversity. We describe the method. The strategy can be applied to other GPCRs.

The role of gene duplication and unconstrained selective pressures in the melanopsin gene family evolution and vertebrate circadian rhythm regulation

Melanopsin is a photosensitive cell protein involved in regulating circadian rhythms and other non-visual responses to light. The melanopsin gene family is represented by two paralogs, OPN4x and OPN4m, which originated through gene duplication early in the emergence of vertebrates. Here we studied the melanopsin gene family using an integrated gene/protein evolutionary approach, which revealed that the rhabdomeric urbilaterian ancestor had the same amino acid patterns (DRY motif and the Y and E conterions) as extant vertebrate species, suggesting that the mechanism for light detection and regulation is similar to rhabdomeric rhodopsins. Both OPN4m and OPN4x paralogs are found in vertebrate genomic paralogons, suggesting that they diverged following this duplication event about 600 million years ago, when the complex eye emerged in the vertebrate ancestor. Melanopsins generally evolved under negative selection (ω = 0.171) with some minor episodes of positive selection (proportion of sites = 25%) and functional divergence (θ(I) = 0.349 and θ(II) = 0.126). The OPN4m and OPN4x melanopsin paralogs show evidence of spectral divergence at sites likely involved in melanopsin light absorbance (200F, 273S and 276A). Also, following the teleost lineage-specific whole genome duplication (3R) that prompted the teleost fish radiation, type I divergence (θ(I) = 0.181) and positive selection (affecting 11% of sites) contributed to amino acid variability that we related with the photo-activation stability of melanopsin. The melanopsin intracellular regions had unexpectedly high variability in their coupling specificity of G-proteins and we propose that Gq/11 and Gi/o are the two G-proteins most-likely to mediate the melanopsin phototransduction pathway. The selection signatures were mainly observed on retinal-related sites and the third and second intracellular loops, demonstrating the physiological plasticity of the melanopsin protein group. Our results provide new insights on the phototransduction process and additional tools for disentangling and understanding the links between melanopsin gene evolution and the specializations observed in vertebrates, especially in teleost fish.

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.

New insights for drug design from the X-ray crystallographic structures of G-protein-coupled receptors

Methodological advances in X-ray crystallography have made possible the recent solution of X-ray structures of pharmaceutically important G protein-coupled receptors (GPCRs), including receptors for biogenic amines, peptides, a nucleoside, and a sphingolipid. These high-resolution structures have greatly increased our understanding of ligand recognition and receptor activation. Conformational changes associated with activation common to several receptors entail outward movements of the intracellular side of transmembrane helix 6 (TM6) and movements of TM5 toward TM6. Movements associated with specific agonists or receptors have also been described [e.g., extracellular loop (EL) 3 in the A(2A) adenosine receptor]. The binding sites of different receptors partly overlap but differ significantly in ligand orientation, depth, and breadth of contact areas in TM regions and the involvement of the ELs. A current challenge is how to use this structural information for the rational design of novel potent and selective ligands. For example, new chemotypes were discovered as antagonists of various GPCRs by subjecting chemical libraries to in silico docking in the X-ray structures. The vast majority of GPCR structures and their ligand complexes are still unsolved, and no structures are known outside of family A GPCRs. Molecular modeling, informed by supporting information from site-directed mutagenesis and structure-activity relationships, has been validated as a useful tool to extend structural insights to related GPCRs and to analyze docking of other ligands in already crystallized GPCRs.

Platform for the rapid construction and evaluation of GPCRs for crystallography in Saccharomyces cerevisiae

BACKGROUND

Recent successes in the determination of G-protein coupled receptor (GPCR) structures have relied on the ability of receptor variants to overcome difficulties in expression and purification. Therefore, the quick screening of functionally expressed stable receptor variants is vital.

RESULTS

We developed a platform using Saccharomyces cerevisiae for the rapid construction and evaluation of functional GPCR variants for structural studies. This platform enables us to perform a screening cycle from construction to evaluation of variants within 6-7 days. We firstly confirmed the functional expression of 25 full-length class A GPCRs in this platform. Then, in order to improve the expression level and stability, we generated and evaluated the variants of the four GPCRs (hADRB2, hCHRM2, hHRH1 and hNTSR1). These stabilized receptor variants improved both functional activity and monodispersity. Finally, the expression level of the stabilized hHRH1 in Pichia pastoris was improved up to 65 pmol/mg from negligible expression of the functional full-length receptor in S. cerevisiae at first screening. The stabilized hHRH1 was able to be purified for use in crystallization trials.

CONCLUSIONS

We demonstrated that the S. cerevisiae system should serve as an easy-to-handle and rapid platform for the construction and evaluation of GPCR variants. This platform can be a powerful prescreening method to identify a suitable GPCR variant for crystallography.

Transmembrane helix: simple or complex

Transmembrane helical segments (TMs) can be classified into two groups of so-called ‘simple’ and ‘complex’ TMs. Whereas the first group represents mere hydrophobic anchors with an overrepresentation of aliphatic hydrophobic residues that are likely attributed to convergent evolution in many cases, the complex ones embody ancestral information and tend to have structural and functional roles beyond just membrane immersion. Hence, the sequence homology concept is not applicable on simple TMs. In practice, these simple TMs can attract statistically significant but evolutionarily unrelated hits during similarity searches (whether through BLAST- or HMM-based approaches). This is especially problematic for membrane proteins that contain both globular segments and TMs. As such, we have developed the transmembrane helix: simple or complex (TMSOC) webserver for the identification of simple and complex TMs. By masking simple TM segments in seed sequences prior to sequence similarity searches, the false-discovery rate decreases without sacrificing sensitivity. Therefore, TMSOC is a novel and necessary sequence analytic tool for both the experimentalists and the computational biology community working on membrane proteins. It is freely accessible at http://tmsoc.bii.a-star.edu.sg or available for download.

Molecular dynamics simulations reveal specific interactions of post-translational palmitoyl modifications with rhodopsin in membranes

We present a detailed analysis of the behavior of the highly flexible post-translational lipid modifications of rhodopsin from multiple-microsecond all-atom molecular dynamics simulations. Rhodopsin was studied in a realistic membrane environment that includes cholesterol, as well as saturated and polyunsaturated lipids with phosphocholine and phosphoethanolamine headgroups. The simulation reveals striking differences between the palmitoylations at Cys322 and Cys323 as well as between the palmitoyl chains and the neighboring lipids. Notably the palmitoyl group at Cys322 shows considerably greater contact with helix H1 of rhodopsin, yielding frequent chain upturns with longer reorientational correlation times, and relatively low order parameters. While the palmitoylation at Cys323 makes fewer protein contacts and has increased order compared to Cys322, it nevertheless exhibits greater flexibility with smaller order parameters than the stearoyl chains of the surrounding lipids. The dynamical structure of the palmitoylations-as well as their extensive fluctuations-suggests a complex function for the post-translational modifications in rhodopsin and potentially other G protein-coupled receptors, going beyond their role as membrane anchoring elements. Rather, we propose that the palmitoylation at Cys323 has a potential role as a lipid anchor, whereas the palmitoyl-protein interaction observed for Cys322 suggests a more specific interaction that affects the stability of the dark state of rhodopsin.

G-protein-coupled receptor inactivation by an allosteric inverse-agonist antibody

G protein-coupled receptors (GPCRs) are the largest class of cell-surface receptors, and these membrane proteins exist in equilibrium between inactive and active states.^1-13 Conformational changes induced by extracellular ligands binding to GPCRs result in a cellular response through the activation of G-proteins. The A_2A adenosine receptor (A_2AAR) is responsible for regulating blood flow to the cardiac muscle and is important in the regulation of glutamate and dopamine release in the brain.¹⁴ In this study, we have successfully raised a mouse monoclonal antibody against human A_2AAR that prevents agonist but not antagonist binding to the extracellular ligand-binding pocket. The structure of the A_2AAR-antibody Fab fragment (Fab2838) complex reveals that the fragment, unexpectedly, recognises the intracellular surface of A_2AAR and that its complementarity determining region, CDR-H3, penetrates into the receptor. CDR-H3 is located in a similar position to the G-protein C-terminal fragment in the active opsin structure¹ and to the CDR-3 of the nanobody in the active β₂ adrenergic receptor structure² but locks the A_2AAR in an inactive conformation. These results shed light on a novel strategy to modulate GPCR activity.

Comparative modeling of lipid receptors

Comparative modeling is a powerful technique to generate models of proteins from families already represented by members with experimentally characterized three-dimensional structures. The method is particularly important for modeling membrane-bound receptors in the G Protein-Coupled Receptor (GPCR) family, such as many of the lipid receptors (such as the cannabinoid, prostanoid, lysophosphatidic acid, sphingosine 1-phosphate, and eicosanoid receptor family members), as these represent particularly challenging targets for experimental structural characterization methods. Although challenging modeling targets, these receptors have been linked to therapeutic indications that vary from nociception to cancer, and thus are of interest as therapeutic targets. Accurate models of lipid receptors are therefore valuable tools in the drug discovery and optimization phases of therapeutic development. This chapter describes the construction and evaluation of comparative structural models of lipid receptors beginning with the selection of template structures.

Homology modeling of class a G protein-coupled receptors

G protein-coupled receptors (GPCRs) are a large superfamily of membrane bound signaling proteins that hold great pharmaceutical interest. Since experimentally elucidated structures are available only for a very limited number of receptors, homology modeling has become a widespread technique for the construction of GPCR models intended to study the structure-function relationships of the receptors and aid the discovery and development of ligands capable of modulating their activity. Through this chapter, various aspects involved in the constructions of homology models of the serpentine domain of the largest class of GPCRs, known as class A or rhodopsin family, are illustrated. In particular, the chapter provides suggestions, guidelines, and critical thoughts on some of the most crucial aspect of GPCR modeling, including: collection of candidate templates and a structure-based alignment of their sequences; identification and alignment of the transmembrane helices of the query receptor to the corresponding domains of the candidate templates; selection of one or more templates receptor; election of homology or de novo modeling for the construction of specific extracellular and intracellular domains; construction of the 3D models, with special consideration to extracellular regions, disulfide bridges, and interhelical cavity; validation of the models through controlled virtual screening experiments.

The significance of G protein-coupled receptor crystallography for drug discovery

Crucial as molecular sensors for many vital physiological processes, seven-transmembrane domain G protein-coupled receptors (GPCRs) comprise the largest family of proteins targeted by drug discovery. Together with structures of the prototypical GPCR rhodopsin, solved structures of other liganded GPCRs promise to provide insights into the structural basis of the superfamily's biochemical functions and assist in the development of new therapeutic modalities and drugs. One of the greatest technical and theoretical challenges to elucidating and exploiting structure-function relationships in these systems is the emerging concept of GPCR conformational flexibility and its cause-effect relationship for receptor-receptor and receptor-effector interactions. Such conformational changes can be subtle and triggered by relatively small binding energy effects, leading to full or partial efficacy in the activation or inactivation of the receptor system at large. Pharmacological dogma generally dictates that these changes manifest themselves through kinetic modulation of the receptor's G protein partners. Atomic resolution information derived from increasingly available receptor structures provides an entrée to the understanding of these events and practically applying it to drug design. Supported by structure-activity relationship information arising from empirical screening, a unified structural model of GPCR activation/inactivation promises to both accelerate drug discovery in this field and improve our fundamental understanding of structure-based drug design in general. This review discusses fundamental problems that persist in drug design and GPCR structural determination.

Characterizing and predicting the functional and conformational diversity of seven-transmembrane proteins

The activation of seven-transmembrane receptors (7TMRs) allows cells to sense their environment and convert extracellular signals (like hormone binding) into intracellular signals (through G protein-coupled and/or β arrestin-coupled pathways). A single 7TMR is capable of transducing a wide spectrum of physiological responses inside a cell by coupling to these pathways. This intracellular pleiotropic action is enabled by multiple conformations exhibited by these receptors. Developments in membrane protein structure determination technologies have led to a rapid increase in crystal structures for many 7TMRs. Majority of these receptors have been crystallized in their inactive conformation and, for some, one of the many active conformations has also been crystallized. Given the topological constraints of a lipid bilayer that results in a single fold of seven almost parallel TM helices connected by mostly unstructured loops, these structures exhibit a diversity of conformations not only across the receptors but also across the different functional forms for receptors with structures for one of the functionally active conformations. Here we present a method to characterize this conformational diversity in terms of transmembrane helix topology (TMHTOP) parameters and how to use these helix orientation parameters to predict functionally-distinct multiple conformations for these receptors. The TMHTOP parameters enable a quantification of the structural changes that underlie 7TMR activation and also sheds a unique mechanistic light on the pleiotropic nature of these receptors. It provides a common language to describe the 7TMR activation mechanisms as well as differences across many receptors in terms of visually intuitive structural parameters. Protein structure prediction methods can use these parameters to describe 7TMR conformational ensembles, which coupled to experimental data can be used to develop testable hypotheses for the structural basis of 7TMR functions.

Preparation of an activated rhodopsin/transducin complex using a constitutively active mutant of rhodopsin

The interaction of rhodopsin and transducin has been the focus of study for more than 30 years, but only recently have efforts to purify an activated complex in detergent solution materialized. These efforts have used native rhodopsin isolated from bovine retina and employed either sucrose density gradient centrifugation or size exclusion chromatography to purify the complex. While there is general agreement on most properties of the activated complex, subunit stoichiometry is not yet settled, with rhodopsin/transducin molar ratios of both 2/1 and 1/1 reported. In this report, we introduce methods for preparation of the complex that include use of recombinant rhodopsin, so as to take advantage of mutations that confer constitutive activity and enhanced thermal stability on the protein, and immunoaffinity chromatography for purification of the complex. We show that chromatography on ConA-Sepharose can substitute for the immunoaffinity column and that bicelles can be used instead of detergent solution. We demonstrate the following: that rhodopsin has a covalently bound all-trans-retinal chromophore and therefore corresponds to the active metarhodopin II state; that transducin has an empty nucleotide-binding pocket; that the isolated complex is active and dissociates upon addition of guanine nucleotide; and finally that the stoichiometry corresponds reproducibly to a 1/1 molar ratio of rhodopsin to transducin.

Diversity and modularity of G protein-coupled receptor structures

G protein-coupled receptors (GPCRs) comprise the most 'prolific' family of cell membrane proteins, which share a general mechanism of signal transduction, but greatly vary in ligand recognition and function. Crystal structures are now available for rhodopsin, adrenergic, and adenosine receptors in both inactive and activated forms, as well as for chemokine, dopamine, and histamine receptors in inactive conformations. Here we review common structural features, outline the scope of structural diversity of GPCRs at different levels of homology, and briefly discuss the impact of the structures on drug discovery. Given the current set of GPCR crystal structures, a distinct modularity is now being observed between the extracellular (ligand-binding) and intracellular (signaling) regions. The rapidly expanding repertoire of GPCR structures provides a solid framework for experimental and molecular modeling studies, and helps to chart a roadmap for comprehensive structural coverage of the whole superfamily and an understanding of GPCR biological and therapeutic mechanisms.

Collation and data-mining of literature bioactivity data for drug discovery

The challenge of translating the huge amount of genomic and biochemical data into new drugs is a costly and challenging task. Historically, there has been comparatively little focus on linking the biochemical and chemical worlds. To address this need, we have developed ChEMBL, an online resource of small-molecule SAR (structure–activity relationship) data, which can be used to support chemical biology, lead discovery and target selection in drug discovery. The database contains the abstracted structures, properties and biological activities for over 700000 distinct compounds and in excess of more than 3 million bioactivity records abstracted from over 40000 publications. Additional public domain resources can be readily integrated into the same data model (e.g. PubChem BioAssay data). The compounds in ChEMBL are largely extracted from the primary medicinal chemistry literature, and are therefore usually ‘drug-like’ or ‘lead-like’ small molecules with full experimental context. The data cover a significant fraction of the discovery of modern drugs, and are useful in a wide range of drug design and discovery tasks. In addition to the compound data, ChEMBL also contains information for over 8000 protein, cell line and whole-organism ‘targets’, with over 4000 of those being proteins linked to their underlying genes. The database is searchable both chemically, using an interactive compound sketch tool, protein sequences, family hierarchies, SMILES strings, compound research codes and key words, and biologically, using a variety of gene identifiers, protein sequence similarity and protein families. The information retrieved can then be readily filtered and downloaded into various formats. ChEMBL can be accessed online at https://www.ebi.ac.uk/chembldb.