ALBERT

Notes: Abstract The partial differential equation of the random walk problem with persistence of direction and external bias is derived. By persistence of direction or internal bias we mean that the probability a particle will travel in a given direction need not be the same for all directions, but depends solely upon the particle's previous direction of motion. The external bias arises from an anisotropy of the medium or an external force on the particle. The problem is treated by considering that the net displacement of a particle arises from two factors, namely, that neither the probability of the particle traveling in any direction after turning nor the distance the particle travels in a given direction need be the same for all directions. A modified Fokker-Planck equation is first obtained using the assumptions that the particles have a distribution of travel times and speeds and that the average time of travel between turns need not be zero. The fional equation incopporating the assumption of a persistence of direction and an external bias is then derived. Applications to the study of diffusion and to long-chain polymers are then made.

Notes: Abstract The transmission of some information or behavior pattern is treated as a flow of “particles” which execute random motions over a population of individuals and which may multiply or disappear. Equations are derived for the number density of these “particles” and from this is calculated the number of individuals through which the “particles” have passed. The results are applied to a number of situations such as 1) uniform spatial distribution with multiplication factor decreasing with time because of loss of interest or confusion of the information, 2) multiplication factor constant but the rate of spreal decreasing with multiple hearings, 3) one-dimensional region with a small starting region with or without an absorbing barrier 4) two-dimensional region with absorbing barrier, 5) continous sources of information within a small region in one dimension, 6) uniform spatial distribution in which individuals do not respond to more than one hearing.

Notes: Abstract A situation is considered in which a fluid containing a substance flows through a vessel at a constant rate, the substance being permeable to the vessel wall. In the region outside the vessel there is supposed to be rapid mixing in the direction perpendicular to the axis of the vessel but no mixing longitudinally. The solution for the spatial distribution at any time is given for the case of an arbitrary initial distribution along the vessel length in the absence of an input. The solution is also given for the case of a single impulsive input, the concentration being initially zero everywhere.

Notes: Abstract Organisms orient themselves to a stimulus by two general methods. One method is by directed orientation (taxis); the other is by undirected locomotory reaction (kinesis). An equation, and the methods for finding the necessary parameters of this equation, is derived for the distribution of organisms within a container, with the following limitations: (1) the organisms have no accommodation, (2) they are always active, and (3) the stimulus changes slowly with position. Necessary modifications of the equation are then derived, so that the last two limitations may be eliminated. The equation cannot be solved excatly because of its complexity; hence an approximation method must be used. This method is discussed, an approximate solution is found, and a time constant for equilibrium to be established is derived. Applications tovarious experiments in the literature are then made with fairly satisfactory results. A new interpretation of the theory of klino-kinesis with accommodation is found upon application of the equations developed to experimental work. Further limitations and uses of these equations are then discussed.

Notes: Abstract Certain parameters are defined which roughly characterize the internal structure of networks. A given network structure uniquely determines the values of the parameters, but the reverse is not true. The parameters therefore define certain classes of networks. One of the parameters, thedispersion D(S) gives an indication of the “compactness” of the internal structure. Addition theorems and inequalities are derived relating the dispersions of sub-systems to the dispersion of the complete structure.

Notes: Abstract A mathematical theory is developed which permits the determination of certain parameters of an inhomogenous tissue, such as a nerve trunk without its epineurium. The parameters are the permeability coefficients for entrance into an exit of a substance from the nerve fibers, and the diffusion coefficient of the interstitial material. The experimental data required are the dimensions of the cross-section, the average diameter of the fibers, and the ratio of the cross-sectional are of the fibers to the total cross-section, as well as the time course of the decrease of the fraction of the substance left in the nerve trunk, when the trunk is immersed in a bathing solution containing none of it.

Notes: Abstract A model based on enzyme localization is developed which gives rise to an apparent active transport of a metabolite into or out of cells. The model is applied to three simple situations, using Fick's equation and the Rashevsky approximation. It is shown that the apparent efficiency can be made as large as desired if, for constant reaction, the outer cell region is made sufficiently small, or, for autocatalytic reaction, if the metabolite concentration in the outer region is sufficiently small. The physical limitations imposed by this mechanism are developed for all three situations.

