Millisecond Timescale Dynamics of Human Liver Fatty Acid Binding Protein: Testing of Its Relevance to the Ligand Entry Process

For over a decade, scientists have been attempting to know more about the conformational dynamics of fatty acid binding proteins (FABPs), to answer the puzzling question of how ligands could access the internalized binding site(s). Conformational exchange of FABPs on the microsecond to millisecond timescales has been found in many FABPs and offers an important hypothesis for the ligand entry mechanism. Despite the potential significance, the validity of this hypothesis has not been verified yet. In this study, the slow dynamics of human liver fatty acid binding protein (hLFABP) that was shown previously to be highly flexible on millisecond timescales was quantitatively characterized in detail. In addition, the interaction between hLFABP and 1,8-ANS was studied using NMR spectroscopy, and the kinetic rate of ANS association to hLFABP was measured. We believe the current result excludes the possibility that the intrinsic millisecond dynamics of hLFABP represents a critical conformational reorganization process required for ligand entry, but implies that it may represent the exchange between the apo-state and a state resembling the singly-bound conformation. Furthermore, we suggest these results show that the ligand-entry related functional dynamics could occur on the microsecond/submicrosecond timescales, highly encouraging future computational studies on this topic.     

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.     

Molecular Dynamics Simulation of Ligand Dissociation from Liver Fatty Acid Binding Protein

The mechanisms of how ligands enter and leave the binding cavity of fatty acid binding proteins (FABPs) have been a puzzling question over decades. Liver fatty acid binding protein (LFABP) is a unique family member which accommodates two molecules of fatty acids in its cavity and exhibits the capability of interacting with a variety of ligands with different chemical structures and properties. Investigating the ligand dissociation processes of LFABP is thus a quite interesting topic, which however is rather difficult for both experimental approaches and ordinary simulation strategies. In the current study, random expulsion molecular dynamics simulation, which accelerates ligand motions for rapid dissociation, was used to explore the potential egress routes of ligands from LFABP. The results showed that the previously hypothesized “portal region” could be readily used for the dissociation of ligands at both the low affinity site and the high affinity site. Besides, one alternative portal was shown to be highly favorable for ligand egress from the high affinity site and be related to the unique structural feature of LFABP. This result lends strong support to the hypothesis from the previous NMR exchange studies, which in turn indicates an important role for this alternative portal. Another less favored potential portal located near the N-terminal end was also identified. Identification of the dissociation pathways will allow further mechanistic understanding of fatty acid uptake and release by computational and/or experimental techniques.     

A Thermoacidophile-Specific Protein Family,DUF3211, Functions as a Fatty Acid Carrier with Novel Binding Mode

STK_08120 is a member of the thermoacidophile-specific DUF3211 protein family from Sulfolobus tokodaii strain 7. Its molecular function remains obscure, and sequence similarities for obtaining functional remarks are not available. In this study, the crystal structure of STK_08120 was determined at 1.79-Å resolution to predict its probable function using structure similarity searches. The structure adopts an α/β structure of a helix-grip fold, which is found in the START domain proteins with cavities for hydrophobic substrates or ligands. The detailed structural features implied that fatty acids are the primary ligand candidates for STK_08120, and binding assays revealed that the protein bound long-chain saturated fatty acids (>C₁₄) and their trans-unsaturated types with an affinity equal to that for major fatty acid binding proteins in mammals and plants. Moreover, the structure of an STK_08120-myristic acid complex revealed a unique binding mode among fatty acid binding proteins. These results suggest that the thermoacidophile-specific protein family DUF3211 functions as a fatty acid carrier with a novel binding mode.     

Expression of rat L-FABP in mouse fibroblasts: role in fat absorption

