Stabilization of G protein-coupled receptors by point mutations

The path to the G protein-coupled receptor structural landscape: Major milestones and future directions

G protein-coupled receptors (GPCRs) play a crucial role in cell function by transducing signals from the extracellular environment to the inside of the cell. They mediate the effects of various stimuli, including hormones, neurotransmitters, ions, photons, food tastants and odorants, and are renowned drug targets. Advancements in structural biology techniques, including X-ray crystallography and cryo-electron microscopy (cryo-EM), have driven the elucidation of an increasing number of GPCR structures. These structures reveal novel features that shed light on receptor activation, dimerization and oligomerization, dichotomy between orthosteric and allosteric modulation, and the intricate interactions underlying signal transduction, providing insights into diverse ligand-binding modes and signalling pathways. However, a substantial portion of the GPCR repertoire and their activation states remain structurally unexplored. Future efforts should prioritize capturing the full structural diversity of GPCRs across multiple dimensions. To do so, the integration of structural biology with biophysical and computational techniques will be essential. We describe in this review the progress of nuclear magnetic resonance (NMR) to examine GPCR plasticity and conformational dynamics, of atomic force microscopy (AFM) to explore the spatial-temporal dynamics and kinetic aspects of GPCRs, and the recent breakthroughs in artificial intelligence for protein structure prediction to characterize the structures of the entire GPCRome. In summary, the journey through GPCR structural biology provided in this review illustrates how far we have come in decoding these essential proteins architecture and function. Looking ahead, integrating cutting-edge biophysics and computational tools offers a path to navigating the GPCR structural landscape, ultimately advancing GPCR-based applications.

The application of artificial intelligence to accelerate G protein-coupled receptor drug discovery

The application of artificial intelligence (AI) approaches to drug discovery for G protein-coupled receptors (GPCRs) is a rapidly expanding area. Artificial intelligence can be used at multiple stages during the drug discovery process, from aiding our understanding of the fundamental actions of GPCRs to the discovery of new ligand-GPCR interactions or the prediction of clinical responses. Here, we provide an overview of the concepts behind artificial intelligence, including the subfields of machine learning and deep learning. We summarise the published applications of artificial intelligence to different stages of the GPCR drug discovery process. Finally, we reflect on the benefits and limitations of artificial intelligence and share our vision for the exciting potential for further development of applications to aid GPCR drug discovery. In addition to making the drug discovery process "faster, smarter and cheaper," we anticipate that the application of artificial intelligence will create exciting new opportunities for GPCR drug discovery. LINKED ARTICLES: This article is part of a themed issue Therapeutic Targeting of G Protein-Coupled Receptors: hot topics from the Australasian Society of Clinical and Experimental Pharmacologists and Toxicologists 2021 Virtual Annual Scientific Meeting. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v181.14/issuetoc.

Structural and dynamic insights into the activation of the μ-opioid receptor by an allosteric modulator

G-protein-coupled receptors (GPCRs) play pivotal roles in various physiological processes. These receptors are activated to different extents by diverse orthosteric ligands and allosteric modulators. However, the mechanisms underlying these variations in signaling activity by allosteric modulators remain largely elusive. Here, we determine the three-dimensional structure of the μ-opioid receptor (MOR), a class A GPCR, in complex with the Gi protein and an allosteric modulator, BMS-986122, using cryogenic electron microscopy. Our results reveal that BMS-986122 binding induces changes in the map densities corresponding to R1673.50 and Y2545.58, key residues in the structural motifs conserved among class A GPCRs. Nuclear magnetic resonance analyses of MOR in the absence of the Gi protein reveal that BMS-986122 binding enhances the formation of the interaction between R1673.50 and Y2545.58, thus stabilizing the fully-activated conformation, where the intracellular half of TM6 is outward-shifted to allow for interaction with the Gi protein. These findings illuminate that allosteric modulators like BMS-986122 can potentiate receptor activation through alterations in the conformational dynamics in the core region of GPCRs. Together, our results demonstrate the regulatory mechanisms of GPCRs, providing insights into the rational development of therapeutics targeting GPCRs.

High-Temperature Tolerance Protein Engineering through Deep Evolution

Protein engineering aimed at increasing temperature tolerance through iterative mutagenesis and high-throughput screening is often labor-intensive. Here, we developed a deep evolution (DeepEvo) strategy to engineer protein high-temperature tolerance by generating and selecting functional sequences using deep learning models. Drawing inspiration from the concept of evolution, we constructed a high-temperature tolerance selector based on a protein language model, acting as selective pressure in the high-dimensional latent spaces of protein sequences to enrich those with high-temperature tolerance. Simultaneously, we developed a variant generator using a generative adversarial network to produce protein sequence variants containing the desired function. Afterward, the iterative process involving the generator and selector was executed to accumulate high-temperature tolerance traits. We experimentally tested this approach on the model protein glyceraldehyde 3-phosphate dehydrogenase, obtaining 8 variants with high-temperature tolerance from just 30 generated sequences, achieving a success rate of over 26%, demonstrating the high efficiency of DeepEvo in engineering protein high-temperature tolerance.

Scaling up Functional Analyses of the G Protein-Coupled Receptor Rhodopsin

Eukaryotic cells use G protein-coupled receptors (GPCRs) to convert external stimuli into internal signals to elicit cellular responses. However, how mutations in GPCR-coding genes affect GPCR activation and downstream signaling pathways remain poorly understood. Approaches such as deep mutational scanning show promise in investigations of GPCRs, but a high-throughput method to measure rhodopsin activation has yet to be achieved. Here, we scale up a fluorescent reporter assay in budding yeast that we engineered to study rhodopsin's light-activated signal transduction. Using this approach, we measured the mutational effects of over 1200 individual human rhodopsin mutants, generated by low-frequency random mutagenesis of the GPCR rhodopsin (RHO) gene. Analysis of the data in the context of rhodopsin's three-dimensional structure reveals that transmembrane helices are generally less tolerant to mutations compared to flanking helices that face the lipid bilayer, which suggest that mutational tolerance is contingent on both the local environment surrounding specific residues and the specific position of these residues in the protein structure. Comparison of functional scores from our screen to clinically identified rhodopsin disease variants found many pathogenic mutants to be loss of function. Lastly, functional scores from our assay were consistent with a complex counterion mechanism involved in ligand-binding and rhodopsin activation. Our results demonstrate that deep mutational scanning is possible for rhodopsin activation and can be an effective method for revealing properties of mutational tolerance that may be generalizable to other transmembrane proteins.

Molecular simulations of SSTR2 dynamics and interaction with ligands

