Bifurcation and singularity analysis of a molecular network for the induction of long-term memory.

Withdrawal reflexes of the mollusk Aplysia exhibit sensitization, a simple form of long-term memory (LTM). Sensitization is due, in part, to long-term facilitation (LTF) of sensorimotor neuron synapses. LTF is induced by the modulatory actions of serotonin (5-HT). Pettigrew et al. developed a computational model of the nonlinear intracellular signaling and gene network that underlies the induction of 5-HT-induced LTF. The model simulated empirical observations that repeated applications of 5-HT induce persistent activation of protein kinase A (PKA) and that this persistent activation requires a suprathreshold exposure of 5-HT. This study extends the analysis of the Pettigrew model by applying bifurcation analysis, singularity theory, and numerical simulation. Using singularity theory, classification diagrams of parameter space were constructed, identifying regions with qualitatively different steady-state behaviors. The graphical representation of these regions illustrates the robustness of these regions to changes in model parameters. Because persistent protein kinase A (PKA) activity correlates with Aplysia LTM, the analysis focuses on a positive feedback loop in the model that tends to maintain PKA activity. In this loop, PKA phosphorylates a transcription factor (TF-1), thereby increasing the expression of an ubiquitin hydrolase (Ap-Uch). Ap-Uch then acts to increase PKA activity, closing the loop. This positive feedback loop manifests multiple, coexisting steady states, or multiplicity, which provides a mechanism for a bistable switch in PKA activity. After the removal of 5-HT, the PKA activity either returns to its basal level (reversible switch) or remains at a high level (irreversible switch). Such an irreversible switch might be a mechanism that contributes to the persistence of LTM. The classification diagrams also identify parameters and processes that might be manipulated, perhaps pharmacologically, to enhance the induction of memory. Rational drug design, to affect complex processes such as memory formation, can benefit from this type of analysis.

Analysis of the mechanisms that maintain coordinate production of $\alpha$- and $\beta$-tubulin in mammalian cells

Alpha and beta tubulin are essential proteins in all eukaryotic cells. To study how cells maintain coordinate levels of these two interacting proteins, we have used PCR to add a 9 amino acid epitope from influenza hemagglutinin protein onto the carboxyl terminus of $\alpha$1 and $\beta$1-tubulin. The chimeric tubulin genes (HA$\alpha$1 and HA$\beta$1) were transfected into CHO cells and cell lines that stably express each gene were selected. Cells transfected with HA-tubulin do not exhibit any gross changes in growth or morphology. Immunofluorescence analysis demonstrated that HA-tubulins incorporate into both cytoplasmic and spindle microtubules. A quantitative biochemical assay was used to show that HA-tubulins incorporate into microtubules to a normal extent and do not alter the steady state distribution of endogenous tubulin between monomer and polymer pools. Two-dimensional gel analysis of pulse-labeled cells indicated that when HA$\beta$1-tubulin is expressed at high levels, it slightly represses the synthesis of the endogenous $\beta$-tubulin but produces a small increase in the synthesis of $\alpha$-tubulin. Analysis of cells labeled to steady state showed that HA$\beta$1-tubulin accumulates to a similar level as the wild-type gene product, but together these polypeptides produce only a small increase in total tubulin content consistent with the increased synthesis of $\alpha$-tubulin. It thus appears that HA$\beta$1-tubulin successfully competes with endogenous $\beta$-tubulin for heterodimer formation and that free $\beta$-tubulin subunits (endogenous and HA$\beta$1) are selectively degraded to maintain coordinate amounts of $\alpha$- and $\beta$-tubulin. In addition, the increased synthesis of $\alpha$-tubulin suggested the existence of a mechanism to ensure coordinate synthesis of $\alpha$- and $\beta$-tubulin subunits. To analyze whether reciprocal changes in endogenous tubulin synthesis occur when $\alpha$-tubulin is overexpressed, stably transfected CHO cell lines were isolated in which HA$\alpha$1-tubulin represents 50% of the total $\alpha$-tubulin, and its relative abundance can be further increased to 85-90% by treatment with sodium butyrate. In contrast with results obtained using HA$\beta$1-tubulin, transfection of HA$\alpha$1-tubulin decreased the synthesis of endogenous $\alpha$-tubulin to 60% of normal with little or no change in $\beta$-tubulin synthesis. When the transfected cells were treated with sodium butyrate to further increase HA$\beta$1-tubulin production, a larger decrease in the synthesis of endogenous $\alpha$-tubulin (to 30% of normal) was observed. The repression on the synthesis of endogenous $\alpha$-tubulin polypeptide was found to be directly proportional to the expression of HA$\alpha$1-tubulin indicating the existence of an autoregulatory loop, where $\alpha$-tubulin inhibits its own synthesis. To determine whether overproduction of HA$\alpha$1-tubulin affected the transcription, message stability or translation of endogenous $\alpha$-tubulin, the steady state levels of $\alpha$-tubulin mRNA were analyzed by ribonuclease protection assays. The results showed that the steady state level of $\alpha$-tubulin mRNA is not affected by the overexpression of HA$\alpha$1-tubulin, indicating that the repression is translational. The results are compatible with a model in which $\beta$-tubulin synthesis is largely unperturbed by overexpression of other tubulin subunits, and excess $\beta$-tubulin subunits are rapidly degraded to maintain coordinate $\alpha$- and $\beta$-tubulin levels at steady state. In contrast, free $\alpha$-tubulin represses its own synthesis at the translational level, suggesting that its level of production may be controlled by the amount of $\beta$-tubulin available for heterodimer formation. ^

