ALBERT

Notes: Abstract By studying the behavior of various tracer species in the lungs, one can assess many important characteristics which distinguish normal and abnormal function. Quantitative evaluation of function depends on the use of an appropriate model in conjunction with experimental data. A multi-compartment model is derived from mass balances to describe dynamic as well as (breath-averaged) steady-state transport processes between the environment and pulmonary capillary blood. The breathing cycle is divided into three time periods (inspiration, expiration, and pause) so that the model equations are discrete in time. No other model of tracer species transport in the lungs deals simultaneously with species dynamics, variable breathing pattern, distribution inhomogeneities, and non-equilibrium between alveolar gas and capillary blood. Models currently in the literature are shown to be special cases of the model presented here.

Notes: Abstract Local stability seems to imply global stability for population models. To investigate this claim, we formally define apopulation model. This definition seems to include the one-dimensional discrete models now in use. We derive a necessary and sufficient condition for the global stability of our defined class of models. We derive an easily testable sufficient condition for local stability to imply global stability. We also show that if a discrete model is majorized by one of these stable population models, then the discrete model is globally stable. We demonstrate the utility of these theorems by using them to prove that the regions of local and global stability coincide for six models from the literature. We close by arguing that these theorems give a method for demonstrating global stability that is simpler and easier to apply than the usual method of Liapunov functions.

Notes: Abstract The asymptotic behaviour of a logistic equation with diffusion on a bounded region and a diffusionally coupled delay is investigated. An equivelent parabolic system is derived for certain types of delays. Using a Layapunov functional, sufficient conditions for the global asymptotic stability of the constant steady state are obtained. When the global stability is lost, using Hopf's bifurcation theory, existence of travelling waves is shown for ring-like and periodic one dimensional habitats.

Notes: Abstract It was hypothesized in an earlier work that sensory perception can occur only when the perceiving system is uncertain about the nature of the event being perceived. In the absence of any uncertainty, perception will not take place. The response of the sensory afferent neuron (impulse transmission rate) was calculated using Shannon's measure of uncertainty or entropy. It will now be shown that when the event being perceived is the position and momentum of a particle, Shannon's measure of uncertainty leads to the Heisenberg Uncertainty relationship.

Notes: Abstract It is not unusual for several classifications to be given for the same collection of objects. We present a method, called majority rule, which can be used to define a consensus of these classifications. We also discuss some mathematical properties of this consensus tree.

Notes: Abstract The dependence of the spatial concentration profiles of morphogens on a characteristic dimension is obtained by continuation techniques for Gierer and Meinhardt's activator-inhibitor model of morphogenesis. The study of the behaviour of the system during growth, where the linear and exponential increase of the characteristic dimension is considered, revealed that more complex patterns of morphogen spatial concentrations appear regularly in a reproducible way.

Notes: Abstract Computer models have been used by various authors to simulate both the growth of normal cellular tissue and the development of cancerous cells within normal tissue. As these models were the result of considerable idealization, the purpose of the present paper is to propose a model in which the degree of simplification is relaxed: the features of simultaneous growth, and cell growth whose rate depends on the free absorbing periphery of the cell are introduced. Simulation experiments have been conducted using the model, and the results are presented.

Notes: Abstract The theory of complementary variational principles is used to obtain maximum and minimum principles for a nonlinear model of heat conduction in the human head. Accurate variational solutions are obtained in illustrative calculations. The effect of nonlinearity is seen to be significant from a comparison with the linearized model.

Notes: Abstract A model for the nerve impulse due to Zeeman (1972) and based on catastrophe theory is compared with alternative models and criticisms of Zeeman's model by Sussmann and Zahler (1977, 1978) are assessed. The criticisms of Zeeman's motivation for his model are found to carry some weight. Sussmann and Zahler (1977, 1978) list numerous features of Zeeman's model which, they state, are not in agreement with experiment. These statements as they stand are largely erroneous, and the model still remains to be tested by a critical series of experiments. However, a detailed analysis reveals defects in Zeeman's model, not among those claimed by Sussmann and Zahler, showing that the explicit equations of the model cannot be correct. The possibility of a modified approach along similar lines and its ultimate adoption remains open.

