RNA Polymerase Binding Site Research Articles

Fidaxomicin resistance in Clostridioides difficile: a systematic review and predictive modeling with RNA polymerase binding sites.

Fidaxomicin (FDX), an RNA polymerase (RNAP) inhibitor antibiotic, is a guideline-recommended therapy for Clostridioides difficile infection. Mutations associated with reduced FDX minimum inhibitory concentrations (MICs) have been identified. However, the molecular characterization of these mutations on FDX binding and the development of FDX resistance have not been studied. The purpose of this systematic review was to identify FDX resistance in C. difficile isolates and determine whether single nucleotide polymorphisms associated with increased FDX MIC aligned with the RNAP binding pocket interacting residues. A systematic literature search was done in PubMed (1991-2023) with identified articles and their bibliographies searched for papers that included C. difficile genetic mutations and increased FDX MIC. Visualization of FDX-RNAP interactions was performed on Schrödinger Maestro using the publicly available C. difficile RNAP with fidaxomicin sequence (code 7L7B) on the Protein Data Bank. Seven articles were identified after applying inclusion and exclusion criteria. The most common mutation in clinical and laboratory isolates was at position V1143 of the β subunit, which accounted for approximately 50% of the identified mutations. Most other mutations occurred within the β' subunit of RNAP. Approximately one-third of the identified mutation aligned directly with FDX interacting residues with C. difficile RNAP (7/20) with most of the remainder occurring within 5 Å of the binding residues. C. difficile strains with elevated FDX MIC align closely with the known RNAP binding residues. These data demonstrate the potential to identify genomic methods to identify emerging FDX resistance.

Visualizing Viral RNA Packaging Signals in Action

Here we confirm, using genome-scale RNA fragments in assembly competition assays, that multiple sub-sites (Packaging Signals, PSs) across the 5′ two-thirds of the gRNA of Satellite Tobacco Necrosis Virus-1 make sequence-specific contacts to the viral CPs helping to nucleate formation of its T = 1 virus-like particle (VLP). These contacts explain why natural virions only package their positive-sense genomes. Asymmetric cryo-EM reconstructions of these VLPs suggest that interactions occur between amino acid residues in the N-terminal ends of the CP subunits and the gRNA PS loop sequences. The base-paired stems of PSs also act non-sequence-specifically by electrostatically promoting the assembly of CP trimers. Importantly, alterations in PS-CP affinity result in an asymmetric distribution of bound PSs inside VLPs, with fuller occupation of the higher affinity 5′ PS RNAs around one vertex, decreasing to an RNA-free opposite vertex within the VLP shell. This distribution suggests that gRNA folding regulates cytoplasmic genome extrusion so that the weakly bound 3′ end of the gRNA, containing the RNA polymerase binding site, extrudes first. This probably occurs after cation-loss induced swelling of the CP-shell, weakening contacts between CP subunits. These data reveal for the first time in any virus how differential PS folding propensity and CP affinities support the multiple roles genomes play in virion assembly and infection. The high degree of conservation between the CP fold of STNV-1 and those of the CPs of many other viruses suggests that these aspects of genome function will be widely shared.

Implementation of computational strategies to combat drug resistance in Mycobacterium tuberculosis

Mycobacterium tuberculosis (MTB) is an infectious pathogen that causes tuberculosis (TB). Even though it is curable, the first-line anti-TB medication rifampicin has been a source of rising concern for human health across the globe. However, rifampicin resistance results from mutations in the RNA polymerase's binding site of the target protein. This infectious pathogen can survive in macrophages and can induce a long-lasting latent infection. The prolonged survival of mycobacterial species was achieved with the aid of serine/threonine protein kinase G (PknG) which performs a preventive role in this situation. Thus, PknG is believed to be an essential target to prevent the pathogen from entering the latency stage. The present investigation was oriented towards molecular-based virtual screening carried out using AutoDock for the FDA-approved and nutraceutical subsets of compounds from the DrugBank repository. Initially, a total of 3035 compounds were screened based on their absorption, distribution, metabolism, excretion and toxicity (ADMET) properties. The screened molecules were taken for docking purposes after which post-docking analysis was performed by calculating the random forest (RF) score. Further, the structural characterization of hit compounds was validated using interaction studies and in silico bioactivity prediction studies. Finally, the anti-myotubular activity prediction was carried out in the search for potential inhibitors. Importantly, four potential ligands, namely DB00080, DB00807, DB06249 and DB06274 were preferred after a thorough investigation. These ligands showed significant structural and protein-ligand interactions, suggesting that they might be potentially used as therapeutic candidates to combat MTB.

