## Abstract

The DNA molecule is the most evolved and most complex molecule created by nature. The primary role of DNA in medicine is long-term storage of genetic information. Genetic modifying is one of the most critical challenges that scientists face. On the other hand, it is said that under the influence of acoustic, electromagnetic, and scalar waves, the genetic code of DNA can be read or rewritten. In this article, the most accurate and comprehensive dynamic model will be presented for DNA. Each of the two strands is modeled with an out of plane curved beam and then by doubling this two strands with springs, consider the hydrogen bond strength between this two strands. Beams are traditionally descriptions of mechanical engineering structural elements or building. However, any structure such as automotive automobile frames, aircraft components, machine frames, and other mechanical or structural systems contain beam structures that are designed to carry lateral loads are analyzed similarly. Also, in this model, the mass of the nucleobases in the DNA structure, the effects of the fluid surrounding the DNA (nucleoplasm) and the effects of temperature changes are also considered. Finally, by deriving governing equations from Hamilton’s principle method and solving these equations with the generalized differential quadrature method (GDQM), the frequency and mode shape of the DNA is obtained for the first time. In the end, validation of the obtained results from solving the governing equations of mathematical model compared to the obtained results from the COMSOL software is confirmed. By the help of these results, a conceptual idea for controlling cancer with using the DNA resonance frequency is presented. This idea will be presented to stop the cancerous cell’s protein synthesis and modifying DNA sequence and genetic manipulation of the cell. On the other hand, by the presented DNA model and by obtaining DNA frequency, experimental studies of the effects of waves on DNA such as phantom effect or DNA teleportation can also be studied scientifically and precisely.

## Introduction

Deoxyribonucleic acid, more commonly known as DNA, is an evolved molecule that contains the genetic code of organisms. Every living thing has DNA within their cells. It is important for inheritance, coding for proteins and the genetic instruction guide for life and its processes. DNA holds the instructions for an organism or each cell development processes, reproduction and ultimately death. Over the past decades, empirical discussions have been proposed to modify the genes in the DNA. These changes have been much discussed in the medical field by drawing on applications such as treatment, preventing the development of cancer, or erupting an organ (for example, a tooth). On the other hand, it should be noted that only 3% of DNA capacity is considered in medical fields. In the last two decades, a topic called “wave genome” has been raised by Russian scientists, which states that 97% of other DNA is not only inapplicable but also has a more significant role; because DNA can be affected by acoustic, electromagnetic, and scalar waves. Under the influence of these waves, the genetic code can be read or rewritten. Another claim of Russian scientists is that DNA is a biological network that binds all humans. About impressionability of DNA from the wave frequency, many experimental research studies have been carried out which have opened a new branch in science, called wave genome. Konstantin Meyl adapted the scalar waves described by Nicola Tesla to biology and proposed the relationship between the scalar waves and DNA^{1}. Greg Braddon and colleagues in 3 experiments investigated the impressionability of DNA from human emotions^{2}. Rein and Mccraty studied the impact of music on the DNA^{3,4,5}. Another study was carried out on the effect of sound waves on the synthesis and genes of chrysanthemum^{6}. Peter Garjajev and his research group proved that DNA can be reprogrammed by words and using the correct resonant frequencies of DNA^{7,8}. Russian quantum biologist Poponin tried to prove that human DNA has a direct effect on the physical world using some experiments^{9}. Also, he found out that our DNA can cause disturbing patterns in the vacuum, thus producing magnetized microscopic wormholes^{10}. Nobel Prize-winning scientist Luc Montagnier known for his study on HIV and AIDS, claims to have demonstrated that DNA can be generated by teleportation through quantum imprint and also showed that DNA emits electromagnetic signals that teleport the DNA to other places, such as water molecules^{11,12}.

