Electrostatic energy analysis of 8-oxoguanine DNA lesion-molecular dynamics study

One nanosecond molecular dynamics (MD) simulation was performed for two DNA segments each composed of 30 base pairs. In one DNA segment the native guanines at nucleotides positions 17 and 19 were replaced with two 8-oxoguanines (8-oxoG) (8-oxoG is mutagenic DNA oxo-lesion). The analysis of results was focused on the electrostatic energy that is supposed to be significant factor causing the disruption of DNA base stacking in DNA duplex and may also serve as a signal toward the repair enzyme informing the presence of the lesion. The repulsive interaction between 8-oxoG and the entire DNA molecule was observed, which caused the extrahelical position of 8-oxoG (position 19). The repulsive electrostatic interaction between both 8-oxoG lesions contributed to the flipping out of one 8-oxoG and to the local instability of the lesioned DNA region. The electrostatic potential at the surface of DNA close to the lesions has more negative value than the same region on the native DNA. This electrostatic potential may signal presence of the lesion to the repair enzyme. In the simulation of native DNA segment, no significant structural changes were observed and B-DNA structure was well preserved throughout the MD simulation.     

Molecular dynamics simulation of clustered DNA damage sites containing 8-oxoguanine and abasic site

Clustered DNA damage sites induced by ionizing radiation have been suggested to have serious consequences to organisms, such as cancer, due to their reduced probability to be repaired by the enzymatic repair machinery of the cell. Although experimental results have revealed that clustered DNA damage sites effectively retard the efficient function of repair enzymes, it remains unclear as to what particular factors influence this retardation. In this study, approaches based on molecular dynamics (MD) simulation have been applied to examine conformational changes and energetic properties of DNA molecules containing clustered damage sites consisting of two lesioned sites, namely 7,8-dihydro-8-oxoguanine (8-oxoG) and apurinic/apyrimidinic (AP) site, located within a few base pairs of each other. After 1 ns of MD simulation, one of the six DNA molecules containing a clustered damage site develops specific characteristic features: sharp bending at the lesioned site and weakening or complete loss of electrostatic interaction energy between 8-oxoG and bases located on the complementary strand. From these results it is suggested that these changes would make it difficult for the repair enzyme to bind to the lesions within the clustered damage site and thereby result in a reduction of its repair capacity.     

8-oxoguanine lesioned B-DNA molecule complexed with repair enzyme hOGG1: a molecular dynamics study

The molecular dynamics (MD) simulation of DNA mutagenic oxidative lesion, 7,8-dihydro-8-oxoguanine (8-oxoG), complexed with the repair enzyme, human oxoguanine glycosylase 1 (hOGG1), was performed for 1 nanosecond (ns) in order to describe the dynamical process of DNA-enzyme complex formation. After 900 picoseconds of MD the lesioned DNA and enzyme formed a complex that lasted until the end of the simulation at 1 ns. The complex was mainly represented by the overlapping van der Waals surfaces of DNA and enzyme molecules. The amino group of arginine 324 was located close to the phosphodiester bond of the nucleotide with 8-oxoG enabling chemical reactions between amino acid and lesion. The broken hydrogen bonds resulting in locally collapsed B-DNA structure were observed at the lesion site. The phosphodiester bond at C5' of 8-oxoG was displaced to the position close to the amino group of arginine 324. The water-mediated hydrogen bond network was formed in each contact area between DNA and enzyme, further enhancing the stability of the complex. In the background simulation of the identical molecular system with the native DNA, neither the complex nor the water- mediated hydrogen bond network was observed.     

8-oxoguanine and 8-oxodeoxyguanosine Biomarkers of Oxidative DNA Damage: A Review on HPLC–ECD Determination

Reactive oxygen species (ROS) are continuously produced in living cells due to metabolic and biochemical reactions and due to exposure to physical, chemical and biological agents. Excessive ROS cause oxidative stress and lead to oxidative DNA damage. Within ROS-mediated DNA lesions, 8-oxoguanine (8-oxoG) and its nucleotide 8-oxo-2′-deoxyguanosine (8-oxodG)—the guanine and deoxyguanosine oxidation products, respectively, are regarded as the most significant biomarkers for oxidative DNA damage. The quantification of 8-oxoG and 8-oxodG in urine, blood, tissue and saliva is essential, being employed to determine the overall effects of oxidative stress and to assess the risk, diagnose, and evaluate the treatment of autoimmune, inflammatory, neurodegenerative and cardiovascular diseases, diabetes, cancer and other age-related diseases. High-performance liquid chromatography with electrochemical detection (HPLC–ECD) is largely employed for 8-oxoG and 8-oxodG determination in biological samples due to its high selectivity and sensitivity, down to the femtomolar range. This review seeks to provide an exhaustive analysis of the most recent reports on the HPLC–ECD determination of 8-oxoG and 8-oxodG in cellular DNA and body fluids, which is relevant for health research.     