High-resolution landscape of an antibiotic binding site

Antibiotic binding sites are located in important domains of essential enzymes and have been extensively studied in the context of resistance mutations; however, their study is limited by positive selection. Using multiplex genome engineering1 to overcome this constraint, we generate and characterize a collection of 760 single-residue mutants encompassing the entire rifampicin binding site of Escherichiacoli RNA polymerase (RNAP). By genetically mapping drug–enzyme interactions, we identify an alpha helix where mutations considerably enhance or disrupt rifampicin binding. We find mutations in this region that prolong antibiotic binding, converting rifampicin from a bacteriostatic to bactericidal drug by inducing lethal DNA breaks. The latter are replication dependent, indicating that rifampicin kills by causing detrimental transcription–replication conflicts at promoters. We also identify additional binding site mutations that greatly increase the speed of RNAP.Fast RNAP depletes the cell of nucleotides, alters cell sensitivity to different antibiotics and provides a cold growth advantage. Finally, by mapping natural rpoB sequence diversity, we discover that functional rifampicin binding site mutations that alter RNAP properties or confer drug resistance occur frequently in nature.

Identifying Innate Resistance Hotspots for SARS-CoV-2 Antivirals Using In Silico Protein Techniques.

The development and approval of antivirals against SARS-CoV-2 has further equipped clinicians with treatment strategies against the COVID-19 pandemic, reducing deaths post-infection. Extensive clinical use of antivirals, however, can impart additional selective pressure, leading to the emergence of antiviral resistance. While we have previously characterized possible effects of circulating SARS-CoV-2 missense mutations on proteome function and stability, their direct effects on the novel antivirals remains unexplored. To address this, we have computationally calculated the consequences of mutations in the antiviral targets: RNA-dependent RNA polymerase and main protease, on target stability and interactions with their antiviral, nucleic acids, and other proteins. By analyzing circulating variants prior to antiviral approval, this work highlighted the inherent resistance potential of different genome regions. Namely, within the main protease binding site, missense mutations imparted a lower fitness cost, while the opposite was noted for the RNA-dependent RNA polymerase binding site. This suggests that resistance to nirmatrelvir/ritonavir combination treatment is more likely to occur and proliferate than that to molnupiravir. These insights are crucial both clinically in drug stewardship, and preclinically in the identification of less mutable targets for novel therapeutic design.

A CRISPRi-based investigation reveals that multiple promoter elements drive gene expression from the genome of mycobacteriophage D29.

A unique feature found in the genomes of mycobacteriophages such as L5 belonging to the A cluster is the presence of multiple dispersed repeated elements known as stoperators. The phage repressor binds these repeat elements, shutting off transcription globally and thereby promoting lysogeny. Interestingly, the sequence of these stoperators closely matches that of the consensus -35 region of prokaryotic promoters, leading us to propose that they may have a role to play in the initiation of transcription by serving as RNA polymerase binding sites. Mycobacteriophage D29 is closely related to phage L5, and their genome organizations are very similar. As in L5, there are multiple stoperators in the genome of D29. The positions occupied by the stoperators in the two genomes are almost identical. The significant difference between the two phages is that D29 lacks the gene encoding the equivalent of the L5 repressor. Since phage D29 does not produce a repressor, we considered it to be a suitable model for testing our hypothesis that the stoperators function as promoters in the absence of the repressor. To prove our point, we targeted CRISPR guide RNAs against six stoperators. In the case of five out of the six, we found a significant reduction in downstream gene expression and phage growth. Based on this observation and primer extension assays, we conclude that promoting gene expression is likely to be the primary function of stoperators.

