Determinants of cooperativity and site selectivity in human ileal bile acid binding protein

Human ileal bile acid binding protein (I-BABP) is a member of the family of intracellular lipid-binding proteins and is thought to play a role in the enterohepatic circulation of bile salts. Our group has previously shown that human I-BABP binds two molecules of glycocholate (GCA) with low intrinsic affinity but an extraordinary high degree of positive cooperativity. Besides the strong positive cooperativity, human I-BABP exhibits a high degree of site selectivity in its interactions with GCA and glycochenodeoxycholate (GCDA), the two major bile salts in humans. In this study, on the basis of our first generation nuclear magnetic resonance (NMR) structure of the ternary complex of human I-BABP with GCA and GCDA, we introduced single-residue mutations at certain key positions in the binding pocket that might disrupt a hydrogen-bonding network, a likely way of energetic communication between the two sites. Macroscopic binding parameters were determined using isothermal titration calorimetry, and site selectivity was monitored by NMR spectroscopy of isotopically enriched bile salts. According to our results, cooperativity and site selectivity are not linked in human I-BABP. While cooperativity is governed by a subtle interplay of entropic and enthalpic contributions, site selectivity appears to be determined by more localized enthalpic effects. Possible communication pathways between the two binding sites are discussed.     

Kinetic mechanism of ligand binding in human ileal bile acid binding protein as determined by stopped-flow fluorescence analysis

Cooperative ligand binding to human ileal bile acid binding protein (I-BABP) was studied using the stopped-flow fluorescence technique. The kinetic data obtained for wild-type protein are in agreement with a four-step mechanism where after a fast conformational change on the millisecond time scale, the ligands bind in a sequential manner, followed by another, slow conformational change on the time scale of seconds. This last step is more pronounced in the case of glycocholate (GCA), the bile salt that binds with high positive cooperativity and is absent in mutant I-BABP proteins that lack positive cooperativity in their bile salt binding. These results suggest that positive cooperativity in human I-BABP is related to a slow conformational change of the protein, which occurs after the second binding step. Analogous to that in the intestinal fatty acid binding protein (I-FABP), we hypothesize that ligand binding in I-BABP is linked to a disorder-order transition between an open and a closed form of the protein.     

Conformational and dynamics changes induced by bile acids binding to chicken liver bile acid binding protein

The correlation between protein motions and function is a central problem in protein science. Several studies have demonstrated that ligand binding and protein dynamics are strongly correlated in intracellular lipid binding proteins (iLBPs), in which the high degree of flexibility, principally occurring at the level of helix-II, CD, and EF loops (the so-called portal area), is significantly reduced upon ligand binding. We have recently investigated by NMR the dynamic properties of a member of the iLBP family, chicken liver bile acid binding protein (cL-BABP), in its apo and holo form, as a complex with two bile salts molecules. Binding was found to be regulated by a dynamic process and a conformational rearrangement was associated with this event. We report here the results of molecular dynamics (MD) simulations performed on apo and holo cL-BABP with the aim of further characterizing the protein regions involved in motion propagation and of evaluating the main molecular interactions stabilizing bound ligands. Upon binding, the root mean square fluctuation values substantially decrease for CD and EF loops while increase for the helix-loop-helix region, thus indicating that the portal area is the region mostly affected by complex formation. These results nicely correlate with backbone dynamics data derived from NMR experiments. Essential dynamics analysis of the MD trajectories indicates that the major concerted motions involve the three contiguous structural elements of the portal area, which however are dynamically coupled in different ways whether in the presence or in the absence of the ligands. Motions of the EF loop and of the helical region are part of the essential space of both apo and holo-BABP and sample a much wider conformational space in the apo form. Together with NMR results, these data support the view that, in the apo protein, the flexible EF loop visits many conformational states including those typical of the holo state and that the ligand acts stabilizing one of these pre-existing conformations. The present results, in agreement with data reported for other iLBPs, sharpen our knowledge on the binding mechanism for this protein family.     

Steroid ring hydroxylation patterns govern cooperativity in human bile acid binding protein

