Mice with GNAO1 R209H Movement Disorder Variant Display Hyperlocomotion Alleviated by Risperidone.

Neurodevelopmental disorder with involuntary movements (Online Mendelian Inheritance in Man: 617493) is a severe, early onset neurologic condition characterized by a delay in psychomotor development, hypotonia, and hyperkinetic involuntary movements. Heterozygous de novo mutations in the GNAO1 gene cause neurodevelopmental disorder with involuntary movements. Gα o, the gene product of GNAO1, is the alpha subunit of Go, a member of the heterotrimeric Gi/o family of G proteins. Go is found abundantly throughout the brain, but the pathophysiological mechanisms linking Gα o functions to clinical manifestations of GNAO1-related disorders are still poorly understood. One of the most common mutant alleles among the GNAO1 encephalopathies is the c.626G>A or p.Arg209His (R209H) mutation. We developed heterozygous knock-in Gnao1 +/R209H mutant mice using CRISPR/Cas9 methodology to assess whether a mouse model could replicate aspects of the neurodevelopmental disorder with involuntary movements clinical pattern. Mice carrying the R209H mutation exhibited increased locomotor activity and a modest gait abnormality at 8-12 weeks. In contrast to mice carrying other mutations in Gnao1, the Gnao1 +/R209H mice did not show enhanced seizure susceptibility. Levels of protein expression in multiple brain regions were unchanged from wild-type (WT) mice, but the nucleotide exchange rate of mutant R209H Gα o was 6.2× faster than WT. The atypical neuroleptic risperidone has shown efficacy in a patient with the R209H mutation. It also alleviated the hyperlocomotion phenotype observed in our mouse model but suppressed locomotion in WT mice as well. In this study, we show that Gnao1 +/R209H mice mirror elements of the patient phenotype and respond to an approved pharmacological agent. SIGNIFICANCE STATEMENT: Children with de novo mutations in the GNAO1 gene may present with movement disorders with limited effective therapeutic options. The most common mutant variant seen in children with GNAO1-associated movement disorder is R209H. Here we show, using a novel Gnao1 +/R209H mouse, that there is a clear behavioral phenotype that is suppressed by risperidone. However, risperidone also affects wild-type mouse activity, so its effects are not selective for the GNAO1-associated movement disorder.

Loss-of-Function Mutations in Human Regulator of G Protein Signaling RGS2 Differentially Regulate Pharmacological Reactivity of Resistance Vasculature.

Regulator of G protein signaling 2 (RGS2) plays a role in reducing vascular contraction and promoting relaxation due to its GTPase accelerating protein activity toward Gαq. Previously, we identified four human loss-of-function (LOF) mutations in RGS2 (Q2L, D40Y, R44H, and R188H). This study aimed to investigate whether those RGS2 LOF mutations disrupt the ability of RGS2 to regulate vascular reactivity. Isolated mesenteric arteries (MAs) from RGS2-/- mice showed an elevated contractile response to 5 nM angiotensin II and a loss of acetylcholine (ACh)-mediated vasodilation. Reintroduction of a wild-type (WT) RGS2-GFP plasmid into RGS2-/- MAs suppressed the vasoconstrictor response to angiotensin II. RGS2 LOF mutants failed to suppress the angiotensin II constriction response compared with RGS2 WT. In contrast, ACh-mediated vasoconstriction was restored by expression of RGS2 WT, D40Y, and R44H but not by RGS2 Q2L or R188H. Phosphorylation of RGS2 D40Y and R44H by protein kinase G (PKG) may explain their maintained function to support relaxation in MAs. This is supported by phosphomimetic mutants and suppression of vasorelaxation mediated by RGS2 D40Y by a PKG inhibitor. These results demonstrate that RGS2 attenuates vasoconstriction in MAs and that RGS2 LOF mutations cannot carry out this effect. Among them, the Q2L and R188H mutants supported less relaxation to ACh, whereas relaxation mediated by the D40Y and R44H mutant proteins was equal to that with WT protein. Phosphorylation of RGS2 by PKG appears to contribute to this vasorelaxation. These results provide insights for precision medicine targeting the rare individuals carrying these RGS2 mutations. SIGNIFICANCE STATEMENT: Regulator of G protein signaling 2 (RGS2) has been implicated in the control of blood pressure; rare mutations in the RGS2 gene have been identified in large-scale human gene sequencing studies. Four human mutations in RGS2 that cause loss of function (LOF) in cell-based assays were examined in isolated mouse arteries for effects on both vasoconstriction and vasodilation. All mutants showed the expected LOF effects in suppressing vasoconstriction. Surprisingly, the D40Y and R44H mutant RGS2 showed normal control of vasodilation. We propose that this is due to rescue of the mislocalization phenotype of these two mutants by nitric oxide-mediated/protein kinase G-dependent phosphorylation. These mechanisms may guide drug discovery or drug repurposing efforts for hypertension by enhancing RGS2 function.