The cyclic peptide hormone somatostatin regulates physiological processes involved in growth and metabolism, through its binding to G-protein coupled somatostatin receptors. The isoform 2 (SSTR2) is of particular relevance for the therapy of neuroendocrine tumours for which different analogues to somatostatin are currently in clinical use. We present an extensive and systematic computational study on the dynamics of SSTR2 in three different states: active agonist-bound, inactive antagonist-bound and apo inactive. We exploited the recent burst of SSTR2 experimental structures to perform μs-long multi-copy molecular dynamics simulations to sample conformational changes of the receptor and rationalize its binding to different ligands (the agonists somatostatin and octreotide, and the antagonist CYN154806). Our findings suggest that the apo form is more flexible compared to the holo ones, and confirm that the extracellular loop 2 closes upon the agonist octreotide but not upon the antagonist CYN154806. Based on interaction fingerprint analyses and free energy calculations, we found that all peptides similarly interact with residues buried into the binding pocket. Conversely, specific patterns of interactions are found with residues located in the external portion of the pocket, at the basis of the extracellular loops, particularly distinguishing the agonists from the antagonist. This study will help in the design of new somatostatin-based compounds for theranostics of neuroendocrine tumours.

Enhancing the stability and homogeneity of non-ionic polymer nanodiscs by tuning electrostatic interactions

The nanodisc technology is increasingly used for structural studies on membrane proteins and drug delivery. The development of synthetic polymer nanodiscs and the recent discovery of non-ionic inulin-based polymers have significantly broadened the scope of nanodiscs. While the lipid exchange and size flexibility properties of the self-assembled polymer-based nanodiscs are valuable for various applications, the non-ionic polymer nanodiscs are remarkably unique in that they enable the reconstitution of any protein, protein-protein complexes, or drugs irrespective of their charge. However, the non-ionic nature of the belt could influence the stability and size homogeneity of inulin-based polymer nanodiscs. In this study, we investigate the size stability and homogeneity of nanodiscs formed by non-ionic lipid-solubilizing polymers using different biophysical methods. Polymer nanodiscs containing zwitterionic DMPC and different ratios of DMPC:DMPG lipids were made using anionic SMA-EA or non-ionic pentyl-inulin polymers. Non-ionic polymer nanodiscs made using zwitterionic DMPC lipids produced a very broad elution profile on SEC due to their instability in the column, thus affecting sample monodispersity which was confirmed by DLS experiments that showed multiple peaks. However, the inclusion of anionic DMPG lipids improved the stability as observed from SEC and DLS profiles, which was further confirmed by TEM images. Whereas, anionic SMA-EA-based DMPC-nanodiscs showed excellent stability and size homogeneity when solubilizing zwitterionic lipids. The stability of DMPC:DMPG non-ionic polymer nanodiscs is attributed to the inter-nanodisc repulsion by the anionic-DMPG that prevents the uncontrolled collision and fusion of nanodiscs. Thus, the reported results demonstrate the use of electrostatic interactions to tune the solubility, stability, and size homogeneity of non-ionic polymer nanodiscs which are important features for enabling functional and atomic-resolution structural studies of membrane proteins, other lipid-binding molecules, and water-soluble biomolecules including cytosolic proteins, nucleic acids and metabolites.

Rapid One-Step Capturing of Native, Cell-Free Synthesized and Membrane-Embedded GLP-1R

G protein-coupled receptors (GPCRs) are of outstanding pharmacological interest as they are abundant in cell membranes where they perform diverse functions that are closely related to the vitality of cells. The analysis of GPCRs in natural membranes is laborious, as established methods are almost exclusively cell culture-based and only a few methods for immobilization in a natural membrane outside the cell are known. Within this study, we present a one-step, fast and robust immobilization strategy of the GPCR glucagon-like peptide 1 receptor (GLP-1R). GLP-1R was synthesized in eukaryotic lysates harboring endogenous endoplasmic reticulum-derived microsomes enabling the embedment of GLP-1R in a natural membrane. Interestingly, we found that these microsomes spontaneously adsorbed to magnetic Neutravidin beads thus providing immobilized membrane protein preparations which required no additional manipulation of the target receptor or its supporting membrane. The accessibility of the extracellular domain of membrane-embedded and bead-immobilized GLP-1R was demonstrated by bead-based enzyme-linked immunosorbent assay (ELISA) using GLP-1R-specific monoclonal antibodies. In addition, ligand binding of immobilized GLP-1R was verified in a radioligand binding assay. In summary, we present an easy and straightforward synthesis and immobilization methodology of an active GPCR which can be beneficial for studying membrane proteins in general.

Fluorescence-Detection Size-Exclusion Chromatography-Based Thermostability Assay for Membrane Proteins

Green fluorescent proteins (GFPs) have lightened up almost every aspect of biological research including protein sciences. In the field of membrane protein structural biology, GFPs have been used widely to monitor membrane protein localization, expression level, the purification process and yield, and the stability inside the cells and in the test tube. Of particular interest is the fluorescence-detector size-exclusion chromatography-based thermostability assay (FSEC-TS). By simple heating and FSEC, the generally applicable method allows rapid assessment of the thermostability of GFP-fused membrane proteins without purification. Here we describe the experimental details and some typical results for the FSEC-TS method.

New insights into the structure and function of class B1 GPCRs

G protein-coupled receptors (GPCRs) are the largest family of cell surface receptors. Class B1 GPCRs constitute a subfamily of 15 receptors that characteristically contain large extracellular domains (ECDs) and respond to long polypeptide hormones. Class B1 GPCRs are critical regulators of homeostasis, and as such, many are important drug targets. While most transmembrane proteins, including GPCRs, are recalcitrant to crystallization, recent advances in electron cryo-microscopy (cryo-EM) have facilitated a rapid expansion of the structural understanding of membrane proteins. As a testament to this success, structures for all the class B1 receptors bound to G proteins have been determined by cryo-EM in the past five years. Further advances in cryo-EM have uncovered dynamics of these receptors, ligands, and signalling partners. Here, we examine the recent structural underpinnings of the class B1 GPCRs with an emphasis on structure-function relationships.

Molecular recognition of morphine and fentanyl by the human μ-opioid receptor

Morphine and fentanyl are among the most used opioid drugs that confer analgesia and unwanted side effects through both G protein and arrestin signaling pathways of μ-opioid receptor (μOR). Here, we report structures of the human μOR-G protein complexes bound to morphine and fentanyl, which uncover key differences in how they bind the receptor. We also report structures of μOR bound to TRV130, PZM21, and SR17018, which reveal preferential interactions of these agonists with TM3 side of the ligand-binding pocket rather than TM6/7 side. In contrast, morphine and fentanyl form dual interactions with both TM3 and TM6/7 regions. Mutations at the TM6/7 interface abolish arrestin recruitment of μOR promoted by morphine and fentanyl. Ligands designed to reduce TM6/7 interactions display preferential G protein signaling. Our results provide crucial insights into fentanyl recognition and signaling of μOR, which may facilitate rational design of next-generation analgesics.

Filling of a water-free void explains the allosteric regulation of the β1-adrenergic receptor by cholesterol