Human ileal bile acid binding protein (I-BABP) is a member of the intracellular lipid binding protein family. This protein is thought to function in the transcellular transport and enterohepatic circulation of bile salts. Human I-BABP binds two molecules of glycocholate, the physiologically most abundant bile salt, with modest intrinsic affinity but a remarkably high degree of positive cooperativity. Here we report a calorimetric analysis for the binding of a broad panel of bile salts to human I-BABP. The interaction of I-BABP with nine physiologically relevant derivatives of cholic acid, chenodeoxycholic acid, and deoxycholic acid in their conjugated (glycine and taurine) and unconjugated forms was monitored by isothermal titration calorimetry. All bile salts bound to I-BABP with a 2:1 stoichiometry and similar overall affinity, but the derivatives of cholic acid displayed much higher Hill coefficients, a measure of macroscopic positive cooperativity. To test whether the cooperativity was dependent on individual structural features of the bile salt side chain, a series of side-chain-extended bile salts that lacked a hydrogen bond donor or acceptor at C-24 were chemically synthesized. These synthetic variants exhibited the same energetic and cooperativity profile as the naturally occurring bile salts. Our findings indicate that cooperativity in bile salt-I-BABP recognition is governed by the pattern of steroid B- and C-ring hydroxylation and not the presence or type of side-chain conjugation.     

Structural requirements for cooperativity in ileal bile acid-binding proteins

Ileal bile acid-binding proteins (I-BABP), belonging to the family of intracellular lipid-binding proteins, control bile acid trafficking in enterocytes and participate in regulating the homeostasis of these cholesterol-derived metabolites. I-BABP orthologues share the same structural fold and are able to host up to two ligands in their large internal cavities. However variations in the primary sequences determine differences in binding properties such as the degree of binding cooperativity. To investigate the molecular requirements for cooperativity we adopted a gain-of-function approach, exploring the possibility to turn the noncooperative chicken I-BABP (cI-BABP) into a cooperative mutant protein. To this aim we first solved the solution structure of cI-BABP in complex with two molecules of the physiological ligand glycochenodeoxycholate. A comparative structural analysis with closely related members of the same protein family provided the basis to design a double mutant (H99Q/A101S cI-BABP) capable of establishing a cooperative binding mechanism. Molecular dynamics simulation studies of the wild type and mutant complexes and essential dynamics analysis of the trajectories supported the role of the identified amino acid residues as hot spot mediators of communication between binding sites. The emerging picture is consistent with a binding mechanism that can be described as an extended conformational selection model.     

Mechanisms of allosteric gene regulation by NMR quantification of microsecond-millisecond protein dynamics

The trp RNA-binding attenuation protein (TRAP) is a paradigmatic allosteric protein that regulates the tryptophan biosynthetic genes associated with the trp operon in bacilli. The ring-shaped 11-mer TRAP is activated for recognition of a specific trp-mRNA target by binding up to 11 tryptophan molecules. To characterize the mechanisms of tryptophan-induced TRAP activation, we have performed methyl relaxation dispersion (MRD) nuclear magnetic resonance (NMR) experiments that probe the time-dependent structure of TRAP in the microsecond-to-millisecond "chemical exchange" time window. We find significant side chain flexibility localized to the RNA and tryptophan binding sites of the apo protein and that these dynamics are dramatically reduced upon ligand binding. Analysis of the MRD NMR data provides insights into the structural nature of transiently populated conformations sampled in solution by apo TRAP. The MRD data are inconsistent with global two-state exchange, indicating that conformational sampling in apo TRAP is asynchronous. These findings imply a temporally heterogeneous population of structures that are incompatible with RNA binding and substantiate the study of TRAP as a paradigm for probing and understanding essential dynamics in allosteric, regulatory proteins.     

Guanine nucleotides differentially modulate backbone dynamics of the STAS domain of the SulP/SLC26 transport protein Rv1739c of Mycobacterium tuberculosis

