Conformational dynamics of a regulator of G-protein signaling protein reveals a mechanism of allosteric inhibition by a small molecule

The roles of RGS proteins in cardiometabolic disease

G protein-coupled receptors (GPCRs) are the most prominent receptors on the surface of the cell and play a central role in the regulation of cardiac and metabolic functions. GPCRs transmit extracellular stimuli to the interior of the cells by activating one or more heterotrimeric G proteins. The duration and intensity of G protein-mediated signalling are tightly controlled by a large array of intracellular mediators, including the regulator of G protein signalling (RGS) proteins. RGS proteins selectively promote the GTPase activity of a subset of Gα subunits, thus serving as negative regulators in a pathway-dependent manner. In the current review, we summarise the involvement of RGS proteins in cardiometabolic function with a focus on their tissue distribution, mechanisms of action and dysregulation under various disease conditions. We also discuss the potential therapeutic applications for targeting RGS proteins in treating cardiometabolic conditions and current progress in developing RGS modulators. LINKED ARTICLES: This article is part of a themed issue Therapeutic Targeting of G Protein-Coupled Receptors: hot topics from the Australasian Society of Clinical and Experimental Pharmacologists and Toxicologists 2021 Virtual Annual Scientific Meeting. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v181.14/issuetoc.

Identification of Core Allosteric Sites through Temperature- and Nucleus-Invariant Chemical Shift Covariance

Allosteric regulation is essential to control biological function. In addition, allosteric sites offer a promising venue for selective drug targeting. However, accurate mapping of allosteric sites remains challenging since allostery relies on often subtle, yet functionally relevant, structural and dynamical changes. A viable approach proposed to overcome such challenge is the chemical shift covariance analysis (CHESCA). Although CHESCA offers an exhaustive map of allosteric networks, it is critical to define the core allosteric sites to be prioritized in subsequent functional studies or the design of allosteric drugs. Here we propose two new CHESCA-based methodologies, called temperature CHESCA (T-CHESCA) and CLASS-CHESCA, aimed at narrowing down allosteric maps to the core allosteric residues. Both T- and CLASS-CHESCAs rely on the invariance of core inter-residue correlations to changes in the chemical shifts of the active and inactive conformations interconverting in fast exchange. In the T-CHESCA the chemical shifts of such states are modulated through temperature changes, while in the CLASS-CHESCA through variations in the spin-active nuclei involved in pairwise correlations. The T- and CLASS-CHESCAs as well as complete-linkage CHESCA were applied to the cAMP-binding domain of the exchange protein directly activated by cAMP (EPAC), which serves as a prototypical allosteric switch. Residues consistently identified by the three CHESCA methods were found in previously identified EPAC allosteric core sites. Hence, the T-, CLASS- and CL-CHESCA provide a toolset to establish allosteric site hierarchy and triage allosteric sites to be further analyzed by mutations and functional assays. Furthermore, the core allosteric networks selectively revealed through T- and CLASS-CHESCA are expected to facilitate the mechanistic understanding of disease-related mutations and the design of selective allosteric modulators.

Mixed-solvent molecular dynamics simulation-based discovery of a putative allosteric site on regulator of G protein signaling 4

Regulator of G protein signaling 4 (RGS4) is an intracellular protein that binds to the G_α subunit ofheterotrimeric G proteins and aids in terminating G protein coupled receptor signaling. RGS4 has been implicated in pain, schizophrenia, and the control of cardiac contractility. Inhibitors of RGS4 have been developed but bind covalently to cysteine residues on the protein. Therefore, we sought to identify alternative druggable sites on RGS4 using mixed-solvent molecular dynamics simulations, which employ low concentrations of organic probes to identify druggable hotspots on the protein. Pseudo-ligands were placed in consensus hotspots, and perturbation with normal mode analysis led to the identification and characterization of a putative allosteric site, which would be invaluable for structure-based drug design of non-covalent, small molecule inhibitors. Future studies on the mechanism of this allostery will aid in the development of novel therapeutics targeting RGS4.

Allosteric Pathways Originating at Cysteine Residues in Regulators of G-Protein Signaling Proteins

