Sequence-function relationships in E. coli transcriptional regulation

Understanding how transcriptional regulatory sequence maps to regulatory function remains a difficult problem in regulatory biology. Given a particular DNA sequence for a bacterial promoter region, we would like to be able to say which transcription factors bind there, how strongly they bind, and whether they interact with each other and/or RNA polymerase, with the ultimate objective of integrating knowledge of these parameters into a prediction of gene expression levels. The theoretical framework of statistical thermodynamics provides a useful framework for doing so, enabling us to predict how gene expression levels depend on transcription factor binding energies and concentrations. We used thermodynamic models, coupled with models of the sequence-dependent binding energies of transcription factors and RNAP, to construct a genotype to phenotype map for the level of repression exhibited by the lac promoter, and tested it experimentally using a set of promoter variants from E. coli strains isolated from different natural environments. For this work, we sought to ``reverse engineer'' naturally occurring promoter sequences to understand how variations in promoter sequence affects gene expression. The natural inverse of this approach is to ``forward engineer'' promoter sequences to obtain targeted levels of gene expression. We used a high precision model of RNAP-DNA sequence dependent binding energy, coupled with a thermodynamic model relating binding energy to gene expression, to predictively design and verify a suite of synthetic E. coli promoters whose expression varied over nearly three orders of magnitude.

However, although thermodynamic models enable predictions of mean levels of gene expression, it has become evident that cell-to-cell variability or ``noise'' in gene expression can also play a biologically important role. In order to address this aspect of gene regulation, we developed models based on the chemical master equation framework and used them to explore the noise properties of a number of common E. coli regulatory motifs; these properties included the dependence of the noise on parameters such as transcription factor binding strength and copy number. We then performed experiments in which these parameters were systematically varied and measured the level of variability using mRNA FISH. The results showed a clear dependence of the noise on these parameters, in accord with model predictions.

Finally, one shortcoming of the preceding modeling frameworks is that their applicability is largely limited to systems that are already well-characterized, such as the lac promoter. Motivated by this fact, we used a high throughput promoter mutagenesis assay called Sort-Seq to explore the completely uncharacterized transcriptional regulatory DNA of the E. coli mechanosensitive channel of large conductance (MscL). We identified several candidate transcription factor binding sites, and work is continuing to identify the associated proteins.

Regulation of the fitness of Atlantic salmon (Salmo salar) by intra-specific competition amongst the juveniles