Targeting Transcriptional and Translational Hindrances in a Modular T7RNAP Expression System in Engineered Pseudomonas putida.

The T7 RNA polymerase is considered one of the most popular tools for heterologous gene expression in the gold standard biotechnological host Escherichia coli. However, the exploitation of this tool in other prospective hosts, such as the biotechnologically relevant bacterium Pseudomonas putida, is still very scarce. The majority of the existing T7-based systems in P.putida show low expression strengths and possess only weak controllability. A fundamental understanding of these systems is necessary in order to design robust and predictable biotechnological processes. To fill this gap, we established and characterized a modular T7 RNA polymerase-based system for heterologous protein production in P.putida, using the enhanced Green Fluorescent Protein (eGFP) as an easy-to-quantify reporter protein. We have effectively targeted the limitations associated with the initial genetic setup of the system, such as slow growth and low protein production rates. By replacing the T7 phage-inherent TΦ terminator downstream of the heterologous gene with the synthetic tZ terminator, growth and protein production rates improved drastically, and the T7 RNA polymerase system reached a productivity level comparable to that of an intrinsic RNA polymerase-based system. Furthermore, we were able to show that the system was saturated with T7 RNA polymerase by applying a T7 RNA polymerase ribosome binding site library to tune heterologous protein production. This saturation indicates an essential role for the ribosome binding sites of the T7 RNA polymerase since, in an oversaturated system, cellular resources are lost to the synthesis of unnecessary T7 RNA polymerase. Eventually, we combined the experimental data into a model that can predict the eGFP production rate with respect to the relative strength of the ribosome binding sites upstream of the T7 gene.

Structural and Functional Analysis of Toxin and Small RNA Gene Promoter Regions in Bacillus anthracis.

It was previously demonstrated that anthrax toxin activator (AtxA) binds directly to the σA-like promoter region of pagA (encoding protective antigen, PA) immediately upstream of the RNA polymerase binding site. In this study, using electrophoretic mobility shift assays and in vivo analyses, we identified AtxA-binding sites in the promoter regions of the lef and cya genes (encoding lethal and edema factors, respectively) and of two Bacillus anthracis small RNAs (XrrA and XrrB). Activities of all four newly studied promoters were enhanced in the presence of CO2/bicarbonate and AtxA, as previously seen for the pagA promoter. Notably, the cya promoter was less activated by AtxA and CO2/bicarbonate conditions. The putative promoter of a recently described third small RNA, XrrC, showed a negligible response to AtxA and CO2/bicarbonate. RNA polymerase binding sites of the newly studied promoters show no consensus and differ from the σA-like promoter region of pagA. In silico analysis of the probable AtxA binding sites in the studied promoters revealed several palindromes. All the analyzed palindromes showed very little overlap with the σA-like pagA promoter. It remains unclear as to how AtxA and DNA-dependent RNA-polymerase identify such diverse DNA-sequences and differentially regulate promoter activation of the studied genes. IMPORTANCE Anthrax toxin activator (AtxA) is the major virulence regulator of Bacillus anthracis, the causative agent of anthrax. Understanding AtxA's mechanism of regulation could facilitate the development of therapeutics for B. anthracis infection. We provide evidence that AtxA binds to the promoters of the cya, lef, xrrA, and xrrB genes. In vivo assays confirmed the activities of all four promoters were enhanced in the presence of AtxA and CO2/bicarbonate, as previously seen for the pagA promoter. The cya and lef genes encode important toxin components. The xrrA and xrrB genes encode sRNAs with a suggested function as cell physiology regulators. Our data provides further evidence for the direct regulatory role of AtxA that was previously shown with the pagA promoter.

