ALBERT

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 The parameter domain for which the quasi-steady state assumption is valid can be considerably extended merely by a simple change of variable. This is demonstrated for a variety of biologically significant examples taken from enzyme kinetics, immunology and ecology.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Notes: Abstract We have studied an ecological system of two species, which we denotestrong andweak, respectively, that compete for a single food resource. This system is modelled as a three component reaction-diffusion process. In the presence of a solitary pulse of increased resources, the weaker competitor can diffuse toward this surplus, gaining a competitive advantage and hence persisting in contraposition with the classical Lotka-Volterra result. An exact analytical solution has been found through a quantum mechanical analogy. A stability analysis of this solution against changes in different parameters has been carried out.

Notes: Abstract Collapsible-tube flow with self-excited oscillations has been extensively investigated. Though physiologically relevant, forced oscillation coupled with self-excited oscillation has received little attention in this context. Based on an ODE model of collapsible-tube flow, the present study applies modern dynamics methods to investigate numerically the responses of forced oscillation to a limit-cycle oscillation which has topological characteristics discovered in previous unforced experiments. A devil's staircase and period-doubling cascades are presented with forcing frequency and amplitude as control parameters. In both cases, details are provided in a bifurcation diagram. Poincaré sections, a frequency spectrum and the largest Lyapunov exponents verify the existence of chaos in some circumstances. The thin fractal structure found in the strange attractors is believed to be a result of high damping and low stiffness in such systems.

Notes: Abstract In this work, we show that a one-dimensional model of the blood flow across the lungs can reproduce the evolution of a bolus versus the time. Solving the differential equation governing the bolus concentration in the framework of this model, we determine the solution which fulfills Gaussian initial boundary conditions. An effective parameter related to the ratio of a diffusion coefficient to the square of the mean speed of the flow is defined. The determination of its numerical values following a semi-empirical approach enables us to know accurately the mean transit time and the cardiac output. The results have been compared to other methods, and were found in good agreement. Such an approach could be of interest in all studies where the knowledge of flow—including micro-circulation—is needed.

Notes: Abstract How two species interact during and after colonization influences which of them will be present in each stage of succession. In the tolerance model of ecological succession in a patchy environment, empty patches can be colonized by any species, but the ability to tolerate reduced resource levels determines which species will exclude the other. Here, we analyze a meta-population model of the possible roles of competition in colonization and succession, using non-linear Markov chains as a mathematical framework. Different kinds of competition affect the final equilibrial, abundances of the species involved in qualitatively different ways. An explicit criterion is given to determine which interactions have stronger effects on the final equilibrial levels of the weaker, species. Precise conditions are stated for the co-existence of both species. Both species are more likely to co-exist in the presence of an intermediate disturbance frequency.

Notes: Abstract We study boundary value problems for the time-harmonic form of the Maxwell equations, as well as for other related systems of equations, on arbitrary Lipschitz domains in the three-dimensional Euclidean space. The main goal is to develop the corresponding theory for Lp-integrable bounday data for optimal values of p's. We also discuss a number of relevant applications in electromagnetic scattering.

Notes: Abstract Beurling's algebra $A^*=\{f:\sum_{k=0}^{\infty} \sup_{k\le |m|} |\hat f (m)| 〈 \infty \}$ is considered. A* arises quite naturally in problems of summability of the Fourier series at Lebesgue points, whereas Wiener's algebra A of functions with absolutely convergent Fourier series arises when studying the norm convergence of linear means. Certainly, both algebras are used in some other areas. A* has many properties similar to those of A, but there are certain essential distinctions. A* is a regular Banach algebra, its space of maximal ideals coincides with $[-\pi,\pi],$ and its dual space is indicated. Analogs of Herz's and Wiener-Ditkin's theorems hold. Quantitative parameters in an analog of the Beurling-Pollard theorem differ from those for A. Several inclusion results comparing the algebra A* with certain Banach spaces of smooth functions are given. Some special properties of the analogous space for Fourier transforms on the real axis are presented. The paper ends with a summary of some open problems.

Notes: Abstract We construct an algebra of left-invariant pseudodifferential operators on SU(2). We require only that the symbols be homogeneous and C2. For Fourier-bandlimited symbols, we derive the expected formulae for composition and commutators and construct an orthonormal basis of common approximate eigenvectors that could be used to study spectral theory. Some remarks on applications to matrices of operators are made.

