Hiroaki Nakamura, Hisanori Miyanishi, Takuo Yasunaga, Susumu Fujiwara, Tomoko Mizuguchi, Ayako Nakata, Tsuyoshi Miyazaki, Takao Otsuka, Takahiro Kenmotsu, Yuji Hatano, and Shinji Saito
Using molecular dynamics (MD) simulation, we simulate the structural change of a telomeric DNA by β-decay of substituted tritium to helium-3. The configuration of the telomeric DNA is obtained by removing TRF2 protein from the TRF2-Dbd-DNA complex (Protein Data Bank ID is 3SJM). We assume that hydrogens (H) of guanines in the telomeric DNA are replaced to helium-3. Since this replacement of the H atoms to the 3He atoms changes the charge distribution significantly, the charge distribution used in the MD simulation for the modified guanine is obtained by the density functional theory calculations. We adopt, as the MD simulation, nanoscale molecular dynamics code with CHARMM36 force field using Langevin thermostat and Nosé–Hoover Langevin piston to control the temperature and pressure of the system, respectively. Moreover, changing both the number of replaced guanine N and the temperature of the system T, we calculate the root mean square deviation RMSD to quantify the dependence of the durability of the telomeric DNA on the β-decays. From the MD simulation, it is found that as N or T becomes larger, the RMSD of the DNA becomes also larger. Namely, it denotes that as the intensity of the β-decays becomes larger or as the temperature is increased, the DNA structure becomes more fragile. [ABSTRACT FROM AUTHOR]