About Mathematical Modeling of DNA, three eminent mathematical models have been proposed for describing DNA. The most famous model presented is the PBD (Peyrard-Bishop-Dauxois) model^{13}, although this model has been updated by various researchers^{14,15,16,17,18}, the models presented in these studies usually have at least three significant weaknesses such as being discrete, not being outside the plane, not being spiral, and not considering the position of nucleobases*.* Examples of the PBD model are displayed in Fig. 1. The other model is a rod model looking at DNA on a larger scale^{19,20,21}, with significant weaknesses, including the lack of attention to the nucleobase positions and the hydrogen bond, as well as considering DNA with one strand (Fig. 2). There is another model called SIDD (stress-induced DNA duplex destabilization) that is entirely mathematical and applied in the field of molecular dynamics^{22,23}. These models mentioned were mostly designed to investigate the vibration of DNA, and there are several models available in other fields such as DNA’s entropic elasticity^{24} and bending of a DNA^{25,26}. A beam is a structural element that primarily resists loads applied laterally to the beam axis. Its mode of deflection is primarily by bending. The Timoshenko beam theory takes the shear deformation and rotational bending effects into consideration for describing the behavior of thick beams. On the other hand, all previous studies related to curved beam vibrations focus on out of plane vibration of curved beams (with inline coordinates) and do not study out of plane vibration of the curved beams. A.Y.T. Leung is the only reference derived from the governing equations for a helical beam with rectangular cross-sections with pre-twist^{27}.

According to the above contents based on the necessity of DNA vibration analysis and the weaknesses of the previously proposed models, the necessity of carrying out this research becomes more evident. The dynamic model presented for DNA here has been named GMDM (Ghadiri Marvi DNA Model), it is provided by connecting two out of plane nano curved beams with spring and a damper. Each of the beams is a model for one of the two DNA strands (sugar-phosphate backbone). Also, spring and damper is a model for hydrogen bonds between nitrogen-containing nucleobases. The effects of the nucleobases (cytosine, guanine, adenine, and thymine) are also considered with their mass (Fig. 3). Also, the effects of DNA surrounding fluid (nucleoplasm) have been applied by using the Navier Stokes equations. Effects of temperature change on DNA are also applied to equations with external work. Finally, by using the relations of all the effects that mentioned above and using the Hamilton principle, the DNA equations will be extracted.

It should be noted that with the help of the theory of nonlocal, the effects of size were considered. By solving these equations, DNA natural frequency will be obtained for the first time. Numerical method can solve the equations derived from Hamilton’s principle method. The generalized differential quadrature method (GDQM) is one of the most numerical methods can be used for solving governing differential equations.

The idea brought up in this study is that if DNA is influenced by wave frequencies as much as DNA natural frequency then resonance occurs and DNA vibrates with large amplitude oscillations. With many shakes, DNA strands go up and down and at this moment a nucleobase in the nucleotide of one of the DNA strands establishes a hydrogen bond with the nucleobases lower or upper and the nucleobase aligns itself on the other DNA strands. This mechanism and idea which is presented schematically in Fig. 4, can cause disorganization in a sequence of DNA, and finally, with the help of a restriction enzyme (endonuclease), DNA in the cancerous cell loss of the ability to biosynthesis of proteins and the cancer is controlled just like the CRISPR/CAS9 technology.

## Nonlocal Elasticity Theory

In accordance with the nonlocal elasticity theory, the stress state at a reference point in an elastic continuum depends not only on the strain components at the same position but also a function of strains of all points in the neighbor regions. Therefore, the nonlocal stress tensor’s at point *x* is expressed as:

where *T(x)* is the classical, macroscopic stress tensor at point x related to strain by Hooke’s law with Eq. (1)*. C*(x) is the fourth-order elasticity tensor, ε(*x*) is the strain tensor and K(|*x*′ − *x*|, *τ*) denotes nonlocal modulus. *|x*′ − *x*| represents the distance and *τ* is a material constant that depends on the internal and external characteristic length defined as τ=e0aLτ=e0aL where *e*_{0} is a constant appropriate to each material, *a* is an internal characteristics length (e.g., bonds length) and *L* is an external characteristics length (e.g., wavelength).

Solving the integral constitutive relation in Eq. (1) is difficult. Thus, equivalent relation in a differential form was proposed by Eringen^{28,29} as follows:

∇^{2} is the Laplacian operator

### Beam theory and displacement of a double helical nanobeam

An out of plane curved beam is the beam having a twist and curvature around its central axis. The geometry situation of this model needs to choose the coordinate system that every moment changes its vector location (Fig. 5). Thus, for the analysis of the out of plane curved beam, Frenet triad must be used.

- Article
- Open Access
- Published:

# A Mathematical Model for Vibration Behavior Analysis of DNA and Using a Resonant Frequency of DNA for Genome Engineering

*Scientific Reports* **volume 10**, Article number: 3439 (2020)

## Abstract

