Supramolecular switches based on the guanine-cytosine (GC) Watson-Crick pair: effect of neutral and ionic substituents

We have theoretically analyzed Watson-Crick guanine-cytosine (GC) base pairs in which purine-C8 and/or pyrimidine-C6 positions carry a substituent X = NH(-), NH(2), NH(3) (+) (N series), O(-), OH, or OH(2) (+) (O series), using the generalized gradient approximation (GGA) of density functional theory at the BP86/TZ2P level. The purpose is to study the effects on structure and hydrogen-bond strength if X= H is substituted by an anionic, neutral, or cationic substituent. We found that replacing X = H by a neutral substituent has relatively small effects. Introducing a charged substituent, on the other hand, led to substantial and characteristic changes in hydrogen-bond lengths, strengths, and hydrogen-bonding mechanism. In general, introducing an anionic substituent reduces the hydrogen-bond-donating and increases the hydrogen-bond-accepting capabilities of a DNA base, and vice versa for a cationic substituent. Thus, along both the N and O series of substituents, the geometric shape and bond strength of our DNA base pair can be chemically switched between three states, thus yielding a chemically controlled supramolecular switch. Interestingly, the orbital-interaction component in some of these hydrogen bonds was found to contribute to more than 49 % of the attractive interactions and is thus virtually equal in magnitude to the electrostatic component, which provides the other (somewhat less than) 51 % of the attraction.     

Hydrogen bonding in mimics of Watson-Crick base pairs involving C-H proton donor and F proton acceptor groups: a theoretical study.

We have theoretically analyzed mimics of Watson-Crick adenine-thymine (AT) and guanine-cytosine (GC) base pairs in which N-H...O and N...H-N hydrogen bonds are replaced by N-H...F and N...H-C, respectively, by using the generalized gradient approximation of density functional theory at BP86/TZ2P. The general effect of the above substitutions is an elongation and weakening of the hydrogen bonds that hold together the base pairs. However, the precise effects depend on how many and, in particular, on which hydrogen bonds are substituted in AT and GC. Another purpose of this work is to clarify the relative importance of electrostatic attraction versus orbital interaction in the weak hydrogen bonds involved in the mimics, by using a quantitative bond-energy decomposition scheme. At variance with widespread believe, the orbital interaction component in these weak hydrogen bonds is found to contribute 34-42% of the attractive interactions and is thus of the same order of magnitude as the electrostatic component, which provides the remaining attraction component.     

Watson-crick base pairs with thiocarbonyl groups: How sulfur changes the hydrogen bonds in DNA

We have theoretically analyzed mimics of Watson-Crick AT and GC base pairs in which N-H···O hydrogen bonds are replaced by N-H···S, using the generalized gradient approximation (GGA) of density functional theory at BP86/TZ2P level. The general effect of the above substitutions is an elongation and a slight weakening of the hydrogen bonds that hold together the base pairs. However, the precise effects depend on how many, and in particular, on which hydrogen bonds AT and GC are substituted.. Another purpose of this work is to clarify the relative importance of electrostatic attraction versus orbital interaction in the hydrogen bonds involved in the mimics, using a quantitative bond energy decomposition scheme. At variance with widespread believe, the orbital interaction component in these hydrogen bonds is found to contribute more than 40% of the attractive interactions and is thus of the same order of magnitude as the electrostatic component, which provides the remaining attraction.      

Controlling Emissive Properties by Intramolecular Hydrogen Bonds: Alkyl and Aryl meso-Substituted Porphycenes

