Studies on the role of histones in the structure and function of chromatin

Studies on the dissociation of histones from chromatin by increasing concentrations of sodium deoxycholate (DOC) have shown that histrone II is removed at lowest concentrations of DOC, while slightly higher concentrations remove histones III and IV. Still higher concentrations remove histone I.

The complete separation of chromatin and ¹⁴C-DOC by sucrose sedimentation indicated that the binding of DOC to chromatin is readily and completely reversible.

The dissociation of histones from chromatin by increasing concentrations of related cholanic acids and some of their conjugated derivatives were studied. The results suggested that the driving force for the interaction between the cholanic acid anion and histones is the lowering of the activity coefficient of the cholanic acid anion which occurs when it is partially removed from solution by interaction with hydrophobic regions of the positively charged histones.

The role of histones in the structure of chromatin has been studied by comparing the effects of selective removal of histones from chromatin by increasing concentrations of DOC with those caused by NaCl (removes histone I at lowest concentrations, while higher concentrations remove histones II, III, and IV). Properties studied included thermal denaturation, sedimentation velocity, flow dichroism, relaxation times of molecules oriented in a flow field, and the irreversible disruption of a 130 S, cross-linked component of sheared chromatin. The data indicated that none of the structural or chemical parameters with which these properties are correlated show a dependence on the presence of one particular histone fraction.

The template activity (ability to prime a 0.2 M KC1 DNA-dependent RNA synthesis system catalyzed by E. coli RNA polymerase) increases from that of native chromatin (approximately 25 per cent of that pure DNA) to that of pure DNA in a fashion which shows a nearly linear relationship to the amount of histone coverage of the template. The precipitability of partially dehistonized chromatin samples in 0.15 M NaCl shows a large dependence on the presence of histone I.

Functional genomic studies of the structure and regulation of eukaryotic transcriptomes

The main focus of this thesis is the use of high-throughput sequencing technologies in functional genomics (in particular in the form of ChIP-seq, chromatin immunoprecipitation coupled with sequencing, and RNA-seq) and the study of the structure and regulation of transcriptomes. Some parts of it are of a more methodological nature while others describe the application of these functional genomic tools to address various biological problems. A significant part of the research presented here was conducted as part of the ENCODE (ENCyclopedia Of DNA Elements) Project.

The first part of the thesis focuses on the structure and diversity of the human transcriptome. Chapter 1 contains an analysis of the diversity of the human polyadenylated transcriptome based on RNA-seq data generated for the ENCODE Project. Chapter 2 presents a simulation-based examination of the performance of some of the most popular computational tools used to assemble and quantify transcriptomes. Chapter 3 includes a study of variation in gene expression, alternative splicing and allelic expression bias on the single-cell level and on a genome-wide scale in human lymphoblastoid cells; it also brings forward a number of critical to the practice of single-cell RNA-seq measurements methodological considerations.

The second part presents several studies applying functional genomic tools to the study of the regulatory biology of organellar genomes, primarily in mammals but also in plants. Chapter 5 contains an analysis of the occupancy of the human mitochondrial genome by TFAM, an important structural and regulatory protein in mitochondria, using ChIP-seq. In Chapter 6, the mitochondrial DNA occupancy of the TFB2M transcriptional regulator, the MTERF termination factor, and the mitochondrial RNA and DNA polymerases is characterized. Chapter 7 consists of an investigation into the curious phenomenon of the physical association of nuclear transcription factors with mitochondrial DNA, based on the diverse collections of transcription factor ChIP-seq datasets generated by the ENCODE, mouseENCODE and modENCODE consortia. In Chapter 8 this line of research is further extended to existing publicly available ChIP-seq datasets in plants and their mitochondrial and plastid genomes.

The third part is dedicated to the analytical and experimental practice of ChIP-seq. As part of the ENCODE Project, a set of metrics for assessing the quality of ChIP-seq experiments was developed, and the results of this activity are presented in Chapter 9. These metrics were later used to carry out a global analysis of ChIP-seq quality in the published literature (Chapter 10). In Chapter 11, the development and initial application of an automated robotic ChIP-seq (in which these metrics also played a major role) is presented.

The fourth part presents the results of some additional projects the author has been involved in, including the study of the role of the Piwi protein in the transcriptional regulation of transposon expression in Drosophila (Chapter 12), and the use of single-cell RNA-seq to characterize the heterogeneity of gene expression during cellular reprogramming (Chapter 13).

