Simulations of fatty acid-binding proteins suggest sites important for function. I. Stearic acid

Concerted dynamic motions of an FABP4 model and its ligands revealed by microsecond molecular dynamics simulations

In this work, we investigate the dynamic motions of fatty acid binding protein 4 (FABP4) in the absence and presence of a ligand by explicitly solvated all-atom molecular dynamics simulations. The dynamics of one ligand-free FABP4 and four ligand-bound FABP4s is compared via multiple 1.2 μs simulations. In our simulations, the protein interconverts between the open and closed states. Ligand-free FABP4 prefers the closed state, whereas ligand binding induces a conformational transition to the open state. Coupled with opening and closing of FABP4, the ligand adopts distinct binding modes, which are identified and compared with crystal structures. The concerted dynamics of protein and ligand suggests that there may exist multiple FABP4–ligand binding conformations. Thus, this work provides details about how ligand binding affects the conformational preference of FABP4 and how ligand binding is coupled with a conformational change of FABP4 at an atomic level.

Dissection of the critical binding determinants of cellular retinoic acid binding protein II by mutagenesis and fluorescence binding assay

The binding of retinoic acid to mutants of Cellular Retinoic Acid Binding Protein II (CRABPII) was evaluated to better understand the importance of the direct protein/ligand interactions. The important role of Arg111 for the correct structure and function of the protein was verified and other residues that directly affect retinoic acid binding have been identified. Furthermore, retinoic acid binding to CRABPII mutants that lack all previously identified interacting amino acids was rescued by providing a carboxylic acid dimer partner in the form of a Glu residue.

Molecular dynamics simulations of palmitate entry into the hydrophobic pocket of the fatty acid binding protein

The entry of substrate into the active site is the first event in any enzymatic reaction. However, due to the short time interval between the encounter and the formation of the stable complex, the detailed steps are experimentally unobserved. In the present study, we report a molecular dynamics simulation of the encounter between palmitate molecule and the Toad Liver fatty acid binding protein, ending with the formation of a stable complex resemblance in structure of other proteins of this family. The forces operating on the system leading to the formation of the tight complex are discussed.

Release of albumin from oligoester plastic matrices: effect of magnesium oxide and bivalent stearates

Biodegradable implantable matrices containing bovine serum albumin were prepared from oligoesters by melting, and subsequently tested on in vitro albumin release. The linear poly (DL-lactic acid) and the branched terpolymer of DL-lactic acid, glycolic acid, and mannitol were synthesized. Products were of similar molecular weight and possessed different thermal and swelling characteristics. Oligoesters were loaded with 4% albumin and plasticized by 30% triacetin. Other additives added into the matrices as albumin stabilizers were divalent stearates and magnesium oxide. The influences of oligomer molecules constitution, divalent ion stearates or magnesium oxide addition, and triacetin concentration on the albumin release were quantified. SDS-PAGE revealed protein hydrolysis during the dissolution tests.

The pH-dependent unfolding mechanism of P2 myelin protein: an experimental and computational study

The P2 protein is a small, extrinsic protein of the myelin membrane in the peripheral nervous system that structurally belongs to the fatty acid binding proteins (FABPs) family, sharing with them a 10 strands beta-barrel structure. FABPs appear to be involved in cellular fatty acid transport, but very little is known about the role of P2 in the metabolism of peripheral myelin lipids. Study of protein conformation at different pHs is a useful tool for the characterization of the unfolding mechanisms and the intrinsic conformational properties of the protein, and may give insight into factors that guide protein folding pathways. In particular, low pH conditions have been shown to induce partially folded states in several proteins. In this paper, the acidic unfolding of purified P2 protein was studied with both spectroscopic techniques and molecular dynamics simulation. Both experimental and computational results indicate the presence of a partly folded state at low pH, which shows structural changes mainly involving the lid that is formed by the helix-turn-helix domain. The opening of the lid, together with a barrel relaxation, could regulate the ligand exchanges near the cell membrane, supporting the hypothesis that the P2 protein may transport fatty acids between Schwann cells and peripheral myelin.

Fatty acid binding proteins: same structure but different binding mechanisms? Molecular dynamics simulations of intestinal fatty acid binding protein

