From signal perception to signal transduction: ligand‐induced dimeric switch of DctB sensory domain in solution

Sinorhizobium meliloti DctB is a typical transmembrane sensory histidine kinase, which senses C₄‐dicarboxylic acids (DCA) and regulates the expression of DctA, the DCA transporter. We previously reported the crystal structures of its periplasmic sensory domain (DctBp) in apo and succinate‐bound states, and these structures showed dramatic conformational changes at dimeric level. Here we show a ligand‐induced dimeric switch in solution and a strong correlation between DctBp's dimerization states and the in vivo activities of DctB. Using site‐directed mutagenesis, we identify important determinants for signal perception and transduction. Specifically, we show that the ligand‐binding pocket is essential for DCA‐induced ‘on’ activity of DctB. Mutations at different sections of DctBp's dimerization interface can lock full‐length DctB at either ‘on’ or ‘off’ state, independent of ligand binding. Taken together, these results suggest that DctBp's signal perception and transduction occur through a ‘ligand‐induced dimeric switch’, in which the changes in the dimeric conformations upon ligand binding are responsible for the signal transduction in DctB.     

Tryptophan fluorescence reveals conformational changes in the acetylcholine binding protein

The recent characterization of an acetylcholine binding protein (AChBP) from the fresh water snail, Lymnaea stagnalis, shows it to be a structural homolog of the extracellular domain of the nicotinic acetylcholine receptor (nAChR). To ascertain whether the AChBP exhibits the recognition properties and functional states of the nAChR, we have expressed the protein in milligram quantities from a synthetic cDNA transfected into human embryonic kidney (HEK) cells. The protein secreted into the medium shows a pentameric rosette structure with ligand stoichiometry approximating five sites per pentamer. Surprisingly, binding of acetylcholine, selective agonists, and antagonists ranging from small alkaloids to larger peptides results in substantial quenching of the intrinsic tryptophan fluorescence. Using stopped-flow techniques, we demonstrate rapid rates of association and dissociation of agonists and slow rates for the alpha-neurotoxins. Since agonist binding occurs in millisecond time frames, and the alpha-neurotoxins may induce a distinct conformational state for the AChBP-toxin complex, the snail protein shows many of the properties expected for receptor recognition of interacting ligands. Thus, the marked tryptophan quenching not only documents the importance of aromatic residues in ligand recognition, but establishes that the AChBP will be a useful functional as well as structural surrogate of the nicotinic receptor.     

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.     

Microsecond simulations of mdm2 and its complex with p53 yield insight into force field accuracy and conformational dynamics

The p53‐MDM2 complex is both a major target for cancer drug development and a valuable model system for computational predictions of protein‐ligand binding. To investigate the accuracy of molecular simulations of MDM2 and its complex with p53, we performed a number of long (200 ns to 1 µs) explicit‐solvent simulations using a range of force fields. We systematically compared nine popular force fields (AMBER ff03, ff12sb, ff14sb, ff99sb, ff99sb‐ildn, ff99sb‐ildn‐nmr, ff99sb‐ildn‐phi, CHARMM22*, and CHARMM36) against experimental chemical shift data, and found similarly accurate results, with microsecond simulations achieving better agreement compared to 200‐ns trajectories. Although the experimentally determined apo structure has a closed binding cleft, simulations in all force fields suggest the apo state of MDM2 is highly flexible, and able to sample holo‐like conformations, consistent with a conformational selection model. Initial structuring of the MDM2 lid region, known to competitively bind the binding cleft, is also observed in long simulations. Taken together, these results show molecular simulations can accurately sample conformations relevant for ligand binding. We expect this study to inform future computational work on folding and binding of MDM2 ligands. Proteins 2015; 83:1665–1676. © 2015 Wiley Periodicals, Inc.     

Influence of ligand binding on structure and thermostability of human α1‐acid glycoprotein

