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.     

Mathematical modeling of capillary formation and development in tumor angiogenesis: penetration into the stroma

The purpose of this paper is to present a mathematical model for the tumor vascularization theory of tumor growth proposed by Judah Folkman in the early 1970s and subsequently established experimentally by him and his coworkers [Ausprunk, D. H. and J. Folkman (1977) Migration and proliferation of endothelial cells in performed and newly formed blood vessels during tumor angiogenesis, Microvasc Res., 14, 53–65; Brem, S., B. A. Preis, ScD. Langer, B. A. Brem and J. Folkman (1997) Inhibition of neovascularization by an extract derived from vitreous Am. J. Opthalmol., 84, 323–328; Folkman, J. (1976) The vascularization of tumors, Sci. Am., 234, 58–64; Gimbrone, M. A. Jr, R. S. Cotran, S. B. Leapman and J. Folkman (1974) Tumor growth and neovascularization: an experimental model using the rabbit cornea, J. Nat. Cancer Inst., 52, 413–419]. In the simplest version of this model, an avascular tumor secretes a tumor growth factor (TGF) which is transported across an extracellular matrix (ECM) to a neighboring vasculature where it stimulates endothelial cells to produce a protease that acts as a catalyst to degrade the fibronectin of the capillary wall and the ECM. The endothelial cells then move up the TGF gradient back to the tumor, proliferating and forming a new capillary network. In the model presented here, we include two mechanisms for the action of angiostatin. In the first mechanism, substantiated experimentally, the angiostatin acts as a protease inhibitor. A second mechanism for the production of protease inhibitor from angiostatin by endothelial cells is proposed to be of Michaelis-Menten type. Mathematically, this mechanism includes the former as a subcase. Our model is different from other attempts to model the process of tumor angiogenesis in that it focuses (1) on the biochemistry of the process at the level of the cell; (2) the movement of the cells is based on the theory of reinforced random walks; (3) standard transport equations for the diffusion of molecular species in porous media. One consequence of our numerical simulations is that we obtain very good computational agreement with the time of the onset of vascularization and the rate of capillary tip growth observed in rabbit cornea experiments [Ausprunk, D. H. and J. Folkman (1977) Migration and proliferation of endothelial cells in performed and newly formed blood vessels during tumor angiogenesis, Microvasc Res., 14, 73–65; Brem, S., B. A. Preis, ScD. Langer, B. A. Brem and J. Folkman (1997) Inhibition of neovascularization by an extract derived from vitreous Am. J. Opthalmol., 84, 323–328; Folkman, J. (1976) The vascularization of tumors, Sci. Am., 234, 58–64; Gimbrone, M. A. Jr, R. S. Cotran, S. B. Leapman and J. Folkman (1974) Tumor growth and neovascularization: An experimental model using the rabbit cornea, J. Nat. Cancer Inst., 52, 413–419]. Furthermore, our numerical experiments agree with the observation that the tip of a growing capillary accelerates as it approaches the tumor [Folkman, J. (1976) The vascularization of tumors, Sci. Am., 234, 58–64]. An erratum to this article is available at .     

Homeostatic control of thyroxin concentration expressed by a set of linear differential equations

The purpose of this work is to express current concepts on the relationship between the rates of secretion of thyroxin and of thyroid stimulating hormone (TSH) by a set of linear differential equations (two attempts have been made previously in this direction; cf. Roston,Bull. Math. Biophysics,21, 271–282, 1959; Danziger and Elmergreen,Bull. Math. Biophysics,16, 15–21, 1954), and to show that the solutions to these equations fulfill two criteria: that they correctly express the previously observed behavior of thyroxin and TSH, and that they allow certain predictions to be made which are amenable to experimental verification or disproval by currently existing techniques. This mathematical model is necessarily only an approximation of reality.     

Chemotaxis and chemokinesis in eukaryotic cells: The Keller-Segel equations as an approximation to a detailed model

