Bioactive Properties of Venoms Isolated from Whiptail Stingrays and the Search for Molecular Mechanisms and Targets.

The venom-containing barb attached to their 'whip-like' tail provides stingrays a defensive mechanism for evading predators such as sharks. From human encounters, dermal stingray envenomation is characterized by intense pain often followed by tissue necrosis occurring over several days to weeks. The bioactive components in stingray venoms (SRVs) and their molecular targets and mechanisms that mediate these complex responses are not well understood. Given the utility of venom-derived proteins from other venomous species for biomedical and pharmaceutical applications, we set out to characterize the bioactivity of SRV extracts from three local species that belong to the Dasyatoidea 'whiptail' superfamily. Multiple cell-based assays were used to quantify and compare the in vitro effects of these SRVs on different cell lines. All three SRVs demonstrated concentration-dependent growth-inhibitory effects on three different human cell lines tested. In contrast, a mouse fibrosarcoma cell line was markedly resistant to all three SRVs, indicating the molecular target(s) for mediating the SRV effects are not expressed on these cells. The multifunctional SRV responses were characterized by an acute disruption of cell adhesion leading to apoptosis. These findings aim to guide future investigations of individual SRV proteins and their molecular targets for potential use in biomedical applications.

Structural Determinants Mediating Tertiapin Block of Neuronal Kir3.2 Channels.

Tertiapin (TPN) is a 21 amino acid venom peptide from Apis mellifera that inhibits certain members of the inward rectifier potassium (Kir) channel family at a nanomolar affinity with limited specificity. Structure-based computational simulations predict that TPN behaves as a pore blocker; however, the molecular determinants mediating block of neuronal Kir3 channels have been inconclusive and unvalidated. Here, using molecular docking and molecular dynamics (MD) simulations with 'potential of mean force' (PMF) calculations, we investigated the energetically most favored interaction of TPN with several Kir3.x channel structures. The resulting binding model for Kir3.2-TPN complexes was then tested by targeted mutagenesis of the predicted contact sites, and their impact on the functional channel block was measured electrophysiologically. Together, our findings indicate that a high-affinity TPN block of Kir3.2 channels involves a pore-inserting lysine side chain requiring (1) hydrophobic interactions at a phenylalanine ring surrounding the channel pore and (2) electrostatic interactions with two adjacent Kir3.2 turret regions. Together, these interactions collectively stabilize high-affinity toxin binding to the Kir3.2 outer vestibule, which orients the ε-amino group of TPN-K21 to occupy the outermost K+ binding site of the selectivity filter. The structural determinants for the TPN block described here also revealed a favored subunit arrangement for assembled Kir3.x heteromeric channels, in addition to a multimodal binding capacity of TPN variants consistent with the functional dyad model for polybasic peptide pore blockers. These novel findings will aid efforts in re-engineering the TPN pharmacophore to develop peptide variants having unique and distinct Kir channel blocking properties.

Identification of Aethina tumida Kir Channels as Putative Targets of the Bee Venom Peptide Tertiapin Using Structure-Based Virtual Screening Methods.

Venoms are comprised of diverse mixtures of proteins, peptides, and small molecules. Identifying individual venom components and their target(s) with mechanism of action is now attainable to understand comprehensively the effectiveness of venom cocktails and how they collectively function in the defense and predation of an organism. Here, structure-based computational methods were used with bioinformatics tools to screen and identify potential biological targets of tertiapin (TPN), a venom peptide from Apis mellifera (European honey bee). The small hive beetle (Aethina tumida (A. tumida)) is a natural predator of the honey bee colony and was found to possess multiple inwardly rectifying K+ (Kir) channel subunit genes from a genomic BLAST search analysis. Structure-based virtual screening of homology modelled A. tumida Kir (atKir) channels found TPN to interact with a docking profile and interface "footprint" equivalent to known TPN-sensitive mammalian Kir channels. The results support the hypothesis that atKir channels, and perhaps other insect Kir channels, are natural biological targets of TPN that help defend the bee colony from infestations by blocking K+ transport via atKir channels. From these in silico findings, this hypothesis can now be subsequently tested in vitro by validating atKir channel block as well as in vivo TPN toxicity towards A. tumida. This study highlights the utility and potential benefits of screening in virtual space for venom peptide interactions and their biological targets, which otherwise would not be feasible.

