The dynamics of the G protein-coupled neuropeptide Y2 receptor in monounsaturated membranes investigated by solid-state NMR spectroscopy

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

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

Multiplexing experiments in NMR and multi-nuclear MRI

Multiplexing NMR experiments by direct detection of multiple free induction decays (FIDs) in a single experiment offers a dramatic increase in the spectral information content and often yields significant improvement in sensitivity per unit time. Experiments with multi-FID detection have been designed with both homonuclear and multinuclear acquisition, and the advent of multiple receivers on commercial spectrometers opens up new possibilities for recording spectra from different nuclear species in parallel. Here we provide an extensive overview of such techniques, designed for applications in liquid- and solid-state NMR as well as in hyperpolarized samples. A brief overview of multinuclear MRI is also provided, to stimulate cross fertilization of ideas between the two areas of research (NMR and MRI). It is shown how such techniques enable the design of experiments that allow structure elucidation of small molecules from a single measurement. Likewise, in biomolecular NMR experiments multi-FID detection allows complete resonance assignment in proteins. Probes with multiple RF microcoils routed to multiple NMR receivers provide an alternative way of increasing the throughput of modern NMR systems, effectively reducing the cost of NMR analysis and increasing the information content at the same time. Solid-state NMR experiments have also benefited immensely from both parallel and sequential multi-FID detection in a variety of multi-dimensional pulse schemes. We are confident that multi-FID detection will become an essential component of future NMR methodologies, effectively increasing the sensitivity and information content of NMR measurements.

The Conformational Equilibrium of the Neuropeptide Y2 Receptor in Bilayer Membranes

Dynamic structural transitions within the seven-transmembrane bundle represent the mechanism by which G-protein-coupled receptors convert an extracellular chemical signal into an intracellular biological function. Here, the conformational dynamics of the neuropeptide Y receptor type 2 (Y2R) during activation was investigated. The apo, full agonist-, and arrestin-bound states of Y2R were prepared by cell-free expression, functional refolding, and reconstitution into lipid membranes. To study conformational transitions between these states, all six tryptophans of Y2R were 13 C-labeled. NMR-signal assignment was achieved by dynamic-nuclear-polarization enhancement and the individual functional states of the receptor were characterized by monitoring 13 C NMR chemical shifts. Activation of Y2R is mediated by molecular switches involving the toggle switch residue Trp2816.48 of the highly conserved SWLP motif and Trp3277.55 adjacent to the NPxxY motif. Furthermore, a conformationally preserved "cysteine lock"-Trp11623.50 was identified.

The Dynamics of the Neuropeptide Y Receptor Type 1 Investigated by Solid-State NMR and Molecular Dynamics Simulation

We report data on the structural dynamics of the neuropeptide Y (NPY) G-protein-coupled receptor (GPCR) type 1 (Y1R), a typical representative of class A peptide ligand GPCRs, using a combination of solid-state NMR and molecular dynamics (MD) simulation. First, the equilibrium dynamics of Y1R were studied using 15N-NMR and quantitative determination of 1H-13C order parameters through the measurement of dipolar couplings in separated-local-field NMR experiments. Order parameters reporting the amplitudes of the molecular motions of the C-H bond vectors of Y1R in DMPC membranes are 0.57 for the Cα sites and lower in the side chains (0.37 for the CH2 and 0.18 for the CH3 groups). Different NMR excitation schemes identify relatively rigid and also dynamic segments of the molecule. In monounsaturated membranes composed of longer lipid chains, Y1R is more rigid, attributed to a higher hydrophobic thickness of the lipid membrane. The presence of an antagonist or NPY has little influence on the amplitude of motions, whereas the addition of agonist and arrestin led to a pronounced rigidization. To investigate Y1R dynamics with site resolution, we conducted extensive all-atom MD simulations of the apo and antagonist-bound state. In each state, three replicas with a length of 20 μs (with one exception, where the trajectory length was 10 μs) were conducted. In these simulations, order parameters of each residue were determined and showed high values in the transmembrane helices, whereas the loops and termini exhibit much lower order. The extracellular helix segments undergo larger amplitude motions than their intracellular counterparts, whereas the opposite is observed for the loops, Helix 8, and termini. Only minor differences in order were observed between the apo and antagonist-bound state, whereas the time scale of the motions is shorter for the apo state. Although these relatively fast motions occurring with correlation times of ns up to a few µs have no direct relevance for receptor activation, it is believed that they represent the prerequisite for larger conformational transitions in proteins.

Integration of Cell-Free Expression and Solid-State NMR to Investigate the Dynamic Properties of Different Sites of the Growth Hormone Secretagogue Receptor

