Next Article in Journal
Expression Analysis of PIN Genes in Root Tips and Nodules of Medicago truncatula
Previous Article in Journal
Comparative Analysis for the Presence of IgG Anti-Aquaporin-1 in Patients with NMO-Spectrum Disorders
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Institute of Physical Electrons, University of Electronic Science and Technology of China, Chengdu 610054, China
2
Department of Physics, East China Normal University, Shanghai 200062, China
3
State Key Laboratory of Electrical Insulation and Power Equipment, Xi’an Jiaotong University, Xi’an 710049, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2016, 17(8), 1130; https://doi.org/10.3390/ijms17081130
Submission received: 16 January 2016 / Revised: 25 June 2016 / Accepted: 27 June 2016 / Published: 25 July 2016
(This article belongs to the Section Physical Chemistry, Theoretical and Computational Chemistry)

Abstract

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

1. 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 [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23]. 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.

2. Theories of Energy Transport in Protein Molecules and the Mechanism of Influence of EFs on Transport

2.1. Davydov’s Theory of Energy Transport and Its Features

Energy in living systems comes mainly from the hydrolysis reaction of ATP molecules, which is denoted [9,24] by
ATP 4 + H 2 O ADP 3 + HPO 4 2 + H + + 0.43 e V
where ADP is adenosine diphosphate. The energy of 0.43 eV released from this chemical reaction is used in many life activities and processes, such as muscle contraction, neuroelectric pulse transfer along the neurolemma, DNA reduplication, and calcium and sodium pumping in cell membranes, which are some basic life processes and activities in life systems. This function suggests clearly that the energy released in this reaction and its transport is related closely to the growth and development of human beings and animals. Thus, most life processes require the release of energy by ATP hydrolysis, with protein molecules playing a key role in energy transport.
As it is known, the amino acid molecule consists of an amino group (NH2), a carboxyl group (COOH), and a side group or chain functional group attached to an α-carbon atom, they polymerize to form long chains of residues of …H–N–C=O…H–N–C=O…H–N–C=O…H–N–C=O…, where the dotted lines indicate the hydrogen bond, that constitute the chain structure of protein molecules. When the polypeptide chains have been formed, they can further fold into an α-helix, a β-sheet, and a globular conformation. In an α-helix structure, three chains of hydrogen-bonded peptide groups can be represented approximately along the horizontal direction of the sequence. In this process, the distributions of the positive and negative charges in each amino acid residue change, and most of the amino acid residues are polarized and have certain electric dipole moments. Figure 1 shows a typical structure of an α-helix protein molecule [24,25,26,27,28,29,30,31,32,33], in which the forms of the three polypeptide chains built by hydrogen bonds and amide-I bonds are clearly exhibited.
However, questions remain as to how bio-energy is transported into life systems. Although several decades of research have been amassed, in 1971 Davydov proposed the soliton theory. According to this theory, bio-energy transports along the α-helical protein molecules, for which the structure is denoted in Figure 1 [24,25,26,27,28,29,30,31,32,33]. The mechanism of the bio-energy transport can be described as follows. The stretching vibration of C=O bonds(Amide Is) in the amino acid residues in the proteins is affected by energy released from ATP hydrolysis, in which the vibrational quanta, called excitons, and the deformation of amino acid residues also occur simultaneously and influence each other in this case. Thus, the nonlinear coupling interaction between the excitons and the deformation of amino acid residues, that occur in this process, make the excitons self-trap as a soliton, which transports further along the protein molecules in virtue of the dipole–dipole interaction between neighboring amino acid residues. This is a well-known soliton mechanism of energy transport in life systems. Davydov’s theory works in α-helix protein molecules with three channels in Figure 1 [24,25,26,27,28,29,30,31,32,33]. In this case, a biological process of energy transport always occurs from the produced place to the absorbed place in the living system, which is carried out by means of protein molecules such as collagen, myosin, myoglobin, and actin.
According to the mechanism of bio-energy transport, Davydov established the theory of bio-energy transport in protein molecules [24,25,26,27,28,29,30,31,32,33], in which Davydov gave the Hamiltonian of the protein molecules, which was represented by
H D = n [ ε 0 B n + B n J ( B n + B n + B n B n + ) ] + n [ P n 2 2 M + 1 2 w ( u n u n 1 ) 2 ] + n [ χ 1 ( u n + 1 u n 1 ) B n + B n ] = H e x + H p h + H int
where B n + ( B n ) is the creation (annihilation) operator for an Amide I quantum (exciton) in the site n, un is the displacement operator of amino acid residue at site n, Pn is its conjugate momentum operator, M is the mass of an amino acid residue and M = 1.17 × 10−25 kg in single-protein molecules or 5.73 × 10−25 kg in α-helix proteins with three channels, w is the elastic constant of the protein molecular chains and w = (13–19.5) N/m for single-protein molecules or (39–58.5) N/m for α-helix proteins, χ 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 r 0 = 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
D 2 ( t ) = | φ D | β ( t ) = n φ n ( t ) B n + exp { i n [ β n ( t ) P n π n ( t ) u n ] } | 0
where | 0 = | 0 e x | 0 p h , | 0 e x and | 0 p h 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
φ D ( x , t ) = ( u D 2 ) 1 / 2  seach [ μ D r 0 ( x x 0 v t ) ] exp { i [ v 2 J r 0 2 ( x x 0 ) E v t / ] }
If, | φ D ( t ) = n φ n ( t ) B n + | 0 e x in Equation (2) is an eigenstate of the number operator, N = n B n + B n , then the soliton contains only one exciton because N = φ D ( t ) | N | φ D ( t ) = 1 , i.e., N | φ D ( t ) = 1 | φ D ( t ) . 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 / ( 1 s 2 ) w J , s 2 = v 2 / v 0 2 , and v 0 = r 0 ( w / M ) 1 / 2 ; where v 0 is the sound speed in the protein molecular chains; v is the velocity of the soliton; and G D = 4 J μ 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 B D = χ 1 4 / 3 J w 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 [34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73]. 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 [34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71]. 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−12–10−13 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.

