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.

Olfactory transduction is intrinsically noisy.

The sources of noise that limit olfactory signal detection were investigated in dissociated rat olfactory receptor cells. Near-threshold odorant-evoked currents exhibited large random fluctuation. However, similar fluctuations were observed in the absence of applied odorants when currents were induced by elevating the intracellular cyclic AMP concentration. This suggests that the fluctuations reflect noise intrinsic to the transduction mechanism, rather than the quantal nature of an odorant stimulus. For many odorants, this intrinsic noise may preclude the reliable detection of single odorant molecules.

Simultaneous determination of helical unwinding angles and intrinsic association constants in ligand–DNA complexes: The interaction between DNA and calichearubicin B

We present a helical unwinding assay for reversibly binding DNA ligands that uses closed circular DNA, topoisomerase I (Topo I), and two-dimensional agarose gel electrophoresis. Serially diluted Topo I relaxation reactions at constant DNA/ligand ratio are performed, and the resulting apparent unwinding of the closed circular DNA is used to calculate both ligand unwinding angle (φ) and intrinsic association constant (Ka). Mathematical treatment of apparent unwinding is formally analogous to that of apparent extinction coefficient data for optical binding titrations. Extrapolation to infinite DNA concentration yields the true unwinding angle of a given ligand and its association constant under Topo I relaxation conditions. Thus this assay delivers simultaneous structural and thermodynamic information describing the ligand–DNA complex. The utility of this assay has been demonstrated by using calichearubicin B (CRB), a synthetic hybrid molecule containing the anthraquinone chromophore of (DA) and the carbohydrate domain of calicheamicin γ1I. The unwinding angle for CRB calculated by this method is −5.3 ± 0.5°. Its Ka value is 0.20 × 106 M−1. For comparison, the unwinding angles of ethidium bromide and DA have been independently calculated, and the results are in agreement with canonical values for these compounds. Although a stronger binder to selected sites, CRB is a less potent unwinder than its parent compound DA. The assay requires only small amounts of ligand and offers an attractive option for analysis of DNA binding by synthetic and natural compounds.

Yeast HOS3 forms a novel trichostatin A-insensitive homodimer with intrinsic histone deacetylase activity

Histone deacetylases such as human HDAC1 and yeast RPD3 are trichostatin A (TSA)-sensitive enzymes that are members of large, multiprotein complexes. These contain specialized subunits that help target the catalytic protein to histones at the appropriate DNA regulatory element, where the enzyme represses transcription. To date, no deacetylase catalytic subunits have been shown to have intrinsic activity, suggesting that noncatalytic subunits of the deacetylase complex are required for their enzymatic function. In this paper we describe a novel yeast histone deacetylase HOS3 that is relatively insensitive to the histone deacetylase inhibitor TSA, forms a homodimer when expressed ectopically both in yeast and Escherichia coli, and has intrinsic activity when produced in the bacterium. Most HOS3 protein can be found associated with a larger complex in partially purified yeast nuclear extracts, arguing that the HOS3 homodimer may be dissociated from a very large nuclear structure during purification. We also demonstrate, using a combination of mass spectrometry, tandem mass spectrometry, and proteolytic digestion, that recombinant HOS3 has a distinct specificity in vitro for histone H4 sites K5 and K8, H3 sites K14 and K23, H2A site K7, and H2B site K11. We propose that while factors that interact with HOS3 may sequester the catalytic subunit at specific cellular sites, they are not required for HOS3 histone deacetylase activity.

Use of intrinsic binding energy for catalysis by an RNA enzyme

The contribution of several individual ribozyme⋅substrate base pairs to binding and catalysis has been investigated using hammerhead ribozyme substrates that were truncated at their 3′ or 5′ ends. The base pairs at positions 1.1–2.1 and 15.2–16.2, which flank the conserved core, each contribute 104-fold in the chemical step, without affecting substrate binding. In contrast, base pairs distal to the core contribute to substrate binding but have no effect on the chemical step. These results suggest a “fraying model” in which each ribozyme⋅substrate helix can exist in either an unpaired (“open”) state or a helical (“closed”) state, with the closed state required for catalysis. The base pairs directly adjacent to the conserved core contribute to catalysis by allowing the closed state to form. Once the number of base pairs is sufficient to ensure that the closed helical state predominates, additional residues provide stabilization of the helix, and therefore increase binding, but have no further effect on the chemical step. Remarkably, the >5 kcal/mol free energy contribution to catalysis from each of the internal base pairs is considerably greater than the free energy expected for formation of a base pair. It is suggested that this unusually large energetic contribution arises because free energy that is typically lost in constraining residues within a base pair is expressed in the transition state, where it is used for positioning. This extends the concept of “intrinsic binding energy” from protein to RNA enzymes, suggesting that intrinsic binding energy is a fundamental feature of biological catalysis.

