Structural Modeling of Cytokine-Receptor-JAK2 Signaling Complexes Using AlphaFold Multimer

Homodimeric class 1 cytokine receptors include the erythropoietin (EPOR), thrombopoietin (TPOR), granulocyte colony-stimulating factor 3 (CSF3R), growth hormone (GHR), and prolactin receptors (PRLR). These cell-surface single-pass transmembrane (TM) glycoproteins regulate cell growth, proliferation, and differentiation and induce oncogenesis. An active TM signaling complex consists of a receptor homodimer, one or two ligands bound to the receptor extracellular domains, and two molecules of Janus Kinase 2 (JAK2) constitutively associated with the receptor intracellular domains. Although crystal structures of soluble extracellular domains with ligands have been obtained for all of the receptors except TPOR, little is known about the structure and dynamics of the complete TM complexes that activate the downstream JAK-STAT signaling pathway. Three-dimensional models of five human receptor complexes with cytokines and JAK2 were generated here by using AlphaFold Multimer. Given the large size of the complexes (from 3220 to 4074 residues), the modeling required a stepwise assembly from smaller parts, with selection and validation of the models through comparisons with published experimental data. The modeling of active and inactive complexes supports a general activation mechanism that involves ligand binding to a monomeric receptor followed by receptor dimerization and rotational movement of the receptor TM α-helices, causing proximity, dimerization, and activation of associated JAK2 subunits. The binding mode of two eltrombopag molecules to the TM α-helices of the active TPOR dimer was proposed. The models also help elucidate the molecular basis of oncogenic mutations that may involve a noncanonical activation route. Models equilibrated in explicit lipids of the plasma membrane are publicly available.

Structural modeling of cytokine-receptor-JAK2 signaling complexes using AlphaFold Multimer

Homodimeric class 1 cytokine receptors include the erythropoietin (EPOR), thrombopoietin (TPOR), granulocyte colony-stimulating factor 3 (CSF3R), growth hormone (GHR), and prolactin receptors (PRLR). They are cell-surface single-pass transmembrane (TM) glycoproteins that regulate cell growth, proliferation, and differentiation and induce oncogenesis. An active TM signaling complex consists of a receptor homodimer, one or two ligands bound to the receptor extracellular domains and two molecules of Janus Kinase 2 (JAK2) constitutively associated with the receptor intracellular domains. Although crystal structures of soluble extracellular domains with ligands have been obtained for all the receptors except TPOR, little is known about the structure and dynamics of the complete TM complexes that activate the downstream JAK-STAT signaling pathway. Three-dimensional models of five human receptor complexes with cytokines and JAK2 were generated using AlphaFold Multimer. Given the large size of the complexes (from 3220 to 4074 residues), the modeling required a stepwise assembly from smaller parts with selection and validation of the models through comparisons with published experimental data. The modeling of active and inactive complexes supports a general activation mechanism that involves ligand binding to a monomeric receptor followed by receptor dimerization and rotational movement of the receptor TM α-helices causing proximity, dimerization, and activation of associated JAK2 subunits. The binding mode of two eltrombopag molecules to TM α-helices of the active TPOR dimer was proposed. The models also help elucidating the molecular basis of oncogenic mutations that may involve non-canonical activation route. Models equilibrated in explicit lipids of the plasma membrane are publicly available.

Effects of calreticulin mutations on cell transformation and immunity

Myeloproliferative neoplasms (MPNs) are cancers involving dysregulated production and function of myeloid lineage hematopoietic cells. Among MPNs, Essential thrombocythemia (ET), Polycythemia Vera (PV) and Myelofibrosis (MF), are driven by mutations that activate the JAK-STAT signalling pathway. Somatic mutations of calreticulin (CRT), an endoplasmic reticulum (ER)-localized lectin chaperone, are driver mutations in approximately 25% of ET and 35% of MF patients. The MPN-linked mutant CRT proteins have novel frameshifted carboxy-domain sequences and lack an ER retention motif, resulting in their secretion. Wild type CRT is a regulator of ER calcium homeostasis and plays a key role in the assembly of major histocompatibility complex (MHC) class I molecules, which are the ligands for antigen receptors of CD8+ T cells. Mutant CRT-linked oncogenesis results from the dysregulation of calcium signalling in cells and the formation of stable complexes of mutant CRT with myeloproliferative leukemia (MPL) protein, followed by downstream activation of the JAK-STAT signalling pathway. The intricate participation of CRT in ER protein folding, calcium homeostasis and immunity suggests the involvement of multiple mechanisms of mutant CRT-linked oncogenesis. In this review, we highlight recent findings related to the role of MPN-linked CRT mutations in the dysregulation of calcium homeostasis, MPL activation and immunity.

Membranome 3.0: Database of single-pass membrane proteins with AlphaFold models

The Membranome database provides comprehensive structural information on single‐pass (i.e., bitopic) membrane proteins from six evolutionarily distant organisms, including protein–protein interactions, complexes, mutations, experimental structures, and models of transmembrane α‐helical dimers. We present a new version of this database, Membranome 3.0, which was significantly updated by revising the set of 5,758 bitopic proteins and incorporating models generated by AlphaFold 2 in the database. The AlphaFold models were parsed into structural domains located at the different membrane sides, modified to exclude low‐confidence unstructured terminal regions and signal sequences, validated through comparison with available experimental structures, and positioned with respect to membrane boundaries. Membranome 3.0 was re‐developed to facilitate visualization and comparative analysis of multiple 3D structures of proteins that belong to a specified family, complex, biological pathway, or membrane type. New tools for advanced search and analysis of proteins, their interactions, complexes, and mutations were included. The database is freely accessible at https://membranome.org.

Comparative Molecular Dynamics Simulation Studies of Realistic Eukaryotic, Prokaryotic, and Archaeal Membranes

We present a comparative all-atom molecular dynamics simulation study of 18 biomembrane systems with lipid compositions corresponding to eukaryotic, bacterial, and archaebacterial membranes together with three single-component lipid bilayers. A total of 105 lipid types used in this study include diverse sterols and glycerol-based lipids with acyl chains of various lengths, unsaturation degrees, and branched or cyclic moieties. Our comparative analysis provides deeper insight into the influences of sterols and lipid unsaturation on the structural and mechanical properties of these biomembranes, including water permeation into the membrane hydrocarbon core. For sterol-containing membranes, sterol fraction is correlated with the membrane thickness, the area compressibility modulus, and lipid order but anticorrelated with the area per lipid and sterol tilt angles. Similarly, for all 18 biomembranes, lipid order is correlated with the membrane thickness and area compressibility modulus. Sterols and lipid unsaturation produce opposite effects on membrane thickness, but only sterols influence water permeation into the membrane. All membrane systems are accessible for public use in CHARMM-GUI Archive. They can be used as templates to expedite future modeling of realistic cell membranes with transmembrane and peripheral membrane proteins to study their structure, dynamics, molecular interactions, and function in a nativelike membrane environment.

