Time-resolved synchrotron X-ray "footprinting", a new approach to the study of nucleic acid structure and function: application to protein-DNA interactions and RNA folding

Characterization of tertiary folding of RNA by circular dichroism and urea

CD spectroscopy can be used to monitor RNA tertiary folding transitions that may not be observable by absorbance spectroscopy. With the use of computer-controlled titrators, data can be acquired rapidly, and accurate thermodynamic properties can be obtained over a wide variety of conditions. Thus, CD spectroscopy provides a useful complement to site-resolved or chemical modification methods.

Quantitative studies of ribosome conformational dynamics

The ribosome is a dynamic machine that undergoes many conformational rearrangements during the initiation of protein synthesis. Significant differences exist between the process of protein synthesis initiation in eubacteria and eukaryotes. In particular, the initiation of eukaryotic protein synthesis requires roughly an order of magnitude more initiation factors to promote efficient mRNA recruitment and ribosomal recognition of the start codon than are needed for eubacterial initiation. The mechanisms by which these initiation factors promote ribosome conformational changes during stages of initiation have been studied using cross-linking, footprinting, site-directed probing, cryo-electron microscopy, X-ray crystallography, fluorescence spectroscopy and single-molecule techniques. Here, we review how the results of these different approaches have begun to converge to yield a detailed molecular understanding of the dynamic motions that the eukaryotic ribosome cycles through during the initiation of protein synthesis.

Following molecular transitions with single residue spatial and millisecond time resolution

"Footprinting" describes assays in which ligand binding or structure formation protects polymers such as nucleic acids and proteins from either cleavage or modification; footprinting allows the accessibility of individual residues to be mapped in solution. Equilibrium and time-dependent footprinting links site-specific structural information with thermodynamic and kinetic transitions, respectively. The hydroxyl radical (*OH) is a uniquely insightful footprinting probe by virtue of it being among the most reactive chemical oxidants; it reports the solvent accessibility of reactive sites on macromolecules with as fine as a single residue resolution. A novel method of millisecond time-resolved *OH footprinting is presented based on the Fenton reaction, Fe(II) + H(2)O(2) --> Fe(III) + *OH + OH(-). It is implemented using a standard three-syringe quench-flow mixer. The utility of this method is demonstrated by its application to the studies on RNA folding. Its applicability to a broad range of biological questions involving the function of DNA, RNA, and proteins is discussed.

RNA misfolding and the action of chaperones

RNA folds to a myriad of three-dimensional structures and performs an equally diverse set of functions. The ability of RNA to fold and function in vivo is all the more remarkable because, in vitro, RNA has been shown to have a strong propensity to adopt misfolded, non-functional conformations. A principal factor underlying the dominance of RNA misfolding is that local RNA structure can be quite stable even in the absence of enforcing global tertiary structure. This property allows non-native structure to persist, and it also allows native structure to form and stabilize non-native contacts or non-native topology. In recent years it has become clear that one of the central reasons for the apparent disconnect between the capabilities of RNA in vivo and its in vitro folding properties is the presence of RNA chaperones, which facilitate conformational transitions of RNA and therefore mitigate the deleterious effects of RNA misfolding. Over the past two decades, it has been demonstrated that several classes of non-specific RNA binding proteins possess profound RNA chaperone activity in vitro and when overexpressed in vivo, and at least some of these proteins appear to function as chaperones in vivo. More recently, it has been shown that certain DExD/H-box proteins function as general chaperones to facilitate folding of group I and group II introns. These proteins are RNA-dependent ATPases and have RNA helicase activity, and are proposed to function by using energy from ATP binding and hydrolysis to disrupt RNA structure and/or to displace proteins from RNA-protein complexes. This review outlines experimental studies that have led to our current understanding of the range of misfolded RNA structures, the physical origins of RNA misfolding, and the functions and mechanisms of putative RNA chaperone proteins.

Folding of a universal ribozyme: the ribonuclease P RNA

Ribonuclease P is among the first ribozymes discovered, and is the only ubiquitously occurring ribozyme besides the ribosome. The bacterial RNase P RNA is catalytically active without its protein subunit and has been studied for over two decades as a model system for RNA catalysis, structure and folding. This review focuses on the thermodynamic, kinetic and structural frameworks derived from the folding studies of bacterial RNase P RNA.

Time-resolved NMR methods resolving ligand-induced RNA folding at atomic resolution

Structural transitions of RNA between alternate conformations with similar stabilities are associated with important aspects of cellular function. Few techniques presently exist that are capable of monitoring such transitions and thereby provide insight into RNA dynamics and function at atomic resolution. Riboswitches are found in the 5'-UTR of mRNA and control gene expression through structural transitions after ligand recognition. A time-resolved NMR strategy was established in conjunction with laser-triggered release of the ligand from a photocaged derivative in situ to monitor the hypoxanthine-induced folding of the guanine-sensing riboswitch aptamer domain of the Bacillus subtilis xpt-pbuX operon at atomic resolution. Combining selective isotope labeling of the RNA with NMR filter techniques resulted in significant spectral resolution and allowed kinetic analysis of the buildup rates for individual nucleotides in real time. Three distinct kinetic steps associated with the ligand-induced folding were delineated. After initial complex encounter the ligand-binding pocket is formed and results in subsequent stabilization of a remote long-range loop-loop interaction. Incorporation of NMR data into experimentally restrained molecular dynamics simulations provided insight into the RNA structural ensembles involved during the conformational transition.

Communication between RNA folding domains revealed by folding of circularly permuted ribozymes

To study the role of sequence and topology in RNA folding, we determined the kinetic folding pathways of two circularly permuted variants of the Tetrahymena group I ribozyme, using time-resolved hydroxyl radical footprinting. Circular permutation changes the distance between interacting residues in the primary sequence, without changing the native structure of the RNA. In the natural ribozyme, tertiary interactions in the P4-P6 domain form in 1 s, while interactions in the P3-P9 form in 1-3 min at 42 degrees C. Permutation of the 5' end to G111 in the P4 helix allowed the stable P4-P6 domain to fold in 200 ms at 30 degrees C, five times faster than in the wild-type RNA, while the other domains folded five times more slowly (5-8 min). By contrast, circular permutation of the 5' end to G303 in J8/7 decreased the folding rate of the P4-P6 domain. In this permuted RNA, regions joining P2, P3 and P4 were protected in 500 ms, while the P3-P9 domain was 60-80% folded within 30 s. RNase T(1) digestion and FMN photocleavage showed that circular permutation of the RNA sequence alters the initial ensemble of secondary structures, thereby changing the tertiary folding pathways. Our results show that the natural 5'-to-3' order of the structural domains in group I ribozymes optimizes structural communication between tertiary domains and promotes self-assembly of the catalytic center.