Venom-derived peptides inhibiting Kir channels: Past, present, and future.

Inwardly rectifying K+ (Kir) channels play a significant role in vertebrate and invertebrate biology by regulating the movement of K+ ions involved in membrane transport and excitability. Yet unlike other ion channels including their ancestral K+-selective homologs, there are very few venom toxins known to target and inhibit Kir channels with the potency and selectivity found for the Ca2+-activated and voltage-gated K+ channel families. It is unclear whether this is simply due to a lack of discovery, or instead a consequence of the evolutionary processes that drive the development of venom components towards their targets based on a collective efficacy to 1) elicit pain for defensive purposes, 2) promote paralysis for prey capture, or 3) facilitate delivery of venom components into the circulation. The past two decades of venom screening has yielded three venom peptides with inhibitory activity towards mammalian Kir channels, including the discovery of tertiapin, a high-affinity pore blocker from the venom of the European honey bee Apis mellifera. Venomics and structure-based computational approaches represent exciting new frontiers for venom peptide development, where re-engineering peptide 'scaffolds' such as tertiapin may aid in the quest to expand the palette of potent and selective Kir channel blockers for future research and potentially new therapeutics. This article is part of the Special Issue entitled 'Venom-derived Peptides as Pharmacological Tools.'

RGS Redundancy and Implications in GPCR-GIRK Signaling.

Regulators of G protein signaling (RGS proteins) are key components of GPCR complexes, interacting directly with G protein α-subunits to enhance their intrinsic GTPase activity. The functional consequence is an accelerated termination of G protein effectors including certain ion channels. RGS proteins have a profound impact on the membrane-delimited gating behavior of G-protein-activated inwardly rectifying K(+) (GIRK) channels as demonstrated in reconstitution assays and recent RGS knockout mice studies. Akin to GPCRs and G protein αßγ subunits, multiple RGS isoforms are expressed within single GIRK-expressing neurons, suggesting functional redundancy and/or specificity in GPCR-GIRK channel signaling. The extent and impact of RGS redundancy in neuronal GPCR-GIRK channel signaling is currently not fully appreciated; however, recent studies from RGS knockout mice are providing important new clues on the impact of individual endogenous RGS proteins and the extent of RGS functional redundancy. Incorporating "tools" such as engineered RGS-resistant Gαi/o subunits provide an important assessment method for determining the impact of all endogenous RGS proteins on a given GPCR response and an accounting benchmark to assess the impact of individual RGS knockouts on overall RGS redundancy within a given neuron. Elucidating the degree of regulation attributable to specific RGS proteins in GIRK channel function will aid in the assessment of individual RGS proteins as viable therapeutic targets in epilepsy, ataxia's, memory disorders, and a growing list of neurological disorders.

A computational design approach for virtual screening of peptide interactions across K(+) channel families.

