Evolution of the family of intracellular lipid binding proteins in vertebrates

Members of the family of intracellular lipid binding proteins (iLBPs) have been implicated in cytoplasmic transport of lipophilic ligands, such as long-chain fatty acids and retinoids. iLBPs are low molecular mass proteins (14–16 kDa) sharing a common structural fold. The iLBP family likely arose through duplication and diversification of an ancestral iLBP gene. Phylogenetic analysis undertaken in the present study indicates that the ancestral iLBP gene arose after divergence of animals from fungi and plants. The first gene duplication was dated around 930 millions of years ago, and subsequent duplications in the succeeding 550 millions of years gave rise to the 16 iLBP types currently recognized in vertebrates. Four clusters of proteins, each binding a characteristic range of ligands, are evident from the phylogenetic tree. Evolution of different binding properties probably allowed cytoplasmic trafficking of distinct ligands. It is speculated that recruitment of an iLBP during evolution of animals enabled the mitochondrial oxidation of long-chain fatty acids.     

New insights into intracellular lipid binding proteins: The role of buried water

The crystal structures of most intracellular lipid binding proteins (LBPs) show between 5 and 20 internally bound water molecules, depending on the presence or the absence of ligand inside the protein cavity. The structural and functional significance of these waters has been discussed for several LBPs based on studies that used various biophysical techniques. The present work focuses on two very different LBPs, heart-type fatty acid binding protein (H-FABP) and ileal lipid binding protein (ILBP). Using high-resolution nuclear magnetic resonance spectroscopy, certain resonances belonging to side-chain protons that are located inside the water-filled lipid binding cavity were observed. In the case of H-FABP, the pH- and temperature-dependent behavior of selected side-chain resonances (Ser82 OgH and the imidazole ring protons of His93) indicated an unusually slow exchange with the solvent, implying that the intricate hydrogen-bonding network of amino-acid side-chains and water molecules in the protein interior is very rigid. In addition, holo H-FABP appeared to display a reversible self-aggregation at physiological pH. For ILBP, on the other hand, a more solvent-accessible protein cavity was deduced based on the pH titration behavior of its histidine residues. Comparison with data from other LBPs implies that the evolutionary specialization of LBPs for certain ligand types was not only because of mutations of residues directly involved in ligand binding but also to a refinement of the internal water scaffold.     

A comparative study of the backbone dynamics of two closely related lipid binding proteins: Bovine heart fatty acid binding protein and porcine ileal lipid binding protein

The backbone dynamics of bovine heart fatty acid binding protein (H-FABP) and porcine ileal lipid binding protein (ILBP) were studied by 15N NMR relaxation (T1 and T2) and steady state heteronuclear 15N{1H} NOE measurements. The microdynamic parameters characterizing the backbone mobility were determined using the model-free approach. For H-FABP, the non-terminal backbone amide groups display a rather compact protein structure of low flexibility. In contrast, for ILBP an increased number of backbone amide groups display unusually high internal mobility. Furthermore, the data indicate a higher degree of conformational exchange processes in the sec-msec time range for ILBP compared to H-FABP. These results suggest significant differences in the conformational stability for these two structurally highly homologous members of the fatty acid binding protein family.     

Effects of ligand binding on the dynamics of rice nonspecific lipid transfer protein 1: a model from molecular simulations

