The kinetics of mammalian gene expression

Identifying Key In Silico Knockout for Enhancement of Limonene Yield Through Dynamic Metabolic Modelling

Living cells display dynamic and complex behaviors. To understand their response and to infer novel insights not possible with traditional reductionist approaches, over the last few decades various computational modelling methodologies have been developed. In this chapter, we focus on modelling the dynamic metabolic response, using linear and nonlinear ordinary differential equations, of an engineered Escherichia coli MG1655 strain with plasmid pJBEI-6409 that produces limonene. We show the systems biology steps involved from collecting time-series data of living cells, to dynamic model creation and fitting the model with experimental responses using COPASI software.

Nuclear deformation regulates YAP dynamics in cancer associated fibroblasts

Cells cultured on stiff 2D substrates exert high intracellular force, resulting in mechanical deformation of their nuclei. This nuclear deformation (ND) plays a crucial role in the transport of Yes Associated Protein (YAP) from the cytoplasm to the nucleus. However, cells in vivo are in soft 3D environment with potentially much lower intracellular forces. Whether and how cells may deform their nuclei in 3D for YAP localization remains unclear. Here, by culturing human colon cancer associated fibroblasts (CAFs) on 2D, 2.5D, and 3D substrates, we differentiated the effects of stiffness, force, and ND on YAP localization. We found that nuclear translocation of YAP depends on the degree of ND irrespective of dimensionality, stiffness and total force. ND induced by the perinuclear force, not the total force, and nuclear membrane curvature correlate strongly with YAP activation. Immunostained slices of human tumors further supported the association between ND and YAP nuclear localization, suggesting ND as a potential biomarker for YAP activation in tumors. Additionally, we conducted quantitative analysis of the force dynamics of CAFs on 2D substrates to construct a stochastic model of YAP kinetics. This model revealed that the probability of YAP nuclear translocation, as well as the residence time in the nucleus follow a power law. This study provides valuable insights into the regulatory mechanisms governing YAP dynamics and highlights the significance of threshold activation in YAP localization. STATEMENT OF SIGNIFICANCE: Yes Associated Protein (YAP), a transcription cofactor, has been identified as one of the drivers of cancer progression. High tumor stiffness is attributed to driving YAP to the nucleus, wherein it activates pro-metastatic genes. Here we show, using cancer associated fibroblasts, that YAP translocation to the nucleus depends on the degree of nuclear deformation, irrespective of stiffness. We also identified that perinuclear force induced membrane curvature correlates strongly with YAP nuclear transport. A novel stochastic model of YAP kinetics unveiled a power law relationship between the activation threshold and persistence time of YAP in the nucleus. Overall, this study provides novel insights into the regulatory mechanisms governing YAP dynamics and the probability of activation that is of immense clinical significance.

Principles of Computation by Competitive Protein Dimerization Networks

Many biological signaling pathways employ proteins that competitively dimerize in diverse combinations. These dimerization networks can perform biochemical computations, in which the concentrations of monomers (inputs) determine the concentrations of dimers (outputs). Despite their prevalence, little is known about the range of input-output computations that dimerization networks can perform (their "expressivity") and how it depends on network size and connectivity. Using a systematic computational approach, we demonstrate that even small dimerization networks (3-6 monomers) are expressive, performing diverse multi-input computations. Further, dimerization networks are versatile, performing different computations when their protein components are expressed at different levels, such as in different cell types. Remarkably, individual networks with random interaction affinities, when large enough (≥8 proteins), can perform nearly all (~90%) potential one-input network computations merely by tuning their monomer expression levels. Thus, even the simple process of competitive dimerization provides a powerful architecture for multi-input, cell-type-specific signal processing.

Phototoxic effects of nonlinear optical microscopy on cell cycle, oxidative states, and gene expression

Nonlinear optical imaging modalities, such as stimulated Raman scattering (SRS) microscopy, use pulsed-laser excitation with high peak intensity that can perturb the native state of cells. In this study, we used bulk RNA sequencing, quantitative measurement of cell proliferation, and fluorescent measurement of the generation of reactive oxygen species to assess phototoxic effects of near-IR pulsed laser radiation, at different time scales, for laser excitation settings relevant to SRS imaging. We define a range of laser excitation settings for which there was no significant change in mouse Neuro2A cells after laser exposure. This study provides guidance for imaging parameters that minimize photo-induced perturbations in SRS microscopy to ensure accurate interpretation of experiments with time-lapse imaging or with paired measurements of imaging and sequencing on the same cells.

Expression of calpain system transcripts responds inversely to beef tenderization after vitamin D3 supplementation in Nellore cattle

ABSTRACT The calpain system is the central player for meat tenderization and the calpastatin activity plays an important role in beef tenderness of Bos taurus indicus cattle. This study investigated whether dietary vitamin D3 induced gene expression of μ-calpain (CAPN1), total calpastatin (CAST T), and their variants (CAST I and II). Animals received none or 2 × 106 IU of vitamin D3 for either 2 or 8 days before slaughter and were submitted to different conditions during feedlot: sun exposure or artificial shade. Steaks from Longissimus lumborum were fabricated, aged for 1, 7, and 21 days post-mortem, and later used for the analyses of shear force and the myofibrillar fragmentation index. Vitamin D3 did not influence mRNA abundance; however, it induced a greater CAST II transcript in animals supplemented 8 days before slaughter. There was a negative association between CAST II abundance and the shear force, which revealed an important modulation of the calpain system expression due to vitamin D supplementation. This result is an important tool for strategies to improve beef tenderness.

Transcriptional Riboswitches Integrate Timescales for Bacterial Gene Expression Control

Transcriptional riboswitches involve RNA aptamers that are typically found in the 5' untranslated regions (UTRs) of bacterial mRNAs and form alternative secondary structures upon binding to cognate ligands. Alteration of the riboswitch's secondary structure results in perturbations of an adjacent expression platform that controls transcription elongation and termination, thus turning downstream gene expression "on" or "off." Riboswitch ligands are typically small metabolites, divalent cations, anions, signaling molecules, or other RNAs, and can be part of larger signaling cascades. The interconnectedness of ligand binding, RNA folding, RNA transcription, and gene expression empowers riboswitches to integrate cellular processes and environmental conditions across multiple timescales. For a successful response to an environmental cue that may determine a bacterium's chance of survival, a coordinated coupling of timescales from microseconds to minutes must be achieved. This review focuses on recent advances in our understanding of how riboswitches affect such critical gene expression control across time.

Development and Analysis of a Quantitative Mathematical Model of Bistability in the Cross Repression System Between APT and SLBO Within the JAK/STAT Signaling Pathway

