Genetic susceptibility, cigarette smoking, and wood dust exposure in lung cancer: Case control analyses

This case control study was conducted to assess the association between lung cancer risk, mutagen sensitivity (a marker of cancer susceptibility), and a putative lung carcinogen, wood dust exposure. There were 165 cases (98 African-Americans, 67 Mexican-Americans) with newly diagnosed, previously untreated lung cancer, and 239 controls, frequency-matched on age, sex, and ethnicity.^ Mutagen sensitivity ($\ge$1 break/cell) was associated with a statistically significant elevated risk for lung cancer (odds ratio (OR) = 4.1, 95% confidence limits (CL) = 2.3,7.2). Wood dust exposure was also a significant predictor of risk (OR = 2.8, 95% CL = 1.2,6.6) after controlling for smoking and mutagen sensitivity. When stratified by ethnicity, wood dust exposure was a significant risk factor for African-Americans (OR = 4.0, 95% CL = 1.4,11.5), but not for Mexican-Americans (OR = 1.5, 95% CL = 0.3,7.1). Stratified analysis suggested a greater than multiplicative interaction between wood dust exposure and both mutagen sensitivity and smoking.^ The cases had significantly more breaks on chromosomes 4 and 5 than the controls did with ORs of 4.9 (95% CL = 2.0, 11.7) and 3.9 (95% CL = 1.6, 9.3), respectively. Breaks at 4p14, 4q27, 4q31, 5q21-22, 5q31, and 5q33 were significantly more common in lung cancer patients than in controls. Lung cancer risk had a dose-response relationship with breaks on chromosomes 4 and 5. Cigarette smoking had a strong interaction with breaks on chromosomes 2, 4, and 5.^ In a molecular cytogenetic study, using chromosome painting and G-banding, we showed that: (1) the proportion of chromosome 5 abnormalities surviving as chromosome-type aberrations remained significantly higher in cells of lung cancer cases (14%) than in controls (5%) (P $<$ 0.001). However, no significant differences were detected in chromosome 4 abnormalities between cases and controls; (2) the proportion of chromosome 5q13-22 abnormalities was 5.3% in the cases and 0.7% in the controls (P $<$ 0.001). 5q13-22 regions represented 40% of all abnormalities on chromosome 5 in the cases and only 14% in the controls.^ This study suggests that mutagen sensitivity, wood dust exposure, and cigarette smoking were independent risk factors for lung cancer, and the susceptibility of particular chromosome loci to mutagenic damage may be a genetic marker for specific types of lung cancer. ^

Studies of the expression of DNA repair related genes

This dissertation examines the biological functions and the regulation of expression of DNA ligase I by studying its expression under different conditions.^ The gene expression of DNA ligase I was induced two- to four-fold in S-phase lymphoblastoid cells but was decreased to 15% of control after administration of a DNA damaging agent, 4-nitroquinoline-1-oxide. When cells were induced into differentiation, the expression level of DNA ligase I was decreased to less than 15% of that of the control cells. When the gene of DNA ligase I was examined for tissue specific expression in adult rats, high levels of DNA ligase I mRNA were observed in testis (8-fold), intermediate levels in ovary and brain (4-fold), and low levels were found in intestine, spleen, and liver (1- to 2-fold).^ In confluent cells of normal skin fibroblasts, UV irradiation induced the gene expression of DNA ligase I at 24 and 48 h. The induction of DNA ligase I gene expression requires active p53 protein. Introducing a vector containing the wild type p53 protein in the cells caused an induction of the DNA ligase I protein 24 h after the treatment.^ Our results indicate that, in addition to the regulation by phosphorylation/dephosphorylation, cellular DNA ligase I activity can be regulated at the gene transcription level, and the p53 tumor suppresser is one of the transcription factors for the DNA ligase I gene. Also, our results suggest that DNA ligase I is involved in DNA repair as well as in DNA replication.^ Also, as an early attempt to clone the human homolog of the yeast CDC9 gene which has been shown to be involved in DNA replication, DNA repair, and DNA recombination, we have identified a human gene with mRNA of 1.7 kb. This dissertation studies the gene regulation and the possible biological functions of this new human gene by examining its expression at different stages of the cell cycle, during cell differentiation, and in cellular response to DNA damage.^ The new gene that we recently identified from human cells is highly expressed in brain and reproductive organs (BRE). This BRE gene encodes an mRNA of 1.7-1.9 kb, with an open reading frame of 1,149 bp, and gives rise to a deduced polypeptide of 383 amino acid residues. No extensive homology was found between BRE and sequences from the EMBL-Gene Banks. BRE showed tissue-specific expression in adult rats. The steady state mRNA levels were high in testis (5-6 fold), ovary and brain (3-4 fold) compared to the spleen level, but low in intestine and liver (1-2 fold). The expression of this gene is responsive to DNA damage and/or retinoic acid (RA) treatment. Treatment of fibroblast cells with UV irradiation and 4-nitroquinoline-1-oxide caused more than 90% and 50% decreases in BRE mRNA, respectively. Similar decreases in BRE expression were observed after treatment of the brain glioma cell line U-251 and the promyelocytic cell line HL-60 with retinoic acid. (Abstract shortened by UMI). ^

