Effects of ligand binding on dynamics of fatty acid binding protein and interactions with membranes

Intracellular transport of fatty acids involves binding of ligands to their carrier fatty acid binding proteins (FABPs) and interactions of ligand-free and -bound FABPs with membranes. Previous studies focused on ligand-free FABPs. Here, our amide hydrogen exchange data showed that oleic acid binding to human intestinal FABP (hIFABP) stabilizes the protein, most likely through enhancing the hydrogen-bonding network, and induces rearrangement of sidechains even far away from the ligand binding site. Using NMR relaxation techniques, we found that the ligand binding affects not only conformational exchanges between major and minor states but also the affinity of hIFABP to nanodiscs. Analyses of the relaxation and amide exchange data suggested that two minor native-like states existing in both ligand-free and -bound hIFABPs originate from global “breathing” motions, while one minor native-like state comes from local motions. The amide hydrogen exchange data also indicated that helix αII undergoes local unfolding through which ligands can exit from the binding cavity.     

Molecular Dynamics Simulation of Ligand Dissociation from Liver Fatty Acid Binding Protein

The mechanisms of how ligands enter and leave the binding cavity of fatty acid binding proteins (FABPs) have been a puzzling question over decades. Liver fatty acid binding protein (LFABP) is a unique family member which accommodates two molecules of fatty acids in its cavity and exhibits the capability of interacting with a variety of ligands with different chemical structures and properties. Investigating the ligand dissociation processes of LFABP is thus a quite interesting topic, which however is rather difficult for both experimental approaches and ordinary simulation strategies. In the current study, random expulsion molecular dynamics simulation, which accelerates ligand motions for rapid dissociation, was used to explore the potential egress routes of ligands from LFABP. The results showed that the previously hypothesized “portal region” could be readily used for the dissociation of ligands at both the low affinity site and the high affinity site. Besides, one alternative portal was shown to be highly favorable for ligand egress from the high affinity site and be related to the unique structural feature of LFABP. This result lends strong support to the hypothesis from the previous NMR exchange studies, which in turn indicates an important role for this alternative portal. Another less favored potential portal located near the N-terminal end was also identified. Identification of the dissociation pathways will allow further mechanistic understanding of fatty acid uptake and release by computational and/or experimental techniques.     

Molecular simulations of β‐lactoglobulin complexed with fatty acids reveal the structural basis of ligand affinity to internal and possible external binding sites

The interaction of saturated fatty acids of different length (C8:0 to C18:0) with β‐lactoglobulin (βLG) was investigated by molecular dynamics simulation and docking approaches. The results show that the presence of such ligands in the hydrophobic central cavity of βLG, known as the protein calyx, determines an enhancement of atomic fluctuations compared with the unliganded form, especially for loops at the entrance of the binding site. Concerted motions are evidenced for protein regions that could favor the binding of ligands. The mechanism of anchoring of fatty acids of different length is similar for the carboxylate head‐group, through electrostatic interactions with the side chains of Lys60/Lys69. The key protein residues to secure the hydrocarbon chain are Phe105/Met107, which adapt their conformation upon ligand binding. In particular, Phe105 provides an additional hydrophobic clamp only for the tail of the two fatty acids with the longest chains, palmitic, and stearic acid, which are known to bind βLG with a high affinity. The search of additional external binding sites for fatty acids, distinct from the calyx, was also carried out for palmitic acid. Two external sites with a lower affinity were identified as secondary sites, one consisting in a hydrophobic cavity allowing two distinct binding modes for the fatty acid, and the other corresponding to a surface crevice close to the protein α‐helix. The overall results provide a comprehensive picture of the dynamical behavior of βLG in complex with fatty acids, and elucidate the structural basis of the binding of these physiological ligands. Proteins 2014; 82:2609–2619. © 2014 Wiley Periodicals, Inc.     

Titration and exchange studies of liver fatty acid-binding protein with 13C-labeled long-chain fatty acids

