Molecular dynamics reveals the stabilizing role of loop closing residues in kissing interactions: comparison between TAR-TAR* and TAR-aptamer

Strong and selective interactions of palmatine with G-rich sequences in TRF2 promoter; experimental and computational studies

G-rich sequences have the potential to fold into G-quadruplexes (GQs). G-quadruplexes, particularly those positioned in the regulatory regions of proto-oncogenes, have recently garnered attention in anti-cancer drug design. A thermal FRET assay was employed to conduct preliminary screening of various alkaloids, aiming to identify stronger interactions with a specific set of G-rich double-labeled oligonucleotides in both K + and Na + buffers. These oligonucleotides were derived from regions associated with Kit, Myc, Ceb, Bcl2, human telomeres, and potential G-quadruplex forming sequences found in the Nrf2 and Trf2 promoters. Palmatine generally increased the stability of different G-rich sequences into their folded GQ structures, more or less in a concentration dependent manner. The thermal stability and interaction of palmatine was further studied using transition FRET (t-FRET), CD and UV-visible spectroscopy and molecular dynamics simulation methods. Palmatine showed the strongest interaction with T RF2 in both K+ and Na+ buffers even at equimolar concentration ratio. T-FRET studies revealed that palmatine has the potential to disrupt double-strand formation by the T RF2 sequence in the presence of its complementary strand. Palmatine exhibits a stronger interaction with G-rich strand DNA, promoting its folding into G-quadruplex structures. It is noteworthy that palmatine exhibits the strongest interaction with T RF2, which is the shortest sequence among the G-rich oligonucleotides studied, featuring only one nucleotide for two of its loops. Palmatine represents a suitable structure for drug design to develop more specific ligands targeting G-quadruplexes. Whether palmatine can also affect the expression of the T RF2 gene requires further studies.Communicated by Ramaswamy H. Sarma.

Engineering Light-Up Aptamers for the Detection of RNA Hairpins through Kissing Interaction

Aptasensors are biosensors that include aptamers for detecting a target of interest. We engineered signaling aptasensors for the detection of RNA hairpins from the previously described malachite green (MG) RNA aptamer. The top part of this imperfect hairpin aptamer was modified in such a way that it can engage loop-loop (so-called kissing) interactions with RNA hairpins displaying partly complementary apical loops. These newly derived oligonucleotides named malaswitches bind their cognate fluorogenic ligand (MG) exclusively when RNA-RNA kissing complexes are formed, whereas MG does not bind to malaswitches alone. Consequently, the formation of the ternary target RNA-malaswitch RNA-MG complex results in fluorescence emission, and malaswitches constitute sensors for detecting RNA hairpins. Malaswitches were designed that specifically detect precursors of microRNAs let7b and miR-206.

Triggering nucleic acid nanostructure assembly by conditional kissing interactions

Nucleic acids are biomolecules of amazing versatility. Beyond their function for information storage they can be used for building nano-objects. We took advantage of loop–loop or kissing interactions between hairpin building blocks displaying complementary loops for driving the assembly of nucleic acid nano-architectures. It is of interest to make the interaction between elementary units dependent on an external trigger, thus allowing the control of the scaffold formation. To this end we exploited the binding properties of structure-switching aptamers (aptaswitch). Aptaswitches are stem–loop structured oligonucleotides that engage a kissing complex with an RNA hairpin in response to ligand-induced aptaswitch folding. We demonstrated the potential of this approach by conditionally assembling oligonucleotide nanorods in response to the addition of adenosine.

Electrostatics Explains the Position-Dependent Effect of G⋅U Wobble Base Pairs on the Affinity of RNA Kissing Complexes

In the RNA realm, non-Watson-Crick base pairs are abundant and can affect both the RNA 3D structure and its function. Here, we investigated the formation of RNA kissing complexes in which the loop-loop interaction is modulated by non-Watson-Crick pairs. Mass spectrometry, surface plasmon resonance, and UV-melting experiments show that the G⋅U wobble base pair favors kissing complex formation only when placed at specific positions. We tried to rationalize this effect by molecular modeling, including molecular mechanics Poisson-Boltzmann surface area (MMPBSA) thermodynamics calculations and PBSA calculations of the electrostatic potential surfaces. Modeling reveals that the G⋅U stabilization is due to a specific electrostatic environment defined by the base pairs of the entire loop-loop region. The loop is not symmetric, and therefore the identity and position of each base pair matters. Predicting and visualizing the electrostatic environment created by a given sequence can help to design specific kissing complexes with high affinity, for potential therapeutic, nanotechnology or analytical applications.

Single-molecule observations of RNA-RNA kissing interactions in a DNA nanostructure

