On systems thinking, systems biology and the in Silico Plant

The recent summary report of a Department of Energy Workshop on Plant Systems Biology (P.V. Minorsky [2003] Plant Physiol 132: 404-409) offered a welcomed advocacy for systems analysis as essential in understanding plant development, growth, and production. The goal of the Workshop was to consider methods for relating the results of molecular research to real-world challenges in plant production for increased food supplies, alternative energy sources, and environmental improvement. The rather surprising feature of this report, however, was that the Workshop largely overlooked the rich history of plant systems analysis extending over nearly 40 years (Sinclair and Seligman, 1996) that has considered exactly those challenges targeted by the Workshop. Past systems research has explored and incorporated biochemical and physiological knowledge into plant simulation models from a number of perspectives. The research has resulted in considerable understanding and insight about how to simulate plant systems and the relative contribution of various factors in influencing plant production. These past activities have contributed directly to research focused on solving the problems of increasing biomass production and crop yields. These modeling approaches are also now providing an avenue to enhance integration of molecular genetic technologies in plant improvement (Hammer et al., 2002).

Advanced computing for systems biology

Systems biology is based on computational modelling and simulation of large networks of interacting components. Models may be intended to capture processes, mechanisms, components and interactions at different levels of fidelity. Input data are often large and geographically disperse, and may require the computation to be moved to the data, not vice versa. In addition, complex system-level problems require collaboration across institutions and disciplines. Grid computing can offer robust, scaleable solutions for distributed data, compute and expertise. We illustrate some of the range of computational and data requirements in systems biology with three case studies: one requiring large computation but small data (orthologue mapping in comparative genomics), a second involving complex terabyte data (the Visible Cell project) and a third that is both computationally and data-intensive (simulations at multiple temporal and spatial scales). Authentication, authorisation and audit systems are currently not well scalable and may present bottlenecks for distributed collaboration particularly where outcomes may be commercialised. Challenges remain in providing lightweight standards to facilitate the penetration of robust, scalable grid-type computing into diverse user communities to meet the evolving demands of systems biology.

Inherent size constraints on prokaryote gene networks due to "accelerating" growth

Networks exhibiting accelerating growth have total link numbers growing faster than linearly with network size and either reach a limit or exhibit graduated transitions from nonstationary-to-stationary statistics and from random to scale-free to regular statistics as the network size grows. However, if for any reason the network cannot tolerate such gross structural changes then accelerating networks are constrained to have sizes below some critical value. This is of interest as the regulatory gene networks of single-celled prokaryotes are characterized by an accelerating quadratic growth and are size constrained to be less than about 10,000 genes encoded in DNA sequence of less than about 10 megabases. This paper presents a probabilistic accelerating network model for prokaryotic gene regulation which closely matches observed statistics by employing two classes of network nodes (regulatory and non-regulatory) and directed links whose inbound heads are exponentially distributed over all nodes and whose outbound tails are preferentially attached to regulatory nodes and described by a scale-free distribution. This model explains the observed quadratic growth in regulator number with gene number and predicts an upper prokaryote size limit closely approximating the observed value. (c) 2005 Elsevier GmbH. All rights reserved.

Accelerated leap methods for simulating discrete stochastic chemical kinetics

Biologists are increasingly conscious of the critical role that noise plays in cellular functions such as genetic regulation, often in connection with fluctuations in small numbers of key regulatory molecules. This has inspired the development of models that capture this fundamentally discrete and stochastic nature of cellular biology - most notably the Gillespie stochastic simulation algorithm (SSA). The SSA simulates a temporally homogeneous, discrete-state, continuous-time Markov process, and of course the corresponding probabilities and numbers of each molecular species must all remain positive. While accurately serving this purpose, the SSA can be computationally inefficient due to very small time stepping so faster approximations such as the Poisson and Binomial τ-leap methods have been suggested. This work places these leap methods in the context of numerical methods for the solution of stochastic differential equations (SDEs) driven by Poisson noise. This allows analogues of Euler-Maruyuma, Milstein and even higher order methods to be developed through the Itô-Taylor expansions as well as similar derivative-free Runge-Kutta approaches. Numerical results demonstrate that these novel methods compare favourably with existing techniques for simulating biochemical reactions by more accurately capturing crucial properties such as the mean and variance than existing methods.

Analysis and application of an equilibrium model for in vitro bioassay systems with three components: Receptor, hormone and hormone-binding-protein

A simple theoretical framework is presented for bioassay studies using three component in vitro systems. An equilibrium model is used to derive equations useful for predicting changes in biological response after addition of hormone-binding-protein or as a consequence of increased hormone affinity. Sets of possible solutions for receptor occupancy and binding protein occupancy are found for typical values of receptor and binding protein affinity constants. Unique equilibrium solutions are dictated by the initial condition of total hormone concentration. According to the occupancy theory of drug action, increasing the affinity of a hormone for its receptor will result in a proportional increase in biological potency. However, the three component model predicts that the magnitude of increase in biological potency will be a small fraction of the proportional increase in affinity. With typical initial conditions a two-fold increase in hormone affinity for its receptor is predicted to result in only a 33% increase in biological response. Under the same conditions an Ii-fold increase in hormone affinity for receptor would be needed to produce a two-fold increase in biological potency. Some currently used bioassay systems may be unrecognized three component systems and gross errors in biopotency estimates will result if the effect of binding protein is not calculated. An algorithm derived from the three component model is used to predict changes in biological response after addition of binding protein to in vitro systems. The algorithm is tested by application to a published data set from an experimental study in an in vitro system (Lim et al., 1990, Endocrinology 127, 1287-1291). Predicted changes show good agreement (within 8%) with experimental observations. (C) 1998 Academic Press Limited.

Visual biology of Hawaiian coral reef fishes. II. Colors of Hawaiian coral reef fish

The colors of 51 species of Hawaiian reef fish have been measured using a spectrometer and therefore can be described in objective terms that are not influenced by the human visual experience. In common with other known reef fish populations, the colors of Hawaiian reef fish occupy spectral positions from 300-800nm; yellow or orange with blue, yellow with black, and black with white are the most frequently combined colors; and there is no link between possession of ultraviolet (UV) reflectance and UV visual sensitivity or the potential for UV visual sensitivity. In contrast to other reef systems, blue, yellow, and orange appear more frequently in Hawaiian reef fish. Based on spectral quality of reflections from fish skin, trends in fish colors can be seen that are indicative of both visually driven selective pressures and chemical or physical constraints on the design of colors. UV-reflecting colors can function as semiprivate communication signals. White or yellow with black form highly contrasting patterns that transmit well through clear water. Labroid fishes display uniquely complex colors but lack the ability to see the UV component that is common in their pigments. Step-shaped spectral curves are usually long-wavelength colors such as yellow or red, and colors with a peak-shaped spectral curves are green, blue, violet, and UV.

Muscle-specific regulation of tropomyosin gene expression and myofibrillogenesis differs among muscle systems examined at metamorphosis of the gastropod Haliotis rufescens

The spatial and temporal association of muscle-specific tropomyosin gene expression, and myofibril assembly and degradation during metamorphosis is analyzed in the gastropod mollusc. Haliotis rufescens. Metamorphosis of tile planktonic larva to the benthic juvenile includes rearrangement and atrophy of specific larval muscles, and biogenesis of the new juvenile muscle system. The major muscle of the larva - the larval retractor muscle - reorganizes at metamorphosis, with two suites of cells having different fates. The ventral cells degenerate, while the dorsal cells become part of the developing juvenile mantle musculature. Prior to these changes in myofibrillar structure, tropomyosin mRNA prevalence declines until undetectable in the ventral cells, while increasing markedly in the dorsal cells. In the foot muscle and right shell muscle, tropomyosin mRNA levels remain relatively stable, even trough myofibril content increases. In a population of median mesoderm cells destined to form de novo the major muscle of the juvenile and adult (the columellar muscle), tropomyosin expression is initiated at 45 h after induction of metamorphosis. Myofibrillar filamentous actin is not detected in these cells until about 7 days later. Given that patterns of tropomyosin mRNA accumulation in relation to myofibril assembly and disassembly differ significantly among the four major muscle systems examined, we suggest that different regulatory mechanisms, probably operating at both transcriptional and post-transcriptional levels, control the biogenesis and atrophy of different larval and postlarval muscles at metamorphosis.

Differential perturbation of neuronal and glial glutamate transport systems in retinal ischaemia

Glutamate is the major excitatory neurotransmitter in the retina and is removed from the extracellular space by an energy-dependent process involving neuronal and glial cell transporters. The radial glial Muller cells express the glutamate transporter, GLAST, and preferentially accumulate glutamate. However, during an ischaemic episode, extracellular glutamate concentrations may rise to excitotoxic levels. Is this catastrophic rise in extracellular glutamate due to a failure of GLAST? Using immunocytochemistry, we monitored the transport of the glutamate transporter substrate, D-aspartate, in the retina under normal and ischaemic conditions. Two models of compromised retinal perfusion were compared: (1) Anaesthetised rats had their carotid arteries occluded for 7 days to produce a chronic reduction in retinal blood flow. Retinal function was assessed by electroretinography. D-aspartate was injected into the eye for 45 min, Following euthanasia, the retina was processed for D-aspartate. GLAST and glutamate immunocytochemistry. Although reduced retinal perfusion suppresses the electroretinogram b-wave, neither retinal histology, GLAST expression, nor the ability of Muller cells to uptake D-aspartate is affected. As this insult does not appear to cause excitotoxic neuronal damage, these data suggest that GLAST function and glutamate clearance are maintained during periods of reduced retinal perfusion. (2) Occlusion of the central retinal artery for 60 min abolishes retinal perfusion, inducing histological damage and electroretinogram suppression. Although GLAST expression appears to be normal. its ability to transport D-aspartate into Muller cells is greatly reduced. Interestingly, D-aspartate is transported into neuronal cells, i.e. photoreceptors, bipolar and ganglion cells. This suggests that while GLAST is vitally important for the clearance of excess extracellular glutamate, its capability to sustain inward transport is particularly susceptible to an acute ischaemic attack. Manipulation of GLAST function could alleviate the degeneration and blindness that result from ischaemic retinal disease. (C) 2001 Elsevier Science Ltd, All rights reserved.

Lipid, sugar and liposaccharide based delivery systems

Although there are formidable barriers to the oral delivery of biologically active drugs, considerable progress in the field has been made, using both physical and chemical strategies of absorption enhancement. A possible method to enhance oral absorption is to exploit the phenomenon of lipophilic modification and mono and oligosaccharide conjugation. Depending on the uptake mechanism targeted, different modifications can be employed. To target passive diffusion, lipid modification has been used, whereas the targeting of sugar transport systems has been achieved through drugs conjugated with sugars. These drug delivery units can be specifically tailored to transport a wide variety of poorly absorbed drugs through the skin, and across the barriers that normally inhibit absorption from the gut or into the brain. The delivery system can be conjugated to the drug in such a way as to release the active compound after it has been absorbed (i.e. the drug becomes a prodrug), or to form a biologically stable and active molecule (i.e. the conjugate becomes a new drug moiety). Examples where lipid, sugar and lipid-sugar conjugates have resulted in enhanced drug delivery will be highlighted in this review.

Phase variable restriction-modification systems in Moraxella catarrhalis

A repetitive DNA motif was used as a marker to identify novel genes in the mucosal pathogen Moraxella catarrhalis. There is a high prevalence of such repetitive motifs in virulence genes that display phase variable expression. Two repeat containing loci were identified using a digoxigenin-labelled 5'-(CAAC)(6)-3' oligonucleotide probe. The repeats are located in the methylase components of two distinct type III restriction-modification (R-M) systems. We suggest that the phase variable nature of these R-M systems indicates that they have an important role in the biology of M. catarrhalis. (C) 2002 Published by Elsevier Science B.V. on behalf of the Federation of European Microbiological Societies.

Partial automation of database processing of simulation outputs from L-systems models of plant morphogenesis

Models of plant architecture allow us to explore how genotype environment interactions effect the development of plant phenotypes. Such models generate masses of data organised in complex hierarchies. This paper presents a generic system for creating and automatically populating a relational database from data generated by the widely used L-system approach to modelling plant morphogenesis. Techniques from compiler technology are applied to generate attributes (new fields) in the database, to simplify query development for the recursively-structured branching relationship. Use of biological terminology in an interactive query builder contributes towards making the system biologist-friendly. (C) 2002 Elsevier Science Ireland Ltd. All rights reserved.