Cell-free expression represents an attractive method to produce large quantities of selectively labeled protein for NMR applications. Here, cell-free expression was used to label specific regions of the growth hormone secretagogue receptor (GHSR) with NMR-active isotopes. The GHSR is a member of the class A family of G protein-coupled receptors. A cell-free expression system was established to produce the GHSR in the precipitated form. The solubilized receptor was refolded in vitro and reconstituted into DMPC lipid membranes. Methionines, arginines, and histidines were chosen for ¹³C-labeling as they are representative for the transmembrane domains, the loops and flanking regions of the transmembrane α-helices, and the C-terminus of the receptor, respectively. The dynamics of the isotopically labeled residues was characterized by solid-state NMR measuring motionally averaged ¹H-¹³C dipolar couplings, which were converted into molecular order parameters. Separated local field DIPSHIFT experiments under magic-angle spinning conditions using either varying cross polarization contact times or direct excitation provided order parameters for these residues showing that the C-terminus was the segment with the highest motional amplitude. The loop regions and helix ends as well as the transmembrane regions of the GHSR represent relatively rigid segments in the overall very flexible receptor molecule. Although no site resolution could be achieved in the experiments, the previously reported highly dynamic character of the receptor concluded from uniformly ¹³C labeled receptor samples could be further specified by this segmental labeling approach, leading to a more diversified understanding of the receptor dynamics under equilibrium conditions.

Serotonin Alters the Phase Equilibrium of a Ternary Mixture of Phospholipids and Cholesterol

Unsaturated and saturated phospholipids tend to laterally segregate, especially in the presence of cholesterol. Small molecules such as neurotransmitters, toxins, drugs etc. possibly modulate this lateral segregation. The small aromatic neurotransmitter serotonin (5-HT) has been found to bind to membranes. We studied the lipid structure and packing of a ternary membrane mixture consisting of palmitoyl-oleoyl-phosphatidylcholine, palmitoyl-sphingomyelin, and cholesterol at a molar ratio of 4/4/2 in the absence and in the presence of 5-HT, using a combination of solid-state 2H NMR, atomic force microscopy, and atomistic molecular dynamics (MD) simulations. Both NMR and MD report formation of a liquid ordered (L o ) and a liquid disordered (L d ) phase coexistence with small domains. Lipid exchange between the domains was fast such that single component 2H NMR spectra are detected over a wide temperature range. A drastic restructuring of the domains was induced when 5-HT is added to the membranes at a 9 mol% concentration relative to the lipids. 2H NMR spectra of all components of the mixture showed two prominent contributions indicative of molecules of the same kind residing both in the disordered and the ordered phase. Compared to the data in the absence of 5-HT, the lipid chain order in the disordered phase was further decreased in the presence of 5-HT. Likewise, addition of serotonin increased lipid chain order within the ordered phase. These characteristic lipid chain order changes were confirmed by MD simulations. The 5-HT-induced larger difference in lipid chain order results in more pronounced differences in the hydrophobic thickness of the individual membrane domains. The correspondingly enlarged hydrophobic mismatch between ordered and disordered phases is assumed to increase the line tension at the domain boundary, which drives the system into formation of larger size domains. These results not only demonstrate that small membrane binding molecules such as neurotransmitters have a profound impact on essential membrane properties. It also suggests a mechanism by which the interaction of small molecules with membranes can influence the function of membrane proteins and non-cognate receptors. Altered membrane properties may modify lateral sorting of membrane protein, membrane protein conformation, and thus influence their function as suspected for neurotransmitters, local anesthetics, and other small drug molecules.

Capturing Peptide-GPCR Interactions and Their Dynamics

Many biological functions of peptides are mediated through G protein-coupled receptors (GPCRs). Upon ligand binding, GPCRs undergo conformational changes that facilitate the binding and activation of multiple effectors. GPCRs regulate nearly all physiological processes and are a favorite pharmacological target. In particular, drugs are sought after that elicit the recruitment of selected effectors only (biased ligands). Understanding how ligands bind to GPCRs and which conformational changes they induce is a fundamental step toward the development of more efficient and specific drugs. Moreover, it is emerging that the dynamic of the ligand–receptor interaction contributes to the specificity of both ligand recognition and effector recruitment, an aspect that is missing in structural snapshots from crystallography. We describe here biochemical and biophysical techniques to address ligand–receptor interactions in their structural and dynamic aspects, which include mutagenesis, crosslinking, spectroscopic techniques, and mass-spectrometry profiling. With a main focus on peptide receptors, we present methods to unveil the ligand–receptor contact interface and methods that address conformational changes both in the ligand and the GPCR. The presented studies highlight a wide structural heterogeneity among peptide receptors, reveal distinct structural changes occurring during ligand binding and a surprisingly high dynamics of the ligand–GPCR complexes.

