ALBERT

Notes: Abstract Biological tree-like structures, such as mammalian tracheobronchial airways, are complicated branching systems. One problem in modeling such systems is the reassignment of the number of segments at a given generation in the model being constructed. A hypothesis is proposed which has successfully been used in modeling mammalian lung airways.

Notes: Abstract This paper is a sequel to a paper by the author entitled “Restricted Transition Probabilities and Their Applications to Some Problems in the Dynamics of Biological Populations” (Bull. Math. Biophysics, 1966,28, 315–331). The paper is divided into two parts. In part one some aspects of the maximum size attained by the population during a finite time interval are studied for the case the stochastic process underlying the evolution of the population is a birth process. Two interesting by-products emerge from the study presented in part one; namely a combinatorial method of finding solutions to the Kolmogorov differential equations in special cases, and secondly, a set of criteria for the optimum allocation of genotypes in the host population of a host-pathogen system. The optimum allocation of genotypes in the host population is a problem of practical importance in controlling plant pathogens. In part two the theory of restricted transition probabilities developed in the companion paper is applied in finding the distribution of the time to the appearance of the first mutation for the case of a two dimensional birth process. The distribution of the time to the appearance of the first mutation is of importance in understanding the role mutation plays in the evolution of a population, particularly in the pathogen population of a host-pathogen system.

Notes: Abstract The complex arrangement of the muscle fibers in the ventricular wall and the nonsymmetric contraction and expansion of the ventricle preclude the writing of a differential equation of motion for the ventricle as a whole. We can, however, describe the motion of the ventricle by describing the motion of the dimensional parameters length and diameter; the radius, circumference, cross-sectional area, and volume following naturally from these. The ventricle is assumed to be an ellipsoid of revolution and the dimensional parameters to be periodic functions of time. Each of the parameters is expressed as a Fourier series.

Notes: Abstract Le modèle de Nelson peut-être considéré comme une approximation du modèle de Hodgkin-Huxley. Moins précis, il est plus maniable. Le modèle de Nelson peut également être considéré comme une généralisation du modèle de Hodgkin-Huxley. En effet, il introduit des liaisons synaptiques localisées ou diffusantes, et un processus de facilitation. Le mécanisme des liaisons synaptiques ne se traduit pas facilement dans le langage de Hodgkin-Huxley. Par contre, le processus de facilitation s'interprète facilement. Nelson's model can be taken as an approximation of Hodgkin-Huxley's model. Its precision is lesser, but it is more usable. Nelson's model can also be taken as a generalization of Hodgkin-Huxley's one; for it introduces localized or diffusing synaptic connexions and a facilitating process. The mechanism of synaptic connexions cannot be easily translated into Hodgkin-Huxley's language. On the contrary, the facilitating process is easily interpreted.

Notes: Abstract A spherical model for the human left ventricle with two different types of aneurysms, circular and rectangular-square, is proposed and meaningful relations are derived between the parameters of the aneurysms and ventricle. Such ventricular parameters as stroke volume, end-diastolic volume, and end-systolic volume are given normal human values to compute values for end-systolic radius and percentage shortening of muscle for various sized circular and rectangular-square aneurysms.

Notes: Abstract The Volterra theory of two competing populations is extended to the contemporary social problem of crime control. Domains of stability for the time dependence of the numbers in the criminal and enforcement groups are exposed by a numerical example. Both augmentation and reduction of enforcement can produce a stable system. Average values of the ratio of members in each group show great sensitivity to the control policies adopted by the remaining sector of the total population.

Notes: Abstract The paper deals with interactions of substances via an enzymatic reaction (Bull. Math. Biophysics,25, 141–154, 1963). The substances are the activators, inhibitors and/or substrates of the reaction. Due to the bimolecularity of the processes in the reaction, the quantitative relation between the steady state amount of complexes and the amounts of the substances assumes a typical form. In multiple enzymatic reactions this form is more complicated, though basically similar. Because the substances may influence the steady state amounts of the complexes in opposite directions, the compensation and blocking effects are the properties of enzymatic reactions. The substances with the same direction of influence may potentiate each other. In the enzymatic reaction here considered, the potentiation is always non-negative.

Notes: Abstract Growth-rate functions in analytic form have been obtained for cell cultures in which the doubling times follow the Gaussian and Poisson distributions. The growth-rate functions are calculated by using Laplace transforms to solve an integral equation previously presented. Oscillatory solutions result if a substantial fraction of the cells in a culture are synchronized to divide at some particular time. The synchrony and, hence, the oscillatory character of the growth-rate function eventually disappear because of the non-zero variance of the doubling-time distribution. If their variances are sufficiently small, the Gaussian and Poisson doubling-time distributions lead to growth-rate functions that become identical in the limit of large time.

Notes: Abstract IfN(t) is the expected number of cells in a culture at timet, $$\dot N(t)$$ the corresponding time derivative, andf(t−τ)dt the probability that a cell of aget−τ at timet will divide in the succeeding time intervaldt, then according to Hirsch and Engelberg (this issue) there obtains the integral equation $$\dot N(t) = 2\int_{ - \infty }^t {f(t - \tau )\dot N(\tau )d\tau }$$ for describing the dynamics of the cell population. It is the purpose of this note to give two alternative derivations of this equation, one based on the age density equation of Von Foerster, and the other based on a generalized form of the Harris-Bellman equation describing the first moment of an age dependent, branching process. In addition, a probability model is posed from which the Von Foerster equation and, hence, the Hirsch-Engelberg equation readily follows.

