Structural features of β2 adrenergic receptor: crystal structures and beyond

Early Events in β2AR Dimer Dynamics Mediated by Activation-Related Microswitches

G-Protein-Coupled Receptors (GPCRs) make up around 3-4% of the human genome and are the targets of one-third of FDA-approved drugs. GPCRs typically exist as monomers but also aggregate to form higher-order oligomers, including dimers. β2AR, a pharmacologically relevant GPCR, is known to be targeted for the treatment of asthma and cardiovascular diseases. The activation of β2AR at the dimer level remains under-explored. In the current study, molecular dynamics (MD) simulations have been performed to understand activation-related structural changes in β2AR at the dimer level. The transition from inactive to active and vice versa has been studied by starting the simulations in the apo, agonist-bound, and inverse agonist-bound β2AR dimers for PDB ID: 2RH1 and PDB ID: 3P0G, respectively. A cumulative total of around 21-μs simulations were performed. Residue-based distances, RMSD, and PCA calculations suggested that either of the one monomer attained activation-related features for the apo and agonist-bound β2AR dimers. The TM5 and TM6 helices within the two monomers were observed to be in significant variation in all the simulations. TM5 bulge and proximity of TM2 and TM7 helices may be contributing to one of the early events in activation. The dimeric interface between TM1 and helix 8 were observed to be well maintained in the apo and agonist-bound simulations. The presence of inverse agonists favored inactive features in both the monomers. These key features of activation known for monomers were observed to have an impact on β2AR dimers, thereby providing an insight into the oligomerization mechanism of GPCRs.

Differential Behavior of Conformational Dynamics in Active and Inactive States of Cannabinoid Receptor 1

Cannabinoid receptor 1 (CB1) is a G protein-coupled receptor that regulates critical physiological processes including pain, appetite, and cognition. Understanding the conformational dynamics of CB1 associated with transitions between inactive and active signaling states is imperative for developing targeted modulators. Using microsecond-level all-atom molecular dynamics simulations, we identified marked differences in the conformational ensembles of inactive and active CB1 in apo. The inactive state exhibited substantially increased structural heterogeneity and plasticity compared to the more rigidified active state in the absence of stabilizing ligands. Transmembrane helices TM3 and TM7 were identified as distinguishing factors modulating the state-dependent dynamics. TM7 displayed amplified fluctuations selectively in the inactive state simulations attributed to disruption of conserved electrostatic contacts anchoring it to surrounding helices in the active state. Additionally, we identified significant reorganizations in key salt bridge and hydrogen bond networks contributing to the CB1 activation/inactivation. For instance, D213-Y224 hydrogen bond and D184-K192 salt bridge showed marked rearrangements between the states. Collectively, these findings reveal the specialized role of TM7 in directing state-dependent CB1 dynamics through electrostatic switch mechanisms. By elucidating the intrinsic enhanced flexibility of inactive CB1, this study provides valuable insights into the conformational landscape enabling functional transitions. Our perspective advances understanding of CB1 activation mechanisms and offers opportunities for structure-based drug discovery targeting the state-specific conformational dynamics of this receptor.

Screening for promising multi-target bioactive components from Cortex Mori Radicis for the treatment of chronic cor pulmonale based on immobilized beta1-adrenergic receptor and beta2-adrenergic receptor chromatography

Cortex Morin Radicis (CMR) is the dried root bark of Morus alba. L. It has a variety of effects such as antibacterial, anti-tumour, treatment of cardiovascular diseases or upper respiratory tract disease and so on. The pursuit for drugs selected from Cortex Mori Radicis having improved therapeutic efficacy necessitates increasing research on new assays for screening bioactive compounds with multi-targets. In this work, we applied immobilized β1-AR and β2-AR as the stationary phase in chromatographic column to screen bioactive compounds from Cortex Morin Radicis. Specific ligands of the two receptors (e.g. esmolol, metoprolol, atenolol, salbutamol, methoxyphenamine, tulobuterol and clorprenaline) were utilized to characterize the specificity and bioactivity of the columns. We used high performance affinity chromatography coupled with ESI-MS to screen targeted compounds of Cortex Morin Radicis. By zonal elution, we identified morin as a bioactive compound simultaneously binding to β1-AR and β2-AR. The compound exhibited the association constants of 3.10 × 104 and 2.60 × 104 M-1 on the β1-AR and β2-AR column. On these sites, the dissociation rate constants were calculated to be 0.131 and 0.097 s-1. Molecular docking indicated that the binding of morin to the two receptors occurred on Asp200, Asp121, and Val122 of β1-AR, Asn312, Thr110, Asp113, Tyr316, Gly90, Phe193, Ile309, and Trp109 of β2-AR. Likewise, mulberroside C was identified as the bioactive compound binding to β2-AR. The association constants and dissociation rate constants were calculated to be 1.08 × 104 M-1 and 0.900 s-1. Molecular docking also indicated that mulberroside C could bind to β2-AR receptor on its agonist site. Taking together, we demonstrated that the chromatographic strategy to identify bioactive natural products based on the β1-AR and β2-AR immobilization, has potential for screening bioactive compounds with multi-targets from complex matrices including traditional Chinese medicines.

Engineered mini-G proteins block the internalization of cognate GPCRs and disrupt downstream intracellular signaling

Mini-G proteins are engineered, thermostable variants of Gα subunits designed to stabilize G protein-coupled receptors (GPCRs) in their active conformations. Because of their small size and ease of use, they are popular tools for assessing GPCR behaviors in cells, both as reporters of receptor coupling to Gα subtypes and for cellular assays to quantify compartmentalized signaling at various subcellular locations. Here, we report that overexpression of mini-G proteins with their cognate GPCRs disrupted GPCR endocytic trafficking and associated intracellular signaling. In cells expressing the Gαs-coupled GPCR glucagon-like peptide 1 receptor (GLP-1R), coexpression of mini-Gs, a mini-G protein derived from Gαs, blocked β-arrestin 2 recruitment and receptor internalization and disrupted endosomal GLP-1R signaling. These effects did not involve changes in receptor phosphorylation or lipid nanodomain segregation. Moreover, we found that mini-G proteins derived from Gαi and Gαq also inhibited the internalization of GPCRs that couple to them. Finally, we developed an alternative intracellular signaling assay for GLP-1R using a nanobody specific for active Gαs:GPCR complexes (Nb37) that did not affect GLP-1R internalization. Our results have important implications for designing methods to assess intracellular GPCR signaling.

