Ab initio modeling and molecular dynamics simulation of the alpha 1b-adrenergic receptor activation

Disruption of TAB1/p38α interaction using a cell-permeable peptide limits myocardial ischemia/reperfusion injury

Targeting the adaptor protein (transforming growth factor-β (TGF-β)-activated protein kinase 1 (TAK1)-binding protein 1) (TAB1)-mediated non-canonical activation of p38α to limit ischemia/reperfusion (I/R) injury after an acute myocardial infarction seems to be attractive since TAB1/p38α interaction occurs specifically in very limited circumstances and possesses unique structural basis. However, so far no TAB1/p38α interaction inhibitor has been reported due to the limited knowledge about the interfaces. In this study, we sought to identify key amino acids essential for the unique mode of interaction with computer-guided molecular simulations and molecular docking. After validation of the predicted three-dimensional (3-D) structure of TAB1/p38α complex, we designed several peptides and evaluated whether they could block TAB1/p38α interaction with selectivity. We found that a cell-permeable peptide worked as a selective TAB1/p38α interaction inhibitor and decreased myocardial I/R injury. To our knowledge, this is the first TAB1/p38α interaction inhibitor.

GPCR Conformations: Implications for Rational Drug Design

G protein-coupled receptors (GPCRs) comprise a large class of transmembrane proteins that play critical roles in both normal physiology and pathophysiology. These critical roles offer targets for therapeutic intervention, as exemplified by the substantial fraction of current pharmaceutical agents that target members of this family. Tremendous contributions to our understanding of GPCR structure and dynamics have come from both indirect and direct structural characterization techniques. Key features of GPCR conformations derived from both types of characterization techniques are reviewed.

Constitutive activity and inverse agonism at the α(₁a) and α(₁b) adrenergic receptor subtypes

The α(1b)-adrenergic receptor (AR) was, after rhodopsin, the first G protein-coupled receptor (GPCR) in which point mutations were shown to trigger constitutive (agonist-independent) activity. Constitutively activating mutations have been found in other AR subtypes as well as in several GPCRs. This chapter briefly summarizes the main findings on constitutively active mutants of the α(1a)- and α(1b)-AR subtypes and the methods used to predict activating mutations, to measure constitutive activity of Gq-coupled receptors and to investigate inverse agonism. In addition, it highlights the implications of studies on constitutively active AR mutants on elucidating the molecular mechanisms of receptor activation and drug action.

Glutamic acid decarboxylase-derived epitopes with specific domains expand CD4(+)CD25(+) regulatory T cells

Background

CD4⁺CD25⁺ regulatory T cell (Treg)-based immunotherapy is considered a promising regimen for controlling the progression of autoimmune diabetes. In this study, we tested the hypothesis that the therapeutic effects of Tregs in response to the antigenic epitope stimulation depend on the structural properties of the epitopes used.

Methodology/Principal Findings

Splenic lymphocytes from nonobese diabetic (NOD) mice were stimulated with different glutamic acid decarboxylase (GAD)-derived epitopes for 7–10 days and the frequency and function of Tregs was analyzed. We found that, although all expanded Tregs showed suppressive functions in vitro, only p524 (GAD524–538)-expanded CD4⁺CD25⁺ T cells inhibited diabetes development in the co-transfer models, while p509 (GAD509–528)- or p530 (GAD530–543)-expanded CD4⁺CD25⁺ T cells had no such effects. Using computer-guided molecular modeling and docking methods, the differences in structural characteristics of these epitopes and the interaction mode (including binding energy and identified domains in the epitopes) between the above-mentioned epitopes and MHC class II I-A^g7 were analyzed. The theoretical results showed that the epitope p524, which induced protective Tregs, possessed negative surface-electrostatic potential and bound two chains of MHC class II I-A^g7, while the epitopes p509 and p530 which had no such ability exhibited positive surface-electrostatic potential and bound one chain of I-A^g7. Furthermore, p524 bound to I-A^g7 more stably than p509 and p530. Of importance, we hypothesized and subsequently confirmed experimentally that the epitope (GAD570–585, p570), which displayed similar characteristics to p524, was a protective epitope by showing that p570-expanded CD4⁺CD25⁺ T cells suppressed the onset of diabetes in NOD mice.

