Environmental effects on the enhancement in natural and damaged DNA nucleobase acidity because of discrete hydrogen-bonding interactions

The present study uses density functional theory to carefully consider the effects of the environment on the enhancement in (natural and damaged) DNA nucleobase acidities because of multiple hydrogen-bonding interactions. Although interactions with one small molecule can increase the acidity of the nucleobases by up to 60 kJ mol-1 in the gas phase, the maximum increase in enzymatic-like environments is expected to be approximately 40 kJ mol-1, which reduces to approximately 30 kJ mol-1 in water. Furthermore, the calculated (simultaneous) effects of two, three, or four molecules are increasingly less than the sum of the individual (additive) effects with an increase in the number and acidity of the small molecules bound or the dielectric constant of the solvent. Regardless of these trends, our calculations reveal that additional hydrogen-bonding interactions will have a significant effect on nucleobase acidity in a variety of environments, where the exact magnitude of the effect depends on the properties of the small molecule bound, the nucleobase binding site, and the solvent. The maximum increase in nucleobase acidity because of interactions with up to four small molecules is approximately 80 kJ mol-1 in enzymatic-like environments (or 65 kJ mol-1 in water). These results suggest that hydrogen-bonding interactions likely play an important role in many biological processes by changing the physical and chemical properties of the nucleobases.     

Direct on-strip analysis of size- and time-resolved aerosol impactor samples using laser induced fluorescence spectra excited at 263 and 351 nm

We report a novel atmospheric aerosol characterization technique, in which dual wavelength UV laser induced fluorescence (LIF) spectrometry marries an eight-stage rotating drum impactor (RDI), namely UV-LIF-RDI, to achieve size- and time-resolved analysis of aerosol particles on-strip. The UV-LIF-RDI technique measured LIF spectra via direct laser beam illumination onto the particles that were impacted on a RDI strip with a spatial resolution of 1.2 mm, equivalent to an averaged time resolution in the aerosol sampling of 3.6 h. Excited by a 263 nm or 351 nm laser, more than 2000 LIF spectra within a 3-week aerosol collection time period were obtained from the eight individual RDI strips that collected particles in eight different sizes ranging from 0.09 to 10 μm in Djibouti. Based on the known fluorescence database from atmospheric aerosols in the US, the LIF spectra obtained from the Djibouti aerosol samples were found to be dominated by fluorescence clusters 2, 5, and 8 (peaked at 330, 370, and 475 nm) when excited at 263 nm and by fluorescence clusters 1, 2, 5, and 6 (peaked at 390 and 460 nm) when excited at 351 nm. Size- and time-dependent variations of the fluorescence spectra revealed some size and time evolution behavior of organic and biological aerosols from the atmosphere in Djibouti. Moreover, this analytical technique could locate the possible sources and chemical compositions contributing to these fluorescence clusters. Advantages, limitations, and future developments of this new aerosol analysis technique are also discussed.     

The hydrogen bonding properties of cytosine: a computational study of cytosine complexed with hydrogen fluoride, water, and ammonia

Density functional theory is used to study the hydrogen bonding pattern in cytosine, which does not contain alternating proton donor and acceptor sites and therefore is unique compared with the other pyrimidines. Complexes between various small molecules (HF, H(2)O, and NH(3)) and four main binding sites in (neutral and (N1) anionic) cytosine are considered. Two complexes (O2(N1) and N3(N4)) involve neighboring cytosine proton acceptor and donor sites, which leads to cooperative interactions and bidendate hydrogen bonds. The third (less stable) complex (N4) involves a single cytosine donor. The final (O2-N3) complex involves two cytosine proton acceptors, which leads to an anticooperative hydrogen bonding pattern for H(2)O and NH(3). On the neutral surface, the anticooperative O2-N3 complex is less stable than those involving bidentate hydrogen bonds, and the H(2)O complex cannot be characterized when diffuse functions are included in the (6-31G(d,p)) basis set. On the contrary, the anionic O2-N3 structure is the most stable complex, while the HF and H(2)O N3(N4) complexes cannot be characterized with diffuse functions. B3LYP and MP2 potential energy surface scans are used to consider the relationship between the water N3(N4) and O2-N3 complexes. These calculations reveal that diffuse functions reduce the conversion barrier between the two complexes on both the neutral and anionic surfaces, where the reduction leads to a (O2-N3) energy plateau on the neutral surface and complete (N3(N4)) complex destabilization on the anionic surface. From these complexes, the effects of hydrogen bonds on the (N1) acidity of cytosine are determined, and it is found that the trends in the effects of hydrogen bonds on the (N1) acidity are similar for all pyrimidines.     

