Two remarks on the mathematical biology of automobile driving

In connection with a series of previous papers by this author (Bulletin of Mathematical Biophysics,21, 299–308, 375–385;22, 257–262, 263–267;23, 19–29;24, 319–325) results obtained by A. Crawford (Economics 5, 417–428) on the effects of irrelevant lights on reaction times toward a given light stimulus are discussed. The conclusions from a previous paper of this author (Bulletin of Mathematical Biophysics,23, 19–29) are elaborated.     

Probability distributions associated with distinct hits on targets

Some probability distributions connected with distinct hits on targets, using two different firing schemes, are developed. It is assumed that any shot has a probabilityp, not necessarily unity, of hitting the target at which it was aimed. The development uses a well-known expression for the probability that exactlyt ofN possible events occur simultaneously. Some of the formulae developed here include as special cases the probabilities derived separately and by more complicated arguments in papers by N. Rashevsky. (Bull. Math. Biophysics,17, 45–50, 1955) and A. Rapoport (Bull. Math. Biophysics,13, 133–38, 1951).     

Analysis of tracer experiments: IV. The kinetics of generalN compartment systems

The methods of C. W. Sheppard and A. S. Householder (Jour. App. Physcis,22, 510–20, 1951), H. D. Landahl (Bull. Math. Biophysics,16, 151–54, 1954) and H. E. Hart (Bull. Math. Biophysics,17, 87–94, 1955;ibid.,19, 61–72, 1957;ibid.,20, 281–87, 1958) are employed in studying the kinetics of generalN compartment systems. It is shown that the nature of the transfer processes occurring in fluid flow systems and the chemical processes occurring in quadratic systems and in catalyzed quadratic systems can in principle be completely determined for all polynomial dependencies. Systems involving three-body and higher-order interactions can be completely solved, however, only if supplementary information is available. Research supported by the Atomic Energy Commission, Contract AT (30-1)-1551.     

A note on epilepsies and similar disorders

This paper compares two previously published neural models for epilepsies (Bull. Math. Biophysics,33, 539–553, 1971;34, 71–78, 1972). The second model is developed in more detail and an attempt is made to bring it more in line with established neurological findings. The question of classification of some epilepsies is briefly discussed.     

Contributions to relational biology

The principle of biotopological mapping (Rashevsky, 1954,Bull. Math. Biophysics,16, 317–48) is given a generalized formulation, as the principle of relational epimorphism in biology. The connection between this principle and Robert Rosen’s representation of organisms by means of categories (1958,Bull. Math. Biophysics,20, 317–41) is studied. Rosen’s theory of (M,R)-systems, (1958,Bull. Math. Biophysics,20, 245–60) is generalized by dropping the assumption that only terminalM _i components are sending inputs into theR _i components. It is shown that, if the primordial organism is an (M,R)-system, then the higher organisms, obtained by a construction well discussed previously (1958,Bull. Math. Biophysics,20, 71–93), are also (M,R)-systems. Several theorems about such derived (M,R)-systems are demonstrated. It is shown that Rosen’s concept of an organism as a set of mappings throws light on phenomena of synesthesia and also leads to the conclusion that Gestalt phenomena must occur not only in the fields of visual and auditory perception but in perceptions of any modality.     

Analysis of chemotactic bacterial distributions in population migration assays using a mathematical model applicable to steep or shallow attractant gradients

The mathematical model developed by Riveroet al. (1989,Chem. Engng Sci. 44, 2881–2897) is applied to literature data measuring chemotactic bacterial population distributions in response to steep as well as shallow attractant gradients. This model is based on a fundamental picture of the sensing and response mechanisms of individual bacterial cells, and thus relates individual cell properties such as swimming speed and tumbling frequency to population parameters such as the random motility coefficient and the chemotactic sensitivity coefficient. Numerical solution of the model equations generates predicted bacterial density and attractant concentration profiles for any given experimental assay. We have previously validated the mathematical model from experimental work involving a step-change in the attractant gradient (Fordet al., 1991Biotechnol. Bioengng.37, 647–660; For and Lauffenburger, 1991,Biotechnol. Bioengng,37, 661–672). Within the context of this experimental assay, effects of attractant diffusion and consumption, random motility, and chemotactic sensitivity on the shape of the profiles are explored to enhance our understanding of this complex phenomenon. We have applied this model to various other types of gradients with successful intepretation of data reported by Dalquistet al. (1972,Nature New Biol. 236, 120–123) forSalmonella typhimurum validating the mathematical model and supportin the involvement of high and low affinity receptors for serine chemotaxis by these cells.     

