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.

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.

A reduced model clarifies the role of feedback loops and time delays in the Drosophila circadian oscillator.

Although several detailed models of molecular processes essential for circadian oscillations have been developed, their complexity makes intuitive understanding of the oscillation mechanism difficult. The goal of the present study was to reduce a previously developed, detailed model to a minimal representation of the transcriptional regulation essential for circadian rhythmicity in Drosophila. The reduced model contains only two differential equations, each with time delays. A negative feedback loop is included, in which PER protein represses per transcription by binding the dCLOCK transcription factor. A positive feedback loop is also included, in which dCLOCK indirectly enhances its own formation. The model simulated circadian oscillations, light entrainment, and a phase-response curve with qualitative similarities to experiment. Time delays were found to be essential for simulation of circadian oscillations with this model. To examine the robustness of the simplified model to fluctuations in molecule numbers, a stochastic variant was constructed. Robust circadian oscillations and entrainment to light pulses were simulated with fewer than 80 molecules of each gene product present on average. Circadian oscillations persisted when the positive feedback loop was removed. Moreover, elimination of positive feedback did not decrease the robustness of oscillations to stochastic fluctuations or to variations in parameter values. Such reduced models can aid understanding of the oscillation mechanisms in Drosophila and in other organisms in which feedback regulation of transcription may play an important role.

Negative regulation of Drosophila Transformer2 by SR splicing factors

The Drosophila Transformer-2 (Tra2) protein activates the splicing of doublesex and fruitless pre-mRNA and represses M1 intron splicing in its own RNA in male germline. The M1 retention is part of negative feedback mechanism that controls Tra2 protein synthesis. However it is not known how the M1 intron is repressed or why Tra2 activates splicing of some RNAs while repressing splicing in others. Here we show that Tra2 and SR protein Rbp1 function together to specifically repress M1 splicing in vitro through the same intronic silencer by binding independently to distinct sites. The role of Rbp1 in M1 repression in vivo was validated by the finding that increased expression of Rbp1 in S2 cells promotes M1 retention. Furthermore, Tra2 blocks prespliceosomal A complex formation, a step corresponding to U2 snRNP recruitment to the branchpoint. High levels of Tra2 repression require an upstream enhancer. Together, we propose that the complex formed by Tra2 and Rbp1 on the silencer achieves splicing repression by blocking the recognition of the branchpoint or antagonizing enhancer function. ^ In addition, both splicing regulatory activities of Tra2 are essential developmental events, doublesex splicing is the key for Drosophila sex determination in the soma, while M1 retention occurs in the male germline and is necessary for spermatogenesis. However, active Tra2 is expressed ubiquitously. So another issue we have studied is how Tra2 accomplishes negative and positive splicing regulation in a tissue-specific fashion. Surprisingly, we found that nuclear extract from somatically-derived S2 cells support M1 repression in vitro. This led us to hypothesize that no germline specific factor is required and that high levels of Tra2 expression in the male germline is sufficient to trigger M1 retention. To test it, I examined whether increased expression of Tra2 could promote M1 retention in cells outside male germline. My results show that increased Tra2 expression promotes M1 retention in somatically-derived S2 cells as well as in the somatic tissues of living flies. These results show that somatic tissues are capable of supporting M1 repression but do not normally do so because the low levels of Tra2 do not trigger negative feedback regulation. ^

Transcriptional regulation of the {\it neu\/} oncogene and its negative regulation by the c-{\it myc\/} oncogene