Cell migration is a key component in development, homeostasis, immune function, and pathology. It is important to understand the molecular activity that allows some cells to migrate. Drosophila melanogaster is a useful model system because its genes are largely conserved with humans and it is straightforward to study biologically. The well-conserved transcriptional regulator Signal Transducer and Activator of Transcription (STAT) promotes cell migration, but its signaling is modulated by downstream targets Apontic (APT) and Slow Border Cells (SLBO). Inhibition of STAT activity by APT and cross-repression of APT and SLBO determines whether an epithelial cell in the Drosophila egg chamber becomes motile or remains stationary. Through mathematical modeling and analysis, we examine how the interaction of STAT, APT, and SLBO creates bistability in the Janus Kinase (JAK)/STAT signaling pathway. In this paper, we update and analyze earlier models to represent mechanistically the processes of the JAK/STAT pathway. We utilize parameter, bifurcation, and phase portrait analyses, and make reductions to the system to produce a minimal three-variable quantitative model. We analyze the manifold between migratory and stationary steady states in this minimal model and show that when the initial conditions of our model are near this manifold, cell migration can be delayed.

Network inference performance complexity: a consequence of topological, experimental and algorithmic determinants

MOTIVATION

Network inference algorithms aim to uncover key regulatory interactions governing cellular decision-making, disease progression and therapeutic interventions. Having an accurate blueprint of this regulation is essential for understanding and controlling cell behavior. However, the utility and impact of these approaches are limited because the ways in which various factors shape inference outcomes remain largely unknown.

RESULTS

We identify and systematically evaluate determinants of performance-including network properties, experimental design choices and data processing-by developing new metrics that quantify confidence across algorithms in comparable terms. We conducted a multifactorial analysis that demonstrates how stimulus target, regulatory kinetics, induction and resolution dynamics, and noise differentially impact widely used algorithms in significant and previously unrecognized ways. The results show how even if high-quality data are paired with high-performing algorithms, inferred models are sometimes susceptible to giving misleading conclusions. Lastly, we validate these findings and the utility of the confidence metrics using realistic in silico gene regulatory networks. This new characterization approach provides a way to more rigorously interpret how algorithms infer regulation from biological datasets.

AVAILABILITY AND IMPLEMENTATION

Code is available at http://github.com/bagherilab/networkinference/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Nascent transcript analysis of glucocorticoid crosstalk with TNF defines primary and cooperative inflammatory repression

The glucocorticoid receptor (NR3C1, also known as GR) binds to specific DNA sequences and directly induces transcription of anti-inflammatory genes that contribute to cytokine repression, frequently in cooperation with NF-kB. Whether inflammatory repression also occurs through local interactions between GR and inflammatory gene regulatory elements has been controversial. Here, using global run-on sequencing (GRO-seq) in human airway epithelial cells, we show that glucocorticoid signaling represses transcription within 10 min. Many repressed regulatory regions reside within “hyper-ChIPable” genomic regions that are subject to dynamic, yet nonspecific, interactions with some antibodies. When this artifact was accounted for, we determined that transcriptional repression does not require local GR occupancy. Instead, widespread transcriptional induction through canonical GR binding sites is associated with reciprocal repression of distal TNF-regulated enhancers through a chromatin-dependent process, as evidenced by chromatin accessibility and motif displacement analysis. Simultaneously, transcriptional induction of key anti-inflammatory effectors is decoupled from primary repression through cooperation between GR and NF-kB at a subset of regulatory regions. Thus, glucocorticoids exert bimodal restraints on inflammation characterized by rapid primary transcriptional repression without local GR occupancy and secondary anti-inflammatory effects resulting from transcriptional cooperation between GR and NF-kB.

Homeostasis: How Neurons Achieve Temperature Invariance

Temperature influences physiological processes and can corrupt nervous system function. A modelling study shows how regulation of ion channel expression can establish an acute temperature invariance of neuronal responses despite temperature-dependent and variable ionic conductances.

Acute Noise Exposure Is Associated With Intrinsic Apoptosis in Murine Central Auditory Pathway

Noise that is capable of inducing the hearing loss (NIHL) has a strong impact on the inner ear structures and causes early and most obvious pathophysiological changes in the auditory periphery. Several studies indicated that intrinsic apoptotic cell death mechanisms are the key factors inducing cellular degeneration immediately after noise exposure and are maintained for days or even weeks. In addition, studies demonstrated several changes in the central auditory system following noise exposure, consistent with early apoptosis-related pathologies. To clarify the underlying mechanisms, the present study focused on the noise-induced gene and protein expression of the pro-apoptotic protease activating factor-1 (APAF1) and the anti-apoptotic B-cell lymphoma 2 related protein a1a (BCL2A1A) in the cochlear nucleus (CN), inferior colliculus (IC) and auditory cortex (AC) of the murine central auditory pathway. The expression of Bcl2a1a mRNA was upregulated immediately after trauma in all tissues investigated, whereas the protein levels were significantly reduced at least in the auditory brainstem. Conversely, acute noise has decreased the expression of Apaf1 gene along the auditory pathway. The changes in APAF1 protein level were not statistically significant. It is tempting to speculate that the acoustic overstimulation leads to mitochondrial dysfunction and induction of apoptosis by regulation of proapoptotic and antiapoptotic proteins. The inverse expression pattern on the mRNA level of both genes might reflect a protective response to decrease cellular damage. Our results indicate the immediate presence of intrinsic apoptosis following noise trauma. This, in turn, may significantly contribute to the development of central structural deficits. Auditory pathway-specific inhibition of intrinsic apoptosis could be a therapeutic approach for the treatment of acute (noise-induced) hearing loss to prevent irreversible neuronal injury in auditory brain structures and to avoid profound deficits in complex auditory processing.

Determining Genome-wide Transcript Decay Rates in Proliferating and Quiescent Human Fibroblasts

Quiescence is a temporary, reversible state in which cells have ceased cell division, but retain the capacity to proliferate. Multiple studies, including ours, have demonstrated that quiescence is associated with widespread changes in gene expression. Some of these changes occur through changes in the level or activity of proliferation-associated transcription factors, such as E2F and MYC. We have demonstrated that mRNA decay can also contribute to changes in gene expression between proliferating and quiescent cells. In this protocol, we describe the procedure for establishing proliferating and quiescent cultures of human dermal foreskin fibroblasts. We then describe the procedures for inhibiting new transcription in proliferating and quiescent cells with Actinomycin D (ActD). ActD treatment represents a straightforward and reproducible approach to dissociating new transcription from transcript decay. A disadvantage of ActD treatment is that the time course must be limited to a short time frame because ActD affects cell viability. Transcript levels are monitored over time to determine transcript decay rates. This procedure allows for the identification of genes and isoforms that exhibit differential decay in proliferating versus quiescent fibroblasts.

Synergistic Activation of RD29A Via Integration of Salinity Stress and Abscisic Acid in Arabidopsis thaliana

