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.

MOLECULAR AND GENOMIC BASED INSIGHTS INTO THE EVOLUTION OF ENTEROCOCCUS FAECIUM FROM COMMENSAL TO HOSPITAL-ADAPTED PATHOGEN

The basis for the recent transition of Enterococcus faecium from a primarily commensal organism to one of the leading causes of hospital-acquired infections in the United States is not yet understood. To address this, the first part of my project assessed isolates from early outbreaks in the USA and South America using sequence analysis, colony hybridizations, and minimal inhibitory concentrations (MICs) which showed clinical isolates possess virulence and antibiotic resistance determinants that are less abundant or lacking in community isolates. I also revealed that the level of ampicillin resistance increased over time in clinical strains. By sequencing the pbp5 gene, I demonstrated an ~5% difference in the pbp5 gene between strains with MICs <4ug/ml and those with MICs >4µg/ml, but no specific sequence changes correlated with increases in MICs within the latter group. A 3-10% nucleotide difference was also seen in three other genes analyzed, which suggested the existence of two distinct subpopulations of E. faecium. This led to the second part of my project analyzing concatenated core gene sequences, SNPs, the 16S rRNA, and phylogenetics of 21 E. faecium genomes confirming two distinct clades; a community-associated (CA) clade and hospital-associated (HA) clade. Molecular clock calculations indicate that these two clades likely diverged ~ 300,000 to > 1 million years ago, long before the modern antibiotic era. Genomic analysis also showed that, in addition to core genomic differences, HA E. faecium harbor specific accessory genetic elements that may confer selection advantages over CA E. faecium. The third part of my project discovered 6 E. faecium genes with the newly identified “WxL” domain. My analyses, using RT-PCR, western blots, patient sera, whole-cell ELISA, and immunogold electron microscopy, indicated that E. faecium WxL genes exist in operons, encode bacterial cell surface localized proteins, that WxL proteins are antigenic in humans, and are more exposed on the surface of clinical isolates versus community isolates (even though they are ubiquitous in both clades). ELISAs and BIAcore analyses also showed that proteins encoded by these operons bind several different host extracellular matrix proteins, as well as to each other, suggesting a novel cell-surface complex. In summary, my studies provide new insights into the evolution of E. faecium by showing that there are two distantly related clades; one being more successful in the hospital setting. My studies also identified operons encoding WxL proteins whose characteristics could also contribute to colonization and virulence within this species.

Mechanisms Underlying the Heterogeneous Sensitivities of Cancer Cells to Proteasome Inhibitors

The mechanisms underlying cellular response to proteasome inhibitors have not been clearly elucidated in solid tumor models. Evidence suggests that the ability of a cell to manage the amount of proteotoxic stress following proteasome inhibition dictates survival. In this study using the FDA-approved proteasome inhibitor bortezomib (Velcade®) in solid tumor cells, we demonstrated that perhaps the most critical response to proteasome inhibition is repression of global protein synthesis by phosphorylation of the eukaryotic initiation factor 2-α subunit (eIF2α). In a panel of 10 distinct human pancreatic cancer cells, we showed marked heterogeneity in the ability of cancer cells to induce eIF2α phosphorylation upon stress (eIF2α-P); lack of inducible eIF2α-P led to excessive accumulation of aggregated proteins, reactive oxygen species, and ultimately cell death. In addition, we examined complementary cytoprotective mechanisms involving the activation of the heat shock response (HSR), and found that induction of heat shock protein 70 kDa (Hsp72) protected against proteasome inhibitor-induced cell death in human bladder cancer cells. Finally, investigation of a novel histone deacetylase 6 (HDAC6)-selective inhibitor suggested that the cytoprotective role of the cytoplasmic histone deacetylase 6 (HDAC6) in response to proteasome inhibition may have been previously overestimated.

Interview with Dianna Milewicz

An oral interview with Dianna Milewicz, M.D., Ph.D., Professor and Vice Chairman of the Department of Internal Medicine at The University of Texas Medical School at Houston and Director of the M.D./Ph.D. Program and co-Director of the Biomedical Engineering Center. Her research interests include the genetic basis of cardiovascular diseases, and understanding the effect of identified mutations on protein function. She has recently established a genetic core laboratory to provide molecular biology and genetic expertise to clinicians who want to initiate genetic studies on their patient populations.