Spatial arrangement of proteins in planar and curved membranes by PPM 3.0

Cellular protrusions, invaginations, and many intracellular organelles have strongly curved membrane regions. Transmembrane and peripheral membrane proteins that induce, sense, or stabilize such regions cannot be properly fitted into a single flat bilayer. To treat such proteins, we developed a new method and a web tool, PPM 3.0, for positioning proteins in curved or planar, single or multiple membranes. This method determines the energetically optimal spatial position, the hydrophobic thickness, and the radius of intrinsic curvature of a membrane-deforming protein structure by arranging it in a single or several sphere-shaped or planar membrane sections. In addition, it can define the lipid-embedded regions of a protein that simultaneously spans several membranes or determine the optimal position of a peptide in a spherical micelle. The PPM 3.0 web server operates with 17 types of biological membranes and 4 types of artificial bilayers. It is publicly available at https://opm.phar.umich.edu/ppm_server3. PPM 3.0 was applied to identify and characterize arrangements in membranes of 128 proteins with a significant intrinsic curvature, such as BAR domains, annexins, Piezo, and MscS mechanosensitive channels, cation-chloride cotransporters, as well as mitochondrial ATP synthases, calcium uniporters, and TOM complexes. These proteins form large complexes that are mainly localized in mitochondria, plasma membranes, and endosomes. Structures of bacterial drug efflux pumps, AcrAB-TolC, MexAB-OrpM, and MacAB-TolC, were positioned in both membranes of the bacterial cell envelop, while structures of multimeric gap-junction channels were arranged in two opposed cellular membranes.

Thermodynamics-Based Molecular Modeling of α-Helices in Membranes and Micelles