Regulators of G-protein signaling (RGS) proteins play a central role in modulating signaling via G-protein coupled receptors (GPCRs). Specifically, RGS proteins bind to activated Gα subunits in G-proteins, accelerate the GTP hydrolysis, and thereby rapidly dampen GPCR signaling. Therefore, covalent molecules targeting conserved cysteine residues among RGS proteins have emerged as potential candidates to inhibit the RGS/Gα protein-protein interaction and enhance GPCR signaling. Although these inhibitors bind to conserved cysteine residues among RGS proteins, we have previously suggested [J. Am. Chem. Soc. 2018;140:3454–3460] that their potencies and specificities are related to differential protein dynamics among RGS proteins. Using data from all-atom molecular dynamics simulations, we reveal these differences in dynamics of RGS proteins by partitioning the protein structural space into a network of communities that allow allosteric signals to propagate along unique pathways originating at inhibitor binding sites and terminating at the RGS/Gα protein-protein interface.

Engineering and application of a biosensor with focused ligand specificity

Cell factories converting bio-based precursors to chemicals present an attractive avenue to a sustainable economy, yet screening of genetically diverse strain libraries to identify the best-performing whole-cell biocatalysts is a low-throughput endeavor. For this reason, transcriptional biosensors attract attention as they allow the screening of vast libraries when used in combination with fluorescence-activated cell sorting (FACS). However, broad ligand specificity of transcriptional regulators (TRs) often prohibits the development of such ultra-high-throughput screens. Here, we solve the structure of the TR LysG of Corynebacterium glutamicum, which detects all three basic amino acids. Based on this information, we follow a semi-rational engineering approach using a FACS-based screening/counterscreening strategy to generate an l-lysine insensitive LysG-based biosensor. This biosensor can be used to isolate l-histidine-producing strains by FACS, showing that TR engineering towards a more focused ligand spectrum can expand the scope of application of such metabolite sensors.

Transcriptional biosensors represent powerful tools for the screening of vast strain libraries, but the broad ligand specificity of some transcriptional regulators (TRs) can prohibit such applications. Here authors present the engineering of a LysG-based biosensor with a focused ligand specificity to isolate L-histidine-producing strains.

Regulator of G-protein signaling (RGS) proteins as drug targets: Progress and future potentials

G protein-coupled receptors (GPCRs) play critical roles in regulating processes such as cellular homeostasis, responses to stimuli, and cell signaling. Accordingly, GPCRs have long served as extraordinarily successful drug targets. It is therefore not surprising that the discovery in the mid-1990s of a family of proteins that regulate processes downstream of GPCRs generated great excitement in the field. This finding enhanced the understanding of these critical signaling pathways and provided potentially new targets for pharmacological intervention. These regulators of G-protein signaling (RGS) proteins were viewed by many as nodes downstream of GPCRs that could be targeted with small molecules to tune signaling processes. In this review, we provide a brief overview of the discovery of RGS proteins and of the gradual and continuing discovery of their roles in disease states, focusing particularly on cancer and neurological disorders. We also discuss high-throughput screening efforts that have led to the discovery first of peptide-based and then of small-molecule inhibitors targeting a subset of the RGS proteins. We explore the unique mechanisms of RGS inhibition these chemical tools have revealed and highlight the most up-to-date studies using these tools in animal experiments. Finally, we discuss the future opportunities in the field, as there are clearly more avenues left to be explored and potentials to be realized.

Structural Analysis and Dynamic Processes of the Transmembrane Segment Inside Different Micellar Environments-Implications for the TM4 Fragment of the Bilitranslocase Protein

The transmembrane (TM) proteins are gateways for molecular transport across the cell membrane that are often selected as potential targets for drug design. The bilitranslocase (BTL) protein facilitates the uptake of various anions, such as bilirubin, from the blood into the liver cells. As previously established, there are four hydrophobic transmembrane segments (TM1-TM4), which constitute the structure of the transmembrane channel of the BTL protein. In our previous studies, the 3D high-resolution structure of the TM2 and TM3 transmembrane fragments of the BTL in sodium dodecyl sulfate (SDS) micellar media were solved using Nuclear Magnetic Resonance (NMR) spectroscopy and molecular dynamics simulations (MD). The high-resolution 3D structure of the fourth transmembrane region (TM4) of the BTL was evaluated using NMR spectroscopy in two different micellar media, anionic SDS and zwitterionic DPC (dodecylphosphocholine). The presented experimental data revealed the existence of an α -helical conformation in the central part of the TM4 in both micellar media. In the case of SDS surfactant, the α -helical conformation is observed for the Pro258-Asn269 region. The use of the zwitterionic DPC micelle leads to the formation of an amphipathic α -helix, which is characterized by the extension of the central α -helix in the TM4 fragment to Phe257-Thr271. The complex character of the dynamic processes in the TM4 peptide within both surfactants was analyzed based on the relaxation data acquired on 15 N and 31 P isotopes. Contrary to previously published and present observations in the SDS micelle, the zwitterionic DPC environment leads to intensive low-frequency molecular dynamic processes in the TM4 fragment.