Notes: Abstract A neuron is assumed to receive synaptic input of both excitatory and inhibitory natures from a large number of neighboring neurons; it is also assumed that a large number of such impulses are required to raise the neuron’s transmembrane potential to its threshold potential, at which it “fires” or “spikes”. The model is similar to one of Gerstein and Mandelbrot, except that in the absence of input an exponential decay of potential toward a resting level is introduced. Computational methods of determining the spike timeinterval distribution are discussed, along with the inverse problem of estimating the parameters of the system from observed spike time-interval data.

Notes: Abstract Moens-Korteweg relations are developed for the velocity of wave propagation through an orthotropic elastic tube based on the three dimensional equations of elasticity. Numerical examples are presented for the femoral artery of a dog and several other orthotropic materials. These results are compared with those obtained from the equations of motion for an orthotropic elastic medium.

Notes: Abstract Up to the present time, the following property of the product component in the reversible one substrate-one intermediate-one product enzymic mechanism has been taken only as anassumption, viz., during the course of the reaction, the time-rate of change of product concentration is never negative and the product concentration never exceeds its equilibrium value. Applying the methods of the geometric theory of ordinary differential equations it is shown that this result follows as a direct deduction from the differential equations governing the mechanism together with the initial conditions. Further, the nature of the equilibrium point as a stable node is established.

Notes: Abstract An integro-differential equation treatment of multi-compartment systems is developed which permits formal analysis of the incomplete data which is available from partly accessible, partly injectable systems. New transport functions are defined which can be obtained directly from the experimental data. These functions serve to characterize the communication and topology between different accessible compartments and also the reentrant contributions from inaccessible sites. The method gives solutions consistent with those of the differential equation approach when the system is uniformly contiguous and accessible, more complete solutions than those of the integral equation approach when all measured compartments are injectable, and in addition provides complete or partial solutions for certain otherwise analytically intractable systems. Detailed numerical illustrations of the method are given.

Notes: Abstract In continuation of previous work (Rashevsky,Some Medical Aspects of Mathematical Biology, Springfield, Ill.: Charles C. Thomas, 1964, Chap. 23 and Appendix 14), the study of the effects of the physical parameters of the cells of endocrine glands on the onset of sustained periodical oscillations in the interaction between the anterior pituitary and the thyroid hormones is generalized to include the possible effect of the intercellular fluid and of the degree of vascularization. Some conclusions of the previous study remain valid although some modifications must be made. A decreased relative volume of the intercellular fluid and an increased vascularization favor the conditions for sustained oscillations. The permeability of the cells and the permeability of the capillaries appear explicitly in the expressions which show the conditions for sustained periodicities.

Notes: Abstract A general model of a large 2m-ploid breeding population, withr loci ands h alleles at the h th locus is considered. It is assumed that the population is bisexual, non-overlapping and breeds by random mating. The genotypic structure of the population is presented as a bilinear form in the gametic output vectors where the genotype distribution is in the matrix form. Using the concept of the segregation distribution, the genotype proportions in the (n+1)st generation are given. An equilibrium condition for random chromosome segregation is obtained in terms of gene frequencies.

Notes: Abstract The rigidity of the skull and the inertial characteristics and incompressibility of its contents cause the elastic cerebral arteries and veins to act over brief periods of time like rigid tubes of relatively small diameter. Poiseuille's law is applicable to their behavior. The use of this law, in combination with the fact that, during brief intervals, the total volume of the cerebral arteries and veins remains constant, permits derivation of a mathematical expression for the average arterial flow in terms of an average arterial radius. The differentiated equation has five positive roots which represent maxima and minima of the average flow in terms of the average arterial radius. The theoretical results have physiological implications and potential clinical usefulness, which are discussed.

Notes: Abstract Equations are derived for the flow of an anisotropic fluid in a tube. It is argued that these provide a model for arterial blood flow. Particular attention is paid to the effect of radial differences in hematocrit. Sequels to this paper (Bull. Math. Biophysics,29, 565–574; forthcoming, 1967) will respectively demonstrate possible wall-directed forces on the erythrocyte and enlarge on the physiological consequence of hematocrit variations. The present article develops the basic equations and explores the possible role of anisotropic effects in blood flow.

Notes: Abstract Concentric circular lines are present on X-ray diffraction patterns of DNA. They cannot be accounted for by the double helix model. It is suggested that they may indicate the presence of double helical side chains.

