Identification of Position-Specific Correlations between DNA-Binding Domains and Their Binding Sites. Application to the MerR Family of Transcription Factors

Development of a highly sensitive PbrR-based biosensor via directed evolution and its application for lead detection

The Whole-cell biosensor (WCB) is a convenient and practical assay that can monitor bioavailable lead (Pb) contamination. However, existing Pb-responsive WCB struggle to meet practical detection needs due to the lack of sensitivity, specificity, and stability. In this study, we developed a Pb WCB using the Pb resistance transcriptional regulatory factor (PbrR) and green fluorescent protein (GFP), and improved its performance by directed evolution in conjunction with fluorescence-activated cell sorting (FACS). After 3 rounds of screening, we acquired a biosensor mutant (PbrR-E3). The evolved biosensor exhibited an approximately 11-fold increase in maximum fluorescence output signal compared to the non-evolved biosensor, resulting in an improvement of its sensitivity and specificity. This biosensor demonstrated a limit of detection (LOD) of 0.045 μg/L. Furthermore, the evolved biosensor showcased outstanding performance in the detection of Pb(II) in tea infusion and also demonstrated good stability in tests with spiked real water samples. These results highlight the potential of the evolved WCB as a viable approach for monitoring Pb.

Predicting bacterial transcription factor binding sites through machine learning and structural characterization based on DNA duplex stability

Transcriptional factors (TFs) in bacteria play a crucial role in gene regulation by binding to specific DNA sequences, thereby assisting in the activation or repression of genes. Despite their central role, deciphering shape recognition of bacterial TFs-DNA interactions remains an intricate challenge. A deeper understanding of DNA secondary structures could greatly enhance our knowledge of how TFs recognize and interact with DNA, thereby elucidating their biological function. In this study, we employed machine learning algorithms to predict transcription factor binding sites (TFBS) and classify them as directed-repeat (DR) or inverted-repeat (IR). To accomplish this, we divided the set of TFBS nucleotide sequences by size, ranging from 8 to 20 base pairs, and converted them into thermodynamic data known as DNA duplex stability (DDS). Our results demonstrate that the Random Forest algorithm accurately predicts TFBS with an average accuracy of over 82% and effectively distinguishes between IR and DR with an accuracy of 89%. Interestingly, upon converting the base pairs of several TFBS-IR into DDS values, we observed a symmetric profile typical of the palindromic structure associated with these architectures. This study presents a novel TFBS prediction model based on a DDS characteristic that may indicate how respective proteins interact with base pairs, thus providing insights into molecular mechanisms underlying bacterial TFs-DNA interaction.

ResR/McdR-regulated protein translation machinery contributes to drug resilience in Mycobacterium tuberculosis

Survival response of the human tuberculosis pathogen, Mycobacterium tuberculosis (Mtb) to a diverse environmental cues is governed through its versatile transcription regulatory mechanisms with the help of a large pool of transcription regulators (TRs). Rv1830 is one such conserved TR, which remains uncharacterized in Mtb. It was named as McdR based on an effect on cell division upon its overexpression in Mycobacterium smegmatis. Recently, it has been implicated in antibiotic resilience in Mtb and reannotated as ResR. While Rv1830 affects cell division by modulating the expression of M. smegmatis whiB2, the underlying cause of its essentiality and regulation of drug resilience in Mtb is yet to be deciphered. Here we show that ResR/McdR, encoded by ERDMAN_2020 in virulent Mtb Erdman, is pivotal for bacterial proliferation and crucial metabolic activities. Importantly, ResR/McdR directly regulates ribosomal gene expression and protein synthesis, requiring distinct disordered N-terminal sequence. Compared to control, bacteria depleted with resR/mcdR exhibit delayed recovery post-antibiotic treatment. A similar effect upon knockdown of rplN operon genes further implicates ResR/McdR-regulated protein translation machinery in attributing drug resilience in Mtb. Overall, findings from this study suggest that chemical inhibitors of ResR/McdR may be proven effective as adjunctive therapy for shortening the duration of TB treatment.

Chimeric MerR-Family Regulators and Logic Elements for the Design of Metal Sensitive Genetic Circuits in Bacillus subtilis

