Receptor-bound States Research Articles

Host receptor binding induces conformational changes in a PfEMP1 protein

In this issue of Structure, Raghavan etal. present the cryo-EM structures of the malaria-associated group A PfEMP1 HB3VAR03 head in both host receptor-free and receptor-bound states. The structures are complemented by biophysical analysis and introduce an innovative model in which host-receptor binding induces conformational changes in a PfEMP1 protein.

Structural models of viral insulin-like peptides and their analogs.

The insulin receptor (IR), the insulin-like growth factor-1 receptor (IGF1R), and the insulin/IGF1 hybrid receptors (hybR) are homologous transmembrane receptors. The peptide ligands, insulin and IGF1, exhibit significant structural homology and can bind to each receptor via site-1 and site-2 residues with distinct affinities. The variants of the Iridoviridae virus family show capability in expressing single-chain insulin/IGF1 like proteins, termed viral insulin-like peptides (VILPs), which can stimulate receptors from the insulin family. The sequences of VILPs lacking the central C-domain (dcVILPs) are known, but their structures in unbound and receptor-bound states have not been resolved to date. We report all-atom structural models of three dcVILPs (dcGIV, dcSGIV, and dcLCDV1) and their complexes with the receptors (μIR, μIGF1R, and μhybR), and probed the peptide/receptor interactions in each system using all-atom molecular dynamics (MD) simulations. Based on the nonbonded interaction energies computed between each residue of peptides (insulin and dcVILPs) and the receptors, we provide details on residues establishing significant interactions. The observed site-1 insulin/μIR interactions are consistent with previous experimental studies, and a residue-level comparison of interactions of peptides (insulin and dcVILPs) with the receptors revealed that, due to sequence differences, dcVILPs also establish some interactions distinct from those between insulin and IR. We also designed insulin analogs and report enhanced interactions between some analogs and the receptors.

Decoding the Host-Parasite Protein Interactions Involved in Cerebral Malaria Through Glares of Molecular Dynamics Simulations.

Malaria causes millions of deaths every year. The malaria parasite spends a substantial part of its life cycle inside human erythrocytes. Inside erythrocytes, it synthesizes and displays various proteins onto the erythrocyte surface, such as Plasmodium falciparum erythrocytic membrane protein-1 (PfEMP1). This protein contains cysteine-rich interdomain region (CIDR) domains which have many subtypes based on sequence diversity and can cross-talk with host molecules. The CIDRα1.4 subtype can attach host endothelial protein C receptor (EPCR). This interaction facilitates infected erythrocyte adherence to brain endothelium and subsequent development of cerebral malaria. Through molecular dynamics simulations in conjunction with the molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) method, we explored the mechanism of interaction in the CIDRα1-EPCR complex. We examined the structural behavior of two CIDRα1 molecules (encoded by HB3-isolate var03-gene and IT4-isolate var07-gene) with EPCR unbound and bound (complex) forms. HB3var03CIDRα1 in apo and complexed with EPCR was comparatively more stable than IT4var07CIDRα1. Both of the complexes adopted two distinct conformational energy states. The hydrophobic residues played a crucial role in the binding of both complexes. For HB3var03CIDRα1-EPCR, the dominant energetic components were total polar interactions, while in IT4var07CIDRα1-EPCR, the primary interaction was van der Waals and nonpolar solvation energy. The study also revealed details such as correlated conformational motions and secondary structure evolution. Further, it elucidated various hotspot residues involved in protein-protein recognition. Overall, our study provides additional information on the structural behavior of CIDR molecules in unbound and receptor-bound states, which will help to design potent inhibitors.

The structural basis for the phospholipid remodeling by lysophosphatidylcholine acyltransferase 3

