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

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

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.     

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.     

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)     

Dynamics of a de novo designed three-helix bundle protein studied by 15N, 13C, and 2H NMR relaxation methods

Understanding how the amino acid sequence of a polypeptide chain specifies a unique, functional three-dimensional structure remains an important goal, especially in the context of the emerging discipline of de novo protein design. Alpha3D is a single chain protein of 73 amino acids resulting from a de novo design effort. Previous solution nuclear magnetic resonance studies of alpha3D confirm that the protein adopts the designed structure of a three-helix bundle. Furthermore, alpha3D has been previously shown to possess all of the major thermodynamic and structural characteristics of natural proteins, though it shares no sequence homology to any protein sequence in the database. In this work, the backbone and side-chain dynamics of alpha3D were investigated using 15N, 13C, and 2H nuclear magnetic resonance relaxation methods with the aim of assessing the character of the internal motions of this native-like protein of de novo design. At the backbone level, both 15N and 13C(alpha) relaxation studies indicate highly restrictive motion on the picosecond to nanosecond time scale in the alpha-helical regions of alpha3D, with increasing mobility at the ends of the alpha-helices and in the two loop regions. This is largely consistent with what is seen in proteins of natural origin. Overall, the view provided by both 2H and 13C methyl relaxation methods suggest that the side chains of alpha3D are more dynamic compared to natural proteins. Regions of relative flexibility bound clusters of rigid methyl-bearing side-chain groups that are interspersed with aromatic and beta-branched amino acids. The time scale of motions associated with methyl-bearing side chains of alpha3D are significantly longer than that seen in natural proteins. These results indicate that the strategies underlying the design of alpha3D have largely, but not completely, captured both the structural and dynamic character of natural proteins.     

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.     

Monitoring aromatic picosecond to nanosecond dynamics in proteins via 13C relaxation: expanding perturbation mapping of the rigidifying core mutation, V54A, in eglin c

Long-range effects, such as allostery, have evolved in proteins as a means of regulating function via communication between distal sites. An NMR-based perturbation mapping approach was used to more completely probe the dynamic response of the core mutation V54A in the protein eglin c by monitoring changes in picosecond to nanosecond aromatic side-chain dynamics and H/D exchange stabilities. Previous side-chain dynamics studies on this mutant were limited to methyl-bearing residues, most of which were found to rigidify on the picosecond to nanosecond time scale in the form of a contiguous "network". Here, high precision (13)C relaxation data from 13 aromatic side chains were acquired by applying canonical relaxation experiments to a newly developed carbon labeling scheme [Teilum et al. (2006) J. Am. Chem. Soc. 128, 2506-2507]. The fitting of model-free parameters yielded S (2) variability which is intermediate with respect to backbone and methyl-bearing side-chain variability and tau e values that are approximately 1 ns. Inclusion of the aromatic dynamic response results in an expanded network of dynamically coupled residues, with some aromatics showing increases in flexibility, which partially offsets the rigidification in methyl side chains. Using amide hydrogen exchange, dynamic propagation on a slower time scale was probed in response to the V54A perturbation. Surprisingly, regional stabilization (slowed exchange) 10-12 A from the site of mutation was observed despite a global destabilization of 1.5 kcal x mol (-1). Furthermore, this unlikely pocket of stabilized residues colocalizes with increases in aromatic flexibility on the faster time scale. Because the converse is also true (destabilized residues colocalize with rigidification on the fast time scale), a plausible entropy-driven mechanism is discussed for relating colocalization of opposing dynamic trends on vastly different time scales.     

Conformational exchange of aromatic side chains characterized by L-optimized TROSY-selected 13C CPMG relaxation dispersion

