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

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

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

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

Single-stranded DNA binding of the cold-shock protein CspB from Bacillus subtilis: NMR mapping and mutational characterization

Cold-shock proteins (CSPs) bind to single-stranded nucleic acids, thereby acting as a "RNA chaperone." To gain deeper insights into the rather unspecific nature of ssDNA/RNA binding, we characterized the binding interface of CspB from Bacillus subtilis to a 25-mer of ssDNA (Y-Box25) using heteronuclear 2D NMR spectroscopy. Seventeen residues, including eight out of nine aromatic amino acids, are directly involved in the Y-Box25 interaction and were identified by extreme line broadening of their cross-peaks. Eight residues belong to the earlier proposed RNP binding motifs. A second set of seven backbone amides becomes evident by major chemical shift perturbations reporting remote conformational rearrangements upon binding. These residues are located in loop beta3-beta4 and loopbeta4-beta5, and include Ile18. The individual contributions of the so-identified residues were examined by fluorescence titration experiments of 15 CspB variants. Phenylalanine substitutions in- and outside the RNP motifs significantly reduce the binding affinity. Unrestricted possible backbone conformations of loop beta3-beta4 also markedly contribute to binding. Stopped-flow fluorescence kinetics revealed that the different binding affinities of CspB variants are determined by the dissociation rate, whereas the association rate remains unchanged. This might be of importance for the "RNA chaperone" activity of CspB.     

Identification of Slow Motions in the Reduced Recombinant High-potential Iron Sulfur Protein I (HiPIP I) from Ectothiorhodospira Halophila via 15N Rotating-frame NMR Relaxation Measurements

Rotating-frame 15N relaxation rate (R1) NMR experiments have been performed in order to study the dynamic behavior of the reduced recombinant high-potential iron-sulfur protein iso I (HiPIP I) from Ectothiorhodospira halophila, in the s to ms time range. Measurements of R1 were performed as a function of the effective spin-lock magnetic field amplitude by using both on and off-resonance radio frequency irradiation. The two data sets provided consistent results and were fit globally in order to identify possible exchange processes in an external loop of the reduced HiPIP I. The loop consists of residues 43-45 and the correlation time of the exchange process was determined to be 50 ± 8 s for the backbone nitrogen of Gln 44.     

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

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

Microsecond time scale dynamics in the RXR DNA-binding domain from a combination of spin-echo and off-resonance rotating frame relaxation measurements

Slow protein dynamics can be studied by 15N spin-echo (CPMG) and off-resonance rotating frame relaxation through the effective field dependence of the exchange-mediated relaxation contribution. It is shown that, by a combination of these complementary techniques, a more extended sampling of the microsecond time scale processes is achieved than by either method alone. 15N R2 and improved off-resonance R1 experiments [Mulder et al. (1998) J. Magn. Reson., 131, 351–357] were applied to the 9- cis-retinoic acid receptor DNA-binding domain and allowed the identification of 14 residues exhibiting microsecond time scale dynamics. Assuming exchange between two conformational substates, average lifetimes ranging from 37 to 416 s, and chemical shift differences of up to 3 ppm were obtained. The largest perturbation of tertiary structure was observed for the second zinc finger region, which was found to be disordered in the solution structure [Holmbeck et al. (1998) J. Mol. Biol., 281, 271–284]. Since this zinc-coordinating domain comprises the principal dimerization interface for RXR in a wide repertoire of complexes with different hormone receptors to their cognate response elements, this finding has important implications for our understanding of nuclear receptor assembly on DNA direct repeats. The flexibility observed for the dimerization domain may explain how RXR, through the ability to adaptively interact with a wide variety of highly homologous partner molecules, demonstrates such a versatile DNA-binding repertoire.