Whole-cell biosensors are emerging as promising tools for monitoring environmental pollutants such as heavy metals. These sensors constitute a genetic circuit comprising a sensing module and an output module, such that a detectable signal is produced in the presence of the desired analyte. The MerR family of metal-responsive regulators offers great potential for the construction of metal sensing circuits, due to their high sensitivity, tight transcription control, and large diversity in metal-specificity. However, the sensing diversity is broadest in Gram-negative systems, while chassis organisms are often selected from Gram-positive species, particularly sporulating bacilli. This can be problematic, because Gram-negative biological parts, such as promoters, are frequently observed to be nonfunctional in Gram-positive hosts. Herein, we combined construction of synthetic genetic circuits and chimeric MerR regulators, supported by structure-guided design, to generate metal-sensitive biosensor modules that are functional in the biotechnological work-horse species Bacillus subtilis. These chimeras consist of a constant Gram-positive derived DNA-binding domain fused to variable metal binding domains of Gram-negative origins. To improve the specificity of the whole-cell biosensor, we developed a modular "AND gate" logic system based on the B. subtilis two-subunit σ-factor, SigO-RsoA, designed to maximize future use for synthetic biology applications in B. subtilis. This work provides insights into the use of modular regulators, such as the MerR family, in the design of synthetic circuits for the detection of heavy metals, with potentially wider applicability of the approach to other systems and genetic backgrounds.

Comparative Analysis of the IclR-Family of Bacterial Transcription Factors and Their DNA-Binding Motifs: Structure, Positioning, Co-Evolution, Regulon Content

The IclR-family is a large group of transcription factors (TFs) regulating various biological processes in diverse bacteria. Using comparative genomics techniques, we have identified binding motifs of IclR-family TFs, reconstructed regulons and analyzed their content, finding co-occurrences between the regulated COGs (clusters of orthologous genes), useful for future functional characterizations of TFs and their regulated genes. We describe two main types of IclR-family motifs, similar in sequence but different in the arrangement of the half-sites (boxes), with GKTYCRYW_3-4RYGRAMC and TGRAACAN_1-2TGTTYCA consensuses, and also predict that TFs in 32 orthologous groups have binding sites comprised of three boxes with alternating direction, which implies two possible alternative modes of dimerization of TFs. We identified trends in site positioning relative to the translational gene start, and show that TFs in 94 orthologous groups bind tandem sites with 18-22 nucleotides between their centers. We predict protein-DNA contacts via the correlation analysis of nucleotides in binding sites and amino acids of the DNA-binding domain of TFs, and show that the majority of interacting positions and predicted contacts are similar for both types of motifs and conform well both to available experimental data and to general protein-DNA interaction trends.

A genome-wide scan for correlated mutations detects macromolecular and chromatin interactions in Arabidopsis thaliana

The concept of exploiting correlated mutations has been introduced and applied successfully to identify interactions within and between biological macromolecules. Its rationale lies in the preservation of physical interactions via compensatory mutations. With the massive increase of available sequence information, approaches based on correlated mutations have regained considerable attention. We analyzed a set of 10 707 430 single nucleotide polymorphisms detected in 1135 accessions of the plant Arabidopsis thaliana. To measure their covariance and to reveal the global genome-wide sequence correlation structure of the Arabidopsis genome, the adjusted mutual information has been estimated for each possible pair of polymorphic sites. We developed a series of filtering steps to account for genetic linkage and lineage relations between Arabidopsis accessions, as well as transitive covariance as possible confounding factors. We show that upon appropriate filtering, correlated mutations prove indeed informative with regard to molecular interactions, and furthermore, appear to reflect on chromosomal interactions. Our study demonstrates that the concept of correlated mutations can also be applied successfully to within-species sequence variation and establishes a promising approach to help unravel the complex molecular interactions in A. thaliana and other species with broad sequence information.

Deciphering the protein-DNA code of bacterial winged helix-turn-helix transcription factors

Background

