Experimental and Theoretical Determination of Nucleic Acid Magnetic Susceptibility: Importance for the Study of Dynamics by Field-Induced Residual Dipolar Couplings

Magnetic control of self-assembly and disassembly in organic materials

Because organic molecules and materials are generally insensitive or weakly sensitive to magnetic fields, a certain means to enhance their magnetic responsiveness needs to be exploited. Here we show a strategy to amplify the magnetic responsiveness of self-assembled peptide nanostructures by synergistically combining the concepts of perfect α-helix and rod-coil supramolecular building blocks. Firstly, we develop a monomeric, nonpolar, and perfect α-helix (MNP-helix). Then, we employ the MNP-helix as the rod block of rod-coil amphiphiles (rod-coils) because rod-coils are well-suited for fabricating responsive assemblies. We show that the self-assembly processes of the designed rod-coils and disassembly of rod-coil/DNA complexes can be controlled in a magnetically responsive manner using the relatively weak magnetic field provided by the ordinary neodymium magnet [0.07 ~ 0.25 Tesla (T)]. These results demonstrate that magnetically responsive organic assemblies usable under practical conditions can be realized by using rod-coil supramolecular building blocks containing constructively organized diamagnetic moieties.

HERMES - A Software Tool for the Prediction and Analysis of Magnetic-Field-Induced Residual Dipolar Couplings in Nucleic Acids

Field-Induced Residual Dipolar Couplings (fiRDC) are a valuable source of long-range information on structure of nucleic acids (NA) in solution. A web application (HERMES) was developed for structure-based prediction and analysis of the (fiRDCs) in NA. fiRDC prediction is based on input 3D model structure(s) of NA and a built-in library of nucleobase-specific magnetic susceptibility tensors and reference geometries. HERMES allows three basic applications: (i) the prediction of fiRDCs for a given structural model of NAs, (ii) the validation of experimental or modeled NA structures using experimentally derived fiRDCs, and (iii) assessment of the oligomeric state of the NA fragment and/or the identification of a molecular NA model that is consistent with experimentally derived fiRDC data. Additionally, the program's built-in routine for rigid body modeling allows the evaluation of relative orientation of domains within NA that is in agreement with experimental fiRDCs.

Impact of nucleic acid self-alignment in a strong magnetic field on the interpretation of indirect spin-spin interactions

Heteronuclear and homonuclear direct (D) and indirect (J) spin-spin interactions are important sources of structural information about nucleic acids (NAs). The Hamiltonians for the D and J interactions have the same functional form; thus, the experimentally measured apparent spin-spin coupling constant corresponds to a sum of J and D. In biomolecular NMR studies, it is commonly presumed that the dipolar contributions to Js are effectively canceled due to random molecular tumbling. However, in strong magnetic fields, such as those employed for NMR analysis, the tumbling of NA fragments is anisotropic because the inherent magnetic susceptibility of NAs causes an interaction with the external magnetic field. This motional anisotropy is responsible for non-zero D contributions to Js. Here, we calculated the field-induced D contributions to 33 structurally relevant scalar coupling constants as a function of magnetic field strength, temperature and NA fragment size. We identified two classes of Js, namely (1)JCH and (3)JHH couplings, whose quantitative interpretation is notably biased by NA motional anisotropy. For these couplings, the magnetic field-induced dipolar contributions were found to exceed the typical experimental error in J-coupling determinations by a factor of two or more and to produce considerable over- or under-estimations of the J coupling-related torsion angles, especially at magnetic field strengths >12 T and for NA fragments longer than 12 bp. We show that if the non-zero D contributions to J are not properly accounted for, they might cause structural artifacts/bias in NA studies that use solution NMR spectroscopy.

Base pair openings and temperature dependence of DNA flexibility

The relationship of base pair openings to DNA flexibility is examined. Published experimental data on the temperature dependence of the persistence length by two different groups are well described in terms of an inhomogeneous Kratky-Porot model with soft and hard joints, corresponding to open and closed base pairs, and sequence-dependent statistical information about the state of each pair provided by a Peyrard-Bishop-Dauxois (PBD) model calculation with no freely adjustable parameters.

Extending the NMR spatial resolution limit for RNA by motional couplings

Experimental resolution of distinct dynamical processes in molecules can prove impossible when they are correlated to one another. In nuclear magnetic resonance (NMR) spectroscopy, couplings between internal and overall motions lead to intractable complexity, depriving insights into functionally important motions. Here we demonstrate that motional couplings can be used to anchor NMR frames of reference onto different parts of an RNA molecule, thus extending the spatial resolution limit for dynamical characterization.