Uniformly (13)C-labeled long-chain fatty acids were used to probe ligand binding to rat liver fatty acid-binding protein (LFABP), an atypical member of the fatty acid-binding protein (FABP) family that binds more than one molecule of long-chain fatty acid, accommodates a variety of diverse ligands, and exhibits diffusion-mediated lipid transport to membranes. Two sets of (1)H-(13)C resonances were found in a titration series of NMR spectra for oleate-LFABP complexes, indicating that two molecules of the fatty acid are situated in the protein cavity. However, no distinct resonances were observed for the excess fatty acid in solution, suggesting that at least one ligand undergoes rapid exchange with oleate in the bulk solution. An exchange rate of 54 +/- 6 s(-1) between the two sets of resonances was measured directly using (13)C z,z-exchange spectroscopy. In light of these NMR measurements, possible molecular mechanisms for the ligand-exchange process are evaluated and implications for the anomalous fatty acid transport mechanism of LFABP are discussed.     

Fluorine-19 NMR studies on the acid state of the intestinal fatty acid binding protein

The intestinal fatty acid binding protein (IFABP) is composed of two beta-sheets with a large hydrophobic cavity into which ligands bind. After eight 4-(19)F-phenylalanines were incorporated into the protein, the acid state of both apo- and holo-IFABP (at pH 2.8 and 2.3) was characterized by means of (1)H NMR diffusion measurements, circular dichroism, and (19)F NMR. Diffusion measurements show a moderately increased hydrodynamic radius while near- and far-UV CD measurements suggest that the acid state has substantial secondary structure as well as persistent tertiary interactions. At pH 2.8, these tertiary interactions have been further characterized by (19)F NMR and show an NOE cross-peak between residues that are located on different beta-strands. Side chain conformational heterogeneity on the millisecond time scale was captured by phase-sensitive (19)F-(19)F NOESY. At pH 2.3, native NMR peaks are mostly gone, but the protein can still bind fatty acid to form the holoprotein. An exchange cross-peak of one phenylalanine in the holoprotein is attributed to increased motional freedom of the fatty acid backbone caused by the slight opening of the binding pocket at pH 2.8. In the acid environment Phe128 and Phe17 show dramatic line broadening and chemical shift changes, reflecting greater degrees of motion around these residues. We propose that there is a separation of specific regions of the protein that gives rise to the larger radius of hydration. Temperature and urea unfolding studies indicate that persistent hydrophobic clusters are nativelike and may account for the ability of ligand to bind and induce nativelike structure, even at pH 2.3.     

Crystal structure of the HNF4 alpha ligand binding domain in complex with endogenous fatty acid ligand

HNF4 alpha is an orphan member of the nuclear receptor family with prominent functions in liver, gut, kidney and pancreatic beta cells. We have solved the x-ray crystal structure of the HNF4 alpha ligand binding domain, which adopts a canonical fold. Two conformational states are present within each homodimer: an open form with alpha helix 12 (alpha 12) extended and collinear with alpha 10 and a closed form with alpha 12 folded against the body of the domain. Although the protein was crystallized without added ligands, the ligand binding pockets of both closed and open forms contain fatty acids. The carboxylic acid headgroup of the fatty acid ion pairs with the guanidinium group of Arg(226) at one end of the ligand binding pocket, while the aliphatic chain fills a long, narrow channel that is lined with hydrophobic residues. These findings suggest that fatty acids are endogenous ligands for HNF4 alpha and establish a framework for understanding how HNF4 alpha activity is enhanced by ligand binding and diminished by MODY1 mutations.     

Dynamic conformational equilibria in the physiological function of the Bombyx mori pheromone-binding protein