Ligand binding of neutral progesterone, basic propranolol, and acidic warfarin to human α₁‐acid glycoprotein (AGP) was investigated by Raman spectroscopy. The binding itself is characterized by a uniform conformational shift in which a tryptophan residue is involved. Slight differences corresponding to different contacts of the individual ligands inside the β‐barrel are described. Results are compared with in silico ligand docking into the available crystal structure of deglycosylated AGP using quantum/molecular mechanics. Calculated binding energies are ?18.2, ?14.5, and ?11.5 kcal/mol for warfarin, propranolol, and progesterone, respectively. These calculations are consistent with Raman difference spectroscopy; nevertheless, minor discrepancies in the precise positions of the ligands point to structural differences between deglycosylated and native AGP. Thermal dynamics of AGP with/without bounded warfarin was followed by Raman spectroscopy in a temperature range of 10–95 °C and analyzed by principal component analysis. With increasing temperature, a slight decrease of α‐helical content is observed that coincides with an increase in β‐sheet content. Above 45 °C, also β‐strands tend to unfold, and the observed decrease in β‐sheet coincides with an increase of β‐turns accompanied by a conformational shift of the nearby disulfide bridge from high‐energy trans‐gauche‐trans to more relaxed gauche‐gauche‐trans. This major rearrangement in the vicinity of the bridge is not only characterized by unfolding of the β‐sheet but also by subsequent ligand release. Hereby, ligand binding alters the protein dynamics, and the more rigid protein–ligand complex shows an improved thermal stability, a finding that contributes to the reported chaperone‐like function of AGP. Copyright © 2015 John Wiley & Sons, Ltd.     

Structural determinants for interaction of partial agonists with acetylcholine binding protein and neuronal α7 nicotinic acetylcholine receptor

The pentameric acetylcholine‐binding protein (AChBP) is a soluble surrogate of the ligand binding domain of nicotinic acetylcholine receptors. Agonists bind within a nest of aromatic side chains contributed by loops C and F on opposing faces of each subunit interface. Crystal structures of Aplysia AChBP bound with the agonist anabaseine, two partial agonists selectively activating the α7 receptor, 3‐(2,4‐dimethoxybenzylidene)‐anabaseine and its 4‐hydroxy metabolite, and an indole‐containing partial agonist, tropisetron, were solved at 2.7–1.75 Å resolution. All structures identify the Trp 147 carbonyl oxygen as the hydrogen bond acceptor for the agonist‐protonated nitrogen. In the partial agonist complexes, the benzylidene and indole substituent positions, dictated by tight interactions with loop F, preclude loop C from adopting the closed conformation seen for full agonists. Fluctuation in loop C position and duality in ligand binding orientations suggest molecular bases for partial agonism at full‐length receptors. This study, while pointing to loop F as a major determinant of receptor subtype selectivity, also identifies a new template region for designing α7‐selective partial agonists to treat cognitive deficits in mental and neurodegenerative disorders.     

Ligand binding to a high‐energy partially unfolded protein

The conformational energy landscape of a protein determines populations of all possible conformations of the protein and also determines the kinetics of the conversion between the conformations. Interaction with ligands influences the conformational energy landscapes of proteins and shifts populations of proteins in different conformational states. To investigate the effect of ligand binding on partial unfolding of a protein, we use Escherichia coli dihydrofolate reductase (DHFR) and its functional ligand NADP⁺ as a model system. We previously identified a partially unfolded form of DHFR that is populated under native conditions. In this report, we determined the free energy for partial unfolding of DHFR at varying concentrations of NADP⁺ and found that NADP⁺ binds to the partially unfolded form as well as the native form. DHFR unfolds partially without releasing the ligand, though the binding affinity for NADP⁺ is diminished upon partial unfolding. Based on known crystallographic structures of NADP⁺‐bound DHFR and the model of the partially unfolded protein we previously determined, we propose that the adenosine‐binding domain of DHFR remains folded in the partially unfolded form and interacts with the adenosine moiety of NADP⁺. Our result demonstrates that ligand binding may affect the conformational free energy of not only native forms but also high‐energy non‐native forms.     