Notes: Abstract A model of the dissolution kinetics of powdered enamel is developed based on the kinetic rate termq, $$q = K'H - k'C \cdot P_1$$ , whereH=[H +],C=[Ca ++] andP 1=[HPO 4 = ]. The differential equations describing the rate of mineral dissolution (and the linearly related rate of appearance of calcium and phosphate in solution) have been derived and solved for three basic cases: (1) when thepH of the solution and surface area of the enamel are considered constant, (2) when thepH is assumed constant, but the reduction in surface area during dissolution is considered, and (3) when the rise ofpH resulting from the buffering effect of the dissolved enamel is considered along with the change in surface area. Analytical solutions have been obtained for cases (1) and (2), while a numerical solution has been found for case (3). Data from a study on enamel dissolution are presented that agree with the theory of case (3), and it is noted that apH rise as large as 0.5 can occur, as has been shown elsewhere in the literature.

Notes: Abstract On the basis of Landahl's theory of two-choice learning it is shown that application of punishment for wrong responses, without giving award for correct ones, does not lead to complete learning, no matter how many trials are used. If initially a “wrong response” was learned, then an attempt to inhibit it by punishment alone will in a class of cases lead only to a 50% suppression of that wrong response. Possible connection with the problem of effectiveness of punishment as a deterrent for crime is mentioned.

Notes: Abstract A set of characteristic parameters is given for electrophoresis accompanied by diffusion, followed by a method of simplification of the transport equation. The concept of electrophoretic similarity is introduced in connection with the presentation of solutions and the final section contains some dimensional considerations of the potential equation.