Interview with Anna Steinberger

An oral interview with Dr. Anna Steinberger, who taught and conducted basic research in Reproductive Biology and served as Assistant Dean for Faculty Affairs at UT Medical School-Houston. Her research yielded over 250 scientific articles, books, and book chapters for which she received numerous awards and recognitions in the USA and abroad.

GENOME-WIDE PROFILING UNVEILS CRITICIAL FUNCTIONS OF p53 IN HUMAN EMBRYONIC STEM CELLS

Embryonic stem cells (ESCs) possess two unique characteristics: infinite self-renewal and the potential to differentiate into almost every cell type (pluripotency). Recently, global expression analyses of metastatic breast and lung cancers revealed an ESC-like expression program or signature, specifically for cancers that are mutant for p53 function. Surprisingly, although p53 is widely recognized as the guardian of the genome, due to its roles in cell cycle checkpoints, programmed cell death or senescence, relatively little is known about p53 functions in normal cells, especially in ESCs. My hypothesis is that p53 has specific transcription regulatory functions in human ESCs (hESCs) that a) oppose pluripotency and b) protect the stem cell genome in response to DNA damage and stress signaling. In mouse ESCs, these roles are believed to coincide, as p53 promotes differentiation in response to DNA damage, but this is unexplored in hESCs. To determine the biological roles of p53, specifically in hESCs, we mapped genome-wide chromatin interactions of p53 by chromatin immunoprecipitation and massively parallel tag sequencing (ChIP-Seq), and did so under three VIdifferent conditions of hESC status: pluripotency, differentiation-initiated and DNA-damage-induced. ChIP-Seq showed that p53 is enriched at distinct, induction-specific gene loci during each of these different conditions. Microarray gene expression analysis and functional annotation of the distinct p53-target genes revealed that p53 regulates specific genes encoding developmental regulators, which are expressed in differentiation-initiated but not DNA- damaged hESCs. We further discovered that, in response to differentiation signaling, p53 binds regions of chromatin that are repressed but also poised for rapid activation by core pluripotency factors OCT4 and NANOG in pluripotent hESCs. In response to DNA damage, genes associated with migration and motility are targeted by p53; whereas, the prime targets of p53 in control of cell death are conserved for p53 regulation in both differentiation and DNA damage. Our genome-wide profiling and bioinformatics analyses show that p53 occupies a special set of developmental regulatory genes during early differentiation of hESCs and functions in an induction-specific manner. In conclusion, our research unveiled previously unknown functions of p53 in ESC biology, which augments our understanding of one of the most deregulated proteins in human cancers.

INTERACTION BETWEEN BRK AND HER2 IN BREAST CANCER