Notes: Abstract The irreversible Michaelis-Menten reaction is studied by the use of the method of multiple scales. Three stages of the reaction are identified, one of which is studied in detail. The results are compared with those of two earlier analyses.

Notes: Abstract A numerical study of the coupled nerve fibre problem is given which verifies and extends the perturbation theory of Luzader. Pulses on adjacent fibres can couple together with two possible stable pulse separations.

Notes: Abstract A possible mechanism for effects of microwave radiation on the auditory system is the generation of field-induced forces at interfaces that divide materials of dissimilar electrical properties. A general expression for these “Maxwell stresses” is derived and then used to calculate the approximate magnitude of field-induced force within the organ of Corti during microwave exposure. Comparison of the results with data on the force needed to excite cochlear hair cells indicates auditory responses could be evoked by this mechanism at power densities near the threshold of rf hearing sensations.

Notes: Abstract A definition of homogeneity for neural networks is given which permits their construction as group quotients. The significance of this for neural dynamics is discussed.

Notes: Abstract The left ventricle is represented as a cylinder contracting both radially and longitudinally. A simple method is indicated to derive an expression for the rate of change of the kinetic energy of this three-dimensional model, which quantity can be used as an index for the study of the contractile behaviour of the myocardium. An application to the study of muscle mechanics is also indicated.

Notes: Abstract A perturbation method is proposed to calculate approximately the limit cycle type nonequilibrium steady-state resulting from periodic perturbation of coefficients of stable population systems; the two species Lotka-Volterra competition system is explicity studied and the results are formulated for general multi-species population systems. Avoidance of competitive or other types of exclusion of species in a periodic environment is indicated.

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 The theory of computational complexity and certain explicitly-stated hypotheses imply limitations on the information processing power of biological systems. Parallelism, special purpose organization, and analog mechanisms may provide speedup critical for life processes, but have little power in the face of exponential growth. We show that “polynomially simulatable” biological systems cannot exhibit dynamic behavior which produces the solution of an intractable problem. The argument implies that parallelism does not allow biological systems to defeat the exponential explosion, but rather is important because it allows polynomial time algorithms to be used more efficiently.

Notes: Abstract A correlation matrix analysis is applied to the base sequence of MS2 and ϕX174 in comparison with sets of simulated sequences with different degrees of constaint Significant differences between a codified sequence, and a statistical one in terms of the “correlation matrix” for sets of different length cannot be found. This result is analysed in terms of nucleotide sequences with different levels of informational content.

Notes: Abstract A method of calculating the volume of a tree distal to a cut at the origin of a branch, using branching, diameter and length ratios, has been developed. The method was applied to bronchial tree casts from human, dog, sheep, hamster, and rat lungs. It was found that the exponenta in the equation weight=k×diameter a is approximately equal to 3.0 in sheep lung casts, as found by Hooper (1977), but it is greater than 3.0 in casts from the other four species.

Notes: Abstract In this note we examine a continuous time version of a compartmental model introduced in a discrete time setting by S. R. Bernard. The model allows for more than one particle to leave the system at any time. This introduces additional randomness into the system, over the pure death system and this is reflected in the variance function.

Notes: Abstract The mean first passage time for free diffusion can be derived directly by solving a simple analogue steady state problem. In this problem the diffusion starting region is considered as a time independent source of diffusing particles and the diffusion target assumes the behaviour of a perfectly absorbing sink. It is shown here that the transit time between the source and the sink, which in this particular problem is equal to the ratio between the holdup of the system and the total flux, is identical to the Brownian movement concept of the mean first passage time for free diffusion. This established identity considerably facilitates the derivation and investigation of the timing of diffusion in complicated structures such as those commonly found in living organisms.

Notes: Abstract In this paper three stochastic models are developed for a class of two-compartment systems to analyse the randomness of the leaving process of the particles in the system. Results in closed form for the distribution of the leaving process of the particles in the system are given both for general and exponential sojourn time distributions and also in association with forward recurrence time distributions with and without Poisson input.