Notes: Abstract We show that when we represent (ℓ, ℛ)-systems with fixed genome as automata (sequential machines), we get automata with output-dependent states. This yields a short proof that ((ℓ, ℛ)-systems from a subcategory of automata—and with more homomorphisms than previously exhibited. We show how ((ℓ, ℛ)-systems with variable genetic structure may be represented as automata and use this embedding to set up a larger subcategory of the category of automata. An analogy with dynamical systems is briefly discussed. This paper presents a formal exploration and extension of some of the ideas presented by Rosen (Bull. Math. Biophyss,26, 103–111, 1964;28, 141–148;28 149–151). We refer the reader to these papers, and references cited therein, for a discussion of the relevance of this material to relational biology.

Notes: Abstract A two-dimensional nonlinear integro-differential equation with time-varying coefficients describing the behavior of the fluttering wing-body systems typical of natural flight mechanisms has been deduced from the Navier-Stokes equation which generalizes local pressure and velocity distributions in the externally oscillating air field. The resulting equation for the wing forces is combined with an analogous expression for the forces of gravitation and acceleration associated with the body. The air acceleration force, not previously considered in bio-physical models of insect and bird flight, is shown to arise from a formal analysis of unsteady or time-varying contributions to the velocity field, while the square form of the conventional steady state aerodynamic forces is derived from the intertial terms in the Navier-Stokes equation with the aid of the approximations of Newtonian impact theory. Previous calculations (Houghton, 1964) have indicated that the contribution to gravitational stability of air acceleration and aerodynamic life are roughly in the ratio of 3:1.

Notes: Abstract Certain types of cortical electrical events are non-propagated so that the associated electric fields must have standing wave characteristics. However, cortical electric events typically are generated by neurone populations which cannot be activated simultaneously on impulse driving. Hence the sum of the standing wave fields due to asynchronous activation of adjoining regions of cortical neurones must give the appearance of a traveling wave. Analysis of cortical waveforms is further complicated by curvature in cortical surfaces. A model is presented that shows the effects of curvature and time lag in activation on the form of the potential at points in space around a laminar array of elements simulating a population of cortical neurones. The results are compared with waveforms evoked by single-shock stimulation of the prepyriform cortex in cats.

Notes: Abstract A continuity equation for cell-number density in a population of cells is derived, and a system of equations for eliminating parameters between the general solution and the initial distribution obtained.

Notes: Abstract Use of an electrical model of the left ventricle of the heart and the arterial system permits analysis of the changes which take place as the capacity of the myocardium for generation of force decreases. The model is simple in structure, and its construction and practical testing would not be difficult. It demonstrates that, as the heart muscle weakens, the peak of intracardiac force occurs later in systole, and the difference between the intracardiac pressure and the aortic pressure in the second half of systole is much greater than for the normal heart. The feedback mechanisms which are proposed to affect myocardial contractility would affect this compensation for cardiac weakening. Indices to categorize the behavior of the normal, compensated though weakened, and decompensated myocardium are proposed.

Notes: Abstract The frequency distribution in a population of cells of the quantityCD (defined as the amount of some chromosomal substance in a cell which containsC chromosomes) is calculated using the distribution in the population of the amount per chromosome,D, and the distribution of chromosome number,C.

Notes: Abstract The rate of chromosomal DNA synthesis in an exponentially growing population of cells having chromosome-number dispersion is calculated using DNA histogram data, chromosome-number distribution data, and the assumptions that the synthesis rate is constant and DNA double exactly.

Notes: Abstract An estimate is made of the rate coefficient for linear DNA synthesis with exact doubling in an exponentially growing population of Ehrlich ascites tumor cells having chromosome-number dispersion. Comparison of calculated and experimental results suggest that the assumptions used in the calculation are tenable, but further experimental evidence is needed to prove this.

Notes: Abstract The paper develops further some suggestions made previously (Bulletin of Mathematical Biophysics,28, 283–308, 1966) that certain biological phenomena may be more easily interpreted from a “sociological” point of view by considering the organism as a social aggregate of cells and a cell as a social aggregate of genes. In this light the problems of origin of life on earth, of aging, and of parasitism and symbiosis are discussed. The notion of social aggregates of different orders is introduced.

Notes: Abstract A theoretical and experimental study was made of the mechanical behavior of the cornea. The theoretical analysis included an analytical solution for the symmetrical constraint of a thin, shallow, spherical shell by a rigid indenter. The experimental study investigated the rheology of the cornea with particular emphasis on its compliance with the requirements of the Boltzmann Superposition Principle. Representative results of tests on twenty enucleated hog eyes and two human eyes have been reported. The corneas of the human and hog eyes behaved as linear viscoelastic solids; the human eyes differed from the hog eyes in having a long term creep component. Several eyes were tested at the site of procurement, six to seven minutes after the animal's death, and it was established that creep is not an artifact due to aging or enucleation. The analytical and experimental results were combined to study some instruments used to detect the level of pressure in the eye. The theoretical analysis predicted that a type of elastic instability occurs during the process of flattening a small portion of the cornea; this is discussed with reference to the Goldmann and Mackay-Marg tonometers. The role of corneal creep was considered with reference to the response of the Schiøtz indentation tonometer during the time dependent process known as tonography.

Notes: Abstract Following previous studies, differential equations are established which determine the variation of the stimulus towards a corrective turn of the steering wheel and its effect on the excitation of the centers in the brain which results in the production of the corrective turn. The equations are derived under the highly oversimplified assumption that all excitation thresholds are so small that they can be neglected. Under these assumptions it is found that the tracking curve of a car is a sinusoid with negative damping, that is, with an ever increasing amplitude. Driving under these assumptions is imposible since the car will always eventually jump off the road. The possible effects of the threshold as well as stimuli towards corrective turns other than the distance from the edge of the lane are very briefly discussed. In spite of the negative results of the paper, its interest lies in the circumstance that with the complication of the model, we find that driving depends not only on the reaction times as the only “purely biological” parameter, but on three other neurobiophysical constants. In a subsequent paper (Rashevsky, 1967) it is shown how the introduction of one or more purely biological parameters of the driver makes a stable driving regime possible.

Notes: Abstract A vast number of biologically important processes are based upon bimolecular systems. In these systems intermediate complexes are formed. Bimolecular systems in which no complex-complex interactions occur are called linear systems of complexes. A definition and some characteristic properties of these systems are given here. There may exist a contradiction of Onsager's principle of detailed balancing in these systems; however, no principal differences are found between the steady state behavior of an open system and that of a closed system. It is shown that the steady state behavior of a linear system of complexes of arbitrary complexity has some similarities with the steady state behavior of a simple bimolecular system, e.g., Michaelis-Menten enzymatic reaction. Multiplicity of action of the substances participating in biomolecular processes may produce some qualitative differences in the steady state behavior of the system.

Notes: Abstract A time-dependent DNA histogram is calculated for an irradiated population of cells under the limiting assumption that the cells cannot pass through prophase due to the effects of the radiation. The population is assumed to increase exponentially prior to irradiation, but after irradiation to neither gain nor lose cells. Chromosome-number dispersion is taken into account in the calculation. The qualitative behavior of the calculated and experimental histograms are in reasonable agreement. The quantitative agreement between the two is relatively good at short post-irradiation times but is poor at long post-irradiation times (say, greater than half the doubling time). This suggests that recovery phenomena cannot be neglected at long post-irradiation times.

Notes: Abstract It is pointed out that the three different stimuli for a corrective turn, namely the distance from the edge of the lane, the rate of approach to the edge, and the angle between the direction of the car and the direction of the lane (Bull. Math. Biophysics,28, 645–654, 1966,29, 181–186, 1967) may act all three simultaneously. It is found that in that case the tracking curve of the car is stable below a critical speed and becomes unstable above it.

Notes: Abstract Continuing a previous study (Bull. Math. Biophysics, 28, 645–654, 1966), the biophysical mechanism of a corrective turn is investigated for the case where the stimulus for the corrective turn is produced not only by the perception of the nearness of an edge of the lane, but also by the rate of approach of the car towards the edge. In that case it is found that the tracking curve of the car may consist of a series of damped sinusoids and safe driving would be possible at any speed if it were not for the endogenous fluctuation in the driver's central nervous system. If the effect of the rate of approach increases sufficiently rapidly as the distance to the edge of the lane decreases, then a stable undamped oscillating tracking curve is possible. The case is also studied where the driver makes a corrective turn in response to a direct perception of the angle between the direction of the lane and the longitudinal axis of the car.

Notes: Abstract The principle of minimal work requires that the conducting airways of the human lung should have a maximum radius for minimal resistance to gas flow. At the same time there is a requirement that the airways should have a minimal volume for economy of space. These two opposing requirements have been investigated mathematically, and a method for calculating the angle of branching which produces minimal volume has been derived. The relationship of the radii of the parent and daughter branches to produce minimal resistance has been similarly defined. By measurement of a bronchial cast from a human lung the extent to which the predicted optimum structure is realized in practice has been shown. The change in structure associated with change of function at the transition from conducting airway to diffusion zone has been demonstrated.

Notes: Abstract This paper considers a class of set-theoretical entities, calledn-rank Linnaean structures, which are intended as abstract models of the taxonomic classificatory systems of biology. In the first part, devoted to formalism, finite Linnaean structures are discussed in complete generality; but, in addition, eight distinct subclasses are noted and some of the properties of their elements are explored. In the second part, concerned with applications, it is shown that taxonomic systems may be recast in the form of finite Linnaean structures, and an effort is made to show that some undesirable features of earlier models are avoided without artificiality and without abandoning extensional mathematics.

Notes: Abstract Generalizing an idea of M. Richardson (Fundamentals of Mathematics, New York: Macmillan Co., 1958), an APS on a given populationP is a non-empty collection of non-empty subsets ofP such that ifA is in the collection andA⊆B, thenB is in the collection. From a structure of this kind a partial ordering ofP, called therelated bumping order, is derived. The question is raised as to what kinds of partial orderings can be so obtained. For structures determined by voting weights of the members of the population, a complete characterization of all possible bumping orders is obtained.

Notes: Abstract The “second method” of Liapunov is used to perform a stability analysis of a mathematical model of the neuron. This analysis is based on the hypothesis that the firing of the neuron coincides with a temporary state of instability of the system, and that the initiation of all-or-none process depends on the magnitude of membrane depolarization and its first time derivative. It is found that the stability (and hence the possibility of a second firing) is restored approximately when the rate of membrane repolarization is at a maximum. This result predicts that the duration of the period of absolute refractoriness in neurons would be about 75 per cent of the spike duration, and thus shorter than the value usually obtained from experimental measurements.

Notes: Abstract Some aspects of masking phenomena are considered in terms of the simplest possible model of two-factor neural elements. The effect of a number of variables can be accounted for, but the introduction of an internuncial element results in a masking function which need not be symmetric about zero delay interval. As an illustration, the results for a special case are compared with available data. In general, such a model results in a masking function which depends on the intensity, area, and duration of the stimuli, as well as on the temporal and spatial separation between them.

Notes: Abstract The general equations are discussed describing two species in competition or in symbiosis or feeding one on the other.

Notes: Abstract It is shown that the principle of biological epimorphism (Rashevsky,Mathematical Principles in Biology and Their Applications, Springfield, Ill.: Charles Thomas, 1960) is contained in the theory of organismic sets (Bull. Math. Biophysics,29, 139–152, 1967) if an additional postulate not directly connected to mappings is made.

Notes: Abstract This paper deals with bimolecular systems in which also complex-complex interactions occur. Because of the complexity of the problem, an approximation in a form of coupled linear systems of complexes (Bull. Math. Biophysics,29, 1–16, 1967) is considered. Two types of couplings, serial and parallel, are studied. In the serial coupling the nonlinear system of complexes has the same behavior as its subsystems. An entity, initial sensitivity, has interesting properties: in serial coupling it is at most equal to the product and in parallel coupling, at most equal to the sum of partial initial sensitivities.

Notes: Abstract Amplification effect in the catalytic bimolecular systems is a consequence of the kinetic characteristic of the catalyst. Two types of the coefficient of amplification are defined. The applicability of these definitions is given by the type of the bimolecular system. In a simple example it is shown that the concept of amplification is meaningful in these systems. Furthermore, two rules, analogous to those for a coupling of amplifiers, are derived for the two basic modes of coupling of catalytic systems. Thus, in biological systems the catalytic reactions may be regarded as biologically effective amplifiers.

Notes: Abstract It is postulated that cell hydration is governed by adsorption of water on cell proteins in accord with the Bradley adsorption isotherm, and that the action of a solute in the surrounding solution is to lower the vapor pressure of the solution so that cell water adsorption is decreased by moving down the Bradley isotherm. From these concepts, it is derived that cell volume (V) should be related to solute concentration (x) by the equationV=−E log10 x+F whereE andF are constants which are independent of type of solute. For a non-adsorbed solute this agrees well with experimental data. For solutes which are adsorbed by cell proteins, a correction in the above equation may be necessary at higher solute concentrations, which is shown to be compatible with various experimental data. The types of experiments which are generally used to support the osmotic pressure theory of cell hydration agree equally well with the adsorption theory. The virtue of the adsorption theory is that, unlike the osmotic pressure theory of cell swelling, it is compatible with permeability of the cell membrane to solutes, which has been experimentally observed for various solutes.

Notes: Abstract Conditions under which a time varying electromagnetic field problem (such as arises in electrophysiology, electrocardiography, etc.) can be reduced to the conventional quasistatic problem are summarized. These conditions are discussed for typical physiological parameters.

Notes: Abstract A compartmental lung model with any number of synchronously filling and emptying functional chambers and a common dead space or conducting region is considered. It is shown that the model gives rise to an output, in an open circuit washout determination, which is a weighted sum of exponentials. From estimates of these weights and exponential components, estimates of the model parameters can be recovered. Relations giving the unique correspondence between the output parameters and the model parameters are derived and the existence and uniqueness of solutions established.

Notes: Abstract A physical model that incorporates all the experimental information on the formation of the visual pigment rhodopsin is presented. The visual pigments consist of a chromophore bound to an appropriate protein. Thus rhodopsin (λm 505 mμ) is formed by a Schiff’s base linkage C19H27CH=NH+-opsin (λm 440 mμ) between 11-cis retinal (λm 380 mμ) and the protein opsin (λm 280 mμ). It is found that there exists a red shift in the spectrum of rhodopsin from the Schiff’s base. The model brings an explanation for this red shift. It is shown that such a shift may be due to a charge transfer process (R. S. Mulliken,J. Am. Chem. Soc.,74, 811–824, 1952) between an electron at the double bond of carbons C11−C12 and an atomic orbital of the sulphur present in cysteine. This provides an explanation of the presence of SH-groups in the protein after the absorption of light. A one-electron approximation is used and the dipole momentμ NV ; hence, the oscillator strengthf of the transitionNV is estimated and compared with the experimentally determined extinction coefficient ∈m by mixing 3.5×10−3 M of 11-cis retinal with 8.3×10−5 M of cysteine at pH ranges 6 through 8. Reasonable agreement is found. Solvent, concentration and temperature dependence are shown also.

Notes: Abstract By assigning coordinates to the environmental function space comprising all physical and mental stimuli, mathematical interpretations can be based on such terms as adaptability, and reactivity which relate to individuals interacting with their environment within a society. These psychometric concepts are incorporated into a framework of functional analysis, which permits the optimization of social change by maximizing the satisfaction integral through the use of variational or dynamic programming methods in conjunction with some optimal social policy. The approach provides a mathematical connection between psychology and sociology, and further demonstrates that existing forms of government are simulated by differential equations belonging to the same general class. The synthesis of new classes of functional equations describing social progress is visualized as a legitimate objective for abstract mathematical sociology.

Notes: Abstract The theory of imitative behavior as developed hitherto by the author was based on the assumption that each individual has a natural preference for one of the two mutually exclusive behaviors. The endogenous fluctuations in the central nervous system then result in the individual’s exhibiting the two behaviors alternately with a relative frequency determined by the natural preference. Imitation shifts the natural preference towards one or the other of the two mutually exclusive behaviors. In the present approach it is suggested that the relative frequency of the two mutually exclusive behaviors exhibited alternately is determined by maximizing the “satisfaction function” of the individual, that is by hedonistic factors rather than by purely random fluctuations. Corresponding equations are developed. It is shown that in certain cases, even when the imitation effect is absent, a sort of “pseudoimitation” may occur. Another situation leads, in the case of two individuals only, to a complete “division of labor” between them, with respect to the two behaviors. Each one exhibits only one behavior. After that imitation is introduced explicitly by assuming that imitation by one individual or another increases the satisfaction function of the imitating individual. Results thus obtained show similarities to the results of the old theory.

Notes: Abstract In a previous paper (Bull. Math. Biophysics,29, 565–574, 1967) the author developed equations to represent velocity and hematocrit profiles in quasi-Poiseuille flow of blood. It was assumed that energy dissipation was minimized and that the viscosity depended on hematocrit and shear rate according to the Casson formula. These equations are simplified considerably, placed in a form more suitable for numerical solution and shown to depend on a single dimensionless parameter. Typicalin vivo values for this parameter are calculated.

Notes: Abstract In this paper, discrete models of reproduction are studied. In part one, definitions are given, particularly on order of the reproduction; part two concerns the growth of the population; part three, the phenomena of delay or acceleration; and part four, the consequences of mortality.

Notes: Abstract A model of the regulation of thyroid hormone in the bloodstream of living systems is formulated and analyzed. The portion of this model defined as theregulator includes components representing the thyroid, anterior pituitary and hypothalamic organs and their intercommunicating channels, that is, the peripheral plasma and hypophysial portal circulations and certain neuro-secretory connections. The loss of hormones from the plasma in the living system associated with physiological mechanisms within the peripheral tissue space and the excretory pathways is represented in the model by a lumpedload on the regulator. The model is reduced to a system of differential equations involving eleven parameters and variables, all of which are identified with certain physiological structures and states. Five of these are currently observable by available laboratory techniques and two others are computable explicity from the equations of the model; the remaining four can be computed in the same way to within a multiplicative constant. Procedires for carrying out ten of these measurements and calculations are suggested. On the basis of the equations and parameters of the model, a discussion of the normal behavior and the response of this system to certain types of disturbances is presented. A systematic effort has been made in the development of this model to include all relevant physiological data and relationships reported in the biological literature. A summary of this literature, reflecting the views and interpretations made by the authors of this paper, is included for completeness and ease of reference.

Notes: Abstract In this paper an expression is derived which describes the transient overall uptake of an inert solute by a section of tissue excised with parallel faces and placed upon an impermeable base. The approach diverges from the conventional analyses for perfused tissue (Morales and Smith,Bull. Math. Biophysics,6, 125–141, 1944;7, 47–99, 1945) because the extravascular zone is regarded as a heterogeneous diffusion medium. Account for this is taken by regarding tissue as effectively composed of two phases—a continuous (extracellular) phase similar to water, and a dispersed phase comprising cells of irregular profile. In both phases the relevant mode of uptake is taken as bulk diffusion rather than surface permeation, thus emphasizing the influence of the internal geometry of the tissue upon its overall exchange response.

Notes: Abstract A method for the identification of flow systems by frequency domain analysis has been extended to include systems with recirculation and truncated data curves. Application of the technique to clinical indicator-dilution curves indicates that the method may be useful in the quantitation of intracardiac shunts. A number of numerical examples which demonstrate the accuracy of the method are included.

Notes: Abstract By assigning time-varying coordinates to all environmental stimuli, it has been possible to axiomatize psychoanalytic theory on the five principles of multiple causation, growth-aging influence, genetic influence, historic influence and conscious-unconscious activity. The theorems of summation of response and the inevitability of conscious-unconscious conflict with their corollaries follow directly from the axiomatic foundations, as does the existence of an adaptation-defense mechanism. The interpretation of the defense mechanism in terms of an ego-id feedback system provides the basis for the structural existence of conscious-conscious and unconscious-unconscious conflict.

Notes: Abstract In previous studies of (M,R) (Rosen, 1961; Demetrius, 1966), it was assumed that changes in the structure of (M,R) which were induced by environmental alternations occurred without error. Here, the effect of both “genetic” and “metabolic” malfunctions on the behavior of (M,R) is examined and a subclass of these systems whose behavior is invulnerable to such errors is specified.

Notes: Abstract The definition of an (M,R) is formulated in a way that emphasizes its mathematical properties. Neglecting interactions between the components, it is shown that: (1) An (M,R) contains only one non-reestablishable component. (2) If an (M,R) contains only one non-reestablishable component, then that component is central. Examples are given to illustrate the biological significance of these two results. The notion of “lag-independence” is introduced, and it is shown that if a system possesses only one non-reestablishable component which is “lag-independent” then all components are lag-independent. The concepts of reestablishability, centrality and lag-independence are applied in order to suggest various criteria for optimal organization of (M,R).

Notes: Abstract A mathematical representation for the analysis of control mechanisms in biochemical reactions is presented. First, the theoretical concept of concentration in biological systems is developed. Then a system consisting of two functions λ and τ is constructed as a network of single output automata. The range of λ is taken to be formed by a set of twostates qualitatively different from the “repair function” Φ f of a mappingf: A→B in the stimulated Φ1 and unstimulated state Φ0. Likewise, the range of τ is formed by the set δ={f o ,f 1} wheref 1 means the mappingf in its stimulated state andf o in the unstimulated one. It is demonstrated that the mathematical structure described acts as a control mechanism over thef and Φ f , so that two biochemical components,A→B, are transformed at a controlled rate. Some of the biological applications of this model are briefly examined. The Jacob-Monod model, the enzymatic adaptation phenomenon, and the “rheon unit” hypothesis are discussed within our framework. Eventually, a concrete model for the RNA-polymerase mechanism, based on the above discussion, is presented.

Notes: Abstract The application of Rashevsky’s transformationT to a primordial graph yields a set of graphs corresponding to different stages in the development of the organism. However, sinceT is multiple-valued the graphs obtained are not ordered. To obtain an ordering, it is first shown that the set of graphs under consideration is equivalent to a well defined setO (for “organism”) ofn-tuples. A metric is then introduced which is based on a biological consideration discussed by Rashevsky (Bull. Math. Biophysics,16, 317–348, 1954). Since a metric implies an ordering of the setO, with a knowledge of the structure of the primordial, one can obtain the developmental sequence. Unfortunately, at present, the structure of the primordial graph is unknown which makes the direct application of the above principle impossible. Consequently, an indirect approach which makes use of more accessible biological phenomena is discussed as well. The hypothesis thatrate of development decreases exponentially and the implications this has with regard to the metric onO are discussed. It is shown that if the hypothesis is accepted the search for the developmental sequence is narrowed.

Notes: Abstract The chromosomal theory of inbreeding based on a gametic interaction system lead us to define a depression coefficientD. Comparison of random, sib and half-sib matings (with inbreeding coefficientF=0, 1/4 and 1/8) shows thatD depends on the structure of the starting population and on values of the model parameters. This result accounts for responses of lines whose depression does not depend directly on the inbreeding coefficient and which theories of inbreeding based on increasing homozygosity fail to explain.

Notes: Abstract An idealization of chemical combination is formulated as a model of computability, and it is shown that this model has universal computational power just in case assembly has at least two-dimensional space in which to occur. It is also shown that this model, under reinterpretation, corresponds to a cellular automaton in which growth occurs by differentiation only (i.e., the state into which any cell is born is thereadfter fixed). Hence this latter model of growth is also computationally universal.

Notes: Abstract Kinetics of biological light emission processes do not mean what they seem to mean, because measured light intensity is not proportional to reactant concentration but to reaction rate. Therefore, the differential equation for light decay is usually different from that of concentration decay, so that mass action interpretations cannot be applied directly to light intensity decay. An observed second order light decay for Chlorella at 6.5°C, implies Elovich solid state reaction kinetics, which agrees with other evidence for solid state processes in photosynthesis. An observed 1.5 order light decay for Cholorella at 28°C implies second order liquid or solid state reaction kinetics. First ordere light decay implies first order reaction kinetics.

Notes: Abstract Analysis based on the integration of differential inequalities is employed to derive upper and lower bounds on the total populationN(t) = ∫ R θ(x 1,x 2,t) dx 1 dx 2 of a biological species with an area-density distribution function θ=θ(x 1,x 2,t) (≥0) governed by a reaction-diffusion equation of the form ∂θ/∂t =D∇2θ +fθ −gθ n+1 whereD (〉0),n (〉0),f andg are constant parameters, θ=0 at all points on the boundary ∂R of an (arbitrary) two-dimensional regionR, and the initial distribution (θ(x 1,x 2, 0) is such thatN(0) is finite. Forg≥0 withR the entire two-dimensional Euclidean space, a lower bound onN(t) is obtained, showing in particular thatN(∞) is bounded below by a finite positive quantity forf≥0 andn〉1. An upper bound onN(t) is obtained for arbitrary bounded or unbounded)R withn=1,f andg negative, and ∫ R θ(x 1,x 2, 0)2 dx 1 dx 2 sufficiently small in magnitude, implying that the population goes to extinction with increasing values of the time,N(∞)=0. Forg≥0 andR of finite area, the analysis yields upper bounds onN(t), predicting eventual extinction of the population if eitherf≤0 or if the area ofR is less than a certain grouping of the parameters in cases for whichf is positive. These results are directly applicable to biological species with distributions satisfying the Fisher equation in two spatial dimensions and to species governed by certain specialized population models.

Notes: Abstract The equilibrium probability distribution of the process level is studied for a general class of reversible stochastic reactions. A calculationally convenient approximation for equilibrium probabilities is derived and its accuracy is investigated over a range of values of the equilibrium constant. A method of approximating the equilibrium means and variance is developed and illustrated forQ th-order processes.

Notes: Abstract Beside the concept of material inputs and outputs of components of the representation of biological systems given to us by Rosen, the concept of energy is incorporated. The interaction of material and energy is represented by a cartesian product; and separate material and energetical mappings are considered as the new representation of components. These developments generate aMα category, and it is shown thatMα is isomorphic to theM category of previous developments.

Notes: Abstract This paper discusses the solution of a generaln-compartment system with time dependent transition probabilities utilizing the technique described by Cardenas and Matis (1975) (hereafter abbreviated (CM)). In addition, the cumulant generating function is derived for a special class of reversiblen-compartment systems where the time-dependent intensity coefficients corresponding to the migration and death rates are some multiple of each other. The immigration rates can be any integrable function of time. The moments are also obtained and the solution to the two-compartment system is presented explicitly. The solution is illustrated with a linear and a periodic function which forms have been widely reported in the literature.

Notes: Abstract The relations (inflow) = (dose)/(area under indicator curve), and (volume of distribution) = (throughflow) × (mean transit time) are derived by a matrix method for a system of interconnected subsystems, within which spatial indicator activity gradients may exist, and for compartments, within which the indicator activity is spatially uniform. The inflow theorem, is different from the outflow theorem. Equivalent labeling of multi-input systems reduces them formally to single input systems. Foreign indicator flow-volume kinetics are more general than, and include as a special case, tracer flux-mass (metabolic) kinetics. Volume of distribution in the indicator steady state may be different from the equilibrium volume of distribution.

Notes: Abstract Perturbation methods are applied to a differential equation predator-prey model to find the approximate amplitudes and period of limit cycle solutions. In the model the feeding rate per unit predator per unit prey decreases as the prey become scare. The rigorous applicability of the perturbation technique depends on the assumptions that the limit cycle amplitude is relatively small and that near the equilibrium point the growth rate of each species is most sensitive to changes in the density of the other species. The second assumption is usually roughly satisfied in practice and examples are considered which suggest that the first assumption can be greatly relaxed.

Notes: Abstract Signal Detection Theory can be used to provide a mathematical model describing the choice of a predator trying to distinguish between a model and a Batesian mimic. The mathematical model yields a number of a deductions, in particular that it may or may not assist the mimic population if mimics more closely resemble their models. The assumptions underlying the analysis are discussed in some detail.

Notes: Abstract A new type of physical transition, denotedS→S *, has been detected in irradiated organic molecules (λ=546 nm) through their interaction with specific biological macromolecules. In a specific enzyme-substrate interaction, a clear enhancement of the reaction rate is observed, when the substrate is irradiated with sharply well defined times. These “efficient irradiation times” are always of the 5k sec type (k=1, 2, 3, …). They have been consistently revealed in a great number of specific biological interactions. The present note demonstrates an important property, i.e. that forevery irradiation time aS→S * transition is induced in organic molecules. It is shown that for any irradiation times different from the 5k sec type (k=1, 2, 3, …) states of theS * type may occur, but the biological macromolecules may “detect” only theS * states induced by irradiations of the 5k sec type.

Notes: Abstract A semi-empirical model applicable to the flow of blood and other particulate suspensions through narrow tubes has been developed. It envisages a central core of blood surrounded by a wall layer of reduced hematocrit. With the help of this model the wall layer thickness and extent of plug flow may be calculated using pressure drop, flow rate and hematocrit reduction data. It has been found from the available data in the literature that for a given sample of blood the extent of plug flow increases with decreasing tube diameter. Also for a flow through a given tube it increases with hematocrit. The wall layer thickness is found to decrease with increase in blood hematocrit. A comparison between the results of rigid particulate suspensions and blood reveals that the thicker wall layer and smaller plug flow radius in the case of blood may be attributed to the deformability of the erythrocytes.

Notes: Abstract Three dimensional laminar, viscid flow is developed for Newtonian fluids which provides absolute values for axial, radial and tangential velocity fields everywhere if the dimensions of the vessel are known and two simultaneous axial velocities e.g. on and off the central axis in the same plane, and the central axis axial velocity gradient are measured. In addition, normal and shear stresses are determinable. The equation set satisfies geometric and other known flow limiting conditions such as no slip at surfaces etc. and are amenable for inclusion in general, dynamic flow expressions. Alternatively they may be used alone for certain problems involving gradients and secondary flows. A range of illustrations are shown for a distorting vessel with elliptic cross-section and small axial taper (analogous to the pulmonary trunk during ejection).

Notes: Abstract A regulated left ventricular dynamics model is presented which involves interaction of the dynamics of the left ventricular and circulatory systems and their regulation by the central nervous system. On-line human parametric simulation (parameter estimation) and consequential prognostic implications (based on parametric values) are demonstrated. Model responses to simulated physiologic stresses help delineate tolerances of subjects. In order to have an estimate of the reliability of the model, the sensitivity of the model's responses to changes in the values of its intrinsic parameters is assessed. Also determined is the extent to which errors in measuring the pressure affect the calculated values of the model's simulation parameters and subsequently influence the values of other diagnostically useful variables (such as contractility, oxygen consumption rate, heart rate), when the model is used to determine the limiting physiological stress sustainable by the subject. A comparison of the model's composition with those of other similar cardio-circulatory models is included.

Notes: Abstract Then-stage harvesting strategy of Elizarov and Svirezhev is examined. As a result, some important new features appear. A discussion is presented on whether or not one should harvest a species at one time stage or wait until a later time. The paper is concerned with contributions which are primarily mathematical formulations and results for continuous, as well as discrete time, logistic growth of a single species being harvested. Age class structure is ignored.

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 For certain environments, the Darwinian model allows unique prediction of a function that any surviving system adapted to such an environment has to perform. This is the case for those environments that determine a “survival functional” of position in space-time of known shape. Purely temporal survival functionals can be distinguished from spatial and mixed ones. In each case, there exists an optimum path in combined physical and (reduced) metabolic space. Dependent on the admissible error, approximate solutions of different complexity are sufficient. All solutions possess an afferent, a central, and an efferent part. Within this general frame, specific, “probably simplest”, solutions are proposed for adaptive chemotaxis, insect locomotion, lower vertebrates locomotion, higher vertebrates locomotion, chronobiological systems, and immune systems, respectively—or rather, for the underlying functionals.

Notes: Abstract Mathematical models afford a procedure of unifying concepts and hypotheses by expressing quantitative relationships between observables. The model presented indicates the roles of both insulin and glucagon as regulators of blood glucose, albeit in different ranges of the blood glucose concentrations. Insulin secretion is induced during hyperglycemia, while glucagon secretion results during hypoglycemia. These are demonstrated by simulations of a mathematical model conformed to data from the oral glucose tolerance test and the insulin infusion test in normal control subjects and stable and unstable diabetic patients. The model studies suggest the parameters could prove of value in quantifying the diabetic condition by indicating the degree of instability.

Notes: Abstract A model for the dynamics of a single-species population whose birth rate depends on densities of previous generations is introduced. A difference equation formulation is proposed and the solutions classified for the various parameter values. Data from an experimental population of mice growing in limited space is cited and compared with the model predictions.

Notes: Abstract All soft tissues are modeled as either one-dimensionalstrings, two-dimensionalmembranes, or three-dimensionalsolids. Attention is restricted to tissues in which one of the principal stress components is large and positive in comparison with the other negligible components. Results indicate the following: (1) If a deformed string isconstrained to lie on a surface and is free of tangential pressure, the tension is carried by rays which are geodesics of the surface. If a string or membrane isfree to deform in space without normal pressure, the tension rays are straight lines. If a membrane deforms without tangential surface loads, the tension rays are always geodesics on the deformed surface. If a solid deforms without body forces, the tension rays are straight lines. (2) The stress in a string is a constant if the string is free of tangential pressure and has constant cross-sectional area. The stress in flat tension fields free of tangential surface loads decays inversely with distance along a tension ray from the edge of regression. The stress in a spherically symmetric tension field free of body forces decays inversely with the square of the distance from the center of the sphere. (3) Stress singularities can occur in soft tissues, such as at the corners of a closed rectangular hole in a flat membrane strip. (4) The tension rays in the torsion of soft annular membranes are more steeply inclined from the radial direction than the tension rays for hard metals equally displaced.