Importance of model size in quantum mechanical studies of DNA intercalation

Fast and Accurate Quantum Mechanical Modeling of Large Molecular Systems Using Small Basis Set Hartree-Fock Methods Corrected with Atom-Centered Potentials

There has been significant interest in developing fast and accurate quantum mechanical methods for modeling large molecular systems. In this work, by utilizing a machine learning regression technique, we have developed new low-cost quantum mechanical approaches to model large molecular systems. The developed approaches rely on using one-electron Gaussian-type functions called atom-centered potentials (ACPs) to correct for the basis set incompleteness and the lack of correlation effects in the underlying minimal or small basis set Hartree-Fock (HF) methods. In particular, ACPs are proposed for ten elements common in organic and bioorganic chemistry (H, B, C, N, O, F, Si, P, S, and Cl) and four different base methods: two minimal basis sets (MINIs and MINIX) plus a double-ζ basis set (6-31G*) in combination with dispersion-corrected HF (HF-D3/MINIs, HF-D3/MINIX, HF-D3/6-31G*) and the HF-3c method. The new ACPs are trained on a very large set (73 832 data points) of noncovalent properties (interaction and conformational energies) and validated additionally on a set of 32 048 data points. All reference data are of complete basis set coupled-cluster quality, mostly CCSD(T)/CBS. The proposed ACP-corrected methods are shown to give errors in the tenths of a kcal/mol range for noncovalent interaction energies and up to 2 kcal/mol for molecular conformational energies. More importantly, the average errors are similar in the training and validation sets, confirming the robustness and applicability of these methods outside the boundaries of the training set. In addition, the performance of the new ACP-corrected methods is similar to complete basis set density functional theory (DFT) but at a cost that is orders of magnitude lower, and the proposed ACPs can be used in any computational chemistry program that supports effective-core potentials without modification. It is also shown that ACPs improve the description of covalent and noncovalent bond geometries of the underlying methods and that the improvement brought about by the application of the ACPs is directly related to the number of atoms to which they are applied, allowing the treatment of systems containing some atoms for which ACPs are not available. Overall, the ACP-corrected methods proposed in this work constitute an alternative accurate, economical, and reliable quantum mechanical approach to describe the geometries, interaction energies, and conformational energies of systems with hundreds to thousands of atoms.

Theoretical Investigation of the 4,5-Dibromorodamine Methyl Ester (TH9402) Photosensitizer Used in Photodynamic Therapy: Photophysics, Reactions in the Excited State, and Interactions with DNA

Photosensitizer (PS) molecules play a critical role in photodynamic therapy of cancer and the understanding of the molecular mechanism involved in the photophysics of these compounds, and their reactions in the excited state are, therefore, of great interest for the development of this technique. In this article, the photophysics of the cationic PS 4,5-dibromorodamine methyl ester (TH9402), its electron- and energy-transfer reactions in the excited triplet state, with molecular oxygen, nitric oxide, guanosine-5'-monophosphate (GMP), and guanine, and the interaction with DNA were evaluated. Time-dependent density functional theory calculations at the TPSSh/Def2-TZVP//B3LYP/Def2-TZVP level of theory in water solution reveals that the PS has a bright S₁ state 2.33 eV above the ground state that produces a fluorescent rate constant of 5.40 × 10⁷ s^-1, calculated using Fermi's golden rule within a path integral formalism. Once excited to the bright state, the main intersystem crossing (ISC) channel involves the coupling with the T₂ state just below S₁ (S₁ → T₂ → T₁) with an overall ISC rate constant of 10.1 × 10⁷ s^-1, in good agreement with the experimental data. Excited-state reaction thermodynamics, computed at the M06-2X/Def2-TZVP//B3LYP/Def2-TZVP level of theory in water, showed that from all the excited-state electron-transfer reactions studied, only the transfer from GMP to the PS is thermodynamically favorable, independent of the protonation state of guanosine, which indicates a possible DNA photo-oxidation mechanism for the PS. Triplet-triplet energy-transfer reactions from TH9402 to molecular oxygen, producing reactive singlet oxygen, and to the deprotonated guanosine, producing ³GMP^2-, are also thermodynamically favorable, with ΔG = -2.0 and -24.0 kcal//mol, respectively. However, the energy transfer to the monoprotonated guanosine is not favorable, (ΔG = 36.1), suggesting that in the DNA double-strand environment, this energy-transfer process may not be observed. The results show that the PS can act through electron transfer and triplet-triplet energy-transfer reactions involved in mechanism types I and II in photodynamic therapy. Interactions of TH9402 with the d(AGACGTCT)₂ octanucleotide revealed that the PS can intercalate between the d(GpC)-d(CpG) base pairs in three different orientations and, upon intercalation, the π → π* transition of the PS shows a bathochromic shift up to 90 nm and up to 60% decrease in intensity. Interactions through groove binding showed a smaller bathochromic shift of 52.2 nm and a 56% decrease in intensity of the main transition band.