Search

Four uracil-containing intrastrand cross-links have been detected in human cells upon UV irradiation of 5-bromouracil-containing DNA, namely 5'-G[8-5]U-3', 5'-U[5-8]G-3', 5'-A[8-5]U-3', and 5'-A[2-5]U-3'. These lesions feature unique composition and connectivity compared with other intrastrand cross-links reported in the literature. For the first time, structural information obtained using molecular dynamics (MD) simulations reveal that all four lesions distort the DNA helix, which can involve an extrahelical location of the cross-link, changes in the helical interactions of the complementary nucleotides, or disruption of hydrogen bonding in the flanking base pairs up to two positions from the cross-linked site; however, the degree of distortion varies between the cross-links, being affected by the sequence, nucleobase-nucleobase connectivity, and the purine involved. Most importantly, the relative distortion of the damaged DNA provides the first structural explanation for the observed abundances of the four uracil-containing cross-links. Furthermore, the highly distorted conformations suggest that these lesions will likely have severe implications for DNA replication and repair processes in cells.

The magnitudes of select π–π interactions are studied with MP2, as well as the fragment molecular orbital (FMO) method, which evaluates the total energy as a sum of polarized monomer energies and interfragment interaction energies. The monomer polarization is of particular interest and is found to increase in magnitude with system size and the electrostatic contribution to the interaction. As a result, differences will exist in the interaction energy reported from FMO calculations, and those calculated using a traditional supermolecular approach.

Previous computational studies have examined pi-pi and pi(+)-pi stacking and T-shaped interactions in nucleobase-amino acid dimers, yet it is important to investigate how additional amino acids affect these interactions since simultaneous contacts often appear in nature. Therefore, this paper investigates the geometries and binding strengths of amino acid-nucleobase-amino acid trimers, which are compared to the corresponding nucleobase-amino acid dimer interactions. We concentrate on systems containing the natural nucleobase adenine or its (cationic) damaged counterpart, 3-methyladenine, and the aromatic amino acid histidine, in both the neutral and protonated forms. This choice of molecules provides information about pi-pi and pi(+)-pi stacking and T-shaped interactions in asymmetric, biologically relevant systems. We determined that both stacked and T-shaped interactions, as well as both pi-pi and pi(+)-pi interactions, exhibit geometric additivity. To investigate the energetic additivity in our trimers, the synergy (E(syn)) and the additivity (E(add)) energy were examined. E(add) reveals that it is important to consider the interaction between the two amino acids when examining the additivity of nucleobase-amino acid interactions. Additionally, E(syn) and E(add) indicate that pi(+)-pi interactions are quite different from pi-pi interactions. The magnitude of E(add) is generally less than 2 kJ mol(-1), which suggests that these interactions are additive. However, the interaction energy analysis does not provide information about the individual interactions in the trimers. Therefore, the quantum theory of atoms in molecules (QTAIM) was implemented. We find inconsistent conclusions from our QTAIM analysis and interaction energy evaluation. However, the magnitudes of the differences between the dimer and trimer critical point properties are extremely small and therefore may not be able to yield conclusive descriptions of differences (if any) between the dimer and trimer interactions. We hypothesize that, due to the limited number of investigations of this type, it is currently unclear how QTAIM can improve our understanding of pi-pi and pi(+)-pi dimers and trimers. Therefore, future work must systematically alter the pi-pi or pi(+)-pi system to definitively determine how the geometry, symmetry, and system size alter the QTAIM analysis, which can then be used to understand biologically relevant complexes.

Changes in the magnitude of (M06-2X/6-31+G(d,p)) π-π stacking and T-shaped (nucleobase-edge and amino acid-edge) interactions between (neutral or protonated) histidine (His) and adenine (A) dimers upon microsolvation with up to four discrete water molecules were determined. A variety of histidine-water interactions were considered including conventional (N-H···O, N···H-O, C-H···O) hydrogen bonding and nonconventional (X-H···π (neutral His) or lone-pair···π (protonated His)) contacts. Overall, the effects of discrete His-H(2)O interactions on the neutral histidine-adenine π-π contacts are negligible (3 kJ mol(-1) or 15%) regardless of the type of water binding, the number of water molecules bound, or the His-A dimer (stacked or (amino acid- or nucleobase-edge) T-shaped) configuration. This suggests that previously reported gas-phase binding strengths for a variety of neutral amino acid-nucleobase dimers are likely relevant for a wide variety of (microsolvated) environments. In contrast, the presence of water decreases the histidine-adenine π(+)-π interaction by up to 15 kJ mol(-1) (or 30%) for all water binding modes and orientations, as well as different stacked and T-shaped His(+)-A dimers. Regardless of the larger effect of discrete histidine-water interactions on the magnitude of the π(+)-π compared with π-π interactions, the π(+)-π binding strengths remain substantially larger than the corresponding π-π contacts. These findings emphasize the distinct nature of π(+)-π and π-π interactions and suggest that π(+)-π contacts can provide significant stabilization in biological systems relative to π-π contacts under many different environmental conditions.

