Stimulus amplification, efficacy, and the operational model. Part II--ternary complex occupancy mechanisms.

To understand how deviations from simple binary occupancy affect measures of efficacy, the generalized stimulus function developed in Part I was used to examine the actions of drugs in systems where occupation of the receptor was modeled using a two-state, ternary complex, or combination of mechanisms.Amplification of drug responses can occur during formation of an active agonist-receptor complex, during generation of the initial stimulus, and during signal transduction. Expressions were derived to characterize the separate contributions of these three phases. Ideally, comparison of relative intrinsic efficacy measures differences in the ability of the agonists to convert active complex into an active stimulus. In practice, differences in the ability of the drugs to form the stimulus-generating complex may also contribute to the efficacy ratio and must be taken into consideration. Failure to adequately account for differences in occupancy can result in overestimation of the efficacy ratio. The magnitude of the difference between true and experimental measures of intrinsic efficacy may be affected by G protein concentration, by the affinity between the G protein and receptor, and (in some models) by the receptor activation constant. Provided that the dissociation constant between the G protein and receptor is of the same order of magnitude as, or lower than the receptor concentration, however, experimental estimates should provide reasonably accurate estimates of the true efficacy ratio.In agreement with previously published experimental data, total G protein level was found capable of influencing agonist maximal response, Emax, and EC50 values in all four ternary complex models. The magnitude of the changes in Emax and EC50 appear to be dependent upon the efficacy of the agonist as well as characteristics of the post-receptor stimulus sequence.Additionally, the concentration-response relations for all four ternary complex models could be reduced to a modified operational format in which the apparent dissociation constant Kapp replaced the true KA, and an apparent operational efficacy, tauapp, replaced tau. tauapp can be estimated experimentally from measurements of the Kapp and EC50, while the operational maximum, Em, may be found from the calculated tauapp and the measured Emax of the response curve. These findings support the use of direct operational model-fitting in a variety of systems, regardless of the mechanisms underlying occupancy.Values of Kapp calculated using the exact formula for [ARG] displayed an anomalous rise or discontinuity where the concentrations of total G protein equaled that of the receptor protein. This discontinuity is not observed in the estimates based on approximations to [ARG], and may explain practical difficulties in evaluating the dissociation constant under these conditions.     

Inverse, protean, and ligand-selective agonism: matters of receptor conformation.

Concepts regarding the mechanisms by which drugs activate receptors to produce physiological response have progressed beyond considering the receptor as a simple on-off switch. Current evidence suggests that the idea that agonists produce only varying degrees of receptor activation is obsolete and must be reconciled with data to show that agonist efficacy has texture as well as magnitude. Thus, agonists can block system constitutive response (inverse agonists), behave as positive and inverse agonists on the same receptor (protean agonists), and differ in the stimulus pattern they produce in physiological systems (ligand-selective agonists). The molecular mechanism for this seemingly diverse array of activities is the same, namely, the selective microaffinity of ligands for different conformational states of the receptor. This paper reviews evidence for the existence of the various types of agonism and the potential therapeutic utility of different agonist types.-Kenakin, T. Inverse, protean, and ligand-selective agonism: matters of receptor conformation.     

Reduction of beta-adrenoceptor function by oxidative stress in the heart

