Nucleosome Histone Tail Conformation and Dynamics: Impacts of Lysine Acetylation and a Nearby Minor Groove Benzo[a]pyrene-Derived Lesion.

Histone tails in nucleosomes play critical roles in regulation of many biological processes, including chromatin compaction, transcription, and DNA repair. Moreover, post-translational modifications, notably lysine acetylation, are crucial to these functions. While the tails have been intensively studied, how the structures and dynamics of tails are impacted by the presence of a nearby bulky DNA lesion is a frontier research area, and how these properties are impacted by tail lysine acetylation remains unexplored. To obtain molecular insight, we have utilized all atom 3 µs molecular dynamics simulations of nucleosome core particles (NCPs) to determine the impact of a nearby DNA lesion, 10S (+)-trans-anti-B[a]P-N2-dG-the major adduct derived from the procarcinogen benzo[a]pyrene-on H2B tail behavior in unacetylated and acetylated states. We similarly studied lesion-free NCPs to investigate the normal properties of the H2B tail in both states. In the lesion-free NCPs, charge neutralization upon lysine acetylation causes release of the tail from the DNA. When the lesion is present, it stably engulfs part of the nearby tail, impairing the interactions between DNA and tail. With the tail in an acetylated state, the lesion still interacts with part of it, although unstably. The lesion's partial entrapment of the tail should hinder the tail from interacting with other nucleosomes, and other proteins such as acetylases, deacetylases, and acetyl-lysine binding proteins, and thus disrupt critical tail-governed processes. Hence, the lesion would impede tail functions modulated by acetylation or deacetylation, causing aberrant chromatin structures and impaired biological transactions such as transcription and DNA repair.

Entrapment of a Histone Tail by a DNA Lesion in a Nucleosome Suggests the Lesion Impacts Epigenetic Marking: A Molecular Dynamics Study.

Errors in epigenetic markings are associated with human diseases, including cancer. We have used molecular dynamics simulations of a nucleosome containing the 10S (+)-trans-anti-B[a]P-N(2)-dG lesion, derived from the environmental pro-carcinogen benzo[a]pyrene, to elucidate the impact of the lesion on the structure and dynamics of a nearby histone N-terminal tail. Our results show that a lysine-containing part of this H2B tail that is subject to post-translational modification is engulfed by the enlarged DNA minor groove imposed by the lesion. The tail entrapment suggests that epigenetic markings could be hampered by this lesion, thereby impacting critical cellular functions, including transcription and repair.

Structural basis for the recognition of diastereomeric 5',8-cyclo-2'-deoxypurine lesions by the human nucleotide excision repair system.

The hydroxyl radical is a powerful oxidant that generates DNA lesions including the stereoisomeric R and S 5',8-cyclo-2'-deoxyadenosine (cdA) and 5',8-cyclo-2'-deoxyguanosine (cdG) pairs that have been detected in cellular DNA. Unlike some other oxidatively generated DNA lesions, cdG and cdA are repaired by the human nucleotide excision repair (NER) apparatus. The relative NER efficiencies of all four cyclopurines were measured and compared in identical human HeLa cell extracts for the first time under identical conditions, using identical sequence contexts. The cdA and cdG lesions were excised with similar efficiencies, but the efficiencies for both 5'R cyclopurines were greater by a factor of ∼2 than for the 5'S lesions. Molecular modeling and dynamics simulations have revealed structural and energetic origins of this difference in NER-incision efficiencies. These lesions cause greater DNA backbone distortions and dynamics relative to unmodified DNA in 5'R than in 5'S stereoisomers, producing greater impairment in van der Waals stacking interaction energies in the 5'R cases. The locally impaired stacking interaction energies correlate with relative NER incision efficiencies, and explain these results on a structural basis in terms of differences in dynamic perturbations of the DNA backbone imposed by the R and S covalent 5',8 bonds.

Differences in the Access of Lesions to the Nucleotide Excision Repair Machinery in Nucleosomes.

In nucleosomes, the access of DNA lesions to nucleotide excision repair is hindered by histone proteins. However, evidence that the nature of the DNA lesions may play a role in facilitating access is emerging, but these phenomena are not well-understood. We have used molecular dynamics simulations to elucidate the structural, dynamic, and energetic properties of the R and S 5'-8-cyclo-2'-dG and the (+)-cis-anti-B[a]P-dG lesions in a nucleosome. Our results show that the (+)-cis-anti-B[a]P-dG adduct is more dynamic and more destabilizing than the smaller and more constrained 5',8-cyclo-2'-dG lesions, suggesting more facile access to the more bulky (+)-cis-anti-B[a]P-dG lesion.

