Stereochemical, structural, and thermodynamic origins of stability differences between stereoisomeric benzo[a]pyrene diol epoxide deoxyadenosine adducts in a DNA mutational hot spot sequence

Carcinogen-induced DNA structural distortion differences in the RAS gene isoforms; the importance of local sequence

BACKGROUND

Local sequence context is known to have an impact on the mutational pattern seen in cancer. The RAS genes and a smoking carcinogen, Benzo[a]pyrene diol epoxide (BPDE), have been utilised to explore these context effects. BPDE is known to form an adduct at the guanines in a number of RAS gene sites, KRAS codons 12, 13 and 14, NRAS codon 12, and HRAS codons 12 and 14.

RESULTS

Molecular modelling techniques, along with multivariate analysis, have been utilised to determine the sequence influenced differences between BPDE-adducted RAS gene sequences as well as the local distortion caused by the adducts.

CONCLUSIONS

We conclude that G:C > T:A mutations at KRAS codon 12 in the tumours of lung cancer patients (who smoke), proposed to be predominantly caused by BPDE, are due to the effect of the interaction methyl group at the C5 position of the thymine base in the KRAS sequence with the BPDE carcinogen investigated causing increased distortion. We further suggest methylated cytosine would have a similar effect, showing the importance of methylation in cancer development.

Towards a Better Understanding of Computational Models for Predicting DNA Methylation Effects at the Molecular Level

Human DNA is a very sensitive macromolecule and slight changes in the structure of DNA can have disastrous effects on the organism. When nucleotides are modified, or changed, the resulting DNA sequence can lose its information, if it is part of a gene, or it can become a problem for replication and repair. Human cells can regulate themselves by using a process known as DNA methylation. This methylation is vitally important in cell differentiation and expression of genes. When the methylation is uncontrolled, however, or does not occur in the right place, serious pathophysiological consequences may result. Excess methylation causes changes in the conformation of the DNA double helix. The secondary structure of DNA is highly dependent upon the sequence. Therefore, if the sequence changes slightly the secondary structure can change as well. These slight changes will then cause the doublestranded DNA to be more open and available in some places where large adductions can come in and react with the DNA base pairs. Computer models have been used to simulate a variety of biological processes including protein function and binding, and there is a growing body of evidence that in silico methods can shed light on DNA methylation. Understanding the anomeric effect that contributes to the structural and conformational flexibility of furanose rings through a combination of quantum mechanical and experimental studies is critical for successful molecular dynamic simulations.

Repair-Resistant DNA Lesions

The eukaryotic global genomic nucleotide excision repair (GG-NER) pathway is the major mechanism that removes most bulky and some nonbulky lesions from cellular DNA. There is growing evidence that certain DNA lesions are repaired slowly or are entirely resistant to repair in cells, tissues, and in cell extract model assay systems. It is well established that the eukaryotic DNA lesion-sensing proteins do not detect the damaged nucleotide, but recognize the distortions/destabilizations in the native DNA structure caused by the damaged nucleotides. In this article, the nature of the structural features of certain bulky DNA lesions that render them resistant to NER, or cause them to be repaired slowly, is compared to that of those that are good-to-excellent NER substrates. Understanding the structural features that distinguish NER-resistant DNA lesions from good NER substrates may be useful for interpreting the biological significance of biomarkers of exposure of human populations to genotoxic environmental chemicals. NER-resistant lesions can survive to replication and cause mutations that can initiate cancer and other diseases. Furthermore, NER diminishes the efficacy of certain chemotherapeutic drugs, and the design of more potent pharmaceuticals that resist repair can be advanced through a better understanding of the structural properties of DNA lesions that engender repair-resistance.

Base damage, local sequence context and TP53 mutation hotspots: a molecular dynamics study of benzo[a]pyrene induced DNA distortion and mutability

The mutational pattern for the TP53 tumour suppressor gene in lung tumours differs to other cancer types by having a higher frequency of G:C>T:A transversions. The aetiology of this differing mutation pattern is still unknown. Benzo[a]pyrene,diol epoxide (BPDE) is a potent cigarette smoke carcinogen that forms guanine adducts at TP53 CpG mutation hotspot sites including codons 157, 158, 245, 248 and 273. We performed molecular modelling of BPDE-adducted TP53 duplex sequences to determine the degree of local distortion caused by adducts which could influence the ability of nucleotide excision repair. We show that BPDE adducted codon 157 has greater structural distortion than other TP53 G:C>T:A hotspot sites and that sequence context more distal to adjacent bases must influence local distortion. Using TP53 trinucleotide mutation signatures for lung cancer in smokers and non-smokers we further show that codons 157 and 273 have the highest mutation probability in smokers. Combining this information with adduct structural data we predict that G:C>T:A mutations at codon 157 in lung tumours of smokers are predominantly caused by BPDE. Our results provide insight into how different DNA sequence contexts show variability in DNA distortion at mutagen adduct sites that could compromise DNA repair at well characterized cancer related mutation hotspots.

Adenine versus guanine DNA adducts of aristolochic acids: role of the carcinogen-purine linkage in the differential global genomic repair propensity

