Efficient analysis of macromolecular rotational diffusion from heteronuclear relaxation data

A novel program has been developed for the interpretation of ¹⁵N relaxation rates in terms of macromolecular anisotropic rotational diffusion. The program is based on a highly efficient simulated annealing/minimization algorithm, designed specifically to search the parametric space described by the isotropic, axially symmetric and fully anisotropic rotational diffusion tensor models. The high efficiency of this algorithm allows extensive noise-based Monte Carlo error analysis. Relevant statistical tests are systematically applied to provide confidence limits for the proposed tensorial models. The program is illustrated here using the example of the cytochrome c from Rhodobacter capsulatus, a four-helix bundle heme protein, for which data at three different field strengths were independently analysed and compared.     

Backbone dynamics of the oligomerization domain of p53 determined from 15N NMR relaxation measurements.

The backbone dynamics of the tetrameric p53 oligomerization domain (residues 319-360) have been investigated by two-dimensional inverse detected heteronuclear 1H-15N NMR spectroscopy at 500 and 600 MHz. 15N T1, T2, and heteronuclear NOEs were measured for 39 of 40 non-proline backbone NH vectors at both field strengths. The overall correlation time for the tetramer, calculated from the T1/T2 ratios, was found to be 14.8 ns at 35 degrees C. The correlation times and amplitudes of the internal motions were extracted from the relaxation data using the model-free formalism (Lipari G, Szabo A, 1982, J Am Chem Soc 104:4546-4559). The internal dynamics of the structural core of the p53 oligomerization domain are uniform and fairly rigid, with residues 327-354 exhibiting an average generalized order parameter (S2) of 0.88 +/- 0.08. The N- and C-termini exhibit substantial mobility and are unstructured in the solution structure of p53. Residues located at the N- and C-termini, in the beta-sheet, in the turn between the alpha-helix and beta-sheet, and at the C-terminal end of the alpha-helix display two distinct internal motions that are faster than the overall correlation time. Fast internal motions (< or = 20 ps) are within the extreme narrowing limit and are of uniform amplitude. The slower motions (0.6-2.2 ns) are outside the extreme narrowing limit and vary in amplitude.(ABSTRACT TRUNCATED AT 250 WORDS)     

Characterization of the overall and local dynamics of a protein with intermediate rotational anisotropy: Differentiating between conformational exchange and anisotropic diffusion in the B3 domain of protein G

Because the overall tumbling provides a major contribution to protein spectral densities measured in solution, the choice of a proper model for this motion is critical for accurate analysis of protein dynamics. Here we study the overall and backbone dynamics of the B3 domain of protein G using ¹⁵N relaxation measurements and show that the picture of local motions is markedly dependent on the model of overall tumbling. The main difference is in the interpretation of the elevated R ₂ values in the -helix: the isotropic model results in conformational exchange throughout the entire helix, whereas no exchange is predicted by anisotropic models that place the longitudinal axis of diffusion tensor almost parallel to the helix axis. Due to small size (fast tumbling) of the protein, the T ₁ values have low sensitivity to NH bond orientation. The diffusion tensor derived from orientation dependence of R ₂/R ₁ is anisotropic (D _par/D _perp=1.4), with a small rhombic component. In order to distinguish the correct picture of motion, we apply model-independent methods that are sensitive to conformational exchange and do not require knowledge of protein structure or assumptions about its dynamics. A comparison of the CSA/dipolar cross-correlation rate constants with ¹⁵N relaxation rates and the estimation of R _ex terms from relaxation data at 9.4 and 14.1 T indicate no conformational exchange in the helix, in support of the anisotropic models. The experimentally derived diffusion tensor is in excellent agreement with theoretical predictions from hydrodynamic calculations; a detailed comparison with various hydrodynamic models revealed optimal parameters for hydrodynamic calculations.     

Anisotropic rotational diffusion in model-free analysis for a ternary DHFR complex