The first part of my research involved the characterization of the neu gene promoter. I subcloned a 2.2-kb sequence located upstream to the extreme 5$\sp\prime$ end of the neu gene, in front of the bacterial reporter gene, chloramphenicol acetyltransferase (CAT). Transfection of this construct into different cell lines and subsequent CAT assays demonstrated that this 2.2-kb fragment was functional as a promoter. A series of deletion constructs was engineered to study the contribution of different fragments to transcription. Subcloning of individual fragments was followed by a cotransfection competition experiment, which demonstrated the involvement of protein factors interacting with the promoter. A gel retardation assay was also performed to show the physical binding of protein factors to the promoter. The combined results suggested that both positively and negatively acting protein factors are involved in interacting with different regions of the promoter, contributing to the overall transcription activity. My findings provide an insight into the regulation of neu gene expression, which in turn provides the tools to understand the molecular mechanisms of overexpression of the neu gene in some breast cancer and ovarian cancer cell lines.^ In the second part of my research, I discovered that another oncogene, c-myc, was able to reverse the transformed morphology that was induced by the neu oncogene. Utilizing the promoter constructs that I made, I was able to show that the c-myc oncogene has a negative regulatory effect on the expression of the neu oncogene. Further studies suggested that c-myc is able to lower the effective concentration of a positive factor(s) that interact with a 139-bp fragment of the neu gene promoter. These findings may provide a direct evidence of the long suspected role of the c-myc gene in transcriptional regulation. The neu gene may very well be the first identified mammalian target gene that is regulated by the c-myc oncogene. Since c-myc is known to be stimulated by various mitogenic signals and the neu gene is likely to be a growth factor receptor, it is possible that c-myc, when stimulated by the signal transduction pathway of the neu gene, would function as a negative feedback regulator on the neu gene receptor. (Abstract shortened with permission of author.) ^

Simulation of Drosophila circadian oscillations, mutations, and light responses by a model with VRI, PDP-1, and CLK.

A model of Drosophila circadian rhythm generation was developed to represent feedback loops based on transcriptional regulation of per, Clk (dclock), Pdp-1, and vri (vrille). The model postulates that histone acetylation kinetics make transcriptional activation a nonlinear function of [CLK]. Such a nonlinearity is essential to simulate robust circadian oscillations of transcription in our model and in previous models. Simulations suggest that two positive feedback loops involving Clk are not essential for oscillations, because oscillations of [PER] were preserved when Clk, vri, or Pdp-1 expression was fixed. However, eliminating positive feedback by fixing vri expression altered the oscillation period. Eliminating the negative feedback loop in which PER represses per expression abolished oscillations. Simulations of per or Clk null mutations, of per overexpression, and of vri, Clk, or Pdp-1 heterozygous null mutations altered model behavior in ways similar to experimental data. The model simulated a photic phase-response curve resembling experimental curves, and oscillations entrained to simulated light-dark cycles. Temperature compensation of oscillation period could be simulated if temperature elevation slowed PER nuclear entry or PER phosphorylation. The model makes experimental predictions, some of which could be tested in transgenic Drosophila.

Role of protein kinase C in short-term sensitization of {\it Aplysia}