Notes: Abstract An extension of an earlier model simulating the effects of light on the drosophila eclosion rhythm is presented. The effects of variable light intensity are described. This allows not only the simulation of certain experiments not covered by the earlier model, but also it permits an extension of the model to other organisms. By changing only its sensitivity to light the model simulates the phase response curves of certain mammals as well as Aschoff’srule.

Notes: Abstract Previous derivations of the Stevens Power Law relating loudness to intensity assume, in addition to level invariance, that the relation must be analytic or at least differentiable. This last condition is replaced here by the weaker one of requiring only continuity.

Notes: Abstract A global kinetic analysis of a model consisting of an autocatalytic zymogen-activation process, in which an irreversible inhibitor competes with the zymogen for the active site of the proteinase, and a monitoring coupled reaction, in which the enzyme acts upon one of its substrates, is presented. This analysis is based on the progress curves of any of the two products released in the monitoring reaction. The general solution is applied to an important particular case in which rapid equilibrium conditions prevail. Finally, we suggest a procedure to predict whether the inhibition or activation route dominates in the steady state of the system. These results generalize our previous analysis of simpler mechanisms.

Notes: Abstract In this work, treating the artery as a thick-walled cylindrical shell made of an incompressible, isotropic and elastic solid, utilizing the large deformation theory and the stress-strain relation proposed by Demiray (1976b,Trans. ASME Ser. E, J. Appl. Mech.,98, 194–197), an explicit expression for the pulse speed is obtained and the effect of lumen pressure and the axial stretch on wave speed is discussed. Numerical results indicate that the wave speed increases with lumen pressure but decreases with the axial stretch. The results of the present model are compared with our previous work (Demiray, 1988,J. Biomech. 21, 55–58) on the same subject.

Notes: Abstract For a one-locus selection model, Svirezhev introduced an integral variational principle by defining a Lagrangian which remained stationary on the trajectory followed by the population undergoing selection. It is shown here (i) that this principle can be extended to multiple loci in some simple cases and (ii) that the Lagrangian is defined by a straightforward generalization of the one-locus case, but (iii) that in two-locus or more general models there is no straightforward extension of this principle if linkage and epistasis are present. The population trajectories can be constructed as trajectories of steepest ascent in a Riemannian metric space. A general method is formulated to find the metric tensor and the surface in the metric space on which the trajectories, which characterize the variations in the gene structure of the population, lie. The local optimality principle holds good in such a space. In the special case when all possible linkage disequilibria are zero, the phase point of then-locus genetic system moves on the surface of the product space ofn higher dimensional unit spheres in a certain Riemannian metric space of gene frequencies so that the rate of change of mean fitness is maximum along the trajectory. In the two-locus case the corresponding surface is a hyper-torus.

Notes: Abstract Most recent models of the immune network are based upon a phenomenological log bell-shaped interaction function. This function depends on a single parameter, the “field”, which is the sum of all ligand concentrations weighted by their respective affinities. The typical behavior of these models is dominated by percolation, a phenomenon in which a local stimulus spreads globally throughout the network. The usual reason for employing a log bell-shaped interaction function is that B cells are activated by cross-linking of their surface immunoglobulin receptors. Here we formally derive a new phenomenological log bell-shaped function from the chemistry of receptor cross-linking by bivalent ligand. Specifying how this new function depends on the ligand concentrations requires two fields: a binding field and a cross-linking field. When we compare the activation functions for ligand-receptor pairs with different affinities, the one-field and the two-field functions differ markedly. In the case of the one-field activation function, its graph is shifted to increasingly higher concentration as the affinity decreases but keeps its width and height. In the case of the two-field activation function, the graph of a low-affinity interaction is nested within the graphs of all higher-affinity interactions. We show that this difference in the relations among activation functions for different affinities radically changes the network behavior. In models that described B cell proliferation using the one-field activation function, network behavior was dominated by low-affinity interactions. Conversely, in our new model, the high-affinity interactions are the most significant. As a consequence, percolation is no longer the only typical network behavior.

Notes: Abstract The dynamics of a microbial community consisting of a eucaryotic ciliateTetrahymena pyriformis and procaryoticEscherichia coli in a batch culture is explored by employing an individual-based approach. In this portion of the article, Part I, population models are presented. Because both models are individual-based, models of individual organisms are developed prior to construction of the population models. The individual models use an energy budget method in which growth depends on energy gain from feeding and energy sinks such as maintenance and reproduction. These models are not limited by simplifying assumptions about constant yield, constant energy sinks and Monod growth kinetics as are traditional models of microbal organisms. Population models are generated from individual models by creating distinct individual types and assigning to each type the number of real individuals they represent. A population is a compilation of individual types that vary in a phase of cell cycle and physiological parameters such as filtering rate for ciliates and maximum anabolic rate for bacteria. An advantage of the developed models is that they realistically describe the growth of the individual cells feeding on resource which varies in density and composition. Part II, the core of the project, integrates models into a dynamic microbial community and provides model analysis based upon available data.