Ligand-induced conformational changes in the acetylcholine-binding protein analyzed by hydrogen-deuterium exchange mass spectrometry

Recent x-ray crystallographic studies of the acetylcholine-binding protein (AChBP) suggest that loop C, found at the circumference of the pentameric molecule, shows distinctive conformational changes upon antagonist and agonist occupation. We have employed hydrogen-deuterium exchange mass spectrometry to examine the influence of bound ligands on solvent exposure of AChBP. Quantitative measurements of deuterium incorporation are possible for approximately 56% of the Lymnaea AChBP sequence, covering primarily the outer surface of AChBP. In the apoprotein, two regions flanking the ligand occupation site at the subunit interface, loop C (residues 175-193) and loop F (residues 164-171), show greater extents of solvent exchange than other regions of the protein including the N- and C-terminal regions. Occupation by nicotinic agonists, epibatidine and lobeline, and nicotinic antagonists, methyllycaconitine, alpha-bungarotoxin, and alpha-cobratoxin, markedly restricts the exchange of loop C amide protons, influencing both the rates and degrees of exchange. Solvent exposure of loop C and its protection by ligand suggest that in the apoprotein, loop C exhibits rapid fluctuations in an open conformation. Bound agonists restrict solvent exposure through loop closure, whereas the larger antagonists restrict solvent exposure largely through occlusion of solvent. Loop F, found on the complementary subunit surface at the interface, also reveals ligand selective changes in amide proton exchange rates. Agonists do not affect solvent accessibility of loop F, whereas certain antagonists cause subtle accessibility changes. These results reveal dynamic states and fluctuating movements in the vicinity of the binding site for unligated AChBP that can be influenced selectively by ligands.     

Structural properties of amyloid β(1‐40) dimer explored by replica exchange molecular dynamics simulations

Replica exchange molecular dynamics simulations (300 ns) were used to study the dimerization of amyloid β(1‐40) (Aβ(1‐40)) polypeptide. Configurational entropy calculations revealed that at physiological temperature (310 K, 37°C) dynamic dimers are formed by randomly docked monomers. Free energy of binding of the two chains to each other was ?93.56 ± 6.341 kJ mol^?1. Prevalence of random coil conformations was found for both chains with the exceptions of increased β‐sheet content from residues 16‐21 and 29‐32 of chain A and residues 15‐21 and 30‐33 of chain B with β‐turn/β‐bend conformations in both chains from residues 1‐16, 21‐29 of chain A, 1‐16, and 21‐29 of chain B. There is a mixed β‐turn/β‐sheet region from residues 33‐38 of both chains. Analysis of intra‐ and interchain residue distances shows that, although the individual chains are highly flexible, the dimer system stays in a loosely packed antiparallel β‐sheet configuration with contacts between residues 17‐21 of chain A with residues 17‐21 and 31‐36 of chain B as well as residues 31‐36 of chain A with residues 17‐21 and 31‐36 of chain B. Based on dihedral principal component analysis, the antiparallel β‐sheet‐loop‐β‐sheet conformational motif is favored for many low energy sampled conformations. Our results show that Aβ(1‐40) can form dynamic dimers in aqueous solution that have significant conformational flexibility and are stabilized by collapse of the central and C‐terminal hydrophobic cores with the expected β‐sheet‐loop‐β‐sheet conformational motif. Proteins 2017; 85:1024–1045. © 2017 Wiley Periodicals, Inc.     

Crystal Structure of Lymnaea stagnalis AChBP Complexed with the Potent nAChR Antagonist DHβE Suggests a Unique Mode of Antagonism