Plasticity at the connections between sensory neurons and their follower cells in Aplysia has been used extensively as a model system to examine mechanisms of simple forms of learning, such as sensitization. Sensitization is induced, at least in part, by the transmitter serotonin (5-HT) and expressed in several forms, including facilitation of sensorimotor connections. Spike broadening has been believed to be a key mechanism underlying facilitation of nondepressed synapses. Previously, this broadening was believed to be dependent primarily on cAMP/protein kinase A (PKA)-mediated reduction of a noninactivating, relatively voltage-independent K$\sp{+}$ current termed the S-K$\sp+$ current (I$\sb{\rm K{,}S}$). Recent evidence, however, suggests that 5-HT-induced somatic spike broadening is composed of at least two components: a cAMP-dependent, rapidly developing component and a cAMP-independent, slowly developing component.^ Phorbol esters, activators of protein kinase C (PKC), mimicked the cAMP-independent component of 5-HT-induced broadening. Staurosporine, which inhibits PKC, had little effect on the rapidly developing component of 5-HT-induced broadening, but inhibited significantly the slowly developing component. These results suggest that PKC is involved in the cAMP-independent component of 5-HT-induced broadening. The membrane currents responsible for the slowly developing component of broadening were examined. Activation of PKC mimicked, and partially occluded, 5-HT-induced modulation of membrane currents above 0 mV, where a voltage-dependent K$\sp+$ current (I$\sb{\rm K{,}V}$) is significantly activated. This modulation was complex because it was associated with a reduction in the magnitude of I$\sb{\rm K{,}V}$, as well as a slowing of both activation and inactivation kinetics of I$\sb{\rm K{,}V}$. These results support the hypothesis that PKC modulates I$\sb{\rm K{,}V}$ and that this modulation contributes to the slowly developing component of 5-HT-induced broadening. Based on these results and others, a new scheme for 5-HT-induced spike broadening is proposed in which the modulatory effects are mediated via two second messenger/protein kinase systems converging and diverging on multiple ionic conductances.^ The relationship between spike broadening and synaptic facilitation was also examined. Pharmacological reduction of I$\sb{\rm K{,}V}$ by low concentrations of 4-aminopyridine (4-AP) led to spike broadening and facilitation of the nondepressed sensorimotor connections, indicating that spike broadening via the reduction of I$\sc{K,V}$ can facilitate the synaptic connection. Further analyses, however, revealed that 4-AP-induced facilitation has qualitative differences from 5-HT- and PKC-induced facilitation. These results suggest that 5-HT- and PKC-induced facilitation of nondepressed synapses is mediated, at least in part, by spike-duration independent (SDI) processes. Under certain conditions, the PKC inhibitor, staurosporine, significantly inhibited the 5-HT-induced facilitation of sensorimotor connections.^ Finally, it was found that activation of PKC increased a basal level of cAMP and that PKC caused desensitization of the 5-HT receptor, which may be a possible negative feedback mechanism through which an extracellular ligand, 5-HT, is regulated. These results suggest that these two second messenger/protein kinase pathways can interact in the sensory neuron. Thus, neuronal plasticity that may contribute to learning and memory appears to involve several complex and interactive processes. ^

Functional architecture of mammalian horizontal cells

In the rabbit retina, there are two kinds of horizontal cells (HCs). The A-type HC is a large axonless cell which contacts cones exclusively. The B-type HC is an axon bearing cell. While the somatic dendrites of B-type HCs also contact cones, the axon expands into an elaborately branched structure, the axon terminal (AT), which contacts a large number of rods. It is difficult to label the different HCs selectively by immunochemical methods. Therefore, we developed dye injection methods to label each type of HC. Then it was possible, (1) to describe the detailed structure of the AT (2) to identify the glutamate receptors mediating cone input to A and B-type HCs and rod input to ATs and (3) to test the hypothesis that the B-type HCs are coupled via Cx57 gap junctions. ^ To obtain well filled examples of single HCs, it was necessary to block gap junction coupling to stop the spread of Neurobiotin through the network. We used dye coupling in A-type HCs to screen a series of potential gap junction antagonists. One of these compounds, meclofenamic acid (MFA), was potent, water soluble and easily reversible. This compound may be a useful tool to manipulate gap junction coupling. ^ In the presence of MFA, Neurobiotin passed down the axon of B-type HCs to reveal the detailed structure of the AT. We observed that only one AT ending entered each rod spherule invagination. This observation was confirmed by calculation and two dye injections. ^ Glutamate is the neurotransmitter used by both rods and cones. AMPA receptors were colocalized with the dendrites of A and B-type HCs at each cone pedicle. In addition, AMPA receptors were located on the AT ending at each rod spherule. Thus rod and cone input to HCs is mediated by AMPA receptors. ^ A-type and B-type HCs may express different connexins because they have different dye-coupling properties. Recently, we found that connexin50 (Cx50) is expressed by A-type HCs. B-type HCs and B-type ATs are also independently coupled. Cx57 was expressed in the OPL and double label studies showed that Cx 57 was colocalized with the AT matrix but not with the somatic dendrites of B-type HCs. ^ In summary, we have identified a useful gap junction antagonist, MFA. There is one AT ending at each rod spherule, rods inputs to ATs is mediated by AMPA receptors and coupling in the AT matrix is mediated by Cx57. This confirms that HCs with different properties use distinct connexins. The properties of ATs described in this research are consistent. The connections and properties reported here suggest that ATs functions as rod HCs and provide a negative feedback signal to rods. ^

