Molecular dynamics simulations of the effect of the G-protein and diffusible ligands on the β2-adrenergic receptor

G-protein-coupled receptors have extraordinary therapeutic potential as targets for a broad spectrum of diseases. Understanding their function at the molecular level is therefore essential. A variety of crystal structures have made the investigation of the inactive receptor state possible. Recently released X-ray structures of opsin and the β₂-adrenergic receptor (β₂AR) have provided insight into the active receptor state. In addition, we have contributed to the crystal structure of an irreversible agonist-β₂ adrenoceptor complex. These extensive studies and biophysical investigations have revealed that agonist binding leads to a low-affinity conformation of the active state that is suggested to facilitate G-protein binding. The high-affinity receptor state, which promotes signal transduction, is only formed in the presence of both agonist and G-protein. Despite numerous crystal structures, it is not yet clear how ligands tune receptor dynamics and G-protein binding. We have now used molecular dynamics simulations to elucidate the distinct impact of agonist and inverse agonist on receptor conformation and G-protein binding by investigating the influence of the ligands on the structure and dynamics of a complex composed of β₂AR and the C-terminal end of the Gα_s subunit (GαCT). The simulations clearly showed that the agonist isoprenaline and the inverse agonist carazolol influence the ligand-binding site and the interaction between β₂AR and GαCT differently. Isoprenaline induced an inward motion of helix 5, whereas carazolol blocked the rearrangement of the extracellular part of the receptor. Moreover, in the presence of isoprenaline, β₂AR and GαCT form a stable interaction that is destabilized by carazolol.     

The unique ligand-binding pocket for the human prostacyclin receptor. Site-directed mutagenesis and molecular modeling

The human prostacyclin receptor is a seven-transmembrane alpha-helical G-protein coupled receptor, which plays important roles in both vascular smooth muscle relaxation as well as prevention of blood coagulation. The position of the native ligand-binding pocket for prostacyclin as well as other derivatives of the 20-carbon eicosanoid, arachidonic acid, has yet to be determined. Through the use of prostanoid receptor sequence alignments, site-directed mutagenesis, and the 2.8-A x-ray crystallographic structure of bovine rhodopsin, we have developed a three-dimensional model of the agonist-binding pocket within the seven-transmembrane (TM) domains of the human prostacyclin receptor. Upon mutation to alanine, 11 of 29 candidate residues within TM domains II, III, IV, V, and VII exhibited a marked decrease in agonist binding. Of this group, four amino acids, Arg-279 (TMVII), Phe-278 (TMVII), Tyr-75 (TMII), and Phe-95 (TMIII), were identified (via receptor amino acid sequence alignment, ligand structural comparison, and computer-assisted homology modeling) as having direct molecular interactions with ligand side-chain constituents. This binding pocket is distinct from that of the biogenic amine receptors and rhodopsin where the native ligands (also composed of a carbon ring and a carbon chain) are accommodated in an opposing direction. These findings should assist in the development of novel and highly specific ligands including selective antagonists for further molecular pharmacogenetic studies of the human prostacyclin receptor.     

Active state-like conformational elements in the beta2-AR and a photoactivated intermediate of rhodopsin identified by dynamic properties of GPCRs

G-Protein-coupled receptors (GPCRs) adopt various functionally relevant conformational states in cell signaling processes. Recently determined crystal structures of rhodopsin and the beta 2-adrenergic receptor (beta 2-AR) offer insight into previously uncharacterized active conformations, but the molecular states of these GPCRs are likely to contain both inactive and active-like conformational elements. We have identified conformational rearrangements in the dynamics of the TM7-HX8 segment that relate to the properties of the conserved NPxxY(x)5,6F motif and show that they can be used to identify active state-like conformational elements in the corresponding regions of the new structures of rhodopsin and the beta 2-AR.     

Exploring the binding site structure of the PPAR gamma ligand-binding domain by computational solvent mapping

