Regulation of the heat shock response by thiol-reactive compounds in the yeast Saccharomyces cerevisiae

Cells govern their activities and modulate their interactions with the environment to achieve homeostasis. The heat shock response (HSR) is one of the most well studied fundamental cellular responses to environmental and physiological challenges, resulting in rapid synthesis of heat shock proteins (HSPs), which serve to protect cellular constituents from the deleterious effects of stress. In addition to its role in cytoprotection, the HSR also influences lifespan and is associated with a variety of human diseases including cancer, aging and neurodegenerative disorders. In most eukaryotes, the HSR is primarily mediated by the highly conserved transcription factor HSF1, which recognizes target hsp genes by binding to heat shock elements (HSEs) in their promoters. In recent years, significant efforts have been made to identify small molecules as potential pharmacological activators of HSF1 that could be used for therapeutic benefit in the treatment of human diseases relevant to protein conformation. However, the detailed mechanisms through which these molecules drive HSR activation remain unclear. In this work, I utilized the baker's yeast Saccharomyces cerevisiae as a model system to identify a group of thiol-reactive molecules including oxidants, transition metals and metalloids, and electrophiles, as potent activators of yeast Hsf1. Using an artificial HSE-lacZ reporter and the glucocorticoid receptor system (GR), these diverse thiol-reactive compounds are shown to activate Hsf1 and inhibit Hsp90 chaperone complex activity in a reciprocal, dose-dependent manner. To further understand whether cells sense these reactive compounds through accumulation of unfolded proteins, the proline analog azetidine-2-carboxylic acid (AZC) and protein cross-linker dithiobis(succinimidyl propionate) (DSP) were used to force misfolding of nascent polypeptides and existing cytosolic proteins, respectively. Both unfolding reagents display kinetic HSP induction profiles dissimilar to those generated by thiol-reactive compounds. Moreover, AZC treatment leads to significant cytotoxicity, which is not observed in the presence of the thiol-reactive compounds at the concentrations sufficient to induce Hsf1. Additionally, DSP treatment has little to no effect on Hsp90 functions. Together with the ultracentrifugation analysis of cell lysates that detected no insoluble protein aggregates, my data suggest that at concentrations sufficient to induce Hsf1, thiol-reactive compounds do not induce the HSR via a mechanism based on accumulation of unfolded cytosolic proteins. Another possibility is that thiol-reactive compounds may influence aspects of the protein quality control system such as the ubiquitin-proteasome system (UPS). To address this hypothesis, β-galactosidase reporter fusions were used as model substrates to demonstrate that thiol-reactive compounds do not inhibit ubiquitin activating enzymes (E1) or proteasome activity. Therefore, thiol-reactive compounds do not activate the HSR by inhibiting UPS-dependent protein degradation. I therefore hypothesized that these molecules may directly inactivate protein chaperones, known as repressors of Hsf1. To address this possibility, a thiol-reactive biotin probe was used to demonstrate in vitro that the yeast cytosolic Hsp70 Ssa1, which partners with Hsp90 to repress Hsf1, is specifically modified. Strikingly, mutation of conserved cysteine residues in Ssa1 renders cells insensitive to Hsf1 activation by cadmium and celastrol but not by heat shock. Conversely, substitution with the sulfinic acid and steric bulk mimic aspartic acid led to constitutive activation of Hsf1. Cysteine 303, located in the nucleotide-binding/ATPase domain of Ssa1, was shown to be modified in vivo by a model organic electrophile using Click chemistry technology, verifying that Ssa1 is a direct target for thiol-reactive compounds through adduct formation. Consistently, cadmium pretreatment promoted cells thermotolerance, which is abolished in cells carrying SSA1 cysteine mutant alleles. Taken together, these findings demonstrate that Hsp70 acts as a sensor to induce the cytoprotective heat shock response in response to environmental or endogenously produced thiol-reactive molecules and can discriminate between two distinct environmental stressors.

Functional Analysis of Drosophila Integrator Complex in snRNA 3' End Processing