Nuclear magnetic resonance studies of an N2-guanine adduct derived from the tumorigen dibenzo[a,l]pyrene in DNA: impact of adduct stereochemistry, size, and local DNA sequence on solution conformations.

The dimensions and arrangements of aromatic rings (topology) in adducts derived from the reactions of polycyclic aromatic hydrocarbon (PAH) diol epoxide metabolites with DNA influence the distortions and stabilities of double-stranded DNA, and hence their recognition and processing by the human nucleotide excision repair (NER) system. Dibenzo[a,l]pyrene (DB[a,l]P) is a highly tumorigenic six-ring PAH, which contains a nonplanar and aromatic fjord region that is absent in the structurally related bay region five-ring PAH benzo[a]pyrene (B[a]P). The PAH diol epoxide-DNA adducts formed include the stereoisomeric 14S and 14R trans-anti-DB[a,l]P-N(2)-dG and the stereochemically analogous 10S- and 10R-B[a]P-N(2)-dG (B[a]P-dG) guanine adducts. However, nuclear magnetic resonance (NMR) solution studies of the 14S-DB[a,l]P-N(2)-dG adduct in DNA have not yet been presented. Here we have investigated the 14S-DB[a,l]P-N(2)-dG adduct in two different sequence contexts using NMR methods with distance-restrained molecular dynamics simulations. In duplexes with dC opposite the adduct deleted, a well-resolved base-displaced intercalative adduct conformation can be observed. In full duplexes, in contrast to the intercalated 14R stereoisomeric adduct, the bulky DB[a,l]P residue in the 14S adduct is positioned in a greatly widened and distorted minor groove, with significant disruptions and distortions of base pairing at the lesion site and two 5'-side adjacent base pairs. These unique structural features are significantly different from those of the stereochemically analogous but smaller B[a]P-dG adduct. The greater size and different topology of the DB[a,l]P aromatic ring system lead to greater structurally destabilizing DNA distortions that are partially compensated by stabilizing DB[a,l]P-DNA van der Waals interactions, whose combined effects impact the NER response to the adduct. These structural results broaden our understanding of the structure-function relationship in NER.

Structural, energetic and dynamic properties of guanine(C8)-thymine(N3) cross-links in DNA provide insights on susceptibility to nucleotide excision repair.

The one-electron oxidation of guanine in DNA by carbonate radical anions, a decomposition product of peroxynitrosocarbonate which is associated with the inflammatory response, can lead to the formation of intrastrand cross-links between guanine and thymine bases [Crean et al. (Oxidation of single-stranded oligonucleotides by carbonate radical anions: generating intrastrand cross-links between guanine and thymine bases separated by cytosines. Nucleic Acids Res. 2008; 36: 742-755.)]. These involve covalent bonds between the C8 positions of guanine (G*) and N3 of thymine (T*) in 5'-d(…G*pT*…) and 5'-d(…G*pCpT*…) sequence contexts. We have performed nucleotide excision repair (NER) experiments in human HeLa cell extracts which show that the G*CT* intrastrand cross-link is excised with approximately four times greater efficiency than the G*T* cross-link embedded in 135-mer DNA duplexes. In addition, thermal melting studies reveal that both lesions significantly destabilize duplex DNA, and that the destabilization induced by the G*CT* cross-link is considerably greater. Consistent with this difference in NER, our computations show that both lesions dynamically distort and destabilize duplex DNA. They disturb Watson-Crick base-pairing and base-stacking interactions, and cause untwisting and minor groove opening. These structural perturbations are much more pronounced in the G*CT* than in the G*T* cross-link. Our combined experimental and computational studies provide structural and thermodynamic understanding of the features of the damaged duplexes that produce the most robust NER response.

Role of structural and energetic factors in regulating repair of a bulky DNA lesion with different opposite partner bases.

Extensive molecular modeling with molecular dynamics simulations and van der Waals energy analyses were used to elucidate the striking finding that a mutagenic benzo[a]pyrene-derived DNA lesion, the base-displaced intercalated 10R-(+)-cis-anti-B[a]P-N(2)-dG (G*), manifests large differences in nucleotide excision repair (NER) efficiencies in DNA duplexes, which depend on the identities of the partner base opposite G*. The nature of the partner base causes marked differences in the extent of its major groove extrusion and dynamics, as well as energetic stability of the intercalation pocket that parallels the relative NER efficiencies.

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.

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.