The DNA molecule is the most evolved and most complex molecule created by nature. The primary role of DNA in medicine is long-term storage of genetic information. Genetic modifying is one of the most critical challenges that scientists face. On the other hand, it is said that under the influence of acoustic, electromagnetic, and scalar waves, the genetic code of DNA can be read or rewritten. In this article, the most accurate and comprehensive dynamic model will be presented for DNA. Each of the two strands is modeled with an out of plane curved beam and then by doubling this two strands with springs, consider the hydrogen bond strength between this two strands. Beams are traditionally descriptions of mechanical engineering structural elements or building. However, any structure such as automotive automobile frames, aircraft components, machine frames, and other mechanical or structural systems contain beam structures that are designed to carry lateral loads are analyzed similarly. Also, in this model, the mass of the nucleobases in the DNA structure, the effects of the fluid surrounding the DNA (nucleoplasm) and the effects of temperature changes are also considered. Finally, by deriving governing equations from Hamilton’s principle method and solving these equations with the generalized differential quadrature method (GDQM), the frequency and mode shape of the DNA is obtained for the first time. In the end, validation of the obtained results from solving the governing equations of mathematical model compared to the obtained results from the COMSOL software is confirmed. By the help of these results, a conceptual idea for controlling cancer with using the DNA resonance frequency is presented. This idea will be presented to stop the cancerous cell’s protein synthesis and modifying DNA sequence and genetic manipulation of the cell. On the other hand, by the presented DNA model and by obtaining DNA frequency, experimental studies of the effects of waves on DNA such as phantom effect or DNA teleportation can also be studied scientifically and precisely.

## Introduction

Deoxyribonucleic acid, more commonly known as DNA, is an evolved molecule that contains the genetic code of organisms. Every living thing has DNA within their cells. It is important for inheritance, coding for proteins and the genetic instruction guide for life and its processes. DNA holds the instructions for an organism or each cell development processes, reproduction and ultimately death. Over the past decades, empirical discussions have been proposed to modify the genes in the DNA. These changes have been much discussed in the medical field by drawing on applications such as treatment, preventing the development of cancer, or erupting an organ (for example, a tooth). On the other hand, it should be noted that only 3% of DNA capacity is considered in medical fields. In the last two decades, a topic called “wave genome” has been raised by Russian scientists, which states that 97% of other DNA is not only inapplicable but also has a more significant role; because DNA can be affected by acoustic, electromagnetic, and scalar waves. Under the influence of these waves, the genetic code can be read or rewritten. Another claim of Russian scientists is that DNA is a biological network that binds all humans. About impressionability of DNA from the wave frequency, many experimental research studies have been carried out which have opened a new branch in science, called wave genome. Konstantin Meyl adapted the scalar waves described by Nicola Tesla to biology and proposed the relationship between the scalar waves and DNA^{1}. Greg Braddon and colleagues in 3 experiments investigated the impressionability of DNA from human emotions^{2}. Rein and Mccraty studied the impact of music on the DNA^{3,4,5}. Another study was carried out on the effect of sound waves on the synthesis and genes of chrysanthemum^{6}. Peter Garjajev and his research group proved that DNA can be reprogrammed by words and using the correct resonant frequencies of DNA^{7,8}. Russian quantum biologist Poponin tried to prove that human DNA has a direct effect on the physical world using some experiments^{9}. Also, he found out that our DNA can cause disturbing patterns in the vacuum, thus producing magnetized microscopic wormholes^{10}. Nobel Prize-winning scientist Luc Montagnier known for his study on HIV and AIDS, claims to have demonstrated that DNA can be generated by teleportation through quantum imprint and also showed that DNA emits electromagnetic signals that teleport the DNA to other places, such as water molecules^{11,12}.

