From cell death to embryo arrest: Mathematical models of human preimplantation embryo development

Human preimplantation embryos exhibit high levels of apoptotic cells and high rates of developmental arrest during the first week in vitro. The relation between the two is unclear and difficult to determine by conventional experimental approaches, partly because of limited numbers of embryos. We apply a mixture of experiment and mathematical modeling to show that observed levels of cell death can be reconciled with the high levels of embryo arrest seen in the human only if the developmental competence of embryos is already established at the zygote stage, and environmental factors merely modulate this. This suggests that research on improving in vitro fertilization success rates should move from its current concentration on optimizing culture media to focus more on the generation of a healthy zygote and on understanding the mechanisms that cause chromosomal and other abnormalities during early cleavage stages.

Failure of programmed cell death and differentiation as causes of tumors: some simple mathematical models.

Most models of tumorigenesis assume that the tumor grows by increased cell division. In these models, it is generally supposed that daughter cells behave as do their parents, and cell numbers have clear potential for exponential growth. We have constructed simple mathematical models of tumorigenesis through failure of programmed cell death (PCD) or differentiation. These models do not assume that descendant cells behave as their parents do. The models predict that exponential growth in cell numbers does sometimes occur, usually when stem cells fail to die or differentiate. At other times, exponential growth does not occur: instead, the number of cells in the population reaches a new, higher equilibrium. This behavior is predicted when fully differentiated cells fail to undergo PCD. When cells of intermediate differentiation fail to die or to differentiate further, the values of growth parameters determine whether growth is exponential or leads to a new equilibrium. The predictions of the model are sensitive to small differences in growth parameters. Failure of PCD and differentiation, leading to a new equilibrium number of cells, may explain many aspects of tumor behavior--for example, early premalignant lesions such as cervical intraepithelial neoplasia, the fact that some tumors very rarely become malignant, the observation of plateaux in the growth of some solid tumors, and, finally, long lag phases of growth until mutations arise that eventually result in exponential growth.

Topoisomerase II drives DNA transport by hydrolyzing one ATP

DNA topoisomerase II is a homodimeric molecular machine that couples ATP usage to the transport of one DNA segment through a transient break in another segment. In the presence of a nonhydrolyzable ATP analog, the enzyme is known to promote a single turnover of DNA transport. Current models for the enzyme’s mechanism based on this result have hydrolysis of two ATPs as the last step, used only to reset the enzyme for another round of reaction. Using rapid-quench techniques, topoisomerase II recently was shown to hydrolyze its two bound ATPs in a strictly sequential manner. This result is incongruous with the models based on the nonhydrolyzable ATP analog data. Here we present evidence that hydrolysis of one ATP by topoisomerase II precedes, and accelerates, DNA transport. These results indicate that important features of this enzyme’s mechanism previously have been overlooked because of the reliance on nonhydrolyzable analogs for studying a single reaction turnover. A model for the mechanism of topoisomerase II is presented to show how hydrolysis of one ATP could drive DNA transport.

Models of immune memory: On the role of cross-reactive stimulation, competition, and homeostasis in maintaining immune memory

There has been much debate on the contribution of processes such as the persistence of antigens, cross-reactive stimulation, homeostasis, competition between different lineages of lymphocytes, and the rate of cell turnover on the duration of immune memory and the maintenance of the immune repertoire. We use simple mathematical models to investigate the contributions of these various processes to the longevity of immune memory (defined as the rate of decline of the population of antigen-specific memory cells). The models we develop incorporate a large repertoire of immune cells, each lineage having distinct antigenic specificities, and describe the dynamics of the individual lineages and total population of cells. Our results suggest that, if homeostatic control regulates the total population of memory cells, then, for a wide range of parameters, immune memory will be long-lived in the absence of persistent antigen (T1/2 > 1 year). We also show that the longevity of memory in this situation will be insensitive to the relative rates of cross-reactive stimulation, the rate of turnover of immune cells, and the functional form of the term for the maintenance of homeostasis.

Variation in per capita interaction strength: Thresholds due to nonlinear dynamics and nonequilibrium conditions

I measured the strength of interaction between a marine herbivore and its growing resource over a realistic range of absolute and relative abundances. The herbivores (hermit crabs: Pagurus spp.) have slow and/or weak functional and numerical responses to epiphytic diatoms (Isthmia nervosa), which show logistic growth in the absence of consumers. By isolating this interaction in containers in the field, I mimicked many of the physical and biological variables characteristic of the intertidal while controlling the densities of focal species. The per capita effects of consumers on the population dynamics of their resource (i.e., interaction strength) were defined by using the relationship between hermit crab density and proportional change in the resource. When this relationship is fit by a Weibull function, a single parameter distinguishes constant interaction strength from one that varies as a function of density. Constant interaction strength causes the proportion of diatoms to fall linearly or proportionally as hermit crab density increases whereas per capita effects that increase with density cause an accelerating decline. Although many mathematical models of species interactions assume linear dynamics and invariant parameters, at least near equilibrium, the per capita effects of hermit crabs on diatoms varied substantially, apparently crossing a threshold from weak to strong when consumption exceeded resource production. This threshold separates a domain of coexistence from one of local extinction of the resource. Such thresholds may help explain trophic cascades, resource compensation, and context-dependent interaction strengths, while indicating a way to predict trophic effects, despite nonlinearities, as a function of vital rates.

