Structure and dynamics of a constitutively active neurotensin receptor

G Protein-Coupled Receptor-Ligand Pose and Functional Class Prediction

G protein-coupled receptor (GPCR) transmembrane protein family members play essential roles in physiology. Numerous pharmaceuticals target GPCRs, and many drug discovery programs utilize virtual screening (VS) against GPCR targets. Improvements in the accuracy of predicting new molecules that bind to and either activate or inhibit GPCR function would accelerate such drug discovery programs. This work addresses two significant research questions. First, do ligand interaction fingerprints provide a substantial advantage over automated methods of binding site selection for classical docking? Second, can the functional status of prospective screening candidates be predicted from ligand interaction fingerprints using a random forest classifier? Ligand interaction fingerprints were found to offer modest advantages in sampling accurate poses, but no substantial advantage in the final set of top-ranked poses after scoring, and, thus, were not used in the generation of the ligand-receptor complexes used to train and test the random forest classifier. A binary classifier which treated agonists, antagonists, and inverse agonists as active and all other ligands as inactive proved highly effective in ligand function prediction in an external test set of GPR31 and TAAR2 candidate ligands with a hit rate of 82.6% actual actives within the set of predicted actives.

Bayesian network models identify cooperative GPCR:G protein interactions that contribute to G protein coupling

Cooperative interactions in protein-protein interfaces demonstrate the interdependency or the linked network-like behavior and their effect on the coupling of proteins. Cooperative interactions also could cause ripple or allosteric effects at a distance in protein-protein interfaces. Although they are critically important in protein-protein interfaces, it is challenging to determine which amino acid pair interactions are cooperative. In this work, we have used Bayesian network modeling, an interpretable machine learning method, combined with molecular dynamics trajectories to identify the residue pairs that show high cooperativity and their allosteric effect in the interface of G protein-coupled receptor (GPCR) complexes with Gα subunits. Our results reveal six GPCR:Gα contacts that are common to the different Gα subtypes and show strong cooperativity in the formation of interface. Both the C terminus helix5 and the core of the G protein are codependent entities and play an important role in GPCR coupling. We show that a promiscuous GPCR coupling to different Gα subtypes, makes all the GPCR:Gα contacts that are specific to each Gα subtype (Gαs, Gαi, and Gαq). This work underscores the potential of data-driven Bayesian network modeling in elucidating the intricate dependencies and selectivity determinants in GPCR:G protein complexes, offering valuable insights into the dynamic nature of these essential cellular signaling components.

Insights on the G protein-coupled receptor helix 8 solution structure and orientation using a neurotensin receptor 1 peptide

G-protein coupled receptors (GPCRs) are the largest class of membrane proteins encoded in the human genome with high pharmaceutical relevance and implications to human health. These receptors share a prevalent architecture of seven transmembrane helices followed by an intracellular, amphipathic helix 8 (H8) and a disordered C-terminal tail (Ctail). Technological advancements have led to over 1000 receptor structures in the last two decades, yet frequently H8 and the Ctail are conformationally heterogeneous or altogether absent. Here we synthesize a peptide comprising the neurotensin receptor 1 (NTS1) H8 and Ctail (H8-Ctail) to investigate its structural stability, conformational dynamics, and orientation in the presence of detergent and phospholipid micelles, which mimic the membrane. Circular dichroism (CD) and nuclear magnetic resonance (NMR) measurements confirm that zwitterionic 1,2-diheptanoyl-sn-glycero-3-phosphocholine is a potent stabilizer of H8 structure, whereas the commonly-used branched detergent lauryl maltose neopentyl glycol (LMNG) is unable to completely stabilize the helix - even at amounts four orders of magnitude greater than its critical micellar concentration. We then used NMR spectroscopy to assign the backbone chemical shifts. A series of temperature and lipid titrations were used to define the H8 boundaries as F376-R392 from chemical shift perturbations, changes in resonance intensity, and chemical-shift-derived phi/psi angles. Finally, the H8 azimuthal and tilt angles, defining the helix orientation relative of the membrane normal were measured using paramagnetic relaxation enhancement NMR. Taken together, our studies reveal the H8-Ctail region is sensitive to membrane physicochemical properties and is capable of more adaptive behavior than previously suggested by static structural techniques.

A peroxidase-derived ligand that induces Fusarium graminearum Ste2 receptor-dependent chemotropism

Introduction

The fungal G protein-coupled receptors Ste2 and Ste3 are vital in mediating directional hyphal growth of the agricultural pathogen Fusarium graminearum towards wheat plants. This chemotropism is induced by a catalytic product of peroxidases secreted by the wheat. Currently, the identity of this product, and the substrate it is generated from, are not known.

Methods and results

We provide evidence that a peroxidase substrate is derived from F. graminearum conidia and report a simple method to extract and purify the FgSte2-activating ligand for analyses by mass spectrometry. The mass spectra arising from t he ligand extract are characteristic of a 400 Da carbohydrate moiety. Consistent with this type of molecule, glycosidase treatment of F. graminearum conidia prior to peroxidase treatment significantly reduced the amount of ligand extracted. Interestingly, availability of the peroxidase substrate appears to depend on the presence of both FgSte2 and FgSte3, as knockout of one or the other reduces the chemotropism-inducing effect of the extracts.

Conclusions

While further characterization is necessary, identification of the F. graminearum-derived peroxidase substrate and the FgSte2-activating ligand will unearth deeper insights into the intricate mechanisms that underlie fungal pathogenesis in cereal crops, unveiling novel avenues for inhibitory interventions.

Single-molecule visualization of human A2A adenosine receptor activation by a G protein and constitutively activating mutations

Mutations that constitutively activate G protein-coupled receptors (GPCRs), known as constitutively activating mutations (CAMs), modify cell signaling and interfere with drugs, resulting in diseases with limited treatment options. We utilize fluorescence imaging at the single-molecule level to visualize the dynamic process of CAM-mediated activation of the human A2A adenosine receptor (A2AAR) in real time. We observe an active-state population for all CAMs without agonist stimulation. Importantly, activating mutations significantly increase the population of an intermediate state crucial for receptor activation, notably distinct from the addition of a partner G protein. Activation kinetics show that while CAMs increase the frequency of transitions to the intermediate state, mutations altering sodium sensitivity increase transitions away from it. These findings indicate changes in GPCR function caused by mutations may be predicted based on whether they favor or disfavor formation of an intermediate state, providing a framework for designing receptors with altered functions or therapies that target intermediate states.

