How Can Mutations Thermostabilize G-Protein-Coupled Receptors?

Investigation of F508del CFTR unfolding and a search for stabilizing small molecules

Mutation of phenylalanine at position 508 in the cystic fibrosis transmembrane conductance regulator (F508del CFTR) yields a protein unstable at physiological temperatures that is rapidly degraded in the cell. This mutation is present in about 90% of cystic fibrosis patients, hence there is great interest in compounds reversing its instability. We have previously reported the expression of the mutated protein at low temperature and its purification in detergent. Here we describe the use of the protein to screen compounds present in a library of Federal Drug Administration (FDA) - approved drugs and also in a small natural product library. The kinetics of unfolding of F508del CFTR at 37 °C were probed by the increase in solvent-exposed cysteine residues accessible to a fluorescent reporter molecule. This occurred in a bi-exponential manner with a major (≈60%) component of half-life around 5 min and a minor component of around 60 min. The faster kinetics match those observed for loss of channel activity of F508del CFTR in cells at 37 °C. Most compounds tested had no effect on the fluorescence increase, but some were identified that significantly slowed the kinetics. The general properties of these compounds, and any likely mechanisms for inducing stability in purified CFTR are discussed. These experimental data may be useful for artificial intelligence - aided design of CFTR-specific drugs and in the identification of stabilizing additives for membrane proteins (in general).

Engineering G protein-coupled receptors for stabilization

G protein-coupled receptors (GPCRs) are one of the most important families of targets for drug discovery. One of the limiting steps in the study of GPCRs has been their stability, with significant and time-consuming protein engineering often used to stabilize GPCRs for structural characterization and drug screening. Unfortunately, computational methods developed using globular soluble proteins have translated poorly to the rational engineering of GPCRs. To fill this gap, we propose GPCR-tm, a novel and personalized structurally driven web-based machine learning tool to study the impacts of mutations on GPCR stability. We show that GPCR-tm performs as well as or better than alternative methods, and that it can accurately rank the stability changes of a wide range of mutations occurring in various types of class A GPCRs. GPCR-tm achieved Pearson's correlation coefficients of 0.74 and 0.46 on 10-fold cross-validation and blind test sets, respectively. We observed that the (structural) graph-based signatures were the most important set of features for predicting destabilizing mutations, which points out that these signatures properly describe the changes in the environment where the mutations occur. More specifically, GPCR-tm was able to accurately rank mutations based on their effect on protein stability, guiding their rational stabilization. GPCR-tm is available through a user-friendly web server at https://biosig.lab.uq.edu.au/gpcr_tm/.

Industrial and academic approaches to the search for alternative melatonin receptor ligands: An historical survey

The search for melatonin receptor agonists formed the main part of melatonin medicinal chemistry programs for the last three decades. In this short review, we summarize the two main aspects of these programs: the development of all the necessary tools to characterize the newly synthesized ligands at the two melatonin receptors MT1 and MT2, and the medicinal chemist's approaches to find chemically diverse ligands at these receptors. Both strategies are described. It turns out that the main source of tools were industrial laboratories, while the medicinal chemistry was mainly carried out in academia. Such complete accounts are interesting, as they delineate the spirits in which the teams were working demonstrating their strength and innovative character. Most of the programs were focused on nonselective agonists and few of them reached the market. In contrast, discovery of MT1-selective agonists and melatonergic antagonists with proven in vivo activity and MT1 or MT2-selectivity is still in its infancy, despite the considerable interest that subtype selective compounds may bring in the domain, as the physiological respective roles of the two subtypes of melatonin receptors, is still poorly understood. Poly-pharmacology applications and multitarget ligands have also been considered.

Rational Design of Drugs Targeting G-Protein-Coupled Receptors: A Structural Biology Perspective

G protein-coupled receptors (GPCRs) play a key role in the transduction of extracellular signals to cells and regulation of many biological processes, which makes these membrane proteins one of the most important targets for pharmacological agents. A significant increase in the number of resolved atomic structures of GPCRs has opened the possibility of developing pharmaceuticals targeting these receptors via structure-based drug design (SBDD). SBDD employs information on the structure of receptor-ligand complexes to search for selective ligands without the need for an extensive high-throughput experimental ligand screening and can significantly expand the chemical space for ligand search. In this review, we describe the process of deciphering GPCR structures using X-ray diffraction analysis and cryoelectron microscopy as an important stage in the rational design of drugs targeting this receptor class. Our main goal was to present modern developments and key features of experimental methods used in SBDD of GPCR-targeting agents to a wide range of specialists.

Characterization of RNA polymerase II trigger loop mutations using molecular dynamics simulations and machine learning

