X-ray crystallographic structures of adipocyte lipid-binding protein complexed with palmitate and hexadecanesulfonic acid. Properties of cavity binding sites

Heart-Type Fatty Acid Binding Protein Binds Long-Chain Acylcarnitines and Protects against Lipotoxicity

Heart-type fatty-acid binding protein (FABP3) is an essential cytosolic lipid transport protein found in cardiomyocytes. FABP3 binds fatty acids (FAs) reversibly and with high affinity. Acylcarnitines (ACs) are an esterified form of FAs that play an important role in cellular energy metabolism. However, an increased concentration of ACs can exert detrimental effects on cardiac mitochondria and lead to severe cardiac damage. In the present study, we evaluated the ability of FABP3 to bind long-chain ACs (LCACs) and protect cells from their harmful effects. We characterized the novel binding mechanism between FABP3 and LCACs by a cytotoxicity assay, nuclear magnetic resonance, and isothermal titration calorimetry. Our data demonstrate that FABP3 is capable of binding both FAs and LCACs as well as decreasing the cytotoxicity of LCACs. Our findings reveal that LCACs and FAs compete for the binding site of FABP3. Thus, the protective mechanism of FABP3 is found to be concentration dependent.

Unveiling the Role of the Fatty Acid Binding Protein 4 in the Metabolic-Associated Fatty Liver Disease

Metabolic-associated fatty liver disease (MAFLD), the main cause of chronic liver disease worldwide, is a progressive disease ranging from fatty liver to steatohepatitis (metabolic-associated steatohepatitis; MASH). Nevertheless, it remains underdiagnosed due to the lack of effective non-invasive methods for its diagnosis and staging. Although MAFLD has been found in lean individuals, it is closely associated with obesity-related conditions. Adipose tissue is the main source of liver triglycerides and adipocytes act as endocrine organs releasing a large number of adipokines and pro-inflammatory mediators involved in MAFLD progression into bloodstream. Among the adipocyte-derived molecules, fatty acid binding protein 4 (FABP4) has been recently associated with fatty liver and additional features of advanced stages of MAFLD. Additionally, emerging data from preclinical studies propose FABP4 as a causal actor involved in the disease progression, rather than a mere biomarker for the disease. Therefore, the FABP4 regulation could be considered as a potential therapeutic strategy to MAFLD. Here, we review the current knowledge of FABP4 in MAFLD, as well as its potential role as a therapeutic target for this disease.

Concise Modular Synthesis and NMR Structural Determination of Gallium Mycobactin T

A modular synthesis of mycobactin T and its N-acetyl analogue is reported in a route that facilitates permutation of the lipid tails. A key feature is the generation of N(α)-Cbz-N(ε)-benzyloxy-N(ε)-Boc-lysine (A4) with methyl(trifluoromethyl)dioxirane in 59% yield. Selective hydroxamate N-acylation was achieved with acyl fluorides, enabling installation of lipids tails in the final step. O-Benzyl-dehydrocobactin T (B4) was prepared by modifying a known five-step sequence with an overall yield of 49%. 2-Hydroxyphenyl-4-carboxyloxazoline (C3) was prepared from 2-hydroxybenzoic acid and l-serine methyl ester in three steps with an overall yield of 55%. Ester coupling of A4 and B4 with EDCI afforded MbI-1 in 73% yield. Catalytic hydrogenation with Pd/BaSO₄ and 50 psi of H₂ simultaneously effected alkene reduction and debenzylation to afford MbI-2 in 96% yield. Fragment C3 was converted into acyl fluoride C4, which coupled with MbI-2 to afford MbI-3 in 51% yield. Finally, Boc-removal with HCl/EtOAc and treatment of the resultant hydroxylamine with stearyl fluoride furnished mycobactin T in 65% yield. Overall, the yield is 4% over 14 steps. The gallium mycobactin T-N-acetyl derivative (GaMbT-NAc) structure was determined by ¹H NMR. The structure shows an octahedral Ga and two internal hydrogen bonds between peptidic N-Hs and two of the oxygen atoms coordinating Ga.

Molecular mechanisms of inhibitor bindings to A-FABP deciphered by using molecular dynamics simulations and calculations of MM-GBSA

Adipocyte fatty-acid binding protein (A-FABP) plays a central role in many aspects of metabolic diseases. It is an important target in drug design for treatment of FABP-related diseases. In this study, molecular dynamics (MD) simulations followed by calculations of molecular mechanics generalized Born surface area (MM-GBSA) and principal components analysis (PCA) were implemented to decipher molecular mechanism correlating with binding of inhibitors 57Q, 57P and L96 to A-FABP. The results show that van der Waals interactions are the leading factors to control associations of 57Q, 57P, and L96 with A-FABP, which reveals an energetic basis for designing of clinically available inhibitors towards A-FABP. The information from PCA and cross-correlation analysis rationally unveils that inhibitor bindings affect conformational changes of A-FABP and change relative movements between residues. Decomposition of binding affinity into contributions of individual residues not only detects hot spots of inhibitor/A-FABP binding but also shows that polar interactions of the positively charged residue Arg126 with three inhibitors provide a significant contribution for stabilization of the inhibitor/A-FABP bindings. Furthermore, the binding strength of L96 to residues Ser55, Phe57 and Lys58 are stronger than that of inhibitors 57Q and 57P to these residues.

Identification of Cryptic Binding Sites Using MixMD with Standard and Accelerated Molecular Dynamics

Protein dynamics play an important role in small molecule binding and can pose a significant challenge in the identification of potential binding sites. Cryptic binding sites have been defined as sites which require significant rearrangement of the protein structure to become physically accessible to a ligand. Mixed-solvent MD (MixMD) is a computational protocol which maps the surface of the protein using molecular dynamics (MD) of the unbound protein solvated in a 5% box of probe molecules with explicit water. This method has successfully identified known active and allosteric sites which did not require reorganization. In this study, we apply the MixMD protocol to identify known cryptic sites of 12 proteins characterized by a wide range of conformational changes. Of these 12 proteins, three require reorganization of side chains, five require loop movements, and four require movement of more significant structures such as whole helices. In five cases, we find that standard MixMD simulations are able to map the cryptic binding sites with at least one probe type. In two cases (guanylate kinase and TIE-2), accelerated MD, which increases sampling of torsional angles, was necessary to achieve mapping of portions of the cryptic binding site missed by standard MixMD. For more complex systems where movement of a helix or domain is necessary, MixMD was unable to map the binding site even with accelerated dynamics, possibly due to the limited timescale (100 ns for individual simulations). In general, similar conformational dynamics are observed in water-only simulations and those with probe molecules. This could imply that the probes are not driving opening events but rather take advantage of mapping sites that spontaneously open as part of their inherent conformational behavior. Finally, we show that docking to an ensemble of conformations from the standard MixMD simulations performs better than docking the apo crystal structure in nine cases and even better than half of the bound crystal structures. Poorer performance was seen in docking to ensembles of conformations from the accelerated MixMD simulations.

