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.

Crystal structure of the Anabaena sensory rhodopsin transducer.

We present crystal structures of the Anabaena sensory rhodopsin transducer (ASRT), a soluble cytoplasmic protein that interacts with the first structurally characterized eubacterial retinylidene photoreceptor Anabaena sensory rhodopsin (ASR). Four crystal structures of ASRT from three different spacegroups were obtained, in all of which ASRT is present as a planar (C4) tetramer, consistent with our characterization of ASRT as a tetramer in solution. The ASRT tetramer is tightly packed, with large interfaces where the well-structured beta-sandwich portion of the monomers provides the bulk of the tetramer-forming interactions, and forms a flat, stable surface on one side of the tetramer (the beta-face). Only one of our four different ASRT crystals reveals a C-terminal alpha-helix in the otherwise all-beta protein, together with a large loop from each monomer on the opposite face of the tetramer (the alpha-face), which is flexible and largely disordered in the other three crystal forms. Gel-filtration chromatography demonstrated that ASRT forms stable tetramers in solution and isothermal microcalorimetry showed that the ASRT tetramer binds to ASR with a stoichiometry of one ASRT tetramer per one ASR photoreceptor with a K(d) of 8 microM in the highest affinity measurements. Possible mechanisms for the interaction of this transducer tetramer with the ASR photoreceptor via its flexible alpha-face to mediate transduction of the light signal are discussed.

Crystal structure of the alpha(6)beta(6) holoenzyme of propionyl-coenzyme A carboxylase.

Propionyl-coenzyme A carboxylase (PCC), a mitochondrial biotin-dependent enzyme, is essential for the catabolism of the amino acids Thr, Val, Ile and Met, cholesterol and fatty acids with an odd number of carbon atoms. Deficiencies in PCC activity in humans are linked to the disease propionic acidaemia, an autosomal recessive disorder that can be fatal in infants. The holoenzyme of PCC is an alpha(6)beta(6) dodecamer, with a molecular mass of 750 kDa. The alpha-subunit contains the biotin carboxylase (BC) and biotin carboxyl carrier protein (BCCP) domains, whereas the beta-subunit supplies the carboxyltransferase (CT) activity. Here we report the crystal structure at 3.2-A resolution of a bacterial PCC alpha(6)beta(6) holoenzyme as well as cryo-electron microscopy (cryo-EM) reconstruction at 15-A resolution demonstrating a similar structure for human PCC. The structure defines the overall architecture of PCC and reveals unexpectedly that the alpha-subunits are arranged as monomers in the holoenzyme, decorating a central beta(6) hexamer. A hitherto unrecognized domain in the alpha-subunit, formed by residues between the BC and BCCP domains, is crucial for interactions with the beta-subunit. We have named it the BT domain. The structure reveals for the first time the relative positions of the BC and CT active sites in the holoenzyme. They are separated by approximately 55 A, indicating that the entire BCCP domain must translocate during catalysis. The BCCP domain is located in the active site of the beta-subunit in the current structure, providing insight for its involvement in the CT reaction. The structural information establishes a molecular basis for understanding the large collection of disease-causing mutations in PCC and is relevant for the holoenzymes of other biotin-dependent carboxylases, including 3-methylcrotonyl-CoA carboxylase (MCC) and eukaryotic acetyl-CoA carboxylase (ACC).

Rubisco in complex with Rubisco large subunit methyltransferase.

SET domain protein lysine methyltransferases (PKMT) are a structurally unique class of enzymes that catalyze the specific methylation of lysine residues in a number of different substrates. Especially histone-specific SET domain PKMTs have received widespread attention because of their roles in the regulation of epigenetic gene expression and the development of some cancers. Rubisco large subunit methyltransferase (RLSMT) is a chloroplast-localized SET domain PKMT responsible for the formation of trimethyl-lysine-14 in the large subunit of Rubisco, an essential photosynthetic enzyme. Here, we have used cryoelectron microscopy to produce an 11-A density map of the Rubisco-RLSMT complex. The atomic model of the complex, obtained by fitting crystal structures of Rubisco and RLSMT into the density map, shows that the extensive contact regions between the 2 proteins are mainly mediated by hydrophobic residues and leucine-rich repeats. It further provides insights into potential conformational changes that may occur during substrate binding and catalysis. This study presents the first structural analysis of a SET domain PKMT in complex with its intact polypeptide substrate.

Crystal structures of progressive Ca2+ binding states of the Ca2+ sensor Ca2+ binding domain 1 (CBD1) from the CALX Na+/Ca2+ exchanger reveal incremental conformational transitions.

Na(+)/Ca(2+) exchangers (NCX) constitute a major Ca(2+) export system that facilitates the re-establishment of cytosolic Ca(2+) levels in many tissues. Ca(2+) interactions at its Ca(2+) binding domains (CBD1 and CBD2) are essential for the allosteric regulation of Na(+)/Ca(2+) exchange activity. The structure of the Ca(2+)-bound form of CBD1, the primary Ca(2+) sensor from canine NCX1, but not the Ca(2+)-free form, has been reported, although the molecular mechanism of Ca(2+) regulation remains unclear. Here, we report crystal structures for three distinct Ca(2+) binding states of CBD1 from CALX, a Na(+)/Ca(2+) exchanger found in Drosophila sensory neurons. The fully Ca(2+)-bound CALX-CBD1 structure shows that four Ca(2+) atoms bind at identical Ca(2+) binding sites as those found in NCX1 and that the partial Ca(2+) occupancy and apoform structures exhibit progressive conformational transitions, indicating incremental regulation of CALX exchange by successive Ca(2+) binding at CBD1. The structures also predict that the primary Ca(2+) pair plays the main role in triggering functional conformational changes. Confirming this prediction, mutagenesis of Glu(455), which coordinates the primary Ca(2+) pair, produces dramatic reductions of the regulatory Ca(2+) affinity for exchange current, whereas mutagenesis of Glu(520), which coordinates the secondary Ca(2+) pair, has much smaller effects. Furthermore, our structures indicate that Ca(2+) binding only enhances the stability of the Ca(2+) binding site of CBD1 near the hinge region while the overall structure of CBD1 remains largely unaffected, implying that the Ca(2+) regulatory function of CBD1, and possibly that for the entire NCX family, is mediated through domain interactions between CBD1 and the adjacent CBD2 at this hinge.

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.