Conformational and dynamics changes induced by bile acids binding to chicken liver bile acid binding protein

The correlation between protein motions and function is a central problem in protein science. Several studies have demonstrated that ligand binding and protein dynamics are strongly correlated in intracellular lipid binding proteins (iLBPs), in which the high degree of flexibility, principally occurring at the level of helix-II, CD, and EF loops (the so-called portal area), is significantly reduced upon ligand binding. We have recently investigated by NMR the dynamic properties of a member of the iLBP family, chicken liver bile acid binding protein (cL-BABP), in its apo and holo form, as a complex with two bile salts molecules. Binding was found to be regulated by a dynamic process and a conformational rearrangement was associated with this event. We report here the results of molecular dynamics (MD) simulations performed on apo and holo cL-BABP with the aim of further characterizing the protein regions involved in motion propagation and of evaluating the main molecular interactions stabilizing bound ligands. Upon binding, the root mean square fluctuation values substantially decrease for CD and EF loops while increase for the helix-loop-helix region, thus indicating that the portal area is the region mostly affected by complex formation. These results nicely correlate with backbone dynamics data derived from NMR experiments. Essential dynamics analysis of the MD trajectories indicates that the major concerted motions involve the three contiguous structural elements of the portal area, which however are dynamically coupled in different ways whether in the presence or in the absence of the ligands. Motions of the EF loop and of the helical region are part of the essential space of both apo and holo-BABP and sample a much wider conformational space in the apo form. Together with NMR results, these data support the view that, in the apo protein, the flexible EF loop visits many conformational states including those typical of the holo state and that the ligand acts stabilizing one of these pre-existing conformations. The present results, in agreement with data reported for other iLBPs, sharpen our knowledge on the binding mechanism for this protein family.     

Molecular dynamics simulations of apo and holo forms of fatty acid binding protein 5 and cellular retinoic acid binding protein II reveal highly mobile protein,retinoic acid ligand,and water molecules

Structural and dynamic properties from a series of 300 ns molecular dynamics, MD, simulations of two intracellular lipid binding proteins, iLBPs, (Fatty Acid Binding Protein 5, FABP5, and Cellular Retinoic Acid Binding Protein II, CRABP-II) in both the apo form and when bound with retinoic acid reveal a high degree of protein and ligand flexibility. The ratio of FABP5 to CRABP-II in a cell may determine whether it undergoes natural apoptosis or unrestricted cell growth in the presence of retinoic acid. As a result, FABP5 is a promising target for cancer therapy. The MD simulations presented here reveal distinct differences in the two proteins and provide insight into the binding mechanism. CRABP-II is a much larger, more flexible protein that closes upon ligand binding, where FABP5 transitions to an open state in the holo form. The traditional understanding obtained from crystal structures of the gap between two β-sheets of the β-barrel common to iLBPs and the α-helix cap that forms the portal to the binding pocket is insufficient for describing protein conformation (open vs. closed) or ligand entry and exit. When the high degree of mobility between multiple conformations of both the ligand and protein are examined via MD simulation, a new mode of ligand motion that improves understanding of binding dynamics is revealed.     

Opening and closing motions in the periplasmic vitamin B12 binding protein BtuF

BtuF is the periplasmic binding protein (PBP) in the vitamin B(12) uptake system in Escherichia coli where it is associated with the ABC transporter BtuCD. When the ligand binds, PBPs generally display large conformational changes, commonly termed the Venus flytrap mechanism. BtuF belongs to a group of PBPs that, on the basis of crystal structures, does not appear to display such behavior. Using 480 ns multicopy molecular dynamics simulations of apo and holo forms of the protein, we investigate the dynamics of BtuF. We find BtuF to be more flexible than previously assumed, displaying clear opening and closing motions which are more pronounced in the apo form. The protein behavior is compatible with a PBP functional model that postulates a closed conformation for the ligand-bound state, whereas the empty form fluctuates between open and closed conformations. Elastic network normal-mode analysis suggests that all BtuF-like PBPs are capable of similar opening and closing motions. It also makes the typical Venus flytrap domain motions a likely common means of how PBP-ABC transporter interaction could occur.     

Molecular motions and conformational changes of HPPK