Resistance of bulky DNA lesions to nucleotide excision repair can result from extensive aromatic lesion-base stacking interactions.

The molecular basis of resistance to nucleotide excision repair (NER) of certain bulky DNA lesions is poorly understood. To address this issue, we have studied NER in human HeLa cell extracts of two topologically distinct lesions, one derived from benzo[a]pyrene (10R-(+)-cis-anti-B[a]P-N(2)-dG), and one from the food mutagen 2-amino-1-methyl-6-phenylimidazo[4,5-b]pyridine (C8-dG-PhIP), embedded in either full or 'deletion' duplexes (the partner nucleotide opposite the lesion is missing). All lesions adopt base-displaced intercalated conformations. Both full duplexes are thermodynamically destabilized and are excellent substrates of NER. However, the identical 10R-(+)-cis-anti-B[a]P-N(2)-dG adduct in the deletion duplex dramatically enhances the thermal stability of this duplex, and is completely resistant to NER. Molecular dynamics simulations show that B[a]P lesion-induced distortion/destabilization is compensated by stabilizing aromatic ring system-base stacking interactions. In the C8-dG-PhIP-deletion duplex, the smaller size of the aromatic ring system and the mobile phenyl ring are less stabilizing and yield moderate NER efficiency. Thus, a partner nucleotide opposite the lesion is not an absolute requirement for the successful initiation of NER. Our observations are consistent with the hypothesis that carcinogen-base stacking interactions, which contribute to the local DNA stability, can prevent the successful insertion of an XPC ß-hairpin into the duplex and the normal recruitment of other downstream NER factors.

A CT-Based Radiomics Nomogram Combined with Clinic-Radiological Characteristics for Preoperative Prediction of the Novel IASLC Grading of Invasive Pulmonary Adenocarcinoma.

RATIONALE AND OBJECTIVES: The novel International Association for the Study of Lung Cancer (IASLC) grading system of invasive lung adenocarcinoma (ADC) demonstrated a remarkable prognostic effect and enabled numerous patients to benefit from adjuvant chemotherapy. We sought to build a CT-based nomogram for preoperative prediction of the IASLC grading. MATERIALS AND METHODS: This work retrospectively analyzed the CT images and clinical data of 303 patients with pathologically confirmed invasive ADC. The histological subtypes and radiological characteristics of the patients were re-evaluated. Radiomics features were extracted, and the optimal subset of features was established by ANOVA, spearman correlation analysis, and the least absolute shrinkage and selection operator (LASSO). Univariate and multivariate analyses identified the independent clinical and radiological variables. Finally, multivariate logistic regression analysis incorporated clinical, radiological, and optimal radiomics features into the nomogram. Receiver operating characteristic (ROC) curve, and accuracy were applied to assess the model's performance. Decision curve analysis (DCA), and calibration curve were applied to assess the clinical usefulness. RESULTS: Nine selected CT image features were used to develop the radiomics model. The accuracy, precision, sensitivity, and specificity of the radiomics model outperformed the clinic-radiological model in the training and testing sets. Integrating Radscore with independent radiological characteristics showed higher prediction performance than clinic-radiological characteristics alone in the training (AUC, 0.915 vs. 0.882; DeLong, p < 0.05) and testing (AUC, 0.838 vs. 0.782; DeLong, p < 0.05) sets. Good calibration and decision curve analysis demonstrated the clinical usefulness of the nomogram. CONCLUSION: Radiomics features effectively predict high-grade ADC. The combined nomogram may facilitate selecting patients who benefit from adjuvant treatment.

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.

Nuclear magnetic resonance solution structure of an N(2)-guanine DNA adduct derived from the potent tumorigen dibenzo[a,l]pyrene: intercalation from the minor groove with ruptured Watson-Crick base pairing.

