Ligand-induced Conformational Changes Research Articles

C-terminal Movement during Gating in Cyclic Nucleotide-modulated Channels

Activation of cyclic nucleotide-modulated channels such as CNG and HCN channels is promoted by ligand-induced conformational changes in their C-terminal regions. The primary intersubunit interface of these C termini includes two salt bridges per subunit, formed between three residues (one positively charged and two negatively charged amino acids) that we term the SB triad. We previously hypothesized that the SB triad is formed in the closed channel and breaks when the channel opens. Here we tested this hypothesis by dynamically manipulating the SB triad in functioning CNGA1 channels. Reversing the charge at positions Arg-431 and Glu-462, two of the SB triad residues, by either mutation or application of charged reagents increased the favorability of channel opening. To determine how a charge reversal mutation in the SB triad structurally affects the channel, we solved the crystal structure of the HCN2 C-terminal region with the equivalent E462R mutation. The backbone structure of this mutant was very similar to that of wild type, but the SB triad was rearranged such that both salt bridges did not always form simultaneously, suggesting a mechanism for the increased ease of opening of the mutant channels. To prevent movement in the SB triad, we tethered two components of the SB triad region together with cysteine-reactive cross-linkers. Preventing normal movement of the SB triad region with short cross-linkers inhibited channel opening, whereas longer cross-linkers did not. These results support our hypothesis that the SB triad forms in the closed channel and indicate that this region expands as the channel opens.

Structural Insight into the Mechanism of Activation of the Toll Receptor by the Dimeric Ligand Spätzle

The Drosophila Toll receptor, which functions in both embryonic patterning and innate immunity to fungi and Gram-positive bacteria, is activated by a dimeric cytokine ligand, Spätzle (Spz). Previous studies have suggested that one Spz cross-links two Toll receptor molecules to form an activated complex. Here we report electron microscopy structures of the Toll ectodomain in the absence and presence of Spz. Contrary to expectations, Spz does not directly cross-link two Toll ectodomains. Instead, Spz binding at the N-terminal end of Toll predominantly induces the formation of a 2:2 complex, with two sites of interaction between the ectodomain chains, one located near to the N terminus of the solenoid and the other between the C-terminal juxtamembrane sequences. Moreover, Toll undergoes a ligand-induced conformational change, becoming more tightly curved than in the apo form. The unexpected 2:2 complex was confirmed by mass spectrometry under native conditions. These results suggest that activation of Toll is an allosteric mechanism induced by an end-on binding mode of its ligand Spz.

Molecular dynamics simulations of ligand-induced backbone conformational changes in the binding site of the periplasmic lysine-, arginine-, ornithine-binding protein

The periplasmic lysine-, arginine-, ornithine-binding protein (LAOBP) traps its ligands by a large hinge bending movement between two globular domains. The overall geometry of the binding site remains largely unchanged between the open (unliganded) and closed (liganded) forms, with only a small number of residues exhibiting limited movement of their side chains. However, in the case of the ornithine-bound structure, the backbone peptide bond between Asp11 and Thr12 undergoes a large rotation. Molecular dynamics simulations have been used to investigate the origin and mechanism of this backbone movement. Simulations allowing flexibility of a limited region and of the whole binding site, with and without bound ligands, suggest that this conformational change is induced by the binding of ornithine, leading to the stabilisation of an energetically favourable alternative conformation.

Structure-Function Relationships of the C-Terminal End of the Saccharomyces cerevisiae ADP/ATP Carrier Isoform 2

The adenine nucleotide carrier (Ancp) catalyzes the transport of ADP and ATP across the mitochondrial inner membrane, thus playing an essential role in the cellular energy metabolism. Two regions of Anc2p from Saccharomyces cerevisiae are specifically photolabeled using a photoactivable ADP derivative; they are the central matrix loop, m2, and the C-terminal end. To get more insights into the structure-function relationships of the C-terminal region during nucleotide transport, we have developed two independent approaches. In the first we have deleted the last eight amino acids of Anc2p (Anc2pDeltaCter) and demonstrated that the C-terminal end of Anc2p plays an essential role in yeast growth on a non-fermentable carbon source. This resulted from impaired nucleotide binding properties of the Anc2pDeltaCter variant in line with conversion of ADP binding sites from high to low affinity. In the second we probed the ligand-induced conformational changes of Anc2p C-terminal end (i) by assessing its accessibility to anti-C-terminal antibodies and (ii) by measuring intrinsic fluorescence changes of an Anc2p mutant containing only one tryptophan residue located at its C-terminal end (Anc2p3Y-u). We show that the C-terminal region is no further accessible to antibodies when Anc2p binds non-transportable analogues of ADP. Besides, Trp-316 fluorescence is highly increased upon ligand binding, suggesting large conformational changes. Taken together, our results highlight the involvement of the Anc2p C-terminal region in nucleotide recognition, binding, and transport.

