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.     

The quantitative patterns of the modified one-center bioreceptor model

A modified one site model of the bioreceptor have been used to estimate quantitatively the phenomenon of the full agonism. Threshold phenomenon, spare receptors and linear dependence of the biological effect on concentration of complex agonist-receptor has been determinate by the general correlation equations. The equations of one site model provides a good fit to the experimental curves "dose-response" for the full agonists and allows to calculate the value of space receptors. The model includes occupancy Clark's theory and law "all or nothing". The interaction of acetylcholine and aklyltrimethylammonium salts with muscarinic acetylholine receptors is analysed as an example of use of this equation.     

Partial Agonism: Mechanisms Based on Ligand-Receptor Interactions and on Stimulus-Response Coupling

Abstract

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 (EC₅₀), or similar response sensitivity parameters, are not included in the definition of partial agonism. However, tight correlations exist between maximal response and EC₅₀ in many, but not all, generic groups of agonistically acting substances. These relationships are frequently linear; intercepts and slopes of these ‘E, K_E 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, K_E 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.     

Quantifying agonist activity at G protein-coupled receptors

When an agonist activates a population of G protein-coupled receptors (GPCRs), it elicits a signaling pathway that culminates in the response of the cell or tissue. This process can be analyzed at the level of a single receptor, a population of receptors, or a downstream response. Here we describe how to analyze the downstream response to obtain an estimate of the agonist affinity constant for the active state of single receptors. Receptors behave as quantal switches that alternate between active and inactive states (Figure 1). The active state interacts with specific G proteins or other signaling partners. In the absence of ligands, the inactive state predominates. The binding of agonist increases the probability that the receptor will switch into the active state because its affinity constant for the active state (K(b)) is much greater than that for the inactive state (K(a)). The summation of the random outputs of all of the receptors in the population yields a constant level of receptor activation in time. The reciprocal of the concentration of agonist eliciting half-maximal receptor activation is equivalent to the observed affinity constant (K(obs)), and the fraction of agonist-receptor complexes in the active state is defined as efficacy (ε) (Figure 2). Methods for analyzing the downstream responses of GPCRs have been developed that enable the estimation of the K(obs) and relative efficacy of an agonist. In this report, we show how to modify this analysis to estimate the agonist K(b) value relative to that of another agonist. For assays that exhibit constitutive activity, we show how to estimate K(b) in absolute units of M(-1). Our method of analyzing agonist concentration-response curves consists of global nonlinear regression using the operational model. We describe a procedure using the software application, Prism (GraphPad Software, Inc., San Diego, CA). The analysis yields an estimate of the product of K(obs) and a parameter proportional to efficacy (τ). The estimate of τK(obs) of one agonist, divided by that of another, is a relative measure of K(b) (RA(i)). For any receptor exhibiting constitutive activity, it is possible to estimate a parameter proportional to the efficacy of the free receptor complex (τ(sys)). In this case, the K(b) value of an agonist is equivalent to τK(obs)/τ(sys). Our method is useful for determining the selectivity of an agonist for receptor subtypes and for quantifying agonist-receptor signaling through different G proteins.     

Potential errors in agonist dissociation constant estimation caused by desensitization

A simple model of receptor desensitization is developed and analysed to predict the consequences of acute receptor loss on the pharmacological quantification of agonist action. The model incorporates the desensitization scheme of Katz and co-workers (1957) into the operational model of agonism (Black & Leff, 1983) and, therefore, it assumes the occupancy theory of agonist action. Pharmacological effect-time profiles are simulated which illustrate (i) the extent to which overt fade may be detectable under different conditions and (ii) the extent to which measured pharmacological effects deviate from those which would be measured in the absence of desensitization. It is shown that the resulting agonist concentration-effect curves may be displaced rightwards from their equilibrium positions and that the agonist dissociation constants estimated from them may be overestimated. Such errors are predicted to occur regardless of whether or not fade is detectable in the effect-time profiles. Estimates of agonist efficacy appears to be unaffected by desensitization. The results of this analysis are discussed with respect to their implications for receptor classification using agonist dissociation constant estimates and for the development of agonist drugs.     

