Molecular dynamics simulations of adipocyte lipid-binding protein: effect of electrostatics and acyl chain unsaturation

A FABP4-PPARγ signaling axis regulates human monocyte responses to electrophilic fatty acid nitroalkenes

Nitro-fatty acids (NO2-FA) are electrophilic lipid mediators derived from unsaturated fatty acid nitration. These species are produced endogenously by metabolic and inflammatory reactions and mediate anti-oxidative and anti-inflammatory responses. NO2-FA have been postulated as partial agonists of the Peroxisome Proliferator-Activated Receptor gamma (PPARγ), which is predominantly expressed in adipocytes and myeloid cells. Herein, we explored molecular and cellular events associated with PPARγ activation by NO2-FA in monocytes and macrophages. NO2-FA induced the expression of two PPARγ reporter genes, Fatty Acid Binding Protein 4 (FABP4) and the scavenger receptor CD36, at early stages of monocyte differentiation into macrophages. These responses were inhibited by the specific PPARγ inhibitor GW9662. Attenuated NO2-FA effects on PPARγ signaling were observed once cells were differentiated into macrophages, with a significant but lower FABP4 upregulation, and no induction of CD36. Using in vitro and in silico approaches, we demonstrated that NO2-FA bind to FABP4. Furthermore, the inhibition of monocyte FA binding by FABP4 diminished NO2-FA-induced upregulation of reporter genes that are transcriptionally regulated by PPARγ, Keap1/Nrf2 and HSF1, indicating that FABP4 inhibition mitigates NO2-FA signaling actions. Overall, our results affirm that NO2-FA activate PPARγ in monocytes and upregulate FABP4 expression, thus promoting a positive amplification loop for the downstream signaling actions of this mediator.

Fatty Acid and Retinol-Binding Protein: Unusual Protein Conformational and Cavity Changes Dictated by Ligand Fluctuations

Lipid-binding proteins (LBPs) are soluble proteins responsible for the uptake, transport, and storage of a large variety of hydrophobic lipophilic molecules including fatty acids, steroids, and other lipids in the cellular environment. Among the LBPs, fatty acid binding proteins (FABPs) present preferential binding affinities for long-chain fatty acids. While most of FABPs in vertebrates and invertebrates present similar β-barrel structures with ligands accommodated in their central cavity, parasitic nematode worms exhibit additional unusual α-helix rich fatty acid- and retinol-binding proteins (FAR). Herein, we report the comparison of extended molecular dynamics (MD) simulations performed on the ligand-free and palmitic acid-bond states of the Necator americanus FAR-1 (Na-FAR-1) with respect to other classical β-barrel FABPs. Principal component analysis (PCA) has been used to identify the different conformations adopted by each system during MD simulations. The α-helix fold encompasses a complex internal ligand-binding cavity with a remarkable conformational plasticity that allows reversible switching between distinct states in the holo-Na-FAR-1. The cavity can change up to one-third of its size affected by conformational changes of the protein-ligand complex. Besides, the ligand inside the cavity is not fixed but experiences large conformational changes between bent and stretched conformations. These changes in the ligand conformation follow changes in the cavity size dictated by the transient protein conformation. On the contrary, protein-ligand complex in β-barrel FABPs fluctuates around a unique conformation. The significantly more flexible holo-Na-FAR-1 ligand-cavity explains its larger ligand multiplicity respect to β-barrel FABPs.

Atomistic Insights into the Functional Instability of the Second Helix of Fatty Acid Binding Protein

Structural dynamics of fatty acid binding proteins (FABPs), which accommodate poorly soluble ligands in the internalized binding cavities, are intimately related to their function. Recently, local unfolding of the α-helical cap in a variant of human intestinal FABP (IFABP) has been shown to correlate with the kinetics of ligand association, shedding light on the nature of the critical conformational reorganization. Yet, the physical origin and mechanism of the functionally relevant transient unfolding remain elusive. Here, we investigate the intrinsic structural instability of the second helix (αII) of IFABP in comparison with other segments of the protein using hydrogen-exchange NMR spectroscopy, microsecond molecular dynamics simulations, and enhanced sampling techniques. Although tertiary interactions positively contribute to the stability of helices in IFABP, the intrinsic unfolding tendency of αII is encoded in its primary sequence and can be described by the Lifson-Roig theory in the absence of tertiary interactions. The unfolding pathway of αII in intact proteins involves an on-pathway intermediate state that is characterized with the fraying of the last helical turn, captured by independent enhanced sampling methods. The simulations in this work, combined with hydrogen-exchange NMR data, provide new, to our knowledge, atomistic insights into the functional local unfolding of FABPs.

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.

Structural and functional characterization of human peripheral nervous system myelin protein P2

