Chromatin is the complex of genomic DNA and proteins that fills the nucleus of eukaryotic cells. The packaging of chromatin within the nucleus involves hierarchical folding at several levels. The first two levels produce nucleosomes and the 30-nm chromatin fiber. Additional levels give rise to the higher order organization of chromatin. Although the details of this organization remain unclear, their elucidation has recently been facilitated by chromosome conformation capture (3C) and related experimental methods, such as Hi-C, which detect physical interactions between different loci of genomic DNA. Knowledge of these interactions may be used to infer the conformation of the chromatin fiber. Motivated by this possibility, the present dissertation proposes an overall plan to develop and validate computational methods for the analysis of DNA interactions from Hi-C experiments. Within this plan, solutions for three specific problems are described. The first problem involves deducing an ensemble of chromatin conformations consistent with a given set of reference contact probabilities (CPs). The proposed solution consists of matching the given reference CPs with CPs estimated from simulations of a bead-chain polymer model representing the chromatin fiber. The second problem concerns accelerating the estimation of CPs from the simulated conformation ensembles. The proposed solution relies on fitting the extended generalized lambda distribution (EGLD) to distributions of simulated internal distances. The third problem involves quantifying the specificity of DNA cleavages in Hi-C experiments. The proposed solution consists of expressing an empirical local site distribution, obtained from Hi-C data, as a linear combination of conditional local site distributions, obtained from a reference genomic sequence, and then inferring cleavage fractions as weights of the linear combination. The above computational methods are promising steps toward the analysis of DNA interactions to investigate the higher order spatial organization of chromatin. This dissertation also includes a review on knots in biophysics, and describes a computational approach to identify efficient target sites for trans- splicing ribozymes

Recruiting endogenous adenosine deaminases using exogenous guide RNAs to edit cellular RNAs is a promising therapeutic strategy, but editing efficiency and durability remain low using current guide RNA designs. In this study, we engineered circular ADAR-recruiting guide RNAs (cadRNAs) to enable more efficient programmable adenosine-to-inosine RNA editing without requiring co-delivery of any exogenous proteins. Using these cadRNAs, we observed robust and durable RNA editing across multiple sites and cell lines, in both untranslated and coding regions of RNAs, and high transcriptome-wide specificity. Additionally, we increased transcript-level specificity for the target adenosine by incorporating interspersed loops in the antisense domains, reducing bystander editing. In vivo delivery of cadRNAs via adeno-associated viruses enabled 53% RNA editing of the mPCSK9 transcript in C57BL/6J mice livers and 12% UAG-to-UGG RNA correction of the amber nonsense mutation in the IDUA-W392X mouse model of mucopolysaccharidosis type I-Hurler syndrome. cadRNAs enable efficient programmable RNA editing in vivo with diverse protein modulation and gene therapeutic applications.

Adenosine deaminases acting on RNA (ADARs) can be repurposed to enable programmable RNA editing, however their exogenous delivery leads to transcriptome-wide off-targeting, and additionally, enzymatic activity on certain RNA motifs, especially those flanked by a 5' guanosine is very low thus limiting their utility as a transcriptome engineering toolset. Towards addressing these issues, we first performed a novel deep mutational scan of the ADAR2 deaminase domain, directly measuring the impact of every amino acid substitution across 261 residues, on RNA editing. This enabled us to create a domain-wide mutagenesis map while also revealing a novel hyperactive variant with improved enzymatic activity at 5'-GAN-3' motifs. As overexpression of ADAR enzymes, especially hyperactive variants, can lead to significant transcriptome-wide off-targeting, we next engineered a split-ADAR2 deaminase which resulted in >100-fold more specific RNA editing as compared to full-length deaminase overexpression. Taken together, we anticipate this systematic engineering of the ADAR2 deaminase domain will enable broader utility of the ADAR toolset for RNA biotechnology applications.

The eukaryotic genome is partitioned into distinct topological domains separated by boundary elements. Emerging data support the concept that several well-established nuclear compartments are ribonucleoprotein condensates assembled through the physical process of phase separation. Here, based on our demonstration that chemical disruption of nuclear condensate assembly weakens the insulation properties of a specific subset (∼20%) of topologically associated domain (TAD) boundaries, we report that the disrupted boundaries are characterized by a high level of transcription and striking spatial clustering. These topological boundary regions tend to be spatially associated, even interchromosomally, segregate with nuclear speckles, and harbor a specific subset of "housekeeping" genes widely expressed in diverse cell types. These observations reveal a previously unappreciated mode of genome organization mediated by conserved boundary elements harboring highly and widely expressed transcription units and associated transcriptional condensates.

We present in vivo sequence-specific RNA base editing via adenosine deaminases acting on RNA (ADAR) enzymes with associated ADAR guide RNAs (adRNAs). To achieve this, we systematically engineered adRNAs to harness ADARs, and comprehensively evaluated the specificity and activity of the toolsets in vitro and in vivo via two mouse models of human disease. We anticipate that this platform will enable tunable and reversible engineering of cellular RNAs for diverse applications.

Cell-cycle control is accomplished by cyclin-dependent kinases (CDKs), motivating extensive research into CDK targeting small-molecule drugs as cancer therapeutics. Here we use combinatorial CRISPR/Cas9 perturbations to uncover an extensive network of functional interdependencies among CDKs and related factors, identifying 43 synthetic-lethal and 12 synergistic interactions. We dissect CDK perturbations using single-cell RNAseq, for which we develop a novel computational framework to precisely quantify cell-cycle effects and diverse cell states orchestrated by specific CDKs. While pairwise disruption of CDK4/6 is synthetic-lethal, only CDK6 is required for normal cell-cycle progression and transcriptional activation. Multiple CDKs (CDK1/7/9/12) are synthetic-lethal in combination with PRMT5, independent of cell-cycle control. In-depth analysis of mRNA expression and splicing patterns provides multiple lines of evidence that the CDK-PRMT5 dependency is due to aberrant transcriptional regulation resulting in premature termination. These inter-dependencies translate to drug-drug synergies, with therapeutic implications in cancer and other diseases.

A crucial feature of differentiated cells is the rapid activation of enhancer-driven transcriptional programs in response to signals. The potential contributions of physicochemical properties of enhancer assembly in signaling events remain poorly understood. Here we report that in human breast cancer cells, the acute 17β-estradiol-dependent activation of functional enhancers requires assembly of an enhancer RNA-dependent ribonucleoprotein (eRNP) complex exhibiting properties of phase-separated condensates. Unexpectedly, while acute ligand-dependent assembly of eRNPs resulted in enhancer activation sensitive to chemical disruption of phase separation, chronically activated enhancers proved resistant to such disruption, with progressive maturation of eRNPs to a more gel-like state. Acute, but not chronic, stimulation resulted in ligand-induced, condensin-dependent changes in spatial chromatin conformation based on homotypic enhancer association, resulting in cooperative enhancer-activation events. Thus, distinct physicochemical properties of eRNP condensates on enhancers serve as determinants of rapid ligand-dependent alterations in chromosomal architecture and cooperative enhancer activation.