Single ion channel's view of classical receptor theory.

The nature of drug agonism has been the central mystery of two conceptually different approaches: classical receptor theory, which does not require any knowledge of mechanism, and mechanistic theories, which do. Ligand-activated ion channel macromolecules that contain both the agonist receptor site and molecular machinery to generate a response present a unique experimental system with which to explore agonism and antagonism. Electrical recordings from one channel at a time offer phenomena and a perspective quite different from that usually encountered in studies of the drug-receptor interaction. This review describes patch-clamp recordings that illustrate the ligand-evoked behavior that gives rise to classical phenomena. A comparison of the channel currents recorded in the presence of different agonists reveals how these drugs act as full or partial agonists. At higher concentrations of agonist, the conformational and kinetic transitions that underlie desensitization can be observed. Receptor conformational changes induced by agonist and antagonist binding further expand our ideas about what these drugs do, and contribute to the growth of concepts that will further our understanding of drug agonism.     

CRH-induced electrical activity and calcium signalling in pituitary corticotrophs

Pituitary corticotroph cells generate repetitive action potentials and associated Ca2+ transients in response to the agonist corticotropin releasing hormone (CRH). There is indirect evidence suggesting that the agonist, by way of complex intracellular mechanisms, modulates the voltage sensitivity of the L-type Ca2+ channels embedded in the plasma membrane. We have previously constructed a Hodgkin-Huxley-type model of this process, which indicated that an increase in the L-type Ca2+ current is sufficient to generate repetitive action potentials (LeBeau et al. (1997). Biophys. J.73, 1263-1275). CRH is also believed to inhibit an inwardly rectifying K+ current. In this paper, we have found that a CRH-induced inhibition of the inwardly rectifying K+ current increases the model action potential firing frequency, [Ca2+]i transients and membrane excitability. This dual modulatory action of CRH on inward rectifier and voltage-gated Ca2+ channels better describes the observed CRH-induced effects. This structural alteration to the model along with parameter changes bring the model firing frequency in line with experimental data. We also show that the model exhibits experimentally observed bursting behaviour, where the depolarization spike is followed by small oscillations in the membrane potential.     

Allosteric agonist modulators

This article discusses a model to describe the effects of molecules that bind to a site on the receptor separate from that of the endogenous agonist to actively produce receptor signals (direct agonism). In addition, these molecules also modify the biological responses of the endogenous agonist (either potentiation or antagonism). The effects of such compounds in high-throughput screening assays are described as well as their effects on the dose-response curves to conventional agonists.     

Allosteric Agonist Modulators

This article discusses a model to describe the effects of molecules that bind to a site on the receptor separate from that of the endogenous agonist to actively produce receptor signals (direct agonism). In addition, these molecules also modify the biological responses of the endogenous agonist (either potentiation or antagonism). The effects of such compounds in high-throughput screening assays are described as well as their effects on the dose-response curves to conventional agonists.     

Receptor-agonist interactions in service-theoretic perspective, effects of molecular timing on the shape of dose-response curves