Recent high-pressure NMR results indicate that the preactive conformation of the β₁-adrenergic receptor (β₁AR) harbours completely empty cavities of ~100 ų volume, which disappear in the active conformation of the receptor. Here we have localized these cavities using X-ray crystallography of xenon-derivatized β₁AR crystals. One of the cavities is in direct contact with the cholesterol-binding pocket. Solution NMR shows that addition of the cholesterol analogue cholesteryl hemisuccinate impedes the formation of the active conformation of detergent-solubilized β₁AR by blocking conserved G protein-coupled receptor microswitches, concomitant with an affinity reduction of both isoprenaline and G protein-mimicking nanobody Nb80 for β₁AR detected by isothermal titration calorimetry. This wedge-like action explains the function of cholesterol as a negative allosteric modulator of β₁AR. A detailed understanding of G protein-coupled receptor regulation by cholesterol by filling of a dry void and the easy scouting for such voids by xenon may provide new routes for the development of allosteric drugs.

Computational and experimental approaches to probe GPCR activation and signaling

G protein-coupled receptors (GPCRs) regulate different physiological functions, e.g., sensation, growth, digestion, reproductivity, nervous and immune systems response, and many others. In eukaryotes, they are also responsible for intercellular communication in response to pathogens. The major primary messengers binding to these cell-surface receptors constitute small-molecule or peptide hormones and neurotransmitters, nucleotides, lipids as well as small proteins. The simplicity of the way how GPCR signaling can be regulated by their endogenous agonists prompted the usage of GPCRs as major drug targets in modern pharmacology. Drugs targeting GPCRs inhibit pathological processes at the very beginning. This enables to significantly reduce the occurrence of morphological changes caused by diseases. Until recently, X-ray crystallography was the method of the first choice to obtain high-resolution structural information about GPCRs. Following X-ray crystallography, cryo-EM gained attention in GPCR studies as a quick and low-cost alternative. FRET microscopy is also widely used for GPCRs in the analysis of protein-protein interactions (PPIs) in intact cells as well as for screening purposes. Regarding computational methods, molecular dynamics (MD) for many years has proven its usefulness in studying the GPCR activation. MODELLER and Rosetta were widely used to generate preliminary homology models of GPCRs for MD simulation systems. Apart from the conventional all-atom approach with explicitly defined solvent, also other techniques have been applied to GPCRs, e.g., MARTINI or hybrid methods involving the coarse-grained representation, less demanding regarding computational resources, and thus offering much larger simulation timescales.

Yeast-based directed-evolution for high-throughput structural stabilization of G protein-coupled receptors (GPCRs)

The immense potential of G protein-coupled receptors (GPCRs) as targets for drug discovery is not fully realized due to the enormous difficulties associated with structure elucidation of these profoundly unstable membrane proteins. The existing methods of GPCR stability-engineering are cumbersome and low-throughput; in addition, the scope of GPCRs that could benefit from these techniques is limited. Here, we present a yeast-based screening platform for a single-step isolation of GRCR variants stable in the presence of short-chain detergents, a feature essential for their successful crystallization using vapor diffusion method. The yeast detergent-resistant cell wall presents a unique opportunity for compartmentalization, to physically link the receptor's phenotype to its encoding DNA, and thus enable discovery of stable GPCR variants with unprecedent efficiency. The scope of mutations identified by the method reveals a surprising amenability of the GPCR scaffold to stabilization, and suggests an intriguing possibility of amending the stability properties of GPCR by varying the structural status of the C-terminus.

Angiotensin and Endothelin Receptor Structures With Implications for Signaling Regulation and Pharmacological Targeting

In conjunction with the endothelin (ET) type A (ETAR) and type B (ETBR) receptors, angiotensin (AT) type 1 (AT1R) and type 2 (AT2R) receptors, are peptide-binding class A G-protein-coupled receptors (GPCRs) acting in a physiologically overlapping context. Angiotensin receptors (ATRs) are involved in regulating cell proliferation, as well as cardiovascular, renal, neurological, and endothelial functions. They are important therapeutic targets for several diseases or pathological conditions, such as hypertrophy, vascular inflammation, atherosclerosis, angiogenesis, and cancer. Endothelin receptors (ETRs) are expressed primarily in blood vessels, but also in the central nervous system or epithelial cells. They regulate blood pressure and cardiovascular homeostasis. Pathogenic conditions associated with ETR dysfunctions include cancer and pulmonary hypertension. While both receptor groups are activated by their respective peptide agonists, pathogenic autoantibodies (auto-Abs) can also activate the AT1R and ETAR accompanied by respective clinical conditions. To date, the exact mechanisms and differences in binding and receptor-activation mediated by auto-Abs as opposed to endogenous ligands are not well understood. Further, several questions regarding signaling regulation in these receptors remain open. In the last decade, several receptor structures in the apo- and ligand-bound states were determined with protein X-ray crystallography using conventional synchrotrons or X-ray Free-Electron Lasers (XFEL). These inactive and active complexes provide detailed information on ligand binding, signal induction or inhibition, as well as signal transduction, which is fundamental for understanding properties of different activity states. They are also supportive in the development of pharmacological strategies against dysfunctions at the receptors or in the associated signaling axis. Here, we summarize current structural information for the AT1R, AT2R, and ETBR to provide an improved molecular understanding.

Sequence coevolution and structure stabilization modulate olfactory receptor expression

Olfactory receptors (ORs) belong to class A G-protein coupled receptors (GPCRs) and are activated by a variety of odorants. To date, there is no three-dimensional structure of an OR available. One of the major bottlenecks in obtaining purified protein for structural studies of ORs is their poor expression in heterologous cells. To design mutants that enhance expression and thereby enable protein purification, we first identified computable physical properties that recapitulate OR and class A GPCR expression and further conducted an iterative computational prediction-experimental test cycle and generated human OR mutants that express as high as biogenic amine receptors for which structures have been solved. In the process of developing the computational method to recapitulate the expression of ORs in membranes, we identified properties, such as amino acid sequence coevolution, and the strength of the interactions between intracellular loop 1 (ICL1) and the helix 8 region of ORs, to enhance their heterologous expression. We identified mutations that are directly located in these regions as well as other mutations not located in these regions but allosterically strengthen the ICL1-helix 8 enhance expression. These mutants also showed functional responses to known odorants. This method to enhance heterologous expression of mammalian ORs will facilitate high-throughput "deorphanization" of ORs, and enable OR purification for biochemical and structural studies to understand odorant-OR interactions.

Magic angle spinning NMR of G protein-coupled receptors

G protein-coupled receptors (GPCRs) have a simple seven transmembrane helix architecture which has evolved to recognize a diverse number of chemical signals. The more than 800 GPCRs encoded in the human genome function as receptors for vision, smell and taste, and mediate key physiological processes. Consequently, these receptors are a major target for pharmaceuticals. Protein crystallography and electron cryo-microscopy have provided high resolution structures of many GPCRs in both active and inactive conformations. However, these structures have not sparked a surge in rational drug design, in part because GPCRs are inherently dynamic and the structural changes induced by ligand or drug binding to stabilize inactive or active conformations are often subtle rearrangements in packing or hydrogen-bonding interactions. NMR spectroscopy provides a sensitive probe of local structure and dynamics at specific sites within these receptors as well as global changes in receptor structure and dynamics. These methods can also capture intermediate states and conformations with low populations that provide insights into the activation pathways. We review the use of solid-state magic angle spinning NMR to address the structure and activation mechanisms of GPCRs. The focus is on the large and diverse class A family of receptors. We highlight three specific class A GPCRs in order to illustrate how solid-state, as well as solution-state, NMR spectroscopy can answer questions in the field involving how different GPCR classes and subfamilies are activated by their associated ligands, and how small molecule drugs can modulate GPCR activation.

