Crystal structure of a protein associated with cell division from Mycoplasma pneumoniae (GI: 13508053): a novel fold with a conserved sequence motif

MraZ Transcriptionally Controls the Critical Level of FtsL Required for Focusing Z-Rings and Kickstarting Septation in Bacillus subtilis

The bacterial division and cell wall (dcw) cluster is a highly conserved region of the genome which encodes several essential cell division factors, including the central divisome protein FtsZ. Understanding the regulation of this region is key to our overall understanding of the division process. mraZ is found at the 5' end of the dcw cluster, and previous studies have described MraZ as a sequence-specific DNA binding protein. In this article, we investigate MraZ to elucidate its role in Bacillus subtilis. Through our investigation, we demonstrate that increased levels of MraZ result in lethal filamentation due to repression of its own operon (mraZ-mraW-ftsL-pbpB). We observed rescue of filamentation upon decoupling ftsL expression, but not other genes in the operon, from MraZ control. Our data suggest that regulation of the mra operon may be an alternative way for cells to quickly arrest cytokinesis, potentially during entry into the stationary phase and in the event of DNA replication arrest. Furthermore, through time-lapse microscopy, we were able to identify that overexpression of mraZ or depletion of FtsL results in decondensation of the FtsZ ring (Z-ring). Using fluorescent d-amino acid labeling, we also observed that coordinated peptidoglycan insertion at the division site is dysregulated in the absence of FtsL. Thus, we reveal that the precise role of FtsL is in Z-ring maturation and focusing septal peptidoglycan synthesis. IMPORTANCE MraZ is a highly conserved protein found in a diverse range of bacteria, including genome-reduced Mycoplasma. We investigated the role of MraZ in Bacillus subtilis and found that overproduction of MraZ is toxic due to cell division inhibition. Upon further analysis, we observed that MraZ is a repressor of its own operon, which includes genes that encode the essential cell division factors FtsL and PBP2B. We noted that decoupling of ftsL alone was sufficient to abolish MraZ-mediated cell division inhibition. Using time-lapse microscopy, we showed that under conditions where the FtsL level is depleted, the cell division machinery is unable to initiate cytokinesis. Thus, our results pinpoint that the precise role of FtsL is in concentrating septal cell wall synthesis to facilitate cell division.

Binding site of MraZ transcription factor in Mollicutes

Mollicutes (mycoplasmas) feature a significant loss of known regulators of gene expression. Here, we identified the recognition site of the MraZ-family regulator of Mycoplasma gallisepticum, which is conserved in many species of different clades within class Mollicutes. The MraZ binding site is AAAGTG[T/G], in the promoter of mraZ gene it forms a series of direct repeats with a structure (AAAGTG[T/G]N3)k, where k = 3 most frequently. MraZ binds to a single repeat as an octamer complex. MraZ can also bind a single binding site or a series of repeats with different spacer lengths (2-4 nt); thus, it may play a role in the regulation of multiple operons in Mollicutes. In M. gallisepticum, MraZ acts as a transcriptional activator. The overexpression of MraZ leads to moderate filamentation of cells and the formation of aggregates, likely as a result of incomplete cytokinesis.

The highly conserved MraZ protein is a transcriptional regulator in Escherichia coli

