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

Development of enhanced conformational sampling methods to probe the activation landscape of GPCRs

G protein-coupled receptors (GPCRs) make up the largest superfamily of integral membrane proteins and play critical signal transduction roles in many physiological processes. Developments in molecular biology, biophysical, biochemical, pharmacological, and computational techniques aimed at these important therapeutic targets are beginning to provide unprecedented details on the structural as well as functional basis of their pleiotropic signaling mediated by G proteins, β arrestins, and other transducers. This pleiotropy presents a pharmacological challenge as the same ligand-receptor interaction can cause a therapeutic effect as well as an undesirable on-target side-effect through different downstream pathways. GPCRs don't function as simple binary on-off switches but as finely tuned shape-shifting machines described by conformational ensembles, where unique subsets of conformations may be responsible for specific signaling cascades. X-ray crystallography and more recently cryo-electron microscopy are providing snapshots of some of these functionally-important receptor conformations bound to ligands and/or transducers, which are being utilized by computational methods to describe the dynamic conformational energy landscape of GPCRs. In this chapter, we review the progress in computational conformational sampling methods based on molecular dynamics and discrete sampling approaches that have been successful in complementing biophysical and biochemical studies on these receptors in terms of their activation mechanisms, allosteric effects, actions of biased ligands, and effects of pathological mutations. Some of the sampled simulation time scales are beginning to approach receptor activation time scales. The list of conformational sampling methods and example uses discussed is not exhaustive but includes representative examples that have pushed the limits of classical molecular dynamics and discrete sampling methods to describe the activation energy landscape of GPCRs.

Transmembrane Helix Integrity versus Fraying To Expose Hydrogen Bonds at a Membrane-Water Interface

Transmembrane helices dominate the landscape for many membrane proteins. Often flanked by interfacial aromatic residues, these transmembrane helices also contain loops and interhelix segments, which could help in stabilizing a transmembrane orientation. Using ²H nuclear magnetic resonance spectroscopy to monitor bilayer-incorporated model GWALP23 family peptides, we address systematically the issue of helix fraying in relation to the dynamics and orientation of highly similar individual transmembrane helices. We inserted aromatic (Phe, Trp, Tyr, and His) or non-aromatic residues (Ala and Gly) into positions 4 and 5 adjacent to a core transmembrane helix to examine the side-chain dependency of the transmembrane orientation, dynamics, and helix integrity (extent and location of unraveling). Incorporation of [²H]alanine labels enables one to assess the helicity of the core sequence and the peptide termini. For most of the helices, we observed substantial unwinding involving at least three residues at both ends. For the unique case of histidine at positions 4 and 5, an extended N-terminal unwinding was observed up to residue 7. For further investigation of the onset of fraying, we employed A^4,5GWALP23 with ²H labels at residues 4 and 5 and found that the number of terminal residues involved in the unwinding depends on bilayer thicknesses and helps to govern the helix dynamics. The combined results enable us to compare and contrast the extent of fraying for each related helix, as reflected by the deviation of experimental ²H quadrupolar splitting magnitudes of juxta-terminal alanines A3 and A21 from those represented by an ideal helix geometry.

Novel and recurrent variants in AVPR2 in 19 families with X-linked congenital nephrogenic diabetes insipidus

Congenital nephrogenic diabetes insipidus (CNDI) is characterized by the reduced ability of renal collecting duct cells to reabsorb water in response to the antidiuretic effect of vasopressin. Chronic polyuria and polydipsia are the hallmarks of the disease. Approximately 90% of all patients with CNDI have X-linked inherited disease caused by variants in the arginine vasopressin receptor 2 (AVPR2) gene. We present genetic findings in 34 individuals from 19 kindreds including one or more family members with CNDI. Coding regions of AVPR2 were sequenced bi-directionally. We identified eight novel disease-causing variants in AVPR2, p.Arg68Alafs*124, p.Ser171Arg, p.Gln174Pro, p.Trp200Arg, p.Gly201Cys, p.Gly220Arg, p.Val226Glu, and p.Gln291Pro in nine kindreds. In all three families with more than one affected individual, the novel variants segregated with the disease. We also identified eight recurrent disease-causing variants, p.Val88Met, p.Leu111Valfs*80, p.Arg113Trp, p.Tyr124*, p.Ser167Leu, p.Thr207Asn, p.Arg247Alafs*12, and p.Arg337* in ten kindreds. Our findings contribute to the growing list of AVPR2 variants causing X-linked CNDI.