The last part of the thesis provides a review of the results of the ENCODE Project and the interpretation of the complexity of the biochemical activity exhibited by mammalian genomes that they have revealed (Chapters 15 and 16), an overview of the expected in the near future technical developments and their impact on the field of functional genomics (Chapter 14), and a discussion of some so far insufficiently explored research areas, the future study of which will, in the opinion of the author, provide deep insights into many fundamental but not yet completely answered questions about the transcriptional biology of eukaryotes and its regulation.

Towards understanding the mixing characteristics of turbulent buoyant flows

This work proposes a new simulation methodology in which variable density turbulent flows can be studied in the context of a mixing layer with or without the presence of gravity. Specifically, this methodology is developed to probe the nature of non-buoyantly-driven (i.e. isotropically-driven) or buoyantly-driven mixing deep inside a mixing layer. Numerical forcing methods are incorporated into both the velocity and scalar fields, which extends the length of time over which mixing physics can be studied. The simulation framework is designed to allow for independent variation of four non-dimensional parameters, including the Reynolds, Richardson, Atwood, and Schmidt numbers. Additionally, the governing equations are integrated in such a way to allow for the relative magnitude of buoyant energy production and non-buoyant energy production to be varied.

The computational requirements needed to implement the proposed configuration are presented. They are justified in terms of grid resolution, order of accuracy, and transport scheme. Canonical features of turbulent buoyant flows are reproduced as validation of the proposed methodology. These features include the recovery of isotropic Kolmogorov scales under buoyant and non-buoyant conditions, the recovery of anisotropic one-dimensional energy spectra under buoyant conditions, and the preservation of known statistical distributions in the scalar field, as found in other DNS studies.

This simulation methodology is used to perform a parametric study of turbulent buoyant flows to discern the effects of varying the Reynolds, Richardson, and Atwood numbers on the resulting state of mixing. The effects of the Reynolds and Atwood numbers are isolated by looking at two energy dissipation rate conditions under non-buoyant (variable density) and constant density conditions. The effects of Richardson number are isolated by varying the ratio of buoyant energy production to total energy production from zero (non-buoyant) to one (entirely buoyant) under constant Atwood number, Schmidt number, and energy dissipation rate conditions. It is found that the major differences between non-buoyant and buoyant turbulent flows are contained in the transfer spectrum and longitudinal structure functions, while all other metrics are largely similar (e.g. energy spectra, alignment characteristics of the strain-rate tensor). Also, despite the differences noted between fully buoyant and non-buoyant turbulent fields, the scalar field, in all cases, is unchanged by these. The mixing dynamics in the scalar field are found to be insensitive to the source of turbulent kinetic energy production (non-buoyant vs. buoyant).

Velocity resolved - scalar modeled simulations of high Schmidt number turbulent transport

The objective of this thesis is to develop a framework to conduct velocity resolved - scalar modeled (VR-SM) simulations, which will enable accurate simulations at higher Reynolds and Schmidt (Sc) numbers than are currently feasible. The framework established will serve as a first step to enable future simulation studies for practical applications. To achieve this goal, in-depth analyses of the physical, numerical, and modeling aspects related to Sc>>1 are presented, specifically when modeling in the viscous-convective subrange. Transport characteristics are scrutinized by examining scalar-velocity Fourier mode interactions in Direct Numerical Simulation (DNS) datasets and suggest that scalar modes in the viscous-convective subrange do not directly affect large-scale transport for high Sc. Further observations confirm that discretization errors inherent in numerical schemes can be sufficiently large to wipe out any meaningful contribution from subfilter models. This provides strong incentive to develop more effective numerical schemes to support high Sc simulations. To lower numerical dissipation while maintaining physically and mathematically appropriate scalar bounds during the convection step, a novel method of enforcing bounds is formulated, specifically for use with cubic Hermite polynomials. Boundedness of the scalar being transported is effected by applying derivative limiting techniques, and physically plausible single sub-cell extrema are allowed to exist to help minimize numerical dissipation. The proposed bounding algorithm results in significant performance gain in DNS of turbulent mixing layers and of homogeneous isotropic turbulence. Next, the combined physical/mathematical behavior of the subfilter scalar-flux vector is analyzed in homogeneous isotropic turbulence, by examining vector orientation in the strain-rate eigenframe. The results indicate no discernible dependence on the modeled scalar field, and lead to the identification of the tensor-diffusivity model as a good representation of the subfilter flux. Velocity resolved - scalar modeled simulations of homogeneous isotropic turbulence are conducted to confirm the behavior theorized in these a priori analyses, and suggest that the tensor-diffusivity model is ideal for use in the viscous-convective subrange. Simulations of a turbulent mixing layer are also discussed, with the partial objective of analyzing Schmidt number dependence of a variety of scalar statistics. Large-scale statistics are confirmed to be relatively independent of the Schmidt number for Sc>>1, which is explained by the dominance of subfilter dissipation over resolved molecular dissipation in the simulations. Overall, the VR-SM framework presented is quite effective in predicting large-scale transport characteristics of high Schmidt number scalars, however, it is determined that prediction of subfilter quantities would entail additional modeling intended specifically for this purpose. The VR-SM simulations presented in this thesis provide us with the opportunity to overlap with experimental studies, while at the same time creating an assortment of baseline datasets for future validation of LES models, thereby satisfying the objectives outlined for this work.

