Algorithm independent properties of RNA secondary structure predictions

Algorithms predicting RNA secondary structures based on different folding criteria – minimum free energies (mfe), kinetic folding (kin), maximum matching (mm) – and different parameter sets are studied systematically. Two base pairing alphabets were used: the binary GC and the natural four-letter AUGC alphabet. Computed structures and free energies depend strongly on both the algorithm and the parameter set. Statistical properties, such as mean number of base pairs, mean numbers of stacks, mean loop sizes, etc., are much less sensitive to the choice of parameter set and even of algorithm. Some features of RNA secondary structures, such as structure correlation functions, shape space covering and neutral networks, seem to depend only on the base pairing logic (GC or AUGC alphabet). Received: 16 May 1996 / Accepted: 10 July 1996     

A novel family of linear plasmids with homology to plasmid pAL2-1 of Podospora anserina

Three recently isolated wild-type strains of the ascomycete Podospora anserina were analyzed for the presence of linear mitochondrial plasmids. In one of these strains, designated Wa6, at least 12 distinct plasmid-like elements were identified. From molecular analyses a minimum number of 78 individual linear molecules with proteins bound to their 5′ ends was estimated. In addition, the different members of this family of typical linear plasmids were shown to possess a common central region and terminal sequences which differ from one plasmid to another due to the presence of different numbers of a 2.4 kb sequence module. Finally, the pWa6 plasmids share a high degree of sequence similarity with pAL2-1, a linear plasmid previously identified in mitochondria of a long-lived mutant of P.anserina. A mechanism is proposed which explains the generation of these distinct, closely related extrachromosomal genetic traits.     

Statistics of RNA melting kinetics

We present and study the behavior of a simple kinetic model for the melting of RNA secondary structures, given that those structures are known. The model is then used as a map that. assigns structure dependent overall rate constants of melting (or refolding) to a sequence. This induces a landscape of reaction rates, or activation energies, over the space of sequences with fixed length. We study the distribution and the correlation structure of these activation energies. Correspondence to: P. Schuster     

RNA based evolutionary optimization

The notion of an RNA world has been introduced for a prebiotic scenario that is dominated by RNA molecules and their properties, in particular their capabilities to act as templates for reproduction and as catalysts for several cleavage and ligation reactions of polynucleotides and polypeptides. This notion is used here also for simple experimental assays which are well suited to study evolution in the test tube. In molecular evolution experiments fitness is determined in essence by the molecular structures of RNA molecules. Evidence is presented for adaptation to environment in cell-free media. RNA based molecular evolution experiments have led to interesting spin-offs in biotechnology, commonly called applied molecular evolution, which make use of Darwinian trial-and-error strategies in order to synthesize new pharmacological compounds and other advanced materials on a biological basis.Error-propagation in RNA replication leads to formation of mutant spectra called quasispecies. An increase in the error rate broadens the mutant spectrum. There exists a sharply defined threshold beyond which heredity breaks down and evolutionary adaptation becomes impossible. Almost all RNA viruses studied so far operate at conditions close to this error threshold. Quasispecies and error thresholds are important for an understanding of RNA virus evolution, and they may help to develop novel antiviral strategies.Evolution of RNA molecules can be studied and interpreted by considering secondary structures. The notion of sequence space introduces a distance between pairs of RNA sequences which is tantamount to counting the minimal number of point mutations required to convert the sequences into each other. The mean sensitivity of RNA secondary structures to mutation depends strongly on the base pairing alphabet: structures from sequences which contain only one base pair (GC or AU are much less stable against mutation than those derived from the natural (AUGC) sequences. Evolutionary optimization of two-letter sequences in thus more difficult than optimization in the world of natural RNA sequences with four bases. This fact might explain the usage of four bases in the genetic language of nature.Finally we study the mapping from RNA sequences into secondary structures and explore the topology of RNA shape space. We find that neutral paths connecting neighbouring sequences with identical structures go very frequently through entire sequence space. Sequences folding into common structures are found everywhere in sequence space. Hence, evolution can migrate to almost every part of sequence space without hill climbing and only small fractions of the entire number of sequences have to be searched in order to find suitable structures.     

Polysulfide as a possible substrate for sulfur-reducing bacteria

Because of its low solubility it is unlikely that elemental sulfur serves as the direct substrate for sulfur-reducing bacteria. To test the hypothesis that polysulfide may represent a soluble intermediate of sulfur reduction, the maximal polysulfide concentrations formed from elemental sulfur in aqueous sulfide solutions were measured at near neutral pH and at temperatures up to 90°C. The saturation concentrations decreased by two orders of magnitude when the pH was lowered from 7 to 6 at a given temperature, and increased about tenfold when the temperature was raised from 37°C to 90°C at a given pH. The dissolution of 0.1 mM zerovalent sulfur in 1 mM sulfide (H₂S+HS^–) required a pH of 7.5 at 20°C and of only 6.1 at 100°C. A comparison with the growth optima of sulfur-reducers suggests that polysulfide is present at sufficient concentration at the growth conditions of the Bacteria and the moderately acidophilic Archaea. Polysulfide is apparently not available at the growth conditions of the extremely acidophilic Archaea. Alternative mechanisms for the sulfur utilization under these conditions are discussed.Abbreviations MOPS Morpholinopropanesulfonate - PIPES 1,4 piperazine-N,N-bis(2-ethanesulfonate) - HEPES N-2-hydroxy-ethylpiperazine-N-ethanesulfonate     

Refined algorithm and computer program for calculating all non-negative fluxes admissible in steady states of biochemical reaction systems with or without some flux rates fixed

A refined algorithm together with a computer procedure for determiningthe complete set of non–negative, steady–state fluxesin biochemical reaction systems of any complexily, with or withoutsome flux rates fixed, is given. It is shown that this set isa convex polyhedron, which may or may not be bounded. The algorithmis illustrated by several examples; one of them concerns intermediarymetabolism. A computer cade in standard C is presented.     

Reaction intensification for biocatalytic production of polyphenolic natural product di-C-β-glucosides

Polyphenolic aglycones featuring two sugars individually attached via C-glycosidic linkage (di-C-glycosides) represent a rare class of plant natural products with unique physicochemical properties and biological activities. Natural scarcity of such di-C-glycosides limits their use-inspired exploration as pharmaceutical ingredients. Here, we show a biocatalytic process technology for reaction-intensified production of the di-C-β-glucosides of two representative phenol substrates, phloretin (a natural flavonoid) and phenyl-trihydroxyacetophenone (a phenolic synthon for synthesis), from sucrose. The synthesis proceeds via an iterative two-fold C-glycosylation of the respective aglycone, supplied as inclusion complex with 2-hydroxypropyl β-cyclodextrin for enhanced water solubility of up to 50 mmol/L, catalyzed by a kumquat di-C-glycosyltransferase (di-CGT), and it uses UDP-Glc provided in situ from sucrose by a soybean sucrose synthase, with catalytic amounts (≤3 mol%) of UDP added. Time course analysis reveals the second C-glycosylation as rate-limiting (0.4–0.5 mmol/L/min) for the di-C-glucoside production. With internal supply from sucrose keeping the UDP-Glc at a constant steady-state concentration (≥50% of the UDP added) during the reaction, the di-C-glycosylation is driven to completion (≥95% yield). Contrary to the mono-C-glucoside intermediate which is stable, the di-C-glucoside requires the addition of reducing agent (10 mmol/L 2-mercaptoethanol) to prevent its decomposition during the synthesis. Both di-C-glucosides are isolated from the reaction mixtures in excellent purity (≥95%), and their expected structures are confirmed by NMR. Collectively, this study demonstrates efficient glycosyltransferase cascade reaction for flexible use in natural product di-C-β-glucoside synthesis from expedient substrates.     

The influence of prey consumption and demographic stochasticity on population growth rate of Isle Royale wolves Canis lupus

The relationship between the rates of prey capture and predator population growth is a fundamental aspect of predation, yet it is rarely measured for vertebrate predators. For the isolated wolf population on Isle Royale, annual variation in kill rate explains 22% of the variation in wolf population growth rate. From the slope of this relationship, we estimate that the production efficiency (ratio of production to respiration) of wolves is between 0.5% and 1.5%. More generally, we assess the relative extent to which wolf population growth rate is affected by density dependence, prey availability (moose, Alces alces ), winter weather, and demographic stochasticity. Prey availability explains the most variation in wolf growth rate (42%), but this is only recognized after accounting for the influence of a disease-induced population crash and age structure of the prey population (i.e. number of vulnerable moose, >9 years of age). Demographic stochasticity accounts for approximately 30% of the variation in wolf growth rate. This recognition is important, but not surprising, given that the average population size of Isle Royale wolves is 22. Previous work indicates that the effect of winter climate, as mediated through prey vulnerability and kill rates, is substantial. This work indicates that the direct effect of winter climate is weak, and explains only about 4% of the variation in wolf growth rate (P=0.10).