ALBERT

Notes: Abstract The entropy budget of a white-tailed deer (50kg) on a maintenance diet and a full-feed diet in a standing posture in an open field under clear nocturnal skies with an air temperature of −20°C is investigated based on the energetics given by Moen. Entropy inflow into a white-tailed deer due to infra-red radiation and entropy outflows from a deer due to infra-red radiation, convection, evaporation of water and conduction to ingested food are calculated. Also the entropy production due to metabolic heat production is estimated. Net entropy flow into a deer from its environment becomes negative. On the assumption that a white-tailed deer is in a steady state in entropy, the total entropy production in a deer on a maintenance diet becomes +0.46 J/sec/K. Positiveness of the entropy production shows that the Second Law of Thermodynamics certainly holds in a white-tailed deer. The entropy production per effective radiating surface area of a deer on a maintenance diet is 0.32×10−4 J/cm2/sec/K. On the other hand, the entropy production in a deer on a full-feed diet is 0.59 J/sec/K and that per effective surface area is 0.41×10−4 J/cm2/sec/K. Uptake of 1 g of food produces 22 J/K of entropy within the body of a white-tailed deer. Comparison is made with the results for entropy production in a lizard and in plant leaves.

Notes: Abstract Biological adaptability has been proved to be analysable by means of the Maximum Entropy Formalism (MAXENT) in some cases of non-interacting systems. This formalism is extended to the biomass statistical structures of populations exhibiting internal interactions (i.e. predatorprey effects).

Notes: Abstract The temporal behaviours of the nonlinear substructure of a self-organized compartmental model of calcium metabolism were investigated. The order-two autocatalytic process included in this simple two-dimensional model is compared to some secondary nucleation mechanisms which should take place at the extracellular fluid-bone interface. The model gives rise to complex dynamic behaviours, and multistability properties, involving up to two stable periodic regimes (birhythmicity), were established in different topological configurations. The bifurcations occurring on the boundaries between regions of different qualitative behaviour have been determined. These properties are discussed in relation to the dynamical behaviour of other two-variable models, especially those including the same nonlinearity.

Notes: Abstract A linearized oscillation theorem due to Kulenović, Ladas and Meimaridou (1987,Quart. appl. Math. XLV, 155–164) and an extension of it are applied to obtain the oscillation of solutions of several equations which have appeared in population dynamics. They include the logistic equation with several delays, Nicholson's blowflies model as described by Gurney, Blythe and Nisbet (1980,Nature, Lond. 287, 17–21) and the Lasota-Wazewska model of the red blood cell supply in an animal. We also developed a linearized oscillation result for difference equations and applied it to several equations taken from the biological literature.

Notes: Abstract A theoretical approach to the explanation of the structural design of metabolic pathway is presented. It is based on the hypothesis that due to natural selection during evolution the cellular metabolism of present-day organisms may be characterized by optimal properties. Two cardinal terms enter the theory: (i) the efficiency of a metabolic pathway and (ii) the evolutionary effort for the change of the kinetic parameters of enzymes by mutations of the corresponding genes. For both quantities simple mathematical expressions are proposed. While the efficiency is related to the reaction rates of the enzymes constituting the metabolic pathway, the evolutionary effort is considered to be a monotonically increasing function of the parameter values. By maximizing the efficiency under the constraint of a fixed evolutionary effort the theory allows the calculation of the optimal parameter distribution as the outcome of evolution processes. The methods developed are applied to the following systems: (a) linear reaction sequences with very low affinities of the enzymes towards substrates, (b) linear sequences consisting of saturable enzymatic reactions, (c) branched metabolic pathways consisting of segments of linear chains and (d) glycolysis of erythrocytes. The conclusion is derived that the optimal distribution of kinetic constants depends strongly on the equilibrium constants of the reactions as well as on the total osmolarity of the metabolic intermediates. Without osmotic constraints the evolutionary effort is mainly spent on the enzymes at the beginning of the chain. Using Michaelis-Menten equations the optimal state is characterized by a decrease of the maximal activities of the enzymes towards the end of the chain. These results are modified if osmotic constraints are taken into account. At the investigation of branched pathways the following results were obtained: firstly, if a certain end product may be synthesized along different pathways those which are thermodynamically more unfavourable (e.g. characterized by a small change of free energy) are eliminated in the course of evolution; secondly, if a branched pathway leads to several important end products those reaction segments which are thermodynamically unfavourable are characterized by a higher evolutionary effort. The application of the theory to a realistic model of glycolysis of erythrocytes leads to a correct description of various functionally important properties of the system, such as the ratio between fluxes through different branches and the ATP/ADP ratio, whereas the theory cannot predict the strong separation of time constants observed in the real glycolytic system. It is concluded that the improvement of the predictive power of the theory necessitates the use of more complex functionals for the efficiency which take into account not only the fluxes but also other system properties such as the stability of the pathway or homoeostatic effects.

Notes: Abstract In order to evaluate the effect of anatomic asymmetries on the gas concentration distribution in the pulmonary airways, a Monte Carlo simulation of combined bulk flow and molecular diffusion was carried out in a realistic distal airway model (Parkeret al., 1971). This airway model, composed of branches distal to the 0.5-ram diameter airways, contained an upper symmetric segment consisting of four generations of conducting airways and a lower asymmetric segment of alveolar ducts and sacs arranged in five transport paths of varying lengths. In accounting for the volume increases of these ducts and sacs occurring during normal respiration, uniform alveolar filling rates and a fixed length-to-diameter ratio of all airways were assumed. For a pulse injection of inert tracer gas, the simulation was employed to determine the longitudinal concentration profiles in the conducting airways. In the alveolated airways, not only were the longitudinal profiles determined along each path, but radial transport from the core to the periphery of the airways was considered. The results of the simulations indicate that geometric asymmetries alone contribute substantially to regional concentration variations in the distal airways. For example, when a gas bolus is injected at mid*inspiration, there are concentration differences as great as 40% between two points along different transport paths located equi-distant from the proximal end of the model. As viewed from the terminal end of the model (acinus), average concentration differences as large as 6-to-1 exist between the longest and shortest transport paths respectively for gas boli introduced near the end of inspiration. The results further indicate because of large radial diffusion rates, no significant concentration differences exist between the periphery a-ld the central core of alveolated airways. Simulation of the expired concentration profiles indicate that boll injected very late during inspiration exhibit a sloping tail, unlike the earlier injected boll whose tails are virtually horizontal. Through the use of superposition teehniqnes, it was found that these sloping tails correspond to an alveolar slope of 1.5 vol% between 750 and 1250 ml expired for a continuous washing of tracer. This result is in disagreement with other transport analyses which did not directly account for the effect of geometric asymmetries.

