Energy localization, via modulation instability, is addressed in a modified twistopening model of DNA with solvent interactions. The Fourier expansion method is used to reduce the complex roto-torsional equations of the system to a set of discrete coupled nonlinear Schrödinger equations, which are used to perform the analytical investigation of modulation instability. We find that the instability criterion is highly influenced by the solvent parameters. Direct numerical simulations, performed on the generic model, further confirm our analytical predictions, as solvent interactions bring about highly localized energy patterns. These patterns are also shown to be robust under thermal fluctuations.