Model-free analysis has been extensively used to extract information on motions in proteins over a wide range of timescales from NMR relaxation data. We present a detailed analysis of the effects of rotational anisotropy on the model-free analysis of a ternary complex for dihydrofolate reductase (DHFR). Our findings show that the small degree of anisotropy exhibited by DHFR (D_||/D=1.18) introduces erroneous motional models, mostly exchange terms, to over 50% of the NH spins analyzed when isotropic tumbling is assumed. Moreover, there is a systematic change in S², as large as 0.08 for some residues. The significant effects of anisotropic rotational diffusion on model-free motional parameters are in marked contrast to previous studies and are accentuated by lowering of the effective correlation time using isotropic tumbling methods. This is caused by the preponderance of NH vectors aligned perpendicular to the principal diffusion tensor axis and is readily detected because of the high quality of the relaxation data. A novel procedure, COPED (COmparison of Predicted and Experimental Diffusion tensors) is presented for distinguishing genuine motions from the effects of anisotropy by comparing experimental relaxation data and data predicted from hydrodynamic analyses. The procedure shows excellent agreement with the slow motions detected from the axially symmetric model-free analysis and represents an independent procedure for determining rotational diffusion and slow motions that can confirm or refute established procedures that rely on relaxation data. Our findings show that neglect of even small degrees of rotational diffusion anisotropy can introduce significant errors in model-free analysis when the data is of high quality. These errors can hinder our understanding of the role of internal motions in protein function.     

A suite of Mathematica notebooks for the analysis of protein main chain 15N NMR relaxation data

A suite of Mathematica notebooks has been designed to ease the analysis of protein main chain ¹⁵N NMR relaxation data collected at a single magnetic field strength. Individual notebooks were developed to perform the following tasks: nonlinear fitting of ¹⁵N-T ₁ and -T ₂ relaxation decays to a two parameter exponential decay, calculation of the principal components of the inertia tensor from protein structural coordinates, nonlinear optimization of the principal components and orientation of the axially symmetric rotational diffusion tensor, model-free analysis of ¹⁵N-T ₁, -T ₂, and {¹H}–¹⁵N NOE data, and reduced spectral density analysis of the relaxation data. The principle features of the notebooks include use of a minimal number of input files, integrated notebook data management, ease of use, cross-platform compatibility, automatic visualization of results and generation of high-quality graphics, and output of analyses in text format.L. Spyracopoulos is an AHFMR Medical Research Senior Scholar     

Investigation of the backbone dynamics of the IgG-binding domain of streptococcal protein G by heteronuclear two-dimensional 1H-15N nuclear magnetic resonance spectroscopy.

The backbone dynamics of the immunoglobulin-binding domain (B1) of streptococcal protein G, uniformly labeled with 15N, have been investigated by two-dimensional inverse detected heteronuclear 1H-15N NMR spectroscopy at 500 and 600 MHz. 15N T1, T2, and nuclear Overhauser enhancement data were obtained for all 55 backbone NH vectors of the B1 domain at both field strengths. The overall correlation time obtained from an analysis of the T1/T2 ratios was 3.3 ns at 26 degrees C. Overall, the B1 domain is a relatively rigid protein, consistent with the fact that over 95% of the residues participate in secondary structure, comprising a four-stranded sheet arranged in a -1, +3x, -1 topology, on top of which lies a single helix. Residues in the turns and loops connecting the elements of secondary structure tend to exhibit a higher degree of mobility on the picosecond time scale, as manifested by lower values of the overall order parameter. A number of residues at the ends of the secondary structure elements display two distinct internal motions that are faster than the overall rotational correlation time: one is fast (< 20 ps) and lies in the extreme narrowing limit, whereas the other is one to two orders of magnitude slower (1-3 ns) and lies outside the extreme narrowing limit. The slower motion can be explained by large-amplitude (20-40 degrees) jumps in the N-H vectors between states with well-defined orientations that are stabilized by hydrogen bonds.(ABSTRACT TRUNCATED AT 250 WORDS)     

Interpretation of 15N NMR relaxation data of globular proteins using hydrodynamic calculations with HYDRONMR

HYDRONMR is an implementation of state of the art hydrodynamic modeling to calculate the spectral density functions for NH or C-H vectors in a rigid protein structure starting from an atomic level representation. Thus HYDRONMR can be used to predict NMR relaxation times from a rigid model and to compare them with the experimental results. HYDRONMR contains a single adjustable parameter, the atomic element radius. A protocol to determine the value that gives the best agreement between calculated and experimental T₁/T₂values is described. For most proteins, the value of the atomic element radius ranges between 2.8 Å and 3.8 Å with a distribution centered at 3.3 Å. Deviations from the usual range towards larger values are associated to aggregation in several proteins. Deviations to lower values may be related to large-scale motions or inappropriate model structures.If the average structure is correct, deviations between experimental T₁/T₂values and those calculated with HYDRONMR can be used to distinguish residues affected by anisotropic motion from those that are involved in chemical exchange.     