Plant nonspecific lipid transfer proteins (nsLTPs) are small, basic proteins constituted mainly of alpha-helices and stabilized by four conserved disulfide bridges. They are characterized by the presence of a tunnel-like hydrophobic cavity, capable of transferring various lipid molecules between lipid bilayers in vitro. In this study, molecular dynamics (MD) simulations were performed at room temperature to investigate the effects of lipid binding on the dynamic properties of rice nsLTP1. Rice nsLTP1, either in the free form or complexed with one or two lipids was subjected to MD simulations. The C-terminal loop was very flexible both before and after lipid binding, as revealed by calculating the root-mean-square fluctuation. After lipid binding, the flexibility of some residues that were not in direct contact with lipid molecules increased significantly, indicating an increase of entropy in the region distal from the binding site. Essential dynamics analysis revealed clear differences in motion between unliganded and liganded rice nsLTP1s. In the free form of rice nsLTP1, loop1 exhibited the largest directional motion. This specific essential motion mode diminished after binding one or two lipid molecules. To verify the origin of the essential motion observed in the free form of rice nsLTP1, we performed multiple sequence alignments to probe the intrinsic motion encoded in the primary sequence. We found that the amino acid sequence of loop1 is highly conserved among plant nsLTP1s, thus revealing its functional importance during evolution. Furthermore, the sequence of loop1 is composed mainly of amino acids with short side chains. In this study, we show that MD simulations, together with essential dynamics analysis, can be used to determine structural and dynamic differences of rice nsLTP1 upon lipid binding.     

Aluminum binding to phosphatidylcholine lipid bilayer membranes: 27Al and 31P NMR spectroscopic studies

27Al and 31P nuclear magnetic resonance (NMR) spectroscopies were used to investigate aluminum interactions at pH 3.4 with model membranes composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC). A solution state 27Al NMR difference assay was developed to quantify aluminum binding to POPC multilamellar vesicles (MLVs). Corresponding one-dimensional (1D) fast magic angle spinning (MAS) 31P NMR spectra showed that aluminum induced the appearance of two new isotropic resonances for POPC shifted to -6.4 ppm and -9.6 ppm upfield relative to, and in slow exchange with, the control resonance at -0.6 ppm. Correlation of the (27)Al and (31)P NMR binding data revealed a 1:2 aluminum:phospholipid stoichiometry in the aluminum-bound complex at -9.6 ppm and a 1:1 aluminum:phospholipid stoichiometry in that at -6.4 ppm. Slow MAS 31P NMR spectra demonstrated shifts in the anisotropic chemical shift tensor components of the aluminum-bound POPC consistent with a close coordination of aluminum with phosphorus. A model of the aluminum-bis-phospholipid complex is proposed on the basis of these findings.     

Interdomain dynamics and ligand binding: molecular dynamics simulations of glutamine binding protein

Periplasmic binding proteins from Gram-negative bacteria possess a common architecture, comprised of two domains linked by a hinge region, a fold which they share with the neurotransmitter-binding domains of ionotropic glutamate receptors (GluRs). Glutamine-binding protein (GlnBP) is one such protein, whose crystal structure has been solved in both open and closed forms. Multi-nanosecond molecular dynamics simulations have been used to explore motions about the hinge region and how they are altered by ligand binding. Glutamine binding is seen to significantly reduce inter-domain motions about the hinge region. Essential dynamics analysis of inter-domain motion revealed the presence of both hinge-bending and twisting motions, as has been reported for a related sugar-binding protein. Significantly, the influence of the ligand on GlnBP dynamics is similar to that previously observed in simulations of rat glutamate receptor (GluR2) ligand-binding domain. The essential dynamics analysis of GlnBP also revealed a third class of motion which suggests a mechanism for signal transmission in GluRs.     

The liver fatty acid binding protein - comparison of cavity properties of intracellular lipid-binding proteins

The crystal and solution structures of all of the intracellular lipid binding proteins (iLBPs) reveal a common -barrel framework with only small local perturbations. All existing evidence points to the binding cavity and a poorly delimited portal region as defining the function of each family member. The importance of local structure within the cavity appears to be its influence on binding affinity and specificity for the lipid. The portal region appears to be involved in the regulation of ligand exchange. Within the iLBP family, liver fatty acid binding protein or LFABP, has the unique property of binding two fatty acids within its internalized binding cavity rather than the commonly observed stoichiometry of one. Furthermore, LFABP will bind hydrophobic molecules larger than the ligands which will associate with other iLBPs. The crystal structure of LFABP contains two bound oleate molecules and provides the explanation for its unusual stoichiometry. One of the bound fatty acids is completely internalized and has its carboxylate interacting with an arginine and two serines. The second oleate represents an entirely new binding mode with the carboxylate on the surface of LFABP. The two oleates also interact with each other. Because of this interaction and its inner location, it appears the first oleate must be present before the second more external molecule is bound.     

