Electron Cryotomography of Bacterial Cells

While much is already known about the basic metabolism of bacterial cells, many fundamental questions are still surprisingly unanswered, including for instance how they generate and maintain specific cell shapes, establish polarity, segregate their genomes, and divide. In order to understand these phenomena, imaging technologies are needed that bridge the resolution gap between fluorescence light microscopy and higher-resolution methods such as X-ray crystallography and NMR spectroscopy.Electron cryotomography (ECT) is an emerging technology that does just this, allowing the ultrastructure of cells to be visualized in a near-native state, in three dimensions (3D), with "macromolecular" resolution (~4nm).^{1, 2} In ECT, cells are imaged in a vitreous, "frozen-hydrated" state in a cryo transmission electron microscope (cryoTEM) at low temperature (< -180°C). For slender cells (up to ~500 nm in thickness³), intact cells are plunge-frozen within media across EM grids in cryogens such as ethane or ethane/propane mixtures. Thicker cells and biofilms can also be imaged in a vitreous state by first "high-pressure freezing" and then, "cryo-sectioning" them. A series of two-dimensional projection images are then collected through the sample as it is incrementally tilted along one or two axes. A three-dimensional reconstruction, or "tomogram" can then be calculated from the images. While ECT requires expensive instrumentation, in recent years, it has been used in a few labs to reveal the structures of various external appendages, the structures of different cell envelopes, the positions and structures of cytoskeletal filaments, and the locations and architectures of large macromolecular assemblies such as flagellar motors, internal compartments and chemoreceptor arrays.^{1, 2}In this video article we illustrate how to image cells with ECT, including the processes of sample preparation, data collection, tomogram reconstruction, and interpretation of the results through segmentation and in some cases correlation with light microscopy.     

Visualization and interpretation of protein networks in Mycobacterium tuberculosis based on hierarchical clustering of genome-wide functional linkage maps

Genome-wide functional linkages among proteins in cellular complexes and metabolic pathways can be inferred from high throughput experimentation, such as DNA microarrays, or from bioinformatic analyses. Here we describe a method for the visualization and interpretation of genome-wide functional linkages inferred by the Rosetta Stone, Phylogenetic Profile, Operon and Conserved Gene Neighbor computational methods. This method involves the construction of a genome-wide functional linkage map, where each significant functional linkage between a pair of proteins is displayed on a two-dimensional scatter-plot, organized according to the order of genes along the chromosome. Subsequent hierarchical clustering of the map reveals clusters of genes with similar functional linkage profiles and facilitates the inference of protein function and the discovery of functionally linked gene clusters throughout the genome. We illustrate this method by applying it to the genome of the pathogenic bacterium Mycobacterium tuberculosis, assigning cellular functions to previously uncharacterized proteins involved in cell wall biosynthesis, signal transduction, chaperone activity, energy metabolism and polysaccharide biosynthesis.     

Mitochondrial gene order is not conserved in arthropods: prostriate and metastriate tick mitochondrial genomes

The entire mitochondrial genome was sequenced in a prostriate tick, Ixodes hexagonus, and a metastriate tick, Rhipicephalus sanguineus. Both genomes encode 22 tRNAs, 13 proteins, and two ribosomal RNAs. Prostriate ticks are basal members of Ixodidae and have the same gene order as Limulus polyphemus. In contrast, in R. sanguineus, a block of genes encoding NADH dehydrogenase subunit 1 (ND1), tRNA(Leu)(UUR), tRNA(Leu)(CUN), 16S rDNA, tRNA(Val), 12S rDNA, the control region, and the tRNA(Ile) and tRNA(Gln) have translocated to a position between the tRNA(Glu) and tRNA(Phe) genes. The tRNA(Cys) gene has translocated between the control region and the tRNA(Met) gene, and the tRNA(Leu)(CUN) gene has translocated between the tRNA(Ser)(UCN) gene and the control region. Furthermore, the control region is duplicated, and both copies undergo concerted evolution. Primers that flank these rearrangements confirm that this gene order is conserved in all metastriate ticks examined. Correspondence analysis of amino acid and codon use in the two ticks and in nine other arthropod mitochondrial genomes indicate a strong bias in R. sanguineus towards amino acids encoded by AT-rich codons.      

Nucleotide polymorphism at the alcohol dehydrogenase locus of pocket gophers, genus Geomys