The Bombyx mori pheromone-binding protein (BmorPBP) undergoes a pH-dependent conformational transition from a form at basic pH, which contains an open cavity suitable for ligand binding (BmorPBP^B), to a form at pH 4.5, where this cavity is occupied by an additional helix (BmorPBP^A). This helix α7 is formed by the C-terminal dodecapeptide 131-142, which is flexibly disordered on the protein surface in BmorPBP^B and in its complex with the pheromone bombykol. Previous work showed that the ligand-binding cavity cannot accommodate both bombykol and helix α7. Here we further investigated mechanistic aspects of the physiologically crucial ejection of the ligand at lower pH values by solution NMR studies of the variant protein BmorPBP(1-128), where the C-terminal helix-forming tetradecapeptide is removed. The NMR structure of the truncated protein at pH 6.5 corresponds closely to BmorPBP^B. At pH 4.5, BmorPBP(1-128) maintains a B-type structure that is in a slow equilibrium, on the NMR chemical shift timescale, with a low-pH conformation for which a discrete set of ¹⁵N-¹H correlation peaks is NMR unobservable. The full NMR spectrum was recovered upon readjusting the pH of the protein solution to 6.5. These data reveal dual roles for the C-terminal tetradecapeptide of BmorPBP in the mechanism of reversible pheromone binding and transport, where it governs dynamic equilibria between two locally different protein conformations at acidic pH and competes with the ligand for binding to the interior cavity.     

A nuclear magnetic resonance-based structural rationale for contrasting stoichiometry and ligand binding site(s) in fatty acid-binding proteins

Liver fatty acid-binding protein (LFABP) is a 14 kDa cytosolic polypeptide, differing from other family members in the number of ligand binding sites, the diversity of bound ligands, and the transfer of fatty acid(s) to membranes primarily via aqueous diffusion rather than direct collisional interactions. Distinct two-dimensional (1)H-(15)N nuclear magnetic resonance (NMR) signals indicative of slowly exchanging LFABP assemblies formed during stepwise ligand titration were exploited, without determining the protein-ligand complex structures, to yield the stoichiometries for the bound ligands, their locations within the protein binding cavity, the sequence of ligand occupation, and the corresponding protein structural accommodations. Chemical shifts were monitored for wild-type LFABP and an R122L/S124A mutant in which electrostatic interactions viewed as being essential to fatty acid binding were removed. For wild-type LFABP, the results compared favorably with the data for previous tertiary structures of oleate-bound wild-type LFABP in crystals and in solution: there are two oleates, one U-shaped ligand that positions the long hydrophobic chain deep within the cavity and another extended structure with the hydrophobic chain facing the cavity and the carboxylate group lying close to the protein surface. The NMR titration validated a prior hypothesis that the first oleate to enter the cavity occupies the internal protein site. In contrast, (1)H and (15)N chemical shift changes supported only one liganded oleate for R122L/S124A LFABP, at an intermediate location within the protein cavity. A rationale based on protein sequence and electrostatics was developed to explain the stoichiometry and binding site trends for LFABPs and to put these findings into context within the larger protein family.     

Structural basis of non-specific lipid binding in maize lipid-transfer protein complexes revealed by high-resolution X-ray crystallography

Non-specific lipid-transfer proteins (nsLTPs) are involved in the movement of phospholipids, glycolipids, fatty acids, and steroids between membranes. Several structures of plant nsLTPs have been determined both by X-ray crystallography and nuclear magnetic resonance. However, the detailed structural basis of the non-specific binding of hydrophobic ligands by nsLTPs is still poorly understood. In order to gain a better understanding of the structural basis of the non-specific binding of hydrophobic ligands by nsLTPs and to investigate the plasticity of the fatty acid binding cavity in nsLTPs, seven high-resolution (between 1.3 A and 1.9 A) crystal structures have been determined. These depict the nsLTP from maize seedlings in complex with an array of fatty acids.A detailed comparison of the structures of maize nsLTP in complex with various ligands reveals a new binding mode in an nsLTP-oleate complex which has not been seen before. Furthermore, in the caprate complex, the ligand binds to the protein cavity in two orientations with equal occupancy. The volume of the hydrophobic cavity in the nsLTP from maize shows some variation depending on the size of the bound ligands.The structural plasticity of the ligand binding cavity and the predominant involvement of non-specific van der Waals interactions with the hydrophobic tail of the ligands provide a structural explanation for the non-specificity of maize nsLTP. The hydrophobic cavity accommodates various ligands from C10 to C18. The C18:1 ricinoleate with its hydroxyl group hydrogen bonding to Ala68 possibly mimics cutin monomer binding which is of biological importance. Some of the myristate binding sites in human serum albumin resemble the maize nsLTP, implying the importance of a helical bundle in accommodating the non-specific binding of fatty acids.     

Two fatty acid‐binding proteins expressed in the intestine interact differently with endocannabinoids

