Structured population on two patches: modeling dispersal and delay

We derive from the age-structured model a system of delay differential equations to describe the interaction of spatial dispersal (over two patches) and time delay (arising from the maturation period). Our model analysis shows that varying the immature death rate can alter the behavior of the homogeneous equilibria, leading to transient oscillations around an intermediate equilibrium and complicated dynamics (in the form of the coexistence of possibly stable synchronized periodic oscillations and unstable phase-locked oscillations) near the largest equilibrium.     

Oscillations in a refractory neural net

A functional differential equation that arises from the classic theory of neural networks is considered. As the length of the absolute refractory period is varied, there is, as shown here, a super-critical Hopf bifurcation. As the ratio of the refractory period to the time constant of the network increases, a novel relaxation oscillation occurs. Some approximations are made and the period of this oscillation is computed.     

Inverse bifurcation problem, singular Wiener-Hopf equations, and mathematical models in ecology

A single-species population dynamics with dispersal in a spatially heterogeneous environment is modeled by a nonlinear reaction-diffusion equation with a potential term. To each nonlinear kinetics there corresponds a bifurcation curve that describes the relation between the growth rate and the central density of a steady-state population distribution. Our main concern is an inverse problem for this correspondence. The existence of nonlinear kinetics realizing a prescribed bifurcation curve is established. It is shown that the freedom of such kinetics is of degree finite and even, depending only on the heterogeneity of the environment, and conversely that any nonnegative even integer occurs as the degree of freedom in some environments. A discussion is also made on under what kind of environment the degree is equal to zero or is positive. The mathematical analysis involves the development of a general theory for singular multiplicative Wiener-Hopf integral equations.     

Mathematical analysis of delay differential equation models of HIV-1 infection

Models of HIV-1 infection that include intracellular delays are more accurate representations of the biology and change the estimated values of kinetic parameters when compared to models without delays. We develop and analyze a set of models that include intracellular delays, combination antiretroviral therapy, and the dynamics of both infected and uninfected T cells. We show that when the drug efficacy is less than perfect the estimated value of the loss rate of productively infected T cells, delta, is increased when data is fit with delay models compared to the values estimated with a non-delay model. We provide a mathematical justification for this increased value of delta. We also provide some general results on the stability of non-linear delay differential equation infection models.     

The effect of dispersal on single-species nonautonomous dispersal models with delays

In this paper, single-species nonautonomous dispersal models with delays are considered. An interesting result on the effect of dispersal for persistence and extinction is obtained. That is, if the species is persistent in a patch then it is also persistent in all other patches; if the species is permanent in a patch then it is also permanent in all other patches; if the species is extinct in a patch then it is also extinct in all other patches. Furthermore, some new sufficient conditions for the permanence and extinction of the species in a patch are established. The existence of positive periodic solutions is obtained in the periodic case by employing Teng and Chen's results on the existence of positive periodic solutions for functional differential equations. Received: 26 June 2000 / Revised version: 6 October 2000 / Published online: 10 April 2001     

Global analysis of the Michaelis-Menten-type ratio-dependent predator-prey system

The recent broad interest on ratio-dependent based predator functional response calls for detailed qualitative study on ratio-dependent predator-prey differential systems. A first such attempt is documented in the recent work of Kuang and Beretta(1998), where Michaelis-Menten-type ratio-dependent model is studied systematically. Their paper, while contains many new and significant results, is far from complete in answering the many subtle mathematical questions on the global qualitative behavior of solutions of the model. Indeed, many of such important open questions are mentioned in the discussion section of their paper. Through a simple change of variable, we transform the Michaelis-Menten-type ratio-dependent model to a better studied Gause-type predator-prey system. As a result, we can obtain a complete classification of the asymptotic behavior of the solutions of the Michaelis-Menten-type ratio-dependent model. In some cases we can determine how the outcomes depend on the initial conditions. In particular, open questions on the global stability of all equilibria in various cases and the uniqueness of limit cycles are resolved. Biological implications of our results are also presented.     

An alternative stochastic model of generation of oligodendrocytes in cell culture

