ALBERT

Notizen: Abstract In the theory of organismic sets (Bull. Math. Biophysics,31, 159–198, 1969) we considered organisms as sets endowed with certain “activities,” the latter’s resulting in a set of “products.” Those products may be of a material nature, like a hormone secreted by a cell, or of a non-material nature, like a feeling or an attitude. In the present paper aggressiveness and submissiveness are considered as such non-material products of the activities of the brain cells. A general description of aggressiveness and submissiveness is given in terms of organismic sets. Cycles in “peck order” are thus naturally explained.

Notizen: Abstract In line with previous studies on organismic sets, the division of all organismic sets intogeneral autotrophic and heterotrophic is introduced. The first produce their food themselves from some external source of energy, which in general may be an energy of any kind. The others use other organismic sets as the source of their food and energy. On earth we know only one kind of generalgeneral autotrophic organismic sets, namely, the autotrophic plants which use solar radiation as their source of energy and for production of their own food. It is shown why autotrophic animals do not exist on earth except as microorganisms like, e.g.,Euglena. A rigorous proof of the previously derived theorem that in an organismic set of ordern〉1 no element can be completely specialized is given. It requires the introduction of new postulates. Finally, in considering the organic world as a whole, the notion of organismic sets ofmixed order is introduced.

Notizen: Abstract It appears to be axiomatic that termolecular and higher order reactions occur relatively rarely. The basis for this judgment seems to lie in the supposition that successful 3-Body collisions of 3 interactive species of molecules cannot occur frequently enought to account for chemical or biochemical transformation. In order to provide a more complete mathematical framework than now exists for examining this hypothesis the probability of effective termolecular “δ-collisions” as a function of time is derived. This amounts to adding to the class of reactions for which stochastic models are now available the termolecular reaction. In common with the unimolecular and bimolecular cases this process is seen to satisfy the criterion of consistency-in-the-mean with respect to deterministic formulations. It is planned next to use the termolecular process and the lower order processes in computer-assistedin numero experimental studies aimed at comparing alternative mechanisms of reaction.

Notizen: Abstract Small sample properties of the maximum likelihood estimator for the rate constant of a stochastic first order reaction are investigated. The approximate bias and variance of the maximum likelihood estimator are derived and tabulated. If observations of the system are made at timesiτ,i=1, 2, ...,N; τ〉0, the observational spacing τ which minimizes the approximate variance of the maximum likelihood estimator is found. The non-applicability of large sample theory to confidence interval derivation is demonstrated by examination of the relative likelihood. Bartlett’s method is employed to derive approximate confidence limits, and is illustrated by using simulated kinetic runs.

Notizen: Abstract The representation of biological systems by means of organismic supercategories, developed in previous papers (Bull. Math. Biophysics,30, 625–636;31, 59–71;32, 539–561), is further discussed. The different approaches to relational biology, developed by Rashevsky, Rosen and by Băianu and Marinescu, are compared with Qualitative Dynamics of Systems which was initiated by Henri Poincaré (1881). On the basis of this comparison some concrete result concerning dynamics of genetic system, development, fertilization, regeneration, analogies, and oncogenesis are derived.

Notizen: Abstract Some years ago (Rosen 1958a, b; 1959) we described a class of metaphorical, relational paradigms for cellular activity which we termed (M, R)-systems. A sizable amount of subsequent work, to be itemized below, has been devoted to an exploration of some of the properties of these systems. The main purpose of the present paper is to put this class of paradigms into a general system-theoretic perspective, with a particular view to appraising the relation between the type of system description embodied in the (M, R)-system and other kinds of physical and mathematical descriptions of cellular systems. Thus, the principal aim is to establish the relationships and connections between the global relational formalism embodied in the (M, R)-systems and the empirical descriptions which still represent the bulk of our biological knowledge.

Notizen: Abstract After giving a brief review of the theory of organismic sets (Bull. Math. Biophysics,29, 139–152, 1967;31, 159–198, 1969), in which the concept of relational forces, introduced earlier (Bull. Math. Biophysics,28, 283–308, 1966a) plays a fundamental role, the author discusses examples of possible different structures produced by relational forces. For biological organisms the different structures found theoretically are in general agreement with observation. For societies, which are also organismic sets as discussed in the above references, the structures can be described only in an abstract space, the nature of which is discussed. Different isomorphisms between anatomical structures, as described in ordinary Euclidean space, and the sociological structures described in an abstract space are noted, as should be expected from the theory of organismic sets.

Notizen: Abstract Current psychological research into the inference (diagnostic) process is briefly reviewed, using as a vehicle an investigation of the prediction of the probability of success of hypothetical applicants to a graduate program in biology. Brunswik’s lens model and multiple regression analysis are used, as is a Bayesian approach. Four judges’ (biologists’) predictions are analyzed. Some general conclusions about inference, drawn from the current data in psychology, are presented.

Notizen: Abstract A mathematical model has been developed to simulate the glucose-insulin interaction following a glucose load such as occurs in an IVGTT. This model differs from earlier models in that the insulin response to glucose loading is a recurring all or none threshold response. The model has been simulated on a digital computer using the digital analog simulation language CSMP.

Notizen: Abstract The composite nature of bone dictates the use of a model for bone which is transversely isotropic. We solve the associated sets of partial differential equations governing the dynamic elastic behavoor of a two-layered cylindrical-shaped bone. The solution is analyzed for long, short, and intermediate length waves. The special case of compact bone is treated for long and short wave lengths and a numerical example is worked out to determine the wave speeds (for short wave lengths) given a set of elastic constants, determined by ultrasonic methods, and the bone density, wave frequency, and radius.

Notizen: Abstract A somewhat different approach to the principle of biotopological mapping, discussed in previous publications, is given. The organism is considered as a set of properties, each of which is in its turn a set of numerous subproperties which are logically included in the corresponding properties. Topology is introduced by an appropriate definition of neighborhoods, and four postulates are stated which concern the mapping of the spaces corresponding to higher organisms on those of lower ones. A number of conclusions are drawn from the postulates. Some of them correspond to well-known facts. For example, in man and some higher organisms appropriate emotional stimuli should produce gastrointestinal or cardiovascular disturbances; or some microorganisms should produce substances harmful to other microorganisms (antibiotics). Some other conclusions are still awaiting verification. One of them is, for example, that there must exist unicellular organisms which produce antibodies to appropriate antigens.

