Development and implementation of a metaphase DNA model for ionizing radiation induced DNA damage calculation

Objective. To develop a metaphase chromosome model representing the complete genome of a human lymphocyte cell to support microscopic Monte Carlo (MMC) simulation-based radiation-induced DNA damage studies.Approach. We first employed coarse-grained polymer physics simulation to obtain a rod-shaped chromatid segment of 730 nm in diameter and 460 nm in height to match Hi-C data. We then voxelized the segment with a voxel size of 11 nm per side and connected the chromatid with 30 types of pre-constructed nucleosomes and 6 types of linker DNAs in base pair (bp) resolutions. Afterward, we piled different numbers of voxelized chromatid segments to create 23 pairs of chromosomes of 1-5μm long. Finally, we arranged the chromosomes at the cell metaphase plate of 5.5μm in radius to create the complete set of metaphase chromosomes. We implemented the model in gMicroMC simulation by denoting the DNA structure in a four-level hierarchical tree: nucleotide pairs, nucleosomes and linker DNAs, chromatid segments, and chromosomes. We applied the model to compute DNA damage under different radiation conditions and compared the results to those obtained with G0/G1 model and experimental measurements. We also performed uncertainty analysis for relevant simulation parameters.Main results. The chromatid segment was successfully voxelized and connected in bps resolution, containing 26.8 mega bps (Mbps) of DNA. With 466 segments, we obtained the metaphase chromosome containing 12.5 Gbps of DNA. Applying it to compute the radiation-induced DNA damage, the obtained results were self-consistent and agreed with experimental measurements. Through the parameter uncertainty study, we found that the DNA damage ratio between metaphase and G0/G1 phase models was not sensitive to the chemical simulation time. The damage was also not sensitive to the specific parameter settings in the polymer physics simulation, as long as the produced metaphase model followed a similar contact map distribution.Significance. Experimental data reveal that ionizing radiation induced DNA damage is cell cycle dependent. Yet, DNA chromosome models, except for the G0/G1 phase, are not available in the state-of-the-art MMC simulation. For the first time, we successfully built a metaphase chromosome model and implemented it into MMC simulation for radiation-induced DNA damage computation.

Modeling of DNA and protein organization levels with Cn3D software

The molecular structure of living organisms and the complex interactions amongst its components are the basis for the diversity observed at the macroscopic level. Proteins and nucleic acids are some of the major molecular components, and play a key role in several biological functions, such as those of development and evolution. This article presents an educational, bioinformatics-based process, designed to enhance a better understanding of the structure of the above molecules. In addition, by using selected protein molecules, it attempts to clarify protein organization levels and how these are related with the structural and functional diversity, which define the biodiversity of living organisms at the macroscopic level. In the framework of this project, molecular modeling has been performed using the Cn3D software, created by the US National Center for Biotechnology Information (NCBI). Cn3D is a user friendly application, which is easy for students to get familiar with quickly. Our suggested process may be easily enriched by a multitude of protein, nucleic acid, or other molecule structures, which are freely accessible at NCBI website. The described process has been implemented by students (n = 225) of the 5th high school grade, in two high schools in Greece. © 2016 by The International Union of Biochemistry and Molecular Biology, 45(2):126-129, 2017.

A twist opening model for DNA

The real mechanisms of several biological processes involving DNA are not yet understood. We discuss here some aspects of the initiation of transcription, in particular the formation of the open complex and the activation mechanism associated to enhancer binding proteins. Transcription activation seems to be governed by underlying dynamical mechanisms related to several distortions of the double chain structure: a dynamical approach on a mesoscopic description level could then allow a deeper understanding of this complex process. Starting from the Peyrard Bishop (PB) model, that considers only the hydrogen bond stretching of each base pair, we describe here an extended DNA model, proposed in [1], that allows a rather good representation of the double helix geometry and of its structural features by the introduction of angular variables related to the twist angle. Using a generalized multiple scale expansion for the case of vectorial lattices derived elsewhere [2], we derive analytically small amplitude approximate solutions of the model which are movable and spatially localized: we present here the results of this calculation and show how the special shape of the solutions is in good agreement with what can be expected for coupled angular radial distortions in the real molecule.