RNA molecules uniquely form a complex through specific hairpin loops, called a kissing complex. The kissing complex is widely investigated and used for the construction of RNA nanostructures. Molecular switches have also been created by combining a kissing loop and a ligand-binding aptamer to control the interactions of RNA molecules. In this study, we incorporated two kinds of RNA molecules into a DNA origami structure and used atomic force microscopy to observe their ligand-responsive interactions at the single-molecule level. We used a designed RNA aptamer called GTPswitch, which has a guanosine triphosphate (GTP) responsive domain and can bind to the target RNA hairpin named Aptakiss in the presence of GTP. We observed shape changes of the DNA/RNA strands in the DNA origami, which are induced by the GTPswitch, into two different shapes in the absence and presence of GTP, respectively. We also found that the switching function in the nanospace could be improved by using a cover strand over the kissing loop of the GTPswitch or by deleting one base from this kissing loop. These newly designed ligand-responsive aptamers can be used for the controlled assembly of the various DNA and RNA nanostructures.

On the conformational stability of the smallest RNA kissing complexes maintained through two G·C base pairs

Two identical 5'GACG3' tetra-loop motifs with different stem sequences (called H2 and H3) are found in the 5' end region of Moloney Murine Leukemia Virus (MMLV) genomic RNA. They play important roles in RNA dimerization and encapsidation through two identical tetra-loops (5'GACG3') forming a loop-to-loop kissing complex, the smallest RNA kissing complex ever found in nature. We examined the effects of a loop-closing base pair as well as a stem sequence on the conformational stability of the kissing complex. UV melting analysis and gel electrophoresis were performed on eight RNA sequences mimicking the H2 and H3 hairpin tetra-loops with variation in loop-closing base pairs. Our results show that changing the loop-closing base pair from the wildtype (5'A·U3' for H3, 5'U·A3' for H2) to 5'G·C3'/5'C·G3' has significant effect on the stability of the kissing complexes: the substitution to 5'C·G3' significantly decreases both thermal and mechanical stability, while switching to the 5'G·C3' significantly increases the mechanical stability only. The kissing complexes with the wildtype loop-closing base pairs (5'A·U3' for H3 and 5'U·A3' for H2) show different stability when attached to a different stem sequence (H2 stem vs. H3 stem). This suggests that not only the loop-closing base pair itself, but also the stem sequence, affects the conformational stability of the RNA kissing complex.

Riboswitches based on kissing complexes for the detection of small ligands

Biosensors derived from aptamers were designed for which folding into a hairpin shape is triggered by binding of the cognate ligand. These aptamers (termed aptaswitches) thus switch between folded and unfolded states in the presence and absence of the ligand, respectively. The apical loop of the folded aptaswitch is recognized by a second hairpin called the aptakiss through loop-loop or kissing interactions, whereas the aptakiss does not bind the unfolded aptaswitch. Therefore, the formation of a kissing complex signals the presence of the ligand. Aptaswitches were designed that enable the detection of GTP and adenosine in a specific and quantitative manner by surface plasmon resonance when using a grafted aptakiss or in solution by anisotropy measurement with a fluorescently labeled aptakiss. This approach is generic and can potentially be extended to the detection of any molecule for which hairpin aptamers have been identified, as long as the apical loop is not involved in ligand binding.

HIV-1 DIS stem loop forms an obligatory bent kissing intermediate in the dimerization pathway

The HIV-1 dimerization initiation sequence (DIS) is a conserved palindrome in the apical loop of a conserved hairpin motif in the 5′-untranslated region of its RNA genome. DIS hairpin plays an important role in genome dimerization by forming a ‘kissing complex’ between two complementary hairpins. Understanding the kinetics of this interaction is key to exploiting DIS as a possible human immunodeficiency virus (HIV) drug target. Here, we present a single-molecule Förster resonance energy transfer (smFRET) study of the dimerization reaction kinetics. Our data show the real-time formation and dissociation dynamics of individual kissing complexes, as well as the formation of the mature extended duplex complex that is ultimately required for virion packaging. Interestingly, the single-molecule trajectories reveal the presence of a previously unobserved bent intermediate required for extended duplex formation. The universally conserved A272 is essential for the formation of this intermediate, which is stabilized by Mg²⁺, but not by K⁺ cations. We propose a 3D model of a possible bent intermediate and a minimal dimerization pathway consisting of three steps with two obligatory intermediates (kissing complex and bent intermediate) and driven by Mg²⁺ ions.

Unmasking the information encoded as structural motifs of viral RNA genomes: a potential antiviral target