The Folding of Membrane-Associated Peptides (FMAP) method was developed for modeling α-helix formation by linear peptides in micelles and lipid bilayers. FMAP 2.0 identifies locations of α-helices in the amino acid sequence, generates their three-dimensional models in planar bilayers or spherical micelles, and estimates their thermodynamic stabilities and tilt angles, depending on temperature and pH. The method was tested for 723 peptides (926 data points) experimentally studied in different environments and for 170 single-pass transmembrane (TM) proteins with available crystal structures. FMAP 2.0 detected more than 95% of experimentally observed α-helices with an average error in helix end determination of around 2, 3, 4, and 5 residues per helix for peptides in water, micelles, bilayers, and TM proteins, respectively. Helical and nonhelical residue states were predicted with an accuracy from 0.86 to 0.96, and the Matthews correlation coefficient was from 0.64 to 0.88 depending on the environment. Experimental micelle- and membrane-binding energies and tilt angles of peptides were reproduced with a root-mean-square deviation of around 2 kcal/mol and 7°, respectively. The TM and non-TM states of hydrophobic and pH-triggered α-helical peptides in various lipid bilayers were reproduced in more than 95% of cases. The FMAP 2.0 web server (https://membranome.org/fmap) is publicly available to explore the structural polymorphism of antimicrobial, cell-penetrating, fusion, and other membrane-binding peptides, which is important for understanding the mechanisms of their biological activities.

Mechanism of mutant calreticulin-mediated activation of the thrombopoietin receptor in cancers

Myeloproliferative neoplasms (MPNs) are frequently driven by mutations within the C-terminal domain (C-domain) of calreticulin (CRT). CRT_Del52 and CRT_Ins5 are recurrent mutations. Oncogenic transformation requires both mutated CRT and the thrombopoietin receptor (Mpl), but the molecular mechanism of CRT-mediated constitutive activation of Mpl is unknown. We show that the acquired C-domain of CRT_Del52 mediates both Mpl binding and disulfide-linked CRT_Del52 dimerization. Cysteine mutations within the novel C-domain (C400A and C404A) and the conserved N-terminal domain (N-domain; C163A) of CRT_Del52 are required to reduce disulfide-mediated dimers and multimers of CRT_Del52. Based on these data and published structures of CRT oligomers, we identify an N-domain dimerization interface relevant to both WT CRT and CRT_Del52. Elimination of disulfide bonds and ionic interactions at both N-domain and C-domain dimerization interfaces is required to abrogate the ability of CRT_Del52 to mediate cell proliferation via Mpl. Thus, MPNs exploit a natural dimerization interface of CRT combined with C-domain gain of function to achieve cell transformation.

Aromatic-Amine Pendants Produce Highly Potent and Efficacious Mixed Efficacy μ-Opioid Receptor (MOR)/δ-Opioid Receptor (DOR) Peptidomimetics with Enhanced Metabolic Stability

We previously reported a novel SAR campaign that converted a metabolically unstable series of μ-opioid receptor (MOR) agonist/δ-opioid receptor (DOR) antagonist bicyclic core peptidomimetics with promising analgesic activity and reduced abuse liabilities into a more stable series of benzylic core analogues. Herein, we expanded the SAR of that campaign and determined that the incorporation of amines into the benzylic pendant produces enhanced MOR-efficacy in this series, whereas the reincorporation of an aromatic ring into the pendant enhanced MOR-potency. Two compounds, which contain a piperidine (14) or an isoindoline (17) pendant, retained the desired opioid profile in vitro, possessed metabolic half-lives of greater than 1 h in mouse liver microsomes (MLMs), and were active antinociceptive agents in the acetic acid stretch assay (AASA) at subcutaneous doses of 1 mg/kg.

Incorporation of Agouti-Related Protein (AgRP) Human Single Nucleotide Polymorphisms (SNPs) in the AgRP-Derived Macrocyclic Scaffold c[Pro-Arg-Phe-Phe-Asn-Ala-Phe-dPro] Decreases Melanocortin-4 Receptor Antagonist Potency and Results in the Discovery of Melanocortin-5 Receptor Antagonists

While the melanocortin receptors (MCRs) are known to be involved in numerous biological pathways, the potential roles of the MC5R have not been clearly elucidated in humans. Agouti-related protein (AgRP), an MC3R/MC4R antagonist and MC4R inverse agonist, contains an exposed β-hairpin loop composed of six residues (Arg-Phe-Phe-Asn-Ala-Phe) that is imperative for binding and function. Within this active loop of AgRP, four human missense polymorphisms were deposited into the NIH Variation Viewer database. These polymorphisms, Arg111Cys, Arg111His, Phe112Tyr, and Ala115Val (AgRP full-length numbering), were incorporated into the peptide macrocycles c[Pro¹-Arg²-Phe³-Phe⁴-Xaa⁵-Ala⁶-Phe⁷-dPro⁸], where Xaa was Dap⁵ or Asn⁵, to explore the functional effects of these naturally occurring substitutions in a simplified AgRP scaffold. All peptides lowered potency at least 10-fold in a cAMP accumulation assay compared to the parent sequences at the MC4Rs. Compounds MDE 6-82-3c, ZMK 2-82, MDE 6-82-1c, ZMK 2-85, and ZMK 2-112 are also the first AgRP-based chemotypes that antagonize the MC5R.

Phospholipid Component Defines Pharmacokinetic and Pharmacodynamic Properties of Synthetic High-Density Lipoproteins

Synthetic high-density lipoprotein (sHDL) nanoparticles composed of apolipoprotein A-I mimetic peptide and phospholipids have been shown to reduce atherosclerosis in animal models. Cholesterol is mobilized from atheroma macrophages by sHDL into the blood compartment and delivered to the liver for elimination. Historically, sHDL drug discovery efforts were focused on optimizing peptide sequences for interaction with cholesterol cellular transporters rather than understanding how both sHDL components, peptide and lipid, influence its pharmacokinetic and pharmacodynamic profiles. We designed two sets of sHDL having either identical phospholipid but variable peptide sequences with different plasma stability or identical peptide and phospholipids with variable fatty acid chain length and saturation. We found that sHDL prepared with proteolytically stable 22A-P peptide had 2-fold longer circulation half-time relative to the less stable 22A peptide. Yet, longer half-life did not translate into any improvement in cholesterol mobilization. In contrast, sHDL with variable phospholipid compositions showed significant differences in phospholipid PK, with distearoyl phosphatidylcholine-based sHDL demonstrating the longest half-life of 6.0 hours relative to 1.0 hour for palmitoyl-oleoyl phosphatidylcholine-based sHDL. This increase in half-life corresponded to an approx. 6.5-fold increase in the area under the curve for the mobilized cholesterol. Therefore, the phospholipid component in sHDL plays a major role in cholesterol mobilization in vivo and should not be overlooked in the design of future sHDL. SIGNIFICANCE STATEMENT: The phospholipid composition in sHDL plays a critical role in determining half-life and cholesterol mobilization in vivo.

Novel Dimethyltyrosine-Tetrahydroisoquinoline Peptidomimetics with Aromatic Tetrahydroisoquinoline Substitutions Show in Vitro Kappa and Mu Opioid Receptor Agonism

The dimethyltyrosine-tetrahydroisoquinoline (Dmt-Tiq) scaffold was originally developed in the production of selective delta opioid receptor (DOR) antagonists. Installation of a 7-benzyl pendant on the tetrahydroisoquinoline core of this classic opioid scaffold introduced kappa opioid receptor (KOR) agonism. Further modification of this pendant resulted in retention of KOR agonism and the addition of mu opioid receptor (MOR) partial agonism, a bifunctional profile with potential to be used in the treatment of cocaine addiction.

Membranome 2.0: database for proteome-wide profiling of bitopic proteins and their dimers

Motivation

Structural studies of TM domains of single-spanning (bitopic) membrane proteins are impeded by their instability, flexibility and heterogeneity. The new computational method TMDOCK allows reliable modeling of homodimers of transmembrane (TM) α-helices on a proteomic scale.

Results

3D models of 2129 parallel homodimers formed by TM α-helices of bitopic proteins from six evolutionarily distant organisms were modeled by TMDOCK, verified through experimental data available for nearly 600 proteins, and included in the Membranome database (v.2.0) along with related information to facilitate structural and evolutionary analysis of bitopic proteins.

Availability and implementation

http://membranome.org.

Contact

almz@umich.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

PerMM: A Web Tool and Database for Analysis of Passive Membrane Permeability and Translocation Pathways of Bioactive Molecules

The PerMM web server and database were developed for quantitative analysis and visualization of passive translocation of bioactive molecules across lipid membranes. The server is the first physics-based web tool that calculates membrane binding energies and permeability coefficients of diverse molecules through artificial and natural membranes (phospholipid bilayers, PAMPA-DS, blood-brain barrier, and Caco-2/MDCK cell membranes). It also visualizes the transmembrane translocation pathway as a sequence of translational and rotational positions of a permeant as it moves across the lipid bilayer, along with the corresponding changes in solvation energy. The server can be applied for prediction of permeability coefficients of compounds with diverse chemical scaffolds to facilitate selection and optimization of potential drug leads. The complementary PerMM database allows comparison of computationally and experimentally determined permeability coefficients for more than 500 compounds in different membrane systems. The website and database are freely accessible at https://permm.phar.umich.edu/ .

Physics-Based Method for Modeling Passive Membrane Permeability and Translocation Pathways of Bioactive Molecules

Assessment of permeability is a critical step in the drug development process for selection of drug candidates with favorable ADME properties. We have developed a novel physics-based method for fast computational modeling of passive permeation of diverse classes of molecules across lipid membranes. The method is based on heterogeneous solubility-diffusion theory and operates with all-atom 3D structures of solutes and the anisotropic solvent model of the lipid bilayer characterized by transbilayer profiles of dielectric and hydrogen bonding capacity parameters. The optimal translocation pathway of a solute is determined by moving an ensemble of representative conformations of the molecule through the dioleoyl-phosphatidylcholine (DOPC) bilayer and optimizing their rotational orientations in every point of the transmembrane trajectory. The method calculates (1) the membrane-bound state of the solute molecule; (2) free energy profile of the solute along the permeation pathway; and (3) the permeability coefficient obtained by integration over the transbilayer energy profile and assuming a constant size-dependent diffusivity along the membrane normal. The accuracy of the predictions was evaluated against experimental permeability coefficients measured in pure lipid membranes (for 78 compounds, R² was 0.88 and rmse was 1.15 log units), PAMPA-DS (for 280 compounds, R² was 0.75 and rmse was 1.59 log units), BBB (for 182 compounds, R² was 0.69 and rmse was 0.87 log units), and Caco-2/MDCK assays (for 165 compounds, R² was 0.52 and rmse was 0.89 log units).

Placement of Hydroxy Moiety on Pendant of Peptidomimetic Scaffold Modulates Mu and Kappa Opioid Receptor Efficacy

In an effort to expand the structure–activity relationship (SAR) studies of a series of mixed-efficacy opioid ligands, peptidomimetics that incorporate methoxy and hydroxy groups around a benzyl or 2-methylindanyl pendant on a tetrahydroquinoline (THQ) core of the peptidomimetics were evaluated. Compounds containing a methoxy or hydroxy moiety in the o- or m-positions increased binding affinity to the kappa opioid receptor (KOR), whereas compounds containing methoxy or hydroxy groups in the p-position decreased KOR affinity and reduced or eliminated efficacy at the mu opioid receptor (MOR). The results from a substituted 2-methylindanyl series aligned with the findings from the substituted benzyl series. Our studies culminated in the development of 8c, a mixed-efficacy MOR agonist/KOR agonist with subnanomolar binding affinity for both MOR and KOR.

Membranome: a database for proteome-wide analysis of single-pass membrane proteins

The Membranome database was developed to assist analysis and computational modeling of single-pass (bitopic) transmembrane (TM) proteins and their complexes by providing structural information about these proteins on a genomic scale. The database currently collects data on >6000 bitopic proteins from Homo sapiens, Arabidopsis thaliana, Dictyostelium discoideum, Saccharomyces cerevisiae, Escherichia coli and Methanocaldococcus jannaschii It presents the following data: (i) hierarchical classification of bitopic proteins into 15 functional classes, 689 structural superfamilies and 1404 families; (ii) 446 complexes of bitopic proteins with known three-dimensional (3D) structures classified into 129 families; (iii) computationally generated three-dimensional models of TM α-helices positioned in membranes; (iv) amino acid sequences, domain architecture, functional annotation and available experimental structures of bitopic proteins; (v) TM topology and intracellular localization, (vi) physical interactions between proteins from the database along with links to other resources. The database is freely accessible at http://membranome.org There is a variety of options for browsing, sorting, searching and retrieval of the content, including downloadable coordinate files of TM domains with calculated membrane boundaries.

Effects of N-Substitutions on the Tetrahydroquinoline (THQ) Core of Mixed-Efficacy μ-Opioid Receptor (MOR)/δ-Opioid Receptor (DOR) Ligands

N-Acetylation of the tetrahydroquinoline (THQ) core of a series of μ-opioid receptor (MOR) agonist/δ-opioid receptor (DOR) antagonist ligands increases DOR affinity, resulting in ligands with balanced MOR and DOR affinities. We report a series of N-substituted THQ analogues that incorporate various carbonyl-containing moieties to maintain DOR affinity and define the steric and electronic requirements of the binding pocket across the opioid receptors. 4h produced in vivo antinociception (ip) for 1 h at 10 mg/kg.

Gramicidin A disassembles large conductive clusters of its lysine-substituted derivatives in lipid membranes

N-terminally substituted lysine derivatives of gramicidin A (gA), [Lys1]gA and [Lys3]gA, but not glutamate- or aspartate-substituted peptides have been previously shown to cause the leakage of carboxyfluorescein from liposomes. Here, the leakage induction was also observed for [Arg1]gA and [Arg3]gA, while [His1]gA and [His3]gA were inactive at neutral pH. The Lys3-containing analogue with all tryptophans replaced by isoleucines did not induce liposome leakage, similar to gA. This suggests that the presence of both tryptophans and N-terminal cationic residues is critical for pore formation. Remarkably, the addition of gA blocked the leakage induced by [Lys3]gA. By examining with fluorescence correlation spectroscopy the peptide-induced leakage of fluorescent markers from liposomes, we estimated the diameter of pores responsible for the leakage to be about 1.6 nm. Transmission electron cryo-microscopy imaging of liposomes with [Lys3]gA showed that the liposomal membranes contained high electron density particles with a size of about 40 Å, suggesting the formation of peptide clusters. No such clusterization was observed in liposomes incorporating gA or a mixture of gA with [Lys3]gA. Three-dimensional reconstruction of the clusters was compatible with their pentameric arrangement. Based on experimental data and computational modeling, we suggest that the large pore formed by [Lys3]gA represents a barrel-stave oligomeric cluster formed by antiparallel double-stranded helical dimers (DH). In a tentative model, the pentamer of dimers may be stabilized by aromatic Trp-Trp and cation-π Trp-Lys interactions between the neighboring DHs. The inhibiting effect of gA on the [Lys3]gA-induced leakage can be attributed to breaking of cation-π interactions, which prevents peptide clusterization and pore formation.

Further Optimization and Evaluation of Bioavailable, Mixed-Efficacy μ-Opioid Receptor (MOR) Agonists/δ-Opioid Receptor (DOR) Antagonists: Balancing MOR and DOR Affinities

In a previously described peptidomimetic series, we reported the development of bifunctional μ-opioid receptor (MOR) agonist and δ-opioid receptor (DOR) antagonist ligands with a lead compound that produced antinociception for 1 h after intraperitoneal administration in mice. In this paper, we expand on our original series by presenting two modifications, both of which were designed with the following objectives: (1) probing bioavailability and improving metabolic stability, (2) balancing affinities between MOR and DOR while reducing affinity and efficacy at the κ-opioid receptor (KOR), and (3) improving in vivo efficacy. Here, we establish that, through N-acetylation of our original peptidomimetic series, we are able to improve DOR affinity and increase selectivity relative to KOR while maintaining the desired MOR agonist/DOR antagonist profile. From initial in vivo studies, one compound (14a) was found to produce dose-dependent antinociception after peripheral administration with an improved duration of action of longer than 3 h.

Use of Functional Polymorphisms To Elucidate the Peptide Binding Site of TAP Complexes

TAP1/TAP2 complexes translocate peptides from the cytosol to the endoplasmic reticulum lumen to enable immune surveillance by CD8(+) T cells. Peptide transport is preceded by peptide binding to a cytosol-accessible surface of TAP1/TAP2 complexes, but the location of the TAP peptide-binding pocket remains unknown. Guided by the known contributions of polymorphic TAP variants to peptide selection, we combined homology modeling of TAP with experimental measurements to identify several TAP residues that interact with peptides. Models for peptide-TAP complexes were generated, which indicate bent conformation for peptides. The peptide binding site of TAP is located at the hydrophobic boundary of the cytosolic membrane leaflet, with striking parallels to the glutathione binding site of NaAtm1, a transporter that functions in bacterial heavy metal detoxification. These studies illustrate the conservation of the ligand recognition modes of bacterial and mammalians transporters involved in peptide-guided cellular surveillance.

Life at the border: adaptation of proteins to anisotropic membrane environment

This review discusses main features of transmembrane (TM) proteins which distinguish them from water-soluble proteins and allow their adaptation to the anisotropic membrane environment. We overview the structural limitations on membrane protein architecture, spatial arrangement of proteins in membranes and their intrinsic hydrophobic thickness, co-translational and post-translational folding and insertion into lipid bilayers, topogenesis, high propensity to form oligomers, and large-scale conformational transitions during membrane insertion and transport function. Special attention is paid to the polarity of TM protein surfaces described by profiles of dipolarity/polarizability and hydrogen-bonding capacity parameters that match polarity of the lipid environment. Analysis of distributions of Trp resides on surfaces of TM proteins from different biological membranes indicates that interfacial membrane regions with preferential accumulation of Trp indole rings correspond to the outer part of the lipid acyl chain region-between double bonds and carbonyl groups of lipids. These "midpolar" regions are not always symmetric in proteins from natural membranes. We also examined the hydrophobic effect that drives insertion of proteins into lipid bilayer and different free energy contributions to TM protein stability, including attractive van der Waals forces and hydrogen bonds, side-chain conformational entropy, the hydrophobic mismatch, membrane deformations, and specific protein-lipid binding.

Translation of structure-activity relationships from cyclic mixed efficacy opioid peptides to linear analogues

Most opioid analgesics used in the treatment of pain are mu opioid receptor (MOR) agonists. While effective, there are significant drawbacks to opioid use, including the development of tolerance and dependence. However, the coadministration of a MOR agonist with a delta opioid receptor (DOR) antagonist slows the development of MOR-related side effects, while maintaining analgesia. We have previously reported a series of cyclic mixed efficacy MOR agonist/DOR antagonist ligands. Here we describe the transfer of key features from these cyclic analogs to linear sequences. Using the linear MOR/DOR agonist, Tyr-DThr-Gly-Phe-Leu-Ser-NH2 (DTLES), as a lead scaffold, we replaced Phe(4) with bulkier and/or constrained aromatic residues shown to confer DOR antagonism in our cyclic ligands. These replacements failed to confer DOR antagonism in the DTLES analogs, presumably because the more flexible linear ligands can adopt binding poses that will fit in the narrow binding pocket of the active conformations of both MOR and DOR. Nonetheless, the pharmacological profile observed in this series, high affinity and efficacy for MOR and DOR with selectivity relative to KOR, has also been shown to reduce the development of unwanted side effects. We further modified our lead MOR/DOR agonist with a C-terminal glucoserine to improve bioavailability. The resulting ligand displayed high efficacy and potency at both MOR and DOR and no efficacy at KOR.

Structural adaptations of proteins to different biological membranes

To gain insight into adaptations of proteins to their membranes, intrinsic hydrophobic thicknesses, distributions of different chemical groups and profiles of hydrogen-bonding capacities (α and β) and the dipolarity/polarizability parameter (π*) were calculated for lipid-facing surfaces of 460 integral α-helical, β-barrel and peripheral proteins from eight types of biomembranes. For comparison, polarity profiles were also calculated for ten artificial lipid bilayers that have been previously studied by neutron and X-ray scattering. Estimated hydrophobic thicknesses are 30-31Å for proteins from endoplasmic reticulum, thylakoid, and various bacterial plasma membranes, but differ for proteins from outer bacterial, inner mitochondrial and eukaryotic plasma membranes (23.9, 28.6 and 33.5Å, respectively). Protein and lipid polarity parameters abruptly change in the lipid carbonyl zone that matches the calculated hydrophobic boundaries. Maxima of positively charged protein groups correspond to the location of lipid phosphates at 20-22Å distances from the membrane center. Locations of Tyr atoms coincide with hydrophobic boundaries, while distributions maxima of Trp rings are shifted by 3-4Å toward the membrane center. Distributions of Trp atoms indicate the presence of two 5-8Å-wide midpolar regions with intermediate π* values within the hydrocarbon core, whose size and symmetry depend on the lipid composition of membrane leaflets. Midpolar regions are especially asymmetric in outer bacterial membranes and cell membranes of mesophilic but not hyperthermophilic archaebacteria, indicating the larger width of the central nonpolar region in the later case. In artificial lipid bilayers, midpolar regions are observed up to the level of acyl chain double bonds.

Solvation models and computational prediction of orientations of peptides and proteins in membranes

Membrane-associated peptides and proteins function in the highly heterogeneous environment of the lipid bilayer whose physico-chemical properties change non-monotonically along the bilayer normal. To simulate insertion of peptides and proteins into membranes and correctly reproduce the energetics of this process, an adequate solvation model and physically realistic representation of the lipid bilayer should be employed. We present a brief overview of the existing solvation models and their application for prediction of binding affinities and orientations of proteins in membranes. Particular emphasis is placed on the recently proposed PPM method, the corresponding web server, and the OPM database that were designed for positioning in membranes of integral and peripheral proteins with known three-dimensional structures.

Modulation of opioid receptor ligand affinity and efficacy using active and inactive state receptor models

Mu opioid receptor (MOR) agonists are widely used for the treatment of pain; however, chronic use results in the development of tolerance and dependence. It has been demonstrated that coadministration of a MOR agonist with a delta opioid receptor (DOR) antagonist maintains the analgesia associated with MOR agonists, but with reduced negative side-effects. Using our newly refined opioid receptor models for structure-based ligand design, we have synthesized several pentapeptides with tailored affinity and efficacy profiles. In particular, we have obtained pentapeptides 8, Tyr-c(S-S)[DCys-1Nal-Nle-Cys]NH(2), and 12, Tyr-c(S-S)[DCys-1Nal-Nle-Cys]OH, which demonstrates high affinity and full agonist behavior at MOR, high affinity but very low efficacy for DOR, and minimal affinity for the kappa opioid receptor (KOR). Functional properties of these peptides as MOR agonists/DOR antagonists lacking undesired KOR activity make them promising candidates for future in vivo studies of MOR/DOR interactions. Subtle structural variation of 12, by substituting D-Cys(5) for L-Cys(5), generated analog 13, which maintains low nanomolar MOR and DOR affinity, but which displays no efficacy at either receptor. These results demonstrate the power and utility of accurate receptor models for structure-based ligand design, as well as the profound sensitivity of ligand function on its structure.

Rescue of misrouted GnRHR mutants reveals its constitutive activity

G protein-coupled receptors (GPCR) play central roles in almost all physiological functions, and mutations in GPCR are responsible for over 30 hereditary diseases associated with loss or gain of receptor function. Gain of function mutants are frequently described as having constitutive activity (CA), that is, they activate effectors in the absence of agonist occupancy. Although many GPCR have mutants with CA, the GnRH receptor (GnRHR) was not, until 2010, associated with any CA mutants. The explanation for the failure to observe CA appears to be that the quality control system of the cell recognizes CA mutants of GnRHR as misfolded and retains them in the endoplasmic reticulum. In the present study, we identified several human (h)GnRHR mutants with substitutions in transmembrane helix 6 (F(272)K, F(272)Q, Y(284)F, C(279)A, and C(279)S) that demonstrate varying levels of CA after being rescued by pharmacoperones from different chemical classes and/or deletion of residue K(191), a modification that increases trafficking to the plasma membrane. The movement of the mutants from the endoplasmic reticulum (unrescued) to the plasma membrane (after rescue) is supported by confocal microscopy. Judging from the receptor-stimulated inositol phosphate production, mutants F(272)K and F(272)Q, after rescue, display the largest level of CA, an amount that is comparable with agonist-stimulated activation. Because mutations in other GPCR are, like the hGnRHR, scrutinized by the quality control system, this general approach may reveal CA in receptor mutants from other systems. A computer model of the hGnRHR and these mutants was used to evaluate the conformation associated with CA.

Development and in vitro characterization of a novel bifunctional μ-agonist/δ-antagonist opioid tetrapeptide

The development of tolerance to and dependence on opioid analgesics greatly reduces their long-term usefulness. Previous studies have demonstrated that co-administration of a μ-opioid receptor (MOR) agonist and δ-opioid receptor (DOR) antagonist can decrease MOR agonist-induced tolerance and dependence development after chronic exposure. Clinically, a single ligand displaying multiple efficacies (e.g., MOR agonism concurrently with DOR antagonism) would be of increased value over two drugs administered simultaneously. Guided by modeling of receptor-ligand complexes we have developed a series of potent non-selective opioid tetrapeptides that have differing efficacy at MOR and DOR. In particular, our lead peptide (KSK-103) binds with equal affinity to MOR and DOR but acts as a MOR agonist with similar efficacy but greater potency than morphine and a DOR antagonist in cellular assays measuring both G protein stimulation and adenylyl cyclase inhibition.

OPM database and PPM web server: resources for positioning of proteins in membranes

The Orientations of Proteins in Membranes (OPM) database is a curated web resource that provides spatial positions of membrane-bound peptides and proteins of known three-dimensional structure in the lipid bilayer, together with their structural classification, topology and intracellular localization. OPM currently contains more than 1200 transmembrane and peripheral proteins and peptides from approximately 350 organisms that represent approximately 3800 Protein Data Bank entries. Proteins are classified into classes, superfamilies and families and assigned to 21 distinct membrane types. Spatial positions of proteins with respect to the lipid bilayer are optimized by the PPM 2.0 method that accounts for the hydrophobic, hydrogen bonding and electrostatic interactions of the proteins with the anisotropic water-lipid environment described by the dielectric constant and hydrogen-bonding profiles. The OPM database is freely accessible at http://opm.phar.umich.edu. Data can be sorted, searched or retrieved using the hierarchical classification, source organism, localization in different types of membranes. The database offers downloadable coordinates of proteins and peptides with membrane boundaries. A gallery of protein images and several visualization tools are provided. The database is supplemented by the PPM server (http://opm.phar.umich.edu/server.php) which can be used for calculating spatial positions in membranes of newly determined proteins structures or theoretical models.

Salt bridges overlapping the gonadotropin-releasing hormone receptor agonist binding site reveal a coincidence detector for G protein-coupled receptor activation

G protein-coupled receptors (GPCRs) play central roles in most physiological functions, and mutations in them cause heritable diseases. Whereas crystal structures provide details about the structure of GPCRs, there is little information that identifies structural features that permit receptors to pass the cellular quality control system or are involved in transition from the ground state to the ligand-activated state. The gonadotropin-releasing hormone receptor (GnRHR), because of its small size among GPCRs, is amenable to molecular biological approaches and to computer modeling. These techniques and interspecies comparisons are used to identify structural features that are important for both intracellular trafficking and GnRHR activation yet distinguish between these processes. Our model features two salt (Arg(38)-Asp(98) and Glu(90)-Lys(121)) and two disulfide (Cys(14)-Cys(200) and Cys(114)-Cys(196)) bridges, all of which are required for the human GnRHR to traffic to the plasma membrane. This study reveals that both constitutive and ligand-induced activation are associated with a "coincidence detector" that occurs when an agonist binds. The observed constitutive activation of receptors lacking Glu(90)-Lys(121), but not Arg(38)-Asp(98) ionic bridge, suggests that the role of the former connection is holding the receptor in the inactive conformation. Both the aromatic ring and hydroxyl group of Tyr(284) and the hydrogen bonding of Ser(217) are important for efficient receptor activation. Our modeling results, supported by the observed influence of Lys(191) from extracellular loop 2 (EL2) and a four-residue motif surrounding this loop on ligand binding and receptor activation, suggest that the positioning of EL2 within the seven-α-helical bundle regulates receptor stability, proper trafficking, and function.

Anisotropic solvent model of the lipid bilayer. 2. Energetics of insertion of small molecules, peptides, and proteins in membranes

A new computational approach to calculating binding energies and spatial positions of small molecules, peptides, and proteins in the lipid bilayer has been developed. The method combines an anisotropic solvent representation of the lipid bilayer and universal solvation model, which predicts transfer energies of molecules from water to an arbitrary medium with defined polarity properties. The universal solvation model accounts for hydrophobic, van der Waals, hydrogen-bonding, and electrostatic solute-solvent interactions. The lipid bilayer is represented as a fluid anisotropic environment described by profiles of dielectric constant (ε), solvatochromic dipolarity parameter (π*), and hydrogen bonding acidity and basicity parameters (α and β). The polarity profiles were calculated using published distributions of quasi-molecular segments of lipids determined by neutron and X-ray scattering for DOPC bilayer and spin-labeling data that define concentration of water in the lipid acyl chain region. The model also accounts for the preferential solvation of charges and polar groups by water and includes the effect of the hydrophobic mismatch for transmembrane proteins. The method was tested on calculations of binding energies and preferential positions in membranes for small-molecules, peptides and peripheral membrane proteins that have been experimentally studied. The new theoretical approach was implemented in a new version (2.0) of our PPM program and applied for the large-scale calculations of spatial positions in membranes of more than 1000 peripheral and integral proteins. The results of calculations are deposited in the updated OPM database ( http://opm.phar.umich.edu ).

Pharmacological chaperones restore function to MC4R mutants responsible for severe early-onset obesity

Heterozygous null mutations in the melanocortin-4 receptor (MC4R) cause early-onset obesity in humans, indicating that metabolic homeostasis is sensitive to quantitative variation in MC4R function. Most of the obesity-causing MC4R mutations functionally characterized so far lead to intracellular retention of receptors by the cell's quality control system. Thus, recovering cell surface expression of mutant MC4Rs could have a beneficial therapeutic value. We tested a pharmacological chaperone approach to restore cell surface expression and function of 10 different mutant forms of human melanocortin-4 receptor found in obese patients. Five cell-permeant MC4R-selective ligands were tested and displayed pharmacological chaperone activities, restoring cell surface targeting and function of the receptors with distinct efficacy profiles for the different mutations. Such mutation-specific efficacies suggested a structure-activity relationship between compounds and mutant receptor conformations that may open a path toward personalized therapy. In addition, one of the five pharmacological chaperones restored function to most of the mutant receptors tested. Combined with its ability to reach the central nervous system and its selectivity for the MC4R, this pharmacological chaperone may represent a candidate for the development of a targeted therapy suitable for a large subset of patients with MC4R-deficient obesity.

Pentapeptides displaying mu opioid receptor agonist and delta opioid receptor partial agonist/antagonist properties

Chronic use of mu-opioid agonists has been shown to cause neurochemical adaptations resulting in tolerance and dependence. While the analgesic effects of these drugs are mediated by mu-opioid receptors (MOR), several studies have shown that antagonism or knockdown of delta-opioid receptors (DOR) can lessen or prevent development of tolerance and dependence. On the basis of computational modeling of putative active and inactive conformations of MOR and DOR, we have synthesized a series of pentapeptides with the goal of developing a MOR agonist/DOR antagonist peptide with similar affinity at both receptors as a tool to probe functional opioid receptor interaction(s). The eight resulting naphthylalanine-substituted cyclic pentapeptides displayed variable mixed-efficacy profiles. The most promising peptide (9; Tyr-c(S-CH(2)-S)[D-Cys-Phe-2-Nal-Cys]NH(2)) displayed a MOR agonist and DOR partial agonist/antagonist profile and bound with equipotent affinity (K(i) approximately 0.5 nM) to both receptors, but also showed kappa opioid receptor (KOR) agonist activity.

Functional characterization and structural modeling of obesity associated mutations in the melanocortin 4 receptor

Mutations in the melanocortin 4 receptor (MC4R) gene are the most common known cause of monogenic human obesity. The MC4R gene was sequenced in 2000 subjects with severe early-onset obesity. We detected seven different nonsense and 19 nonsynonymous mutations in a total of 94 probands, some of which have been reported previously by others. We functionally characterized the 11 novel obesity associated missense mutations. Seven of these mutants (L54P, E61K, I69T, S136P, M161T, T162I, and I269N) showed impaired cell surface trafficking, reduced level of maximal binding of the radioligand [125I]NDP-MSH, and reduced ability to generate cAMP in response to ligand. Four mutant MC4Rs (G55V, G55D, S136F, and A303T) displayed cell surface expression and agonist binding similar to the wild-type receptor but showed impaired cAMP production, suggesting that these residues are likely to be critical for conformational rearrangement essential for receptor activation. Homology modeling of these mutants using a model of MC4R based on the crystal structure of the beta2-adrenoreceptor was used to provide insights into the possible structural basis for receptor dysfunction. Transmembrane (TM) domains 1, 3, 6, 7, and peripheral helix 8 appear to participate in the agonist-induced conformational rearrangement necessary for coupling of ligand binding to signaling. We conclude that G55V, G55D, S136F, and A303T mutations are likely to strengthen helix-helix interactions between TM1 and TM2, TM3 and TM6, and TM7 and helix 8, respectively, preventing relative movement of these helices during receptor activation. The combination of functional studies and structural modeling of naturally occurring pathogenic mutations in MC4R can provide valuable information regarding the molecular mechanism of MC4R activation and its dysfunction in human disease.

Melanocortin tetrapeptide Ac-His-DPhe-Arg-Trp-NH2 modified at the para position of the benzyl side chain (DPhe): importance for mouse melanocortin-3 receptor agonist versus antagonist activity

The melanocortin-3 and -4 receptors (MC3R, MC4R) have been implicated in energy homeostasis and obesity. Whereas the physiological role of the MC4R is extensively studied, little is known about the MC3R. One caveat is the limited availability of ligands that are selective for the MC3R. Previous studies identified Ac-His-DPhe(p-I)-Arg-Trp-NH 2, which possessed partial agonist/antagonist pharmacology at the mMC3R while retaining full nanomolar agonist pharmacology at the mMC4R. These data allowed for the hypothesis that the DPhe position in melanocortin tetrapeptides can be used to examine ligand side-chain determinants important for differentiation of mMC3R agonist versus antagonist activity. A series of 15 DPhe (7) modified Ac-His-DPhe (7)-Arg-Trp-NH 2 tetrapeptides has been synthesized and pharmacologically characterized. Most notable results include the identification of modifications that resulted in potent antagonists/partial agonists at the mMC3R and full, potent agonists at the mMC4R. These SAR studies provide experimental evidence that the molecular mechanism of antagonism at the mMC3R differentiates this subtype from the mMC4R.

Novel peptide ligands of RGS4 from a focused one-bead, one-compound library

Regulators of G protein signaling accelerate GTP hydrolysis by G alpha subunits and profoundly inhibit signaling by G protein-coupled receptors. The distinct expression patterns and pathophysiologic regulation of regulators of G protein signaling proteins suggest that inhibitors may have therapeutic potential. We previously reported the design, mechanistic evaluation, and structure-activity relationships of a disulfide-containing cyclic peptide inhibitor of RGS4, YJ34 (Ac-Val-Lys-c[Cys-Thr-Gly-Ile-Cys]-Glu-NH(2), S-S) (Roof et al., Chem Biol Drug Des, 67, 2006, 266). Using a focused one-bead, one-compound peptide library that contains features known to be necessary for the activity of YJ34, we now identify peptides that bind to RGS4. Six peptides showed confirmed binding to RGS4 by flow cytometry. Two analogs of peptide 2 (Gly-Thr-c[Cys-Phe-Gly-Thr-Cys]-Trp-NH(2), S-S with a free or acetylated N-terminus) inhibited RGS4-stimulated G alpha(o) GTPase activity at 25-50 microM. They selectively inhibit RGS4 but not RGS7, RGS16, and RGS19. Their inhibition of RGS4 does not depend on cysteine-modification of RGS4, as they do not lose activity when all cysteines are removed from RGS4. Peptide 2 has been modeled to fit in the same binding pocket predicted for YJ34 but in the reverse orientation.

The role of hydrophobic interactions in positioning of peripheral proteins in membranes

BACKGROUND

Three-dimensional (3D) structures of numerous peripheral membrane proteins have been determined. Biological activity, stability, and conformations of these proteins depend on their spatial positions with respect to the lipid bilayer. However, these positions are usually undetermined.

RESULTS

We report the first large-scale computational study of monotopic/peripheral proteins with known 3D structures. The optimal translational and rotational positions of 476 proteins are determined by minimizing energy of protein transfer from water to the lipid bilayer, which is approximated by a hydrocarbon slab with a decadiene-like polarity and interfacial regions characterized by water-permeation profiles. Predicted membrane-binding sites, protein tilt angles and membrane penetration depths are consistent with spin-labeling, chemical modification, fluorescence, NMR, mutagenesis, and other experimental studies of 53 peripheral proteins and peptides. Experimental membrane binding affinities of peripheral proteins were reproduced in cases that did not involve a helix-coil transition, specific binding of lipids, or a predominantly electrostatic association. Coordinates of all examined peripheral proteins and peptides with the calculated hydrophobic membrane boundaries, subcellular localization, topology, structural classification, and experimental references are available through the Orientations of Proteins in Membranes (OPM) database.

CONCLUSION

Positions of diverse peripheral proteins and peptides in the lipid bilayer can be accurately predicted using their 3D structures that represent a proper membrane-bound conformation and oligomeric state, and have membrane binding elements present. The success of the implicit solvation model suggests that hydrophobic interactions are usually sufficient to determine the spatial position of a protein in the membrane, even when electrostatic interactions or specific binding of lipids are substantial. Our results demonstrate that most peripheral proteins not only interact with the membrane surface, but penetrate through the interfacial region and reach the hydrocarbon interior, which is consistent with published experimental studies.

Peptide and small molecules rescue the functional activity and agonist potency of dysfunctional human melanocortin-4 receptor polymorphisms

The melanocortin pathway, specifically the melanocortin-4 receptor and the cognate endogenous agonist and antagonist ligands, have been strongly implicated in the regulation of energy homeostasis and satiety. Genetic studies of morbidly obese human patients and normal weight control patients have resulted in the discovery of over 70 human melanocortin-4 receptor (MC4R) polymorphisms observed as both heterozygous and homozygous forms. A number of laboratories have been studying these hMC4R polymorphisms attempting to understand the molecular mechanism(s) that might explain the obese human phenotype. Herein, we have studied 13 polymorphic hMC4Rs that have been identified to possess statistically significant decreased endogenous agonist potency with synthetic peptides and small molecules attempting to identify ligands that can pharmacologically rescue the hMC4R polymorphic agonist response. The ligands examined in this study include NDP-MSH, MTII, Ac-His-DPhe-Arg-Trp-NH2 (JRH887-9), Ac-Anc-DPhe-Arg-Trp-NH2 (amino-2-naphtylcarboxylic acid, Anc, JRH420-12), Ac-His-(pI)DPhe-Arg-Trp-NH2 (JRH322-18), chimeric AGRP-melanocortin based ligands (Tyr-c[Cys-His-DPhe-Arg-Trp-Asn-Ala-Phe-Cys]-Tyr-NH2, AMW3-130 and Ac-mini-(His-DPhe-Arg-Trp)-hAGRP-NH2, AMW3-106), and the small molecules JB25 and THIQ. The hMC4R polymorphisms included in this study are S58C, N97D, I102S, L106P, S127L, T150I, R165Q, R165W, L250Q, G252S, C271Y, Y287Stop, and I301T. These studies resulted in the NDP-MSH, MTII, AMW3-130, THIQ, and AMW3-106 ligands possessing nanomolar to subnanomolar agonist potency at the hMC4R polymorphisms examined in this study. Thus, these ligands could generically rescue the potency and stimulatory response of the abnormally functioning hMC4Rs studied and may provide tools to further clarify the molecular mechanism(s) involving these receptor modifications.

Molecular mechanism of the constitutive activation of the L250Q human melanocortin-4 receptor polymorphism

The Melanocortin-4 Receptor is a G-protein coupled receptor that has been physiologically linked to participate in the regulation of energy homeostasis. The Melanocortin-4 Receptor is stimulated by endogenous melanocortin agonists derived from the pro-opiomelanocortin gene transcript and antagonized by the endogenous antagonist agouti-related protein. Central administration of melanocortin agonists has been demonstrated to decrease food intake and conversely, treatment with antagonists resulted in increased food intake. Deletion of the Melanocortin-4 Receptor gene from the mouse genome results in an obese and hyperphagic phenotype. Polymorphisms of the human Melanocortin-4-Receptor have been found in severely obese individuals, suggesting that Melanocortin-4 Receptor malfunction might be involved in human obesity and obesity-associated diabetes. Herein, we have performed experiments to understand the molecular mechanisms associated with the L250Q human Melanocortin-4-Receptor polymorphism discovered in an extremely obese woman. This L250Q human Melanocortin-4-Receptor has been pharmacologically characterized to result in a constitutively active receptor. The fact that a constitutively active human Melanocortin-4-Receptor mutation was found in an obese person is a physiologic contradiction, as chronic activation of the human Melanocortin-4-Receptor and subsequently high cyclic adenosine monophosphate levels should theoretically result in a normal or lean phenotype. In this study, we demonstrated that agouti-related protein acts as an inverse agonist at this constitutively active receptor, and we propose a mechanism by which agouti-related protein might contribute to the obese phenotype in the L250Q patient. In addition, using receptor mutagenesis, pharmacology, and computer modeling approaches, we investigated the molecular mechanism by which modification of the L250 residue results in constitutive activation of the human Melanocortin-4-Receptor.

Positioning of proteins in membranes: a computational approach

A new computational approach has been developed to determine the spatial arrangement of proteins in membranes by minimizing their transfer energies from water to the lipid bilayer. The membrane hydrocarbon core was approximated as a planar slab of adjustable thickness with decadiene-like interior and interfacial polarity profiles derived from published EPR studies. Applicability and accuracy of the method was verified for a set of 24 transmembrane proteins whose orientations in membranes have been studied by spin-labeling, chemical modification, fluorescence, ATR FTIR, NMR, cryo-microscopy, and neutron diffraction. Subsequently, the optimal rotational and translational positions were calculated for 109 transmembrane, five integral monotopic and 27 peripheral protein complexes with known 3D structures. This method can reliably distinguish transmembrane and integral monotopic proteins from water-soluble proteins based on their transfer energies and membrane penetration depths. The accuracies of calculated hydrophobic thicknesses and tilt angles were approximately 1 A and 2 degrees, respectively, judging from their deviations in different crystal forms of the same proteins. The hydrophobic thicknesses of transmembrane proteins ranged from 21.1 to 43.8 A depending on the type of biological membrane, while their tilt angles with respect to the bilayer normal varied from zero in symmetric complexes to 26 degrees in asymmetric structures. Calculated hydrophobic boundaries of proteins are located approximately 5 A lower than lipid phosphates and correspond to the zero membrane depth parameter of spin-labeled residues. Coordinates of all studied proteins with their membrane boundaries can be found in the Orientations of Proteins in Membranes (OPM) database:http://opm.phar.umich.edu/.

Design of high affinity cyclic pentapeptide ligands for kappa-opioid receptors

Using results from our previously reported cyclic opioid peptide series and reliable models for mu-, delta-, and kappa-opioid receptors (MOR, DOR, and KOR, respectively) and their complexes with peptide ligands, we have designed and synthesized a series of cyclic pentapeptides of structure Tyr-C[D-Cys-Phe-Phe-X]-NH2, cyclized via disulfide, methylene, or ethylene dithioethers, and where X = D- or L-Cys; or D- or L-penicillamine (Pen; beta,beta-dimethylcysteine). Determination of binding affinities to MOR, DOR, and KOR revealed that members of this series with X = D- or L-Cys display KOR affinities in the low nanomolar range, demonstrating that a 'DPDPE-like' tetrapeptide scaffold is suitable not only for DOR and MOR ligands, but also for KOR ligands. The cyclic pentapeptides reported here are not, however, selective for KOR, rather they display significant selectivity and high affinity for MOR. Indeed, peptide 8, Tyr-C[D-Cys-Phe-Phe-Cys]-NH2-cyclized via a methylene dithioether, shows picomolar binding affinity for MOR ( = 16 pm) with more than 100-fold selectivity for MOR vs. DOR or KOR, and may be of interest as a high affinity, high selectivity MOR ligand. Nonetheless, the high affinity KOR peptides in this series represent excellent leads for the development of structurally related, selective KOR ligands designed to exploit structurally specific features of KOR, MOR, and DOR.

OPM: Orientations of Proteins in Membranes database

SUMMARY

The Orientations of Proteins in Membranes (OPM) database provides a collection of transmembrane, monotopic and peripheral proteins from the Protein Data Bank whose spatial arrangements in the lipid bilayer have been calculated theoretically and compared with experimental data. The database allows analysis, sorting and searching of membrane proteins based on their structural classification, species, destination membrane, numbers of transmembrane segments and subunits, numbers of secondary structures and the calculated hydrophobic thickness or tilt angle with respect to the bilayer normal. All coordinate files with the calculated membrane boundaries are available for downloading.

AVAILABILITY

http://opm.phar.umich.edu.