DETERMINING THE GENOTYPE-PHENOTYPE CONNECTION IN SYNTHETIC INDUCIBLE GENE EXPRESSION SYSTEMS

Introduction Gene expression is an important process whereby the genotype controls an individual cell’s phenotype. However, even genetically identical cells display a variety of phenotypes, which may be attributed to differences in their environment. Yet, even after controlling for these two factors, individual phenotypes still diverge due to noisy gene expression. Synthetic gene expression systems allow investigators to isolate, control, and measure the effects of noise on cell phenotypes. I used mathematical and computational methods to design, study, and predict the behavior of synthetic gene expression systems in S. cerevisiae, which were affected by noise. Methods I created probabilistic biochemical reaction models from known behaviors of the tetR and rtTA genes, gene products, and their gene architectures. I then simplified these models to account for essential behaviors of gene expression systems. Finally, I used these models to predict behaviors of modified gene expression systems, which were experimentally verified. Results Cell growth, which is often ignored when formulating chemical kinetics models, was essential for understanding gene expression behavior. Models incorporating growth effects were used to explain unexpected reductions in gene expression noise, design a set of gene expression systems with “linear” dose-responses, and quantify the speed with which cells explored their fitness landscapes due to noisy gene expression. Conclusions Models incorporating noisy gene expression and cell division were necessary to design, understand, and predict the behaviors of synthetic gene expression systems. The methods and models developed here will allow investigators to more efficiently design new gene expression systems, and infer gene expression properties of TetR based systems.

Altered cholesterol metabolism in a hypercholesterolemia -resistant rabbit colony

A colony of rabbits has been developed at the University of Texas Medical School at Houston that is resistant to dietary-induced hypercholesterolemia. The liver of resistant rabbits had higher levels of ($\sp{125}$I) $\beta$-VLDL binding and 3-hydroxy-3-methylglutaryl (HMGCoA) reductase activity, but lower acyl coenzyme A:cholesterol acyltransferase (ACAT) activity than normal rabbits. Direct quantitation of intracellular cholesterol content of the liver revealed that the resistant rabbits had $<$10% of the intracellular free cholesterol present in normal rabbits. Fibroblasts isolated from normal and resistant rabbits exhibited differences in ($\sp{125}$I) LDL binding, HMGCoA reductase activity and ACAT activity that were similar to those found in the liver. No structural differences were found in the LDL receptor of normal and resistant fibroblasts that would account for the increased binding capacity of the resistant cells. The regulation of LDL receptor levels by exogenous oxygenated sterols was similar in normal and resistant fibroblasts. The regulation of LDL receptor binding capacity by LDL was attenuated in the resistant compared to normal fibroblasts, suggesting that the resistant fibroblasts have an alternate pathway for processing lipoprotein-derived cholesterol. Sterol-balance studies revealed that the cholesterol-fed resistant rabbits increased lithocholic acid excretion compared to the basal state, and had higher levels of deoxycholic acid excretion than cholesterol-fed normal rabbits. In addition, the specific activity and mRNA levels of cholesterol 7$\alpha$-hydroxylase (C7$\alpha$H) were higher in resistant rabbits than normal rabbits, suggesting that the increased bile acid excretion was due to an increase in bile acid synthesis. Increased clearance of cholesterol relieves the negative feedback inhibition cholesterol exerts on expression of the LDL receptor. The number of cell surface LDL receptors is then increased in resistant rabbits and allows rapid clearance of lipoproteins from the plasma compartment, thereby reducing plasma cholesterol levels. The low intracellular cholesterol level also relieves the negative feedback inhibition cholesterol exerts on HMGCoA reductase. Increased synthesis of cholesterol from acetate provides cells with cholesterol for bile acid synthesis and/or homeostasis. The activity of ACAT is then decreased due to the flux of cholesterol through the bile acid synthetic pathways. ^