Distinct contribution of electrostatics, initial conformational ensemble, and macromolecular stability in RNA folding

We distinguish the contribution of the electrostatic environment, initial conformational ensemble, and macromolecular stability on the folding mechanism of a large RNA using a combination of time-resolved "Fast Fenton" hydroxyl radical footprinting and exhaustive kinetic modeling. This integrated approach allows us to define the folding landscape of the L-21 Tetrahymena thermophila group I intron structurally and kinetically from its earliest steps with unprecedented accuracy. Distinct parallel pathways leading the RNA to its native form upon its Mg(2+)-induced folding are observed. The structures of the intermediates populating the pathways are not affected by variation of the concentration and type of background monovalent ions (electrostatic environment) but are altered by a mutation that destabilizes one domain of the ribozyme. Experiments starting from different conformational ensembles but folding under identical conditions show that whereas the electrostatic environment modulates molecular flux through different pathways, the initial conformational ensemble determines the partitioning of the flux. This study showcases a robust approach for the development of kinetic models from collections of local structural probes.

Electrophoretic mobility shift assay (EMSA) for detecting protein-nucleic acid interactions

The gel electrophoresis mobility shift assay (EMSA) is used to detect protein complexes with nucleic acids. It is the core technology underlying a wide range of qualitative and quantitative analyses for the characterization of interacting systems. In the classical assay, solutions of protein and nucleic acid are combined and the resulting mixtures are subjected to electrophoresis under native conditions through polyacrylamide or agarose gel. After electrophoresis, the distribution of species containing nucleic acid is determined, usually by autoradiography of 32P-labeled nucleic acid. In general, protein-nucleic acid complexes migrate more slowly than the corresponding free nucleic acid. In this protocol, we identify the most important factors that determine the stabilities and electrophoretic mobilities of complexes under assay conditions. A representative protocol is provided and commonly used variants are discussed. Expected outcomes are briefly described. References to extensions of the method and a troubleshooting guide are provided.

Unusual metal specificity and structure of the group I ribozyme from Chlamydomonas reinhardtii 23S rRNA

Group I intron ribozymes require cations for folding and catalysis, and the current literature indicates that a number of cations can promote folding, but only Mg2+ and Mn2+ support both processes. However, some group I introns are active only with Mg2+, e.g. three of the five group I introns in Chlamydomonas reinhardtii. We have investigated one of these ribozymes, an intron from the 23S LSU rRNA gene of Chlamydomonas reinhardtii (Cr.LSU), by determining if the inhibition by Mn2+ involves catalysis, folding, or both. Kinetic analysis of guanosine-dependent cleavage by a Cr.LSU ribozyme, 23S.5 Delta Gb, that lacks the 3' exon and intron-terminal G shows that Mn2+ does not affect guanosine binding or catalysis, but instead promotes misfolding of the ribozyme. Surprisingly, ribozyme misfolding induced by Mn2+ is highly cooperative, with a Hill coefficient larger than that of native folding induced by Mg2+. At lower Mn2+ concentrations, metal inhibition is largely alleviated by the guanosine cosubstrate (GMP). The concentration dependence of guanosine cosubstrate-induced folding suggests that it functions by interacting with the G binding site, perhaps by displacing an inhibitory Mn2+. Because of these and other properties of Cr.LSU, the tertiary structure of the intron from 23S.5 Delta Gb was examined using Fe2+-EDTA cleavage. The ground-state structure shows evidence of an unusually open ribozyme core: the catalytic P3-P7 domain and the nucleotides that connect it to the P4-P5-P6 domain are exposed to solvent. The implications of this structure for the in vitro and in vivo properties of this intron ribozyme are discussed.

Mg2+-RNA interaction free energies and their relationship to the folding of RNA tertiary structures

Mg2+ ions are very effective at stabilizing tertiary structures in RNAs. In most cases, folding of an RNA is so strongly coupled to its interactions with Mg2+ that it is difficult to separate free energies of Mg2+-RNA interactions from the intrinsic free energy of RNA folding. To devise quantitative models accounting for this phenomenon of Mg2+-induced RNA folding, it is necessary to independently determine Mg2+-RNA interaction free energies for folded and unfolded RNA forms. In this work, the energetics of Mg2+-RNA interactions are derived from an assay that measures the effective concentration of Mg2+ in the presence of RNA. These measurements are used with other measures of RNA stability to develop an overall picture of the energetics of Mg2+-induced RNA folding. Two different RNAs are discussed, a pseudoknot and an rRNA fragment. Both RNAs interact strongly with Mg2+ when partially unfolded, but the two folded RNAs differ dramatically in their inherent stability in the absence of Mg2+ and in the free energy of their interactions with Mg2+. From these results, it appears that any comprehensive framework for understanding Mg2+-induced stabilization of RNA will have to (i) take into account the interactions of ions with the partially unfolded RNAs and (ii) identify factors responsible for the widely different strengths with which folded tertiary structures interact with Mg2+.

Fast Fenton footprinting: a laboratory-based method for the time-resolved analysis of DNA, RNA and proteins

'Footprinting' describes assays in which ligand binding or structure formation protects polymers such as nucleic acids and proteins from either cleavage or modification; footprinting allows the accessibility of individual residues to be mapped in solution. Equilibrium and time-dependent footprinting links site-specific structural information with thermodynamic and kinetic transitions. The hydroxyl radical (*OH) is a particularly valuable footprinting probe by virtue of it being among the most reactive of chemical oxidants; it reports the solvent accessibility of reactive sites on macromolecules with as fine as a single residue resolution. A novel method of millisecond time-resolved .OH footprinting has been developed based on the Fenton reaction, Fe(II) + H2O2 --> Fe(III) + *OH + OH-. This method can be implemented in laboratories using widely available three-syringe quench flow mixers and inexpensive reagents to study local changes in the solvent accessibility of DNA, RNA and proteins associated with their biological function.

Local kinetic measures of macromolecular structure reveal partitioning among multiple parallel pathways from the earliest steps in the folding of a large RNA molecule

