Monte carlo simulations of protein assembly, disassembly, and linear motion on DNA

We use Monte Carlo simulations to analyze the simultaneous interactions of multiple proteins to a long DNA molecule. We study the time dependence of protein organization on DNA for different regimes that comprise (non)cooperative sequence-independent protein assembly, dissociation, and linear motion. A range of different behaviors is observed for the dynamics, final coverage, and cluster size distributions. We observe that the DNA substrate is almost never completely covered by protein when taking into account only (non)cooperative binding, because gaps remain on the substrate that are smaller than the binding site size of the protein. Due to these gaps, the apparent binding size of a protein during noncooperative binding can be overestimated by up to 30%. During dissociation of cooperatively bound proteins, the dissociation curve can be exponentially shaped even when allowing only end-dependent dissociation. We discuss the potential of our method for the analysis of a number of single-molecule experiments, for example, the binding of the DNA-repair proteins RecA and Rad51 to DNA.