The mraZ and mraW genes are highly conserved in bacteria, both in sequence and in their position at the head of the division and cell wall (dcw) gene cluster. Located directly upstream of the mraZ gene, the Pmra promoter drives the transcription of mraZ and mraW, as well as many essential cell division and cell wall genes, but no regulator of Pmra has been found to date. Although MraZ has structural similarity to the AbrB transition state regulator and the MazE antitoxin and MraW is known to methylate the 16S rRNA, mraZ and mraW null mutants have no detectable phenotypes. Here we show that overproduction of Escherichia coli MraZ inhibited cell division and was lethal in rich medium at high induction levels and in minimal medium at low induction levels. Co-overproduction of MraW suppressed MraZ toxicity, and loss of MraW enhanced MraZ toxicity, suggesting that MraZ and MraW have antagonistic functions. MraZ-green fluorescent protein localized to the nucleoid, suggesting that it binds DNA. Consistent with this idea, purified MraZ directly bound a region of DNA containing three direct repeats between Pmra and the mraZ gene. Excess MraZ reduced the expression of an mraZ-lacZ reporter, suggesting that MraZ acts as a repressor of Pmra, whereas a DNA-binding mutant form of MraZ failed to repress expression. Transcriptome sequencing (RNA-seq) analysis suggested that MraZ also regulates the expression of genes outside the dcw cluster. In support of this, purified MraZ could directly bind to a putative operator site upstream of mioC, one of the repressed genes identified by RNA-seq.

Quantification of mRNA and protein and integration with protein turnover in a bacterium

Determination of the average cellular copy number of 400 proteins under different growth conditions and integration with protein turnover and absolute mRNA levels reveals the dynamics of protein expression in the genome-reduced bacterium Mycoplasma pneumoniae.

Our study provides a fine-grained, quantitative picture to unprecedented detail in an established model organism for systems-wide studies.

Our integrative approach reveals a novel, dynamic view on the processes, interactions and regulations underlying the central dogma pathway and the composition of protein complexes.

Simulations using our quantitative data on mRNA, protein and turnover show how an organism copes with stochastic noise in gene expression in vivo.

Our data serve as an important resource for colleagues both within our field of research and in related disciplines.

A hallmark of Systems Biology is the integration of diverse, large quantitative data sets with the aim to gain novel insights into how biological processes work. We measured individual mRNA and protein abundances as well as protein turnover in the bacterium Mycoplasma pneumoniae. This human pathogen is an ideal model organism for organism-wide studies. It can be readily cultured under laboratory conditions and it has a very small genome with only 690 protein-coding genes. This comparably low complexity allows for the exhaustive analysis of major cellular biomolecules avoiding constrains introduced by limitations of available analysis techniques.

Using a recently developed mass spectrometry-based approach, we determined the average cellular copy number for over 400 individual proteins under different growth and stress conditions. The 20 most abundant proteins, including Elongation factor Tu, cellular chaperones, and proteins involved in metabolizing glucose, the major energy source of M. pneumoniae account for nearly 44% of the total cellular protein mass. We observed abundance changes of many expected and several unexpected proteins in response to cellular stress, such as heat shock, DNA damage and osmotic stress, as well as along batch culture growth over 4 days.

Integration of the protein abundance data with quantitative mRNA measurements revealed a modest correlation between these two classes of biomolecules. However, for several classical stress-induced proteins, we observed a correlated induction of mRNA and protein in response to heat shock. A focused analysis of mRNA–protein abundance dynamics during batch culture growth suggested that the regulation of gene expression is largely decoupled from protein dynamics in M. pneumoniae, indicating extensive post-transcriptional and post-translational regulation influencing the cellular mRNA–protein ratios.

To investigate the factors influencing the cellular protein abundance, we measured individual protein turnover rates by mass spectrometry using a label-chase approach involving stable isotope-labelled amino acids. The average half-life of a protein in M. pneumoniae is 23 h. Based on the measured quantitative mRNA data, the protein abundances and their half-lives, we established an ordinary differential equations model for the estimation of individual in vivo protein degradation and translation efficiency rates. We found out that translation efficiency rather than protein turnover is the dominating factor influencing protein abundance. Using our abundance and turnover data, we additionally performed stochastic simulations of gene expression. We observed that long protein half-life and low translational efficiency buffers gene expression noise propagating from low cellular mRNA levels in vivo.

