Cellular dynamic simulator: an event driven molecular simulation environment for cellular physiology.

In this paper, we present the Cellular Dynamic Simulator (CDS) for simulating diffusion and chemical reactions within crowded molecular environments. CDS is based on a novel event driven algorithm specifically designed for precise calculation of the timing of collisions, reactions and other events for each individual molecule in the environment. Generic mesh based compartments allow the creation / importation of very simple or detailed cellular structures that exist in a 3D environment. Multiple levels of compartments and static obstacles can be used to create a dense environment to mimic cellular boundaries and the intracellular space. The CDS algorithm takes into account volume exclusion and molecular crowding that may impact signaling cascades in small sub-cellular compartments such as dendritic spines. With the CDS, we can simulate simple enzyme reactions; aggregation, channel transport, as well as highly complicated chemical reaction networks of both freely diffusing and membrane bound multi-protein complexes. Components of the CDS are generally defined such that the simulator can be applied to a wide range of environments in terms of scale and level of detail. Through an initialization GUI, a simple simulation environment can be created and populated within minutes yet is powerful enough to design complex 3D cellular architecture. The initialization tool allows visual confirmation of the environment construction prior to execution by the simulator. This paper describes the CDS algorithm, design implementation, and provides an overview of the types of features available and the utility of those features are highlighted in demonstrations.

A Cbfa1-dependent genetic pathway controls bone formation beyond embryonic development.

The molecular mechanisms controlling bone extracellular matrix (ECM) deposition by differentiated osteoblasts in postnatal life, called hereafter bone formation, are unknown. This contrasts with the growing knowledge about the genetic control of osteoblast differentiation during embryonic development. Cbfa1, a transcriptional activator of osteoblast differentiation during embryonic development, is also expressed in differentiated osteoblasts postnatally. The perinatal lethality occurring in Cbfa1-deficient mice has prevented so far the study of its function after birth. To determine if Cbfa1 plays a role during bone formation we generated transgenic mice overexpressing Cbfa1 DNA-binding domain (DeltaCbfa1) in differentiated osteoblasts only postnatally. DeltaCbfa1 has a higher affinity for DNA than Cbfa1 itself, has no transcriptional activity on its own, and can act in a dominant-negative manner in DNA cotransfection assays. DeltaCbfa1-expressing mice have a normal skeleton at birth but develop an osteopenic phenotype thereafter. Dynamic histomorphometric studies show that this phenotype is caused by a major decrease in the bone formation rate in the face of a normal number of osteoblasts thus indicating that once osteoblasts are differentiated Cbfa1 regulates their function. Molecular analyses reveal that the expression of the genes expressed in osteoblasts and encoding bone ECM proteins is nearly abolished in transgenic mice, and ex vivo assays demonstrated that DeltaCbfa1-expressing osteoblasts were less active than wild-type osteoblasts. We also show that Cbfa1 regulates positively the activity of its own promoter, which has the highest affinity Cbfa1-binding sites characterized. This study demonstrates that beyond its differentiation function Cbfa1 is the first transcriptional activator of bone formation identified to date and illustrates that developmentally important genes control physiological processes postnatally.

Significance of Increased Tissue Transglutaminase in Hormone Refractory Prostate Cancer

The progression of hormone responsive to hormone refractory prostate cancer poses a major clinical challenge in the successful treatment of prostate cancer. The hormone refractory prostate cancer cells exhibit resistance not only to castrate levels of testosterone, but also to other therapeutic modalities and hence become lethal. Currently, there is no effective treatment available for managing this cancer. These observations underscore the urgency to investigate mechanism(s) that contribute to the progression of hormone-responsive to hormone-refractory prostate cancer and to target them for improved clinical outcomes. Tissue transglutaminase (TG2) is a multifunctional pro-inflammatory protein involved in diverse physiological processes such as inflammation, tissue repair, and wound healing. Its expression is also implicated in pathological conditions such as cancer and fibrosis. Interestingly, we found that the androgen-independent prostate cancer cell lines, which lacked androgen receptor (AR) expression, contained high basal levels of tissue transglutaminase. Inversely, the cell lines that expressed androgen receptor lacked transglutaminase expression. This attracted our attention to investigate the possible role this protein may play in the progression of prostate cancer, especially in view of recent observations that its expression is linked with increased invasion, metastasis, and drug resistance in multiple cancer cell types. The results we obtained were rather surprising and revealed that stable expression of tissue transglutaminase in androgen-sensitive LNCaP prostate cancer cells rendered these cells independent of androgen for growth and survival by silencing the AR expression. The AR silencing in TG2 expressing cells (TG2-infected LNCaP and PC-3 cells) was due to TG2-induced activation of the inflammatory nuclear transcription factor-kB (NF-kB). Thus, TG2 induced NF-kB was found to directly bind to the AR promoter. Importantly, TG2 protein was specifically recruited to the AR promoter in complex with the p65 subunit of NF-kB. Moreover, TG2 expressing LNCaP and PC-3 cells exhibited epithelial-to-mesenchymal transition, as evidenced by gain of mesenchymal (such as fibronectin, vimentin, etc.) and loss of epithelial markers (such as E-cadherin, b-catenin). Taken together, these results suggested a new function for TG2 and revealed a novel mechanism that is responsible for the progression of prostate cancer to the aggressive hormone-refractory phenotype.