A Study on the Robustness and Stability of Explainable Deep Learning in an Imbalanced Setting: The Exploration of the Conformational Space of G Protein-Coupled Receptors

G-protein coupled receptors (GPCRs) are transmembrane proteins that transmit signals from the extracellular environment to the inside of the cells. Their ability to adopt various conformational states, which influence their function, makes them crucial in pharmacoproteomic studies. While many drugs target specific GPCR states to exert their effects-thereby regulating the protein's activity-unraveling the activation pathway remains challenging due to the multitude of intermediate transformations occurring throughout this process, and intrinsically influencing the dynamics of the receptors. In this context, computational modeling, particularly molecular dynamics (MD) simulations, may offer valuable insights into the dynamics and energetics of GPCR transformations, especially when combined with machine learning (ML) methods and techniques for achieving model interpretability for knowledge generation. The current study builds upon previous work in which the layer relevance propagation (LRP) technique was employed to interpret the predictions in a multi-class classification problem concerning the conformational states of the β2-adrenergic (β2AR) receptor from MD simulations. Here, we address the challenges posed by class imbalance and extend previous analyses by evaluating the robustness and stability of deep learning (DL)-based predictions under different imbalance mitigation techniques. By meticulously evaluating explainability and imbalance strategies, we aim to produce reliable and robust insights.

Differential Behavior of Conformational Dynamics in Active and Inactive States of Cannabinoid Receptor 1

The cannabinoid receptor CB1 is a G protein-coupled receptor that regulates critical physiological processes including pain, appetite, and cognition. Understanding the conformational dynamics of CB1 associated with transitions between inactive and active signaling states is imperative for developing targeted modulators. Using microsecond-level all-atom molecular dynamics simulations, we identified marked differences in the conformational ensembles of inactive and active CB1 states in apo conditions. The inactive state exhibited substantially increased structural heterogeneity and plasticity compared to the more rigidified active state in the absence of stabilizing ligands. Transmembrane helices TM3 and TM7 were identified as distinguishing factors modulating the state-dependent dynamics. TM7 displayed amplified fluctuations selectively in the inactive state simulations attributed to disruption of conserved electrostatic contacts anchoring it to surrounding helices in the active state. Additionally, we identified significant reorganization of key salt bridge and hydrogen bond networks known to control CB1 activation between states. For instance, a conserved D213-Y224 hydrogen bond and D184-K192 salt bridge interactions showed marked rearrangements between the states. Collectively, these findings reveal the specialized role of TM7 in directing state-dependent CB1 dynamics through electrostatic switch mechanisms. By elucidating the intrinsic enhanced flexibility of inactive CB1, this study provides valuable insights into the conformational landscape enabling functional transitions. Our perspective advances understanding of CB1 activation mechanisms and offers opportunities for structure-based drug discovery targeting the state-specific conformational dynamics of this receptor.

Coevolution-Driven Method for Efficiently Simulating Conformational Changes in Proteins Reveals Molecular Details of Ligand Effects in the β2AR Receptor

With the advent of AI-powered structure prediction, the scientific community is inching closer to solving protein folding. An unresolved enigma, however, is to accurately, reliably, and deterministically predict alternative conformational states that are crucial for the function of, e.g., transporters, receptors, or ion channels where conformational cycling is innately coupled to protein function. Accurately discovering and exploring all conformational states of membrane proteins has been challenging due to the need to retain atomistic detail while enhancing the sampling along interesting degrees of freedom. The challenges include but are not limited to finding which degrees of freedom are relevant, how to accelerate the sampling along them, and then quantifying the populations of each micro- and macrostate. In this work, we present a methodology that finds relevant degrees of freedom by combining evolution and physics through machine learning and apply it to the conformational sampling of the β2 adrenergic receptor. In addition to predicting new conformations that are beyond the training set, we have computed free energy surfaces associated with the protein's conformational landscape. We then show that the methodology is able to quantitatively predict the effect of an array of ligands on the β2 adrenergic receptor activation through the discovery of new metastable states not present in the training set. Lastly, we also stake out the structural determinants of activation and inactivation pathway signaling through different ligands and compare them to functional experiments to validate our methodology and potentially gain further insights into the activation mechanism of the β2 adrenergic receptor.

β2-adrenergic signaling promotes higher-affinity B cells and antibodies

Stress-induced β2-adrenergic receptor (β2AR) activation in B cells increases IgG secretion; however, the impact of this activation on antibody affinity and the underlying mechanisms remains unclear. In the current study, we demonstrate that stress in mice following ovalbumin (OVA) or SARS-CoV-2 RBD immunization significantly increases both serum and surface-expressed IgG binding to the immunogen, while concurrently reducing surface IgG expression and B cell clonal expansion. These effects were abolished by pharmacological β2AR blocking or when the experiments were conducted in β2AR -/- mice. In the second part of our study, we used single B cell sorting to characterize the monoclonal antibodies (mAbs) generated following β2AR activation in cultured RBD-stimulated B cells from convalescent SARS-CoV-2 donors. Ex vivo β2AR activation increased the affinities of the produced anti-RBD mAbs by 100-fold compared to mAbs produced by the same donor control cultures. Consistent with the mouse experiments, β2AR activation reduced both surface IgG levels and the frequency of expanded clones. mRNA sequencing revealed a β2AR-dependent upregulation of the PI3K pathway and B cell receptor (BCR) signaling through AKT phosphorylation, as well as an increased B cell motility. Overall, our study demonstrates that stress-mediated β2AR activation drives changes in B cells associated with BCR activation and higher affinity antibodies.

Stepwise Frontal Analysis Coupled with Affinity Chromatography: A Fast and Reliable Method for Potential Ligand Isolation and Evaluation from Mahuang-Fuzi-Xixin Decoction