Notizen: Zusammenfassung Für die Praxis der Pflanzenernährung ist es wichtig, zu wissen, in welcher Weise die Ertragsbildung von der Konzentration eines mineralischen Nährstoffes in der Umgebung der Pflanze abhängt. Da nur diejenigen Nährstoffmengen das physiologische Geschehen in der Pflanze unmittelbar zu beeinflussen vermögen, die sich in der Pflanze befinden, wird angenommen, dass das Wachstum zum Zeitpunktt, d.h. die Geschwindigkeit der Trockensubstanzzunahme zu diesem Zeitpunkt, eine Funktion der zur Zeitt in der Pflanze enthaltenen Nährstoffmenge ist. Diese Nährstoffmenge wird natürlich im Intervall vor dem Zeitpunktt aufgenommen. Deshalb und auch noch aus anderen Gründen hängt das Wachstum zur Zeitt davon ab, wie die in der Umgebung der Pflanze herrschende Konzentration des betrachteten Nährstoffes in demjenigen Zeitintervall verläuft, das sich von der Aussaat bis zum Zeitpunktt erstreckt. Die angegebene Annahme fürhrt zusammen mit einigen weiteren naheliegenden Annahmen zu einem Ansatz, der Ergebnisse liefert, die in verschiedener Hinsicht gut mit der Erfahrung übereinstimmen. Jedoch gibt es auch noch Widersprüche zwischen Theorie und Erfahrung. Durch weitere Ausgestaltung der Theorie lassen sich diese Widersprüche beseitigen. Es wird angeregt, Versuche durchzuführen, deren Resultate Hinweise für die weitere Ausgestaltung der Theorie liefern.

Notizen: Abstract The dynamics of cell multiplication and differentiation in tissues in asteady state and the kinetics of isotope incorporation into the DNA have been theoretically analyzed. Equations have been derived, with the aid of which thegeneration time, thelife span, and the distribution or rate of death of the cells can be obtained if the tissue is in asteady state, i.e., if the number of cells is maintained constant by constant, equal rates of cell division and cell death and if the mean DNA content per cell is also constant. An equation has also been derived which gives thegeneration time in the case of logarithmic multiplication of cells. Two special cases have been analyzed: InCase 1, the isotope is considered as being introduced into the metabolic system at zero time only; inCase 2, the specific activity of the DNA precursor is considered as being maintained constant. The use of the method has been illustrated by an example in which thegeneration time and themean, themedian, and themode life span, as well as the curve of the rate of death of leukocytes in a patient with chronic leukemic granulocytic leukemia, have been obtained from the rate of P32 incorporation into the DNA. The merits and the limitations of the method are discussed.

Notizen: Abstract A stochastic model is developed for a compartment with a single time-dependent input, and generalized to include inputs from several sources. With the number of particles of a given molecular species in the compartment as the random variable, the mean, variance and third central moment of this variable are calculated from its generating function, and compared with previous results. The behavior of the calculated moments is discussed, and the possibility of applying the model to chemical and biological systems is considered.

Notizen: Abstract By means of the Laplace transform, the behavior of a simplified model of the cardiovascular system is mathematically formulated. This formulation allows mathematical expression of the periodicity of the cardiac output and the systemic response. With the cardiac output represented as half of a sine function cycle, the systolic aortic pressure becomes the sum of a sine term and exponential terms, while the sum of the exponential terms alone represents the diastolic pressure. The characteristics of the mathematical expressions for systole and diastole are analyzed, and some relationships of potentially practical value are derived. Variation in the parameters of the system yields mathematical results consistent with the expected physical ones.

Notizen: Abstract A generalization of Landahl's approximation method (H. D. Landahl,Bull. Math. Biophysics,15, 49–61, 1953) for non-linear diffusion problems is suggested. The method is applied to sorption, desorption, and free diffusion problems involving concentration-dependent diffusion coefficients. With some limitations, the results compare favorably with those obtained by numerical methods.

Notizen: Abstract Recently a theorem for representing current generators in a volume conductor by the superposition of a central dipole, quadrupole, octopole, etc., has been established by G. C. K. Yeh, J. Martinek, and H. de Beaumont (Bull. Math. Biophysics,20, 203–16, 1958). This theorem makes possible the representation of any discrete or line, surface- or volume-distributed current source by a unique model which can be determined for each given case by surface potential measurements and closed form analysis. In this paper the multipole representations of an eccentric dipole and an eccentric double-layer are obtained in terms of the various parameters of the assumed singularities, and the contributions to surface potentials due to each of the multipoles are compared. Certain numerical results corresponding to those of E. Frank (Amer. Heart J.,46, 364–78, 1953) are carried out and compared. Furthermore, the multipole representation of a partially damaged double-layer is also determined and compared with that of an undamaged one. It is concluded that within the range of parameters corresponding to human subjects the higher-order multipoles can contribute significantly to the surface potentials compared with the dipole.

Notizen: Abstract In line with a recent suggestion by the author (Bull. Math. Biophysics,20, 267–73, September, 1958) that not only does the organism as a whole map on the primordial, but that each organ can also be thus mapped, it is shown that the previously introduced abstract spaces, which represent an organism, contain subspaces which map continuously on the space of the primordial. Several theorems about those subspaces are proven.

Notizen: Abstract The DNA-protein coding problem is given a general algebraic formulation, the consequences of which are then explored by standard mathematical methods. To keep the treatment self-contained, the mathematical techniques to be used are explained in detail. It is demonstrated that there exista priori a countably infinite number of different abstract DNA-protein codes, thereby showing that inductive attempts to construct such a code will most likely be fruitless. A notion of ergodicity is then introduced, which imposes a number of restrictions on the admissible codes, and, in fact, these considerations enable us toderive a small portion of a code which, if our hypothesis of ergodicity is correct, must occur in nature. Finally, we discuss briefly the problem as to whether there can exist more than one DNA-protein code in nature.

Notizen: Abstract The present-day practices of electrocardiography and vectorardiography are based upon the theory that the surface potential differences can be assumed to be due to a single dipole inside the body. It is shown in this paper that a dipole cannot account for all the surface potentials due to realistic current generators, and hence the determination of the current generator from surface potential measurements based upon such a theory will lead to inconsistent representations of the heart for one and the same subject. To demonstrate this point two eccentric dipoles of different strengths and locations representing two muscle fibers are taken to be the current generator in a homogeneous spherical conductor. The exact surface potentials are then expressed by means of the “interior sphere theorem” of the authors. With these expressions the magnitude, direction, and location of the resultant dipole are determined by the method of D. Gabor and C. V. Nelson (J. App. Physics,25, 413–16, 1954). The surface potentials due to this resultant dipole are again exactly expressed by means of the “interior sphere theorem” and compared with those due to the eccentric dipoles assumed. It can be seen that the differences can be considerable. It is suggested that the multipole model of the authors (Bull. Math. Biophysics,20, 203–16, 1958) be used as a more accurate and the only unique representation of the heart.

Notizen: Abstract In a preceding paper (Bull. Math. Biophysics 20, 71–93, 1958) the principle of biotopological mapping was formulated in terms of a continuous mapping of an abstract space, made from the set of biological properties which characterize the organism, by an appropriate definition of neighborhoods. In this paper it is shown that we may consider directly the mappings of the different sets of properties which characterize different organisms without taking recourse to abstract spaces. All the verificable conclusions made in the preceding paper remain valid. A serious difficulty mentioned previously is, however, avoided and the possibility of more general predictions is established.

