Multiscale QM/MM molecular dynamics study on the first steps of guanine damage by free hydroxyl radicals in solution

Photophysics of a nucleic acid-protein crosslinking model strongly depends on solvation dynamics: an experimental and theoretical study

We present a combined experimental and theoretical study of the photophysics of 5-benzyluracil (5BU) in methanol, which is a model system for interactions between nucleic acids and proteins. A molecular dynamics study of 5BU in solution through efficient DFT-based hybrid ab initio potentials revealed a remarkable conformational flexibility - allowing the population of two main conformers - as well as specific solute-solvent interactions, which both appear as relevant factors for the observed 5BU optical absorption properties. The simulated absorption spectrum, calculated on such an ensemble, enabled a molecular interpretation of the experimental UV-Vis lowest energy band, which is also involved in the induced photo-reactivity upon irradiation. In particular, the first two excited states (mainly involving the uracil moiety) both contribute to the 5BU lowest energy absorption. Moreover, as a key finding, the nature and brightness of such electronic transitions are strongly influenced by 5BU conformation and the microsolvation of its heteroatoms.

Influence of hydrogen bonds on the reaction of guanine and hydroxyl radical: DFT calculations in C(H+)GC motif

A comprehensive theoretical investigation was performed to illuminate the influence of hydrogen bonds (H-bonds) on the obscure reaction of a hydroxyl radical (HO˙) and guanine (G) by selecting the building block of parallel triplex DNA, C(H+)GC, as the model. By mapping the energy profiles for addition and hydrogen abstraction reactions, the favorable pathway is predicted. The results reveal that in the C(H+)GC context, barrierless hydrogen abstraction from N2 of G leading to a neutral radical G(N2-H)˙ appears to become significant, but electrophilic attack by HO˙ on C8 of G resulting in 8-oxoG is the most thermodynamically favorable course. This shows a strong structural dependence due to the context constrained by the H-bond, which is dramatically different from the situation in unencumbered G. More interestingly, it proves that the stability order of resulting adduct radicals is not altered by H-bonding, but the activity for possible sites of the hydroxylation reaction changes. The significant influence of the H-bond on elementary reactions involved in the reaction is emphasized in the C(H+)GC context but is not restricted to the H-abstraction reaction. It is greatly anticipated that the present study could provide thoughtful insights into the vague hydroxyl radical-induced oxidation chemistry.

A Theoretical Study of Hydrogen Abstraction Reactions in Guanosine and Uridine

All practically possible hydrogen abstraction reactions for guanosine and uridine have been investigated through quantum chemical calculations of energy barriers and rate constants. This was done at the level of density functional theory (DFT) with the ωB97X-D functional and the 6-311++G(2df,2pd) Pople basis set. Transition state theory with the Eckart tunneling correction was used to calculate the rate constants. The results show that the reaction involving the hydrogen labelled C4' in the ribofuranose part has the largest rate constant for guanosine with the value 5.097×1010 L mol-1s-1 and the largest for uridine with the value 1.62×1010 L mol-1s-1. Based on the results for these two nucleosides, there is a noticeable similarity between the rate constants in the ribofuranose part of the molecule, even though they are bound to two entirely different nucleobases.

Towards the ionizing radiation induced bond dissociation mechanism in oxygen, water, guanine and DNA fragmentation: a density functional theory simulation

The radiation-induced damages in bio-molecules are ubiquitous processes in radiotherapy and radio-biology, and critical to space projects. In this study, we present a precise quantification of the fragmentation mechanisms of deoxyribonucleic acid (DNA) and the molecules surrounding DNA such as oxygen and water under non-equilibrium conditions using the first-principle calculations based on density functional theory (DFT). Our results reveal the structural stability of DNA bases and backbone that withstand up to a combined threshold of charge and hydrogen abstraction owing to simultaneously direct and indirect ionization processes. We show the hydrogen contents of the molecules significantly control the stability in the presence of radiation. This study provides comprehensive information on the impact of the direct and indirect induced bond dissociations and DNA damage and introduces a systematic methodology for fine-tuning the input parameters necessary for the large-scale Monte Carlo simulations of radio-biological responses and mitigation of detrimental effects of ionizing radiation.

Single charge exchange in collisions of energetic nuclei with biomolecules of interest to ion therapy

