Mutations affecting beta-tubulin folding and degradation.

Revertants of a colcemid-resistant Chinese hamster ovary cell line with an altered (D45Y) beta-tubulin have allowed the identification of four cis-acting mutations (L187R, Y398C, a 12-amino acid in-frame deletion, and a C-terminal truncation) that act by destabilizing the mutant tubulin and preventing it from incorporating into microtubules. These unstable beta-tubulins fail to form heterodimers and are predominantly found in association with the chaperonin CCT, suggesting that they cannot undergo productive folding. In agreement with these in vivo observations, we show that the defective beta-tubulins do not stably interact with cofactors involved in the tubulin folding pathway and, hence, fail to exchange with beta-tubulin in purified alphabeta heterodimers. Treatment of cells with MG132 causes an accumulation of the aberrant tubulins, indicating that improperly folded beta-tubulin is degraded by the proteasome. Rapid degradation of the mutant tubulin does not elicit compensatory changes in wild-type tubulin synthesis or assembly. Instead, loss of beta-tubulin from the mutant allele causes a 30-40% decrease in cellular tubulin content with no obvious effect on cell growth or survival.

STABILITY AND FOLDING OF MUTANT FORMS OF ISO-1 CYTOCHROME-C FROM SACCHAROMYCES CEREVISIAE

Partially functional forms of iso-1-cytochrome c from Saccharomyces cerevisiae were obtained by replacements of the evolutionarily conserved proline 71 with valine, isoleucine and threonine (Ernst et.al.,1985). Pro-71 lies at the juncture of two short helical regions and is believed to be important for proper local polypeptide chain folding within the iso-1-cytochrome c structure.^ To study folding in the absence of intermolecular disulfide dimer formation the free sulfhydryl group of Cys-102 was modified in both wild type and mutant proteins with an alkylating reagent, methyl methanethiosulfonate. Spectral analysis of the wild type and mutant proteins shows that the native-like functional (or partially functional) folded structure of cytochrome c is retained in the chemically modified derivatives. The replacement of Pro-71 with valine, isoleucine or threonine reduces the intensity of the 696 nm absorbance band which is an indicator of the Met-80 ligation to the heme. Thermal stability and guanidine hydrochloride unfolding studies of the mutant proteins shows a destabilization of the protein as a result of mutation. The degree of destabilization depends on the chemical nature of the substituent amino acid in the mutant protiens.^ Kinetics of folding/unfolding reactions of the proteins were monitored by fluorescence changes using stopped flow mixing to obtain guanidine hydrochloride concentration jumps ending below, within, and above the transition zone. The replacement of Pro-71 alters the rate on one of the fastest phases, $\tau\sb3$, while the two other phases, $\tau\sb1$ & $\tau\sb2$, remain the same.^ Slow refolding kinetic studies indicate that replacement of Pro-71 does not completely eliminate the absorbance or fluorescence detected slow phases leading to the conclusion that Pro-71 is not involved in the generation of the slow phases in the folding kinetics of iso-1-cytochrome c.^ The alkaline conformational change involving the disappearance of the 696 nm absorbance band occurs with increasing pH in the alkaline pH region (Davis et al., 1974). The apparent pK of this conformational change in mutant proteins is shifted as much as two pH units compared to wild type. The equilibrium and kinetic data of alkaline transition for the wild type follows a simple mechanism proposed by Davis et al., (1974) for horse heart cytochrome c. A more complex mechanism is proposed for the behavior of the mutant proteins. ^

MECHANISM-BASED STRATEGIES TO ENHANCE THE ACTIONS OF A