At the heart of the RNA folding problem is the number, structures, and relationships among the intermediates that populate the folding pathways of most large RNA molecules. Unique insight into the structural dynamics of these intermediates can be gleaned from the time-dependent changes in local probes of macromolecular conformation (e.g. reports on individual nucleotide solvent accessibility offered by hydroxyl radical (()OH) footprinting). Local measures distributed around a macromolecule individually illuminate the ensemble of separate changes that constitute a folding reaction. Folding pathway reconstruction from a multitude of these individual measures is daunting due to the combinatorial explosion of possible kinetic models as the number of independent local measures increases. Fortunately, clustering of time progress curves sufficiently reduces the dimensionality of the data so as to make reconstruction computationally tractable. The most likely folding topology and intermediates can then be identified by exhaustively enumerating all possible kinetic models on a super-computer grid. The folding pathways and measures of the relative flux through them were determined for Mg(2+) and Na(+)-mediated folding of the Tetrahymena thermophila group I intron using this combined experimental and computational approach. The flux during Mg(2+)-mediated folding is divided among numerous parallel pathways. In contrast, the flux during the Na(+)-mediated reaction is predominantly restricted through three pathways, one of which is without detectable passage through intermediates. Under both conditions, the folding reaction is highly parallel with no single pathway accounting for more than 50% of the molecular flux. This suggests that RNA folding is non-sequential under a variety of different experimental conditions even at the earliest stages of folding. This study provides a template for the systematic analysis of the time-evolution of RNA structure from ensembles of local measures that will illuminate the chemical and physical characteristics of each step in the process. The applicability of this analysis approach to other macromolecules is discussed.

Mg2+-induced compaction of single RNA molecules monitored by tethered particle microscopy

We have applied tethered particle microscopy (TPM) as a single molecule analysis tool to studies of the conformational dynamics of poly-uridine(U) messenger (m)RNA and 16S ribosomal (r)RNA molecules. Using stroboscopic total internal reflection illumination and rigorous selection criteria to distinguish from nonspecific tethering, we have tracked the nanometer-scale Brownian motion of RNA-tethered fluorescent microspheres in all three dimensions at pH 7.5, 22 degrees C, in 10 mM or 100 mM NaCl in the absence or presence of 10 mM MgCl(2). The addition of Mg(2+) to low-ionic strength buffer results in significant compaction and stiffening of poly(U) mRNA, but not of 16S rRNA. Furthermore, the motion of poly(U)-tethered microspheres is more heterogeneous than that of 16S rRNA-tethered microspheres. Analysis of in-plane bead motion suggests that poly(U) RNA, but less so 16S rRNA, can be modeled both in the presence and absence of Mg(2+) by a statistical Gaussian polymer model. We attribute these differences to the Mg(2+)-induced compaction of the relatively weakly structured and structurally disperse poly(U) mRNA, in contrast to Mg(2+)-induced reinforcement of existing secondary and tertiary structure contacts in the highly structured 16S rRNA. Both effects are nonspecific, however, as they are dampened in the presence of higher concentrations of monovalent cations.

Perturbation of the hierarchical folding of a large RNA by the destabilization of its Scaffold's tertiary structure

The P4-P6 domain serves as a scaffold against which the periphery and catalytic core organize and fold during Mg2+-mediated folding of the Tetrahymena thermophila ribozyme. The most prominent structural motif of the P4-P6 domain is the tetraloop-tetraloop receptor interaction which "clamps" the distal parts of its hairpin-like structure. Destabilization of the tertiary structure of the P4-P6 domain by perturbation of the tetraloop-tetraloop receptor interaction alters the Mg2+-mediated folding pathway. The folding hierarchy of P5c approximately P4-P6 > periphery > catalytic core that is a striking attribute of the folding of the wild-type RNA is abolished. The initial steps in folding of the mutant RNA are > or =50-fold faster than those of the wild-type ribozyme with the earliest observed tertiary contacts forming around regions known to specifically bind Mg2+. The interaction between the mutant tetraloop and the tetraloop receptor appears coincidently with slowly forming catalytic core tertiary contacts. Thus, the stability conferred upon the P4-P6 domain by the tetraloop-tetraloop receptor interaction dictates the preferred folding pathway by stabilizing an early intermediate. A sub-denaturing concentration of urea diminishes the early barrier to folding the wild-type ribozyme along with complex effects on the subsequent steps of folding the wild-type and mutant RNA.

Mapping nucleic acid structure by hydroxyl radical cleavage

Hydroxyl radical footprinting is a widely used method for following the folding of RNA molecules in solution. This method has the unique ability to provide experimental information on the solvent accessibility of each nucleotide in an RNA molecule, so that the folding of all domains of the RNA species can be followed simultaneously at single-nucleotide resolution. In recent work, hydroxyl radical footprinting has been used, often in combination with other global measures of structure, to work out detailed folding pathways and three-dimensional structures for increasingly large and complicated RNA molecules. These include synthetic ribozymes, and group I and group II ribozymes, from yeast, the Azoarcus cyanobacterium and Tetrahymena thermophila. Advances have been made in methods for analysis of hydroxyl radical data, so that the large datasets that result from kinetic folding experiments can be analyzed in a semi-automated and quantitative manner.

SAFA: semi-automated footprinting analysis software for high-throughput quantification of nucleic acid footprinting experiments

Footprinting is a powerful and widely used tool for characterizing the structure, thermodynamics, and kinetics of nucleic acid folding and ligand binding reactions. However, quantitative analysis of the gel images produced by footprinting experiments is tedious and time-consuming, due to the absence of informatics tools specifically designed for footprinting analysis. We have developed SAFA, a semi-automated footprinting analysis software package that achieves accurate gel quantification while reducing the time to analyze a gel from several hours to 15 min or less. The increase in analysis speed is achieved through a graphical user interface that implements a novel methodology for lane and band assignment, called "gel rectification," and an optimized band deconvolution algorithm. The SAFA software yields results that are consistent with published methodologies and reduces the investigator-dependent variability compared to less automated methods. These software developments simplify the analysis procedure for a footprinting gel and can therefore facilitate the use of quantitative footprinting techniques in nucleic acid laboratories that otherwise might not have considered their use. Further, the increased throughput provided by SAFA may allow a more comprehensive understanding of molecular interactions. The software and documentation are freely available for download at http://safa.stanford.edu.

Monovalent ion-mediated folding of the Tetrahymena thermophila ribozyme

The time-course of monovalent cation-induced folding of the L-21 Sca1 Tetrahymena thermophila ribozyme and a selected mutant was quantitatively followed using synchrotron X-ray (.OH) footprinting. Initiating folding by increasing the concentration of either Na+ or K+ to 1.5M from an initial condition of approximately 0.008 M Na+ at 42 degrees C resulted in the complete formation of tertiary contacts within the P5abc subdomain and between the peripheral helices within the dead time of our measurements (k>50 s(-1)). These results contrast with folding rates of 2-0.2 s(-1) previously observed for formation of these contacts in 10mM Mg2+ from the same initial condition. Thus, the initial formation of native tertiary contacts is inhibited by divalent but not monovalent cations. The native contacts within the catalytic core form without a detectable burst phase at rates of 0.4-1.0 s(-1) in a manner reminiscent of the Mg2+-dependent folding behavior, although tenfold faster. The tertiary interactions stabilizing the catalytic core interaction with P4-P6 and P2.1, as well as one of the protections internal for the P4-P6 domain, display progress curves with appreciable burst amplitudes and a phase comparable in rate to that of the catalytic core. That the slow folding of the ribozyme's core is a consequence of the alt-P3 secondary structure is shown by the 100% burst phase amplitudes that are observed for folding of the U273A mutant ribozyme within which the native secondary structure (P3) is strengthened. Thus, formation of a misfolded intermediate(s) resulting from the alt-P3 secondary structure is independent of ion valency while the rate at which the respective intermediates are resolved is sensitive to ion valency. The overall portrait painted by these results is that ion valency differentially affects steps in the folding process and that folding in monovalent ion alone for the U273A mutant Tetrahymena ribozyme is fast and direct.

Semi-automated, single-band peak-fitting analysis of hydroxyl radical nucleic acid footprint autoradiograms for the quantitative analysis of transitions

Hydroxyl radical footprinting can probe the solvent accessibility of the ribose moiety of the individual nucleotides of DNA and RNA. Semi-automated analytical tools are presented for the quantitative analyses of nucleic acid footprint transitions in which processes such as folding or ligand binding are followed as a function of time or ligand concentration. Efficient quantitation of the intensities of the electrophoretic bands comprising the footprinting reaction products is achieved by fitting a series of Lorentzian curves to line profiles obtained from gels utilizing sequentially relaxed constraints consistent with electrophoretic mobility. An automated process of data 'standardization' has been developed that corrects for differences in the loading amounts in the electrophoresis. This process enhances the accuracy of the derived transitions and makes generating them easier. Together with visualization of the processed footprinting in false-color two-dimensional maps, DNA and RNA footprinting data can be accurately, precisely and efficiently processed allowing transitions to be objectively and comprehensively analyzed. The utility of this new analysis approach is illustrated by its application to the ion-meditated folding of a large RNA molecule.

Analysis of Chemical and Enzymatic Cleavage Frequencies in Supercoiled DNA

Chemical and enzymatic probing methods are powerful techniques for examining details of sequence-dependent structure in DNA and RNA. Reagents that cleave nucleic acid molecules in a structure-specific, but relatively sequence-non-specific manner, such as hydroxyl radical or DNase I, have been used widely to probe helical geometry in nucleic acid structures, nucleic acid-drug complexes, and in nucleoprotein assemblies. Application of cleavage-based techniques to structures present in superhelical DNA has been hindered by the fact that the cleavage pattern attributable to supercoiling-dependent structures is heavily mixed with non-specific cleavage signals that are inevitable products of multiple cleavage events. We present a rigorous mathematical procedure for extracting the cleavage pattern specific to supercoiled DNA and use this method to investigate the hydroxyl radical cleavage pattern in a cruciform DNA structure formed by a 60 bp inverted repeat sequence embedded in a negatively supercoiled plasmid. Our results support the presence of a stem-loop structure in the expected location and suggest that the helical geometry of the cruciform stem differs from that of the normal duplex form.

Multiple monovalent ion-dependent pathways for the folding of the L-21 Tetrahymena thermophila ribozyme

Synchrotron hydroxyl radical (*OH) footprinting is a technique that monitors the local changes in solvent accessibility of the RNA backbone on milliseconds to minutes time-scales. The Mg(2+)-dependent folding of the L-21 Sca 1 Tetrahymena thermophila ribozyme has been followed using this technique at an elevated concentration of monovalent ion (200 mM NaCl) and as a function of the initial annealing conditions and substrate. Previous studies conducted at low concentrations of monovalent ion displayed sequential folding of the P4-P6 domain, the peripheral helices and the catalytic core, with each protection displaying monophasic kinetics. For ribozyme annealed in buffer containing 200 mM NaCl and folded by the addition of 10 mM MgCl(2), multiple kinetic phases are observed for *OH protections throughout the ribozyme. The independently folding P4-P6 domain is the first to fold with its protections displaying 50-90% burst phase amplitudes. That the folding of P4-P6 within the ribozyme does not display the 100% burst phase of isolated P4-P6 at 200 mM NaCl shows that interactions with the remainder of the ribozyme impede this domain's folding. In addition, *OH protections constituting each side of a tertiary contact are not coincident in some cases, consistent with the formation of transient non-native interactions. While the peripheral contacts and triple helical scaffold exhibit substantial burst phases, the slowest protection to appear is J8/7 in the catalytic core, which displays a minimal burst amplitude and whose formation is coincident with the recovery of catalytic activity. The number of kinetic phases as well as their amplitudes and rates are different when the ribozyme is annealed in low-salt buffer and folded by the concomitant addition of monovalent and divalent cations. Annealed substrate changes the partitioning of the ribozyme among the multiple folding populations. These results provide a map of the early steps in the ribozyme's folding landscape and the degree to which the preferred pathways are dependent upon the initial reaction conditions.

Structural changes of tRNA and 5S rRNA induced with magnesium and visualized with synchrotron mediated hydroxyl radical cleavage

