Agonists in the Extended Conformation Stabilize the Active State of β-Adrenoceptors

In this study, we conducted a comparative analysis of the structure of agonists and antagonists of transmembrane (TM) β-adrenoceptors (β-ARs) and their interactions with the β-ARs and proposed the mechanism of receptor activation. A characteristic feature of agonist and antagonist molecules is the presence of a hydrophobic head (most often, one or two aromatic rings) and a tail with a positively charged amino group. All β-adrenergic agonists have two carbon atoms between the aromatic ring of the head and the nitrogen atom of the amino group. In antagonist molecules, this fragment can be either reduced or increased to four atoms due to the additional carbon and oxygen atoms. The agonist head, as a rule, has two H-bond donors or acceptors in the para- and meta-positions of the aromatic rings, while in the antagonist heads, these donors/acceptors are absent or located in other positions. Analysis of known three-dimensional structures of β-AR complexes with agonists showed that the agonist head forms two H-bonds with the TM5 helix, and the tail forms an ionic bond with the D3.32 residue of the TM3 helix and one or two H-bonds with the TM7 helix. The tail of the antagonist can form similar bonds, but the interaction between the head and the TM5 helix is much weaker. As a result of these interactions, the agonist molecule acquires an extended “strained string” conformation, in contrast to the antagonist molecule, which has a longer, bended, and flexible tail. The “strained string” of the agonist interacts with the TM6 helix (primarily with the W6.48 residue) and turns it, which leads to the opening of the G protein-binding site on the intracellular side of the receptor, while flexible and larger antagonist molecules do not have the same effect on the receptor.

INTRODUCTION G protein-coupled receptors (GPCRs) constitute the largest family of cell membrane receptors that includes over 800 human proteins targeted by at least 30% of current medicines (see reviews [1,2]). This might be the main reason why GPCRs have been extensively studied in several recent decades, resulting in the elucidation of many aspects of their biology, biochemistry, and pharma-cology. Originally, the behavior of GPCRs was described in terms of a simple two-state model; however, there is growing body of evidence that GPCRs are not molecular switches, but rather molecular relays, i.e., dynamic proteins with multiple states between active and inactive conformations [4][5][6][7][8].
Recent crystallography data have provided the snapshots of both active and inactive functional states of GP-CRs [9,10]. GPCR structures solved so far share the same overall fold -a bundle of seven transmembrane (TM) α-helices with three extra-and three intracellular loops. The extracellular interface is responsible for the ligand binding, while the intracellular portion interacts with G proteins, β-arrestins, and other downstream effectors. Analysis of known GPCR structures indicates that receptor activation is associated with subtle changes in the extracellular portion of the protein and extensive rearrangement of the TM helices on the cytoplasmic side [11,12]. Agonist binding at the GPCR extracellular interface results in the opening of the intracellular part for the G protein binding, which promotes G protein activation (GDP release) and initiates the signaling cascade.
The progress in membrane protein crystallography and related techniques in the past decade [7,13] has allowed to elucidate many aspects of GPCR structure, activation, and physiology; however, some details of ligand recognition and receptor activation remain poorly understood. One of the main aims of this study was comparative analysis of the structures of β-adrenoceptor (β-AR) ligands. Based on the results of this analysis, we proposed that agonist molecules adopt an extended ("strained string") conformation and stabilize the active state of β-AR, in contrast to antagonist molecules, which have longer and fl exible tail and fail to produce the same eff ect on the receptors.

