Search

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., ((c) 2008 Wiley-Liss, Inc.)

It has been proposed that intracellular carrier proteins mediate active transport of the bile acids within hepatocytes and ileocytes, during the enterohepatic circulation. In mammalian species only ileal bile acid binding proteins have been so far identified, while liver cytosolic carriers have never been found. On the contrary, in non-mammalian vertebrates, only liver, and not ileal, bile acid binding proteins were reported. The aim of the present work is to find the missing cytosolic transport proteins. A bioinformatic search allowed us to identify a non-mammalian putative bile acid binding protein in the chicken ileum (cI-BABP), which we recombinantly expressed and purified. The protein exhibits the capability, tested by in vitro NMR experiments, of binding bile acids. Furthermore, strong NMR evidence reported that the human liver fatty acid binding protein (hL-FABP) can also bind bile acids. Taken together, these data strongly suggest that both cI-BABP and hL-FABP have a bile acid binding function in the two organisms, and support a previous hypothesis on the role of hL-FABP in regulating bile acid metabolism and determining bile acid pool size., ((c) 2007 Wiley-Liss, Inc.)

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., (2007 Wiley-Liss, Inc.)

A new energy decomposition approach, aimed at identifying residues playing a folding key role, has been applied here to three homologous proteins, belonging to the calycin superfamily, namely bovine and porcine beta-lactoglobulins and Liver basic fatty acid binding protein, sharing the same beta-barrel fold and different degree of sequence identities. All-atom, explicit solvent molecular dynamics simulations around the native conformation were used to generate, for each of the three proteins, energy maps which were further simplified through eigenvalue decomposition. Analysis of the components of the eigenvector associated with the lowest eigenvalue singled out those residues (hot sites) behaving as strongly interacting and possible nucleation centers. The results fit well with experimental folding data and, especially, with the analysis of side chain-side chain interaction conservation., (Copyright 2005 Wiley-Liss, Inc.)

The pH-driven opening and closure of beta-lactoglobulin EF loop, acting as a lid and closing the internal cavity of the protein, has been studied by molecular dynamics (MD) simulations and free energy calculations based on molecular mechanics/Poisson-Boltzmann (PB) solvent-accessible surface area (MM/PBSA) methodology. The forms above and below the transition pH differ presumably only in the protonation state of residue Glu89. MM/PBSA calculations are able to reproduce qualitatively the thermodynamics of the transition. The analysis of MD simulations using a combination of MM/PBSA methodology and the colony energy approach is able to highlight the driving forces implied in the transition. The analysis suggests that global rearrangements take place before the equilibrium local conformation is reached. This conclusion may bear general relevance to conformational transitions in all lipocalins and proteins in general., ((c) 2005 Wiley-Liss, Inc.)

The chemical unfolding behavior of porcine beta-lactoglobulin (PLG) has been followed at pH 2 and 6 in the presence of guanidinium hydrochloride. The PLG unfolding transition, monitored by tryptophan fluorescence, far and near UV circular dichroism and 1D-NMR, can be described by a three-state transition suggesting the presence of at least one intermediate state that appears to display an excess of non-native alpha-helical structures. The thermodynamic parameters, as determined through a global analysis fitting procedure, give estimates of the free energy differences of the transitions connecting the native, the intermediate and the unfolded state: DeltaG(NI) (0) = 2.8 +/- 0.7 kcal mol(-1) (pH 2) and 4.2 +/- 0.5 kcal mol(-1) (pH 6) and DeltaG(NU) (0) = 7.2 +/- 0.6 kcal mol(-1) (pH 2) and 6.9 +/- 0.6 kcal mol(-1) (pH 6). CD unfolding data of the bovine species (BLG) have been collected here under the same experimental conditions of PLG to allow a careful comparison of the two beta-lactoglobulins. Intermediates with different characteristics have been identified for BLG and PLG, and their nature has been discussed on a structural analysis basis. The thermodynamic data reported here for PLG and BLG and the comparative analysis with data reported for equine beta lactoglobulin, show that homologous beta-barrel proteins, belonging to the same family and displaying high sequence identity (52-64%) populate unfolding intermediates to different extents, even though a common tendency to the formation of non-native alpha-helical intermediates, can be envisaged. The present results provide a prerequisite foundation of knowledge for the design and interpretation of future folding kinetic studies., ((c) 2004 Wiley-Liss, Inc.)

Bovine beta-Lactoglobulin (BLG) has been studied for many decades, but only recently structural data have been obtained, making it possible to simulate its molecular properties. In the present study, electrostatic properties of BLG are investigated theoretically using Poisson-Boltzmann calculations and experimentally following pH titration via NMR. Electrostatic properties are determined for several structural models, including an ensemble of NMR structures obtained at low pH. The changes in electrostatic forces upon changes in ionic strength, solvent dielectric constant, and pH are calculated and compared with experiments. pK(a)s are computed for all titratable sites and compared with NMR titration data. The analysis of theoretical and experimental results suggests that (1) there may be more than one binding sites for negatively charged ligands; (2) at low pH the core of the molecule is more compact than observed in the structures obtained via restrained molecular dynamics from NMR data, but loop and terminal regions must be disordered.