Fatty acid binding to the human transport proteins FABP3, FABP4, and FABP5 from a Ligand's perspective.

Fatty acid binding proteins (FABPs) are a family of amphiphilic transport proteins with high diversity in terms of their amino acid sequences and binding preferences. Beyond their main biological role as cytosolic fatty acid transporters, many aspects regarding their binding mechanism and functional specializations in human cells remain unclear. In this work, the binding properties and thermodynamics of FABP3, FABP4, and FABP5 were analyzed under various physical conditions. For this purpose, the FABPs were loaded with fatty acids bearing fluorescence or spin probes as model ligands, comparing their binding affinities via microscale thermophoresis (MST) and continuous-wave electron paramagnetic resonance (CW EPR) spectroscopy. The CW EPR spectra of non-covalently bound 5- and 16-DOXYL stearic acid (5/16-DSA) deliver in-depth information about the dynamics and chemical environments of ligands inside the binding pockets of the FABPs. EPR spectral simulations allow the construction of binding curves, revealing two different binding states ('intermediately' and 'strongly' bound). The proportion of bound 5/16-DSA depends strongly on the FABP concentration and the temperature but with remarkable differences between the three isoforms. Additionally, the more dynamic state ('intermediately bound') seems to dominate at body temperature with thermodynamic preference. The ligand binding studies were supplemented by aggregation studies via dynamic light scattering and bioinformatic analyses. Beyond the remarkably fine-tuned binding properties exhibited by each FABP, which were discernible with our EPR-centered approach, the results of this work attest to the power of simple spectroscopic experiments to provide new insights into the ligand binding mechanisms of proteins in general on a molecular level.

Integrating Explicit and Implicit Fullerene Models into UNRES Force Field for Protein Interaction Studies.

Fullerenes, particularly C60, exhibit unique properties that make them promising candidates for various applications, including drug delivery and nanomedicine. However, their interactions with biomolecules, especially proteins, remain not fully understood. This study implements both explicit and implicit C60 models into the UNRES coarse-grained force field, enabling the investigation of fullerene-protein interactions without the need for restraints to stabilize protein structures. The UNRES force field offers computational efficiency, allowing for longer timescale simulations while maintaining accuracy. Five model proteins were studied: FK506 binding protein, HIV-1 protease, intestinal fatty acid binding protein, PCB-binding protein, and hen egg-white lysozyme. Molecular dynamics simulations were performed with and without C60 to assess protein stability and investigate the impact of fullerene interactions. Analysis of contact probabilities reveals distinct interaction patterns for each protein. FK506 binding protein (1FKF) shows specific binding sites, while intestinal fatty acid binding protein (1ICN) and uteroglobin (1UTR) exhibit more generalized interactions. The explicit C60 model shows good agreement with all-atom simulations in predicting protein flexibility, the position of C60 in the binding pocket, and the estimation of effective binding energies. The integration of explicit and implicit C60 models into the UNRES force field, coupled with recent advances in coarse-grained modeling and multiscale approaches, provides a powerful framework for investigating protein-nanoparticle interactions at biologically relevant scales without the need to use restraints stabilizing the protein, thus allowing for large conformational changes to occur. These computational tools, in synergy with experimental techniques, can aid in understanding the mechanisms and consequences of nanoparticle-biomolecule interactions, guiding the design of nanomaterials for biomedical applications.

A combined computational-biophysical approach to understanding fatty acid binding to FABP7.

Members of the fatty acid binding protein (FABP) family function as intracellular transporters of long-chain fatty acids and other hydrophobic molecules to different cellular compartments. Brain FABP (FABP7) exhibits ligand-directed differences in cellular transport. For example, when FABP7 binds to docosahexaenoic acid (DHA), the complex relocates to the nucleus and influences transcriptional activity, whereas FABP7 bound with monosaturated fatty acids remains in the cytosol. Preferential binding of FABP7 to polyunsaturated fatty acids like DHA has been previously observed and is thought to play a role in differential localization. However, we find that at 37°C, FABP7 does not display strong selectivity, suggesting that the conformational ensemble of FABP7 and its perturbation upon binding may be important. We use molecular dynamics simulations, NMR, and a variety of biophysical techniques to better understand the conformational ensemble of FABP7, how it is perturbed by fatty acid binding, and how this may be related to ligand-directed transport. We find that FABP7 has high degree of conformational heterogeneity that is substantially reduced upon ligand binding. We also observe substantial heterogeneity in ligand binding poses, which is consistent with our finding that ligand binding is resistant to mutations in key polar residues in the binding pocket. Our NMR experiments show that DHA binding leads to chemical shift perturbations in residues near the nuclear localization signal, which may point toward a mechanism of differential transport.