Factors affecting the fitness of juvenile salmon are discussed. Although fitness from the genetic point of view is defined as the relative capacity of carriers of a given genotype to transmit their genes to the gene pool of the following generations, growth and survival of individuals are also components of fitness, and are influenced by responses to competition, which is the major topic of this article including implications for management. In order to better understand the relationships of density-dependent survival in Newfoundland, egg depositions were manipulated experimentally in the Freshwater River. Figures demonstrate the relationship between stock (number of eggs per 100 m2 of river) and recruitment (number of smolts per l00 m2 of Atlantic salmon, and also the percentage survival from egg to smolt stage related to potential egg depositions.

Molecular genetic studies on voltage-gated ion channels

Several different methods have been employed in the study of voltage-gated ion channels. Electrophysiological studies on excitable cells in vertebrates and molluscs have shown that many different voltage-gated potassium (K⁺) channels and sodium channels may coexist in the same organism. Parallel genetic studies in Drosophila have identified mutations in several genes that alter the properties of specific subsets of physiologically identified ion channels. Chapter 2 describes molecular studies that identify two Drosophila homologs of vertebrate sodium-channel genes. Mutations in one of these Drosophila sodium-channel genes are shown to be responsible for the temperature-dependent paralysis of a behavioural mutant para^ts. Evolutionary arguments, based on the partial sequences of the two Drosophila genes, suggest that subfamilies of voltage-gated sodium channels in vertebrates remain to be identified.

In Drosophila, diverse voltage-gated K⁺ channels arise from alternatively spliced mRNAs generated at the Shaker locus. Chapter 3 and the Appendices describe the isolation and characterization of several human K⁺-channel genes, similar in sequence to Shaker. Each of these human genes has a highly conserved homolog in rodents; thus, this K⁺-channel gene family probably diversified prior to the mammalian radiation. Functional K⁺ channels encoded by these genes have been expressed in Xenopus oocytes and their properties have been analyzed by electrophysiological methods. These studies demonstrate that both transient and noninactivating voltage-gated K⁺ channels may be encoded by mammalian genes closely related to Shaker. In addition, results presented in Appendix 3 clearly demonstrate that independent gene products from two K⁺-channel genes may efficiently co-assemble into heterooligomeric K⁺ channels with properties distinct from either homomultimeric channel. This finding suggests yet another molecular mechanism for the generation of K⁺-channel diversity.

The Mechanical Genome in Regulation and Infection

Biological information storage and retrieval is a dynamic process that requires the genome to undergo dramatic structural rearrangements. Recent advances in single-molecule techniques have allowed precise quantification of the nano-mechanical properties of DNA [1, 2], and direct in vivo observation of molecules in action [3]. In this work, we will examine elasticity in protein-mediated DNA looping, whose structural rearrangement is essential for transcriptional regulation in both prokaryotes and eukaryotes. We will look at hydrodynamics in the process of viral DNA ejection, which mediates information transfer and exchange and has prominent implications in evolution. As in the case of Kepler's laws of planetary motion leading to Newton's gravitational theory, and the allometric scaling laws in biology revealing the organizing principles of complex networks [4], experimental data collapse in these biological phenomena has guided much of our studies and urged us to find the underlying physical principles.

Investigation of competitive antagonist binding to the nicotinic acetylcholine receptor using voltage-jump and light-flash techniques

1. The effect of 2,2’-bis-[α-(trimethylammonium)methyl]azobenzene (2BQ), a photoisomerizable competitive antagonist, was studied at the nicotinic acetycholine receptor of Electrophorus electroplaques using voltage-jump and light-flash techniques.

2. 2BQ, at concentrations below 3 μΜ, reduced the amplitude of voltage-jump relaxations but had little effect on the voltage-jump relaxation time constants under all experimental conditions. At higher concentrations and voltages more negative than -150 mV, 2BQ caused significant open channel blockade.

3. Dose-ratio studies showed that the cis and trans isomers of 2BQ have equilibrium binding constants (K _ᵢ) of .33 and 1.0 μΜ, respectively. The binding constants determined for both isomers are independent of temperature, voltage, agonist concentration, and the nature of the agonist.

4. In a solution of predominantly cis-2BQ, visible-light flashes led to a net cis→trans isomerization and caused an increase in the agonist-induced current. This increase had at least two exponential components; the larger amplitude component had the same time constant as a subsequent voltage-jump relaxation; the smaller amplitude component was investigated using ultraviolet light flashes.

5. In a solution of predominantly trans-2BQ, UV-light flashes led to a net trans→cis isomerization and caused a net decrease in the agonist-induced current. This effect had at least two exponential components. The smaller and faster component was an increase in agonist-induced current and had a similar time constant to the voltage-jump relaxation. The larger component was a slow decrease in the agonist-induced current with rate constant approximately an order of magnitude less than that of the voltage-jump relaxation. This slow component provided a measure of the rate constant for dissociation of cis-2BQ (k_ = 60/s at 20°C). Simple modelling of the slope of the dose-rate curves yields an association rate constant of 1.6 x 10⁸/M/s. This agrees with the association rate constant of 1.8 x 10⁸/M/s estimated from the binding constant (K_i). The Q₁₀ of the dissociation rate constant of cis-2BQ was 3.3 between 6° and 20°C. The rate constants for association and dissociation of cis-28Q at receptors are independent of voltage, agonist concentration, and the nature of the agonist.

6. We have measured the molecular rate constants of a competitive antagonist which has roughly the same K _ᵢ as d-tubocurarine but interacts more slowly with the receptor. This leads to the conclusion that curare itself has an association rate constant of 4 x 10⁹/M/s or roughly as fast as possible for an encounter-limited reaction.

The Regulation of Sleep and Circadian Rhythms: The Role of Melatonin and Adenosine in Zebrafish

Sleep is a highly conserved behavioral state whose regulation is still unclear. In this thesis I initially briefly introduce the known sleep circuitry and regulation in vertebrates, and why zebrafish is seen as a good model to study sleep-regulation. I describe the existing two-process model of sleep regulation, which posits that the two processes C (circadian) and S (homeostatic) control timing of sleep-wake behavior. I then study the role melatonin plays in the circadian regulation of sleep using zebrafish. Firstly, we find that the absence of melatonin results in a reduction of sleep at night, establishing that endogenous melatonin is required for sleep at night. Secondly, melatonin mutants show a reduction in sleep in animals with no functional behavioral rhythms suggesting that melatonin does not require intact circadian rhythms for its effect on sleep. Thirdly, melatonin mutants do not exhibit any changes in circadian rhythms, suggesting that the circadian clock does not require melatonin for its function. Fourthly, we find that in the absence of melatonin, there is no rhythmic expression of sleep, suggesting that melatonin is the output molecule of process C. Lastly, we describe a connection between adenosine signaling (output molecules of process S), and melatonin. Following this we proceed to study the role adenosine signaling plays in sleep-wake behavior. We find that firstly, adenosine receptor A1 and A2 are involved in sleep- wake behavior in zebrafish, based on agonist/antagonist behavioral results. Secondly, we find that several brain regions such as PACAP cells in the rostral midbrain, GABAergic cells in the forebrain and hindbrain, Dopamine and serotonin cells in the caudal hypothalamus and sox2 cells lining the hindbrain ventricle are activated in response to the A1 antagonist and VMAT positive cells are activated in response to the A2A agonist, suggesting these areas are involved in adenosine signaling in zebrafish. Thirdly, we find that knocking out the zebrafish adenosine receptors has no effect on sleep architecture. Lastly, we find that while the A1 agonist phenotype requires the zfAdora1a receptor, the antagonist and the A2A agonist behavioral phenotypes are not mediated by the zfAdora1a, zfAdora1b and zfAdoraA2Aa, zfAdora2Ab receptors respectively.

Part I. Regulation of development in the cellular slime mold Dictyostelium discoideum. Part II. Polysomes and RNA synthesis during early development of the surf clam Spisula solidissima

Part I. The cellular slime mold Dictyostelium discoideum is a simple eukaryote which undergoes a multi-cellular developmental process. Single cell myxamoebae divide vegetatively in the presence of a food source. When the food is depleted or removed, the cells aggregate, forming a migrating pseudoplasmodium which differentiates into a fruiting body containing stalk and spore cells. I have shown that during the developmental cycle glycogen phosphorylase, aminopeptidase, and alanine transaminase are developmentally regulated, that is their specific activities increased at a specific time in the developmental cycle. Phosphorylase activity is undetectable in developing cells until mid-aggregation whereupon it increases and reaches a maximum at mid-culmination. Thereafter the enzyme disappears. Actinomycin D and cycloheximide studies as well as studies with morphologically aberrant and temporally deranged mutants indicate that prior RNA and concomitant protein synthesis are necessary for the rise and decrease in activity and support the view that the appearance of the enzyme is regulated at the transcriptional level. Aminopeptidase and alanine transaminase increase 3 fold starting at starvation and reach maximum activity at 18 and 5 hours respectively.

The cellular DNA s of D. discoideum were characterized by CsC1 buoyant density gradient centrifugation and by renaturation kinetics. Whole cell DNA exhibits three bands in CsCl: ρ = 1.676 g/cc (nuclear main band), 1.687 (nuclear satellite), and 1.682 (mitochondrial). Reassociation kinetics at a criterion of Tm -23°C indicates that the nuclear reiterated sequences make up 30% of the genome (Cot_1/2 (pure) 0.28) and the single-copy DNA 70% (Cot_1/2(pure) 70). The complexity of the nuclear genome is 30 x 10⁹ daltons and that of the mitochondrial DNA is 35-40 x 10⁶ daltons (Cot_1/2 0.15). rRNA cistrons constitute 2.2% of nuclear DNA and have a ρ = 1.682.

RNA extracted from 4 stages during developmental cycle of Dictyostelium was hybridized with purified single-copy nuclear DNA. The hybrids had properties indicative of single-copy DNA-RNA hybrids. These studies indicate that there are, during development, qualitative and quantitative changes in the portion of the single-copy of the genome transcribed. Overall, 56% of the genome is represented by transcripts between the amoeba and mid-culmination stages. Some 19% are sequences which are represented at all stages while 37% of the genome consists of stage specific sequences.

Part II. RNA and protein synthesis and polysome formation were studied during early development of the surf clam Spisula solidissima embryos. The oocyte has a small number of polysomes and a low but measurable rate of protein synthesis (leucine-³H incorporation). After fertilization, there is a continual increase in the percentage of ribosomes sedimenting in the polysome region. Newly synthesized RNA (uridine-5-³H incorporation) was found in polysomes as early as the 2-cell stage. During cleavage, the newly formed RNA is associated mainly with the light polysomes.

RNA extracted from polysomes labeled at the 4-cell stage is polydisperse, nonribosomal, and non-4 S. Actinomycin D causes a reduction of about 30% of the polysomes formed between fertilization and the 16-cell stage.

In the early cleavage stages the light polysomes are mostly affected by actinomycin.