New Schemes of Dynamic Preservation of Diversity: Remarks on Stability and Topology.

We address the biological dynamics problem of the persistence of several species in conditions of non-existence of an equilibrium, including an example of stabilization by predation and the very controversial "competitive exclusion" (which depends on the precise definition of persistence). We give normal forms for various examples of such (essentially dynamical) persistence and comments on the involved topology, which implies the presence of exceptional heteroclinic connections binding equilibria on the boundary.

A Mathematical Model for Alternation of Polygamy and Parthenogenesis: Stability Versus Efficiency and Analogy with Parasitism.

The present work is a contribution to the understanding of the sempiternal problem of the "burden of factor two" implied by sexual reproduction versus asexual one, as males are energy consumers not contributing to the production of offspring. We construct a deterministic mathematical model in population dynamics where a species enjoys both sexual and parthenogenetic capabilities of reproduction and lives on a limited resource. We then show how polygamy implies instability of a parthenogenetic population with a small number of sexually born males. This instability implies evolution of the system towards an attractor involving both (sexual and asexual) populations (which does not imply optimality of the population). We also exhibit the analogy with a parasite/host system.

On a minimal model for hemodynamics and metabolism of lactate: application to low grade glioma and therapeutic strategies.

WHO II low grade glioma evolves inevitably to anaplastic transformation. Magnetic resonance imaging is a good non-invasive way to watch it, by hemodynamic and metabolic modifications, thanks to multinuclear spectroscopy (1)H/(31)P. In this work we study a multi-scale minimal model of hemodynamics and metabolism applied to the study of gliomas. This mathematical analysis leads us to a fast-slow system. The control of the position of the stationary point brings to the concept of domain of viability. Starting from this system, the equations bring to light the parameters that push glioma cells out of their domain of viability. Four fundamental factors are highlighted. The first two are cerebral blood flow and the rate of lactate transport through monocarboxylate transporters, which must be reduced in order to push glioma out of its domain of viability. Another factor is the intra arterial lactate, which must be increased. The last factor is pH, indeed a decrease of intra cellular pH could interfere with glioma growth. These reflections suggest that these four parameters could lead to new therapeutic strategies for the management of low grade gliomas.

Mathematical modeling of energy metabolism and hemodynamics of WHO grade II gliomas using in vivo MR data.

Therapeutic management of low-grade gliomas (LGG) is a challenge because they have undergone anaplastic transformation with variable delay. Today, only progressive volume growth on successive MRI allows an in vivo monitoring of this evolution. On the other hand, multinuclear spectroscopy and perfusion available during MRI may also provide assessment of metabolic changes underlying morphological modifications. To overcome this drawback, we developed a mathematical model of the metabolism and the hemodynamic of gliomas, based on a physiological model previously published, and including the MR parameters. This allows us to suggest that some specific profiles of metabolic and hemodynamic changes would be good indicators of potential anaplastic transformation.

Propagation of bursting oscillations.

We investigate a system of partial differential equations of reaction-diffusion type which displays propagation of bursting oscillations. This system represents the time evolution of an assembly of cells constituted by a small nucleus of bursting cells near the origin immersed in the middle of excitable cells. We show that this system displays a global attractor in an appropriated functional space. Numerical simulations show the existence in this attractor of recurrent solutions which are waves propagating from the central source. The propagation seems possible if the excitability of the neighbouring cells is above some threshold.

SAPHIR: a physiome core model of body fluid homeostasis and blood pressure regulation.

We present the current state of the development of the SAPHIR project (a Systems Approach for PHysiological Integration of Renal, cardiac and respiratory function). The aim is to provide an open-source multi-resolution modelling environment that will permit, at a practical level, a plug-and-play construction of integrated systems models using lumped-parameter components at the organ/tissue level while also allowing focus on cellular- or molecular-level detailed sub-models embedded in the larger core model. Thus, an in silico exploration of gene-to-organ-to-organism scenarios will be possible, while keeping computation time manageable. As a first prototype implementation in this environment, we describe a core model of human physiology targeting the short- and long-term regulation of blood pressure, body fluids and homeostasis of the major solutes. In tandem with the development of the core models, the project involves database implementation and ontology development.

SAPHIR - a multi-scale, multi-resolution modeling environment targeting blood pressure regulation and fluid homeostasis.

We present progress on a comprehensive, modular, interactive modeling environment centered on overall regulation of blood pressure and body fluid homeostasis. We call the project SAPHIR, for "a Systems Approach for PHysiological Integration of Renal, cardiac, and respiratory functions". The project uses state-of-the-art multi-scale simulation methods. The basic core model will give succinct input-output (reduced-dimension) descriptions of all relevant organ systems and regulatory processes, and it will be modular, multi-resolution, and extensible, in the sense that detailed submodules of any process(es) can be "plugged-in" to the basic model in order to explore, eg. system-level implications of local perturbations. The goal is to keep the basic core model compact enough to insure fast execution time (in view of eventual use in the clinic) and yet to allow elaborate detailed modules of target tissues or organs in order to focus on the problem area while maintaining the system-level regulatory compensations.

Approximation for limit cycles and their isochrons.

Local analysis of trajectories of dynamical systems near an attractive periodic orbit displays the notion of asymptotic phase and isochrons. These notions are quite useful in applications to biosciences. In this note, we give an expression for the first approximation of equations of isochrons in the setting of perturbations of polynomial Hamiltonian systems. This method can be generalized to perturbations of systems that have a polynomial integral factor (like the Lotka-Volterra equation).

A differential geometry approach for biomedical image processing.

We show in this paper how simple considerations about bio-arrays images lead to a peak segmentation allowing the genes activity analysis. Bio-arrays images have a particular structure and the aim of the paper is to present a mathematical method allowing their automatic processing. The differential geometry approach used here can be also employed for other types of images presenting grey level peaks corresponding to a functional activity or to a chemical concentration. The mathematical method is based on elementary techniques of differential geometry and dynamical systems theory and provides a simple efficient algorithm when the peaks to segment are isolated.