Involvement of Asn-293 in stereospecific agonist recognition and in activation of the beta 2-adrenergic receptor

Can the heptapeptide ASSIVSF of the β2-adrenoceptor recognize ephedrine and pseudoephedrine epimers in a complex system?

Epimer separation is crucial in the field of analytical chemistry, separation science, and the pharmaceutical industry. No reported methods could separate simultaneously epimers or even isomers and remove other unwanted, co-existing, interfering substances from complex systems like herbal extracts. Herein, we prepared a heptapeptide-modified stationary phase for the separation of 1R,2S-(-)-ephedrine [(-)-Ephe] and 1S,2S-(+)-pseudoephedrine [(+)-Pse] epimers from Ephedra sinica Stapf extract and blood samples. The heptapeptide stationary phase was comprehensively characterized by scanning electron microscopy, X-ray photoelectron spectroscopy, and Fourier transform infrared spectroscopy. The separation efficiency of the heptapeptide column was compared with an affinity column packed with full-length β2-AR functionalized silica gel (β2-AR column). The binding affinity of the heptapeptide with (+)-Pse was 3-fold greater than that with (-)-Ephe. Their binding mechanisms were extensively characterized by chromatographic analysis, ultraviolet spectra, circular dichroism analysis, isothermal titration calorimetry, and molecule docking. An enhanced hydrogen bonding was clearly observed in the heptapeptide-(+)-Pse complex. Such results demonstrated that the heptapeptide can recognize (+)-Pse and (-)-Ephe epimers in a complex system. This work, we believe, was the first report to simultaneously separate epimers and remove non-specific interfering substances from complex samples. The method was potentially applicable to more challenging sample separation, such as chiral separation from complex systems.

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

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

Stereoselectivity in the Membrane Transport of Phenylethylamine Derivatives by Human Monoamine Transporters and Organic Cation Transporters 1, 2, and 3

Stereoselectivity is well known and very pronounced in drug metabolism and receptor binding. However, much less is known about stereoselectivity in drug membrane transport. Here, we characterized the stereoselective cell uptake of chiral phenylethylamine derivatives by human monoamine transporters (NET, DAT, and SERT) and organic cation transporters (OCT1, OCT2, and OCT3). Stereoselectivity differed extensively between closely related transporters. High-affinity monoamine transporters (MATs) showed up to 2.4-fold stereoselective uptake of norepinephrine and epinephrine as well as of numerous analogs. While NET and DAT preferentially transported (S)-norepinephrine, SERT preferred the (R)-enantiomer. In contrast, NET and DAT showed higher transport for (R)-epinephrine and SERT for (S)-epinephrine. Generally, MAT stereoselectivity was lower than expected from their high affinity to several catecholamines and from the high stereoselectivity of some inhibitors used as antidepressants. Additionally, the OCTs differed strongly in their stereoselectivity. While OCT1 showed almost no stereoselective uptake, OCT2 was characterized by a roughly 2-fold preference for most (R)-enantiomers of the phenylethylamines. In contrast, OCT3 transported norphenylephrine and phenylephrine with 3.9-fold and 3.3-fold preference for their (R)-enantiomers, respectively, while the para-hydroxylated octopamine and synephrine showed no stereoselective OCT3 transport. Altogether, our data demonstrate that stereoselectivity is highly transporter-to-substrate specific and highly diverse even between homologous transporters.

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.

How Do Modulators Affect the Orthosteric and Allosteric Binding Pockets?

Allosteric modulators (AMs) that bind allosteric sites can exhibit greater selectivity than the orthosteric ligands and can either enhance agonist-induced receptor activity (termed positive allosteric modulator or PAM), inhibit agonist-induced activity (negative AM or NAM), or have no effect on activity (silent AM or SAM). Until now, it is not clear what the exact effects of AMs are on the orthosteric active site or the allosteric binding pocket(s). In the present work, we collected both the three-dimensional (3D) structures of receptor-orthosteric ligand and receptor-orthosteric ligand-AM complexes of a specific target protein. Using our novel algorithm toolset, molecular complex characterizing system (MCCS), we were able to quantify the key residues in both the orthosteric and allosteric binding sites along with potential changes of the binding pockets. After analyzing 21 pairs of 3D crystal or cryo-electron microscopy (cryo-EM) complexes, including 4 pairs of GPCRs, 5 pairs of ion channels, 11 pairs of enzymes, and 1 pair of transcription factors, we found that the binding of AMs had little impact on both the orthosteric and allosteric binding pockets. In return, given the accurately predicted allosteric binding pocket(s) of a drug target of medicinal interest, we can confidently conduct the virtual screening or lead optimization without concern that the huge conformational change of the pocket could lead to the low accuracy of virtual screening.

Unraveling the hidden temporal range of fast β2-adrenergic receptor mobility by time-resolved fluorescence

G-protein-coupled receptors (GPCRs) are hypothesized to possess molecular mobility over a wide temporal range. Until now the temporal range has not been fully accessible due to the crucially limited temporal range of available methods. This in turn, may lead relevant dynamic constants to remain masked. Here, we expand this dynamic range by combining fluorescent techniques using a spot confocal setup. We decipher mobility constants of β₂-adrenergic receptor over a wide time range (nanosecond to second). Particularly, a translational mobility (10 µm²/s), one order of magnitude faster than membrane associated lateral mobility that explains membrane protein turnover and suggests a wider picture of the GPCR availability on the plasma membrane. And a so far elusive rotational mobility (1-200 µs) which depicts a previously overlooked dynamic component that, despite all complexity, behaves largely as predicted by the Saffman-Delbrück model.

The mobility of the β2-adrenergic receptor, from the nanosecond to the second range, is revealed by combining several fluorescent spectroscopy techniques. These data also show a previously hidden mobility of this GPCR.

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.

Unleashing the shape of L-DOPA at last

Herein, we report the first rotational study of neutral L-DOPA, an extensively used supramolecular synthon and an amino acid precursor of the neurotransmitters dopamine, norepinephrine (noradrenaline), and epinephrine (adrenaline) using broadband and narrowband Fourier transform microwave spectroscopies coupled with a laser ablation vaporization system. The spectroscopic parameters derived from the analysis of the rotational spectrum conclusively identify the existence of four distinct conformers of L-DOPA in the supersonic jet, further rejecting the previously reported catechol ring-induced conformational restriction. The analysis of the ¹⁴N nuclear quadrupole coupling hyperfine structure further revealed the orientation of the N-bearing functional group, proving the existence of stabilizing N-H⋯π interactions for the observed structures.

Identification and characterization of an atypical Gαs-biased β2AR agonist that fails to evoke airway smooth muscle cell tachyphylaxis