According to our previous model, oligodendrocyte – type 2 (O-2A) astrocyte progenitor cells become competent for differentiation in vitro after they complete a certain number of critical mitotic cycles. After attaining the competency to differentiate, progenitor cells divide with fixed probability p in subsequent cycles. The number of critical cycles is random; analysis of data suggests that it varies from zero to two. The present paper presents an alternative model in which there are no critical cycles, and the probability that a progenitor cell will divide again decreases gradually to a plateau value as the number of completed mitotic cycles increases. In particular all progenitor cells have the ability to differentiate from the time of plating. The Kiefer-Wolfowitz procedure is used to fit the new model to experimental data on the clonal growth of purified O-2A progenitor cells obtained from the optic nerves of 7 day old rats. The new model is shown to fit the experimental data well, indicating that it is not possible to determine whether critical cycles exist on the basis of these experimental data. In contrast to the fit of the previous model, which suggested that the addition of thyroid hormone increased the limiting probability of differentiation as the number of mitotic cycles increases, the fit of the new model suggests that the addition of thyroid hormone has almost no effect on the limiting probability of differentiation. Received: 6 March 2000 / Revised version: 18 September 2000 / Published online: 30 April 2001     

Equilibrium structure and stability in a frequency-dependent,two-population diploid model

We investigate the equilibrium structure for an evolutionary genetic model in discrete time involving two monoecious populations subject to intraspecific and interspecific random pairwise interactions. A characterization for local stability of an equilibrium is found, related to the proximity of this equilibrium with evolutionarily stable strategies (ESS). This extends to a multi-population framework a principle initially proposed for single populations, which states that the mean population strategy at a locally stable equilibrium is as close as possible to an ESS.     

Parametric analysis of the ratio-dependent predator–prey model

We present a complete parametric analysis of stability properties and dynamic regimes of an ODE model in which the functional response is a function of the ratio of prey and predator abundances. We show the existence of eight qualitatively different types of system behaviors realized for various parameter values. In particular, there exist areas of coexistence (which may be steady or oscillating), areas in which both populations become extinct, and areas of "conditional coexistence" depending on the initial values. One of the main mathematical features of ratio-dependent models, distinguishing this class from other predator-prey models, is that the Origin is a complicated equilibrium point, whose characteristics crucially determine the main properties of the model. This is the first demonstration of this phenomenon in an ecological model. The model is investigated with methods of the qualitative theory of ODEs and the theory of bifurcations. The biological relevance of the mathematical results is discussed both regarding conservation issues (for which coexistence is desired) and biological control (for which extinction is desired).     

The influence of drug treatment on the maintenance of schistosome genetic diversity

Drug treatment of patients with schistosomiasis may select for drug-resistant parasites. In this article, we formulate a deterministic model with multiple strains of schistosomes (helminth parasites with a two-host life cycles) in order to explore the role of drug treatment in the maintenance of a polymorphism of parasite strains that differ in their resistance levels. The basic reproductive numbers for all strains are computed, and are shown to determine the stabilities of equilibria of the model and consequently the distribution of parasite phenotypes with different levels of drug tolerance. Analysis of our model shows that the likelihood that resistant strains will increase in frequency depends on the interplay between their relative fitness, the cost of resistance, and the degree of selection pressure exerted by the drug treatments.     

Persistent solutions for age-dependent pair-formation models

Conditions on the vital rates and the mating function are derived which imply existence or nonexistence of exponentially growing persistent age-distributions for age-dependent pair-formation models.     

Models of central pattern generators for quadruped locomotion I. Primary gaits

In this paper we continue the analysis of a network of symmetrically coupled cells modeling central pattern generators for quadruped locomotion proposed by Golubitsky, Stewart, Buono, and Collins. By a cell we mean a system of ordinary differential equations and by a coupled cell system we mean a network of identical cells with coupling terms. We have three main results in this paper. First, we show that the proposed network is the simplest one modeling the common quadruped gaits of walk, trot, and pace. In doing so we prove a general theorem classifying spatio-temporal symmetries of periodic solutions to equivariant systems of differential equations. We also specialize this theorem to coupled cell systems. Second, this paper focuses on primary gaits; that is, gaits that are modeled by output signals from the central pattern generator where each cell emits the same waveform along with exact phase shifts between cells. Our previous work showed that the network is capable of producing six primary gaits. Here, we show that under mild assumptions on the cells and the coupling of the network, primary gaits can be produced from Hopf bifurcation by varying only coupling strengths of the network. Third, we discuss the stability of primary gaits and exhibit these solutions by performing numerical simulations using the dimensionless Morris-Lecar equations for the cell dynamics.     