Differences in electrostatic potential around DNA fragments containing guanine and 8-oxo-guanine

Changes of electrostatic potential (EP) around the DNA molecule resulting from chemical modifications of nucleotides may play a role in enzymatic recognition of damaged sites. Effects of chemical modifications of nucleotides on the structure of DNA have been characterized through large-scale density functional theory computations. Quantum mechanical structural optimizations of DNA fragments with three pairs of nucleotides and accompanying counteractions were performed with a B3LYP exchange-correlation functional and 6–31G** basis sets. The “intact” DNA fragment contained guanine in the middle layer, while the “damaged” fragment had the guanine replaced with 8-oxo-guanine. The electrostatic potential around these DNA fragments was projected on a surface around the double helix. The 2D maps of EP of intact and damaged DNA fragments were analyzed to identify these modifications of EP that result from the occurrence of 8-oxo-guanine. It was found that distortions of the phosphate groups and displacements of the accompanying countercations are clearly reflected in the EP maps. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Modelling CreA protein–DNA recognition determinants. A molecular dynamics study of fully charged CreA–DNA model in water

CreA, the negative regulator mediating carbon catabolism repression in Aspergillus nidulans, is a protein that contains a DNA-binding domain comprising two zinc finger motifs. A 3D model for the CreA–G4 (5′-GCGGGGGCGT-3′) complex is constructed on the basis of the structure of the Zif268–DNA crystal complex [Science 252 (1991) 809] and using similarity analysis and computer assisted modelling techniques. The CreA–G4 model was then subjected to a set of molecular dynamics (MD) studies. Based on our previous nano second long Zif268–DNA MD simulation, a 170 ps long trajectory was deemed sufficient to test possible DNA–protein interactions. A screening of the static model and of the trajectory was performed for protein amino acids, nucleotide bases, phosphates backbone and water molecules mediating protein–DNA contacts. Energy, root mean square deviation (RMSD), principal inertial moment and distances were analysed. For this time span, the stability shown in fluctuation patterns reveals the presence of a complex behaving in a manner similar to Zif268–DNA.

An unambiguous characterisation of the amino acids involved in DNA-binding was obtained. These results could contribute towards the establishment of a code of protein–DNA recognition for this class of DNA-binding motifs.     

Reactions of the OOH radical with guanine: Mechanisms of formation of 8-oxoguanine and other products

The mutagenic product 8-oxoguanine (8-oxoGua) is formed due to intermediacy of peroxyl (OOR) radicals in lipid peroxidation and protein oxidation-induced DNA damage. The mechanisms of these reactions are not yet understood properly. Therefore, in the present study, the mechanisms of formation of 8-oxoGua and other related products due to the reaction of the guanine base of DNA with the hydroperoxyl radical (OOH) were investigated theoretically employing the B3LYP and BHandHLYP hybrid functionals of density functional theory and the polarizable continuum model for solvation. It is found that the reaction of the OOH radical with guanine can occur following seven different mechanisms leading to the formation of various products including 8-oxoGua, its radicals, 5-hydroxy-8-oxoguanine and CO₂. The mechanism that yields 8-oxoGua as an intermediate and 5-hydroxy-8-oxoGua as the final product was found to be energetically most favorable.     

毛细管电泳法高效筛选8-氧代鸟嘌呤DNA糖基化酶的核酸适配体

