NMR (13)C-relaxation study of base and sugar dynamics in GCAA RNA hairpin tetraloop

Intramolecular dynamics of a 14-mer RNA hairpin including GCAA tetraloop was investigated by (13)C NMR relaxation. R(1) and R(1rho) relaxation rates were measured for all protonated base carbons as well as for C1' carbons of ribose sugars at several magnetic field strengths. The data has been interpreted in the framework of modelfree analysis [G. Lipari and A. Szabo. J Am Chem Soc 104, 4546-4559 (1982); G. Lipari and A. Szabo. J Am Chem Soc 104, 4559-4570 (1982)] characterizing the internal dynamics of the molecule by order parameters and correlation times for fast motions on picosecond to nanosecond time scale and by contributions of the chemical exchange. The fast dynamics reveals a rather rigid stem and a significantly more flexible loop. The cytosine and the last adenine bases in the loop as well as all the loop sugars exhibit a significant contribution of conformational equilibrium on microsecond to millisecond time scale. The high R(1rho) values detected on both base and sugar moieties of the loop indicate coordinated motions in this region. A semiquantitative analysis of the conformational equilibrium suggests the exchange rates on the order of 10(4) s(-1). The results are in general agreement with dynamics studies of GAAA loops by NMR relaxation and fluorescent spectroscopy and support the data on the GCAA loop dynamics obtained by MD simulations.     

Backbone dynamics in dihydrofolate reductase complexes: role of loop flexibility in the catalytic mechanism

To elucidate the influence of local motion of the polypeptide chain on the catalytic mechanism of an enzyme, we have measured (15)N relaxation data for Escherichia coli dihydrofolate reductase in three different complexes, representing different stages in the catalytic cycle of the enzyme. NMR relaxation data were analyzed by the model-free approach, corrected for rotational anisotropy, to provide insights into the backbone dynamics. There are significant differences in the backbone dynamics in the different complexes. Complexes in which the cofactor binding site is occluded by the Met20 loop display large amplitude motions on the picosecond/nanosecond time scale for residues in the Met20 loop, the adjacent betaF-betaG loop and for residues 67-69 in the adenosine binding loop. Formation of the closed Met20 loop conformation in the ternary complex with folate and NADP(+), results in attenuation of the motions in the Met20 loop and the betaF-betaG loop but leads to increased flexibility in the adenosine binding loop. New fluctuations on a microsecond/millisecond time scale are observed in the closed E:folate:NADP(+) complex in regions that form hydrogen bonds between the Met20 and the betaF-betaG loops. The data provide insights into the changes in backbone dynamics during the catalytic cycle and point to an important role of the Met20 and betaF-betaG loops in controlling access to the active site. The high flexibility of these loops in the occluded conformation is expected to promote tetrahydrofolate-assisted product release and facilitate binding of the nicotinamide ring to form the Michaelis complex. The backbone fluctuations in the Met20 loop become attenuated once it closes over the active site, thereby stabilizing the nicotinamide ring in a geometry conducive to hydride transfer. Finally, the relaxation data provide evidence for long-range motional coupling between the adenosine binding loop and distant regions of the protein.     

Tyr74 is essential for the formation, stability and function of Plasmodium falciparum triosephosphate isomerase dimer