Computational modeling is employed to provide a plausible structural explanation for the experimentally-observed differential global genome repair (GGR) propensity of the ALII-N(2)-dG and ALII-N(6)-dA DNA adducts of aristolochic acid II. Our modeling studies suggest that an intrinsic twist at the carcinogen-purine linkage of ALII-N(2)-dG induces lesion site structural perturbations and conformational heterogeneity of damaged DNA. These structural characteristics correlate with the relative repair propensities of AA-adducts, where GGR recognition occurs for ALII-N(2)-dG, but is evaded for intrinsically planar ALII-N(6)-dA that minimally distorts DNA and restricts the conformational flexibility of the damaged duplex. The present analysis on the ALII adduct model systems will inspire future experimental studies on these adducts, and thereby may extend the list of structural factors that directly correlate with the propensity for GGR recognition.

Inhibition of HIV-1 reverse transcriptase-catalyzed synthesis by intercalated DNA Benzo[a]Pyrene 7,8-Dihydrodiol-9,10-Epoxide adducts

To aid in the characterization of the relationship of structure and function for human immunodeficiency virus type-1 reverse transcriptase (HIV-1 RT), this investigation utilized DNAs containing benzo[a]pyrene-7,8-dihydrodiol-9,10-epoxide (BPDE)-modified primers and templates as a probe of the architecture of this complex. BPDE lesions that differed in their stereochemistry around the C₁₀ position were covalently linked to N⁶-adenine and positioned in either the primer or template strand of a duplex template-primer. HIV-1 RT exhibited a stereoisomer-specific and strand-specific difference in replication when the BPDE-lesion was placed in the template versus the primer strand. When the C₁₀R-BPDE adduct was positioned in the primer strand in duplex DNA, 5 nucleotides from the 3΄ end of the primer terminus, HIV-1 RT could not fully replicate the template, producing truncated products; this block to further synthesis did not affect rates of dissociation or DNA binding affinity. Additionally, when the adducts were in the same relative position, but located in the template strand, similar truncated products were observed with both the C₁₀R and C₁₀S BPDE adducts. These data suggest that the presence of covalently-linked intercalative DNA adducts distant from the active site can lead to termination of DNA synthesis catalyzed by HIV-1 RT.

Free energy profiles of base flipping in intercalative polycyclic aromatic hydrocarbon-damaged DNA duplexes: energetic and structural relationships to nucleotide excision repair susceptibility

The crystal structure of Rad4/Rad23, the yeast homolog of the human nucleotide excision repair (NER) lesion recognition factor XPC-RAD23B ( Min , J. H. and Pavletich , N. P. ( 2007 ) Nature 449 , 570 - 575 ) reveals that the lesion-partner base is flipped out of the helix and binds to amino acids of the protein. This suggests the hypothesis that the flipping of this partner base must overcome a free energy barrier, which constitutes one element contributing to changes in the thermodynamic properties induced by the DNA damage and sensed by the recognition protein. We explored this hypothesis by computing complete flipping free energy profiles for two lesions derived from the procarcinogenic polycyclic aromatic hydrocarbons (PAHs), dibenzo[a,l]pyrene (DB[a,l]P) and benzo[a]pyrene (B[a]P), R-trans-anti-DB[a,l]P-N(6)-dA (R-DB[a,l]P-dA) and R-trans-anti-B[a]P-N(6)-dA (R-B[a]P-dA), and the corresponding unmodified duplex. The DB[a,l]P and B[a]P adducts differ in number and organization of their aromatic rings. We integrate these results with prior profiles for the R-trans-anti-DB[a,l]P-dG adduct ( Zheng , H. et al. ( 2010 ) Chem. Res. Toxicol. 23 , 1868 - 1870 ). All adopt conformational themes involving intercalation of the PAH aromatic ring system into the DNA duplex; however, R-DB[a,l]P-dA and R-B[a]P-dA intercalate from the major groove, while R-DB[a,l]P-dG intercalates from the minor groove. These structural differences produce different computed van der Waals stacking interaction energies between the flipping partner base with the lesion aromatic ring system and adjacent bases; we find that the better the stacking, the higher the relative flipping free energy barrier and hence lower flipping probability. The better relative NER susceptibilities correlate with greater ease of flipping in these three differently intercalated lesions. In addition to partner base flipping, the Rad4/Rad23 crystal structure shows that a protein-β-hairpin, BHD3, intrudes from the major groove side between the DNA strands at the lesion site. We present a molecular modeling study for the R-DB[a,l]P-dG lesion in Rad4/Rad23 showing BHD3 β-hairpin intrusion with lesion eviction, and we hypothesize that lesion steric effects play a role in the recognition of intercalated adducts.

Adenine-DNA adducts derived from the highly tumorigenic Dibenzo[a,l]pyrene are resistant to nucleotide excision repair while guanine adducts are not