Sequence-specific binding by transcription factors (TFs) plays a significant role in the selection and regulation of target genes. At the protein:DNA interface, amino acid side-chains construct a diverse physicochemical network of specific and non-specific interactions, and seemingly subtle changes in amino acid identity at certain positions may dramatically impact TF:DNA binding. Variation of these specificity-determining residues (SDRs) is a major mechanism of functional divergence between TFs with strong structural or sequence homology.

Methods

In this study, we employed a combination of high-throughput specificity profiling by SELEX and Spec-seq, structural modeling, and evolutionary analysis to probe the binding preferences of winged helix-turn-helix TFs belonging to the OmpR sub-family in Escherichia coli.

Results

We found that E. coli OmpR paralogs recognize tandem, variably spaced repeats composed of "GT-A" or "GCT"-containing half-sites. Some divergent sequence preferences observed within the "GT-A" mode correlate with amino acid similarity; conversely, "GCT"-based motifs were observed for a subset of paralogs with low sequence homology. Direct specificity profiling of a subset of OmpR homologues (CpxR, RstA, and OmpR) as well as predicted "SDR-swap" variants revealed that individual SDRs may impact sequence preferences locally through direct contact with DNA bases or distally via the DNA backbone.

Conclusions

Overall, our work provides evidence for a common structural code for sequence-specific wHTH:DNA interactions, and demonstrates that surprisingly modest residue changes can enable recognition of highly divergent sequence motifs. Further examination of SDR predictions will likely reveal additional mechanisms controlling the evolutionary divergence of this important class of transcriptional regulators.

Scoring Targets of Transcription in Bacteria Rather than Focusing on Individual Binding Sites

Reliable identification of targets of bacterial regulators is necessary to understand bacterial gene expression regulation. These targets are commonly predicted by searching for high-scoring binding sites in the upstream genomic regions, which typically leads to a large number of false positives. In contrast to the common approach, here we propose a novel concept, where overrepresentation of the scoring distribution that corresponds to the entire searched region is assessed, as opposed to predicting individual binding sites. We explore two implementations of this concept, based on Kolmogorov-Smirnov (KS) and Anderson-Darling (AD) tests, which both provide straightforward P-value estimates for predicted targets. This approach is implemented for pleiotropic bacterial regulators, including σ70 (bacterial housekeeping σ factor) target predictions, which is a classical bioinformatics problem characterized by low specificity. We show that KS based approach is both faster and more accurate, departing from the current paradigm of AD being slower, but more accurate. Moreover, KS approach leads to a significant increase in the search accuracy compared to the standard approach, while at the same time straightforwardly assigning well established P-values to each potential target. Consequently, the new KS based method proposed here, which assigns P-values to fixed length upstream regions, provides a fast and accurate approach for predicting bacterial transcription targets.

Two-domain analysis of JmjN-JmjC and PHD-JmjC lysine demethylases: Detecting an inter-domain evolutionary stress

Residues at different positions of a multiple sequence alignment sometimes evolve together, due to a correlated structural or functional stress at these positions. Co-evolution has thus been evidenced computationally in multiple proteins or protein domains. Here, we wish to study whether an evolutionary stress is exerted on a sequence alignment across protein domains, i.e., on longer sequence separations than within a single protein domain. JmjC-containing lysine demethylases were chosen for analysis, as a follow-up to previous studies; these proteins are important multidomain epigenetic regulators. In these proteins, the JmjC domain is responsible for the demethylase activity, and surrounding domains interact with histones, DNA or partner proteins. This family of enzymes was analyzed at the sequence level, in order to determine whether the sequence of JmjC-domains was affected by the presence of a neighboring JmjN domain or PHD finger in the protein. Multiple positions within JmjC sequences were shown to have their residue distributions significantly altered by the presence of the second domain. Structural considerations confirmed the relevance of the analysis for JmjN-JmjC proteins, while among PHD-JmjC proteins, the length of the linker region could be correlated to the residues observed at the most affected positions. The correlation of domain architecture with residue types at certain positions, as well as that of overall architecture with protein function, is discussed. The present results thus evidence the existence of an across-domain evolutionary stress in JmjC-containing demethylases, and provide further insights into the overall domain architecture of JmjC domain-containing proteins.