More than 20 years after its proposal, Keller and Segel's model (1971,J. theor. Biol.,30, 235–248) remains by far the most popular model for chemical control of cell movement. However, before the Keller-Segel equations can be applied to a particular system, appropriate functional forms must be specified for the dependence on chemical concentration of the cell transport coefficients and the chemical degradation rate. In the vast majority of applications, these functional forms have been chosen using simple intuitive criteria. We focus on the particular case of eukaryotic cell movement, and derive an approximation to the detailed model of Sherrattet al. (1993,J. theor. Biol.,162, 23–40). The approximation consists of the Keller-Segel equations, with specific forms predicted for the cell transport coefficients and chemical degradation rate. Moreover, the parameter values in these functional forms can be directly measured experimentally. In the case of the much studied neutrophil-peptide system, we test our approximation using both the Boyden chamber and under-agarose assays. Finally, we show that for other cell-chemical interactions, a simple comparison of time scales provides a rapid check on the validity of our Keller-Segel approximation.     

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.     

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.     

A neural model for conditioned avoidance learning

A mathematical model for learning of a conditioned avoidance behavior is presented. An identification of the net excitation of a neural model (Rashevsky, N., 1960.Mathematical Biophysics. Vol. II. New York: Dover Publications, Inc.) with the instantaneous probability of response is introduced and its usefulness in discussing block-trial learning performances in the conditioned avoidance situation is outlined for normal and brain-operated animals, using experimental data collected by the author. Later, the model is applied to consecutive trial learning and connection is made with the approach of H. D. Landahl (1964. “An Avoidance Learning Situation. A Neural Net Model.”Bull. Math. Biophysics,26, 83–89; and 1965, “A Neural Net Model for Escape Learning.”Bull. Math. Biophysics,27, Special Edition, 317–328) wherein lie further data with which the model can be compared.     

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.     

A suggestion of a new approach to the theory of some biological periodicities

Previous studies of L. Danziger and G. Elmergreen (Bull. Math. Biophysics,16, 15–21, 1954;18, 1–13, 1956) of possible biochemical periodicities in organisms assumed non-linear biochemical interaction between different metabolites, because linear systems do not lead to undamped ocsillations. They treated homogeneous systems. Later N. Rashevsky generalized their results to a more realistic case where the non-homogeneity due to the histological structure is considered. (Some Medical Aspects of Mathematical Biology, Springfield, Illinois: Charles C. Thomas, Publisher, 1964;Bull. Math. Biophysics,29, 389–393, 1967.) As long as the histological structure remains constant, the existence of sustained periodicities requires the assumption of non-linearity of biochemical interactions. If, however, the secretions of an endocrine gland affect the histological structure of the target organ, notably as in the menstrual cycle, and if there is a feed-back, the equations become non-linear and may admit sustained periodic solutions even if the purely biochemical interactions are linear.     

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.     

Modeling of recombinant bacteria fermentation for enhanced productivity

Accurate kinetic models are necessary in order to scale biochemical processes which utilize recombinant organisms. In this work, the performance of a kinetic model proposed by Miao and Kompala (Biotechnol. Bioeng.,40, 787–796, 1992) has been tested against experimental data at different concentrations of substrates. The experimental results from the induction of the batch cultures with IPTG (isopropyl-β-D-thiogalactopyranoside) prove that the model equations are reliable in predicting the biomass and foreign protein concentrations as well as the optimum induction time.     

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     

Mathematical theory of biological periodicities: Formulation of then-body case

In a series of papers, L. Danziger and G. Elmergreen (Bull. Math. Biophysics,16, 15–21, 1954;18, 1–13, 1956;19, 9–18, 1957) showed that a non-linear biochemical interaction between the anterior pituary gland and the thyroid gland may result under certain conditions in sustained periodical oscillations of the rates of production and of the blood level of the thyrotropic and of the thyroid hormone. They treated the systems, however, as a homogeneous one. N. Rashevsky (Some Medical Aspects of Mathematical Biology, Springfield, Illinois: Charles C. Thomas, Publisher, 1965;Bull. Math. Biophysics,29, 395–401, 1967) generalized the above results by taking into account the histological structures of the two glands as well as the diffusion coefficients and permeabilities of cells involved. The present paper is the first step toward the theory of interaction of any numbern of glands or, more generally,n components. The differential equations which govern the behavior of such a system represent a system of2n ²+n non-linear first order ordinary equations and involve a total of 7 n ²+3n parameters of partly histological, partly biochemical nature. The requirements of the existence of sustained oscillations demand 4n ²+2n+2 inequalities between those 7n ²+3n parameters.     

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.     

