Temperature dependence of fast dynamics in proteins

The temperature dependence of the internal dynamics of recombinant human ubiquitin has been measured using solution NMR relaxation techniques. Nitrogen-15 relaxation has been employed to obtain a measure of the amplitude of subnanosecond motion at amide N-H sites in the protein. Deuterium relaxation has been used to obtain a measure of the amplitude of motion of methyl-groups in amino-acid side chains. Data was obtained between 5 and 55 degrees C. The majority of amide N-H and methyl groups show a roughly linear (R(2)>0.75) temperature dependence of the associated Lipari-Szabo model-free squared generalized-order parameter (O(2)) describing the amplitude of motion. Interestingly, for those sites showing a linear response, the temperature dependence of the backbone is distinct from that of the methyl-bearing side chains with the former being characterized by a significantly larger Lambda-value, where Lambda is defined as d ln(1 - O)/d lnT. These results are comparable to the sole previous such study of the temperature dependence of protein motion obtained for a calmodulin-peptide complex. This suggests that the distinction between the main chain and methyl-bearing side chains may be general. Insight into the temperature dependence is gathered from a simple two-state step potential model.     

Characterization of the backbone and side chain dynamics of the CaM-CaMKIp complex reveals microscopic contributions to protein conformational entropy

Calmodulin is a central mediator of calcium-dependent signal transduction pathways and regulates the activity of a large number of diverse targets. Calcium-dependent interactions of calmodulin with regulated proteins are of generally high affinity but of quite variable thermodynamic origins. Here we investigate the influence of the binding of the calmodulin-binding domain of calmodulin kinase I on the fast internal dynamics of calcium-saturated calmodulin. NMR relaxation was used to probe motion on the backbone (viewed through the backbone amide NH group) and the side chains (viewed through methyl groups). The distribution of the amplitudes of side chain dynamics is trimodal. The microscopic details of side chain motion are compared with those of a thermodynamically and structurally similar complex of calmodulin with the calmodulin-binding domain of the smooth muscle myosin light chain kinase. While there are no significant differences in backbone dynamics and no net change in methyl-bearing side chain dynamics, a large redistribution of the amplitude of methyl dynamics is observed between the two complexes. The variation in dynamics was largely localized to the heterogeneously dynamic target-binding interface, suggesting that differential dynamics of the binding surface plays a functional role in the high-affinity binding interactions of calmodulin. These results begin to reveal a fundamental role for residual protein entropy in molecular recognition by calmodulin.     

An increase in side chain entropy facilitates effector binding: NMR characterization of the side chain methyl group dynamics in Cdc42Hs

Cdc42Hs is a signal transduction protein that is involved in cytoskeletal growth and organization. We describe here the methyl side chain dynamics of three forms of (2)H,(13)C,(15)N-Cdc42Hs [GDP-bound (inactive), GMPPCP-bound (active), and GMPPCP/PBD46-bound (effector-bound)] from (13)C-(1)H NMR measurements of deuterium T(1) and T(1 rho) relaxation times. A wide variation in flexibility was observed throughout the protein, with methyl axis order parameters (S(2)(axis)) ranging from 0.2 to 0.4 (highly disordered) in regions near the PBD46 binding site to 0.8--1.0 (highly ordered) in some helices. The side chain dynamics of the GDP and GMPPCP forms are similar, with methyl groups on the PBD46 binding surface experiencing significantly greater mobility (lower S(2)(axis)) than those not on the binding surface. Binding of PBD46 results in a significant increase in the disorder and a corresponding increase in entropy for the majority of methyl groups. Many of the methyl groups that experience an increase in mobility are found in residues that are not part of the PBD46 binding interface. This entropy gain represents a favorable contribution to the overall entropy of effector binding and partially offsets unfavorable entropy losses such as those that occur in the backbone.     

Temperature dependence of the internal dynamics of a calmodulin-peptide complex

