X-ray crystallographic observation of "in-line" and "adjacent" conformations in a bulged self-cleaving RNA/DNA hybrid

Dioxaphosphorinane-constrained nucleic Acid dinucleotides as tools for structural tuning of nucleic acids

We describe a rational approach devoted to modulate the sugar-phosphate backbone geometry of nucleic acids. Constraints were generated by connecting one oxygen of the phosphate group to a carbon of the sugar moiety. The so-called dioxaphosphorinane rings were introduced at key positions along the sugar-phosphate backbone allowing the control of the six-torsion angles α to ζ defining the polymer structure. The syntheses of all the members of the D-CNA family are described, and we emphasize the effect on secondary structure stabilization of a couple of diastereoisomers of α,β-D-CNA exhibiting wether B-type canonical values or not.

Hydrolysis of bulged nucleotides in hybrids formed by RNA and imidazole-derivatized oligo-2'-O-methylribonucleotides

In order to enhance the efficacy of small antisense molecules, we examined a series of antisense oligonucleotides derivatized with functional groups designed to enable them to hydrolyze their RNA target. Solid phase synthetic methods were used to prepare imidazole-derivatized antisense oligo-2'-O-methylribonucleotides. Upon binding, these oligonucleotides create internal bulged bases in the target RNA that serve as sites for hydrolysis. We observed that an oligonucleotide derivatized with a side chain containing two imidazole groups was capable of hydrolyzing 58% of its RNA target when incubated with the target for 48 hours at 37°C and physiological pH.

Use of nucleic Acid analogs for the study of nucleic Acid interactions

Unnatural nucleosides have been explored to expand the properties and the applications of oligonucleotides. This paper briefly summarizes nucleic acid analogs in which the base is modified or replaced by an unnatural stacking group for the study of nucleic acid interactions. We also describe the nucleoside analogs of a base pair-mimic structure that we have examined. Although the base pair-mimic nucleosides possess a simplified stacking moiety of a phenyl or naphthyl group, they can be used as a structural analog of Watson-Crick base pairs. Remarkably, they can adopt two different conformations responding to their interaction energies, and one of them is the stacking conformation of the nonpolar aromatic group causing the site-selective flipping of the opposite base in a DNA double helix. The base pair-mimic nucleosides can be used to study the mechanism responsible for the base stacking and the flipping of bases out of a nucleic acid duplex.

CHARMM force field parameters for simulation of reactive intermediates in native and thio-substituted ribozymes

Force field parameters specifically optimized for residues important in the study of RNA catalysis are derived from density-functional calculations, in a fashion consistent with the CHARMM27 all-atom empirical force field. Parameters are presented for residues that model reactive RNA intermediates and transition state analogs, thio-substituted phosphates and phosphoranes, and bound Mg(2+) and di-metal bridge complexes. Target data was generated via density-functional calculations at the B3LYP/6-311++G(3df,2p)// B3LYP/6-31++G(d,p) level. Partial atomic charges were initially derived from CHelpG electrostatic potential fitting and subsequently adjusted to be consistent with the CHARMM27 charges. Lennard-Jones parameters were determined to reproduce interaction energies with water molecules. Bond, angle, and torsion parameters were derived from the density-functional calculations and renormalized to maintain compatibility with the existing CHARMM27 parameters for standard residues. The extension of the CHARMM27 force field parameters for the nonstandard biological residues presented here will have considerable use in simulations of ribozymes, including the study of freeze-trapped catalytic intermediates, metal ion binding and occupation, and thio effects.

Base-tetrad swapping results in dimerization of RNA quadruplexes: implications for formation of the i-motif RNA octaplex

