A theory for controlling cell cycle dynamics using a reversibly binding inhibitor

We demonstrate, by using mathematical modeling of cell division cycle (CDC) dynamics, a potential mechanism for precisely controlling the frequency of cell division and regulating the size of a dividing cell. Control of the cell cycle is achieved by artificially expressing a protein that reversibly binds and inactivates any one of the CDC proteins. In the simplest case, such as the checkpoint-free situation encountered in early amphibian embryos, the frequency of CDC oscillations can be increased or decreased by regulating the rate of synthesis, the binding rate, or the equilibrium constant of the binding protein. In a more complex model of cell division, where size-control checkpoints are included, we show that the same reversible binding reaction can alter the mean cell mass in a continuously dividing cell. Because this control scheme is general and requires only the expression of a single protein, it provides a practical means for tuning the characteristics of the cell cycle in vivo.

Molecular dynamics of a grid-mounted molecular dipolar rotor in a rotating electric field

Classical molecular dynamics is applied to the rotation of a dipolar molecular rotor mounted on a square grid and driven by rotating electric field E(ν) at T ≃ 150 K. The rotor is a complex of Re with two substituted o-phenanthrolines, one positively and one negatively charged, attached to an axial position of Rh\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{2}^{4+}}}\end{equation*}\end{document} in a [2]staffanedicarboxylate grid through 2-(3-cyanobicyclo[1.1.1]pent-1-yl)malonic dialdehyde. Four regimes are characterized by a, the average lag per turn: (i) synchronous (a < 1/e) at E(ν) = |E(ν)| > Ec(ν) [Ec(ν) is the critical field strength], (ii) asynchronous (1/e < a < 1) at Ec(ν) > E(ν) > Ebo(ν) > kT/μ, [Ebo(ν) is the break-off field strength], (iii) random driven (a ≃ 1) at Ebo(ν) > E(ν) > kT/μ, and (iv) random thermal (a ≃ 1) at kT/μ > E(ν). A fifth regime, (v) strongly hindered, W > kT, Eμ, (W is the rotational barrier), has not been examined. We find Ebo(ν)/kVcm−1 ≃ (kT/μ)/kVcm−1 + 0.13(ν/GHz)1.9 and Ec(ν)/kVcm−1 ≃ (2.3kT/μ)/kVcm−1 + 0.87(ν/GHz)1.6. For ν > 40 GHz, the rotor behaves as a macroscopic body with a friction constant proportional to frequency, η/eVps ≃ 1.14 ν/THz, and for ν < 20 GHz, it exhibits a uniquely molecular behavior.

From residue matching patterns to protein folding topographies: General model and bovine pancreatic trypsin inhibitor

A coarse-grained model for protein-folding dynamics is introduced based on a discretized representation of torsional modes. The model, based on the Ramachandran map of the local torsional potential surface and the class (hydrophobic/polar/neutral) of each residue, recognizes patterns of both torsional conformations and hydrophobic-polar contacts, with tolerance for imperfect patterns. It incorporates empirical rates for formation of secondary and tertiary structure. The method yields a topological representation of the evolving local torsional configuration of the folding protein, modulo the basins of the Ramachandran map. The folding process is modeled as a sequence of transitions from one contact pattern to another, as the torsional patterns evolve. We test the model by applying it to the folding process of bovine pancreatic trypsin inhibitor, obtaining a kinetic description of the transitions between the contact patterns visited by the protein along the dominant folding pathway. The kinetics and detailed balance make it possible to invert the result to obtain a coarse topographic description of the potential energy surface along the dominant folding pathway, in effect to go backward or forward between a topological representation of the chain conformation and a topographical description of the potential energy surface governing the folding process. As a result, the strong structure-seeking character of bovine pancreatic trypsin inhibitor and the principal features of its folding pathway are reproduced in a reasonably quantitative way.

Molecular dynamics of MHC genesis unraveled by sequence analysis of the 1,796,938-bp HLA class I region

The intensely studied MHC has become the paradigm for understanding the architectural evolution of vertebrate multigene families. The 4-Mb human MHC (also known as the HLA complex) encodes genes critically involved in the immune response, graft rejection, and disease susceptibility. Here we report the continuous 1,796,938-bp genomic sequence of the HLA class I region, linking genes between MICB and HLA-F. A total of 127 genes or potentially coding sequences were recognized within the analyzed sequence, establishing a high gene density of one per every 14.1 kb. The identification of 758 microsatellite provides tools for high-resolution mapping of HLA class I-associated disease genes. Most importantly, we establish that the repeated duplication and subsequent diversification of a minimal building block, MIC-HCGIX-3.8–1-P5-HCGIV-HLA class I-HCGII, engendered the present-day MHC. That the currently nonessential HLA-F and MICE genes have acted as progenitors to today’s immune-competent HLA-ABC and MICA/B genes provides experimental evidence for evolution by “birth and death,” which has general relevance to our understanding of the evolutionary forces driving vertebrate multigene families.