Notes: Abstract Mathematical models are developed for the chemotherapy of AIDS. The models are systems of differential equations describing the interaction of the HIV infected immune system with AZT chemotherapy. The models produce the three types of qualitative clinical behavior: anuninfected steady state, aninfected steady state (latency) and aprogression to AIDS state. The effect of treatment is to perturb the system from progression to AIDS back to latency. Simulation of treatment schedules is provided for the consideration of treatment regimes. The following issues of chemotherapy are addressed: (i) daily frequency of treatment, (ii) early versus late initiation of treatment and (iii) intermittent treatment with intervals of no treatment. The simulations suggest the following properties of AZT chemotherapy: (i) the daily period of treatment does not affect the outcome of the treatment, (ii) treatment should not begin until after the final decline of T cells begins (not until the T cell population falls below approximately 300 mm−3) and then, it should be administered immediately and (iii) a possible strategy for treatment which may cope with side effects and/or resistance, is to treat intermittently with chemotherapy followed by interruptions in the treatment during which either a different drug or no treatment is administered. These properties are revealed in the simulations, as the model equations incorporate AZT chemotherapy as a weakly effective treatment process. We incorporate into the model the fact that AZT treatment does not eliminate HIV, but only restrains its progress. The mathematical model, although greatly simplified as a description of an extremely complex process, offers a means to pose hypotheses concerning treatment protocols, simulate alternative strategies and guide the qualitative understanding of AIDS chemotherapy.

Notes: Abstract A stochastic cellular automata model for the population dynamics of the army antEciton burchelli on Barro Colorado Island in Panama is set up. It is simulated on the computer and shown to give good agreement with biological data. It is analysed using two approximations akin to the mean field approximation in statistical mechanics, and good agreement with the simulations is obtained. Finally, the role of distance between successive statary phase bivouacs is discussed with regard to the rate of colony growth. There are two aspects of the biological system studied here that make it of general importance. First, the population is structured, since the size of each colony of army ants is crucial. Second, the spatial behaviour of the population, as in many others, is not diffusion-like, although it is random. This has implications for the kind of model that is chosen.

Notes: Abstract An analytical model for thermal damage of retinal tissue due to absorption of laser energy by finite-sized melanin granules is developed. Since melanin is the primary absorber of visible and near-IR light in the skin and in the retina, bulk heating of tissue can be determined by superposition of individual melanin granule effects. Granules are modeled as absorbing spheres surrounded by an infinite medium of water. Analytical solutions to the heat equation result in computations that are quick and accurate. Moreover, the model does not rely on symmetric beam profiles, and so arbitrary images can be studied. The important contribution of this model is to provide a more accurate biological description of submillisecond pulse exposures than previous retinal models, while achieving agreement for longer pulses. This model can also be naturally extended into the sub-microsecond domain by including vaporization as a damage mechanism. It therefore represents the beginning of a model which can be applied across the entire pulse duration domain.

Notes: Abstract Most of the elastic tubes found in the mammalian body will collapse from a distended circular cross section and when collapsed may undergo flow-induced oscillations. A mathematical model describing fluid flow in a collapsible tube is analysed using the software package AUTO-86. AUTO-86 is used for continuation and bifurcation problems in systems of non-linear ordinary differential equations. The model is a third-order lumped-parameter type and is based on the classical “Starling resistor”; it describes the unsteady flow behaviour and, in particular, the experimentally observed self-excited oscillations, in a way which is simple enough to give physical understanding, yet still firmly based on fluid mechanical principles. Some of the bifurcation types found in this model bear close resemblance to the types suggested by experimental observations of self-excited oscillations in collapsible tubes; they thus shed some light on the various topological changes which occur in practice, particularly in view of the fact that some of the points found numerically are diffcult to achieve experimentally, while the existence of others can only be inferred indirectly and uncertainly from experiment.