The temperature dependence of the fast internal dynamics of calcium-saturated calmodulin in complex with a peptide corresponding to the calmodulin-binding domain of the smooth muscle myosin light chain kinase is examined using 15N and 2H NMR relaxation methods. NMR relaxation studies of the complex were carried out at 13 temperatures that span 288-346 K. The dynamics of the backbone and over four dozen methyl-bearing side chains, distributed throughout the calmodulin molecule, were probed. The side chains show a much more variable and often considerably larger response to temperature than the backbone. A significant variation in the temperature dependence of the amplitude of motion of individual side chains is seen. The amplitude of motion of some side chains is essentially temperature-independent while many show a simple roughly linear temperature dependence. In a few cases, angular order increases with temperature, which is interpreted as arising from interactions with neighboring residues. In addition, a number of side chains display a nonlinear temperature dependence. The significance of these and other results is illuminated by several simple interpretative models. Importantly, analysis of these models indicates that changes in generalized order parameters can be robustly related to corresponding changes in residual entropy. A simple cluster model that incorporates features of cooperative or conditional motion reproduces many of the unusual features of the experimentally observed temperature dependence and illustrates that side chain interactions result in a dynamically changing environment that significantly influences the motion of internal side chains. This model also suggests that the intrinsic entropy of interacting clusters of side chains is only modestly reduced from that of independent side chain motion. Finally, estimates of protein heat capacity support the view that the major contribution to the heat capacity of protein solutions largely arises from local bond vibrations and solvent interactions and not from torsional oscillations of side chains.     

Fast internal main-chain dynamics of human ubiquitin.

The fast internal dynamics of human ubiquitin have been studied by the analysis of 15N relaxation of backbone amide nitrogens. The amide 15N resonances have been assigned by use of heteronuclear multiple-quantum spectroscopy. Spin lattice relaxation times at 60.8 and 30.4 MHz and the steady-state nuclear Overhauser effect at 60.8 MHz have been determined for 67 amide 15N sites in the protein using two-dimensional spectroscopy. These data have been analyzed in terms of the model free treatment of Lipari and Szabo [Lipari, G., & Szabo, A. (1982) J. Am. Chem. Soc. 104, 4546-4559]. The global motion of the protein is shown to be isotropic and is characterized by a correlation time of 4.1 ns rad-1. The generalized order parameters (S2) of backbone amide N-H vectors in the globular region of the protein range from 0.5 to 0.95. No apparent correlation between secondary structure and generalized order parameters is observed. There is, however, a strong correlation between the magnitude of the generalized order parameters of a given N-H vector and the presence of hydrogen bonding of the amide hydrogen or its peptide bond associated carbonyl. Using a chemical shift tensor breadth of 160 ppm, the N-H vectors of peptide linkages participating in one or more hydrogen bonds to the main chain show an average generalized order parameter of 0.80 (SD 0.06), while those amide NH of peptide linkages free of hydrogen-bonding interactions with the main chain show an average order parameter of 0.69 (SD 0.06).(ABSTRACT TRUNCATED AT 250 WORDS)     

Mapping backbone dynamics in solution with site-directed spin labeling: GCN4-58 bZip free and bound to DNA

In site-directed spin labeling, a nitroxide-containing side chain is introduced at selected sites in a protein. The EPR spectrum of the labeled protein encodes information about the motion of the nitroxide on the nanosecond time scale, which has contributions from the rotary diffusion of the protein, from internal motions in the side chain, and from backbone fluctuations. In the simplest model for the motion of noninteracting (surface) side chains, the contribution from the internal motion is sequence independent, as is that from protein rotary diffusion. Hence, differences in backbone motions should be revealed by comparing the sequence-dependent motions of nitroxides at structurally homologous sites. To examine this model, nitroxide side chains were introduced, one at a time, along the GCN4-58 bZip sequence, for which NMR (15)N relaxation experiments have identified a striking gradient of backbone mobility along the DNA-binding region [Bracken et al. (1999) J. Mol. Biol. 285, 2133]. Spectral simulation techniques and a simple line width measure were used to extract dynamical parameters from the EPR spectra, and the results reveal a mobility gradient similar to that observed in NMR relaxation, indicating that side chain motions mirror backbone motions. In addition, the sequence-dependent side chain dynamics were analyzed in the DNA/protein complex, which has not been previously investigated by NMR relaxation methods. As anticipated, the backbone motions are damped in the DNA-bound state, although a gradient of motion persists with residues at the DNA-binding site being the most highly ordered, similar to those of helices on globular proteins.     

Hydration dependence of backbone and side chain polylysine dynamics: a 13C solid-state NMR and IR spectroscopy study

