Analytical description of finite size effects for RNA secondary structures

The ensemble of RNA secondary structures of uniform sequences is studied analytically. We calculate the partition function for very long sequences and discuss how the cross-over length, beyond which asymptotic scaling laws apply, depends on thermodynamic parameters. For realistic choices of parameters this length can be much longer than natural RNA molecules. This has to be taken into account when applying asymptotic theory to interpret experiments or numerical results.Comment: 10 pages, 13 figures, published in Phys. Rev.

Qualification Tests of 474 Photomultiplier Tubes for the Inner Detector of the Double Chooz Experiment

The hemispherical 10" photomultiplier tube (PMT) R7081 from Hamamatsu Photonics K.K. (HPK) is used in various experiments in particle and astroparticle physics. We describe the test and calibration of 474 PMTs for the reactor antineutrino experiment Double Chooz. The unique test setup at Max-Planck-Institut f\"ur Kernphysik Heidelberg (MPIK) allows one to calibrate 30 PMTs simultaneously and to characterize the single photo electron response, transit time spread, linear behaviour and saturation effects, photon detection efficiency and high voltage calibration

Qualification Tests of 474 Photomultiplier Tubes for the Inner Detector of the Double Chooz Experiment

The hemispherical 10" photomultiplier tube (PMT) R7081 from Hamamatsu Photonics K.K. (HPK) is used in various experiments in particle and astroparticle physics. We describe the test and calibration of 474 PMTs for the reactor antineutrino experiment Double Chooz. The unique test setup at Max-Planck-Institut f\"ur Kernphysik Heidelberg (MPIK) allows one to calibrate 30 PMTs simultaneously and to characterize the single photo electron response, transit time spread, linear behaviour and saturation effects, photon detection efficiency and high voltage calibration

Qualification Tests of 474 Photomultiplier Tubes for the Inner Detector of the Double Chooz Experiment

The hemispherical 10" photomultiplier tube (PMT) R7081 from Hamamatsu Photonics K.K. (HPK) is used in various experiments in particle and astroparticle physics. We describe the test and calibration of 474 PMTs for the reactor antineutrino experiment Double Chooz. The unique test setup at Max-Planck-Institut f\"ur Kernphysik Heidelberg (MPIK) allows one to calibrate 30 PMTs simultaneously and to characterize the single photo electron response, transit time spread, linear behaviour and saturation effects, photon detection efficiency and high voltage calibration

RNA denaturation: excluded volume, pseudoknots and transition scenarios

A lattice model of RNA denaturation which fully accounts for the excluded volume effects among nucleotides is proposed. A numerical study shows that interactions forming pseudoknots must be included in order to get a sharp continuous transition. Otherwise a smooth crossover occurs from the swollen linear polymer behavior to highly ramified, almost compact conformations with secondary structures. In the latter scenario, which is appropriate when these structures are much more stable than pseudoknot links, probability distributions for the lengths of both loops and main branches obey scaling with nonclassical exponents.Comment: 4 pages 3 figure

Translocation of structured polynucleotides through nanopores

We investigate theoretically the translocation of structured RNA/DNA molecules through narrow pores which allow single but not double strands to pass. The unzipping of basepaired regions within the molecules presents significant kinetic barriers for the translocation process. We show that this circumstance may be exploited to determine the full basepairing pattern of polynucleotides, including RNA pseudoknots. The crucial requirement is that the translocation dynamics (i.e., the length of the translocated molecular segment) needs to be recorded as a function of time with a spatial resolution of a few nucleotides. This could be achieved, for instance, by applying a mechanical driving force for translocation and recording force-extension curves (FEC's) with a device such as an atomic force microscope or optical tweezers. Our analysis suggests that with this added spatial resolution, nanopores could be transformed into a powerful experimental tool to study the folding of nucleic acids.Comment: 9 pages, 5 figure

Statistical mechanics of RNA folding: importance of alphabet size

We construct a minimalist model of RNA secondary-structure formation and use it to study the mapping from sequence to structure. There are strong, qualitative differences between two-letter and four or six-letter alphabets. With only two kinds of bases, there are many alternate folding configurations, yielding thermodynamically stable ground-states only for a small set of structures of high designability, i.e., total number of associated sequences. In contrast, sequences made from four bases, as found in nature, or six bases have far fewer competing folding configurations, resulting in a much greater average stability of the ground state.Comment: 7 figures; uses revtex

Zero Temperature Properties of RNA Secondary Structures

We analyze different microscopic RNA models at zero temperature. We discuss both the most simple model, that suffers a large degeneracy of the ground state, and models in which the degeneracy has been remove, in a more or less severe manner. We calculate low-energy density of states using a coupling perturbing method, where the ground state of a modified Hamiltonian, that repels the original ground state, is determined. We evaluate scaling exponents starting from measurements of overlaps and energy differences. In the case of models without accidental degeneracy of the ground state we are able to clearly establish the existence of a glassy phase with $\theta \simeq 1/3$.Comment: 20 pages including 9 eps figure

Statistical mechanics of secondary structures formed by random RNA sequences

The formation of secondary structures by a random RNA sequence is studied as a model system for the sequence-structure problem omnipresent in biopolymers. Several toy energy models are introduced to allow detailed analytical and numerical studies. First, a two-replica calculation is performed. By mapping the two-replica problem to the denaturation of a single homogeneous RNA in 6-dimensional embedding space, we show that sequence disorder is perturbatively irrelevant, i.e., an RNA molecule with weak sequence disorder is in a molten phase where many secondary structures with comparable total energy coexist. A numerical study of various models at high temperature reproduces behaviors characteristic of the molten phase. On the other hand, a scaling argument based on the extremal statistics of rare regions can be constructed to show that the low temperature phase is unstable to sequence disorder. We performed a detailed numerical study of the low temperature phase using the droplet theory as a guide, and characterized the statistics of large-scale, low-energy excitations of the secondary structures from the ground state structure. We find the excitation energy to grow very slowly (i.e., logarithmically) with the length scale of the excitation, suggesting the existence of a marginal glass phase. The transition between the low temperature glass phase and the high temperature molten phase is also characterized numerically. It is revealed by a change in the coefficient of the logarithmic excitation energy, from being disorder dominated to entropy dominated.Comment: 24 pages, 16 figure

RNA secondary structure prediction from multi-aligned sequences

It has been well accepted that the RNA secondary structures of most functional non-coding RNAs (ncRNAs) are closely related to their functions and are conserved during evolution. Hence, prediction of conserved secondary structures from evolutionarily related sequences is one important task in RNA bioinformatics; the methods are useful not only to further functional analyses of ncRNAs but also to improve the accuracy of secondary structure predictions and to find novel functional RNAs from the genome. In this review, I focus on common secondary structure prediction from a given aligned RNA sequence, in which one secondary structure whose length is equal to that of the input alignment is predicted. I systematically review and classify existing tools and algorithms for the problem, by utilizing the information employed in the tools and by adopting a unified viewpoint based on maximum expected gain (MEG) estimators. I believe that this classification will allow a deeper understanding of each tool and provide users with useful information for selecting tools for common secondary structure predictions.Comment: A preprint of an invited review manuscript that will be published in a chapter of the book `Methods in Molecular Biology'. Note that this version of the manuscript may differ from the published versio