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

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.

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.

Structural and Dynamic Characterizations Highlight the Deleterious Role of SULT1A1 R213H Polymorphism in Substrate Binding

Sulfotransferase 1A1 (SULT1A1) is responsible for catalyzing various types of endogenous and exogenous compounds. Accumulating data indicates that the polymorphism rs9282861 (R213H) is responsible for inefficient enzymatic activity and associated with cancer progression. To characterize the detailed functional consequences of this mutation behind the loss-of-function of SULT1A1, the present study deployed molecular dynamics simulation to get insights into changes in the conformation and binding energy. The dynamics scenario of SULT1A1 in both wild and mutated types as well as with and without ligand showed that R213H induced local conformational changes, especially in the substrate-binding loop rather than impairing overall stability of the protein structure. The higher conformational changes were observed in the loop3 (residues, 235-263), turning loop conformation to A-helix and B-bridge, which ultimately disrupted the plasticity of the active site. This alteration reduced the binding site volume and hydrophobicity to decrease the binding affinity of the enzyme to substrates, which was highlighted by the MM-PBSA binding energy analysis. These findings highlight the key insights of structural consequences caused by R213H mutation, which would enrich the understanding regarding the role of SULT1A1 mutation in cancer development and also xenobiotics management to individuals in the different treatment stages.

Structural insights into ligand binding features of dual FABP4/5 inhibitors by molecular dynamics simulations

The fatty acid binding protein (FABP) 4 and 5 have been considered as potential targets for the treatment of metabolic diseases. A compensatory upregulation of FABP5 due to the gene ablation of FABP4 in adipocytes indicated the importance of dual FABP4/5 inhibitors. A few compounds have been discovered as dual FABP4/5 inhibitors. However, none exhibited equivalent inhibitory activity against both FABP4 and FABP5, and almost all compounds showed weaker inhibition against FABP5. To provide a better structural understanding for the design of potent dual FABP4/5 inhibitors, molecular dynamics simulations have been performed for 100 ns to disclose the ligand binding features in FABP4 and FABP5 using Amber14, respectively. Key residues were identified by analysis of close contact, hydrogen bond occupancy, binding free energy and alanine scanning mutagenesis. In addition, induced-fit effects have been observed upon ligand binding in the process of simulations. The shifted alkyl chain of ligand in FABP4 was significantly different from that in FABP5 due to the corresponding residues (Phe58^FABP4 and Leu60^FABP5). Thus, to avoid different steric effects made by these two residues, hydrophobic groups of suitable size should be taken into account. Besides, electrostatic and steric effects with Arg107^FABP4 and Arg109^FABP5 should be paid more attention to. The results will facilitate the rational design of dual FABP4/5 inhibitors.

Molecular Dynamics Exploration of Selectivity of Dual Inhibitors 5M7, 65X, and 65Z toward Fatty Acid Binding Proteins 4 and 5

Designing highly selective inhibitors of fatty acid binding proteins 4 and 5 (FABP4 and FABP5) is of importance for treatment of some diseases related with inflammation, metabolism, and tumor growth. In this study, molecular dynamics (MD) simulations combined with molecular mechanics generalized Born surface area (MM-GBSA) method were performed to probe binding selectivity of three inhibitors (5M7, 65X, and 65Z) to FABP4/FABP5 with K_i values of 0.022/0.50 μM, 0.011/0.086 μM, and 0.016/0.12 μM, respectively. The results not only suggest that all inhibitors associate more tightly with FABP4 than FABP5, but also prove that the main forces driving the selective bindings of inhibitors to FABP4 and FABP5 stem from the difference in the van der Waals interactions and polar interactions of inhibitors with two proteins. Meanwhile, a residue-based free energy decomposition method was applied to reveal molecular basis that inhibitors selectively interact with individual residues of two different proteins. The calculated results show that the binding difference of inhibitors to the residues (Phe16, Phe19), (Ala33, Gly36), (Phe57, Leu60), (Ala75, Ala78), (Arg126, Arg129), and (Tyr128, Tyr131) in (FABP4, FABP5) drive the selectivity of inhibitors toward FABP4 and FABP5. This study will provide great help for further design of effective drugs to protect against a series of metabolic diseases, arteriosclerosis, and inflammation.

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.

Exploration of gated ligand binding recognizes an allosteric site for blocking FABP4-protein interaction

Fatty acid binding protein 4 (FABP4), reversibly binding to fatty acids and other lipids with high affinities, is a potential target for treatment of cancers. The binding site of FABP4 is buried in an interior cavity and thereby ligand binding/unbinding is coupled with opening/closing of FABP4. It is a difficult task both experimentally and computationally to illuminate the entry or exit pathway, especially with the conformational gating. In this report we combine extensive computer simulations, clustering analysis, and the Markov state model to investigate the binding mechanism of FABP4 and troglitazone. Our simulations capture spontaneous binding and unbinding events as well as the conformational transition of FABP4 between the open and closed states. An allosteric binding site on the protein surface is recognized for the development of novel FABP4 inhibitors. The binding affinity is calculated and compared with the experimental value. The kinetic analysis suggests that ligand residence on the protein surface may delay the binding process. Overall, our results provide a comprehensive picture of ligand diffusion on the protein surface, ligand migration into the buried cavity, and the conformational change of FABP4 at an atomic level.