De Novo Design of the ArsR Regulated Pars Promoter Enables a Highly Sensitive Whole-Cell Biosensor for Arsenic Contamination.

Whole-cell biosensors for arsenic contamination are typically designed based on natural bacterial sensing systems, which are often limited by their poor performance for precisely tuning the genetic response to environmental stimuli. Promoter design remains one of the most important approaches to address such issues. Here, we use the arsenic-responsive ArsR-Pars regulation system from Escherichia coli MG1655 as the sensing element and coupled gfp or lacZ as the reporter gene to construct the genetic circuit for characterizing the refactored promoters. We first analyzed the ArsR binding site and a library of RNA polymerase binding sites to mine potential promoter sequences. A set of tightly regulated Pars promoters by ArsR was designed by placing the ArsR binding sites into the promoter’s core region, and a novel promoter with maximal repression efficiency and optimal fold change was obtained. The fluorescence sensor PlacV-ParsOC2 constructed with the optimized ParsOC2 promoter showed a fold change of up to 63.80-fold (with green fluorescence visible to the naked eye) at 9.38 ppb arsenic, and the limit of detection was as low as 0.24 ppb. Further, the optimized colorimetric sensor PlacV-ParsOC2-lacZ with a linear response between 0 and 5 ppb was used to perform colorimetric reactions in 24-well plates combined with a smartphone application for the quantification of the arsenic level in groundwater. This study offers a new approach to improve the performance of bacterial sensing promoters and will facilitate the on-site application of arsenic whole-cell biosensors.

A detailed investigation to study the pattern of the interplay of Cyclic AMP Receptor Protein (CRP) of E. coli with its different classes of promoters

The activity of most of the promoters in Escherichia coli, involved in the metabolism of sugars other than glucose, is controlled by a CRP (cAMP receptor protein) or CAP (catabolite activator protein). CRP-dependent promoters are differentiated into various classes (Class I, Class II, and Class III) based on its cognate binding site’s position on DNA. The promoters regulated by CAP are differentially regulated by this transcriptional factor and it is also imperative to mention that these promoters vary greatly in respect to the binding site of CAP to its cognate binding site, it has also been reported that either it overlaps with the binding site of RNA polymerase or it present upstream to it. In Class I CAP-dependent promoters, a particular CAP molecule makes protein-protein interaction for the start of transcription. In Class II CAP-dependent promoters, a particular CAP molecule makes multiple interactions for the start of transcription. At last, in Class III-CAP dependent promoters, more than one CAP molecule is involved and activation of transcription is done synergistically. It has also been documented that CAP shows a kind of biphasic behavior in some promoters. So, the main focus of this work is to find out whether this biphasic behavior is true for other E. coli promoters as well. Experiments have been performed to know more about this biphasic nature and the various patterns of interactions of catabolite activator protein (CAP) of E. coli with its different classes of promoters.

Inhibition of the RNA-dependent RNA Polymerase of the SARS-CoV-2 by Short Peptide Inhibitors

The rapid proliferation of SARS-CoV-2 in COVID-19 patients has become detrimental to their lives. However, blocking the replication cycle of SARS-CoV-2 will help in suppressing the viral loads in patients, which would ultimately help in the early recovery. To discover such drugs, molecular docking, MD-simulations, and MM/GBSA approaches have been used herein to examine the role of several short ionic peptides in inhibiting the RNA binding site of the RNA-dependent RNA polymerase (RdRp). Out of the 49 tri- and tetrapeptide inhibitors studied, 8 inhibitors were found to bind RdRp strongly as revealed by the docking studies. Among these inhibitors, the Ala1-Arg2-Lys3-Asp4 and Ala1-Lys2-Lys3-Asp4 are found to make the most stable complexes with RdRp and possess the ΔGbind of -17.41 and -14.21 kcal/mol respectively as revealed by the MD and MM/GBSA studies. Hence these peptide inhibitors would be highly potent in inhibiting the activities of RdRp. It is further found that these inhibitors can occupy the positions of the nucleotide triphosphate (NTP) insertion site, thereby inhibiting the replication of the viral genome by obstructing the synthesis of new nucleotides. Structural and energetic comparisons of these inhibitors with Remdesivir and similar nucleotide drugs show that these peptides would be more specific and hence may act as promiscuous antiviral agents against RdRp.