Enzymatic catalysis and protein signaling are dynamic processes that involve local and/or global conformational changes occurring across a broad range of time scales. (1) H-(15) N relaxation NMR provides a comprehensive understanding of protein backbone dynamics both in the apo (unliganded) and ligand-bound conformations, enabling both fast and slow internal motions of individual amino acid residues to be observed. We recently reported the structure and nucleotide binding properties of the sulfate transporter and anti-sigma factor antagonist (STAS) domain of Rv1739c, a SulP anion transporter protein of Mycobacterium tuberculosis. In the present study, we report (1) H-(15) N NMR backbone dynamics measurements [longitudinal (T(1) ), transverse (T(2) ) and steady-state ({(1) H}-(15) N) heteronuclear NOE] of the Rv1739c STAS domain, in the absence and presence of saturating concentrations of GTP and GDP. Analysis of measured relaxation data and estimated dynamic parameters indicated distinct features differentiating the binding of GTP and GDP to Rv1739c STAS. The 9.55 ns overall rotational correlation time of Rv1739c STAS increased to 10.48 ns in the presence of GTP, and to 13.25 ns in the presence of GDP, indicating significant nucleotide-induced conformational changes. These conformational changes were accompanied by slow time scale (μs to ms) motions in discrete regions of the protein, as reflected by guanine nucleotide-induced changes in relaxation parameters. The observed nucleotide-specific alterations in the relaxation properties of individual STAS residues reflect an increased molecular anisotropy and/or the emergence of conformational equilibria governing functional properties of the STAS domain.     

Dynamics of cellular retinoic acid binding protein I on multiple time scales with implications for ligand binding

Cellular retinoic acid binding protein I (CRABPI) belongs to the family of intracellular lipid binding proteins (iLBPs), all of which bind a hydrophobic ligand within an internal cavity. The structures of several iLBPs reveal minimal structural differences between the apo (ligand-free) and holo (ligand-bound) forms, suggesting that dynamics must play an important role in the ligand recognition and binding processes. Here, a variety of nuclear magnetic resonance (NMR) spectroscopy methods were used to systematically study the dynamics of both apo and holo CRABPI at various time scales. Translational and rotational diffusion constant measurements were used to study the overall motions of the proteins. Both apo and holo forms of CRABPI tend to self-associate at high (1.2 mM) concentrations, while at low concentrations (0.2 mM), they are predominantly monomeric. Rapid amide exchange rate and laboratory frame relaxation rate measurements at two spectrometer field strengths (500 and 600 MHz) were used to probe the internal motions of the individual residues. Several residues in the apo form, notably within the ligand recognition region, exhibit millisecond time scale motions that are significantly arrested in the holo form. In contrast, no significant differences in the high-frequency motions were observed between the two forms. These results provide direct experimental evidence for dynamics-induced ligand recognition and binding at a specifically defined time scale. They also exemplify the importance of dynamics in providing a more comprehensive understanding of how a protein functions.     

The X-ray structure of zebrafish (Danio rerio) ileal bile acid-binding protein reveals the presence of binding sites on the surface of the protein molecule

The ileal bile acid-binding proteins (I-BABPs), also called ileal lipid-binding proteins or gastrotropins, belong to the family of the fatty acid-binding proteins and play an important role in the solubilization and transport of bile acids in the enterocyte. This article describes the expression, purification, crystallization, and three-dimensional structure determination of zebrafish (Danio rerio) I-BABP both in its apo form and bound to cholic acid. This is the first X-ray structure of an I-BABP. The structure of the apoprotein was determined to a resolution of 1.6 Å, and two different monoclinic crystal forms of the holoprotein were solved and refined to 2.2 Å resolution. Three protein molecules are present in the asymmetric unit of one of the co-crystal forms and two in the other, and therefore, the results of this study refer to observations made on five different protein molecules in the crystalline state. In every case, two cholate ligands were found bound in approximately the same position in the internal cavity of the protein molecules, but an unexpected result is the presence of clear and unambiguous electron density for several cholate molecules bound on hydrophobic patches on the surface of all the five independent protein molecules examined. Isothermal titration calorimetry was used for the thermodynamic characterization of the binding mechanism and has yielded results that are consistent with the X-ray data. Ligand binding is described in detail, and the conformational changes undergone by the protein molecule in the apo-to-holo transition are examined by superposition of the apo- and holoprotein models. The structure of the holoprotein is also compared with that of the liver BABP from the same species and those of other I-BABPs determined by NMR.     