Mahuang-Fuzi-Xixin Decoction (MFXD) is widely used in the treatment of asthma, however, the functional components in the decoction targeting beta2-adrenoceptor (β₂ -AR) remain unclear. Herein, we immobilized the haloalkane dehalogenase (Halo)-tagged β₂ -AR on the 6-chlorocaproic acid-modified microspheres. Using the affinity stationary phase, the interactions of four ligands with the receptor were analyzed by stepwise frontal analysis. The association constants were (4.75±0.28)×10⁴  M^-1 for salbutamol, (2.93±0.15)×10⁴  M^-1 for terbutaline, (1.23±0.03)×10⁴  M^-1 for methoxyphenamine, (5.67±0.38)×10⁴  M^-1 for clorprenaline at high-affinity binding site, and (2.73±0.05)×10³  M^-1 at low-affinity binding site. These association constants showed the same rank order as the radioligand binding assay, demonstrating that immobilized β₂ -AR had capacity to screen bioactive compounds binding to the receptor while stepwise frontal analysis could predict their binding affinities. Application of the immobilized receptor in analysis of MFXD by chromatographic method revealed that ephedrine, aconifine, karakoline, and chasmanine were the bioactive compounds targeting β₂ -AR. Among them, ephedrine and chasmanine exhibited association constants of (2.94±0.02)×10⁴ M^-1 and (4.60±0.15)×10⁴  M^-1 to the receptor by stepwise frontal analysis. Molecular docking analysis demonstrated that ephedrine, chasmanine, and the other two compounds interact with β₂ -AR through the same pocket involving the key amino acids such as Asn312, Asp113, Phe289, Trp286, Tyr316, and Val114. As such, we reasoned that the four compounds dominate the therapeutic effect of MFXD against asthma through β₂ -AR mediating pathway. This work shed light on the potential of immobilized β₂ -AR for drug discovery and provided a valuable methodology for rapid screening.

Mass spectrometry captures biased signalling and allosteric modulation of a G-protein-coupled receptor

G-protein-coupled receptors signal through cognate G proteins. Despite the widespread importance of these receptors, their regulatory mechanisms for G-protein selectivity are not fully understood. Here we present a native mass spectrometry-based approach to interrogate both biased signalling and allosteric modulation of the β1-adrenergic receptor in response to various ligands. By simultaneously capturing the effects of ligand binding and receptor coupling to different G proteins, we probed the relative importance of specific interactions with the receptor through systematic changes in 14 ligands, including isoprenaline derivatives, full and partial agonists, and antagonists. We observed enhanced dynamics of the intracellular loop 3 in the presence of isoprenaline, which is capable of acting as a biased agonist. We also show here that endogenous zinc ions augment the binding in receptor-Gs complexes and propose a zinc ion-binding hotspot at the TM5/TM6 intracellular interface of the receptor-Gs complex. Further interrogation led us to propose a mechanism in which zinc ions facilitate a structural transition of the intermediate complex towards the stable state.

L-DOPA and Droxidopa: From Force Field Development to Molecular Docking into Human β2-Adrenergic Receptor

The increasing interest in the molecular mechanism of the binding of different agonists and antagonists to β₂-adrenergic receptor (β₂AR) inactive and active states has led us to investigate protein–ligand interactions using molecular docking calculations. To perform this study, the 3.2 Å X-ray crystal structure of the active conformation of human β₂AR in the complex with the endogenous agonist adrenaline has been used as a template for investigating the binding of two exogenous catecholamines to this adrenergic receptor. Here, we show the derivation of L-DOPA and Droxidopa OPLS all atom (AA) force field (FF) parameters via quantum mechanical (QM) calculations, molecular dynamics (MD) simulations in aqueous solutions of the two catecholamines and the molecular docking of both ligands into rigid and flexible β₂AR models. We observe that both ligands share with adrenaline similar experimentally observed binding anchor sites, which are constituted by Asp113/Asn312 and Ser203/Ser204/Ser207 side chains. Moreover, both L-DOPA and Droxidopa molecules exhibit binding affinities comparable to that predicted for adrenaline, which is in good agreement with previous experimental and computational results. L-DOPA and Droxidopa OPLS AA FFs have also been tested by performing MD simulations of these ligands docked into β₂AR proteins embedded in lipid membranes. Both hydrogen bonds and hydrophobic interaction networks observed over the 1 μs MD simulation are comparable with those derived from molecular docking calculations and MD simulations performed with the CHARMM FF.

Evaluation of Neurotropic Activity and Molecular Docking Study of New Derivatives of pyrano[4″,3″:4',5']pyrido[3',2':4,5]thieno[3,2-d]pyrimidines on the Basis of pyrano[3,4-c]pyridines

Background: Heterocyclic compounds and their fused analogs, which contain pharmacophore fragments such as pyridine, thiophene and pyrimidine rings, are of great interest due to their broad spectrum of biological activity. Chemical compounds containing two or more pharmacophore groups due to additional interactions with active receptor centers usually enhance biological activity and can even lead to a new type of activity. The search for new effective neurotropic drugs in the series of derivatives of heterocycles containing pharmacophore groups in organic, bioorganic and medical chemistry is a serious problem. Methods: Modern methodology of drugs involves synthesis, physicochemical study, molecular modeling and selection of active compounds through virtual screening and experimental evaluation of the biological activity of new chimeric compounds with pharmacophore fragments. For the synthesis of new compounds, classical organic methods were used and developed. For the evaluation of neurotropic activity of new synthesized compounds, some biological methods were used according to indicators characterizing anticonvulsant, sedative and antianxiety activity as well as side effects. For docking analysis, various soft ware packages and methods were used. Results: As a result of multistep reactions, 11 new, tri- and tetracyclic heterocyclic systems were obtained. The studied compounds exhibit protection against pentylenetetrazole (PTZ) seizures as well as some psychotropic effects. The biological assays evidenced that nine of the eleven studied compounds showed a high anticonvulsant activity by antagonism with pentylenetetrazole. The toxicity of the compounds is low, and they do not induce muscle relaxation in the studied doses. According to the study of psychotropic activity, it was found that the selected compounds have an activating behavior and anxiolytic effects on the “open field” and “elevated plus maze” (EPM) models. The data obtained indicate the anxiolytic (antianxiety) activity of the derivatives of tricyclic thieno[2,3-b]pyridines and tetracyclic pyridothieno[3,2-d]pyrimidin-8-ones, especially pronounced in compounds 3b–f and 4e. The studied compounds increase the latent time of first immobilization on the “forced swimming” (FS) model and exhibit antidepressant effects; compounds 3e and 3f especially exhibit these effects, similarly to diazepam. Docking studies revealed that compounds 3c and 4b bound tightly in the active site of γ-aminobutyric acid type A (GABA_A) receptors with a value of the scoring function that estimates free energy of binding (∆G) at −10.0 ± 5 kcal/mol. Compound 4e showed the best affinity ((∆G) at −11.0 ± 0.54 kcal/mol) and seems to be an inhibitor of serotonin (SERT) transporter. Compounds 3c–f and 4e practically bound with the groove of T4L of 5HT_1A and blocked it completely, while the best affinity observed was in compound 3f ((∆G) at −9.3 ± 0.46 kcal/mol). Conclusions: The selected compounds have an anticonvulsant, activating behavior and anxiolytic effects and at the same time exhibit antidepressant effects.

