ALBERT

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 A stochastic approximation problem for polynomic operators is formulated. The performance of polynomic operators of various degrees is compared. The effect of causality constraints is also examined.

Notes: Abstract A recently developed algebraic approach to the feedback system design problem is reviewed via the derivation of the theory in the single-variate case. This allows the simple algebraic nature of the theory to be brought to the fore while simultaneously minimizing the complexities of the presentation. Rather than simply giving a single solution to the prescribed design problem we endeavor to give a complete parameterization of the set of compensators which meet specifications. Although this might at first seem to complicate our theory it, in fact, opens the way for a sequential approach to the design problem in which one parameterizes the subset of those compensators which meet the second specification...etc. Specific problems investigated include feedback system stabilization, the tracking and disturbance rejection problem, robust design, transfer function design, pole placement, simultaneous stabilization, and stable stabilization.

Notes: Abstract We describe a recently developed framework for exploring the structure of linear time-invariant models of large systems, and for constructing interpretable or physically-based, reduced-order models that reproduce selected modes of the original systems to a desired accuracy. Application of this framework to constructing lumped approximations for interconnections of lumped and distributed systems is briefly explored.

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.

Notes: Abstract In a previous paper, in which a possible mathematical approach to history was outlined, it was shown that urbanization plays an important part in the propagation of new ideas. The rate of such propagation influences the rate of historical developments. The present paper deals in more detail with possible mechanisms of formation of earliest cities. Equations are derived which give the limiting size of such cities and their rate of growth. Of particular importance for the spread of new ideas is the spread of information. The latter largely depends on the fraction of individuals who travel between city and country. Expressions for this quantity are derived. An approach is outlined to the mathematical study of the earliest social classes, which may have been formed as a result of military, religious, or economic stratifications.

Notes: Abstract The velocity of propagation of a disturbance wave in a liquid flowing in a distensible tube is computed. The mathematical model is more general than those used in previous analyses: the tube wall properties are realistic; the convective part of the axial inertia forces is taken into account; radial inertia forces of both the fluid and tube wall are present; viscous stresses are present. Four parameters influencing the velocity of propagation are obtained and discussed. Curves are plotted illustrating the effects of the parameters. Contrary to the results of previous analyses, viscous effects are shown to be appreciable in blood flow. It is also shown that radial inertia effects can be important in laboratory set-ups.

Notes: Abstract A general theory of the drying of frozen tissue is developed and applied to the measurement of the drying rate of frozen guinea pig liver. It is shown that for a given temperature of the subliming ice crystals the mininum drying time of a piece of guinea pig liver is greater than the minimum sublimation time of a piece of ice of the same size and shape by a factor of the order of one thousand. This fact has many implications in the design of freeze-dry apparatus which will be developed in a following paper.

Notes: Abstract The equations governing the time course of the exchange of substances between the blood in the capillaries and the extracellular space are solved for the case of substances which do not penetrate the cells. The equations given relate the time course of the exchange process to the various tissue and circulation parameters such as the specific capillary wall area, the pore area, the inter-capillary distance, the size of the extra-vascular, extra-cellular space, the diffusion coefficient in this space, and the velocity of blood in the capillaries. Some experimental work on capillary exchange is discussed in relation to the theory and estimates are made of the relative importance of the various tissue and circulation parameters in the exchange of substances in different tissues.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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