Solvent mapping moves molecular probes, small organic molecules containing various functional groups, around the protein surface, finds favorable positions, clusters the conformations, and ranks the clusters based on the average free energy. Using at least six different solvents as probes, the probes cluster in major pockets of the functional site, providing detailed and reliable information on the amino acid residues that are important for ligand binding. Solvent mapping was applied to 12 structures of the peroxisome proliferator activated receptor gamma (PPARgamma) ligand-binding domain (LBD), including 2 structures without a ligand, 2 structures with a partial agonist, and 8 structures with a PPAR agonist bound. The analysis revealed 10 binding "hot spots", 4 in the ligand-binding pocket, 2 in the coactivator-binding region, 1 in the dimerization domain, 2 around the ligand entrance site, and 1 minor site without a known function. Mapping is a major source of information on the role and cooperativity of these sites. It shows that large portions of the ligand-binding site are already formed in the PPARgamma apostructure, but an important pocket near the AF-2 transactivation domain becomes accessible only in structures that are cocrystallized with strong agonists. Conformational changes were seen in several other sites, including one involved in the stabilization of the LBD and two others at the region of the coactivator binding. The number of probe clusters retained by these sites depends on the properties of the bound agonist, providing information on the origin of correlations between ligand and coactivator binding.     

Constitutive activity and ligand-dependent activation of the nuclear receptor CAR-insights from molecular dynamics simulations

The constitutive androstane receptor (CAR) possesses, unlike most other nuclear receptors, a pronounced basal activity in vitro whose structural basis is still not fully understood. Using comparative molecular dynamics simulations of CAR X-ray crystal structures, we evaluated the molecular basis for constitutive activity and ligand-dependent receptor activation. Our results suggest that a combination of van der Waals interactions and hydrogen bonds is required to maintain the activation helix in the active conformation also in absence of a ligand. Furthermore, we identified conformational rearrangements within the ligand-binding pocket upon agonist binding and an influence of CAR inducers pregnanedione and CITCO on the helical conformation of the activation helix. Based on the results a model for ligand-dependent CAR activation is suggested.     

Solution structure and backbone dynamics of human Raf-1 kinase inhibitor protein

Human Raf-1 kinase inhibitor protein (hRKIP) is a small multi-functional protein of 187 residues. It contains a conserved pocket, which binds a wide range of ligands from various small molecules to distinct proteins. To provide a structural basis for the ligand diversity of RKIP, we herein determined the solution structure of hRKIP, and analyzed its structural dynamics. In solution, hRKIP mainly comprises two antiparallel β sheets, two α helices and two 3₁₀ helices. NMR dynamic analysis reveals that the overall structure of hRKIP is rigid, but its C-terminal helix which is close to the ligand-binding site is mobile. In addition, residues around the ligand-binding pocket exhibit significant conformational exchange on the μs–ms timescale. Conformational flexibility may allow the ligand-binding pocket and the C-terminal helix to adopt various conformations to interact with different substrates. This work may shed light on the underlying molecular mechanisms of how hRKIP recognizes and binds diverse substrate ligands.     

Targeted molecular dynamics study of C-loop closure and channel gating in nicotinic receptors

The initial coupling between ligand binding and channel gating in the human α7 nicotinic acetylcholine receptor (nAChR) has been investigated with targeted molecular dynamics (TMD) simulation. During the simulation, eight residues at the tip of the C-loop in two alternating subunits were forced to move toward a ligand-bound conformation as captured in the crystallographic structure of acetylcholine binding protein (AChBP) in complex with carbamoylcholine. Comparison of apo- and ligand-bound AChBP structures shows only minor rearrangements distal from the ligand-binding site. In contrast, comparison of apo and TMD simulation structures of the nAChR reveals significant changes toward the bottom of the ligand-binding domain. These structural rearrangements are subsequently translated to the pore domain, leading to a partly open channel within 4 ns of TMD simulation. Furthermore, we confirmed that two highly conserved residue pairs, one located near the ligand-binding pocket (Lys145 and Tyr188), and the other located toward the bottom of the ligand-binding domain (Arg206 and Glu45), are likely to play important roles in coupling agonist binding to channel gating. Overall, our simulations suggest that gating movements of the α7 receptor may involve relatively small structural changes within the ligand-binding domain, implying that the gating transition is energy-efficient and can be easily modulated by agonist binding/unbinding.     