Conclusions/Significance

These data suggest that molecular modeling-based structural analysis of epitopes may be an instrumental tool for prediction of protective epitopes to expand functional Tregs.

Computational study of the heterodimerization between mu and delta receptors

A growing body of evidence indicated that the G protein coupled receptors exist as homo- or hetero-dimers in the living cell. The heterodimerization between mu and delta opioid receptors has attracted researchers' particular interests, it is reported to display novel pharmacological and signalling regulation properties. In this study, we construct the full-length 3D-model of mu and delta opioid receptors using the homology modelling method. Threading program was used to predict the possible templates for the N- and C-terminus domains. Then, a 30 ns molecular dynamics simulations was performed with each receptor embedded in an explicit membrane-water environment to refine and explore the conformational space. Based on the structures extracted from the molecular dynamics, the likely interface of mu-delta heterodimer was investigated through the analysis of protein-protein docking, cluster, shape complementary and interaction energy. The computational modelling works revealed that the most likely interface of heterodimer was formed between the transmembrane1,7 (TM1,7) domains of mu receptor and the TM(4,5) domains of delta receptor, with emphasis on mu-TM1 and delta-TM4, the next likely interface was mu(TM6,7)-delta(TM4,5), with emphasis on mu-TM6 and delta-TM4. Our results were consistent with previous reports.

Inactive and active states and supramolecular organization of GPCRs: insights from computational modeling

Herein we make an overview of the results of our computational experiments aimed at gaining insight into the molecular mechanisms of GPCR functioning either in their normal conditions or when hit by gain-of-function or loss-of-function mutations. Molecular simulations of a number of GPCRs in their wild type and mutated as well as free and ligand-bound forms were instrumental in inferring the structural features, which differentiate the mutation- and ligand-induced active from the inactive states. These features essentially reside in the interaction pattern of the E/DRY arginine and in the degree of solvent exposure of selected cytosolic domains. Indeed, the active states differ from the inactive ones in the weakening of the interactions made by the highly conserved arginine and in the increase in solvent accessibility of the cytosolic interface between helices 3 and 6. Where possible, the structural hallmarks of the active and inactive receptor states are translated into molecular descriptors useful for in silico functional screening of novel receptor mutants or ligands. Computational modeling of the supramolecular organization of GPCRs and their intracellular partners is the current challenge toward a deep understanding of their functioning mechanisms.

A two model receptor system of the alpha1D adrenergic receptor to describe interactions with epinephrine and BMY7378

In this study, we have developed a two receptor model system to describe the R and R states of G-protein coupled receptors, specifically the alpha(1D) adrenergic receptor. The two models interact with agonist (epinephrine) and antagonist (BMY7378) differently. The active model has increased interactions with epinephrine. The inactive model has increased interactions with BMY7378. We also explored the protonation state of the ligands. When the most basic amine was protonated, we found increased hydrogen bonding and increased aromatic interactions. Protonated epinephrine hydrogen bonds with Asp176 and has aromatic residues Trp172, Trp235, Trp361, and Phe388 within 3 Angstroms. Protonated BMY7378 hydrogen bonds with Trp172 and Lys236 and has aromatic residues Trp172, Trp254, Phe364, Phe384, and Phe388 within 3 Angstroms. We conclude that the two model system is required to represent the two states of the receptor and that protonation of the ligand is also critical.

In vivo identification of the interaction site of ErbB2 extracellular domain with its autoinhibitor

