“Coarse” stability and bifurcation analysis using time-steppers: A reaction-diffusion example

Evolutionary, pattern forming partial differential equations (PDEs) are often derived as limiting descriptions of microscopic, kinetic theory-based models of molecular processes (e.g., reaction and diffusion). The PDE dynamic behavior can be probed through direct simulation (time integration) or, more systematically, through stability/bifurcation calculations; time-stepper-based approaches, like the Recursive Projection Method [Shroff, G. M. & Keller, H. B. (1993) SIAM J. Numer. Anal. 30, 1099–1120] provide an attractive framework for the latter. We demonstrate an adaptation of this approach that allows for a direct, effective (“coarse”) bifurcation analysis of microscopic, kinetic-based models; this is illustrated through a comparative study of the FitzHugh-Nagumo PDE and of a corresponding Lattice–Boltzmann model.

High-sensitivity array analysis of gene expression for the early detection of disseminated breast tumor cells in peripheral blood

Early detection is an effective means of reducing cancer mortality. Here, we describe a highly sensitive high-throughput screen that can identify panels of markers for the early detection of solid tumor cells disseminated in peripheral blood. The method is a two-step combination of differential display and high-sensitivity cDNA arrays. In a primary screen, differential display identified 170 candidate marker genes differentially expressed between breast tumor cells and normal breast epithelial cells. In a secondary screen, high-sensitivity arrays assessed expression levels of these genes in 48 blood samples, 22 from healthy volunteers and 26 from breast cancer patients. Cluster analysis identified a group of 12 genes that were elevated in the blood of cancer patients. Permutation analysis of individual genes defined five core genes (P ≤ 0.05, permax test). As a group, the 12 genes generally distinguished accurately between healthy volunteers and patients with breast cancer. Mean expression levels of the 12 genes were elevated in 77% (10 of 13) untreated invasive cancer patients, whereas cluster analysis correctly classified volunteers and patients (P = 0.0022, Fisher's exact test). Quantitative real-time PCR confirmed array results and indicated that the sensitivity of the assay (1:2 × 108 transcripts) was sufficient to detect disseminated solid tumor cells in blood. Expression-based blood assays developed with the screening approach described here have the potential to detect and classify solid tumor cells originating from virtually any primary site in the body.

PEA3 sites within the progression elevated gene-3 (PEG-3) promoter and mitogen-activated protein kinase contribute to differential PEG-3 expression in Ha-ras and v-raf oncogene transformed rat embryo cells

Transformation of normal cloned rat embryo fibroblast (CREF) cells with cellular oncogenes results in acquisition of anchorage-independent growth and oncogenic potential in nude mice. These cellular changes correlate with an induction in the expression of a cancer progression-promoting gene, progression elevated gene-3 (PEG-3). To define the mechanism of activation of PEG-3 as a function of transformation by the Ha-ras and v-raf oncogenes, evaluations of the signaling and transcriptional regulation of the ~2.0 kb promoter region of the PEG-3 gene, PEG-Prom, was undertaken. The full-length and various mutated regions of the PEG-Prom were linked to a luciferase reporter construct and tested for promoter activity in CREF and oncogene-transformed CREF cells. An analysis was also performed using CREF cells doubly transformed with Ha-ras and the Ha-ras specific suppressor gene Krev-1, which inhibits the transformed phenotype in vitro. These assays document an association between expression of the transcription regulator PEA3 and PEG-3. The levels of PEA3 and PEG-3 RNA and proteins are elevated in the oncogenically transformed CREF cells, and reduced in transformation and tumorigenic suppressed Ha-ras/Krev-1 doubly transformed CREF cells. Enhanced tumorigenic behavior, PEG-3 promoter function and PEG-3 expression in Ha-ras transformed cells were all dependent upon increased activity within the mitogen-activated protein kinase (MAPK) pathway. Electrophoretic mobility shift assays and DNase I footprinting experiments indicate that PEA3 binds to sites within the PEG-Prom in transformed rodent cells in an area adjacent to the TATA box in a MAPK-dependent fashion. These findings demonstrate an association between Ha-ras and v-raf transformation of CREF cells with elevated PEA3 and PEG-3 expression, and they implicate MAPK signaling via PEA3 as a signaling cascade involved in activation of the PEG-Prom.

Proteomic analysis of the bacterial cell cycle

A global approach was used to analyze protein synthesis and stability during the cell cycle of the bacterium Caulobacter crescentus. Approximately one-fourth (979) of the estimated C. crescentus gene products were detected by two-dimensional gel electrophoresis, 144 of which showed differential cell cycle expression patterns. Eighty-one of these proteins were identified by mass spectrometry and were assigned to a wide variety of functional groups. Pattern analysis revealed that coexpression groups were functionally clustered. A total of 48 proteins were rapidly degraded in the course of one cell cycle. More than half of these unstable proteins were also found to be synthesized in a cell cycle-dependent manner, establishing a strong correlation between rapid protein turnover and the periodicity of the bacterial cell cycle. This is, to our knowledge, the first evidence for a global role of proteolysis in bacterial cell cycle control.

Genome-wide expression analysis reveals dysregulation of myelination-related genes in chronic schizophrenia

Neuropathological and brain imaging studies suggest that schizophrenia may result from neurodevelopmental defects. Cytoarchitectural studies indicate cellular abnormalities suggestive of a disruption in neuronal connectivity in schizophrenia, particularly in the dorsolateral prefrontal cortex. Yet, the molecular mechanisms underlying these findings remain unclear. To identify molecular substrates associated with schizophrenia, DNA microarray analysis was used to assay gene expression levels in postmortem dorsolateral prefrontal cortex of schizophrenic and control patients. Genes determined to have altered expression levels in schizophrenics relative to controls are involved in a number of biological processes, including synaptic plasticity, neuronal development, neurotransmission, and signal transduction. Most notable was the differential expression of myelination-related genes suggesting a disruption in oligodendrocyte function in schizophrenia.

Molecular Analysis of (R)-(+)-Mandelonitrile Lyase Microheterogeneity in Black Cherry1

The flavoprotein (R)-(+)-mandelonitrile lyase (MDL; EC 4.1.2.10), which plays a key role in cyanogenesis in rosaceous stone fruits, occurs in black cherry (Prunus serotina Ehrh.) homogenates as several closely related isoforms. Biochemical and molecular biological methods were used to investigate MDL microheterogeneity and function in this species. Three novel MDL cDNAs of high sequence identity (designated MDL2, MDL4, and MDL5) were isolated. Like MDL1 and MDL3 cDNAs (Z. Hu, J.E. Poulton [1997] Plant Physiol 115: 1359–1369), they had open reading frames that predicted a flavin adenine dinucleotide-binding site, multiple N-glycosylation sites, and an N-terminal signal sequence. The N terminus of an MDL isoform purified from seedlings matched the derived amino acid sequence of the MDL4 cDNA. Genomic sequences corresponding to the MDL1, MDL2, and MDL4 cDNAs were obtained by polymerase chain reaction amplification of genomic DNA. Like the previously reported mdl3 gene, these genes are interrupted at identical positions by three short, conserved introns. Given their overall similarity, we conclude that the genes mdl1, mdl2, mdl3, mdl4, and mdl5 are derived from a common ancestral gene and constitute members of a gene family. Genomic Southern-blot analysis showed that this family has approximately eight members. Northern-blot analysis using gene-specific probes revealed differential expression of the genes mdl1, mdl2, mdl3, mdl4, and mdl5.

Differential Expression of Three Members of the 1-Aminocyclopropane-1-Carboxylate Synthase Gene Family in Carnation1

We investigated the expression patterns of three 1-aminocyclopropane-1-carboxylate (ACC) synthase genes in carnation (Dianthus caryophyllus cv White Sim) under conditions previously shown to induce ethylene biosynthesis. These included treatment of flowers with 2,4-dichlorophenoxyacetic acid, ethylene, LiCl, cycloheximide, and natural and pollination-induced flower senescence. Accumulation of ACC synthase transcripts in leaves following mechanical wounding and treatment with 2,4-dichlorophenoxyacetic acid or LiCl was also determined by RNA gel-blot analysis. As in other species, the carnation ACC synthase genes were found to be differentially regulated in a tissue-specific manner. DCACS2 and DCACS3 were preferentially expressed in styles, whereas DCACS1 mRNA was most abundant in petals. Cycloheximide did not induce increased accumulation of ACC synthase transcripts in carnation flowers, whereas the expression of ACC synthase was up-regulated by auxin, ethylene, LiCl, pollination, and senescence in a floral-organ-specific manner. Expression of the three ACC synthases identified in carnation did not correspond to elevated ethylene biosynthesis from wounded or auxin-treated leaves, and there are likely additional members of the carnation ACC synthase gene family responsible for ACC synthase expression in vegetative tissues.

Differential fluorescence labeling of cysteinyl clusters uncovers high tissue levels of thionein

The isolation of thionein (T) from tissues has not been reported heretofore. T contains 20 cysteinyl residues that react with 7-fluorobenz-2-oxa-1,3-diazole-4-sulfonamide to form fluorescent adducts. In metallothionein (MT) the cysteinyl residues, which are bound to zinc, do not react. However, they do react in the presence of a chelating agent such as EDTA. The resultant difference in chemical reactivity provides a means to measure T in the absence of EDTA, (MT + T) in its presence, and, of course, MT by difference. The 7-fluorobenz-2-oxa-1,3-diazole-4-sulfonamide derivative of T can be isolated from tissue homogenates by HPLC and quantified fluorimetrically with a detection limit in the femtomolar range and a linear response over 3 orders of magnitude. Analysis of liver, kidney, and brain of rats reveals almost as much T as MT. Moreover, in contrast to earlier views, MT in tissue extracts appears to be less stable than T. The existence of T in tissues under normal physiological conditions has important implications for its function both in zinc metabolism and the redox balance of the cell.

Differential Regulation of Enolase during Anaerobiosis in Maize1

It was reported previously that enolase enzyme activity and ENO1 transcript levels are induced by anaerobic stress in maize (Zea mays). Here we show that not all isoforms of maize enolase are anaerobically induced. We cloned and sequenced a second enolase cDNA clone (pENO2) from maize. Sequence analysis showed that pENO2 shares 75.6% nucleotide and 89.5% deduced amino acid sequence identity with pENO1 and is encoded by a distinct gene. Expression of ENO2 is constitutive under aerobic conditions, whereas ENO1 levels are induced 10-fold in maize roots after 24 h of anaerobic treatment. Western-blot analysis and N-terminal sequencing of in vivo-labeled maize roots identified two major proteins selectively synthesized upon anaerobic stress as isozymes of enolase. We describe the expression of enolase in maize roots under anaerobic stress.

Differential Regulation of Sugar-Sensitive Sucrose Synthases by Hypoxia and Anoxia Indicate Complementary Transcriptional and Posttranscriptional Responses1

The goal of this research was to resolve the hypoxic and anoxic responses of maize (Zea mays) sucrose (Suc) synthases known to differ in their sugar regulation. The two maize Suc synthase genes, Sus1 and Sh1, both respond to sugar and O2, and recent work suggests commonalities between these signaling systems. Maize seedlings (NK508 hybrid, W22 inbred, and an isogenic sh1-null mutant) were exposed to anoxic, hypoxic, and aerobic conditions (0, 3, and 21% O2, respectively), when primary roots had reached approximately 5 cm. One-centimeter tips were excised for analysis during the 48-h treatments. At the mRNA level, Sus1 was rapidly up-regulated by hypoxia (approximately 5-fold in 6 h), whereas anoxia had less effect. In contrast, Sh1 mRNA abundance increased strongly under anoxia (approximately 5-fold in 24 h) and was much less affected by hypoxia. At the enzyme level, total Suc synthase activity rose rapidly under hypoxia but showed little significant change during anoxia. The contributions of SUS1 and SH1 activities to these responses were dissected over time by comparing the sh1-null mutant with the isogenic wild type (Sus+, Sh1+). Sh1-dependent activity contributed most markedly to a rapid protein-level response consistently observed in the first 3 h, and, subsequently, to a long-term change mediated at the level of mRNA accumulation at 48 h. A complementary midterm rise in SUS1 activity varied in duration with genetic background. These data highlight the involvement of distinctly different genes and probable signal mechanisms under hypoxia and anoxia, and together with earlier work, show parallel induction of “feast and famine” Suc synthase genes by hypoxia and anoxia, respectively. In addition, complementary modes of transcriptional and posttranscriptional regulation are implicated by these data, and provide a mechanism for sequential contributions from the Sus1 and Sh1 genes during progressive onset of naturally occurring low-O2 events.

4-Coumarate:Coenzyme A Ligase in Hybrid Poplar1 : Properties of Native Enzymes, cDNA Cloning, and Analysis of Recombinant Enzymes

The enzyme 4-coumarate:coenzyme A ligase (4CL) is important in providing activated thioester substrates for phenylpropanoid natural product biosynthesis. We tested different hybrid poplar (Populus trichocarpa × Populus deltoides) tissues for the presence of 4CL isoforms by fast-protein liquid chromatography and detected a minimum of three 4CL isoforms. These isoforms shared similar hydroxycinnamic acid substrate-utilization profiles and were all inactive against sinapic acid, but instability of the native forms precluded extensive further analysis. 4CL cDNA clones were isolated and grouped into two major classes, the predicted amino acid sequences of which were 86% identical. Genomic Southern blots showed that the cDNA classes represent two poplar 4CL genes, and northern blots provided evidence for their differential expression. Recombinant enzymes corresponding to the two genes were expressed using a baculovirus system. The two recombinant proteins had substrate utilization profiles similar to each other and to the native poplar 4CL isoforms (4-coumaric acid > ferulic acid > caffeic acid; there was no conversion of sinapic acid), except that both had relatively high activity toward cinnamic acid. These results are discussed with respect to the role of 4CL in the partitioning of carbon in phenylpropanoid metabolism.

Two Isoforms of NADPH:Cytochrome P450 Reductase in Arabidopsis thaliana : Gene Structure, Heterologous Expression in Insect Cells, and Differential Regulation

We have investigated two NADPH-cytochrome (Cyt) P450 reductase isoforms encoded by separate genes (AR1 and AR2) in Arabidopsis thaliana. We isolated AR1 and AR2 cDNAs using a mung bean (Phaseolus aureus L.) NADPH-Cyt P450 reductase cDNA as a probe. The recombinant AR1 and AR2 proteins produced using a baculovirus expression system showed similar Km values for Cyt c and NADPH, respectively. In the reconstitution system with a recombinant cinnamate 4-hydroxylase (CYP73A5), the recombinant AR1 and AR2 proteins gave the same level of cinnamate 4-hydroxylase activity (about 70 nmol min−1 nmol−1 P450). The AR2 gene expression was transiently induced by 4- and 3-fold within 1 h of wounding and light treatments, respectively, and the induction time course preceded those of CYP73A5 and a phenylalanine ammonia-lyase (PAL1) gene. On the contrary, the AR1 expression level did not change during the treatments. Analysis of the AR1 and AR2 gene structure revealed that only the AR2 promoter contained three putative sequence motifs (boxes P, A, and L), which are involved in the coordinated expression of CYP73A5 and other phenylpropanoid pathway genes. These results suggest the possibility that AR2 transcription may be functionally linked to the induced levels of phenylpropanoid pathway enzymes.

Coordinating DNA replication initiation with cell growth: differential roles for DnaA and SeqA proteins.

We describe here the development of a new approach to the analysis of Escherichia coli replication control. Cells were grown at low growth rates, in which case the bacterial cell cycle approximates that of eukaryotic cells with G1, S, and G2 phases: cell division is followed sequentially by a gap period without DNA replication, replication of the single chromosome, another gap period, and finally the next cell division. Flow cytometry of such slowly growing cells reveals the timing of replication initiation as a function of cell mass. The data show that initiation is normally coupled to cell physiology extremely tightly: the distribution of individual cell masses at the time of initiation in wild-type cells is very narrow, with a coefficient of variation of less than 9%. Furthermore, a comparison between wild-type and seqA mutant cells shows that initiation occurs at a 10-20% lower mass in the seqA mutant, providing direct evidence that SeqA is a bona fide negative regulator of replication initiation. In dnaA (Ts) mutants the opposite is found: the mass at initiation is dramatically increased and the variability in cell mass at initiation is much higher than that for wild-type cells. In contrast to wild-type and dnaA(Ts) cells, seqA mutant cells frequently go through two initiation events per cell division cycle, and all the origins present in each cell are not initiated in synchrony. The implications for the complex interplay amongst growth, cell division, and DNA replication are discussed.

Parallel human genome analysis: microarray-based expression monitoring of 1000 genes.

Microarrays containing 1046 human cDNAs of unknown sequence were printed on glass with high-speed robotics. These 1.0-cm2 DNA "chips" were used to quantitatively monitor differential expression of the cognate human genes using a highly sensitive two-color hybridization assay. Array elements that displayed differential expression patterns under given experimental conditions were characterized by sequencing. The identification of known and novel heat shock and phorbol ester-regulated genes in human T cells demonstrates the sensitivity of the assay. Parallel gene analysis with microarrays provides a rapid and efficient method for large-scale human gene discovery.

Differential functions of the two Src homology 2 domains in protein tyrosine phosphatase SH-PTP1.

SH-PTP1 (also known as PTP1C, HCP, and SHP) is a non-transmembrane protein tyrosine phosphatase (PTPase) containing two tandem Src homology 2 (SH2) domains. We show here that the two SH2 (N-SH2 and C-SH2) domains in SH-PTP1 have different functions in regulation of the PTPase domain and thereby signal transduction. While the N-terminal SH2 domain is both necessary and sufficient for autoinhibition through an intramolecular association with the PTPase domain, truncation of the C-SH2 domain [SH-PTP1 (delta CSH2) construct] has little effect on SH-PTP1 activity. A synthetic phosphotyrosine residue (pY) peptide derived from the erythropoietin receptor (EpoR pY429) binds to the N-SH2 domain and activates both wild-type SH-PTP1 and SH-PTP1 (delta CSH2) 60- to 80-fold. Another pY peptide corresponding to a phosphorylation site on the IgG Fc receptor (Fc gamma RIIB1 pY309) associates with both the C-SH2 domain (Kd = 2.8 microM and the N-SH2 domain (Kd = 15.0 microM) and also activates SH-PTP1 12-fold. By analysis of the effect of the Fc gamma RIIB1 pY309 peptide on SH-PTP1 (delta CSH2), SH-PTP1 (R30K/R33E), SH-PTP1 (R30K/R136K), and SH-PTP1 (R136K) mutants in which the function of either the N- or C-SH2 domain has been impaired, we have determined that both synthetic pY peptides stimulate SH-PTP1 by binding to its N-SH2 domain; binding of pY ligand to the C-SH2 domain has no effect on SH-PTP1 activity. We propose that the N-terminal SH2 domain serves both as a regulatory domain and as a recruiting unit, whereas the C-terminal SH2 domain acts merely as a recruiting unit.