This study is on the theory of single electron capture by fast nuclei from a variety of molecular targets of biological significance with high relevance to ion therapy of deep-seated tumors. The adopted theoretical framework is that of the first principles of quantum physics. As such, no free, adjustable parameters are used. This is in sharp contrast to the associated existing cross section input data to Monte Carlo simulations that all abound with empirical/phenomenological formulae. The present theory has the well-established track of its predictive power. This means that the computed cross sections can confidently be used in the cases for which no experimental data exist. These cross sections are from the full continuum distorted wave method (CDW). We first compute atomic cross sections in the independent electron model and then generate the corresponding molecular cross sections. The latter follow from the former within the independent atom model accompanied by the Bragg additivity rule. The investigated atomic targets are from the backbone of DNA and/or RNA molecules. These are atomic hydrogen, carbon, nitrogen and oxygen (H, C, N, O). Neon is also added to this sequence of targets as an isoelectronic atomic counterpart of water vapor, methane and ammonia molecules. The studied molecular targets are H₂O (water vapor), CO (carbon-monoxide), CO₂ (carbon-dioxide), CH₄ (methane), C₂H₄ (ethylene), C₂H₆ (ethane), C₄H₁₀ (butane) as well as the DNA/RNA nucleobases C₄H₄N₂O₂ (uracil), C₅H₅N₅ (adenine), C₅H₅N₅O (guanine), C₅H₆N₂O₂ (thymine) and C₄H₅N₃O (cytosine). The obtained total cross sections for any electronic target shell are compared with the available experimental data and overall favorable agreement is recorded at intermediate and high impact energies, which is the validity domain of the CDW method.

Targeted DNA oxidation and trajectory of radical DNA using DFT based QM/MM dynamics

Molecular insight into electronic rearrangements and structural trajectories arising from oxidative damages to DNA backbone is of crucial importance in understanding the effect of ionizing radiation, developing DNA biosensors and designing effective DNA cleaving molecules. Employing a Density Functional Theory based multi-scale Quantum-Mechanical-Molecular-Mechanical (QM/MM) simulation and a suitable partitioning of the Hamiltonian on solvated nucleotide, and single-, and double-stranded DNA, we mimic hydrogen transfer reactions from the backbone by OH radicals and report structural trajectories arising from on-the-fly electronic charge- and spin-density redistribution in these three different structural topologies of DNA. Trajectories reveal that H4′ abstraction can disrupt the deoxyribose moiety through the formation of C4′=O4′ ketone and a π-bond with base at C1′-N9 in a nucleotide versus only partial ketone formation in single- and double-stranded DNA, where the orientation of the base is topologically restrained. However, H5′ abstraction can lead DNA cleavage at 5′ end through the formation of C5′=O5′ ketone and breakage of P-O5′ bond. Results demonstrate that structural damages from oxidative reactions are restrained by base stacking and base-pair hydrogen bonding. The methodology can be suitably used to study targeted DNA and RNA damages from radicals and radiomimetic drugs to design DNA cleaving molecules for chemotherapy.

Multi-level Quantum Mechanics and Molecular Mechanics Study of Ring Opening Process of Guanine Damage by Hydroxyl Radical in Aqueous Solution

Combining multi-level quantum mechanics theories and molecular mechanics with an explicit water model, we investigated the ring opening process of guanine damage by hydroxyl radical in aqueous solution. The detailed, atomic-level ring-opening mechanism along the reaction pathway was revealed in aqueous solution at the CCSD(T)/MM levels of theory. The potentials of mean force in aqueous solution were calculated at both the DFT/MM and CCSD(T)/MM levels of the theory. Our study found that the aqueous solution has a significant effect on this reaction in solution. In particular, by comparing the geometries of the stationary points between in gas phase and in aqueous solution, we found that the aqueous solution has a tremendous impact on the torsion angles much more than on the bond lengths and bending angles. Our calculated free-energy barrier height 31.6 kcal/mol at the CCSD(T)/MM level of theory agrees well with the one obtained based on gas-phase reaction profile and free energies of solvation. In addition, the reaction path in gas phase was also mapped using multi-level quantum mechanics theories, which shows a reaction barrier at 19.2 kcal/mol at the CCSD(T) level of theory, agreeing very well with a recent ab initio calculation result at 20.8 kcal/mol.

Free energy profiles for two ubiquitous damaging agents: methylation and hydroxylation of guanine in B-DNA

DNA methylation and hydroxylation are two ubiquitous reactions in DNA damage induction, yet insights are scarce concerning the free energy of activation within B-DNA. We resort to multiscale simulations to investigate the attack of a hydroxyl radical and of the primary diazonium onto a guanine embedded in a solvated dodecamer. Reaction free energy profiles characterize two strongly exergonic processes, yet allow unprecedented quantification of the barrier towards this damage reaction, not higher than 6 kcal mol^-1 and sometimes inexistent, and of the exergonicities. In the case of the [G(C8)-OH]˙ intermediate, we challenge the functional dependence of such simulations: recently-proposed functionals, such as M06-2X and LC-BLYP, agree on a ∼4 kcal mol^-1 barrier, whereas the hybrid GGA B3LYP functional predicts a barrier-less pathway. In the long term, multiscale approaches can help build up a unified panorama of DNA lesion induction. These results stress the importance of DFT/MM-MD simulations involving new functionals towards the sound modelling of biomolecule damage even in the ground state.

Understanding DNA under oxidative stress and sensitization: the role of molecular modeling

