Two hypotheses--one answer. Sequence comparison does not support an evolutionary link between halobacterial retinal proteins including bacteriorhodopsin and eukaryotic G-protein-coupled receptors

Automatic Rhodopsin Modeling with Multiple Protonation Microstates

Automatic Rhodopsin Modeling (ARM) is a simulation protocol providing QM/MM models of rhodopsins capable of reproducing experimental electronic absorption and emission trends. Currently, ARM is restricted to a single protonation microstate for each rhodopsin model. Herein, we incorporate an extension of the minimal electrostatic model (MEM) into the ARM protocol to account for all relevant protonation microstates at a given pH. The new ARM+MEM protocol determines the most important microstates contributing to the description of the absorption spectrum. As a test case, we have applied this methodology to simulate the pH-dependent absorption spectrum of a toy model, showing that the single-microstate picture breaks down at certain pH values. Subsequently, we applied ARM+MEM toAnabaenasensory rhodopsin, confirming an improved description of its absorption spectrum when the titration of several key residues is considered.

Convergent evolution of animal and microbial rhodopsins

Rhodopsins, a family of photoreceptive membrane proteins, contain retinal as a chromophore and were firstly identified as reddish pigments from frog retina in 1876. Since then, rhodopsin-like proteins have been identified mainly from animal eyes. In 1971, a rhodopsin-like pigment was discovered from the archaeon Halobacterium salinarum and named bacteriorhodopsin. While it was believed that rhodopsin- and bacteriorhodopsin-like proteins were expressed only in animal eyes and archaea, respectively, before the 1990s, a variety of rhodopsin-like proteins (called animal rhodopsins or opsins) and bacteriorhodopsin-like proteins (called microbial rhodopsins) have been progressively identified from various tissues of animals and microorganisms, respectively. Here, we comprehensively introduce the research conducted on animal and microbial rhodopsins. Recent analysis has revealed that the two rhodopsin families have common molecular properties, such as the protein structure (i.e., 7-transmembrane structure), retinal structure (i.e., binding ability to cis- and trans-retinal), color sensitivity (i.e., UV- and visible-light sensitivities), and photoreaction (i.e., triggering structural changes by light and heat), more than what was expected at the early stages of rhodopsin research. Contrastingly, their molecular functions are distinctively different (e.g., G protein-coupled receptors and photoisomerases for animal rhodopsins and ion transporters and phototaxis sensors for microbial rhodopsins). Therefore, based on their similarities and dissimilarities, we propose that animal and microbial rhodopsins have convergently evolved from their distinctive origins as multi-colored retinal-binding membrane proteins whose activities are regulated by light and heat but independently evolved for different molecular and physiological functions in the cognate organism.

Superconserved receptors expressed in the brain: Expression, function, motifs and evolution of an orphan receptor family

GPR27, GPR85 and GPR173 constitute a small family of G protein-coupled receptors (GPCR) that share the distinctive characteristics of being highly conserved throughout vertebrate evolution and predominantly expressed in the brain. Accordingly, they have been coined as "Superconserved Receptors Expressed in the Brain" (SREB), although their expression profile is more complex than what was originally thought. SREBs have no known validated endogenous ligands and are thus labeled as "orphan" receptors. The investigation of this particular category of uncharacterized receptors holds great promise both in terms of physiology and drug development. In the largest GPCR family, the Rhodopsin-like or Class A, around 100 receptors are considered orphans. Because GPCRs are the most successful source of drug targets, the discovery of a novel function or ligand most likely will lead to significant breakthroughs for the discovery of innovative therapies. The high level of conservation is one of the characteristic features of the SREBs. We propose herein a detailed analysis of the putative evolutionary origin of this family. We highlight the properties that distinguish SREBs from other rhodopsin-like GPCRs. We present the current evidence for these receptors downstream signaling pathways and functions. We discuss the pharmacological challenge for the identification of natural or synthetic ligands of orphan receptors like SREBs. The different SREB-related scientific questions are presented with a highlight on what should be addressed in the near future, including the confirmation of published evidence and their validation as drug targets. In particular, we discuss in which pathological conditions these receptors may be of great relevance to solve unmet medical needs.

Expansion of the "Sodium World" through Evolutionary Time and Taxonomic Space