The structural origins of differences in susceptibilities of various DNA lesions to nucleotide excision repair (NER) are poorly understood. Here we compared, in the same sequence context, the relative NER dual incision efficiencies elicited by two stereochemically distinct pairs of guanine (N(2)-dG) and adenine (N(6)-dA) DNA lesions, derived from enantiomeric genotoxic diol epoxides of the highly tumorigenic fjord region polycyclic aromatic hydrocarbon dibenzo[a,l]pyrene (DB[a,l]P). Remarkably, in cell-free HeLa cell extracts, the guanine adduct with R absolute chemistry at the N(2)-dG linkage site is ∼35 times more susceptible to NER dual incisions than the stereochemically identical N(6)-dA adduct. For the guanine and adenine adducts with S stereochemistry, a similar but somewhat smaller effect (factor of ∼15) is observed. The striking resistance of the bulky N(6)-dA in contrast to the modest to good susceptibilities of the N(2)-dG adducts to NER is interpreted in terms of the balance between lesion-induced DNA distorting and DNA stabilizing van der Waals interactions in their structures, that are partly reflected in the overall thermal stabilities of the modified duplexes. Our results are consistent with the hypothesis that the high genotoxic activity of DB[a,l]P is related to the formation of NER-resistant and persistent DB[a,l]P-derived adenine adducts in cellular DNA.

Nucleotide excision repair efficiencies of bulky carcinogen-DNA adducts are governed by a balance between stabilizing and destabilizing interactions

The nucleotide excision repair (NER) machinery, the primary defense against cancer-causing bulky DNA lesions, is surprisingly inefficient in recognizing certain mutagenic DNA adducts and other forms of DNA damage. However, the biochemical basis of resistance to repair remains poorly understood. To address this problem, we have investigated a series of intercalated DNA-adenine lesions derived from carcinogenic polycyclic aromatic hydrocarbon (PAH) diol epoxide metabolites that differ in their response to the mammalian NER apparatus. These stereoisomeric PAH-derived adenine lesions represent ideal model systems for elucidating the effects of structural, dynamic, and thermodynamic properties that determine the recognition of these bulky DNA lesions by NER factors. The objective of this work was to gain a systematic understanding of the relation between aromatic ring topology and adduct stereochemistry with existing experimental NER efficiencies and known thermodynamic stabilities of the damaged DNA duplexes. For this purpose, we performed 100 ns molecular dynamics studies of the lesions embedded in identical double-stranded 11-mer sequences. Our studies show that, depending on topology and stereochemistry, stabilizing PAH-DNA base van der Waals stacking interactions can compensate for destabilizing distortions caused by these lesions that can, in turn, cause resistance to NER. The results suggest that the balance between helix stabilizing and destabilizing interactions between the adduct and nearby DNA residues can account for the variability of NER efficiencies observed in this class of PAH-DNA lesions.

Structure and stability of DNA containing an aristolactam II-dA lesion: implications for the NER recognition of bulky adducts

Aristolochic acids I and II are prevalent plant toxicants found in the Aristolochiaceae plant family. Metabolic activation of the aristolochic acids leads to the formation of a cyclic N-hydroxylactam product that can react with the peripheral amino group of purine bases generating bulky DNA adducts. These lesions are mutagenic and established human carcinogens. Interestingly, although AL-dG adducts progressively disappear from the DNA of laboratory animals, AL-dA lesions has lasting persistence in the genome. We describe here NMR structural studies of an undecameric duplex damaged at its center by the presence of an ALII-dA adduct. Our data establish a locally perturbed double helical structure that accommodates the bulky adduct by displacing the counter residue into the major groove and stacking the ALII moiety between flanking bases. The presence of the ALII-dA perturbs the conformation of the 5′-side flanking base pair, but all other pairs of the duplex adopt standard conformations. Thermodynamic studies reveal that the lesion slightly decreases the energy of duplex formation in a sequence-dependent manner. We discuss our results in terms of its implications for the repair of ALII-dA adducts in mammalian cells.

Mechanism of epoxide hydrolysis in microsolvated nucleotide bases adenine, guanine and cytosine: a DFT study

Six water molecules have been used for microsolvation to outline a hydrogen bonded network around complexes of ethylene epoxide with nucleotide bases adenine (EAw), guanine (EGw) and cytosine (ECw). These models have been developed with the MPWB1K-PCM/6-311++G(3df,2p)//MPWB1K/6-31+G(d,p) level of DFT method and calculated S(N)2 type ring opening of the epoxide due to amino group of the nucleotide bases, viz. the N6 position of adenine, N2 position of guanine and N4 position of cytosine. Activation energy (E(act)) for the ring opening was found to be 28.06, 28.64, and 28.37 kcal mol(-1) respectively for EAw, EGw and ECw. If water molecules were not used, the reactions occurred at considerably high value of E(act), viz. 53.51 kcal mol(-1) for EA, 55.76 kcal mol(-1) for EG and 56.93 kcal mol(-1) for EC. The ring opening led to accumulation of negative charge on the developing alkoxide moiety and the water molecules around the charge localized regions showed strong hydrogen bond interactions to provide stability to the intermediate systems EAw-1, EGw-1 and ECw-1. This led to an easy migration of a proton from an activated water molecule to the alkoxide moiety to generate a hydroxide. Almost simultaneously, a proton transfer chain reaction occurred through the hydrogen bonded network of water molecules and resulted in the rupture of one of the N-H bonds of the quaternized amino group. The highest value of E(act) for the proton transfer step of the reaction was 2.17 kcal mol(-1) for EAw, 2.93 kcal mol(-1) for EGw and 0.02 kcal mol(-1) for ECw. Further, the overall reaction was exothermic by 17.99, 22.49 and 13.18 kcal mol(-1) for EAw, EGw and ECw, respectively, suggesting that the reaction is irreversible. Based on geometric features of the epoxide-nucleotide base complexes and the energetics, the highest reactivity is assigned for adenine followed by cytosine and guanine. Epoxide-mediated damage of DNA is reported in the literature and the present results suggest that hydrated DNA bases become highly S(N)2 active on epoxide systems and the occurrence of such reactions can inflict permanent damage to the DNA.