Targeting ADRB2 enhances sensitivity of non-small cell lung cancer to VEGFR2 tyrosine kinase inhibitors

Vascular Endothelial Growth Factor Receptor 2 (VEGFR2) tyrosine kinase inhibitors (TKIs) have achieved remarkable clinical progress in the treatment of non-small-cell lung cancer; however, resistance has limited their therapeutic efficacy. Therefore, understanding the mechanisms of VEGF-TKI and ICI resistance will help to develop effective treatment strategies for patients with advanced NSCLC. Our results suggested that treatment with VEGFR2-TKIs upregulated ADRB2 expression in NSCLC cells. Propranolol, a common ADRB2 antagonist, significantly enhanced the therapeutic effect of VEGFR2-TKIs by inhibiting the ADRB2 signaling pathway in NSCLC cells in vitro and in vivo. Mechanically, the treatment-induced ADRB2 upregulation and the enhancement of ADRB2/VEGFR2 interaction caused resistance to VEGFR2-TKIs in NSCLC. And the inhibition of the ADRB2/CREB/PSAT1 signaling pathway sensitized cells to VEGFR2-TKIs. We demonstrated that ADRB2 signaling is crucial in mediating resistance to VEGFR2-TKIs and provided a novel promising combinatory approach to enhance the antitumor effect of VEGFR2-TKIs in NSCLC combining with propranolol.

Allosteric Effect of Nanobody Binding on Ligand-Specific Active States of the β2 Adrenergic Receptor

Nanobody binding stabilizes G-protein-coupled receptors (GPCR) in a fully active state and modulates their affinity for bound ligands. However, the atomic-level basis for this allosteric regulation remains elusive. Here, we investigate the conformational changes induced by the binding of a nanobody (Nb80) on the active-like β2 adrenergic receptor (β2AR) via enhanced sampling molecular dynamics simulations. Dimensionality reduction analysis shows that Nb80 stabilizes structural features of the β2AR with an ∼14 Å outward movement of transmembrane helix 6 and a close proximity of transmembrane (TM) helices 5 and 7, and favors the fully active-like conformation of the receptor, independent of ligand binding, in contrast to the conditions under which no intracellular binding partner is bound, in which case the receptor is only stabilized in an intermediate-active state. This activation is supported by the residues located at hotspots located on TMs 5, 6, and 7, as shown by supervised machine learning methods. Besides, ligand-specific subtle differences in the conformations assumed by intracellular loop 2 and extracellular loop 2 are captured from the trajectories of various ligand-bound receptors in the presence of Nb80. Dynamic network analysis further reveals that Nb80 binding triggers tighter and stronger local communication networks between the Nb80 and the ligand-binding sites, primarily involving residues around ICL2 and the intracellular end of TM3, TM5, TM6, as well as ECL2, ECL3, and the extracellular ends of TM6 and TM7. In particular, we identify unique allosteric signal transmission mechanisms between the Nb80-binding site and the extracellular domains in conformations modulated by a full agonist, BI167107, and a G-protein-biased partial agonist, salmeterol, involving mainly TM1 and TM2, and TM5, respectively. Altogether, our results provide insights into the effect of intracellular binding partners on the GPCR activation mechanism, which should be taken into account in structure-based drug discovery.

Analysis of L-DOPA and droxidopa binding to human β2-adrenergic receptor

Over the last two decades, an increasing number of studies has been devoted to a deeper understanding of the molecular process involved in the binding of various agonists and antagonists to active and inactive conformations of β2-adrenergic receptor (β2AR). The 3.2 Å x-ray crystal structure of human β2AR active state in combination with the endogenous low affinity agonist adrenaline offers an ideal starting structure for studying the binding of various catecholamines to adrenergic receptors. We show that molecular docking of levodopa (L-DOPA) and droxidopa into rigid and flexible β2AR models leads for both ligands to binding anchor sites comparable to those experimentally reported for adrenaline, namely D113/N312 and S203/S204/S207 side chains. Both ligands have a hydrogen bond network that is extremely similar to those of noradrenaline and dopamine. Interestingly, redocking neutral and protonated versions of adrenaline to rigid and flexible β2AR models results in binding poses that are more energetically stable and distinct from the x-ray crystal structure. Similarly, lowest energy conformations of noradrenaline and dopamine generated by docking into flexible β2AR models had binding free energies lower than those of best poses in rigid receptor models. Furthermore, our findings show that L-DOPA and droxidopa molecules have binding affinities comparable to those predicted for adrenaline, noradrenaline, and dopamine, which are consistent with previous experimental and computational findings and supported by the molecular dynamics simulations of β2AR-ligand complexes performed here.

Functional solubilization of the β2-adrenoceptor using diisobutylene maleic acid

The β2-adrenoceptor (β2AR) is a well-established target in asthma and a prototypical G protein-coupled receptor for biophysical studies. Solubilization of membrane proteins has classically involved the use of detergents. However, the detergent environment differs from the native membrane environment and often destabilizes membrane proteins. Use of amphiphilic copolymers is a promising strategy to solubilize membrane proteins within their native lipid environment in the complete absence of detergents. Here we show the isolation of the β₂AR in the polymer diisobutylene maleic acid (DIBMA). We demonstrate that β₂AR remains functional in the DIBMA lipid particle and shows improved thermal stability compared with the n-dodecyl-β-D-maltopyranoside detergent-solubilized β₂AR. This unique method of extracting β₂AR offers significant advantages over previous methods routinely employed such as the introduction of thermostabilizing mutations and the use of detergents, particularly for functional biophysical studies.