RNA viruses show enormous capacity to evolve and adapt to new cellular and molecular contexts, a consequence of mutations arising from errors made by viral RNA-dependent RNA polymerase during replication. Sequence variation must occur, however, without compromising functions essential for the completion of the viral cycle. RNA viruses are safeguarded in this respect by their genome carrying conserved information that does not code only for proteins but also for the formation of structurally conserved RNA domains that directly perform these critical functions. Functional RNA domains can interact with other regions of the viral genome and/or proteins to direct viral translation, replication and encapsidation. They are therefore potential targets for novel therapeutic strategies. This review summarises our knowledge of the functional RNA domains of human RNA viruses and examines the achievements made in the design of antiviral compounds that interfere with their folding and therefore their function.

Inhibition of HIV Replication by Cyclic and Hairpin PNAs Targeting the HIV-1 TAR RNA Loop

Human immunodeficiency virus-1 (HIV-1) replication and gene expression entails specific interaction of the viral protein Tat with its transactivation responsive element (TAR), to form a highly stable stem-bulge-loop structure. Previously, we described triphenylphosphonium (TPP) cation-based vectors that efficiently deliver nucleotide analogs (PNAs) into the cytoplasm of cells. In particular, we showed that the TPP conjugate of a linear 16-mer PNA targeting the apical stem-loop region of TAR impedes Tat-mediated transactivation of the HIV-1 LTR in vitro and also in cell culture systems. In this communication, we conjugated TPP to cyclic and hairpin PNAs targeting the loop region of HIV-1 TAR and evaluated their antiviral efficacy in a cell culture system. We found that TPP-cyclic PNAs containing only 8 residues, showed higher antiviral potency compared to hairpin PNAs of 12 or 16 residues. We further noted that the TPP-conjugates of the 8-mer cyclic PNA as well as the 16-mer linear PNA displayed similar antiviral efficacy. However, cyclic PNAs were shown to be highly specific to their target sequences. This communication emphasizes on the importance of small constrained cyclic PNAs over both linear and hairpin structures for targeting biologically relevant RNA hairpins.

In vitro selection of RNA aptamers derived from a genomic human library against the TAR RNA element of HIV-1

The transactivating responsive (TAR) element is a RNA hairpin located in the 5' untranslated region of HIV-1 mRNA. It is essential for full-length transcription of the retroviral genome and therefore for HIV-1 replication. Hairpin aptamers that generate highly stable and specific complexes with TAR were previously identified, thus decreasing the level of TAR-dependent expression in cultured cells [Kolb, G., et al. (2006) RNA Biol. 3, 150-156]. We performed genomic SELEX against TAR using a human RNA library to identify human transcripts that might interact with the retroviral genome through loop-loop interactions and potentially contribute to the regulation of TAR-mediated processes. We identified a genomic aptamer termed a1 that folds as a hairpin with an apical loop complementary to five nucleotides of the TAR hexanucleotide loop. Surface plasmon resonance experiments performed on a truncated or mutated version of the a1 aptamer, in the presence of the Rop protein of Escherichia coli, indicate the formation of a highly stable a1-TAR kissing complex. The 5' ACCCAG loop of a1 constitutes a new motif of interaction with the TAR loop.

Exploring TAR-RNA aptamer loop-loop interaction by X-ray crystallography, UV spectroscopy and surface plasmon resonance

In HIV-1, trans-activation of transcription of the viral genome is regulated by an imperfect hairpin, the trans-activating responsive (TAR) RNA element, located at the 5′ untranslated end of all viral transcripts. TAR acts as a binding site for viral and cellular proteins. In an attempt to identify RNA ligands that would interfere with the virus life-cycle by interacting with TAR, an in vitro selection was previously carried out. RNA hairpins that formed kissing-loop dimers with TAR were selected [Ducongé F. and Toulmé JJ (1999) RNA, 5:1605–1614]. We describe here the crystal structure of TAR bound to a high-affinity RNA aptamer. The two hairpins form a kissing complex and interact through six Watson–Crick base pairs. The complex adopts an overall conformation with an inter-helix angle of 28.1°, thus contrasting with previously reported solution and modelling studies. Structural analysis reveals that inter-backbone hydrogen bonds between ribose 2′ hydroxyl and phosphate oxygens at the stem-loop junctions can be formed. Thermal denaturation and surface plasmon resonance experiments with chemically modified 2′-O-methyl incorporated into both hairpins at key positions, clearly demonstrate the involvement of this intermolecular network of hydrogen bonds in complex stability.

Real time analysis of the RNAI-RNAII-Rop complex by surface plasmon resonance: from a decaying surface to a standard kinetic analysis

