Intrinsic noise in gene regulatory networks

Cells are intrinsically noisy biochemical reactors: low reactant numbers can lead to significant statistical fluctuations in molecule numbers and reaction rates. Here we use an analytic model to investigate the emergent noise properties of genetic systems. We find for a single gene that noise is essentially determined at the translational level, and that the mean and variance of protein concentration can be independently controlled. The noise strength immediately following single gene induction is almost twice the final steady-state value. We find that fluctuations in the concentrations of a regulatory protein can propagate through a genetic cascade; translational noise control could explain the inefficient translation rates observed for genes encoding such regulatory proteins. For an autoregulatory protein, we demonstrate that negative feedback efficiently decreases system noise. The model can be used to predict the noise characteristics of networks of arbitrary connectivity. The general procedure is further illustrated for an autocatalytic protein and a bistable genetic switch. The analysis of intrinsic noise reveals biological roles of gene network structures and can lead to a deeper understanding of their evolutionary origin.

RecA binding to a single double-stranded DNA molecule: A possible role of DNA conformational fluctuations

Most genetic regulatory mechanisms involve protein–DNA interactions. In these processes, the classical Watson–Crick DNA structure sometimes is distorted severely, which in turn enables the precise recognition of the specific sites by the protein. Despite its key importance, very little is known about such deformation processes. To address this general question, we have studied a model system, namely, RecA binding to double-stranded DNA. Results from micromanipulation experiments indicate that RecA binds strongly to stretched DNA; based on this observation, we propose that spontaneous thermal stretching fluctuations may play a role in the binding of RecA to DNA. This has fundamental implications for the protein–DNA binding mechanism, which must therefore rely in part on a combination of flexibility and thermal fluctuations of the DNA structure. We also show that this mechanism is sequence sensitive. Theoretical simulations support this interpretation of our experimental results, and it is argued that this is of broad relevance to DNA–protein interactions.

DHR3, an ecdysone-inducible early-late gene encoding a Drosophila nuclear receptor, is required for embryogenesis

Response to the steroid hormone ecdysone in Drosophila is controlled by genetic regulatory hierarchies that include eight members of the nuclear receptor protein family. The DHR3 gene, located within the 46F early-late ecdysone-inducible chromosome puff, encodes an orphan nuclear receptor that recently has been shown to exert both positive and negative regulatory effects in the ecdysone-induced genetic hierarchies at metamorphosis. We used a reverse genetics approach to identify 11 DHR3 mutants from a pool of lethal mutations in the 46F region on the second chromosome. Two DHR3 mutations result in amino acid substitutions within the conserved DNA binding domain. Analysis of DHR3 mutants reveals that DHR3 function is required to complete embryogenesis. All DHR3 alleles examined result in nervous system defects in the embryo.

The Ath5 proneural genes function upstream of Brn3 POU domain transcription factor genes to promote retinal ganglion cell development

During retinogenesis, the Xenopus basic helix–loop–helix transcription factor Xath5 has been shown to promote a ganglion cell fate. In the developing mouse and chicken retinas, gene targeting and overexpression studies have demonstrated critical roles for the Brn3 POU domain transcription factor genes in the promotion of ganglion cell differentiation. However, the genetic relationship between Ath5 and Brn3 genes is unknown. To understand the genetic regulatory network(s) that controls retinal ganglion cell development, we analyzed the relationship between Ath5 and Brn3 genes by using a gain-of-function approach in the chicken embryo. We found that during retinogenesis, the chicken Ath5 gene (Cath5) is expressed in retinal progenitors and in differentiating ganglion cells but is absent in terminally differentiated ganglion cells. Forced expression of both Cath5 and the mouse Ath5 gene (Math5) in retinal progenitors activates the expression of cBrn3c following central-to-peripheral and temporal-to-nasal gradients. As a result, similar to the Xath5 protein, both Cath5 and Math5 proteins have the ability to promote the development of ganglion cells. Moreover, we found that forced expression of all three Brn3 genes also can stimulate the expression of cBrn3c. We further found that Ath5 and Brn3 proteins are capable of transactivating a Brn3b promoter. Thus, these data suggest that the expression of cBrn3c in the chicken and Brn3b in the mouse is initially activated by Ath5 factors in newly generated ganglion cells and later maintained by a feedback loop of Brn3 factors in the differentiated ganglion cells.