8-氧代鸟嘌呤DNA糖基化酶(OGG1)是人体中重要的功能蛋白,在修复DNA氧化性损伤过程中起关键作用。氧化应激等引起的氧化损伤易导致炎症反应的发生,对OGG1的抑制可以一定程度上起到缓解作用;对癌细胞OGG1的抑制有望作为癌症治疗的新方法。目前的研究多集中于小分子对OGG1功能的影响和调控,而OGG1的适配体筛选尚未见报道。作为功能配体,适配体具有合成简单、高亲和力及高特异性等优点。该文筛选了OGG1的核酸适配体,结合毛细管电泳高效快速的优点建立了两种基于毛细管电泳-指数富集进化(CE-SELEX)技术的筛选方法:同步竞争法和多轮筛选法。同步竞争法利用单链结合蛋白(SSB)与核酸库中单链核酸的强结合能力,与目标蛋白OGG1组成竞争体系,并通过增加SSB浓度来增加竞争筛选压力,以去除与OGG1弱结合的核酸序列,一步筛选即可获得与OGG1强结合的核酸序列。多轮筛选法在相同孵育条件和电泳条件下,经3轮筛选获得OGG1的核酸适配体。比较两种筛选方法的筛选结果,筛选结果中频次最高的3条候选核酸适配体序列一致,其解离常数(K_D)值在1.71~2.64 μmol/L之间。分子对接分析结果表明候选适配体1(Apt 1)可能与OGG1中具有修复氧化性损伤功能的活性口袋结合。通过对两种筛选方法的对比,证明同步竞争法更加快速高效,对其他蛋白核酸适配体筛选方法的选择具有一定的指导意义。得到的适配体有望用于OGG1功能调控,以抑制其修复功能。     

Efficient calculation of many stacking and pairing free energies in DNA from a few molecular dynamics simulations

Through the use of the one-step perturbation approach, 130 free energies of base stacking and 1024 free energies of base pairing in DNA have been calculated from only five simulations of a nonphysical reference state. From analysis of a diverse set of 23 natural and unnatural bases, it appears that stacking free energies and stacking conformations play an important role in pairing of DNA nucleotides. On the one hand, favourable pairing free energies were found for bases that do not have the possibility to form canonical hydrogen bonds, while on the other hand, good hydrogen-bonding possibilities do not guarantee a favourable pairing free energy if the stacking of the bases dictates an unfavourable conformation. In this application, the one-step perturbation approach yields a wealth of both energetic and structural information at minimal computational cost.     

In silico study directed towards identification of the key structural features of GyrB inhibitors targeting MTB DNA gyrase: HQSAR,CoMSIA and molecular dynamics simulations

ABSTRACT

Mycobacterium tuberculosis DNA gyrase subunit B (GyrB) has been identified as a promising target for rational drug design against fluoroquinolone drug-resistant tuberculosis. In this study, we attempted to identify the key structural feature for highly potent GyrB inhibitors through 2D-QSAR using HQSAR, 3D-QSAR using CoMSIA and molecular dynamics (MD) simulations approaches on a series of thiazole urea core derivatives. The best HQSAR and CoMSIA models based on IC₅₀ and MIC displayed the structural basis required for good activity against both GyrB enzyme and mycobacterial cell. MD simulations and binding free energy analysis using MM-GBSA and waterswap calculations revealed that the urea core of inhibitors has the strongest interaction with Asp79 via hydrogen bond interactions. In addition, cation-pi interaction and hydrophobic interactions of the R₂ substituent with Arg82 and Arg141 help to enhance the binding affinity in the GyrB ATPase binding site. Thus, the present study provides crucial structural features and a structural concept for rational design of novel DNA gyrase inhibitors with improved biological activities against both enzyme and mycobacterial cell, and with good pharmacokinetic properties and drug safety profiles.     

Twist‐dependent stacking energy of base‐pair steps in B‐DNA geometry: A density functional theory approach

Stacking energy of all the 10 unique DNA base‐pair steps (bp step) are calculated using density functional theory within the ultrasoft pseudopotential plane wave method and local density approximation for the exchange‐correlation functional. We have studied the dependence of stacking energy on twist angle, an aspect found difficult to explain using classical theory. We have found that the twist angle for different bp steps at stacking energy minimum matches extremely well with the values of average twist obtained from B‐DNA crystal structure data. This indicates that the use of a proper quantum chemical method to calculate the π‐π electronic interactions may explain stacking energy without incorporating hydrophobic interaction through solvent or effect of backbone through pseudobond. From the twist angle‐dependent stacking energy profile, we have also generated the probability distributions of twist for all the bp steps and calculated the variance of the distribution. Our calculated variances show similar trend to that of the experimental data for which sufficient numbers of data are available. The TA, AT, and CG doublets show large variances among the 10 possible bp steps, indicating their maximum flexibility. This might be the case of unusual deformation observed at the TATA‐box while binding to TBP protein. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

DNA strand break: structural and electrostatic properties studied by molecular dynamics simulation

