Brownian-dynamics simulations of metal-ion binding to four-way junctions

Geometric deep learning for the prediction of magnesium-binding sites in RNA structures

Magnesium ions (Mg2+) are essential for the folding, functional expression, and structural stability of RNA molecules. However, predicting Mg2+-binding sites in RNA molecules based solely on RNA structures is still challenging. The molecular surface, characterized by a continuous shape with geometric and chemical properties, is important for RNA modelling and carries essential information for understanding the interactions between RNAs and Mg2+ ions. Here, we propose an approach named RNA-magnesium ion surface interaction fingerprinting (RMSIF), a geometric deep learning-based conceptual framework to predict magnesium ion binding sites in RNA structures. To evaluate the performance of RMSIF, we systematically enumerated decoy Mg2+ ions across a full-space grid within the range of 2 to 10 Å from the RNA molecule and made predictions accordingly. Visualization techniques were used to validate the prediction results and calculate success rates. Comparative assessments against state-of-the-art methods like MetalionRNA, MgNet, and Metal3DRNA revealed that RMSIF achieved superior success rates and accuracy in predicting Mg2+-binding sites. Additionally, in terms of the spatial distribution of Mg2+ ions within the RNA structures, a majority were situated in the deep grooves, while a minority occupied the shallow grooves. Collectively, the conceptual framework developed in this study holds promise for advancing insights into drug design, RNA co-transcriptional folding, and structure prediction.

Molecular Mechanisms of DNA Replication and Repair Machinery: Insights from Microscopic Simulations

Reproduction, the hallmark of biological activity, requires making an accurate copy of the genetic material to allow the progeny to inherit parental traits. In all living cells, the process of DNA replication is carried out by a concerted action of multiple protein species forming a loose protein-nucleic acid complex, the replisome. Proofreading and error correction generally accompany replication but also occur independently, safeguarding genetic information through all phases of the cell cycle. Advances in biochemical characterization of intracellular processes, proteomics and the advent of single-molecule biophysics have brought about a treasure trove of information awaiting to be assembled into an accurate mechanistic model of the DNA replication process. In this review, we describe recent efforts to model elements of DNA replication and repair processes using computer simulations, an approach that has gained immense popularity in many areas of molecular biophysics but has yet to become mainstream in the DNA metabolism community. We highlight the use of diverse computational methods to address specific problems of the fields and discuss unexplored possibilities that lie ahead for the computational approaches in these areas.

Holliday Junction Thermodynamics and Structure: Coarse-Grained Simulations and Experiments

Holliday junctions play a central role in genetic recombination, DNA repair and other cellular processes. We combine simulations and experiments to evaluate the ability of the 3SPN.2 model, a coarse-grained representation designed to mimic B-DNA, to predict the properties of DNA Holliday junctions. The model reproduces many experimentally determined aspects of junction structure and stability, including the temperature dependence of melting on salt concentration, the bias between open and stacked conformations, the relative populations of conformers at high salt concentration, and the inter-duplex angle (IDA) between arms. We also obtain a close correspondence between the junction structure evaluated by all-atom and coarse-grained simulations. We predict that, for salt concentrations at physiological and higher levels, the populations of the stacked conformers are independent of salt concentration, and directly observe proposed tetrahedral intermediate sub-states implicated in conformational transitions. Our findings demonstrate that the 3SPN.2 model captures junction properties that are inaccessible to all-atom studies, opening the possibility to simulate complex aspects of junction behavior.

DNA self-assembly: from chirality to evolution

Transient or long-term DNA self-assembly participates in essential genetic functions. The present review focuses on tight DNA-DNA interactions that have recently been found to play important roles in both controlling DNA higher-order structures and their topology. Due to their chirality, double helices are tightly packed into stable right-handed crossovers. Simple packing rules that are imposed by DNA geometry and sequence dictate the overall architecture of higher order DNA structures. Close DNA-DNA interactions also provide the missing link between local interactions and DNA topology, thus explaining how type II DNA topoisomerases may sense locally the global topology. Finally this paper proposes that through its influence on DNA self-assembled structures, DNA chirality played a critical role during the early steps of evolution.

DNA-directed base pair opening