Notes: Abstract Assuming a spherical geometry for the left ventricle, passive elastic stiffness-stress relations have been obtained on the basis of linear elasticity theory and large deformation theory. Employing pressure-volume aata taken from rat hearts of various age groups, it is shown that young rat heart muscle (1 month) is stiffer than either adult (7 months) or old rat heart muscle (17 months). Although the qualitative results are similar for both elasticity theories, the large deformation theory gave results in closer agreement with those obtained from papillary muscle studies. These results imply that stiffness of muscleper se can be assessed from left ventricular pressure-volume data.

Notes: Abstract Deliberate evaluation of the quantum theory of nerve excitation is made by comparing it with Hill's theory in fitting the experimental data on threshold-frequency relation, optimum frequency (v0) for nerve excitation and strength-duration relation. Decrease of v0 and increase of all the time constants (Hill's λ andk, Wei'sT 2 and spike durationw) with decreasing temperature are interpreted on the basis of the dipole relaxation timeT 2 but inexplicable from Hill's theory or any other existing theory. The closeness ofk,T 2 andw values is explained. A variety of experimental results obtained by others is discussed. Finally, a comparison is made between the Hodgkin-Huxley equations and the quantum theory. Most of the facts (electrical and non-electrical) tend to support the thesis that nerve excitation is a macroscopic expression of quantum transitions of dipoles between energy states.

Notes: Abstract In the periodic Leslie model the asymptotic period of total population is a divisor of the asymptotic period of the population vector. Under reasonable circumstances these periods are identical.

Notes: Abstract A number of biological branching systems, such as the bronchial and pulmonary arterial trees, are being investigated in an ongoing study in order to define their physiological properties. The technique involves the description of branching trees by the use of hierarchical systems of ordering, especially those described by Horsfield and by Strahler. During this work some mathematical properties of branching trees were demonstrated and these are described in this paper.

Notes: Abstract Luteinizing hormone (LH) is secreted continuously from the anterior pituitary gland. The concentration in the blood of this gonadotropic hormone plays a regulatory role in the development of puberty in both sexes, in the induction of ovulation in females, and in the production of testosterone in males. The secretion of LH is in turn controlled by luteinizing hormone releasing hormone (LHRH) secreted by the hypothalamus. LH and LHRH are removed from the blood by degradation and excretion. This hormonal system is modelled by a system of ordinary differential equations based upon specific physiological and biochemical assumptions current among experimentalists in this field. The one exception is the assumption that LHRH may bind reversibly to a serum protein; an analysis of the data shows that this or a similar mechanism is a crucial specification. Data on the serum levels of LH and LHRH in two human subjects were fitted using the model. The data consist of the transients and subsequent decays created by a bolus intravenous injection of LHRH.

Notes: Abstract A thick-wall incompressible, elastic sphere was used as a model for the diastolic rat left ventricle. A model for myocardial nonhomogeneity was derived assuming that fiber (circumferential) stress was independent of position in the ventricular wall. The theoretical implications of the resulting constitutive relations together with the spherical model were analyzed in the context of large deformation elasticity theory. It was found that muscle stiffness at a given level of uniaxial stress increased monotonically from the endocardium to the epicardium. In addition, fiber stress was found to be essentially a linear function of transmural pressure above a pressure of 6 g/cm2. It was also shown theoretically that neglecting the nonhomogeneity of the myocardium resulted in a state of stress which differed significantly from that predicted by the nonhomogeneous model. For example, at a transmural pressure of 14 g/cm2, fiber stress in the nonhomogenous model was equal to 17 g/cm2 while fiber stress in the homogeneous model varied between 100 g/cm2 at the endocardial surface and 2 g/cm2 at the epicardial surface. The change in muscle stiffness with position which characterized the nonhomogeneous model also tended to linearize the highly curvilinear radial stress distribution predicted by the homogeneous model at a given transmural pressure.

Notes: Abstract The phenomenon of axonal transport has been well documented (Ochs, 971; Lasek, 1970; and Grafstein, 1967). In a previous paper, we showed how diffusion alone could not account for this process. In this report we show that convection or convection with diffusion can account for the observed build-up of material. By including a first-order catabolic sequestration term, we are able to offer an understanding of the several apparent rates of transport with the same underlying velocity and variable sequestration.

Notes: Abstract It is known that the Lotka-Volterra coupled nonlinear differential equations for a two-species prey-predator ecosystem possess a periodic solution, although its exact form is not yet obtained analytically. The conventional linearization approximation for solving these nonlinear equations leads to a harmonic oscillator whose frequency depends only on the intraspecific coefficients. We propose here a prescription for obtaining nonlinear correction to the linear frequency by using the Hamilton-Jacobi canonical formalism of classical mechanics. It is found that the first-order correction, which also involves interspecific parameters, exhibits the basic qualitative features of the nonlinearity.

Notes: Abstract Environmental safety testing typically requires procedures for extrapolating from the relatively high experimental to the very low use doses of potentially harmful substances. In the present paper, a stochastic mammillary compartmental model for environmental safety testing is proposed and extrapolation procedures based on its dose-response relationship are developed. The proposed model is a direct generalization of one of the basic safety models, the one-hit model, in that a harmful reaction is assumed to occur if at any time any of the peripheral compartments attains a specified threshold of particles. Consideration of a closed model yields an upper bound on the probability of attaining a certain threshold level, thus providing a conservative procedure for extrapolating to a low dose, while a lower bound obtained from a related open model provides a useful monitoring device as to the sharpness of the upper, bound. The extrapolation procedure is illustrated with simulated data and approximations for initial values are developed.