The structure of native yeast tRNA(Phe) and wheat germ ribosomal 5S RNA induced by different magnesium ion concentrations was studied in solution with a synchrotron mediated hydroxyl radical RNA cleavage reaction. We showed that very small amounts of Mg+2 can induce significant changes in the hydroxyl radical cleavage pattern of tRNA(Phe). It also turned out that a reactivity of tRNAz(Phe) towards *OH coincides with the strong metal binding sites. Because of the Mg ions are heavily hydrated one can suggest the strong correlation of the observed nucleosides reactivity in vicinity of Mg2+ binding sites with availability of water molecules as a source of hydroxyl radical. On the other hand the structure of wheat germ 5S rRNA is less sensitive to the hydroxyl radical reaction than tRNA(Phe) although some changes are visible at 4 mM Mg ions. It is probably due to the lack of strong Mg+2 binding sites in that molecule. The reactivity of nucleotides in loops C and D of 5S rRNA is not effected, what suggests their flexibility or involvement in higher order structure formation. There is different effect of magnesium on tRNA and 5S rRNA folding. We found that nucleotides forming strong binding sites for magnesium are very sensitive to X-ray generated hydroxyl radical and can be mapped with *OH. The results show, that guanine nucleotides are preferentially hydrated. X-ray footprinting mediated hydroxyl radical RNA cleavage is a very powerful method and has been applied to studies of stable RNAs for the first time.

Mispaired P3 region in the hierarchical folding pathway of the Tetrahymena ribozyme

BACKGROUND

The Tetrahymena group I ribozyme folds into a complex three-dimensional structure for performing catalytic reactions. The catalysis depends on its catalytic core consisting of two helical domains, P4-P6 and P3-P7, connected by single stranded regions. In the folding process, most of this ribozyme folds in a hierarchical manner in which a kinetically stable intermediate determines the overall folding rate.

RESULTS

Although the nature of this intermediate has not yet been elucidated, a mispaired P3 stem (alt-P3) appears a likely candidate. To examine the effects of the alt-P3 structure on the kinetic and thermodynamic properties of the active structure of the ribozyme or its P3-P7 domain formation, we prepared and analysed variant ribozymes in which relative stabilities of the original P3 and alt-P3 structure were altered systematically.

CONCLUSION

The results indicate that the alt-P3 structure is not the major rate-limiting factor in the folding process.

Concerted binding and bending of DNA by Escherichia coli integration host factor

Integration host factor (IHF) is a heterodimeric Escherichia coli protein that plays essential roles in a variety of cellular processes including site-specific recombination, transcription, and DNA replication. The IHF-DNA interface extends over three helical turns and includes sequential minor groove contacts that present strong, sequence specific protection patterns against hydroxyl radical cleavage. Synchrotron X-ray footprinting has been used to follow the kinetics of formation of DNA-protein contacts in the IHF-DNA complex with single base-pair spatial, and millisecond time, resolution. The three sites of IHF protection on the DNA develop with similar time-dependence, indicating that sequence specific binding and bending occur concertedly. Two distinct phases are observed in the association process. The first "burst" phase is characterized by a rate that is greater than diffusion limited (>10(10) s(-1) M(-1)) and the second phase is on the order of diffusion controlled (approximately 10(8) M(-1) s(-1)). The overall kinetics of association become faster with increasing IHF concentration showing that complex formation is second-order with protein. The rate of association is maximal between 100 and 200 mM KCl decreasing at higher and lower concentrations. The rate of IHF dissociation from site-specifically bound DNA increases monotonically as KCl concentration is increased. The dissociation progress curves are biphasic with the amplitude of the first phase dependent upon competitor DNA concentration. These results are the first analysis by synchrotron footprinting of the fast kinetics of a protein-DNA interaction and suggest that IHF binds its specific site through a multiple-step mechanism in which the first step is facilitated diffusion along the length of the duplex followed by subsequent binding and bending of the DNA in a concerted manner.

Productive folding to the native state by a group II intron ribozyme

Group II introns are large catalytic RNA molecules that fold into compact structures essential for the catalysis of splicing and intron mobility reactions. Despite a growing body of information on the folded state of group II introns at equilibrium, there is currently no information on the folding pathway and little information on the ionic requirements for folding. Folding isotherms were determined by hydroxyl radical footprinting for the 32 individual protections that are distributed throughout a group II intron ribozyme derived from intron ai5gamma. The isotherms span a similar range of Mg(2+) concentrations and share a similar index of cooperativity. Time-resolved hydroxyl radical footprinting studies show that all regions of the ribozyme fold slowly and with remarkable synchrony into a single catalytically active structure at a rate comparable to those of other ribozymes studied thus far. The rate constants for the formation of tertiary contacts and recovery of catalytic activity are identical within experimental error. Catalytic activity analyses in the presence of urea provide no evidence that the slow folding of the ai5gamma intron is attributable to the presence of unproductive kinetic traps along the folding pathway. Taken together, the data suggest that the rate-limiting step for folding of group II intron ai5gamma occurs early along the reaction pathway. We propose that this behavior resembles protein folding that is limited in rate by high contact order, or the need to form key tertiary interactions from partners that are located far apart in the primary or secondary structure.

Studies of biomolecular conformations and conformational dynamics by mass spectrometry

In the post-genomic era, a wealth of structural information has been amassed for proteins from NMR and crystallography. However, static protein structures alone are not a sufficient description: knowledge of the dynamic nature of proteins is essential to understand their wide range of functions and behavior during the life cycle from synthesis to degradation. Furthermore, few proteins have the ability to act alone in the crowded cellular environment. Assemblies of multiple proteins governed by complex signaling pathways are often required for the tasks of target recognition, binding, transport, and function. Mass spectrometry has emerged over the past several years as a powerful tool to address many of these questions. Recent improvements in "soft" ionization techniques have enabled researchers to study proteins and biomolecular complexes, both directly and indirectly. Likewise, continuous improvements in instrumental design in recent years have resulted in a dramatic expansion of the m/z range and resolution, enabling observation of large multi-protein assemblies whose structures are retained in the gas phase. In this article, we discuss some of the mass spectrometric techniques applied to investigate the nature of the conformations and dynamical properties that govern protein function.

Nucleic acid fragmentation on the millisecond timescale using a conventional X-ray rotating anode source: application to protein-DNA footprinting

Nucleic acid fragmentation (footprinting) by *OH radicals is used often as a tool to probe nucleic acid structure and nucleic acid-protein interactions. This method has proven valuable because it provides structural information with single base pair resolution. Recent developments in the field introduced the 'synchrotron X-ray footprinting' method, which uses a high-flux X-ray source to produce single base pair fragmentation of nucleic acid in tens of milliseconds. We developed a complementary method that utilizes X-rays generated from a conventional rotating anode machine in which nucleic acid footprints can be generated by X-ray exposures as short as 100-300 ms. Our theoretical and experimental studies indicate that efficient cleavage of nucleic acids by X-rays depends upon sample preparation, energy of the X-ray source and the beam intensity. In addition, using this experimental set up, we demonstrated the feasibility of conducting X-ray footprinting to produce protein-DNA protection portraits at sub-second timescales.