The (gas-phase) MP2/6-31G*(0.25) π···π stacking interactions between the five natural bases and the aromatic amino acids calculated using (truncated) monomers composed of conjugated rings and/or (extended) monomers containing the biological backbone (either the protein backbone or deoxyribose sugar) were previously compared. Although preliminary energetic results indicated that the protein backbone strengthens, while the deoxyribose sugar either strengthens or weakens, the interaction calculated using truncated models, the reasons for these effects were unknown. The present work explains these observations by dissecting the interaction energy of the extended complexes into individual backbone···π and π···π components. Our calculations reveal that the total interaction energy of the extended complex can be predicted as a sum of the backbone···π and π···π components, which indicates that the biological backbone does not significantly affect the ring system through π-polarization. Instead, we find that the backbone can indirectly affect the magnitude of the π···π contribution by changing the relative ring orientations in extended dimers compared with truncated dimers. Furthermore, the strengths of the individual backbone···π contributions are determined to be significant (up to 18 kJ mol(-1)). Therefore, the origin of the energetic change upon model extension is found to result from a balance between an additional (attractive) backbone···π component and differences in the strength of the π···π interaction. In addition, to understand the effects of the biological backbone on the stacking interactions at DNA-protein interfaces in nature, we analyzed the stacking interactions found in select DNA-protein crystal structures, and verified that an additive approach can be used to examine the strength of these interactions in biological complexes. Interestingly, although the presence of attractive backbone···π contacts is qualitatively confirmed using the quantum theory of atoms in molecules (QTAIM), QTAIM electron density analysis is unable to quantitatively predict the additive relationship of these interactions. Most importantly, this work reveals that both the backbone···π and π···π components must be carefully considered to accurately determine the overall stability of DNA-protein assemblies.

Bulky DNA addition products (adducts) formed through attack at the C8 site of guanine can adopt the syn orientation about the glycosidic bond due to changes in conformational stability or hydrogen-bonding preferences directly arising from the bulky group. Indeed, the bulky substituent may improve the stability of (non-native) Hoogsteen pairs. Therefore, such adducts often result in mutations upon DNA replication. This work examines the hydrogen-bonded pairs between the Watson-Crick and Hoogsteen faces of the ortho or para C8-phenoxyl-2'-deoxyguanosine adduct and each natural (undamaged) nucleobase with the goal to clarify the conformational preference of this type of damage, as well as provide insight into the likelihood of subsequent mutation events. B3LYP/6-311+G(2df,p)//B3LYP/6-31G(d) hydrogen-bond strengths were determined using both nucleobase and nucleoside models for adduct pairs, as well as the corresponding complexes involving natural 2'-deoxyguanosine. In addition to the magnitude of the binding strengths, the R(C1'···C1') distances and ∠(N9C1'C1') angles, as well as the degree of propeller-twist and buckle distortions, were carefully compared to the values observed in natural DNA strands. Due to structural changes in the adduct monomer upon inclusion of the sugar moiety, the monomer deformation energy significantly affects the relative hydrogen-bond strengths calculated with the nucleobase and nucleoside models. Therefore, we recommend the use of at least a nucleoside model to accurately evaluate hydrogen-bond strengths of base pairs involving flexible, bulky nucleobase adducts. Our results also emphasize the importance of considering both the magnitude of the hydrogen-bond strength and the structure of the base pair when predicting the preferential binding patterns of nucleobases. Using our best models, we conclude that the Watson-Crick face of the ortho phenoxyl adduct forms significantly more stable complexes than the Hoogsteen face, which implies that the anti orientation of the damaged base will be favored by hydrogen bonding in DNA helices. Additionally, regardless of the hydrogen-bonding face involved, cytosine forms the most stable base pair with the ortho adduct, which implies that misincorporation due to this type of damage is unlikely. Similarly, cytosine is the preferred binding partner for the Watson-Crick face of the para adduct. However, Hoogsteen interactions with the para adduct are stronger than those with natural 2'-deoxyguanosine or the ortho adduct, and this form of damage binds with nearly equal stability to both cytosine and guanine in the Hoogsteen orientation. Therefore, the para adduct may adopt multiple orientations in DNA helices and potentially cause mutations by forming pairs with different natural bases. Models of oligonucleotide duplexes must be used in future work to further evaluate other factors (stacking, major groove contacts) that may influence the conformation and binding preference of these adducts in DNA helices.