Notes: Abstract A solution algorithm yielding the pressure and flow-rate distributions for steady flow in an arbitrary, tree-like network is provided. Given the tree topology, the conductance of each segment and the pressure distribution at the boundary nodes, the solution is obtained from a simple recursion based on perfect Gauss elimination. An iterative solution method using this algorithm is suggested to solve for the pressure and flow-rate distributions in an arbitrary diverging-converging (arterial-venous) network consisting of two tree-like networks which are connected to each other at the capillary nodes. A number of special solutions for tree-like networks are obtained for which the general algorithm is either simplified or can be replaced by closed form solutions of the pressure and flow-rate distributions. These special solutions can also be obtained for each tree of diverging-converging networks having particular topologies and conductance distributions. Sample numerical results are provided.

Notes: Abstract Previous game-theoretic models of reciprocity have assumed that populations are large and organisms effectively sessile. This paper analyzes an iterated prisoner's dilemma among non-sessile organisms in a finite population, on the assumption that an individual's chance of remaining in one place is not influenced by a partner's behavior. This mode of interaction is suitable for analyzing potentially cooperative behaviors that are secondary to the advantage of group formation, e.g. allogrooming among social mammals. The analysis yields necessary conditions for stable reciprocity in terms of three parameters, namely, a benefit/cost ratio, the probability of further interaction and the probability of partner retention. The results suggest that, in highly mobile organisms such as fish, birds and mammals, reciprocity may be stable only if the population is small and the relative benefit and future interaction probability are both large.

Notes: Abstract We consider the kinetics of an autocatalytic reaction network in which replication and catalytic actions are separated by a translation step. We find that the behaviour of such a system is closely related to second-order replicator equations, which describe the kinetics of autocatalytic reaction networks in which the replicators act also as catalysts. In fact, the qualitative dynamics seems to be described almost entirely be the second-order reaction rates of the replication step. For two species we recover the qualitative dynamics of the replicator equations. Larger networks show some deviations, however. A hypercyclic system consisting of three interacting species can converge toward a stable limit cycle in contrast to the replicator equation case. A singular perturbation analysis shows that the replication-translation system reduces to a second-order replicator equation if translation is fast. The influence of mutations on replication-translation networks is also very similar to the behavior of selection-mutation equations.

Notes: Abstract In this work, we studied the propagation of non-linear waves in a pre-stressed thin elastic tube filled with an inviscid fluid. In the analysis, analogous to the physiological conditions of the arteries, the tube is assumed to be subject to a uniform pressureP 0 and a constant axial stretch ratio λz. In the course of blood flow it is assumed that a large dynamic displacement is superimposed on this static field. Furthermore, assuming that the displacement gradient in the axial direction is small, the non-linear equation of motion of the tube is obtained. Using the reductive perturbation technique, the propagation of weakly non-linear waves in the long-wave approximation is investigated. It is shown that the governing equations reduce to the Korteweg-deVries equation which admits a solitary wave solution. The result is discussed for some elastic materials existing in the literature.

Notes: Abstract The present paper deals with the bifurcation analysis of a simple food chain model consisting of components like detritus, nutrients, microorganisms, phytoplankton and zooplankton in an aquatic environment. The food chain model is described by a system of differential equations. If the length of the food chain (LFCH) is equal to 3 or 4, then an asymptotically stable equilibrium exists. For LFCH=5 or 6 the non-trivial equilibrium is unstable and the food-chain model has periodic orbits.

Notes: Abstract We develop a macroscopic model for delivering drug to brain tumors. The model accounts for bulk convective and diffusive transport across the blood-brain barrier and through the interstitial space. Through mathematical analysis and simulations, we assess the effects of changing parameters (within physiological bounds) on drug delivery. We find that there is an optimal treatment for convective drug delivery to the center of the tumor. We interpret this phenomenon in terms of traffic flow. The implications of our analyses on existing chemotherapeutic protocols are discussed.

Notes: Abstract The concept of shape space proposed by Perelson and Oster (1979,J. Theor. Biol. 81, 645–670) has been a useful tool for theoretical immunologists, who have invoked it to model idiotypic binding, which plays a significant role in mathematical models of immune networks. The actual construction of such a space from its definition requires specialized experimental information, which is not completely available. In this article, we discuss, with illustrative examples, how graphical representations similar to the idea of shape space can be derived by analyzing real affinity matrices, and the relative merits of such representations to approximations that might be obtained by the approach of Perelson and Oster. We also give directions for future research with a view toward applications.