Functional analysis of GNA -1, a Galphai superfamily member found in the filamentous fungus Neurospora crassa

Heterotrimeric GTP-binding proteins, G proteins, are integral components of eukaryotic signaling systems linking extracellular signals to intracellular responses. Through coupling to seven-transmembrane helix receptors, G proteins convey primary signaling events into multi-leveled cascades of intracellular activity by regulating downstream enzymes, collectively called effectors. The effector enzymes regulated by G proteins include adenylyl cyclase, cAMP phosphodiesterase, phospolipase C-β, mitogen-activated protein kinases, and ion channels. ^ Neurospora crassa is a multicellular, filamentous fungus that is capable of both asexual and sexual reproduction by elaboration of specialized, developmentally controlled structures that give rise to either asexual or sexual spores, respectively. N. crassa possesses at least three heterotrimeric Gα proteins (GNA-1–3) and one Gβ subunit (GNB-1). GNA-1 was the first microbial protein that could be classified in the Gαi superfamily based on its amino acid identity and demonstration that it is a substrate for ADP-ribosylation by pertussis toxin. ^ Experiments were designed to identify the signal transduction pathways and the effector enzymes regulated by GNA-1. Targeted gene-replacement of gna-1 revealed that GNA-1 controls multiple developmental pathways including both asexual and sexual reproduction, maintenance of growth, and resistance to osmotic stress. The Gαi and Gαz members of the Gαi superfamily negatively regulate adenylyl cyclase activity in mammalian cells; therefore, adenylyl cyclase and cAMP levels were measured in Δgna-1 strains and also in strains that were deleted for both gna-1 and gna-2, a second Gα in N. crassa shown to have overlapping functions with GNA-1. Direct measurements of adenylyl cyclase activity revealed that GNA-1, but not GNA-2, was responsible for GTP-stimulated adenylyl cyclase activity in N. crassa. Furthermore, anti-GNA-1 IgG could specifically inhibit GTP-stimulated adenylyl cyclase activity in wild-type strain extracts. These studies also provided evidence that N. crassa possesses feedback mechanisms that control steady-state cAMP levels through indirect regulation of cAMP-phosphodiesterase activity; mutations in gna-1 and gna-2 were additive in their effect on lowering cAMP-phosphodiesterase activity under growth conditions where steady-state cAMP levels were normal but GTP-stimulated adenylyl cyclase activity was reduced 90% in comparison to control strains. ^ Genetic and biochemical epistasis experiments utilizing a Δ gna-1 cr-1 mutant suggest that GNA-1 is essential for female fertility in a cAMP-independent pathway. Furthermore, deletion of gna-1 in a cr-1 background exacerbated many of the defects already observed in the cr-1 strain including more severe growth restriction and developmental defects. However, deletion of gna-1 had no effect on the increased thermotolerance of cr-1, which has been attributed to loss of cAMP. cr-1 possesses GNA-1 protein, and crude membrane fractions from this strain reconstituted GTP-stimulated adenylyl cyclase activity in Δgna-1 membrane fractions. These studies provide direct evidence for the involvement of Gα proteins in the regulation of adenylyl cyclase activity in eukaryotic microbes. ^

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.

