Bulge Loop Research Articles

Effect of Internal and Bulge Loops on the Thermal Stability of Small DNA Duplexes.

The thermal stabilities of DNA duplexes analogous to the let-7 microRNA:lin-41 mRNA complex from Caenorhabditis elegans have been measured by free solution capillary electrophoresis. DNA duplexes with the same stems but different types of internal or bulge loops and a control with no loop have also been studied. The melting temperatures of the DNA derivatives increased linearly with the logarithm of the Na+ or K+ ion concentration in the solution. Peaks in the electropherograms corresponding to duplexes with internal or bulge loops exhibited extensive tailing at high temperatures, suggesting that denaturation occurred by slow exchange between the duplexes and their component single strands. The single strands did not separate completely from the duplexes in aqueous solutions; instead, they appeared as small subpeaks on the tails of the duplex peaks. However, complete separation of the duplexes from their component single strands was observed at 20 °C in solutions containing 300 mM tetrapropylammonium ions. In addition, counterion condensation appears to be significantly reduced in DNA duplexes containing internal or bulge loops.

Entangled World of DNA Quadruplex Folds.

DNA quadruplexes participate in many biological functions. It takes up a variety of folds based on the sequence and environment. Here, a meticulous analysis of experimentally determined 437 quadruplex structures (433 PDBs) deposited in the PDB is carried out. The analysis reveals the modular representation of the quadruplex folds. Forty-eight unique quadruplex motifs (whose diversity arises out of the propeller, bulge, diagonal, and lateral loops that connect the quartets) are identified, leading to simple to complex inter/intramolecular quadruplex folds. The two-layered structural motifs are further classified into 33 continuous and 15 discontinuous motifs. While the continuous motifs can directly be extended to a quadruplex fold, the discontinuous motif requires an additional loop(s) to complete a fold, as illustrated here with examples. Similarly, higher-order quadruplex folds can also be represented by continuous or discontinuous motifs or their combinations. Such a modular representation of the quadruplex folds may assist in custom engineering of quadruplexes, designing motif-based drugs, and the prediction of the quadruplex structure. Furthermore, it could facilitate understanding of the role of quadruplexes in biological functions and diseases.

The Programmable Catalytic Core of 8-17 DNAzymes.

8-17 DNAzymes (8-17, 17E, Mg5, and 17EV1) are in vitro-selected catalytic DNA molecules that are capable of cleaving complementary RNAs. The conserved residues in their similar catalytic cores, together with the metal ions, were suggested to contribute to the catalytic reaction. Based on the contribution of the less conserved residues in the bulge loop residues (W12, A15, A15.0) and the internal stem, new catalytic cores of 8-17 DNAzymes were programmed. The internal stem CTC-GAG seems to be more favorable for the DNAzymes than CCG-GGC, while an extra W12.0 led to a significant loss of activity of DNAzymes, which is contrary to the positive effect of A15.0, by which a new active DNAzyme 17EM was derived. It conducts a faster reaction than 17E. It is most active in the presence of Pb2+, with the metal ion preference of Pb2+ >> Zn2+ > Mn2+ > Ca2+ ≈ Mg2+. In the Pb2+ and Zn2+-mediated reactions of 17EM and 17E, the same Na+- and pH dependence were also observed as what was observed for 17E and other 8-17 DNAzymes. Therefore, 17EM is another member of the 8-17 DNAzymes, and it could be applied as a potential biosensor for RNA and metal ions.

Solid Phase Synthesis and TAR RNA-Binding Activity of Nucleopeptides Containing Nucleobases Linked to the Side Chains via 1,4-Linked-1,2,3-triazole.

