Different conformational responses of the β2-adrenergic receptor-Gs complex upon binding of the partial agonist salbutamol or the full agonist isoprenaline

G Protein-Coupled Receptor-Ligand Pose and Functional Class Prediction

G protein-coupled receptor (GPCR) transmembrane protein family members play essential roles in physiology. Numerous pharmaceuticals target GPCRs, and many drug discovery programs utilize virtual screening (VS) against GPCR targets. Improvements in the accuracy of predicting new molecules that bind to and either activate or inhibit GPCR function would accelerate such drug discovery programs. This work addresses two significant research questions. First, do ligand interaction fingerprints provide a substantial advantage over automated methods of binding site selection for classical docking? Second, can the functional status of prospective screening candidates be predicted from ligand interaction fingerprints using a random forest classifier? Ligand interaction fingerprints were found to offer modest advantages in sampling accurate poses, but no substantial advantage in the final set of top-ranked poses after scoring, and, thus, were not used in the generation of the ligand-receptor complexes used to train and test the random forest classifier. A binary classifier which treated agonists, antagonists, and inverse agonists as active and all other ligands as inactive proved highly effective in ligand function prediction in an external test set of GPR31 and TAAR2 candidate ligands with a hit rate of 82.6% actual actives within the set of predicted actives.

State-dependent dynamics of extramembrane domains in the β2 -adrenergic receptor

G protein-coupled receptors (GPCRs) are membrane-bound signaling proteins that play an essential role in cellular signaling processes. Due to their intrinsic function of transmitting internal signals in response to external cues, these receptors are adapted to be highly dynamic in nature. The β2 -adrenergic receptor (β2 AR) is a representative member of the family that has been extensively analyzed in terms of its structure and activation. Although the structure of the transmembrane domain has been characterized in the different functional states of the receptor, the conformational dynamics of the extramembrane domains, especially the intrinsically disordered regions are still emerging. In this study, we analyze the state-dependent dynamics of extramembrane domains of β2 AR using atomistic molecular dynamics simulations. We introduce a parameter, the residue excess dynamics that allows us to better quantify receptor dynamics. Using this measure, we show that the dynamics of the extramembrane domains are sensitive to the receptor state. Interestingly, the ligand-bound intermediateR ' state shows the maximal dynamics compared to either the active R*G or inactive R states. Ligand binding appears to be correlated with high residue excess dynamics that are dampened upon G protein coupling. The intracellular loop-3 (ICL3) domain has a tendency to flip towards the membrane upon ligand binding, which could contribute to receptor "priming." We highlight an important ICL1-helix-8 interplay that is broken in the ligand-bound state but is retained in the active state. Overall, our study highlights the importance of characterizing the functional dynamics of the GPCR loop domains.

Predicting the duration of action of β2-adrenergic receptor agonists: Ligand and structure-based approaches

Agonists of the β2 adrenergic receptor (ADRB2) are an important class of medications used for the treatment of respiratory diseases. They can be classified as short acting (SABA) or long acting (LABA), with each class playing a different role in patient management. In this work we explored both ligand-based and structure-based high-throughput approaches to classify β2-agonists based on their duration of action. A completely in-silico prediction pipeline using an AlphaFold generated structure was used for structure-based modelling. Our analysis identified the ligands' 3D structure and lipophilicity as the most relevant features for the prediction of the duration of action. Interaction-based methods were also able to select ligands with the desired duration of action, incorporating the bias directly in the structure-based drug discovery pipeline without the need for further processing.

Activation/Deactivation Free-Energy Profiles for the β2-Adrenergic Receptor: Ligand Modes of Action

We use enhanced-sampling simulations with an effective collective variable to study the activation of the β2-adrenergic receptor in the presence of ligands with different efficacy. The free-energy profiles are computed for the ligand-free (apo) receptor and binary (apo-receptor + G-protein α-subunit and receptor + ligand) and ternary complexes. The results are not only compatible with available experiments but also allow unprecedented structural insight into the nature of GPCR conformations along the activation pathway and their role in the activation mechanism. In particular, the simulations reveal an unexpected mode of action of partial agonists such as salmeterol and salbutamol that arises already in the binary complex without the G-protein. Specific differences in the polar interactions with residues in TM5, which are required to stabilize an optimal TM6 conformation that facilitates G-protein binding and receptor activation, play a major role in differentiating them from full agonists.

