Temperature dependence for the rate of hole transfer in DNA: Nonadiabatic regime

The positive charge transfer in DNA is investigated, using the first principle treatment of the electron-vibrational interaction. We show that rearrangements of atoms belonging to base pairs induced by charge transfer are essentially quantum mechanical in nature. Particularly at room temperature, around half of the rearrangements occur via quantum tunneling, while the other half takes place via thermally activated transitions. This effect reduces activation energies for charge transfer between both AT and GC pairs by a factor of two compared to their classical values. These behaviors are described within small polaron theory for the non-adiabatic charge transfer and compared to the experimental data and previous theoretical studies.     

Charge hopping in DNA.

The efficiency of charge migration through stacked Watson-Crick base pairs is analyzed for coherent hole motion interrupted by localization on guanine (G) bases. Our analysis rests on recent experiments, which demonstrate the competition of hole hopping transitions between nearest neighbor G bases and a chemical reaction of the cation G(+) with water. In addition, it has been assumed that the presence of units with several adjacent stacked G bases on the same strand leads to the additional vibronic relaxation process (G(+)G...G) --> (GG...G)(+). The latter may also compete with the hole transfer from (G(+)G...G) to a single G site, depending on the relative positions of energy levels for G(+) and (G(+)G...G). A hopping model is proposed to take the competition of these three rate steps into account. It is shown that the model includes two important limits. One corresponds to the situation where the charge relaxation inside a multiple guanine unit is faster than hopping. In this case hopping is terminated by several adjacent G bases located on the same strand, as has been observed for the GGG triple. In the opposite, slow relaxation limit the GG...G unit allows a hole to migrate further in accord with experiments on strand cleavage exploiting GG pairs. We demonstrate that for base pair sequences with only the GGG triple, the fast relaxation limit of our model yields practically the same sequence- and distance dependencies as measurements, without invoking adjustable parameters. For sequences with a certain number of repeating adenine:thymine pairs between neighboring G bases, our analysis predicts that the hole transfer efficiency varies in inverse proportion to the sequence length for short sequences, with change to slow exponential decay for longer sequences. Calculations performed within the slow relaxation limit enable us to specify parameters that provide a reasonable fit of our numerical results to the hole migration efficiency deduced from experiments with sequences containing GG pairs. The relation of the results obtained to other theoretical and experimental studies of charge transfer in DNA is discussed. We propose experiments to gain a deeper insight into complicated kinetics of charge-transfer hopping in DNA.     

Electronic coupling mediated by stacked [Thymine-Hg-Thymine] base pairs

Very recently it has been shown that stable metal-mediated base pairs [Thymine-Hg-Thymine] can form in DNA. To estimate the effect of such pairs on the efficiency of charge transfer through DNA, we carry out quantum mechanical calculations of double-stranded pi-stacks GXG, GXXG, and GXXXG, where X = [Thymine-Hg-Thymine] and stacks GT(n)G of canonical base pairs. The charge-transfer efficiency in short duplexes GXG and GTG is found to be similar. However, the donor-acceptor coupling in GXXG and GXXXG is stronger by a factor of 2.5-3.0 than that in GT(n)G (n = 2 and 3), respectively. It is shown that the valence orbitals of Hg atoms do not essentially participate in mediating the electronic coupling for hole transfer; however, they may play an important role in excess electron transfer.     

Reorganization energy, activation energy, and mechanism of hole transfer process in DNA: a theoretical study

The density functional calculations with aug-cc-pVDZ basis sets on cationic guanine-cytosine (GC(+)) and adenine-thymine (AT(+)) base pairs suggest that the cationic charge is almost entirely localized on the G and A units with significant changes in the N-H and N...O distances around the H-bonded area. While the calculated intramolecular reorganization energy (lambda(v)) for a GC base pair (0.75 eV) is remarkably larger than that for an isolated G base (0.49 eV), for the AT base pairs these values (0.44 and 0.40 eV) are almost the same. The gas phase activation energies (E(a)) for GC(+)GC-->GCGC(+), AT(+)AT-->ATAT(+), and GC(+)AT-->GCAT(+) hole transfer processes are 0.19, 0.11, and 0.73 eV with rate constants of 1.69 x 10(11), 3.15 x 10(11), and 4.61(0.168) s(-1), respectively, at 298 K. An alternative mechanism of hole transfer has been proposed on the basis of energy barriers.     

Absolute rates of hole transfer in DNA

