ThermoBRET: A Ligand-Engagement Nanoscale Thermostability Assay Applied to GPCRs.

Measurements of membrane protein thermostability reflect ligand binding. Current thermostability assays often require protein purification or rely on pre-existing radiolabelled or fluorescent ligands, limiting their application to established targets. Alternative methods, such as fluorescence-detection size exclusion chromatography thermal shift, detect protein aggregation but are not amenable to high-throughput screening. Here, we present a ThermoBRET method to quantify the relative thermostability of G protein coupled receptors (GPCRs), using cannabinoid receptors (CB1 and CB2 ) and the ß2 -adrenoceptor (ß2 AR) as model systems. ThermoBRET reports receptor unfolding, does not need labelled ligands and can be used with non-purified proteins. It uses Bioluminescence Resonance Energy Transfer (BRET) between Nanoluciferase (Nluc) and a thiol-reactive fluorescent dye that binds cysteines exposed by unfolding. We demonstrate that the melting point (Tm ) of Nluc-fused GPCRs can be determined in non-purified detergent solubilised membrane preparations or solubilised whole cells, revealing differences in thermostability for different solubilising conditions and in the presence of stabilising ligands. We extended the range of the assay by developing the thermostable tsNLuc by incorporating mutations from the fragments of split-Nluc (Tm of 87 °C versus 59 °C). ThermoBRET allows the determination of GPCR thermostability, which is useful for protein purification optimisation and drug discovery screening.

Structural Basis of the Negative Allosteric Modulation of 5-BDBD at Human P2X4 Receptors.

The P2X4 receptor is a ligand-gated ion channel activated by extracellular ATP. P2X4 activity is associated with neuropathic pain, vasodilation, and pulmonary secretion and is therefore of therapeutic interest. The structure-activity relationship of P2X4 antagonists is poorly understood. Here we elucidate the structure-activity of 5-(3-bromophenyl)-1,3-dihydro-2H-benzofuro[3,2-e]-1,4-diazepin-2-one (5-BDBD) at human P2X4 by combining pharmacology, electrophysiology, molecular modeling, and medicinal chemistry. 5-BDBD antagonized P2X4 in a noncompetitive manner but lacked effect at human P2X2. Molecular modeling and site-directed mutagenesis suggested an allosteric binding site for 5-BDBD located between two subunits in the body region of P2X4, with M109, F178, Y300, and I312 on one subunit and R301 on the neighboring subunit as key residues involved in antagonist binding. The bromine group of 5-BDBD was redundant for the antagonist activity of 5-BDBD, although an interaction between the carbonyl group of 5-BDBD and R301 in P2X4 was associated with 5-BDBD activity. 5-BDBD could inhibit the closed channel but poorly inhibited the channel in the open/desensitizing state. We hypothesize that this is due to constriction of the allosteric site after transition from closed to open channel state. We propose that M109, F178, Y300, R301, and I312 are key residues for 5-BDBD binding; provide a structural explanation of how they contribute to 5-BDBD antagonism; and highlight that the limited action of 5-BDBD on open versus closed channels is due to a conformational change in the allosteric site. SIGNIFICANCE STATEMENT: Activity of P2X4 receptor is associated with neuropathic pain, inflammation, and vasodilatation. Molecular information regarding small-molecule interaction with P2X4 is very limited. Here, this study provides a structural explanation for the action of the small-molecule antagonist 5-BDBD at the human P2X4 receptor.

Mechano-sensitivity of ß2-adrenoceptors enhances constitutive activation of cAMP generation that is inhibited by inverse agonists.