Colocalization of cell division proteins FtsZ and FtsA to cytoskeletal structures in living Escherichia coli cells by using green fluorescent protein

In the current model for bacterial cell division, FtsZ protein forms a ring that marks the division plane, creating a cytoskeletal framework for the subsequent action of other proteins such as FtsA. This putative protein complex ultimately generates the division septum. Herein we report that FtsZ and FtsA proteins tagged with green fluorescent protein (GFP) colocalize to division-site ring-like structures in living bacterial cells in a visible space between the segregated nucleoids. Cells with higher levels of FtsZ–GFP or with FtsA–GFP plus excess wild-type FtsZ were inhibited for cell division and often exhibited bright fluorescent spiral tubules that spanned the length of the filamentous cells. This suggests that FtsZ may switch from a septation-competent localized ring to an unlocalized spiral under some conditions and that FtsA can bind to FtsZ in both conformations. FtsZ–GFP also formed nonproductive but localized aggregates at a higher concentration that could represent FtsZ nucleation sites. The general domain structure of FtsZ–GFP resembles that of tubulin, since the C terminus of FtsZ is not required for polymerization but may regulate polymerization state. The N-terminal portion of Rhizobium FtsZ polymerized in Escherichia coli and appeared to copolymerize with E. coli FtsZ, suggesting a degree of interspecies functional conservation. Analysis of several deletions of FtsA–GFP suggests that multiple segments of FtsA are important for its localization to the FtsZ ring.

Antigen-driven CD4+ T cell and HIV-1 dynamics: Residual viral replication under highly active antiretroviral therapy

Antigen-induced stimulation of the immune system can generate heterogeneity in CD4+ T cell division rates capable of explaining the temporal patterns seen in the decay of HIV-1 plasma RNA levels during highly active antiretroviral therapy. Posttreatment increases in peripheral CD4+ T cell counts are consistent with a mathematical model in which host cell redistribution between lymph nodes and peripheral blood is a function of viral burden. Model fits to patient data suggest that, although therapy reduces HIV replication below replacement levels, substantial residual replication continues. This residual replication has important consequences for long-term therapy and the evolution of drug resistance and represents a challenge for future treatment strategies.

In vivo telomere dynamics of human hematopoietic stem cells

Aging in vivo and cell division in vitro are associated with telomere shortening. Several lines of evidence suggest that telomere length may be a good predictor of the long term replicative capacity of cells. To investigate the natural fate of chromosome telomeres of hematopoietic stem cells in vivo, we measured the telomere length of peripheral blood granulocytes from 11 fully engrafted bone marrow transplant recipients and from their respective donors. In 10 of 11 donor–recipient pairs, the telomere length was significantly reduced in the recipient and the extent of reduction correlated inversely with the number of nucleated cells infused. These data provide internally controlled in vivo evidence that, concomitantly with their proliferation, hematopoietic stem cells lose telomere length; it is possible that, as a result, their proliferative potential is reduced. These findings must be taken into account when developing new protocols in which few stem cells are used for bone marrow transplantation or for gene therapy.

In vivo estimates of division and death rates of human T lymphocytes.

We present data on the decay, after radiotherapy, of naive and memory human T lymphocytes with stable chromosome damage. These data are analyzed in conjunction with existing data on the decay of naive and memory T lymphocytes with unstable chromosome damage and older data on unsorted lymphocytes. The analyses yield in vivo estimates for some life-history parameters of human T lymphocytes. Best estimates of proliferation rates have naive lymphocytes dividing once every 3.5 years and memory lymphocytes dividing once every 22 weeks. It appears that memory lymphocytes can revert to the naive phenotype, but only, on average, after 3.5 years in the memory class. The lymphocytes with stable chromosome damage decay very slowly, yielding surprisingly low estimates of their death rate. The estimated parameters are used in a simple mathematical model of the population dynamics of undamaged naive and memory lymphocytes. We use this model to illustrate that it is possible for the unprimed subset of a constantly stimulated clone to stay small, even when there is a large population of specific primed cells reverting to the unprimed state.