Search for additional replication terminators in the Bacillus subtilis 168 chromosome

An affinity-structure database of helix-turn-helix: DNA complexes with a universal coordinate system

Background

Molecular interactions between proteins and DNA molecules underlie many cellular processes, including transcriptional regulation, chromosome replication, and nucleosome positioning. Computational analyses of protein-DNA interactions rely on experimental data characterizing known protein-DNA interactions structurally and biochemically. While many databases exist that contain either structural or biochemical data, few integrate these two data sources in a unified fashion. Such integration is becoming increasingly critical with the rapid growth of structural and biochemical data, and the emergence of algorithms that rely on the synthesis of multiple data types to derive computational models of molecular interactions.

Description

We have developed an integrated affinity-structure database in which the experimental and quantitative DNA binding affinities of helix-turn-helix proteins are mapped onto the crystal structures of the corresponding protein-DNA complexes. This database provides access to: (i) protein-DNA structures, (ii) quantitative summaries of protein-DNA binding affinities using position weight matrices, and (iii) raw experimental data of protein-DNA binding instances. Critically, this database establishes a correspondence between experimental structural data and quantitative binding affinity data at the single basepair level. Furthermore, we present a novel alignment algorithm that structurally aligns the protein-DNA complexes in the database and creates a unified residue-level coordinate system for comparing the physico-chemical environments at the interface between complexes. Using this unified coordinate system, we compute the statistics of atomic interactions at the protein-DNA interface of helix-turn-helix proteins. We provide an interactive website for visualization, querying, and analyzing this database, and a downloadable version to facilitate programmatic analysis.

Conclusions

This database will facilitate the analysis of protein-DNA interactions and the development of programmatic computational methods that capitalize on integration of structural and biochemical datasets. The database can be accessed at http://ProteinDNA.hms.harvard.edu.

Three enhancements to the inference of statistical protein-DNA potentials

The energetics of protein-DNA interactions are often modeled using so-called statistical potentials, that is, energy models derived from the atomic structures of protein-DNA complexes. Many statistical protein-DNA potentials based on differing theoretical assumptions have been investigated, but little attention has been paid to the types of data and the parameter estimation process used in deriving the statistical potentials. We describe three enhancements to statistical potential inference that significantly improve the accuracy of predicted protein-DNA interactions: (i) incorporation of binding energy data of protein-DNA complexes, in conjunction with their X-ray crystal structures, (ii) use of spatially-aware parameter fitting, and (iii) use of ensemble-based parameter fitting. We apply these enhancements to three widely-used statistical potentials and use the resulting enhanced potentials in a structure-based prediction of the DNA binding sites of proteins. These enhancements are directly applicable to all statistical potentials used in protein-DNA modeling, and we show that they can improve the accuracy of predicted DNA binding sites by up to 21%.

Sequence-directed DNA export guides chromosome translocation during sporulation in Bacillus subtilis

In prokaryotes, the transfer of DNA between cellular compartments is essential for the segregation and exchange of genetic material. SpoIIIE and FtsK are AAA+ ATPases responsible for intercompartmental chromosome translocation in bacteria. Despite functional and sequence similarities, these motors were proposed to use drastically different mechanisms: SpoIIIE was suggested to be a unidirectional DNA transporter that exports DNA from the compartment in which it assembles, whereas FtsK was shown to establish translocation directionality by interacting with highly skewed chromosomal sequences. Here we use a combination of single-molecule, bioinformatics and in vivo fluorescence methodologies to study the properties of DNA translocation by SpoIIIE in vitro and in vivo. These data allow us to propose a sequence-directed DNA exporter model that reconciles previously proposed models for SpoIIIE and FtsK, constituting a unified model for directional DNA transport by the SpoIIIE/FtsK family of AAA+ ring ATPases.

Replication terminator protein-based replication fork-arrest systems in various Bacillus species

The replication terminator protein (RTP) of Bacillus subtilis interacts with its cognate DNA terminators to cause replication fork arrest, thereby ensuring that the forks approaching one another at the conclusion of a round of replication meet within a restricted terminus region. A similar situation exists in Escherichia coli, but it appears that the fork-arrest systems in these two organisms have evolved independently of one another. In the present work, RTP homologs in four species closely related to B. subtilis (B. atrophaeus, B. amyloliquefaciens, B. mojavensis, and B. vallismortis) have been identified and characterized. An RTP homolog could not be identified in another closely related species, B. licheniformis. The nucleotide and amino acid changes from B. subtilis among the four homologs are consistent with the recently established phylogenetic tree for these species. The GC contents of the rtp genes raise the possibility that these organisms arose within this branch of the tree by horizontal transfer into a common ancestor after their divergence from B. licheniformis. Only 5 amino acid residue positions were changed among the four homologs, despite an up to 17.2% change in the nucleotide sequence, a finding that highlights the importance of the precise folded structure to the functioning of RTP. The absence of any significant change in the proposed DNA-binding region of RTP emphasizes the importance of its high affinity for the DNA terminator in its functioning. By coincidence, the single change (E30K) found in the B. mojavensis RTP corresponds exactly to that purposefully introduced by others into B. subtilis RTP to implicate a crucial role for E30 in the fork-arrest mechanism. The natural occurrence of this variant is difficult to reconcile with such an implication, and it was shown directly that RTP.E30K functions normally in fork arrest in B. subtilis in vivo. Additional DNA terminators were identified in the new RTP homolog-containing strains, allowing the definition of a Bacillus terminator consensus and identification of two more terminators in the B. subtilis 168 genome sequence to bring the total to nine.

Analyzing genomes with cumulative skew diagrams

A novel method of cumulative diagrams shows that the nucleotide composition of a microbial chromosome changes at two points separated by about a half of its length. These points coincide with sites of replication origin and terminus for all bacteria where such sites are known. The leading strand is found to contain more guanine than cytosine residues. This fact is used to predict origin and terminus locations in other bacterial and archaeal genomes. Local changes, visible as diagram distortions, may represent recent genome rearrangements, as demonstrated for two strains of Escherichia coli . Analysis of the diagrams of viral and mitochondrial genomes suggests a link between the base composition bias and the time spent by DNA in a single stranded state during replication.