Leonardo Pardo

Crystal structure of the µ-opioid receptor bound to a morphinan antagonist

Opium is one of the world’s oldest drugs, and its derivatives morphine and codeine are among the most used clinical drugs to relieve severe pain. These prototypical opioids produce analgesia as well as many undesirable side effects (sedation, apnoea and dependence) by binding to and activating the G-protein-coupled µ-opioid receptor (µ-OR) in the central nervous system. Here we describe the 2.8 Å crystal structure of the mouse µ-OR in complex with an irreversible morphinan antagonist. Compared to the buried binding pocket observed in most G-protein-coupled receptors published so far, the morphinan ligand binds deeply within a large solvent-exposed pocket. Of particular interest, the µ-OR crystallizes as a two-fold symmetrical dimer through a four-helix bundle motif formed by transmembrane segments 5 and 6. These high-resolution insights into opioid receptor structure will enable the application of structure-based approaches to develop better drugs for the management of pain and addiction.

Ligand-specific regulation of the extracellular surface of a G-protein-coupled receptor

G-protein-coupled receptors (GPCRs) are seven-transmembrane proteins that mediate most cellular responses to hormones and neurotransmitters. They are the largest group of therapeutic targets for a broad spectrum of diseases. Recent crystal structures of GPCRs have revealed structural conservation extending from the orthosteric ligand-binding site in the transmembrane core to the cytoplasmic G-protein-coupling domains. In contrast, the extracellular surface (ECS) of GPCRs is remarkably diverse and is therefore an ideal target for the discovery of subtype-selective drugs. However, little is known about the functional role of the ECS in receptor activation, or about conformational coupling of this surface to the native ligand-binding pocket. Here we use NMR spectroscopy to investigate ligand-specific conformational changes around a central structural feature in the ECS of the β2 adrenergic receptor: a salt bridge linking extracellular loops 2 and 3. Small-molecule drugs that bind within the transmembrane core and exhibit different efficacies towards G-protein activation (agonist, neutral antagonist and inverse agonist) also stabilize distinct conformations of the ECS. We thereby demonstrate conformational coupling between the ECS and the orthosteric binding site, showing that drugs targeting this diverse surface could function as allosteric modulators with high subtype selectivity. Moreover, these studies provide a new insight into the dynamic behaviour of GPCRs not addressable by static, inactive-state crystal structures.

Glycoprotein hormone receptors: determinants in leucine-rich repeats responsible for ligand specificity

Glycoprotein hormone receptors [thyrotropin (TSHr), luteinizing hormone/chorionic gonadotropin (LH/CGr), follicle stimulating hormone (FSHr)] are rhodopsin‐like G protein‐coupled receptors with a large extracellular N‐terminal portion responsible for hormone recognition and binding. In structural models, this ectodomain is composed of two cysteine clusters flanking nine leucine‐rich repeats (LRRs). The LRRs form a succession of β‐strands and α‐helices organized into a horseshoe‐shaped structure. It has been proposed that glycoprotein hormones interact with residues of the β‐strands making the concave surface of the horseshoe. Gain‐of‐function homology scanning of the β‐strands of glycoprotein hormone receptors allowed identification of the critical residues responsible for the specificity towards human chorionic gonadotropin (hCG). Substitution of eight or two residues of the LH/CGr into the TSHr or FSHr, respectively, resulted in constructs displaying almost the same affinity and sensitivity for hCG as wild‐type LH/CGr. Molecular dynamics simulations and additional site‐directed mutagenesis provided a structural rationale for the evolution of binding specificity in this duplicated gene family.

Ligand-regulated oligomerization of β2-adrenoceptors in a model lipid bilayer

The β2‐adrenoceptor (β2AR) was one of the first Family A G protein‐coupled receptors (GPCRs) shown to form oligomers in cellular membranes, yet we still know little about the number and arrangement of protomers in oligomers, the influence of ligands on the organization or stability of oligomers, or the requirement for other proteins to promote oligomerization. We used fluorescence resonance energy transfer (FRET) to characterize the oligomerization of purified β2AR site‐specifically labelled at three different positions with fluorophores and reconstituted into a model lipid bilayer. Our results suggest that the β2AR is predominantly tetrameric following reconstitution into phospholipid vesicles. Agonists and antagonists have little effect on the relative orientation of protomers in oligomeric complexes. In contrast, binding of inverse agonists leads to significant increases in FRET efficiencies for most labelling pairs, suggesting that this class of ligand promotes tighter packing of protomers and/or the formation of more complex oligomers by reducing conformational fluctuations in individual protomers. The results provide new structural insights into β2AR oligomerization and suggest a possible mechanism for the functional effects of inverse agonists.

Serine and Threonine Residues Bend α-Helices in the χ1 = g− Conformation

