A peptide profile of the Bacillus subtilis genome

Three-helix bundle and SH3-type barrels: autonomously stable structural motifs in small and large proteins

In this study, we investigated two variants of a three-helix bundle and SH3-type barrel, compact in space, present in small and large proteins of various living organisms. Using a neural graph network, proteins with three-helix bundle (n = 1377) and SH3-type barrels (n = 1914) spatial folds were selected. Molecular experiments were performed for small proteins with these folds, and motifs were studied autonomously outside the protein environment at 300, 340, and 370 K. A comparative analysis of the main parameters of the structures in the course of the experiment was performed, including gyration radius, area accessible to the solvent, number of hydrophobic and hydrogen bonds, and root-mean-square deviation of atomic positions (RMSD). We exhibited an autonomous stability of the studied folds outside the protein environment in an aquatic medium. We aimed to demonstrate the possibility of analyzing three-helix bundle and SH3-type barrels autonomously outside the protein globule, thereby reducing the computational time and increasing performance without significant loss of information.Communicated by Ramaswamy H. Sarma.

Prophage-encoded small protein YqaH counteracts the activities of the replication initiator DnaA in Bacillus subtilis

Bacterial genomes harbour cryptic prophages that are mostly transcriptionally silent with many unannotated genes. Still, cryptic prophages may contribute to their host fitness and phenotypes. In Bacillus subtilis, the yqaF-yqaN operon belongs to the prophage element skin, and is tightly repressed by the Xre-like repressor SknR. This operon contains several small ORFs (smORFs) potentially encoding small-sized proteins. The smORF-encoded peptide YqaH was previously reported to bind to the replication initiator DnaA. Here, using a yeast two-hybrid assay, we found that YqaH binds to the DNA binding domain IV of DnaA and interacts with Spo0A, a master regulator of sporulation. We isolated single amino acid substitutions in YqaH that abolished the interaction with DnaA but not with Spo0A. Then, using a plasmid-based inducible system to overexpress yqaH WT and mutant derivatives, we studied in B. subtilis the phenotypes associated with the specific loss-of-interaction with DnaA (DnaA_LOI). We found that expression of yqaH carrying DnaA_LOI mutations abolished the deleterious effects of yqaH WT expression on chromosome segregation, replication initiation and DnaA-regulated transcription. When YqaH was induced after vegetative growth, DnaA_LOI mutations abolished the drastic effects of YqaH WT on sporulation and biofilm formation. Thus, YqaH inhibits replication, sporulation and biofilm formation mainly by antagonizing DnaA in a manner that is independent of the cell cycle checkpoint Sda.

Sterilization effects of UV laser irradiation on Bacillus atrophaeus spore viability, structure, and proteins

Bacillus spores are highly resistant to toxic chemicals and extreme environments. Because some Bacillus species threaten public health, spore inactivation techniques have been intensively investigated. We exposed Bacillus atrophaeus spores to a 266 nm Nd:YVO₄ laser at a laser power of 1 W and various numbers of scans. As a result, the UV laser reduced the viability of Bacillus atrophaeus spores. Although the outer coat of spores remained intact after UV laser irradiation of 720 scans, damage inside the spores was observed. Spore proteins were identified by matrix-assisted laser desorption/ionization Fourier transform ion cyclotron resonance mass spectrometry during the course of UV laser irradiation. Photochemical and photothermal processes are believed to be involved in the UV laser sterilization of Bacillus spores. Our findings suggest that a UV laser is capable of sterilizing Bacillus atrophaeus spores.

Small but Smart: On the Diverse Role of Small Proteins in the Regulation of Cyanobacterial Metabolism

Over the past few decades, bioengineered cyanobacteria have become a major focus of research for the production of energy carriers and high value chemical compounds. Besides improvements in cultivation routines and reactor technology, the integral understanding of the regulation of metabolic fluxes is the key to designing production strains that are able to compete with established industrial processes. In cyanobacteria, many enzymes and metabolic pathways are regulated differently compared to other bacteria. For instance, while glutamine synthetase in proteobacteria is mainly regulated by covalent enzyme modifications, the same enzyme in cyanobacteria is controlled by the interaction with unique small proteins. Other prominent examples, such as the small protein CP12 which controls the Calvin-Benson cycle, indicate that the regulation of enzymes and/or pathways via the attachment of small proteins might be a widespread mechanism in cyanobacteria. Accordingly, this review highlights the diverse role of small proteins in the control of cyanobacterial metabolism, focusing on well-studied examples as well as those most recently described. Moreover, it will discuss their potential to implement metabolic engineering strategies in order to make cyanobacteria more definable for biotechnological applications.

