Backbone and side chain dynamics of uncomplexed human adipocyte and muscle fatty acid-binding proteins

Investigating protein-membrane interactions using native reverse micelles constructed from naturally sourced lipids

Advancing the study of membrane associated proteins and their interactions is dependent on accurate membrane models. While a variety of membrane models for high-resolution membrane protein study exist, most do not reflect the diversity of lipids found within biological membranes. In this work, we have developed native reverse micelles (nRMs) formulated with lipids from multiple eukaryotic sources, which encapsulate proteins and enable them to interact as they would with a biological membrane. Diverse formulations of nRMs using soy lecithin, porcine brain lipids, or bovine heart lipids combined with n-dodecylphosphocholine were developed and characterized by dynamic light scattering and 31 P-NMR. To optimize protein encapsulation, ubiquitin was used as a standard and protein NMR verified minimal changes to its structure. Peripheral membrane proteins, which bind reversibly to membranes, were encapsulated and include glutathione peroxidase 4 (GPx4), phosphatidylethanolamine-binding protein 1 (PEBP1), and fatty acid binding protein 4 (FABP4). All three proteins showed anticipated interactions with the membrane-like inner surface of the nRMs as assessed by protein NMR. The nRM formulations developed here allow for efficient, high-resolution study of membrane interacting proteins up to and beyond ~21 kDa, in a more biologically relevant context compared to other non-native membrane models. The approach outlined here may be applied to a wide range of lipid extracts, allowing study of a variety of membrane associated proteins in their specific biological context.

Structural and Dynamic Determinants of Molecular Recognition in Bile Acid-Binding Proteins

Disorders in bile acid transport and metabolism have been related to a number of metabolic disease states, atherosclerosis, type-II diabetes, and cancer. Bile acid-binding proteins (BABPs), a subfamily of intracellular lipid-binding proteins (iLBPs), have a key role in the cellular trafficking and metabolic targeting of bile salts. Within the family of iLBPs, BABPs exhibit unique binding properties including positive binding cooperativity and site-selectivity, which in different tissues and organisms appears to be tailored to the local bile salt pool. Structural and biophysical studies of the past two decades have shed light on the mechanism of bile salt binding at the atomic level, providing us with a mechanistic picture of ligand entry and release, and the communication between the binding sites. In this review, we discuss the emerging view of bile salt recognition in intestinal- and liver-BABPs, with examples from both mammalian and non-mammalian species. The structural and dynamic determinants of the BABP-bile–salt interaction reviewed herein set the basis for the design and development of drug candidates targeting the transcellular traffic of bile salts in enterocytes and hepatocytes.

Dynamics of proteins in solution

The dynamics of proteins in solution includes a variety of processes, such as backbone and side-chain fluctuations, interdomain motions, as well as global rotational and translational (i.e. center of mass) diffusion. Since protein dynamics is related to protein function and essential transport processes, a detailed mechanistic understanding and monitoring of protein dynamics in solution is highly desirable. The hierarchical character of protein dynamics requires experimental tools addressing a broad range of time- and length scales. We discuss how different techniques contribute to a comprehensive picture of protein dynamics, and focus in particular on results from neutron spectroscopy. We outline the underlying principles and review available instrumentation as well as related analysis frameworks.

Protocols and pitfalls in obtaining fatty acid-binding proteins for biophysical studies of ligand-protein and protein-protein interactions

Adipocyte fatty acid-binding protein (AFABP: FABP4) is a member of the intracellular lipid-binding protein family that is thought to target long-chain fatty acids to nuclear receptors such as peroxisome proliferator-activated receptor gamma (PPARγ), which in turn plays roles in insulin resistance and obesity. A molecular understanding of AFABP function requires robust isolation of the protein in liganded and free forms as well as characterization of its oligomerization state(s) under physiological conditions. We report development of a protocol to optimize the production of members of this protein family in pure form, including removal of their bound lipids by mixing with hydrophobically functionalized hydroxypropyl dextran beads and validation by two-dimensional NMR spectroscopy. The formation of self-associated or covalently bonded protein dimers was evaluated critically using gel filtration chromatography, revealing conditions that promote or prevent formation of disulfide-linked homodimers. The resulting scheme provides a solid foundation for future investigations of AFABP interactions with key ligand and protein partners involved in lipid metabolism.

On the ability of molecular dynamics force fields to recapitulate NMR derived protein side chain order parameters

Molecular dynamics (MD) simulations have become a central tool for investigating various biophysical questions with atomistic detail. While many different proxies are used to qualify MD force fields, most are based on largely structural parameters such as the root mean square deviation from experimental coordinates or nuclear magnetic resonance (NMR) chemical shifts and residual dipolar couplings. NMR derived Lipari-Szabo squared generalized order parameter (O(2) ) values of amide NH bond vectors of the polypeptide chain were also often employed for refinement and validation. However, with a few exceptions, side chain methyl symmetry axis order parameters have not been incorporated into experimental reference sets. Using a test set of five diverse proteins, the performance of several force fields implemented in the NAMDD simulation package was examined. It was found that simulations employing explicit water implemented using the TIP3 model generally performed significantly better than those using implicit water in reproducing experimental methyl symmetry axis O(2) values. Overall the CHARMM27 force field performs nominally better than two implementations of the Amber force field. It appeared that recent quantum mechanics modifications to side chain torsional angles of leucine and isoleucine in the Amber force field have significantly hindered proper motional modeling for these residues. There remained significant room for improvement as even the best correlations of experimental and simulated methyl group Lipari-Szabo generalized order parameters fall below an R(2) of 0.8.

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

In this work, we investigate the dynamic motions of fatty acid binding protein 4 (FABP4) in the absence and presence of a ligand by explicitly solvated all-atom molecular dynamics simulations. The dynamics of one ligand-free FABP4 and four ligand-bound FABP4s is compared via multiple 1.2 μs simulations. In our simulations, the protein interconverts between the open and closed states. Ligand-free FABP4 prefers the closed state, whereas ligand binding induces a conformational transition to the open state. Coupled with opening and closing of FABP4, the ligand adopts distinct binding modes, which are identified and compared with crystal structures. The concerted dynamics of protein and ligand suggests that there may exist multiple FABP4–ligand binding conformations. Thus, this work provides details about how ligand binding affects the conformational preference of FABP4 and how ligand binding is coupled with a conformational change of FABP4 at an atomic level.

The role of dynamics in modulating ligand exchange in intracellular lipid binding proteins