RNA loop-loop complexes are motifs that regulate biological functions in both prokaryotic and eukaryotic organisms. In E. coli, RNAI, an antisense RNA encoded by the ColE1 plasmid, regulates the plasmid replication by recognizing through loop-loop interactions RNAII, the RNA primer that binds to the plasmidic DNA to initiate the replication. Rop, a plasmid-encoded homodimeric protein interacts with this transient RNAI-RNAII kissing complex. A surface plasmon resonance (SPR)-based biosensor was used to investigate this protein-nucleic acid ternary complex, at 5 degrees C, in experimental conditions such as the protein binds either to the loop-loop complex while it dissociates or to a saturated stable RNAI-RNAII surface. The results show that RNAI hairpin dissociates from the RNAII surface 110 times slower in the presence of Rop than in its absence. Rop binds to RNAI-RNAII with an on-rate of 3.6 x 10(6) M(-1) s(-1) and an off-rate of 0.11 s(-1), resulting in a binding equilibrium constant equal to 31 nM. A Scatchard-plot analysis of the interaction monitored by SPR confirms a 1:1 complex of Rop and RNAI-RNAII as observed for non-natural Rop-loop-loop complexes.

Liquid-crystal NMR structure of HIV TAR RNA bound to its SELEX RNA aptamer reveals the origins of the high stability of the complex

Transactivation-response element (TAR) is a stable stem-loop structure of HIV RNA, which plays a crucial role during the life cycle of the virus. The apical loop of TAR acts as a binding site for essential cellular cofactors required for the replication of HIV. High-affinity aptamers directed against the apical loop of TAR have been identified by the SELEX approach. The RNA aptamers with the highest affinity for TAR fold as hairpins and form kissing complexes with the targeted RNA through loop-loop interactions. The aptamers with the strongest binding properties all possess a GA base pair combination at the loop-closing position. Using liquid-crystal NMR methodology, we have obtained a structural model in solution of a TAR-aptamer kissing complex with an unprecedented accuracy. This high-resolution structure reveals that the GA base pair is unilaterally shifted toward the 5' strand and is stabilized by a network of intersugar hydrogen bonds. This specific conformation of the GA base pair allows for the formation of two supplementary stable base-pair interactions. By systematic permutations of the loop-closing base pair, we establish that the identified atomic interactions, which form the basis for the high stability of the complex, are maintained in several other kissing complexes. This study rationalizes the stabilizing role of the loop-closing GA base pairs in kissing complexes and may help the development or improvement of drugs against RNA loops of viruses or pathogens as well as the conception of biochemical tools targeting RNA hairpins involved in important biological functions.

Self-association of adenine-dependent hairpin ribozymes

Hairpin ribozymes are flexible molecules that catalyse reversible self-cleavage after the docking of two independently folded internal loops, A and B. The activities, self-association and structures in solution of two 85 base adenine-dependent hairpin ribozymes (ADHR1 and ADHR2) were studied by native gel electrophoresis, analytical centrifugation, and small angle neutron scattering. Bi-molecular RNA interactions such as linear-linear, loop-loop, loop-linear or kissing interactions have been found to be important in the control of various biological functions, and hairpin loops present rich potential for establishing both intra- and intermolecular interactions through standard Watson-Crick base pairing or non-canonical interactions. Similar results were obtained for ADHR1 and ADHR2. At room temperature, they indicated end-to-end self-association of the ribozymes in rod-like structures with a cross-section corresponding to two double strands side-by-side. Dimers, which predominate at low concentration ( approximately 0.1 mg/ml), associate into longer rods, with increasing concentration ( approximately 1 mg/ml). Above 65 degrees C, the dimers and rods dissociated into compact monomers, with a radius of gyration similar to that of tRNA (about 70 bases). The dimers were non-active for catalysis, which suggests that dimer formation, probably by preventing the correct docking of loops A and B, could act as an inhibition mechanism for the regulation of hairpin ribozyme catalysis.

Antiviral aptamers

Aptamers are rare nucleic acid ligands, which can be concocted in the laboratory from the randomized pool of molecules by affinity and amplification processes. Aptamers have several properties as they can be applied complementarily to antibodies and have several advantages over antibodies. In the past, several aptamers have been selected with a view to develop antiviral agents for therapeutic applications. This review summarizes potent antiviral aptamers and their strategies to prevent the viral replication.

NMR structure of a kissing complex formed between the TAR RNA element of HIV-1 and a LNA-modified aptamer