The relationship between the Ser, Thr, and Cys side-chain conformation (chi(1) = g(-), t, g(+)) and the main-chain conformation (phi and psi angles) has been studied in a selection of protein structures that contain alpha-helices. The statistical results show that the g(-) conformation of both Ser and Thr residues decreases their phi angles and increases their psi angles relative to Ala, used as a control. The additional hydrogen bond formed between the O(gamma) atom of Ser and Thr and the i-3 or i-4 peptide carbonyl oxygen induces or stabilizes a bending angle in the helix 3-4 degrees larger than for Ala. This is of particular significance for membrane proteins. Incorporation of this small bending angle in the transmembrane alpha-helix at one side of the cell membrane results in a significant displacement of the residues located at the other side of the membrane. We hypothesize that local alterations of the rotamer configurations of these Ser and Thr residues may result in significant conformational changes across transmembrane helices, and thus participate in the molecular mechanisms underlying transmembrane signaling. This finding has provided the structural basis to understand the experimentally observed influence of Ser residues on the conformational equilibrium between inactive and active states of the receptor, in the neurotransmitter subfamily of G protein-coupled receptors.

An activation switch in the rhodopsin family of G protein coupled receptors: The thyrotropin receptor

We aimed at understanding molecular events involved in the activation of a member of the G protein-coupled receptor family, the thyrotropin receptor. We have focused on the transmembrane region and in particular on a network of polar interactions between highly conserved residues. Using molecular dynamics simulations and site-directed mutagenesis techniques we have identified residue Asn-7.49, of the NPxxY motif of TM 7, as a molecular switch in the mechanism of thyrotropin receptor (TSHr) activation. Asn-7.49 appears to adopt two different conformations in the inactive and active states. These two states are characterized by specific interactions between this Asn and polar residues in the transmembrane domain. The inactive gauche+ conformation is maintained by interactions with residues Thr-6.43 and Asp-6.44. Mutation of these residues into Ala increases the constitutive activity of the receptor by factors of ∼14 and ∼10 relative to wild type TSHr, respectively. Upon receptor activation Asn-7.49 adopts the trans conformation to interact with Asp-2.50 and a putatively charged residue that remains to be identified. In addition, the conserved Leu-2.46 of the (N/S)LxxxD motif also plays a significant role in restraining the receptor in the inactive state because the L2.46A mutation increases constitutive activity by a factor of ∼13 relative to wild type TSHr. As residues Leu-2.46, Asp-2.50, and Asn-7.49 are strongly conserved, this molecular mechanism of TSHr activation can be extended to other members of the rhodopsin-like family of G protein-coupled receptors.

The Role of Internal Water Molecules in the Structure and Function of the Rhodopsin Family of G Protein‐Coupled Receptors

Membrane receptors coupled to guanine nucleotide-binding regulatory proteins (commonly known as G protein-coupled ACHTUNGTRENNUNGreceptors, GPCRs) constitute one of the most attractive pharmaceutical targets, as around 40% of clinically prescribed drugs and 25% of the top-selling drugs act at these receptors. GPCRs are receptors for sensory signals of external origin such as odors, pheromones, or tastes; and for endogenous signals such as neurotransmitters, (neuro)peptides, divalent cations, proteases, glycoprotein hormones, and purine ligands. Phylogenetic analyses of the human genome have permitted GPCR sequences to be classified into five main families: rhodopsin (class A or family 1), secretin (class B or family 2), glutamate (class C or family 3), adhesion, and frizzled/ taste2. Specialized databases of GPCRs can be found at http://www.gpcr.org/7tm, http://gris.ulb.ac.be/, and http:// www.iuphar-db.org. Due to the low natural abundance of GPCRs and the difficulty in producing and purifying recombinant protein, only one member of this family, rhodopsin, the photoreceptor protein of rod cells, has been crystallized so far. Five structural models of inactive rhodopsin are available at the Protein Data Bank, at resolutions of 2.8 A (PDB IDs: 1F88 and 1HZX), 2.65 A (1GZM), 2.6 A (1L9H), and 2.2 A (1U19). Structural models of rhodopsin photointermediates such as bathorhodopsin (2G87), lumirhodopsin (2HPY), metarhodopsin I, and a photoactivated deprotonated intermediate reminiscent of metarhodopsin II (2I37) are also available. Rhodopsin is formed by an extracellular N terminus of four b-strands, seven transmembrane helices (TM1 to TM7) connected by alternating intracellular (I1 to I3) and extracellular (E1 to E3) hydrophilic loops, a disulfide bridge between E2 and TM3, and a cytoplasmic C terminus containing an a-helix (Hx8) parallel to the cell membrane. Statistical analysis of the residues forming the TM helices of the rhodopsin family of GPCRs shows a large number of conserved sequence patterns; this suggests a common TM structure. Thus, the availability of the rhodopsin structure allows the use of homology modeling techniques to build three-dimensional models of other homologous GPCRs. The putative structural homology between rhodopsin and other GPCRs probably does not extend to the extracellular domain, since the extracellular N terminus and loop fragments are highly variable in length and amino acid content. The class A family of GPCRs contains highly conserved Pro residues in the middle of TMs 5 (P5.50, conserved in 77% of the sequences), 6 (P6.50, 100%), and 7 (P7.50, 96%; residues are identified by the generic numbering scheme of Ballesteros and Weinstein, which allows easy comparison among residues in the 7TM segments of different receptors). Pro residues are normally observed in the TM helices of membrane proteins where they usually induce a significant distortion named a “Pro-kink”. This break arises in order to avoid a steric clash between the pyrrolidine ring of the Pro side chain (at position i) and the carbonyl oxygen of the residue in the preceding turn (position i 4) and leads to a distortion of the helical structure. However, TM segments of rhodopsin, either with or without Pro residues in their sequence are far from being standard Pro-kinked or ideal helices, respectively. Their distortions are energetically stabilized through complementary intraand interhelical interactions involving polar side chains, backbone carbonyls, and, in some cases, specific structural and functional water molecules embedded in the TM bundle. Here we review the role of these water molecules in the structure and function of GPCRs and in building computergenerated homology models of class A GPCRs. We propose that water molecules present in the vicinity of highly conserved motifs are most likely present in the rhodopsin family of GPCRs, being another conserved structural element in the family.