Fatty acid binding proteins (FABPs) carry fatty acids (FAs) and other lipids in the cellular environment, and are thus involved in processes such as FA uptake, transport, and oxidation. These proteins bind either one or two ligands in a binding site, which appears to be inaccessible from the bulk. Thus, the entry of the substrate necessitates a conformational change, whose nature is still unknown. A possible description of the ligand binding process is given by the portal hypothesis, which suggests that the FA enters the protein through a dynamic area known as the portal region. On the other hand, recent simulations of the adipocyte lipid binding protein (ALBP) suggested a different entry site (the alternative portal). In this article, we discuss molecular dynamics simulations of the apo-intestinal-FABP (I-FABP) in the presence of palmitate molecule(s) in the simulation box. The simulations were carried out to study whether the FA can enter the protein during the simulations (as in the ALBP) and where the ligand entry site is (the portal region, the alternative portal or a different domain). The analysis of the simulations revealed a clear difference between the ALBP and the I-FABP. In the latter case, the palmitate preferentially adsorbed to the portal region, which was more mobile than the rest of the protein. However, no ligand entry was observed in the multi-nanosecond-long simulations, in contrast to ALBP. These findings suggest that, although the main structural motif of the FABPs is common, the fine details of each individual protein structure grossly modulate its reactivity.

Water dynamics in the large cavity of three lipid-binding proteins monitored by (17)O magnetic relaxation dispersion

Intracellular lipid-binding proteins contain a large binding cavity filled with water molecules. The role played by these water molecules in ligand binding is not well understood, but their energetic and dynamic properties must be important for protein function. Here, we use the magnetic relaxation dispersion (MRD) of the water 17O resonance to investigate the water molecules in the binding cavity of three different lipid-binding proteins: heart fatty acid-binding protein (H-FABP), ileal lipid-binding protein (I-LBP) and intestinal fatty acid-binding protein (I-FABP). Whereas about half of the crystallographically visible water molecules appear to be expelled by the ligand, we find that ligand binding actually increases the number of water molecules within the cavity. At 300 K, the water molecules in the cavity exchange positions on a time-scale of about 1ns and exchange with external water on longer time-scales (0.01-1 micros). Exchange of water molecules among hydration sites within the cavity should be strongly coupled to ligand motion. Whereas a recent MD simulation indicates that the structure of the cavity water resembles a bulk water droplet, the present MRD results show that its dynamics is more than two orders of magnitude slower than in the bulk. These findings may have significant implications for the strength, specificity and kinetics of lipid binding.

The crystal structure of Echinococcus granulosus fatty-acid-binding protein 1

We describe the 1.6 A crystal structure of the fatty-acid-binding protein EgFABP1 from the parasitic platyhelminth Echinococcus granulosus. E. granulosus causes hydatid disease, which is a major zoonosis. EgFABP1 has been implicated in the acquisition, storage, and transport of lipids, and may be important to the organism since it is incapable of synthesising most of its lipids de novo. Moreover, EgFABP1 is a promising candidate for a vaccine against hydatid disease. The crystal structure reveals that EgFABP1 has the expected 10-stranded beta-barrel fold typical of the family of intracellular lipid-binding proteins, and that it is structurally most similar to P2 myelin protein. We describe the comparison of the crystal structure of EgFABP1 with these proteins and with an older homology model for EgFABP1. The electron density reveals the presence of a bound ligand inside the cavity, which we have interpreted as palmitic acid. The carboxylate group of the fatty acid interacts with the protein's P2 motif, consisting of a conserved triad R em leader R-x-Y. The hydrophobic tail of the ligand assumes a fairly flat, U-shaped conformation and has relatively few interactions with the protein.We discuss some of the structural implications of the crystal structure of EgFABP1 for related platyhelminthic FABPs.

Simulations of apo and holo-fatty acid binding protein: structure and dynamics of protein, ligand and internal water

Two molecular dynamics simulations of 5 ns each have been carried out for rat intestinal fatty acid binding protein, in apo-form and with bound palmitate. The fatty acid and a number of water molecules are encapsulated in a large interior cavity of the barrel-shaped protein. The simulations are compared to experimental data and analyzed in terms of root mean square deviations, atomic B-factors, secondary structure elements, hydrogen bond patterns, and distance constraints derived from nuclear Overhauser experiments. Excellent agreement is found between simulated and experimental solution structures of holo-FABP, but a number of differences are observed for the apo-form. The ligand in holo-FABP shows considerable displacement after about 1.5 ns and displays significant configurational entropy. A novel computational approach has been employed to identify internal water and to capture exchange pathways. Orifices in the portal and gap regions of the protein, discussed in the experimental literature, have been confirmed as major openings for solvent exchange between the internal cavity and bulk water. A third opening on the opposite side of the barrel experiences significant exchange but it does not provide a pathway for further passage to the central cavity. Internal water is characterized in terms of density distributions, interaction energies, mobility, protein contact times, and water molecule coordination. A number of differences are observed between the apo and holo-forms and related to differences in the protein structure. Solvent inside apo-FABP, for example, shows characteristics of a water droplet, while solvent in holo-FABP benefits from interactions with the ligand headgroup and slightly stronger interactions with protein residues.