The trans-activating responsive (TAR) RNA element located in the 5′ untranslated region of the HIV-1 genome is a 57-nt imperfect stem-loop essential for the viral replication. TAR regulates transcription by interacting with both viral and cellular proteins. RNA hairpin aptamers specific for TAR were previously identified by in vitro selection [Ducongé,F. and Toulmé,J.J. (1999) In vitro selection identifies key determinants for loop-loop interactions: RNA aptamers selective for the TAR RNA element of HIV-1. RNA, 5, 1605–1614]. These aptamers display a 5′-GUCCCAGA-3′ consensus apical loop, partially complementary to the TAR one, leading to the formation of a TAR–aptamer kissing complex. The conserved GA combination (underlined in the consensus sequence) has been shown to be crucial for the formation of a highly stable complex. To improve the nuclease resistance of the aptamer and to increase its affinity for TAR, locked nucleic acid (LNA) nucleotides were introduced in the aptamer apical loop. LNA are nucleic acids analogues that contain a 2′-O,4′-C methylene linkage and that raise the thermostablity of duplexes. We solved the NMR solution structure of the TAR–LNA-modified aptamer kissing complex. Structural analysis revealed the formation of a non-canonical G•A pair leading to increased stacking at the stem-loop junction. Our data also showed that the introduction of LNA residues provides an enhanced stability while maintaining a normal Watson–Crick base pairing with a loop–loop conformation close to an A-type.

Determination of refractive index increment ratios for protein-nucleic acid complexes by surface plasmon resonance

Three nucleic acid-protein complexes of 1:1 stoichiometry were analyzed by surface plasmon resonance on a Biacore biosensor to test whether or not proteins and nucleic acids yielded similar refractive index increments on binding. The expected maximum response in resonance units, (RU(exp))(max), and the observed one, (RU(obs))(max), on saturation of immobilized targets by interacting partners were compared to determine the ratio of (deltan/deltaC)(protein) to (deltan/deltaC)(nucleic acid), where n is the refractive index at the surface and C is the concentration of one partner. Our results suggest that proteins and nucleic acids behave similarly and that the discrepancy between the expected and observed maximum responses for such complexes reflects inaccurate evaluation of the binding responses. Therefore, no correction of the instrument response is required for protein and nucleic acid interaction studies on a Biacore biosensor.

Coupling of fast and slow modes in the reaction pathway of the minimal hammerhead ribozyme cleavage

By employing classical molecular dynamics, correlation analysis of coupling between slow and fast dynamical modes, and free energy (umbrella) sampling using classical as well as mixed quantum mechanics molecular mechanics force fields, we uncover a possible pathway for phosphoryl transfer in the self-cleaving reaction of the minimal hammerhead ribozyme. The significance of this pathway is that it initiates from the minimal hammerhead crystal structure and describes the reaction landscape as a conformational rearrangement followed by a covalent transformation. The delineated mechanism is catalyzed by two metal (Mg(2+)) ions, proceeds via an in-line-attack by CYT 17 O2' on the scissile phosphorous (ADE 1.1 P), and is therefore consistent with the experimentally observed inversion configuration. According to the delineated mechanism, the coupling between slow modes involving the hammerhead backbone with fast modes in the cleavage site appears to be crucial for setting up the in-line nucleophilic attack.

Conformational changes in the expression domain of the Escherichia coli thiM riboswitch

The thiM riboswitch contains an aptamer domain that adaptively binds the coenzyme thiamine pyrophosphate (TPP). The binding of TPP to the aptamer domain induces structural rearrangements that are relayed to a second domain, the so-called expression domain, thereby interfering with gene expression. The recently solved crystal structures of the aptamer domains of the thiM riboswitches in complex with TPP revealed how TPP stabilizes secondary and tertiary structures in the RNA ligand complex. To understand the global modes of reorganization between the two domains upon metabolite binding the structure of the entire riboswitch in presence and absence of TPP needs to be determined. Here we report the secondary structure of the entire thiM riboswitch from Escherichia coli in its TPP-free form and its transition into the TPP-bound variant, thereby depicting domains of the riboswitch that serve as communication links between the aptamer and the expression domain. Furthermore, structural probing provides an explanation for the lack of genetic control exerted by a riboswitch variant with mutations in the expression domain that still binds TPP.

Aptamers in the virologists' toolkit

Aptamers are artificial nucleic acid ligands that can be generated in vitro against a wide range of molecules, including the gene products of viruses. Aptamers are isolated from complex libraries of synthetic nucleic acids by an iterative, cell-free process that involves repetitively reducing the complexity of the library by partitioning on the basis of selective binding to the target molecule, followed by reamplification. For virologists, aptamers have potential uses as tools to help to analyse the molecular biology of virus replication, as a complement to the more familiar monoclonal antibodies. They also have potential applications as diagnostic biosensors and in the development of antiviral agents. In recent years, these two promising avenues have been explored increasingly by virologists; here, the progress that has been made is reviewed.

LNA derivatives of a kissing aptamer targeted to the trans-activating responsive RNA element of HIV-1