Nicotinic acetylcholine receptors (nAChRs) are pentameric ligand-gated ion channels that belong to the Cys-loop receptor superfamily. These receptors are allosteric proteins that exist in different conformational states, including resting (closed), activated (open), and desensitized (closed) states. The acetylcholine binding protein (AChBP) is a structural homologue of the extracellular ligand-binding domain of nAChRs. In previous studies, the degree of the C-loop radial extension of AChBP has been assigned to different conformational states of nAChRs. It has been suggested that a closed C-loop is preferred for the active conformation of nAChRs in complex with agonists whereas an open C-loop reflects an antagonist-bound (closed) state. In this work, we have determined the crystal structure of AChBP from the water snail Lymnaea stagnalis (Ls) in complex with dihydro-β-erythroidine (DHβE), which is a potent competitive antagonist of nAChRs. The structure reveals that binding of DHβE to AChBP imposes closure of the C-loop as agonists, but also a shift perpendicular to previously observed C-loop movements. These observations suggest that DHβE may antagonize the receptor via a different mechanism compared to prototypical antagonists and toxins.     

Conformational selection in protein binding and function

Protein binding and function often involves conformational changes. Advanced nuclear magnetic resonance (NMR) experiments indicate that these conformational changes can occur in the absence of ligand molecules (or with bound ligands), and that the ligands may “select” protein conformations for binding (or unbinding). In this review, we argue that this conformational selection requires transition times for ligand binding and unbinding that are small compared to the dwell times of proteins in different conformations, which is plausible for small ligand molecules. Such a separation of timescales leads to a decoupling and temporal ordering of binding/unbinding events and conformational changes. We propose that conformational‐selection and induced‐change processes (such as induced fit) are two sides of the same coin, because the temporal ordering is reversed in binding and unbinding direction. Conformational‐selection processes can be characterized by a conformational excitation that occurs prior to a binding or unbinding event, while induced‐change processes exhibit a characteristic conformational relaxation that occurs after a binding or unbinding event. We discuss how the ordering of events can be determined from relaxation rates and effective on‐ and off‐rates determined in mixing experiments, and from the conformational exchange rates measured in advanced NMR or single‐molecule fluorescence resonance energy transfer experiments. For larger ligand molecules such as peptides, conformational changes and binding events can be intricately coupled and exhibit aspects of conformational‐selection and induced‐change processes in both binding and unbinding direction.     

Galanthamine and non-competitive inhibitor binding to ACh-binding protein: evidence for a binding site on non-alpha-subunit interfaces of heteromeric neuronal nicotinic receptors

Rapid neurotransmission is mediated through a superfamily of Cys-loop receptors that includes the nicotinic acetylcholine (nAChR), gamma-aminobutyric acid (GABA(A)), serotonin (5-HT(3)) and glycine receptors. A class of ligands, including galanthamine, local anesthetics and certain toxins, interact with nAChRs non-competitively. Suggested modes of action include blockade of the ion channel, modulation from undefined extracellular sites, stabilization of desensitized states, and association with annular or boundary lipid. Alignment of mammalian Cys-loop receptors shows aromatic residues, found in the acetylcholine or ligand-binding pocket of nAChRs, are conserved in all subunit interfaces of neuronal nAChRs, including those that are not formed by alpha subunits on the principal side of the transmitter binding site. The amino-terminal domain containing the ligand recognition site is homologous to the soluble acetylcholine-binding protein (AChBP) from mollusks, an established structural and functional surrogate. We assess ligand specificity and employ X-ray crystallography with AChBP to demonstrate ligand interactions at subunit interfaces lacking vicinal cysteines (i.e. the non-alpha subunit interfaces in nAChRs). Non-competitive nicotinic ligands bind AChBP with high affinity (K(d) 0.015-6 microM). We mutated the vicinal cysteine residues in loop C of AChBP to mimic the non-alpha subunit interfaces of neuronal nAChRs and other Cys loop receptors. Classical nicotinic agonists show a 10-40-fold reduction in binding affinity, whereas binding of ligands known to be non-competitive are not affected. X-ray structures of cocaine and galanthamine bound to AChBP (1.8 A and 2.9 A resolution, respectively) reveal interactions deep within the subunit interface and the absence of a contact surface with the tip of loop C. Hence, in addition to channel blocking, non-competitive interactions with heteromeric neuronal nAChR appear to occur at the non-alpha subunit interface, a site presumed to be similar to that of modulating benzodiazepines on GABA(A) receptors.     