Catalysis and fidelity of multisubunit RNA polymerases rely on a highly conserved active site domain called the trigger loop (TL), which achieves roles in transcription through conformational changes and interaction with NTP substrates. The mutations of TL residues cause distinct effects on catalysis including hypo- and hyperactivity and altered fidelity. We applied molecular dynamics simulation (MD) and machine learning (ML) techniques to characterize TL mutations in the Saccharomyces cerevisiae RNA Polymerase II (Pol II) system. We did so to determine relationships between individual mutations and phenotypes and to associate phenotypes with MD simulated structural alterations. Using fitness values of mutants under various stress conditions, we modeled phenotypes along a spectrum of continual values. We found that ML could predict the phenotypes with 0.68 R2 correlation from amino acid sequences alone. It was more difficult to incorporate MD data to improve predictions from machine learning, presumably because MD data is too noisy and possibly incomplete to directly infer functional phenotypes. However, a variational auto-encoder model based on the MD data allowed the clustering of mutants with different phenotypes based on structural details. Overall, we found that a subset of loss-of-function (LOF) and lethal mutations tended to increase distances of TL residues to the NTP substrate, while another subset of LOF and lethal substitutions tended to confer an increase in distances between TL and bridge helix (BH). In contrast, some of the gain-of-function (GOF) mutants appear to cause disruption of hydrophobic contacts among TL and nearby helices.

Effect of Ligands and Transducers on the Neurotensin Receptor 1 Conformational Ensemble

Using a discrete, intracellular ¹⁹F nuclear magnetic resonance (NMR) probe on transmembrane helix 6 of the neurotensin receptor 1 (NTS1), we aim to understand how ligands and transducers modulate the receptor's structural ensemble in a solution. For apo NTS1, ¹⁹F NMR spectra reveal an ensemble of at least three conformational substates (one inactive and two active-like) in equilibrium that exchange on the millisecond to second timescale. Dynamic NMR experiments reveal that these substates follow a linear three-site exchange process that is both thermodynamically and kinetically remodeled by orthosteric ligands. As previously observed in other G protein-coupled receptors (GPCRs), the full agonist is insufficient to completely stabilize the active-like state. The inactive substate is abolished upon coupling to β-arrestin-1 (βArr1) or the C-terminal helix of Gα_q, which comprises ≳60% of the GPCR/G protein interface surface area. Whereas βArr1 exclusively selects for pre-existing active-like substates, the Gα_q peptide induces a new substate. Both transducer molecules promote substantial line broadening of active-like states, suggesting contributions from additional microsecond to millisecond exchange processes. Together, our study suggests that (i) the NTS1 allosteric activation mechanism may be alternatively dominated by induced fit or conformational selection depending on the coupled transducer, and (ii) the available static structures do not represent the entire conformational ensemble observed in a solution.

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.

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.

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.

[Activation Mechanism of Prostanoid Receptors -X-ray Crystallography of EP3 Receptor]

Prostanoids [prostaglandins (PGs) and thromboxanes (TXs)] are a series of bioactive lipid metabolites that function in an autacoid manner via activation of cognate G protein-coupled receptors (GPCRs). The nine subtypes of prostanoid receptors (DP1, DP2, EP1, EP2, EP3, EP4, FP, IP, TP) are involved in a wide range of functions, including inflammation, immune response, reproduction, and homeostasis of the intestinal mucosa and cardiovascular system. Among the prostanoid receptors, the structure of antagonist-bound DP2, which belongs to the chemoattractant receptor family, was previously determined. However, the mechanisms of prostanoid recognition and receptor activation remained elusive. To address this issue, we determined the crystal structures of antagonist-bound EP4 and PGE₂-bound EP3. The EP3-PGE₂ complex exhibits an active-like conformation, including outward movement of the cytoplasmic end of transmembrane (TM) 6 relative to the cytoplasmic end of TM6 of the EP4 complex. The carboxyl moiety of PGE₂ is recognized through three hydrogen bonds formed by highly conserved residues: Y114^_2.65, T206^{Extracelluar loop 2 (ECL2)}, and R333^7.40 (superscripts denote Ballesteros-Weinstein numbering). In addition, the ω-chain of PGE₂ orients toward TM6, which appears to contribute to receptor activation. The structure reveals important insights into the activation mechanism of prostanoid receptors and provides a molecular basis for the binding modes of endogenous ligands. These findings should facilitate the development of subtype-selective and non-PG-like ligands.

Fine Sampling of Sequence Space for Membrane Protein Structural Biology

We describe an enhancement of traditional genomics-based approaches to improve the success of structure determination of membrane proteins. Following a broad screen of sequence space to identify initial expression-positive targets, we employ a second step to select orthologs with closely related sequences to these hits. We demonstrate that a greater percentage of these latter targets express well and are stable in detergent, increasing the likelihood of identifying candidates that will ultimately yield structural information.