The Transcription Factor Rv1453 Regulates the Expression of qor and Confers Resistant to Clofazimine in Mycobacterium tuberculosis

ObjectiveClofazimine plays an important role in the treatment of drug-resistant tuberculosis. However, the mechanism of clofazimine resistance remains unclear. In order to slow down the occurrence of clofazimine resistance, it is necessary to study its resistance mechanism.MethodsIn this study, we constructed Rv1453 knockout, complementary and overexpressed strain. The minimum inhibitory concentration (MIC) of clofazimine against Mycobacterium tuberculosis was detected by microplate alamar blue assay (MABA). The transcription levels of Rv1453 and its adjacent genes were detected by quantitative reverse transcriptase PCR. The purified Rv1453 protein was used for electrophoretic mobility shift assay (EMSA) to identify the binding site of Rv1453 protein.ResultsThe minimum inhibitory concentration (MIC) of clofazimine increased about 4-fold for the Rv1453 knockout strain and decreased about 4-fold for the Rv1453 overexpressed strain compared with Mycobacterium tuberculosis H37Rv. Further analysis showed that Rv1453 protein, as a regulatory protein, binds to the RNA polymerase binding site of qor and blocks the transcription process.ConclusionThis study preliminarily revealed that Rv1453 protein of Mycobacterium tuberculosis affects its susceptibility to clofazimine by regulating the transcription level of qor, which is shedding a new light on the mechanism of clofazimine resistance.

Effects of Synonymous Mutations beyond Codon Bias: The Evidence for Adaptive Synonymous Substitutions from Microbial Evolution Experiments.

Synonymous mutations are often assumed to be neutral with respect to fitness because they do not alter the encoded amino acid and so cannot be “seen” by natural selection. Yet a growing body of evidence suggests that synonymous mutations can have fitness effects that drive adaptive evolution through their impacts on gene expression and protein folding. Here, we review what microbial experiments have taught us about the contribution of synonymous mutations to adaptation. A survey of site-directed mutagenesis experiments reveals the distributions of fitness effects for nonsynonymous and synonymous mutations are more similar, especially for beneficial mutations, than expected if all synonymous mutations were neutral, suggesting they should drive adaptive evolution more often than is typically observed. A review of experimental evolution studies where synonymous mutations have contributed to adaptation shows they can impact fitness through a range of mechanisms including the creation of illicit RNA polymerase binding sites impacting transcription and changes to mRNA folding stability that modulate translation. We suggest that clonal interference in evolving microbial populations may be the reason synonymous mutations play a smaller role in adaptive evolution than expected based on their observed fitness effects. We finish by discussing the impacts of falsely assuming synonymous mutations are neutral and discuss directions for future work exploring the role of synonymous mutations in adaptive evolution.

Multiplexed characterization of rationally designed promoter architectures deconstructs combinatorial logic for IPTG-inducible systems

A crucial step towards engineering biological systems is the ability to precisely tune the genetic response to environmental stimuli. In the case of Escherichia coli inducible promoters, our incomplete understanding of the relationship between sequence composition and gene expression hinders our ability to predictably control transcriptional responses. Here, we profile the expression dynamics of 8269 rationally designed, IPTG-inducible promoters that collectively explore the individual and combinatorial effects of RNA polymerase and LacI repressor binding site strengths. We then fit a statistical mechanics model to measured expression that accurately models gene expression and reveals properties of theoretically optimal inducible promoters. Furthermore, we characterize three alternative promoter architectures and show that repositioning binding sites within promoters influences the types of combinatorial effects observed between promoter elements. In total, this approach enables us to deconstruct relationships between inducible promoter elements and discover practical insights for engineering inducible promoters with desirable characteristics.

