Evolution in vitro of an RNA enzyme with altered metal dependence

Switch-like compaction of poly(ADP-ribose) upon cation binding

Poly(ADP-ribose) (PAR) is a homopolymer of adenosine diphosphate ribose that is added to proteins as a posttranslational modification to regulate numerous cellular processes. PAR also serves as a scaffold for protein binding in macromolecular complexes, including biomolecular condensates. It remains unclear how PAR achieves specific molecular recognition. Here, we use single-molecule fluorescence resonance energy transfer (smFRET) to evaluate PAR flexibility under different cation conditions. We demonstrate that, compared to RNA and DNA, PAR has a longer persistence length and undergoes a sharper transition from extended to compact states in physiologically relevant concentrations of various cations (Na+, Mg2+, Ca2+, and spermine4+). We show that the degree of PAR compaction depends on the concentration and valency of cations. Furthermore, the intrinsically disordered protein FUS also served as a macromolecular cation to compact PAR. Taken together, our study reveals the inherent stiffness of PAR molecules, which undergo switch-like compaction in response to cation binding. This study indicates that a cationic environment may drive recognition specificity of PAR.

Switch-like Compaction of Poly(ADP-ribose) Upon Cation Binding

Poly(ADP-ribose) (PAR) is a homopolymer of adenosine diphosphate ribose that is added to proteins as a post-translational modification to regulate numerous cellular processes. PAR also serves as a scaffold for protein binding in macromolecular complexes, including biomolecular condensates. It remains unclear how PAR achieves specific molecular recognition. Here, we use single-molecule fluorescence resonance energy transfer (smFRET) to evaluate PAR flexibility under different cation conditions. We demonstrate that, compared to RNA and DNA, PAR has a longer persistence length and undergoes a sharper transition from extended to compact states in physiologically relevant concentrations of various cations (Na ⁺ , Mg ²⁺ , Ca ²⁺ , and spermine). We show that the degree of PAR compaction depends on the concentration and valency of cations. Furthermore, the intrinsically disordered protein FUS also served as a macromolecular cation to compact PAR. Taken together, our study reveals the inherent stiffness of PAR molecules, which undergo switch-like compaction in response to cation binding. This study indicates that a cationic environment may drive recognition specificity of PAR.

Significance

Poly(ADP-ribose) (PAR) is an RNA-like homopolymer that regulates DNA repair, RNA metabolism, and biomolecular condensate formation. Dysregulation of PAR results in cancer and neurodegeneration. Although discovered in 1963, fundamental properties of this therapeutically important polymer remain largely unknown. Biophysical and structural analyses of PAR have been exceptionally challenging due to the dynamic and repetitive nature. Here, we present the first single-molecule biophysical characterization of PAR. We show that PAR is stiffer than DNA and RNA per unit length. Unlike DNA and RNA which undergoes gradual compaction, PAR exhibits an abrupt switch-like bending as a function of salt concentration and by protein binding. Our findings points to unique physical properties of PAR that may drive recognition specificity for its function.

Extended magnesium and calcium force field parameters for accurate ion-nucleic acid interactions in biomolecular simulations

Magnesium and calcium play an essential role in the folding and function of nucleic acids. To correctly describe their interactions with DNA and RNA in biomolecular simulations, an accurate parameterization is crucial. In most cases, the ion parameters are optimized based on a set of experimental solution properties such as solvation free energies, radial distribution functions, water exchange rates, and activity coefficient derivatives. However, the transferability of such bulk-optimized ion parameters to quantitatively describe biomolecular systems is limited. Here, we extend the applicability of our previous bulk-optimized parameters by including experimental binding affinities toward the phosphate oxygen on nucleic acids. In particular, we systematically adjust the combination rules that are an integral part of the pairwise interaction potentials of classical force fields. This allows us to quantitatively describe specific ion binding to nucleic acids without changing the solution properties in the most simple and efficient way. We show the advancement of the optimized Lorentz combination rule for two representative nucleic acid systems. For double-stranded DNA, the optimized combination rule for Ca²⁺ significantly improves the agreement with experiments, while the standard combination rule leads to unrealistically distorted DNA structures. For the add A-riboswitch, the optimized combination rule for Mg²⁺ improves the structure of two specifically bound Mg²⁺ ions as judged by the experimental distance to the binding site. Including experimental binding affinities toward specific ion binding sites on biomolecules, therefore, provides a promising perspective to develop a more accurate description of metal cations for biomolecular simulations.

Evolution of the Standard Genetic Code

A near-universal Standard Genetic Code (SGC) implies a single origin for present Earth life. To study this unique event, I compute paths to the SGC, comparing different plausible histories. Notably, SGC-like coding emerges from traditional evolutionary mechanisms, and a superior route can be identified. To objectively measure evolution, progress values from 0 (random coding) to 1 (SGC-like) are defined: these measure fractions of random-code-to-SGC distance. Progress types are spacing/distance/delta Polar Requirement, detecting space between identical assignments/mutational distance to the SGC/chemical order, respectively. The coding system is based on selected RNAs performing aminoacyl-RNA synthetase reactions. Acceptor RNAs exhibit SGC-like Crick wobble; alternatively, non-wobbling triplets uniquely encode 20 amino acids/start/stop. Triplets acquire 22 functions by stereochemistry, selection, coevolution, or at random. Assignments also propagate to an assigned triplet’s neighborhood via single mutations, but can also decay. A vast code universe makes futile evolutionary paths plentiful. Thus, SGC evolution is critically sensitive to disorder from random assignments. Evolution also inevitably slows near coding completion. The SGC likely avoided these difficulties, and two suitable paths are compared. In late wobble, a majority of non-wobble assignments are made before wobble is adopted. In continuous wobble, a uniquely advantageous early intermediate yields an ordered SGC. Revised coding evolution (limited randomness, late wobble, concentration on amino acid encoding, chemically conservative coevolution with a chemically ordered elite) produces varied full codes with excellent joint progress values. A population of only 600 independent coding tables includes SGC-like members; a Bayesian path toward more accurate SGC evolution is available.

