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.

3.88 A structure of cytoplasmic polyhedrosis virus by cryo-electron microscopy.

Cytoplasmic polyhedrosis virus (CPV) is unique within the Reoviridae family in having a turreted single-layer capsid contained within polyhedrin inclusion bodies, yet being fully capable of cell entry and endogenous RNA transcription. Biochemical data have shown that the amino-terminal 79 residues of the CPV turret protein (TP) is sufficient to bring CPV or engineered proteins into the polyhedrin matrix for micro-encapsulation. Here we report the three-dimensional structure of CPV at 3.88 A resolution using single-particle cryo-electron microscopy. Our map clearly shows the turns and deep grooves of alpha-helices, the strand separation in beta-sheets, and densities for loops and many bulky side chains; thus permitting atomic model-building effort from cryo-electron microscopy maps. We observed a helix-to-beta-hairpin conformational change between the two conformational states of the capsid shell protein in the region directly interacting with genomic RNA. We have also discovered a messenger RNA release hole coupled with the mRNA capping machinery unique to CPV. Furthermore, we have identified the polyhedrin-binding domain, a structure that has potential in nanobiotechnology applications.

High Resolution Vascular Imaging of the Rat Spine using Liposomal Blood Pool MR Agent

High resolution, vascular magnetic resonance imaging of the spine region in small animals poses several challenges. The small anatomical features, extravascular diffusion, and the low signal-to-noise ratio limit the use of conventional contrast agents. We hypothesize that a long circulating, intravascular liposomal-encapsulated MR contrast agent (liposomal-Gd) would facilitate visualization of small anatomical features of the perispinal vasculature not visible with conventional contrast agent (Gd-DTPA).

High-resolution vascular imaging of the rat spine using liposomal blood pool MR agent.

BACKGROUND AND PURPOSE: High-resolution, vascular MR imaging of the spine region in small animals poses several challenges. The small anatomic features, extravascular diffusion, and low signal-to-noise ratio limit the use of conventional contrast agents. We hypothesize that a long-circulating, intravascular liposomal-encapsulated MR contrast agent (liposomal-Gd) would facilitate visualization of small anatomic features of the perispinal vasculature not visible with conventional contrast agent (gadolinium-diethylene-triaminepentaacetic acid [Gd-DTPA]). METHODS: In this study, high-resolution MR angiography of the spine region was performed in a rat model using a liposomal-Gd, which is known to remain within the blood pool for an extended period. The imaging characteristics of this agent were compared with those of a conventional contrast agent, Gd-DTPA. RESULTS: The liposomal-Gd enabled acquisition of high quality angiograms with high signal-to-noise ratio. Several important vascular features, such as radicular arteries, posterior spinal vein, and epidural venous plexus were visualized in the angiograms obtained with the liposomal agent. The MR angiograms obtained with conventional Gd-DTPA did not demonstrate these vessels clearly because of marked extravascular soft-tissue enhancement that obscured the vasculature. CONCLUSIONS: This study demonstrates the potential benefit of long-circulating liposomal-Gd as a MR contrast agent for high-resolution vascular imaging applications.

Structural studies of proteinase entrapment by human alpha 2-macroglobulin using 3-D electron microscopy

Human a2 -macroglobulin ( a2 M; homotetramer, Mr 720 kDa) is an essential scavenger of proteinases in the serum. Each of its four subunits has a ‘bait region’, with cleavage sequences for almost all endo-proteinases, an unusual thiol ester moiety and a receptor-binding domain (RBD). Bait region cleavage in native a2 M ( a2 M-N) by a proteinase results in rapid thiol ester breakage, with a large-scale structural transformation, in which a2 M uniquely entraps the proteinase in a cage-like structure and exposes receptor-binding domains for rapid endocytosis. Transformed a2 M ( a2 M-TR) contains up to two proteinases, which remain active to small substrates. 3-D electron microscopy is optimally suited to study this unusual structural change at resolutions near (1/30) Å−1. ^ The structural importance of the thiol esters was demonstrated by a genetically-engineered a2 M, with the cysteines involved in thiol ester formation mutated to serines, which appeared structurally homologous to a2 M-TR. This demonstrates that the four highly labile thiol esters alone maintain the a2 M-N structure, while the ‘closed trap’ formed by a2 M-TR is a more stable structural form. ^ Half-transformed a2 M ( a2 M-HT), with cleaved bait regions and thiol esters in only two of its four subunits, provides an important structural link between a2 M-N and a2 M-TR. A comparison with a2 M-N showed the two proteinase-entrapping domains were above and below the plane bisecting the long axis. Both a2 M-N and a2 M-TR consist of two dense, oppositely twisted strands with significant interconnections, indicating that the structural change involves a rotation of these strands. In a2 M-HT these strands were partially untwisted with large central openings, revealing the manner in which the proteinase enters the internal cavity of a2 M. ^ In reconstructions of a2 M-N, a2 M-HT and a2 M-TR labeled with a monoclonal Fab, the Fabs were located on distal ends of each constitutive strand, demonstrating an anti-parallel arrangement of the subunits. Separation between the top and bottom pairs of Fabs was nearly the same on all structures, but the pairs were rotated about the long axis. Taken together, these results indicate that upon proteinase cleavage the two strands in a2 M-N separate. The proteinase enters the structure, while the strands re-twist to encage it. In a2 M-TR, which displays receptor-binding arms, more than two subunits are transformed as strands in the transformed half of a2 M-HT were not separated. ^