The most potent tumorigen identified among the polycyclic aromatic hydrocarbons (PAH) is the nonplanar fjord region dibenzo[a,l]pyrene (DB[a,l]P). It is metabolically activated in vivo through the widely studied diol epoxide (DE) pathway to form covalent adducts with DNA bases, predominantly guanine and adenine. The (+)-11S,12R,13R,14S DE enantiomer forms adducts via its C14 position with the exocyclic amino group of guanine. Here, we present the first nuclear magnetic resonance solution structure of a DB[a,l]P-derived adduct, the 14R-(+)-trans-anti-DB[a,l]P-N(2)-dG (DB[a,l]P-dG) lesion in double-stranded DNA. In contrast to the stereochemically identical benzo[a]pyrene-derived N(2)-dG adduct (B[a]P-dG) in which the B[a]P rings reside in the B-DNA minor groove on the 3'-side of the modifed deoxyguanosine, in the DB[a,l]P-derived adduct the DB[a,l]P rings intercalate into the duplex on the 3'-side of the modified base from the sterically crowded minor groove. Watson-Crick base pairing of the modified guanine with the partner cytosine is broken, but these bases retain some stacking with the bulky DB[a,l]P ring system. This new theme in PAH DE-DNA adduct conformation differs from (1) the classical intercalation motif in which Watson-Crick base pairing is intact at the lesion site and (2) the base-displaced intercalation motif in which the damaged base and its partner are extruded from the helix. The structural considerations that lead to the intercalated conformation of the DB[a,l]P-dG lesion in contrast to the minor groove alignment of the B[a]P-dG adduct, and the implications of the DB[a,l]P-dG conformational motif for the recognition of such DNA lesions by the human nucleotide excision repair apparatus, are discussed.

Extraction of Treatment Information From Electronic Health Records and Evaluation of Testosterone Recovery in Patients With Prostate Cancer.

PURPOSE: Data quality and standardization remain a challenge when analyzing real-world clinical data. We built a clinical research database, using machine learning and natural learning processing, and investigated factors influencing testosterone recovery (T-recovery) in patients with localized prostate cancer (LPC) after initial androgen deprivation therapy (ADT). METHODS: Medication and treatment-associated dates missing in structured tables were extracted from patient notes using ConceptMapper, an automated data extraction tool, standardized and curated in Sema4 clinical research database. ADT usage duration was evaluated, and T-recovery in patients with LPC was analyzed by the Kaplan-Meier method and multivariable Cox proportional hazards models. We assessed the prognostic value of post-ADT T-recovery with prostate-specific antigen progression-free survival and failure-free survival. RESULTS: In total, 4,125 of 30,832 (13.4%) patients with prostate cancer had medication exclusively from notes with high precision and recall, F1 score ≥ 0.95. Association of dates with medication usage had a F1 score of 0.76. ADT duration estimation had higher accuracy combining information from notes to tables from electronic medical record (70% v 45%). Baseline testosterone was the strongest predictor of T-recovery in these patients. Patients with a baseline testosterone ≥ 300 ng/dL recovered in 9.79 versus 38 months for patients with baseline testosterone < 300 ng/dL (P < .0001). Shorter prostate-specific antigen progression-free interval was observed for patients with T-recovery (≥ 300 ng/dL) at 6 months after ADT cessation compared with patients without T-recovery (< 300 ng/dL; 13.7 v 25.1 months; P = .055). CONCLUSION: We augmented structured electronic medical record data with data extracted from notes and improved the accuracy of medication information for patients. ADT exposure and T-recovery in patients with LPC produced results consistent with the literature and clinical experience and illustrates the power of applying machine learning methods to enhance the quality of real-world evidence in answering clinically relevant questions.

Intercalative conformations of the 14R (+)- and 14S (-)-trans-anti-DB[a,l]P-N6-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-N6-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-N6-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-N6-dA adducts. For the DB[a,l]P-N6-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-N6-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.

[Research on the arrangement of hospital beds and empirical analysis].

This paper was aimed to find an efficient plan for the arrangement of the beds in hospitals so that the bed-operating efficiency could be raised and the requirements of patients could be perfectly met. After analyzing the data using the method of mathematical statistics, we defined the concept of "Waste Index". In addition to the concept, we realized the auto-simulation of the model using the dynamic programming model, and design the program by MATLAB. We got the characteristics of bed-occupation conditions by data analysis. We obtained the result of optimizing model by programming and the most efficient plan which was obviously better than the traditional FCFS (First Come First Serve) plan. It can be well concluded that the mathematical model should be of great maneuverability in the practical situation to ameliorate the management efficiency of hospital beds considerably.

A bulky DNA lesion derived from a highly potent polycyclic aromatic tumorigen stabilizes nucleosome core particle structure.

The impact of a bulky DNA lesion on the structure and dynamics of a nucleosome core particle (NCP) containing a lesion derived from the unusually potent tumorigen dibenzo[a,l]pyrene that resists nucleotide excision repair (NER) in free DNA was investigated using 65 ns molecular dynamics simulations. Our results reveal that, relative to unmodified NCP, the lesion stabilizes the nucleosome via stacking interactions, improved Watson-Crick base pairing, hydrogen bonding between DNA and histones, and damped dynamics. These findings suggest that such lesions should be as resistant to NER in the nucleosome environment as they are in free DNA.