Service-theoretic concepts and methods, widely used in other fields (e.g., telecommunication and operations research), are useful also in a biochemical setting because the treatment of biocatalysts (enzymes, receptors) as servers and their ligands as customers, based on the established formal methods of service or queuing theory, may lead to insights and results unobtainable by conventional, mass-action-law-based theories. In this article, we apply the service-theoretic approach to receptor-agonist systems and show how by changing the stochastic time pattern of "operationally relevant" point events (e.g., instants of agonist arrival, instants of post-climax agonist departure) a great variety of dose-response curves may be generated, even in very simple reaction schemes, which, according to mass action kinetics, invariably lead to hyperbolic r(A) curves (r and A stand for response and agonist concentration, respectively). The molecular timing inherent to a hyperbolic response system is not optimal: for instance, at the agonist concentration A(50), half of the agonist molecules are rejected ("lost") because of unfortunate timing of the arrival events. The fraction of lost arrivers can be diminished considerably if the arrivals are better timed: "sub-Poisson" arrivals improve the timing and, thus, convert hyperbolic r(A) curves into "lifted" nonhyperbolic ones. Conversely, "super-Poisson" arrivals make the non-optimal timing in hyperbolic response systems even worse and, thus, convert hyperbolic r(A) curves into "depressed" nonhyperbolic ones. Furthermore, under special timing conditions, nonhyperbolic r(A) curves can be generated, which are partly lifted, partly depressed relative to the reference hyperbola, and which resemble in shape well-known nonhyperbolic forms of enzyme and receptor kinetics (negatively cooperative, positively cooperative, and sigmoidal kinetics). In addition unusual (undulatory and sawtooth-like) r(A) curves can be generated solely by changing the temporal pattern of arrival and service completion instants. Virtually any shape of dose-response curves may be obtained by allowing for probability distributions whose characteristic shape varies with their mean; we call such distributions "variomorphic" and apply them to the arrival process of agonist molecules.     

Identification of an efficacy switch region in the ghrelin receptor responsible for interchange between agonism and inverse agonism

The carboxyamidated wFwLL peptide was used as a core ligand to probe the structural basis for agonism versus inverse agonism in the constitutively active ghrelin receptor. In the ligand, an efficacy switch could be built at the N terminus, as exemplified by AwFwLL, which functioned as a high potency agonist, whereas KwFwLL was an equally high potency inverse agonist. The wFw-containing peptides, agonists as well as inverse agonists, were affected by receptor mutations covering the whole main ligand-binding pocket with key interaction sites being an aromatic cluster in transmembrane (TM)-VI and -VII and residues on the opposing face of TM-III. Gain-of-function in respect of either increased agonist or inverse agonist potency or swap between high potency versions of these properties was obtained by substitutions at a number of positions covering a broad area of the binding pocket on TM-III, -IV, and -V. However, in particular, space-generating substitutions at position III:04 shifted the efficacy of the ligands from inverse agonism toward agonism, whereas similar substitutions at position III: 08, one helical turn below, shifted the efficacy from agonism toward inverse agonism. It is suggested that the relative position of the ligand in the binding pocket between this "efficacy shift region" on TM-III and the opposing aromatic cluster on TM-VI and TM-VII leads either to agonism, i.e. in a superficial binding mode, or it leads to inverse agonism, i.e. in a more profound binding mode. This relationship between different binding modes and opposite efficacy is in accordance with the Global Toggle Switch model for 7TM receptor activation.     

Effect of asymmetry of concentration-response curves on the results obtained by the receptorial responsiveness method (RRM): an in silico study

The receptorial responsiveness method (RRM) was proposed to estimate changes in the concentration of an agonist in the microenvironment of its receptor. Usually, this is done by providing the equieffective concentration of another agonist for the same receptor or for a largely overlapping postreceptorial signaling ("test agonist"). The RRM is a special nonlinear regression algorithm to analyze a concentration-response (E/c) curve that represents the simultaneous actions of a single agonist concentration to be estimated and of increasing concentrations of the test agonist. The aim of this study was to explore whether asymmetry of the E/c curve to be analyzed influences the reliability of the RRM. For this purpose, computer simulation was performed by constructing symmetric and asymmetric E/c curves using the operational model of agonism, and then these curves were analyzed with the RRM. To perform the RRM, 2 types of equations were used: one involving the Hill equation, the simplest model of the E/c relationship, and one containing the Richards equation, an advanced model properly handling E/c curve asymmetry. Results of this study indicate that E/c curve asymmetry does not significantly influence the accuracy of the estimates provided by the RRM. Thus, when using the RRM, it is not necessary to replace the Hill equation with the Richards equation to obtain useful estimates. Furthermore, it was found that estimation of a high concentration of a high-efficacy agonist can fail when the RRM is performed with a low-efficacy test agonist in a system characterized by a small operational slope factor.     