CONCLUSION

Being a rapid diagnostic tool for CNDI, direct sequencing of AVPR2 should be encouraged in newborns with familial predisposition to CNDI. What is Known: • Disease-causing variants in AVPR2 cause X-linked congenital nephrogenic diabetes insipidus (CNDI). • DNA sequencing of AVPR2 is rapid, facilitates differential diagnosis, early intervention, and genetic diagnosis thus reducing morbidity in CNDI. What is New: • We identified eight novel disease-causing variants in AVPR2: p.Arg68Alafs*124, p.Ser171Arg, p.Gln174Pro, p.Trp200Arg, p.Gly201Cys, p.Gly220Arg, p.Val226Glu, and p.Gln291Pro, thereby adding to the growing list of AVPR2 disease-causing variants and emphasizing the importance of genetic testing in CNDI.

Optimization and 13CH3 methionine labeling of a signaling competent neurotensin receptor 1 variant for NMR studies

Neurotensin is a 13-residue peptide that acts as a neuromodulator of classical neurotransmitters such as dopamine and glutamate in the mammalian central nervous system, mainly by activating the G protein-coupled receptor (GPCR), neurotensin receptor 1 (NTS₁). Agonist binding to GPCRs shifts the conformational equilibrium of the transmembrane helices towards distinct, thermodynamically favorable conformations that favor effector protein interactions and promotes cell signaling. The introduction of site specific labels for NMR spectroscopy has proven useful for investigating this dynamic process, but the low expression levels and poor stability of GPCRs is a hindrance to solution NMR experiments. Several thermostabilized mutants of NTS₁ have been engineered to circumvent this, with the crystal structures of four of these published. The conformational dynamics of NTS₁ however, has not been thoroughly investigated with NMR. It is generally accepted that stabilized GPCRs exhibit attenuated signaling, thus we thoroughly characterized the signaling characteristics of several thermostabilized NTS₁ variants to identify an optimal variant for protein NMR studies. A variant termed enNTS₁ exhibited the best combination of signaling capability and stability upon solubilization with detergents. enNTS₁ was subsequently labeled with ¹³CH₃-methionine in E. coli and purified to homogeneity in the absence of bound ligands. Using solution NMR spectroscopy we observed several well dispersed ¹³CH₃-methionine resonances, many of which exhibited chemical shift changes upon the addition of the high affinity agonist peptide, NT8-13. Thus, enNTS₁ represents a novel tool for investigating ligand induced conformational changes in NTS₁ to gain insights into the molecular mechanisms underlying neurotensin signaling.

Structure-Based Sequence Alignment of the Transmembrane Domains of All Human GPCRs: Phylogenetic, Structural and Functional Implications