A new explicit stability criterion for human periodic breathing<--RID="*"--><--ID="*" This research was supported by grants from DRED (96-920).-->

The aim of this paper is to carry out a stability analysis for periodic breathing in humans that incorporates the dynamic characteristics of ventilation control. A simple CO₂ model that takes into account the main elements of the respiratory system, i.e. the lungs and the ventilatory controller with its dynamic properties, is presented. This model results in a three-dimensional non-linear delay differential system for which there exists a unique equilibrium point. Our stability analysis of this equilibrium point leads to the definition of a new explicit stability criterion and to the demonstration of the existence of a Hopf bifurcation. Numerical simulations illustrate the influence of physiological parameters on the stability of ventilation, and particularly the major role of the dynamic characteristics of the respiratory controller. Received: 2 February 1999 / Revised version: 18 June 1999 / Published online: 23 October 2000     

Mathematical modelling of the intravenous glucose tolerance test

 Several attempts at building a satisfactory model of the glucose-insulin system are recorded in the literature. The minimal model, which is the model currently mostly used in physiological research on the metabolism of glucose, was proposed in the early eighties for the interpretation of the glucose and insulin plasma concentrations following the intravenous glucose tolerance test. It is composed of two parts: the first consists of two differential equations and describes the glucose plasma concentration time-course treating insulin plasma concentration as a known forcing function; the second consists of a single equation and describes the time course of plasma insulin concentration treating glucose plasma concentration as a known forcing function. The two parts are to be separately estimated on the available data. In order to study glucose-insulin homeostasis as a single dynamical system, a unified model would be desirable. To this end, the simple coupling of the original two parts of the minimal model is not appropriate, since it can be shown that, for commonly observed combinations of parameter values, the coupled model would not admit an equilibrium and the concentration of active insulin in the “distant” compartment would be predicted to increase without bounds. For comparison, a simple delay-differential model is introduced, is demonstrated to be globally asymptotically stable around a unique equilibrium point corresponding to the pre-bolus conditions, and is shown to have positive and bounded solutions for all times. The results of fitting the delay-differential model to experimental data from ten healthy volunteers are also shown. It is concluded that a global unified model is both theoretically desirable and practically usable, and that any such model ought to undergo formal analysis to establish its appropriateness and to exclude conflicts with accepted physiological notions. Received: 22 June 1998 / Revised version: 24 February 1999     

A mathematical model for drug administration by using the phagocytosis of red blood cells

 A mathematical model for the delivery of drug directly to the macrophages by using the phagocytosis of senescent red blood cells is proposed. The model is based on the following assumption: At time t=0 a preassigned red blood cell population n(0, a)=φ(a), a>0, loaded by the drug, is injected in the blood circulation. Among the cells of that population only those with an age a≧ā (ā=120 days) will be phagocytosed by macrophages. Of course, the lifetime of the drug must be higher than ā. Within the red blood cells it cannot be metabolized, neither can it diffuse through their membranes. The emphasis of the paper is on the mathematical properties and on the formulation of the control problem. Received 15 December 1994; received in revised form 20 July 1995     

Approximating the long-term behaviour of a model for parasitic infection

In a companion paper two stochastic models, useful for the initial behaviour of a parasitic infection, were introduced. Now we analyse the long term behaviour. First a law of large numbers is proved which allows us to analyse the deterministic analogues of the stochastic models. The behaviour of the deterministic models is analogous to the stochastic models in that again three basic reproduction ratios are necessary to fully describe the information needed to separate growth from extinction. The existence of stationary solutions is shown in the deterministic models, which can be used as a justification for simulation of quasi-equilibria in the stochastic models. Host-mortality is included in all models. The proofs involve martingale and coupling methods.     

Metabolite analysis in positron emission tomography studies: examples from food sciences