The molecular dynamics of solid poly-L-lysine has been studied by the following natural abundance (13)C-NMR relaxation methods: measurements of the relaxation times T(1) at two resonance frequencies, off-resonance T(1rho) at two spin-lock frequencies, and proton-decoupled T(1rho). Experiments were performed at different temperatures and hydration levels (up to 17% H(2)O by weight). The natural abundance (13)C-CPMAS spectrum of polylysine provides spectral resolution of all types of backbone and side chain carbons and thus, dynamic parameters could be determined separately for each of them. At the same time, the conformational properties of polylysine were investigated by Fourier transform infrared spectroscopy. The data obtained from the different NMR experiments were simultaneously analyzed using the correlation function formalism and model-free approach. The results indicate that in dry polylysine both backbone and side chains take part in two low amplitude motions with correlation times of the order of 10(-4) s and 10(-9) s. Upon hydration, the dynamic parameters of the backbone remain almost constant except for the amplitude of the slower process that increases moderately. The side chain dynamics reveals a much stronger hydration response: the amplitudes of both slow and fast motions increase significantly and the correlation time of the slow motion shortens by about five orders of magnitude, and at hydration levels of more than 10% H(2)O fast and slow side chain motions are experimentally indistinguishable. These changes in the molecular dynamics cannot be ascribed to any hydration-dependent conformational transitions of polylysine because IR spectra reveal almost no hydration dependence in either backbone or side chain absorption domains. The physical nature of the fast and slow motions, their correlation time distributions, and hydration dependence of microdynamic parameters are discussed.     

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²) 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² 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² of 0.8.     

Increased rigidity of eglin c at acidic pH: evidence from NMR spin relaxation and MD simulations

To gain physical insights into how proteins respond to changes in pH, the picosecond to nanosecond time scale dynamics of the small serine protease inhibitor eglin c have been studied by NMR spin relaxation experiments and MD simulations under two pH solution conditions, pH 7 and 3. Like many proteins, eglin c is destabilized by a lowering of the pH, although it retains enough stability to maintain its native conformation at pH 3. Backbone (15)N relaxation results show comparable global tumbling times (tau(m)) and model-free order parameters (S(2)) under the two pH conditions, indicating that the molecule maintains its overall molecular shape and structure at low pH, although the backbone rigidity is slightly increased (/ = 0.6%). In contrast, the side-chain methyl dynamics, as measured from (2)H relaxation experiments, show a substantial increase in rigidity at lower pH (/ = 14.8%). Molecular dynamics simulations performed at these pH states produce results consistent with NMR measurements, showing that the two methods are in qualitative agreement. Although a full accounting of the physical basis for the concurrent conformational rigidification and destabilization at low pH requires further investigation, the high level of detail in the MD simulations provides a potential molecular mechanism: the breaking of the hydrogen bond between the side chains of Asp46 and Arg53, and changes in electrostatic interactions, appear to allow the binding loop to move closer to the core part of the protein, resulting in a more compact structure at low pH. This more compact structure may be responsible for the increased level of restriction of molecular motion. As these findings show, the stability of a molecular structure is distinct from its conformational rigidity, and the two can even change in opposite directions, against na?ve expectation.     

Side chain and backbone dynamics of phospholamban in phospholipid bilayers utilizing 2H and 15N solid-state NMR spectroscopy

2H and 15N solid-state NMR spectroscopic techniques were used to investigate both the side chain and backbone dynamics of wild-type phospholamban (WT-PLB) and its phosphorylated form (P-PLB) incorporated into 1-palmitoyl-2-oleoyl-sn-glycerophosphocholine (POPC) phospholipid bilayers. 2H NMR spectra of site-specific CD3-labeled WT-PLB (at Leu51, Ala24, and Ala15) in POPC bilayers were similar under frozen conditions (-25 degrees C). However, significant differences in the line shapes of the 2H NMR spectra were observed in the liquid crystalline phase at and above 0 degrees C. The 2H NMR spectra indicate that Leu51, located toward the lower end of the transmembrane (TM) helix, shows restricted side chain motion, implying that it is embedded inside the POPC lipid bilayer. Additionally, the line shape of the 2H NMR spectrum of CD3-Ala24 reveals more side chain dynamics, indicating that this residue (located in the upper end of the TM helix) has additional backbone and internal side chain motions. 2H NMR spectra of both WT-PLB and P-PLB with CD3-Ala15 exhibit strong isotropic spectral line shapes. The dynamic isotropic nature of the 2H peak can be attributed to side chain and backbone motions to residues located in an aqueous environment outside the membrane. Also, the spectra of 15N-labeled amide WT-PLB at Leu51 and Leu42 residues showed only a single powder pattern component indicating that these two 15N-labeled residues located in the TM helix are motionally restricted at 25 degrees C. Conversely, 15N-labeled amide WT-PLB at Ala11 located in the cytoplasmic domain showed both powder and isotropic components at 25 degrees C. Upon phosphorylation, the mobile component contribution increases at Ala11. The 2H and 15N NMR data indicate significant backbone motion for the cytoplasmic domain of WT-PLB when compared to the transmembrane section.     