A mathematical analysis of factors determining the time-course of membrane conductance change evoked by electrophoretic drug application.

The paper discusses a mathematical model in which an attempt was made to simulate electrophoretic application of drugs to the chemosensitive area of muscle fibre and examines the resulting conductance change.The simulated conductance change of the chemosensitive membrane was studied as the function of receptor distribution, drug-receptor cooperative interaction, the duration of agonist release pulse and the position of the microelectrode tip. It is shown that the time-course of conductance change does not depend only on the kinetics of agonist-receptor interaction and co-operativity but also on the geometry of receptor distribution and the position of the microelectrode tip. These factors should be also considered in future electrophoretic experiments designed to study agonist-receptor interaction.     

Receptor-Agonist Interactions in Service-Theoretic Perspective,Effects of Molecular Timing on the Shape of Dose-Response Curves

Service-theoretic concepts and methods, widely used in other fields (e.g., telecommunication and operations research), are useful also in a biochemical setting because the treatment of biocatalysts (enzymes, receptors) as servers and their ligands as customers, based on the established formal methods of service or queuing theory, may lead to insights and results unobtainable by conventional, mass-action-law-based theories. In this article, we apply the service-theoretic approach to receptor-agonist systems and show how by changing the stochastic time pattern of “operationally relevant” point events (e.g., instants of agonist arrival, instants of postclimax agonist departure) a great variety of dose-response curves may be generated, even in very simple reaction schemes, which, according to mass action kinetics, invariably lead to hyperbolic r(A) curves (r and A stand for response and agonist concentration, respectively). The molecular timing inherent to a hyperbolic response system is not optimal: for instance, at the agonist concentration A₅₀, half of the agonist molecules are rejected (“lost”) because of unfortunate timing of the arrival events. The fraction of lost arrivers can be diminished considerably if the arrivals are better timed: “sub-Poisson” arrivals improve the timing and, thus, convert hyperbolic r(A) curves into “lifted” nonhyperbolic ones. Conversely, “super-Poisson” arrivals make the nonoptimal timing in hyperbolic response systems even worse and, thus, convert hyperbolic r(A) curves into “depressed” nonhyperbolic ones. Furthermore, under special timing conditions, nonhyperbolic r(A) curves can be generated, which are partly lifted, partly depressed relative to the reference hyperbola, and which resemble in shape well-known nonhyperbolic forms of enzyme and receptor kinetics (negatively cooperative, positively cooperative, and sigmoidal kinetics). In addition unusual (undulatory and sawtooth-like) r(A) curves can be generated solely by changing the temporal pattern of arrival and service completion instants. Virtually any shape of dose-response curves may be obtained by allowing for probability distributions whose characteristic shape varies with their mean; we call such distributions “variomorphic” and apply them to the arrival process of agonist molecules.     

Interpretation of agonist affinity estimations: the question of distributed receptor states

In the receptor-transducer model of pharmacological agonism, rejection of the traditional assumption that receptor molecules are in vast excess of transducer molecules permits the receptors to become distributed among unbound, bound and complexed states. Under these conditions, agonist affinities are liable to be overestimated when the method of irreversible receptor antagonism is used. Graphical tests have been developed to detect distribution, and these were applied to experimental data from the interaction between 5-HT and phenoxybenzamine on aortic tissue. Significant receptor distribution was not detected by the method. However, in the model it was assumed that there was a linear relation between the concentration of ternary complex and pharmacological effect. If this relation was replaced with a saturable one the effect of receptor distribution would be masked. The implications for pharmacologists and medicinal chemists are discussed.     

Synthesis and SAR of analogues of the M1 allosteric agonist TBPB. Part I: Exploration of alternative benzyl and privileged structure moieties