A combinatorial method to isolate short ribozymes from complex ribozyme libraries

In vitro selections are the only known methods to generate catalytic RNAs (ribozymes) that do not exist in nature. Such new ribozymes are used as biochemical tools, or to address questions on early stages of life. In both cases, it is helpful to identify the shortest possible ribozymes since they are easier to deploy as a tool, and because they are more likely to have emerged in a prebiotic environment. One of our previous selection experiments led to a library containing hundreds of different ribozyme clusters that catalyze the triphosphorylation of their 5'-terminus. This selection showed that RNA systems can use the prebiotically plausible molecule cyclic trimetaphosphate as an energy source. From this selected ribozyme library, the shortest ribozyme that was previously identified had a length of 67 nucleotides. Here we describe a combinatorial method to identify short ribozymes from libraries containing many ribozymes. Using this protocol on the library of triphosphorylation ribozymes, we identified a 17-nucleotide sequence motif embedded in a 44-nucleotide pseudoknot structure. The described combinatorial approach can be used to analyze libraries obtained by different in vitro selection experiments.

Construction and characterization of an improved DNA-launched infectious clone of duck hepatitis a virus type 1

Background

DNA-launched infectious system is a useful tool with high rescue efficiency that allows the introduction of mutations in specific positions to investigate the function of an individual viral element. Rescued virus particles could be harvested by directly transfecting the DNA-launched recombinant plasmid to the host cells, which will reduce labor and experimental cost by skipping the in vitro transcription assay.

Methods

A total of four overlapping fragments covering the entire viral genome were amplified and then were assembled into a transformation vector based on pIRES2-EGFP to establish the DNA-launched infectious system of duck hepatitis A virus type 1 (DHAV-1), named pIR-DHAV-1. Reverse transcription polymerase chain reaction (RT-PCR) detection, quantitative real-time polymerase chain reaction (qRT-PCR), western blotting assay and indirect immunofluorescence (IFA) were conducted for rescued virus identification. A total of 4.0 μg of recombinant plasmid of pIR-DHAV-1 and in vitro transcribed product of 4.0 μg of RNA-launched infectious clone named pR-DHAV-1 were transfected into BHK-21 cells to analyze the rescue efficiency. Following that, tissue tropism of rescued virus (rDHAV-1) and parental virus (pDHAV-1) were assayed for virulence testing in 1-day-old ducklings.

Results

Rescued virus particles carry the designed genetic marker which could be harvested by directly transfecting pIR-DHAV-1 to BHK-21 cells. The qRT-PCR and western blotting results indicated that rDHAV-1 shared similar growth characteristics with pDHAV-1. Furthermore, DNA-launched infectious system possessed much higher rescue efficiency assay compared to RNA-launched infectious system. The mutation at position 3042 from T to C has no impact on viral replication and tissue tropism. From 1 h post infection (hpi) to 48 hpi, the viral RNA copies of rDHAV-1 in liver were the highest among the six tested tissues (with an exception of thymus at 6 hpi), while the viral RNA copy numbers in heart and kidney were alternately the lowest.

Conclusion

We have constructed a genetically stable and highly pathogenic DNA-launched infectious clone, from which the rescued virus could be harvested by direct transfection with recombinant plasmids. rDHAV-1 shared similar growth characteristics and tissue tropism with pDHAV-1. The DNA-launched infectious system of DHAV-1 possessed higher rescue efficiency compared to the traditional RNA-launched infectious system.

Analysis of in vitro evolution reveals the underlying distribution of catalytic activity among random sequences

The emergence of catalytic RNA is believed to have been a key event during the origin of life. Understanding how catalytic activity is distributed across random sequences is fundamental to estimating the probability that catalytic sequences would emerge. Here, we analyze the in vitro evolution of triphosphorylating ribozymes and translate their fitnesses into absolute estimates of catalytic activity for hundreds of ribozyme families. The analysis efficiently identified highly active ribozymes and estimated catalytic activity with good accuracy. The evolutionary dynamics follow Fisher's Fundamental Theorem of Natural Selection and a corollary, permitting retrospective inference of the distribution of fitness and activity in the random sequence pool for the first time. The frequency distribution of rate constants appears to be log-normal, with a surprisingly steep dropoff at higher activity, consistent with a mechanism for the emergence of activity as the product of many independent contributions.

A divalent cation-dependent variant of the glmS ribozyme with stringent Ca2+ selectivity co-opts a preexisting nonspecific metal ion-binding site

Ribozymes use divalent cations for structural stabilization, as catalytic cofactors, or both. Because of the prominent role of Ca²⁺ in intracellular signaling, engineered ribozymes with stringent Ca²⁺ selectivity would be important in biotechnology. The wild-type glmS ribozyme (glmS^WT) requires glucosamine-6-phosphate (GlcN6P) as a catalytic cofactor. Previously, a glmS ribozyme variant with three adenosine mutations (glmS^AAA) was identified, which dispenses with GlcN6P and instead uses, with little selectivity, divalent cations as cofactors for site-specific RNA cleavage. We now report a Ca²⁺-specific ribozyme (glmS^Ca) evolved from glmS^AAA that is >10,000 times more active in Ca²⁺ than Mg²⁺, is inactive in even 100 mM Mg²⁺, and is not responsive to GlcN6P. This stringent selectivity, reminiscent of the protein nuclease from Staphylococcus, allows rapid and selective ribozyme inactivation using a Ca²⁺ chelator such as EGTA. Because glmS^Ca functions in physiologically relevant Ca²⁺ concentrations, it can form the basis for intracellular sensors that couple Ca²⁺ levels to RNA cleavage. Biochemical analysis of glmS^Ca reveals that it has co-opted for selective Ca²⁺ binding a nonspecific cation-binding site responsible for structural stabilization in glmS^WT and glmS^AAA Fine-tuning of the selectivity of the cation site allows repurposing of this preexisting molecular feature.

The X-ray Structures of Six Octameric RNA Duplexes in the Presence of Different Di- and Trivalent Cations