Differential dynamic properties of scleroderma fibroblasts in response to perturbation of environmental stimuli.

Diseases are believed to arise from dysregulation of biological systems (pathways) perturbed by environmental triggers. Biological systems as a whole are not just the sum of their components, rather ever-changing, complex and dynamic systems over time in response to internal and external perturbation. In the past, biologists have mainly focused on studying either functions of isolated genes or steady-states of small biological pathways. However, it is systems dynamics that play an essential role in giving rise to cellular function/dysfunction which cause diseases, such as growth, differentiation, division and apoptosis. Biological phenomena of the entire organism are not only determined by steady-state characteristics of the biological systems, but also by intrinsic dynamic properties of biological systems, including stability, transient-response, and controllability, which determine how the systems maintain their functions and performance under a broad range of random internal and external perturbations. As a proof of principle, we examine signal transduction pathways and genetic regulatory pathways as biological systems. We employ widely used state-space equations in systems science to model biological systems, and use expectation-maximization (EM) algorithms and Kalman filter to estimate the parameters in the models. We apply the developed state-space models to human fibroblasts obtained from the autoimmune fibrosing disease, scleroderma, and then perform dynamic analysis of partial TGF-beta pathway in both normal and scleroderma fibroblasts stimulated by silica. We find that TGF-beta pathway under perturbation of silica shows significant differences in dynamic properties between normal and scleroderma fibroblasts. Our findings may open a new avenue in exploring the functions of cells and mechanism operative in disease development.

Differential dynamic properties of scleroderma fibroblasts in response to perturbation of environmental stimuli.

Diseases are believed to arise from dysregulation of biological systems (pathways) perturbed by environmental triggers. Biological systems as a whole are not just the sum of their components, rather ever-changing, complex and dynamic systems over time in response to internal and external perturbation. In the past, biologists have mainly focused on studying either functions of isolated genes or steady-states of small biological pathways. However, it is systems dynamics that play an essential role in giving rise to cellular function/dysfunction which cause diseases, such as growth, differentiation, division and apoptosis. Biological phenomena of the entire organism are not only determined by steady-state characteristics of the biological systems, but also by intrinsic dynamic properties of biological systems, including stability, transient-response, and controllability, which determine how the systems maintain their functions and performance under a broad range of random internal and external perturbations. As a proof of principle, we examine signal transduction pathways and genetic regulatory pathways as biological systems. We employ widely used state-space equations in systems science to model biological systems, and use expectation-maximization (EM) algorithms and Kalman filter to estimate the parameters in the models. We apply the developed state-space models to human fibroblasts obtained from the autoimmune fibrosing disease, scleroderma, and then perform dynamic analysis of partial TGF-beta pathway in both normal and scleroderma fibroblasts stimulated by silica. We find that TGF-beta pathway under perturbation of silica shows significant differences in dynamic properties between normal and scleroderma fibroblasts. Our findings may open a new avenue in exploring the functions of cells and mechanism operative in disease development.

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.

Patient Characteristics and Outcomes of Gastrointestinal Stromal Tumor Patients Undergoing Imatinib Plasma Level Testing