Molecular Dynamics Simulations for the Determination of the Characteristic Structural Differences between Inactive and Active States of Wild Type and Mutants of the Orexin2 Receptor

The orexin2 receptor (OX2R), which is classified as a class A G protein-coupled receptor (GPCR), is the target of our study. We performed over 20 several-microsecond-scale molecular dynamics simulations of the wild type and mutants of OX2R to extract the characteristics of the structural changes taking place in the active state. We introduced mutations that exhibited the stable inactive state and the constitutively active state in class A GPCRs. In these simulations, significant characteristic structural changes were observed in the V309^6.40Y mutant, which corresponded to a constitutively active mutant. These conformational changes include the outward movement of the transmembrane helix 6 (TM6) and the inward movement of TM7, which are common structural changes in the activation of GPCRs. In addition, we extracted a suitable index for the quantitative evaluation of the active and inactive states of GPCRs, namely, the inter-atomic distance of Cα atoms between x(3.46) and Y(7.53). The structures of the inactive and active states solved by X-ray crystallography and cryo-electron microscopy can be classified using the inter-atomic distance. Furthermore, we clarified that the inward movement of TM7 requires the swapping of M305^6.36 on TM6 and L367^7.56 on TM7. Finally, we discussed the structural advantages of TM7 inward movement for GPCR activation.

A universal allosteric mechanism for G protein activation

G proteins play a central role in signal transduction and pharmacology. Signaling is initiated by cell-surface receptors, which promote guanosine triphosphate (GTP) binding and dissociation of Gα from the Gβγ subunits. Structural studies have revealed the molecular basis of subunit association with receptors, RGS proteins, and downstream effectors. In contrast, the mechanism of subunit dissociation is poorly understood. We use cell signaling assays, molecular dynamics (MD) simulations, and biochemistry and structural analyses to identify a conserved network of amino acids that dictates subunit release. In the presence of the terminal phosphate of GTP, a glycine forms a polar network with an arginine and glutamate, putting torsional strain on the subunit binding interface. This "G-R-E motif" secures GTP and, through an allosteric link, discharges the Gβγ dimer. Replacement of network residues prevents subunit dissociation regardless of agonist or GTP binding. These findings reveal the molecular basis of the final committed step of G protein activation.

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.

Loss of a water-mediated network results in reduced agonist affinity in a β2-adrenergic receptor clinical variant

The β₂-adrenergic receptor (β₂AR) is a member of the G protein-coupled receptor (GPCR) family that is an important drug target for asthma and COPD. Clinical studies coupled with biochemical data have identified a critical receptor variant, Thr164Ile, to have a reduced response to agonist-based therapy, although the molecular mechanism underlying this seemingly "non-deleterious" substitution is not clear. Here, we couple molecular dynamics simulations with network analysis and free-energy calculations to identify the molecular determinants underlying the differential drug response. We are able to identify hydration sites in the transmembrane domain that are essential to maintain the integrity of the binding site but are absent in the variant. The loss of these hydration sites in the variant correlates with perturbations in the intra-protein interaction network and rearrangements in the orthosteric ligand binding site. In conjunction, we observe an altered binding and reduced free energy of a series of agonists, in line with experimental trends. Our work identifies a functional allosteric pathway connected by specific hydration sites in β₂AR that has not been reported before and provides insight into water-mediated networks in GPCRs in general. Overall, the work is one of the first step towards developing variant-specific potent and selective agonists.

Fake It 'Till You Make It-The Pursuit of Suitable Membrane Mimetics for Membrane Protein Biophysics

Membrane proteins evolved to reside in the hydrophobic lipid bilayers of cellular membranes. Therefore, membrane proteins bridge the different aqueous compartments separated by the membrane, and furthermore, dynamically interact with their surrounding lipid environment. The latter not only stabilizes membrane proteins, but directly impacts their folding, structure and function. In order to be characterized with biophysical and structural biological methods, membrane proteins are typically extracted and subsequently purified from their native lipid environment. This approach requires that lipid membranes are replaced by suitable surrogates, which ideally closely mimic the native bilayer, in order to maintain the membrane proteins structural and functional integrity. In this review, we survey the currently available membrane mimetic environments ranging from detergent micelles to bicelles, nanodiscs, lipidic-cubic phase (LCP), liposomes, and polymersomes. We discuss their respective advantages and disadvantages as well as their suitability for downstream biophysical and structural characterization. Finally, we take a look at ongoing methodological developments, which aim for direct in-situ characterization of membrane proteins within native membranes instead of relying on membrane mimetics.

Structure and Dynamics of GPCRs in Lipid Membranes: Physical Principles and Experimental Approaches