Role of the fatty acid-binding protein 4 in heart failure and cardiovascular disease

Obesity and ectopic fat accumulation in non-adipose tissues are major contributors to heart failure (HF) and cardiovascular disease (CVD). Adipocytes act as endocrine organs by releasing a large number of bioactive molecules into the bloodstream, which participate in a communication network between white adipose tissue and other organs, including the heart. Among these molecules, fatty acid-binding protein 4 (FABP4) has recently been shown to increase cardiometabolic risk. Both clinical and experimental evidence have identified FABP4 as a relevant player in atherosclerosis and coronary artery disease, and it has been directly related to cardiac alterations such as left ventricular hypertrophy (LVH) and both systolic and diastolic cardiac dysfunction. The available interventional studies preclude the establishment of a direct causal role of this molecule in CVD and HF and propose FABP4 as a biomarker rather than as an aetiological factor. However, several experimental reports have suggested that FABP4 may act as a direct contributor to cardiac metabolism and physiopathology, and the pharmacological targeting of FABP4 may restore some of the metabolic alterations that are conducive to CVD and HF. Here, we review the current knowledge regarding FABP4 in the context of HF and CVD as well as the molecular basis by which this protein participates in the regulation of cardiac function.

Interaction Analysis of FABP4 Inhibitors by X-ray Crystallography and Fragment Molecular Orbital Analysis

X-ray crystal structural determination of FABP4 in complex with four inhibitors revealed the complex binding modes, and the resulting observations led to improvement of the inhibitory potency of FABP4 inhibitors. However, the detailed structure-activity relationship (SAR) could not be explained from these structural observations. For a more detailed understanding of the interactions between FABP4 and inhibitors, fragment molecular orbital analyses were performed. These analyses revealed that the total interfragment interaction energies of FABP4 and each inhibitor correlated with the ranking of the K i value for the four inhibitors. Furthermore, interactions between each inhibitor and amino acid residues in FABP4 were identified. The oxygen atom of Lys58 in FABP4 was found to be very important for strong interactions with FABP4. These results might provide useful information for the development of novel potent FABP4 inhibitors.

High-resolution neutron and X-ray diffraction room-temperature studies of an H-FABP-oleic acid complex: study of the internal water cluster and ligand binding by a transferred multipolar electron-density distribution

Crystal diffraction data of heart fatty acid binding protein (H-FABP) in complex with oleic acid were measured at room temperature with high-resolution X-ray and neutron protein crystallography (0.98 and 1.90 Å resolution, respectively). These data provided very detailed information about the cluster of water molecules and the bound oleic acid in the H-FABP large internal cavity. The jointly refined X-ray/neutron structure of H-FABP was complemented by a transferred multipolar electron-density distribution using the parameters of the ELMAMII library. The resulting electron density allowed a precise determination of the electrostatic potential in the fatty acid (FA) binding pocket. Bader's quantum theory of atoms in molecules was then used to study interactions involving the internal water molecules, the FA and the protein. This approach showed H⋯H contacts of the FA with highly conserved hydrophobic residues known to play a role in the stabilization of long-chain FAs in the binding cavity. The determination of water hydrogen (deuterium) positions allowed the analysis of the orientation and electrostatic properties of the water molecules in the very ordered cluster. As a result, a significant alignment of the permanent dipoles of the water molecules with the protein electrostatic field was observed. This can be related to the dielectric properties of hydration layers around proteins, where the shielding of electrostatic interactions depends directly on the rotational degrees of freedom of the water molecules in the interface.

Manganese-mefenamic acid complexes exhibit high lipoxygenase inhibitory activity

The coordination of non-steroidal anti-inflammatory drugs (NSAIDs) to metal ions could improve the pharmaceutical efficacy of NSAIDs due to the unique characteristics of metal complexes.

Structure of the human-heart fatty-acid-binding protein 3 in complex with the fluorescent probe 1-anilinonaphthalene-8-sulphonic acid

Heart-type fatty-acid-binding protein (FABP3), which is a cytosolic protein abundantly found in cardiomyocytes, plays a role in trafficking fatty acids throughout cellular compartments by reversibly binding intracellular fatty acids with relatively high affinity. The fluorescent probe 1-anilinonaphthalene-8-sulfonate (ANS) is extensively utilized for examining the interaction of ligands with fatty-acid-binding proteins. The X-ray structure of FABP3 was determined in the presence of ANS and revealed the detailed ANS-binding mechanism. Furthermore, four water molecules were clearly identified in the binding cavity. Through these water molecules, the bound ANS molecule forms indirect hydrogen-bond interactions with FABP3. The adipocyte-type fatty-acid-binding protein (FABP4) exhibits 67% sequence identity with FABP3 and its crystal structure is almost the same as that of FABP3. However, FABP4 can bind with a higher affinity to ANS than FABP3. To understand the difference in their ligand specificities, a structural comparison was performed between FABP3-ANS and FABP4-ANS complexes. The result revealed that the orientation of ANS binding to FABP3 is completely opposite to that of ANS binding to FABP4, and the substitution of valine in FABP4 to leucine in FABP3 may result in greater steric hindrance between the side-chain of Leu115 and the aniline ring of ANS.

Adipocyte fatty acid binding protein: a novel adipokine involved in the pathogenesis of metabolic and vascular disease?

Adipocyte fatty acid binding protein (AFABP, also known as aP2 and FABP4) has recently been introduced as a novel fat-derived circulating protein. AFABP serum concentrations are positively correlated with markers of the metabolic syndrome and vascular disease in various cross-sectional and interventional studies. Furthermore, a small set of prospective studies indicates that high AFABP serum levels at baseline predict the risk for metabolic and vascular morbidity and mortality. Studies in Afabp (also known as Fabp4) knockout mice and AFABP inhibitor-treated animals suggest that total AFABP promotes insulin resistance, hypertriacylglycerolaemia and atherosclerosis by ligand/ligand delivery, as well as ligand-independent mechanisms. In contrast, the pathophysiological significance of circulating AFABP and the mechanisms leading to its release remain to be established. The current review summarises recent findings on the regulation and potential role of AFABP in metabolic and vascular disease.

