ALBERT

Notes: Abstract The theory of complementary variational principles is used to obtain maximum and minimum principles for diffusion problems with Michaelis-Menten kinetics. In an illustrative calculation we obtain an extremely accurate variational solution in good agreement with the numerical solution of McElwain (1978).

Notes: Abstract Necessary and sufficient conditions for primitivity of a product of two Leslie matrices are given. Such a product could be used in modeling the growth of a population governed alternately by two different sets of fertility and survival parameters.

Notes: Abstract Zadeh's transfer function method for linear time-variable systems is used to apply frequency-domain analysis to a periodically time-varying elastance model of the left ventricle. Left ventricular pressure computed from the system function of the time-varying elastance and the phasors of aortic flow shows a typical waveform of the measured ventricular pressure.

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 A theory of noise fluctuations is developed which is applicable to systems of any size in which unimolecular or bimolecular reactions are occurring. The main difference between small and large reacting systems is that in the former the probability of finding a particle in a particular state does not obey a Gaussian distribution, but satisfies a distribution which reflects the mechanism of the chemical reaction. This difference is reflected in the main result of the theory: an autocorrelation function that is expressible as a sum of exponentials, the amplitudes of which are explicit functions of the moments of the distribution. Thus, by using small systems, the autocorrelation function,in principle, allows the elucidation of reaction mechanisms. Numerical simulations indicate that for reacting systems having ten or fewer particles, the deviation of the autocorrelation function from a single exponential should be easily detectable, and that estimates of the first four moments of the distribution should be possible. Accurate inference of the distribution, however, will require further mathematical and experimental advances.

Notes: Abstract The recent mathematical formalization of the concepts of matter and extrinsical energy, which are used for the relational representation of biological systems, is employed in the analysis of the important experimental discoveries of Comorosanet al. related to low energy electromagnetic irradiations on enzyme substrates. By means of the present analysis one of the properties inherent to the experimental phenomena is more precisely exposed, and theoretical developments corresponding to “energetical evolutions” in a biological system (Leguizamón, 1976) may now have an experimental basis. Important limitations are introduced for the validity of the commutativity and associativity of cartesian product of sets, when they represent matter and its linked extrinsical energy. In connection with this last aspect, new important knowledge is obtained for the relational mathematical representation of biological systems.

Notes: Abstract The structure of solutions to a simple spatially dependent population model involving growth and death is investigated. Two forms of motility of the population are considered: (1) random motion only modeled by a Fickian law, and (2) a directed component of motion (chemotaxis), included in addition to the random motion. Under certain growth conditions a traveling wave of constant speed is approached. This speed can be increased by the addition of the chemotaxis with a corresponding increase in the asymptotic population. Development of initial conditions into a wave is illustrated numerically.

Notes: Abstract This paper describes mechanisms of intracellular and intercellular adaptation that are due to spatial or temporal factors. The spatial mechanisms support self-regulating pattern formation that is capable of directing self-organization in a large class of systems, including examples of directed intercellular growth, transmitter production, and intracellular conductance changes. A balance between intracellular flows and counterflows causes adaptation. This balance can be shifted by environmental inputs. The decrease in Ca2+-modulated outward K+ conductance in certain molluscan nerve cells is a likely example. Examples wherein Ca2+ acts as a second messenger that shunts receptor sensitivity can also be discussed from this perspective. The systems differ in basic ways from recent diffusion models. Chemical transducers driven by membrane-bound intracellular signals can establish long-range intercellular interactions that compensate for variable intercellular distances and are invariant under developmental size changes; diffusional signals do not. The intracellular adaptational mechanisms are formally analogous to intercellular mechanisms that include cellular properties which are omitted in recent reaction-diffusion models of pattern formation. The cellular models use these properties to compute size-invariant properties despite wide variations in their intercellular signals. Mechanisms of temporal adaptation can be derived from the simplest laws of chemical transduction by using a correspondence principle. These mechanisms lead to such properties of intercellular signals as transient overshoot, antagonistic rebound, and an inverted U in sensitivity as intracellular signals or adaptation levels shift. Such effects are implicated in studies of behavioral, reinforcement, motor control, and cognitive coding.