Plasmodium falciparum triosephosphate isomerase (PfTIM) is known to be functional only as a homodimer. Although many studies have shown that the interface Cys13 plays a major role in the stability of the dimer, a few reports have demonstrated that structurally conserved Tyr74 may be essential for the stability of PfTIM dimer. To understand the role of Tyr74, we have performed molecular dynamics (MD) simulations of monomeric and dimeric PfTIM mutated to glycine and cysteine at position 74. Simulations of the monomer revealed that mutant Tyr74Gly does not produce changes in folding and stability of the monomer. Interestingly, comparison of the flexibility of Tyr74 in the monomer and dimer revealed that this residue possesses an intrinsic restricted mobility, indicating that Tyr74 is an anchor residue required for homodimerization. Tyr74 also appears to play an important role in binding by facilitating the disorder-to-order transitions of loops 1 and 3, which allows Cys13 to form favorable interactions with loop 3 and Lys12 to be locked in a favorable position for catalysis. High-temperature MD simulations of the wild-type and Tyr74Gly PfTIM dimers showed that the aromatic moiety of Tyr74 is necessary to preserve the geometry and native contacts between loops 1 and 3 at the interface of the dimer. Disulfide cross-linking between mutant Tyr74Cys and Cys13 further revealed that Tyr74 stabilizes the geometry of loop 1 (which contains the catalytic residue Lys12) and the interactions between loops 1 and 3 via aromatic-aromatic interactions with residues Phe69, Tyr101, and Phe102. Principal component analysis showed that Tyr74 is also necessary to preserve the collective motions in the dimer that contribute to the catalytic efficiency of PfTIM dimer. We conclude that Tyr74 not only plays a role in the stability of the dimer, but also participates in the dimerization process and collective motions via coupled disorder-to-order transitions of intrinsically disordered regions, necessary for efficiency in the catalytic function of PfTIM.     

Modelling the dynamics of an antigenic peptide using NMR relaxation data.

The internal dynamics of a cyclic peptide which was designed to mimic an antigenic loop of the haemagglutinin, is studied through heteronuclear relaxation along the 13C alpha-1H alpha vectors and through homonuclear relaxation along the 1H alpha-1HN and 1H beta-1H beta' vectors. Order parameters are extracted from the longitudinal and cross-relaxation data. Molecular dynamics simulations are performed and the order parameters are calculated in different ways from the trajectories. The simulation, which is performed in vacuo, gives smaller order parameters (vector motions of larger amplitude) than the experimental results. However, the general features of the experimental order parameters are reproduced by the molecular dynamics simulation. The flexibility of the molecule can then be investigated from the results of the molecular dynamics. It shows that the mobility observed through the order parameters is due to motions in flanking regions, remote from the observed vectors.     

Dynamics of the active site loops in catalyzing aminoacylation reaction in seryl and histidyl tRNA synthetases

Aminoacylation reaction is the first step of protein biosynthesis. The catalytic reorganization at the active site of aminoacyl tRNA synthetases (aaRSs) is driven by the loop motions. There remain lacunae of understanding concerning the catalytic loop dynamics in aaRSs. We analyzed the functional loop dynamics in seryl tRNA synthetase from Methanopyrus kandleri (^mkSerRS) and histidyl tRNA synthetases from Thermus thermophilus (^ttHisRS), respectively, using molecular dynamics. Results confirm that the motif 2 loop and other active site loops are flexible spots within the catalytic domain. Catalytic residues of the loops form a network of interaction with the substrates to form a reactive state. The loops undergo transitions between closed state and open state and the relaxation of the constituent residues occurs in femtosecond to nanosecond time scale. Order parameters are higher for constituent catalytic residues which form a specific network of interaction with the substrates to form a reactive state compared to the Gly residues within the loop. The development of interaction is supported from mutation studies where the catalytic domain with mutated loop exhibits unfavorable binding energy with the substrates. During the open-close motion of the loops, the catalytic residues make relaxation by ultrafast librational motion as well as fast diffusive motion and subsequently relax rather slowly via slower diffusive motion. The Gly residues act as a hinge to facilitate the loop closing and opening by their faster relaxation behavior. The role of bound water is analyzed by comparing implicit solvent-based and explicit solvent-based simulations. Loops fail to form catalytically competent geometry in absence of water. The present result, for the first time reveals the nature of the active site loop dynamics in aaRS and their influence on catalysis.     

Simultaneous processing of solid-state NMR relaxation and 1D-MAS exchange data: the backbone dynamics of free vs. binase-bound barstar

