Ligand-induced conformational changes in a thermophilic ribose-binding protein

Background

Members of the periplasmic binding protein (PBP) superfamily are involved in transport and signaling processes in both prokaryotes and eukaryotes. Biological responses are typically mediated by ligand-induced conformational changes in which the binding event is coupled to a hinge-bending motion that brings together two domains in a closed form. In all PBP-mediated biological processes, downstream partners recognize the closed form of the protein. This motion has also been exploited in protein engineering experiments to construct biosensors that transduce ligand binding to a variety of physical signals. Understanding the mechanistic details of PBP conformational changes, both global (hinge bending, twisting, shear movements) and local (rotamer changes, backbone motion), therefore is not only important for understanding their biological function but also for protein engineering experiments.

Results

Here we present biochemical characterization and crystal structure determination of the periplasmic ribose-binding protein (RBP) from the hyperthermophile Thermotoga maritima in its ribose-bound and unliganded state. The T. maritima RBP (tmRBP) has 39% sequence identity and is considerably more resistant to thermal denaturation ( ^app T _m value is 108°C) than the mesophilic Escherichia coli homolog (ecRBP) ( ^app T _m value is 56°C). Polar ligand interactions and ligand-induced global conformational changes are conserved among ecRBP and tmRBP; however local structural rearrangements involving side-chain motions in the ligand-binding site are not conserved.

Conclusion

Although the large-scale ligand-induced changes are mediated through similar regions, and are produced by similar backbone movements in tmRBP and ecRBP, the small-scale ligand-induced structural rearrangements differentiate the mesophile and thermophile. This suggests there are mechanistic differences in the manner by which these two proteins bind their ligands and are an example of how two structurally similar proteins utilize different mechanisms to form a ligand-bound state.     

Ligand-induced conformational changes within a hexameric Acyl-CoA thioesterase

Acyl-coenzyme A (acyl-CoA) thioesterases play a crucial role in the metabolism of activated fatty acids, coenzyme A, and other metabolic precursor molecules including arachidonic acid and palmitic acid. These enzymes hydrolyze coenzyme A from acyl-CoA esters to mediate a range of cellular functions including β-oxidation, lipid biosynthesis, and signal transduction. Here, we present the crystal structure of a hexameric hot-dog domain-containing acyl-CoA thioesterase from Bacillus halodurans in the apo-form and provide structural and comparative analyses to the coenzyme A-bound form to identify key conformational changes induced upon ligand binding. We observed dramatic ligand-induced changes at both the hot-dog dimer and the trimer-of-dimer interfaces; the dimer interfaces in the apo-structure differ by over 20% and decrease to about half the size in the ligand-bound state. We also assessed the specificity of the enzyme against a range of fatty acyl-CoA substrates and have identified a preference for short-chain fatty acyl-CoAs. Coenzyme A was shown both to negatively regulate enzyme activity, representing a direct inhibitory feedback, and consistent with the structural data, to destabilize the quaternary structure of the enzyme. Coenzyme A-induced conformational changes in the C-terminal helices of enzyme were assessed through mutational analysis and shown to play a role in regulating enzyme activity. The conformational changes are likely to be conserved from bacteria through to humans and provide a greater understanding, particularly at a structural level, of thioesterase function and regulation.     

Ligand-induced conformational and structural dynamics changes in Escherichia coli cyclic AMP receptor protein