The effect of oxidative stress on beta-adrenoceptor function in the heart was determined. To this end ventricle membranes, field-stimulated rat left atria and field-stimulated rat right ventricle strips were exposed to 0.1 mM cumene hydroperoxide for 20 min. It was found that oxidative stress increased beta-adrenoceptor number and reduced c-AMP formation in the ventricle membranes. In the rat left atria and rat right ventricle strips the efficacy of beta-adrenoceptor agonists was reduced to approximately 30% of the control value, whereas maximal beta-adrenoceptor-mediated response was reduced to 50%. Using membranes from control atria and from atria exposed to oxidative stress, it was found that oxidative stress had no effect on beta-adrenoceptor density, nor on the affinity of (-)isoproterenol for the receptor. c-AMP production in membranes prepared from atria exposed to oxidative stress was reduced to approximately 30% of the c-AMP production in membranes prepared of control atria. In addition, it was found that the shape of the function that transduces the stimulus which is generated by receptor activation into an effect, is not altered by oxidative stress. It was concluded that the reduction of the efficacy of beta-adrenoceptor agonists by oxidative stress is probably caused by the reduction of c-AMP formation. Because the efficacy of forskolin and of dibutyryl c-AMP was not affected by oxidative stress, the reduced c-AMP formation is probably caused by an impaired coupling between the receptor and adenylate cyclase. The reduction of maximal beta-adrenoceptor-mediated response might be the result of cytotoxic aldehydes that are produced during oxidative stress. In ischemia, catecholamine release and subsequent beta-adrenoceptor hyperstimulation lead to cardiotoxicity. As shown in the present study, oxidative stress reduces beta-adrenoceptor function. This might represent a protective physiological feedback mechanism that protects the heart against excessive beta-adrenoceptor stimulation.     

Coding "what" and "when" in the Archer fish retina

Traditionally, the information content of the neural response is quantified using statistics of the responses relative to stimulus onset time with the assumption that the brain uses onset time to infer stimulus identity. However, stimulus onset time must also be estimated by the brain, making the utility of such an approach questionable. How can stimulus onset be estimated from the neural responses with sufficient accuracy to ensure reliable stimulus identification? We address this question using the framework of colour coding by the archer fish retinal ganglion cell. We found that stimulus identity, “what”, can be estimated from the responses of best single cells with an accuracy comparable to that of the animal''s psychophysical estimation. However, to extract this information, an accurate estimation of stimulus onset is essential. We show that stimulus onset time, “when”, can be estimated using a linear-nonlinear readout mechanism that requires the response of a population of 100 cells. Thus, stimulus onset time can be estimated using a relatively simple readout. However, large nerve cell populations are required to achieve sufficient accuracy.

Authors Summary

In our interaction with the environment we are flooded with a stream of numerous objects and events. Our brain needs to understand the nature of these complex and rich stimuli in order to react. Research has shown ways in which a ‘what’ stimulus was presented can be encoded by the neural responses. However, to understand ‘what was the nature of the stimulus’ the brain needs to know ‘when’ the stimulus was presented. Here, we investigated how the onset of visual stimulus can be signalled by the retina to higher brain regions. We used archer fish as a framework to test the notion that the answer to the question of ‘when’ something has been presented lies within the larger cell population, whereas the answer to the question of ‘what’ has been presented may be found at the single-neuron level. The utility of the archer fish as model animal stems from its remarkable ability to shoot down insects settling on the foliage above the water level, and its ability to distinguish between artificial targets. Thus, the archer fish can provide the fish equivalent of a monkey or a human that can report psychophysical decisions.     

Quantification and statistical verification of neuronal stimulus responses from noisy spike train data

Usually neuronal responses to short-lasting stimuli are displayed as peri-stimulus time histogram. The function estimated by such a histogram allows to obtain informations about stimulus-induced postsynaptic events as long as the interpretation is restricted to the first response component after the stimulus. The interpretation of secondary response components is much more difficult, as they may be either due to stimulus effects or represent an echo of the primary response. In the present paper two output functions are developed that do not show such an echoing of responses. The first one, the interspike interval change function, represents an ideal way to quantify a neuronal stimulus response as its amplitude was found to be almost independent of the stimulation strategy used during acquisition of the spike train data. The other function, the displaced impulses function, allows to verify the statistical significance of an observed response component. Both functions may be estimated from stimulus-correlated spike train data, even if the neuron under investigation shows considerable interspike-interval variability in the absence of stimulation. The concepts underlying these neuronal output functions are developed on simulated responses of a Hodgkin-Huxley-type model for a mammalian neuron at body temperature that is exposed to a transient excitatory conductance increase. Additionally, estimation of these output functions is also demonstrated on responses of human soleus motoneurons that were exposed to electrical stimuli of the tibial nerve in the popliteal fossa.     