MATERIALS AND METHODS
As the main research approaches in this study, we used stereochemical analysis of the known three-dimensional structures of β-AR complexes with the corresponding ligands and comparative analysis of the chemical structure of agonists and antagonists and their conformations in the complexes. For this, we created a database of such complexes that included 64 structures determined by crystallography so far and a database of β-AR sequences from the Swiss-Prot UniProt [14]. The atomic coordinates of the complexes were taken from the Protein Data Bank (PDB, URL: https://www.rcsb.org) [15]. The three-dimensional structures of the receptors and their ligands were analyzed using RasMol [16] and PyMOL Molecular Graphics System Version 1.4.1 Schrödinger, LLC. β-AR subtypes designated using the nomenclature recommended by the NC-IUPHAR Subcommittee on Adrenergic Receptors. Amino acid residues in β-ARs were designated according to the Ballesteros and Weinstein nomenclature [17]. Multiple sequence alignment was performed with the Clustal Omega program (1.2.4) included in the UniProt resource [14]. Images of ligand molecules were taken from the Drug Informational Portal, ChEBI, and ChemSpider.

RESULTS AND DISCUSSION
Comparative analysis of chemical structures of β-AR ligands. In terms of receptor activation (intrinsic effi ca-cy), GPCR ligands can be categorized into four groups: (i and ii) full and partial agonists that produce the maximal or sub-maximal functional response, respectively; (iii) inverse agonists that decrease the basal receptor activity (activity in the absence of ligand); and (iv) neutral antagonists, that compete with other ligands for the orthosteric binding site, but their interaction with the receptor does not result in the G protein binding. Table 1 shows chemical structures of β-AR agonists that have been co-crystallized with the receptor. A characteristic feature of these molecules is the presence of an aromatic head and a tail with a positively charged amine. The tail consists of ethanolamine core and various substitutes connected to the amine group. Hydroxyl groups in the paraand meta-positions of the catechol moieties and in the para-position in aromatic rings of non-catechol agonists can form hydrogen bonds with the receptor helices.
Aromatic heterocyclic heads of non-catechol agonists can have other donors and acceptors of hydrogen bonds that might be involved in the binding with the receptor. Hence, GPCR agonists have two centers of polar interactions with the receptor (donors/acceptors of H-bonds of the head and donors/acceptors and positively charged amine of the tail). One may ask if the distance between these centers has any infl uence on the specifi city of receptor-ligand interactions. As seen in Table 1, the length of the tail fragment between the N-atom and the aromatic ring is the same in all agonists and comprises three covalent bonds. In other words, the N-atom is separated by 2 carbons from the substituted benzene or other aromatic ring. Moreover, in all agonists with six-membered aromatic rings, the O-atoms of hydroxyls located in the paraand meta-positions and the N-atoms of the tail are separated by seven and six covalent bonds, respectively.
tFor comparison, Table 2 shows partial agonists, antagonists, and inverse agonists (referred to as β-blockers after Emtage et al. [28]. These compounds also have aromatic heads and tails with positively charged amine groups. However, compared to agonists, these ligands have either longer (Table 2) or shorter (e.g., doxepin and bretylium tosylate not shown here) tail fragments between the amine N-atom and the aromatic ring. In most cases, the N-atom of the amine group is separated from the aromatic ring by 4 atoms (often one of them is O). Moreover, antagonists and inverse agonists presented in Table 2 do not contain donors or acceptors of H-bonds in the para-position of the aromatic head or even lack them at all. Note that the NH-groups of cyanopindolol and carazolol are located in the meta-position relative to the tail. Another feature of β-blockers is that they typically have larger heads consisting of two and three rings as compared to the agonists.
Two centers of polar interactions between β-ARs and their agonists. The orthosteric binding pocket of β-ARs is located within the TM region and is primarily composed of the extracellular fragments of the TM helices 3, 5, 6,      and 7 (Fig. 1). Multiple biochemical and mutagenesis studies, as well as analysis of crystal structures of aminergic GPCRs, have allowed to locate critical amino acid residues in the TM helices [45]. It was demonstrated that the charged amine of the ligand interacts with Asp residue D3.32. GPCRs with D3.32 also have a Tyr residue at position Y7.43 and contain Asn at position N7.39, which suggests that their side chains are involved in the interactions with the positively charged moieties in the ligands. There-fore, these key amino acids form the center for the binding of amino groups and other polar groups of the ligand tail.
In the second center of polar interactions, all catecholamine receptors have Ser residues at positions S5.42 and S5.46; most of them also have Ser at position S5.43, so that these residues can form H-bonds with the donors and/or acceptors of the aromatic head. Recent data [10] indicate that Thr at position T3.37 also can interact with polar groups of the ligand head.  Figure 2 shows the alignment of amino acid sequences of TM helices 3, 5, 6, and 7 in β-AR subtypes β 1 , β 2 , and β 3 from diff erent animal species. All key amino acid residues mentioned above are highly conserved in all β-ARs. Figure 3 shows crystal structures of the agonist isoprenaline (panel a) and neutral antagonist cyanopindolol (panel b) bound in the main binding pocket of β 1 -Ars [18,33]. The fi gure shows the side chains of the amino acid residues forming the two centers of polar interactions as well as ligand conformations. As seen in Fig. 3b, the O-C-C-C group of atoms in cyanopindolol tail acquires a gauche-conformation (i.e., forms a kink in the tail), while the C-C-C-N group of atoms in the isoprenaline tail has a trans-conformation (Fig. 3a). Figure 4 shows the superimposed structures from Fig. 3. It is clearly seen that the O-C-C-C group of atoms forms a kink in the cyanopindolol tail, and the isoprenaline tail has an extended conformation. Analysis shows that the same picture is observed in other complexes of agonists and antagonists with the receptors (Tables 1 and 2).

Comparison of conformations of agonists, antagonists, and inverse agonists bound to β-ARs.
The amine groups of both ligands are located in the amine-binding center (D3.32, N7.39, Y7.43), while catechol hydroxyls of isoprenaline and NH-and cyanogroups of cyanopindolol are situated in the other center of polar interactions (S5.42, S5.43, S5.46). However, the overall geometries of the agonist and antagonist molecules are diff erent. In isoprenaline, the tail has an extended (trans-) conformation, and the molecule resembles a "strained string" that stabilizes the active state of the receptor. The tail of cyanopindolol is bent between the amine N-atom and the aromatic ring, resulting in the zigzag-like conformation (Figs. 3 and 4). Under other conditions (for example, in a native membrane, and not in a crystal), the tail of the antagonist can be transformed into an extended conformation and back, i.e., TM helices in the receptor complexes with antagonists can have a greater dynamic mobility than in the complexes with agonists (see Nygaard et al. (2013) [6]).

Fig. 2.
Multiple sequence alignment of TM3, TM5, TM6, and TM7 domains of β-ARs for diff erent animal species. Domain information and protein sequences were taken from Swiss-Prot UniProt Knowledgebase (35 proteins) [14]. Key amino acids in the sequences are shown with a gray background.  We have examined other agonists bound to the corresponding aminergic GPCRs and, indeed, the tail fragment between the amine N-atom and the aromatic ring in these complexes also has an extended conformation (PDB ID: 2Y00, 2Y02, 2Y04, 3P0G, 3PDS, 3SN6, and other structures presented in Table 1). On the other hand, antagonist molecules are bent in the tail fragments between the amine N-atom and the aromatic ring (PDB ID: 2VT4, 2YCW, 2RH1, 3D4S, 3NY8, 3NY9, 3NYA, 3PBL, 3RZE, and other structures presented in Table 2). These data suggest that the extended conformation of the agonist tail is of particular importance in the active state stabilization in β-ARs and other aminergic GPCRs. We can also speculate that the distance between the two centers of polar interactions with the agonists should be equal (or close) to the distance between the corresponding centers in the activated GPCR.
The role of the agonist "strained string" conformation in the β-AR activation. Simple geometry considerations outlined above suggest that the agonist molecule can act as a "strained string" that stabilizes the arrangement of TM helices 3, 5, 6, and 7 corresponding to the receptor active state (Fig. 1). It appears that the interaction with the agonist bring the TM5 helix closer to the TM3 and TM7 helices. In particular, this "strained string" interacts with Trp at position W6. 48 and Phe residues at positions F6.51 and F6.52. These interactions are likely to be responsible for the rearrangement of TM6, i.e., for its rotation and/or vertical see-saw movement around a pivot in the middle of the membrane, which results in the opening of the intracellular portion for the G protein binding [46][47][48]. The interactions with W6.48 result in a subtle rotation of TM6 in the extracellular portion, which is amplifi ed towards the cytoplasmic side by the characteristic kink in the helix introduced by Pro residue P6.50 [11]. Unlike an agonist, an antagonist or an inverse agonist has a longer tail fragments with a loose conformation, a kink in the middle of the tail, and weak polar interactions of the head with the TM5 helix. Therefore, they cannot act as a "strained string", but instead occupy the active site of the receptor due to the polar and hydrophobic interactions.
Based on these fi ndings, we propose several predictions on how minor diff erences in the ligand structure can infl uence its functional characteristics: 1. Elongation of the agonist tail between the amine N-atom and the aromatic ring could result in the agonist conversion into antagonist.
2. Shortening of the tail fragment in the antagonist or inverse agonist with polar moieties in the corresponding positions of their heads may result in their transformation into agonists or partial agonists.
3. Removal of catechol hydroxyls or the corresponding polar substitutes in other aromatic heads could reduce the agonist activity. Modifi cations of these polar substitutes with aromatic or aliphatic groups substitutes, especially bulky ones, are likely to have a similar eff ect.
4. Modifi cation of amine groups of ligands (both agonists and antagonists) with bulky substitutes reduce or even prevent ligand binding to GPCRs.
5. The trans-conformation (extended conformation) of the agonist tail can be transformed into a gaucheor cis-conformation by chemical modifi cation, which can result in the reduction of agonist activity or even transform it into an antagonist.
Funding. This work was supported by the Russian Foundation for Basic Research (project no. 20-04-00453).
Ethics declarations. The authors declare no confl icts of interest. This article does not contain descriptions of studies involving humans or animals as study subjects.
Open access. This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution, and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third-party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.