Emerging roles of allosteric modulators in the regulation of protein-protein interactions (PPIs): A new paradigm for PPI drug discovery

Protein-protein interactions (PPIs) are closely implicated in various types of cellular activities and are thus pivotal to health and disease states. Given their fundamental roles in a wide range of biological processes, the modulation of PPIs has enormous potential in drug discovery. However, owing to the general properties of large, flat, and featureless interfaces of PPIs, previous attempts have demonstrated that the generation of therapeutic agents targeting PPI interfaces is challenging, rendering them almost "undruggable" for decades. To date, rapid progress in chemical and structural biology techniques has promoted the exploitation of allostery as a novel approach in drug discovery. By attaching to allosteric sites that are topologically and spatially distinct from PPI interfaces, allosteric modulators can achieve improved physiochemical properties. Thus, allosteric modulators may represent an alternative strategy to target intractable PPIs and have attracted intense pharmaceutical interest. In this review, we first briefly introduce the characteristics of PPIs and then present different approaches for investigating PPIs, as well as the latest methods for modulating PPIs. Importantly, we comprehensively review the recent progress in the development of allosteric modulators to inhibit or stabilize PPIs. Finally, we conclude with future perspectives on the discovery of allosteric PPI modulators, especially the application of computational methods to aid in allosteric PPI drug discovery.

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.

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.

Screen Targeting Lung and Prostate Cancer Oncogene Identifies Novel Inhibitors of RGS17 and Problematic Chemical Substructures

Regulator of G protein signaling (RGS) proteins temporally regulate heterotrimeric G protein signaling cascades elicited by G protein-coupled receptor activation and thus are essential for cell homeostasis. The dysregulation of RGS protein expression has been linked to several pathologies, spurring discovery efforts to identify small-molecule inhibitors of these proteins. Presented here are the results of a high-throughput screening (HTS) campaign targeting RGS17, an RGS protein reported to be inappropriately upregulated in several cancers. A screen of over 60,000 small molecules led to the identification of five hit compounds that inhibit the RGS17-Gα_o protein-protein interaction. Chemical and biochemical characterization demonstrated that three of these hits inhibited the interaction through the decomposition of parent compound into reactive products under normal chemical library storage/usage conditions. Compound substructures susceptible to decomposition are reported and the decomposition process characterized, adding to the armamentarium of tools available to the screening field, allowing for the conservation of resources in follow-up efforts and more efficient identification of potentially decomposed compounds. Finally, analogues of one hit compound were tested, and the results establish the first ever structure-activity relationship (SAR) profile for a small-molecule inhibitor of RGS17.

Natural Products Discovered in a High-Throughput Screen Identified as Inhibitors of RGS17 and as Cytostatic and Cytotoxic Agents for Lung and Prostate Cancer Cell Lines

Regulator of G Protein Signaling (RGS) 17 is an overexpressed promoter of cancer survival in lung and prostate tumors, the knockdown of which results in decreased tumor cell proliferation in vitro. Identification of drug-like molecules inhibiting this protein could ameliorate the RGS17's pro-tumorigenic effect. Using high-throughput screening, a chemical library containing natural products was interrogated for inhibition of the RGS17-Gα_o interaction. Initial hits were verified in control and counter screens. Leads were characterized via biochemical, mass spectrometric, Western blot, microscopic, and cytotoxicity measures. Four known compounds (1-4) were identified with IC₅₀ values ranging from high nanomolar to low micromolar. Three compounds were extensively characterized biologically, demonstrating cellular activity determined by confocal microscopy, and two compounds were assessed via ITC exhibiting high nanomolar to low micromolar dissociation constants. The compounds were found to have a cysteine-dependent mechanism of binding, verified through site-directed mutagenesis and cysteine reactivity assessment. Two compounds, sanguinarine (1) and celastrol (2), were found to be cytostatic against lung and prostate cancer cell lines and cytotoxic against prostate cancer cell lines in vitro, although the dependence of RGS17 on these phenomena remains elusive, a result that is perhaps not surprising given the multimodal cytostatic and cytotoxic activities of many natural products.