Due to the polyanionic nature of RNA, the principles of charge neutralization and electrostatic condensation require that cations help to overcome the repulsive forces in order for RNA to adopt a three-dimensional structure. A precise structural knowledge of RNA-metal ion interactions is crucial to understand the mechanism of metal ions in the catalytic or regulatory activity of RNA. We solved the crystal structure of an octameric RNA duplex in the presence of the di- and trivalent metal ions Ca(2+), Mn(2+), Co(2+), Cu(2+), Sr(2+), and Tb(3+). The detailed investigation reveals a unique innersphere interaction to uracil and extends the knowledge of the influence of metal ions for conformational changes in RNA structure. Furthermore, we could demonstrate that an accurate localization of the metal ions in the X-ray structures require the consideration of several crystallographic and geometrical parameters as well as the anomalous difference map.

In vitro evolution of distinct self-cleaving ribozymes in diverse environments

In vitro evolution experiments have long been used to evaluate the roles of RNA in both modern and ancient biology, and as a tool for biotechnology applications. The conditions under which these experiments have been conducted, however, do not reflect the range of cellular environments in modern biology or our understanding of chemical environments on the early earth, when the atmosphere and oceans were largely anoxic and soluble Fe²⁺ was abundant. To test the impact of environmental factors relevant to RNA's potential role in the earliest forms of life, we evolved populations of self-cleaving ribozymes in an anoxic atmosphere with varying pH in the presence of either Fe²⁺ or Mg²⁺. Populations evolved under these different conditions are dominated by different sequences and secondary structures, demonstrating global differences in the underlying fitness landscapes. Comparisons between evolutionary outcomes and catalytic activities also indicate that Mg²⁺ can readily take the place of Fe²⁺ in supporting the catalysis of RNA cleavage at neutral pH, but not at lower pH. These results highlight the importance of considering the specific environments in which functional biopolymers evolve when evaluating their potential roles in the origin of life, extant biology, or biotechnology.

From formamide to RNA, the path is tenuous but continuous

Reactions of formamide (NH2COH) in the presence of catalysts of both terrestrial and meteoritic origin yield, in plausible and variegated conditions, a large panel of precursors of (pre)genetic and (pre)metabolic interest. Formamide chemistry potentially satisfies all of the steps from the very initial precursors to RNA. Water chemistry enters the scene in RNA non-enzymatic synthesis and recombination.

Empirical demonstration of environmental sensing in catalytic RNA: evolution of interpretive behavior at the origins of life

Background

The origins of life on the Earth required chemical entities to interact with their environments in ways that could respond to natural selection. The concept of interpretation, where biotic entities use signs in their environment as proxy for the existence of other items of selective value in their environment, has been proposed on theoretical grounds to be relevant to the origins and early evolution of life. However this concept has not been demonstrated empirically.

Results

Here, we present data that certain catalytic RNA sequences have properties that would enable interpretation of divalent cation levels in their environment. By assaying the responsiveness of two variants of the Tetrahymena ribozyme to the Ca²⁺ ion as a sign for the more catalytically useful Mg²⁺ ion, we show an empirical proof-of-principle that interpretation can be an evolvable trait in RNA, often suggested as a model system for early life. In particular we demonstrate that in vitro, the wild-type version of the Tetrahymena ribozyme is not interpretive, in that it cannot use Ca²⁺ as a sign for Mg²⁺. Yet a variant of this sequence containing five mutations that alter its ability to utilize the Ca²⁺ ion engenders a strong interpretive characteristic in this RNA.

Conclusions

We have shown that RNA molecules in a test tube can meet the minimum criteria for the evolution of interpretive behaviour in regards to their responses to divalent metal ion concentrations in their environment. Interpretation in RNA molecules provides a property entirely dependent on natural physico-chemical interactions, but capable of shaping the evolutionary trajectory of macromolecules, especially in the earliest stages of life’s history.

Aptamers as both drugs and drug-carriers

Aptamers are short nucleic acid oligos. They may serve as both drugs and drug-carriers. Their use as diagnostic tools is also evident. They can be generated using various experimental, theoretical, and computational techniques. The systematic evolution of ligands by exponential enrichment which uses iterative screening of nucleic acid libraries is a popular experimental technique. Theory inspired methodology entropy-based seed-and-grow strategy that designs aptamer templates to bind specifically to targets is another one. Aptamers are predicted to be highly useful in producing general drugs and theranostic drugs occasionally for certain diseases like cancer, Alzheimer's disease, and so on. They bind to various targets like lipids, nucleic acids, proteins, small organic compounds, and even entire organisms. Aptamers may also serve as drug-carriers or nanoparticles helping drugs to get released in specific target regions. Due to better target specific physical binding properties aptamers cause less off-target toxicity effects. Therefore, search for aptamer based drugs, drug-carriers, and even diagnostic tools is expanding fast. The biophysical properties in relation to the target specific binding phenomena of aptamers, energetics behind the aptamer transport of drugs, and the consequent biological implications will be discussed. This review will open up avenues leading to novel drug discovery and drug delivery.

Calcium inhibition of ribonuclease H1 two-metal ion catalysis

Most phosphate-processing enzymes require Mg(2+) as a cofactor to catalyze nucleotide cleavage and transfer reactions. Ca(2+) ions inhibit many of these enzymatic activities, despite Ca(2+) and Mg(2+) having comparable binding affinities and overall biological abundances. Here we study the molecular details of the calcium inhibition mechanism for phosphodiester cleavage, an essential reaction in the metabolism of nucleic acids and nucleotides, by comparing Ca(2+)- and Mg(2+) catalyzed reactions. We study the functional roles of the specific metal ion sites A and B in enabling the catalytic cleavage of an RNA/DNA hybrid substrate by B. halodurans ribonuclease (RNase) H1 using hybrid quantum-mechanics/molecular mechanics (QM/MM) free energy calculations. We find that Ca(2+) substitution of either of the two active-site Mg(2+) ions substantially increases the height of the reaction barrier and thereby abolishes the catalytic activity. Remarkably, Ca(2+) at the A site is inactive also in Mg(2+)-optimized active-site structures along the reaction path, whereas Mg(2+) substitution recovers activity in Ca(2+)-optimized structures. Geometric changes resulting from Ca(2+) substitution at metal ion site A may thus be a secondary factor in the loss of catalytic activity. By contrast, at metal ion site B geometry plays a more important role, with only a partial recovery of activity after Mg(2+) substitution in Ca(2+)-optimized structures. Ca(2+)-substitution also leads to a change in mechanism, with deprotonation of the water nucleophile requiring a closer approach to the scissile phosphate, which in turn increases the barrier. As a result, Ca(2+) is less efficient in activating the water. As a likely cause for the different reactivities of Mg(2+) and Ca(2+) ions in site A, we identify differences in charge transfer to the ions and the associated decrease in the pKa of the oxygen nucleophile attacking the phosphate group.