Direct interference with the transforming potential of ErbB2 has become a subject of great interest. Disruption of critical ErbB2 ectodomain interactions may lead to novel therapeutic approaches for the treatment of various tumors. The ErbB receptor signaling can be inhibited by rationally designed peptide mimetics based on the subdomains of ErbB ectodomain. The mimetics can bind to the ErbB receptor specifically and block inter-receptor interactions, resulting in the growth inhibition of ErbB2-overexpressing cells in vitro. In this study, three-dimensional structure of herstatin, an autoinhibitor of ErbB2 and ErbB2 ectodomain complex was constructed by computer-aided molecular modeling. The binding site on ErbB2 ectodomain for herstatin was determined at S1 domain. The mutants of ErbB2 ectodomain were constructed. The interactions of ErbB2 ectodomain and its mutants with herstatin were analyzed for the first time in living cells that coexpressed herstatin and ErbB2 ectodomain or the mutants. The S1 domain in ErbB2 ectodomain was verified as the interaction site with herstatin by immunoprecipitation, confocal microscopy, and fluorescence resonance energy transfer (FRET). The binding region of herstatin on ErbB2 ectodomain might be a potential target region for the drug design.

Atypical and typical antipsychotic drug interactions with the dopamine D2 receptor

A model of the dopamine D2 receptor was used to study the receptor interactions of dopamine, the typical antipsychotics haloperidol and loxapine, and the atypical antipsychotics clozapine and melperone. The atypical antipsychotics interacted with the halogen atom of the ring system in the direction of the transmembrane helices (TMHs) 2, 3 and 7, while the typical had the corresponding halogen atom in the direction of TMH5. Molecular dynamics simulations indicated that the average helical displacement upon binding increased in the order: typical < atypical < dopamine. Upon binding, the atypical induced larger displacements into TMH5 than did the typical. The typical had stronger non-bonded interactions with the receptor than had the atypical, which is in agreement with the experimental observation that the atypical antipsychotic drugs dissociate faster from the receptor than the typical antipsychotic drugs.

Molecular model of the neural dopamine transporter

The dopamine transporter (DAT) regulates the action of dopamine by reuptake of the neurotransmitter into presynaptic neurons, and is the main molecular target of amphetamines and cocaine. DAT and the Na+/H+ antiporter (NhaA) are secondary transporter proteins that carry small molecules across a cell membrane against a concentration gradient, using ion gradients as energy source. A 3-dimensional projection map of the E. coli NhaA has confirmed a topology of 12 membrane spanning domains, and was previously used to construct a 3-dimensional NhaA model with 12 trans-membrane alpha-helices (TMHs). The NhaA model, and site directed mutagenesis data on DAT, were used to construct a detailed 3-dimensional DAT model using interactive molecular graphics and empiric force field calculations. The model proposes a dopamine transport mechanism involving TMHs 1, 3, 4, 5, 7 and 11. Asp79, Tyr252 and Tyr274 were the primary cocaine binding residues. Binding of cocaine or its analogue, (-)-2beta-carbomethoxy-3beta-(4-fluorophenyl)tropane (CFT), seemed to lock the transporter in an inactive state, and thus inhibit dopamine transport. The present model may be used to design further experimental studies of the molecular structure and mechanisms of DAT and other secondary transporter proteins.

Intrinsic activity and comparative molecular dynamics of buspirone analogues at the 5-HT1A receptors

In CNS, the 5-hydroxytryptamine(1A) (5-HT(1A)) receptors exist in two different populations with different behavioural and physiological effects: (1) somatodendritic autoreceptors located pre-synaptically of 5-HT containing neurons and (2) receptors located post-synaptic to 5-HT containing neurons. Clinical studies have shown that 5-HT(1A) partial agonists have anxiolytic properties, while antagonists of pre-synaptical autoreceptors shorten the onset time of selective serotonin reuptake inhibitors (SSRIs). In the present study, the pre- and post-synaptic activity of structural analogues of buspirone was evaluated in animal models. A three dimensional model of the 5-HT(1A) receptor was used to study their interaction modes and helical displacements upon receptor binding. The predicted receptor-ligand interactions indicated similarities in the receptor binding modes for all buspirone analogues, and no clear relationship between receptor contact residues and activity at pre- and post-synaptic receptors. Comparative molecular dynamics (MD) simulations for 650ps indicated that pre-synaptic antagonistic behaviour is connected to large displacements of transmembrane helix (TMH) 7 upon binding, while pre-synaptic agonistic behaviour is connected to large displacements of TMH2 and small displacements of TMH7. Post-synaptic partial agonist behaviour is connected to large displacements of TMH4 and TMH5 upon binding, while post-synaptic antagonists only slightly displace these helices.