About Mathematical Modeling of DNA, three eminent mathematical models have been proposed for describing DNA. The most famous model presented is the PBD (Peyrard-Bishop-Dauxois) model^{13}, although this model has been updated by various researchers^{14,15,16,17,18}, the models presented in these studies usually have at least three significant weaknesses such as being discrete, not being outside the plane, not being spiral, and not considering the position of nucleobases*.* Examples of the PBD model are displayed in Fig. 1. The other model is a rod model looking at DNA on a larger scale^{19,20,21}, with significant weaknesses, including the lack of attention to the nucleobase positions and the hydrogen bond, as well as considering DNA with one strand (Fig. 2). There is another model called SIDD (stress-induced DNA duplex destabilization) that is entirely mathematical and applied in the field of molecular dynamics^{22,23}. These models mentioned were mostly designed to investigate the vibration of DNA, and there are several models available in other fields such as DNA’s entropic elasticity^{24} and bending of a DNA^{25,26}. A beam is a structural element that primarily resists loads applied laterally to the beam axis. Its mode of deflection is primarily by bending. The Timoshenko beam theory takes the shear deformation and rotational bending effects into consideration for describing the behavior of thick beams. On the other hand, all previous studies related to curved beam vibrations focus on out of plane vibration of curved beams (with inline coordinates) and do not study out of plane vibration of the curved beams. A.Y.T. Leung is the only reference derived from the governing equations for a helical beam with rectangular cross-sections with pre-twist^{27}.

According to the above contents based on the necessity of DNA vibration analysis and the weaknesses of the previously proposed models, the necessity of carrying out this research becomes more evident. The dynamic model presented for DNA here has been named GMDM (Ghadiri Marvi DNA Model), it is provided by connecting two out of plane nano curved beams with spring and a damper. Each of the beams is a model for one of the two DNA strands (sugar-phosphate backbone). Also, spring and damper is a model for hydrogen bonds between nitrogen-containing nucleobases. The effects of the nucleobases (cytosine, guanine, adenine, and thymine) are also considered with their mass (Fig. 3). Also, the effects of DNA surrounding fluid (nucleoplasm) have been applied by using the Navier Stokes equations. Effects of temperature change on DNA are also applied to equations with external work. Finally, by using the relations of all the effects that mentioned above and using the Hamilton principle, the DNA equations will be extracted.

It should be noted that with the help of the theory of nonlocal, the effects of size were considered. By solving these equations, DNA natural frequency will be obtained for the first time. Numerical method can solve the equations derived from Hamilton’s principle method. The generalized differential quadrature method (GDQM) is one of the most numerical methods can be used for solving governing differential equations.

The idea brought up in this study is that if DNA is influenced by wave frequencies as much as DNA natural frequency then resonance occurs and DNA vibrates with large amplitude oscillations. With many shakes, DNA strands go up and down and at this moment a nucleobase in the nucleotide of one of the DNA strands establishes a hydrogen bond with the nucleobases lower or upper and the nucleobase aligns itself on the other DNA strands. This mechanism and idea which is presented schematically in Fig. 4, can cause disorganization in a sequence of DNA, and finally, with the help of a restriction enzyme (endonuclease), DNA in the cancerous cell loss of the ability to biosynthesis of proteins and the cancer is controlled just like the CRISPR/CAS9 technology.

## Nonlocal Elasticity Theory

In accordance with the nonlocal elasticity theory, the stress state at a reference point in an elastic continuum depends not only on the strain components at the same position but also a function of strains of all points in the neighbor regions. Therefore, the nonlocal stress tensor’s at point *x* is expressed as:

where *T(x)* is the classical, macroscopic stress tensor at point x related to strain by Hooke’s law with Eq. (1)*. C*(x) is the fourth-order elasticity tensor, ε(*x*) is the strain tensor and K(|*x*′ − *x*|, *τ*) denotes nonlocal modulus. *|x*′ − *x*| represents the distance and *τ* is a material constant that depends on the internal and external characteristic length defined as τ=e0aLτ=e0aL where *e*_{0} is a constant appropriate to each material, *a* is an internal characteristics length (e.g., bonds length) and *L* is an external characteristics length (e.g., wavelength).

Solving the integral constitutive relation in Eq. (1) is difficult. Thus, equivalent relation in a differential form was proposed by Eringen^{28,29} as follows:

∇^{2} is the Laplacian operator

### Beam theory and displacement of a double helical nanobeam

An out of plane curved beam is the beam having a twist and curvature around its central axis. The geometry situation of this model needs to choose the coordinate system that every moment changes its vector location (Fig. 5). Thus, for the analysis of the out of plane curved beam, Frenet triad must be used.

The base beam theory used to model an out of plane curved beam is Timoshenko beam theory. By following Timoshenko’s assumption, the displacements *u* are defined as consisting of two parts, part one is the displacements at the centerline along the local axes *v* and part 2 is the rotations of the cross-section *θ*. This two parts according to Timoshenko’s assumption, only defined at the *x*_{3} axis.