Lipids are essential for many biological processes and crucial in the pathogenesis of several diseases. Intracellular lipid-binding proteins (iLBPs) provide mobile hydrophobic binding sites that allow hydrophobic or amphipathic lipid molecules to penetrate into and across aqueous layers. Thus iLBPs mediate the lipid transport within the cell and participate to a spectrum of tissue-specific pathways involved in lipid homeostasis. Structural studies have shown that iLBPs' binding sites are inaccessible from the bulk, implying that substrate binding should involve a conformational change able to produce a ligand entry portal. Many studies have been reported in the last two decades on iLBPs indicating that their dynamics play a pivotal role in regulating ligand binding and targeted release. The ensemble of reported data has not been reviewed until today. This review is thus intended to summarize and possibly generalize the results up to now described, providing a picture which could help to identify the missing notions necessary to improve our understanding of the role of dynamics in iLBPs' molecular recognition. Such notions would clarify the chemistry of lipid binding to iLBPs and set the basis for the development of new drugs.

Microscopic insights into the NMR relaxation-based protein conformational entropy meter

Conformational entropy is a potentially important thermodynamic parameter contributing to protein function. Quantitative measures of conformational entropy are necessary for an understanding of its role but have been difficult to obtain. An empirical method that utilizes changes in conformational dynamics as a proxy for changes in conformational entropy has recently been introduced. Here we probe the microscopic origins of the link between conformational dynamics and conformational entropy using molecular dynamics simulations. Simulation of seven proteins gave an excellent correlation with measures of side-chain motion derived from NMR relaxation. The simulations show that the motion of methyl-bearing side chains are sufficiently coupled to that of other side chains to serve as excellent reporters of the overall side-chain conformational entropy. These results tend to validate the use of experimentally accessible measures of methyl motion--the NMR-derived generalized order parameters--as a proxy from which to derive changes in protein conformational entropy.

A simple method to measure protein side-chain mobility using NMR chemical shifts

Protein side-chain motions are involved in many important biological processes including enzymatic catalysis, allosteric regulation, and the mediation of protein-protein, protein-DNA, protein-RNA, and protein-cofactor interactions. NMR spectroscopy has long been used to provide insights into the motions of side-chain groups. Currently, the method of choice for studying side-chain dynamics by NMR is the measurement of methyl (2)H autorelaxation. Methyl (2)H autorelaxation exhibits simple relaxation mechanisms and can be straightforwardly converted to meaningful dynamic parameters. However, methyl groups can only be found in 6 of 19 side-chain bearing amino acids. Consequently, only a sparse assessment of protein side-chain dynamics is possible. Therefore, there is a significant interest in developing novel methods of studying side-chain motions that can be applied to all types of side-chains. Here, we show how side-chain chemical shifts can be used to determine the magnitude of fast side-chain motions in proteins. The chemical shift method is applicable to all side-chain bearing residues and does not require any additional measurements beyond standard NMR experiments for backbone and side-chain assignments.

The description of protein internal motions aids selection of ligand binding poses by the INPHARMA method

Protein internal motions influence observables of NMR experiments. The effect of internal motions occurring at the sub-nanosecond timescale can be described by NMR order parameters. Here, we report that the use of order parameters derived from Molecular Dynamics (MD) simulations of two holo-structures of Protein Kinase A increase the discrimination power of INPHARMA, an NMR based methodology that selects docked ligand orientations by maximizing the correlation of back-calculated to experimental data. By including internal motion in the back-calculation of the INPHARMA transfer, we obtain a more realistic description of the system, which better represents the experimental data. Furthermore, we propose a set of generic order parameters, derived from MD simulations of globular proteins, which can be used in the back-calculation of INPHARMA NOEs for any protein-ligand complex, thus by-passing the need of obtaining system-specific order parameters for new protein-ligand complexes.

Solution structure and backbone dynamics of human liver fatty acid binding protein: fatty acid binding revisited

Liver fatty acid binding protein (L-FABP), a cytosolic protein most abundant in liver, is associated with intracellular transport of fatty acids, nuclear signaling, and regulation of intracellular lipolysis. Among the members of the intracellular lipid binding protein family, L-FABP is of particular interest as it can i), bind two fatty acid molecules simultaneously and ii), accommodate a variety of bulkier physiological ligands such as bilirubin and fatty acyl CoA. To better understand the promiscuous binding and transport properties of L-FABP, we investigated structure and dynamics of human L-FABP with and without bound ligands by means of heteronuclear NMR. The overall conformation of human L-FABP shows the typical β-clam motif. Binding of two oleic acid (OA) molecules does not alter the protein conformation substantially, but perturbs the chemical shift of certain backbone and side-chain protons that are involved in OA binding according to the structure of the human L-FABP/OA complex. Comparison of the human apo and holo L-FABP structures revealed no evidence for an "open-cap" conformation or a "swivel-back" mechanism of the K90 side chain upon ligand binding, as proposed for rat L-FABP. Instead, we postulate that the lipid binding process in L-FABP is associated with backbone dynamics.

Relative stability of the open and closed conformations of the active site loop of streptavidin

The eight-residue surface loop, 45-52 (Ser, Ala, Val, Gly, Asn, Ala, Glu, Ser), of the homotetrameric protein streptavidin has a "closed" conformation in the streptavidin-biotin complex, where the corresponding binding affinity is one of the strongest found in nature (ΔG ∼ -18 kcal∕mol). However, in most of the crystal structures of apo (unbound) streptavidin, the loop conformation is "open" and typically exhibits partial disorder and high B-factors. Thus, it is plausible to assume that the loop structure is changed from open to closed upon binding of biotin, and the corresponding difference in free energy, ΔF = F(open) - F(closed) in the unbound protein, should therefore be considered in the total absolute free energy of binding. ΔF (which has generally been neglected) is calculated here using our "hypothetical scanning molecular-dynamics" (HSMD) method. We use a protein model in which only the atoms closest to the loop are considered (the "template") and they are fixed in the x-ray coordinates of the free protein; the x-ray conformation of the closed loop is attached to the same (unbound) template and both systems are capped with the same sphere of TIP3P water. Using the force field of the assisted model building with energy refinement (AMBER), we carry out two separate MD simulations (at temperature T = 300 K), starting from the open and closed conformations, where only the atoms of the loop and water are allowed to move (the template-water and template-loop interactions are considered). The absolute F(open) and F(closed) (of loop + water) are calculated from these trajectories, where the loop and water contributions are obtained by HSMD and a thermodynamic integration (TI) process, respectively. The combined HSMD-TI procedure leads to total (loop + water) ΔF = -27.1 ± 2.0 kcal∕mol, where the entropy TΔS constitutes 34% of ΔF, meaning that the effect of S is significant and should not be ignored. Also, ΔS is positive, in accord with the high flexibility of the open loop observed in crystal structures, while the energy ΔE is unexpectedly negative, thus also adding to the stability of the open loop. The loop and the 250 capped water molecules are the largest system studied thus far, which constitutes a test for the efficiency of HSMD-TI; this efficiency and technical issues related to the implementation of the method are also discussed. Finally, the result for ΔF is a prediction that will be considered in the calculation of the absolute free energy of binding of biotin to streptavidin, which constitutes our next project.

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.

