Adenine Bulge Research Articles

Structure of a (3+1) hybrid G-quadruplex in the PARP1 promoter.

Poly (ADP-ribose) polymerase 1 (PARP1) has emerged as an attractive target for cancer therapy due to its key role in DNA repair processes. Inhibition of PARP1 in BRCA-mutated cancers has been observed to be clinically beneficial. Recent genome-mapping experiments have identified a non-canonical G-quadruplex-forming sequence containing bulges within the PARP1 promoter. Structural features, like bulges, provide opportunities for selective chemical targeting of the non-canonical G-quadruplex structure within the PARP1 promoter, which could serve as an alternative therapeutic approach for the regulation of PARP1 expression. Here we report the G-quadruplex structure formed by a 23-nucleotide G-rich sequence in the PARP1 promoter. Our study revealed a three-layered intramolecular (3+1) hybrid G-quadruplex scaffold, in which three strands are oriented in one direction and the fourth in the opposite direction. This structure exhibits unique structural features such as an adenine bulge and a G·G·T base triple capping structure formed between the central edgewise loop, propeller loop and 5′ flanking terminal. Given the highly important role of PARP1 in DNA repair and cancer intervention, this structure presents an attractive opportunity to explore the therapeutic potential of PARP1 inhibition via G-quadruplex DNA targeting.

Structural features of a 3' splice site in influenza a.

Influenza A is an RNA virus with a genome of eight negative sense segments. Segment 7 mRNA contains a 3′ splice site for alternative splicing to encode the essential M2 protein. On the basis of sequence alignment and chemical mapping experiments, the secondary structure surrounding the 3′ splice site has an internal loop, adenine bulge, and hairpin loop when it is in the hairpin conformation that exposes the 3′ splice site. We report structural features of a three-dimensional model of the hairpin derived from nuclear magnetic resonance spectra and simulated annealing with restrained molecular dynamics. Additional insight was provided by modeling based on 1H chemical shifts. The internal loop containing the 3′ splice site has a dynamic guanosine and a stable imino (cis Watson–Crick/Watson–Crick) GA pair. The adenine bulge also appears to be dynamic with the A either stacked in the stem or forming a base triple with a Watson–Crick GC pair. The hairpin loop is a GAAA tetraloop closed by an AC pair.

Binding of pixantrone to DNA at CpA dinucleotide sequences and bulge structures.

The binding of the anti-cancer drug pixantrone to three oligonucleotide sequences, d(TCATATGA)2, d(CCGAGAATTCCGG)2 {double bulge = DB} and the non-self complementary d(TACGATGAGTA) : d(TACCATCGTA) {single bulge = SB}, has been studied by NMR spectroscopy and molecular modelling. The upfield shifts observed for the aromatic resonances of pixantrone upon addition of the drug to each oligonucleotide confirmed the drug bound by intercalation. For the duplex sequence d(TCATATGA)2, NOEs were observed from the pixantrone aromatic H7/8 and aliphatic Ha/Hb protons to the H6/H8 and H1' protons of the C2, A3, T6 and G7 nucleotides, demonstrating that pixantrone preferentially binds at the symmetric CpA sites. However, weaker NOEs observed to various protons from the T4 and A5 residues indicated alternative minor binding sites. NOEs from the H7/H8 and Ha/Hb protons to both major (H6/H8) and minor groove (H1') protons indicated approximately equal proportions of intercalation was from the major and minor groove at the CpA sites. Intermolecular NOEs were observed between the H7/H8 and H4 protons of pixantrone and the A4H1' and G3H1' protons of the oligonucleotide that contains two symmetrically related bulge sites (DB), indicative of binding at the adenine bulge sites. For the oligonucleotide that only contains a single bulge site (SB), NOEs were observed from pixantrone protons to the SB G7H1', A8H1' and G9H1' protons, confirming that the drug bound selectively at the adenine bulge site. A molecular model of pixantrone-bound SB could be constructed with the drug bound from the minor groove at the A8pG9 site that was consistent with the observed NMR data. The results demonstrate that pixantrone preferentially intercalates at adenine bulge sites, compared to duplex DNA, and predominantly from the minor groove.

