ban operon of bacteriophage P1. Mutational analysis of the c1 repressor-controlled operator

The repressor of bacteriophage P1, encoded by the c1 gene, represses the phage lytic functions and is responsible for maintaining the P1 prophage in the lysogenic state. The c1 repressor interacts with at least 11 binding sites or operators widely scattered over the P1 genome. From these operators, a 17 base-pair asymmetric consensus sequence, ATTGCTCTAATAAATTT, was derived. Here, we show that the operator, Op72 of the P1ban operon consists of two overlapping 17 base-pair sequences a and b forming an incomplete palindrome. Op72a matches the consensus sequence, whereas Op72b contains two mismatches. The evidence is based on the sequence analysis of 27 operator mutants constitutive for ban expression. They were identified as single-base substitutions at positions 2 to 10 of Op72a (26 mutants) and at position 8 of Op72b (one mutant). We conclude from gel retardation and footprinting studies that two repressor molecules bind to the operator and that positions 4, 5 and 7 to 10 of the operator play an essential role in repressor recognition.     

Dynamics of small autocatalytic reaction networks—I. bifurcations,permanence and exclusion

Catalysis in replication networks has become an important issue in biophysics and other areas of biology. Examples are RNA catalysis, idiotype recognition in the immune response and dynamical models of Maynard-Smith games in sociobiology. Chemical reaction networks describing catalysed, template-induced reproduction of three species are analysed in full generality. The nine-dimensional parameter space is reduced to three relevant angular coordinates which determine completely the phase portraits (PPs) and the bifurcation patterns. All cases are classified and all generic as well as most of the nongeneric transitions are listed and described. This paper has been reproduced directly from disc using a LA-T_EX system.     

Alveolar slope and dead space of He and SF6 in dogs: comparison of airway and venous loading

Series (Fowler) dead space (VD) and slope of the alveolar plateau of two inert gases (He and SF6) with similar blood-gas partition coefficients (approximately 0.01) but different diffusivities were analyzed in 10 anesthetized paralyzed mechanically ventilated dogs (mean body wt 20 kg). Single-breath constant-flow expirograms were simultaneously recorded in two conditions: 1) after equilibration of lung gas with the inert gases at tracer concentrations [airway loading (AL)] and 2) during steady-state elimination of the inert gases continuously introduced into venous blood by a membrane oxygenator and partial arteriovenous bypass [venous loading (VL)]. VD was consistently larger for SF6 than for He, but there was no difference between AL and VL. The relative alveolar slope, defined as increment of partial pressure per increment of expired volume and normalized to mixed expired-inspired partial pressure difference, was larger by a factor of two in VL than in AL for both He and SF6. The He-to-SF6 ratio of relative alveolar slope was generally smaller than unity in both VL and AL. Whereas unequal ventilation-volume distribution combined with sequential emptying of parallel lung regions appears to be responsible for the sloping alveolar plateau during AL, the steeper slope during VL is attributed to the combined effects of continuing gas exchange and ventilation-perfusion inequality coupled with sequential emptying. The differences between He and SF6 point at the contributing role of diffusion-dependent mechanisms in intrapulmonary gas mixing.     

Transfer of rps19 to the nucleus involves the gain of an RNP-binding motif which may functionally replace RPS13 in Arabidopsis mitochondria.

The discovery of disrupted rps19 genes in Arabidopsis mitochondria prompted speculation about the transfer to the nuclear compartment. We here describe the functional gene transfer of rps19 into the nucleus of Arabidopsis. Molecular cloning and sequence analysis of rps19 show that the nuclear gene encodes a long N-terminal extension. Import studies of the precursor protein indicate that only a small part of this extension is cleaved off during import. The larger part of the extension, which shows high similarity to conserved RNA-binding domains of the RNP-CS type, became part of the S19 protein. In the Escherichia coli ribosome S19 forms an RNA-binding complex as heterodimer with S13. By using immuno-analysis and import studies we show that a eubacterial-like S13 protein is absent from Arabidopsis mitochondria, and is not substituted by either a chloroplastic or a cytosolic homologue of this ribosomal protein. We therefore propose that either a highly diverged or missing RPS13 has been functionally replaced by an RNP domain that most likely derived from a glycine-rich RNA-binding protein. These results represent the first case of a functional replacement of a ribosomal protein by a common RNA-binding domain and offer a new view on the flexibility of biological systems in using well-adapted functional domains for different jobs.     

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     

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.     

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.