The structural basis for agonist and partial agonist action on a β(1)-adrenergic receptor

Revealing the graded activation mechanism of neurotensin receptor 1

Graded activation contributes to the precise regulation of GPCR activity, presenting new opportunities for drug design. In this work, a total of 10 μs enhanced-sampling simulations are performed to provide molecular insights into the binding dynamics differences of the neurotensin receptor 1 (NTSR1) to the full agonist SRI-9829, partial agonist RTI-3a and inverse agonist SR48692. The possible graded activation mechanism of NTSR1 is revealed by an integrated analysis utilizing the reweighted potential of mean force (PMF), deep learning (DL) and transfer entropy (TE). Specifically, the orthosteric pocket is observed to undergo expansion and contraction, with the G-protein-binding site experiencing interconversions among the inactive, intermediate and active-like states. Detailed structural comparisons capture subtle conformational differences arising from ligand binding in allosteric signaling, which can well explain the graded activation. Critical microswitches that contribute to graded activation are efficiently identified with the DL model. TE calculations enable the visualization of allosteric communication networks within the receptor, elucidating the driver-responder relationships associated with signal transduction. Fortunately, the dissociation of the full agonist from the orthosteric pocket is observed. The current findings systematically reveal the mechanism of NTSR1 graded activation, and also provide implications for structure-based drug design.

Computational Methods for the Discovery and Optimization of TAAR1 and TAAR5 Ligands

G-protein-coupled receptors (GPCRs) represent a family of druggable targets when treating several diseases and continue to be a leading part of the drug discovery process. Trace amine-associated receptors (TAARs) are GPCRs involved in many physiological functions with TAAR1 having important roles within the central nervous system (CNS). By using homology modeling methods, the responsiveness of TAAR1 to endogenous and synthetic ligands has been explored. In addition, the discovery of different chemo-types as selective murine and/or human TAAR1 ligands has helped in the understanding of the species-specificity preferences. The availability of TAAR1-ligand complexes sheds light on how different ligands bind TAAR1. TAAR5 is considered an olfactory receptor but has specific involvement in some brain functions. In this case, the drug discovery effort has been limited. Here, we review the successful computational efforts developed in the search for novel TAAR1 and TAAR5 ligands. A specific focus on applying structure-based and/or ligand-based methods has been done. We also give a perspective of the experimental data available to guide the future drug design of new ligands, probing species-specificity preferences towards more selective ligands. Hints for applying repositioning approaches are also discussed.

Deconvoluting Monomer- and Dimer-Specific Distance Distributions between Spin Labels in a Monomer/Dimer Mixture Using T1-Edited DEER EPR Spectroscopy

Double electron-electron resonance (DEER) EPR is a powerful tool in structural biology, providing distances between pairs of spin labels. When the sample consists of a mixture of oligomeric species (e.g., monomer and dimer), the question arises as to how to assign the peaks in the DEER-derived probability distance distribution to the individual species. Here, we propose incorporating an EPR longitudinal electron relaxation (T1) inversion recovery experiment within a DEER pulse sequence to resolve this problem. The apparent T1 between dipolar coupled electron spins measured from the inversion recovery time (τinv) dependence of the peak intensities in the T1-edited DEER-derived probability P(r) distance distribution will be affected by the number of nitroxide labels attached to the biomolecule of interest, for example, two for a monomer and four for a dimer. We show that global fitting of all the T1-edited DEER echo curves, recorded over a range of τinv values, permits the deconvolution of distances between spin labels originating from monomeric (longer T1) and dimeric (shorter T1) species. This is especially useful when the trapping of spin labels in different conformational states during freezing gives rise to complex P(r) distance distributions. The utility of this approach is demonstrated for two systems, the β1 adrenergic receptor and a construct of the huntingtin exon-1 protein fused to the immunoglobulin domain of protein G, both of which exist in a monomer-dimer equilibrium.

The Discovery of Novel α2a Adrenergic Receptor Agonists Only Coupling to Gαi/O Proteins by Virtual Screening

Most α2-AR agonists derived from dexmedetomidine have few structural differences between them and have no selectivity for α2A/2B-AR or Gi/Gs, which can lead to side effects in drugs. To obtain novel and potent α2A-AR agonists, we performed virtual screening for human α2A-AR and α2B-AR to find α2A-AR agonists with higher selectivity. Compound P300-2342 and its three analogs significantly decreased the locomotor activity of mice (p < 0.05). Furthermore, P300-2342 and its three analogs inhibited the binding of [3H] Rauwolscine with IC50 values of 7.72 ± 0.76 and 12.23 ± 0.11 μM, respectively, to α2A-AR and α2B-AR. In α2A-AR-HEK293 cells, P300-2342 decreased forskolin-stimulated cAMP production without increasing cAMP production, which indicated that P300-2342 activated α2A-AR with coupling to the Gαi/o pathway but without Gαs coupling. P300-2342 exhibited no agonist but slight antagonist activities in α2B-AR. Similar results were obtained for the analogs of P300-2342. The docking results showed that P300-2342 formed π-hydrogen bonds with Y394, V114 in α2A-AR, and V93 in α2B-AR. Three analogs of P300-2342 formed several π-hydrogen bonds with V114, Y196, F390 in α2A-AR, and V93 in α2B-AR. We believe that these molecules can serve as leads for the further optimization of α2A-AR agonists with potentially few side effects.

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.

Orphan receptor GPR88 as a potential therapeutic target for CNS disorders - an in silico approach

The G-protein-coupled receptors are a part of the largest and most physiologically relevant family of membrane proteins. One-third of the medications, now on the market, target the GPCR receptor family, which is one of the most important therapeutic targets for many disorders. In the reported work, we have focussed on orphan GPR88 receptor which is a part of the GPCR protein family and a potential target for central nervous system disorders. GPR88 is known to show the highest expression in the striatum, which is a key region in motor control and cognitive functions. Recent studies have reported that GPR88 is activated by two agonists, 2-PCCA and RTI-13951-33. In this study, we have predicted the three-dimensional protein structure for the orphan GPR88 using the homology modeling approach. We then used shape-based screening techniques based on known agonists and structure-based virtual screening methods employing docking to uncover novel GPR88 ligands. The screened GPR88-ligand complexes were further subjected to molecular dynamics simulation studies. The selected ligands could fasten the development of novel treatments for the vast list of movement and central nervous system disorders.Communicated by Ramaswamy H. Sarma.

Structural Insights into Partial Activation of the Prototypic G Protein-Coupled Adenosine A2A Receptor