Notes: Abstract A previously derived iteration formula for a random net was applied to some data on the spread of information through a population. It was found that if the axon density (the only free parameter in the formula) is determined by the first pair of experimental values, the predicted spread is much more rapid than the observed one. If the successive values of the “apparent axon density” are calculated from the successive experimental values, it is noticed that this quantity at first suffers a sharp drop from an initial high value to its lowest value and then gradually “recovers”. An attempt is made to account for this behavior of the apparent axon density in terms of the “assumption of transitivity”, based on a certain socio-structural bias, namely, that the likely contacts of two individuals who themselves have been in contact are expected to be strongly overlapping. The assumption of transitivity leads to a drop in the apparent axon density from an arbitrary initial value to the vicinity of unity (if the actual axon density is not too small). However, the “recovery” is not accounted for, and thus the predicted spread turns out to beslower than the observed.

Notes: Abstract The assumption of transitivity treated in part I is modified in various ways to describe an information-diffusion process, in which a certain amount of randomness of contact does occur. In one model a parameter is introduced which is indicative of a tendency to go beyond one's immediate vicinity to spread the information as the vicinity becomes saturated with knowers. In another model the randomness appears in the assumption that new knowers are uniformly distributed among the knowers. Two of the equations thus derived, each with two free parameters are in good agreement with experimental results.

Notes: Abstract It is observed that a dynamical continuity equation for biomass distribution yields the asymptotic steady-state exponential dependencen=A exp( $$ - m/\bar m$$ ) exhibited by certain fishery data, wherem is the biomass of an individual,n is the number of individuals per unit biomass interval, andA, $$\bar m$$ are positive constants. This dynamical approach to biomass distribution is an alternative to the global maximization principle proposed recently by Lurié and Wagensberg.

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 Linear birth and death processes are used to derive simple expressions for sequential extinction times and gene fixation probabilities in asexual populations.

Notes: Abstract We are here concerned with the functionf which assigns to each pointP of an object the numberf(P) which is the shortest distance fromP to the border. This function appears in various guises in diverse biological studies. The functionf(P) is itself a measure of shape—or more precisely, an infinite set of measures, one for each point (and hence, in view of its geometric definition, usually in a form inconvenient for use). Thus in this paper we sought a reasonable representative of this infinite set of measures, namely themean of the numbersf(P) asP ranges over all points of the entity. Computability studies are developed for various classes of shapes. For example, (1) the mean for a lamina bounded by a polygon circumscribable about a circle of radiusr isr/3; (2) the mean for a domain bounded by a polyhedron circumscribable about a sphere of radiusr isr/4. The transition from pointwise to piecewisef(P), especially in the non-convex case, requires working with inequalities.

Notes: Abstract A mathematical model for traveling bands of motile and chemotactic bacteria in the presence of cell growth and death is examined. It is found that asymptotic traveling wave solutions exist in the absence of chemotaxis, due to the balance of growth, death and random motility. Thus random motility confers the ecological advantage of population propagation through migration into nutrient-rich regions. The presence of chemotaxis amplifies this advantage by moving more cells into higher nutrient concentration regions, resulting in larger and faster bands. Therefore there seem to be two types of traveling bands that can be attained by chemotactic bacteria in the presence of growth and death: (1) these growth/death/motility bands; and (2) pure chemotactic ‘Keller-Segel'-type bands. Comparison to experimental observations by Chapman in 1973 indicate that the latter seem to be formed. The relationship between these two types of solution is at present uncertain. The growth/death/motility bands may have relevance on longer time or distance scales characteristic of microbial ecological systems.

Notes: Abstract Based on the principle of minimum power, a mathematical model of the pathological functional state of the oxygen transport system is presented. The model is used to determine the optimal functional parameters of the oxygen transport system in hyperthyroidism, anemia and hypertension. Theoretical results are compared with clinical data.

Notes: Abstract A modified SIRS model is developed as a framework for the study of epizootiological dynamics in an insect-pathogen system. Linearized stability analysis reveals that the system with one immune and one susceptible host class can exhibit stable, periodic or unstable behavior depending on model parameters. In general, high pathogenicity, short pathogen propagule lifespan and high host reproductive rate are stabilizing influences. Pathogen transmissibility and propagule production/host do not influence local stability. The effect of seasonal host reproduction is studied because most insect hosts are seasonal in temperate climates. The basic stability dependence on model parameters holds except as modified by the length of the reproduction interval. The results of this study are compared with the recent work of Anderson and May.

Notes: Abstract The random predator-prey type interactions of the population species in a random varying environment have been investigated. The Fokker-Planck equation for the transition probability, corresponding to the stochastic differential equations established from Lotka-Volterra equations by the introduction of randomness and variability, has been integrated in the form of a path integral. The transition probabilities for extinction or survival of one or several species have been approximately evaluated and investigated.

