An ab initio theoretical study of electronic structure and properties of 2'-deoxyguanosine in gas phase and aqueous media

Drosophila Larval Brain Neoplasms Present Tumour-Type Dependent Genome Instability

Single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) are found at different rates in human cancer. To determine if these genetic lesions appear in Drosophila tumors we have sequenced the genomes of 17 malignant neoplasms caused by mutations in l(3)mbt, brat, aurA, or lgl We have found CNVs and SNPs in all the tumors. Tumor-linked CNVs range between 11 and 80 per sample, affecting between 92 and 1546 coding sequences. CNVs are in average less frequent in l(3)mbt than in brat lines. Nearly half of the CNVs fall within the 10 to 100Kb range, all tumor samples contain CNVs larger that 100 Kb and some have CNVs larger than 1Mb. The rates of tumor-linked SNPs change more than 20-fold depending on the tumor type: at late time points brat, l(3)mbt, and aurA and lgl lines present median values of SNPs/Mb of exome of 0.16, 0.48, and 3.6, respectively. Higher SNP rates are mostly accounted for by C > A transversions, which likely reflect enhanced oxidative stress conditions in the affected tumors. Both CNVs and SNPs turn over rapidly. We found no evidence for selection of a gene signature affected by CNVs or SNPs in the cohort. Altogether, our results show that the rates of CNVs and SNPs, as well as the distribution of CNV sizes in this cohort of Drosophila tumors are well within the range of those reported for human cancer. Genome instability is therefore inherent to Drosophila malignant neoplastic growth at a variable extent that is tumor type dependent.

Mechanisms and energetics for N-glycosidic bond cleavage of protonated 2'-deoxyguanosine and guanosine

Experimental and theoretical investigations suggest that hydrolysis of N-glycosidic bonds generally involves a concerted SN2 or a stepwise SN1 mechanism. While theoretical investigations have provided estimates for the intrinsic activation energies associated with N-glycosidic bond cleavage reactions, experimental measurements to validate the theoretical studies remain elusive. Here we report experimental investigations for N-glycosidic bond cleavage of the protonated guanine nucleosides, [dGuo+H](+) and [Guo+H](+), using threshold collision-induced dissociation (TCID) techniques. Two major dissociation pathways involving N-glycosidic bond cleavage, resulting in production of protonated guanine or the elimination of neutral guanine are observed in competition for both [dGuo+H](+) and [Guo+H](+). The detailed mechanistic pathways for the N-glycosidic bond cleavage reactions observed are mapped via electronic structure calculations. Excellent agreement between the measured and B3LYP calculated activation energies and reaction enthalpies for N-glycosidic bond cleavage of [dGuo+H](+) and [Guo+H](+) in the gas phase is found indicating that these dissociation pathways involve stepwise E1 mechanisms in analogy to the SN1 mechanisms that occur in the condensed phase. In contrast, MP2 is found to significantly overestimate the activation energies and slightly overestimate the reaction enthalpies. The 2'-hydroxyl substituent is found to stabilize the N-glycosidic bond such that [Guo+H](+) requires ∼25 kJ mol(-1) more than [dGuo+H](+) to activate the glycosidic bond.

Quantum theoretical study of cleavage of the glycosidic bond of 2'-deoxyadenosine: base excision-repair mechanism of DNA by MutY

The enzyme adenine DNA glycosylase, also called MutY, is known to catalyze base excision repair by removal of adenine from the abnormal 2'-deoxyadenosine:8-oxo-2'-deoxyguanosine pair in DNA. The active site of the enzyme was considered to consist of a glutamic acid residue along with two water molecules. The relevant reaction mechanism involving different barrier energies was studied theoretically. Molecular geometries of the various molecules and complexes involved in the reaction, e.g., the reactant, intermediate, and product complexes as well as transition states, were optimized employing density functional theory at the B3LYP/6-31G(d,p) level in the gas phase. It was followed by single-point energy calculations at the B3LYP/AUG-cc-pVDZ, BHandHLYP/AUG-cc-pVDZ, and MP2/AUG-cc-pVDZ levels in the gas phase. Single-point energy calculations were also carried out at the B3LYP/AUG-cc-pVDZ and BHandHLYP/AUG-cc-pVDZ levels in aqueous media as well as in the solvents chlorobenzene and dichloroethane. For the solvation calculations, the integral equation formalism of the polarizable continuum model (IEF-PCM) was employed. It is found that glutamic acid along with two water molecules would effectively cleave the glycosidic bond of adenosine by a new two-step reaction mechanism proposed here which is different from the three-step mechanism proposed by other authors earlier regarding the working mechanism of MutY.