Over the past decade, the vast amount of information generated through structural and biophysical studies of GPCRs has provided unprecedented mechanistic insight into the complex signalling behaviour of these receptors. With this recent information surge, it has also become increasingly apparent that in order to reproduce the various effects that lipids and membranes exert on the biological function for these allosteric receptors, in vitro studies of GPCRs need to be conducted under conditions that adequately approximate the native lipid bilayer environment. In the first part of this review, we assess some of the more general effects that a membrane environment exerts on lipid bilayer-embedded proteins such as GPCRs. This is then followed by the consideration of more specific effects, including stoichiometric interactions with specific lipid subtypes. In the final section, we survey a range of different membrane mimetics that are currently used for in vitro studies, with a focus on NMR applications.

IMPROvER: the Integral Membrane Protein Stability Selector

Identifying stabilising variants of membrane protein targets is often required for structure determination. Our new computational pipeline, the Integral Membrane Protein Stability Selector (IMPROvER) provides a rational approach to variant selection by employing three independent approaches: deep-sequence, model-based and data-driven. In silico tests using known stability data, and in vitro tests using three membrane protein targets with 7, 11 and 16 transmembrane helices provided measures of success. In vitro, individual approaches alone all identified stabilising variants at a rate better than expected by random selection. Low numbers of overlapping predictions between approaches meant a greater success rate was achieved (fourfold better than random) when approaches were combined and selections restricted to the highest ranked sites. The mix of information IMPROvER uses can be extracted for any helical membrane protein. We have developed the first general-purpose tool for selecting stabilising variants of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\upalpha$$\end{document}α-helical membrane proteins, increasing efficiency and reducing workload. IMPROvER can be accessed at http://improver.ddns.net/IMPROvER/.

Releasing the technical 'shackles' on GPCR drug discovery: opportunities enabled by detergent-free polymer lipid particle (PoLiPa) purification

G-protein-coupled receptor (GPCR) drug research is presently hindered by the technical challenges associated with generating purified receptors. Consequently, the application of critical modern discovery technologies has been limited, and the vast untapped opportunity for new GPCR-directed medicines is not being realised. A simple but transformative solution is to purify receptors without removing them from their native phospholipid environment by using polymer lipid particle (PoLiPa) technology, with reagents such as styrene-maleic acid co-polymer (SMA). Compared with contemporary detergent-based and stabilising mutagenesis methods, the PoLiPa approach is simple and generic and, therefore, offers huge advantages, with the potential to revolutionise GPCR research by facilitating the availability of the purified receptors that are required for structural biology, biophysical, and panning technologies.

Cryo-EM structure of an activated VIP1 receptor-G protein complex revealed by a NanoBiT tethering strategy

Vasoactive intestinal polypeptide receptor (VIP1R) is a widely expressed class B G protein-coupled receptor and a drug target for the treatment of neuronal, metabolic, and inflammatory diseases. However, our understanding of its mechanism of action and the potential of drug discovery targeting this receptor is limited by the lack of structural information of VIP1R. Here we report a cryo-electron microscopy structure of human VIP1R bound to PACAP27 and Gs heterotrimer, whose complex assembly is stabilized by a NanoBiT tethering strategy. Comparison with other class B GPCR structures reveals that PACAP27 engages VIP1R with its N-terminus inserting into the ligand binding pocket at the transmembrane bundle of the receptor, which subsequently couples to the G protein in a receptor-specific manner. This structure has provided insights into the molecular basis of PACAP27 binding and VIP receptor activation. The methodology of the NanoBiT tethering may help to provide structural information of unstable complexes.

Vasoactive intestinal polypeptide receptor (VIP1R) is a widely expressed class B G protein-coupled receptor and a drug target for the treatment of inflammatory diseases. Here authors report a cryoelectron microscopy structure of human VIP1R bound to PACAP27 and Gs heterotrimer, which provides insights into PACAP27 binding and VIP receptor activation.

Improving the Accuracy of Protein Thermostability Predictions for Single Point Mutations

Accurately predicting the protein thermostability changes upon single point mutations in silico is a challenge that has implications for understanding diseases as well as industrial applications of protein engineering. Free energy perturbation (FEP) has been applied to predict the effect of single point mutations on protein stability for over 40 years and emerged as a potentially reliable prediction method with reasonable throughput. However, applications of FEP in protein stability calculations in industrial settings have been hindered by a number of limitations, including the inability to model mutations to and from prolines in which the bonded topology of the backbone is modified and the complexity in modeling charge-changing mutations. In this study, we have extended the FEP+ protocol to enable the accurate modeling of the effects on protein stability from proline mutations and from charge-changing mutations. We also evaluated the influence of the unfolded model in the stability calculations using increasingly longer peptides with native sequence and conformations. With the abovementioned improvements, the accuracy of FEP predictions of protein stability over a data set of 87 mutations on five different proteins has drastically improved compared with previous studies, with a mean unsigned error of 0.86 kcal/mol and root mean square error of 1.11 kcal/mol, comparable with the accuracy of previously published state-of-the-art small-molecule relative binding affinity calculations, which have been shown to be capable of driving discovery projects.

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.