The understanding of G-protein coupled receptors (GPCRs) is undergoing a revolution due to increased information about their signaling and the experimental determination of structures for more than 25 receptors. The availability of at least one receptor structure for each of the GPCR classes, well separated in sequence space, enables an integrated superfamily-wide analysis to identify signatures involving the role of conserved residues, conserved contacts, and downstream signaling in the context of receptor structures. In this study, we align the transmembrane (TM) domains of all experimental GPCR structures to maximize the conserved inter-helical contacts. The resulting superfamily-wide GpcR Sequence-Structure (GRoSS) alignment of the TM domains for all human GPCR sequences is sufficient to generate a phylogenetic tree that correctly distinguishes all different GPCR classes, suggesting that the class-level differences in the GPCR superfamily are encoded at least partly in the TM domains. The inter-helical contacts conserved across all GPCR classes describe the evolutionarily conserved GPCR structural fold. The corresponding structural alignment of the inactive and active conformations, available for a few GPCRs, identifies activation hot-spot residues in the TM domains that get rewired upon activation. Many GPCR mutations, known to alter receptor signaling and cause disease, are located at these conserved contact and activation hot-spot residue positions. The GRoSS alignment places the chemosensory receptor subfamilies for bitter taste (TAS2R) and pheromones (Vomeronasal, VN1R) in the rhodopsin family, known to contain the chemosensory olfactory receptor subfamily. The GRoSS alignment also enables the quantification of the structural variability in the TM regions of experimental structures, useful for homology modeling and structure prediction of receptors. Furthermore, this alignment identifies structurally and functionally important residues in all human GPCRs. These residues can be used to make testable hypotheses about the structural basis of receptor function and about the molecular basis of disease-associated single nucleotide polymorphisms.

G-protein coupled receptors (GPCRs) are a large superfamily of integral membrane proteins that share a characteristic 7 transmembrane helix fold. They detect various molecules outside of the cell and signal their presence to the inside of the cell. At least half of the 800 human GPCRs are potential drug targets, so understanding their structure and function is critical. Experimental structures are now available for at least one receptor from each GPCR class. The structure of the 7 helix fold is highly conserved even for receptors with very low sequence similarity. We analyze the available experimental structures and compare the common inter-helical contacts. Our analysis leads to a unified sequence-structure alignment of the GPCR superfamily that can then be used as the starting point for structure prediction of all other GPCRs. A key result of our analysis is a list of conserved contact residues and activation “hot-spots” residues that are critical for GPCR folding and function. We propose that mutations and natural variants of amino acids at these locations in the GPCRs can dramatically influence their activation state and alter intracellular signaling. This provides hypotheses for the molecular mechanisms underlying disease causing mutants for any GPCR.

Predicted 3D structures of olfactory receptors with details of odorant binding to OR1G1

Olfactory receptors (ORs) are responsible for mediating the sense of smell; they allow humans to recognize an enormous number of odors but the connection between binding and perception is not known. We predict the ensemble of low energy structures for the human OR1G1 (hOR1G1) and also for six other diverse ORs, using the G protein-coupled receptor Ensemble of Structures in Membrane BiLayer Environment complete sampling method that samples 13 trillion different rotations and tilts using four different templates to predict the 24 structures likely to be important in binding and activation. Our predicted most stable structures of hOR1G1 have a salt-bridge between the conserved D3.49 and K6.30 in the D(E)RY region, that we expect to be associated with an inactive form. The hOR1G1 structure also has specific interaction in transmembrane domains (TMD) 3-6 (E3.39 and H6.40), which is likely an important conformational feature for all hORs because of the ~94 to 98 % conservation among all hOR sequences. Of the five ligands studied (nonanal, 9-decen-1-ol, 1-nonanol, camphor, and n-butanal), we find that the 4 expected to bind lead to similar binding energies with nonanol the strongest.

Ligand- and mutation-induced conformational selection in the CCR5 chemokine G protein-coupled receptor