Structure and Dynamics of the Human beta 2 Adrenoceptor

G protein coupled receptors (GPCRs) are remarkably versatile signaling molecules. The beta 2 adrenoceptor (beta2AR) is a prototypical Family A GPCR that mediates physiologic responses to adrenaline and noradrenaline. The function of the beta2AR can be modulated by a spectrum of synthetic ligands ranging from full agonists to inverse agonists. We have used crystallography to determine the three‐dimensional structure of the beta2AR, and fluorescence spectroscopy to map ligand‐induced conformational changes and characterize the structure of beta2AR dimers. I will discuss what we these studies have taught us about the structural basis of beta2ARs function.

Ligand-Stabilized Conformational States of Human β2 Adrenergic Receptor: Insight into G-Protein-Coupled Receptor Activation

G-protein-coupled receptors (GPCRs) are known to exist in dynamic equilibrium between inactive- and several active-state conformations, even in the absence of a ligand. Recent experimental studies on the β 2 adrenergic receptor ( β 2AR) indicate that structurally different ligands with varying efficacies trigger distinct conformational changes and stabilize different receptor conformations. We have developed a computational method to study the ligand-induced rotational orientation changes in the transmembrane helices of GPCRs. This method involves a systematic spanning of the rotational orientation of the transmembrane helices (TMs) that are in the vicinity of the ligand for predicting the helical rotations that occur on ligand binding. The predicted ligand-stabilized receptor conformations are characterized by a simultaneous lowering of the ligand binding energy and a significant gain in interhelical and receptor-ligand hydrogen bonds. Using the β 2AR as a model, we show that the receptor conformational state depends on the structure and efficacy of the ligand for a given signaling pathway. We have studied the ligand-stabilized receptor conformations of five different ligands, a full agonist, norepinephrine; a partial agonist, salbutamol; a weak partial agonist, dopamine; a very weak agonist, catechol; and an inverse agonist, ICI-115881. The predicted ligand-stabilized receptor models correlate well with the experimentally observed conformational switches in β 2AR, namely, the breaking of the ionic lock between R131 3.50 at the intracellular end of TM3 (part of the DRY motif) and E268 6.30 on TM6, and the rotamer toggle switch on W286 6.48 on TM6. In agreement with trp-bimane quenching experiments, we found that norepinephrine and dopamine break the ionic lock and engage the rotamer toggle switch, whereas salbutamol, a noncatechol partial agonist only breaks the ionic lock, and the weak agonist catechol only engages the rotamer toggle switch. Norepinephrine and dopamine occupy the same binding region, between TM3, TM5, and TM6, whereas the binding site of salbutamol is shifted toward TM4. Catechol binds deeper into the protein cavity compared to the other ligands, making contact with TM5 and TM6. A part of the catechol binding site overlaps with those of dopamine and norepinephrine but not with that of salbutamol. Virtual ligand screening on 10,060 ligands on the norepinephrine-stabilized receptor conformation shows an enrichment of 38% compared to ligand unbound receptor conformation. These results show that ligand-induced conformational changes are important for developing functionally specific drugs that will stabilize a particular receptor conformation. These studies represent the first step toward a more universally applicable computational method for studying ligand efficacy and GPCR activation.

Complex Ligand-Induced Conformational Changes in tRNAAsp Revealed by Single-Nucleotide Resolution SHAPE Chemistry

