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.

Structural conservations and diversities of dsRNA viruses: Structural studies of the cytoplasmic polyhedrosis virus by electron cryomicroscopy and 3D reconstruction

The Reoviridae virus family is a group of economically and pathologically important viruses that have either single-, double-, or triple-shelled protein layers enclosing a segmented double stranded RNA genome. Each virus particle in this family has its own viral RNA dependent RNA polymerase and the enzymatic activities necessary for the mature RNA synthesis. Based on the structure of the inner most cores of the viruses, the Reoviridae viruses can be divided into two major groups. One group of viruses has a smooth surfaced inner core, surrounded by complete outer shells of one or two protein layers. The other group has an inner core decorated with turrets on the five-fold vertices, and could either completely lack or have incomplete outer protein layers. The structural difference is one of the determinant factors for their biological differences during the infection. ^ Cytoplasmic polyhedrosis virus (CPV) is a single-shelled, turreted virus and the structurally simplest member in Reoviridae. It causes specific chronic infections in the insect gut epithelial cells. Due to its wide range of insect hosts, CPV has been engineered as a potential insecticide for use in fruit and vegetable farming. Its unique structural simplicity, unparalleled capsid stability and ease of purification make CPV an ideal model system for studying the structural basis of dsRNA virus assembly at the highest possible resolution by electron cryomicroscopy (cryoEM) and three-dimensional (3D) reconstruction. ^ In this thesis work, I determined the first 3D structure of CPV capsids using 100 kV cryoEM. At an effective resolution of 17 Å, the full capsid reveals a 600-Å diameter, T = 1 icosahedral shell decorated with A and B spikes at the 5-fold vertices. The internal space of the empty CPV is unoccupied except for 12 mushroom-shaped densities that are attributed to the transcriptional enzyme complexes. The inside of the full capsid is packed with icosahedrally-ordered viral genomic RNA. The interactions of viral RNA with the transcriptional enzyme complexes and other capsid proteins suggest a mechanism for RNA transcription and subsequent release. ^ Second, the interactions between the turret proteins (TPs) and the major capsid shell protein (CSPs) have been identified through 3D structural comparisons of the intact CPV capsids with the spikeless CPV capsids, which were generated by chemical treatments. The differential effects of these chemical treatment experiments also indicated that CPV has a significantly stronger structural integrity than other dsRNA viruses, such as the orthoreovirus subcores, which are normally enclosed within outer protein shells. ^ Finally, we have reconstructed the intact CPV to an unprecendented 8 Å resolution from several thousand of 400kV cryoEM images. The 8 Å structure reveals interactions among the 120 molecules of each of the capsid shell protein (CSP), the large protrusion protein (LPP), and 60 molecules of the turret protein (TP). A total of 1980 α-helices and 720 β-sheets have been identified in these capsid proteins. The CSP structure is largely conserved, with the majority of the secondary structures homologous to those observed in the x-ray structures of corresponding proteins of other reoviruses, such as orthoreovirus and bluetongue virus. The three domains of TP are well positioned to play multifunctional roles during viral transcription. The completely non-equivalent interactions between LPP and CSP and those between the anchoring domain of TP and CSP account for the unparalleled stability of this structurally simplest member of the Reoviridae. ^

Geochemical analyses of secondary minerals from ODP Leg 115 basalts (Tables 1+2)

Basalts recovered along the Reunion hotspot track on Ocean Drilling Program (ODP) Leg 115 range in age from 34 Ma at Site 706 to 64 Ma at Site 707. They have undergone various degrees of secondary alteration. Within single holes the amount of alteration can vary from a few percent to near complete replacement of phenocrysts and groundmass by secondary minerals. Olivine appears to be the most susceptible to alteration and in some sections it is the only mineral altered. In other sections, olivine, pyroxene and plagioclase phenocrysts, and groundmass have been completely replaced by secondary minerals. Clays are the predominant form of secondary mineralization. In addition to replacing olivine, pyroxene, glass, and groundmass, clays have filled veins, vesicles, and voids. Minor amounts of calcite, zeolites, and K-feldspar were also detected. The clays that filled vesicles and veins often show color zonations of dark, opaque bands near the edges that grade into tan or green transparent regions in the centers of the veins. The electron microprobe was used to obtain chemical analyses of these veins as well as to characterize isolated clays that replaced specific minerals and filled voids and vesicles.

Mineralogy and chemistry of basalts and secondary minerals from DSDP Site 417, Atlantic Ocean

Basalts from DSDP Site 417 (109 Ma) exhibit the effects of several stages of alteration reflecting the evolution of seawater-derived solution compositions and control by the structure and permeability of the crust. Characteristic secondary mineral assemblages occur in often superimposed alteration zones within individual basalt fragments. By combining bulk rock and single phase chemical analyses with detailed mineralogic and petrographic studies, chemical changes have been determined for most of the alteration stages identified in the basalts. 1) Minor amounts of saponite, chlorite, and pyrite formed locally in coarse grained portions of massive units, possibly at high temperatures during initial cooling of the basalts. No chemical changes could be determined for this stage. 2) Possible mixing of cooled hydrothermal fluids with seawater resulted in the formation of celadonite-nontronite and Fe-hydroxide-rich black halos around cracks and pillow rims. Gains of K, Rb, H20, increase of Fe 3 +/FeT and possibly some losses of Ca and Mg occurred during this stage. 3a) Extensive circulation of oxygenated seawater resulted in the formation of various smectites, K-feldspar, and Fe-hydroxides in brown and light grey alteration zones around formerly exposed surfaces. K, Rb, H20, and occasionally P were added to the rocks, Fe3+/FeT increased, and Ca, Mg, Si and occasionally Al and Na were lost. 3 b) Anoxic alteration occurred during reaction of basalt with seawater at low water-rock ratios, or with seawater that had previously reacted with basalt. Saponite-rich dark grey alteration zones formed which exhibit very little chemical change: generally only slight increases in Fe 3 +/FeT and H20 occurred. 4) Zeolites and calcite formed from seawater-derived fluids modified by previous reactions with basalt. Chemical changes involved increases of Ca, Na, H20 , and CO2 in the rocks. 5) A late stage of anoxic conditions resulted in the formation of minor amounts of Mn-calcites and secondary sulfides in previously oxidized rocks. No chemical changes were determined for this stage. Recognition of such alteration sequences is important in understanding the evolution of submarine hydrothermal systems and in interpreting chemical exchange due to seawater-basalt reactions.

Mineral chemistry of primary and secondary phases in basaltic rocks

Primary magmatic phases (spinel, olivine, plagioclase, clinopyroxene, amphibole, and biotite) and secondary phyllosilicates (smectite, chlorite-smectite, and celadonite) were analyzed by electron microprobe in alkalic and tholeiitic dolerites and basalts from Ocean Drilling Program Sites 800, 801, and 802. Aphyric alkalic dolerite sills (Hole 800A) and basalt flows (Holes 801B and 801C) share common mineralogical features: matrix feldspars are strongly zoned from labradorite cores to discrete sodic rims of alkali feldspar with a high Or component, which overlaps that of quench microlites in glassy mesostasis; little fractionated clinopyroxenes are Ti-rich diopsides and augites (with marked aegirine-augite rims at Site 801); rare, brown, Fe**3+-rich amphibole is winchite; and late biotites exhibit variable Ti contents. Alkalic rims to feldspars probably developed at the same time as quenched mesostasis feldspars and late-stage magmatic biotite, and represent the buildup of K-rich hydrous fluids during crystallization. Phenocryst phases in primitive mid-ocean ridge tholeiites from Hole 801C (Mg numbers about 70) have extreme compositions with chrome spinel (Cr/Cr + Al ratios about 0.2-0.4), Ni-rich olivine (Fo90), and highly calcic plagioclase (An90). Later glomerophyric clumps of plagioclase (An75-80) and clinopyroxene (diopside-augite) are strongly zoned and probably reflect rapidly changing melt conditions during upward transport, prior to seafloor quenching. In contrast, phenocryst phases (olivine, plagioclase, and clinopyroxene) in the Hole 802A tholeiites show limited variation and do not have such primitive compositions, reflecting the uniform and different chemical composition of all the bulk rocks. Replacive phyllosilicates in both alkalic and tholeiitic basalts include various colored smectites (Fe-, Mg-, and Al-saponites), chlorite-smectite and celadonite. Smectite compositions typically reflect the replaced host composition; glass is replaced by brown Fe-saponites (variable Fe/Mg ratios) and olivine by greenish Mg-saponites (or Al-rich chlorite-smectite).

Geochemistry of primary and secondary mineral phases in ODP Leg 209 sediments