The myelin sheath is a tightly packed multilayered membrane structure insulating selected axons in the central and the peripheral nervous systems. Myelin is a biochemically unique membrane, containing a specific set of proteins. In this study, we expressed and purified recombinant human myelin P2 protein and determined its crystal structure to a resolution of 1.85 A. A fatty acid molecule, modeled as palmitate based on the electron density, was bound inside the barrel-shaped protein. Solution studies using synchrotron radiation indicate that the crystal structure is similar to the structure of the protein in solution. Docking experiments using the high-resolution crystal structure identified cholesterol, one of the most abundant lipids in myelin, as a possible ligand for P2, a hypothesis that was proven by fluorescence spectroscopy. In addition, electrostatic potential surface calculations supported a structural role for P2 inside the myelin membrane. The potential membrane-binding properties of P2 and a peptide derived from its N terminus were studied. Our results provide an enhanced view into the structure and function of the P2 protein from human myelin, which is able to bind both monomeric lipids inside its cavity and membrane surfaces.

Modeling fatty acid delivery from intestinal fatty acid binding protein to a membrane

Intestinal fatty acid binding protein (IFABP) interacts with biological membranes and delivers fatty acid (FA) into them via a collisional mechanism. However, the membrane-bound structure of the protein and the pathway of FA transfer are not precisely known. We used molecular dynamics (MD) simulations with an implicit membrane model to determine the optimal orientation of apo- and holo-IFABP (bound with palmitate) on an anionic membrane. In this orientation, the helical portal region, delimited by the alphaII helix and the betaC-betaD and betaE-betaF turns, is oriented toward the membrane whereas the putative beta-strand portal, delimited by the betaB-betaC, betaF-betaG, betaH-betaI turns and the N terminus, is exposed to solvent. Starting from the MD structure of holo-IFABP in the optimal orientation relative to the membrane, we examined the release of palmitate via both pathways. Although the domains can widen enough to allow the passage of palmitate, fatty acid release through the helical portal region incurs smaller conformational changes and a lower energetic cost.

Glycoimmunomics of human cancer: current concepts and future perspectives

Future strategies for the treatment of human cancer require a full appreciation of the intracellular and extracellular changes that accompany neoplastic transformation. The changes may involve a variety of micro- and macro-molecules, including, but not restricted to, peptides, proteins (with sugar and/or lipid moieties), oligosaccharides, glycolipids (neutral or acidic, e.g., gangliosides), ceramides, fatty acids and other lipids. Although several therapeutic approaches have been well developed in recent years, most of the reported studies focus on proteins and peptides. Glycoantigens and lipoantigens have been neglected. Elucidation of the profiles and properties of all molecules associated with tumor progression is required to develop a successful strategy to treat human cancer. This review describes the unique immunomics of tumor-associated glycoantigens and explains why the field of glycoimmunomics may yield clinically important biomarkers and treatments for the management of human cancer.

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.

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.

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.

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

Molecular dynamics simulations of two structurally similar fatty acid-binding proteins interacting with stearic acid are described. The calculations relate to recent ligand binding measurements and suggest similarities and differences between the two systems. Charged and neutral forms of the fatty acid were examined. The charged forms led to rapid trajectory divergence, whereas the protonated forms remained stable over the length of their 1-ns production trajectories. The two protein systems showed similar sets of total interaction energies with the ligand. However, the strengths of individual amino acids interacting with the ligand differ. Furthermore, covariance analysis of the ligand with both protein and water suggests that the stearic acid in the adipocyte fatty acid-binding protein is coupled more strongly to the water than to the protein. The stearic acid in the muscle fatty acid-binding protein is seen to be coupled differentially along the length of the chain to the protein. These differences could help to rationalize the stronger binding affinity for stearic acid in the human muscle fatty acid-binding protein. An importance scale, based on both covariance and interaction energy with the ligand, is proposed to identify residues that may be important for binding function.

Ligand specificity of brain lipid-binding protein

Brain lipid-binding protein (BLBP) is a member of the fatty acid-binding protein (FABP) family. Although BLBP expression in the developing central nervous system is complex, a close correlation between its expression and radial glial differentiation has been observed. Furthermore, antibodies to BLBP can block glial cell differentiation in mixed primary cell cultures. Here we describe the ligand binding properties of BLBP. The binding affinities of BLBP for oleic acid (Kd approximately 0.44 microM) and arachidonic acid (Kd approximately 0.25 microM) are similar to those reported for other FABPs, but BLBP does not bind to palmitic acid or arachidinic acid. These and other experiments establish that BLBP has a strong preference for binding long chain polyunsaturated fatty acids. A probable in vivo ligand for BLBP is docosahexaenoic acid (DHA), since its binding affinity (Kd approximately 10 nM) is the highest yet reported for an FABP/ligand interaction, exceeding even the affinity of retinoic acid for its binding proteins. Furthermore, the requirement of DHA for nervous system development and the coincident expression of BLBP during these developmental stages suggest that the physiologic role of BLBP may involve DHA utilization. Finally, we present a structural model of BLBP/DHA interaction that provides insight into both the structural characteristics important for ligand binding and the effects of specific mutations upon BLBP/ligand interactions.