Nucleic acids adopt different multistranded helical architectures to perform various biological functions. Here, we report a crystal structure of an RNA quadruplex containing "base-tetrad swapping" and bulged nucleotide at 2.1-Angstroms resolution. The base-tetrad swapping results in a dimer of quadruplexes with an intercalated octaplex fragment at the 5' end junction. The intercalated base tetrads provide the basic repeat unit for constructing a model of intercalated RNA octaplex. The model we obtained shows fundamentally different characteristics from duplex, triplex, and quadruplex. We also observed two different orientations of bulged uridine residues that are related to the interaction with surroundings. This structural evidence reflects the conformational flexibility of bulged nucleotides in RNA quadruplexes and implies the potential roles of bulged nucleotides as recognition and interaction sites in RNA-protein and RNA-RNA interactions.

Conformational transitions in RNA single uridine and adenosine bulge structures: a molecular dynamics free energy simulation study

Extra unmatched nucleotides (single base bulges) are common structural motifs in folded RNA molecules and can participate in RNA-ligand binding and RNA tertiary structure formation. Often these processes are associated with conformational transitions in the bulge region such as flipping out of the bulge base from an intrahelical stacked toward a looped out state. Knowledge of the flexibility of bulge structures and energetics of conformational transitions is an important prerequisite to better understand the function of this RNA motif. Molecular dynamics simulations were performed on single uridine and adenosine bulge nucleotides at the center of eight basepair RNA molecules and indicated larger flexibility of the bulge bases compared to basepaired regions. The umbrella sampling method was applied to study the bulge base looping out process and accompanying conformational and free energy changes. Looping out toward the major groove resulted in partial disruption of adjacent basepairs and was found to be less favorable compared to looping out toward the minor groove. For both uridine and adenosine bulges, a positive free energy change for full looping out was obtained which was approximately 1.5 kcal mol-1 higher in the case of the adenosine compared to the uridine bulge system. The simulations also indicated stable partially looped out states with the bulge bases located in the RNA minor groove and forming base triples with 5'-neighboring basepairs. In the case of the uridine bulge this state was more stable than the intrahelical stacked bulge structure. Induced looping out toward the minor groove involved crossing of an energy barrier of approximately 3.5 kcal mol-1 before reaching the base triple state. A continuum solvent analysis of intermediate bulge states indicated that electrostatic interactions stabilize looped out and base triple states, whereas van der Waals interactions and nonpolar contributions favor the stacked bulge conformation.

Site-Selective RNA Cleavage by DNA Bearing a Base Pair-Mimic Nucleoside

We have synthesized the deoxyadenosine derivative tethering a phenyl group (X), which mimics the Watson-Crick A/T base pair. The RNA/DNA hybrid duplexes containing X in the middle of the DNA sequence showed a similar thermal stability regardless of the ribonucleotide species (A, G, C, or U) opposite to X, probably because of the phenyl group stacking inside of the duplex accompanied by the opposite ribonucleotide base flipped in an extrahelical position. The RNA strand hybridized with the DNA strand bearing X was cleaved on the 3'-side of the ribonucleotide opposite to X in the presence of MgCl2, and the RNA sequence to be cleaved was not restricted. The site-specific RNA hydrolysis suggests that the DNA strand bearing X has the advantage of the site-selective base flipping in the target sequence and the development of a "universal deoxyribozyme" to exclusively cleave a target RNA sequence.

Molecular dynamics simulation study of oriented polyamine- and Na-DNA: sequence specific interactions and effects on DNA structure