We predicted the structural basis for pleiotropic signaling of the C-C chemokine type 5 (CCR5) G protein-coupled receptor (GPCR) by predicting the binding of several ligands to the lower-energy conformations of the CCR5 receptor and 11 mutants. For each case, we predicted the ∼ 20 most stable conformations for the receptor along with the binding sites for four anti-HIV ligands. We found that none of the ligands bind to the lowest-energy apo-receptor conformation. The three ligands with a similar pharmacophore (Maraviroc, PF-232798, and Aplaviroc) bind to a specific higher-energy receptor conformation whereas TAK-779 (with a different pharmacophore) binds to a different high-energy conformation. This result is in agreement with the very different binding-site profiles for these ligands obtained by us and others. The predicted Maraviroc binding site agrees with the recent structure of CCR5 receptor cocrystallized with Maraviroc. We performed 11 site-directed mutagenesis experiments to validate the predicted binding sites. Here, we independently predicted the lowest 10 mutant protein conformations for each of the 11 mutants and then docked the ligands to these lowest conformations. We found the predicted binding energies to be in excellent agreement with our mutagenesis experiments. These results show that, for GPCRs, each ligand can stabilize a different protein conformation, complicating the use of cocrystallized structures for ligand screening. Moreover, these results show that a single-point mutation in a GPCR can dramatically alter the available low-energy conformations, which in turn alters the binding site, potentially altering downstream signaling events. These studies validate the conformational selection paradigm for the pleiotropic function and structural plasticity of GPCRs.

Modelling three-dimensional protein structures for applications in drug design

A structural perspective of drug target and anti-target proteins, and their molecular interactions with biologically active molecules, largely advances many areas of drug discovery, including target validation, hit and lead finding and lead optimisation. In the absence of experimental 3D structures, protein structure prediction often offers a suitable alternative to facilitate structure-based studies. This review outlines recent methodical advances in homology modelling, with a focus on those techniques that necessitate consideration of ligand binding. In this context, model quality estimation deserves special attention because the accuracy and reliability of different structure prediction techniques vary considerably, and the quality of a model ultimately determines its usefulness for structure-based drug discovery. Examples of G-protein-coupled receptors (GPCRs) and ADMET-related proteins were selected to illustrate recent progress and current limitations of protein structure prediction. Basic guidelines for good modelling practice are also provided.

Recent developments in biased agonism

The classic paradigm of G protein-coupled receptor (GPCR) activation was based on the understanding that agonist binding to a receptor induces or stabilizes a conformational change to an 'active' conformation. In the past decade, however, it has been appreciated that ligands can induce distinct 'active' receptor conformations with unique downstream functional signaling profiles. Building on the initial recognition of the existence of such 'biased ligands', recent years have witnessed significant developments in several areas of GPCR biology. These include increased understanding of structural and biophysical mechanisms underlying biased agonism, improvements in characterization and quantification of ligand efficacy, as well as clinical development of these novel ligands. Here we review recent major developments in these areas over the past several years.

Structural dynamics and regulation of the mammalian SLC9A family of Na⁺/H⁺ exchangers

Mammalian Na⁺/H⁺ exchangers of the SLC9A family are widely expressed and involved in numerous essential physiological processes. Their primary function is to mediate the 1:1 exchange of Na⁺ for H⁺ across the membrane in which they reside, and they play central roles in regulation of body, cellular, and organellar pH. Their function is tightly regulated through mechanisms involving interactions with multiple protein and lipid-binding partners, phosphorylations, and other posttranslational modifications. Biochemical and mutational analyses indicate that the SLC9As have a short intracellular N-terminus, 12 transmembrane (TM) helices necessary and sufficient for ion transport, and a C-terminal cytoplasmic tail region with essential regulatory roles. No high-resolution structures of the SLC9As exist; however, models based on crystal structures of the bacterial NhaAs support the 12 TM organization and suggest that TMIV and XI may form a central part of the ion-translocation pathway, whereas pH sensing may involve TMII, TMIX, and several intracellular loops. Similar to most ion transporters studied, SLC9As likely exist as coupled dimers in the membrane, and this appears to be important for the well-studied cooperativity of H⁺ binding. The aim of this work is to summarize and critically discuss the currently available evidence on the structural dynamics, regulation, and binding partner interactions of SLC9As, focusing in particular on the most widely studied isoform, SLC9A1/NHE1. Further, novel bioinformatic and structural analyses are provided that to some extent challenge the existing paradigm on how ions are transported by mammalian SLC9As.