Ion channels represent a large family of membrane proteins with many being well established targets in pharmacotherapy. The 'druggability' of heteromeric channels comprised of different subunits remains obscure, due largely to a lack of channel-specific probes necessary to delineate their therapeutic potential in vivo. Our initial studies reported here, investigated the family of inwardly rectifying potassium (Kir) channels given the availability of high resolution crystal structures for the eukaryotic constitutively active Kir2.2 channel. We describe a 'limited' homology modeling approach that can yield chimeric Kir channels having an outer vestibule structure representing nearly any known vertebrate or invertebrate channel. These computationally-derived channel structures were tested ""in silico for 'docking' to NMR structures of tertiapin (TPN), a 21 amino acid peptide found in bee venom. TPN is a highly selective and potent blocker for the epithelial rat Kir1.1 channel, but does not block human or zebrafish Kir1.1 channel isoforms. Our Kir1.1 channel-TPN docking experiments recapitulated published in vitro ""findings for TPN-sensitive and TPN-insensitive channels. Additionally, in silico site-directed mutagenesis identified 'hot spots' within the channel outer vestibule that mediate energetically favorable docking scores and correlate with sites previously identified with in vitro thermodynamic mutant-cycle analysis. These 'proof-of-principle' results establish a framework for virtual screening of re-engineered peptide toxins for interactions with computationally derived Kir channels that currently lack channel-specific blockers. When coupled with electrophysiological validation, this virtual screening approach may accelerate the drug discovery process, and can be readily applied to other ion channels families where high resolution structures are available.

GPCR-Kir channel signaling complexes: defining rules of engagement.

Ion channels and G protein-coupled receptors (GPCRs) are integral transmembrane proteins vital to a multitude of cell signaling and physiological functions. Members of these large protein families are known to interact directly with various intracellular protein partners in a dynamic and isoform-dependent manner, ultimately shaping their life cycle and signal output. The family of G protein-gated inwardly rectifying potassium channels (Kir3 or GIRK) expressed in brain, heart, and endocrine tissues were recently shown to stably associate with several different GPCRs, forming the basis of a macromolecular ion channel-GPCR signaling complex. The molecular determinants that mediate and maintain GPCR-Kir3 channel complexes are currently not well understood. Recent findings and emerging hypotheses on the assembly and stability of multiprotein GPCR-Kir channel signaling complexes are discussed, highlighting distinct mechanisms used by different Kir channel families. These protein-protein interaction processes are crucial in determining both the synaptic response times and the extent of GPCR "cross-talk" in Kir3-mediated inhibitory synaptic transmission.

RGS3 and RGS4 differentially associate with G protein-coupled receptor-Kir3 channel signaling complexes revealing two modes of RGS modulation. Precoupling and collision coupling.

RGS3 and RGS4 are GTPase-activating proteins expressed in the brain and heart that accelerate the termination of G(i/o)- and G(q)-mediated signaling. We report here the determinants mediating selective association of RGS4 with several G protein-coupled receptors (GPCRs) that form macromolecular complexes with neuronal G protein-gated inwardly rectifying potassium (Kir3 or GIRK) channels. Kir3 channels are instrumental in regulating neuronal firing in the central and peripheral nervous system and pacemaker activity in the heart. By using an epitope-tagged degradation-resistant RGS4 mutant, RGS4(C2V), immunoprecipitation of several hemagglutinin-tagged G(i/o)-coupled and G(q)-coupled receptors expressed in Chinese hamster ovary (CHO-K1) cells readily co-precipitated both Kir3.1/Kir3.2a channels and RGS4(C2V). In contrast to RGS4(C2V), the closely related and functionally active RGS3 "short" isoform (RGS3s) did not interact with any of the GPCR-Kir3 channel complexes examined. Deletion and chimeric RGS constructs indicate both the N-terminal domain and the RGS domain of RGS4(C2V) are necessary for association with m2 receptor-Kir3.1/Kir3.2a channel complexes, where the GPCR was found to be the major target for RGS4(C2V) interaction. The functional impact of RGS4(C2V) "precoupling" to the GPCR-Kir3 channel complex versus RGS3s "collision coupling" was a 100-fold greater potency in the acceleration of G protein-dependent Kir3 channel-gating kinetics with no attenuation in current amplitude. These findings demonstrate that RGS4, a highly regulated modulator and susceptibility gene for schizophrenia, can directly associate with multiple GPCR-Kir3 channel complexes and may affect a wide range of neurotransmitter-mediated inhibitory and excitatory events in the nervous and cardiovascular systems.

R7BP augments the function of RGS7*Gbeta5 complexes by a plasma membrane-targeting mechanism.

The RGS7 (R7) family of G protein regulators, Gbeta5, and R7BP form heterotrimeric complexes that potently regulate the kinetics of G protein-coupled receptor signaling. Reversible palmitoylation of R7BP regulates plasma membrane/nuclear shuttling of R7*Gbeta5*R7BP heterotrimers. Here we have investigated mechanisms whereby R7BP controls the function of the R7 family. We show that unpalmitoylated R7BP undergoes nuclear/cytoplasmic shuttling and that a C-terminal polybasic motif proximal to the palmitoylation acceptor sites of R7BP mediates nuclear localization, palmitoylation, and plasma membrane targeting. These results suggest a novel mechanism whereby palmitoyltransferases and nuclear import receptors both utilize the C-terminal domain of R7BP to determine the trafficking fate of R7*Gbeta5*R7BP heterotrimers. Analogous mechanisms may regulate other signaling proteins whose distribution between the plasma membrane and nucleus is controlled by palmitoylation. Lastly, we show that cytoplasmic RGS7*Gbeta5*R7BP heterotrimers and RGS7*Gbeta5 heterodimers are equivalently inefficient regulators of G protein-coupled receptor signaling relative to plasma membrane-bound heterotrimers bearing palmitoylated R7BP. Therefore, R7BP augments the function of the complex by a palmitoylation-regulated plasma membrane-targeting mechanism.

Neuronal Kir3.1/Kir3.2a channels coupled to serotonin 1A and muscarinic m2 receptors are differentially modulated by the "short" RGS3 isoform.

"Regulators of G protein signaling" (RGS proteins) have profound effects on ion channels regulated by G protein-coupled receptor (GPCR) signaling, including the G protein-gated inwardly rectifying K+ (GIRK) channels that inhibit excitability of neuronal, endocrine, and cardiac cells. Here we describe the effects of an alternatively spliced "short" RGS3 isoform (RGS3s) in comparison to RGS4, on the temporal and steady-state gating properties of neuronal GIRK channels (Kir3.1/Kir3.2a) activated by either serotonin 1A (5-HT(1A)) receptors or muscarinic m2 receptors expressed in Chinese hamster ovary (CHO-K1) cells. RGS3s is abundantly expressed in brain and contains a unique short N-terminus via alternative splicing that distinguishes it from other RGS3 isoforms as well as other members of the B/R4 RGS gene subfamily. Our results indicate that RGS3s and RGS4 similarly affect the temporal and steady-state gating properties of 5-HT(1A) receptor-coupled Kir3.1/Kir3.2a channels, but differentially modulate muscarinic m2 receptor-coupled channels. RGS3s caused a significant approximately 45% reduction in the maximal acetylcholine (ACh)-evoked GIRK current amplitude and a marked shift in the steady-state ACh dose-response relation indicative of a reduction in functionally coupled m2 receptor-GIRK channel complexes. Yet RGS3s still accelerated the m2 receptor-dependent GIRK activation, deactivation, and acute desensitization time course consistent with the RGS-enhanced GAP activity that was also observed with RGS4. Several mechanisms that may contribute to the receptor-dependent effects of RGS3s are discussed with particular attention to the role of the distinct N-terminal domain. Our findings highlight the potential impact of selective RGS-GPCR interactions on neuronal GIRK channel function that may affect the properties of inhibitory postsynaptic potentials activated by different GPCR-GIRK channel complexes.

Palmitoylation regulates plasma membrane-nuclear shuttling of R7BP, a novel membrane anchor for the RGS7 family.