Understanding FABP7 binding to fatty acid micelles and membranes.

Fatty acid-binding proteins (FABPs) are chaperones that facilitate the transport of long-chain fatty acids within the cell and can provide cargo-dependent localization to specific cellular compartments. Understanding the nature of this transport is important because lipid signaling functions are associated with metabolic pathways impacting disease pathologies including cancer, autism, and schizophrenia. FABPs often associate with cell membranes to acquire and deliver their bound cargo as part of transport. We focus on brain FABP (FABP7), which demonstrates localization to the cytoplasm and nucleus, influencing transcription and fatty acid metabolism. We use a combined biophysical-computational approach to elucidate the interaction between FABP7 and model membranes. Specifically, we use multiple experiments to demonstrate that FABP7 can bind oleic acid and docosahexaenoic acid micelles. Data from NMR and multiscale molecular dynamics simulations reveal that the interaction with micelles is through FABP7's portal region residues. Simulations suggest that binding to membranes occurs through the same residues as micelles. Simulations also capture binding events where fatty acids dissociate from the membrane and enter FABP7's binding pocket. Overall, our data shed light on the interactions between FABP7 and OA or DHA micelles and provide insight into the transport of long-chain fatty acids.

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.

On the Dynamic Mechanism of Long-Flexible Fatty Acid Binding to Fatty Acid Binding Protein: Resolving the Long-Standing Debate.

Fatty acids (FAs) are one of the essential energy sources for physiological processes, and they play a vital role in regulating immune and inflammatory responses, promoting cell differentiation and apoptosis, and inhibiting tumor growth. These functions are carried out by FA binding proteins (FABPs) that recognize and transport FAs. Although the crystal structure of the FA-FABPs complex has long been characterized, the mechanism behind FA binding and dissociation from FABP remains unclear. This study employed conventional MD simulations and enhanced sampling technologies to investigate the atomic-scale complexes of heart fatty acid binding proteins and stearic acid (SA). The results revealed two primary pathways for the binding or dissociation of the flexible long-chain ligand, with the orientation of the SA carboxyl head during dissociation determining the chosen path. Conformational changes in the portal region of FABP during the ligand binding/unbinding were found to be trivial, and the overturn of the ″cap″ or the unfolding of the α2 helix was not required. This study resolves the long-standing debate on the binding mechanism of SA with the long-flexible tail to FABP, which significantly improves the understanding of the transport mechanism of FABPs and the development of related therapeutic agents.

Disclosing Environmental Ligands of L-FABP and PPARγ: Should We Re-evaluate the Chemical Safety of Hydrocarbon Surfactants?

Chemical contaminants can cause adverse effects by binding to the liver-fatty acid binding protein (L-FABP) and peroxisome proliferator-activated nuclear receptor γ (PPARγ), which are vital in lipid metabolism. However, the presence of numerous compounds in the environment has hindered the identification of their ligands, and thus only a small portion have been discovered to date. In this study, protein Affinity Purification with Nontargeted Analysis (APNA) was employed to identify the ligands of L-FABP and PPARγ in indoor dust and sewage sludge. A total of 83 nonredundant features were pulled-out by His-tagged L-FABP as putative ligands, among which 13 were assigned as fatty acids and hydrocarbon surfactants. In contrast, only six features were isolated when His-tagged PPARγ LBD was used as the protein bait. The binding of hydrocarbon surfactants to L-FABP and PPARγ was confirmed using both recombinant proteins and reporter cells. These hydrocarbon surfactants, along with >50 homologues and isomers, were detected in dust and sludge at high concentrations. Fatty acids and hydrocarbon surfactants explained the majority of L-FABP (57.7 ± 32.9%) and PPARγ (66.0 ± 27.1%) activities in the sludge. This study revealed hydrocarbon surfactants as the predominant synthetic ligands of L-FABP and PPARγ, highlighting the importance of re-evaluating their chemical safety.

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.