Notes: Abstract If two microbial populations compete for a single resource in a homogeneous environment with time invariant inputs they cannot coexist indefinitely if the resource competed for is not renewed by biological activity within the system. Mathematical studies have shown that in a predator-prey system, where the resource (prey) is self-renewing, the two competitors (predators) can coexist in a limit cycle. This suggests that if the resource competed for is renewed by biological activity within the system coexistence can occur in any microbial system provided that it exhibits the same features as, but without being, a predator-prey one. A food chain involving commensalism, competition and amensalism is presented here. Two subcases are considered. It is only when maintenance effects are taken into account that coexistence, in limit cycles, can occur for this system. Limit cycle solutions for the system are demonstrated with the help of computer simulations. Some necessary conditions for coexistence are presented, as are some speculations regarding the possible physical explanations of the results.

Notes: Abstract The branching structure of the mammalian arterial tree has been known to be close to that of an optimal conduit system of the minimum work model characterized as the branch system of constant wall shear rate. The physiological mechanism producing such construction was considered to be based on the local response of arterial caliber induced by the wall shear stress (shear rate × blood viscosity) and thereby maintaining this stress constant, which was previously observed at the canine common carotid artery shunted to the external jugular vein. The stress levels at various parts of the arterial system estimated from available data fell within ±50% of the mean (15 dyn/cm2), which was consistent with the value predicted from the model. Theoretical analyses on the cost function of the model indicated that the suspected variation of shear rate levels in the arterial tree due to the anomalous changes in blood viscosity which might bring about 3- to 4-fold differences between the minimum and maximum shear rates would cause less than 10% increase in the total energy cost. It was concluded that a local adaptive response to wall shear stress is the mechanism which effectively optimizes the design of the arterial tree.

Notes: Abstract The structure of the genetic code is related to a Gray code, which is a plausible theoretical model for an amino acid code. The proposed model implies that the most important factor in shaping the code was the effects of mistakes in translation, not effects of mutations. Another possible implication is that the preservation of stiffness and flexibility at appropriate places in a protein chain is as important in protein structure as the appropriate placement of hydrophilic (external) and hydrophobic (internal) residues. Other results are a simple conceptualization of the relationships among the 20 amino acids and their relations to their codons. The detailed relationships are summarized in the following ‘similarity alphabet’: ala, thr, gly, pro, ser; asp, asn, glu, gln, lys; his, arg, trp, tyr, phe; leu, met, ile, val, cys; (ATGPS DNEQK HRWYF LMIVC in the one-letter code). This alphabet falls into four groups of amino acids: small, external, large, internal. The approximate relation of the groups to their codons is expressed as: the first base of a codon controls size—a purine means a small amino acid, a pyrimidine means large; the middle base controls cloisterednes—purine means external, pyrimidine means internal. These relationships express the minimum change principle upon which the code appears to be founded.

Notes: Abstract A theoretical study of the Brusselator model with non-uniform distribution of component A and a concentration-dependent diffusion coefficient has been performed. Numerical simulation reveals that a variable diffusion coefficient alters the bifurcation pattern and the stability properties of the steady-state as well as periodic solutions. A simple approximate method, based on one-point collocation, has been proposed to analyze the bifurcation phenomena for the case of fixed boundary conditions and low system size.

Notes: Abstract In this paper the effects of the occurrence of cut trees in the topological analysis of branching patterns have been studied. It is assumed that branches are removed at random from the trees. We prove that, for both the segmental and terminal growth models, the probability distributions of the cut trees are identical to those of complete trees.

Notes: Abstract The goal of this work is an examination of capillary exchange models as mathematical operators. The concentration function relations for the Krogh cylinder of a single capillary, basic to many organ models, are studied via the theory of operators on the Lebesgue normed spacesL p[0,∞], (1〈-p〈-∞). A discussion is included of theL p -normsvis-à-vis the coefficient of variation currently used in finding capillary parameters and evaluating parameter searches. The capillary model determines two operators on the space of locally integrable functions: O K (relating extravascular concentration to intravascular) and K a, k (relating intravascular concentration to input), wherek is the ratio of permeabilitysurface area (PS) to extravascular volume, and α is the ratio of PS to flow. These operators are shown to induce contractive (‖O K ‖ p 〈-1, ‖K a, k ‖ p 〈-1), isotone, linear operators onL p . The uniform convergence relation $$K_{a,k} = \mathop {\lim _{(p)} }\limits_{N \to \infty } \left( {\sum\limits_{n = 0}^N {P_n (a)O_k^n } } \right)$$ (as operators onL p) is derived, whereP n (a) is the Poisson probabilitye −a a n /n!. For the important special cases ofp=∞, 1, 2 the norms are found (‖Ok‖=‖Ka,k‖p=1). Consideration is also given to the norms and operators when the functions involved are limited to a finite interval of time.

