Stochastic Chlamydia Dynamics and Optimal Spread.

Chlamydia trachomatis is an important bacterial pathogen that has an unusual developmental switch from a dividing form (reticulate body or RB) to an infectious form (elementary body or EB). RBs replicate by binary fission within an infected host cell, but there is a delay before RBs convert into EBs for spread to a new host cell. We developed stochastic optimal control models of the Chlamydia developmental cycle to examine factors that control the number of EBs produced. These factors included the probability and timing of conversion, and the duration of the developmental cycle before the host cell lyses. Our mathematical analysis shows that the observed delay in RB-to-EB conversion is important for maximizing EB production by the end of the intracellular infection.

Regulatory feedback on receptor and non-receptor synthesis for robust signaling.

Elaborate regulatory feedback processes are thought to make biological development robust, that is, resistant to changes induced by genetic or environmental perturbations. How this might be done is still not completely understood. Previous numerical simulations on reaction-diffusion models of Dpp gradients in Drosophila wing imaginal disc have showed that feedback (of the Hill function type) on (signaling) receptors and/or non-(signaling) receptors are of limited effectiveness in promoting robustness. Spatial nonuniformity of the feedback processes has also been shown theoretically to lead to serious shape distortion and a principal cause for ineffectiveness. Through mathematical modeling and analysis, the present article shows that spatially uniform nonlocal feedback mechanisms typically modify gradient shape through a shape parameter (that does not change with location). This in turn enables us to uncover new multi-feedback instrument for effective promotion of robust signaling gradients.

Fastest time to cancer by loss of tumor suppressor genes.

Genetic instability promotes cancer progression (by increasing the probability of cancerous mutations) as well as hinders it (by imposing a higher cell death rate for cells susceptible to cancerous mutation). With the loss of tumor suppressor gene function known to be responsible for a high percentage of breast and colorectal cancer (and a good fraction of lung cancer and other types as well), it is important to understand how genetic instability can be orchestrated toward carcinogenesis. In this context, this paper gives a complete characterization of the optimal (time-varying) cell mutation rate for the fastest time to a target cancerous cell population through the loss of both copies of a tumor suppressor gene. Similar to the (one-step) oncogene activation model previously analyzed, the optimal mutation rate of the present two-step model changes qualitatively with the convexity of the (mutation rate-dependent) cell death rate. However, the structure of the Hamiltonian for the new model differs significantly and intrinsically from that of the one-step model, and a completely new approach is needed for the solution of the present two-step problem. Considerable insight into the biology of optimal switching (between corner controls) is extracted from numerical results for cases with nonconvex death rates.

ROBUSTNESS OF MORPHOGEN GRADIENTS WITH "BUCKET BRIGADE" TRANSPORT THROUGH MEMBRANE-ASSOCIATED NON-RECEPTORS.

Robust multiple-fate morphogen gradients are essential for embryo development. Here, we analyze mathematically a model of morphogen gradient (such as Dpp in Drosophila wing imaginal disc) formation in the presence of non-receptors with both diffusion of free morphogens and the movement of morphogens bound to non-receptors. Under the assumption of rapid degradation of unbound morphogen, we introduce a method of functional boundary value problem and prove the existence, uniqueness and linear stability of a biologically acceptable steady-state solution. Next, we investigate the robustness of this steady-state solution with respect to significant changes in the morphogen synthesis rate. We prove that the model is able to produce robust biological morphogen gradients when production and degradation rates of morphogens are large enough and non-receptors are abundant. Our results provide mathematical and biological insight to a mechanism of achieving stable robust long distance morphogen gradients. Key elements of this mechanism are rapid turnover of morphogen to non-receptors of neighoring cells resulting in significant degradation and transport of non-receptor-morphogen complexes, the latter moving downstream through a "bucket brigade" process.

Cell lineages and the logic of proliferative control.

