A model of the roles of essential kinases in the induction and expression of late long-term potentiation.

The induction of late long-term potentiation (L-LTP) involves complex interactions among second-messenger cascades. To gain insights into these interactions, a mathematical model was developed for L-LTP induction in the CA1 region of the hippocampus. The differential equation-based model represents actions of protein kinase A (PKA), MAP kinase (MAPK), and CaM kinase II (CAMKII) in the vicinity of the synapse, and activation of transcription by CaM kinase IV (CAMKIV) and MAPK. L-LTP is represented by increases in a synaptic weight. Simulations suggest that steep, supralinear stimulus-response relationships between stimuli (e.g., elevations in [Ca(2+)]) and kinase activation are essential for translating brief stimuli into long-lasting gene activation and synaptic weight increases. Convergence of multiple kinase activities to induce L-LTP helps to generate a threshold whereby the amount of L-LTP varies steeply with the number of brief (tetanic) electrical stimuli. The model simulates tetanic, -burst, pairing-induced, and chemical L-LTP, as well as L-LTP due to synaptic tagging. The model also simulates inhibition of L-LTP by inhibition of MAPK, CAMKII, PKA, or CAMKIV. The model predicts results of experiments to delineate mechanisms underlying L-LTP induction and expression. For example, the cAMP antagonist RpcAMPs, which inhibits L-LTP induction, is predicted to inhibit ERK activation. The model also appears useful to clarify similarities and differences between hippocampal L-LTP and long-term synaptic strengthening in other systems.

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.

Loss of DNA polymerase zeta enhances spontaneous tumorigenesis.

Mammalian genomes encode at least 15 distinct DNA polymerases, functioning as specialists in DNA replication, DNA repair, recombination, or bypass of DNA damage. Although the DNA polymerase zeta (polzeta) catalytic subunit REV3L is important in defense against genotoxins, little is known of its biological function. This is because REV3L is essential during embryogenesis, unlike other translesion DNA polymerases. Outstanding questions include whether any adult cells are viable in the absence of polzeta and whether polzeta status influences tumorigenesis. REV3L-deficient cells have properties that could influence the development of neoplasia in opposing ways: markedly reduced damage-induced point mutagenesis and extensive chromosome instability. To answer these questions, Rev3L was conditionally deleted from tissues of adult mice using MMTV-Cre. Loss of REV3L was tolerated in epithelial tissues but not in the hematopoietic lineage. Thymic lymphomas in Tp53(-/-) Rev3L conditional mice occurred with decreased latency and higher incidence. The lymphomas were populated predominantly by Rev3L-null T cells, showing that loss of Rev3L can promote tumorigenesis. Remarkably, the tumors were frequently oligoclonal, consistent with accelerated genetic changes in the absence of Rev3L. Mammary tumors could also arise from Rev3L-deleted cells in both Tp53(+/+) and Tp53(+/-) backgrounds. Mammary tumors in Tp53(+/-) mice deleting Rev3L formed months earlier than mammary tumors in Tp53(+/-) control mice. Prominent preneoplastic changes in glandular tissue adjacent to these tumors occurred only in mice deleting Rev3L and were associated with increased tumor multiplicity. Polzeta is the only specialized DNA polymerase yet identified that inhibits spontaneous tumor 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.

Differential gene expression in multistep tumorigenesis using an {\it in vitro\/} human cell model