Why Are We Still Cloning Melatonin Receptors? A Commentary

Cloning may seem to be a view from the past. The time before software, computers and AI were invented. It seems to us worth discussing these points in view of our favorite target: the melatoninergic system. In a few stances, it might be important to point out that even in the new era of dry science, there is still a need to experiment and to prove at the bench that our in silico assertions are right. Most of the living animals express to some extend the melatonin receptors. Some of these animal genomes were completely or partially sequenced, and it is tempting to extract from this huge information the sequence(s) of our favorite genes (MLT receptors). Then, why bother cloning, as opposed to simply built the gene and express it in a host cell? Because the genetic boundaries of the expressed sequence(s) are not 100% sure. Because the melatonin receptor gene(s) comprise a first exon 25,000 base pair far from the second one and the limits between this Ex1 and In1-as between In1 and Ex2-are subject to changes that might have a huge impact on the biochemical properties of the receptor, once expressed. Because a receptor is a biochemical entity with characteristics that are important for the functioning of this particular pathway, and more generally, for the functioning of life.

The Novel Application of Geometric Morphometrics with Principal Component Analysis to Existing G Protein-Coupled Receptor (GPCR) Structures

The G protein-coupled receptor (GPCR) superfamily is a large group of membrane proteins which, because of their vast involvement in cell signalling pathways, are implicated in a plethora of disease states and are therefore considered to be key drug targets. Despite advances in techniques to study these receptors, current prophylaxis is often limited due to the challenging nature of their dynamic, complex structures. Greater knowledge and understanding of their intricate structural rearrangements will therefore undoubtedly aid structure-based drug design against GPCRs. Disciplines such as anthropology and palaeontology often use geometric morphometrics to measure variation between shapes and we have therefore applied this technique to analyse GPCR structures in a three-dimensional manner, using principal component analysis. Our aim was to create a novel system able to discriminate between GPCR structures and discover variation between them, correlated with a variety of receptor characteristics. This was conducted by assessing shape changes at the extra- and intracellular faces of the transmembrane helix bundle, analysing the XYZ coordinates of the amino acids at those positions. We have demonstrated that GPCR structures can be classified based on characteristics such as activation state, bound ligands and fusion proteins, with the most significant results focussed at the intracellular face. Conversely, our analyses provide evidence that thermostabilising mutations do not cause significant differences when compared to non-mutated GPCRs. We believe that this is the first time geometric morphometrics has been applied to membrane proteins on this scale, and believe it can be used as a future tool in sense-checking newly resolved structures and planning experimental design.

Evidence that specific interactions play a role in the cholesterol sensitivity of G protein-coupled receptors

G protein-coupled receptors (GPCRs) are known to be modulated by membrane cholesterol levels, but whether or not the effects are caused by specific receptor-cholesterol interactions or cholesterol's general effects on the membrane is not well-understood. We performed coarse-grained molecular dynamics (CGMD) simulations coupled with structural bioinformatics approaches on the β2-adrenergic receptor (β2AR) and the cholecystokinin (CCK) receptor subfamily. The β2AR has been shown to be sensitive to membrane cholesterol and cholesterol molecules have been clearly resolved in numerous β2AR crystal structures. The two CCK receptors are highly homologous and preserve similar cholesterol recognition motifs but despite their homology, CCK1R shows functional sensitivity to membrane cholesterol while CCK2R does not. Our results offer new insights into how cholesterol modulates GPCR function by showing cholesterol interactions with β2AR that agree with previously published data; additionally, we observe differential and specific cholesterol binding in the CCK receptor subfamily while revealing a previously unreported Cholesterol Recognition Amino-acid Consensus (CRAC) sequence that is also conserved across 38% of class A GPCRs. A thermal denaturation assay (LCP-Tm) shows that mutation of a conserved CRAC sequence on TM7 of the β2AR affects cholesterol stabilization of the receptor in a lipid bilayer. The results of this study provide a better understanding of receptor-cholesterol interactions that can contribute to novel and improved therapeutics for a variety of diseases.

Clinical, Biochemical, and Genetic Heterogeneity in Glutaric Aciduria Type II Patients

The variants of electron transfer flavoprotein (ETFA, ETFB) and ETF dehydrogenase (ETFDH) are the leading cause of glutaric aciduria type II (GA-II). In this study, we identified 13 patients harboring six variants of two genes associated with GA-II. Out of the six variants, four were missense, and two were frameshift mutations. A missense variant (ETFDH:p.Gln269His) was observed in a homozygous state in nine patients. Among nine patients, three had experienced metabolic crises with recurrent vomiting, abdominal pain, and nausea. In one patient with persistent metabolic acidosis, hypoglycemia, and a high anion gap, the ETFDH:p.Gly472Arg, and ETFB:p.Pro94Thrfs*8 variants were identified in a homozygous, and heterozygous state, respectively. A missense variant ETFDH:p.Ser442Leu was detected in a homozygous state in one patient with metabolic acidosis, hypoglycemia, hyperammonemia and liver dysfunction. The ETFDH:p.Arg41Leu, and ETFB:p.Ile346Phefs*19 variants were observed in a homozygous state in one patient each. Both these variants have not been reported so far. In silico approaches were used to evaluate the pathogenicity and structural changes linked with these six variants. Overall, the results indicate the importance of a newborn screening program and genetic investigations for patients with GA-II. Moreover, careful interpretation and correlation of variants of uncertain significance with clinical and biochemical findings are needed to confirm the pathogenicity of such variants.

Optimizing the Expression of Human Dopamine Receptors in Escherichia coli

The human dopamine receptors D_2S and D₃ belong to the group of G protein-coupled receptors (GPCRs) and are important drug targets. Structural analyses and development of new receptor subtype specific drugs have been impeded by low expression yields or receptor instability. Fusing the T4 lysozyme into the intracellular loop 3 improves crystallization but complicates conformational studies. To circumvent these problems, we expressed the human D_2S and D₃ receptors in Escherichia coli using different N- and C-terminal fusion proteins and thermostabilizing mutations. We optimized expression times and used radioligand binding assays with whole cells and membrane homogenates to evaluate K_D-values and the number of receptors in the cell membrane. We show that the presence but not the type of a C-terminal fusion protein is important. Bacteria expressing receptors capable of ligand binding can be selected using FACS analysis and a fluorescently labeled ligand. Improved receptor variants can thus be generated using error-prone PCR. Subsequent analysis of clones showed the distribution of mutations over the whole gene. Repeated cycles of PCR and FACS can be applied for selecting highly expressing receptor variants with high affinity ligand binding, which in the future can be used for analytical studies.

Membrane Protein Stabilization Strategies for Structural and Functional Studies