Cyclic AMP receptor protein (CRP) regulates the expression of a large number of genes in E. coli. It is activated by cAMP binding, which leads to some yet undefined conformational changes. These changes do not involve significant redistribution of secondary structures. A potential mechanism of activation is a ligand-induced change in structural dynamics. Hence, the cAMP-mediated conformational and structural dynamics changes in the wild-type CRP were investigated using hydrogen-deuterium exchange and Fourier transform infrared spectroscopy. Upon cAMP binding, the two functional domains within the wild-type CRP undergo conformational and structural dynamics changes in two opposite directions. While the smaller DNA-binding domain becomes more flexible, the larger cAMP-binding domain shifts to a less dynamic conformation, evidenced by a faster and a slower amide H-D exchange, respectively. To a lesser extent, binding of cGMP, a nonfunctional analogue of cAMP, also stabilizes the cAMP-binding domain, but it fails to mimic the relaxation effect of cAMP on the DNA-binding domain. Despite changes in the conformation and structural dynamics, cAMP binding does not alter significantly the secondary structural composition of the wild-type CRP. The apparent difference between functional and nonfunctional analogues of cAMP is the ability of cAMP to effect an increase in the dynamic motions of the DNA binding domain.     

Ligand-induced conformational changes in cytosolic protein kinase C

The changes in intrinsic spectral properties of protein kinase C were monitored upon association with its divalent cation and lipid activators in a model membrane system. The enzyme demonstrated changes in both its intrinsic fluorescence and far ultraviolet circular dichroism spectra upon association with lipid vesicles in the absence of calcium. The acidic phospholipid, phosphatidylserine, significantly quenched the intrinsic tryptophan fluorescence and was also the most potent lipid support for the phosphorylating activity of the enzyme. The enzyme was fully activated by a number of Ca2(+)-lipid combinations which correlated with maximal fluorescence quenching (40-50%) of available tryptophan residues in hydrophobic domains. The circular dichroism structure of the associated active-protein Ca2(+)-lipid complexes suggested different active enzyme secondary structures. However, the Ca2(+)-dependent changes in fluorescence and circular dichroism spectra were observed only after the enzyme associated with the lipid vesicles. These data suggest that protein kinase C has the properties of a complex multidomain protein and provides an additional perspective into the mechanism of protein kinase C activation.     

Ligand-induced changes in dynamics in the RT loop of the C-terminal SH3 domain of Sem-5 indicate cooperative conformational coupling

We report the effects of peptide binding on the (15)N relaxation rates and chemical shifts of the C-SH3 of Sem-5. (15)N spin-lattice relaxation time (T(1)), spin-spin relaxation time (T(2)), and ((1)H)-(15)N NOE were obtained from heteronuclear 2D NMR experiments. These parameters were then analyzed using the Lipari-Szabo model free formalism to obtain parameters that describe the internal motions of the protein. High-order parameters (S(2) > 0.8) are found in elements of regular secondary structure, whereas some residues in the loop regions show relatively low-order parameters, notably the RT loop. Peptide binding is characterized by a significant decrease in the (15)N relaxation in the RT loop. Concomitant with the change in dynamics is a cooperative change in chemical shifts. The agreement between the binding constants calculated from chemical shift differences and that obtained from ITC indicates that the binding of Sem-5 C-SH3 to its putative peptide ligand is coupled to a cooperative conformational change in which a portion of the binding site undergoes a significant reduction in conformational heterogeneity.     

Ligand-induced conformational changes and a mechanism for domain closure in Aspergillus nidulans dehydroquinate synthase

In order to investigate systematically substrate and cofactor-induced conformational changes in the enzyme dehydroquinate synthase (DHQS), eight structures representing a series of differently liganded states have been determined in a total of six crystal forms. DHQS in the absence of the substrate analogue carbaphosphonate, either unliganded or in the presence of NAD or ADP, is in an open form where a relative rotation of 11-13 degrees between N and C-terminal domains occurs.Analysis of torsion angle difference plots between sets of structures reveals eight rearrangements that appear relevant to domain closure and a further six related to crystal packing. Overlapping 21 different copies of the individual N and C-terminal DHQS domains further reveals a series of pivot points about which these movements occur and illustrates the way in which widely separated secondary structure elements are mechanically inter-linked to form "composite elements", which propagate structural changes across large distances.This analysis has provided insight into the basis of DHQS ligand-initiated domain closure and gives rise to the proposal of an ordered sequence of events involving substrate binding, and local rearrangements within the active site that are propagated to the hinge regions, leading to closure of the active-site cleft.     

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.     