Using a human terato-carcinoma cell line, PA-1, the functional role of the oncogenes and tumor suppressor gene involved in the multistep process of carcinogenesis have been analyzed. The expression of AP-2 was strongly correlated with the susceptibility to ras transformation. The differential responsiveness to growth factors between stage 1 ras resistant cells and stage 2 ras susceptible cells was observed, indicating that the ability of stage 2 cells to respond to the mutated ras oncogenes in transformation correlated with the ability to be stimulated by certain growth factors. Using differential screening of cDNA libraries, a number of differentially expressed cDNA clones was isolated. One of those, clone 12, is overexpressed in ras transformed stage 3 cells. The amino acid sequence of clone 12 is almost identical to a mouse LLrep3 gene that was growth-regulated, and 78% similar to a yeast ribosomal protein S4. These results suggest that the S4 gene may be involved in regulation of growth. Clone 9 is expressed in stage 1 ras resistant cells (3.5-kb and 3.0-kb transcripts) but the expression of this clone in stage 2 ras susceptible cells and stage 3 ras-transformed cells is greatly diminished. The expression of this cDNA clone was increased to at least five fold in ras resistant cells and nontumorigenic hybrids treated with retinoic acid but not increased in retinoic acid treated ras susceptible cells, ras transformed cells and the tumorigenic segregants. Partial sequence of this clone showed no homology to the sequences in Genbank. These findings suggest that clone 9 could be a suppressor gene or the genes that are involved in the biochemical pathway of tumor suppression or neurogenic differentiation. The apparent pleiotropic effect of the loss of this suppressor gene function support Harris' proposal that tumor suppressor genes regulate differentiation. The tumor suppressor gene may act as negative regulator of tumor growth by controlling gene expression in differentiation. ^

Function and formation of $\beta$1,4-galactosyltransferase-specific adhesions during growth and differentiation of F9 embryonal carcinoma cells

Galactosyltransferase (GalTase) is localized in the Golgi, where it functions in oligosaccharide synthesis, as well as on the cell surface where it serves as a cell adhesion molecule. GalTase-specific adhesions are functional in a number of important biological events, including F9 embryonal carcinoma (EC) cell adhesions. GalTase-based adhesions are formed by recognition and binding to terminal N-acetylglucosamine (GlcNAc) residues on its glycoprotein counterpart on adjacent cell surfaces. The object of this work has been to investigate the formation and function of GalTase-specific adhesions during F9 cell growth and differentiation. We initially investigated GalTase synthesis during differentiation and found that the increase in GalTase activity was specific for the Golgi compartment; surface GalTase levels remained constant during differentiation. These data indicated that the increase in cell adhesions expected with increased cell-matrix interaction in differentiated F9 cells is not the consequence of increased surface GalTase expression and, more interestingly, that the two pools of GalTase are under differential regulation. Synthesis and recognition of the consociate glycoprotein component was next investigated. Surface GalTase recognized several surface glycoproteins in a pattern that changes with differentiation. Uvomorulin, lysosome-associated membrane protein-1 (LAMP-1), and laminin were recognized by surface GalTase and are, therefore, potential components in GalTase-specific adhesions. Furthermore, these interactions were aberrant in an adhesion-defective F9 cell line that results, at least in part, from abnormal oligosaccharide synthesis. The function played by surface GalTase in growth and induction of differentiation was examined. Inhibition of surface GalTase function by a panel of reagents inhibited F9 cell growth. GalTase expression at both the transcription and protein levels were differentially regulated during the cell cycle, with surface expression greatest in the G1 phase. Disruption of GalTase adhesion by exposure to anti-GalTase antibodies during this period resulted in extension of the G2 phase, a result similar to that seen with agents known to inhibit growth and induce differentiation. Finally, other studies have suggested that a subset of cell adhesion molecules have the capability to induce differentiation in EC cells systems. We have determined in F9 cells that dissociating GalTase adhesion by galactosylation of and release of the consociate glycoproteins induces differentiation, as defined by increased laminin synthesis. The ability to induce differentiation by surface galactosylation was greatest in cells grown in cultures promoting cell-cell adhesions, relative to cultures with minimal cell-cell interactions. ^

Transcriptional analysis of liver-specific expression and differential interleukin-6 response of rat kininogen genes