In vivo analysis of interactions between trans-acting factors and their target genes

The investigations presented in this thesis use various in vivo techniques to understand how trans-acting factors control gene expression. The first part addresses the transcriptional regulation of muscle creatine kinase (MCK). MCK expression is activated during the course of development and is found only in differentiated muscle. Several in vivo footprints are observed at the enhancer of this gene, but all of these interactions are limited to cell types that express MCK. This is interesting because two of the footprints appear to represent muscle specific use of general transcription factors, while the other two correspond to sites that can bind the myogenic regulator, MyoD1, in vitro. MyoD1 and these general factors are present in myoblasts, but can bind to the enhancer only in myocytes. This suggests that either the factors themselves are post-translationally modified (phosphorylation or protein:protein interactions), or the accessibility of the enhancer to the factors is limited (changes in chromatin structure). The in vivo footprinting study of MCK was performed with a new ligation mediated, single-sided PCR (polymerase chain reaction) technique that I have developed.

The second half of the thesis concerns the regulation of mouse metallothionein (MT). Metallothioneins are a family of highly conserved housekeeping genes whose expression can be induced by heavy metals, steroids, and other stresses. By adapting a primer extension method of genomic sequencing to in vivo footprinting, I've observed both metal inducible and noninducible interactions at the promoter of MT-I. From these results I've been able to limit the possible mechanisms by which metal responsive trans-acting factors induce transcription. These interpretations correlate with a second line of experiments involving the stable titration of positive acting factors necessary for induction of MT. I've amplified the promoter of MT to 10^2-10^3 copies per cell by fusing the 5' and 3' ends of the MT gene to the coding region of DHFR and selecting cells for methotrexate resistance. In these cells, there is a metal-specific titration effect, and although it acts at the level of transcription, it appears to be independent of direct DNA binding factors.

DNA-mediated oxidation of transcription factor p53

Transcription factor p53 is the most commonly altered gene in human cancer. As a redox-active protein in direct contact with DNA, p53 can directly sense oxidative stress through DNA-mediated charge transport. Electron hole transport occurs with a shallow distance dependence over long distances through the π-stacked DNA bases, leading to the oxidation and dissociation of DNA-bound p53. The extent of p53 dissociation depends upon the redox potential of the response element DNA in direct contact with each p53 monomer. The DNA sequence dependence of p53 oxidative dissociation was examined by electrophoretic mobility shift assays using radiolabeled oligonucleotides containing both synthetic and human p53 response elements with an appended anthraquinone photooxidant. Greater p53 dissociation is observed from DNA sequences containing low redox potential purine regions, particularly guanine triplets, within the p53 response element. Using denaturing polyacrylamide gel electrophoresis of irradiated anthraquinone-modified DNA, the DNA damage sites, which correspond to locations of preferred electron hole localization, were determined. The resulting DNA damage preferentially localizes to guanine doublets and triplets within the response element. Oxidative DNA damage is inhibited in the presence of p53, however, only at DNA sites within the response element, and therefore in direct contact with p53. From these data, predictions about the sensitivity of human p53-binding sites to oxidative stress, as well as possible biological implications, have been made. On the basis of our data, the guanine pattern within the purine region of each p53-binding site determines the response of p53 to DNA-mediated oxidation, yielding for some sequences the oxidative dissociation of p53 from a distance and thereby providing another potential role for DNA charge transport chemistry within the cell.