Two different members of the fatty acid‐binding protein (FABP) family are found in enterocyte cells of the gastrointestinal system, namely liver‐type and intestinal fatty acid‐binding proteins (LFABP and IFABP, also called FABP1 and FABP2, respectively). Striking phenotypic differences have been observed in knockout mice for either protein, for example, high fat‐fed IFABP‐null mice remained lean, whereas LFABP‐null mice were obese, correlating with differences in food intake. This finding prompted us to investigate the role each protein plays in directing the specificity of binding to ligands involved in appetite regulation, such as fatty acid ethanolamides and related endocannabinoids. We determined the binding affinities for nine structurally related ligands using a fluorescence competition assay, revealing tighter binding to IFABP than LFABP for all ligands tested. We found that the head group of the ligand had more impact on binding affinity than the alkyl chain, with the strongest binding observed for the carboxyl group, followed by the amide, and then the glycerol ester. These trends were confirmed using two‐dimensional ¹H–¹⁵N nuclear magnetic resonance (NMR) to monitor chemical shift perturbation of the protein backbone resonances upon titration with ligand. Interestingly, the NMR data revealed that different residues of IFABP were involved in the coordination of endocannabinoids than those implicated for fatty acids, whereas the same residues of LFABP were involved for both classes of ligand. In addition, we identified residues that are uniquely affected by binding of all types of ligand to IFABP, suggesting a rationale for its tighter binding affinity compared with LFABP.     

The ligand-mediated affinity of brain-type fatty acid-binding protein for membranes determines the directionality of lipophilic cargo transport

The intracellular transport of lipophilic cargoes is a highly dynamic process. In eukaryotic cells, the uptake and release of long-chain fatty acids (LCFAs) are executed by fatty-acid binding proteins. However, how these carriers control the directionality of cargo trafficking remains unclear. Here, we revealed that the unliganded archetypal Drosophila brain-type fatty acid-binding protein (dFABP) possesses a stronger binding affinity than its liganded counterpart for empty nanodiscs (ND). Titrating unliganded dFABP and nanodiscs with LCFAs rescued the broadening of FABP cross-peak intensities in HSQC spectra from a weakened protein–membrane interaction. Two out of the 3 strongest LCFA contacting residues in dFABP identified by NMR HSQC chemical shift perturbation (CSP) are also part of the 30 ND-contacting residues (out of the total 130 residues in dFABP), revealed by attenuated TROSY signal in the presence of lipid ND to apo-like dFABP. Our crystallographic temperature factor data suggest enhanced αII helix dynamics upon LCFA binding, compensating for the entropic loss in the βC-D/βE-F loops. The aliphatic tail of bound LCFA impedes the charge-charge interaction between dFABP and the head groups of the membrane, and dFABP is prone to dissociate from the membrane upon ligand binding. We therefore conclude that lipophilic ligands participate directly in the control of the functionally required membrane association and dissociation of FABPs.     

Model of β-Sheet of Muscle Fatty Acid Binding Protein of Locusta migratoria Displays Characteristic Topology

The β-sheet of muscle fatty acid binding protein of Locusta migratoria (Lm-FABP) was modeled by employing 2-D NMR data and the Rigid Body Assembly method. The model shows the β-sheet to comprise ten β-strands arranged anti-parallel to each other. There is a β-bulge between Ser 13 and Gln 14 which is a difference from the published structure of β-sheet of bovine heart Fatty Acid Binding Protein. Also, a hydrophobic patch consisting of Ile 45, Phe 51, Phe 64 and Phe 66 is present on the surface which is characteristic of most Fatty Acid Binding Proteins. A “gap” is present between βD and βE that provides evidence for the presence of a portal or opening between the polypeptide chains which allows ligand fatty acids to enter the protein cavity and bind to the protein.     

Conformational exchange of fatty acid binding protein induced by protein-nanodisc interactions

