Traveling waves in a model of influenza A drift

Between major pandemics, the influenza A virus changes its antigenic properties by accumulating point mutations (drift) mainly in the RNA genes that code for the surface proteins hemagglutinin (HA) and neuraminidase (NA). The successful strain (variant) that will cause the next epidemic is selected from a reduced number of progenies that possess relatively high transmissibility and the ability to escape from the immune surveillance of the host. In this paper, we analyse a one-dimensional model of influenza A drift (Z. Angew. Math. Mech. 76 (2) (1996) 421) that generalizes the classical SIR model by including mutation as a diffusion process in a phenotype space of variants. The model exhibits traveling wave solutions with an asymptotic wave speed that matches well those obtained from numerical simulations. As exact solutions for these waves are not available, asymptotic estimates for the amplitudes of infected and recovered classes are provided through an exponential approximation based on the smallness of the diffusion constant. Through this approximation, we find simple scaling properties to several parameters of relevance to the epidemiology of the disease.     

Traveling waves in the discrete fast buffered bistable system

We study the existence and uniqueness of traveling wave solutions of the discrete buffered bistable equation. Buffered excitable systems are used to model, among other things, the propagation of waves of increased calcium concentration, and discrete models are often used to describe the propagation of such waves across multiple cells. We derive necessary conditions for the existence of waves, and, under some restrictive technical assumptions, we derive sufficient conditions. When the wave exists it is unique and stable.      

Traveling Waves and Spread Rates for a West Nile Virus Model

A reaction–diffusion model for the spatial spread of West Nile virus is developed and analysed. Infection dynamics are based on a modified version of a model for cross infection between birds and mosquitoes (Wonham et al., 2004, An epidemiological model for West-Nile virus: Invasion analysis and control application. Proc. R. Soc. Lond. B 271), and diffusion terms describe movement of birds and mosquitoes. Working with a simplified version of the model, the cooperative nature of cross-infection dynamics is utilized to prove the existence of traveling waves and to calculate the spatial spread rate of infection. Comparison theorem results are used to show that the spread rate of the simplified model may provide an upper bound for the spread rate of a more realistic and complex version of the model.     

The hypercycle,traveling waves,and Wright's equation

A formal relation between the hypercycle equation and the delay differential equation of E. M. Wright is exhibited using a traveling waves approach. Several unsolved questions in either problem can be related and interpreted, in particular new motivation for the study of Wright's equation is obtained.     

An asymptotic solution for traveling waves of a nonlinear-diffusion Fisher's equation

We examine traveling-wave solutions for a generalized nonlinear-diffusion Fisher equation studied by Hayes [J. Math. Biol. 29, 531–537 (1991)]. The density-dependent diffusion coefficient used is motivated by certain polymer diffusion and population dispersal problems. Approximate solutions are constructed using asymptotic expansions. We find that the solution will have a corner layer (a shock in the derivative) as the diffusion coefficient approaches a step function. The corner layer at z = 0 is matched to an outer solution for z < 0 and a boundary layer for z > 0 to produce a complete solution. We show that this model also admits a new class of nonphysical solutions and obtain conditions that restrict the set of valid traveling-wave solutions. Supported by a National Science Foundation graduate fellowship. This work was performed under National Science Foundation grant DMS-9024963 and Air Force Office of Scientific Research grant AFOSR-F49620-94-1-0044.     

Traveling wave solutions of a reaction diffusion model for competing pioneer and climax species

Presented is a reaction-diffusion model for the interaction of pioneer and climax species. For certain parameters the system exhibits bistability and traveling wave solutions. Specifically, we show that when the climax species diffuses at a slow rate there are traveling wave solutions which correspond to extinction waves of either the pioneer or climax species. A leading order analysis is used in the one-dimensional spatial case to estimate the wave speed sign that determines which species becomes extinct. Results of these analyses are then compared to numerical simulations of wave front propagation for the model on one and two-dimensional spatial domains. A simple mechanism for harvesting is also introduced.     

Traveling waves in a piecewise-linear reaction-diffusion model of excitable medium

One-dimensional autowaves (traveling waves) in excitable medium described by a piecewise-linear reaction-diffusion system have been investigated. Two main types of wave have been considered: a single pulse and a periodic sequence of pulses (wave trains). In a two-component system, oscillations are due to the second component of the reaction-diffusion system, while in a one-component system, they are caused by external periodic excitation (forcing). Using semianalytical solutions for the wave profile, the shape and velocity of autowaves have been found. It is shown that the dispersion relation for oscillating sequences of pulses has an anomalous character.     