6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) belongs to a class of catalytic enzymes involved in phosphoryl transfer and is a new target for the development of novel antimicrobial agents. In the present study, the fundamental consideration is to view the overall structure of HPPK as a network of interacting residues and to extract the most cooperative collective motions that define its global dynamics. A coarse-grained model, harmonically constrained according to HPPK's crystal structure is used. Four crystal structures of HPPK (one apo and three holo forms with different nucleotide and pterin analogs) are studied with the goal of providing insights about the function-dynamic correlation and ligand induced conformational changes. The dynamic differences are examined between HPPK's apo- and holo-forms, because they are involved in the catalytic reaction steps. Our results indicate that the palm-like structure of HPPK is nearly rigid, whereas the two flexible loops: L2 (residues 43-53) and L3 (residues 82-92) exhibit the most concerted motions for ligand recognition and presumably, catalysis. These two flexible loops are involved in the recognition of HPPKs nucleotide and pterin ligands, whereas the rigid palm region is associated with binding of these cognate ligands. Six domains of collective motions are identified, comprised of structurally close but not necessarily sequential residues. Two of these domains correspond to the flexible loops (L2 and L3), whereas the remaining domains correspond to the rigid part of the molecule.     

Prediction of interaction sites from apo 3D structures when the holo conformation is different

The predictability of catalytic and binding sites from apo structures is addressed for proteins that undergo significant conformational change upon binding. Theoretical microscopic titration curves (THEMATICS), an electrostatics-based method for the prediction of functional sites, is performed on a test set of 24 proteins with both apo and holo structures available. For 23 of these 24 proteins (96%), THEMATICS predicts the correct catalytic or binding site for both the apo and holo forms. For only one of the 24 proteins, THEMATICS makes the correct prediction for the holo structure but fails for the apo structure. The metrics used by THEMATICS to identify functional residues generally are larger in absolute value for the functional residues in the holo forms compared to the corresponding residues in the apo forms. However, even in the apo forms, these identifying metrics are still statistically significantly larger for functional residues than for residues not involved in catalysis or binding. This indicates that some of the unusual electrostatic properties of functional residues are preserved in the apo conformation. Evidence is presented that certain residues immediately surrounding the active catalytic and binding residues impart functionally important chemical and electrostatic properties to the active residues. At least parts of these microenvironments exist in the unbound conformations, such that THEMATICS is able to distinguish the functional residues even in the apo structures.     

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.     

Structure and dynamics of the fatty acid binding cavity in apo rat intestinal fatty acid binding protein.

The structure and dynamics of the fatty acid binding cavity in I-FABP (rat intestinal fatty acid binding protein) were analyzed. In the crystal structure of apo I-FABP, the probe occupied cavity volume and surface are 539+/-8 A3 and 428 A2, respectively (1.4 A probe). A total of 31 residues contact the cavity with their side chains. The side-chain cavity surface is partitioned according to the residue type as follows: 36-39% hydrophobic, 21-25% hydrophilic, and 37-43% neutral or ambivalent. Thus, the cavity surface is neither like a typical protein interior core, nor is like a typical protein external surface. All hydrophilic residues that contact the cavity-with the exception of Asp74-are clustered on the one side of the cavity. The cavity appears to expand its hydrophobic surface upon fatty acid binding on the side opposite to this hydrophilic patch. In holo I-FABP the fatty acid chain interactions with the hydrophilic side chains are mediated by water molecules. Molecular dynamics (MD) simulation of fully solvated apo I-FABP showed global conformational changes of I-FABP, which resulted in a large, but seemingly transient, exposure of the cavity to the external solvent. The packing density of the side chains lining the cavity, studied by Voronoi volumes, showed the presence of two distinctive small hydrophobic cores. The MD simulation predicts significant structural perturbations of the cavity on the subnanosecond time scale, which are capable of facilitating exchange of I-FABP internal water.     

Molecular dynamics simulations of the bacterial periplasmic heme binding proteins ShuT and PhuT