High-resolution microarray analysis of chromosome 20q in human colon cancer metastasis model systems

Amplification of human chromosome 20q DNA is the most frequently occurring chromosomal abnormality detected in sporadic colorectal carcinomas and shows significant correlation with liver metastases. Through comprehensive high-resolution microarray comparative genomic hybridization and microarray gene expression profiling, we have characterized chromosome 20q amplicon genes associated with human colorectal cancer metastasis in two in vitro metastasis model systems. The results revealed increasing complexity of the 20q genomic profile from the primary tumor-derived cell lines to the lymph node and liver metastasis derived cell lines. Expression analysis of chromosome 20q revealed a subset of over expressed genes residing within the regions of genomic copy number gain in all the tumor cell lines, suggesting these are Chromosome 20q copy number responsive genes. Bases on their preferential expression levels in the model system cell lines and known biological function, four of the over expressed genes mapping to the common intervals of genomic copy gain were considered the most promising candidate colorectal metastasis-associated genes. Validation of genomic copy number and expression array data was carried out on these genes, with one gene, DNMT3B, standing out as expressed at a relatively higher levels in the metastasis-derived cell lines compared with their primary-derived counterparts in both the models systems analyzed. The data provide evidence for the role of chromosome 20q genes with low copy gain and elevated expression in the clonal evolution of metastatic cells and suggests that such genes may serve as early biomarkers of metastatic potential. The data also support the utility of the combined microarray comparative genomic hybridization and expression array analysis for identifying copy number responsive genes in areas of low DNA copy gain in cancer cells. ^

Intercellular adhesion and membrane polarity in rat hepatocytes: Localization of cell -CAM 105 and other domain-specific membrane proteins {\it in situ\/} and {\it in vitro\/} by electron microscopy

Cell-CAM 105 has been identified as a cell adhesion molecule (CAM) based on the ability of monospecific and monovalent anti-cell-CAM 105 antibodies to inhibit the reaggregation of rat hepatocytes. Although one would expect to find CAMs concentrated in the lateral membrane domain where adhesive interactions predominate, immunofluorescence analysis of rat liver frozen sections revealed that cell-CAM 105 was present exclusively in the bile canalicular (BC) domain of the hepatocyte. To more precisely define the in situ localization of cell-CAM 105, immunoperoxidase and electron microscopy were used to analyze intact and mechanically dissociated fixed liver tissue. Results indicate that although cell-CAM 105 is apparently restricted to the BC domain in situ, it can be detected in the pericanalicular region of the lateral membranes when accessibility to lateral membranes is provided by mechanical dissociation. In contrast, when hepatocytes were labeled following incubation in vitro under conditions used during adhesion assays, cell-CAM 105 had redistributed to all areas of the plasma membrane. Immunofluorescence analysis of primary hepatocyte cultures revealed that cell-CAM 105 and two other BC proteins were localized in discrete domains reminscent of BC while cell-CAM 105 was also present in regions of intercellular contact. These results indicate that the distribution of cell-CAM 105 under the experimental conditions used for cell adhesion assays differs from that in situ and raises the possibility that its adhesive function may be modulated by its cell surface distribution. The implications of these and other findings are discussed with regard to a model for BC formation.^ Analysis of molecular events involved in BC formation would be accelerated if an in vitro model system were available. Although BC formation in culture has previously been observed, repolarization of cell-CAM 105 and two other domain-specific membrane proteins was incomplete. Since DMSO had been used by Isom et al. to maintain liver-specific gene expression in vitro, the effect of this differentiation system on the polarity of these membrane proteins was examined. Based on findings presented here, DMSO apparently prolongs the expression and facilitates polarization of hepatocyte membrane proteins in vitro. ^