RNA conformation is both highly dependent on and sensitive to the presence of charged ligands. Mono- and divalent ions stabilize the native fold of RNA, whereas other polyvalent cationic ligands can act to both stabilize or disrupt native RNA structure. In this work, we analyze the effects of two ligands, Mg (2+) and tobramycin, on the folding of S. cerevisiae tRNA (Asp) transcripts using single nucleotide resolution SHAPE chemistry. Surprisingly, reducing the Mg (2+) concentration favors a structural rearrangement in which the D- and variable loops pair. The tobramycin polycation binds to loops in tRNA (Asp) and induces RNA unfolding in two distinct transitions: the loss of tertiary interactions between the T- and D-loops followed by complete unfolding of the D-stem. Although Mg (2+) and tobramycin are relatively simple ligands, both modulate tRNA (Asp) folding in unanticipatedly complex ways, neither of which is consistent with simple hierarchical folding or unfolding of this RNA. Monitoring the structural consequences of ligand binding to RNA at single nucleotide resolution makes it possible to define intermediate structures that contribute to the complex energy landscapes often observed for RNA folding processes and lays the groundwork for a significantly improved understanding of the interactions between RNA and its solution environment.

Molecular modelling of human CYP2D6 and molecular docking of a series of ajmalicine- and quinidine-like inhibitors

3D-models were created and refined for CYP2D6 and for its complexes with ajmalicine and quinidine. The influence of the conformation of the enzyme active site on its interaction with ligands was evaluated by performing three series of molecular docking on selected ajmalicine- and quinidine-like inhibitors. The results suggested that the experimental binding values of ajmalicine- and quinidine-like inhibitors better fit with the energetic terms derived from their interaction with structures of CYP2D6 obtained by, respectively, optimizing the ajmalicine/CYP2D6 and the quinidine/CYP2D6 complexes, rather than exploiting the 3D-strucure of the enzyme not subjected to a ligand-induced conformational change. It suggests the relevance of induced-fit phenomena in the biological system of interest.

Modeling and Selection of Flexible Proteins for Structure‐Based Drug Design: Backbone and Side Chain Movements in p38 MAPK

Receptor rearrangement upon ligand binding (induced fit) is a major stumbling block in docking and virtual screening. Even though numerous studies have stressed the importance of including protein flexibility in ligand docking, currently available methods provide only a partial solution to the problem. Most of these methods, being computer intensive, are often impractical to use in actual drug discovery settings. We had earlier shown that ligand-induced receptor side-chain conformational changes could be modeled statistically using data on known receptor-ligand complexes. In this paper, we show that a similar approach can be used to model more complex changes like backbone flips and loop movements. We have used p38 MAPK as a test case and have shown that a few simple structural features of ligands are sufficient to predict the induced variation in receptor conformations. Rigorous validation, both by internal resampling methods and on an external test set, corroborates this finding and demonstrates the robustness of the models. We have also compared our results with those from an earlier molecular dynamics simulation study on DFG loop conformations of p38 MAPK, and found that the results matched in the two cases. Our statistical approach enables one to predict the final ligand-induced conformation of the active site of a protein, based on a few ligand properties, prior to docking the ligand. We can do this without having to trace the step-by-step process by which this state is arrived at (as in molecular dynamics simulations), thereby drastically reducing computational effort.

Gα-Subunits Differentially Alter the Conformation and Agonist Affinity of κ-Opioid Receptors

Although ligand-induced conformational changes in G protein-coupled receptors (GPCRs) are well-documented, there is little direct evidence for G protein-induced changes in GPCR conformation. To investigate this possibility, the effects of overexpressing Galpha-subunits (Galpha16 or Galphai2) with the kappa-opioid receptor (KOR) were examined. The changes in KOR conformation were subequently examined via the substituted cysteine accessibility method (SCAM) in transmembrane domains 6 (TM6) and 7 (TM7) and extracellular loop 2 (EL2). Significant conformational changes were observed on TM7, the extracellular portion of TM6, and EL2. Seven SCAM-sensitive residues (S3107.33, F3147.37, and I3167.39 to Y3207.43) on TM7 presented a cluster pattern when the KOR was exposed to baseline amounts of G protein, and additional residues became sensitive upon overexpression of various G proteins. In TM7, S3117.34 and N3267.49 were found to be sensitive in Galpha16-overexpressed cells and Y3137.36, N3227.45, S3237.46, and L3297.52 in Galphai2-overexpressed cells. In addition, the degree of sensitivity for various TM7 residues was augmented, especially in Galphai2-overexpressed cells. A similar phenomenon was also observed for residues in TM6 and EL2. In addition to an enhanced sensitivity of certain residues, our findings also indicated that a slight rotation was predicted to occur in the upper part of TM7 upon G protein overexpression. These relatively modest conformational changes engendered by G protein overexpression had both profound and differential effects on the abilities of agonists to bind to KOR. These data are significant because they demonstrate that Galpha-subunits differentially modulate the conformation and agonist affinity of a prototypical GPCR.