Development of a bimolecular luminescence complementation assay for RGS: G protein interactions in cells

Cell based assessment tools and screening platforms are the preferred paradigm for small molecule identification and validation due to selectively identifying molecules with cellular activity and validation of compound activity against target proteins in their native environment. With respect to Regulator of G Protein Signaling (RGS) proteins, current cell based methodologies are either low throughput or monitor downstream signaling consequences. The increasing number of reports indicating RGS function in various disease pathogeneses highlights the need for a robust RGS inhibitor discovery and characterization paradigm. Promega's NanoBit Protein Complementation Assay utilizes NanoLuc, an engineered luciferase with enhanced luminescence characteristics which allow for both robust and kinetic assessment of protein interaction formation and disruption. Here we characterized 15 separate RGS: G protein interactions using this system. The binding profile of RGS: Gα interactions correlates to prior published biochemical binding profiles of these proteins. Additionally, we demonstrated this system is suitable for high throughput screening efforts via calculation of Z-factors for three of the interactions and demonstrated that a known small molecule inhibitor of RGS4 disrupts the RGS4: Gα_i1 protein-protein interaction. In conclusion, the NanoBit Protein Complementation Assay holds promise as a robust platform for discovery and characterization of RGS inhibitors.

Building Graphs To Describe Dynamics, Kinetics, and Energetics in the d-ALa:d-Lac Ligase VanA

The d-Ala:d-Lac ligase, VanA, plays a critical role in the resistance of vancomycin. Indeed, it is involved in the synthesis of a peptidoglycan precursor, to which vancomycin cannot bind. The reaction catalyzed by VanA requires the opening of the so-called “ω-loop”, so that the substrates can enter the active site. Here, the conformational landscape of VanA is explored by an enhanced sampling approach: the temperature-accelerated molecular dynamics (TAMD). Analysis of the molecular dynamics (MD) and TAMD trajectories recorded on VanA permits a graphical description of the structural and kinetics aspects of the conformational space of VanA, where the internal mobility and various opening modes of the ω-loop play a major role. The other important feature is the correlation of the ω-loop motion with the movements of the opposite domain, defined as containing the residues A149–Q208. Conformational and kinetic clusters have been determined and a path describing the ω-loop opening was extracted from these clusters. The determination of this opening path, as well as the relative importance of hydrogen bonds along the path, permit one to propose some key residue interactions for the kinetics of the ω-loop opening.

Advances in NMR Methods To Map Allosteric Sites: From Models to Translation

The last five years have witnessed major developments in the understanding of the allosteric phenomenon, broadly defined as coupling between remote molecular sites. Such advances have been driven not only by new theoretical models and pharmacological applications of allostery, but also by progress in the experimental approaches designed to map allosteric sites and transitions. Among these techniques, NMR spectroscopy has played a major role given its unique near-atomic resolution and sensitivity to the dynamics that underlie allosteric couplings. Here, we highlight recent progress in the NMR methods tailored to investigate allostery with the goal of offering an overview of which NMR approaches are best suited for which allosterically relevant questions. The picture of the allosteric "NMR toolbox" is provided starting from one of the simplest models of allostery (i.e., the four-state thermodynamic cycle) and continuing to more complex multistate mechanisms. We also review how such an "NMR toolbox" has assisted the elucidation of the allosteric molecular basis for disease-related mutations and the discovery of novel leads for allosteric drugs. From this overview, it is clear that NMR plays a central role not only in experimentally validating transformative theories of allostery, but also in tapping the full translational potential of allosteric systems.

NMR Methods to Study Dynamic Allostery

Nuclear magnetic resonance (NMR) spectroscopy provides a unique toolbox of experimental probes for studying dynamic processes on a wide range of timescales, ranging from picoseconds to milliseconds and beyond. Along with NMR hardware developments, recent methodological advancements have enabled the characterization of allosteric proteins at unprecedented detail, revealing intriguing aspects of allosteric mechanisms and increasing the proportion of the conformational ensemble that can be observed by experiment. Here, we present an overview of NMR spectroscopic methods for characterizing equilibrium fluctuations in free and bound states of allosteric proteins that have been most influential in the field. By combining NMR experimental approaches with molecular simulations, atomistic-level descriptions of the mechanisms by which allosteric phenomena take place are now within reach.

Free energy and hidden barriers of the β-sheet structure of prion protein