Different population dynamics of human T cell lymphotropic virus type II in intravenous drug users compared with endemically infected tribes

The phylogeny of human T cell lymphotropic virus type II (HTLV-II) was investigated by using strains isolated from Amerindian and Pygmy tribes, in which the virus is maintained primarily through mother-to-child transmission via breast-feeding, and strains from intravenous drug users (IDUs), in which spread is mainly blood-borne via needle sharing. Molecular clock analysis showed that HTLV-II has two different evolutionary rates with the molecular clock for the virus in IDUs ticking 150–350 times faster than the one in endemically infected tribes: 2.7 × 10−4 compared with 1.71/7.31 × 10−7 nucleotide substitutions per site per year in the long terminal repeat region. This dramatic acceleration of the evolutionary rate seems to be related with the mode of transmission. Mathematical models showed the correlation of these two molecular clocks with an endemic spread of HTLV-II in infected tribes compared with the epidemic spread in IDUs. We also noted a sharp increase in the population size of the virus among IDUs during the last decades probably caused by the worldwide increase in intravenous drug use.

Specific therapy regimes could lead to long-term immunological control of HIV

We use mathematical models to study the relationship between HIV and the immune system during the natural course of infection and in the context of different antiviral treatment regimes. The models suggest that an efficient cytotoxic T lymphocyte (CTL) memory response is required to control the virus. We define CTL memory as long-term persistence of CTL precursors in the absence of antigen. Infection and depletion of CD4+ T helper cells interfere with CTL memory generation, resulting in persistent viral replication and disease progression. We find that antiviral drug therapy during primary infection can enable the development of CTL memory. In chronically infected patients, specific treatment schedules, either including deliberate drug holidays or antigenic boosts of the immune system, can lead to a re-establishment of CTL memory. Whether such treatment regimes would lead to long-term immunologic control deserves investigation under carefully controlled conditions.

Cytotoxic T lymphocytes and viral turnover in HIV type 1 infection

To understand the role of the immune system in limiting HIV type 1 replication, it is critical to know to what extent the rapid turnover of productively infected cells is caused by viral cytopathicity or by immune-mediated lysis. We show that uncultured peripheral blood mononuclear cells of many patients contain cytotoxic T lymphocytes (CTL) that lyse target cells—at plausible peripheral blood mononuclear cell-to-target ratios—with half-lives of less than 1 day. In 23 patients with CD4 counts ranging from 10 to 900 per μl, the average rate of CTL-mediated lysis corresponds to a target cell half-life of 0.7 day. We develop mathematical models to calculate the turnover rate of infected cells subjected to immune-mediated lysis and viral cytopathicity and to estimate the fraction of cells that are killed by CTL as opposed to virus. The models provide new interpretations of drug treatment dynamics and explain why the observed rate of virus decline is roughly constant for different patients. We conclude that in HIV type 1 infection, CTL-mediated lysis can reduce virus load by limiting virus production, with small effects on the half-life of infected cells.

Frontiers of finance: Evolution and efficient markets

In this review article, we explore several recent advances in the quantitative modeling of financial markets. We begin with the Efficient Markets Hypothesis and describe how this controversial idea has stimulated a number of new directions of research, some focusing on more elaborate mathematical models that are capable of rationalizing the empirical facts, others taking a completely different tack in rejecting rationality altogether. One of the most promising directions is to view financial markets from a biological perspective and, specifically, within an evolutionary framework in which markets, instruments, institutions, and investors interact and evolve dynamically according to the “law” of economic selection. Under this view, financial agents compete and adapt, but they do not necessarily do so in an optimal fashion. Evolutionary and ecological models of financial markets is truly a new frontier whose exploration has just begun.

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.

Antigenic variation and the within-host dynamics of parasites.

Many parasites exhibit antigenic variation within their hosts. We use mathematical models to investigate the dynamical interaction between an antigenically varying parasite and the host's immune system. The models incorporate antigenic variation in the parasite population and the generation of immune responses directed against (i) antigens specific to individual parasite variants and (ii) antigens common to all the parasite variants. Analysis of the models allows us to evaluate the relative importance of variant-specific and cross-reactive immune responses in controlling the parasite. Early in the course of infection within the host, when parasite diversity is below a defined threshold value (the value is determined by the biological properties of the parasite and of the host's immune response), the variant-specific immune responses are predominant. Later, when the parasite diversity is high, the cross-reactive immune response is largely responsible for controlling the parasitemia. It is argued that increasing antigenic diversity leads to a switch from variant-specific to cross-reactive immune responses. These simple models mimic various features of observed infections recorded in the experimental literature, including an initial peak in parasitemia, a long and variable duration of infection with fluctuating parasitemia that ends with either the clearance of the parasite or persistent infection.

Complexity, contingency, and criticality.

Complexity originates from the tendency of large dynamical systems to organize themselves into a critical state, with avalanches or "punctuations" of all sizes. In the critical state, events which would otherwise be uncoupled become correlated. The apparent, historical contingency in many sciences, including geology, biology, and economics, finds a natural interpretation as a self-organized critical phenomenon. These ideas are discussed in the context of simple mathematical models of sandpiles and biological evolution. Insights are gained not only from numerical simulations but also from rigorous mathematical analysis.