Metabotropic receptor activation, desensitization and sequestration-I: modelling calcium and inositol 1,4,5-trisphosphate dynamics following receptor activation

Modeling P2Y receptor-Ca2+ response coupling in taste cells

Here we elaborated an analytical approach for the simulation of dose-response curves mediated by cellular receptors coupled to PLC and Ca(2+) mobilization. Based on a mathematical model of purinergic Ca(2+) signaling in taste cells, the analysis of taste cells responsiveness to nucleotides was carried out. Consistently with the expression of P2Y(2) and P2Y(4) receptors in taste cells, saturating ATP and UTP equipotently mobilized intracellular Ca(2+). Cellular responses versus concentration of BzATP, a P2Y(2) agonist and a P2Y(4) antagonist, implicated high and low affinity BzATP receptors. Suramin modified the BzATP dose-response curve in a manner that suggested the low affinity receptor to be weakly sensitive to this P2Y antagonist. Given that solely P2Y(2) and P2Y(11) are BzATP receptors, their high sensitivity to suramin is poorly consistent with the suramin effects on BzATP responses. We simulated a variety of dose-response curves for different P2Y receptor sets and found that the appropriate fit of the overall pharmacological data was achievable only with dimeric receptors modeled as P2Y(2)/P2Y(4) homo- and heterodimers. Our computations and analytical analysis of experimental dose-response curves raise the possibility that ATP responsiveness of mouse taste cells is mediated by P2Y(2) and P2Y(4) receptors operative mostly in the dimeric form.

Inverse drug screens: a rapid and inexpensive method for implicating molecular targets

Identification of gene products that function in some specific process of interest is a common goal in developmental biology. Although use of drug compounds to probe biological systems has a very long history in teratology and toxicology, systematic hierarchical drug screening has not been capitalized upon by the developmental biology community. This "chemical genetics" approach can greatly benefit the study of embryonic and regenerative systems, and we have formalized a strategy for using known pharmacological compounds to implicate specific molecular candidates in any chosen biological phenomenon. Taking advantage of a hierarchical structure that can be imposed on drug reagents in a number of fields such as ion transport, neurotransmitter function, metabolism, and cytoskeleton, any assay can be carried out as a binary search algorithm. This inverse drug screen methodology is much more efficient than exhaustive testing of large numbers of drugs, and reveals the identity of a manageable number of specific molecular candidates that can then be validated and targeted using more expensive and specific molecular reagents. Here, we describe the process of this loss-of-function screen and illustrate its use in uncovering novel bioelectrical and serotonergic mechanisms in embryonic patterning. This technique is an inexpensive and rapid complement to existing molecular screening strategies. Moreover, it is applicable to maternal proteins, and model species in which traditional genetic screens are not feasible, significantly extending the opportunities to identify key endogenous players in biological processes.

Purinergic junctional transmission and propagation of calcium waves in spinal cord astrocyte networks

Micro-photolithographic methods have been employed to form discrete patterns of spinal cord astrocytes that allow quantitative measurements of Ca(2+) wave propagation. Astrocytes were confined to lanes 20-100 microm wide and Ca(2+) waves propagated from a point of mechanical stimulation or of application of adenosine triphosphate; all Ca(2+) wave propagation was blocked by simultaneous application of purinergic P2Y(1) and P2Y(2) antagonists. Stimulation of an astrocyte at one end of a lane, followed by further stimulation of this astrocyte, gave rise to Ca(2+) transients in the same astrocytes; however, if the second stimulation was applied to an astrocyte at the other end of the lane, then this gave rise to a different but overlapping set of astrocytes generating a Ca(2+) signal. Both the amplitude and velocity of the Ca(2+) wave decreased over 270 microm from the point of initiation, and thereafter remained, on average, constant with random variations for at least a further 350 microm. Also, the percentage of astrocytes that gave a Ca(2+) transient decreased with distance along lanes. All the above observations were quantitatively predicted by our recent theoretical model of purinergic junctional transmission, as was the Ca(2+) wave propagation along and between parallel lanes of astrocytes different distances apart. These observations show that a model in which the main determinants are the diffusion of adenosine triphosphates regeneratively released from a stimulated astrocyte, together with differences in the properties and density of the purinergic P2Y receptors on astrocytes, is adequate to predict a wide range of Ca(2+) wave transmission and propagation phenomena.

A quantitative model of purinergic junctional transmission of calcium waves in astrocyte networks

A principal means of transmitting intracellular calcium (Ca2+) waves at junctions between astrocytes involves the release of the chemical transmitter adenosine triphosphate (ATP). A model of this process is presented in which activation of purinergic P2Y receptors by ATP triggers the release of ATP, in an autocrine manner, as well as concomitantly increasing intracellular Ca2+. The dependence of the temporal characteristics of the Ca2+ wave are shown to critically depend on the dissociation constant (K(R)) for ATP binding to the P2Y receptor type. Incorporating this model astrocyte into networks of these cells successfully accounts for many of the properties of propagating Ca2+ waves, such as the dependence of velocity on the type of P2Y receptor and the time-lag of the Ca2+ wave behind the ATP wave. In addition, the conditions under which Ca2+ waves may jump from one set of astrocytes across an astrocyte-free lane to another set of astrocytes are quantitatively accounted for by the model. The properties of purinergic transmission at astrocyte junctions may determine many of the characteristics of Ca2+ propagation in networks of these cells.

Functional selectivity of G protein signaling by agonist peptides and thrombin for the protease-activated receptor-1

Thrombin activates protease-activated receptor-1 (PAR-1) by cleavage of the amino terminus to unmask a tethered ligand. Although peptide analogs can activate PAR-1, we show that the functional responses mediated via PAR-1 differ between the agonists. Thrombin caused endothelial monolayer permeability and mobilized intracellular calcium with EC(50) values of 0.1 and 1.7 nm, respectively. The opposite order of activation was observed for agonist peptide (SFLLRN-CONH(2) or TFLLRNKPDK) activation. The addition of inactivated thrombin did not affect agonist peptide signaling, suggesting that the differences in activation mechanisms are intramolecular in origin. Although activation of PAR-1 or PAR-2 by agonist peptides induced calcium mobilization, only PAR-1 activation affected barrier function. Induced barrier permeability is likely to be Galpha(12/13)-mediated as chelation of Galpha(q)-mediated intracellular calcium with BAPTA-AM, pertussis toxin inhibition of Galpha(i/o), or GM6001 inhibition of matrix metalloproteinase had no effect, whereas Y-27632 inhibition of the Galpha(12/13)-mediated Rho kinase abrogated the response. Similarly, calcium mobilization is Galpha(q)-mediated and independent of Galpha(i/o) and Galpha(12/13) because pertussis toxin Y-27632 and had no effect, whereas U-73122 inhibition of phospholipase C-beta blocked the response. It is therefore likely that changes in permeability reflect Galpha(12/13) activation, and changes in calcium reflect Galpha(q) activation, implying that the pharmacological differences between agonists are likely caused by the ability of the receptor to activate Galpha(12/13) or Galpha(q). This functional selectivity was characterized quantitatively by a mathematical model describing each step leading to Rho activation and/or calcium mobilization. This model provides an estimate that peptide activation alters receptor/G protein binding to favor Galpha(q) activation over Galpha(12/13) by approximately 800-fold.

A signal transduction pathway model prototype I: From agonist to cellular endpoint

The postgenomic era is providing a wealth of information about the genes involved in many cellular processes. However, the ability to apply this information to understanding cellular signal transduction is limited by the lack of tools that quantitatively describe cellular signaling processes. The objective of the current studies is to provide a framework for modeling cellular signaling processes beginning at a plasma membrane receptor and ending with a measurable endpoint in the signaling process. Agonist-induced Ca(2+) mobilization coupled to down stream phosphorylation events was modeled using knowledge of in vitro and in vivo process parameters. The simulation process includes several modules that describe cellular processes involving receptor activation phosphoinositide metabolism, Ca(2+)-release, and activation of a calmodulin-dependent protein kinase. A Virtual Cell-based simulation was formulated using available literature data and compared to new and existing experimental results. The model provides a new approach to facilitate hypothesis-driven investigation and experimental design based upon simulation results. These investigations may be directed at the timing of multiple phosphorylation/dephosphorylation events affecting key enzymatic activities in the signaling pathway.

Calcium mobilization and spontaneous transient outward current characteristics upon agonist activation of P2Y2 receptors in smooth muscle cells