Notes: Abstract Large radiation doses to the lung can cause early death from cardiopulmonary insufficiency resulting from radiation pneumonitis and pulmonary fibrosis. A model for early death following inhalation of insoluble radioactive particles is propose. The model is based on three assumptions: (1) early death results from damage to a cluster of cells from a large number of cell clusters at risk, (2) the dose that causes early death depends on how the radiation is delivered in time and (3) the cell clusters at risk to damage are equally sensitive ro radiation. Results from asymptotic theory of extreme values, along with biophysical considerations, suggest that the cumultive distribution function for the absorbed radiation dose to the production of pulmonary injury sufficient to cause early death is best estimated by the third asymptotic distribution without a threshold. This distribution function is identical to the Weibull cumulative distribution function. Data for Beagle dogs after inhaling relatively insoluble forms of alpha- or beta-gamma-emitting particles are shown to support the Weibull model.

Notes: Abstract Models of the human respiratory tract were developed based on detailed morphometric measurements of a silicone rubber cast of the human tracheobronchial airways. Emphasis was placed on the “Typical Path Lung Model” which used one typical pathway to represent a portion of the lung, such as a lobe, or to represent the whole lung. The models contain geometrical parameters, including airway segment diameters, lengths, branching angles and angles of inclination to gravity, which are needed for estimating inhaled particle deposition. Aerosol depositions for various breathing patterns and particle sizes were calculated using these lung models and the modified Findeisen-Landahl computational scheme. The results agree reasonably well with recent experimental data. Regional deposition, including lobar deposition fractions, are also calculated and compared with results based on the ICRP lung deposition model.

Notes: Abstract The completely symmetrical system is defined as having identical transfer coefficients between pairs of compartments and the same loss coefficient for each compartment. The eigenvalues and eigenvector are explicitly found along with the inverses of the system matrix and the matrix of eigenvectors. Many properties, special instances of more general theorems, can be seen at once from the explicit analytic solution of the initial value, washout and washin problems. The system serves as a known case for testing estimation procedures, algorithms for solutions of linear systems, eigenvalue-eigenvector and inversion routines and is of considerable tutorial value.

Notes: Abstract The mathematical structures underlying the theories of organismic sets, (M, R)-systems and molecular sets are shown to be transformed naturally within the theory of categories and functors. Their natural transformations allow the comparison of distinct entities, as well as the modelling of dynamics in “organismic” structures.

Notes: Abstract To explain the sodium conductance change using Wei's dipole model (Wei, 1969), we may expect that during depolarization the dipole's population difference, ΔN, is first reduced and then returns more slowly to its resting value. This paper shows that the experimental results of gating currents support this idea. Such time course of ΔN, however, is not a usual relaxation process. To account for the unusual behavior of ΔN, we propose two additional assumptions: (1) there exists a special coupling system (probably the intramolecular vibrations) whose coupling strength with the dipoles is much stronger than with the thermal bath (intermolecular vibrations), and (2) there also exist “traps” for the dipole's excitation energy so that this energy is transformed into other energy forms at a rate increasing with the increase of depolarization. Experiments suggest that the traps are proteins located at the inner membrane surface.

Notes: Abstract Current research into the dynamics of iterative ecological and biological models has lead to a number of theorems concerning the existence of various types of iterative dynamical behavior. In particular, much study has been done on the dynamical behavior of the “simplest dynamical system”f b(x)=bx(1−x), which is just the canonical discrete form of logistic growth equations found in ecology, sociobiology, and population biology. In this paper, we make use of some of the techniques and concepts of topological dynamics to construct a number of generalized conjugacy theorems. These theorems are then used to demonstrate that the mappingf b has a number of conjugacy classes in which the dynamics of the iterates is equivalent to within a change of variables. The concepts of fitness and survival in logistic equations are then shown to be independent, if we follow certain intuitive definitions for these concepts. This conclusion follows from a comparison of the conjugacy classes of the functionf b and the extinction sets off b.