It is widely accepted that the growth and regeneration of tissues and organs is tightly controlled. Although experimental studies are beginning to reveal molecular mechanisms underlying such control, there is still very little known about the control strategies themselves. Here, we consider how secreted negative feedback factors ("chalones") may be used to control the output of multistage cell lineages, as exemplified by the actions of GDF11 and activin in a self-renewing neural tissue, the mammalian olfactory epithelium (OE). We begin by specifying performance objectives-what, precisely, is being controlled, and to what degree-and go on to calculate how well different types of feedback configurations, feedback sensitivities, and tissue architectures achieve control. Ultimately, we show that many features of the OE-the number of feedback loops, the cellular processes targeted by feedback, even the location of progenitor cells within the tissue-fit with expectations for the best possible control. In so doing, we also show that certain distinctions that are commonly drawn among cells and molecules-such as whether a cell is a stem cell or transit-amplifying cell, or whether a molecule is a growth inhibitor or stimulator-may be the consequences of control, and not a reflection of intrinsic differences in cellular or molecular character.

ROBUSTNESS OF SIGNALING GRADIENT IN DROSOPHILA WING IMAGINAL DISC.

Quasi-stable gradients of signaling protein molecules (known as morphogens or ligands) bound to cell receptors are known to be responsible for differential cell signaling and gene expressions. From these follow different stable cell fates and visually patterned tissues in biological development. Recent studies have shown that the relevant basic biological processes yield gradients that are sensitive to small changes in system characteristics (such as expression level of morphogens or receptors) or environmental conditions (such as temperature changes). Additional biological activities must play an important role in the high level of robustness observed in embryonic patterning for example. It is natural to attribute observed robustness to various type of feedback control mechanisms. However, our own simulation studies have shown that feedback control is neither necessary nor sufficient for robustness of the morphogen decapentaplegic (Dpp) gradient in wing imaginal disc of Drosophilas. Furthermore, robustness can be achieved by substantial binding of the signaling morphogen Dpp with nonsignaling cell surface bound molecules (such as heparan sulfate proteoglygans) and degrading the resulting complexes at a sufficiently rapid rate. The present work provides a theoretical basis for the results of our numerical simulation studies.

Spatial dynamics of multistage cell lineages in tissue stratification.

In developing and self-renewing tissues, terminally differentiated (TD) cell types are typically specified through the actions of multistage cell lineages. Such lineages commonly include a stem cell and multiple progenitor (transit-amplifying) cell stages, which ultimately give rise to TD cells. As the tissue reaches a tightly controlled steady-state size, cells at different lineage stages assume distinct spatial locations within the tissue. Although tissue stratification appears to be genetically specified, the underlying mechanisms that direct tissue lamination are not yet completely understood. Herein, we use modeling and simulations to explore several potential mechanisms that can be utilized to create stratification during developmental or regenerative growth in general systems and in the model system, the olfactory epithelium of mouse. Our results show that tissue stratification can be generated and maintained through controlling spatial distribution of diffusive signaling molecules that regulate the proliferation of each cell type within the lineage. The ability of feedback molecules to stratify a tissue is dependent on a low TD death rate: high death rates decrease tissue lamination. Regulation of the cell cycle lengths of stem cells by feedback signals can lead to transient accumulation of stem cells near the base and apex of tissue.

Formation of the BMP activity gradient in the Drosophila embryo.

The dorsoventral axis of the Drosophila embryo is patterned by a gradient of bone morphogenetic protein (BMP) ligands. In a process requiring at least three additional extracellular proteins, a broad domain of weak signaling forms and then abruptly sharpens into a narrow dorsal midline peak. Using experimental and computational approaches, we investigate how the interactions of a multiprotein network create the unusual shape and dynamics of formation of this gradient. Starting from observations suggesting that receptor-mediated BMP degradation is an important driving force in gradient dynamics, we develop a general model that is capable of capturing both subtle aspects of gradient behavior and a level of robustness that agrees with in vivo results.

Do morphogen gradients arise by diffusion?

Many patterns of cell and tissue organization are specified during development by gradients of morphogens, substances that assign different cell fates at different concentrations. Gradients form by morphogen transport from a localized site, but whether this occurs by simple diffusion or by more elaborate mechanisms is unclear. We attempt to resolve this controversy by analyzing recent data in ways that appropriately capture the complexity of systems in which transport, receptor interaction, endo- and exocytosis, and degradation occur together. We find that diffusive mechanisms of morphogen transport are much more plausible-and nondiffusive mechanisms much less plausible-than has generally been argued. Moreover, we show that a class of experiments, endocytic blockade, thought to effectively distinguish between diffusive and nondiffusive transport models actually fails to draw useful distinctions.