As the major component of cell membranes, phosphatidylcholine (PC) is synthesized de novo in the Kennedy pathway and then undergoes extensive deacylation-reacylation remodeling via Lands’ cycle. The re-acylation is catalyzed by lysophosphatidylcholine acyltransferase (LPCAT) and among the four LPCAT members in human, the LPCAT3 preferentially introduces polyunsaturated acyl onto the sn-2 position of lysophosphatidylcholine, thereby modulating the membrane fluidity and membrane protein functions therein. Combining the x-ray crystallography and the cryo-electron microscopy, we determined the structures of LPCAT3 in apo-, acyl donor-bound, and acyl receptor-bound states. A reaction chamber was revealed in the LPCAT3 structure where the lysophosphatidylcholine and arachidonoyl-CoA were positioned in two tunnels connected near to the catalytic center. A side pocket was found expanding the tunnel for the arachidonoyl CoA and holding the main body of arachidonoyl. The structural and functional analysis provides the basis for the re-acylation of lysophosphatidylcholine and the substrate preference during the reactions.

Dynamics and electrostatics define an allosteric druggable site within the receptor-binding domain of SARS-CoV-2 spike protein.

The pathogenesis of the SARS‐CoV‐2 virus initiates through recognition of the angiotensin‐converting enzyme 2 (ACE2) receptor of the host cells by the receptor‐binding domain (RBD) located at the spikes of the virus. Here, using molecular dynamics simulations, we have demonstrated the allosteric crosstalk within the RBD in the apo‐ and the ACE2 receptor‐bound states, revealing the contribution of the dynamics‐based correlated motions and the electrostatic energy perturbations to this crosstalk. While allostery, based on correlated motions, dominates inherent distal communication in the apo‐RBD, the electrostatic energy perturbations determine favorable pairwise crosstalk within the RBD residues upon binding to ACE2. Interestingly, the allosteric path is composed of residues which are evolutionarily conserved within closely related coronaviruses, pointing toward the biological relevance of the communication and its potential as a target for drug development.

Conformational selection of vasopressin upon V1a receptor binding

The neuropeptide vasopressin (VP) and its three G protein-coupled receptors (V1aR, V1bR and V2R) are of high interest in a wide array of drug discovery programs. V1aR is of particular importance due to its cardiovascular functions and diverse roles in the central nervous system. The structure–activity relationships underpinning ligand-receptor interactions remain however largely unclear, hindering rational drug design. This is not least due to the high structural flexibility of VP in its free as well as receptor-bound states. In this work, we developed a novel approach to reveal features of conformational selectivity upon VP-V1aR complex formation. We employed virtual screening strategies to probe VP’s conformational space for transiently adopted structures that favor binding to V1aR. To this end, we dissected the VP conformational space into three sub-ensembles, each containing distinct structural sets for VP’s three-residue C-terminal tail. We validated the computational results with experimental nuclear magnetic resonance (NMR) data and docked each sub-ensemble to V1aR. We observed that the conformation of VP’s three-residue tail significantly modulated the complex dissociation constants. Solvent-exposed and proline trans-configured VP tail conformations bound to the receptor with three-fold enhanced affinities compared to compacted or cis-configured conformations. The solvent-exposed and more flexible structures facilitated unique interaction patterns between VP and V1aR transmembrane helices 3, 4, and 6 which led to high binding energies. The presented “virtual conformational space screening” approach, integrated with NMR spectroscopy, thus enabled identification and characterization of a conformational selection-type complex formation mechanism that confers novel perspectives on targeting the VP-V1aR interactions at the level of the encounter complex – an aspect that opens novel research avenues for understanding the functionality of the evolutionary selected conformational properties of VP, as well as guidance for ligand design strategies to provide more potent and selective VP analogues.

Forbidden Coherence Transfer of 19F Nuclei to Quantitatively Measure the Dynamics of a CF₃-Containing Ligand in Receptor-Bound States.

