Monte Carlo simulation of melting transition on DNA nanocompartment

DNA nanocompartment is a typical DNA-based machine whose function is dependent of molecular collective effect. Fundamental properties of the device have been addressed via electrochemical analysis, fluorescent microscopy, and atomic force microscopy. Interesting and novel phenomena emerged during the switching of the device. We have found that DNAs in this system exhibit a much steep melting transition compared to ones in bulk solution or conventional DNA array. To achieve an understanding to this discrepancy, we introduced DNA-DNA interaction potential to the conventional Ising-like Zimm-Bragg theory and Peyrard-Bishop model of DNA melting. To avoid unrealistic numerical calculation caused by modification of the Peyrard-Bishop nonlinear Hamiltonian with the DNA-DNA interaction, we established coarse-gained Monte Carlo recursion relations by elucidation of five components of energy change during melting transition. The result suggests that DNA-DNA interaction potential accounts for the observed steep transition.