Although gastrointestinal stromal tumor (GIST) is effectively treated with imatinib, there are a number of clinical challenges in the optimal treatment of these patients. The plasma steady-state trough level of imatinib has been proposed to correlate with clinical outcome. Plasma imatinib level may be affected by a number of patient characteristics. Additionally, the ideal plasma trough concentration of imatinib is likely to vary based on the KIT genotype (genotype determines imatinib binding affinity) of the individual patient. Patients’ genotype or plasma imatinib level may influence the type and duration of response that is appreciable by clinical evaluation. The objectives of this study were to determine effects of genotype on the type of response appreciable by current imaging criteria, to determine the distribution of plasma imatinib levels in patients with GIST, to determine factors that correlate with plasma imatinib level, to determine the incremental effects of imatinib dose escalation; and to explore the median plasma levels and outcomes of patients with various KIT mutations. We therefore obtained KIT mutation information and analyzed CT response for size and density measurement of GISTs at baseline and within the first four moths of imatinib treatment. In 126 patients with metastatic/unresectable disease, the KIT genotype of patients’ tumor was significantly associated with unique response characteristics measurable by CT. Furthermore, hepatic and peritoneal metastases differed in their response characteristics. A subgroup of patients with KIT exon 9 mutation, who received higher doses of imatinib and experienced higher trough imatinib levels, experienced improved progression-free survival similar to that of KIT exon 11 patients. Therefore, we have found that imatinib plasma levels were higher in patients with elevated Aspartate amino transferase, were women, were older, or were being treated concomitantly with CYP450 substrate drugs. As expected, CYP450 inducers correlated with a lower plasma imatinib levels in GIST patients. Renal metabolism of imatinib accounts for <10%, so it was not included in the analysis but may affect covariates. Interestingly, there was a trend for low imatinib levels and inferior progression-free survival in patients who had undergone complete gastrectomy. Patients with KIT exon 9 mutation in our cohort received higher imatinib doses, experienced higher trough imatinib levels, and experienced a PFS similar to that of KIT exon 11 patients. In conclusion, imatinib plasma levels are influenced by a number of patient characteristics. The optimal imatinib plasma level for individual patients is not known but is an area of intense investigation. Our study confirms patients with KIT exon 9 mutations benefit from high-dose imatinib and higher trough imatinib levels.

Insight into the strand exchange reaction promoted by the RecA protein of {\it Escherichia coli\/}

In vitro, RecA protein catalyses the exchange of single strands of DNA between different DNA molecules with sequence complementarity. In order to gain insight into this complex reaction and the roles of ATP binding and hydrolysis, two different approaches have been taken. The first is to use short single-stranded deoxyoligonucleotides as the ssDNA in strand exchange. These were used to determine the signal for hydrolysis and the structure of the RecA-DNA complex that hydrolyses ATP. I present a defined kinetic analysis of the nucleotide triphosphatase activity of RecA protein using short oligonucleotides as ssDNA cofactor. I compare the effects of both homopolymers and mixed base composition oligomers on the ATPase activity of RecA protein. I examine the steady state kinetic parameters of the ATPase reaction using these oligonucleotides as ssDNA cofactor, and show that although RecA can both bind to, and utilise, oligonucleotides 7 to 20 residues in length to support the repressor cleavage activity of RecA, these oligonucleotides are unable to efficiently stimulate the ATPase activity of RecA protein. I show that the K$\sb{\rm m}\sp{\rm ATP}$, the Hill coefficient for ATP binding, the extent of reaction, and k$\sb{\rm cat}$ are all a function of ssDNA chain length and that secondary structure may also play a role in determining the effects of a particular chain length on the ATPase activity of RecA protein.^ The second approach is to utilise one of the many mutants of RecA to gain insight into this complex reaction. The mutant selected was RecA1332. Surprisingly, in vitro, this mutant possesses a DNA-dependent ATPase activity. The K$\sb{\rm m}\sp{\rm ATP}$, Hill coefficient for ATP binding, and K$\sb{\rm m}\sp{\rm DNA}$ are similar to that of wild type. k$\sb{\rm cat}$ for the ATPase activity is reduced 3 to 12-fold, however. RecA1332 is unable to use deoxyoligonucleotides as DNA cofactors in the ATPase reaction, and demonstrates an increased sensitivity to inhibition by monovalent ions. It is able to perform strand exchange with ATP and ATP$\lbrack\gamma\rbrack$S but not with UTP, whereas the wild type protein is able to use all three nucleotide triphosphates. RecA1332 appears to be slowed in its ability to form intermediates and to convert these intermediates to products. (Abstract shortened by UMI.) ^

The unstable (CTG)$\sb{\rm n}$ trinucleotide repeat associated with myotonic dystrophy: What it does and where it came from