Structure-seq2 probing of RNA structure upon amino acid starvation reveals both known and novel RNA switches in Bacillus subtilis

RNA structure influences numerous processes in all organisms. In bacteria, these processes include transcription termination and attenuation, small RNA and protein binding, translation initiation, and mRNA stability, and can be regulated via metabolite availability and other stresses. Here we use Structure-seq2 to probe the in vivo RNA structurome of Bacillus subtilis grown in the presence and absence of amino acids. Our results reveal that amino acid starvation results in lower overall dimethyl sulfate (DMS) reactivity of the transcriptome, indicating enhanced protection owing to protein binding or RNA structure. Starvation-induced changes in DMS reactivity correlated inversely with transcript abundance changes. This correlation was particularly pronounced in genes associated with the stringent response and CodY regulons, which are involved in adaptation to nutritional stress, suggesting that RNA structure contributes to transcript abundance change in regulons involved in amino acid metabolism. Structure-seq2 accurately reported on four known amino acid-responsive riboswitches: T-box, SAM, glycine, and lysine riboswitches. Additionally, we discovered a transcription attenuation mechanism that reduces yfmG expression when amino acids are added to the growth medium. We also found that translation of a leader peptide (YfmH) encoded just upstream of yfmG regulates yfmG expression. Our results are consistent with a model in which a slow rate of yfmH translation caused by limitation of the amino acids encoded in YfmH prevents transcription termination in the yfmG leader region by favoring formation of an overlapping antiterminator structure. This novel RNA switch offers a way to simultaneously monitor the levels of multiple amino acids.

An integrative strategy to identify the entire protein coding potential of prokaryotic genomes by proteogenomics

Accurate annotation of all protein-coding sequences (CDSs) is an essential prerequisite to fully exploit the rapidly growing repertoire of completely sequenced prokaryotic genomes. However, large discrepancies among the number of CDSs annotated by different resources, missed functional short open reading frames (sORFs), and overprediction of spurious ORFs represent serious limitations. Our strategy toward accurate and complete genome annotation consolidates CDSs from multiple reference annotation resources, ab initio gene prediction algorithms and in silico ORFs (a modified six-frame translation considering alternative start codons) in an integrated proteogenomics database (iPtgxDB) that covers the entire protein-coding potential of a prokaryotic genome. By extending the PeptideClassifier concept of unambiguous peptides for prokaryotes, close to 95% of the identifiable peptides imply one distinct protein, largely simplifying downstream analysis. Searching a comprehensive Bartonella henselae proteomics data set against such an iPtgxDB allowed us to unambiguously identify novel ORFs uniquely predicted by each resource, including lipoproteins, differentially expressed and membrane-localized proteins, novel start sites and wrongly annotated pseudogenes. Most novelties were confirmed by targeted, parallel reaction monitoring mass spectrometry, including unique ORFs and single amino acid variations (SAAVs) identified in a re-sequenced laboratory strain that are not present in its reference genome. We demonstrate the general applicability of our strategy for genomes with varying GC content and distinct taxonomic origin. We release iPtgxDBs for B. henselae, Bradyrhizobium diazoefficiens and Escherichia coli and the software to generate both proteogenomics search databases and integrated annotation files that can be viewed in a genome browser for any prokaryote.

Common and phylogenetically widespread coding for peptides by bacterial small RNAs

Background

While eukaryotic noncoding RNAs have recently received intense scrutiny, it is becoming clear that bacterial transcription is at least as pervasive. Bacterial small RNAs and antisense RNAs (sRNAs) are often assumed to be noncoding, due to their lack of long open reading frames (ORFs). However, there are numerous examples of sRNAs encoding for small proteins, whether or not they also have a regulatory role at the RNA level.

Methods