Integration of multiple instructive cues by neural crest stem cells reveals cell-intrinsic biases in relative growth factor responsiveness

Growth factors can influence lineage determination of neural crest stem cells (NCSCs) in an instructive manner, in vitro. Because NCSCs are likely exposed to multiple signals in vivo, these findings raise the question of how stem cells would integrate such combined influences. Bone morphogenetic protein 2 (BMP2) promotes neuronal differentiation and glial growth factor 2 (GGF2) promotes glial differentiation; if NCSCs are exposed to saturating concentrations of both factors, BMP2 appears dominant. By contrast, if the cells are exposed to saturating concentrations of both BMP2 and transforming growth factor β1 (which promotes smooth muscle differentiation), the two factors appear codominant. Sequential addition experiments indicate that NCSCs require 48–96 hrs in GGF2 before they commit to a glial fate, whereas the cells commit to a smooth muscle fate within 24 hr in transforming growth factor β1. The delayed response to GGF2 does not reflect a lack of functional receptors; however, because the growth factor induces rapid mitogen-activated protein kinase phosphorylation in naive cells. Furthermore, GGF2 can attenuate induction of the neurogenic transcription factor mammalian achaete-scute homolog 1, by low doses of BMP2. This short-term antineurogenic influence of GGF2 is not sufficient for glial lineage commitment, however. These data imply that NCSCs exhibit cell-intrinsic biases in the timing and relative dosage sensitivity of their responses to instructive factors that influence the outcome of lineage decisions in the presence of multiple factors. The relative delay in glial lineage commitment, moreover, apparently reflects successive short-term and longer-term actions of GGF2. Such a delay may help to explain why glia normally differentiate after neurons, in vivo.

Transformation with human dihydrofolate reductase renders malaria parasites insensitive to WR99210 but does not affect the intrinsic activity of proguanil

Increasing resistance of Plasmodium falciparum malaria parasites to chloroquine and the dihydrofolate reductase (DHFR) inhibitors pyrimethamine and cycloguanil have sparked renewed interest in the antimalarial drugs WR99210 and proguanil, the cycloguanil precursor. To investigate suggestions that WR99210 and proguanil act against a target other than the reductase moiety of the P. falciparum bifunctional DHFR–thymidylate synthase enzyme, we have transformed P. falciparum with a variant form of human DHFR selectable by methotrexate. Human DHFR was found to fully negate the antiparasitic effect of WR99210, thus demonstrating that the only significant action of WR99210 is against parasite DHFR. Although the human enzyme also resulted in greater resistance to cycloguanil, no decrease was found in the level of susceptibility of transformed parasites to proguanil, thus providing evidence of intrinsic activity of this parent compound against a target other than DHFR. The transformation system described here has the advantage that P. falciparum drug-resistant lines are uniformly sensitive to methotrexate and will complement transformation with existing pyrimethamine-resistance markers in functional studies of P. falciparum genes. This system also provides an approach for screening and identifying novel DHFR inhibitors that will be important in combined chemotherapeutic formulations against malaria.

Use of intrinsic modes in biology: Examples of indicial response of pulmonary blood pressure to ± step hypoxia

Recently, a new method to analyze biological nonstationary stochastic variables has been presented. The method is especially suitable to analyze the variation of one biological variable with respect to changes of another variable. Here, it is illustrated by the change of the pulmonary blood pressure in response to a step change of oxygen concentration in the gas that an animal breathes. The pressure signal is resolved into the sum of a set of oscillatory intrinsic mode functions, which have zero “local mean,” and a final nonoscillatory mode. With this device, we obtain a set of “mean trends,” each of which represents a “mean” in a definitive sense, and together they represent the mean trend systematically with different degrees of oscillatory content. Correspondingly, the oscillatory content of the signal about any mean trend can be represented by a set of partial sums of intrinsic mode functions. When the concept of “indicial response function” is used to describe the change of one variable in response to a step change of another variable, we now have a set of indicial response functions of the mean trends and another set of indicial response functions to describe the energy or intensity of oscillations about each mean trend. Each of these can be represented by an analytic function whose coefficients can be determined by a least-squares curve-fitting procedure. In this way, experimental results are stated sharply by analytic functions.