On-the-fly free-energy parametrization is a new collective variable biasing approach akin to metadynamics with one important distinction: rather than acquiring an accelerated distribution via a history-dependent bias potential, sampling on this distribution is achieved from the beginning of the simulation using temperature-accelerated molecular dynamics. In the present work, we compare the performance of both approaches to compute the free-energy profile along a scalar collective variable measuring the H-bond registry of the β-sheet structure of the mouse Prion protein. Both methods agree on the location of the free-energy minimum, but free-energy profiles from well-tempered metadynamics are subject to a much higher degree of statistical noise due to hidden barriers. The sensitivity of metadynamics to hidden barriers is shown to be a consequence of the history dependence of the bias potential, and we detail the nature of these barriers for the prion β-sheet. In contrast, on-the-fly parametrization is much less sensitive to these barriers and thus displays improved convergence behavior relative to that of metadynamics. While hidden barriers are a frequent and central issue in free-energy methods, on-the-fly free-energy parametrization appears to be a robust and preferable method to confront this issue.

Extended Phase-Space Methods for Enhanced Sampling in Molecular Simulations: A Review

Molecular Dynamics simulations are a powerful approach to study biomolecular conformational changes or protein-ligand, protein-protein, and protein-DNA/RNA interactions. Straightforward applications, however, are often hampered by incomplete sampling, since in a typical simulated trajectory the system will spend most of its time trapped by high energy barriers in restricted regions of the configuration space. Over the years, several techniques have been designed to overcome this problem and enhance space sampling. Here, we review a class of methods that rely on the idea of extending the set of dynamical variables of the system by adding extra ones associated to functions describing the process under study. In particular, we illustrate the Temperature Accelerated Molecular Dynamics (TAMD), Logarithmic Mean Force Dynamics (LogMFD), and Multiscale Enhanced Sampling (MSES) algorithms. We also discuss combinations with techniques for searching reaction paths. We show the advantages presented by this approach and how it allows to quickly sample important regions of the free-energy landscape via automatic exploration.

Direct and Propagated Effects of Small Molecules on Protein-Protein Interaction Networks

Networks of protein–protein interactions (PPIs) link all aspects of cellular biology. Dysfunction in the assembly or dynamics of PPI networks is a hallmark of human disease, and as such, there is growing interest in the discovery of small molecules that either promote or inhibit PPIs. PPIs were once considered undruggable because of their relatively large buried surface areas and difficult topologies. Despite these challenges, recent advances in chemical screening methodologies, combined with improvements in structural and computational biology have made some of these targets more tractable. In this review, we highlight developments that have opened the door to potent chemical modulators. We focus on how allostery is being used to produce surprisingly robust changes in PPIs, even for the most challenging targets. We also discuss how interfering with one PPI can propagate changes through the broader web of interactions. Through this analysis, it is becoming clear that a combination of direct and propagated effects on PPI networks is ultimately how small molecules re-shape biology.

Structural Model of the Bilitranslocase Transmembrane Domain Supported by NMR and FRET Data

We present a 3D model of the four transmembrane (TM) helical regions of bilitranslocase (BTL), a structurally uncharacterized protein that transports organic anions across the cell membrane. The model was computed by considering helix-helix interactions as primary constraints, using Monte Carlo simulations. The interactions between the TM2 and TM3 segments have been confirmed by Förster resonance energy transfer (FRET) spectroscopy and nuclear magnetic resonance (NMR) spectroscopy, increasing our confidence in the model. Several insights into the BTL transport mechanism were obtained by analyzing the model. For example, the observed cis-trans Leu-Pro peptide bond isomerization in the TM3 fragment may indicate a key conformational change during anion transport by BTL. Our structural model of BTL may facilitate further studies, including drug discovery.

Selectivity and anti-Parkinson's potential of thiadiazolidinone RGS4 inhibitors