Selective pressures for and against genetic instability in cancer: a time-dependent problem.

Genetic instability in cancer is a two-edge sword. It can both increase the rate of cancer progression (by increasing the probability of cancerous mutations) and decrease the rate of cancer growth (by imposing a large death toll on dividing cells). Two of the many selective pressures acting upon a tumour, the need for variability and the need to minimize deleterious mutations, affect the tumour's 'choice' of a stable or unstable 'strategy'. As cancer progresses, the balance of the two pressures will change. In this paper, we examine how the optimal strategy of cancerous cells is shaped by the changing selective pressures. We consider the two most common patterns in multistage carcinogenesis: the activation of an oncogene (a one-step process) and an inactivation of a tumour-suppressor gene (a two-step process). For these, we formulate an optimal control problem for the mutation rate in cancer cells. We then develop a method to find optimal time-dependent strategies. It turns out that for a wide range of parameters, the most successful strategy is to start with a high rate of mutations and then switch to stability. This agrees with the growing biological evidence that genetic instability, prevalent in early cancers, turns into stability later on in the progression. We also identify parameter regimes where it is advantageous to keep stable (or unstable) constantly throughout the growth.

Replication-dependent size reduction precedes differentiation in Chlamydia trachomatis.

Chlamydia trachomatis is the most common cause of bacterial sexually transmitted infection. It produces an unusual intracellular infection in which a vegetative form, called the reticulate body (RB), replicates and then converts into an elementary body (EB), which is the infectious form. Here we use quantitative three-dimensional electron microscopy (3D EM) to show that C. trachomatis RBs divide by binary fission and undergo a sixfold reduction in size as the population expands. Conversion only occurs after at least six rounds of replication, and correlates with smaller RB size. These results suggest that RBs only convert into EBs below a size threshold, reached by repeatedly dividing before doubling in size. A stochastic mathematical model shows how replication-dependent RB size reduction produces delayed and asynchronous conversion, which are hallmarks of the Chlamydia developmental cycle. Our findings support a model in which RB size controls the timing of RB-to-EB conversion without the need for an external signal.

Optimal Proliferation and Differentiation of Chlamydia Trachomatis.

Chlamydia trachomatis is a bacterium that causes eye infection and blindness in humans. It has an unusual life cycle involving two developmental forms. Within a cytoplasmic inclusion, the reticulate body (RB) repeatedly divides by binary fission and asynchronously differentiates into the infectious elementary body (EB). Upon the death of the mammalian cell that host many such inclusions, only the EB form of the bacteria survive and proceed to infect other cells. Given the bacteria's fast spreading infection, conventional wisdom would have the few initial EB turn into RB, divide and proliferate first, and then eventually start converting in order to maximize the terminal EB population upon host cell lysis. Several biological processes are seen as possible mechanisms for implementing such a conversion strategy. However, the optimality of an instinctual strategy with a period of proliferate without conversion prior to the onset of differentiation has never been substantiated theoretically or justified mathematically. This paper formulates three relatively simple models that capture the essential features of the Chlamydia life cycle. When the initial infection is caused by the endocytosis of a small EB population well below the carrying capacity of the host cell, the Maximum Principle requires for these models an optimal conversion strategy that confirms and rigorously justifies the prevailing view of no conversion at the early stage of the host cell infection. However, the conventional supposition is found to be inappropriate for an initial EB (-to-RB) population near or above the carrying capacity. Previously suggested and new biological mechanisms are examined for their role in implementing the different optimal conversion strategies associated with models investigated herein.

TRANSIENT FEEDBACK AND ROBUST SIGNALING GRADIENTS.