Repair of O6-methylguanine to guanine by cysteine in the absence and presence of histidine and by cysteine thiolate anion: a quantum chemical study

O6-methylguanine (O6mG) is known to be a potential mutagenic modification of guanine as it mispairs with thymine in DNA and causes GC to AT transversion mutation. It is experimentally known that O6mG can be repaired to guanine by the protein O6-alkylguanine-DNA alkyltransferase (AGT), a cysteine residue being the main active site. In the present work, the mechanisms of repair of cis-O6-methylguanine (O6mG) to guanine due to its reaction with cysteine in the absence and presence of histidine and with cysteine thiolate anion were investigated theoretically using the B3LYP hybrid functional of density functional theory and the second order Møller-Plesset perturbation (MP2) theory. Reactant, intermediate and product complexes as well as transition states involved in these reactions were fully optimized at the B3LYP/6-31 + G* level of theory in the gas phase. The solvent effect of water was treated using the polarizable continuum model (PCM). Single point energy calculations were performed at the B3LYP/AUG-cc-pVDZ and MP2/6-31 + G* levels of theory in the gas phase and aqueous media. It is found that cysteine alone can repair the cis-O6mG to guanine, but the involvement of histidine along with cysteine lowers down the barrier energy significantly. However, when cysteine thiolate anion is used in place of cysteine, the barrier energy is strongly reduced. These results broadly support the suggestions based on experimental studies.

A quantum chemical study of repair of O6-methylguanine to guanine by tyrosine: evaluation of the winged helix-turn-helix model

The winged helix-turn-helix model for the repair of O6-MeG to guanine involving the reaction of O6-MeG with a tyrosine residue of the protein O6-alkylguanine-DNA alkyltransferase (AGT) was examined by studying the reaction mechanism and barrier energies. Molecular geometries of the species and complexes involved in the reaction, i.e. the reactant, intermediate and product complexes as well as transition states, were optimized employing density functional theory in gas phase. It was followed by single point energy calculations using density functional theory along with a higher basis set and second order M(phi)ller-Plesset perturbation theory (MP2) along with two different basis sets in gas phase and aqueous media. For the solvation calculations in aqueous media, the integral equation formalism of the polarizable continuum model (IEF-PCM) was employed. Vibrational frequency analysis was performed for each optimized structure and genuineness of transition states was ensured by visualizing the vibrational modes. It is found that tyrosine can repair O6-MeG to guanine by a two-step reaction. The present results have been compared with those obtained considering the helix-turn-helix model where the repair reaction primarily involves cysteine and occurs in a single-step. It is concluded that the repair through tyrosine envisaged in the winged helix-turn-helix model would be less efficient than that through cysteine envisaged in the helix-turn-helix model.

Analysis of actinomycin D-DNA model complexes using a quantum-chemical criterion: Mulliken overlap populations

The binding of the antitumoral drug actinomycin D to single- and double-stranded DNA was investigated using molecular modeling in the frame of MM+ molecular mechanics and AM1 semi-empirical method. Two other programs, especially conceived to analyze hydrogen-bonding patterns in biological macromolecules, HBexplore, based on geometrical criteria and SHB_interactions, based on quantum-chemical criteria (Mulliken overlap populations), were also used. The results account for the non-cooperative intercalative binding process previously investigated, and outline the contribution of specific hydrogen bonding as well as CH...O(N) and other atom-atom intermolecular interactions to the stabilization of the actinomycin D-DNA complexes. They also support the hemi-intercalation model proposed in literature for the actinomycin D-ssDNA complex.

Significant effects of phosphorylation on relative stabilities of DNA and RNA sugar radicals: remarkably high susceptibility of h-2' abstraction in RNA