Many current therapies target G protein coupled receptors (GPCR), transporters, or ion channels. In addition to directly targeting these proteins, disrupting the protein-protein interactions that localize or regulate their function could enhance selectivity and provide unique pharmacologic actions. Regulators of G protein signaling (RGS) proteins, especially RGS4, play significant roles in epilepsy and Parkinson's disease. Thiadiazolidinone (TDZD) inhibitors of RGS4 are nanomolar potency blockers of the biochemical actions of RGS4 in vitro. Here, we demonstrate the substantial selectivity (8- to >5000-fold) of CCG-203769 for RGS4 over other RGS proteins. It is also 300-fold selective for RGS4 over GSK-3β, another target of this class of chemical scaffolds. It does not inhibit the cysteine protease papain at 100 μM. CCG-203769 enhances Gαq-dependent cellular Ca(2+) signaling in an RGS4-dependent manner. TDZD inhibitors also enhance Gαi-dependent δ-OR inhibition of cAMP production in SH-SY-5Y cells, which express endogenous receptors and RGS4. Importantly, CCG-203769 potentiates the known RGS4 mechanism of Gαi-dependent muscarinic bradycardia in vivo. Furthermore, it reverses raclopride-induced akinesia and bradykinesia in mice, a model of some aspects of the movement disorder in Parkinson's disease. A broad assessment of compound effects revealed minimal off-target effects at concentrations necessary for cellular RGS4 inhibition. These results expand our understanding of the mechanism and specificity of TDZD RGS inhibitors and support the potential for therapeutic targeting of RGS proteins in Parkinson's disease and other neural disorders.

Temperature Accelerated Molecular Dynamics with Soft-Ratcheting Criterion Orients Enhanced Sampling by Low-Resolution Information

Many proteins exhibit an equilibrium between multiple conformations, some of them being characterized only by low-resolution information. Visiting all conformations is a demanding task for computational techniques performing enhanced but unfocused exploration of collective variable (CV) space. Otherwise, pulling a structure toward a target condition biases the exploration in a way difficult to assess. To address this problem, we introduce here the soft-ratcheting temperature-accelerated molecular dynamics (sr-TAMD), where the exploration of CV space by TAMD is coupled to a soft-ratcheting algorithm that filters the evolving CV values according to a predefined criterion. Any low resolution or even qualitative information can be used to orient the exploration. We validate this technique by exploring the conformational space of the inactive state of the catalytic domain of the adenyl cyclase AC from Bordetella pertussis. The domain AC gets activated by association with calmodulin (CaM), and the available crystal structure shows that in the complex the protein has an elongated shape. High-resolution data are not available for the inactive, CaM-free protein state, but hydrodynamic measurements have shown that the inactive AC displays a more globular conformation. Here, using as CVs several geometric centers, we use sr-TAMD to enhance CV space sampling while filtering for CV values that correspond to centers moving close to each other, and we thus rapidly visit regions of conformational space that correspond to globular structures. The set of conformations sampled using sr-TAMD provides the most extensive description of the inactive state of AC up to now, consistent with available experimental information.

Theoretical and computational studies of peptides and receptors of the insulin family

Synergistic interactions among peptides and receptors of the insulin family are required for glucose homeostasis, normal cellular growth and development, proliferation, differentiation and other metabolic processes. The peptides of the insulin family are disulfide-linked single or dual-chain proteins, while receptors are ligand-activated transmembrane glycoproteins of the receptor tyrosine kinase (RTK) superfamily. Binding of ligands to the extracellular domains of receptors is known to initiate signaling via activation of intracellular kinase domains. While the structure of insulin has been known since 1969, recent decades have seen remarkable progress on the structural biology of apo and liganded receptor fragments. Here, we review how this useful structural information (on ligands and receptors) has enabled large-scale atomically-resolved simulations to elucidate the conformational dynamics of these biomolecules. Particularly, applications of molecular dynamics (MD) and Monte Carlo (MC) simulation methods are discussed in various contexts, including studies of isolated ligands, apo-receptors, ligand/receptor complexes and intracellular kinase domains. The review concludes with a brief overview and future outlook for modeling and computational studies in this family of proteins.

Two pathways mediate interdomain allosteric regulation in pin1

Allostery is an essential means for regulating biomolecular functions and provides unique opportunities for drug design, yet our ability to elucidate allosteric mechanisms remains limited. Here, based on extensive molecular dynamics simulations, we present an atomistic picture of the pathways mediating the allosteric regulation of the PPIase domain of Pin1 by its WW domain. Two pathways jointly propagate the action of substrate-WW binding to produce closure and rigidification of three PPIase catalytic-site loops. One pathway preexists in the apo protein, but remains dormant until substrate-WW binding completes the second. The reduction in conformational entropy and preorganization of the catalytic-site loops observed here may explain why substrate-WW binding enhances ligand affinity and catalytic activity of the PPIase domain and suggest a combination drug therapy for Pin1-related diseases. Whereas the traditional view of allostery has emphasized conformational transition, our study identifies a distinct role of conformational dynamics in eliciting allostery.