Phosphodiesterase 4 mediates interleukin-8-induced heterologous desensitization of the β2 -adrenergic receptor

Acute respiratory distress syndrome (ARDS) is a life-threatening illness characterized by decreased alveolar-capillary barrier function, pulmonary edema consisting of proteinaceous fluid, and inhibition of net alveolar fluid transport responsible for resolution of pulmonary edema. There is currently no pharmacotherapy that has proven useful to prevent or treat ARDS, and two trials using beta-agonist therapy to treat ARDS demonstrated no effect. Prior studies indicated that IL-8-induced heterologous desensitization of the beta2-adrenergic receptor (β₂ -AR) led to decreased beta-agonist-induced mobilization of cyclic adenosine monophosphate (cAMP). Interestingly, phosphodiesterase (PDE) 4 inhibitors have been used in human airway diseases characterized by low intracellular cAMP levels and increases in specific cAMP hydrolyzing activity. Therefore, we hypothesized that PDE4 would mediate IL-8-induced heterologous internalization of the β₂ -AR and that PDE4 inhibition would restore beta-agonist-induced functions. We determined that CINC-1 (a functional IL-8 analog in rats) induces internalization of β₂ -AR from the cell surface, and arrestin-2, PDE4, and β₂ -AR form a complex during this process. Furthermore, we determined that cAMP associated with the plasma membrane was adversely affected by β₂ -AR heterologous desensitization. Additionally, we determined that rolipram, a PDE4 inhibitor, reversed CINC-1-induced derangements of cAMP and also caused β₂ -AR to successfully recycle back to the cell surface. Finally, we demonstrated that rolipram could reverse CINC-1-mediated inhibition of beta-agonist-induced alveolar fluid clearance in a murine model of trauma-shock. These results indicate that PDE4 plays a role in CINC-1-induced heterologous internalization of the β₂ -AR; PDE4 inhibition reverses these effects and may be a useful adjunct in particular ARDS patients.

Predicting aptamer sequences that interact with target proteins using an aptamer-protein interaction classifier and a Monte Carlo tree search approach

Oligonucleotide-based aptamers, which have a three-dimensional structure with a single-stranded fragment, feature various characteristics with respect to size, toxicity, and permeability. Accordingly, aptamers are advantageous in terms of diagnosis and treatment and are materials that can be produced through relatively simple experiments. Systematic evolution of ligands by exponential enrichment (SELEX) is one of the most widely used experimental methods for generating aptamers; however, it is highly expensive and time-consuming. To reduce the related costs, recent studies have used in silico approaches, such as aptamer-protein interaction (API) classifiers that use sequence patterns to determine the binding affinity between RNA aptamers and proteins. Some of these methods generate candidate RNA aptamer sequences that bind to a target protein, but they are limited to producing candidates of a specific size. In this study, we present a machine learning approach for selecting candidate sequences of various sizes that have a high binding affinity for a specific sequence of a target protein. We applied the Monte Carlo tree search (MCTS) algorithm for generating the candidate sequences using a score function based on an API classifier. The tree structure that we designed with MCTS enables nucleotide sequence sampling, and the obtained sequences are potential aptamer candidates. We performed a quality assessment using the scores of docking simulations. Our validation datasets revealed that our model showed similar or better docking scores in ZDOCK docking simulations than the known aptamers. We expect that our method, which is size-independent and easy to use, can provide insights into searching for an appropriate aptamer sequence for a target protein during the simulation step of SELEX.

Optimized iLID Membrane Anchors for Local Optogenetic Protein Recruitment

Optogenetic protein dimerization systems are powerful tools to investigate the biochemical networks that cells use to make decisions and coordinate their activities. These tools, including the improved Light-Inducible Dimer (iLID) system, offer the ability to selectively recruit components to subcellular locations, such as micron-scale regions of the plasma membrane. In this way, the role of individual proteins within signaling networks can be examined with high spatiotemporal resolution. Currently, consistent recruitment is limited by heterogeneous optogenetic component expression, and spatial precision is diminished by protein diffusion, especially over long time scales. Here, we address these challenges within the iLID system with alternative membrane anchoring domains and fusion configurations. Using live cell imaging and mathematical modeling, we demonstrate that the anchoring strategy affects both component expression and diffusion, which in turn impact recruitment strength, kinetics, and spatial dynamics. Compared to the commonly used C-terminal iLID fusion, fusion proteins with large N-terminal anchors show stronger local recruitment, slower diffusion of recruited components, efficient recruitment over wider gene expression ranges, and improved spatial control over signaling outputs. We also define guidelines for component expression regimes for optimal recruitment for both cell-wide and subcellular recruitment strategies. Our findings highlight key sources of imprecision within light-inducible dimer systems and provide tools that allow greater control of subcellular protein localization across diverse cell biological applications.

The role of ADRB2 gene polymorphisms in malignancies

Beta-2-adrenergic receptor is a member of the G protein-coupled receptor superfamily, which is highly expressed in most malignancies. There is increasing evidence showing that beta-2-adrenergic receptors are associated with carcinogenesis, proliferation, immune regulation, invasion, angiogenesis, clinical prognosis and treatment resistance in malignancies. Polymorphisms of the ADRB2 gene have been confirmed to be associated with transcriptional activity, mRNA translation, and beta-2-adrenergic receptor expression and sensitivity. This review discusses clinically relevant examples of single nucleotide polymorphisms of ADRB2 in malignancies and the effects of these polymorphisms on cancer susceptibility, prognosis and treatment response of cancer patients.

Chirality of β2-agonists. An overview of pharmacological activity, stereoselective analysis, and synthesis