We previously identified an RNA aptamer targeted to the trans-activating responsive (TAR) element of the HIV-1 genome [F. Ducongé, J.J. Toulmé, In vitro selection identifies key determinants for loop--loop interactions: RNA aptamers selective for the TAR RNA element of HIV--1. RNA 5 (1999) 1605--1614]. This hairpin aptamer binds to its target through loop-loop interactions. We derived chemically modified R06 aptamers that show improved nuclease resistance and affinity for TAR. We review here the results obtained with chimeric aptamers containing locked nucleic acid (LNA) residues. Chimeras containing 2 to 4 LNA residues in an RNA or 2'-O-methyl,RNA context display binding properties of interest and compete with the viral protein Tat for binding to TAR. NMR studies have shown that these properties are modulated by the conformation of the loop-loop helix depending on the presence of LNA residues.

Molecular dynamics simulations of RNA: an in silico single molecule approach

RNA molecules are now known to be involved in the processing of genetic information at all levels, taking on a wide variety of central roles in the cell. Understanding how RNA molecules carry out their biological functions will require an understanding of structure and dynamics at the atomistic level, which can be significantly improved by combining computational simulation with experiment. This review provides a critical survey of the state of molecular dynamics (MD) simulations of RNA, including a discussion of important current limitations of the technique and examples of its successful application. Several types of simulations are discussed in detail, including those of structured RNA molecules and their interactions with the surrounding solvent and ions, catalytic RNAs, and RNA-small molecule and RNA-protein complexes. Increased cooperation between theorists and experimentalists will allow expanded judicious use of MD simulations to complement conceptually related single molecule experiments. Such cooperation will open the door to a fundamental understanding of the structure-function relationships in diverse and complex RNA molecules. .

Stabilities of HIV-1 DIS type RNA loop-loop interactions in vitro and in vivo

RNA loop–loop interactions are a prevalent motif in the formation of tertiary structure and are well suited to trigger molecular recognition between RNA molecules. We determined the stabilities of several loop–loop interactions with a constant 6 bp core sequence and varying unpaired flanking nucleotides and found that the flanking bases have a strong influence on the stability and ion dependence of the kissing complex. In general, the stabilities determined in 1 M Na⁺ are equivalent to those in the presence of near physiological Mg²⁺ concentrations. Therefore we further tested whether the stabilities determined in vitro and within yeast cells correlate, using a recently developed yeast RNA-hybrid system. For the majority of the loop types analyzed here, the melting temperatures determined in vitro are in good agreement with the relative β-galactosidase activity in yeast cells, showing that data derived from in vitro measurements reflect in vivo properties. The most stable interactions are the naturally occurring HIV-1 DIS MAL and LAI derived loops with the motif (5′ A^A/_GN₆A 3′), emphasizing the crucial role of stable kissing complexes in HIV genome dimerization.

A cyclic PNA-based compound targeting domain IV of HCV IRES RNA inhibits in vitro IRES-dependent translation

A cyclic molecule 1 constituted by a hepta-peptide nucleic acid sequence complementary to the apical loop of domain IV of hepatitis C virus (HCV) internal ribosome entry site (IRES) RNA has been prepared via a 'mixed' liquid-phase strategy, which relies on easily available protected PNA and poly(2-aminoethylglycinamide) building blocks. This compound 1 has been elaborated to mimic 'loop-loop' interactions. For comparison, its linear analog has also been investigated. Although preliminary biological assays have revealed the ability of 1 to inhibit in vitro the HCV IRES-dependent translation in a dose-dependent manner, the linear analog has shown a slightly higher activity.

Structural dynamics of precursor and product of the RNA enzyme from the hepatitis delta virus as revealed by molecular dynamics simulations

The hepatitis delta virus (HDV) ribozyme is a self-cleaving RNA enzyme involved in the replication of a human pathogen, the hepatitis delta virus. Recent crystal structures of the precursor and product of self-cleavage, together with detailed kinetic analyses, have led to hypotheses on the catalytic strategies employed by the HDV ribozyme. We report molecular dynamics (MD) simulations (approximately 120 ns total simulation time) to test the plausibility that specific conformational rearrangements are involved in catalysis. Site-specific self-cleavage requires cytidine in position 75 (C75). A precursor simulation with unprotonated C75 reveals a rather weak dynamic binding of C75 in the catalytic pocket with spontaneous, transient formation of a H-bond between U-1(O2') and C75(N3). This H-bond would be required for C75 to act as the general base. Upon protonation in the precursor, C75H+ has a tendency to move towards its product location and establish a firm H-bonding network within the catalytic pocket. However, a C75H+(N3)-G1(O5') H-bond, which would be expected if C75 acted as a general acid catalyst, is not observed on the present simulation timescale. The adjacent loop L3 is relatively dynamic and may serve as a flexible structural element, possibly gated by the closing U20.G25 base-pair, to facilitate a conformational switch induced by a protonated C75H+. L3 also controls the electrostatic environment of the catalytic core, which in turn may modulate C75 base strength and metal ion binding. We find that a distant RNA tertiary interaction involving a protonated cytidine (C41) becomes unstable when left unprotonated, leading to disruptive conformational rearrangements adjacent to the catalytic core. A Na ion temporarily compensates for the loss of the protonated hydrogen bond, which is strikingly consistent with the experimentally observed synergy between low pH and high Na+ concentrations in mediating residual self-cleavage of the HDV ribozyme in the absence of divalents.