The dynamic property of a ligand in the receptor-bound state is an important metric to characterize the interactions in the ligand–receptor interface, and the development of an experimental strategy to quantify the amplitude of motions in the bound state is of importance to introduce the dynamic aspect into structure-guided drug development (SGDD). Fluorine modifications are frequently introduced at the hit-to-lead optimization stage to enhance the binding potency and other characteristics of a ligand. However, the effects of fluorine modifications are generally difficult to predict, owing to the pleiotropic nature of the interactions. In this study, we report an NMR-based approach to experimentally evaluate the local dynamics of trifluoromethyl (CF3)-containing ligands in the receptor-bound states. For this purpose, the forbidden coherence transfer (FCT) analysis, which has been used to study the dynamics of methyl moieties in proteins, was extended to the 19F nuclei of CF3-containing ligands. By applying this CF3–FCT analysis to a model interaction system consisting of a ligand, AST-487, and a receptor, p38α, we successfully quantified the amplitude of the CF3 dynamics in the p38α-bound state. The strategy would bring the CF3-containing ligands within the scope of dynamic SGDD to improve the affinity and specificity for the drug-target receptors.

Probing Gαi1 protein activation at single-amino acid resolution.

We present comprehensive maps at single-amino acid resolution of the residues stabilizing the human Gαi1 subunit in nucleotide- and receptor-bound states. We generated these maps by measuring the effects of alanine mutations on the stability of Gαi1 and the rhodopsin-Gαi1 complex. We identified stabilization clusters in the GTPase and helical domains responsible for structural integrity and the conformational changes associated with activation. In activation cluster I, helices α1 and α5 pack against strands β1-β3 to stabilize the nucleotide-bound states. In the receptor-bound state, these interactions are replaced by interactions between α5 and strands β4-β6. Key residues in this cluster are Y320, which is crucial for the stabilization of the receptor-bound state, and F336, which stabilizes nucleotide-bound states. Destabilization of helix α1, caused by rearrangement of this activation cluster, leads to the weakening of the interdomain interface and release of GDP.

A survey of conformational and energetic changes in G protein signaling

Cell signaling is a fundamental process for all living organisms. G protein-coupled receptors (GPCRs) are a large and diverse group of transmembrane receptors which convert extracellular signals into intracellular responses primarily via coupling to heterotrimeric G proteins. In order to integrate the range of very diverse extracellular signals into a message the cell can recognize and respond to, conformational changes occur that rewire the interactions between the receptor and heterotrimer in a specific and coordinated manner. By interrogating the energetics of these interactions within the individual proteins and across protein-protein interfaces, a communication network between amino acids involved in conformational changes for signaling, is created. To construct this mapping of pairwise interactions in silico, we analyzed the Rhodopsin GPCR coupled to a Gαi1β1γ1 heterotrimer. The structure of this G protein complex was modeled in the receptor-bound and unbound heterotrimeric states as well as the activated, monomeric Gα(GTP) state. From these tertiary structural models, we computed the average pairwise residue-residue interactions and interface energies across ten models of each state using the ROSETTA modeling software suite. Here we disseminate a comprehensive survey of all critical interactions and create intra-protein network communication maps. These networks represent nodes of interaction necessary for G protein activation.

Identification of Receptor Binding-induced Conformational Changes in Non-visual Arrestins

The non-visual arrestins, arrestin-2 and arrestin-3, belong to a small family of multifunctional cytosolic proteins. Non-visual arrestins interact with hundreds of G protein-coupled receptors (GPCRs) and regulate GPCR desensitization by binding active phosphorylated GPCRs and uncoupling them from heterotrimeric G proteins. Recently, non-visual arrestins have been shown to mediate G protein-independent signaling by serving as adaptors and scaffolds that assemble multiprotein complexes. By recruiting various partners, including trafficking and signaling proteins, directly to GPCRs, non-visual arrestins connect activated receptors to diverse signaling pathways. To investigate arrestin-mediated signaling, a structural understanding of arrestin activation and interaction with GPCRs is essential. Here we identified global and local conformational changes in the non-visual arrestins upon binding to the model GPCR rhodopsin. To detect conformational changes, pairs of spin labels were introduced into arrestin-2 and arrestin-3, and the interspin distances in the absence and presence of the receptor were measured by double electron electron resonance spectroscopy. Our data indicate that both non-visual arrestins undergo several conformational changes similar to arrestin-1, including the finger loop moving toward the predicted location of the receptor in the complex as well as the C-tail release upon receptor binding. The arrestin-2 results also suggest that there is no clam shell-like closure of the N- and C-domains and that the loop containing residue 136 (homolog of 139 in arrestin-1) has high flexibility in both free and receptor-bound states.