Interpretation of relative potencies, relative efficacies and apparent affinity constants of agonist drugs estimated from concentration-response curves

Differences in the relative potencies of agonists have been used successfully in the past to classify receptors. Such use of agonists can be justified on the basis of ideas and equations developed using the occupancy model of drug action. However the occupancy model makes no allowance for possible complications which may arise when the drug-receptor complex interacts with a transducer-effector system. For some receptor-effector systems use of an equilibrium ternary complex model may be better than use of the occupancy model but the former still does not take into account the possible effect of guanosine-5'-triphosphate on the system. A steady-state version of the ternary complex model has therefore been analysed to explore possible interpretations of relative potencies, relative efficacies and apparent affinity constants estimated from concentration-response curves. It is concluded that for agonists which act on receptors which function through G-proteins these pharmacological parameters may depend on the concentration of the relevant G-protein in the cell membranes and on the intracellular concentrations of guanosine-5'-triphosphate and guanosine-5'-diphosphate. If these concentrations vary appreciably between tissues then the parameters are also likely to vary, even for a single receptor-transducer system. It follows that the use of such agonist parameters to classify receptors or receptor-transducer systems is not likely to be totally dependable. It is also possible that agonists which interact with only one receptor-transducer system may show selectivity between tissues with different concentrations of G-proteins and of guanine nucleotides.     

Quantitative analysis of ultrasensitive responses

Ultrasensitive responses are common in cellular information transfer because they allow cells to decode extracellular stimuli in an all-or-none manner. Biochemical responses are usually analyzed by fitting the Hill equation, and the estimated Hill coefficient is taken as a measure of sensitivity. However, this approach is not appropriate if the response under consideration significantly deviates from the best-fit Hill equation. In addition, Hill coefficients greater than unity do not necessarily imply ultrasensitive behaviour if basal activation is significant. In order to circumvent these problems we propose a general method for the quantitative analysis of sensitivity, the relative amplification plot, which is based on the response coefficient defined in metabolic control analysis. To quantify sensitivity globally (i.e. over the whole stimulus range) we introduce the integral-based relative amplification coefficient. Our relative amplification approach can easily be extended to monotonically decreasing, bell-shaped or nonsaturated responses.     

Partial agonism: mechanisms based on ligand-receptor interactions and on stimulus-response coupling

Substances eliciting, at very high concentrations, a lower maximal response of a particular biological system than a defined standard, are defined as partial agonists. The convention rests on the definition of a standard substance that achieves a 'full' maximal response; partial agonism being, therefore, relative. Various mechanisms lie behind this phenomenon: 1. Receptor-related mechanisms: the agonist-receptor complex exists in several conformational states from which only one, or only a few, activate the cell signaling pathway. This may occur when the receptor itself, or the agonist, exists in multiple states (e.g., in the form of enantiomers or stereoisomers), or when the agonist-receptor complex changes its conformation (receptor switch: two-state model of receptor activation). Furthermore, a steric hindrance by a 'wrong-way binding' of a part of the agonist's molecules may prevent the full 'correct' occupancy of receptors. 2. Mechanisms based on the efficacy of the stimulus-response coupling. The efficacy is then proportional to the sum of probabilities that receptors in individual states activate the cell-signaling pathway. Doses (concentrations) eliciting the half maximal response (EC50), or similar response sensitivity parameters, are not included in the definition of partial agonism. However, tight correlations exist between maximal response and EC50 in many, but not all, generic groups of agonistically acting substances. These relationships are frequently linear; intercepts and slopes of these 'E, KE plots' are characteristic for individual, putative mechanisms. Dose-response curves of partial agonists are akin to those obtained for a response to a full agonist after a stepwise partial inactivation of receptors by an irreversible inhibitor. Also, the E, KE plots obtained in these instances are similar to those of partial agonists. The receptor reserve, rather vaguely defined in early reports, is therefore closely linked to the phenomenon of partial agonism.     

Synergistic responses in human platelets. Comparison between aggregation, secretion and cytosolic Ca2+ concentration