An approach to therapeutic agents through selective targeting of destabilised nucleic acid duplex sequences

The binding of ΔΔ/ΛΛ-[{Ru(phen)(2)}(2)(μ-bb(n))](4+) {where phen = 1,10-phenanthroline, bb(n) = 1,n-bis[4(4'-methyl-2,2'-bipyridyl)]-alkane (ΔΔ/ΛΛ-Rubb(n))} to the non-self complementary oligonucleotide 5'-d(CGCGATAAGCCGC·5'-GCGGCATTACGCG) (3-DB) has been examined using a 4',6-diamidino-2-phenylindole dihydrochloride (DAPI) displacement assay. The 3-DB oligonucleotide contains two single adenine bulge nucleotides that are separated by three base pairs. (1)H NMR spectroscopy data demonstrated that the adenine bases are intra-helical and that the segment containing the two bulge nucleotides and the three A·T base pairs between the bulges forms a destabilised segment within the stable duplex oligonucleotide. The DAPI displacement assay demonstrated that ΔΔ-Rubb(7)-bound 3-DB with higher affinity than the other members of the ΔΔ/ΛΛ-Rubb(n) series. Molecular models suggested that the seven-carbon chain length in ΔΔ-Rubb(7) was ideal to span the distance between the two bulge sites. The binding of ΔΔ-Rubb(7) to 3-DB was also studied by (1)H NMR spectroscopy and molecular modelling. The selective changes in chemical shifts for the resonances from 3-DB upon addition of ΔΔ-Rubb(7) suggested that the metal complex specifically bound at the destabilised segment between A(5) and A(19). Observation in NOESY spectra of NOE cross peaks between 3-DB and ΔΔ-Rubb(7) confirmed that one of the ruthenium centres bound at the A(5) bulge site, with the other metal centre positioned at the A(19) bulge. In addition, ΔΔ-Rubb(7) was found to bind chromosomal DNA extracted from a suspension of Staphylococcus aureus that had been incubated with the ruthenium(ii) complex. As inert dinuclear ruthenium(ii) complexes are capable of being transported into a bacterial cell and bind chromosomal DNA, it is possible that they could be developed into anti-microbial agents that specifically target destabilised segments of DNA that are recognised by essential DNA-binding proteins.

NMR structure of the let-7 miRNA interacting with the site LCS1 of lin-41 mRNA from Caenorhabditis elegans

We have determined the 3D structure of a 34-nt RNA construct, herein named LCS1co, which mimics the interaction of let-7 microRNA (miRNA) to one of its complementary binding sites, LCS1, in the 3′-untranslated region of lin-41 mRNA by solution-state NMR spectroscopy. let-7 miRNAs control the timing of development of the nematode Caenorhabditis elegans and are highly conserved in mammals. The sequence and structure of the two conserved let-7 complementary sites, LCS1 and LCS2, in the 3′-untranslated region of lin-41 mRNA are important for a proper downregulation of lin-41. The high-resolution NMR structure reveals details of the binding of let-7 miRNA to lin-41 mRNA which involves formation of a complex with non-canonical structural elements within the seed region. LCS1co exhibits a stem-loop structure with two stems, an asymmetric internal loop and an adenine bulge. Comparison with the NMR solution-state structure of the let-7:lin-41 complex involving the LCS2-binding site shows that conformational freedom of the asymmetric internal loop of LCS1co correlates with a smaller bend between the upper and lower stems in comparison to the well-defined asymmetric loop of LCS2co.

Binding of a Flexibly-linked Dinuclear Ruthenium(II) Complex to Adenine-bulged DNA Duplexes