We compared the abundance ratios of proteins associating into complexes in vivo with their expected functional stoichiometries. We observed that for stable protein complexes, such as the GroEL/ES chaperonin or DNA gyrase, our measured abundance ratios reflected the expected subunit stoichiometries. More dynamic protein complexes, such as the DnaK/J/GrpE chaperone system or RNA polymerase, showed several unusual subunit ratios, pointing towards transient interaction of sub-stoichiometric subunits for function. A detailed, quantitative analysis of the ribosome, the largest cellular protein complex, revealed large abundance differences of the 51 subunits. This observation indicates a multi-functionality for several, abundant ribosomal proteins.

Finally, a comparison of the determined average cellular protein abundances with a different pathogenic bacterium, Leptospira interrogans, revealed that cellular protein abundances closely reflect their respective lifestyles.

Our study represents an organism-wide, quantitative analysis of cellular protein abundances. Integrating our proteomics data with determined mRNA levels and protein turnover rates reveals insights into the dynamic interplay and regulation of mRNA and proteins, the central biomolecules of a cell.

Biological function and cellular responses to environmental perturbations are regulated by a complex interplay of DNA, RNA, proteins and metabolites inside cells. To understand these central processes in living systems at the molecular level, we integrated experimentally determined abundance data for mRNA, proteins, as well as individual protein half-lives from the genome-reduced bacterium Mycoplasma pneumoniae. We provide a fine-grained, quantitative analysis of basic intracellular processes under various external conditions. Proteome composition changes in response to cellular perturbations reveal specific stress response strategies. The regulation of gene expression is largely decoupled from protein dynamics and translation efficiency has a higher regulatory impact on protein abundance than protein turnover. Stochastic simulations using in vivo data show how low translation efficiency and long protein half-lives effectively reduce biological noise in gene expression. Protein abundances are regulated in functional units, such as complexes or pathways, and reflect cellular lifestyles. Our study provides a detailed integrative analysis of average cellular protein abundances and the dynamic interplay of mRNA and proteins, the central biomolecules of a cell.

Rapid model building of alpha-helices in electron-density maps

A method for the identification of alpha-helices in electron-density maps at low resolution followed by interpretation at moderate to high resolution is presented. Rapid identification is achieved at low resolution, where alpha-helices appear as tubes of density. The positioning and direction of the alpha-helices is obtained at moderate to high resolution, where the positions of side chains can be seen. The method was tested on a set of 42 experimental electron-density maps at resolutions ranging from 1.5 to 3.8 A. An average of 63% of the alpha-helical residues in these proteins were built and an average of 76% of the residues built matched helical residues in the refined models of the proteins. The overall average r.m.s.d. between main-chain atoms in the modeled alpha-helices and the nearest atom with the same name in the refined models of the proteins was 1.3 A.

Rapid model building of beta-sheets in electron-density maps

A method for rapidly building beta-sheets into electron-density maps is presented. beta-Strands are identified as tubes of high density adjacent to and nearly parallel to other tubes of density. The alignment and direction of each strand are identified from the pattern of high density corresponding to carbonyl and C(beta) atoms along the strand averaged over all repeats present in the strand. The beta-strands obtained are then assembled into a single atomic model of the beta-sheet regions. The method was tested on a set of 42 experimental electron-density maps at resolutions ranging from 1.5 to 3.8 A. The beta-sheet regions were nearly completely built in all but two cases, the exceptions being one structure at 2.5 A resolution in which a third of the residues in beta-sheets were built and a structure at 3.8 A in which under 10% were built. The overall average r.m.s.d. of main-chain atoms in the residues built using this method compared with refined models of the structures was 1.5 A.

Rapid chain tracing of polypeptide backbones in electron-density maps

A method for rapid chain tracing of polypeptide backbones at moderate resolution is presented.