Self-consistent residual dipolar coupling based model-free analysis for the robust determination of nanosecond to microsecond protein dynamics

Residual dipolar couplings (RDCs) provide information about the dynamic average orientation of inter-nuclear vectors and amplitudes of motion up to milliseconds. They complement relaxation methods, especially on a time-scale window that we have called supra-tau(c) (tau(c) < supra-tau(c) < 50 micros). Here we present a robust approach called Self-Consistent RDC-based Model-free analysis (SCRM) that delivers RDC-based order parameters-independent of the details of the structure used for alignment tensor calculation-as well as the dynamic average orientation of the inter-nuclear vectors in the protein structure in a self-consistent manner. For ubiquitin, the SCRM analysis yields an average RDC-derived order parameter of the NH vectors 0.72 +/- 0.02 compared to = 0.778 +/- 0.003 for the Lipari-Szabo order parameters, indicating that the inclusion of the supra-tau(c) window increases the averaged amplitude of mobility observed in the sub-supra-tau(c) window by about 34%. For the beta-strand spanned by residues Lys48 to Leu50, an alternating pattern of backbone NH RDC order parameter S2(rdc)(NH) = (0.59, 0.72, 0.59) was extracted. The backbone of Lys48, whose side chain is known to be involved in the poly-ubiquitylation process that leads to protein degradation, is very mobile on the supra-tau(c) time scale (S2(rdc)(NH) = 0.59 +/- 0.03), while it is inconspicuous (S2(LS)(NH)= 0.82) on the sub-tau(c) as well as on micros-ms relaxation dispersion time scales. The results of this work differ from previous RDC dynamics studies of ubiquitin in the sense that the results are essentially independent of structural noise providing a much more robust assessment of dynamic effects that underlie the RDC data.     

Comparison of backbone dynamics of oxidized and reduced putidaredoxin by 15N NMR relaxation measurements.

The backbone dynamics of uniformly 15N-labeled reduced and oxidized putidaredoxin (Pdx) have been studied by 2D 15N NMR relaxation measurements. 15N T1 and T2 values and 1H-15N NOEs have been measured for the diamagnetic region of the protein. These data were analyzed by using a model-free dynamics formalism to determine the generalized order parameters (S2), the effective correlation time for internal motions (tau e), and the 15N exchange broadening contributions (Rex) for each residue, as well as the overall correlation time (tau(m)). Order parameters for the reduced Pdx are generally higher than for the oxidized Pdx, and there is increased mobility on the microsecond to millisecond time scale for the oxidized Pdx, in comparison with the reduced Pdx. These results clearly indicate that the oxidized protein exhibits higher mobility than the reduced one, which is in agreement with the recently published redox-dependent dynamics studied by amide proton exchange. In addition, we observed very high T1/T2 ratios for residues 33 and 34, giving rise to a large Rex contribution. Residue 34 is believed to be involved in the binding of Pdx to cytochrome P450cam (CYP101). The differences in the backbone dynamics are discussed in relation to the oxidation states of Pdx, and their impact on electron transfer. The entropy change occurring on oxidation of reduced Pdx has been calculated from the order parameters of the two forms.     

Backbone dynamics measurements on leukemia inhibitory factor, a rigid four-helical bundle cytokine