ABCG2 requires a single aromatic amino acid to "clamp" substrates and inhibitors into the binding pocket

ATP-binding cassette sub-family G member 2 (ABCG2) is a homodimeric ATP-binding cassette (ABC) transporter that not only has a key role in helping cancer cells to evade the cytotoxic effects of chemotherapy, but also in protecting organisms from multiple xeno- and endobiotics. Structural studies indicate that substrate and inhibitor (ligands) binding to ABCG2 can be differentiated quantitatively by the number of amino acid contacts, with inhibitors displaying more contacts. Although binding is the obligate initial step in the transport cycle, there is no empirical evidence for one amino acid being primarily responsible for ligand binding. By mutagenesis and biochemical studies, we demonstrated that the phylogenetically conserved amino acid residue, F439, was critical for both transport and the binding of multiple substrates and inhibitors. Structural modeling implied that the π-π interactions from each F439 monomer mediated the binding of a surprisingly diverse array of structurally unrelated substrates and inhibitors and that this symmetrical π-π interaction "clamps" the ligand into the binding pocket. Key molecular features of diverse ABCG2 ligands using the π-π clamp along with structural studies created a pharmacophore model. These novel findings have important therapeutic implications because key properties of ligands interacting with ABCG2 have been disovered. Furthermore, mechanistic insights have been revealed by demonstrating that for ABCG2 a single amino acid is essential for engaging and initiating transport of multiple drugs and xenobiotics.

Bridging non-overlapping reads illuminates high-order epistasis between distal protein sites in a GPCR

Epistasis emerges when the effects of an amino acid depend on the identities of interacting residues. This phenomenon shapes fitness landscapes, which have the power to reveal evolutionary paths and inform evolution of desired functions. However, there is a need for easily implemented, high-throughput methods to capture epistasis particularly at distal sites. Here, we combine deep mutational scanning (DMS) with a straightforward data processing step to bridge reads in distal sites within genes (BRIDGE). We use BRIDGE, which matches non-overlapping reads to their cognate templates, to uncover prevalent epistasis within the binding pocket of a human G protein-coupled receptor (GPCR) yielding variants with 4-fold greater affinity to a target ligand. The greatest functional improvements in our screen result from distal substitutions and substitutions that are deleterious alone. Our results corroborate findings of mutational tolerance in GPCRs, even in conserved motifs, but reveal inherent constraints restricting tolerated substitutions due to epistasis.

Epistasis effects among amino acids at distal sites within binding pockets can have important impacts on protein fitness landscapes. Here the authors present BRIDGE, which matches non-overlapping sequence reads with their cognate DNA templates.

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.

Machine Learning for Prioritization of Thermostabilizing Mutations for G-Protein Coupled Receptors

Although the three-dimensional structures of G-protein coupled receptors (GPCRs), the largest superfamily of drug targets, have enabled structure-based drug design, there are no structures available for 87% of GPCRs. This is due to the stiff challenge in purifying the inherently flexible GPCRs. Identifying thermostabilized mutant GPCRs via systematic alanine scanning mutations has been a successful strategy in stabilizing GPCRs, but it remains a daunting task for each GPCR. We developed a computational method that combines sequence-, structure-, and dynamics-based molecular properties of GPCRs that recapitulate GPCR stability, with four different machine learning methods to predict thermostable mutations ahead of experiments. This method has been trained on thermostability data for 1231 mutants, the largest publicly available data set. A blind prediction for thermostable mutations of the complement factor C5a receptor 1 retrieved 36% of the thermostable mutants in the top 50 prioritized mutants compared to 3% in the first 50 attempts using systematic alanine scanning.

Prediction of Conformation Specific Thermostabilizing Mutations for Class A G Protein-Coupled Receptors

G protein-coupled receptors (GPCRs) are highly flexible and prone to denaturation during protein extraction in detergents and purification. This poses a huge challenge to purify a conformationally homogeneous solution of GPCRs. Thermostabilizing mutations have been used widely to purify and obtain crystal structures of several GPCRs. However, identifying thermostabilizing mutations for GPCRs remains a tedious and expensive task as they are not transferable even among closely related GPCRs. Additionally, the mutations stabilizing one conformational state of a GPCR do not always stabilize other conformational state(s) of the same GPCR. Previously we developed a computational method, LiticonDesign, for rapid prediction of thermostabilizing mutations for a specific GPCR conformation. In this study, we have used LiticonDesign to predict thermostabilizing mutations for the agonist bound active-intermediate state of the human adenosine receptor (A2AR) using the structure of the inactive state of the same GPCR and vice versa. Our study shows that the thermostable mutation predictions using LiticonDesign, for an active-intermediate state of a GPCR (A2AR in our case), requires a homology model that is derived from an active/active-intermediate state GPCR structure as a template. Similarly, the homology models derived from inactive state GPCR conformations are better in predicting the thermostable mutations for the inactive state of A2AR. Overall, LiticonDesign method is not only efficient in predicting thermostabilizing mutations for a given GPCR sequence but also can recover conformation specific mutations for a state of interest, if a suitable starting structure of desired conformation is chosen.

Selection of Biophysical Methods for Characterisation of Membrane Proteins

Over the years, there have been many developments and advances in the field of integral membrane protein research. As important pharmaceutical targets, it is paramount to understand the mechanisms of action that govern their structure-function relationships. However, the study of integral membrane proteins is still incredibly challenging, mostly due to their low expression and instability once extracted from the native biological membrane. Nevertheless, milligrams of pure, stable, and functional protein are always required for biochemical and structural studies. Many modern biophysical tools are available today that provide critical information regarding to the characterisation and behaviour of integral membrane proteins in solution. These biophysical approaches play an important role in both basic research and in early-stage drug discovery processes. In this review, it is not our objective to present a comprehensive list of all existing biophysical methods, but a selection of the most useful and easily applied to basic integral membrane protein research.

Folding and Misfolding of Human Membrane Proteins in Health and Disease: From Single Molecules to Cellular Proteostasis

Advances over the past 25 years have revealed much about how the structural properties of membranes and associated proteins are linked to the thermodynamics and kinetics of membrane protein (MP) folding. At the same time biochemical progress has outlined how cellular proteostasis networks mediate MP folding and manage misfolding in the cell. When combined with results from genomic sequencing, these studies have established paradigms for how MP folding and misfolding are linked to the molecular etiologies of a variety of diseases. This emerging framework has paved the way for the development of a new class of small molecule "pharmacological chaperones" that bind to and stabilize misfolded MP variants, some of which are now in clinical use. In this review, we comprehensively outline current perspectives on the folding and misfolding of integral MPs as well as the mechanisms of cellular MP quality control. Based on these perspectives, we highlight new opportunities for innovations that bridge our molecular understanding of the energetics of MP folding with the nuanced complexity of biological systems. Given the many linkages between MP misfolding and human disease, we also examine some of the exciting opportunities to leverage these advances to address emerging challenges in the development of therapeutics and precision medicine.

The role of NMR spectroscopy in mapping the conformational landscape of GPCRs

Over recent years, nuclear magnetic resonance (NMR) spectroscopy has developed into a powerful mechanistic tool for the investigation of G protein-coupled receptors (GPCRs). NMR provides insights which underpin the dynamic nature of these important receptors and reveals experimental evidence for a complex conformational energy landscape that is explored during receptor activation resulting in signalling. NMR studies have highlighted both the dynamic properties of different receptor states as well as the exchange pathways and intermediates formed during activation, extending the static view of GPCRs obtained from other techniques. NMR studies can be undertaken in realistic membrane-like phospholipid environments and an ever-increasing choice of labelling strategies provides comprehensive, receptor-wide information. Combined with other structural methods, NMR is contributing to our understanding of allosteric signal propagation and the interaction of GPCRs with intracellular binding partners (IBP), crucial to explaining cellular signalling.

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.

Solution NMR spectroscopy of GPCRs: Residue-specific labeling strategies with a focus on 13C-methyl methionine labeling of the atypical chemokine receptor ACKR3

The past decade has witnessed remarkable progress in the determination of G protein-coupled receptor (GPCR) structures, profoundly expanding our understanding of how GPCRs recognize ligands, become activated, and interact with intracellular signaling components. In recent years, numerous studies have used solution nuclear magnetic resonance (NMR) spectroscopy to investigate GPCRs, providing fundamental insights into GPCR conformational changes, allostery, dynamics, and other facets of GPCR function are challenging to study using other structural techniques. Despite these advantages, NMR-based studies of GPCRs are few relative to the number of published structures, due in part to the challenges and limitations of NMR for the characterization of large membrane proteins. Several studies have circumvented these challenges using a variety of isotopic labeling strategies, including side chain derivatization and metabolic incorporation of NMR-active nuclei. In this chapter, we provide an overview of different isotopic labeling strategies and describe an in-depth protocol for the expression, purification, and NMR studies of the chemokine GPCR atypical chemokine receptor 3 (ACKR3) via ¹³CH₃-methionine incorporation. The goal of this chapter is to provide a resource to the GPCR community for those interested in pursuing NMR studies of GPCRs.