Traveling wave solutions for a one-dimensional crawling nematode sperm cell model

In this paper, we proved that the one-dimensional crawling nematode sperm cell model proposed by Mogilner and Verzi (2003) supports traveling wave solutions if there is no disassembly of unbundled filaments in the cell. Uniqueness of traveling wave is established under additional assumptions and numerical examples are also given in the paper. Mathematical methods used include dynamical system techniques, implicit function theorem and global bifurcation theory.Revised version: 16 September 2003     

A chemotactic model of trunk neural crest cell migration

Trunk neural crest cells follow a common ventral migratory pathway but are distributed into two distinct locations to form discrete sympathetic and dorsal root ganglia along the vertebrate axis. Although fluorescent cell labeling and time‐lapse studies have recorded complex trunk neural crest cell migratory behaviors, the signals that underlie this dynamic patterning remain unclear. The absence of molecular information has led to a number of mechanistic hypotheses for trunk neural crest cell migration. Here, we review recent data in support of three distinct mechanisms of trunk neural crest cell migration and develop and simulate a computational model based on chemotactic signaling. We show that by integrating the timing and spatial location of multiple chemotactic signals, trunk neural crest cells may be accurately positioned into two distinct targets that correspond to the sympathetic and dorsal root ganglia. In doing so, we honor the contributions of Wilhelm His to his identification of the neural crest and extend the observations of His and others to better understand a complex question in neural crest cell biology.     

Traveling waves in a simple population model involving growth and death

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.     

Existence of traveling waves for integral recursions with nonmonotone growth functions

A class of integral recursion models for the growth and spread of a synchronized single-species population is studied. It is well known that if there is no overcompensation in the fecundity function, the recursion has an asymptotic spreading speed c*, and that this speed can be characterized as the speed of the slowest non-constant traveling wave solution. A class of integral recursions with overcompensation which still have asymptotic spreading speeds can be found by using the ideas introduced by Thieme (J Reine Angew Math 306:94–121, 1979) for the study of space-time integral equation models for epidemics. The present work gives a large subclass of these models with overcompensation for which the spreading speed can still be characterized as the slowest speed of a non-constant traveling wave. To illustrate our results, we numerically simulate a series of traveling waves. The simulations indicate that, depending on the properties of the fecundity function, the tails of the waves may approach the carrying capacity monotonically, may approach the carrying capacity in an oscillatory manner, or may oscillate continually about the carrying capacity, with its values bounded above and below by computable positive numbers. B. Li’s research was partially supported by the National Science Foundation under Grant DMS-616445. M. A. Lewis research was supported by “The Canada Research Chairs program,” and a grant from the Natural Sciences and Engineering Research Council of Canada.     

A gradually slowing travelling band of chemotactic bacteria

A model for describing the motion of chemotactic bacteria in a capillary tube containing substrate is treated. Chemotactic substrate threshold effects are included in the chemotactic response coefficient. The ratio of the substrate threshold, s _T, to the substrate level far ahead of the travelling band, s , is used as a small parameter in developing an asymptotic solution of near travelling wave form.Also at: Membrane Filtration Technology, Kiryat Weizmann, Rehovot, Israel.Also at: Department of Mathematical Sciences, Rensselaer Polytechnic Institute, Troy, NY 12181, USA.     

Microwave-induced pressure waves in a model of muscle tissue

Microwave-induced mechanical stress waves were studied in simulated muscle tissue. Pulsed microwave energy at 5.655 GH_z induced pressure waves that were recorded with a hydrophone transducer. Each pulse produced a peak power density greater than 1.5 kW/cm². Microwave absorption measurements within the model showed energy deposition to be mostly confined to a region within 2 cm of the irradiated surface. The average specific absorption rate (SAR) at the surface of the sample was about 100 W/kg. The microwave-induced stress wave propagated at a velocity of 1,600 m/sec with peak pressures of approximately 300 pascals and was detectable after having traveled a total distance of 0.61 m on a path that included two reflections at model-container interfaces.     

A chemotactic method for the isolation of Actinoplanaceae

A simple chemotactic method for the isolation of Actinoplanaceae from soil is described. The method is based on a combination of the aerotactic behavior of the spores and the attraction to chloride ions. A simple isolation chamber is described. The method is simpler and less time-consuming than the current baiting techniques.     