Notes: Abstract For chemical reactions not at equilibrium but proceeding in the forward direction in the steady state, a result found by a method first introduced by H. G. Britton (1963, 1965) is generalized to prove that if $${{\vec J} \mathord{\left/ {\vphantom {{\vec J} {\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\leftarrow}$}}{J} }}} \right. \kern-\nulldelimiterspace} {\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\leftarrow}$}}{J} }}$$ is the unidirectional flux ratio, $${{\vec J} \mathord{\left/ {\vphantom {{\vec J} {\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\leftarrow}$}}{J} }}} \right. \kern-\nulldelimiterspace} {\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\leftarrow}$}}{J} }}$$ exp (−ΔG/RT). The conditions under which the equality or inequality applies are discussed. If the unidirectional fluxes are not in the steady state, the unidirectional flux ratio is time invariant in certain specific situations. One such important case is for chemical reaction systems with an ordered sequence of reactions. For systems with more than one pathway, $${{\vec J} \mathord{\left/ {\vphantom {{\vec J} {\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\leftarrow}$}}{J} }}} \right. \kern-\nulldelimiterspace} {\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\leftarrow}$}}{J} }}$$ is not constant except for special cases. These results also apply to diffusional and active transport systems.

Notes: Abstract The nonlinear second-order difference equationx n+1=axn(1-xn−1), where 0≦x nX≦1 anda ≧1, is examined from varying points of view, analytical, numerical and geometrical. An analytic expression is obtained for an invariant attracting curveC ∞ (a) in phase space, which becomes the central object of study. This basic curve, which replaces the simple parabolic shape typical of many analogous first-order models, may have a complicated geometrical structure. As the parametera increases,C ∞(a) undergoes transformations characterized by the dynamical descriptions: stable node→stable focus→stable limit cycle →chaotic attractor. Although the limited characterization ofchaos by the appearance of nonperiodic solutions and solutions of arbitrarily large period is relied upon, this appears to be only a simplified approximation of the real behavior of solutions. Trajectories (x n, xn+1),n=0,1,…, are calculated using the related nonlinear planar mapT a(x,y)=(y,ay(1−x)), and regions of persistence and escape are described for characteristic values ofa. The study of persistence, of even more fundamental interest than the associated problems of periodicity and stability, receives special attention. We introduce a geometrical model, similar in many respects to that for the well-known analoguex n+1=axn(1−x n), but having several new and important features. It appears that as the parametera increases in the chaotic regime there are infinitely many intermittent bursts of increase in the probability that any initial point (x 0, x1) will persist in the unit square under successive iterations of the mappingT a, an unexpected property that should be of interest for applications. A discussion of the applicability of these results to population dynamics theory is given, and it is suggested that such equations might find useful application to problems in developmental biology as well.

Notes: Abstract Based on the principle of minimum power, a mathematical model of the functional state of the oxygen transport system is presented. The optimization model minimizes the power expenditure of the heart, bone marrow, lung and other tissues. The model is used to determine the functional parameters of the oxygen transport system in man under both normal and varying barometric pressures. Theoretical results are compared with experimental data.

Notes: Abstract A quantitative model of ion binding and molecular interactions in the lipid bilayer membrane is proposed and found to be useful in examining the factors underlying such membrane characteristics as shape, sidedness, stability and vesicle size at various cation concentrations. The lipid membrane behaves as a bilayer couple whose preferential radius of curvature depends on the expansion or contraction of one monolayer relative to the other. It is proposed that molecular packing may be altered by electrostatic repulsion of adjacent like-charged phospholipid headgroups, or by bringing two headgroups closer together by divalent cation crossbridging. The surface concentrations of each type of cation-phospholipid complex can be described by simple binding equilibria and the Gouy-Chapman-Stern formulation for the surface potential in a diffuse double layer. The asymmetric distribution of acidic phospholipids in most biological membranes can account for the differential effects of identical ionic environments on either side of the bilayer. The fraction of vesicle material which tends to have a right-side-out orientation may be approximated by a normal distribution about the mean curvature. The theory generates vesicle sidedness distributions that, when fitted to experimental results from human erythrocyte membranes, provide an alternative method of estimating intrinsic cationphospholipid dissociation constants and other molecular parameters of the bilayer. The results also corroborate earlier suggestions that the Gouy-Chapman theory tends to overestimate free counter-ion concentrations at the surface under large surface potentials.