Strand separation is a fundamental molecular process essential for the reading of the genetic information during DNA replication, transcription and recombination. However, DNA melting in physiological conditions in which the double helix is expected to be stable represents a challenging problem. Current models propose that negative supercoiling destabilizes the double helix and promotes the spontaneous, sequence-dependent DNA melting. The present review examines an alternative view and reveals how DNA compaction may trigger the sequence dependent opening of the base pairs. This analysis shows that in DNA crystals, tight DNA-DNA interactions destabilize the double helices at various degrees, from the alteration of the base-stacking to the opening of the base-pairs. The electrostatic repulsion generated by the DNA close approach of the negatively charged sugar phosphate backbones may therefore provide a potential source of the energy required for DNA melting. These observations suggest a new molecular mechanism for the initial steps of strand separation in which the coupling of the DNA tertiary and secondary interactions both actively triggers the base pair opening and stabilizes the intermediate states during the melting pathway.

Molecular dynamics of a DNA Holliday junction: the inverted repeat sequence d(CCGGTACCGG)₄

All-atom molecular dynamics (MD) computer simulations have been applied successfully to duplex DNA structures in solution for some years and found to give close accord with observed results. However, the MD force fields have generally not been parameterized against unusual DNA structures, and their use to obtain dynamical models for this class of systems needs to be independently validated. The four-way junction (4WJ), or Holliday junction, is a dynamic DNA structure involved in central cellular processes of homologous replication and double strand break repair. Two conformations are observed in solution: a planar open-X form (OPN) with a mobile center and four duplex arms, and an immobile stacked-X (STX) form with two continuous strands and two crossover strands, stabilized by high salt conditions. To characterize the accuracy of MD modeling on 4WJ, we report a set of explicit solvent MD simulations of ∼100 ns on the repeat sequence d(CCGGTACCGG)(4) starting from the STX structure (PDB code 1dcw), and an OPN structure built for the same sequence. All 4WJ MD simulations converged to a stable STX structure in close accord with the crystal structure. Our MD beginning in the OPN form converts to the STX form spontaneously at both high and low salt conditions, providing a model for the conformational transition. Thus, these simulations provide a successful account of the dynamical structure of the STX form of d(CCGGTACCGG)(4) in solution, and provide new, to our knowledge, information on the conformational stability of the junction and distribution of counterions in the junction interior.

Structure of d(CGGGTACCCG)4 as a four-way Holliday junction

The crystal structure of the decamer sequence d(CGGGTACCCG)(4) as a four-way Holliday junction has been determined at 2.35 Å resolution. The sequence was designed in order to understand the principles that govern the relationship between sequence and branching structure. It crystallized as a four-way junction structure with an overall geometry similar to those of previously determined Holliday junction structures.

Dissecting electrostatic screening, specific ion binding, and ligand binding in an energetic model for glycine riboswitch folding

Riboswitches are gene-regulating RNAs that are usually found in the 5'-untranslated regions of messenger RNA. As the sugar-phosphate backbone of RNA is highly negatively charged, the folding and ligand-binding interactions of riboswitches are strongly dependent on the presence of cations. Using small angle X-ray scattering (SAXS) and hydroxyl radical footprinting, we examined the cation dependence of the different folding stages of the glycine-binding riboswitch from Vibrio cholerae. We found that the partial folding of the tandem aptamer of this riboswitch in the absence of glycine is supported by all tested mono- and divalent ions, suggesting that this transition is mediated by nonspecific electrostatic screening. Poisson-Boltzmann calculations using SAXS-derived low-resolution structural models allowed us to perform an energetic dissection of this process. The results showed that a model with a constant favorable contribution to folding that is opposed by an unfavorable electrostatic term that varies with ion concentration and valency provides a reasonable quantitative description of the observed folding behavior. Glycine binding, on the other hand, requires specific divalent ions binding based on the observation that Mg(2+), Ca(2+), and Mn(2+) facilitated glycine binding, whereas other divalent cations did not. The results provide a case study of how ion-dependent electrostatic relaxation, specific ion binding, and ligand binding can be coupled to shape the energetic landscape of a riboswitch and can begin to be quantitatively dissected.

Differential stability of DNA crossovers in solution mediated by divalent cations