Ligand-induced changes in the structure and dynamics of a human class Mu glutathione S-transferase

Glutathione transferases are detoxification enzymes that catalyze the addition of glutathione (GSH) to a wide variety of hydrophobic compounds. Although this group of enzymes has been extensively characterized by crystallographic studies, little is known about their dynamic properties. This study investigates the role of protein dynamics in the mechanism of a human class mu enzyme (GSTM2-2) by characterizing the motional properties of the unliganded enzyme, the enzyme-substrate (GSH) complex, an enzyme-product complex [S-(2,4-dinitrobenzyl)glutathione, GSDNB], and an enzyme-inhibitor complex (S-1-hexylglutathione, GSHEX). The kinetic on- and off-rates for these ligands are 10-20-fold lower than the diffusion limit, suggesting dynamic conformational heterogeneity of the active site. The off-rate of GSDNB is similar to the turnover number for its enzymatic formation, suggesting that product release is rate-limiting when 1-chloro-2,4-dinitrobenzene is the substrate. The dynamic properties of GSTM2-2 were investigated over a wide range of time scales using (15)N nuclear spin relaxation, residual dipolar couplings, and amide hydrogen-deuterium exchange rates. These data show that the majority of the protein backbone is rigid on the nanosecond to picosecond time scale for all forms of the enzyme. The presence of motion on the millisecond to microsecond time scale was detected for a small number of residues within the active site. These motions are likely to play a role in facilitating substrate binding and product release. The residual dipolar couplings also show that the conformation of the active site region is more open in solution than in the crystalline environment, further enhancing ligand accessibility to the active site. Amide hydrogen-deuterium exchange rates indicate a reduction in the dynamic properties of several residues near the active site due to the binding of ligand. GSH binding reduces the exchange rate of a number of residues in proximity to its binding site, while GSHEX causes a reduction in amide-exchange rates throughout the entire active site region. The location of the dinitrobenzene (DNB) ring in the GSDNB-GSTM2-2 complex was modeled using chemical shift changes that occur when GSDNB binds to the enzyme. The DNB ring makes a number of contacts with hydrophobic residues in the active site, including Met108. Replacement of Met108 with Ala increases the turnover number of the enzyme by a factor of 1.7.     

NMR detection of slow conformational dynamics in an endonuclease toxin

The cytotoxic activity of the secreted bacterial toxin colicin E9 is due to a non-specific DNase housed in the C-terminus of the protein. Double-resonance and triple-resonance NMR studies of the 134-amino acid¹⁵ N- and ¹³C/¹⁵N-labelled DNase domain are presented. Extensive conformational heterogeneity was evident from the presence of far more resonances than expected based on the amino acid sequence of the DNase, and from the appearance of chemical exchange cross-peaks in TOCSY and NOESY spectra. EXSY spectra were recorded to confirm that slow chemical exchange was occurring. Unambiguous sequence-specific resonance assignments are presented for one region of the protein, Pro⁶⁵-Asn⁷², which exists in two slowly exchanging conformers based on the identification of chemical exchange cross-peaks in 3D ¹H-¹H-¹⁵N EXSY-HSQC, NOESY-HSQC and TOCSY-HSQC spectra, together with C and C chemical shifts measured in triple-resonance spectra and sequential NH NOEs. The rates of conformational exchange for backbone amide resonances in this stretch of amino acids, and for the indole NH of either Trp²² or Trp⁵⁸, were determined from the intensity variation of the appropriate diagonal and chemical exchange cross-peaks recorded in 3D¹ H-¹H-¹⁵N NOESY-HSQC spectra. The data fitted a model in which this region of the DNase has two conformers, N_A and N_B, which interchange at 15 °C with a forward rate constant of 1.61 ± 0.5 s^-1 and a backward rate constant of 1.05 ± 0.5 s^-1. Demonstration of this conformational equilibrium has led to a reappraisal of a previously proposed kinetic scheme describing the interaction of E9 DNase with immunity proteins [Wallis et al. (1995) Biochemistry, 34, 13743–13750 and 13751–13759]. The revised scheme is consistent with the specific inhibitor protein for the E9 DNase, Im9, associating with both the N_A and N_B conformers of the DNase and with binding only to the N_B conformer detected because the rate of dissociation of the complex of Im9 and the N_A conformer, N_AI, is extremely rapid. In this model stoichiometric amounts of Im9 convert, the E9 DNase is converted wholly into the N_BI form. The possibility that cis–trans isomerisation of peptide bonds preceding proline residues is the cause of the conformational heterogeneity is discussed. E9 DNase contains 10 prolines, with two bracketing the stretch of amino acids that have allowed the N_A N_B interconversion to be identified, Pro⁶⁵ and Pro⁷³. The model assumes that one or both of these can exist in either the cis or trans form with strong Im9 binding possible to only one form.     