Uridine-rich small nuclear RNAs (U snRNAs) play essential roles in eukaryotic gene expression by facilitating the removal of introns from mRNA precursors and the processing of the replication-dependent histone pre-mRNAs. Formation of the 3’ end of these snRNAs is carried out by a poorly characterized, twelve-membered protein complex named Integrator Complex. In the effort to understand Integrator Complex function in the formation of the snRNA 3’ end, we performed a functional RNAi screen in Drosophila S2 cells to identify protein factors required for snRNA 3’ end formation. This screen was conducted by using a fluorescence-based reporter that elicits GFP expression in response to a deficiency in snRNA processing. Besides scoring the known Integrator subunits, we identified Asunder and CG4785 as additional core members of the Integrator Complex. Additionally, we also found a conserved requirement for Cyclin C and Cdk8 in both fly and human snRNA 3’ end processing. We have further demonstrated that the kinase activity of Cdk8 is critical for snRNA 3’ end processing and is likely to function independent of its well-documented function within the Mediator Cdk8 module. Taken together, this work functionally defines the Drosophila Integrator Complex and demonstrates a novel function for Cyclin C/Cdk8 in snRNA 3’ end formation. This thesis work has also characterized an important functional interaction mediated by a microdomain within Integrator subunit 12 (IntS12) and IntS1 that is required for the activity of the Integrator Complex in processing the snRNA 3’ end. Through the development of a reporter-based functional RNAi-rescue assay in Drosophila S2 cells, we analyzed domains within IntS12 required for snRNA 3’ end formation. This analysis unexpectedly revealed that an N-terminal 30 amino acid region and not the highly conserved central PHD finger domain, is required for snRNA 3’ end cleavage. The IntS12 microdomain (1-45) functions autonomously, and is sufficient to interact and stabilize the putative scaffold protein IntS1. Our findings provide more details of the Integrator Complex for understanding the molecular mechanism of snRNA 3’ end processing. Moreover, these results lay the foundation for future studies of the complex through the identification of a novel functional domain within one subunit and the identification of additional subunits.

Characterization of a Novel Role for the Tousled- like Kinase in Kinetochore Assembly and Function in Caenorhabditis elegans

Chromosome segregation is a critical step during cell division to avoid aneuploidy and promote proper organismal development. Correct sister chromatid positioning and separation during mitosis helps to achieve faithful transmission of genetic material to daughter cells. This prevents improper chromosome partitioning that can potentially result in extrachromosomal fragments, increasing the tumorigenic potential of the cells. The kinetochore is a protenaicious structure responsible for the initiation and orchestration of chromosome movement during mitosis. This highly conserved structure among eukaryotes is required for chromosome attachment to the mitotic spindle and failure to assemble the kinetochore results in aberrant chromosome segregation. Thus elucidating the mechanism of kinetochore assembly is important to have a better understanding of the regulation that controls chromosome segregation. Our previous work identified the C. elegans Tousled-like kinase (TLK-1) as a mitotic kinase and depletion of TLK-1 results in embryonic lethality, characterized by nuclei displaying poor mitotic chromosome alignment, lagging chromosome, and chromosome bridges during anaphase. Additionally, previous studies from our group revealed that TLK-1 is phosphorylated independently by Aurora B at serine 634, and by CHK-1 at threonine T610. The research presented herein reveals that both phosphorylated forms of TLK-1 associate with the kinetochore during mitosis. Moreover, by systematic depletion of kinetochore proteins, I uncovered that pTLK-1 is bona fide kinetochore component that is located at the outer kinetochore layer, influencing the microtubule-binding interface. I also demonstrated that TLK-1 is necessary for the kinetochore localization of the microtubule interacting proteins CLS-2 and LIS-1 and I show that embryos depleted of TLK-1 presented an aberrant twisted kinetochore pattern. Furthermore, I established that the inner kinetochore protein KNL-2 is an in vitro substrate of TLK-1 indicating a possible role of TLK-1 in regulating centromeric assembly. Collectively, these results suggest a novel role for the Tousled-like kinase in regulation of kinetochore assembly and microtubule dynamics and demonstrate the necessity of TLK-1 for proper chromosome segregation in C. elegans.

IDENTIFYING GENETIC VARIANTS AND CHARACTERIZING THEIR ROLE IN CLUBFOOT

Clubfoot is a common, complex birth defect affecting 4,000 newborns in the United States and 135,000 world-wide each year. The clubfoot deformity is characterized by inward and rigid downward displacement of one or both feet, along with persistent calf muscle hypoplasia. Despite strong evidence for a genetic liability, there is a limited understanding of the genetic and environmental factors contributing to the etiology of clubfoot. The studies described in this dissertation were performed to identify variants and/or genes associated with clubfoot. Genome-wide linkage scan performed on ten multiplex clubfoot families identified seven new chromosomal regions that provide new areas to search for clubfoot genes. Troponin C (TNNC2) the strongest candidate gene, located in 20q12-q13.11, is involved in muscle contraction. Exon sequencing of TNNC2 did not identify any novel coding variants. Interrogation of fifteen muscle contraction genes found strong associations with SNPs located in potential regulatory regions of TPM1 (rs4075583 and rs3805965), TPM2 (rs2025126 and rs2145925) and TNNC2 (rs383112 and rs437122). In previous studies, a strong association was found with rs3801776 located in the basal promoter of HOXA9, a gene also involved in muscle development and patterning. Altogether, this data suggests that SNPs located in potential regulatory regions of genes involved in muscle development and function could alter transcription factor binding leading to changes in gene expression. Functional analysis of 3801776/HOXA9, rs2025126/TPM2 and rs2145925/TPM2 showed altered protein binding, which significantly influenced promoter activity. Although the ancestral allele (G) of rs4075583/TPM1 creates a DNA-protein complex, it did not affect TPM1 promoter activity. However and importantly, in the context of a haplotype, rs4075583/G significantly decreased TPM1 promoter activity. These results suggest dysregulation of multiple skeletal muscle genes, TPM1, TPM2, TNNC2 and HOXA9, working in concert may contribute to clubfoot. However, specific allelic combinations involving these four regulatory SNPs did not confer a significantly higher risk for clubfoot. Other combinations of these variants are being evaluated. Moreover, these variants may interact with yet to be discovered variants in other genes to confer a higher clubfoot risk. Collectively, we show novel evidence for the role of skeletal muscle genes in clubfoot indicating that there are multiple genetic factors contributing to this complex birth defect.