Notes: Abstract By using chromosome images as a framework, algorithms for finding most dissimilar images are presented and illustrated by examples. In terms of angles, a chromosome image consists of two exterior biangles and two interior biangles. Biangles are defined and classified into 180° biangles, 〉180° biangles and 〈180° biangles. The dissimilarity of biangles and its geometric interpretation together with various properties of biangles are also presented. The results may have useful applications in pattern recognition, scene analysis, information storage and retrieval, artificial intelligence and fuzzy set theory.

Notes: Abstract The Volterra equations which represent competitions between two species are utilized to examine the phenomenon of boundary formation between two species of plants. The set of stable stationary points for these equations is determined and is illustrated in a product space of parameters and dynamical variables. The stages of boundary appearance and succession are visualized by considering slow changes of the parameters as functions of time and space.

Notes: Abstract Bayesian image processing formalisms which incorporatea priori information about valued-uncorrelated and valued-correlated (patterned) source distributions are introduced and the corresponding iterative algorithms are derived using the EM technique. Striking improvement in image processing is demonstrated when applying these algorithms to Poisson and Gaussian randomized data in one-dimensional cases.

Notes: Abstract For a broad class of interconnected nonlinear systems, this paper develops a complete design methodology for decentralized variable structure control. Specifically, the paper sets forth design schemes for local switching surfaces and the related local switched feedback gains which together force the original nonlinear interconnected system to behave as a reduced order interconnected equivalent system having a desired response such as stability, tracking, or prespecified eigenvalues. Also developed is a numerical algorithm for constructing the switched local feedback gains. A simple nonlinear example illustrates the control strategy.

Notes: Abstract A new stabilization method of a large-scale dynamic system, consisting of a set of interconnected subsystems, is presented in this paper. The topology of the interconnected subsystems is given as a network containing nodes with only one ingoing link, and none, one, or more outgoing links. Here, when the notion “node” is used a subsystem is assumed, and the links stand for the subsystem interconnections. The stabilization method is made only by the use of local linear state feedback around each subsystem, in order to satisfy constraints given in the problem. The interconnections among the subsystems are assumed to be nonlinear, time-varying. According to the topology of the large-scale system, the method of stabilization is hierarchic, one proceeds from node to node, and is applicable from a computer standpoint. A design algorithm follows directly, and can be made using the Generate and Test method for each subsystem independently, thus enabling designers to use a computer which has a video terminal as a peripheral unit and providing a possibility for interactive applications.

Notes: Abstract The estimation of covariance matrices which are structured, for example, of Toeplitz type, from measurement data is considered. The problem is considered in the context of array beamforming, and various methods of estimation are derived and compared, such comparison including consideration of the behavior of the estimate in beamforming applications.

Notes: Abstract The planar least-squares inverse (PLSI) polynomials are used for stabilization of two-dimensional unstable recursive filters. In order to obtain the PLSI polynomials, the main work involved consists in forming a set of linear equations and then solving them. In this paper we present an efficient and simple method to form the necessary set of linear equations (i.e., the required coefficient matrix) for a chosen pattern and order of the desired PLSI polynomial, starting from the denominator polynomial of a two-dimensional unstable recursive filter.

Notes: Abstract This paper shows how to use orthogonal functions to invert singular (i.e., generalized state-space) systems. The approach is to express the inverse system itself as a singular system, and then to apply the theory of orthogonal functions to convert that differential-algebraic system to a purely algebraicgeneralized Lyapunov equation whose solution yields the input of the original system given its output. Both left and right inversion are treated. Necessary and sufficient conditions for the existence and uniqueness of the generalized Lyapunov equation are derived, and a generalizedQZ algorithm is given for its efficient solution. It is also shown that the coefficients in the Walsh function expansion may be approximately found using an FFT-type butterfly network. These results provide both an extension in theory, by investigating the properties of a new Lyapunov equation, and an extension in the implementation of system inversion, by providing a scheme which applies to generalized state-space systems and uses an unconventional approach which may prove to be a useful contribution.

Notes: Abstract In digital communication networks, a special class of complex biquad recursive digital filters called orthogonal filters is increasingly being used. The separate effects of the overflow and quantization nonlinearities on these orthogonal filters' responses have been investigated [5], [6]. In this paper we examine the zero-input stability properties of the actual orthogonal filter having both overflow and quantization nonlinearities. The overflow nonlinearities considered include saturation, bit-by-bit inversion, zeroing, and modulo 2 arithmetic. The quantization techniques used may be roundoff, magnitude, or value truncation. An example demonstrates the adverse coupling effect between the overflow and quantization nonlinearities. Two criteria are therefore derived to ensure asymptotic overflowstability of the filter in the presence of quantization. These criteria have been translated to the coefficient plane; various regions corresponding to different minimum wordlengths required to ensure decoupling of the overflow and quantization phenomena have been derived.

Notes: Abstract The identification problem for electromagnetic objects excited by transients is discussed. Several classes of models are reviewed, and an output error model is selected. An algorithm for solving the transient identification problem using this model is presented, and some of the issues connected with its use are considered. Examples of the application of this algorithm to electromagnetic data are given.

Notes: Abstract A stochastic model for the population regulated by logistic growth and spreading in a given region of two-or three-dimensional space has been introduced. For many-species population the interactions among the species have also been icorporated in this model. From the random variables that describe stochastic processes of a Wiener type the space-dependent random population densities have been formed and shown to satisfy the Langevin equations. The Fokker-Planck equation corresponding to these Langevin equations has been approximately solved for the transition probability of the population spreading and it has been found that such approximate expressions of the transition probability depend on the solutions of the deterministic equations of the diffusion model with logistic growth and interactions. Also, the stationary or equilibrium solutions of the Fokker-Planck equation together with the special discussion on the pattern of single-species population spreading have been made.