yDNA versus xDNA pyrimidine nucleobases: computational evidence for dependence of duplex stability on spacer location

The structural and binding properties of the natural and x- and y-pyrimidines were compared using computational methods. Our calculations show that although the x-pyrimidines favor different orientations about the glycosidic bond compared to the natural pyrimidines, which could have implications for the formation and resulting stability of xDNA duplexes and jeopardize the selectivity of expanded nucleobases, y-pyrimidines have rotational profiles more similar to the natural bases. Increasing the pyrimidine size using a benzene spacer leads to relatively minor changes in the hydrogen-bond strength of isolated Watson-Crick base pairs. However, differences in the anomeric carbon distances in pairs composed of x- or y-pyrimidines suggest yDNA may yield a more optimal expanded structure. By stacking two monomers via their centers of mass, we find that the expanded nucleobases stack much stronger than the natural bases. Additionally, although replacing xT by yT changes the stacking energy by less than 5 kJ mol (-1), replacing xC by yC significantly strengthens complexes with the natural nucleobases (by up to 30%). Calculations on larger duplex models composed of four nucleobases reveal that x- and y-pyrimidines can increase duplex stability of natural helices by strengthening both the intra and interstrand stacking interactions. Furthermore, when the total stability (sum of all hydrogen-bonding and (intrastrand and interstrand) stacking interactions) of the larger models is considered, y-pyrimidines lead to more stable complexes than x-pyrimidines for all but three duplex sequences. Thus, through analysis of a variety of properties, our calculations suggest that the location of the benzene spacer affects the properties of expanded nucleobases and the stability of expanded duplexes, and therefore should be carefully considered when designing future expanded analogues.     

How do size-expanded DNA nucleobases enhance duplex stability? Computational analysis of the hydrogen-bonding and stacking ability of xDNA bases

Computational chemistry (B3LYP, MP2) is used to study the properties of size-expanded DNA nucleobases generated by inserting a benzene spacer into the natural nucleobases. Although the addition of the spacer does not significantly affect the hydrogen-bonding properties of natural nucleobases, the orientation of the base about the glycosidic bond necessary for Watson-Crick binding is destabilized, which could have implications for the selectivity of expanded bases, as well as the stability of expanded duplexes. Consideration of the (stacked) binding energies in the preferred relative orientation of natural and expanded nucleobases aligned according to their centers of mass reveals that the stacking within natural dimers can be increased by up to 50% upon expansion of one nucleobase and up to 90% upon expansion of two nucleobases. The implications of these findings to the stability of expanded duplexes were revealed by considering simplified models of natural and mixed duplexes composed of four nucleobases. Although intra- and interstrand interactions within double helices are typically less than those predicted when nucleobases are stacked according to their centers of mass, some nucleobases utilize their full stacking potential within double helices, where both intra- and interstrand interactions can be significant. Most importantly, increasing the size of nucleobases within the duplex significantly increases both intra- and interstrand stacking interactions. Specifically, some interactions are double the magnitude of the corresponding intrastrand interactions in natural helices, and even greater increases in interstrand interactions are sometimes found. Thus, our work suggests that mixed duplexes composed of natural bases hydrogen bound to expanded bases may exploit the increase in the inherent stacking ability of the expanded bases in more than one way and thereby afford duplexes with greater stability than natural DNA.     

The kinetic effects of decreasing catalytic activity on a simple rate expression

When kinetics are studied for a catalyzed reaction, the active catalyst surface available affects the results. This article illustrates the effect of a progressively decreasing available surface on the heterogeneous kinetics of a gas‐phase reaction in a closed system. The present contribution focuses on the effect of simple mth order surface deactivation on a simple nth order kinetic expression. The basic analysis and general results are unchanged if more complex equations are used. It is shown that there are certain common anomalous characteristics of kinetic expressions involving deactivation. In particular, the apparent rate constant and the apparent order are usually dependent on both the actual and initial pressures of the reactant, and the reaction may not go to completion. © 2000 John Wiley & Sons, Inc. Int J Chem Kinet 32: 7–16, 2000     

An efficient method for large scale configuration interaction calculations