Here, we apply flexible machine learning techniques based on sequence features and comparative genomics to quantify the prevalence of sRNA ORFs under natural selection to maintain protein-coding function in 14 phylogenetically diverse bacteria. Importantly, we quantify uncertainty in our predictions, and follow up on them using mass spectrometry proteomics and comparison to datasets including ribosome profiling.

Results

A majority of annotated sRNAs have at least one ORF between 10 and 50 amino acids long, and we conservatively predict that 409±191.7 unannotated sRNA ORFs are under selection to maintain coding (mean estimate and 95% confidence interval), an average of 29 per species considered here. This implies that overall at least 10.3±0.5% of sRNAs have a coding ORF, and in some species around 20% do. 165±69 of these novel coding ORFs have some antisense overlap to annotated ORFs. As experimental validation, many of our predictions are translated in published ribosome profiling data and are identified via mass spectrometry shotgun proteomics. B. subtilis sRNAs with coding ORFs are enriched for high expression in biofilms and confluent growth, and S. pneumoniae sRNAs with coding ORFs are involved in virulence. sRNA coding ORFs are enriched for transmembrane domains and many are predicted novel components of type I toxin/antitoxin systems.

Conclusions

We predict over two dozen new protein-coding genes per bacterial species, but crucially also quantified the uncertainty in this estimate. Our predictions for sRNA coding ORFs, along with predicted novel type I toxins and tools for sorting and visualizing genomic context, are freely available in a user-friendly format at http://disco-bac.web.pasteur.fr. We expect these easily-accessible predictions to be a valuable tool for the study not only of bacterial sRNAs and type I toxin-antitoxin systems, but also of bacterial genetics and genomics.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3932-y) contains supplementary material, which is available to authorized users.

Small proteins: untapped area of potential biological importance

Polypeptides containing ≤100 amino acid residues (AAs) are generally considered to be small proteins (SPs). Many studies have shown that some SPs are involved in important biological processes, including cell signaling, metabolism, and growth. SP generally has a simple domain and has an advantage to be used as model system to overcome folding speed limits in protein folding simulation and drug design. But SPs were once thought to be trivial molecules in biological processes compared to large proteins. Because of the constraints of experimental methods and bioinformatics analysis, many genome projects have used a length threshold of 100 amino acid residues to minimize erroneous predictions and SPs are relatively under-represented in earlier studies. The general protein discovery methods have potential problems to predict and validate SPs, and very few effective tools and algorithms were developed specially for SPs identification. In this review, we mainly consider the diverse strategies applied to SPs prediction and discuss the challenge for differentiate SP coding genes from artifacts. We also summarize current large-scale discovery of SPs in species at the genome level. In addition, we present an overview of SPs with regard to biological significance, structural application, and evolution characterization in an effort to gain insight into the significance of SPs.

Contributions of the σ(W) , σ(M) and σ(X) regulons to the lantibiotic resistome of Bacillus subtilis

In Bacillus subtilis, the extracytoplasmic function (ECF) σ factors σ(M) , σ(W) and σ(X) all contribute to resistance against lantibiotics. Nisin, a model lantibiotic, has a dual mode of action: it inhibits cell wall synthesis by binding lipid II, and this complex also forms pores in the cytoplasmic membrane. These activities can be separated in a nisin hinge-region variant (N20P M21P) that binds lipid II, but no longer permeabilizes membranes. The major contribution of σ(M) to nisin resistance is expression of ltaSa, encoding a stress-activated lipoteichoic acid synthase, and σ(X) functions primarily by activation of the dlt operon controlling d-alanylation of teichoic acids. Together, σ(M) and σ(X) regulate cell envelope structure to decrease access of nisin to its lipid II target. In contrast, σ(W) is principally involved in protection against membrane permeabilization as it provides little protection against the nisin hinge region variant. σ(W) contributes to nisin resistance by regulation of a signal peptide peptidase (SppA), phage shock proteins (PspA and YvlC, a PspC homologue) and tellurite resistance related proteins (YceGHI). These defensive mechanisms are also effective against other lantibiotics such as mersacidin, gallidermin and subtilin and comprise an important subset of the intrinsic antibiotic resistome of B. subtilis.