Synchrotron protein footprinting: a technique to investigate protein-protein interactions

Traditional approaches for macromolecular structure elucidation, including NMR, crystallography and cryo-EM have made significant progress in defining the structures of protein-protein complexes. A substantial number of macromolecular structures, however, have not been examined with atomic detail due to sample size and heterogeneity, or resolution limitations of the technique; therefore, the general applicability of each method is greatly reduced. Synchrotron footprinting attempts to bridge the gap in these methods by monitoring changes in accessible surface areas of discrete macromolecular moieties. As evidenced by our previous studies on RNA folding and DNA-protein interactions, the three-dimensional structure is probed by examining the reactions of these moieties with hydroxyl radicals generated by synchrotron X-rays. Here we report the application of synchrotron footprinting to the investigation of protein- protein interactions, as the novel technique has been utilized to successfully map the contact sites of gelsolin segment-1 in the gelsolin segment 1/actin complex. Footprinting results demonstrate that phenylalanine 104, located on the actin binding helix of gelsolin segment 1, is protected from hydroxyl radical modification in the presence of actin. This change in reactivity results from the specific protection of gelsolin segment-1, consistent with the substantial decrease in solvent accessibility of F104 upon actin binding, as calculated from the crystal structural of the gelsolin segment 1/actin complex. The results presented here establish synchrotron footprinting as a broadly applicable method to probe structural features of macromolecular complexes that are not amenable to conventional approaches.

Rapid magnesium chelation as a method to study real-time tertiary unfolding of RNA

This unit describes a method to measure the unfolding of RNA tertiary structure on a millisecond time scale. A stopped-flow spectrophotometer is used to measure the rate of unfolding induced by the addition of EDTA to an RNA whose tertiary structure has been stabilized in the presence of magnesium ions. Using this methodology, rate constants for unfolding of tertiary or secondary structure can be obtained over a range of temperatures, and these values can be used to construct Arrhenius and Eyring plots, from which activation energy, Arrhenius pre-exponential factor, and enthalpy and entropy of activation can be obtained. These data provide information about the energy of the transition state and the energy barriers between secondary and tertiary structure, which is necessary for predicting RNA tertiary structure from secondary structure.

Radical approaches to probe protein structure, folding, and interactions by mass spectrometry

This review describes mass spectrometry-based strategies and investigations to determine protein structure, folding dynamics, and protein-protein interactions in solution through the use of radical reagents. The radicals are generated in high flux within microseconds from synchrotron radiation and discharge sources, and react with proteins on time scales that are less than those often attributed to structural reorganization and folding. The oxygen-based radicals generated in aqueous solution react with proteins to effect limited oxidation at specific amino acids throughout the sequence of the protein. The extent of oxidation at these residue markers is highly influenced by the accessibility of the reaction site to the bulk solvent. The extent of oxidation allows protection levels to be measured based on the degree to which a reaction occurs. A map of a protein's three-dimensional structure is subsequently assembled as in a footprinting experiment. Protein solutions that contain various concentrations of substrates that either promote or disrupt dynamic structural transitions can be investigated to facilitate site-specific equilibrium and time-resolved studies of protein folding. The radical-based strategies can also be employed in the study of protein-protein interactions to provide a new avenue for investigating protein complexes and assemblies with high structural resolution. The urea-induced unfolding of apomyoglobin and the binding of gelsolin to actin are among the systems presented to illustrate the approach.

Mechanisms of photoinitiated cleavage of DNA by 1,8-naphthalimide derivatives

Using water-soluble 1,8-naphthalimide derivatives, the mechanisms of photosensitized DNA damage have been elucidated. Specifically, a comparison of rate constants for the photoinduced relaxation of supercoiled to circular DNA, as a function of dissolved halide, oxygen and naphthalimide concentration, has been carried out. The singlet excited states of the naphthalimide derivatives were quenched by chloride, bromide and iodide. In all cases the quenching products were naphthalimide triplet states, produced by induced intersystem crossing within the collision complex. Similarly, the halides were found to quench the triplet excited state of the 1,8-naphthalimide derivatives by an electron transfer mechanism. Bimolecular rate constants were < 10(5) M-1 s-1 for quenching by bromide and chloride. As expected from thermodynamic considerations quenching by iodide was 6.7 x 10(9) and 8.8 x 10(9) M-1 s-1 for the two 1,8-naphthalimide derivatives employed. At sufficiently high ground-state concentration self-quenching of the naphthalimide triplet excited state also occurs. The photosensitized conversion of supercoiled to circular DNA is fastest when self-quenching reactions are favored. The results suggest that, in the case of 1,8-naphthalimide derivatives, radicals derived from quenching of the triplet state by ground-state chromophores are more effective in cleaving DNA than reactive oxygen species or radicals derived from halogen atoms.

Probing the folding landscape of the Tetrahymena ribozyme: commitment to form the native conformation is late in the folding pathway

Large, structured RNAs traverse folding landscapes in which intermediates and long-lived misfolded states are common. To obtain a comprehensive description of the folding landscape for a structured RNA, it is necessary to understand the connections between productive folding pathways and pathways to these misfolded states. The Tetrahymena group I ribozyme partitions between folding to the native state and to a long-lived misfolded conformation. Here, we show that the observed rate constant for commitment to fold to the native or misfolded states is 1.9 min(-1) (37 degrees C, 10 mM Mg(2+)), the same within error as the rate constant for overall folding to the native state. Thus, the commitment to alternative folding pathways is made late in the folding process, concomitant with or after the rate-limiting step for overall folding. The ribozyme forms much of its tertiary structure significantly faster than it reaches this commitment point and the tertiary structure is expected to be stable, suggesting that the commitment to fold along pathways to the native or misfolded states is made from a partially structured intermediate. These results allow the misfolded conformation to be incorporated into a folding framework that reconciles previous data and gives quantitative information about the energetic topology of the folding landscape for this RNA.

Quantitative analysis of the isolated GAAA tetraloop/receptor interaction in solution: a site-directed spin labeling study