FXRE can function as an LXRE in the promoter of human ileal bile acid-binding protein (I-BABP) gene

Ileal bile acid-binding protein (I-BABP) is a 14 kDa cytosolic protein which binds bile acids with a high affinity. It is thought to be implicated in the enterohepatic circulation of bile acids and, hence, in cholesterol homeostasis. Using a combination of in vivo and in vitro experiments, we have recently shown that I-BABP gene expression can be indirectly up-regulated by cholesterol through the activation of sterol-responsive element-binding protein 1c (SREBP1c) by liver X-receptor (LXR). We report here that I-BABP can be also a direct target for LXR. I-BABP regulation by LXR is maintained when the SREBP binding site is deleted in the I-BABP promoter and occurs, in the absence of conventional LXRE sequences, through an IR1 sequence previously identified as a farnesoid X-receptor-responsive element (FXRE). Electrophoretic mobility shift assays demonstrated that the LXR/RXR heterodimer specifically recognizes the FXRE. Collectively, these data strongly suggest that LXR can regulate the I-BABP gene by both direct and indirect mechanisms.     

Dynamics of backbone conformational heterogeneity in Bacillus subtilis ribonuclease P protein

Interconversion of protein conformations is imperative to function, as evidenced by conformational changes associated with enzyme catalytic cycles, ligand binding and post-translational modifications. In this study, we used 15N NMR relaxation experiments to probe the fast (i.e., ps-ns) and slow (i.e., micros-ms) conformational dynamics of Bacillus subtilis ribonuclease P protein (P protein) in its folded state, bound to two sulfate anions. Using the Lipari-Szabo mapping method [Andrec, M., Montelione, G. T., and Levy, R. M. (2000) J. Biomol. NMR 18, 83-100] to interpret the data, we find evidence for P protein dynamics on the mus-ms time scale in the ensemble. The residues that exhibit these slow internal motions are found in regions that have been previously identified as part of the P protein-P RNA interface. These results suggest that structural flexibility within the P protein ensemble may be important for proper RNase P holoenzyme assembly and/or catalysis.     

NMR-based modeling and binding studies of a ternary complex between chicken liver bile acid binding protein and bile acids

Chicken liver bile acid binding protein (cL-BABP) is involved in bile acid transport in the liver cytosol. A detailed study of the mechanism of binding and selectivity of bile acids binding proteins towards the physiological pool of bile salts is a key issue for the complete understanding of the role of these proteins and their involvement in cholesterol homeostasis. In the present study, we modeled the ternary complex of cL-BABP with two molecules of bile salts using the data driven docking program HADDOCK on the basis of NMR and mass spectrometry data. Docking resulted in good 3D models, satisfying the majority of experimental restraints. The docking procedure represents a necessary step to help in the structure determination and in functional analysis of such systems, in view of the high complexity of the 3D structure determination of a ternary complex with two identical ligands. HADDOCK models show that residues involved in binding are mainly located in the C-terminal end of the protein, with two loops, CD and EF, playing a major role in ligand binding. A spine, comprising polarresidues pointing toward the protein interior and involved in motion communication, has a prominent role in ligand interaction. The modeling approach has been complemented with NMR interaction and competition studies of cL-BABP with chenodeoxycholic and cholic acids. A higher affinity for chenodeoxycholic acid was observed and a Kd upper limit estimate was obtained. The binding is highly cooperative and no site selectivity was detected for the different bile salts, thus indicating that site selectivity and cooperativity are not correlated. Differences in physiological pathways and bile salt pools in different species is discussed in light of the binding results thus enlarging the body of knowledge of BABPs biological functions.     

Modeling conformational ensembles of slow functional motions in Pin1-WW