Single transmembrane peptide DinQ modulates membrane-dependent activities

The functions of several SOS regulated genes in Escherichia coli are still unknown, including dinQ. In this work we characterize dinQ and two small RNAs, agrA and agrB, with antisense complementarity to dinQ. Northern analysis revealed five dinQ transcripts, but only one transcript (+44) is actively translated. The +44 dinQ transcript translates into a toxic single transmembrane peptide localized in the inner membrane. AgrB regulates dinQ RNA by RNA interference to counteract DinQ toxicity. Thus the dinQ-agr locus shows the classical features of a type I TA system and has many similarities to the tisB-istR locus. DinQ overexpression depolarizes the cell membrane and decreases the intracellular ATP concentration, demonstrating that DinQ can modulate membrane-dependent processes. Augmented DinQ strongly inhibits marker transfer by Hfr conjugation, indicating a role in recombination. Furthermore, DinQ affects transformation of nucleoid morphology in response to UV damage. We hypothesize that DinQ is a transmembrane peptide that modulates membrane-dependent activities such as nucleoid compaction and recombination.

Exposure of the bacterium Escherichia coli to DNA damaging agents induces the SOS response, which up-regulates gene functions involved in numerous cellular processes such as DNA repair, cell division, and replication. Most of the SOS regulated genes in E. coli have been characterized, but still there are several genes of unknown function. One of these uncharacterized genes is dinQ. In this work we characterize dinQ and two novel small RNAs, agrA and agrB, that regulate expression of dinQ. The DinQ peptide is localized in the inner membrane as a single transmembrane peptide of 27 amino acids. Small proteins of less than 50 amino acids are important in cellular processes such as regulation, signalling, and antibacterial action. Here we demonstrate that DinQ modulates recombination and transformation of nucleoid morphology in response to UV damage. Our results provide new insights into small hydrophobic peptides that could regulate important DNA metabolic processes dependent on the inner membrane of the cell.

Hydrophobic peptides: novel regulators within bacterial membrane

Identification of short coding sequences is challenging, both experimentally and in silico, and functional natural peptides (< 50 amino acids) have to a large extent been overlooked in Gram-negative bacteria. Recent results have converged to highlight the role of hydrophobic peptides that form a novel class of active molecules in Escherichia coli and Salmonella enterica serovar Typhimurium. These peptides can play a regulatory role by interacting with protein partners at the inner membrane and by modulating protein partner activity or stability. Genome-wide analyses in both bacterial species have identified several conserved short open reading frames encoding a single transmembrane segment. We discuss the known and predicted membrane-associated peptides and the tools for their identification. Besides the identification of novel regulatory networks, characterization of peptides with a single transmembrane helix segment and proteins that interact with them provides a powerful opportunity to study interactions between alpha helices within biological membranes. In addition, some bioactive membrane peptides could provide a basis for engineering membrane protein antagonists.

A genome-wide survey of short coding sequences in streptococci

Identification of short genes that encode peptides of fewer than 60 aa is challenging, both experimentally and in silico. As a consequence, the universe of these short coding sequences (CDSs) remains largely unknown, although some are acknowledged to play important roles in cell-cell communication, particularly in Gram-positive bacteria. This paper reports a thorough search for short CDSs across streptococcal genomes. Our bioinformatic approach relied on a combination of advanced intrinsic and extrinsic methods. In the first step, intrinsic sequence information (nucleotide composition and presence of RBSs) served to identify new short putative CDSs (spCDSs) and to eliminate the differences between annotation policies. In the second step, pseudogene fragments and false predictions were filtered out. The last step consisted of screening the remaining spCDSs for lines of extrinsic evidence involving sequence and gene-context comparisons. A total of 789 spCDSs across 20 complete genomes (19 Streptococcus and one Enterococcus) received the support of at least one line of extrinsic evidence, which corresponds to an average of 20 short CDSs per million base pairs. Most of these had no known function, and a significant fraction (31%) are not even annotated as hypothetical genes in GenBank records. As an illustration of the value of this list, we describe a new family of CDSs, encoding very short hydrophobic peptides (20-23 aa) situated just upstream of some of the positive transcriptional regulators of the Rgg family. The expression of seven other short CDSs from Streptococcus thermophilus CNRZ1066 that encode peptides ranging in length from 41 to 56 aa was confirmed by real-time quantitative RT-PCR and revealed a variety of expression patterns. Finally, one peptide from this list, encoded by a gene that is not annotated in GenBank, was identified in a cell-envelope-enriched fraction of S. thermophilus CNRZ1066.