A method for the rapid tracing of polypeptide backbones has been developed. The method creates an approximate chain tracing that is useful for visual evaluation of whether a structure has been solved and for use in scoring the quality of electron-density maps. The essence of the method is to (i) sample candidate C^α positions at spacings of approximately 0.6 Å along ridgelines of high electron density, (ii) list all possible nonapeptides that satisfy simple geometric and density criteria using these candidate C^α positions, (iii) score the nonapeptides and choose the highest scoring ones, and (iv) find the longest chains that can be made by connecting nonamers. An indexing and storage scheme that allows a single calculation of most distances and density values is used to speed up the process. The method was applied to 42 density-modified electron-density maps at resolutions from 1.5 to 3.8 Å. A total of 21 428 residues in these maps were traced in 24 CPU min with an overall r.m.s.d. of 1.61 Å for C^α atoms compared with the known refined structures. The method appears to be suitable for rapid evaluation of electron-density map quality.

Decision-making in structure solution using Bayesian estimates of map quality: the PHENIX AutoSol wizard

Ten measures of experimental electron-density-map quality are examined and the skewness of electron density is found to be the best indicator of actual map quality. A Bayesian approach to estimating map quality is developed and used in the PHENIX AutoSol wizard to make decisions during automated structure solution.

Estimates of the quality of experimental maps are important in many stages of structure determination of macromolecules. Map quality is defined here as the correlation between a map and the corresponding map obtained using phases from the final refined model. Here, ten different measures of experimental map quality were examined using a set of 1359 maps calculated by re-analysis of 246 solved MAD, SAD and MIR data sets. A simple Bayesian approach to estimation of map quality from one or more measures is presented. It was found that a Bayesian estimator based on the skewness of the density values in an electron-density map is the most accurate of the ten individual Bayesian estimators of map quality examined, with a correlation between estimated and actual map quality of 0.90. A combination of the skewness of electron density with the local correlation of r.m.s. density gives a further improvement in estimating map quality, with an overall correlation coefficient of 0.92. The PHENIX AutoSol wizard carries out automated structure solution based on any combination of SAD, MAD, SIR or MIR data sets. The wizard is based on tools from the PHENIX package and uses the Bayesian estimates of map quality described here to choose the highest quality solutions after experimental phasing.

Common evolutionary origin of swapped-hairpin and double-psi beta barrels

The core of swapped-hairpin and double-psi beta barrels is formed by duplication of a conserved betaalphabeta element, suggesting a common evolutionary origin. The path connecting the two folds is unclear as the two barrels are not interconvertible by a simple topological modification, such as circular permutation. We have identified a protein family whose sequence properties are intermediate to the two folds. The structure of one of these proteins, Pyrococcus horikoshii PhS018, is also built by duplication of the conserved betaalphabeta element but shows yet a third topology, which we name the RIFT barrel. This topology is widespread in the structure database and spans three folds of the SCOP classification, including the middle domain of EF-Tu and the N domain of F1-ATPase. We propose that swapped-hairpin beta barrels arose from an ancestral RIFT barrel by strand invasion and double-psi beta barrels by a strand swap. We group the three barrel types into a metafold, the cradle-loop barrels.

Mycoplasmas: a distinct cytoskeleton for wall-less bacteria

The bacterial genus Mycoplasma includes a large number of highly genomically-reduced species which in nature are associated with hosts either commensally or pathogenically. Several Mycoplasma species, including Mycoplasma pneumoniae, feature a multifunctional polar structure, the terminal organelle. Essential for colonization of the host and for gliding motility, the terminal organelle is associated with an internal cytoskeleton crucial to its assembly and function. This cytoskeleton is structurally and compositionally novel as compared with the cytoskeletons of other organisms, including other bacteria, is also involved in the cell division process. In this review we discuss the cytoskeletal structures and protein components of the attachment organelle and how they might interact and contribute to its various functions.

Revised structure of the AbrB N-terminal domain unifies a diverse superfamily of putative DNA-binding proteins

