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.

High resolution study of sediment core GeoB3606-1

Based on a high-resolution analysis of the diatom signal and biogenic bulk components at site GeoB3606-1 (25°S, off Namibia), we describe rapid palaeoceanographic changes in the Benguela Upwelling System (BUS) from early MIS 3 through to the early Holocene (55 000 to 7 000 14C yr BP). Coastal upwelling strongly varied at 25°S from MIS 3 through to MIS 2. The abrupt decrease in the accumulation rate of biogenic silica and diatoms from MIS 3 into MIS 2 records rapid oceanographic changes in the BUS off Namibia. During MIS 3, leakage of excess H4SiO4 acid from the Southern Ocean into low-latitude surface waters, as indicated by the occurrence of Antarctic diatoms, enhanced the production of spores of Chaetoceros at the expense of calcareous phytoplankton. Furthermore, shallower Antarctic Intermediate Water (AAIW) would have enriched the thermocline off Namibia with silicate transported from the Southern Ocean. The strong decrease of the siliceous signal throughout MIS 2 represents a decrease in the nutrient input to the BUS, even though the diatom assemblage is still dominated by spores of the upwelling-associated diatom genus Chaetoceros. Depletion of silicate in the thermocline from the onset of MIS 2 through to the early Holocene reflects the shutdown of AAIW injection from the Southern Ocean into the BUS, causing upwelled waters to become reduced in silicate, hence less favourable for diatom production. The deglaciation and early Holocene are characterised by the replacement of the upwelling-associated flora by a non-upwelling-related diatom community, reflecting weakened upwelling, retraction of the seaward extension of the chlorophyll filament off Lüderitz, and dominance of warmer waters.

High-resolution isotope stratigraphy of ODP Site 184-1144

Narrow-spaced oxygen and carbon stable isotope records of the planktonic foraminifer Globigerinoides ruber (white) were obtained at Ocean Drilling Program Leg 184 Site 1144 to establish a first record of high-resolution Pleistocene monsoon variability on orbital to centennial timescales in the northern South China Sea. The new records extend from the Holocene back to marine isotope Stage (MIS) 34 (1.1 Ma). Sedimentation rates average 0.56 m/k.y. for the upper Matuyama and Brunhes Chrons and increase to 1.8 m/k.y. over the last 100 k.y. Stable isotope records thus reach an average time resolution of 270-500 yr for the last 375 k.y. and 570 yr further back to 700 ka. On the other hand, major stratigraphic gaps were identified for peak warm Stages 5.5, 7.5 (down to 8.4), 11.3, and 15.5. These gaps probably resulted from short-lasting events of contour current erosion induced by short-term enhanced incursions of Upper Pacific Deep Water near the end of glacial terminations. A further major hiatus extends from MIS 34 to MIS 73(?). The long-term variations in monsoon climate were largely dominated by the 100-k.y. eccentricity cycle. Planktonic delta13C values culminated near 30, 480, and 1035 ka and reflect an overlying 450-k.y. eccentricity cycle of minimum nutrient concentrations in the surface ocean. Superimposed on the orbital variations, millennial-scale cycles were prominent throughout the last 700 k.y., mainly controlled by short-term changes in monsoon-driven precipitation and freshwater input from mainland China. During the last 110 k.y. these short-lasting oscillations closely match the record of 1500-yr Dansgaard-Oeschger climate cycles in the Greenland ice core record.

High-resolution IP25 of three sediment cores in the western North Pacific and Bering Sea

Due to its strong influence on heat and moisture exchange between the ocean and the atmosphere, sea ice is an essential component of the global climate system. In the context of its alarming decrease in terms of concentration, thickness and duration, understanding the processes controlling sea-ice variability and reconstructing paleo-sea-ice extent in polar regions have become of great interest for the scientific community. In this study, for the first time, IP25, a recently developed biomarker sea-ice proxy, was used for a high-resolution reconstruction of the sea-ice extent and its variability in the western North Pacific and western Bering Sea during the past 18,000 years. To identify mechanisms controlling the sea-ice variability, IP25 data were associated with published sea-surface temperature as well as diatom and biogenic opal data. The results indicate that a seasonal sea-ice cover existed during cold periods (Heinrich Stadial 1 and Younger Dryas), whereas during warmer intervals (Bølling-Allerød and Holocene) reduced sea ice or ice-free conditions prevailed in the study area. The variability in sea-ice extent seems to be linked to climate anomalies and sea-level changes controlling the oceanographic circulation between the subarctic Pacific and the Bering Sea, especially the Alaskan Stream injection though the Aleutian passes.