In Eq. (3) *r* is the vector of the main displacement functions.

In Fig. 6, a small cut of a cross section of a beam has been displayed.

To find the strain vector of out of plane curved beams must be differentiated the displacement concerning the arc length *x*_{3}.

With using the Frenet triad and its special differential can be derived the displacement vector gradient of an out of plane curved beam as follows^{27}:

Parameters τ,κ,μτ,κ,μ used in the matrix Eq. (4) respectively represent tortuosity, curvature and pre-twist rate, in which, C=cos(μs)C=cos(μs), S=sin(μs)S=sin(μs) and μsμs is twist per length.

The non-vanishing strain matrix based on the Timoshenko’s assumption are derived as:

For simplicity in the use of nonlocal theory, strain matrix is separated into three matrices:

In which:

Now considering the nonlocal Eringen equation and replacing the strain term with the Eq. (6).

By defining axial force to form N=∫σdAN=∫σdA and bended moment as M=∫xσdAM=∫xσdA, we will have:

Where I1=∫x22dAI1=∫x22dA and I2=∫x12dAI2=∫x12dA and {M1,M2}=∫{x1,x2}σdA{M1,M2}=∫{x1,x2}σdA

For convenience, Eqs. (8–10) are brought into a matrix.

The GMDM model consists of two out of plane curved beams that are used to model a DNA. These two beams are connected with springs and dampers. By considering the mass of the nucleobases (Fig. 3b) and then by writing the strain energy and kinetic energy equations of the components, and the use of Navier-Stokes equations to apply the effects of a nucleoplasm, the effects of the temperature increase as an external work, and putting all of these equations in the Hamilton equation, the governing equation for the DNA model will be derived.

To derive the governing equations in the first step, we need to find the strain energy. Strain energy for two beams and also damper and spring connecting two beams will be as follows.

where:

Finally, the strain energy of all the components is brought together and the strain energy of the whole system is obtained as

For simplicity, Eqs. (12) and (13) can be demonstrated in the form of a matrix as the following:

It should be noted that in order to prepare the conditions for using variational method form, the matrix [e0e1e2]9∗12[e0e1e2]9∗12 was separated and rearranged to this matrix [e01e11e21e02e12e22]9∗(6+6)[e01e02e11e12e21e22]9∗(6+6), in which:

The kinetic energy of two beams and nucleobases with vibration frequency ωω is given by:

where ρρ is the density of strands of DNA and:

It should be noted that in Fig. 3 and GMDM model, the mass of nucleobases have been considered. Adenine always provides hydrogen bond with Thymine and also Guanine always makes hydrogen bond with Cytosine. Therefore, by the averages of weight of Adenine and Thymine, also weight of Guanine and Cytosine, attach this average mass to the both beams. Depending on the genetic sequence in the DNA code, it is determined that which one of the nucleobases kinetic energy (T1AT,T1CG,T2AT,T2CGT1AT,T1CG,T2AT,T2CG) should be used.

The effects of the nucleoplasm on the vibration of DNA are also considered as an external work, which the Navier Stokes equations are also used to apply this effect. The effects of the nucleoplasm for one beam and for two directions x1x1 and x2x2 that are perpendicular to the cross-section are shown below:

In which *ρ*_{n}, *P and μ*_{v} are density of nucleoplasm, intracellular pressure and viscosity of nucleoplasm respectively.

The effects of fluid are considered as an external work.

In which:

Regardless of gravity and nonlinear term, the effects of fluid on DNA strands will be as follows:

where W1S1W1S1 represents the external work applied to number 1 strand in x2x2 directions and W2S1W2S1 represents the external work applied to number 2 strand in x1x1 directions. Also, the effects of temperature increase on DNA as linear in term of thickness of the beam will also be applied as an external work as below (since strands are larger than nucleobases, we neglected the effect of temperature on nucleobases).

In which:

By replacing the Eq. (24) in to the Eq. (23), we will have:

The Eq. (25) can be written as follows:

In which:

By substituting the strain energy (Eq. (16)) kinetic energy (Eq. (19)) external works (Eqs. (21, 26)) and also using variation method and put them in Hamilton relation, we will have:

By merging Eq. (11) in the Eq. (27), the governing equations of DNA with considering the effects of the fluid and temperature change will be as: