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

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-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-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-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].

Figure 1

Schematic of PBD model (a) and its updates (b–d). These images have been drawn with Adobe Photoshop CC 2018.

Figure 2

A rod model that looks at DNA on a larger scale. This image has been designed with Adobe Photoshop CC 2018.

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.

Figure 3

(a) Imagined shape for DNA[43] © 2013 Nature Education Adapted from Pray, L. (2008) Discovery of DNA structure and function: Watson and Crick. Nature Education 1(1):100. All rights reserved. (b) The GMDM Mathematical model, which is presented in this paper for dynamics investigations of DNA. This image has been modeled with COMSOL Multiphysics 5.3a and edited with Adobe Photoshop CC 2018.

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.

Figure 4

(a) DNA’s shape before reaching resonant frequency. (b) DNA’s shape after reaching resonant frequency. (c) Establishes a hydrogen bond between the nucleobases upper with nucleobases lower. (d) Removal of additional and remaining nucleotides with the help of restriction enzyme. (e) New DNA’s shape and disorganization in a sequence of DNA. These images have been modeled with COMSOL Multiphysics 5.3a and edited with Adobe Photoshop CC 2018.

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 where e0 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.

Figure 5

An out of plane curved beam. This image has been modeled with CATIA V5 and edited with Adobe Photoshop CC 2018.

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 x3 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.

Figure 6

A small cut of a cross section of a beam. This image has been modeled with Adobe Photoshop CC 2018.

To find the strain vector of out of plane curved beams must be differentiated the displacement concerning the arc length x3.

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, , and 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 and bended moment as , we will have:

Where and and

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 was separated and rearranged to this matrix , 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 () 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 and that are perpendicular to the cross-section are shown below:

In which ρ, P and μ 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 represents the external work applied to number 1 strand in directions and represents the external work applied to number 2 strand in 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:

This model is the most comprehensive model to investigate the dynamic behavior of DNA considering mass of nucleobases, their hydrogen bond and effect of the surrounding fluid of DNA, is extracted to 12 PDE equations.

Also boundary condition equation of this system is as below:

The special boundary conditions that was used in are defined as clamped-clamped (at x = 0 & x = L), and simply supported (at x = 0 & x = L) as follows:

Clamped:

Simply supported:

Solution with using GDQM

Different numerical techniques can be used to solve the governing equations and related boundary conditions. In this paper the GDQ method that is introduced by Shu[30] is used. This technique has been successfully employed to solve a variety of problems in vibration analysis and dynamical systems. The GDQM is a powerful method that can be used to solve partial differential equations extended and generalized with high accuracy and convergence and performance.

According to GDQM the order derivative of function of with respect to x at is:Where N is the number of total grid points and C(n) is the respective weighting coefficients matrix that can be obtained as followswhere M(x) is defined as below:

The weighting coefficients for order are derived from:

The Chebyshev-Gauss-Lobatto technique can be used for the distribution of grid points in domain:

By writing the governing equations as a matrix form:whereas is the mass matrix and is stiffness matrix

With coupling this six matrix equations and using computer software, natural frequency is found.

Geometry and mechanical properties

Dimensions and mechanical properties of DNA are depicted in Tables 1 and 2 according to various scientific references.

Table 1

mechanical properties of DNA.

Mechanical or thermal properties	Symbol	Quantity	Reference
Young modulus	E	0.3 (GPa)	[44]
Poisson’s ratio	ν	0.5	[45]
Shear modulus	G	0.1 (GPa)	G=E2(1+ν)
Density	ρ	1.7 (g/m3)	[46]
Mass of Adenine	M	226 * 10−27 (Kg)	[47]
Mass of Thymine	M	211 * 10−27 (Kg)
Mass of Guanine	M	252 * 10−27 (Kg)
Mass of Cytosine	M	185 * 10−27 (Kg)
Hydrogen bond strength	Between Adenine-Thymine (kAT)	19.5 (N/m)	[18]
	Between Guanine-Cytosine (KGC)	56.3 (N/m)
Damper constant	C	0.05 (N s/m)	[48]
Density of nucleoplasm	ρn	0.14 (g/cm3)	[49]
Viscosity of nucleoplasm	μn	1.35 (cP) or 0.135 * 10−2 (N.sm^2)	[50]
Osmotic pressure	P	4 (atm)	[50]
Thermal conductivity	kt	k=150Wmk	[51]
Specific heat capacity	cp	C_P=40calmolk or2.56∗10^−4Kjgrk	[52]
Thermal diffusivity	α	α=3.44∗10^−8m^2s	α=kρc_P

Table 2

Geometry properties of DNA.