The human fatty acid-binding protein family: evolutionary divergences and functions

Fatty acid-binding proteins (FABPs) are members of the intracellular lipid-binding protein (iLBP) family and are involved in reversibly binding intracellular hydrophobic ligands and trafficking them throughout cellular compartments, including the peroxisomes, mitochondria, endoplasmic reticulum and nucleus. FABPs are small, structurally conserved cytosolic proteins consisting of a water-filled, interior-binding pocket surrounded by ten anti-parallel beta sheets, forming a beta barrel. At the superior surface, two alpha-helices cap the pocket and are thought to regulate binding. FABPs have broad specificity, including the ability to bind long-chain (C16-C20) fatty acids, eicosanoids, bile salts and peroxisome proliferators. FABPs demonstrate strong evolutionary conservation and are present in a spectrum of species including Drosophila melanogaster, Caenorhabditis elegans, mouse and human. The human genome consists of nine putatively functional protein-coding FABP genes. The most recently identified family member, FABP12, has been less studied.

X-ray crystallographic analysis of adipocyte fatty acid binding protein (aP2) modified with 4-hydroxy-2-nonenal

Fatty acid binding proteins (FABP) have been characterized as facilitating the intracellular solubilization and transport of long-chain fatty acyl carboxylates via noncovalent interactions. More recent work has shown that the adipocyte FABP is also covalently modified in vivo on Cys117 with 4-hydroxy-2-nonenal (4-HNE), a bioactive aldehyde linked to oxidative stress and inflammation. To evaluate 4-HNE binding and modification, the crystal structures of adipocyte FABP covalently and noncovalently bound to 4-HNE have been solved to 1.9 A and 2.3 A resolution, respectively. While the 4-HNE in the noncovalently modified protein is coordinated similarly to a carboxylate of a fatty acid, the covalent form show a novel coordination through a water molecule at the polar end of the lipid. Other defining features between the two structures with 4-HNE and previously solved structures of the protein include a peptide flip between residues Ala36 and Lys37 and the rotation of the side chain of Phe57 into its closed conformation. Representing the first structure of an endogenous target protein covalently modified by 4-HNE, these results define a new class of in vivo ligands for FABPs and extend their physiological substrates to include bioactive aldehydes.

Identification and characterization of a small molecule inhibitor of Fatty Acid binding proteins

Molecular disruption of the lipid carrier AFABP/aP2 in mice results in improved insulin sensitivity and protection from atherosclerosis. Because small molecule inhibitors may be efficacious in defining the mechanism(s) of AFABP/aP2 action, a chemical library was screened and identified 1 (HTS01037) as a pharmacologic ligand capable of displacing the fluorophore 1-anilinonaphthalene 8-sulfonic acid from the lipid binding cavity. The X-ray crystal structure of 1 bound to AFABP/aP2 revealed that the ligand binds at a structurally similar position to a long-chain fatty acid. Similar to AFABP/aP2 knockout mice, 1 inhibits lipolysis in 3T3-L1 adipocytes and reduces LPS-stimulated inflammation in cultured macrophages. 1 acts as an antagonist of the protein-protein interaction between AFABP/aP2 and hormone sensitive lipase but does not activate PPARgamma in macrophage or CV-1 cells. These results identify 1 as an inhibitor of fatty acid binding and a competitive antagonist of protein-protein interactions mediated by AFABP/aP2.

Binding of small-molecule ligands to proteins: "what you see" is not always "what you get"

We review insights from computational studies of affinities of ligands binding to proteins. The power of structural biology is in translating knowledge of protein structures into insights about their forces, binding, and mechanisms. However, the complementary power of computer modeling is in showing "the rest of the story" (i.e., how motions and ensembles and alternative conformers and the entropies and forces that cannot be seen in single molecular structures also contribute to binding affinities). Upon binding to a protein, a ligand can bind in multiple orientations; the protein or ligand can be deformed by the binding event; waters, ions, or cofactors can have unexpected involvement; and conformational or solvation entropies can sometimes play large and otherwise unpredictable roles. Computer modeling is helping to elucidate these factors.

Molecular dynamics simulations of palmitate entry into the hydrophobic pocket of the fatty acid binding protein

The entry of substrate into the active site is the first event in any enzymatic reaction. However, due to the short time interval between the encounter and the formation of the stable complex, the detailed steps are experimentally unobserved. In the present study, we report a molecular dynamics simulation of the encounter between palmitate molecule and the Toad Liver fatty acid binding protein, ending with the formation of a stable complex resemblance in structure of other proteins of this family. The forces operating on the system leading to the formation of the tight complex are discussed.

The consequences of scoring docked ligand conformations using free energy correlations

Ligands from a set of 19 protein-ligand complexes were re-docked with AutoDock, GOLD and FlexX using the scoring algorithms native to these programs supplemented by analysis using the HINT free energy force field. A HINT scoring function was calibrated for this data set using a simple linear regression of total HINT score for crystal-structure complexes vs. measured free energy of binding. This function had an r(2) of 0.84 and a standard error of +/-0.42 kcal mol(-1). The free energies of binding were calculated for the best poses using the AutoDock, GOLD and FlexX scoring functions. The AutoDock and GoldScore algorithms estimated more than half of the binding free energies within the reported calibration standard errors for these functions, while that of FlexX did not. In contrast, the calibrated HINT scoring function identified optimized poses with standard errors near +/-0.5 kcal mol(-1). When the metric of success is minimum RMSD (vs. crystallographic coordinates) the three docking programs were more successful, with mean RMSDs for the top-ranking poses in the 19 complexes of 3.38, 2.52 and 2.62 A for AutoDock, GOLD and FlexX, respectively. Two key observations in this study have general relevance for computational medicinal chemistry: first, while optimizing RMSD with docking score functions is clearly of value, these functions may be less well optimized for free energy of binding, which has broader applicability in virtual screening and drug discovery than RMSD; second, scoring functions uniquely calibrated for the data set or sets under study should nearly always be preferable to universal scoring functions. Due to these advantages, the poses selected by the HINT score also required less post-docking structure optimization to produce usable molecular models. Most of these features may be achievable with other scoring functions.

