Model for melting transition of twisted DNA in a thermal bath

We investigated the melting transition of deoxyribonucleic acid (DNA) embedded in a Langevin fluctuation–dissipation thermal bath. Torsional effects were taken into consideration by introducing a twist angle φ between neighboring base pairs stacked along the molecule backbone. We use the Barbi–Cocco–Payrard model to numerically study the impact of the twist angle on the melting temperature, considering four different sequences composed of 69 base pairs. According to the outcomes of our simulation, for all heterogeneous sequences, an increase in twist angle leads to a linear rise in melting temperature with a positive slope. For angles greater than the so-called equilibrium angle, the DNA chain becomes very rigid against opening and accordingly high temperatures are required to initiate the melting process. We also investigate the opening probability of bubbles, the bubble lifetime profiles and bubble length along the different DNA sequences. Graphical abstract