Molecular dynamics simulations of apo and holo forms of fatty acid binding protein 5 and cellular retinoic acid binding protein II reveal highly mobile protein,retinoic acid ligand,and water molecules

Structural and dynamic properties from a series of 300 ns molecular dynamics, MD, simulations of two intracellular lipid binding proteins, iLBPs, (Fatty Acid Binding Protein 5, FABP5, and Cellular Retinoic Acid Binding Protein II, CRABP-II) in both the apo form and when bound with retinoic acid reveal a high degree of protein and ligand flexibility. The ratio of FABP5 to CRABP-II in a cell may determine whether it undergoes natural apoptosis or unrestricted cell growth in the presence of retinoic acid. As a result, FABP5 is a promising target for cancer therapy. The MD simulations presented here reveal distinct differences in the two proteins and provide insight into the binding mechanism. CRABP-II is a much larger, more flexible protein that closes upon ligand binding, where FABP5 transitions to an open state in the holo form. The traditional understanding obtained from crystal structures of the gap between two β-sheets of the β-barrel common to iLBPs and the α-helix cap that forms the portal to the binding pocket is insufficient for describing protein conformation (open vs. closed) or ligand entry and exit. When the high degree of mobility between multiple conformations of both the ligand and protein are examined via MD simulation, a new mode of ligand motion that improves understanding of binding dynamics is revealed.     

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.     

Obtaining ligand geometries from paramagnetic shifts in low-spin haem proteins

b and c with His, Met, and cyanide ligands. Variations in the electronic structure of the haem and the magnetic susceptibility tensors have been shown to depend primarily on the axial ligand geometry, and the shifts of haem substituents have been used to obtain the first structural information for several cytochromes. Recently, the database of assigned spectra for bis-His haems has been extended sufficiently for an empirical equation to be produced for treating ¹H NMR data from haem methyl groups at 298 K. However, the database used contains large systematic deviations and the form of the equation leads to systematic errors in the ligand geometries. This article describes the link with the semi-empirical methods used previously and provides a set of corrected empirical parameters as well as an improved equation. The possibilities for generalising the empirical method to account for ligands other than His and temperatures other than 298 K are discussed. Received: 19 October 1999 / Accepted: 21 February 2000     

The influence of proteins and peptides on the phase properties of lipids

This paper reviews model membrane studies on the modulation of the macroscopic structure of lipids by lipid-protein interactions, with particular emphasis on the gramicidin molecule. This hydrophobic peptide has three main effects on lipid polymorphism: (1) in lysophosphatidylcholine it triggers a micellar to bilayer transition, (2) in phosphatidylethanolamine it lowers the bilayer to hexagonal H_II phase transition temperature and (3) in phosphatidylcholine and other bilayer preferring lipids it is able to induce the formation of an H_II phase. From experiments in which the gramicidin molecule was chemically modified it can be concluded that the tryptophan residues play a determining role in the peptide-induced changes in polymorphism. The experimental data lead to the proposal that gramicidin molecules have a tendency to self-associate, possibly mediated by tryptophan-tryptophan interactions and organize into tubular structures such as found in the H_II phase.     

The effect of heterogeneous structural diffusion on ligand binding to heme proteins

The kinetics of ligand binding to heme proteins studied by flash photolysis display an algebraic time dependence at low temperatures in contrast exponential recombination observed under physiological conditions. This result shows that protein structures should be viewed as a time average of interconverting microstates which are frozen in at low temperatures. We propose a quasi-one-dimensional model of heterogeneous structural diffusion coupled to ligand binding which describes freezing transition as an inherent property of protein fluctuations. The structural hopping rates are derived from a temperature invariant spectrum of activation energies. The model predicts power law kinetics of the form t ^- at long times. The exponent is constant (0.5) at high temperatures but decreases below a critical temperature in the frozen regime. These results are compared to experiments performed with myoglobin and -chains of hemoglobin.     