Some considerations on relational equilibria

A previous study (Bull. Math. Biophysics,31, 417–427, 1969) on the definitions of stability of equilibria in organismic sets determined byQ relations is continued. An attempt is made to bring this definition into a form as similar as possible to that used in physical systems determined byF-relations. With examples taken from physics, biology and sociology, it is shown that a definition of equilibria forQ-relational systems similar to the definitions used in physics can be obtained, provided the concept of stable or unstable structures of a system determined byQ-relations is considered in a probabilistic manner. This offers an illustration of “fuzzy categories,” a notion introduced by I. Bąianu and M. Marinescu (Bull. Math. Biophysics,30, 625–635, 1968), in their paper on organismic supercategories, which is designed to provide a mathematical formalism for Rashevsky's theory of Organismic Sets (Bull. Math. Biophysics,29, 389–393, 1967;30, 163–174, 1968;31, 159–198, 1969). A suggestion is made for a method of mapping the abstract discrete space ofQ-relations on a continuum of variables ofF-relations. Problems of polymorphism and metamorphosis, both in biological and social organisms, are discussed in the light of the theory.     

Linearized oscillations in population dynamics

A linearized oscillation theorem due to Kulenović, Ladas and Meimaridou (1987,Quart. appl. Math. XLV, 155–164) and an extension of it are applied to obtain the oscillation of solutions of several equations which have appeared in population dynamics. They include the logistic equation with several delays, Nicholson's blowflies model as described by Gurney, Blythe and Nisbet (1980,Nature, Lond. 287, 17–21) and the Lasota-Wazewska model of the red blood cell supply in an animal. We also developed a linearized oscillation result for difference equations and applied it to several equations taken from the biological literature.     

Kinetic study of DNA binding of cisplatin and of a bicycloalkyl-substituted (ethylenediamine)dichloroplatinum(II) complex

 Salmon sperm DNA platination has been conducted under strictly pseudo-first-order conditions with cisplatin (1) and rac-{(1S,2S,4S)-exo-2-(aminomethyl)-2-amino-7-bicyclo[2.2.1]heptane}dichloroplatinum(II) (2). An aquation step first occurs for both complexes, with the rate constants k ₁ = 1.12(0.02)×10^–4 s^–1 and 1.47(0.02)×10^–4 s^–1 respectively for 1 and 2 at 37  °C, values in agreement with those previously reported. It is followed by the actual platination step whose second-order rate constant has been determined for the first time by physicochemical techniques. The values for 1 and 2 respectively are: k ₂ = 2.08(0.07) M^–1 s^–1 and 3.9(0.4) M^–1 s^–1. These kinetic data are discussed in the context of a comparison of several biological properties of the two complexes. Received: 15 May 1998 / Accepted: 26 June 1998     

On relations between sets: III

The theory of relations between sets, proposed and outlined in previous publications (Bull. Math. Biophysics,23, 233–235, 1961;28, 117–124, 1966;28, 309–313, 1966), is tentatively expanded and generalized with a view to biological applications.     

A note on the development of organismic sets

It is suggested that the development of organismic sets is governed not by the maximalization of the integral survival value, as suggested previously (Bull. Math. Biophysics,28, 283–308, 1966;29, 139–152, 1967;30, 163–174, 1968), but by maximizing the number of new relations which appear as an organismic set develops.     

On the number of balanced signed graphs

The classical enumeration theorem of Pólya (Acta Math.,68, 145–254, 1937) is applied to a modified version of Harary’s (Pacific J. Math.,8, 743–755, 1958) generating functions for counting bicolored graphs to derive a counting function for the number of balanced signed graphs. Methods for computing these counting polynomial functions are discussed.     

Sporadic sustained nonperiodic oscillations in the activities of organismic sets

In combining the author's theories of organismic sets (Rashevsky,Bull. Math. Biophysics,31, 159–198, 1969a) and Robert Rosen's theory of (M, R)-systems (Bull. Math. Biophysics,20, 245–265, 1958), a conclusion is reached that the number of either normal or pathological phenomena in organismic sets may occur. Those phenomena are characterized by occurring spontaneously once in a while but are not exactly periodic. Some epilepsies are an example of such pathological phenomena in the brain.     

Models for reserpine effects on normal and brain-operated animals on a continuous reinforcement bar pressing schedule