Carbonylation of adipose proteins in obesity and insulin resistance: identification of adipocyte fatty acid-binding protein as a cellular target of 4-hydroxynonenal

Obesity is a state of mild inflammation correlated with increased oxidative stress. In general, pro-oxidative conditions lead to production of reactive aldehydes such as trans-4-hydroxy-2-nonenal (4-HNE) and trans-4-oxo-2-nonenal implicated in the development of a variety of metabolic diseases. To investigate protein modification by 4-HNE as a consequence of obesity and its potential relationship to the development of insulin resistance, proteomics technologies were utilized to identify aldehyde-modified proteins in adipose tissue. Adipose proteins from lean insulin-sensitive and obese insulin-resistant C57Bl/6J mice were incubated with biotin hydrazide and detected using horseradish peroxidase-conjugated streptavidin. High carbohydrate, high fat feeding of mice resulted in a approximately 2-3-fold increase in total adipose protein carbonylation. Consistent with an increase in oxidative stress in obesity, the abundance of glutathione S-transferase A4 (GSTA4), a key enzyme responsible for metabolizing 4-HNE, was decreased approximately 3-4-fold in adipose tissue of obese mice. To identify specific carbonylated proteins, biotin hydrazide-modified adipose proteins from obese mice were captured using avidin-Sepharose affinity chromatography, proteolytically digested, and subjected to LC-ESI MS/MS. Interestingly enzymes involved in cellular stress response, lipotoxicity, and insulin signaling such as glutathione S-transferase M1, peroxiredoxin 1, glutathione peroxidase 1, eukaryotic elongation factor 1alpha-1 (eEF1alpha1), and filamin A were identified. The adipocyte fatty acid-binding protein, a protein implicated in the regulation of insulin resistance, was found to be carbonylated in vivo with 4-HNE. In vitro modification of adipocyte fatty acid-binding protein with 4-HNE was mapped to Cys-117, occurred equivalently using either the R or S enantiomer of 4-HNE, and reduced the affinity of the protein for fatty acids approximately 10-fold. These results indicate that obesity is accompanied by an increase in the carbonylation of a number of adipose-regulatory proteins that may serve as a mechanistic link between increased oxidative stress and the development of insulin resistance.

Potent and selective biphenyl azole inhibitors of adipocyte fatty acid binding protein (aFABP)

Herein we report the first disclosure of biphenyl azoles that are nanomolar binders of adipocyte fatty acid binding protein (aFABP or aP2) with up to thousand-fold selectivity against muscle fatty acid binding protein and epidermal fatty acid binding protein. In addition a new radio-ligand to determine binding against the three fatty acid binding proteins was also synthesized.

Expression, purification, crystallization and structure of human adipocyte lipid-binding protein (aP2)

Human adipocyte lipid-binding protein (aP2) belongs to a family of intracellular lipid-binding proteins involved in the transport and storage of lipids. Here, the crystal structure of human aP2 with a bound palmitate is described at 1.5 A resolution. Unlike the known crystal structure of murine aP2 in complex with palmitate, this structure shows that the fatty acid is in a folded conformation and that the loop containing Phe57 acts as a lid to regulate ligand binding by excluding solvent exposure to the central binding cavity.

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.

Fully quantum mechanical energy optimization for protein-ligand structure

We present a quantum mechanical approach to study protein-ligand binding structure with application to a Adipocyte lipid-binding protein complexed with Propanoic Acid. The present approach employs a recently develop molecular fractionation with a conjugate caps (MFCC) method to compute protein-ligand interaction energy and performs energy optimization using the quasi-Newton method. The MFCC method enables us to compute fully quantum mechanical ab initio protein-ligand interaction energy and its gradients that are used in energy minimization. This quantum optimization approach is applied to study the Adipocyte lipid-binding protein complexed with Propanoic Acid system, a complex system consisting of a 2057-atom protein and a 10-atom ligand. The MFCC calculation is carried out at the Hartree-Fock level with a 3-21G basis set. The quantum optimized structure of this complex is in good agreement with the experimental crystal structure. The quantum energy calculation is implemented in a parallel program that dramatically speeds up the MFCC calculation for the protein-ligand system. Similarly good agreement between MFCC optimized structure and the experimental structure is also obtained for the streptavidin-biotin complex. Due to heavy computational cost, the quantum energy minimization is carried out in a six-dimensional space that corresponds to the rigid-body protein-ligand interaction.

Specificity determinants for lipids bound to beta-barrel proteins

The family of proteins accountable for the intracellular movement of lipids is characterized by a 10-stranded beta-barrel that forms an internalized cavity varying in size and binding preferences. The loop connecting beta-strands E and F (the fifth and sixth strands) is the most striking conformational difference between adipocyte lipid binding protein (ALBP; fatty acids) and cellular retinoic acid binding protein type I (CRABP I). A three-residue mutation was made in wild-type (WT)-ALBP [ALBP with a three-residue mutation (EF-ALBP)] to mimic CRABP I. Crystal structures of ligand-free and EF-ALBP with bound oleic acid were solved to resolutions of 1.5 A and 1.7 A, respectively, and compared with previous studies of WT-ALBP. The changes in three residues of one loop of the protein appear to have altered the positioning of the C18 fatty acid, as observed in the electron density of EF-ALBP. The crystallographic studies made it possible to compare the protein conformation and ligand positioning with those found in the WT protein. Although the cavity binding sites in both the retinoid and fatty acid binding proteins are irregular, the ligand atoms appear to favor a relatively planar region of the cavities. Preliminary chemical characterization of the mutant protein indicated changes in some binding properties and overall protein stability.

The crystal structure of Echinococcus granulosus fatty-acid-binding protein 1

We describe the 1.6 A crystal structure of the fatty-acid-binding protein EgFABP1 from the parasitic platyhelminth Echinococcus granulosus. E. granulosus causes hydatid disease, which is a major zoonosis. EgFABP1 has been implicated in the acquisition, storage, and transport of lipids, and may be important to the organism since it is incapable of synthesising most of its lipids de novo. Moreover, EgFABP1 is a promising candidate for a vaccine against hydatid disease. The crystal structure reveals that EgFABP1 has the expected 10-stranded beta-barrel fold typical of the family of intracellular lipid-binding proteins, and that it is structurally most similar to P2 myelin protein. We describe the comparison of the crystal structure of EgFABP1 with these proteins and with an older homology model for EgFABP1. The electron density reveals the presence of a bound ligand inside the cavity, which we have interpreted as palmitic acid. The carboxylate group of the fatty acid interacts with the protein's P2 motif, consisting of a conserved triad R em leader R-x-Y. The hydrophobic tail of the ligand assumes a fairly flat, U-shaped conformation and has relatively few interactions with the protein.We discuss some of the structural implications of the crystal structure of EgFABP1 for related platyhelminthic FABPs.