Protein-protein interactions are often mediated by flexible loops that experience conformational dynamics on the microsecond to millisecond time scales. NMR relaxation studies can map these dynamics. However, defining the network of inter-converting conformers that underlie the relaxation data remains generally challenging. Here, we combine NMR relaxation experiments with simulation to visualize networks of inter-converting conformers. We demonstrate our approach with the apo Pin1-WW domain, for which NMR has revealed conformational dynamics of a flexible loop in the millisecond range. We sample and cluster the free energy landscape using Markov State Models (MSM) with major and minor exchange states with high correlation with the NMR relaxation data and low NOE violations. These MSM are hierarchical ensembles of slowly interconverting, metastable macrostates and rapidly interconverting microstates. We found a low population state that consists primarily of holo-like conformations and is a "hub" visited by most pathways between macrostates. These results suggest that conformational equilibria between holo-like and alternative conformers pre-exist in the intrinsic dynamics of apo Pin1-WW. Analysis using MutInf, a mutual information method for quantifying correlated motions, reveals that WW dynamics not only play a role in substrate recognition, but also may help couple the substrate binding site on the WW domain to the one on the catalytic domain. Our work represents an important step towards building networks of inter-converting conformational states and is generally applicable.     

Characterization of a novel 14 kDa bile acid-binding protein from rat ileal cytosol

A 14 kDa polypeptide in rat ileal cytosol has been identified as the major intestinal cytosolic bile acid-binding protein (I-BABP) by photoaffinity labeling with the radiolabeled 7,7-azo derivative of taurocholate (7,7-azo-TC). To further characterize I-BABP, the protein was purified by lysylglycocholate Sepharose 4B affinity and DE-52 anion-exchange chromatography. The purified I-BABP contained a single 14 kDa band on SDS-PAGE. The 14 kDa protein showed a 26-fold increase in binding affinity for [3H]7,7-azo-TC compared to cytosolic protein. Immunoblotting of protein fractions separated by affinity chromatography showed that neither liver fatty acid binding protein (L-FABP) nor intestinal fatty acid binding protein (I-FABP) bind to the affinity column and that the 14 kDa protein which bound to the column and was subsequently eluted with detergent did not cross-react with anti-L-FABP or anti-I-FABP. The 14 kDa protein labeled with [3H]7,7-azo-TC was radioimmunoprecipitated from cytosol by rabbit antiserum raised against purified I-BABP. I-BABP was shown to have a blocked N-terminus; however, its mixed internal sequence generated from cyanogen bromide-cleaved protein and amino acid composition indicated that it was related to (although clearly distinct from) both I-FABP and L-FABP. These studies have isolated a 14 kDa bile acid-binding protein from rat ileal cytosol which is immunologically and biochemically distinct from I-FABP and L-FABP.     

Temperature-dependent spectral density analysis applied to monitoring backbone dynamics of major urinary protein-I complexed with the pheromone 2-sec-butyl-4,5-dihydrothiazole*

Backbone dynamics of mouse major urinary protein I (MUP-I) was studied by (15)N NMR relaxation. Data were collected at multiple temperatures for a complex of MUP-I with its natural pheromonal ligand, 2- sec -4,5-dihydrothiazole, and for the free protein. The measured relaxation rates were analyzed using the reduced spectral density mapping. Graphical analysis of the spectral density values provided an unbiased qualitative picture of the internal motions. Varying temperature greatly increased the range of analyzed spectral density values and therefore improved reliability of the analysis. Quantitative parameters describing the dynamics on picosecond to nanosecond time scale were obtained using a novel method of simultaneous data fitting at multiple temperatures. Both methods showed that the backbone flexibility on the fast time scale is slightly increased upon pheromone binding, in accordance with the previously reported results. Zero-frequency spectral density values revealed conformational changes on the microsecond to millisecond time scale. Measurements at different temperatures allowed to monitor temperature dependence of the motional parameters.     

NMR investigations of allosteric processes in a two-domain Thermus thermophilus Hsp70 molecular chaperone