The adenosine A2A receptor (A2AAR) belongs to the rhodopsin-like G protein-coupled receptor (GPCR) family, which constitutes the largest class of GPCRs. Partial agonists show reduced efficacy as compared to physiological agonists and can even act as antagonists in the presence of a full agonist. Here, we determined an X-ray crystal structure of the partial A2AAR agonist 2-amino-6-[(1H-imidazol-2-ylmethyl)sulfanyl]-4-p-hydroxyphenyl-3,5-pyridinedicarbonitrile (LUF5834) in complex with the A2AAR construct A2A-PSB2-bRIL, stabilized in its inactive conformation and being devoid of any mutations in the ligand binding pocket. The determined high-resolution structure (2.43 Å) resolved water networks and crucial binding pocket interactions. A direct hydrogen bond of the p-hydroxy group of LUF5834 with T883.36 was observed, an amino acid that was mutated to alanine in the most frequently used A2AAR crystallization constructs thus preventing the discovery of its interactions in most of the previous A2AAR co-crystal structures. G protein dissociation studies confirmed partial agonistic activity of LUF5834 as compared to that of the full agonist N-ethylcarboxamidoadenosine (NECA). In contrast to NECA, the partial agonist was still able to bind to the receptor construct locked in its inactive conformation by an S913.39K mutation, although with an affinity lower than that at the native receptor. This could explain the compound's partial agonistic activity: while full A2AAR agonists bind exclusively to the active conformation, likely following conformational selection, partial agonists bind to active as well as inactive conformations, showing higher affinity for the active conformation. This might be a general mechanism of partial agonism also applicable to other GPCRs.

GPCRs in the round: SMA-like copolymers and SMALPs as a platform for investigating GPCRs

G-protein-coupled receptors (GPCRs) are the largest family of membrane proteins, regulate a plethora of physiological responses and are the therapeutic target for 30-40% of clinically-prescribed drugs. They are integral membrane proteins deeply embedded in the plasma membrane where they activate intracellular signalling via coupling to G-proteins and β-arrestin. GPCRs are in intimate association with the bilayer lipids and that lipid environment regulates the signalling functions of GPCRs. This complex lipid 'landscape' is both heterogeneous and dynamic. GPCR function is modulated by bulk membrane properties including membrane fluidity, microdomains, curvature, thickness and asymmetry but GPCRs are also regulated by specific lipid:GPCR binding, including cholesterol and anionic lipids. Understanding the molecular mechanisms whereby GPCR signalling is regulated by lipids is a very active area of research currently. A major advance in membrane protein research in recent years was the application of poly(styrene-co-maleic acid) (SMA) copolymers. These spontaneously generate SMA lipid particles (SMALPs) encapsulating membrane protein in a nano-scale disc of cell membrane, thereby removing the historical need for detergent and preserving lipid:GPCR interaction. The focus of this review is how GPCR-SMALPs are increasing our understanding of GPCR structure and function at the molecular level. Furthermore, an increasing number of 'second generation' SMA-like copolymers have been reported recently. These are reviewed from the context of increasing our understanding of GPCR molecular mechanisms. Moreover, their potential as a novel platform for downstream biophysical and structural analyses is assessed and looking ahead, the translational application of SMA-like copolymers to GPCR drug discovery programmes in the future is considered.

Determinants of Neutral Antagonism and Inverse Agonism in the β2-Adrenergic Receptor

Free-energy profiles for the activation/deactivation of the β2-adrenergic receptor (ADRB2) with neutral antagonist and inverse agonist ligands have been determined with well-tempered multiple-walker (MW) metadynamics simulations. The inverse agonists carazolol and ICI118551 clearly favor single inactive conformational minima in both the binary and ternary ligand-receptor-G-protein complexes, in accord with the inverse-agonist activity of the ligands. The behavior of neutral antagonists is more complex, as they seem also to affect the recruitment of the G-protein. The results are analyzed in terms of the conformational states of the well-known microswitches that have been proposed as indicators of receptor activity.

"Photo-Adrenalines": Photoswitchable β2 -Adrenergic Receptor Agonists as Molecular Probes for the Study of Spatiotemporal Adrenergic Signaling

β2 -adrenergic receptor (β2 -AR) agonists are used for the treatment of asthma and chronic obstructive pulmonary disease, but also play a role in other complex disorders including cancer, diabetes and heart diseases. As the cellular and molecular mechanisms in various cells and tissues of the β2 -AR remain vastly elusive, we developed tools for this investigation with high temporal and spatial resolution. Several photoswitchable β2 -AR agonists with nanomolar activity were synthesized. The most potent agonist for β2 -AR with reasonable switching is a one-digit nanomolar active, trans-on arylazopyrazole-based adrenaline derivative and comprises valuable photopharmacological properties for further biological studies with high structural accordance to the native ligand adrenaline.

Binding kinetics drive G protein subtype selectivity at the β1-adrenergic receptor

G protein-coupled receptors (GPCRs) bind to different G protein α-subtypes with varying degrees of selectivity. The mechanism by which GPCRs achieve this selectivity is still unclear. Using 13C methyl methionine and 19F NMR, we investigate the agonist-bound active state of β1AR and its ternary complexes with different G proteins in solution. We find the receptor in the ternary complexes adopts very similar conformations. In contrast, the full agonist-bound receptor active state assumes a conformation differing from previously characterised activation intermediates or from β1AR in ternary complexes. Assessing the kinetics of binding for the agonist-bound receptor with different G proteins, we find the increased affinity of β1AR for Gs results from its much faster association with the receptor. Consequently, we suggest a kinetic-driven selectivity gate between canonical and secondary coupling which arises from differential favourability of G protein binding to the agonist-bound receptor active state.

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.

IUPHAR themed review: Opioid efficacy, bias, and selectivity

Drugs acting at the opioid receptor family are clinically used to treat chronic and acute pain, though they represent the second line of treatment behind GABA analogs, antidepressants and SSRI's. Within the opioid family mu and kappa opioid receptor are commonly targeted. However, activation of the mu opioid receptor has side effects of constipation, tolerance, dependence, euphoria, and respiratory depression; activation of the kappa opioid receptor leads to dysphoria and sedation. The side effects of mu opioid receptor activation have led to mu receptor drugs being widely abused with great overdose risk. For these reasons, newer safer opioid analgesics are in high demand. For many years a focus within the opioid field was finding drugs that activated the G protein pathway at mu opioid receptor, without activating the β-arrestin pathway, known as biased agonism. Recent advances have shown that this may not be the way forward to develop safer analgesics at mu opioid receptor, though there is still some promise at the kappa opioid receptor. Here we discuss recent novel approaches to develop safer opioid drugs including efficacy vs bias and fine-tuning receptor activation by targeting sub-pockets in the orthosteric site, we explore recent works on the structural basis of bias, and we put forward the suggestion that Gα subtype selectivity may be an exciting new area of interest.

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.

