Pattern formation in an immobilized bienzyme system. A morphogenetic model

A reaction-diffusion model of ROS-induced ROS release in a mitochondrial network

Loss of mitochondrial function is a fundamental determinant of cell injury and death. In heart cells under metabolic stress, we have previously described how the abrupt collapse or oscillation of the mitochondrial energy state is synchronized across the mitochondrial network by local interactions dependent upon reactive oxygen species (ROS). Here, we develop a mathematical model of ROS-induced ROS release (RIRR) based on reaction-diffusion (RD-RIRR) in one- and two-dimensional mitochondrial networks. The nodes of the RD-RIRR network are comprised of models of individual mitochondria that include a mechanism of ROS-dependent oscillation based on the interplay between ROS production, transport, and scavenging; and incorporating the tricarboxylic acid (TCA) cycle, oxidative phosphorylation, and Ca(2+) handling. Local mitochondrial interaction is mediated by superoxide (O2.-) diffusion and the O2.(-)-dependent activation of an inner membrane anion channel (IMAC). In a 2D network composed of 500 mitochondria, model simulations reveal DeltaPsi(m) depolarization waves similar to those observed when isolated guinea pig cardiomyocytes are subjected to a localized laser-flash or antioxidant depletion. The sensitivity of the propagation rate of the depolarization wave to O(2.-) diffusion, production, and scavenging in the reaction-diffusion model is similar to that observed experimentally. In addition, we present novel experimental evidence, obtained in permeabilized cardiomyocytes, confirming that DeltaPsi(m) depolarization is mediated specifically by O2.-). The present work demonstrates that the observed emergent macroscopic properties of the mitochondrial network can be reproduced in a reaction-diffusion model of RIRR. Moreover, the findings have uncovered a novel aspect of the synchronization mechanism, which is that clusters of mitochondria that are oscillating can entrain mitochondria that would otherwise display stable dynamics. The work identifies the fundamental mechanisms leading from the failure of individual organelles to the whole cell, thus it has important implications for understanding cell death during the progression of heart disease.

Biological processes in organised media

Embedding a simple Michaelis-Menten enzyme in a gel slice may allow the catalysis of not only scalar processes but also vectorial ones, including uphill transport of a substrate between two compartments, and may make it seem as if two enzymes or transporters are present or as if an allosterically controlled enzyme/transporter is operating. The values of kinetic parameters of an enzyme in a partially hydrophobic environment are usually different from those actually measured in a homogeneous aqueous solution. This implies that fitting kinetic data (expressed in reciprocal co-ordinates) from in vivo studies of enzymes or transporters to two straight lines or a sigmoidal curve does not prove the existence of two different membrane mechanisms or allosteric control. In the artificial transport systems described here, a functional asymmetry was sufficient to induce uphill transport, therefore, although the active transport systems characterised so far correspond to proteins asymmetrically anchored in a membrane, the past or present existence of structurally symmetrical systems of transport in vivo cannot be excluded. The fact that oscillations can be induced in studies of the maintenance of the electrical potential of frog skin by addition of lithium allowed evaluation of several parameters fundamental to the functioning of the system in vivo (e.g., relative volumes of internal compartments, characteristic times of ionic exchanges between compartments). Hence, under conditions that approach real biological complexity, increasing the complexity of the behaviour of the system may provide information that cannot be obtained by a conventional, reductionist approach.

Chemical or biological activity in open chaotic flows

We investigate the evolution of particle ensembles in open chaotic hydrodynamical flows. Active processes of the type A+B-->2B and A+B-->2C are considered in the limit of weak diffusion. As an illustrative advection dynamics we consider a model of the von Kármán vortex street, a time-periodic two-dimensional flow of a viscous fluid around a cylinder. We show that a fractal unstable manifold acts as a catalyst for the process, and the products cover fattened-up copies of this manifold. This may account for the observed filamental intensification of activity in environmental flows. The reaction equations valid in the wake are derived either in the form of dissipative maps or differential equations depending on the regime under consideration. They contain terms that are not present in the traditional reaction equations of the same active process: the decay of the products is slower while the productivity is much faster than in homogeneous flows. Both effects appear as a consequence of underlying fractal structures. In the long time limit, the system locks itself in a dynamic equilibrium state synchronized to the flow for both types of reactions. For particles of finite size an emptying transition might also occur leading to no products left in the wake.

An allometric interpretation of the spatio-temporal organization of molecular and cellular processes

Different levels of organization distinguished by characteristics spatial dimensions, Ec, and relaxation times, Tr, of biological processes ranging from electron transport in energy transduction to growth of microbial and plant cells, are shown to be related through a relation that may be interpreted as allometric and characterized by two different slopes. Processes, at levels of organization occurring in spatial dimensions of micrometers and relaxing in the order of minutes, delimit a 'transition point' between the two curves, that we interpret as a limit for the emergence of macroscopic coherence. The characteristic spatial dimension, Ec, and the relaxation time, Tr, contain dynamical information about the processes occurring at a given level of organization. When a steady state of a biological process at a certain level of organization becomes unstable, the system undergoes a transition to another level of organization. To exemplify the appearance of macroscopic order at levels of organization further from the 'transition point' we present in this report various experimental systems involving many levels of organization allometrically related that exhibit different kinds of self-organized behavior, i.e. bi-stability, oscillations, changes in (a)symmetry.