Biomarker circadian rhythm profiles in critically ill mechanically ventilated patients

Objective: To explore the natural trajectory of core body temperature (CBT) and cortisol (CORT) circadian rhythms in mechanically ventilated intensive care unit (MV ICU) patients. ^ Design: Prospective, observational, time-series pilot study. ^ Setting: Medical-surgical and pulmonary ICUs in a tertiary care hospital. ^ Sample: Nine (F = 3, M = 6) adults who were mechanically ventilated within 12 hrs of ICU admission with mean ± SD age of 65.2 ± 14 years old. ^ Measurements: Core body temperature and environmental measures of light, sound, temperature, and relative humidity were logged in 1-min intervals. Hourly urine specimens and 2-hr interval blood specimens were collected for up to 7 consecutive days for CORT assay. Mechanical ventilation days, ICU length of stay, and ICU mortality were documented. Acute Physiology and Chronic Health Evaluation (APACHE) II scores were computed for each study day. The data of each biologic and environmental variable were analyzed using single cosinor analysis of 24-hr serial segments. One patient did not complete the study because mortality occurred within 8 hrs of enrollment. Nine ICU patients completed the study in 1.6 to 7.0 days. ^ Results: No normal circadian rhythm pattern was found when the cosinor-derived parameters of amplitude (one-half the peak-trough variability) and acrophase (peak time) were compared with cosinor-derived parameter reference ranges of healthy, diurnally active humans, although 83% of patient-day CBT segments showed statistically significant (p ≤ .05) and biologically meaningful (R2≥ 0.30) 24-hr rhythms with abnormal cosinor parameters. Cosinor parameters of the environmental temporal profiles showed 27% of light, 76% of ambient temperature, and 78% of relative humidity serial segments had a significant and meaningful 24-hr diurnal pattern. Average daily light intensity varied from 34 to 187 lx with a maximum light exposure of 1877 lx. No sound measurement segment had a statistically significant cosine pattern, and numerous 1-minute interval peaks ≥ 60 dB occurred around the clock. Average daily ambient temperature and relative humidity varied from 19 to 24°C and from 25% to 61%, respectively. There was no statistically significant association between CBT or clinical outcomes and cosinor-derived parameters of the environmental variables. Circadian rhythms of urine and plasma CORT were deferred for later analysis. ^ Conclusions: The natural trajectory of the CBT circadian rhythm in MV ICU patients demonstrated persistent cosinor parameter alteration, even when a significant and meaningful 24-hr rhythm was present. The ICU environmental measures showed erratic light and sound exposures. Room temperature and relative humidity data produced the highest rate of significant and meaningful diurnal 24-hr patterns. Additional research is needed to clarify relations among the CBT biomarker of the circadian clock and environmental variables of MV ICU patients. ^

A case-control study of the association between corneal arcus and death by cardiovascular disease among cornea donors

The cornea of the human eye can develop deposits of lipids in the periphery known as corneal arcus. [2, 10] For over a century, these deposits have been of interest as possible indicators of the accumulation of lipids in arterial walls of the heart and body with implications for heart disease. [2, 10, 11, 29] Heart disease is currently the leading cause of death in this country. [5, 29] There have been several publications suggesting an association between the development of atherosclerotic lesions and corneal arcus. [2, 12, 29] Investigators have differed in their interpretation of the relevance of corneal arcus to coronary heart disease or cardiovascular disease. However, there is widespread consensus that the presence of corneal arcus in patients under the age of 50 should prompt physicians to further investigate for dyslipidemia or heart disease. [2, 3, 6, 8, 19] Earlier studies have often suffered from difficulty in determining the presence or severity of atherosclerosis and from inconsistencies in evaluating corneal arcus. This study involves the review of mortality data, medical and social history and standardized slit lamp examination of corneal tissue donors to evaluate the prevalence of corneal arcus in relation to death by CHD or CVD. The prevalence of arcus, odds ratio, and logistic regression was utilized for statistical analysis.^