Notes: Abstract The “yellow strips” on the cuticle of the Oriental Hornet (Vespa orientalis, Hymenoptera, Vespinae), present photoelectric properties. A mathematical model for the relative changes in resistance as a photoconductive process conforms to the general model for a semiconductor with traps.

Notes: Abstract Damped nonlinear oscillations in biological and biochemical systems are investigated by the extended Krylov-Bogoliubov-Mitropolskii (KBM) method. A review on the extension made by Popov to the KBM method is given and also further improvements are presented. Applications are made to models of oscillating chemical reactions (Lefever and Nicolis, 1971), FitzHugh (1961) equations, and population dynamics (Gatto and Rinaldi, 1977). Comparison to damped oscillating physical and engineering systems is made.

Notes: Abstract The conditions that will allow the lumping together of several age classes in the Leslie model are investigated. We show that if the lumping is to be valid for all population distributions, then the parameters of the model must be periodic. Lumping is valid when the population is in equilibrium, but equilibrium should be tested before the model is lumped.

Notes: Abstract Catastrophe theory is a mathematical theory which, allied with a new and controversial methodology, has claimed wide application, particularly in the biological and the social sciences. These claims have recently been heatedly opposed. This article describes the debate and assesses the merits of the different arguments advanced.

Notes: Abstract Estimates of capillary tracer permeability calculated using multiple indicator data depend upon the particular model adopted to describe blood tissue exchange. The model proposed by Crone (1963) is appropriate when some of the injected tracer diffuses into the tissue but does not return appreciably to the bloodstream before data collection is terminated. Under these conditions extraction of tracer by the tissue depends on a single dimensionless parameter, αcap, defined as the ratio of capillary permeability surface area to water flow. The effects of finite red cell tracer permeability on the Crone model estimate of capillary permeability are examined in the present study. The results indicate that even when back diffusion from the extravascular space is negligible, significant errors in the Crone model estimate can be expected when capillary permeability is relatively high and the ratio of red cell to capillary permeability is less than unity. However, when an aliquot of blood is equilibrated with tracer prior to injection and the dimensionless capillary permeability is relatively low (i.e. αcap ≦ 0.25 for a haematocrit≦50%), the whole blood Crone model estimate of αcap will be within 10% of the actual value, irrespective of red cell permeability. Red cell-plasma exchange for commonly used tracer-organ combinations should not significantly affect Crone estimates of capillary permeability under normal physiological conditions, but may be important in low flow situations.

Notes: Abstract Assuming truncated ellipsoidal geometry for the right and left ventricles, a model is developed for the myocardium enabling biventricular mechanical behavior to be studied. Employing pressure-volume data taken from normal dog hearts and from hearts in which the pulmonary artery has been banded over periods of 2–40 weeks, it is shown that: (a) right ventricular wall stresses are higher than left ventricular stresses; (b) right ventricular wall stress increases initially to a maximum after 3–4 weeks followed by a decline to normal and even subnormal levels, attaining a minimum value at 32–33 weeks; (c) left ventricular stresses behave in a similar manner, attaining their maximum and minimum levels after 7–8 weeks and 32–33 weeks respectively. These results suggest that surgical or medical therapy in patients with hypertrophied ventricles might be more appropriate during the period of wall stress reduction.