Notes: Abstract A novel intact circular dsDNA supercoil is proposed as an alternative to the conventional DNA supercoil, so that the two complementary strands of ssDNA circles are separable without any covalent bond breakage. This new structure can be visualized by using two tubings: one black and one clear. Twist the black tubing a number of times and connect its two ends. Do the same for the clear tubing. Then wrap the two tubings together. This forms the separable or novel supercoil. On the other hand, the conventional supercoil can be modeled by twisting the black and clear tubings together and then connect their respective ends, so that the two tubings are not separable unless one of them is cut. Experimentally, in the absence of any enzyme, many intact plasmid dsDNA circles give two bands on agarose gel electrophoresis under a certain given condition, while the same plasmid molecules after cutting once by a restriction enzyme give only one band under the same, condition. In the case of intact pUC19 plasmids, these two bands can then be, recovered and sequenced separately, using two primers in opposite directions. Each band gives mostly one sequence which is complementary to that of the other band. The combination of the above theoretical model and experimental results strongly suggests that there is an alternative structure of DNA which does not have the usual difficulty of unwinding, rewinding and requiring numerous covalent bond breakages and ligations during semiconservative replication.

Notes: Abstract A method of analyzing thymidine labeling in a population of cells is formulated. The formulation establishes a unique relation between a specific set of labeling data and a specific set of cells in the population, viz. that set of cells having a particular chromosome number. The analysis employs a cell-state variable, i.e., a quantity which specifies the progress of a cell through its lifecycle. This variable is defined in terms of the nucleo-protein content and configuration of the chromosomes. The relation mentioned above leads immediately to an expression for the number of cells present at a particular time following labeling which have a given amount of label per cell and a given chromosome number.

Notes: Abstract An equation relating radiation-induced metaphase delay to the dose-rate and duration of irradiation is obtained. The equation is derived from a model specifying the effects of radiation on the normal chromosome coiling process. The basic assumptions of the model are (1) that normal coiling proceeds by contractile protein acting on segements of a viscoelastic chromosomal fiber; (2) that radiation causes cross-linking of adjacent chromosomal fibers which hinders the coiling process.

Notes: Abstract Normal micturition is controlled primarily by a neural system. Certain physical effects become evident when neural control is destroyed, and the automatic or autonomous bladder phenomena occur. It is shown in this paper that a physical system simulating the alternating periods of continence and voiding of the automatic bladder may comprise only passive elastic components, and that periodic voiding does not per se imply neural control.

Notes: Abstract A number of inaccuracies in previous papers are pointed out and amended, and some implications of the correct situation are outlined.

Notes: Abstract We investigate a model of the renal medulla in which active NaCl transport is restricted to the thick ascending limb of Henle's loop. The model contains a vas rectum, a loop of Henle, salt, and water. The model generates interstitial osmolality curves consonant with the known functioning of the kidney in water diuresis. Using data from the white rat and the curves generated by the model, one can predict the permeability of the thin limb of Henle's loop to NaCl and the percentage of total renal blood flow entering the inner medulla. In this model interstitial osmolality at the papilla can be about twice plasma osmolality, so that NaCl transport restricted to the outer medulla can contribute significantly to the work required in producing a hypertonic urine. However, the interstitial osmolality monotonically decreases proceeding from the junction of the outer and inner medulla to the papilla, and the maximum interstitial osmolality in the outer medulla is greater than the maximum interstitial osmolality in the inner medulla. Thus we infer that a source of active transport located in the inner medulla is needed to explain the high osmolalities observed in hydropenia. A sketch of an alternative model, a “lineal multiplication mechanism”, for the renal concentrating process is presented in which active transport in the inner medulla is restricted to active salt transport by the collecting duct. The lineal multiplication mechanism makes no use of counter-current multipliers in the inner medulla.

Notes: Abstract Starting from the basic flux equation, it is possible to obtain an integral form relating the current componentsI i at an arbitrary pointr 2 to the distribution of mobilities and concentrationsc i, potential forces $$\bar \mu $$ , and chemical productivityp i without any restrictive assumptions such as constant mobilities, constant field, steady state, or electrical neutrality. The equation is $$\begin{gathered} I_i (r_2 ) = G_i (r_2 )\left[ {\Delta \bar \mu _i - \int_{r_1 }^{r_2 } {z_i } FA\left( {p_i - dc_i /dt} \right)\left( {\frac{1}{{G_i (r)}}} \right)dr} \right]; \hfill \\ G_i (r) = 1/\int_{r_1 }^r {\frac{{dr}}{{z_i^2 F^2 c_i u_i }}.} \hfill \\ \end{gathered} $$ On the basis of this equation, it is possible to give a more general and systematic development of the basic equations of electrophysiology which clarifies a number of questions concerning the physical interpretation of and the necessary and sufficient conditions for the applicability of some of the standard equations and gives their proper extensions to more general conditions. It is found that the relation between the current components and chemical reactions present arises in a very natural way via the continuity equation and enables one to discuss the incorporation of the metabolic and active transport parameters by assuming a very general physical condition. On the basis of this general integration technique one may then compare the physical interpretation of the differential conductance, the chord conductance and the integral conductance.

