Influences of Electromagnetic Energy on Bio-Energy Transport through Protein Molecules in Living Systems and Its Experimental Evidence

The influences of electromagnetic fields (EMFs) on bio-energy transport and its mechanism of changes are investigated through analytic and numerical simulation and experimentation. Bio-energy transport along protein molecules is performed by soliton movement caused by the dipole–dipole electric interactions between neighboring amino acid residues. As such, EMFs can affect the structure of protein molecules and change the properties of the bio-energy transported in living systems. This mechanism of biological effect from EMFs involves the amino acid residues in protein molecules. To study and reveal this mechanism, we simulated numerically the features of the movement of solitons along protein molecules with both a single chain and with three channels by using the Runge–Kutta method and Pang’s soliton model under the action of EMFs with the strengths of 25,500, 51,000, 76,500, and 102,000 V/m in the single-chain protein, as well as 17,000, 25,500, and 34,000 V/m in the three-chain protein, respectively. Results indicate that electric fields (EFs) depress the binding energy of the soliton, decrease its amplitude, and change its wave form. Also, the soliton disperses at 102,000 V/m in a single-chain protein and at 25,500 and 34,000 V/m in three-chain proteins. These findings signify that the influence of EMFs on the bio-energy transport cannot be neglected; however, these variations depend on both the strength and the direction of the EF in the EMF. This direction influences the biological effects of EMF, which decrease with increases in the angle between the direction of the EF and that of the dipole moment of amino acid residues; however, randomness at the macroscopic level remains. Lastly, we experimentally confirm the existence of a soliton and the validity of our conclusion by using the infrared spectra of absorption of the collagens, which is activated by another type of EF. Thus, we can affirm that both the described mechanism and the corresponding theory are correct and that EMFs or EFs can influence the features of energy transport in living systems and thus have certain biological effects.