Notes: Abstract The notion of an evolutive hierarchical system proposed in this paper is a mathematical model for systems, like organisms, with more or less complex objects. This model, based on category theory, retains the following characteristics of natural systems: they have an internal organization consisting of components with interrelations; they maintain their organization in time though their components are changing; their components are divided into several levels corresponding to the increasing complexity of their own organization, and the system may be studied at any of these levels (e.g. molecular, cellular...). The state of the system at a given instant is modeled by a category whose objects are its components, the state transition by a functor, a complex object by the (direct) limit of a pattern of linked objects (which describes its internal organization). The properties of limits in a category make it possible to ‘measure’ the emergence of properties for a complex object with respect to its components, and to reduce the study of a hierarchical system to that of its components of the lowest degree and their links. Categorical constructions describe the formation of a hierarchical evolutive system stepwise, by means of the operations: absorption of external objects, destruction of some components, formation of new complex objects.

Notes: Abstract An algorithmic formulation is presented for the inference procedure concerning lineage models. The problem is to find lineage rules from observed sequences of tree structures under the assumption that no interactions take place in the course of development and that sufficiently frequent observations are available at equal time intervals. The underlying structural pattern is taken to be a OL system, and the goal is to find propagating and deterministic OL schemes with minimal properties satifsying certain biological reliance criteria. Upper bounds have been found for the complexity of the inference algorithms.

Notes: Abstract A statistical theory of non-equilibrium fluctuation in Volterra-Lotka systems has been presented on the basis of the technique of statistical linearization of non-linear coupled stochastic differential equations.

Notes: Abstract Under certain basic assumptions the branching pattern of dendrites can be modeled as a Galton-Watson process in varying environment. Using results from graph theory we compute the probability distributions, expectations and variances for biologically significant variables such as the number of (intermediate and terminal) branches, the maximum number of orders, etc., together with the limit behavior of these quantities. Furthermore, the probability measure induced by the Galton-Watson process on the set of all trees is calculated. The measure assigns to any set of branching patterns the probability that it is realized by a certain process, which is completely described through the bifurcation probabilities.

Notes: Abstract Theoretical models for DNA repair mechanisms are constructed. Reliability studies considering the living cell as a repairable system are done. The DNA repair process is discussed along with applications and comparison with available experimental data.

Notes: Abstract A mathematical model has been developed to simulatein vivo transmural accumulation of an intravenously injected tracer in the aortic wall of experimental animals. Parameters have been included to represent the following processes that affect tracer distribution: permeation of the blood-tissue interface, diffusion through the layers of the artery wall,convective solute drag through the same, and degradation. Of particular interest for thein vivo situation situation is the inclusion of boundary conditions that account for the variation in the plasma concentration of injected tracer as a function of time. Two analytical solutions are presented. The first describes a system in which two boundaries must be delineated; it pertains if the tracer is allowed to circulate until it enters the avascular media of the artery wall both across its luminal boundary and from the capillaries in its outer layer. The second applies to shorter duration experiments in which entry across only the luminal boundary is considered. A limiting case of the solution for short circulation times is presented, compared with a previously published solution, and examined for its potential utility in parameter estimation. Because of its treatment of time-dependent boundary conditions, the model has unique application toin vivo experiments related to macromolecular transport in atherosclerosis that may otherwise elude proper interpretation.

Notes: Abstract Several current reaction-diffusion mechanisms have been proposed as models for morphogenesis in the Turing (1952,Phil. Trans. R. Soc. Lond. B 237, 37–72) sense. We introduce and exploit a quantity, we have termed heterogeneity, which allows us to elaborate the differences between the various models with regard to spatial pattern formation. It is shown that this quantity provides a concise view for the comparison of theoretical models with experimental observations. Two model mechanisms are treated explicitly both for linear and for biased diffusion.

Notes: Abstract The optimum number of total capillaries in the whole human body was estimated from the analysis of the efficiency for oxygen (O2) transport in the vascular-tissue system. We used a tissue model composed of uniform spheres in which O2 diffuses from the capillary located at the centre of each sphere towards the surrounding tissue consuming O2 at a constant rate. The tissue mass supplied by a single capillary was estimated as the area of positive O2 concentration under a given condition of capillary flow and O2 consumption rate. Total tissue mass was determined as the function of the capillary numbern and the total blood flow. On the other hand, the energy cost required to maintain the vascular system withn terminals was assessed by using the minimum volume model (Kamiya and Togawa,Bull. math. Biophys. 34, 431–438, 1972). The efficiency of the entire vascular-tissue system was evaluated by calculating the ratio of total tissue mass/cost function. The result of the calculation using physiological data of cardiac output and O2 consumption for an average human adult during a heavy exercise revealed the maximum efficiency occurring at the capillary number 3.7×1010 which coincided well with its normal range of physiological estimates (3.2×1010–4.2×1010). We concluded that the entire vascular-tissue system is constructed so as to attain the highest efficiency in O2 transport at its maximum activity.

Notes: Abstract A kinetic model involving synthesis of proinsulin in the rough endoplasmic reticulum, maturation through the Golgi apparatus and granules, with conversion to insulin, is proposed to account for data on the amount of insulin and of proinsulin both secreted during various time intervals and remaining in islets. Introducing three compartments for granules makes it possible to account for the measurement of both hot (pulse labeled with tritiated leucine) and cold proinsulin and insulin over a period of 21/2 hr under constant glucose. Data from islets from animals pretreated with tolbutamide are also presented and modeled. The model is then expanded so that it can be successfully applied to available data on the effects of a period of glucose deprivation on secretion of both hot and cold hormone. Parameters have essentially the same values, where they overlap, as were obtained (Landahl and Grodsky, 1982Bull. math. Biol. 44, 399–410) from insulin secretion by perfused rat pancreas stimulated by a variety of temporal patterns of glucose concentration.

Notes: Abstract Two models of binary tree growth are examined in terms of the Strahler order branching ratio (Rb) and the types of vertex produced during growth, and their inter-relationship. The sequential growth model is that described by Van Pelt and Verwer (1985,Bull. math. Biol. 47, 323–336) in which random growth occurs according to attributed probabilities on terminal or internal segments, one branch at a time. This model generates values ofRb≥3. The synchronous growth model is new and permits more than one segment to branch at a time, again randomly with attributed probabilities. This model generates values ofRb≥2 and in particular, when only terminal branching is permitted, gives 2≤Rb〈3. Such a model might explain the branching in the human bronchial tree, in which 2.5≤Rb≤2.8. Our synchronous model is an alternative to the centrifugal-order-dependent sequential model of Van Pelt and Verwer.