Ligand Binding Induces a Conformational Change in ifnar1 that Is Propagated to Its Membrane-Proximal Domain

The type I interferon (IFN) receptor plays a key role in innate immunity against viral and bacterial infections. Here, we show by intramolecular Förster resonance energy transfer spectroscopy that ligand binding induces substantial conformational changes in the ectodomain of ifnar1 (ifnar1-EC). Binding of IFNα2 and IFNβ induce very similar conformations of ifnar1, which were confirmed by single-particle electron microscopy analysis of the ternary complexes formed by IFNα2 or IFNβ with the two receptor subunits ifnar1-EC and ifnar2-EC. Photo-induced electron-transfer-based fluorescence quenching and single-molecule fluorescence lifetime measurements revealed that the ligand-induced conformational change in the membrane-distal domains of ifnar1-EC is propagated to its membrane-proximal domain, which is not involved in ligand recognition but is essential for signal activation. Temperature-dependent ligand binding studies as well as stopped-flow fluorescence experiments corroborated a multistep conformational change in ifnar1 upon ligand binding. Our results thus suggest that the relatively intricate architecture of the type I IFN receptor complex is designed to propagate the ligand binding event to and possibly even across the membrane by conformational changes.

Predominance of activated EGFR higher-order oligomers on the cell surface

The epidermal growth factor receptor (EGFR) kinase is generally considered to be activated by either ligand-induced dimerisation or a ligand-induced conformational change within pre-formed dimers. We report the relationship between ligand-induced higher-order EGFR oligomerization and EGFR phosphorylation on the surface of intact cells. We have combined lifetime-detected Förster resonance energy transfer, as a probe of the receptor phosphorylation state and image correlation spectroscopy, to extract the relative association state of activated versus unactivated EGFR, to determine the ratio of the average number of receptors for active (phosphorylated) and inactive clusters. There are at least four times as many receptors in the ligand-induced active clusters than inactive clusters. Contrary to the prevailing view that the EGFR dimer is the predominant, active form, our data determine that higher-order EGFR oligomers are the dominant species associated with the ligand activated EGFR tyrosine kinase.

The short length of the extracellular domain of ζ is crucial for T cell antigen receptor function

The T cell antigen receptor (TCR-CD3) consists of the pMHC-binding TCRαβ heterodimer and the signalling dimers CD3δɛ, CD3γɛ and ζζ. The very short length of the extracellular domain (EC) of the ζ chain is preserved through evolution, however a rational explanation for this observation has not been elucidated. Here, we show that TCR-CD3 assembly is clearly defective when the murine ζ EC domain is artificially enlarged. Under these conditions, the TCR-CD3 complex is super-competent in transducing activation signals upon engagement. Furthermore, the TCR-CD3 complexes containing enlarged ζ EC domains underwent ligand-induced conformation changes with higher efficiency than TCR-CD3 complexes with an unmodified ζ EC domain. Together these data suggest that a short ζ EC domain is needed to correctly assemble the TCR-CD3 complex. When this domain is enlarged, the resulting TCR-CD3 complex is distorted leading to a hyperactive phenotype and enhanced T cell activation.

Ligand-induced Conformational Changes and Conformational Dynamics in the Solution Structure of the Lactose Repressor Protein