In 1986, Vladimir Skulachev and his colleagues coined the term "Sodium World" for the group of diverse organisms with sodium (Na)-based bioenergetics. Albeit only few such organisms had been discovered by that time, the authors insightfully noted that "the great taxonomic variety of organisms employing the Na-cycle points to the ubiquitous distribution of this novel type of membrane-linked energy transductions". Here we used tools of bioinformatics to follow expansion of the Sodium World through the evolutionary time and taxonomic space. We searched for those membrane protein families in prokaryotic genomes that correlate with the use of the Na-potential for ATP synthesis by different organisms. In addition to the known Na-translocators, we found a plethora of uncharacterized protein families; most of them show no homology with studied proteins. In addition, we traced the presence of Na-based energetics in many novel archaeal and bacterial clades, which were recently identified by metagenomic techniques. The data obtained support the view that the Na-based energetics preceded the proton-dependent energetics in evolution and prevailed during the first two billion years of the Earth history before the oxygenation of atmosphere. Hence, the full capacity of Na-based energetics in prokaryotes remains largely unexplored. The Sodium World expanded owing to the acquisition of new functions by Na-translocating systems. Specifically, most classes of G-protein-coupled receptors (GPCRs), which are targeted by almost half of the known drugs, appear to evolve from the Na-translocating microbial rhodopsins. Thereby the GPCRs of class A, with 700 representatives in human genome, retained the Na-binding site in the center of the transmembrane heptahelical bundle together with the capacity of Na-translocation. Mathematical modeling showed that the class A GPCRs could use the energy of transmembrane Na-potential for increasing both their sensitivity and selectivity. Thus, GPCRs, the largest protein family coded by human genome, stem from the Sodium World, which encourages exploration of other Na-dependent enzymes of eukaryotes.

Recent Advances and Applications of Molecular Docking to G Protein-Coupled Receptors

The growing number of studies on G protein-coupled receptors (GPCRs) family are a source of noticeable improvement in our understanding of the functioning of these proteins. GPCRs are responsible for a vast part of signaling in vertebrates and, as such, invariably remain in the spotlight of medicinal chemistry. A deeper insight into the underlying mechanisms of interesting phenomena observed in GPCRs, such as biased signaling or allosteric modulation, can be gained with experimental and computational studies. The latter play an important role in this process, since they allow for observations on scales inaccessible for most other methods. One of the key steps in such studies is proper computational reconstruction of actual ligand-receptor or protein-protein interactions, a process called molecular docking. A number of improvements and innovative applications of this method were documented recently. In this review, we focus particularly on innovations in docking to GPCRs.

Eukaryotic G protein-coupled receptors as descendants of prokaryotic sodium-translocating rhodopsins

Abstract

Microbial rhodopsins and G-protein coupled receptors (GPCRs, which include animal rhodopsins) are two distinct (super) families of heptahelical (7TM) membrane proteins that share obvious structural similarities but no significant sequence similarity. Comparison of the recently solved high-resolution structures of the sodium-translocating bacterial rhodopsin and various Na⁺-binding GPCRs revealed striking similarity of their sodium-binding sites. This similarity allowed us to construct a structure-guided sequence alignment for the two (super)families, which highlighted their evolutionary relatedness. Our analysis supports a common underlying molecular mechanism for both families that involves a highly conserved aromatic residue playing a pivotal role in rotation of the 6th transmembrane helix.

Reviewers

This article was reviewed by Oded Beja, G. P. S. Raghava and L. Aravind.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-015-0091-4) contains supplementary material, which is available to authorized users.

An empirical test of convergent evolution in rhodopsins

Rhodopsins are photochemically reactive membrane proteins that covalently bind retinal chromophores. Type I rhodopsins are found in both prokaryotes and eukaryotic microbes, whereas type II rhodopsins function as photoactivated G-protein coupled receptors (GPCRs) in animal vision. Both rhodopsin families share the seven transmembrane α-helix GPCR fold and a Schiff base linkage from a conserved lysine to retinal in helix G. Nevertheless, rhodopsins are widely cited as a striking example of evolutionary convergence, largely because the two families lack detectable sequence similarity and differ in many structural and mechanistic details. Convergence entails that the shared rhodopsin fold is so especially suited to photosensitive function that proteins from separate origins were selected for this architecture twice. Here we show, however, that the rhodopsin fold is not required for photosensitive activity. We engineered functional bacteriorhodopsin variants with novel folds, including radical noncircular permutations of the α-helices, circular permutations of an eight-helix construct, and retinal linkages relocated to other helices. These results contradict a key prediction of convergence and thereby provide an experimental attack on one of the most intractable problems in molecular evolution: how to establish structural homology for proteins devoid of discernible sequence similarity.