Structural features of the inactive and active states of the melanin-concentrating hormone receptors: Insights from molecular simulations

Comparative molecular dynamics simulations of both subtypes 1 and 2 of the melanin-concentrating hormone receptor (MCHR1 and MCHR2, respectively) in their free and hormone-bound forms have been carried out. The hormone has been used in its full-length and truncated forms, as well as in 16 mutated forms. Moreover, MCHR1 has been simulated in complex with T-226296, a novel orally active and selective antagonist. The comparative analysis of an extended number of receptor configurations suggests that the differences between inactive (i.e., free and antagonist-bound) and active (i.e., agonist-bound) states of MCHRs involve the receptor portions close to the E/DRY and NPxxY motifs, with prominence to the cytosolic extensions of helices 2, 3, 6, and 7. In fact, the active forms of these receptors share the release of selected intramolecular interactions found in the inactive forms, such as that between R3.50 of the E/DRY motif and D2.40, and that between Y7.53 of the NPxxY motif and F7.60. Another feature of the active forms of both MCHRs is the approach of "helix 8" to the cytosolic extension of helix 3. These features of the active forms are concurrent with the opening of a cleft at the cytosolic end of the helix bundle. For both MCHRs, the agonist-induced chemical information transfer from the extracellular to the cytosolic domains is mediated by a cluster of aromatic amino acids in helix 6, following the ligand interaction with selected amino acids in the extracellular half of the receptor.

G Protein–Coupled Receptor Oligomerization

The cardiovascular system is richly endowed with G protein–coupled receptors (GPCRs), members of the largest family of plasma membrane-localized receptors. During the last 10 years, it has become increasingly clear that many, if not all, GPCRs function in oligomeric complexes, as either homo- or hetero-oligomers. This review explores the mechanistic implications of GPCR dimerization and/or oligomerization on receptor activation and interactions with G proteins. The effects of GPCR oligomerization on receptor pharmacology, GPCR-mediated signaling, and potential contributions to GPCR crosstalk will be considered in the context of receptors important in the cardiovascular system. Our evolving understanding of the structural and functional consequences of GPCR oligomerization may provide novel and more selective sites for pharmacological tuning of cardiovascular function.

Sequence analysis reveals how G protein-coupled receptors transduce the signal to the G protein

Sequence entropy-variability plots based on alignments of very large numbers of sequences-can indicate the location in proteins of the main active site and modulator sites. In the previous article in this issue, we applied this observation to a series of well-studied proteins and concluded that it was possible to detect most of the residues with a known functional role. Here, we apply the method to rhodopsin-like G protein-coupled receptors. Our conclusion is that G protein binding is the main evolutionary constraint on these receptors, and that other ligands, such as agonists, act as modulators. The activation of the receptors can be described as a simple, two-step process, and the residues involved in signal transduction can be identified.

Molecular Dynamics Simulations of the Ligand-Induced Chemical Information Transfer in the 5-HT1A Receptor

