Dynamics of a minimal model of interlocked positive and negative feedback loops of transcriptional regulation by cAMP-response element binding proteins.

cAMP-response element binding (CREB) proteins are involved in transcriptional regulation in a number of cellular processes (e.g., neural plasticity and circadian rhythms). The CREB family contains activators and repressors that may interact through positive and negative feedback loops. These loops can be generated by auto- and cross-regulation of expression of CREB proteins, via CRE elements in or near their genes. Experiments suggest that such feedback loops may operate in several systems (e.g., Aplysia and rat). To understand the functional implications of such feedback loops, which are interlocked via cross-regulation of transcription, a minimal model with a positive and negative loop was developed and investigated using bifurcation analysis. Bifurcation analysis revealed diverse nonlinear dynamics (e.g., bistability and oscillations). The stability of steady states or oscillations could be changed by time delays in the synthesis of the activator (CREB1) or the repressor (CREB2). Investigation of stochastic fluctuations due to small numbers of molecules of CREB1 and CREB2 revealed a bimodal distribution of CREB molecules in the bistability region. The robustness of the stable HIGH and LOW states of CREB expression to stochastic noise differs, and a critical number of molecules was required to sustain the HIGH state for days or longer. Increasing positive feedback or decreasing negative feedback also increased the lifetime of the HIGH state, and persistence of this state may correlate with long-term memory formation. A critical number of molecules was also required to sustain robust oscillations of CREB expression. If a steady state was near a deterministic Hopf bifurcation point, stochastic resonance could induce oscillations. This comparative analysis of deterministic and stochastic dynamics not only provides insights into the possible dynamics of CREB regulatory motifs, but also demonstrates a framework for understanding other regulatory processes with similar network architecture.

Simulation of Drosophila circadian oscillations, mutations, and light responses by a model with VRI, PDP-1, and CLK.

A model of Drosophila circadian rhythm generation was developed to represent feedback loops based on transcriptional regulation of per, Clk (dclock), Pdp-1, and vri (vrille). The model postulates that histone acetylation kinetics make transcriptional activation a nonlinear function of [CLK]. Such a nonlinearity is essential to simulate robust circadian oscillations of transcription in our model and in previous models. Simulations suggest that two positive feedback loops involving Clk are not essential for oscillations, because oscillations of [PER] were preserved when Clk, vri, or Pdp-1 expression was fixed. However, eliminating positive feedback by fixing vri expression altered the oscillation period. Eliminating the negative feedback loop in which PER represses per expression abolished oscillations. Simulations of per or Clk null mutations, of per overexpression, and of vri, Clk, or Pdp-1 heterozygous null mutations altered model behavior in ways similar to experimental data. The model simulated a photic phase-response curve resembling experimental curves, and oscillations entrained to simulated light-dark cycles. Temperature compensation of oscillation period could be simulated if temperature elevation slowed PER nuclear entry or PER phosphorylation. The model makes experimental predictions, some of which could be tested in transgenic Drosophila.

A reduced model clarifies the role of feedback loops and time delays in the Drosophila circadian oscillator.

Although several detailed models of molecular processes essential for circadian oscillations have been developed, their complexity makes intuitive understanding of the oscillation mechanism difficult. The goal of the present study was to reduce a previously developed, detailed model to a minimal representation of the transcriptional regulation essential for circadian rhythmicity in Drosophila. The reduced model contains only two differential equations, each with time delays. A negative feedback loop is included, in which PER protein represses per transcription by binding the dCLOCK transcription factor. A positive feedback loop is also included, in which dCLOCK indirectly enhances its own formation. The model simulated circadian oscillations, light entrainment, and a phase-response curve with qualitative similarities to experiment. Time delays were found to be essential for simulation of circadian oscillations with this model. To examine the robustness of the simplified model to fluctuations in molecule numbers, a stochastic variant was constructed. Robust circadian oscillations and entrainment to light pulses were simulated with fewer than 80 molecules of each gene product present on average. Circadian oscillations persisted when the positive feedback loop was removed. Moreover, elimination of positive feedback did not decrease the robustness of oscillations to stochastic fluctuations or to variations in parameter values. Such reduced models can aid understanding of the oscillation mechanisms in Drosophila and in other organisms in which feedback regulation of transcription may play an important role.