Ligand Entry into Fatty Acid Binding Protein via Local Unfolding Instead of Gap Widening.

Fatty acid binding proteins play an important role in the transportation of fatty acids. Despite intensive studies, how fatty acids enter the protein cavity for binding is still controversial. Here, a gap-closed variant of human intestinal fatty acid binding protein was generated by mutagenesis, in which the gap is locked by a disulfide bridge. According to its structure determined here by NMR, this variant has no obvious openings as the ligand entrance and the gap cannot be widened by internal dynamics. Nevertheless, it still takes up fatty acids and other ligands. NMR relaxation dispersion, chemical exchange saturation transfer, and hydrogen-deuterium exchange experiments show that the variant exists in a major native state, two minor native-like states, and two locally unfolded states in aqueous solution. Local unfolding of either ßB-ßD or helix 2 can generate an opening large enough for ligands to enter the protein cavity, but only the fast local unfolding of helix 2 is relevant to the ligand entry process.

A ligand-induced structural change in fatty acid-binding protein 1 is associated with potentiation of peroxisome proliferator-activated receptor α agonists.

Peroxisome proliferator-activated receptor α (PPARα) is a transcriptional regulator of lipid metabolism. GW7647 is a potent PPARα agonist that must reach the nucleus to activate this receptor. In cells expressing human fatty acid-binding protein 1 (FABP1), GW7647 treatment increases FABP1's nuclear localization and potentiates GW7647-mediated PPARα activation; GW7647 is less effective in cells that do not express FABP1. To elucidate the underlying mechanism, here we substituted residues in FABP1 known to dictate lipid signaling by other intracellular lipid-binding proteins. Substitutions of Lys-20 and Lys-31 to Ala in the FABP1 helical cap affected neither its nuclear localization nor PPARα activation. In contrast, Ala substitution of Lys-57, Glu-77, and Lys-96, located in the loops adjacent to the ligand-binding portal region, abolished both FABP1 nuclear localization and GW7647-induced PPARα activation but had little effect on GW7647-FABP1 binding affinity. Using solution NMR spectroscopy, we determined the WT FABP1 structure and analyzed the dynamics in the apo and GW7647-bound structures of both the WT and the K57A/E77A/K96A triple mutant. We found that GW7647 binding causes little change in the FABP1 backbone, but solvent exposes several residues in the loops around the portal region, including Lys-57, Glu-77, and Lys-96. These residues also become more solvent-exposed upon binding of FABP1 with the endogenous PPARα agonist oleic acid. Together with previous observations, our findings suggest that GW7647 binding stabilizes a FABP1 conformation that promotes its interaction with PPARα. We conclude that full PPARα agonist activity of GW7647 requires FABP1-dependent transport and nuclear localization processes.

Acyl-chain selectivity and physiological roles of Staphylococcus aureus fatty acid-binding proteins.

Fatty acid (FA) kinase produces acyl-phosphate for the synthesis of membrane phospholipids in Gram-positive bacterial pathogens. FA kinase consists of a kinase protein (FakA) that phosphorylates an FA substrate bound to a second module, an FA-binding protein (FakB). Staphylococcus aureus expresses two distinct, but related, FakBs with different FA selectivities. Here, we report the structures of FakB1 bound to four saturated FAs at 1.6-1.93 Å resolution. We observed that the different FA structures are accommodated within a slightly curved hydrophobic cavity whose length is governed by the conformation of an isoleucine side chain at the end of the tunnel. The hydrophobic tunnel in FakB1 prevents the binding of cis-unsaturated FAs, which are instead accommodated by the kinked tunnel within the FakB2 protein. The differences in the FakB interiors are not propagated to the proteins' surfaces, preserving the protein-protein interactions with their three common partners, FakA, PlsX, and PlsY. Using cellular thermal shift analyses, we found that FakB1 binds FA in vivo, whereas a significant proportion of FakB2 does not. Incorporation of exogenous FA into phospholipid in ΔfakB1 and ΔfakB2 S. aureus knockout strains revealed that FakB1 does not efficiently activate unsaturated FAs. FakB2 preferred unsaturated FAs, but also allowed the incorporation of saturated FAs. These results are consistent with a model in which FakB1 primarily functions in the recycling of the saturated FAs produced by S. aureus metabolism, whereas FakB2 activates host-derived oleate, which S. aureus does not produce but is abundant at infection sites.