Two types of dynamic solid-state NMR experiments-relaxation and 1D-MAS exchange-were combined for the investigation of the backbone dynamics of a 15% randomly 15N-enriched protein barstar in both free and binase-bound states. The main novelty of this work is a simultaneous quantitative processing of the results of these two types of experiments that we call Simultaneous Relaxation and Exchange Data Analysis (SREDA) approach. It extends the well-known model-free approach such that it permits to discriminate between various motional models (jumps between different sites, wobbling in a cone, etc.). This objective cannot be achieved by analyzing the relaxation or exchange data separately. The SREDA approach was applied to probe a modification of the average backbone dynamics of barstar upon forming a complex with another protein binase. T(1) and off-resonance T(1rho) relaxation times of 15N backbone nuclei were measured at three temperatures between 0 and 45 degrees C, 1D-MAS exchange (CODEX) data were obtained at room temperature within the mixing time range from 0.3 to 200 ms. It has been found that the barstar backbone participates in two molecular processes with correlation times in the 10(-9)-10(-7) and 10(-3)-10(-2) s ranges. Forming the complex with binase results in a significant decrease of the amplitudes of both motions, suggesting that the complex is a more rigid and stable structure than free barstar.     

A molecular dynamics analysis of protein structural elements

The relation between protein secondary structure and internal motions was examined by using molecular dynamics to calculate positional fluctuations of individual helix, beta-sheet, and loop structural elements in free and substrate-bound hen egg-white lysozyme. The time development of the fluctuations revealed a general correspondence between structure and dynamics; the fluctuations of the helices and beta-sheets converged within the 101 psec period of the simulation and were lower than average in magnitude, while the fluctuations of the loop regions were not converged and were mostly larger than average in magnitude. Notable exceptions to this pattern occurred in the substrate-bound simulation. A loop region (residues 101-107) of the active site cleft had significantly reduced motion due to interactions with the substrate. Moreover, part of a loop and a 3(10) helix (residues of 67-88) not in contact with the substrate showed a marked increase in fluctuations. That these differences in dynamics of free and substrate-bound lysozyme did not result simply from sampling errors was established by an analysis of the variations in the fluctuations of the two halves of the 101 psec simulation of free lysozyme. Concerted transitions of four to five mainchain phi and psi angles between dihedral wells were shown to be responsible for large coordinate shifts in the loops. These transitions displaced six or fewer residues and took place either abruptly, in 1 psec or less, or with a diffusive character over 5-10 psec. Displacements of rigid secondary structures involved longer timescale motions in bound lysozyme; a 0.5 A rms change in the position of a helix occurred over the 55 psec simulation period. This helix reorientation within the protein appears to be a response to substrate binding. There was little correlation between the solvent accessible surface area and the dynamics of the different structural elements.     

Global and local motions in ribonuclease A: a molecular dynamics study

The understanding of protein dynamics is one of the major goals of structural biology. A direct link between protein dynamics and function has been provided by x-ray studies performed on ribonuclease A (RNase A) (B. F. Rasmussen et al., Nature, 1992, Vol. 357, pp. 423-424; L. Vitagliano et al., Proteins: Structure, Function, and Genetics, 2002, Vol. 46, pp. 97-104). Here we report a 3 ns molecular dynamics simulation of RNase A in water aimed at characterizing the dynamical behavior of the enzyme. The analysis of local and global motions provides interesting insight on the dynamics/function relationship of RNase A. In agreement with previous crystallographic reports, the present study confirms that the RNase A active site is constituted by rigid (His12, Asn44, Thr45) and flexible (Lys41, Asp83, His119, Asp121) residues. The analysis of the global motions, performed using essential dynamics, shows that the two beta-sheet regions of RNase A move coherently in opposite directions, thus modifying solvent accessibility of the active site, and that the mixed alpha/3(10)-helix (residues 50-60) behaves as a mechanical hinge during the breathing motion of the protein. These data demonstrate that this motion, essential for RNase A substrate binding and release, is an intrinsic dynamical property of the ligand-free enzyme.     