Notes: Abstract Previous papers by F. M. Snell (Jour. Theor. Biol.,8, 469–479, 1965) and M. A. Fox and H. D. Landahl (Bull. Math. Biophysics,27, Spec. Issue, 183–190, 1965) have found that the formulation by previous authors for the oxygen flow rates through hemoglobin solution as a function of pressure determined by E. Hemmingsen and P. F. Scholander (Science,132, 1379–1381, 1960) did not give a satisfactory quantitative fit of the curve for constant pressure difference. The suggestion of Fox and Landahl that the Bohr effect involving the shift in acidity accompanying the oxidation of Hb should give rise to voltage and pH differences in oxyhemoglobin transport is examined in more detail. In this paper, the previous expressions for the total oxygen flow rate in terms of the end point concentrations are extended to include the effects of the electrical field. Estimates of the potential difference shows it to be negligible. A derivation of a voltage-pH relation shows that the Nernst relation does not apply and a negligible voltage difference does not preclude a pH shift which is the more probable explanation of the discrepancies observed. Several other predictions suitable for experimental testing are made.

Notes: Abstract A model for the human eye-movement mechanism is derived. The derivation is based on a literature search directed toward identifying and mathematically describing each component through physiological and anatomical considerations. It is felt that although certain parameter values may not be exactly correct (for the data were taken from a wide variety of animals), we can place a great deal of confidence in the configuration.

Notes: Abstract The urethra as seen on X-ray films may show alternate regions of constriction and distension. That these regions do not necessarily correspond to high and low tensions in the circumferential muscle sheath is shown by calculated stable configurations under uniform tension.

Notes: Abstract The discussion as to whether societies are organisms andvice versa has been going on for a long time. The question is meaningless unless a clear definition of the term “organism” is made. Once such a definition is made, the question may be answered by studying whether there exists any relational isomorphism between what the biologist calls an organism and what the sociologist calls society. Such a study should also include animal societies studied by ecologists. Both human and animal societies are sets of individuals together with certain other objects which are the products of their activities. A multicellular organism is a set of cells together with some products of their activities. A cell itself may be regarded as a set of genes together with the products of their activities because every component of the cell is either directly or indirectly the result of the activities of the genes. Thus it is natural to define both biological and social organisms as special kinds of sets. A number of definitions are given in this paper which define what we call here organismic sets. Postulates are introduced which characterize such sets, and a number of conclusions are drawn. It is shown that an organismic set, as defined here, does represent some basic relational aspects of both biological organisms and societies. In particular a clarification and a sharpening of the Postulate of Relational Forces given previously (Bull. Math. Biophysics,28, 283–308, 1966) is presented. It is shown that from the basic definitions and postulates of the theory of organismic sets, it folows that only such elements of those sets will aggregate spontaneously, which are not completely “specialized” in the performance of only one activity. It is further shown that such “non-specialized” elements undergo a process of specialization, and as a result of it their spontaneous aggregation into organismic sets becomes impossible. This throws light on the problem of the origin of life on Earth and the present absence of the appearance of life by spontaneous generation. Some applications to problems of ontogenesis and philogenesis are made. Finally the relation between physics, biology, and sociology is discussed in the light of the theory of organismic sets.

Notes: Abstract Ureteral obstruction has been known to reduce the renal concentration gradients of chloride, sodium, urea, and osmolality. A quantitative analysis of the factors responsible for the decrease in the gradients has been performed. By applying the equation of conservation of matter to data obtained in this laboratory it is concluded that:a. at least, 62.5 per cent of the decrease in concentration gradients is due to drainage by circulatory vessels;b. at most, 25 per cent of the decrease is due to increase in water content of medulla; andc. at most, 12.5 per cent of the decrease is due to diffusion of solutes through the interstitium.