Detailed (gas-phase) MP2/6-31G*(0.25) potential energy surface scans and CCSD(T) energy calculations at the complete basis set (CBS) limit were used to analyze the (face-to-face) stacking and (edge-to-face) T-shaped interactions between histidine (modeled as imidazole) and the DNA nucleobases. For the first time, a variety of relative monomer arrangements between both neutral and protonated histidine and the natural nucleobases were considered to determine the effects of charge on the optimum dimer geometry and binding strength. Our results reveal that protonation of histidine changes the preferred relative orientations of the monomers and propose that these geometric differences may be combined with experimental crystal structures to assess the protonation state of histidine in different environments. It is also found that protonation affects the nucleobase binding preference, as well as the magnitude of the stacking and T-shaped interactions. Indeed, the maximum possible stacking and T-shaped interactions involving the neutral histidine range between approximately 20 and 45 kJ mol(-1), while this range increases to 40-105 kJ mol(-1) upon protonation, which represents an up to 330% enhancement. Although an increase in the interaction energies upon protonation of histidine is expected, the present work provides a measure of the magnitude of this enhancement in the gas phase and reveals that the amplification is almost entirely due to larger electrostatic contributions. The relative strengthening of different classifications of dimers upon protonation leads to stronger T-shaped interactions than stacking energies for protonated histidine, while the stacking and T-shaped interactions involving neutral histidine are of comparable magnitude. Thus, there is a significant difference in the nature of the pi(cation)-pi interactions involving protonated histidine and the pi-pi interactions involving neutral histidine. The calculated strengths of the interactions studied in the present work suggest that both neutral and cationic histidine contacts will provide significant stabilization to DNA-protein complexes. Although solvation effects will decrease the magnitude of the reported interactions, our results are applicable to a variety of low-polarity, biologically-relevant environments such as nonpolar enzyme active sites. Therefore, our calculations suggest that these interactions may also be important for many biological processes. The proposed significance of these interactions is supported by the large number of histidine-nucleobase contacts that appear in experimental crystal structures. The highly accurate (MP2/6-31G*(0.25)) preferred structures and (CCSD(T)/CBS) binding strengths reported in the present work can be used as benchmarks to analyze the performance of existing, or to develop new, molecular mechanics force fields for use in large-scale molecular dynamics (MD) studies of DNA-protein complexes.

The ability of a dinucleoside monophosphate to mimic the conformation of B-DNA was tested using a combination of different phosphate models (anionic, neutral, counterion), environments (gas, water), and density functionals (B3LYP, MPWB1K, M06-2X) with the 6-31G(d,p) basis set. Three sequences (5'-GX(Py)-3', where X(Py) = T, U or (Br)U) were considered, which vary in the (natural or modified) 3' pyrimidine nucleobase (X(Py)). These bases were selected due to their presence in natural DNA, structural similarity to T and/or applications in radical-initiated anti-tumour therapies. The accuracy of each of the 54 model, method and sequence combinations was judged based on the ability to reproduce key structural features of natural B-DNA, including the stacked base-base orientation and important backbone torsion angles. B3LYP yields distorted or tilted relative base-base orientations, while many computational challenges were encountered for MPWB1K. Despite wide use in computational studies of DNA, the neutral (protonated) phosphate model could not consistently predict a stacked arrangement for all sequences. In contrast, stacked base-base arrangements were obtained for all sequences with M06-2X in conjunction with both the anionic and (sodium) counterion phosphate models. However, comparison of calculated and experimental backbone conformations reveals the charge-neutralized counterion phosphate model best mimics B-DNA. Structures optimized with implicit solvent (water) are comparable to gas-phase structures, which suggests similar results should be obtained in an environment of intermediate polarity. We recommend the use of either gas-phase or water M06-2X optimizations with the counterion phosphate model to study the structure and/or reactivity of natural or modified DNA with a dinucleoside monophosphate.

The pi-pi stacking (face-to-face) interactions between the five natural DNA or RNA nucleobases and the four aromatic amino acids were compared using three different types of dimers: (1) a truncated nucleoside (nucleobase) stacked with a truncated amino acid; (2) a truncated nucleoside (nucleobase) stacked with an extended amino acid; and (3) a nucleoside (extended nucleobase) stacked with a truncated amino acid. Systematic (MP2/6-31G*(0.25)) potential energy surface scans reveal important information about the effects of the deoxyribose sugar and protein backbone on the structure and binding energy between truncated nucleobase and amino acid models that are typically implemented in the literature. Most notably, electrostatic and steric interactions arising from the bulkiness of the biological backbones can change the preferred relative orientations of DNA and protein pi-systems. More importantly, the protein backbone can strengthen the stacking energy (by up to 10 kJ mol(-1)), while the deoxyribose moiety can strengthen or weaken the stacking interaction depending on the positioning of the amino acid relative to the sugar residue. These effects are likely due to additional interactions between the amino acid or nucleobase ring and the backbone in the extended monomer rather than significant changes in the properties of the biological pi-systems upon model extension. Since the present work reveals that all calculated DNA-protein stacking interactions are significant and approach the strength of other noncovalent interactions between biomolecules, both pi-pi and backbone-pi interactions must be considered when attempting to gain a complete picture of DNA-protein binding.