Comparative molecular dynamics simulations of the 5-HT(1A) receptor in its empty as well as agonist- (i.e. active) and antagonist-bound (i.e. nonactive) forms have been carried out. The agonists 5-HT and (R)-8-OH-DPAT as well as the antagonist WAY100635 have been employed. The results of this study strengthen the hypothesis that the receptor portions close to the E/DRY/W motif, with prominence to the cytosolic extensions of helices 3 and 6, are particularly susceptible to undergo structural modification in response to agonist binding. Despite the differences in the structural/dynamics behavior of the two agonists when docked into the 5-HT(1A) receptor, they both exert a destabilization of the intrahelical and interhelical interactions found in the empty and antagonist-bound receptor forms between the arginine of the E/DRY sequence and both D133(3.49) and E340(6.30). For both agonists, the chemical information transfer from the extracellular to the cytosolic domains is mediated by a cluster of aromatic amino acids in helix 6, following the ligand interaction with selected amino acids in the extracellular half of the receptor, such as D116(3.32), S199(5.42), Y195(5.38), and F361(6.51). A significant reduction in the bend at P360(6.50), as compared to the empty and the antagonist-bound receptor forms, is one of the features of the agonist-bound forms that is related to the breakage of the interhelical salt bridge between the E/DRY arginine and E340(6.30). Another structural feature, shared by the agonist-bound receptor forms and not by the empty and antagonist-bound forms, is the detachment of helices 2 and 4, as marked by the movement of W161(4.50) away from helix 2, toward the membrane space.

Mutagenesis and modelling of the alpha(1b)-adrenergic receptor highlight the role of the helix 3/helix 6 interface in receptor activation

Computer simulations on a new model of the alpha1b-adrenergic receptor based on the crystal structure of rhodopsin have been combined with experimental mutagenesis to investigate the role of residues in the cytosolic half of helix 6 in receptor activation. Our results support the hypothesis that a salt bridge between the highly conserved arginine (R143(3.50)) of the E/DRY motif of helix 3 and a conserved glutamate (E289(6.30)) on helix 6 constrains the alpha1b-AR in the inactive state. In fact, mutations of E289(6.30) that weakened the R143(3.50)-E289(6.30) interaction constitutively activated the receptor. The functional effect of mutating other amino acids on helix 6 (F286(6.27), A292(6.33), L296(6.37), V299(6.40,) V300(6.41), and F303(6.44)) correlates with the extent of their interaction with helix 3 and in particular with R143(3.50) of the E/DRY sequence.

Interaction of 1,2,4-substituted piperazines, new serotonin receptor ligands, with 5-HT1A and 5-HT2A receptors

In the present paper, we describe affinities to 5-HT1A and 5-HT2A receptors of several new 1,2,4-trisubstituted piperazine derivatives. The affinities were compared with those described earlier for 1,4-disubstituted piperazines and the influence of the third (methyl) substituent on the affinity to both receptors is discussed. The difference between two- and three-substituted derivatives was rationalised in terms of molecular modelling of the respective ligand-receptor complexes. Additionally, the functional activity of some 1,2,4-trisubstituted piperazines for 5-HT1A receptor was examined in behavioural and biochemical models. The obtained results have shown that some trisubstituted compounds exhibited a higher affinity to 5-HT2A receptors than their respective disubstituted analogues (with the affinity to 5-HT1A receptors remaining the same or somewhat improving). The molecular dynamics simulations suggested that the presence of the third substituent in the piperazine ring of those compounds may induce stabilising effect on the ligand-receptor complexes. The results of the in vivo studies have shown that some of the examined trisubstituted piperazines (10-13, 16, 17) exhibited properties of postsynaptic 5-HT1A partial agonists. Moreover, compounds 13 and 16 exhibited features of 5-HT1A presynaptic agonists in in vitro test, and compound 16 also in in vivo tests.

The alpha1b-adrenergic receptor subtype: molecular properties and physiological implications

The aim of this review is to summarize some of the main findings from our laboratory as well as from others concerning the biochemical, molecular, and functional properties of the alpha1b-adrenergic receptor. Experimental and computational mutagenesis of the alpha1b-adrenergic receptor have been instrumental in elucidating some of the molecular mechanisms underlying receptor activation and receptor coupling to Gq. The knockout mouse model lacking the alpha1b-adrenergic receptor has highlighted the potential implication of this receptor subtype in variety of functions including the regulation of blood pressure, glucose homeostasis, and the rewarding response to drugs of abuse.

Mutational and computational analysis of the alpha(1b)-adrenergic receptor. Involvement of basic and hydrophobic residues in receptor activation and G protein coupling