Phylogenetic footprinting of transcription factor binding sites in proteobacterial genomes

Toward the goal of identifying complete sets of transcription factor (TF)-binding sites in the genomes of several gamma proteobacteria, and hence describing their transcription regulatory networks, we present a phylogenetic footprinting method for identifying these sites. Probable transcription regulatory sites upstream of Escherichia coli genes were identified by cross-species comparison using an extended Gibbs sampling algorithm. Close examination of a study set of 184 genes with documented transcription regulatory sites revealed that when orthologous data were available from at least two other gamma proteobacterial species, 81% of our predictions corresponded with the documented sites, and 67% corresponded when data from only one other species were available. That the remaining predictions included bona fide TF-binding sites was proven by affinity purification of a putative transcription factor (YijC) bound to such a site upstream of the fabA gene. Predicted regulatory sites for 2097 E.coli genes are available at http://www.wadsworth.org/resnres/bioinfo/.

Targeted deletion of the mouse POU domain gene Brn-3a causes selective loss of neurons in the brainstem and trigeminal ganglion, uncoordinated limb movement, and impaired suckling.

The Brn-3 subfamily of POU domain genes are expressed in sensory neurons and in select brainstem nuclei. Earlier work has shown that targeted deletion of the Brn-3b and Brn-3c genes produce, respectively, defects in the retina and in the inner ear. We show herein that targeted deletion of the Brn-3a gene results in defective suckling and in uncoordinated limb and trunk movements, leading to early postnatal death. Brn-3a (-/-) mice show a loss of neurons in the trigeminal ganglia, the medial habenula, the red nucleus, and the caudal region of the inferior olivary nucleus but not in the retina and dorsal root ganglia. In the trigeminal and dorsal root ganglia, but not in the retina, there is a marked decrease in the frequency of neurons expressing Brn-3b and Brn-3c, suggesting that Brn-3a positively regulates Brn-3b and Brn-3c expression in somatosensory neurons. Thus, Brn-3a exerts its major developmental effects in somatosensory neurons and in brainstem nuclei involved in motor control. The pheno-types of Brn-3a, Brn-3b, and Brn-3c mutant mice indicate that individual Brn-3 genes have evolved to control development in the auditory, visual, or somatosensory systems and that despite differences between these systems in transduction mechanisms, sensory organ structures, and central information processing, there may be fundamental homologies in the genetic regulatory events that control their development.

Inhibition of the association of RNA polymerase II with the preinitiation complex by a viral transcriptional repressor.

Transcriptional repression is an important component of regulatory networks that govern gene expression. In this report, we have characterized the mechanisms by which the immediate early protein 2 (IE2 or IE86), a master transcriptional regulator of human cytomegalovirus, down-regulates its own expression. In vitro transcription and DNA binding experiments demonstrate that IE2 blocks specifically the association of RNA polymerase II with the preinitiation complex. Although, to our knowledge, this is the first report to describe a eukaryotic transcriptional repressor that selectively impedes RNA polymerase II recruitment, we present data that suggest that this type of repression might be widely used in the control of transcription by RNA polymerase II.

The underlying pathway structure of biochemical reaction networks