ShuT and PhuT are two periplasmic heme binding proteins that shuttle heme between the outer and inner membranes of the Gram-negative bacteria. Periplasmic binding proteins (PBPs) generally exhibit considerable conformational changes during the ligand binding process, whereas ShuT and PhuT belong to a class of PBPs that do not show such behavior based on their apo and holo crystal structures. By employing a series of molecular dynamic simulations on the ShuT and the PhuT, the dynamics and functions of the two PBPs were investigated. Through monitoring the distance changes between the two conserved glutamates of ShuT and PhuT, it was found the two PBPs were more flexible than previously assumed, exhibiting obvious opening-closing motions which were more remarkable in the apo runs of ShuT. Based on the results of the domain motion analysis, large scale conformational transitions were found in all apo runs of ShuT and PhuT, hinting that the domain motions of the two PBPs may be intrinsic. On the basis of the results of the principle component analysis, distinct opening-closing and twisting motion tendencies were observed not only in the apo, but also in the holo simulations of the two PBPs. The Gaussian network model was applied in order to analyze the hinge bending regions. The most important bending regions of ShuT and PhuT are located around the midpoints of their respective connecting helixes. Finally, the flexibilities and the details of the simulations of ShuT and PhuT were discussed. Characterized by the remarkably large flexibilities, the loop constituted by Ala 169, Gly170 and Gly171 of ShuT and the beta-turn constituted by Ala176, Gly177 and Gly178 of PhuT may be important for the functions of the two PBPs. Furthermore, the Asn254 of ShuT and the Arg228 of PhuT may be indispensable for the binding or unbinding of heme, since it is involved in the important hydrogen bonding to the propionate side-chains of heme.     

Internal motions and exchange processes in human ileal bile acid binding protein as studied by backbone (15)N nuclear magnetic resonance spectroscopy

Human ileal bile acid binding protein (I-BABP), a member of the family of intracellular lipid binding proteins, is thought to play a role in the enterohepatic circulation of bile salts. Previously, we have shown by stopped-flow fluorescence analysis that positive binding cooperativity exhibited by I-BABP in its interactions with glycocholate (GCA) and glycochenodeoxycholate (GCDA), the two primary bile salts in humans, is related to a slow conformational change in the protein. In this study, we used backbone (15)N relaxation nuclear magnetic resonance (NMR) techniques to obtain residue-specific information about the internal dynamics of apo I-BABP and the doubly ligated I-BABP:GCA:GCDA complex on various time scales. According to our NMR data, bile salt binding is accompanied by a slight rigidification of the (15)N-(1)H bond vectors on the picosecond to nanosecond time scale, with most pronounced changes occurring in the C-D region. In contrast to the minor effects of ligation on fast motions, relaxation dispersion NMR experiments indicate a marked difference between the two protein states on the microsecond to millisecond time scale. In the apo form, an extensive network of conformational fluctuations is detected throughout segments of the EFGHIJ β-strands and the C-D loop, which cease upon complexation. Our NMR data are in agreement with a conformational selection model we proposed earlier for I-BABP and support the hypothesis of an allosteric mechanism of ligand binding. According to the NMR measurements, the helical cap region may have a less crucial role in mediating ligand entry and release than what has been indicated for fatty acid binding proteins.     

Direct observation of protein-ligand interaction kinetics

Internal dynamics on the micro- to millisecond time scale have a strong influence on the affinity and specificity with which a protein binds ligands. This time scale is accessible through relaxation dispersion measurements using NMR. By studying the dynamics of a protein with different concentrations of a ligand, one can determine the dynamic effects induced by the ligand. Here we have studied slow internal dynamics of the N-terminal src homology 2 domain of phosphatidylinositide 3-kinase to probe the role of individual residues for the interaction with a tyrosine-phosphorylated binding sequence from polyoma middle T antigen. While slow dynamic motion was restricted to a few residues in the free SH2 and in the SH2 complex, motion was significantly enhanced by adding even small amounts of ligand. Kinetic rates induced by ligand binding varied between 300 and 2000 s(-1). High rates reflected direct interactions with the ligand or rearrangements caused by ligand binding. Large differences in rates were observed for residues adjacent in the primary sequence reflecting their individual roles in ligand interaction. However, rates were similar for residues involved in the same side chain interactions, reflecting concerted motions during ligand binding. For a subset of residues, exchange must involve structural intermediates which play a crucial role in high-affinity ligand binding. This analysis supports a new view of the dynamics of individual sites of a protein during ligand interaction.     

Protein function annotation by local binding site surface similarity

