Predicting the onset of DNA supercoiling using a non-linear hemitropic elastic rod

Elastic rod models provide a means to interpret single molecule DNA experiments as well as predict DNA behavior under physiological conditions. Here we use an elastic rod model to predict the stability boundary (critical torque vs. applied tension) for single molecule DNA experiments in which the molecule is subjected to applied tension and twist. We discuss the shortcomings of the usual isotropic rod model. We then derive a consistent non-linear material law from the general representation for a hemitropic (chiral) rod. Finally, we present results of a standard bifurcation analysis predicting the stability boundary. We find results from the non-linear hemitropic rod to match the data closely.