Bovine beta-lactoglobulin: interaction studies with palmitic acid

Bovine beta-lactoglobulin (BLG) in vivo has been found complexed with fatty acids, especially palmitic and oleic acid. To elucidate the still unknown structure-function relationship in this protein, the interactions between 13C enriched palmitic acid (PA) and BLG were investigated by means of one-, two-, and three-dimensional NMR spectroscopy in the pH range 8.4-2.1. The NMR spectra revealed that at neutral pH the ligand is bound within the central cavity of BLG, with the methyl end deeply buried within the protein. The analysis of 13C spectra of the holo protein revealed the presence of conformational variability of bound PA carboxyl end in the pH range 8.4-5.9, related to the Tanford transition. The release of PA starts at pH lower than 6.0, and it is nearly complete at acidic pH. This finding is relevant in relation to the widely reported hypothesis that this protein can act as a transporter through the acidic gastric tract. Ligand binding and release is shown to be completely reversible over the entire pH range examined, differently from other fatty acid binding proteins whose behavior is analyzed throughout the paper. The mode of interaction of BLG is compatible with the proposed function of facilitating the digestion of milk fat during the neonatal period of calves.

Differences between apo and three holo forms of the intestinal fatty acid binding protein seen by molecular dynamics computer calculations

It is commonly believed that binding affinity can be estimated by consideration of local changes of ligand and protein. This paper discusses a set of molecular dynamics simulations of intestinal fatty acid binding protein addressing the protein's response to presence or absence of different ligands. A 5-ns simulation was performed of the protein without a ligand, and three simulations (one 5-ns and two 2-ns) were performed with different fatty acids bound. The results indicate that, although the basic protein structure is unchanged by the presence of the ligand, other properties are significantly affected by ligand binding. For example, zero-time covariance patterns between protein, bound waters, and ligand vary between the different simulations. Moreover, the interaction energies between ligand and specific residues indicate that different ligands are stabilized in different ways. In sum, the results suggest that binding thermodynamics within this system will need to be calculated not from a subset of nearby protein:ligand interactions, but will depend on a knowledge of the motions coupling together water, protein, and ligand.

Muscle fatty acid-binding protein

Muscle or heart fatty acid-binding protein is a low molecular weight protein that binds long-chain fatty acids in the cytosol of muscle tissues. The three-dimensional structure of the human, bovine and insect proteins are known, either via X-ray or NMR techniques. The folding of the protein closely resembles that of the other FABPs: ten anti-parallel beta-strands are arranged to form a clam shell, closed at one end by two alpha-helices. This arrangement allows the formation of an internal cavity where the fatty acid can be accommodated, protected and isolated from the external environment. The fatty acid in the protein interior is stabilized by electrostatic and hydrogen bond interactions of its carboxylic head with charged or polar residues of the protein and by interactions of its tail with hydrophobic residues. The three-dimensional structure of different fatty acid-protein complexes along with molecular dynamics simulations are now providing insight into the molecular details of the specificity of the ligand binding.

Structural properties of the adipocyte lipid binding protein

The adipocyte lipid binding protein, ALBP (also adipocyte fatty acid binding protein, A-FABP, 422 protein, aP2, and p15 protein), is one of the most studied of the intracellular lipid binding protein family. Here we sequentially compare the different sources of ALBP and describe the idea that one-third of the amino acid side chains near the N-terminal end appear to play a major role in conformational dynamics and in ligand transfer. Crystallographic data for mouse ALBP are summarized and the ligand binding cavity analyzed in terms of the overall surface and conformational dynamics. The region of the proposed ligand portal is described. Amino acid side chains critical to cavity formation and fatty acid interactions are analyzed by comparing known crystal structures containing a series of different hydrophobic ligands. Finally, we address ALBP ligand binding affinity and thermodynamic studies.

Structure and dynamics of the fatty acid binding cavity in apo rat intestinal fatty acid binding protein