Anisotropic rotational diffusion of perdeuterated HIV protease from 15N NMR relaxation measurements at two magnetic fields

Summary ¹⁵N NMR relaxation times in perdeuterated HIV-1 protease, complexed with the sub-nanomolar inhibitor DMP323, have been measured at 600 and 360 MHz ¹H frequency. The relative magnitudes of the principal components of the inertia tensor, calculated from the X-ray coordinates of the protein-drug complex, are 1.0:0.85:0.44. The relation between the T₁/T₂ ratios observed for the individual backbone amides and their N-H orientation within the 3D structure of the protease dimer yields a rotational diffusion tensor oriented nearly collinear to the inertia tensor. The relative magnitudes of its principal components (1.00:1.11:1.42) are also in good agreement with hydrodynamic modeling results. The orientation and magnitude of the diffusion tensors derived from relaxation data obtained at 360 and 600 MHz are nearly identical. The anisotropic nature of the rotational diffusion has little influence on the order parameters derived from the ¹⁵N T₁ and T₂ relaxation times; however, if anisotropy is ignored, this can result in erroneous identification of either exchange broadening or internal motions on a nanosecond time scale. The average ratio of the T₁ values measured at 360 and 600 MHz is 0.50±0.015, which is slightly larger than the value of 0.466 expected for an isotropic rigid rotor with _c = 10.7 ns. The average ratio of the T₂ values measured at 360 and 600 MHz is 1.14±0.04, which is also slightly larger than the expected ratio of 1.11. This magnetic field dependence of the T₁ and T₂ relaxation times suggests that the spectral density contribution from fast internal motions is not negligible, and that the chemical shift anisotropy of peptide backbone amides, on average, is larger than the 160 ppm value commonly used in ¹⁵N relaxation studies of proteins.     

Structural and dynamic characterization of the urea denatured state of the immunoglobulin binding domain of streptococcal protein G by multidimensional heteronuclear NMR spectroscopy.

The structure and dynamics of the urea-denatured B1 immunoglobulin binding domain of streptococcal protein G (GB1) has been investigated by multidimensional heteronuclear NMR spectroscopy. Complete 1H, 15N, and 13C assignments are obtained by means of sequential through-bond correlations. The nuclear Overhauser enhancement, chemical shift, and 3JHN alpha coupling constant data provide no evidence for the existence of any significant population of residual native or nonnative ordered structure. 15N relaxation measurements at 500 and 600 MHz, however, provide evidence for conformationally restricted motions in three regions of the polypeptide that correspond to the second beta-hairpin, the N-terminus of the alpha-helix, and the middle of the alpha-helix in the native protein. The time scale of these motions is longer than the apparent overall correlation time (approximately 3 ns) and could range from about 6 ns in the case of one model to between 4 microseconds and 2 ms in another; it is not possible to distinguish between these two cases with certainty because the dynamics are highly complex and hence the analysis of the time scale of this slower motion is highly model dependent. It is suggested that these three regions may correspond to nucleation sites for the folding of the GB1 domain. With the exception of the N- and C-termini, where end effects predominate, the amplitude of the subnanosecond motions, on the other hand, are fairly uniform and model independent, with an overall order parameter S2 ranging from 0.4 to 0.5.     

Internal motional amplitudes and correlated bond rotations in an alpha-helical peptide derived from 13C and 15N NMR relaxation