Intercalative conformations of the 14R (+)- and 14S (-)-trans-anti-DB[a,l]P-N⁶-dA adducts: molecular modeling and MD simulations

Among the polycyclic aromatic hydrocarbon class of chemical carcinogens, dibenzo[a,l]pyrene (DB[a,l]P) is the most potent tumorigen that has been identified to date. Structurally, it is bulky with six aromatic rings, and it contains the nonplanar fjord-region. The conformational properties of DB[a,l]P-derived DNA adducts responsible for its extraordinary carcinogenicity are hence of great interest. We have carried out molecular modeling and MD simulations for the 14R (+)- and 14S (-)-trans-anti-DB[a,l]P-N⁶-dA adducts derived from the reactions of the DB[a,l]P diol epoxides with adenine in double-stranded DNA. The structures are based on the classically intercalated NMR solution structures of the analogous fjord-region benzo[c]phenanthrene-derived-N⁶-dA adducts. One objective was to gain insight on the impact of the more bulky DB[a,l]P ring system on the structural characteristics of the intercalative adduct conformations. A further objective was to elucidate the effect of the flexible twist associated with the sterically hindered aromatic ring in the fjord-region on the intercalated conformations, for comparison with the intercalated but planar bay-region benzo[a]pyrene-derived-N⁶-dA adducts. For the DB[a,l]P-N⁶-dA adducts, our results show that the 14R (+)-adduct is more favorably intercalated on the 5'-side of the modified adenine than the stereoisomeric 14S (-)-adduct, intercalated on its 3'-side. The 14R (+)-adduct manifests better van der Waals stacking interactions with flanking base pairs, less perturbed Watson-Crick hydrogen bonding, less local groove enlargement, less unwinding, and a lower solvent exposure than the 14S (-)-adduct. These structural findings are consistent with observed thermodynamic melting data, UV absorption properties, and fluorescence quenching studies. By contrast, the NMR solution structures for the analogous but less bulky B[c]Ph-derived adducts reveal no such stereoisomeric effect, while the planar bay-region benzo[a]pyrene-derived-N⁶-dA adducts do. Differences in nucleotide excision repair susceptibilities of the fjord and bay region adducts stem from distinctions in their intercalative conformations, produced by the intrinsic topological variations in their polycyclic aromatic ring systems.

Novel enzymatic function of DNA polymerase nu in translesion DNA synthesis past major groove DNA-peptide and DNA-DNA cross-links

DNA polymerase ν (POLN or pol ν) is a newly discovered A family polymerase that generates a high error rate when incorporating nucleotides opposite dG; its translesion DNA synthesis (TLS) capability has only been demonstrated for high fidelity replication bypass of thymine glycol lesions. In the current investigation, we describe a novel TLS substrate specificity of pol ν, demonstrating that it is able to bypass exceptionally large DNA lesions whose linkages are through the DNA major groove. Specifically, pol ν catalyzed efficient and high fidelity TLS past peptides linked to N⁶-dA via a reduced Schiff base linkage with a γ-hydroxypropano-dA. Additionally, pol ν could bypass DNA interstrand cross-links with linkage between N⁶-dAs in complementary DNA strands. However, the chemically identical DNA−peptide and DNA interstrand cross-links completely blocked pol ν when they were located in the minor groove via a N²-dG linkage. Furthermore, we showed that pol ν incorporated a nucleotide opposite the 1,N⁶-etheno-dA (εdA) in an error-free manner and (+)-trans-anti-benzo[a]pyrene-7,8-dihydrodiol 9,10-epoxide-dA [(+)-BPDE-dA] in an error-prone manner, albeit with a greatly reduced capability. Collectively, these data suggest that although pol ν bypass capacity cannot be generalized to all major groove DNA adducts, this polymerase could be involved in TLS when genomic replication is blocked by extremely large major groove DNA lesions. In view of the recent observation that pol ν may have a role in cellular tolerance to DNA cross-linking agents, our findings provide biochemical evidence for the potential functioning of this polymerase in the bypass of some DNA−protein and DNA−DNA cross-links.

H2AX phosphorylation in A549 cells induced by the bulky and stable DNA adducts of benzo[a]pyrene and dibenzo[a,l]pyrene diol epoxides

Early events in the cellular response to DNA damage, such as double strand breaks, rely on lesion recognition and activation of proteins involved in maintenance of genomic stability. One important component of this process is the phosphorylation of the histone variant H2AX. To investigate factors explaining the variation in carcinogenic potency between different categories of polycyclic aromatic hydrocarbons (PAHs), we have studied the phosphorylation of H2AX (H2AXgamma). A549 cells were exposed to benzo[a]pyrene diol epoxide [(+)-anti-BPDE] (a bay-region PAH) and dibenzo[a,l]pyrene diol epoxide [(-)-anti-DBPDE] (a fjord-region PAH) and H2AXgamma was studied using immunocytochemistry and Western blot. Hydrogen peroxide (H(2)O(2)) was used to induce oxidative DNA damage and strand breaks. As showed with single cell gel electrophoresis, neither of the diol epoxides resulted in DNA strand breaks relative to H(2)O(2). Visualisation of H2AXgamma formation demonstrated that the proportion of cells exhibiting H2AXgamma staining at 1h differed between BPDE, 40% followed by a decline, and DBPDE, <10% followed by an increase. With H(2)O(2) treatment, almost all cells demonstrated H2AXgamma at 1h. Western blot analysis of the H2AXgamma formation also showed concentration and time-dependent response patterns. The kinetics of H2AXgamma formation correlated with the previously observed kinetics of elimination of BPDE and DBPDE adducts. Thus, the extent of H2AXgamma formation and persistence was related to both the number of adducts and their structural features.