Heat shock protein 90 (HSP90) is an abundant molecular chaperone that regulates the functional stability of client oncoproteins, such as STAT3, Raf-1 and Akt, which play a role in the survival of malignant cells. The chaperone function of HSP90 is driven by the binding and hydrolysis of ATP. The geldanamycin analog, 17-AAG, binds to the ATP pocket of HSP90 leading to the degradation of client proteins. However, treatment with 17-AAG results in the elevation of the levels of antiapoptotic proteins HSP70 and HSP27, which may lead to cell death resistance. The increase in HSP70 and HSP27 protein levels is due to the activation of the transcription factor HSF-1 binding to the promoter region of HSP70 and HSP27 genes. HSF-1 binding subsequently promotes HSP70 and HSP27 gene expression. Based on this, I hypothesized that inhibition of transcription/translation of HSP or client proteins would enhance 17-AAG-mediated cytotoxicity. Multiple myeloma (MM) cell lines MM.1S, RPMI-8226, and U266 were used as a model. To test this hypothesis, two different strategies were used. For the first approach, a transcription inhibitor was combined with 17-AAG. The established transcription inhibitor Actinomycin D (Act D), used in the clinic, intercalates into DNA and blocks RNA elongation. Stress inducible (HSP90á, HSP70 and HSP27) and constitutive (HSP90â and HSC70) mRNA and protein levels were measured using real time RT-PCR and immunoblot assays. Treatment with 0.5 µM 17-AAG for 8 hours resulted in the induction of all HSP transcript and protein levels in the MM cell lines. This induction of HSP mRNA levels was diminished by 0.05 µg/mL Act D for 12 hours in the combination treatment, except for HSP70. At the protein level, Act D abrogated the 17-AAG-mediated induction of all HSP expression levels, including HSP70. Cytotoxic evaluation (Annexin V/7-AAD assay) of Act D in combination with 17-AAG suggested additive or more than additive interactions. For the second strategy, an agent that affected bioenergy production in addition to targeting transcription and translation was used. Since ATP is necessary for the proper folding and maturation of client proteins by HSP90, ATP depletion should lead to a decrease in client protein levels. The transcription and translation inhibitor 8-Chloro-Adenosine (8-Cl-Ado), currently in clinical trials, is metabolized into its cytotoxic form 8-Cl-ATP causing a parallel decrease of the cellular ATP pool. Treatment with 0.5 µM 17-AAG for 8 hours resulted in the induction of all HSP transcript and protein levels in the three MM cell lines evaluated. In the combination treatment, 10 µM 8-Cl-Ado for 20 hours did not abrogate the induction of HSP mRNA or protein levels. Since cellular bioenergy is necessary for the stabilization of oncoproteins by HSP90, immunoblot assays analyzing for expression levels of client proteins such as STAT3, Raf-1, and Akt were performed. Immunoblot assays detecting for the phosphorylation status of the translation repressor 4E-BP1, whose activity is modulated by upstream kinases sensitive to changes in ATP levels, were also performed. The hypophosphorylated state of 4E-BP1 leads to translation repression. Data indicated that treatment with 17-AAG alone resulted in a minor (<10%) change in STAT3, Raf-1, and Akt protein levels, while no change was observed for 4E-BP1. The combination treatment resulted in more than 50% decrease of the client protein levels and hypophosphorylation of 4E-BP1 in all MM cell lines. Treatment with 8-Cl-Ado alone resulted in less than 30% decrease in client protein levels as well as a decrease in 4E-BP1 phosphorylation. Cytotoxic evaluation of 8-Cl-Ado in combination with 17-AAG resulted in more than additive cytotoxicity when drugs were combined in a sequential manner. In summary, these data suggest that the mechanism-based combination of agents that target transcription, translation, or decrease cellular bioenergy with 17-AAG results in increase cytotoxicity when compared to the single agents. Such combination strategies may be applied in the clinic since these drugs are established chemotherapeutic agents or currently in clinical trials.

Identification and phenotypic characterization of a second collagen adhesin, Scm, and genome-based identification and analysis of 13 other predicted MSCRAMMs, including four distinct pilus loci, in Enterococcus faecium.

Attention has recently been drawn to Enterococcus faecium because of an increasing number of nosocomial infections caused by this species and its resistance to multiple antibacterial agents. However, relatively little is known about the pathogenic determinants of this organism. We have previously identified a cell-wall-anchored collagen adhesin, Acm, produced by some isolates of E. faecium, and a secreted antigen, SagA, exhibiting broad-spectrum binding to extracellular matrix proteins. Here, we analysed the draft genome of strain TX0016 for potential microbial surface components recognizing adhesive matrix molecules (MSCRAMMs). Genome-based bioinformatics identified 22 predicted cell-wall-anchored E. faecium surface proteins (Fms), of which 15 (including Acm) had characteristics typical of MSCRAMMs, including predicted folding into a modular architecture with multiple immunoglobulin-like domains. Functional characterization of one [Fms10; redesignated second collagen adhesin of E. faecium (Scm)] revealed that recombinant Scm(65) (A- and B-domains) and Scm(36) (A-domain) bound to collagen type V efficiently in a concentration-dependent manner, bound considerably less to collagen type I and fibrinogen, and differed from Acm in their binding specificities to collagen types IV and V. Results from far-UV circular dichroism measurements of recombinant Scm(36) and of Acm(37) indicated that these proteins were rich in beta-sheets, supporting our folding predictions. Whole-cell ELISA and FACS analyses unambiguously demonstrated surface expression of Scm in most E. faecium isolates. Strikingly, 11 of the 15 predicted MSCRAMMs clustered in four loci, each with a class C sortase gene; nine of these showed similarity to Enterococcus faecalis Ebp pilus subunits and also contained motifs essential for pilus assembly. Antibodies against one of the predicted major pilus proteins, Fms9 (redesignated EbpC(fm)), detected a 'ladder' pattern of high-molecular-mass protein bands in a Western blot analysis of cell surface extracts from E. faecium, suggesting that EbpC(fm) is polymerized into a pilus structure. Further analysis of the transcripts of the corresponding gene cluster indicated that fms1 (ebpA(fm)), fms5 (ebpB(fm)) and ebpC(fm) are co-transcribed, a result consistent with those for pilus-encoding gene clusters of other Gram-positive bacteria. All 15 genes occurred frequently in 30 clinically derived diverse E. faecium isolates tested. The common occurrence of MSCRAMM- and pilus-encoding genes and the presence of a second collagen-binding protein may have important implications for our understanding of this emerging pathogen.