This Letter describes the first account of the synthesis and SAR, developed through an iterative analogue library approach, of analogues of the highly selective M1 allosteric agonist TBPB. With slight structural changes, mAChR selectivity was maintained, but the degree of partial M1 agonism varied considerably.     

On ballistic parameters of less lethal projectiles influencing the severity of thoracic blunt impacts

The development and safety certification of less lethal projectiles require an understanding of the influence of projectile parameters on projectile–chest interaction and on the resulting terminal effect. Several energy-based criteria have been developed for chest injury assessment. Many studies consider kinetic energy (KE) or energy density as the only projectile parameter influencing terminal effect. In a common KE range (100–160 J), analysis of the firing tests of two 40 mm projectiles of different masses on animal surrogates has been made in order to investigate the severity of the injuries in the thoracic region. Experimental results have shown that KE and calibre are not sufficient to discriminate between the two projectiles as regards their injury potential. Parameters, such as momentum, shape and impedance, influence the projectile–chest interaction and terminal effect. A simplified finite element model of projectile–structure interaction confirms the experimental tendencies. Within the range of ballistic parameters used, it has been demonstrated that maximum thoracic deflection is a useful parameter to predict the skeletal level of injury, and it largely depends on the projectile pre-impact momentum. However, numerical simulations show that these results are merely valid for the experimental conditions used and cannot be generalised. Nevertheless, the transmitted impulse seems to be a more general factor governing the thorax deflection.     

Membrane permeability and receptor function.

This paper presents a model of receptor function for receptors that mediate selective increases in ionic permeabilities (e.g. receptors that are sensitive to acetylcholine and gamma-aminobutyric acid). The model distinguishes itself from previous models by the qualitative and quantitative application of current biophysical approaches to ion permeation. Explicit consideration of receptor-controlled permeation sites and membrane permeability results in a specific and testable theory into which complex receptor phenomena can easily be incorporated.Examination of agonist concentration-conductance data from a well-studied GABAnergic receptor allows invalidation of a standard assumption for determining kinetic constants of agonist-receptor interactions. Specifically, calculations are presented which indicate that the size of an individual receptor (elementary) conductance event can vary with agonist concentration. The following section develops the expressions relating agonist concentration to total membrane conductance for possible receptor models which incorporate variations in the amplitude of elementary conductance events.In the case of desensitization, the “inactivated” receptor is suggested to result from occupation and inactivation, not at the agonist-binding (activation) site of the receptor, but at the ionophore. The ionophore occupation model provides a parsimonious explanation of different types of desensitization. In addition, incorporation of this model into the cyclic reaction model of Katz &; Thesleff (1957) provides increased specification for the cyclic model and results in new experimental predictions.     

Influence of adenine nucleotides on oligomycin inhibition of energy-transducing reactions in intact rat-liver mitochondria.

The inhibitory effect of oligomycin was investigated in intact mitochondria through oxidative phosphorylation and uncoupler induced ATPase activity. Results show that oligomycin inhibition curves can be either sigmoidal or hyperbolic depending on experimental conditions and chiefly on the metabolic state of mitochondria with regard to the distribution of mitochondrial endogenous adenine-nucleotides. Active respiration and uncoupler-induced ATPse activity produce sigmoidal titration curves for a high initial ATP : ADP ratio and hyperbolic curves for a low ATP : ADP ratio. Time-dependent inhibitions are observed for the two reactions. The maximal inhibitory action for low concentrations of the inhibitor is delayed by the initial presence of ATP or the possibility of generating from inorganic phosphate before adding oligomycin. Results presented here show that the initial adenine-nucleotide distribution is important for oligomycin sensitivity of energy-linked reactions. Although a limited conformational change of the oligomycin-sensitivity to the inhibitor, it is more likely that a gross structural change of the inner membrane induced by adenine-nucleotides modifies membrane permeability to oligomycin.