NMR and molecular modeling study, as complementary techniques to capillary electrophoresis method to elucidate the separation mechanism of linezolid enantiomers

NMR study and molecular modeling were performed to improve the level of understanding of the chiral recognition process occurring between linezolid and anionic single-isomer cyclodextrin-heptakis-(2,3-diacetyl-6-sulfato)-beta-cyclodextrin (HDAS-beta-CD). NMR spectrometry allowed to estimate the stoichiometry of the complexes between HDAS-beta-CD and S- or R-linezolid and to determine the binding constants. The 1:1 complex stoichiometry was detected in millimolar concentrations and the mode of binding was proposed. The binding constants Ka of the complexes were of the order of 30-80 M(-1). Molecular dynamic simulations of 40ns for four complexes and calculations of binding free energies were performed. These calculations allowed determining the mode of binding of linezolid to HDAS-beta-CD and explaining the binding enantioselectivity.

Dynamics of a benzo[a]pyrene-derived guanine DNA lesion in TGT and CGC sequence contexts: enhanced mobility in TGT explains conformational heterogeneity, flexible bending, and greater susceptibility to nucleotide excision repair

The nucleotide excision repair (NER) machinery excises a variety of bulky DNA lesions, but with varying efficiencies. The structural features of the DNA lesions that govern these differences are not well understood. An intriguing model system for studying structure-function relationships in NER is the major adduct derived from the reaction of the highly tumorigenic metabolite of benzo[a]pyrene, (+)-anti-benzo[a]pyrene diol epoxide, with the exocyclic amino group of guanine ((+)-trans-anti-[BP]-N(2)-dG, or G*). The rates of incision of the stereochemically identical lesions catalyzed by the prokaryotic UvrABC system was shown to be greater by a factor of 2.3+/-0.3 in the TG*T than in the CG*C sequence context [Biochemistry 46 (2007) 7006-7015]. Here we employ molecular dynamics simulations to elucidate the origin of the greater excision efficiency in the TG*T case and, more broadly, to delineate structural parameters that enhance NER. Our results show that the BP aromatic ring system is 5'-directed along the modified strand in the B-DNA minor groove in both sequence contexts. However, the TG*T modified duplex is much more dynamically flexible, featuring more perturbed and mobile Watson-Crick hydrogen bonding adjacent to the lesion, a greater impairment in stacking interactions, more dynamic local roll/bending, and more minor groove flexibility. These characteristics explain a number of experimental observations concerning the (+)-trans-anti-[BP]-N(2)-dG adduct in double-stranded DNA with the TG*T sequence context: its conformational heterogeneity in NMR solution studies, its highly flexible bend, and its lower thermal stability. By contrast, the CG*C modified duplex is characterized by a single BP conformation and a rigid bend. While current recognition models of bulky lesions by NER factors have stressed the importance of impaired Watson-Crick pairing/stacking and bending, our results highlight the likelihood of an important role for the local dynamics in the vicinity of the lesion.

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.

Multiple binding modes of the camptothecin family to DNA oligomers

The binding constants of camptothecin, topotecan and its lactone ring-opened carboxylate derivative to DNA octamers were measured by UV and NMR spectroscopy. The self-association of topotecan (TPT) was also measured. The carboxylate form of TPT binds in the same way as the lactone, but more weakly. Titration of TPT into d(GCGATCGC)2 shows a preferred location stacked onto the terminal G1 base. However, the intermolecular NOEs cannot be reconciled with a single conformation of the complex, and suggest a model of a limited number of conformations in fast exchange. MD calculations on four pairs of starting structures with TPT stacked onto the G1-C8 base pair in different orientations were therefore performed. The use of selected experimental "docking" restraints yielded ten MD trajectories covering a wide conformational space. From a combination of calculated free energies, NOEs and chemical shifts, some of the structures produced could be eliminated, and it is concluded that the data are consistent with two major families of conformations in fast exchange. One of these is the conformation found in a crystal of a TPT/DNA/topoisomerase I ternary complex [Proc. Natl. Acad. Sci. USA 2002, 99, 15 387-15 392].

Oxetane Locked Thymidine in the Dickerson-Drew Dodecamer Causes Local Base Pairing Distortions—An NMR Structure and Hydration Study