On the binding determinants of the glutamate agonist with the glutamate receptor ligand binding domain

Ionotropic glutamate receptors (GluRs) are ligand-gated membrane channel proteins found in the central neural system that mediate a fast excitatory response of neurons. In this paper, we report theoretical analysis of the ligand-protein interactions in the binding pocket of the S1S2 (ligand binding) domain of the GluR2 receptor in the closed conformation. By utilizing several theoretical methods ranging from continuum electrostatics to all-atom molecular dynamics simulations and quantum chemical calculations, we were able to characterize in detail glutamate agonist binding to the wild-type and E705D mutant proteins. A theoretical model of the protein-ligand interactions is validated via direct comparison of theoretical and Fourier transform infrared spectroscopy (FTIR) measured frequency shifts of the ligand's carboxylate group vibrations [Jayaraman et al. (2000) Biochemistry 39, 8693-8697; Cheng et al. (2002) Biochemistry 41, 1602-1608]. A detailed picture of the interactions in the binding site is inferred by analyzing contributions to vibrational frequencies produced by protein residues forming the ligand-binding pocket. The role of mobility and hydrogen-bonding network of water in the ligand-binding pocket and the contribution of protein residues exposed in the binding pocket to the binding and selectivity of the ligand are discussed. It is demonstrated that the molecular surface of the protein in the ligand-free state has mainly positive electrostatic potential attractive to the negatively charged ligand, and the potential produced by the protein in the ligand-binding pocket in the closed state is complementary to the distribution of the electrostatic potential produced by the ligand itself. Such charge complementarity ensures specificity to the unique charge distribution of the ligand.     

Light-driven activation of beta 2-adrenergic receptor signaling by a chimeric rhodopsin containing the beta 2-adrenergic receptor cytoplasmic loops

Structure-function studies of rhodopsin indicate that both intradiscal and transmembrane (TM) domains are required for retinal binding and subsequent light-induced structural changes in the cytoplasmic domain. Further, a hypothesis involving a common mechanism for activation of G-protein-coupled receptor (GPCR) has been proposed. To test this hypothesis, chimeric receptors were required in which the cytoplasmic domains of rhodopsin were replaced with those of the beta(2)-adrenergic receptor (beta(2)-AR). Their preparation required identification of the boundaries between the TM domain of rhodopsin and the cytoplasmic domain of the beta(2)-AR necessary for formation of the rhodopsin chromophore and its activation by light and subsequent optimal activation of beta(2)-AR signaling. Chimeric receptors were constructed in which the cytoplasmic loops of rhodopsin were replaced one at a time and in combination. In these replacements, size of the third cytoplasmic (EF) loop critically determined the extent of chromophore formation, its stability, and subsequent signal transduction specificity. All the EF loop replacements showed significant decreases in transducin activation, while only minor effects were observed by replacements of the CD and AB loops. Light-dependent activation of beta(2)-AR leading to Galphas signaling was observed only for the EF2 chimera, and its activation was further enhanced by replacements of the other loops. The results demonstrate coupling between light-induced conformational changes occurring in the transmembrane domain of rhodopsin and the cytoplasmic domain of the beta(2)-AR.     

GluR2 ligand-binding core complexes: importance of the isoxazolol moiety and 5-substituent for the binding mode of AMPA-type agonists

X-ray structures of the GluR2 ligand-binding core in complex with (S)-Des-Me-AMPA and in the presence and absence of zinc ions have been determined. (S)-Des-Me-AMPA, which is devoid of a substituent in the 5-position of the isoxazolol ring, only has limited interactions with the partly hydrophobic pocket of the ligand-binding site, and adopts an AMPA-like binding mode. The structures, in comparison with other agonist complex structures, disclose the relative importance of the isoxazolol ring and of the substituent in the 5-position for the mode of binding. A relationship appears to exist between the extent of interaction of the ligand with the hydrophobic pocket and the affinity of the ligand.     

Protein dynamics and conformational selection in bidirectional signal transduction