Dimerization affects collective dynamics of triosephosphate isomerase

Molecular dynamics simulations (30-60 ns runs) are performed on free/apo triosephosphate isomerase (TIM) to determine any correlation between collective motions and loop 6 dynamics. Native TIM is reported to be active only as a homodimer even though cooperativity has not been observed between the two identical subunits. Both dimeric and monomeric (isolated from dimer) forms of TIM are simulated in explicit water at 300 K and 1 bar to inspect any differences between the structures in terms of fluctuation dynamics and functionally important loop 6 dynamics/closure. Significant cross-correlations between residue fluctuations are observed in the dimer, which result from the global counter-rotations of the two identical subunits in the essential modes of the dimer. Specifically, the first essential mode contributing to 34% of overall motion of the dimer is strongly coupled to the loop 6's closure over the active site. In contrast, such significant correlations cannot be observed in the monomeric structure, which maintains relatively localized motions of the loops in the essential modes. Thus, the onset of collective motions at ns time scale due to dimerization has functional implications as to the coordination of loop 6 closure.     

Role of loop dynamics in thermal stability of mesophilic and thermophilic adenylosuccinate synthetase: a molecular dynamics and normal mode analysis study

Enzymes from thermophiles are poorly active at temperatures at which their mesophilic homologs exhibit high activity and attain corresponding active states at high temperatures. In this study, comparative molecular dynamics (MD) simulations, supplemented by normal mode analysis, have been performed on an enzyme Adenylosuccinate synthetase (AdSS) from E. coli (mesophilic) and P. horikoshii (thermophilic) systems to understand the effects of loop dynamics on thermal stability of AdSS. In mesophilic AdSS, both ligand binding and catalysis are facilitated through the coordinated movement of five loops on the protein. The simulation results suggest that thermophilic P. horikoshii preserves structure and catalytic function at high temperatures by using the movement of only a subset of loops (two out of five) for ligand binding and catalysis unlike its mesophilic counterpart in E. coli. The pre-arrangement of the catalytic residues in P. horikoshii is well-preserved and salt bridges remain stable at high temperature (363K). The simulations suggest a general mechanism (including pre-arrangement of catalytic residues, increased polar residue content, stable salt bridges, increased rigidity, and fewer loop movements) used by thermophilic enzymes to preserve structure and be catalytically active at elevated temperatures.     

Induction of apoptosis of tumor cells by binase

An induction of apoptosis by RNase from Bacillus intermedius (binase) and its mutants characterized with low catalytic activity (Lys26Ala and His101Glu) in human myelogenic erythroleukemia K562 cells, human lung carcinoma A549 cells and human peripheral blood mononuclear cells was studied. For the first time selective apoptogenic effects of binase toward leukemic blood cells was determined. Neither antiproliferative nor apoptotic effects of binase were detected in normal human peripheral blood mononuclear cells. Formation of low molecular weight oligonucleosomal DNA fragments (less than 50 Kb) which are an early marks of apoptosis was registered in solid tumor cells treated by binase. Using mutant RNases it was shown that decrease of catalytic activity to 2.5% of wild type enzyme activity leads to the loss of apoptogenic properties of enzyme. Selective apoptogenicity of binase found towards malignant cells confirmed that antitumor agents based on bacterial RNases could be considered as an alternative to standard chemotherapeutic drugs.     

Backbone dynamics of apocytochrome b5 in its native, partially folded state

