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.

Interlinked dual-time feedback loops can enhance robustness to stochasticity and persistence of memory.

Multiple interlinked positive feedback loops shape the stimulus responses of various biochemical systems, such as the cell cycle or intracellular Ca2+ release. Recent studies with simplified models have identified two advantages of coupling fast and slow feedback loops. This dual-time structure enables a fast response while enhancing resistances of responses and bistability to stimulus noise. We now find that (1) the dual-time structure similarly confers resistance to internal noise due to molecule number fluctuations, and (2) model variants with altered coupling, which better represent some specific biochemical systems, share all the above advantages. We also develop a similar bistable model with coupling of a fast autoactivation loop to a slow loop. This model's topology was suggested by positive feedback proposed to play a role in long-term synaptic potentiation (LTP). The advantages of fast response and noise resistance are also present in this autoactivation model. Empirically, LTP develops resistance to reversal over approximately 1h . The model suggests this resistance may result from increased amounts of synaptic kinases involved in positive feedback.

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.

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.

MULTIPLE POSTTRANSCRIPTIONAL REGULATORY FEATURES CONTROL EXPRESSION OF ETHANOLAMINE UTILIZATION GENES IN ENTEROCOCCUS FAECALIS

Enterococcus faecalis is a Gram-positive bacterium that lives as a commensal organism in the mammalian gastrointestinal tract, but can behave as an opportunistic pathogen. Our lab discovered that mutation of the eutK gene attenuates virulence of E. faecalis in the C. elegans model host. eutK is part of the ethanolamine metabolic pathway which was previously unknown in E. faecalis. I discovered the presence of two unique posttranscriptional regulatory features that control expression of eut locus genes. The first feature I found is an AdoCBL riboswitch, a cis-acting RNA regulatory element that acts as a positive regulator of gene expression. The second feature I discovered is a unique two-component system, EutVW. The EutV response regulator contains an ANTAR family domain, which binds RNA to trigger transcriptional antitermination. I determined that induction of expression of several genes in the eut locus is dependent on ethanolamine, AdoCBL and the two-component system. AdoCBL and ethanolamine are both required for induction of eut locus gene expression. Additionally, I discovered eutG is regulated by a unique mechanism of antitermination. Both the AdoCBL riboswitch and EutV response regulator control the expression of the downstream gene eutG. EutV potentially acts through a novel antitermination mechanism in which a dimer of EutV binds to a pair of mRNA stem loops forming an antitermination complex. My data show a unique mechanism by which two environmental signals are integrated by two different posttranscriptional regulators to regulate a single locus.

Translational regulation of nuclear gene COX4 expression by mitochondrial content of phosphatidylglycerol and cardiolipin in Saccharomyces cerevisiae.

Previous results indicated that translation of four mitochondrion-encoded genes and one nucleus-encoded gene (COX4) is repressed in mutants (pgs1Delta) of Saccharomyces cerevisiae lacking phosphatidylglycerol and cardiolipin. COX4 translation was studied here using a mitochondrially targeted green fluorescence protein (mtGFP) fused to the COX4 promoter and its 5' and 3' untranslated regions (UTRs). Lack of mtGFP expression independent of carbon source and strain background was established to be at the translational level. The translational defect was not due to deficiency of mitochondrial respiratory function but was rather caused directly by the lack of phosphatidylglycerol and cardiolipin in mitochondrial membranes. Reintroduction of a functional PGS1 gene under control of the ADH1 promoter restored phosphatidylglycerol synthesis and expression of mtGFP. Deletion analysis of the 5' UTR(COX4) revealed the presence of a 50-nucleotide fragment with two stem-loops as a cis-element inhibiting COX4 translation. Binding of a protein factor(s) specifically to this sequence was observed with cytoplasm from pgs1Delta but not PGS1 cells. Using HIS3 and lacZ as reporters, extragenic spontaneous recessive mutations that allowed expression of His3p and beta-galactosidase were isolated, which appeared to be loss-of-function mutations, suggesting that the genes mutated may encode the trans factors that bind to the cis element in pgs1Delta cells.

3.88 A structure of cytoplasmic polyhedrosis virus by cryo-electron microscopy.

Cytoplasmic polyhedrosis virus (CPV) is unique within the Reoviridae family in having a turreted single-layer capsid contained within polyhedrin inclusion bodies, yet being fully capable of cell entry and endogenous RNA transcription. Biochemical data have shown that the amino-terminal 79 residues of the CPV turret protein (TP) is sufficient to bring CPV or engineered proteins into the polyhedrin matrix for micro-encapsulation. Here we report the three-dimensional structure of CPV at 3.88 A resolution using single-particle cryo-electron microscopy. Our map clearly shows the turns and deep grooves of alpha-helices, the strand separation in beta-sheets, and densities for loops and many bulky side chains; thus permitting atomic model-building effort from cryo-electron microscopy maps. We observed a helix-to-beta-hairpin conformational change between the two conformational states of the capsid shell protein in the region directly interacting with genomic RNA. We have also discovered a messenger RNA release hole coupled with the mRNA capping machinery unique to CPV. Furthermore, we have identified the polyhedrin-binding domain, a structure that has potential in nanobiotechnology applications.

LOCALIZATION OF DNA REPAIR FUNCTIONS TO THE NUCLEAR MATRIX

