Binding of cationic and neutral phenanthridine intercalators to a DNA oligomer is controlled by dispersion energy: quantum chemical calculations and molecular mechanics simulations.

Correlated ab initio as well as semiempirical quantum chemical calculations and molecular dynamics simulations were used to study the intercalation of cationic ethidium, cationic 5-ethyl-6-phenylphenanthridinium and uncharged 3,8-diamino-6-phenylphenanthridine to DNA. The stabilization energy of the cationic intercalators is considerably larger than that of the uncharged one. The dominant energy contribution with all intercalators is represented by dispersion energy. In the case of the cationic intercalators, the electrostatic and charge-transfer terms are also important. The DeltaG of ethidium intercalation to DNA was estimated at -4.5 kcal mol(-1) and this value agrees well with the experimental result. Of six contributions to the final free energy, the interaction energy value is crucial. The intercalation process is governed by the non-covalent stacking (including charge-transfer) interaction while the hydrogen bonding between the ethidium amino groups and the DNA backbone is less important. This is confirmed by the evaluation of the interaction energy as well as by the calculation of the free energy change. The intercalation affects the macroscopic properties of DNA in terms of its flexibility. This explains the easier entry of another intercalator molecule in the vicinity of an existing intercalation site.

Correlated ab initio study of nucleic acid bases and their tautomers in the gas phase, in a microhydrated environment and in aqueous solution. Part 4. Uracil and thymine.

Altogether 13 keto and enol tautomers of uracil and 13 keto and enol tautomers of thymine were studied theoretically in the gas phase, in a microhydrated environment (1 and 2 water molecules) and in a water environment. Bulk water was described using the thermodynamic integration method, Conductor-like polarizable continuum model (C-PCM, COSMO) and hybrid model (C-PCM + 1-2 explicit water molecules). The structures of various tautomers were determined at the RI-MP2 level using the TZVPP basis set while relative energies were determined at the CCSD(T) level. The relative free energies at 298 K were based on the relative energies mentioned above and zero-point vibration energies, and temperature dependent enthalpy terms and entropies evaluated at the MP2/6-31G** level. The effect of bulk solvent on the relative stability of uracil and thymine tautomers was studied using molecular dynamics free energy calculations by means of the thermodynamic integration method and self-consistent reaction field. Despite the completely different nature of these methods they provide comparable solvation free energies. Besides theoretical investigation, experimental detection of uracil and thymine tautomers was performed by means of steady-state fluorescence. We conclude that it is impossible to utilize the method used by Suwaiyan and Morsy (M. A. Morsy, A. M. Al-Somali and A. Suwaiyan, J. Phys. Chem. B, 1999, 103(50), 11205) for tautomer detection, even if a very sensitive fluorimeter is used. Theoretical relative energies and free energies for isolated uracil and thymine tautomers support the existence of the canonical form only. The microhydrated environment and bulk solvent stabilize enol forms more than the canonical keto one, but gas phase destabilization of these enol forms is too high. Population of rare enol forms of uracil and thymine in bulk water will thus be very low and canonical structure will also be dominant in this phase.

New parameterization of the Cornell et al. empirical force field covering amino group nonplanarity in nucleic acid bases.

The Cornell et al. empirical potential (Cornell et al., J Am Chem Soc 1995, 117, 5197) was modified by the introduction of nonplanarity of the amino group in guanine, adenine, and cytosine. Reparameterization was performed for 12 bond parameters of the amino group (three valence angles (C-N-H(1), C-N-H(2), and H(1)-N-H(2)) and the improper dihedral angle at nitrogen, four dihedral angles (X-C-N-H(1), Y-C-N-H(2), X-X-C-N, Y-Y-C-N), three valence angles (X-C-N, Y-C-N, X-C-Y), and the improper dihedral angle at the adjacent carbon), and was based on correlated ab initio potential energy surfaces. Calculations were performed using the resolution of identity MP2 (RIMP2) method with SVP (3s2p1d/2s1p), TZVP (5s3p1d/3s1p), TZVPP (5s3p2d1f/3s2p1d), and augTZVPP (6s4p3d2f/4s3p2d) basis sets. Results obtained on the latter two levels are practically identical with the literature reference data [MP2/6-311++G(2df,p)] and the RIMP2/augTZVPP data can be considered a new reference set. The potential energy surface used for reparameterization was evaluated at the RIMP2/SVP level because respective geometry as well as energy data were close to the reference ones, and the computational time was very favorable. The modified potential was tested for 19 H-bonded and eight stacked nucleic acid base pairs. Interaction energies as well as geometries were described by the modified potential considerably better than by the original one and significant improvement resulted in the description of the nonplanar H-bonded and stacked complexes.

Correlated ab initio study of nucleic acid bases and their tautomers in the gas phase, in a microhydrated environment and in aqueous solution. guanine: surprising stabilization of rare tautomers in aqueous solution.