The assembly of DNA duplexes into higher-order structures plays a major role in many vital cellular functions such as recombination, chromatin packaging and gene regulation. However, little is currently known about the molecular structure and stability of direct DNA–DNA interactions that are required for such functions. In nature, DNA helices minimize electrostatic repulsion between double helices in several ways. Within crystals, B-DNA forms either right-handed crossovers by groove–backbone interaction or left-handed crossovers by groove–groove juxtaposition. We evaluated the stability of such crossovers at various ionic concentrations using large-scale atomistic molecular dynamics simulations. Our results show that right-handed DNA crossovers are thermodynamically stable in solution in the presence of divalent cations. Attractive forces at short-range stabilize such crossover structures with inter-axial separation of helices less than 20 Å. Right-handed crossovers, however, dissociate swiftly in the presence of monovalent ions only. Surprisingly, left-handed crossovers, assembled by sequence-independent juxtaposition of the helices, appear unstable even at the highest concentration of Mg²⁺studied here. Our study provides new molecular insights into chiral association of DNA duplexes and highlights the unique role divalent cations play in differential stabilization of crossover structures. These results may serve as a rational basis to understand the role DNA crossovers play in biological processes.

Preparation of selective and segmentally labeled single-stranded DNA for NMR by self-primed PCR and asymmetrical endonuclease double digestion

We demonstrate a new, efficient and easy-to-use method for enzymatic synthesis of (stereo-)specific and segmental (13)C/(15)N/(2)H isotope-labeled single-stranded DNA in amounts sufficient for NMR, based on the highly efficient self-primed PCR. To achieve this, new approaches are introduced and combined. (i) Asymmetric endonuclease double digestion of tandem-repeated PCR product. (ii) T4 DNA ligase mediated ligation of two ssDNA segments. (iii) In vitro dNTP synthesis, consisting of in vitro rNTP synthesis followed by enzymatic stereo-selective reduction of the C2' of the rNTP, and a one-pot add-up synthesis of dTTP from dUTP. The method is demonstrated on two ssDNAs: (i) a 36-nt three-way junction, selectively (13)C(9)/(15)N(3)/(2)H((1',2'',3',4',5',5''))-dC labeled and (ii) a 39-nt triple-repeat three-way junction, selectively (13)C(9)/(15)N(3)/(2)H((1',2'',3',4',5',5''))-dC and (13)C(9)/(15)N(2)/(2)H((1',2'',3',4',5',5''))-dT labeled in segment C20-C39. Their NMR spectra show the spectral simplification, while the stereo-selective (2)H-labeling in the deoxyribose of the dC-residues, straightforwardly provided assignment of their C1'-H2' and C2'-H2' resonances. The labeling protocols can be extended to larger ssDNA molecules and to more than two segments.

The stability and characteristics of a DNA Holliday junction switch

A Holliday junction (HJ) consists of four DNA double helices, with a branch point discontinuity at the intersection of the component strands. At low ionic strength, the HJ adopts an open conformation, with four widely spaced arms, primarily due to strong electrostatic repulsion between the phosphate groups on the backbones. At high ionic strength, screening of this repulsion induces a switch to a more compact (closed) junction conformation. Fluorescent labelling with dyes placed on the HJ arms allows this conformational switch to be detected optically using fluorescence resonance energy transfer (FRET), producing a sensitive fluorescent output of the switch state. This paper presents a systematic and quantitative survey of the switch characteristics of such a labelled HJ. A short HJ (arm length 8 bp) is shown to be prone to dissociation at low switching ion concentration, whereas an HJ of arm length 12 bp is shown to be stable over all switching ion concentrations studied. The switching characteristics of this HJ have been systematically and quantitatively studied for a variety of switching ions, by measuring the required ion concentration, the sharpness of the switching transition and the fluorescent output intensity of the open and closed states. This stable HJ is shown to have favourable switch characteristics for a number of inorganic switching ions, making it a promising candidate for use in nanoscale biomolecular switch devices.

Assessment of adenyl residue reactivity within model nucleic acids by surface enhanced Raman spectroscopy

We rank the reactivity of the adenyl residues (A) of model DNA and RNA molecules with electropositive subnano size [Ag]n+ sites as a function of nucleic acid primary sequences and secondary structures and in the presence of biological amounts of Cl- and Na+ or Mg2+ ions. In these conditions A is markedly more reactive than any other nucleic acid bases. A reactivity is higher in ribo (r) than in deoxyribo (d) species [pA>pdA and (pA)n>>(pdA)n]. Base pairing decreases A reactivity in corresponding duplexes but much less in r than in d. In linear single and paired dCAG or dGAC loci, base stacking inhibits A reactivity even if A is bulged or mispaired (A.A). dA tracts are highly reactive only when dilution prevents self-association and duplex structures. In d hairpins the solvent-exposed A residues are reactive in CAG and GAC triloops and even more in ATC loops. Among the eight rG1N2R3A4 loops, those bearing a single A (A4) are the least reactive. The solvent-exposed A2 is reactive, but synergistic structural transitions make the initially stacked A residues of any rGNAA loop much more reactive. Mg2+ cross-bridging single strands via phosphates may screen A reactivity. In contrast d duplexes cross-bridging enables "A flipping" much more in rA.U pairs than in dA.T. Mg2+ promotes A reactivity in unpaired strands. For hairpins Mg2+ binding stabilizes the stems, but according to A position in the loops, A reactivity may be abolished, reduced, or enhanced. It is emphasized that not only accessibility but also local flexibility, concerted docking, and cation and anion binding control A reactivity.