Entropy and Free Energy of a Mobile Loop Based on the Crystal Structures of the Free and Bound Proteins

A mobile loop changes its conformation from "open" (free enzyme) to "closed" upon ligand binding. The difference in the Helmholtz free energy, ΔF(loop) between these states sheds light on the mechanism of binding. With our "hypothetical scanning molecular dynamics" (HSMD-TI) method ΔF(loop) = F(free) - F(bound) where F(free) and F(bound) are calculated from two MD samples of the free and bound loop states; the contribution of water is obtained by a thermodynamic integration (TI) procedure. In previous work the free and bound loop structures were both attached to the same "template" which was "cut" from the crystal structure of the free protein. Our results for loop 287-290 of AcetylCholineEsterase agree with the experiment, ΔF(loop)~ -4 kcal/mol if the density of the TIP3P water molecules capping the loop is close to that of bulk water, i.e., N(water) = 140 - 180 waters in a sphere of a 18 Å radius. Here we calculate ΔF(loop) for the more realistic case, where two templates are "cut" from the crystal structures, 2dfp.pdb (bound) and 2ace.pdb (free), where N(water) = 40 - 160; this requires adding a computationally more demanding (second) TI procedure. While the results for N(water) ≤ 140 are computationally sound, ΔF(loop) is always positive (18 ± 2 kcal/mol for N(water) = 140). These (disagreeing) results are attributed to the large average B-factor, 41.6 of 2dfp (23.4 Å(2) for 2ace). While this conformational uncertainty is an inherent difficulty, the (unstable) results for N(water) = 160 suggest that it might be alleviated by applying different (initial) structural optimizations to each template.

Methods for calculating the absolute entropy and free energy of biological systems based on ideas from polymer physics

The commonly used simulation techniques, Metropolis Monte Carlo (MC) and molecular dynamics (MD) are of a dynamical type which enables one to sample system configurations i correctly with the Boltzmann probability, P(i)(B), while the value of P(i)(B) is not provided directly; therefore, it is difficult to obtain the absolute entropy, S approximately -ln P(i)(B), and the Helmholtz free energy, F. With a different simulation approach developed in polymer physics, a chain is grown step-by-step with transition probabilities (TPs), and thus their product is the value of the construction probability; therefore, the entropy is known. Because all exact simulation methods are equivalent, i.e. they lead to the same averages and fluctuations of physical properties, one can treat an MC or MD sample as if its members have rather been generated step-by-step. Thus, each configuration i of the sample can be reconstructed (from nothing) by calculating the TPs with which it could have been constructed. This idea applies also to bulk systems such as fluids or magnets. This approach has led earlier to the "local states" (LS) and the "hypothetical scanning" (HS) methods, which are approximate in nature. A recent development is the hypothetical scanning Monte Carlo (HSMC) (or molecular dynamics, HSMD) method which is based on stochastic TPs where all interactions are taken into account. In this respect, HSMC(D) can be viewed as exact and the only approximation involved is due to insufficient MC(MD) sampling for calculating the TPs. The validity of HSMC has been established by applying it first to liquid argon, TIP3P water, self-avoiding walks (SAW), and polyglycine models, where the results for F were found to agree with those obtained by other methods. Subsequently, HSMD was applied to mobile loops of the enzymes porcine pancreatic alpha-amylase and acetylcholinesterase in explicit water, where the difference in F between the bound and free states of the loop was calculated. Currently, HSMD is being extended for calculating the absolute and relative free energies of ligand-enzyme binding. We describe the whole approach and discuss future directions.

Entropy and free energy of a mobile protein loop in explicit water

Estimation of the energy from a given Boltzmann sample is straightforward since one just has to average the contribution of the individual configurations. On the other hand, calculation of the absolute entropy, S (hence the absolute free energy F) is difficult because it depends on the entire (unknown) ensemble. We have developed a new method called "the hypothetical scanning molecular dynamics" (HSMD) for calculating the absolute S from a given sample (generated by any simulation technique). In other words, S (like the energy) is "written" on the sample configurations, where HSMD provides a prescription of how to "read" it. In practice, each sample conformation, i, is reconstructed with transition probabilities, and their product leads to the probability of i, hence to the entropy. HSMD is an exact method where all interactions are considered, and the only approximation is due to insufficient sampling. In previous studies HSMD (and HS Monte CarloHSMC) has been extended systematically to systems of increasing complexity, where the most recent is the seven-residue mobile loop, 304-310 (Gly-His-Gly-Ala-Gly-Gly-Ser) of the enzyme porcine pancreatic alpha-amylase modeled by the AMBER force field and AMBER with the implicit solvation GB/SA (paper I, Cheluvaraja, S.; Meirovitch, H. J. Chem. Theory Comput. 2008, 4, 192). In the present paper we make a step further and extend HSMD to the same loop capped with TIP3P explicit water at 300 K. As in paper I, we are mainly interested in entropy and free energy differences between the free and bound microstates of the loop, which are obtained from two separate MD samples of these microstates. The contribution of the loop to S and F is calculated by HSMD and that of water by a particular thermodynamic integration procedure. As expected, the free microstate is more stable than the bound microstate by a total free energy difference, Ffree-Fbound=-4.8+/-1, as compared to -25.5 kcal/mol obtained with GB/SA. We find that relatively large systematic errors in the loop entropies, Sfree(loop) and Sbound(loop) are cancelled in their difference which is thus obtained efficiently and with high accuracy, i.e., with a statistical error of 0.1 kcal/mol. This cancellation, which has been observed in previous HSMD studies, is in accord with theoretical arguments given in paper I.

A simple model of backbone flexibility improves modeling of side-chain conformational variability