β2-Agonists (β2-adrenergic agonists, bronchodilatants, and sympathomimetic drugs) are a group of drugs that are mainly used in asthma and obstructive pulmonary diseases. In practice, the substances used to contain one or more stereogenic centers in their structure and their enantiomers exhibit different pharmacological properties. In terms of bronchodilatory activity, (R)-enantiomers showed higher activity. The investigation of stereoselectivity in action and disposition of chiral drugs together with the preparation of pure enantiomer drugs calls for efficient stereoselective analytical methods. The overview focuses on the stereoselectivity in pharmacodynamics and pharmacokinetics of β2-agonists and summarizes the stereoselective analytical methods for the enantioseparation of racemic beta-agonists (HPLC, LC-MS, GC, TLC, CE). Some methods of the stereoselective synthesis for β2-agonists preparation are also presented.

Preferential Coupling of Dopamine D2S and D2L Receptor Isoforms with Gi1 and Gi2 Proteins-In Silico Study

The dopamine D₂ receptor belongs to rhodopsin-like G protein-coupled receptors (GPCRs) and it is an important molecular target for the treatment of many disorders, including schizophrenia and Parkinson's disease. Here, computational methods were used to construct the full models of the dopamine D₂ receptor short (D_2S) and long (D_2L) isoforms (differing with 29 amino acids insertion in the third intracellular loop, ICL3) and to study their coupling with G_i1 and G_i2 proteins. It was found that the D_2L isoform preferentially couples with the G_i2 protein and D_2S isoform with the G_i1 protein, which is in accordance with experimental data. Our findings give mechanistic insight into the interplay between isoforms of dopamine D₂ receptors and G_i proteins subtypes, which is important to understand signaling by these receptors and their mediation by pharmaceuticals, in particular psychotic and antipsychotic agents.

The β2-adrenergic receptor-ROS signaling axis: An overlooked component of β2AR function?

β2-Adrenergic receptor (β2AR) agonists are clinically used to elicit rapid bronchodilation for the treatment of bronchospasms in pulmonary diseases such as asthma and COPD, both of which exhibit characteristically high levels of reactive oxygen species (ROS); likely secondary to over-expression of ROS generating enzymes and chronically heightened inflammation. Interestingly, β2AR has long-been linked to ROS, yet the involvement of ROS in β2AR function has not been as vigorously studied as other aspects of β2AR signaling. Herein, we discuss the existing body of evidence linking β2AR activation to intracellular ROS generation and importantly, the role of ROS in regulating β2AR function. The reciprocal interplay of the β2AR and ROS appear to endow this receptor with the ability to self-regulate signaling efficacy and ligand binding, hereby unveiling a redox-axis that may be unfavorably altered in pathological states contributing to both disease progression and therapeutic drug responses.

Detailed insight on β-adrenoceptors as therapeutic targets

Human G protein-coupled receptors (GPCRs), especially adrenoceptors, play a crucial role in maintaining important physiological activities including cardiovascular and pulmonary functions. Among all adrenoceptors, β-adrenoceptors are the best characterized GPCRs and possess distinctive features as drug targets. Similarly, ligands that activate/deactivate β-adrenoceptors also hold a significant position in the field of biomarker identification and drug discovery. Several studies regarding molecular characterization of the β-adrenoceptor ligands have revealed that ligands with abilities to inhibit basal or intrinsic receptor activity or prevent receptor desensitization are particularly important to efficiently manage detrimental health conditions, including chronic heart failure, asthma, chronic obstructive pulmonary disease, obesity, and diabetes. Given the importance of β-adrenoceptors as molecular targets for many pathological conditions, this review aims to provide a detailed insight on the structural and functional aspects of β-adrenoceptors, with a particular emphasis on their importance as biomarkers and therapeutic targets.

PM2.5-induced ADRB2 hypermethylation contributed to cardiac dysfunction through cardiomyocytes apoptosis via PI3K/Akt pathway

BACKGROUND

Long-term exposure to fine particulate matter (PM_2.5) can causally contribute to progression of atherosclerosis, risk of ischemic heart disease and death, but the underlying mechanism is little known. Since DNA methylation impacts the process of heart disease, it might be useful in exploring potential mechanistic pathways linking PM_2.5 exposure and heart disease.

OBJECTIVES

Here, we investigated the PM_2.5-induced ADRB2 hypermethylation and the involving epigenetic mechanism of PM_2.5-induced cardiomyocytes apoptosis and cardiac dysfunction.

METHODS AND RESULTS

In vitro, PM_2.5 markedly augmented cardiotoxicity including oxidative damage and apoptosis in cardiomyocytes AC16 as well as epigenetic alteration. DNA methylation profiling revealed a significant gene-ADRB2 was involved in the cardiac relative GO and KEGG pathways. Methylation chip and Bisulfite Sequencing PCR (BSP) both identified the hypermethylation status of ADRB2 which encodes β2-Adrenergic receptor (β2AR). Mechanistic study showed ADRB2 hypermethylation-induced down-regulation of β2AR inhibited PI3K/Akt and then activated Bcl-2/BAX and p53 pathway in AC16. The transgenic cell lines showed over-expression of ADRB2 weakened the PM_2.5-induced cardiomyocytes apoptosis in opposite way, but was augmented by PI3K inhibitor (LY294002). In vivo, echocardiography showed the heart contractile function was decreased after SD rats intratracheal instillation of PM_2.5 for 30 days. The myocardial interstitial edema, myocardial gap expansion and myofibril disorder in PM_2.5 treated group were observed in rats heart tissue. What's more, basal expression of β2AR and VEGFR2 decreased in heart tissue as the dosage of PM_2.5 increasing, meanwhile PM_2.5 markedly attenuated PI3K/Akt pathway followed by augmented Bcl-2/BAX and p53 pathway, thus caused a greater number of TUNEL positive cardiomyocytes resulted in cardiac dysfunction in vivo.

CONCLUSIONS

PM_2.5 exposure could cause the myocardial ADRB2 hypermethylation and activate the β2AR/PI3K/Akt pathway, resulted in PM_2.5-induced cardiomyocytes apoptosis and cardiac dysfunction. Our study suggested that the ADRB2 demethylation or ADRB2/β2AR activation may serve as a potential pathway to prevent cardiac dysfunction induced by PM_2.5 exposure.