Primary and secondary mineral phases from Holes 1268A (11 samples), 1272A (9 samples), and 1274A (12 samples) were analyzed by electron microprobe in Bonn and Cologne (Germany). Bulk rock powders of these samples were also analyzed geochemically, including major and trace elements (Paulick et al., 2006, doi:10.1016/j.chemgeo.2006.04.011). Ocean Drilling Program (ODP) Leg 209 Holes 1268A, 1272A, and 1274A differ remarkably in alteration intensity and mineralogy, and details regarding their lithologic characteristics are presented in Bach et al. (2004, doi:10.1029/2004GC000744) and Shipboard Scientific Party (2004, doi:10.2973/odp.proc.ir.209.101.2004). Because of the least altered character of peridotite in Hole 1274A, abundant clinopyroxene, orthopyroxene, olivine, and spinel were analyzed at this site. In Hole 1272A, primary silicates are rare and analyses were restricted to some samples that contain traces of olivine and orthopyroxene. Because of the intensity of alteration, Hole 1268A is devoid of primary phases except spinel. Commonly, alteration is pseudomorphic and serpentinization of olivine and orthopyroxene can be distinguished. Accordingly, compositional variations of the alteration minerals with regard to the precursor minerals are one of the issues investigated in this data report.

Secondary minerals and geochemistry at DSDP Legs 68 and 69

Basalt samples recovered during DSDP Legs 68, 69, and 70 from a 550-meter-thick section in two holes near the Costa Rica Rift (Holes 501 and 504B) were found to contain the following secondary minerals: trioctahedral and dioctahedral smectite, chlorite, mixed-layer clays, talc, hematite, pyrite, foujasite, phillipsite, analcime, natrolite, thomsonite, gyrolite, aragonite, calcite, anhydrite, chalcocite, Fe-hydrosilicate, okenite, apophyllite, actinolite, cristobalite, quartz, and magnesite. A less positive identification of bismutite was made. A mineral rich in Mn and minerals with strong reflections at 12.9 Å and 3.20 Å remain unidentified. Trioctahedral smectite replaces glass and olivine in the basalt groundmass. The other secondary minerals occur in veins. The distribution of the secondary minerals in the basalt section shows both hydrothermal and oxidizing-nonoxidizing zonation. Most of the secondary minerals formed under alkaline, nonoxidizing conditions at temperatures up to 120° C. An acidic regime probably existed in the lowest portion of basalt. Oxidative diagenesis followed nonoxidative diagenesis in the upper part of the section. Oxidative diagenesis is characterized by the absence of celadonite, rare occurrences of dioctahedral smectite, and widespread hematite and phillipsite.

(Table 7) Representative electron microprobe analyses for natroalunite and larger FeOx spheroids of experiment ADSU6

Acid-sulfate alteration of basalt by SO2-bearing volcanic vapors has been proposed as one possible origin for sulfate-rich deposits on Mars. To better define mineralogical signatures of acid-sulfate alteration, laboratory experiments were performed to investigate alteration pathways and geochemical processes during reaction of basalt with sulfuric acid. Pyroclastic cinders composed of phenocrysts including plagioclase, olivine, and augite embedded in glass were reacted with sulfuric acid at 145 °C for up to 137 days at a range of fluid : rock ratios. During the experiments, the phenocrysts reacted rapidly to form secondary products, while the glass was unreactive. Major products included amorphous silica, anhydrite, and Fe-rich natroalunite, along with minor iron oxides/oxyhydroxides (probably hematite) and trace levels of other sulfates. At the lowest fluid : rock ratio, hexahydrite and an unidentified Fe-silicate phase also occurred as major products. Reaction-path models indicated that formation of the products required both slow dissolution of glass and kinetic inhibitions to precipitation of a number of minerals including phyllosilicates and other aluminosilicates as well as Al- and Fe-oxides/oxyhydroxides. Similar models performed for Martian basalt compositions predict that the initial stages of acid-sulfate alteration of pyroclastic deposits on Mars should result in formation of amorphous silica, anhydrite, Fe-bearing natroalunite, and kieserite, along with relict basaltic glass. In addition, analysis of the experimental products indicates that Fe-bearing natroalunite produces a Mössbauer spectrum closely resembling that of jarosite, suggesting that it should be considered an alternative to the component in sulfate-rich bedrocks at Meridiani Planum that has previously been identified as jarosite.