The nonsteroidal anti-inflammatory drug indomethacin induces heterogeneity in lipid membranes: potential implication for its diverse biological action.

BACKGROUND: The nonsteroidal anti-inflammatory drug (NSAID), indomethacin (Indo), has a large number of divergent biological effects, the molecular mechanism(s) for which have yet to be fully elucidated. Interestingly, Indo is highly amphiphilic and associates strongly with lipid membranes, which influence localization, structure and function of membrane-associating proteins and actively regulate cell signaling events. Thus, it is possible that Indo regulates diverse cell functions by altering micro-environments within the membrane. Here we explored the effect of Indo on the nature of the segregated domains in a mixed model membrane composed of dipalmitoyl phosphatidyl-choline (di16:0 PC, or DPPC) and dioleoyl phosphatidyl-choline (di18:1 PC or DOPC) and cholesterol that mimics biomembranes. METHODOLOGY/PRINCIPAL FINDINGS: Using a series of fluorescent probes in a fluorescence resonance energy transfer (FRET) study, we found that Indo induced separation between gel domains and fluid domains in the mixed model membrane, possibly by enhancing the formation of gel-phase domains. This effect originated from the ability of Indo to specifically target the ordered domains in the mixed membrane. These findings were further confirmed by measuring the ability of Indo to affect the fluidity-dependent fluorescence quenching and the level of detergent resistance of membranes. CONCLUSION/SIGNIFICANCE: Because the tested lipids are the main lipid constituents in cell membranes, the observed formation of gel phase domains induced by Indo potentially occurs in biomembranes. This marked Indo-induced change in phase behavior potentially alters membrane protein functions, which contribute to the wide variety of biological activities of Indo and other NSAIDs.

Dynamics of a minimal model of interlocked positive and negative feedback loops of transcriptional regulation by cAMP-response element binding proteins.

cAMP-response element binding (CREB) proteins are involved in transcriptional regulation in a number of cellular processes (e.g., neural plasticity and circadian rhythms). The CREB family contains activators and repressors that may interact through positive and negative feedback loops. These loops can be generated by auto- and cross-regulation of expression of CREB proteins, via CRE elements in or near their genes. Experiments suggest that such feedback loops may operate in several systems (e.g., Aplysia and rat). To understand the functional implications of such feedback loops, which are interlocked via cross-regulation of transcription, a minimal model with a positive and negative loop was developed and investigated using bifurcation analysis. Bifurcation analysis revealed diverse nonlinear dynamics (e.g., bistability and oscillations). The stability of steady states or oscillations could be changed by time delays in the synthesis of the activator (CREB1) or the repressor (CREB2). Investigation of stochastic fluctuations due to small numbers of molecules of CREB1 and CREB2 revealed a bimodal distribution of CREB molecules in the bistability region. The robustness of the stable HIGH and LOW states of CREB expression to stochastic noise differs, and a critical number of molecules was required to sustain the HIGH state for days or longer. Increasing positive feedback or decreasing negative feedback also increased the lifetime of the HIGH state, and persistence of this state may correlate with long-term memory formation. A critical number of molecules was also required to sustain robust oscillations of CREB expression. If a steady state was near a deterministic Hopf bifurcation point, stochastic resonance could induce oscillations. This comparative analysis of deterministic and stochastic dynamics not only provides insights into the possible dynamics of CREB regulatory motifs, but also demonstrates a framework for understanding other regulatory processes with similar network architecture.

High-resolution microarray analysis of chromosome 20q in human colon cancer metastasis model systems