Sulfatides Partition Disabled-2 in Response to Platelet Activation

BackgroundPlatelets contact each other at the site of vascular injury to stop bleeding. One negative regulator of platelet aggregation is Disabled-2 (Dab2), which is released to the extracellular surface upon platelet activation. Dab2 inhibits platelet aggregation through its phosphotyrosine-binding (PTB) domain by competing with fibrinogen for αIIbβ3 integrin receptor binding by an unknown mechanism.Methodology/Principal FindingsUsing protein-lipid overlay and liposome-binding assays, we identified that the N-terminal region of Dab2, including its PTB domain (N-PTB), specifically interacts with sulfatides. Moreover, we determined that such interaction is mediated by two conserved basic motifs with a dissociation constant (Kd) of 0.6 µM as estimated by surface plasmon resonance (SPR) analysis. In addition, liposome-binding assays combined with mass spectroscopy studies revealed that thrombin, a strong platelet agonist, cleaved N-PTB at a site located between the basic motifs, a region that becomes protected from thrombin cleavage when bound to sulfatides. Sulfatides on the platelet surface interact with coagulation proteins, playing a major role in haemostasis. Our results show that sulfatides recruit N-PTB to the platelet surface, sequestering it from integrin receptor binding during platelet activation. This is a transient recruitment that follows N-PTB internalization by an actin-dependent process.Conclusions/SignificanceOur experimental data support a model where two pools of Dab2 co-exist at the platelet surface, in both sulfatide- and integrin receptor-bound states, and their balance controls the extent of the clotting response.

Crystal Structure of a Prolactin Receptor Antagonist Bound to the Extracellular Domain of the Prolactin Receptor

The crystal structure of the complex between an N-terminally truncated G129R human prolactin (PRL) variant and the extracellular domain of the human prolactin receptor (PRLR) was determined at 2.5A resolution by x-ray crystallography. This structure represents the first experimental structure reported for a PRL variant bound to its cognate receptor. The binding of PRL variants to the PRLR extracellular domain was furthermore characterized by the solution state techniques, hydrogen exchange mass spectrometry, and NMR spectroscopy. Compared with the binding interface derived from mutagenesis studies, the structural data imply that the definition of PRL binding site 1 should be extended to include residues situated in the N-terminal part of loop 1 and in the C terminus. Comparison of the structure of the receptor-bound PRL variant with the structure reported for the unbound form of a similar analogue ( Jomain, J. B., Tallet, E., Broutin, I., Hoos, S., van Agthoven, J., Ducruix, A., Kelly, P. A., Kragelund, B. B., England, P., and Goffin, V. (2007) J. Biol. Chem. 282, 33118-33131 ) demonstrates that receptor-induced changes in the backbone of the four-helix bundle are subtle, whereas large scale rearrangements and structuring occur in the flexible N-terminal part of loop 1. Hydrogen exchange mass spectrometry data imply that the dynamics of the four-helix bundle in solution generally become stabilized upon receptor interaction at binding site 1.

New probes of ligand flexibility in drug design: transferred (13)C CSA-dipolar cross-correlated relaxation at natural abundance.

