STUDIES INTO THE MOLECULAR MECHANISMS REGULATING MUSCULAR ATROPHY RESULTING FROM HINDLIMB IMMOBILIZATION IN THE RAT (PROTEIN, SYNTHESIS, RNA)

The purpose of the work performed in this dissertation was to examine some of the possible regulatory mechanisms involved in the initiation of muscular atrophy during periods of decreased muscle utilization resulting from hindlimb immobilization in the rat. A 37% decrease in the rate of total muscle protein synthesis which has been observed to occur in the first 6 h of immobilization contributes significantly to the observed loss of protein during immobilization.^ The rates of cytochrome c and actin synthesis were determined in adult rat red vastus lateralis and gastrocnemius muscles, respectively, by the constant infusion and incorporation of ('3)H-tyrosine into protein. The fractional synthesis rates of both actin and cytochrome c were significantly decreased (P < 0.05) in the 6th h of hindlimb immobilization.^ RHA was extracted from adult rat gastrocnemius muscle by modification of the phenol: chloroform: SDS extraction procedures commonly used for preparation of RNA for hybridization analysis from other mammalian tissues. RNA content of rat gastrocnemius muscle, as determined by this method of extraction and its subsequent quantification by UV absorbance and orcinol assay, was significantly greater than the RNA content previously determined for adult rat gastrocnemius by other commonly employed methods.^ RNA extracted by this method from gastrocnemius muscles of control and 6h immobilized rats was subjected to "dot blot" hybridization to ('32)P-labelled probe from plasmid p749, containing a cDNA sequence complementary to (alpha)-actin mRNA and from rat skeletal muscle. (alpha)-Actin specific mRNA content as estimated by this procedure is not significantly decreased in rat gastrocnemius following 6h or hindlimb immobilization. However, (alpha)-actin specific mRNA content is significantly decreased (P < 0.05) in adult rat gastrocnemius (alpha)-actin specific mRNA is not decreased in adult rat gastrocnemius muscle following 6h of immobilization, a time when actin synthesis is significantly decreased, it is concluded that a change in (alpha)-actin specific mRNA content is not the initiating event responsible for the early decrease in actin synthesis observed in the 6th h of immobilization. ^

Mapping and characterization of the Xlsirt P11 RNA cis-acting localization element

In many organisms, polarity of the oocyte is established post-transcriptionally via subcellular RNA localization. Many RNAs are localized during oogenesis in Xenopus laevis, including Xlsirts ( Xenopus laevis short interspersed repeat transcripts) [Kloc, 1993]. Xlsirts constitute a large family defined by highly homologous repeat units 79–81 nucleotides in length. Endogenous Xlsirt RNAs use the METRO (Message Transport Organizer) pathway of localization, where RNAs are transported from the nucleus to the mitochondrial cloud in stage I oocytes. Secondly, RNAs anchor at the vegetal pole in stage II oocytes. Exogenous Xlsirt RNAs can also utilize the Late pathway of localization, which involves localization to the vegetal cortex during stage III of oogenesis and results in RNAs anchored in the cortex of the entire vegetal hemisphere. ^ The Xlsirts localization signal is contained within the repeat region. This study was designed to test the hypothesis that there are cis -acting localization elements in Xlsirts, and that higher order structure plays a role. Results of experiments on Xlsirt P11, a 1700 basepair (bp) family member, led to the conclusion that a 137-bp fragment of the repetitive region is necessary and sufficient for METRO and Late pathway localization. This analysis definitively demonstrates that the Xlsirt localization signal for the METRO and Late pathways reside within the repetitive region and not within the flanking regions. Analysis of Xlsirt linker scanning mutations revealed two METRO-pathway specific subelements, and one Late-pathway specific subelement. Functional, computer, and biochemical evidence relates the higher order structure of this element to its ability to function as a localization element. ^ Xlsirt 137 is 99% identical to the Xlsirt consensus sequence identified in this study, suggesting that it is the localization element for all localized Xlsirt family members. The repeat unit was reframed based on function, rather than arbitrarily based on sequence. This work supports the hypothesis presented in 1981 by George Spohr, who originally isolated the Xlsirts, which stated that the highly conserved repetitive elements must be constrained from variability due to some unknown function of the repeats themselves. These studies shed light on the mechanism of RNA localization, linking structure and function. ^