Illuminating the Energy Landscape of GPCRs: The Key Contribution of Solution-State NMR Associated with Escherichia coli as an Expression Host

Conformational dynamics of GPCRs are central to their function but are difficult to explore at the atomic scale. Solution-state NMR has provided the major contribution in that area of study during the past decade, despite nonoptimized labeling schemes due to the use of insect cells and, to a lesser extent, yeast as the main expression hosts. Indeed, the most efficient isotope-labeling scheme ever to address energy landscape issues for large proteins or protein complexes relies on the use of ¹³CH₃ probes immersed in a perdeuterated dipolar environment, which is essentially out of reach of eukaryotic expression systems. In contrast, although its contribution has been underestimated because of technical issues, Escherichia coli is by far the best-adapted host for such labeling. As it is now tightly controlled, we show in this review that bacterial expression can provide an NMR spectral resolution never achieved in the GPCR field.

Helix formation and stability in membranes

In this article we review current understanding of basic principles for the folding of membrane proteins, focusing on the more abundant alpha-helical class. Membrane proteins, vital to many biological functions and implicated in numerous diseases, fold into their active conformations in the complex environment of the cell bilayer membrane. While many membrane proteins rely on the translocon and chaperone proteins to fold correctly, others can achieve their functional form in the absence of any translation apparatus or other aides. Nevertheless, the spontaneous folding process is not well understood at the molecular level. Recent findings suggest that helix fraying and loop formation may be important for overall structure, dynamics and regulation of function. Several types of membrane helices with ionizable amino acids change their topology with pH. Additionally we note that some peptides, including many that are rich in arginine, and a particular analogue of gramicidin, are able passively to translocate across cell membranes. The findings indicate that a final protein structure in a lipid-bilayer membrane is sequence-based, with lipids contributing to stability and regulation. While much progress has been made toward understanding the folding process for alpha-helical membrane proteins, it remains a work in progress. This article is part of a Special Issue entitled: Emergence of Complex Behavior in Biomembranes edited by Marjorie Longo.

Improved in Vitro Folding of the Y2 G Protein-Coupled Receptor into Bicelles

Prerequisite for structural studies on G protein-coupled receptors is the preparation of highly concentrated, stable, and biologically active receptor samples in milligram amounts of protein. Here, we present an improved protocol for Escherichia coli expression, functional refolding, and reconstitution into bicelles of the human neuropeptide Y receptor type 2 (Y₂R) for solution and solid-state NMR experiments. The isotopically labeled receptor is expressed in inclusion bodies and purified using SDS. We studied the details of an improved preparation protocol including the in vitro folding of the receptor, e.g., the native disulfide bridge formation, the exchange of the denaturating detergent SDS, and the functional reconstitution into bicelle environments of varying size. Full pharmacological functionality of the Y₂R preparation was shown by a ligand affinity of 4 nM and G-protein activation. Further, simple NMR experiments are used to test sample quality in high micromolar concentration.

Expression, Functional Characterization, and Solid-State NMR Investigation of the G Protein-Coupled GHS Receptor in Bilayer Membranes

The expression, functional reconstitution and first NMR characterization of the human growth hormone secretagogue (GHS) receptor reconstituted into either DMPC or POPC membranes is described. The receptor was expressed in E. coli. refolded, and reconstituted into bilayer membranes. The molecule was characterized by ¹⁵N and ¹³C solid-state NMR spectroscopy in the absence and in the presence of its natural agonist ghrelin or an inverse agonist. Static ¹⁵N NMR spectra of the uniformly labeled receptor are indicative of axially symmetric rotational diffusion of the G protein-coupled receptor in the membrane. In addition, about 25% of the ¹⁵N sites undergo large amplitude motions giving rise to very narrow spectral components. For an initial quantitative assessment of the receptor mobility, ¹H-¹³C dipolar coupling values, which are scaled by molecular motions, were determined quantitatively. From these values, average order parameters, reporting the motional amplitudes of the individual receptor segments can be derived. Average backbone order parameters were determined with values between 0.56 and 0.69, corresponding to average motional amplitudes of 40–50° of these segments. Differences between the receptor dynamics in DMPC or POPC membranes were within experimental error. Furthermore, agonist or inverse agonist binding only insignificantly influenced the average molecular dynamics of the receptor.

Structure and dynamics of the aliphatic cholesterol side chain in membranes as studied by (2)H NMR spectroscopy and molecular dynamics simulation