Synergistic interaction between ADP, adrenaline, 5-hydroxytryptamine (5HT) and [8-arginine]vasopressin is not observed for the aggregatory response of aspirin-treated human platelets when this response is estimated directly from the decrease in the number of single platelets in the suspension. This finding is in marked contrast with prior reports of synergistic interaction between these agonists when the rate and extent of the aggregometer response is estimated from the increase in the light transmittance of the suspension, using a platelet aggregometer. We propose that the apparent synergistic response detected using the aggregometer results from the inability of this instrument to respond during the initial phase of aggregation. Significant synergistic interaction is observed for the increase in cytosolic [Ca2+] induced by addition of the ADP/5HT and, to a lesser extent, of the ADP/vasopressin agonist pairs as compared with that caused by addition of the individual agonists. This effect is not, however, typical of the system since increases in cytosolic [Ca2+] induced by addition of the ADP/thrombin or 5HT/vasopressin agonist pairs are no greater than the sum of the responses to these agonists added separately. Addition of collagen prior to ADP or 11,9-epoxymethanoprostaglandin H2 (U46619) fails to enhance the increase in cytosolic [Ca2+] induced by these latter agonists. Adrenaline, when added prior to non-saturating concentrations of U46619, thrombin, vasopressin or ADP, significantly enhances the increase in cytosolic [Ca2+] induced by these agonists in platelets suspended in media containing less than 0.1 microM or 1 mM Ca2+. However, adrenaline fails to enhance the increase in cytosolic [Ca2+] induced by the divalent cation ionophore, ionomycin. Enhancement by adrenaline of Ca2+ influx induced by U46619, thrombin and ADP has been shown by using Mn2+ as probe. Adrenaline also enhances the extent of [3H]5HT secretion induced by U46619, thrombin and vasopressin but fails to increase that induced by ADP in this aspirin-treated preparation. These results are in part consistent with the postulate that adrenaline, acting via an alpha 2-adrenoceptor, modulates receptor--phospholipase-C coupling. However, such modulation does not appear to involve inhibition of adenylate cyclase.(ABSTRACT TRUNCATED AT 400 WORDS)     

Similar acute molecular responses to equivalent volumes of isometric, lengthening, or shortening mode resistance exercise.

The present study was undertaken to test the hypothesis that the contraction mode of action [static-isometric (Iso), shortening-concentric (Con), or lengthening-eccentric (Ecc)] used to stress the muscle provides a differential mechanical stimulus eliciting greater or lesser degrees of anabolic response at the initiation of a resistance training program. We performed an acute resistance training study in which different groups of rodents completed four training sessions in either the Iso, Con, or Ecc mode of contraction under conditions of activation and movement specifically designed to elicit equivalent volumes of force accumulation. The results of this experiment indicate that the three modes of contraction produced nearly identical cell signaling, indicative of an anabolic response involving factors such as increased levels of mRNA for IGF-I, procollagen III alpha1, decreased myostatin mRNA, and increased total RNA concentration. The resulting profiles collectively provide evidence that pure mode of muscle action, in and of itself, does not appear to be a primary variable in determining the efficacy of increased loading paradigms with regard to the initiation of selected muscle anabolic responses.     

Systems analysis of digoxin kinetics and inotropic response in the rat heart: effects of calcium and KB-R7943

To gain more insight into the mechanistic processes controlling the kinetics of inotropic response of digoxin in the perfused whole heart, an integrated kinetic model was developed incorporating digoxin uptake, receptor binding (Na(+)-K(+)-ATPase inhibition), and cellular events linking receptor occupation and response. The model was applied to data obtained in the single-pass Langendorff-perfused rat heart for external [Ca(2+)] of 0.5 and 1.5 mM under control conditions and in the presence of the reverse-mode Na(+)/Ca(2+) exchange inhibitor KB-R7943 (0.1 microM) in perfusate. Outflow concentration and left ventricular developed pressure data measured for three consecutive doses (15, 30, and 45 microg) in each heart were analyzed simultaneously. While disposition kinetics of digoxin was determined by interaction with a heterogeneous receptor population consisting of a high-affinity/low-capacity and a low-affinity/high- capacity binding site, response generation was >80% mediated by binding to the high-affinity receptor. Digoxin sensitivity increased at lower external [Ca(2+)] due to higher stimulus amplification. Coadministration of KB-R7943 significantly reduced the positive inotropic effect of digoxin at higher doses (30 and 45 microg) and led to a saturated and delayed receptor occupancy-response relationship in the cellular effectuation model. The results provide further evidence for the functional heterogeneity of the Na(+)-K(+)-ATPase and suggest that in the presence of KB-R7943 a reduction of the Ca(2+) influx rate via the reverse mode Na(+)/Ca(2+) exchanger might become the limiting factor in digoxin response generation.     