Plants perceive information from the surroundings and elicit appropriate molecular responses. How plants dynamically respond to combinations of external inputs is yet to be revealed, despite the detailed current knowledge of intracellular signaling pathways. We measured dynamics of Response-to-Dehydration 29A (RD29A) expression induced by single or combined NaCl and ABA treatments in Arabidopsis thaliana. RD29A expression in response to a combination of NaCl and ABA leads to unique dynamic behavior that cannot be explained by the sum of responses to individual NaCl and ABA. To explore the potential mechanisms responsible for the observed synergistic response, we developed a mathematical model of the DREB2 and AREB pathways based on existing knowledge, where NaCl and ABA act as the cognate inputs, respectively, and examined various system structures with cross-input modulation, where non-cognate input affects expression of the genes involved in adjacent signaling pathways. The results from the analysis of system structures, combined with the insights from microarray expression profiles and model-guided experiments, predicted that synergistic activation of RD29A originates from enhancement of DREB2 activity by ABA. Our analysis of RD29A expression profiles demonstrates that a simple mathematical model can be used to extract information from temporal dynamics induced by combinatorial stimuli and produce experimentally testable hypotheses.

Non-invasive prediction of implantation window in controlled hyperstimulation cycles: Can the time from the menstrual day at embryo transfer to expected menstrual cycle give a clue?

Objective:

The aim of this study was to assess whether the time from the menstrual day at embryo transfer to expected menstrual cycle (TETEMC) is associated with the implantation in women with regular cycles or not.

Materials and Methods:

Forty women with successful implantation and forty women without implantation with regular cycles were randomly selected from prospectively collected database of assisted reproductive technology clinic of Zeynep Kamil Women And Children’s Health Training and Research Hospital. TETEMC was calculated for each case to assess relationship with the successful implantation.

Results:

Comparison of groups revealed significant differences with regard to TETEMC and the menstrual period (p<0.05).

In ROC analyses both the TETEMC (AUC=0.824, p<0.001) and the menstrual period (AUC=0.797, p<0.001) were significant predictors for clinical pregnancy. Cut off value for the menstrual period was found to be 27.5 days with 82.6% sensitivity and 65% specificity. Cut off value for TETEMC was 11.5 days with 75% sensitivity and 63.2% specificity.

Conclusion:

Longer menstrual cycle and the TETEMC seem to be associated with the implantation failure.

Nucleotide excision repair is not induced in human embryonic lung fibroblasts treated with environmental pollutants

The cellular response to genotoxic treatment depends on the cell line used. Although tumor cell lines are widely used for genotoxicity tests, the interpretation of the results may be potentially hampered by changes in cellular processes caused by malignant transformation. In our study we used normal human embryonic lung fibroblasts (HEL12469 cells) and tested their response to treatment with benzo[a]pyrene (B[a]P) and extractable organic matter (EOM) from ambient air particles <2.5 µm (PM2.5) collected in two Czech cities differing in levels and sources of air pollution. We analyzed multiple endpoints associated with exposure to polycyclic aromatic hydrocarbons (PAHs) including the levels of bulky DNA adducts and the nucleotide excision repair (NER) response [expression of XPE, XPC and XPA genes on the level of mRNA and proteins, unscheduled DNA synthesis (UDS)]. EOMs were collected in the winter and summer of 2011 in two Czech cities with different levels and sources of air pollution. The effects of the studied compounds were analyzed in the presence (+S9) and absence (–S9) of the rat liver microsomal S9 fraction. The levels of bulky DNA adducts were highest after treatment with B[a]P, followed by winter EOMs; their induction by summer EOMs was weak. The induction of both mRNA and protein expression was observed, with the most pronounced effects after treatment with B[a]P (–S9); the response induced by EOMs from both cities and seasons was substantially weaker. The expression of DNA repair genes was not accompanied by the induction of UDS activity. In summary, our results indicate that the tested compounds induced low levels of DNA damage and affected the expression of NER genes; however, nucleotide excision repair was not induced.

Short-term, high-fat diet increases the expression of key intestinal genes involved in lipoprotein metabolism in healthy men

BACKGROUND

The modulation of cholesterol and fatty acid homeostasis by dietary fatty acids is thought to be mediated by changes in the expression of key intestinal genes involved in lipoprotein metabolism. However, the short-term effect of dietary fat intake on the expression of these genes has not been fully investigated in humans.

OBJECTIVE

To test whether short-term changes in dietary fatty acid intake affect the expression of key intestinal genes involved in lipoprotein metabolism, we conducted a randomized, double-blind, crossover study in 12 nonobese, healthy men with normal plasma lipid profiles.

DESIGN

Participants were subjected to the following 2 intensive 3-d dietary interventions under isocaloric conditions: 1) a high-fat diet (37% of energy from fat and 50% of energy from carbohydrates) and 2) a low-fat diet (25% of energy from fat and 62% of energy from carbohydrates). Expressions of key genes involved in lipoprotein metabolism were compared by using real-time polymerase chain reaction quantification on duodenal biopsy specimens obtained in a fasting state after each diet.

RESULTS

After the 3-d high-fat diet, plasma cholesterol, LDL cholesterol, and HDL cholesterol concentrations were significantly higher than concentrations observed after the low-fat diet was consumed. The high-fat diet also resulted in significant increases in the intestinal messenger RNA expression of several key genes involved in lipoprotein metabolism. Plasma triglycerides and apolipoprotein B-48 concentrations were significantly lower after the high-fat diet than after the low-fat diet.

CONCLUSION

These findings suggest that short-term exposure to a high-fat diet upregulates the expression of key genes involved in lipid and lipoprotein metabolism at the enterocyte level. This trial was registered at clinicaltrials.gov as NCT01806441.

Surface contamination of dental implants assessed by gene expression analysis in a whole-blood in vitro assay: a preliminary study

AIM

We aimed at evaluating pyrogen contamination of dental implants made of titanium and zirconia by using gene expression analysis in a whole-blood in vitro assay.

MATERIAL AND METHODS

Titanium and zirconia implants (five each) were incubated in human whole blood. Samples were assayed for gene expression levels of toll-like receptor 4 (TLR4), TLR9, interleukin (IL)-1β, nuclear factor 'kappa-light-chain-enhancer' of activated B-cells (NF-kB), tumour necrosis factor (TNF)-α, and Fas-associated protein with death domain (FADD) as indicators of surface contamination resulting in lipopolysaccharides (LPS)-stimulated TLR- or TNF-mediated immune responses. Gene expression was assayed using real-time quantitative polymerase chain reaction (RT-qPCR). Non-stimulated blood from the same donor served as a negative control, and blood stimulated with LPS served as a positive control. After dry-heat treatment with dry heat, all implants were re-analysed as described above.

RESULTS

Both implant systems contained surface contaminants evoking a pro-inflammatory response similar to that induced by LPS. After dry-heat treatment, gene expression was significantly decreased to levels similar to those of negative control samples.

CONCLUSIONS

The results demonstrated LPS-like surface-bound contaminants in both tested implant systems. Depyrogenation with dry heat seems to be an effective means of reducing such contamination in dental implants.

Bridging time scales in cellular decision making with a stochastic bistable switch

Background

Cellular transformations which involve a significant phenotypical change of the cell's state use bistable biochemical switches as underlying decision systems. Some of these transformations act over a very long time scale on the cell population level, up to the entire lifespan of the organism.

Results