An Interhelical Salt Bridge Controls Flexibility and Inhibitor Potency for Regulators of G-protein Signaling Proteins 4, 8, and 19.

Regulators of G-protein signaling (RGS) proteins modulate receptor signaling by binding to activated G-protein α-subunits, accelerating GTP hydrolysis. Selective inhibition of RGS proteins increases G-protein activity and may provide unique tissue specificity. Thiadiazolidinones (TDZDs) are covalent inhibitors that act on cysteine residues to inhibit RGS4, RGS8, and RGS19. There is a correlation between protein flexibility and potency of inhibition by the TDZD 4-[(4- fluorophenyl)methyl]-2-(4-methylphenyl)-1,2,4-thiadiazolidine-3,5-dione (CCG-50014). In the context of a single conserved cysteine residue on the α 4 helix, RGS19 is the most flexible and most potently inhibited by CCG-50014, followed by RGS4 and RGS8. In this work, we identify residues responsible for differences in both flexibility and potency of inhibition among RGS isoforms. RGS19 lacks a charged residue on the α 4 helix that is present in RGS4 and RGS8. Introducing a negative charge at this position (L118D) increased the thermal stability of RGS19 and decreased the potency of inhibition of CCG-50014 by 8-fold. Mutations eliminating salt bridge formation in RGS8 and RGS4 decreased thermal stability in RGS8 and increased potency of inhibition of both RGS4 and RGS8 by 4- and 2-fold, respectively. Molecular dynamics simulations with an added salt bridge in RGS19 (L118D) showed reduced RGS19 flexibility. Hydrogen-deuterium exchange studies showed striking differences in flexibility in the α 4 helix of RGS4, 8, and 19 with salt bridge-modifying mutations. These results show that the α 4 salt bridge-forming residue controls flexibility in several RGS isoforms and supports a causal relationship between RGS flexibility and the potency of TDZD inhibitors. SIGNIFICANCE STATEMENT: Inhibitor potency is often viewed in relation to the static structure of a target protein binding pocket. Using both experimental and computation studies we assess determinants of dynamics and inhibitor potency for three different RGS proteins. A single salt bridge-forming residue determines differences in flexibility between RGS isoforms; mutations either increase or decrease protein motion with correlated alterations in inhibitor potency. This strongly suggests a causal relationship between RGS protein flexibility and covalent inhibitor potency.

Interplay of cysteine exposure and global protein dynamics in small-molecule recognition by a regulator of G-protein signaling protein.

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.

Interpreting Hydrogen-Deuterium Exchange Events in Proteins Using Atomistic Simulations: Case Studies on Regulators of G-Protein Signaling Proteins.