Notes: Abstract Two simple models are proposed and analysed, in which it is shown that the formation of a new polymer, resulting from an “error” in the template action mechanism of production of an old polymer, may compromise the stability of the initial system under specific conditions, in the context of prebiotic evolution. Autocatalysis is shown to be a “selective advantage”, enabling the “mutant” to dominate in concentration and even replace the initial polymer. The addition of a third molecule playing the role of a catalyst causes hysteresis effects.

Notes: Abstract The problem of extinction of the prey population in a microbial predator-prey interaction in a chemostat has been examined. Usual deterministic lumped parameter models were used for the dynamics of the chemostat for large numbers of the two populations; the generalized birth and death stochastic process was employed for the description of the random variations at small prey numbers. Extinction probabilities of the prey population were calculated for different holding times and chemostat volumes, and their dependence upon the growth parameters of the two populations was studied. It was found that extinction was possible when the Monod model was used for the specific growth rate of the predators as a function of the prey number density. On the other hand, the decrease of the feeding activity of the predators at low prey densities predicted by the multiple saturation model acts as a regulatory factor that prevents extinction of the prey. In view of the fact that extinction of the prey has never been observed in the laboratory, the latter model seems more appropriate to describe the dynamics of microbial predation.

Notes: Abstract During exposures of the eye to light, the choroidal circulation may have a regulatory influence on the retinal temperature. This is investigated using a mathematical model and a finite-difference technique. It is predicted that the choroidal blood flow a small effect on retinal temperature, which may be important.

Notes: Abstract A probabilistic model of a spatially localized, mutually exitatory (inhibitory) population of neurons is formulated to help explain average evoked potential and post-stimulus time histogram measurements. The model is based on the stochastic activity of single neurons within interactive masses of neurons which exhibit co-operative behavior. Macrostate variables corresponding to the above measurements are related through the model to features of neural operation at the individual and ensemble level. Steady-state solution are obtained and their physiological implications are discussed.

Notes: Abstract We consider a one-compartment system with stochastic transfer rate characterized either by Gaussian or by two-level jump process and study the time evolution of the (statistical) moments of the (random) amount of the substance present in the system. The effect of the coloured as well as of the white noise is investigated and it is found that the presence of stochasticity in the transfer rate parameter increases the relaxation time of the system. Finally, we obtain the conditions for the stability of the system in the moment sense.

Notes: Abstract A model is described in which damage to a single intracellular locus can lead to a tumorigenic transformation. Assuming a large number of independent intracellular loci to be at risk and assuming that damage to a locus sufficient to cause a tumorigenic transformation occurs with probability greater than zero for all doses greater than zero, leads to the use of the Weibull distribution to characterize the probability of a nonspontaneous tumorigenic cellular transformation occurring after exposure to a given dose of carcinogen. The excess lifetime tumor incidence (i.e., the proportion of tumor bearers) above the spontaneous incidence is used as an estimate of the non-spontaneous incidence and is characterized by a tumor incidence function that represents the probability of occurrence of one or more non-spontaneous tumorigenic cellular transformations amongN(D) independent surviving cells per individual, after exposure to a doseD of carcinogen. The tumor incidence function is fitted to published data for the excess tumor incidence after exposure of animals or humans to ionizing radiation and after exposure of animals to chemical carcinogens.

Notes: Abstract This paper deals with a stochasticn-compartment irreversible system with a non-homogeneous Poisson input and arbitrary residence time for each of the compartments. Results relating to the number of particles present in each of the compartments as well as the total number of particles present in the system at any time are derived. Further, explicit expressions for the auto covariance function for each compartment and the cross-covariance function between any two compartments with a given time lag are obtained. As a particular case, then-compartment irreversible system is analyzed with homogeneous Poisson input and exponential residence time distribution for each of the compartments. The possible applications of the model are discussed.

Notes: Abstract This paper deals with the pulsatile blood flow in the lung alveolar sheets by idealizing each of them as a channel covered by porous media. As the blood flow in the lung is of low Reynolds number, a creeping flow is assumed in the channel. The analytical and numerical results for the velocity and pressure distribution in the porous medium are presented. The effect of an imposed slip condition is also studied. Comparisons with the corresponding results for the steady-state case are made at the end.