Travelling waves in a tissue interaction model for skin pattern formation

Tissue interaction plays a major role in many morphogenetic processes, particularly those associated with skin organ primordia. We examine travelling wave solutions in a tissue interaction model for skin pattern formation which is firmly based on the known biology. From a phase space analysis we conjecture the existence of travelling waves with specific wave speeds. Subsequently, analytical approximations to the wave profiles are derived using perturbation methods. We then show numerically that such travelling wave solutions do exist and that they are in good agreement with our analytical results. Finally, the biological implications of our analysis are discussed.     

Cell balance equation for chemotactic bacteria with a biphasic tumbling frequency

Alts three-dimensional cell balance equation characterizing the chemotactic bacteria was analyzed under the presence of one-dimensional spatial chemoattractant gradients. Our work differs from that of others who have developed rather general models for chemotaxis in the use of a non-smooth anisotropic tumbling frequency function that responds biphasically to the combined temporal and spatial chemoattractant gradients. General three-dimensional expressions for the bacterial transport parameters were derived for chemotactic bacteria, followed by a perturbation analysis under the planar geometry. The bacterial random motility and chemotaxis were summarized by a motility tensor and a chemotactic velocity vector, respectively. The consequence of invoking the diffusion-approximation assumption and using intrinsic one-dimensional models with modified cellular swimming speeds was investigated by numerical simulations. Characterizing the bacterial random orientation after tumbles by a turn angle probability distribution function, we found that only the first-order angular moment of this turn angle probability distribution is important in influencing the bacterial long-term transport. Mathematics Subject Classification (2000):60G05, 60J60, 82A70     

Bias in the gradient-sensing response of chemotactic cells

We apply linear stability theory and perform perturbation studies to better characterize, and to generate new experimental predictions from, a model of chemotactic gradient sensing in eukaryotic cells. The model uses reaction-diffusion equations to describe 3(') phosphoinositide signaling and its regulation at the plasma membrane. It demonstrates a range of possible gradient-sensing mechanisms and captures such characteristic behaviors as strong polarization in response to static gradients, adaptation to differing mean levels of stimulus, and plasticity in response to changing gradients. An analysis of the stability of polarized steady-state solutions indicates that the model is most sensitive to off-axis perturbations. This biased sensitivity is also reflected in responses to localized external stimuli, and leads to a clear experimental prediction, namely, that a cell which is polarized in a background gradient will be most sensitive to transient point-source stimuli lying within a range of angles that are oblique with respect to the polarization axis. Stimuli at angles below this range will elicit responses whose directions overshoot the stimulus angle, while responses to stimuli applied at larger angles will undershoot the stimulus angle. We argue that such a bias is likely to be a general feature of gradient sensing in highly motile cells, particularly if they are optimized to respond to small gradients. Finally, an angular bias in gradient sensing might lead to preferred turn angles and zigzag movements of cells moving up chemotactic gradients, as has been noted under certain experimental conditions.     

Spatial spreading of West Nile Virus described by traveling waves

In this work, we propose a spatial model to analyze the West Nile Virus propagation across the USA, from east to west. West Nile Virus is an arthropod-borne flavivirus that appeared for the first time in New York City in the summer of 1999 and then spread prolifically among birds. Mammals, such as humans and horses, do not develop sufficiently high bloodstream titers to play a significant role in the transmission, which is the reason to consider the mosquito-bird cycle. The model aims to study this propagation based on a system of partial differential reaction-diffusion equations taking the mosquito and the avian populations into account. Diffusion and advection movements are allowed for both populations, being greater in the avian than in the mosquito population. The traveling wave solutions of the model are studied to determine the speed of disease dissemination. This wave speed is obtained as a function of the model's parameters, in order to assess the control strategies. The propagation of West Nile Virus from New York City to California state is established as a consequence of the diffusion and advection movements of birds. Mosquito movements do not play an important role in the disease dissemination, while bird advection becomes an important factor for lower mosquito biting rates.     

Dynamics of bacterial population waves as dependent on the type of the chemotactic response function

A study was made of the time of formation and the rate of propagation of bacterial population waves as dependent on substrate concentration. Experimentally and theoretically it is shown that the time of bacterial wave formation grows with increasing concentration of the nutritive component. At the same time, the wave propagation rate drops at relatively high concentrations of the nutritive component. In numerical experiments with a mathematical model, comparative analysis has been made of the dynamics of the behavior of bacterial waves at different chemotactic responses.