Bayesian network models identify co-operative GPCR:G protein interactions that contribute to G protein coupling

Cooperative interactions in protein-protein interfaces demonstrate the interdependency or the linked network-like behavior of interface interactions and their effect on the coupling of proteins. Cooperative interactions also could cause ripple or allosteric effects at a distance in protein-protein interfaces. Although they are critically important in protein-protein interfaces it is challenging to determine which amino acid pair interactions are cooperative. In this work we have used Bayesian network modeling, an interpretable machine learning method, combined with molecular dynamics trajectories to identify the residue pairs that show high cooperativity and their allosteric effect in the interface of G protein-coupled receptor (GPCR) complexes with G proteins. Our results reveal a strong co-dependency in the formation of interface GPCR:G protein contacts. This observation indicates that cooperativity of GPCR:G protein interactions is necessary for the coupling and selectivity of G proteins and is thus critical for receptor function. We have identified subnetworks containing polar and hydrophobic interactions that are common among multiple GPCRs coupling to different G protein subtypes (Gs, Gi and Gq). These common subnetworks along with G protein-specific subnetworks together confer selectivity to the G protein coupling. This work underscores the potential of data-driven Bayesian network modeling in elucidating the intricate dependencies and selectivity determinants in GPCR:G protein complexes, offering valuable insights into the dynamic nature of these essential cellular signaling components.

NMR sample optimization and backbone assignment of a stabilized neurotensin receptor

G protein-coupled receptors (GPCRs) are involved in a multitude of cellular signaling cascades and consequently are a prominent target for pharmaceutical drugs. In the past decades, a growing number of high-resolution structures of GPCRs has been solved, providing unprecedented insights into their mode of action. However, knowledge on the dynamical nature of GPCRs is equally important for a better functional understanding, which can be obtained by NMR spectroscopy. Here, we employed a combination of size exclusion chromatography, thermal stability measurements and 2D-NMR experiments for the NMR sample optimization of the stabilized neurotensin receptor type 1 (NTR1) variant HTGH4 bound to the agonist neurotensin. We identified the short-chain lipid di-heptanoyl-glycero-phosphocholine (DH7PC) as a promising membrane mimetic for high resolution NMR experiments and obtained a partial NMR backbone resonance assignment. However, internal membrane-incorporated parts of the protein were not visible due to lacking amide proton back-exchange. Nevertheless, NMR and hydrogen deuterium exchange (HDX) mass spectrometry experiments could be used to probe structural changes at the orthosteric ligand binding site in the agonist and antagonist bound states. To enhance amide proton exchange we partially unfolded HTGH4 and observed additional NMR signals in the transmembrane region. However, this procedure led to a higher sample heterogeneity, suggesting that other strategies need to be applied to obtain high-quality NMR spectra of the entire protein. In summary, the herein reported NMR characterization is an essential step toward a more complete resonance assignment of NTR1 and for probing its structural and dynamical features in different functional states.

Screening of Membrane Protein Production by Comparison of Transient Expression in Insect and Mammalian Cells

Membrane proteins are difficult biomolecules to express and purify. In this paper, we compare the small-scale production of six selected eukaryotic integral membrane proteins in insect and mammalian cell expression systems using different techniques for gene delivery. The target proteins were C terminally fused to the green fluorescent marker protein GFP to enable sensitive monitoring. We show that the choice of expression systems makes a considerable difference to the yield and quality of the six selected membrane proteins. Virus-free transient gene expression (TGE) in insect High Five cells combined with solubilization in dodecylmaltoside plus cholesteryl hemisuccinate generated the most homogeneous samples for all six targets. Further, the affinity purification of the solubilized proteins using the Twin-Strep^® tag improved protein quality in terms of yield and homogeneity compared to His-tag purification. TGE in High Five insect cells offers a fast and economically attractive alternative to the established methods that require either baculovirus construction and the infection of the insect cells or relatively expensive transient gene expression in mammalian cells for the production of integral membrane proteins.

Neurotensin(8-13) analogs as dual NTS1 and NTS2 receptor ligands with enhanced effects on a mouse model of Parkinson's disease

The modulatory interactions between neurotensin (NT) and the dopaminergic neurotransmitter system in the brain suggest that NT may be associated with the progression of Parkinson's disease (PD). NT exerts its neurophysiological effects by interactions with the human NT receptors type 1 (hNTS1) and 2 (hNTS2). Therefore, both receptor subtypes are promising targets for the development of novel NT-based analogs for the treatment of PD. In this study, we used a virtually guided molecular modeling approach to predict the activity of NT(8-13) analogs by investigating the docking models of ligands designed for binding to the human NTS1 and NTS2 receptors. The importance of the residues at positions 8 and/or 9 for hNTS1 and hNTS2 receptor binding affinity was experimentally confirmed by radioligand binding assays. Further in vitro ADME profiling and in vivo studies revealed that, compared to the parent peptide NT(8-13), compound 10 exhibited improved stability and BBB permeability combined with a significant enhancement of the motor function and memory in a mouse model of PD. The herein reported NTS1/NTS2 dual-specific NT(8-13) analogs represent an attractive tool for the development of therapeutic strategies against PD and potentially other CNS disorders.

Protein Design Strategies for the Structural–Functional Studies of G Protein-Coupled Receptors

G protein-coupled receptors (GPCRs) are an important family of membrane proteins responsible for many physiological functions in human body. High resolution GPCR structures are required to understand their molecular mechanisms and perform rational drug design, as GPCRs play a crucial role in a variety of diseases. That is difficult to obtain for the wild-type proteins because of their low stability. In this review, we discuss how this problem can be solved by using protein design strategies developed to obtain homogeneous stabilized GPCR samples for crystallization and cryoelectron microscopy.

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.

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.

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.

Structure-Based Virtual Screening for Ligands of G Protein-Coupled Receptors: What Can Molecular Docking Do for You?