Genetic Variability of Tick-Borne Encephalitis Virus Genome 5'-UTR from Northern Eurasia

This paper reports the analysis of the nucleotide sequences of the 5'-untranslated region (5'-UTR) of tick-borne encephalitis virus (TBEV) genomic RNA isolated from 39 individual taiga ticks collected in several regions of Northern Eurasia. The sequences of 5'-UTRs of the Siberian and Far East TBEV genotypes were 89% and 95% identical to the prototype strains (Zausaev and 205), respectively. The detected nucleotide substitutions were typical for these two TBEV genotypes, which made possible unambiguous identification. Both conservative and variable motifs were detected in the 5'-UTR RNA. The B2, C1, and C2 elements of the Y-shaped 5'-UTR structure and the presumable viral RNA-dependent RNA-polymerase binding site were the most variable. The A2, CS A, CS В elements as well as the start codon were conservative. Interestingly, five substitutions in the 5'-UTR C1 variable element of the TBEVs isolated in different geographical regions were strictly conservative, while 11 different substitutions were detected in this element among the laboratory TBEV variants. A little less that a third of all nucleotide substitutions were mapped outside the main elements of the Y-shaped structure. In general, nucleotide substitutions were localized to stem structures, not being found in the hairpin regions of the TBEV 5'-UTR. The results indicated significant variability of the genomic RNA 5'-UTR in the TBEV laboratory strains and field isolates obtained from different geographical regions. It has been suggested that genetic variability of 5'-UTR is characteristic of the TBEV genome 5'-UTR organization and may serve as a structural basis for virus efficient replication in various avian, mammalian, and ixodic tick cells.

Deep learning approach for predicting functional Z-DNA regions using omics data

Computational methods to predict Z-DNA regions are in high demand to understand the functional role of Z-DNA. The previous state-of-the-art method Z-Hunt is based on statistical mechanical and energy considerations about B- to Z-DNA transition using sequence information. Z-DNA CHiP-seq experiment results showed little overlap with Z-Hunt predictions implying that sequence information only is not sufficient to explain emergence of Z-DNA at different genomic locations. Adding epigenetic and other functional genomic mark-ups to DNA sequence level can help revealing the functional Z-DNA sites. Here we take advantage of the deep learning approach that can analyze and extract information from large volumes of molecular biology data. We developed a machine learning approach DeepZ that aggregates information from genome-wide maps of epigenetic markers, transcription factor and RNA polymerase binding sites, and chromosome accessibility maps. With the developed model we not only verify the experimental Z-DNA predictions, but also generate the whole-genome annotation, introducing new possible Z-DNA regions, which have not yet been found in experiments and can be of interest to the researchers from various fields.

A highly rifampicin resistant Mycobacterium tuberculosis strain emerging in Southern Brazil

Here we described phenotypical, molecular and epidemiological features of a highly rifampicin-resistant Mycobacterium tuberculosis strain emerging in Southern Brazil, that carries an uncommon insertion of 12 nucleotides at the codon 435 in the rpoB gene. Employing a whole-genome sequencing-based study on drug-resistant Mycobacterium tuberculosis strains, we identified this emergent strain in 16 (9.19%) from 174 rifampicin-resistant clinical strains, all of them belonging to LAM RD115 sublineage. Nine of these 16 strains were available to minimum inhibitory concentration determination and for all of them was found a high rifampicin-resistance level (≥to 32 mg/L). This high resistance level could be explained by structural changes into the RIF binding site of RNA polymerase caused by the insertions, and consequent low-affinity interaction with rifampicin complex confirmed through protein modeling and molecular docking simulations. Epidemiological investigation showed that most of the individuals (56.25%) infected by the studied strains were prison inmate individuals or that spent some time in prison. The phylogenomic approach revealed that strains carrying on insertion belonged to same genomic cluster, evidencing a communal transmission chain involving inmate individuals and community. We stress the importance of tuberculosis genomic surveillance and introduction of measures to interrupt Mycobacterium tuberculosis transmission chain in this region.