The backbone dynamics of the four-helical bundle cytokine leukemia inhibitory factor (LIF) have been investigated using 15N NMR relaxation and amide proton exchange measurements on a murine-human chimera, MH35-LIF. For rapid backbone motions (on a time scale of 10 ps to 100 ns), as probed by 15N relaxation measurements, the dynamics parameters were calculated using the model-free formalism incorporating the model selection approach. The principal components of the inertia tensor of MH35-LIF, as calculated from its NMR structure, were 1:0.98:0.38. The global rotational motion of the molecule was, therefore, assumed to be axially symmetric in the analysis of its relaxation data. This yielded a diffusion anisotropy D(parallel)/D(perpendicular) of 1.31 and an effective correlation time (4D(perpendicular) + 2D(parallel))(-1) of 8.9 ns. The average values of the order parameters (S2) for the four helices, the long interhelical loops, and the N-terminus were 0.91, 0.84, and 0.65, respectively, indicating that LIF is fairly rigid in solution, except at the N-terminus. The S2 values for the long interhelical loops of MH35-LIF were higher than those of their counterparts in short-chain members of the four-helical bundle cytokine family. Residues involved in LIF receptor binding showed no consistent pattern of backbone mobilities, with S2 values ranging from 0.71 to 0.95, but residues contributing to receptor binding site III had relatively lower S2 values, implying higher amplitude motions than for the backbone of sites I and II. In the relatively slow motion regime, backbone amide exchange measurements showed that a number of amides from the helical bundle exchanged extremely slowly, persisting for several months in 2H2O at 37 degrees C. Evidence for local unfolding was considered, and correlations among various structure-related parameters and the backbone amide exchange rates were examined. Both sets of data concur in showing that LIF is one of the most rigid four-helical bundle cytokines.     

Dynamics of methyl groups in proteins as studied by proton-detected 13C NMR spectroscopy. Application to the leucine residues of staphylococcal nuclease.