The RGS7 (R7) family of RGS proteins bound to the divergent Gbeta subunit Gbeta5 is a crucial regulator of G protein-coupled receptor (GPCR) signaling in the visual and nervous systems. Here, we identify R7BP, a novel neuronally expressed protein that binds R7-Gbeta5 complexes and shuttles them between the plasma membrane and nucleus. Regional expression of R7BP, Gbeta5, and R7 isoforms in brain is highly coincident. R7BP is palmitoylated near its COOH terminus, which targets the protein to the plasma membrane. Depalmitoylation of R7BP translocates R7BP-R7-Gbeta5 complexes from the plasma membrane to the nucleus. Compared with nonpalmitoylated R7BP, palmitoylated R7BP greatly augments the ability of RGS7 to attenuate GPCR-mediated G protein-regulated inward rectifying potassium channel activation. Thus, by controlling plasma membrane nuclear-shuttling of R7BP-R7-Gbeta5 complexes, reversible palmitoylation of R7BP provides a novel mechanism that regulates GPCR signaling and potentially transduces signals directly from the plasma membrane to the nucleus.

Measuring the modulatory effects of RGS proteins on GIRK channels.

Discovery of "regulators of G-protein signaling" (RGS) as GTPase-activating proteins for heterotrimeric G proteins has provided a highly sought "missing link," reconciling past discrepancies between the in vitro GTPase activity of purified G proteins and the kinetics of physiological responses mediated by G-protein signaling in vivo. With the number of RGS genes in the mammalian genome at more than 30, associating specific RGS proteins to specific G-protein-coupled receptor (GPCR) signaling events has become a focus of RGS investigators. The ubiquitous expression of multiple RGS proteins has complicated this effort, yet the outlook has been encouraged with the identification of RGS9 as the determinant mediating rapid recovery of the transducin-dependent photoresponse. G-protein-gated inwardly rectifying potassium (GIRK) channels that mediate inhibitory synaptic transmission via GPCR activation of pertussis toxin-sensitive G proteins are similarly accelerated by RGS proteins when reconstituted in heterologous cell expression systems and fully reproduce the gating properties of native GIRK channels in neurons and cardiomyocytes. The endogenous neuronal and cardiac RGS protein(s) that accelerate GPCR-->GIRK channel-gating kinetics are currently not known. This article describes methods used to measure the receptor-dependent GIRK channel-gating parameters reconstituted in Chinese hamster ovary (CHO-K1) cells and Xenopus oocytes, as well as rat atrial myocytes and rat cerebellar granule neurons as model cells with native GPCR-->GIRK channel signaling. Applications of these methods for structure-function-based studies of RGS proteins, G proteins, and GPCRs are discussed. We also describe single cell reverse transcriptase polymerase chain reaction methods developed to profile atrial myocyte and neuronal RGS expression to identify specific RGS proteins for targeted knockdown or knockout.

Gbetagamma-activated inwardly rectifying K(+) (GIRK) channel activation kinetics via Galphai and Galphao-coupled receptors are determined by Galpha-specific interdomain interactions that affect GDP release rates.

Gbetagamma-activated inwardly rectifying K(+) (GIRK) channels have distinct gating properties when activated by receptors coupled specifically to Galpha(o) versus Galpha(i) subunit isoforms, with Galpha(o)-coupled currents having approximately 3-fold faster agonist-evoked activation kinetics. To identify the molecular determinants in Galpha subunits mediating these kinetic differences, chimeras were constructed using pertussis toxin (PTX)-insensitive Galpha(oA) and Galpha(i2) mutant subunits (Galpha(oA(C351G)) and Galpha(i2(C352G))) and examined in PTX-treated Xenopus oocytes expressing muscarinic m2 receptors and Kir3.1/3.2a channels. These experiments revealed that the alpha-helical N-terminal region (amino acids 1-161) and the switch regions of Galpha(i2) (amino acids 162-262) both partially contribute to slowing the GIRK activation time course when compared with the Galpha(oA(C351G))-coupled response. When present together, they fully reproduce Galpha(i2(C352G))-coupled GIRK kinetics. The Galpha(i2) C-terminal region (amino acids 263-355) had no significant effect on GIRK kinetics. Complementary responses were observed with chimeras substituting the Galpha(o) switch regions into the Galpha(i2(C352G)) subunit, which partially accelerated the GIRK activation rate. The Galpha(oA)/Galpha(i2) chimera results led us to examine an interaction between the alpha-helical domain and the Ras-like domain previously implicated in mediating a 4-fold slower in vitro basal GDP release rate in Galpha(i1) compared with Galpha(o). Mutations disrupting the interdomain contact in Galpha(i2(C352G)) at either the alphaD-alphaE loop (R145A) or the switch III loop (L233Q/A236H/E240T/M241T), significantly accelerated the GIRK activation kinetics consistent with the Galpha(i2) interdomain interface regulating receptor-catalyzed GDP release rates in vivo. We propose that differences in Galpha(i) versus Galpha(o)-coupled GIRK activation kinetics are due to intrinsic differences in receptor-catalyzed GDP release that rate-limit Gbetagamma production and is attributed to heterogeneity in Galpha(i) and Galpha(o) interdomain contacts.

