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.

Catalytic mechanism of the adenylyl and guanylyl cyclases: Modeling and mutational analysis

The adenylyl and guanylyl cyclases catalyze the formation of 3′,5′-cyclic adenosine or guanosine monophosphate from the corresponding nucleoside 5′-triphosphate. The guanylyl cyclases, the mammalian adenylyl cyclases, and their microbial homologues function as pairs of homologous catalytic domains. The crystal structure of the rat type II adenylyl cyclase C2 catalytic domain was used to model by homology a mammalian adenylyl cyclase C1-C2 domain pair, a homodimeric adenylyl cyclase of Dictyostelium discoideum, a heterodimeric soluble guanylyl cyclase, and a homodimeric membrane guanylyl cyclase. Mg2+ATP or Mg2+GTP were docked into the active sites based on known stereochemical constraints on their conformation. The models are consistent with the activities of seven active-site mutants. Asp-310 and Glu-432 of type I adenylyl cyclase coordinate a Mg2+ ion. The D310S and D310A mutants have 10-fold reduced Vmax and altered [Mg2+] dependence. The NTP purine moieties bind in mostly hydrophobic pockets. Specificity is conferred by a Lys and an Asp in adenylyl cyclase, and a Glu, an Arg, and a Cys in guanylyl cyclase. The models predict that an Asp from one domain is a general base in the reaction, and that the transition state is stabilized by a conserved Asn-Arg pair on the other domain.

Ultrafast detection and control of molecular dynamics

Many elementary chemical and physical processes such as the breaking of a chemical bond or the vibrational motion of atoms within a molecule take place on a femtosecond (fs = 10−15 s) or picosecond (ps = 10−12 s) time scale. It is now possible to monitor these events as a function of time with temporal resolution well below 100 fs. This capability is based on the pump-probe technique where one optical pulse triggers a reaction and a second delayed optical pulse probes the changes that ensue. To illustrate this capability, the dynamics of ligand motion within a protein are presented. Moving beyond casual observation of a reaction to active control of its outcome requires additional experimental and theoretical effort. To illustrate the concept of control, the effect of optical pulse duration on the vibrational dynamics of a tri-atomic molecule are discussed. The experimental and theoretical resources currently available are poised to make the dream of reaction control a reality for certain molecular systems.

The neural development and organization of letter recognition: Evidence from functional neuroimaging, computational modeling, and behavioral studies

Although much of the brain’s functional organization is genetically predetermined, it appears that some noninnate functions can come to depend on dedicated and segregated neural tissue. In this paper, we describe a series of experiments that have investigated the neural development and organization of one such noninnate function: letter recognition. Functional neuroimaging demonstrates that letter and digit recognition depend on different neural substrates in some literate adults. How could the processing of two stimulus categories that are distinguished solely by cultural conventions become segregated in the brain? One possibility is that correlation-based learning in the brain leads to a spatial organization in cortex that reflects the temporal and spatial clustering of letters with letters in the environment. Simulations confirm that environmental co-occurrence does indeed lead to spatial localization in a neural network that uses correlation-based learning. Furthermore, behavioral studies confirm one critical prediction of this co-occurrence hypothesis, namely, that subjects exposed to a visual environment in which letters and digits occur together rather than separately (postal workers who process letters and digits together in Canadian postal codes) do indeed show less behavioral evidence for segregated letter and digit processing.

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.

Auxin-Growth Relationships in Maize Coleoptiles and Pea Internodes and Control by Auxin of the Tissue Sensitivity to Auxin

Growth of a zone of maize (Zea mays L.) coleoptiles and pea (Pisum sativum L.) internodes was greatly suppressed when the organ was decapitated or ringed at an upper position with the auxin transport inhibitor N-1-naphthylphthalamic acid (NPA) mixed with lanolin. The transport of apically applied 3H-labeled indole-3-acetic acid (IAA) was similarly inhibited by NPA. The growth suppressed by NPA or decapitation was restored by the IAA mixed with lanolin and applied directly to the zone, and the maximal capacity to respond to IAA did not change after NPA treatment, although it declined slightly after decapitation. The growth rate at IAA saturation was greater than the rate in intact, nontreated plants. It was concluded that growth is limited and controlled by auxin supplied from the apical region. In maize coleoptiles the sensitivity to IAA increased more than 3 times when the auxin level was reduced over a few hours with NPA treatment. This result, together with our previous result that the maximal capacity to respond to IAA declines in pea internodes when the IAA level is enhanced for a few hours, indicates that the IAA concentration-response relationship is subject to relatively slow adaptive regulation by IAA itself. The spontaneous growth recovery observed in decapitated maize coleoptiles was prevented by an NPA ring placed at an upper position of the stump, supporting the view that recovery is due to regenerated auxin-producing activity. The sensitivity increase also appeared to participate in an early recovery phase, causing a growth rate greater than in intact plants.

Independence of metastatic ability and extravasation: metastatic ras-transformed and control fibroblasts extravasate equally well.

Escape of cancer cells from the circulation (extravasation) is thought to be a major rate-limiting step in metastasis, with few cells being able to extravasate. Furthermore, highly metastatic cells are believed to extravasate more readily than poorly metastatic cells. We assessed in vivo the extravasation ability of highly metastatic ras-transformed NIH 3T3 cells (PAP2) versus control nontumorigenic nontransformed NIH 3T3 cells and primary mouse embryo fibroblasts. Fluorescently labeled cells were injected intravenously into chicken embryo chorioallantoic membrane and analyzed by intravital videomicroscopy. The chorioallantoic membrane is an appropriate model for studying extravasation, since, at the embryonic stage used, the microvasculature exhibits a continuous basement membrane and adult permeability properties. The kinetics of extravasation were assessed by determining whether individual cells (n = 1481) were intravascular, extravascular, or in the process of extravasation, at 3, 6, and 24 h after injection. Contrary to expectations, our results showed that all three cell types extravasated with the same kinetics. By 24 h after injection > 89% of observed cells had completed extravasation from the capillary plexus. After extravasation, individual fibroblasts of all cell types demonstrated preferential migration within the mesenchymal layer toward arterioles, not to venules or lymphatics. Thus in this model and for these cells, extravasation is independent of metastatic ability. This suggests that the ability to extravasate in vivo is not necessarily predictive of subsequent metastasis formation, and that postextravasation events may be key determinants in metastasis.

Modeling and optimization of populations subject to time-dependent mutation.

It has become clear that many organisms possess the ability to regulate their mutation rate in response to environmental conditions. So the question of finding an optimal mutation rate must be replaced by that of finding an optimal mutation schedule. We show that this task cannot be accomplished with standard population-dynamic models. We then develop a "hybrid" model for populations experiencing time-dependent mutation that treats population growth as deterministic but the time of first appearance of new variants as stochastic. We show that the hybrid model agrees well with a Monte Carlo simulation. From this model, we derive a deterministic approximation, a "threshold" model, that is similar to standard population dynamic models but differs in the initial rate of generation of new mutants. We use these techniques to model antibody affinity maturation by somatic hypermutation. We had previously shown that the optimal mutation schedule for the deterministic threshold model is phasic, with periods of mutation between intervals of mutation-free growth. To establish the validity of this schedule, we now show that the phasic schedule that optimizes the deterministic threshold model significantly improves upon the best constant-rate schedule for the hybrid and Monte Carlo models.