Relocating the active-site lysine in rhodopsin and implications for evolution of retinylidene proteins

Type I and type II rhodopsins share several structural features including a G protein-coupled receptor fold and a highly conserved active-site Lys residue in the seventh transmembrane segment of the protein. However, the two families lack significant sequence similarity that would indicate common ancestry. Consequently, the rhodopsin fold and conserved Lys are widely thought to have arisen from functional constraints during convergent evolution. To test for the existence of such a constraint, we asked whether it were possible to relocate the highly conserved Lys296 in the visual pigment bovine rhodopsin. We show here that the Lys can be moved to three other locations in the protein while maintaining the ability to form a pigment with 11-cis-retinal and activate the G protein transducin in a light-dependent manner. These results contradict the convergent hypothesis and support the homology of type I and type II rhodopsins by divergent evolution from a common ancestral protein.

Structure determination of membrane proteins by nuclear magnetic resonance spectroscopy

Many biological membranes consist of 50% or more (by weight) membrane proteins, which constitute approximately one-third of all proteins expressed in biological organisms. Helical membrane proteins function as receptors, enzymes, and transporters, among other unique cellular roles. Additionally, most drugs have membrane proteins as their receptors, notably the superfamily of G protein-coupled receptors with seven transmembrane helices. Determining the structures of membrane proteins is a daunting task because of the effects of the membrane environment; specifically, it has been difficult to combine biologically compatible environments with the requirements for the established methods of structure determination. There is strong motivation to determine the structures in their native phospholipid bilayer environment so that perturbations from nonnatural lipids and phases do not have to be taken into account. At present, the only method that can work with proteins in liquid crystalline phospholipid bilayers is solid-state NMR spectroscopy.

Evolution of GPCR: change and continuity

Once introduced into the very early eukaryotic blueprint, seven-transmembrane receptors soon became the central and versatile components of the evolutionary highly successful G protein-coupled transmembrane signaling mechanism. In contrast to all other components of this signal transduction pathway, G protein-coupled receptors (GPCR) evolved in various structural families, eventually comprising hundreds of members in vertebrate genomes. Their functional diversity is in contrast to the conserved transmembrane core and the invariant set of intracellular signaling mechanisms, and it may be the interplay of these properties that is the key to the evolutionary success of GPCR. The GPCR repertoires retrieved from extant vertebrate genomes are the recent endpoints of this long evolutionary process. But the shaping of the fine structure and the repertoire of GPCR is still ongoing, and signatures of recent selection acting on GPCR genes can be made visible by modern population genetic methods. The very dynamic evolution of GPCR can be analyzed from different perspectives: at the levels of sequence comparisons between species from different families, orders and classes, and at the level of populations within a species. Here, we summarize the main conclusions from studies at these different levels with a specific focus on the more recent evolutionary dynamics of GPCR.

The predictability of evolution: glimpses into a post-Darwinian world

The very success of the Darwinian explanation, in not only demonstrating evolution from multiple lines of evidence but also in providing some plausible explanations, paradoxically seems to have served to have stifled explorations into other areas of investigation. The fact of evolution is now almost universally yoked to the assumption that its outcomes are random, trends are little more than drunkard's walks, and most evolutionary products are masterpieces of improvisation and far from perfect. But is this correct? Let us consider some alternatives. Is there evidence that evolution could in anyway be predictable? Can we identify alternative forms of biological organizations and if so how viable are they? Why are some molecules so extraordinarily versatile, while others can be spoken of as "molecules of choice"? How fortuitous are the major transitions in the history of life? What implications might this have for the Tree of Life? To what extent is evolutionary diversification constrained or facilitated by prior states? Are evolutionary outcomes merely sufficient or alternatively are they highly efficient, even superb? Here I argue that in sharp contradistinction to an orthodox Darwinian view, not only is evolution much more predictable than generally assumed but also investigation of its organizational substrates, including those of sensory systems, which indicates that it is possible to identify a predictability to the process and outcomes of evolution. If correct, the implications may be of some significance, not least in separating the unexceptional Darwinian mechanisms from underlying organizational principles, which may indicate evolutionary inevitabilities.