Bioinformatics is yielding extensive, and in some cases complete, genetic and biochemical information about individual cell types and cellular processes, providing the composition of living cells and the molecular structure of its components. These components together perform integrated cellular functions that now need to be analyzed. In particular, the functional definition of biochemical pathways and their role in the context of the whole cell is lacking. In this study, we show how the mass balance constraints that govern the function of biochemical reaction networks lead to the translation of this problem into the realm of linear algebra. The functional capabilities of biochemical reaction networks, and thus the choices that cells can make, are reflected in the null space of their stoichiometric matrix. The null space is spanned by a finite number of basis vectors. We present an algorithm for the synthesis of a set of basis vectors for spanning the null space of the stoichiometric matrix, in which these basis vectors represent the underlying biochemical pathways that are fundamental to the corresponding biochemical reaction network. In other words, all possible flux distributions achievable by a defined set of biochemical reactions are represented by a linear combination of these basis pathways. These basis pathways thus represent the underlying pathway structure of the defined biochemical reaction network. This development is significant from a fundamental and conceptual standpoint because it yields a holistic definition of biochemical pathways in contrast to definitions that have arisen from the historical development of our knowledge about biochemical processes. Additionally, this new conceptual framework will be important in defining, characterizing, and studying biochemical pathways from the rapidly growing information on cellular function.

A Late Mitotic Regulatory Network Controlling Cyclin Destruction in Saccharomyces cerevisiae

Exit from mitosis requires the inactivation of mitotic cyclin-dependent kinase–cyclin complexes, primarily by ubiquitin-dependent cyclin proteolysis. Cyclin destruction is regulated by a ubiquitin ligase known as the anaphase-promoting complex (APC). In the budding yeast Saccharomyces cerevisiae, members of a large class of late mitotic mutants, including cdc15, cdc5, cdc14, dbf2, and tem1, arrest in anaphase with a phenotype similar to that of cells expressing nondegradable forms of mitotic cyclins. We addressed the possibility that the products of these genes are components of a regulatory network that governs cyclin proteolysis. We identified a complex array of genetic interactions among these mutants and found that the growth defect in most of the mutants is suppressed by overexpression of SPO12, YAK1, and SIC1 and is exacerbated by overproduction of the mitotic cyclin Clb2. When arrested in late mitosis, the mutants exhibit a defect in cyclin-specific APC activity that is accompanied by high Clb2 levels and low levels of the anaphase inhibitor Pds1. Mutant cells arrested in G1 contain normal APC activity. We conclude that Cdc15, Cdc5, Cdc14, Dbf2, and Tem1 cooperate in the activation of the APC in late mitosis but are not required for maintenance of that activity in G1.

Altered recognition mutants of the response regulator PhoB: A new genetic strategy for studying protein–protein interactions

Two-component regulatory systems require highly specific interactions between histidine kinase (transmitter) and response regulator (receiver) proteins. We have developed a novel genetic strategy that is based on tightly regulated synthesis of a given protein to identify domains and residues of an interacting protein that are critical for interactions between them. Using a reporter strain synthesizing the nonpartner kinase VanS under tight arabinose control and carrying a promoter-lacZ fusion activated by phospho-PhoB, we isolated altered recognition (AR) mutants of PhoB showing enhanced activation (phosphorylation) by VanS as arabinose-dependent Lac+ mutants. Changes in the PhoBAR mutants cluster in a “patch” near the proposed helix 4 of PhoB based on the CheY crystal structure (a homolog of the PhoB receiver domain) providing further evidence that helix 4 lies in the kinase-regulator interface. Based on the CheY structure, one mutant has an additional change in a region that may propagate a conformational change to helix 4. The overall genetic strategy described here may also be useful for studying interactions of other components of the vancomycin resistance and Pi signal transduction pathways, other two-component regulatory systems, and other interacting proteins. Conditionally replicative oriRR6Kγ attP “genome targeting” suicide plasmids carrying mutagenized phoB coding regions were integrated into the chromosome of a reporter strain to create mutant libraries; plasmids encoding mutant PhoB proteins were subsequently retrieved by P1-Int-Xis cloning. Finally, the use of similar genome targeting plasmids and P1-Int-Xis cloning should be generally useful for constructing genomic libraries from a wide array of organisms.

