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.

Palaeoceanographic data of sediment core MD01-2416

We present an SiF4 separation line, coupled to a laser fluorination system, which allows for an efficient combined silica d18O and d30Si analysis (50 min per sample). The required sample weight of 1.5-2.0 mg allows for high-resolution isotope studies on biogenic opal. Besides analytical tests, the new instrumentation set-up was used to analyse two marine diatom fractions (>63 µm, 10-20 µm) with different diatom species compositions extracted from a Bølling/Allerød-Holocene core section [MD01-2416, North-West (NW) Pacific] to evaluate the palaeoceanographic significance of the diatom isotopic signals and to address isotopic effects related to contamination and species-related isotope effects (vital and environmental effects). While d30Si offsets between the two fractions were not discernible, supporting the absence of species-related silicon isotope effects, systematic offsets occur between the d18O records. Although small, these offsets point to species-related isotope effects, as bias by contamination can be discarded. The new records strengthen the palaeoceanographic history during the last deglaciation in the NW Pacific characterized by a sequence of events with varying surface water structure and biological productivity. With such palaeoceanographic evolution it becomes unlikely that the observed systematic d18O offsets signal seasonal temperature variability. This calls for reconsideration of vital effects, generally excluded to affect d18O measurements.

A multiproxy analysis of late MIS 6 to early MIS 5d OSP Site 161-975 sediments

A multiproxy analysis based on planktic foraminiferal abundances, derived SSTs, and stable planktic isotopes measurements together with alkenone abundances and Uk'37 SSTs was performed on late MIS 6 to early MIS 5d sediment recovered from Site 975 (ODP Leg 161) in the South Balearic Islands basin (Western Mediterranean) with emphasis on reconstructing the climate progression of the last interglacial period. A number of abrupt climate changes related to alternative influence of nutrient rich northern and oligotrophic southern water masses were revealed. Heinrich event 11 and cooling events C27, C26, C25, C24, C23, which have been previously described in the North Atlantic, were recognized. However, in comparison to the eastern North Atlantic mid-latitude region, events C27 and C26 at Site 975 seem to be significantly more pronounced. Together with evidence of a two-phase climate optimum with maximum SSTs reached during its later phase, this implies a close similarity in climate dynamics between the Western Mediterranean and the Nordic seas. We propose that postglacial effects in the Nordic seas had an influence on the western Mediterranean climate via atmospheric circulation and that these effects competed with the insolation force.

Age determination and sedimentation rates of the subpolar northeast Atlantic

In order to monitor the evolution of the British-Irish Ice Sheet (BIIS) and its influence in surface ocean structure during marine isotopic stages (MIS) 2 and 3, we have analyzed the sediments recovered in core MD04-2829CQ (Rosemary Bank, north Rockall Trough, northeast Atlantic) dated between ~41 and ~18 ka B.P. Ice-rafted debris flux and composition, 40Ar/39Ar ages of individual hornblende grains, multispecies planktonic stable isotope records, planktonic foraminifera assemblage data and faunal-based sea surface temperatures (SSTs) demonstrate a close interaction between BIIS dynamics and surface ocean structure and water properties in this region. The core location lies beneath the North Atlantic Current (NAC) and is ideal for monitoring the shifts in the position of its associated oceanic fronts, as recorded by faunal changes. These data reveal a succession of BIIS-sourced iceberg calving events related to low SST, usually synchronous with dramatic changes in the composition of the planktonic foraminifera assemblage and with variations in the stable isotope records of the taxa Neogloboquadrina pachyderma (sinistral coiling) and Globigerina bulloides. The pacing of the calving events, from typically Dansgaard-Oeschger millennial timescales during late MIS 3 to multicentennial cyclicity from ~28 ka B.P., represents the build-up of the BIIS and its growing instability toward Heinrich Event (HE) 2 and the Last Glacial Maximum. Our data confirm the strong coupling between BIIS instabilities and the temperature and salinity of surface waters in the adjacent northeast Atlantic and demonstrate the BIIS's ability to modify the NAC on its flow toward the Nordic Seas. In contrast, subsurface water masses were less affected except during the Greenland stadials that contain HEs, when most intense water column reorganizations occurred simultaneously with the deposition of cream-colored carbonate sourced from the Laurentide Ice Sheet.

Radiative properties for warm and hot dense matter

We will present calculations of opacities for matter under LTE conditions. Opacities are needed in radiation transport codes to study processes like Inertial Confinement Fusion and plasma amplifiers in X-ray secondary sources. For the calculations we use the code BiGBART, with either a hydrogenic approximation with j-splitting or self-consistent data generated with the atomic physics code FAC. We calculate the atomic structure, oscillator strengths, radiative transition energies, including UTA computations, and photoionization cross-sections. A DCA model determines the configurations considered in the computation of the opacities. The opacities obtained with these two models are compared with experimental measurements.

Icosahedral Ni nanowires formed from nanocontacts breaking: Identification and characterization by Molecular Dynamics