We present here the results of a series of small-angle X-ray scattering studies aimed at understanding the role of conformational changes and structural flexibility in DNA binding and allosteric signaling in a bacterial transcription regulator, lactose repressor protein (LacI). Experiments were designed to detect possible conformational changes that occur when LacI binds either DNA or the inducer IPTG, or both. Our studies included the native LacI dimer of homodimers and a dimeric variant (R3), enabling us to probe conformational changes within the homodimers and distinguish them from those involving changes in the homodimer–homodimer relationships. The scattering data indicate that removal of operator DNA (oDNA) from R3 results in an unfolding and extension of the hinge helix that connects the LacI regulatory and DNA-binding domains. In contrast, only very subtle conformational changes occur in the R3 dimer–oDNA complex upon IPTG binding, indicative of small adjustments in the orientations of domains and/or subdomains within the structure. The binding of IPTG to native (tetrameric) LacI–oDNA complexes also appears to facilitate a modest change in the average homodimer–homodimer disposition. Notably, the crystal structure of the native LacI–oDNA complex differs significantly from the average solution conformation. The solution scattering data are best fit by an ensemble of structures that includes (1) ∼60% of the V-shaped dimer of homodimers observed in the crystal structure and (2) ∼40% of molecules with more “open” forms, such as those generated when the homodimers move with respect to each other about the tetramerization domain. In gene regulation, such a flexible LacI would be beneficial for the interaction of its two DNA-binding domains, positioned at the tips of the V, with the required two of three LacI operators needed for full repression.

A Domain Linking the AMPA Receptor Agonist Binding Site to the Ion Pore Controls Gating and CauseslurcherProperties when Mutated

Ionotropic, AMPA-type glutamate receptors (GluRs) critically shape excitatory synaptic signals in the CNS. Ligand binding induces conformational changes in the glutamate-binding domain of the receptors that are converted into opening of the channel pore via three short linker sequences, a process referred to as gating. Although crystallization of the glutamate-binding domain and structural models of the ion pore advanced our understanding of ligand-binding dynamics and pore movements, the allosteric coupling of both events by the short linkers has not been described in detail. To study the role of the linkers in gating GluR1, we transplanted them between different GluRs and examined the electrophysiological properties of the resulting chimeric receptors in Xenopus laevis oocytes and HEK293 cells. We found that all three linkers decisively affect receptor functionality, agonist potency, and desensitization. One linker chimera was nondesensitizing and exhibited strongly increased agonist potencies, while fluxing ions even in the absence of agonist, similar to properties reported for the GluR1 lurcher mutation. Combining this new lurcher-like linker chimera with the original lurcher mutation allowed us to reassess the effect of lurcher on GluR1 gating properties. The observed differential but interdependent influence of linker and lurcher mutations on receptor properties suggests that the linkers are part of a fine-tuned structural element that normally stabilizes the closed ion pore. We propose that lurcher-like mutations act by disrupting this element such that ligand-induced conformational changes are not necessarily required to gate the channel.

Sodium-Dependent Reorganization of the Sugar-Binding Site of SGLT1

The sodium-dependent glucose cotransporter SGLT1 undergoes a series of voltage- and ligand-induced conformational changes that underlie the cotransport mechanism. In this study we describe how the binding of external Na changes the conformation of the sugar-binding domain, exposing residues that are involved in sugar recognition to the external environment. We constructed 15 individual Cys mutants in the four transmembrane helices (TMHs) that form the sugar binding and translocation domain. Each mutant was functionally characterized for transport kinetics and substrate specificity. Identification of interactions between mutated residues and hydroxyls on the pyranose ring was assessed by comparing the affinities of deoxy sugars to those of glucose. We determined conformation-dependent accessibility to the mutated residues by both a traditional substituted cysteine accessibility method (SCAM) and a new fluorescence binding assay. These data were integrated to orient the helices and construct a framework of residues that comprise the external sugar binding site. We present evidence that R499, Q457, and T460 play a direct role in sugar recognition and that five other residues are indirectly involved in transport. Arranging the four TMHs to account for Na-dependent accessibility and potential for sugar interaction allows us to propose a testable model for the SGLT1 sugar binding site.

Conformational selection in silico: Loop latching motions and ligand binding in enzymes