New insights into intracellular lipid binding proteins: The role of buried water

The crystal structures of most intracellular lipid binding proteins (LBPs) show between 5 and 20 internally bound water molecules, depending on the presence or the absence of ligand inside the protein cavity. The structural and functional significance of these waters has been discussed for several LBPs based on studies that used various biophysical techniques. The present work focuses on two very different LBPs, heart-type fatty acid binding protein (H-FABP) and ileal lipid binding protein (ILBP). Using high-resolution nuclear magnetic resonance spectroscopy, certain resonances belonging to side-chain protons that are located inside the water-filled lipid binding cavity were observed. In the case of H-FABP, the pH- and temperature-dependent behavior of selected side-chain resonances (Ser82 OgH and the imidazole ring protons of His93) indicated an unusually slow exchange with the solvent, implying that the intricate hydrogen-bonding network of amino-acid side-chains and water molecules in the protein interior is very rigid. In addition, holo H-FABP appeared to display a reversible self-aggregation at physiological pH. For ILBP, on the other hand, a more solvent-accessible protein cavity was deduced based on the pH titration behavior of its histidine residues. Comparison with data from other LBPs implies that the evolutionary specialization of LBPs for certain ligand types was not only because of mutations of residues directly involved in ligand binding but also to a refinement of the internal water scaffold.

Structure and backbone dynamics of Apo- and holo-cellular retinol-binding protein in solution

Retinoid-binding proteins play an important role in regulating transport, storage, and metabolism of vitamin A and its derivatives. The solution structure and backbone dynamics of rat cellular retinol-binding protein type I (CRBP) in the apo- and holo-form have been determined and compared using multidimensional high resolution NMR spectroscopy. The global fold of the protein is consistent with the common motif described for members of the intracellular lipid-binding protein family. The most relevant difference between the NMR structure ensembles of apo- and holoCRBP is the higher backbone disorder, in the ligand-free form, of some segments that frame the putative entrance to the ligand-binding site. These comprise alpha-helix II, the subsequent linker to beta-strand B, the hairpin turn between beta-strands C and D, and the betaE-betaF turn. The internal backbone dynamics, obtained from 15N relaxation data (T1, T2, and heteronuclear nuclear Overhauser effect) at two different fields, indicate several regions with significantly higher backbone mobility in the apoprotein, including the betaC-betaD and betaE-betaF turns. Although apoCRBP contains a binding cavity more shielded than that of any other retinoid carrier, conformational flexibility in the portal region may assist retinol uptake. The stiffening of the backbone in the holoprotein guarantees the stability of the complex during retinol transport and suggests that targeted retinol release requires a transiently open state that is likely to be promoted by the acceptor or the local environment.

Simple, intuitive calculations of free energy of binding for protein-ligand complexes. 1. Models without explicit constrained water

The prediction of the binding affinity between a protein and ligands is one of the most challenging issues for computational biochemistry and drug discovery. While the enthalpic contribution to binding is routinely available with molecular mechanics methods, the entropic contribution is more difficult to estimate. We describe and apply a relatively simple and intuitive calculation procedure for estimating the free energy of binding for 53 protein-ligand complexes formed by 17 proteins of known three-dimensional structure and characterized by different active site polarity. HINT, a software model based on experimental LogP(o/w) values for small organic molecules, was used to evaluate and score all atom-atom hydropathic interactions between the protein and the ligands. These total scores (H(TOTAL)), which have been previously shown to correlate with DeltaG(interaction) for protein-protein interactions, correlate with DeltaG(binding) for protein-ligand complexes in the present study with a standard error of +/-2.6 kcal mol(-1) from the equation DeltaG(binding) = -0.001 95 H(TOTAL) - 5.543. A more sophisticated model, utilizing categorized (by interaction class) HINT scores, produces a superior standard error of +/-1.8 kcal mol(-1). It is shown that within families of ligands for the same protein binding site, better models can be obtained with standard errors approaching +/-1.0 kcal mol(-1). Standardized methods for preparing crystallographic models for hydropathic analysis are also described. Particular attention is paid to the relationship between the ionization state of the ligands and the pH conditions under which the binding measurements are made. Sources and potential remedies of experimental and modeling errors affecting prediction of DeltaG(binding) are discussed.

One site fits both: a model for the ternary complex of folate + NADPH in R67 dihydrofolate reductase, a D2 symmetric enzyme

R67 dihydrofolate reductase (DHFR) is a novel enzyme that confers resistance to the antibiotic trimethoprim. The crystal structure of R67 DHFR displays a toroidal structure with a central active-site pore. This homotetrameric protein exhibits 222 symmetry, with only a few residues from each chain contributing to the active site, so related sites must be used to bind both substrate (dihydrofolate) and cofactor (NADPH) in the productive R67 DHFR.NADPH.dihydrofolate complex. Whereas the site of folate binding has been partially resolved crystallographically, an interesting question remains: how can the highly symmetrical active site also bind and orient NADPH for catalysis? To model this ternary complex, we employed DOCK and SLIDE, two methods for docking flexible ligands into proteins using quite different algorithms. The bound pteridine ring of folate (Fol I) from the crystal structure of R67 DHFR was used as the basis for docking the nicotinamide-ribose-Pi (NMN) moiety of NADPH. NMN was positioned by both DOCK and SLIDE on the opposite side of the pore from Fol I, where it interacts with Fol I at the pore's center. Numerous residues serve dual roles in binding. For example, Gln 67 from both the B and D subunits has several contacts with the pteridine ring, while the same residue from the A and C subunits has several contacts with the nicotinamide ring. The residues involved in dual roles are generally amphipathic, allowing them to make both hydrophobic and hydrophilic contacts with the ligands. The result is a 'hot spot' binding surface allowing the same residues to co-optimize the binding of two ligands, and orient them for catalysis.

Ligand-receptor docking with the Mining Minima optimizer