Myotonic dystrophy (DM), an autosomal dominant disorder mapping to human chromosome 19q13.3, is the most common neuromuscular disease in human adults.^ Following the identification of the mutation underlying the DM phenotype, an unstable (CTG)$\sb{n}$ trinucleotide repeat in the 3$\prime$ untranslated region (UTR) of a gene encoding a ser/thr protein kinase named DM protein kinase (DMPK), the study was targeted at two questions: (1) the identification of the disease-causing mechanism(s) of the unstable repeat, and at a more basic level, (2) the identification of the origin and the mechanism(s) involved in repeat instability. The first goal was to identify the pathophysiological mechanisms of the (CTG)$\sb{n}$ repeat.^ The normal repeat is transcribed but not translated; therefore, initial studies centered on the effect on RNA transcript levels. The vast majority of DM affecteds are heterozygous for the mutant expansion, so that the normal allele interferes with the analysis of the mutant allele. A quantitative allele-specific RT-PCR procedure was developed and applied to a spectrum of patient tissue samples and cell lines. Equal levels of unprocessed pre-mRNA were determined for the wild type (+) and disease (DM) alleles in skeletal muscle and cell lines of heterozygous DM patients, indicating that any nucleosome binding has no effect at the level of transcriptional initiation and transcription of the mutant DMPK locus. In contrast, processed mRNA levels from the DM allele were reduced relative to the + allele as the size of the expansion increased. The unstable repeat, therefore, impairs post-transcriptional processing of DM allele transcripts. This phenomenon has profound effects on overall DMPK locus steady-state transcript levels in cells missing a wild type allele and does not appear to be mediated by imprinting, decreased mRNA stability, generation of aberrant splice forms, or absence of polyadenylation of the mutant allele.^ In Caucasian DM subjects, the unstable repeat is in complete linkage disequlibrium with a single haplotype composed of nine alleles within and flanking DMPK over a physical distance of 30 kb. A detailed haplotype analysis of the DM region was conducted on a Nigerian (Yoruba) DM family, the only indigenous sub-Saharan DM case reported to date. Each affected member of this family had an expanded (CTG)$\sb{n}$ repeat in one of their DMPK alleles. However, unlike all other DM populations studied thus far, disassociation of the (CTG)$\sb{n}$ repeat expansion from other alleles of the putative predisposing haplotype was found. Thus, the expanded (CTG)$\sb{n}$ repeat in this family was the result of an independent mutational event. Consequently, the origin of DM is unlikely the result of a single mutational event, and the hypothesis that a single ancestral haplotype predisposes to repeat expansion is not compelling. (Abstract shortened by UMI.) ^

Regulatory factors and control of anthrax toxin gene expression

Bacillus anthracis plasmid pXO1 carries genes for three anthrax toxin proteins, pag (protective antigen), cya (edema factor), and lef (lethal factor). Expression of the toxin genes is enhanced by two signals: CO$\sb2$/bicarbonate and temperature. The CO$\sb2$/bicarbonate effect requires the presence of pXO1. I hypothesized that pXO1 harbors a trans-acting regulatory gene(s) required for CO$\sb2$/bicarbonate-enhanced expression of the toxin genes. Characterization of such a gene(s) will lead to increased understanding of the mechanisms by which B. anthracis senses and responds to host environments.^ A regulatory gene (atxA) on pXO1 was identified. Transcription of all three toxin genes is decreased in an atxA-null mutant. There are two transcriptional start sites for pag. Transcription from the major site, P1, is enhanced in elevated CO$\sb2$. Only P1 transcripts are significantly decreased in the atxA mutant. Deletion analysis of the pag upstream region indicates that the 111-bp region upstream of the P1 site is sufficient for atxA-mediated increase of this transcript. The cya and lef genes each have one apparent transcriptional start site. The cya and lef transcripts are significantly decreased in the atxA mutant. The atxA mutant is avirulent in mice. The antibody response to all three toxin proteins is significantly decreased in atxA mutant-infected mice. These data suggest that the atxA gene product activates expression of the toxin genes and is essential for virulence.^ Since expression of the toxin genes is dependent on atxA, whether increased toxin gene expression in response to CO$\sb2$/bicarbonate and temperature is associated with increased atxA expression was investigated. I monitored steady state levels of atxA mRNA and AtxA protein in different growth conditions. The results indicate that expression of atxA is not influenced by CO$\sb2$/bicarbonate. Steady state levels of atxA mRNA and AtxA protein are higher at 37$\sp\circ$C than 28$\sp\circ$C. However, increased pag expression at high temperature can not be attributed directly to increased atxA expression.^ There is evidence that an additional factor(s) may be involved in regulation of pag. Expression of pag in strains overproducing AtxA is significantly decreased compared to the wildtype strain. A specific interaction of tagged-AtxA with the pag upstream DNA has not been demonstrated. Furthermore, four proteins in B. anthracis extract can be co-immunoprecipitated with tagged-AtxA. Amino-terminal sequence of one protein has been determined and found highly homologous to chaperonins of GroEL family. Studies are under way to determine if this GroEL-like protein interactions with AtxA and plays any role in atxA-mediated activation of toxin genes. ^