The introduction of a North-type sugar conformation constrained oxetane T block, 1-(1',3'-O-anhydro-beta-D-psicofuranosyl) thymine, at the T(7) position of the self-complementary Dickerson-Drew dodecamer, d[(5'-C(1)G(2)C(3)G(4)A(5)A(6)T(7)T(8)C(9)G(10)C(11)G(12)-3')](2), considerably perturbs the conformation of the four central base pairs, reducing the stability of the structure. UV spectroscopy and 1D NMR display a drop in melting temperature of approximately 10 degrees C per modification for the T(7) oxetane modified duplex, where the T(7) block has been introduced in both strands, compared to the native Dickerson-Drew dodecamer. The three dimensional structure has been determined by NMR spectroscopy and has subsequently been compared with the results of 2.4 ns MD simulations of the native and the T(7) oxetane modified duplexes. The modified T(7) residue is found to maintain its constrained sugar- and the related glycosyl torsion conformations in the duplex, resulting in staggered and stretched T(7).A(6) and A(6).T(7) non-linear base pairs. The stacking is less perturbed, but there is an increased roll between the two central residues compared to the native counterpart, which is compensated by tilts of the neighboring base steps. The one dimensional melting profile of base protons of the T(7) and T(8) residues reveals that the introduction of the North-type sugar constrained thymine destabilizes the core of the modified duplex, promoting melting to start simultaneously from the center as well as from the ends. Temperature dependent hydration studies by NMR demonstrate that the central T(7).A(6)/A(6).T(7) base pairs of the T(7) oxetane modified Dickerson-Drew dodecamer have at least one order of magnitude higher water exchange rates (correlated to the opening rate of the base pair) than the corresponding base pairs in the native duplex.

Deformations of promoter DNA bound to carcinogens help interpret effects on TATA-element structure and activity

The TATA-box binding protein (TBP) is required by eukaryotic RNA polymerases for correct transcription initiation. TBP binds to the minor groove of an 8 base pair (bp) DNA-promoter element known as the TATA box and severely bends the TATA box. The promoter-DNA substrate can be damaged by components present in the cell or the environment to produce covalent carcinogen-DNA adducts. These may lead to transcription blockage or unfaithful transcription. Benzo[a]pyrene (BP) is a widespread environmental chemical carcinogen which can be metabolically converted to DNA-reactive enantiomeric (+) and (-)-anti-benzo[a]pyrene diol epoxides (BPDEs). Recent experimental studies of a pair of stereoisomeric adenine adducts, derived from (+) and (-)-anti-BPDEs, have revealed how these lesions influence the complexation of TBP with the TATA box. Depending on the adduct's location in the TATA box and its stereochemistry, the stability of monomeric TATA-TBP complexes was found to increase or decrease relative to the unmodified DNA. We report here analyses of molecular-dynamics simulations to interpret these findings. Structural analyses of 12 DNA-protein systems representing different combinations of adduct stereoisomer type and placement within the promoter reveal that the location of the adduct within the TATA octamer determines whether the stability of TATA-TBP complexes is increased or decreased. The effect on binding stability can be interpreted in terms of conformational freedom and major-groove space available to BP due to the hydrogen bonds and inserted phenylalanines of the TATA-TBP complex; that is, depending on the position of the adenine to which BP is covalently bound, BP can be accommodated in an intercalated or major-groove orientation with ease or with difficulty (due to interference with TATA-TBP interactions). The unravelled structures and interactions thus reveal the effect of different adduct locations on TATA-TBP complex formation and suggest how transcription initiation may be affected by the presence of a bulky BP.

Relative free energy of binding and binding mode calculations of HIV-1 RT inhibitors based on dock-MM-PB/GS

Tetrahydroimidazo-[4,5,l-jk][1,4]-benzodiazepin-2-(1H)-one (TIBO) derivatives are important nonnucleoside human immunodeficiency virus-1 reverse transcriptase inhibitors (NNRTI). Several TIBO derivatives have shown high potency to inhibit reverse transcriptase (RT) and one (Tivirapine) has entered into clinical trials. The free energy of binding (FEB) is a numerical way to express the binding affinity of a ligand to its receptor and has been applied in screening candidates in rational drug design. In this work, the FEB of 42 TIBOs in RT was studied. Relative FEB is expressed in the form of a linear combination of vdW, electrostatic, solvation, and nonpolar solvation energy terms. The predicted FEB activity of the TIBOs studied has a good correlation (r(2) = 0.8680, q(2) = 0.8298) with respect to the experimental activity (pIC(50)). Based on the data reported here, the Finite Difference Poisson Boltzmann with a Gaussian Smooth Dielectric Constant Function method (PB/GS) solvation energy term is very important in predicting the binding affinity of TIBOs in RT. In summary, the Dock-Molecular Mechanics (MM)-PB/GS method is a promising technique in predicting ligand/receptor binding affinity and it can be used to screen relatively large sets of molecules in a reasonable amount of computer time.

Dynamic conformational heterogeneities of carcinogen-DNA adducts and their mutagenic relevance

Arylamines and polycyclic aromatic hydrocarbons (PAHs), which are known as "bulky" carcinogens, have been studied extensively and upon activation in vivo, react with cellular DNA to form DNA-adducts. The available structure data accumulated thus far has revealed that conformational heterogeneity is a common theme among duplex DNA modified with these carcinogens. Several conformationally diverse structures have been elucidated and found to be in equilibrium in certain cases. The dynamics of the heterogeneity appear to be modulated by the nature of the adduct structure and the base sequences neighboring the lesion site. These can be termed as "adduct- and sequence-induced conformational heterogeneities," respectively. Due to the small energy differences, the population levels of these conformers could readily be altered within the active sites of repair or replicate enzymes. Thus, the complex role of "enzyme-induced conformational heterogeneity" must also be taken into consideration for the establishment of a functional structure-mutation relationship. Ultimately, a major challenge in mutation structural biology is to carry out adduct- and site-specific experiments in a conformationally specific manner within biologically relevant environments. Results from such experiments should provide an accurate account of how a single chemically homogenous adduct gives rise to complex multiple mutations, the earliest step in the induction of cancer.