Using the strictly neutral model as a null hypothesis, we tested for deviations from expected levels of nucleotide polymorphism at the alcohol dehydrogenase locus (Adh-1) within and among four species of pocket gophers (Geomys bursarius major, G. knoxjonesi, G. texensis llanensis, and G. attwateri). The complete protein-encoding region was examined, and 10 unique alleles, representing both electromorphic and cryptic alleles, were used to test hypotheses (e.g., the neutral model) concerning the maintenance of genetic variation. Nineteen variable sites were identified among the 10 alleles examined, including 9 segregating sites occurring in synonymous positions and 10 that were nonsynonymous. Several statistical methods, including those that test for within-species variation as well as those that examine variation within and among species, failed to reject the null hypothesis that variation (both within and between species of Geomys) at the Adh locus is consistent with the neutral theory. However, there was significant heterogeneity in the ratio of polymorphism to divergence across the gene, with polymorphisms clustered in the first half of the coding region and fixed differences clustered in the second half of the gene. Two alternative hypotheses are discussed as possible explanations for this heterogeneity: an old balanced polymorphism in the first half of the gene or a recent selective sweep in the second half of the gene.      

The genomics of disulfide bonding and protein stabilization in thermophiles

Thermophilic organisms flourish in varied high-temperature environmental niches that are deadly to other organisms. Recently, genomic evidence has implicated a critical role for disulfide bonds in the structural stabilization of intracellular proteins from certain of these organisms, contrary to the conventional view that structural disulfide bonds are exclusively extracellular. Here both computational and structural data are presented to explore the occurrence of disulfide bonds as a protein-stabilization method across many thermophilic prokaryotes. Based on computational studies, disulfide-bond richness is found to be widespread, with thermophiles containing the highest levels. Interestingly, only a distinct subset of thermophiles exhibit this property. A computational search for proteins matching this target phylogenetic profile singles out a specific protein, known as protein disulfide oxidoreductase, as a potential key player in thermophilic intracellular disulfide-bond formation. Finally, biochemical support in the form of a new crystal structure of a thermophilic protein with three disulfide bonds is presented together with a survey of known structures from the literature. Together, the results provide insight into biochemical specialization and the diversity of methods employed by organisms to stabilize their proteins in exotic environments. The findings also motivate continued efforts to sequence genomes from divergent organisms.     

Escherichia coli Peptidoglycan Structure and Mechanics as Predicted by Atomic-Scale Simulations

Bacteria face the challenging requirement to maintain their shape and avoid rupture due to the high internal turgor pressure, but simultaneously permit the import and export of nutrients, chemical signals, and virulence factors. The bacterial cell wall, a mesh-like structure composed of cross-linked strands of peptidoglycan, fulfills both needs by being semi-rigid, yet sufficiently porous to allow diffusion through it. How the mechanical properties of the cell wall are determined by the molecular features and the spatial arrangement of the relatively thin strands in the larger cellular-scale structure is not known. To examine this issue, we have developed and simulated atomic-scale models of Escherichia coli cell walls in a disordered circumferential arrangement. The cell-wall models are found to possess an anisotropic elasticity, as known experimentally, arising from the orthogonal orientation of the glycan strands and of the peptide cross-links. Other features such as thickness, pore size, and disorder are also found to generally agree with experiments, further supporting the disordered circumferential model of peptidoglycan. The validated constructs illustrate how mesoscopic structure and behavior emerge naturally from the underlying atomic-scale properties and, furthermore, demonstrate the ability of all-atom simulations to reproduce a range of macroscopic observables for extended polymer meshes.     

Activated chemoreceptor arrays remain intact and hexagonally packed

Bacterial chemoreceptors cluster into exquisitively sensitive, tunable, highly ordered, polar arrays. While these arrays serve as paradigms of cell signalling in general, it remains unclear what conformational changes transduce signals from the periplasmic tips, where attractants and repellents bind, to the cytoplasmic signalling domains. Conflicting reports support and contest the hypothesis that activation causes large changes in the packing arrangement of the arrays, up to and including their complete disassembly. Using electron cryotomography, here we show that in Caulobacter crescentus, chemoreceptor arrays in cells grown in different media and immediately after exposure to the attractant galactose all exhibit the same 12 nm hexagonal packing arrangement, array size and other structural parameters. ΔcheB and ΔcheR mutants mimicking attractant- or repellent-bound states prior to adaptation also show the same lattice structure. We conclude that signal transduction and amplification must be accomplished through only small, nanoscale conformational changes.