We sought β₂AR agonists for treating obstructive lung diseases such as asthma, in which this receptor relaxes airway smooth muscle (ASM) cells and opens airways. Agonists favoring Gs coupling (leads to airway relaxation) compared with activating β-arrestin (limits effectiveness due to receptor desensitization) were pursued in a 40-million-compound screening library. Of several agonists identified, one was apparently biased away from β-arrestin. Agonist–receptor–G protein modeling revealed different receptor interactions compared with other agonists. The favorable effects of the apparent biasing with this agonist were demonstrated in a physiologic system (ASM relaxation). These studies point to a different structural class of β-agonists that might be used to treat obstructive lung diseases without the adverse effects associated with tachyphylaxis.

G protein–coupled receptors display multifunctional signaling, offering the potential for agonist structures to promote conformational selectivity for biased outputs. For β₂-adrenergic receptors (β₂AR), unbiased agonists stabilize conformation(s) that evoke coupling to Gαs (cyclic adenosine monophosphate [cAMP] production/human airway smooth muscle [HASM] cell relaxation) and β-arrestin engagement, the latter acting to quench Gαs signaling, contributing to receptor desensitization/tachyphylaxis. We screened a 40-million-compound scaffold ranking library, revealing unanticipated agonists with dihydroimidazolyl-butyl-cyclic urea scaffolds. The S-stereoisomer of compound C1 shows no detectable β-arrestin engagement/signaling by four methods. However, C1-S retained Gαs signaling—a divergence of the outputs favorable for treating asthma. Functional studies with two models confirmed the biasing: β₂AR-mediated cAMP signaling underwent desensitization to the unbiased agonist albuterol but not to C1-S, and desensitization of HASM cell relaxation was observed with albuterol but not with C1-S. These HASM results indicate biologically pertinent biasing of C1-S, in the context of the relevant physiologic response, in the human cell type of interest. Thus, C1-S was apparently strongly biased away from β-arrestin, in contrast to albuterol and C5-S. C1-S structural modeling and simulations revealed binding differences compared with unbiased epinephrine at transmembrane (TM) segments 3,5,6,7 and ECL2. C1-S (R2 = cyclohexane) was repositioned in the pocket such that it lost a TM6 interaction and gained a TM7 interaction compared with the analogous unbiased C5-S (R2 = benzene group), which appears to contribute to C1-S biasing away from β-arrestin. Thus, an agnostic large chemical-space library identified agonists with receptor interactions that resulted in relevant signal splitting of β₂AR actions favorable for treating obstructive lung disease.

An Intracellular Tripeptide Arg-His-Trp of Serum Origin Detected in MCF-7 Cells Is A Possible Agonist to β2 Adrenoceptor

BACKGROUND

The need of agonists and antagonists of β2 adrenoceptor (β2AR) is warranted in various human disease conditions including cancer, cardiovascular and other metabolic disorders. However, the sources of agonists of β2AR are diverse in nature. Interestingly, there is a complete gap in the exploration of agonists of β2AR from serum that is a well-known component of culture media which supports growth and proliferation of normal and cancer cells in vitro.

METHODS

In this paper, we employed a novel vertical tube gel electrophoresis (VTGE)-assisted purification of intracellular metabolites of MCF-7 cells grown in vitro in complete media with fetal bovine serum (FBS). Intracellular metabolites of MCF-7 cells were then analyzed by LC-HRMS. Identified intracellular tripeptides of FBS origin were evaluated for their molecular interactions with various extracellular and intracellular receptors including β2AR (PDB ID: 2RH1) by employing molecular docking and molecular dynamics simulations (MDS). A known agonist of β2AR, isoproterenol was used as a positive control in molecular docking and MDS analyses.

RESULTS

We report here identification of a few novel intracellular tripeptides, namely Arg-His-Trp, (PubChem CID-145453842), Pro-Ile-Glu, (PubChem CID-145457492), Cys-Gln-Gln, (PubChem CID-71471965), Glu-Glu-Lys, (PubChem CID-11441068) and Gly-Cys-Leu (PubChem CID-145455600) of FBS origin in MCF-7 cells. Molecular docking and MDS analyses revealed that among these molecules, the tripeptide Arg-His-Trp shows a favorable binding affinity with β2AR (-9.8 Kcal/mol). The agonistic effect of Arg-His-Trp is significant and comparable with that of a known agonist of β2AR, isoproterenol.

CONCLUSION

In conclusion, we identified a unique Arg-His-Trp tripeptide of FBS origin in MCF-7 cells by employing a novel approach. This unique tripeptide Arg-His-Trp is suggested to be a potential agonist of β2AR and it may have applications in the context of various human diseases like bronchial asthma and chronic obstructive pulmonary disease (COPD).

Ligand Pose Predictions for Human G Protein-Coupled Receptors: Insights from the Amber-Based Hybrid Molecular Mechanics/Coarse-Grained Approach

Human G protein-coupled receptors (hGPCRs) are the most frequent targets of Food and Drug Administration (FDA)-approved drugs. Structural bioinformatics, along with molecular simulation, can support structure-based drug design targeting hGPCRs. In this context, several years ago, we developed a hybrid molecular mechanics (MM)/coarse-grained (CG) approach to predict ligand poses in low-resolution hGPCR models. The approach was based on the GROMOS96 43A1 and PRODRG united-atom force fields for the MM part. Here, we present a new MM/CG implementation using, instead, the Amber 14SB and GAFF all-atom potentials for proteins and ligands, respectively. The new implementation outperforms the previous one, as shown by a variety of applications on models of hGPCR/ligand complexes at different resolutions, and it is also more user-friendly. Thus, it emerges as a useful tool to predict poses in low-resolution models and provides insights into ligand binding similarly to all-atom molecular dynamics, albeit at a lower computational cost.

Can the Partial Peptide SIVSF of the β2-Adrenergic Receptor Recognize Chirality of the Epinephrine Neurotransmitter?

Chirality plays an essential role in biological molecular recognition, such as neurotransmission. Here, we applied electrospray-cold ion trap spectroscopy to complexes of a partial binding motif SIVSF of a β₂-adrenergic receptor pocket with L- and D-epinephrine AdH⁺. The ultraviolet spectrum of the SIVSF-AdH⁺ complex changed drastically when L-AdH⁺ was replaced by its enantiomer. The isomer-selected infrared spectra revealed that D-AdH⁺ was bound to SIVSF by its protonated amino-group or a single catechol OH and induced nonhelical secondary structures of SIVSF. This is in sharp contrast to the helical SIVSF complex with L-AdH⁺, which is close to the natural binding structure with two catechol OHs binding in the receptor. This shows that a short pentapeptide SIVSF can distinguish the chirality of the ligand AdH⁺ as well as the receptor. This stereoselectivity is suggested to arise from an additional interaction involving the hydroxyl group on the chiral carbon.

Aminergic GPCR-Ligand Interactions: A Chemical and Structural Map of Receptor Mutation Data