The backbone dynamics in the native state of apocytochrome b5 were studied using 15N nuclear magnetic spin relaxation measurements. The field (11.7 and 14.1 T) and temperature (10-25 degrees C) dependence of the relaxation parameters (R1, R2, and R1rho) and the 1H-15N NOE established that the protein undergoes multiple time scale internal motions related to the secondary structure. The relaxation data were analyzed with the reduced spectral density mapping approach and within the extended model-free framework. The apoprotein was confirmed to contain a disordered heme-binding loop of approximately 30 residues with dynamics on the sub-nanosecond time scale (0.6 < S2 < 0.7, 100 ps < taue < 500 ps). This loop is attached to a structured hydrophobic core, rigid on the picosecond time scale (S2 > 0.75, taue < 50 ps). The inability to fit the data for several residues with the model-free protocol revealed the presence of correlated motion. An exchange contribution was detected in the transverse relaxation rate (R2) of all residues. The differential temperature response of R2 along the backbone supported slower exchange rates for residues in the loop (tauex > 300 micros) than for the folded polypeptide chain (tauex < 150 micros). The distribution of the reduced spectral densities at the 1H and 15N frequencies followed the dynamic trend and predicted the slowing of the internal motions at 10 degrees C. Comparison of the dynamics with those of the holoprotein [Dangi, B., Sarma, S., Yan, C., Banville, D. L., and Guiles, R. D. (1998) Biochemistry 37, 8289-8302] demonstrated that binding of the heme alters the time scale of motions both in the heme-binding loop and in the structured hydrophobic core.     

Global and local molecular dynamics of a bacterial carboxylesterase provide insight into its catalytic mechanism

Carboxylesterases (CEs) are ubiquitous enzymes responsible for the detoxification of xenobiotics. In humans, substrates for these enzymes are far-ranging, and include the street drug heroin and the anticancer agent irinotecan. Hence, their ability to bind and metabolize substrates is of broad interest to biomedical science. In this study, we focused our attention on dynamic motions of a CE from B. subtilis (pnbCE), with emphasis on the question of what individual domains of the enzyme might contribute to its catalytic activity. We used a 10 ns all-atom molecular dynamics simulation, normal mode calculations, and enzyme kinetics to understand catalytic consequences of structural changes within this enzyme. Our results shed light on how molecular motions are coupled with catalysis. During molecular dynamics, we observed a distinct C-C bond rotation between two conformations of Glu310. Such a bond rotation would alternately facilitate and impede protonation of the active site His399 and act as a mechanism by which the enzyme alternates between its active and inactive conformation. Our normal mode results demonstrate that the distinct low-frequency motions of two loops in pnbCE, coil_5 and coil_21, are important in substrate conversion and seal the active site. Mutant CEs lacking these external loops show significantly reduced rates of substrate conversion, suggesting this sealing motion prevents escape of substrate. Overall, the results of our studies give new insight into the structure-function relationship of CEs and have implications for the entire family of α/β fold family of hydrolases, of which this CE is a member.     

Surface dynamics of bacteriorhodopsin as revealed by (13)C NMR studies on [(13)C]Ala-labeled proteins: detection of millisecond or microsecond motions in interhelical loops and C-terminal alpha-helix

We have recorded (13)C NMR spectra of [2-(13)C]-, [1-(13)C]-, [3-(13)C],- and [1,2,3-(13)C(3)]Ala-labeled bacteriorhodopsin (bR), and its mutants, A196G, A160G, and A103C, by means of cross polarization-magic angle spinning (CP-MAS) and dipolar decoupled-magic angle spinning (DD-MAS) techniques, to reveal the conformation and dynamics of bR, with emphasis on the loop and C-terminus structures. The (13)C NMR signals of the loop (C-D, E-F, and F-G) regions were almost completely suppressed from [2-(13)C]-, [1-(13)C]Ala-, and [1-(13)C]Gly-labeled bR, due to the presence of conformational fluctuation with correlation times of 10(-4) s that interfered with the peak-narrowing by magic angle spinning. The observation of such suppressed peaks for specific residues provides a unique means of detecting intermediate frequency motions on the time scale of ms or micros in the surface loops of membrane proteins. Instead, the three well-resolved (13)C CP-MAS NMR signals of [2-(13)C]Ala-bR, at 50.38, 49.90, and 47.96 ppm, were ascribed to the C-terminal alpha-helix previously proposed from the data for [3-(13)C]Ala-bR: the former two peaks were assigned to Ala 232 and 238, in view of the results of successive proteolysis experiments, while the highest-field peak was ascribed to Ala 235 prior to Pro 236. Even such (13)C NMR signals were substantially broadened when (13)C NMR spectra of fully labeled [1,2,3-(13)C]Ala-bR were recorded, because the broadening and splitting of peaks due to the accelerated transverse relaxation rate caused by the increased number of relaxation pathways through a number of (13)C-(13)C homo-nuclear dipolar interactions and scalar J couplings, respectively, are dominant among (13)C-labeled nuclei. In addition, approximate correlation times for local conformational fluctuations of different domains, including the C-terminal tail, C-terminal alpha-helix, loops, and transmembrane alpha-helices, were estimated by measurement of the spin-lattice relaxation times in the laboratory frame and spin-spin relaxation times under the conditions of cross-polarization-magic angle spinning, and comparative study of suppressed specific peaks between the CP-MAS and DD-MAS experiments.     

