Recent advances in scientific computing enable the full atomistic simulation of DNA molecules. However, there exists length and time scale limitations in molecular dynamics (MD) simulation for large DNA molecules. In this work, a two-level homogenization of DNA molecules is proposed. A wavelet projection method is first introduced to form a coarse-grained DNA molecule represented with superatoms. The coarsened MD model offers a simplified molecular structure for the continuum description of DNA molecules. The coarsened DNA molecular structure is then homogenized into a three-dimensional beam with embedded molecular properties. The methods to determine the elasticity constants in the continuum model are also presented. The proposed continuum model is adopted for the study of mechanical behavior of DNA loop.