Notes: Abstract In this paper we describe a mathematical model of the oscillations of the diaphragm which limits the vitreous body from the anterior segment of the human eye after the lens has been removed in a cataract operation. We study the motion of this diaphragm driven by movements of the eye. Firstly, a mathematical statement of the problem is given and then we solve the problem exactly for a given class of eye movements. From the analysis we deduce that significant oscillations of the membrane are driven by saccades and that it is the angular acceleration of the eye which causes these types of oscillations. A numerical example is given.

Notes: Abstract The Lotka-Volterra system of prey-predator equations is considered with a special type of continuous time delay. In the case of equal diffusion coefficients Hopf’s bifurcation technique is used to show the existence of travelling wave train solutions for the prey-predator system.

Notes: Abstract A mathematical model of prothrombin activation is being proposed which includes the feedback mechanism of thrombin and the alteration of factor V by thrombin. This model is in good agreement with experimental data for the dependence of the rate of thrombin formation on the concentrations of factors V and X a . In particular, it correctly predicts the existence and location of a maximum in both of these cases.

Notes: Abstract A new model of the upper tracheobronchial tree is proposed to account for the three-dimensional nature of the airway system. In addition to the tube length, the tube diameter, and the branching angle, the model includes information on the orientation angle of each tube relative to its parent tube. The orientation angle, defined as the angle between two successive bifurcations, is useful for calculating the gravitational inclination of each tube. The information on orientation angle is further used to construct a binary coding system for identifying individual tubes in the airway tree. The proposed model is asymmetrical, but the same principles can be readily used to construct a symmetrical one.

Notes: Abstract In any control system for which the number of independent controls is smaller than the number of degrees of freedom to be controlled, our choice of control in any state is restricted to a submanifold of smaller dimension than the tangent space. This simple fact has a number of important consequences for questions of biological import; we consider its implications for adaptation, for senescent phenomena and for the determination of tertiary structures of polypeptides through control of certain average properties. We also formulate the Pontryagin Maximum Principle of Optimal control theory in such a way as to inquire whether specific biodynamic systems can be regarded as optimal with respect to rate of accumulation of particular quantities of the system. We find that if this is possible, the quantity in question must play the role of a clock.

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 computational model to predict deposition of a wide variety of particulate pollutants in several species of mammals is presented. The model incorporates breathing pattern and detailed anatomical models of the respiratory tract based on extensive morphometric measurements of individual airways. The predicted deposition from this general model is in close agreement with observed deposition of monodisperse aerosols in rats. Particle size and density and respiratory breathing patterns are the critical parameters affecting regional deposition.

Notes: Abstract A theory of antigen-antibody induced particulate aggregation is developed by investigating the stability of model systems of particles. Conditions for the formation of large aggregates are derived by imposing the requirement that at equilibrium a statistically significant number of redundant bonds would occur in a reduced monomer-dimer model system. A relationship is obtained which predicts the fractional agglutination in the reduced dimer system as a function of the antigen, antibody and particulate concentrations: $$\frac{g}{{2f c_0 (1 - g)^{2^ - } }} = \frac{{s_1 }}{r} + \frac{{s_1 s_2 }}{{2!r^2 }} + ... + \frac{{s_1 s_2 ...s_j }}{{j!r^j }},$$ wherec 0 is the initial concentration of monomer,f is a proximity factor,g is the fractional agglutination,s i is the average rate of formation of theith bond from an (i−1)th bound dimer, andr is the average rate of dissociation of a single antibody-antigen bond.

Notes: Abstract The roles of the concentrations of the three interacting constituents in the aggregation process (antibodies, antigens and particulates) are analyzed in detail. It is shown that the basic equation derived in Part I is consistent over a broad range of conditions with experimental findings previously reported.

Notes: Abstract A general mathematical model describing the biochemical interactions of the hormones luteinizing hormone releasing hormone (LHRH), luteinizing hormone (LH) and testosterone (T) in the male is presented. The model structure consists of a negative feedback system of three ordinary differential equations, in which the qualitative behavior is either a stable constant equilibrium solution or oscillatory solutions. A specific realization of the model is used to describe the experimental observations of pulsatile hormone release, its experimental suppression, the onset of puberty, the effects of castration, and several other qualitative and quantitative results. This model is presented as a first step in understanding the physicochemical interactions of the hypothalamic-pituitary-gonadal axis.