Conformational model of the Holliday junction transition deduced from molecular dynamics simulations

Homologous recombination plays a key role in the restart of stalled replication forks and in the generation of genetic diversity. During this process, two homologous DNA molecules undergo strand exchange to form a four-way DNA (Holliday) junction. In the presence of metal ions, the Holliday junction folds into the stacked-X structure that has two alternative conformers. Experiments have revealed the spontaneous transitions between these conformers, but their detailed pathways are not known. Here, we report a series of molecular dynamics simulations of the Holliday junction at physiological and elevated (400 K) temperatures. The simulations reveal new tetrahedral intermediates and suggest a schematic framework for conformer transitions. The tetrahedral intermediates bear resemblance to the junction conformation in complex with a junction-resolving enzyme, T7 endonuclease I, and indeed, one intermediate forms a stable complex with the enzyme as demonstrated in one simulation. We also describe free energy minima for various states of the Holliday junction system, which arise during conformer transitions. The results show that magnesium ions stabilize the stacked-X form and destabilize the open and tetrahedral intermediates. Overall, our study provides a detailed dynamic model of the Holliday junction undergoing a conformer transition.

Electrostatic Interactions and the Folding of the Four-way DNA Junction: Analysis by Selective Methyl Phosphonate Substitution

The structure and dynamics of the four-way (Holliday) junction are strongly dependent on the presence of metal ions. In this study, the importance of phosphate charge in and around the point of strand exchange has been explored by selective replacement with electrically neutral methyl phosphonate groups, guided by crystal structures of the junction in the folded, stacked X conformation. Junction conformation has been analysed by comparative gel electrophoresis and fluorescence resonance energy transfer (FRET). Three of sets of phosphate groups on the exchanging strands have been analysed; those at the point of strand exchange and those to their 3' and 5' sides. The exchanging and 3' phosphate groups form a box of negatively charged groups on the minor groove face of the junction, while the 5' phosphate groups face each other on the major groove side, with their proR oxygen atoms directed at one another. The largest effects are observed on substitution of the exchanging phosphate groups; replacement of both groups leads to the loss of the requirement for addition of metal ions to allow junction folding. When the equivalent phosphate groups on the continuous strands were substituted, a proportion of the junction folded into the alternative conformer so as to bring these phosphate groups onto the exchanging strands. These species did not interconvert, and thus this is likely to result from the alternative diasteromeric forms of the methyl phosphonate group. This shows that some of the conformational effects result from more than purely electrostatic interactions. Smaller but significant effects were observed on substitution of the flanking phosphate groups. All methyl phosphonate substitutions at these positions allowed folding to proceed at a reduced concentration of magnesium ions, with double substitutions more effective than single substitutions. Substitution of 5' phosphates resulted in a greater degree of folding at a given ionic concentration compared to the corresponding 3' phosphate substitutions. These results show that the phosphate groups at the point of strand exchange exert the largest electrostatic effect on junction folding, but a number of phosphate groups in the vicinity of the exchange region contribute to the overall effects.

Global structure of a DNA three-way junction by solution NMR: towards prediction of 3H fold