Protein dynamics on the millisecond time scale commonly reflect conformational transitions between distinct functional states. NMR relaxation dispersion experiments have provided important insights into biologically relevant dynamics with site-specific resolution, primarily targeting the protein backbone and methyl-bearing side chains. Aromatic side chains represent attractive probes of protein dynamics because they are over-represented in protein binding interfaces, play critical roles in enzyme catalysis, and form an important part of the core. Here we introduce a method to characterize millisecond conformational exchange of aromatic side chains in selectively (13)C labeled proteins by means of longitudinal- and transverse-relaxation optimized CPMG relaxation dispersion. By monitoring (13)C relaxation in a spin-state selective manner, significant sensitivity enhancement can be achieved in terms of both signal intensity and the relative exchange contribution to transverse relaxation. Further signal enhancement results from optimizing the longitudinal relaxation recovery of the covalently attached (1)H spins. We validated the L-TROSY-CPMG experiment by measuring fast folding-unfolding kinetics of the small protein CspB under native conditions. The determined unfolding rate matches perfectly with previous results from stopped-flow kinetics. The CPMG-derived chemical shift differences between the folded and unfolded states are in excellent agreement with those obtained by urea-dependent chemical shift analysis. The present method enables characterization of conformational exchange involving aromatic side chains and should serve as a valuable complement to methods developed for other types of protein side chains.     

Backbone and side chain dynamics of mutant calmodulin-peptide complexes

The mechanism of long-range coupling of allosteric sites in calcium-saturated calmodulin (CaM) has been explored by characterizing structural and dynamics effects of mutants of calmodulin in complex with a peptide corresponding to the smooth muscle myosin light chain kinase calmodulin-binding domain (smMLCKp). Four CaM mutants were examined: D95N and D58N, located in Ca2+-binding loops; and M124L and E84K, located in the target domain-binding site of CaM. Three of these mutants have altered allosteric coupling either between Ca2+-binding sites (D58N and D95N) or between the target- and Ca2+-binding sites (E84K). The structure and dynamics of the mutant calmodulins in complex with smMLCKp were characterized using solution NMR. Analysis of chemical shift perturbations was employed to detect largely structural perturbations. 15N and 2H relaxation was employed to detect perturbations of the dynamics of the backbone and methyl-bearing side chains of calmodulin. The least median squares method was found to be robust in the detection of perturbed sites. The main chain dynamics of calmodulin are found to be largely unresponsive to the mutations. Three mutants show significantly perturbed dynamics of methyl-bearing side chains. Despite the pseudosymmetric location of Ca2+-binding loop mutations D58N and D95N, the dynamic response of CaM is asymmetric, producing long-range perturbation in D58N and almost none in D95N. The mutations located at the target domain-binding site have quite different effects. For M124L, a local perturbation of the methyl dynamics is observed, while the E84K mutation produces a long-range propagation of dynamic perturbations along the target domain-binding site.     

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.     

The unusual internal motion of the villin headpiece subdomain

The thermostable 36‐residue subdomain of the villin headpiece (HP36) is the smallest known cooperatively folding protein. Although the folding and internal dynamics of HP36 and close variants have been extensively studied, there has not been a comprehensive investigation of side‐chain motion in this protein. Here, the fast motion of methyl‐bearing amino acid side chains is explored over a range of temperatures using site‐resolved solution nuclear magnetic resonance deuterium relaxation. The squared generalized order parameters of methyl groups extensively spatially segregate according to motional classes. This has not been observed before in any protein studied using this methodology. The class segregation is preserved from 275 to 305 K. Motions detected in Helix 3 suggest a fast timescale of conformational heterogeneity that has not been previously observed but is consistent with a range of folding and dynamics studies. Finally, a comparison between the order parameters in solution with previous results based on solid‐state nuclear magnetic resonance deuterium line shape analysis of HP36 in partially hydrated powders shows a clear disagreement for half of the sites. This result has significant implications for the interpretation of data derived from a variety of approaches that rely on partially hydrated protein samples.     

Dynamic transition associated with the thermal denaturation of a small Beta protein

We studied the temperature dependence of the picosecond internal dynamics of an all-beta protein, neocarzinostatin, by incoherent quasielastic neutron scattering. Measurements were made between 20 degrees C and 71 degrees C in heavy water solution. At 20 degrees C, only 33% of the nonexchanged hydrogen atoms show detectable dynamics, a number very close to the fraction of protons involved in the side chains of random coil structures, therefore suggesting a rigid structure in which the only detectable diffusive movements are those involving the side chains of random coil structures. At 61.8 degrees C, although the protein structure is still native, slight dynamic changes are detected that could reflect enhanced backbone and beta-sheet side-chain motions at this higher temperature. Conversely, all internal dynamics parameters (amplitude of diffusive motions, fraction of immobile scatterers, mean-squared vibration amplitude) rapidly change during heat-induced unfolding, indicating a major loss of rigidity of the beta-sandwich structure. The number of protons with diffusive motion increases markedly, whereas the volume occupied by the diffusive motion of protons is reduced. At the half-transition temperature (T = 71 degrees C) most of backbone and beta-sheet side-chain hydrogen atoms are involved in picosecond dynamics.     