Using 1H NMR spectroscopy and molecular modelling, the DNA binding of a chiral dinuclear ruthenium(ii) complex {Δ,Δ-[{Ru(phen)2}2(μ-bb7)]4+; phen = 1,10-phenanthroline, bb7 = 1,7-bis[4(4′-methyl-2,2′-bipyridyl)]-heptane} involving a bridging ligand containing a flexible aliphatic chain has been studied. The binding of the ruthenium(ii) complex was examined with the non-self-complementary duplexes d(CCGAGAATCGGCC):d(GGCCGATTCCGG) (containing a single adenine bulge: designated SB) and d(CCGAGCCGTGCC):d(GGCACGAGCCGG) (containing two adenine bulge sites separated by two base-pairs: designated DB). The NMR data indicated that the ruthenium(ii) complex bound at the bulge site of SB, with one ruthenium centre located at the bulge site with the second metal centre binding with lower affinity and selectivity in the duplex region adjacent to the bulge site. Less specific binding is inferred from chemical shift changes of nucleotide protons two to five base pairs from the single adenine bulge. The ruthenium(ii) complex selectively bound the DB duplex with one metal centre located at each bulge site. The NMR results also suggested that the metal complex binding induced greater changes to the structure of the SB duplex, compared with the DB duplex. Modelling indicates the bridging ligand allowed each ruthenium(ii) metal centre to bind one adenine bulge of the doubly-bulged duplex without disrupting the DNA structure, using the additional torsional flexibility conferred by the aliphatic bridging ligand. However, the second ruthenium(ii) metal centre is not able to bind in the minor groove of the singly-bulged duplex without disrupting the structure, as the metal centre is too bulky. The results of this study suggest dinuclear ruthenium(ii) complexes have considerable potential as probes for DNA and RNA sequences that contain two bulge sites separated by a small number of base-pairs.

Bulged Adenosine Influence on the RNA Duplex Conformation in Solution

The RNA single bulge motif is an unpaired residue within a strand of several complementary base pairs. To gain insight into structural changes induced by the presence of the adenosine bulge on RNA duplex, the solution structures of RNA duplex containing a single adenine bulge (5'-GCAGAAGAGCG-3'/5'-CGCUCUCUGC-3') and a reference duplex with all Watson-Crick base pairs (5'-GCAGAGAGCG-3'/5'-CGCUCUCUGC-3') have been determined by NMR spectroscopy. The reference duplex structure is a regular right-handed helix with all of the attributes of an A-type helix. In the bulged duplex, single adenine bulge stacks into the helix, and the bulge region forms a well-defined structure. Both structures were analyzed by the use of calculated helical parameters. Distortions induced by the accommodation of unpaired residue into the helical structure propagate over the entire structure and are manifested as the reduced base pairs inclination and x-displacement. Intrahelical position of bulged adenine A5 is stabilized by efficient stacking with 5'-neighboring residues G4.

Dinuclear ruthenium(ii) complexes as potential probes for RNA bulge sites