The GNRA (N: any nucleotide; R: purine) tetraloop/receptor interaction is believed to be one of the most frequently occurring tertiary interaction motifs in RNAs, but an isolated tetraloop/receptor complex has not been identified in solution. In the present work, site-directed spin labeling is applied to detect tetraloop/receptor complex formation and estimate the free energy of interaction. For this purpose, the GAAA tetraloop/receptor interaction was chosen as a model system. A method was developed to place nitroxide labels at specific backbone locations in an RNA hairpin containing the GAAA tetraloop. Formation of the tetraloop/receptor complex was monitored through changes in the rotational correlation time of the tetraloop and the attached nitroxide. Results show that a hairpin containing the GAAA tetraloop forms a complex with an RNA containing the 11-nucleotide GAAA tetraloop receptor motif with an apparent Kd that is strongly dependent on Mg2+. At 125 mM MgCl2, Kd = 0.40 +/- 0.05 mM. The corresponding standard free energy of complex formation is -4.6 kcal/mol, representing the energetics of the tetraloop/receptor interaction in the absence of other tertiary constraints. The experimental strategy presented here should have broad utility in quantifying weak interactions that would otherwise be undetectable, for both nucleic acids and nucleic acid-protein complexes.

Visualizing the solvent-inaccessible core of a group II intron ribozyme

Group II introns are well recognized for their remarkable catalytic capabilities, but little is known about their three-dimensional structures. In order to obtain a global view of an active enzyme, hydroxyl radical cleavage was used to define the solvent accessibility along the backbone of a ribozyme derived from group II intron ai5gamma. These studies show that a highly homogeneous ribozyme population folds into a catalytically compact structure with an extensively internalized catalytic core. In parallel, a model of the intron core was built based on known tertiary contacts. Although constructed independently of the footprinting data, the model implicates the same elements for involvement in the catalytic core of the intron.

A deteriorated triple-helical scaffold accelerates formation of the Tetrahymena ribozyme active structure

The Tetrahymena group I ribozyme requires a hierarchical folding process to form its correct three-dimensional structure. Ribozyme activity depends on the catalytic core consisting of two domains, P4-P6 and P3-P7, connected by a triple-helical scaffold. The folding proceeds in the following order: (i) fast folding of the P4-P6 domain, (ii) slow folding of the P3-P7 domain, and (iii) structure rearrangement to form the active ribozyme structure. The third step is believed to directly determine the conformation of the active catalytic domain, but as yet the precise mechanisms remain to be elucidated. To investigate the folding kinetics of this step, we analyzed mutant ribozymes having base substitution(s) in the triple-helical scaffold and found that disruption of the scaffold at A105G results in modest slowing of the P3-P7 folding (1.9-fold) and acceleration of step (iii) by 5.9-fold. These results suggest that disruption or destabilization of the scaffold is a normal component in the formation process of the active structure of the wild type ribozyme.

Role of counterion condensation in folding of the Tetrahymena ribozyme. I. Equilibrium stabilization by cations

Folding of RNA into an ordered, compact structure requires substantial neutralization of the negatively charged backbone by positively charged counterions. Using a native gel electrophoresis assay, we have examined the effects of counterion condensation upon the equilibrium folding of the Tetrahymena ribozyme. Incubation of the ribozyme in the presence of mono-, di- and trivalent ions induces a conformational state that is capable of rapidly forming the native structure upon brief exposure to Mg2+. The cation concentration dependence of this transition is directly correlated with the charge of the counterion used to induce folding. Substrate cleavage assays confirm the rapid onset of catalytic activity under these conditions. These results are discussed in terms of classical counterion condensation theory. A model for folding is proposed which predicts effects of charge, ionic radius and temperature on counterion-induced RNA folding transitions.

Time-resolved hydroxyl-radical footprinting of RNA using Fe(II)-EDTA

Chemical footprinting methods have been used extensively to probe the structures of biologically important RNAs at nucleotide resolution. One of these methods, hydroxyl-radical footprinting, has recently been employed to study the kinetics of RNA folding. Hydroxyl radicals can be generated by a number of different methods, including Fe(II)-EDTA complexes, synchrotron radiation, and peroxynitrous acid disproportionation. The latter two methods have been used for kinetic studies of RNA folding. We have taken advantage of rapid hydroxyl-radical generation by Fe(II)-EDTA-hydrogen peroxide solutions to develop a benchtop method to study folding kinetics of RNA complexes. This technique can be performed using commercially available chemicals, and can be used to accurately define RNA folding rate constants slower than 6 min(-1). Here we report the method and an example of time-resolved footprinting on the hairpin ribozyme, a small endoribonuclease and RNA ligase.

Determination of macromolecular folding and structure by synchrotron x-ray radiolysis techniques

Radiolysis of water by synchrotron X-rays generates oxygen-containing radicals that undergo reactions with solvent accessible sites of macromolecules inducing stable covalent modifications or cleavage on millisecond time scales. The extent and site of these reactions are determined by gel electrophoresis and mass spectrometry analysis. These data are used to construct a high-resolution map of solvent accessibility at individual reactive sites. The experiments can be performed in a time-resolved manner to provide kinetic rate constants for dynamic events occurring at individual sites within macromolecules or can provide equilibrium parameters of binding and thermodynamics of folding processes. The application of this synchrotron radiolysis technique to the study of lysozyme protein structure and the equilibrium urea induced unfolding of apomyoglobin are described. The Mg2+-induced folding of Tetrahymena thermophila group I ribozyme shows the capability of the method to study kinetics of folding.

Concerted kinetic folding of a multidomain ribozyme with a disrupted loop-receptor interaction

The free energy landscape for the folding of large, multidomain RNAs is rugged, and kinetically trapped, misfolded intermediates are a hallmark of RNA folding reactions. Here, we examine the role of a native loop-receptor interaction in determining the ruggedness of the energy landscape for folding of the Tetrahymena ribozyme. We demonstrate a progressive smoothing of the energy landscape for ribozyme folding as the strength of the loop-receptor interaction is reduced. Remarkably, with the most severe mutation, global folding is more rapid than for the wild-type ribozyme and proceeds in a concerted fashion without the accumulation of long-lived kinetic intermediates. The results demonstrate that a complex interplay between native tertiary interactions, divalent ion concentration, and non-native secondary structure determines the ruggedness of the energy landscape. Furthermore, the results suggest that kinetic folding transitions involving large regions of highly structured RNAs can proceed in a concerted fashion, in the absence of significant stable, preorganized tertiary structure.

Kinetic analysis of the M1 RNA folding pathway