Polyamine detergents tailored for native mass spectrometry studies of membrane proteins

Native mass spectrometry (MS) is a powerful technique for interrogating membrane protein complexes and their interactions with other molecules. A key aspect of the technique is the ability to preserve native-like structures and noncovalent interactions, which can be challenging depending on the choice of detergent. Different strategies have been employed to reduce charge on protein complexes to minimize activation and preserve non-covalent interactions. Here, we report the synthesis of a class of polyamine detergents tailored for native MS studies of membrane proteins. These detergents, a series of spermine covalently attached to various alkyl tails, are exceptional charge-reducing molecules, exhibiting a ten-fold enhanced potency over spermine. Addition of polyamine detergents to proteins solubilized in maltoside detergents results in improved, charge-reduced native mass spectra and reduced dissociation of subunits. Polyamine detergents open new opportunities to investigate membrane proteins in different detergent environments that have thwarted previous native MS studies.

A Novel "Activation Switch" Motif Common to All Aminergic Receptors

Aminergic receptors are G protein-coupled receptors (GPCRs) that transduce signals from small endogenous biogenic amines to regulate intracellular signaling pathways. Agonist binding in the ligand binding pocket on the extracellular side opens and prepares a cavity on the intracellular face of the receptors to interact with and activate G proteins and β-arrestins. Here, by reviewing and analyzing all available aminergic receptor structures, we seek to identify activation-related conformational changes that are independent of the specific scaffold of the bound agonist, which we define as "activation conformational changes" (ACCs). While some common intracellular ACCs have been well-documented, identifying common extracellular ACCs, including those in the ligand binding pocket, is complicated by local adjustments to different ligand scaffolds. Our analysis shows no common ACCs at the extracellular ends of the transmembrane helices. Furthermore, the restricted access to the ligand binding pocket identified previously in some receptors is not universal. Notably, the Trp6.48 toggle switch and the Pro5.50-Ile3.40-Phe6.44 (PIF) motif at the bottom of the ligand binding pocket have previously been proposed to mediate the conformational consequences of ligand binding to the intracellular side of the receptors. Our analysis shows that common ACCs in the ligand binding pocket are associated with the PIF motif and nearby residues, including Trp6.48, but fails to support a shared rotamer toggle associated with activation. However, we identify two common rearrangements between the extracellular and middle subsegments, and propose a novel "activation switch" motif common to all aminergic receptors. This motif includes the middle subsegments of transmembrane helices 3, 5, and 6 and integrates both the PIF motif and Trp6.48.

Genetic Variants of Adrenoceptors

Adrenoceptors are class A G-protein-coupled receptors grouped into three families (α1-, α2-, and β-adrenoceptors), each one including three members. All nine corresponding adrenoceptor genes display genetic variation in their coding and adjacent non-coding genomic region. Coding variants, i.e., nucleotide exchanges within the transcribed and translated receptor sequence, may result in a difference in amino acid sequence thus altering receptor function and signaling. Such variants have been intensely studied in vitro in overexpression systems and addressed in candidate-gene studies for distinct clinical parameters. In recent years, large cohorts were analyzed in genome-wide association studies (GWAS), where variants are detected as significant in context with specific traits. These studies identified two of the in-depth characterized 18 coding variants in adrenoceptors as repeatedly statistically significant genetic risk factors - p.Arg389Gly in the β1- and p.Thr164Ile in the β2-adrenoceptor, along with 56 variants in the non-coding regions adjacent to the adrenoceptor gene loci, the functional role of which is largely unknown at present. This chapter summarizes current knowledge on the two coding variants in adrenoceptors that have been consistently validated in GWAS and provides a prospective overview on the numerous non-coding variants more recently attributed to adrenoceptor gene loci.

Molecular recognition of an aversive odorant by the murine trace amine-associated receptor TAAR7f

There are two main families of G protein-coupled receptors that detect odours in humans, the odorant receptors (ORs) and the trace amine-associated receptors (TAARs). Their amino acid sequences are distinct, with the TAARs being most similar to the aminergic receptors such as those activated by adrenaline, serotonin and histamine. To elucidate the structural determinants of ligand recognition by TAARs, we have determined the cryo-EM structure of a murine receptor, mTAAR7f, coupled to the heterotrimeric G protein Gs and bound to the odorant N,N-dimethylcyclohexylamine (DMCH) to an overall resolution of 2.9 Å. DMCH is bound in a hydrophobic orthosteric binding site primarily through van der Waals interactions and a strong charge-charge interaction between the tertiary amine of the ligand and an aspartic acid residue. This site is distinct and non-overlapping with the binding site for the odorant propionate in the odorant receptor OR51E2. The structure, in combination with mutagenesis data and molecular dynamics simulations suggests that the activation of the receptor follows a similar pathway to that of the β-adrenoceptors, with the significant difference that DMCH interacts directly with one of the main activation microswitch residues.

The dynamics of agonist-β2-adrenergic receptor activation induced by binding of GDP-bound Gs protein

There is considerable uncertainty about the mechanism by which the β₂-adrenergic receptor (β₂AR) is activated. Here we use molecular metadynamics computations to predict the mechanism by which an agonist induces the activation of the β₂AR and its cognate Gs protein. We found that binding agonist alone to the inactive β₂AR does not break the ionic lock and hence does not drive the β₂AR towards the activated conformation. However, we found that attaching the inactive Gs protein to the agonist-bound inactive β₂AR (containing the ionic lock) leads to partial insertion of Gαs-α5 into the core of β₂AR, which breaks the ionic lock, leading to activation of the Gs protein coupled to β₂AR. Upon activation, the Gαs protein undergoes a remarkable opening of the GDP binding pocket, making the GDP available for exchange or release. Concomitantly, Gαs-α5 undergoes a remarkable expansion in the β₂AR cytoplasmic region after the ionic lock is broken, inducing TM6 to displace outward by ~5 Å from TM3.

Dietary Protective Potential of Fucoxanthin as an Active Food Component on Neurological Disorders

The prevalence of neurodegenerative, cerebrovascular, and psychiatric diseases and other neurological disorders has increased dramatically worldwide. Fucoxanthin is an algal pigment with many biological functions, and there is rising evidence that fucoxanthin plays a preventive and therapeutic role in neurological disorders. This review focuses on the metabolism, bioavailability, and blood-brain barrier penetration of fucoxanthin. Furthermore, the neuroprotective potential of fucoxanthin in neurodegenerative diseases, cerebrovascular diseases, and psychiatric diseases as well as other neurological disorders such as epilepsy, neuropathic pain, and brain tumors by acting on multiple targets will be summarized. The multiple targets include regulating apoptosis, reducing oxidative stress, activating the autophagy pathway, inhibiting Aβ aggregation, improving dopamine secretion, reducing α-synuclein aggregation, attenuating neuroinflammation, modulating gut microbiota, and activating brain-derived neurotrophic factor, etc. Additionally, we look forward to brain-targeted oral transport systems due to the low bioavailability and blood-brain barrier permeability of fucoxanthin. We also propose exploring the systemic mechanisms of fucoxanthin metabolism and transport through the gut-brain process and envision new therapeutic targets for fucoxanthin to act on the central nervous system. Finally, we propose dietary fucoxanthin delivery interventions to achieve preventive effects on neurological disorders. This review provides a reference for the application of fucoxanthin in the neural field.

