Designed architectural proteins that tune DNA looping in bacteria

Coarse-grained modeling reveals the impact of supercoiling and loop length in DNA looping kinetics

Measurements of protein-mediated DNA looping reveal that in vivo conditions favor the formation of loops shorter than those that occur in vitro, yet the precise physical mechanisms underlying this shift remain unclear. To understand the extent to which in vivo supercoiling may explain these shifts, we develop a theoretical model based on coarse-grained molecular simulation and analytical transition state theory, enabling us to map out looping energetics and kinetics as a function of two key biophysical parameters: superhelical density and loop length. We show that loops on the scale of a persistence length respond to supercoiling over a much wider range of superhelical densities and to a larger extent than longer loops. This effect arises from a tendency for loops to be centered on the plectonemic end region, which bends progressively more tightly with superhelical density. This trend reveals a mechanism by which supercoiling favors shorter loop lengths. In addition, our model predicts a complex kinetic response to supercoiling for a given loop length, governed by a competition between an enhanced rate of looping due to torsional buckling and a reduction in looping rate due to chain straightening as the plectoneme tightens at higher superhelical densities. Together, these effects lead to a flattening of the kinetic response to supercoiling within the physiological range for all but the shortest loops. Using experimental estimates for in vivo superhelical densities, we discuss our model's ability to explain available looping data, highlighting both the importance of supercoiling as a regulatory force in genetics and the additional complexities of looping phenomena in vivo.

emDNA - A Tool for Modeling Protein-decorated DNA Loops and Minicircles at the Base-pair Step Level

Computational modeling of nucleic acids plays an important role in molecular biology, enhancing our general understanding of the relationship between structure and function. Biophysical studies have provided a wealth of information on how double-helical DNA responds to proteins and other molecules in its local environment but far less understanding of the larger scale structural responses found in protein-decorated loops and minicircles. Current computational models of DNA range from detailed all-atom molecular dynamics studies, which produce rich time and spatially dependent depictions of small DNA fragments, to coarse-grained simulations, which sacrifice detailed physical and chemical information to treat larger scale systems. The treatment of DNA used here, at the base-pair step level with rigid-body parameters, allows one to develop quality models hundreds of base pairs long from local, sequence-specific features found from experiment. The emDNA software takes advantage of this framework, producing elastically optimized models of DNA at thermal equilibrium with built-in or user-generated elastic models. This versatile program, in combination with case studies included in this article, allows users of any skill level to develop and investigate mesoscale models of their own design. The functionality of emDNA includes a tool to incorporate experiment-specific configurations, e.g. protein-bound and/or melted DNA from known high-resolution structures, within higher-order 3D models by fixing the orientation and position of user-specified base pairs. The software provides a new avenue into multiscale genetic modeling, giving a wide range of users a deeper understanding of DNA mesoscale organization and the opportunity to pose new questions in genetic research. The publicly available emDNA software, including build instructions and usage information, is available on GitHub (https://nicocvn.github.io/emDNA/).