Next-generation sequencing reveals how RNA catalysts evolve from random space

Catalytic RNAs are attractive objects for studying molecular evolution. To understand how RNA libraries can evolve from randomness toward highly active catalysts, we analyze the original samples that led to the discovery of Diels-Alderase ribozymes by next-generation sequencing. Known structure-activity relationships are used to correlate abundance with catalytic performance. We find that efficient catalysts arose not just from selection for reactivity among the members of the starting library, but from improvement of less potent precursors by mutations. We observe changes in the ribozyme population in response to increasing selection pressure. Surprisingly, even after many rounds of enrichment, the libraries are highly diverse, suggesting that potential catalysts are more abundant in random space than generally thought. To highlight the use of next-generation sequencing as a tool for in vitro selections, we also apply this technique to a recent, less characterized ribozyme selection. Making use of the correlation between sequence evolution and catalytic activity, we predict mutations that improve ribozyme activity and validate them biochemically. Our study reveals principles underlying ribozyme in vitro selections and provides guidelines to render future selections more efficient, as well as to predict the conservation of key structural elements, allowing the rational improvement of catalysts.

Low selection pressure aids the evolution of cooperative ribozyme mutations in cells

Understanding the evolution of functional RNA molecules is important for our molecular understanding of biology. Here we tested experimentally how two evolutionary parameters, selection pressure and recombination, influenced the evolution of an evolving RNA population. This was done using four parallel evolution experiments that employed low or gradually increasing selection pressure, and recombination events either at the end or dispersed throughout the evolution. As model system, a trans-splicing group I intron ribozyme was evolved in Escherichia coli cells over 12 rounds of selection and amplification, including mutagenesis and recombination. The low selection pressure resulted in higher efficiency of the evolved ribozyme populations, whereas differences in recombination did not have a strong effect. Five mutations were responsible for the highest efficiency. The first mutation swept quickly through all four evolving populations, whereas the remaining four mutations accumulated later and more efficiently under low selection pressure. To determine why low selection pressure aided this evolution, all evolutionary intermediates between the wild type and the 5-mutation variant were constructed, and their activities at three different selection pressures were determined. The resulting fitness profiles showed a high cooperativity among the four late mutations, which can explain why high selection pressure led to inefficient evolution. These results show experimentally how low selection pressure can benefit the evolution of cooperative mutations in functional RNAs.

Enhanced group II intron retrohoming in magnesium-deficient Escherichia coli via selection of mutations in the ribozyme core

Mobile group II introns are bacterial retrotransposons thought to be evolutionary ancestors of spliceosomal introns and retroelements in eukaryotes. They consist of a catalytically active intron RNA ("ribozyme") and an intron-encoded reverse transcriptase, which function together to promote RNA splicing and intron mobility via reverse splicing of the intron RNA into new DNA sites ("retrohoming"). Although group II introns are active in bacteria, their natural hosts, they function inefficiently in eukaryotes, where lower free Mg(2+) concentrations decrease their ribozyme activity and constitute a natural barrier to group II intron proliferation within nuclear genomes. Here, we show that retrohoming of the Ll.LtrB group II intron is strongly inhibited in an Escherichia coli mutant lacking the Mg(2+) transporter MgtA, and we use this system to select mutations in catalytic core domain V (DV) that partially rescue retrohoming at low Mg(2+) concentrations. We thus identified mutations in the distal stem of DV that increase retrohoming efficiency in the MgtA mutant up to 22-fold. Biochemical assays of splicing and reverse splicing indicate that the mutations increase the fraction of intron RNA that folds into an active conformation at low Mg(2+) concentrations, and terbium-cleavage assays suggest that this increase is due to enhanced Mg(2+) binding to the distal stem of DV. Our findings indicate that DV is involved in a critical Mg(2+)-dependent RNA folding step in group II introns and demonstrate the feasibility of selecting intron variants that function more efficiently at low Mg(2+) concentrations, with implications for evolution and potential applications in gene targeting.

Complex formation of cadmium with sugar residues, nucleobases, phosphates, nucleotides, and nucleic acids