Notes: Abstract One particular kind of structure offers possible explanations, for long-term memory, efficient consolidation of stored information from the environment, clustering of data strings and multimodal functioning. It is a possible model for pieces of neural structure and its use offers a uniform method for both studying and constructing an extensive class of mechanisms.

Notes: Abstract Levenshtein dissimilarity measures are used to compare sequences in application areas including coding theory, computer science and macromolecular biology. In general, they measure sequence dissimilarity by the length of a shortest weighted sequence of insertions, deletions and substitutions required, to transform one sequence into another. Those Levenshtein dissimilarity measures based on insertions and deletions are analyzed by a model involving valuations on a partially ordered set. The model reveals structural relationships among poset, valuation and dissimilarity measure. As a consequence, certain Levenshtein dissimilarity measures are shown to be metrics characterized by betweenness properties and computable in terms of well-known measures of sequence similarity.

Notes: Abstract It is pointed out that the methane flux measured experimentally for certain ponds and swamps is quantitatively consistent with a commensal dependence of Methanobacteria on O2-chemotactic motile aerobic bacteria. The Methano species is thereby shielded from oxygen and provided with carbon dioxide for the anaerobic production of methane.

Notes: Abstract A sufficient condition is given for stochastic boundedness persistence of a top predator in generalized Lotka-Volterra-type stochastic food web models in arbitrary bounded regions of state space. The main result indicates that persistence in the corresponding deterministic system is preserved in the stochastic system if the intensities of the random fluctuations are not too large.

Notes: Abstract One of Bobisud's models for the evolution of cannibalism is reanalyzed by applying the method of finding evolutionarily stable strategies (or ESS's). It is demonstrated that ‘no cannibalism’ never will be an ESS if the initial rate of cannibalism is too large. It is further demonstrated that individual selection may even result in the evolution of cannibalism during food abundance. Some empirical case studies are briefly discussed in relation to this model.

Notes: Abstract A computer algorithm is presented which equiprobably generates any member of the set of all directed trees withk labeled terminal nodes and unlabeled interior nodes. The algorithm requires roughlyk 2 /2 storage locations. The one-time initialization requiresO(k 2 ) time, while generating each tree requiresO(k) time.

Notes: Abstract The comparison of several sequences is central to many problems of molecular biology. Finding consensus patterns that define genetic control regions or that determine structural or functional themes are examples of these problems. Previously proposed methods, such as dynamic programming, are not adequate for solving problems of realistic size. This paper gives a new and practical solution for finding unknown patterns that occur imperfectly above a preset frequency. Algorithms for finding the patterns are given as well as estimates of statistical significance.

Notes: Abstract A new development is introduced here in the use of dynamic programming in finding pattern similarities in genetic sequences, as was first done by Needleman and Wunsch (1969). A condition of pattern similarity is defined and an algorithm is given which scans any set of similarities and screens out those which fail to meet the condition. When the set to be scanned contains every pair of segments, one from each of two given sequences of lengthsm andn (i.e. every possible location for a pattern similarity), then it completes the scan in a number of computational steps proportional tom·n, leaving those pairs of segments which satisfy the similarity condition. The algorithm is based on the concept of match density, as suggested by Goad and Kanehisa (1982).

Notes: Abstract This paper concerns sequences of letters in which certain “distinguished” words are of interest. Such sequences arise as data in numerous fields including genetics and neuroscience. A probability distribution is given for the number of occurrences of a chosen word in a randomized sequence of letters. Such words are considered “favored” if they occur more than expected at random. Favored words have been discovered in nerve impulse trains and may reflect a neural coding scheme.

Notes: Abstract As the volume of protein sequence data grows, rapid methods for searching the protein sequence database become of primary importance. Rigorous comparison of sequences is obtained with the well-known dynamic programming algorithms. However, these algorithms are not rapid enough to use for routinely searching the entire database. In this paper we discuss some methods that can be used for rapid searches.

Notes: Abstract We give algorithms for computing the extent of similarity between two or three sequences of letters. The similarity measures we consider include a penalty for inserting gaps within the sequence in order to enhance similarity. The magnitude of the penalty for gaps is assumed to be independent of their size in order to accommodate certain biological applications. Our algorithm for three sequence comparisons, which is based on solving a system of recursive equations, improves upon the efficiency of existing methods. Although the system of recursive equations utilized by the algorithm is quite complicated as it stands, it has none the less been simplified by appeal to combinatorial considerations.

