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.

Identification and analysis of novel mitotic inhibitors of the Caenorhabditis elegans Aurora B kinase

Proper execution of mitosis requires the accurate segregation of replicated DNA into each daughter cell. The highly conserved mitotic kinase AIR-2/Aurora B is a dynamic protein that interacts with subsets of cofactors and substrates to coordinate chromosome segregation and cytokinesis in Caenorhabdiris elegans. To identify components of the AIR-2 regulatory pathway, a genome-wide RNAi-based screen for suppressors of air-2 temperature-sensitive mutant lethality was conducted. Here, I present evidence that two classes of suppressors identified in this screen are bona fide regulators of the AIR-2 kinase. The strongest suppressor cdc-48.3, encodes an Afg2/Spaf-related Cdc48-like AAA+ ATPase that regulates AIR-2 kinase activity and stability during C. elegans embryogenesis. Loss of CDC-48.3 suppresses the lethality of air-2 mutant embryos, marked by the restoration of the dynamic behavior of AIR-2 and rescue of chromosome segregation and cytokinesis defects. Loss of CDC-48.3 leads to mitotic delays and abnormal accumulation of AIR-2 during late telophase/mitotic exit. In addition, AIR-2 kinase activity is significantly upregulated from metaphase through mitotic exit in CDC-48.3 depleted embryos. Inhibition of the AIR-2 kinase is dependent on (1) a direct physical interaction between CDC-48.3 and AIR-2, and (2) CDC-48.3 ATPase activity. Importantly, the increase in AIR-2 kinase activity does not correlate with the stabilization of AIR-2 in late mitosis. Hence, CDC-48.3 is a bi-functional inhibitor of AIR-2 that is likely to act via distinct mechanisms. The second class of suppressors consists of psy-2/smk-1 and pph-4.1, which encode two components of the conserved PP4 phosphatase complex that is essential for spindle assembly, chromosome segregation, and overall mitotic progression. AIR-2 and its substrates are likely to be targets of this complex since mitotic AIR-2 kinase activity is significantly increased during mitosis when either PSY-2/SMK-1 or PPH-4.l is depleted. Altogether, this study demonstrates that during the C. elegans embryonic cell cycle, regulators including the CDC-48.3 ATPase and PP4 phosphatase complex interact with and control the kinase activity, targeting behavior and protein stability of the Aurora B kinase to ensure accurate and timely progression of mitosis. ^

EphA2 expression is regulated by EGFR and KRAS and promotes non-small cell lung cancer progression

Lung cancer is the leading cause of cancer deaths worldwide. The development of improved systemic therapy is needed for the most common form of the disease, non-small cell lung cancer (NSCLC). This will depend on the identification of valid molecular targets. Recent studies point to the receptor tyrosine kinase EphA2 as a novel therapeutic target. Overexpression of EphA2 has been demonstrated in a number of epithelial cancers, and its expression has been associated with more severe disease. Regulation of EphA2 in cancer is poorly understood. Recently, regulation of EphA2 by EGFR and KRAS has been reported in a number of in vitro models, but no examination of this relationship has been undertaken in patient tumors. Because of the established importance of EGFR and KRAS in NSCLC, we have investigated the relationship between these mutations and EphA2 in NSCLC patient tissues and cell lines. The significance of Epha2 expression was further examined by testing for correlation with survival, metastases, histology, and smoking status in patient tissues, and tumor cell proliferation and migration in vitro. EphA2 expression was analyzed in by immunohistochemistry in tissue microarray (TMA) format utilizing surgically resected lung cancer specimens. EGFR and KRAS mutation status was determined for the majority of specimens. EphA2 expression was detected in >90% of NSCLC tumors. High EphA2 expression was associated with decreased time to recurrence and metastases, and predicted poorer progression free and overall survival. Expression of EphA2 was positively correlated with activated EGFR and with KRAS mutation. Expression of EphA2 was also positively correlated with a history of smoking. There was no association between gender or histology and EphA2 expression. In H322 cells, activation of EGFR or KRAS resulted in an increase in EphA2 protein expression. Downregulation of EphA2 resulted in decreased proliferation in a clonal growth assay, and inhibited migration in a wound healing assay, in a panel of cell lines. The decrease in proliferation correlated with a transient decrease in the levels of phospho-ERK, a downstream effector of EGFR and KRAS. Based on these data, the potential of EphA2 as a therapeutic target for NSCLC should be further investigated. ^