Notes: Abstract Based upon the transition rate equation of dipoles in the membrane, we deal with two important aspects of interaction of nerve signals: (1) conditions for nerve excitation and (2) frequency spectrum analysis of nerve impulse. Interrelations between signal amplitudes and frequencies are formulated in detail. There are several important conclusions which can be drawn from our calculations. First, toexcite the nerve, low frequencies are generally more effective than high frequencies. Second, tosedate the nerve (i.e. to suppress undesired activities), high frequencies would suit better. Third, harmonics produced through interactions of nerve signals are not necessarily weaker than the fundamental frequencies. The great significance of our theory is that it indicates in principle the feasibility to alter or rewrite the information contents of a nerve message in our body by applying stimulations of appropriate strengths and frequencies. Thus, the theory provides a physical basis and hence some understanding for a new branch of medicine—neuro therapy such as Nogier's auriculotherapy, Lamy's phonophoresis, Voll's electroacupuncture and the fast rising TENS (transcutaneous electro-neuro stimulation).

Notes: Abstract A new physical property, called resonance of the B-type is hypothetically attached to the λ =546 nm irradiated crystalline (small) molecules. In this respect an up or down configuration is assumed for those states obtained through irradiation times that are multiples of 5 sec. With these assumptions, the cellular receptors that may detect these states appear to belong to three classes: the up, down and alternatively mixed up-down. Using the classic formalism of eigenvectors and eigenvalues, a simple spectroscopic type of formula is derived, through which all the possible states of the above characteristic may be obtained.

Notes: Abstract A model of ecosystems with migration is proposed from the viewpoint of flow. This model explains the following two points: (1) How the density-dependent terms in population dynamics arise as a consequence of migration. (2) How the ecosystem exhibits a hierarchy in energy per unit biomass.

Notes: Abstract For precise boundary conditions of biological relevance, it is proved that the steadily propagating plane-wave solution to the Fisher equation requires the unique (eigenvalue) velocity of advance 2(Df)1/2, whereD is the diffusivity of the mutant species andf is the frequency of selection in favor of the mutant. This rigorous result shows that a so-called “wrong equation”, i.e. one which differs from Fisher's by a term that is seemingly inconsequential for certain initial conditions, cannot be employed readily to obtain approximate solutions to Fisher's, for the two equations will often have qualitatively different manifolds of exact solutions. It is noted that the Fisher equation itself may be inappropriate in certain biological contexts owing to the manifest instability of the lowerconcentration uniform equilibrium state (UES). Depicting the persistence of a mutantdeficient spatial pocket, an exact steady-state solution to the Fisher equation is presented. As an alternative and perhaps more faithful model equation for the propagation of certain species properties through a homogeneous population, we consider a reaction-diffusion equation that features a cubic-polynomial rate expression in the species concentration, with two stable UES and one intermediate unstable UES. This equation admits a remarkably simple exact analytical solution to the steadily propagating plane-wave eigenvalue problem. In the latter solution, the sign of the eigenvelocity is such that the wave propagates to yield the “preferred” stable UES (namely, the one further removed from the unstable intermediate UES) at all spatial points ast→∞. The cubic-polynomial equation also admits an exact steady-state solution for a mutant-deficient or mutant-isolated spatial pocket. Finally, the perpetuating growth of a mutant population from an arbitrary localized initial distribution, a mathematical problem analogous to that for ignition in laminar flame theory, is studied by applying differential inequality analysis, and rigorous sufficient conditions for extinction are derived here.

Notes: Abstract The binding of mono-, di- and trivalent cations to negatively charged surfaces is studied within the framework of a modified Gouy-Chapman equation. For any given combination of ions of the above valences, the existence and uniqueness of the solution for the surface potential is shown. The treatment provides the surface potential and charge density. For a system containing only monovalent and divalent ions, analytical solutions are given. When trivalent ions are also present, a procedure based on numerical integration is described. The distance dependence of the electrostatic potential for planar surfaces is given. The calculations provide the amount of cations tightly bound and the amount trapped in the double layer region. The competition between cations for binding to surfaces is elucidated.