Due to their lethal consequences and a relatively high probability of introduction of repair errors and mutations, single and double strand breaks are among the most important and dangerous DNA lesions. However, the mechanisms of their recognition and repair processes are only poorly known at present. This work defines and analyzes a DNA with single strand break as a template study for future complex analyses of biologically serious double strand break damage and its enzymatic repair mechanisms. Besides a non-damaged DNA serving as a reference system with no surprising results, system with open valences of the atoms at the strand break ends as well as a system with filled valences were simulated. In both cases during the first few nanoseconds the broken ends of strand breaks are significantly exposed to the outside of the molecule. However, with increasing time, the system with single strand break with open valences is partially disrupted. On the contrary, the system with filled valences shows stable conformation with newly created hydrogen bond between the two strand break endings. Moreover, these endings are steadily situated in the inner part of the molecule, thus making the recognition and docking process of a repair enzyme more complicated in the case of filled valences.     

7,8-二氢-8-氧鸟嘌呤碱基对的量子化学研究

遗传信息的完整性不断受氧化基因的威胁,7,8-二氢-8-氧鸟嘌呤（8-oxo-G）是氧化DNA损伤最常见的产物. 氧化碱基会引起基因突变、癌变及衰老等. 应用量子化学方法分析得出：鸟嘌呤（G）被氧化为8-oxo-G后,其电荷分布、氢键的供体和受体位点的数目和位置随之改变,N7和O6原子所带的电荷变得更负,使得它们作为氢键供体的能力增强. 从而G被误认为其他碱基,与正常碱基形成多种氢键复合物. 可将8-oxo-G划分为3个作用位点与正常碱基相互作用. 与正常的单体相比,碱基对中形成氢键的受体原子上所带电荷平均变负0.05e,占原电荷的8%; 供体H原子所带电荷平均变正0.02e,占原电荷的4%. 1位点与正常碱基作用形成的氢键复合物更稳定,2位点和3位点性质相似,水溶剂使碱基对的结合能力减弱,其中与C作用形成氢键复合物的结合能减弱程度最大,且使碱基对结合能力的次序改变. 在8-oxo-G导致的GC→TA突变中,亲核反应位点从G所在链转到A（C）所在链,影响酶对碱基的识别,从而产生基因突变.     

Electrostatic complex of didodecyldimethylammonium and DNA: Langmuir monolayer, Langmuir–Blodgett film and dye recognition at air/water interface

DNA–didodecyldimethylammonium (DNA–DDDA) electrostatic complex was prepared and characterized through Fourier transformation infrared (FT-IR), ¹H NMR and circular dichroism (CD) spectroscopy. When the dye molecule aqueous solutions were used as the subphase, the interaction between three dye molecules, acridine orange (AO), ethidium bromide (EB) and 5,10,15,20-tetrakis(4-N-methylpyridyl)porphine tetra(p-toluenesulfonate) (TMPyP) and the complex at air/solution interface were investigated through the surface pressure–area (π–A) isotherms, Brewster angle microscopy and UV-Vis spectroscopy, respectively. Our investigation indicates that the interaction capabilities of the three dyes to DNA–DDDA complex are different and present an order of TMPyP>AO>EB. For the interaction forms, we believe that TMPyP intercalates into the double helix of DNA, and AO adsorbs onto the surface of the DNA. As for EB, the measured signal is too weak to give a definite interaction form in the present experiment.     

A molecular dynamics study on permeability of gases through parylene AF8 membranes

The amorphous and crystalline poly‐p‐xylylene (parylene) AF8 membranes are constructed by a novel computational technique, that is, a combined method of NVT + NPT‐molecular dynamics (MD) and gradually reduce the size. The related permeation properties are calculated using Grand Canonical Monte Carlo (GCMC) and NVT‐MD methods. The results show that amorphous and crystalline parylene AF8 membranes have different permeation properties. Compared with amorphous parylene AF8 area, crystalline parylene AF8 membranes provide less walking paths of gases. In parylene AF8 membranes, gases walk front and back among the stable sites for holding the minimal energy. Furthermore, gases walk just along the y axis of crystalline cell, however, randomly in amorphous area. The corresponding permeability coefficients approach the experimental data. Copyright © 2011 John Wiley & Sons, Ltd.     

Denaturation of Drew-Dickerson DNA in a high salt concentration medium: molecular dynamics simulations