Altogether eight keto and enol tautomers of guanine were studied theoretically in the gas phase, in a microhydrated environment (1 and 2 water molecules) and in bulk water. The structures of isolated, as well as mono- and dihydrated tautomers were determined by means of the RI-MP2 method using the extended TZVPP (5s3p2d1f/3s2p1d) basis set. The relative energies of isolated tautomers included the correction to higher correlation energy terms evaluated at the CCSD(T)/aug-cc-pVDZ level. The relative enthalpies at 0 K and relative free energies at 298 K were based on the above-mentioned relative energies and zero-point vibration energies, temperature-dependent enthalpy terms and entropies evaluated at the MP2/6-31G level. The keto form having hydrogen atom at N7 is the global minimum while the canonical form having hydrogen atom at N9 represents the first local minimum at all theoretical levels in vacuo and in the presence of 1 and 2 water molecules. All three unusual rare tautomers having hydrogens at N3 and N7, at N3 and N9, and also at N9 and N7 are systematically considerably less stable and can be hardly detected in the gas phase. The theoretical predictions fully agree with existing theoretical as well as experimental results. The effect of bulk solvent on the relative stability of guanine tautomers was studied by self-consistent reaction field and molecular dynamics free energy calculations using the thermodynamic integration method. Bulk solvent, surprisingly, strongly favored these three rare tautomers over all remaining low-energy tautomers and probably only these forms can exist in water phase. The global minimum (tautomer with hydrogens at N3 and N7) is by 13 kcal/mol more stable than the canonical form (3rd local minimum). Addition of one or two water molecules does not change the relative stability order of isolated guanine tautomers but the respective trend clearly supports the surprising stabilization of three rare forms.

Molecular dynamics simulations and thermodynamics analysis of DNA-drug complexes. Minor groove binding between 4',6-diamidino-2-phenylindole and DNA duplexes in solution.

An extended set of nanosecond-scale molecular dynamics simulations of DNA duplex sequences in explicit solvent interacting with the minor groove binding drug 4',6-diamidino-2-phenylindole (DAPI) are investigated for four different and sequence specific binding modes. Force fields for DAPI have been parametrized to properly reflect its internal nonplanarity. Sequences investigated include the binding modes observed experimentally, that is, AATT in d(CGCGAATTCGCG)(2) and ATTG in d(GGCCAATTGG)(2) and alternative shifted binding modes ATTC and AATT, respectively. In each case, stable MD simulations are obtained, well reproducing specific hydration patterns seen in the experiments. In contrast to the 2.4 A d(CGCGAATTCGCG)(2) crystal structure, the DAPI is nonplanar, consistent with its gas-phase geometry and the higher resolution crystal structure. The simulations also suggest that the DAPI molecule is able to adopt different conformational substates accompanied by specific hydration patterns that include long-residing waters. The MM_PBSA technology for estimating relative free energies was utilized. The most consistent free energy results were obtained with an approach that uses a single trajectory of the DNA-DAPI complex to estimate all free energy terms. It is demonstrated that explicit inclusion of a subset of bound water molecules shifts the calculated relative binding free energies in favor of both crystallographically observed binding modes, underlining the importance of structured hydration.

Intercalators. 1. Nature of stacking interactions between intercalators (ethidium, daunomycin, ellipticine, and 4',6-diaminide-2-phenylindole) and DNA base pairs. Ab initio quantum chemical, density functional theory, and empirical potential study.

Properties of isolated intercalators (ethidium (E), daunomycin (D), ellipticine (EL), and 4,6'-diaminide-2-phenylindole (DAPI)) and their stacking interactions with adenine...thymine (AT) and guanine...cytosine (GC) nucleic acid base pairs were investigated by means of a nonempirical correlated ab initio method. All intercalators exhibit large charge delocalization, and none of them (including the DAPI dication) exhibits a site with dominant charge. All intercalators have large polarizability and are good electron acceptors, while base pairs are good electron donors. MP2/6-31G*(0.25) stabilization energies of intercalator...base pair complexes are large (E...AT, 22.4 kcal/mol; D...GC, 17.8 kcal/mol; EL...GC, 18.2 kcal/mol; DAPI...GC, 21.1 kcal/mol) and are well reproduced by modified AMBER potential (van der Waals radii of intercalator atoms are enlarged and their energy depths are increased). Standard AMBER potential underestimates binding, especially for DAPI-containing complexes. Because the DAPI dication is the best electron acceptor (among all intercalators studied), this difference is explained by the importance of the charge-transfer term, which is not included in the AMBER potential. For the neutral EL molecule, the standard AMBER force field provides correct results. The Hartree-Fock and DFT/B3LYP methods, not covering the dispersion energy, fail completely to reveal any energy minimum at the potential energy curve of the E...AT complex, and these methods thus cannot be recommended for a study of intercalation process. On the other hand, an approximate version of the DFT method, which was extended to cover London dispersion energy, yields for all complexes very good stabilization energies that are well comparable with referenced ab initio data. Besides the vertical dependence of the interaction, an energy twist dependence of the interaction energy was also investigated by a reference correlated ab initio method and empirical potentials. It is concluded that, despite the cationic (E +1, D +1, DAPI +2) or polar (EL) character of the intercalators investigated, it is the dispersion energy which predominantly contributes to the stability of intercalator...base pair complexes. Any procedure which does not cover dispersion energy is thus not suitable for studying the process of intercalation.