Notes: Abstract This paper deals with a model describing the behavior of barium-treatedApalysia neurons. The model is represented by a dynamical system, so-called “complete system”, defined in R4 and depending on a small parameter. The study of this system under zero membrane current conditions was performed with the use of the qualitative theory of singular perturbations. We show that this system has a stable periodic solution of the discontinuous type when the small parameter tends to 0+. A reduced system defined in R3, associated to the complete system was also studied: it corresponds to a constant activation of the inward current. We demonstrate that the corresponding hypothetical cell remains silent under zero current conditions.

Notes: Abstract In the study of chemical modification of enzymes and other biologically active proteins, plots of fractional residual activity as a function of number of groups modified per enzyme molecule are often used to establish a correlation between the chemical modification and enzyme inactivation reactions and to determine the stoichiometry of the modification reaction. This paper presents a critical examination of the underlying theoretical framework of such graphs. Whereas these plots are usually presented as linear functions, it is shown here that the general equation describing the relationship between inactivation and modification contains an exponential term; therefore, in the general case, the plot is actually a curve. It is suggested that caution be exercised in the interpretation of such plots and that equations such as those derived in the text be used to fit theoretical curves to the data, in order to maximize the information gained from chemical modification experiments.

Notes: Abstract This communiction argues that so-called “hermaphroditic” tracer systems, which are neither open nor closed, do not exist physically. The argument is based on the assumption that any observable (possibly nonhomogeneous) macroscopic compartment can be approximated by a compartmentC with a finite number of entry points for the tracer, each associated with an abstract subcompartment ofC. It is shown that the “hermaphroditic” property requires that the mean waiting time be infinite in at least one of the subcompartments, or in a subcompartment elsewhere in the system. A subcompartment with infinite mean waiting time must have some sort of memory, of infinite duration, which knows how long a given particle has been retained, however long that is, and thereby determines its probability of departure. Assuming, as seems likely, that no physical basis exists for such an infinite memory, it follows that “hermaphroditic” systems do not exist.

Notes: Abstract The implication of state space structure on the existence of a repeatable experimentE designed to determine if a states∈L has propertyP or notP is investigated. It is shown that if a state spaceL is connected, then no experimentE is repeatable. This formalism is used to demonstrate that if a propertyP has an associated set of points inL which is dense with dense complement inL, then there exists no repeatable experimentE which can be used to test whethers has propertyP or notP. Other consequences of this formalization are discussed.

Notes: Abstract It has been shown that the resistance of flow and the wall shear increase with the size of the stenosis but these increases are comparatively small due to non-Newtonian behaviour of the blood indicating the usefulness of its rheological character in the functioning of the diseased arterial circulation.

Notes: Abstract The modern theory of generalized Hamiltonian systems is used to construct a unified canonical description of the linear Lagrangian biodynamics introduced by Kerner.

Notes: Abstract The purpose of this paper is to justify an asymptotic method developed for the study of peristaltic transport in a tube of arbitrary cross section. Within the framework of long wave approximation, the three-dimensional nonlinear Navier-Stokes equations are reduced to a sequence of two-dimensional linear boundary value problems of Laplace and biharmonic operators. It is shown that, if a Reynolds number is less than some constant, the solution of the approximate equations is indeed an asymptotic approximation to the exact solution of the problem as the ratio of the maximum radius of the tube to the wave length of the peristaltic motion of the wall tends to zero, and the error estimates are expressed inL 2 norms. Furthermore, under the same condition the exact solution is shown to be unique and stable under arbitrary perturbation of spatially periodic disturbance. Application of the stability condition to peristaltic transport in a tube of circular cross section is given.