DNA is constantly exposed to damaging threats coming from oxidative stress, i.e., from the presence of free radicals and reactive oxygen species. Sensitization from exogenous and endogenous compounds that strongly enhance the frequency of light-induced lesions also plays an important role. The experimental determination of DNA lesions, though a difficult subject, is somehow well established and allows to elucidate even extremely rare DNA lesions. In parallel, molecular modeling has become fundamental to clearly understand the fine mechanisms related to DNA defects induction. Indeed, it offers an unprecedented possibility to get access to an atomistic or even electronic resolution. Ab initio molecular dynamics may also describe the time-evolution of the molecular system and its reactivity. Yet the modeling of DNA (photo-)reactions does necessitate elaborate multi-scale methodologies to tackle a damage induction reactivity that takes place in a complex environment. The double-stranded DNA environment is first characterized by a very high flexibility, but also a strongly inhomogeneous electrostatic embedding. Additionally, one aims at capturing more subtle effects, such as the sequence selectivity which is of critical important for DNA damage. The structure and dynamics of the DNA/sensitizers complexes, as well as the photo-induced electron- and energy-transfer phenomena taking place upon sensitization, should be carefully modeled. Finally the factors inducing different repair ratios for different lesions should also be rationalized. In this review we will critically analyze the different computational strategies used to model DNA lesions. A clear picture of the complex interplay between reactivity and structural factors will be sketched. The use of proper multi-scale modeling leads to the in-depth comprehension of DNA lesions mechanisms and also to the rational design of new chemo-therapeutic agents.

Influence of reactive species on the modification of biomolecules generated from the soft plasma

Plasma medicine is an upcoming research area that has attracted the scientists to explore more deeply the utility of plasma. So, apart from the treating biomaterials and tissues with plasma, we have studied the effect of soft plasma with different feeding gases such as Air, N2 and Ar on modification of biomolecules. Hence, in this work we have used the soft plasma on biomolecules such as proteins ((Hemoglobin (Hb) and Myoglobin (Mb)), calf thymus DNA and amino acids. The structural changes or structural modification of proteins and DNA have been studied using circular dichroism (CD), fluorescence spectroscopy, protein oxidation test, gel electrophoresis, UV-vis spectroscopy, dynamic light scattering (DLS) and 1D NMR, while Liquid Chromatograph/Capillary Electrophoresis-Mass Spectrometer (LC/CE-MS) based on qualitative and quantitative bio-analysis have been used to study the modification of amino acids. Further, the thermal analysis of the protein has been studied with differential scanning calorimetry (DSC) and CD. Additionally, we have performed docking studies of H2O2 with Hb and Mb, which reveals that H2O2 molecules preferably attack the amino acids near heme group. We have also shown that N2 gas plasma has strong deformation action on biomolecules and compared to other gases plasma.

Effect of nanosecond-pulsed plasma on the structural modification of biomolecules

Nanosecond-pulsed plasma (NPP) action on biomolecules modification.

Attacking mechanism of hydroxyl radical to DNA base-pair: density functional study in vacuum and in water

Recently, the influence of radiation on human body has been recognized as a serious problem. In particular, highly reactive hydroxyl radicals *OH produced by the radiation react with DNA, resulting in a great damage on its structure and electronic properties. It is thus important to investigate the reaction mechanism of *OH to DNA for elucidating the initial damage in DNA induced by the radiation. In the present study, we search for transition states (TS) of the reaction between G-C/A-T base-pair and [Formula: see text] in vacuum and in water, by the density functional theory (DFT) calculations. At first, we obtain the stable structures for the dehydrogenated G-C and A-T, in which the hydrogen atom of NH2 group of G or A base is abstracted by [Formula: see text]. From the structures of the dehydrogenated as well as the natural base-pairs, the TS between these structures is searched for and the activation free energy (AFE) is estimated for the reaction. In vacuum, AFEs for the G-C and A-T are almost the same each other, while the stabilization energy by the reaction for G-C is about 4.9 kcal/mol larger than that for A-T, indicating that the population of the dehydrogenated G-C is remarkably larger than that of the dehydrogenated A-T in vacuum. On the other hand, in water approximated by the continuum solvation model, the AFE for A-T is 2.6 kcal/mol smaller than that for G-C, indicating that the reaction dehydrogenated by [Formula: see text] occurs more frequently for the solvated A-T base-pair than G-C.

Molecular modeling of nucleic acid structure: energy and sampling

An overview of computer simulation techniques as applied to nucleic acid systems is presented. This unit discusses methods used to treat energy and to sample representative configurations. Emphasis is placed on molecular mechanics and empirical force fields.

Molecular modeling of nucleic acid structure: energy and sampling

An overview of computer simulation techniques as applied to nucleic acid systems is presented. This unit expands an accompanying overview unit (UNIT 7.5) by discussing methods used to treat the energy and sample representative configurations. Emphasis is placed on molecular mechanics and empirical force fields.