Ligand binding frequently induces significant conformational changes in a protein receptor. Understanding and predicting such conformational changes represent an important challenge for computational biology, including applications to structure-based drug design. We describe an approach to this problem based on the assumption that the holo state is at least transiently populated in the absence of a ligand; this hypothesis has been referred to as "conformational selection." Here, we apply a method that tests this hypothesis on a challenging class of ligand-induced conformational changes, which we refer to as loop latching: the closing of a loop around an active site that sequesters the ligand from solvent. The method uses a combination of replica exchange molecular dynamics and a loop prediction algorithm to generate low-energy loop structures, and docking to select the conformation appropriate for binding a particular ligand. On a test set of six proteins, it yields loop structures including hololike conformations, generally below 2 A RMSD from the liganded structure, for loops that span up to 15 residues. Docking serves as a stringent test of the predictions. In five of the six cases, the predicted loop conformations improve the ranks of cognate ligands relative to using the apo structure, although the results remain, in most cases, significantly worse than using a holo structure. The poses of the cognate ligands are correct in four of the six test cases, while they are correct for five of the six using a holo structure.

Complex, unusual conformational changes in kidney betaine aldehyde dehydrogenase suggested by chemical modification with disulfiram

The NAD +-dependent animal betaine aldehyde dehydrogenases participate in the biosynthesis of glycine betaine and carnitine, as well as in polyamines catabolism. We studied the kinetics of inactivation of the porcine kidney enzyme (pkBADH) by the drug disulfiram, a thiol-reagent, with the double aim of exploring the enzyme dynamics and investigating whether it could be an in vivo target of disulfiram. Both inactivation by disulfiram and reactivation by reductants were biphasic processes with equal limiting amplitudes. Under certain conditions half of the enzyme activity became resistant to disulfiram inactivation. NAD + protected almost 100% at 10 μM but only 50% at 5 mM, and vice versa if the enzyme was pre-incubated with NAD + before the chemical modification. NADH, betaine aldehyde, and glycine betaine also afforded greater protection after pre-incubation with the enzyme than without pre-incubation. Together, these findings suggest two kinds of active sites in this seemingly homotetrameric enzyme, and complex, unusual ligand-induced conformational changes. In addition, they indicate that, in vivo, pkBADH is most likely protected against disulfiram inactivation.

Structure‐based design of robust glucose biosensors using a Thermotoga maritima periplasmic glucose‐binding protein

We report the design and engineering of a robust, reagentless fluorescent glucose biosensor based on the periplasmic glucose-binding protein obtained from Thermotoga maritima (tmGBP). The gene for this protein was cloned from genomic DNA and overexpressed in Escherichia coli, the identity of its cognate sugar was confirmed, ligand binding was studied, and the structure of its glucose complex was solved to 1.7 Angstrom resolution by X-ray crystallography. TmGBP is specific for glucose and exhibits high thermostability (midpoint of thermal denaturation is 119 +/- 1 degrees C and 144 +/- 2 degrees C in the absence and presence of 1 mM glucose, respectively). A series of fluorescent conjugates was constructed by coupling single, environmentally sensitive fluorophores to unique cysteines introduced by site-specific mutagenesis at positions predicted to be responsive to ligand-induced conformational changes based on the structure. These conjugates were screened to identify engineered tmGBPs that function as reagentless fluorescent glucose biosensors. The Y13C*Cy5 conjugate is bright, gives a large response to glucose over concentration ranges appropriate for in vivo monitoring of blood glucose levels (1-30 mM), and can be immobilized in an orientation-specific manner in microtiter plates to give a reversible response to glucose. The immobilized protein retains its response after long-term storage at room temperature.

Distinct Structural Changes in a G Protein-coupled Receptor Caused by Different Classes of Agonist Ligands

The activity of G protein-coupled receptors can be modulated by different classes of ligands, including agonists that promote receptor signaling and inverse agonists that reduce basal receptor activity. The conformational changes in receptor structure induced by different agonist ligands are not well understood at present. In this study, we employed an in situ disulfide cross-linking strategy to monitor ligand-induced conformational changes in a series of cysteine-substituted mutant M(3) muscarinic acetylcholine receptors. The observed disulfide cross-linking patterns indicated that muscarinic agonists trigger a separation of the N-terminal segment of the cytoplasmic tail (helix 8) from the cytoplasmic end of transmembrane domain I. In contrast, inverse muscarinic agonists were found to increase the proximity between these two receptor regions. These findings provide a structural basis for the opposing biological effects of muscarinic agonists and inverse agonists. This study also provides the first piece of direct structural information as to how the conformations induced by these two functionally different classes of ligands differ at the molecular level. Given the high degree of structural homology found among most G protein-coupled receptors, our findings should be of broad general relevance.