Notes: Abstract A mathematical analysis, including existence and uniqueness, is given for some boundary value problems which model the flow of a fluid-solute mixture in a tube which is placed in an interstitium. The model permits an interchange of fluid and solute across the tube walls.

Notes: Abstract The vertebrate nervous system has topographic interconnections in many parts, known for example as retinotopy, somatotopy, etc. It is plausible that modifiable synapses play an important role in forming and refining these connections together with the sensory experiences. To elucidate the mechanism of topographic organization, we propose a simple model consisting of two nerve fields connected by modifiable excitatory synapses. The model also includes modifiable inhibitory synapses. The behavior of the model is described by a set of simultaneous non-linear integro-differential equations. By analyzing the equations, we obtain the equilibrium solution of topographic connections. It is also proved that a part of the presynaptic field which is frequently stimulated comes to be mapped on a large area of the postsynaptic field so that it has a good resolution.

Notes: Abstract The steady-state solution of the equations governing substrate exchange between vascular and extravascular compartments separated by a membrane with finite, symmetrical substrate permeability is presented. Substrate removal from the extravascular compartment by Michaelis-Menten saturation type kinetics with negligible diffusion in the axial and instantaneous diffusion in the transverse directions in both compartments are assumed. It is shown that the solution degenerates into known expressions for special linearized and asymptotic cases. The method of solution is also applied to an extension of the original model incorporating autoregulatory feedback effects upon the process responsible for substrate removal.

Notes: Abstract Global asymptotic stability and equilibrium coexistence is established in two species Lotka-Volterra-type competition when there are time delays in interspecific interaction terms and the intraspecies competition is stronger than the interspecies competition.

Notes: Abstract The realized (observed) value of Landau’s dominance hierarchy index is examined. Under a model of constant pairwise dominance probabilities, the observed index is shown to be a strongly consistent estimator of the underlying (true) index. However, a large number of encounters between animals is shown to be required in order to reduce bias and variance to practical levels except when the pairwise dominance probabilities are near one.

Notes: Abstract The steady state potassium conductance as a function of measured membrane potential difference (p.d.) ϕ of the squid giant axon is corrected for the effect of accumulation of potassium in the periaxonal space. This correction is made on the assumption that several mathematical models of the axon are valid. These are (i) the McIlroy (1975), McIlroy-Hahn (1978) model of membrane conductanceg i(i=K, Na) which is a detailed model of passive transport of ions across the axonal membrane with the aid of mobile, negatively-charged carriers, (ii) the Adelmanet al. (1973) compartmental model of the periaxonal and external bathing-solution spaces, (iii) the enzymatic theory of nervous conduction due to McIlroy (1970 a, b, c), (iv) the Wien dissociative effect of the axolemmic electric field on the weak membrane buffer proposed by Bass and Moore (1968) as a trigger mechanism in nervous excitation and (v) the model (McIlroy, 1979) of the interfacial double-layer p.d.s. which are assumed to exist at the membrane’s surfaces because of the presence of a fixed surface charge. From the correctedg k (ϕ) curves the values of the double-layer p.d.s. of model (v) are deduced and these are shown to lead to a consistent, physically reasonable solution for the distance (approx. 6.8Å) between the fixed surface charges and for the dissociation constants of these sites in their interactions with the ions of the extra-membrane electrolytes. Assuming that the selectivity coefficint of the potassium conducting system for the squid giant axon is approx. 52 it is deduced that the potassium permeability,P ks , of the periaxonal barrier ≈1.37(±0.5)×10−4 cm sec−1 and the thickness of the periaxonal space ≈451±159Å.

Notes: Abstract The effects of peripheral layer viscosity on physiological characteristics of blood flow through the artery with mild stenosis have been studied. It has been shown that the resistance to flow and the wall shear decrease as the peripheral layer viscosity decreases.

Notes: Abstract A study was made of Higgins’ model of glycolysis incorporating molecular diffusion of intermediates, utilizing an earlier conjecture due to Landau. Conditions for the existence of asymptotically stable spatio-temporal periodic solutions are obtained.

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.