Using Lindenmayer systems to model morphogenesis in a tropical pasture legume Stylosanthes scabra

This paper summarizes the processes involved in designing a mathematical model of a growing pasture plant, Stylosanthes scabra Vog. cv. Fitzroy. The model is based on the mathematical formalism of Lindenmayer systems and yields realistic computer-generated images of progressive plant geometry through time. The processes involved in attaining growth data, retrieving useful growth rules, and constructing a virtual plant model are outlined. Progressive output morphological data proved useful for predicting total leaf area and allowed for easier quantification of plant canopy size in terms of biomass and total leaf area.

New integrable model of correlated electrons with off-diagonal long-range order from so(5) symmetry

We present a new integrable model for correlated electrons which is based on so(5) symmetry. By using an eta-pairing realization we construct eigenstates of the Hamiltonian with off-diagonal long-range order. It is also shown that these states lie in the ground state sector. We exactly solve the model on a one-dimensional lattice by the Bethe ansatz.

Protein sequence threading, the alignment problem, and a two-step strategy

Conventionally, protein structure prediction via threading relies on some nonoptimal method to align a protein sequence to each member of a library of known structures. We show how a score function (force field) can be modified so as to allow the direct application of a dynamic programming algorithm to the problem. This involves an approximation whose damage can be minimized by an optimization process during score function parameter determination. The method is compared to sequence to structure alignments using a more conventional pair-wise score function and the frozen approximation. The new method produces results comparable to the frozen approximation, but is faster and has fewer adjustable parameters. It is also free of memory of the template's original amino acid sequence, and does not suffer from a problem of nonconvergence, which can be shown to occur with the frozen approximation. Alignments generated by the simplified score function can then be ranked using a second score function with the approximations removed. (C) 1999 John Wiley & Sons, Inc.

Ameloblast apoptosis and IGF-1 receptor expression in the continuously erupting rat incisor model

Enamel-producing cells (ameloblasts) pass through several phenotypic and functional stages during enamel formation. In the transition between secretory and maturation stages, about one quarter of the ameloblasts suddenly undergo apoptosis. We have studied this phenomenon using the continuously erupting rat incisor model. A special feature of this model is that all stages of ameloblast differentiation are presented within a single longitudinal section of the developing tooth. This permits investigation of the temporal sequence of gene and growth factor receptor expression during ameloblast differentiation and apoptosis. We describe the light and electron microscopic morphology of ameloblast apoptosis and the pattern of insulin-like growth factor-1 receptor expression by ameloblasts in the continuously erupting rat incisor model. In the developing rat incisor, ameloblast apoptosis is associated with downregulated expression of the insulin-like growth factor-1 receptor. These data are consistent with the hypothesis that ameloblasts are hard wired for apoptosis and that insulin-like growth factor-1 receptor expression is required to block the default apoptotic pathway. Possible mechanisms of insulin-like growth factor-1 inhibition of ameloblast apoptosis are presented. The rat incisor model may be useful in studies of physiological apoptosis as it presents apoptosis in a predictable pattern in adult tissues.

RWMODEL II: Computer simulation of the Rescorla-Wagner model of Pavlovian conditioning

RWMODEL II simulates the Rescorla-Wagner model of Pavlovian conditioning. It is written in Delphi and runs under Windows 3.1 and Windows 95. The program was designed for novice and expert users and can be employed in teaching, as well as in research. It is user friendly and requires a minimal level of computer literacy but is sufficiently flexible to permit a wide range of simulations. It allows the display of empirical data, against which predictions from the model can be validated.

Class attainment among British men - A multivariate extension of the CASMIN model of intergenerational class mobility

The CASMIN Project is arguably the most influential contemporary study of class mobility in the world. However, CASMIN results with respect to weak vertical status effects on class mobility have been extensively criticized. Drawing on arguments about how to model vertical mobility, Hout and Hauser (1992) show that class mobility is strongly determined by vertical socioeconomic differences. This paper extends these arguments by estimating the CASMIN model while explicitly controlling for individual determinants of socioeconomic attainment. Using the 1972 Oxford Mobility Data and the 1979 and 1983 British Election Studies, the paper employs mixed legit models to show how individual socioeconomic factors and categorical differences between classes shape intergenerational mobility. The findings highlight the multidimensionality of class mobility and its irreducibility to vertical movement up and down a stratification hierarchy.

Determining forest structural attributes using an inverted geometric-optical model in mixed eucalypt forests, Southeast Queensland, Australia

The Montreal Process indicators are intended to provide a common framework for assessing and reviewing progress toward sustainable forest management. The potential of a combined geometrical-optical/spectral mixture analysis model was assessed for mapping the Montreal Process age class and successional age indicators at a regional scale using Landsat Thematic data. The project location is an area of eucalyptus forest in Emu Creek State Forest, Southeast Queensland, Australia. A quantitative model relating the spectral reflectance of a forest to the illumination geometry, slope, and aspect of the terrain surface and the size, shape, and density, and canopy size. Inversion of this model necessitated the use of spectral mixture analysis to recover subpixel information on the fractional extent of ground scene elements (such as sunlit canopy, shaded canopy, sunlit background, and shaded background). Results obtained fron a sensitivity analysis allowed improved allocation of resources to maximize the predictive accuracy of the model. It was found that modeled estimates of crown cover projection, canopy size, and tree densities had significant agreement with field and air photo-interpreted estimates. However, the accuracy of the successional stage classification was limited. The results obtained highlight the potential for future integration of high and moderate spatial resolution-imaging sensors for monitoring forest structure and condition. (C) Elsevier Science Inc., 2000.

Marine protected areas for spatially structured exploited stocks

Many harvested marine and terrestrial populations have segments of their range protected in areas free from exploitation. Reasons for areas being protected from harvesting include conservation, tourism, research, protection of breeding grounds, stock recovery, harvest regulation, or habitat that is uneconomical to exploit. In this paper we consider the problem of optimally exploiting a single species local population that is connected by dispersing larvae to an unharvested local population. We define a spatially-explicit population dynamics model and apply dynamic optimization techniques to determine policies for harvesting the exploited patch. We then consider how reservation affects yield and spawning stock abundance when compared to policies that have not recognised the spatial structure of the metapopulation. Comparisons of harvest strategies between an exploited metapopulation with and without a harvest refuge are also made. Results show that in a 2 local population metapopulation with unidirectional larval transfer, the optimal exploitation of the harvested population should be conducted as if it were independent of the reserved population. Numerical examples suggest that relative source populations should be exploited if the objective is to maximise spawning stock abundance within a harvested metapopulation that includes a protected local population. However, this strategy can markedly reduce yield over a sink harvested reserve system and may require strict regulation for conservation goals to be realised. If exchange rates are high, results indicate that spawning stock abundance can be less in a reserve system than in a fully exploited metapopulation. In order to maximise economic gain in the reserve system, results indicate that relative sink populations should be harvested. Depending on transfer levels, loss in harvest through reservation can be minimal, and is likely to be compensated by the potential environmental and economic benefits of the reserve.

A model for assembly and activation of the GM-CSF, IL-3 and IL-5 receptors: Insights from activated mutants of the common beta subunit

Granulocyte-macrophage colony stimulating factor (GM-CSF), Interleukin-3 (IL-3) and Interleukin-5 (IL-5) have overlapping, pleiotropic effects on hematopoietic cells, including neutrophils, eosinophils, monocytes and early progenitor cells. The high-affinity receptors for human GM-CSF, IL-3, and IL-5 share a common beta-subunit (h beta(c)), which is essential for signalling and plays a major role in recruiting intracellular signalling molecules. While activation of the cytoplasmic tyrosine kinase JAK2 appears to be the initiating event for signalling, the immediate events that trigger this are still unclear. We have isolated a number of activated mutants of h beta(c), which can be grouped into classes defined by their state of receptor phosphorylation, their requirement for alpha subunit as a cofactor, and their activities in primary cells and cell lines. We discuss these findings with regard to the stoichiometry, activation, and signalling of the normal GM-CSF/IL-3/IL-5 receptor complexes. Specifically, this work has implications for the role of the ligand-specific alpha-subunits in initiating the signalling through the beta-subunit, the role of beta subunit dimerization as a receptor trigger, and the function of receptor tyrosine phosphorylation in generating growth and survival signals. Based on the properties of the activated mutants and the recent structures of erythropoietin receptor (Epo-R) complexes, we propose a model in which (1) activation of h beta(c) can occur via alternative states that differ with respect to stoichiometry and subunit assembly, but which all mediate proliferative responses, and (2) each of the different classes of activated mutants mimics one of these alternative states. (C) 2000 International Society for Experimental Hematology. Published by Elsevier Science Inc.

Optimal release strategies for biological control agents: an application of stochastic dynamic programming to population management

1. Establishing biological control agents in the field is a major step in any classical biocontrol programme, yet there are few general guidelines to help the practitioner decide what factors might enhance the establishment of such agents. 2. A stochastic dynamic programming (SDP) approach, linked to a metapopulation model, was used to find optimal release strategies (number and size of releases), given constraints on time and the number of biocontrol agents available. By modelling within a decision-making framework we derived rules of thumb that will enable biocontrol workers to choose between management options, depending on the current state of the system. 3. When there are few well-established sites, making a few large releases is the optimal strategy. For other states of the system, the optimal strategy ranges from a few large releases, through a mixed strategy (a variety of release sizes), to many small releases, as the probability of establishment of smaller inocula increases. 4. Given that the probability of establishment is rarely a known entity, we also strongly recommend a mixed strategy in the early stages of a release programme, to accelerate learning and improve the chances of finding the optimal approach.

Finite element modelling of dissipative structures for nonequilibrium chemical reactions in fluid-saturated porous media

We use the finite element method to model and predict the dissipative structures of chemical species for a nonequilibrium chemical reaction system in a fluid-saturated porous medium. In particular, we explore the conditions under which dissipative structures of the species may exist in the Brusselator type of nonequilibrium chemical reaction. Since this is the first time the finite element method and related strategies have been used to study the chemical instability problems in a fluid-saturated porous medium, it is essential to validate the method and strategies before they are put into application. For this purpose, we have rigorously derived the analytical solutions for dissipative structures of chemical species in a benchmark problem, which geometrically is a square. Comparison of the numerical solutions with the analytical ones demonstrates that the proposed numerical method and strategy are robust enough to solve chemical instability problems in a fluid-saturated porous medium. Finally, the related numerical results from two application examples indicate that both the regime and the magnitude of pore-fluid flow have significant effects on the nature of the dissipative structures that developed for a nonequilibrium chemical reaction system in a fluid-saturated porous medium. The motivation for this study is that self-organization under conditions of pore-fluid flow in a porous medium is a potential mechanism of the orebody formation and mineralization in the upper crust of the Earth. (C) 2000 Elsevier Science S.A. All rights reserved.

The role of mobile belts for the longevity of deep cratonic lithosphere: The crumple zone model

Many Archean cratons are surrounded by Proterozoic mobile belts that have experienced episodes of tectonic re-activation over their lifetimes. This suggests that mobile belt lithosphere may be associated with long lived, inherited weakness. It is proposed that the proximity of this weakness can increase the longevity of deep Archean lithosphere by buffering Archean cratons from mantle derived stresses. The physical plausibility of this idea is explored through numerical simulations of mantle convection that include continents and allow for material rheologies that model the combined brittle and ductile behavior of the lithosphere. Within the simulations, the longevity of deep cratonic lithosphere does increase if it is buffered by mobile belts that can fail at relatively low stress levels.

Modeling of hepatic elimination and organ distribution kinetics with the extended convection-dispersion model

The conventional convection-dispersion (also called axial dispersion) model is widely used to interrelate hepatic availability (F) and clearance (Cl) with the morphology and physiology of the liver and to predict effects such as changes in liver blood flow on F and Cl. An extended form of the convection-dispersion model has been developed to adequately describe the outflow concentration-time profiles for vascular markers at both short and long times after bolus injections into perfused livers. The model, based on flux concentration and a convolution of catheters and large vessels, assumes that solute elimination in hepatocytes follows either fast distribution into or radial diffusion in hepatocytes. The model includes a secondary vascular compartment, postulated to be interconnecting sinusoids. Analysis of the mean hepatic transit time (MTT) and normalized variance (CV2) of solutes with extraction showed that the discrepancy between the predictions of MTT and CV2 for the extended and conventional models are essentially identical irrespective of the magnitude of rate constants representing permeability, volume, and clearance parameters, providing that there is significant hepatic extraction. In conclusion, the application of a newly developed extended convection-dispersion model has shown that the unweighted conventional convection-dispersion model can be used to describe the disposition of extracted solutes and, in particular, to estimate hepatic availability and clearance in booth experimental and clinical situations.

The Fornax spectroscopic survey I. Survey strategy and preliminary results on the redshift distribution of a complete sample of stars and galaxies

The Fornax Spectroscopic Survey will use the Two degree Field spectrograph (2dF) of the Angle-Australian Telescope to obtain spectra for a complete sample of all 14000 objects with 16.5 less than or equal to b(j) less than or equal to 19.7 in a 12 square degree area centred on the Fornax Cluster. The aims of this project include the study of dwarf galaxies in the cluster (both known low surface brightness objects and putative normal surface brightness dwarfs) and a comparison sample of background field galaxies. We will also measure quasars and other active galaxies, any previously unrecognised compact galaxies and a large sample of Galactic stars. By selecting all objects-both stars and galaxies-independent of morphology, we cover a much larger range of surface brightness and scale size than previous surveys. In this paper we first describe the design of the survey. Our targets are selected from UK Schmidt Telescope sky survey plates digitised by the Automated Plate Measuring (APM) facility. We then describe the photometric and astrometric calibration of these data and show that the APM astrometry is accurate enough for use with the 2dF. We also describe a general approach to object identification using cross-correlations which allows us to identify and classify both stellar and galaxy spectra. We present results from the first 2dF field. Redshift distributions and velocity structures are shown for all observed objects in the direction of Fornax, including Galactic stars? galaxies in and around the Fornax Cluster, and for the background galaxy population. The velocity data for the stars show the contributions from the different Galactic components, plus a small tail to high velocities. We find no galaxies in the foreground to the cluster in our 2dF field. The Fornax Cluster is clearly defined kinematically. The mean velocity from the 26 cluster members having reliable redshifts is 1560 +/- 80 km s(-1). They show a velocity dispersion of 380 +/- 50 km s(-1). Large-scale structure can be traced behind the cluster to a redshift beyond z = 0.3. Background compact galaxies and low surface brightness galaxies are found to follow the general galaxy distribution.