Hybrid stochastic simulations of intracellular reaction-diffusion systems.

With the observation that stochasticity is important in biological systems, chemical kinetics have begun to receive wider interest. While the use of Monte Carlo discrete event simulations most accurately capture the variability of molecular species, they become computationally costly for complex reaction-diffusion systems with large populations of molecules. On the other hand, continuous time models are computationally efficient but they fail to capture any variability in the molecular species. In this study a hybrid stochastic approach is introduced for simulating reaction-diffusion systems. We developed an adaptive partitioning strategy in which processes with high frequency are simulated with deterministic rate-based equations, and those with low frequency using the exact stochastic algorithm of Gillespie. Therefore the stochastic behavior of cellular pathways is preserved while being able to apply it to large populations of molecules. We describe our method and demonstrate its accuracy and efficiency compared with the Gillespie algorithm for two different systems. First, a model of intracellular viral kinetics with two steady states and second, a compartmental model of the postsynaptic spine head for studying the dynamics of Ca+2 and NMDA receptors.

THEORETICAL STUDIES ON THE METHODS OF RECONSTRUCTING PHYLOGENETIC TREES FROM DNA SEQUENCE DATA (MOLECULAR EVOLUTION, HOMINOID EVOLUTION, COMPUTER SIMULATION)

(1) A mathematical theory for computing the probabilities of various nucleotide configurations is developed, and the probability of obtaining the correct phylogenetic tree (model tree) from sequence data is evaluated for six phylogenetic tree-making methods (UPGMA, distance Wagner method, transformed distance method, Fitch-Margoliash's method, maximum parsimony method, and compatibility method). The number of nucleotides (m*) necessary to obtain the correct tree with a probability of 95% is estimated with special reference to the human, chimpanzee, and gorilla divergence. m* is at least 4,200, but the availability of outgroup species greatly reduces m* for all methods except UPGMA. m* increases if transitions occur more frequently than transversions as in the case of mitochondrial DNA. (2) A new tree-making method called the neighbor-joining method is proposed. This method is applicable either for distance data or character state data. Computer simulation has shown that the neighbor-joining method is generally better than UPGMA, Farris' method, Li's method, and modified Farris method on recovering the true topology when distance data are used. A related method, the simultaneous partitioning method, is also discussed. (3) The maximum likelihood (ML) method for phylogeny reconstruction under the assumption of both constant and varying evolutionary rates is studied, and a new algorithm for obtaining the ML tree is presented. This method gives a tree similar to that obtained by UPGMA when constant evolutionary rate is assumed, whereas it gives a tree similar to that obtained by the maximum parsimony tree and the neighbor-joining method when varying evolutionary rate is assumed. ^

Mechanisms by which nonsense codons downregulate T-cell receptor transcripts

Nonsense-mediated mRNA decay (NMD) is a quality control mechanism that degrades aberrant mRNAs harboring premature termination codons (PTCs). Two out of three T-cell receptor β (TCRβ) transcripts carry PTCs as a result of error-prone programmed rearrangements that occur at this locus during lymphocyte maturation. PTCs decrease TCRβ mRNA levels to a much greater extent than mRNAs transcribed from non-rearranging genes. This robust decrease in TCRβ mRNA levels is not a unique characteristic of the T-cell environment or the TCRβ promoter. The simplest explanation for this is that PTC-bearing TCRβ mRNAs elicit a stronger NMD response. An alternative explanation is NMD collaborates with another mechanism to dramatically decrease PTC-bearing TCRβ mRNA levels. ^ In my dissertation, I investigated the molecular mechanism behind the strong decrease in TCRβ mRNA levels triggered by PTCs. To determine the location of this response, I performed mRNA half-life analysis and found that PTCs elicited more rapid TCRβ mRNA decay in the nuclear fraction, not the cytoplasmic fraction. Although decay was restricted to the nuclear fraction, PTC-bearing TCRβ transcript levels were extremely low in the cytoplasm, a phenomenon that I named the nonsense-codon induced partitioning shift (NIPS). I established that NIPS shares several qualities with NMD, including its dependence on translation and NMD factors. Several lines of evidence suggested that NIPS results from PTCs eliciting retention of TCRβ transcripts in the nuclear fraction. This retention, as well as rapid TCRβ mRNA decay, most likely occurs in either the nucleoplasm or the outer nuclear membrane, based on analysis of nuclear and cytoplasmic markers in the highly purified nuclei I used for my studies. To further address the location of decay, I asked whether nuclear or cytoplasmic RNA decay factors mediated the destruction of PTC-bearing mRNAs. My results suggested that a nuclear component of the 3'-to-5' exosome, as well as an endonucleolytic activity, are involved in the destruction of PTC-containing TCRβ mRNAs. Individual endogenous NMD substrates had differential requirements for nuclear and cytoplasmic exonucleases. In summary, my results provide evidence that PTCs trigger multiple mechanisms involving multiple decay factors to remove and regulate mRNAs in mammalian cells. ^

Characterization of a Novel Role for the Tousled- like Kinase in Kinetochore Assembly and Function in Caenorhabditis elegans

Chromosome segregation is a critical step during cell division to avoid aneuploidy and promote proper organismal development. Correct sister chromatid positioning and separation during mitosis helps to achieve faithful transmission of genetic material to daughter cells. This prevents improper chromosome partitioning that can potentially result in extrachromosomal fragments, increasing the tumorigenic potential of the cells. The kinetochore is a protenaicious structure responsible for the initiation and orchestration of chromosome movement during mitosis. This highly conserved structure among eukaryotes is required for chromosome attachment to the mitotic spindle and failure to assemble the kinetochore results in aberrant chromosome segregation. Thus elucidating the mechanism of kinetochore assembly is important to have a better understanding of the regulation that controls chromosome segregation. Our previous work identified the C. elegans Tousled-like kinase (TLK-1) as a mitotic kinase and depletion of TLK-1 results in embryonic lethality, characterized by nuclei displaying poor mitotic chromosome alignment, lagging chromosome, and chromosome bridges during anaphase. Additionally, previous studies from our group revealed that TLK-1 is phosphorylated independently by Aurora B at serine 634, and by CHK-1 at threonine T610. The research presented herein reveals that both phosphorylated forms of TLK-1 associate with the kinetochore during mitosis. Moreover, by systematic depletion of kinetochore proteins, I uncovered that pTLK-1 is bona fide kinetochore component that is located at the outer kinetochore layer, influencing the microtubule-binding interface. I also demonstrated that TLK-1 is necessary for the kinetochore localization of the microtubule interacting proteins CLS-2 and LIS-1 and I show that embryos depleted of TLK-1 presented an aberrant twisted kinetochore pattern. Furthermore, I established that the inner kinetochore protein KNL-2 is an in vitro substrate of TLK-1 indicating a possible role of TLK-1 in regulating centromeric assembly. Collectively, these results suggest a novel role for the Tousled-like kinase in regulation of kinetochore assembly and microtubule dynamics and demonstrate the necessity of TLK-1 for proper chromosome segregation in C. elegans.

Factors associated with HPV16 status and survival among patients with squamous cell carcinoma of the oropharynx

Background: Squamous cell carcinoma of the oropharynx (SCCOP) is characterized by local tumor aggressiveness, high recurrence rates, a high incidence of second primary tumors, and medical comorbidities. Significant trends in demographic and clinical characteristics as well as survival among SCCOP patients have been observed over time, likely owing to the changing etiology of the disease. Human papillomavirus type 16 (HPV16) infection is associated with a significant proportion of these cancers. Biomarkers that may aid in identifying patients that are at higher risk of recurrence and death are important so that these patients may be followed more closely to improve their quality of life. ^ Study population and methods: The retrospective review (Specific Aim 2) included 3891 newly diagnosed, previously untreated patients presenting to our institution between 1955 and 2004. A total of 2299 patients treated at our institution were included in survival and recursive partitioning analysis. The prospective cohort study (Specific Aim 3) included 266 patients presenting to our institution between January 2006 and September 2009. ^ Results: The results from the retrospective review showed that over time, patients presented at younger ages and were more likely to have base of tongue/tonsil tumors and to be never/former smokers and moreover survival improved significantly over time. In survival and recursive partitioning analyses, the TNM staging system was efficient in prognosticating patients prior to 1995. However, in the recent decade, the TNM staging system was completely inadequate. The factors having the greatest positive effect on overall survival since 1995 were those common to HPV-associated SCCOP. The results from the prospective cohort study indicate that patients with high nodal stage and those with late stage disease have increased levels of pretreatment serum HPV DNA. ^ Conclusions: We saw a distinct improvement in survival among SCCOP patients over the past 50 years at our institution. The main factors contributing to this were changes in clinical characteristics, in particular surrogates for HPV status. The current TNM staging system for SCCOP is inadequate and incorporation of HPV status (and perhaps smoking status) is encouraged. Furthermore, although pretreatment circulating levels of HPV DNA was associated with higher N category and overall disease stage, it has limited utility as a marker for recurrence among SCCOP patients.^