In this work, we aim at linking cellular decisions taking place on a time scale of years to decades with the biochemical dynamics in signal transduction and gene regulation, occuring on a time scale of minutes to hours. We show that a stochastic bistable switch forms a viable biochemical mechanism to implement decision processes on long time scales. As a case study, the mechanism is applied to model the initiation of follicle growth in mammalian ovaries, where the physiological time scale of follicle pool depletion is on the order of the organism's lifespan. We construct a simple mathematical model for this process based on experimental evidence for the involved genetic mechanisms.

Conclusions

Despite the underlying stochasticity, the proposed mechanism turns out to yield reliable behavior in large populations of cells subject to the considered decision process. Our model explains how the physiological time constant may emerge from the intrinsic stochasticity of the underlying gene regulatory network. Apart from ovarian follicles, the proposed mechanism may also be of relevance for other physiological systems where cells take binary decisions over a long time scale.

Toward a system-level understanding of microRNA pathway via mathematical modeling

The microRNA (miRNA) pathway plays multiple roles in regulating mechanisms controlling both physiological and pathological processes such as the cell proliferation and cancers. But little is known about the dynamic properties, key rate-limiting steps as well as the stochastic noise in this pathway. Presently, a system-theoretic approach was presented to analyze and quantitative modeling of a generic miRNA pathway, which can be implemented deterministically and stochastically. Our results show that the inferred dynamic properties obtained from the mathematical models of the miRNA pathway are well consistent with previous experimental observations. By sensitivity analysis, the key steps in this pathway were found to be the miRNA gene transcription, RISC decay and mRNA formation. In addition, the results of quantified noise strength along the pathway demonstrate that the pathway can reduce the ingress noise and reveal the noise robustness property. Our findings also present testable hypothesis for experimental biologists to further investigate miRNA's increasing functional roles in regulating various cellular processes.

Small RNAs establish delays and temporal thresholds in gene expression

Noncoding RNAs are crucial regulators of gene expression in prokaryotes and eukaryotes, but how they affect the dynamics of transcriptional networks remains poorly understood. We analyzed the temporal characteristics of the cyanobacterial iron stress response by mathematical modeling and quantitative experimental analyses and focused on the role of a recently discovered small noncoding RNA, IsrR. We found that IsrR is responsible for a pronounced delay in the accumulation of isiA mRNA encoding the late-phase stress protein, IsiA, and that it ensures a rapid decline in isiA levels once external stress triggers are removed. These kinetic properties allow the system to selectively respond to sustained (as opposed to transient) stimuli and thus establish a temporal threshold, which prevents energetically costly IsiA accumulation under short-term stress conditions. Biological information is frequently encoded in the quantitative aspects of intracellular signals (e.g., amplitude and duration). Our simulations reveal that competitive inhibition and regulated degradation allow intracellular regulatory networks to efficiently discriminate between transient and sustained inputs.

Phenotype prediction in regulated metabolic networks

Background

Due to the growing amount of biological knowledge that is incorporated into metabolic network models, their analysis has become more and more challenging. Here, we examine the capabilities of the recently introduced chemical organization theory (OT) to ease this task. Considering only network stoichiometry, the theory allows the prediction of all potentially persistent species sets and therewith rigorously relates the structure of a network to its potential dynamics. By this, the phenotypes implied by a metabolic network can be predicted without the need for explicit knowledge of the detailed reaction kinetics.

Results

We propose an approach to deal with regulation – and especially inhibitory interactions – in chemical organization theory. One advantage of this approach is that the metabolic network and its regulation are represented in an integrated way as one reaction network. To demonstrate the feasibility of this approach we examine a model by Covert and Palsson (J Biol Chem, 277(31), 2002) of the central metabolism of E. coli that incorporates the regulation of all involved genes. Our method correctly predicts the known growth phenotypes on 16 different substrates. Without specific assumptions, organization theory correctly predicts the lethality of knockout experiments in 101 out of 116 cases. Taking into account the same model specific assumptions as in the regulatory flux balance analysis (rFBA) by Covert and Palsson, the same performance is achieved (106 correctly predicted cases). Two model specific assumptions had to be considered: first, we have to assume that secreted molecules do not influence the regulatory system, and second, that metabolites with increasing concentrations indicate a lethal state.

Conclusion

The introduced approach to model a metabolic network and its regulation in an integrated way as one reaction network makes organization analysis a universal technique to study the potential behavior of biological network models. Applying multiple methods like OT and rFBA is shown to be valuable to uncover critical assumptions and helps to improve model coherence.

Comprehensive transcriptional analysis of the oxidative response in yeast

The oxidative stress response in Saccharomyces cerevisiae has been analyzed by parallel determination of mRNA levels and transcription rates for the entire genome. A mathematical algorithm has been adapted for a dynamic situation such as the response to stress, to calculate theoretical mRNA decay rates from the experimental data. Yeast genes have been grouped into 25 clusters according to mRNA level and transcription rate kinetics, and average mRNA decay rates have been calculated for each cluster. In most of the genes, changes in one or both experimentally determined parameters occur during the stress response. 24% of the genes are transcriptionally induced without an increase in mRNA levels. The lack of parallelism between the evolution of the mRNA amount and transcription rate predicts changes in mRNA stability during stress. Genes for ribosomal proteins and rRNA processing enzymes are abundant among those whose mRNAs are predicted to destabilize. The number of genes whose mRNAs are predicted to stabilize is lower, although some protein folding or proteasomal genes are among the latter. We have confirmed the mathematical predictions for several genes pertaining to different clusters by experimentally determining mRNA decay rates using the regulatable tetO promoter in transcriptional expression conditions not affected by the oxidative stress. This study indicates that the oxidative stress response in yeast cells is not only conditioned by gene transcription but also by the mRNA decay dynamics and that this complex response may be particularly relevant to explain the temporary down-regulation of protein synthesis occurring during stress.

Sharp developmental thresholds defined through bistability by antagonistic gradients of retinoic acid and FGF signaling

The establishment of thresholds along morphogen gradients in the embryo is poorly understood. Using mathematical modeling, we show that mutually inhibitory gradients can generate and position sharp morphogen thresholds in the embryonic space. Taking vertebrate segmentation as a paradigm, we demonstrate that the antagonistic gradients of retinoic acid (RA) and Fibroblast Growth Factor (FGF) along the presomitic mesoderm (PSM) may lead to the coexistence of two stable steady states. Here, we propose that this bistability is associated with abrupt switches in the levels of FGF and RA signaling, which permit the synchronized activation of segmentation genes, such as mesp2, in successive cohorts of PSM cells in response to the segmentation clock, thereby defining the future segments. Bistability resulting from mutual inhibition of RA and FGF provides a molecular mechanism for the all-or-none transitions assumed in the "clock and wavefront" somitogenesis model. Given that mutually antagonistic signaling gradients are common in development, such bistable switches could represent an important principle underlying embryonic patterning.

Genomics and gene transcription kinetics in yeast

