A simple method for computing exact probabilities of mutation numbers

We describe a method for the recursive computation of exact probability distributions for the number of neutral mutations segregating in samples of arbitrary size and configuration. Construction of the recursions requires only characterization of evolutionary changes as a Markov process and determination of one-step transition matrices. We address the pattern of nucleotide diversity at a neutral marker locus linked to a determinant of mating type. Under a reformulation of parameters, the method also applies directly to metapopulation models with island migration among demes. Characterization of complete probability distributions facilitates parameter estimation and hypothesis testing by likelihood- as well as moment-based approaches.     

Partition function and base pairing probabilities of RNA heterodimers

Background  

RNA has been recognized as a key player in cellular regulation in recent years. In many cases, non-coding RNAs exert their function by binding to other nucleic acids, as in the case of microRNAs and snoRNAs. The specificity of these interactions derives from the stability of inter-molecular base pairing. The accurate computational treatment of RNA-RNA binding therefore lies at the heart of target prediction algorithms.     

Prediction of RNA base pairing probabilities on massively parallel computers.

We present an implementation of McCaskill's algorithm for computing the base pair probabilities of an RNA molecule for massively parallel message passing architectures. The program can be used to routinely fold RNA sequences of more than 10,000 nucleotides. Applications to complete viral genomes are discussed.     

GUUGle: a utility for fast exact matching under RNA complementary rules including G-U base pairing

MOTIVATION: RNA secondary structure analysis often requires searching for potential helices in large sequence data. RESULTS: We present a utility program GUUGle that efficiently locates potential helical regions under RNA base pairing rules, which include Watson-Crick as well as G-U pairs. It accepts a positive and a negative set of sequences, and determines all exact matches under RNA rules between positive and negative sequences that exceed a specified length. The GUUGle algorithm can also be adapted to use a precomputed suffix array of the positive sequence set. We show how this program can be effectively used as a filter preceding a more computationally expensive task such as miRNA target prediction. AVAILABILITY: GUUGle is available via the Bielefeld Bioinformatics Server at http://bibiserv.techfak.uni-bielefeld.de/guugle     

Mutation as an error in base pairing

Summary The model of mutation by transitional change (Freese 1959) predicts that a heritable change in genotype is established when two replications of DNA succeed the initial incorporation of an analogue. The model was tested in populations ofSalmonella typhimurium strainstryD-10 andtryD-79 whose division had been synchronized by fractional filtration. Mutation from auxotrophy to prototrophy (try ^–try ⁺) induced by 5-bromodeoxyuridine (BUDR) and 2-aminopurine (AP) occurred in accordance with DNA replication. Two subsequent DNA replications were necessary to obtain BUDR-induced prototrophs inD-79, one subsequent DNA replication was required for AP-induced prototrophs inD-79, while no subsequent DNA replication was necessary for AP-induced prototrophs inD-10. This was observed whether the mutagens were present continuously or during only the first replication and also when the cells were allowed to replicate their DNA without cell division in the presence of inhibitory concentrations of the base analogue or when protein synthesis was blocked in the presence of chloramphenicol. A statistical analysis of the patterns of mutant increase observed for six mutant strains was used to distinguish between errors in replication and errors in incorporation induced by the base analogues and thereby the base pair at the mutant site was identified.With 10 Figures in the TextSupported in part by grants from the American Cancer Society the U.S. Public Health Service and the National Science Foundation administered by ProfessorF. J. Ryan.     

Secondary structure of the 3'-noncoding region of flavivirus genomes: comparative analysis of base pairing probabilities.

The prediction of the complete matrix of base pairing probabilities was applied to the 3'' noncoding region (NCR) of flavivirus genomes. This approach identifies not only well-defined secondary structure elements, but also regions of high structural flexibility. Flaviviruses, many of which are important human pathogens, have a common genomic organization, but exhibit a significant degree of RNA sequence diversity in the functionally important 3''-NCR. We demonstrate the presence of secondary structures shared by all flaviviruses, as well as structural features that are characteristic for groups of viruses within the genus reflecting the established classification scheme. The significance of most of the predicted structures is corroborated by compensatory mutations. The availability of infectious clones for several flaviviruses will allow the assessment of these structural elements in processes of the viral life cycle, such as replication and assembly.     

Informational parameters of an exact DNA base sequence

A series of hierarchical chemical reactivity calculations have been performed to elucidate the alkylation properties of a methyldiazonium ion toward DNA base sites. Both MINDO/3 and CNDO/2 approximate methods have been employed. For the isolated bases the O6 of guanine is predicted to be the most reactive site. This prediction may also be relevant to single-stranded DNA chains containing guanine. For base-pairs, the N7 and O6 sites on guanine are about equally favored for alkylation. The previous study of aziridinium ion alkylation gave about the same results with N7 guanine modestly favored as the preferred site of alkylation for base-pairs. In composite we conclude that N7 guanine and/or O6 guanine will be the preferred sites for alkylation by a methyldiazonium ion but cannot distinguish between these two in terms of chemical specificity.     