A functional NMR for membrane proteins: dynamics,ligand binding,and allosteric modulation

By nature of conducting ions, transporting substrates and transducing signals, membrane channels, transporters and receptors are expected to exhibit intrinsic conformational dynamics. It is therefore of great interest and importance to understand the various properties of conformational dynamics acquired by these proteins, for example, the relative population of states, exchange rate, conformations of multiple states, and how small molecule ligands modulate the conformational exchange. Because small molecule binding to membrane proteins can be weak and/or dynamic, structural characterization of these effects is very challenging. This review describes several NMR studies of membrane protein dynamics, ligand‐induced conformational rearrangements, and the effect of ligand binding on the equilibrium of conformational exchange. The functional significance of the observed phenomena is discussed.     

Plant lipid binding proteins: properties and applications

Lipid transfer proteins (LTP) and puroindolines are abundant lipid binding proteins of plant seeds. While LTP are ubiquitous plant proteins, puroindolines are only found in the seeds of plants from the Triticae and Avenae tribes. These proteins display a similar overall folding pattern but different lipid binding properties. The unique and diverse biological and technological functions of LTPs and puroindolines are closely related to their structural and lipid binding properties. These proteins are attractive to improve the agronomic performances and food quality of crops. Heterologous expression and genetic engineering should allow industrial production and enlarge applications of these lipid binding proteins.     

Insights on pH-dependent conformational changes of mosquito odorant binding proteins by molecular dynamics simulations

Chemical recognition plays an important role for the survival and reproduction of many insect species. Odorant binding proteins (OBPs) are the primary components of the insect olfactory mechanism and have been documented to play an important role in the host-seeking mechanism of mosquitoes. They are “transport proteins” believed to transport odorant molecules from the external environment to their respective membrane targets, the olfactory receptors. The mechanism by which this transport occurs in mosquitoes remains a conundrum in this field. Nevertheless, OBPs have proved to be amenable to conformational changes mediated by a pH change in other insect species. In this paper, the effect of pH on the conformational flexibility of mosquito OBPs is assessed computationally using molecular dynamics simulations of a mosquito OBP “CquiOBP1” bound to its pheromone 3OG (PDB ID: 3OGN). Conformational twist of a loop, driven by a set of well-characterized changes in intramolecular interactions of the loop, is demonstrated. The concomitant (i) closure of what is believed to be the entrance of the binding pocket, (ii) expansion of what could be an exit site, and (iii) migration of the ligand towards this putative exit site provide preliminary insights into the mechanism of ligand binding and release of these proteins in mosquitoes. The correlation of our results with previous experimental observations based on NMR studies help us provide a cardinal illustration on one of the probable dynamics and mechanism by which certain mosquito OBPs could deliver their ligand to their membrane-bound receptors at specific pH conditions.     

Hydrogen exchange and ligand binding: ligand-dependent and ligand-independent protection in the Src SH3 domain

Amide hydrogen-deuterium exchange has proven to be a powerful tool for detecting and characterizing high-energy conformations in protein ensembles. Since interactions with ligands can modulate these high-energy conformations, hydrogen exchange appears to be an ideal experimental probe of the physical mechanisms underlying processes like allosteric regulation. The chemical mechanism of hydrogen exchange, however, can complicate such studies. Here, we examine hydrogen exchange rates in a simple model system, the c-Src SH3 domain interacting with a short peptide ligand. Addition of ligand slows the rates of hydrogen exchange at nearly every amide for which we can obtain data. Careful analysis, however, reveals that this slowing is due primarily to a reduction in the population of free protein in the system, and not to any specific property of the complex. We present a method to separate the contributions of free and bound protein to the exchange kinetics that has allowed us to identify the subset of amides where exchange arises directly from the complex. These results demonstrate that the slowing of hydrogen exchange induced by ligand interactions should be interpreted with caution, and more extensive experiments are required to correlate changes in hydrogen exchange with changes in structure or internal dynamics.     