As an adaptive response to new conditions, mRNA concentrations in eukaryotes are readjusted after any environmental change. Although mRNA concentrations can be modified by altering synthesis and/or degradation rates, the rapidity of the transition to a new concentration depends on the regulation of mRNA stability. There are several plausible transcriptional strategies following environmental change, reflecting different degrees of compromise between speed of response and cost of synthesis. The recent development of genomic techniques now enables researchers to determine simultaneously (either directly or indirectly) the transcription rates and mRNA half-lifes, together with mRNA concentrations, corresponding to all yeast genes. Such experiments could provide a new picture of the transcriptional response, by enabling us to characterize the kinetic strategies that are used by different genes under given environmental conditions.

Radiation-induced posttranscriptional control of M6P/IGF2r expression in breast cancer cell lines

The mannose 6-phosphate/insulin-like growth factor 2 receptor (M6P/IGF2r), a member of the IGF axis of growth factors, is a negative regulator of cell growth and a putative tumor suppressor gene. Regulation of M6P/IGF2r levels is critical in breast physiology; low expression is associated with various aspects of breast cancer. We have found that ionizing radiation induces the rapid expression of M6P/IGF2r in a dose-dependent manner in MCF7 human breast cancer cells. We show that this increase is mediated, at least in part, by a stabilization of M6P/IGF2r transcripts by radiation in both ER positive (MCF7 and T47D) and ER negative (MDA-MB-231) breast cancer cell lines. It is probable, therefore, that posttranscriptional dysregulation of M6P/IGF2r is a contributing mechanism in breast cancer development and breast cancer response to therapy. This is a novel find that underscores the importance of posttranscriptional control of radiation-induced gene expression-a phenomenon that has often been paradigmatically attributed to transcriptional control.

A Dose-Dependent Decrease in the Fraction of Cases Harboring M6P/IGF2R Mutations in Hepatocellular Carcinomas from the Atomic Bomb Survivors

The risk for hepatocellular carcinoma (HCC) development is significantly heightened in the atomic bomb survivors, but the mechanism is unclear. We have previously reported finding a radiation dose-dependent increase in HCCs with TP53 mutations from the survivors. We now show that, in the same HCC samples, the frequency of 3'-untranslated region (3'UTR) mutations in M6P/IGF2R, a candidate HCC tumor suppressor gene, decreases with dose (P = 0.0091), implying a radiation dose-dependent negative selection of cells harboring such mutations. The fact that they were in the 3'UTR implicates changes in transcript stability rather than in protein function as the mechanism. Moreover, these M6P/IGF2R 3'UTR mutations and the TP53 mutations detected previously were mutually exclusive in most of the tumors, suggesting two independent pathways to HCC development, with the TP53 pathway being more favored with increasing radiation dose than the M6P/IGF2R pathway. These results suggest that tumors attributable to radiation may be genotypically different from tumors of other etiologies and hence may provide a way of distinguishing radiation-induced cancers from "background" cancers--a shift from the current paradigm.

In vivo detection of Egr2 binding to target genes during peripheral nerve myelination

Egr2/Krox20 is a zinc finger transactivator that regulates a diverse array of genes required for peripheral nerve myelination. Although several studies have elucidated the Egr2-regulated gene network, it is not clear if Egr2 regulates its target genes directly or indirectly through induction of other transactivators. Moreover, very few Egr2 binding sites have been identified in regulatory elements of myelin genes. To address this issue, we have successfully adapted chromatin immunoprecipitation assays to test if Egr2 binds directly to target genes in myelinating rat sciatic nerve. These experiments demonstrate direct binding of Egr2 to previously described binding sites within the Schwann cell enhancer of the myelin basic protein gene. Furthermore, we show Egr2 binding to a conserved site within the myelin-associated glycoprotein gene. Finally, our experiments provide the first evidence that Egr2 directly regulates expression of desert hedgehog, which is critically involved in development, maintenance and regeneration of multiple nerve elements including myelinated fibers. Surprisingly, this analysis has identified an apparent preponderance of Egr2 binding sites within conserved intron sequences of several myelin genes. Application of chromatin immunoprecipitation analysis to myelination in vivo will prove to be a valuable asset in assaying transcription factor binding and chromatin modifications during activation of myelin genes.

Endothelin-1 expression is strongly repressed by AU-rich elements in the 3'-untranslated region of the gene

The regulation of the synthesis of the endothelial-derived vasoconstrictor ET-1 (endothelin-1) is a complex process that occurs mainly at the mRNA level. Transcription of the gene accounts for an important part of the regulation of expression, as already described for different modulators such as the cytokine TGF-beta (transforming growth factor-beta). However, very little is known about mechanisms governing ET-1 expression at the post-transcriptional level. The aim of the present study was to investigate the regulation of the ET-1 expression at this level. Since the 3'-UTR (3'-untranslated region) of mRNAs commonly contains genetic determinants for the post-transcriptional control of gene expression, we focused on the potential role of the 3'-UTR of ET-1 mRNA. Experiments performed with luciferase reporter constructs containing the 3'-UTR showed that this region exerts a potent destabilizing effect. Deletional analyses allowed us to locate this activity within a region at positions 924-1127. Some (but not all) of the AREs (AU-rich elements) present in this region were found to be essential for this mRNA-destabilizing activity. We also present evidence that cytosolic proteins from endothelial cells interact specifically with these RNA elements, and that a close correlation exists between the ability of the AREs to destabilize ET-1 mRNA and the binding of proteins to these elements. Our results are compatible with the existence of a strong repressional control of ET-1 expression mediated by destabilization of the mRNA exerted through the interaction of specific cytosolic proteins with AREs present in the 3'-UTR of the gene.

Developing a biosensor for estrogens in water samples: Study of the real-time response of live cells of the estrogen-sensitive yeast strain RMY/ER-ERE using fluorescence microscopy

Using a fluorescein di-beta-D-galactopyranoside (FDG) substrate we show that in live cells of an estrogen-sensitive yeast strain RMY/ER-ERE with human estrogen receptor (ERalpha) gene and the lacZ gene which encodes beta-galactosidase, the uptake of 17beta-estradiol (E2) and the subsequent production of beta-galactosidase enzyme occur quite rapidly, with maximal enzyme-catalyzed product formation evident after about 30 min of exposure to E2. This finding which agrees with the well-known rates of enzyme-catalyzed reactions could have implications for shortening the duration of environmental sample screening and monitoring regimes using yeast-based estrogen assays, and the development of biosensors for environmental estrogens to complement quantification methods.

A thermodynamic model of transcriptome formation

The genome supplies information on both the quality and quantity of the transcriptome. However, as it remains unknown how a cell determines transcript levels from the genome sequences, despite comprehensive knowledge of the cellular components involved, the quantity information held by the genome cannot as yet be derived from nucleotide sequences. The model presented here explains on a thermodynamic basis how the components decode the genome to form and maintain the transcriptome. The model describes the level of a transcript as a pseudo-equilibrium between velocities of synthesis and degradation, both of which are controlled by sequence-specific interactions between protein factors and nucleic acids. Each of the transcript levels can be described by a single equation expressing a function of the activity concentrations of the protein factors. Quantitative information in the genome can thus be transformed into constants determined from the nucleotide sequences. Using this model, the transcriptome can be traced back to the protein factors and the state of chromosome packaging. The total description of transcript levels allows the model to be verified through comparison of derived hypotheses with comprehensive measurements of the transcriptome. The hypotheses thus derived in the present study are well supported by experimental microarray data, confirming the appropriateness of the model.