We have performed molecular dynamics simulation on B-DNA duplex (CGCGAATTGCGC) at different temperatures. The DNA was immerged in a salt-water medium with 1 M NaCl concentration to investigate salt effect on the denaturation process. At each temperature, configurational entropy is estimated using the covariance matrix of atom-positional fluctuations, from which the melting temperature (T(m)) was found to be 349 K. The calculated configuration entropy for different bases shows that the melting process involves more peeling (including fraying from the ends) conformations, and therefore the untwisting of the duplex and peeling states form the transition state of the denaturation process. There is a narrow minor groove in the AATT sequence that becomes wider by increasing temperature which disappears at high temperatures, especially above the melting temperature. We have also calculated the fraction of denatured base pairs, f-curve, from which T(m) was found to be 340 K, close to experimental value of 341 K. We found that DNA at high salt concentrations has few hydrogen bonds even at temperatures higher than the T(m). Our calculations show the fact that adding salt leads to increase of T(m) and stabilization of DNA.     

Structure,thermodynamics and diffusion in asymmetric binary mixtures: a molecular dynamics simulation study

Structural and thermodynamic properties as well as diffusion coefficients of binary fluid mixtures with asymmetry in mass, size, charge and their combinations have been studied using classical molecular dynamics simulations. The fluid mixture is modelled as spherical particles interacting via the Weeks–Chandler–Andersen and Coulomb potential. The diameter, charge and mass of the fluid particles are in the range 6–60 Å, 1–10e and 1—500 amu, respectively. Systematic variations in pair-correlation functions, thermodynamic properties as well as the self-diffusion coefficient are found with the size, charge and mass ratio of the particles. The self-diffusion coefficient for systems having more than one type of asymmetry is calculated and expressed in terms of diffusion coefficients of systems with only one type of asymmetry.     

LET dependence of the formation of oxidative damage 8-hydroxy-2′-deoxyguanosine (8-OHdG) in 2′-deoxyguanosine aqueous solution irradiated with heavy ions

To evaluate the contribution of indirect action mediated by OH radicals in biological effects of high LET radiations, we examined the production of 8-hydroxy-2′-deoxyguanosine (8-OHdG) in 2′-deoxyguanosine aqueous solution using various ion species with an LET range from 20 to 300 keV/μm. The 8-OHdG yield decreased with increasing LET. In the hypoxic irradiation condition, the yield showed constant or rather increasing tendency above about 100 keV/μm, which is consistent with an oxygen-in-the-track hypothesis to explain the diminishment of oxygen effect.     

Molecular dynamics study of DNA binding by INT‐DBD under a polarized force field

The DNA binding domain of transposon Tn916 integrase (INT‐DBD) binds to DNA target site by positioning the face of a three‐stranded antiparallel β‐sheet within the major groove. As the negatively charged DNA directly interacts with the positively charged residues (such as Arg and Lys) of INT‐DBD, the electrostatic interaction is expected to play an important role in the dynamical stability of the protein–DNA binding complex. In the current work, the combined use of quantum‐based polarized protein‐specific charge (PPC) for protein and polarized nucleic acid‐specific charge (PNC) for DNA were employed in molecular dynamics simulation to study the interaction dynamics between INT‐DBD and DNA. Our study shows that the protein–DNA structure is stabilized by polarization and the calculated protein–DNA binding free energy is in good agreement with the experimental data. Furthermore, our study revealed a positive correlation between the measured binding energy difference in alanine mutation and the occupancy of the corresponding residue's hydrogen bond. This correlation relation directly relates the contribution of a specific residue to protein–DNA binding energy to the strength of the hydrogen bond formed between the specific residue and DNA. © 2013 Wiley Periodicals, Inc.     

The Electrostatic Origin of Low‐Hydration Polymorphism in DNA

In recent years, significant progress has been made towards uncovering the physical mechanisms of low‐hydration polymorphism in double‐helical DNA. The effect appears to be mechanistically similar in different biological systems, and it is due to the ability of water to form spanning H‐bonded networks around biomacromolecules via a quasi‐two‐dimensional percolation transition. In the case of DNA, disintegration of the spanning H‐bonded network leads to electrostatic condensation of DNA strands because, below the percolation threshold, water loses its high dielectric permittivity, whereas the concentration of neutralizing counterions becomes high. In this Concept article arguments propose that this simple electrostatic mechanism represents the universal origin of low‐hydration polymorphism in DNA.