Notes: Abstract From an energy budget of a deciduous plant leaf in moderate conditions, entropy fluxes into or out of the leaf due to solar radiation, infrared radiation, evaporation of water and heat conduction are calculated. Net entropy flow into the leaf is negative. On the assumption that the entropy in the leaf is in a steady state, the entropy production in the typical deciduous leaf in moderate conditions [the solar energy absorbed by both sides of the leaf isE solar=0.0602 (J cm−2 s−1)] becomesS prod=1.8×10−4 (J cm−2 s−1 K−1). The positiveness of the entropy production shows that the Second Law of Thermodynamics certainly holds in the plant leaf. Entropy productions in other conditions are also calculated. The entropy production in the leafS prod becomes a linear function of the solar energy absorbed by the leafE solar:S prod≈-(29.5E solar)×10−4. A theorem is presented: the entropy production in plant leaves oscillates during the period of one day, paralleling the daily solar energy absorbed by leaves.

Notes: Abstract Molecular biologists strive to infer evolutionary relationships from quantitative macromolecular comparisons obtained by immunological, DNA hybridization, electrophoretic or amino acid sequencing techniques. The problem is to find unrooted phylogenies that best approximate a given dissimilarity matrix according to a goodness-of-fit measure, for example the least-squares-fit criterion or Farris'sf statistic. Computational costs of known algorithms guaranteeing optimal solutions to these problems increase exponentially with problem size; practical computational considerations limit the algorithms to analyzing small problems. It is established here that problems of phylogenetic inference based on the least-squares-fit criterion and thef statistic are NP-complete and thus are so difficult computationally that efficient optimal algorithms are unlikely to exist for them.

Notes: Abstract A method, based on symmetry, is suggested for determining the information content of systems. A comparison has been made between the information for symmetry, topology, and chemical composition. The new information measure increases when the asymmetry of the molecules and the number of atoms in the latter increases. It can distinguish between different molecular conformations, and give a linear correlation with the absolute entropy for homologous series of chemical compounds.

Notes: Abstract The micromorphic theory of Eringen is applied to study the tube flow of blood. The blood is considered to be a deformable suspension, with constitutive relations of the form of those of simple microfluids. By means of energy consideration, a relation is established between the local concentration parameter and the measure of rotationality involving both macro-and micromotions. The tube flow problem is then solved with some analyses on viscosity coefficients and boundary conditions. The results obtained indicate an integrated explanation of various important physical phenomena associated with blood flow, such as the tube size dependence of the apparent viscosity and the non-uniform concentration distribution over a tube cross section.

Notes: Abstract By observing that the n-tuple of rate functionsQ(c) is orthogonal to the c-space gradients of each of the (n - 1) constants of the motion Φ v (c), a generic canonical expression for the rate functions is given in terms of the exterior product of the gradients of the (n - 1) Φ v 's. For models withQ so prescribed from the outset, an analytical general solution is obtainable directly for the system of autonomous ordinary differential equations dc/dt =Q(c). Thus, the generic canonical expression for the rate functions can be utilized to construct analytically solvable models for interacting biological species, as ilIus~rated by examples here.

Notes: Abstract A model for the extraocular plant of the human visual eye tracking mechanisms is discussed. Its sensitivity to variation of controller signal nervous activity is studied in order to determine the type of activity that yields realistic simulations characteristic of typical saccadic eye movements.

Notes: Abstract Mathematical models of predator-prey systems in which the prey species has a three-stage life cycle are studied. Certain stages of the prey life history are allowed to use younger stages as food. It is shown that sufficiently restricted cannibalism can result in an increase in the numbers of adult prey on a sustained basis when cannibalism decreases the vulnerability of a stage subject to predation or increases overall productivity.

Notes: Abstract A general analysis is presented for the thermal behavior of a biological tissue. Energy transport by the circulatory system is assumed to be represented by a modified Fick's law. General boundary conditions are assumed for the two-dimensional model and solutions are obtained for rectangular, cylindrical, and spherical geometries. The effects of blood perfusion rate, metabolic rate, arterial temperature and heat exchange with the environment are considered. Results indicate a region of almost constant temperature in the deeper layers of the tissue and reaffirm the important role which blood flow plays in maintaining homeostasis.

Notes: Abstract The oscillatory aspect in a system having two steady states is studied theoretically using a model of excitable nerve membrane. The condition for the occurrence of oscillatory instability is discussed on the basis of the kinetic picture of nerve excitation in consideration of the non-Markoffian effect caused by ion transport in the system. Small oscillations around a steady state as well as a giant fluctuation between two states are obtained. Results are compared with experiments carried out with squid giant axons perfused intracellularly.

Notes: Abstract An expression for the variance in birth volumes during balanced growth of a cell population is derived. The requirement of this expression being positive and finite allows a discussion of some of the requirements imposed on the mechanisms of growth and division.

Notes: Abstract The phenomenon of axonal transport of material has been well documented (Ochs, 1971; Lasek, 1970; and Grafstein, 1967). This report seeks to establish the role of diffusion—if any—in such a transport process. We report that diffusion cannot account for the observed build-up of material as reported in the literature.

Notes: Abstract The question of how to fit a general cubic model of a multicomponent, interactive growth system to observed data is addressed. A multidimensional-polynomial type of regression analysis is used, with a least-squares criterion. By testing the scheme on a problem with known solution, the way in which the accuracy of the results varies with the number of datum points used is investigated in an heuristic manner.

Notes: Abstract The describing function method is used as a guide to the behaviour of the solutions of the equations of Danziger and Elmergreen, proposed as a model of periodic catatonia. The method suggests that whenever the equilibrium point is unstable it is surrounded by a stable closed periodic orbit. This is confirmed in specific cases by computation.

Notes: Abstract A theory of ambiguous pattern perception is formulated. This theory proposes a feature selector (field of attention) based on the time-sequential discrete property of the attention, a short-term memory for storage of the selected features, and a displayer (perception) to display the consecutively stored features. Since the selected features continuously enter, and since the features can only be stored in the short-term memory for a short period, the features which can be displayed in the displayer vary with time. When all the essential features belonging to one pattern happen to be in the displayer, the picture is perceived to be that pattern; when all the essential features belonging to another pattern happen to be in the displayer, then the picture is perceived to be the other pattern. Thus the picture appears to vary with time and alternate between two patterns. A numerical calculation is presented.