Analyses of rat T1 kininogen gene/chloramphenicol acetyltransferase (T1K/CAT) constructs revealed two regions important for tissue-specific and induced regulation of T1 kininogen.^ Although the T1 kininogen gene is inducible by inflammatory cytokines, a highly homologous K kininogen gene is minimally responsive. Moreover, the basal expression of a KK/CAT construct was 5- to 7-fold higher than that of the analogous T1K/CAT construct. To examine the molecular basis of this differential regulation, a series of promoter swapping experiments was carried out. Our transfection results showed that at least two regions in the K kininogen gene are important for its high basal expression: a distal 19-bp region (C box) constituted a binding site for CCAAT/enhancer binding protein (C/EBP) family proteins and a proximal 66-bp region contained two adjacent binding sites for hepatocyte nuclear factor-3 (HNF-3). The distal HNF-3 binding site from the K kininogen promoter demonstrated a stronger affinity than that from the T1 kininogen promoter. Since C/EBP and HNF-3 are highly enriched in the liver and known to enhance transcription of liver-specific genes, differential binding affinities of these factors accounted for the higher basal expression of the K kininogen gene.^ In contrast to the K kininogen C box, the T1 kininogen C box does not bind C/EBP presumably due to their two-nucleotide divergence. This sequence divergence, however, converts it to a consensus binding sequence for two IL-6-inducible transcription factors--IL-6 response element binding protein and acute-phase response factor. To functionally determine whether C box sequences are important for their differential acute-phase response, T1 and K kininogen C boxes were swapped and analyzed after transfection into Hep3B cells. Our results showed that the T1 kininogen C box is indeed one of the IL-6 response elements in T1 kininogen promoter. Furthermore, its function can be modulated by a 5$\sp\prime$-adjacent C/EBP-binding site (B box) whose mutation significantly reduced the overall induced activity. Moreover, this B box is the target site for binding and transactivation of another IL-6 inducible transcription factor C/EBP$\delta.$ Evolutionary divergence of a few critical nucleotides can either lead to subtle changes in the binding affinities of a given transcription factor or convert a binding sequence for a constitutive factor to a site recognized by an inducible factor. (Abstract shortened by UMI.) ^

The differential expression of the $\beta\sb2$-adrenergic receptor modifies agonist-dependent adenylylcyclase activation

There have been multiple reports which indicate that variations in $\beta$AR expression affect the V$\sb{\rm max}$ observed for the agonist-dependent activation of adenylylcyclase. This observation has been ignored by most researchers when V$\sb{\rm max}$ values obtained for wild type and mutant receptors are compared. Such an imprecise analysis may lead to erroneous conclusions concerning the ability of a receptor to activate adenylylcyclase. Equations were derived from the Cassel-Selinger model of GTPase activity and Tolkovsky and Levitzki's Collision Coupling model which predict that the EC$\sb{50}$ and V$\sb{\rm max}$ for the activation of adenylylcyclase are a function of receptor number. Experimental results for L cell clones in which either hamster or human $\beta$AR were transfected at varying levels showed that EC$\sb{50}$ decreases and V$\sb{\rm max}$ increases as receptor number increases. Comparison of these results with simulations obtained from the equations describing EC$\sb{50}$ and V$\sb{\rm max}$ showed a close correlation. This documents that the kinetic parameters of adenylylcyclase activation change with the level of receptor expression and relates this phenomenon to a theoretical framework concerning the mechanisms involved in $\beta$AR signal transduction.^ One of the terms used in the equations which expressed the EC$\sb{50}$ and V$\sb{\rm max}$ as a function of receptor number is coupling efficiency, defined as $\rm k\sb1/k\sb{-1}$. Calculation of $\rm k\sb1/k\sb{-1}$ can be accomplished for wild type receptors with the easily measured experimental values of agonist K$\sb{\rm d}$, EC$\sb{50}$ and receptor number. This was demonstrated for hamster $\beta$AR which yielded a coupling efficiency of 0.15 $\pm$ 0.003 and human $\beta$AR which yielded a coupling efficiency of 0.90 $\pm$ 0.031. $\rm k\sb1/k\sb{-1}$ replaces the traditional qualitative evaluation of the ability to activate adenylylcyclase, which utilizes V$\sb{\rm max}$ without correction for variation in receptor number, with a quantitative definition that more accurately describes the ability of $\beta$AR to couple to G$\sb{\rm s}$.^ The equations which express the EC$\sb{50}$ and V$\sb{\rm max}$ for adenylylcyclase activation as a function of receptor number and coupling efficiency were tested to determine whether they could accurately simulate the changes seen in these parameters during desensitization. Data from original desensitization experiments and data from the literature (24,25,52,54,83) were compared to simulated changes in EC$\sb{50}$ and V$\sb{\rm max}$. In a variety of systems the predictions of the equations were consistent with the changes observed in EC$\sb{50}$ and V$\sb{\rm max}$. In addition reductions in the calculated value of $\rm k\sb1/k\sb{-1}$ was shown to correlate well with $\beta$AR phosphorylation and to be minimally affected by sequestration and down-regulation. ^