RNA binding characteristics and overall topology of the vaccinia virus polyadenylation complex

mRNA 3′ polyadenylation is central to mRNA biogenesis in prokaryotes and eukaryotes, and is implicated in numerous aspects of mRNA metabolism, including efficiency of mRNA export from the nucleus, message stability, and initiation of translation. However, due to the great complexity of the eukaryotic polyadenylation apparatus, the mechanisms of RNA 3 ′ end processing have remained elusive. Although the RNA processing reactions leading to polyadenylated messenger RNA have been studied in many systems, and much progress has been made, a complete understanding of the biochemistry of the poly(A) polymerase enzyme is still lacking. My research uses Vaccinia virus as a model system to gain a better understanding of this complicated polyadenylation process, which consist of RNA binding, catalysis and polymerase translocation. ^ Vaccinia virus replicates in the cytoplasm of its host cell, so it must employ its own poly(A) polymerase (PAP), a heterodimer of two virus encoded proteins, VP55 and VP39. VP55 is the catalytic subunit, adding 30 adenylates to a non-polyadenylated RNA in a rapid processive manner before abruptly changing to a slow, non-processive mode of adenylate addition and dissociating from the RNA. VP39 is the stimulatory subunit. It has no polyadenylation catalytic activity by itself, but when associated with VP55 it facilitates the semi-processive synthesis of tails several hundred adenylates in length. ^ Oligonucleotide selection and competition studies have shown that the heterodimer binds a minimal motif of (rU)2 (N)25 U, the “heterodimer binding motif”, within an oligonucleotide, and its primer selection for polyadenylation is base-type specific. ^ Crosslinking studies using photosensitive uridylate analogs show that within a VP55-VP39-primer ternary complex, VP55 comes into contact with all three required uridylates, while VP39 only contacts the downstream uridylate. Further studies, using a backbone-anchored photosensitive crosslinker show that both PAP subunits are in close proximity to the downstream −10 to −21 region of 50mer model primers containing the heterodimer binding motif. This equal crosslinking to both subunits suggests that the dimerization of VP55 and VP39 creates either a cleft or a channel between the two subunits through which this region of RNA passes. ^ Peptide mapping studies of VP39 covalently crosslinked to the oligonucleotide have identified residue R107 as the amino acid in close proximity to the −10 uridylate. This helps us project a conceptual model onto the known physical surface of this subunit. In the absence of any tertiary structural data for VP55, we have used a series of oligonucleotide selection assays, as well as crosslinking, nucleotide transfer assays, and gel shift assays to gain insight into the requirements for binding, polyadenylation and translocation. Collectively, these data allow us to put together a comprehensive model of the structure and function of the polyadenylation ternary complex consisting of VP39, VP55 and RNA. ^

Correction of aberrant FGFR1 RNA splicing through blocking intronic regulatory elements

Analysis of the human genome has revealed that more than 74% of human genes undergo alternative RNA splicing. Aberrations in alternative RNA splicing have been associated with several human disorders, including cancer. ^ We studied the aberrant expression of alternative RNA splicing isoforms of the Fibroblast Growth Factor Receptor 1 (FGFR1) gene in a human glioblastoma cancer model. Normal glial cells express the FGFR1α, which contains three extracellular domains. In tumors the most abundant isoform is the FGFR1β, which lacks the first extracellular domain due to the skipping of a single exon, termed alpha. The skipping of the α-exon is regulated by two intronic silencing sequences within the precursor mRNA. Since we observed no mutations on these elements in tumor cells, we hypothesized that the over-expression of regulatory proteins that recognize these sequences is responsible for the aberrant expression of splicing isoforms. Hence, we blocked the formation of protein complexes on the ISS using antisense RNA oligonucleotides in vitro. We also evaluated the impact of the ISS antisense oligonucleotides on the endogenous FGFR1 splicing, in a glioblastoma cell model. By targeting intronic regulatory elements we were able to increase the level of alpha exon inclusion up to 90% in glioblastoma cells. The effect was dose dependent, sequence specific and reproducible in glioblastoma and other cancer cells, which also exhibit an alpha exon skipping phenotype. Targeting FGFR1 endogenous ISS1 and ISS2 sequences did not have an additive or synergistic effect, which suggest a regulatory splicing mechanism that requires the interaction of complexes formed on these elements. An increase in the levels of the FGFR1α isoform resulted in a reduction in cell invasiveness. Also, a significant increase in the levels of caspase 3/7 activities, which is indicative of an elevation in apoptosis levels, suggests that expression of FGFR1β might be relevant for tumor survival. These studies demonstrate that it is possible to prevent aberrant expression of exon skipping events through the targeting of intronic regulatory elements, providing an important new therapeutic tool for the correction of human disease caused by alternative RNA splicing. ^