Peptide GFSKAELAKARAAKRGGY folds in an alpha-helical conformation that is stabilized by formation of a hydrophobic staple motif and an N-terminal capping box (Munoz V. Blanco FJ, Serrano L, 1995, Struct Biol 2:380-385). To investigate backbone and side-chain internal motions within the helix and hydrophobic staple, residues F2, A5, L7, A8, and A10 were selectively 13C- and 15N-enriched and NMR relaxation experiments were performed in water and in water/trifluoroethanol (TFE) solution at four Larmor frequencies (62.5, 125, 150, and 200 MHz for 13C). Relaxation data were analyzed using the model free approach and an anisotropic diffusion model. In water, angular variances of motional vectors range from 10 to 20 degrees and backbone phi,psi bond rotations for helix residues A5, L7, A8, and A10 are correlated indicating the presence of Calpha-H, Calpha-Cbeta, and N-H rocking-type motions along the helix dipole axis. L7 side-chain CbetaH2 and CgammaH motions are also correlated and as motionally restricted as backbone CalphaH, suggesting considerable steric hindrance with neighboring groups. In TFE which stabilizes the fold, internal motional amplitudes are attenuated and rotational correlations are increased. For the side chain of hydrophobic staple residue F2, wobbling-in-a-cone type motions dominate in water, whereas in TFE, the Cbeta-Cgamma bond and phenyl ring fluctuate more simply about the Calpha-Cbeta bond. These data support the Daragan-Mayo model of correlated bond rotations (Daragan VA, Mayo KH, 1996, J Phys Chem 100:8378-8388) and contribute to a general understanding of internal motions in peptides and proteins.     

Backbone dynamics of free barnase and its complex with barstar determined by 15N NMR relaxation study

Backbone dynamics of uniformly ¹⁵N-labeled free barnase and its complex with unlabelled barstar have been studied at 40°C, pH 6.6, using ¹⁵N relaxation data obtained from proton-detected 2D {¹H}-¹⁵N NMR spectroscopy. ¹⁵N spin-lattice relaxation rate constants (R₁), spin-spin relaxation rate constants (R₂), and steady-state heteronuclear {¹H}-¹⁵N NOEs have been measured at a magnetic field strength of 14.1 Tesla for 91 residues of free barnase and for 90 residues out of a total of 106 in the complex (excluding three prolines and the N-terminal residue) backbone amide ¹⁵N sites of barnase. The primary relaxation data for both the cases have been analyzed in the framework of the model-free formalism using both isotropic and axially symmetric models of the rotational diffusion tensor. As per the latter, the overall rotational correlation times (_m) are 5.0 and 9.5 ns for the free and complexed barnase, respectively. The average order parameter is found to be 0.80 for free barnase and 0.86 for the complex. However, the changes are not uniform along the backbone and for about 5 residues near the binding interface there is actually a significant decrease in the order parameters on complex formation. These residues are not involved in the actual binding. For the residues where the order parameter increases, the magnitudes vary significantly. It is observed that the complex has much less internal mobility, compared to free barnase. From the changes in the order parameters, the entropic contribution of NH bond vector motion to the free energy of complex formation has been calculated. It is apparent that these motions cause significant unfavorable contributions and therefore must be compensated by many other favorable contributions to effect tight complex formation. The observed variations in the motion and their different locations with regard to the binding interface may have important implications for remote effects and regulation of the enzyme action.     

Determination of protein rotational correlation time from NMR relaxation data at various solvent viscosities

An accurate determination of the overall rotation of a protein plays a crucial role in the investigation of its internal motions by NMR. In the present work, an innovative approach to the determination of the protein rotational correlation time _R from the heteronuclear relaxation data is proposed. The approach is based on a joint fit of relaxation data acquired at several viscosities of a protein solution. The method has been tested on computer simulated relaxation data as compared to the traditional _R determination method from T₁/T₂ ratio. The approach has been applied to ribonuclease barnase from Bacillus amyloliquefaciens dissolved in an aqueous solution and deuterated glycerol as a viscous component. The resulting rotational correlation time of 5.56 ± 0.01 ns and other rotational diffusion tensor parameters are in good agreement with those determined from T₁/T₂ ratio.     

Mobility of NH bonds in DNA-binding protein HU of shape Bacillus stearothermophilus from reduced spectral density mapping analysis at multiple NMR fields