Nucleopeptides (NPs) represent synthetic polymers created by attaching nucleobases to the side chains of amino acid residues within peptides. These compounds amalgamate the characteristics of peptides and nucleic acids, showcasing a unique ability to recognize RNA structures. In this study, we present the design and synthesis of Fmoc-protected nucleobase amino acids (1,4-TzlNBAs) and a new class of NPs, where canonical nucleobases are affixed to the side chain of L-homoalanine (Hal) through a 1,4-linked-1,2,3-triazole (HalTzl). Fmoc-protected 1,4-TzlNBAs suitable for HalTzl synthesis were obtained via Cu(I)-catalyzed azide-alkyne cycloaddition (CuAAC) conjugation of Fmoc-L-azidohomoalanine (Fmoc-Aha) and N1- or N9-propargylated nucleobases or their derivatives. Following this, two trinucleopeptides, HalTzlAAA and HalTzlAGA, and the hexanucleopeptide HalTzlTCCCAG, designed to complement bulge and outer loop structures of TAR (trans-activation response element) RNA HIV-1, were synthesized using the classical solid-phase peptide synthesis (SPPS) protocol. The binding between HalTzls and fluorescently labeled 5'-(FAM(6))-TAR UCU and UUU mutant was characterized using circular dichroism (CD) and fluorescence spectroscopy. CD results confirmed the binding of HalTzls to TAR RNA, which was evident by a decrease in ellipticity band intensity around 265 nm during complexation. CD thermal denaturation studies indicated a relatively modest effect of complexation on the stability of TAR RNA structure. The binding of HalTzls at an equimolar ratio only marginally increased the melting temperature (Tm) of the TAR RNA structure, with an increment of less than 2 °C in most cases. Fluorescence spectroscopy revealed that HalTzlAAA and HalTzlAGA, complementary to UUU or UCU bulges, respectively, exhibited disparate affinities for the TAR RNA structure (with Kd ≈ 30 and 256 µM, respectively). Hexamer HalTzlTCCCAG, binding to the outer loop of TARUCU, demonstrated a moderate affinity with Kd ≈ 38 µM. This study demonstrates that newly designed HalTzls effectively bind the TAR RNA structure, presenting a potential new class of RNA binders and may be a promising scaffold for the development of a new class of antiviral drugs.

How sequence alterations enhance the stability and delay expansion of DNA triplet repeat domains

Abstract DNA sequence alterations within DNA repeat domains inexplicably enhance the stability and delay the expansion of interrupted repeat domains. Here we propose mechanisms that rationalise such unanticipated outcomes. Specifically, we describe how interruption of a DNA repeat domain restricts the ensemble space available to dynamic, slip out, repeat bulge loops by introducing energetic barriers to loop migration. We explain how such barriers arise because some possible loop isomers result in energetically costly mismatches in the duplex portion of the repeat domain. We propose that the reduced ensemble space is the causative feature for the observed delay in repeat DNA expansion. We further posit that the observed loss of the interrupting repeat in some expanded DNAs reflects the transient occupation of loop isomer positions that result in a mismatch in the duplex stem due to ‘leakiness’ in the energy barrier. We propose that if the lifetime of such a low probability event allows for recognition by the mismatch repair system, then ‘repair’ of the repeat interruption can occur; thereby rationalising the absence of the interruption in the final expanded DNA ‘product.’ Our proposed mechanistic pathways provide reasoned explanations for what have been described as ‘puzzling’ observations, while also yielding insights into a biomedically important set of coupled genotypic phenomena that map the linkage between DNA origami thermodynamics and phenotypic disease states.

Adjacent and Nonadjacent Dipyrimidine Photoproducts as Intrinsic Probes of DNA Secondary and Tertiary Structure†.

While the photochemistry of duplex DNA has been extensively studied, the photochemistry of nonduplex DNA structures is largely unexplored. Because the structure and stereochemistry of DNA photoproducts depend on the secondary structure and conformation of the DNA precursor, they can serve as intrinsic probes of DNA structure. This review focuses on the structures and stereoisomers of pyrimidine dimer photoproducts arising from adjacent and nonadjacent pyrimidines in A, B and denatured DNA, bulge loops, G-quadruplexes and reverse Hoogsteen hairpins and methods for their detection.

Evaluation of weak interactions of proteins and organic cations with DNA duplex structures

