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.

Kinetic coupling between electron and proton transfer in cytochrome c oxidase: simultaneous measurements of conductance and absorbance changes.

Bovine heart cytochrome c oxidase is an electron-current driven proton pump. To investigate the mechanism by which this pump operates it is important to study individual electron- and proton-transfer reactions in the enzyme, and key reactions in which they are kinetically and thermodynamically coupled. In this work, we have simultaneously measured absorbance changes associated with electron-transfer reactions and conductance changes associated with protonation reactions following pulsed illumination of the photolabile complex of partly reduced bovine cytochrome c oxidase and carbon monoxide. Following CO dissociation, several kinetic phases in the absorbance changes were observed with time constants ranging from approximately 3 microseconds to several milliseconds, reflecting internal electron-transfer reactions within the enzyme. The data show that the rate of one of these electron-transfer reactions, from cytochrome a3 to a on a millisecond time scale, is controlled by a proton-transfer reaction. These results are discussed in terms of a model in which cytochrome a3 interacts electrostatically with a protonatable group, L, in the vicinity of the binuclear center, in equilibrium with the bulk through a proton-conducting pathway, which determines the rate of proton transfer (and indirectly also of electron transfer). The interaction energy of cytochrome a3 with L was determined independently from the pH dependence of the extent of the millisecond-electron transfer and the number of protons released, as determined from the conductance measurements. The magnitude of the interaction energy, 70 meV (1 eV = 1.602 x 10(-19) J), is consistent with a distance of 5-10 A between cytochrome a3 and L. Based on the recently determined high-resolution x-ray structures of bovine and a bacterial cytochrome c oxidase, possible candidates for L and a physiological role for L are discussed.

Intrinsic compressibility and volume compression in solvated proteins by molecular dynamics simulation at high pressure.