Evaluation of Polymer Gel Dosimeters for Measurements of Dose and LET in Proton Beams

This project assessed the effectiveness of polymer gel dosimeters as tools for measuring the dose deposited by and LET of a proton beam. A total of three BANG® dosimeter formulations were evaluated: BANG®-3-Pro-2 BANGkits™ for dose measurement and two BANG®-3 variants, the LET-Baseline and LET-Meter dosimeters, for LET measurement. All dosimeters were read out using an OCT scanner. The basic characteristics of the BANGkits™ were assessed in a series of photon and electron irradiations. The dose-response relationship was found to be sigmoidal with a threshold for response of approximately 15 cGy. The active region of the dosimeter, the volume in which dosimeter response is not inhibited by oxygen, was found to make up roughly one fourth of the total dosimeter volume. Delivering a dose across multiple fractions was found to yield a greater response than delivering the same dose in a single irradiation. The dosimeter was found to accurately measure a dose distribution produced by overlapping photon fields, yielding gamma pass rates of 95.4% and 93.1% from two planar gamma analyses. Proton irradiations were performed for measurements of proton dose and LET. Initial irradiations performed through the side of a dosimeter led to OCT artifacts. Gamma pass rates of 85.7% and 89.9% were observed in two planar gamma analyses. In irradiations performed through the base of a dosimeter, gel response was found to increase with height in the dosimeter, even in areas of constant dose. After a correction was applied, gamma pass rates of 94.6% and 99.3% were observed in two planar gamma analyses. Absolute dose measurements were substantially higher (33%-100%) than the delivered doses for proton irradiations. Issues encountered while calibrating the LET-Meter gel restricted analysis of the LET measurement data to the SOBP of a proton beam. LET-Meter overresponse was found to increase linearly with track-average LET across the LET range that could be investigated (1.5 keV/micron – 3.5 keV/micron).

CHARACTERIZATION OF SELECTIVE INHIBITION OF ADENYLYL CYCLASE ACTIVITY BY SMALL MOLECULE INHIBITOR NKY80

The nine membrane-bound isoforms of adenylyl cyclase (AC), via synthesis of the signaling molecule cyclic AMP (cAMP), are involved in many isoform specific physiological functions. Decreasing AC5 activity has been shown to have potential therapeutic benefit, including reduced stress on the heart, pain relief, and attenuation of morphine dependence and withdrawal behaviors. However, AC structure is well conserved, and there are currently no isoform selective AC inhibitors in clinical use. P-site inhibitors inhibit AC directly at the catalytic site, but with an uncompetitive or noncompetitive mechanism. Due to this mechanism and nanomolar potency in cell-free systems, attempts at ligand-based drug design of novel AC inhibitors frequently use P-site inhibitors as a starting template. One small molecule inhibitor designed through this process, NKY80, is described as an AC5 selective inhibitor with low micromolar potency in vitro. P-site inhibitors reveal important ligand binding “pockets” in the AC catalytic site, but specific interactions that give NKY80 selectivity are unclear. Identifying and characterizing unique interactions between NKY80 and AC isoforms would significantly aid the development of isoform selective AC inhibitors. I hypothesized that NKY80’s selective inhibition is conferred by AC isoform specific interactions with the compound within the catalytic site. A structure-based virtual screen of the AC catalytic site was used to identify novel small molecule AC inhibitors. Identified novel inhibitors are isoform selective, supporting the catalytic site as a region capable of more potent isoform selective inhibition. Although NKY80 is touted commercially as an AC5 selective inhibitor, its characterization suggests strong inhibition of both AC5 and the closely related AC6. NKY80 was also virtually docked to AC to determine how NKY80 binds to the catalytic site. My results show a difference between NKY80 binding and the conformation of classic P-site inhibitors. The selectivity and notable differences in NKY80 binding to the AC catalytic site suggest a catalytic subregion more flexible in AC5 and AC6 that can be targeted by selective small molecule inhibitors.