Accounting for nearly two-thirds of known druggable targets, membrane proteins are highly relevant for cell physiology and pharmacology. In this regard, the structural determination of pharmacologically relevant targets would facilitate the intelligent design of new drugs. The structural biology of membrane proteins is a field experiencing significant growth as a result of the development of new strategies for structure determination. However, membrane protein preparation for structural studies continues to be a limiting step in many cases due to the inherent instability of these molecules in non-native membrane environments. This review describes the approaches that have been developed to improve membrane protein stability. Membrane protein mutagenesis, detergent selection, lipid membrane mimics, antibodies, and ligands are described in this review as approaches to facilitate the production of purified and stable membrane proteins of interest for structural and functional studies.

GPCR Activation States Induced by Nanobodies and Mini-G Proteins Compared by NMR Spectroscopy

In this work, we examine methyl nuclear magnetic resonance (NMR) spectra of the methionine ε-[¹³CH₃] labelled thermostabilized β₁ adrenergic receptor from turkey in association with a variety of different effectors, including mini-G_s and nanobody 60 (Nb60), which have not been previously studied in complex with β₁ adrenergic receptor (β₁AR) by NMR. Complexes with pindolol and Nb60 induce highly similar inactive states of the receptor, closely resembling the resting state conformational ensemble. We show that, upon binding of mini-Gs or nanobody 80 (Nb80), large allosteric changes throughout the receptor take place. The conformation of tβ₁AR stabilized by the native-like mini-G_s protein is highly similar to the conformation induced by the currently used surrogate Nb80. Interestingly, in both cases residual dynamics are present, which were not observed in the resting states. Finally, we reproduce a pharmaceutically relevant situation, where an antagonist abolishes the interaction of the receptor with the mini-G protein in a competitive manner, validating the functional integrity of our preparation. The presented system is therefore well suited for reproducing the individual steps of the activation cycle of a G protein-coupled receptor (GPCR) in vitro and serves as a basis for functional and pharmacological characterizations of more native-like systems in the future.

An improved fluorescent tag and its nanobodies for membrane protein expression, stability assay, and purification

Green fluorescent proteins (GFPs) are widely used to monitor membrane protein expression, purification, and stability. An ideal reporter should be stable itself and provide high sensitivity and yield. Here, we demonstrate that a coral (Galaxea fascicularis) thermostable GFP (TGP) is by such reasons an improved tag compared to the conventional jellyfish GFPs. TGP faithfully reports membrane protein stability at temperatures near 90 °C (20-min heating). By contrast, the limit for the two popular GFPs is 64 °C and 74 °C. Replacing GFPs with TGP increases yield for all four test membrane proteins in four expression systems. To establish TGP as an affinity tag for membrane protein purification, several high-affinity synthetic nanobodies (sybodies), including a non-competing pair, are generated, and the crystal structure of one complex is solved. Given these advantages, we anticipate that TGP becomes a widely used tool for membrane protein structural studies.

Ligand-Specific Factors Influencing GLP-1 Receptor Post-Endocytic Trafficking and Degradation in Pancreatic Beta Cells

The glucagon-like peptide-1 receptor (GLP-1R) is an important regulator of blood glucose homeostasis. Ligand-specific differences in membrane trafficking of the GLP-1R influence its signalling properties and therapeutic potential in type 2 diabetes. Here, we have evaluated how different factors combine to control the post-endocytic trafficking of GLP-1R to recycling versus degradative pathways. Experiments were performed in primary islet cells, INS-1 832/3 clonal beta cells and HEK293 cells, using biorthogonal labelling of GLP-1R to determine its localisation and degradation after treatment with GLP-1, exendin-4 and several further GLP-1R agonist peptides. We also characterised the effect of a rare GLP1R coding variant, T149M, and the role of endosomal peptidase endothelin-converting enzyme-1 (ECE-1), in GLP1R trafficking. Our data reveal how treatment with GLP-1 versus exendin-4 is associated with preferential GLP-1R targeting towards a recycling pathway. GLP-1, but not exendin-4, is a substrate for ECE-1, and the resultant propensity to intra-endosomal degradation, in conjunction with differences in binding affinity, contributes to alterations in GLP-1R trafficking behaviours and degradation. The T149M GLP-1R variant shows reduced signalling and internalisation responses, which is likely to be due to disruption of the cytoplasmic region that couples to intracellular effectors. These observations provide insights into how ligand- and genotype-specific factors can influence GLP-1R trafficking.

Thermostability of a recombinant G protein-coupled receptor expressed at high level in mammalian cell culture

Rational design of pharmaceutical drugs targeting integral membrane G protein-coupled receptors (GPCR) requires thorough understanding of ligand binding and mechanism of activation through high resolution structural studies of purified proteins. Due to inherent conformational flexibility of GPCR, stabilization of these proteins solubilized from cell membranes into detergents is a challenging task. Here, we take advantage of naturally occurring post-translational modifications for stabilization of purified GPCR in detergent micelles. The recombinant cannabinoid CB₂ receptor was expressed at high yield in Expi293F mammalian cell cultures, solubilized and purified in Façade detergent. We report superior stability of the mammalian cell-expressed receptor compared to its E.coli-expressed counterpart, due to contributions from glycosylation of the N terminus and palmitoylation of the C terminus of CB₂. Finally, we demonstrate that the mammalian Expi293F amino acid labelling kit is suitable for preparation of multi-milligram quantities of high quality, selectively stable isotope-labeled GPCR for studies by nuclear magnetic resonance.

MPTherm-pred: Analysis and Prediction of Thermal Stability Changes upon Mutations in Transmembrane Proteins

The stability of membrane proteins differs from globular proteins due to the presence of nonpolar membrane-spanning regions. Using a dataset of 929 membrane protein mutations whose effects on thermal stability (ΔT_m) were experimentally determined, we found that the average ΔT_m due to 190 stabilizing and 232 destabilizing mutations occurring in membrane-spanning regions are 2.43(3.1) °C and -5.48(5.5) °C, respectively. The ΔT_m values for mutations occurring in solvent-exposed regions are 2.56(2.82) and - 6.8(7.2) °C. We have systematically analyzed the factors influencing the stability of mutants and observed that changes in hydrophobicity, number of contacts between Cα atoms and frequency of aliphatic residues are important determinants of the stability change induced by mutations occurring in membrane-spanning regions. We have developed structure- and sequence-based machine learning predictors of ΔT_m due to mutations specifically for membrane proteins. They showed a correlation and mean absolute error (MAE) of 0.72 and 2.85 °C, respectively, between experimental and predicted ΔT_m for mutations in membrane-spanning regions on 10-fold group-wise cross-validation. The average correlation and MAE for mutations in aqueous regions are 0.73 and 3.7 °C, respectively. These MAE values are about 50% lower than standard deviations from the mean ΔT_m values. The reliability of the method was affirmed on a test set of mutations occurring in evolutionary independent protein sequences. The developed MPTherm-pred server for predicting thermal stability changes upon mutations in membrane proteins is available at https://web.iitm.ac.in/bioinfo2/mpthermpred/. Our results provide insights into factors influencing the stability of membrane proteins and can aid in designing mutants that are more resistant to thermal stress.