The crystal structure of ORP3 reveals the conservative PI4P binding pattern.

Oxysterol-binding protein (OSBP) and its related protein (ORP) constitute a conserved family of lipid transfer proteins (LTPs). ORPs have been implicated as intracellular lipid exchanger and sensor in recent years, which regulate the lipid homeostasis and signal pathway. OSBP-related protein 3 plays key role in controlling cell adhesion and migration and could be developed as the drug target for cancer therapy. Here, we report the crystal structures of human ORP3 ORD to 2.1 Å and ORD-PI4P complex to 3.2 Å. The binding assay in vitro confirms the ORP3 has the capability of PI4P binding. This study further verifies that the PI4P is the common ligand of all ORPs and ORPs should be the lipid exchanger in membrane contact sites(MCS).

Proteomic profiling of fatty acid binding proteins in muscular dystrophy.

Introduction: Duchenne muscular dystrophy is a neuromuscular disorder, which is caused by abnormalities in the DMD gene that encodes the membrane cytoskeletal protein dystrophin. Besides progressive skeletal muscle wasting, dystrophinopathy also affects non-skeletal muscle tissues, including cells in the cardio-respiratory system, the central nervous system, the liver and the kidney.Areas covered: This review summarizes the proteomic characterization of a key class of lipid chaperones, the large family of fatty acid binding proteins, and their potential role in muscular dystrophy. Recent proteomic surveys using animal models and patient specimens are reviewed. Pathobiochemical changes in specific proteoforms of fatty acid binding protein in the multi-system pathology of dystrophinopathy are discussed.Expert opinion: The mass spectrometric identification of distinct changes in fatty acid binding proteins in muscle, heart, liver, kidney and serum demonstrates that considerable alterations occur in key steps of metabolite transport and fat metabolism in muscular dystrophy. These new findings might be helpful to further develop a comprehensive biomarker signature of metabolic changes in X-linked muscular dystrophy, which should improve (i) our understanding of complex pathobiochemical changes due to dystrophin deficiency, (ii) the identification of novel therapeutic targets, and (iii) the design of differential diagnostic, prognostic and therapy-monitoring approaches.

Evolution of chalcone isomerase from a noncatalytic ancestor.

The emergence of catalysis in a noncatalytic protein scaffold is a rare, unexplored event. Chalcone isomerase (CHI), a key enzyme in plant flavonoid biosynthesis, is presumed to have evolved from a nonenzymatic ancestor related to the widely distributed fatty-acid binding proteins (FAPs) and a plant protein family with no isomerase activity (CHILs). Ancestral inference supported the evolution of CHI from a protein lacking isomerase activity. Further, we identified four alternative founder mutations, i.e., mutations that individually instated activity, including a mutation that is not phylogenetically traceable. Despite strong epistasis in other cases of protein evolution, CHI's laboratory reconstructed mutational trajectory shows weak epistasis. Thus, enantioselective CHI activity could readily emerge despite a catalytically inactive starting point. Accordingly, X-ray crystallography, NMR, and molecular dynamics simulations reveal reshaping of the active site toward a productive substrate-binding mode and repositioning of the catalytic arginine that was inherited from the ancestral fatty-acid binding proteins.

Potential Antidiabetic Fumiquinazoline Alkaloids from the Marine-Derived Fungus Scedosporium apiospermum F41-1.

Fumiquinazoline alkaloids have attracted much attention from medicinal and natural product chemists due to their interesting structures and biological potential. In this study, three new and 12 known fumiquinazoline alkaloids were isolated and characterized from the marine fungus Scedosporium apiospermum F41-1. The structures of the new compounds and their absolute configurations were determined using NMR spectroscopy, ECD, and OR calculations. The compounds were evaluated for their antidiabetic potential by determining their triglyceride-promoting activity using 3T3-L1 adipocytes. One of the new compounds, scequinadoline J (14), as well as scequinadolines D (9) and E (10), was found to promote triglyceride accumulation in 3T3-L1 cells. Scequinadoline D (9) demonstrated the most potent activity, with an EC50 value of 0.27 ± 0.03 µM. Quantitative polymerase chain reaction experiments suggested that scequinadoline D (9) acts through activation of the PPARγ pathway. It stimulated the mRNA expression of PPARγ, AMPKα, C/EBPα, LXRα, SCD-1, and FABP4. In addition, its triglyceride-promoting efficacy could be blocked by a double dose of the PPARγ antagonist GW9662. These results indicated that scequinadoline D (9) is a potent insulin sensitizer that targets adipocytes and may be useful for the treatment of type 2 diabetes mellitus after further investigation.

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.