Fatty acid-binding proteins (FABP) are abundant cytosolic proteins whose level is responsive to nutritional, endocrine, and a variety of pathological states. Although FABPs have been investigatedin vitro for several decades, little is known of their physiological function. Liver L-FABP binds both fatty acids and cholesterol. Competitive binding analysis and molecular modeling studies of L-FABP indicate the presence of two ligand binding pockets that accomodate one fatty acid each. One fatty acid binding site is identical to the cholesterol binding site. To test whether these observations obtainedin vitro were physiologically relevant, the cDNA encoding L-FABP was transfected into L-cells, a cell line with very low endogenous FABP and sterol carrier proteins. Uptake of both ligands did not differ between control cells and low expression clones. In contrast, both fatty acid uptake and cholesterol uptake were stimulated in the high expression cells. In high expression cells, uptake of fluorescent cis-parinaric acid was enhanced more than that of trans-parinaric acid. This is consistent with the preferential binding of cis-fatty acids to L-FABP but in contrast to the preferential binding of trans-parinaric acid to the L-cell plasma membrane fatty acid transporter (PMFABP). These data show that the level of cytosolic fatty acids in intact cells can regulate both the extent and specificity of fatty acid uptake. Last, sphingomyelinase treatment of L-cells released cholesterol from the plasma membrane to the cytoplasm and stimulated microsomal acyl-CoA: cholesteryl acyl transferase (ACAT). This process was accelerated in high expression cells. These observations show for the first time in intact cells that L-FABP, a protein most prevalent in liver and intestine where much fat absorption takes place, may have a role in fatty acid and cholesterol absorption.Abbreviations FABP fatty acid-binding protein - L-FABP liver fatty acid-binding protein - I-FABP intestinal fatty acid-binding protein - H-FABP heart fatty acid-binding protein - A-FABP adipocyte fatty acid-binding protein - PMFABP plasma membrane fatty acid-binding protein - SCP-2 sterol carrier protein-2 - Dehydroergosterol (DHE) d-5,7,9(11),22-ergostatetraene-3b-ol - cis-parinaric acid-9Z, 11E, 13E, 15Z-octatetraenoic acid - trans parinaric acid, 9E, 11E, 13E, 14E-octatetraenoic acid - BSA bovine serum albumin - KRH Krebs-Ringer-Henseleit buffer     

Decoding Structural Properties of a Partially Unfolded Protein Substrate: En Route to Chaperone Binding

Many proteins comprising of complex topologies require molecular chaperones to achieve their unique three-dimensional folded structure. The E.coli chaperone, GroEL binds with a large number of unfolded and partially folded proteins, to facilitate proper folding and prevent misfolding and aggregation. Although the major structural components of GroEL are well defined, scaffolds of the non-native substrates that determine chaperone-mediated folding have been difficult to recognize. Here we performed all-atomistic and replica-exchange molecular dynamics simulations to dissect non-native ensemble of an obligate GroEL folder, DapA. Thermodynamics analyses of unfolding simulations revealed populated intermediates with distinct structural characteristics. We found that surface exposed hydrophobic patches are significantly increased, primarily contributed from native and non-native β-sheet elements. We validate the structural properties of these conformers using experimental data, including circular dichroism (CD), 1-anilinonaphthalene-8-sulfonic acid (ANS) binding measurements and previously reported hydrogen-deutrium exchange coupled to mass spectrometry (HDX-MS). Further, we constructed network graphs to elucidate long-range intra-protein connectivity of native and intermediate topologies, demonstrating regions that serve as central “hubs”. Overall, our results implicate that genomic variations (or mutations) in the distinct regions of protein structures might disrupt these topological signatures disabling chaperone-mediated folding, leading to formation of aggregates.     

脂肪细胞型脂肪酸结合蛋白的研究进展

脂肪细胞型脂肪酸结合蛋白(adipocyte fatty acid binding protein,AFABP/aP2)作为脂肪酸结合蛋白(FABPS)超家族成员之一,广泛存在于各种正常的组织细胞中,参与脂肪酸贮存,运输与降解等过程。近年来,对脂肪细胞型脂肪酸结合蛋白的研究已成为热点,本文就其主要特征及其与各类疾病的关系作一简要综述。     

Analyzing Protein Dynamics Using Hydrogen Exchange Mass Spectrometry