The structure and dynamics of the fatty acid binding cavity in I-FABP (rat intestinal fatty acid binding protein) were analyzed. In the crystal structure of apo I-FABP, the probe occupied cavity volume and surface are 539+/-8 A3 and 428 A2, respectively (1.4 A probe). A total of 31 residues contact the cavity with their side chains. The side-chain cavity surface is partitioned according to the residue type as follows: 36-39% hydrophobic, 21-25% hydrophilic, and 37-43% neutral or ambivalent. Thus, the cavity surface is neither like a typical protein interior core, nor is like a typical protein external surface. All hydrophilic residues that contact the cavity-with the exception of Asp74-are clustered on the one side of the cavity. The cavity appears to expand its hydrophobic surface upon fatty acid binding on the side opposite to this hydrophilic patch. In holo I-FABP the fatty acid chain interactions with the hydrophilic side chains are mediated by water molecules. Molecular dynamics (MD) simulation of fully solvated apo I-FABP showed global conformational changes of I-FABP, which resulted in a large, but seemingly transient, exposure of the cavity to the external solvent. The packing density of the side chains lining the cavity, studied by Voronoi volumes, showed the presence of two distinctive small hydrophobic cores. The MD simulation predicts significant structural perturbations of the cavity on the subnanosecond time scale, which are capable of facilitating exchange of I-FABP internal water.

Water molecules in the binding cavity of intestinal fatty acid binding protein: dynamic characterization by water 17O and 2H magnetic relaxation dispersion

The hydration of intestinal fatty acid binding protein (IFABP) in apo-form and complexed with palmitate, oleate, and 1-anilino-8-naphthalene sulfonate (ANS) has been studied by water 17O and 2H magnetic relaxation dispersion (MRD) measurements. These ligands bind in a large internal cavity, displacing most of the crystallographically identified cavity water molecules. Unlike most other proteins, IFABP gives rise to MRD profiles with two dispersion steps. The low-frequency dispersion yields a correlation time of 7 ns at 300 K, matching the known tumbling time of IFABP. The dispersion amplitude requires only three (apo) or four (holo) long-lived and ordered water molecules (residence time 0.01-4 microseconds at 300 K). Comparison of MRD profiles from the different complexes indicates that the displaced cavity water molecules are short-lived. The few long-lived (>10 ns) water molecules required by the MRD data are tentatively assigned to crystallographic hydration sites on the basis of accessibility, positional order, and H-bonding. The amplitude of the high-frequency dispersion corresponds to 10-20 moderately ordered water molecules, with a correlation time of ca. 1 ns that may reflect a transient opening of the cavity required for exchange with external water.

The third leg: molecular dynamics simulations of lipid binding proteins

Molecular dynamics computer simulations can provide a third leg which balances the contributions of both structural biology and binding studies performed on the lipid binding protein family. In this context, these calculations help to establish a dialogue between all three communities, by relating experimental observables with details of structure. Working towards this connection is important, since experience has shown the difficulty of inferring thermodynamic properties from a single static conformation. The challenge is exemplified by ongoing attempts to interpret the impact of mutagenesis on structure and function (i.e. binding). A detailed atomic-level understanding of this system could be achieved with the support of all three legs, paving the way towards rational design of proteins with novel specificities. This paper provides an outline of the connections possible between experiment and theory concerning lipid binding proteins.

Simulations of fatty acid-binding proteins. II. Sites for discrimination of monounsaturated ligands

Fatty acid binding proteins (FABPs) can discriminate between saturated and unsaturated fatty acids via molecular mechanisms that are not understood. Molecular dynamics computer calculations are used to suggest the relationship between tertiary structure and binding specificity. Three separate 1-ns simulations, with explicit solvent, are presented: 1) oleic acid (C18:1 cis) bound to adipocyte FABP, 2) oleic acid bound to human muscle FABP, and 3) elaidic acid (C18:1 trans) bound to human muscle FABP. The average structural, dynamic, and energetic properties of the trajectory were analyzed, as were the motional correlations. The molecular dynamics trajectories reveal intriguing differences among all three systems. For example, the two proteins have different strengths of interaction energy with the ligand and different motional coupling, as seen with covariance analysis. This suggests distinctive molecular behavior of monounsaturated fatty acids in the two similar proteins. An importance scale, based on motional correlation and interaction energy between protein and ligand, is proposed, to help identify amino acids involved with the discrimination of ligand saturation state or geometric isomerization.