The nucleus of a eukaryotic cell contains both structural and functional elements that contribute to the controlled operation of the cell. In this context, functional components refers to those nuclear constituents that perform metabolic activities such as DNA replication and RNA transcription. Structural nuclear components, designated nuclear matrix, organize the DNA into loops or domains and appear to provide a framework for nuclear DNA organization. However, the boundary between structural and functional components is not clear cut as evinced by reports of associations between metabolic functions and the nuclear matrix. The studies reported here attempt to determine the relationship of another nuclear function, DNA repair, to the nuclear matrix.^ One objective of these studies was to study the initiation of DNA repair by directly measuring the UV-incision activities in human cells and determine the influence of various extractable nuclear components on these activities. The assay for incision activities required the development of a nuclear isolation protocol that produced nuclei with intact DNA; the conformation of the nuclear DNA and its physical characteristics in response to denaturing conditions were determined.^ The nuclei produced with this protocol were then used as substrates for endogenous UV-specific nuclease activities. The isolated nuclei were shown to contain activities that cause breaks in nuclear DNA in response to UV-irradiation. These UV-responsive activities were tightly associated with nuclear components, being unextractable with salt concentration of up to 0.6 M.^ The tight association of the incision activities with salt-extracted nuclei suggested that other repair function might also be associated with salt-stable components of the nucleus. The site of unscheduled DNA synthesis (UDS) was determined in salt-extracted nuclei (nucleoids) using autoradiography and fluorescent microscopy. UDS was found to occur in association with the nuclear matrix following low-doses (2.55 J/M('2)) of ultraviolet light, but the association became looser after higher doses of ultraviolet light (10-30 J/m('2)). ^

Molecular mechanisms of prostacyclin receptor signaling: The intracellular domains in G protein coupling

To better understand the mechanisms of how the human prostacyclin receptor (1P) mediates vasodilation and platelet anti-aggregation through Gs protein coupling, a strategy integrating multiple approaches including high resolution NMR experiments, synthetic peptide, fluorescence spectroscopy, molecular modeling, and recombinant protein was developed and used to characterize the structure/function relationship of important segments and residues of the IP receptor and the α-subunit of the Gs protein (Gαs). The first (iLP1) and third (iLP3) intracellular loops of the IP receptor, as well as the Gαs C-terminal domain, relevant to the Gs-mediated IP receptor signaling, were first identified by observation of the effects of the mini gene-expressed corresponding protein segments in HEK293 cells which co-expressed the receptor and Gαs. Evidence of the IP iLP1 domain interacted with the Gαs C-terminal domain was observed by fluorescence and NMR spectroscopic studies using a constrained synthetic peptide, which mimicked the IP iLP1 domain, and the synthetic peptide, which mimicked Gαs C-terminal domain. The solution structural models and the peptide-peptide interaction of the two synthetic protein segments were determined by high resolution NMR spectroscopy. The important residues in the corresponding domains of the IP receptor and the Gαs predicted by NMR chemical shift mapping were used to guide the identification of their protein-protein interaction in cells. A profile of the residues Arg42 - Ala48 of the IP iLP1 domain and the three residues Glu392 ∼ Leu394 of the Gαs C-terminal domain involved in the IP/Gs protein coupling were confirmed by recombinant proteins. The data revealed an intriguing speculation on the mechanisms of how the signal of the ligand-activated IP receptor is transmitted to the Gs protein in regulating vascular functions and homeostasis, and also provided substantial insights into other prostanoid receptor signaling. ^

Lipid-dependent topogenesis and function of membrane proteins

Membranes are essential for the integrity and function of the cell. The collective property of the lipid bilayer is critical in providing an optimal functioning environment for membrane proteins. The simple yet well-characterized bacterium Escherichia coli serves an ideal model system to study the function of specific lipids since its lipid content can be easily manipulated. The most abundant lipid in E. coli membrane is phosphatidylethanolamine (PE, 70-80%). A PE-lacking E. coli mutant displays a complex mixture of deficient phenotypes, suggesting a profound role for PE in different aspects of cell function. A novel role of PE as a topological and functional determinant for membrane proteins has been established using lactose permease (LacY) as a model protein. PE is found to be required for energy-dependent uphill transport process of LacY. In PE-lacking membranes, LacY undergoes a dramatic conformational change, and the first half of the protein adopts an inverted topology with respect to the bilayer plane. ^ The work reported here was initiated to understand the molecular properties of lipids that enable their function as topological and functional determinants for membrane proteins. A glycolipid, monoglucosyldiacylglycerol (MGlcDAG) which shares physicochemical similarities with PE, was introduced to PE-lacking E. coli membranes. The introduction of MGlcDAG suppresses many of the PE-deficient phenotypes, and in particular supports the function and native topology of LacY. ^ The lipid-sensitive topogenic signals encoded in the amino acid sequence of LacY were also identified. Native LacY adopts an inverted topology when synthesized without PE, but mutation of specific acidic residues in the cytoplasmic extra-membrane domains can prevent this inversion and supports a native topological organization of LacY in PE-lacking membranes. These results suggest that it is the interplay between the collective charge properties of the lipid bilayer and extra-membrane loops of protein that determines the final orientation of transmembrane domains. By comparing the similarities as well as differences between these two lipids, we established how specific physical and chemical properties of lipids influence various cell functions and elucidated the molecular basis for the novel role of lipids in determining membrane protein topology. ^