Sodium is a negative allosteric regulator of the ghrelin receptor.

The functional properties of G protein-coupled receptors (GPCRs) are intimately associated with the different components in their cellular environment. Among them, sodium ions have been proposed to play a substantial role as endogenous allosteric modulators of GPCR-mediated signaling. However, this sodium effect and the underlying mechanisms are still unclear for most GPCRs. Here, we identified sodium as a negative allosteric modulator of the ghrelin receptor GHSR (growth hormone secretagogue receptor). Combining 23Na-nuclear magnetic resonance (NMR), molecular dynamics, and mutagenesis, we provide evidence that, in GHSR, sodium binds to the allosteric site conserved in class A GPCRs. We further leveraged spectroscopic and functional assays to show that sodium binding shifts the conformational equilibrium toward the GHSR-inactive ensemble, thereby decreasing basal and agonist-induced receptor-catalyzed G protein activation. All together, these data point to sodium as an allosteric modulator of GHSR, making this ion an integral component of the ghrelin signaling machinery.

MDexciteR: Enhanced Sampling Molecular Dynamics by Excited Normal Modes or Principal Components Obtained from Experiments.

Molecular dynamics with excited normal modes (MDeNM) is an enhanced sampling method for exploring conformational changes in proteins with minimal biases. The excitation corresponds to injecting kinetic energy along normal modes describing intrinsic collective motions. Herein, we developed a new automated open-source implementation, MDexciteR (https://github.com/mcosta27/MDexciteR), enabling the integration of MDeNM with two commonly used simulation programs with GPU support. Second, we generalized the method to include the excitation of principal components calculated from experimental ensembles. Finally, we evaluated whether the use of coarse-grained normal modes calculated with elastic network representations preserved the performance and accuracy of the method. The advantages and limitations of these new approaches are discussed based on results obtained for three different protein test cases: two globular and a protein/membrane system.

Concerted conformational dynamics and water movements in the ghrelin G protein-coupled receptor.

There is increasing support for water molecules playing a role in signal propagation through G protein-coupled receptors (GPCRs). However, exploration of the hydration features of GPCRs is still in its infancy. Here, we combined site-specific labeling with unnatural amino acids to molecular dynamics to delineate how local hydration of the ghrelin receptor growth hormone secretagogue receptor (GHSR) is rearranged upon activation. We found that GHSR is characterized by a specific hydration pattern that is selectively remodeled by pharmacologically distinct ligands and by the lipid environment. This process is directly related to the concerted movements of the transmembrane domains of the receptor. These results demonstrate that the conformational dynamics of GHSR are tightly coupled to the movements of internal water molecules, further enhancing our understanding of the molecular bases of GPCR-mediated signaling.

Allosteric modulation of ghrelin receptor signaling by lipids.

The membrane is an integral component of the G protein-coupled receptor signaling machinery. Here we demonstrate that lipids regulate the signaling efficacy and selectivity of the ghrelin receptor GHSR through specific interactions and bulk effects. We find that PIP2 shifts the conformational equilibrium of GHSR away from its inactive state, favoring basal and agonist-induced G protein activation. This occurs because of a preferential binding of PIP2 to specific intracellular sites in the receptor active state. Another lipid, GM3, also binds GHSR and favors G protein activation, but mostly in a ghrelin-dependent manner. Finally, we find that not only selective interactions but also the thickness of the bilayer reshapes the conformational repertoire of GHSR, with direct consequences on G protein selectivity. Taken together, this data illuminates the multifaceted role of the membrane components as allosteric modulators of how ghrelin signal could be propagated.

Cryo-electron microscopy structure of the antidiuretic hormone arginine-vasopressin V2 receptor signaling complex.

The antidiuretic hormone arginine-vasopressin (AVP) forms a signaling complex with the V2 receptor (V2R) and the Gs protein, promoting kidney water reabsorption. Molecular mechanisms underlying activation of this critical G protein-coupled receptor (GPCR) signaling system are still unknown. To fill this gap of knowledge, we report here the cryo-electron microscopy structure of the AVP-V2R-Gs complex. Single-particle analysis revealed the presence of three different states. The two best maps were combined with computational and nuclear magnetic resonance spectroscopy constraints to reconstruct two structures of the ternary complex. These structures differ in AVP and Gs binding modes. They reveal an original receptor-Gs interface in which the Gαs subunit penetrates deep into the active V2R. The structures help to explain how V2R R137H or R137L/C variants can lead to two severe genetic diseases. Our study provides important structural insights into the function of this clinically relevant GPCR signaling complex.

Structure and dynamics of G protein-coupled receptor-bound ghrelin reveal the critical role of the octanoyl chain.

Ghrelin plays a central role in controlling major biological processes. As for other G protein-coupled receptor (GPCR) peptide agonists, the structure and dynamics of ghrelin bound to its receptor remain obscure. Using a combination of solution-state NMR and molecular modeling, we demonstrate that binding to the growth hormone secretagogue receptor is accompanied by a conformational change in ghrelin that structures its central region, involving the formation of a well-defined hydrophobic core. By comparing its acylated and nonacylated forms, we conclude that the ghrelin octanoyl chain is essential to form the hydrophobic core and promote access of ghrelin to the receptor ligand-binding pocket. The combination of coarse-grained molecular dynamics studies and NMR should prove useful in improving our mechanistic understanding of the complex conformational space explored by a natural peptide agonist when binding to its GPCR. Such information should also facilitate the design of new ghrelin receptor-selective drugs.

Molecular Dynamics Simulations of the Allosteric Modulation of the Adenosine A2A Receptor by a Mini-G Protein.

Through their coupling to G proteins, G Protein-Coupled Receptors (GPCRs) trigger cellular responses to various signals. Some recent experiments have interestingly demonstrated that the G protein can also act on the receptor by favoring a closed conformation of its orthosteric site, even in the absence of a bound agonist. In this work, we explored such an allosteric modulation by performing extensive molecular dynamics simulations on the adenosine A2 receptor (A2AR) coupled to the Mini-Gs protein. In the presence of the Mini-Gs, we confirmed a restriction of the receptor's agonist binding site that can be explained by a modulation of the intrinsic network of contacts of the receptor. Of interest, we observed similar effects with the C-terminal helix of the Mini-Gs, showing that the observed effect on the binding pocket results from direct local contacts with the bound protein partner that cause a rewiring of the whole receptor's interaction network.

GHSR-D2R heteromerization modulates dopamine signaling through an effect on G protein conformation.

The growth hormone secretagogue receptor (GHSR) and dopamine receptor (D2R) have been shown to oligomerize in hypothalamic neurons with a significant effect on dopamine signaling, but the molecular processes underlying this effect are still obscure. We used here the purified GHSR and D2R to establish that these two receptors assemble in a lipid environment as a tetrameric complex composed of two each of the receptors. This complex further recruits G proteins to give rise to an assembly with only two G protein trimers bound to a receptor tetramer. We further demonstrate that receptor heteromerization directly impacts on dopamine-mediated Gi protein activation by modulating the conformation of its α-subunit. Indeed, association to the purified GHSR:D2R heteromer triggers a different active conformation of Gαi that is linked to a higher rate of GTP binding and a faster dissociation from the heteromeric receptor. This is an additional mechanism to expand the repertoire of GPCR signaling modulation that could have implications for the control of dopamine signaling in normal and physiopathological conditions.

Use of Molecular Modeling to Design Selective NTS2 Neurotensin Analogues.

Neurotensin exerts potent analgesia by acting at both NTS1 and NTS2 receptors, whereas NTS1 activation also results in other physiological effects such as hypotension and hypothermia. Here, we used molecular modeling approach to design highly selective NTS2 ligands by investigating the docking of novel NT[8-13] compounds at both NTS1 and NTS2 sites. Molecular dynamics simulations revealed an interaction of the Tyr11 residue of NT[8-13] with an acidic residue (Glu179) located in the ECL2 of hNTS2 or with a basic residue (Arg212) at the same position in hNTS1. The importance of the residue at position 11 for NTS1/NTS2 selectivity was further demonstrated by the design of new NT analogues bearing basic (Lys, Orn) or acid (Asp or Glu) function. As predicted by the molecular dynamics simulations, binding of NT[8-13] analogues harboring a Lys11 exhibited higher affinity toward the hNTS1-R212E mutant receptor, in which Arg212 was substituted by the negatively charged Glu residue.

Coarse-Grained Prediction of Peptide Binding to G-Protein Coupled Receptors.

In this study, we used the Martini Coarse-Grained model with no applied restraints to predict the binding mode of some peptides to G-Protein Coupled Receptors (GPCRs). Both the Neurotensin-1 and the chemokine CXCR4 receptors were used as test cases. Their ligands, NTS8-13 and CVX15 peptides, respectively, were initially positioned in the surrounding water box. Using a protocol based on Replica Exchange Molecular Dynamics (REMD), both opening of the receptors and entry of the peptides into their dedicated pockets were observed on the µs time-scale. After clustering, the most statistically representative orientations were closely related to the X-ray structures of reference, sharing both RMSD lower than 3 Å and most of the native contacts. These results demonstrate that such a model, that does not require access to tremendous computational facilities, can be helpful in predicting peptide binding to GPCRs as well as some of the receptor's conformational changes required for this key step. We also discuss how such an approach can now help to predict, de novo, the interactions of GPCRs with other intra- or extracellular peptide/protein partners.

Structural and Molecular Determinants of Membrane Binding by the HIV-1 Matrix Protein.

Assembly of HIV-1 particles is initiated by the trafficking of viral Gag polyproteins from the cytoplasm to the plasma membrane, where they co-localize and bud to form immature particles. Membrane targeting is mediated by the N-terminally myristoylated matrix (MA) domain of Gag and is dependent on the plasma membrane marker phosphatidylinositol-4,5-bisphosphate [PI(4,5)P2]. Recent studies revealed that PI(4,5)P2 molecules containing truncated acyl chains [tr-PI(4,5)P2] are capable of binding MA in an "extended lipid" conformation and promoting myristoyl exposure. Here we report that tr-PI(4,5)P2 molecules also readily bind to non-membrane proteins, including HIV-1 capsid, which prompted us to re-examine MA-PI(4,5)P2 interactions using native lipids and membrane mimetic liposomes and bicelles. Liposome binding trends observed using a recently developed NMR approach paralleled results of flotation assays, although the affinities measured under the equilibrium conditions of NMR experiments were significantly higher. Native PI(4,5)P2 enhanced MA binding to liposomes designed to mimic non-raft-like regions of the membrane, suggesting the possibility that binding of the protein to disordered domains may precede Gag association with, or nucleation of, rafts. Studies with bicelles revealed a subset of surface and myr-associated MA residues that are sensitive to native PI(4,5)P2, but cleft residues that interact with the 2'-acyl chains of tr-PI(4,5)P2 molecules in aqueous solution were insensitive to native PI(4,5)P2 in bicelles. Our findings call to question extended-lipid MA:membrane binding models, and instead support a model put forward from coarse-grained simulations indicating that binding is mediated predominantly by dynamic, electrostatic interactions between conserved basic residues of MA and multiple PI(4,5)P2 and phosphatidylserine molecules.

Conformational Equilibrium of CDK/Cyclin Complexes by Molecular Dynamics with Excited Normal Modes.

Cyclin-dependent kinases (CDKs) and their associated regulatory cyclins are central for timely regulation of cell-cycle progression. They constitute attractive pharmacological targets for development of anticancer therapeutics, since they are frequently deregulated in human cancers and contribute to sustained, uncontrolled tumor proliferation. Characterization of their structural/dynamic features is essential to gain in-depth insight into structure-activity relationships. In addition, the identification of druggable pockets or key intermediate conformations yields potential targets for the development of novel classes of inhibitors. Structural studies of CDK2/cyclin A have provided a wealth of information concerning monomeric/heterodimeric forms of this kinase. There is, however, much less structural information for other CDK/cyclin complexes, including CDK4/cyclin D1, which displays an alternative (open) position of the cyclin partner relative to CDK, contrasting with the closed CDK2/cyclin A conformation. In this study, we carried out normal-mode analysis and enhanced sampling simulations with our recently developed method, molecular dynamics with excited normal modes, to understand the conformational equilibrium on these complexes. Interestingly, the lowest-frequency normal mode computed for each complex described the transition between the open and closed conformations. Exploration of these motions with an explicit-solvent representation using molecular dynamics with excited normal modes confirmed that the closed conformation is the most stable for the CDK2/cyclin A complex, in agreement with their experimentally available structures. On the other hand, we clearly show that an open↔closed equilibrium may exist in CDK4/cyclin D1, with closed conformations resembling that captured for CDK2/cyclin A. Such conformational preferences may result from the distinct distributions of frustrated contacts in each complex. Using the same approach, the putative roles of the Thr(160) phosphoryl group and the T-loop conformation were investigated. These results provide a dynamic view of CDKs revealing intermediate conformations not yet characterized for CDK members other than CDK2, which will be useful for the design of inhibitors targeting critical conformational transitions.

Identification of TAX2 peptide as a new unpredicted anti-cancer agent.

The multi-modular glycoprotein thrombospondin-1 (TSP-1) is considered as a key actor within the tumor microenvironment. Besides, TSP-1 binding to CD47 is widely reported to regulate cardiovascular function as it promotes vasoconstriction and angiogenesis limitation. Therefore, many studies focused on targeting TSP-1:CD47 interaction, aiming for up-regulation of physiological angiogenesis to enhance post-ischemia recovery or to facilitate engraftment. Thus, we sought to identify an innovative selective antagonist for TSP-1:CD47 interaction. Protein-protein docking and molecular dynamics simulations were conducted to design a novel CD47-derived peptide, called TAX2. TAX2 binds TSP-1 to prevent TSP-1:CD47 interaction, as revealed by ELISA and co-immunoprecipitation experiments. Unexpectedly, TAX2 inhibits in vitro and ex vivo angiogenesis features in a TSP-1-dependent manner. Consistently, our data highlighted that TAX2 promotes TSP-1 binding to CD36-containing complexes, leading to disruption of VEGFR2 activation and downstream NO signaling. Such unpredicted results prompted us to investigate TAX2 potential in tumor pathology. A multimodal imaging approach was conducted combining histopathological staining, MVD, MRI analysis and µCT monitoring for tumor angiography longitudinal follow-up and 3D quantification. TAX2 in vivo administrations highly disturb syngeneic melanoma tumor vascularization inducing extensive tumor necrosis and strongly inhibit growth rate and vascularization of human pancreatic carcinoma xenografts in nude mice.

Ghrelin receptor conformational dynamics regulate the transition from a preassembled to an active receptor:Gq complex.

How G protein-coupled receptor conformational dynamics control G protein coupling to trigger signaling is a key but still open question. We addressed this question with a model system composed of the purified ghrelin receptor assembled into lipid discs. Combining receptor labeling through genetic incorporation of unnatural amino acids, lanthanide resonance energy transfer, and normal mode analyses, we directly demonstrate the occurrence of two distinct receptor:Gq assemblies with different geometries whose relative populations parallel the activation state of the receptor. The first of these assemblies is a preassembled complex with the receptor in its basal conformation. This complex is specific of Gq and is not observed with Gi. The second one is an active assembly in which the receptor in its active conformation triggers G protein activation. The active complex is present even in the absence of agonist, in a direct relationship with the high constitutive activity of the ghrelin receptor. These data provide direct evidence of a mechanism for ghrelin receptor-mediated Gq signaling in which transition of the receptor from an inactive to an active conformation is accompanied by a rearrangement of a preassembled receptor:G protein complex, ultimately leading to G protein activation and signaling.

Coarse-grained simulations of the HIV-1 matrix protein anchoring: revisiting its assembly on membrane domains.

In the accepted model for human immunodeficiency virus preassembly in infected host cells, the anchoring to the intracellular leaflet of the membrane of the matrix domain (MA) that lies at the N-terminus of the viral Gag protein precursor appears to be one of the crucial steps for particle assembly. In this study, we simulated the membrane anchoring of human immunodeficiency virus-1 myristoylated MA protein using a coarse-grained representation of both the protein and the membrane. Our calculations first suggest that the myristoyl group could spontaneously release from its initial hydrophobic pocket before MA protein interacts with the lipid membrane. All-atom simulations confirmed this possibility with a related energy cost estimated to be ~5 kcal.mol(-1). The phosphatidylinositol (4,5) bisphosphate (PI(4,5)P2) head binds preferentially to the MA highly basic region as described in available NMR data, but interestingly without flipping of its 2' acyl chain into the MA protein. Moreover, MA was able to confine PI(4,5)P2 lipids all around its molecular surface after having found a stable orientation at the membrane surface. Our results suggest that this orientation is dependent on Myr anchoring and that this confinement induces a lateral segregation of PI(4,5)P2 in domains. This is consistent with a PI(4,5)P2 enrichment of the virus envelope as compared to the host cell membrane.

Conformational restriction of G-proteins Coupled Receptors (GPCRs) upon complexation to G-proteins: a putative activation mode of GPCRs?

GPCRs undergo large conformational changes during their activation. Starting from existing X-ray structures, we used Normal Modes Analyses to study the collective motions of the agonist-bound ß2-adrenergic receptor both in its isolated "uncoupled" and G-protein "coupled" conformations. We interestingly observed that the receptor was able to adopt only one major motion in the protein:protein complex. This motion corresponded to an anti-symmetric rotation of both its extra- and intra-cellular parts, with a key role of previously identified highly conserved proline residues. Because this motion was also retrieved when performing NMA on 7 other GPCRs which structures were available, it is strongly suspected to possess a significant biological role, possibly being the "activation mode" of a GPCR when coupled to G-proteins.

Interaction between the elastin peptide VGVAPG and human elastin binding protein.

The elastin binding protein (EBP), a spliced variant of lysosomal ß-galactosidase, is the primary receptor of elastin peptides that have been linked to emphysema, aneurysm and cancer progression. The sequences recognized by EBP share the XGXXPG consensus pattern found in numerous matrix proteins, notably in elastin where the VGVAPG motif is repeated. To delineate the elastin binding site of human EBP, we built a homology model of this protein and docked VGVAPG on its surface. Analysis of this model suggested that Gln-97 and Asp-98 were required for interaction with VGVAPG because they contribute to the definition of a pocket thought to represent the elastin binding site of EBP. Additionally, we proposed that Leu-103, Arg-107, and Glu-137 were essential residues because they could interact with VGVAPG itself. Site-directed mutagenesis experiments at these key positions validated our model. This work therefore provides the first structural data concerning the interaction of the VGVAPG with its cognate receptor. The present structural data should now allow the development of EBP-specific antagonists.

Dissociation of membrane-anchored heterotrimeric G-protein induced by G(α) subunit binding to GTP.

Heterotrimeric G-proteins' activation on the intracellular side of the cell membrane is initiated by stimulation of the G-Protein Coupled Receptors (GPCRs) extra-cellular part. This two-step activation mechanism includes (1) an exchange between GDP and GTP molecules in the G(α) subunit and (2) a dissociation of the whole G(αßγ) complex into two membrane-anchored blocks, namely the isolated G(α) and G(ßγ) subunits. Although X-ray data are available for both inactive G(αßγ):GDP and active G(α):GTP complexes, intermediate steps involved in the molecular mechanism of the dissociation have not yet been addressed at the molecular level. In this study, we first built a membrane-anchored intermediate G(iαßγ):GTP complex. This model was then equilibrated by molecular dynamics simulations before the Targeted Molecular Dynamics (TMD) technique was used to force the G(α) subunit to evolve from its inactive (GDP-bound) to its active (GTP-bound) conformations, as described by available X-ray data. The TMD constraint was applied only to the G(α) subunit so that the resulting global rearrangements acting on the whole G(αßγ):GTP heterotrimer could be analyzed. We showed how these mainly local conformational changes of G(α) could initiate large domain:domain motions of the whole complex, the G(ßγ) behaving as an almost quasi-rigid block. This separation of the two G(α):GTP and G(ßγ) subunits required the loss of several interactions at the G(α):G(ßγ) interface that were reported. This study provided an atomistic view of the crucial intermediate step of the G-proteins activation, e.g., the dissociation, that could hardly be elucidated by the experiment.

GDP release preferentially occurs on the phosphate side in heterotrimeric G-proteins.

After extra-cellular stimulation of G-Protein Coupled Receptors (GPCRs), GDP/GTP exchange appears as the key, rate limiting step of the intracellular activation cycle of heterotrimeric G-proteins. Despite the availability of a large number of X-ray structures, the mechanism of GDP release out of heterotrimeric G-proteins still remains unknown at the molecular level. Starting from the available X-ray structure, extensive unconstrained/constrained molecular dynamics simulations were performed on the complete membrane-anchored Gi heterotrimer complexed to GDP, for a total simulation time overcoming 500 ns. By combining Targeted Molecular Dynamics (TMD) and free energy profiles reconstruction by umbrella sampling, our data suggest that the release of GDP was much more favored on its phosphate side. Interestingly, upon the forced extraction of GDP on this side, the whole protein encountered large, collective motions in perfect agreement with those we described previously including a domain to domain motion between the two ras-like and helical sub-domains of G(α).

Ligands and signaling proteins govern the conformational landscape explored by a G protein-coupled receptor.

The dynamic character of G protein-coupled receptors is essential to their function. However, the details of how ligands stabilize a particular conformation to selectively activate a signaling pathway and how signaling proteins affect this conformational repertoire remain unclear. Using a prototypical peptide-activated class A G protein-coupled receptor (GPCR), the ghrelin receptor, reconstituted as a monomer into lipid discs and labeled with a fluorescent conformational reporter, we demonstrate that ligand efficacy and functional selectivity are directly related to different receptor conformations. Of importance, our data bring direct evidence that distinct effector proteins affect the conformational landscape of the ghrelin receptor in different ways. Whereas G proteins affect the balance between active and inactive receptor substates in favor of the active state, agonist-induced arrestin recruitment is accompanied by a marked change in the structural features of the receptor that adopt a conformation different from that observed in the absence of arrestin. In contrast to G proteins and arrestins, µ-AP2 has no significant effect on the organization of the transmembrane core of the receptor. Such a modulation of a GPCR conformational landscape by pharmacologically distinct ligands and effectors provides insights into the structural bases that decisively affect ligand efficacy and subsequent biological responses. This is also likely to have major implications for the design of drugs activating specific GPCR-associated signaling pathways.