Characterization of osteopontin charge forms produced by osteoblasts

Osteopontin (OPN) is a highly-phosphorylated extracellular matrix protein localized in bone, kidney, placenta, T-lymphocytes, macrophages, smooth muscle of the vascular system, milk, urine, and plasma. In ROS 17/2.8 osteoblast-like osteosarcoma cells, 1,25-dihydroxyvitamin D3 [1,25(OH)2D 3] regulates OPN at the transcriptional level resulting in increased steady state mRNA levels and increased production of OPN protein, maximal at 48 hours. Using ROS 17/2.8 cells as an osteoblast model, OPN was purified from culture medium after three hour treatments of either vehicle (ethanol) or 1,25(OH)2D3 via barium citrate precipitation followed by immunoaffinity chromatography. ^ Here, further evidence of regulation of OPN by 1,25(OH)2D 3 at the posttranslational level is presented. Prior to the up-regulation of OPN at the transcriptional level, 1,25(OH)2D3 induces a shift in OPN isoelectric point (pI) detected on two-dimensional gels from pI 4.6 to pI 5.1. Loading equal amounts of [32P]-labeled OPN recovered from ROS 17/2.8 cells exposed to 1,25(OH)2D3 or vehicle alone for three hours reveals that the shift from pI 4.6 to 5.1 is the result of reduced phosphorylation. Using structural analogs to 1,25(OH) 2D3, analog AT [25-(OH)-16-ene-23-yne-D3], which triggers Ca2+ influx through voltage sensitive Ca2+ channels but does not bind to the vitamin D receptor, mimicked the OPN pI shift while analog BT [1,25(OH)2-22-ene-24-cyclopropyl-D 3], which binds to the vitamin D receptor but does not allow Ca 2+ influx, did not. Inclusion of the Ca2+ channel blocker nifedipine also blocks the charge shift conversion of OPN. Further analysis of the signaling pathway initiated by 1,25(OH)2D3 reveals that inhibition of the cyclic 3′,5′ -adenosine monophosphate-dependent kinase, protein kinase A, or inhibition of the cyclic 3′,5′-guanine monophosphate-dependent kinase, protein kinase G, also prevents the charge shift conversion. ^ Isolation of OPN from rat femurs and tibiae provides evidence for the existence of these two OPN charge forms in vivo, evidenced by differential migration on isoelectric focusing gels and sodium dodecyl sulfate-polyacrylamide gels. Peptide sequencing of rat long bone fractions revealed the presence of a presumed dentin specific protein, dentin matrix protein-1 (DMP-1). Western blot analysis confirmed the existence of DMP-1 in these fractions. ^ Using the OPN charge forms in functional assays, it was determined that the charge forms have differential roles in both cell surface and mineralization functions. In cell attachment assays and Ca2+ influx assays using PC-3 prostate cancer cells, the pI 5.1 charge form of OPN was found to permit binding and increase intracellular Ca2+ concentrations of PC-3 cells. The increase in intracellular Ca2+ concentration was found to be integrin αvβ3-dependent. In mineralization assays, the pI 4.6 charge form of OPN promoted hydroxyapatite formation, while the pI 5.1 charge form had improved Ca2+ binding ability. ^ In conclusion, these findings suggest that 1,25(OH) 2D3 regulates OPN not only at the transcriptional level, but also plays a role in determination of the OPN phosphorylation state. The latter involves a short term (less than three hours) treatment and is associated with membrane-initiated Ca2+ influx. Functional assays utilizing the two OPN charge forms reveal the dependence of OPN post-translational state on its function. ^