RNA polymerase II large subunit degradation induced by ecteinascidin 743: Molecular characterization and subsequent rational investigation of antitumor mechanisms in cancers with clinical response

Ecteinascidin 743 (Et-743), which is a novel DNA minor groove alkylator with a unique spectrum of antitumor activity, is currently being evaluated in phase II/III clinical trials. Although the precise molecular mechanisms responsible for the observed antitumor activity are poorly understood, recent data suggests that post-translational modifications of RNA polymerase II Large Subunit (RNAPII LS) may play a central role in the cellular response to this promising anticancer agent. The stalling of an actively transcribing RNAPII LS at Et-743-DNA adducts is the initial cellular signal for transcription-coupled nucleotide excision repair (TC-NER). In this manner, Et-743 poisons TC-NER and produces DNA single strand breaks. Et-743 also inhibits the transcription and RNAPII LS-mediated expression of selected genes. Because the poisoning of TC-NER and transcription inhibition are critical components of the molecular response to Et-743 treatment, we have investigated if changes in RNAPII LS contribute to the disruption of these two cellular pathways. In addition, we have studied changes in RNAPII LS in two tumors for which clinical responses were reported in phase I/II clinical trials: renal cell carcinoma and Ewing's sarcoma. Our results demonstrate that Et-743 induces degradation of the RNAPII LS that is dependent on active transcription, a functional 26S proteasome, and requires functional TC-NER, but not global genome repair. Additionally, we have provided the first experimental data indicating that degradation of RNAPII LS might lead to the inhibition of activated gene transcription. A set of studies performed in isogenic renal carcinoma cells deficient in von Hippel-Lindau protein, which is a ubiquitin-E3-ligase for RNAPII LS, confirmed the central role of RNAPII LS degradation in the sensitivity to Et-743. Finally, we have shown that RNAPII LS is also degraded in Ewing's sarcoma tumors following Et-743 treatment and provide data to suggest that this event plays a role in decreased expression of the Ewing's sarcoma oncoprotein, EWS-Fli1. Altogether, these data implicate degradation of RNAPII LS as a critical event following Et-743 exposure and suggest that the clinical activity observed in renal carcinoma and Ewing's sarcoma may be mediated by disruption of molecular pathways requiring a fully functional RNAPII LS. ^

Regulation of RNA splicing by nonsense mutations

Translation termination as a result of premature nonsense codon-incorporation in a RNA transcript can lead to the production of aberrant proteins with gain-of-function or dominant negative properties that could have deletrious effects on the cell. T-cell Receptor (TCR) genes acquire premature termination codons two-thirds of the time as a result of the error-prone programmed rearrangement events that normally occur during T-cell development. My studies have focused on the fate of TCR precursor mRNAs in response to in-frame nonsense mutations. ^ Previous published studies from our laboratory have shown that TCR precursor mRNAs are subject to nonsense mediated upregulation of pre-mRNA (NMUP). In this dissertation, I performed substitution and deletion analysis to characterize specific regions of TCR which are required to elicit NMUP. I performed frame- and factor-dependence studies to determine its relationship with other nonsense codon induced responses using several approaches including (i) translation dependence studies (ii) deletion and mutational analysis, as well as (iii) siRNA mediated knockdown of proteins involved. I also addressed the underlying molecular mechanism for this pre-mRNA upregulation by (i) RNA half-life studies using a c-fos inducible promoter, and (ii) a variety of assays to determine pre-mRNA splicing efficiency. ^ Using these approaches, I have identified a region of TCR that is both necessary and sufficient to elicit (NMUP). I have also found that neither cytoplasmic translation machinery nor the protein UPF1 are involved in eliciting this nuclear event. I have shown that the NMUP can be induced not only by nonsense and frameshift mutations, but also missense mutations that disrupt a cis splicing element in the exon that contains the mutation. However, the effect of nonsense mutations on pre-mRNA is unique and distinguishable from that of missense mutations in that nonsense mutations can upregulate pre-mRNA in a frame-dependent manner. Lastly, I provide evidence that NMUP occurs by a mechanism in which nonsense mutations inhibit the splicing of introns. In summary, I have found that TCR precursor mRNAs are subject to multiple forces involving both RNA splicing and translation that can either increase or decrease the levels of these precursor mRNAs. ^

