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.

Neocortical hyperexcitability defect in a mutant mouse model of spike-wave epilepsy, {\it stargazer}

Single-locus mutations in mice can express epileptic phenotypes and provide critical insights into the naturally occurring defects that alter excitability and mediate synchronization in the central nervous system (CNS). One such recessive mutation (on chromosome (Chr) 15), stargazer(stg/stg) expresses frequent bilateral 6-7 cycles per second (c/sec) spike-wave seizures associated with behavioral arrest, and provides a valuable opportunity to examine the inherited lesion associated with spike-wave synchronization.^ The existence of distinct and heterogeneous defects mediating spike-wave discharge (SWD) generation has been demonstrated by the presence of multiple genetic loci expressing generalized spike-wave activity and the differential effects of pharmacological agents on SWDs in different spike-wave epilepsy models. Attempts at understanding the different basic mechanisms underlying spike-wave synchronization have focused on $\gamma$-aminobutyric acid (GABA) receptor-, low threshold T-type Ca$\sp{2+}$ channel-, and N-methyl-D-aspartate receptor (NMDA-R)-mediated transmission. It is believed that defects in these modes of transmission can mediate the conversion of normal oscillations in a trisynaptic circuit, which includes the neocortex, reticular nucleus and thalamus, into spike-wave activity. However, the underlying lesions involved in spike-wave synchronization have not been clearly identified.^ The purpose of this research project was to locate and characterize a distinct neuronal hyperexcitability defect favoring spike-wave synchronization in the stargazer brain. One experimental approach for anatomically locating areas of synchronization and hyperexcitability involved an attempt to map patterns of hypersynchronous activity with antibodies to activity-induced proteins.^ A second approach to characterizing the neuronal defect involved examining the neuronal responses in the mutant following application of pharmacological agents with well known sites of action.^ In order to test the hypothesis that an NMDA receptor mediated hyperexcitability defect exists in stargazer neocortex, extracellular field recordings were used to examine the effects of CPP and MK-801 on coronal neocortical brain slices of stargazer and wild type perfused with 0 Mg$\sp{2+}$ artificial cerebral spinal fluid (aCSF).^ To study how NMDA receptor antagonists might promote increased excitability in stargazer neocortex, two basic hypotheses were tested: (1) NMDA receptor antagonists directly activate deep layer principal pyramidal cells in the neocortex of stargazer, presumably by opening NMDA receptor channels altered by the stg mutation; and (2) NMDA receptor antagonists disinhibit the neocortical network by blocking recurrent excitatory synaptic inputs onto inhibitory interneurons in the deep layers of stargazer neocortex.^ In order to test whether CPP might disinhibit the 0 Mg$\sp{2+}$ bursting network in the mutant by acting on inhibitory interneurons, the inhibitory inputs were pharmacologically removed by application of GABA receptor antagonists to the cortical network, and the effects of CPP under 0 Mg$\sp{2+}$aCSF perfusion in layer V of stg/stg were then compared with those found in +/+ neocortex using in vitro extracellular field recordings. (Abstract shortened by UMI.) ^

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. ^

Mean dynamic topography and oceanographic parameters estimated from an inverse model and satellite geodesy, with link to model result in one single NetCDF file (392 MB), including the inverse data error covariance