Model systems for understanding DNA base pairing

The fact that nucleic acid bases recognize each other to form pairs is a canonical part of the dogma of biology. However, they do not recognize each other well enough in water to account for the selectivity and efficiency that is needed in the transmission of biological information through a cell. Thus proteins assist in this recognition in multiple ways, and recent data suggest that these mechanisms of recognition can vary widely with context. To probe how the chemical differences of the four nucleobases are defined in various biological contexts, chemists and biochemists have developed modified versions that differ in their polarity, shape, size, and functional groups. This brief review covers recent advances in this field of research.     

Comparison of an exact and a simulation method for calculating gene extinction probabilities in pedigrees

A comparison is made between a much used simulation method, commonly called gene dropping, and the exact computational technique of peeling. These methods are illustrated using the problem of finding the distribution of the number of distinct ancestral genes surviving at an autosomal locus. Each method is used on several real zoo pedigrees, of varying size and complexity, and the results are compared. Gene dropping is found to be a good approximation to peeling, but for all but the most complex pedigrees surveyed, peeling is preferable. The relationship between heterozygosity and allelic variability is investigated.     

Oxetane locked thymidine in the Dickerson-Drew dodecamer causes local base pairing distortions -- an NMR structure and hydration study

The introduction of a North-type sugar conformation constrained oxetane T block, 1-(1',3'-O-anhydro-beta-D-psicofuranosyl) thymine, at the T(7) position of the self-complementary Dickerson-Drew dodecamer, d[(5'-C(1)G(2)C(3)G(4)A(5)A(6)T(7)T(8)C(9)G(10)C(11)G(12)-3')](2), considerably perturbs the conformation of the four central base pairs, reducing the stability of the structure. UV spectroscopy and 1D NMR display a drop in melting temperature of approximately 10 degrees C per modification for the T(7) oxetane modified duplex, where the T(7) block has been introduced in both strands, compared to the native Dickerson-Drew dodecamer. The three dimensional structure has been determined by NMR spectroscopy and has subsequently been compared with the results of 2.4 ns MD simulations of the native and the T(7) oxetane modified duplexes. The modified T(7) residue is found to maintain its constrained sugar- and the related glycosyl torsion conformations in the duplex, resulting in staggered and stretched T(7).A(6) and A(6).T(7) non-linear base pairs. The stacking is less perturbed, but there is an increased roll between the two central residues compared to the native counterpart, which is compensated by tilts of the neighboring base steps. The one dimensional melting profile of base protons of the T(7) and T(8) residues reveals that the introduction of the North-type sugar constrained thymine destabilizes the core of the modified duplex, promoting melting to start simultaneously from the center as well as from the ends. Temperature dependent hydration studies by NMR demonstrate that the central T(7).A(6)/A(6).T(7) base pairs of the T(7) oxetane modified Dickerson-Drew dodecamer have at least one order of magnitude higher water exchange rates (correlated to the opening rate of the base pair) than the corresponding base pairs in the native duplex.     

RNA structure and dynamics: A base pairing perspective

RNA is now known to possess various structural, regulatory and enzymatic functions for survival of cellular organisms. Functional RNA structures are generally created by three-dimensional organization of small structural motifs, formed by base pairing between self-complementary sequences from different parts of the RNA chain. In addition to the canonical Watson–Crick or wobble base pairs, several non-canonical base pairs are found to be crucial to the structural organization of RNA molecules. They appear within different structural motifs and are found to stabilize the molecule through long-range intra-molecular interactions between basic structural motifs like double helices and loops. These base pairs also impart functional variation to the minor groove of A-form RNA helices, thus forming anchoring site for metabolites and ligands. Non-canonical base pairs are formed by edge-to-edge hydrogen bonding interactions between the bases. A large number of theoretical studies have been done to detect and analyze these non-canonical base pairs within crystal or NMR derived structures of different functional RNA. Theoretical studies of these isolated base pairs using ab initio quantum chemical methods as well as molecular dynamics simulations of larger fragments have also established that many of these non-canonical base pairs are as stable as the canonical Watson–Crick base pairs. This review focuses on the various structural aspects of non-canonical base pairs in the organization of RNA molecules and the possible applications of these base pairs in predicting RNA structures with more accuracy.     

A rRNA-mRNA base pairing model for UGA-dependent termination.

A series of site-directed mutations has been constructed in E coli 16S rRNA and shown to suppress UGA-dependent translational termination. With the exception of the C726 to G base change, all were constructed in helix 34. Characterization of these mutations is reviewed here and from these data and mRNA-rRNA base pairing model for the termination event is presented. The interaction functions via antiparallel base pairing between either 1 of the 2 UCA motifs in helix 34 and the complementary UGA stop codon on the message, thus forming a quasicontinuous A-type helical structure that is further stabilized by stacking enthalpy. Finally, rRNA motifs potentially required for UAA and UAG-dependent translational termination are discussed.     

A fast algorithm for computing multilocus recombination