Amplification of human chromosome 20q DNA is the most frequently occurring chromosomal abnormality detected in sporadic colorectal carcinomas and shows significant correlation with liver metastases. Through comprehensive high-resolution microarray comparative genomic hybridization and microarray gene expression profiling, we have characterized chromosome 20q amplicon genes associated with human colorectal cancer metastasis in two in vitro metastasis model systems. The results revealed increasing complexity of the 20q genomic profile from the primary tumor-derived cell lines to the lymph node and liver metastasis derived cell lines. Expression analysis of chromosome 20q revealed a subset of over expressed genes residing within the regions of genomic copy number gain in all the tumor cell lines, suggesting these are Chromosome 20q copy number responsive genes. Bases on their preferential expression levels in the model system cell lines and known biological function, four of the over expressed genes mapping to the common intervals of genomic copy gain were considered the most promising candidate colorectal metastasis-associated genes. Validation of genomic copy number and expression array data was carried out on these genes, with one gene, DNMT3B, standing out as expressed at a relatively higher levels in the metastasis-derived cell lines compared with their primary-derived counterparts in both the models systems analyzed. The data provide evidence for the role of chromosome 20q genes with low copy gain and elevated expression in the clonal evolution of metastatic cells and suggests that such genes may serve as early biomarkers of metastatic potential. The data also support the utility of the combined microarray comparative genomic hybridization and expression array analysis for identifying copy number responsive genes in areas of low DNA copy gain in cancer cells. ^

Concomitant Mycobacterium tuberculosis infection and type 2 diabetes: A descriptive model with global health implications

The type 2 diabetes (diabetes) pandemic is recognized as a threat to tuberculosis (TB) control worldwide. This secondary data analysis project estimated the contribution of diabetes to TB in a binational community on the Texas-Mexico border where both diseases occur. Newly-diagnosed TB patients > 20 years of age were prospectively enrolled at Texas-Mexico border clinics between January 2006 and November 2008. Upon enrollment, information regarding social, demographic, and medical risks for TB was collected at interview, including self-reported diabetes. In addition, self-reported diabetes was supported by blood-confirmation according to guidelines published by the American Diabetes Association (ADA). For this project, data was compared to existing statistics for TB incidence and diabetes prevalence from the corresponding general populations of each study site to estimate the relative and attributable risks of diabetes to TB. In concordance with historical sociodemographic data provided for TB patients with self-reported diabetes, our TB patients with diabetes also lacked the risk factors traditionally associated with TB (alcohol abuse, drug abuse, history of incarceration, and HIV infection); instead, the majority of our TB patients with diabetes were characterized by overweight/obesity, chronic hyperglycemia, and older median age. In addition, diabetes prevalence among our TB patients was significantly higher than in the corresponding general populations. Findings of this study will help accurately characterize TB patients with diabetes, thus aiding in the timely recognition and diagnosis of TB in a population not traditionally viewed as at-risk. We provide epidemiological and biological evidence that diabetes continues to be an increasingly important risk factor for TB.^

BIOLOGICAL MECHANISMS AND CLINICAL IMPLICATIONS OF BCR-ABL-INDUCED MITOCHONDRIAL OXIDATIVE STRESS AND CELL SURVIVAL IN CHRONIC MYELOID LEUKEMIA

Chronic myeloid leukemia (CML), a myeloproliferative disorder, represents approximately 15-20% of all adult leukemia. The development of CML is clearly linked to the constitutively active protein-tyrosine kinase BCR-ABL, which is encoded by BCR-ABL fusion gene as the result of chromosome 9/22 translocation (Philadelphia chromosome). Previous studies have demonstrated that oxidative stress-associated genetic, metabolic and biological alterations contribute to CML cell survival and drug refractory. Mitochondria and NAD(P)H oxidase (NOX) are the major sources of BCR-ABL-induced cellular reactive oxygen species (ROS) production. However, it is still unknown how CML cells maintain the altered redox status, while escaping from the persistent oxidative stress-induced cell death. Therefore, elucidation of the mechanisms by which CML cells cope with oxidative stress will provide new insights into CML leukemogenesis. The major goal of this study is to identify the survival factors protecting CML cells against oxidative stress and develop novel therapeutic strategies to overcome drug resistance. Several experimental models were used to test CML cell redox status and cellular sensitivity to oxidative stress, including BCR-ABL inducible cell lines, BCR-ABL stably transformed cell lines and BCR-ABL-expressing CML blast crisis cells with differential BCL-XL/BCL-2 expressions. Additionally, an artificial CML cell model with heterogenic BCL-XL/BCL-2 expression was established to assess the correlation between differential survival factor expression patterns and cell sensitivity to Imatinib and oxidative stress. In this study, BCL-XL and GSH have been identified as the major survival factors responsive to BCR-ABL-promoted cellular oxidative stress and play a dominant role in regulating the threshold of oxidative stress-induced apoptosis. Cell survival factors BCL-XL and BCL-2 differentially protect mitochondria under oxidative stress. BCL-XL is an essential survival factor in preventing excessive ROS-induced cell death while BCL-2 seems to play a relatively minor role. Furthermore, the redox modulating reagent β-phenethyl isothiocyanate (PEITC) has been found to efficiently deplete GSH and induce potent cell killing effects in drug-resistant CML cells. Combination of PEITC with BCL-XL/BCL2 inhibitor ABT737 or suppression of BCL-XL by BCR-ABL inhibitor Gleevec dramatically sensitizes CML cells to apoptosis. These results have suggested that elevation of BCL-XL and cellular GSH are important for the development of CML, and that redox-directed therapy is worthy of further clinical investigations in CML.