Structural equation modeling of the medical outcome study: Short Form 36 (SF-36)

The factorial validity of the SF-36 was evaluated using confirmatory factor analysis (CFA) methods, structural equation modeling (SEM), and multigroup structural equation modeling (MSEM). First, the measurement and structural model of the hypothesized SF-36 was explicated. Second, the model was tested for the validity of a second-order factorial structure, upon evidence of model misfit, determined the best-fitting model, and tested the validity of the best-fitting model on a second random sample from the same population. Third, the best-fitting model was tested for invariance of the factorial structure across race, age, and educational subgroups using MSEM.^ The findings support the second-order factorial structure of the SF-36 as proposed by Ware and Sherbourne (1992). However, the results suggest that: (a) Mental Health and Physical Health covary; (b) general mental health cross-loads onto Physical Health; (c) general health perception loads onto Mental Health instead of Physical Health; (d) many of the error terms are correlated; and (e) the physical function scale is not reliable across these two samples. This hierarchical factor pattern was replicated across both samples of health care workers, suggesting that the post hoc model fitting was not data specific. Subgroup analysis suggests that the physical function scale is not reliable across the "age" or "education" subgroups and that the general mental health scale path from Mental Health is not reliable across the "white/nonwhite" or "education" subgroups.^ The importance of this study is in the use of SEM and MSEM in evaluating sample data from the use of the SF-36. These methods are uniquely suited to the analysis of latent variable structures and are widely used in other fields. The use of latent variable models for self reported outcome measures has become widespread, and should now be applied to medical outcomes research. Invariance testing is superior to mean scores or summary scores when evaluating differences between groups. From a practical, as well as, psychometric perspective, it seems imperative that construct validity research related to the SF-36 establish whether this same hierarchical structure and invariance holds for other populations.^ This project is presented as three articles to be submitted for publication. ^

THE EFFECT OF GABAMIMETICS ON NEUROTRANSMITTER RECEPTOR BINDING AND FUNCTION IN RAT BRAIN