Protein conformational dynamics simultaneously allow promiscuity and specificity in binding. The multiple conformations of the free EphA4 ligand-binding domain observed in two new EphA4 crystal structures provide a unique insight into the conformational dynamics of EphA4 and its signaling pathways. The heterogeneous ensemble and loop dynamics explain how the EphA4 receptor is able to bind multiple A- and B-ephrin ligands and small molecules via conformational selection, which helps to fine-tune cellular signal response in both receptor and ligand cells.     

Structures of Aplysia AChBP complexes with nicotinic agonists and antagonists reveal distinctive binding interfaces and conformations

Upon ligand binding at the subunit interfaces, the extracellular domain of the nicotinic acetylcholine receptor undergoes conformational changes, and agonist binding allosterically triggers opening of the ion channel. The soluble acetylcholine-binding protein (AChBP) from snail has been shown to be a structural and functional surrogate of the ligand-binding domain (LBD) of the receptor. Yet, individual AChBP species display disparate affinities for nicotinic ligands. The crystal structure of AChBP from Aplysia californica in the apo form reveals a more open loop C and distinctive positions for other surface loops, compared with previous structures. Analysis of Aplysia AChBP complexes with nicotinic ligands shows that loop C, which does not significantly change conformation upon binding of the antagonist, methyllycaconitine, further opens to accommodate the peptidic antagonist, alpha-conotoxin ImI, but wraps around the agonists lobeline and epibatidine. The structures also reveal extended and nonoverlapping interaction surfaces for the two antagonists, outside the binding loci for agonists. This comprehensive set of structures reflects a dynamic template for delineating further conformational changes of the LBD of the nicotinic receptor.     

Differential hydrogen/deuterium exchange mass spectrometry analysis of protein-ligand interactions

Functional regulation of ligand-activated receptors is driven by alterations in the conformational dynamics of the protein upon ligand binding. Differential hydrogen/deuterium exchange (HDX) coupled with mass spectrometry has emerged as a rapid and sensitive approach for characterization of perturbations in conformational dynamics of proteins following ligand binding. While this technique is sensitive to detecting ligand interactions and alterations in receptor dynamics, it also can provide important mechanistic insights into ligand regulation. For example, HDX has been used to determine a novel mechanism of ligand activation of the nuclear receptor peroxisome proliferator activated receptor-γ, perform detailed analyses of binding modes of ligands within the ligand-binding pocket of two estrogen receptor isoforms, providing insight into selectivity, and helped classify different types of estrogen receptor-α ligands by correlating their pharmacology with the way they interact with the receptor based solely on hierarchical clustering of receptor HDX signatures. Beyond small-molecule-receptor interactions, this technique has also been applied to study protein-protein complexes, such as mapping antibody-antigen interactions. In this article, we summarize the current state of the differential HDX approaches and the future outlook. We summarize how HDX analysis of protein-ligand interactions has had an impact on biology and drug discovery.     

Structural insights into agonist-induced activation of G-protein-coupled receptors

Recent years have seen tremendous breakthroughs in structure determination of G-protein-coupled receptors (GPCRs). In 2011, two agonist-bound active-state structures of rhodopsin have been published. Together with structures of several rhodopsin activation intermediates and a wealth of biochemical and spectroscopic information, they provide a unique structural framework on which to understand GPCR activation. Here we use this framework to compare the recent crystal structures of the agonist-bound active states of the β(2) adrenergic receptor (β(2)AR) and the A(2A) adenosine receptor (A(2A)AR). While activation of these three GPCRs results in rearrangements of TM5 and TM6, the extent of this conformational change varies considerably. Displacements of the cytoplasmic side of TM6 ranges between 3 and 8? depending on whether selective stabilizers of the active conformation are used (i.e. a G-protein peptide in the case of rhodopsin or a conformationally selective nanobody in the case of the β(2)AR) or not (A(2A)AR). The agonist-induced conformational changes in the ligand-binding pocket are largely receptor specific due to the different chemical nature of the agonists. However, several similarities can be observed, including a relocation of conserved residues W6.48 and F6.44 towards L5.51 and P5.50, and of I/L3.40 away from P5.50. This transmission switch links agonist binding to the movement of TM5 and TM6 through the rearrangement of the TM3-TM5-TM6 interface, and possibly constitutes a common theme of GPCR activation.     

