“…In recent years, MD simulations have been widely used to study the flexibility of DNAs, as DNA structural flexibility is closely associated with many biological processes involving the storage or encoding of genetic information [ 65 , 66 ]. Although many results from single-molecule experiments can be well-described by the commonly accepted WLC models [ 27 , 28 , 29 ], atomistic MD simulations are extensively used to obtain microscopic descriptions of DNA flexibility, such as the width and depth of the major/minor grooves and the distances/twist angles between neighbor base pairs [ 67 , 68 , 69 ]. For example, to explain the experimental results that short DNAs consisting of tens of base pairs (bps) may have seemingly higher flexibility than those of kilobase pairs, Wu et al performed MD simulations for short dsDNAs with a finite-length of 5–50 bps using the Amber parmbsc0 force field.…”