Theory for Spin−Lattice Relaxation of Spin Probes on Weakly Deformable DNA

Direct simulation of magnetic resonance relaxation rates and line shapes from molecular trajectories

We simulate spin relaxation processes, which may be measured by either continuous wave or pulsed magnetic resonance techniques, using trajectory-based simulation methodologies. The spin-lattice relaxation rates are extracted numerically from the relaxation simulations. The rates obtained from the numerical fitting of the relaxation curves are compared to those obtained by direct simulation from the relaxation Bloch-Wangsness-Abragam-Redfield theory (BWART). We have restricted our study to anisotropic rigid-body rotational processes, and to the chemical shift anisotropy (CSA) and a single spin-spin dipolar (END) coupling mechanisms. Examples using electron paramagnetic resonance (EPR) nitroxide and nuclear magnetic resonance (NMR) deuterium quadrupolar systems are provided. The objective is to compare those rates obtained by numerical simulations with the rates obtained by BWART. There is excellent agreement between the simulated and BWART rates for a Hamiltonian describing a single spin (an electron) interacting with the bath through the chemical shift anisotropy (CSA) mechanism undergoing anisotropic rotational diffusion. In contrast, when the Hamiltonian contains both the chemical shift anisotropy (CSA) and the spin-spin dipolar (END) mechanisms, the decay rate of a single exponential fit of the simulated spin-lattice relaxation rate is up to a factor of 0.2 smaller than that predicted by BWART. When the relaxation curves are fit to a double exponential, the slow and fast rates extracted from the decay curves bound the BWART prediction. An extended BWART theory, in the literature, includes the need for multiple relaxation rates and indicates that the multiexponential decay is due to the combined effects of direct and cross-relaxation mechanisms.

Conformation and dynamics of nucleotides in bulges and symmetric internal loops in duplex DNA studied by EPR and fluorescence spectroscopies

The dynamics and conformation of base bulges and internal loops in duplex DNA were studied using the bifunctional spectroscopic probe Ç, which becomes fluorescent (Ç(f)) upon reduction of the nitroxide functional group, along with EPR and fluorescence spectroscopies. A one-base bulge was in a conformational equilibrium between looped-out and stacked states, the former favored at higher temperature and the latter at lower temperature. Stacking of bulge bases was favored in two- and three-base bulges, independent of temperature, resulting in DNA bending as evidenced by increased fluorescence of Ç(f). EPR spectra of Ç-labeled three-, four- and five-base symmetrical interior DNA bulges at 20 °C showed low mobility, indicating that the spin-label was stacked within the loop. The spin-label mobility at 37 °C increased as the loops became larger. A considerable variation in fluorescence between different loops was observed, as well as a temperature-dependence within constructs. Fluorescence unexpectedly increased as the size of the loop decreased at 2 °C. Fluorescence of the smallest loops, where a single T·T mismatch was located between the stem region and the probe, was even larger than for the single strand, indicating a considerable local structural deformation of these loops from regular B-DNA. These results show the value of combining EPR and fluorescence spectroscopy to study non-helical regions of nucleic acids.

Conformational equilibria of bulged sites in duplex DNA studied by EPR spectroscopy

Conformational flexibility in nucleic acids provides a basis for complex structures, binding, and signaling. One-base bulges directly neighboring single-base mismatches in nucleic acids can be present in a minimum of two distinct conformations, complicating the examination of the thermodynamics by calorimetry or UV-monitored melting techniques. To provide additional information about such structures, we demonstrate how electron paramagnetic resonance (EPR) active spin-labeled base analogues, base-specifically incorporated into the DNA, are monitors of the superposition of different bulge-mismatch conformations. EPR spectra provide information about the dynamic environments of the probe. This information is cast in terms of "dynamic signatures" that have an underlying basis in structural variations. By examining the changes in the equilibrium of the different states across a range of temperatures, the enthalpy and entropy of the interconversion among possible conformations can be determined. The DNA constructs with a single bulge neighboring a single-base mismatch ("bulge-mismatches") may be approximately modeled as an equilibrium between two possible conformations. This structural information provides insight into the local composition of the bulge-mismatch sequences. Experiments on the bulge-mismatches show that basepairing across the helix can be understood in terms of purine and pyrimidine interactions, rather than specific bases. Measurements of the enthalpy and entropy of formation for the bulge-mismatches by differential scanning calorimetry and UV-monitored melting confirm that the formation of bulge-mismatches is in fact more complicated than a simple two-state process, consistent with the base-specific spectral data that bulge-mismatches exist in multiple conformations in the premelting temperature region. We find that the calculations with the nearest-neighbor (NN) model for the two likely conformations do not correlate well with the populations of structures and thermodynamic parameters inferred from the base-specific EPR dynamics probe. We report that the base-specific spin probes are able to identify a bistable, temperature dependent, switching between conformations for a particular complex bulged construct.

Probing the dynamics of the P1 helix within the Tetrahymena group I intron

RNA conformational transformations are integral to RNA's biological functions. Further, structured RNA molecules exist as a series of dynamic intermediates in the course of folding or complexation with proteins. Thus, an understanding of RNA folding and function will require deep and incisive understanding of its dynamic behavior. However, existing tools to investigate RNA dynamics are limited. Here, we introduce a powerful fluorescence polarization anisotropy approach that utilizes a rare base analogue that retains substantial fluorescence when incorporated into helices. We show that 6-methylisoxanthopterin (6-MI) can be used to follow the nanosecond dynamics of individual helices. We then use 6-MI to probe the dynamics of an individual helix, referred to as P1, within the 400nt Tetrahymena group I ribozyme. Comparisons of the dynamics of the P1 helix in wild type and mutant ribozymes and in model constructs reveal a highly immobilized docked state of the P1 helix, as expected, and a relatively mobile "open complex" or undocked state. This latter result rules out a model in which slow docking of the P1 helix into its cognate tertiary interactions arises from a stable alternatively docked conformer. The results are consistent with a model in which stacking and tertiary interactions of the A(3) tether connecting the P1 helix to the body of the ribozyme limit P1 mobility and slow its docking, and this model is supported by cross-linking results. The ability to isolate the nanosecond motions of individual helices within complex RNAs and RNA/protein complexes will be valuable in distinguishing between functional models and in discerning the fundamental behavior of important biological species.

Preparation, resonance assignment, and preliminary dynamics characterization of residue specific 13C/15N-labeled elongated DNA for the study of sequence-directed dynamics by NMR

DNA is a highly flexible molecule that undergoes functionally important structural transitions in response to external cellular stimuli. Atomic level spin relaxation NMR studies of DNA dynamics have been limited to short duplexes in which sensitivity to biologically relevant fluctuations occurring at nanosecond timescales is often inadequate. Here, we introduce a method for preparing residue-specific (13)C/(15)N-labeled elongated DNA along with a strategy for establishing resonance assignments and apply the approach to probe fast inter-helical bending motions induced by an adenine tract. Preliminary results suggest the presence of elevated A-tract independent end-fraying internal motions occurring at nanosecond timescales, which evade detection in short DNA constructs and that penetrate deep (7 bp) within the DNA helix and gradually fade away towards the helix interior.