Notes: Abstract Well-known dynamic programming algorithms exist for comparing two finite sequences inO(N 2) time and storage, whereN is the common sequence length. Extensions to the comparison ofM finite sequences requireO((2N) M) time and storage, making such algorithms difficult even forM=3. A simple generalization of the sequences makes it possible to obtain some results about the geometry of sequence alignments. These ideas suggest heuristic approaches to problems of comparing several sequences. IfM sequences are known to be related by a binary tree, they can be aligned inO(MN 2) time andO(N 2+NM) storage.

Notes: Abstract Protein sequences of the Dayhoff databank of 1984 have been analyzed to evaluate the occurrences of the 400 dipeptides and 8000 tripeptides. Expected values and standard deviations for the di- and tripeptides were determined by Monte Carlo and binomial approximation. A condensed format containing this information, labeled a uniqueness diagram, is presented and made available in the form of a microfiche.

Notes: Abstract In two-state sliding filament models of muscle contraction a partial differential equation must be solved to find the cross-bridge distribution functionn(x, t). In this paper the analytical form of this function is obtained by integration along the characteristic line and special cases are presented in which the explicit expression forn(x, t) can be completely determined. These analytical solutions provide a direct mathematical connection between the microscopic contraction parameters contained in the kinetic theories and macroscopic muscle dynamics and are thus used to investigate what parameters influence the transient contractile tension in typical experimental conditions. The results of this investigation are consistent with relevant aspects of muscle physiology.

Notes: Abstract Experiments on random binary, ternary, etc. (P=2, 3,…, 10) switching nets are reported. Behavioral cycle lengths are examined as functions of output variety,P, input connectance,K, and net size,N. Overall, output variety appears an influential, well-behaved net property. Strong, but well-behaved interactions appear among net variables. In high connectance nets, median cycle length grows approx. asP N/2. Other factors constant, one-connected nets show the shortest cycles, and connectance effects appear to converge asymptotically aroundN. Data for cycle length as a function of net size suggest a concavity not compatible with the Kauffman “square root law” (Kauffman, 1969). Evidence of a positive relationship between cycle length and run-in length is found in two-input nets; weaker evidence is obtained that in higher connectance nets this relationship becomes negative in sign. The “modular complexity” ofP〉2 nets is examined briefly.

Notes: Abstract The influence of topographical situation on the spread of infection is studied. The investigation is based on a multigroup model. The population under consideration is thought to be divided into subpopulations living in regions that are separated from each other by natural barriers (mountains). Infection is carried from one region to another by migrating infectives. Migration is possible only along the river system so that the structure of the epidemiological network is that of a symmetric tree. The results allow comparison of the velocity of propagation of the epidemic for different geographical situations and allow quantification of the “channel-effect”, according to which mountainous regions are channels rather than barriers to the spread of an epidemic.

Notes: Abstract The impressed field, “Window Field” (WF), due to a half-wave action potential on a muscle fiber, has been calculated on the basis of potential theory. It has been shown that in spite of the small intensity of the field, its integrated action can transfer the energy needed to induce, contraction from the membrane to the interior of the fiber. The energy of polarization has been found to be sufficient to exceed the energy of, thermal agitation on that length of fiber, which can be identified as the length of a sarcomere. The changes of ion concentration, caused by the WF, if calculated on the assumption of the semipermeability of theZ membranes, was found to be equal to the changes necessary to induce contraction of actomyosinin vitro.

Notes: Abstract Some general properties of the solution of the diffusion equation are deduced for the steady-state, spherically symmetric system. On the basis of these developments some results of N. Rashevsky (Bull. Math. Biophysics,11, 15, 1949) are discussed and the results of a previous investigation (Hearon,Bull. Math. Biophysics,12, 135, 1950b) are extended to more general conditions. In particular these extensions apply to the flow of a soluteagainst its concentration gradient, the nonzero gradient of an inert metabolite, and theaccumulation or exclusion of an inert metabolite in a metabolic system.

Notes: Abstract The approximation method of N. Rashevsky is discussed and reviewed. It is shown that in addition to theexplicit assumptions and approximations there is involved the assumption that the rate of metabolism is the same at every point in the cell and that theaverage rate of metabolism is different from zero. An expression is given for the error in the approximate method when the rate of metabolism is any function of the concentration. It is also shown that a solution in theform of that obtained by the approximate method is not possible if the generalized laws of diffusion are assumed to apply.