Hundreds of protein crystal structures exist for proteins whose function cannot be confidently determined from sequence similarity. Surflex‐PSIM, a previously reported surface‐based protein similarity algorithm, provides an alternative method for hypothesizing function for such proteins. The method now supports fully automatic binding site detection and is fast enough to screen comprehensive databases of protein binding sites. The binding site detection methodology was validated on apo/holo cognate protein pairs, correctly identifying 91% of ligand binding sites in holo structures and 88% in apo structures where corresponding sites existed. For correctly detected apo binding sites, the cognate holo site was the most similar binding site 87% of the time. PSIM was used to screen a set of proteins that had poorly characterized functions at the time of crystallization, but were later biochemically annotated. Using a fully automated protocol, this set of 8 proteins was screened against ～60,000 ligand binding sites from the PDB. PSIM correctly identified functional matches that predated query protein biochemical annotation for five out of the eight query proteins. A panel of 12 currently unannotated proteins was also screened, resulting in a large number of statistically significant binding site matches, some of which suggest likely functions for the podorly characterized proteins. Proteins 2014; 82:679–694. © 2013 Wiley Periodicals, Inc.     

Substrate binding modifies the hinge bending characteristics of human 3‐phosphoglycerate kinase: A molecular dynamics study

3‐Phosphogycerate kinase (PGK) is a two domain enzyme, with a binding site of the 1,3‐bisphosphoglycerate on the N‐domain and of the ADP on the C‐domain. To transfer a phosphate group the enzyme has to undergo a hinge bending motion from open to closed conformation to bring the substrates to close proximity. Molecular dynamics simulation was used to elucidate the effect of ligand binding onto the domain motions of this enzyme. The simulation results of the apo form indicate a hinge bending motion in the ns timescale while the time period of the hinge bending motion of the complex form is clearly over the 20 ns simulation time. The apo form exhibits several hinge points that contribute to the hinge bending motion while upon binding the ligands, the hinge bending becomes strictly restrained with one dominant hinge point in the vicinity of the substrates. At the same time, ligand binding results in an enhanced correlation of internal domain motions. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Molecular dynamics simulations reveal that apo‐HisJ can sample a closed conformation

The Escherichia coli histidine binding protein HisJ is a type II periplasmic binding protein (PBP) that preferentially binds histidine and interacts with its cytoplasmic membrane ABC transporter, HisQMP₂, to initiate histidine transport. HisJ is a bilobal protein where the larger Domain 1 is connected to the smaller Domain 2 via two linking strands. Type II PBPs are thought to undergo “Venus flytrap” movements where the protein is able to reversibly transition from an open to a closed conformation. To explore the accessibility of the closed conformation to the apo state of the protein, we performed a set of all‐atom molecular dynamics simulations of HisJ starting from four different conformations: apo‐open, apo‐closed, apo‐semiopen, and holo‐closed. The simulations reveal that the closed conformation is less dynamic than the open one. HisJ experienced closing motions and explored semiopen conformations that reverted to closed forms resembling those found in the holo‐closed state. Essential dynamics analysis of the simulations identified domain closing/opening and twisting as main motions. The formation of specific inter‐hinge strand and interdomain polar interactions contributed to the adoption of the closed apo‐conformations although they are up to 2.5‐fold less prevalent compared with the holo‐closed simulations. The overall sampling of the closed form by apo‐HisJ provides a rationale for the binding of unliganded PBPs with their cytoplasmic membrane ABC transporters. Proteins 2014; 82:386–398. © 2013 Wiley Periodicals, Inc.     

Molecular dynamics simulations of the holo and apo forms of retinol binding protein. Structural and dynamical changes induced by retinol removal