Understanding the impact of molecular flexibility remains an important outstanding problem in rational drug design. Toward this end, we present new NMR relaxation methods that describe ligand flexibility at the atomic level. Specifically, we measure natural abundance (13)C cross-correlated relaxation parameters for ligands in rapid exchange between the free and receptor-bound states. The rapid exchange transfers the bound state relaxation parameters to the free state, such that a comparison of relaxation rates in the absence and presence of protein receptor yields site-specific information concerning the bound ligand flexibility. We perform these measurements for aromatic carbons, which are highly prevalent in drug-like molecules and demonstrate significant cross-correlated relaxation between the (13)C-(1)H dipole-dipole (DD) and (13)C chemical shift anisotropy (CSA) relaxation mechanisms. Our use of natural abundance measurements addresses the practical difficulties of obtaining isotope-labeled ligands in pharmaceutical research settings. We demonstrate our methods on a small ligand of the 42 kDa kinase domain of the p38 MAP kinase. We show that exchange-transferred cross-correlated relaxation measurements are not only sensitive probes of bound ligand flexibility but also offer complementary advantages over standard R(1) = 1/T(1) and R(2) = 1/T(2) measurements. The ligand flexibility profiles obtained from the relaxation data can help assess the influence of dynamics on ligand potency or pharmacokinetic properties or both, and thereby include inherent molecular flexibility in drug design.

Peptide conformational changes induced by tryptophan–phosphocholine interactions in a micelle

Sodium dodecylsulfate (SDS) and dodecylphosphocholine (DPC) micelles are often used to mimic the membrane- or receptor-bound states of peptides in NMR studies. From the present examination of a 26-residue analog of exendin-4 (TrEX4) by NMR and CD in water, aqueous 30% trifluoroethanol (TFE), and bound to both SDS and DPC micelles, it is clear that these two lipid micelles can yield very different peptide structures. The Trp-cage fold (also observed in 30% TFE) is present when TrEX4 is bound to SDS micelles; however, tertiary structure is absent in the presence of DPC micelles. The loss of tertiary structure is attributed to an energetically favorable interaction (estimated as 2-3 kcal/mol) of the tryptophan side chain with the phosphocholine head groups. These dramatic structural differences suggest that care must be taken when using either SDS or DPC to mimic the membrane- or receptor-bound states.

NMR solution conformation of an antitoxic analogue of alpha-conotoxin GI: identification of a common nicotinic acetylcholine receptor alpha 1-subunit binding surface for small ligands and alpha-conotoxins.

The three-dimensional solution conformation of an 11-residue antitoxic analogue of alpha-conotoxin GI, des-Glu1-[Cys3Ala]-des-Cys13-conotoxin GI (CANPACGRHYS-NH(2), designated "GI-15" henceforth), has been determined using two-dimensional (1)H NMR spectroscopy. The disulfide loop region (1C-6C) and the C-terminal tail (8R-11S) are connected by a flexible hinge formed near 7G, and the pairwise backbone rmsds for the former and the latter are 0.58 and 0.65 A, respectively. Superpositioning GI-15 with the structure of alpha-conotoxin GI shows that the two share an essentially identical fold in the common first disulfide loop region (1C-6C). However, the absence of the second disulfide loop in GI-15 results in segmental motion of the C-terminal half, causing the key receptor subtype selectivity residue 8R (Arg9 in alpha-conotoxin GI) to lose its native spatial orientation. The combined features of structural equivalence in the disulfide loop and a mobile C-terminal tail appear to be responsible for the activity of GI-15 as a competitive antagonist against native toxin. Electrostatic surface potential comparisons of the first disulfide region of GI-15 with other alpha-conotoxins or receptor-bound states of acetylcholine and d-tubocurarine show a common protruding surface that may serve as the minimal binding determinant for the neuromuscular acetylcholine receptor alpha 1-subunit. On the basis of the original "Conus toxin macrosite model" [Olivera, B. M., Rivier, J., Scott, J. K., Hillyard, D. R., and Cruz, L. J. (1991) J. Biol. Chem. 266, 1923-1936], we propose a revised binding model which incorporates these results.