Notizen: Abstract The general Theory of Categories is applied to the study of the (M, R)-systems previously defined. A set of axioms is provided which characterize “abstract (M, R)-systems”, defined in terms of the Theory of Categories. It is shown that the replication of the repair components of these systems may be accounted for in a natural way within this framework, thereby obviating the need for anad hoc postulation of a replication mechanism. A time-lag structure is introduced into these abstract (M, R)-systems. In order to apply this structure to a discussion of the “morphology” of these systems, it is necessary to make certain assumptions which relate the morphology to the time lags. By so doing, a system of abstract biology is in effect constructed. In particular, a formulation of a general Principle of Optimal Design is proposed for these systems. It is shown under what conditions the repair mechanism of the system will be localized into a spherical region, suggestive of the nuclear arrangements in cells. The possibility of placing an abstract (M, R)-system into optimal form in more than one way is then investigated, and a necessary and sufficient condition for this occurrence is obtained. Some further implications of the above assumptions are then discussed.

Notizen: Abstract The transient stage of the random dispersal of logistic populations is investigated, using a Sturm-Liouville series leading to an infinite system of non-linear integral equations. These equations are then solved via a successive approximation scheme. R. A. Fisher's (steady-state) velocity of advance paradox is discussed. An illustrative example is worked to the second order of approximation.

Notizen: Abstract An approximation method using a sine function is used to solve the second degree growth equation for the case in which an organism may simultaneously become dispersed throughout a uniform region. The resulting expression for a special case is compared with the expression obtained by R. Barakat (1959,Bull. Math. Biophysics,21, 141–51), giving the first two terms, by an iterative, procedure. The agreement is satisfactory.

Notizen: Abstract Diffusion through a flat pore into a large open region is proportional to the linear dimension of the pore and not to its area. This was first explained by Brown and Escombe (1900) for a circular pore and is here generalized, by means of a dimensional argument, to include any type of regular opening. The problem is further generalized to include diffusion through pores of finite thickness, finite distance apart, and into finite regions. Since this problem cannot be solved exactly, an approximation method is introduced. Reasons for the credibility of the approximation are presented. It is then shown, by means of the approximation method, that the diffusive flow through a pore is equal to the total concentration difference divided by the resistance of the system. The resistance, in turn, is the sum of the resistances of all portions of the system, each of which is calculated. The result is compared with results which have been calculated exactly for limiting cases and found to agree very well. The results are then applied to a standard method of computing pore size in membranes, and it is shown that the correction factor is negligible.

Notizen: Abstract V. S. Ivlev [Experimental Ecology of Nutrition of Fishes, 1955, Moscow (in Russia)] has shown that the food uptake by fishes during a fixed interval of time is an exponential function of the concentration of food. Ivlev's equation is derived here, and it is shown that it can hold only for non-stationary conditions, such as prevailed in Ivlev's experiments. For a stationary state, the rate of food uptake should tend asymptotically to a limiting value as the concentration increases, but the variation is not exponential. Different other aspects of the problem are investigated, and definite new experimental procedures suggested. The implications of Ivlev's findings on the effect of non-uniformity of food distribution upon the rate of food consumption are studied from a mathematical point of view. The conclusion is reached that whereas a fish does not, in the process of eating, move directly to an individual food particle which it perceives, it does move more or less directly to large aggregates of particles, if the latter are distributed nonuniformly.

Notizen: Abstract Some relational aspects of the property of self-reproduction of biological systems are studied. If in addition to the requirement of the property of self-reproduction we add also the requirement of adaptability of the organism to changing environment, this imposes certain conditions on the topology of the graphs which represent such systems. A further study of the relational properties of such systems seems to offer the possibility of deriving the principle of biological mapping from the requirement of self-reproduction and adaptability. An examination of the problem of the original formation of such self-reproducing systems in connection with the established fact of impossibility of spontaneous generation leads to the conclusion that an organism must inhibit such processes which, in the absence of organisms, would lead to spontaneous generation.

Notizen: Abstract In the human, the antagonistic, extensor-flexor system of the leg is an example of a common type of neurophysiological feedback system. After a brief introduction to the neuroanatomy and physiology of this feedback system, the paper formulates transfer functions from temporal response data available in the literature. A feedback stability analysis, based on the extension of Nyquist's stability criteria to multiple-loop systems and utilizing flow-graph techniques, demonstrates the stable behavior of the system. Expressions are given relating the sensitivity of the system to variations in muscle response and Golgi tendon organ (tension receptor) response. By considering the events taking place at synapses and end-plates during “isometric tension-small knee angle excursion” conditions as stationary stochastic processes, an external “noise” input to the system is given, whose spectrum is derived from the statistics of a shot-process representation of these events. The paper concludes with some correlations between the analytical results and clinical syndromes.

Notizen: Abstract In this continuation of a previous report it is shown how the Volterra population dynamics, which underlies the statistical theory, can be based on a variational principle; how the dynamics can be generalized as regards both the behavior of total populations and migration phenomena; and how many directly observable data, such as amplitudes and frequencies of oscillation of a population, fit into the statistical theory and can test it. Such a test is carried out in some detail using the fox-catch data of Elton, with a clear indication that the theory is capable of comprehending the major statistical properties of population-time curves. A final section sketches an extension of the theory to cover secular variations of external conditions such as temperature of the environment.

Notizen: Abstract It is shown that from the definition of organismic sets (Rashevsky,Organismic Sets. Some Reflections on the Nature of Life and Society, Holland, Michigan, Mathematical Biology, Inc. and Grosse Pointe, Michigan, J. M. Richards Laboratory) a complete sensory deprivation of an organismic set of ordern=2 should result in malfunctioning of the set. A generalization to higher order sets is suggested.

Notizen: Abstract Optimality in branching structure of the vascular tree was studied. Analysis on its physiological roles as the duct system for blood supply to the capillaries predicted that the vascular tree should be constructed with minimum volume under restriction of determinant pressure, flow and location at the origin and the terminals. Mathematical derivations of this conditional extremum problem yielded some equations expressing the relations between the radii of the branches and their branching angles, which provided numerical solutions for branching points of bi- and poli-terminal minimum volume trees. Comparison of the peritoneal vascular tree in a dog with the minimum volume one computed under the same restrictive conditions showed good agreement in their branching structure.

Notizen: Abstract In order to determine the kinetics of passage of a substance through an organ containing a tangle of vessels, we study the response of a tube to various inflows (perfusion, brief injection, ...). The introduction of the catabolic terms and of the spatial dependence between bulk concentration and surface concentration allows one to account for the difference of arteriovenous concentrations observed experimentally for many metabolites. The relationships between the physico-chemical parameters of the organ and the operational parameters of the model demonstrate the importance of the transit time through the considered vessels. If one considers the different pathways as independent, the introduction of the transit time distribution for an inert substance enables one to compute the response of the organ analytically or by recurrence, using convolution. The parameters of the model can be obtained by the moments method.

Notizen: Abstract The solution of the diffusion equation in the gas phase of the human lung is very difficult because of the structure of the bronchial tree. It is shown by means of physical arguments, how one can reduce the diffusion equation to a simple one-dimensional form. The solution is then obtained by a stochastic simulation, which is easily realized on a digital computer.

Notizen: Abstract A new mathematical model of the oscillatory behavior of the respiratory center has been developed based upon published records of neuronal activity during respiration in the pons and medulla. In contrast with a previous model, four, rather than two, networks are assumed to interact in the respiratory center so as to produce the respiratory oscillation. A mathematical description of this interaction, in the form of a set of four first-order, nonlinear, coupled differential equations, is derived; the behavior of the solutions of this system is studied qualitatively, and expressions for the durations of the inspiratory and expiratory phases are obtained in terms of some parameters. It is found that central and chemical influences drive the medullar neurons to a position somewhere between saturation and full cutoff, and the pontine neurons deeply into cutoff. The control of the duration of the different phases by these chemical and central means is discussed. In order to effect a decrease in the magnitude of the various times, the neurons have to be driven towards operating points of higher central facilitation.

Notizen: Abstract In vivo control of calcium is analysed under the assumption that hormonal influences via plasma levels of parathormone and calcitonin are of prime (but not absolutely dominating) importance. A brief review concerning the physiological significance of body calcium and the mode of action of these two hormones is presented as an introduction to the basic philosophy of the study. A theoretical quasi-linear lumped-parameter model is developed to describe variations in ionic calcium, parathormone and calcitonin plasma concentrations to specific input stimuli. Formal evaluation of the system response requires the determination of ten constants, together with quantitation of ingested calcium entry into the plasma compartment which isindependent of hormonal influences. Values for various parameters are deduced from published data and experimental procedures are outlined to facilitate determination of the remaining unknowns. It is suggested that the proposed model should prove useful for investigations concerning general hormonal actions on calcium homeostatic mechanisms in both normal and diseased states, with particular reference to calcitonin.

Notizen: Abstract Ann species predator-prey chain is analyzed to determine what oscillations occur in population sizes. It is found that only the populations of the first and second species in the chain must necessarily oscillate around the point of equilibrium if they do not come to equilibrium. The other species may or may not oscillate.

Notizen: Abstract Sangren and Sheppard developed a mathematical model for first-order processes taking place in the regional circulation, applicable—for example—to tracer studies of potassium transport. It permits calculation of specific activity at any point along a “tube of flow” or in the cuff of tissue surrounding it as a function of time following a spike injection of tracer. In efforts to relate to the exchange a rate curves obtained within vivo counters pointed at the region of interest, we developed a compartment-system model of the process. In investigating the properties of the Sangren and Sheppard model integrated over an entire circulatory bed, as thein vivo counter would see it, we found that when the distribution of transit times of the “tubes of flow” can be approximated by an exponential sum, the solution reduces to that of the compartment system model. This results in an important simplification in the calculation, and insight into the assumptions underlying the two different models. A curve-fitting computer program for the compartment model has been written and applied to double-isotope studies of potassium transport in the hind leg of the dog.

Notizen: Abstract Two chemicals,A andB, are allowed to diffuse together and a reaction described by $$A + B\mathop \rightleftharpoons \limits_{K_{ - 1} }^{K_1 } C$$ is allowed to proceed. This system is described mathematically by a system of partial differential equations. A numerical procedure is presented to find the rate constants ofK 1 andK −1. A systematic analysis of the effects of errors is also presented.

Notizen: Abstract Equations are developed to describe the energy expenditure of the human heart. As well as the external potential and kinetic energy terms, general consideration is given to other possible avenues of energy consumption. Emphasis is placed upon using mathematical variables which are readily available for experimental verification. The errors involved in assuming that mean values for the physiological parameters give reasonable estimations for the external mechanical performance are examined, and a theoretical estimation for the discrepancy in the kinetic component is presented. Logical extension of the mathematical derivation leads to a determination of cardiac external mechanical efficiency and clearly demonstrates the significance of the ventricular pressure-volume loop in this context. Finally, experimental procedures are suggested to clarify further some of the conclusions reached through the theoretical analysis.

Notizen: Abstract The question is discussed as to the reason why some animal societies, such as bees or ants, are sexually differentiated, that is, onlysome of its members are exhibiting reproducing activities. It is indicated that human society may be on its way to such a sexual differentiation which may eventually come.

Notizen: Abstract X-ray diffraction patterns obtained experimentally for fibers, together with their chemical structures, can be analyzed theoretically in terms of an integral equation. The partially unknown electron density function can be solved by iteration. This mathematical technique has been applied with success to study the secondary structures of DNA fibers.

Notizen: Abstract For an environmental system described by a system of nonlinear first-order differential equations, the problem of achieving specified terminal conditions in a given time with a minimum expenditure of resources is considered. The initial conditions and the minimum value are found numerically in a particular example.

Notizen: Abstract A kinetic model of neural systems is introduced and discussed with statistical mechanics techniques. It is assumed that, for a macroscopic description of the model, it suffices to consider only the distribution for the velocity and position of the impulses, and the distribution for the excitation and position of the neurons, at any timet. Making use of Boltzmann's method for the study of a dilute gas, coupled differential equations for the rate of change with time of the distributions have been constructed.

Notizen: Abstract Creeping flow of a Newtonian fluid through a rigid permeable tube is considered and the transmural seepage is assumed to obey Darcy's law. Closed-form solutions for the pressure and velocity fields are presented and equations describing the axial variation of the mean cross-sectional pressure, the axial volumetric flow and the transmural fluid flux are derived. Approximate solutions for small seepage rates are given and are applied to the flow in the proximal renal tubule. Probable values for the epithelium permeability and the intraluminal hydrostatic pressure drop are obtained.

Notizen: Abstract The paper demonstrates that it is possible to construct memory models where the information inserted is stored in disseminated form, using sequential coding, the changes in the units forming the models being determined by their geometrical connections and by the incoming stream of information. The models are shown to have large storage capacity and their efficiency can be made insensitive to loss of or damage to a large fraction of their units. The satisfactory verification by computer simulation of the analysis and results described in the present paper will be the subject of a future paper.

Notizen: Abstract We examine in detail Edward Kerner’s method for linearizing the equations of enzyme kinetics. Our main result is the determination of canonical forms for systems which can be linearized by the method. This is done both in general and in the special cases of two and three dimensions where complete results are obtained. The practical problem of identifying linearizable systems is also considered and computable necessary criteria are presented.

Notizen: Abstract A general method for determination of the volume of a space in a non steady state condition, in case diffusion might be significant, is developed. Instantaneous mixing of indicators with native fluid is assumed in this first stage of investigation. Theoretical expressions are obtained for the volume of the space and the diffusion coefficient as a function of time. An analysis of feasibility of the method is also included.

Notizen: Abstract The theory of transfer of low-molecular nonelectrolytes across deformable semipermeable membranes of large curvature developed in Part I (Rubinstein, 1974) is used to describe the dynamics of swelling and shrinking of a muscle fiber at the influx and efflux of low-molecular nonelectrolytes. A large set of computations showed that the theory explains the experiments described in the literature.

Notizen: Abstract Green's function for heat and matter transport is calculated for an infinite medium in which a convection field v(r,t) makes a contribution to the total heat and matter current. It is given by a uniformly and absolutely convergent series in which every term is calculated from the preceding one merely by integration. The solution procedure is interpreted physically and illustrated by a simple problem in which v(r,t)=const. in space and time. Since the solution contains no intrinsic spatial symmetry, it can serve as a starting point for a theory of heat and mass transport in perfused biological tissue.

Notizen: Abstract The hydrodynamics of a microorganism swimming in a channel is investigated. The microorganism is modeled as a two-dimensional sheet swimming at low Reynolds numbers between two rigid walls. The wavelengths of the propulsive waves passing down the sheet are assummed to be very large compared to the channel spacing, but the amplitude of the propulsive waves is arbitrary. Explicit analytical solutions for the propulsive velocity and the rate of energy dissipated in terms of the wave amplitude, channel spacing, wave number, and wave speeds are given.

Notizen: Abstract It is shown from the statistical-mechanical overview of Volterra's ecological model how to reckon the fluctuations of collective variables such as the total population of a genus: and that these fluctuations are much decreased (or that the collective populationsteadiness is enhanced) as the speciation is increased. (A niching of species in time, or phase-niching, is entailed here.) Secondly, it is shown how Preston's log-normal distribution describing the species-abundance relationship, as well as a generalization of such distributions, come forth simply and naturally from the statistical-Volterra-dynamics.

Notizen: Abstract In the framework of the neural network theory effects similar to hypnotic displays are constructed. They are based on the associative paradigm involving non-linear interaction of excitatory and inhibitory channels with synaptic memory. The non-linearity of long-term memorizing processes may cause effects exhibited by blind spots, which are interpreted as the first stage of hypnosis. More complicated phenomena are discussed in terms of a two-layer network.

Notizen: Abstract Mutation is introduced into autocatalytic reaction networks. The differential equations obtained are neither of repliator-type nor can they be transformed straightway into a linear equation. Examples of low dimensional dynamical systems —n=2, 3 and 4 — are discussed and complete qualitative analysis is presented. Error thresholds known from simple replication-mutation kinetics with frequency independent replication rates occur here as well. Instead of cooperative transitions or higher order phase transitions the thresholds appear here as supercritical or subcritical bifurcations being analogous to first-order phase transitions.

Notizen: Abstract The non-linear behavior of a differential equations-based predator-prey model, incorporating a spatial refuge protecting a consant proportion of prey and with temperature-dependent parameters chosen appropriately for a mite interaction on fruit trees, is examined using the numerical bifurcation code AUTO 86. The most significant result of this analysis is the existence of a temperature interval in which increasing the amount of refuge dynamically destabilizes the system; and on part of this interval the interaction is less likely to persist in that predator and prey minimum population densities are lower than when no refuge is available. It is also shown that increasing the amount of refuge can lead to population outbreaks due to the presence of multiple stable states. The ecological implications of a refuge are discussed with respect to the biological control of mite pests.

Notizen: Abstract In many applications of control theory on plant growth models biomass maximization is postulated to avoid analytically unsolvable problems while fruit maximization is commonly considered to be a more realistic criterion. In a special case, we are able to compare these criteria. Iwasa and Roughgarden (1984,Theor. Pop. Biol. 25, 78–105) have investigated a certain class of plant growth models using a fruit maximization criterion. They proved that, in the vegetative growth period, the organs follow a certain path of balanced growth. We show that this path remains optimal when biomass maximization is postulated. This underlines the importance of the balanced growth path found by Iwasa and Roughgarden. Furthermore, our result suggests that in the vegetative growth period the biomass maximization criterion is a good approximation of fruit maximization. In another theoretical control investigation, Schultzeet al. (1983,Oecologia 58, 169–177) derived a different type of balanced growth path. We apply the theory of Iwasa and Roughgarden to an improved version of the model of Schulzeet al. This leads to a new description of balanced growth between root and shoot that reflects non-linearities in the water uptake process and constitutes an interesting hypothesis for further experimental testing.

Notizen: Abstract In this paper the effects of changing the ion concentration in and around a sample of soft tissue are investigated. The triphasic theory developed by Laiet al. (1990,Biomechanics of Diarthrodial Joints, Vol. 1, Berlin, Springer-Verlag) is reduced to two coupled partial differential equations involving fluid ion concentration and tissue solid deformation. These equations are given in general form for Cartesian, cylindrical and spherical geometries. After solving the two equations quantities such as fluid velocity, fluid pressure, chemical potentials and chemical expansion stress may be easily calculated. In the Cartesian geometry comparison is made with the experimental and theoretical work of Myerset al. (1984,ASME J. biomech. Engng,106, 151–158). This dealt with changing the ion concentration of a salt shower on a strip of bovine articular cartilage. Results were obtained in both free swelling and isometric tension states, using an empirical formula to acount for ion induced deformation. The present theory predicts lower ion concentrations inside the tissue than this earlier work. A spherical sample of tissue subjected to a change in salt bath ion concentration is also considered. Numerical results are obtained for both hypertonic and hypotonic bathing solutions. Of particular interest is the finding that tissue may contract internally before reaching a final swollen equilibrium state or swell internally before finally contracting. By considering the relative magnitude, and also variation throughout the time course of terms in the governing equations, an even simpler system is deduced. As well as being linear the concentration equation in the new system is uncoupled. Results obtained from the linear system compare well with those from the spherical section. Thus, biological swelling situations may be modelled by a simple system of equations with the possibility, of approximate analytic solutions in certain cases.

Notizen: Abstract Many models of immune networks have been proposed since the original work of Jerne [1974,Ann. Immun. (Inst. Pasteur) 125C, 373–389]. Recently, a limited class of models (Weisbuchet al., 1990,J. theor. Biol. 146, 483–499) have been shown to maintain immunological memory by idiotypic network interactions. We examine generalizations of these models when the networks are both large and highly connected to study their memory capacity, i.e. their ability to account for immunization to a large number of random antigens. Our calculations show that in these minimal models, random connectivities with continuously distributed affinities reduce the memory capacity to essentially nil.

Notizen: Abstract A kinetic model is proposed to delineate the factors that determine the coronary reactive hyperemic response (RHR) to transient ischemia. The model comprises of myocardial-interstitial (M) and vascular (V) compartments. Vasodilator metabolites (VM) are produced in the M compartment during the interval of coronary occlusion. The rate of VM production is dependent on the flow rate during the ischemic period, the ratio of excess flow above the control level (R) to the loss of flow during occlusion period (D), the amount of oxygen stored and the degree of vasodilation in the V compartment prior to occlusion. Following a complete release of occlusion, VM are transported from the M to V compartment and are washed out or degraded with time. The time course of RHR is determined by the coronary patency which is proportional to VM concentration in the V compartment. Based on a set of numerical constants, the model is tested by simulating RHR to the various occlusion manoeuvres: a pair of 10 sec occlusions separated by brief release, a 15 sec release followed by a second brief occlusion, a brief release of an occlusion followed by restriced inflow and a period of restricted inflow after occlusion. The simulated results fit the experimental R/D and RH durations data of canine hearts. Factors that determine the impairment of RH capacity in coronary stenosis are suggested in terms of the model scheme.

Notizen: Abstract In the present paper a kinetic study is made of the behaviour of a Michaelis-Menten enzyme-catalysed reaction in the presence of irreversible inhibitors rendered unstable in the medium by their reaction with the product of enzymatic catalysis. A general mechanism involving competitive, non-competitive, uncompetitive and mixed irreversible inhibition with one or two steps has been analysed. The differential equation that describes the kinetics of the reaction is non-linear and computer simulations of its dynamic behaviour are presented. The results obtained show that the systems studied here present kinetic co-operativity for a target enzyme that follows the simple Michaelis-Menten mechanism in its action on the substrate, except in the case of an uncompetitive-type inhibitor.

Notizen: Abstract The relative contributions of mitochondrial β-oxidation and peroxisomal β-oxidation and peroxisomal ω-oxidation to the oxidation of a given fatty acidin vivo can be quantitated by an isotopic method. The approach requires infusion of a fatty acid labelled on two specific carbon atoms (e.g. [1-14C] and [11-14C] palmitate) to an isotopic steady state, with subsequent isolation and degradation of an acetylated conjugate as a product of the liver cytosolic acetyl CoA pool and of ketone bodies as a product of the liver mitochondrial acetyl CoA pool.

Notizen: Abstract Pancreatic β-cells in intact islets of Langerhans perfused with various glucose concentrations exhibit periodic bursting electrical activity (BEA) consisting of active and silent phases. The fraction of the time spent in the active phase is called the plateau fraction and appears to be strongly correlated with the rate of release of insulin from islets as glucose concentration is varied. Here this correlation is quantified and a theoretical development is presented in detail. Experimental rates of insulin release are correlated with “effective” plateau fractions over a range of glucose concentrations. There are a number of different models for BEA in pancreatic β-cells and a method is developed here to quantify the dependence of a glucose dependent parameter on glucose concentration. As an example, the plateau fractions computed from the Sherman-Rinzel-Keizer model are matched with experimental plateau fractions to obtain a relationship between the model's glucose-dependent parameter, β, and glucose concentration. Knowledge of the relationships between β and glucose concentration and between experimental measurements of rates of insulin release and plateau fractions permits the determination of theoretical rates of insulin release from the model.

Notizen: Abstract When a suspension of bacterial cells of the speciesBacillus subtilis is placed in a chamber with its upper surface open to the atmosphere complex bioconvection patterns are observed. These arise because the cells: (1) are denser than water; and (2) usually swim upwards, so that the density of an initially uniform suspension becomes greater at the top than the bottom. When the vertical density gradient becomes large enough, an overturning instability occurs which ultimately evolves into the observed patterns. The reason that the cells swim upwards is that they are aerotactic, i.e. they swim up gradients of oxygen, and they consume oxygen. These properties are incorporated in conservation equations for the cell (N) and oxygen (C) concentrations, and these are solved in the pre-instability phase of development whenN andC depend only on the vertical coordinate and time. Numerical results are obtained for both shallow- and deep-layer chambers, which are intrinsically different and require different mathematical and numerical treatments. It is found that, for both shallow and deep chambers, a thin boundary layer, densely packed with cells, forms near the surface. Beneath this layer the suspension becomes severely depleted of cells. Furthermore, in the deep chamber cases, a discontinuity in the cell concentration arises between this cell-depleted region and a cell-rich region further below, where no significant oxygen concentration gradients develop before the oxygen is fully consumed. The results obtained from the model are in good qualitative agreement with the experimental observations.

Notizen: Abstract We describe a classification scheme for bursting oscillations which encompasses many of those found in the literature on bursting in excitable media. This is an extension of the scheme of Rinzel (inMathematical Topics in Population Biology, Springer, Berlin, 1987), put in the context of a sequence of horizontal cuts through a two-parameter bifurcation diagram. We use this to describe the phenomenological character of different types of bursting, addressing the issue of how well the bursting can be characterized given the limited amount of information often available in experimental settings.

Notizen: Abstract To ensure its sustained growth, a tumour may secrete chemical compounds which cause neighbouring capillaries to form sprouts which then migrate towards it, furnishing the tumour with an increased supply of nutrients. In this paper a mathematical model is presented which describes the migration of capillary sprouts in response to a chemoattractant field set up by a tumour-released angiogenic factor, sometimes termed a tumour angiogenesis factor (TAF). The resulting model admits travelling wave solutions which correspond either to successful neovascularization of the tumour or failure of the tumour to secure a vascular network, and which exhibit many of the characteristic features of angiogenesis. For example, the increasing speed of the vascular front, and the evolution of an increasingly developed vascular network behind the leading capillary tip front (the brush-border effect) are both discernible from the numerical simulations. Through the development and analysis of a simplified caricature model, valuable insight is gained into how the balance between chemotaxis, tip proliferation and tip death affects the tumour's ability to induce a vascular response from neighbouring blood vessels. In particular, it is possible to define the success of angiogenesis in terms of known parameters, thereby providing a potential framework for assessing the viability of tumour neovascularization in terms of measurable quantities.

Notizen: Abstract The parameter domain for which the quasi-steady state assumption is valid can be considerably extended merely by a simple change of variable. This is demonstrated for a variety of biologically significant examples taken from enzyme kinetics, immunology and ecology.

Notizen: Abstract The ability of random fluctuations in selection to maintain genetic diversity is greatly increased when generations overlap. This result has been derived previously using genetic models with very special assumptions about the population age structure. Here we explore its robustness in more realistic population models, with very general age structure or physiological structure. For a range of genetic models (haploid, diploid, single and multilocus) we find that the condition for maintaining genetic diversity generalizes almost without change. Genetic diversity is maintained by selection if a product of the form (generation overlap)×(selection intensity)×(variability in the selection regime) is sufficiently large, where the generation overlap is measured in units of Fisher's reproductive value. This conclusion is based on a local evolutionary stability analysis, which differs from the standard “protected polymorphism” criterion for the maintenance of genetic diversity. Simulation results match the predictions from the local stability analysis, but not those from the protected polymorphism criterion. The condition obtained here for maintaining genetic diversity requires fitness fluctuations that are substantial but well within the range observed in many studies of natural populations.

Notizen: Abstract Premised on relatively simple assumptions, mathematical models like those of Monod, Pirt or Droop inadequately explain the complex transient behavior of microbial populations. In particular, these models fail to explain many aspects of the dynamics of aTetrahymena pyriformis-Escherichia coli community. In this study an alternative approach, an individual-based model, is employed to investigate the growth and interactions ofTetrahymena pyriformis andE. coli in a batch culture. Due to improved representation of physiological processes, the model provides a better agreement with experimental data of bacterial density and ciliate biomass than previous modeling studies. It predicts a much larger coexistence domain than rudimentary models, dependence of biomass dynamics on initial conditions (bacteria to ciliate biomasses ratio) and appropriate timing of minimal bacteria density. Moreover, it is found that accumulation ofE. coli sized particles andE. coli toxic metabolites has a stabilizing effect on the system.

Notizen: Abstract At the core of contemporarymorphometrics—the quantitative study of biological shape variation—is a synthesis of two originally divergent methodological styles. One contributory tradition is the multivariate analysis of covariance matrices originally developed as biometrics and now dominant across a broad expanse of applied statistics. This approach, couched solely in the linear geometry of covariance structures, ignores biomathematical aspects of the original measurements. The other tributary emphasizes the direct visualization of changes in biological form. However, making objective the biological meaning of the features seen in those diagrams was always problematical; also, the representation of variation, as distinct from pairwise difference, proved infeasible. To combine these two variants of biomathematical modeling into a valid praxis for quantitative studies of biological shape was a goal earnestly sought though most of this century. That goal was finally achieved in the 1980s when techniques from mathematical statistics, multivariate biometrics, non-Euclidean geometry and computer graphics were combined in a coherent new system of tools for the complete regionalized quantitative analysis oflandmark points together with the biomedical images in which they are seen. In this morphometric synthesis, correspondence of landmarks (biologically labeled geometric points, like “bridge of the nose”) across specimens is taken as a biomathematical primitive. The shapes of configurations of landmarks are defined as equivalence classes with respect to the Euclidean similarity group and then represented as single points in David Kendall'sshape space, a Riemannian manifold with Procrustes distance as metric. All conventional multivariate strategies carry over to the study of shape variation and covariation when shapes are interpreted in the tangent space to the shape manifold at an average shape. For biomathematical interpretation of such analyses, one needs a basis for the tangent space compatible with the reality of local biotheoretical processes and explanations at many different geometric scales, and one needs graphics for visualizing average shape differences and other statistical contrasts there. Both of these needs are managed by thethin-plate spline, a deformation function that has an unusually helpful linear algebra. The spline also links the biometrics of landmarks to deformation analysis of the images from which the landmarks originally arose. This article reviews the history and principal tools of this synthesis in their biomathematical and biometrical context and demonstrates their usefulness in a study of focal neuroanatomical anomalies in schizophrenia.

Notizen: Abstract A competition model describing tumor-normal cell interaction with the added effects of periodically pulsed chemotherapy is discussed. The model describes parameter conditions needed to prevent relapse following attempts to remove the tumor or tumor metastasis. The effects of resistant tumor subpopulations are also investigated and recurrence prevention strategies are explored.

Notizen: Abstract Increasing attention is being paid to the configuration and development of vascular structures and their possible correlations with physiological events. The study of angiogenesis in normal and pathological states as well as in the embryo and adult has provided new insights into the mechanism of vessel growth and organization of the vasculature. Various mathematical branching models have been developed. These constructions are mainly geometrical and only involve a branching phenomenon. We propose the use of a deterministic non-linear model based on physiological laws and hydrodynamics. Growth, branching and anastomosis, the three actual main events occurring in vascular growth, are included in this model. Space growth, including cells and vessels, is defined by a decreasing transformation. Space density and the length of new sprouts are controlled by a set of parameters. The conditions on these parameters are well established, which allows the production of realistic patterns.

Notizen: Abstract The quasi-stationary distribution of a population within a system of interacting populations is approximated by a stochastic logistic process. The parameters of this process can be expressed in the parameters of the full system. Using the diffusion approximation, an expression for the expected extinction time is derived from this logistic process. Since the expected extinction time is expressed in the parameters of the full system, the effect of these parameters on the extinction risk can be easily evaluated, which may be of use for studies in ecology, conservation biology and epidemiology. The outcome is compared with simulation results for the case of a prey-predator system.

Notizen: Abstract The cytokines are the information superhighway of the immune system. They are an important component of the integrated behavior of the system. In order to be able to have a good understanding of the immune system, we must be able to model the effect of cytokines and their combined effect. This work is a step in that direction. We study the combined effect of two cytokines: interleukin-2 (IL-2) and interleukin-4 (IL-4) on some cells of the immune system. Interleukin-2 and interleukin-4 are important growth and differentiation factors for B and T cells. Interleukin-4 antagonizes the effect of interleukin-2 on B cells and some T cells while it synergizes with interleukin-2 on other T cells. We build a mathematical model of the interaction of both cytokines on T and B cells as a building block toward a model of the Th1/Th2 cross-regulation. The response of a given cell to the combination of interleukin-2 and interleukin-4 is shown to involve competing dynamical effects which can lead to either antagnostic or synergistic combined effect.

Notizen: Abstract We propose a mathematical approach to the modelling of self-organizing hierarchies in animal societies. This approach relies on a basic positive feedback mechanism that reinforces the ability of a given individual to win or to lose in a hierarchical interaction, depending on how many times it won or lost in previous interactions. Motivated by experiments carried out on primitively eusocial waspsPolistes, the model, is based on coupled differential equations supplemented with a small stochastic term. Numerical integrations allow many different hierarchical profiles to be obtained depending on the model parameters: (1) the particular form of the probability for an individual to win or lose a fight given its history, (2) the probability of interaction between two individuals, (3) the forgetting strength, which determines the rate at which events in the past are forgotten and no longer influence the force of an individual and (4) two individual recognition parameters, which set the contribution of individual recognition in the process of hierarchical genesis. We compare the results, expressed in terms of a hierarchical index or of the Landau number that describes the degree of linearity of the hierarchy, with various experimental results.

Notizen: Abstract The normal process of dermal wound healing fails in some cases, due to fibro-proliferative disorders such as keloid and hypertrophic scars. These types of abnormal healing may be regarded as pathologically excessive responses to wounding in terms of fibroblastic cell profiles and their inflammatory growth-factor mediators. Biologically, these conditions are poorly understood and current medical treatments are thus unreliable. In this paper, the authors apply an existing deterministic mathematical model for fibroplasia and wound contraction in adult mammalian dermis (Olsenet al., J. theor. Biol. 177, 113–128, 1995) to investigate key clinical problems concerning these healing disorders. A caricature model is proposed which retains the fundamental cellular and chemical components of the full model, in order to analyse the spatiotemporal dynamics of the initiation, progression, cessation and regression of fibro-contractive diseases in relation to normal healing. This model accounts for fibroblastic cell migration, proliferation and death and growth-factor diffusion, production by cells and tissue removal/decay. Explicit results are obtained in terms of the model processes and parameters. The rate of cellular production of the chemical is shown to be critical to the development of a stable pathological state. Further, cessation and/or regression of the disease depend on appropriate spatiotemporally varying forms for this production rate, which can be understood in terms of the bistability of the normal dermal and pathological steady states—a central property of the model, which is evident from stability and bifurcation analyses. The work predicts novel, biologically realistic and testable pathogenic and control mechanisms, the understanding of which will lead toward more effective strategies for clinical therapy of fibro-proliferative disorders.

Notizen: Abstract Populations often exhibit abrupt changes in abundance associated with a smooth, continuous change in some component of their environment, with the abruptness usually attributed to inter-specific interactions or physical extremes. This paper presents a spatially explicit single-species population model in which intra-specific interactions alone are responsible for such an abrupt change. The essential mechanism involves cooperation in both colonization (through enhanced recruitment near other individuals) and mortality (protection through a “safety-in-numbers” interaction). Large fluctuations in population density would likely be observable near the transition region.

Notizen: Abstract A mathematical model of the nitrogen transformation cycle in an aquatic environment is studied. Using Pontryagin's maximum principle, a preferential utilization of ammonium to nitrate by phytoplankton is explained and verified by experimental data. A multiparameter bifurcation is given. The model was found to have four types of equilibrium sets. It is shown that a Hopf bifurcation may occur.

Notizen: Abstract Parallel computation employing a domain decomposition method was used to calculate precisely without approximations the spatio-temporal distribution of Ca2+ in nerve terminals. The results showed, contrary to expectations, that for equal admitted Ca2+ currents at low (one channel open) and high (four channels open) depolarization, the average Ca2+ concentration at the release area is higher at the low depolarization. These calculations provide additional support for the Ca2+-voltage hypothesis for neurotransmitter release.

Notizen: Abstract Type I hypersensitivity, which functions to protect the organism from parasites, is caused by binding of antigen to IgE antibodies pre-attached to the cell surface of tissue mast cells and their circulating counterparts, the basophils. In “allergy,” type I hypersensitivity is inappropriately induced by protein-based foreign substances (such as pollen) or protein components of insect stings, which in the normal course of events would be cleared from the organism without causing any damage. Paradoxically, a successful clinical treatment of allergy involves repeated immunization of allergic persons with low doses of the allergen—immunotherapy. Investigation of the available experimental evidence leads to the conclusion that the phenomena of immunotherapy are best addressed in terms of the interplay among the mechanism(s) of immune memory—Th1/Th2 cross-regulation—and the physical compart-mentalization of the immune system. These conclusions are illustrated with a numerical simulation.

Notizen: Abstract We consider a host-solitary parasitoid system with three categories of individuals: parasitoids, healthy hosts and parasitized hosts. Parasitoids are assumed to discriminate perfectly between the two kinds of hosts and they can reject those which are already parasitized. If parasitoids systematically accept or reject superparasitism or behave randomly, the system is always unstable. Using an optimal foraging model, we determine the behavior of parasitoids which leads to maximization of the instantaneous reproductive rate. When following this adaptive decision rule, parasitoids accept or refuse superparasitism according to the densities of both healthy and parasitized hosts. We study the dynamics of the system when parasitoids follow the optimal rule and show that under certain conditions it possesses a locally stable equilibrium point. In addition, our model predicts that at equilibrium parasitoids show partial preferences for superparasitism.

Notizen: Abstract A system of differential equations for the control of tumor growth cells in a cycle nonspecific chemotherapy is analyzed. Spontaneously acquired drug resistance is taken into account, and a criterion for the selection of chemotherapeutic treatment is used. This criterion purports to describe the possibility of improvement of the patient's health when treatment is discontinued. Contrary to our early results which also take drug resistance into account, in this context strategies of continuous chemotherapy in which rest periods take part may be better than maximum drug concentration throughout the treatment (which appears to be in accordance with clinical practice). This bears out our previous conjecture that when drug resistance is accounted for, the imperfections in the usual modelling of treatment criteria, which in general do not allow for patient recuperation, ruled out the possibility of rest periods in optimal continuous chemotherapy.

Notizen: Abstract A logistic density-dependent matrix model is developed in which the matrices contain only parameters and recruitment is a function of adult population density. The model was applied to simulate introductions of white-tailed deer into an area; the fitted model predicted a carrying capacity of 215 deer, which was close to the observed carrying capacity of 220 deer. The rate of population increase depends on the dominant eigenvalue of the Leslie matrix, and the age structure of the simulated population approaches a stable age distribution at the carrying capacity, which was similar to that generated by the Leslie matrix. The logistic equation has been applied to study many phenomena, and the matrix model can be applied to these same processes. For example, random variation can be added to life history parameters, and population abundances generated with random effects on fecundity show both the affect of annual variation in fecundity and a longer-term pattern resulting from the age structure.

Notizen: Abstract In the presence of seasonal forcing, predator-prey models with quadratic interaction terms and weak dissipation can exhibit infinite numbers of coexisting periodic attractors corresponding to cycles of different magnitude and frequency. These motions are best understood with reference to the conservative case, for which the degree of dissipation is, by definition, zero. Here one observes the familiar mix of “regular” (neutrally stable orbits and tori) and chaotic motion typical of non-integrable Hamiltonian systems. Perturbing away from the conservative limit, the chaos becomes transitory. In addition, the invariant tori are destroyed and the neutrally stable periodic orbits becomes stable limit cycles, the basins of attraction of which are intertwined in a complicated fashion. As a result, stochastic perturbations can bounce the system from one basin to another with consequent changes in system behavior. Biologically, weak dissipation corresponds to the case in which predators are able to regulate the density of their prey well below carrying capacity.

Notizen: Abstract The standard method for measuringin vitro antibiotic efficacy is based on a point observation of bacterial activity 18 hours after inoculation. The method, while simple, forgoes significant information by ignoring the dynamics of the interations between antibiotic and bacteria. This paper proposes a simple dynamic model describing these interactions. The model consists of two non-linear differential equations of the S-system type. Its parameter values are estimated, through the minimization of residual errors, from data on the effect of the carbapenem antibiotic imipenem onPseudomonas aeruginosa. The model adequately describes the dynamic behavior of the bacterial populations in the presence of the antibiotic: beginning with drug administration, then through the decline of the bacterial population and possibly ending with bacterial resurgence.