To determine whether the change in p53 response element occupancy observed in vitro also correlates in cellulo, chromatin immunoprecipition (ChIP) and quantitative PCR (qPCR) were used to directly quantify p53 binding to certain response elements in HCT116N cells. The HCT116N cells containing a wild type p53 were treated with the photooxidant [Rh(phi)2bpy]³⁺, Nutlin-3 to upregulate p53, and subsequently irradiated to induce oxidative genomic stress. To covalently tether p53 interacting with DNA, the cells were fixed with disuccinimidyl glutarate and formaldehyde. The nuclei of the harvested cells were isolated, sonicated, and immunoprecipitated using magnetic beads conjugated with a monoclonal p53 antibody. The purified immounoprecipiated DNA was then quantified via qPCR and genomic sequencing. Overall, the ChIP results were significantly varied over ten experimental trials, but one trend is observed overall: greater variation of p53 occupancy is observed in response elements from which oxidative dissociation would be expected, while significantly less change in p53 occupancy occurs for response elements from which oxidative dissociation would not be anticipated.

The chemical oxidation of transcription factor p53 via DNA CT was also investigated with respect to the protein at the amino acid level. Transcription factor p53 plays a critical role in the cellular response to stress stimuli, which may be modulated through the redox modulation of conserved cysteine residues within the DNA-binding domain. Residues within p53 that enable oxidative dissociation are herein investigated. Of the 8 mutants studied by electrophoretic mobility shift assay (EMSA), only the C275S mutation significantly decreased the protein affinity (KD) for the Gadd45 response element. EMSA assays of p53 oxidative dissociation promoted by photoexcitation of anthraquinone-tethered Gadd45 oligonucleotides were used to determine the influence of p53 mutations on oxidative dissociation; mutation to C275S severely attenuates oxidative dissociation while C277S substantially attenuates dissociation. Differential thiol labeling was used to determine the oxidation states of cysteine residues within p53 after DNA-mediated oxidation. Reduced cysteines were iodoacetamide labeled, while oxidized cysteines participating in disulfide bonds were ¹³C₂D₂-iodoacetamide labeled. Intensities of respective iodoacetamide-modified peptide fragments were analyzed using a QTRAP 6500 LC-MS/MS system, quantified with Skyline, and directly compared. A distinct shift in peptide labeling toward ¹³C₂D₂-iodoacetamide labeled cysteines is observed in oxidized samples as compared to the respective controls. All of the observable cysteine residues trend toward the heavy label under conditions of DNA CT, indicating the formation of multiple disulfide bonds potentially among the C124, C135, C141, C182, C275, and C277. Based on these data it is proposed that disulfide formation involving C275 is critical for inducing oxidative dissociation of p53 from DNA.

Mechanisms of Transposable Element Repression by Piwi Proteins in the piRNA Pathway of Drosophila Germ Cells

The ability to reproduce is a defining characteristic of all living organisms. During reproduction, the integrity of genetic material transferred from one generation to the next is of utmost importance. Organisms have diverse strategies to ensure the fidelity of genomic information inherited between generations of individuals. In sexually reproducing animals, the piRNA pathway is an RNA-interference (RNAi) mechanism that protects the genomes of germ cells from the replication of ‘selfish’ genetic sequences called transposable elements (TE). When left unabated, the replication of TE sequences can cause gene disruption, double-stranded DNA breaks, and germ cell death that results in sterility of the organism. In Drosophila, the piRNA pathway is divided into a cytoplasmic and nuclear branch that involves the functions of three Piwi-clade Argonaute proteins—Piwi, Aubergine (Aub) and Argonaute-3 (Ago3)—which bind piwi-interacting RNA (piRNA) to form the effector complexes that represses deleterious TE sequences.