Structural insights into the subtype-selective antagonist binding to the M2 muscarinic receptor

Human muscarinic receptor, M₂ is one of the five subtypes of muscarinic receptors belonging to the family of G protein-coupled receptors. Muscarinic receptors are targets for multiple neurodegenerative diseases. The challenge has been designing subtype selective ligands against one of the five muscarinic receptors. We report high resolution structures of a thermostabilized mutant M₂ receptor bound to a subtype selective antagonist AF-DX 384 and a non-selective antagonist NMS. The thermostabilizing mutation S110R in M₂ was predicted using a theoretical strategy previously developed in our group. Comparison of the crystal structures and pharmacological properties of the M₂ receptor shows that the Arg in the S110R mutant mimics the stabilizing role of the sodium cation, that is known to allosterically stabilize inactive state(s) of class A GPCRs. Molecular Dynamics simulations reveal that tightening of the ligand-residue contacts in M₂ receptor compared to M₃ receptor leads to subtype selectivity of AF-DX 384.

Engineering Salt Bridge Networks between Transmembrane Helices Confers Thermostability in G-Protein-Coupled Receptors

Introduction of specific point mutations has been an effective strategy in enhancing the thermostability of G-protein-coupled receptors (GPCRs). Our previous work showed that a specific residue position on transmembrane helix 6 (TM6) in class A GPCRs consistently yields thermostable mutants. The crystal structure of human chemokine receptor CCR5 also showed increased thermostability upon mutation of two positions, A233D^6.33 and K303E^7.59. With the goal of testing the transferability of these two thermostabilizing mutations in other chemokine receptors, we tested the mutations A237D^6.33 and R307E^7.59 in human CCR3 for thermostability and aggregation properties in detergent solution. Interestingly, the double mutant exhibited a 6-10-fold decrease in the aggregation propensity of the wild-type protein. This is in stark contrast to the two single mutants whose aggregation properties resemble the wild type (WT). Moreover, unlike in CCR5, the two single mutants separately showed no increase in thermostability compared to the wild-type CCR3, while the double-mutant A237D^6.33/R307E^7.59 confers an increase of 2.6 °C in the melting temperature compared to the WT. Extensive all-atom molecular dynamics (MD) simulations in detergent micelles show that a salt bridge network between transmembrane helices TM3, TM6, and TM7 that is absent in the two single mutants confers stability in the double mutant. The free energy surface of the double mutant shows conformational homogeneity compared to the single mutants. An annular n-dodecyl maltoside detergent layer packs tighter to the hydrophobic surface of the double-mutant CCR3 compared to the single mutants providing additional stability. The purification of other C-C chemokine receptors lacking such stabilizing residues may benefit from the incorporation of these two point mutations.

Membrane protein engineering to the rescue

The inherent hydrophobicity of membrane proteins is a major barrier to membrane protein research and understanding. Their low stability and solubility in aqueous environments coupled with poor expression levels make them a challenging area of research. For many years, the only way of working with membrane proteins was to optimise the environment to suit the protein, through the use of different detergents, solubilising additives, and other adaptations. However, with innovative protein engineering methodologies, the membrane proteins themselves are now being adapted to suit the environment. This mini-review looks at the types of adaptations which are applied to membrane proteins from a variety of different fields, including water solubilising fusion tags, thermostabilising mutation screening, scaffold proteins, stabilising protein chimeras, and isolating water-soluble domains.

An engineered thermal-shift screen reveals specific lipid preferences of eukaryotic and prokaryotic membrane proteins

Membrane bilayers are made up of a myriad of different lipids that regulate the functional activity, stability, and oligomerization of many membrane proteins. Despite their importance, screening the structural and functional impact of lipid–protein interactions to identify specific lipid requirements remains a major challenge. Here, we use the FSEC-TS assay to show cardiolipin-dependent stabilization of the dimeric sodium/proton antiporter NhaA, demonstrating its ability to detect specific protein-lipid interactions. Based on the principle of FSEC-TS, we then engineer a simple thermal-shift assay (GFP-TS), which facilitates the high-throughput screening of lipid- and ligand- interactions with membrane proteins. By comparing the thermostability of medically relevant eukaryotic membrane proteins and a selection of bacterial counterparts, we reveal that eukaryotic proteins appear to have evolved to be more dependent to the presence of specific lipids.

Membrane bilayers are made up of a myriad of different lipids that affect membrane proteins, but identifying those specific lipid requirements remains a challenge. Here authors present an engineered thermal-shift screen which reveals specific lipid preferences of eukaryotic and prokaryotic membrane proteins.

Transporters through the looking glass. An insight into the mechanisms of ion-coupled transport and methods that help reveal them