Robust development of biological organisms in the presence of genetic and epi-genetic perturbations is important for time spans short relative to evolutionary time. Gradients of receptor bound signaling morphogens are responsible for patterning formation and development. A variety of inhibitors for reducing ectopic signaling activities are known to exist and their specific role in down-regulating the undesirable ectopic activities reasonably well understood. However, how a developing organism manages to adjust inhibition/stimulation in response to genetic and/or environmental changes remains to be uncovered. The need to adjust for ectopic signaling activities requires the presence of one or more feedback mechanisms to stimulate the needed adjustment. As the ultimate effect of many inhibitors (including those of the nonreceptor type) is to reduce the availability of signaling morphogens for binding with signaling receptors, a negative feedback on signaling morphogen synthesis rate based on a root-mean-square measure of the spatial distribution of signaling concentration offers a simple approach to robusness and has been demonstrated to be effective in a proof-of-concept implementation. In this paper, we complement the previous investigation of feedback in steady state by examining the effect of one or more feedback adjustments during the transient phase of the biological development.

Robust and precise morphogen-mediated patterning: trade-offs, constraints and mechanisms.

The patterning of many developing tissues is organized by morphogens. Genetic and environmental perturbations of gene expression, protein synthesis and ligand binding are among the sources of unreliability that limit the accuracy and precision of morphogen-mediated patterning. While it has been found that the robustness of morphogen gradients to the perturbation of morphogen synthesis can be enhanced by particular mechanisms, how such mechanisms affect robustness to other perturbations, such as to receptor synthesis for the same morphogen, has been little explored. Here, we investigate the interplay between the robustness of patterning to the changes in receptor synthesis and morphogen synthesis and to the effects of cell-to-cell variability. Our analysis elucidates the trade-offs and constraints that arise as a result of achieving these three performance objectives simultaneously in the context of simple, steady-state morphogen gradients formed by diffusion and receptor-mediated uptake. Analysis of the interdependence between length scales of patterning and these performance objectives reveals several potential mechanisms for mitigating such trade-offs and constraints. One involves downregulation of receptor synthesis in the morphogen source, while another involves the presence of non-signalling cell-surface morphogen-binding molecules. Both of these mechanisms occur in Drosophila wing discs during their patterning. We computationally elucidate how these mechanisms improve the robustness and precision of morphogen-mediated patterning.

On the behavior of solutions in viral dynamical models.

We consider simple mathematical models for the early population dynamics of the human immunodefficiency type 1 virus (HIV-1). Although these systems of differential equations may be solved by numerical methods, few general theoretical results are available due to nonlinearities. We analyze a model whose components are plasma densities of uninfected CD4+ T-cells and infected cells (assumed in this model to be proportional to virion density). In addition to analyzing the nature of the equilibrium points, we show that there are no periodic or limit-cycle solutions. Depending on the values of the parameters, solutions either tend without oscillation to an equilibrium point with zero virion density or to an equilibrium point in which there are a nonzero number of virions. In the latter case the approach to equilibrium may be through damped oscillations or without oscillation.

Cell-Surface Bound Nonreceptors and Signaling Morphogen Gradients.

The patterning of many developing tissues is orchestrated by gradients of signaling morphogens. Included among the molecular events that drive the formation of morphogen gradients are a variety of elaborate regulatory interactions. Such interactions are thought to make gradients robust, i.e. insensitive to change in the face of genetic or environmental perturbations. But just how this is accomplished is a major unanswered question. Recently extensive numerical simulations suggest that robustness of signaling gradients can be achieved through morphogen degradation mediated by cell surface bound non-signaling receptor molecules (or nonreceptors for short) such as heparan sulfate proteoglycans (HSPG). The present paper provides a mathematical validation of the results from the aforementioned numerical experiments. Extension of a basic extracellular model to include reversible binding with nonreceptors synthesized at a prescribed rate and mediated morphogen degradation shows that the signaling gradient diminishes with increasing concentration of cell-surface nonreceptors. Perturbation and asymptotic solutions obtained for i) low (receptor and nonreceptor) occupancy, and ii) high nonreceptor concntration permit more explicit delineation of the effects of nonreceptors on signaling gradients and facilitate the identification of scenarios in which the presence of nonreceptors may or may not be effective in promoting robustness.