Quantitative analysis of sterol-modulated monomer-dimer equilibrium of the β1-adrenergic receptor by DEER spectroscopy

G protein-coupled receptors (GPCR) activate numerous intracellular signaling pathways. The oligomerization properties of GPCRs, and hence their cellular functions, may be modulated by various components within the cell membrane (such as the presence of cholesterol). Modulation may occur directly via specific interaction with the GPCR or indirectly by affecting the physical properties of the membrane. Here, we use pulsed Q-band double electron-electron resonance (DEER) spectroscopy to probe distances between R1 nitroxide spin labels attached to Cys163 and Cys344 of the β₁-adrenergic receptor (β₁AR) in n-dodecyl-β-D-maltoside micelles upon titration with two soluble cholesterol analogs, cholesteryl hemisuccinate (CHS) and sodium cholate. The former, like cholesterol, inserts itself into the lipid membrane, parallel to the phospholipid chains; the latter is aligned parallel to the surface of membranes. Global quantitative analysis of DEER echo curves upon titration of spin-labeled β₁AR with CHS and sodium cholate reveal the following: CHS binds specifically to the β₁AR monomer at a site close to the Cys163-R1 spin label with an equilibrium dissociation constant [Formula: see text] ~1.4 ± 0.4 mM. While no direct binding of sodium cholate to the β₁AR receptor was observed by DEER, sodium cholate induces specific β₁AR dimerization ([Formula: see text] ~35 ± 6 mM and a Hill coefficient n ~ 2.5 ± 0.4) with intersubunit contacts between transmembrane helices 1 and 2 and helix 8. Analysis of the DEER data obtained upon the addition of CHS to the β₁AR dimer in the presence of excess cholate results in dimer dissociation with species occupancies as predicted from the individual K_D values.

Comparative Study of Receptor-, Receptor State-, and Membrane-Dependent Cholesterol Binding Sites in A2A and A1 Adenosine Receptors Using Coarse-Grained Molecular Dynamics Simulations

We used coarse-grained molecular dynamics (CG MD) simulations to study protein-cholesterol interactions for different activation states of the A_2A adenosine receptor (A_2AR) and the A₁ adenosine receptor (A₁R) and predict new cholesterol binding sites indicating amino acid residues with a high residence time in three biologically relevant membranes. Compared to 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC)-cholesterol and POPC-phosphatidylinositol-bisphosphate (PIP₂)-cholesterol, the plasma mimetic membrane best described the cholesterol binding sites previously detected for the inactive state of A_2AR and revealed the binding sites with long-lasting amino acid residues. We observed that using the plasma mimetic membrane and plotting residues with cholesterol residence time ≥2 μs, our CG MD simulations captured most obviously the cholesterol-protein interactions. For the inactive A_2AR, we identified one more binding site in which cholesterol is bound to residues with a long residence time compared to the previously detected, for the active A₁R, three binding sites, and for the inactive A₁R, two binding sites. We calculated that for the active states, cholesterol binds to residues with a much longer residence time compared to the inactive state for both A_2AR and A₁R. The stability of the identified binding sites to A₁R or A_2AR with CG MD simulations was additionally investigated with potential of mean force calculations using umbrella sampling. We observed that the binding sites with residues to which cholesterol has a long residence time in A_2AR have shallow binding free energy minima compared to the related binding sites in A₁R, suggesting a stronger binding for cholesterol to A₁R. The differences in binding sites in which cholesterol is stabilized and interacts with residues with a long residence time between active and inactive states of A₁R and A_2AR can be important for differences in functional activity and orthosteric agonist or antagonist affinity and can be used for the design of allosteric modulators, which can bind through lipid pathways. We observed a stronger binding for cholesterol to A₁R (i.e., generally higher association rates) compared to A_2AR, which remains to be demonstrated. For the active states, cholesterol binds to residues with much longer residence times compared to the inactive state for both A_2AR and A₁R. Taken together, binding sites of active A₁R may be considered as promising allosteric targets.

Rational design, synthesis, and pharmacological evaluation of a cohort of novel beta-adrenergic receptors ligands enables an assessment of structure-activity relationships

Biomedical applications of molecules that are able to modulate β-adrenergic signaling have become increasingly attractive over the last decade, revealing that β-adrenergic receptors (β-ARs) are key targets for a plethora of therapeutic interventions, including cancer. Despite successes in β-AR drug discovery, identification of β-AR ligands that are useful as selective chemical tools in pharmacological studies of the three β-AR subtypes, or lead compounds for drug development is still a highly challenging task. This is mainly due to the intrinsic plasticity of β-ARs as G protein-coupled receptors in conjunction with the requirement for functional receptor subtype selectivity, tissue specificity and minimal off-target effects. With the aim to provide insight into structure-activity relationships for the three β-AR subtypes, we have synthesized and obtained the pharmacological profile of a series of structurally diverse compounds (named MC) that were designed based on the aryloxy-propanolamine scaffold of SR59230A. Comparative analysis of their predicted binding mode within the active and inactive states of the receptors in combination with their pharmacological profile revealed key structural elements that control their activity as agonists or antagonists, in addition to clues about substituents that mediate selectivity for one receptor subtype over the others. We anticipate that these results will facilitate selective β-AR drug development efforts.

Enhanced peroxidase-like activity of 2(3), 9(10), 16(17), 23(24)-octamethoxyphthalocyanine modified CoFe LDH for a sensor array for reducing substances with catechol structure