Involvement of lipocalin 2 in leukemia and breast cancer

Human lipocalin 2 is described as the neutrophil gelatinase-associated lipocalin (NGAL). The lipocalin 2 gene encodes a small, secreted glycoprotein that possesses a variety of functions, of which the best characterized function is organic iron binding activity. Elevated NGAL expression has been observed in many human cancers including breast, colorectal, pancreatic and ovarian cancers. I focused on the characterization of NGAL function in chronic myelogenous leukemia (CML) and breast cancer. Using the leukemic xenograft mouse model, we demonstrated that over-expression of NGAL in K562 cells, a leukemic cell line, led to a higher apoptotic rate and an atrophy phenotype in the spleen of inoculated mice compared to K562 cells alone. These results indicate that NGAL plays a primary role in suppressing hematopoiesis by inducing apoptosis within normal hematopoietic cells. In the breast cancer project, we analyzed two microarray data sets of breast cancer cell lines ( n = 54) and primary breast cancer samples (n = 318), and demonstrated that high NGAL expression is significantly correlated with several tumor characteristics, including negative estrogen receptor (ER) status, positive HER2 status, high tumor grade, and lymph node metastasis. Ectopic NGAL expression in non-aggressive (ZR75.1 and MCF7) cells led to aggressive tumor phenotypes in vitro and in vivo. Conversely, knockdown of NGAL expression in various breast cancer cell lines by shRNA lentiviral infection significantly decreased migration, invasion, and metastasis activities of tumor cells both in vitro and in vivo . It has been previously reported that transgenic mice with a mutation in the region of trans-membrane domain (V664E) of HER2 develop mammary tumors that progress to lung metastasis. However, we observed that genetic deletion of the 24p3 gene, a mouse homolog of NGAL, in HER2 transgenic mice by breeding with 24p3-null mice resulted in a significant delay of mammary tumor formation and reduction of lung metastasis. Strikingly, we also found that treatment with affinity purified 24p3 antibodies in the 4T1 breast cancer mice strongly reduced lung metastasis. Our studies provide evidence that NGAL plays a critical role in breast cancer development and progression, and thus NGAL has potential as a new therapeutic target in breast cancer.^

Novel Imaging-Based Techniques Reveal a Role for PD-1/PD-L1 in Tumor Immune Surveillance in the Lung

The binding of immune inhibitory receptor Programmed Death 1 (PD-1) on T cells to its ligand PD-L1 has been implicated as a major contributor to tumor induced immune suppression. Clinical trials of PD-L1 blockade have proven effective in unleashing therapeutic anti-tumor immune responses in a subset of patients with advanced melanoma, yet current response rates are low for reasons that remain unclear. Hypothesizing that the PD-1/PD-L1 pathway regulates T cell surveillance within the tumor microenvironment, we employed intravital microscopy to investigate the in vivo impact of PD-L1 blocking antibody upon tumor-associated immune cell migration. However, current analytical methods of intravital dynamic microscopy data lack the ability to identify cellular targets of T cell interactions in vivo, a crucial means for discovering which interactions are modulated by therapeutic intervention. By developing novel imaging techniques that allowed us to better analyze tumor progression and T cell dynamics in the microenvironment; we were able to explore the impact of PD-L1 blockade upon the migratory properties of tumor-associated immune cells, including T cells and antigen presenting cells, in lung tumor progression. Our results demonstrate that early changes in tumor morphology may be indicative of responsiveness to anti-PD-L1 therapy. We show that immune cells in the tumor microenvironment as well as tumors themselves express PD-L1, but immune phenotype alone is not a predictive marker of effective anti-tumor responses. Through a novel method in which we quantify T cell interactions, we show that T cells are largely engaged in interactions with dendritic cells in the tumor microenvironment. Additionally, we show that during PD-L1 blockade, non-activated T cells are recruited in greater numbers into the tumor microenvironment and engage more preferentially with dendritic cells. We further show that during PD-L1 blockade, activated T cells engage in more confined, immune synapse-like interactions with dendritic cells, as opposed to more dynamic, kinapse-like interactions with dendritic cells when PD-L1 is free to bind its receptor. By advancing the contextual analysis of anti-tumor immune surveillance in vivo, this study implicates the interaction between T cells and tumor-associated dendritic cells as a possible modulator in targeting PD-L1 for anti-tumor immunotherapy.