Gender-dependent effects of selenite on the perfused rat heart

Gender differences are related to the manner in which the heart responds to chronic and acute stress conditions of physiological and pathological nature. Depending on dose, sodium selenite acts as an antioxidant proven to have beneficial effects in several pathological conditions G. Drasch, J. Schopfer, and G. N. Schrauzer, Selenium/cadmium ratios in human prostates: indicators of prostate cancer risk of smokers and non-smokers, and relevance to the cancer protective effects of selenium,Biol. Trace Element Res. 103(2), 103–107 (2005); R. G. Kasseroller and G. N. Schrauzer, Treatment of secondary lymphedema of the arm with physical decongestive therapy and sodium selenite: a review,Am. J. Ther. 7(4), 273–279 (2000); G. N. Schrauzer, Anticarcinogenic effects of selenium,Cell. Mol. Life Sci. 57(13–14), 1864–1873 (2000); I. S. Palmer and O. E. Olson, Relative toxicities of selenite and selenate in the drinking water of rats,J. Nutr. 104(3), 306–314 (1974). To date, little is known about the gender-dependent direct effects of toxic doses of selenite on electrophysiology of the cardiovascular system H. A. Schroeder and M. Mitchener, Selenium and tellurium in rats: effect on growth, survival and tumors,J. Nutr. 101(11), 1531–1540 (1971); G. N. Schrauzer, The nutritional significance, metabolism and toxicology of selenomethionine,Adv. Food Nutr. Res. 47, 73–112 (2003). In the present study, the effects of in vitro toxic concentrations of sodium selenite ranging from 10^-6 M to 10^-3 M were tested on both male and female rat heart preparations. The toxic effects seen in an electrocardiogram and left ventricular pressure were dose and sex dependent at most of the tested concentrations. The present study clearly shows that at toxic doses, stress conditions are induced by selenite, resulting in gender-dependent modifications of the heart function. This modification is more pronounced in the contraction cascade of female rats. Males, on the other hand, had been much more affected in excitation-related parameters.     

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.     

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.     

Erythrocyte distribution in arterial blood flow: I. Basic equations and anisotropic effects

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

Comparison of surface potentials due to several singularity representations of the human heart

In electrocardiography the electrical potentials due to the heart actions can be analyzed by assuming the human body to be a conductor of homogeneous medium and the heart to be a combination of singularities within it. For a spherical conductor the “interior sphere theorem” of G. Ludford, J. Martinek, and G. Yeh (Proc. Cambridge Phil. Soc.,51, 389–93, 1955) renders potential expressions due to any singularity. For a conductor of prolate spheroidal shape the potential expressions due to a source-sink pair and a general dipole have been given by J. R. Wait (Jour. App. Physics,24, 496–97, 1953) and the authors (paper at the Conference on the Electrophysiology of the Heart, Feb. 16–17, 1956, in New York, to appear in theAnn. N. Y. Acad. Sciences) respectively. (A theorem which applies to any singularity inside a prolate or oblate spheroid will be published by the authors soon). This paper presents numerical and graphical results of potentials on the surfaces of a prolate spheroid and a sphere due to source-sink pairs and dipoles of several locations and directions and compares the various representations. A discussion on the judicious choice of heart models concludes the paper. This investigation was supported by The National Heart Institute under a research grant H-2263.     

A mathematical model for the inhibition of the multidrug resistance-associated P-glycoprotein pump

An extension of an earlier model of the p170 glycoprotein pump is presented. In an earlier work (Michelson and Slate,Bull. math. Biol. 54, 1023–1038, 1992), the pump was modeled using an energy-dependent model of facilitated diffusion. In this paper we add an inhibitor to the model. New equations are derived which represent either competitive or non-competitive inhibition in the pumping action of the glycoprotein. Numerical simulations were run which provide a response surface (initial loading concentration of inhibitor and its ability to compete with an ideal anti-cancer drug vs a summary measure of cytoplasmic exposure) for each scenario. The importance of the exposure profile, how it is related to ultimate tumor cell survival, and the binding requirements for developing multidrug resistance inhibitors are discussed.