Benzyl and benzoyl benzoic acid inhibitors of bacterial RNA polymerase-sigma factor interaction.

Transcription is an essential biological process in bacteria requiring a core enzyme, RNA polymerase (RNAP). Bacterial RNAP is catalytically active but requires sigma (σ) factors for transcription of natural DNA templates. σ factor binds to RNAP to form a holoenzyme which specifically recognizes a promoter, melts the DNA duplex, and commences RNA synthesis. Inhibiting the binding of σ to RNAP is expected to inhibit bacterial transcription and growth. We previously identified a triaryl hit compound that mimics σ at its major binding site of RNAP, thereby inhibiting the RNAP holoenzyme formation. In this study, we modified this scaffold to provide a series of benzyl and benzoyl benzoic acid derivatives possessing improved antimicrobial activity. A representative compound demonstrated excellent activity against Staphylococcus epidermidis with minimum inhibitory concentrations reduced to 0.5μg/mL, matching that of vancomycin. The molecular mechanism of inhibition was confirmed using biochemical and cellular assays. Low cytotoxicity and metabolic stability of compounds demonstrated the potential for further studies.

Ciliate mitoribosome illuminates evolutionary steps of mitochondrial translation.

To understand the steps involved in the evolution of translation, we used Tetrahymena thermophila, a ciliate with high coding capacity of the mitochondrial genome, as the model organism and characterized its mitochondrial ribosome (mitoribosome) using cryo-EM. The structure of the mitoribosome reveals an assembly of 94-ribosomal proteins and four-rRNAs with an additional protein mass of ~700 kDa on the small subunit, while the large subunit lacks 5S rRNA. The structure also shows that the small subunit head is constrained, tRNA binding sites are formed by mitochondria-specific protein elements, conserved protein bS1 is excluded, and bacterial RNA polymerase binding site is blocked. We provide evidence for anintrinsic protein targeting system through visualization of mitochondria-specific mL105 by the exit tunnel that would facilitate the recruitment of a nascent polypeptide. Functional protein uS3m is encoded by three complementary genes from the nucleus and mitochondrion, establishing a link between genetic drift and mitochondrial translation. Finally, we reannotated nine open reading frames in the mitochondrial genome that code for mitoribosomal proteins.

Spo0A can efficiently enhance the expression of the alkaline protease gene aprE in Bacillus licheniformis by specifically binding to its regulatory region

The expression of enzymes in Bacillus licheniformis, such as the valuable extracellular alkaline protease AprE, is highly regulated by a complex transcriptional regulation mechanism. Here, we found that the transcript abundance of aprE varies >343-fold in response to the supply of nutrients or to environmental challenges. To identify the underlying regulatory mechanism, the core promoter of aprE and several important upstream regulatory regions outside the promoter were firstly confirmed by 5′-RACE and mutagenesis experiments. The specific proteins that bind to the identified sequences were subsequently captured by DNA pull-down experiments, which yielded the transcriptional factors (TFs) Spo0A, CggR, FruR, YhcZ, as well as fragments of functionally unassigned proteins. Further electrophoretic mobility shift assay (EMSA) and DNase I foot-printing experiments indicated that Spo0A can directly bind to the region from −92 to −118 nucleotides upstream of the transcription start site, and the deletion of this specific region drastically decreased the production of AprE. Taken together, these results indicated that the expression of aprE was mainly regulated by the interplay between Spo0A and its cognate DNA sequence, which was successfully applied to overproduce AprE in a genetically modified host harboring three aprE expression cassettes. The DNA binding proteins may serve to increase the efficiency of transcription by creating an additional binding site for RNA polymerase. The discovery of this mechanism significantly increases our understanding of the aprE transcription mechanism, which is of great importance for AprE overproduction.