The biological activity of large RNAs is dependent on the formation of complex folded structures that determine function. Typically the creation of such structures requires divalent magnesium and in many cases the folding process takes place over the course of several minutes. It has been proposed that the folding paths of large RNAs proceed through discrete intermediates but the nature of these intermediates is not known in most cases. Here, we describe our studies on the folding of the M1 RNA sub-unit of Escherichia coli RNase P. We performed kinetic footprinting studies of M1 RNA folding with the chemical footprinting reagent peroxynitrous acid to provide a detailed description of the folding pathway of RNase P RNA. Our results indicate that, in contrast to the Group I ribozyme, the M1 RNA folds into its catalytically active structure through the formation of two separately folded domains and that the folding of each proceeds through a discrete series of intermediates. Similar rates of folding were observed for regions believed to form the interface between the two domains. This observation is consistent with a kinetic trap which occurs by interaction of the domains during folding.

Methyl-RNA: an evolutionary bridge between RNA and DNA?

Given the apparent limitation of double-stranded RNA (dsRNA) genomes to about 30 kb, together with the complexity of DNA synthesis, it appears difficult for a dsRNA genome to encode all the information required before the transition from an RNA to a DNA genome. Ribonucleotide reductase itself, which synthesizes deoxyribonucleotides from ribonucleotides, requires complex protein radical chemistry, and RNA world genomes may have reached their limits of coding capacity well before such complex enzymes had evolved. The transition from RNA to DNA thus appears to require intermediate steps, and we suggest that the naturally occurring 2'-O-methylated RNA, with chemical properties intermediate between RNA and DNA, is a suitable candidate.

A detailed interpretation of OH radical footprints in a TBP-DNA complex reveals the role of dynamics in the mechanism of sequence-specific binding

The hydroxyl radical footprint of the TATA-binding protein (TBP) bound to the high-affinity sequence TATAAAAG of the adenovirus 2 major late promoter has been quantitatively compared to a 2 ns molecular dynamics simulation of the complex in aqueous solution at room temperature using the CHARMM23 potential. The nucleotide-by-nucleotide analysis of the TBP-TATA hydroxyl radical footprint correlates with the solvent-accessible surface calculated from the dynamics simulation. The results suggest that local reactivity towards OH radicals results from the interplay between the local DNA geometry imposed by TBP binding, and the dynamics of the side-chains contacting the sugar hydrogen atoms. Analysis of the dynamics suggests that, over time, TBP forms stable interactions with the sugar-phosphate backbone through multiple contacts to different partners. This mechanism results in an enthalpic advantage to complex formation at a low entropic cost.

Multiple folding pathways for the P4-P6 RNA domain

We recently described site-specific pyrene labeling of RNA to monitor Mg(2+)-dependent equilibrium formation of tertiary structure. Here we extend these studies to follow the folding kinetics of the 160-nucleotide P4-P6 domain of the Tetrahymena group I intron RNA, using stopped-flow fluorescence with approximately 1 ms time resolution. Pyrene-labeled P4-P6 was prepared using a new phosphoramidite that allows high-yield automated synthesis of oligoribonucleotides with pyrene incorporated at a specific 2'-amino-2'-deoxyuridine residue. P4-P6 forms its higher-order tertiary structure rapidly, with k(obs) = 15-31 s(-1) (t(1/2) approximately 20-50 ms) at 35 degrees C and [Mg(2+)] approximately 10 mM in Tris-borate (TB) buffer. The folding rate increases strongly with temperature from 4 to 45 degrees C, demonstrating a large activation enthalpy DeltaH(double dagger) approximately 26 kcal/mol; the activation entropy DeltaS(double dagger) is large and positive. In low ionic strength 10 mM sodium cacodylate buffer at 35 degrees C, a slow (t(1/2) approximately 1 s) folding component is also observed. The folding kinetics are both ionic strength- and temperature-dependent; the slow phase vanishes upon increasing [Na(+)] in the cacodylate buffer, and the kinetics switch completely from fast at 30 degrees C to slow at 40 degrees C. Using synchrotron hydroxyl radical footprinting, we confirm that fluorescence monitors the same kinetic events as hydroxyl radical cleavage, and we show that the previously reported slow P4-P6 folding kinetics apply only to low ionic strength conditions. One model to explain the fast and slow folding kinetics postulates that some tertiary interactions are present even without Mg(2+) in the initial state. The fast kinetic phase reflects folding that is facilitated by these interactions, whereas the slow kinetics are observed when these interactions are disrupted at lower ionic strength and higher temperature.

Metal-phosphate interactions in the hammerhead ribozyme observed by 31P NMR and phosphorothioate substitutions

The hammerhead ribozyme is a catalytic RNA that requires divalent metal cations for activity under moderate ionic strength. Two important sites that are proposed to bind metal ions in the hammerhead ribozyme are the A9/G10.1 site, located at the junction between stem II and the conserved core, and the scissile phosphate (P1.1). (31)P NMR spectroscopy in conjunction with phosphorothioate substitutions is used in this study to investigate these putative metal sites. The (31)P NMR feature of a phosphorothioate appears in a unique spectral window and can be monitored for changes upon addition of metals. Addition of 1-2 equiv of Cd(2+) to the hammerhead with an A9-S(Rp) or A9-S(S)(Rp) substitution results in a 2-3 ppm upfield shift of the (31)P NMR resonance. In contrast, the P1.1-S(Rp) and P1.1-S(Sp) (31)P NMR features shift slightly and in opposite directions, with a total change in delta of </=0.6 ppm with addition of up to 10 equiv of Cd(2+). No significant shifts are observed for an RNA.RNA duplex with a single, internal phosphorothioate modification upon addition of Cd(2+). Data obtained using model compounds including diethyl phosphate/thiophosphate, AMP, and AMPS, show that a Cd(2+)-S interaction yields an upfield shift for the (31)P NMR resonance, even in the case of a weak coordination such as with diethyl thiophosphate. Taken together, these data predict that Cd(2+) has a high affinity for the A9 site and suggest that there is flexibility in metal coordination within the binding pocket. Cd(2+) interactions with the cleavage site P1.1-S positions are weaker and appear to be stereospecific. These data have implications for mechanisms that have been proposed to explain the influence of metal ions on hammerhead ribozyme activity. These experiments also show the potential utility of (31)P NMR spectroscopy in conjunction with phosphorothioates as a probe for metal binding sites in nucleic acids.