SOCIAL AND BIOLOGICAL DETERMINANTS OF SHORT-TERM GROWTH VELOCITY AMONG PRESCHOOLERS OF RURAL GUATEMALA: A PATH ANALYSIS APPROACH (WEIGHT-GAIN, IRON DEFICIENCY, GASTROINTESTINAL SYMPTOMS)

This dissertation develops and tests through path analysis a theoretical model to explain how socioeconomic, socioenvironmental, and biologic risk factors simultaneously influence each other to further produce short-term, depressed growth in preschoolers. Three areas of risk factors were identified: child's proximal environment, maturational stage, and biological vulnerability. The theoretical model represented both the conceptual framework and the nature and direction of the hypotheses. Original research completed in 1978-80 and in 1982 provided the background data. It was analyzed first by nested-analysis of variance, followed by path analysis. The study provided evidence of mild iron deficiency and gastrointestinal symptomatology in the etiology of depressed, short-term weight gain. Also, there was evidence suggesting that family resources for material and social survival significantly contribute to the variability of short-term, age-adjusted growth velocity. These results challenge current views of unifocal intervention, whether for prevention or control. For policy formulations, though, the mechanisms underlying any set of interlaced relationships must be decoded. Theoretical formulations here proposed should be reassessed under a more extensive research design. It is suggested that studies should be undertaken where social changes are actually in progress; otherwise, nutritional epidemiology in developing countries operates somewhere between social reality and research concepts, with little grasp of its real potential. The study stresses that there is a connection between substantive theory, empirical observation, and policy issues. ^

Development and Characterization of an in vitro Four-species Anaerobic Dental Biofilm Model

Dental caries is the most common chronic disease worldwide. It is characterized by the demineralization of tooth enamel caused by acid produced by cariogenic dental bacteria growing on tooth surfaces, termed bacterial biofilms. Cariogenesis is a complex biological process that is influence by multiple factors and is not attributed to a sole causative agent. Instead, caries is associated with multispecies microbial biofilm communities composed of some bacterial species that directly influence the development of a caries lesion and other species that are seemingly benign but must contribute to the community in an uncharacterized way. Clinical analysis of dental caries and its microbial populations is challenging due to many factors including low sensitivity of clinical measurement tools, variability in saliva chemistry, and variation in the microbiota. Our laboratory has developed an in vitro anaerobic biofilm model for dental carries to facilitate both clinical and basic research-based analyses of the multispecies dynamics and individual factors that contribute to cariogenicity. The rational for development of this system was to improve upon the current models that lack key elements. This model places an emphasis on physiological relevance and ease of maintenance and reproducibility. The uniqueness of the model is based on integrating four critical elements: 1) a biofilm community composed of four distinct and representative species typically associated with dental caries, 2) a semi-defined synthetic growth medium designed to mimic saliva, 3) physiologically relevant biofilm growth substrates, and 4) a novel biofilm reactor device designed to facilitate the maintenance and analysis. Specifically, human tooth sections or hydroxyapatite discs embedded into poly(methyl methacrylate) (PMMA) discs are incubated for an initial 24 hr in a static inverted removable substrate (SIRS) biofilm reactor at 37°C under anaerobic conditions in artificial saliva (CAMM) without sucrose in the presence of 1 X 106 cells/ml of each Actinomyces odontolyticus, Fusobacterium nucleatum, Streptococcus mutans, and Veillonella dispar. During days 2 and 3 the samples are maintained continually in CAMM with various exposures to 0.2% sucrose; all of the discs are transferred into fresh medium every 24 hr. To validate that this model is an appropriate in vitro representation of a caries-associated multispecies biofilm, research aims were designed to test the following overarching hypothesis: an in vitro anaerobic biofilm composed of four species (S. mutans, V. dispar, A. odontolyticus, and F. nucleatum) will form a stable biofilm with a community profile that changes in response to environmental conditions and exhibits a cariogenic potential. For these experiments the biofilms as described above were exposed on days 2 and 3 to either CAMM lacking sucrose (no sucrose), CAMM with 0.2% sucrose (constant sucrose), or were transferred twice a day for 1 hr each time into 0.2% sucrose (intermittent sucrose). Four types of analysis were performed: 1) fluorescence microscopy of biofilms stained with Syto 9 and hexidium idodine to determine the biofilm architecture, 2) quantitative PCR (qPCR) to determine the cell number of each species per cm2, 3) vertical scanning interferometry (VSI) to determine the cariogenic potential of the biofilms, and 4) tomographic pH imaging using radiometric fluorescence microscopy after exposure to pH sensitive nanoparticles to measure the micro-environmental pH. The qualitative and quantitative results reveal the expected dynamics of the community profile when exposed to different sucrose conditions and the cariogenic potential of this in vitro four-species anaerobic biofilm model, thus confirming its usefulness for future analysis of primary and secondary dental caries.