Characterization of a small PlcR-regulated gene co-expressed with cereolysin O

Background

In the human pathogen Bacillus cereus, the expression of most extracellular virulence factors is controlled by the transcriptional activator PlcR. Among these virulence factors, cereolysin O (Clo) is an haemolysin belonging to the cholesterol-dependant cytolysins, a protein family extensively studied in Gram-positive bacteria.

Results

In the genomes of bacteria belonging to the B. cereus group, including Bacillus anthracis and Bacillus thuringiensis, a small gene encoding a 26-amino acid peptide was present in multicopy. One copy was always found upstream from the gene encoding Clo. In B. cereus ATCC 14579, the small gene and the clo gene are co-transcribed. Transcriptional fusions showed that the three paralogues identified in this strain were expressed in a PlcR-dependent manner. We propose to name these peptides Spp for small PlcR-regulated peptides. We show that a synthetic peptide corresponding to the deduced product of the spp genes displayed antibacterial activity.

Conclusion

The co-expression of spp, a small PlcR-regulated multicopy gene with clo suggests a yet unidentified relationship between Spp and the cholesterol-dependent cytolysin in bacteria belonging to the B.cereus group.

Conservation of genes and processes controlled by the quorum response in bacteria: characterization of genes controlled by the quorum-sensing transcription factor ComA in Bacillus subtilis

Quorum or diffusion responses in bacteria are mediated by secreted signalling molecules that accumulate extracellularly as cultures grow to high density. The regulatory response to these signalling molecules can result in dramatic changes in gene expression. In Bacillus subtilis, a quorum response is mediated by a secreted 10-amino-acid modified peptide (ComX pheromone) that activates a receptor histidine kinase (ComP) that activates a response regulator transcription factor (ComA). We have used DNA microarrays to identify genes controlled by the ComX-ComP-ComA quorum-sensing pathway. We found that ComX, ComP and ComA affect the same set of genes, indicating that the kinase ComP is the only receptor for the signalling molecule ComX, and that ComA is the only transcription factor activated directly by ComP, under the conditions tested. Expression of over 20 genes appears to be controlled directly by this signalling pathway, and expression of over 150 additional genes, including those involved in competence development, appears to be controlled indirectly. The genes affected appear to have three general functions: (i) to co-ordinate physiological changes involved in developmental pathways, (ii) to produce extracellular products under conditions in which high concentrations of product are needed to be effective and (iii) to enhance survival, growth and colonization under conditions of crowding or limited diffusion. Many of the genes and processes controlled by the quorum response in B. subtilis are also regulated by quorum sensing in Gram-positive and Gram-negative bacteria. The quorum-sensing signalling molecules and regulatory proteins are quite different between Gram-positives and Gram-negatives and the convergent physiological regulation of similar genes and processes indicate the important and conserved nature of the quorum response.

Genes involved in SkfA killing factor production protect a Bacillus subtilis lipase against proteolysis

Small lipases of Bacillus species, such as LipA from Bacillus subtilis, have a high potential for industrial applications. Recent studies showed that deletion of six AT-rich islands from the B. subtilis genome results in reduced amounts of extracellular LipA. Here we demonstrate that the reduced LipA levels are due to the absence of four genes, skfABCD, located in the prophage 1 region. Intact skfABCD genes are required not only for LipA production at wild-type levels by B. subtilis 168 but also under conditions of LipA overproduction. Notably, SkfA has bactericidal activity and, probably, requires the SkfB to SkfD proteins for its production. The present results show that LipA is more prone to proteolytic degradation in the absence of SkfA and that high-level LipA production can be improved significantly by employing multiple protease-deficient B. subtilis strains. In conclusion, our findings imply that SkfA protects LipA, directly or indirectly, against proteolytic degradation. Conceivably, SkfA could act as a modulator in LipA folding or as a protease inhibitor.