The work presented in this thesis examines the function and regulation of Piwi proteins in Drosophila germ cells. Chapter 1 presents an introduction to piRNA biogenesis and to the essential roles occupied by each Piwi protein in the repression of TE. We discuss the architecture and function of germ granules as the cellular compartments where much of the piRNA pathway operates. In Chapter 2, we present how Piwi in the nucleus co-transcriptionally targets genomic loci expressing TE sequences to direct the deposition of repressive chromatin marks. Chapter 3 examines the cytoplasmic function of the piRNA pathway, where we find that the protein Krimper coordinates Aub and Ago3 in the piRNA ping-pong pathway to adaptively target and destroy TE transcripts. Chapter 4 explores how interactions of Piwis with associated proteins are modulated by arginine methylation modifications. Lastly, in Chapter 5 I present evidence that the cytoplasmic branch of the piRNA pathway can potentially ‘cross-talk’ with the nuclear branch to transfer sequence information to better target and co-transcriptionally silence the genomic loci coding active TE sequences. Overall, the work presented in this thesis constitutes a part of the first steps in understanding the molecular mechanisms that protect germ cells from invasion by TE sequences.

Characterization and modeling of premixed turbulent n-heptane flames in the thin reaction zone regime

n-heptane/air premixed turbulent flames in the high-Karlovitz portion of the thin reaction zone regime are characterized and modeled in this thesis using Direct Numerical Simulations (DNS) with detailed chemistry. In order to perform these simulations, a time-integration scheme that can efficiently handle the stiffness of the equations solved is developed first. A first simulation with unity Lewis number is considered in order to assess the effect of turbulence on the flame in the absence of differential diffusion. A second simulation with non-unity Lewis numbers is considered to study how turbulence affects differential diffusion. In the absence of differential diffusion, minimal departure from the 1D unstretched flame structure (species vs. temperature profiles) is observed. In the non-unity Lewis number case, the flame structure lies between that of 1D unstretched flames with "laminar" non-unity Lewis numbers and unity Lewis number. This is attributed to effective Lewis numbers resulting from intense turbulent mixing and a first model is proposed. The reaction zone is shown to be thin for both flames, yet large chemical source term fluctuations are observed. The fuel consumption rate is found to be only weakly correlated with stretch, although local extinctions in the non-unity Lewis number case are well correlated with high curvature. These results explain the apparent turbulent flame speeds. Other variables that better correlate with this fuel burning rate are identified through a coordinate transformation. It is shown that the unity Lewis number turbulent flames can be accurately described by a set of 1D (in progress variable space) flamelet equations parameterized by the dissipation rate of the progress variable. In the non-unity Lewis number flames, the flamelet equations suggest a dependence on a second parameter, the diffusion of the progress variable. A new tabulation approach is proposed for the simulation of such flames with these dimensionally-reduced manifolds.

Investigations of Pyrrole-Imidazole Polyamide Effects on DNA Replication

Pyrrole–Imidazole polyamides are programmable, cell-permeable small molecules that bind in the minor groove of double-stranded DNA sequence-specifically. Polyamide binding has been shown to alter the local helical structure of DNA, disrupt protein-DNA interactions, and modulate endogenous gene expression. Py–Im polyamides targeted to the androgen receptor-DNA interface have been observed to decrease expression of androgen-regulated genes, upregulate p53, and induce apoptosis in a hormone-sensitive prostate cancer cell line. Here we report that androgen response element (ARE)-targeted polyamides induced DNA replication stress in a hormone-insensitive prostate cancer cell line. The ATR checkpoint kinase was activated in response to this stress, causing phosphorylation of MCM2, and FANCD2 was monoubiquitinated. Surprisingly, little single-stranded DNA was exhibited, and the ATR targets RPA2 and Chk1 were not phosphorylated. We conclude that polyamide induces relatively low level replication stress, and suggest inhibition of the replicative helicase as a putative mechanism based on in vitro assays. We also demonstrate polyamide-induced inhibition of DNA replication in cell free extracts from X. laevis oocytes. In this system, inhibition of chromatin decondensation is observed, preventing DNA replication initiation. Finally, we show that Py-Im polyamides targeted to the ARE and ETS binding sequence downregulate AR- and ERG-driven signaling in a prostate cancer cell line harboring the TMPRSS2-ERG fusion. In a mouse xenograft model, ARE-targeted polyamide treatment reduced growth of the tumor.

High resolution electron microscope studies of chromosomal fibers

Techniques are described for mounting and visualizing biological macromolecules for high resolution electron microscopy. Standard techniques are included in a discussion of new methods designed to provide the highest structural resolution. Methods are also discussed for handling samples on the grid, for making accurate size measurements at the 20 Å level, and for photographically enhancing image contrast.