Molecular interactions and reactions in living cells occur with high background concentrations of organic compounds including proteins. Uncharged water-soluble polymers are commonly used cosolutes in studies on molecular crowding, and most studies argue about the effects of intracellular crowding based on results obtained using polymer cosolutes. Further investigations using protein crowders and organic cations are important in understanding the effects of cellular environments on nucleic acids with negatively charged surfaces. We assessed the effects of using model globular proteins, serum proteins, histone proteins, structurally flexible polypeptides, di- and polyamines, and uncharged polymers. Thermal stability analysis of DNA oligonucleotide structures revealed that unlike conventional polymer cosolutes, basic globular proteins (lysozyme and cytochrome c) at high concentrations stabilized long internal and bulge loop structures but not fully matched duplexes. The selective stabilization of long loop structures suggests preferential binding to unpaired nucleotides in loops through weak electrostatic interactions. Furthermore, the ability of the proteins to stabilize the loop structures was enhanced under macromolecular crowding conditions. Remarkably, the effects of basic proteins on the stability of fully matched duplexes were dissimilar to those of basic amino-acid-rich polypeptides and polyamines. This study provides new insights into the interaction of nucleic acid structures with organic cations.

Bulge-Forming miRNases Cleave Oncogenic miRNAs at the Central Loop Region in a Sequence-Specific Manner.

The selective degradation of disease-associated microRNA is promising for the development of new therapeutic approaches. In this study, we engineered a series of bulge-loop-forming oligonucleotides conjugated with catalytic peptide [(LeuArg)2Gly]2 (BC–miRNases) capable of recognizing and destroying oncogenic miR-17 and miR-21. The principle behind the design of BC–miRNase is the cleavage of miRNA at a three-nucleotide bulge loop that forms in the central loop region, which is essential for the biological competence of miRNA. A thorough study of mono- and bis-BC–miRNases (containing one or two catalytic peptides, respectively) revealed that: (i) the sequence of miRNA bulge loops and neighbouring motifs are of fundamental importance for efficient miRNA cleavage (i.e., motifs containing repeating pyrimidine–A bonds are more susceptible to cleavage); (ii) the incorporation of the second catalytic peptide in the same molecular scaffold increases the potency of BC–miRNase, providing a complete degradation of miR-17 within 72 h; (iii) the synergetic co-operation of BC–miRNases with RNase H accelerates the rate of miRNA catalytic cleavage by both the conjugate and the enzyme. Such synergy allows the rapid destruction of constantly emerging miRNA to maintain sufficient knockdown and achieve a desired therapeutic effect.

A Test and Refinement of Folding Free Energy Nearest Neighbor Parameters for RNA Including N6-Methyladenosine

RNA folding free energy change parameters are widely used to predict RNA secondary structure and to design RNA sequences. These parameters include terms for the folding free energies of helices and loops. Although the full set of parameters has only been traditionally available for the four common bases and backbone, it is well known that covalent modifications of nucleotides are widespread in natural RNAs. Covalent modifications are also widely used in engineered sequences. We recently derived a full set of nearest neighbor terms for RNA that includes N6-methyladenosine (m6A). In this work, we test the model using 98 optical melting experiments, matching duplexes with or without N6-methylation of A. Most experiments place RRACH, the consensus site of N6-methylation, in a variety of contexts, including helices, bulge loops, internal loops, dangling ends, and terminal mismatches. For matched sets of experiments that include either A or m6A in the same context, we find that the parameters for m6A are as accurate as those for A. Across all experiments, the root mean squared deviation between estimated and experimental free energy changes is 0.67 kcal/mol. We used the new experimental data to refine the set of nearest neighbor parameter terms for m6A. These parameters enable prediction of RNA secondary structures including m6A, which can be used to model how N6-methylation of A affects RNA structure.

In Vitro Evaluation of Bis-3-Chloropiperidines as RNA Modulators Targeting TAR and TAR-Protein Interaction.

After a long limbo, RNA has gained its credibility as a druggable target, fully earning its deserved role in the next generation of pharmaceutical R&D. We have recently probed the trans-activation response (TAR) element, an RNA stem–bulge–loop domain of the HIV-1 genome with bis-3-chloropiperidines (B-CePs), and revealed the compounds unique behavior in stabilizing TAR structure, thus impairing in vitro the chaperone activity of the HIV-1 nucleocapsid (NC) protein. Seeking to elucidate the determinants of B-CePs inhibition, we have further characterized here their effects on the target TAR and its NC recognition, while developing quantitative analytical approaches for the study of multicomponent RNA-based interactions.

"Bind, cleave and leave": multiple turnover catalysis of RNA cleavage by bulge-loop inducing supramolecular conjugates.