On the origins of arrestin and rhodopsin

Background

G protein coupled receptors (GPCRs) are the most numerous proteins in mammalian genomes, and the most common targets of clinical drugs. However, their evolution remains enigmatic. GPCRs are intimately associated with trimeric G proteins, G protein receptor kinases, and arrestins. We conducted phylogenetic studies to reconstruct the history of arrestins. Those findings, in turn, led us to investigate the origin of the photosensory GPCR rhodopsin.

Results

We found that the arrestin clan is comprised of the Spo0M protein family in archaea and bacteria, and the arrestin and Vps26 families in eukaryotes. The previously known animal arrestins are members of the visual/beta subfamily, which branched from the founding "alpha" arrestins relatively recently. Curiously, we identified both the oldest visual/beta arrestin and opsin genes in Cnidaria (but not in sponges). The arrestin clan has 14 human members: 6 alphas, 4 visual/betas, and 4 Vps26 genes. Others recently showed that the 3D structure of mammalian Vps26 and the biochemical function of the yeast alpha arrestin PalF are similar to those of beta arrestins. We note that only alpha arrestins have PY motifs (known to bind WW domains) in their C-terminal tails, and only visual/betas have helix I in the Arrestin N domain.

Conclusion

We identified ciliary opsins in Cnidaria and propose this subfamily is ancestral to all previously known animal opsins. That finding is consistent with Darwin's theory that eyes evolved once, and lends some support to Parker's hypothesis that vision triggered the Cambrian explosion of life forms. Our arrestin findings have implications on the evolution of GPCR signaling, and on the biological roles of human alpha arrestins.

The tachykinin tale: molecular recognition in a historical perspective

Crystallography, mutational mapping and crosslinking are but a few of the experimental techniques that have helped to elucidate the underlying principles of molecular recognition between macromolecules and to improve our understanding of the evolution of the structure-activity relationship (SAR). While this development has been particularly successful for small and rigid ligands and substrates that bind to larger hydrophilic biomolecules, our understanding of membrane-embedded proteins is still rather limited. This review uses the example of the neuropeptide family of tachykinins and their G-protein coupled receptors (GPCR) to present how complementary experimental strategies over the past decades have nourished and modified conceptual models of the structural requisites of molecular recognition and function. Given the little we know, the pertinent question is how we proceed from here.

Three-dimensional models for beta-adrenergic receptor complexes with agonists and antagonists

Molecular modeling methods have been used to construct three-dimensional models for agonist and antagonist complexes with beta-adrenergic receptors. The recent rhodopsin crystal structure was used as a template in standard homology modeling methods. The rhodopsin-based homology models were assessed for agreement with experimental results for beta-adrenergic receptors, and compared with receptor models developed using de novo modeling techniques. While the de novo and homology-derived receptor models are generally quite similar, there are some localized structural differences that impact the putative ligand-binding site significantly. The de novo receptor models appear to provide much better agreement with experimental data, particularly for receptor complexes with agonist ligands. The de novo receptor models also yield some interesting and testable hypotheses for the structural basis of beta-adrenergic receptor subtype ligand selectivity.

New hypotheses for the GPCR 3D arrangement based on a molecular model of the human sweet-taste receptor

A molecular model of the human sweet-taste receptor has been inferred from superpositions of 3D maps of sweetener interaction sites (themselves previously deduced from extensive structure-activity relationship studies on highly potent sweeteners) onto three well-known G protein-coupled receptors (GPCRs)-rhodopsin, beta(2)- and alpha(2A)-adrenergic receptors-assumed to be linked by common evolutionary origins. The model gives new answers to old questions on the GPCR 3D structure, such as on the orientation and arrangement of the binding helices, their interaxial distances, radial orientations and relative heights. The model should be useful as a new approach to the rational design of drugs.

Evidence for kinship between diverse G-protein coupled receptors