SuperBiHelix method for predicting the pleiotropic ensemble of G-protein-coupled receptor conformations

There is overwhelming evidence that G-protein-coupled receptors (GPCRs) exhibit several distinct low-energy conformations, each of which might favor binding to different ligands and/or lead to different downstream functions. Understanding the function of such proteins requires knowledge of the ensemble of low-energy configurations that might play a role in this pleiotropic functionality. We earlier reported the BiHelix method for efficiently sampling the (12)(7) = 35 million conformations resulting from 30° rotations about the axis (η) of all seven transmembrane helices (TMHs), showing that the experimental structure is reliably selected as the best conformation from this ensemble. However, various GPCRs differ sufficiently in the tilts of the TMHs that this method need not predict the optimum conformation starting from any other template. In this paper, we introduce the SuperBiHelix method in which the tilt angles (θ, ϕ) are optimized simultaneously with rotations (η) efficiently enough that it is practical and sufficient to sample (5 × 3 × 5)(7) = 13 trillion configurations. This method can correctly identify the optimum structure of a GPCR starting with the template from a different GPCR. We have validated this method by predicting known crystal structure conformations starting from the template of a different protein structure. We find that the SuperBiHelix conformational ensemble includes the higher energy conformations associated with the active protein in addition to those associated with the more stable inactive protein. This methodology was then applied to design and experimentally confirm structures of three mutants of the CB1 cannabinoid receptor associated with different functions.

How broadly tuned olfactory receptors equally recognize their agonists. Human OR1G1 as a test case

The molecular features that dominate the binding mode of agonists by a broadly tuned olfactory receptor are analyzed through a joint approach combining cell biology, calcium imaging, and molecular modeling. The odorant/receptor affinities, estimated through statistics accrued during molecular dynamics simulations, are in accordance with the experimental ranking. Although in many systems receptors recognize their target through a network of oriented interactions, such as H-bonding, the binding by broadly tuned olfactory receptors is dominated by non-polar terms. We show how such a feature allows chemicals belonging to different chemical families to similarly activate the receptors through compensations of interactions within the binding site.

Identification and characterization of an activating F229V substitution in the V2 vasopressin receptor in an infant with NSIAD

Gain-of-function mutations in the gene encoding the V2 vasopressin receptor (V2R) cause nephrogenic syndrome of inappropriate antidiuresis. To date, reported mutations lead to the substitution of arginine 137 by either a cysteine or leucine (R137C/L). Here, we describe a 3-month-old hyponatremic infant found to have a phenylalanine 229 to valine (F229V) substitution in V2R. Characterization of this substitution in vitro revealed that it leads to high constitutive activity of the receptor, compatible with spontaneous antidiuresis. In contrast to R137C/L mutant receptors, F229V receptors do not undergo spontaneous desensitization, which results in sustained, high basal activity. Notably, the V2R-selective inverse agonists tolvaptan and satavaptan completely silenced the constitutive signaling activity of the F229V mutant receptor, indicating that this substitution does not lock the receptor in an irreversible active state. Thus, inverse agonists might prove to be effective therapies for treating patients with this or other spontaneously activating mutations that do not lock the V2R in its active state. These results emphasize the importance of genetic testing and the functional characterization of mutant receptors for patients with nephrogenic syndrome of inappropriate antidiuresis because the results might inform treatment decisions.

Genetically encoded photo-cross-linkers map the binding site of an allosteric drug on a G protein-coupled receptor

G protein-coupled receptors (GPCRs) are dynamic membrane proteins that bind extracellular molecules to transduce signals. Although GPCRs represent the largest class of therapeutic targets, only a small percentage of their ligand-binding sites are precisely defined. Here we describe the novel application of targeted photo-cross-linking using unnatural amino acids to obtain structural information about the allosteric binding site of a small molecule drug, the CCR5-targeted HIV-1 co-receptor blocker maraviroc.