Ser88 phosphorylation regulates dynein light chain 1 (DLC1) function in the mouse mammary gland

Dynein light chain 1 (DLC1) is a highly conserved and ubiquitously expressed protein which might have critical cellular function as total loss of DLC1 caused Drosophila embryonic death. Despite many proteins and RNAs interaction with it identified, DLC1's function(s) and regulation are largely unknown. Recently, DLC1 was identified as a physiological substrate of P21-activate kinase 1(Pak1) kinase from a human mammary cDNA library in a yeast-2-hybridization screening assay. Studies in primary human tumors and cell culture implicated that DLC1 could promote mammary cancerous phenotypes, and more importantly, Ser88 phosphorylation of DLC1by Pak1 kinase was found to be essential for DLC1's tumorigenic activities. Based on the above tissue culture studies, we hypothesized that Ser88 phosphorylation regulates DLC1. ^ To test this hypothesis, we generated two transgenic mouse models: MMTV-DLC1 and MMTV-DLC1-S88A mice with mammary specific expression of the DLC1 and DLC1-S88A cDNAs. Both of the transgenic mice mammary glands showed rare tumor incidence which indicated DLC1 alone may not be sufficient for tumorigenesis in vivo. However, these mice showed a significant alteration of mammary development. Mammary glands from the MMTV-DLC1 mice had hyperbranching and alveolar hyperplasia, with elevated cell proliferation. Intriguingly, these phenotypes were not seen in the mammary glands from the MMTV-S88A mice. Furthermore, while MMTV-DLC1 glands were normal during involution, MMTV-S88A mice showed accelerated mammary involution with increase apoptosis and altered expression of involution-associated genes. Further analysis of the MMTV-S88A glands showed they had increased steady state level of Bim protein which might be responsible for the early involution. Finally, our in vitro data showed that Ser88 phosphorylation abolished DLC1 dimer and consequently might disturb its interaction with Bim and destabilize Bim. ^ Collectively, our findings provided in vivo evidence that Ser88 phosphorylation of DLC1 can regulate DLC1's function. In addition, Ser88 phosphorylation might be critical for DLC1 dimer-monomer transition. ^

Control of the master virulence regulatory gene atxA in Bacillus anthracis

