The
Regulator of G protein signaling (RGS) proteins play a pivotal role in regulation of G protein‐coupled receptor (GPCR) signaling and are therefore becoming an increasingly important therapeutic target. Recently discovered thiadiazolidinone (TDZD) compounds that target cysteine residues have shown different levels of specificities and potencies for the RGS4 protein, thereby suggesting intrinsic differences in dynamics of this protein upon binding of these compounds. In this work, we investigated using atomistic molecular dynamics (MD) simulations the effect of binding of several small‐molecule inhibitors on perturbations and dynamical motions in RGS4. Specifically, we studied two conformational models of RGS4 in which a buried cysteine residue is solvent‐exposed due to side‐chain motions or due to flexibility in neighboring helices. We found that TDZD compounds with aromatic functional groups perturb the RGS4 structure more than compounds with aliphatic functional groups. Moreover, small‐molecules with aromatic functional groups but lacking sulfur atoms only transiently reside within the protein and spontaneously dissociate to the solvent. We further measured inhibitory effects of TDZD compounds using a protein–protein interaction assay on a single‐cysteine RGS4 protein showing trends in potencies of compounds consistent with our simulation studies. Thermodynamic analyses of RGS4 conformations in the apo‐state and on binding to TDZD compounds revealed links between both conformational models of RGS4. The exposure of cysteine side‐chains appears to facilitate initial binding of TDZD compounds followed by migration of the compound into a bundle of four helices, thereby causing allosteric perturbations in the RGS/Gα protein–protein interface.
more » « less- Award ID(s):
- 1508595
- NSF-PAR ID:
- 10462636
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Proteins: Structure, Function, and Bioinformatics
- Volume:
- 87
- Issue:
- 2
- ISSN:
- 0887-3585
- Page Range / eLocation ID:
- p. 146-156
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
mer operon in bacteria encodes a set of proteins and enzymes that impart resistance to environmental mercury toxicity by importing Hg2+and reducing it to volatile Hg(0). Because the reduction occurs in the cytoplasm, mercuric ions must first be transported across the cytoplasmic membrane by one of a few known transporters. MerF is the smallest of these, containing only two transmembrane helices and two pairs of vicinal cysteines that coordinate mercuric ions. In this work, we use molecular dynamics simulations to characterize the dynamics of MerF in its apo and Hg2+‐bound states. We find that the apo state positions one of the cysteine pairs closer to the periplasmic side of the membrane, while in the bound state the same pair approaches the cytoplasmic side. This finding is consistent with the functional requirement of accepting Hg2+from the periplasmic space, sequestering it on acceptance, and transferring it to the cytoplasm. Conformational changes in the TM helices facilitate the functional interaction of the two cysteine pairs. Free‐energy calculations provide a barrier of 16 kcal/mol for the association of the periplasmic Hg2+‐bound protein MerP with MerF and 7 kcal/mol for the subsequent association of MerF's two cysteine pairs. Despite the significant conformational changes required to move the binding site across the membrane, coarse‐grained simulations of multiple copies of MerF support the expectation that it functions as a monomer. Our results demonstrate how conformational changes and binding thermodynamics could lead to such a small membrane protein acting as an ion transporter. Published 2019. This article is a U.S. Government work and is in the public domain in the USA. -
G protein-coupled receptors (GPCRs) represent the largest group of membrane receptors for transmembrane signal transduction. Ligand-induced activation of GPCRs triggers G protein activation followed by various signaling cascades. Understanding the structural and energetic determinants of ligand binding to GPCRs and GPCRs to G proteins is crucial to the design of pharmacological treatments targeting specific conformations of these proteins to precisely control their signaling properties. In this study, we focused on interactions of a prototypical GPCR, beta-2 adrenergic receptor (β 2 AR), with its endogenous agonist, norepinephrine (NE), and the stimulatory G protein (G s ). Using molecular dynamics (MD) simulations, we demonstrated the stabilization of cationic NE, NE(+), binding to β 2 AR by G s protein recruitment, in line with experimental observations. We also captured the partial dissociation of the ligand from β 2 AR and the conformational interconversions of G s between closed and open conformations in the NE(+)–β 2 AR–G s ternary complex while it is still bound to the receptor. The variation of NE(+) binding poses was found to alter G s α subunit (G s α) conformational transitions. Our simulations showed that the interdomain movement and the stacking of G s α α1 and α5 helices are significant for increasing the distance between the G s α and β 2 AR, which may indicate a partial dissociation of G s α The distance increase commences when G s α is predominantly in an open state and can be triggered by the intracellular loop 3 (ICL3) of β 2 AR interacting with G s α, causing conformational changes of the α5 helix. Our results help explain molecular mechanisms of ligand and GPCR-mediated modulation of G protein activation.more » « less
-
The conformational heterogeneity and dynamics of protein side chains contribute to function, but investigating exactly how is hindered by experimental challenges arising from the fast timescales involved and the spatial heterogeneity of protein structures. The potential of two-dimensional infrared (2D IR) spectroscopy for measuring conformational heterogeneity and dynamics with unprecedented spatial and temporal resolution has motivated extensive effort to develop amino acids with functional groups that have frequency-resolved absorptions to serve as probes of their protein microenvironments. We demonstrate the full advantage of the approach by selective incorporation of the probe p -cyanophenylalanine at six distinct sites in a Src homology 3 domain and the application of 2D IR spectroscopy to site-specifically characterize heterogeneity and dynamics and their contribution to cognate ligand binding. The approach revealed a wide range of microenvironments and distinct responses to ligand binding, including at the three adjacent, conserved aromatic residues that form the recognition surface of the protein. Molecular dynamics simulations performed for all the labeled proteins provide insight into the underlying heterogeneity and dynamics. Similar application of 2D IR spectroscopy and site-selective probe incorporation will allow for the characterization of heterogeneity and dynamics of other proteins, how heterogeneity and dynamics are affected by solvation and local structure, and how they might contribute to biological function.more » « less
-
Abstract G-proteins are molecular on–off switches that are involved in transmitting a variety of extracellular signals to their intracellular targets. In animal and yeast systems, the switch property is encoded through nucleotides: a GDP-bound state is the “off-state” and the GTP-bound state is the “on-state”. The G-protein cycle consists of the switch turning on through nucleotide exchange facilitated by a G-protein coupled receptor and the switch turning off through hydrolysis of GTP back to GDP, facilitated by a protein designated REGULATOR OF G SIGNALING 1 (RGS). In plants, G-protein signaling dramatically differs from that in animals and yeast. Despite stringent conservation of the nucleotide binding and catalytic structures over the 1.6 billion years that separate the evolution of plants and animals, genetic and biochemical data indicate that nucleotide exchange is less critical for this switch to operate in plants. Also, the loss of the single RGS protein in Arabidopsis (Arabidopsis thaliana) confers unexpectedly weaker phenotypes consistent with a diminished role for the G cycle, at least under static conditions. However, under dynamic conditions, genetic ablation of RGS in Arabidopsis results in a strong phenotype. We explore explanations to this conundrum by formulating a mathematical model that takes into account the accruing evidence for the indispensable role of phosphorylation in G-protein signaling in plants and that the G-protein cycle is needed to process dynamic signal inputs. We speculate that the plant G-protein cycle and its attendant components evolved to process dynamic signals through signaling modulation rather than through on–off, switch-like regulation of signaling. This so-called change detection may impart greater fitness for plants due to their sessility in a dynamic light, temperature, and pest environment.more » « less
-
more » « less
Abstract Heterotrimeric G-protein complexes comprising Gα-, Gβ-, and Gγ-subunits and the regulator of G-protein signaling (RGS) are conserved across most eukaryotic lineages. Signaling pathways mediated by these proteins influence overall growth, development, and physiology. In plants, this protein complex has been characterized primarily from angiosperms with the exception of spreading-leaved earth moss (Physcomitrium patens) and Chara braunii (charophytic algae). Even within angiosperms, specific G-protein components are missing in certain species, whereas unique plant-specific variants—the extra-large Gα (XLGα) and the cysteine-rich Gγ proteins—also exist. The distribution and evolutionary history of G-proteins and their function in nonangiosperm lineages remain mostly unknown. We explored this using the wealth of available sequence data spanning algae to angiosperms representing extant species that diverged approximately 1,500 million years ago, using BLAST, synteny analysis, and custom-built Hidden Markov Model profile searches. We show that a minimal set of components forming the XLGαβγ trimer exists in the entire land plant lineage, but their presence is sporadic in algae. Additionally, individual components have distinct evolutionary histories. The XLGα exhibits many lineage-specific gene duplications, whereas Gα and RGS show several instances of gene loss. Similarly, Gβ remained constant in both number and structure, but Gγ diverged before the emergence of land plants and underwent changes in protein domains, which led to three distinct subtypes. These results highlight the evolutionary oddities and summarize the phyletic patterns of this conserved signaling pathway in plants. They also provide a framework to formulate pertinent questions on plant G-protein signaling within an evolutionary context.