G protein-coupled receptors (GPCRs) constitute the largest family of membrane proteins in the human genome and are important therapeutic targets. During the last decade, the number of atomic-resolution structures of GPCRs has increased rapidly, providing insights into drug binding at the molecular level. These breakthroughs have created excitement regarding the potential of using structural information in ligand design and initiated a new era of rational drug discovery for GPCRs. The molecular docking method is now widely applied to model the three-dimensional structures of GPCR-ligand complexes and screen for chemical probes in large compound libraries. In this review article, we first summarize the current structural coverage of the GPCR superfamily and the understanding of receptor-ligand interactions at atomic resolution. We then present the general workflow of structure-based virtual screening and strategies to discover GPCR ligands in chemical libraries. We assess the state of the art of this research field by summarizing prospective applications of virtual screening based on experimental structures. Strategies to identify compounds with specific efficacy and selectivity profiles are discussed, illustrating the opportunities and limitations of the molecular docking method. Our overview shows that structure-based virtual screening can discover novel leads and will be essential in pursuing the next generation of GPCR drugs. SIGNIFICANCE STATEMENT: Extraordinary advances in the structural biology of G protein-coupled receptors have revealed the molecular details of ligand recognition by this large family of therapeutic targets, providing novel avenues for rational drug design. Structure-based docking is an efficient computational approach to identify novel chemical probes from large compound libraries, which has the potential to accelerate the development of drug candidates.

Nonantimicrobial Actions of Macrolides: Overview and Perspectives for Future Development

Macrolides are among the most widely prescribed broad spectrum antibacterials, particularly for respiratory infections. It is now recognized that these drugs, in particular azithromycin, also exert time-dependent immunomodulatory actions that contribute to their therapeutic benefit in both infectious and other chronic inflammatory diseases. Their increased chronic use in airway inflammation and, more recently, of azithromycin in COVID-19, however, has led to a rise in bacterial resistance. An additional crucial aspect of chronic airway inflammation, such as chronic obstructive pulmonary disease, as well as other inflammatory disorders, is the loss of epithelial barrier protection against pathogens and pollutants. In recent years, azithromycin has been shown with time to enhance the barrier properties of airway epithelial cells, an action that makes an important contribution to its therapeutic efficacy. In this article, we review the background and evidence for various immunomodulatory and time-dependent actions of macrolides on inflammatory processes and on the epithelium and highlight novel nonantibacterial macrolides that are being studied for immunomodulatory and barrier-strengthening properties to circumvent the risk of bacterial resistance that occurs with macrolide antibacterials. We also briefly review the clinical effects of macrolides in respiratory and other inflammatory diseases associated with epithelial injury and propose that the beneficial epithelial effects of nonantibacterial azithromycin derivatives in chronic inflammation, even given prophylactically, are likely to gain increasing attention in the future. SIGNIFICANCE STATEMENT: Based on its immunomodulatory properties and ability to enhance the protective role of the lung epithelium against pathogens, azithromycin has proven superior to other macrolides in treating chronic respiratory inflammation. A nonantibiotic azithromycin derivative is likely to offer prophylactic benefits against inflammation and epithelial damage of differing causes while preserving the use of macrolides as antibiotics.

Motions around conserved helical weak spots facilitate GPCR activation

G protein‐coupled receptors (GPCRs) participate in most physiological processes and are important drug targets in many therapeutic areas. Recently, many GPCR X‐ray structures became available, facilitating detailed studies of their sequence‐structure‐mobility‐function relations. We show that the functional role of many conserved GPCR sequence motifs is to create weak spots in the transmembrane helices that provide the structural plasticity necessary for ligand binding and signaling. Different receptor families use different conserved sequence motifs to obtain similar helix irregularities that allow for the same motions upon GPCR activation. These conserved motions come together to facilitate the timely release of the conserved sodium ion to the cytosol. Most GPCR crystal structures could be determined only after stabilization of the transmembrane helices by mutations that remove weak spots. These mutations often lead to diminished binding of agonists, but not antagonists, which logically agrees with the fact that large helix rearrangements occur only upon agonist binding. Upon activation, six of the seven TM helices in GPCRs undergo helix motions and/or deformations facilitated by weak spots in these helices. The location of these weak spots is much more conserved than the sequence motifs that cause them. Knowledge about these weak spots helps understand the activation process of GPCRs and thus helps design medicines.

Allosteric communication regulates ligand-specific GPCR activity

G protein-coupled receptors (GPCRs) are membrane-bound proteins that are ubiquitously expressed in many cell types and take part in mediating multiple signaling pathways. GPCRs are dynamic proteins and exist in an equilibrium between an ensemble of conformational states such as inactive and fully active states. This dynamic nature of GPCRs is one of the factors that confers their basal activity even in the absence of any ligand-mediated activation. Ligands selectively bind and stabilize a subset of the conformations from the ensemble leading to a shift in the equilibrium toward the inactive or the active state depending on the nature of the ligand. This ligand-selective effect is achieved through allosteric communication between the ligand binding site and G protein or β-arrestin coupling site. Similarly, the G protein coupling to the receptor exerts the allosteric effect on the ligand binding region leading to increased binding affinity for agonists and decreased affinity for antagonists or inverse agonists. In this review, we enumerate the current state of our understanding of the mechanism of allosteric communication in GPCRs with a specific focus on the critical role of computational methods in delineating the residues involved in allosteric communication. Analyzing allosteric communication mechanism using molecular dynamics simulations has revealed (a) a structurally conserved mechanism of allosteric communication that regulates the G protein coupling, (b) a rational structure-based approach to designing selective ligands, and (c) an approach to designing allosteric GPCR mutants that are either ligand and G protein or β-arrestin selective.

Engineering of Challenging G Protein-Coupled Receptors for Structure Determination and Biophysical Studies

Membrane proteins such as G protein-coupled receptors (GPCRs) exert fundamental biological functions and are involved in a multitude of physiological responses, making these receptors ideal drug targets. Drug discovery programs targeting GPCRs have been greatly facilitated by the emergence of high-resolution structures and the resulting opportunities to identify new chemical entities through structure-based drug design. To enable the determination of high-resolution structures of GPCRs, most receptors have to be engineered to overcome intrinsic hurdles such as their poor stability and low expression levels. In recent years, multiple engineering approaches have been developed to specifically address the technical difficulties of working with GPCRs, which are now beginning to make more challenging receptors accessible to detailed studies. Importantly, successfully engineered GPCRs are not only valuable in X-ray crystallography, but further enable biophysical studies with nuclear magnetic resonance spectroscopy, surface plasmon resonance, native mass spectrometry, and fluorescence anisotropy measurements, all of which are important for the detailed mechanistic understanding, which is the prerequisite for successful drug design. Here, we summarize engineering strategies based on directed evolution to reduce workload and enable biophysical experiments of particularly challenging GPCRs.