Summary. Substances of various chemical structures can be labelled with appropriate positron emitting isotopes and applied as tracer compounds in PET examinations. Using dynamic data acquisition protocols, time-activity curves of radioactivity uptake in organs can be derived and the measurements of tissue tracer concentrations can be translated into quantitative values of tissue function. However, analysis of metabolites of these tracers regarding their nature and distribution in the living organism is an essential need for the quantitative analysis of PET measurements. In addition, metabolite analysis contributes to the interpretation of the images obtained as well as to the identification of pathological changes in metabolic pathways. This paper reports on representative examples of radiolabelled compounds which might be of importance in food science (e.g., amino acids, polyphenols, and model compounds for advanced glycation end products (AGEs)). Typical procedures of analysis (radio-HPLC, radio-TLC) including pre-analytical sample preparation are described. Specific challenges of the method, e.g., trace amounts of radiolabelled compounds and the influence of the often very short half-lives of positron-emitting nuclides used are highlighted. Representative results of analyses of plasma, urine, and tissue samples are presented and discussed in terms of the metabolic fate of the tracers.     

Stochastic models of a parasitic infection, exhibiting three basic reproduction ratios

Two closely related stochastic models of parasitic infection are investigated: a non-linear model, where density dependent constraints are included, and a linear model appropriate to the initial behaviour of an epidemic. Host-mortality is included in both models. These models are appropriate to transmission between homogeneously mixing hosts, where the amount of infection which is transferred from one host to another at a single contact depends on the number of parasites in the infecting host. In both models, the basic reproduction ratio R0 can be defined to be the lifetime expected number of offspring of an adult parasite under ideal conditions, but it does not necessarily contain the information needed to separate growth from extinction of infection. In fact we find three regions for a certain parameter where different combinations of parameters determine the behavior of the models. The proofs involve martingale and coupling methods.     

Stability analysis of the partial selfing selection model

We undertake a detailed study of the one-locus two-allele partial selfing selection model. We show that a polymorphic equilibrium can exist only in the cases of overdominance and underdominance and only for a certain range of selfing rates. Furthermore, when it exists, we show that the polymorphic equilibrium is unique. The local stability of the polymorphic equilibrium is investigated and exact analytical conditions are presented. We also carry out an analysis of local stability of the fixation states and then conclude that only overdominance can maintain polymorphism in the population. When the linear local analysis is inconclusive, a quadratic analysis is performed. For some sets of selective values, we demonstrate global convergence. Finally, we compare and discuss results under the partial selfing model and the random mating model.     

Near infrared spectroscopy, cluster and multivariate analysis hyphenated to thin layer chromatography for the analysis of amino acids

Summary. A method based on near-infrared spectroscopy (NIRS) was developed for the rapid and non-destructive determination and quantification of solid and dissolved amino acids. The statistical results obtained after optimisation of measurement conditions were evaluated on the basis of statistical parameters, Q-value (quality of calibrations), R², standard error of estimation (SEE), standard error of prediction (SEP), BIAS applying cluster and different multivariate analytical procedures. Experimental optimisation comprised the selection of the highest suitable optical thin-layer (0.5, 1.0, 1.5, 2.0, 2.5, 3.0 mm), sample temperature (10–30 °C), measurement option (light fibre, 0.5 mm optical thin-layer; boiling point tube; different types of cuvettes) and sample concentration in the range between 100 and 500 ppm. Applying the optimised conditions and a 115-QS Suprasil^? cuvette (V = 400 μl), the established qualitative model enabled to distinguish between different dissolved amino acids with a Q-value of 0.9555. Solid amino acids were investigated in the transflectance mode, allowing to differentiate them with a Q-value of 0.9155. For the qualitative and quantitative analysis of amino acids in complex matrices NIRS was established as a detection system directly onto the plate after prior separation on cellulose based thin-layer chromatography (TLC) sheets employing n-butanol, acetic acid and distilled water at a ratio of 8:4:2 (v/v/v) as an optimised mobile phase. Due to the prior separation step, the established calibration curve was found to be more stable than the one calculated from the dissolved amino acids. The found lower limit of detection was 0.01 mg/ml. Finally, this optimised TLC-NIRS method was successfully applied for the qualitative and quantitative analysis of L-lysine in apple juice. NIRS is shown not only to offer a fast, non-destructive detection tool but also to provide an easy-to-use alternative to more complicated detection methods such as mass spectrometry (MS) for qualitative and quantitative TLC analysis of amino acids in crude samples.