Glioma-specific alternative RNA splicing: A role for polypyrimidine tract binding protein-1

Alternative RNA splicing plays an integral role in cell fate determination and function, especially in the cells of the brain. Errors in RNA processing contribute to diseases such as cancer, where it leads to the production of oncogenic proteins or the loss of tumor suppressors. In silica mining suggests that hundreds of splice isoforms are misexpressed in the glial cell-derived glioma. However, there is little experimental evidence of the prevalence and contribution of these changes and whether they contribute to the formation and progression of this devastating malignancy. To determine the frequency of these aberrant events, global profiling of alternative RNA splice patterns in glioma and nontumor brain was conducted using an exon array. Most splicing changes were less than 5-fold in magnitude and 14 cassette exon events were validated, including 7 previously published events. To determine the possible causes of missplicing, the differential expression levels of splicing factors in these two tissues were also analyzed. Six RNA splicing factors had greater than 2-fold changes in expression. The highest differentially expressed factor was polypyrimidine tract binding protein-1 (PTB). Evaluation by immunohistochemistry determined that this factor was elevated in both early and late stages of glioma. Glial cell-specific PTB expression in the adult brain led me to examine the role of PTB in gliomagenesis. Downregulation of PTB slowed glioma cell proliferation and migration and enhanced cell adhesion to fibronectin and vitronectin. To determine whether PTB was affecting these processes through splicing, genome-wide exon expression levels were correlated with PTB levels. Surprisingly, previously reported PTB target transcripts were insensitive to changes in PTB levels in both patient samples and PTB-depleted glioma cells. Only one validated glioma-specific splice target, RTN4/Nogo, had a significant PTB-mediated splicing change. Downregulation of PTB enhanced inclusion of its alternative exon 3, which encodes an auxiliary domain within a neurite inhibitor protein. Overexpression of this splice isoform in glioma cells slowed proliferation in a manner similar to that observed in PTB knockdown cells. In summary, aberrant expression of splicing factors such as PTB in glioma may elicit changes in splicing patterns that enhance tumorigenesis. ^

Micro RNA expression alterations in patients with glioblastoma

Glioblastoma, also known as glioblastoma multiform or GBM, is the most common and most malignant primary brain tumor. The clinical history of patients with glioblastoma is short, usually less than 3 months in more than 50% of cases after diagnosis. Currently, the methods of glioblastoma treatment are chemotherapy, radiotherapy and surgery. Even with the more effective treatment options, patients with glioblastoma most likely have a median survival time of 10 to 12 months. It is necessary to seek other treatment methods, including gene-targeted treatment. The success of gene-targeted treatment depends critically on the knowledge of genes that may be the cause of, or contribute to disease. To establish a correlate between glioblastoma survival timeline and micro RNA expression alteration, a study of 91 glioblastoma patients was conducted at the University of Texas M. D. Anderson Cancer Center. These 91 glioblastoma patients were newly diagnosed from 2002 to 2007. Statistical analysis was conducted to test the association of miRNA expression alteration between long-term survival and short-term survival glioblastoma. The completion of this proposed study will provide a better understanding of the regulatory role of miRNA in glioblastoma progression.^

Post-transcriptional Regulation of Mammalian Gene Expression in Non-coding Region of Target RNA