To investigate their role in receptor coupling to G(q), we mutated all basic amino acids and some conserved hydrophobic residues of the cytosolic surface of the alpha(1b)-adrenergic receptor (AR). The wild type and mutated receptors were expressed in COS-7 cells and characterized for their ligand binding properties and ability to increase inositol phosphate accumulation. The experimental results have been interpreted in the context of both an ab initio model of the alpha(1b)-AR and of a new homology model built on the recently solved crystal structure of rhodopsin. Among the twenty-three basic amino acids mutated only mutations of three, Arg(254) and Lys(258) in the third intracellular loop and Lys(291) at the cytosolic extension of helix 6, markedly impaired the receptor-mediated inositol phosphate production. Additionally, mutations of two conserved hydrophobic residues, Val(147) and Leu(151) in the second intracellular loop had significant effects on receptor function. The functional analysis of the receptor mutants in conjunction with the predictions of molecular modeling supports the hypothesis that Arg(254), Lys(258), as well as Leu(151) are directly involved in receptor-G protein interaction and/or receptor-mediated activation of the G protein. In contrast, the residues belonging to the cytosolic extensions of helices 3 and 6 play a predominant role in the activation process of the alpha(1b)-AR. These findings contribute to the delineation of the molecular determinants of the alpha(1b)-AR/G(q) interface.

GABA(B2) is essential for g-protein coupling of the GABA(B) receptor heterodimer

GABA(B) receptors are unique among G-protein-coupled receptors (GPCRs) in their requirement for heterodimerization between two homologous subunits, GABA(B1) and GABA(B2), for functional expression. Whereas GABA(B1) is capable of binding receptor agonists and antagonists, the role of each GABA(B) subunit in receptor signaling is unknown. Here we identified amino acid residues within the second intracellular domain of GABA(B2) that are critical for the coupling of GABA(B) receptor heterodimers to their downstream effector systems. Our results provide strong evidence for a functional role of the GABA(B2) subunit in G-protein coupling of the GABA(B) receptor heterodimer. In addition, they provide evidence for a novel "sequential" GPCR signaling mechanism in which ligand binding to one heterodimer subunit can induce signal transduction through the second partner of a heteromeric complex.

Molecular dynamics of 5-HT1A and 5-HT2A serotonin receptors with methylated buspirone analogues

In the present study experimentally determined ligand selectivity of three methylated buspirone analogues (denoted as MM2, MM5 and P55) towards 5-HT1A and 5-HT2A serotonin receptors was theoretically investigated on a molecular level. The relationships between the ligand structure and 5-HT1A and 5-HT2A receptor affinities were studied and the results were found to be in agreement with the available site-directed mutagenesis and binding affinity data. Molecular dynamics (MD) simulations of ligand-receptor complexes were performed for each investigated analogue, docked twice into the central cavity of 5-HT1A/5-HT2A, each time in a different orientation. Present results were compared with our previous theoretical results, obtained for buspirone and its non-methylated analogues. It was found that due to the presence of the methyl group in the piperazine ring the ligand position alters and the structure of the ligand-receptor complex is modified. Further, the positions of derivatives with pyrimidinyl aromatic moiety and quinolinyl moiety are significantly different at the 5-HT2A receptor. Thus, methylation of such derivatives alters the 3D structures of ligand-receptor complexes in different ways. The ligand-induced changes of the receptor structures were also analysed. The obtained results suggest, that helical domains of both receptors have different dynamical behaviour. Moreover, both location and topography of putative binding sites for buspirone analogues are different at 5-HT1A and 5-HT2A receptors.

Ligand induced conformational states of the 5-HT(1A) receptor