Molecular modeling of the major benzo[a]pyrene N2-dG adduct in cases where mutagenesis results are known in double stranded DNA

The potent mutagen/carcinogen benzo[a]pyrene (B[a]P) is metabolically activated to (+)-anti-B[a]PDE, which induces a full spectrum of mutations (e.g. GC-->TA, GC-->AT, etc.). One hypothesis for this complexity is that different mutations are induced by different conformations of its major adduct [+ta]-B[a]P-N2-dG when bypassed during DNA replication (probably by different DNA polymerases). Previous molecular modeling studies suggested that B[a]P-N2-dG adducts can in principle adopt at least 16 potential conformational classes in ds-DNA. Herein we report on molecular modeling studies with the eight conformations most likely to be relevant to base substitution mutagenesis in 10 cases where mutagenesis has been studied in ds-DNA plasmids in E. coli with B[a]P-N2-dG adducts of differing stereoisomers and DNA sequence contexts, as well as in five cases where the conformation is known by NMR. Of the approximately 11,000 structures generated in this study, the computed lowest energy structures are reported for 120 cases (i.e. eight conformations and 15 examples), and their conformations compared. Of the eight conformations, four are virtually always computed to be high in energy. The remaining four lower energy conformations include two with the BP moiety in the minor groove (designated: BPmi5 and BPmi3), and two base-displaced conformations, one with the dG moiety in the major groove (designated: Gma5) and one with the dG in the minor groove (designated: Gmi3). Interestingly, these four are the only conformations that have been observed for B[a]P-N2-dG adducts in NMR studies. Independent of sequence contexts and adduct stereochemistry, BPmi5 structures tend to look reasonably similar, as do BPmi3 structures, while the base-displaced structures Gma5 and BPmi3 tend to show greater variability in structure. A correlation was sought between modeling and mutagenesis results in the case of the low energy conformations BPmi5, BPmi3, Gma5 and Gma3. Plots of log[(G-->T)/(G-->A)] versus energy[(conformation X)-(conformation Y)] were constructed for all six pairwise combinations of these four conformations, and the only plot giving a straight line involved Gma5 and Gmi3. While this finding is striking, its significance is unclear (as discussed).

Understanding discrimination by the ribosome: stability testing and groove measurement of codon-anticodon pairs

The ribosome must discriminate between correct and incorrect tRNAs with sufficient speed and accuracy to sustain an adequate rate of cell growth. Here, we report the results of explicit solvent molecular dynamics simulations, which address the mechanism of discrimination by the ribosome. The universally conserved 16S rRNA base A1493 and the kink in mRNA between A and P sites amplify differences in stability between cognate and near-cognate codon-anticodon pairs. Destabilization by the mRNA kink also provides a geometric explanation for the higher error rates observed for mismatches in the first codon position relative to mismatches in the second codon position. For more stable near-cognates, the repositioning of the universally conserved bases A1492 and G530 results in increased solvent exposure and an uncompensated loss of hydrogen bonds, preventing correct codon-anticodon-ribosome interactions from forming.

Simulating structural and thermodynamic properties of carcinogen-damaged DNA

A pair of stereoisomeric covalent adducts to guanine in double-stranded DNA, derived from the reaction of mutagenic and tumorigenic metabolites of benzo[a]pyrene, have been well characterized structurally and thermodynamically. Both high-resolution NMR solution structures and an array of thermodynamic data are available for these 10S (+)- and 10R (-)-trans-anti -[BP]-N(2)-dG adducts in double-stranded deoxyoligonucleotides. The availability of experimentally well-characterized duplexes containing these two stereoisomeric guanine adducts provides an opportunity for evaluating the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method for computing thermodynamic properties from molecular dynamics ensembles. We have carried out 3-ns molecular dynamics simulations, using NMR solution structures as the starting models for the 10S (+)- and 10R (-)-trans-anti-dG adducts in a DNA duplex 11-mer using AMBER 6.0. We employed the MM-PBSA method to compute the free energies, enthalpies, and entropies of the two adducts. Our complete thermodynamic analysis agrees quite well with the full experimental thermodynamic characterization of these adducts, showing essentially equal stabilities of the two adducts. We also calculated the nuclear Overhauser effect (NOE) distances from the molecular dynamics trajectories, and compared them against the experimental NMR-derived NOE distances. Our results showed that the simulated structures are in good agreement with the NMR experimental NOE data. Furthermore, the molecular dynamics simulations provided new structural and biological insights. Specifically, the puzzling observation that the BP aromatic ring system in the 10S (+)-trans-anti-dG adduct is more exposed to the aqueous solvent than the 10R (-)-trans-anti-dG adduct, is rationalized in terms of the adduct structures. The structural and thermodynamic features of these stereoisomeric adducts are also discussed in relation to their reported low susceptibilities to nucleotide excision repair.

Theoretical methods for the simulation of nucleic acids