A quantitative model is provided that links the process of metabotropic receptor activation and sequestration to the generation of inositol 1,4,5-trisphosphate, the subsequent release of calcium from the central sarcoplasmic reticulum, and the consequent release of calcium from subsarcolemma sarcoplasmic reticulum that acts on large-conductance potassium channels to generate spontaneous transient outward currents (STOCs). This model is applied to the case of STOC generation in vascular A7r5 smooth muscle cells that have been transfected with a chimera of the P2Y(2) metabotropic receptor and green fluorescent protein (P2Y(2)-GFP) and exposed to the P2Y(2) receptor agonist uridine 5'-triphosphate. The extent of P2Y(2)-GFP sequestration from the membrane on exposure to uridine 5'-triphosphate, the ensuing changes in cytosolic calcium concentration, as well as the interval between STOCs that are subsequently generated, are used to determine parameter values in the model. With these values, the model gives a good quantitative prediction of the dynamic changes in STOC amplitude observed upon activation of metabotropic P2Y(2) receptors in the vascular smooth muscle cell line.

A quantitative description of the contraction of blood vessels following the release of noradrenaline from sympathetic varicosities

A model is presented that highlights the principal factors determining the form and extent of contraction in arteries upon stimulation of their sympathetic nerve supply. This model incorporates a previous quantitative model of the process of noradrenaline (NAd) diffusion into the vascular media and reuptake into sympathetic varicosities during nerve stimulation (J. Theor. Biol. 226 (2004) 359). It is also dependent on a model of how the subsequent activation of metabotropic receptors initiates a G-protein cascade, resulting in the production of inositol trisphosphate (IP3) and an increase in intracellular calcium concentration, [Ca2+]i, in the smooth muscle cells (J. Theor. Biol. 223 (2003) 93). In the present work we couple this rise in [Ca2+]i to the increase in phosphorylated myosin bound to actin in the cells and hence determine the force development in arteries due to nerve stimulation. The model accounts for force development as a function of [Ca2+]i and for the rate of change of force as a function of the rate of change of [Ca2+]i in single smooth muscle cells. It also accounts for the characteristic time course of the force developed by the media of the rat-tail artery upon nerve stimulation. This consists of a rapid rise to a transient peak followed by a sustained plateau of contraction during the stimulation period, after which the contraction slowly decays back to baseline at a rate dependent on the strength of the stimulation. The model indicates that the transient peak is primarily due to the partial block of the IP3 receptor by the rise in [Ca2+]i and that the main determinant of the equilibrium condition indicated by the plateau phase is the rate of pumping of calcium into the sarcoplasmic reticulum. The relatively slow decline of contraction at the end of nerve stimulation is primarily a consequence of the slow rates of removal of NAd from the media by diffusion and reuptake into the sympathetic varicosities. The model thus provides a quantitative account of vascular smooth muscle contraction upon sympathetic nerve stimulation.

Metabotropic receptor activation, desensitization and sequestration-II: modelling the dynamics of the pleckstrin homology domain

Recent observations have been made regarding the generation of inositol 1,4,5-trisphosphate (IP(3)), using chimeras of green fluorescent protein and the pleckstrin homology domain of phospholipase C-delta. In this paper a model is presented giving the quantitative relations between the green fluorescent protein-pleckstrin homology domain (GFP-PHD) construct and membrane phosphatidylinositol 4,5-bisphosphate (PIP(2)) levels as well as the concentration of IP(3), the product of hydrolysis of PIP(2). The model can correctly reproduce the dependence of cytosolic GFP-PHD fluorescence on IP(3) concentration. This model extends a previous one (Metabotropic receptor activation, desensitization and sequestration-I: modelling calcium and inositol 1,4,5-trisphosphate dynamics following receptor activation, in this issue) dealing with the processes governing the production of IP(3) and the subsequent calcium (Ca2+) changes in cells following activation of metabotropic receptors. This model is applied to the case of purinergic P(2)Y(2) receptor activation in Madin-Darby Canine Kidney (MDCK) cells with adenosine triphosphate (ATP) (Science 284 (1999) 1527). It is shown that it can correctly reproduce the dependence of GFP-PHD fluorescence on the concentration of P(2)Y(2) receptor ligand, as well as the temporal changes of GFP-PHD fluorescence following application of ligand.