It has been shown that G-protein coupled receptors have seven transmembrane alpha-helices, but the structural changes occurring in a G-protein coupled receptor as a response on agonist stimulus and the molecular events leading to blockade of the signal transduction by antagonists are not well understood. In the present study, the AMBER 5.0 force field was used for comparative molecular dynamics simulations of a 5-HT(1A) receptor model in the absence of ligand, in complex with a 5-HT(1A) receptor agonist (R)-8-hydroxy-2-(di-n-propylamino)tetralin [(R)-8-OH-DPAT], in complex with a selective 5-HT(1A) receptor antagonist (S)-N-tert-butyl-3-[4-(2-methoxyphenyl)piperazin-1-yl ]-2-phenylpropanamide [(S)-WAY100135], and in complex with the partial agonist, buspirone. In the simulations, the agonist induced larger conformational changes into transmembrane helix 3 and 6 than into the other helices, while the main conformational differences between the agonist bound receptor and the antagonist bound receptor were in transmembrane helix 5 and 6. During the simulations, all the three ligands constrained the helical movements compared to those observed in the receptor without any ligand.

Switching agonist/antagonist properties of opiate alkaloids at the delta opioid receptor using mutations based on the structure of the orphanin FQ receptor

In an earlier study, we have demonstrated that by mutating five amino acid residues to those conserved in the opioid receptors, the OFQ receptor could be converted to a functional receptor that bound many opioid alkaloids with nanomolar affinities. Surprisingly, when the reciprocal mutations, Lys-214 --> Ala (TM5), Ile-277 --> Val/His-278 --> Gln/Ile-279 --> Val (TM6), and Ile-304 --> Thr (TM7), are introduced in the delta receptor, neither the individual mutations nor their various combinations significantly reduce the binding affinities of opioid alkaloids tested. However, these mutations cause profound alterations in the functional characteristics of the mutant receptors as measured in guanosine 5'-3-O-(thio)triphosphate binding assays. Some agonists become antagonists at some constructs as they lose their ability to activate them. Some alkaloid antagonists are transformed into agonists at other constructs, but their agonistic effects can still be blocked by the peptide antagonist TIPP. Even the delta inverse agonist 7-benzylidenenaltrexone becomes an agonist at the mutant containing both the Ile-277 --> Val/His-278 --> Gln/Ile-279 --> Val and Ile-304 --> Thr mutations. Thus, although the mutated residues are thought to be part of the binding pocket, they are critically involved in the control of the delta receptor activation process. These findings shed light on some of the structural bases of ligand efficacy. They are also compatible with the hypothesis that a ligand may achieve high affinity binding in several different ways, each having different effects on receptor activation.

Mutational analysis of the interaction of the N- and C-terminal ends of angiotensin II with the rat AT(1A) receptor

1. The role of different residues of the rat AT(1A) receptor in the interaction with the N- and C-terminal ends of angiotensin II (AngII) was studied by determining ligand binding and production of inositol phosphates (IP) in COS-7 cells transiently expressing the following AT(1A) mutants: T88H, Y92H, G196I, G196W and D278E. 2. G196W and G196I retained significant binding and IP-production properties, indicating that bulky substituents in position 196 did not affect the interaction of AngII's C-terminal carboxyl with Lys(199) located three residues below. 3. Although the T88A mutation did not affect binding, the T88H mutant had greatly decreased affinity for AngII, suggesting that substitution of Thr(88) by His might hinder binding through an indirect effect. 4. The Y92H mutation caused loss of affinity for AngII that was much less pronounced than that reported for Y92A, indicating that His in that position can fulfil part of the requirements for binding. 5. Replacing Asp(278) by Glu caused a much smaller reduction in affinity than replacing it by Ala, indicating the importance of Asp's beta-carboxyl group for AngII binding. 6. Mutations in residues Thr(88), Tyr(92) and Asp(278) greatly reduced affinity for AngII but not for Sar(1) Leu(8)-AngII, suggesting unfavourable interactions between these residues and AngII's aspartic acid side-chain or N-terminal amino group, which might account for the proposed role of the N-terminal amino group of AngII in the agonist-induced desensitization (tachyphylaxis) of smooth muscles.

Theoretical study on mutation-induced activation of the luteinizing hormone receptor