Constant pressure and temperature molecular dynamics techniques have been employed to investigate the changes in structure and volumes of two globular proteins, superoxide dismutase and lysozyme, under pressure. Compression (the relative changes in the proteins' volumes), computed with the Voronoi technique, is closely related with the so-called protein intrinsic compressibility, estimated by sound velocity measurements. In particular, compression computed with Voronoi volumes predicts, in agreement with experimental estimates, a negative bound water contribution to the apparent protein compression. While the use of van der Waals and molecular volumes underestimates the intrinsic compressibilities of proteins, Voronoi volumes produce results closer to experimental estimates. Remarkably, for two globular proteins of very different secondary structures, we compute identical (within statistical error) protein intrinsic compressions, as predicted by recent experimental studies. Changes in the protein interatomic distances under compression are also investigated. It is found that, on average, short distances compress less than longer ones. This nonuniform contraction underlines the peculiar nature of the structural changes due to pressure in contrast with temperature effects, which instead produce spatially uniform changes in proteins. The structural effects observed in the simulations at high pressure can explain protein compressibility measurements carried out by fluorimetric and hole burning techniques. Finally, the calculation of the proteins static structure factor shows significant shifts in the peaks at short wavenumber as pressure changes. These effects might provide an alternative way to obtain information concerning compressibilities of selected protein regions.

ATPase activity of Escherichia coli Rep helicase crosslinked to single-stranded DNA: implications for ATP driven helicase translocation.

To examine the coupling of ATP hydrolysis to helicase translocation along DNA, we have purified and characterized complexes of the Escherichia coli Rep protein, a dimeric DNA helicase, covalently crosslinked to a single-stranded hexadecameric oligodeoxynucleotide (S). Crosslinked Rep monomers (PS) as well as singly ligated (P2S) and doubly ligated (P2S2) Rep dimers were characterized. The equilibrium and kinetic constants for Rep dimerization as well as the steady-state ATPase activities of both PS and P2S crosslinked complexes were identical to the values determined for un-crosslinked Rep complexes formed with dT16. Therefore, ATP hydrolysis by both PS and P2S complexes are not coupled to DNA dissociation. This also rules out a strictly unidirectional sliding mechanism for ATP-driven translocation along single-stranded DNA by either PS or the P2S dimer. However, ATP hydrolysis by the doubly ligated P2S2 Rep dimer is coupled to single-stranded DNA dissociation from one subunit of the dimer, although loosely (low efficiency). These results suggest that ATP hydrolysis can drive translocation of the dimeric Rep helicase along DNA by a "rolling" mechanism where the two DNA binding sites of the dimer alternately bind and release DNA. Such a mechanism is biologically important when one subunit binds duplex DNA, followed by subsequent unwinding.

An intrinsic genetic program for autonomous differentiation of muscle cells in the ascidian embryo.

The B-line presumptive muscle cells of ascidian embryos have extensive potential for self-differentiation dependent on determinants prelocalized in the myoplasm of fertilized eggs. Ascidian larval muscle cells therefore provide an experimental system with which to explore an intrinsic genetic program for autonomous specification of embryonic cells. Experiments with egg fragments suggested that maternal mRNAs are one of the components of muscle determinants. Expression of larval muscle actin genes begins as early as the 32-cell stage, prior to the developmental fate restriction of the cells. The timing of initiation of the actin gene expression proceeds the expression of an ascidian homologue of vertebrate MyoD by a few hours. Mutations in the proximal E-box of the 5' flanking region of the actin genes did not alter the promoter activity for muscle-specific expression of reporter gene. These results, together with results of deletion constructs of fusion genes, suggest that muscle determinants regulate directly, or indirectly via regulatory factors other than MyoD, the transcription of muscle-specific structural genes leading to the terminal differentiation.

Subthreshold facilitation and suppression in primary visual cortex revealed by intrinsic signal imaging.

Neurons in primary visual cortex (area 17) respond vigorously to oriented stimuli within their receptive fields; however, stimuli presented outside the suprathreshold receptive field can also influence their responses. Here we describe a fundamental feature of the spatial interaction between suprathreshold center and subthreshold surround. By optical imaging of intrinsic signals in area 17 in response to a stimulus border, we show that a given stimulus generates activity primarily in iso-orientation domains, which extend for several millimeters across the cortical surface in a manner consistent with the architecture of long-range horizontal connections in area 17. By mapping the receptive fields of single neurons and imaging responses from the same cortex to stimuli that include or exclude the aggregate suprathreshold receptive field, we show that intrinsic signals strongly reveal the subthreshold surround contribution. Optical imaging and single-unit recording both demonstrate that the relative contrast of center and surround stimuli regulates whether surround interactions are facilitative or suppressive: the same surround stimulus facilitates responses when center contrast is low, but suppresses responses when center contrast is high. Such spatial interactions in area 17 are ideally suited to contribute to phenomena commonly regarded as part of "higher-level" visual processing, such as perceptual "popout" and "filling-in."

Phosphorylation by protein kinase C of serine-23 of the alpha-1 subunit of rat Na+,K(+)-ATPase affects its conformational equilibrium.

Phosphorylation of the alpha-1 subunit of rat Na+,K(+)-ATPase by protein kinase C has been shown previously to decrease the activity of the enzyme in vitro. We have now undertaken an investigation of the mechanism by which this inhibition occurs. Analysis of the phosphorylation of recombinant glutathione S-transferase fusion proteins containing putative cytoplasmic domains of the protein, site-directed mutagenesis, and two-dimensional peptide mapping indicated that protein kinase C phosphorylated the alpha-1 subunit of the rat Na+,K(+)-ATPase within the extreme NH2-terminal domain, on serine-23. The phosphorylation of this residue resulted in a shift in the equilibrium toward the E1 form, as measured by eosin fluorescence studies, and this was associated with a decrease in the apparent K+ affinity of the enzyme, as measured by ATPase activity assays. The rate of transition from E2 to E1 was apparently unaffected by phosphorylation by protein kinase C. These results, together with previous studies that examined the effects of tryptic digestion of Na+,K(+)-ATPase, suggest that the NH2-terminal domain of the alpha-1 subunit, including serine-23, is involved in regulating the activity of the enzyme.

Mutations in the rotated abdomen locus affect muscle development and reveal an intrinsic asymmetry in Drosophila.

In bilateral animals, the left and right sides of the body usually present asymmetric structures, the genetic bases of whose generation are still largely unknown [CIBA Foundation (1991) Biological Asymmetry and Handedness, CIBA Foundation Symposium 162 (Wiley, New York), pp. 1-327]. In Drosophila melanogaster, mutations in the rotated abdomen (rt) locus cause a clockwise helical rotation of the body. Even null alleles are viable but exhibit defects in embryonic muscle development, rotation of the whole larval body, and helical staggering of cuticular patterns in abdominal segments of the adult. rotated abdomen is expressed in the embryonic mesoderm and midgut but not in the ectoderm; it encodes a putative integral membrane glycoprotein (homologous to key yeast mannosyltransferases). Mesodermal cells defective in O-glycosylation lead to an impaired larval muscular system. We propose that the staggering of the adult abdominal segments would be a consequence of the relaxation of intrinsic rotational torque of muscle architecture, preventing the colateral alignment of the segmental histoblast cells during their proliferation at metamorphosis.

Expression of corticotropin-releasing factor in inflamed tissue is required for intrinsic peripheral opioid analgesia.

Immune cell-derived opioid peptides can activate opioid receptors on peripheral sensory nerves to inhibit inflammatory pain. The intrinsic mechanisms triggering this neuroimmune interaction are unknown. This study investigates the involvement of endogenous corticotropin-releasing factor (CRF) and interleukin-1beta (IL-1). A specific stress paradigm, cold water swim (CWS), produces potent opioid receptor-specific antinociception in inflamed paws of rats. This effect is dose-dependently attenuated by intraplantar but not by intravenous alpha-helical CRF. IL-1 receptor antagonist is ineffective. Similarly, local injection of antiserum against CRF, but not to IL-1, dose-dependently reverses this effect. Intravenous anti-CRF is only inhibitory at 10(4)-fold higher concentrations and intravenous CRF does not produce analgesia. Pretreatment of inflamed paws with an 18-mer 3'-3'-end inverted CRF-antisense oligodeoxynucleotide abolishes CWS-induced antinociception. The same treatment significantly reduces the amount of CRF extracted from inflamed paws and the number of CRF-immunostained cells without affecting gross inflammatory signs. A mismatch oligodeoxynucleotide alters neither the CWS effect nor CRF immunoreactivity. These findings identify locally expressed CRF as the predominant agent to trigger opioid release within inflamed tissue. Endogenous IL-1, circulating CRF or antiinflammatory effects, are not involved. Thus, an intact immune system plays an essential role in pain control, which is important for the understanding of pain in immunosuppressed patients with cancer or AIDS.

Determination of the transiently lowered pKa of the retinal Schiff base during the photocycle of bacteriorhodopsin.

Reprotonation of the transiently deprotonated retinal Schiff base in the bacteriorhodopsin photocycle is greatly slowed when the proton donor Asp-96 is removed with site-specific mutagenesis, but its rate is restored upon adding azide or other weak acids such as formate and cyanate. As expected, between pH 3 and 7 the rate of Schiff base protonation in the photocycle of the D96N mutant correlates with the concentrations of the acid forms of these agents. Dissection of the rates in the biexponential reprotonation kinetics of the Schiff base between pH 7 and 9 yielded calculated rate constants for the protonation equilibrium. Their dependencies on pH and azide or cyanate concentrations are consistent with both earlier suggested mechanisms: (i) azide and other weak acids may function as proton carriers in the protonation equilibrium of the Schiff base, or (ii) the binding of their anionic forms may catalyze proton conduction to and from the Schiff base. The measured rate constants allow the calculation of the pKa of the Schiff base during its reprotonation in the photocycle of D96N. It is 8.2-8.3, a value much below the pKa determined earlier in unphotolyzed bacteriorhodopsin.

On the nucleation and growth of amyloid beta-protein fibrils: detection of nuclei and quantitation of rate constants.

We have studied the fibrillogenesis of synthetic amyloid beta-protein-(1-40) fragment (A beta) in 0.1 M HCl. At low pH, A beta formed fibrils at a rate amenable to detailed monitoring by quasi-elastic light-scattering spectroscopy. Examination of the fibrils with circular dichroism spectroscopy and electron microscopy showed them to be highly similar to those found in amyloid plaques. We determined the hydrodynamic radii of A beta aggregates during the entire process of fibril nucleation and growth. Above an A beta concentration of approximately 0.1 mM, the initial rate of elongation and the final size of fibrils were independent of A beta concentration. Below an A beta concentration of 0.1 mM, the initial elongation rate was proportional to the peptide concentration, and the resulting fibrils were significantly longer than those formed at higher concentration. We also found that the surfactant n-dodecylhexaoxyethylene glycol monoether (C12E6) slowed nucleation and elongation of fibrils in a concentration-dependent manner. Our observations are consistent with a model of A beta fibrillogenesis that includes the following key steps: (i) peptide micelles form above a certain critical A beta concentration, (ii) fibrils nucleate within these micelles or on heterogeneous nuclei (seeds), and (iii) fibrils grow by irreversible binding of monomers to fibril ends. Interpretation of our data enabled us to determine the sizes of fibril nuclei and A beta micelles and the rates of fibril nucleation (from micelles) and fibril elongation. Our approach provides a powerful means for the quantitative assay of A beta fibrillogenesis.

The intrinsic ability of ribosomes to bind to endoplasmic reticulum membranes is regulated by signal recognition particle and nascent-polypeptide-associated complex.

Signal peptides direct the cotranslational targeting of nascent polypeptides to the endoplasmic reticulum (ER). It is currently believed that the signal recognition particle (SRP) mediates this targeting by first binding to signal peptides and then by directing the ribosome/nascent chain/SRP complex to the SRP receptor at the ER. We show that ribosomes can mediate targeting by directly binding to translocation sites. When purified away from cytosolic factors, including SRP and nascent-polypeptide-associated complex (NAC), in vitro assembled translation intermediates representing ribosome/nascent-chain complexes efficiently bound to microsomal membranes, and their nascent polypeptides could subsequently be efficiently translocated. Because removal of cytosolic factors from the ribosome/nascent-chain complexes also resulted in mistargeting of signalless nascent polypeptides, we previously investigated whether readdition of cytosolic factors, such as NAC and SRP, could restore fidelity to targeting. Without SRP, NAC prevented all nascent-chain-containing ribosomes from binding to the ER membrane. Furthermore, SRP prevented NAC from blocking ribosome-membrane association only when the nascent polypeptide contained a signal. Thus, NAC is a global ribosome-binding prevention factor regulated in activity by signal-peptide-directed SRP binding. A model presents ribosomes as the targeting vectors for delivering nascent polypeptides to translocation sites. In conjunction with signal peptides, SRP and NAC contribute to this specificity of ribosomal function by regulating exposure of a ribosomal membrane attachment site that binds to receptors in the ER membrane.

Transcription termination at intrinsic terminators: the role of the RNA hairpin.

Intrinsic termination of transcription in Escherichia coli involves the formation of an RNA hairpin in the nascent RNA. This hairpin plays a central role in the release of the transcript and polymerase at intrinsic termination sites on the DNA template. We have created variants of the lambda tR2 terminator hairpin and examined the relationship between the structure and stability of this hairpin and the template positions and efficiencies of termination. The results were used to test the simple nucleic acid destabilization model of Yager and von Hippel and showed that this model must be modified to provide a distinct role for the rU-rich sequence in the nascent RNA, since a perfect palindromic sequence that is sufficiently long to form an RNA hairpin that could destabilize the entire putative 12-bp RNA-DNA hybrid does not trigger termination at the expected positions. Rather, our results show that both a stable terminator hairpin and the run of 6-8 rU residues that immediately follows are required for effective intrinsic termination and that termination occurs at specific and invariant template positions relative to these two components. Possible structural or kinetic modifications of the simple model are proposed in the light of these findings and of recent results implicating "inchworming" and possible conformational heterogeneity of transcription complexes in intrinsic termination. Thus, these findings argue that the structure and dimensions of the hairpin are important determinants of the termination-elongation decision and suggest that a complete mechanism is likely to involve specific interactions of the polymerase, the RNA terminator hairpin, and, perhaps, the dT-rich template sequence that codes for the run of rU residues at the 3' end of the nascent transcript.

Intrinsic changes in developing retinal neurons result in regenerative failure of their axons.

The failure of mature mammalian central nervous system axons to regenerate after transection is usually attributed to influences of the extraneuronal milieu. Using explant cocultures of retina and midbrain tectum from hamsters, we have found evidence that these influences account for failure of regrowth of only a small minority of retinal axons. For most of the axons, there is a programmed loss of ability to elongate in the central nervous system. We show that there is a precipitous decline in the ability of retinal axons to reinnervate tectal targets when the retina is derived from pups on or after postnatal day 2, even when the target is embryonic. By contrast, embryonic retinal axons can regrow into tectum of any age, overcoming growth-inhibiting influences of glial factors.

Analysis of the interaction of ZAP-70 and syk protein-tyrosine kinases with the T-cell antigen receptor by plasmon resonance.

Tyrosine phosphorylation of a 17-amino acid immunoreceptor tyrosine-based activation motif (ITAM), conserved in each of the signaling subunits of the T-cell antigen receptor (TCR), mediates the recruitment of ZAP-70 and syk protein-tyrosine kinases (PTKs) to the activated receptor. The interaction between the two tandemly arranged Src-homology 2 (SH2) domains of this family of PTKs and each of the phosphotyrosine-containing ITAMs was examined by real-time measurements of kinetic parameters. The association rate and equilibrium binding constants for the ZAP-70 and syk SH2 domains were determined for the CD3 epsilon ITAM. Both PTKs bound with ka and Kd values of 5 x 10(6) M-1.sec-1 and approximately 25 nM, respectively. Bindings to the other TCR ITAMs (zeta 1, zeta 2, gamma, and delta ITAMs) were comparable, although the zeta 3 ITAM bound approximately 2.5-fold less well. Studies of the affinity of a single functional SH2 domain of ZAP-70 provided evidence for the cooperative nature of binding of the dual SH2 domains. Mutation of either single SH2 domain decreased the Kd by > 100-fold. Finally, the critical features of the ITAM for syk binding were found to be similar to those required for ZAP-70 binding. These data provide insight into the mechanism by which the multisubunit TCR interacts with downstream effector molecules.