Cholesterol is an evolutionarily highly optimized molecule particularly known for its ability to condense the phospholipids in cellular membranes. Until recently, the accompanying increase in the chain order of the surrounding phospholipids was attributed to the planar and rigid tetracyclic ring structure of cholesterol. However, detailed investigations of cholesterol's aliphatic side chain demonstrated that this side chain is responsible for approximately half of the condensation effect. Therefore, we investigated the structure and dynamics of the aliphatic side chain of cholesterol using (2)H solid-state nuclear magnetic resonance (NMR) spectroscopy and microsecond timescale all-atom molecular dynamics (MD) simulations in four different model membranes: POPC, DPPC, PSM, and POPC/PSM (1 : 1 mol/mol) and at three different temperatures: 5 °C, 37 °C, and 50 °C. A cholesterol variant, in which 11 hydrogens of the aliphatic side chain were exchanged for deuterium, was used and the respective (2)H NMR spectra confirmed the axially asymmetric rotational diffusion of cholesterol in DPPC and PSM. Furthermore, NMR spectra indicated that some hydrogens showed an unexpected magnetic inequivalency. This finding was confirmed by all-atom molecular dynamics simulations and detailed analysis revealed that the hydrogens of the methylene groups at C22, C23, and C24 are magnetically inequivalent. This inequivalency is caused by steric clashes of the aliphatic side chain with the ring structure of cholesterol as well as the branched C21 methyl group. These excluded volume effects result in reduced conformational flexibility of the aliphatic side chain of cholesterol and explain its high order (order parameter of 0.78 for chain motions) and large contribution to the condensation effect. Additionally, the motional pattern of the side chain becomes highly anisotropic such that it shows larger fluctuations perpendicular to the ring plane of cholesterol with a biaxiality of the distribution of 0.046. Overall, our results shed light on the mechanism how the aliphatic side chain is able to contribute about half of the condensation effect of cholesterol.

Studying Dynamics by Magic-Angle Spinning Solid-State NMR Spectroscopy: Principles and Applications to Biomolecules

Magic-angle spinning solid-state NMR spectroscopy is an important technique to study molecular structure, dynamics and interactions, and is rapidly gaining importance in biomolecular sciences. Here we provide an overview of experimental approaches to study molecular dynamics by MAS solid-state NMR, with an emphasis on the underlying theoretical concepts and differences of MAS solid-state NMR compared to solution-state NMR. The theoretical foundations of nuclear spin relaxation are revisited, focusing on the particularities of spin relaxation in solid samples under magic-angle spinning. We discuss the range of validity of Redfield theory, as well as the inherent multi-exponential behavior of relaxation in solids. Experimental challenges for measuring relaxation parameters in MAS solid-state NMR and a few recently proposed relaxation approaches are discussed, which provide information about time scales and amplitudes of motions ranging from picoseconds to milliseconds. We also discuss the theoretical basis and experimental measurements of anisotropic interactions (chemical-shift anisotropies, dipolar and quadrupolar couplings), which give direct information about the amplitude of motions. The potential of combining relaxation data with such measurements of dynamically-averaged anisotropic interactions is discussed. Although the focus of this review is on the theoretical foundations of dynamics studies rather than their application, we close by discussing a small number of recent dynamics studies, where the dynamic properties of proteins in crystals are compared to those in solution.

Juxta-terminal Helix Unwinding as a Stabilizing Factor to Modulate the Dynamics of Transmembrane Helices

Transmembrane helices of integral membrane proteins often are flanked by interfacial aromatic residues that can serve as anchors to aid the stabilization of a tilted transmembrane orientation. Yet, physical factors that govern the orientation or dynamic averaging of individual transmembrane helices are not well understood and have not been adequately explained. Using solid-state (2) H NMR spectroscopy to examine lipid bilayer-incorporated model peptides of the GWALP23 (acetyl-GGALW(LA)6 LWLAGA-amide) family, we observed substantial unwinding at the terminals of several tilted helices spanning the membranes of DLPC, DMPC, or DOPC lipid bilayers. The fraying of helix ends might be vital for defining the dynamics and orientations of transmembrane helices in lipid bilayer membranes.

Identification of a Conformational Equilibrium That Determines the Efficacy and Functional Selectivity of the μ-Opioid Receptor

G-protein-coupled receptor (GPCR) ligands impart differing degrees of signaling in the G-protein and arrestin pathways, in phenomena called "biased signaling". However, the mechanism underlying the biased signaling of GPCRs is still unclear, although crystal structures of GPCRs bound to the G protein or arrestin are available. In this study, we observed the NMR signals from methionine residues of the μ-opioid receptor (μOR) in the balanced- and biased-ligand-bound states. We found that the intracellular cavity of μOR exists in an equilibrium between closed and multiple open conformations with coupled conformational changes on the transmembrane helices 3, 5, 6, and 7, and that the population of each open conformation determines the G-protein- and arrestin-mediated signaling levels in each ligand-bound state. These findings provide insight into the biased signaling of GPCRs and will be helpful for development of analgesics that stimulate μOR with reduced tolerance and dependence.