Notes: Abstract Rate equations for the enzymatic oxidation of succinic acid are derived on the assumption that when a single molecule of substrate combines with an enzyme molecule, it can do so with either one or two sites on the enzyme, and that oxidation occurs only in the second case. In addition it is assumed that the product of the reaction, fumaric acid, combines reversibly with the enzyme. With certain enzyme preparations the data fitted such an equation satisfactorily. In others the rate was that of a first-order reaction, but addition of cytochrome changed it to the former type. It was concluded that the transfer of hydrogen to oxygen was a first-order reaction and dominated the whole rate when enzyme preparations were used which had been washed relatively free of cytochrome. When the limiting factor was succino-dehydrogenase the rates followed the new equation. Criteria for recognizing noncompetitive inhibition are given, and inhibition by di-tertiary butyl peroxide was shown to be of this type.

Notes: Abstract This paper deals with the following question: Which distributions of radiosensitivity in a population can lead to an exponential survival curve? The problem is solved exactly, with statistical fluctuations in dose fully accounted for. It is shown that only an exponential distribution of sensitivities can give rise to an exponential survival curve.

Notes: Abstract An approximation method is introduced which enables a number of diffusion-type problems to be solved in an approximate but simple manner. Many cases require only the solution of a simple first-order differential equation. The method is applied to a number of cases in which the exact solutions are available. A comparison shows that the method is quite satisfactory in these cases. The method is applied to diffusion problems with rate of consumption proportional to concentration or to the square of the concentration. In the latter case, the result obtained is essentially the same as that found by H. G. Landau (1950) after elaborate calculations.

Notes: Abstract It is shown that a slight modification of a model of excitatory phenomena in irritable tissues, which has been treated before, exhibits spontaneous oscillations. The frequency of these oscillations and the time-course of the potential across the model membrane have been determined, together with the dependence of some of their characteristics on some important parameters, particularly (Ca++).

Notes: Abstract It is demonstrated that, if the variations of viscosity throughout a cell are considered, swelling stresses may produce elongation and division. To do this it is necessary to generalize Betti's theorem to cover systems containing viscosity gradients and such a generalization is presented. On the basis of two special assumptions it is shown that most of the results of the diffusion drag theory of cell division may be duplicated by the present theory.

Notes: Abstract The theory of imitative behavior, developed previously, is applied to the case of two social groups which are separated spatially. If the information of each group as to the behavior of the other is complete, the case reduces to that of a single group. When any information is lacking at all, the two groups are independent. If we have two mutually exclusive behaviorsA andB, all four combinationsAA, AB, BA, andBB are possible. If the mutual information gradually increases from zero, then for a certain value of it, the group which is more informed about the behavior of the other will change to that behavior if it did not already exhibit it. If for constant information the size of the group increases, then above a certain threshold value, the larger group imposes its behavior on the smaller.

Notes: Abstract Local regulation of blood flow as determined by capillary diameter and the number of open capillaries in a region is considered. The local changes in capillary diameter and in the number of open capillaries are assumed to be due to concentration changes of a diffusible, nonspecified metabolite. This metabolite is produced in the tissue and carried away by the blood stream. Using these assumptions and applying pertinent data on capillaries, deductions are made concerning: (a) the law of blood flow as a function of temperature and capillary radius for the hyperemia of high temperature, (b) high flow as it depends on metabolism during strenuous exercise of muscle, and (c) a first approximation to the time duration of occlusion hyperemia.

Notes: Abstract The necessary and sufficient condition is given forn integers to be the score structure of a society with a dominance relation. A proof is also given for a theorem showing that there are members who dominate every other member either directly or indirectly through a single intermediate member.

Notes: Abstract On the basis of a previous general formulation (Bull. Math. Biophysics,15, 21–29, 1953a) a discussion is given of the error in the approximation method of N. Rashevsky. This error, inherent in the method when the metabolic rate is different at each point in the cell, is discussed in detail and numerical values are presented for two particular cases: the rate proportional to the concentration and the rate a prescribed function of the spatial coordinates. It is shown that the formulation for the first case also applies to several other cases, that the error is negligible provided the rate is sufficiently small, and that the error is fairly sensitive to the cell size. If the rate depends upon the coordinatesalone a small rate is not sufficient to insure a negligible error. The relations between the exact method, the standard approximate method, an earlier approximate method (Physics,7 260, 1936), and a more recent refinement (Bull. Math. Biophysics,10, 201, 1948) of the standard method are discussed.