Benchmarking GPCR homology model template selection in combination with de novo loop generation

G protein-coupled receptors (GPCR) comprise the largest family of membrane proteins and are of considerable interest as targets for drug development. However, many GPCR structures remain unsolved. To address the structural ambiguity of these receptors, computational tools such as homology modeling and loop modeling are often employed to generate predictive receptor structures. Here we combined both methods to benchmark a protocol incorporating homology modeling based on a locally selected template and extracellular loop modeling that additionally evaluates the presence of template ligands during these modeling steps. Ligands were also docked using three docking methods and two pose selection methods to elucidate an optimal ligand pose selection method. Results suggest that local template-based homology models followed by loop modeling produce more accurate and predictive receptor models than models produced without loop modeling, with decreases in average receptor and ligand RMSD of 0.54 Å and 2.91 Å, respectively. Ligand docking results showcased the ability of MOE induced fit docking to produce ligand poses with atom root-mean-square deviation (RMSD) values at least 0.20 Å lower (on average) than the other two methods benchmarked in this study. In addition, pose selection methods (software-based scoring, ligand complementation) selected lower RMSD poses with MOE induced fit docking than either of the other methods (averaging at least 1.57 Å lower), indicating that MOE induced fit docking is most suited for docking into GPCR homology models in our hands. In addition, target receptor models produced with a template ligand present throughout the modeling process most often produced target ligand poses with RMSD values ≤ 4.5 Å and Tanimoto coefficients > 0.6 after selection based on ligand complementation than target receptor models produced in the absence of template ligands. Overall, the findings produced by this study support the use of local template homology modeling in combination with de novo ECL2 modeling in the presence of a ligand from the template crystal structure to generate GPCR models intended to study ligand binding interactions.

Docking of peptides to GPCRs using a combination of CABS-dock with FlexPepDock refinement

The structural description of peptide ligands bound to G protein-coupled receptors (GPCRs) is important for the discovery of new drugs and deeper understanding of the molecular mechanisms of life. Here we describe a three-stage protocol for the molecular docking of peptides to GPCRs using a set of different programs: (1) CABS-dock for docking fully flexible peptides; (2) PD2 method for the reconstruction of atomistic structures from C-alpha traces provided by CABS-dock and (3) Rosetta FlexPepDock for the refinement of protein–peptide complex structures and model scoring. We evaluated the proposed protocol on the set of seven different GPCR–peptide complexes (including one containing a cyclic peptide), for which crystallographic structures are available. We show that CABS-dock produces high resolution models in the sets of top-scored models. These sets of models, after reconstruction to all-atom representation, can be further improved by Rosetta high-resolution refinement and/or minimization, leading in most of the cases to sub-Angstrom accuracy in terms of interface root-mean-square-deviation measure.

Computational Investigations on the Binding Mode of Ligands for the Cannabinoid-Activated G Protein-Coupled Receptor GPR18

GPR18 is an orphan G protein-coupled receptor (GPCR) expressed in cells of the immune system. It is activated by the cannabinoid receptor (CB) agonist ∆⁹-tetrahydrocannabinol (THC). Several further lipids have been proposed to act as GPR18 agonists, but these results still require unambiguous confirmation. In the present study, we constructed a homology model of the human GPR18 based on an ensemble of three GPCR crystal structures to investigate the binding modes of the agonist THC and the recently reported antagonists which feature an imidazothiazinone core to which a (substituted) phenyl ring is connected via a lipophilic linker. Docking and molecular dynamics simulation studies were performed. As a result, a hydrophobic binding pocket is predicted to accommodate the imidazothiazinone core, while the terminal phenyl ring projects towards an aromatic pocket. Hydrophobic interaction of Cys251 with substituents on the phenyl ring could explain the high potency of the most potent derivatives. Molecular dynamics simulation studies suggest that the binding of imidazothiazinone antagonists stabilizes transmembrane regions TM1, TM6 and TM7 of the receptor through a salt bridge between Asp118 and Lys133. The agonist THC is presumed to bind differently to GPR18 than to the distantly related CB receptors. This study provides insights into the binding mode of GPR18 agonists and antagonists which will facilitate future drug design for this promising potential drug target.

Characterizing Interhelical Interactions of G-Protein Coupled Receptors with the Fragment Molecular Orbital Method

G-protein coupled receptors (GPCRs) are the largest superfamily of membrane proteins, regulating almost every aspect of cellular activity and serving as key targets for drug discovery. We have identified an accurate and reliable computational method to characterize the strength and chemical nature of the interhelical interactions between the residues of transmembrane (TM) domains during different receptor activation states, something that cannot be characterized solely by visual inspection of structural information. Using the fragment molecular orbital (FMO) quantum mechanics method to analyze 35 crystal structures representing different branches of the class A GPCR family, we have identified 69 topologically equivalent TM residues that form a consensus network of 51 inter-TM interactions, providing novel results that are consistent with and help to rationalize experimental data. This discovery establishes a comprehensive picture of how defined molecular forces govern specific interhelical interactions which, in turn, support the structural stability, ligand binding, and activation of GPCRs.

Thermostabilization of Membrane Proteins by Consensus Mutation: A Case Study for a Fungal Δ8-7 Sterol Isomerase

Membrane proteins are generally challenging to work with because of their notorious instability. Protein engineering has been used increasingly to thermostabilize labile membrane proteins such as G-protein-coupled receptors for structural and functional studies in recent years. Two major strategies exist. Scanning mutagenesis systematically eliminates destabilizing residues, whereas the consensus approach assembles mutants with the most frequent residues among selected homologs, bridging sequence conservation with stability. Here, we applied the consensus concept to stabilize a fungal homolog of the human sterol Δ8-7 isomerase, a 26.4 kDa protein with five transmembrane helices. The isomerase is also called emopamil-binding protein (EBP), as it binds this anti-ischemic drug with high affinity. The wild-type had an apparent melting temperature (T_m) of 35.9 °C as measured by the fluorescence-detection size-exclusion chromatography-based thermostability assay. A total of 87 consensus mutations sourced from 22 homologs gained expression level and thermostability, increasing the apparent T_m to 69.9 °C at the cost of partial function loss. Assessing the stability and activity of several systematic chimeric constructs identified a construct with an apparent T_m of 79.8 °C and two regions for function rescue. Further back-mutations of the chimeric construct in the two target regions yielded the final construct with similar apparent activity to the wild-type and an elevated T_m of 88.8 °C, totaling an increase of 52.9 °C. The consensus approach is effective and efficient because it involves fewer constructs compared with scanning mutagenesis. Our results should encourage more use of the consensus strategy for membrane protein thermostabilization.

Structural insights into melatonin receptors