Decoding cortical neuronal signals: Network models,information estimation and spatial tuning

We have studied the encoding of spatial pattern information by complex cells in the primary visual cortex of awake monkeys. Three models for the conditional probabilities of different stimuli, given the neuronal response, were fit and compared using cross-validation. For our data, a feed-forward neural network proved to be the best of these models.The information carried by a cell about a stimulus set can be calculated from the estimated conditional probabilities. We performed a spatial spectroscopy of the encoding, examining how the transmitted information varies with both the average coarseness of the stimulus set and the coarseness differences within it. We find that each neuron encodes information about many features at multiple scales. Our data do not appear to allow a characterization of these variations in terms of the detection of simple single features such as oriented bars.     

An analysis of the interaction of a platinum complex and radiation with CHO cells using the molecular theory of cell survival.

A mathematical analysis is developed from the concepts of the molecular theory of cell survival to explain the cytotoxic action of a platinum complex on CHO cells and its synergistic interaction with radiation. In the analysis, it is assumed that both inter- and intra-strand cross-links induced by the Pt complex in the DNA contribute to cell killing and that the synergistic effect arises from the interaction between an intra-strand cross-link and a radiation-induced single-strand break on complementary strands of the DNA. The analysis is shown to be compatible with experimental results. Further, a mechanistic model is developed in an attempt to explain in more detail the processes which appear to be involved in the expression of the cytological damage.     

Mechanisms of ligand binding and efficacy at the human D2(short) dopamine receptor

Mechanisms of ligand binding and receptor activation for the human D2(short) dopamine receptor have been probed using two homologous series of monohydroxylated and dihydroxylated agonists (phenylethylamines and 2-dipropylaminotetralins). In ligand binding studies, the majority of compounds exhibited competition curves versus [3H]spiperone that were best fitted using a two site binding model. The compounds had different abilities (potencies and maximal effects) to stimulate [35S]GTPgammaS binding and to inhibit forskolin-stimulated cAMP accumulation. From the data it can be concluded that: (i) the ability of an agonist to stabilize receptor/G protein coupling can be used to predict agonist efficacy for some groups of compounds (2-dipropylaminotetralins) but not for others (phenylethylamines); (ii) the receptor may be activated by unhydroxylated compounds; (iii) single hydroxyl groups or pairs of hydroxyl groups on the agonist may contribute to binding affinity, potency and efficacy; and (iv) for the 2-dipropylaminotetralin series two modes of agonist/receptor interaction have been identified associated with different relative efficacy.     

The influence of affinity, efficacy, and slope factor on the estimates obtained by the receptorial responsiveness method (RRM): a computer simulation study

The receptorial responsiveness method (RRM) was proposed to characterize changes in the concentration of degradable agonists in the microenvironment of their receptors. The characterization is done by providing concentrations of a stable agonist for the same receptor that is equieffective with the change in concentration to be characterized. RRM is based on the analysis of concentration-effect (E/c) curves reflecting the simultaneous action of the degradable and the stable agonist. In the present study, we investigated whether dissimilar affinity and (or) efficacy of the coacting agonists as well as the steepness of the E/c curves influence the reliability of RRM. E/c curves were simulated based on the operational model and then analyzed with RRM. We found that dissimilarity in affinity of the coacting agonists did not affect the accuracy of RRM estimates. In contrast, accuracy of the estimation depended on the magnitude of the concentration to be assessed, the operational slope factor, and the operational efficacy ratio of the coacting agonists. However, our results suggest that proper choice of a stable agonist for a degradable one can help to ensure reliable results, since information about the change in concentration of a degradable agonist is otherwise difficult to obtain.     