An efficient method of handling large scale configuration interaction calculations is developed and applied to the H₂O molecule as a test case. The method, which is based upon matrix partitioning, is shown to be capable of calculating the ¹B₁ spectrum of H₂O to an accuracy level of 0.1 eV for each state with very moderate computational effort.     

Modeling the chemical step utilized by human alkyladenine DNA glycosylase: a concerted mechanism AIDS in selectively excising damaged purines

Human alkyladenine DNA glycosylase (AAG) initiates the repair of a wide variety of (neutral or cationic) alkylated and deaminated purines by flipping damaged nucleotides out of the DNA helix and catalyzing the hydrolytic N-glycosidic bond cleavage. Unfortunately, the limited number of studies on the catalytic pathway has left many unanswered questions about the hydrolysis mechanism. Therefore, detailed ONIOM(M06-2X/6-31G(d):AMBER) reaction potential energy surface scans are used to gain the first atomistic perspective of the repair pathway used by AAG. The lowest barrier for neutral 1,N(6)-ethenoadenine (εA) and cationic N(3)-methyladenine (3MeA) excision corresponds to a concerted (A(N)D(N)) mechanism, where our calculated ΔG(?) = 87.3 kJ mol(-1) for εA cleavage is consistent with recent kinetic data. The use of a concerted mechanism supports previous speculations that AAG uses a nonspecific strategy to excise both neutral (εA) and cationic (3MeA) lesions. We find that AAG uses nonspecific active site DNA-protein π-π interactions to catalyze the removal of inherently more difficult to excise neutral lesions, and strongly bind to cationic lesions, which comes at the expense of raising the excision barrier for cationic substrates. Although proton transfer from the recently proposed general acid (protein-bound water) to neutral substrates does not occur, hydrogen-bond donation lowers the catalytic barrier, which clarifies the role of a general acid in the excision of neutral lesions. Finally, our work shows that the natural base adenine (A) is further inserted into the AAG active site than the damaged substrates, which results in the loss of a hydrogen bond with Y127 and misaligns the general base (E125) and water nucleophile to lead to poor nucleophile activation. Therefore, our work proposes how AAG discriminates against the natural purines in the chemical step and may also explain why some damaged pyrimidines are bound but are not excised by this enzyme.     

Computational and experimental evidence for the structural preference of phenolic C-8 purine adducts

The structural and spectral properties of (ortho and para) C8-aryl-purine adducts formed from carbon attachment by phenolic toxins were investigated through DFT calculations and UV-vis absorbance and emission studies. The global minima of both the deoxyadenosine (dA) and deoxyguanosine (dG) adducts adopted a syn conformation about the glycosidic bond due to the presence of an O5'-H...N3 hydrogen bond, where the anti minima are 20-30 kJ mol-1 higher in energy. While the nucleobase adducts are planar, the presence of the deoxyribose sugar induces a twist about the carbon-carbon bond connecting the phenol and nucleobase rings. ortho-Phenolic adducts are less twisted than the corresponding para adducts due to stabilization provided by an intramolecular O-H...N7 bond. Solvation calculations, in combination with UV-vis studies, demonstrate that the structural preference is solvent dependent, where solvents with hydrogen-bonding abilities disrupt the intramolecular O-H...N7 hydrogen bond such that a greater degree of twist is observed, and less polar solvents stabilize the planar structure. Indeed, the ratio of twisted to planar conformers is estimated to be as large as 50:50 in some aprotic solvents. Thus, the combined experimental and computational approach has provided a greater understanding of the structure of the ortho- and para-dA and dG C-bonded phenoxyl adducts as the first step to understanding the biological consequences of this form of DNA damage.     

The degenerate properties of modified purine and pyrimidine DNA bases: a density functional study

Density-functional theory was used to study the properties (binding geometries and affinities for the natural DNA bases) of various degenerate nucleobases, which bind without discrimination to the purines or pyrimidines. The data for purine mimics (Z and K) indicates that although stronger binding strengths are calculated for pairs with cytosine compared with thymine, cytosine binds to a less stable tautomer of the nucleobase mimic. Indeed, the energy differences between the binding strengths and the tautomers effectively cancel and thereby provide a possible explanation for the observed degenerate properties of these molecules. Similar trends are found for the pyrimidine mimics (M and P); however, the energy differences do not cancel, even upon inclusion of environmental effects.