The considerable flexibility of side-chains in folded proteins is important for protein stability and function, and may have a role in mediating allosteric interactions. While sampling side-chain degrees of freedom has been an integral part of several successful computational protein design methods, the predictions of these approaches have not been directly compared to experimental measurements of side-chain motional amplitudes. In addition, protein design methods frequently keep the backbone fixed, an approximation that may substantially limit the ability to accurately model side-chain flexibility. Here, we describe a Monte Carlo approach to modeling side-chain conformational variability and validate our method against a large dataset of methyl relaxation order parameters derived from nuclear magnetic resonance (NMR) experiments (17 proteins and a total of 530 data points). We also evaluate a model of backbone flexibility based on Backrub motions, a type of conformational change frequently observed in ultra-high-resolution X-ray structures that accounts for correlated side-chain backbone movements. The fixed-backbone model performs reasonably well with an overall rmsd between computed and predicted side-chain order parameters of 0.26. Notably, including backbone flexibility leads to significant improvements in modeling side-chain order parameters for ten of the 17 proteins in the set. Greater accuracy of the flexible backbone model results from both increases and decreases in side-chain flexibility relative to the fixed-backbone model. This simple flexible-backbone model should be useful for a variety of protein design applications, including improved modeling of protein-protein interactions, design of proteins with desired flexibility or rigidity, and prediction of correlated motions within proteins.

Protein-ligand NOE matching: a high-throughput method for binding pose evaluation that does not require protein NMR resonance assignments

Given the three-dimensional (3D) structure of a protein, the binding pose of a ligand can be determined using distance restraints derived from assigned intra-ligand and protein-ligand nuclear Overhauser effects (NOEs). A primary limitation of this approach is the need for resonance assignments of the ligand-bound protein. We have developed an approach that utilizes data from 3D 13C-edited, 13C/15N-filtered HSQC-NOESY spectra for evaluating ligand binding poses without requiring protein NMR resonance assignments. Only the 1H NMR assignments of the bound ligand are essential. Trial ligand binding poses are generated by any suitable method (e.g., computational docking). For each trial binding pose, the 3D 13C-edited, 13C/15N-filtered HSQC-NOESY spectrum is predicted, and the predicted and observed patterns of protein-ligand NOEs are matched and scored using a fast, deterministic bipartite graph matching algorithm. The best scoring (lowest "cost") poses are identified. Our method can incorporate any explicit restraints or protein assignment data that are available, and many extensions of the basic procedure are feasible. Only a single sample is required, and the method can be applied to both slowly and rapidly exchanging ligands. The method was applied to three test cases: one complex involving muscle fatty acid-binding protein (mFABP) and two complexes involving the leukocyte function-associated antigen 1 (LFA-1) I-domain. Without using experimental protein NMR assignments, the method identified the known binding poses with good accuracy. The addition of experimental protein NMR assignments improves the results. Our "NOE matching" approach is expected to be widely applicable; i.e., it does not appear to depend on a fortuitous distribution of binding pocket residues.

Conformational plasticity of the lipid transfer protein SCP2

The nonspecific lipid transfer protein sterol carrier protein 2 (SCP2) is involved in organellar fatty acid metabolism. A hydrophobic cavity in the structure of SCP2 accommodates a wide variety of apolar ligands such as cholesterol derivatives or fatty acyl-coenzyme A (CoA) conjugates. The properties of this nonspecific lipid binding pocket are explored using NMR chemical shift perturbations, paramagnetic relaxation enhancement, amide hydrogen exchange, and 15N relaxation measurements. A common binding cavity shared by different physiological ligands is identified. NMR relaxation measurements reveal that residues in the three C-terminal alpha-helices within the lipid binding region exhibit mobility at fast (picosecond to nanosecond) and slow (microsecond to millisecond) time scales. Ligand binding is associated with a considerable loss of peptide backbone mobility. The observed conformational dynamics in SCP2 may play a role for the access of hydrophobic ligands to an occluded binding pocket. The C-terminal peroxisomal targeting signal of SCP2 is specifically recognized by the Pex5p receptor protein, which conducts cargo proteins toward the peroxisomal organelle. Neither the C-terminal targeting signal nor the N-terminal precursor sequence interferes with lipid binding by SCP2. The alpha-helices involved in lipid binding also mediate a secondary interaction interface with the Pex5p receptor. Silencing of conformational dynamics of the peptide backbone in these helices upon either lipid or Pex5p binding might communicate the loading state of the cargo protein to the targeting receptor.

Recent developments in methodologies for calculating the entropy and free energy of biological systems by computer simulation

The Helmholtz free energy, F, plays an important role in proteins because of their rugged potential energy surface, which is 'decorated' with a tremendous number of local wells (denoted microstates, m). F governs protein folding, whereas differences DeltaF(mn) determine the relative populations of microstates that are visited by a flexible cyclic peptide or a flexible protein segment (e.g. a surface loop). Recently developed methodologies for calculating DeltaF(mn) (and entropy differences, DeltaS(mn)) mainly use thermodynamic integration and calculation of the absolute F; interesting new approaches in these categories are the adaptive integration method and the hypothetical scanning molecular dynamics method, respectively.

Calculation of the entropy and free energy of peptides by molecular dynamics simulations using the hypothetical scanning molecular dynamics method

Hypothetical scanning (HS) is a method for calculating the absolute entropy S and free energy F from a sample generated by any simulation technique. With this approach each sample configuration is reconstructed with the help of transition probabilities (TPs) and their product leads to the configuration's probability, hence to the entropy. Recently a new way for calculating the TPs by Monte Carlo (MC) simulations has been suggested, where all system interactions are taken into account. Therefore, this method--called HSMC--is in principle exact where the only approximation is due to insufficient sampling. HSMC has been applied very successfully to liquid argon, TIP3P water, self-avoiding walks on a lattice, and peptides. Because molecular dynamics (MD) is considered to be significantly more efficient than MC for a compact polymer chain, in this paper HSMC is extended to MD simulations as applied to peptides. Like before, we study decaglycine in vacuum but for the first time also a peptide with side chains, (Val)(2)(Gly)(6)(Val)(2). The transition from MC to MD requires implementing essential changes in the reconstruction process of HSMD. Results are calculated for three microstates, helix, extended, and hairpin. HSMD leads to very stable differences in entropy TDeltaS between these microstates with small errors of 0.1-0.2 kcal/mol (T=100 K) for a wide range of calculation parameters with extremely high efficiency. Various aspects of HSMD and plans for future work are discussed.

NMR Structure of a Potent Small Molecule Inhibitor Bound to Human Keratinocyte Fatty Acid-Binding Protein

The NMR structure is presented for compound 1 (BMS-480404) (Ki = 33 (+/-2) nM) bound to keratinocyte fatty acid-binding protein. This article describes interactions between a high affinity drug-like compound and a member of the fatty acid-binding protein family. A benzyl group ortho to the mandelic acid in 1 occupies an area of the protein that fatty acids do not normally contact. Similar to that in the kFABP-palmitic acid structure, the acid moiety in 1 is proximal to R129 and Y131. Computational modeling indicates that the acid moiety in 1 interacts indirectly via a modeled water molecule to R109.