Porphycene, a porphyrin isomer, is an efficient fluorophore. However, four-fold meso substitution with alkyl groups decreases the fluorescence quantum yield by orders of magnitude. For aryl substituents, this effect is small. To explain this difference, we have synthesized and studied a mixed aryl-alkyl-substituted compound, 9,20-diphenyl-10,19-dimethylporphycene, as well as the 9,20-diphenyl and 9,20-dimethyl derivatives. Analysis of the structural, spectroscopic, and photophysical data of the six porphycenes, combined with quantum chemical calculations, shows a clear correlation between the strength of the intramolecular NH⋅⋅⋅N hydrogen bonds and the efficiency of the radiationless depopulation of the lowest-excited singlet state. This result led us to propose a model in which the delocalization of the inner protons in the cavity of the macrocycle is responsible for the nonradiative deactivation channel. The applicability of the model is confirmed by the literature data for other alkyl- or aryl-substituted porphycenes. The finding of a correlation between structural and emissive characteristics enables a rational design of porphycenes with desired photophysical properties.     

Substituent Effects on Hydrogen Bonding in Watson–Crick Base Pairs. A Theoretical Study

We have theoretically analyzed Watson–Crick AT and GC base pairs in which purine C8 and/or pyrimidine C6 positions carry a substituent X = H, F, Cl or Br, using the generalized gradient approximation (GGA) of density functional theory at BP86/TZ2P. The purpose is to study the effects on structure and hydrogen bond strength if X = H is substituted by a halogen atom. Furthermore, we wish to explore the relative importance of electrostatic attraction versus orbital interaction in the above multiply hydrogen-bonded systems, using a quantitative bond energy decomposition scheme. We find that replacing X = H by a halogen atom has relatively small yet characteristic effects on hydrogen bond lengths, strengths and bonding mechanism. In general, it reduces the hydrogen-bond-accepting- and increases the hydrogen-bond-donating capabilities of a DNA base. The orbital interaction component in these hydrogen bonds is found for all substituents (X = H, F, Cl, and Br) to contribute about 41% of the attractive interactions and is thus of the same order of magnitude as the electrostatic component, which provides the remaining 59% of the attraction.     

Kinetic Stabilization and Reactivity of π Single‐Bonded Species: Effect of the Alkoxy Group on the Lifetime of Singlet 2,2‐Dialkoxy‐1,3‐diphenyloctahydropentalene‐1,3‐diyls

Kinetic stabilization and reactivity of π single‐bonded species have been investigated in detail by generating a series of singlet 2,2‐dialkoxy‐1,3‐diphenyloctahydropentalene‐1,3‐diyls ( DR s). The lifetime at 293 K in benzene was found to increase when the carbon chain length of the alkoxy groups was increased; 292 ns ( DRb ; OR=OR′=OCH₃) <880 ns ( DRc ; OR=OR′=OC₂H₅) <1899 ns ( DRd ; OR=OR′=OC₃H₇) ≈2292 ns ( DRe ; OR=OR′=OC₆H₁₃) ≈2146 ns ( DRf ; OR=OR′=OC₁₀H₂₁). DRh (OR=OC₃H₇, OR′=OCH₃; 935 ns) with the mixed‐acetal moiety is a longer‐lived species than another diastereomer DRg (OR=OCH₃, OR′=OC₃H₇; 516 ns). Activation parameters determined for the first‐order decay process reveal that the enthalpy factor plays a crucial role in determining the energy barrier of the ring‐closing reaction, that is, from the π‐bonding to the σ‐bonding compounds. Computational studies using density functional theory provided more insight into the structures of the singlet species with π single‐bonded character and the transition states for the ring‐closing reaction, thereby clarifying the role of the alkoxy group on the lifetime and the stereoselectivity of the ring‐closing reaction.     

Characterization of the monovalent ion position and hydrogen-bond network in guanine quartets by DFT calculations of NMR parameters

