In the present study, the effect of 2H/1H isotopic exchange in hydrogen bonds between nitrogenous base pairs on occurrence and open states zones dynamics is investigated. These processes are studied using mathematical modeling, taking into account the number of open states between base pairs. The calculations of the probability of occurrence of open states in different parts of the gene were done depending on the localization of the deuterium atom. The mathematical modeling study demonstrated significant inequality (dependent on single 2H/1H replacement in DNA) among three parts of the gene similar in length of the frequency of occurrence of the open states. In this paper, the new convenient approach of the analysis of the abnormal frequency of open states in different parts of the gene encoding interferon alpha 17 was presented, which took into account both rising and decreasing of them that allowed to make a prediction of the functional instability of the specific DNA regions. One advantage of the new algorithm is diminishing the number of both false positive and false negative results in data filtered by this approach compared to the pure fractile methods, such as deciles or quartiles.

The formation and dynamics of the open states in a double-stranded DNA molecule are largely determined by its mechanical parameters. The main one is the torque. However, the experimental study of DNA dynamics and the occurrence of open states is limited by the spatial resolution of available biophysical instruments. Therefore, in this work, on the basis of a mechanical mathematical model of DNA, calculations of the torque effect on the process of occurrence and dynamics of open states were carried out for the interferon alpha 17 gene. It was shown that torsion action leads to the occurrence of rotational movements of nitrogenous bases. This influence is nonlinear, and an increase in the amplitude of the torsion action does not lead to an automatic increase in the amplitude of rotational movements and an increase in the zones’ open states. Calculations with a constant torsion moment demonstrate that open states zones are more often formed at the boundaries of the gen and in regions with a predominance of A–T pairs. It is shown, that for the occurrence of open states in the part of the gene that contains a small number of A–T pairs, a large amount of torque is required. When the torque is applied to a certain region of the gene, the probability of the formation of the open state depends on the content of A–T pairs in this region, the size of this region, and on the exposure time. For this mathematical model, open states zones can be closed when the torsion action stops. The simulation results showed that the values of the torsion moment required for the appearance of open states zones, in some cases, are close to experimentally measured (13–15 pN·nm).

The sensitivity of DNA to electromagnetic radiation in different ranges differs depending on various factors. The aim of this study was to examine the molecular dynamics of DNA under the influence of external periodic influences with different frequencies. In the present paper, within the framework of a mechanical model without simplifications, we investigated the effect of various frequencies of external periodic action in the range from 1011 s−1 to 108 s−1 on the dynamics of a DNA molecule. It was shown that under the influence of an external periodic force, a DNA molecule can perform oscillatory movements with a specific frequency characteristic of this molecule, which differs from the frequency of the external influence ω. It was found that the frequency of such specific vibrations of a DNA molecule depends on the sequence of nucleotides. Using the developed mathematical model describing the rotational motion of the nitrogenous bases around the sugar–phosphate chain, it is possible to calculate the frequency and amplitude of the oscillations of an individual DNA area. Such calculations can find application in the field of molecular nanotechnology.

The role of single protium-deuterium exchanges in formation of open-state bubbles is investigated using mathematical modeling methods. It is shown that the probability for formation of bubbles with a certain length (12 to 27 nucleotides) depends on localization of the deuterium atom in the DNA molecule and may be appreciably different from the general probability for occurrence of open states.

The effect of single substitutions of protium for deuterium in hydrogen bonds between pairs of nitrogenous bases on the open states occurrence probability at high critical breaking energies of these bonds has been studied. The study was carried out using numerical methods based on the angular mathematical model of DNA. The IFNA17 gene was divided into three approximately equal parts. A comparison of the open states occurrence probability in these parts of the gene was done. To improve the accuracy of the results, a special data processing algorithm was developed. The developed methods have shown their suitability for taking into account the occurrence of open states in the entire range of high critical energies. It has been established that single 2H/1H substitutions in certain nitrogenous bases can be a mechanism for maintaining the vital activity of IFNA17 under critical conditions. In general, the developed method of the mathematical modeling provide unprecedented insight into the DNA behavior under the highest critical energy range, which greatly expands scientific understanding of nucleobases interaction.