Mathematical models are presented for the relationship between changes in the bar pressing rate of albino rats on a continuous reinforcement schedule and alterations of the brain amines norepinephrine and serotonin induced by electrolytic lesions in the hypothalamus and by the drug, reserpine. The discussion utilizes central concepts from some of the more prominent hypotheses on the relation of the amines to behavior and for the potentiation effects of the drug by the lesion (A. Carlsson, 1964. “Functional Significance of Drug-induced Changes in Brain Monoamine Levels.” InBiogenic Amines, eds. H. E. Himwich and W. A. Himwich, New York: Elsevier, 9–27; B. B. Brodie and E. Costa, 1962. “Some Current Views on Brain Monoamines.”Psychopharm. Service Centre Bull.,2, 1–25.) Laboratory data (F. Grabarits and J. Harvey, 1966. “The Effects of Reserpine on Behavior and on Brain Concentrations of Serotonin and Norepinephrine in Control Rats and Rats with Hypothalamic Lesions.”J. Pharm.,153, 403–411) are interpreted, and suggestions for further experimental and theoretical work are made.     

Evolutionary tradeoff and equilibrium in an aquatic predator-prey system

Due to the conventional distinction between ecological (rapid) and evolutionary (slow) timescales, ecological and population models have typically ignored the effects of evolution. Yet the potential for rapid evolutionary change has been recently established and may be critical to understanding how populations persist in changing environments. In this paper we examine the relationship between ecological and evolutionary dynamics, focusing on a well-studied experimental aquatic predator-prey system (Fussmann et al., 2000, Science, 290, 1358–1360; Shertzer et al., 2002, J. Anim. Ecol., 71, 802–815; Yoshida et al., 2003, Nature, 424, 303–306). Major properties of predator-prey cycles in this system are determined by ongoing evolutionary dynamics in the prey population. Under some conditions, however, the populations tend to apparently stable steady-state densities. These are the subject of the present paper. We examine a previously developed model for the system, to determine how evolution shapes properties of the equilibria, in particular the number and identity of coexisting prey genotypes. We then apply these results to explore how evolutionary dynamics can shape the responses of the system to ‘management’: externally imposed alterations in conditions. Specifically, we compare the behavior of the system including evolutionary dynamics, with predictions that would be made if the potential for rapid evolutionary change is neglected. Finally, we posit some simple experiments to verify our prediction that evolution can have significant qualitative effects on observed population-level responses to changing conditions.     

Physics,biology and sociology: II. Suggestion for a synthesis

In continuation of previous studies (Bull. Math. Biophysics,28, 283–308; 655–661, 1966;29, 139–152, 1967) it is shown that the difference between the “metric” aspects of physics and the “relational” aspects of biological and social sciences disappear by accepting the broader definition of “relation”, such as that given in mathematics and logic. A conceptual superstructure then becomes possible from which all three branches of knowledge may be derived, though none of them can be derived from the others.     

Relational forces and structures produced by them

After giving a brief review of the theory of organismic sets (Bull. Math. Biophysics,29, 139–152, 1967;31, 159–198, 1969), in which the concept of relational forces, introduced earlier (Bull. Math. Biophysics,28, 283–308, 1966a) plays a fundamental role, the author discusses examples of possible different structures produced by relational forces. For biological organisms the different structures found theoretically are in general agreement with observation. For societies, which are also organismic sets as discussed in the above references, the structures can be described only in an abstract space, the nature of which is discussed. Different isomorphisms between anatomical structures, as described in ordinary Euclidean space, and the sociological structures described in an abstract space are noted, as should be expected from the theory of organismic sets.     

Human trachea primary epithelial cells express both sialyl(α2-3)Gal receptor for human parainfluenza virus type 1 and avian influenza viruses, and sialyl(α2-6)Gal receptor for human influenza viruses