Three-way junctions (3H) are the simplest and most commonly occurring branched nucleic acids. They consist of three double helical arms (A to C), connected at the junction point, with or without a number of unpaired bases in one or more of the three different strands. Three-way junctions with two unpaired bases in one strand (3HS2) have a high tendency to adopt either of two alternative stacked conformations in which two of the three arms A, B and C are coaxially stacked, i.e. A/B-stacked or A/C-stacked. Empirical stacking rules, which successfully predict for DNA 3HS2 A/B-stacking preference from sequence, have been extended to A/C-stacked conformations. Three novel DNA 3HS2 sequences were designed to test the validity of these extended stacking rules and their conformational behavior was studied by solution NMR. All three show the predicted A/C-stacking preference even in the absence of multivalent cations. The stacking preference for both classes of DNA 3HS2 can thus be predicted from sequence. The high-resolution NMR solution structure for one of the stacked 3HS2 is also reported. It shows a well-defined local and global structure defined by an extensive set of classical NMR restraints and residual dipolar couplings. Analysis of its global conformation and that of other representatives of the 3H family, shows that the relative orientations of the stacked and non-stacked arms, are restricted to narrow regions of conformational space, which can be understood from geometric considerations. Together, these findings open up the possibility of full prediction of 3HS2 conformation (stacking and global fold) directly from sequence.

Microenvironment analysis and identification of magnesium binding sites in RNA

Interactions with magnesium (Mg2+) ions are essential for RNA folding and function. The locations and function of bound Mg2+ ions are difficult to characterize both experimentally and computationally. In particular, the P456 domain of the Tetrahymena thermophila group I intron, and a 58 nt 23s rRNA from Escherichia coli have been important systems for studying the role of Mg2+ binding in RNA, but characteristics of all the binding sites remain unclear. We therefore investigated the Mg2+ binding capabilities of these RNA systems using a computational approach to identify and further characterize their Mg2+ binding sites. The approach is based on the FEATURE algorithm, reported previously for microenvironment analysis of protein functional sites. We have determined novel physicochemical descriptions of site-bound and diffusely bound Mg2+ ions in RNA that are useful for prediction. Electrostatic calculations using the Non-Linear Poisson Boltzmann (NLPB) equation provided further evidence for the locations of site-bound ions. We confirmed the locations of experimentally determined sites and further differentiated between classes of ion binding. We also identified potentially important, high scoring sites in the group I intron that are not currently annotated as Mg2+ binding sites. We note their potential function and believe they deserve experimental follow-up.

RNA folding: models and perspectives

Intrinsic events during RNA folding include conformational search and metal ion binding. Several experimentally testable models have been proposed to explain how large ribozymes accomplish folding. Future challenges include the validation of these models, and the correlation of experimental results and theoretical simulations.

Type III intermediate filament proteins interact with four-way junction DNA and facilitate its cleavage by the junction-resolving enzyme T7 endonuclease I

The isolation from proliferating mouse and human embryo fibroblasts of SDS-stable crosslinkage products of vimentin with DNA fragments containing inverted repeats capable of cruciform formation under superhelical stress and the competitive effect of a synthetic Holliday junction on the binding of cytoplasmic intermediate filament (cIF) proteins to supercoiled DNA prompted a detailed investigation of the proteins' capacity to associate with four-way junction DNA and to influence its processing by junction-resolving endonucleases. Electrophoretic mobility shift analysis of reaction products obtained from vimentin and Holliday junctions under varying ionic conditions revealed efficient complex formation of the filament protein not only with the unstacked, square-planar configuration of the junctions but also with their coaxially stacked X-conformation. Glial fibrillary acidic protein (GFAP) was less efficient and desmin virtually inactive in complex formation. Electron microscopy showed binding of vimentin tetramers or octamers almost exclusively to the branchpoint of the Holliday junctions under physiological ionic conditions. Even at several hundredfold molar excess, sequence-related single- and double-stranded DNAs were unable to chase Holliday junctions from their complexes with vimentin. Vimentin also stimulated bacteriophage T7 endonuclease I in introducing single-strand cuts diametrically across the branchpoint and thus in the resolution of the Holliday junctions. This effect is very likely due to vimentin-induced structural distortion of the branchpoint, as suggested by the results of hydroxyl radical footprinting of Holliday junctions in the absence and the presence of vimentin. Moreover, vimentin, and to a lesser extent GFAP and desmin, interacted with the cruciform structures of inverted repeats inserted into a supercoiled vector plasmid, thereby changing their configuration via branch migration and sensibilizing them to processing by T7 endonuclease I. This refers to both plasmid relaxation caused by unilateral scission and, particularly, linearization via bilateral scission at primary and cIF protein-induced secondary cruciform branchpoints that were identified by T7 endonuclease I footprinting. cIF proteins share these activities with a variety of other architectural proteins interacting with and structurally modulating four-way DNA junctions. In view of the known and hypothetical functions of four-way DNA junctions and associated protein factors in DNA metabolism, cIF proteins as complementary nuclear matrix proteins may play important roles in such nuclear matrix-associated processes as DNA replication, recombination, repair, and transcription, with special emphasis on both the preservation and evolution of the genome.