Identification of Fatty Acid Binding Protein 5 Inhibitors Through Similarity-Based Screening.

Fatty acid binding protein 5 (FABP5) is a promising target for development of inhibitors to help control pain and inflammation. In this work, computer-based docking (DOCK6 program) was employed to screen ∼2 M commercially available compounds to FABP5 based on an X-ray structure complexed with the small molecule inhibitor SBFI-26 previously identified by our group (also through virtual screening). The goal was discovery of additional chemotypes. The screen resulted in the purchase of 78 candidates, which led to the identification of a new inhibitor scaffold (STK-0) with micromolar affinity and apparent selectivity for FABP5 over FABP3. A second similarity-based screen resulted in three additional hits (STK-15, STK-21, STK-22) from which preliminary SAR could be derived. Notably, STK-15 showed comparable activity to the SBFI-26 reference under the same assay conditions (1.40 vs 0.86 µM). Additional molecular dynamics simulations, free energy calculations, and structural analysis (starting from DOCK-generated poses) revealed that R enantiomers (dihydropyrrole scaffold) of STK-15 and STK-22 have a more optimal composition of functional groups to facilitate additional H-bonds with Arg109 of FABP5. This observation suggests enantiomerically pure compounds could show enhanced activity. Overall, our study highlights the utility of using similarity-based screening methods to discover new inhibitor chemotypes, and the identified FABP5 hits provide a strong starting point for future efforts geared to improve activity.

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.

Cannabinoid Interactions with Proteins: Insights from Structural Studies.

Cannabinoids have been widely used for recreational and medicinal purposes. The increasing legalization of cannabinoid use and the growing success in Medicinal Chemistry of cannabinoids have fueled recent interest in cannabinoid-sensing sites in receptor proteins. Here, we review structural data from high-resolution cryo-EM and crystallography studies that depict phytocannabinoid, endocannabinoid, and synthetic cannabinoid molecules bound to various proteins. The latter include antigen-binding fragment (Fab), cellular retinol binding protein 2 (CRBP2), fatty acid-binding protein 5 (FABP5), peroxisome proliferator-activated receptor γ (PPAR γ), and cannabinoid receptor types 1 and 2 (CB1 and CB2). Cannabinoid-protein complexes reveal the complex design of cannabinoid binding sites that are usually presented by conventional ligand-binding pockets on respective proteins. However, subtle differences in cannabinoid interaction with amino acids within the binding pocket often result in diverse consequences for protein function. The rapid increase in available structural data on cannabinoid-protein interactions will ultimately direct drug design efforts toward rendering highly potent cannabinoid-related pharmacotherapies that are devoid of side effects.

The Observation of Ligand-Binding-Relevant Open States of Fatty Acid Binding Protein by Molecular Dynamics Simulations and a Markov State Model.

As a member of the fatty acids transporter family, the heart fatty acid binding proteins (HFABPs) are responsible for many important biological activities. The binding mechanism of fatty acid with FABP is critical to the understanding of FABP functions. The uncovering of binding-relevant intermediate states and interactions would greatly increase our knowledge of the binding process. In this work, all-atom molecular dynamics (MD) simulations were performed to characterize the structural properties of nativelike intermediate states. Based on multiple 6 µs MD simulations and Markov state model (MSM) analysis, several "open" intermediate states were observed. The transition rates between these states and the native closed state are in good agreement with the experimental measurements, which indicates that these intermediate states are binding relevant. As a common property in the open states, the partially unfolded α2 helix generates a larger portal and provides the driving force to facilitate ligand binding. On the other side, there are two kinds of open states for the ligand-binding HFABP: one has the partially unfolded α2 helix, and the other has the looser ß-barrel with disjointing ßD-ßE strands. Our results provide atomic-level descriptions of the binding-relevant intermediate states and could improve our understanding of the binding mechanism.