δ-Tonoplast intrinsic protein defines unique plant vacuole functions

Plant cell vacuoles may have either storage or degradative functions. Vegetative storage proteins (VSPs) are synthesized in response to wounding and to developmental switches that affect carbon and nitrogen sinks. Here we show that VSPs are stored in a unique type of vacuole that is derived from degradative central vacuoles coincident with insertion of a new tonoplast intrinsic protein (TIP), δ-TIP, into their membranes. This finding demonstrates a tight coupling between the presence of δ-TIP and acquisition of a specialized storage function and indicates that TIP isoforms may determine vacuole identity.

Intrinsic responses to Borna disease virus infection of the central nervous system

Immune cells invading the central nervous system (CNS) in response to Borna disease virus (BDV) antigens are central to the pathogenesis of Borna disease (BD). We speculate that the response of the resident cells of the brain to infection may be involved in the sensitization and recruitment of these inflammatory cells. To separate the responses of resident cells from those of cells infiltrating from the periphery, we used dexamethasone to inhibit inflammatory reactions in BD. Treatment with dexamethasone prevented the development of clinical signs of BD, and the brains of treated animals showed no neuropathological lesions and a virtual absence of markers of inflammation, cell infiltration, or activation normally seen in the CNS of BDV-infected rats. In contrast, treatment with dexamethasone exacerbated the expression of BDV RNA, which was paralleled by a similarly elevated expression of mRNAs for egr-1, c-fos, and c-jun. Furthermore, dexamethasone failed to inhibit the increase in expression of mRNAs for tumor necrosis factor α, macrophage inflammatory protein 1β, interleukin 6, and mob-1, which occurs in the CNS of animals infected with BDV. Our findings suggest that these genes, encoding transcription factors, chemokines, and proinflammatory cytokines, might be directly activated in CNS resident cells by BDV. This result supports the hypothesis that the initial phase of the inflammatory response to BDV infection in the brain may be dependent upon virus-induced activation of CNS resident cells.

Role of intrinsic synaptic circuitry in collicular sensorimotor integration

The superficial gray layer of the superior colliculus contains a map that represents the visual field, whereas the underlying intermediate gray layer contains a vector map of the saccades that shift the direction of gaze. These two maps are aligned so that a particular region of the visual field is represented directly above the neurons that orient the highest acuity area of the retina toward that region. Although it has been proposed that the transmission of information from the visuosensory to the motor map plays an important role in the generation of visually guided saccades, experiments have failed to demonstrate any functional linkage between the two layers. We examined synaptic transmission between these layers in vitro by stimulating the superficial layer while using whole-cell patch-clamp methods to measure the responses of intermediate layer neurons. Stimulation of superficial layer neurons evoked excitatory postsynaptic currents in premotor cells. This synaptic input was columnar in organization, indicating that the connections between the layers link corresponding regions of the visuosensory and motor maps. Excitatory postsynaptic currents were large enough to evoke action potentials and often occurred in clusters similar in duration to the bursts of action potentials that premotor cells use to command saccades. Our results indicate the presence of functional connections between the superficial and intermediate layers and show that such connections could play a significant role in the generation of visually guided saccades.

The intrinsic radioresistance of glioblastoma-derived cell lines is associated with a failure of p53 to induce p21BAX expression