Notes: Abstract The thermodynamics of irreversible processes is derived from the principles of dynamical field theory independently of all elements of thermostatics, in particular the assumption of local equilibrium. Field thermodynamics proceeds from the premise that all driving forces experienced by the molecules in a continuum are conservative and arise from scalar potential functions. Dynamically the temperature potentialT is no different from the pressure potentialp. A field is converted to a force upon multiplication by a scale factor. A potential is converted to potential energy by the same scale factor. To scale the field −∇p to the force per mole of molecular speciesk, the partial molar, volume $$\bar V_k $$ is the scale factor. Similarly the partial molar entropy, $$\bar S_k $$ , scales the temperature field. The transition from the scale factors (which are physical parameters) to the systemic variables, for example $$\bar S_k \to s\left( {x,y,z;t} \right)$$ , is not trivial. From the dynamics and the structure of the derived potential energy function are inducted the conjugate variables such as (p, V I) and (T, s). The meta-mechanical properties of the thermal variables (T, s) are discovered via the local First Law of Thermodynamics, which relates internal energy, thermal flux, and work, and from the local Second Law, which prescribes, the possible partitions of internal energy between kinetic, potential, and thermal energies. From the form of the potential energy come Maxwell's relationships. From the energy partition comes the equation of continuity for entropy, with its important source term. In contrast to earlier theories of irreversible thermodynamics, the dissipation function does not include the stress tensor, a constitutive parameter.

Notes: Abstract The transfer of solute through a membrane separating two aqueous solutions is studied with the time-dependent diffusion equation for composite media. By introducing new independent and dependent variables it is shown that the differential equations and boundary conditions can be transformed into a dimensionless form which does not explicitly depend on the diffusivities of the media. Laplace transforms are used to derive explicit solutions for the solute concentration as a function of position and time. It is shown that at large time the concentration approaches the equilibrium distribution exponentially. Explicit results are given for the decay time as a function of the parameters of the system. In addition, an accurate and simplified expression is derived for the decay time for the case of small membrane permeability. The accuracy of the analytic solutions for the concentration profiles is tested by comparing them with numerical results obtained by solving the diffusion equations by the method of finite differences. Excellent agreement is found.

Notes: Abstract Physiological systems are often modelled by a set of compartments. Alternatively they can be described by the diffusion-convection-reaction equations governing distributed systems. The problem considered here is that of identifying a continuously changing input of some metabolite )tracee), endogenous to the system and hence inaccessible, when a nonlinear or time-varying component is also introduced into the loss parameter, as for example through feedback mechanisms. A tracer is used to determine the steady-state impulse response under time-invariant, linear conditions. A known input of tracer is also administered when the system is driven out of steady state. The integral equations developed utilize the predetermined impulse response, the measured concentrations of both tracer and tracee (output) in some region of the system to estimate the changing loss parameter and the unknown input in a continuous fashion.

Notes: Abstract The equilibrium distribution for a generalQth-order multivariate reaction system is studied. The state transition intensity matrix is developed and examples are given for small numbers of reaction components. A closed-form expression for the equilibrium distribution for systems which are symmetric with respect to the order of component reactions is presented. Numerical examples for three component systems are discussed.

Notes: Abstract The model recently proposed by Dreitlein and Smoes for oscillatory kinetic systems is studied. Diffusion of the oscillating species is taken into account, and bounds on the total number of individuals of each species are determined for both two- and three-dimensional finite regions with various boundary conditons applied. It is found that in general the effect of diffusion on the system behavior is to reduce the maximum possible radius of limit cycles. In particular, in some cases global limit cycle behavior is precluded.

Notes: Abstract The Kedem-Katchalsky equations for fluid flux across membranes may not be adequate for large solvent flows. In particular, for an example of two membranes in series, it is argued that they would predict physically unreasonable behavior. An alternate equation for solute flow is proposed for a simple sieving membrane. For the same example, this equation predicts more physically reasonable results.

Notes: Abstract In order to understand the abnormal flow conditions of blood in a locally constricted blood vessel, the analytical results are obtained for the oscillatory flow of blood which behaves as a Newtonian fluid. It is here assumed that the surface roughness is cosine-shaped and the maximum height of the roughness is very small compared with the radius of the unconstricted tube. Numerical solutions are presented for the instantaneous flow rate, resistive impedance, wall shear stress and phase lag.

Notes: Abstract In order to better understand the effect of initial stress in blood flow in arteries, a theoretical analysis of wave propagation in an initially inflated and axially stretched cylindrical thick shell is investigated. For simplicity in the mathematical analysis, the blood is assumed to be an incompressible inviscid fluid while the arterial wall is taken to be an isotropic, homogeneous and incompressible elastic material. Employing the theory of small deformations superimposed on a large initial field the governing differential equations of perturbed solid motions are obtained in cylindrical polar coordinates. Considering the difficulty in obtaining a closed form solution for the field equations, an approximate power series method is utilized. The dispersion relations for the most general case of this approximation and for the thin tube case are thoroughly discussed. The speeds of waves propagating in an unstressed tube are obtained as a special case of our general treatment. It is observed that the speeds of both waves increase with increasing inner pressure and axial stretch.

Notes: Abstract The paper introduces a basic mathematical form, which is characteristic of a number of linear one-dimensional diffusion equations with coefficients being represented as simple polynomials in the spatial coordinate. A number of particular diffusion equations are introduced and their corresponding exact mathematical solutions are obtained.

Notes: Abstract In this paper we examine a set of nonlinear rate equations (devised by M. Eigen (1971)) which describe the process of selection in a collection of self-reproducing, macromolecular information carriers. We construct exact solutions to the equations for the case of constant rate parameters and constant error distributions. The solutions allow the direct assessment of the effect of mutations on the “selective value” parameters discussed by Eigen as well as the distribution of the molecular species selected in steady state. In addition we show that the selection process may be characterized by an extremal principle.