The aminergic family of G protein-coupled receptors (GPCRs) plays an important role in various diseases and represents a major drug discovery target class. Structure determination of all major aminergic subfamilies has enabled structure-based ligand design for these receptors. Site-directed mutagenesis data provides an invaluable complementary source of information for elucidating the structural determinants of binding of different ligand chemotypes. The current study provides a comparative analysis of 6692 mutation data points on 34 aminergic GPCR subtypes, covering the chemical space of 540 unique ligands from mutagenesis experiments and information from experimentally determined structures of 52 distinct aminergic receptor-ligand complexes. The integrated analysis enables detailed investigation of structural receptor-ligand interactions and assessment of the transferability of combined binding mode and mutation data across ligand chemotypes and receptor subtypes. An overview is provided of the possibilities and limitations of using mutation data to guide the design of novel aminergic receptor ligands.

Structural insights into binding specificity, efficacy and bias of a β2AR partial agonist

Salmeterol is a partial agonist for the β2 adrenergic receptor (β2AR) and the first long-acting β2AR agonist to be widely used clinically for the treatment of asthma and chronic obstructive pulmonary disease. Salmeterol's safety and mechanism of action have both been controversial. To understand its unusual pharmacological action and partial agonism, we obtained the crystal structure of salmeterol-bound β2AR in complex with an active-state-stabilizing nanobody. The structure reveals the location of the salmeterol exosite, where sequence differences between β1AR and β2AR explain the high receptor-subtype selectivity. A structural comparison with the β2AR bound to the full agonist epinephrine reveals differences in the hydrogen-bond network involving residues Ser2045.43 and Asn2936.55. Mutagenesis and biophysical studies suggested that these interactions lead to a distinct active-state conformation that is responsible for the partial efficacy of G-protein activation and the limited β-arrestin recruitment for salmeterol.

Implications for human odor sensing revealed from the statistics of odorant-receptor interactions

Binding of odorants to olfactory receptors (ORs) elicits downstream chemical and neural signals, which are further processed to odor perception in the brain. Recently, Mainland and colleagues have measured more than 500 pairs of odorant-OR interaction by a high-throughput screening assay method, opening a new avenue to understanding the principles of human odor coding. Here, using a recently developed minimal model for OR activation kinetics, we characterize the statistics of OR activation by odorants in terms of three empirical parameters: the half-maximum effective concentration EC₅₀, the efficacy, and the basal activity. While the data size of odorants is still limited, the statistics offer meaningful information on the breadth and optimality of the tuning of human ORs to odorants, and allow us to relate the three parameters with the microscopic rate constants and binding affinities that define the OR activation kinetics. Despite the stochastic nature of the response expected at individual OR-odorant level, we assess that the confluence of signals in a neuron released from the multitude of ORs is effectively free of noise and deterministic with respect to changes in odorant concentration. Thus, setting a threshold to the fraction of activated OR copy number for neural spiking binarizes the electrophysiological signal of olfactory sensory neuron, thereby making an information theoretic approach a viable tool in studying the principles of odor perception.

Despite the decades of research, quantitative details of human olfaction have remained largely unexplored. However, a high-throughput measurement has recently been carried out to produce dose-response data between a set of odorants and a repertoire of human olfactory receptors. We characterized each pair of odorant-receptor interaction in terms of EC50, efficacy, and basal level, a strategy often adopted in biochemical, pharmacological sciences to describe the response of receptors to cognate ligands. The distributions of EC50 values and efficacies acquired from the analysis provide glimpses into how human olfactory receptors are tuned to odorants. Specifically, the response of human ORs is optimized around ∼ 100μM of odorant. Next, the efficacies of OR responses to odorants are bi-exponentially distributed, which indicates that the strength of odorant-OR interaction is classified into strong and weak subgroups. By showing that the stochastic response of individual receptor to odorant can effectively be binarized at cellular level through olfactory processes, we also provide a theoretical basis for an information theoretical approach in studying the principles of odor perception.

Stereoselective binding of agonists to the β2-adrenergic receptor: insights into molecular details and thermodynamics from molecular dynamics simulations

Enhanced-sampling molecular dynamics simulations show the molecular details of stereospecific binding of agonists to the β₂-adrenergic receptor.

Investigation of allosteric coupling in human β2-adrenergic receptor in the presence of intracellular loop 3

Background

This study investigates the allosteric coupling that exists between the intra- and extracellular parts of human β₂-adrenergic receptor (β₂-AR), in the presence of the intracellular loop 3 (ICL3), which is missing in all crystallographic experiments and most of the simulation studies reported so far. Our recent 1 μs long MD run has revealed a transition to the so-called very inactive state of the receptor, in which ICL3 packed under the G protein’s binding cavity and completely blocked its accessibility to G protein. Simultaneously, an outward tilt of transmembrane helix 5 (TM5) caused an expansion of the extracellular ligand-binding site. In the current study, we performed independent runs with a total duration of 4 μs to further investigate the very inactive state with packed ICL3 and the allosteric coupling event (three unrestrained runs and five runs with bond restraints at the ligand-binding site).

Results

In all three independent unrestrained runs (each 500 ns long), ICL3 preserved its initially packed/closed conformation within the studied time frame, suggesting an inhibition of the receptor’s activity. Specific bond restraints were later imposed between some key residues at the ligand-binding site, which have been experimentally determined to interact with the ligand. Restraining the binding site region to an open state facilitated ICL3 closure, whereas a relatively constrained/closed binding site hindered ICL3 packing. However, the reverse operation, i.e. opening of the packed ICL3, could not be realized by restraining the binding site region to a closed state. Thus, any attempt failed to free the ICL3 from its locked state due to the presence of persistent hydrogen bonds.

Conclusions

Overall, our simulations indicated that starting with very inactive states, the receptor stayed almost irreversibly inhibited, which in turn decreased the overall mobility of the receptor. Bond restraints which represented the geometric restrictions caused by ligands of various sizes when bound at the ligand-binding site, induced the expected conformational changes in TM5, TM6 and consequently, ICL3. Still, once ICL3 was packed, the allosteric coupling became ineffective due to strong hydrogen bonds connecting ICL3 to the core of the receptor.

Electronic supplementary material

The online version of this article (doi:10.1186/s12900-016-0061-9) contains supplementary material, which is available to authorized users.

How To Minimize Artifacts in Atomistic Simulations of Membrane Proteins, Whose Crystal Structure Is Heavily Engineered: β₂-Adrenergic Receptor in the Spotlight

Atomistic molecular dynamics (MD) simulations are used extensively to elucidate membrane protein properties. These simulations are based on three-dimensional protein structures that in turn are often based on crystallography. The protein structures resolved in crystallographic studies typically do not correspond to pristine proteins, however. Instead the crystallized proteins are commonly engineered, including structural modifications (mutations, replacement of protein sequences by antibodies, bound ligands, etc.) whose impact on protein structure and dynamics is largely unknown. Here we explore this issue through atomistic MD simulations (∼5 μs in total), focusing on the β2-adrenergic receptor (β2AR) that is one of the most studied members of the G-protein coupled receptor superfamily. Starting from an inactive-state crystal structure of β2AR, we remove the many modifications in β2AR systematically one at a time, in six consecutive steps. After each step, we equilibrate the system and simulate it quite extensively. The results of this step-by-step approach highlight that the structural modifications used in crystallization can affect ligand and G-protein binding sites, packing at the transmembrane-helix interface region, and the dynamics of connecting loops in β2AR. When the results of the systematic step-by-step approach are compared to an all-at-once technique where all modifications done on β2AR are removed instantaneously at the same time, it turns out that the step-by-step method provides results that are superior in terms of maintaining protein structural stability. The results provide compelling evidence that for membrane proteins whose 3D structure is based on structural engineering, the preparation of protein structure for atomistic MD simulations is a delicate and sensitive process. The results show that most valid results are found when the structural modifications are reverted slowly, one at a time.