Insights into signaling from the beta2-adrenergic receptor structure

With more than 800 members, the G protein-coupled receptor family constitutes the largest group of membrane proteins involved in signal transduction. Until the end of last year, high-resolution three-dimensional structures were available for only one of them--the light receptor rhodopsin. Recently the structure of the beta(2)-adrenergic receptor has been obtained, and it revealed interesting differences with the structure of rhodopsin. Analyses of these differences raise important questions about the binding modes of diffusible ligands in the receptor and allow formulation of testable hypotheses about the structural determinants linking drug binding to specific signaling responses. The three-dimensional structure derived from the beta(2)-adrenergic receptor crystal has been used to virtually dock ligands with distinct activities. The different binding modes of these ligands, which correlated with their reported efficacy profiles, suggest that it could be possible to predict the structural determinants of drug signaling efficacies.     

Model structures of the N-methyl-D-aspartate receptor subunit NR1 explain the molecular recognition of agonist and antagonist ligands

Molecular models of the ligand-binding domain of N-methyl-d-aspartate subunit R1 (NR1) were made using the published crystal structures of rat glutamate receptor B (GluRB), the bacterial glutamate receptor (GluR0), and the glutamine-binding protein (QBP) of Escherichia coli. Separate models of NR1 were built to represent the ligand-binding conformation for agonist (glycine, d- and l-isomers of serine and alanine, and the partial agonist ligand d-cycloserine) and antagonist (5,7-dichloro-4-oxo-1,4-dihydroquinoline-2-carboxylic acid (DCKA) and E-3-(2-phenyl-2-carboxyethenyl)-4,6-dichloro-1-H-indole-2-carboxylic acid (MDL 105,519)) ligands. Side-chain conformations of residues within the NR1 ligand-binding site were selected that optimized the hydrophobic packing and hydrogen bonding among residues, while taking into account published data comparing receptor mutants with wild-type NR1. Ligands docked to the model structures provide a rational explanation for the observed differences in binding affinity and receptor activation among agonist and antagonist ligands. NR1 prefers smaller ligands (glycine, serine, and alanine) in comparison with GluRB and GluR0 that bind l-glutamate: the bulky side chain of W731 in NR1 dramatically reduces the size of the ligand-binding site, functioning to selectively restrict recognition to glycine and the d-isomers of serine and alanine. Nevertheless, many of the interactions seen for ligands bound to GluRB, GluR0, and periplasmic-binding proteins are present for the ligands docked to the model structures of NR1.     

GPCRs through the keyhole: the role of protein flexibility in ligand binding to β-adrenoceptors

G protein-coupled receptors (GPCRs) are proteins of pharmaceutical importance, with over 30% of all drugs in clinical use targeting them. Increasing numbers of X-ray crystal (XRC) structures of GPCRs offer a wealth of data relating to ligand binding. For the β-adrenoceptors (β-ARs), XRC structures are available for human β₂- and turkey β₁-subtypes, in complexes with a range of ligands. While these structures provide insight into the origins of ligand structure-activity relationships (SARs), questions remain. The ligands in all published complexed XRC structures lack extensive substitution, with no obvious way the ligand-binding site can accommodate β₁-AR-selective antagonists with extended side-chains para- to the common aryloxypropanolamine pharmacophore. Using standard computational docking tools with such ligands generally returns poses that fail to explain known SARs. Application of our Active Site Pressurisation modelling method to β-AR XRC structures and homology models, however, reveals a dynamic area in the ligand-binding pocket that, through minor changes in amino acid side chain orientations, opens a fissure between transmembrane helices H4 and H5, exposing intra-membrane space. This fissure, which we term the “keyhole”, is ideally located to accommodate extended moieties present in many high-affinity β₁-AR-selective ligands, allowing the rest of the ligand structure to adopt a canonical pose in the orthosteric binding site. We propose the keyhole may be a feature of both β₁- and β₂-ARs, but that subtle structural differences exist between the two, contributing to subtype-selectivity. This has consequences for the rational design of future generations of subtype-selective ligands for these therapeutically important targets.