Localization of ligand binding site in proteins identified in silico

Knowledge-based models for protein folding assume that the early-stage structural form of a polypeptide is determined by the backbone conformation, followed by hydrophobic collapse. Side chain–side chain interactions, mostly of hydrophobic character, lead to the formation of the hydrophobic core, which seems to stabilize the structure of the protein in its natural environment. The fuzzy-oil-drop model is employed to represent the idealized hydrophobicity distribution in the protein molecule. Comparing it with the one empirically observed in the protein molecule reveals that they are not in agreement. It is shown in this study that the irregularity of hydrophobic distributions is aim-oriented. The character and strength of these irregularities in the organization of the hydrophobic core point to the specificity of a particular protein’s structure/function. When the location of these irregularities is determined versus the idealized fuzzy-oil-drop, function-related areas in the protein molecule can be identified. The presented model can also be used to identify ways in which protein–protein complexes can possibly be created. Active sites can be predicted for any protein structure according to the presented model with the free prediction server at . The implication based on the model presented in this work suggests the necessity of active presence of ligand during the protein folding process simulation. Figure Fuzzy-oil-drop model applied to identify the ligation site in lysozyme complexed with N-acetylglucosamine (PDB ID:1LMQ) in form of hydrophobicity deficiency (ΔH) profile and three-dimensional distribution of on protein surface     

Computational investigation of pressure profiles in lipid bilayers with embedded proteins

The distribution of surface tension within a lipid bilayer, also referred to as the lateral pressure profile, has been the subject of theoretical scrutiny recently due to its potential to radically alter the function of biomedically important membrane proteins. Experimental measurements of the pressure profile are still hard to come by, leaving first-principles all-atom calculations of the profile as an important investigative tool. We describe and validate an efficient implementation of pressure profile calculations in the molecular dynamics package NAMD, capable of distinguishing between internal, bonded and nonbonded contributions as well as those of selected atom groups. The new implementation can also be used in conjunction with Ewald summation for long-range electrostatics, improving the accuracy and reproducibility of the calculated profiles. We then describe results of the calculation of a pressure profile for a simple protein–lipid system consisting of melittin embedded in a DMPC bilayer. While the lateral pressure in the protein–lipid system is nearly the same as that of the bilayer alone, partitioning of the lateral pressure by atom type revealed substantial perturbation of the pressure profile and surface tension in an asymmetric manner.     

Structural and dynamic roles of permanent water molecules in ligand molecular recognition by chicken liver bile acid binding protein

Chicken liver bile acid binding protein (cL-BABP) crystallizes with water molecules in its binding site. To obtain insights on the role of internal water, we performed two 100 ns molecular dynamics (MD) simulations in explicit solvent for cL-BABP, as apo form and as a complex with two molecules of cholic acid, and analyzed in detail the dynamics properties of all water molecules. The diffusion coefficients of the more persistent internal water molecules are significantly different from the bulk, but similar between the two protein forms. A different number of molecules and a different organization are observed for apo- and holo-cL-BABP. Most water molecules identified in the binding site of the apo-crystal diffuse to the bulk during the simulation. In contrast, almost all the internal waters of the holo-crystal maintain the same interactions with internal sidechains and ligands, which suggests they have a relevant role in protein-ligand molecular recognition. Only in the presence of these water molecules we were able to reproduce, by a classical molecular docking approach, the structure of the complex cL-BABP::cholic acid with a low ligand root mean square deviation (RMSD) with respect to its reference positioning. Literature data reported a conserved pattern of hydrogen bonds between a single water molecule and three amino acid residues of the binding site in a series of crystallized FABP. In cL-BABP, the interactions between this conserved water molecule and the three residues are present in the crystal of both apo- and holo-cL-BABP but are lost immediately after the start of molecular dynamics. Copyright (c) 2008 John Wiley & Sons, Ltd.