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     

Approaching the multifaceted nature of energy metabolism: inactivation of the cytosolic creatine kinases via homologous recombination in mouse embryonic stem cells

To study the physiological role of the creatine kinase/phosphocreatine (CK/PCr) system in cells and tissues with a high and fluctuating energy demand we have concentrated on the site-directed inactivation of the B- and M-CK genes encoding the cytosolic CK protein subunits. In our approach we used homologous recombination in mouse embryonic stem (ES) cells from strain 129/Sv. Using targeting constructs based on strain 129/Sv isogenic DNA we managed to ablate the essential exons of the B-CK and M-CK genes at reasonably high frequencies. ES clones with fully disrupted B-CK and two types of M-CK gene mutations, a null (M-CK^–) and leaky (M-CK¹) mutation, were used to generate chimaeric mutant mice via injection in strain C57BL/6 derived blastocysts. Chimaeras with the B-CK null mutation have no overt abnormalities but failed to transmit the mutation to their offspring. For the M-CK^– and M-CK¹ mutations successful transmission was achieved and heterozygous and homozygous mutant mice were bred. Animals deficient in MM-CK are phenotypically normal but lack muscular burst activity. Fluxes through the CK reaction in skeletal muscle are highly impaired and fast fibres show adaptation in cellular architecture and storage of glycogen. Mice homozygous for the leaky M-CK allele, which have 3-fold reduced MM-CK activity, show normal fast fibres but CK fluxes and burst activity are still not restored to wildtype levels.     

CENTRAL MOMENTS AND PROBABILITY DISTRIBUTION OF COLLESS'S COEFFICIENT OF TREE IMBALANCE

The great increase in the number of phylogenetic studies of a wide variety of organisms in recent decades has focused considerable attention on the balance of phylogenetic trees—the degree to which sister clades within a tree tend to be of equal size—for at least two reasons: (1) the degree of balance of a tree may affect the accuracy of estimates of it; (2) the degree of balance, or imbalance, of a tree may reveal something about the macroevolutionary processes that produced it. In particular, variation among lineages in rates of speciation or extinction is expected to produce trees that are less balanced than those that result from phylogenetic evolution in which each extant species of a group has the same probability of speciation or extinction. Several coefficients for measuring the balance or imbalance of phylogenetic trees have been proposed. I focused on Colless's coefficient of imbalance (7) for its mathematical tractability and ease of interpretation. Earlier work on this statistic produced exact methods only for calculating the expected value. In those studies, the variance and confidence limits, which are necessary for testing the departure of observed values of I from the expected, were estimated by Monte Carlo simulation. I developed recursion equations that allow exact calculation of the mean, variance, skewness, and complete probability distribution of I for two different probability-generating models for bifurcating tree shapes. The Equal-Rates Markov (ERM) model assumes that trees grow by the random speciation and extinction of extant species, with all species that are extant at a given time having the same probability of speciation or extinction. The Equal Probability (EP) model assumes that all possible labeled trees for a given number of terminal taxa have the same probability of occurring. Examples illustrate how these theoretically derived probabilities and parameters may be used to test whether the evolution of a monophyletic group or set of monophyletic groups has proceeded according to a Markov model with equal rates of speciation and extinction among species, that is, whether there has been significant variation among lineages in expected rates of speciation or extinction.     

Energy metabolism of Chironomus anthracinus (Diptera: Chironomidae) from the profundal zone of Lake Esrom,Denmark, as a function of body size,temperature and oxygen concentration

The profundal zone of Lake Esrom, Denmark has a dense population of Chironomus anthracinus, which survives 2–4 months of oxygen depletion each summer during stratification. The metabolism of 3rd and 4th instar larvae was examined in regard to variation in biomass and temperature. Respiration at air saturation was described by a curvilinear multiple regression relating oxygen consumption to individual AFDW and temperature. At 10 °C and varying oxygen regimes the O₂ consumption and CO₂ production of 4th instar larvae were almost unaltered from saturation to about 3 mg O₂ l^–1, but decreased steeply below this level. The respiratory quotient increased from 0.82 at saturation to about 3.4 at oxygen concentrations near 0.5 mg O₂ l^–1. This implied a shift from aerobic to partially anaerobic metabolism. At 0.5 mg O₂ l^–1 the total energy production equalled 20% of the rate at saturation of which more than one third was accounted for by anaerobic degradation of glycogen. This corresponded to a daily loss of 12 µg mg AFDW^–1 or approximately 5% of the body reserves. At unchanged metabolic rate the glycogen store would last three weeks, but long term oxygen deficiency causes a further suppression of the energy metabolism in C. anthracinus.     

Transpiration measurements on apple trees with an improved stem heat balance method

Since the late eighties a handy and user-friendly sap flow meter (Dynagage®) is on the market which can quantify 0205 the sap flow through intact plant stems, based on the stem heat balance method. The documentation about its accuracy and reliability, however, is still too limited to use it as a standard method in field experiments with apple trees. We therefore tested this commercial system on potted apple trees (Malus domestica L.; cv. Red Elstar and Jonagold; on rootstock M9 vf) with stem diameters of 1.8 to 4 cm. The measured sap flow was compared with mass loss measured by an automated balance, supposing the total mass loss of the trees was equal to the water loss by transpiration. The results revealed three major problems:

The Wageningen Rhizolab — a facility to study soil-root-shoot-atmosphere interactions in crops

A research facility is described for the integrated study of soil-root-shoot-atmosphere relationships in crops. The Wageningen Rhizolab has been in use since 1990, and consists of two rows, each with eight below-ground compartments aligned along a corridor. A rain shelter automatically covers the experimental area at the start of rainfall. Compartments are 125 cm × 125 cm and 200 cm deep. Each compartment has a separate drip irrigation system. Crop canopy photosynthesis, respiration, and transpiration can be measured simultaneously and continuously on four out of eight compartments at a time. Each compartment can be filled with a selected soil material (repacked soil) and is accessible from the corridor over its full depth. Multiple sensors for measuring soil moisture status, electrical conductivity, temperature, soil respiration, trace gases and oxygen are installed in spatial patterns in accordance with the requirements of the experiments. Sensors are connected to control and data-acquisition devices. Likewise, provisions have been made to sample manually the soil solution and soil atmosphere. Root observation tubes (minirhizotrons) are installed horizontally at depth intervals ranging from 5 cm (upper soil layers) to 25 cm (below 1 m). The facility is at present in use to study growth and development of vegetation (crops) in relation to drought, nutrient status, soil-borne diseases, and underground root competition. One important application is the study of elevated CO₂ concentration and climate change and the way they affect crops and their carbon economy. Growth and development of field grown vegetables and winter cover crops are also evaluated. The common aspect of those studies is to gain a better understanding of crop growth under varying environmental conditions, and to collect datasets that may help to improve mechanistic crop growth simulation models that can address suboptimal growth conditions.     

Is there an energy conservation “system” in brain that protects against the consequences of energy depletion?

A poorly understood marked decrease (circa 50% of control) in local cerebral glucose utilization is caused by sublethal doses of NaCN. The decrease is global, occurring in essentially all brain regions and is entirely reversible within hours, leaving no obvious pathology. This event is not unique to NaCN in so far as a strikingly similar pattern of decreased glucose utilization occus with some other toxins. Nor can it be attributed to a direct action of NaCN since local application by microdialysis to the striatum produces a global depression. These results imply that some widely distributed system or substance is involved. We speculate the existence of a system possibly related to the reticular activating system that senses a fall in energy production and acts globally to make cells quiescent and thus would give some protection from excitotoxic driven damage.Special issue dedicated to Dr. Sidney Ochs.     

Analytical study of the effects of some soil and plant parameters on root growth due to absorption of one mobile ion: A free-boundary model

The object of our study is: a model for root growth through a free-boundary problem and the effects resulting from differences in nutrient availability and transport of only one mobile nutrient between the root surface and the rhizosphere produced by an absorption Michaelis-Menten for low and high concentrations. The model equations are solved by two methods: the quasi-stationary method and the balance integral method. The numerical solutions are used to compute radial root growth. Curves of nutrient concentration at the root-soil interface, curve as a function of root radius as well as curves representing root radius as a function of time are plotted. The parameters which are varied are the root absorption power, flux velocity at the root surface, efflux, rhizosphere radius, diffusion coefficient, buffer power, and maximum influx. The two methods show the theoretical results for radial root growth in the range of low and high concentrations. The balance integral method provides more detailed information.     

The sequence of electron carriers in the reaction of cytochromec oxidase with oxygen

Kinetic studies of the electron transfer processes performed by cytochrome oxidase have assigned rates of electron transfer between the metal centers involved in the oxidation of ferrocytochromec by molecular oxygen. Transient-state studies of the reaction with oxygen have led to the proposal of a sequence of carriers from cytochromec, to Cu_A, to cytochromea, and then to the binuclear (i.e., cytochromea ₃-Cu_B) center. Electron exchange rates between these centers agree with relative center-to-center distances as follows; cytochromec to Cu_A 5–7 Å, cytochromec to cytochromea 20–25 Å, Cu_A to cytochromea 14–16 Å and cytochromea to cytochrome a₃-Cu_B 8–10 Å. It is proposed that the step from cytochromea to the binuclear center is the key control point in the reaction and that this step is one of the major points of energy transduction in the reaction cycle.