Design, Structural Optimization, and Characterization of the First Selective Macrocyclic Neurotensin Receptor Type 2 Non-opioid Analgesic

Neurotensin (NT) receptor type 2 (NTS2) represents an attractive target for the development of new NT-based analgesics. Here, we report the synthesis and functional in vivo characterization of the first constrained NTS2-selective macrocyclic NT analog. While most chemical optimization studies rely on the NT(8-13) fragment, we focused on NT(7-12) as a scaffold to design NTS2-selective macrocyclic peptides. Replacement of Ile12 by Leu, and Pro7/Pro10 by allylglycine residues followed by cyclization via ring-closing metathesis led to macrocycle 4, which exhibits good affinity for NTS2 (50 nM), high selectivity over NTS1 (>100 μM), and improved stability compared to NT(8-13). In vivo profiling in rats reveals that macrocycle 4 produces potent analgesia in three distinct rodent pain models, without causing the undesired effects associated with NTS1 activation. We further provide evidence of its non-opioid antinociceptive activity, therefore highlighting the strong therapeutic potential of NTS2-selective analogs for the management of acute and chronic pain.

Complexes of the neurotensin receptor 1 with small-molecule ligands reveal structural determinants of full, partial, and inverse agonism

Neurotensin receptor 1 (NTSR1) and related G protein-coupled receptors of the ghrelin family are clinically unexploited, and several mechanistic aspects of their activation and inactivation have remained unclear. Enabled by a new crystallization design, we present five new structures: apo-state NTSR1 as well as complexes with nonpeptide inverse agonists SR48692 and SR142948A, partial agonist RTI-3a, and the novel full agonist SRI-9829, providing structural rationales on how ligands modulate NTSR1. The inverse agonists favor a large extracellular opening of helices VI and VII, undescribed so far for NTSR1, causing a constriction of the intracellular portion. In contrast, the full and partial agonists induce a binding site contraction, and their efficacy correlates with the ability to mimic the binding mode of the endogenous agonist neurotensin. Providing evidence of helical and side-chain rearrangements modulating receptor activation, our structural and functional data expand the mechanistic understanding of NTSR1 and potentially other peptidergic receptors.

Preventing the Interaction between Coronaviruses Spike Protein and Angiotensin I Converting Enzyme 2: An In Silico Mechanistic Case Study on Emodin as a Potential Model Compound

Emodin, a widespread natural anthraquinone, has many biological activities including health-protective and adverse effects. Amongst beneficial effects, potential antiviral activity against coronavirus responsible for the severe acute respiratory syndrome outbreak in 2002–2003 has been described associated with the inhibition of the host cells target receptors recognition by the viral Spike protein. However, the inhibition mechanisms have not been fully characterized, hindering the rational use of emodin as a model compound to develop more effective analogues. This work investigates emodin interaction with the Spike protein to provide a mechanistic explanation of such inhibition. A 3D molecular modeling approach consisting of docking simulations, pharmacophoric analysis and molecular dynamics was used. The plausible mechanism is described as an interaction of emodin at the protein–protein interface which destabilizes the viral protein-target receptor complex. This analysis has been extended to the Spike protein of the coronavirus responsible for the current pandemic hypothesizing emodin’s functional conservation. This solid knowledge-based foothold provides a possible mechanistic rationale of the antiviral activity of emodin as a future basis for the potential development of efficient antiviral cognate compounds. Data gaps and future work on emodin-related adverse effects in parallel to its antiviral pharmacology are explored.

Conformational Changes in Tyrosine 11 of Neurotensin Are Required to Activate the Neurotensin Receptor 1

Cell-cell communication via endogenous peptides and their receptors is vital for controlling all aspects of human physiology and most peptides signal through G protein-coupled receptors (GPCRs). Disordered peptides bind GPCRs through complex modes for which there are few representative crystal structures. The disordered peptide neurotensin (NT) is a neuromodulator of classical neurotransmitters such as dopamine and glutamate, through activation of neurotensin receptor 1 (NTS1). While several experimental structures show how NT binds NTS1, details about the structural dynamics of NT during and after binding NTS1, or the role of peptide dynamics on receptor activation, remain obscure. Here saturation transfer difference (STD) NMR revealed that the binding mode of NT fragment NT10-13 is heterogeneous. Epitope maps of NT10-13 at NTS1 suggested that tyrosine 11 (Y11) samples other conformations to those observed in crystal structures of NT-bound NTS1. Molecular dynamics (MD) simulations confirmed that when NT is bound to NTS1, residue Y11 can exist in two χ1 rotameric states, gauche plus (g+) or gauche minus (g-). Since only the g+ Y11 state is observed in all the structures solved to date, we asked if the g- state is important for receptor activation. NT analogues with Y11 replaced with 7-OH-Tic were synthesized to restrain the dynamics of the side chain. P(OH-TIC)IL bound NTS1 with the same affinity as NT10-13 but did not activate NTS1, instead acted as an antagonist. This study highlights that flexibility of Y11 in NT may be required for NT activation of NTS1.

Olfactory marker protein elevates basal cAMP concentration

Olfactory marker protein (OMP), which is expressed abundantly in mature olfactory receptor neurons, operates as a cAMP-binding protein. OMP captures phasic cAMP surges induced by sensory stimuli and punctuates the downstream signalling in the cilia. On the other hand, OMP is also abundant in the soma. At equilibrium, OMP should exhibit association/dissociation reactions with cAMP. To examine the steady-state function of OMP, we expressed OMP in an HEK293 heterologous expression system and measured the activity of cAMP-dependent protein kinase (PKA) using a cAMP response element/luciferase reporter assay. In the presence of OMP, the basal activity level of PKA was elevated to approximately twice as much as that in the absence of OMP. Upon tonic stimulation by membrane-permeable cAMP, the PKA activity increased in a dose-dependent manner and was greater in the presence of OMP at all doses until saturation. These results indicate that OMP, a cytosolic cAMP-binding protein, operates as a cAMP reservoir by increases the basal cAMP concentration and enhances tonic cAMP actions. Together with the previous finding that OMP acutely sequesters cAMP-related responses, these results indicate that OMP can buffer acute surges in cAMP and tonic production, which stabilizes the basal cAMP pool in the long run.