Fatty acid binding proteins (FABPs) can facilitate the transfer of long-chain fatty acids between intracellular membranes across considerable distances. The transfer process involves fatty acids, their donor membrane and acceptor membrane, and FABPs, implying that potential protein-membrane interactions exist. Despite intensive studies on FABP-membrane interactions, the interaction mode remains elusive, and the protein-membrane association and dissociation rates are inconsistent. In this study, we used nanodiscs (NDs) as mimetic membranes to investigate FABP-membrane interactions. Our NMR experiments showed that human intestinal FABP interacts weakly with both negatively charged and neutral membranes, but it prefers the negatively charged one. Through simultaneous analysis of NMR relaxation in the rotating-frame (R_1ρ), relaxation dispersion, chemical exchange saturation transfer, and dark-state exchange saturation transfer data, we estimated the affinity of the protein to negatively charged NDs, the dissociation rate, and apparent association rate. We further showed that the protein in the ND-bound state adopts a conformation different from the native structure and the second helix is very likely involved in interactions with NDs. We also found a membrane-induced FABP conformational state that exists only in the presence of NDs. This state is native-like, different from other conformational states in structure, unbound to NDs, and in dynamic equilibrium with the ND-bound state.     

The crystal structure of odorant binding protein 7 from Anopheles gambiae exhibits an outstanding adaptability of its binding site

Anopheles gambiae (Agam) targets human and animals by using its olfactory system, leading to the spread of Plasmodium falciparum, the malaria vector. Odorant binding proteins (OBPs) participate to the first event in odorant recognition and constitute an interesting target for insect control. OBPs interact with olfactory receptors to which they deliver the odorant molecule. We have undertaken a large-scale study of proteins belonging to the olfactory system of Agam with in mind of designing strong olfactory repellants. Here, we report the expression, three-dimensional structures and binding properties of AgamOBP07, a member of a new structural class of OBPs, characterized by the occurrence of eight cysteines. We showed that AgamOBP07 possesses seven α-helices and four disulfide bridges, instead of six α-helices and three disulfide bridges in classical OBPs. The extra seventh helix is located at the surface of the protein, locked by the fourth disulfide bridge, and forms a wall of the internal cavity. The binding site of the protein is mainly hydrophobic, elongated and open and is able to accommodate elongated ligands, linear or polycyclic, as suggested also by binding experiments. An elongated electron density was observed in the internal cavity of the purified protein, belonging to a serendipitous ligand. The structure of AgamOBP07 in complex with an azo-bicyclic model compound reveals that a large conformational change in the protein has reshaped its binding site, provoking helix 4 unfolding and doubling of the cavity volume.     

Testing of the portal hypothesis: analysis of a V32G, F57G, K58G mutant of the fatty acid binding protein of the murine adipocyte

The portal region of fatty acid binding proteins is hypothesized to function as a dynamic aperture, controlling accessibility of external ligands to the internal fatty acid binding cavity. To test this hypothesis, a triple mutant of the murine FABP4 has been developed (V32G, F57G, K58G, referred to as the portal mutant) that is predicted to constitutively enlarge the opening due to a reduction in the molecular dimensions of the side chains of key portal amino acids. The portal mutant was purified from expressing Escherichia coli, its stability was evaluated, and the thermodynamics and kinetics of ligand binding were compared to that of wild-type protein. Introduction of the three amino acid substitutions caused no significant change in the stability of the protein with a free energy of unfolding of 13.7 kJ/mol as compared to 14.0 kJ/mol for the wild-type protein. The portal mutant exhibited a modest decrease (4-fold) in ligand binding affinity using the fluorescent probe 1-anilinonaphthalene-8-sulfonic acid (1,8-ANS) as a surrogate ligand. 1,8-ANS displacement assays revealed that the binding affinity for oleate increased from a K0.5 of 196 +/- 15 nM for the wild-type protein to 165 +/- 8 for the portal mutant, while that for arachidonate decreased from the wild type of 186 +/- 11 nM to 418 +/- 26 nM for the portal mutant. To evaluate cavity accessibility, rate of 1,8-ANS binding was assessed between the portal and wild-type protein. Using equimolar amounts of ligand and protein at 4 degrees, 1,8-ANS bound within the cavity to 95% saturation (t0.95) in 750 ms, while the mutant protein was fully modified in less than 1.4 ms. To independently evaluate cavity accessibility, modification of the sole protein cysteine residue, C117 residing within the cavity near C2-C4 of the bound ligand, was monitored using 5,5'-dithio-bis(2-nitrobenzoic acid) modification. The half time for modification (t0.5) for the wild-type protein was approximately 20 s, while that for V32G F57G K58G occurred in less than a second. As such, enlargement of the portal region of FABP4 markedly increased the accessibility of ligands to the cavity while having only modest effects on ligand affinity. Taken together, these data provide support for the portal region hypothesis and suggest dynamic fluctuations in this region regulate cavity access, but not ligand affinity or selectivity.     