In an earlier publication we described similarities at the primary sequence level between the first probable plant G-protein coupled receptor (GPCR) and three GPCR families (families A, B and F according to Kolakowski's classification) that were previously considered evolutionarily unrelated. Here we analyze further the relatedness among different GPCR families. By using PSI-BLAST, which is a search algorithm that is more potent in detecting weak similarities, one finds additional similarities between GPCR families that have not previously been described. Based on these comparisons, it is possible to divide all the GPCR families into one large clade and two smaller ones. The large clade includes the rhodopsin family (family A), the glucagon receptor family (family B), cyclic AMP receptors (family F), an Arabidopsis thaliana receptor, the Frizzled family and probably also the STE3 pheromone receptors (family E) and vomeronasal receptors type 1. The smaller clades consist of, in one case, BOSS and the GABA-B family of receptors (family C), and in the other the STE2 pheromone receptors (family D) alone. Although our findings are likely to reflect a common ancestry within each of these clades, whether or not two or all three of the clades also share an even more ancient ancestor between them remains an open question that cannot be answered from our present data.

Structural organization of G-protein-coupled receptors

Atomic-resolution structures of the transmembrane 7-alpha-helical domains of 26 G-protein-coupled receptors (GPCRs) (including opsins, cationic amine, melatonin, purine, chemokine, opioid, and glycoprotein hormone receptors and two related proteins, retinochrome and Duffy erythrocyte antigen) were calculated by distance geometry using interhelical hydrogen bonds formed by various proteins from the family and collectively applied as distance constraints, as described previously [Pogozheva et al., Biophys. J., 70 (1997) 1963]. The main structural features of the calculated GPCR models are described and illustrated by examples. Some of the features reflect physical interactions that are responsible for the structural stability of the transmembrane alpha-bundle: the formation of extensive networks of interhelical H-bonds and sulfur-aromatic clusters that are spatially organized as 'polarity gradients'; the close packing of side-chains throughout the transmembrane domain; and the formation of interhelical disulfide bonds in some receptors and a plausible Zn2+ binding center in retinochrome. Other features of the models are related to biological function and evolution of GPCRs: the formation of a common 'minicore' of 43 evolutionarily conserved residues; a multitude of correlated replacements throughout the transmembrane domain; an Na(+)-binding site in some receptors, and excellent complementarity of receptor binding pockets to many structurally dissimilar, conformationally constrained ligands, such as retinal, cyclic opioid peptides, and cationic amine ligands. The calculated models are in good agreement with numerous experimental data.

Evolution of the archaeal rhodopsins: evolution rate changes by gene duplication and functional differentiation

The amino acid sequences of 25 archaeal retinal proteins from 13 different strains of extreme halophiles were analyzed to establish their molecular phylogenetic relationship. On the basis of amino acid sequence similarity, these proteins apparently formed a distinct family designated as the archaeal rhodopsin family (ARF), which was not related to other known proteins, including G protein-coupled receptors. The archaeal rhodopsin family was further divided into four clusters with different functions; H+ pump (bacteriorhodopsin), Cl- pump (halorhodopsin), and two kinds of sensor (sensory rhodopsin and phoborhodopsin). These four rhodopsin clusters seemed to have occurred by gene duplication(s) before the generic speciation of halophilic archaea, based on phylogenetic analysis. Therefore, the degrees of differences in amino acid sequences within each cluster simply reflected the divergent evolution of halophilic archaea. By comparing the branch lengths after speciation points of the reconstituted tree, we calculated the relative evolution rates of the four archaeal rhodopsins bacteriorhodopsin:halorhodopsin:sensory rhodopsin: phoborhodopsin to be 5:4:3:10. From these values, the degrees of functional and structural restriction of each protein can be inferred. The branching topology of four clusters grouped bacteriorhodopsin and halorhodopsin versus sensory rhodopsin and phoborhodopsin by likelihood mapping. Using bacteriorhodopsin (and halorhodopsin) as an outgroup, the gene duplication point of sensory rhodopsin/phoborhodopsin was determined. By calculating the branch lengths between the gene duplication point and each halophilic archaea speciation point, we could speculate upon the relative evolution rate of pre-sensory rhodopsin and pre-phoborhodopsin. The evolution rate of pre-sensory rhodopsin was fivefold faster than that of pre-phoborhodopsin, which suggests that the original function of the ancestral sensor was similar to that of phoborhodopsin, and that sensory rhodopsin evolved from pre-sensory rhodopsin by the accumulation of mutations. The changes in evolution rate by gene duplication and functional differentiation were demonstrated in the archaeal rhodopsin family using the gene duplication date and halobacterial speciation date as common time stamps.

Towards structural models of molecular recognition in olfactory receptors

The G protein coupled receptors (GPCR) are an important class of proteins that act as signal transducers through the cytoplasmic membrane. Understanding the structure and activation mechanism of these proteins is crucial for understanding many different aspects of cellular signalling. The olfactory receptors correspond to the largest family of GPCRs. Very little is known about how the structures of the receptors govern the specificity of interaction which enables identification of particular odorant molecules. In this paper, we review recent developments in two areas of molecular modelling: methods for modelling the configuration of trans-membrane helices and methods for automatic docking of ligands into receptor structures. We then show how a subset of these methods can be combined to construct a model of a rat odorant receptor interacting with lyral for which experimental data are available. This modelling can help us make progress towards elucidating the specificity of interactions between receptors and odorant molecules.

Evolutionary relationships among proteins probed by an iterative neighborhood cluster analysis (INCA). Alignment of bacteriorhodopsins with the yeast sequence YRO2

PURPOSE

Searching the existing databases for homologous sequences is essential to understanding a protein's structure and function. For a query sequence, its nearest neighbors can be identified by BLAST (basic local alignment search tool). However, a single query sequence is sufficient to define the entire neighborhood of related sequences, and multiple BLAST queries are needed. We describe here a program which permits automated and iterative BLAST analysis of an entire neighborhood of sequences and apply this to search for homologs of the bacteriorhodopsins outside the archaea phylum.

METHODS

We have developed a Java program, 'Iterative Neighborhood Cluster Analysis' (INCA), which performs iterative BLAST searches, beginning with a single starter sequence, and proceeding with any other sequence achieving a predefined minimum alignment score. This results in a cluster of sequences where each sequence is related to at least one other sequence by the cutoff score, additional lists of more distantly related sequences for each member of cluster.

RESULTS

Bacteriorhodopsins had not been previously aligned with any other protein family with scores indicative of probable homology. Using INCA, we identified a probable homolog in yeast, YRO2_YEAST, also containing seven putative transmembrane domains. A finding of probable homology was supported by additional alignment strategies.

CONCLUSIONS

INCA is a useful tool to assess complete protein neighborhoods. With an increasing database, INCA can serve to detect the emergence of evolutionary links between even the most distantly related protein families. Identifying a homolog of the bacteriorhodopsins in yeast illustrates this approach but at the same time highlights the vast evolutionary distances between polytopic membrane proteins, such as the bacteriorhodopsins.

The transmembrane 7-alpha-bundle of rhodopsin: distance geometry calculations with hydrogen bonding constraints

A 3D model of the transmembrane 7-alpha-bundle of rhodopsin-like G-protein-coupled receptors (GPCRs) was calculated using an iterative distance geometry refinement with an evolving system of hydrogen bonds, formed by intramembrane polar side chains in various proteins of the family and collectively applied as distance constraints. The alpha-bundle structure thus obtained provides H bonding of nearly all buried polar side chains simultaneously in the 410 GPCRs considered. Forty evolutionarily conserved GPCR residues form a single continuous domain, with an aliphatic "core" surrounded by six clusters of polar and aromatic side chains. The 7-alpha-bundle of a specific GPCR can be calculated using its own set of H bonds as distance constraints and the common "average" model to restrain positions of the helices. The bovine rhodopsin model thus determined is closely packed, but has a few small polar cavities, presumably filled by water, and has a binding pocket that is complementary to 11-cis (6-s-cis, 12-s-trans, C = N anti)-retinal or to all-trans-retinal, depending on conformations of the Lys296 and Trp265 side chains. A suggested mechanism of rhodopsin photoactivation, triggered by the cis-trans isomerization of retinal, involves rotations of Glu134, Tyr223, Trp265, Lys296, and Tyr306 side chains and rearrangement of their H bonds. The model is in agreement with published electron cryomicroscopy, mutagenesis, chemical modification, cross-linking, Fourier transform infrared spectroscopy, Raman spectroscopy, electron paramagnetic resonance spectroscopy, NMR, and optical spectroscopy data. The rhodopsin model and the published structure of bacteriorhodopsin have very similar retinal-binding pockets.

Three-dimensional models for agonist and antagonist complexes with beta 2 adrenergic receptor

Computer-modeling techniques have been used to generate docked complexes for a series of beta adrenergic agonists and antagonists with a three-dimensional model of the beta 2 adrenergic receptor. For all ligands tested, it proved possible to dock low-energy conformers in the receptor model, with sensible electrostatic, steric, and hydrogen-bonding interactions, many of which are supported by experimental studies of the beta 2 receptor. Our results illustrate the power of molecular modeling techniques, when coupled with appropriate experimental methods and data, to investigate structure-function properties of integral membrane receptor proteins that cannot yet be studied by direct structural methods.

Involvement of Asn-293 in stereospecific agonist recognition and in activation of the beta 2-adrenergic receptor

To investigate the molecular mechanism for stereospecific binding of agonists to beta 2-adrenergic receptors we used receptor models to identify potential binding sites for the beta-OH-group of the ligand, which defines the chiral center. Ser-165, located in transmembrane helix IV, and Asn-293, situated in the upper half of transmembrane helix VI, were identified as potential binding sites. Mutation of Ser-165 to Ala did not change the binding of either isoproterenol isomer as revealed after transient expression in human embryonic kidney (HEK)-293 cells. In contrast, a receptor mutant in which Asn-293 was replaced by Leu showed substantial loss of stereospecific isoproterenol binding. Adenylyl cyclase stimulation by this mutant after stable expression in CHO cells confirmed the substantial loss of stereospecificity for isoproterenol. In a series of agonists the loss of affinity in the Leu-293 mutant receptor was strongly correlated with the intrinsic activity of the compounds. Full agonists showed a 10-30-fold affinity loss, whereas partial agonists had almost the same affinity for both receptors. Stereospecific recognition of antagonists was unaltered in the Leu-293 mutant receptor. These data indicate a relationship between stereospecificity and intrinsic activity of agonists and suggest that Asn-293 is important for both properties of the agonist-receptor interaction.

An analysis of the conserved residues between halobacterial retinal proteins and G-protein coupled receptors: implications for GPCR modeling

An alignment of the transmembrane domains of halobacterial retinal proteins (including bacteriorhodopsin) and G-protein coupled receptors (GPCRs) is presented based on the commonality of conserved residues between families. Due to the limited sequence homology displayed by these proteins, an alternative strategy is proposed for sequence alignment that correlates residues within secondary structure elements. The nonsequential alignment developed identifies three proline and two aspartates residues that share common positions and, in the former case, similar functions in the transmembrane domain. The alignment is further applied to model the packing of transmembrane helices 5 and 6 of the beta-adrenergic receptor based on the backbone coordinates of bacteriorhodopsin helices 3 and 2, respectively. Unlike models derived from standard sequential alignments, the approach developed here allows the key structural features conferred by the proline residues to be captured during model building. The structure described is also compared with available site directed mutagenesis results as well as existing GPCR models. In addition to the implications to model building, the commonality observed suggests a potential relationship among the GPCRs and retinal proteins.

Automated method for modeling seven-helix transmembrane receptors from experimental data

A rule-based automated method is presented for modeling the structures of the seven transmembrane helices of G-protein-coupled receptors. The structures are generated by using a simulated annealing Monte Carlo procedure that positions and orients rigid helices to satisfy structural restraints. The restraints are derived from analysis of experimental information from biophysical studies on native and mutant proteins, from analysis of the sequences of related proteins, and from theoretical considerations of protein structure. Calculations are presented for two systems. The method was validated through calculations using appropriate experimental information for bacteriorhodopsin, which produced a model structure with a root mean square (rms) deviation of 1.87 A from the structure determined by electron microscopy. Calculations are also presented using experimental and theoretical information available for bovine rhodopsin to assign the helices to a projection density map and to produce a model of bovine rhodopsin that can be used as a template for modeling other G-protein-coupled receptors.