Dynamics and Stability of GCAA Tetraloops with 2-Aminopurine and Purine Substitutions

The contributions of various interactions in the GGCGCAAGCC hairpin containing a GCAA tetraloop were studied by computer simulations using the substitutions of functional groups. The guanosine (G) in the first tetraloop position or in the C-G closing base pair was replaced by 2-aminopurine (AP), and the individual tetraloop's adenosines (A) were replaced by purine (PUR). These substitutions eliminated particular hydrogen bonds thought to stabilize the GCAA tetraloop. For each substitution, molecular dynamics (MD) simulations were carried out in an aqueous solution with sodium counterions, using the CHARMM27 force field. The MD simulations showed that the substitutions in the first (G-->AP) and the third (A-->PUR) position of the GCAA tetraloop did not significantly influence the conformation of the hairpin. A long-lived bridging water molecule observed in the GCAA loop was present in both modified loops. The substitutions made in the last loop position (A-->PUR) or in the C-G base pair closing the tetraloop (G-->AP) to some extent influenced the loop structure and dynamics. These loops did not display the long-lived bridging water molecules. When the second A in the GCAA loop was replaced by PUR, the first A in the loop was observed in the anti or in the syn orientation about the glycosyl bond. The G to AP substitution in C-G base pair led to a change of their arrangement from the Watson-Crick to wobble. The MD simulations of the hairpin with C-AP wobble closing base pair showed increased conformational dynamics of the hairpin. The changes of hairpin formation free energy associated with the substitutions of individual bases were calculated by the free energy perturbation method. Our theoretical estimates suggest a larger destabilization for the G to AP substitutions in GCAA loop than for the substitutions of individual A's by PUR, which is in accordance with experimental tendency. The calculations predicted a similar free energy change for G to AP substitutions in the GCAA tetraloop and in the C-G closing base pair.

Closing loop base pairs in RNA loop-loop complexes: structural behavior, interaction energy and solvation analysis through molecular dynamics simulations

Nanosecond molecular dynamics using the Ewald summation method have been performed to elucidate the structural and energetic role of the closing base pair in loop-loop RNA duplexes neutralized by Mg2+ counterions in aqueous phases. Mismatches GA, CU and Watson-Crick GC base pairs have been considered for closing the loop of an RNA in complementary interaction with HIV-1 TAR. The simulations reveal that the mismatch GA base, mediated by a water molecule, leads to a complex that presents the best compromise between flexibility and energetic contributions. The mismatch CU base pair, in spite of the presence of an inserted water molecule, is too short to achieve a tight interaction at the closing-loop junction and seems to force TAR to reorganize upon binding. An energetic analysis has allowed us to quantify the strength of the interactions of the closing and the loop-loop pairs throughout the simulations. Although the water-mediated GA closing base pair presents an interaction energy similar to that found on fully geometry-optimized structure, the water-mediated CU closing base pair energy interaction reaches less than half the optimal value.

Determinants of apical loop–internal loop RNA–RNA interactions involving the HCV IRES

Domain II of the hepatitis C virus internal ribosome entry site is a major RNA structure involved in the viral mRNA translation. It comprises four different structural domains. We performed in vitro selection against the apical loop of the domain II and we identified RNA aptamers folding as an imperfect hairpin with an internal loop of interacting with the apical loop of the domain II. This RNA-RNA interaction creates apical loop-internal loop complex. The aptamer binds the target with an apparent K(d) of 35nM. In this study, the main structural elements of the target and the aptamer involved in the formation of the complex are characterized by mutation, deletion, and RNase probing analysis. We demonstrate that a complementary loop flanked by G,C rich upper and lower stems are crucial for such RNA-RNA interactions.

Molecular dynamics simulations of RNA kissing-loop motifs reveal structural dynamics and formation of cation-binding pockets