Transcription of the Bacillus anthracis structural genes for the anthrax toxin proteins and biosynthetic operon for capsule are positively regulated by AtxA, a transcription regulator with unique properties. Consistent with the role of atxA in virulence factor expression, a B. anthracis atxA-null mutant is avirulent in a murine model for anthrax. In batch culture, multiple signals impact atxA transcript levels, and the timing and steady state level of atxA expression is critical for optimal toxin and capsule synthesis. Despite the apparent complex control of atxA transcription, only one trans-acting protein, the transition state regulator AbrB, has been demonstrated to directly interact with the atxA promoter. The AbrB-binding site has been described, but additional cis-acting control sequences have not been defined. Using transcriptional lacZ fusions, electrophoretic mobility shift assays, and Western blot analysis, the cis-acting elements and trans-acting factors involved in regulation of atxA in B. anthracis strains containing either both virulence plasmids, pXO1 and pXO2, or only one plasmid, pXO1, were studied. This work demonstrates that atxA transcription from the major start site P1 is dependent upon a consensus sequence for the housekeeping sigma factor SigA, and an A+T-rich upstream element (UP-element) for RNA polymerase (RNAP). In addition, the data show that a trans-acting protein(s) other than AbrB negatively impacts atxA transcription when it binds specifically to a 9-bp palindrome within atxA promoter sequences located downstream of P1. Mutation of the palindrome prevents binding of the trans-acting protein(s) and results in a corresponding increase in AtxA and anthrax toxin production in a strain- and culture-dependent manner. The identity of the trans-acting repressor protein(s) remains elusive; however, phenotypes associated with mutation of the repressor binding site have revealed that the trans-acting repressor protein(s) indirectly controls B. anthracis development. Mutation of the repressor binding site results in misregulation and overexpression of AtxA in conditions conducive for development, leading to a marked sporulation defect that is both atxA- and pXO2-61-dependent. pXO2-61 is homologous to the sensor domain of sporulation sensor histidine kinases and is proposed to titrate an activating signal away from the sporulation phosphorelay when overexpressed by AtxA. These results indicate that AtxA is not only a master virulence regulator, but also a modulator of proper B. anthracis development. Also demonstrated in this work is the impact of the developmental regulators AbrB, Spo0A, and SigH on atxA expression and anthrax toxin production in a genetically incomplete (pXO1+, pXO2-) and genetically complete (pXO1+, pXO2+) strain background. AtxA and anthrax toxin production resulting from deletion of the developmental regulators are strain-dependent suggesting that factors on pXO2 are involved in control of atxA. The only developmental deletion mutant that resulted in a prominent and consistent strain-independent increase in AtxA protein levels was an abrB-null mutant. As a result of increased AtxA levels, there is early and increased production of anthrax toxins in an abrB-null mutant. In addition, the abrB-null mutant exhibited an increase in virulence in a murine model for anthrax. In contrast, virulence of the atxA promoter mutant was unaffected in a murine model for anthrax despite the production of 5-fold more AtxA than the abrB-null mutant. These results imply that AtxA is not the only factor impacting pathogenesis in an abrB-null mutant. Overall, this work highlights the complex regulatory network that governs expression of atxA and provides an additional role for AtxA in B. anthracis development.

External assessment of myocardial glucose metabolism with $\sp{18}$F-2 -deoxy-2 -fluoro-D-glucose

Despite the popularity of the positron emitting glucose analog, ($\sp{18}$F) -2-deoxy-2-fluoro-D-glucose (2FDG), for the noninvasive "metabolic imaging" of organs with positron emission tomography (PET), the physiological basis for the tracer has not been tested, and the potential of 2FDG for the rapid kinetic analysis of altered glucose metabolism in the intact heart has not been fully exploited. We, therefore, developed a quantitative method to characterize metabolic changes of myocardial glucose metabolism noninvasively and with high temporal resolution.^ The first objective of the work was to provide direct evidence that the initial steps in the metabolism of 2FDG are the same as for glucose and that 2FDG is retained by the tissue in proportion to the rate of glucose utilization. The second objective was to characterize the kinetic changes in myocardial glucose transport and phosphorylation in response to changes in work load, competing substrates, acute ischemia and reperfusion, and the addition of insulin. To assess changes in myocardial glucose metabolism isolated working rat hearts were perfused with glucose and 2FDG. Tissue uptake of 2FDG and the input function were measured on-line by external detection. The steady state rate of 2FDG phosphorylation was determined by graphical analysis of 2FDG time-activity curves.^ The rate of 2FDG uptake was linear with time and the tracer was retained in its phosphorylated form. Tissue accumulation of 2FDG decreased within seconds with a reduction in work load, in the presence of competing substrates, and during reperfusion after global ischemia. Thus, most interventions known to alter glucose metabolism induced rapid parallel changes in 2FDG uptake. By contrast, insulin caused a significant increase in 2FDG accumulation only in hearts from fasted animals when perfused at a sub-physiological work load. The mechanism for this phenomenon is not known but may be related to the existence of two different glucose transporter systems and/or glycogen metabolism in the myocardial cell.^ It is concluded that (1) 2FDG traces glucose uptake and phosphorylation in the isolated working rat heart; and (2) early and transient kinetic changes in glucose metabolism can be monitored with high temporal resolution with 2FDG and a simple positron coincidence counting system. The new method has revealed transients of myocardial glucose metabolism, which would have remained unnoticed with conventional methods. These transients are not only important for the interpretation of glucose metabolic PET scans, but also provide insights into mechanisms of glucose transport and phosphorylation in heart muscle. ^