Ligands selectively tune the local and global motions of neurotensin receptor 1 (NTS(1))

Advances in uncovering the mechanisms of macromolecular conformational entropy

During protein folding, proteins transition from a disordered polymer into a globular structure, markedly decreasing their conformational degrees of freedom, leading to a substantial reduction in entropy. Nonetheless, folded proteins retain substantial entropy as they fluctuate between the conformations that make up their native state. This residual entropy contributes to crucial functions like binding and catalysis, supported by growing evidence primarily from NMR and simulation studies. Here, we propose three major ways that macromolecules use conformational entropy to perform their functions; first, prepaying entropic cost through ordering of the ground state; second, spatially redistributing entropy, in which a decrease in entropy in one area is reciprocated by an increase in entropy elsewhere; third, populating catalytically competent ensembles, in which conformational entropy within the enzymatic scaffold aids in lowering transition state barriers. We also provide our perspective on how solving the current challenge of structurally defining the ensembles encoding conformational entropy will lead to new possibilities for controlling binding, catalysis and allostery.

CXC Chemokine Ligand 12 Facilitates Gi Protein Binding to CXC Chemokine Receptor 4 by Stabilizing Packing of the Proline-Isoleucine-Phenylalanine Motif: Insights from Automated Path Searching

The activation of G protein-coupled receptors (GPCRs) is a complex multibody multievent process involving agonist binding, receptor activation, G protein coupling, and subsequent G protein activation. The order and energetics of these events, though crucial for the rational design of selective GPCR drugs, are challenging to characterize and remain largely underexplored. Here, we employed molecular dynamics simulations and our recently developed traveling salesman-based automated path searching (TAPS) algorithm to efficiently locate the minimum free-energy paths for the coupling events of the CXC chemokine receptor 4 (CXCR4) with its endogenous ligand CXC chemokine ligand 12 (CXCL12) and Gi protein. We show that, after overcoming three low energy barriers (3.24-6.89 kcal/mol), Gi alone can precouple with CXCR4 even without CXCL12, consistent with previous reports on the existence of the apo CXCR4-Gi complex and our NanoBiT experiments. The highest barrier of 6.89 kcal/mol in this process corresponds to the packing of the proline-isoleucine-phenylalanine (PIF) motif of CXCR4. Interestingly, without Gi, CXCL12 alone cannot activate CXCR4 (high barrier of 18.89 kcal/mol). Instead, it can enhance Gi coupling by circumventing the energy barrier of PIF packing. Shedding new light on the activation mechanism of CXCR4, these results present TAPS as a promising tool for uncovering complete activation pathways of GPCRs and the corresponding agonist design.

Deciphering the allosteric dynamics of GPR120-fatty acid interactions within a bilayer nanogold electrochemical receptor biosensor: the impact of replacing tryptophan 198 with proline

GPR120 is a free fatty acid receptor capable of signalling excess fatty acids. GPR120 can be activated by various types of free fatty acids, causing intracellular signal transduction and exerting energy regulation, immune homeostasis, and neuronal functions. It has been suggested that Trp198 may be an important residue in the recognition and activation of GPR120 by fatty acid ligands, but direct experimental evidence is lacking. In this study, a GPR120-based bilayer gold nanoparticle biosensor (Trp198→Pro) was constructed by genetically manipulating Trp198 on GPR120 by replacing it with proline for the determination of linkage variability between 14 naturally occurring fatty acid ligands and mutant receptors. The results showed that both before and after amino acid substitution the GPR120 bilayer nanogold receptor sensor responded to all 14 natural fatty acid ligands. And the linkage transformation constants of crotonic acid, dodecanoic acid, oleic acid, linoleic acid, α-linolenic acid, and DHA decreased after Trp198 was replaced by Pro. To further reveal its molecular recognition mechanism, molecular simulation docking experiments were performed on GPR120 and 14 fatty acid ligand compounds before and after amino acid substitutions, respectively. The results showed that before and after the amino acid substitutions, the binding conformational affinity values of GPR120 docked with the ligands were negative, implying that these fatty acid ligands can spontaneously bind to the active pocket of GPR120 without absorbing external energy. Upon replacement of Trp198 by Pro, the active pocket of GPR120 at the optimal docking site with the fatty acid ligand is altered, leading to changes in the amino acid residues that exert the interaction. The above results demonstrate that Trp198 indeed plays an important role in the recognition of fatty acid ligands on GPR120. The present study provides direct quantitative evidence for the roles played by different amino acid residues in receptor-ligand recognition and interaction. At the same time, it provides new ideas for the study of other receptor-ligand-linked metastable mechanisms and kinetic laws.

Large scale investigation of GPCR molecular dynamics data uncovers allosteric sites and lateral gateways