This paper describes the application of recently developed nuclear magnetic resonance (NMR) pulse sequences to obtain information about the internal dynamics of isotopically enriched hydrophobic side chains in proteins. The two-dimensional spectra provided by the pulse sequences enable one to make accurate measurements of nuclear Overhauser effects (NOE) and longitudinal (T1) and transverse (T2) relaxation times of enriched methyl carbons in proteins. Herein, these techniques are used to investigate the internal dynamics of the 11 leucine side chains of staphylococcal nuclease (SNase), a small enzyme having Mr = 16.8K, in the absence and presence of ligands thymidine 3',5'-bisphosphate (pdTp) and Ca2+. We report the synthesis of [5,5'-13C2]leucine, the preparation of SNase containing the labeled leucine, the sequential assignment of the leucine methyl carbons and protons in the liganded and unliganded proteins, and the measurement of the 13C T1, T2, and NOE values for the SNase leucine methyl carbons. Analysis of the relaxation parameters using the formalism of Lipari and Szabo shows that the internal motions of the leucine methyl carbons are characterized by effective correlation times tau f (5-80 ps) and tau s (less than 2 ns). The fast motion is identified with the rapid rotation of the methyl group about the C gamma-C delta bond axis, while the slow motion is associated with reorientation of the C gamma-C delta bond axis itself. The mean squared order parameters associated with the latter motion, Ss2, lie in the range 0.34-0.92. The values of Ss2 correlate reasonably well with the temperature factors of the leucine methyl carbons obtained from the crystal structures, but some are smaller than anticipated on the basis of the fact that nearly all leucine methyl carbons are buried and have temperature factors no larger than that of the leucine backbone atoms. Five leucine residues in liganded SNase and eight in unliganded SNase have values of Ss2 less than 0.71. These order parameters correspond to large amplitude motions (angular excursions of 27-67 degrees) of the C gamma-C delta bond axis. These results indicate that, in solution, the internal motions of the leucine side chains of SNase are significantly larger than suggested by the X-ray structures or by qualitative analysis of NOESY spectra. Comparison of Ss2 values obtained from liganded and unliganded SNase reveals a strong correlation between delta Ss2 and distance between the leucine methyl carbon and the ligands.(ABSTRACT TRUNCATED AT 400 WORDS)     

Characterization of threonine side chain dynamics in an antifreeze protein using natural abundance 13C NMR spectroscopy

The dynamics of threonine side chains of the Tenebrio molitor antifreeze protein (TmAFP) were investigated using natural abundance (13)C NMR. In TmAFP, the array of threonine residues on one face of the protein is responsible for conferring its ability to bind crystalline ice and inhibit its growth. Heteronuclear longitudinal and transverse relaxation rates and the [(1)H]-(13)C NOE were determined in this study. The C alpha H relaxation measurements were compared to the previously measured (15)N backbone parameters and these are found to be in agreement. For the analysis of the threonine side chain motions, the model of restricted rotational diffusion about the chi(1) dihedral angle was employed [London and Avitabile (1978) J. Am. Chem. Soc., 100, 7159-7165]. We demonstrate that the motion experienced by the ice binding threonine side chains is highly restricted, with an approximate upper limit of less than +/-25 degrees.     

Main chain and side chain dynamics of the ubiquitin conjugating enzyme variant human Mms2 in the free and ubiquitin-bound States

Protein ubiquitination involves a cascade of enzymatic steps where ubiquitin (Ub) is sequentially transferred as a thiolester intermediate from an E1 enzyme to an E2 enzyme and finally to the protein target with the help of a Ub-protein ligase. Protein ubiquitination brought about by the Ubc13-Mms2 (E2-E2) complex has a unique role in the cell, unrelated to protein degradation. The Mms2-Ubc13 heterodimer links Ub molecules to one another through an isopeptide bond between its own C-terminus and Lys-63 on another Ub. The role of Mms2 is to orient a target-bound Ub molecule such that its Lys-63 is proximal to the C-terminus of the Ub molecule that is covalently linked to the active site of Ubc13. To gain insight into the influence of protein dynamics on the affinity of Ub for Mms2, we have determined pico- to nanosecond time scale fluctuations of the main chain and methyl side chains of human Mms2 in the free and Ub-bound states using solution state (15)N and (2)H nuclear magnetic resonance relaxation measurements. Analysis of the relaxation data allows for a semiquantitative estimation of the conformational entropy change (TDeltaS(NMR)) for the main chain and side chain methyl groups of Mms2 upon binding Ub. The value of TDeltaS(NMR) for the main chain and side chain methyl groups of Mms2 is -8 +/- 2 and -2 +/- 2 kcal mol(-)(1), respectively. The experimental DeltaG(binding) for the Mms2.Ub complex is -6 kcal mol(-)(1). Estimation of DeltaG(binding) using an empirical structure-based approach that does not account for changes in main chain entropy yields a value of -17 +/- 2 kcal mol(-)(1). However, inclusion of TDeltaS(NMR) for the main chain of Mms2 increases the estimated DeltaG(binding) to -9 +/- 3 kcal mol(-)(1). Assuming that changes in Ub main chain dynamics contribute to TDeltaS(NMR) to the same extent as Mms2, the estimated DeltaG(binding) is further reduced to -1 +/- 4 kcal mol(-)(1), a value close to the experimental DeltaG(binding).     

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.     

Structure of the oxidized long-chain flavodoxin from Anabaena 7120 at 2 A resolution.

The structure of the long-chain flavodoxin from the photosynthetic cyanobacterium Anabaena 7120 has been determined at 2 A resolution by the molecular replacement method using the atomic coordinates of the long-chain flavodoxin from Anacystis nidulans. The structure of a third long-chain flavodoxin from Chondrus crispus has recently been reported. Crystals of oxidized A. 7120 flavodoxin belong to the monoclinic space group P2(1) with a = 48.0, b = 32.0, c = 51.6 A, and beta = 92 degrees, and one molecule in the asymmetric unit. The 2 A intensity data were collected with oscillation films at the CHESS synchrotron source and processed to yield 9,795 independent intensities with Rmerg of 0.07. Of these, 8,493 reflections had I > 2 sigma and were used in the analysis. The model obtained by molecular replacement was initially refined by simulated annealing using the XPLOR program. Repeated refitting into omit maps and several rounds of conjugate gradient refinement led to an R-value of 0.185 for a model containing atoms for protein residues 2-169, flavin mononucleotide (FMN), and 104 solvent molecules. The FMN shows many interactions with the protein with the isoalloxazine ring, ribityl sugar, and the 5'-phosphate. The flavin ring has its pyrimidine end buried into the protein, and the functional dimethyl benzene edge is accessible to solvent. The FMN interactions in all three long-chain structures are similar except for the O4' of the ribityl chain, which interacts with the hydroxyl group of Thr 88 side chain in A. 7120, while with a water molecule in the other two. The phosphate group interacts with the atoms of the 9-15 loop as well as with NE1 of Trp 57. The N5 atom of flavin interacts with the amide NH of Ile 59 in A. 7120, whereas in A. nidulans it interacts with the amide NH of Val 59 in a similar manner. In C. crispus flavodoxin, N5 forms a hydrogen bond with the side chain hydroxyl group of the equivalent Thr 58. The hydrogen bond distances to the backbone NH groups in the first two flavodoxins are 3.6 A and 3.5 A, respectively, whereas in the third flavodoxin the distance is 3.1 A, close to the normal value. Even though the hydrogen bond distances are long in the first two cases, still they might have significant energy because their microenvironment in the protein is not accessible to solvent.(ABSTRACT TRUNCATED AT 400 WORDS)