A functional genetic screen identifies regions at the C-terminal tail and death-domain of death-associated protein kinase that are critical for its proapoptotic activity

Death-associated protein kinase (DAP-kinase) is a Ca+2/calmodulin-regulated serine/threonine kinase with a multidomain structure that participates in apoptosis induced by a variety of signals. To identify regions in this protein that are critical for its proapoptotic activity, we performed a genetic screen on the basis of functional selection of short DAP-kinase-derived fragments that could protect cells from apoptosis by acting in a dominant-negative manner. We expressed a library of randomly fragmented DAP-kinase cDNA in HeLa cells and treated these cells with IFN-γ to induce apoptosis. Functional cDNA fragments were recovered from cells that survived the selection, and those in the sense orientation were examined further in a secondary screen for their ability to protect cells from DAP-kinase-dependent tumor necrosis factor-α-induced apoptosis. We isolated four biologically active peptides that mapped to the ankyrin repeats, the “linker” region, the death domain, and the C-terminal tail of DAP-kinase. Molecular modeling of the complete death domain provided a structural basis for the function of the death-domain-derived fragment by suggesting that the protective fragment constitutes a distinct substructure. The last fragment, spanning the C-terminal serine-rich tail, defined a new regulatory region. Ectopic expression of the tail peptide (17 amino acids) inhibited the function of DAP-kinase, whereas removal of this region from the complete protein caused enhancement of the killing activity, indicating that the C-terminal tail normally plays a negative regulatory role. Altogether, this unbiased screen highlighted functionally important regions in the protein and revealed an additional level of regulation of DAP-kinase apoptotic function that does not affect the catalytic activity.

Manipulating the genetic identity and biochemical surface properties of individual cells with electric-field-induced fusion

A method for cell–cell and cell–liposome fusion at the single-cell level is described. Individual cells or liposomes were first selected and manipulated either by optical trapping or by adhesion to a micromanipulator-controlled ultramicroelectrode. Spatially selective fusion of the cell–cell or cell–liposome pair was achieved by the application of a highly focused electric field through a pair of 5-μm o.d. carbon-fiber ultramicroelectrodes. The ability to fuse together single cells opens new possibilities in the manipulation of the genetic and cellular makeup of individual cells in a controlled manner. In the study of cellular networks, for example, the alteration of the biochemical identity of a selected cell can have a profound effect on the behavior of the entire network. Fusion of a single liposome with a target cell allows the introduction of the liposomal content into the cell interior as well as the addition of lipids and membrane proteins onto the cell surface. This cell–liposome fusion represents an approach to the manipulation of the cytoplasmic contents and surface properties of single cells. As an example, we have introduced a membrane protein (γ-glutamyltransferase) reconstituted in liposomes into the cell plasma membrane.

A yeast genetic system for selecting small molecule inhibitors of protein–protein interactions in nanodroplets

Cellular processes are mediated by complex networks of molecular interactions. Dissection of their role most commonly is achieved by using genetic mutations that alter, for example, protein–protein interactions. Small molecules that accomplish the same result would provide a powerful complement to the genetic approach, but it generally is believed that such molecules are rare. There are several natural products, however, that illustrate the feasibility of this approach. Split-pool synthesis now provides a simple mechanical means to prepare vast numbers of complex, even natural product-like, molecules individually attached to cell-sized polymer beads. Here, we describe a genetic system compatible with split-pool synthesis that allows the detection of cell-permeable, small molecule inhibitors of protein–protein interactions in 100- to 200-nl cell culture droplets, prepared by a recently described technique that arrays large numbers of such droplets. These “nanodroplets” contain defined media, cells, and one or more beads containing ≈100 pmol of a photoreleasable small molecule and a controlled number of cells. The engineered Saccharomyces cerevisiae cells used in this study express two interacting proteins after induction with galactose whose interaction results in cell death in the presence of 5-fluoroorotic acid (inducible reverse two-hybrid assay). Disruption of the interaction by a small molecule allows growth, and the small molecule can be introduced into the system hours before induction of the toxic interaction. We demonstrate that the interaction between the activin receptor R1 and the immunophilin protein FKBP12 can be disrupted by the small molecule FK506 at nanomolar concentrations in nanodroplets. This system should provide a general method for selecting cell-permeable ligands that can be used to study the relevance of protein–protein interactions in living cells or organisms.