Notes: Abstract The vector equation for the general motion of a body in an inertial system is used to analyze the accelerations in the semicircular canals of the cat when the head undergoes rotation about a vertical axis only, rotation about the naso-occipital axis only, and both rotations simultaneously. The corresponding mean forces and mean pressures in the endolymph are calculated by means of a closed line integral along each canal circumference. The importance of the area of the semicircular canal and of its orientation in space become evident. One can see through this mathematical analysis that the input pattern received by the labyrinthine system depends on a set of well-specified geometrical and mechanical conditions, which must be precisely evaluated in order to interpret the nystagmic outputs.

Notes: Abstract This paper describes some analytical models for the system which regulates the daily eclosion rhythm of the drosophila. A general topological model is described which can simulate practically all the known experimental results about the behavior of the system under various light stimuli. From that a more specific model is proposed which can shortly be described as follows: The system contains a basic oscillator whose output is a substances. This is produced in the presence of an enzymer. During part of the cycler is deactivated ands dissipates until it reaches a lower level whenr is reactivated again. Light has the effect of deactivatingr immediately. The substances causes the production of a second substanceq which triggers a series of reactions leading to eclosion when it exceeds a threshold. The main oscillator (s—r) is temperature-compensated, but the production ofq is accelerated in the presence of light or higher temperature.

Notes: Abstract Some theoretical results obtained in a previous publication (Bull. Math. Biophysics,28: 25–50, 1966) are studied from the numerical point of view. Possible medical interpretations are suggested.

Notes: Abstract The findings of molecular biology concerning biosynthesis of macromolecules are applied to the deduction of the kinetics of mass and volume growth in individual cells between divisions. The time course of increase of all macromolecules and of the total dry mass is found to be linear, in agreement with the available data; the corresponding volume growth curves are either quadratic, or exponential with a linear asymptote, depending on the relative contributions of metabolism and transport to cell water. A self-limiting mass and volume kinetics is derived by including repression among the other molecular mechanisms.

Notes: Abstract An error in a previous analysis of a simple model for myelinated nerve propagation is pointed out and rectified. The model previously chosen does not lend itself to analysis because it has a region ofinfinite negative conductivity which leads to an instability in steady propagation. A model which assumes afinite negative conductivity is examined in detail and a more general results is derived.

Notes: Abstract In a deterministic model for the spread of news in a closed homogeneously mixing population of individuals who never forget (or, of an epidemic without recovery), it is shown that the fractions πi of the population first hearing the news (contracting the disease)i th hand are given by the terms of a truncated Poisson distribution.

Notes: Abstract A mathematical model of the left ventricle is constructed. This model is a spherical thinwalled chamber with a radius and wall density approximately that of the left ventricle and enclosing a fluid whose pressure is that noted when heart sounds occur. The solution for the natural frequency of vibration is of the same value as that noted physiologically. Substitution of the values for the right ventricle yields similar results. The lowest natural frequencies of the idealized cardiac chamber possessing physical properties similar to a real cardiac chamber are of the same order of magnitude as the lowest significant observed cardiac frequencies. Such observation reinforces the likelihood that cardiac sounds are generated by the vibrating cardiac walls.

Notes: Abstract From a quantum mechanical standpoint, electron tunneling may occur in a nerve axon because the nerve membrane (60 to 100 angstroms) is thin enough for the tunnel effect and because the external solution and the axoplasma are redox systems which can provide electrons. Calculations and diagrams of the density-of-states are given to predict theN-shaped current-voltage characteristics which have been observed in the studies of squid giant axons, of the reconstituted liquid membranes and of the marine algae.

Notes: Abstract A theory of growth of a cell which takes up nutrients by diffusion or active transport is discussed. The main conclusion is that the volume should grow at least as fast as the third power of the time. Existing experimental evidence is not a conclusive test of the theory, and further experiments to test it are proposed.

Notes: Abstract This is an analytic study of mucous flow caused by ciliary motion. The computed flow data may be compared with that already found exprimentally. The effects of mucous density, viscosity and layer depth on flow phenomena are investigated. The effects of cilia diameter, length, spacing and oscillation frequency are determined from the equations governing the flow of the mucous blanket. A pertinent finding of the analysis is that the mucous flow in the airway tubes can satisfy physical constraints only through the assumption of a variable viscosity in the covering mucous blanket. The mucous viscosity must increase considerably from the low value at the cillium layer to a much higher value at the air-mucus interface.