Ligand binding to the membrane-bound acetylcholine receptor from Torpedo marmorata: a complete mathematical analysis.

We have studied by means of equilibrium binding and kinetic experiments the interaction of the membrane-bound nicotinic acetylcholine receptor (nACHR) from Torpedo marmorata with [3H]acetylcholine and the fluorescent agonist NBD-5-acylcholine. In agreement with previous studies by others, we observed the preexistence, in the absence of ligand, of an equilibrium between two states of the nAChR, one with high affinity and the other with low affinity for agonist. As additional requirements for a minimal reaction scheme, we recognized (i) the existence of two ligand-binding sites, each of which may exist in two conformational states when occupied, and (ii) ligand-induced transitions between these conformations. Employing a special form of the allosteric model which considers these requirements, we then developed a suitable algorithm in order to simultaneously fit the whole set of equilibrium binding and kinetic data obtained for the two ligands. In this way we determined for a minimal model of the mechanism of action of the nAChR the complete set of rate constants and KD values involved. With these values available, we were able to simulate the rise and fall in the concentrations of individual receptor-ligand complexes and conformations occurring in the course of excitatory events at the electrocyte synapse. The membrane environment of the nAChR plays a decisive role with respect to the rates of conformational change of the nAChR occurring in the course of ligand interaction. Thus, artificial changes in membrane structure and composition can speed up by several orders of magnitude the rate of conformational change ("desensitization"). A proper structure of the surrounding membrane hence is a prerequisite for the physiological function of the membrane-embedded nAChR.     

Structural studies on 10-hydroxygeraniol dehydrogenase: A novel linear substrate-specific dehydrogenase from Catharanthus roseus

Conversion of 10-hydroxygeraniol to 10-oxogeranial is a crucial step in iridoid biosynthesis. This reaction is catalyzed by a zinc-dependent alcohol dehydrogenase, 10-hydroxygeraniol dehydrogenase, belonging to the family of medium-chain dehydrogenase/reductase (MDR). Here, we report the crystal structures of a novel 10-hydroxygeraniol dehydrogenase from Catharanthus roseus in its apo and nicotinamide adenine dinucleotide phosphate (NADP⁺) bound forms. Structural analysis and docking studies reveal how subtle conformational differences of loops L1, L2, L3, and helix α9' at the orifice of the catalytic site confer differential activity of the enzyme toward various substrates, by modulating the binding pocket shape and volume. The present study, first of its kind, provides insights into the structural basis of substrate specificity of MDRs specific to linear substrates. Furthermore, comparison of apo and NADP⁺ bound structures suggests that the enzyme adopts open and closed states to facilitate cofactor binding.     

Cofactor‐induced reversible folding of Flavodoxin‐4 from Lactobacillus acidophilus

Flavodoxins in combination with the flavin mononucleotide (FMN) cofactor play important roles for electron transport in prokaryotes. Here, novel insights into the FMN‐binding mechanism to flavodoxins‐4 were obtained from the NMR structures of the apo‐protein from Lactobacillus acidophilus (YP_193882.1) and comparison of its complex with FMN. Extensive reversible conformational changes were observed upon FMN binding and release. The NMR structure of the FMN complex is in agreement with the crystal structure (PDB ID: 3EDO ) and exhibits the characteristic flavodoxin fold, with a central five‐stranded parallel β–sheet and five α‐helices forming an α/β‐sandwich architecture. The structure differs from other flavoproteins in that helix α2 is oriented perpendicular to the β‐sheet and covers the FMN‐binding site. This helix reversibly unfolds upon removal of the FMN ligand, which represents a unique structural rearrangement among flavodoxins.     

Molecular modeling and NMR studies of benzyl substituted mannosyl trisaccharide binding to two mannose‐specific lectins: Allium sativam agglutinin I and Concanavalin A