Notes: Abstract The relationships that define the structure of a given ecosystem, social system, or even a physiological function can only exist if certain parameters are confined to a certain range of values. As the values change and exceed this given range the relationships are forced to change, and so produce a new pattern of relationships. The concept of a dynamic structure captures this potential for structural change in relation to a set of parameters. The precise definition of structure and allowable transformation constitutes the definition of a category. The total range of parameters associated with all the relevant structures provides a parameter space which is assumed to be a manifold. Maps with extra structure from the manifold to the category define dynamic structures. The domain of differential dynamic systems is the manifold, and a flow or movement across the manifold is associated with a series of structural transformations in the category. In some cases a structure outruns its parameter range, to be faced with an obstruction—an absence of possible transformations. Ways of studying such “obstructions” are considered along with the related problem of extending a dynamic structure beyond a previously given set of parameters. The cost or resistance of transformations is also studied. The concepts of dynamic structures are illustrated by the structural change of food webs and they are used in a necessarily qualitative fashion to study dominance structures of social orders and finally to speculate on the qualitative nature of evolutionary change of functional aspects of organisms.

Notes: Abstract The preceding paper (Thorn, 1981) has shown that in a linear pharmacokinetic system with a multimodal impulse response the peak drug level may sometimes be smaller with slower rates of injection. This paper presents two theorems on this paradoxical injection rate effect where the injection is a constant infusion of finite duration. The first theorem establishes a graphical method for determining whether a given impulse response will give a paradoxical injection rate effect; and the second establishes that the maximum paradoxical increase in peak drug level is by a factor of two. It is further shown that in order to approach this maximum paradoxical increase the impulse response must contain two isolated, sharp, narrow pulses of approximately equal area. Some examples of bimodal arterial dye-dilution curves from the literature are discussed as impulse responses; and there is also a discussion of the behavior of drug level maxima and minima at different injection rates.

Notes: Abstract This paper presents three theorems on the peak drug levels that result from injection into a linear pharmacokinetic system. As a preliminary, the “rate of injection” is defined in terms of time expansion or time contraction of the injection function (input). The first theorem then states that the peak drug level will not be greater when the rate of injection is slow than when it is fast, if the impulse response is unimodal. The second theorem sets limits for the time of the maximum drug level, in relation to the time of the maximum of the (unimodal) impulse response and the duration of the input. The third theorem defines conditions which assure a definitely lower peak drug level if the rate of injection is slower. A graphical method is suggested for determining the times and magnitudes of the peak drug levels that result from constant infusions of a fixed dose at different rates. An example is provided to show that if the impulse response is multimodal then the peak drug level may sometimes increase with a decrease in the injection rate.

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 A two urn Polya-type scheme is considered in whichr black balls (corresponding to the stable form of an element) are added to urn one at every stage and the same number of balls are removed at random at every stage from the same urn. In between these two operations, which form a stage or iteration, a fixed number of balls is exchanged at random between urns one and two. Urn one has a given initial number of white balls (corresponding to a radioactive form of the same element). The problem of interest is to study the stochastic aspect of the number of white balls remaining in urn one (and/or urn two) aftern iterations.

Notes: Abstract We obtain within the action-angle variable approach new expressions, involving the Dirac delta function, for time periods and time averages of dynamical variables which are useful for nonlinear biological oscillator problems. We combine these with Laplace transformation techniques for evaluating the required perturbation expansions. The radii of convergence of these series are determined through a complex variable approach. The method is powerful enough to yield explicit results for such systems as the two species Volterra model, Goodwin's model of protein synthesis etc. and as an illustration, is applied here to Cowan's model of neuroelectric activity. We also point out the usefulness of the action integral in the case where parameters occurring in dynamics have slow time variations.