Cell membranes, despite providing a barrier to protect intracellular constituents, require selective gating for influx of important metabolites including ions, sugars, amino acids, neurotransmitters and efflux of toxins and metabolic end-products. The machinery involved in carrying out this gating process comprises of integral membrane proteins that use ionic electrochemical gradients or ATP hydrolysis, to drive concentrative uptake or efflux. The mechanism through which ion-coupled transporters function is referred to as alternating-access. In the recent past, discrete modes of alternating-access have been described with the elucidation of new transporter structures and their snapshots in altered conformational states. Despite X-ray structures being the primary sources of mechanistic information, other biophysical methods provide information related to the structural dynamics of these transporters. Methods including EPR and smFRET, have extensively helped validate or clarify ion-coupled transport mechanisms, in a near-native environment. This review seeks to highlight the mechanistic details of ion-coupled transport and delve into the biophysical tools and methods that help in understanding these fascinating molecules.

Screening of Detergents for Stabilization of Functional Membrane Proteins

Membrane protein studies usually require use of detergents to extract and isolate proteins from membranes and manipulate them in a soluble context for their functional or structural characterization. However, solubilization with detergent may interfere with MP stability and may directly affect MP function or structure. Moreover, detergent properties can be affected such as critical micellar concentration (CMC) can be affected by the experimental conditions. Consequently, the experimenter must pay attention to both the protein and the behavior of the detergent. This article provides a convenient protocol for estimating the CMC of detergents in given experimental conditions. Then, it presents two protocols aimed at monitoring the function of a membrane protein in the presence of detergent. Such experiments may help to test various detergents for their inactivating or stabilizing effects on long incubation times, ranging from few hours to some days. © 2018 by John Wiley & Sons, Inc.

Reprogramming G protein coupled receptor structure and function

The prominence of G protein-coupled receptors (GPCRs) in human physiology and disease has resulted in their intense study in various fields of research ranging from neuroscience to structural biology. With over 800 members in the human genome and their involvement in a myriad of diseases, GPCRs are the single largest family of drug targets, and an ever-present interest exists in further drug discovery and structural characterization efforts. However, low GPCR expression and stability outside the natural lipid environments have challenged these efforts. In vivo functional studies of GPCR signaling are complicated not only by the need for specific spatiotemporal activation, but also by downstream effector promiscuity. In this review, we summarize the present and emerging GPCR engineering methods that have been employed to overcome the challenges involved in receptor characterization, and to better understand the functional role of these receptors.

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.

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.

Membrane properties that shape the evolution of membrane enzymes

Spectacular recent progress in structural biology has led to determination of the structures of many integral membrane enzymes that catalyze reactions in which at least one substrate also is membrane bound. A pattern of results seems to be emerging in which the active site chemistry of these enzymes is usually found to be analogous to what is observed for water soluble enzymes catalyzing the same reaction types. However, in light of the chemical, structural, and physical complexity of cellular membranes plus the presence of transmembrane gradients and potentials, these enzymes may be subject to membrane-specific regulatory mechanisms that are only now beginning to be uncovered. We review the membrane-specific environmental traits that shape the evolution of membrane-embedded biocatalysts.

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 Basis for G Protein-Coupled Receptor Signaling

G protein-coupled receptors (GPCRs), which mediate processes as diverse as olfaction and maintenance of metabolic homeostasis, have become the single most effective class of therapeutic drug targets. As a result, understanding the molecular basis for their activity is of paramount importance. Recent technological advances have made GPCR structural biology increasingly tractable, offering views of these receptors in unprecedented atomic detail. Structural and biophysical data have shown that GPCRs function as complex allosteric machines, communicating ligand-binding events through conformational change. Changes in receptor conformation lead to activation of effector proteins, such as G proteins and arrestins, which are themselves conformational switches. Here, we review how structural biology has illuminated the agonist-induced cascade of conformational changes that culminate in a cellular response to GPCR activation.

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.

Structural Basis for G Protein-Coupled Receptor Activation

G protein-coupled receptors (GPCRs) are critical regulators of human physiology and make up the largest single class of therapeutic drug targets. Although GPCRs regulate highly diverse physiology, they share a common signaling mechanism whereby extracellular stimuli induce conformational changes in the receptor that enable activation of heterotrimeric G proteins and other intracellular effectors. Advances in GPCR structural biology have made it possible to examine ligand-induced GPCR activation at an unprecedented level of detail. Here, we review the structural basis for family A GPCR activation, with a focus on GPCRs for which structures are available in both active or active-like states and inactive states. Crystallographic and other biophysical data show how chemically diverse ligands stabilize highly conserved conformational changes on the intracellular side of the receptors, allowing many different extracellular stimuli to utilize shared downstream signaling molecules. Finally, we discuss the remaining challenges in understanding GPCR activation and signaling and highlight new technologies that may allow unanswered questions to be resolved.