Improving the catalytic activity of artificial nanozymes to realize the real-time detection of small molecules becomes an important task. Herein, a highly active nanozyme, 2(3), 9(10), 16(17), 23(24)-octamethoxyphthalocyanine (Pc(OH)₈) modified CoFe LDH microspheres (Pc(OH)₈-CoFe LDH) have been prepared by the two-step hydrothermal method. The 3,3',5,5'-tetramylbenzidine (TMB), a chromogenic substrate, was fast oxidized into blue oxTMB by H₂O₂ in the presence of Pc(OH)₈-CoFe LDH, indicating that Pc(OH)₈-CoFe LDH possesses high peroxidase-like activity rather than pure CoFe LDH. The enhancement peroxidase-like activity of the Pc(OH)₈-CoFe LDH is ascribed to the synergistic action between Pc(OH)₈ and CoFe LDH. Experimental results of radical scavenger and fluorescence probe verify that superoxide radical (•O₂^-) plays an important role during the catalytic reaction. Interestingly, the absorption intensity of reaction system has been enhanced largely, due to adding of the reducing substances containing catechol structure. Based on this, the three reducing substances (dopamine, procyanidin B2, catechins) containing catechol structure were distinguished from other reducing substances without catechol structure. Thus, a colorimetric array has been constructed using reaction time as the sensing element to realize the sensitive and selective recognition of catechol structures at a certain concentration.

Protein Design Strategies for the Structural–Functional Studies of G Protein-Coupled Receptors

G protein-coupled receptors (GPCRs) are an important family of membrane proteins responsible for many physiological functions in human body. High resolution GPCR structures are required to understand their molecular mechanisms and perform rational drug design, as GPCRs play a crucial role in a variety of diseases. That is difficult to obtain for the wild-type proteins because of their low stability. In this review, we discuss how this problem can be solved by using protein design strategies developed to obtain homogeneous stabilized GPCR samples for crystallization and cryoelectron microscopy.

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.

Structure of the GOLD-domain seven-transmembrane helix protein family member TMEM87A

TMEM87s are eukaryotic transmembrane proteins with two members (TMEM87A and TMEM87B) in humans. TMEM87s have proposed roles in protein transport to and from the Golgi, as mechanosensitive ion channels, and in developmental signaling. TMEM87 disruption has been implicated in cancers and developmental disorders. To better understand TMEM87 structure and function, we determined a cryo-EM structure of human TMEM87A in lipid nanodiscs. TMEM87A consists of a Golgi-dynamics (GOLD) domain atop a membrane-spanning seven-transmembrane helix domain with a large cavity open to solution and the membrane outer leaflet. Structural and functional analyses suggest TMEM87A may not function as an ion channel or G-protein coupled receptor. We find TMEM87A shares its characteristic domain arrangement with seven other proteins in humans; three that had been identified as evolutionary related (TMEM87B, GPR107, and GPR108) and four previously unrecognized homologs (GPR180, TMEM145, TMEM181, and WLS). Among these structurally related GOLD domain seven-transmembrane helix (GOST) proteins, WLS is best characterized as a membrane trafficking and secretion chaperone for lipidated Wnt signaling proteins. We find key structural determinants for WLS function are conserved in TMEM87A. We propose TMEM87A and structurally homologous GOST proteins could serve a common role in trafficking membrane-associated cargo.

Best Practices in Constant pH MD Simulations: Accuracy and Sampling

Various approaches have been proposed to include the effect of pH in molecular dynamics (MD) simulations. Among these, the λ-dynamics approach proposed by Brooks and co-workers [Kong, X.; Brooks III, C. L. J. Chem. Phys.1996, 105, 2414−2423] can be performed with little computational overhead and hfor each typeence be used to routinely perform MD simulations at microsecond time scales, as shown in the accompanying paper [Aho, N. et al. J. Chem. Theory Comput.2022, DOI: 10.1021/acs.jctc.2c00516]. At such time scales, however, the accuracy of the molecular mechanics force field and the parametrization becomes critical. Here, we address these issues and provide the community with guidelines on how to set up and perform long time scale constant pH MD simulations. We found that barriers associated with the torsions of side chains in the CHARMM36m force field are too high for reaching convergence in constant pH MD simulations on microsecond time scales. To avoid the high computational cost of extending the sampling, we propose small modifications to the force field to selectively reduce the torsional barriers. We demonstrate that with such modifications we obtain converged distributions of both protonation and torsional degrees of freedom and hence consistent pK_a estimates, while the sampling of the overall configurational space accessible to proteins is unaffected as compared to normal MD simulations. We also show that the results of constant pH MD depend on the accuracy of the correction potentials. While these potentials are typically obtained by fitting a low-order polynomial to calculated free energy profiles, we find that higher order fits are essential to provide accurate and consistent results. By resolving problems in accuracy and sampling, the work described in this and the accompanying paper paves the way to the widespread application of constant pH MD beyond pK_a prediction.

Finding the Perfect Fit: Conformational Biosensors to Determine the Efficacy of GPCR Ligands

G protein-coupled receptors (GPCRs) are highly druggable targets that adopt numerous conformations. A ligand's ability to stabilize specific conformation(s) of its cognate receptor determines its efficacy or ability to produce a biological response. Identifying ligands that produce different receptor conformations and potentially discrete pharmacological effects (e.g., biased agonists, partial agonists, antagonists, allosteric modulators) is a major goal in drug discovery and necessary to develop drugs with better effectiveness and fewer side effects. Fortunately, direct measurements of ligand efficacy, via receptor conformational changes are possible with the recent development of conformational biosensors. In this review, we discuss classical efficacy models, including the two-state model, the ternary-complex model, and multistate models. We describe how nanobody-, transducer-, and receptor-based conformational biosensors detect and/or stabilize specific GPCR conformations to identify ligands with different levels of efficacy. In particular, conformational biosensors provide the potential to identify and/or characterize therapeutically desirable but often difficult to measure conformations of receptors faster and better than current methods. For drug discovery/development, several recent proof-of-principle studies have optimized conformational biosensors for high-throughput screening (HTS) platforms. However, their widespread use is limited by the fact that few sensors are reliably capable of detecting low-frequency conformations and technically demanding assay conditions. Nonetheless, conformational biosensors do help identify desirable ligands such as allosteric modulators, biased ligands, or partial agonists in a single assay, representing a distinct advantage over classical methods.

19F NMR: A promising tool for dynamic conformational studies of G protein-coupled receptors

Advances in X-ray crystallography and cryoelectron microscopy enabled unprecedented insights into the activation processes of G protein-coupled receptors (GPCRs). However, these static receptor structures provide limited information about dynamics and conformational transitions that play pivotal roles in mediating signaling diversity through the multifaceted interactions between ligands, receptors, and transducers. Developing NMR approaches to probe the dynamics of conformational transitions will push the frontier of receptor science toward a more comprehensive understanding of these signaling processes. Although much progress has been made during the last decades, it remains challenging to delineate receptor conformational states and interrogate the functions of the individual states at a quantitative level. Here we cover the progress of ¹⁹F NMR applications in GPCR conformational and dynamic studies during the past 20 years. Current challenges and limitations of ¹⁹F NMR for studying GPCR dynamics are also discussed, along with experimental strategies that will drive this field forward.

Filling of a water-free void explains the allosteric regulation of the β1-adrenergic receptor by cholesterol