Cadmium(II), commonly classified as a relatively soft metal ion, prefers indeed aromatic-nitrogen sites (e.g., N7 of purines) over oxygen sites (like sugar-hydroxyl groups). However, matters are not that simple, though it is true that the affinity of Cd(2+) towards ribose-hydroxyl groups is very small; yet, a correct orientation brought about by a suitable primary binding site and a reduced solvent polarity, as it is expected to occur in a folded nucleic acid, may facilitate metal ion-hydroxyl group binding very effectively. Cd(2+) prefers the guanine(N7) over the adenine(N7), mainly because of the steric hindrance of the (C6)NH(2) group in the adenine residue. This Cd(2+)-(N7) interaction in a guanine moiety leads to a significant acidification of the (N1)H meaning that the deprotonation reaction occurs now in the physiological pH range. N3 of the cytosine residue, together with the neighboring (C2)O, is also a remarkable Cd(2+) binding site, though replacement of (C2)O by (C2)S enhances the affinity towards Cd(2+) dramatically, giving in addition rise to the deprotonation of the (C4)NH(2) group. The phosphodiester bridge is only a weak binding site but the affinity increases further from the mono- to the di- and the triphosphate. The same also holds for the corresponding nucleotides. Complex stability of the pyrimidine-nucleotides is solely determined by the coordination tendency of the phosphate group(s), whereas in the case of purine-nucleotides macrochelate formation takes place by the interaction of the phosphate-coordinated Cd(2+) with N7. The extents of the formation degrees of these chelates are summarized and the effect of a non-bridging sulfur atom in a thiophosphate group (versus a normal phosphate group) is considered. Mixed ligand complexes containing a nucleotide and a further mono- or bidentate ligand are covered and it is concluded that in these species N7 is released from the coordination sphere of Cd(2+). In the case that the other ligand contains an aromatic residue (e.g., 2,2'-bipyridine or the indole ring of tryptophanate) intramolecular stack formation takes place. With buffers like Tris or Bistris mixed ligand complexes are formed. Cd(2+) coordination to dinucleotides and to dinucleoside monophosphates provides some insights regarding the interaction between Cd(2+) and nucleic acids. Cd(2+) binding to oligonucleotides follows the principles of coordination to its units. The available crystal studies reveal that N7 of purines is the prominent binding site followed by phosphate oxygens and other heteroatoms in nucleic acids. Due to its high thiophilicity, Cd(2+) is regularly used in so-called thiorescue experiments, which lead to the identification of a direct involvement of divalent metal ions in ribozyme catalysis.

Characterization of metal ion-nucleic acid interactions in solution

Metal ions are inextricably involved with nucleic acids due to their polyanionic nature. In order to understand the structure and function of RNAs and DNAs, one needs to have detailed pictures on the structural, thermodynamic, and kinetic properties of metal ion interactions with these biomacromolecules. In this review we first compile the physicochemical properties of metal ions found and used in combination with nucleic acids in solution. The main part then describes the various methods developed over the past decades to investigate metal ion binding by nucleic acids in solution. This includes for example hydrolytic and radical cleavage experiments, mutational approaches, as well as kinetic isotope effects. In addition, spectroscopic techniques like EPR, lanthanide(III) luminescence, IR and Raman as well as various NMR methods are summarized. Aside from gaining knowledge about the thermodynamic properties on the metal ion-nucleic acid interactions, especially NMR can be used to extract information on the kinetics of ligand exchange rates of the metal ions applied. The final section deals with the influence of anions, buffers, and the solvent permittivity on the binding equilibria between metal ions and nucleic acids. Little is known on some of these aspects, but it is clear that these three factors have a large influence on the interaction between metal ions and nucleic acids.

Cryptic genetic variation promotes rapid evolutionary adaptation in an RNA enzyme

Cryptic variation is caused by the robustness of phenotypes to mutations. Cryptic variation has no effect on phenotypes in a given genetic or environmental background, but it can have effects after mutations or environmental change. Because evolutionary adaptation by natural selection requires phenotypic variation, phenotypically revealed cryptic genetic variation may facilitate evolutionary adaptation. This is possible if the cryptic variation happens to be pre-adapted, or "exapted", to a new environment, and is thus advantageous once revealed. However, this facilitating role for cryptic variation has not been proven, partly because most pertinent work focuses on complex phenotypes of whole organisms whose genetic basis is incompletely understood. Here we show that populations of RNA enzymes with accumulated cryptic variation adapt more rapidly to a new substrate than a population without cryptic variation. A detailed analysis of our evolving RNA populations in genotype space shows that cryptic variation allows a population to explore new genotypes that become adaptive only in a new environment. Our observations show that cryptic variation contains new genotypes pre-adapted to a changed environment. Our results highlight the positive role that robustness and epistasis can have in adaptive evolution.

HAPIscreen, a method for high-throughput aptamer identification

Background

Aptamers are oligonucleotides displaying specific binding properties for a predetermined target. They are selected from libraries of randomly synthesized candidates through an in vitro selection process termed SELEX (Systematic Evolution of Ligands by EXponential enrichment) alternating selection and amplification steps. SELEX is followed by cloning and sequencing of the enriched pool of oligonucleotides to enable comparison of the selected sequences. The most represented candidates are then synthesized and their binding properties are individually evaluated thus leading to the identification of aptamers. These post-selection steps are time consuming and introduce a bias to the expense of poorly amplified binders that might be of high affinity and are consequently underrepresented. A method that would circumvent these limitations would be highly valuable.

Results

We describe a novel homogeneous solution-based method for screening large populations of oligonucleotide candidates generated from SELEX. This approach, based on the AlphaScreen^® technology, is carried out on the exclusive basis of the binding properties of the selected candidates without the needs of performing a priori sequencing. It therefore enables the functional identification of high affinity aptamers. We validated the HAPIscreen (High throughput APtamer Identification screen) methodology using aptamers targeted to RNA hairpins, previously identified in our laboratory. We then screened pools of candidates issued from SELEX rounds in a 384 well microplate format and identify new RNA aptamers to pre-microRNAs.

Conclusions

HAPIscreen, an Alphascreen^®-based methodology for the identification of aptamers is faster and less biased than current procedures based on sequence comparison of selected oligonucleotides and sampling binding properties of few individuals. Moreover this methodology allows for screening larger number of candidates. Used here for selecting anti-premiR aptamers, HAPIscreen can be adapted to any type of tagged target and is fully amenable to automation.

In vitro selection of RNA aptamers to a protein target by filter immobilization

This unit describes the selection of aptamers from a pool of single-stranded RNA by binding to a protein target. Aptamers generated from this selection experiment can potentially act as protein function inhibitors, and may find applications as therapeutic or diagnostic reagents. A pool of dsDNA is used to generate an ssRNA pool, which is mixed with the protein target. Bound complexes are separated from unbound reagents by filtration, and the RNA:protein complexes are amplified by a combination of reverse transcription, PCR, and in vitro transcription.

The continuous evolution in vitro technique

In vitro experimentation techniques were developed in response to the necessity of exploring new molecular structures and functions and to better understand evolutionary phenomena that shape organismal and molecular populations. The advancement of these techniques has allowed further exploration of more complicated evolutionary dynamics. One such technique is the continuous evolution in vitro (CE) method, to which this unit is devoted. The CE method is characterized by continuous cycles of amplification of RNA molecules that occur without much participation of the researcher. This feature allows us to evolve lineages in which the evolutionary phenomena occurring at the molecular level more closely mimic what happens in organismal populations in the present, or what may have happened in RNA populations during the RNA world stage of life.

Natural and artificial RNAs occupy the same restricted region of sequence space

Different chemical and mutational processes within genomes give rise to sequences with different compositions and perhaps different capacities for evolution. The evolution of functional RNAs may occur on a "neutral network" in which sequences with any given function can easily mutate to sequences with any other. This neutral network hypothesis is more likely if there is a particular region of composition that contains sequences that are functional in general, and if many different functions are possible within this preferred region of composition. We show that sequence preferences in active sites recovered by in vitro selection combine with biophysical folding rules to support the neutral network hypothesis. These simple active-site specifications and folding preferences obtained by artificial selection experiments recapture the previously observed purine bias and specific spread along the GC axis of naturally occurring aptamers and ribozymes isolated from organisms, although other types of RNAs, such as miRNA precursors and spliceosomal RNAs, that act primarily through complementarity to other amino acids do not share these preferences. These universal evolved sequence features are therefore intrinsic in RNA molecules that bind small-molecule targets or catalyze reactions.

In vitro selection of RNA aptamers to a protein target by filter immobilization

This unit describes the selection of aptamers from a pool of single-stranded RNA by binding to a protein target. Aptamers generated from this selection experiment can potentially act as protein function inhibitors, and may find applications as therapeutic or diagnostic reagents. A pool of dsDNA is used to generate an ssRNA pool, which is mixed with the protein target. Bound complexes are separated from unbound reagents by filtration, and the RNA:protein complexes are amplified by a combination of reverse transcription, PCR, and in vitro transcription.

A genotype-to-phenotype map of in vitro selected RNA-cleaving DNAzymes: implications for accessing the target phenotype

Herein, we describe a case study into the population dynamics of in vitro selection, using RNA-cleaving DNAzymes as a model system. We sought to understand how the composition of the population can change over time in response to different levels of selection pressure, and how well these changes are correlated with selection of the target phenotype. The model population is composed of 857 DNAzyme clones representing 215 discrete sequence classes, which had previously been identified from two parallel selection experiments, conducted under an increasingly stringent, or permissive and constant selection time pressure. In this report, we determined the principal phenotypic properties (i.e. k_obs, maximum cleavage yield and PCR efficiency) from a sample of 58 clones representing 46 different major and minor sequence classes from various rounds of each selection experiment. Interestingly, a positive correlation between the catalytic rate constant and the corresponding frequency and temporal position of a given DNAzyme was not consistently observed; however, the strength of the correlation was qualitatively higher under conditions of more stringent selection time pressure. These results suggest that the selective sampling paradigm on which in vitro selection is based, may underestimate the true functional capacity of any given random-sequence library.

In vitro selection of RNA aptamers to a protein target by filter immobilization

This unit describes the selection of aptamers from a pool of single-stranded RNA by binding to a protein target. Aptamers generated from this selection experiment can potentially function as protein inhibitors, and may find applications as therapeutic or diagnostic reagents. A pool of dsDNA is used to generate a ssRNA pool, which is mixed with the protein target. Bound complexes are separated from unbound reagents by filtration, and the RNA:protein complexes are amplified by a combination of reverse transcription, PCR, and in vitro transcription.

In vitro selection of RNA aptamers to a protein target by filter immobilization

In vitro selection of RNA aptamers that bind to a protein target is detailed in the protocols presented in this unit. Aptamers generated from these types of selection experiments can potentially function as protein inhibitors, and are often used as diagnostic or therapeutic reagents.

Accommodation of Ca(II) ions for catalytic activity by a group I ribozyme

The wildtype Tetrahymena ribozyme cannot catalyze detectable levels of phosphotransfer activity in vitro on an exogenous RNA substrate oligonucleotide when calcium(II) is supplied as the only available divalent ion. Nevertheless, low-error mutants of this ribozyme have been acquired through directed evolution that do have activity in 10mM CaCl(2). The mechanisms for such Ca(II) accommodation are not known. Here, we assayed the entire molecule in an effort to identify the roles of the mutations in allowing catalytic activity in Ca(II). We used four biochemical probing techniques - native-gel electrophoresis, hydroxyl radical footprinting, terbium(III) cleavage footprinting, and phosphorothioate interference mapping - to compare the solution structure of the wildtype ribozyme with that of a Ca(II)-active five-site mutant. We compared the gross folding patterns and specific metal-binding sites in both MgCl(2) and CaCl(2) solutions. We detected no large-scale folding differences between the two RNAs in either metal. However, we did discover a limited number of local folding differences, involving regions of the RNA affected by positions 42, 188, and 270. These data support the notion that Ca(II) is accommodated by the Tetrahymena ribozyme by a slight breathing at the active site, but that alterations at, near to, and distal from the active site can all contribute to Ca(II)-based activity.

Selection of fluorescent aptamer beacons that light up in the presence of zinc

In order to generate nucleic acid biosensors that could undergo a reversible conformation change in the presence of the metal zinc, a random sequence pool of single-stranded DNA was immobilized on an oligonucleotide affinity column. In the presence of zinc, those species that underwent a conformational change were released from the column, collected, and amplified. A series of negative and positive selections refined the metal specificity of the selected aptamer beacons. Since the aptamer beacons contained a fluorophore, while the bound oligonucleotide contained a quencher, zinc binding also resulted in an increase in fluorescence. One of the selected beacons, Zn-6m2, bound zinc in the low micromolar range, gave a dose-dependent fluorescence signal, and showed an approximately sixfold increase in fluorescence on zinc binding. While some cross-reactivity with cadmium was observed, it should nonetheless prove possible to use the novel selection method to generate and tune the specificity of a variety of reversible metal biosensors. Such biosensors could potentially be used for continuous monitoring of metals in environmental samples.

Evidence that binding of C5 protein to P RNA enhances ribozyme catalysis by influencing active site metal ion affinity

The RNA subunit (P RNA) of the bacterial RNase P ribonucleoprotein is a ribozyme that catalyzes the Mg-dependent hydrolysis of pre-tRNA, but it requires an essential protein cofactor (P protein) in vivo that enhances substrate binding affinities and catalytic rates in a substrate dependent manner. Previous studies of Bacillus subtilis RNase P, containing a Type B RNA subunit, showed that its cognate protein subunit increases the affinity of metal ions important for catalysis, but the functional role of these ions is unknown. Here, we demonstrate that the Mg2+ dependence of the catalytic step for Escherichia coli RNase P, which contains a more common Type A RNA subunit, is also modulated by its cognate protein subunit (C5), indicating that this property is fundamental to P protein. To monitor specifically the binding of active site metal ions, we analyzed quantitatively the rescue by Cd2+ of an inhibitory Rp phosphorothioate modification at the pre-tRNA cleavage site. The results show that binding of C5 protein increases the apparent affinity of the rescuing Cd2+, providing evidence that C5 protein enhances metal ion affinity in the active site, and thus is likely to contribute significantly to rate enhancement at physiological metal ion concentrations.

Evolution of High-Branching Deoxyribozymes from a Catalytic DNA with a Three-Way Junction

Here, we report the evolution of two star-shaped (five-way junction) deoxyribozymes from a catalytic DNA containing a three-way junction scaffold. The transition was shown to be a switch rather than a gradual progression. The star-shaped motifs, surprisingly, only took five selection cycles to be detected, and another four to dominate the evolving population. Chemical probing experiments indicated that the two deoxyribozymes belong to the same family despite noticeable variations in both the primary sequence and the secondary structure. Our findings not only describe the evolution of high-branching nucleic acid structures from a low-branching catalytic module, but they also illustrate the idea of deriving a rare structural motif by sampling the sequence variants of a given functional nucleic acid.

Unusual metal specificity and structure of the group I ribozyme from Chlamydomonas reinhardtii 23S rRNA

Group I intron ribozymes require cations for folding and catalysis, and the current literature indicates that a number of cations can promote folding, but only Mg2+ and Mn2+ support both processes. However, some group I introns are active only with Mg2+, e.g. three of the five group I introns in Chlamydomonas reinhardtii. We have investigated one of these ribozymes, an intron from the 23S LSU rRNA gene of Chlamydomonas reinhardtii (Cr.LSU), by determining if the inhibition by Mn2+ involves catalysis, folding, or both. Kinetic analysis of guanosine-dependent cleavage by a Cr.LSU ribozyme, 23S.5 Delta Gb, that lacks the 3' exon and intron-terminal G shows that Mn2+ does not affect guanosine binding or catalysis, but instead promotes misfolding of the ribozyme. Surprisingly, ribozyme misfolding induced by Mn2+ is highly cooperative, with a Hill coefficient larger than that of native folding induced by Mg2+. At lower Mn2+ concentrations, metal inhibition is largely alleviated by the guanosine cosubstrate (GMP). The concentration dependence of guanosine cosubstrate-induced folding suggests that it functions by interacting with the G binding site, perhaps by displacing an inhibitory Mn2+. Because of these and other properties of Cr.LSU, the tertiary structure of the intron from 23S.5 Delta Gb was examined using Fe2+-EDTA cleavage. The ground-state structure shows evidence of an unusually open ribozyme core: the catalytic P3-P7 domain and the nucleotides that connect it to the P4-P5-P6 domain are exposed to solvent. The implications of this structure for the in vitro and in vivo properties of this intron ribozyme are discussed.

Calcium(II)-dependent catalytic activity of the Azoarcus ribozyme: testing the limits of resolution for in vitro selection

Group I intron ribozymes isolated from natural sources have a strict dependence on the divalent metal cations Mg(II) or Mn(II) for catalytic activity. However, mutant versions of the Tetrahymena ribozyme have been previously isolated in the laboratory that show demonstrable activity in 10 mM CaCl(2) as the only supplied salt. Here, we sought to discover similar variants of another group I intron that is likely more evolutionarily specialized. We used in vitro selection to isolate a Ca(II)-dependent variant of the naturally-occurring form of the Azoarcus ribozyme, which is half the size of the Tetrahymena ribozyme and possesses an extremely high G+C content of 71%. A mutation of G to A at position 118 was selected in multiple independent trials. Activity of the mutant is very poor in Ca(II) and can only be observed after RT-PCR, highlighting the power of in vitro selection to isolate molecules with rare and low-level activities. The mutation likely confers an alternate but rare folded conformation that permits accommodation of Ca(II) ions and catalysis. This work also serves to caution that although a selection may be successful, isolates may not be catalytically proficient enough to provide useful levels of activity.

Detection of high levels of recombination generated during PCR amplification of RNA templates

Recombination during the PCR amplification of DNA templates can be a serious problem for those seeking to genotype heterogeneous populations, yet a boon to those seeking to enhance variation during in vitro evolution. Here, the extent to which PCR generates chimeric full-length products was estimated using a powerful restriction fragment-length polymorphism (RFLP) assay involving the use of fluorescently labeled PCR primers. Three different RNA-encoding DNA templates were assayed: (i) one for a group I ribozyme, (ii) one for a 16S ribosomal RNA (rRNA), and (iii) one for a messenger RNA (mRNA). In all cases, the observed frequency of chimeric PCR products exceeded 20%, and longer templates appear to produce more chimeric products. Although two of these templates have the potential to form secondary structures during the PCR, this tendency does not seem to heighten recombination frequency. These results corroborate previous studies that show that the production of chimeras can be best attenuated to a certain extent by varying the extension times in PCR.

[Aptamers: selection and scope of applications]

Aptamers are short oligonucleotides selected from large combinatorial pools of sequences for their capacity to bind to many different targets ranging from small molecules (amino acids, antibiotics...) to proteins or nucleic acid structures. Aptamers present the same high specificity and affinity for their targets as antibodies. In addition to efficient binding, aptamers have been shown in many cases to display an inhibitory activity against their targets. Many aptamers are now being developed against biomedical relevant targets, and one aptamer that inhibits the human VEGF165 already received approval for the treatment of age-related macular degeneration. Here we discuss the principles and the practical way of selecting aptamers (SELEX technology) as well as the structural basis for their performance as ligands. A wide scope of applications is described - aptamers have been used as tools for studying nucleic acids/proteins interactions, detecting, purifying or imaging target molecules, regulating gene expression - and includes recent developments of aptamers for therapy and diagnosis.

Diverse Evolutionary Trajectories Characterize a Community of RNA-Cleaving Deoxyribozymes: A Case Study into the Population Dynamics of In Vitro Selection

Two parallel in vitro selections (denoted Selection A and Selection B) were conducted under different selection-pressure regimes, yielding a diverse community of RNA-cleaving deoxyribozymes. In Selection A, the reaction time was reduced four times (from 5 h to 5 s) over the course of 24 generations, while in Selection B the reaction time was maintained at 5 h for 30 rounds of selective amplification. Sequence alignment was conducted on more than 800 clones assembled from 18 generations that span both selections. Many prominent catalytic sequence classes, including some that extend across both selections, were identified and used to construct fitness landscapes depicting their rise and fall over time. The landscapes from both selections exhibit similar global trends despite differences in population dynamics. Some deoxyribozymes were predominant in the early rounds of selection but gave way to other species that dominated in the middle rounds. Ultimately, these middle classes disappeared from the landscape in favor of new and presumably more fit deoxyribozyme sequence classes. The shape of these landscapes alludes to the presence of many latent deoxyribozymes in the initial library, which can only be accessed by changes in the selection pressure and/or by adaptive mutations. Basic computer simulations provide theoretical corroboration of the experimentally observed pattern of staggered sequence-class transitions across the fitness landscapes. These simulations model the influence of one or more contributing factors, including catalytic rate, folding efficiency, PCR amplification efficiency, and random mutagenesis. This is the first study which thoroughly documents the topography of a deoxyribozyme fitness landscape over many generations of in vitro selection.

Expanded divalent metal-ion tolerance of evolved ligase ribozymes

Class I ligases are artificial ribozymes that catalyze the joining of two single-stranded RNAs. These ribozymes are between 120 and 160 nucleotides in length, making them intermediate in size for catalytic RNAs. Previous characterization of the b1-207 ribozyme suggests that it behaves similar to larger ribozymes in terms of divalent metal-ion dependence. This molecule displays a strong preference for magnesium for catalysis, and is inactive in any other metal except manganese, which actually inhibits its operation in magnesium. Here, we sought to examine the metal-ion usages of two ligases that were obtained through continuous evolution in vitro from the b1-207 sequence framework. We found an expanded catalytic range for the E(100)(#3) and B(16)(#19) ribozymes, as they are both catalytically active in calcium and strontium, and less inhibited by manganese. Though not selected for activity in these salts, the evolved ribozymes exhibit several adaptations to in vitro catalysis, and their ability to accommodate metals other than magnesium can be viewed as an example of a molecular exaptation.

Eutectic phase polymerization of activated ribonucleotide mixtures yields quasi-equimolar incorporation of purine and pyrimidine nucleobases

The RNA world hypothesis requires a plausible mechanism by which RNA itself (or precursor RNA-like polymers) can be synthesized nonenzymatically from the corresponding building blocks. Simulation experiments have exploited chemically reactive mononucleotides as monomers. Solutions of such monomers in the prebiotic environment were likely to be very dilute, but in experimental simulations of polymerization reactions dilute solutions of activated mononucleotides in the millimolar range hydrolyze extensively, and only trace amounts of dimers and trimers are formed. We report here that random medium-size RNA analogues with mixed sequences (5- to 17-mers with traces of longer products) can be synthesized in ice eutectic phases that are produced when dilute solutions of activated monomers and catalysts (Mg(II) and Pb(II)) are frozen and maintained at -18 degrees C for periods up to 38 days. Under these conditions, the monomers are concentrated as eutectics in an ice matrix. Hydrolysis of the activated mononucleotides was suppressed at low-temperature ranges, and polymerization was enhanced with yields up to 90%. Analysis of the mixed oligomers established that incorporation of both purine and pyrimidine bases proceeded at comparable rates and yields. These results suggest that ice deposits on the early Earth could have facilitated the synthesis of short- and medium-size random sequence RNA analogues and thereby provided a microenvironment suitable for the formation of biopolymers or their precursors.

Assessing the likelihood of recurrence during RNA evolution in vitro

Recurrence is the possibility of resulting in the same endpoint multiple times when a living system is allowed to evolve repeatedly starting from a given initial point. This concept is of concern to both evolutionary theoreticians and molecular biologists who use nucleic acid selection techniques to mimic biotic and computorial processes in the test tube. Using the continuous in vitro evolution methodology, many replicate experimental evolutionary lineages with populations of catalytic RNA were performed to gain insight into the parameters that could affect recurrence. The likelihood that the same genotype will result in parallel trials of an evolution experiment in vitro depends on several factors, including the phenotype under selection, the size and composition of the initial diverse pool of nucleic acids used in the experiment, the degree of mutation possible during the experiment, the shape of the fitness landscape through which the population evolves, and the strategies used to invoke selection and to search the landscape, among others. By considering these factors, it can be predicted that recurrence is more likely when a small, wild-type-based starting pool is used with efficient selection and search strategies involving little online mutagenesis within a rugged adaptive landscape with a strong local optimum. The recurrence experiments performed here on the 150-nucleotide ligase ribozyme demonstrate that it repeatedly jumps from one peak in a fitness landscape to another, apparently hurdling a deep fitness valley. These predictions can and should be tested by additional multiple replicates of actual evolution experiments in the laboratory.