Comparison of alignment tensors generated for native tRNA(Val) using magnetic fields and liquid crystalline media

Residual dipolar couplings (RDCs) complement standard NOE distance and J-coupling torsion angle data to improve the local and global structure of biomolecules in solution. One powerful application of RDCs is for domain orientation studies, which are especially valuable for structural studies of nucleic acids, where the local structure of a double helix is readily modeled and the orientations of the helical domains can then be determined from RDC data. However, RDCs obtained from only one alignment media generally result in degenerate solutions for the orientation of multiple domains. In protein systems, different alignment media are typically used to eliminate this orientational degeneracy, where the combination of RDCs from two (or more) independent alignment tensors can be used to overcome this degeneracy. It is demonstrated here for native E. coli tRNA(Val) that many of the commonly used liquid crystalline alignment media result in very similar alignment tensors, which do not eliminate the 4-fold degeneracy for orienting the two helical domains in tRNA. The intrinsic magnetic susceptibility anisotropy (MSA) of the nucleobases in tRNA(Val) was also used to obtain RDCs for magnetic alignment at 800 and 900 MHz. While these RDCs yield a different alignment tensor, the specific orientation of this tensor combined with the high rhombicity for the tensors in the liquid crystalline media only eliminates two of the four degenerate orientations for tRNA(Val). Simulations are used to show that, in optimal cases, the combination of RDCs obtained from liquid crystalline medium and MSA-induced alignment can be used to obtain a unique orientation for the two helical domains in tRNA(Val).

NMR methods for studying quadruplex nucleic acids

Solution NMR spectroscopy has traditionally played a central role in examining quadruplex structure, dynamics, and interactions. Here, an overview is given of the methods currently applied to structural, dynamics, thermodynamics, and kinetics studies of nucleic acid quadruplexes and associated cations.

Magnetic field induced residual dipolar couplings of imino groups in nucleic acids from measurements at a single magnetic field

For base-paired nucleic acids, variations in (1) J (NH) and the imino (1)H chemical shift are both dominated by hydrogen bond length. In the absence of molecular alignment, the (1) J (NH) coupling for the imino proton then can be approximated by (1) J (NH) = (1.21Hz/ppm)delta(H) - 103.5 +/- 0.6 Hz, where delta(H) represents the chemical shift of the imino proton in ppm. This relation permits imino residual dipolar couplings (RDCs) resulting from magnetic susceptibility anisotropy (MSA) to be extracted from measurement of ((1) J (NH) + RDC) splittings at a single magnetic field strength. Magnetic field-induced RDCs were measured for tRNA(Val) and the alignment tensor determined from magnetic-field alignment of tRNA(Val) agrees well with the tensor calculated by summation of the MSA tensors of the individual nucleobases.

NMR studies of RNA dynamics and structural plasticity using NMR residual dipolar couplings

An increasing number of RNAs are being discovered that perform their functions by undergoing large changes in conformation in response to a variety of cellular signals, including recognition of proteins and small molecular targets, changes in temperature, and RNA synthesis itself. The measurement of NMR residual dipolar couplings (RDCs) in partially aligned systems is providing new insights into the structural plasticity of RNA through combined characterization of large-amplitude collective helix motions and local flexibility in noncanonical regions over a wide window of biologically relevant timescales (<milliseconds). Here, we review RDC methodology for studying RNA structural dynamics and survey what has been learnt thus far from application of these methods. Future methodological challenges are also identified.

Chemical shift tensors of protonated base carbons in helical RNA and DNA from NMR relaxation and liquid crystal measurements

Knowledge of (13)C chemical shift anisotropy (CSA) tensors in nucleotide bases is important for interpretation of NMR relaxation data in terms of local dynamic properties of nucleic acids and for analysis of residual chemical shift anisotropy (RCSA) resulting from weak alignment. CSA tensors for protonated nucleic acid base carbons have been derived from measurements on a uniformly (13)C-enriched helical A-form RNA segment and a helical B-form DNA dodecamer at natural (13)C abundance. The magnitudes of the derived CSA principal values are tightly restricted by the magnetic field dependencies of the (13)C transverse relaxation rates, whereas the tensor orientation and asymmetry follow from quantitative measurements of interference between (13)C-{(1)H} dipolar and (13)C CSA relaxation mechanisms. Changes in the chemical shift between the isotropic and aligned states, Deltadelta, complement these measurements and permit cross-validation. The CSA tensors are determined from the experimental Deltadelta values and relaxation rates, under the assumption that the CSA tensor of any specific carbon in a given type of base is independent of the base position in either the RNA or DNA helix. However, the experimental data indicate that for pyrimidine C(6) carbons in A-form RNA the CSA magnitude is considerably larger than in B-form DNA. This result is supported by quantum chemical calculations and is attributed in part to the close proximity between intranucleotide C(6)H and O(5)' atoms in RNA. The magnitudes of the measured CSA tensors, on average, agree better with previous solid-state NMR results obtained on powdered nucleosides than with prior results from quantum chemical calculations on isolated bases, which depend rather strongly on the level of theory at which the calculations are carried out. In contrast, previously computed orientations of the chemical shift tensors agree well with the present experimental results and exhibit less dependence on the level of theory at which the computations are performed.