Recent high-pressure NMR results indicate that the preactive conformation of the β₁-adrenergic receptor (β₁AR) harbours completely empty cavities of ~100 ų volume, which disappear in the active conformation of the receptor. Here we have localized these cavities using X-ray crystallography of xenon-derivatized β₁AR crystals. One of the cavities is in direct contact with the cholesterol-binding pocket. Solution NMR shows that addition of the cholesterol analogue cholesteryl hemisuccinate impedes the formation of the active conformation of detergent-solubilized β₁AR by blocking conserved G protein-coupled receptor microswitches, concomitant with an affinity reduction of both isoprenaline and G protein-mimicking nanobody Nb80 for β₁AR detected by isothermal titration calorimetry. This wedge-like action explains the function of cholesterol as a negative allosteric modulator of β₁AR. A detailed understanding of G protein-coupled receptor regulation by cholesterol by filling of a dry void and the easy scouting for such voids by xenon may provide new routes for the development of allosteric drugs.

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.

A cholesterol analog stabilizes the human β2-adrenergic receptor nonlinearly with temperature

In cell membranes, G protein-coupled receptors (GPCRs) interact with cholesterol, which modulates their assembly, stability, and conformation. Previous studies have shown how cholesterol modulates the structural properties of GPCRs at ambient temperature. Here, we characterized the mechanical, kinetic, and energetic properties of the human β2-adrenergic receptor (β2AR) in the presence and absence of the cholesterol analog cholesteryl hemisuccinate (CHS) at room temperature (25°C), at physiological temperature (37°C), and at high temperature (42°C). We found that CHS stabilized various structural regions of β2AR differentially, which changed nonlinearly with temperature. Thereby, the strongest effects were observed for structural regions that are important for receptor signaling. Moreover, at 37°C, but not at 25° or 42°C, CHS caused β2AR to increase and stabilize conformational substates to adopt to basal activity. These findings indicate that the nonlinear, temperature-dependent action of CHS in modulating the structural and functional properties of this GPCR is optimized for 37°C.

Beta-2 Adrenergic Receptor Gene Expression in HER2-Positive Early-Stage Breast Cancer Patients: A Post-hoc Analysis of the NCCTG-N9831 (Alliance) Trial

BACKGROUND

Beta-2 adrenergic receptor (ß2AR) modulates immune activation and may enhance trastuzumab activity. We assessed the impact of ß2AR gene (ADRB2) expression on the outcomes of patients with HER2-positive early-stage breast cancer enrolled on the NCCTG-N9831 trial.

PATIENTS AND METHODS

This is a post-hoc analysis of the NCCTG-N9831 trial, which compared chemotherapy (arm A) versus chemotherapy plus trastuzumab (arms B&C) as adjuvant treatment of patients with HER2-positive early-stage breast cancer, with disease-free survival (DFS) as primary endpoint. Gene expression levels retrieved by DASL assay were used to classify patients as ADRB2-high or ADRB2-low. Hazard ratios (HRs) were calculated by a Cox proportional model adjusted for prognostic variables and ADRB2 expression. Correlations between ADRB2 expression and stromal tumor-infiltrating lymphocyte (TIL) levels were assessed with Pearson coefficient. A multivariable Cox regression model with interaction term was performed to assess the interaction between ADRB2 expression and treatment arm; and ADRB2 expression and a 8-gene signature previously shown to predict trastuzumab benefit.

RESULTS

Overall, 1,282 patients were included (ADRB2-high [N = 944] / ADRB2-low [N = 338]). A high expression of ADRB2 was associated with a longer DFS (P = .01) in the overall population. The addition of trastuzumab to chemotherapy improved DFS only in patients with ADRB2-high tumors (P < .01). ADRB2 expression was correlated with TIL levels (r = 0.24, P < .001). No association between ADRB2 expression and the 8-gene trastuzumab benefit signature was observed (P = .32).

CONCLUSION

Our findings suggest that a high ADRB2 expression is a favorable prognostic factor and may identify patients with HER2-positive early-stage breast cancer who benefit from adjuvant trastuzumab.

TRIAL REGISTRATION

clinicaltrials.gov NCT00005970.

Cryo-EM structures of the β3 adrenergic receptor bound to solabegron and isoproterenol

The β₃-adrenergic receptor (β₃AR) is the most essential drug target for overactive bladder and has therapeutic potentials for the treatments of type 2 diabetes and obesity. Here, we report the cryo-electron microscopy structures of the β₃AR-G_s signaling complexes with the selective agonist, solabegron and the nonselective agonist, isoproterenol. Comparison of the isoproterenol-, mirabegron-, and solabegron-bound β₃AR structures revealed that the extracellular loop 2 changes its conformation depending on the bound agonist and plays an essential role in solabegron binding. Moreover, β₃AR has an intrinsically narrow exosite, regardless of the agonist type. This structural feature clearly explains why β₃AR prefers mirabegron and solabegron, as the narrow exosite is suitable for binding with agonists with elongated shapes. Our study deepens the understanding of the binding characteristics of β₃AR agonists and may pave the way for developing β₃AR-selective drugs.

Structural insights into ligand recognition, activation, and signaling of the α2A adrenergic receptor

The α_2A adrenergic receptor (α_2AAR) is a G protein (heterotrimeric guanine nucleotide-binding protein)-coupled receptor that mediates important physiological functions in response to the endogenous neurotransmitters norepinephrine and epinephrine, as well as numerous chemically distinct drugs. However, the molecular mechanisms of drug actions remain poorly understood. Here, we report the cryo-electron microscopy structures of the human α_2AAR-GoA complex bound to norepinephrine and three imidazoline derivatives (brimonidine, dexmedetomidine, and oxymetazoline). Together with mutagenesis and functional data, these structures provide important insights into the molecular basis of ligand recognition, activation, and signaling at the α_2AAR. Further structural analyses uncover different molecular determinants between α_2AAR and βARs for recognition of norepinephrine and key regions that determine the G protein coupling selectivity. Overall, our studies provide a framework for understanding the signal transduction of the adrenergic system at the atomic level, which will facilitate rational structure-based discovery of safer and more effective medications for α_2AAR.

Simplified quality assessment for small-molecule ligands in the Protein Data Bank

More than 70% of the experimentally determined macromolecular structures in the Protein Data Bank (PDB) contain small-molecule ligands. Quality indicators of ∼643,000 ligands present in ∼106,000 PDB X-ray crystal structures have been analyzed. Ligand quality varies greatly with regard to goodness of fit between ligand structure and experimental data, deviations in bond lengths and angles from known chemical structures, and inappropriate interatomic clashes between the ligand and its surroundings. Based on principal component analysis, correlated quality indicators of ligand structure have been aggregated into two largely orthogonal composite indicators measuring goodness of fit to experimental data and deviation from ideal chemical structure. Ranking of the composite quality indicators across the PDB archive enabled construction of uniformly distributed composite ranking score. This score is implemented at RCSB.org to compare chemically identical ligands in distinct PDB structures with easy-to-interpret two-dimensional ligand quality plots, allowing PDB users to quickly assess ligand structure quality and select the best exemplars.