Notes: Abstract The electromagnetic molecular electronic resonance (EMER) frequencies of the molecular chains of α-chymotrypsin are calculated. The chain length relations and coupling positions suggest a possible energy transfer, at the EMER frequencies, from one chain to the other. Photon enzyme activation data indicate that the energy corresponding to the EMER frequencies of its molecular chains is used by α-chymotrypsin for its enzyme function.

Notes: Abstract Based on Wei's dipole flip-flop model and with the assumption that the dipole is coupled to the membrane matrix, the cathode-make-excitation, the anode-break-excitation and the cathode-gap-excitation can be explained in a systematic way. The istrength-duration relations for these three processes are derived.

Notes: Abstract It is usually stated that only systems with multiple steady states can exhibit hysteresis. For protein conformations, this would violate the idea that, in fixed environments, primary structure uniquely determines tertiary structure. It is shown that hysteresis-like phenomena can be exhibited by systems possessing only a single steady-state configuration. This property is placed in a more general theoretical setting of recognition and classification systems, and some implications for processes such as memory, learning and pattern generation (morphogenesis) are discussed.

Notes: Abstract The natural historical literature contains a considerable body of work which indicates that harvesting (or predation) can alleviate competitive instabilities. In order to arrive at an understanding of this, the appropriate bifurcation structure for a rather general family of two-dimensional competitive systems is here investigated. The results of this analysis suggest that, in more complicated ecosystems with many competing species, (1) there is a good chance that harvesting at moderate rates will increase species diversity if one species is dominant in the unharvested system, while an increase in diversity is not likely to result from harvesting from a system with no dominant species, (2) whenever harvesting does increase species diversity, maximal diversity will occur at moderate harvesting rates, with less diversity at both very high and very low harvesting rates.