Introduction
Prevalent in our environment, electromagnetic fields (EMFs) or waves (EMWs) of varying frequencies and strengths are generated by such sources as the irradiations of high-voltage transmission lines, electrical appliances, microwave stations, and radio equipment. As such, it is necessary to understand whether these externally applied EMFs, EMWs, or electric fields (EFs) have any biological effects on animals and humans [1,2]. Much research has focused on determining any relationship between EMFs and biological processes . Li et al. [10][11][12][13][14][15] studied the biologic effects of environmental EMFs, observing the variations of the dielectric constant and conductivity and absorbed strength of infrared light of bio-tissues as well as of the rheology features of the blood in animals under the influences of different electromagnetic fields relative to those without electromagnetic fields. Zhang et al. [16,17] measured the expression of matrix metalloproteinase and tight-junction proteins in rats by following a pulse EMF at high frequencies, which induced blood-brain barrier permeability changes, and also studied the effects of a pulse EMF with high frequencies on the fluorescence spectrum of serum in the rats. Beebe et al. [18] studied the nanosecond pulsed electric field effect on cell growth and the development of bio-tissues, as well as apoptosis induction and tumor growth inhibition. Lazzaro et al. [19] and Apollonio et al. [20] discussed the feasibility of microwave energy affecting biological systems via nonthermal mechanisms and offered systematic approaches, respectively. Hofmann et al. [21] investigated the possibility of electroporation therapy for head and neck cancer. Kekez et al. [22] researched contributions to the biophysics of lethal effects of EFs on microorganisms. Pang et al. [23] investigated and discussed the mechanism and properties of biological thermal effects of microwaves and noted that microwaves can influence the proliferations of cells and tissues. These results indicated that EMFs and EFs have certain biological effects. Although animals or humans are not measurably affected by EMFs unless they come in close contact with high-voltage transmission lines, research should focus more on the proliferation of microwaves and radio waves and their influences on health. Hence, the biological effects of these influences on the health of human beings merit examination.
Only a few conclusions supporting the theory that EMFs and EFs containing microwave radio waves have biological effects have been accepted. Because of varying types of research and instrumentation, experimental results for the same questions range considerably. Thus, a widely accepted and unified conclusion is difficult to acquire, and the mechanism of influence of EMFs on life activity has not been revealed or elucidated [1][2][3]. However, epidemiological investigations have shown that EMFs, including microwaves, radio waves, and those from high-voltage transmission lines, have always had adverse effects on health [3][4][5][6][7][8][9][10][11][12][13][14][15][16]. Therefore, the influence of EMFs on the health of human beings and animals cannot be dismissed. The rapidly increasing usage of electrical appliances, microwave instruments, cell phones, and other electromagnetic devices has led to expansive distribution of EMFs in our environment. In this case, it is necessary to study and release the biological effects of EMFs and their influences on the health of human beings and animals. This need suggests that investigations should focus on the biological effects of EMFs-in particular, the mechanisms in depth through the latest ideas and methods in theoretical analyses and experimental measurements.
One key issue would be to analyze and investigate the electromagnetic properties of biotissues and biomacromolecules and the distributions and properties of charged groups, or atoms and molecules, as well as any variations under the influence of externally applied EMFs. Examination of these charged particles, molecules, or species, as well as their properties of movement, could lead to a better understanding of biological mechanisms under the effect of EMFs. In practice, these charged particles, molecules, or species exist widely in cells and biomacromolecules, such as protein and DNA [9,24], upon which the externally applied EMF can influence their biological processes and properties, and any biological effects can also be exhibited clearly. In this paper, we first seek any mechanism of biological effect from EMFs and EFs. Energy transport, released from the hydrolysis reaction of adenosine triphosphate (ATP) molecules along the protein molecules, is a typical target that exemplifies this described mechanism; i.e., EMFs can disturb and influence energy transport considerably through its interaction with the amino acid residues with certain electric dipole moments.
Here we reveal the mechanism of influence of EMFs on energy transport in protein molecules and further study the properties of this mechanism.
H D " ř n rε 0 Bǹ B n´J pBǹ B n`Bn Bǹ qs`ř n r P 2 n 2M`1 2 wpu n´un´1 q 2 s ř n rχ 1 pu n`1´un´1 qBǹ B n s " H ex`Hph`Hint (1) where Bǹ pB n q is the creation (annihilation) operator for an Amide I quantum (exciton) in the site n, u n is the displacement operator of amino acid residue at site n, P n is its conjugate momentum 1 χ = 6.2 × 10 −11 N is a nonlinear coupling parameter and represents the size of the exciton-phonon interaction in this process, 0 ε = 0.205 eV is the energy of the Amide I quantum (exciton), J is the dipole-dipole interaction energy between neighboring amino acid residues, J = 1.55 × 10 −22 or J = 9.68 × 10 −4 eV, and 0 r = 4.5 × 10 −10 m is the average distance between the neighboring amino acid residues [24][25][26][27][28][29][30][31][32][33]. The wave function of the systems proposed by Davydov [24][25][26][27][28][29] is the form of  (1) and (2) in the semiclassical limit and using the continuum approximation [24][25][26][27][28][29], has the form of If,  The wave function of the systems proposed by Davydov [24][25][26][27][28][29] is the form of D 2 ptq " |ϕ D y |βptqy " ÿ n ϕ n ptqBǹ expt´i ÿ n rβ n ptqP n´πn ptqu n su |0y (2) where |0y " |0y ex |0y ph , |0y ex and |0y ph are the ground states of the exciton and phonon, respectively, Davydov's soliton, obtained from Equations (1) and (2) in the semiclassical limit and using the continuum approximation [24][25][26][27][28][29], has the form of If, |ϕ D ptqy " ř n ϕ n ptqBǹ |0y ex in Equation (2) is an eigenstate of the number operator, N " ř n Bǹ B n , then the soliton contains only one exciton because N " xϕ D ptq|N |ϕ D ptqy " 1, i.e., N |ϕ D ptqy " 1 |ϕ D ptqy. This finding indicates that the Davydov soliton contains only one excitation, which corresponds to the excited state of a single particle and is localized over a scale of r 0 {µ D ; where µ D " χ 2 2 {p1´s 2 qwJ, s 2 " ν 2 {ν 2 0 , and ν 0 " r 0 pw{Mq 1{2 ; where ν 0 is the sound speed in the protein molecular chains; ν is the velocity of the soliton; and G D " 4Jµ D is its nonlinear interaction energy accepted in this process. This finding shows that Davydov's soliton is formed through the self-trapping of one exciton and that it has a binding energy E BD "´χ 4 1 {3Jw 2 . These results clearly exhibit that the energy released from the hydrolysis reaction of the ATP molecule promoted the form of a soliton and its transported along the protein molecules in a bell-type solitary wave or a soliton with an invariable amplitude and velocity, as given in Equation (3) [9,27,34]. This finding indicates that the energy cannot be damped and dissipated in the transport process and is significant for biological processes because energy is retained in the transport process due to the feature and essence of the soliton. Thus, these life activities clearly require the soliton.
Davydov's idea yields a compelling picture of the mechanism of energy transport in living systems and has led to extensive research in biophysics . However, the issues related to Davydov's model, including its foundation, accuracy quantum and classical properties, thermal stability, and lifetimes, have been the focus of much research . These investigations and discussions focus mainly on the validity of the Davydov theory and the thermal stability of the Davydov soliton at the biological temperature 300 K. Some numerical simulations have indicated that the Davydov soliton is not stable at this temperature [36][37][38][39][40][41][42][43][44][45][46][47]. At the same time, Monte Carlo numerical calculations indicate that the correlation characteristic of soliton-like quasiparticles occurs only at low temperatures [42][43][44], approximately T < 10 K, for widely accepted parameter values. This finding is consistent at a qualitative level with the results of Cottingham et al. [46] and Schweitzer [47]. The latter is a straightforward quantum-mechanical perturbation calculation, in which the lifetime of the Davydov's soliton is approximately 10´1 2 -10´1 3 s at 300 K, in which the soliton can transport only over approximately 10 amino acid residues. Therefore, Davydov's theory is not suitable for protein molecules. In addition, Förner's investigations showed that Davydov's soliton is stable only at 40 K and that it disperses completely at higher temperatures [39][40][41][42][43]. These results indicate clearly that Davydov's soliton is not a real carrier of the energy transport in protein molecules; thus, Davydov's theory is not appropriate to the systems. This finding demonstrates the necessity of developing new theories of energy transport in living systems.

Pang's Theory of Energy Transport and Its Properties
On the basis of the difficulties described concerning Davydov's theory and the results researched by Cruzeiro-Hansson [37,38] and Förner et al. , we improve Davydov's model by changing simultaneously the Hamiltonian and the wave function of the systems, in which we added a new coupling interaction of the excitons with the displacement of amino acid residues into the Hamiltonian in Equation (1) and replaced further the Davydov's wave function of the one-quantum (exciton)excited state in Equation (2) by a quasi-coherent two-quantum state . In this case, the representations in Equations (1) and (2) for the single-channel protein molecules are replaced by |Φptqy " |αptqy |βptqy " 1 λ r1`ř n α n ptqBǹ`1 2! p ř n α n ptqB 2 b q 2 s |0y exê xpt´i ř n rβ n ptqP n´πn u n |0y ph (4) and H " H ex`Hph`Hint " ř n rε 0 Bǹ B n´J pBǹ B n`Bn Bǹ qs`ř n r P 2 n 2M`1 2 wpu n´un´1 q 2 sř n rχ 1 pu n`1´un´1 qBǹ B n`χ2 pu n`1´un qpBǹ`1B n`Bǹ B n`1 qs respectively; where the created and annihilated operators of the exciton are represented by Bǹ and B n , respectively; |0y ex and |0y ph are the ground states of the exciton and phonon, respectively; and u n and P n are the displacement and momentum operators of amino acid residue at the site n, respectively. The β n ptq " xΦptq| u n |Φptqy, and π n ptq " xΦptq| P n |Φptqy are two sets of unknown functions; and λ is a normalization constant. Present nonlinear coupling constants are χ 1 and χ 2 = (10´15)ˆ10´1 1 N, which represent the modulations of the on-site energy and dipole-dipole interaction energy of excitons due to the variations of displacements of amino acid residue in the protein molecules, respectively. Other parameters are same as those in the Davydov's model mentioned previously. Using Equations (4) and (5) and working from the Schrödinger equation and the Heisenberg equation, we obtained i . α n ptq " ε 0 α n ptq´Jrα n`1 ptq`α n´1 ptqs`χ 1 rq n`1 ptq´q n´1 ptqsα n ptqχ 2 rq n`1 ptq´q n ptqsrα n`1 ptq`α n´1 ptqs5 q n ptq " Wrq n`1 ptq´2q n ptq`q n´1 ptqs`2χ 1 r|α n`1 | 2´| α n´1 ptq|s2 χ 2 αn ptqrα n`1 ptq´α n´1 ptqs`α n ptqrαn`1ptq´αn´1ptqs ( In the continuum approximation, we get from Equations (6) and (7) and where Rptq "" ε 0´2 J`5 2 tWptq´1 2 ř m r . β n ptqπ m ptq´. π n ptqβptqsu.
The soliton solution of Equation (8) is denoted by With µ P " 2pχ 1`χ2 q 2 wp1´s 2 qJ and the nonlinear interaction en G P " 8pχ 1`χ2 q 2 wp1´s 2 q and s " ν{ν 0 . In this case, the energy of the soliton in Equation (10), or the energy transported by the soliton, is obtained by The static (rest) energy of the soliton is where W " r2pχ 1`χ2 q 4 s{3w 2 J is the energy of deformation of the amino acid residues. The effective mass of the soliton is M sol " 2m ex`8 pχ 1`χ2 q 4 p9s 2`2´3 s 4 q 3w 2 Jp1´s 2 q In such a case, the binding or forming energy of the soliton in Pang's theory is denoted by These mathematical models clearly indicate that the properties of Pang's soliton  differ from those of Davydov's. Through concrete calculations, we found that Pang's soliton is thermally stable and has a sufficient lifetime of 10´9-10´1 0 s at the physiological temperature of 300 K, which is approximately 800ˆthat of Davydov's soliton (10´1 2´1 0´1 3 s), in which Pang's soliton can transport more than several hundreds of amino acid residues. Thus, Pang's soliton plays an important role in the biological processes. The comparison of Pang's soliton with Davydov's soliton is given in Table 1, which clearly indicates that Pang's soliton (or theory) differs considerably from that of Davydov's soliton (or theory). Since then, the results of Pang's theory were also confirmed by other studies [75][76][77][78][79][80][81][82][83][84][85][86]. Therefore, Pang's soliton is a carrier of energy transport, and the Pang's model can be applied to energy transport in protein molecules [97,[101][102][103][104][105][106]. Thus, we here applied the Pang's model to study the influences of EFs and EMFs on properties of energy transport in α-helical protein molecules ( Figure 1).

Mechanism of Influence of Electric Field on Energy Transport in Protein Molecules
Equations (11)- (14) indicate that the properties of the energy transported by the soliton involving its velocity (v), effective mass (M sol ), energy (E), and rest energy (E 0 ), which are determined by the features of the exciton, the amino acid residues, and the structure of the protein molecule, including the effective mass of the exciton (m ex ), the mass of the amino acid residue (M), the elastic constant of the protein molecular chains (w), the nonlinear exciton-phonon interaction coefficients (χ 1 and χ 2 ), the Amide I vibrational energy (ε 0 ), the dipole-dipole interaction energy between the neighboring amino acid residues (J), and the distance between the neighboring amino acid residues (r 0 ). This clearly indicates that the properties of the energy transport or the soliton will change once these characteristic parameters of the exciton, amino acid residues, and protein molecule structures are altered under the influences of externally applied electromagnetic or light fields, as well as temperatures. On the basis of these characters and features, we seek the mechanism that influences the EMFs or the EFs on energy transport in protein molecules.
As described, energy released from the hydrolysis reaction of the ATP molecules is transported along protein molecules in the solitons by virtue of the dipole-dipole interaction between neighboring amino acid residues with certain electric dipole moments. The externally applied EFs or EMFs can interact with these amino acid residues according to the EMF theory; thus, the properties of the energy transport or movement of soliton will vary, along with the externally applied EFs or EMFs, because the latter can vary the sizes and directions of the dipole-dipole interaction between the neighboring amino acid residues in such a case. Thus, some biological effects arising from the variations of the biological energy in protein molecule will occur correspondingly, which is the influencing mechanism of EMF or EF on the energy transport. In other words, this is a mechanism of biological effect through EMFs or EFs.
If we assume that the strength of the EF in the EMF is denoted by Ñ E and that the electric dipole moments of the amino acid residue in the protein molecules is denoted by interaction of the EF or the EMF with the amino acid residues can be denoted by Ñ E . Ñ p , according to the EMF theory [82,87,99]. In such a case, the dipole-dipole interaction energy between neighboring amino acid residues in the protein molecules will be changed from J to J`Ñ E . Ñ p . Then we can affirm that the properties of the energy transported by the soliton will be altered notably because the soliton in Equations (10)- (14) is very sensitive to variations of the dipole-dipole interaction between neighboring amino acid residues. This conclusion will be confirmed by the following results.

Analytic Results for Changes of Properties of the Soliton Transporting the Energy under the Influence of EFs
As mentioned previously, when protein molecules are exposed in an EMF or an EF, the dipole-dipole interaction between the neighboring amino acid residues in the protein molecules is changed into J`Ñ E . Ñ p . Then the variation of features of the solitons in Equations (10)- (14) in Pang's model will appear in this case [75][76][77][78][79][80][81][82][83][84][85][86]. These variations in the soliton from Equations (10)- (14) can be obtained, which are described as follows: (1) The amplitude (µ P ) effective mass (M sol ), energy (E), rest energy (E 0 ), and binding energy E BP of the soliton are decreased and depressed after EF is applied because the physical parameters in Equations (11)- (14) are inversely proportional to the dipole-dipole interactional energy J. This means that the capability and value of the energy transported by the soliton is decreased because of the increases of dipole-dipole interaction between neighboring amino acid residues. If the strength of an EF in an EMF is very high, then the capability and value of the energy transport will depress considerably. In such a case, we could infer and suppose that some new biological effects will occur and that these new biological effects have just arisen from the EF (2) The EF varies the form and outline of the soliton in Equation (10) Therefore, the EF has a stronger biological effect on the former and no biological effect on the latter. This result indicates clearly that the biological effect depends on the direction of externally applied EF with respect to that of electric dipole moments of the amino acid residues. If the directions of EF are different, although their strengths are the same, then their biological effects are also different. This finding implies that the biological effects of EFs depend on both the strength and the direction with respect of the dipole moment of amino acid residues [81,86,97].
However, it is worth noting the direction from which the biological effects arise from the EF at the macroscopic level. The orientations of many proteins and amino acid residues within the proteins over any macroscopic region tend to be random. However, this condition depends on the biotissue structure because some biotissues are well organized. Clearly, the orientation of the EF to the structure will be important in this case. In other cases, the overall effect being calculated will not depend on the direction of the field.

Results in Single-Protein Chains
We now investigate the influences of EFs in EMFs on the energy transported by the soliton in protein molecules with single channels and three channels by using the numerical simulation method [107,108], respectively. In this case, the variation of the dipole-dipole interaction is represented by pJ`Ñ E . Ñ p q´J "ˇˇˇˇÑ EˇˇˇˇˇˇÑ pˇˇcosθ, which is related to the strengths of the EF in the EMF. Therefore, our purpose is to determinate the changed features of the movement of the soliton transporting the energy with varying EFs, which can be obtained by the fourth-order Runge-Kutta method [107,108]. In such a case, we must establish the dynamic equations of the soliton in this numerical simulation. Obviously, the dynamic equations of the soliton can be obtained from Equations (6) and (7) in Pang's model. Their derivations are depicted as follows.
. y " Wpq n`1´2 q n`qn´1 q`2χ 1 pαr 2 n`1`α i 2 n`1´α r 2 n´1´α i 2 n´1 q4 χ 2 rαr n pαr n`1´α r n´1 q`αi n pαi n`1´α i n´1 qs (18) |α n | 2 " |αr n | 2`| αi n | 2 (19) where ar n and ai n are real and imaginary parts of a n , respectively. Equations (15)- (19) are discretely coupled and can be used to determine the dynamic features and states of the soliton by using the fourth-order Runge-Kutta method in the numerical simulation. Obviously, in this case there are four equations for one amino acid residue from Equations (15)- (18). Thus, we should find solutions of 4 N associated equations for the protein molecules constructed by N amino acid residues. However, when the fourth-order Runge-Kutta method [107,108] is used to find the solutions of Equations (15)-(18), we should discretize these equations further. Thus, the n is replaced by j, the time is denoted by n, and the step length of the space variable is denoted by h in the previous equations in this simulation. Then we can get the representations of the solutions from the given equations. The concrete equations for finding the solutions are described in Appendix A.
From Equations (15)- (19) using the given values of parameters, we can determine the changes of the solution of Equations (15)- (19) numerically, aptqr n and aptqi n , and of the corresponding |a n ptq| 2 , which denotes the probability of appearance of the soliton related to its amplitude, with varying time and space by the fourth-order Runge-Kutta method [107,108]. In this calculation, we choose the time step ∆t to be between 0.2 and 0.5pM{wq 1{2 and the initial excitation of a n p0q " Asechrpnń 0 qˆpχ 1`χ2 q{4Jws (where A is normalization constant) and the initial states of amino acid residues, q n (0) = π n (0) = 0 at the site n; where the molecular chain is fixed, the number of amino acid residues, N, is chosen to be 400 and a time step of 0.0195 is used. This means that 800 coupled differential equations with first-order-of-time derivation had to be solved in such a case. At the same time, the system of units eV for energy, 0 A for length, and ps for time have been used to be suitable for the numerical solutions of Equations (15)- (19). Otherwise, by using numerical simulations obtained by the fourth-order Runge-Kutta method [107,108], we could also satisfy the total energy of the soliton, E ď Φ(t)|H|Φ(t) ě constant, in which a possible imaginary part of the energy, which could occur because of numerical inaccuracies, is zero. Its accuracy is 0.001 feV; when the soliton is in motion, its probability of appearance must be normalized at all times, i.e., the number of particles for the system must be conservative, ř n |a n ptq| 2 " ř n |a n p0q| 2 " 1 up to 0.3 ppm. Meanwhile, the initial excitation conditions described previously are also required in this calculation. Total numerical simulation is performed through data-paralleling algorithms and MATLAB language. On the basis of the methods and conditions given, we obtain the representation for the change of |a n | 2 =|a j | 2 , which is the probability of the soliton occurred at the nth (or jth) amino acid molecule, with varying positions and times under the actions of different EFs.
We first calculate the solutions of Equations (15)- (19) numerically in the uniform and periodic proteins with a single chain using the above-average values for these parameters, i.e., M = 1.17ˆ10´2 5 kg; J = 1.55ˆ10´2 2 J; or = 9.68ˆ10´4 eV; ε 0 = 0.20 5 eV; χ 1 = 6.2ˆ10´1 1 N; w = (13´19.5) N/m; and χ 2 " p11´18qˆ10´1 2 N, which are widely accepted in investigations of energy transport in α-helical protein molecules [24][25][26][27][28][29][30][31][32][33][65][66][67][68][69][70][71][72][73][74]. If the fourth-order Runge-Kutta method is used, then the results for the state of transport of the soliton in the time-place are shown in Figures 2  and 3, respectively, in which we show the movements or progress of the soliton transporting the energy (in other words, its probability |a n qtq| 2 ) with varying time and space, where the symbols of three axes and their physical significances are marked clearly; where n is the number of amino acid residues the soliton has transported; and t is time of movement of the soliton. If the probabilities or amplitudes of the soliton, obtained from this calculation, are always constant in the movement process, then we can affirm that the solution of the dynamic equations is a soliton. Otherwise, these equations do not provide the soliton solution, which indicates that the soliton was dispersed or damped and that the energy transport was also dispersed or dissipated due to the influence of an externally applied field, according to the soliton theory [9,27,36]. Therefore, we can study and confirm the influences of EFs or EMFs on the soliton or energy transport in the protein molecules through the results obtained from the numerical simulations. Figures 2 and 3 indicate that this solution is a soliton because its probabilities or amplitudes are invariable in the propagation processes of the long-time period of 250 ps and the long distances of 400 amino acid residues, which are shown in Figure 2. Figure 3 exhibits the collision property of two solitons and also shows that the probabilities or amplitudes of the two solitons are invariable after the collision, resembling a feature of classical particles. Therefore, Figures 2 and 3 clearly show that the numerical solutions of Equations (15)- (19) are reliable and indicate the presence of a soliton according to the soliton theory [9,27,34].Thus, we can affirm that Equations (15)- (19) have the soliton solution, which is exactly the carrier of the energy transport in uniform and periodic protein molecules. are shown in Figure 2. Figure 3 exhibits the collision property of two solitons and also shows that the probabilities or amplitudes of the two solitons are invariable after the collision, resembling a feature of classical particles. Therefore, Figures 2 and 3 clearly show that the numerical solutions of Equations (15)- (19) are reliable and indicate the presence of a soliton according to the soliton theory [9,27,34].Thus, we can affirm that Equations (15)- (19) have the soliton solution, which is exactly the carrier of the energy transport in uniform and periodic protein molecules.  When protein molecules are exposed to EFs, the states of the soliton will change because of variations in the dipole-dipole interaction J, which would be replaced by at 0 0 = θ , as shown previously. As such, studies should explore further the properties of various states of the soliton. We studied these variations by using Equations (15)- (19) and the fourth-order processes of the long-time period of 250 ps and the long distances of 400 amino acid residues, which are shown in Figure 2. Figure 3 exhibits the collision property of two solitons and also shows that the probabilities or amplitudes of the two solitons are invariable after the collision, resembling a feature of classical particles. Therefore, Figures 2 and 3 clearly show that the numerical solutions of Equations (15)- (19) are reliable and indicate the presence of a soliton according to the soliton theory [9,27,34].Thus, we can affirm that Equations (15)- (19) have the soliton solution, which is exactly the carrier of the energy transport in uniform and periodic protein molecules.  When protein molecules are exposed to EFs, the states of the soliton will change because of variations in the dipole-dipole interaction J, which would be replaced by at 0 0 = θ , as shown previously. As such, studies should explore further the properties of various states of the soliton. We studied these variations by using Equations (15)- (19) and the fourth-order When protein molecules are exposed to EFs, the states of the soliton will change because of variations in the dipole-dipole interaction J, which would be replaced by J`Ñ E . Ñ p " J`ˇˇˇˇÑ EˇˇˇˇˇˇÑ pˇǎ t θ " 0 0 , as shown previously. As such, studies should explore further the properties of various states of the soliton. We studied these variations by using Equations (15)- (19) and the fourth-order Runge-Kutta method [107,108] under the influences of different strengths of EFs withˇˇˇˇÑ Eˇˇˇˇ= 25,500, 51,000, 76,500, and 102,000 V/m, which correspond to the variations of the dipole-dipole interaction energy of ∆J " J´J = 5%, 10%, 15%, and 20%, respectively, between the neighboring amino acid residues arising from the EF, in which the electric dipole momentˇˇÑ pˇˇ= 10 debyes is used, in which other parameters used by their average values , that are M, w, ε 0 , χ 1 , and χ 2 . Figure 4 shows the states of movement of the soliton under these conditions. For instance, the range of variations increases with the increase of strength of the EFs, such that the soliton is stable atˇˇˇˇÑ Eˇˇˇˇ= 25,500 and 51,000 V/m because its amplitudes have not been altered during the transport processes. However, the soliton disperses atˇˇˇˇÑ Eˇˇˇˇ= 76,500 V/m, as shown in Figure 4c, and disperses significantly atˇˇˇˇÑ Eˇˇˇˇ= 102,000 V/m, as shown in Figure 4d. Also, Figure 4 shows that their amplitudes are changed or reduced in these cases and that small ripples occur around the solitons in Figure 4c,d. This occurrence indicates that the states of the soliton, that transport the bio-energy, are changed because of influences by the EFs. This finding suggests that the features of bio-energy are highly sensitive to changes in the dipole-dipole interaction or the states and properties of externally applied EFs.
increases with the increase of strength of the EFs, such that the soliton is stable at E = 25,500 and 51,000 V/m because its amplitudes have not been altered during the transport processes. However, the soliton disperses at E = 76,500 V/m, as shown in Figure 4c, and disperses significantly at E = 102,000 V/m, as shown in Figure 4d. Also, Figure 4 shows that their amplitudes are changed or reduced in these cases and that small ripples occur around the solitons in Figure 4c,d. This occurrence indicates that the states of the soliton, that transport the bio-energy, are changed because of influences by the EFs. This finding suggests that the features of bio-energy are highly sensitive to changes in the dipole-dipole interaction or the states and properties of externally applied EFs.
In this numerical simulation, we use the fourth-order Runge-Kutta method [107,108] with the initial condition of a nα p0q " A 1 sechrpnα´n 0 αqˆpχ 1`χ2 q{4Jws, where A' is the normalization factor. Meanwhile, we use M = 5.73ˆ10´2 5 kg, w = 39 N/m, ε 0 = 0.2055 eV, J = 9.68ˆ10´4 eV, χ 1 = 6.2ˆ10´1 1 N, χ 2 = (10´18)ˆ10´1 2 N, and L= 1.5 meV for the α-helix protein molecules in the simulation, which are widely accepted in energy-transport investigations [24][25][26][27][28][29]. The numerical solutions of Equations (22) and (23) and their variations for the protein molecules with three channels are shown in Figure 5a, which shows the motion behaviors of the solution, in which the previously given initial conditions are motivated simultaneously on the first ends of the three chains. From this figure, we see that their solution can always retain a bell shape while moving over a long distance of the spacing of 400 amino acid residues in the period of 40 ps along the molecular chains without the dispersion. This feature is similar to the numerical results in Figure 2 for the single-chain proteins. On the basis of the soliton theory [9,27,36], we can also affirm from this result that this solution of Equations (22) and (23) also indicate the existence of a soliton in this case. Therefore, Equations (22) and (23) can be used to simulate the dynamic properties of the soliton excited in the protein molecules with three channels under the action of EF.
We studied the collision property of two solitons, setting up from opposite ends of the three channels in the protein molecules by using the fourth-order Runge-Kutta method. The result is shown in Figure 5b, in which the initial conditions given simultaneously motivate the opposite ends of the three channels, where the initial two solitons separating the 100 amino acid spaces in each channel collide with each other after approximately 17 ps. After this collision, the two solitons in each channel go through each other without scattering to propagate toward and separately along the three chains, which is also similar to the rules of collision for two solitons in the numerical results shown in Figure 3 for single-chain proteins. In accordance with the soliton theory [9,27,34], we can also confirm from the two results in Figure 5 that this solution to Equations (22) and (23) is still a soliton. Thus, by using this method, we can study the changes to the properties of energy transported by the solitons in Equations (22) and (23)  shown in Figure 5b, in which the initial conditions given simultaneously motivate the opposite ends of the three channels, where the initial two solitons separating the 100 amino acid spaces in each channel collide with each other after approximately 17 ps. After this collision, the two solitons in each channel go through each other without scattering to propagate toward and separately along the three chains, which is also similar to the rules of collision for two solitons in the numerical results shown in Figure 3 for single-chain proteins. In accordance with the soliton theory [9,27,34], we can also confirm from the two results in Figure 5 that this solution to Equations (22) and (23) is still a soliton. Thus, by using this method, we can study the changes to the properties of energy transported by the solitons in Equations (22) and (23)  When the protein molecules with the three channels are acted on by an EF, the states of the soliton will also change because of variations in the dipole-dipole interaction J, as noted previously. We studied the changing states of movement of the soliton under the effects of EFs of different strengths by using Equations (22) and (23) and the fourth-order Runge-Kutta method [107,108], in When the protein molecules with the three channels are acted on by an EF, the states of the soliton will also change because of variations in the dipole-dipole interaction J, as noted previously. We studied the changing states of movement of the soliton under the effects of EFs of different strengths by using Equations (22) and (23) and the fourth-order Runge-Kutta method [107,108], in which the dipole-dipole interaction J is replaced by J`Ñ E . Ñ p " J`ˇˇˇˇÑ EˇˇˇˇˇˇÑ pˇˇat θ " 90 0 ; where J = 9.68ˆ10´4 eV andˇˇÑ pˇˇ" 10 debyes; other parameters are used as their average values, M, w, ε 0 , ε 0 , χ 1 , and χ 2 , as described. Figure 6 shows the states of movement of the solitons in the cases ofˇˇˇˇÑ Eˇˇˇˇ= 17,000, 25,500, and 34,000 V/m, which correspond to variations of the dipole-dipole interaction energy of ∆J " J´J = 3%, 5%, and 7% between the neighboring amino acid residues, respectively. From these figures, we see clearly that the states of the soliton vary with increasing EF strength. The results show that the soliton is stable atˇˇˇˇÑ Eˇˇˇˇ= 17,000 V/m but that it disperses atˇˇˇˇÑ Eˇˇˇˇ= 25,500 V/m and damps at 34,000 V/m. The stability of solitons in three-channel α-helix protein molecules is reduced with respect to the solitons in single-chain protein molecules under the influence of an EF through the existence of the disperse effect, which is caused by the chain-chain interactions among the three channels. Thus, we can surmise that the state of a soliton transporting energy is highly sensitive to the EFs in EMFs.
indicate the relevance of our investigations of the biological effects of EFs or EMFs because the strengths of electric field, as described previously, can generate corresponding electrically interactional energies to the amino acid residues within a few nanometers, which is approximately one order of magnitude smaller than that of the dipole-dipole interactional energy J between them. Therefore, these estimated values are basically the same as those in the calculations previously described, which, we believe, are basically correct.

Experimental Evidences for this Theory
The theoretical results and mechanisms of the biological effects of EMFs mentioned previously require further experimental confirmation. Experimental evidence should relate mainly to the affirmation of the nonlinear excitation or the soliton, which provides the bio-energy transport in an α-helix protein molecule, and its changes with varying EF in EMF. Investigations indicate that the nonlinear interaction, or the soliton, can be confirmed from the infrared spectrum of absorption from the α-helix protein molecule [113]. The essence of the experimental evidence is to confirm the real existence of the soliton in transporting the bio-energy within the α-helix protein molecules. Therefore, the first step of the investigation is to determine the existence of soliton excitation in the protein molecules with α-helix structures. Subsequently, investigation should focus on the exact influence of the EF on the soliton. By confirming the existence of a soliton and indicating that EMFs can influence or change the properties of the soliton, we can confirm the correctness of the previously given mechanics for the biological effects of EFs, that is, EMFs can influence the bio-energy transport in protein molecules. Careri et al. [113][114][115][116][117], Scott [30,31,[118][119][120][121], and Alexander et al. [122,123] measured the infrared spectrum of absorption for acetanilide (ACN) [113], a molecular structure of like-protein molecules. From these investigations, they confirmed that the 1666-and 1650-cm −1 bands that occurred in the infrared spectra represent the excited states of amide-I mode in ACN, which corresponds to the excitations of the exciton and the soliton, respectively. Their infrared absorption strengths change linearly and exponentially by increasing the temperatures in ACN, respectively. Marracino et al. [112] showed the effects of EFs in EMFs on chemical reactions and determined that a 100,000 V/m electric field can have a significant effect on chemical reaction in micelles. In such a case, we think that this result might serve somewhat as a basis of comparison. These comparisons indicate the relevance of our investigations of the biological effects of EFs or EMFs because the strengths of electric field, as described previously, can generate corresponding electrically interactional energies to the amino acid residues within a few nanometers, which is approximately one order of magnitude smaller than that of the dipole-dipole interactional energy J between them. Therefore, these estimated values are basically the same as those in the calculations previously described, which, we believe, are basically correct.

Experimental Evidences for this Theory
The theoretical results and mechanisms of the biological effects of EMFs mentioned previously require further experimental confirmation. Experimental evidence should relate mainly to the affirmation of the nonlinear excitation or the soliton, which provides the bio-energy transport in an α-helix protein molecule, and its changes with varying EF in EMF. Investigations indicate that the nonlinear interaction, or the soliton, can be confirmed from the infrared spectrum of absorption from the α-helix protein molecule [113]. The essence of the experimental evidence is to confirm the real existence of the soliton in transporting the bio-energy within the α-helix protein molecules. Therefore, the first step of the investigation is to determine the existence of soliton excitation in the protein molecules with α-helix structures. Subsequently, investigation should focus on the exact influence of the EF on the soliton. By confirming the existence of a soliton and indicating that EMFs can influence or change the properties of the soliton, we can confirm the correctness of the previously given mechanics for the biological effects of EFs, that is, EMFs can influence the bio-energy transport in protein molecules. Careri et al. [113][114][115][116][117], Scott [30,31,[118][119][120][121], and Alexander et al. [122,123] measured the infrared spectrum of absorption for acetanilide (ACN) [113], a molecular structure of like-protein molecules. From these investigations, they confirmed that the 1666-and 1650-cm´1 bands that occurred in the infrared spectra represent the excited states of amide-I mode in ACN, which corresponds to the excitations of the exciton and the soliton, respectively. Their infrared absorption strengths change linearly and exponentially by increasing the temperatures in ACN, respectively. Thus, the existence of the soliton in ACN was affirmed [115][116][117][118][119]. We chose a collagen protein with α-helix conformation to confirm the existence of the soliton and study the influence of EFs, and to affirm further the validity of the mechanism and theory previously described. These investigations are summarized as follows.

Experimental Evidence of the Existence of Solitons in Protein Molecules
Because collagen is a common biomacromolecule that widely exists in a solid-like or soft condensed state from 0 to 95˝C in living systems, it is a basic component of surface-tissue musculature, e.g., smooth muscle [113][114][115][116][117][118][119][120][121][122][123][124][125]. Tropocollagen, a kind of collagen, is an α-superhelical biopolymer with three channels (Figure 7), in which each α-peptide chain contains 1050 amino acid residues and a sugary side chain. The structure consists of approximately 35% glycine, 10% proline, and 9% hydroproline, as well as some alumine and hydrolysin. Its molecular structure might be described as follows. Its primary conformation is a chain of (Gly-X-Y) n , where X and Y are the proline and the hydroproline, respectively. Its secondary structure consists of main chains that are regularly folded into an α-helical chain with a left-hand spin, in which the glycine is at the center of the helix and the proline and hydroxylproline are at its exterior. If the three chains are again folded into a multi-helix structure, a tertiary conformation is constructed in which each chain contains several helices. When the three chains are assembled and wind further into a right-hand spin to construct a long fiber, in virtue of the hydrogen bonds and the C=O . . . NH groups, its quaternary conformation is obtained. Thus, the linkage between the two lysine residue side groups in the peptide chains occurs through covalent bonding, which is formed and produced by the oxidoreductase reaction of lysine, as represented by -CO-CH(NH)-CH 2 -CH 2 -CH 2 -CH=N-CH 2 -CH 2 -CH 2 -CH 2 -CH(NH)-CO-. Collagen is, therefore, a kind of α-helical protein molecule with an α-helix conformation (Figure 7), in which the hydrogen bonds and the covalent bonds are arranged alternately, which plays an important role in stabilizing the molecular structure of the collagen and completing its biological functions (such as the transport of energy and information) and enhancing its tensile strength. We used the features of the infrared spectrum to confirm the existence of the soliton and the influences of EMF and could thus verify further the validity of the previously given mechanism and theory.
in which the hydrogen bonds and the covalent bonds are arranged alternately, which plays an important role in stabilizing the molecular structure of the collagen and completing its biological functions (such as the transport of energy and information) and enhancing its tensile strength. We used the features of the infrared spectrum to confirm the existence of the soliton and the influences of EMF and could thus verify further the validity of the previously given mechanism and theory. We measured the infrared spectrum of absorption of the collagen, which was supplied by Sigma-Aldrich (St. Louis, MO, USA), using the spectrum GXFFIR spectrometer with a DTGS detector and a variable-temperature bath with a reported accuracy of ±1 °C, in which infrared silicon-carbon bars served as a light source, was provided by Perkin Elmer (Waltham, MA, USA), and a 670 Nicolet FT-IR spectrometer with a 4-cm −1 resolution and ATR full-reflective device with 16 successive scans were supplied by Nicolet Nexus (Ramsey, MN, USA), respectively. In this experiment we measured and collected mainly the infrared absorption spectrum of the collagen in the region of 1000-2000 cm −1 by using, first, the Perkin Elmer spectrum GXFFIR spectrometer [108][109][110][111][112][113][114][115][116][117][118][119][120][121][124][125][126] and, again, the 670 Nicolet FT-IR and the GXFFIR spectrometers with resolutions of 4 cm −1 to check the results. The measured samples of collagen were sandwiched between KBr windows at temperatures controlled by a variable-temperature bath. Throughout the experiment, we inspected the infrared absorption spectrum of the collagen in the range of 1000-4000 cm −1 . To obtain an acceptable signal-to-noise ratio, we took 16 scans of these infrared spectra. When the temperatures of the samples changed from 15 to 95 °C, variations of the infrared spectra of the collagen were also measured and inspected by a spectrometer. In this experiment, the water and the water vapor in these samples were completely drained from the tested samples so as not to influence results.
The infrared spectrum of collagen at 25 °C in the range of 1540-1710 cm −1 is shown in Figure 8. Its spectrum at high frequencies is dominated by the amide spectrum, in which three clear Amide I vibrational modes are present: 1666.11, 1680.38, and 1695.32 cm −1 , where an important feature in the spectrum is the appearance of a new band at 1650.01 cm −1 . Other amide bands, such as the We measured the infrared spectrum of absorption of the collagen, which was supplied by Sigma-Aldrich (St. Louis, MO, USA), using the spectrum GXFFIR spectrometer with a DTGS detector and a variable-temperature bath with a reported accuracy of˘1˝C, in which infrared silicon-carbon bars served as a light source, was provided by Perkin Elmer (Waltham, MA, USA), and a 670 Nicolet FT-IR spectrometer with a 4-cm´1 resolution and ATR full-reflective device with 16 successive scans were supplied by Nicolet Nexus (Ramsey, MN, USA), respectively. In this experiment we measured and collected mainly the infrared absorption spectrum of the collagen in the region of 1000-2000 cm´1 by using, first, the Perkin Elmer spectrum GXFFIR spectrometer [108][109][110][111][112][113][114][115][116][117][118][119][120][121][124][125][126] and, again, the 670 Nicolet FT-IR and the GXFFIR spectrometers with resolutions of 4 cm´1 to check the results. The measured samples of collagen were sandwiched between KBr windows at temperatures controlled by a variable-temperature bath. Throughout the experiment, we inspected the infrared absorption spectrum of the collagen in the range of 1000-4000 cm´1. To obtain an acceptable signal-to-noise ratio, we took 16 scans of these infrared spectra. When the temperatures of the samples changed from 15 to 95˝C, variations of the infrared spectra of the collagen were also measured and inspected by a spectrometer. In this experiment, the water and the water vapor in these samples were completely drained from the tested samples so as not to influence results.
The infrared spectrum of collagen at 25˝C in the range of 1540-1710 cm´1 is shown in Figure 8. Its spectrum at high frequencies is dominated by the amide spectrum, in which three clear Amide I vibrational modes are present: 1666.11, 1680.38, and 1695.32 cm´1, where an important feature in the spectrum is the appearance of a new band at 1650.01 cm´1. Other amide bands, such as the Amide II at 1624.9 cm´1, the Amide III at 1553.1, and so on are also observed and obtained from the previously designated.
The infrared spectrum of collagen denoted in Figure 8 has the following properties: (1) Amide II, III bands, in addition to an Amide I band, exist in the infrared absorption spectrum of the collagen; (2) a new band of 1650 cm´1 appears besides the conventional 1666 cm´1 band. These results are very interesting and important because from this experiment we confirm the appearances of 1650 cm´1 and 1666 cm´1,which correspond the excitations of soliton and exciton in the collagen in this case, respectively, in accordance with Careri et al. [113][114][115][116][117] and Scott [30,31,[118][119][120][121] researched conclusions. Thus we may use the method and material to investigate further the influences of EF in EMF on the soliton, or energy transport in collagen, which is described in detail as follows. Amide II at 1624.9 cm −1 , the Amide III at 1553.1, and so on are also observed and obtained from the previously designated. The infrared spectrum of collagen denoted in Figure 8 has the following properties: (1) Amide II, III bands, in addition to an Amide I band, exist in the infrared absorption spectrum of the collagen; (2) a new band of 1650 cm −1 appears besides the conventional 1666 cm −1 band. These results are very interesting and important because from this experiment we confirm the appearances of 1650 cm −1 and 1666 cm −1 ,which correspond the excitations of soliton and exciton in the collagen in this case, respectively, in accordance with Careri et al. [113][114][115][116][117] and Scott [30,31,[118][119][120][121] researched conclusions. Thus we may use the method and material to investigate further the influences of EF in EMF on the soliton, or energy transport in collagen, which is described in detail as follows.
At the same time, the infrared spectra of absorption of collagen at different temperatures are carefully recorded, where the temperatures of the samples varied from 15 to 95 °C, with the intervals of 10 °C and a reported accuracy of ±1 °C, which is controlled by variable-temperature equipment. The temperature-dependent intensity of infrared absorption at temperatures ranging from 15 to 95 °C is shown in Figure 9. Figure 9 clearly indicates that the intensity of the 1650-cm −1 band increases with the decrease of temperature, without any apparent change in frequency and shape, but is weakened at 95 °C. However, the strength of the Amide I infrared absorption of the 1666-cm −1 band decreases with a decrease in temperature.  At the same time, the infrared spectra of absorption of collagen at different temperatures are carefully recorded, where the temperatures of the samples varied from 15 to 95˝C, with the intervals of 10˝C and a reported accuracy of˘1˝C, which is controlled by variable-temperature equipment. The temperature-dependent intensity of infrared absorption at temperatures ranging from 15 to 95˝C is shown in Figure 9. Figure 9 clearly indicates that the intensity of the 1650-cm´1 band increases with the decrease of temperature, without any apparent change in frequency and shape, but is weakened at 95˝C. However, the strength of the Amide I infrared absorption of the 1666-cm´1 band decreases with a decrease in temperature. Amide II at 1624.9 cm −1 , the Amide III at 1553.1, and so on are also observed and obtained from the previously designated. The infrared spectrum of collagen denoted in Figure 8 has the following properties: (1) Amide II, III bands, in addition to an Amide I band, exist in the infrared absorption spectrum of the collagen; (2) a new band of 1650 cm −1 appears besides the conventional 1666 cm −1 band. These results are very interesting and important because from this experiment we confirm the appearances of 1650 cm −1 and 1666 cm −1 ,which correspond the excitations of soliton and exciton in the collagen in this case, respectively, in accordance with Careri et al. [113][114][115][116][117] and Scott [30,31,[118][119][120][121] researched conclusions. Thus we may use the method and material to investigate further the influences of EF in EMF on the soliton, or energy transport in collagen, which is described in detail as follows.
At the same time, the infrared spectra of absorption of collagen at different temperatures are carefully recorded, where the temperatures of the samples varied from 15 to 95 °C, with the intervals of 10 °C and a reported accuracy of ±1 °C, which is controlled by variable-temperature equipment. The temperature-dependent intensity of infrared absorption at temperatures ranging from 15 to 95 °C is shown in Figure 9. Figure 9 clearly indicates that the intensity of the 1650-cm −1 band increases with the decrease of temperature, without any apparent change in frequency and shape, but is weakened at 95 °C. However, the strength of the Amide I infrared absorption of the 1666-cm −1 band decreases with a decrease in temperature.  The peak strengths of the 1650-and 1666-cm´1 bands versus temperature are clearly shown in Figure 9. Although the temperature dependences of strengths for the 1650-and 1666-cm´1 bands differ, the absorption strength of the 1666-cm´1 band increases linearly with increasing temperature; however, the strength of the 1650-cm´1 band decreases exponentially with an increase in temperature This relationship can be approximately simulated in the exponential form of I " I 0 exp t´r0.437`8.987ˆ10 -6 pT{˝Cq 2 su (24) where T is represented in Celsius temperature (˝C); and I 0 is a constant related to the initial strength. Figure 10 shows the experimental data of the relationship between the logarithm of relative strength, Ln (I/I 0 ), and the temperature, (T/˝C) 2 , for the 1650-cm´1 peak in the collagen, and this relationship can be simulated by Ln (I/I 0 ) =´(0.437 + 8.987ˆ10´6 (T/˝C) 2 ). It denotes that the strengths of the 1650-cm´1 peak in the infrared spectrum of absorption in the collagen decreased exponentially with increasing temperatures. Figure 10b gives the linear changed relation of relative strength, I/I 0 , for the 1666-cm´1 peak changing with increasing temperatures in the infrared spectrum of the absorption of the collagen molecule, a novel and interesting finding. The peak strengths of the 1650-and 1666-cm −1 bands versus temperature are clearly shown in Figure 9. Although the temperature dependences of strengths for the 1650-and 1666-cm −1 bands differ, the absorption strength of the 1666-cm −1 band increases linearly with increasing temperature; however, the strength of the 1650-cm −1 band decreases exponentially with an increase in temperature This relationship can be approximately simulated in the exponential form of I = I0 exp {−[0.437 + 8.987 × 10 −6 (T/°C) 2 ]} (24) where T is represented in Celsius temperature (°C); and I0 is a constant related to the initial strength. Figure 10 shows the experimental data of the relationship between the logarithm of relative strength, Ln (I/I0), and the temperature, (T/°C) 2 , for the 1650-cm −1 peak in the collagen, and this relationship can be simulated by Ln (I/I0) = −(0.437 + 8.987 × 10 −6 (T/°C) 2 ). It denotes that the strengths of the 1650-cm −1 peak in the infrared spectrum of absorption in the collagen decreased exponentially with increasing temperatures. Figure 10b gives the linear changed relation of relative strength, I/I0, for the 1666-cm −1 peak changing with increasing temperatures in the infrared spectrum of the absorption of the collagen molecule, a novel and interesting finding. On the basis of the experimental results in Figures 8-10, we can affirm again the existence of solitons in collagen according to the theory and conclusions of Careri et al. [113][114][115][116][117], Scott [30,31,[108][109][110][111][112][113][114][115][116][117][118][119][120][121], and Alexander et al. [122,123]. Taking into account these results, we studied further the influences of EMFs on solitons, because the outline and features of solitons depend closely on the dipole-dipole interaction between the neighboring amino acid residues with certain electric dipole moments in the collagen, as mentioned previously.

Experimental Evidence of the Influence of EFs on Solitons in Protein Molecules
We investigated and measured further the influence of EMFs on the properties of the 1650-cm −1 peak in collagen, which is related to the soliton excited, by using the 670 Nicolet FT-IR spectrometer based on the above results. If we obtain the changes of features of the peak of 1650 cm −1 by varying the externally applied electric voltage, we can affirm that the EF varies the states and features of the soliton, or the bio-energy transport, in the collagen because of variations of the electric dipole moments of amino acids in the protein.
In this experiment the electric voltage, which is applied on the collagen, changed from 15,000 to 20,000 V, and our results are shown in Figures 11-14. Figure 11 shows the results of the infrared spectrum of absorption of collagen in 480-2000 cm −1 , and it is apparent that the 1650 cm −1 peak occurs in this case. Hence, the soliton exists within the collagen in this case. In this experiment, the externally applied voltages are linked on the two sides of the thin plates, where the measured collagen samples are first made into powders, which are again inserted into the center of the thin plates and were kept On the basis of the experimental results in Figures 8-10, we can affirm again the existence of solitons in collagen according to the theory and conclusions of Careri et al. [113][114][115][116][117], Scott [30,31,[108][109][110][111][112][113][114][115][116][117][118][119][120][121], and Alexander et al. [122,123]. Taking into account these results, we studied further the influences of EMFs on solitons, because the outline and features of solitons depend closely on the dipole-dipole interaction between the neighboring amino acid residues with certain electric dipole moments in the collagen, as mentioned previously.

Experimental Evidence of the Influence of EFs on Solitons in Protein Molecules
We investigated and measured further the influence of EMFs on the properties of the 1650-cm´1 peak in collagen, which is related to the soliton excited, by using the 670 Nicolet FT-IR spectrometer based on the above results. If we obtain the changes of features of the peak of 1650 cm´1 by varying the externally applied electric voltage, we can affirm that the EF varies the states and features of the soliton, or the bio-energy transport, in the collagen because of variations of the electric dipole moments of amino acids in the protein.
In this experiment the electric voltage, which is applied on the collagen, changed from 15,000 to 20,000 V, and our results are shown in Figures 11-14. Figure 11 shows the results of the infrared spectrum of absorption of collagen in 480-2000 cm´1, and it is apparent that the 1650 cm´1 peak occurs in this case. Hence, the soliton exists within the collagen in this case. In this experiment, the externally applied voltages are linked on the two sides of the thin plates, where the measured collagen samples are first made into powders, which are again inserted into the center of the thin plates and were kept in tight contact with each other through extrusion and pressurization. Subsequently, they were placed in the sample bath in the FT-IR spectrometer. When the electric circuit built was connected, the collagen samples were exposed to the externally applied EF. We could then measure and record the infrared spectra of the collagen samples and their variations (Figures 12-14). Figure 12 shows the decreases in the peak height with increases in EF strength up to 30,000 V/m, and the inverse beyond 30,000 V/m. Figures 13 and 14 indicate changes in the position of the 1650-cm´1 peak and its half-peak width with increases to the externally applied EF strength. Figure 13 clearly shows that the position of this peak lifts with increases to the EF strength below 30,000 V/m but lowers sluggishly after 30,000 V/m. For its half-peak width, we find that the half-peak width decreases with increases to the EF strength up to 30,000 V/m, but it increases sluggishly after 30,000 V/m ( Figure 14). Therefore, Figures 12-14 exhibit clearly that the features of the 1650-cm´1 peak vary under the influences of externally applied electric voltages, which are basically consistent with those in Figure 6. Thus, we can affirm and judge from these experiments that the externally applied EFs influence and change the states and features of the soliton and the bio-energy transport in the collagen (Figures 12-14). Hence, we verified the validity of these results and the mechanism of influence of EMFs and EFs on life bodies through changes of the bio-energy transport or in the features of the soliton caused by changes to the electric properties of amino acids in protein molecules.
Int. J. Mol. Sci. 2016, 17, 1130 in tight contact with each other through extrusion and pressurization. Subsequently, they were placed in the sample bath in the FT-IR spectrometer. When the electric circuit built was connected, the collagen samples were exposed to the externally applied EF. We could then measure and record the infrared spectra of the collagen samples and their variations (Figures 12-14). Figure 12 shows the decreases in the peak height with increases in EF strength up to 30,000 V/m, and the inverse beyond 30,000 V/m. Figures 13 and 14 indicate changes in the position of the 1650-cm −1 peak and its half-peak width with increases to the externally applied EF strength. Figure 13 clearly shows that the position of this peak lifts with increases to the EF strength below 30,000 V/m but lowers sluggishly after 30,000 V/m. For its half-peak width, we find that the half-peak width decreases with increases to the EF strength up to 30,000 V/m, but it increases sluggishly after 30,000 V/m ( Figure 14). Therefore, Figures  12-14 exhibit clearly that the features of the 1650-cm −1 peak vary under the influences of externally applied electric voltages, which are basically consistent with those in Figure 6. Thus, we can affirm and judge from these experiments that the externally applied EFs influence and change the states and features of the soliton and the bio-energy transport in the collagen (Figures 12-14). Hence, we verified the validity of these results and the mechanism of influence of EMFs and EFs on life bodies through changes of the bio-energy transport or in the features of the soliton caused by changes to the electric properties of amino acids in protein molecules.   in tight contact with each other through extrusion and pressurization. Subsequently, they were placed in the sample bath in the FT-IR spectrometer. When the electric circuit built was connected, the collagen samples were exposed to the externally applied EF. We could then measure and record the infrared spectra of the collagen samples and their variations (Figures 12-14). Figure 12 shows the decreases in the peak height with increases in EF strength up to 30,000 V/m, and the inverse beyond 30,000 V/m. Figures 13 and 14 indicate changes in the position of the 1650-cm −1 peak and its half-peak width with increases to the externally applied EF strength. Figure 13 clearly shows that the position of this peak lifts with increases to the EF strength below 30,000 V/m but lowers sluggishly after 30,000 V/m. For its half-peak width, we find that the half-peak width decreases with increases to the EF strength up to 30,000 V/m, but it increases sluggishly after 30,000 V/m ( Figure 14). Therefore, Figures  12-14 exhibit clearly that the features of the 1650-cm −1 peak vary under the influences of externally applied electric voltages, which are basically consistent with those in Figure 6. Thus, we can affirm and judge from these experiments that the externally applied EFs influence and change the states and features of the soliton and the bio-energy transport in the collagen (Figures 12-14). Hence, we verified the validity of these results and the mechanism of influence of EMFs and EFs on life bodies through changes of the bio-energy transport or in the features of the soliton caused by changes to the electric properties of amino acids in protein molecules.     Changes to the width of a half-peak for a 1650-cm −1 peak with increases to the externally applied voltage.

Conclusions
The mechanism of influence of an externally applied EF or EMF on bio-energy transport along a protein molecule and their properties are investigated by using analytic methods and numerical simulation as well as experimental measurements. Energy, released by the hydrolysis reaction of an ATP molecule, and its transport are the basic biological activities in the life system, including muscle contraction, DNA duplication, neuroelectric pulse transfer along the neurolemma, and calcium and sodium pumping. This transport is conducted through movement of the soliton along the protein molecules and the dipole-dipole interaction between neighbor amino acid residues. Thus, we affirmed that externally applied EFs or EMFs can directly change the strength and direction of electric dipole moments of amino acid residues in protein molecules, which results directly in variations of property of the bio-energy transported by the soliton, thus a series of new biological effects can occur in living systems. In this investigation we obtained the following results.
The first result confirms one mechanism of the biological effect of EMFs and EFs, that is, the EFs and EMFs changed the features of bio energy transport through variations of electric dipole moments of the amino acid residue or the dipole-dipole interactional energy between the neighboring amino acid residues. Thus, we determined first that the targets of the biological effects of EMFs and EFs are the amino acid residues in protein molecules.
The second result establishes a theory of the biological effects of EMFs and EFs through the study of soliton mechanisms and bio-energy transport in protein molecules. In this theory, the  Changes to the width of a half-peak for a 1650-cm −1 peak with increases to the externally applied voltage.

Conclusions
The mechanism of influence of an externally applied EF or EMF on bio-energy transport along a protein molecule and their properties are investigated by using analytic methods and numerical simulation as well as experimental measurements. Energy, released by the hydrolysis reaction of an ATP molecule, and its transport are the basic biological activities in the life system, including muscle contraction, DNA duplication, neuroelectric pulse transfer along the neurolemma, and calcium and sodium pumping. This transport is conducted through movement of the soliton along the protein molecules and the dipole-dipole interaction between neighbor amino acid residues. Thus, we affirmed that externally applied EFs or EMFs can directly change the strength and direction of electric dipole moments of amino acid residues in protein molecules, which results directly in variations of property of the bio-energy transported by the soliton, thus a series of new biological effects can occur in living systems. In this investigation we obtained the following results.
The first result confirms one mechanism of the biological effect of EMFs and EFs, that is, the EFs and EMFs changed the features of bio energy transport through variations of electric dipole moments of the amino acid residue or the dipole-dipole interactional energy between the neighboring amino acid residues. Thus, we determined first that the targets of the biological effects of EMFs and EFs are the amino acid residues in protein molecules.
The second result establishes a theory of the biological effects of EMFs and EFs through the study of soliton mechanisms and bio-energy transport in protein molecules. In this theory, the

Conclusions
The mechanism of influence of an externally applied EF or EMF on bio-energy transport along a protein molecule and their properties are investigated by using analytic methods and numerical simulation as well as experimental measurements. Energy, released by the hydrolysis reaction of an ATP molecule, and its transport are the basic biological activities in the life system, including muscle contraction, DNA duplication, neuroelectric pulse transfer along the neurolemma, and calcium and sodium pumping. This transport is conducted through movement of the soliton along the protein molecules and the dipole-dipole interaction between neighbor amino acid residues. Thus, we affirmed that externally applied EFs or EMFs can directly change the strength and direction of electric dipole moments of amino acid residues in protein molecules, which results directly in variations of property of the bio-energy transported by the soliton, thus a series of new biological effects can occur in living systems. In this investigation we obtained the following results.
The first result confirms one mechanism of the biological effect of EMFs and EFs, that is, the EFs and EMFs changed the features of bio energy transport through variations of electric dipole moments of the amino acid residue or the dipole-dipole interactional energy between the neighboring amino acid residues. Thus, we determined first that the targets of the biological effects of EMFs and EFs are the amino acid residues in protein molecules.
The second result establishes a theory of the biological effects of EMFs and EFs through the study of soliton mechanisms and bio-energy transport in protein molecules. In this theory, the dipole-dipole interactional energy J between neighboring amino acid residues is replaced by J`Ñ E . Ñ p " J`ˇˇˇˇÑ EˇˇˇˇˇˇÑ pˇˇcosθ. We found changes to the features of bio-energy transport in the Pang's soliton under different EFs using analytical, numerical, and experimental methods.
Our analytical investigations indicated that EMFs and EFs influenced both the amplitude of the soliton as well as its form and outline. Its degree of effect depended on the strength of the EF, as well as its direction with respect to the electric dipole moments of amino acid residues. Thus, the biological effects of EFs and EMFs decrease with increases to the angle θ between them. If θ " 0 0 , then Ñ E . Ñ p "ˇˇˇˇÑ EˇˇˇˇˇˇÑ pˇˇ, but Ñ E . Ñ p " 0 at θ " 90 0 , therefore, EMFs and EFs have no biological effect in this case, even at high strengths. Hence, different EMFs and EFs have different biological effects as well, whereas the same EMFs and EFs also have different biological effects because their effects are determined by both the strength and the direction of the EMF or the EF. However, this effect has obvious randomness at the macroscopic level.
In the numerical simulation, we found that EMFs and EFs change the states and features of the movement of the soliton by using the fourth-order Runge-Kutta numerical simulation method in Pang's model [89][90][91][92][93][94][95][96][97][98][99][100][101][102][103][104][105]. By applying this simulation, we found that the solitons are stable at EF strengths of 25,500 and 51,000 V/m, but the solitons disperse at 76,500 V/m and significantly so at 102,000 V/m for protein molecules with single chains. However, the soliton is stable only at the EF strength of 17,000 V/m. It disperses at 25,500 V/m and damps at 34,000 V/m in protein molecules with three channels. This finding implies that the stability of a soliton in a three-channel α-helix protein molecule is reduced compared with that in a single-chain protein molecule because of the influence of the dispersed effect of chain-chain interactions among the three channels under the actions of EFs and EMFs. At the same time, we found that the capability of the bio-energy transported by the soliton is depressed considerably if the strength of the EMF or the EF is very high. These results verified that EMFs and EFs can vary the properties of the energy transport and the states of the soliton and that other biological effects could occur in living systems.
Finally we inspected and confirmed experimentally the real existence of the soliton, which corresponds to the excitation of a 1650-cm´1 peak, by using the infrared spectra of absorption of collagen with an α-helical structure and its change of features with varying external EFs, in which the height, position, and half-peak width of the soliton vary with different external EFs.
However, we should note the following problems associated with these investigations.
(1) We have not studied concretely the macroscopic biological effects arising from changes to energy transport under the influences of EMFs and EFs. Thus, further research into the molecular and cellular biology must be pursued. Thus, we cannot confirm that the concrete biological effects arising from variations in the bio-energy transport in protein molecules in the presence of EMFs or EFs are advantageous or harmful to the health of humans and animals. (2) The mechanism of the biological effects of EFs and EMFs that we propose require additional experimental confirmation, i.e., through instrumentation and novel methods, taking direct measurements of any changed features in the electric dipole moments of α-helix proteins or of the amino acid residues in them arising from EMFs and EFs, which entails many challenges. (3) Because the strength of EMF frequency is altered in a real EMF, we studied the biological effects of static EFs or DC fields only. Therefore, further investigations would focus on the biological effects of AC fields or altered EFs. and considering further the neighboring interactions among the three channels we can determine i . a nα ptq " ε 0 a nα ptq´Jra n`1α ptq`a n´1α ptqs`χ 1 rq n`1α ptq´q n´1α ptqsa nα ptq χ 2 rq n`1α ptq´q n´1α ptqsra n`1α ptq`a n´1α ptqs q nα " Wrq n`1α ptq´2q nα ptq`q n´1α ptqs`2χ 1 r|a n`1α ptq| 2´| a n´1α ptq| 2 s 2χ 2 tan α ptqra n`1α ptq´a n´1α ptqs`a nα ptqra˚n`1 α ptq´a˚n´1 α ptqsu.