Conformational dynamics of a G protein-coupled receptor helix 8 in lipid membranes

G protein-coupled receptors (GPCRs) are the largest and pharmaceutically most important class of membrane proteins encoded in the human genome, characterized by a seven-transmembrane helix architecture and a C-terminal amphipathic helix 8 (H8). In a minority of GPCR structures solved to date, H8 either is absent or adopts an unusual conformation. The controversial existence of H8 of the class A GPCR neurotensin receptor 1 (NTS1) has been examined here for the nonthermostabilized receptor in a functionally supporting membrane environment using electron paramagnetic resonance, molecular dynamics simulations, and circular dichroism. Lipid-protein interactions with phosphatidylserine and phosphatidylethanolamine lipids, in particular, stabilize the residues 374 to 390 of NTS1 into forming a helix. Furthermore, introduction of a helix-breaking proline residue in H8 elicited an increase in ß-arrestin-NTS1 interactions observed in pull-down assays, suggesting that the structure and/or dynamics of H8 might play an important role in GPCR signaling.

Insightful Backbone Modifications Preventing Proteolytic Degradation of Neurotensin Analogs Improve NTS1-Induced Protective Hypothermia

Therapeutic hypothermia represents a brain-protective strategy for multiple emergency situations, such as stroke or traumatic injury. Neurotensin (NT), which exerts its effects through activation of two G protein-coupled receptors, namely NTS1 and NTS2, induces a strong and long-lasting decrease in core body temperature after its central administration. Growing evidence demonstrates that NTS1 is the receptor subtype mediating the hypothermic action of NT. As such, potent NTS1 agonists designed on the basis of the minimal C-terminal NT(8-13) bioactive fragment have been shown to produce mild hypothermia and exert neuroprotective effects under various clinically relevant conditions. The high susceptibility of NT(8-13) to protease degradation (half-life <2 min) represents, however, a serious limitation for its use in pharmacological therapy. In light of this, we report here a structure-activity relationship study in which pairs of NT(8-13) analogs have been developed, based on the incorporation of a reduced Lys⁸-Lys⁹ bond. To further stabilize the peptide bonds, a panel of backbone modifications was also inserted along the peptide sequence, including Sip¹⁰, D-Trp¹¹, Dmt¹¹, Tle¹², and TMSAla¹³. Our results revealed that the combination of appropriate chemical modifications leads to compounds exhibiting improved resistance to proteolytic cleavages (>24 h; 16). Among them, the NT(8-13) analogs harboring the reduced amine bond combined with the unnatural amino acids TMSAla¹³ (4) and Sip¹⁰ (6) or the di-substitution Lys¹¹ - TMSAla¹³ (12), D-Trp¹¹-TMSAla¹³ (14), and Dmt¹¹-Tle¹² (16) produced sustained hypothermic effects (−3°C for at least 1 h). Importantly, we observed that hypothermia was mainly driven by the increased stability of the NT(8-13) derivatives, instead of the high binding-affinity at NTS1. Altogether, these results reveal the importance of the reduced amine bond in optimizing the metabolic properties of the NT(8-13) peptide and support the development of stable NTS1 agonists as first drug candidate in neuroprotective hypothermia.

A Structural Study on the Listeria Monocytogenes Internalin A-Human E-cadherin Interaction: A Molecular Tool to Investigate the Effects of Missense Mutations

Listeria monocytogenes is a widespread foodborne pathogen of high concern and internalin A is an important virulence factor that mediates cell invasion upon the interaction with the host protein E-cadherin. Nonsense mutations of internalin A are known to reduce virulence. Although missense mutations are largely overlooked, they need to be investigated in respect to their effects in cell invasion processes. This work presented a computational workflow to early characterize internalin A missense mutations. The method reliably estimated the effects of a set of engineered missense mutations in terms of their effects on internalin A–E-cadherin interaction. Then, the effects of mutations of an internalin A variant from a L. monocytogenes isolate were calculated. Mutations showed impairing effects on complex stability providing a mechanistic explanation of the low cells invasion capacity previously observed. Overall, our results provided a rational approach to explain the effects of internalin A missense mutations. Moreover, our findings highlighted that the strength of interaction may not directly relate to the cell invasion capacity reflecting the non-exclusive role of internalin A in determining the virulence of L. monocytogenes. The workflow could be extended to other virulence factors providing a promising platform to support a better molecular understanding of L. monocytogenes epidemiology.

Design of Bimodal Ligands of Neurotensin Receptor 1 for Positron Emission Tomography Imaging and Fluorescence-Guided Surgery of Pancreatic Cancer

Neurotensin receptor 1 (NTSR1) is overexpressed in most human pancreatic ductal adenocarcinomas. It makes it an attractive target for the development of pancreatic cancer imaging agents. In this study, we sought to develop a bimodal positron emission tomography (PET)/fluorescent imaging agent capable of specifically targeting these receptors. Starting from the structure of a known NTSR1 agonist, a series of tracers were synthesized, radiometalated with gallium-68, and evaluated in vitro and in vivo, in mice bearing an AsPC-1 xenograft. PET imaging allowed us to identify the compound [⁶⁸Ga]Ga-NODAGA-Lys(Cy5**)-AEEAc-[Me-Arg⁸,Tle¹²]-NT(7-13) as the one with the most promising biodistribution profile, characterized by high tumor uptake (2.56 ± 0.97%ID/g, 1 h post-injection) and rapid elimination from nontargeted organs, through urinary excretion. Fluorescence imaging gave similar results. On this basis, fluorescence-guided resection of tumor masses was successfully carried out on a preclinical model.

A complex structure of arrestin-2 bound to a G protein-coupled receptor