All cellular processes depend on the functionality of proteins. Although the functionality of a given protein is the direct consequence of its unique amino acid sequence, it is only realized by the folding of the polypeptide chain into a single defined three-dimensional arrangement or more commonly into an ensemble of interconverting conformations. Investigating the connection between protein conformation and its function is therefore essential for a complete understanding of how proteins are able to fulfill their great variety of tasks. One possibility to study conformational changes a protein undergoes while progressing through its functional cycle is hydrogen-¹H/²H-exchange in combination with high-resolution mass spectrometry (HX-MS). HX-MS is a versatile and robust method that adds a new dimension to structural information obtained by e.g. crystallography. It is used to study protein folding and unfolding, binding of small molecule ligands, protein-protein interactions, conformational changes linked to enzyme catalysis, and allostery. In addition, HX-MS is often used when the amount of protein is very limited or crystallization of the protein is not feasible. Here we provide a general protocol for studying protein dynamics with HX-MS and describe as an example how to reveal the interaction interface of two proteins in a complex.        

脂肪细胞型脂肪酸结合蛋白的研究进展

脂肪细胞型脂肪酸结合蛋白（adipocyte fatty acid binding protein,AFABP/aP2）作为脂肪酸结合蛋白（FABPS）超家族成员之一,广泛存在于各种正常的组织细胞中,参与脂肪酸贮存,运输与降解等过程。近年来,对脂肪细胞型脂肪酸结合蛋白的研究已成为热点,本文就其主要特征及其与各类疾病的关系作一简要综述。     

SPRY Domain-Containing SOCS Box Protein 2: Crystal Structure and Residues Critical for Protein Binding

The four mammalian SPRY (a sequence repeat in dual-specificity kinase splA and ryanodine receptors) domain-containing suppressor of cytokine signalling (SOCS) box proteins (SSB-1 to -4) are characterised by a C-terminal SOCS box and a central SPRY domain. The latter is a protein interaction module found in over 1600 proteins, with more than 70 encoded in the human genome. Here we report the crystal structure of the SPRY domain of murine SSB-2 and compare it with the SSB-2 solution structure and crystal structures of other B30.2/SPRY domain-containing family proteins. The structure is a bent β-sandwich, consisting of two seven-stranded β-sheets wrapped around a long loop that extends from the centre strands of the inner or concave β-sheet; it closely matches those of GUSTAVUS and SSB-4. The structure is also similar to those of two recently determined Neuralized homology repeat (NHR) domains (also known as NEUZ domains), with detailed comparisons, suggesting that the NEUZ/NHR domains form a subclass of SPRY domains. The binding site on SSB-2 for the prostate apoptosis response-4 (Par-4) protein has been mapped in finer detail using mutational analyses. Moreover, SSB-1 was shown to have a Par-4 binding surface similar to that identified for SSB-2. Structural perturbations of SSB-2 induced by mutations affecting its interaction with Par-4 and/or c-Met have been characterised by NMR. These comparisons, in conjunction with previously published dynamics data from NMR relaxation studies and coarse-grained dynamics simulation using normal mode analysis, further refine our understanding of the structural basis for protein recognition of SPRY domain-containing proteins.     

From Nonspecific DNA–Protein Encounter Complexes to the Prediction of DNA–Protein Interactions

DNA–protein interactions are involved in many essential biological activities. Because there is no simple mapping code between DNA base pairs and protein amino acids, the prediction of DNA–protein interactions is a challenging problem. Here, we present a novel computational approach for predicting DNA-binding protein residues and DNA–protein interaction modes without knowing its specific DNA target sequence. Given the structure of a DNA-binding protein, the method first generates an ensemble of complex structures obtained by rigid-body docking with a nonspecific canonical B-DNA. Representative models are subsequently selected through clustering and ranking by their DNA–protein interfacial energy. Analysis of these encounter complex models suggests that the recognition sites for specific DNA binding are usually favorable interaction sites for the nonspecific DNA probe and that nonspecific DNA–protein interaction modes exhibit some similarity to specific DNA–protein binding modes. Although the method requires as input the knowledge that the protein binds DNA, in benchmark tests, it achieves better performance in identifying DNA-binding sites than three previously established methods, which are based on sophisticated machine-learning techniques. We further apply our method to protein structures predicted through modeling and demonstrate that our method performs satisfactorily on protein models whose root-mean-square Cα deviation from native is up to 5 Å from their native structures. This study provides valuable structural insights into how a specific DNA-binding protein interacts with a nonspecific DNA sequence. The similarity between the specific DNA–protein interaction mode and nonspecific interaction modes may reflect an important sampling step in search of its specific DNA targets by a DNA-binding protein.     