Gating properties of GIRK channels activated by Galpha(o)- and Galpha(i)-coupled muscarinic m2 receptors in Xenopus oocytes: the role of receptor precoupling in RGS modulation.

'Regulators of G protein Signalling' (RGSs) accelerate the activation and deactivation kinetics of G protein-gated inwardly rectifying K(+) (GIRK) channels. In an apparent paradox, RGSs do not reduce steady-state GIRK current amplitudes as expected from the accelerated rate of deactivation when reconstituted in Xenopus oocytes. We present evidence here that this kinetic anomaly is dependent on the degree of G protein-coupled receptor (GPCR) precoupling, which varies with different Galpha(i/o)-RGS complexes. The gating properties of GIRK channels (Kir3.1/Kir3.2a) activated by muscarinic m2 receptors at varying levels of G protein expression were examined with or without the co-expression of either RGS4 or RGS7 in Xenopus oocytes. Different levels of specific m2 receptor-Galpha coupling were established by uncoupling endogenous pertussis toxin (PTX)-sensitive Galpha(i/o) subunits with PTX, while expressing varying amounts of a single PTX-insensitive subunit (Galpha(i1(C351G)), Galpha(i2(C352G)), Galpha(i3(C351G)), Galpha(oA(C351G)), or Galpha(oB(C351G))). Co-expression of each of the PTX-insensitive Galpha(i/o) subunits rescued acetylcholine (ACh)-elicited GIRK currents (I(K,ACh)) in a concentration-dependent manner, with Galpha(o) isoforms being more effective than Galpha(i) isoforms. Receptor-independent 'basal' GIRK currents (I(K,basal)) were reduced with increasing expression of PTX-insensitive Galpha subunits and were accompanied by a parallel rise in I(K,ACh). These effects together are indicative of increased Gbetagamma scavenging by the expressed Galpha subunit and the subsequent formation of functionally coupled m2 receptor-G protein heterotrimers (Galpha((GDP))betagamma). Co-expression of RGS4 accelerated all the PTX-insensitive Galpha(i/o)-coupled GIRK currents to a similar extent, yet reduced I(K,ACh) amplitudes 60-90 % under conditions of low Galpha(i/o) coupling. Kinetic analysis indicated the RGS4-dependent reduction in steady-state GIRK current was fully explained by the accelerated deactivation rate. Thus kinetic inconsistencies associated with RGS4-accelerated GIRK currents occur at a critical threshold of G protein coupling. In contrast to RGS4, RGS7 selectively accelerated Galpha(o)-coupled GIRK currents. Co-expression of Gbeta5, in addition to enhancing the kinetic effects of RGS7, caused a significant reduction (70-85 %) in steady-state GIRK currents indicating RGS7-Gbeta5 complexes disrupt Galpha(o) coupling. Altogether these results provide further evidence for a GPCR-Galphabetagamma-GIRK signalling complex that is revealed by the modulatory affects of RGS proteins on GIRK channel gating. Our functional experiments demonstrate that the formation of this signalling complex is markedly dependent on the concentration and composition of G protein-RGS complexes.