Proton nuclear magnetic resonance study of an active pentapeptide fragment of ubiquitin

The aqueous solution conformation of Tyr-Asn-Ile-Gln-Lys (UB5) corresponding to positions 59-63 of the polypeptide, ubiquitin, has been investigated by proton NMR. Like the parent protein, UB5 induces nonspecifically both T and B lymphocyte differentiation. The various NH and CH resonances of this pentapeptide have been assigned, and its solution conformation has been probed through a study of chemical shift variations with pH, temperature dependence of amide hydrogen chemical shifts, vicinal NH--C alpha H and C alpha H--C beta H2 coupling constant data, and amide hydrogen-exchange rates. The latter were measured in H2O by using a combination of transfer of solvent saturation and saturation recovery NMR experiments. The data are compatible with the assumption of a highly motile dynamic equilibrium among different conformations for this peptide. The various secondary amide hydrogens remain essentially exposed to the solvent. The temperature-dependence study of the amide hydrogen chemical shifts also did not reveal any strong internal hydrogen bonds. A rotamer population analysis of tyrosine and asparagine side chains suggests that two of the rotomers are predominantly populated for each of these residues. From these results, a picture emerges of the dynamic conformation of UB5 in aqueous solution.     

13C NMR and fluorescence analysis of tryptophan dynamics in wild-type and two single-Trp variants of Escherichia coli thioredoxin.

The rotational motion of tryptophan side chains in oxidized and reduced wild-type (WT) Escherichia coli thioredoxin and in two single-tryptophan variants of E. coli thioredoxin was studied in solution in the temperature range 20-50 degrees C from 13C-NMR relaxation rate measurements at 75.4 and 125.7 MHz and at 20 degrees C from steady-state and time-resolved trp fluorescence anisotropy measurements. Tryptophan enriched with 13C at the delta 1 and epsilon 3 sites of the indole ring was incorporated into WT thioredoxin and into two single-trp mutants, W31F and W28F, in which trp-28 or trp-31 of WT thioredoxin was replaced, respectively, with phenylalanine. The NMR relaxation data were interpreted using the Lipari and Szabo "model-free" approach (G. Lipari and A. Szabo. 1982. J. Amer. Chem. Soc. 104:4546-4559) with trp steady-state anisotropy data included for the variants at 20 degrees C. Values for the correlation time for the overall rotational motion (tau m) from NMR of oxidized and reduced WT thioredoxin at 35 degrees C agree well with those given by Stone et al. (Stone, M. J., K. Chandrasekhar, A. Holmgren, P. E. Wright, and H. J. Dyson. 1993. Biochemistry. 32:426-435) from 15N NMR relaxation rates, and the dependence of tau m on viscosity and temperature was in accord with the Stokes-Einstein relationship. Order parameters (S2) near 1 were obtained for the trp side chains in the WT proteins even at 50 degrees C. A slight increase in the amplitude of motion (decrease in S2) of trp-31, which is near the protein surface, but not of trp-28, which is partially buried in the protein matrix, was observed in reduced relative to oxidized WT thioredoxin. For trp-28 in W31F, order parameters near 1 (S2 > or = 0.8) at 20 degrees C were found, whereas trp-31 in W28F yielded the smallest order parameters (S2 approximately 0.6) of any of the cases. Analysis of time-resolved anisotropy decays in W28F and W31F yielded S2 values in good agreement with NMR, but gave tau m values about 60% smaller. Generally, values of tau e, the effective correlation time for the internal motion, were < or = 60 ps from NMR, whereas somewhat longer times were obtained from fluorescence. The ability of NMR and fluorescence techniques to detect subnanosecond motions in proteins reliably is examined.     

Methodological advancements for characterising protein side chains by NMR spectroscopy

