Controllability of selection-mutation systems

The well-known Fisher type selection-mutation model is studied from the point of view of mathematical systems theory. Mutation rates are considered as control functions. Based on a general sufficient condition for local controllability of non-linear systems with invariant manifold, a method is proposed to guarantee the controllability of the considered population into a polymorphic equilibrium.     

Controllability of non-linear biochemical systems

Mathematical methods of biochemical pathway analysis are rapidly maturing to a point where it is possible to provide objective rationale for the natural design of metabolic systems and where it is becoming feasible to manipulate these systems based on model predictions, for instance, with the goal of optimizing the yield of a desired microbial product. So far, theory-based metabolic optimization techniques have mostly been applied to steady-state conditions or the minimization of transition time, using either linear stoichiometric models or fully kinetic models within biochemical systems theory (BST). This article addresses the related problem of controllability, where the task is to steer a non-linear biochemical system, within a given time period, from an initial state to some target state, which may or may not be a steady state. For this purpose, BST models in S-system form are transformed into affine non-linear control systems, which are subjected to an exact feedback linearization that permits controllability through independent variables. The method is exemplified with a small glycolytic-glycogenolytic pathway that had been analyzed previously by several other authors in different contexts.     

Three-dimensional waves in excitable reaction-diffusion systems

The eikonal equation [5] for excitable media is generalised to three dimensional systems. The main result of the investigation is the demonstration of the existence of toroidal and twisted toroidal scroll waves in the limit of large values of the major radius of the torus. The existence of a helical wave near the z-axis follows from the eidonal equation but its connection with the twisted toroidal scroll remains to be demonstrated. The eikonal equation also predicts a non-uniform rate of rotation of the cross-sectional spiral wave near the toroidal axis. The notion of geometrical stability is introduced for the case of an expanding sphere; in particular it is shown that a discussion of stability of solutions of the eikonal equation must take into account the possible shift in the origin of the coordinate systems with respect to which patterns are defined.On leave of absence from: The Department of Mathematics, Glasgow College of Technology, Cowcaddens Road, Glasgow G4 0BA, UK     

Negative refractoriness in excitable systems with cross-diffusion

The results of numerical experiments with mathematical models of excitable systems with cross-diffusion are presented. It was shown that the refractoriness in such systems may be negative. The effects of negative refractoriness on the propagation and interaction of waves are demonstrated.     

Subharmonic resonance and chaos in forced excitable systems

 Forced excitable systems arise in a number of biological and physiological applications and have been studied analytically and computationally by numerous authors. Existence and stability of harmonic and subharmonic solutions of a forced piecewise-linear Fitzhugh-Nagumo-like system were studied in Othmer ad Watanabe (1994) and in Xie et al. (1996). The results of those papers were for small and moderate amplitude forcing. In this paper we study the existence of subharmonic solutions of this system under large-amplitude forcing. As in the case of intermediate-amplitude forcing, bistability between 1 : 1 and 2 : 1 solutions is possible for some parameters. In the case of large-amplitude forcing, bistability between 2 : 2 and 2 : 1 solutions, which does not occur in the case of intermediate-amplitude forcing, is also possible for some parameters. We identify several new canonical return maps for a singular system, and we show that chaotic dynamics can occur in some regions of parameter space. We also prove that there is a direct transition from 2 : 2 phase-locking to chaos after the first period-doubling bifurcation, rather than via the infinite sequence of period doublings seen in a smooth quadratic interval map. Coexistence of chaotic dynamics and stable phase-locking can also occur. Received: 6 July 1998 / Revised version: 2 October 1998     

Mathematical models for excitable systems in biology and medicine.

The excitable systems play a very important role in Biology and Medicine. Phenomena such as the transmission of impulses between neurons, the cardiac arrhythmia, the aggregation of amoebas, the appearance of organized structures in the cortex of egg cells, all derive from the activity of excitable media. In the first part of this work a general definition of excitable system is given; we then analyze some cases of excitability, distinguishing between electrical and chemical excitability and comparing experimental observations with simulations carried out by appropriate mathematical models. Such models are almost always formulated by partial differential equations of "reaction-diffusion" type and they have the characteristic to describe propagations of electrical waves (neurons, pacemaker cardiac cells, pancreatic b-cells) or chemical and mechanical waves (propagation of Ca++ waves or mechanical waves in the endoplasmic reticulum). The aim is to put in evidence that the biological systems can show not only excitability of electrical type, but also excitability of chemical nature, which can be observed in the first steps of development of egg cells or, for example, in the formation of pigments in vertebrate skin or in clam shells.     

Hopf bifurcation and bursting synchronization in an excitable systems with chemical delayed coupling

In this paper we consider the Hopf bifurcation and synchronization in the two coupled Hindmarsh–Rose excitable systems with chemical coupling and time-delay. We surveyed the conditions for Hopf bifurcations by means of dynamical bifurcation analysis and numerical simulation. The results show that the coupled excitable systems with no delay have supercritical Hopf bifurcation, while the delayed system undergoes Hopf bifurcations at critical time delays when coupling strength lies in a particular region. We also investigated the effect of the delay on the transition of bursting synchronization in the coupled system. The results are helpful for us to better understand the dynamical properties of excitable systems and the biological mechanism of information encoding and cognitive activity.     

The excitable fluid mosaic