The application of these techniques to the study of the binding of DNA polymerase to DNA is described. It is shown that the electron micrographs of this material are in agreement with the model proposed by Dr. Arthur Kornberg. A model is described which locates several active sites on the enzyme.

The chromosomal material of the protozoan tetrahymena has been isolated and characterized by biochemical techniques and by electron microscopy. This material is shown to be typical of chromatin of higher creatures.

Comparison with other chromatins discloses that the genome of tetrahymena is highly template active and has a relatively simple genetic construction.

High resolution electron microscope procedures developed in this work have been combined with standard biochemical techniques to give a comprehensive picture of the structure of interphase chromosome fibers. The distribution of the chromosomal proteins along its DNA is discussed.

The interactions of basic proteins and DNA

I. ELECTROPHORESIS OF THE NUCLEIC ACIDS

A zone electrophoresis apparatus using ultraviolet optics has been constructed to study nucleic acids at concentrations less than 0.004%. Native DNA has a mobility about 15% higher than denatured DNA over a range of conditions. Otherwise, the electrophoretic mobility is independent of molecular weight, base composition or source. DNA mobilities change in the expected way with pH but the fractional change in mobility is less than the calculated change in charge. A small decrease in mobility accompanies an increase in ionic strength. RNA’s from various sources have mobilities slightly lower than denatured DNA except for s-RNA which travels slightly faster. The important considerations governing the mobility of nucleic acids appear to be the nature of the hydrodynamic segment, and the binding of counterions. The differences between electrophoresis and sedimentation stem from the fact that all random coil polyelectrolytes are fundamentally free draining in electrophoresis.

II. THE CYTOCHROME C/DNA COMPLEX

The basic protein, cytochrome c, has been complexed to DNA. Up to a cytochrome:DNA mass ratio of 2, a single type of complex is formed. Dissociation of this complex occurs between 0.05F and 0.1F NaCl. The complexing of cytochrome to DNA causes a slight increase in the melting temperature of the DNA, and a reduction of the electrophoretic mobility proportional to the decrease in net charge. Above a cytochrome:DNA mass ratio of 2.5, a different type of complex is formed. The results suggest that complexes such as are formed in the Kleinschmidt technique of electron microscopy would not exist in bulk solution and are exclusively film phenomena.

III. STUDIES OF THE ELECTROPHORESIS AND MELTING BEHAVIOUR OF NUCLEOHISTONES

Electrophoresis studies on reconstituted nucleohistones indicate that the electrophoretic mobility for these complexes is a function of the net charge of the complex. The mobility is therefore dependent on the charge density of the histone complexing the DNA, as well as on the histone/DNA ratio. It is found that the different histones affect the transition from native to denatured DNA in different ways. It appears that histone I is exchanging quite rapidly between DNA molecules in 0.01 F salt, while histone II is irreversibly bound. Histone III-IV enhances the capacity of non-strand separated denatured DNA to reanneal. Studies on native nucleoproteins indicate that there are no gene-sized uncomplexed DNA regions in any preparations studied.

IV. THE DISSOCIATION OF HISTONE FROM CALF THYMUS CROMATIN

Calf thymus nucleoprotein was treated with varying concentrations of NaCl. The identity of the histones associated and dissociated from the DNA at each salt concentration was determined by gel electrophoresis. It was found that there is no appreciable histone dissociation below 0.4 F NaCl. The lysine rich histones dissociate between 0.4 and 0.5 F NaCl. Their dissociation is accompanies by a marked increase in the solubility of the chromatin. The moderately lysine rich histones dissociate mainly between 0.8 and 1.1 F NaCl. There are two arginine rich histone components: the first dissociates between 0.8 F and 1.1 F NaCl, but the second class is the very last to be dissociated from the DNA (dissociation beginning at 1.0 F NaCl). By 2.0 F NaCl, essentially all the histones are dissociated.

The properties of the extracted nucleoprotein were studied. The electrophoretic mobility increases and the melting temperature decreases as more histones are dissociated from the DNA. A comparison with the dissociation of histones from DNA in NaClO₄ shows that to dissociate the same class of histones, the concentration of NaCl required is twice that of NaClO₄.

Small scale turbulence in high Karlovitz number premixed flames