The measure of success: constraints, objectives, and tradeoffs in morphogen-mediated patterning.

A large, diverse, and growing number of strategies have been proposed to explain how morphogen gradients achieve robustness and precision. We argue that, to be useful, the evaluation of such strategies must take into account the constraints imposed by competing objectives and performance tradeoffs. This point is illustrated through a mathematical and computational analysis of the strategy of self-enhanced morphogen clearance. The results suggest that the usefulness of this strategy comes less from its ability to increase robustness to morphogen source fluctuations per se, than from its ability to overcome specific kinds of noise, and to increase the fraction of a morphogen gradient within which robust threshold positions may be established. This work also provides new insights into the longstanding question of why morphogen gradients show a maximum range in vivo.

Feedback regulation in multistage cell lineages.

Studies of developing and self-renewing tissues have shown that differentiated cell types are typically specified through the actions of multistage cell lineages. Such lineages commonly include a stem cell and multiple progenitor (transit amplifying; TA) cell stages, which ultimately give rise to terminally differentiated (TD) cells. In several cases, self-renewal and differentiation of stem and progenitor cells within such lineages have been shown to be under feedback regulation. Together, the existence of multiple cell stages within a lineage and complex feedback regulation are thought to confer upon a tissue the ability to autoregulate development and regeneration, in terms of both cell number (total tissue volume) and cell identity (the proportions of different cell types, especially TD cells, within the tissue). In this paper, we model neurogenesis in the olfactory epithelium (OE) of the mouse, a system in which the lineage stages and mediators of feedback regulation that govern the generation of terminally differentiated olfactory receptor neurons (ORNs) have been the subject of much experimental work. Here we report on the existence and uniqueness of steady states in this system, as well as local and global stability of these steady states. In particular, we identify parameter conditions for the stability of the system when negative feedback loops are represented either as Hill functions, or in more general terms. Our results suggest that two factors -- autoregulation of the proliferation of transit amplifying (TA) progenitor cells, and a low death rate of TD cells -- enhance the stability of this system.

Compact integration factor methods in high spatial dimensions.

The dominant cost for integration factor (IF) or exponential time differencing (ETD) methods is the repeated vector-matrix multiplications involving exponentials of discretization matrices of differential operators. Although the discretization matrices usually are sparse, their exponentials are not, unless the discretization matrices are diagonal. For example, a two-dimensional system of N × N spatial points, the exponential matrix is of a size of N(2) × N(2) based on direct representations. The vector-matrix multiplication is of O(N(4)), and the storage of such matrix is usually prohibitive even for a moderate size N. In this paper, we introduce a compact representation of the discretized differential operators for the IF and ETD methods in both two- and three-dimensions. In this approach, the storage and CPU cost are significantly reduced for both IF and ETD methods such that the use of this type of methods becomes possible and attractive for two- or three-dimensional systems. For the case of two-dimensional systems, the required storage and CPU cost are reduced to O(N(2)) and O(N(3)), respectively. The improvement on three-dimensional systems is even more significant. We analyze and apply this technique to a class of semi-implicit integration factor method recently developed for stiff reaction-diffusion equations. Direct simulations on test equations along with applications to a morphogen system in two-dimensions and an intra-cellular signaling system in three-dimensions demonstrate an excellent efficiency of the new approach.

The role of feedback in the formation of morphogen territories.

In this paper, we consider a mathematical model for the formation of spatial morphogen territories of two key morphogens: Wingless (Wg) and Decapentaplegic (DPP), involved in leg development of Drosophila. We define a gene regulatory network (GRN) that utilizes autoactivation and cros-sinhibition (modeled by Hill equations) to establish and maintain stable boundaries of gene expression. By computational analysis we find that in the presence of a general activator, neither autoactivation, nor cross-inhibition alone are sufficient to maintain stable sharp boundaries of morphogen production in the leg disc. The minimal requirements for a self-organizing system are a coupled system of two morphogens in which the autoactivation and cross-inhibition have Hill coefficients strictly greater than one. In addition, the GRN modeled here describes the regenerative responses to genetic manipulations of positional identity in the leg disc.