Denaturation and unfolding of human anaphylatoxin C3a: an unusually low covalent stability of its native disulfide bonds.

The complement C3a anaphylatoxin is a major molecular mediator of innate immunity. It is a potent activator of mast cells, basophils and eosinophils and causes smooth muscle contraction. Structurally, C3a is a relatively small protein (77 amino acids) comprising a N-terminal domain connected by 3 native disulfide bonds and a helical C-terminal segment. The structural stability of C3a has been investigated here using three different methods: Disulfide scrambling; Differential CD spectroscopy; and Reductive unfolding. Two uncommon features regarding the stability of C3a and the structure of denatured C3a have been observed in this study. (a) There is an unusual disconnection between the conformational stability of C3a and the covalent stability of its three native disulfide bonds that is not seen with other disulfide proteins. As measured by both methods of disulfide scrambling and differential CD spectroscopy, the native C3a exhibits a global conformational stability that is comparable to numerous proteins with similar size and disulfide content, all with mid-point denaturation of [GdmCl](1/2) at 3.4-5M. These proteins include hirudin, tick anticoagulant protein and leech carboxypeptidase inhibitor. However, the native disulfide bonds of C3a is 150-1000 fold less stable than those proteins as evaluated by the method of reductive unfolding. The 3 native disulfide bonds of C3a can be collectively and quantitatively reduced with as low as 1mM of dithiothreitol within 5 min. The fragility of the native disulfide bonds of C3a has not yet been observed with other native disulfide proteins. (b) Using the method of disulfide scrambling, denatured C3a was shown to consist of diverse isomers adopting varied extent of unfolding. Among them, the most extensively unfolded isomer of denatured C3a is found to assume beads-form disulfide pattern, comprising Cys(36)-Cys(49) and two disulfide bonds formed by two pair of consecutive cysteines, Cys(22)-Cys(23) and Cys(56)-Cys(57), a unique disulfide structure of polypeptide that has not been documented previously.

Activation of DegP chaperone-protease via formation of large cage-like oligomers upon binding to substrate proteins.

Cells use molecular chaperones and proteases to implement the essential quality control mechanism of proteins. The DegP (HtrA) protein, essential for the survival of Escherichia coli cells at elevated temperatures with homologues found in almost all organisms uniquely has both functions. Here we report a mechanism for DegP to activate both functions via formation of large cage-like 12- and 24-mers after binding to substrate proteins. Cryo-electron microscopic and biochemical studies revealed that both oligomers are consistently assembled by blocks of DegP trimers, via pairwise PDZ1-PDZ2 interactions between neighboring trimers. Such interactions simultaneously eliminate the inhibitory effects of the PDZ2 domain. Additionally, both DegP oligomers were also observed in extracts of E. coli cells, strongly implicating their physiological importance.

Structure of the Hsp110:Hsc70 nucleotide exchange machine.

Hsp70s mediate protein folding, translocation, and macromolecular complex remodeling reactions. Their activities are regulated by proteins that exchange ADP for ATP from the nucleotide-binding domain (NBD) of the Hsp70. These nucleotide exchange factors (NEFs) include the Hsp110s, which are themselves members of the Hsp70 family. We report the structure of an Hsp110:Hsc70 nucleotide exchange complex. The complex is characterized by extensive protein:protein interactions and symmetric bridging interactions between the nucleotides bound in each partner protein's NBD. An electropositive pore allows nucleotides to enter and exit the complex. The role of nucleotides in complex formation and dissociation, and the effects of the protein:protein interactions on nucleotide exchange, can be understood in terms of the coupled effects of the nucleotides and protein:protein interactions on the open-closed isomerization of the NBDs. The symmetrical interactions in the complex may model other Hsp70 family heterodimers in which two Hsp70s reciprocally act as NEFs.