The dynamics of the backbone NH bonds of protein HU from Bacillus stearothermophilus (HUBst) have been characterized using measurements of cross-relaxation, longitudinal and transverse relaxation rates at 11.7, 14.1 and 17.6 T. Linear regression of the values with the squared Larmor frequency _N ² has revealed global exchange processes, which contributed on the order of 0.5–5.0 s^-1to the transverse relaxation rate. Subsequently, the experimental values were corrected for these exchange contributions. A reduced spectral density mapping procedure has been employed with the experimental relaxation rates and seven values of the spectral density function J() have been extracted. These spectral densities have been fitted within the framework of the model-free approach. The densities agree well with an axially symmetric rotational diffusion tensor with a diffusion anisotropy D_/D_ of 1.15, indicating that the flexible arms of HUBst do not significantly contribute to the rotational diffusion. The overall correlation time is 8.9 ± 0.6 ns/rad. The fast internal motions of most of the NH bonds in the core display order parameters ranging between 0.74 and 0.83 and internal correlation times between 1 and 20 ps. For the residues in the DNA-binding -arms, an extended version of the model function has been used. The slow internal motions show correlation times of 1–2 ns. The concomitant order parameters (0.3–0.6) are lower than those observed on the fast time scale, indicating that the flexibility of the -arms is mainly determined by the slower internal motions. A substantial decrease of the generalized order parameters in the -arms starting at residues Arg⁵⁵ and Ser⁷⁴, opposite on both strands of the -ribbon arms, has been explained as a hinge motion. A comparison of the order parameters for free and DNA-bound protein has demonstrated that the slow hinge motions largely disappear when HU binds DNA.     

The relative orientation of the fibronectin 6F11F2 module pair: A 15N NMR relaxation study

The structure of a pair of modules (⁶F1¹F2), that forms part of the collagen-binding region of fibronectin, is refined using heteronuclear relaxation data. A structure of the pair was previously derived from ¹H-¹H NOE and ³ J _HHN data [Bocquier et al. (1999) Structure, 7, 1451–1460] and a weak module–module interface, comprising Leu19 and Leu28, in ⁶F1, and Tyr68 in ²F1, was identified. In this study, the definition of the average relative orientation of the two modules is improved using the dependence of ¹⁵N relaxation on rotational diffusion anisotropy. This structure refinement is based on the selection of a subset of structures from sets calculated with NOE and ³ J _HHN data alone, using the quality of the fits to the relaxation data as the selection criterion. This simple approach is compared to a refinement strategy where ¹⁵N relaxation data are included in the force field as additional restraints [Tjandra et al. (1997) Nat. Struct. Biol., 4, 443–449].     

Measurement of 15N relaxation in the detergent-solubilized tetrameric KcsA potassium channel

A set of TROSY-HNCO (tHNCO)-based 3D experiments is presented for measuring ¹⁵N relaxation parameters in large, membrane-associated proteins, characterized by slow tumbling times and significant spectral overlap. Measurement of backbone ¹⁵N R ₁, R _1ρ, ¹⁵N–{¹H} NOE, and ¹⁵N CSA/dipolar cross correlation is demonstrated and applied to study the dynamic behavior of the homotetrameric KcsA potassium channel in SDS micelles under conditions where this channel is in the closed state. The micelle-encapsulated transmembrane domain, KcsA^TM, exhibits a high degree of order, tumbling as an oblate ellipsoid with a global rotational correlation time, τ_c = 38 ± 2.5 ns, at 50 °C and a diffusion anisotropy, , corresponding to an aspect ratio a/b ≥ 1.4. The N- and C-terminal intracellular segments of KcsA exhibit considerable internal dynamics (S ² values in the 0.2–0.45 range), but are distinctly more ordered than what has been observed for unstructured random coils. Relaxation behavior in these domains confirms the position of the C-terminal helix, and indicates that in SDS micelles, this amphiphilic helix does not associate into a stable homotetrameric helical bundle. The relaxation data indicate the absence of elevated backbone dynamics on the ps–ns time scale for the 5-residue selectivity filter, which selects K⁺ ions to enter the channel. Electronic Supplementary Material Supplementary material is available to authorised users in the online version of this article at . An erratum to this article can be found at     

15N relaxation study of the cold shock protein CspB at various solvent viscosities