Arrestins comprise a family of signal regulators of G-protein-coupled receptors (GPCRs), which include arrestins 1 to 4. While arrestins 1 and 4 are visual arrestins dedicated to rhodopsin, arrestins 2 and 3 (Arr2 and Arr3) are β-arrestins known to regulate many nonvisual GPCRs. The dynamic and promiscuous coupling of Arr2 to nonvisual GPCRs has posed technical challenges to tackle the basis of arrestin binding to GPCRs. Here we report the structure of Arr2 in complex with neurotensin receptor 1 (NTSR1), which reveals an overall assembly that is strikingly different from the visual arrestin-rhodopsin complex by a 90° rotation of Arr2 relative to the receptor. In this new configuration, intracellular loop 3 (ICL3) and transmembrane helix 6 (TM6) of the receptor are oriented toward the N-terminal domain of the arrestin, making it possible for GPCRs that lack the C-terminal tail to couple Arr2 through their ICL3. Molecular dynamics simulation and crosslinking data further support the assembly of the Arr2‒NTSR1 complex. Sequence analysis and homology modeling suggest that the Arr2‒NTSR1 complex structure may provide an alternative template for modeling arrestin-GPCR interactions.

A3 adenosine receptor activation mechanisms: molecular dynamics analysis of inactive, active, and fully active states

We investigated the Gi-coupled A₃ adenosine receptor (A₃AR) activation mechanism by running 7.2 µs of molecular dynamics (MD) simulations. Based on homology to G protein-coupled receptor (GPCR) structures, three constitutively active mutant (CAM) and the wild-type (WT) A₃ARs in the apo form were modeled. Conformational signatures associated with three different receptor states (inactive R, active R*, and bound to Gi protein mimic) were predicted by analyzing and comparing the CAMs with WT receptor and by considering site-directed mutagenesis data. Detected signatures that were correlated with receptor state included: Persistent salt-bridges involving key charged residues for activation (including a novel, putative ionic lock), rotameric state of conserved W^6.48, and Na⁺ ions and water molecules present. Active-coupled state signatures similar to the X-ray structures of β₂ adrenergic receptor-Gs protein and A_2AAR-mini-Gs and the recently solved cryo-EM A₁AR-Gi complexes were found. Our MD analysis suggests that constitutive activation might arise from the D107^3.49-R108^3.50 ionic lock destabilization in R and the D107^3.49-R111^3.53 ionic lock stabilization in R* that presumably lowers the energy barrier associated with an R to R* transition. This study provides new opportunities to understand the underlying interactions of different receptor states of other Gi protein-coupled GPCRs.

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.

The nature of efficacy at G protein-coupled receptors

G protein-coupled receptors (GPCRs) participate in many pathophysiological processes as well as almost all aspects of normal physiology. They are present at the surface of all cell types making them amenable and attractive targets for pharmaceutical therapeutics. GPCRs possess complex pharmacology with the ability to be turned on to various extents, have their constitutive activity suppressed and even switch between signaling pathways to which they couple. Underlying this complex pharmacology is GPCR signaling efficacy, and differences in efficacy promoted by alternative ligands and in different tissues is of great interest to biology in general and also the pharmaceutical industry. In this review we hope to discuss what the molecular foundations of efficacy are and whether a new approach utilizing a rate-dependent model may provide new insights into this phenomenon.

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.

Conformational transitions of a neurotensin receptor 1-Gi1 complex

Neurotensin receptor 1 (NTSR1) is a G-protein-coupled receptor (GPCR) that engages multiple subtypes of G protein, and is involved in the regulation of blood pressure, body temperature, weight and the response to pain. Here we present structures of human NTSR1 in complex with the agonist JMV449 and the heterotrimeric Gi1 protein, at a resolution of 3 Å. We identify two conformations: a canonical-state complex that is similar to recently reported GPCR-Gi/o complexes (in which the nucleotide-binding pocket adopts more flexible conformations that may facilitate nucleotide exchange), and a non-canonical state in which the G protein is rotated by about 45 degrees relative to the receptor and exhibits a more rigid nucleotide-binding pocket. In the non-canonical state, NTSR1 exhibits features of both active and inactive conformations, which suggests that the structure may represent an intermediate form along the activation pathway of G proteins. This structural information, complemented by molecular dynamics simulations and functional studies, provides insights into the complex process of G-protein activation.

A benchmark study of loop modeling methods applied to G protein-coupled receptors

G protein-coupled receptors (GPCR) are important drug discovery targets. Despite progress, many GPCR structures have not yet been solved. For these targets, comparative modeling is used in virtual ligand screening to prioritize experimental efforts. However, the structure of extracellular loop 2 (ECL2) is often poorly predicted. This is significant due to involvement of ECL2 in ligand binding for many Class A GPCR. Here we examine the performance of loop modeling protocols available in the Rosetta (cyclic coordinate descent [CCD], KIC with fragments [KICF] and next generation KIC [NGK]) and Molecular Operating Environment (MOE) software suites (de novo search). ECL2 from GPCR crystal structures served as the structure prediction targets and were divided into four sets depending on loop length. Results suggest that KICF and NGK sampled and scored more loop models with sub-angstrom and near-atomic accuracy than CCD or de novo search for loops of 24 or fewer residues. None of the methods were able to sample loop conformations with near-atomic accuracy for the longest targets ranging from 25 to 32 residues based on 1000 models generated. For these long loop targets, increased conformational sampling is necessary. The strongly conserved disulfide bond between Cys3.25 and Cys45.50 in ECL2 proved an effective filter. Setting an upper limit of 5.1 Å on the S-S distance improved the lowest RMSD model included in the top 10 scored structures in Groups 1-4 on average between 0.33 and 1.27 Å. Disulfide bond formation and geometry optimization of ECL2 provided an additional incremental benefit in structure quality.

In vivo assembly and large-scale purification of a GPCR - Gα fusion with Gβγ, and characterization of the active complex