A fast algorithm for computing recombination is developed for model organisms with selection on haploids. Haplotype frequencies are transformed to marginal frequencies; random mating and recombination are computed; marginal frequencies are transformed back to haplotype frequencies. With L diallelic loci, this algorithm is theoretically a factor of a constant times (3/8)^L faster than standard computations with selection on diploids, and up to 16 recombining loci have been computed. This algorithm is then applied to model the opposing evolutionary forces of multilocus epistatic selection and recombination. Selection is assumed to favor haplotypes with specific alleles either all present or all absent. When the number of linked loci exceeds a critical value, a jump bifurcation occurs in the two-dimensional parameter space of the selection coefficient s and the recombination frequency r. The equilibrium solution jumps from high to low mean fitness with increasing r or decreasing s. These numerical results display an unexpected and dramatic nonlinear effect occurring in linkage models with a large number of loci.     

'Extremotaxis': computing with a bacterial-inspired algorithm

We present a general-purpose optimization algorithm inspired by "run-and-tumble", the biased random walk chemotactic swimming strategy used by the bacterium Escherichia coli to locate regions of high nutrient concentration The method uses particles (corresponding to bacteria) that swim through the variable space (corresponding to the attractant concentration profile). By constantly performing temporal comparisons, the particles drift towards the minimum or maximum of the function of interest. We illustrate the use of our method with four examples. We also present a discrete version of the algorithm. The new algorithm is expected to be useful in combinatorial optimization problems involving many variables, where the functional landscape is apparently stochastic and has local minima, but preserves some derivative structure at intermediate scales.     

BALSA: Bayesian algorithm for local sequence alignment

The Smith–Waterman algorithm yields a single alignment, which, albeit optimal, can be strongly affected by the choice of the scoring matrix and the gap penalties. Additionally, the scores obtained are dependent upon the lengths of the aligned sequences, requiring a post-analysis conversion. To overcome some of these shortcomings, we developed a Bayesian algorithm for local sequence alignment (BALSA), that takes into account the uncertainty associated with all unknown variables by incorporating in its forward sums a series of scoring matrices, gap parameters and all possible alignments. The algorithm can return both the joint and the marginal optimal alignments, samples of alignments drawn from the posterior distribution and the posterior probabilities of gap penalties and scoring matrices. Furthermore, it automatically adjusts for variations in sequence lengths. BALSA was compared with SSEARCH, to date the best performing dynamic programming algorithm in the detection of structural neighbors. Using the SCOP databases PDB40D-B and PDB90D-B, BALSA detected 19.8 and 41.3% of remote homologs whereas SSEARCH detected 18.4 and 38% at an error rate of 1% errors per query over the databases, respectively.     

Nuclear Overhauser effect study of yeast tRNAVal 1: evidence for uridine-pseudouridine base pairing.

The proton NMR spectrum of yeast tRNAVal 1 has been studied using nuclear Overhauser effect (NOE), including comparison of NOE patterns between purine C8 deuterated and nondeuterated samples. Studies of the downfield region enable us to reliably assign many resonances in the acceptor and D stems. Prominent among these reliable assignments is that of the unusual base pair U psi, which is made here for the first time. Other identifications include GU2, U8-A14, the three AU base pairs of the acceptor stem, and N1 and N3 protons of psi 55.     

Hoogsteen G-G base pairing is dispensable for telomere healing in yeast.

The G-rich strands of most eukaryotic telomeres are capable of forming highly folded structures in vitro, mediated, in part, through Hoogsteen G-G base pairing. The ability of most telomeres to form these structures has led to the suggestion that they play an important role in telomere addition. I have investigated this possibility in the yeast Saccharomyces cerevisiae through the use of an in vivo assay that measures healing via poly(G1-3T) addition onto plasmid substrates containing synthetic telomeres. Synthetic telomere healing is a highly size- and sequence-specific process that allows the discrimination of telomeres of differing efficiency. Plasmids containing synthetic telomeres with differing abilities to form secondary structures were tested in this assay for healing in vivo. The results of this study demonstrate that telomeres incapable of forming Hoogsteen base pairs nonetheless serve as efficient substrates for poly(G1-3T) addition, indicating that intramolecular Hoogsteen G-G base pairing is not essential for this process.     

Parallel nucleic acid helices with hoogsteen base pairing: Symmetry and structure

Molecular structures for parallel DNA and RNA double helices with Hoogsteen pairing are proposed for the first time. The DNA helices have sugars in the C2′-endo region and the phosphodiester conformations are (trans, gauche^?), and the RNA helices have sugars in the C3′-endo region and the phosphodiester conformations are (gauche^?, gauche^?). A pseudorotational symmetry relates the two parallel strands of DNA helices and a screw symmetry relates the two strands of RNA helices, which have an associated tilt of the The conformational space of parallel helices with Hoogsteen base pairing, unlike the Watson-Crick duplex, is highly restricted due to the unique positioning of the symmetry axis in the former case. The features of the parallel double helix with Hoogsteen pairing are compared with the Watson-Crick duplex and the corresponding triple helix. © 1994 John Wiley & Sons, Inc.