Antisense sequence-specific knockdown of pathogenic RNA offers opportunities to find new solutions for therapeutic treatments. However, to gain a desired therapeutic effect, the multiple turnover catalysis is critical to inactivate many copies of emerging RNA sequences, which is difficult to achieve without sacrificing the sequence-specificity of cleavage. Here, engineering two or three catalytic peptides into the bulge–loop inducing molecular framework of antisense oligonucleotides achieved catalytic turnover of targeted RNA. Different supramolecular configurations revealed that cleavage of the RNA backbone upon sequence-specific hybridization with the catalyst accelerated with increase in the number of catalytic guanidinium groups, with almost complete demolition of target RNA in 24 h. Multiple sequence-specific cuts at different locations within and around the bulge–loop facilitated release of the catalyst for subsequent attacks of at least 10 further RNA substrate copies, such that delivery of only a few catalytic molecules could be sufficient to maintain knockdown of typical RNA copy numbers. We have developed fluorescent assay and kinetic simulation tools to characterise how the limited availability of different targets and catalysts had restrained catalytic reaction progress considerably, and to inform how to accelerate the catalytic destruction of shorter linear and larger RNAs even further.

A least square support vector machine approach based on bvRNA-GA for modeling photovoltaic systems

Accurate model plays an important role in designing, assessing, and controlling photovoltaic (PV) systems. In this work, the least-squares support vector machine (LSSVM) is adopted to model the current–voltage (V–I) characteristic curves of different PV systems. A novel RNA genetic algorithm (bvRNA-GA) is proposed to determine the parameters of LSSVM. The bvRNA-GA is featured by designing the bulge loop crossover operator and the virus-induced mutation operator, they are employed to balance the exploration and exploitation capacities. Different experiments with 10 benchmark functions are conducted to show that the search efficiency of bvRNA-GA is better than the other four state-of-art algorithms. The outputs of bvRNA-GA optimized LSSVM models can better agree with the real outputs of different PV systems, the modeling results demonstrate the effectiveness of bvRNA-GA in solving real-world problems.

Differential repair enzyme-substrate selection within dynamic DNA energy landscapes.

We demonstrate that reshaping of the dynamic, bulged-loop energy landscape of DNA triplet repeat ensembles by the presence of an abasic site alters repair outcomes by the APE1 enzyme. This phenomenon depends on the structural context of the lesion, despite the abasic site always having the same neighbors in sequence space. We employ this lesion-induced redistribution of DNA states and a kinetic trap to monitor different occupancies of the DNA bulge loop states. We show how such dynamic redistribution and associated differential occupancies of DNA states impact APE1 repair outcomes and APE1 induced interconversions. We correlate the differential biophysical properties of the dynamic, DNA ensemble states, with their ability to be recognized and processed as substrates by the APE1 DNA repair enzyme. Enzymatic digestions and biophysical characterizations reveal that APE1 cuts a fraction (10-12%) of the dynamic, rollameric substrates within the initial kinetic distribution. APE1 interactions also 'induce' rollamer redistribution from a kinetically trapped distribution to an equilibrium distribution, the latter not containing viable APE1 substrates. We distinguish between kinetically controlled ensemble (re)distributions of potential DNA substrates, versus thermodynamically controlled ensemble (re)distribution; features of importance to DNA regulation. We conclude that APE1 activity catalyzes/induces ensembles that represent the thermodynamically optimal loop distribution, yet which also are nonviable substrate states for abasic site cleavage by APE1. We propose that by inducing substrate redistributions in a dynamic energy landscape, the enzyme actually reduces the available substrate competent species for it to process, reflective of a regulatory mechanism for enzymatic self-repression. If this is a general phenomenon, such a consequence would have a profound impact on slowing down and/or misdirecting DNA repair within dynamic energy landscapes, as exemplified here within triplet repeat domains. In short, APE1-instigated redistribution of potential substrates induces a preferred pathway to an equilibrium ensemble of enzymatically incompetent states.

Vfold2D-MC: A Physics-Based Hybrid Model for Predicting RNA Secondary Structure Folding.