New relationships found in the process of updating the structural classification of proteins (SCOP) database resulted in the revision of the structure of the N-terminal, DNA-binding domain of the transition state regulator AbrB. The dimeric AbrB domain shares a common fold with the addiction antidote MazE and the subunit of uncharacterized protein MraZ implicated in cell division and cell envelope formation. It has a detectable sequence similarity to both MazE and MraZ thus providing an evolutionary link between the two proteins. The putative DNA-binding site of AbrB is found on the same face as the DNA-binding site of MazE and appears similar, both in structure and sequence, to the exposed conserved region of MraZ. This strongly suggests that MraZ also binds DNA and allows for a consensus model of DNA recognition by the members of this novel protein superfamily.

Transcriptional analysis of the conserved ftsZ gene cluster in Mycoplasma genitalium and Mycoplasma pneumoniae

Several experimental approaches were used to construct a detailed transcriptional profile of the phylogenetically conserved ftsZ cell division gene cluster in both Mycoplasma genitalium and its closest relative, Mycoplasma pneumoniae. We determined initiation and termination points for the cluster, as well as an absolute steady-state RNA level for each gene. Transcription of this cluster in both these organisms was shown to be highly strand specific. While the four genes in this cluster are cotranscribed, their transcription unit also includes two genes of close proximity yet disparate function. A transcription initiation point immediately upstream of these two genes was detected in M. genitalium but not M. pneumoniae. In M. pneumoniae, transcription of the six genes terminates at a poly(U)-tailed hairpin. In M. genitalium, this transcription terminates at two closely spaced points by an unknown mechanism. Real-time reverse transcription-PCR analysis of this cluster in M. pneumoniae shows that mRNA levels for all six genes vary at most fivefold and form a gradient of decreasing quantity with increasing distance from the promoter at the beginning of the cluster. mRNA from coding regions was approximately 20- to 100-fold more abundant than that from intergenic regions. We estimated the most abundant mRNA we detected at 0.6 copy per cell. We conclude that groups of functionally related genes in M. genitalium and M. pneumoniae are often preceded by promoters but rarely followed by terminators. This causes functionally unrelated genes to be commonly cotranscribed in these organisms.

AbrB-like Transcription Factors Assume a Swapped Hairpin Fold that Is Evolutionarily Related to Double-Psi β Barrels

AbrB is a key transition-state regulator of Bacillus subtilis. Based on the conservation of a betaalphabeta structural unit, we proposed a beta barrel fold for its DNA binding domain, similar to, but topologically distinct from, double-psi beta barrels. However, the NMR structure revealed a novel fold, the "looped-hinge helix." To understand this discrepancy, we undertook a bioinformatics study of AbrB and its homologs; these form a large superfamily, which includes SpoVT, PrlF, MraZ, addiction module antidotes (PemI, MazE), plasmid maintenance proteins (VagC, VapB), and archaeal PhoU homologs. MazE and MraZ form swapped-hairpin beta barrels. We therefore reexamined the fold of AbrB by NMR spectroscopy and found that it also forms a swapped-hairpin barrel. The conservation of the core betaalphabeta element supports a common evolutionary origin for swapped-hairpin and double-psi barrels, which we group into a higher-order class, the cradle-loop barrels, based on the peculiar shape of their ligand binding site.

MraZ from Escherichia coli: cloning, purification, crystallization and preliminary X-ray analysis

The MraZ family of proteins, also referred to as the UPF0040 family, are highly conserved in bacteria and are thought to play a role in cell-wall biosynthesis and cell division. The murein region A (mra) gene cluster encodes MraZ proteins along with a number of other proteins involved in this complex process. To date, there has been no clear functional assignment provided for MraZ proteins and the structure of a homologue from Mycoplasma pneumoniae, MPN314, failed to suggest a molecular function. The b0081 gene from Escherichia coli that encodes the MraZ protein was cloned and the protein was overexpressed, purified and crystallized. This data is presented along with evidence that the E. coli homologue exists in a different oligomeric state to the MPN314 protein.