Principles and Methods in Computational Membrane Protein Design

After decades of progress in computational protein design, the design of proteins folding and functioning in lipid membranes appears today as the next frontier. Some notable successes in the de novo design of simplified model membrane protein systems have helped articulate fundamental principles of protein folding, architecture and interaction in the hydrophobic lipid environment. These principles are reviewed here, together with the computational methods and approaches that were used to identify them. We provide an overview of the methodological innovations in the generation of new protein structures and functions and in the development of membrane-specific energy functions. We highlight the opportunities offered by new machine learning approaches applied to protein design, and by new experimental characterization techniques applied to membrane proteins. Although membrane protein design is in its infancy, it appears more reachable than previously thought.     

Energy Landscapes of Protein Aggregation and Conformation Switching in Intrinsically Disordered Proteins

The protein folding problem was apparently solved recently by the advent of a deep learning method for protein structure prediction called AlphaFold. However, this program is not able to make predictions about the protein folding pathways. Moreover, it only treats about half of the human proteome, as the remaining proteins are intrinsically disordered or contain disordered regions. By definition these proteins differ from natively folded proteins and do not adopt a properly folded structure in solution. However these intrinsically disordered proteins (IDPs) also systematically differ in amino acid composition and uniquely often become folded upon binding to an interaction partner. These factors preclude solving IDP structures by current machine-learning methods like AlphaFold, which also cannot solve the protein aggregation problem, since this meta-folding process can give rise to different aggregate sizes and structures. An alternative computational method is provided by molecular dynamics simulations that already successfully explored the energy landscapes of IDP conformational switching and protein aggregation in multiple cases. These energy landscapes are very different from those of ‘simple’ protein folding, where one energy funnel leads to a unique protein structure. Instead, the energy landscapes of IDP conformational switching and protein aggregation feature a number of minima for different competing low-energy structures. In this review, I discuss the characteristics of these multifunneled energy landscapes in detail, illustrated by molecular dynamics simulations that elucidated the underlying conformational transitions and aggregation processes.     

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

Fatty acid binding proteins (FABPs), are evolutionarily conserved small cytoplasmic proteins that occur in many tissue-specific types. One of their primary functions is to facilitate the clearance of the cytoplasmic matrix from free fatty acids and of other detergent-like compounds. Crystallographic studies of FABP proteins have revealed a well defined binding site located deep inside their β-clam structure that is hardly exposed to the bulk solution. However, NMR measurements revealed that, when the protein is equilibrated with its ligands, residues that are clearly located on the outer surface of the protein do interact with the ligand. To clarify this apparent contradiction we applied molecular dynamics simulations to follow the initial steps associated with the FABP–fatty acid interaction using, as a model, the interaction of toad liver basic FABP, or chicken liver bile acid binding protein, with a physiological concentration of palmitate ions. The simulations (~200 ns of accumulated time) show that fatty acid molecules interact, unevenly, with various loci on the protein surface, with the favored regions being the portal and the anti-portal domains. Random encounters with palmitate at these regions led to lasting adsorption to the surface, while encounters at the outer surface of the β-clam were transient. Therefore, we suggest that the protein surface is capable of sequestering free fatty acids from solution, where brief encounters evolve into adsorbed states, which later mature by migration of the ligand into a more specific binding site.     

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.     

Dynamics of the BH3-Only Protein Binding Interface of Bcl-xL

The balance and interplay between pro-death and pro-survival members of the B-cell lymphoma-2 (Bcl-2) family proteins play key roles in regulation of the mitochondrial pathway of programmed cell death. Recent NMR and biochemical studies have revealed that binding of the proapoptotic BH3-only protein PUMA induces significant unfolding of antiapoptotic Bcl-xL at the interface, which in turn disrupts the Bcl-xL/p53 interaction to activate apoptosis. However, the molecular mechanism of such regulated unfolding of Bcl-xL is not fully understood. Analysis of the existing Protein Data Bank structures of Bcl-xL in both bound and unbound states reveal substantial intrinsic heterogeneity at its BH3-only protein binding interface. Large-scale atomistic simulations were performed in explicit solvent for six representative structures to further investigate the intrinsic conformational dynamics of Bcl-xL. The results support that the BH3-only protein binding interface of Bcl-xL is much more dynamic compared to the rest of the protein, both unbound and when bound to various BH3-only proteins. Such intrinsic interfacial conformational dynamics likely provides a physical basis that allows Bcl-xL to respond sensitively to detailed biophysical properties of the ligand. The ability of Bcl-xL to retain or even enhance dynamics at the interface in bound states could further facilitate the regulation of its interactions with various BH3-only proteins such as through posttranslational modifications.     