Gamma-aminobutyric acid (GABA) is a major inhibitory neurotransmitter in the central nervous system and alterations in central GABAergic transmission may contribute to the symptoms of a number of neurological and psychiatric disorders. Because of this relationship, numerous laboratories are attempting to develop agents which will selectively enhance GABA neurotransmission in brain. Due to these efforts, several promising compounds have recently been discovered. Should these drugs prove to be clinically effective, they will be used to treat chronic neuropsychiatric disabilities and, therefore, will be administered for long periods of time. Accordingly, the present investigation was undertaken to determine the neurochemical consequences of chronic activation of brain GABA systems in order to better define the therapeutic potential and possible side-effect liability of GABAmimetic compounds.^ Chronic (15 day) administration to rats of low doses of amino-oxyacetic acid (AOAA, 10 mg/kg, once daily), isonicotinic acid hydrazide (20 mg/kg, b.i.d.), two non-specific inhibitors of GABA-T, the enzyme which catabolizes GABA in brain, or (gamma)-acetylenic GABA (10 mg/kg, b.i.d.) a catalytic inhibitor of this enzyme, resulted in a significant elevation of brain and CSF GABA content throughout the course of treatment. In addition, chronic administration of these drugs, as well as the direct acting GABA receptor agonists THIP (8 mg/kg, b.i.d.) or kojic amine (18 mg/kg, b.i.d.) resulted in a significant increase in dopamine receptor number and a significant decrease in GABA receptor number in the corpus striatum of treated animals as determined by standard in vitro receptor binding techniques. Changes in the GABA receptor were limited to the corpus striatum and occurred more rapidly than did alterations in the dopamine receptor. The finding that dopamine-mediated stereotypic behavior was enhanced in animals treated chronically with AOAA suggested that the receptor binding changes noted in vitro have some functional consequence in vitro.^ Coadministration of atropine (a muscarinic cholinergic receptor antagonist) blocked the GABA-T inhibitor-induced increase in striatal dopamine receptors but was without effect on receptor alterations seen following chronic administration of direct acting GABA receptor agonists. Atropine administration failed to influence the drug-induced decreases in striatal GABA receptors.^ Other findings included the discovery that synaptosomal high affinity ('3)H-choline uptake, an index of cholinergic neuronal activity, was significantly increased in the corpus striatum of animals treated acutely, but not chronically, with GABAmimetics.^ It is suggested that the dopamine receptor supersensitivity observed in the corpus striatum of animals following long-term treatment with GABAmimetics is a result of the chronic inhibition of the nigrostriatal dopamine system by these drugs. Changes in the GABA receptor, on the other hand, are more likely due to a homospecific regulation of these receptors. An hypothesis based on the different sites of action of GABA-T inhibitors vis-a-vis the direct acting GABA receptor agonists is proposed to account for the differential effect of atropine on the response to these drugs.^ The results of this investigation provide new insights into the functional interrelationships that exist in the basal ganglia and suggest that chronic treatment with GABAmimetics may produce extrapyramidal side-effects in man. In addition, the constellation of neurochemical changes observed following administration of these drugs may be a useful guide for determining the GABAmimetic properties of neuropharmacological agents. ^

14-3-3 ZETA OVEREXPRESSION SERVES AS A NOVEL MOLECULAR SWITCH TURNING TGF-BETA FROM TUMOR SUPPRESSOR TO TUMOR PROMOTER

TGF-β plays an important role in differentiation and tissue morphogenesis as well as cancer progression. However, the role of TGF-β in cancer is complicate. TGF-β has primarily been recognized as tumor suppressor, because it can directly inhibit cell proliferation of normal and premalignant epithelial cell. However, in the last stage of tumor progression, TGF-β functions as tumor promoter to enhance tumor cells metastatic dissemination and expands metastatic colonies. Currently, the mechanism of how TGF-β switches its role from tumor suppressor to promoter still remains elusive. Here we identify that overexpression of 14-3-3ζ inhibits TGF-β’s cell cytostatic program through destabilizing p53 in non-transformed human mammary epithelial cells. Mechanistically, we found that 14-3-3ζ overexpression leads to 14-3-3σ downregulation, thereby activates PI3K/Akt signaling pathway and degrades p53, and further inhibits TGF-β induced p21 expression and cell cytostatic function. In addition, we found that overexpression of 14-3-3ζ promotes TGF-β induced breast cancer cells bone metastatic colonization through stabilizing Gli2, which is an important co-transcriptional factor for p-smad2 to activate PTHrP expression and bone osteolytic effect. Taken together, we reveal a novel mechanism that 14-3-3ζ dictates the tumor suppressor or metastases promoter activities of TGF-β signaling pathway through switching p-smad2 binding partner from p53 to Gli2. The expected results will not only provide us the better understanding of the important role of 14-3-3ζ in the early stage of breast cancer development, but also deeply impact our knowledge of signaling mechanisms underlying the complex roles of TGF-β in cancer, which will give us a more accurate strategy to determine when and how anti-TGF-β targeted therapy might be effective.