Mammalian Eps15 homology domain 1 potentiates angiogenesis of non-small cell lung cancer by regulating β2AR signaling

BACKGROUND

Non-small cell lung cancer (NSCLC) is a devastating disease with a heterogeneous prognosis, and the molecular mechanisms underlying tumor progression remain elusive. Mammalian Eps15 homology domain 1 (EHD1) plays a promotive role in tumor progression, but its role in cancer angiogenesis remains unknown. This study thus explored the role of EHD1 in angiogenesis in NSCLC.

METHODS

The changes in angiogenesis were evaluated through human umbilical vein endothelial cell (HUVEC) proliferation, migration and tube formation assays. The impact of EHD1 on β2-adrenoceptor (β2AR) signaling was evaluated by Western blotting, quantitative real-time polymerase chain reaction (qRT-PCR) analysis, and enzyme-linked immunosorbent assay (ELISA). The interaction between EHD1 and β2AR was confirmed by immunofluorescence (IF) and coimmunoprecipitation (Co-IP) experiments, and confocal microscopy immunofluorescence studies revealed that β2AR colocalized with the recycling endosome marker Rab11, which indicated β2AR endocytosis. Xenograft tumor models were used to investigate the role of EHD1 in NSCLC tumor growth.

RESULTS

The microarray analysis revealed that EHD1 was significantly correlated with tumor angiogenesis, and loss- and gain-of-function experiments demonstrated that EHD1 potentiates HUVEC proliferation, migration and tube formation. EHD1 knockdown inhibited β2AR signaling activity, and EHD1 upregulation promoted vascular endothelial growth factor A (VEGFA) and β2AR expression. Interestingly, EHD1 interacted with β2AR and played a novel and critical role in β2AR endocytic recycling to prevent receptor degradation. Aberrant VEGFA or β2AR expression significantly affected EHD1-mediated tumor angiogenesis. The proangiogenic role of EHD1 was confirmed in xenograft tumor models, and immunohistochemistry (IHC) analysis confirmed that EHD1 expression was positively correlated with VEGFA expression, microvessel density (MVD) and β2AR expression in patient specimens.

CONCLUSION

Collectively, the data obtained in this study suggest that EHD1 plays a critical role in NSCLC angiogenesis via β2AR signaling and highlight a potential target for antiangiogenic therapy.

Born This Way: Using Intrinsic Disorder to Map the Connections between SLITRKs, TSHR, and Male Sexual Orientation

Recently, genome-wide association study reveals a significant association between specific single nucleotide polymorphisms (SNPs) in men and their sexual orientation. These SNPs (rs9547443 and rs1035144) reside in the intergenic region between the SLITRK5 and SLITRK6 genes and in the intronic region of the TSHR gene and might affect functionality of SLITRK5, SLITRK6, and TSHR proteins that are engaged in tight control of key developmental processes, such as neurite outgrowth and modulation, cellular differentiation, and hormonal regulation. SLITRK5 and SLITRK6 are single-pass transmembrane proteins, whereas TSHR is a heptahelical G protein-coupled receptor (GPCR). Mutations in these proteins are associated with various diseases and are linked to phenotypes found at a higher rate in homosexual men. A bioinformatics analysis of SLITRK5, SLITRK6, and TSHR proteins is conducted to look at their structure, protein interaction networks, and propensity for intrinsic disorder. It is assumed that this information might improve understanding of the roles that SLITRK5, SLITRK6, and TSHR play within neuronal and thyroidal tissues and give insight into the phenotypes associated with male homosexuality.

Biophysical and functional characterization of Norrin signaling through Frizzled4

Wnt signaling is initiated by Wnt ligand binding to the extracellular ligand binding domain, called the cysteine-rich domain (CRD), of a Frizzled (Fzd) receptor. Norrin, an atypical Fzd ligand, specifically interacts with Fzd4 to activate β-catenin-dependent canonical Wnt signaling. Much of the molecular basis that confers Norrin selectivity in binding to Fzd4 was revealed through the structural study of the Fzd4_CRD-Norrin complex. However, how the ligand interaction, seemingly localized at the CRD, is transmitted across full-length Fzd4 to the cytoplasm remains largely unknown. Here, we show that a flexible linker domain, which connects the CRD to the transmembrane domain, plays an important role in Norrin signaling. The linker domain directly contributes to the high-affinity interaction between Fzd4 and Norrin as shown by ∼10-fold higher binding affinity of Fzd4_CRD to Norrin in the presence of the linker. Swapping the Fzd4 linker with the Fzd5 linker resulted in the loss of Norrin signaling, suggesting the importance of the linker in ligand-specific cellular response. In addition, structural dynamics of Fzd4 associated with Norrin binding investigated by hydrogen/deuterium exchange MS revealed Norrin-induced conformational changes on the linker domain and the intracellular loop 3 (ICL3) region of Fzd4. Cell-based functional assays showed that linker deletion, L430A and L433A mutations at ICL3, and C-terminal tail truncation displayed reduced β-catenin-dependent signaling activity, indicating the functional significance of these sites. Together, our results provide functional and biochemical dissection of Fzd4 in Norrin signaling.

PtdIns(4,5)P2 stabilizes active states of GPCRs and enhances selectivity of G-protein coupling

G-protein-coupled receptors (GPCRs) are involved in many physiological processes and are therefore key drug targets1. Although detailed structural information is available for GPCRs, the effects of lipids on the receptors, and on downstream coupling of GPCRs to G proteins are largely unknown. Here we use native mass spectrometry to identify endogenous lipids bound to three class A GPCRs. We observed preferential binding of phosphatidylinositol-4,5-bisphosphate (PtdIns(4,5)P2) over related lipids and confirm that the intracellular surface of the receptors contain hotspots for PtdIns(4,5)P2 binding. Endogenous lipids were also observed bound directly to the trimeric Gαsβγ protein complex of the adenosine A2A receptor (A2AR) in the gas phase. Using engineered Gα subunits (mini-Gαs, mini-Gαi and mini-Gα12)2, we demonstrate that the complex of mini-Gαs with the β1 adrenergic receptor (β1AR) is stabilized by the binding of two PtdIns(4,5)P2 molecules. By contrast, PtdIns(4,5)P2 does not stabilize coupling between β1AR and other Gα subunits (mini-Gαi or mini-Gα12) or a high-affinity nanobody. Other endogenous lipids that bind to these receptors have no effect on coupling, highlighting the specificity of PtdIns(4,5)P2. Calculations of potential of mean force and increased GTP turnover by the activated neurotensin receptor when coupled to trimeric Gαiβγ complex in the presence of PtdIns(4,5)P2 provide further evidence for a specific effect of PtdIns(4,5)P2 on coupling. We identify key residues on cognate Gα subunits through which PtdIns(4,5)P2 forms bridging interactions with basic residues on class A GPCRs. These modulating effects of lipids on receptors suggest consequences for understanding function, G-protein selectivity and drug targeting of class A GPCRs.