Genomic and genetic dissection of an archaeal regulon

The extremely halophilic archaeon Halobacterium sp. NRC-1 can grow phototrophically by means of light-driven proton pumping by bacteriorhodopsin in the purple membrane. Here, we show by genetic analysis of the wild type, and insertion and double-frame shift mutants of Bat that this transcriptional regulator coordinates synthesis of a structural protein and a chromophore for purple membrane biogenesis in response to both light and oxygen. Analysis of the complete Halobacterium sp. NRC-1 genome sequence showed that the regulatory site, upstream activator sequence (UAS), the putative binding site for Bat upstream of the bacterio-opsin gene (bop), is also present upstream to the other Bat-regulated genes. The transcription regulator Bat contains a photoresponsive cGMP-binding (GAF) domain, and a bacterial AraC type helix–turn–helix DNA binding motif. We also provide evidence for involvement of the PAS/PAC domain of Bat in redox-sensing activity by genetic analysis of a purple membrane overproducer. Five additional Bat-like putative regulatory genes were found, which together are likely to be responsible for orchestrating the complex response of this archaeon to light and oxygen. Similarities of the bop-like UAS and transcription factors in diverse organisms, including a plant and a γ-proteobacterium, suggest an ancient origin for this regulon capable of coordinating light and oxygen responses in the three major branches of the evolutionary tree of life. Finally, sensitivity of four of five regulon genes to DNA supercoiling is demonstrated and correlated to presence of alternating purine–pyrimidine sequences (RY boxes) near the regulated promoters.

Methylation of the Protein Phosphatase 2A Catalytic Subunit Is Essential for Association of Bα Regulatory Subunit But Not SG2NA, Striatin, or Polyomavirus Middle Tumor Antigen

Binding of different regulatory subunits and methylation of the catalytic (C) subunit carboxy-terminal leucine 309 are two important mechanisms by which protein phosphatase 2A (PP2A) can be regulated. In this study, both genetic and biochemical approaches were used to investigate regulation of regulatory subunit binding by C subunit methylation. Monoclonal antibodies selectively recognizing unmethylated C subunit were used to quantitate the methylation status of wild-type and mutant C subunits. Analysis of 13 C subunit mutants showed that both carboxy-terminal and active site residues are important for maintaining methylation in vivo. Severe impairment of methylation invariably led to a dramatic decrease in Bα subunit binding but not of striatin, SG2NA, or polyomavirus middle tumor antigen (MT) binding. In fact, most unmethylated C subunit mutants showed enhanced binding to striatin and SG2NA. Certain carboxy-terminal mutations decreased Bα subunit binding without greatly affecting methylation, indicating that Bα subunit binding is not required for a high steady-state level of C subunit methylation. Demethylation of PP2A in cell lysates with recombinant PP2A methylesterase greatly decreased the amount of C subunit that could be coimmunoprecipitated via the Bα subunit but not the amount that could be coimmunoprecipitated with Aα subunit or MT. When C subunit methylation levels were greatly reduced in vivo, Bα subunits were found complexed exclusively to methylated C subunits, whereas striatin and SG2NA in the same cells bound both methylated and unmethylated C subunits. Thus, C subunit methylation is critical for assembly of PP2A heterotrimers containing Bα subunit but not for formation of heterotrimers containing MT, striatin, or SG2NA. These findings suggest that methylation may be able to selectively regulate the association of certain regulatory subunits with the A/C heterodimer.