Ligand-induced conformational changes in the capping subdomain of a bacterial old yellow enzyme homologue and conserved sequence fingerprints provide new insights into substrate binding

We have recently reported that Shewanella oneidensis, a Gram-negative gamma-proteobacterium with a rich arsenal of redox proteins, possesses four old yellow enzyme (OYE) homologues. Here, we report a series of high resolution crystal structures for one of these OYEs, Shewanella yellow enzyme 1 (SYE1), in its oxidized form at 1.4A resolution, which binds a molecule of PEG 400 in the active site, and in its NADH-reduced and p-hydroxybenzaldehyde- and p-hydroxyacetophenone-bound forms at 1.7A resolution. Although the overall structure of SYE1 reveals a monomeric enzyme based on the alpha(8)beta(8) barrel scaffold observed for other OYEs, the active site exhibits a unique combination of features: a strongly butterfly-bent FMN cofactor both in the oxidized and NADH-reduced forms, a collapsed and narrow active site tunnel, and a novel combination of conserved residues involved in the binding of phenolic ligands. Furthermore, we identify a second p-hydroxybenzaldehyde-binding site in a hydrophobic cleft next to the entry of the active site tunnel in the capping subdomain, formed by a restructuring of Loop 3 to an "open" conformation. This constitutes the first evidence to date for the entire family of OYEs that Loop 3 may indeed play a dynamic role in ligand binding and thus provides insights into the elusive NADH complex and into substrate binding in general. Structure-based sequence alignments indicate that the novelties we observe in SYE1 are supported by conserved residues in a number of structurally uncharacterized OYEs from the beta- and gamma-proteobacteria, suggesting that SYE1 represents a new subfamily of bacterial OYEs.     

Ligand-induced conformational changes in lactose repressor: a phosphorescence and ODMR study of single-tryptophan mutants.

Phosphorescence and optically detected magnetic resonance (ODMR) measurements are reported on four single-tryptophan mutants of lac repressor protein from Escherichia coli: H74W/Wless, W201Y, Y273W/Wless, and F293W/Wless, where Wless represents a protein background containing the double mutation W201Y/W220Y. The single-tryptophan residues are located in the protein core region, either in the monomer-monomer interface of the tetrameric protein or in the region of the inducer binding cleft. Inducer binding elicits large changes in the energy (0,0-band wavelength shifts) and zero-field splitting energies (ZFS) of the triplet states for each of the mutant proteins except W201Y which exhibits more modest effects. F293W/Wless exists in two distinguishable conformations, only one of which appears to be sensitive to the presence of inducer. These effects of inducer binding can be attributed to a conformational change that alters specific polar interactions that occur at each affected tryptophan site. Changes in the tryptophan triplet state indicator depend on the existence of specific polar interactions that are altered by local atomic relocations.     