G protein coupled receptors (GPCRs) are central players in recognizing a variety of stimuli to mediate diverse cellular responses. This myriad of functions is accomplished by their modular interactions with downstream intracellular transducers, such as heterotrimeric G proteins and arrestins. Assembling a specific GPCR-G protein pair as a purified complex for their structural and functional investigations remains a challenging task, however, because of the low affinity of the interaction. Here, we optimized fusion constructs of the Gα subunit of the heterotrimeric G protein and engineered versions of rat Neurotensin receptor 1 (NTR1), coexpressed and assembled in vivo with Gβ and Gγ. This was achieved by using the baculovirus-based MultiBac system. We thus generated a functional receptor-G protein fusion complex, which can be efficiently purified using ligand-based affinity chromatography on large scales. Additionally, we utilized a purification method based on a designed ankyrin repeat protein tightly binding to Green Fluorescent Protein (GFP-DARPin) that may be used as a generic approach for a large-scale purification of GPCR-G protein fusion complexes for which no ligands column can be generated. The purification methods described herein will support future studies that aim to understand the structural and functional framework of GPCR activation and signaling.

Activation Dynamics of the Neurotensin G Protein-Coupled Receptor 1

A replica-exchange protocol remarkably enhances the sampling of the activation dynamics of the neurotensin receptor type 1, a G protein-coupled receptor (GPCR) and important drug target. Our work highlights the dynamic communication between conformational changes of the agonist and the G protein-binding site, via contraction-oscillation of the orthosteric pocket. It also gives insights into the mechanism by which certain mutations diminish or stimulate activation. The replica-exchange protocol effectively enhances barrier crossing where standard brute-force molecular dynamics simulations fail. It is readily applicable to other GPCRs and represents a promising approach for virtual ligand screening, using the typical features of receptor activation as a benchmark.

Identifying G protein-coupled receptor dimers from crystal packings

Dimers of G protein-coupled receptors (GPCRs) are believed to be important for signaling with their associated G proteins. Low-resolution electron microscopy has shown rhodopsin dimers in native retinal membranes, and CXCR4 dimers have been found in several different crystal structures. Evidence for dimers of other GPCRs is more indirect. An alternative to computational modeling studies is to search for parallel dimers in the packing environments of the reported crystal structures of GPCRs. Two major structural types of GPCR dimers exist (as predicted by others), but there is considerable structural variation within each cluster. The different structural variants described here might reflect different functional properties and should provide a range of model structures for computational and experimental examination.

Dynamic tuneable G protein-coupled receptor monomer-dimer populations

G protein-coupled receptors (GPCRs) are the largest class of membrane receptors, playing a key role in the regulation of processes as varied as neurotransmission and immune response. Evidence for GPCR oligomerisation has been accumulating that challenges the idea that GPCRs function solely as monomeric receptors; however, GPCR oligomerisation remains controversial primarily due to the difficulties in comparing evidence from very different types of structural and dynamic data. Using a combination of single-molecule and ensemble FRET, double electron–electron resonance spectroscopy, and simulations, we show that dimerisation of the GPCR neurotensin receptor 1 is regulated by receptor density and is dynamically tuneable over the physiological range. We propose a “rolling dimer” interface model in which multiple dimer conformations co-exist and interconvert. These findings unite previous seemingly conflicting observations, provide a compelling mechanism for regulating receptor signalling, and act as a guide for future physiological studies.

Evidence suggests oligomerisation of G protein-coupled receptors in membranes, but this is controversial. Here, authors use single-molecule and ensemble FRET, and spectroscopy to show that the neurotensin receptor 1 forms multiple dimer conformations that interconvert - “rolling” interfaces.

Towards monitoring conformational changes of the GPCR neurotensin receptor 1 by single-molecule FRET

Neurotensin receptor 1 (NTSR1) is a G protein-coupled receptor that is important for signaling in the brain and the gut. Its agonist ligand neurotensin (NTS), a 13-amino-acid peptide, binds with nanomolar affinity from the extracellular side to NTSR1 and induces conformational changes that trigger intracellular signaling processes. Our goal is to monitor the conformational dynamics of single fluorescently labeled NTSR1. For this, we fused the fluorescent protein mNeonGreen to the C terminus of NTSR1, purified the receptor fusion protein from E. coli membranes, and reconstituted NTSR1 into liposomes with E. coli polar lipids. Using single-molecule anisotropy measurements, NTSR1 was found to be monomeric in liposomes, with a small fraction being dimeric and oligomeric, showing homoFRET. Similar results were obtained for NTSR1 in detergent solution. Furthermore, we demonstrated agonist binding to NTSR1 by time-resolved single-molecule Förster resonance energy transfer (smFRET), using neurotensin labeled with the fluorophore ATTO594.

In Silico Studies Targeting G-protein Coupled Receptors for Drug Research Against Parkinson's Disease

Parkinson's Disease (PD) is a long-term neurodegenerative brain disorder that mainly affects the motor system. The causes are still unknown, and even though currently there is no cure, several therapeutic options are available to manage its symptoms. The development of novel antiparkinsonian agents and an understanding of their proper and optimal use are, indeed, highly demanding. For the last decades, L-3,4-DihydrOxyPhenylAlanine or levodopa (L-DOPA) has been the gold-standard therapy for the symptomatic treatment of motor dysfunctions associated to PD. However, the development of dyskinesias and motor fluctuations (wearing-off and on-off phenomena) associated with long-term L-DOPA replacement therapy have limited its antiparkinsonian efficacy. The investigation for non-dopaminergic therapies has been largely explored as an attempt to counteract the motor side effects associated with dopamine replacement therapy. Being one of the largest cell membrane protein families, G-Protein-Coupled Receptors (GPCRs) have become a relevant target for drug discovery focused on a wide range of therapeutic areas, including Central Nervous System (CNS) diseases. The modulation of specific GPCRs potentially implicated in PD, excluding dopamine receptors, may provide promising non-dopaminergic therapeutic alternatives for symptomatic treatment of PD. In this review, we focused on the impact of specific GPCR subclasses, including dopamine receptors, adenosine receptors, muscarinic acetylcholine receptors, metabotropic glutamate receptors, and 5-hydroxytryptamine receptors, on the pathophysiology of PD and the importance of structure- and ligand-based in silico approaches for the development of small molecules to target these receptors.

G protein-coupled receptors: the evolution of structural insight