Molecular dynamics (MD) computer simulations have been carried out on four systems that correspond to an infinite array of parallel ordered B-DNA, mimicking the state in oriented DNA fibers and also being relevant for crystals of B-DNA oligonucleotides. The systems were all comprised of a periodical hexagonal cell with three identical DNA decamers, 15 water molecules per nucleotide, and counterions balancing the DNA charges. The sequence of the double helical DNA decamer was d(5'-ATGCAGTCAG)xd(5'-TGACTGCATC). The counterions were the two natural polyamines spermidine(3+) (Spd(3+)) and putrescine(2+) (Put(2+)), the synthetic polyamine diaminopropane(2+) (DAP(2+)), and the simple monovalent cation Na(+). This work compares the specific structures of the polyamine- and Na-DNA systems and how they are affected by counterion interactions. It also describes sequence-specific hydration and interaction of the cations with DNA. The local DNA structure is dependent on the nature of the counterion. Even the very similar polyamines, Put(2+) and DAP(2+), show clear differences in binding to DNA and in effect on hydration and local structure. Generally, the polyamines disorder the hydration of the DNA around their binding sites whereas Na(+) being bound to DNA attracts and organizes water in its vicinity. Cation binding at the selected sites in the minor and in the major groove is compared for the different polyamines and Na(+). We conclude that the synthetic polyamine (DAP(2+)) binds specifically to several structural and sequence-specific motifs on B-DNA, unlike the natural polyamines, Spd(3+) and Put(2+). This specificity of DAP(2+) compared to the more dynamic behavior of Spd(3+) and Put(2+) may explain why the latter polyamines are naturally occurring in cells.

Mechanics of DNA flexibility visualized by selective 2'-amine acylation at nucleotide bulges

We used selective acylation of 2'-amine-substituted nucleotides to visualize local backbone conformations that occur preferentially at bulged sites in DNA duplexes. 2'-Amine acylation reports local nucleotide flexibility because unconstrained 2'-amino nucleotides more readily reach a reactive conformation in which the amide-forming transition state is stabilized by interactions between the amine nucleophile and the adjacent 3'-phosphodiester group. Bulged 2'-amine-substituted cytidine nucleotides react approximately 20-fold more rapidly than nucleotides constrained by base-pairing at 35 degrees C. In contrast, base-paired 2'-amine-substituted nucleotides flanked by a 5' or 3' bulge react two- or six-fold more rapidly, respectively, than the perfectly paired duplex. The relative lack of 2'-amine reactivity for nucleotides adjacent to a DNA bulge emphasizes, first, that structural perturbations do not extend significantly into the flanking duplex structure. Second, the exquisite sensitivity towards very local perturbations in nucleic acid structure suggests that 2'-amine acylation can be used to chemically interrogate deletion mutations in DNA. Finally, these data support the mechanical interpretation that the reactive ribose conformation for 2'-amine acylation requires that the base lies out of the helix and in the major groove, a mechanistic insight useful for designing 2'-amine-based sensors.

A molecular dynamics simulation study of oriented DNA with polyamine and sodium counterions: diffusion and averaged binding of water and cations

Four different molecular dynamics (MD) simulations have been performed for ordered DNA decamers, d(5'-ATGCAGTCAG)*d(5'-TGACTGCATC). The counterions were the two natural polyamines spermidine3+ (Spd3+) and putrescine2+ (Put2+), the synthetic polyamine diaminopropane2+ (DAP2+) and Na+. The simulation set-up corresponds to an infinite array of parallel DNA mimicking the state in oriented DNA fibers or crystals. This work describes general properties of polyamine and Na+ binding to DNA. Simulated diffusion coefficients show satisfactory agreement with experimental NMR diffusion data of comparable systems. The interaction of the polyamines with DNA is dynamic in character and the cations mostly form short-lived contacts with the electronegative binding sites of DNA. Polyamines, Na+ and water interact most frequently with the charged phosphate atoms with preference for association from the minor groove side with O1P over O2P. There is a strong anti-correlation in the cation binding to the electronegative groups of DNA, i.e. the presence of a cation near one of the DNA sites repels other cations from binding to this and to the other sites separated by <7.5 A from each other. In contrast to the other polyamines, DAP2+ is able to form 'bridges' connecting neighboring phosphate groups along the DNA strand. A small fraction of DAP2+ and Put2+ can be found in the major grooves, while Spd3+ is absent there. The results of the MD simulations reveal principal differences in the polyamine-DNA interactions between the natural (Spd3+, Put2+ and spermine4+) and synthetic (DAP2+) polyamines.

Molecular modeling and dynamics studies of HIV-1 kissing loop structures