Absolute rates of hole transfer between guanine nucleobases separated by one or two A:T base pairs in stilbenedicarboxamide-linked DNA hairpins were obtained by improved kinetic analysis of experimental data. The charge-transfer rates in four different DNA sequences were calculated using a density-functional-based tight-binding model and a semiclassical superexchange model. Site energies and charge-transfer integrals were calculated directly as the diagonal and off-diagonal matrix elements of the Kohn-Sham Hamiltonian, respectively, for all possible combinations of nucleobases. Taking into account the Coulomb interaction between the negative charge on the stilbenedicarboxamide linker and the hole on the DNA strand as well as effects of base pair twisting, the relative order of the experimental rates for hole transfer in different hairpins could be reproduced by tight-binding calculations. To reproduce quantitatively the absolute values of the measured rate constants, the effect of the reorganization energy was taken into account within the semiclassical superexchange model for charge transfer. The experimental rates could be reproduced with reorganization energies near 1 eV. The quantum chemical data obtained were used to discuss charge carrier mobility and hole-transport equilibria in DNA.     

Protonic delocalization and quantum correlations in the H-bond dynamics of G–C and κ–π DNA base pairs

A new base pair (called κ–π) of Watson–Crick type, with a H -bond pattern different from that in A –T and G –C base pairs, has been recently synthesized and shown to be stable and incorporable into duplex DNA and RNA by polymerases. This new basepair, which contains three H -bonds, is compared with G –C , in the framework of modern dynamical theory of quantum nonlocality and quantum correlations. Connection with the traditional treatment of proton transfer in DNA base pairs, which uses the adiabatic approximation, is explicitly made. As a result, the dynamics of the H -bond pattern of G –C is shown to exhibit a specific quantum mechanical phase stability, which is clearly missing in the case of κ–π. This finding is discussed and illustrated, also in connection with recent quantum chemical calculations of proton transfers in DNA base pairs. Additionally, certain speculations concerning the “evolutionary advantage” of G –C with respect to κ–π are briefly considered. © 1993 John Wiley & Sons, Inc.     

Preparation of supramolecular chromophoric assemblies using a DNA duplex

Organization of supramolecular assemblies of chromophores with precisely-controlled orientation and sequence remains challenging. Nucleic acids with complementary base sequences spontaneously form double-helical structures. Therefore, covalent attachment of chromophores to DNA or RNA can be used to control assembly and orientation of chromophores. In this perspective, we first review our recent work on the assemblies of fluorophores (pyrene and perylene) by using natural base pairs. The interaction between dyes can be strictly controlled by means of cluster and interstrand wedge motifs. We then discuss novel artificial base pairs that can suppress the interaction between fluorophores and nucleobases. We incorporated a cyclohexane moiety into DNA, and showed that these artificial base pairs suppressed the electron-hole transfer between fluorophores and nucleobases and enhanced the quantum yields of fluorophores. These base pairs can potentially be used to accumulate fluorophores inside DNA duplexes without decreasing quantum yields.     

Elementary steps for charge transport in DNA: thermal activation vs. tunneling

Using stacks of Watson–Crick base pairs as an important example of multichromophoric molecular assemblies, we studied charge migration in DNA with special emphasis on the mechanism of hole hopping between neighboring guanines (G) connected by the adenine–thymine (AT) bridge. The tight-binding model proposed for this elementary step shows that for short AT bridges, hole transfer between two G bases proceeds via quantum mechanical tunneling. By contrast, hopping over long bridges requires thermal activation. The condition for crossover between tunneling and thermal activation near room temperature is specified and applies to the analysis of experimental data. We show that thermal activation dominates, if the bridge between two G bases contains more than three AT pairs. Our theoretical findings predict that the replacement of AT base pairs by GC pairs increases the efficiency of hole transport only in the case of short base pair sequences. For long sequences, however, the opposite effect is expected.     

Quasiclassical trajectory study of the collision-induced dissociation dynamics of Ar + CH3SH+ using an ab initio interpolated potential energy surface

Classical trajectory calculations have been performed to investigate the collision-induced dissociation (CID) of the CH(3)SH(+) cation with Ar atoms. A new intramolecular potential energy surface for the CH(3)SH(+) cation is evaluated by interpolation of 3000 ab initio data points calculated at the MP2/6-311G(d,p) level of theory. The new potential energy surface includes seven accessible dissociation channels of the cation. The present QCT calculations show that migration of hydrogen atoms, leading to the rearrangement CH(3)SH(+) <--> CH(2)SH(2)(+), is significant at the collision energies considered (6.5-34.7 eV) and that the formation of CH(3)(+), CH(3)S(+), and CH(2)(+) cations takes place primarily by a "shattering" mechanism in which the products are formed just after the collision. The theoretical product abundances are found to be in qualitative agreement with the experimental data. However, at high collision energies, the calculated total cross sections for the formation of CH(3)(+) and CH(2)SH(+) cations are noticeably larger than the experimental determinations. Several features of the dynamics of the CID processes are discussed.     