2.2. 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. [39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74], 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 [75,76,77,78,79,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,100,101,102]. In this case, the representations in Equations (1) and (2) for the single-channel protein molecules are replaced by
| Φ ( t ) = | α ( t ) | β ( t ) = 1 λ [ 1 + n α n ( t ) B n + + 1 2 ! ( n α n ( t ) B b 2 ) 2 ] | 0 e x × exp { i n [ β n ( t ) P n π n u n | 0 p h
and
H = H e x + H p h + H int = n [ ε 0 B n + B n J ( B n + B n + B n B n + ) ] + n [ P n 2 2 M + 1 2 w ( u n u n 1 ) 2 ] + n [ χ 1 ( u n + 1 u n 1 ) B n + B n + χ 2 ( u n + 1 u n ) ( B n + 1 + B n + B n + B n + 1 ) ]
respectively; where the created and annihilated operators of the exciton are represented by B n + and B n , respectively; | 0 e x and | 0 p h are the ground states of the exciton and phonon, respectively; and un and Pn are the displacement and momentum operators of amino acid residue at the site n, respectively. The β n ( t ) = Φ ( t ) | u n | Φ ( t ) , and π n ( t ) = Φ ( t ) | P n | Φ ( t ) are two sets of unknown functions; and λ is a normalization constant. Present nonlinear coupling constants are χ 1 and χ 2 = (10−15) × 10−11 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 ( t ) = ε 0 α n ( t ) J [ α n + 1 ( t ) + α n 1 ( t ) ] + χ 1 [ q n + 1 ( t ) q n 1 ( t ) ] α n ( t ) + χ 2 [ q n + 1 ( t ) q n ( t ) ] [ α n + 1 ( t ) + α n 1 ( t ) ] + 5 2 { w ( t ) 1 2 m q m ( t ) π m ( t ) π . m ( t ) q . m ( t ) ] } α n ( t )
M q .. n ( t ) = W [ q n + 1 ( t ) 2 q n ( t ) + q n 1 ( t ) ] + 2 χ 1 [ | α n + 1 | 2 | α n 1 ( t ) | ] + 2 χ 2 { α n * ( t ) [ α n + 1 ( t ) α n 1 ( t ) ] + α n ( t ) [ α n + 1 * ( t ) α n 1 * ( t ) ] }
In the continuum approximation, we get from Equations (6) and (7) [73,74,75,76,77,78,79,80,81,82,83,84]
i t α ( x , t ) = R ( t ) α ( x , t ) J r 0 2 2 x 2 α ( x , t ) G p | α ( x , t ) | 2 α ( x , t )
and
M 2 β ( x , t ) t 2 w r o 2 2 β ( x , t ) x 2 = 4 ( χ 1 + χ 1 ) r 0 x | α ( x , t ) | 2
where R ( t ) = = ε 0 2 J + 5 2 { W ( t ) 1 2 m [ β . n ( t ) π m ( t ) π . n ( t ) β ( t ) ] } .
The soliton solution of Equation (8) is denoted by
α ( x , t ) = ( μ p 2 ) 1 / 2 sec h [ ( μ p / r 0 ) ( x x 0 v t ) ] × exp { i [ v 2 J r 0 2 ( x x 0 ) E v t ] }
With μ P = 2 ( χ 1 + χ 2 ) 2 w ( 1 s 2 ) J and the nonlinear interaction en G P = 8 ( χ 1 + χ 2 ) 2 w ( 1 s 2 ) and s = v / v 0 .
In this case, the energy of the soliton in Equation (10), or the energy transported by the soliton, is obtained by
E = < Φ ( t ) | H | Φ ( t ) > = 1 r 0 2 [ J r 0 2 ( α x ) 2 + R | α ( x , t ) | 2 G p | α ( x , t ) | 4 d x + 1 r 0 1 2 [ M ( β ( x , t ) t ) 2 + w r 0 ( β ( x , t ) x ) 2 ] d x = E 0 + 1 2 M s o l v 2 .
The static (rest) energy of the soliton is
E 0 = 2 ( ε 0 2 J ) 8 ( χ 1 + χ 2 ) 4 3 w 3 J = E s 0 + W
where W = [ 2 ( χ 1 + χ 2 ) 4 ] / 3 w 2 J is the energy of deformation of the amino acid residues. The effective mass of the soliton is
M s o l = 2 m e x + 8 ( χ 1 + χ 2 ) 4 ( 9 s 2 + 2 3 s 4 ) 3 w 2 J ( 1 s 2 ) 3 v 0 2
In such a case, the binding or forming energy of the soliton in Pang’s theory is denoted by
E B P = 8 ( χ 1 + χ 2 ) 4 3 J w 2
These mathematical models clearly indicate that the properties of Pang’s soliton [75,76,77,78,79,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,100,101] 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−10 s at the physiological temperature of 300 K, which is approximately 800× that of Davydov’s soliton (10−12−10−13 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).

2.3. 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 s o l ), energy (E), and rest energy (E0), 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 (mex), 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 p , then the energy of 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.

3. Variation of Properties of Energy Transport Arising from the EF in Protein Molecules

3.1. 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 (Msol), energy (E), rest energy (E0), and binding energy EBP 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) because the form and outline of the solitary wave, S e c h [ ( μ P / r 0 ) ( x x 0 v t ) ] × exp { i [ v 2 J r 0 2 ( χ 1 + χ 2 ) E v t ] } , in Equation (10) and its amplitude of envelope, μ P / r 0 , as well as its phase of the carrier wave, [ v 2 J r 0 2 ( x x 0 ) E v t ] , are all changed with the variation of dipole–dipole interaction J under the influenceof EF. These changes arising from EFs will also affect the proliferation of cells and life bodies because of variations in the bio-energy they obtain.
(3)
The biological effects of EFs closely depend on both strength and direction with respect to the dipole moment of amino acid residues because the externally applied electric-field, E and dipole moments of the amino acid residue, p , all possess a certain strength and direction. In this case, we should consider the direction and strength of the EF, E , and the dipole moments of the amino acid residue as well as their relationships. Thus, the electromagnetic energy of interaction between them should be represented by E . p = | E | | p | cos θ , where θ is the angle between the two vectors, E and p . This implies that the variations of the dipole–dipole interaction between neighboring amino acid residues caused by the EF should be expressed by ( J + E . p ) J = | E | | p | cos θ , which decreases when the angle (θ) increases. If θ = 0, then E . p = | E | | p | . If θ = 90 0 , then E . p = 0 . 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.

3.2. Results of Numerical Simulation for Changes of Property of the Energy Transport Resulting from EF

3.2.1. 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 ( J + E . p ) 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.
First, we use first the transformation a n ( t ) a n ( t ) exp [ i ε 0 t / ] to eliminate the term ε 0 a n ( t ) in Equation (6). Then, we make the transformation according to the Runge–Kutta method [107,108]: a n ( t ) = a n ( t ) r n + i a ( t ) i n . Thus, Equations (6) and (7) become
α . r n = J ( α i n + 1 + α i n 1 ) + χ 1 ( q n + 1 q n 1 ) α i n + χ 2 ( q n + 1 q n ) ( a i n + 1 + α i n 1 )
α . i n = J ( α r n + 1 + α r n 1 ) + χ 1 ( q n + 1 q n 1 ) α i r n + χ 2 ( q n + 1 q n ) ( a r n + 1 + α r n 1 )
q . n = y / M
y . = W ( q n + 1 2 q n + q n 1 ) + 2 χ 1 ( α r n + 1 2 + α i n + 1 2 α r n 1 2 α i n 1 2 ) + 4 χ 2 [ α r n ( α r n + 1 α r n 1 ) + α i n ( α i n + 1 α i n 1 ) ]
| α n | 2 = | α r n | 2 + | α i n | 2
where arn and ain are real and imaginary parts of an, 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.
In concrete calculations, we must determine the values of the parameters in Equations (15)–(18). The values for the parameters M, J, W, ε 0 , x1, and x2 are known in Pang’s theory of the energy transport in protein molecules. The widely accepted values of these parameters are J = 9.68 × 10−4 eV, ε 0 = 0.205 eV, χ 1 = 6.2 × 10 11 N , w = 13  N / m , χ 2 = ( 10 18 ) × 10 12  N and M = 1.17 × 10−25 kg in single-protein molecules or 5.73 × 10−25 kg in α-helix proteins [24,25,26,27,28,29,30,31,32,33,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88,89,90]. However, we need to confirm the value of the electric dipole moment | p | of the amino acid residues, which is a key parameter for finding the relationship to varying solitons transporting the bio-energy with changing EFs applied externally to protein molecules. Many researchers have calculated and measured the values of | p | of the amino acid residues in protein molecules [24,109,110,111]. Results have shown the electric dipole moments in the region of a few −10 debyes (1 debye = 3.33 × 10−30 cm). Therefore, it is reasonable to choose | p | = 10 debyes in the calculation [24,109,110,111].
From Equations (15)–(19) using the given values of parameters, we can determine the changes of the solution of Equations (15)–(19) numerically, a ( t ) r n and a ( t ) i n , and of the corresponding | a n ( t ) | 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.5 ( M ¯ / w ¯ ) 1 / 2 and the initial excitation of a n ( 0 ) = A sec h [ ( n n 0 ) × ( χ ¯ 1 + χ ¯ 2 ) / 4 J ¯ w ¯ ] (where A is normalization constant) and the initial states of amino acid residues, qn(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, A 0 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 ( t ) | 2 = n | a n ( 0 ) | 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 |an|2=|aj|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−25 kg; J ¯ = 1.55 × 10−22 J; or = 9.68 × 10−4 eV; ε 0 ¯ = 0.20 5 eV; χ 1 ¯ = 6.2 × 10−11 N; w ¯ = (13 − 19.5) N/m; and χ 2 ¯ = ( 11 18 ) × 10 12  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 Figure 2 and Figure 3, respectively, in which we show the movements or progress of the soliton transporting the energy (in other words, its probability | a n ) t ) | 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. Figure 2 and Figure 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, Figure 2 and Figure 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 J + E . p = J + | E | | p | 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 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.