Contractile Torque as a Steering Mechanism for Orientation of Adherent Cells

It is well established that adherent cells change their orientation in response to non-uniform substrate stretching. Most observations indicate that cells orient away from the direction of the maximal substrate strain, whereas in some cases cells also align with the direction of the maximal strain. Previous studies suggest that orientation and steering of the cell may be closely tied to cytoskeletal contractile stress but they could not explain the mechanisms that direct cell reorientation. This led us to develop a simple, mechanistic theoretical model that could predict a direction of cell orientation in response to mechanical nonuniformities of the substrate. The model leads to a simple physical mechanism -- namely the contractile torque -- that directs the cell toward a new orientation in response to anisotropic substrate stretching or substrate material anisotropy. A direction of the torque is determined by a dependence of the contractile stress on substrate strain. Model predictions are tested in the case of simple elongation of the substrate and found to be consistent with experimental data from the literature.     

Guanine nucleotide regulation of agonist binding to muscarinic cholinergic receptors. Relation to efficacy of agonists for stimulation of phosphoinositide breakdown and Ca2+ mobilization.

The efficacies of a series of six muscarinic cholinergic receptor agonists for stimulation of phosphoinositide breakdown and unidirectional efflux of 45Ca2+ in 1321N1 human astrocytoma cells were compared with the relative capacity of these agonists for formation of a GTP-sensitive high-affinity binding state in washed membranes. Carbachol and methacholine were 'full' agonists as regards phosphoinositide breakdown and Ca2+ mobilization, whereas bethanechol, arecoline and oxotremorine were 'partial' agonists for these two responses. Pilocarpine was the least efficacious of the six drugs tested. Except for pilocarpine, competition curves generated with the agonists and [3H]quinuclidinyl benzilate did not follow the Law of Mass Action for ligand interaction at a single site. Non-linear regression analyses of these data indicated that the data significantly better fit a two-, rather than a single-, site model with a high- and a low-affinity binding component. Competition curves generated in the presence of GTP were shifted to the right, and the extent of receptors in the high-affinity agonist-binding state was decreased. The relative efficacies of the six agonists for stimulation of phosphoinositide breakdown and Ca2+ mobilization were significantly correlated with the difference in affinities (KL/KH) between the two affinity states for each agonist. The relative efficacy of the agonists for stimulation of Ca2+ mobilization also was significantly correlated with the extent of receptors in the high-affinity state (%H) for each agonist. The results suggest that interaction with an as-yet unidentified guanine nucleotide regulatory protein is important in the mechanism whereby muscarinic receptors stimulate phosphoinositide breakdown in 1321N1 astrocytoma cells.     

Data-driven analysis in drug discovery

In the process of drug discovery for new chemical entities, application of appropriate pharmacological models often is not possible because the molecular mechanism of the compound is not yet elucidated. Therefore, a data-driven approach using generic tools designed to quantify characteristic patterns of concentration-response curves is required. This article outlines the options available for quantifying agonist and antagonist activity. Specifically, for agonists, the use of the Operational model for the determination of functional effects (equimolar potency ratios for full agonists, calculation of relative efficacy) is described. For antagonists, the measurement of pKB (-log of the equilibrium dissociation constant of the antagonist-receptor complex) for orthosteric antagonists that do not alter basal response (simple competitive antagonists), increase basal response (partial agonists), and decrease basal response (in constitutively active systems; inverse agonists) is discussed. In addition, this article considers methods to discern orthosteric receptor antagonism from allosteric antagonism whereby the agonist and antagonist bind to separate sites and interact through a conformational change in the receptor. Methods for the measurement of the pKB for allosteric modulators as well as co-operativity constants for these modulators is described.