Base flipping free energy profiles for damaged and undamaged DNA.

Lesion-induced thermodynamic destabilization is believed to facilitate ß-hairpin intrusion by the human XPC/hHR23B nucleotide excision repair (NER) recognition factor, accompanied by partner-base flipping, as suggested by the crystal structure of the yeast orthologue (Min, J. H., and Pavletich, N. P. (2007) Nature 449, 570-575). To investigate this proposed mechanism, we employed the umbrella sampling method to compute partner base flipping free energies for the repair susceptible 14R (+)-trans-anti-DB[a,l]P-N(2)-dG modified duplex 11-mer, derived from the fjord region polycyclic aromatic hydrocarbon dibenzo[a,l]pyrene, and for the undamaged duplex. Our flipping free energy profiles show that the adduct has a lower flipping barrier by ∼7.7 kcal/mol, consistent with its thermally destabilizing impact on the damaged DNA duplex and its susceptibility to NER.

Variable impact of conformationally distinct DNA lesions on nucleosome structure and dynamics: Implications for nucleotide excision repair.

The packaging of DNA in nucleosomes presents a barrier for biological transactions including replication, transcription and repair. However, despite years of research, how the DNA is freed from the histone proteins and thereby allows the molecular machines to access the DNA remains poorly understood. We are interested in global genomic nucleotide excision repair (GG-NER). It is established that the histones are obstacles to this process, and DNA lesions are repaired less efficiently in nucleosomes than in free DNA. In the present study, we utilized molecular dynamics simulations to elucidate the nature of the distortions and dynamics imposed in the nucleosome by a set of three structually different lesions that vary in GG-NER efficiencies in free DNA, and in nucleosomes [Shafirovich, Geacintov, et. al, 2019]. Two of these are bulky lesions derived from metabolic activation of the environmental carcinogen benzo[a]pyrene, the 10R (+)-cis-anti-B[a]P-N2-dG and the stereoisomeric 10S (+)-trans-anti-B[a]P-N2-dG, which respectively adopt base-displaced/intercalated and minor groove-aligned conformations in DNA. The third is a non-bulky lesion, the 5'R-8-cyclo-2'-deoxyguanosine cross-link, produced by reactive oxygen and nitrogen species; cyclopurine lesions are highly mutagenic. These adducts are placed near the dyad axis, and rotationally with the lesion-containing strand facing towards or away from the histones. While each lesion has distinct conformational characteristics that are retained in the nucleosome, a spectrum of structural and dynamic disturbances, from slight to substantial, are displayed that depend on the lesion's structure and position in the nucleosome. We hypothesize that these intrinsic structural and dynamic distinctions provide different signals to initiate the cascade of chromatin-opening processes, including acetylation and other post translational modifications, remodeling by ATP-dependent complexes and spontaneous unwrapping that regulate the rate of access to the lesion; this may translate ultimately into varying GG-NER efficiencies, including repair resistance when signals for access are too weak.

Exocyclic amino groups of flanking guanines govern sequence-dependent adduct conformations and local structural distortions for minor groove-aligned benzo[a]pyrenyl-guanine lesions in a GG mutation hotspot context.

The environmental carcinogen benzo[a]pyrene (BP) is metabolized to reactive diol epoxides that bind to cellular DNA by predominantly forming N2-guanine adducts (G*). Mutation hotspots for these adducts are frequently found in 5'- ... GG ... dinucleotide sequences, but their origins are poorly understood. Here we used high resolution NMR and molecular dynamics simulations to investigate differences in G* adduct conformations in 5'- ... CG*GC ... and 5'- ... CGG* C... sequence contexts in otherwise identical 12-mer duplexes. The BP rings are positioned 5' along the modified strand in the minor groove in both cases. However, subtle orientational differences cause strong distinctions in structural distortions of the DNA duplexes, because the exocyclic amino groups of flanking guanines on both strands compete for space with the BP rings in the minor groove, acting as guideposts for placement of the BP. In the 5'- ... CGG* C ... case, the 5'-flanking G . C base pair is severely untwisted, concomitant with a bend deduced from electrophoretic mobility. In the 5'- ... CG*GC ... context, there is no untwisting, but there is significant destabilization of the 5'-flanking Watson-Crick base pair. The minor groove width opens near the lesion in both cases, but more for 5'- ... CGG*C.... Differential sequence-dependent removal rates of this lesion result and may contribute to the mutation hotspot phenomenon.