The optimizer developed for the Mining Minima algorithm, which uses ideas from Genetic Algorithms, the Global Underestimator Method, and Poling, has been adapted for use in ligand-receptor docking. The present study describes the resulting methodology and evaluates its accuracy and speed for 27 test systems. The performance of the new docking algorithm appears to be competitive with that of previously published methods. The energy model, an empirical force field with a distance-dependent dielectric treatment of solvation, is adequate for a number of test cases, although incorrect low-energy conformations begin to compete with the correct conformation for larger sampling volumes and for highly solvent-exposed binding sites that impose little steric constraint on the ligand.

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.

Studies of the ligand binding reaction of adipocyte lipid binding protein using the fluorescent probe 1, 8-anilinonaphthalene-8-sulfonate

The fluorescent probe anilinonaphthalene-8-sulfonate binds to adipocyte lipid binding protein at a site that competes with normal physiological ligands, such as fatty acids. Binding to the protein is accompanied by a relatively large increase in fluorescent intensity. To correlate the major change in optical properties and to determine the mechanism of competitive inhibition with fatty acids, the crystal structure of the protein with the bound fluorophore has been determined. In addition, the thermodynamic contributions to the binding reaction have been studied by titration calorimetry. Because the binding site is in a relatively internal position, kinetic studies have also been carried out to determine k(on). The results indicate that binding is not accompanied by any major conformational change. However, the negatively charged sulfonate moiety is not positioned the same as the carboxylate of fatty acid ligands as determined in previous studies. Nonetheless, the binding reaction is still driven by enthalpic effects. As judged by the crystallographic structure, a significant amount of the surface of the fluorophore is no longer exposed to water in the bound state.

Surface properties of adipocyte lipid-binding protein: Response to lipid binding, and comparison with homologous proteins

Adipocyte lipid-binding protein (ALBP) is one of a family of intracellular lipid-binding proteins (iLBPs) that bind fatty acids, retinoids, and other hydrophobic ligands. The different members of this family exhibit a highly conserved three-dimensional structure; and where structures have been determined both with (holo) and without (apo) bound lipid, observed conformational changes are extremely small (Banaszak, et al., 1994, Adv. Prot. Chem. 45, 89; Bernlohr, et al., 1997, Annu. Rev. Nutr. 17, 277). We have examined the electrostatic, hydrophobic, and water accessible surfaces of ALBP in the apo form and of holo forms with a variety of bound ligands. These calculations reveal a number of previously unrecognized changes between apo and holo ALBP, including: 1) an increase in the overall protein surface area when ligand binds, 2) expansion of the binding cavity when ligand is bound, 3) clustering of individual residue exposure increases in the area surrounding the proposed ligand entry portal, and 4) ligand-binding dependent variation in the topology of the electrostatic potential in the area surrounding the ligand entry portal. These focused analyses of the crystallographic structures thus reveal a number of subtle but consistent conformational and surface changes that might serve as markers for differential targeting of protein-lipid complexes within the cell. Most changes are consistent from ligand to ligand, however there are some ligand-specific changes. Comparable calculations with intestinal fatty-acid-binding protein and other vertebrate iLBPs show differences in the electrostatic topology, hydrophobic topology, and in localized changes in solvent exposure near the ligand entry portal. These results provide a basis toward understanding the functional and mechanistic differences among these highly structurally homologous proteins. Further, they suggest that iLBPs from different tissues exhibit one of two predominant end-state structural distributions of the ligand entry portal.

Thermodynamics of fatty acid binding to engineered mutants of the adipocyte and intestinal fatty acid-binding proteins

We constructed 18 single amino acid mutants of the adipocyte fatty acid-binding protein (A-FABP) and 17 of the intestinal fatty acid-binding protein (I-FABP), at locations in the fatty acid (FA) binding sites. For each mutant protein, we measured thermodynamic parameters that characterize FA binding. Binding affinities ranged from about 200-fold smaller to 30-fold larger than the wild type (WT) proteins. Thermodynamic parameters revealed that binding affinities often inaccurately reported changes in protein-FA interactions because changes in the binding entropy and enthalpy were usually compensatory and larger than the binding free energy. FA-FABP interactions were quite different for I-FABP and A-FABP proteins. Binding affinities were larger and decreased to a greater degree with increasing FA solubility for most of the I-FABP as compared with the A-FABP proteins, consistent with a more hydrophobic binding site for the I-FABP proteins. In A-FABP, Ala substitutions for Arg106 and Arg126, which interact with the FA carboxylate, reduce affinities by about 100-fold, but in I-FABP, R106A increases affinities up to 30-fold. Moreover, in A-FABP, the thermodynamic parameters predict that the FA carboxylate location switches from the 126-position in R106A to the 106 position in R126A. Finally, the A-FABP proteins, in contrast to the I-FABP proteins, reveal significant heat capacity changes (DeltaCp) upon FA binding, and substitutions at residues Arg106 and Arg126 reduce the magnitude of DeltaCp.

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.

Mutants of rat intestinal fatty acid-binding protein illustrate the critical role played by enthalpy-entropy compensation in ligand binding

Site-specific variants of rat intestinal fatty acid-binding protein were constructed to identify the molecular interactions that are important for binding to fatty acids (FAs). Several variants displayed affinities that appeared incompatible with the crystal structure of the protein-FA complex. Thermodynamic measurements provided an explanation for these apparent inconsistencies and revealed that binding affinities often inaccurately reported changes in protein-FA interactions because changes in the binding entropy and enthalpy were usually compensatory. These results demonstrate that understanding the effects of amino acid replacements on ligand binding requires measurements of enthalpy and entropy, in addition to affinity.

Biochemical and crystallographic analyses of a portal mutant of the adipocyte lipid-binding protein

