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.

drp, A novel cell density regulated protein

Growing cells are continuously processing signals of all varieties and responding to these signals by changes in cellular gene expression. One signal that cells in close proximity relay to each other is cell-cell contact. Non-transformed cells respond to cell-cell contact by arrest of growth and entry into G$\sb0,$ a process known as contact inhibition. Transformed cells do not respond to contact inhibition and continue to grow to high cell density, forming foci when in cell culture and tumors in the living organism. The events surrounding the generation, transduction, and response to cellular contact are poorly understood. In the present study, a novel gene product, drp, is shown to be expressed at high levels in cultured cells at high cell density. This density regulated protein, drp, has an apparent molecular weight of 70 kDa. Northern analysis shows drp to be highly expressed in cardiac and skeletal muscle and least abundant in lung and kidney tissues. By homology to two independently derived sequence tagged sites (STSs) used in the human genome project, drp or a closely related sequence maps to human chromosome 12. Density-dependent increases in drp expression have been demonstrated in six different cell lines including NIH 3T3, Hela and a human teratocarcinoma cell line, PA-1. Cells exhibit increased drp expression both when they are plated at increasing concentrations per unit area, or plated at low density and allowed to grow naturally to higher cell density. Cells at high density can exhibit several phenotypes including growth arrest, accumulation of soluble factors in the media, and increased numbers of cell contacts. Growth arrest by serum starvation or TGF-$\beta$ treatment fails to produce an increase in drp expression. Similarly, treatment of low density cells with conditioned media from high density cells fails to elicit drp expression. These results argue that neither soluble factors accumulated or expressed at high density nor simple exit from the cell cycle is sufficient to produce an increase in drp expression. The expression of drp appears to be uniquely regulated by cell density alone. ^

REGULATION OF TOXIN SYNTHESIS BY CLOSTRIDIUM DIFFICILE

Clostridium difficile is the leading definable cause of nosocomial diarrhea worldwide due to its virulence, multi-drug resistance, spore-forming ability, and environmental persistence. The incidence of C. difficile infection (CDI) has been increasing exponentially in the last decade. Virulent strains of C. difficile produce either toxin A and/or toxin B, which are essential for the pathogenesis of this bacterium. Current methods for diagnosing CDI are mostly qualitative tests that detect the bacterium, the toxins, or the toxin genes. These methods do not differentiate virulent C. difficile strains that produce active toxins from non-virulent strains that do not produce toxins or produce inactive toxins. Based on the knowledge that C. difficile toxins A and B cleave a substrate that is stereochemically similar to the native substrate of the toxins, uridine diphosphoglucose, a quantitative, cost-efficient assay, the Cdifftox activity assay, was developed to measure C. difficile toxin activity. The concept behind the activity assay was modified to develop a novel, rapid, sensitive, and specific assay for C. difficile toxins in the form of a selective and differential agar plate culture medium, the Cdifftox Plate assay (CDPA). This assay combines in a single step the specific identification of C. difficile strains and the detection of active toxin(s). The CDPA was determined to be extremely accurate (99.8% effective) at detecting toxin-producing strains based on the analysis of 528 C. difficile isolates selected from 50 tissue culture cytotoxicity assay-positive clinical stool samples. This new assay advances and improves the culture methodology in that only C. difficile strains will grow on this medium and virulent strains producing active toxins can be differentiated from non-virulent strains. This new method reduces the time and effort required to isolate and confirm toxin-producing C. difficile strains and provides a clinical isolate for antibiotic susceptibility testing and strain typing. The Cdifftox activity assay was used to screen for inhibitors of toxin activity. Physiological levels of the common human conjugated bile salt, taurocholate, was found to inhibit toxin A and B in vitro activities. When co-incubated ex vivo with purified toxin B, taurocholate protected Caco-2 colonic epithelial cells from the damaging effects of the toxin. Furthermore, using a caspase-3 detection assay, taurocholate reduced the extent of toxin B-induced Caco-2 cell apoptosis. These results suggest that bile salts can be effective in protecting the gut epithelium from C. difficile toxin damage, thus, the delivery of physiologic amounts of taurocholate to the colon, where it is normally in low concentration, could be useful in CDI treatment. These findings may help to explain why bile rich small intestine is spared damage in CDI, while the bile salt poor colon is vulnerable in CDI. Toxin synthesis in C. difficile occurs during the stationary phase, but little is known about the regulation of these toxins. It was hypothesized that C. difficile toxin synthesis is regulated by a quorum sensing mechanism. Two lines of evidence supported this hypothesis. First, a small (KDa), diffusible, heat-stable toxin-inducing activity accumulates in the medium of high-density C. difficile cells. This conditioned medium when incubated with low-density log-phase cells causes them to produce toxin early (2-4 hrs instead of 12-16 hrs) and at elevated levels when compared with cells grown in fresh medium. These data suggested that C. difficile cells extracellularly release an inducing molecule during growth that is able to activate toxin synthesis prematurely and demonstrates for the first time that toxin synthesis in C. difficile is regulated by quorum signaling. Second, this toxin-inducing activity was partially purified from high-density stationary-phase culture supernatant fluid by HPLC and confirmed to induce early toxin synthesis, even in C. difficile virulent strains that over-produce the toxins. Mass spectrometry analysis of the purified toxin-inducing fraction from HPLC revealed a cyclic compound with a mass of 655.8 Da. It is anticipated that identification of this toxin-inducing compound will advance our understanding of the mechanism involved in the quorum-dependent regulation of C. difficile toxin synthesis. This finding should lead to the development of even more sensitive tests to diagnose CDI and may lead to the discovery of promising novel therapeutic targets that could be harnessed for the treatment C. difficile infections.

The LPS O -antigen is required for Myxococcus xanthus social motility and multicellular development

The gliding bacterium Myxococcus xanthus aggregates to form spore-filled fruiting bodies when starved at high density. All of the identified M. xanthus lipopolysaccharide (LPS) O-antigen biosynthesis mutants exhibit defective motility and fruiting-body development. To determine the cause of these phenotypes, the cell-surface properties of the LPS O-antigen mutants were compared to wild-type cells. The binding characteristics of wild-type and LPS O-antigen-defective strains to cationic resin indicate that the mutant cell surfaces are more electronegative. Antibiotic sensitivity and hexadecane adhesion assays indicate that the wild-type M. xanthus cell surface is hydrophobic, supporting the idea that phospholipids are present in the outer leaflet of the outer membrane. The absence of the LPS O-antigen appears to expose charges associated with phospholipids and LPS core/lipid A, resulting in a dramatic alteration of the cell-surface organization and charge. These differences may affect the interaction of the LPS O-antigen mutants with their substratum and neighboring cells, leading to defects in social and single-cell gliding motility and thus, deficiencies in fruiting body formation. ^ The LPS O-antigen biosynthetic mutations also bypass the requirement of 4521 gene expression for the cell-density signal, A signal. The 4521 gene is overexpressed in these mutants. This 4521 overexpression is dependent on the sensor kinase SasS. Co-development with wild-type cells, or the addition of crude polysaccharides or membrane vesicles restores the ability of LPS O-antigen mutants to form fruiting bodies and lowers 4521 developmental gene expression to wild-type levels. Wild-type vesicles may attach or incorporate into the outer membrane of the mutants that lack LPS O-antigen, restoring a wild-type periplasmic status and allowing for normal levels of 4521 activity and fruiting body formation. We propose that the LPS composition and the configuration of the outer membrane are important elements for the complex behavioral response of M. xanthus fruiting body development. ^