G protein-coupled receptors (GPCRs) constitute a functionally diverse protein family and are targets for a broad spectrum of pharmaceuticals. Technological progress in X-ray crystallography and cryogenic electron microscopy has enabled extensive, high-resolution structural characterisation of GPCRs in different conformational states. However, as highly dynamic events underlie GPCR signalling, a complete understanding of GPCR functionality requires insights into their conformational dynamics. Here, we present a large dataset of molecular dynamics simulations covering 60% of currently available GPCR structures. Our analysis reveals extensive local "breathing" motions of the receptor on a nano- to microsecond timescale and provides access to numerous previously unexplored receptor conformational states. Furthermore, we reveal that receptor flexibility impacts the shape of allosteric drug binding sites, which frequently adopt partially or completely closed states in the absence of a molecular modulator. We demonstrate that exploring membrane lipid dynamics and their interaction with GPCRs is an efficient approach to expose such hidden allosteric sites and even lateral ligand entrance gateways. The obtained insights and generated dataset on conformations, allosteric sites and lateral entrance gates in GPCRs allows us to better understand the functionality of these receptors and opens new therapeutic avenues for drug-targeting strategies.

Biased activation of the vasopressin V2 receptor probed by molecular dynamics simulations, NMR and pharmacological studies

G protein-coupled receptors (GPCRs) control critical cell signaling. Their response to extracellular stimuli involves conformational changes to convey signals to intracellular effectors, among which the most important are G proteins and β-arrestins (βArrs). Biased activation of one pathway is a field of intense research in GPCR pharmacology. Combining NMR, site-directed mutagenesis, molecular pharmacology, and molecular dynamics (MD) simulations, we studied the conformational diversity of the vasopressin V2 receptor (V2R) bound to different types of ligands: the antagonist Tolvaptan, the endogenous unbiased agonist arginine-vasopressin, and MCF14, a partial Gs protein-biased agonist. A double-labeling NMR scheme was developed to study the receptor conformational changes and ligand binding: V2R was subjected to lysine 13CH3 methylation for complementary NMR studies, whereas the agonists were tagged with a paramagnetic probe. Paramagnetic relaxation enhancements and site-directed mutagenesis validated the ligand binding modes in the MD simulations. We found that the bias for the Gs protein over the βArr pathway involves interactions between the conserved NPxxY motif in the transmembrane helix 7 (TM7) and TM3, compacting helix 8 (H8) toward TM1 and likely inhibiting βArr signaling. A similar mechanism was elicited for the pathogenic mutation I130N, which constitutively activates the Gs proteins without concomitant βArr recruitment. The findings suggest common patterns of biased signaling in class A GPCRs, as well as a rationale for the design of G protein-biased V2R agonists.

A Fast Method to Monitor Tyrosine Kinase Inhibitor Mechanisms

Methionine residues within the kinase domain of Src serve as unique NMR probes capable of distinguishing between distinct conformational states of full-length Src, including alternative drug-inhibited forms. This approach offers a rapid method to differentiate between various inhibition mechanisms at any stage of drug development, eliminating the need to resolve the structure of Src-drug complexes. Using selectively 13C-methyl-enriched methionine, spectra can be acquired in under an hour, while natural abundance spectra with comparable information are achievable within a few hours.

Revealing the graded activation mechanism of neurotensin receptor 1

Graded activation contributes to the precise regulation of GPCR activity, presenting new opportunities for drug design. In this work, a total of 10 μs enhanced-sampling simulations are performed to provide molecular insights into the binding dynamics differences of the neurotensin receptor 1 (NTSR1) to the full agonist SRI-9829, partial agonist RTI-3a and inverse agonist SR48692. The possible graded activation mechanism of NTSR1 is revealed by an integrated analysis utilizing the reweighted potential of mean force (PMF), deep learning (DL) and transfer entropy (TE). Specifically, the orthosteric pocket is observed to undergo expansion and contraction, with the G-protein-binding site experiencing interconversions among the inactive, intermediate and active-like states. Detailed structural comparisons capture subtle conformational differences arising from ligand binding in allosteric signaling, which can well explain the graded activation. Critical microswitches that contribute to graded activation are efficiently identified with the DL model. TE calculations enable the visualization of allosteric communication networks within the receptor, elucidating the driver-responder relationships associated with signal transduction. Fortunately, the dissociation of the full agonist from the orthosteric pocket is observed. The current findings systematically reveal the mechanism of NTSR1 graded activation, and also provide implications for structure-based drug design.

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.

Persistence of Methionine Side Chain Mobility at Low Temperatures in a Nine-Residue Low Complexity Peptide, as Probed by 2 H Solid-State NMR

Methionine side chains are flexible entities which play important roles in defining hydrophobic interfaces. We utilize deuterium static solid-state NMR to assess rotameric inter-conversions and other dynamic modes of the methionine in the context of a nine-residue random-coil peptide (RC9) with the low-complexity sequence GGKGMGFGL. The measurements in the temperature range of 313 to 161 K demonstrate that the rotameric interconversions in the hydrated solid powder state persist to temperatures below 200 K. Removal of solvation significantly reduces the rate of the rotameric motions. We employed 2 H NMR line shape analysis, longitudinal and rotation frame relaxation, and chemical exchange saturation transfer methods and found that the combination of multiple techniques creates a significantly more refined model in comparison with a single technique. Further, we compare the most essential features of the dynamics in RC9 to two different methionine-containing systems, characterized previously. Namely, the M35 of hydrated amyloid-β1-40 in the three-fold symmetric polymorph as well as Fluorenylmethyloxycarbonyl (FMOC)-methionine amino acid with the bulky hydrophobic group. The comparison suggests that the driving force for the enhanced methionine side chain mobility in RC9 is the thermodynamic factor stemming from distributions of rotameric populations, rather than the increase in the rate constant.