Fourier transform microwave spectroscopy of Ac-Ser-NH2: the role of side chain interactions in peptide folding

Serine capped dipeptide N-acetyl-l-serinamide (Ac-Ser-NH2) has been investigated using Fourier transform microwave spectroscopic techniques combined with laser ablation sources. Spectral signatures originating from one dominant species have been detected in the supersonic expansion. Rotational and nuclear quadrupole coupling constants of the two (14)N nuclei have been used in the characterization of a C/γ-turn structure, which is stabilized by a CO∙∙∙HN intramolecular hydrogen bond closing a seven-membered ring. Two extra hydrogen bonds involving the polar side chain (-CH2OH) further stabilize the structure. The non-observation of C5 species, attributed to the presence of the polar side chain, is in contrast with the previous gas phase observation of the related dipeptides containing glycine or alanine residues. The A-E splitting pattern arising from the internal rotation of the methyl group has been analyzed and the internal rotation barrier has been determined.

Structural signatures of DRD4 mutants revealed using molecular dynamics simulations: Implications for drug targeting

Human Dopamine Receptor D4 (DRD4) orchestrates several neurological functions and represents a target for many psychological disorders. Here, we examined two rare variants in DRD4; V194G and R237L, which elicit functional alterations leading to disruption of ligand binding and G protein coupling, respectively. Using atomistic molecular dynamics (MD) simulations, we provide in-depth analysis to reveal structural signatures of wild and mutant complexes with their bound agonist and antagonist ligands. We constructed intra-protein network graphs to discriminate the global conformational changes induced by mutations. The simulations also allowed us to elucidate the local side-chain dynamical variations in ligand-bound mutant receptors. The data suggest that the mutation in transmembrane V (V194G) drastically disrupts the organization of ligand binding site and causes disorder in the native helical arrangement. Interestingly, the R237L mutation leads to significant rewiring of side-chain contacts in the intracellular loop 3 (site of mutation) and also affects the distant transmembrane topology. Additionally, these mutations lead to compact ICL3 region compared to the wild type, indicating that the receptor would be inaccessible for G protein coupling. Our findings thus reveal unreported structural determinants of the mutated DRD4 receptor and provide a robust framework for design of effective novel drugs.

Multiresidue Method for Analysis of β Agonists in Swine Urine by Enzyme Linked Receptor Assay Based on β2 Adrenergic Receptor Expressed in HEK293 Cells

A novel enzyme-linked receptor assay (ELRA) based on β₂-adrenergic receptor (β₂-AR) has been developed for rapid and high-throughput detection of β-adrenergic agonists (β-agonists) in urine. Human embryonic kidney cells (HEK293) were introduced as the expression system to enhance the functionality of the recombinant β₂-AR, and the attempt to detect β-agonists in swine urine using such approaches was accomplished unprecedentedly. In this article, a recombinant porcine β₂-AR was produced in the inner membrane of HEK293 cells and purified from crude membrane protein by nickel-nitrilotriacetic acid affinity chromatography. After activity identification, the recombinant receptor was used in the development of direct competitive ELRA. Several parameters such as blocking buffer and blocking process were optimized and the performance of the system was determined. The IC50 concentrations of clenbuterol, salbutamol, and ractopamine were 34, 53 and 63 μg/L, and the average recovery rates were 68.2%, 60.3% and 65.5%, respectively. ELRA based on β₂-AR shows a series of advantages such as safety, easy operation, and high efficiency, making it promising for the rapid screening of β-agonists in animal urine.

Principal component analysis of binding energies for single-point mutants of hT2R16 bound to an agonist correlate with experimental mutant cell response

Directed evolution is a technique that enables the identification of mutants of a particular protein that carry a desired property by successive rounds of random mutagenesis, screening, and selection. This technique has many applications, including the development of G protein-coupled receptor-based biosensors and designer drugs for personalized medicine. Although effective, directed evolution is not without challenges and can greatly benefit from the development of computational techniques to predict the functional outcome of single-point amino acid substitutions. In this article, we describe a molecular dynamics-based approach to predict the effects of single amino acid substitutions on agonist binding (salicin) to a human bitter taste receptor (hT2R16). An experimentally determined functional map of single-point amino acid substitutions was used to validate the whole-protein molecular dynamics-based predictive functions. Molecular docking was used to construct a wild-type agonist-receptor complex, providing a starting structure for single-point substitution simulations. The effects of each single amino acid substitution in the functional response of the receptor to its agonist were estimated using three binding energy schemes with increasing inclusion of solvation effects. We show that molecular docking combined with molecular mechanics simulations of single-point mutants of the agonist-receptor complex accurately predicts the functional outcome of single amino acid substitutions in a human bitter taste receptor.

Agonist binding by the β2-adrenergic receptor: an effect of receptor conformation on ligand association-dissociation characteristics

The β₂-adrenergic receptor (β₂-AR), a G protein-coupled receptor (GPCR), is a physiologically important transmembrane protein that is a target for drugs used for treatment of asthma and cardiovascular diseases. Study of the first steps of ligand recognition and the molecular basis of ligand binding to the orthosteric site is essential for understanding the pharmacological properties of the receptor. In this work we investigated the characteristic features of the agonist association–dissociation process to and from the different conformational forms of β₂-AR by use of advanced molecular modeling techniques. The investigation was focused on estimating the free energy profiles (FEPs) corresponding to the process of a full agonist ((R,R)-fenoterol) and an inverse agonist (carazolol) binding and unbinding to and from β₂-AR. The two different conformational forms of β₂-AR, i.e. active β₂-AR–PDB: 3P0G and inactive β₂-AR–PDB: 2RH1 were included in this stage of the study. We revealed several significant qualitative differences between FEPs characteristic of both conformational forms. Both FEPs suggest the existence of three transient binding sites in the extracellular domain of β₂-AR. Comparison of the residues surrounding these transient binding sites in both β₂-AR states revealed the importance of the aromatic residues F194, H93^2.64, H296^6.58, and H178 (extracellular part of β₂-AR) in the early stages of the binding process. In addition, slightly different exit and entry paths are preferred by the ligand molecule in the extracellular part of β₂-AR, depending on the conformation of the receptor.

Understanding the effects on constitutive activation and drug binding of a D130N mutation in the β2 adrenergic receptor via molecular dynamics simulation

