Diurnal water storage in the stems of Picea sitchensis (Bong.) Carr.

Abstract. Two models of the relationship between diurnal variation in shoot water potential and transpiration in 14-year-old Picea sitchensis (Bong.) Carr. were compared. The first model was a physiologically based resistance-capacitance (R-C) analogue with its associated differential equations. The second was a non-physiological discrete-difference (D-D) or stochastic transfer function model. The RC model included only the effect of water storage in the phloem and bark while the D-D model implicity included all storage mechanisms. The R-C and D-D models explained similar fractions (62% and 68% respectively) of the variation in shoot water potential due to diurnal changes in transpiration rate. However, the D-D model had fewer parameters than the R-C model. The results from the D-D model showed that the resistance to flow from soil to shoots along the trunk, (RT), was 5 × 10³ MPa kg^-1s and the capacitance of the phloem and bark treated as a single store, (C_s), was 1.6 kg MPa^-1. It is suggested that the resistance to flow into storage (R_s) is much greater than R_T and can be disregarded. A non-linear version of the D-D model suggested [hat resistance to flow in the trunk increases with increasing transpiration rate.     

On the solution of mathematical models of herd immunity in human helminth infections

The general solution of the mathematical model of herd immunity to human helminth infections recently proposed by Anderson and May [3] is obtained. The numerical solution of a more accurate biological model is indistinguishable from the corresponding exact solution of a more tractable mathematical model. Computer simulations of some particular cases of this model support the notion that both ecological and immunological factors determine the observed convex patterns of age-prevalence and age-intensity curves of human helminth infections.This work was made thanks to the advise and support of Dr. Robert M. May while the author was Postdoctoral Fellow at Princeton University     

Relationships among recent models for insect population dynamics with variable rates of development

A classification scheme for those population models which allow variation in development rates is proposed, based on two ways of modifying standard age-structured models. The resulting classes of models are termed development index models and sojourn time models. General formulations for the two classes of models are developed from two basic balance equations, and numerous specific models from the literature are shown to fit into the scheme. Concepts from competing risks theory are shown to be important in understanding the interplay between mortality and maturation. Relationships among the classes are investigated both for the most general forms of the models and for the simpler forms often used. The scheme can provide guidance in developing appropriate insect population models for specific modelling situations.Contribution 3878871     

A simulation of the process of sedimentation of suspended solids in the Yoshii River estuary

A vertical plane 2D model has been applied to the Yoshii River estuary, Japan, as a means of simulating the process of sedimentation. The current velocity and isohaline distributions were well reproduced in the estuary and the model simulated the concentration distribution of suspended solids for 8 particle size classes by using a transport equation based only on physical processes. The changes of the particle size distribution in surface water were well reproduced by the model. The vertical velocity component plays an important role in the behaviour of suspended solids.     

Models for contact-mediated pattern formation: cells that form parallel arrays

Kinetic continuum models are derived for cells that crawl over a 2D substrate, undergo random reorientation, and turn in response to contact with a neighbor. The integro-partial differential equations account for changes in the distribution of orientations in the population. It is found that behavior depends on parameters such as total mass, random motility, adherence, and sloughing rates, as well as on broad aspects of the contact response. Linear stability analysis, and numerical, and cellular automata simulations reveal that as parameters are varied, a bifurcation leads to loss of stability of a uniform (isotropic) steady state, in favor of an (anisotropic) patterned state in which cells are aligned in parallel arrays.     

Conformational and helicoidal analysis of the molecular dynamics of proteins: "curves," dials and windows for a 50 psec dynamic trajectory of BPTI

A new procedure for the graphic analysis of molecular dynamics (MD) simulations on proteins is introduced, in which comprehensive visualization of results and pattern recognition is greatly facilitated. The method involves determining the conformational and helicoidal parameters for each structure entering the analysis via the method "Curves," developed for proteins by Sklenar, Etchebest, and Lavery (Proteins: Structure, Function Genet. 6:46-60, 1989) followed by a novel computer graphic display of the results. The graphic display is organized systematically using conformation wheels ("dials") for each torsional parameter and "windows" on the range values assumed by the linear and angular helicoidal parameters, and is present in a form isomorphous with the primary structure per se. The complete time evolution of dynamic structure can then be depicted in a set of four composite figures. Dynamic aspects of secondary and tertiary structure are also provided. The procedure is illustrated with an analysis of a 50 psec in vacuo simulation on the 58 residue protein, bovine pancreatic trypsin inhibitor (BPTI), in the vicinity of the local minimum on the energy surface corresponding to a high resolution crystal structure. The time evolution of 272 conformational and 788 helicoidal parameters for BPTI is analyzed. A number of interesting features can be discerned in the analysis, including the dynamic range of conformational and helicoidal motions, the dynamic extent of 2 degrees structure motifs, and the calculated fluctuations in the helix axis. This approach is expected to be useful for a critical analysis of the effects of various assumptions about force field parameters, truncation of potentials, solvation, and electrostatic effects, and can thus contribute to the development of more reliable simulation protocols for proteins. Extensions of the analysis to present differential changes in conformational and helicoidal parameters is expected to be valuable in MD studies of protein complexes with substrates, inhibitors, and effectors and in determining the nature of structural changes in protein-protein interactions.     

Total root length as estimated from small sub-samples

A method is proposed for estimating the total length of a root system from sub-samples. This method is based on the measurement of the length and diameter of small pieces of roots, and on the measurements of the bulk density of root sub-samples. It is assumed that roots are cylinders with a given bulk density. The length and diameter of small root pieces are measured by image analysis. A weighted quadratic mean (W.Q.M.) root diameter is then calculated and used in estimating the root length. This W.Q.M. diameter is defined as the real mean diameter of an equivalent single root with the same length and volume as the tested root system. The accuracy of prediction is demonstrated for one theoretical root system. The standard deviation of estimation can be calculated using sampling simulations.     

Correlation of immunomodulatory and therapeutic activities of interferon and interferon inducers in metastatic disease

The mechanism of therapeutic activity of recombinant murine interferon-gamma (rMu IFN-gamma) and the IFN inducer polyinosinic-polycytidylic acid solubilized with poly-L-lysine in carboxy methyl cellulose (pICLC) in treating metastatic disease was investigated by comparing effector cell augmentation with therapeutic activity in mice bearing experimental lung metastases (B16-BL6 melanoma). Effector cell functions in spleen, peripheral blood, and lung (the organ with tumor) were tested after 1 and 3 weeks of rMu IFN-gamma or pICLC administration (intravenous, three times a week). In these studies, natural killer (NK), lymphokine-activated killer (LAK), cytolytic T lymphocytes (CTL) (against specific and nonspecific targets), and macrophage tumoricidal and tumoristatic activities were measured. rM IFN-gamma and pICLC had therapeutic activity and immunomodulatory activity in most assays of immune function examined. Specific CTL activity of pulmonary parenchymal mononuclear cells (PPMC), but not in splenocytes or peripheral blood lymphocytes (PBL), during week 3 and not during week 1, correlated with the therapeutic activity of rMu IFN-gamma and of pICLC. Macrophage tumoricidal activity in PPMC, but not in alveolar macrophages, also correlated with the therapeutic activity of rMu IFN-gamma, but the opposite was true for the therapeutic activity of pICLC. NK activity of PPMC, but not of splenocytes or PBL, during week 1 correlated with the therapeutic activity of pICLC; in contrast, NK activity at any site did not correlate with the therapeutic activity of rMu IFN-gamma. LAK activity at any site did not correlate with the therapeutic activity of either agent.     

The effect of time-varying mortality and carbon assimilation on models of carbon allocation in annual plants

Summary Models of optimal carbon allocation schedules have influenced the way plant ecologists think about life history evolution, particularly for annual plants. The present study asks (1) how, within the framework of these models, are their predictions affected by within-season variation in mortality and carbon assimilation rates?; and (2) what are the consequences of these prediction changes for empirical tests of the models? A companion paper examines the basic assumptions of the models themselves. I conducted a series of numerical experiments with a simple carbon allocation model. Results suggest that both qualitative and quantitative predictions can sometimes be sensitive to parameter values for net assimilation rate and mortality: for some parameter values, both the time and size at onset of reproduction, as well as the number of reproductive intervals, vary considerably as a result of small variations in these parameters. For other parameter values, small variations in the parameters result in only small changes in predicted phenotype, but these have very large fitness consequences. Satisfactory empirical tests are thus likely to require much accuracy in parameter estimates. The effort required for parameter estimation imposes a practical constraint on empirical tests, making large multipopulation comparisons impractical. It may be most practical to compare the predicted and observed fitness consequences of variation in the timing of onset of reproduction.     

Scaling physiological measurements for individuals of different body size

This paper examines how selected physiological performance variables, such as maximal oxygen uptake, strength and power, might best be scaled for subject differences in body size. The apparent dilemma between using either ratio standards or a linear adjustment method to scale was investigated by considering how maximal oxygen uptake (l.min-1), peak and mean power output (W) might best be adjusted for differences in body mass (kg). A curvilinear power function model was shown to be theoretically, physiologically and empirically superior to the linear models. Based on the fitted power functions, the best method of scaling maximum oxygen uptake, peak and mean power output, required these variables to be divided by body mass, recorded in the units kg 2/3. Hence, the power function ratio standards (ml.kg-2/3.min-1) and (W.kg-2/3) were best able to describe a wide range of subjects in terms of their physiological capacity, i.e. their ability to utilise oxygen or record power maximally, independent of body size. The simple ratio standards (ml.kg-1.min-1) and (W.kg-1) were found to best describe the same subjects according to their performance capacities or ability to run which are highly dependent on body size. The appropriate model to explain the experimental design effects on such ratio standards was shown to be log-normal rather than normal. Simply by taking logarithms of the power function ratio standard, identical solutions for the design effects are obtained using either ANOVA or, by taking the unscaled physiological variable as the dependent variable and the body size variable as the covariate, ANCOVA methods.