Accurate prediction of RNA structure and folding stability has a far-reaching impact on our understanding of RNA functions. Here we develop Vfold2D-MC, a new physics-based model, to predict RNA structure and folding thermodynamics from the sequence. The model employs virtual bond-based coarse-graining of RNA backbone conformation and generates RNA conformations through Monte Carlo sampling of the bond angles and torsional angles of the virtual bonds. Using a coarse-grained statistical potential derived from the known structures, we assign each conformation with a statistical weight. The weighted average over the conformational ensemble gives the entropy and free energy parameters for the hairpin, bulge, and internal loops, and multiway junctions. From the thermodynamic parameters, we predict RNA structures, melting curves, and structural changes from the sequence. Theory-experiment comparisons indicate that Vfold2D-MC not only gives improved structure predictions but also enables the interpretation of thermodynamic results for different RNA structures, including multibranched junctions. This new model sets a promising framework to treat more complicated RNA structures, such as pseudoknotted and intramolecular kissing loops, for which experimental thermodynamic parameters are often unavailable.

Thermodynamic Examination of P3a Bulge Loop of a Glycine Riboswitch

The glycine riboswitch (GRS) is a non-coding RNA which regulates expression of glycine cleavage and transport genes through transcription attenuation. The riboswitch contains a tandem aptamer which binds glycine to cause structural changes in the downstream expression platform. The 4×3-nt P3a bulge loop motif is implicated in glycine binding in both aptamers. Glycine hydrogen bonds to U78 in the bulge loop along with at least one magnesium ion in the crystal structure. We designed RNA constructs containing a 4×3-nt bulge loop to measure thermodynamics of the P3a bulge loop motif seen in aptamer 1 of V. cholerae. Our goalwasto examine if the bulge loop would bind to glycine in the absence of the larger RNA. We measured the binding of glycine using isothermal calorimetry (ITC) in 10 mM MgCl2 and 1 M KCl buffers at pH 7. A U78A modified construct was also analyzed. Although both structures formed in solution, the wild-type construct did not bind glycine, suggesting that the bulge loop motif does not independently support glycine binding within the riboswitch. Initial data analysis on hybridization of the wild-type RNA structure showed an increase in stability of ~3 kcal/mol in the presence of 10 mM MgCl2 over 1 M KCl. The thermodynamics of the U78A construct was not affected by the presence of magnesium, with only a 0.03 kcal/mol increase in stability in 10 mM MgCl2 buffer.

Bis-3-Chloropiperidines Targeting TAR RNA as A Novel Strategy to Impair the HIV-1 Nucleocapsid Protein.

Specific RNA sequences regulate functions essential to life. The Trans-Activation Response element (TAR) is an RNA stem–bulge–loop structure involved in several steps of HIV-1 replication. In this work, we show how RNA targeting can inhibit HIV-1 nucleocapsid (NC), a highly conserved protein known to catalyze nucleic acid melting and strand transfers during reverse transcription. Our RNA targeting strategy consists of the employment of bis-3-chloropiperidines (B-CePs) to impair RNA melting through bifunctional alkylation. Specific interactions between B-CePs and TAR RNA were analytically investigated by gel electrophoresis and mass spectrometry, allowing the elucidation of B-CePs’ recognition of TAR, and highlighting an RNA-directed mechanism of protein inhibition. We propose that B-CePs can freeze TAR tridimensional conformation, impairing NC-induced dynamics and finally inhibiting its functions in vitro.

Multipolar electrostatics for hairpin and pseudoknots in RNA: Improving the accuracy of force field potential energy function.

Molecular dynamics (MD) simulations that rely on force field methods has been widely used to explore the structure and function of RNAs. However, the current commonly used force fields are limited by the electrostatic description offered by atomic charge, dipole and at most quadrupole moments, failing to capture the anisotropic picture of electronic features. Actually, the distribution of electrons around atomic nuclei is not spherically symmetric but is geometry dependent. A multipolar electrostatic model based on high rank multipole moments is described in this work, which allows us to combine polarizability and anisotropy of electron density. RNA secondary structure was taken as a research system, and its substructures including stem, loops (hairpin loop, bulge loop, internal loop, and multi-branch loop), and pseudoknots (H-type and K-type) were investigated, respectively, as well as the hairpin. First, the atom-atom electrostatic properties derived from one chain of a duplex RNA 2MVY in our previous work (Ref. 58) were measured by the pilot RNA systems of hairpin, hairpin loop, stem, and H-type pseudoknot, respectively. The prediction results were not satisfactory. Consequently, to obtain a general set of electrostatic parameters for RNA force fields, the convergence behavior of the atom-atom electrostatic interactions in the pilot RNA systems was explored using high rank atomic multipole moments. The pilot RNA systems were cut into four types of different-sized molecular fragments, and the single nucleotide fragment and nucleotide-paired fragment proved to be the most reasonable systems for base-unpairing regions and base-pairing regions to investigate the convergence behavior of all types of atom-atom electrostatic interactions, respectively. Transferability of the electrostatic properties drawn from the pilot RNA systems to the corresponding test systems was also investigated. Furthermore, the convergence behavior of atomic electrostatic interactions in other substructures including bulge loop, internal loop, multi-branch loop, and K-type pseudoknot was expected to be modeled via the hairpin.