Tumor Suppressor Candidate 2 (TUSC2) is a novel tumor suppressor gene located in the human chromosome 3p21.3 region. TUSC2 mRNA transcripts could be detected on Northern blots in both normal lung and some lung cancer cell lines, but no endogenous TUSC2 protein could be detected in a majority of lung cancer cell lines. Mechanisms regulating TUSC2 protein expression and its inactivation in primary lung cancer cells are largely unknown. We investigated the role of the 5’- and 3’-untranslated regions (UTRs) of the TUSC2 gene in the regulation of TUSC2 protein expression. We found that two small upstream open-reading frames (uORFs) in the 5’UTR of TUSC2 could markedly inhibit the translational initiation of TUSC2 protein by interfering with the “scanning” of the ribosome initiation complexes. Site-specific stem-loop array reverse transcription-polymerase chain reaction (SLA-RT-PCR) verified several micoRNAs (miRNAs) targeted at 3’UTR and directed TUSC2 cleavage and degradation. In addition, we used the established let-7-targeted high mobility group A2 (Hmga2) mRNA as a model system to study the mechanism of regulation of target mRNA by miRNAs in mammalian cells under physiological conditions. There have been no evidence of direct link between mRNA downregulation and mRNA cleavages mediated by miRNAs. Here we showed that the endonucleolytic cleavages on mRNAs were initiated by mammalian miRNA in seed pairing style. Let-7 directed cleavage activities among the eight predicted potential target sites have varied efficiency, which are influenced by the positional and the structural contexts in the UTR. The 5’ cleaved RNA fragments were mostly oligouridylated at their 3’-termini and accumulated for delayed 5’–3’ degradation. RNA fragment oligouridylation played important roles in marking RNA fragments for delayed bulk degradation and in converting RNA degradation mode from 3’–5’ to 5’–3’ with cooperative efforts from both endonucleolytic and non-catalytic miRNA-induced silencing complex (miRISC). Our findings point to a mammalian miRNA-mediated mechanism for the regulation of mRNA that miRNA can decrease target mRNA through target mRNA cleavage and uridine addition

GENETIC INTERACTIONS OF FACTORS THAT REGULATE ALTERNATIVE RNA SPLICING IN THE MALE GERM LINE OF DROSOPHILA

Alternative RNA splicing is a critical process that contributes variety to protein functions, and further controls cell differentiation and normal development. Although it is known that most eukaryotic genes produce multiple transcripts in which splice site selection is regulated, how RNA binding proteins cooperate to activate and repress specific splice sites is still poorly understood. In addition how the regulation of alternative splicing affects germ cell development is also not well known. In this study, Drosophila Transformer 2 (Tra2) was used as a model to explore both the mechanism of its repressive function on its own pre-mRNA splicing, and the effect of the splicing regulation on spermatogenesis in testis. Half-pint (Hfp), a protein known as splicing activator, was identified in an S2 cell-based RNAi screen as a co-repressor that functions in combination with Tra2 in the splicing repression of the M1 intron. Its repressive splicing function is found to be sequence specific and is dependent on both the weak 3’ splice site and an intronic splicing silencer within the M1 intron. In addition we found that in vivo, two forms of Hfp are expressed in a cell type specific manner. These alternative forms differ at their amino terminus affecting the presence of a region with four RS dipeptides. Using assays in Drosophila S2 cells, we determined that the alternative N terminal domain is necessary in repression. This difference is probably due to differential localization of the two isoforms in the nucleus and cytoplasm. Our in vivo studies show that both Hfp and Tra2 are required for normal spermatogenesis and cooperate in repression of M1 splicing in spermatocytes. But interestingly, Tra2 and Hfp antagonize each other’s function in regulating germline specific alternative splicing of Taf1 (TBP associated factor 1). Genetic and cytological studies showed that mutants of Hfp and Taf1 both cause similar defects in meiosis and spermatogenesis. These results suggest Hfp regulates normal spermatogenesis partially through the regulation of taf1 splicing. These observations indicate that Hfp regulates tra2 and taf1 activity and play an important role in germ cell differentiation of male flies.

Modeling Techniques for the High-Resolution Interpretation of Cryo-Electron Microscopy Reconstructions