Steroid hormones regulate gene expression posttranscriptionally by altering the stabilities of messenger RNAs

Hormones exert powerful effects on reproductive physiology by regulating gene expression. Recent discoveries in hormone action emphasize that regulation of gene expression is not restricted to their alterations of the rate of gene transcription. On the contrary, hormonal effects on the stability of a specific mRNA can profoundly alter its steady-state concentration. The mRNAs encoding hormone receptors are commonly regulated by their own hormones to create autoregulatory feedback loops. Negative and positive autoregulatory feedback loops serve to limit or augment hormonal responses, respectively. After introducing the topics of mRNA degradation and regulated stability, this review focuses on steroid hormone effects on mRNA stabilities. Autoregulation of the mRNAs encoding estrogen, progesterone, androgen, and glucocorticoid receptors by the steroid hormones in reproductive tissues is discussed. In addition, steroid hormone effects on the stabilities of many other mRNAs that are important to reproductive biology are reviewed. These include mRNAs that encode gonadotropin hormones, integrins, growth factors, and inflammatory response proteins. Through these posttranscriptional effects, steroid hormones impact the expression of a large population of genes. Studies of the molecular mechanisms of hormonally regulated mRNA stabilities continue to identify critical mRNA sequence elements and their interactions with proteins. Increased understanding of how hormones affect mRNA stability may yield novel approaches to the therapeutic control of hormone effects, including those essential to reproductive physiology in animals.

Continuous-time identification of gene expression models

One objective of systems biology is to create predictive, quantitative models of the transcriptional regulation networks that govern numerous cellular processes. Gene expression measurements, as provided by microarrays, are commonly used in studies that attempt to infer the regulation underlying these processes. At present, most gene expression models that have been derived from microarray data are based in discrete-time, which have limited applicability to common biological data sets, and may impede the integration of gene expression models with other models of biological processes that are formulated as ordinary differential equations (ODEs). To overcome these difficulties, a continuous-time approach for process identification to identify gene expression models based in ODEs was developed. The approach utilizes the modulating functions method of parameter identification. The method was applied to three simulated systems: (1) a linear gene expression model, (2) an autoregulatory gene expression model, and (3) simulated microarray data from a nonlinear transcriptional network. In general, the approach was well suited for identifying models of gene expression dynamics, capable of accurately identifying parameters for small numbers of data samples in the presence of modest experimental noise. Additionally, numerous insights about gene expression modeling were revealed by the case studies.

Modeling feedback loops of the Mammalian circadian oscillator

The suprachiasmatic nucleus governs daily variations of physiology and behavior in mammals. Within single neurons, interlocked transcriptional/translational feedback loops generate circadian rhythms on the molecular level. We present a mathematical model that reflects the essential features of the mammalian circadian oscillator to characterize the differential roles of negative and positive feedback loops. The oscillations that are obtained have a 24-h period and are robust toward parameter variations even when the positive feedback is replaced by a constantly expressed activator. This demonstrates the crucial role of the negative feedback for rhythm generation. Moreover, it explains the rhythmic phenotype of Rev-erbalpha-/- mutant mice, where a positive feedback is missing. The interplay of negative and positive feedback reveals a complex dynamics. In particular, the model explains the unexpected rescue of circadian oscillations in Per2Brdm1/Cry2-/- double-mutant mice (Per2Brdm1 single-mutant mice are arrhythmic). Here, a decrease of positive feedback strength associated with mutating the Per2 gene is compensated by the Cry2-/- mutation that simultaneously decreases the negative feedback strength. Finally, this model leads us to a testable prediction of a molecular and behavioral phenotype: circadian oscillations should be rescued when arrhythmic Per2Brdm1 mutant mice are crossed with Rev- erbalpha -/- mutant mice.

Regulation of metabolic genes in human skeletal muscle by short-term exercise and diet manipulation

Changes in dietary macronutrient intake alter muscle and blood substrate availability and are important for regulating gene expression. However, few studies have examined the effects of diet manipulation on gene expression in human skeletal muscle. The aim of this study was to quantify the extent to which altering substrate availability impacts on subsequent mRNA abundance of a subset of carbohydrate (CHO)- and fat-related genes. Seven subjects consumed either a low- (LOW; 0.7 g/kg body mass CHO) or high- (HIGH; 10 g/kg body mass CHO) CHO diet for 48 h after performing an exhaustive exercise bout to deplete muscle glycogen stores. After intervention, resting muscle and blood samples were taken. Muscle was analyzed for the gene abundances of GLUT4, glycogenin, pyruvate dehydrogenase kinase-4 (PDK-4), fatty acid translocase (FAT/CD36), carnitine palmitoyltransferase I (CPT I), hormone-sensitive lipase (HSL), beta-hydroxyacyl-CoA dehydrogenase (beta-HAD), and uncoupling binding protein-3 (UCP3), and blood samples for glucose, insulin, and free fatty acid (FFA) concentrations. Glycogen-depleting exercise and HIGH-CHO resulted in a 300% increase in muscle glycogen content (P < 0.001) relative to the LOW-CHO condition. FFA concentrations were twofold higher after LOW- vs. HIGH-CHO (P < 0.05). The exercise-diet manipulation exerted a significant effect on transcription of all carbohydrate-related genes, with an increase in GLUT4 and glycogenin mRNA abundance and a reduction in PDK-4 transcription after HIGH-CHO (all P < 0.05). FAT/CD36 (P < 0.05) and UCP3 (P < 0.01) gene transcriptions were increased following LOW-CHO. We conclude that 1) there was a rapid capacity for a short-term exercise and diet intervention to exert coordinated changes in the mRNA transcription of metabolic related genes, and 2) genes involved in glucose regulation are increased following a high-carbohydrate diet.

Homocysteine increases the expression of vascular endothelial growth factor by a mechanism involving endoplasmic reticulum stress and transcription factor ATF4

Vascular endothelial growth factor (VEGF) plays a key role in the development and progression of diabetic retinopathy. We previously demonstrated that amino acid deprivation and other inducers of endoplasmic reticulum-stress (ER stress) up-regulate the expression of VEGF in the retinal-pigmented epithelial cell line ARPE-19. Because homocysteine causes ER stress, we hypothesized that VEGF expression is increased by ambient homocysteine. dl-Homocysteine-induced VEGF expression was investigated in confluent ARPE-19 cultures. Northern analysis showed that homocysteine increased steady state VEGF mRNA levels 4.4-fold. Other thiol-containing compounds, including l-homocysteine thiolactone and DTT, induced VEGF expression 7.9- and 8.8-fold. Transcriptional run-on assays and mRNA decay studies demonstrated that the increase in VEGF mRNA levels was caused by increased transcription rather than mRNA stabilization. VEGF mRNA induction paralleled that of the ER-stress gene GRP78. Homocysteine treatment caused transient phosphorylation of eIF2alpha and an increase in ATF4 protein level. Overexpression of a dominant-negative ATF4 abolished the VEGF response to homocysteine treatment and to amino acid deprivation. VEGF mRNA expression by ATF4-/- MEF did not respond to homocysteine treatment and the response was restored with expression of wild-type ATF4. These studies indicate that expression of the pro-angiogenic factor VEGF is increased by homocysteine and other thiol-containing reductive compounds via ATF4-dependent activation of VEGF transcription.