3.2.2. Results in α-Helix Protein Molecules with Three Channels

We also numerically simulated the dynamic states of the soliton that occurred in the α-helix protein molecules with three channels under the influence of EFs. The Hamiltonian and the wave function of the α-helix protein molecules in Pang’s model [86,87,88,89,90,91,92,93,94,95,96,97,98,99,100,101] are, respectively, represented by
H = n α [ ε 0 B n α + B n α J ( B n α + B n + 1 α + B n α B n + 1 α + ) ] + n α [ P n α 2 2 M + 1 2 w ( q n α q n 1 α ) 2 ] + [ χ 1 ( q n + 1 α q n 1 α ) × B n α + B n α + χ 2 ( q n + 1 α q n 1 α ) ( B n + 1 α + B n α + B n α + B n + 1 α ) + L ( B n α + B n + 1 α + B n α + B n 1 α ) ] , α = 1 , 2 , 3.
and
| Φ ( t ) = | α ( t ) | β ( t ) = 1 λ [ 1 + n α α n α ( t ) B n α + + 1 2 ! ( n α α n α ( t ) B n α + ) 2 ] | 0 e x × exp { i n [ q n α ( t ) P n α π n α ( t ) u n α | 0 p h
where L is the coefficient of the chain–chain interactions among the three channels, and α = 1, 2, 3 denote the three chains. According to this above method, from Equations (20) and (21) we can get the dynamic equations of the soliton in the discrete case, which are
i α n α ( t ) . = ε 0 α n α ( t ) J [ ( α n + 1 α ( t ) + + α n 1 α ( t ) ) ] + χ 1 [ ( q n + 1 α ( t ) q n 1 α ( t ) ] α n α ( t ) + χ 2 [ q n + 1 α ( t ) q n 1 α ( t ) ] × [ ( α n + 1 α ( t ) + α n 1 α ( t ) ] + 5 2 { W ( t ) 1 2 m [ q n α ( t ) π m α ( t ) π . n α ( t ) q . m α ( t ) } α n α ( t ) + L [ α n α + 1 ( t ) + α n α 1 ( t ) ]
M q .. n α ( t ) . = w [ q n + 1 α ( t ) 2 q n α ( t ) + q n 1 α ( t ) ] + 2 χ 1 [ | α n + 1 α ( t ) | 2 | α n 1 α ( t ) | 2 ] + 2 χ 2 { α n α * ( t ) [ α n + 1 α ( t ) α n 1 α ( t ) ] + α n α ( t ) [ α n + 1 α * α n 1 α * ( t ) ] }
By means of Equations (22) and (23) and the fourth-order Runge–Kutta method [107,108] mentioned previously, we can also numerically simulate the features of the solutions of Equations (22) and (23), which are denoted in Appendix B.
In this numerical simulation, we use the fourth-order Runge–Kutta method [107,108] with the initial condition of a n α ( 0 ) = A ' sec h [ ( n α n 0 α ) × ( χ ¯ 1 + χ ¯ 2 ) / 4 J ¯ w ¯ ] , where A’ is the normalization factor. Meanwhile, we use M ¯ = 5.73 × 10−25 kg, w ¯ = 39 N/m, ε 0 ¯ = 0.2055 eV, J ¯ = 9.68 × 10−4 eV, χ 1 ¯ = 6.2 × 10−11 N, χ 2 ¯ = (10 − 18) × 10−12 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,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,100,101]. 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) under the influences of EFs.
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.
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.

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

4.1. 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)–CH2–CH2–CH2–CH=N–CH2–CH2–CH2–CH2–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.
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.
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 = I0 exp {−[0.437 + 8.987 × 10−6 (T/°C)2]}
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 Figure 8, Figure 9 and Figure 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.

4.2. 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 Figure 11, Figure 12, Figure 13 and Figure 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 (Figure 12, Figure 13 and Figure 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. Figure 13 and Figure 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, Figure 12, Figure 13 and Figure 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 (Figure 12, Figure 13 and Figure 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.

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

Acknowledgments

The authors would like to acknowledge the National Basic Research Program of China under Grant No. 2011CB503701 for its financial support.

Author Contributions

Xiaofeng Pang, Shude Chen, Xianghui Wang and Lisheng Zhong finished together this researched work. This paper is mainly written by Xiaofeng Pang.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Solutions to Equations (17)–(20)

Note that these are base equations derived from numerical computer simulations and the Runge–Kutta method.
α r j , n + 1 = α r j n + h 6 K 1 j + 2 K 2 j ) + 3 K 3 j + K 4 j ) , α i j , n + 1 = α i j n + h 6 ( L 1 j + 2 L 2 j + 2 L 3 j + L 4 j ) , q j , n + 1 = q j n + h 6 ( M 1 + 2 M 2 j + 2 M 3 j + M 4 j ) , q n + 1 q n 1 ) α i n + χ 2 ( q n + 1 q n ) ( a i n + 1 + α i n 1 ) , y j , n + 1 = y j n + h 6 ( N 1 j + 2 N 3 j + N 4 j ) ,
where
K 1 = J ( α i j + 1 + α i j 1 n + χ 1 ( q j + 1 n q j 1 n ) α i + χ 2 ( q j + 1 n q j n ) ( α i j + 1 n α i j 1 n ) , K 2 j = K 1 j J h 2 ( L 1 j + 2 L 2 j + 2 L 3 j + L 4 j ) + χ 1 h 2 ( M 1 j + 1 M 1 j 1 ) ( a i j n + h 2 L 1 j ) + χ 2 h 2 ( q j + 1 , n q j n ) ( L 1 j + 1 + L 1 j 1 ) , M 1 j = y j n / M ; M 2 j = M 1 j + h 2 M N 1 j ; M 3 j = M 1 j + h 2 M N 2 j ; M 4 j = M 1 j + h M N 3 j ; N 1 j = W ( q j + 1 , n 2 q j n + q j 1 , n ) + 2 χ 1 ( a r j + 1 , n 2 + a i j + 1 , n 2 a r j 1 , n 2 a i j 1 , n 2 ) + 4 χ 2 [ a r j n ( a r j + 1 , n a r j 1 , n ) + a i j n ( a i j + 1 , n a i j 1 , n ) ] ; N 2 j = W [ q j + 1 , n 2 q j n + q j 1 , n + h 2 ( M 1 j + 1 2 M 1 j + M 1 j 1 ) ] + 2 χ 1 [ ( a r j + 1 , n + h 2 K 1 j + 1 ) 2 + ( a i j + 1 , n + h 2 L 1 j + 1 ) 2 ( a r j 1 , n + h 2 K 1 j 1 ) 2 ( a i j 1 , n + h 2 L 1 j 1 ) 2 ] + 4 χ 2 { ( a r j n + h 2 K 1 j ) [ a r j + 1 , n a r j 1 , n + h 2 ( K 1 j + 1 K 1 j 1 ) ] + ( a i j n + h 2 L 1 j ) [ a i j + 1 , n a i j 1 , n + h 2 ( L 1 j + 1 L 1 j 1 ) ] } ; N 3 j = W [ q j + 1 , n 2 q j n + q j 1 , n + h 2 ( M 2 j + 1 2 M 2 j + M 2 j 1 ) ] + 2 χ 1 [ ( a r j + 1 , n + h 2 K 2 j + 1 ) 2 + ( a i j + 1 , n + h 2 L 2 j + 1 ) 2 ( a r j 1 , n + h 2 K 2 j 1 ) 2 ( a i j 1 , n + h 2 L 2 j 1 ) 2 ] + 4 χ 2 { ( a r j n + h 2 K 2 j ) [ a r j + 1 , n a r 1 , n + h 2 ( K 2 j + 1 K 2 j 1 ) ] + ( a i j n + h 2 L 2 j ) [ a i j + 1 , n a i j 1 , n + h 2 ( L 2 j + 1 L 2 j 1 ) ] } ; N 4 j = W [ q j + 1 , n 2 q j n + q j 1 , n + h ( M 3 j + 1 2 M 3 j + M 3 j 1 ) ] + 2 χ 1 [ ( a r j + 1 , n + h K 3 j + 1 ) 2 + ( a i j + 1 , n + h L 3 j + 1 ) 2 ( a r j 1 , n + h K 3 j 1 ) 2 ( a i j 1 , n h L 3 j 1 ) 2 ] + 4 χ 2 { ( a r j n + h K 3 j ) [ a r j + 1 , n a r j 1 , n + h ( K 3 j + 1 K 3 j 1 ) ] + ( a i j n + h L 3 j ) [ a i j + 1 , n a i j 1 , n + h ( L 3 j + 1 L 3 j 1 ) ] }

Appendix B. Solutions to Equations (17)–(20)

From the time-dependent Schrödinger equation
H | Φ > = i t | Φ >
with the Hamiltonian Equation (3) and the time-dependent wave function Equation (4) and
i   t < |   Φ ( t ) | u n | Φ ( t ) > = < Φ ( t ) | [ u n , H ]   | Φ ( t ) >
i   t < Φ ( t ) | P n | Φ ( t ) = < Φ ( t ) |   [ P n , H ]   |   Φ ( t ) >
and considering further the neighboring interactions among the three channels we can determine
i a ˙ n α ( t ) = ε 0 a n α ( t ) J [ a n + 1 α ( t ) + a n 1 α ( t ) ] + χ 1 [ q n + 1 α ( t ) q n 1 α ( t ) ] a n α ( t )   + χ 2 [ q n + 1 α ( t ) q n 1 α ( t ) ] [ a n + 1 α ( t ) + a n 1 α ( t ) ]   + 5 2 { w ( t ) 1 2 m [ q m α ( t ) π m α ( t ) π ˙ m α ( t ) q ˙ m α ( t ) ] } a n α ( t ) + L [ a n α + 1 ( t ) + a n α 1 ( t ) ]
M q ¨ n α = W [ q n + 1 α ( t ) 2 q n α ( t ) + q n 1 α ( t ) ] + 2 χ 1 [ | a n + 1 α ( t ) | 2 | a n 1 α ( t ) | 2 ]   + 2 χ 2 { a n α * ( t ) [ a n + 1 α ( t ) a n 1 α ( t ) ] + a n α ( t ) [ a * n + 1 α ( t ) a * n 1 α ( t ) ] } .
Also, we can eliminate the term containing ε0 in Equation (B4) through the following transformation:
φ n α ( t ) = a n α ( t ) exp [ i ε 0 t / ]
Because an(t) in Equations (B4) and (B5) is a complex function, we can make the following transformation
a n α ( t ) = a r n α ( t ) + i a i n α ( t )  with  | a n | 2 = | a r n | 2 + | a i n | 2
Thus Equations (B4) and (B5) change as
a ˙ r n α = J ( a i n + 1 α + a i n 1 α ) + χ 1 ( q n + 1 α q n 1 α ) a i n α + χ 2 ( q n + 1 q n 1 α ) ( a i n + 1 α + a i n 1 α ) + L [ a i n α + 1 ( t ) + a i n α 1 ( t ) ]
a ˙ i n α = J ( a r n + 1 α + a r n 1 α ) + χ 1 ( q n + 1 α q n 1 α ) a r n α + χ 2 ( q n + 1 α q n 1 α ) ( a r n + 1 α + a r n 1 α ) + L [ a i n α + 1 ( t ) + a i n α 1 ( t ) ]
q ˙ n α = y n α M
y ˙ n α = W [ q n + 1 α 2 q n α + q n 1 α ] + 2 χ 1 [ a r n + 1 α 2 + a i n + 1 α 2 a i n 1 α 2 a i n 1 α 2 ]   + 4 χ 2 [ a r n α ( a r n + 1 α a r n 1 α ) + a i n α ( a i n + 1 α a i n 1 α ) ] ,
where arn and ain are real and imaginary parts of an(t), respectively. Therefore, Equations (B8)–(B11) are dynamic in Equations (17)–(20), which can be used to simulate and calculate the dynamic properties of particles. Therefore, the dynamic properties of the peptide groups or amino acid molecules mentioned previously can be calculated and studied through Equations (B8)–(B11), as detailed in Appendix A for single-chain protein molecules.

References

  1. Pang, X.F. Bio-Electromagnetism, 1st ed.; The Defense Industry Press of China: Beijing, China, 2009; pp. 14–52. [Google Scholar]
  2. Fröhlich, H. Biological Coherence and Response to External Stimuli, 1st ed.; Springer-Verlag: New York, NY, USA, 1988; pp. 5–64. [Google Scholar]
  3. Ahlbom, A.; Day, N.; Feychting, M.; Skinner, R.E.; Dockerty, J.; Linet, J.; McBride, M.; Michaelis, M.; Olsen, J.; Tynes, T.; et al. A pooled analysis of magnetic fields and childhood leukaemia. Br. J. Cancer 2000, 83, 692–698. [Google Scholar] [CrossRef] [PubMed]
  4. Dowson, D.I.; Lewith, G.; Campbell, M.; Mullee, M.A.; Brewster, L.A. Overhead high-voltage cables and recurrent headache and depressions. Practitioner 1988, 232, 435–436. [Google Scholar] [PubMed]
  5. Greenland, S.; Sheppard, A.R.; Kaune, W.T.; Poole, C.; Kelsh, M.A. A Pooled analysis of magnetic fields, wire codes, and childhood leukemia. Epidemiology 2000, 11, 624–634. [Google Scholar] [CrossRef] [PubMed]
  6. Perry, S.; Pearl, L.; Binns, R. Power frequency magnetic field: Depressive Illness and myocardial infarction. Public Health 1989, 103, 177–180. [Google Scholar] [CrossRef]
  7. Perry, F.S.; Reichmanis, M.; Marino, A.A.; Becker, R.O. Environmental power-frequency magnetic fields and suicide. Health Phys. 1981, 41, 267–277. [Google Scholar] [CrossRef] [PubMed]
  8. Bromberg, J.; Darnell, J.E. The role of STATs in transcriptional control and their impact on cellular function. Oncogene 2000, 19, 2468–2473. [Google Scholar] [CrossRef] [PubMed]
  9. Pang, X.F. Biophysics; Press of University of Electronic Science and Technology of China: Chengdu, China, 2008; pp. 67–173. [Google Scholar]
  10. Li, G.; Pang, X.F. Biological effects of environmental electromagnetic fields. Adv. Mater. Res. 2011, 183–185, 532–536. [Google Scholar] [CrossRef]
  11. Li, G.; Pang, X.F. The influence of electromagnetic field on the electromagnetic properties of biological tissue. Prog. Biochem. Biophys. (Chin.) 2011, 38, 604–610. [Google Scholar] [CrossRef]
  12. Li, G.; Li, B.; Yang, S.F.; Pang, X.F. The influences of static magnetic field on the electric features of testicular tissue of rats. Aerosp. Med. Med. Eng. 2012, 25, 172–175. [Google Scholar]
  13. Li, G.; Yan, Y.Y.; Han, S.G.; Pang, X.F. The influences of static magnetic field on the brain tissue of the rats using infrared spectrum metho. Laser Infrared 2012, 42, 894–896. [Google Scholar]
  14. Li, G.; Yan, Y.J.; Huan, Y.; Zhou, Y.; Pang, X.F. The influences of electromagnetic field with extremely low frequency on the characteristics of infrared spectrum of the sensitive tissue in rat. Spectrosc. Spectr. Anal. 2012, 32, 1194–1197. [Google Scholar]
  15. Li, G.; Li, B.; Yan, Y.J.; Yang, S.F.; Pang, X.F. Calculation of equivalent electrical parameters in human tissue and its applications. Mater. Rev. 2011, 25, 142–145. [Google Scholar]
  16. Zhang, Y.M.; Zhou, Y.; Pang, X.F. The Investigation of effects of pulse electromagnetic field on the fluorescence spectrum of serum in rat. Spectrosc. Spectr. Anal. 2012, 32, 2162–2165. [Google Scholar]
  17. Zhang, Y.M.; Li, G.; Yan, Y.J.; Pang, X.F. Altered expression of matrix metalloproteinases and tight junction proteins with PEMF-induced BBB permeability change in rats. Biomed. Environ. Sci. 2011, 24, 408–414. [Google Scholar]
  18. Beebe, S.J.; Fox, P.M.; Rec, L.J.; Somers, K.; Stark, R.H.; Schoenbach, K.H. Nanosecond pulsed electric field (nsPEF) effects on cells and tissues: Apoptosis induction and tumor growth inhibition. IEEE Trans. Plasma Sci. 2002, 30, 286–292. [Google Scholar] [CrossRef]
  19. Lazzaro, V.D.; Capone, F.; Apollonio, F.; Borea, P.A.; Cadossi, R.; Fassina, L.; Grassi, C.; Liberti, M.; Paffi, A.; Parazzini, M.; et al. A consensus panel review of central nervous system effects of the exposure to low-intensity extremely low-frequency magnetic fields. Brain Stimul. 2013, 6, 469–476. [Google Scholar] [CrossRef] [PubMed]
  20. Apollonio, F.; Liberti, M.; Paffi, A.; Merla, C.; Marracino, P.; Denzi, A.; Marino, C.; d’Inzeo, G. Feasibility for microwaves energy to affect biological systems via nonthermal mechanisms: A systematic approach. IEEE Trans. Microw. Theory Tech. 2013, 61, 2031–2045. [Google Scholar] [CrossRef]
  21. Hofmann, G.A.; Dev, S.B.; Dimer, S.; Nanda, G.S. Electroporation therapy: A new approach for treatment of head and neck cancer. IEEE Trans. Biomed. Eng. 1999, 46, 752–759. [Google Scholar] [CrossRef] [PubMed]
  22. Kekez, M.M.; Savic, P.; Johnson, B.F. Contribution to the biophysics of the lethal effects of electric field on microorganisms. Biochim. Biophys. Acta 1996, 1278, 79–88. [Google Scholar] [CrossRef]
  23. Pang, X.F.; Zhang, A.Y. Investigations of mechanism and properties of biological thermal effects of the microwave. China J. Atom. Mol. Phys. 2001, 18, 278–281. [Google Scholar]
  24. Davydov, A.S. Biology and Quantum Mechanics; Pergaman: New York, NY, USA, 1982; pp. 161–211. [Google Scholar]
  25. Davydov, A.S. Solitons in Molecular Systems, 2nd ed.; Reidel Publishing Company: Dordrecht, The Netherlands, 1991; pp. 23–79. [Google Scholar]
  26. Davydov, A.S.; Ermakov, V.N. Soliton generation at the boundary of molecular chain. Phys. D 1988, 32, 318–329. [Google Scholar] [CrossRef]
  27. Davydov, A.S. Solitons in molecular systems. Phys. Scr. 1979, 20, 387–394. [Google Scholar] [CrossRef]
  28. Davydov, A.S. The lifetime of molecular solitons. J. Biol. Phys. 1991, 18, 111–125. [Google Scholar] [CrossRef]
  29. Davydov, A.S. Solitons and energy transfer along protein molecules. J. Theor. Biol. 1977, 66, 379–387. [Google Scholar] [CrossRef]
  30. Scott, A.C. Davydov’s soliton. Phys. Rep. 1992, 217, 1–67. [Google Scholar] [CrossRef]
  31. Scott, A.C. Dynamics of Davydov’s soliton. Phys. Rev. A 1982, 26, 578–595. [Google Scholar] [CrossRef]
  32. Brown, D.W.; Ivic, Z. Unification of polaron and soliton theories of exciton transport. Phys. Rev. E 1989, 40, 9876–9887. [Google Scholar] [CrossRef]
  33. Brown, D.W.; Lindenberg, K.; West, B.J. Nonlinear density-matrix equation for the study of finite-temperature soliton dynamics. Phys. Rev. B 1987, 35, 6169–6180. [Google Scholar] [CrossRef]
  34. Guo, B.L.; Pang, X.F. Solitons; Chinese Science Press: Beijing, China, 1987; pp. 4–140. [Google Scholar]
  35. Pang, X.F. The transport of bio-energy in protein molecules. Chin. J. Biochem. Biophys. 1986, 18, 1–8. [Google Scholar]
  36. Christiansen, P.L.; Scott, A.C. Davydov’s Soliton Revisited; Plenum Press: New York, NY, USA, 1990; pp. 34–103. [Google Scholar]
  37. Cruzeiro, L.; Halding, J.; Christiansen, P.L.; Skovgard, O.; Scott, A.C. Temperature effects on Davydov soliton. Phys. Rev. A 1985, 37, 880–887. [Google Scholar] [CrossRef]
  38. Cruzeio-Hansson, L. Mechanism of thermal destabilization of the Davydov soliton. Phys. Rev. A 1992, 45, 4111–4115. [Google Scholar] [CrossRef]
  39. Forner, W. Quantum and disorder effects in Davydov soliton theory. Phys. Rev. A 1991, 44, 2694–2708. [Google Scholar] [CrossRef] [PubMed]
  40. Forner, W. Quantum and temperature effects on Davydov soliton dynamics: Averaged Hamiltonian method. J. Phys. Condens. Matter 1992, 4, 1915–1923. [Google Scholar] [CrossRef]
  41. Forner, W. Davydov soliton dynamics: Two-quantum states and diagonal disorder. J. Phys. Condens. Matter 1991, 3, 3235–3254. [Google Scholar] [CrossRef]
  42. Lomdahl, P.S.; Kerr, W.C. Do Davydov solitons exist at 300 K? Phys. Rev. Lett. 1985, 55, 1235–1239. [Google Scholar] [CrossRef] [PubMed]
  43. Kerr, W.C.; Lomdahl, P.S. Quantum-mechanical derivation of the equations of motion for Davydov solitons. Phys. Rev. B 1987, 35, 3629–3632. [Google Scholar] [CrossRef]
  44. Wang, X.; Brown, D.W.; Lindenberg, K.; West, B. Alternative formulation of Davydov theory of energy transport in biomolecules systems. Phys. Rev. A 1988, 37, 3557–3566. [Google Scholar] [CrossRef]
  45. Wang, X.; Brown, D.W.; Lindenberg, K. Quantum Monte Carlo simulation of Davydov model. Phys. Rev. Lett. 1989, 62, 1796–1799. [Google Scholar] [CrossRef] [PubMed]
  46. Cottingham, J.P.; Schweitzer, J.W. Calculation of the lifetime of a Davydov soliton at finite temperature. Phys. Rev. Lett. 1989, 62, 1792–1795. [Google Scholar] [CrossRef] [PubMed]
  47. Schweitzer, J.W. Lifetime of the Davydov soliton. Phys. Rev. A 1992, 45, 8914–8922. [Google Scholar] [CrossRef] [PubMed]
  48. Mettle, B.; Shaw, P.B. Evolution of a molecular exciton on a Davydov lattice at T = 0. Phys. Rev. B 1988, 38, 3075–3086. [Google Scholar]
  49. Ivic, Z.; Przulj, Z.; Kosti, D. Decay and slowing down of the multiquanta Davydov-like solitons in molecular chains. Phys. Rev. E 2000, 61, 6963–6967. [Google Scholar] [CrossRef]
  50. Teki, J.; Ivic, Z.; Zekovi, S.; Przulj, Z. Kinetic properties of multiquanta Davydov-like solitons inmolecular chains. Phys. Rev. E 1999, 60, 821–825. [Google Scholar] [CrossRef]
  51. Zekovi, S.; Ivic, Z. Damping and modification of the multiquanta Davydov-like solitons in molecular chains. Bioelectrochem. Bioenerg. 1999, 48, 297–300. [Google Scholar] [CrossRef]
  52. Pouthier, V. Two-vibron bound states in alpha-helix proteins: The interplay between the intramolecular anharmonicity and the strong vibron-phonon coupling. Phys. Rev. E 2003, 68. [Google Scholar] [CrossRef] [PubMed]
  53. Pouthier, V.; Falvo, C. Relaxation channels of two-vibron bound states in α-helix proteins. Phys. Rev. E 2004, 69. [Google Scholar] [CrossRef] [PubMed]
  54. Falvo, C.; Pouthier, V. Vibron–polaron in α-helices. I. Single-vibron states. J. Chem. Phys. 2005, 123, 184709–184716. [Google Scholar] [CrossRef] [PubMed]
  55. Falvo, C.; Pouthier, V. Vibron–polaron in α-helices. II. Two-vibron states. J. Chem. Phys. 2005, 123. [Google Scholar] [CrossRef] [PubMed]
  56. Cruzeiro, L. The Davydov/Scott model for energy storage and transport in proteins. J. Biol. Phys. 2009, 35, 43–55. [Google Scholar] [CrossRef] [PubMed]
  57. Silva, P.A.S.; Cruzeiro, L. Dynamics of a nonconserving Davydov monomer. Phys. Rev. E 2006, 74. [Google Scholar] [CrossRef] [PubMed]
  58. Silva, P.A.S.; Cruzeiro-Hansson, L. A reduced set of exact equations of motion for a non-number-conserving Hamiltonian. Phys. Lett. A 2003, 315, 447–451. [Google Scholar] [CrossRef]
  59. Pouthier, V. Energy relaxation of the amide-I mode in hydrogen-bonded peptide units: A route to conformational change. J. Chem. Phys. 2008, 128. [Google Scholar] [CrossRef] [PubMed]
  60. Pouthier, V.; Tsybin, Y.O. Amide-I relaxation-induced hydrogen bond distortion: Anintermediate in electron capture dissociation mass pectrometry of α-helical peptides? J. Chem. Phys. 2008, 129. [Google Scholar] [CrossRef] [PubMed]
  61. Moritsugu, K.; Miyashita, O.; Kidera, K. Vibrational energy transfer in a protein molecule. Phys. Rev. Lett. 2000, 85, 3970–3973. [Google Scholar] [CrossRef] [PubMed]
  62. Fujisaki, H.; Zhang, Y.; Straub, J.E. Time-dependent perturbation theory for vibrational energy relaxation and dephasing in peptides and proteins. J. Chem. Phys. 2006, 124. [Google Scholar] [CrossRef] [PubMed]
  63. Pang, X.F. The properties of collective excitation in organic protein molecular system. J. Phys. Condens. Matter 1990, 2, 9541–9553. [Google Scholar]
  64. Pang, X.F. Comment the thermodynamic properties of α-helix protein: A soliton approach. Phys. Rev. 1994, 49, 4747–4751. [Google Scholar]
  65. Pang, X.F. Nonlinear Quantum Mechanical Theory; Chongqing Press: Chongqing, China, 1994; pp. 356–465. [Google Scholar]
  66. Pang, X.F. Properties of soliton in protein molecules with nonlinear nearest neighbour interaction. Chin. Sci. Bull. 1993, 38, 1572–1581. [Google Scholar]
  67. Pang, X.F. The thermodynamic properties of the solitons excited in the protein molecules. Chin. Sci. Bull. 1993, 38, 1665–1677. [Google Scholar]
  68. Pang, X.F. Quantum-mechamical method for the soliton transported bio-energy in protein. Chin. Phys. Lett. 1993, 10, 437–440. [Google Scholar] [CrossRef]
  69. Pang, X.F. Stability of the soliton excited in protein in the biological temperature range. Chin. Phys. Lett. 1993, 10, 573–576. [Google Scholar] [CrossRef]
  70. Pang, X.F. Influence of the soliton in anharmonic molecular crystals with temperature on Mossbauer effect. Eur. Phys. J. B 1999, 10, 415–423. [Google Scholar] [CrossRef]
  71. Pang, X.F.; Chen, X.R. Distribution of vibrational energy- levels of protein molecular chains. Commun. Theor. Phys. 2001, 35, 323–326. [Google Scholar]
  72. Pang, X.F. The effect of Raman scattering accompanied by the soliton excitation occurred in the molecular crystals. Phys. D 2001, 154, 138–158. [Google Scholar] [CrossRef]
  73. Pang, X.F. An improvement of the Davydov theory of bio-energy etransport in the protein molecular systems. Phys. Rev. E 2000, 62, 6989–6998. [Google Scholar]
  74. Pang, X.F. The lifetime of the soliton in the improved Davydov model at the biological temperature 300 K for protein molecules. Eur. Phys. J. B 2001, 19, 297–316. [Google Scholar]
  75. Pang, X.F.; Zhang, H.W.; Yu, J.F.; Feng, Y.P. States and properties of the soliton transported bio-energy in nonuniform protein molecules at physiological temperature. Phys. Lett. A 2005, 335, 408–416. [Google Scholar] [CrossRef]
  76. Pang, X.F.; Luo, Y.H. Stabilization of the soliton transported bio-energy in protein molecules in the improved model. Commun. Theor. Phys. 2004, 41, 470–476. [Google Scholar]
  77. Pang, X.F.; Feng, Y.P. Quantum Mechanics in Nonlinear Systems; World Science Publisher Co.: Singapore, 2005; pp. 487–556. [Google Scholar]
  78. Pang, X.F.; Liu, M.J. Features of Motion of Soliton Transported Bio-energy in aperiodic alpha-helix protein molecules with three channels. Commun. Theor. Phys. 2009, 51, 170–180. [Google Scholar]
  79. Pang, X.F. The effects of damping and temperature of medium on the soliton excited in α-helix protein molecules with three channels. Mod. Phys. Lett. B 2009, 23, 71–88. [Google Scholar] [CrossRef]
  80. Pang, X.F.; Lui, M.J. The influences of temperature and chain-chain interaction on features of solitons excited in α-Helix protein molecules with three channels. Int. J. Mod. Phys. B 2009, 23, 2303–2322. [Google Scholar] [CrossRef]
  81. Pang, X.F. Influence of structure disorders and temperatures of systems on the bio-energy transport in protein molecules. Front. Phys. China 2008, 3, 457–488. [Google Scholar] [CrossRef]
  82. Pang, X.F.; Yu, J.F.; Lao, Y.H. Combination effects of structure nonuniformity of proteins on the soliton transported bio-energy. Int. J. Mod. Phys. B 2007, 21, 13–42. [Google Scholar] [CrossRef]
  83. Pang, X.F.; Liu, M.J. Properties of soliton-transported bio-energy in α-helix protein molecules with three channels. Commun. Theor. Phys. 2007, 48, 369–376. [Google Scholar]
  84. Pang, X.F. Theory of bio-energy transport in protein molecules and its experimental evidences as well as applications (I). Front. Phys. China 2007, 2, 469–493. [Google Scholar] [CrossRef]
  85. Pang, X.F.; Zhang, H.W. Theoretical investigation of properties of infrared absorption of α-helix protein molecules. Int. J. Infrared Millim. Waves 2006, 27, 735–812. [Google Scholar] [CrossRef]
  86. Pang, X.F.; Zhang, H.W.; Liu, M.J. Influences of heat bath and structure disorder in protein molecules on the soliton transported bio-energy in the improved model. J. Phys. Condens. Matter 2006, 18, 613–627. [Google Scholar] [CrossRef]
  87. Pang, X.F.; Zhang, H.W.; Yu, J.F.; Luo, Y.H. Influences of variations of characteristic parameters arising from the structure nonuniformity of the protein molecules on states of the soliton transported bio-energy in the improved model. Int. J. Mod. Phys. B 2006, 20, 3027–3041. [Google Scholar] [CrossRef]
  88. Pang, X.F.; Chen, X.R. The properties of nonlinear energy-spectra of acetanilide. Int. J. Mod. Phys. 2006, 20, 2505–2518. [Google Scholar] [CrossRef]
  89. Pang, X.F.; Yu, J.F.; Luo, Y.H. Influences of quantum and disorder effects on solitons exited in protein molecules in improved model. Commun. Theor. Phys. 2005, 43, 367–376. [Google Scholar]
  90. Pang, X.F.; Zhang, H.W. Changes of the Mössbauer effect caused by the excitation of the solitons in the organic molecular crystals at finite temperature. J. Phys. Chem. Solids 2005, 66, 963–972. [Google Scholar]
  91. Pang, X.F.; Zhang, H.W.; Yu, J.F.; Luo, Y.H. Thermal stability of the new soliton transported bio-energy under influence of fluctuations of characteristic Parameters at biological temperature in the protein molecules. Int. J. Mod. Phys. B 2005, 19, 4677–4699. [Google Scholar] [CrossRef]
  92. Pang, X.F.; Yu, J.F.; Liu, M.J. Changes of properties of the soliton with temperature under influences of structure disorder in the α-helix protein molecules with three channels. Mol. Phys. 2010, 108, 1297–1315. [Google Scholar] [CrossRef]
  93. Pang, X.F.; Chen, X.R. Calculation of vibrational energy-spectra of α-helical protein molecules and its properties. Commun. Theor. Phys. 2002, 37, 715–722. [Google Scholar]
  94. Pang, X.F.; Chen, X.R. Vibrational energy-spectra and infrared absorption of α-helical protein molecules. Chin. Phys. Lett. 2002, 19, 1096–1099. [Google Scholar]
  95. Pang, X.F. Nonlinear Quantum Mechanics; Chinese Electronic Industry Press: Beijing, China, 2009; pp. 356–412. [Google Scholar]
  96. Pang, X.F.; Xiao, H.L.; Cue, G.P.; Zhang, H.W.; Dong, B. Experiment studies of properties of infrared absorption of biological tissues. Int. J. Infrared Millim. Waves 2010, 31, 521–532. [Google Scholar]
  97. Pang, X.F. The properties of bio-energy transport and influence of structure nonuniformity and temperature of systems on energy transport along polypeptide chains. Prog. Biophys. Mol. Biol. 2012, 108, 1–46. [Google Scholar] [PubMed]
  98. Pang, X.F. The mechanism and properties of bio-photon emission and absorption in protein molecules in living systems. J. Appl. Phys. 2012, 111. [Google Scholar] [CrossRef]
  99. Pang, X.F. The theory of bio-energy transport in the protein molecules and its properties. Phys. Life Rev. 2011, 8, 264–286. [Google Scholar] [CrossRef] [PubMed]
  100. Pang, X.F. Correctness and completeness of the theory of bio-energy transport. Phys. Life Rev. 2011, 8, 302–306. [Google Scholar] [CrossRef]
  101. Pang, X.F. The checkout and verification of theory of bio-energy transport in the protein molecules. Biophys. Rev. Lett. 2014, 9, 1–79. [Google Scholar] [CrossRef]
  102. Pang, X.F. The features of infrared spectrum of bio-polymer and its theoretical investigation. Int. J. Mod. Phys. B 2014, 29. [Google Scholar] [CrossRef]
  103. Su, X.D.; Jin, F.J. Davydov-Pang model: An improved Davydov protein soliton theory. Phys. Life Rev. 2011, 8, 300–301. [Google Scholar] [CrossRef] [PubMed]
  104. Liang, S.D. Physical insights to the bio-energy transport in the protein molecules. Phys. Life Rev. 2011, 8, 287–288. [Google Scholar] [CrossRef] [PubMed]
  105. Tao, S. The function of soliton on bio-energy transport in the protein molecules. Phys. Life Rev. 2011, 8, 291–292. [Google Scholar]
  106. He, N.Y. An important biological theory—Solving the transport of bio-energy in living systems. Phys. Life Rev. 2011, 8, 296–297. [Google Scholar] [CrossRef] [PubMed]
  107. Stiefel, J. Einfuhrung in die Numerische Mathematic; Teubner Verlag: Stuttgart, Germany, 1965; pp. 78–156. [Google Scholar]
  108. Atkinson, K.E. An Introdution to Numerical Analysis; Wiley New York Inc.: New York, NY, USA, 1987; pp. 34–104. [Google Scholar]
  109. Kumar, A.; Singh, G. Measurement of dielectric constant and loss factor of the dielectric material at microwave frequencies. Prog. Electromagn. Res. 2007, 69, 47–54. [Google Scholar] [CrossRef]
  110. Clifford, E.F.; Prilusky, J.; Silman, I.; Sussman, J.L. A server and database for dipole moments of proteins. Nucleic Acids Res. Adv. 2007, 35, W512–W521. [Google Scholar]
  111. Antosiewicz, J. Computation of the dipole moments of proteins. Biophys. J. 1995, 69, 1344–1354. [Google Scholar] [CrossRef]
  112. Marracino, P.; Amadei, A.; Apollonio, F.; d’Inzeo, G.; Liberti, M.; diCrescenzo, A.; Fontana, A.; Zappacosta, R.; Aschi, M. Modeling of chemical reactions in micelle: Water-mediated, Keto-Enol Interconversion as a case study. J. Phys. Chem. 2011, 115, 8102–8111. [Google Scholar] [CrossRef] [PubMed]
  113. Careri, G.; Buontempo, U.; Galluzzi, F.; Scott, A.C.; Gratton, E.; Shyamsunder, E. Spectroscopic evidence for Davydov-like solitons in acetanilide. Phys. Rev. B 1984, 30, 4689–4698. [Google Scholar] [CrossRef]
  114. Careri, G.A.; Gransanti, A.; Rupley, J.A. Critical exponents of protonic percolation in hydrated lysozyme powders. Phys. Rev. A 1988, 37, 2703–2705. [Google Scholar] [CrossRef]
  115. Careri, G.A.; Gratton, E.A.; Shyamsunder, E. Fine structure of the amide-I bond in acetanilide. Phys. Rev. A 1988, 37, 4048–4051. [Google Scholar] [CrossRef]
  116. Careri, G.A.; Buonttempo, U.; Caeta, F.; Scott, A.C. Infrared absorption in acetanilide by solitons. Phys. Rev. Lett. 1983, 51, 304–307. [Google Scholar] [CrossRef]
  117. Careri, G.A.; Eilbeck, J.C. Stability of solutions of the discrete self-trapping equation. Phys. Lett. A 1985, 109, 201–204. [Google Scholar]
  118. Scott, A.C. Davydov’s soliton revisited. Phys. D 1990, 51, 333–342. [Google Scholar] [CrossRef]
  119. Scott, A.C.; Gratton, E.; Shyamsunder, E.; Careri, G. IR overtone spectrum of the vibrational soliton crystalline acetanilide. Phys. Rev. B 1985, 32, 5551–5553. [Google Scholar] [CrossRef]
  120. Scott, A.C.; Bigio, I.J.; Johnston, C.T. Polarons in acetanide. Phys. Rev. B 1989, 39, 12883–12887. [Google Scholar] [CrossRef]
  121. Alexander, D.M.; Krumbansl, J.A. Localized excitation in hydeogen-bonded molecular crystals. Phys. Rev. B 1986, 33, 7172–7185. [Google Scholar] [CrossRef]
  122. Alexander, D.M. Analog of small Holstein polaron in hydrogen-bonded amide systems. Phys. Rev. Lett. 1985, 60, 138–141. [Google Scholar] [CrossRef] [PubMed]
  123. Weaver, R.F. Molecular Biology; McGraw Hill: Boston, MA, USA, 2002; pp. 41–125. [Google Scholar]
  124. Chen, H.L. Structure and Function of Biomacromolecules; Press of Shanghu Medicine University: Shanghu, China, 1998; pp. 113–122. [Google Scholar]
  125. Pang, X.F. Molecular Bio-Physics; Printer of University of Electronic Science and Technology of China: Chengdu, China, 2003; pp. 68–71. [Google Scholar]
  126. Susi, H.; Ard, T.S.; Carroll, J. The infrared spectrum and water binding of collagen as a function of relative humidity. Biopolymers 1971, 10, 1597–1604. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The molecular structure of an α-helical protein molecule.
Figure 1. The molecular structure of an α-helical protein molecule.
Ijms 17 01130 g001
Figure 2. State of the new soliton in the case of a longtime period of 250 ps and long spacing of 400 amino acid residues.
Figure 2. State of the new soliton in the case of a longtime period of 250 ps and long spacing of 400 amino acid residues.
Ijms 17 01130 g002
Figure 3. The collision behavior of two solitons for Equations (18)–(21).
Figure 3. The collision behavior of two solitons for Equations (18)–(21).
Ijms 17 01130 g003
Figure 4. States of Pang’s soliton in the case of (a) | E | = 25,500; (b) 51,000; (c) 76,500; and (d) 102,000 V/m, respectively, where “site k” denoted the number of amino acid residue, ”n”.
Figure 4. States of Pang’s soliton in the case of (a) | E | = 25,500; (b) 51,000; (c) 76,500; and (d) 102,000 V/m, respectively, where “site k” denoted the number of amino acid residue, ”n”.
Ijms 17 01130 g004aIjms 17 01130 g004b
Figure 5. States of (a) movement and collision features of (b) Pang’s solitons in uniform α-helix proteins.
Figure 5. States of (a) movement and collision features of (b) Pang’s solitons in uniform α-helix proteins.
Ijms 17 01130 g005
Figure 6. States of solitons in α-helix protein molecules of which | E | = (a) 17,000; (b) 25,500; and (c) 34,000 V/m, respectively.
Figure 6. States of solitons in α-helix protein molecules of which | E | = (a) 17,000; (b) 25,500; and (c) 34,000 V/m, respectively.
Ijms 17 01130 g006aIjms 17 01130 g006b
Figure 7. Molecular structure of collagen: (a) fundamental structure of right α-helix with three channels; (b) atomic distribution on top section of three channel–axes, where G is carbon atom, is hydrogen bond; (c) one-dimensional structure of gly-pro-hydropro.
Figure 7. Molecular structure of collagen: (a) fundamental structure of right α-helix with three channels; (b) atomic distribution on top section of three channel–axes, where G is carbon atom, is hydrogen bond; (c) one-dimensional structure of gly-pro-hydropro.
Ijms 17 01130 g007
Figure 8. The infrared spectrum of collagen in 1540–1740 cm−1.
Figure 8. The infrared spectrum of collagen in 1540–1740 cm−1.
Ijms 17 01130 g008
Figure 9. Changes of intensity for infrared absorption of collagen in the region at (1) 95 °C; (2) 85 °C; (3) 75; (4) 65 °C; (5) 55 °C; (6) 45 °C; (7) 35 °C; (8) 25 °C; and (9) 15 °C.
Figure 9. Changes of intensity for infrared absorption of collagen in the region at (1) 95 °C; (2) 85 °C; (3) 75; (4) 65 °C; (5) 55 °C; (6) 45 °C; (7) 35 °C; (8) 25 °C; and (9) 15 °C.
Ijms 17 01130 g009
Figure 10. Temperature dependences of intensity of peaks in the infrared spectrum of collagen: (a) Relation of logarithm of relative intensity, Ln (I/I0), versus (T/°C)2 for a 1650-cm−1 peak, where “●” denotes experimental data; and (b) linear temperature dependence of relative intensity, I/I0, for a 1666-cm−1 peak.
Figure 10. Temperature dependences of intensity of peaks in the infrared spectrum of collagen: (a) Relation of logarithm of relative intensity, Ln (I/I0), versus (T/°C)2 for a 1650-cm−1 peak, where “●” denotes experimental data; and (b) linear temperature dependence of relative intensity, I/I0, for a 1666-cm−1 peak.
Ijms 17 01130 g010
Figure 11. Infrared spectrum of absorption of collagen in 480–2000 cm−1.
Figure 11. Infrared spectrum of absorption of collagen in 480–2000 cm−1.
Ijms 17 01130 g011
Figure 12. Changes in the relative height of a 1650-cm−1 peak with increases to externally applied voltage.
Figure 12. Changes in the relative height of a 1650-cm−1 peak with increases to externally applied voltage.
Ijms 17 01130 g012
Figure 13. Changes in the position of a 1650-cm−1 peak with increases to the externally applied voltage.
Figure 13. Changes in the position of a 1650-cm−1 peak with increases to the externally applied voltage.
Ijms 17 01130 g013
Figure 14. Changes to the width of a half-peak for a 1650-cm−1 peak with increases to the externally applied voltage.
Figure 14. Changes to the width of a half-peak for a 1650-cm−1 peak with increases to the externally applied voltage.
Ijms 17 01130 g014
Table 1. Comparison of soliton features in Pang’s model and the Davydov model.
Table 1. Comparison of soliton features in Pang’s model and the Davydov model.
Lifetime at 300 K (S)Critical Temperature (K)Number of Amino Acids Traveled by Soliton In LifetimeNonlinear Interaction G (×10−21 J)Amplitude of SolitonWidth of the Soliton (×10−10 m)Binding Energy of Soliton (×10−21 J)
Pang’s model10−9–10−10320Several hundreds1.180.97414.88−0.188
Davydov model10−12–10−13<200<103.81.724.95−7.8

Share and Cite

MDPI and ACS Style

Pang, X.; Chen, S.; Wang, X.; Zhong, L. Influences of Electromagnetic Energy on Bio-Energy Transport through Protein Molecules in Living Systems and Its Experimental Evidence. Int. J. Mol. Sci. 2016, 17, 1130. https://doi.org/10.3390/ijms17081130

AMA Style

Pang X, Chen S, Wang X, Zhong L. Influences of Electromagnetic Energy on Bio-Energy Transport through Protein Molecules in Living Systems and Its Experimental Evidence. International Journal of Molecular Sciences. 2016; 17(8):1130. https://doi.org/10.3390/ijms17081130

Chicago/Turabian Style

Pang, Xiaofeng, Shude Chen, Xianghui Wang, and Lisheng Zhong. 2016. "Influences of Electromagnetic Energy on Bio-Energy Transport through Protein Molecules in Living Systems and Its Experimental Evidence" International Journal of Molecular Sciences 17, no. 8: 1130. https://doi.org/10.3390/ijms17081130

APA Style

Pang, X., Chen, S., Wang, X., & Zhong, L. (2016). Influences of Electromagnetic Energy on Bio-Energy Transport through Protein Molecules in Living Systems and Its Experimental Evidence. International Journal of Molecular Sciences, 17(8), 1130. https://doi.org/10.3390/ijms17081130

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop