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

Insights into the dynamic interactions at chemokine-receptor interfaces and mechanistic models of chemokine binding

Chemokine receptors are the central signaling hubs of several processes such as cell migration, chemotaxis and cell positioning. In this graphical review, we provide an overview of the structural and mechanistic principles governing chemokine recognition that are currently emerging. Structural models of chemokine-receptor co-complexes with endogenous chemokines, viral chemokines and therapeutics have been resolved that highlight multiple interaction sites, termed as CRS1, CRS1.5 etc. The first site of interaction has been shown to be the N-terminal domain of the receptors (CRS1 site). A large structural flexibility of the N-terminal domain has been reported that was supported by both experimental and simulation studies. Upon chemokine binding, the N-terminal domain appears to show constricted dynamics and opens up to interact with the chemokine via a large interface. The subsequent sites such as CRS1.5 and CRS2 sites have been structurally well resolved although differences arise such as the localization of the N-terminus of the ligand to a major or minor pocket of the orthosteric binding site. Several computational studies have highlighted the dynamic protein-protein interface at the CRS1 site that seemingly appears to resolve the differences in NMR and mutagenesis studies. Interestingly, the differential dynamics at the CRS1 site suggests a mixed model of binding with complex signatures of both conformational selection and induced fit models. Integrative experimental and computational approaches could help unravel the structural basis of promiscuity and specificity in chemokine-receptor binding and open up new avenues of therapeutic design.

Allosteric modulation of the chemokine receptor-chemokine CXCR4-CXCL12 complex by tyrosine sulfation

The chemokine receptor CXCR4 and its cognate ligand CXCL12 mediate pathways that lead to cell migration and chemotaxis. Although the structural details of related receptor-ligand complexes have been resolved, the roles of the N-terminal domain of the receptor and post-translational sulfation that are determinants of ligand selectivity and affinity remain unclear. Here, we analyze the structural dynamics induced by receptor sulfation by combining molecular dynamics, docking and network analysis. The sulfotyrosine residues, 7Ys^N-term, 12Ys^N-term and 21Ys^N-term allow the N-terminal domain of the apo-sulfated receptor to adopt an "open" conformation that appears to facilitate ligand binding. The overall topology of the CXCR4-CXCL12 complex is independent of the sulfation state, but an extensive network of protein-protein interactions characterizes the sulfated receptor, in line with its increased ligand affinity. The altered interactions of sulfotyrosine residues, such as 21Ys^N-term-47R^CXCL12 replacing the 21Y^N-term-13F^CXCL12 interaction, propagate via allosteric pathways towards the receptor lumen. In particular, our results suggest that the experimentally-reported receptor-ligand interactions 262D^6.58-8R^CXCL12 and 277E^7.28-12R^CXCL12 could be dependent on the sulfation state of the receptor and need to be carefully analyzed. Our work is an important step in understanding chemokine-receptor interactions and how post-translational modifications could modulate receptor-ligand complexes.

Conformational plasticity and dynamic interactions of the N-terminal domain of the chemokine receptor CXCR1

The dynamic interactions between G protein-coupled receptors (GPCRs) and their cognate protein partners are central to several cell signaling pathways. For example, the association of CXC chemokine receptor 1 (CXCR1) with its cognate chemokine, interleukin-8 (IL8 or CXCL8) initiates pathways leading to neutrophil-mediated immune responses. The N-terminal domain of chemokine receptors confers ligand selectivity, but unfortunately the conformational dynamics of this intrinsically disordered region remains unresolved. In this work, we have explored the interaction of CXCR1 with IL8 by microsecond time scale coarse-grain simulations, complemented by atomistic models and NMR chemical shift predictions. We show that the conformational plasticity of the apo-receptor N-terminal domain is restricted upon ligand binding, driving it to an open C-shaped conformation. Importantly, we corroborated the dynamic complex sampled in our simulations against chemical shift perturbations reported by previous NMR studies and show that the trends are similar. Our results indicate that chemical shift perturbation is often not a reporter of residue contacts in such dynamic associations. We believe our results represent a step forward in devising a strategy to understand intrinsically disordered regions in GPCRs and how they acquire functionally important conformational ensembles in dynamic protein-protein interfaces.

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 G_s 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-G_s 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 G_s binding modes. They reveal an original receptor-G_s 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.

Perspectives on High-Throughput Ligand/Protein Docking With Martini MD Simulations

Molecular docking is central to rational drug design. Current docking techniques suffer, however, from limitations in protein flexibility and solvation models and by the use of simplified scoring functions. All-atom molecular dynamics simulations, on the other hand, feature a realistic representation of protein flexibility and solvent, but require knowledge of the binding site. Recently we showed that coarse-grained molecular dynamics simulations, based on the most recent version of the Martini force field, can be used to predict protein/ligand binding sites and pathways, without requiring any a priori information, and offer a level of accuracy approaching all-atom simulations. Given the excellent computational efficiency of Martini, this opens the way to high-throughput drug screening based on dynamic docking pipelines. In this opinion article, we sketch the roadmap to achieve this goal.

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.

Recent Advances in Coarse-Grained Models for Biomolecules and Their Applications

Molecular dynamics simulations have emerged as a powerful tool to study biological systems at varied length and timescales. The conventional all-atom molecular dynamics simulations are being used by the wider scientific community in routine to capture the conformational dynamics and local motions. In addition, recent developments in coarse-grained models have opened the way to study the macromolecular complexes for time scales up to milliseconds. In this review, we have discussed the principle, applicability and recent development in coarse-grained models for biological systems. The potential of coarse-grained simulation has been reviewed through state-of-the-art examples of protein folding and structure prediction, self-assembly of complexes, membrane systems and carbohydrates fiber models. The multiscale simulation approaches have also been discussed in the context of their emerging role in unravelling hierarchical level information of biosystems. We conclude this review with the future scope of coarse-grained simulations as a constantly evolving tool to capture the dynamics of biosystems.

Probing the cooperative mechanism of the μ–δ opioid receptor heterodimer by multiscale simulation

The activation-cooperativity of the μ–δ opioid receptor heterodimer was probed by multiscale simulation coupled with a protein structure network.