Abstract The accurate ab initio quantum chemical (QM) method multiconfigurational second-order perturbation (CASSPT2)/complete active space self-consistent field (CASSCF) has been used in conjunction with molecular mechanics (MM) procedures to compute molecular properties and photoinduced reactivity of DNA/RNA nucleobases (NABs) in isolation and within a realistic environment, in which the double helix strand, the aqueous media, and the external counterions are included. It is illustrated that the use of an MM model is helpful both to account for short- and long-range effects of the system surrounding the QM molecular core and to provide the proper structural constraints that allow more accurate QM geometry determinations.