We present and discuss an algorithm to identify and characterize the long icosahedral structures (staggered pentagonal nanowires with 1-5-1-5 atomic structure) that appear in Molecular Dynamics simulations of metallic nanowires of different species subjected to stretching. The use of this algorithm allows the identification of pentagonal rings forming the icosahedral structure as well as the determination of its number np , and the maximum length of the pentagonal nanowire Lpm. The algorithm is tested with some ideal structures to show its ability to discriminate between pentagonal rings and other ring structures. We applied the algorithm to Ni nanowires with temperatures ranging between 4K and 865K, stretched along the [111], [100] and [110] directions. We studied statistically the formation of pentagonal nanowires obtaining the distributions of length Lpm and number of rings np as function of the temperature. The Lpm distribution presents a peaked shape, with peaks located at fixed distances whose separation corresponds to the distance between two consecutive pentagonal rings.

Independent mobility of catalytic and regulatory domains of myosin heads

The recent determination of the myosin head atomic structure has led to a new model of muscle contraction, according to which mechanical torque is generated in the catalytic domain and amplified by the lever arm made of the regulatory domain [Fisher, A. J., Smith, C. A., Thoden, J., Smith, R., Sutoh, K., Holden, H. M. & Rayment, I. (1995) Biochemistry 34, 8960–8972]. A crucial aspect of this model is the ability of the regulatory domain to move independently of the catalytic domain. Saturation transfer–EPR measurements of mobility of these two domains in myosin filaments give strong support for this notion. The catalytic domain of the myosin head was labeled at Cys-707 with indane dione spin label; the regulatory domain was labeled at the single cysteine residue of the essential light chain and exchanged into myosin. The mobility of the regulatory domain in myosin filaments was characterized by an effective rotational correlation time (τR) between 24 and 48 μs. In contrast, the mobility of the catalytic domain was found to be τR = 5–9 μs. This difference in mobility between the two domains existed only in the filament form of myosin. In the monomeric form, or when bound to actin, the mobility of the two domains in myosin was indistinguishable, with τR = 1–4 μs and >1,000 μs, respectively. Therefore, the observed difference in filaments cannot be ascribed to differences in local conformations of the spin-labeled sites. The most straightforward interpretation suggests a flexible hinge between the two domains, which would have to stiffen before force could be generated.

Architecture and mechanism of the light-harvesting apparatus of purple bacteria

Photosynthetic organisms fuel their metabolism with light energy and have developed for this purpose an efficient apparatus for harvesting sunlight. The atomic structure of the apparatus, as it evolved in purple bacteria, has been constructed through a combination of x-ray crystallography, electron microscopy, and modeling. The detailed structure and overall architecture reveals a hierarchical aggregate of pigments that utilizes, as shown through femtosecond spectroscopy and quantum physics, elegant and efficient mechanisms for primary light absorption and transfer of electronic excitation toward the photosynthetic reaction center.

Vacancy-induced magnetism in graphene and graphene ribbons

We address the electronic structure and magnetic properties of vacancies and voids both in graphene and graphene ribbons. By using a mean-field Hubbard model, we study the appearance of magnetic textures associated with removing a single atom (vacancy) and multiple adjacent atoms (voids) as well as the magnetic interactions between them. A simple set of rules, based on the Lieb theorem, link the atomic structure and the spatial arrangement of the defects to the emerging magnetic order. The total spin S of a given defect depends on its sublattice imbalance, but some defects with S=0 can still have local magnetic moments. The sublattice imbalance also determines whether the defects interact ferromagnetically or antiferromagnetically with one another and the range of these magnetic interactions is studied in some simple cases. We find that in semiconducting armchair ribbons and two-dimensional graphene without global sublattice imbalance, there is a maximum defect density above which local magnetization disappears. Interestingly, the electronic properties of semiconducting graphene ribbons with uncoupled local moments are very similar to those of diluted magnetic semiconductors, presenting giant Zeeman splitting.

Proceedings of the Cambridge Philosophical Society.

Mode of access: Internet.

Expression and biochemical analysis of the entire HIV-2 gp41 ectodomain: determinants of stability map to N- and C-terminal sequences outside the 6-helix bundle core

The folding of HIV gp41 into a 6-helix bundle drives virus-cell membrane fusion. To examine the structural relationship between the 6-helix bundle core domain and other regions of gp41, we expressed in Escherichia coli, the entire ectodomain of HIV-2(ST) gp41 as a soluble, trimeric maltose-binding protein (MBP)/gp41 chimera. Limiting proteolysis indicated that the Cys-591-Cys-597 disulfide-bonded region is outside a core domain comprising two peptides, Thr-529-Trp-589 and Val-604-Ser-666. A biochemical examination of MBP/gp41 chimeras encompassing these core peptides; indicated that the N-terminal polar segment, 521-528, and C-terminal membrane-proximal segment, 658-666, cooperate in stabilizing the ectodomain. A functional interaction between sequences outside the gp41 core may contribute energy to membrane fusion. (C) 2004 Published by Elsevier B.V. on behalf of the Federation of European Biochemical Societies.