The Anti-Inflammatory Effect of the β1-Adrenergic Receptor Antagonist Metoprolol on High Glucose Treated Human Microvascular Retinal Endothelial Cells

Hyperglycemia-induced impairment of the blood-retinal barrier represents the main pathological event in diabetic retinopathy that is elicited by a reduced cellular response to an accumulation of reactive oxygen species (ROS) and increased inflammation. The purpose of the study was to evaluate whether the selective β1-adrenoreceptor (β1-AR) antagonist metoprolol could modulate the inflammatory response to hyperglycemic conditions. For this purpose, human retinal endothelial cells (HREC) were treated with normal (5 mM) or high glucose (25 mM, HG) in the presence of metoprolol (10 μM), epinephrine (1 μM), or both compounds. Metoprolol prevented both the HG-induced reduction of cell viability (MTT assays) and the modulation of the angiogenic potential of HREC (tube formation assays) reducing the TNF-α, IL-1β, and VEGF mRNA levels (qRT-PCR). Moreover, metoprolol prevented the increase in phospho-ERK1/2, phospho-cPLA₂, COX2, and protein levels (Western blot) as well as counteracting the translocation of ERK1/2 and cPLA₂ (high-content screening). Metoprolol reduced ROS accumulation in HG-stimulated HREC by activating the anti-oxidative cellular response mediated by the Keap1/Nrf2/HO-1 pathway. In conclusion, metoprolol exerted a dual effect on HG-stimulated HREC, decreasing the activation of the pro-inflammatory ERK1/2/cPLA₂/COX2 axis, and counteracting ROS accumulation by activating the Keap1/Nrf2/HO-1 pathway.

From structure to clinic: Design of a muscarinic M1 receptor agonist with potential to treatment of Alzheimer's disease

Current therapies for Alzheimer's disease seek to correct for defective cholinergic transmission by preventing the breakdown of acetylcholine through inhibition of acetylcholinesterase, these however have limited clinical efficacy. An alternative approach is to directly activate cholinergic receptors responsible for learning and memory. The M1-muscarinic acetylcholine (M1) receptor is the target of choice but has been hampered by adverse effects. Here we aimed to design the drug properties needed for a well-tolerated M1-agonist with the potential to alleviate cognitive loss by taking a stepwise translational approach from atomic structure, cell/tissue-based assays, evaluation in preclinical species, clinical safety testing, and finally establishing activity in memory centers in humans. Through this approach, we rationally designed the optimal properties, including selectivity and partial agonism, into HTL9936-a potential candidate for the treatment of memory loss in Alzheimer's disease. More broadly, this demonstrates a strategy for targeting difficult GPCR targets from structure to clinic.

Structural basis for molecular recognition of G protein-coupled estrogen receptor by selected bisphenols

Complicated ligand-dependent signaling pathways of bisphenol A (BPA) and its analogues involve not only intranuclear estrogen receptor but also membrane receptor G protein-coupled estrogen receptor (GPER). However, the structural basis for molecular recognition of GPER by the environmental chemicals remains unknown. To reveal the structural dependence of GPER recognition by bisphenols, a systematic molecular dynamics simulation study was performed for selected bisphenols with different electron hybrid orbitals and substituents on their C atoms connecting two phenol rings. BPA was used as a control, bisphenol C(BPC) as an example for a connecting C with sp² hybrid orbitals to provide more ligand rigidity, bisphenol E(BPE) and bisphenol F(BPF) for decreased steric hindrance and hydrophobicity around the connecting C, and bisphenol B(BPB) and bisphenol AF(BPAF) for increased hydrophobicity and steric hindrance. All the tested bisphenols can bind with GPER at its classic orthosteric site to obtain GPER-ligand complexes, while van der Waals interactions and direct inter-molecular electrostatic energies provide the driving forces for ligand binding. Bulky substituents and structural rigidity of the connecting C dramatically impair hydrogen bonding between GPER and the bisphenols, which results in decreased contribution of both favorable intermolecular hydrogen bonds and unfavorable polar solvation effect to complex stability of BPB and BPC since decreased number of key residues is expected. Increase in substituent lipophilicity enhances the van der Waals interactions and favorable non-polar solvation effect. The six bisphenols of high structural similarity shared two key recognition residues, Leu137^TM3 and Trp272^TM6, the latter of which was in the highly conserved CWxP motif of TM6 and has been reported as key residue for G protein-coupled receptor activation. Based on the obtained knowledge, GPER affinity and relevant toxicity of BPA alternatives can be easily predicted, and the calculated binding free energies are consistent with the available experimental observations.

Structure determination of GPCRs: cryo-EM compared with X-ray crystallography

G protein-coupled receptors (GPCRs) are the largest single family of cell surface receptors encoded by the human genome and they play pivotal roles in co-ordinating cellular systems throughout the human body, making them ideal drug targets. Structural biology has played a key role in defining how receptors are activated and signal through G proteins and β-arrestins. The application of structure-based drug design (SBDD) is now yielding novel compounds targeting GPCRs. There is thus significant interest from both academia and the pharmaceutical industry in the structural biology of GPCRs as currently only about one quarter of human non-odorant receptors have had their structure determined. Initially, all the structures were determined by X-ray crystallography, but recent advances in electron cryo-microscopy (cryo-EM) now make GPCRs tractable targets for single-particle cryo-EM with comparable resolution to X-ray crystallography. So far this year, 78% of the 99 GPCR structures deposited in the PDB (Jan-Jul 2021) were determined by cryo-EM. Cryo-EM has also opened up new possibilities in GPCR structural biology, such as determining structures of GPCRs embedded in a lipid nanodisc and multiple GPCR conformations from a single preparation. However, X-ray crystallography still has a number of advantages, particularly in the speed of determining many structures of the same receptor bound to different ligands, an essential prerequisite for effective SBDD. We will discuss the relative merits of cryo-EM and X-ray crystallography for the structure determination of GPCRs and the future potential of both techniques.

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 receptors (GPCRs) are responsible for most cytoplasmic signaling in response to extracellular ligands with different efficacy profiles. Various spectroscopic techniques have identified that agonists exhibiting varying efficacies can selectively stabilize a specific conformation of the receptor. However, the structural basis for activation of the GPCR-G protein complex by ligands with different efficacies is incompletely understood. To better understand the structural basis underlying the mechanisms by which ligands with varying efficacies differentially regulate the conformations of receptors and G proteins, we determined the structures of β2AR-Gαs[Formula: see text]γ bound with partial agonist salbutamol or bound with full agonist isoprenaline using single-particle cryo-electron microscopy at resolutions of 3.26 Å and 3.80 Å, respectively. Structural comparisons between the β2AR-Gs-salbutamol and β2AR-Gs-isoprenaline complexes demonstrated that the decreased binding affinity and efficacy of salbutamol compared with those of isoprenaline might be attributed to weakened hydrogen bonding interactions, attenuated hydrophobic interactions in the orthosteric binding pocket and different conformational changes in the rotamer toggle switch in TM6. Moreover, the observed stronger interactions between the intracellular loop 2 or 3 (ICL2 or ICL3) of β2AR and Gαs with binding of salbutamol versus isoprenaline might decrease phosphorylation in the salbutamol-activated β2AR-Gs complex. From the observed structural differences between these complexes of β2AR, a mechanism of β2AR activation by partial and full agonists is proposed to provide structural insights into β2AR desensitization.