Cationic FeII Triplex-Forming Metallohelices as DNA Bulge Binders.

Bulges are essential structural elements in nucleic acids. The detection and targeting of bulged DNA sequences are highly important. Small molecules capable of targeting DNA bulges have attracted considerable attention because they cannot only be used as reagents for bulge recognition, but also as potential therapeutic drugs. Herein, the interactions of DNA duplexes, containing bulges of various sizes and base compositions, with a series of FeII triplex-forming metallohelices are reported. The results obtained, with the aid of molecular biophysics methods, show that the investigated metallohelices prefer to bind to bulged DNA, rather than double-stranded DNA, and that their binding affinities towards bulges differ among individual metallohelices. Moreover, their binding affinities towards bulges strongly depend on the bulge size and the base composition of the bulge loop. The investigated metallohelices can enter eukaryotic cells and accumulate in the cell nucleus, allowing them to interact with nucleic acids. Hence, it is reasonable to suggest that the interaction of metallohelices with nucleic acid bulges might contribute to the mechanism of their biological activity.

Shape-Selective Targeting on RNA Bulges by Peptidomimetic Metallohelices.

RNA bulges represent one of the most common motifs in the RNA secondary structure and serve in a variety of biological functions. Compounds stabilizing RNA bulges are important for probing RNA structure and function and for therapy of some diseases. Here, the ability of a series of enantiomeric pairs of optically pure bimetallic metallohelices with different flexible linkers to target various RNA bulges is investigated. The results show that binding affinities of the metallohelices to bulged RNA differ and strongly depend on the size of the bulge and the base composition of the bulge loop. Notably, the shorter, more compact, and less flexible metallohelices bind to RNA bulges most efficiently and selectively. Interestingly, the ability of the metallohelices to bind to RNA bulges correlates with their previously reported antimicrobial activity, which suggests that the selective recognition of bulged regions in RNA by the metallohelices might also contribute to their biological activity.

Heat Capacity Changes (ΔCp) for Interconversions between Differentially-Ordered DNA States within Physiological Temperature Domains: Implications for Biological Regulatory Switches.

Knowledge of differences in heat capacity changes (ΔCp) between biopolymer states provides essential information about the temperature dependence of the thermodynamic properties of these states, while also revealing insights into the nature of the forces that drive the formation of functional and dysfunctional biopolymer "order." In contrast to proteins, for nucleic acids there is a dearth of direct experimental determination of this information-rich parameter, a deficiency that compromises interpretations of the ever-increasing thermodynamic analyses of nucleic acid properties; particularly as they relate to differential nucleic acid (meta)stability states and their potential biological functions. Here we demonstrate that such heat capacity differences, in fact, exist not only between traditionally measured native to fully unfolded (assumed "random coil") DNA states, but also between competing order-to-order transformations. We illustrate the experimental approach by measuring the heat capacity change between "native"/ordered, sequence homologous, "isomeric" DNA states that differ in conformation but not sequence. Importantly, these heat capacity differences occur within biologically relevant temperature ranges. In short, we describe a new and general method to measure the value of such heat capacity differences anywhere in experimentally accessible conformational and temperature space; in this case, between two metastable bulge loop states, implicated in DNA expansion diseases, and their competing, fully paired, thermodynamically more stable duplex states. This measurement reveals a ΔCp of 61 ± 7 cal molbp -1 K -1. Such heat capacity differences between competing DNA "native" ensemble states must be considered when evaluating equilibria between different DNA "ordered" conformations, including the assessment of the differential stabilizing forces and potential biological functions of competing DNA "structured" motifs.