Regulation of T Helper Cell Type 2 Immune Response by Controlling Beta-2 Adrenergic Receptor in Dendritic Cells of Patients with Allergic Rhinitis

INTRODUCTION

Allergic diseases are mediated by T helper cell type 2 (Th2) cells, which are differentiated by dendritic cells (DCs). Recently, it was reported that cAMP concentration in DCs is important for inducing allergic responses. However, the regulatory function of cAMP in DCs in Th2 immune responses is unclear. It was hypothesized that the regulation of G protein-coupled receptors (GPCRs) to increase cAMP levels in DCs would reduce Th2 immune responses.

METHODS

Human DCs from patients with allergic rhinitis (AR) and from healthy controls were subjected to next-generation sequencing (NGS) to identify potential GPCR. To investigate the functions of GPCR agonists, the in vitro co-culture experiment that THP-1 cells were differentiated into DCs and cultured with human CD4+ T-cells and an AR animal in vivo model were used.

RESULTS

Among the GPCRs, the beta-2 adrenergic receptor (ADRB2) of allergic DCs was significantly increased by NGS analysis. The expression of ADRB2 was also increased in Der p 1-treated DCs, which was reduced by treatment with the ADRB2 agonist salbutamol. Salbutamol treatment induced cAMP production in THP-1 derived DCs. In an in vitro co-culture experiment, salbutamol-treated DCs reduced the secretion of Th2 cytokine. In an in vivo AR animal experiment, salbutamol-administered mice showed reduced allergic behavior and Th2 cytokine expression in the nasal mucosa.

CONCLUSIONS

The regulation of ADRB2 with salbutamol alleviated the allergic response in vitro DC-T cell co-culture and in vivo AR animal models, suggesting that ADRB2 is a therapeutic target for AR and that ADRB2 agonists may be a promising medication for AR.

A comprehensive pharmacological analysis of fenoterol and its derivatives to unravel the role of β2-adrenergic receptor in zebrafish

β-adrenergic receptors (βARs) belong to a key molecular targets that regulate the most important processes occurring in the human organism. Although over the last decades a zebrafish model has been developed as a model complementary to rodents in biomedical research, the role of β₂AR in regulation of pathological and toxicological effects remains to elucidate. Therefore, the study aimed to clarify the role of β₂AR with a particular emphasis on the distinct role of subtypes A and B of zebrafish β₂AR. As model compounds selective β₂AR agonists - (R,R)-fenoterol ((R,R)-Fen) and its new derivatives: (R,R)-4'-methoxyfenoterol ((R,R)-MFen) and (R,R)-4'-methoxy-1-naphtylfenoterol ((R,R)-MNFen) - were tested. We described dose-dependent changes observed after fenoterols exposure in terms of general toxicity, cardiotoxicity and neurobehavioural responses. Subsequently, to better characterise the role of β₂-adrenergic stimulation in zebrafish, we have performed a series of molecular docking simulations. Our results indicate that (R,R)-Fen displays the highest affinity for subtype A of zebrafish β₂AR and β_2AAR might be involved in pigment depletion. (R,R)-MFen shows the lowest affinity for zebrafish β₂ARs out of the tested fenoterols and this might be associated with its cardiotoxic and anxiogenic effects. (R,R)-MNFen displays the highest affinity for subtype B of zebrafish β₂AR and modulation of this receptor might be associated with the development of malformations, increases locomotor activity and induces a negative chronotropic effect. Taken together, the presented data offer insights into the functional responses of the zebrafish β₂ARs confirming their intraspecies conservation, and support the translation of the zebrafish model in pharmacological and toxicological research.

The Nature of Functional Features of Different Classes of G-Protein-Coupled Receptors

