Lac repressor mediated DNA looping: Monte Carlo simulation of constrained DNA molecules complemented with current experimental results

Tethered particle motion (TPM) experiments can be used to detect time-resolved loop formation in a single DNA molecule by measuring changes in the length of a DNA tether. Interpretation of such experiments is greatly aided by computer simulations of DNA looping which allow one to analyze the structure of the looped DNA and estimate DNA-protein binding constants specific for the loop formation process. We here present a new Monte Carlo scheme for accurate simulation of DNA configurations subject to geometric constraints and apply this method to Lac repressor mediated DNA looping, comparing the simulation results with new experimental data obtained by the TPM technique. Our simulations, taking into account the details of attachment of DNA ends and fluctuations of the looped subsegment of the DNA, reveal the origin of the double-peaked distribution of RMS values observed by TPM experiments by showing that the average RMS value for anti-parallel loop types is smaller than that of parallel loop types. The simulations also reveal that the looping probabilities for the anti-parallel loop types are significantly higher than those of the parallel loop types, even for loops of length 600 and 900 base pairs, and that the correct proportion between the heights of the peaks in the distribution can only be attained when loops with flexible Lac repressor conformation are taken into account. Comparison of the in silico and in vitro results yields estimates for the dissociation constants characterizing the binding affinity between O1 and Oid DNA operators and the dimeric arms of the Lac repressor. © 2014 Biton et al

Modulating DNA configuration by interfacial traction: an elastic rod model to characterize DNA folding and unfolding

As a continuum model of DNA, a thin elastic rod subjected to interfacial interactions is used to investigate the equilibrium configuration of DNA in intracellular solution. The interfacial traction between the rod and the solution environment is derived in detail. Kirchhoff’s theory of elastic rods is used to analyze the equilibrium configuration of a DNA segment under the action of the interfacial traction. The influences of the interfacial energy factor and bending stiffness on the toroidal spool formation of the DNA segment are discussed. The results show that the equilibrium configuration of DNA is mainly determined by competition between the interfacial energy and elastic strain energy of the DNA itself, and the interfacial traction is one of the forces that drives DNA folding and unfolding