Structural determinants of cholesterol recognition in helical integral membrane proteins

Cholesterol is an integral component of mammalian membranes. It has been shown to modulate membrane fluidity and dynamics and alter integral membrane protein function. However, understanding the molecular mechanisms of how cholesterol impacts protein function is complicated by limited and conflicting structural data. Because of the nature of the crystallization and cryo-EM structure determination, it is difficult to distinguish between specific and biologically relevant interactions and a nonspecific association. The only widely recognized search algorithm for cholesterol-integral-membrane-protein interaction sites is sequence based, i.e., searching for the so-called "Cholesterol Recognition/interaction Amino acid Consensus" motif. Although these motifs are present in numerous integral membrane proteins, there is inconclusive evidence to support their necessity or sufficiency for cholesterol binding. Here, we leverage the increasing number of experimental cholesterol-integral-membrane-protein structures to systematically analyze putative interaction sites based on their spatial arrangement and evolutionary conservation. This analysis creates three-dimensional representations of general cholesterol interaction sites that form clusters across multiple integral membrane protein classes. We also classify cholesterol-integral-membrane-protein interaction sites as either likely-specific or nonspecific. Information gleaned from our characterization will eventually enable a structure-based approach to predict and design cholesterol-integral-membrane-protein interaction sites.     

Yolk proteins from nematodes, chickens, and frogs bind strongly and preferentially to left-handed Z-DNA

Yolk proteins purified from the nematode Caenorhabditis elegans, from the frog Xenopus laevis, and from chicken eggs all have the unexpected property of binding strongly and preferentially to a left-handed Z-DNA probe, brominated poly(dG-dC). We estimate that the nematode proteins bind to Z-DNA with an association constant of at least 10(4) (M-1) and that this association constant is at least 40-50-fold higher than the association constant to B-DNA. Thus, yolk proteins have a higher Z-DNA specificity than most of the Z-DNA binding proteins previously isolated from other sources. Although yolk protein binding to Z-DNA is poorly competed by a wide variety of nucleic acids, the interaction is strongly competed by the phospholipids cardiolipin and phosphatidic acid (500-1000-fold better than by the same mass of B-DNA). We suggest that Z-DNA interacts with the yolk protein phospholipid binding site. In general, our results emphasize the danger of using physical properties to infer biological function. In particular, our results should raise serious questions about the biological relevance of previously isolated Z-DNA binding proteins.     

Computer-based screening of functional conformers of proteins

A long-standing goal in biology is to establish the link between function, structure, and dynamics of proteins. Considering that protein function at the molecular level is understood by the ability of proteins to bind to other molecules, the limited structural data of proteins in association with other bio-molecules represents a major hurdle to understanding protein function at the structural level. Recent reports show that protein function can be linked to protein structure and dynamics through network centrality analysis, suggesting that the structures of proteins bound to natural ligands may be inferred computationally. In the present work, a new method is described to discriminate protein conformations relevant to the specific recognition of a ligand. The method relies on a scoring system that matches critical residues with central residues in different structures of a given protein. Central residues are the most traversed residues with the same frequency in networks derived from protein structures. We tested our method in a set of 24 different proteins and more than 260,000 structures of these in the absence of a ligand or bound to it. To illustrate the usefulness of our method in the study of the structure/dynamics/function relationship of proteins, we analyzed mutants of the yeast TATA-binding protein with impaired DNA binding. Our results indicate that critical residues for an interaction are preferentially found as central residues of protein structures in complex with a ligand. Thus, our scoring system effectively distinguishes protein conformations relevant to the function of interest.