Unravelling the mechanism of neurotensin recognition by neurotensin receptor 1

The conformational ensembles of G protein-coupled receptors (GPCRs) include inactive and active states. Spectroscopy techniques, including NMR, show that agonists, antagonists and other ligands shift the ensemble toward specific states depending on the pharmacological efficacy of the ligand. How receptors recognize ligands and the kinetic mechanism underlying this population shift is poorly understood. Here, we investigate the kinetic mechanism of neurotensin recognition by neurotensin receptor 1 (NTS1) using 19F-NMR, hydrogen-deuterium exchange mass spectrometry and stopped-flow fluorescence spectroscopy. Our results indicate slow-exchanging conformational heterogeneity on the extracellular surface of ligand-bound NTS1. Numerical analysis of the kinetic data of neurotensin binding to NTS1 shows that ligand recognition follows an induced-fit mechanism, in which conformational changes occur after neurotensin binding. This approach is applicable to other GPCRs to provide insight into the kinetic regulation of ligand recognition by GPCRs.

Interplay of thermodynamics and evolution within the ternary ligand-GPCR-G protein complex

Recent studies on G protein-coupled receptors (GPCRs) dynamics report that GPCRs adopt a wide range of conformations that coexist in equilibrium, with the apo state of a GPCR having a high entropy. The formation of a ligand-GPCR-transducer complex comes with a reduction of conformational space and therefore with an entropic cost. We hypothesize that the availability of binding partners, their binding affinity and the rigidity of the respective binding sites are reflected in a distinct degree of sequence conservation to balance the energetic cost of intra- and extracellular binding events. Here, we outline the current findings in delineating the conformational space and include sequential conservation of many-to-many ligand-receptor systems to discuss the entropic cost that comes with GPCR signal transduction.

Stabilization of pre-existing neurotensin receptor conformational states by β-arrestin-1 and the biased allosteric modulator ML314

The neurotensin receptor 1 (NTS1) is a G protein-coupled receptor (GPCR) with promise as a drug target for the treatment of pain, schizophrenia, obesity, addiction, and various cancers. A detailed picture of the NTS1 structural landscape has been established by X-ray crystallography and cryo-EM and yet, the molecular determinants for why a receptor couples to G protein versus arrestin transducers remain poorly defined. We used 13CεH3-methionine NMR spectroscopy to show that binding of phosphatidylinositol-4,5-bisphosphate (PIP2) to the receptor's intracellular surface allosterically tunes the timescale of motions at the orthosteric pocket and conserved activation motifs - without dramatically altering the structural ensemble. β-arrestin-1 further remodels the receptor ensemble by reducing conformational exchange kinetics for a subset of resonances, whereas G protein coupling has little to no effect on exchange rates. A β-arrestin biased allosteric modulator transforms the NTS1:G protein complex into a concatenation of substates, without triggering transducer dissociation, suggesting that it may function by stabilizing signaling incompetent G protein conformations such as the non-canonical state. Together, our work demonstrates the importance of kinetic information to a complete picture of the GPCR activation landscape.

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.

Neurotensin and Alcohol Use Disorders: Towards a Pharmacological Treatment

Harmful alcohol use is responsible for a group of disorders collectively named alcohol use disorders (AUDs), according to the DSM-5 classification. The damage induced by alcohol depends on the amount, time, and consumption patterns (continuous and heavy episodic drinking). It affects individual global well-being and social and familial environments with variable impact. Alcohol addiction manifests with different degrees of organ and mental health detriment for the individual, exhibiting two main traits: compulsive drinking and negative emotional states occurring at withdrawal, frequently causing relapse episodes. Numerous individual and living conditions, including the concomitant use of other psychoactive substances, lie in the complexity of AUD. Ethanol and its metabolites directly impact the tissues and may cause local damage or alter the homeostasis of brain neurotransmission, immunity scaffolding, or cell repair biochemical pathways. Brain modulator and neurotransmitter-assembled neurocircuitries govern reward, reinforcement, social interaction, and consumption of alcohol behaviors in an intertwined manner. Experimental evidence supports the participation of neurotensin (NT) in preclinical models of alcohol addiction. For example, NT neurons in the central nucleus of the amygdala projecting to the parabrachial nucleus strengthen alcohol consumption and preference. In addition, the levels of NT in the frontal cortex were found to be lower in rats bred to prefer alcohol to water in a free alcohol-water choice compared to wild-type animals. NT receptors 1 and 2 seem to be involved in alcohol consumption and alcohol effects in several models of knockout mice. This review aims to present an updated picture of the role of NT systems in alcohol addiction and the possible use of nonpeptide ligands modulating the activity of the NT system, applied to experimental animal models of harmful drinking behavior mimicking alcohol addiction leading to health ruin in humans.