The interaction of trimannoside, α?benzyl 3, 6‐di‐O‐(α‐D ‐mannopyranosyl)‐α‐D ‐mannopyranoside, 1 with ASAI (Allium sativam agglutinin I, garlic lectin) was studied to reveal the conformational preferences of this ligand in bound‐state and detailed binding mode at atomic level. The binding phenomenon was then compared with another well‐known mannose‐binding lectin, ConA (Concanavalin A). Structural studies of the ligand in free state were done using NMR spectroscopy and Molecular Dynamics simulations. It is found that the substituted‐trimannoside can undergo conformational transitions in solution, with one major and one minor conformation per glycosidic linkage (α 1→3 and α 1→6). On the other hand in the bound‐state only one of the two major conformations was significantly populated. The role of phenyl ring in the binding process was explored. An extended binding site was observed for the trimannoside in ASAI utilizing the aromatic substituent, which is not seen in ConA. Binding data from difference absorption spectroscopy supported this fact that the binding of benzyl‐substituted ligand is tighter with ASAI than ConA. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 952–967, 2010.     

Investigation of the influence of external factors on the conformational dynamics of rhodopsin-like receptors by means of molecular dynamics simulation

Abstract

The study reports about the influence of binding of orthosteric ligands on the conformational dynamics of β-2-adrenoreceptor. Using molecular dynamics (MD) simulation, we found that there was a little fraction of active states of the receptor in its apo (ligand-free) ensemble. Analysis of MD trajectories indicated that such spontaneous activation of the receptor is accompanied by the motion in intracellular part of its alpha-helices. Thus, receptor’s constitutive activity directly results from its conformational dynamics. On the other hand, the binding of a full agonist resulted in a significant shift of the initial equilibrium towards its active state. Finally, the binding of the inverse agonist stabilized the receptor in its inactive state. It is likely that the binding of inverse agonists might be a universal way of constitutive activity inhibition in vivo. Our results indicate that ligand binding redistribute pre-existing conformational degrees of freedom (in accordance to the Monod–Wyman–Changeux Model) of the receptor rather than cause induced fit in it. Therefore, the ensemble of biologically relevant receptor conformations is encoded in its spatial structure, and individual conformations from that ensemble might be used by the cell in conformity with the physiological behavior.     

Dual effects of familial Alzheimer's disease mutations (D7H,D7N,and H6R) on amyloid β peptide: Correlation dynamics and zinc binding

Although the N‐terminal region of Amyloid β (Aβ) peptides plays dual roles as metal‐coordinating sites and conformational modulator, few studies have been performed to explore the effects of mutations at this region on the overall conformational ensemble of Aβ and the binding propensity of metal ions. In this work, we focus on how three familial Alzheimer's disease mutations (D7H, D7N, and H6R) alter the structural characteristics and thermodynamic stabilities of Aβ42 using molecular dynamics simulations. We observe that each mutation displays increased β‐sheet structures in both N and C termini. In particular, both the N terminus and central hydrophobic region of D7H can form stable β‐hairpin structures with its C terminus. The conserved turn structure at Val²⁴–Lys²⁸ in all peptides and Zn²⁺‐bound Aβ42 is confirmed as the common structural motif to nucleate folding of Aβ. Each mutant can significantly increase the solvation free energy and thus enhance the aggregation of Aβ monomers. The correlation dynamics between Aβ(1–16) and Aβ(17–42) fragments are elucidated by linking the domain motions with the corresponding structured conformations. We characterize the different populations of correlated domain motions for each mutant from a more macroscopic perspective, and unexpectedly find that Zn²⁺‐bound Aβ42 ensemble shares the same populations as Aβ42, indicating that the binding of Zn²⁺ to Aβ follows the conformational selection mechanism, and thus is independent of domain motions, even though the structures of Aβ have been modified at a residue level. Proteins 2014; 82:3286–3297. © 2014 Wiley Periodicals, Inc.