Membrane Protein Simulations Using AMBER Force Field and Berger Lipid Parameters.

AMBER force fields are among the most commonly used in molecular dynamics (MD) simulations of proteins. Unfortunately, they lack a specific set of lipid parameters, thus limiting its use in membrane protein simulations. In order to overcome this limitation we assessed whether the widely used united-atom lipid parameters described by Berger and co-workers could be used in conjunction with AMBER force fields in simulations of membrane proteins. Thus, free energies of solvation in water and in cyclohexane, and free energies of water to cyclohexane transfer, were computed by thermodynamic integration procedures for neutral amino acid side-chains employing AMBER99, AMBER03, and OPLS-AA amino acid force fields. In addition, MD simulations of three membrane proteins in a POPC lipid bilayer, the β2 adrenergic G protein-coupled receptor, Aquaporin-1, and the outer membrane protein Omp32, were performed with the aim of comparing the AMBER99SB/Berger combination of force fields with the OPLS-AA/Berger combination. We have shown that AMBER99SB and Berger force fields are compatible, they provide reliable free energy estimations relative to experimental values, and their combination properly describes both membrane and protein structural properties. We then suggest that the AMBER99SB/Berger combination is a reliable choice for the simulation of membrane proteins, which links the easiness of ligand parametrization and the ability to reproduce secondary structure of AMBER99SB force field with the largely validated Berger lipid parameters.

Activation of CCR5 by chemokines involves an aromatic cluster between transmembrane helices 2 and 3

CCR5 is a G protein-coupled receptor responding to four natural agonists, the chemokines RANTES (regulated on activation normal T cell expressed and secreted), macrophage inflammatory protein (MIP)-1α, MIP-1β, and monocyte chemotactic protein (MCP)-2, and is the main co-receptor for the macrophage-tropic human immunodeficiency virus strains. We have previously identified a structural motif in the second transmembrane helix of CCR5, which plays a crucial role in the mechanism of receptor activation. We now report the specific role of aromatic residues in helices 2 and 3 of CCR5 in this mechanism. Using site-directed mutagenesis and molecular modeling in a combined approach, we demonstrate that a cluster of aromatic residues at the extracellular border of these two helices are involved in chemokine-induced activation. These aromatic residues are involved in interhelical interactions that are key for the conformation of the helices and govern the functional response to chemokines in a ligand-specific manner. We therefore suggest that transmembrane helices 2 and 3 contain important structural elements for the activation mechanism of chemokine receptors, and possibly other related receptors as well.

Ser and Thr Residues Modulate the Conformation of Pro-Kinked Transmembrane α-Helices

Functionally required conformational plasticity of transmembrane proteins implies that specific structural motifs have been integrated in transmembrane helices. Surveying a database of transmembrane helices and the large family of G-protein coupled receptors we identified a series of overrepresented motifs associating Pro with either Ser or Thr. Thus, we have studied the conformation of Pro-kinked transmembrane helices containing Ser or Thr residues, in both g+ and g- rotamers, by molecular dynamics simulations in a hydrophobic environment. Analysis of the simulations shows that Ser or Thr can significantly modulate the deformation of the Pro. A series of motifs, such as (S/T)P and (S/T)AP in the g+ rotamer and the TAP and PAA(S/T) motifs in the g- rotamer, induce an increase in bending angle of the helix compared to a standard Pro-kink, apparently due to the additional hydrogen bond formed between the side chain of Ser/Thr and the backbone carbonyl oxygen. In contrast, (S/T)AAP and PA(S/T) motifs, in both g+ and g-, and PAA(S/T) in g+ rotamers decrease the bending angle of the helix by either reducing the steric clash between the pyrrolidine ring of Pro and the helical backbone, or by adding a constrain in the form of a hydrogen bond in the curved-in face of the helix. Together with a number of available experimental data, our results strongly suggest that association of Ser and Thr with Pro is commonly used in transmembrane helices to accommodate the structural needs of specific functions.