The roles of nucleic acid radicals in DNA and RNA damage cannot be properly understood in the absence of knowledge of the C-H bond strengths depicting the energy cost to generate each of these radicals. However, previous theoretical studies on the relative energies of different nucleic acid radicals are not fully convincing mainly because of the use of oversimplified model compounds. In the present study we chose nucleoside 3',5'-bisphosphates as model compounds for DNA and RNA, in which the effects of both the nucleobase and phosphorylation were taken into consideration. Using the newly developed ONIOM-G3B3 methods, we calculated the gas-phase bond dissociation enthalpies and solution-phase bond dissociation free energies of all the carbohydrate C-H bonds in the model compounds. It was found that the monoanionic phosphate group (OPO3H-) was a better radical stabilization group than the OH group by 1.3 kcal/mol, whereas the neutral phosphate group (OPO3H2) was a significantly worse radical stabilization group than OH by 4.4 kcal/mol. Due to these reasons, the relative thermodynamic susceptibility of H-abstraction from deoxyribonucleotides and ribonucleotides varied considerably depending on the phosphorylation state and the charge carried by the phosphate groups. Strikingly, the bond dissociation free energy of C2'-H in ribonucleotides was dramatically lower than that of all the other C-H bonds by 5-6 kcal/mol regardless of the phosphorylation state and the charge carried by the phosphate group. This explained the previous experimental finding that radiation damage of RNA occurs mainly via H-abstraction at H-2'. A model study suggested that the strength of the hydrogen bonding interaction between the 2'-OH and 3-phosphate groups should dramatically increase from ribonucleoside 3',5'-bisphosphate to its C2' radical. The strengthened hydrogen bonding stabilized the C2' radical, rendering the C2'-H bond of RNA extraordinarily vulnerable to H-abstraction.

Flexible 5-guanidino-4-nitroimidazole DNA lesions: structures and thermodynamics

5-Guanidino-4-nitroimidazole (NI), derived from guanine oxidation by reactive oxygen and nitrogen species, contains an unusual flexible ring-opened structure, with nitro and guanidino groups which possess multiple hydrogen bonding capabilities. In vitro primer extension experiments with bacterial and mammalian polymerases show that NI incorporates C as well as A and G opposite the lesion, depending on the polymerase. To elucidate structural and thermodynamic properties of the mutagenic NI lesion, we have investigated the structure of the modified base itself and the NI-containing nucleoside with high-level quantum mechanical calculations and have employed molecular modeling and molecular dynamics simulations in solution for the lesion in B-DNA duplexes, with four partner bases opposite the NI. Our results show that NI adopts a planar structure at the damaged base level. However, in the nucleoside and in DNA duplexes, steric hindrance between the guanidino group and its linked sugar causes NI to be nonplanar. The NI lesion can adopt both syn and anti conformations on the DNA duplex level, with the guanidino group positioned in the DNA major and minor grooves, respectively; the specific preference depends on the partner base. On the basis of hydrogen bonding and stacking interactions, groove dimensions, and bending, we find that the least distorted NI-modified duplex contains partner C, consistent with observed incorporation of C opposite NI. However, hydrogen bonding interactions between NI and partner G or A are also found, which would be compatible with the observed mismatches.

Effects of Hydrogen Bonding on the Acidity of Adenine, Guanine, and Their 8-Oxo Derivatives

Complexes between ammonia, water, or hydrogen fluoride and adenine, guanine, or their 8-oxo derivatives are investigated using density-functional theory. The binding strengths of the neutral and (N9) anionic complexes are considered for a variety of purine binding sites. The effects of hydrogen-bonding interactions on the (N9) acidity of the purine derivatives are considered as a function of the molecule bound and the binding site. It is found that hydrogen-bonding interactions with one molecule can increase the acidity of purine derivatives by up to 60 kJ mol(-1). The (calculated) simultaneous effects of up to four molecules on the acidity of the purine derivatives are also considered. Our data suggest that the effects of more than one molecule on the acidity of the purines are generally less than the sum of the individual (additive) effects, where the magnitude of the deviation from additivity increases with the number, as well as the acidity, of molecules bound. Nevertheless, the increase in the acidity due to additional hydrogen-bonding interactions is significant, where the effect of two, three, or four hydrogen-bonding interactions can be as large as approximately 95, 115, and 130 kJ mol(-1), respectively. The present study provides a greater fundamental understanding of hydrogen-bonding interactions involving the natural purines, as well as those generated through oxidative DNA damage, which may aid the understanding of important biological processes.

Spiroiminodihydantoin lesions derived from guanine oxidation: structures, energetics, and functional implications