H-coupled electron transfer in alkane C-h activations with halogen electrophiles

The mechanisms for the reactions of isobutane and adamantane with polyhalogen electrophiles (HHal(2)(+), Hal(3)(+), Hal(5)(+), and Hal(7)(+), Hal = Cl, Br, or I) were studied computationally at the MP2 and B3LYP levels of theory with the 6-31G (C, H, Cl, Br) and 3-21G (I) basis sets, as well as experimentally for adamantane halogenations in Br(2), Br(2)/HBr, and I(+)Cl(-)/CCl(4). The transition structures for the activation step display almost linear C...H...Hal interactions and are characterized by significant charge transfer to the electrophile; the hydrocarbon moieties resemble the respective radical cation structures. The regiospecificities for polar halogenations of the 3-degree C-H bonds of adamantane, the high experimental kinetic isotope effects (k(H)/k(D) = 3-4), the rate accelerations in the presence of Lewis and proton (HBr) acids, and the high kinetic orders for halogen (7.5 for Br(2)) can only be understood in terms of an H-coupled electron-transfer mechanism. The three centered-two electron (3c-2e) electrophilic mechanistic concept based on the attack of the electrophile on a C-H bond does not apply; electrophilic 3c-2e interactions dominate the C-H activations only with nonoxidizing electrophiles such as carbocations. This was shown by a comparative computational analysis of the electrophilic and H-coupled electron-transfer activation mechanisms for the isobutane reaction with an ambident electrophile, the allyl cation, at the above levels of theory.     

Fluorescence Quenching over Short Range in a Donor‐DNA‐Acceptor System

A new donor‐DNA‐acceptor system has been synthesized containing Nile red‐modified 2′‐deoxyuridine as charge donor and 6‐N,N‐dimethylaminopyrene‐modified 2′‐deoxyuridine as acceptor to investigate the charge transfer in DNA duplexes using fluorescence spectroscopy and time‐resolved femtosecond pump‐probe techniques. Fluorescence quenching experiments revealed that the quenching efficiency of Nile red depends on two components: 1) the presence of a charge acceptor and 2) the number of intervening CG and AT base pairs between donor and acceptor. Surprisingly, the quenching efficiency of two base pairs (73 % for CG and the same for AT) is higher than that for one base pair (68 % for CG and 37 % for AT), while at a separation of three base pairs less than 10 % quenching is observed. A comparison with the results of time‐resolved measurements revealed a correlation between quenching efficiency and the first ultrafast time constant suggesting that quenching proceeds via a charge transfer from the donor to the acceptor. All transients are satisfactorily described with two decays: a rapid charge transfer with 600 fs (～10¹² s^?1) that depends strongly and in a non‐linear fashion on the distance between donor and acceptor, and a slower time constant of a few picoseconds (～10¹¹ s^?1) with weak distance dependence. A third time constant on a nanosecond time scale represents the fluorescence lifetime of the donor molecule. According to these results and time‐dependent density functional theory (TDDFT) calculations a combination of single‐step superexchange and multistep hopping mechanisms can be proposed for this short‐range charge transfer. Furthermore, significantly less quenching efficiency and slower charge transfer rates at very short distances indicate that the direct interaction between donor and acceptor leads to a local structural distortion of DNA duplexes which may provide some uncertainty in identifying the charge transfer rates in short‐range systems.     

A comparative DFT study of some N-based aromatic ligand metal complexes as anticancer agents and analysis of their mode of interaction with DNA base pair