The effects of removing retinol from the X-ray structure of holo-retinol binding protein are studied using the molecular dynamics technique. Structural and dynamical properties emerging from an 80 ps simulation of the apo form, for which no crystallographic structure is available, are compared with the results of a 70 ps trajectory of the holo-protein. Dynamical stationarity is attained after roughly 30 ps, and the resulting average structure is proposed as a reasonable model of the apo-protein. Conformational changes are observed for the loops at the beta-barrel entrance during the non-equilibrium part of the apo-trajectory. Tryptophan labelling experiments and retinoid reconstitution experiments point towards this part of the molecule as being involved in prealbumin binding. Structural changes in this region may therefore explain the differences in prealbumin affinity between the apo and holo forms. Furthermore, a change in the position of the alpha-helix, corresponding to a pivot around its C terminus, is observed for the apo-protein. The resulting conformation of the alpha-helix is found to be similar to that in apo-beta-lactoglobulin, which also can bind retinol and for which a crystal structure exists. The results from the holo simulation are compared to the crystallographic data and show good agreement. The dynamics of the secondary and tertiary structural elements are analysed and compared for the two forms. The beta-barrel is found to be extremely cooperative in its atomic motions in both simulations, and the top and bottom beta-sheets perform collective fluctuations with respect to each other in the low-frequency limit of the simulations. The dynamics of the alpha-helical region presents clear differences between the two forms; while the holo-protein has a well-defined spectrum for the longitudinal stretching mode, the apo form displays a fairly large bending of the alpha-helix at several points of the trajectory.     

Perturbation-Response Scanning Reveals Ligand Entry-Exit Mechanisms of Ferric Binding Protein

We study apo and holo forms of the bacterial ferric binding protein (FBP) which exhibits the so-called ferric transport dilemma: it uptakes iron from the host with remarkable affinity, yet releases it with ease in the cytoplasm for subsequent use. The observations fit the “conformational selection” model whereby the existence of a weakly populated, higher energy conformation that is stabilized in the presence of the ligand is proposed. We introduce a new tool that we term perturbation-response scanning (PRS) for the analysis of remote control strategies utilized. The approach relies on the systematic use of computational perturbation/response techniques based on linear response theory, by sequentially applying directed forces on single-residues along the chain and recording the resulting relative changes in the residue coordinates. We further obtain closed-form expressions for the magnitude and the directionality of the response. Using PRS, we study the ligand release mechanisms of FBP and support the findings by molecular dynamics simulations. We find that the residue-by-residue displacements between the apo and the holo forms, as determined from the X-ray structures, are faithfully reproduced by perturbations applied on the majority of the residues of the apo form. However, once the stabilizing ligand (Fe) is integrated to the system in holo FBP, perturbing only a few select residues successfully reproduces the experimental displacements. Thus, iron uptake by FBP is a favored process in the fluctuating environment of the protein, whereas iron release is controlled by mechanisms including chelation and allostery. The directional analysis that we implement in the PRS methodology implicates the latter mechanism by leading to a few distant, charged, and exposed loop residues. Upon perturbing these, irrespective of the direction of the operating forces, we find that the cap residues involved in iron release are made to operate coherently, facilitating release of the ion.     

Structural and dynamic roles of permanent water molecules in ligand molecular recognition by chicken liver bile acid binding protein

Chicken liver bile acid binding protein (cL-BABP) crystallizes with water molecules in its binding site. To obtain insights on the role of internal water, we performed two 100 ns molecular dynamics (MD) simulations in explicit solvent for cL-BABP, as apo form and as a complex with two molecules of cholic acid, and analyzed in detail the dynamics properties of all water molecules. The diffusion coefficients of the more persistent internal water molecules are significantly different from the bulk, but similar between the two protein forms. A different number of molecules and a different organization are observed for apo- and holo-cL-BABP. Most water molecules identified in the binding site of the apo-crystal diffuse to the bulk during the simulation. In contrast, almost all the internal waters of the holo-crystal maintain the same interactions with internal sidechains and ligands, which suggests they have a relevant role in protein-ligand molecular recognition. Only in the presence of these water molecules we were able to reproduce, by a classical molecular docking approach, the structure of the complex cL-BABP::cholic acid with a low ligand root mean square deviation (RMSD) with respect to its reference positioning. Literature data reported a conserved pattern of hydrogen bonds between a single water molecule and three amino acid residues of the binding site in a series of crystallized FABP. In cL-BABP, the interactions between this conserved water molecule and the three residues are present in the crystal of both apo- and holo-cL-BABP but are lost immediately after the start of molecular dynamics. Copyright (c) 2008 John Wiley & Sons, Ltd.     

Evidence of chemical exchange in recombinant Major Urinary Protein and quenching thereof upon pheromone binding

The internal dynamics of recombinant Major Urinary Protein (rMUP) have been investigated by monitoring transverse nitrogen-15 relaxation using multiple-echo Carr–Purcell–Meiboom–Gill (CPMG) experiments. While the ligand-free protein (APO-rMUP) features extensive evidence of motions on the milliseconds time scale, the complex with 2-methoxy-3-isobutylpyrazine (HOLO-rMUP) appears to be much less mobile on this time scale. At 308 K, exchange rates k _ex = 500–2000 s⁻¹ were typically observed in APO-rMUP for residues located adjacent to a β-turn comprising residues 83–87. These residues occlude an entry to the binding pocket and have been proposed to be a portal for ligand entry in other members of the lipocalin family, such as the retinol binding protein and the human fatty-acid binding protein. Exchange rates and populations are largely uncorrelated, suggesting local ‘breathing’ motions rather than a concerted global conformational change. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Dynamics of the S1S2 glutamate binding domain of GluR2 measured using 19F NMR spectroscopy

Ionotropic glutamate receptors mediate the majority of vertebrate excitatory synaptic transmission. Although the structure of the GluR2 binding domain (S1S2) is well known (agonist binding site between two lobes), little is known about the time scales of conformational transitions or the relationship between dynamics and function. (19)F NMR ((19)F-labeled tryptophan) spectroscopy was used to monitor motions in the S1S2 domain bound to ligands with varying efficacy and in the apo state. One tryptophan (Trp-671) undergoes chemical exchange in some but not all agonists, consistent with mus-ms motion. The dynamics can be correlated to ligand affinity, and a likely source of the motion is a peptide bond capable of transiently forming hydrogen bonds across the lobe interface. Another tryptophan (Trp-767) appears to monitor motions of the relative positions of the lobes and suggests that the relative orientation in the apo- and antagonist-bound forms can exchange between at least two conformations on the ms time scale.     

Conformational selection in silico: loop latching motions and ligand binding in enzymes

Ligand binding frequently induces significant conformational changes in a protein receptor. Understanding and predicting such conformational changes represent an important challenge for computational biology, including applications to structure-based drug design. We describe an approach to this problem based on the assumption that the holo state is at least transiently populated in the absence of a ligand; this hypothesis has been referred to as "conformational selection." Here, we apply a method that tests this hypothesis on a challenging class of ligand-induced conformational changes, which we refer to as loop latching: the closing of a loop around an active site that sequesters the ligand from solvent. The method uses a combination of replica exchange molecular dynamics and a loop prediction algorithm to generate low-energy loop structures, and docking to select the conformation appropriate for binding a particular ligand. On a test set of six proteins, it yields loop structures including hololike conformations, generally below 2 A RMSD from the liganded structure, for loops that span up to 15 residues. Docking serves as a stringent test of the predictions. In five of the six cases, the predicted loop conformations improve the ranks of cognate ligands relative to using the apo structure, although the results remain, in most cases, significantly worse than using a holo structure. The poses of the cognate ligands are correct in four of the six test cases, while they are correct for five of the six using a holo structure.     

Effect of cofactor binding and loop conformation on side chain methyl dynamics in dihydrofolate reductase

Dihydrofolate reductase (DHFR) has several flexible active site loops that facilitate ligand binding and catalysis. Previous studies of backbone dynamics in several complexes of DHFR indicate that the time scale and amplitude of motion depend on the conformation of the active site loops. In this study, information on dynamics is extended to methyl-containing side chains. To understand the role of side chain dynamics in ligand binding and loop conformation, methyl deuterium relaxation rates of Escherichia coli DHFR in binary folate and ternary folate:NADP+ complexes have been measured, together with chi(1) rotamer populations for threonine, isoleucine, and valine residues, determined from measurements of 3J(CgammaCO) and 3J(CgammaN) coupling constants. The results indicate that, in addition to backbone motional restriction in the adenosine-binding site, side chain flexibility in the active site and the surrounding active site loops is diminished upon binding NADP+. Resonances for several methyls in the active site and the surrounding active site loops were severely broadened in the folate:NADP+ ternary complex, suggesting the presence of motion on the chemical shift time scale. The side chains of Ile14 and Ile94, which pack against the nicotinamide and pterin rings of the cofactor and substrate, respectively, exhibit rotamer disorder in the ternary folate:NADP+ complex. Conformational fluctuations of these side chains may play a role in transition state stabilization; the observed line broadening for Ile14 suggests motions on a microsecond/millisecond time scale.