Importance of input perturbations and stochastic gene expression in the reverse engineering of genetic regulatory networks: insights from an identifiability analysis of an in silico network

Gene expression profiles are an increasingly common data source that can yield insights into the functions of cells at a system-wide level. The present work considers the limitations in information content of gene expression data for reverse engineering regulatory networks. An in silico genetic regulatory network was constructed for this purpose. Using the in silico network, a formal identifiability analysis was performed that considered the accuracy with which the parameters in the network could be estimated using gene expression data and prior structural knowledge (which transcription factors regulate which genes) as a function of the input perturbation and stochastic gene expression. The analysis yielded experimentally relevant results. It was observed that, in addition to prior structural knowledge, prior knowledge of kinetic parameters, particularly mRNA degradation rate constants, was necessary for the network to be identifiable. Also, with the exception of cases where the noise due to stochastic gene expression was high, complex perturbations were more favorable for identifying the network than simple ones. Although the results may be specific to the network considered, the present study provides a framework for posing similar questions in other systems.

A Theoretical Model for the Regulation of Sex-lethal, a Gene That Controls Sex Determination and Dosage Compensation in Drosophila melanogaster

Cell fate commitment relies upon making a choice between different developmental pathways and subsequently remembering that choice. Experimental studies have thoroughly investigated this central theme in biology for sex determination. In the somatic cells of Drosophila melanogaster, Sex-lethal (Sxl) is the master regulatory gene that specifies sexual identity. We have developed a theoretical model for the initial sex-specific regulation of Sxl expression. The model is based on the well-documented molecular details of the system and uses a stochastic formulation of transcription. Numerical simulations allow quantitative assessment of the role of different regulatory mechanisms in achieving a robust switch. We establish on a formal basis that the autoregulatory loop involved in the alternative splicing of Sxl primary transcripts generates an all-or-none bistable behavior and constitutes an efficient stabilization and memorization device. The model indicates that production of a small amount of early Sxl proteins leaves the autoregulatory loop in its off state. Numerical simulations of mutant genotypes enable us to reproduce and explain the phenotypic effects of perturbations induced in the dosage of genes whose products participate in the early Sxl promoter activation.

The effects of estradiol and selective estrogen receptor modulators on gene expression and messenger RNA stability in immortalized sheep endometrial stromal cells and human endometrial adenocarcinoma cells

The purpose of this study was to identify an endometrial cell line that maintained the E2 up-regulation of estrogen receptor (ER) mRNA by enhanced message stability and to assess its dependence on ER protein. Estradiol (E2) effects on gene expression were measured in three cell lines: one immortalized from sheep endometrial stroma (ST) and two from human endometrial adenocarcinomas (Ishikawa and ECC-1). E2 up-regulated ER mRNA levels in ST and Ishikawa cells, but down-regulated ER mRNA levels in ECC-1 cells. E2 up-regulated progesterone receptor (PR), glyceraldehyde 3-phosphate dehydrogenase (GAPDH), and transforming growth factor-alpha (TGF-alpha) in both Ishikawa and ECC-1 cells. The selective estrogen receptor modulator ICI 182,780 antagonized the E2-induced up-regulation of ER and/or PR mRNA levels in all three cells, while another, GW 5638, antagonized the up-regulation of PR mRNA in Ishikawa and ECC-1 cells. In mechanistic studies, E2 had no effect on ER mRNA stability in ST cells and it destabilized ER mRNA in ECC-1 cells. Thus, Ishikawa cells appear to be the most physiologically relevant cell line in which to study the up-regulation of ER mRNA levels by enhanced mRNA stability. Its antagonism by ICI 182,780 reveals that ER protein is involved in this E2 response.

A short-term, high-fat diet up-regulates lipid metabolism and gene expression in human skeletal muscle

BACKGROUND

Dietary fatty acids may be important in regulating gene expression. However, little is known about the effect of changes in dietary fatty acids on gene regulation in human skeletal muscle.

OBJECTIVE

The objective was to determine the effect of altered dietary fat intake on the expression of genes encoding proteins necessary for fatty acid transport and beta-oxidation in skeletal muscle.

DESIGN

Fourteen well-trained male cyclists and triathletes with a mean (+/- SE) age of 26.9 +/- 1.7 y, weight of 73.7 +/- 1.7 kg, and peak oxygen uptake of 67.0 +/- 1.3 mL x kg(-1) x min(-1) consumed either a high-fat diet (HFat: > 65% of energy as lipids) or an isoenergetic high-carbohydrate diet (HCho: 70-75% of energy as carbohydrate) for 5 d in a crossover design. On day 1 (baseline) and again after 5 d of dietary intervention, resting muscle and blood samples were taken. Muscle samples were analyzed for gene expression [fatty acid translocase (FAT/CD36), plasma membrane fatty acid binding protein (FABPpm), carnitine palmitoyltransferase I (CPT I), beta-hydroxyacyl-CoA dehydrogenase (beta-HAD), and uncoupling protein 3 (UCP3)] and concentrations of the proteins FAT/CD36 and FABPpm.

RESULTS

The gene expression of FAT/CD36 and beta -HAD and the gene abundance of FAT/CD36 were greater after the HFat than after the HCho diet (P < 0.05). Messenger RNA expression of FABPpm, CPT I, and UCP-3 did not change significantly with either diet.

CONCLUSIONS

A rapid and marked capacity for changes in dietary fatty acid availability to modulate the expression of mRNA-encoding proteins is necessary for fatty acid transport and oxidative metabolism. This finding is evidence of nutrient-gene interactions in human skeletal muscle.

Identification of unstable transcripts in Arabidopsis by cDNA microarray analysis: rapid decay is associated with a group of touch- and specific clock-controlled genes