Biological effects of PEA3 expression in {\it HER-2\/}/{\it neu\/}-overexpressing human cancer cells and the molecular mechanisms of PEA3-mediated transcriptional repression on {\it HER-2\/}/{\it neu\/}

HER-2/neu is a receptor tyrosine kinase highly homologous with epidermal growth factor receptor. Overexpression and/or amplification of HER-2/neu has been implicated in the genesis of a number of human cancers, especially breast and ovarian cancers. Transcriptional upregulation has been shown to contribute significantly to the overexpression of this gene. Studies on the transcriptional regulation of HER-2/neu gene are important for understanding the mechanism of cell transformation and developing the therapeutic strategies to block HER-2/neu-mediated cancers. PEA3 is a DNA binding transcriptional factor and its consensus sequence exists on the HER-2/neu promoter. To examine the role of PEA3 in HER-2/neu expression and cell transformation, we transfected PEA3 into the human breast and ovarian cancer cells that overexpress HER-2/neu and showed that PEA3 dramatically represses HER-2/neu transcription. PEA3 suppresses the oncogenic neu-mediated transformation in mouse fibroblast NIH 3T3 cells. Expression of PEA3 selectively blocks the growth of human cancer cells that overexpress HER-2/neu and inhibits their colony formation. It does not occur in the cancer cells expressing basal level of HER-2/neu. Further studies in the orthotopic ovarian cancer model demonstrated that expression of PEA3 preferentially inhibits growth and tumor development of human cancer cells that overexpress HER-2/neu, the tumor-bearing mice survived significantly longer if treated by injection of the PEA3-liposome complex intraperitoneally. Immunoblotting and immunohistochemical analysis of the tumor tissues indicated that PEA3 mediates the tumor suppression activity through targeting HER-2/neu-p185. Thus, PEA3 is a negative regulator of HER-2/neu gene expression and functions as a tumor suppressor gene in the HER-2/neu-overexpressing human cancer cells.^ The molecular mechanisms of PEA3 mediated transcriptional repression were investigated. PEA3 binds specifically at the PEA3 site on HER-2/neu promoter and this promoter-binding is required for the PEA3 mediated transcriptional repression. Mutation of the PEA3 binding site on HER-2/neu promoter causes decreased transcriptional activity, indicating that the PEA3 binding site is an enhancer-like element in the HER-2/neu-overexpressing cells. We therefore hypothesized that in the HER-2/neu-overexpressing cells, PEA3 competes with a transactivator for binding to the PEA3 site, preventing the putative factor from activating the transcription of HER-2/neu. This hypothesis was supported by the data which demonstrate that PEA3 competes with another nuclear protein for binding to the HER-2/neu promoter in vitro, and expression of a truncated protein which encodes the DNA binding domain of PEA3 is sufficient to repress HER-2/neu transcription in the HER-2/neu-overexpressing human cancer cells. ^