Active site dynamics in the lead-dependent ribozyme

Conformational dynamics are an important property of ribozymes and other RNA molecules but there is currently only limited information on the relationship between dynamics and RNA function. A recent structural study of the lead-dependent ribozyme, known as the leadzyme, showed significant dynamics at the active site and indicated that a structural rearrangement is required for the reaction to proceed from the ground to the transition state. In this work, microsecond-to-millisecond dynamics of the leadzyme are probed by analysis of the power dependence of (13)C NMR relaxation times in the rotating frame (T(1)(rho)). These results revealed a wide range of conformational dynamics for various residues in the leadzyme. For residue A25 in the active site, the power dependence of T(1)(rho) yielded an exchange lifetime similar to that previously measured by line-shape analysis, and provides an important calibration of this T(1)(rho) methodology for probing the dynamics of macromolecules. Strong evidence was also found for a previously suggested dynamic network of hydrogen bonds stabilizing the GAAA tetraloop motif. Within the active site of the leadzyme, internal motions are observed on a wide variety of time scales, suggesting a complex landscape of accessible states, and potential correlations between observed motions and catalytic function are discussed. These results demonstrate that the power dependence of (13)C T(1)(rho) relaxation times provides a valuable method for probing dynamics in nucleic acids.     

Molecular dynamics studies of caspase-3

Caspase-3 is a fundamental target for pharmaceutical interventions against a variety of diseases involving disregulated apoptosis. The enzyme is active as a dimer with two symmetry-related active sites, each featuring a Cys-His catalytic dyad and a selectivity loop, which recognizes the characteristic DEVD pattern of the substrate. Here, a molecular dynamics study of the enzyme in complex with two pentapeptide substrates DEVDG is presented, which provides a characterization of the dynamic properties of the active form in aqueous solution. The mobility of the substrate and that of the catalytic residues are rather low indicating a distinct preorganization effect of the Michaelis complex. An essential mode analysis permits us to identify coupled motions between the two monomers. In particular, it is found that the motions of the two active site loops are correlated and tend to steer the substrate toward the reactive center, suggesting that dimerization has a distinct effect on the dynamic properties of the active site regions. The selectivity loop of one monomer turns out to be correlated with the N-terminal region of the p12 subunit of the other monomer, an interaction that is also found to play a fundamental role in the electrostatic stabilization of the quaternary structure. To further characterize the specific influence of dimerization on the enzyme essential motions, a molecular dynamics analysis is also performed on the isolated monomer.     

Backbone dynamics of calcium-loaded calbindin D9k studied by two-dimensional proton-detected 15N NMR spectroscopy.