The long-anticipated high-resolution structures of the human melatonin G protein-coupled receptors MT₁ and MT₂ , involved in establishing and maintaining circadian rhythm, were obtained in complex with two melatonin analogs and two approved anti-insomnia and antidepression drugs using X-ray free-electron laser serial femtosecond crystallography. The structures shed light on the overall conformation and unusual structural features of melatonin receptors, as well as their ligand binding sites and the melatonergic pharmacophore, thereby providing insights into receptor subtype selectivity. The structures revealed an occluded orthosteric ligand binding site with a membrane-buried channel for ligand entry in both receptors, and an additional putative ligand entry path in MT₂ from the extracellular side. This unexpected ligand entry mode contributes to facilitating the high specificity with which melatonin receptors bind their cognate ligand and exclude structurally similar molecules such as serotonin, the biosynthetic precursor of melatonin. Finally, the MT₂ structure allowed accurate mapping of type 2 diabetes-related single-nucleotide polymorphisms, where a clustering of residues in helices I and II on the protein-membrane interface was observed which could potentially influence receptor oligomerization. The role of receptor oligomerization is further discussed in light of the differential interaction of MT₁ and MT₂ with GPR50, a regulatory melatonin coreceptor. The melatonin receptor structures will facilitate design of selective tool compounds to further dissect the specific physiological function of each receptor subtype as well as provide a structural basis for next-generation sleeping aids and other drugs targeting these receptors with higher specificity and fewer side effects.

NMR investigation of protein-ligand interactions for G-protein coupled receptors

In this review, we report NMR studies of ligand-GPCR interactions, including both ligand-observed and protein-observed NMR experiments. Published studies exemplify how NMR can be used as a powerful tool to design novel GPCR ligands and investigate the ligand-induced conformational changes of GPCRs. The strength of NMR also lies in its capability to explore the diverse signaling pathways and probe the allosteric modulation of these highly dynamic receptors. By offering unique opportunities for the identification, structural and functional characterization of GPCR ligands, NMR will likely play a major role for the generation of novel molecules both as new tools for the understanding of the GPCR function and as therapeutic compounds for a large diversity of pathologies.

Computational design for thermostabilization of GPCRs

GPCR superfamily is the largest clinically relevant family of targets in human genome; however, low thermostability and high conformational plasticity of these integral membrane proteins make them notoriously hard to handle in biochemical, biophysical, and structural experiments. Here, we describe the recent advances in computational approaches to design stabilizing mutations for GPCR that take advantage of the structural and sequence conservation properties of the receptors, and employ machine learning on accumulated mutation data for the superfamily. The fast and effective computational tools can provide a viable alternative to existing experimental mutation screening and are poised for further improvements with expansion of thermostability datasets for training the machine learning models. The rapidly growing practical applications of computational stability design streamline GPCR structure determination and may contribute to more efficient drug discovery.

Computational design of thermostabilizing point mutations for G protein-coupled receptors

Engineering of GPCR constructs with improved thermostability is a key for successful structural and biochemical studies of this transmembrane protein family, targeted by 40% of all therapeutic drugs. Here we introduce a comprehensive computational approach to effective prediction of stabilizing mutations in GPCRs, named CompoMug, which employs sequence-based analysis, structural information, and a derived machine learning predictor. Tested experimentally on the serotonin 5-HT_2C receptor target, CompoMug predictions resulted in 10 new stabilizing mutations, with an apparent thermostability gain ~8.8°C for the best single mutation and ~13°C for a triple mutant. Binding of antagonists confers further stabilization for the triple mutant receptor, with total gains of ~21°C as compared to wild type apo 5-HT_2C. The predicted mutations enabled crystallization and structure determination for the 5-HT_2C receptor complexes in inactive and active-like states. While CompoMug already shows high 25% hit rate and utility in GPCR structural studies, further improvements are expected with accumulation of structural and mutation data.

The trillions of cells in the human body rely on receptors that sit in their cell membranes to communicate with each other. Hundreds of different receptors belong to the G protein-coupled receptor superfamily (called GPCRs for short) and play vital roles in the all organs and bodily systems. Indeed, GPCRs are the targets for almost 40% of therapeutic drugs. As such, deciphering the shape and activity of GPCRs is key to understanding the normal workings of the human biology and could help scientists discover new treatments for various diseases, from depression to high blood pressure to cancer. These receptors, however, are notoriously flimsy and unstable, making them difficult to work with in the laboratory.

Different approaches have been developed to make GPCRs more stable, usually by swapping one or a few of the amino acid building blocks in the protein for other amino acids. Currently, this requires a costly and slow trial-and-error approach in which each amino acid out of 300-400 in the protein is mutated and tested experimentally.

To speed up and reduce the cost of the process, Popov et al. asked if a computer could predict which mutations in the protein would stabilize it, meaning that fewer proteins would actually need to be tested. Four computer algorithms based on four different principles were developed and verified. The first one compares the target GPCR to other closely related receptors, trying to detect variations that cause the instability. The second tries to build in specific stabilizing interactions, or “bridges”, between different parts of the receptor. The third algorithm searches the known structures of other GPCRs for useful mutations. Finally, the fourth one uses accumulated data on the stability of hundreds of mutations in different GPCRs to train a machine learning predictor to recognize stabilizing mutations.

All four algorithms produced useful predictions in a real-life project. Indeed, when combined in one computational tool, named CompoMug, the algorithms made it possible to detect optimal mutations in a human GPCR called 5-HT_2C. This made the protein much easier to work with in the laboratory, and ultimately helped to solve its three-dimensional structure (which was reported in a separate study, published earlier in 2018)

The 5-HT_2C receptor is involved in regulating, among other things, mood and appetite. Details of its structure might therefore help researchers to design new antidepressants and obesity treatments. Moreover, CompoMug is already helping structural biologists to solve the structures of other GPCRs, which will further facilitate many aspects of GPCR drug discovery.

Enabling STD-NMR fragment screening using stabilized native GPCR: A case study of adenosine receptor

Structural studies of integral membrane proteins have been limited by the intrinsic conformational flexibility and the need to stabilize the proteins in solution. Stabilization by mutagenesis was very successful for structural biology of G protein-coupled receptors (GPCRs). However, it requires heavy protein engineering and may introduce structural deviations. Here we describe the use of specific calixarenes-based detergents for native GPCR stabilization. Wild type, full length human adenosine A_2A receptor was used to exemplify the approach. We could stabilize native, glycosylated, non-aggregated and homogenous A_2AR that maintained its ligand binding capacity. The benefit of the preparation for fragment screening, using the Saturation-Transfer Difference nuclear magnetic resonance (STD-NMR) experiment is reported. The binding of the agonist adenosine and the antagonist caffeine were observed and competition experiments with CGS-21680 and ZM241385 were performed, demonstrating the feasibility of the STD-based fragment screening on the native A_2A receptor. Interestingly, adenosine was shown to bind a second binding site in the presence of the agonist CGS-21680 which corroborates published results obtained with molecular dynamics simulation. Fragment-like compounds identified using STD-NMR showed antagonistic effects on A_2AR in the cAMP cellular assay. Taken together, our study shows that stabilization of native GPCRs represents an attractive approach for STD-based fragment screening and drug design.