Geometry	Quantity	Reference
Pitch/turn of helix	34 (Å)	[53]
Rise/bp along axis	3.4 (Å)
Radius	10 (Å)
Rotation/bp	34.3°	[54]
Curvature	0.008 * 1010	κ=radiusradius^2+pitch^2
Twist	0.027 * 1010	τ=pitchradius^2+pitch^2
Straight and open length	7.14 * 10−9	L=(2πr)^2+(pitch)^2
Length of Adenine	2 * 5.8 (Å)	[47]
Length of Thymine	2 * 4.8 (Å)
Length of Guanine	2 * 5.7 (Å)
Length of Cytosine	2 * 4.7 (Å)	[55]
Distance between Adenine-Thymine	2.83 (Å)

Adenine-Thymine (k)

Guanine-Cytosine (K)

0.135 * 10−2

One of the challenges ahead of the present study was to obtain a DNA cross section. This work was carried out in a way that with an average weighted gain from the diameter of the atoms participating in the phosphate and sugar structure[31,32], the thickness of the DNA strands was measured. In addition, by using the length of the atomic bonds in the sugar and phosphate structure[33], and with the help of the angles between bonds[34], the width of the DNA backbone was obtained. The results of these calculations are presented in the Table 3.

Table 3

Cross-section of DNA strands.

Thickness	Width
0.98 (Å)	3.17 (Å)

The sectional average radius of nucleobases is also obtained by an average weighted gain from the diameter of the atoms participating in the nucleobases structure. The sectional radius of nucleobases are shown in Table 4.

Table 4

The sectional radius of nucleobases (Å).

Adenine	Thymine	Guanine	Cytosine
0.53	0.5	0.55	0.45

The second challenge ahead of this study was to obtain a suitable nonlocal parameter for DNA. Parameter is directly related to the aspect ratio of geometry. Depending on the ratio of the length of the DNA to its thickness, and by extrapolating the results presented by Ghavanloo[35], is obtained. Also Parameter a depends on the length of atomic bonds in the target objects, therefore, with an average of length of all atomic bonds in a nucleotide, a will also be obtained. These two parameters have been brought in the Table 5. It should be noted that the calculated value for e0a, applies in the general relation presented for the nonlocal parameter[36,37]. This relationship is in this form: .

Table 5

nonlocal parameter.

e0	a	e0a	e_0aL
17.87	1.348	24.10	0.337

Nuclear DNA does not appear in free linear strands. DNA is wrapped around histones in order to fit within the nucleus and participate in the chromatin structure. According to the cause, the supports in this study are assumed clamped-clamped.

Convergence

The number of grid points affects the natural frequency of DNA. In Table 6 evaluating convergence have been proved with respect to number of elements for the first three frequencies of DNA.

Table 6

Convergence of results.

Mode number	Natural frequency (*109 Hz)
	N = 12	N = 13	N = 14	N = 15	N = 16	N = 17	N = 18	N = 19	N = 20
1	5.0199	5.0199	5.0198	5.0198	5.0198	5.0198	5.0198	5.0198	5.0198
2	7.0388	7.0445	7.0447	7.0443	7.0443	7.0443	7.0444	7.0444	7.0444
3	9.6710	9.6738	9.6967	9.6960	9.6938	9.6938	9.6939	9.6939	9.6939

It is observed that for 20 node later, the answers will converge.

Study of validation

Using simulations in COMSOL software, we obtained the natural frequency of DNA. Then the natural frequency obtained using mathematical modeling is compared with COMSOL simulations in Table 7.

Table 7

Validation of results by COMSOL.

Mode number	Natural frequency from present mathematical model (Hz)	Natural frequency from COMSOL (Hz)	Difference percentage
1	5.0198 * 109	4.6211 * 109	7%
2	7.0444 * 109	7.1816 * 109	2%
3	9.6939 * 109	11.168 * 109	13%

In short articles, DNA natural frequency is mentioned approximately with simulations (Table 8).

Table 8

Approximate DNA frequency obtained in the previous studies.

	[56]	[57]	Present (3 first modes)
Frequency	Above 1 GHz	0.2–10 GHz	5–9 GHz

It should be noted that the primary modes of frequency of DNA obtained in the present study, lie within the range considered in those studies.

Results due to changing parameters

The derived equations, are solved with the mathematical model to find the natural frequencies and mode shapes for DNA and a parametric study is performed on the frequency of the DNA.

The effect of the DNA length on frequency

As mentioned, the length of the DNA varies between the two histones and has 10 and 80 nucleobases. These results which are presented in Table 9 show that with increasing number of nucleobases, the natural frequency decreases. The reason for this change is that by increasing the length of the body, it has more mobility and less rigidity.

Table 9

The effect of the DNA length on frequency.

Number of Nucleobases	Number of turn of helix	Natural frequency (Hz)
		First mode	Second mode	Third mode
10	1	5.0198 * 109	7.0444 * 109	9.6939 * 109
20	2	1.4375 * 109	2.6497 * 109	3.8139 * 109
30	3	8.2032 * 108	1.1533 *109	1.8835 * 109
40	4	6.6750 * 108	6.8998 * 108	1.3232 * 109
50	5	5.3081 * 108	5.4125 * 108	1.0504 * 109
60	6	4.3337 * 108	4.5284 * 108	8.6065 * 108
70	7	3.7059 * 108	3.8425 * 108	7.3720 * 108
80	8	3.2702 * 108	3.3004 * 108	6.5024 * 108