Notes: Abstract It is shown on the basis of (1) conservation of mass, (2) positive concentrations, and (3) the principle of detail balancing that periodic reactions cannot occur in a closed system described bylinear differential equations. The matrix,A, of the rate equations must be such that |A|=0,a ij〉0 fori≠j,a ii〈0, andVAV −1=B, whereV is diagonal andB is symmetric. These properties ofA imply that the latent roots are real and non-positive and that neither catalysis nor inhibition can be described bylinear equations. It is further shown that periodic reactions cannot occur in anopen system for which the matrix associated with the chemical reactions has the above properties and in which thesimple law of diffusion is obeyed. The relation of these results to Onsager's reciprocal relations and to previous work on periodic and cyclic chemical reactions is discussed. The utility of certain of these results for the treatment of isotope kinetics is indicated.

Notes: Abstract It is demonstrated that an explanation of the small radius effect or the so-called sigma phenomenon may be obtained by noting that one of the effects of the presence of suspended particles in a flowing fluid is to increase the velocity of flow near the wall over that existing in the absence of particles. This effect may be considered equivalent to relaxing the boundary conditions at the wall. An expression for the viscosity is compared with data and fit is found to be good.

Notes: Abstract The solution for the spatial distribution of ions in a Donnan equilibrium has been given by J. H. Bartlett and R. A. Kromhout (1952). The present note gives an explicit solution for the case in which the length of the region containing the membrane is large; in biological situations this requires only that the length considered should be greater than a few hundred Ångstrom units. The Donnan equilibrium may be considered to be a special case of a situation in which forces other than electrical act upon the ions; in particular, it represents the case in which only one ion is acted upon and the energy difference on the two sides of the membrane is infinite. An expression is given for the difference in energy of theith in terms of the electrical potential and of the ion concentrations. As an illustration, the results are applied to nerve membrane potentials.

Notes: Abstract A mechanism is described which accounts for the active transport of Na+ ions through a membrane. It is assumed that at one side of the membrane the ion combines with a carrier ion, the resulting carrier compound then diffuses through the membrane and decomposes at the other side of the membrane. The free diffusion of the ions is also taken into account. The time rate of accumulation of the ion in question at the latter side of the membrane is calculated in terms of the concentrations of the ion at both sides of the membrane.

Notes: Abstract The recent extension of the approximation method is applied to enable us to arrive at the time course of the concentrations at both sides of a membrane. From the differential equations which govern these, the steady-state solution is obtained in terms of the parameters, which are determined by the thickness of the diffusion layers, the chemical composition and reactions, and the diffusion constant of the membrane.

Notes: Abstract An equation is derived from the spread of a “state” by contact through a thoroughly mixed population, in which the probability of transmission depends both on the over-all duration of the process and on the time an individual has been in the “state.” Cases in which this probability is a function of only one or the other of the two “times” are worked out. It is shown that in the case of dependence on “private time” alone the asymptotic value of the fraction of the population effected is the same as that derived by the random net approach.

Notes: Abstract The reflection of pressure waves in a fluid enclosed within a tube with an elastic wall is studied for the case of a localized change in diameter of the tube. The concept of impedance is introduced. The relation of the reflection characteristics of the parts of the tube at either side of the change is derived on the basis of the continuity of pressure and mass flow at the site of the change. This relations is used to derive the expression for the ratio of the pressure oscillations measured in front of, and behind, the constriction in terms of the constants of the system. As a result, a method is indicated to locate the coarctation from measurements of the pressures in front of, and behind it.

Notes: Abstract The question is raised concerning the possible causes of abnormally small standard deviations found in counting samples in which particles are distributed at random (e.g., blood cells, fat globules in milk, etc.). The effect of discarding abnormal samples is discounted inasmuch as small standard deviations occur even when all samples are counted. An approximation method is used to calculate the effect of finite particle size, of known repulsive forces between particles and of convection currents. This calculation shows that neither finite size nor the known repulsive forces are sufficient to account for the observed abnormality of standard deviation, but that convection currents can possibly account for it. The possible presence of long-range repulsive forces cannot, however, be excluded.

Notes: Abstract The effect of finite particle size on the standard deviation in sample counts is computed for the one-dimensional case. To a first order of approximation the correction is found to be identifical with that found by H. de Vries (1953) using a general approximation method.