Here, three-dimensional model building and molecular dynamics simulations of the luteinizing hormone receptor have been employed to generate hypotheses about the molecular mechanisms underlying the activation of the receptor induced by naturally occurring activating mutations. The comparative analysis of the wild-type receptor and of 16 constitutively active or inactive mutants has been instrumental in inferring the structural/dynamic features which could characterize the inactive and the active forms of the receptor. These features have been also employed for predicting the functional behavior of new receptor mutants. The results of this study might provide a structural framework to interpret the pathological effects induced by mutations of the luteinizing hormone receptor. In addition, the proposed theoretical model could be useful for engineering new mutations or ligands able to modulate receptor function.

The alpha 1a and alpha 1b-adrenergic receptor subtypes: molecular mechanisms of receptor activation and of drug action

In this chapter we summarize some aspects of the structure-functional relationship of the alpha 1a and alpha 1b-adrenergic receptor subtypes related to the receptor activation process as well as the effect of different alpha-blockers on the constitutive activity of the receptor. Molecular modeling of the alpha 1a and alpha 1b-adrenergic receptor subtypes and computational simulation of receptor dynamics were useful to interpret the experimental findings derived from site directed mutagenesis studies.

Oxytocin receptors and cholesterol: interaction and regulation

Cholesterol affects the ligand binding function of the oxytocin receptor in a highly specific manner. While the structurally-related cholecystokinin receptor shows a strong correlation between the membrane fluidity and its binding function, the oxytocin receptor behaves differently. A stringent and unique requirement of the affinity state of the oxytocin receptor for structural features of the sterol molecule has been found. The molecular requirements differ both from those postulated for sterol-phospholipid interactions and from those known to be necessary for the activity of other proteins. Employing a new detergent-free subcellular fractionation protocol, a two-fold enrichment of the oxytocin receptors (10-15% of total receptors) has been detected in the cholesterol-rich, caveolin-containing membrane domains of the plasma membrane. While most of the properties of the oxytocin receptors were indistinguishable in cholesterol-poor versus cholesterol-rich membrane compartments, high-affinity oxytocin receptors localised in caveolin-enriched low-density membranes showed about a 3-fold higher stability against thermal denaturation at 37 degrees C compared with the oxytocin receptors localised in high-density membranes. Moreover, addition of cholesterol to the cholesterol-poor high-density membranes fully protected the oxytocin receptors against thermal denaturation and partially rescued high-affinity oxytocin binding. Although the membrane fluidity of the caveolin-enriched domains was lower than that in the high-density membranes, there was no correlation between the stability of oxytocin receptors and the fluidity level of the membrane domains. Finally, in a molecular modelling approach a putative cholesterol binding motif on the extracellular surface of the oxytocin receptor was found.

Activation mechanism of human oxytocin receptor: a combined study of experimental and computer-simulated mutagenesis

The aim of this study was to investigate the molecular changes associated with the transition of the human oxytocin receptor from its inactive to its active states. Mutation of the conserved arginine of the glutamate/aspartate-arginine-tyrosine motif located in the second intracellular domain gave rise to the first known constitutively active oxytocin receptor (R137A), whereas mutation of the aspartic acid located in the second transmembrane domain led to an inactive receptor (D85A). The structural features of the constitutively active and inactive receptor mutants were compared with those of the wild type in its free and agonist-bound states. The results suggest that, although differently triggered, the activation process induced by the agonist and the activating mutation are characterized by the opening of a solvent exposed site formed by the 2nd intracellular loop, the cytosolic extension of helix 5, and the 3rd intracellular loop; on the contrary, the D85A mutation prevents oxytocin from triggering the opening of a cytosolic site. On the basis of these findings, we hypothesize that this cytosolic crevice plays an important role in G protein recognition. Finally, comparative analysis of the free- and agonist-bound forms of the wild-type oxytocin receptor and alpha1B adrenergic receptor suggests that the highly conserved polar amino acids and the seven helices play similar mechanistic roles in the different G protein-coupled receptors.