The dynamics and stability of four DNA duplexes are studied by means of molecular dynamics simulations. The four molecules studied are combinations of 4, 15 bases long, single-stranded oligomers, F1, F2, F3, and F4. The sequence of these single strand oligomers are chosen such that F1-F2 and F3-F4 form parallel (ps) DNA double helices, whereas F1-F4 and F2-F3 form antiparallel-stranded (aps) DNA double helices. Simulations were done at low (100 K) and room (300 K) temperatures. At low temperatures the dynamics are quasi-harmonic and the analysis of the trajectories gives good estimates of the low frequency vibrational modes and density of states. These are used to estimate the linear (harmonic) contribution of local fluctuations to the configurational entropy of the systems. Estimates of the differences in enthalpy between ps and aps duplexes show that aps double helices are more stable than the corresponding ps duplexes, in agreement with experiments. At higher temperatures, the distribution of the fluctuations around the average structures are multimodal and estimates of the configurational entropy cannot be obtained. The multi-basin, nonlinear character of the dynamics at 300 K is established using a novel method which extracts large amplitude nonlinear motions from the molecular dynamics trajectories. Our analysis shows that both ps DNA exhibit much larger fluctuations than the two aps DNA. The large fluctuations of ps DNA are explained in terms of correlated transitions in the beta, epsilon, and zeta backbone dihedral angles.