Impact of static and dynamic A-form heterogeneity on the determination of RNA global structural dynamics using NMR residual dipolar couplings

We examined how static and dynamic deviations from the idealized A-form helix propagate into errors in the principal order tensor parameters determined using residual dipolar couplings (rdcs). A 20-ns molecular dynamics (MD) simulation of the HIV-1 transactivation response element (TAR) RNA together with a survey of spin relaxation studies of RNA dynamics reveals that pico-to-nanosecond local motions in non-terminal Watson-Crick base-pairs will uniformly attenuate base and sugar one bond rdcs by approximately 7%. Gaussian distributions were generated for base and sugar torsion angles through statistical comparison of 40 RNA X-ray structures solved to <3.0 A resolution. For a typical number (>or=11) of one bond C-H base and sugar rdcs, these structural deviations together with rdc uncertainty (1.5 Hz) lead to average errors in the magnitude and orientation of the principal axis of order that are <9% and <4 degrees, respectively. The errors decrease to <5% and <4 degrees for >or=17 rdcs. A protocol that allows for estimation of error in A-form order tensors due to both angular deviations and rdc uncertainty (Aform-RDC) is validated using theoretical simulations and used to analyze rdcs measured previously in TAR in the free state and bound to four distinct ligands. Results confirm earlier findings that the two TAR helices undergo large changes in both their mean relative orientation and dynamics upon binding to different targets.

Dynamics-based amplification of RNA function and its characterization by using NMR spectroscopy

The ever-increasing cellular roles ascribed to RNA raise fundamental questions regarding how a biopolymer composed of only four chemically similar building-block nucleotides achieves such functional diversity. Here, I discuss how RNA achieves added mechanistic and chemical complexity by undergoing highly controlled conformational changes in response to a variety of cellular signals. I examine pathways for achieving selectivity in these conformational changes that rely to different extents on the structure and dynamics of RNA. Finally, I review solution-state NMR techniques that can be used to characterize RNA structural dynamics and its relationship to function.

Altered specificity in DNA binding by the lac repressor: a mutant lac headpiece that mimics the gal repressor

Recognition of the lac operator by the lac repressor involves specific interactions between residues in the repressor's recognition helix and bases in the DNA major groove. Tyr17 and Gln18, at positions 1 and 2 in the lac repressor recognition helix, can be exchanged for other amino acids to generate mutant repressors that display altered specificity. We have solved the solution structure of a protein-DNA complex of an altered-specificity mutant lac headpiece in which Tyr17 and Gln18 were exchanged for valine and alanine, respectively, as found in the recognition helix of the gal repressor. As previously described by Lehming et al. (EMBO J. 1987, 6, 3145-3153), this altered-specificity mutant of the lac repressor recognizes a variant lac operator that is similar to the gal operator Oe. The mutant lac headpiece showed the predicted specificity and is also able to mimic the gal repressor by recognizing and bending the natural gal operator Oe. These structural data show that, while most of the anchoring points that help the lac headpiece to assemble on the lac operator were preserved, a different network of protein-DNA interactions connecting Ala17 and Val18 to bases in the DNA major groove drives the specificity towards the altered operator.

NMR methods for studying the structure and dynamics of RNA

Proper functioning of RNAs requires the formation of complex three-dimensional structures combined with the ability to rapidly interconvert between multiple functional states. This review covers recent advances in isotope-labeling strategies and NMR experimental approaches that have promise for facilitating solution structure determinations and dynamics studies of biologically active RNAs. Improved methods for the production of isotopically labeled RNAs combined with new multidimensional heteronuclear NMR experiments make it possible to dramatically reduce spectral crowding and simplify resonance assignments for RNAs. Several novel applications of experiments that directly detect hydrogen-bonding interactions are discussed. These studies demonstrate how NMR spectroscopy can be used to distinguish between possible secondary structures and identify mechanisms of ligand binding in RNAs. A variety of recently developed methods for measuring base and sugar residual dipolar couplings are described. NMR residual dipolar coupling techniques provide valuable data for determining the long-range structure and orientation of helical regions in RNAs. A number of studies are also presented where residual dipolar coupling constraints are used to determine the global structure and dynamics of RNAs. NMR relaxation data can be used to probe the dynamics of macromolecules in solution. The power dependence of transverse rotating-frame relaxation rates was used here to study dynamics in the minimal hammerhead ribozyme. Improved methods for isotopically labeling RNAs combined with new types of structural data obtained from a growing repertoire of NMR experiments are facilitating structural and dynamic studies of larger RNAs.

