Have many estimates of efficacy and affinity been misled? Revisiting the operational model of agonism

A critical re-evaluation of the slope factor of the operational model of agonism: When to exponentiate operational efficacy

Agonist efficacy denoting the "strength" of agonist action is a cornerstone in the proper assessment of agonist selectivity and signalling bias. The simulation models are very accurate but complex and hard to fit experimental data. The parsimonious operational model of agonism (OMA) has become successful in the determination of agonist efficacies and ranking them. In 1983, Black and Leff introduced the slope factor to the OMA to make it more flexible and allow for fitting steep as well as flat concentration-response curves. First, we performed a functional analysis to indicate the potential pitfalls of the OMA. Namely, exponentiation of operational efficacy may break relationships among the OMA parameters. The fitting of the Black & Leff equation to the theoretical curves of several models of functional responses and the experimental data confirmed the fickleness of the exponentiation of operational efficacy affecting estimates of operational efficacy as well as other OMA parameters. In contrast, fitting The OMA based on the Hill equation to the same data led to better estimates of model parameters. In conclusion, Hill equation-based OMA should be preferred over the Black & Leff equation when functional-response curves differ in the slope factor. Otherwise, the Black & Leff equation should be used with extreme caution acknowledging potential pitfalls.

Insights into the operational model of agonism of receptor dimers

INTRODUCTION

Accurate ranking of efficacies and potencies of agonists is essential in the discovery of new selective agonists. For the purpose of system-independent ranking of agonists, the operational model of agonism (OMA) has become a standard. Many receptors function as oligomers which makes functional responses more complex, requiring an extension of the original OMA.

AREAS COVERED

Explicit equations of the operational model of agonism of receptor dimers (OMARD) were derived. The OMARD can be applied to any receptor possessing two orthosteric sites. The behavior of OMARD was analyzed to demonstrate its complexity and relation to experimental data. Properties of OMARD and OMA equations were compared to demonstrate their pros and cons.

EXPERT OPINION

Extension of OMA by slope factors gives simple equations of functional response that are easy to fit experimental data but results may be inaccurate because of exponentiation of operational efficacy. Also, such equations cannot accommodate bell-shaped curves. Explicit equations of OMARD give accurate results but are complex and tedious to fit experimental data. All operational models use inter-dependent parameters that are a hurdle in the fitting. A good understanding of OMARD behavior helps to overcome such obstacles.

Structural Assessment of Agonist Efficacy in the μ-Opioid Receptor: Morphine and Fentanyl Elicit Different Activation Patterns

Over the past two decades, the opioid epidemic in the United States and Canada has evidenced the need for a better understanding of the molecular mechanisms of medications used to fight pain. Morphine and fentanyl are widely used in opiate-mediated analgesia for the treatment of chronic pain. These compounds target the μ-opioid receptor (MOR), a class A G protein-coupled receptor (GPCR). In light of described higher efficacy of fentanyl with respect to morphine, we have performed independent μs-length unbiased molecular dynamics (MD) simulations of MOR complexes with each of these ligands, including the MOR antagonist naltrexone as a negative control. Consequently, MD simulations totaling 58 μs have been conducted to elucidate at the atomic level ligand-specific receptor activity and signal transmission in the MOR. In particular, we have identified stable binding poses of morphine and fentanyl, which interact differently with the MOR. Different ligand-receptor interaction landscapes directly induce sidechain conformational changes of orthosteric pocket residues: Asp149^3.32, Tyr150^3.33, Gln126^2.60, and Lys235^5.39. The induced conformations determine Asp149^3.32-Tyr328^7.43 sidechain-sidechain interactions and Trp295^6.48-Ala242^5.46 sidechain-backbone H-bond formations, as well as Met153^3.36 conformational changes. In addition to differences in ligand binding, different intracellular receptor conformational changes are observed as morphine preferentially activates transmembrane (TM) helices: TM3 and TM5, while fentanyl preferentially activates TM6 and TM7. As conformational changes in TM6 and TM7 are widely described as being the most crucial aspect in GPCR activation, this may contribute to the greater efficacy of fentanyl over morphine. These computationally observed functional differences between fentanyl and morphine may provide new avenues for the design of safer but not weaker opioid drugs because it is desirable to increase the safety of medicines without sacrificing their efficacy.

Quantifying the allosteric interactions within a G-protein-coupled receptor heterodimer

G-protein-coupled receptors are central to signal transduction and cell communication. The possibility that cells use receptor heteromerization to modulate individual receptor pathways is a surmise that cannot be precluded. Given the complexity of these processes, mathematical models contribute to understanding how receptors and their respective ligands regulate signaling. Here, a mathematical model is presented that quantifies the allosteric interactions within a receptor heterodimer. The model is based on the operational model of allosterism including constitutive receptor activity, which provides the pharmacological analysis of heteromerization with well-established and widely used modeling and fitting procedures.