The purpose of this thesis is to characterize the behavior of the smallest turbulent scales in high Karlovitz number (Ka) premixed flames. These scales are particularly important in the two-way coupling between turbulence and chemistry and better understanding of these scales will support future modeling efforts using large eddy simulations (LES). The smallest turbulent scales are studied by considering the vorticity vector, ω, and its transport equation.

Due to the complexity of turbulent combustion introduced by the wide range of length and time scales, the two-dimensional vortex-flame interaction is first studied as a simplified test case. Numerical and analytical techniques are used to discern the dominate transport terms and their effects on vorticity based on the initial size and strength of the vortex. This description of the effects of the flame on a vortex provides a foundation for investigating vorticity in turbulent combustion.

Subsequently, enstrophy, ω² = ω • ω, and its transport equation are investigated in premixed turbulent combustion. For this purpose, a series of direct numerical simulations (DNS) of premixed n-heptane/air flames are performed, the conditions of which span a wide range of unburnt Karlovitz numbers and turbulent Reynolds numbers. Theoretical scaling analysis along with the DNS results support that, at high Karlovitz number, enstrophy transport is controlled by the viscous dissipation and vortex stretching/production terms. As a result, vorticity scales throughout the flame with the inverse of the Kolmogorov time scale, τ_η, just as in homogeneous isotropic turbulence. As τ_η is only a function of the viscosity and dissipation rate, this supports the validity of Kolmogorov’s first similarity hypothesis for sufficiently high Ka numbers (Ka ≳ 100). These conclusions are in contrast to low Karlovitz number behavior, where dilatation and baroclinic torque have a significant impact on vorticity within the flame. Results are unaffected by the transport model, chemical model, turbulent Reynolds number, and lastly the physical configuration.

Next, the isotropy of vorticity is assessed. It is found that given a sufficiently large value of the Karlovitz number (Ka ≳ 100) the vorticity is isotropic. At lower Karlovitz numbers, anisotropy develops due to the effects of the flame on the vortex stretching/production term. In this case, the local dynamics of vorticity in the strain-rate tensor, S, eigenframe are altered by the flame. At sufficiently high Karlovitz numbers, the dynamics of vorticity in this eigenframe resemble that of homogeneous isotropic turbulence.

Combined, the results of this thesis support that both the magnitude and orientation of vorticity resemble the behavior of homogeneous isotropic turbulence, given a sufficiently high Karlovitz number (Ka ≳ 100). This supports the validity of Kolmogorov’s first similarity hypothesis and the hypothesis of local isotropy under these condition. However, dramatically different behavior is found at lower Karlovitz numbers. These conclusions provides/suggests directions for modeling high Karlovitz number premixed flames using LES. With more accurate models, the design of aircraft combustors and other combustion based devices may better mitigate the detrimental effects of combustion, from reducing CO₂ and soot production to increasing engine efficiency.

Investigations of nonhistone chromosomal proteins

The major nonhistone chromosomal proteins (NHC proteins) are a group of 14-20 acidic proteins associated with DNA in eukaryotic chromatin. In comparisons by SDS gel electrophoresis (molecular weight sieving) one observes a high degree of homology among the NHC protein fractions of different tissues from a given species. Tissue-specific protein bands are also observed. The appearance of a new NHC protein, A, in the NHC proteins of rat liver stimulated to divide by partial hepatectomy and of rat ascites cells suggests that this protein may play a role in preparing the cell for division. The NHC proteins of the same tissue from different species are also very similar. Quantitative but not qualitative changes in the NHC proteins of rat uterus are observed on stimulation (in vivo) with estrogen. These observations suggest that the major NHC proteins play a general role in chromatin structure and the regulation of genome expression; several may be enzymes of nucleic acid and histone metabolism and/or structural proteins analogous to histones. One such enzyme, a protease which readily and preferentially degrades histones, can be extracted from chromatin with 0.7 N NaCl.

Although the NHC proteins readily aggregate, they can be separated from histone and fractionated by ion exchange chromatography on Sephadex SE C-25 resin in 10 M urea-25% formic acid (pH 2.5). Following further purification, four fractions of NHC protein are obtained; two of these are single purified proteins, and the other two contain 4-6 and 4-7 different proteins. These NHC proteins show a ratio of acidic to basic amino acids from 2.7 to 1.2 and isoelectric points from apparently less than 3.7 to 8.0. These isolated fractions appear more soluble and easier to work with than any whole NHC protein preparation.