CHEMOSENSITIZATION OF HEPATOCELLULAR CARCINOMA TO GEMCITABINE BY NON-INVASIVE RADIOFREQUENCY FIELD-INDUCED HYPERTHERMIA

Gemcitabine is a potent nucleoside analogue against solid tumors however drug resistance rapidly emerges. Removal of gemcitabine incorporated in the DNA by repair mechanisms could potentially contribute to resistance in chemo-refractory solid tumors. In this study, we evaluated homologous recombination repair of gemcitabine-stalled replication forks as a potential mechanism contributing to resistance. We also studied the effect of hyperthermia on homologous recombination pathway to explain the previously reported synergy between gemcitabine and hyperthermia. We found that hyperthermia degrades and inhibits localization of Mre11 to gemcitabine-stalled replication forks. Furthermore, gemcitabine-treated cells that were also treated with hyperthermia demonstrate a prolonged passage through late S/ G2 phase of cell cycle in comparison to cells treated with gemcitabine alone. This coincides with inhibition of resolution of γH2AX foci. Our findings also demonstrate that thermal sensitization of human hepatocellular carcinoma cell lines to gemcitabine is mediated through an Mre11-dependent homologous recombination repair pathway. Combination of non-invasive radiofrequency field-induced hyperthermia and gemcitabine was superior to either therapy alone (p

Genetic underpinnings of variation in brain ventricular volume: A genome-wide association study

The ventricular system is a critical component of the central nervous system (CNS) that is formed early in the developmental stages and remains functional through the lifetime. Changes in the ventricular system can be easily discerned via neuroimaging procedures and most of the time it reflects changes in the physiology of the CNS. In this study we attempted to identify specific genes associated with variation in ventricular volume in humans. Methods. We conducted a genome wide association (GWA) analysis of the volume of the lateral ventricles among 1605 individuals of European ancestry from two community based cohorts, the Genetics of Microangiopathic Brain Injury (GMBI; N=814) and Atherosclerosis Risk in Communities (ARIC; N=791). Significant findings from the analysis were tested for replication in both the cohorts and then meta-analyzed to get an estimate of overall significance. Results. In our GWA analyses, no single nucleotide polymorphism (SNP) reached a genome-wide significance of p<10−8. There were 25 SNPs in GMBI and 9 SNPs in ARIC that reached a threshold of p<10 −5. However, none of the top SNPs from each cohort were replicated in the other. In the meta-analysis, no SNP reached the genome-wide threshold of 5×10−8, but we identified five novel SNPs associated with variation in ventricular volume at the p<10 −5 level. Strongest association was for rs2112536 in an intergenic region on chromosome 5q33 (Pmeta= 8.46×10−7 ). The remaining four SNPs were located on chromosome 3q23 encompassing the gene for Calsyntenin-2 (CLSTN2). The SNPs with strongest association in this region were rs17338555 (Pmeta= 5.28×10 −6), rs9812091 (Pmeta= 5.89×10−6 ), rs9812283 (Pmeta= 5.97×10−6) and rs9833213 (Pmeta= 6.96×10−6). Conclusions. This GWA study of ventricular volumes in the community-based cohorts of European descent identifies potential locus on chromosomes 3 and 5. Further characterization of these loci may provide insights into pathophysiology of ventricular involvement in various neurological diseases.^

Metabolic syndrome in Mexican American females---prevalence and impact of age, socioeconomic status, and family history of type II diabetes

This thesis presents an analysis of data from Molecular Epidemiology of Type II Diabetes Mellitus in Mexican Americans. The study included 294 families. Among the participating families were 500 Mexican American females aged 19 to 86 who provided information on characteristics such as height, weight, and a variety of biochemical indicators. The research questions for this thesis are: (1) How strong is the association between indicators of the metabolic syndrome in study participants and their family histories of type II diabetes; and (2) How is an individual's family history of type II diabetes, age and socioeconomic status associated with the metabolic syndrome? In this thesis education status of the participants is used as an indicator of socioeconomic status. Answers to these questions are provided through the analysis of women's responses to written questionnaires and biochemical data. ^

Gene-gene interaction and gene-pathway analysis of genome-wide association for schizophrenia

Schizophrenia (SZ) is a complex disorder with high heritability and variable phenotypes that has limited success in finding causal genes associated with the disease development. Pathway-based analysis is an effective approach in investigating the molecular mechanism of susceptible genes associated with complex diseases. The etiology of complex diseases could be a network of genetic factors and within the genes, interaction may occur. In this work we argue that some genes might be of small effect that by itself are neither sufficient nor necessary to cause the disease however, their effect may induce slight changes to the gene expression or affect the protein function, therefore, analyzing the gene-gene interaction mechanism within the disease pathway would play crucial role in dissecting the genetic architecture of complex diseases, making the pathway-based analysis a complementary approach to GWAS technique. ^ In this study, we implemented three novel linkage disequilibrium based statistics, the linear combination, the quadratic, and the decorrelation test statistics, to investigate the interaction between linked and unlinked genes in two independent case-control GWAS datasets for SZ including participants of European (EA) and African (AA) ancestries. The EA population included 1,173 cases and 1,378 controls with 729,454 genotyped SNPs, while the AA population included 219 cases and 288 controls with 845,814 genotyped SNPs. We identified 17,186 interacting gene-sets at significant level in EA dataset, and 12,691 gene-sets in AA dataset using the gene-gene interaction method. We also identified 18,846 genes in EA dataset and 19,431 genes in AA dataset that were in the disease pathways. However, few genes were reported of significant association to SZ. ^ Our research determined the pathways characteristics for schizophrenia through the gene-gene interaction and gene-pathway based approaches. Our findings suggest insightful inferences of our methods in studying the molecular mechanisms of common complex diseases.^

Exploration of genetic susceptibility to pancreatic cancer: A GWAS approach

Pancreatic cancer is the 4th most common cause for cancer death in the United States, accompanied by less than 5% five-year survival rate based on current treatments, particularly because it is usually detected at a late stage. Identifying a high-risk population to launch an effective preventive strategy and intervention to control this highly lethal disease is desperately needed. The genetic etiology of pancreatic cancer has not been well profiled. We hypothesized that unidentified genetic variants by previous genome-wide association study (GWAS) for pancreatic cancer, due to stringent statistical threshold or missing interaction analysis, may be unveiled using alternative approaches. To achieve this aim, we explored genetic susceptibility to pancreatic cancer in terms of marginal associations of pathway and genes, as well as their interactions with risk factors. We conducted pathway- and gene-based analysis using GWAS data from 3141 pancreatic cancer patients and 3367 controls with European ancestry. Using the gene set ridge regression in association studies (GRASS) method, we analyzed 197 pathways from the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Using the logistic kernel machine (LKM) test, we analyzed 17906 genes defined by University of California Santa Cruz (UCSC) database. Using the likelihood ratio test (LRT) in a logistic regression model, we analyzed 177 pathways and 17906 genes for interactions with risk factors in 2028 pancreatic cancer patients and 2109 controls with European ancestry. After adjusting for multiple comparisons, six pathways were marginally associated with risk of pancreatic cancer ( P < 0.00025): Fc epsilon RI signaling, maturity onset diabetes of the young, neuroactive ligand-receptor interaction, long-term depression (Ps < 0.0002), and the olfactory transduction and vascular smooth muscle contraction pathways (P = 0.0002; Nine genes were marginally associated with pancreatic cancer risk (P < 2.62 × 10−5), including five reported genes (ABO, HNF1A, CLPTM1L, SHH and MYC), as well as four novel genes (OR13C4, OR 13C3, KCNA6 and HNF4 G); three pathways significantly interacted with risk factors on modifying the risk of pancreatic cancer (P < 2.82 × 10−4): chemokine signaling pathway with obesity ( P < 1.43 × 10−4), calcium signaling pathway (P < 2.27 × 10−4) and MAPK signaling pathway with diabetes (P < 2.77 × 10−4). However, none of the 17906 genes tested for interactions survived the multiple comparisons corrections. In summary, our current GWAS study unveiled unidentified genetic susceptibility to pancreatic cancer using alternative methods. These novel findings provide new perspectives on genetic susceptibility to and molecular mechanisms of pancreatic cancer, once confirmed, will shed promising light on the prevention and treatment of this disease. ^

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

BIOCHEMICAL CHARACTERIZATION OF BINDING PARTNERS OF TWO HSP70 CO-CHAPERONES IN SACCHAROMYCES CEREVISIAE

Cells are exposed to a variety of environmental and physiological changes including temperature, pH and nutrient availability. These changes cause stress to cells, which results in protein misfolding and altered cellular protein homeostasis. How proteins fold into their three-dimensional functional structure is a fundamental biological process with important relevance to human health. Misfolded and aggregated proteins are linked to multiple neurodegenerative diseases, cardiovascular disease and cystic fibrosis. To combat proteotoxic stress, cells deploy an array of molecular chaperones that assist in the repair or removal of misfolded proteins. Hsp70, an evolutionarily conserved molecular chaperone, promotes protein folding and helps maintain them in a functional state. Requisite co-chaperones, including nucleotide exchange factors (NEFs) strictly regulate and serve to recruit Hsp70 to distinct cellular processes or locations. In yeast and human cells, three structurally non-related cytosolic NEFs are present: Sse1 (Hsp110), Fes1 (HspBP1) and Snl1 (Bag-1). Snl1 is unique among the cytosolic NEFs as it is localized at the ER membrane with its Hsp70 binding (BAG) domain exposed to the cytosol. I discovered that Snl1 distinctly interacts with assembled ribosomes and several lines of evidence indicate that this interaction is both independent of and concurrent with binding to Hsp70 and is not dependent on membrane localization. The ribosome-binding site is identified as a short lysine-rich motif within the amino terminus of the Snl1 BAG domain distinct from the Hsp70 interaction region. In addition, I demonstrate ribosome association with the Snl1 homolog in the pathogenic fungus, Candida albicans and localize this putative NEF to a perinuclear/ER membrane, suggesting functional conservation in fungal BAG domain-containing proteins. As a first step in determining specific domain architecture in fungal BAG proteins, I present the preliminary steps of protein purification and analysis of the minimal Hsp70 binding region in in both S.cerevisiae and C. albicans Snl1. Contrary to previous in vitro evidence which showed the Fes1 NEF to interact with both cytosolic Hsp70s, Ssa and Ssb, Fes1 is shown to interact specifically with Ssa when expressed under normal cellular conditions in S. cerevisiae. This is the first reported evidence of Hsp70 binding selectivity for a cytosolic NEF, and suggests a possible mechanism to achieve specificity in Hsp70-dependent functions. Taken together, the work presented in this dissertation highlights the striking divergence among Hsp70 co-chaperones in selecting binding partners, which may correlate with their specific roles in the cell.

The Role of the Arched Helicases in Exosome-Mediated Function

RNA processing and degradation are two important functions that control gene expression and promote RNA fidelity in the cell. A major ribonuclease complex, called the exosome, is involved in both of these processes. The exosome is composed of ten essential proteins with only one catalytically active subunit, called Rrp44. While the same ten essential subunits make up both the nuclear and cytoplasmic exosome, there are nuclear and cytoplasmic exosome cofactors that promote specific exosome functions in each of the cell compartments. To date, it is unclear how the exosome distinguishes between RNA substrates. We hypothesize that compartment specific cofactors may promote the substrate specificity of the exosome. In this work, I characterize several cofactors of the exosome, both nuclear and cytoplasmic. First, I describe the arch domain, which is a unique domain in a nuclear and a cytoplasmic cofactor of the exosome. Specifically, I show that the arch domain of the nuclear exosome cofactor, Mtr4, is required for specific exosome-mediated activities and overlaps functionally with the exosome-associated exonuclease, Rrp6. Further, I show that the arch domain of Ski2 is required for the degradation of normal and aberrant mRNAs. Additionally, this work describes in detail the Mtr4 domains involved in the physical association with other RNA processing proteins. Further, I characterize the minimal Mtr4-binding region in a third exosome cofactor, Trf5. Understanding how exosome cofactors synergistically promote exosome function will provide us a better understanding of how the exosome complex precisely regulates its catalytic activities. As described here, cofactors play a major role in determining the substrate specificity of the nuclear and cytoplasmic exosome. Moreover, specific accessory domains, which are not involved in the catalytic function of the cofactor, are required for substrate targeting of the eukaryotic RNA exosome.

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.

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.

Interview with Julie Knobil

An oral interview with Juile (Hotchkiss) Knobil, research professor of physiology and then integrative biology at the Medical School, where she lectured on mammalian physiology and perinatal endocrinology.