Backbone dynamics of calcium-loaded calbindin D9k have been investigated by two-dimensional proton-detected heteronuclear nuclear magnetic resonance spectroscopy, using a uniformly 15N enriched protein sample. Spin-lattice relaxation rate constants, spin-spin relaxation rate constants, and steady-state [1H]-15N nuclear Overhauser effects were determined for 71 of the 72 backbone amide 15N nuclei. The relaxation parameters were analyzed using a model-free formalism that incorporates the overall rotational correlation time of the molecule, and a generalized order parameter (S2) and an effective internal correlation time for each amide group. Calbindin D9k contains two helix-loop-helix motifs joined by a linker loop at one end of the protein and a beta-type interaction between the two calcium-binding loops at the other end. The amplitude of motions for the calcium-binding loops and the helices are similar, as judged from the average S2 values of 0.83 +/- 0.05 and 0.85 +/- 0.04, respectively. The linker region joining the two calcium-binding subdomains of the molecule has a significantly higher flexibility, as indicated by a substantially lower average S2 value of 0.59 +/- 0.23. For residues in the linker loop and at the C-terminus, the order parameter is further decomposed into separate order parameters for motional processes on two distinct time scales. The effective correlation times are significantly longer for helices I and IV than for helices II and III or for the calcium-binding loops. Residue by residue comparisons reveal correlations of the order parameters with both the crystallographic B-factors and amide proton exchange rates, despite vast differences in the time scales to which these properties are sensitive. The order parameters are also utilized to distinguish regions of the NMR-derived three-dimensional structure of calbindin D9k that are poorly defined due to inherently high flexibility, from poorly defined regions with average flexibility but a low density of structural constraints.     

Backbone dynamics of barstar: a (15)N NMR relaxation study

Backbone dynamics of uniformly (15)N-labeled barstar have been studied at 32 degrees C, pH 6.7, by using (15)N relaxation data obtained from proton-detected 2D (1)H-(15)N NMR spectroscopy. (15)N spin-lattice relaxation rate constants (R(1)), spin-spin relaxation rate constants (R(2)), and steady-state heteronuclear (1)H-(15)N NOEs have been determined for 69 of the 86 (excluding two prolines and the N-terminal residue) backbone amide (15)N at a magnetic field strength of 14.1 Tesla. The primary relaxation data have been analyzed by using the model-free formalism of molecular dynamics, using both isotropic and axially symmetric diffusion of the molecule, to determine the overall rotational correlation time (tau(m)), the generalized order parameter (S(2)), the effective correlation time for internal motions (tau(e)), and NH exchange broadening contributions (R(ex)) for each residue. As per the axially symmetric diffusion, the ratio of diffusion rates about the unique and perpendicular axes (D( parallel)/D( perpendicular)) is 0.82 +/- 0.03. The two results have only marginal differences. The relaxation data have also been used to map reduced spectral densities for the NH vectors of these residues at three frequencies: 0, omega(H), and omega(N), where omega(H),(N) are proton and nitrogen Larmor frequencies. The value of tau(m) obtained from model-free analysis of the relaxation data is 5.2 ns. The reduced spectral density analysis, however, yields a value of 5.7 ns. The tau(m) determined here is different from that calculated previously from time-resolved fluorescence data (4.1 ns). The order parameter ranges from 0.68 to 0.98, with an average value of 0.85 +/- 0.02. A comparison of the order parameters with the X-ray B-factors for the backbone nitrogens of wild-type barstar does not show any considerable correlation. Model-free analysis of the relaxation data for seven residues required the inclusion of an exchange broadening term, the magnitude of which ranges from 2 to 9.1 s(-1), indicating the presence of conformational averaging motions only for a small subset of residues.     

Preparation of spin-labeled bacterial ribonuclease from Bacillus intermedius 7P for EPR studies of protein dynamics.

The optimal conditions for labeling of binase (a bacterial ribonuclease isolated from Bacillus intermedius 7P) with a bromoacetamide spin label have been determined. The label is suitable for probing the dynamics of the protein by electron paramagnetic resonance (EPR). Binase samples specifically labeled at residue His-101 of the active center were prepared by incubation for 48 h at 20 degrees C in potassium phosphate buffer (pH 5.5) containing the bromoacetamide (1:3 protein/label molar ratio). Fluorescence assay showed that the structure of the labeled binase is indistinguishable from that of the native protein.