G-protein-coupled receptors (GPCRs) are currently one of the largest families of drug targets. The constitutive activation induced by mutation of key GPCR residues is associated closely with various diseases. However, the structural basis underlying such activation and its role in drug binding has remained unclear. Herein, we used all-atom molecular dynamics simulations and free energy calculations to study the effects of a D130N mutation on the structure of β2 adrenergic receptor (β2AR) and its binding of the agonist salbutamol. The results indicate that the mutation caused significant changes in some key helices. In particular, the mutation leads to the departure of transmembrane 3 (TM3) from transmembrane 6 (TM6) and marked changes in the NPxxY region as well as the complete disruption of a key ionic lock, all of which contribute to the observed constitutive activation. In addition, the D130N mutation weakens some important H-bonds, leading to structural changes in these regions. Binding free energy calculations indicate that van der Waals and electrostatic interactions are the main driving forces in binding salbutamol; however, binding strength in the mutant β2AR is significantly enhanced mainly through modifying electrostatic interactions. Further analysis revealed that the increase in binding energy upon mutation stems mainly from the H-bonds formed between the hydroxyl group of salbutamol and the serine residues of TM5. This observation suggests that modifications of the H-bond groups of this drug could significantly influence drug efficacy in the treatment of diseases associated with this mutation.

A structural chemogenomics analysis of aminergic GPCRs: lessons for histamine receptor ligand design

BACKGROUND AND PURPOSE

Chemogenomics focuses on the discovery of new connections between chemical and biological space leading to the discovery of new protein targets and biologically active molecules. G-protein coupled receptors (GPCRs) are a particularly interesting protein family for chemogenomics studies because there is an overwhelming amount of ligand binding affinity data available. The increasing number of aminergic GPCR crystal structures now for the first time allows the integration of chemogenomics studies with high-resolution structural analyses of GPCR-ligand complexes.

EXPERIMENTAL APPROACH

In this study, we have combined ligand affinity data, receptor mutagenesis studies, and amino acid sequence analyses to high-resolution structural analyses of (hist)aminergic GPCR-ligand interactions. This integrated structural chemogenomics analysis is used to more accurately describe the molecular and structural determinants of ligand affinity and selectivity in different key binding regions of the crystallized aminergic GPCRs, and histamine receptors in particular.

KEY RESULTS

Our investigations highlight interesting correlations and differences between ligand similarity and ligand binding site similarity of different aminergic receptors. Apparent discrepancies can be explained by combining detailed analysis of crystallized or predicted protein-ligand binding modes, receptor mutation studies, and ligand structure-selectivity relationships that identify local differences in essential pharmacophore features in the ligand binding sites of different receptors.

CONCLUSIONS AND IMPLICATIONS

We have performed structural chemogenomics studies that identify links between (hist)aminergic receptor ligands and their binding sites and binding modes. This knowledge can be used to identify structure-selectivity relationships that increase our understanding of ligand binding to (hist)aminergic receptors and hence can be used in future GPCR ligand discovery and design.

Exploring the determinants of trace amine-associated receptor 1's functional selectivity for the stereoisomers of amphetamine and methamphetamine

Amphetamines are widely abused drugs that interfere with dopamine transport and storage. Recently, however, another mechanism of action was identified: stereoselective activation of the GαS protein-coupled trace amine-associated receptor 1 (TAAR1). To identify structural determinants of this stereoselectivity, we functionally evaluated six mutant receptors in vitro and then used homology modeling and dynamic simulation to predict drug affinities. Converting Asp102 to Ala rendered mouse and rat TAAR1 (mTAAR1 and rTAAR1, respectively) insensitive to β-phenylethylamine, amphetamine (AMPH), and methamphetamine (METH). Mutating Met268 in rTAAR1 to Thr shifted the concentration-response profiles for AMPH and METH isomers rightward an order of magnitude, whereas replacing Thr268 with Met in mTAAR1 resulted in profiles leftward shifted 10-30-fold. Replacing Asn287 with Tyr in rTAAR1 produced a mouselike receptor, while the reciprocal mTAAR1 mutant was rTAAR1-like. These results confirm TAAR1 is an AMPH/METH receptor in vitro and establish residues 102 (3.32) and 268 (6.55) as major contributors to AMPH/METH binding with residue 287 (7.39) determining species stereoselectivity.

Kinetics and mechanism of G protein-coupled receptor activation

The activation of a G protein-coupled receptor is generally triggered by binding of an agonist to the receptor's binding pocket, or, in the case of rhodopsin, by light-induced changes of the pre-bound retinal. This is followed by a series of a conformational changes towards an active receptor conformation, which is capable of signalling to G proteins and other downstream proteins. In the past few years, a number of new techniques have been employed to analyze the kinetics of this activation process, including X-ray crystallographic three-dimensional structures of receptors in the inactive and the active states, NMR studies of labelled receptors, molecular simulations, and optical analyses with fluorescence resonance energy transfer (FRET). Here we review our current understanding of the activation process of GPCRs as well as open questions in the sequence of events ranging from (sub-)microsecond activation by light or agonist binding to millisecond activation of receptors by soluble ligands and the subsequent generation of an intracellular signal.

Molecular interactions between fenoterol stereoisomers and derivatives and the β₂-adrenergic receptor binding site studied by docking and molecular dynamics simulations

The β2 adrenergic receptor (β2-AR) has become a model system for studying the ligand recognition process and mechanism of the G protein coupled receptors activation. In the present study stereoisomers of fenoterol and some of its derivatives (N = 94 molecules) were used as molecular probes to identify differences in stereo-recognition interactions between β2-AR and structurally similar agonists. The present study aimed at determining the 3D molecular models of the fenoterol derivative-β2-AR complexes. Molecular models of β2-AR have been developed by using the crystal structure of the human β2-AR T4 lysozyme fusion protein with bound (S)-carazolol (PDB ID: 2RH1) and more recently reported structure of a nanobody-stabilized active state of the β2-AR with the bound full agonist BI-167107 (PDB ID: 3P0G). The docking procedure allowed us to study the similarities and differences in the recognition binding site(s) for tested ligands. The agonist molecules occupied the same binding region, between TM III, TM V, TM VI and TM VII. The residues identified by us during docking procedure (Ser203, Ser207, Asp113, Lys305, Asn312, Tyr308, Asp192) were experimentally indicated in functional and biophysical studies as being very important for the agonist-receptor interactions. Moreover, the additional space, an extension of the orthosteric pocket, was identified and described. Furthermore, the molecular dynamics simulations were used to study the molecular mechanism of interaction between ligands ((R,R')- and (S,S')-fenoterol) and β2-AR. Our research offers new insights into the ligand stereoselective interaction with one of the most important GPCR member. This study may also facilitate the design of improved selective medications, which can be used to treat, prevent and control heart failure symptoms.

Determinants involved in subtype-specific functions of rat trace amine-associated receptors 1 and 4

AIMS

The trace amine-associated receptor (Taar) family displays high species- and subtype-specific pharmacology. Several trace amines such as β-phenylethylamine (β-PEA), p-tyramine and tryptamine are agonists at TA(1) but poorly activate rat and mouse Taar4.