Hsp70 chaperones are two-domain proteins that assist in intra-cellular protein (re) folding processes in all species. The protein folding activity of the substrate binding domain of the Hsp70s is regulated by nucleotide binding at the nucleotide-binding domain through an as yet undefined heterotropic allosteric mechanism. The available structures of the isolated domains of Hsp70s have given very limited indications of nucleotide-induced conformational changes that could modulate the affinity for substrate proteins. Here, we present a multi-dimensional NMR study of a prokaryotic Hsp70 homolog, Thermus thermophilus DnaK, using a 54kDa construct containing both nucleotide binding domain and most of the substrate binding domain. It is determined that the nucleotide binding domain and substrate binding domain are closely associated in all ligand states studied. Comparison of the assigned NMR spectra of the two-domain construct with those of the previously studied isolated nucleotide binding domain, allowed the identification of the nucleotide binding domain-substrate binding domain interface. A global three-dimensional structure was obtained for the two-domain construct on the basis of this information and of NMR residual dipolar couplings measurements. This is the first experimental elucidation of the relative positioning of the nucleotide binding domain and substrate binding domain for any Hsp70 chaperone. Comparisons of NMR data between various ligand states including nucleotide-free, ATP, ADP.Pi and ADP.Pi+ peptide bound, identified residues involved in the allosteric inter-domain communication. In particular, peptide binding to the substrate binding domain was found to cause conformational changes in the NBD extending to the nucleotide binding pocket. Detailed analysis suggests that the inter-domain interface becomes tighter in the (nucleotide binding domain ligation/substrate binding domain ligation) order ATP/apo, ADP.Pi/apo ADP.Pi/peptide.     

Conformational change upon ligand binding and dynamics of the PDZ domain from leukemia-associated Rho guanine nucleotide exchange factor

Leukemia-associated Rho guanine nucleotide exchange factor (LARG) is a RhoA-specific guanine nucleotide exchange factor (GEF) that can activate RhoA. The PDZ (PSD-95/Disc-large/ZO-1 homology) domain of LARG interacts with membrane receptors, which can relay extracellular signals to RhoA signal transduction pathways. Until now there is no structural and dynamic information about these interactions. Here we report the NMR structures of the LARG PDZ in the apo form and in complex with the plexin-B1 C-terminal octapeptide. Unobservable resonances of the residues in betaB/betaC and betaE/alphaB loops in apo state were observed in the complex state. A distinct region of the binding groove in the LARG PDZ was found to undergo conformational change compared with other PDZs. Analysis of the (15)N relaxation data using reduced spectral density mapping shows that the apo LARG PDZ (especially its ligand-binding groove) is flexible and exhibits internal motions on both picosecond to nanosecond and microsecond to millisecond timescales. Mutagenesis and thermodynamic studies indicate that the conformation of the betaB/betaC and betaE/alphaB loops affects the PDZ-peptide interaction. It is suggested that the conformational flexibility could facilitate the change of structures upon ligand binding.     

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.     

Ligand-induced conformational change of a protein reproduced by a linear combination of displacement vectors obtained from normal mode analysis

The conformational change of a protein upon ligand binding was examined by normal mode analysis (NMA) based on an elastic-network model (ENM) for a full-atom system using dihedral angles as independent variables. Specifically, we investigated the extent to which conformational change vectors of atoms from an apo form to a holo form of a protein can be represented by a linear combination of the displacement vectors of atoms in the apo form calculated for the lowest-frequency m normal modes (m=1, 2,…, 20). In this analysis, the latter vectors were best fitted to the former ones by the least-squares method. Twenty-two paired proteins in the holo and apo forms, including three dimer pairs, were examined. The results showed that, in most cases, the conformational change vectors were reproduced well by a linear combination of the displacement vectors of a small number of low-frequency normal modes. The conformational change around an active site was reproduced as well as the entire conformational change, except for some proteins that only undergo significant conformational changes around active sites. The weighting factors for 20 normal modes optimized by the least-squares fitting characterize the conformational changes upon ligand binding for these proteins. The conformational changes sampled around the apo form of a protein by the linear combination of the displacement vectors obtained by ENM-based NMA may help solve the flexible-docking problem of a protein with another molecule because the results presented herein suggest that they have a relatively high probability of being involved in an actual conformational change.