205Tl NMR Methods for the Characterization of Monovalent Cation Binding to Nucleic Acids

Monovalent cations play an important role in many biological functions. The guanine rich sequence, d(G4T4G4), requires monovalent cations for formation of the G-quadruplex, d(G4T4G4)2. This requirement can be satisfied by thallium (Tl+), a potassium (K+) surrogate. To verify that the structure of d(G4T4G4)2 in the presence of Tl+ is similar to the K+-form of the G-quadruplex, the solution structure of the Tl+-form of d(G4T4G4)2 was determined. The 10 lowest energy structures have an all atom RMSD of 0.76 +/- 0.16 A. Comparison of this structure to the identical G-quadruplex formed in the presence of K+ validates the isomorphous nature of Tl+ and K+. Using a 1H-205Tl spin-echo difference experiment we show that, in the Tl+-form of d(G4T4G4)2, small scalar couplings (<1 Hz) exist between 205Tl and protons in the G-quadruplex. These data comprise the first 1H-205Tl scalar couplings observed in a biological system and have the potential to provide important constraints for structure determination. These experiments can be applied to any system in which the substituted Tl+ cations are in slow exchange with the bulk ions in solution.

Residue specific ribose and nucleobase dynamics of the cUUCGg RNA tetraloop motif by MNMR 13C relaxation

The dynamics of the nucleobase and the ribose moieties in a 14-nt RNA cUUCGg hairpin-loop uniformly labeled with 13C and 15N were studied by 13C spin relaxation experiments. R1, R1rho and the 13C-[1H] steady-state NOE of C6 and C1' in pyrimidine and C8 and C1' in purine residues were obtained at 298 K. The relaxation data were analyzed by the model-free formalism to yield dynamic information on timescales of pico-, nano- and milli-seconds. An axially symmetric diffusion tensor with an overall rotational correlation time tau(c) of 2.31 +/- 0.13 ns and an axial ratio of 1.35 +/- 0.02 were determined. Both findings are in agreement with hydrodynamic calculations. For the nucleobase carbons, the validity of different reported 13C chemical shift anisotropy values (Stueber, D. and Grant, D. M., 2002 J. Am. Chem. Soc. 124, 10539-10551; Fiala et al., 2000 J. Biomol. NMR 16, 291-302; Sitkoff, D. and Case, D. A., 1998 Prog. NMR Spectroscopy 32, 165-190) is discussed. The resulting dynamics are in agreement with the structural features of the cUUCGg motif in that all residues are mostly rigid (0.82 < S2 < 0.96) in both the nucleobase and the ribose moiety except for the nucleobase of U7, which is protruding into solution (S2 = 0.76). In general, ribose mobility follows nucleobase dynamics, but is less pronounced. Nucleobase dynamics resulting from the analysis of 13C relaxation rates were found to be in agreement with 15N relaxation data derived dynamic information (Akke et al., 1997 RNA 3, 702-709).

Measurement of Ribose Carbon Chemical Shift Tensors for A-form RNA by Liquid Crystal NMR Spectroscopy

Incomplete motional averaging of chemical shift anisotropy upon weak alignment of nucleic acids and proteins in a magnetic field results in small changes in chemical shift. Knowledge of nucleus-specific chemical shift (CS) tensor magnitudes and orientations is necessary to take full advantage of these measurements in biomolecular structure determination. We report the determination by liquid crystal NMR of the CS tensors for all ribose carbons in A-form helical RNA, using a series of novel 3D NMR pulse sequences for accurate and resolved measurement of the ribose (13)C chemical shifts. The orientation of the riboses relative to the rhombic alignment tensor of the molecule studied, a stem-loop sequence corresponding to helix-35 of 23S rRNA, is known from an extensive set of residual dipolar couplings (RDC), previously used to refine its structure. Singular-value-decomposition fits of the chemical shift changes to this structure, or alternatively to a database of helical RNA X-ray structures, provide the CS tensor for each type of carbon. Quantum chemical calculations complement the experimental results and confirm that the most shielded tensor component lies approximately along the local carbon-oxygen bond axis in all cases and that shielding anisotropy for C3' and C4' is much larger than for C1' and C2', with C5' being intermediate.