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

Lipid exchange in crystal-confined fatty acid binding proteins: X-ray evidence and molecular dynamics explanation

Fatty acid binding proteins (FABPs) are responsible for the long-chain fatty acids (FAs) transport inside the cell. However, despite the years, since their structure is known and the many studies published, there is no definitive answer about the stages of the lipid entry-exit mechanism. Their structure forms a β -barrel of 10 anti-parallel strands with a cap in a helix-turn-helix motif, and there is some consensus on the role of the so-called portal region, involving the second α -helix from the cap ( α 2), β C- β D, and β E- β F turns in FAs exchange. To test the idea of a lid that opens, we performed a soaking experiment on an h-FABP crystal in which the cap is part of the packing contacts, and its movement is strongly restricted. Even in these conditions, we observed the replacement of palmitic acid by 2-Bromohexadecanoic acid (Br-palmitic acid). Our MD simulations reveal a two-step lipid entry process: (i) The travel of the lipid head through the cavity in the order of tens of nanoseconds, and (ii) The accommodation of its hydrophobic tail in hundreds to thousands of nanoseconds. We observed this even in the cases in which the FAs enter the cavity by their tail. During this process, the FAs do not follow a single trajectory, but multiple ones through which they get into the protein cavity. Thanks to the complementary views between experiment and simulation, we can give an approach to a mechanistic view of the exchange process.

Molecular dynamics simulation of the heart type fatty acid binding protein in a crystal environment

Crystallographic data comes from a space-time average over all the unit cells within the crystal, so dynamic phenomena do not contribute significantly to the diffraction data. Many efforts have been made to reconstitute the movement of the macromolecules and explore the microstates that the confined proteins can adopt in the crystalline network. We explored different strategies to simulate a heart fatty acid binding protein (H-FABP) crystal by means of Molecular Dynamics (MD) simulations. We evaluate the effect of introducing restraints according to experimental isotropic B-factors and we analyzed the H-FABP motions in the crystal using Principal Component Analysis (PCA), isotropic and anisotropic B-factors. We compared the behavior of the protein simulated in the crystal confinement versus in solution, and we observed the effect of that confinement in the mobility of the protein residues. Restraining one-third of Cα atoms based on experimental B-factors produce lower B-factors than simulations without restraints, showing that the position restraint of the atoms with the lowest experimental B-factor is a good strategy to maintain the geometry of the crystal with an obvious decrease in the degrees of motion of the protein. PCA shows that, as position restraint reduces the conformational space explored by the system, the motion of the crystal is better recovered, for an essential subspace of the same size, in the simulations without restraints. Restraining only one Cα seems to be a good balance between giving flexibility to the system and preserving its structure. Communicated by Ramaswamy H. Sarma.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

Structure and dynamics of a human myelin protein P2 portal region mutant indicate opening of the β barrel in fatty acid binding proteins

Background

Myelin is a multilayered proteolipid sheath wrapped around selected axons in the nervous system. Its constituent proteins play major roles in forming of the highly regular membrane structure. P2 is a myelin-specific protein of the fatty acid binding protein (FABP) superfamily, which is able to stack lipid bilayers together, and it is a target for mutations in the human inherited neuropathy Charcot-Marie-Tooth disease. A conserved residue that has been proposed to participate in membrane and fatty acid binding and conformational changes in FABPs is Phe57. This residue is thought to be a gatekeeper for the opening of the portal region upon ligand entry and egress.

Results

We performed a structural characterization of the F57A mutant of human P2. The mutant protein was crystallized in three crystal forms, all of which showed changes in the portal region and helix α2. In addition, the behaviour of the mutant protein upon lipid bilayer binding suggested more unfolding than previously observed for wild-type P2. On the other hand, membrane binding rendered F57A heat-stable, similarly to wild-type P2. Atomistic molecular dynamics simulations showed opening of the side of the discontinuous β barrel, giving important indications on the mechanism of portal region opening and ligand entry into FABPs. The results suggest a central role for Phe57 in regulating the opening of the portal region in human P2 and other FABPs, and the F57A mutation disturbs dynamic cross-correlation networks in the portal region of P2.

Conclusions

Overall, the F57A variant presents similar properties to the P2 patient mutations recently linked to Charcot-Marie-Tooth disease. Our results identify Phe57 as a residue regulating conformational changes that may accompany membrane surface binding and ligand exchange in P2 and other FABPs.

Electronic supplementary material

The online version of this article (10.1186/s12900-018-0087-2) contains supplementary material, which is available to authorized users.

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.

Fatty acid binding into the highest affinity site of human serum albumin observed in molecular dynamics simulation