Recognition of an RNA loop by another RNA loop is involved in several biological functions. The dimerization of two copies of the HIV-1 genomic RNA is thought to be involved in several steps of the retroviral life cycle. It has been shown that the dimerization of the two HIV-1 RNA genomes is initiated by the so called kissing loop. The 9nt kissing loop consists of a palindromic 6nt sequence that forms Watson-Crick base-pairs at the kissing site in HIV-1. We report the results of our molecular modeling and dynamics studies on two major subtype isolates (MAL and LAI) of HIV-1 kissing loop structures. From our modeling studies, we conclude that the conformation of the loop in the monomer might be closer to the A-RNA-like conformation in order to form an initial kissing structure. This is achieved by the stacking interactions of the bases at the 3' end of the loop and by the intramolecular tertiary interactions of a single linker nucleotide. We discuss the effect of the loop size and the structural limitations on the formation of kissing loop structures. Also, we propose a possible mechanism to convert the kissing loop structure to a stable extended duplex structure without unwinding the stems.

RNA-cleaving DNA enzymes with altered regio- or enantioselectivity

In vitro evolution methods were used to obtain DNA enzymes that cleave either a 2',5'-phosphodiester following a D-ribonucleotide or a 3',5'-phosphodiester following an L-ribonucleotide. Both enzymes can operate in an intermolecular reaction format with multiple turnover. The DNA enzyme that cleaves a 2',5'-phosphodiester exhibits a k(cat) of approximately 0.01 min(-1) and catalytic efficiency, k(cat)/K(m), of approximately 10(8) M(-1) min(-1). The enzyme that cleaves an L-ribonucleotide is about 10-fold slower and has a catalytic efficiency of approximately 4 x 10(5) M(-1) min(-1). Both enzymes require a divalent metal cation for their activity and have optimal catalytic rate at pH 7-8 and 35-50 degrees C. In a comparison of each enzyme's activity with either its corresponding substrate that contains an unnatural ribonucleotide or a substrate that instead contains a standard ribonucleotide, the 2',5'-phosphodiester-cleaving DNA enzyme exhibited a regioselectivity of 6000-fold, while the L-ribonucleotide-cleaving DNA enzyme exhibited an enantioselectivity of 40-fold. These molecules demonstrate how in vitro evolution can be used to obtain regio- and enantioselective catalysts that exhibit specificities for nonnatural analogues of biological compounds.

Crystal structure of an RNA duplex r(gugucgcac)(2) with uridine bulges

The crystal structure of a nonamer RNA duplex with a uridine bulge in each strand, r(gugucgcac)(2), was determined at 1.4 A resolution. The structure was solved by multiple anomalous diffraction phasing method using a three-wavelength data set collected at the Advanced Protein Source and refined to a final R(work)/R(free) of 21.2 %/23.4 % with 33,271 independent reflections (Friedel pairs unmerged). The RNA duplex crystallized in the tetragonal space group P4(1)22 with two independent molecules in the asymmetric unit. The unit cell dimensions are a=b=47.18 A and c=80.04 A. The helical region of the nonamer adopts the A-form conformation. The uridine bulges assume similar conformations, with uracils flipping out and protruding into the minor groove. The presence of the bulge induces very large twist angles (approximately +50 degrees) between the base-pairs flanking the bulges while causing profound kinks in the helix axis at the bulges. This severe twist and the large kink in turn produces a very narrow major groove at the middle of the molecule. The ribose sugars of the guanosines before the bulges adopt the C2'-endo conformation while the rest, including the bulges, are in the C3'-endo conformation. The intrastrand phosphate-phosphate (P-P) distance of the phosphate groups flanking the bulges (approximately 4.4 A) are significantly shorter than the average P-P distance in the duplex (6.0 A). This short distance between the two phosphate groups brings the non-bridging oxygen atoms close to each other where a calcium ion is bound to each strand. The calcium ions in molecule 1 are well defined while the calcium ions in molecule 2 are disordered.