A number of crystallographic studies of the adipocyte lipid-binding protein have established that the fatty acid-binding site is within an internalized water-filled cavity. The same studies have also suggested the existence of a region physically distinct from the fatty acid-binding site which connects the cavity of the protein with the external solvent, hereafter referred to as the portal. In an effort to examine the portal region, we have used site-directed mutagenesis to introduce the mutations V32D/F57H into the murine ALBP cDNA. Mutant protein has been isolated, crystallized, and its stability and binding properties studied by biochemical methods. As assessed by guanidine-HCl denaturation, the mutant form exhibited a slight overall destabilization relative to the wild-type protein under both acid and alkaline conditions. Accessibility to the cavity in both the mutant and wild-type proteins was observed by stopped-flow analysis of the modification of a cavity residue, Cys117, by the sulfhydryl reactive agent 5, 5'-dithiobis(2-nitrobenzoic acid) at pH 8.5. Cys117 of V32D/F57H ALBP was modified 7-fold faster than the wild-type protein. The ligand binding properties of both the V32D/F57H mutant and wild-type proteins were analyzed using a fluorescent probe at pH 6.0 and 8.0. The apparent dissociation constants for 1-anilinonaphthalene-8-sulfonic acid were approximately 9-10-fold greater than the wild-type protein, independent of pH. In addition, there is a 6-fold increase in the Kd for oleic acid for the portal mutant relative to the wild-type at pH 8.0. To study the effect of pH on the double mutant, it was crystallized and analyzed in two distinct space groups at pH 4.5 and 6.4. While in general the differences in the overall main chain conformations are negligible, changes were observed in the crystallographic structures near the site of the mutations. At both pH values, the mutant side chains are positioned somewhat differently than in wild-type protein. To ensure that the mutations had not altered ionic conditions near the binding site, the crystallographic coordinates were used to monitor the electrostatic potentials from the head group site to the positions near the portal region. The differences in the electrostatic potentials were small in all regions, and did not explain the differences in ligand affinity. We present these results within the context of fatty acid binding and suggest lipid association is more complex than that described within a single equilibrium event.

Development and validation of a genetic algorithm for flexible docking

Prediction of small molecule binding modes to macromolecules of known three-dimensional structure is a problem of paramount importance in rational drug design (the "docking" problem). We report the development and validation of the program GOLD (Genetic Optimisation for Ligand Docking). GOLD is an automated ligand docking program that uses a genetic algorithm to explore the full range of ligand conformational flexibility with partial flexibility of the protein, and satisfies the fundamental requirement that the ligand must displace loosely bound water on binding. Numerous enhancements and modifications have been applied to the original technique resulting in a substantial increase in the reliability and the applicability of the algorithm. The advanced algorithm has been tested on a dataset of 100 complexes extracted from the Brookhaven Protein DataBank. When used to dock the ligand back into the binding site, GOLD achieved a 71% success rate in identifying the experimental binding mode.

The crystal structure of the liver fatty acid-binding protein. A complex with two bound oleates

The crystal structure of the recombinant form of rat liver fatty acid-binding protein was completed to 2.3 A and refined to an R factor of 19.0%. The structural solution was obtained by molecular replacement using superimposed polyalanine coordinates of six intracellular lipid-binding proteins as a search probe. The entire amino acid sequence of rat liver fatty acid-binding protein along with an amino-terminal formyl-methionine was modeled in the crystal structure. In addition, the crystal was obtained in the presence of oleic acid, and the initial electron density clearly showed two fatty acid molecules bound within a central cavity. The carboxylate of one fatty acid molecule interacts with arginine 122 and is shielded from free solvent. It has an overall bent conformation. The more solvent-exposed carboxylate of the other oleate is located near the helix-turn-helix that caps one end of the beta-barrel, while the acyl chain lies in the interior. The cavity contains both polar and nonpolar residues but also shows extensive hydrophobic character around the nonpolar atoms of the ligands. The primary and secondary oleate binding sites appear to be totally interdependent, mainly because favorable hydrophobic interactions form between both aliphatic chains.

Intracellular lipid-binding proteins and their genes

Intracellular lipid-binding proteins are a family of low-molecular-weight single-chain polypeptides that form 1:1 complexes with fatty acids, retinoids, or other hydrophobic ligands. These proteins are products of a large multigene family of unlinked loci distributed throughout the genome. Each lipid-binding protein exhibits a distinctive pattern of tissue distribution. Transcriptional control, regulated by a combination of peroxisome proliferator activated receptors and CCAAT/enhancer-binding proteins, allows for a variety of both cell and tissue-specific expression patterns. In some cells, fatty acids increase the expression of the lipid-binding protein genes. Fatty acids, or their metabolites, are activators of the peroxisome proliferator-activated receptor family of transcription factors. Therefore, as the concentration of lipid in the diet increases, the expression of lipid-binding proteins coordinately increases. As revealed by X-ray crystallography, the lipid-binding proteins fold into beta-barrels, forming a large internal water-filled cavity. Fatty acid ligands are bound within the cavity, occupying only about one-third of the accessible volume. The bound fatty acid is stabilized via a combination of enthalpic and entropic forces that govern ligand affinity and selectivity. Cytoplasmic lipid-binding proteins are the intracellular receptors for hydrophobic ligands, delivering them to the appropriate site for use as metabolic fuels and regulatory agents.

Fatty acid binding protein is induced in neurons of the dorsal root ganglia after peripheral nerve injury

Peripheral nerve trauma induces the expression of genes presumed to be involved in the process of nerve degeneration and repair. In the present study, an in vivo paradigm was employed to identify molecules which may have important roles in these processes. A cDNA library was constructed with RNA extracted from rat dorsal root ganglia (DRG) 3 days after a sciatic nerve crush. After differential hybridization to this library, several cDNAs were identified that encoded mRNAs that were upregulated in the DRG ipsilateral to the crush injury, as opposed to the contralateral or naive DRG. Approximately 0.15% of all the clones screened were found to be induced. This report presents the types of induced sequences identified and characterizes one of them, DA11. The 0.7 kb DA11 full length cDNA clone contains a 405 nucleotide open reading frame that encodes a putative protein of 15.2 kDa (135 amino acid residues) and is a member of the family of fatty acid binding proteins (FABP). The DA11 protein differs by one amino acid residue from the sequence of the C-FAPB protein and by eight residues from the sequence of mal1, proteins found in rat and mouse skin, respectively. Northern and Western blot analyses showed that the DA11 mRNA and protein were induced in the injured DRG. Furthermore, studies using antibodies generated against DA11 found that the DA11-like immunoreactivity was more pronounced in the nuclei of neurons located in the DRG ipsilateral to the sciatic cut than those located in the contralateral DRG. The induction of DA11 mRNA and protein in DRG neurons suggests, for the first time, the involvement of a neuronal FABP in the process of degeneration and repair in the nervous system.

Fatty acid binding and conformational stability of mutants of human muscle fatty acid-binding protein