For a detailed NMR study of the dynamics of the cold shock protein CspB from Bacillus subtilis, we determined ¹⁵N transverse and longitudinal relaxation rates and heteronuclear nuclear Overhauser effects at different solvent viscosities. Up to a relative viscosity of 2, which is equivalent to 27% ethylene glycol (EG), the overall correlation time follows the linear Stokes-Einstein equation. At a relative viscosity of 6 (70% EG) the correlation time deviates from linearity by 30%, indicating that CspB tumbles at a higher rate as expected from the solvent viscosity probably due to a preferential binding of water molecules at the protein surface. The corresponding hydrodynamic radii, determined by NMR diffusion experiments, show no variation with viscosity. The amplitudes of intramolecular motions on a sub-nanosecond time scale revealed by an extended Lipari–Szabo analysis were mainly independent of the solvent viscosity. The lower limit of the NMR `observation window' for the internal correlation time shifts above 0.5 ns at 70% EG, which is directly reflected in the experimentally derived internal correlation times. Chemical exchange contributions to the transverse relaxation rates derived from the Lipari-Szabo approach coincide with the experimentally determined values from the transverse ¹H-¹⁵N dipolar/¹⁵N chemical shift anisotropy relaxation interference. These contributions originate from fast protein folding reactions on a millisecond timescale, which get retarded at increased solvent viscosities.     

Backbone dynamics of the regulatory domain of calcium vector protein, studied by 15N relaxation at four fields, reveals unique mobility characteristics of the intermotif linker

CaVP is a calcium-binding protein from amphioxus. It has a modular composition with two domains, but only the two EF-hand motifs localized in the C-terminal domain are functional. We recently determined the solution structure of this regulatory half (C-CaVP) in the Ca(2+)-saturated form and characterized the stepwise ion binding. This paper reports the (15)N nuclear relaxation rates of the Ca(2+)-saturated C-CaVP, measured at four different NMR fields (9.39, 11.74, 14.1, and 18.7 T), which were used to map the spectral density function for the majority of the amide H(N)-N vectors. Fitting the spectral density values at eight frequencies by a model-free approach, we obtained the microdynamic parameters characterizing the global and internal movements of the polypeptide backbone. The two EF-hand motifs, including the ion binding loops, behave like compact structural units with restricted mobility as reflected in the quite uniform order parameter and short internal correlation time (< 20 nsec). Comparative analysis of the two Ca(2+) binding sites shows that site III, having a larger affinity for the metal ion, is generally more rigid, and the amide vector in the second residue of each loop is significantly less restricted. The linker fragment is animated simultaneously by a larger amplitude fast motion and a slow conformational exchange on a microsecond to millisecond time scale. The backbone dynamics of C-CaVP characterized here is discussed in relation with other well-characterized Ca(2+)-binding proteins. Supplemental material: See www.proteinscience.org     

Accessibility of tobacco lipid transfer protein cavity revealed by 15N NMR relaxation studies and molecular dynamics simulations

Plant LTP1 are small helical proteins stabilized by four disulfide bridges and are characterized by the presence of an internal cavity, in which various hydrophobic ligands can be inserted. Recently, we have determined the solution structure of the recombinant tobacco LTP1_1. Unexpectedly, despite a global fold very similar to the structures already known for cereal seed LTP1, its binding properties are different: Tobacco LTP1_1 is able to bind only one monoacylated lipid, whereas cereal LTP1 can bind either one or two. The 3D structure of tobacco LTP1_1 revealed the presence of a hydrophobic cluster, not observed on cereal LTP1 structures, which may hinder one of the two entrances of the cavity defined for wheat LTP1. To better understand the mechanism of lipid entrance for tobacco LTP1_1 and to define the regions of the protein monitoring the accessibility of the cavity, we have complemented our structural data by the study of the internal dynamics of tobacco LTP1_1, using (15)N magnetic relaxation rate data and MD simulations at room and high temperatures. This work allowed us to define two regions of the protein experiencing the largest motions. These two regions delineate a portal that opens up during the simulation constituting a unique entrance of the hydrophobic cavity, in contrast with wheat LTP1 where two routes were detected. The hydrophobic interactions resulting from a few point mutations are strong enough to completely block the second portal so that the accessibility of the cavity is restricted to one entrance, explaining why this particular LTP1 binds only one lipid molecule.     

A comparative study of the backbone dynamics of two closely related lipid binding proteins: Bovine heart fatty acid binding protein and porcine ileal lipid binding protein

The backbone dynamics of bovine heart fatty acid binding protein (H-FABP) and porcine ileal lipid binding protein (ILBP) were studied by 15N NMR relaxation (T1 and T2) and steady state heteronuclear 15N{1H} NOE measurements. The microdynamic parameters characterizing the backbone mobility were determined using the model-free approach. For H-FABP, the non-terminal backbone amide groups display a rather compact protein structure of low flexibility. In contrast, for ILBP an increased number of backbone amide groups display unusually high internal mobility. Furthermore, the data indicate a higher degree of conformational exchange processes in the sec-msec time range for ILBP compared to H-FABP. These results suggest significant differences in the conformational stability for these two structurally highly homologous members of the fatty acid binding protein family.