Serial Femtosecond Crystallography of G Protein-Coupled Receptors

G protein-coupled receptors (GPCRs) represent a large superfamily of membrane proteins that mediate cell signaling and regulate a variety of physiological processes in the human body. Structure-function studies of this superfamily were enabled a decade ago by multiple breakthroughs in technology that included receptor stabilization, crystallization in a membrane environment, and microcrystallography. The recent emergence of X-ray free-electron lasers (XFELs) has further accelerated structural studies of GPCRs and other challenging proteins by overcoming radiation damage and providing access to high-resolution structures and dynamics using micrometer-sized crystals. Here, we summarize key technology advancements and major milestones of GPCR research using XFELs and provide a brief outlook on future developments in the field.

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.

Structural stability of purified human CFTR is systematically improved by mutations in nucleotide binding domain 1

The Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) is an ABC transporter containing two transmembrane domains forming a chloride ion channel, and two nucleotide binding domains (NBD1 and NBD2). CFTR has presented a formidable challenge to obtain monodisperse, biophysically stable protein. Here we report a comprehensive study comparing effects of single and multiple NBD1 mutations on stability of both the NBD1 domain alone and on purified full length human CFTR. Single mutations S492P, A534P, I539T acted additively, and when combined with M470V, S495P, and R555K cumulatively yielded an NBD1 with highly improved structural stability. Strategic combinations of these mutations strongly stabilized the domain to attain a calorimetric T_m > 70 °C. Replica exchange molecular dynamics simulations on the most stable 6SS-NBD1 variant implicated fluctuations, electrostatic interactions and side chain packing as potential contributors to improved stability. Progressive stabilization of NBD1 directly correlated with enhanced structural stability of full-length CFTR protein. Thermal unfolding of the stabilized CFTR mutants, monitored by changes in intrinsic fluorescence, demonstrated that Tm could be shifted as high as 67.4 °C in 6SS-CFTR, more than 20 °C higher than wild-type. H1402S, an NBD2 mutation, conferred CFTR with additional thermal stability, possibly by stabilizing an NBD-dimerized conformation. CFTR variants with NBD1-stabilizing mutations were expressed at the cell surface in mammalian cells, exhibited ATPase and channel activity, and retained these functions to higher temperatures. The capability to produce enzymatically active CFTR with improved structural stability amenable to biophysical and structural studies will advance mechanistic investigations and future cystic fibrosis drug development.

Synergistic Use of GPCR Modeling and SDM Experiments to Understand Ligand Binding

There is a substantial amount of historical ligand binding data available from site-directed mutagenesis (SDM) studies of many different GPCR subtypes. This information was generated prior to the wave of GPCR crystal structure, in an effort to understand ligand binding with a view to drug discovery. Concerted efforts to determine the atomic structure of GPCRs have proven extremely successful and there are now more than 80 GPCR crystal structure in the PDB database, many of which have been obtained in the presence of receptor ligands and associated G proteins. These structural data enable the generation of computational model structures for all GPCRs, including those for which crystal structures do not yet exist. The power of these models in designing novel ligands, especially those with improved residence times, and for better understanding receptor function can be enhanced tremendously by combining them synergistically with historic SDM ligand binding data. Here, we describe a protocol by which historic SDM binding data and receptor models may be used together to identify novel key residues for mutagenesis studies.

High-throughput mutagenesis using a two-fragment PCR approach

Site-directed scanning mutagenesis is a powerful protein engineering technique which allows studies of protein functionality at single amino acid resolution and design of stabilized proteins for structural and biophysical work. However, creating libraries of hundreds of mutants remains a challenging, expensive and time-consuming process. The efficiency of the mutagenesis step is the key for fast and economical generation of such libraries. PCR artefacts such as misannealing and tandem primer repeats are often observed in mutagenesis cloning and reduce the efficiency of mutagenesis. Here we present a high-throughput mutagenesis pipeline based on established methods that significantly reduces PCR artefacts. We combined a two-fragment PCR approach, in which mutagenesis primers are used in two separate PCR reactions, with an in vitro assembly of resulting fragments. We show that this approach, despite being more laborious, is a very efficient pipeline for the creation of large libraries of mutants.

Discovery of functional monoclonal antibodies targeting G-protein-coupled receptors and ion channels

The development of recombinant antibody therapeutics is a significant area of growth in the pharmaceutical industry with almost 50 approved monoclonal antibodies on the market in the US and Europe. Despite this growth, however, certain classes of important molecular targets have remained intractable to therapeutic antibodies due to complexity of the target molecules. These complex target molecules include G-protein-coupled receptors and ion channels which represent a large potential target class for therapeutic intervention with monoclonal antibodies. Although these targets have typically been addressed by small molecule approaches, the exquisite specificity of antibodies provides a significant opportunity to provide selective modulation of these target proteins. Given this opportunity, substantial effort has been applied to address the technical challenges of targeting these complex membrane proteins with monoclonal antibodies. In this review recent progress made in the strategies for discovery of functional monoclonal antibodies for these challenging membrane protein targets is addressed.

Hydrogen/Deuterium Exchange Mass Spectrometry of Human Green Opsin Reveals a Conserved Pro-Pro Motif in Extracellular Loop 2 of Monostable Visual G Protein-Coupled Receptors

Opsins comprise the protein component of light sensitive G protein-coupled receptors (GPCRs) in the retina of the eye that are responsible for the transduction of light into a biochemical signal. Here, we used hydrogen/deuterium (H/D) exchange coupled with mass spectrometry to map conformational changes in green cone opsin upon light activation. We then compared these findings with those reported for rhodopsin. The extent of H/D exchange in green cone opsin was greater than in rhodopsin in the dark and bleached states, suggesting a higher structural heterogeneity for green cone opsin. Further analysis revealed that green cone opsin exists as a dimer in both dark (inactive) and bleached (active) states, and that the predicted glycosylation sites at N32 and N34 are indeed glycosylated. Comparison of deuterium uptake between inactive and active states of green cone opsin also disclosed a reduced solvent accessibility of the extracellular N-terminal region and an increased accessibility of the chromophore binding site. Increased H/D exchange at the extracellular side of transmembrane helix four (TM4) combined with an analysis of sequence alignments revealed a conserved Pro-Pro motif in extracellular loop 2 (EL2) of monostable visual GPCRs. These data present new insights into the locus of chromophore release at the extracellular side of TM4 and TM5 and provide a foundation for future functional evaluation.

Integrating structural and mutagenesis data to elucidate GPCR ligand binding

G protein-coupled receptors (GPCRs) represent the largest family of human membrane proteins, as well as drug targets. A recent boom in GPCR structural biology has provided detailed images of receptor ligand binding sites and interactions on the molecular level. An ever-increasing number of ligands is reported that exhibit activity through multiple receptors, binding in allosteric sites, and bias towards different intracellular signalling pathways. Furthermore, a wealth of single point mutants has accumulated in literature and public databases. Integrating these structural and mutagenesis data will help elucidate new GPCR ligand binding sites, and ultimately design drugs with tailored pharmacological activity.