Notes: Abstract The Weber-Fechner law (response of an organism is a linear function of log of stimulus) has been widely used for description of physiological and psychological data for many years. It is shown here that the Weber-Fechner law is derivable in a simple way from the Elovich equation (−dx/dt=m exp(nx), which is observed experimentally in numerous physiological and biochemical systems, and which has a simple derivation from solid state charge transport across interfaces in these systems. It therefore seems reasonable to interpret data conforming to the Webner-Fechner law to imply that the observed phenomenon is rate-limited by interfacial charge transport in the cell. By a similar analysis, the Loewenstein equation, which may be considered an exact form of the Weber-Fechner law applicable to data over a wider range of values of the variables, is derived from a more exact form of the Elovich equation.

Notes: Abstract A simple similarity transformation of the Leslie matrix renders certain properties obvious. In particular the characteristic polynomial, characteristic vectors and principal vectors can be explicitly written out. Bounds for the dominant root are given.

Notes: Abstract A mathematical description of the process of adherence to glass of the rat thymus lymphocytes is presented. It is based on the result of previous work in which the process was studied in the course of the flow of the51Cr labelled cell suspension through the glass-bead column. The concentration of cell suspension, flow velocity and medium temperature were constant; the experiments were performed with different lengths of the bead bed. The amount of cells captured on the glass beads' surface was calculated as a function of time. Two approaches to the mathematical description of the process are presented. The first one is based on the linear equation of kinetics of cell retention on the elementary thin layer and on the transport equation of the flow of the suspension through the column. In the second one, the differential equation of the adhesion, derived from the experimental data, is discussed.

Notes: Abstract A series of Bayesian image processing algorithms which incorporate various classes ofa priori source information in treating data which obeys Poisson and Gaussian statistics is derived using maximum entropy considerations. The standard maximum likelihood equations are shown to be a special case of Bayesian image processing when thea priori information about a source distribution φ j is solely that a non-vanishing probability for each element value φ j exists only in some finite interval,a j ≤φ j ≤φ j . Bayesian image processing equations for thea priori source information that all φ j are finite -∞〈φ j 〈∞ and each φ j distribution has a defined mean φ j and a defined variance σ j are derived. The Bayesian image processing equations are also derived when thea priori source information is that all φ j ≥0 and that each φ j distribution has a defined mean φ j and a defined variance σ j . The a priori source distribution constraint that a correlation exists among nearby elements is also considered. The results indicate improvement over standard methods.

Notes: Abstract The superficial capillary network of the gastric mucosa can be monitored for red blood cell velocity measurements by a microscopic technique. This network, however, reflects the blood flow in capillaries of more physiological interest, namely those passing by the acid-producing cells and emptying into the superficial network. It is, however, not possible to study these capillaries directly and therefore the problem is to determine in what way and to what degree blood flow measurements in the superficial network reflect the capillary flow of interest. A probabilistic approach where the movements of the red blood cells have been analysed, gives indications of determinable relations between observations on the superficial network flow and the flow passing the acid-producing cells.

Notes: Abstract The present paper describes a calculation technique to determine a value of Young's modulus for the cornea when the intraocular pressure, the shape of the cornea and the thickness variation along the cornea are known fromin vivo measurements. Twenty-eight persons were examined in both eyes and a mean value of Young's modulus for the cornea was calculated to 2.45×104±0.57×104 N/m2 (S.D.).

Notes: Abstract It is pointed out that the asymptotic general solution to the ϕ-model equation for a periodic carrying capacityK(t) andt≳r −1 is identical in form to the generalized logistic equation solution with a built-in developmental time delay τ(≲r −1) and associated parameter ranges of primary biological interest. In the case of the ϕ-model equation, the time delay is a purely dynamical consequence of the nonlinear form featured by the population growth rate.

Notes: Abstract The pulsatile flow of blood through arteries is investigated in this paper by treating the blood vessel as a thin-walled anisotropic, non-linearly viscoelastic, incompressible circular cylindrical shell; nonlinearities of the flow of blood are also paid due consideration. The displacement components of the vessel wall are obtained from the equations of equilibrium which have been linearized by employing the principle of superimposition of a small deformation on a state of known finite deformation. The influence of the wall deformation on the flow properties of blood, has been accounted for by considering suitably formulated continuity conditions. A finitedifference scheme is employed for solving the flow equations together with the boundary and initial conditions by using the locally measured values of pressure and pressure gradient. Numerical results obtained for the velocity profile of blood flowing in a canine middle descending thoracic aorta have been presented through figures.

Notes: Abstract The “protocell” is a mathematical model of a self-maintaining unity based on the dynamics of simple reaction-diffusion processes and a self-controlled dynamics of the surface. In this paper its spatio-temporal behaviour far from the stationary structure is investigated by means of a boundary layer approximation. It is shown in detail how a simplified and mathematically feasible equation can be derived from the original parabolic problem. It turns out that the known instability which is initiated in the linear region around the stationary structure is continued further in the direction to a division by nonlinear dynamics.

Notes: Abstract To describe stationary-state kinetics of solute clearance from an interstitial space with an initial uniform concentration, an explicit solution of the Sangren-Sheppard model and a stratagem for an explicit solution of the Johnson-Wilson model are presented. In both cases expressions for computing the exact intravascular and extravascular concentrations at any time and location are complicated and inconvenient. Simpler and far more convenient formulae for determining upper and lower bounds on solute concentrations and on the pseudo-first-order clearance rate ‘constant’ (k) are derived. For the Johnson-Wilson model, the bounds are so tight thatk, for example, can be estimated with considerable accuracy whenever the capillary blood flow exceeds the permeability-surface area product.

Notes: Abstract Different types of random binary topological trees (like neuronal processes and rivers) occur with relative frequencies that can be explained in terms of growth models. It will be shown how the model parameter determining the mode of growth can be estimated with the maximum likelihood procedure from observed data. Monte Carlo simulations were used to study the distributional properties of this estimator which appeared to have a negligible bias. It is shown that the minimum chi-square procedure yields an estimate that is very close to the maximum likelihood estimate. Moreover, the goodness-of-fit of the growth model can be inferred directly from the chi-square statistic. To illustrate the procedures we examined axonal trees from the goldfish tectum. A notion of complete partition randomness is presented as an alternative to our growth hypotheses.

Notes: Abstract Compartmental models of biological or physical systems are often described by a system of “stiff” differential equations. In this paper an algorithm for solving a system with linear coefficients is presented that employs numerical inversion of the Laplace transform of the model equations. The inversion algorithms and Gear's backward differentiation method are compared for two stiff test problems and a differential system governing a 27-compartment model of bile acid transport and metabolism. The inversion algorithm is reliable, requires modest computation time on a desktop computer and provides better accuracy than Gear's method, especially for the extremely stiff example.

Notes: Abstract Proposed model geometries are described for three types of axisymmetric and four types of flattened sickled erythrocytes. The axisymmetric types are designated as concave-convex, convex-convex and parabolic, while the flattened types carry designations of concave-convex, convex-convex, parabolic and S-shaped. Formulae are provided for the volume and surface area of each type and a figure showing the profile of each is included. Measurements of actual sickled cellsin vitro could be used to find values for the volume and surface area of each type of cell by using parameter values for the appropriate model geometry. These would give close approximations that could be useful in clinical therapies and laboratory investigations for sickle cell anemia. Surface area to volume ratios can also be found to a close approximation for each cell type.

Notes: Abstract In this paper we consider the stability property of single-species patches connected by diffusion with a within-patch dynamics of Volterra type and with continuous time delays. We prove that this system can only have two kinds of equilibria: the positive and the trivial one. By the assumption that the delay kernels are convex combinations of suitable non-negative and normalized functions, the linear chain trick gives an expanded system of O.D.E. with the same stability properties as the original integro-differential system. Homotopy function techniques provide sufficient conditions for the existence of the positive equilibrium and for its global stability. We also prove the local stability of any positive equilibrium and the local instability both of positive and trivial equilibria. The biological meanings of the results obtained are compared with known results from the literature.

Notes: Abstract Currently applied three-copartment models for analyzing kinetic data derived fromin vivo positron emission tomographic (PET) studies of radioligand-neuroreceptor interactions require assumptions which may not be strictly valid. Such assumptions include very rapid kinetics for nonspecific binding and the absence of multiple specific receptors or subtypes. Computer simulations, based on an exact analytical solution of the relevant differential equations, indicate the numerical errors that can arise when the assumptions are invalid. We propose a fourcompartment model which requires fewer assumptions. A simple relationship is derived for expressing the microscopic rate constants of either the three- or four-compartment model as explicit functions of the experimentally-observed macroscopic rate constants. This could eliminate the need for time-consuming, iterative, non-linear, curve-fitting approaches and numerical integration. The usefulness of the four-compartment model is limited, however, by the sensitivity and temporal resolution of current PET imaging devices.

Notes: Abstract Various diffusion processes employed for modelling logistic growth are briefly summarized. A discrete-time, discrete-state space stochastic process for population growth is proposed and analyzed with either Bose-Einstein or Maxwell-Boltzmann statistics for the distribution of offspring in available sites in a restricted region. A diffusion approximation is constructed, which differs from those previously employed. The logistic law is a natural deterministic analog of the diffusion process.

Notes: Abstract In many biophysical and biochemical experiments one observes the decay of some ligand population by an appropriate system of traps. We analyse this decay for a one-dimensional system of randomly distributed traps, and show that one can distinguish three different regimes. The decay starts with a fractional exponential of the form exp[−(t/t 0)1/2], which changes into a fractional exponential of the form exp[−(t/t 1)1/3] for long times, which in its turn changes into a pure exponential time dependence, i.e. exp[−t/t 2] for very long times. With these three regimes, we associate three time scales, related to the average trap density and the diffusion constant characterizing the motion of the ligands.

Notes: Abstract Some key experiments of artificial production ofsitus inversus viscerum are briefly reviewed and a two-step mechanism for the explanation of the systematic asymmetric visceral arrangement in vertebrates is proposed. A two-variable reaction-diffusion system displaying a symmetry-breaking bifurcation is considered, and it is demonstrated that a slight asymmetry of the boundary conditions can give rise to a marked asymmetry in the resulting dissipative structure in both one-and three-dimensional systems. A criterion is formulated allowing classification of reaction-diffusion systems operating in a three-dimensional space with regard to their ability to incorporate slight asymmetries at the boundaries in the form of a chiral dissipative structure.