INTERACTION BETWEEN BRK AND HER2 IN BREAST CANCER Midan Ai, Ph.D. Supervisory Professor: Zhen Fan, M.D. Breast tumor kinase (Brk) is a nonreceptor protein-tyrosine kinase that is highly expressed in approximately two thirds of breast cancers but is not detectable or is expressed at very low levels in normal mammary epithelium. Brk plays important roles in promoting proliferation, survival, invasion, and metastasis of breast cancer cells, but the mechanism(s) of which remain largely unknown. Recent studies showed that Brk is frequently co-overexpressed with human epidermal growth factor receptor-2 (HER2) and is physically associated with HER2 in breast cancer. The mechanism needs to be determined. In my studies, I found that high expression of HER2 is correlated with high expression of Brk in breast cancer cell lines. Silencing HER2 expression via RNA interference in HER2 over-expressed breast cancer cells resulted in Brk protein decrease and overexpression of HER2 in HER2 low-expressed breast cancer cells up-regulated Brk expression. The mechanism study indicated that overexpression of HER2 increased Brk protein stability. Brk was degraded through a Ca2+-dependent protease pathway involving calpain and HER2 stimulated Brk expression via inhibiting calpain activity. Calpastatin is a calpain endogenous inhibitor and the calpain-calpastatin system has been implicated in a number of cell physiological functions. HER2 restrained calpain activation via up-regulating calpastatin expression and HER2 downstream signaling, MAPK pathway, was involved in the regulation. Furthermore, silencing of Brk expression by RNA interference in HER2-overexpressing breast cancer cells decreased HER2-mediated cell proliferation, survival, invasion/metastasis potential and increased cell sensitivity to HER2 kinase inhibitor, lapatinib, treatment, indicating that Brk plays important roles in regulating and mediating the oncogenic functions of HER2. The Stat3 pathway played important roles in Brk mediated cell survival and invasion/metastasis in the context of HER2-overexpressing breast cancer cells. However, transgenic mice with inducible expression of constitutively active Brk (CA) in the mammary epithelium failed to develop malignant change in the mammary glands after Brk induction for 15 months which indicated that expression of Brk protein alone was not sufficiently to induce spontaneous breast tumor. Bitransgenic mice with co-expression of HER2/neu and inducible expression of Brk in the mammary epithelium developed multifocal mammary tumors, but there were no significant difference in the tumor occurring time, tumor size, tumor weight and tumor multiplicity between the mouse group with co-expression of Brk and HER2/neu and the mouse group with HER2/neu expression only.

Acceleration of the PanIN Development in Mice Expressing Oncogenic K-Ras Due to a High Fat Diet

Obesity is postulated to be one of the major risk factors for pancreatic cancer, and recently it was indicated that an elevated body mass index (BMI correlates strongly with a decrease in patient survival. Despite the evident relationship, the molecular mechanisms involved are unclear. Oncogenic mutation of K-Ras is found early and is universal in pancreatic cancer. Extensive evidence indicates oncogenic K-Ras is not entirely active and it requires a triggering event to surpass the activity of Ras beyond the threshold necessary for a Ras-inflammation feed-forward loop. We hypothesize that high fat intake induces a persistent low level inflammatory response triggering increased K-Ras activity and that Cox-2 is essential for this inflammatory reaction. To determine this, LSL-K-Ras mice were crossed with Ela-CreER (Acinar-specific) or Pdx-1-Cre (Pancreas-specific) to “knock-in” oncogenic K-Ras. Additionally, these animals were crossed with Cox-2 conditional knockout mice to access the importance of Cox-2 in the inflammatory loop present. The mice were fed isocaloric diets containing 60% energy or 10% energy from fat. We found that a high fat diet increased K-Ras activity, PanIN formation, and fibrotic stroma significantly compared to a control diet. Genetic deletion of Cox-2 prevented high fat diet induced fibrosis and PanIN formation in oncogenic K-Ras expressing mice. Additionally, long term consumption of high fat diet, increased the progression of PanIN lesions leading to invasive cancer and decreased overall survival rate. These findings indicate that a high fat diet can stimulate the activation of oncogenic K-Ras and initiate an inflammatory feed forward loop requiring Cox-2 leading to inflammation, fibrosis, and PanINs. This mechanism could explain the relationship between a high fat diet and elevated risk for pancreatic cancer.

ROLE OF MIR-19A RELEASED BY INFLAMMATORY BREAST CANCER CELLS IN THE REGULATION OF DENDRITIC CELL FUNCTIONS: IN VITRO MODEL OF CROSSTALK IN THE TUMOR MICROENVIRONMENT OF INFLAMMATORY BREAST CANCER

Inflammatory breast cancer (IBC) is a rare but very aggressive form of locally advanced breast cancer (1-6% of total breast cancer patients in United States), with a 5-year overall survival rate of only 40.5%, compared with 85% of the non-IBC patients. So far, a unique molecular signature for IBC able to explain the dramatic differences in the tumor biology between IBC and non-IBC has not been identified. As immune cells in the tumor microenvironment plays an important role in regulating tumor progression, we hypothesized that tumor-associated dendritic cells (TADC) may be responsible for regulating the development of the aggressive characteristics of IBC. MiRNAs can be released into the extracellular space and mediate the intercellular communication by regulating target gene expression beyond their cells of origin. We hypothesized that miRNAs released by IBC cells can induce an increased activation status, secretion of pro-inflammatory cytokines and migration ability of TADC. In an in vitro model of IBC tumor microenvironment, we found that the co-cultured of the IBC cell line SUM-149 with immature dendritic cells (iDCSUM-149) induced a higher degree of activation and maturation of iDCSUM-149 upon stimulation with lipopolysaccharide (LPS) compared with iDCs co-cultured with the non-IBC cell line SUM-159 (iDCSUM-159), resulting in: increased expression of the costimulatory and activation markers; higher production of pro-inflammatory cytokines (TNF-a, IL-6); and 3) higher migratory ability. These differences were due to the exosome-mediated transfer of miR-19a and miR-146a from SUM-149 and SUM-159, respectively, to iDCs, causing the downregulation of the miR-19a target genes PTEN, SOCS-1 and the miR-146a target genes IRAK1, TRAF6. PTEN, SOCS-1 and IRAK1, TRAF6 are important negative and positive regulator of cytokine- and TLR-mediated activation/maturation signaling pathway in DCs. Increased levels of IL-6 induced the upregulation of miR-19a synthesis in SUM-149 cells that was associated with the induction of CD44+CD24-ALDH1+ cancer stem cells (CSCs) with epithelial-to-mesenchymal transition (EMT) characteristics. In conclusion, in IBC tumor microenvironment IL-6/miR-19a axis can represent a self-sustaining loop able to maintain a pro-inflammatory status of DCs, leading to the development of tumor cells with high metastatic potential (EMT CSCs) responsible of the poor prognosis in IBC patients.