Temperature-accelerated molecular dynamics gives insights into globular conformations sampled in the free state of the AC catalytic domain

The catalytic domain of the adenyl cyclase (AC) toxin from Bordetella pertussis is activated by interaction with calmodulin (CaM), resulting in cAMP overproduction in the infected cell. In the X-ray crystallographic structure of the complex between AC and the C terminal lobe of CaM, the toxin displays a markedly elongated shape. As for the structure of the isolated protein, experimental results support the hypothesis that more globular conformations are sampled, but information at atomic resolution is still lacking. Here, we use temperature-accelerated molecular dynamics (TAMD) simulations to generate putative all-atom models of globular conformations sampled by CaM-free AC. As collective variables, we use centers of mass coordinates of groups of residues selected from the analysis of standard molecular dynamics (MD) simulations. Results show that TAMD allows extended conformational sampling and generates AC conformations that are more globular than in the complexed state. These structures are then refined via energy minimization and further unrestrained MD simulations to optimize inter-domain packing interactions, thus resulting in the identification of a set of hydrogen bonds present in the globular conformations.

Conformational changes in acetylcholine binding protein investigated by temperature accelerated molecular dynamics

Despite the large number of studies available on nicotinic acetylcholine receptors, a complete account of the mechanistic aspects of their gating transition in response to ligand binding still remains elusive. As a first step toward dissecting the transition mechanism by accelerated sampling techniques, we study the ligand-induced conformational changes of the acetylcholine binding protein (AChBP), a widely accepted model for the full receptor extracellular domain. Using unbiased Molecular Dynamics (MD) and Temperature Accelerated Molecular Dynamics (TAMD) simulations we investigate the AChBP transition between the apo and the agonist-bound state. In long standard MD simulations, both conformations of the native protein are stable, while the agonist-bound structure evolves toward the apo one if the orientation of few key sidechains in the orthosteric cavity is modified. Conversely, TAMD simulations initiated from the native conformations are able to produce the spontaneous transition. With respect to the modified conformations, TAMD accelerates the transition by at least a factor 10. The analysis of some specific residue-residue interactions points out that the transition mechanism is based on the disruption/formation of few key hydrogen bonds. Finally, while early events of ligand dissociation are observed already in standard MD, TAMD accelerates the ligand detachment and, at the highest TAMD effective temperature, it is able to produce a complete dissociation path in one AChBP subunit.

RGS6, but not RGS4, is the dominant regulator of G protein signaling (RGS) modulator of the parasympathetic regulation of mouse heart rate

Parasympathetic activity decreases heart rate (HR) by inhibiting pacemaker cells in the sinoatrial node (SAN). Dysregulation of parasympathetic influence has been linked to sinus node dysfunction and arrhythmia. RGS (regulator of G protein signaling) proteins are negative modulators of the parasympathetic regulation of HR and the prototypical M2 muscarinic receptor (M2R)-dependent signaling pathway in the SAN that involves the muscarinic-gated atrial K(+) channel IKACh. Both RGS4 and RGS6-Gβ5 have been implicated in these processes. Here, we used Rgs4(-/-), Rgs6(-/-), and Rgs4(-/-):Rgs6(-/-) mice to compare the relative influence of RGS4 and RGS6 on parasympathetic regulation of HR and M2R-IKACh-dependent signaling in the SAN. In retrogradely perfused hearts, ablation of RGS6, but not RGS4, correlated with decreased resting HR, increased heart rate variability, and enhanced sensitivity to the negative chronotropic effects of the muscarinic agonist carbachol. Similarly, loss of RGS6, but not RGS4, correlated with enhanced sensitivity of the M2R-IKACh signaling pathway in SAN cells to carbachol and a significant slowing of M2R-IKACh deactivation rate. Surprisingly, concurrent genetic ablation of RGS4 partially rescued some deficits observed in Rgs6(-/-) mice. These findings, together with those from an acute pharmacologic approach in SAN cells from Rgs6(-/-) and Gβ5(-/-) mice, suggest that the partial rescue of phenotypes in Rgs4(-/-):Rgs6(-/-) mice is attributable to another R7 RGS protein whose influence on M2R-IKACh signaling is masked by RGS4. Thus, RGS6-Gβ5, but not RGS4, is the primary RGS modulator of parasympathetic HR regulation and SAN M2R-IKACh signaling in mice.