Dynamic regulation of the tryptophan operon: A modeling study and comparison with experimental data

A mathematical model for regulation of the tryptophan operon is presented. This model takes into account repression, feedback enzyme inhibition, and transcriptional attenuation. Special attention is given to model parameter estimation based on experimental data. The model's system of delay differential equations is numerically solved, and the results are compared with experimental data on the temporal evolution of enzyme activity in cultures of Escherichia coli after a nutritional shift (minimal + tryptophan medium to minimal medium). Good agreement is obtained between the numeric simulations and the experimental results for wild-type E. coli, as well as for two different mutant strains.

Blood flow and oxygen delivery to human brain during functional activity: Theoretical modeling and experimental data

Coupling of cerebral blood flow (CBF) and cerebral metabolic rate for oxygen (CMRO2) in physiologically activated brain states remains the subject of debates. Recently it was suggested that CBF is tightly coupled to oxidative metabolism in a nonlinear fashion. As part of this hypothesis, mathematical models of oxygen delivery to the brain have been described in which disproportionately large increases in CBF are necessary to sustain even small increases in CMRO2 during activation. We have explored the coupling of CBF and oxygen delivery by using two complementary methods. First, a more complex mathematical model was tested that differs from those recently described in that no assumptions were made regarding tissue oxygen level. Second, [15O] water CBF positron emission tomography (PET) studies in nine healthy subjects were conducted during states of visual activation and hypoxia to examine the relationship of CBF and oxygen delivery. In contrast to previous reports, our model showed adequate tissue levels of oxygen could be maintained without the need for increased CBF or oxygen delivery. Similarly, the PET studies demonstrated that the regional increase in CBF during visual activation was not affected by hypoxia. These findings strongly indicate that the increase in CBF associated with physiological activation is regulated by factors other than local requirements in oxygen.

Modeling the control of DNA replication in fission yeast

A central event in the eukaryotic cell cycle is the decision to commence DNA replication (S phase). Strict controls normally operate to prevent repeated rounds of DNA replication without intervening mitoses (“endoreplication”) or initiation of mitosis before DNA is fully replicated (“mitotic catastrophe”). Some of the genetic interactions involved in these controls have recently been identified in yeast. From this evidence we propose a molecular mechanism of “Start” control in Schizosaccharomyces pombe. Using established principles of biochemical kinetics, we compare the properties of this model in detail with the observed behavior of various mutant strains of fission yeast: wee1− (size control at Start), cdc13Δ and rum1OP (endoreplication), and wee1− rum1Δ (rapid division cycles of diminishing cell size). We discuss essential features of the mechanism that are responsible for characteristic properties of Start control in fission yeast, to expose our proposal to crucial experimental tests.

Modeling the dynamics of human hair cycles by a follicular automaton

The hair follicle cycle successively goes through the anagen, catagen, telogen, and latency phases, which correspond, respectively, to hair growth, arrest, shedding, and absence before a new anagen phase is initiated. Experimental observations collected over a period of 14 years in a group of 10 male volunteers, alopecic and nonalopecic, allowed us to determine the characteristics of scalp hair follicle cycles. On the basis of these observations, we propose a follicular automaton model to simulate the dynamics of human hair cycles. The automaton model is defined by a set of rules that govern the stochastic transitions of each follicle between the successive states anagen, telogen, and latency, and the subsequent return to anagen. The transitions occur independently for each follicle, after time intervals given stochastically by a distribution characterized by a mean and a variance. The follicular automaton model accounts both for the dynamical transitions observed in a single follicle and for the behavior of an ensemble of independently cycling follicles. Thus, the model successfully reproduces the evolution of the fractions of follicle populations in each of the three phases, which fluctuate around steady-state or slowly drifting values. We apply the follicular automaton model to the study of spatial patterns of follicular growth that result from a spatially heterogeneous distribution of parameters such as the mean duration of anagen phase. When considering that follicles die or miniaturize after going through a critical number of successive cycles, the model can reproduce the evolution to hair patterns similar to well known types of diffuse or androgenetic alopecia.

Singular value decomposition for genome-wide expression data processing and modeling

We describe the use of singular value decomposition in transforming genome-wide expression data from genes × arrays space to reduced diagonalized “eigengenes” × “eigenarrays” space, where the eigengenes (or eigenarrays) are unique orthonormal superpositions of the genes (or arrays). Normalizing the data by filtering out the eigengenes (and eigenarrays) that are inferred to represent noise or experimental artifacts enables meaningful comparison of the expression of different genes across different arrays in different experiments. Sorting the data according to the eigengenes and eigenarrays gives a global picture of the dynamics of gene expression, in which individual genes and arrays appear to be classified into groups of similar regulation and function, or similar cellular state and biological phenotype, respectively. After normalization and sorting, the significant eigengenes and eigenarrays can be associated with observed genome-wide effects of regulators, or with measured samples, in which these regulators are overactive or underactive, respectively.

Modeling back propagating action potential in weakly excitable dendrites of neocortical pyramidal cells.

Simultaneous recordings from the soma and apical dendrite of layer V neocortical pyramidal cells of young rats show that, for any location of current input, an evoked action potential (AP) always starts at the axon and then propagates actively, but decrementally, backward into the dendrites. This back-propagating AP is supported by a low density (-gNa = approximately 4 mS/cm2) of rapidly inactivating voltage-dependent Na+ channels in the soma and the apical dendrite. Investigation of detailed, biophysically constrained, models of reconstructed pyramidal cells shows the following. (i) The initiation of the AP first in the axon cannot be explained solely by morphological considerations; the axon must be more excitable than the soma and dendrites. (ii) The minimal Na+ channel density in the axon that fully accounts for the experimental results is about 20-times that of the soma. If -gNa in the axon hillock and initial segment is the same as in the soma [as recently suggested by Colbert and Johnston [Colbert, C. M. & Johnston, D. (1995) Soc. Neurosci. Abstr. 21, 684.2]], then -gNa in the more distal axonal regions is required to be about 40-times that of the soma. (iii) A backward propagating AP in weakly excitable dendrites can be modulated in a graded manner by background synaptic activity. The functional role of weakly excitable dendrites and a more excitable axon for forward synaptic integration and for backward, global, communication between the axon and the dendrites is discussed.