PRINCIPAL RESULTS

Using rat TA(1) and Taar4 chimera, we identified determinants in transmembrane helices 3 and 6, which, when replaced by the corresponding portion of rat TA(1) , can rescue cell surface expression of rat Taar4. When expressed at the cell surface, rat Taar4 pharmacology was very similar to that of TA(1) and coupled to the Gα(s) -protein/AC pathway. Our data suggest that binding pockets of Taar for surrogate agonists overlap between paralogs.

CONCLUSIONS

This implicates that the repertoire of Taar ensures functional redundancy, tissue- and cell-specific expression and/or different downstream signalling rather than different agonist specificity.

Structure of β-adrenergic receptors

β-Adrenergic receptors (βARs) control key physiological functions by transducing signals encoded in catecholamine hormones and neurotransmitters to activate intracellular signaling pathways. As members of the large family of G protein-coupled receptors (GPCRs), βARs have a seven-transmembrane helix topology and signal via G protein- and arrestin-dependent pathways. Until 2007, three-dimensional structural information of GPCRs activated by diffusible ligands, including βARs, was limited to homology models that used the related photoreceptor rhodopsin as a template. Over many years, several labs have developed strategies that have finally allowed the structures of the turkey β(1)AR and the human β(2)AR to be determined experimentally. The challenges to overcome included heterologous receptor overexpression, design of stabilized and crystallizable modified receptor constructs, ligand-affinity purification of active receptor and the development of novel techniques in crystallization and microcrystallography. The structures of βARs in complex with inverse agonists, antagonists, and agonists have revealed the binding mode of ligands with different efficacies, have allowed to obtain insights into ligand selectivity, and have provided better templates for drug design. Also, the structures of β(2)AR in complex with a G protein and a G protein-mimicking nanobody have provided important insights into the mechanism of receptor activation and G protein coupling. This chapter summarizes the strategies and methods that have been successfully applied to the structural studies of βARs. These are exemplified with detailed protocols toward the structure determination of stabilized turkey β(1)AR-ligand complexes. We also discuss the spectacular insights into adrenergic receptor function that were obtained from the structures.

Identifying ligand binding conformations of the β2-adrenergic receptor by using its agonists as computational probes

Recently available G-protein coupled receptor (GPCR) structures and biophysical studies suggest that the difference between the effects of various agonists and antagonists cannot be explained by single structures alone, but rather that the conformational ensembles of the proteins need to be considered. Here we use an elastic network model-guided molecular dynamics simulation protocol to generate an ensemble of conformers of a prototypical GPCR, β(2)-adrenergic receptor (β(2)AR). The resulting conformers are clustered into groups based on the conformations of the ligand binding site, and distinct conformers from each group are assessed for their binding to known agonists of β(2)AR. We show that the select ligands bind preferentially to different predicted conformers of β(2)AR, and identify a role of β(2)AR extracellular region as an allosteric binding site for larger drugs such as salmeterol. Thus, drugs and ligands can be used as "computational probes" to systematically identify protein conformers with likely biological significance.

The essential role for aromatic cluster in the β3 adrenergic receptor

AIM

To explore the function of the conserved aromatic cluster F213(5.47), F308(6.51), and F309(6.52) in human β3 adrenergic receptor (hβ3AR).

METHODS

Point mutation technology was used to produce plasmid mutations of hβ3AR. HEK-293 cells were transiently co-transfected with the hβ3AR (wild-type or mutant) plasmids and luciferase reporter vector pCRE-luc. The expression levels of hβ3AR in the cells were determined by Western blot analysis. The constitutive signalling and the signalling induced by the β3AR selective agonist, BRL (BRL37344), were then evaluated. To further explore the interaction mechanism between BRL and β3AR, a three-dimensional complex model of β3AR and BRL was constructed by homology modelling and molecular docking.

RESULTS

For F308(6.51), Ala and Leu substitution significantly decreased the constitutive activities of β3AR to approximately 10% of that for the wild-type receptor. However, both the potency and maximal efficacy were unchanged by Ala substitution. In the F308(6.51)L construct, the EC(50) value manifested as a "right shift" of approximately two orders of magnitude with an increased E(max). Impressively, the molecular pharmacological phenotype was similar to the wild-type receptor for the introduction of Tyr at position 308(6.51), though the EC(50) value increased by approximately five-fold for the mutant. For F309(6.52), the constitutive signalling for both F309(6.52)A and F309(6.52)L constructs were strongly impaired. In the F309(6.52)A construct, BRL-stimulated signalling showed a normal E(max) but reduced potency. Leu substitution of F309(6.52) reduced both the E(max) and potency. When F309(6.52) was mutated to Tyr, the constitutive activity was decreased approximately three-fold, and BRL-stimulated signalling was significantly impaired. Furthermore, the double mutant (F308(6.51)A_F309(6.52)A) caused the total loss of β3AR function. The predicted binding mode between β3AR and BRL revealed that both F308(6.51) and F309(6.52) were in the BRL binding pocket of β3AR, while F213(5.47) and W305(6.48) were distant from the binding site.

CONCLUSION

These results revealed that aromatic residues, especially F308(6.51) and F309(6.52), play essential roles in the function of β3AR. Aromatic residues maintained the receptor in a partially activated state and significantly contributed to ligand binding. The results supported the common hypothesis that the aromatic cluster F[Y]5.47/F[Y]6.52/F[Y]6.51 conserved in class A G protein-coupled receptor (GPCR) plays an important role in the structural stability and activation of GPCRs.

A prospective cross-screening study on G-protein-coupled receptors: lessons learned in virtual compound library design

We present the systematic prospective evaluation of a protein-based and a ligand-based virtual screening platform against a set of three G-protein-coupled receptors (GPCRs): the β-2 adrenoreceptor (ADRB2), the adenosine A(2A) receptor (AA2AR), and the sphingosine 1-phosphate receptor (S1PR1). Novel bioactive compounds were identified using a consensus scoring procedure combining ligand-based (frequent substructure ranking) and structure-based (Snooker) tools, and all 900 selected compounds were screened against all three receptors. A striking number of ligands showed affinity/activity for GPCRs other than the intended target, which could be partly attributed to the fuzziness and overlap of protein-based pharmacophore models. Surprisingly, the phosphodiesterase 5 (PDE5) inhibitor sildenafil was found to possess submicromolar affinity for AA2AR. Overall, this is one of the first published prospective chemogenomics studies that demonstrate the identification of novel cross-pharmacology between unrelated protein targets. The lessons learned from this study can be used to guide future virtual ligand design efforts.

Receptor conformations involved in dopamine D(2L) receptor functional selectivity induced by selected transmembrane-5 serine mutations