All in the family: atypical Hsp70 chaperones are conserved modulators of Hsp70 activity.

Divergent relatives of the Hsp70 protein chaperone such as the Hsp110 and Grp170 families have been recognized for some time, yet their biochemical roles remained elusive. Recent work has revealed that these "atypical" Hsp70s exist in stable complexes with classic Hsp70s where they exert a powerful nucleotide-exchange activity that synergizes with Hsp40/DnaJ-type cochaperones to dramatically accelerate Hsp70 nucleotide cycling. This represents a novel evolutionary transition from an independent protein-folding chaperone to what appears to be a dedicated cochaperone. Contributions of the atypical Hsp70s to established cellular roles for Hsp70 now must be deciphered.

The role of Sse1 in the de novo formation and variant determination of the [PSI+] prion.

Yeast prions are a group of non-Mendelian genetic elements transmitted as altered and self-propagating conformations. Extensive studies in the last decade have provided valuable information on the mechanisms responsible for yeast prion propagation. How yeast prions are formed de novo and what cellular factors are required for determining prion "strains" or variants--a single polypeptide capable of existing in multiple conformations to result in distinct heritable phenotypes--continue to defy our understanding. We report here that Sse1, the yeast ortholog of the mammalian heat-shock protein 110 (Hsp110) and a nucleotide exchange factor for Hsp70 proteins, plays an important role in regulating [PSI+] de novo formation and variant determination. Overproduction of the Sse1 chaperone dramatically enhanced [PSI+] formation whereas deletion of SSE1 severely inhibited it. Only an unstable weak [PSI+] variant was formed in SSE1 disrupted cells whereas [PSI+] variants ranging from very strong to very weak were formed in isogenic wild-type cells under identical conditions. Thus, Sse1 is essential for the generation of multiple [PSI+] variants. Mutational analysis further demonstrated that the physical association of Sse1 with Hsp70 but not the ATP hydrolysis activity of Sse1 is required for the formation of multiple [PSI+] variants. Our findings establish a novel role for Sse1 in [PSI+] de novo formation and variant determination, implying that the mammalian Hsp110 may likewise be involved in the etiology of protein-folding diseases.

Molecular cross talk between misfolded proteins in animal models of Alzheimer's and prion diseases.

The central event in protein misfolding disorders (PMDs) is the accumulation of a misfolded form of a naturally expressed protein. Despite the diversity of clinical symptoms associated with different PMDs, many similarities in their mechanism suggest that distinct pathologies may cross talk at the molecular level. The main goal of this study was to analyze the interaction of the protein misfolding processes implicated in Alzheimer's and prion diseases. For this purpose, we inoculated prions in an Alzheimer's transgenic mouse model that develop typical amyloid plaques and followed the progression of pathological changes over time. Our findings show a dramatic acceleration and exacerbation of both pathologies. The onset of prion disease symptoms in transgenic mice appeared significantly faster with a concomitant increase on the level of misfolded prion protein in the brain. A striking increase in amyloid plaque deposition was observed in prion-infected mice compared with their noninoculated counterparts. Histological and biochemical studies showed the association of the two misfolded proteins in the brain and in vitro experiments showed that protein misfolding can be enhanced by a cross-seeding mechanism. These results suggest a profound interaction between Alzheimer's and prion pathologies, indicating that one protein misfolding process may be an important risk factor for the development of a second one. Our findings may have important implications to understand the origin and progression of PMDs.

Plasticity of lipid-protein interactions in the function and topogenesis of the membrane protein lactose permease from Escherichia coli.

Phosphatidylcholine (PC) has been widely used in place of naturally occurring phosphatidylethanolamine (PE) in reconstitution of bacterial membrane proteins. However, PC does not support native structure or function for several reconstituted transport proteins. Lactose permease (LacY) of Escherichia coli, when reconstituted in E. coli phospholipids, exhibits energy-dependent uphill and energy-independent downhill transport function and proper conformation of periplasmic domain P7, which is tightly linked to uphill transport function. LacY expressed in cells lacking PE and containing only anionic phospholipids exhibits only downhill transport and lacks native P7 conformation. Reconstitution of LacY in the presence of E. coli-derived PE, but not dioleoyl-PC, results in uphill transport. We now show that LacY exhibits uphill transport and native conformation of P7 when expressed in a mutant of E. coli in which PC completely replaces PE even though the structure is not completely native. E. coli-derived PC and synthetic PC species containing at least one saturated fatty acid also support the native conformation of P7 dependent on the presence of anionic phospholipids. Our results demonstrate that the different effects of PE and PC species on LacY structure and function cannot be explained by differences in the direct interaction of the lipid head groups with specific amino acid residues alone but are due to more complex effects of the physical and chemical properties of the lipid environment on protein structure. This conclusion is supported by the effect of different lipids on the proper folding of domain P7, which indirectly influences uphill transport function.