PRKCa: Identification of a Novel Downstream Target of WT1

Wilms tumor is a childhood tumor of the kidney arising from the undifferentiated metanephric mesenchyme. Tumorigenesis is attributed to a number of genetic and epigenetic alterations. In 20% of Wilms tumors, Wilms tumor gene 1 (WT1) undergoes inactivating homozygous mutations causing loss of function of the zinc finger transcription factor it encodes. It is hypothesized that mutations in WT1 result in dysregulation of downstream target genes, leading to aberrant kidney development and/or Wilms tumor. These downstream target genes are largely unknown, and identification is important for further understanding Wilms tumor development. Heatmap data of human Wilms tumor protein expression, generated by reverse phase protein assay analysis (RPPA), show significant correlation between WT1 mutation status and low PRKCα expression (p= 0.00013); additionally, p-PRKCα (S657) also shows decreased expression in these samples (p= 0.00373). These data suggest that the WT1 transcription factor regulates PRKCα expression, and that PRKCα plays a potential role in Wilms tumor tumorigenesis. We hypothesize that the WT1 transcription factor directly/indirectly regulates PRKCα and mutations occurring in WT1 lead to decreased expression of PRKCα. Prkcα and Wt1 have been shown to co-localize in E14.5 mesenchymal cells of the developing kidney. siRNA knockdown, in-vivo ablation, and tet-inducible expression of Wt1 each independently confirm regulation of Prkcα expression by Wt1 at both RNA and protein levels, and investigation into possible WT1 binding sites in PRKCα regulatory regions has identified multiple sites to be confirmed by luciferase reporter constructs. With the goal of identifying WT1 and PRKCα downstream targets, RPPA analysis of protein expression in mesenchymal cell culture, following lentiviral delivered shRNA knockdown of Wt1 and shRNA knockdown of Prkcα, will be carried out. Apart from Wilms tumor, WT1 also plays an important role in Acute Myeloid Leukemia (AML). WT1 mutation status has been implicated, controversially, as an independent poor-prognosis factor in leukemia, leading to decreased probability of overall survival, complete remission, and disease free survival. RPPA analysis of AML patient samples showed significant decreases in PRKCα/p-PRKCα protein expression in a subset of patients (Kornblau, personal communication); therefore, the possible role of WT1 and PRKCα in leukemia disease progression is an additional focus of this study. WT1 mutation analysis of diploid leukemia patient samples revealed two patients with mutations predicted to affect WT1 activity; of these two samples, only one corresponded to the low PRKCα expression cohort. Further characterization of the role of WT1 in AML, and further understanding of WT1 regulated PRKCα expression, will be gained following RPPA analysis of protein expression in HL60 leukemia cell lines with lentiviral delivered shRNA knockdown of WT1 and shRNA knockdown of PRKCα.

CONFORMATIONAL DYNAMICS OF K-RAS AND H-RAS PROTEINS: IS THERE FUNCTIONAL SPECIFICITY AT THE CATALYTIC DOMAIN?

Ras proteins serve as crucial signaling modulators in cell proliferation through their ability to hydrolyze GTP and exist in a GTP “on” state and GTP “off” state. There are three different human Ras isoforms: H-ras, N-ras and K-ras (4A and 4B). Although their sequence identity is very high at the catalytic domain, these isoforms differ in their ability to activate different effectors and hence different signaling pathways. Much of the previous work on this topic has attributed this difference to the hyper variable region of Ras proteins, which contains most of the sequence variance among the isoforms and encodes specificity for differential distribution in the membrane. However, we hypothesize that sequence variation on lobe II of Ras catalytic domain alters dynamics and leads to differential preference for different effectors or modulators. In this work, we used all atom molecular dynamics to analyze the dynamics in the catalytic domain of H-ras and K-ras. We have also analyzed the dynamics of a transforming mutant of H-ras and K-ras and further studied the dynamics of an effectorselective mutant of H-ras. Collectively we have determined that wild type K-ras is more dynamic than H-ras and that the structure of the effector binding loop more closely resembles that of the T35S Raf-selective mutant, possibly giving us a new view and insight into the v mode of effector specificity. Furthermore we have determined that specific mutations at the same location perturb the conformational equilibrium differently in H-ras and K-ras and that an enhanced oncogenic potential may arise from different structural perturbations for each point mutation of a specific isoform.

STUDYING AGGREGATE FORMATION BY AMYOTROPHIC LATERAL SCLEROSIS-ASSOCIATED MUTANT SOD1 PROTEIN IN DROSOPHILA MODEL

A common pathological hallmark of most neurodegenerative disorders is the presence of protein aggregates in the brain. Understanding the regulation of aggregate formation is thus important for elucidating disease pathogenic mechanisms and finding effective preventive avenues and cures. Amyotrophic Lateral Sclerosis (ALS), also known as Lou Gehrig’s disease, is a selective neurodegenerative disorder predominantly affecting motor neurons. The majority of ALS cases are sporadic, however, mutations in superoxide dismutase 1 (SOD1) are responsible for about 20% of familial ALS (fALS). Mutated SOD1 proteins are prone to misfold and form protein aggregates, thus representing a good candidate for studying aggregate formation. The long-term goal of this project is to identify regulators of aggregate formation by mutant SOD1 and other ALS-associated disease proteins. The specific aim of this thesis project is to assess the possibility of using the well-established Drosophila model system to study aggregation by human SOD1 (hSOD1) mutants. To this end, using wild type and the three mutant hSOD1 (A4V, G85R and G93A) most commonly found among fALS, I have generated 16 different SOD1 constructs containing either eGFP or mCherry in-frame fluorescent reporters, established and tested both cell- and animal-based Drosophila hSOD1 models. The experimental strategy allows for clear visualization of ectopic hSOD1 expression as well as versatile co-expression schemes to fully investigate protein aggregation specifically by mutant hSOD1. I have performed pilot cell-transfection experiments and verified induced expression of hSOD1 proteins. Using several tissue- or cell type-specific Gal4 lines, I have confirmed the proper expression of hSOD1 from established transgenic fly lines. Interestingly, in both Drosophila S2 cells and different fly tissues including the eye and motor neurons, robust aggregate formation by either wild type or mutant hSOD1 proteins was not observed. These preliminary observations suggest that Drosophila might not be a good experimental organism to study aggregation and toxicity of mutant hSOD1 protein. Nevertheless this preliminary conclusion implies the potential existence of a potent protective mechanism against mutant hSOD1 aggregation and toxicity in Drosophila. Thus, results from my SOD1-ALS project in Drosophila will help future studies on how to best employ this classic model organism to study ALS and other human brain degenerative diseases.

Development of HIF-1α/HIF-1β heterodimerization inhibitors using a novel bioluminescence reporter assay system for in vitro high throughput screening and in vivo imaging

Tumor growth often outpaces its vascularization, leading to development of a hypoxic tumor microenvironment. In response, an intracellular hypoxia survival pathway is initiated by heterodimerization of hypoxia-inducible factor (HIF)-1α and HIF-1β, which subsequently upregulates the expression of several hypoxia-inducible genes, promotes cell survival and stimulates angiogenesis in the oxygen-deprived environment. Hypoxic tumor regions are often associated with resistance to various classes of radio- or chemotherapeutic agents. Therefore, development of HIF-1α/β heterodimerization inhibitors may provide a novel approach to anti-cancer therapy. To this end, a novel approach for imaging HIF-1α/β heterodimerization in vitro and in vivo was developed in this study. Using this screening platform, we identified a promising lead candidate and further chemically derivatized the lead candidate to assess the structure-activity relationship (SAR). The most effective first generation drug inhibitors were selected and their pharmacodynamics and anti-tumor efficacy in vivo were verified by bioluminescence imaging (BLI) of HIF-1α/β heterodimerization in the xenograft tumor model. Furthermore, the first generation drug inhibitors, M-TMCP and D-TMCP, demonstrated efficacy as monotherapies, resulting in tumor growth inhibition via disruption of HIF-1 signaling-mediated tumor stromal neoangiogenesis.

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.

THE MOLECULAR INTERACTION BETWEEN TYPE II DIABETES AND ALZHEIMER’S DISEASE THROUGH CROSS-SEEDING OF PROTEIN MISFOLDING

With the population of the world aging, the prominence of diseases such as Type II Diabetes (T2D) and Alzheimer’s disease (AD) are on the rise. In addition, patients with T2D have an increased risk of developing AD compared to age-matched individuals, and the number of AD patients with T2D is higher than among aged-matched non-AD patients. AD is a chronic and progressive dementia characterized by amyloid-beta (Aβ) plaques, neurofibrillary tangles (NFTs), neuronal loss, brain inflammation, and cognitive impairment. T2D involves the dysfunctional use of pancreatic insulin by the body resulting in insulin resistance, hyperglycemia, hyperinsulinemia, pancreatic beta cell (β-cell) death, and other complications. T2D and AD are considered protein misfolding disorders (PMDs). PMDs are characterized by the presence of misfolded protein aggregates, such as in T2D pancreas (islet amyloid polypeptide - IAPP) and in AD brain (amyloid– Aβ) of affected individuals. The misfolding and accumulation of these proteins follows a seeding-nucleation model where misfolded soluble oligomers act as nuclei to propagate misfolding by recruiting other native proteins. Cross-seeding occurs when oligomers composed by one protein seed the aggregation of a different protein. Our hypothesis is that the pathological interactions between T2D and AD may in part occur through cross-seeding of protein misfolding. To test this hypothesis, we examined how each respective aggregate (Aβ or IAPP) affects the disparate disease pathology through in vitro and in vivo studies. Assaying Aβ aggregates influence on T2D pathology, IAPP+/+/APPSwe+/- double transgenic (DTg) mice exhibited exacerbated T2D-like pathology as seen in elevated hyperglycemia compared to controls; in addition, IAPP levels in the pancreas are highest compared to controls. Moreover, IAPP+/+/APPSwe+/- animals demonstrate abundant plaque formation and greater plaque density in cortical and hippocampal areas in comparison to controls. Indeed, IAPP+/+/APPSwe+/- exhibit a colocalization of both misfolded proteins in cerebral plaques suggesting IAPP may directly interact with Aβ and aggravate AD pathology. In conclusion, these studies suggest that cross-seeding between IAPP and Aβ may occur, and that these protein aggregates exacerbate and accelerate disease pathology, respectively. Further mechanistic studies are necessary to determine how these two proteins interact and aggravate both pancreatic and brain pathologies.