Membrane-facilitated receptor access and binding mechanisms of long-acting β2-adrenergic receptor (β2-AR) agonists

The drugs salmeterol, formoterol, and salbutamol constitute the frontline treatment for asthma and other chronic pulmonary diseases. These drugs activate the β2-adrenergic receptors (β2-AR), a class A G-protein-coupled receptor (GPCR) and differ significantly in their clinical onset and duration of actions. According to the "microkinetic model," the long duration of action of salmeterol and formoterol compared to salbutamol were attributed, at least in part, to their high lipophilicity and increased local concentrations in the membrane near the receptor. However, the structural and molecular bases of how the lipophilic drugs reach the binding site of the receptor from the surrounding membrane remain unknown. Using a variety of classical and enhanced molecular dynamics simulation techniques, we investigated the membrane partitioning characteristics, binding, and unbinding mechanisms of the ligands. The obtained results offer remarkable insight into the functional role of membrane lipids in the ligand association process. Strikingly, salmeterol entered the binding site from the bilayer through transmembrane helices 1 and 7. The entry was preceded by membrane-facilitated rearrangement and presentation of its phenyl-alkoxy-alkyl tail as a passkey to an access route gated by F193, a residue known critical for salmeterol's affinity. Formoterol's access is through the aqueous path shared by other β2-AR agents. We observed a novel secondary path for salbutamol that is distinct from its primary route. Our study offers a mechanistic description for the membrane-facilitated access and binding of ligands to β2-AR and establishes a groundwork for recognizing membrane lipids as an integral component in the molecular recognition process. Significance Statement The cell membrane's functional role behind the duration of action of long-acting β2-adrenergic receptor (β2-AR) agonists such as salmeterol has been a subject of debate for a long time. We investigated the binding and unbinding mechanisms of the three commonly used β2-AR agonists, salmeterol, formoterol, and salbutamol, using advanced simulation techniques. The obtained results offer unprecedented insights into the active role of membrane lipids in facilitating access and binding of the ligands, affecting the molecular recognition process and their pharmacology.

Cryo-EM structure of the β3-adrenergic receptor reveals the molecular basis of subtype selectivity

The β₃-adrenergic receptor (β₃AR) is predominantly expressed in adipose tissue and urinary bladder and has emerged as an attractive drug target for the treatment of type 2 diabetes, obesity, and overactive bladder (OAB). Here, we report the cryogenic electron microscopy structure of the β₃AR-G_s signaling complex with the selective agonist mirabegron, a first-in-class drug for OAB. Comparison of this structure with the previously reported β₁AR and β₂AR structures reveals a receptor activation mechanism upon mirabegron binding to the orthosteric site. Notably, the narrower exosite in β₃AR creates a perpendicular pocket for mirabegron. Mutational analyses suggest that a combination of both the exosite shape and the amino-acid-residue substitutions defines the drug selectivity of the βAR agonists. Our findings provide a molecular basis for βAR subtype selectivity, allowing the design of more-selective agents with fewer adverse effects.

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.

Optimized damping parameters for empirical dispersion corrections to symmetry-adapted perturbation theory

Symmetry-adapted perturbation theory (SAPT) has become an invaluable tool for studying the fundamental nature of non-covalent interactions by directly computing the electrostatics, exchange (steric) repulsion, induction (polarization), and London dispersion contributions to the interaction energy using quantum mechanics. Further application of SAPT is primarily limited by its computational expense, where even its most affordable variant (SAPT0) scales as the fifth power of system size [O(N⁵)] due to the dispersion terms. The algorithmic scaling of SAPT0 is reduced from O(N⁵)→O(N⁴) by replacing these terms with the empirical D3 dispersion correction of Grimme and co-workers, forming a method that may be termed SAPT0-D3. Here, we optimize the damping parameters for the -D3 terms in SAPT0-D3 using a much larger training set than has previously been considered, namely, 8299 interaction energies computed at the complete-basis-set limit of coupled cluster through perturbative triples [CCSD(T)/CBS]. Perhaps surprisingly, with only three fitted parameters, SAPT0-D3 improves on the accuracy of SAPT0, reducing mean absolute errors from 0.61 to 0.49 kcal mol^-1 over the full set of complexes. Additionally, SAPT0-D3 exhibits a nearly 2.5× speedup over conventional SAPT0 for systems with ∼300 atoms and is applied here to systems with up to 459 atoms. Finally, we have also implemented a functional group partitioning of the approach (F-SAPT0-D3) and applied it to determine important contacts in the binding of salbutamol to G-protein coupled β₁-adrenergic receptor in both active and inactive forms. SAPT0-D3 capabilities have been added to the open-source Psi4 software.

Crystal structure of dopamine D1 receptor in complex with G protein and a non-catechol agonist

Dopamine D1 receptor (D1R) is an important drug target implicated in many psychiatric and neurological disorders. Selective agonism of D1R are sought to be the therapeutic strategy for these disorders. Most selective D1R agonists share a dopamine-like catechol moiety in their molecular structure, and their therapeutic potential is therefore limited by poor pharmacological properties in vivo. Recently, a class of non-catechol D1R selective agonists with a distinct scaffold and pharmacological properties were reported. Here, we report the crystal structure of D1R in complex with stimulatory G protein (Gs) and a non-catechol agonist Compound 1 at 3.8 Å resolution. The structure reveals the ligand bound to D1R in an extended conformation, spanning from the orthosteric site to extracellular loop 2 (ECL2). Structural analysis reveals that the unique features of D1R ligand binding pocket explains the remarkable selectivity of this scaffold for D1R over other aminergic receptors, and sheds light on the mechanism for D1R activation by the non-catechol agonist.