Conformational stability of G-quartets found in telomeric DNA quadruplex structures requires the coordination of monovalent ions. Here, an extensive Hartree-Fock and density functional theory analysis of the energetically favored position of Li+, Na+, and K+ ions is presented. The calculations show that at quartet-quartet distances observed in DNA quadruplex structures (3.3 A), the Li+ and Na+ ions favor positions of 0.55 and 0.95 A outside the plane of the G-quartet, respectively. The larger K+ ion prefers a central position between successive G-quartets. The energy barrier separating the minima in the quartet-ion-quartet model are much smaller for the Li+ and Na+ ions compared with the K+ ion; this suggests that K+ ions will not move as freely through the central channel of the DNA quadruplex. Spin-spin coupling constants and isotropic chemical shifts in G-quartets extracted from crystal structures of K+- and Na+-coordinated DNA quadruplexes were calculated with B3LYP/6-311G(d). The results show that the sizes of the trans-hydrogen-bond couplings are influenced primarily by the hydrogen bond geometry and only slightly by the presence of the ion. The calculations show that the R(N2N7) distance of the N2-H2...N7 hydrogen bond is characterized by strong correlations to both the chemical shifts of the donor group atoms and the (h2)J(N2N7) couplings. In contrast, weaker correlations between the (h3)J(N1C6') couplings and single geometric factors related to the N1-H1...O6=C6 hydrogen bond are observed. As such, deriving geometric information on the hydrogen bond through the use of trans-hydrogen-bond couplings and chemical shifts is more complex for the N1-H1...O6=C6 hydrogen bond than for the N2-H2...N7 moiety. The computed trans-hydrogen-bond couplings are shown to correlate with the experimentally determined couplings. However, the experimental values do not show such strong geometric dependencies.     

Molecular Switching Behavior in Isosteric DNA Base Pairs

The structures and proton‐coupled behavior of adenine–thymine (A‐T) and a modified base pair containing a thymine isostere, adenine–difluorotoluene (A‐F), are studied in different solvents by dispersion‐corrected density functional theory. The stability of the canonical Watson–Crick base pair and the mismatched pair in various solvents with low and high dielectric constants is analyzed. It is demonstrated that A‐F base pairing is favored in solvents with low dielectric constant. The stabilization and conformational changes induced by protonation are also analyzed for the natural as well as the mismatched base pair. DNA sequences capable of changing their sequence conformation on protonation are used in the construction of pH‐based molecular switches. An acidic medium has a profound influence in stabilizing the isostere base pair. Such a large gain in stability on protonation leads to an interesting pH‐controlled molecular switch, which can be incorporated in a natural DNA tract.     

Self-assembled organic nanostructures: effect of substituents on the morphology.

Three perylene-3,4;9,10-tetracarboxydiimide (PTCDI) compounds with two dodecyloxy or thiododecyl chains attached at the bay positions of the perylene ring, PTCDIs 1-3, were fabricated into nanoassemblies by a solution injection method. The morphologies of these self-assembled nanostructures were determined by transmission electronic microscopy (TEM), scanning electronic microscopy (SEM), and atomic force microscopy (AFM). PTCDI compound 1, with two dodecyloxy groups, forms long, flexible nanowires with an aspect ratio of over 200, while analogue 3, with two thiododecyl groups, self-assembles into spherical particles. In line with these results, PTCDI 2, with one dodecyloxy group and one thiododecyl group, forms nanorods with an aspect ratio of around 20. Electronic absorption and fluorescence spectroscopy results reveal the formation of H-aggregates in the nanostructures of these PTCDI compounds owing to the pi-pi interaction between the substituted perylene molecules and also suggest a decreasing pi-pi interaction in the order 1>2>3, which corresponds well with the morphology of the corresponding nanoassemblies. On the basis of DFT calculations, the effect of different substituents at the bay positions of the perylene ring on the pi-pi interaction between substituted perylene molecules and the morphology of self-assembled nanostructures is rationalized by the differing degree of twisting of the conjugated perylene system caused by the different substituents and the different bending of the alkoxy and thioalkyl groups with respect to the plane of the naphthalene.     

Ionization-induced proton transfer in model DNA base pairs: a theoretical study of the radical ions of the 7-azaindole dimer.

Proton-transfer reactions of the radical anion and cation of the 7-Azaindole (7AI) dimer were investigated by means of density functional theory (DFT). The calculated results for the dimer anion and cation were very similar. Three equilibrium structures, which correspond to the non-proton-transferred (normal), the single-proton-transferred (SPT) and the double-proton-transferred (tautomeric) forms, were found. The transition states for proton-transfer reactions were also located. The calculations showed that the first proton-transfer reaction (normal-->SPT) is exothermic and almost barrier-free; therefore, it should occur spontaneously in the period of a vibration. In contrast, the second proton-transfer reaction (SPT-->tautomer) was found to be far less-probable in terms of reaction energy and barrier. Hence, it was concluded that both (7Al)2*- and (7Al)2*+ exist in the SPT form. The conclusion was further confirmed by the calculated electron vertical detachment energy (VDE) of the SPT form of (7Al)2*-, 1.33 eV, which is very close to the experimental measurement of 1.35 eV. The calculated VDEs of the normal and tautomeric (7Al)2*- forms were too small compared to the experimental value. The proton transfer process was found to be multidimensional in nature involving not only proton motion but also intermolecular rocking motion. In addition, IR spectra were calculated and reported. The spectra of the three structures showed very different features and, therefore, can be considered as fingerprints for future experimental identifications. The implications of these results to biology and spectroscopy are also briefly discussed.     

Telomere structure and stability: covalency in hydrogen bonds, not resonance assistance, causes cooperativity in guanine quartets

We show that the cooperative reinforcement between hydrogen bonds in guanine quartets is not caused by resonance-assisted hydrogen bonding (RAHB). This follows from extensive computational analyses of guanine quartets (G(4)) and xanthine quartets (X(4)) based on dispersion-corrected density functional theory (DFT-D). Our investigations cover the situation of quartets in the gas phase, in aqueous solution as well as in telomere-like stacks. A new mechanism for cooperativity between hydrogen bonds in guanine quartets emerges from our quantitative Kohn-Sham molecular orbital (MO) and corresponding energy decomposition analyses (EDA). Our analyses reveal that the intriguing cooperativity originates from the charge separation that goes with donor-acceptor orbital interactions in the σ-electron system, and not from the strengthening caused by resonance in the π-electron system. The cooperativity mechanism proposed here is argued to apply, beyond the present model systems, also to other hydrogen bonds that show cooperativity effects.     

Electrostatic calculation of the substituent effect: an efficient test on isolated molecules

The energy of a disubstituted molecule has often been approximated by simple electrostatic formulas that represent the substituents as poles or dipoles. Herein, we test this approach on a new model system that is more direct and more efficient than testing on acid-base properties. The energies of 27 1,4-derivatives of bicyclo[2.2.2]octane were calculated within the framework of the density functional theory at the B3LYP/6-311+G(d,p) level; interaction of the two substituents was evaluated in terms of isodesmic homodesmotic reactions. This interaction energy, checked previously on some experimental gas-phase acidities, was considered to be accurate and served as reference to test the electrostatic approximation. This approximation works well in the qualitative sense as far as the sign and the order of magnitude are concerned: beginning with the strongest interaction between two poles, a weaker interaction between pole and dipole, and the weakest between two dipoles. However, all the electrostatic calculations yield energies that are too small, particularly for weak interaction, and this fundamental defect is not remedied by some possible improvements. In particular, variation of the effective permittivity would require a physically impossible value less than unity. The explanation must lie in a more complex distribution of electron density than anticipated in the electrostatic model. It also follows that possible conclusions about the transmission of substituent effects "through space" have little validity.     

Influence of 2'-deoxy sugar moiety on excited-state protonation equilibrium of adenine and adenosine with acridine inside SDS micelles: a time-resolved study with quantum chemical calculations

The protonation dynamics of the DNA base adenine (Ade) and its nucleoside 2'-deoxyadenosine (d-Ade) are investigated by monitoring the deprotonation kinetics of an N-heterocyclic DNA intercalator, acridine (Acr), in the confined environment of sodium dodecyl sulfate (SDS) micelles. Protonation of acridine (AcrH(+)) occurs at the hydrophilic interface and this species remains in dynamic equilibrium with its deprotonated counterpart (Acr) inside the hydrophobic core of SDS micelles. Quenching of the fluorescence of AcrH(+)* at 478 nm is observed after addition of Ade and d-Ade with Stern-Volmer constant (K(SV)) 298 and 75 M(-1), respectively, with a concomitant increment in Acr* at 425 nm. Time-resolved fluorescence studies reveal quenching in the lifetime of AcrH(+)*. The relative amplitude of AcrH(+)* decreases from 0.97 to 0.51 and 0.97 to 0.89 with equimolar addition of Ade and d-Ade, respectively. These observations are explained by excited-state proton transfer (ESPT) from AcrH(+)* to the bases. The reduced K(SV) value and negligible change in the relative amplitudes of AcrH(+)* with d-Ade infer that ESPT is hindered substantially by the presence of a 2'-deoxy sugar unit. Transient time-resolved absorption spectra of Acr reflect that Ade reduces the absorbance of (3)AcrH(+)*; however, d-Ade keeps it unaltered for more than a time delay of 2 μs. The optimized geometries calculated by quantum chemical methods reflect deprotonation of AcrH(+)* with protonation at the N1 position of Ade, while it remains protonated with d-Ade. The hindered ESPT between AcrH(+)* and d-Ade singles out the significance of the 2'-deoxy sugar moiety in controlling the deprotonation kinetics.     

Influence of Solvent Polarity and Hydrogen Bonding on the Electronic Transition of Coumarin 120: A TDDFT Study

The characteristics of the electronic transition energy of Coumarin 120 (C120) and its H‐bonded complexes in various solvents have been examined by time‐dependent density functional theory (TDDFT) in combination with a polarizable continuum solvent model (PCM). Molecular structures of C120 and its H‐bonded complexes are optimized with the B3LYP method in PCM solution, and the dihedral angle H14? N13? C7? H15 is dependent on solvent polarity and the type of H‐bond. A linear correlation of the absorption maximum of C120 with the solvent polarity function is revealed with the PCM model for all solvents except DMSO. The experimental absorption maximum of C120 in nine solvents is well described by a PCM–TDDFT scheme augmented with explicit inclusion of a few H‐bonded solvent molecules, and quantitative agreement between our calculated results and experimental measurements is obtained with an average error of less than 2 nm. H‐bonding at three different sites shifts the absorption wavelength of C120 either to the blue or to the red, that is, a significant role is played by solvent molecules in the first solvation shell in determining the electronic transition energy of C120. The dependence on the H‐bonding site and solvent polarity is examined by using the Kamlet–Taft equation for solvatochromism.     

Theoretical Study of the Stability of the DNA Duplexes Modified by a Series of Hydrophobic Base Analogues

The geometries of a 13 mer of a DNA double helix (5′‐GCGTAC A CATGCG‐3′) were determined by molecular dynamics simulations using a Cornell et al. empirical force field. The bases in the central base pair (shown in bold) were replaced (one or both) by a series of hydrophobic base analogues (phenyl, biphenyl, phenylnaphathalene, phenylanthracene and phenylphenanthrene). Due to the large fluctuations of the systems, an average geometry could not be determined. The interaction energies of the Model A, which consisted of three central steps of a duplex without a sugar phosphate backbone, taken from molecular dynamics simulations (geometry sampled every 1 ps), were calculated by the self‐consistent charge density functional based tight‐binding (SCC‐DFTB‐D) method and were subsequently averaged. The higher the stability of the systems the higher the aromaticity of the base analogues. To estimate the desolvation energy of the duplex, the COSMO continuum solvent model was used and the calculations were provided on a larger model, Model B (the three central steps of the duplex with a sugar phosphate backbone neutralised by H atoms), taken from molecular dynamics simulations (geometry sampled every 200 ps) and subsequently averaged. The selectivity of the base analogue pairs was ascertained (Model B) by including the desolvation energy and the interaction energy of both strands, as determined by the SCC‐DFTB‐D method. The highest selectivity was found for a phenylphenanthrene. Replacing the nucleic acid bases with a base analogue leads to structural changes of the central pair. Only with the smallest base analogues (phenyl) does the central base pair stay planar. When passing to larger base analogues the central base pair is usually stacked.