Geostrophic surface velocities can be derived from the gradients of the mean dynamic topography-the difference between the mean sea surface and the geoid. Therefore, independently observed mean dynamic topography data are valuable input parameters and constraints for ocean circulation models. For a successful fit to observational dynamic topography data, not only the mean dynamic topography on the particular ocean model grid is required, but also information about its inverse covariance matrix. The calculation of the mean dynamic topography from satellite-based gravity field models and altimetric sea surface height measurements, however, is not straightforward. For this purpose, we previously developed an integrated approach to combining these two different observation groups in a consistent way without using the common filter approaches (Becker et al. in J Geodyn 59(60):99-110, 2012, doi:10.1016/j.jog.2011.07.0069; Becker in Konsistente Kombination von Schwerefeld, Altimetrie und hydrographischen Daten zur Modellierung der dynamischen Ozeantopographie, 2012, http://nbn-resolving.de/nbn:de:hbz:5n-29199). Within this combination method, the full spectral range of the observations is considered. Further, it allows the direct determination of the normal equations (i.e., the inverse of the error covariance matrix) of the mean dynamic topography on arbitrary grids, which is one of the requirements for ocean data assimilation. In this paper, we report progress through selection and improved processing of altimetric data sets. We focus on the preprocessing steps of along-track altimetry data from Jason-1 and Envisat to obtain a mean sea surface profile. During this procedure, a rigorous variance propagation is accomplished, so that, for the first time, the full covariance matrix of the mean sea surface is available. The combination of the mean profile and a combined GRACE/GOCE gravity field model yields a mean dynamic topography model for the North Atlantic Ocean that is characterized by a defined set of assumptions. We show that including the geodetically derived mean dynamic topography with the full error structure in a 3D stationary inverse ocean model improves modeled oceanographic features over previous estimates.

Fission track analyses of ODP Hole 129-800A samples

Four models of fission track annealing in apatite are compared with measured fission track lengths in samples from Site 800 in the East Mariana Basin, Ocean Drilling Program Leg 129, given an independently determined temperature history. The temperature history of Site 800 was calculated using a one-dimensional, compactive, conductive heat flow model assuming two end-member thermal cases: one for cooling of Jurassic ocean crust that has experienced no subsequent heating, and one for cooling of Cretaceous ocean crust. Because the samples analyzed were only shallowly buried and because the tectonic history of the area since sample deposition is simple, resolution of the temperature history is high. The maximum temperature experienced by the sampled bed is between 16°-21°C and occurs at 96 Ma; temperatures since the Cretaceous have dropped in spite of continued pelagic sediment deposition because heat flow has continued to decay exponentially and bottom-water temperatures have dropped. Fission tracks observed within apatite grains from the sampled bed are 14.6 +/- 0.1 µm (1 sigma) long. Given the proposed temperature history of the samples, one unpublished and three published models of fission track annealing predict mean track lengths from 14.8 to 15.9 µm. These models require temperatures as much as 40°C higher than the calculated paleotemperature maximum of the sampled bed to produce the same degree of track annealing. Measured and predicted values are different because annealing models are based on extrapolation of high temperature laboratory data to geologic times. The model that makes the closest prediction is based on the greatest number of experiments performed at low temperature and on an apatite having composition closest to that of the core samples.

Larcform 1 - single column model intercomparison of Arctic air formation, link to model results in NetCDF format

Weather and climate models struggle to represent lower tropospheric temperature and moisture profiles and surface fluxes in Arctic winter, not least because they lack or misrepresent physical processes that are specific to high latitudes. The Arctic boundary layer in winter has been observed to be in either a radiatively clear or cloudy state: The radiatively clear state is characterized by strong surface radiative cooling leading to the build-up of surface-based temperature inversions, whereas the cloudy state occurs when cloud liquid water is present in the atmospheric column, allowing little or no surface radiative cooling and leading to weaker and typically elevated temperature inversions. Many large-scale models have been shown to lack the cloudy state, and some do substantially underestimate stability in the clear state. We here present results from the first Lagrangian ARCtic air FORMation experiment (Larcform 1), a GASS (Global atmospheric system studies) single-column model intercomparison which reproduces these biases of large-scale models in an idealised setup.

Development of an adequate model for verification of design safety-margins of the HTTR nuclear test facility

This work is based on the prototype High Engineering Test Reactor (HTTR) of the Japan Agency of Energy Atomic (JAEA). Its objective is to describe an adequate deterministic model to be used in the assessment of its design safety margins via damage domains. The concept of damage domain is defined and it is shown its relevance in the ongoing effort to apply dynamic risk assessment methods and tools based on the Theory of Stimulated Dynamics (TSD). To illustrate, we present results of an abnormal control rod (CR) withdrawal during subcritical condition and its comparison with results obtained by JAEA. No attempt is made yet to actually assess the detailed scenarios, rather to show how the approach may handle events of its kind

Self-consistent, nonlocal electron heat flux at arbitrary ion charge number

A single, nonlocal expression for the electron heat flux, which closely reproduces known results at high and low ion charge number 2, and “exact” results for the local limit at all 2, is derived by solving the kinetic equation in a narrow, tail-energy range. The solution involves asymptotic expansions of Bessel functions of large argument, and (Z-dependent)order above or below it, corresponding to the possible parabolic or hyperbolic character of the kinetic equation; velocity space diffusion in self-scattering is treated similarly to isotropic thermalization of tail energies in large Z analyses. The scale length H characterizing nonlocal effects varies with Z, suggesting an equal dependence of any ad hoc flux limiter. The model is valid for all H above the mean-free path for thermal electrons.

Influencia del viento en los ensayos de rebase en modelo físico a escala reducida = Wind effect on wave overtopping in physical model test

El artículo que se presenta a continuación recoge la ampliación de una investigación previa sobre los rebases en los espaldones de los diques verticales y en talud. Para ello se han realizado una serie de ensayos en modelo físico a escala reducida sobre la sección vertical del Dique de Levante de Málaga, cuyo objeto principal fue analizar el efecto del viento en el rebase. Los ensayos se han realizado generando oleaje con y sin viento, comparando los resultados obtenidos en cada una de las dos situaciones y se han llevado a cabo en el Canal de Oleaje y Viento de Gran Escala existente en el Laboratorio de Experimentación Marítima del Centro de Estudios de Puertos y Costas del CEDEX. The purpose of the research work as summarised in this article, resulting from diverse work carried out at the CEDEX, is to make an analysis of the influence of wind effects on the wave overtopping of vertical sea-walls. The results obtained in the Málaga´s Levante breakwater tests are presented here. The test was carried out in large sized facilities where waves and wind are generated.