Ligand-induced conformational changes: improved predictions of ligand binding conformations and affinities

A computational docking strategy using multiple conformations of the target protein is discussed and evaluated. A series of low molecular weight, competitive, nonpeptide protein tyrosine phosphatase inhibitors are considered for which the x-ray crystallographic structures in complex with protein tyrosine phosphatase 1B (PTP1B) are known. To obtain a quantitative measure of the impact of conformational changes induced by the inhibitors, these were docked to the active site region of various structures of PTP1B using the docking program FlexX. Firstly, the inhibitors were docked to a PTP1B crystal structure cocrystallized with a hexapeptide. The estimated binding energies for various docking modes as well as the RMS differences between the docked compounds and the crystallographic structure were calculated. In this scenario the estimated binding energies were not predictive inasmuch as docking modes with low estimated binding energies corresponded to relatively large RMS differences when aligned with the corresponding crystal structure. Secondly, the inhibitors were docked to their parent protein structures in which they were cocrystallized. In this case, there was a good correlation between low predicted binding energy and a correct docking mode. Thirdly, to improve the predictability of the docking procedure in the general case, where only a single target protein structure is known, we evaluate an approach which takes possible protein side-chain conformational changes into account. Here, side chains exposed to the active site were considered in their allowed rotamer conformations and protein models containing all possible combinations of side-chain rotamers were generated. To evaluate which of these modeled active sites is the most likely binding site conformation for a certain inhibitor, the inhibitors were docked against all active site models. The receptor rotamer model corresponding to the lowest estimated binding energy is taken as the top candidate. Using this protocol, correct inhibitor binding modes could successfully be discriminated from proposed incorrect binding modes. Moreover, the ranking of the estimated ligand binding energies was in good agreement with experimentally observed binding affinities.     

Ligand-induced conformational change in the minimized insulin receptor

Within the class of insulin and insulin-like growth factor receptors, detailed information about the molecular recognition event at the hormone-receptor interface is limited by the absence of suitable co-crystals. We describe the use of a biologically active insulin derivative labeled with the NBD fluorophore (B29NBD-insulin) to characterize the mechanism of reversible 1:1 complex formation with a fragment of the insulin receptor ectodomain. The accompanying 40 % increase in the fluorescence quantum yield of the label provides the basis for a dynamic study of the hormone-receptor binding event. Stopped-flow fluorescence experiments show that the kinetics of complex formation are biphasic comprising a bimolecular binding event followed by a conformational change. Displacement with excess unlabeled insulin gave monophasic kinetics of dissociation. The rate data are rationalized in terms of available experiments on mutant receptors and the X-ray structure of a non-binding fragment of the receptor of the homologous insulin-like growth factor (IGF-1).     

Ligand-induced conformational changes in the Bacillus subtilis chemoreceptor McpB determined by disulfide crosslinking in vivo

Previously, we characterized the organization of the transmembrane (TM) domain of the Bacillus subtilis chemoreceptor McpB using disulfide crosslinking. Cysteine residues were engineered into serial positions along the two helices through the membrane, TM1 and TM2, as well as double mutants in TM1 and TM2, and the extent of crosslinking determined to characterize the organization of the TM domain. In this study, the organization of the TM domain was studied in the presence and absence of ligand to address what ligand-induced structural changes occur. We found that asparagine caused changes in crosslinking rate on all residues along the TM1-TM1' helical interface, whereas the crosslinking rate for almost all residues along the TM2-TM2' interface did not change. These results indicated that helix TM1 rotated counterclockwise and that TM2 did not move in respect to TM2' in the dimer on binding asparagine. Interestingly, intramolecular crosslinking of paired substitutions in 34/280 and 38/273 were unaffected by asparagine, demonstrating that attractant binding to McpB did not induce a "piston-like" vertical displacement of TM2 as seen for Trg and Tar in Escherichia coli. However, these paired substitutions produced oligomeric forms of receptor in response to ligand. This must be due to a shift of the interface between different receptor dimers, within previously suggested trimers of dimers, or even higher order complexes. Furthermore, the extent of disulfide bond formation in the presence of asparagine was unaffected by the presence of the methyl-modification enzymes, CheB and CheR, or the coupling proteins, CheW and CheV, demonstrating that these proteins must have local structural effects on the cytoplasmic domain that is not translated to the entire receptor. Finally, disulfide bond formation was also unaffected by binding proline to McpC. We conclude that ligand-binding induced a conformational change in the TM domain of McpB dimers as an excitation signal that is likely propagated within the cytoplasmic region of receptors and that subsequent adaptational events do not affect this new TM domain conformation.     