Developmental programming of tumor suppressor gene penetrance

Epidemiological studies have led to the hypothesis that major risk factors for developing diseases such as hypertension, cardiovascular disease and adult-onset diabetes are established during development. This developmental programming hypothesis proposes that exposure to an adverse stimulus or insult at critical, sensitive periods of development can induce permanent alterations in normal physiological processes that lead to increased disease risk later in life. For cancer, inheritance of a tumor suppressor gene defect confers a high relative risk for disease development. However, these defects are rarely 100% penetrant. Traditionally, gene-environment interactions are thought to contribute to the penetrance of tumor suppressor gene defects by facilitating or inhibiting the acquisition of additional somatic mutations required for tumorigenesis. The studies presented herein identify developmental programming as a distinctive type of gene-environment interaction that can enhance the penetrance of a tumor suppressor gene defect in adult life. Using rats predisposed to uterine leiomyoma due to a germ-line defect in one allele of the tuberous sclerosis complex 2 (Tsc-2) tumor suppressor gene, these studies show that early-life exposure to the xenoestrogen, diethylstilbestrol (DES), during development of the uterus increased tumor incidence, multiplicity and size in genetically predisposed animals, but failed to induce tumors in wild-type rats. Uterine leiomyomas are ovarian-hormone dependent tumors that develop from the uterine myometrium. DES exposure was shown to developmentally program the myometrium, causing increased expression of estrogen-responsive genes prior to the onset of tumors. Loss of function of the normal Tsc-2 allele remained the rate-limiting event for tumorigenesis; however, tumors that developed in exposed animals displayed an enhanced proliferative response to ovarian steroid hormones relative to tumors that developed in unexposed animals. Furthermore, the studies presented herein identify developmental periods during which target tissues are maximally susceptible to developmental programming. These data suggest that exposure to environmental factors during critical periods of development can permanently alter normal physiological tissue responses and thus lead to increased disease risk in genetically susceptible individuals. ^

Diffusely adherent Escherichia coli (DAEC) as a cause of acute diarrhea

Acute diarrhea is the most common medical problem in the developing countries. Infectious agents are responsible for a majority of cases of acute diarrhea. Knowing the cause of acute diarrhea is important to developing plans for disease prevention, control and therapy. Acute diarrhea is caused by many viruses, bacteria and parasites. ^ Travelers to developing countries of the world commonly develop diarrhea as a result of eating contaminated food or drinking contaminated water. About 30-50% of travelers who travel from industrialized countries like United States to the developing countries are at risk of developing diarrhea. High risk areas for travelers' diarrhea are Mexico, Latin America and Southeast Asia. Public restaurants are the common sites for exposure to this type of food-borne infectious disease in travelers. Food becomes contaminated when they are handled by people with fecal content on their hands. ^ The importance of Diffusely Adherent Escherichia Coli (DAEC) in travelers to these areas has not been well studied. Some of the studies looking at DAEC have shown the organism to be present in children without symptoms. Other studies have shown a relationship between DAEC infection and presence of symptoms. I have selected this topic because the patho-physiological processes in DAEC infection that allow intestinal and extra-intestinal infections to develop are not fully understood. DAEC related acute diarrhea is a relatively new topic of public health significance. There is a limited number of studies regarding the virulence and pathogenic mechanisms of DAEC. The presumed virulence factor of the organism is diffuse attachment to the intestinal lining of the infected host. However more research needs to be done to identify the pathogenic mechanisms and virulence factors associated with DAEC infection for better treatment planning and diarrhea prevention. ^

A novel mechanism of regulation of phosphatidylinositol 3 -kinase: Tyrosine phosphorylation of p85 residue 688