Different theoretical methods for the description of nucleic acid structures are reviewed. Firstly, we introduce the concept of classical force-field in the context of nucleic acid structures, discussing their accuracy. We then examine theoretical approaches to the description of nucleic acids based on: i) a rigid or quasi-rigid description of the molecule, ii) molecular mechanics optimization, and iii) molecular dynamics. Special emphasis is made ion current state of the art molecular dynamics simulations of nucleic acids structures.

Relating repair susceptibility of carcinogen-damaged DNA with structural distortion and thermodynamic stability

A key issue in the nucleotide excision repair (NER) of bulky carcinogen-DNA adducts is the ability of the NER machinery to recognize and repair certain adducts while failing to repair others. Unrepaired adducts can survive to cause mutations that initiate the carcinogenic process. Benzo[c]phenanthrene (B[c]Ph), a representative fjord region polycyclic aromatic hydrocarbon, can be metabolically activated to the enantiomeric benzo[c]phenanthrene diol epoxides (B[c]PhDEs), (+)-(1S,2R,3R,4S)-3,4- dihydroxy-1,2-epoxy-1,2,3,4-tetrahydrobenzo[c]phenanthrene and the corresponding (-)-(1R,2S,3S,4R) isomer. These react predominantly with adenine residues in DNA to produce the stereoisomeric 1R (+)- and 1S (-)-trans-anti-B[c]Ph-N6-dA adducts. Duplexes containing the 1R (+) or 1S (-) B[c]Ph-dA adduct in codon 61 of the human N-ras mutational hotspot sequence CA*A, with B[c]Ph modification at A*, are not repaired by the human NER system. However, the analogous stereoisomeric DNA adducts of the bay region benzo[a]pyrene diol epoxide (B[a]PDE), 10S (+)- and 10R (-)-trans-anti-B[a]P-N6-dA, are repaired in the same base sequence. In order to elucidate structural and thermodynamic origins of this phenomenon, we have carried out a 2 ns molecular dynamics simulation for the 1R (+)- and 1S (-)-trans-anti-B[c]Ph-N6-dA adducts in an 11mer duplex containing the human N-ras codon 61 sequence, and compared these results with our previous study of the B[a]P-dA adducts in the same sequence. The molecular mechanics Poisson- Boltzmann surface area (MM-PBSA) method was applied to calculate the free energies of the pair of stereoisomeric B[c]Ph-dA adducts, and a detailed structural analysis was carried out. The different repair susceptibilities of the B[a]P-dA adducts and the B[c]Ph-dA adducts can be attributed to different degrees of distortion, stemming from combined effects of differences in the quality of Watson-Crick hydrogen bonding, unwinding, stretching and helix backbone perturbations. These differences are due to the different intrinsic topologies of the rigid, planar bay region adducts versus the twisted, sterically hindered fjord region adducts.

Accurate representation of B-DNA double helical structure with implicit solvent and counterions

High-resolution nuclear magnetic resonance (NMR) and crystallographic data have been taken to refine the force field used in the torsion angle space nucleic acids molecular mechanics program DUPLEX. The population balance deduced from NMR studies of two carcinogen-modified DNA conformers in equilibrium was used to fine tune a sigmoidal, distance-dependent dielectric function so that reasonable relative energies could be obtained. In addition, the base-pair and backbone geometry from high-resolution crystal structures of the Dickerson-Drew dodecamer was used to re-evaluate the deoxyribose pseudorotation profile and the Lennard-Jones nonbonded energy terms. With a modified dielectric function that assumes a very steep distance-dependent form, a deoxyribose pseudorotation profile with reduced energy barriers between C2'- and C3'-endo minima, and a shift of the Lennard-Jones potential energy minimum to a distance approximately 0.4 A greater than the sum of the van der Waals' radii, the sequence-dependent conformational features of the Dickerson-Drew dodecamer in both the solid state and the aqueous liquid crystalline phase are well reproduced. The robust performance of the revised force field, in conjunction with its efficiency through implicit treatment of solvent and counterions, provides a valuable tool for elucidating conformations and structure-function relationships of DNA, including those of molecules modified by carcinogens and other ligands.

Molecular dynamics simulation of nucleic acids: successes, limitations, and promise

In the last five years we have witnessed a significant increase in the number publications describing accurate and reliable all-atom molecular dynamics simulations of nucleic acids. This increase has been facilitated by the development of fast and efficient methods for treating the long-range electrostatic interactions, the availability of faster parallel computers, and the development of well-validated empirical molecular mechanical force fields. With these technologies, it has been demonstrated that simulation is not only capable of consistently reproducing experimental observations of sequence specific fine structure of DNA, but also can give detailed insight into prevalent problems in nucleic acid structure, ion association and specific hydration of nucleic acids, polyadenine tract bending, and the subtle environmental dependence of the A-DNA-B-DNA duplex equilibrium. Despite the advances, there are still issues with the methods that need to be resolved through rigorous controlled testing. In general, these relate to deficiencies of the underlying molecular mechanical potentials or applied methods (such as the imposition of true periodicity in Ewald simulations and the need for energy conservation), and significant limits in effective conformational sampling. In this perspective, we provide an overview of our experiences, provide some cautionary notes, and provide recommendations for further study in molecular dynamics simulation of nucleic acids.