Metal complexed anticancer agents interact with DNA nucleobase pairs (AT and GC) through different types of binding mode such as intercalation, groove binding, covalent binding, etc. Minor and major groove binding mechanism of DNA base pair is the key factor for all kinds of anticancer agent; as metal complexes have a great affinity to bind with DNA nucleobase either through minor or major groove. Ligands in metal complexes also play a vital role during the interaction with DNA base pairs; these ligands directly interact with DNA through different interacting modes. Generally, anticancer agents with less sterically hindered N-based aromatic and planar ligands are the key component for DNA binding; as the structure of such ligands are quite compatible for following intercalation and groove binding mechanism. Since, the experimental investigation for drug-DNA nucleobase complexes are extremely complicated, therefore; quantum mechanical calculations might be very helpful for computing the actual interactions in drug-DNA complexes. Quantum mechanical approaches such as density functional theory (DFT) might be a very important and useful tool to investigate the actual mode of interaction of metal complexed antitumor agents with DNA nucleobase. Herein, we have taken some metal complexes with N-based aromatic ligands as antitumor agents to investigate the proper mode of interaction between drug-DNA complexes.     

Charge separation in acridine- and phenothiazine-modified DNA

The formation of the long-lived, charge-separated state in DNA upon visible light irradiation is of particular interest in molecular-scale optoelectronics, sensor design, and other areas of nanotechnology. However, the efficient generation of the charge-separated state is hampered by fast charge recombination within a contact ion pair, which limits the application of DNA for photoelectrochemical sensors and devices. In this study, a series of protonated 9-alkylamino-6-chloro-2-methoxyacridine (Acr+)- and phenothiazine (Ptz)-modified DNAs were synthesized for the further understanding of the mechanism of charge separation in DNA to generate a long-lived, charge-separated state with a high quantum yield (Phi). The Acr+ serves as a photosensitizer to produce a hole on guanine (G), and the G-C base pairs were used as a hole-transporting pathway to separate a hole from Acr* (the one-electron-reduced form of Acr+) to be trapped at Ptz. Since Acr+ oxides only G upon photoexcitation, the A-T base pair can be used as a spacer between Acr+ and the G-C base pair to avoid the formation of a contact ion pair. The charge injection dynamics was investigated by steady-state fluorescence spectra and fluorescence lifetime measurements, and the Phi and the lifetime of the charge-separated state produced upon photoirradiation were assessed by nanosecond laser flash photolysis of the Acr+- and Ptz-modified DNA. A long-lived, charge-separated state was successfully formed upon visible-light irradiation, and the Phi was the highest for the DNA having a single intervening A-T base pair between Acr+ and the G-C base pair. These results clearly demonstrated that the charge separation process in DNA can be refined by putting a redox-inactive intervening base pair as a spacer between a photosensitizer and the nucleobase to be oxidized to slow down the charge recombination rate.     

A multicoincidence study of fragmentation dynamics in collision of γ-aminobutyric acid with low-energy ions

Fragmentation of the γ-aminobutyric acid molecule (GABA, NH(2)(CH(2))(3)COOH) following collisions with slow O(6+) ions (v≈0.3?a.u.) was studied in the gas phase by a combined experimental and theoretical approach. In the experiments, a multicoincidence detection method was used to deduce the charge state of the GABA molecule before fragmentation. This is essential to unambiguously unravel the different fragmentation pathways. It was found that the molecular cations resulting from the collisions hardly survive the interaction and that the main dissociation channels correspond to formation of NH(2)CH(2)(+), HCNH(+), CH(2)CH(2)(+), and COOH(+) fragments. State-of-the-art quantum chemistry calculations allow different fragmentation mechanisms to be proposed from analysis of the relevant minima and transition states on the computed potential-energy surface. For example, the weak contribution at [M-18](+), where M is the mass of the parent ion, can be interpreted as resulting from H(2)O loss that follows molecular folding of the long carbon chain of the amino acid.     

A molecular dynamics calculations of hole transfer rates in DNA strands

A computational model, which includes both tunneling and thermal hopping mechanisms, has been applied to study the charge transfer in DNA (GC)n and (AT)n strands. The calculations revealed the crucial role played by the A or G NH2-group vibrations in the hole transfer in both types of strands. Charge-transfer rates in the two strands have been determined based on the molecular dynamics calculations. They are in good agreement with the available experimental data. The modeling approach used here may be employed in the theoretical study of the charge transfer in natural and artificial DNA strands containing AT and GC pairs.     

Electronic Structure and Properties of Trihalogen Cations X(3)(+) and XY(2)(+) (X, Y = F, Cl, Br, I)