Cross currents in protein misfolding disorders: interactions and therapy.

Protein Misfolding Disorders (PMDs) are a group of diseases characterized by the accumulation of abnormally folded proteins. Despite the wide range of proteins and tissues involved, PMDs share similar molecular and pathogenic mechanisms. Several epidemiological, clinical and experimental reports have described the co-existence of PMDs, suggesting a possible cross-talk between them. A better knowledge of the molecular basis of PMDs could have important implications for understanding the mechanism by which these diseases appear and progress and ultimately to develop novel strategies for treatment. Due to their similar molecular mechanisms, common therapeutic strategies could be applied for the diseases in this group.

Symmetric inverse consistent nonlinear registration driven by mutual information.

A nonlinear viscoelastic image registration algorithm based on the demons paradigm and incorporating inverse consistent constraint (ICC) is implemented. An inverse consistent and symmetric cost function using mutual information (MI) as a similarity measure is employed. The cost function also includes regularization of transformation and inverse consistent error (ICE). The uncertainties in balancing various terms in the cost function are avoided by alternatively minimizing the similarity measure, the regularization of the transformation, and the ICE terms. The diffeomorphism of registration for preventing folding and/or tearing in the deformation is achieved by the composition scheme. The quality of image registration is first demonstrated by constructing brain atlas from 20 adult brains (age range 30-60). It is shown that with this registration technique: (1) the Jacobian determinant is positive for all voxels and (2) the average ICE is around 0.004 voxels with a maximum value below 0.1 voxels. Further, the deformation-based segmentation on Internet Brain Segmentation Repository, a publicly available dataset, has yielded high Dice similarity index (DSI) of 94.7% for the cerebellum and 74.7% for the hippocampus, attesting to the quality of our registration method.

BIOCHEMICAL CHARACTERIZATION OF BINDING PARTNERS OF TWO HSP70 CO-CHAPERONES IN SACCHAROMYCES CEREVISIAE

Cells are exposed to a variety of environmental and physiological changes including temperature, pH and nutrient availability. These changes cause stress to cells, which results in protein misfolding and altered cellular protein homeostasis. How proteins fold into their three-dimensional functional structure is a fundamental biological process with important relevance to human health. Misfolded and aggregated proteins are linked to multiple neurodegenerative diseases, cardiovascular disease and cystic fibrosis. To combat proteotoxic stress, cells deploy an array of molecular chaperones that assist in the repair or removal of misfolded proteins. Hsp70, an evolutionarily conserved molecular chaperone, promotes protein folding and helps maintain them in a functional state. Requisite co-chaperones, including nucleotide exchange factors (NEFs) strictly regulate and serve to recruit Hsp70 to distinct cellular processes or locations. In yeast and human cells, three structurally non-related cytosolic NEFs are present: Sse1 (Hsp110), Fes1 (HspBP1) and Snl1 (Bag-1). Snl1 is unique among the cytosolic NEFs as it is localized at the ER membrane with its Hsp70 binding (BAG) domain exposed to the cytosol. I discovered that Snl1 distinctly interacts with assembled ribosomes and several lines of evidence indicate that this interaction is both independent of and concurrent with binding to Hsp70 and is not dependent on membrane localization. The ribosome-binding site is identified as a short lysine-rich motif within the amino terminus of the Snl1 BAG domain distinct from the Hsp70 interaction region. In addition, I demonstrate ribosome association with the Snl1 homolog in the pathogenic fungus, Candida albicans and localize this putative NEF to a perinuclear/ER membrane, suggesting functional conservation in fungal BAG domain-containing proteins. As a first step in determining specific domain architecture in fungal BAG proteins, I present the preliminary steps of protein purification and analysis of the minimal Hsp70 binding region in in both S.cerevisiae and C. albicans Snl1. Contrary to previous in vitro evidence which showed the Fes1 NEF to interact with both cytosolic Hsp70s, Ssa and Ssb, Fes1 is shown to interact specifically with Ssa when expressed under normal cellular conditions in S. cerevisiae. This is the first reported evidence of Hsp70 binding selectivity for a cytosolic NEF, and suggests a possible mechanism to achieve specificity in Hsp70-dependent functions. Taken together, the work presented in this dissertation highlights the striking divergence among Hsp70 co-chaperones in selecting binding partners, which may correlate with their specific roles in the cell.

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.