We reported previously that the dominant receptors of influenza A and B viruses, and human and murine respiroviruses, were sialylglycoproteins and gangliosides containing monosialo-lactosamine type I-and II-residues, such as sialic acid-α2-3(6)-Galβ1-3(4)-GlcNAcβ1-. In addition, the Siaα2-3Gal linkage was predominantly recognized by avian and horse influenza viruses, and human parainfluenza virus type 1 (hPIV-1), whereas the Siaα2-6Gal linkage was mainly recognized by human influenza viruses (Paulson JC in “The Receptors' [Conn M Ed] 2, 131–219 (1985); Suzuki Y, Prog Lipid Res 33, 429–57 (1994); Ito T, J Virol 73, 6743–51 (2000); Suzuki Y, J Virol 74, 11825–31 (2000); Suzuki T, J. Virol 75, 4604–4613 (2001); Suzuki Y, Biol. Pharm. Bull. 28, 399–408 (2005)). To clarify the distribution of influenza virus receptors on the human bronchial epithelium cell surface, we investigated a primary culture of normal human bronchial epithelial (NHBE) cells using two types of lectin (MAA and SNA), which recognize sialyl linkages (α2-3 and α2-6), using fluorescence-activated cell-sorting analysis. The results showed that both α2-3- and α2-6-linked Sias were expressed on the surface of primary human bronchial epithelial cells. The cells infected by hPIV-1 bound to MAA, confirming that cells targeted by hPIV-1 have α2-3-linked oligosaccharides. We also compared the ability of hPIV-1 and human influenza A virus to infect primary human bronchial epithelial cells pre-treated with Siaα2-3Gal-specific sialidase from Salmonella typhimurium. No difference was observed in the number of sialidase pre-treated and non-treated cells infected with human influenza A virus, which binds to Siaα2-6Gal-linked oligosaccharides. By contrast, the number of cells infected with hPIV-1 decreased significantly upon sialidase treatment. Thus, cultured NHBE cells showed both α2-3-linked Sias recognized by hPIV-1 and avian influenza virus receptors, and α2-6-linked Sias recognized by human influenza virus receptors.     

Recent Advances in GFP Folding Reporter and Split-GFP Solubility Reporter Technologies. Application to Improving the Folding and Solubility of Recalcitrant Proteins from Mycobacterium tuberculosis

We have improved our green fluorescent protein (GFP) folding reporter technology [Waldo et al., (1999) Nat. Biotechnol. 17, 691–695] to evolve recalcitrant proteins from Mycobacterium tuberculosis. The target protein is inserted into the scaffolding of the GFP, eliminating false-positive artifacts caused by expression of truncated protein variants from internal cryptic ribosome binding sites in the target RNA. In parallel, we have developed a new quantitative fluorescent protein tagging and detection system based on micro-domains of GFP. This split-GFP system, which works both in vivo and in vitro, is amenable to high-throughput assays of protein expression and solubility [Cabantous et al., (2005) Nat. Biotechnol. 23, 102–107]. Together, the GFP folding reporter and split-GFP technologies offer a comprehensive system for manipulating and improving protein folding and solubility.     

Density functional theory study on (LiNH<Subscript>2</Subscript>)<Subscript>n</Subscript> (<Emphasis Type="Italic">n</Emphasis> = 1–5) clusters

Geometrical structures and relative stabilities of (LiNH₂)_n (n = 1–5) clusters were studied using density functional theory (DFT) at the B3LYP/6-31G* and B3LYP/6-31++G* levels. The electronic structures, vibrational properties, N–H bond dissociation energies (BDE), thermodynamic properties, bond properties and ionization potentials were analyzed for the most stable isomers. The calculated results show that the Li–N and Li–Li bonds can be formed more easily than those of the Li–H or N–H bonds in the clusters, in which NH₂ is bound to the framework of Li atomic clusters with fused rings. The average binding energies for each LiNH₂ unit increase gradually from 142 kJ mol⁻¹ up to about 180 kJ mol⁻¹ with increasing n. Natural bond orbital (NBO) analysis suggests that the bonds between Li and NH₂ are of strong ionicity. Three-center–two-electron Li–N–Li bonding exists in the (LiNH₂)₂ dimer. The N–H BDE values indicate that the change in N–H BDE values from the monomer a1 to the singlet-state clusters is small. The N–H bonds in singlet state clusters are stable, while the N–H bonds in triplet clusters dissociate easily. A study of their thermodynamic properties suggests that monomer a1 forms clusters (b1, c1, d2 and e1) easily at low temperature, and clusters with fewer numbers of rings tend to transfer to ones with more rings at low temperature. E _g, E _HOMO and E _av decrease gradually, and become constant. Ring-like (LiNH₂)_3,4 clusters possess higher ionization energy (VIE) and E _g, but lower values of E _HOMO. Ring-like (LiNH₂)_3,4 clusters are more stable than other types. A comparison of structures and spectra between clusters and crystal showed that the NH₂ moiety in clusters has a structure and spectral features similar to those of the crystal.