The concept of agonist-independent signalling that can be attenuated by inverse agonists is a fundamental element of the cubic ternary complex model of G protein-coupled receptor (GPCR) activation. This model shows how a GPCR can exist in two conformational states in the absence of ligands; an inactive R state and an active R* state that differ in their affinities for agonists, inverse agonists, and G-protein alpha subunits. The proportion of R* receptors that exist in the absence of agonists determines the level of constitutive receptor activity. In this study we demonstrate that mechanical stimulation can induce ß2-adrenoceptor agonist-independent Gs-mediated cAMP signalling that is sensitive to inhibition by inverse agonists such as ICI-118551 and propranolol. The size of the mechano-sensitive response is dependent on the cell surface receptor expression level in HEK293G cells, is still observed in a ligand-binding deficient D113A mutant ß2-adrenoceptor and can be attenuated by site-directed mutagenesis of the extracellular N-glycosylation sites on the N-terminus and second extracellular loop of the ß2-adrenoceptor. Similar mechano-sensitive agonist-independent responses are observed in HEK293G cells overexpressing the A2A-adenosine receptor. These data provide new insights into how agonist-independent constitutive receptor activity can be enhanced by mechanical stimulation and regulated by inverse agonists.

Kinetic analysis of endogenous ß2 -adrenoceptor-mediated cAMP GloSensor™ responses in HEK293 cells.

BACKGROUND AND AIM: Standard pharmacological analysis of agonist activity utilises measurements of receptor-mediated responses at a set time-point, or at the peak response level, to characterise ligands. However, the occurrence of non-equilibrium conditions may dramatically impact the properties of the response being measured. Here we have analysed the initial kinetic phases of cAMP responses to ß2 -adrenoceptor agonists in HEK293 cells expressing the endogenous ß2 -adrenoceptor at extremely low levels. EXPERIMENTAL APPROACH: The kinetics of ß2 -adrenoceptor agonist-stimulated cAMP responses were monitored in real-time, in the presence and absence of antagonists, in HEK293 cells expressing the cAMP GloSensor™ biosensor. Potency (EC50 ) and efficacy (Emax ) values were determined at the peak of the agonist GloSensor™ response and compared to kinetic parameters L50 and IRmax values derived from initial response rates. KEY RESULTS: The partial agonists salbutamol and salmeterol displayed reduced relative IRmax values (with respect to isoprenaline) when compared with their Emax values. Except for the fast dissociating bisoprolol, preincubation with ß2 -adrenoceptor antagonists produced a large reduction in the isoprenaline peak response due to a state of hemi-equilibrium in this low receptor reserve system. This effect was exacerbated when IRmax parameters were measured. Furthermore, bisoprolol produced a large reduction in isoprenaline IRmax consistent with its short residence time. CONCLUSIONS AND IMPLICATIONS: Kinetic analysis of real-time signalling data can provide valuable insights into the hemi-equilibria that can occur in low receptor reserve systems with agonist-antagonist interactions, due to incomplete dissociation of antagonist whilst the peak agonist response is developing.

Combined docking and machine learning identify key molecular determinants of ligand pharmacological activity on ß2 adrenoceptor.

G protein-coupled receptors (GPCRs) are valuable therapeutic targets for many diseases. A central question of GPCR drug discovery is to understand what determines the agonism or antagonism of ligands that bind them. Ligands exert their action via the interactions in the ligand binding pocket. We hypothesized that there is a common set of receptor interactions made by ligands of diverse structures that mediate their action and that among a large dataset of different ligands, the functionally important interactions will be over-represented. We computationally docked ~2700 known ß2AR ligands to multiple ß2AR structures, generating ca 75 000 docking poses and predicted all atomic interactions between the receptor and the ligand. We used machine learning (ML) techniques to identify specific interactions that correlate with the agonist or antagonist activity of these ligands. We demonstrate with the application of ML methods that it is possible to identify the key interactions associated with agonism or antagonism of ligands. The most representative interactions for agonist ligands involve K972.68×67 , F194ECL2 , S2035.42×43 , S2045.43×44 , S2075.46×641 , H2966.58×58 , and K3057.32×31 . Meanwhile, the antagonist ligands made interactions with W2866.48×48 and Y3167.43×42 , both residues considered to be important in GPCR activation. The interpretation of ML analysis in human understandable form allowed us to construct an exquisitely detailed structure-activity relationship that identifies small changes to the ligands that invert their pharmacological activity and thus helps to guide the drug discovery process. This approach can be readily applied to any drug target.