1H NMR spectroscopy and molecular modelling have been used to investigate the binding of the DeltaDelta-and LambdaLambda-enantiomers of the dinuclear ruthenium(II) complex [[Ru(Me2bpy)2]2(mu-bpm)]4+ [Me2bpy = 4,4'-dimethyl-2,2'-bipyridine; bpm = 2,2'-bipyrimidine] to an RNA tridecanucleotide duplex containing a single-base bulge [r(CCGAGAAUUCCGG)2]], and the corresponding control dodecanucleotide [r(CCGGAAUUCCGG)2]. Both enantiomers bound the control RNA sequence weakly. From upfield shifts of the metal complex H3 and H3' protons throughout the titration of the control dodecanucleotide with DeltaDelta-[[Ru(Me2bpy)2]2(mu-bpm)]4+, a binding constant of 1 x 10(3) M(-1) was determined. In NOESY spectra of the control sequence with added DeltaDelta-[[Ru(Me2bpy)2]2(mu-bpm)]4+, NOEs were only observed to protons from the terminal base-pair residues. No significant changes in chemical shift were observed for either the metal complex or RNA protons upon addition of the LambdaLambda-enantiomer to the control dodecanucleotide. The DeltaDelta-[[Ru(Me2bpy)2]2(mu-bpm)]4+ complex bound the bulge-containing RNA with a significantly greater affinity (6 x 10(4) M(-1)) than the non-bulge control RNA duplex. Competition binding experiments indicated that the LambdaLambda-isomer bound the tridecanucleotide with similar affinity to the DeltaDelta-enantiomer. Addition of DeltaDelta-[[Ru(Me2bpy)2]2(mu-bpm)]4+ to the bulge-containing tridecanucleotide induced selective changes in chemical shift for the base H8 and sugar H1' resonances from the adenine bulge residue, and resonances from nucleotide residues adjacent to the bulge site. Intermolecular NOEs observed in NOESY spectra of the tridecanucleotide with added DeltaDelta-[[Ru(Me2bpy)2]2(mu-bpm)]4+ confirmed the selective binding of the ruthenium complex at the bulge site. Preliminary binding models, consistent with the NMR data, showed that the ruthenium complex could effectively associate in the RNA minor groove at the bulge site.

Dinuclear ruthenium(ii) complexes with flexible bridges as non-duplex DNA binding agents

The stereoisomers of a series of dinuclear ruthenium(ii) complexes [{Ru(phen)(2)}(2)(micro-BL)](4+) (phen = 1,10-phenanthroline) with flexible bridging ligands (BL) bb2 {1,2-bis[4(4'-methyl-2,2'-bipyridyl)]ethane}, bb5 {1,5-bis[4(4'-methyl-2,2'-bipyridyl)]pentane}, bb7 {1,7-bis[4(4'-methyl-2,2'-bipyridyl)]heptane}, and bb10 {1,10-bis[4(4'-methyl-2,2'-bipyridyl)]decane} have been synthesised. Their binding to a control dodecanucleotide, d(CCGGAATTCCGG)(2), and a tridecanucleotide, d(CCGAGAATTCCGG)(2), which contains a single adenine bulge have been studied using fluorescence displacement assays involving intercalating and groove-binding dyes, equilibrium dialysis and binding affinity chromatography. The fluorescence intercalator displacement (FID) assay indicated that LambdaLambda-[{Ru(phen)(2)}(2)(micro-bb7)](4+) had the greatest binding affinity with all the oligonucleotides, whereas an analogous fluorescence technique using a minor-groove binding dye, equilibrium dialysis and affinity binding chromatography showed that DeltaDelta-[{Ru(phen)(2)}(2)(micro-bb7)](4+) had the strongest binding. An (1)H NMR study of the binding of the DeltaDelta-enantiomer of [{Ru(phen)(2)}(2)(micro-bb7)](4+) to d(CCGAGAATTCCGG)(2) confirmed the selectivity of the metal complex for the bulge site and provided the basis for an energy-minimised binding model of the dinuclear ruthenium complex with the single adenine bulge containing trinucleotide. The binding model demonstrated the ability of the flexibly-linked complex to follow the curvature of the DNA minor groove.

Meso-[{Ru(phen) 2} 2(μ-bpm)] 4+: A high-affinity DNA bulge probe {bpm = 2,2′-bipyrimidine; phen = 1,10-phenanthroline}

The binding of the stereoisomers of [{Ru(Me 2bpy) 2} 2(μ-bpm)] 4+, [{Ru(phen) 2} 2(μ-bpm)] 4+ and [{Ru(Me 2phen) 2} 2(μ-bpm)] 4+ (Me 2bpy = 4,4′-dimethyl-2,2′-bipyridine; bpm = 2,2′-bipyrimidine; phen = 1,10-phenanthroline; Me 2phen = 4,7-dimethyl-1,10-phenanthroline) to a tridecanucleotide d(CCGAGAATTCCGG) 2 which contains a single adenine bulge site, and four control dodecanucleotides, have been studied using a fluorescence intercalator displacement (FID) assay. The meso isomer of [{Ru(phen) 2} 2(μ-bpm)] 4+ showed the strongest binding to the bulge-containing tridecanucleotide. In order to gain a greater understanding of the basis of the higher affinity exhibited by the meso isomer towards the bulge sequence, a 1H NMR study of the binding of the two enantiomers (ΔΔ and ΛΛ) of rac-[{Ru(phen) 2} 2(μ-bpm)] 4+, and the, meso (ΔΛ) diastereoisomer, to the tridecanucleotide d(CCGAGAATTCCGG) 2 was carried out. The NMR results suggest that the meso isomer binds selectively at the bulge site in the tridecanucleotide minor groove, but closer to the 3′-direction and with less structural perturbations of the groove than the ΔΔ and ΛΛ isomers. The results of this study confirm that dinuclear ruthenium complexes have excellent potential as DNA bulge probes, and meso-[{Ru(phen) 2} 2(μ-bpm)] 4+ in particular has a high affinity (1 × 10 6 M −1) and selectivity for a single adenine bulge site.

Total enantioselectivity in the DNA binding of the dinuclear ruthenium(II) complex [[Ru(Me2bpy)2]2(mu-bpm)]4+ [bpm = 2,2'-bipyrimidine; Me2bpy = 4,4'-dimethyl-2,2'-bipyridine

The binding of the three stereoisomers (DeltaDelta-, LambdaLambda- and DeltaLambda-) of the dinuclear ruthenium(II) complex [[Ru(Me2bpy)2]2(mu-bpm)]4+ [Me2bpy = 4,4'-dimethyl-2,2'-bipyridine; bpm = 2,2'-bipyrimidine] to a tridecanucleotide containing a single adenine bulge has been studied by 1H NMR spectroscopy. The addition of the DeltaDelta-isomer to d(CCGAGAATTCCGG)2 induced significant chemical shift changes for the base and sugar resonances of the residues at the bulge site (G3A4G5/C11C10), whereas small shifts were observed upon addition of the enantiomeric LambdaLambda-form. NOESY spectra of the tridecanucleotide bound with the DeltaDelta-isomer revealed intermolecular NOE's between the metal complex and the nucleotide residues at the bulge site, while only weak NOE's were observed to terminal residues to the LambdaLambda-form. Competitive binding studies were performed where both enantiomers were simultaneously added to the tridecanucleotide, and for all ratios of the two stereoisomers the DeltaDelta-isomer remained selectively bound at the bulge site with the LambdaLambda-enantiomer localised at the terminal regions of the tridecanucleotide. The meso-diastereoisomer (DeltaLambda) was found to bind to the tridecanucleotide with characteristics intermediate between the DeltaDelta- and LambdaLambda-enantiomers of the rac form. Two distinct sets of metal complex resonances were observed, with one set having essentially the same shift as the free metal complex, whilst the other set of resonances exhibited significant shifts. The NOE data indicated that the meso-diastereoisomer does not bind as selectively as the DeltaDelta-isomer, with NOE's observed to a greater number of nucleotide residues compared to the DeltaDelta-form. This study provides a rare example of total enantioselectivity in the binding of an inert transition metal complex to DNA, produced by the shape recognition of both ruthenium(II) centres.

Dinuclear ruthenium(ii) complexes as probes for DNA bulge sites

The binding of the dinuclear ruthenium(II) complexes ΔΔ-[{Ru(Me2bpy)2}2(μ-bpm)]4+ and ΔΔ-[{Ru(bpy)2}(μ-bpm){Ru(Me2bpy)2}]4+ {bpy = 2,2′-bipyridine; Me2bpy = 4,4′-dimethyl-2,2′-bipyridine; bpm = 2,2′-bipyrimidine} to a tridecanucleotide containing a single adenine bulge has been studied by 1H NMR spectroscopy. The addition of either complex to d(CCGAGAATTCCGG)2 induced significant chemical shift changes for the base and sugar resonances of the residues at the bulge site (G3A4G5/C11C10). In NOESY spectra of the tridecanucleotide bound with either ruthenium species, NOEs were observed from the H1′ and H4′ protons of the nucleotide residues at the bulge site to the Me2bpy and bpm protons of the metal complex. These results indicate that the dinuclear ruthenium species selectively bind at the adenine bulge site in the minor groove. A simple model was constructed for the binding of the non-symmetrical ΔΔ-[{Ru(bpy)2}(μ-bpm){Ru(Me2bpy)2}]4+ complex, which shows that the minor groove has significantly widened to allow the relatively bulky complex to bind. As the metal complexes specifically bind at an adenine bulge site in a segment of DNA with an affinity considerably greater than that for standard duplex DNA, the results presented here suggest that non-intercalating dinuclear complexes may be excellent diagnostic agents for DNA bulged sequences.

Conformations of an Adenine Bulge in a DNA Octamer and Its Influence on DNA Structure from Molecular Dynamics Simulations

Molecular dynamics simulations have been applied to the DNA octamer d(GCGCA-GAAC) · d(GTTCGCGC), which has an adenine bulge at the center to determine the pathway for interconversion between the stacked and extended forms. These forms are known to be important in the molecular recognition of bulges. From a total of ∼35 ns of simulation time with the most recent CHARMM27 force field a variety of distinct conformations and subconformations are found. Stacked and fully looped-out forms are in excellent agreement with experimental data from NMR and x-ray crystallography. Furthermore, in a number of conformations the bulge base associates with the minor groove to varying degrees. Transitions between many of the conformations are observed in the simulations and used to propose a complete transition pathway between the stacked and fully extended conformations. The effect on the surrounding DNA sequence is investigated and biological implications of the accessible conformational space and the suggested transition pathway are discussed, in particular for the interaction of the MS2 replicase operator RNA with its coat protein.

Crystal structure of an adenine bulge in the RNA chain of a DNA·RNA hybrid, d(CTCCTCTTC)·r(gaagag[formula omitted]gag)

Crystal structure of a DNA·RNA hybrid, d(CTCCTCTTC)·r(gaagagagag), with an adenine bulge in the polypurine RNA strand was determined at 2.3 Å resolution. The structure was solved by the molecular replacement method and refined to a final R-factor of 19.9 % (Rfree 22.2 %). The hybrid duplex crystallized in the space group I222 with unit cell dimens ions, a = 46.66 Å, b = 47.61 Å and c = 54.05 Å, and adopts the A-form conformation. All RNA and DNA sugars are in the C3′-endo conformation, the glycosyl angles in anti conformation and the majority of the C4′-C5′ torsion angles in g+ except two trans angles, in conformity with the C3′-endo rigid nucleotide hypothesis. The adenine bulge is looped out and it is also in the anti C3′-endo conformation. The bulge is involved in a base-triple (C·g)∗a interaction with the end base-pair (C9·g10) in the minor groove of a symmetry-related molecule. The 2′ hydroxyl group of g15 is hydrogen bonded to O2P and O5′ of g17, skipping the bulged adenine a16 and stabilizing the sugar-phosphate backbone of the hybrid. The hydrogen bonding and the backbone conformation at the bulged adenine site is very similar to that found in the crystal structure of a protein-RNA complex.

Conformational Deformability of RNA: A Harmonic Mode Analysis

The harmonic mode analysis method was used to characterize the conformational deformability of regular Watson-Crick paired, mismatch- and bulge-containing RNA. Good agreement between atomic Debye-Waller factors derived from x-ray crystallography of a regular RNA oligonucleotide and calculated atomic fluctuations was obtained. Calculated helical coordinate fluctuations showed a small sequence dependence of up to ∼30–50%. A negative correlation between motions at a given base pair step and neighboring steps was found for most helical coordinates. Only very few calculated modes contribute significantly to global motions such as bending, twisting, and stretching of the RNA molecules. With respect to a local helical description of the RNA helix our calculations suggest that RNA bending is mostly due to a periodic change in the base pair step descriptors slide and roll. The presence of single guanine:uridine or guanine:adenine mismatches had little influence on the calculated RNA flexibility. In contrast, for tandem guanine:adenine base pairs the harmonic mode approach predicts a significantly reduced conformational flexibility in the case of a sheared arrangement and slightly enhanced flexibility for a face-to-face (imino proton) pairing relative to regular RNA. The presence of a single extra adenine bulge nucleotide stacked between flanking sequences resulted in an increased local atomic mobility around the bulge site (∼40%) and a slightly enhanced global bending flexibility. For an adenine bulge nucleotide in a looped-out conformation a strongly enhanced bulge nucleotide mobility but no increased bending flexibility compared to regular RNA was found.

Conformational analysis of single-base bulges in A-form DNA and RNA using a hierarchical approach and energetic evaluation with a continuum solvent model

The analysis and prediction of non-canonical structural motifs in RNA is of great importance for an understanding of the function and design of RNA structures. A hierarchical method has been employed to generate a large variety of sterically possible conformations for a single-base adenine bulge structure in A-form DNA and RNA. A systematic conformational search was performed on the isolated bulge motif and neighboring nucleotides under the constraint to fit into a continuous helical structure. These substructures were recombined with double-stranded DNA or RNA. Energy minimization resulted in more than 300 distinct bulge conformations. Energetic evaluation using a solvation model based on the finite-difference Poisson-Boltzmann method identified three basic classes of low-energy structures. The three classes correspond to conformations with the bulge base stacked between flanking nucleotides (I), location of the bulge base in the minor groove (II) and conformations with a continuous stacking of the flanking helices and a looped out bulge base (III). For the looped out class, two subtypes (IIIa and IIIb) with different backbone geometries at the bulge site could be distinguished. The conformation of lowest calculated energy was a class I structure with backbone torsion angles close to those in standard A-form RNA. Conformations very close to the extra-helical looped out bulge structure determined by X-ray crystallography were also among the low-energy structures. In addition, topologies observed in other experimentally determined bulge structures have been found among low-energy conformers. The implicit solvent model was further tested by comparing an uridine and adenine bulge flanked by guanine:cytosine base-pairs, respectively. In agreement with the experimental observation, a looped out form was found as the energetically most favorable form for the uridine bulge and a stacked conformation in case of the adenine bulge. The inclusion of solvation effects especially electrostatic reaction field contributions turned out to be critically important in order to select realistic low-energy bulge structures from a large number of sterically possible conformations. The results indicate that the approach might be useful to model the three-dimensional structure of non-canonical motifs embedded in double-stranded RNA, in particular, to restrict the number of possible conformations to a manageable number of conformers with energies below a certain threshold.

Analysis of the stability of looped-out and stacked-in conformations of an adenine bulge in DNA using a continuum model for solvent and ions

A combination of conformational search, energy minimization, and energetic evaluation using a continuum solvent treatment has been employed to study the stability of various conformations of the DNA fragment d(CGCAGAA)/d(TTCGCG) containing a single adenine bulge. The extra-helical (looped-out) bulge conformation derived from a published x-ray structure and intra-helical (stacked bulge base) model structures partially based on nuclear magnetic resonance (NMR) data were used as start structures for the conformational search. Solvent-dependent contributions to the stability of the conformations were calculated from the solvent exposed molecular surface area and by using the finite difference Poisson-Boltzmann approach. Three classes (I-III) of bulge conformations with calculated low energies can be distinguished. The lowest-energy conformations were found in class I, corresponding to structures with the bulge base stacked between flanking helices, and class II, composed of structures forming a triplet of the bulge base and a flanking base pair. All extra-helical bulge structures, forming class III, were found to be less stable compared with the lowest energy structures of class I and II. The results are consistent with NMR data on an adenine bulge in the same sequence context indicating an intra-helical or triplet bulge conformation in solution. Although the total energies and total electrostatic energies of the low-energy conformations show only relatively modest variations, the energetic contributions to the stability were found to vary significantly among the classes of bulge structures. All intra-helical bulge structures are stabilized by a more favorable Coulomb charge-charge interaction but destabilized by a larger electrostatic reaction field contribution compared with all extra-helical and most triplet bulge structures. Van der Waals packing interactions and nonpolar surface-area-dependent contributions appear to favor triplet class II structures and to a lesser degree also the intra-helical stacked bulge conformations. The large conformational variation found for class III conformers might add a favorable entropic contribution to the stability of the extra-helical bulge form.