High resolution isotope and Mg/Ca record of sediment core GeoB12615-4

The importance of intermediate water masses in climate change and ocean circulation has been emphasized recently. In particular, Southern Ocean Intermediate Waters (SOIW), such as Antarctic Intermediate Water and Subantarctic Mode Water, are thought to have acted as active interhemispheric transmitter of climate anomalies. Here we reconstruct changes in SOIW signature and spatial and temporal evolution based on a 40 kyr time series of oxygen and carbon isotopes as well as planktic Mg/Ca based thermometry from Site GeoB12615-4 in the western Indian Ocean. Our data suggest that SOIW transmitted Antarctic temperature trends to the equatorial Indian Ocean via the "oceanic tunnel" mechanism. Moreover, our results reveal that deglacial SOIW carried a signature of aged Southern Ocean deep water. We find no evidence of increased formation of intermediate waters during the deglaciation.

Quantifying ice-shelf basal melt and ice-stream dynamics using high-resolution DEM and GPS time series

Thesis (Ph.D.)--University of Washington, 2016-06

High resolution provenancing of long travelled dust deposited on the Southern Alps, New Zealand

On 7 February 2000 an atypical orange discolouration of snowfields in the central Southern Alps, New Zealand occurred following the passage of a cold front. Analysis of snow samples identified fine orangey-brown dust mixed with much coarser grey dust. Air parcel forward trajectories from dust sources in southern and central Australia, where dust storms were reported on 4 February 2000, were computed to identify the deposits source. Geochemical analyses of the dust deposit using 26 trace elements, unaffected by regional air pollution and gravitational sorting, indicate that 20% of the dust was sourced from western New South Wales, with 45% from the eastern Eyre Peninsula of South Australia and the remaining 35% was local New Zealand dust. This provenancing approach provides a spatial resolution of long travelled dust sourcing not previously achieved. (c) 2005 Elsevier B.V. All rights reserved.

A high-resolution sedimentology and age determination of core PS2185 from the central Arctic Ocean

A high-resolution multiparameter stratigraphy allows the identification of late Quaternary glacial and interglacial cycles in a central Arctic Ocean sediment core. Distinct sandy layers in the upper part of the otherwise fine-grained sediment core from the Lomonosov Ridge (lat 87.5°N) correlate to four major glacials since ca. 0.7 Ma. The composition of these ice-rafted terrigenous sediments points to a glaciated northern Siberia as the main source. In contrast, lithic carbonates derived from North America are also present in older sediments and indicate a northern North American glaciation since at least 2.8 Ma. We conclude that large-scale northern Siberian glaciation began much later than other Northern Hemisphere ice sheets.

High resolution study of narrow tailed radio galaxies in clusters

Clusters of galaxies are the most massive and large gravitationally bounded systems in the whole Universe. Their study is of fundamental importance to constrain cosmological parameters and to obtain informations regarding various kind of emission in different wavebands. In particular, in the radio domain, beside the diffuse emission, the study is focused on the radio galaxies emission. Radio galaxies in clusters can have peculiar morphology, since they interact with the intracluster medium (ICM) in which they are embedded. Particularly, in this thesis we focused our attention on the so-called Narrow-Angle Tailed radio galaxies (NAT), which present radio jets that are bent at extreme angle, up to 90 degrees, from their original orientation. Some NAT show a narrow extended structure and the two radio tails are not resolved even with high resolution radio observations. An example is provided by the source IC310, in the Perseus Cluster, whose structure has been recently interpreted as due to Doppler boosting effects of a relativistic jet oriented at a small angle with respect to the line of sight. If the structure is due to relativistic effects, this implies that the jets are relativistic at about 400 kpc from the core, but this is in contrast with unified models, which predict that for low-power radio source (NAT are classified as FRI radio galaxies) the jets decelerate to sub-relativistic speed within a few kpc from the core. To investigate this scientific topic, in this thesis we have analyzed the innermost structure of a sample of eleven radio galaxies showing a very narrow NAT structure. We can conclude that the structure of these radio galaxies is different from that of IC310. These radio galaxies are indeed strongly influenced by environmental effects and are similar to classical NAT sources.