Conformational and hydration effects of site-selective sodium, calcium and strontium ion binding to the DNA Holliday junction structure d(TCGGTACCGA)(4)

The role of metal ions in determining the solution conformation of the Holliday junction is well established, but to date the picture of metal ion binding from structural studies of the four-way DNA junction is very incomplete. Here we present two refined structures of the Holliday junction formed by the sequence d(TCGGTACCGA) in the presence of Na(+) and Ca(2+), and separately with Sr(2+) to resolutions of 1.85A and 1.65A, respectively. This sequence includes the ACC core found to promote spontaneous junction formation, but its structure has not previously been reported. Almost complete hydration spheres can be defined for each metal cation. The Na(+) sites, the most convincing observation of such sites in junctions to date, are one on either face of the junction crossover region, and stabilise the ordered hydration inside the junction arms. The four Ca(2+) sites in the same structure are at the CG/CG steps in the minor groove. The Sr(2+) ions occupy the TC/AG, GG/CC, and TA/TA sites in the minor groove, giving ten positions forming two spines of ions, spiralling through the minor grooves within each arm of the stacked-X structure. The two structures were solved in the two different C2 lattices previously observed, with the Sr(2+) derivative crystallising in the more highly symmetrical form with two-fold symmetry at its centre. Both structures show an opening of the minor groove face of the junction of 8.4 degrees in the Ca(2+) and Na(+) containing structure, and 13.4 degrees in the Sr(2+) containing structure. The crossover angles at the junction are 39.3 degrees and 43.3 degrees, respectively. In addition to this, a relative shift in the base pair stack alignment of the arms of 2.3A is observed for the Sr(2+) containing structure only. Overall these results provide an insight into the so-far elusive stabilising ion structure for the DNA Holliday junction.

Sequence-dependent folding of DNA three-way junctions

Three-way DNA junctions can adopt several different conformers, which differ in the coaxial stacking of the arms. These structural variants are often dominated by one conformer, which is determined by the DNA sequence. In this study we have compared several three-way DNA junctions in order to assess how the arrangement of bases around the branch point affects the conformer distribution. The results show that rearranging the different arms, while retaining their base sequences, can affect the conformer distribution. In some instances this generates a structure that appears to contain parallel coaxially stacked helices rather than the usual anti-parallel arrangement. Although the conformer equilibrium can be affected by the order of purines and pyrimidines around the branch point, this is not sufficient to predict the conformer distribution. We find that the folding of three-way junctions can be separated into two groups of dinucleotide steps. These two groups show distinctive stacking properties in B-DNA, suggesting there is a correlation between B-DNA stacking and coaxial stacking in DNA junctions.

The bulge region of HIV-1 TAR RNA binds metal ions in solution

Binding of Mg2+, Ca2+ and Co(NH3)6(3+) ions to the HIV-1 TAR RNA in solution was analysed by 19F NMR spectroscopy, metal ion-induced RNA cleavages and Brownian dynamics (BD) simulations. Chemically synthesised 29mer oligoribonucleotides of the TAR sequence labelled with 5-fluorouridine (FU) were used for 19F NMR-monitored metal ion titration. The chemical shift changes of fluorine resonances FU-23, FU-25 and FU-40 upon titration with Mg2+ and Ca2+ ions indicated specific, although weak, binding at the bulge region with the dissociation constants (K(d)) of 0.9 +/- 0.6 and 2.7 +/- 1.7 mM, respectively. Argininamide, inducing largest (19)F chemical shifts changes at FU-23, was used as a reference ligand (K(d) = 0.3 +/- 0.1 mM). In the Pb2+-induced TAR RNA cleavage experiment, strong and selective cleavage of the C24-U25 phosphodiester bond was observed, while Mg2+ and Ca2+ induced cuts at all 3-nt residues of the bulge. The inhibition of Pb2+-specific TAR cleavage by di- and trivalent metal ions revealed a binding specificity [in the order Co(NH3)6(3+) > Mg2+ > Ca2+] at the bulge site. A BD simulation search of potential magnesium ion sites within the NMR structure of HIV-1 TAR RNA was conducted on a set of 20 conformers (PDB code 1ANR). For most cases, the bulge region was targeted by magnesium cations.