Multiple molecular dynamics simulations were performed to investigate the association of stearic acid into the highest affinity binding site of human serum albumin. All binding events ended with a rapid (<10 ps) lock-in of the fatty acid due to formation of a hydrogen bond with Tyr401. The kinetics and energetics of the penetration process both depended linearly on the positional shift of the fatty acid, with an average insertion time and free energy reduction of, respectively, 32 ± 20 ps and 0.70 ± 0.15 kcal/mol per methylene group absorbed. Binding events of longer duration (tbind>1 ns) were characterized by a slow exploration of the pocket entry and, frequently, of a nearby protein crevice corresponding to a metastable state along the route to the binding site. Taken all together, these findings reconstruct the following pathway for the binding process of stearic acid: (i) contact with the protein surface, possibly facilitated by the presence of an intermediate location, (ii) probing of the site entry, (iii) insertion into the protein, and (iv) lock-in at the final position. This general description may also apply to other long-chain fatty acids binding into any of the high-affinity sites of albumin, or to specific sites of other lipid-binding proteins.

The simulation approach to lipid-protein interactions

The interactions between lipids and proteins are crucial for a range of biological processes, from the folding and stability of membrane proteins to signaling and metabolism facilitated by lipid-binding proteins. However, high-resolution structural details concerning functional lipid/protein interactions are scarce due to barriers in both experimental isolation of native lipid-bound complexes and subsequent biophysical characterization. The molecular dynamics (MD) simulation approach provides a means to complement available structural data, yielding dynamic, structural, and thermodynamic data for a protein embedded within a physiologically realistic, modelled lipid environment. In this chapter, we provide a guide to current methods for setting up and running simulations of membrane proteins and soluble, lipid-binding proteins, using standard atomistically detailed representations, as well as simplified, coarse-grained models. In addition, we outline recent studies that illustrate the power of the simulation approach in the context of biologically relevant lipid/protein interactions.

Comparative study of the fatty acid binding process of a new FABP from Cherax quadricarinatus by fluorescence intensity, lifetime and anisotropy

Fatty acid-binding proteins (FABPs) are small cytosolic proteins, largely distributed in invertebrates and vertebrates, which accomplish uptake and intracellular transport of hydrophobic ligands such as fatty acids. Although long chain fatty acids play multiple crucial roles in cellular functions (structural, energy metabolism, regulation of gene expression), the precise functions of FABPs, especially those of invertebrate species, remain elusive. Here, we have identified and characterized a novel FABP family member, Cq-FABP, from the hepatopancreas of red claw crayfish Cherax quadricarinatus. We report the characterization of fatty acid-binding affinity of Cq-FABP by four different competitive fluorescence-based assays. In the two first approaches, the fluorescent probe 8-Anilino-1-naphthalenesulfonate (ANS), a binder of internal cavities of protein, was used either by directly monitoring its fluorescence emission or by monitoring the fluorescence resonance energy transfer occurring between the single tryptophan residue of Cq-FABP and ANS. The third and the fourth approaches were based on the measurement of the fluorescence emission intensity of the naturally fluorescent cis-parinaric acid probe or the steady-state fluorescence anisotropy measurements of a fluorescently labeled fatty acid (BODIPY-C16), respectively. The four methodologies displayed consistent equilibrium constants for a given fatty acid but were not equivalent in terms of analysis. Indeed, the two first methods were complicated by the existence of non specific binding modes of ANS while BODIPY-C16 and cis-parinaric acid specifically targeted the fatty acid binding site. We found a relationship between the affinity and the length of the carbon chain, with the highest affinity obtained for the shortest fatty acid, suggesting that steric effects primarily influence the interaction of fatty acids in the binding cavity of Cq-FABP. Moreover, our results show that the binding affinities of several fatty acids closely parallel their prevalences in the hepatopancreas of C. quadricarinatus as measured under specific diet conditions.

Practical considerations for building GROMOS-compatible small-molecule topologies

Molecular dynamics simulations are being applied to increasingly complex systems, including those involving small endogenous compounds and drug molecules. In order to obtain meaningful and accurate data from these simulations, high-quality topologies for small molecules must be generated in a manner that is consistent with the derivation of the force field applied to the system. Often, force fields are designed for use with macromolecules such as proteins, making their transferability to other species challenging. Investigators are increasingly attracted to automated topology generation programs, although the quality of the resulting topologies remains unknown. Here we assess the applicability of the popular PRODRG server that generates small-molecule topologies for use with the GROMOS family of force fields. We find that PRODRG does not reproduce topologies for even the most well-characterized species in the force field due to inconsistent charges and charge groups. We assessed the effects of PRODRG-derived charges on several systems: pure liquids, amino acids at a hydrophobic-hydrophilic interface, and an enzyme-cofactor complex. We found that partial atomic charges generated by PRODRG are largely incompatible with GROMOS force fields, and the behavior of these systems deviates substantially from that of simulations using GROMOS parameters. We conclude by proposing several points as "best practices" for parametrization of small molecules under the GROMOS force fields.