Optimization of the GB/SA solvation model for predicting the structure of surface loops in proteins

Implicit solvation models are commonly optimized with respect to experimental data or Poisson-Boltzmann (PB) results obtained for small molecules, where the force field is sometimes not considered. In previous studies, we have developed an optimization procedure for cyclic peptides and surface loops in proteins based on the entire system studied and the specific force field used. Thus, the loop has been modeled by the simplified solvation function E(tot) = E(FF) (epsilon = 2r) + Sigma(i) sigma(i)A(i), where E(FF) (epsilon = nr) is the AMBER force field energy with a distance-dependent dielectric function, epsilon = nr, A(i) is the solvent accessible surface area of atom i, and sigma(i) is its atomic solvation parameter. During the optimization process, the loop is free to move while the protein template is held fixed in its X-ray structure. To improve on the results of this model, in the present work we apply our optimization procedure to the physically more rigorous solvation model, the generalized Born with surface area (GB/SA) (together with the all-atom AMBER force field) as suggested by Still and co-workers (J. Phys. Chem. A 1997, 101, 3005). The six parameters of the GB/SA model, namely, P(1)-P(5) and the surface area parameter, sigma (programmed in the TINKER package) are reoptimized for a "training" group of nine loops, and a best-fit set is defined from the individual sets of optimized parameters. The best-fit set and Still's original set of parameters (where Lys, Arg, His, Glu, and Asp are charged or neutralized) were applied to the training group as well as to a "test" group of seven loops, and the energy gaps and the corresponding RMSD values were calculated. These GB/SA results based on the three sets of parameters have been found to be comparable; surprisingly, however, they are somewhat inferior (e.g, of larger energy gaps) to those obtained previously from the simplified model described above. We discuss recent results for loops obtained by other solvation models and potential directions for future studies.

Minimalist explicit solvation models for surface loops in proteins

We have performed molecular dynamics simulations of protein surface loops solvated by explicit water, where a prime focus of the study is the small numbers (e.g., ~100) of explicit water molecules employed. The models include only part of the protein (typically 500 - 1000 atoms), and the water molecules are restricted to a region surrounding the loop. In this study, the number of water molecules (N(w)) is systematically varied, and convergence with large N(w) is monitored to reveal N(w)(min), the minimum number required for the loop to exhibit realistic (fully hydrated) behavior. We have also studied protein surface coverage, as well as diffusion and residence times for water molecules as a function of N(w). A number of other modeling parameters are also tested. These include the number of environmental protein atoms explicitly considered in the model, as well as two ways to constrain the water molecules to the vicinity of the loop (where we find one of these methods to perform better when N(w) is small). The results (for RMSD and its fluctuations for four loops) are further compared to much larger, fully solvated systems (using ~10,000 water molecules under periodic boundary conditions and Ewald electrostatics), and to results for the GBSA implicit solvation model. We find that the loop backbone can stabilize with a surprisingly small number of water molecules (as low as 5 molecules per amino acid residue). The side chains of the loop require somewhat larger N(w), where the atomic fluctuations become too small if N(w) is further reduced. Thus, in general, we find adequate hydration to occur at roughly 12 water molecules per residue. This is an important result, because at this hydration level, computational times are comparable to those required for GBSA. Therefore these "minimalist explicit models" can provide a viable and potentially more accurate alternative. The importance of protein loop modeling is discussed in the context of these, and other, loop models, along with other challenges including the relevance of appropriate free energy simulation methodology for assessment of conformational stability.

A new spin probe of protein dynamics: nitrogen relaxation in 15N-2H amide groups

(15)N spin relaxation data have provided a wealth of information on protein dynamics in solution. Standard R(1), R(1)(rho), and NOE experiments aimed at (15)N[(1)H] amide moieties are complemented in this work by HA(CACO)N-type experiments allowing the measurement of nitrogen R(1) and R(1)(rho) rates at deuterated (15)N[(2)D] sites. Difference rates obtained using this approach, R(1)((15)N[(1)H]) - R(1)((15)N[(2)D]) and R(2)((15)N[(1)H]) - R(2)((15)N[(2)D]), depend exclusively on dipolar interactions and are insensitive to (15)N CSA and R(ex) relaxation mechanisms. The methodology has been tested on a sample of peptostreptococcal protein L (63 residues) prepared in 50% H(2)O-50% D(2)O solvent. The results from the new and conventional experiments are found to be consistent, with respect to both local backbone dynamics and overall protein tumbling. Combining several data sets permits evaluation of the spectral density J(omega(D) + omega(N)) for each amide site. This spectral density samples a uniquely low frequency (26 MHz at a 500 MHz field) and, therefore, is expected to be highly useful for characterizing nanosecond time scale local motions. The spectral density mapping demonstrates that, in the case of protein L, J(omega(D) + omega(N)) values are compatible with the Lipari-Szabo interpretation of backbone dynamics based on the conventional (15)N relaxation data.

Implementation and testing of stable, fast implicit solvation in molecular dynamics using the smooth-permittivity finite difference Poisson-Boltzmann method

A fast stable finite difference Poisson-Boltzmann (FDPB) model for implicit solvation in molecular dynamics simulations was developed using the smooth permittivity FDPB method implemented in the OpenEye ZAP libraries. This was interfaced with two widely used molecular dynamics packages, AMBER and CHARMM. Using the CHARMM-ZAP software combination, the implicit solvent model was tested on eight proteins differing in size, structure, and cofactors: calmodulin, horseradish peroxidase (with and without substrate analogue bound), lipid carrier protein, flavodoxin, ubiquitin, cytochrome c, and a de novo designed 3-helix bundle. The stability and accuracy of the implicit solvent simulations was assessed by examining root-mean-squared deviations from crystal structure. This measure was compared with that of a standard explicit water solvent model. In addition we compared experimental and calculated NMR order parameters to obtain a residue level assessment of the accuracy of MD-ZAP for simulating dynamic quantities. Overall, the agreement of the implicit solvent model with experiment was as good as that of explicit water simulations. The implicit solvent method was up to eight times faster than the explicit water simulations, and approximately four times slower than a vacuum simulation (i.e., with no solvent treatment).

The Origin of Protein Sidechain Order Parameter Distributions

Previous work by Wand et al. (Nature 2001, 411, 501-504) showed that the NMR order parameters characterizing the amplitude of motion of protein side chains seemed to form a multimodal distribution. At the time, no detailed explanation of this at the molecular level was offered, yet three "classes" of motion were inferred. We have analyzed a larger published data set and found that, although the distribution is multimodal, the evidence for three classes is weak. More significantly, we have been able to provide a simple physical explanation for the distributions based on the results of molecular dynamics simulations. This result will aid in the interpretation of data from NMR dynamics experiments.

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.

Water dynamics in the large cavity of three lipid-binding proteins monitored by (17)O magnetic relaxation dispersion

Intracellular lipid-binding proteins contain a large binding cavity filled with water molecules. The role played by these water molecules in ligand binding is not well understood, but their energetic and dynamic properties must be important for protein function. Here, we use the magnetic relaxation dispersion (MRD) of the water 17O resonance to investigate the water molecules in the binding cavity of three different lipid-binding proteins: heart fatty acid-binding protein (H-FABP), ileal lipid-binding protein (I-LBP) and intestinal fatty acid-binding protein (I-FABP). Whereas about half of the crystallographically visible water molecules appear to be expelled by the ligand, we find that ligand binding actually increases the number of water molecules within the cavity. At 300 K, the water molecules in the cavity exchange positions on a time-scale of about 1ns and exchange with external water on longer time-scales (0.01-1 micros). Exchange of water molecules among hydration sites within the cavity should be strongly coupled to ligand motion. Whereas a recent MD simulation indicates that the structure of the cavity water resembles a bulk water droplet, the present MRD results show that its dynamics is more than two orders of magnitude slower than in the bulk. These findings may have significant implications for the strength, specificity and kinetics of lipid binding.

Two homologous rat cellular retinol-binding proteins differ in local conformational flexibility

Cellular retinol-binding protein I (CRBP I) and cellular retinol-binding protein II (CRBP II) are closely homologous proteins that play distinct roles in the maintenance of vitamin A homeostasis. The solution structure and dynamics of CRBP I and CRBP II were compared by multidimensional NMR techniques. These studies indicated that differences in the mean backbone structures of CRBP I and CRBP II were localized primarily to the alphaII helix. Intraligand NOE cross-peaks were detected for the hydroxyl proton in the NOESY spectrum of CRBP I-bound retinol, but not for CRBP II-bound retinol, indicating that the conformational dynamics of retinol binding are different for these two proteins. As determined by Lipari-Szabo formalism, both the apo and holo forms of CRBP I and CRBP II are conformationally rigid on the pico- to nanosecond timescale. transverse relaxation optimized spectroscopy-Carr-Purcell-Meiboom-Gill -based 15N relaxation dispersion experiments at both 500 MHz and 600 MHz magnetic fields revealed that 84 and 62 residues for apo-CRBP I and II, respectively, showed detectable conformational exchange on a micro- to millisecond timescale, in contrast to three and seven residues for holo-CRBP I and II, respectively. Thus binding of retinol markedly reduced conformational flexibility in both CRBP I and CRBP II on the micro- to millisecond timescale. The 15N relaxation dispersion curves of apo-CRBP I and II were fit to a two-state conformational exchange model by a global iterative fitting process and by an individual (residue) fitting process. In the process of carrying out the global fit, more than half of the residue sites were eliminated. The individual chemical exchange rates k(ex), and chemical shift differences, Deltadelta, were increased in the putative portal region (alphaII helix and betaC-betaD turn) of apo-CRBP II compared to apo-CRBP I. These differences in conformational flexibility likely contribute to differences in how CRBP I and CRBP II interact with ligands, membranes and retinoid metabolizing enzymes.

Automated nano-electrospray mass spectrometry for protein-ligand screening by noncovalent interaction applied to human H-FABP and A-FABP

A method for ligand screening by automated nano-electrospray ionization mass spectrometry (nano-ESI/MS) is described. The core of the system consisted of a chip-based platform for automated sample delivery from a 96-well plate and subsequent analysis based on noncovalent interactions. Human fatty acid binding protein, H-FABP (heart) and A-FABP (adipose), with small potential ligands was analyzed. The technique has been compared with a previously reported method based on nuclear magnetic resonance (NMR), and excellent correlation with the found hits was obtained. In the current MS screening method, the cycle time per sample was 1.1 min, which is approximately 50 times faster than NMR for single compounds and approximately 5 times faster for compound mixtures. High reproducibility was achieved, and the protein consumption was in the range of 88 to 100 picomoles per sample. Futhermore, a novel protocol for preparation of A-FABP without the natural ligand is presented. The described screening approach is suitable for ligand screening very early in the drug discovery process before conventional high-throughput screens (HTS) are developed and/or used as a secondary screening for ligands identified by HTS.

Structure-based screening as applied to human FABP4: a highly efficient alternative to HTS for hit generation

The time-limiting step in HTS often is the development of an appropriate assay. In addition, hits from HTS fairly often turn out to be false positives and generally display unfavorable properties for further development. Here we describe an alternative process for hit generation, applied to the human adipocyte fatty acid binding protein FABP4. A small molecular ligand for FABP4 that blocks the binding of endogenous ligands may be developed into a drug for the treatment of type-2 diabetes. Using NMR spectroscopy, we screened FABP4 for low-affinity binders in a diversity library consisting of small soluble scaffolds, which yielded 52 initial hits in total. The potencies of these hits were ranked, and crystal structures of FABP4 complexes for two of the hits were obtained. The structural data were subsequently used to direct similarity searches for available analogues, as well as chemical synthesis of 12 novel analogues. In this way, a series of three selective FABP4 ligands with attractive pharmacochemical profiles and potencies of 10 microM or better was obtained.

Solution structure and backbone dynamics of human epidermal-type fatty acid-binding protein (E-FABP)

Human epidermal-type fatty acid-binding protein (E-FABP) belongs to a family of intracellular 14-15 kDa lipid-binding proteins, whose functions have been associated with fatty acid signalling, cell growth, regulation and differentiation. As a contribution to understanding the structure-function relationship, we report in the present study features of its solution structure and backbone dynamics determined by NMR spectroscopy. Applying multi-dimensional high-resolution NMR techniques on unlabelled and 15N-enriched recombinant human E-FABP, the 1H and 15N resonance assignments were completed. On the basis of 2008 distance restraints, the three-dimensional solution structure of human E-FABP was subsequently obtained (backbone atom root-mean-square deviation of 0.92+/-0.11 A; where 1 A=0.1 nm), consisting mainly of 10 anti-parallel beta-strands that form a beta-barrel structure. 15N relaxation experiments (T1, T2 and heteronuclear nuclear Overhauser effects) at 500, 600 and 800 MHz provided information on the internal dynamics of the protein backbone. Nearly all non-terminal backbone amide groups showed order parameters S(2)>0.8, with an average value of 0.88+/-0.04, suggesting a uniformly low backbone mobility in the nanosecond-to-picosecond time range. Moreover, hydrogen/deuterium exchange experiments indicated a direct correlation between the stability of the hydrogen-bonding network in the beta-sheet structure and the conformational exchange in the millisecond-to-microsecond time range. The features of E-FABP backbone dynamics elaborated in the present study differ markedly from those of the phylogenetically closely related heart-type FABP and the more distantly related ileal lipid-binding protein, implying a strong interdependence with the overall protein stability and possibly also with the ligand-binding affinity for members of the lipid-binding protein family.

Deuterium spin probes of side-chain dynamics in proteins. 1. Measurement of five relaxation rates per deuteron in (13)C-labeled and fractionally (2)H-enriched proteins in solution

New pulse sequences are presented for the measurement of the relaxation of deuterium double quantum, quadrupolar order, and transverse antiphase magnetization in (13)CH(2)D methyl groups of (15)N-, (13)C-labeled, fractionally deuterated proteins. Together with previously developed experiments for measuring deuterium longitudinal and transverse decay rates [Muhandiram, D. R.; Yamazaki, T.; Sykes, B. D.; Kay, L. E. J. Am. Chem. Soc. 1995, 117, 11536], these schemes allow measurement of the five unique decay constants of a single deuteron, providing an unprecedented opportunity to investigate side-chain dynamics in proteins. All five deuterium relaxation rates have been measured for deuterons in the methyl groups of the B1 immunoglobulin binding domain of peptostreptococcal protein L and the N-terminal SH3 domain from the protein drk. Since values of the spectral density function at only three different frequencies contribute to the five relaxation rates, the self-consistency of the relaxation data is readily established. Very good agreement is obtained between calculated parameters describing the amplitudes and time scales of motion when different subsets of the relaxation data are employed.

Deuterium spin probes of side-chain dynamics in proteins. 2. Spectral density mapping and identification of nanosecond time-scale side-chain motions

In the previous paper in this issue we have demonstrated that it is possible to measure the five different relaxation rates of a deuteron in (13)CH(2)D methyl groups of (13)C-labeled, fractionally deuterated proteins. The extensive set of data acquired in these experiments provides an opportunity to investigate side-chain dynamics in proteins at a level of detail that heretofore was not possible. The data, acquired on the B1 domain of peptostreptococcal protein L, include 16 (9) relaxation measurements at 4 (2) different magnetic field strengths, 25 degrees C (5 degrees C). These data are shown to be self-consistent and are analyzed using a spectral density mapping procedure which allows extraction of values of the spectral density function at a number of frequencies with no assumptions about the underlying dynamics. Dynamics data from 31 of 35 methyls in the protein for which data could be obtained were well-fitted using the two-parameter Lipari-Szabo model (Lipari, G.; Szabo, A. J. Am. Chem. Soc. 1982, 104, 4546). The data from the remaining 4 methyls can be fitted using a three-parameter version of the Lipari-Szabo model that takes into account, in a simple manner, additional nanosecond time-scale local dynamics. This interpretation is supported by analysis of a molecular dynamics trajectory where spectral density profiles calculated for side-chain methyl sites reflect the influence of slower (nanosecond) time-scale motions involving jumps between rotameric wells. A discussion of the minimum number of relaxation measurements that are necessary to extract the full complement of dynamics information is presented along with an interpretation of the extracted dynamics parameters.

Site-selective screening by NMR spectroscopy with labeled amino acid pairs

