Decoding dissociation of sequence-specific protein-DNA complexes with non-equilibrium simulations

Sequence-specific protein-DNA interactions are crucial in processes such as DNA organization, gene regulation and DNA replication. Obtaining detailed insights into the recognition mechanisms of protein-DNA complexes through experiments is hampered by a lack of resolution in both space and time. Here, we present a molecular simulation approach to quantify the sequence specificity of protein-DNA complexes, that yields results fast, and is generally applicable to any protein-DNA complex. The approach is based on molecular dynamics simulations in combination with a sophisticated steering potential and results in an estimate of the free energy difference of dissociation. We provide predictions of the nucleotide specific binding affinity of the minor groove binding Histone-like Nucleoid Structuring (H-NS) protein, that are in agreement with experimental data. Furthermore, our approach offers mechanistic insight into the process of dissociation. Applying our approach to the major groove binding ETS domain in complex with three different nucleotide sequences identified the high affinity consensus sequence, quantitatively in agreement with experiments. Our protocol facilitates quantitative prediction of protein-DNA complex stability, while also providing high resolution insights into recognition mechanisms. As such, our simulation approach has the potential to yield detailed and quantitative insights into biological processes involving sequence-specific protein-DNA interactions. Graphical Abstract