Phosphatidylinositol 3-kinase (PI3K) phosphorylates membrane constituent phosphatidylinositols, producing second messengers that link membrane bound receptor signals to cellular proliferation and survival. PI3K, a heterodimer consisting of a catalytic p110 subunit and a regulatory p85 subunit, can be activated through induced association with other signaling molecules. The p85 subunit serves to both stabilize and inactivate p110. The inhibitory activity of P85 is relieved by occupancy of the N terminal SH2 domain by phosphorylated tyrosine. PI3K becomes phosphorylated and activated subsequent to a variety of stimuli. Indeed, Src family kinases have been demonstrated to phosphorylate p85 at tyrosine 688, but the role of phosphorylation in PI3K function is unclear. We decided to evaluate the importance of tyrosine phosphorylation to PI3K activity. We demonstrate that tyrosine phosphorylated p85 is associated with a higher specific activity than is non-phosphorylated PI3K. Wild type p85 inhibits PI3K enzyme activity, a process accentuated by mutation of tyrosine 688 to alanine and reversed by mutation to aspartate which functions as a phosphotyrosine mimic in multiple systems. Strikingly, the Y688D mutation completely reverses the p85 inhibitory activity on cell viability and activation of downstream protein NFkB. We demonstrate that tyrosine phosphorylated Y688 or Y688D is sufficient to bind the p85 N terminal SH2 domain, either within full length p85 or in an isolated N terminal SH2 domain, suggesting the possibility of an intramolecular interaction between phosphorylated Y688 and the p85 N terminal SH2 domain that can relieve the p85-induced inhibition of p110. Further, we provide evidence that dephosphorylation of Y688 reduces phosphorylation-induced PI3K activity. We demonstrate that tyrosine phosphatase SHP-1 can physically associate with p85 in a SH2-mediated interaction with the C terminal tail of SHP-1. This association is concomitant with both p85 dephosphorylation and decreased PI3K activity. Altogether, our data suggests the phosphorylation state of p85 is the focal point of a novel mechanism for PI3K activity regulation. As PI3K has been shown to be involved in the vital physiological processes of cell proliferation and apoptosis, a thorough understanding of the regulation of this signaling protein may provide opportunities for the design of novel treatments for cancer. ^

Phenotypic and genotypic analysis of mouse hepatitis virus (JHM) complementation group A temperature-sensitive mutants

Temperature sensitive (ts) mutant viruses have helped elucidate replication processes in many viral systems. Several panels of replication-defective ts mutants in which viral RNA synthesis is abolished at the nonpermissive temperature (RNA$\sp{-})$ have been isolated for Mouse Hepatitis Virus, MHV (Robb et al., 1979; Koolen et al., 1983; Martin et al., 1988; Schaad et al., 1990). However, no one had investigated genetic or phenotypic relationships between these different mutant panels. In order to determine how the panel of MHV-JHM RNA$\sp{-}$ ts mutants (Robb et al., 1979) were genetically related to other described MHV RNA$\sp{-}$ ts mutants, the MHV-JHM mutants were tested for complementation with representatives from two different sets of MHV-A59 ts mutants (Koolen et al., 1983; Schaad et al., 1990). The three ts mutant panels together were found to comprise eight genetically distinct complementation groups. Of these eight complementation groups, three complementation classes are unique to their particular mutant panel; genetically equivalent mutants were not observed within the other two mutant panels. Two complementation groups were common to all three mutant panels. The three remaining complementation groups overlapped two of the three mutant sets. Mutants MHV-JHM tsA204 and MHV-A59 ts261 were shown to be within one of these overlapping complementation groups. The phenotype of the MHV-JHM mutants within this complementation class has been previously characterized (Leibowitz et al., 1982; Leibowitz et al, 1990). When these mutants were grown at the permissive temperature, then shifted up to the nonpermissive temperature at the start of RNA synthesis, genome-length RNA and leader RNA fragments accumulated, but no subgenomic mRNA was synthesized. MHV-A59 ts261 produced leader RNA fragments identical to those observed with MHV-JHM tsA204. Thus, these two MHV RNA$\sp{-}$ ts mutants that were genetically equivalent by complementation testing were phenotypically similar as well. Recombination frequencies obtained from crosses of MHV-A59 ts261 with several of the gene 1 MHV-A59 mutants indicated that the causal mutation(s) of MHV-A59 ts261 was located near the overlapping junction of ORF1a and ORF1b, in the 3$\sp\prime$ end of ORF1a, or the 5$\sp\prime$ end of ORF1b. Sequence analysis of this junction and 1400 nucleotides into the 5$\sp\prime$ end of ORF1b of MHV-A59 ts261 revealed one nucleotide change from the wildtype MHV-A59. This substitution at nucleotide 13,598 (A to G) was a silent mutation in the ORF1a reading frame, but resulted in an amino acid change in ORF1b gene product (I to V). This amino acid change would be expressed only in the readthrough translation product produced upon successful ribosome frameshifting. A revertant of MHV-A59 ts261 (R2) also retained this guanidine residue, but had a second substitution at nucleotide 14,475 in ORF1b. This mutation results in the substitution of valine for an isoleucine.^ The data presented here suggest that the mutation in MHV-A59 ts261 (nucleotide 13,598) would be responsible for the MHV-JHM complementation group A phenotype. A second-site reversion at nucleotide 14,475 may correct this defect in the revertant. Sequencing of gene 1 immediately upstream of nucleotide 13,296 and downstream of nucleotide 15,010 must be conducted to test this hypothesis. ^