mRNA degradation provides a powerful means for controlling gene expression during growth, development, and many physiological transitions in plants and other systems. Rates of decay help define the steady state levels to which transcripts accumulate in the cytoplasm and determine the speed with which these levels change in response to the appropriate signals. When fast responses are to be achieved, rapid decay of mRNAs is necessary. Accordingly, genes with unstable transcripts often encode proteins that play important regulatory roles. Although detailed studies have been carried out on individual genes with unstable transcripts, there is limited knowledge regarding their nature and associations from a genomic perspective, or the physiological significance of rapid mRNA turnover in intact organisms. To address these problems, we have applied cDNA microarray analysis to identify and characterize genes with unstable transcripts in Arabidopsis thaliana (AtGUTs). Our studies showed that at least 1% of the 11,521 clones represented on Arabidopsis Functional Genomics Consortium microarrays correspond to transcripts that are rapidly degraded, with estimated half-lives of less than 60 min. AtGUTs encode proteins that are predicted to participate in a broad range of cellular processes, with transcriptional functions being over-represented relative to the whole Arabidopsis genome annotation. Analysis of public microarray expression data for these genes argues that mRNA instability is of high significance during plant responses to mechanical stimulation and is associated with specific genes controlled by the circadian clock.

Phosphoenolpyruvate carboxykinase is induced in growth-arrested hepatoma cells

Phosphoenolpyruvate carboxykinase (PEPCK) mRNA is elevated in H4IIEC3 rat hepatoma cells cultured at high density, suggesting that PEPCK expression and growth arrest may be coordinately regulated. Induction of growth arrest either by contact inhibition (high culture density) or by serum deprivation correlated with significant increases in PEPCK protein and its mRNA. The observation that PEPCK mRNA was induced by contact inhibition in the presence of serum indicates that the effect of high density is independent of insulin or any other serum component. The magnitudes of the changes in PEPCK expression during growth arrest were greatly enhanced in KRC-7 cells, an H4IIEC3 subclone that is much more sensitive to growth arrest than its parental cell line. Restimulation of proliferation in growth-arrested KRC-7 cells, either by addition of serum or insulin to serum-deprived cells or by replating contact-inhibited cells at low density, caused a rapid decrease in PEPCK expression. However, PEPCK mRNA is not always reduced in proliferating cells since treatment of serum-starved cells with epidermal growth factor stimulated entry into the cell cycle but did not affect PEPCK mRNA levels. Finally, dexamethasone induction of PEPCK mRNA was blunted in cells cultured at high density but was unaffected by the presence or absence of serum. Collectively, these data suggest the possibility of cross-talk between the control of PEPCK expression and growth arrest in KRC-7 cells.

A dynamical systems model to simulate the perturbation kinetics of gene expression by antisense oligonucleotides

Antisense oligonucleotides owe their efficacy to an ability to induce RNase H-dependent suppression of RNA translation, for sufficient time to allow physiological proteolysis. The magnitude and time delay preceding the protein nadir concentration determine the extent and timing of maximum antisense oligonucleotide activity. Antisense oligonucleotide degradation underlies reversal of RNA downregulation. The kinetics of protein downregulation is therefore determined by the complex interaction of both ligand chemistry (nuclease stability, affinity and RNase H activation), and gene expression kinetics. Optimization of antisense oligonucleotide efficacy and experimental design requires understanding of these interactions. The kinetics of protein and RNA downregulation have therefore been simulated by analysing a two-compartment kinetic model incorporating RNase H-dependent transcript degradation. The system of nonlinear differential equations describing this model was solved numerically using Runge-Kutte integration. The timecourse solutions corresponding to the four state variables (RNA, protein, antisense/RNA heteroduplex and antisense oligonucleotide), were determined simultaneously. This allowed systematic in silico examination of the consequences of altering variables such as oligonucleotide concentration, affinity, and stability, or the scheduling of multiple transfections on RNA and protein perturbations. By providing a tool for examining antisense oligonucleotide action theoretically, this heuristic model should facilitate both the rational design and interpretation of antisense experiments.

Identification and cDNA cloning of a novel RNA-binding protein that interacts with the cyclic nucleotide-responsive sequence in the Type-1 plasminogen activator inhibitor mRNA

Incubation of HTC rat hepatoma cells with 8-bromo-cAMP results in a 3-fold increase in the rate of degradation of type-1 plasminogen activator inhibitor (PAI-1) mRNA. We have reported previously that the 3'-most 134 nt of the PAI-1 mRNA is able to confer cyclic nucleotide regulation of message stability onto a heterologous transcript. R-EMSA and UV cross-linking experiments have shown that this 134 nt cyclic nucleotide-responsive sequence (CRS) binds HTC cell cytoplasmic proteins ranging in size from 38 to 76 kDa. Mutations in the A-rich region of the CRS both eliminate cyclic nucleotide regulation of mRNA decay and abolish RN-protein complex formation, suggesting that these RNA-binding proteins may be important regulators of mRNA stability. By sequential R-EMSA and SDS-PAGE we have purified a protein from HTC cell polysomes that binds to the PAI-1 CRS. N-terminal sequence analysis and a search of protein data bases revealed identity with two human sequences of unknown function. We have expressed one of these sequences in E. coli and confirmed that the recombinant protein interacts specifically with the PAI-1 CRS. Mutation of the A-rich portion of the PAI-1 CRS reduces binding by the recombinant PAI-1 RNA-binding protein. The amino acid sequence of this protein includes an RGG box and two arginine-rich regions, but does not include other recognizable RNA binding motifs. Detailed analyses of nucleic acid and protein data bases demonstrate that blocks of this sequence are highly conserved in a number of metazoans, including Arabidopsis, Drosophila, birds, and mammals. Thus, we have described a novel RNA-binding protein that identifies a family of proteins with a previously undefined sequence motif. Our results suggest that this protein, PAI-RBP1, may play a role in regulation of mRNA stability.

Degradation of glial glutamate transporter mRNAs is selectively blocked by inhibition of cellular transcription

Recent studies have demonstrated that the expression of the glial glutamate transporters GLT-1 (glutamate transporter 1) and GLAST (glutamate aspartate transporter) is regulated both in vivo and in vitro. For example, co-culturing with neurons, treatment with N:(6), 2'-O:-dibutyryladenosine 3':5'-cyclic monophosphate (dbcAMP), and treatment with epidermal growth factor all increase the steady-state levels of GLT-1 and GLAST protein in astrocyte cultures. These changes in protein expression are correlated with increased mRNA levels. In the present study, the degradation of GLT-1 and GLAST mRNAs was examined in control and dbcAMP-treated astrocyte cultures after inhibiting transcription with actinomycin D. Although one would predict that inhibition of transcription would cause a decrease in GLT-1 and GLAST mRNAs and that this decrease would depend on the rate of mRNA degradation, the levels of GLT-1 and GLAST mRNAs did not decrease even after 24 h of treatment with actinomycin D. Withdrawal of dbcAMP caused the levels of GLT-1 and GLAST mRNAs to fall to basal levels within 24 h, but this degradation was blocked if actinomycin D was added at the time of dbcAMP withdrawal. Importantly, actinomycin D did not block the degradation of c-fos mRNA also induced by dbcAMP in these cultures. Inhibition of translation with cycloheximide did not stabilize GLT-1 but partially attenuated the degradation of GLAST mRNA. Although the mechanism of this effect remains to be defined, these studies suggest that GLT-1 and GLAST mRNAs belong to a select class of inducible mRNAs stabilized by inhibitors of transcription. The possible relevance of these data to astrocyte differentiation is briefly discussed.