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.

DEVELOPMENTAL AND CELLULAR FUNCTIONS OF DELTA-CATENIN

Catenins have diverse and powerful roles in embryogenesis, homeostasis or disease progression, as best exemplified by the well-known beta-catenin. The less studied delta-catenin likewise contains a central Armadillo-domain. In common with other p120 sub-class members, it acts in a variety of intracellular compartments and modulates cadherin stability, small GTPase activities and gene transcription. In mammals, delta-catenin exhibits neural specific expression, with its knock-out in mice correspondingly producing cognitive defects and synaptic dysfunctions. My work instead employed the amphibian, Xenopus laevis, to explore delta-catenin’s physiological functions in a distinct vertebrate system. Initial isolation and characterization indicated delta-catenin’s expression in Xenopus. Unlike the pattern observed for mammals, delta-catenin was detected in most adult Xenopus tissues, although enriched in embryonic structures of neural fate as visualized using RNA in-situ hybridization. To determine delta-catenin’s requirement in amphibian development, I employed anti-sense morpholinos to knock-down gene products, finding that delta-catenin depletion results in developmental defects in gastrulation, neural crest migration and kidney tubulogenesis, phenotypes that were specific based upon rescue experiments. In biochemical and cellular assays, delta-catenin knock-down reduced cadherin levels and cell adhesion, and impaired activation of RhoA and Rac1, small GTPases that regulate actin dynamics and morphogenetic movements. Indeed, exogenous C-cadherin, or dominant-negative RhoA or dominant-active Rac1, significantly rescued delta-catenin depletion. Thus, my results indicate delta-catenin’s essential roles in Xenopus development, with contributing functional links to cadherins and Rho family small G proteins. In examining delta-catenin’s nuclear roles, I identified delta-catenin as an interacting partner and substrate of the caspase-3 protease, which plays critical roles in apoptotic as well as non-apoptotic processes. Delta-catenin’s interaction with and sensitivity to caspase-3 was confirmed using assays involving its cleavage in vitro, as well as within Xenopus apoptotic extracts or mammalian cell lines. The cleavage site, a highly conserved caspase consensus motif (DELD) within Armadillo-repeat 6 of delta-catenin, was identified through peptide sequencing. Cleavage thus generates an amino- (1-816) and carboxyl-terminal (817-1314) fragment each containing about half of the central Armadillo-domain. I found that cleavage of delta-catenin both abolishes its association with cadherins, and impairs its ability to modulate small GTPases. Interestingly, the carboxyl-terminal fragment (817-1314) possesses a conserved putative nuclear localization signal that I found is needed to facilitate delta-catenin’s nuclear targeting. To probe for novel nuclear roles of delta-catenin, I performed yeast two-hybrid screening of a mouse brain cDNA library, resolving and then validating its interaction with an uncharacterized KRAB family zinc finger protein I named ZIFCAT. My results indicate that ZIFCAT is nuclear, and suggest that it may associate with DNA as a transcriptional repressor. I further determined that other p120 sub-class catenins are similarly cleaved by caspase-3, and likewise bind ZIFCAT. These findings potentially reveal a simple yet novel signaling pathway based upon caspase-3 cleavage of p120 sub-family members, facilitating the coordinate modulation of cadherins, small GTPases and nuclear functions. Together, my work suggested delta-catenin’s essential roles in Xenopus development, and has revealed its novel contributions to cell junctions (via cadherins), cytoskeleton (via small G proteins), and nucleus (via ZIFCAT). Future questions include the larger role and gene targets of delta-catenin in nucleus, and identification of upstream signaling events controlling delta-catenin’s activities in development or disease progression.

Denaturation and unfolding of human anaphylatoxin C3a: an unusually low covalent stability of its native disulfide bonds.

The complement C3a anaphylatoxin is a major molecular mediator of innate immunity. It is a potent activator of mast cells, basophils and eosinophils and causes smooth muscle contraction. Structurally, C3a is a relatively small protein (77 amino acids) comprising a N-terminal domain connected by 3 native disulfide bonds and a helical C-terminal segment. The structural stability of C3a has been investigated here using three different methods: Disulfide scrambling; Differential CD spectroscopy; and Reductive unfolding. Two uncommon features regarding the stability of C3a and the structure of denatured C3a have been observed in this study. (a) There is an unusual disconnection between the conformational stability of C3a and the covalent stability of its three native disulfide bonds that is not seen with other disulfide proteins. As measured by both methods of disulfide scrambling and differential CD spectroscopy, the native C3a exhibits a global conformational stability that is comparable to numerous proteins with similar size and disulfide content, all with mid-point denaturation of [GdmCl](1/2) at 3.4-5M. These proteins include hirudin, tick anticoagulant protein and leech carboxypeptidase inhibitor. However, the native disulfide bonds of C3a is 150-1000 fold less stable than those proteins as evaluated by the method of reductive unfolding. The 3 native disulfide bonds of C3a can be collectively and quantitatively reduced with as low as 1mM of dithiothreitol within 5 min. The fragility of the native disulfide bonds of C3a has not yet been observed with other native disulfide proteins. (b) Using the method of disulfide scrambling, denatured C3a was shown to consist of diverse isomers adopting varied extent of unfolding. Among them, the most extensively unfolded isomer of denatured C3a is found to assume beads-form disulfide pattern, comprising Cys(36)-Cys(49) and two disulfide bonds formed by two pair of consecutive cysteines, Cys(22)-Cys(23) and Cys(56)-Cys(57), a unique disulfide structure of polypeptide that has not been documented previously.

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.

GCN5 maintains genome stability during development

The histone acetyltransferase, GCN5, is essential for survival of mice during embryogenesis. GCN5 null embryos die early during development due to increased apoptosis. We have demonstrated that the increased apoptosis in associated with increased p53 protein levels. Loss of p53 rescues the embryonic apoptosis in the GCN5 null embryos. These results raised the question of what molecular trigger leads to p53 stabilization and cell death in the absence of GCN5. p53 is generally referred to as the gatekeeper of the cell, monitoring cellular responses to DNA damage, genotoxic stress, and other unfavorable conditions in the cell. Therefore, we examined individual cells in wild type and mutant embryos for gross chromosomal aberrations that might trigger a genome integrity checkpoint. Karyotype analysis indicates that approximately 30% of the cells in an E8.5 GCN5 null embryo display chromosomal aberrations, predominantly chromosomal end adhesions and associations. In wild type E8.5 embryos, only 6% of the cells have chromosomal aberrations. Recent data using telomeric FISH demonstrates that cells from GCN5 null embryos have a decreased telomeric signal. Telomere maintenance is essential for maintaining genome integrity. Telomeric defects are associated with loss of chromosomes and chromosomal rearrangements that can lead to detrimental gene fusions involved in many types of cancers. Little is known about the chromatin structures present near the telomeric ends, or whether any of the telomere-associated proteins are subject to post-translational modification such as acetylation. Our results are the first data to demonstrate the involvement of a histone acetyltransferase, GCN5, in maintaining genome integrity through telomere maintenance and/or capping. ^