Although functional selectivity is now widely accepted, the molecular basis is poorly understood. We have studied how aspects of transmembrane region 5 (TM5) of the dopamine D(2L) receptor interacts with three rationally selected rigid ligands (dihydrexidine, dinapsoline, and dinoxyline) and the reference compounds dopamine and quinpirole. As was expected from homology modeling, mutation of three TM5 serine residues to alanine (S5.42A, S5.43A, S5.46A) had little effect on antagonist affinity. All three mutations decreased the affinity of the agonist ligands to different degrees, S5.46A being somewhat less affected. Four functions [adenylate cyclase (AC), extracellular signal-regulated kinase 1/2 phosphorylation (MAPK), arachidonic acid release (AA), and guanosine 5'-O-(3-thio)triphosphate binding (GTPγS)] were assessed. The intrinsic activity (IA) of quinpirole was unaffected by any of the mutations, whereas S5.42A and S5.46A mutations abolished the activity of dopamine and the three rigid ligands, although dihydrexidine retained IA at MAPK function only with S5.42A. Remarkably, S5.43A did not markedly affect IA for AC and MAPK for any of the ligands and eliminated AA activity for dinapsoline and dihydrexidine but not dinoxyline. These data suggest that this mutation did not disrupt the overall conformation or signaling ability of the mutant receptors but differentially affected ligand activation. Computational studies indicate that these D(2) agonists stabilize multiple receptor conformations. This has led to models showing the stabilized conformations and interhelical and receptor-ligand contacts corresponding to the different activation pathways stabilized by various agonists. These data provide a basis for understanding D(2L) functional selectivity and rationally discovering functionally selective D(2) drugs.

Fluorescence/bioluminescence resonance energy transfer techniques to study G-protein-coupled receptor activation and signaling

Fluorescence and bioluminescence resonance energy transfer (FRET and BRET) techniques allow the sensitive monitoring of distances between two labels at the nanometer scale. Depending on the placement of the labels, this permits the analysis of conformational changes within a single protein (for example of a receptor) or the monitoring of protein-protein interactions (for example, between receptors and G-protein subunits). Over the past decade, numerous such techniques have been developed to monitor the activation and signaling of G-protein-coupled receptors (GPCRs) in both the purified, reconstituted state and in intact cells. These techniques span the entire spectrum from ligand binding to the receptors down to intracellular second messengers. They allow the determination and the visualization of signaling processes with high temporal and spatial resolution. With these techniques, it has been demonstrated that GPCR signals may show spatial and temporal patterning. In particular, evidence has been provided for spatial compartmentalization of GPCRs and their signals in intact cells and for distinct physiological consequences of such spatial patterning. We review here the FRET and BRET technologies that have been developed for G-protein-coupled receptors and their signaling proteins (G-proteins, effectors) and the concepts that result from such experiments.

Investigation of D₁ receptor-agonist interactions and D₁/D₂ agonist selectivity using a combination of pharmacophore and receptor homology modeling

The aim of this study was to use a combined structure and pharmacophore modeling approach to extract information regarding dopamine D₁ receptor agonism and D₁/D₂ agonist selectivity. A 3D structure model of the D₁ receptor in its agonist-bound state was constructed with a full D₁ agonist present in the binding site. Two different binding modes were identified using (+)-doxanthrine or SKF89626 in the modeling procedure. The 3D model was further compared with a selective D₁ agonist pharmacophore model. The pharmacophore feature arrangement was found to be in good agreement with the binding site composition of the receptor model, but the excluded volumes did not fully reflect the shape of the agonist binding pocket. A new receptor-based pharmacophore model was developed with forbidden volumes centered on atom positions of amino acids in the binding site. The new pharmacophore model showed a similar ability to discriminate as the previous model. A comparison of the 3D structures and pharmacophore models of D₁ and D₂ receptors revealed differences in shape and ligand-interacting features that determine selectivity of D₁ and D₂ receptor agonists. A hydrogen bond pharmacophoric feature (Ser-TM5) was shown to contribute most to the selectivity. Non-conserved residues in the binding pocket that strongly contribute to D₁/D₂ receptor agonist selectivity were also identified; those were Ser/Cys³·³⁶, Tyr/Phe⁵·³⁸, Ser/Tyr⁵·⁴¹, and Asn/His⁶·⁵⁵ in the transmembrane (TM) helix region, together with Ser/Ile and Leu/Asn in the second extracellular loop (EC2). This work provides useful information for the design of new selective D₁ and D₂ agonists. The combined receptor structure and pharmacophore modeling approach is considered to be general, and could therefore be applied to other ligand-protein interactions for which experimental information is limited.

Investigation of D₂ receptor-agonist interactions using a combination of pharmacophore and receptor homology modeling

A combined modeling approach was used to identify structural factors that underlie the structure–activity relationships (SARs) of full dopamine D₂ receptor agonists and structurally similar inactive compounds. A 3D structural model of the dopamine D₂ receptor was constructed, with the agonist (−)-(R)-2-OH-NPA present in the binding site during the modeling procedure. The 3D model was evaluated and compared with our previously published D₂ agonist pharmacophore model. The comparison revealed an inconsistency between the projected hydrogen bonding feature (Ser-TM5) in the pharmacophore model and the TM5 region in the structure model. A new refined pharmacophore model was developed, guided by the shape of the binding site in the receptor model and with less emphasis on TM5 interactions. The combination of receptor and pharmacophore modeling also identified the importance of His393^6.55 for agonist binding. This convergent 3D pharmacophore and protein structure modeling strategy is considered to be general and can be highly useful in less well-characterized systems to explore ligand–receptor interactions. The strategy has the potential to identify weaknesses in the individual models and thereby provides an opportunity to improve the discriminating predictivity of both pharmacophore searches and structure-based virtual screens.

A novel nonribose agonist, LUF5834, engages residues that are distinct from those of adenosine-like ligands to activate the adenosine A(2a) receptor

The recent publication of both the antagonist- and agonist-bound structures of the adenosine A(2A) receptor have revealed much about how a ligand may bind to a receptor and cause the conformational changes associated with agonist-mediated activation. In particular, the agonist-bound structure revealed key interactions between the ribose group of adenosine-derived agonists and amino acids in the receptor binding pocket that lead to receptor activation. However, agonists without a ribose group also exist, and we wondered whether such compounds occupy the same agonist binding site. Therefore we used a mutagenesis approach in this study to investigate the mode of binding of 2-amino-4-(4-hydroxyphenyl)- 6-(1H-imidazol-2-ylmethylsulfanyl)pyridine-3,5-dicarbonitrile (LUF5834), a potent partial agonist without a ribose moiety, compared with the adenosine-derived reference agonist 2-[p-(2-carboxyethyl)phenyl-ethylamino]-5'-N-ethylcarboxamidoadenosine (CGS21680). Mutation of the orthosteric residue Phe168 to alanine abrogated the function of both agonists. However, mutation to alanine of residues Thr88 and Ser277 shown by the crystal structures to interact with the ribose group of adenosine-like ligands had no effect on the potency of LUF5834. Furthermore, alanine mutation of Asn253, which makes a hydrogen-bonding interaction with the exocyclic nitrogen of the adenine ring, had minimal effect on LUF5834 affinity but removed agonist activity of this ligand. Mutation of other residues, such as the highly conserved Trp246 or Glu13, had significant deleterious effects on the function of CGS21680 but little effect on LUF5834. In summary, our findings suggest that this class of agonist interacts with distinct residues to activate the receptor compared with classic adenosine derived agonists.

GPCR agonist binding revealed by modeling and crystallography

Despite recent progress in structural coverage of the G-protein-coupled receptor (GPCR) family, high plasticity of these membrane proteins poses additional challenges for crystallographic studies of their complexes with different classes of ligands, especially agonists. The ability to predict computationally the binding of natural and clinically relevant agonists and corresponding changes in the receptor pocket, starting from inactive GPCR structures, is therefore of great interest for understanding GPCR biology and drug action. Comparison of computational models published in 2009 and 2010 with recently determined agonist-bound structures of β-adrenergic and adenosine A(2A) receptors reveals high accuracy of the predicted agonist binding poses (0.8 Å and 1.7 Å respectively) and receptor interactions. In the case of the β(2)AR, energy-based models with limited backbone flexibility have also allowed characterization of side-chain rotations and a finite backbone shift in the pocket region as determinants of full, partial or inverse agonism. Development of accurate models of agonist binding for other GPCRs will be instrumental for functional and pharmacological studies, complementing biochemical and crystallographic techniques.

Molecular basis of ligand dissociation in β-adrenergic receptors

The important and diverse biological functions of β-adrenergic receptors (βARs) have promoted the search for compounds to stimulate or inhibit their activity. In this regard, unraveling the molecular basis of ligand binding/unbinding events is essential to understand the pharmacological properties of these G protein-coupled receptors. In this study, we use the steered molecular dynamics simulation method to describe, in atomic detail, the unbinding process of two inverse agonists, which have been recently co-crystallized with β₁ and β₂ARs subtypes, along four different channels. Our results indicate that this type of compounds likely accesses the orthosteric binding site of βARs from the extracellular water environment. Importantly, reconstruction of forces and energies from the simulations of the dissociation process suggests, for the first time, the presence of secondary binding sites located in the extracellular loops 2 and 3 and transmembrane helix 7, where ligands are transiently retained by electrostatic and Van der Waals interactions. Comparison of the residues that form these new transient allosteric binding sites in both βARs subtypes reveals the importance of non-conserved electrostatic interactions as well as conserved aromatic contacts in the early steps of the binding process.

Identification of positions in the human neuropeptide Y/peptide YY receptor Y2 that contribute to pharmacological differences between receptor subtypes

The members of the neuropeptide Y (NPY) family are key players in food-intake regulation. In humans this family consists of NPY, peptide YY (PYY) and pancreatic polypeptide (PP) which interact with distinct preference for the four receptors showing very low sequence identity, i.e. Y1, Y2, Y4 and Y5. The binding of similar peptides to these divergent receptors makes them highly interesting for mutagenesis studies. We present here a site-directed mutagenesis study of four amino acid positions in the human Y2 receptor. T(3.40) was selected based on sequence alignments both between subtypes and between species and G(2.68), L(4.60) and Q(6.55) also on previous binding studies of the corresponding positions in the Y1 receptor. The mutated receptors were characterized pharmacologically with the peptide agonists NPY, PYY, PYY(3-36), NPY(13-36) and the non-peptide antagonist BIIE0246. Interestingly, the affinity of NPY and PYY(3-36) increased for the mutants T(3.40)I and Q(6.55)A. Increased affinity was also observed for PYY to Q(6.55)A. PYY(3-36) displayed decreased affinity for G(2.68)N and L(4.60)A whereas binding of NPY(13-36) was unaffected by all mutations. The antagonist BIIE0246 showed decreased affinity for T(3.40)I, L(4.60)A and Q(6.55)A. Although all positions investigated were found important for interaction with at least one of the tested ligands the corresponding positions in hY1 seem to be of greater importance for ligand binding. Furthermore these data indicate that binding of the agonists and the antagonist differs in their points of interaction. The increase in the binding affinity observed may reflect an indirect effect caused by a conformational change of the receptor. These findings will help to improve the structural models of the human NPY receptors.

Effect of fenoterol stereochemistry on the β2 adrenergic receptor system: ligand-directed chiral recognition

The β(2) adrenergic receptor (β(2)-AR) is a model system for studying the ligand recognition process in G protein-coupled receptors. Fenoterol (FEN) is a β(2)-AR selective agonist that has two centers of chirality and exists as four stereoisomers. Radioligand binding studies determined that stereochemistry greatly influences the binding affinity. Subsequent Van't Hoff analysis shows very different thermodynamics of binding depending on the stereoconfiguration of the molecule. The binding of (S,x')-isomers is almost entirely enthalpy controlled whereas binding of (R,x')-isomers is purely entropy driven. Stereochemistry of FEN molecule also affects the coupling of the receptor to different G proteins. In a rat cardiomyocyte contractility model, (R,R')-FEN was shown to selectively activate G(s) protein signaling while the (S,R')-isomer activated both G(i) and G(s) protein. The overall data demonstrate that the chirality at the two chiral centers of the FEN molecule influences the magnitude of binding affinity, thermodynamics of local interactions within the binding site, and the global mechanism of β(2)-AR activation. Differences in thermodynamic parameters and nonuniform G-protein coupling suggest a mechanism of chiral recognition in which observed enantioselectivities arise from the interaction of the (R,x')-FEN stereoisomers with a different receptor conformation than the one with which the (S,x')-isomer interacts.

Modeling of ligand binding to G protein coupled receptors: cannabinoid CB1, CB2 and adrenergic β 2 AR

Cannabinoid and adrenergic receptors belong to the class A (similar to rhodopsin) G protein coupled receptors. Docking of agonists and antagonists to CB(1) and CB(2) cannabinoid receptors revealed the importance of a centrally located rotamer toggle switch and its possible participation in the mechanism of agonist/antagonist recognition. The switch is composed of two residues, F3.36 and W6.48, located on opposite transmembrane helices TM3 and TM6 in the central part of the membranous domain of cannabinoid receptors. The CB(1) and CB(2) receptor models were constructed based on the adenosine A(2A) receptor template. The two best scored conformations of each receptor were used for the docking procedure. In all poses (ligand-receptor conformations) characterized by the lowest ligand-receptor intermolecular energy and free energy of binding the ligand type matched the state of the rotamer toggle switch: antagonists maintained an inactive state of the switch, whereas agonists changed it. In case of agonists of β(2)AR, the (R,R) and (S,S) stereoisomers of fenoterol, the molecular dynamics simulations provided evidence of different binding modes while preserving the same average position of ligands in the binding site. The (S,S) isomer was much more labile in the binding site and only one stable hydrogen bond was created. Such dynamical binding modes may also be valid for ligands of cannabinoid receptors because of the hydrophobic nature of their ligand-receptor interactions. However, only very long molecular dynamics simulations could verify the validity of such binding modes and how they affect the process of activation.