Reactive oxygen species present in the cell generate DNA damage. One of the major oxidation products of guanine in DNA, 8-oxo-7,8-dihydroguanine, formed by loss of two electrons, is among the most extensively studied base lesions. The further removal of two electrons from this product can yield spiroiminodihydantoin (Sp) R and S stereoisomers. Both in vitro and in vivo experiments have shown that the Sp stereoisomers are highly mutagenic, causing G --> T and G --> C transversions. Hence, they are of interest as examples of endogenous DNA damage that may initiate cancer. To interpret the mutagenic properties of the Sp lesions, an understanding of their structural properties is needed. To elucidate these structural effects, we have carried out computational investigations at the level of the Sp-modified base and nucleoside. At the base level, quantum mechanical geometry optimization studies have revealed exact mirror image symmetry of the R and S stereoisomers, with a near-perpendicular geometry of the two rings. At the nucleoside level, an extensive survey of the potential energy surface by molecular mechanics calculations using AMBER has provided three-dimensional potential energy maps. These maps reveal that the range and flexibility of the glycosidic torsion angles are significantly more restricted in both stereoisomeric adducts than in unmodified 2'-deoxyguanosine. The structural and energetic results suggest that the unusual geometric, steric, and hydrogen bonding properties of these lesions underlie their mutagenicity. In addition, stereoisomer-specific differences indicate the possibility that their processing by cellular replication and repair enzymes may be differentially affected by their absolute configuration.

DNA Nucleosides and Their Radical Anions: Molecular Structures and Electron Affinities

The deoxyribonucleosides have been studied to determine the properties of combinations of 2-deoxyribose with each of the isolated DNA bases for both neutral and anionic species. We have used a carefully calibrated theoretical method [Chem. Rev. 2002, 102, 231], employing the B3LYP hybrid Hartree-Fock/DFT functional with the DZP++ basis set. Predictions are made of the geometric parameters, adiabatic electron affinities, charge distributions based on natural population analysis, and decomposition enthalpy for the neutral and anionic forms of the four 2'-deoxyribonucleosides in DNA: 2'-deoxyriboadenosine (dA), 2'-deoxyribocytidine (dC), 2'-deoxyriboguanosine (dG), and 2'-deoxyribothymidine (dT). Geometric changes in the anions show that the glycosidic bond exhibits little change with excess charge for the guanosine but significant shortening for the adenosine and for the pyrimidines. The zero-point corrected adiabatic electron affinities in eV for each of the 2'-deoxyribonucleosides are as follows: 0.06, dA; 0.09, dG; 0.33, dC; and 0.44, dT. These values are uniformly greater than those of the corresponding isolated bases (-0.28, A; -0.07, G; 0.03, C; 0.20, T) and the isolated 2-deoxyribose (-0.38) at the same level of theory. The vertical detachment energies of dT and dC are substantial, 0.72 and 0.94 eV, and these anions should be observable. A high VDE, 0.91 eV, is also found for dA but its anion is unlikely to be stable due to the small AEA of 0.06 eV. The high VDE reflects the fact that the molecular structures of the anions and the corresponding neutral species are quite different. Valence character is displayed for the SOMOs of dA, dC, and dT, while some component of diffuse character is visible in the SOMO of dG. Further analysis of electronic changes upon electron attachment include an examination of the NPA charges, which show that in the neutral 2'-deoxyribonucleosides the sum of NPA charges for every base is the same, -0.28 with the sum of 2-deoxyribose charges being positive, +0.28. In the anions, the trend in charge division varies based on the nature of the excess electron in the anions. Thermodynamically, the overall enthalpy change for the reaction of water with the neutral nucleosides to give bases and ribose is approximately zero. The analogous decomposition is exothermic by 8 to 11 kcal mol-1 for the anions, indicating possible challenges for anionic gas-phase nucleoside exploration in the presence of water.

Altered structural fluctuations in duplex RNA versus DNA: a conformational switch involving base pair opening

DNA and RNA are known to have different structural properties. In the present study, molecular dynamics (MD) simulations on a series of RNA and DNA duplexes indicate differential structural flexibility for the two classes of oligonucleotides. In duplex RNA, multiple base pairs experienced local opening events into the major groove on the nanosecond time scale, while such events were not observed in the DNA simulations. Three factors are indicated to be responsible for the base opening events in RNA: solvent-base interactions, 2'OH(n)-O4'(n+1) intra-strand hydrogen bonding, and enhanced rigid body motion of RNA at the nucleoside level. Water molecules in the major groove of RNA contribute to initiation of base pair opening. Stabilization of the base pair open state is due to a 'conformational switch' comprised of 2'OH(n)-O4'(n+1) hydrogen bonding and a rigid body motion of the nucleoside moiety in RNA. This rigid body motion is associated with decreased flexibility of the glycosyl linkage and sugar moieties in A-form structures. The observed opening rates in RNA are consistent with the imino proton exchange experiments for AU base pairs, although not for GC base pairs, while structural and flexibility changes associated with the proposed conformational switch are consistent with survey data of RNA and DNA crystal structures. The possible relevance of base pair opening events in RNA to its many biological functions is discussed.