Explicit solvent molecular dynamics (MD) simulations were carried out for three RNA kissing-loop complexes. The theoretical structure of two base pairs (2 bp) complex of H3 stem-loop of Moloney murine leukemia virus agrees with the NMR structure with modest violations of few NMR restraints comparable to violations present in the NMR structure. In contrast to the NMR structure, however, MD shows relaxed intermolecular G-C base pairs. The core region of the kissing complex forms a cation-binding pocket with highly negative electrostatic potential. The pocket shows nanosecond-scale breathing motions coupled with oscillations of the whole molecule. Additional simulations were carried out for 6 bp kissing complexes of the DIS HIV-1 subtypes A and B. The simulated structures agree well with the X-ray data. The subtype B forms a novel four-base stack of bulged-out adenines. Both 6 bp kissing complexes have extended cation-binding pockets in their central parts. While the pocket of subtype A interacts with two hexacoordinated Mg2+ ions and one sodium ion, pocket of subtype B is filled with a string of three delocalized Na+ ions with residency times of individual cations 1-2 ns. The 6 bp complexes show breathing motions of the cation-binding pockets and loop major grooves.

Regulating eukaryotic gene expression with aptamers

Aptamers are RNA or DNA oligonucleotides identified within a randomly synthesized library, through an in vitro selection procedure. The selected candidates display a pre-determined property of interest with respect to a given target. Successful selection has been carried out against targets ranging from small (amino acids, antibiotics) to macro-molecules (proteins, nucleic acids). They generally show an affinity in the nanomolar range and a high specificity of target recognition. Interestingly, aptamers selected against purified targets in the test tube retain their properties within cells. RNA aptamers can be generated in situ from an appropriate DNA construct or delivered as nuclease-resistant oligonucleotide analogues. For example, aptamers recognizing RNA structure through loop-loop interactions modulate the trans-activation of in vitro transcription mediated by the TAR RNA element of human immunodeficiency virus type 1. Consequently, they constitute both exquisite tools for functional genomics analysis and promising prototypes of therapeutic agents. Natural aptameric motifs have been identified within mRNA sequences, which upon binding to a metabolite control the expression of the encoded gene, which is generally involved in the biosynthesis of this particular metabolite.

LNA/DNA chimeric oligomers mimic RNA aptamers targeted to the TAR RNA element of HIV-1

One of the major limitations of the use of phosphodiester oligonucleotides in cells is their rapid degradation by nucleases. To date, several chemical modifications have been employed to overcome this issue but insufficient efficacy and/or specificity have limited their in vivo usefulness. In this work conformationally restricted nucleotides, locked nucleic acid (LNA), were investigated to design nuclease resistant aptamers targeted against the HIV-1 TAR RNA. LNA/DNA chimeras were synthesized from a shortened version of the hairpin RNA aptamer identified by in vitro selection against TAR. The results indicate that these modifications confer good protection towards nuclease digestion. Electrophoretic mobility shift assays, thermal denaturation monitored by UV-spectroscopy and surface plasmon resonance experiments identified LNA/DNA TAR ligands that bind to TAR with a dissociation constant in the low nanomolar range as the parent RNA aptamer. The crucial G, A residues that close the aptamer loop remain a key structural determinant for stable LNA/DNA chimera-TAR complexes. This work provides evidence that LNA modifications alternated with DNA can generate stable structured RNA mimics for interacting with folded RNA targets.

Simulation and modeling of nucleic acid structure, dynamics and interactions

In moving towards the simulation of larger nucleic acid assemblies over longer timescales that include more accurate representations of the environment, we are nearing the end of an era characterized by single nanosecond molecular dynamics simulation of nucleic acids. We are excited by the promise and predictability of the modeling methods, yet remain prudently cautious of sampling and force field limitations. Highlights include the accurate representation of subtle drug-DNA interactions, the detailed study of modified and unusual nucleic acid structures, insight into the influence of dynamics on the structure of DNA, and exploration of the interaction of solvent and ions with nucleic acids.

MD Studies of the DIS/DIS Kissing Complex Solution and X-Ray Structures

As in all retroviruses, human immunodeficiency virus (HIV) genomic RNA is packaged into virions as a dimer. The two copies of the genome are noncovalently linked by their 5'-ends in the dimerization initiating site (DIS), which folds as a hairpin containing an apical autocomplementary sequence. In vitro, DIS is able to dimerize in two conformations: a kissing complex and an extended dimer. Both conformations have been resolved by NMR and x-ray diffraction. Here, we report molecular dynamics (MD) studies of the two available structures for the DIS/DIS kissing complex in aqueous solution and in the presence of sodium counterions. The two structures behave in two different manners. On one hand, the NMR structure displays a very stable behavior, and the simulated structure remains very close to the starting structure. On the other hand, the structure issued from crystallography displays a more dynamic behavior, in which residues A8 and A9 are seen in a new and surprising bulge-in conformation. The transition from the bulge-out to the bulge-in conformation is analyzed, and a new and simple dimerization process is proposed.