Essential biological processes are governed by organized, dynamic interactions between multiple biomolecular systems. Complexes are thus formed to enable the biological function and get dissembled as the process is completed. Examples of such processes include the translation of the messenger RNA into protein by the ribosome, the folding of proteins by chaperonins or the entry of viruses in host cells. Understanding these fundamental processes by characterizing the molecular mechanisms that enable then, would allow the (better) design of therapies and drugs. Such molecular mechanisms may be revealed trough the structural elucidation of the biomolecular assemblies at the core of these processes. Various experimental techniques may be applied to investigate the molecular architecture of biomolecular assemblies. High-resolution techniques, such as X-ray crystallography, may solve the atomic structure of the system, but are typically constrained to biomolecules of reduced flexibility and dimensions. In particular, X-ray crystallography requires the sample to form a three dimensional (3D) crystal lattice which is technically di‑cult, if not impossible, to obtain, especially for large, dynamic systems. Often these techniques solve the structure of the different constituent components within the assembly, but encounter difficulties when investigating the entire system. On the other hand, imaging techniques, such as cryo-electron microscopy (cryo-EM), are able to depict large systems in near-native environment, without requiring the formation of crystals. The structures solved by cryo-EM cover a wide range of resolutions, from very low level of detail where only the overall shape of the system is visible, to high-resolution that approach, but not yet reach, atomic level of detail. In this dissertation, several modeling methods are introduced to either integrate cryo-EM datasets with structural data from X-ray crystallography, or to directly interpret the cryo-EM reconstruction. Such computational techniques were developed with the goal of creating an atomic model for the cryo-EM data. The low-resolution reconstructions lack the level of detail to permit a direct atomic interpretation, i.e. one cannot reliably locate the atoms or amino-acid residues within the structure obtained by cryo-EM. Thereby one needs to consider additional information, for example, structural data from other sources such as X-ray crystallography, in order to enable such a high-resolution interpretation. Modeling techniques are thus developed to integrate the structural data from the different biophysical sources, examples including the work described in the manuscript I and II of this dissertation. At intermediate and high-resolution, cryo-EM reconstructions depict consistent 3D folds such as tubular features which in general correspond to alpha-helices. Such features can be annotated and later on used to build the atomic model of the system, see manuscript III as alternative. Three manuscripts are presented as part of the PhD dissertation, each introducing a computational technique that facilitates the interpretation of cryo-EM reconstructions. The first manuscript is an application paper that describes a heuristics to generate the atomic model for the protein envelope of the Rift Valley fever virus. The second manuscript introduces the evolutionary tabu search strategies to enable the integration of multiple component atomic structures with the cryo-EM map of their assembly. Finally, the third manuscript develops further the latter technique and apply it to annotate consistent 3D patterns in intermediate-resolution cryo-EM reconstructions. The first manuscript, titled An assembly model for Rift Valley fever virus, was submitted for publication in the Journal of Molecular Biology. The cryo-EM structure of the Rift Valley fever virus was previously solved at 27Å-resolution by Dr. Freiberg and collaborators. Such reconstruction shows the overall shape of the virus envelope, yet the reduced level of detail prevents the direct atomic interpretation. High-resolution structures are not yet available for the entire virus nor for the two different component glycoproteins that form its envelope. However, homology models may be generated for these glycoproteins based on similar structures that are available at atomic resolutions. The manuscript presents the steps required to identify an atomic model of the entire virus envelope, based on the low-resolution cryo-EM map of the envelope and the homology models of the two glycoproteins. Starting with the results of the exhaustive search to place the two glycoproteins, the model is built iterative by running multiple multi-body refinements to hierarchically generate models for the different regions of the envelope. The generated atomic model is supported by prior knowledge regarding virus biology and contains valuable information about the molecular architecture of the system. It provides the basis for further investigations seeking to reveal different processes in which the virus is involved such as assembly or fusion. The second manuscript was recently published in the of Journal of Structural Biology (doi:10.1016/j.jsb.2009.12.028) under the title Evolutionary tabu search strategies for the simultaneous registration of multiple atomic structures in cryo-EM reconstructions. This manuscript introduces the evolutionary tabu search strategies applied to enable a multi-body registration. This technique is a hybrid approach that combines a genetic algorithm with a tabu search strategy to promote the proper exploration of the high-dimensional search space. Similar to the Rift Valley fever virus, it is common that the structure of a large multi-component assembly is available at low-resolution from cryo-EM, while high-resolution structures are solved for the different components but lack for the entire system. Evolutionary tabu search strategies enable the building of an atomic model for the entire system by considering simultaneously the different components. Such registration indirectly introduces spatial constrains as all components need to be placed within the assembly, enabling the proper docked in the low-resolution map of the entire assembly. Along with the method description, the manuscript covers the validation, presenting the benefit of the technique in both synthetic and experimental test cases. Such approach successfully docked multiple components up to resolutions of 40Å. The third manuscript is entitled Evolutionary Bidirectional Expansion for the Annotation of Alpha Helices in Electron Cryo-Microscopy Reconstructions and was submitted for publication in the Journal of Structural Biology. The modeling approach described in this manuscript applies the evolutionary tabu search strategies in combination with the bidirectional expansion to annotate secondary structure elements in intermediate resolution cryo-EM reconstructions. In particular, secondary structure elements such as alpha helices show consistent patterns in cryo-EM data, and are visible as rod-like patterns of high density. The evolutionary tabu search strategy is applied to identify the placement of the different alpha helices, while the bidirectional expansion characterizes their length and curvature. The manuscript presents the validation of the approach at resolutions ranging between 6 and 14Å, a level of detail where alpha helices are visible. Up to resolution of 12 Å, the method measures sensitivities between 70-100% as estimated in experimental test cases, i.e. 70-100% of the alpha-helices were correctly predicted in an automatic manner in the experimental data. The three manuscripts presented in this PhD dissertation cover different computation methods for the integration and interpretation of cryo-EM reconstructions. The methods were developed in the molecular modeling software Sculptor (http://sculptor.biomachina.org) and are available for the scientific community interested in the multi-resolution modeling of cryo-EM data. The work spans a wide range of resolution covering multi-body refinement and registration at low-resolution along with annotation of consistent patterns at high-resolution. Such methods are essential for the modeling of cryo-EM data, and may be applied in other fields where similar spatial problems are encountered, such as medical imaging.

RNA-SEQUENCING APPLICATIONS: GENE EXPRESSION QUANTIFICATION AND METHYLATOR PHENOTYPE IDENTIFICATION

My dissertation focuses on two aspects of RNA sequencing technology. The first is the methodology for modeling the overdispersion inherent in RNA-seq data for differential expression analysis. This aspect is addressed in three sections. The second aspect is the application of RNA-seq data to identify the CpG island methylator phenotype (CIMP) by integrating datasets of mRNA expression level and DNA methylation status. Section 1: The cost of DNA sequencing has reduced dramatically in the past decade. Consequently, genomic research increasingly depends on sequencing technology. However it remains elusive how the sequencing capacity influences the accuracy of mRNA expression measurement. We observe that accuracy improves along with the increasing sequencing depth. To model the overdispersion, we use the beta-binomial distribution with a new parameter indicating the dependency between overdispersion and sequencing depth. Our modified beta-binomial model performs better than the binomial or the pure beta-binomial model with a lower false discovery rate. Section 2: Although a number of methods have been proposed in order to accurately analyze differential RNA expression on the gene level, modeling on the base pair level is required. Here, we find that the overdispersion rate decreases as the sequencing depth increases on the base pair level. Also, we propose four models and compare them with each other. As expected, our beta binomial model with a dynamic overdispersion rate is shown to be superior. Section 3: We investigate biases in RNA-seq by exploring the measurement of the external control, spike-in RNA. This study is based on two datasets with spike-in controls obtained from a recent study. We observe an undiscovered bias in the measurement of the spike-in transcripts that arises from the influence of the sample transcripts in RNA-seq. Also, we find that this influence is related to the local sequence of the random hexamer that is used in priming. We suggest a model of the inequality between samples and to correct this type of bias. Section 4: The expression of a gene can be turned off when its promoter is highly methylated. Several studies have reported that a clear threshold effect exists in gene silencing that is mediated by DNA methylation. It is reasonable to assume the thresholds are specific for each gene. It is also intriguing to investigate genes that are largely controlled by DNA methylation. These genes are called “L-shaped” genes. We develop a method to determine the DNA methylation threshold and identify a new CIMP of BRCA. In conclusion, we provide a detailed understanding of the relationship between the overdispersion rate and sequencing depth. And we reveal a new bias in RNA-seq and provide a detailed understanding of the relationship between this new bias and the local sequence. Also we develop a powerful method to dichotomize methylation status and consequently we identify a new CIMP of breast cancer with a distinct classification of molecular characteristics and clinical features.

Estrogen induced desensitization of c -{\it fos\/} messenger RNA expression {\it in vivo\/}

In this thesis, we investigated the regulation of the nuclear proto-oncogene, c-fos by estrogen in vivo. In the uterus, estrogen causes a rapid, dramatic and transient induction of c-fos mRNA and this occurs by transcriptional activation. We have discovered a previously unrecognized regulatory mechanism by which fos becomes desensitized to estrogen following the transient induction. We investigated three aspects of this desensitization: (1) the kinetics and general characteristics of the phenomenon; (2) the molecular mechanism of the desensitization; and (3) the relationship of desensitization to estrogen stimulated DNA synthesis. The desensitization occurs between 3-24 hours after initial hormonal stimulation and is reversible within 72 hours. The desensitization is not species specific, in that it occurs in both the rat and mouse. The desensitization also occurs in at least two estrogen responsive tissues, the uterus and vagina. The desensitization is not unique to c-fos, since both c-myc and c-jun show similar patterns of desensitization. However, the desensitization is not observed with creatine kinase B (CKB), indicating that not all estrogen inducible genes become desensitized. In the second general area, we determined the desensitization is at the transcriptional level. The desensitization is homologous, but not heterologous, since estrogen induction does not desensitize c-fos to other agents. Other studies show that the desensitization is not due to the lack of functional estrogen receptors. Taken together, these findings suggest that the desensitization occurs at the level of the estrogen responsive element. In the third major area, we demonstrated that the desensitization appears to be related to estrogen induced DNA synthesis. Support for this suggestion comes from the observation that short acting estrogens which induce fos, but not DNA synthesis, do not produce desensitization. ^

Inhibitor of cytochrome c messenger RNA -protein interaction in stimulated rat skeletal muscle

The purpose of this work was to examine the possible mechanisms for the regulation of cytochrome c gene expression in response to increased contractile activity in rat skeletal muscle. The working hypothesis was that increased contractile activity enhances cytochrome c gene expression through a cis-element. A 110% increase in cytochrome c mRNA concentration was observed in tibialis anterior (TA) muscle after 9 days of chronic stimulation. Similar difference (120%) exists between soleus (SO) muscle of higher contractile activity and white vastus lateralis (WV) muscle of lower contractile activity. These results suggest that the endogenous cytochrome c gene expression is regulated by contractile activity. Cytochrome c-reporter genes were injected into skeletal muscles to identify the cis-element that is responsible for the regulation. Although the data was inconclusive, part of it suggested the importance of the 3$\sp\prime$-untranslated region (3$\sp\prime$-UTR) in mediating the response to increased contractile activity.^ RNA gel mobility shift (GMSA) and ultraviolet (UV) cross-linking assays revealed specific RNA-protein interaction in a 50-nucleotide region of the 3$\sp\prime$-UTR in unstimulated TA muscle. Computer analysis predicted a stem-loop structure of 17 nucleotides, which provides a structural basis for RNA-protein interaction. These 17 nucleotides are 100% conserved among rat, mouse and human cytochrome c genes and their 13 pseudogenes, suggesting a functional role for this region. The RNA-protein interaction was significantly less in highly active SO muscle than in inactive WV muscle and was dramatically decreased in stimulated TA muscle due to a protein inhibitor(s) associated with ribosome. It is possible that cytochrome c mRNAs undergoing translation are subject to a compartmentalized regulatory influence.^ The conclusion from these results is that increases in contractile activity induce or activate a protein inhibitor(s) associated with ribosome in rat skeletal muscle. The inhibitor decreases RNA-protein interaction in the 3$\sp\prime$-UTR of cytochrome c mRNA, which may result in increased mRNA stability and/or translation. ^