BI and BII conformational substates in the DNA backbone typify canonical B-form DNA. The BI and BII substates are important for structural variation of DNA and have been implicated in protein鈥搉ucleic acid recognition mechanisms. Recent refinements have been made to nucleic acid force fields employed in molecular dynamics simulations that demonstrate a better ability to model the BI and BII states, leading to overall improved modeling of DNA structure and dynamics. These force field improvements have yet to be significantly demonstrated in the context of a protein鈥揇NA system extended to long time scales. Our plan was to run molecular dynamics simulations of a well-studied protein鈥揇NA system (E2鈥揇NA) into the microsecond time scale and determine the ability of the force field to populate BII states in the DNA backbone consistent with dinucleotide steps crystallized in the BII conformation. The results showed that the dinucleotide steps in the E2鈥揇NA complex with the highest BII populations from simulation trajectories corresponded to the dinucleotide steps crystallized in the BII state and that decoy BI and BII states converge to the same results within approximately one microsecond.