Notes: Abstract A critical examination of the “classical” theories of photoreception in view of more recent experimental findings yields the result that these theories do not possess the property to describe all the more significant phenomena of photoreception correctly, and to some extent suffer the lack of more general applicability. The basis for a new and presumably more general theory of photoreception based on dynamical aspects is laid out. Emphasis is put on the time course of afferent and efferent excitation in the photoreception model, consisting of a receptor element, an afferent and an efferent neuron of the one-factor Rashevsky-type, and an effector organ.

Notes: Abstract A mathematical model for the development of human society, beginning with the earliest stages of urban cultures, is outlined. In the early stages of history, behavior was characterized largely by adherence to a number of beliefs and prejudices of diffeirent kinds, which were accepted on faith and not subject to critical rational analysis. Due to psychobiological variability a very small number of individuals spontaneously appear at all times who challenge the accepted beliefs and prejudices and do not follow the accepted patterns of social behavior. The effect of these individuals upon the rest of the society, especially upon the younger generation, depends on the facilities with which information spreads in society. In earliest societies, when modern methods of mass communication were unknown, the channels of communication were practically identical with the channels of economic transport. The latter in its turn depended on the nature of the roads, and especially on the presence of waterway, which facilitated transportation. The sizes of the earliest cities and the distances between them were largely determined by relative ease of transporation. Expressions are derived for the average size of the earliest cities and for the average distance between them. The calculated average populations of the earliest cities are of the order of 103; the distance of the order of 102 km. Both are in agreement with some archaeological findings. An expression for the time spaon required for the development from the earliest stages of urban cultures to the present time is derived and shown to depend on the specific shoreline of the country, that is, the length of the shorline divided by the area of the country. It is pointed out that western Europe's specific shoreline, including land bordering both seas and rivers, is ten times as large as the shoreline area of other parts of the world. It is shown that this greater specific shoreline may account quantitatively for the faster social and technological development of western Europe in the last few centuries. The calculated total span of time of development from earliest urban cultures to our days is found to be of the order of magnitude of ten thousand years. It is shown that the model accounts for the existence at the present time of primitive cultures. A number of suggestions is made in regard to other possible applications of mathematics to history.

Notes: Abstract The assumptions latent in the derivation of the integral equation of Branson are rendered explicit and discussed. It is shown that the equation is valid only for systems in which the substance disappears according to a linear rate law.

Notes: Abstract It is shown that the validity of Branson's integral description of metabolizing systems is subject to severe limitations. The validity is insured only in cases where the reaction is of first order, or quasi of first order. In all other cases Branson's equation has to be modified to insure general applicability. The consequences of a different definition of the metabolizing functionF have also been investigated. With the new definitionF describes the pure effect of metabolization. It is found that in this case the integral equation is only capable of describing first-order reactions. With a slight modification of the integral equation it is possible to describe metabolites “with age”, which do not have reactions of definite order, but which satisfy the superposition principle.

Notes: Abstract A problem in probability is stated with included the problem of the distribution of bacterial mutants as a special case. This problem is solved exactly but since the resulting expressions are too complicated for practical use, various approximate expressions for the distribution are considered, especially for the bacterial mutation case.

Notes: Abstract Simple reaction and discrimination reaction, under the influence of imitation, are considered for the situation in which the stimulus or the stimuli vary slowly with time. The result is analogous to hysteresis under certain conditions. The calculations are facilitated by the solution of $$x = \int_{ - \infty }^{a + \beta x} {g\left( \xi \right)d} \xi ,$$ g(ξ) being the normal error function. Values ofx(α, β) are given in a table.

Notes: Abstract On the basis of simple physical considerations the blood flow in a branching circulatory system is studied. The case of two groups of parallel vessels is treated. The vessels of the same group are supposed to be identical. The resistance of each group is determined by the resistance of each vessel in the group and by the number of vessels in the group. From the dependence of the resistance of each vessel on its radius an expression is obtained for the blood flow through each group of vessels in terms of the numbers and sizes of the vessels in each group. The number of open vessels in an organ and the radius of each of those vessels are assumed to depend on the metabolic rate of that organ. The relations so obtained, together with the expression above, are applied to derive the blood flow through an organ as a function of the metabolic rate of that organ. It is indicated that the relations obtained might describe the shifting of blood from one organ to another if the activity of one of them changes. A way is pointed out to treat neural regulation of this phenomenon.

Notes: Abstract A stochastic model of population growth is treated using the Bellman-Harris theory of agedependent stochastic branching processes. The probability distribution for the population size at any time and the expectation are obtained when it is assumed that there is probability (1−σ), 0≤σ〈1, of the organism dividing into two at the end of its lifetime, and probability σ that division will not take place.