Ligand-induced conformational changes via flexible linkers in the amino-terminal region of the inositol 1,4,5-trisphosphate receptor

Cytoplasmic Ca2+ signals are highly regulated by various ion transporters, including the inositol 1,4,5-trisphosphate (IP(3)) receptor (IP(3)R), which functions as a Ca2+ release channel on the endoplasmic reticulum membrane. Crystal structures of the two N-terminal regulatory regions from type 1 IP(3)R have been reported; those of the IP(3)-binding core (IP(3)R(CORE)) with bound IP(3), and the suppressor domain. This study examines the structural effects of ligand binding on an IP(3)R construct, designated IP(3)R(N), that contains both the IP(3)-binding core and the suppressor domain. Our circular dichroism results reveal that the IP(3)-bound and IP(3)-free states have similar secondary structure content, consistent with preservation of the overall fold within the individual domains. Thermal denaturation data show that, while IP(3) has a large effect on the stability of IP(3)R(CORE), it has little effect on IP(3)R(N), indicating that the suppressor domain is critical to the stability of IP(3)R(N). The NMR data for IP(3)R(N) provide evidence for chemical exchange, which may be due to protein conformational dynamics in both apo and IP(3)-bound states: a conclusion supported by the small-angle X-ray scattering data. Further, the scattering data show that IP(3)R(N) undergoes a change in average conformation in response to IP(3) binding and the presence of Ca2+ in the solution. Taken together, these data lead us to propose that there are two flexible linkers in the N-terminal region of IP(3)R that join stably folded domains and give rise to an equilibrium mixture of conformational sub-states containing compact and more extended structures. IP(3) binding drives the conformational equilibrium toward more compact structures, while the presence of Ca2+ drives it to a more extended set.     

Ligand-induced conformational changes in tissue transglutaminase: Monte Carlo analysis of small-angle scattering data

Small-angle neutron and x-ray scattering experiments have been performed on type 2 tissular transglutaminase to characterize the conformational changes that bring about Ca(2+) activation and guanosine triphosphate (GTP) inhibition. The native and a proteolyzed form of the enzyme, in the presence and in the absence of the two effectors, were considered. To describe the shape of transglutaminase in the different conformations, a Monte Carlo method for calculating small-angle neutron scattering profiles was developed by taking into account the computer-designed structure of the native transglutaminase, the results of the Guinier analysis, and the essential role played by the solvent-exposed peptide loop for the conformational changes of the protein after activation. Although the range of the neutron scattering data is rather limited, by using the Monte Carlo analysis, and because the structure of the native protein is available, the distribution of the protein conformations after ligand interaction was obtained. Calcium activation promotes a rotation of the C-terminal with respect to the N-terminal domain around the solvent-exposed peptide loop that connects the two regions. The psi angle between the longest axes of the two pairs of domains is found to be above 50 degrees, larger than the psi value of 35 degrees calculated for the native transglutaminase. On the other hand, the addition of GTP makes possible conformations characterized by psi angles lower than 34 degrees. These results are in good agreement with the proposed enzyme activity regulation: in the presence of GTP, the catalytic site is shielded by the more compact protein structure, while the conformational changes induced by Ca(2+) make the active site accessible to the substrate.