The Fluid Mosaic Model by Singer & Nicolson proposes that biological membranes consist of a fluid lipid layer into which integral proteins are embedded. The lipid membrane acts as a two-dimensional liquid in which the proteins can diffuse and interact. Until today, this view seems very reasonable and is the predominant picture in the literature. However, there exist broad melting transitions in biomembranes some 10–20 degrees below physiological temperatures that reach up to body temperature. Since they are found below body temperature, Singer & Nicolson did not pay any further attention to the melting process. But this is a valid view only as long as nothing happens. The transition temperature can be influenced by membrane tension, pH, ionic strength and other variables. Therefore, it is not generally correct that the physiological temperature is above this transition. The control over the membrane state by changing the intensive variables renders the membrane as a whole excitable. One expects phase behavior and domain formation that leads to protein sorting and changes in membrane function. Thus, the lipids become an active ingredient of the biological membrane. The melting transition affects the elastic constants of the membrane. This allows for the generation of propagating pulses in nerves and the formation of ion-channel-like pores in the lipid membranes. Here we show that on top of the fluid mosaic concept there exists a wealth of excitable phenomena that go beyond the original picture of Singer & Nicolson.¹     

Electrical properties of an excitable epithelium

The exumbrellar epithelium of the hydromedusa, Euphysa japonica, is composed of a single layer of broad (70 micrometers), thin (1--2 micrometers) cells which are joined by gap junctions and simple appositions. Although the epithelium lacks nerves, it is excitable; electrically stimulating the epithelium initiates a propagated action potential. The average resting potential of the epithelial cells is -46 mV. The action potential, recorded with an intracellular electrode, is an all-or-nothing, positive, overshooting spike. The epithelial cells are electrically coupled. The passive electrical properties of the epithelium were determined from the decrement in membrane hyperpolarization with distance from an intracellular, positive current source. The two-dimensional space constant of the epithelium is 1.3 mm, the internal longitudinal resistivity of the cytoplasm and intercellular junctions is 196 omega cm, and the resistivity of both apical and basal cell membranes is greater than 23 k omega cm2. Although the membrane resistivity is high, the transverse resistivity of the epithelium is quite low (7.5 omega cm2), indicating that the epithelium is leaky with a large, transverse, paracellular shunt.     

Bottom-up excitable models of phytoplankton blooms

A simple nutrient-phytoplankton model is used to explore the dynamics of phytoplankton blooms. The model exhibits excitable behaviour in the sense that a large scale outbreak can only be triggered when a critical nutrient threshold is exceeded. The model takes into account several features often neglected but whose combined effect proves very important: (i) rapid nutrient recycling associated with the microbial loop and patch formation; (ii) self-shading; and (iii) a bottom-up approach, whereby nutrient levels are responsible for both the triggering and the demise of the bloom. Although the literature is replete with studies of ‘top-down’ models in which zooplankton grazing control the triggering and demise of the bloom, bottom-up models are nevertheless appropriate in many circumstances. We provide a full mathematical investigation of the effects of these three different features in an excitable system framework.     

Ionic processes in excitable membranes

The dynamical theory of ionized media is applied to the semi-electrolyte component of an excitable cell membrane, and the adjacent electrolytes. The equations of conservation of charge and momentum for the ions, and Poisson's equation for the electrostatic potential, are applied first to investigate the steady states of the membrane, and then transient effects in the membrane. A dispersion equation is derived, and the characteristic modes of relaxation within the membrane are determined. Among these are oscillating modes whose frequencies and amplitudes are of the correct order of magnitude to explain the observed excitation phenomena.A pair of coupled non-linear equations in the ionic potentials, with action potential solutions, is derived from the time-dependent electrodiffusion equations, and calculations are presented which model the behaviour of the excitable membrane during the voltage clamp. It is not necessary to postulate large changes in the ionic permeabilities in the course of the action potential and the voltage clamp to account for the large transient membrane conductances. It is suggested that the sodium hypothesis be replaced by one which attributes the action potential to non-linear plasma oscillations.     

Gap junctions in excitable cells

Gap junction channels are an integral part of the conduction or propagation of an action potential from cell to cell. Gap junctions have rather unique gating and permeability properties which permit the movement of molecules from cell to cell. These molecules may not be directly linked to action potentials but can alter nonjunctional processes within cells, which in turn can affect conduction velocity. The data described in this review reveal that, for the majority of excitable cells, there are two limiting factors, with respect to gap junctions, that affect the conduction/propagation of action potentials. These are (1) the total number of channels and (2) the selective permeability of the channels. Interestingly, voltage dependence and the time course of voltage inactivation (kinetics) are not rate limiting steps under normal physiological conditions for any of the connexins studied so far. Only specialized rectifying electrical synapses utilize strong voltage dependence and rapid kinetics to permit or deny the continued propagation of an action potential.     

Exact Controllability of Age-Structured Predator-Prey Populations

1IntroductionT5e劝dysisd8贝ami。dpTdatO。preynit。ti。is4iquitousinthebiologicalworldandlsd。id既曲ienit。estan6ipeceto。IOgy．Thetheoretidstudydpredation晌hadahghistoryho。lugwiththe。nalworkofInne．aandVolt。,andcontin。tobeofwid聊readinterest咖y．仇iiithelateof1980’s,most…ator－pteystUdieshavebudonInteractingspedeswithoutsg。sin．－Cture（forexmple,sce［1」,［2」,［3」and［4］）．HOwever,asthelin阿mceof明卜s刎cmreInpoPulaMnhasbecomemorewidelyrecognized,thereisanIn－er。lugnum…