Radiation is the primary modality of therapy for all commonly occurring malignant brain tumors, including medulloblastoma and glioblastoma. These two brain tumors, however, have a distinctly different response to radiation therapy. Medulloblastoma is very sensitive to radiation therapy, whereas glioblastoma is highly resistant, and the long-term survival of medulloblastoma patients exceeds 50%, while there are few long-term survivors among glioblastoma patients. p53-mediated apoptosis is thought to be an important mechanism mediating the cytotoxic response of tumors to radiotherapy. In this study, we compared the response to radiation of five cell lines that have wild-type p53: three derived from glioblastoma and two derived from medulloblastoma. We found that the medulloblastoma-derived cell lines underwent extensive radiation-induced apoptotic cell death, while those from glioblastomas did not exhibit significant radiation-induced apoptosis. p53-mediated induction of p21BAX is thought to be a key component of the pathway mediating apoptosis after the exposure of cells to cytotoxins, and the expression of mRNA encoding p21BAX was correlated with these cell lines undergoing radiation-induced apoptosis. The failure of p53 to induce p21BAX expression in glioblastoma-derived cell lines is likely to be of biologic significance, since inhibition of p21BAX induction in medulloblastoma resulted in a loss of radiation-induced apoptosis, while forced expression of p21BAX in glioblastoma was sufficient to induce apoptosis. The failure of p53 to induce p21BAX in glioblastoma-derived cell lines suggests a distinct mechanism of radioresistance and may represent a critical factor in determining therapeutic responsiveness to radiation in glioblastomas.

Neuronal calcium sparks and intracellular calcium “noise”

Intracellular calcium ions are involved in many forms of cellular function. To accommodate so many control functions, a complex spatiotemporal organization of calcium signaling has developed. In both excitable and nonexcitable cells, calcium signaling was found to fluctuate. Sudden localized increases in the intracellular calcium concentration—or calcium sparks—were found in heart, striated and smooth muscle, Xenopus Laevis oocytes, and HeLa and P12 cells. In the nervous system, intracellular calcium ions were found important in key processes such as transmitter release, repetitive firing, and gene expression. Hence, we examined whether calcium sparks also exist in neurons. Using confocal laser-scanning microscopy and fluorescent probes, we found that calcium sparks exist in two types of neuronal preparations: the presynaptic boutons of the lizard neuromuscular junction and rat hippocampal neurons in cell culture. Control experiments exclude the possibility that these calcium sparks originate from instrumental or biological artifacts. Calcium sparks seem to be just the tip of the iceberg of a more general phenomenon of intracellular calcium “noise.” We speculate that calcium sparks and calcium noise may be of key importance in calcium signaling in the nervous system.

Perceptual learning reflects external noise filtering and internal noise reduction through channel reweighting

To investigate the nature of plasticity in the adult visual system, perceptual learning was measured in a peripheral orientation discrimination task with systematically varying amounts of external (environmental) noise. The signal contrasts required to achieve threshold were reduced by a factor or two or more after training at all levels of external noise. The strong quantitative regularities revealed by this novel paradigm ruled out changes in multiplicative internal noise, changes in transducer nonlinearites, and simple attentional tradeoffs. Instead, the regularities specify the mechanisms of perceptual learning at the behavioral level as a combination of external noise exclusion and stimulus enhancement via additive internal noise reduction. The findings also constrain the neural architecture of perceptual learning. Plasticity in the weights between basic visual channels and decision is sufficient to account for perceptual learning without requiring the retuning of visual mechanisms.

Ubiquitously Expressed Dynamin-II Has a Higher Intrinsic GTPase Activity and a Greater Propensity for Self-assembly Than Neuronal Dynamin-I

To begin to understand mechanistic differences in endocytosis in neurons and nonneuronal cells, we have compared the biochemical properties of the ubiquitously expressed dynamin-II isoform with those of neuron-specific dynamin-I. Like dynamin-I, dynamin-II is specifically localized to and highly concentrated in coated pits on the plasma membrane and can assemble in vitro into rings and helical arrays. As expected, the two closely related isoforms share a similar mechanism for GTP hydrolysis: both are stimulated in vitro by self-assembly and by interaction with microtubules or the SH3 domain-containing protein, grb2. Deletion of the C-terminal proline/arginine-rich domain from either isoform abrogates self-assembly and assembly-dependent increases in GTP hydrolysis. However, dynamin-II exhibits a ∼threefold higher rate of intrinsic GTP hydrolysis and higher affinity for GTP than dynamin-I. Strikingly, the stimulated GTPase activity of dynamin-II can be >40-fold higher than dynamin-I, due principally to its greater propensity for self-assembly and the increased resistance of assembled dynamin-II to GTP-triggered disassembly. These results are consistent with the hypothesis that self-assembly is a major regulator of dynamin GTPase activity and that the intrinsic rate of GTP hydrolysis reflects a dynamic, GTP-dependent equilibrium of assembly and disassembly.