G protein-coupled receptors (GPCR) comprise a diverse superfamily of over 800 proteins that have gained relevance as biological targets for pharmaceutical drug design. Although these receptors have been investigated for decades, three-dimensional structures of GPCR have only recently become available. In this review, we focus on the technological advancements that have facilitated efforts to gain insights into GPCR structure. Progress in these efforts began with the initial crystal structure determination of rhodopsin (PDB: 1F88) in 2000 and has continued to the most recently published structure of the A_1AR (PDB: 5UEN) in 2017. Numerous experimental developments over the past two decades have opened the door for widespread GPCR structural characterization. These efforts have resulted in the determination of three-dimensional structures for over 40 individual GPCR family members. Herein we present a comprehensive list and comparative analysis of over 180 individual GPCR structures. This includes a summary of different GPCR functional states crystallized with agonists, dual agonists, partial agonists, inverse agonists, antagonists, and allosteric modulators.

Structure and Function of Peptide-Binding G Protein-Coupled Receptors

G protein-coupled receptors (GPCRs) are the largest family of cell surface receptors and are important human drug targets. Of the 826 human GPCRs, 118 of them recognize endogenous peptide or protein ligands, and 30 of the 118 are targeted by approved drug molecules, including the very high-profile class B glucagon-like peptide 1 receptor. In this review, we analyze the 21 experimentally determined three-dimensional structures of the known peptide-binding GPCRs in relation to the endogenous peptides and drug molecules that modulate their cell signaling processes. Our integrated analyses reveal that half of the marketed drugs and most of the drugs in clinical trials that interact with peptide GPCRs are small molecules with a wide range of binding modes distinct from those of large peptide ligands. As we continue to collect additional data on these receptors from orthogonal approaches, including nuclear magnetic resonance and electron microscopy, we are beginning to understand how these receptors interact with their ligands at the molecular level and how improving the pharmacology of GPCR signal transduction requires us to study these receptors using multiple biophysical techniques.

New approaches towards the understanding of integral membrane proteins: A structural perspective on G protein-coupled receptors

Three-dimensional structure determination of integral membrane proteins has advanced in unprecedented detail our understanding of mechanistic events of how ion channels, transporters, receptors, and enzymes function. This exciting progress required a tremendous amount of methods development, as exemplified here with G protein-coupled receptors (GPCRs): Optimizing the production of GPCRs in recombinant hosts; increasing the probability of crystal formation using high-affinity ligands, nanobodies, and minimal G proteins for co-crystallization, thus stabilizing receptors into one conformation; using the T4 lysozyme technology and other fusion partners to promote crystal contacts; advancing crystallization methods including the development of novel detergents, and miniaturization and automation of the lipidic cubic phase crystallization method; the concept of conformational thermostabilization of GPCRs; and developing microfocus X-ray synchrotron technologies to analyze small GPCR crystals. However, despite immense progress to explain how GPCRs function, many receptors pose intractable hurdles to structure determination at this time. Three emerging methods, serial femtosecond crystallography, micro electron diffraction, and single particle electron cryo-microscopy, hold promise to overcome current limitations in structural membrane biology.

Structural heterogeneity of the μ-opioid receptor's conformational ensemble in the apo state

G-protein coupled receptors (GPCRs) are the largest and most pharmaceutically relevant family of membrane proteins. Here, fully unbiased, enhanced sampling simulations of a constitutively active mutant (CAM) of a class A GPCR, the μ-opioid receptor (μOR), demonstrates repeated transitions between the inactive (IS) and active-like (AS-L) states. The interconversion features typical activation/inactivation patterns involving established conformational rearrangements of conserved residues. By contrast, wild-type μOR remains in IS during the same course of simulation, consistent with the low basal activity of the protein. The simulations point to an important role of residue W293^6.48 at the "toggle switch" in the mutation-induced constitutive activation. Such role has been already observed for other CAMs of class A GPCRs. We also find a significantly populated intermediate state, rather similar to IS. Based on the remarkable accord between simulations and experiments, we suggest here that this state, which has escaped so far experimental characterization, might constitute an early step in the activation process of the apo μOR CAM.

Gonadotropin-Releasing Hormone (GnRH) Receptor Structure and GnRH Binding

Gonadotropin-releasing hormone (GnRH) regulates reproduction. The human GnRH receptor lacks a cytoplasmic carboxy-terminal tail but has amino acid sequence motifs characteristic of rhodopsin-like, class A, G protein-coupled receptors (GPCRs). This review will consider how recent descriptions of X-ray crystallographic structures of GPCRs in inactive and active conformations may contribute to understanding GnRH receptor structure, mechanism of activation and ligand binding. The structures confirmed that ligands bind to variable extracellular surfaces, whereas the seven membrane-spanning α-helices convey the activation signal to the cytoplasmic receptor surface, which binds and activates heterotrimeric G proteins. Forty non-covalent interactions that bridge topologically equivalent residues in different transmembrane (TM) helices are conserved in class A GPCR structures, regardless of activation state. Conformation-independent interhelical contacts account for a conserved receptor protein structure and their importance in the GnRH receptor structure is supported by decreased expression of receptors with mutations of residues in the network. Many of the GnRH receptor mutations associated with congenital hypogonadotropic hypogonadism, including the Glu^2.53(90) Lys mutation, involve amino acids that constitute the conserved network. Half of the ~250 intramolecular interactions in GPCRs differ between inactive and active structures. Conformation-specific interhelical contacts depend on amino acids changing partners during activation. Conserved inactive conformation-specific contacts prevent receptor activation by stabilizing proximity of TM helices 3 and 6 and a closed G protein-binding site. Mutations of GnRH receptor residues involved in these interactions, such as Arg^3.50(139) of the DRY/S motif or Tyr^7.53(323) of the N/DPxxY motif, increase or decrease receptor expression and efficiency of receptor coupling to G protein signaling, consistent with the native residues stabilizing the inactive GnRH receptor structure. Active conformation-specific interhelical contacts stabilize an open G protein-binding site. Progress in defining the GnRH-binding site has recently slowed, with evidence that Tyr^6.58(290) contacts Tyr⁵ of GnRH, whereas other residues affect recognition of Trp³ and Gly¹⁰NH₂. The surprisingly consistent observations that GnRH receptor mutations that disrupt GnRH binding have less effect on "conformationally constrained" GnRH peptides may now be explained by crystal structures of agonist-bound peptide receptors. Analysis of GPCR structures provides insight into GnRH receptor function.