Electronic structures, charge distributions, geometries, valence force constants, and vibrational frequencies of the homoatomic clusters F(3)(+), Cl(3)(+), Br(3)(+), and I(3)(+) and of the heteroatomic clusters ClF(2)(+), BrF(2)(+), IF(2)(+), BrCl(2)(+), ICl(2)(+), and IBr(2)(+) were determined. The self-consistent field approach extended by MP2-correlation energy or density-functional corrections was applied using various basis sets. It was found that d- and f-type polarization functions play a crucial role as in some other halogen compounds. The MP2 approach yields the most satisfactory results. The effect of the crystalline environment surrounding the Cl(3)(+), Br(3)(+), and I(3)(+) species is successfully simulated by a Madelung potential. Frequencies calculated in the crystal field are in reasonable agreement with the more reliable ones among the experimental results. Coupling force constants were determined. They are not consistent with some empirical rules. Bonding and charge distributions of the formally mixed-valence systems X(+)Y(2)(0) are discussed. X(+) behaves like a divalent chalcogen with high electronegativity. Each of the X(+)-Y bonds in XY(2)(+) is very similar to the bond in X-Y. We predict the experimentally unknown F(3)(+) to be stable in vacuum but not in the solid state. Structures and frequencies of XY(2)(+) species, which are as yet unknown, are also predicted.     

Conceptual Quantum Chemical Analysis of Bonding and Noncovalent Interactions in the Formation of Frustrated Lewis Pairs

The contributions of covalent and noncovalent interactions to the formation of classical adducts of bulky Lewis acids and bases and frustrated Lewis pairs (FLPs) were scrutinized by using various conceptual quantum chemical techniques. Significantly negative complexation energies were calculated for fourteen investigated Lewis pairs containing bases and acids with substituents of various sizes. A Ziegler–Rauk‐type energy decomposition analysis confirmed that two types of Lewis pairs can be distinguished on the basis of the nature of the primary interactions between reactants; dative‐bond formation and concomitant charge transfer from the Lewis base to the acid is the dominant and most stabilizing factor in the formation of Lewis acid–base adducts, whereas weak interactions are the main thermodynamic driving force (>50 %) for FLPs. Moreover, the ease and extent of structural deformation of the monomers appears to be a key component in the formation of the former type of Lewis pairs. A Natural Orbital for Chemical Valence (NOCV) analysis, which was used to visualize and quantify the charge transfer between the base and the acid, clearly showed the importance and lack of this type of interaction for adducts and FLPs, respectively. The Noncovalent Interaction (NCI) method revealed several kinds of weak interactions between the acid and base components, such as dispersion, π–π stacking, C?H ??? π interaction, weak hydrogen bonding, halogen bonding, and weak acid–base interactions, whereas the Quantum Theory of Atoms in Molecules (QTAIM) provided further conceptual insight into strong acid–base interactions.     

Some halogenated anticancer agents,role of halide (-X) and their interaction mechanism with AT/GC base pairs: A computational study

Recently, a new class of halogen-based active anticancer agents have widely been developed which shows effective binding with AT/GC base pairs of DNA nucleobases. Usually intercalation, groove binding and covalent binding mechanisms are the most common drug-DNA binding pathways; but, the groove binding mechanism plays a crucial role in the stability of such drug-DNA complexes. As anticancer agent-DNA nucleobase interactions are very difficult to investigate by using common experimental techniques; therefore, theoretical methods may be quite helpful to analyze the proper mode of interaction for such drug-DNA systems. Past literature reveals that, quantum mechanical (QM) density functional theory (DFT) method is one of the best known tool for analyzing the different binding modes of halogenated anticancer agents with DNA nucleobases. Moreover, the halogen-bonding interaction in any biological system is fundamentally understood by investigating the mechanism of donor-acceptor complex formation between donor halogens and acceptor atoms within a receptor; such study is very competent for exploring the favoured anticancer agent-DNA interaction. In this current work, our main objective is to explore the effect of some intercalating and groove binding halogen-based anticancer agents with DNA nucleobase using computational method.     

Cisplatin cytotoxicity: a theoretical study of induced mutations

We investigate possible mutations in the genetic code induced by cisplatin with an approach combining molecular dynamics (MD) and hybrid quantum mechanics/molecular mechanics (QM/MM) calculations. Specifically, the impact of platination on the natural tautomeric equilibrium in guanine-cytosine (GC) base pairs is assessed to disclose the possible role played by non-canonical forms in anti-tumour activity. To obtain valuable predictions, the main interactions present in a real DNA environment, namely hydration and stacking, are simultaneously taken into account. According to our results, the Pt-DNA adduct promotes a single proton transfer reaction in GC in the DNA sequence AG[combining low line]G[combining low line]C. Such rare tautomers might play an important role in the cisplatin biological activity since they meet the stability requirements necessary to promote a permanent mutation.