G-protein-coupled receptors (GPCRs) are a critical family in the human proteome and are involved in various physiological processes. They are also the most important drug target, with approximately 30% of approved drugs acting on such receptors. The members of the family are divided into six classes based on their structural and functional characteristics. Understanding their structural-functional relationships will benefit us in future drug development. In this article, we investigate the features of protein function, structure, and energy that describe the dynamics of the GPCR activation process between different families. GPCRs straddle the cell membrane and transduce signals from outside the membrane into the cell. During the process, the conformational change in GPCRs that is activated by the binding of signal molecules is essential. During the binding process, different types of signal molecules result in different signal transfer efficiencies. Therefore, the GPCR classes show a variety of structures and activation processes. Based on the experimental crystal structures, we modeled the activation process of the β2 adrenergic receptor (β2AR), glucagon receptor (GCGR), and metabotropic glutamate receptor 2 (mGluR2), which represent class A, B, and C GPCRs, respectively. We calculated their activation free-energy landscapes and analyzed the structure-energy-function relationship. The results show a consistent picture of the activation mechanisms between different types of GPCRs. This could also provide us a way to understand other signal transduction proteins.

Mechanism of enhanced sensitivity of mutated β-adrenergic-like octopamine receptor to amitraz in honeybee Apis mellifera: An insight from MD simulations

BACKGROUND

Amitraz is one of the critical acaricides/insecticides for effective control of pest infestation of Varroa destructor mite, a devastating parasite of Apis mellifera, because of its low toxicity to honeybees. Previous assays verified that a typical G protein-coupled receptor, β-adrenergic-like octopamine receptor (Octβ2R), is the unique target of amitraz, but the honeybee Octβ2R resists to amitraz. However, the underlying molecular mechanism of the enhanced sensitivity or toxicity of amitraz to mutated honeybee Octβ2R^{E208V/I335T/I350V} is not fully understood. Here, molecular dynamics simulations are employed to explore the implied mechanism of the enhanced sensitivity to amitraz in mutant honeybee Octβ2R.

RESULTS

We found that amitraz binding stabilized the structure of Octβ2R, particularly the intracellular loop 3 associated with the Octβ2R signaling. Then, it was further demonstrated that both mutations and ligand binding resulted in a more rigid and compact amitraz binding site, as well as the outward movement of the transmembrane helix 6, which was a prerequisite for G protein coupling and activation. Moreover, mutations were found to promote the binding between Octβ2R and amitraz. Finally, community analysis illuminated that mutations and amitraz strengthened the residue-residue communication within the transmembrane domain, which might facilitate the allosteric signal propagation and activation of Octβ2R.

CONCLUSION

Our results unveiled structural determinants of improved sensitivity in the Octβ2R-amitraz complex and may contribute to further structure-based drug design for safer and less toxic selective insecticides. © 2022 Society of Chemical Industry.

Activation and Speciation Mechanisms in Class A GPCRs

Accurate development of allosteric modulators of GPCRs require a thorough assessment of their sequence, structure, and dynamics, toward gaining insights into their mechanisms of actions shared by family members, as well as dynamic features that distinguish subfamilies. Building on recent progress in the characterization of the signature dynamics of proteins, we analyzed here a dataset of 160 Class A GPCRs to determine their sequence similarities, structural landscape, and dynamic features across different species (human, bovine, mouse, squid, and rat), different activation states (active/inactive), and different subfamilies. The two dominant directions of variability across experimentally resolved structures, identified by principal component analysis of the dataset, shed light to cooperative mechanisms of activation, subfamily differentiation, and speciation of Class A GPCRs. The analysis reveals the functional significance of the conformational flexibilities of specific structural elements, including: the dominant role of the intracellular loop 3 (ICL3) together with the cytoplasmic ends of the adjoining helices TM5 and TM6 in enabling allosteric activation; the role of particular structural motifs at the extracellular loop 2 (ECL2) connecting TM4 and TM5 in binding ligands specific to different subfamilies; or even the differentiation of the N-terminal conformation across different species. Detailed analyses of the modes of motions accessible to the members of the dataset and their variations across members demonstrate how the active and inactive states of GPCRs obey distinct conformational dynamics. The collective fluctuations of the GPCRs are robustly defined in the active state, while the inactive conformers exhibit broad variance among members.

Structures of β1-adrenergic receptor in complex with Gs and ligands of different efficacies

G-protein-coupled receptors (GPCRs) receive signals from ligands with different efficacies, and transduce to heterotrimeric G-proteins to generate different degrees of physiological responses. Previous studies revealed how ligands with different efficacies activate GPCRs. Here, we investigate how a GPCR activates G-proteins upon binding ligands with different efficacies. We report the cryo-EM structures of β₁-adrenergic receptor (β₁-AR) in complex with Gs (Gα_sGβ₁Gγ₂) and a partial agonist or a very weak partial agonist, and compare them to the β₁-AR-Gs structure in complex with a full agonist. Analyses reveal similar overall complex architecture, with local conformational differences. Cellular functional studies with mutations of β₁-AR residues show effects on the cellular signaling from β₁-AR to the cAMP response initiated by the three different ligands, with residue-specific functional differences. Biochemical investigations uncover that the intermediate state complex comprising β₁-AR and nucleotide-free Gs is more stable when binding a full agonist than a partial agonist. Molecular dynamics simulations support the local conformational flexibilities and different stabilities among the three complexes. These data provide insights into the ligand efficacy in the activation of GPCRs and G-proteins.

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.

Evolutionary association of receptor-wide amino acids with G protein-coupling selectivity in aminergic GPCRs

Evolutionary analyses of aminergic G protein–coupled receptors reveal receptor-wide potential determinants of G protein–coupling selectivity.

G protein-coupled receptors (GPCRs) induce signal transduction pathways through coupling to four main subtypes of G proteins (G_s, G_i, G_q, and G_12/13), selectively. However, G protein selective activation mechanisms and residual determinants in GPCRs have remained obscure. Herein, we performed extensive phylogenetic analysis and identified specifically conserved residues for the aminergic receptors having similar coupling profiles. By integrating our methodology of differential evolutionary conservation of G protein–specific amino acids with structural analyses, we identified specific activation networks for G_s, G_i1, G_o, and G_q. To validate that these networks could determine coupling selectivity we further analyzed G_s-specific activation network and its association with G_s selectivity. Through molecular dynamics simulations, we showed that previously uncharacterized Glycine at position 7x41 plays an important role in receptor activation and it may determine G_s coupling selectivity by facilitating a larger TM6 movement. Finally, we gathered our results into a comprehensive model of G protein selectivity called “sequential switches of activation” describing three main molecular switches controlling GPCR activation: ligand binding, G protein selective activation mechanisms, and G protein contact.

Structural insights into the activation of somatostatin receptor 2 by cyclic SST analogues

The endogenous cyclic tetradecapeptide SST14 was reported to stimulate all five somatostatin receptors (SSTR1-5) for hormone release, neurotransmission, cell growth arrest and cancer suppression. Two SST14-derived short cyclic SST analogues (lanreotide or octreotide) with improved stability and longer lifetime were developed as drugs to preferentially activate SSTR2 and treat acromegalia and neuroendocrine tumors. Here, cryo-EM structures of the human SSTR2-Gi complex bound with SST14, octreotide or lanreotide were determined at resolutions of 2.85 Å, 2.97 Å, and 2.87 Å, respectively. Structural and functional analysis revealed that interactions between β-turn residues in SST analogues and transmembrane SSTR2 residues in the ligand-binding pocket are crucial for receptor binding and functional stimulation of the two SST14-derived cyclic octapeptides. Additionally, Q102^2.63, N276^6.55, and F294^7.35 could be responsible for the selectivity of lanreotide or octreotide for SSTR2 over SSTR1 or SSTR4. These results provide valuable insights into further rational development of SST analogue drugs targeting SSTR2.

The Multifaceted Role of GPCRs in Amyotrophic Lateral Sclerosis: A New Therapeutic Perspective?

Amyotrophic lateral sclerosis (ALS) is a degenerating disease involving the motor neurons, which causes a progressive loss of movement ability, usually leading to death within 2 to 5 years from the diagnosis. Much effort has been put into research for an effective therapy for its eradication, but still, no cure is available. The only two drugs approved for this pathology, Riluzole and Edaravone, are onlyable to slow down the inevitable disease progression. As assessed in the literature, drug targets such as protein kinases have already been extensively examined as potential drug targets for ALS, with some molecules already in clinical trials. Here, we focus on the involvement of another very important and studied class of biological entities, G protein-coupled receptors (GPCRs), in the onset and progression of ALS. This workaimsto give an overview of what has been already discovered on the topic, providing useful information and insights that can be used by scientists all around the world who are putting efforts into the fight against this very important neurodegenerating disease.

The Impact of the Secondary Binding Pocket on the Pharmacology of Class A GPCRs

G-protein coupled receptors (GPCRs) are considered important therapeutic targets due to their pathophysiological significance and pharmacological relevance. Class A receptors represent the largest group of GPCRs that gives the highest number of validated drug targets. Endogenous ligands bind to the orthosteric binding pocket (OBP) embedded in the intrahelical space of the receptor. During the last 10 years, however, it has been turned out that in many receptors there is secondary binding pocket (SBP) located in the extracellular vestibule that is much less conserved. In some cases, it serves as a stable allosteric site harbouring allosteric ligands that modulate the pharmacology of orthosteric binders. In other cases it is used by bitopic compounds occupying both the OBP and SBP. In these terms, SBP binding moieties might influence the pharmacology of the bitopic ligands. Together with others, our research group showed that SBP binders contribute significantly to the affinity, selectivity, functional activity, functional selectivity and binding kinetics of bitopic ligands. Based on these observations we developed a structure-based protocol for designing bitopic compounds with desired pharmacological profile.

Universal Properties and Specificities of the β2-Adrenergic Receptor-Gs Protein Complex Activation Mechanism Revealed by All-Atom Molecular Dynamics Simulations

G protein-coupled receptors (GPCRs) are transmembrane proteins of high pharmacological relevance. It has been proposed that their activity is linked to structurally distinct, dynamically interconverting functional states and the process of activation relies on an interconnecting network of conformational switches in the transmembrane domain. However, it is yet to be uncovered how ligands with different extents of functional effect exert their actions. According to our recent hypothesis, based on indirect observations and the literature data, the transmission of the external stimulus to the intracellular surface is accompanied by the shift of macroscopic polarization in the transmembrane domain, furnished by concerted movements of highly conserved polar motifs and the rearrangement of polar species. In this follow-up study, we have examined the β₂-adrenergic receptor (β₂AR) to see if our hypothesis drawn from an extensive study of the μ-opioid receptor (MOP) is fundamental and directly transferable to other class A GPCRs. We have found that there are some general similarities between the two receptors, in agreement with previous studies, and there are some receptor-specific differences that could be associated with different signaling pathways.

Rapidly identifying bioactive compounds from Zhisou oral liquid by immobilized receptor-based high-performance affinity chromatography

The identification of bioactive compounds in complex matrices remains a major challenge due to the lack of highly efficient and specific methods. This work developed an approach based on high-performance affinity chromatography to identify the potential antitussive compounds from Zhisou oral liquid . The main methods include the synthesis of immobilized beta2-adrenoceptor by a one-step method, the screening and identification of the potential bioactive compounds by the receptor column coupled with mass spectrometry, and the binding mechanism analysis of the compounds to the receptor by the in vivo experiment, injection amount dependent method and molecular simulation. We identified the potential bioactive compounds of Zhisou oral liquid as glycyrrhizic acid, platycodin D, tuberostemonine, and hesperidin. In vivo experiment showed that the combinational utilization of the four compounds was possible to present an equivalent antitussive effect to the formula. The docking results demonstrated that hydrogen bonds and Van der Waals forces were the main forces to drive the binding of the four compounds to beta2-adrenoceptor. We concluded that the four compounds are the effective components in Zhisou oral liquid. The proposed strategy is possible to provide an alternative for the development of highly efficient methods to pursue the bioactive compounds of complex matrices.

Ligands of Adrenergic Receptors: A Structural Point of View

Adrenergic receptors are G protein-coupled receptors for epinephrine and norepinephrine. They are targets of many drugs for various conditions, including treatment of hypertension, hypotension, and asthma. Adrenergic receptors are intensively studied in structural biology, displayed for binding poses of different types of ligands. Here, we summarized molecular mechanisms of ligand recognition and receptor activation exhibited by structure. We also reviewed recent advances in structure-based ligand discovery against adrenergic receptors.