A new method for site-selective screening by NMR is presented. The core of the new method is the dual amino acid sequence specific labeling technique. Amino acid X is labeled with (13)C and amino acid Y is labeled with (15)N. Provided only one XY pair occurs in the amino acid sequence, only one signal in the 1D carbonyl (13)C spectrum will display a splitting due to the (1)J(C'N) coupling. Using this labeling strategy it is possible to screen selectively for binding to a selected epitope without the need for sequence specific assignments. An HNCO spectrum (1D or 2D) can be used either directly as a screening experiment or indirectly to identify what signals to monitor in a 2D (1)H-(15)N correlation spectrum. Chemical shift perturbations upon addition of a potential ligand are easily detected even for large proteins due to the reduced spectral complexity resulting from the use of a selectively labeled sample. The new technique is demonstrated on the human adipocyte fatty acid binding protein FABP-4. Due to the reduced spectral complexity, the method should be applicable to larger proteins than are conventional methods.

Main chain and side chain dynamics of oxidized flavodoxin from Cyanobacterium anabaena

Oxidized flavodoxin from Cyanobacterium anabaena PCC 7119 is used as a model system to investigate the fast internal dynamics of a flavin-bearing protein. Virtually complete backbone and side chain resonance NMR assignments of an oxidized flavodoxin point mutant (C55A) have been determined. Backbone and side chain dynamics in flavodoxin (C55A) were investigated using (15)N amide and deuterium methyl NMR relaxation methods. The squared generalized order parameters (S(NH)(2)) for backbone amide N-H bonds are found to be uniformly high (<S(NH)(2)> approximately 0.923 over 109 residues in regular secondary structure), indicating considerable restriction of motion in the backbone of the protein. In contrast, methyl-bearing side chains are considerably heterogeneous in their amplitude of motion, as indicated by obtained symmetry axis squared generalized order parameters (S(axis)(2)). However, in comparison to nonprosthetic group-bearing proteins studied with these NMR relaxation methods, the side chains of oxidized flavodoxin are unusually rigid.

Role of surface lysine residues of adipocyte fatty acid-binding protein in fatty acid transfer to phospholipid vesicles

The tertiary structure of murine adipocyte fatty acid-binding protein (AFABP) is a flattened 10-stranded beta-barrel capped by a helix-turn-helix segment. This helical domain is hypothesized to behave as a "lid" or portal for ligand entry into and exit from the binding cavity. Previously, we demonstrated that anthroyloxy-labeled fatty acid (AOFA) transfer from AFABP to phospholipid membranes occurs by a collisional process, in which ionic interactions between positively charged lysine residues on the protein surface and negatively charged phospholipid headgroups are involved. In the present study, the role of specific lysine residues located in the portal and other regions of AFABP was directly examined using site-directed mutagenesis. The results showed that isoleucine replacement for lysine in the portal region, including the alphaI- and alphaII-helices and the beta C-D turn, resulted in much slower 2-(9-anthroyloxy)palmitate (2AP) transfer rates to acidic membranes than those of native AFABP. An additive effect was found for mutant K22,59I, displaying the slowest rates of FA transfer. Rates of 2AP transfer from "nonportal" mutants on the beta-G and I strands were affected only moderately; however, a lysine --> isoleucine mutation in the nonportal beta-A strand decreased the 2AP transfer rate. These studies suggest that lysines in the helical cap domain are important for governing ionic interactions between AFABP and membranes. Furthermore, it appears that more than one distinct region, including the alphaI-helix, alphaII-helix, beta C-D turn, and the beta-A strand, is involved in these charge-charge interactions.

Microscopic origins of entropy, heat capacity and the glass transition in proteins

Internal motion is central to protein folding, to protein stability through the resulting residual entropy, and to protein function. Despite its importance, the precise nature of the internal motions of protein macromolecules remains a mystery. Here we report a survey of the temperature dependence of the fast dynamics of methyl-bearing side chains in a calmodulin-peptide complex using site-specific deuterium NMR relaxation methods. The amplitudes of motion had a markedly heterogeneous spectrum and segregated into three largely distinct classes. Other proteins studied at single temperatures tend to segregate similarly. Furthermore, a large variability in the degree of thermal activation of the dynamics in the calmodulin complex indicates a heterogeneous distribution of residual entropy and hence reveals the microscopic origins of heat capacity in proteins. These observations also point to an unexpected explanation for the low-temperature 'glass transition' of proteins. It is this transition that has been ascribed to the creation of protein motional modes that are responsible for biological activity.

Main chain and side chain dynamics of a heme protein: 15N and 2H NMR relaxation studies of R. capsulatus ferrocytochrome c2

A detailed characterization of the main chain and side chain dynamics in R. capsulatus ferrocytochrome c(2) derived from (2)H NMR relaxation of methyl group resonances is presented. (15)N relaxation measurements confirm earlier results indicating that R. capsulatus ferrocytochrome c(2) exhibits minor rotational anisotropy in solution. The current study is focused on the use of deuterium relaxation in side chain methyl groups, which has been shown to provide a detailed and accurate measure of internal dynamics. Results obtained indicate that the side chains of ferrocytochrome c(2) exhibit a wide range of motional amplitudes, but are more rigid than generally found in the interior of nonprosthetic group bearing globular proteins. This unusual rigidity is ascribed to the interactions of the protein with the large heme prosthetic group. This observation has significant implications for the potential of the heme-protein interface to modulate the redox properties of the protein and also points to the need for great precision in the design and engineering of heme proteins.

Fluctuations in ion pairs and their stabilities in proteins

This report investigates the effect of systemic protein conformational flexibility on the contribution of ion pairs to protein stability. Toward this goal, we use all NMR conformer ensembles in the Protein Data Bank (1) that contain at least 40 conformers, (2) whose functional form is monomeric, (3) that are nonredundant, and (4) that are large enough. We find 11 proteins adhering to these criteria. Within these proteins, we identify 22 ion pairs that are close enough to be classified as salt bridges. These are identified in the high-resolution crystal structures of the respective proteins or in the minimized average structures (if the crystal structures are unavailable) or, if both are unavailable, in the "most representative" conformer of each of the ensembles. We next calculate the electrostatic contribution of each such ion pair in each of the conformers in the ensembles. This results in a comprehensive study of 1,201 ion pairs, which allows us to look for consistent trends in their electrostatic contributions to protein stability in large sets of conformers. We find that the contributions of ion pairs vary considerably among the conformers of each protein. The vast majority of the ion pairs interconvert between being stabilizing and destabilizing to the structure at least once in the ensembles. These fluctuations reflect the variabilities in the location of the ion pairing residues and in the geometric orientation of these residues, both with respect to each other, and with respect to other charged groups in the remainder of the protein. The higher crystallographic B-factors for the respective side-chains are consistent with these fluctuations. The major conclusion from this study is that salt bridges observed in crystal structure may break, and new salt bridges may be formed. Hence, the overall stabilizing (or, destabilizing) contribution of an ion pair is conformer population dependent.

Optimization of solvation models for predicting the structure of surface loops in proteins

A novel procedure for optimizing the atomic solvation parameters (ASPs) sigma(i) developed recently for cyclic peptides is extended to surface loops in proteins. The loop is free to move, whereas the protein template is held fixed in its X-ray structure. The energy is E(tot) = E(FF)(epsilon = nr) + summation operator sigma(i)A(i), where E(FF)(epsilon = nr) is the force-field energy of the loop-loop and loop-template interactions, epsilon = nr is a distance-dependent dielectric constant, and n is an additional parameter to be optimized. A(i) is the solvent-accessible surface area of atom i. The optimal sigma(i) and n are those for which the loop structure with the global minimum of E(tot)(n, sigma(i)) becomes the experimental X-ray structure. Thus, the ASPs depend on the force field and are optimized in the protein environment, unlike commonly used ASPs such as those of Wesson and Eisenberg (Protein Sci 1992;1:227-235). The latter are based on the free energy of transfer of small molecules from the gas phase to water and have been traditionally combined with various force fields without further calibration. We found that for loops the all-atom AMBER force field performed better than OPLS and CHARMM22. Two sets of ASPs [based on AMBER (n = 2)], optimized independently for loops 64-71 and 89-97 of ribonuclease A, were similar and thus enabled the definition of a best-fit set. All these ASPs were negative (hydrophilic), including those for carbon. Very good (i.e., small) root-mean-square-deviation values from the X-ray loop structure were obtained with the three sets of ASPs, suggesting that the best-fit set would be transferable to loops in other proteins as well. The structure of loop 13-24 is relatively stretched and was insensitive to the effect of the ASPs.

Spin-system heterogeneities indicate a selected-fit mechanism in fatty acid binding to heart-type fatty acid-binding protein (H-FABP)

Recent advances in the characterization of fatty acid-binding proteins (FABPs) by NMR have enabled various research groups to investigate the function of these proteins in aqueous solution. The binding of fatty acid molecules to FABPs, which proceeds through a portal region on the protein surface, is of particular interest. In the present study we have determined the three-dimensional solution structure of human heart-type FABP by multi-dimensional heteronuclear NMR spectroscopy. Subsequently, in combination with data collected on a F57S mutant we have been able to show that different fatty acids induce distinct conformational states of the protein backbone in this portal region, depending on the chain length of the fatty acid ligand. This indicates that during the binding process the protein accommodates the ligand molecule by a "selected-fit" mechanism. In fact, this behaviour appears to be especially pronounced in the heart-type FABP, possibly due to a more rigid backbone structure compared with other FABPs, as suggested by recent NMR relaxation studies. Thus differences in the dynamic behaviours of these proteins may be the key to understanding the variations in ligand affinity and specificity within the FABP family.