DEVELOPMENT OF A BEAM-SPECIFIC PLANNING TARGET VOLUME AND A ROBUST PLAN ANALYSIS TOOL FOR PROTON THERAPY

Proton therapy is growing increasingly popular due to its superior dose characteristics compared to conventional photon therapy. Protons travel a finite range in the patient body and stop, thereby delivering no dose beyond their range. However, because the range of a proton beam is heavily dependent on the tissue density along its beam path, uncertainties in patient setup position and inherent range calculation can degrade thedose distribution significantly. Despite these challenges that are unique to proton therapy, current management of the uncertainties during treatment planning of proton therapy has been similar to that of conventional photon therapy. The goal of this dissertation research was to develop a treatment planning method and a planevaluation method that address proton-specific issues regarding setup and range uncertainties. Treatment plan designing method adapted to proton therapy: Currently, for proton therapy using a scanning beam delivery system, setup uncertainties are largely accounted for by geometrically expanding a clinical target volume (CTV) to a planning target volume (PTV). However, a PTV alone cannot adequately account for range uncertainties coupled to misaligned patient anatomy in the beam path since it does not account for the change in tissue density. In order to remedy this problem, we proposed a beam-specific PTV (bsPTV) that accounts for the change in tissue density along the beam path due to the uncertainties. Our proposed method was successfully implemented, and its superiority over the conventional PTV was shown through a controlled experiment.. Furthermore, we have shown that the bsPTV concept can be incorporated into beam angle optimization for better target coverage and normal tissue sparing for a selected lung cancer patient. Treatment plan evaluation method adapted to proton therapy: The dose-volume histogram of the clinical target volume (CTV) or any other volumes of interest at the time of planning does not represent the most probable dosimetric outcome of a given plan as it does not include the uncertainties mentioned earlier. Currently, the PTV is used as a surrogate of the CTV’s worst case scenario for target dose estimation. However, because proton dose distributions are subject to change under these uncertainties, the validity of the PTV analysis method is questionable. In order to remedy this problem, we proposed the use of statistical parameters to quantify uncertainties on both the dose-volume histogram and dose distribution directly. The robust plan analysis tool was successfully implemented to compute both the expectation value and its standard deviation of dosimetric parameters of a treatment plan under the uncertainties. For 15 lung cancer patients, the proposed method was used to quantify the dosimetric difference between the nominal situation and its expected value under the uncertainties.

Homology Cloning, Heterologous Expression and Characterization of a New Channelrhodopsin