Identification of serotonin 2A receptor as a novel HCV entry factor by a chemical biology strategy

Hepatitis C virus (HCV) is a leading cause of liver disease worldwide. Although several HCV protease/polymerase inhibitors were recently approved by U.S. FDA, the combination of antivirals targeting multiple processes of HCV lifecycle would optimize anti-HCV therapy and against potential drug-resistance. Viral entry is an essential target step for antiviral development, but FDA-approved HCV entry inhibitor remains exclusive. Here we identify serotonin 2A receptor (5-HT_2AR) is a HCV entry factor amendable to therapeutic intervention by a chemical biology strategy. The silencing of 5-HT_2AR and clinically available 5-HT_2AR antagonist suppress cell culture-derived HCV (HCVcc) in different liver cells and primary human hepatocytes at late endocytosis process. The mechanism is related to regulate the correct plasma membrane localization of claudin 1 (CLDN1). Moreover, phenoxybenzamine (PBZ), an FDA-approved 5-HT_2AR antagonist, inhibits all major HCV genotypes in vitro and displays synergy in combination with clinical used anti-HCV drugs. The impact of PBZ on HCV genotype 2a is documented in immune-competent humanized transgenic mice. Our results not only expand the understanding of HCV entry, but also present a promising target for the invention of HCV entry inhibitor.

Structural insights into positive and negative allosteric regulation of a G protein-coupled receptor through protein-lipid interactions

Lipids are becoming known as essential allosteric modulators of G protein-coupled receptor (GPCRs). However, how they exert their effects on GPCR conformation at the atomic level is still unclear. In light of recent experimental data, we have performed several long-timescale molecular dynamics (MD) simulations, totalling 24 μs, to rigorously map allosteric modulation and conformational changes in the β₂ adrenergic receptor (β2AR) that occur as a result of interactions with three different phospholipids. In particular, we identify different sequential mechanisms behind receptor activation and deactivation, respectively, mediated by specific lipid interactions with key receptor regions. We show that net negatively charged lipids stabilize an active-like state of β2AR that is able to dock G_sα protein. Clustering of anionic lipids around the receptor with local distortion of membrane thickness is also apparent. On the other hand, net-neutral zwitterionic lipids inactivate the receptor, generating either fully inactive or intermediate states, with kinetics depending on lipid headgroup charge distribution and hydrophobicity. These chemical differences alter membrane thickness and density, which differentially destabilize the β2AR active state through lateral compression effects.

Ligand binding modes from low resolution GPCR models and mutagenesis: chicken bitter taste receptor as a test-case

Bitter taste is one of the basic taste modalities, warning against consuming potential poisons. Bitter compounds activate members of the bitter taste receptor (Tas2r) subfamily of G protein-coupled receptors (GPCRs). The number of functional Tas2rs is species-dependent. Chickens represent an intriguing minimalistic model, because they detect the bitter taste of structurally different molecules with merely three bitter taste receptor subtypes. We investigated the binding modes of several known agonists of a representative chicken bitter taste receptor, ggTas2r1. Because of low sequence similarity between ggTas2r1 and crystallized GPCRs (~10% identity, ~30% similarity at most), the combination of computational approaches with site-directed mutagenesis was used to characterize the agonist-bound conformation of ggTas2r1 binding site between TMs 3, 5, 6 and 7. We found that the ligand interactions with N93 in TM3 and/or N247 in TM5, combined with hydrophobic contacts, are typically involved in agonist recognition. Next, the ggTas2r1 structural model was successfully used to identify three quinine analogues (epiquinidine, ethylhydrocupreine, quinidine) as new ggTas2r1 agonists. The integrated approach validated here may be applicable to additional cases where the sequence identity of the GPCR of interest and the existing experimental structures is low.

A conformational study of protonated noradrenaline by UV–UV and IR dip double resonance laser spectroscopy combined with an electrospray and a cold ion trap method

In protonated noradrenaline, 3 folded and 2 extended conformers were identified under the ultra-cold condition.

Allergy-specific Phenome-Wide Association Study for Immunogenes in Turkish Children

To dissect the role of immunogenetics in allergy and asthma, we performed a phenome-wide association study in 974 Turkish children selected from a cross-sectional study conducted using ISAAC (International Study of Asthma and Allergies in Children) Phase II tools. We investigated 9 loci involved in different immune functions (ADAM33, ADRB2, CD14, IL13, IL4, IL4R, MS4A2, SERPINE1, and TNF) with respect to 116 traits assessed through blood tests, hypertonic saline challenge tests, questionnaires, and skin prick tests. Multiple associations were observed for ADAM33: rs2280090 was associated with reduced MEF240% (i.e., the ratio of Mean Expiratory Flow after 240s of hypertonic saline inhalation with respect to the age- and ancestry-matched reference value) and with an increased risk of allergic bronchitis (p = 1.77*10⁻⁴ and p = 7.94*10⁻⁴, respectively); rs3918396 was associated with wheezing and eczema comorbidity (p = 3.41*10⁻⁴). IL4 rs2243250 was associated with increased FEV240 (Forced Expiratory Flow Volume after 240s of hypertonic saline inhalation; p = 4.81*10⁻⁴) and CD14 rs2569190 was associated with asthma diagnosis (p = 1.36*10⁻³). ADAM33 and IL4 appeared to play a role in the processes linked to allergic airway inflammation and lung function. Due to its association with wheezing and eczema comorbidity, ADAM33 may also be involved in the atopic march.