The effect of supporting conditions on frequency

It can be seen that clamped-clamped supporting, has the highest frequency (Table 10) because by increasing constraint in supporting (reducing degrees of freedom), rigidity and frequency increase.

Table 10

The effect of supporting conditions on frequency.

Supporting Conditions	Natural frequency (Hz)
	First mode	Second mode	Third mode
Clamped-Clamped	5.0198 * 109	7.0444 * 109	9.6939 * 109
Simply supported-Simply supported	7.5297 * 108	2.8201 * 109	3.0451 * 109

The effect of DNA embedding fluid on frequency

As we know, chromatins and DNA are floating in a fluid called the nucleoplasm. In Table 11, the effects of the presence and absence of this fluid on DNA frequency, have been shown. It is observed that due to the low density of the nucleoplasm, the DNA frequency decreases slightly in the presence of this fluid. The reason for this topic can be found in the absorption of vibrating energy of DNA by the fluid.

Table 11

The effect DNA embedding fluid on Frequency.

Impact of Fluid	Frequency (Hz)
	First mode	Second mode	Third mode
Without considering the effects of the fluid	5.0209 * 109	7.0465 * 109	9.6982 * 109
In the presence of Nucleoplasm	5.0198 * 109	7.0444 * 109	9.6939 * 109

The effects of temperature rise on frequency

Today, the huge effects of temperature increase on the stability of DNA strands, is proven. As DNA close to the melting temperature, displacement of the base pairs increases sharply and hydrogen bound becomes unbound[13,38]. DNA persistence length and internal base pair conformations, strongly depends on temperature[39,40]. Increasing temperature enhances flexibility of the double-stranded DNA chain[41]. As it is said in 354 K, the DNA denatures[42]. In Table 12 the effects of temperature rise on the frequency of DNA, was investigated. It is observed that the DNA frequency decreases with increasing temperature. The reason for this event is that the DNA rigidity is reduced by temperature rise.

Table 12

The effects of temperature rises on frequency.

Temperature increase		Frequency (Hz)
		First mode	Second mode	Third mode
No temperature increase		5.0198 * 109	7.0444 * 109	9.6939 * 109
10° increase in temperature	Mathematical modeling	4.7304 * 109	7.0319 * 109	9.3957 * 109
	COMSOL	4.1906 * 109	6.7811 * 109	10.914 * 109
20° increase in temperature		4.3656 * 109	7.0195 * 109	9.1238 * 109
44° increase in temperature		2.5630 * 109	6.9894 * 109	7.9386 * 109

It is also observed that at 44 °C increase in temperature (reaching 354 K), the DNA natural frequency, drops sharply.

Mode shape of DNA in COMSOL

As mentioned earlier the results were verified with the help of the DNA modeling in COMSOL. In this software, both strands of DNA are modeled with two spiral geometries and the nucleotides are also designed to be cylinders with different radius and height according to their dimensions. It should be noted that the distance between the nucleobases is taken as the length of the hydrogen bond. The nucleotides are arbitrarily arranged in sequence ATCGTAGCATCGTAGCATCG. Hydrogen bonds also have been modeled as springs. Since the DNA strands are inside the nucleobases and in order to apply the properties and forces of fluid in the COMSOL, we used a sphere larger than DNA in which DNA is enclosed by fluid (Fig. 7).

Figure 7

Model of DNA that has been designed with COMSOL Multiphysics 5.3a.

In Fig. 8, the mode shape of the DNA vibrations, obtained in COMSOL software are presented.

Figure 8

First 20 mode shape of DNA’s vibration (3D). These images have been modeled and simulated with COMSOL Multiphysics 5.3a and edited with Adobe Photoshop CC 2018.

Conclusion

The presented study on DNA is the most accurate dynamic model that was named GMDM. This model 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, also spring and damper is a model for hydrogen bonds between nitrogen-containing nucleobases. DNA vibrations and DNA natural frequency is investigated for the first time with considering thermal, mass of nucleobases and fluid effect on GMDM. These effects were applied using the equations like Navier Stokes in to the Hamilton principle. With solving these equations by GDQM, DNA natural frequency will be obtained for the first time. The novelty of this model are as follows:

Being out of plane, spiral with twist and curvature

Being continuous model

Considering the effects of the mass and viscoelastic of the nucleoplasm

Considering the effects of temperature rise

Considering the position of the nucleobases.

As the energy absorption by the nucleoplasm is very insignificant, it can be concluded that during the resonance, the DNA amplitude oscillations range is very large and with severe shakes and using a restriction enzyme, sequence of DNA might be disorganized, and DNA in cancerous cells loses its ability for proteinization and consequently in this way cancer may be controlled.