Channelrhodopsins are phototaxis receptors in the plasma membranes of motile unicellular algae. They function as light-gated cation channels and this channel activity has been exploited to trigger action potentials in neurons with light to control neural circuits (“optogenetics"). Four channelrhodopsins were identified in two algal species, Chlamydomonas reinhardtii and Volvox carteri, with known genome sequences; each species contains 2 channelrhodopsins, one absorbing at longer wavelengths and one at shorter wavelengths, named CrChR1 and CrChR2, respectively. Our goals are to expand knowledge of channelrhodopsin mechanisms and also to identify new channelrhodopsins from various algal species with improved properties for optogenetic use. For these aims we are targeting algae from extreme environments to establish the natural diversity of their properties. We cloned a new channelrhodopsin from the psychrophilic (cold-loving) alga, Chlamydomonas augustae, with degenerate primers based on the 4 known homologs. The new protein is 48% and 52% identical to CrChR1 and CrChR2, respectively. We expressed the channelrhodopsin in HEK293 cells and measured light-induced currents to assess their kinetics and action spectrum. Based on the primary structure, kinetics of light-induced photocurrents in HEK293 cells, and action spectrum maximum of 520 nm near that of the two previously found CrChR1, we named the new channelrhodopsin CaChR1. The properties of robust channel activity at physiological pH, fast on-and-off kinetics, and greatly red-shifted action spectrum maximum from that of CrChR2, make CaChR1 advantageous as an optogenetic tool. To know this new channelrhodopsin better, we expressed His-tagged CaChR1 in Pichia pastoris and the yield is about 6 mg/L. The purified His-tagged CaChR1 exhibited an absorption spectrum identical to the action spectrum of CaChR1-generated photocurrents. The future work will be measurement of the photocycles of CaChR1 by flash photolysis, crystallization of CaChR1 for the structure and mutagenesis of CaChR1 to find the critical amino acids accounting for red-shifted spectra, slow inactivation and rapid on-and-off kinetics. Seven new channelrhodopsins including CaChR1 from different algal species have been cloned in our lab at this time, bringing the total known to 13. The work of cloning of these new channelrhodopsins along with the expression of CaChR1 was published in Photochemistry and Photobiology in January 2012

IMMUNOLOGICAL MECHANISMS OF EXTRACORPOREAL PHOTOPHERESIS IN CUTANEOUS T CELL LYMPHOMA AND GRAFT VERSUS HOST DISEASE

IMMUNOLOGICAL MECHANISMS OF EXTRACORPOREAL PHOTOPHERESIS IN CUTANEOUS T CELL LYMPHOMA AND GRAFT VERSUS HOST DISEASE Publication No.___________ Lisa Harn-Ging Shiue, B.S. Supervisory Professor: Madeleine Duvic, M.D. Extracorporeal photopheresis (ECP) is an effective, low-risk immunomodulating therapy for leukemic cutaneous T cell lymphoma (L-CTCL) and graft versus host disease (GVHD), but whether the mechanism(s) of action in these two diseases is (are) identical or different is unclear. To determine the effects of ECP in vivo, we studied regulatory T cells (T-regs), cytotoxic T lymphocytes (CTLs), and dendritic cells (DCs) by immunofluorescence flow cytometry in 18 L-CTCL and 11 GVHD patients before and after ECP at Day 2, 1 month, 3 months, and 6 months. In this study, ECP was effective in 12/18 L-CTCL patients with a 66.7% overall response rate (ORR) and 6/11 GVHD patients with a 54.5% ORR. Prior to ECP, the percentages of CD4+Foxp3+ T cells in 9 L-CTCL patients were either lower (L-CTCL-Low, n=2) or higher (L-CTCL-High, n=7) than normal. Five of the 7 GVHD patients had high percentages of CD4+Foxp3+ T cells (GVHD-High). Six of 7 L-CTCL-High patients had >80% CD4+Foxp3+ T cells which were correlated with tumor cells, and were responders. Both L-CTCL-High and GVHD-High patients had decreased percentages of CD4+Foxp3+ and CD4+Foxp3+CD25- T cells after 3 months of treatment. CD4+Foxp3+CD25+ T cells increased in GVHD-High patients but decreased in L-CTCL-High patients after 3 months of ECP. In addition, numbers of CTLs were abnormal. We confirmed that numbers of CTLs were low in L-CTCL patients, but high in GVHD patients prior to ECP. After ECP, CTLs increased after 1 month in 4/6 L-CTCL patients whereas CTLs decreased after 6 months in 3/3 GVHD patients. Myeloid (mDCs) and plasmacytoid DCs (pDCs) were also low at baseline in L-CTCL and GVHD patients confirming the DC defect. After 6 months of ECP, numbers and percentages of mDCs and pDCs increased in L-CTCL and GVHD. MDCs were favorably increased in 8/12 L-CTCL responders whereas pDCs were favorably increased in GVHD patients. These data suggest that ECP is favorably modulating the DC subsets. In L-CTCL patients, the mDCs may orchestrate Th1 cell responses to overcome immune suppression and facilitate disease regression. However, in GVHD patients, ECP is favorably down-regulating the immune system and may be facilitating immune tolerance to auto-or allo-antigens. In both L-CTCL and GVHD patients, DCs are modulated, but the T cell responses orchestrated by the DCs are different, suggesting that ECP modulates depending on the immune milieu. _______________

αvβ8 Integrin Interacts with RhoGDI1 to Regulate Rac1 and Cdc42 Activation and Drive Glioblastoma Cell Invasion

As an interface between the circulatory and central nervous systems, the neurovascular unit is vital to the development and survival of tumors. The malignant brain cancer glioblastoma multiforme (GBM) displays invasive growth behaviors that are major impediments to surgical resection and targeted therapies. Adhesion and signaling pathways that drive GBM cell invasion remain largely uncharacterized. Here we have utilized human GBM cell lines, primary patient samples, and pre-clinical mouse models to demonstrate that integrin αvβ8 is a major driver of GBM cell invasion. β8 integrin is overexpressed in many human GBM cells, with higher integrin expression correlating with increased invasion and diminished patient survival. Silencing β8 integrin in human GBM cells leads to impaired tumor cell invasion due to hyperactivation of the Rho GTPases Rac1 and Cdc42. β8 integrin associates with Rho GDP Dissociation Inhibitor 1 (RhoGDI1), an intracellular signaling effector that sequesters Rho GTPases in their inactive GDP-bound states. Silencing RhoGDI1 expression or uncoupling αvβ8 integrin-RhoGDI1 protein interactions blocks GBM cell invasion due to Rho GTPase hyperactivation. These data reveal for the first time that αvβ8 integrin, via interactions with RhoGDI1, suppresses activation of Rho proteins to promote GBM cell invasiveness. Hence, targeting the αvβ8 integrin-RhoGDI1 signaling axis may be an effective strategy for blocking GBM cell invasion.

Molecular Mechanisms of Vascular Disease in Patients with Rare Variants in MYH11

Thoracic aortic aneurysms and dissections (TAAD) are the primary disease affecting the thoracic ascending aorta, with an incidence rate of 10.4/100,000. Although about 20% of patients carry a mutation in a single gene that causes their disease, the remaining 80% of patients may also have genetic factors that increase their risk for developing TAAD. Many of the genes that predispose to TAAD encode proteins involved in smooth muscle cell (SMC) contraction and the disease-causing mutations are predicted to disrupt contractile function. SMCs are the predominant cell type in the ascending aortic wall. Mutations in MYH11, encoding the smooth muscle specific myosin heavy chain, are a rare cause of inherited TAAD. However, rare but recurrent non-synonymous variants in MYH11 are present in the general population but do not cause inherited TAAD. The goal of this study was to assess the potential role of these rare variants in vascular diseases. Two distinct variants were selected: the most commonly seen rare variant, MYH11 R247C, and a duplication of the chromosomal region spanning the MYH11 locus at 16p13.1. Genetic analyses indicated that both of these variants were significantly enriched in patients with TAAD compared with controls. A knock-in mouse model of the Myh11 R247C rare variant was generated, and these mice survive and reproduce normally. They have no structural abnormalities of the aorta or signs of aortic disease, but do have decreased aortic contractility. Myh11R247C/R247C mice also have increased proliferative response to vascular injury in vivo and increased proliferation of SMCs in vitro. Myh11R247C/R247C SMCs have decreased contractile gene and protein expression and are dedifferentiated. In fibroblasts, myosin force generation is required for maturation of focal adhesions, and enhancers of RhoA activity replace enhancers of Rac1 activity as maturation occurs. Consistent with these previous findings, focal adhesions are smaller in Myh11R247C/R247C SMCs, and there is decreased RhoA activation. A RhoA activator (CN03) rescues the dedifferentiated phenotype of Myh11R247C/R247C SMCs. Myh11R247C/R247C mice were bred with an existing murine model of aneurysm formation, the Acta2-/- mouse. Over time, mice carrying the R247C allele in conjunction with heterozygous or homozygous loss of Acta2 had significantly increased aortic diameter, and a more rapid accumulation of pathologic markers. These results suggest that the Myh11 R247C rare variant acts as a modifier gene increasing the risk for and severity of TAAD in mice. In patients with 16p13.1 duplications, aortic MYH11 expression is increased, but there is no corresponding increase in smooth muscle myosin heavy chain protein. Using SMCs that overexpress Myh11, we identified alterations in SMC phenotype leading to excessive protein turnover. All contractile proteins, not just myosin, are affected, and the proteins are turned over by autophagic degradation. Surprisingly, these cells are also more contractile compared with wild-type SMCs. The results described in this dissertation firmly establish that rare variants in MYH11 significantly affect the phenotype of SMCs. Further, the data suggests that these rare variants do increase the risk of TAAD via pathways involving altered SMC phenotype and contraction. Therefore, this study validates that these rare genetic variants alter vascular SMCs and provides model systems to explore the contribution of rare variants to disease.

Radiomics of NSCLC: Quantitative CT Image Feature Characterization and Tumor Shrinkage Prediction

Radiomics is the high-throughput extraction and analysis of quantitative image features. For non-small cell lung cancer (NSCLC) patients, radiomics can be applied to standard of care computed tomography (CT) images to improve tumor diagnosis, staging, and response assessment. The first objective of this work was to show that CT image features extracted from pre-treatment NSCLC tumors could be used to predict tumor shrinkage in response to therapy. This is important since tumor shrinkage is an important cancer treatment endpoint that is correlated with probability of disease progression and overall survival. Accurate prediction of tumor shrinkage could also lead to individually customized treatment plans. To accomplish this objective, 64 stage NSCLC patients with similar treatments were all imaged using the same CT scanner and protocol. Quantitative image features were extracted and principal component regression with simulated annealing subset selection was used to predict shrinkage. Cross validation and permutation tests were used to validate the results. The optimal model gave a strong correlation between the observed and predicted shrinkages with . The second objective of this work was to identify sets of NSCLC CT image features that are reproducible, non-redundant, and informative across multiple machines. Feature sets with these qualities are needed for NSCLC radiomics models to be robust to machine variation and spurious correlation. To accomplish this objective, test-retest CT image pairs were obtained from 56 NSCLC patients imaged on three CT machines from two institutions. For each machine, quantitative image features with concordance correlation coefficient values greater than 0.90 were considered reproducible. Multi-machine reproducible feature sets were created by taking the intersection of individual machine reproducible feature sets. Redundant features were removed through hierarchical clustering. The findings showed that image feature reproducibility and redundancy depended on both the CT machine and the CT image type (average cine 4D-CT imaging vs. end-exhale cine 4D-CT imaging vs. helical inspiratory breath-hold 3D CT). For each image type, a set of cross-machine reproducible, non-redundant, and informative image features was identified. Compared to end-exhale 4D-CT and breath-hold 3D-CT, average 4D-CT derived image features showed superior multi-machine reproducibility and are the best candidates for clinical correlation.

CHARACTERIZATION OF TCL1-Tg:P53-/- MICE THAT RESEMBLE HUMAN CHRONIC LYMPHOCYTIC LEUKEMIA WITH 17P-DELETION

Chronic lymphocytic leukemia (CLL) is the most common adult leukemia in the United Statesand Europe. CLL patients with deletion of chromosome 17p, where the tumor suppressor p53 gene is located, often develop a more aggressive disease with poor clinical outcomes. However, the underlying mechanism remains unclear. In order to understand the underneath mechanism in vivo, I have recently generated mice with Eu-TCL1-Tg:p53-/- genotype and showed that these mice develop aggressive leukemia that resembles human CLL with 17p deletion. The Eu-TCL1-Tg:p53-/- mice developed CLL disease at 3-4 months, significantly earlier than the parental Eu-TCL1-Tg mice that developed CLL disease at 8-12 months. Flow cytometry analysis showed that the CD5+/ IgM+ cell population appeared in the peritoneal cavity, bone marrow, and the spleens of Eu-TCL1-Tg:p53-/- mice significantly earlier than that of the parental Eu-TCL1-Tg mice. Massive infiltration and accumulation of leukemia cells were found in the spleen and peritoneal cavity. In vitro study showed that the leukemia cells isolated from the Eu-TCL1-Tg:p53-/- mice were more resistant to fludarabine treatment than the leukemia cells isolated from spleens of Eu-TCL1-Tg mice. Interestingly, TUNEL assay revealed that there was higher apoptotic cell death found in the Eu-TCL1-Tg spleen tissue compared to the spleens of the Eu-TCL1-Tg:p53-/- mice, suggesting that the loss of p53 compromises the apoptotic process in vivo, and this might in part explain the drug resistant phenotype of CLL cells with 17p-deletion. In the present study, we further demonstrated that the p53 deficiency in the TCL1 transgenic mice resulted in significant down-regulation of microRNAs miR-15a and miR16-1, associated with a substantial up-regulation of Mcl-1, suggesting that the p53-miR15a/16-Mcl-1 axis may play an important role in CLL pathogenesis. Interestingly, we also found that loss of p53 resulted in a significant decrease in expression of the miR-30 family especially miR-30d in leukemia lymphocytes from the Eu-TCL1-Tg:p53-/- mice. Such down-regulation of those microRNAs and up-regulation of Mcl-1 were also found in primary leukemia cells from CLL patients with 17p deletion. To further exam the biological significance of decrease in the miR-30 family in CLL, we investigated the potential involvement of EZH2 (enhancer of zeste homolog 2), a component of the Polycomb repressive complex known to be a downstream target of miR-30d and plays a role in disease progression in several solid cancers. RT-PCR and western blot analyses showed that both EZH2 mRNA transcript and protein levels were significantly increased in the lymphocytes of Eu-TCL1-Tg:p53-/- mice relative to Eu-TCL1-Tg mice. Exposure of leukemia cells isolated from Eu-TCL1-Tg:p53-/- mice to the EZH2 inhibitor 3-deazaneplanocin (DZNep) led to induction of apoptosis, suggesting EZH2 may play a role in promoting CLL cell survival and this may contribute to the aggressive phenotype of CLL with loss of p53. Our study has created a novel CLL mouse model, and suggests that the p53/miR15a/16-Mcl-1 axis & p53/miR30d-EZH2 may contribute to the aggressive phenotype and drug resistance in CLL cells with loss of p53.

OXIDATIVE STRESS BASED STRATEGIES FOR ENHANCING THE EFFICACY OF HISTONE DEACETYLASE INHIBITORS (HDACi)

Histone deacetylase inhibitors (HDACi) are anti-cancer drugs that primarily act upon acetylation of histones, however they also increase levels of intracellular reactive oxygen species (ROS). We hypothesized that agents that cause oxidative stress might enhance the efficacy of HDACi. To test this hypothesis, we treated acute lymphocytic leukemia cells (ALL) with HDACi and adaphostin (ROS generating agent). The combination of two different HDACi (vorinostat or entinostat) with adaphostin synergistically induced apoptosis in ALL. This synergistic effect was blocked when cells were pre-treated with the caspase-9 inhibitor, LEHD. In addition, we showed that loss of the mitochondrial membrane potential is the earliest event observed starting at 12 h. Following this event, we observed increased levels of superoxide at 16 h, and ultimately caspase-3 activation. Pre-treatment with the antioxidant N-acetylcysteine (NAC) blocked ROS generation and reversed the loss of mitochondrial membrane potential for both combinations. Interestingly, DNA fragmentation and caspase-3 activity was only blocked by NAC in cells treated with vorinostat-adaphostin; but not with entinostat-adaphostin. These results suggest that different redox mechanisms are involved in the induction of ROS-mediated apoptosis. To further understand these events, we studied the role of the antioxidants glutathione (GSH) and thioredoxin (Trx). We found that the combination of entinostat-adaphostin induced acetylation of the antioxidant thioredoxin (Trx) and decreased intracellular levels of GSH. However, no effect on Trx activity was observed in either combination. In addition, pre-treatment with GSH ethyl ester, a soluble form of GSH, did not block DNA fragmentation. Together these results suggested that GSH and Trx are not major players in the induction of oxidative stress. Array data examining the expression of genes involved in oxidative stress demonstrated a differential regulation between cells treated with vorinostat-adaphostin and entinostat-adaphostin. Some of the genes differentially expressed between the combinations include aldehyde oxidase 1, glutathione peroxidase-5, -6, peroxiredoxin 6 and myeloperoxidase. Taken together, these experimental results indicate that the synergistic activity of two different HDACi with adaphostin is mediated by distinct redox mechanisms in ALL cells. Understanding the mechanism involved in these combinations will advance scientific knowledge of how the action of HDACi could be augmented in leukemia models. Moreover, this information could be used for the development of effective clinical trials combining HDACi with other anticancer agents.

Improving Cervical Cancer Nodal Boost Radiation Therapy by Quantifying Uncertainties and Exploring Advanced Radiation Therapy Modalities

Radiation therapy for patients with intact cervical cancer is frequently delivered using primary external beam radiation therapy (EBRT) followed by two fractions of intracavitary brachytherapy (ICBT). Although the tumor is the primary radiation target, controlling microscopic disease in the lymph nodes is just as critical to patient treatment outcome. In patients where gross lymphadenopathy is discovered, an extra EBRT boost course is delivered between the two ICBT fractions. Since the nodal boost is an addendum to primary EBRT and ICBT, the prescription and delivery must be performed considering previously delivered dose. This project aims to address the major issues of this complex process for the purpose of improving treatment accuracy while increasing dose sparing to the surrounding normal tissues. Because external beam boosts to involved lymph nodes are given prior to the completion of ICBT, assumptions must be made about dose to positive lymph nodes from future implants. The first aim of this project was to quantify differences in nodal dose contribution between independent ICBT fractions. We retrospectively evaluated differences in the ICBT dose contribution to positive pelvic nodes for ten patients who had previously received external beam nodal boost. Our results indicate that the mean dose to the pelvic nodes differed by up to 1.9 Gy between independent ICBT fractions. The second aim is to develop and validate a volumetric method for summing dose of the normal tissues during prescription of nodal boost. The traditional method of dose summation uses the maximum point dose from each modality, which often only represents the worst case scenario. However, the worst case is often an exaggeration when highly conformal therapy methods such as intensity modulated radiation therapy (IMRT) are used. We used deformable image registration algorithms to volumetrically sum dose for the bladder and rectum and created a voxel-by-voxel validation method. The mean error in deformable image registration results of all voxels within the bladder and rectum were 5 and 6 mm, respectively. Finally, the third aim explored the potential use of proton therapy to reduce normal tissue dose. A major physical advantage of protons over photons is that protons stop after delivering dose in the tumor. Although theoretically superior to photons, proton beams are more sensitive to uncertainties caused by interfractional anatomical variations, and must be accounted for during treatment planning to ensure complete target coverage. We have demonstrated a systematic approach to determine population-based anatomical margin requirements for proton therapy. The observed optimal treatment angles for common iliac nodes were 90° (left lateral) and 180° (posterior-anterior [PA]) with additional 0.8 cm and 0.9 cm margins, respectively. For external iliac nodes, lateral and PA beams required additional 0.4 cm and 0.9 cm margins, respectively. Through this project, we have provided radiation oncologists with additional information about potential differences in nodal dose between independent ICBT insertions and volumetric total dose distribution in the bladder and rectum. We have also determined the margins needed for safe delivery of proton therapy when delivering nodal boosts to patients with cervical cancer.