Notes: Abstract Zwanzig and Mori's projection-operator method is used in order to derive a generalized nonlinear Fokker-Planck equation for one “relevant” species in the many species conservative Volterra model. The deterministic, autonomous, Markovian equations of motion, when averaged over a suitable ensemble of initial conditions in general give rise to a non-autonomous, non-Markovian stochastic process for the evolution of this relevant species. Moreover, this relevant species may show irreversible damping, although self-interaction terms are absent in the many species model.

Notes: Abstract The hydrodynamical problem of flow in proximal renal tubule is investigated by considering axisymmetric flow of a viscous, incompressible fluid through a long narrow tube of varying cross-section with reabsorption at the wall. Two cases for reabsorption have been studied (i) when the bulk flow,Q, decays exponentially with the axial distancex, and (ii) whenQ is an arbitrary function ofx such thatQ-Q 0 can be expressed as a Fourier integral (whereQ 0 is the flux atx=0). The analytic expressions for flow variables have been obtained by applying perturbation method in terms of wall parameter ε. The effects of ε on pressure drop across the tube, radial velocity and wall shear have been studied in the case of exponentially decaying bulk flow and it has been found that the results are in agreement with the existing ones for the renal tubules.

Notes: Abstract Voltage clamp experiments, which determine the kinetic parameters of calcium conductance of cardiac muscle, (d ∞,f ∞, τ d and τ f ) are analyzed with a generally accepted expression for slow inward currentI s=g sdf (E-E R). Activation (d) and inactivation (f) reach the final valuesd ∞ andf ∞ with time constants τ d and τ f respectively. The analysis indicates that the measuredf ∞ agrees with the theoreticalf ∞, but the measuredd ∞ differs from the theoreticald ∞ by a factor which depends on τ d . The peak tension can be made to correlate closely with the theoreticald ∞ after a correction factor is applied to the raw measurements of activation. It can be shown that experiments designed to measure τ f can also be used to determine τ d with greater accuracy.

Notes: Abstract A deterministic model for an SIR epidemic with silent infections is investigated. It is shown for the model studied that the extent to which silent infections are present may be determined from data concerning only those individuals with symptomatic infection.

Notes: Abstract An urn contains balls of different colors. Specified numbers of each color are added and form a reinforcement. The total reinforcement is randomly removed, forming a depletion. The process, not necessarily with the same reinforcements, is performed a number of times. The factorial moment generating function of the urn configurations at any stage is given in terms of multivariate difference operators. Cases when the reinforcement vector is defined as a stochastic variable are considered. The problem is a generalization of an urn model associated with radioactive atoms and stable atoms proposed by S. R. Bernard. The solutions given here have a definite application to the problem of modelling tracers in compartmental systems.

Notes: Abstract A multi-compartmental model with particles producing offspring according to the Markov branching process has been studied. Explicit results are given for the two-compartmental system and for irreversible general multicompartmental systems. The known models in stochastic compartmental analysis are shown to be particular cases of this model and applications are cited.

Notes: Abstract The irreversible Michaelis-Menten scheme may be reduced to a pair of autonomous first-order differential equations. The phase-plane behaviour of these is investigated.

Notes: Abstract A discrete one-dimensional model of convection-diffusion in branching alveolar ducts is described and it is shown that, for a suitable choice of effective axial dispersion, the solution closely approximates that for an axially symmetric representation, at least for Peclet numbers Pe〈1. Following earlier work a composite model of a uniform lung is formed by matching such a respiratory pathway (now having the more convenient one-dimensional form) onto a trumpet representation of the conducting airways. Enhanced mixing due to heart action, and isotropic volume changes of trumpet (in addition to the pathway) during breathing are additional factors included. Calculations are made of O2 concentrations during steady-state breathing and of the concentration of inert gas during single breath wash-out of a gas mixture containing it. Predicted alveolar levels in each case agree extremely well with published data, although no alveolar slope is obtained for the inert gas.

Notes: Abstract The evolution and local stability of a system of two interacting species in a finite two-dimensional habitat is investigated by taking into account the effects of self- and cross-dispersion and convection of the species. In absence of cross-dispersion, an equilibrium state which is stable without dispersion is always stable with dispersion provided that the dispersion coefficients of the two species are equal. However, when the dispersion coefficients of the two species are different, the possibility of self-dispersive instability arises. It is also pointed out that the cross-dispersion of species may lead to stability or instability depending upon the nature and the magnitude of the cross-dispersive interactions in comparison to the self-dispersive interactions. The self-convective movement of species increases the stability of the equilibrium state and can stabilize an otherwise unstable equilibrium state. The effect of cross-convection (in absence of self-dispersion and self-convection) is to stabilize the equilibrium state in a prey-predator model with positive cross-dispersion coefficients for the prey species. Finally, it is shown that if the system is stable under homogeneous boundary conditions it remains so under non-homogeneous boundary conditions.

Notes: Abstract Various observations suggest that sympathetic ganglia act as local integrative centers redistributing preganglionic excitation (i.e. the information issued by the central nervous system) to the postganglionic fibers (and effector organs). In order to support this concept a simple mathematical model of the elementary integration process, treating the case of a single preganglionic compound action potential, has been developed. This quantitative description, based on a few elementary assumptions, shows a possible way of processing preganglionic excitation in the ganglion. It is shown that on a particular ganglion cell the probability distribution of the number of activated synapses obeys hypergeometric distribution and hence, the postganglionic compound action potential is built up of several compound action potentials occurring at different times. The former correspond to different groups of firing cells. The model discloses modes of structural and temporal pattern generation performed by the sympathetic ganglion.

Notes: Abstract Amino acid sequences have already been examined in some detail in order to relate them to structural aspects, homology and gene duplication. This report introduces the concept of internal uniqueness of tripeptides within protein sequences and uses the Monte Carlo method to study this property. Some idea of internal uniqueness may be obtained from such an analysis using only a single sequence if the probability of the random occurrence is about 0.001 or less. This method of analysis is similar to that used in quantitative evaluations of homology. When the probability of the random occurrence is larger than 0.001 a homologous group of sequences is required and the random probabilities may be compared with the real occurrences within the group. From such an examination insulin and cytochrome c are identified as protein sequences with high internal uniqueness. A comparison of data from internal uniqueness and gene duplication analyses shows that these two properties need not be related. Results of the analysis point to internal uniqueness as an additional parameter for inclusion in speculations on why twenty amino acids are coded in protein structure.

Notes: Abstract The concept of a tolerance net formalises simultaneously spatial closeness and nearness of neuronal activity. A method of constructing tolerance group nets is presented, leading to a means of construction of all very homogenous tolerance nets as group quotients. The dihedral group of order eight is taken as an illustrative example.

Notes: Abstract Two methods are described for calculating the value of the exponentx in the equation flow =k×diameter x , as pertaining to a branch of the bronchial tree. In the lungs from three humans, two dogs, one hamster, and one rat mean values ofx between 2.419 and 2.903 were found. They lie within the range of 2.333 to 3.0 predicted by the analysis of Uylings (Bull. Math. Biol. 39, 501–519, 1977).

Notes: Abstract A stochastic analysis of a nonlinear selection model is presented. The model, based on Eigen and Schuster's theory of selection and evolution of biological macromolecules, considers the effects of fluctuations on the individual concentrations of macromolecules as well as the total population numbers in constrained systems. Our analysis shows that one of the models most often treated deterministically (referred to as constant organization in the literature) becomes unstable when fluctuations in the total population number are considered. An alternative model which apparently has built in self-regulating properties is analyzed and proves to be stable except for some special cases of degeneracy.

Notes: Abstract Previous compartmental models have introduced variability either at the particle or at the replicate level. This paper integrates both types of variability through the concept of clustering. The paper develops two different, general clustered models, each with time-dependent hazard rates for the clusters and for the particles within the clusters, and each with random initial number and sizes of clusters. The coefficient of variation of the total number of particles,CV[X(t)], for either model is shown to be bounded below, under very broad conditions, by the coefficient of variation of the initial number of clusters,CV[c(0)]. This high relative variability of the clustered models makes them potentially very useful in kinetic modeling. In many applications, binding and clustering are common phenomena, and two applications of the models to such phenomena are breifly outlined.

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.