Hydrogen-deuterium exchange (HDX) experiments are widely used in studies of protein dynamics. To predict the propensity of amide hydrogens for exchange with deuterium, several models have been reported in which computations of amide-hydrogen protection factors are carried out using molecular dynamics (MD) simulations. Given significant variation in the criteria used in different models, the robustness and broader applicability of these models to other proteins, especially homologous proteins showing distinct amide-exchange patterns, remains unknown. The sensitivity of the predictions when MD simulations are conducted with different force-fields is yet to tested and quantified. Using MD simulations and experimental HDX data on three homologous signaling proteins, we report detailed studies quantifying the performance of seven previously reported models (M1-M7) of two general types: empirical and fractional-population models. We find that the empirical models show inconsistent predictions but predictions of the fractional population models are robust. Contrary to previously reported work, we find that the solvent-accessible surface area of amide hydrogens is a useful metric when combined with a new metric defining the distances of amide hydrogens from the first polar atoms in proteins. On the basis of this, we report two new models, one empirical (M8) and one population-based (M9). We find strong protection of amide hydrogens from solvent exchange both within the stable helical motifs and also in the interhelical loops. We further observe that the exchange-competent states of amide hydrogens occur on the sub 100 ps time-scale via localized fluctuations, and such states among amides of a given protein do not appear to show any cooperativity or allosteric coupling.

Differential Protein Dynamics of Regulators of G-Protein Signaling: Role in Specificity of Small-Molecule Inhibitors.

Small-molecule inhibitor selectivity may be influenced by variation in dynamics among members of a protein family. Regulator of G-protein Signaling (RGS) proteins are a family that plays a key role in G-Protein Coupled Receptor (GPCR) signaling by binding to active Gα subunits and accelerating GTP hydrolysis, thereby terminating activity. Thiadiazolidinones (TDZDs) inhibit the RGS-Gα interaction by covalent modification of cysteine residues in RGS proteins. Some differences in specificity may be explained by differences in the complement of cysteines among RGS proteins. However, key cysteines shared by RGS proteins inhibited by TDZDs are not exposed on the protein surface, and differences in potency exist among RGS proteins containing only buried cysteines. We hypothesize that differential exposure of buried cysteine residues among RGS proteins partially drives TDZD selectivity. Hydrogen-deuterium exchange (HDX) studies and molecular dynamics (MD) simulations were used to probe the dynamics of RGS4, RGS8, and RGS19, three RGS proteins inhibited at a range of potencies by TDZDs. When these proteins were mutated to contain a single, shared cysteine, RGS19 was found to be most potently inhibited. HDX studies revealed differences in α4 and α6 helix flexibility among RGS isoforms, with particularly high flexibility in RGS19. This could cause differences in cysteine exposure and lead to differences in potency of TDZD inhibition. MD simulations of RGS proteins revealed motions that correspond to solvent exposure observed in HDX, providing further evidence for a role of protein dynamics in TDZD selectivity.

The DRY motif and the four corners of the cubic ternary complex model.

Recent structural data on GPCRs using a variety of spectroscopic approaches suggest that GPCRs adopt a dynamic conformational landscape, with ligands stabilizing subsets of these states to activate one or more downstream signaling effectors. A key outstanding question posed by this emerging dynamic structural model of GPCRs is what states, active, inactive, or intermediate are captured by the numerous crystal structures of GPCRs complexed with a variety of agonists, partial agonists, and antagonists. In the early nineties the discovery of inverse agonists and constitutive activity led to the idea that the active receptor state (R⁎) is an intrinsic property of the receptor itself rather than of the RG complex, eventually leading to the formulation of the cubic ternary complex model (CTC). Here, by a careful analysis of a series of data obtained with a number of mutants of the highly conserved E/DRY motif, we show evidences for the existence of all the receptor states theorized by the CTC, four 'uncoupled (R, R⁎ and HR and HR⁎), and, consequently four 'coupled' (RG, R⁎G, HRG and HR⁎G). The E/DRY motif located at the cytosolic end of transmembrane helix III of Class A GPCRs has been widely studied and analyzed because it forms a network of interactions believed to lock receptors in the inactive conformation (R), and, thus, to play a key role in receptor activation. Our conclusions are supported by recent crystal and NMR spectra, as well as by results obtained with two prototypical GPCRs using a new FRET technology that de-couples G protein binding to the receptor from signal transduction. Thus, despite its complexity and limitations, we propose that the CTC is a useful framework to reconcile pharmacological, biochemical and structural data.