The surface of proteins is covered by side chains of polar amino acids that are imperative for modulating protein functionality through the formation of noncovalent intermolecular interactions. However, despite their tremendous importance, the unique structures of protein side chains require tailored approaches for investigation by nuclear magnetic resonance spectroscopy and so have traditionally been understudied compared with the protein backbone. Here, we review substantial recent methodological advancements within nuclear magnetic resonance spectroscopy to address this issue. Specifically, we consider advancements that provide new insight into methyl-bearing side chains, show the potential of using non-natural amino acids and reveal the actions of charged side chains. Combined, the new methods promise unprecedented characterisations of side chains that will further elucidate protein function.     

A paramagnetic probe to localize residues next to carboxylates on protein surfaces

It is shown that the paramagnetic properties of lanthanides can be exploited to obtain information on specific parts of a protein surface. Owing to the high affinity of coordinatively unsaturated lanthanide complexes for oxygen donors, carboxylate groups can be used as preferential targets for the interaction. The DO3A ligand is particularly useful in these studies, as it coordinates lanthanides in a heptadentate fashion, leaving two sites available for exogenous donors. A solution of a (15)N-labeled sample protein, calbindin D(9k) (75 residues), was titrated with up to 200% of Gd(III)-DO3A complex, and an inversion recovery (15)N-(1)H HSQC experiment was used to measure the paramagnetic contributions to the longitudinal relaxation rates of the amide protons. Relaxation data were used as distance constraints to estimate the number of interacting complexes and the occupancies of their binding sites. Four preferential interaction sites on the protein surface are found. Inspection of the various carboxylate side chains on the surface of the protein indicates that Gd(III)-DO3A interacts preferentially with carboxylate-rich regions, rather than with isolated carboxylates, suggesting the possibility of chelation of one Gd(III)-DO3A molecule by two carboxylate groups. Gd(III)-DO3A is thus a valuable semi-selective probe for clusters of negative charges on the protein surface.     

13C relaxation experiments for aromatic side chains employing longitudinal- and transverse-relaxation optimized NMR spectroscopy

Aromatic side chains are prevalent in protein binding sites, perform functional roles in enzymatic catalysis, and form an integral part of the hydrophobic core of proteins. Thus, it is of great interest to probe the conformational dynamics of aromatic side chains and its response to biologically relevant events. Indeed, measurements of (13)C relaxation rates in aromatic moieties have a long history in biomolecular NMR, primarily in the context of samples without isotope enrichment that avoid complications due to the strong coupling between neighboring (13)C spins present in uniformly enriched proteins. Recently established protocols for specific (13)C labeling of aromatic side chains enable measurement of (13)C relaxation that can be analyzed in a straightforward manner. Here we present longitudinal- and transverse-relaxation optimized pulse sequences for measuring R (1), R (2), and {(1)H}-(13)C NOE in specifically (13)C-labeled aromatic side chains. The optimized R (1) and R (2) experiments offer an increase in sensitivity of up to 35 % for medium-sized proteins, and increasingly greater gains are expected with increasing molecular weight and higher static magnetic field strengths. Our results highlight the importance of controlling the magnetizations of water and aliphatic protons during the relaxation period in order to obtain accurate relaxation rate measurements and achieve full sensitivity enhancement. We further demonstrate that potential complications due to residual two-bond (13)C-(13)C scalar couplings or dipolar interactions with neighboring (1)H spins do not significantly affect the experiments. The approach presented here should serve as a valuable complement to methods developed for other types of protein side chains.     

Investigations of the supramolecular structure of individual diphenylalanine nano- and microtubes by polarized Raman microspectroscopy

Polarized Raman microspectroscopy and atomic force microscopy were used to obtain quantitative information regarding the molecular structure of individual diphenylalanine (FF) nano- and microtubes. The frequencies of the Raman spectral bands corresponding to the amide I (1690 cm(-1)) and amide III (1249 cm(-1)) indicated that the FF-molecules interact by hydrogen bonding at the N-H and not at the C═O sites. The calculated mean orientation angles of the principal axes of the Raman tensors (PARTs) obtained from the polarized Raman spectral measurements were 41 ± 4° for the amide I and 59 ± 5° for amide III. On the basis of the orientation of the PART for the amide I mode, it was found that the C═O bond is oriented at an angle of 8 ± 4° to the tube axis. These values did not vary significantly with the diameter of the tubes (range 400-1700 nm) and were in agreement with the molecular structure proposed previously for larger crystalline specimens.