Molecular dynamics simulations of the adipocyte lipid binding protein reveal a novel entry site for the ligand

The adipocyte lipid binding protein (ALBP) binds fatty acids (FA) in a cavity that is inaccessible from the bulk. Therefore, the penetration of the FA necessitates conformational changes whose nature is still unknown. It was suggested that the lipid first enters through a "portal region" which consists of the alphaII helix and the adjacent tight turns. The initial event in the ligand binding must be the interaction of the lipid with the protein surface. To analyze this interaction, we have carried out three molecular dynamics simulations of the apo-ALBP, with a palmitate ion initially located at different positions near the protein surface. The simulation indicated that the ligand could adsorb to the protein in more than one location. Yet, in one case, the ligand managed to penetrate the protein by entering a newly formed cavity some 10 A deep. The entry site is located near the N-terminus, at the junction between the loops connecting the beta-strands. Further progression of the penetration seems to be arrested by the need for desolvation of the COOH end of the headgroup. Evolutionary analysis showed that amino acids in this entry site are well conserved. On the basis of these observations, we suggest that the ligand may enter the protein from a site other than the portal region. Furthermore, the rate-limiting step is proposed to be the desolvation of the FA polar headgroup.     

NMR unfolding studies on a liver bile acid binding protein reveal a global two-state unfolding and localized singular behaviors

The folding properties of a bile acid binding protein, belonging to a subfamily of the fatty acid binding proteins, have been here investigated both by hydrogen exchange measurements, using the SOFAST NMR approach, and urea denaturation experiments. The urea unfolding profiles of individual residues, acting as single probes, were simultaneously analyzed through a global fit, according to a two-state unfolding model. The resulting conformational stability ΔG_U(H₂O) = 7.2 ± 0.25 kcal mol⁻¹ is in good agreement with hydrogen exchange stability ΔG_op. While the majority of protein residues satisfy this model, few amino-acids display a singular behavior, not directly amenable to the presence of a folding intermediate, as reported for other fatty acid binding proteins. These residues are part of a protein patch characterized by enhanced plasticity. To explain this singular behavior a tentative model has been proposed which takes into account the interplay between the dynamic features and the formation of transient aggregates. A functional role for this plasticity, related to translocation across the nuclear membrane, is discussed.     

Insight into the interaction sites between fatty acid binding proteins and their ligands

Fatty acid binding proteins (FABPs), are evolutionarily conserved small cytoplasmic proteins that occur in many tissue-specific types. One of their primary functions is to facilitate the clearance of the cytoplasmic matrix from free fatty acids and of other detergent-like compounds. Crystallographic studies of FABP proteins have revealed a well defined binding site located deep inside their β-clam structure that is hardly exposed to the bulk solution. However, NMR measurements revealed that, when the protein is equilibrated with its ligands, residues that are clearly located on the outer surface of the protein do interact with the ligand. To clarify this apparent contradiction we applied molecular dynamics simulations to follow the initial steps associated with the FABP–fatty acid interaction using, as a model, the interaction of toad liver basic FABP, or chicken liver bile acid binding protein, with a physiological concentration of palmitate ions. The simulations (~200 ns of accumulated time) show that fatty acid molecules interact, unevenly, with various loci on the protein surface, with the favored regions being the portal and the anti-portal domains. Random encounters with palmitate at these regions led to lasting adsorption to the surface, while encounters at the outer surface of the β-clam were transient. Therefore, we suggest that the protein surface is capable of sequestering free fatty acids from solution, where brief encounters evolve into adsorbed states, which later mature by migration of the ligand into a more specific binding site.