Human muscle fatty acid-binding protein (M-FABP) is a 15 kDa cytosolic protein which may be involved in fatty acid transfer and modulation of non-esterified fatty acid concentration in heart, skeletal muscle, kidney and many other tissues. Crystallographic studies have suggested the importance of the amino acids Thr-40, Arg-106, Arg-126 and Tyr-128 for the hydrogen bonding network of the fatty acid carboxylate group. Two phenylalanines at 16 and 57 are positioned to interact with the acyl chain of the fatty acid. We prepared 13 mutant proteins by site-directed mutagenesis and tested them for fatty acid binding and stability. Substitution of amino acids Phe-16, Arg-106 or Arg-126 created proteins which showed a large decrease in or complete loss of oleic acid binding. Substitution of Phe-57 by Ser or Val and of Tyr-128 by Phe had no great effect. The stability of the mutant proteins was tested by denaturation studies on the basis of fatty acid binding or tryptophan fluorescence and compared with that of the wild-type M-FABP. There was no direct relationship between fatty acid-binding activity and stability. Less stable mutants (F57S and Y128F) did not show a marked change in fatty acid-binding activity. Substitution of Arg-126 by Gln or Arg-106 by Thr eliminated binding activity, but the former mutant protein showed wild-type stability, in contrast to the latter. The results are in agreement with crystallographic data.

Hydration in drug design. 2. Influence of local site surface shape on water binding

If water molecules are strongly bound at a protein-ligand interface, they are unlikely to be displaced during ligand binding. Such water molecules can change the shape of the ligand binding site and thus affect strategies for drug design. To understand the nature of water binding, and factors influencing it, water molecules at the ligand binding sites of 26 high-resolution protein-ligand complexes have been examined here. Water molecules bound in deep grooves and cavities between the protein and the ligand are located in the indentations on the protein-site surface, but not in the indentations on the ligand surface. The majority of the water molecules bound in deep indentations on the protein-site surface make multiple polar contacts with the protein surface. This may indicate a strong binding of water molecules in deep indentations on protein-site surfaces. The local shape of the site surface may influence the binding of water molecules that mediate protein-ligand interactions.

Hydration in drug design. 3. Conserved water molecules at the ligand-binding sites of homologous proteins

Water molecules are known to play an important rôle in mediating protein-ligand interactions. If water molecules are conserved at the ligand-binding sites of homologous proteins, such a finding may suggest the structural importance of water molecules in ligand binding. Structurally conserved water molecules change the conventional definition of 'binding sites' by changing the shape and complementarity of these sites. Such conserved water molecules can be important for site-directed ligand/drug design. Therefore, five different sets of homologous protein/protein-ligand complexes have been examined to identify the conserved water molecules at the ligand-binding sites. Our analysis reveals that there are as many as 16 conserved water molecules at the FAD binding site of glutathione reductase between the crystal structures obtained from human and E. coli. In the remaining four sets of high-resolution crystal structures, 2-4 water molecules have been found to be conserved at the ligand-binding sites. The majority of these conserved water molecules are either bound in deep grooves at the protein-ligand interface or completely buried in cavities between the protein and the ligand. All these water molecules, conserved between the protein/protein-ligand complexes from different species, have identical or similar apolar and polar interactions in a given set. The site residues interacting with the conserved water molecules at the ligand-binding sites have been found to be highly conserved among proteins from different species; they are more conserved compared to the other site residues interacting with the ligand. These water molecules, in general, make multiple polar contacts with protein-site residues.

Thermodynamics of water mediating protein-ligand interactions in cytochrome P450cam: a molecular dynamics study

Ordered water molecules are observed by crystallography and nuclear magnetic resonance to mediate protein-ligand interactions. Here, we examine the energetics of hydrating cavities formed at protein-ligand interfaces using molecular dynamics simulations. The free energies of hydrating two cavities in the active site of two liganded complexes of cytochrome P450cam were calculated by multiconfigurational thermodynamic integration. The complex of cytochrome P450cam with 2-phenyl-imidazole contains a crystallographically well defined water molecule mediating hydrogen bonds between the protein and the inhibitor. We calculate that this water molecule is stabilized by a binding free energy of -11.6 +/- kJ/mol. The complex of cytochrome P450cam with its natural substrate, camphor, contains a cavity that is empty in the crystal structure although a water molecule in it could make a hydrogen bond to camphor. Here, solvation of this cavity is calculated to be unfavorable by +15.8 +/- 5.0 kJ/mol. The molecular dynamics simulations can thus distinguish a hydrated interfacial cavity from an empty one. They also provide support for the notion that protein-ligand complexes can accommodate empty interfacial cavities and that such cavities are likely to be unhydrated unless more than one hydrogen bond can be made to a water molecule in the cavity.

Thermodynamics of fatty acid binding to fatty acid-binding proteins and fatty acid partition between water and membranes measured using the fluorescent probe ADIFAB

Using the fluorescent probe ADIFAB (acrylodan-derivatized intestinal fatty acid-binding protein) to determine the equilibrium concentration of the free (unbound) fatty acid (FFA), dissociation constants were measured between 10 and 50 degrees C for the interaction of five different long chain fatty acids (FA) with fatty acid-binding proteins (FABP) from adipocyte, intestine, and heart. Gibbs free energies (delta G) determined from the dissociation constants were between about -9 and -11 kcal/mol at 25 degrees C. Thermodynamic parameters for binding were determined using van't Hoff plots of the dissociation constants, which range, over the entire temperature region, between 2 and 3000 nM. For all the unlabeled FABPs, free energies of binding were dominated by large negative enthalpies that ranged from -7 to -12 kcal/mol, and the enthalpies tended to decrease with increasing FA unsaturation. The entropic contributions (-T detla S) at 25 degrees C between -4 and +2 kcal/mol and tended to increase with increasing FA unsaturation. To assess the role of FA aqueous solubility in FABP binding, measurements of the partition of FA between unilamellar lipid vesicles and water were also done using ADIFAB; the lipid/water partition coefficients (Kp) determined from these measurements were found to be independent of temperature. The binding of FA to FABP is governed by the sum of contributions of various interactions between FA, water, and FABP. An analysis of the individual contributions suggests that the net free energy of binding results from the canceling in part of a number of separate quite large contributions. The entropic contributions sum almost to zero for most FA and FABPs as a result of the canceling of a large increase in bulk solvent entropy by decreases in configurational entropy upon FA binding to FABP. The net, approximately -10 kcal/mol enthalpy of binding, probably results from an increase in FA configurational enthalpy upon binding to FABP plus a large negative enthalpy from the interaction between the FA and the FABP. This large enthalpy of the FA-FABP interaction suggests that in addition to previously identified specific interactions between the carboxylate portion of the FA and charged amino acids within the binding cavity, other significantly larger enthalpic interactions, presumably involving the hydrocarbon portion of the FA, must contribute to the binding energy.