Abstract
An αhelical region of protein molecule was considered in a model of nanotube. The molecule is in conditions of quantum excitations. Such model corresponds to a onedimensional molecular nanocrystal with three molecules in an elementary cell at the presence of excitation. For the analysis of different types of conformational response of the αhelical area of the protein molecule on excitation, the nonlinear response of this area to the intramolecular quantum excitation caused by hydrolysis of adenosine triphosphate (ATP) is taken into account. It has been established that in the simplest case, three types of excitation are realized. As estimates show, each of them ‘serves’ different kinds of protein. The symmetrical type of excitation, most likely, is realized in the reduction of traversalstriped skeletal muscles. It has the highest excitation energy. This well protects from casual actions. Antisymmetric excitations have intermediate energy (between symmetrical and asymmetrical). They, most likely, are realized in membranous and nucleic proteins. It is shown that the conformational response of the αhelical region of the protein is (in angstroms) a quantity of order N_{c}/5, where N_{c} is the number of spiral turns. For the number of turns typical in this case: N_{c} ~ 10, displacement compounds are a quantity of order 2 Å. It qualitatively corresponds to observable values. Asymmetrical excitations have the lowest energy. Therefore, most likely, they are realized in enzymatic proteins. It was shown that at this type of excitation, the bending of the αhelix is formally directed to the opposite side with respect to the antisymmetric excitations. Also, it has a greater value than the antisymmetric case for N_{c} ≤ 14 and smaller for N_{c} > 14.
PACS
92C05
MCS
36.20.Ey
Keywords:
Alphahelix; Protein; NanotubeBackground
Hydrolysis of ATP and amide I excitation
A protein molecule has a rather unique structure not only in the chemicalbiological point of view but also as an interesting physical and mathematical object. If we consider it as a physical object, then such object may be referred to as a nanostructure without any doubt. Thus, the alphahelical region of a protein molecule simultaneously may be considered both as a nanotube and as a nanowire: this depends on the considered level of structure.
Here, the alphahelix is considered at the level of secondary structure where it is a nanotube. It is in the conditions of quantum excitation which is stimulated by reaction of hydrolysis of adenosine triphosphate (ATP). As a result of this reaction, energy in the form of quanta of infrared range is released. It is considered that they are absorbed by a group of energy states known in an alphahelix as amide I, etc. It is considered also that these absorbing states have an internally molecular oscillating nature. The results obtained here allow giving a definite answer to this question, because in the infrared range, absorption can also have the nature of electronic transitions between states with the main quantum number equal to 2.
The alphahelix is interesting as a mathematical object too. Due to the high sensitivity of its ‘crystalline lattice’ in relation to excitation, we are coming to a necessity to solve a nonlinear system of the socalled eigen type, i.e., actually, we are coming to a necessity to search for the eigenvalues and eigenvectors of a nonlinear system of algebraic equations. Such a problem, as it is known to us, is a scantily explored mathematical problem.
Figure 1 shows the alphahelical fragment of a protein molecule. Similar regions in proteins are widespread enough in vivo. The degree of helicity in different proteins varies from 12% to 96%. As can be seen from Figure 1, the alphahelical fragment of protein molecules is structurally a nanotube. The same is true for its physical properties. Therefore, to such regions of protein molecules in their excited states, it is natural to apply methods that are specific for nanotubes.
As a result of hydrolysis of ATP molecule, energy is realized in the range 0.2 to 0.4 eV^{a}. It depends on the charge state of the ATP molecule, in which the composition of the environment influences mainly (pH, etc.). The energy of hydrolysis is absorbed by an alphahelical region of the protein molecule. It takes place due to internal vibrational excitations of the peptide groups (HCNO) in the state amide I. Its energy is also varied within the limits of 0.2 to 0.4 eV. These excitations induce a significant increase of dipole moments of the peptide groups, which is equal to 3.7 D, on 0.29 D[4,5].
There exists another point of view. Excitation of amide I may have an electronic nature. It may correspond to transitions between energy bands with principal quantum numbers that are equal to 2. The physical nature of excitation is inessential for further calculations, but further it will be shown that their nature may be determined experimentally.
Methods
Amide I excitation in the simplest model of alphahelical region of protein
Foremost, we need to determine the model of description of the spatial structure of the alphahelix. Since it is considered as a molecular crystal, the nearest neighbor approximation is used, which is typical for such crystals. However, as seen from Figure 1b, the nearest neighbors for some peptide group with number n are not only group n ± 1 but also group n ± 3.
The simplest model of the spatial structure of the alphahelix is shown in Figure 2. Such simplified model differs from a real molecule only by symmetry. In the model considered, the molecule is independent from each other: translational and axial symmetries. The real molecule has translationalhelical symmetry. Preliminary investigations have already shown that the qualitative picture in terms of types of excitation does not change. Changes will only be quantitative. They will lead only to some displacement of the absorbing states by energy and, may be, to some mixing of states. In this, simplest, model, all turns of the helix closed on itself, although Figure 1 shows that this is not quite so. Each turn of the helix is open for the nearest neighbor. It was previously shown [6] that taking into account open individual cells leads only to quantitative changes. The qualitative picture remains unchanged.
Figure 2. Simplest model of alphahelix as a onedimensional molecular crystal with three molecules per unit cell. Arrows are showing a separate peptide group. They symbolize the dipole moments.
Within the framework of the considered model, every three peptide groups that belong to one turn of the helix grouped into one complex unit cell. We will number these unit cells by indices n, m, etc. The number of such cells is three times less than the number of peptide groups, i.e., N_{0}/3. Peptide groups within a single cell will be enumerated by indices α, β, etc. that may take values 0, 1, 2. The general functional for the alphahelix in this model has the form [7]
w(R_{nα} − R_{mβ}) in this functional is the basic energy of interaction between peptide groups nα and mβ. It is independent on the presence of excitation and exists always. D(R_{nα} − R_{mβ})A_{αn}^{2} is an additional energy to the w(R_{nα} − R_{mβ}) energy of interaction related only to excitation but considerably smaller. Factor A_{αn} is the wave function that describes the excited state of the examined alphahelical region of the protein molecule. It determines the spatialtemporal distribution of excitation in this region. The energy D(R_{nα} − R_{mβ})A_{αn}^{2} leads to the breaking of the equilibrium of the alphahelix and stimulates its conformational response to excitement. Energy is also an additional energy of interaction. However, it is much less than D(R_{nα} − R_{mβ})A_{αn}^{2} but important because it provides the propagation and transfer of excitation along the alphahelix.
As shown in Figure 2, the nearest neighbors for some peptide group nα will only be the peptide groups m = n ± 1, β = α and m = n, β = α ± 1. Taking into account that in the considered model all energy terms depend on the distances between amino acid residues only, the following formulae in the nearest neighbor approximation may be obtained: R_{nα} ≡ R_{n + 1,α} − R_{n,α}, ρ_{nα} ≡ R_{n,α + 1} − R_{n,α}.
Let us take into account that the response of the lattice (Figure 2) on excitation inside of the unit cell is small enough. Thus, it may be neglected in comparison with a similar response between unit cells. In this sense, the equality ρ_{nα} = ρ_{0} is always supposed fulfilled. Factor R_{nα} is the only value that takes into account the response of the alphahelix on excitation. Thus, we will denote its equilibrium value as R_{0}. Values ρ_{0} and R_{0} are shown in Figure 2. Taking into account the normalization condition
the last functional takes the form
Here, w_{⊥} ≡ w(ρ_{0}), D_{⊥} ≡ D(ρ_{0}), M_{⊥} = M(ρ_{0}), and M_{} = M(R_{0}). Obviously, M_{⊥} ≠ M_{}. If resonance interaction has no electronic nature, inequality will be realized: M_{⊥} < M_{}. If excitation has an electronic nature, inequality will be reversed: M_{⊥} > M_{}. This difference may be detected experimentally, and the answer of the question about the physical nature of excitation may be obtained.
New equilibrium values of distances, which actually coincide with the step of alphahelices, are determined using the general condition of minimization: . When interactions between peptide groups are modeled as purely dipole, the step of the alphahelix always decreases and is given by
Next, we must substitute (3) in (2), take into account the condition , designate w(R_{0}) ≡ w_{}, D(R_{0}) ≡ D_{}, , and introduce convenient redesignation: M_{} = −M_{} ≡ −2Λ, M_{⊥} = M_{⊥} ≡ 2Π, which take into account the true signs. Then for the functional (2), finally, the following formula will be obtained:
In Equation 4, E_{осн} = (w_{⊥} + w_{})N_{0} + D_{⊥} + D_{}, and the following is taken into account:
N_{0} is the number of amino acid residues in the alphahelical region of the protein molecule, which is under consideration.
Further, for implementation of the conditional minimization of energy (4) in relation to wave functions A_{αn}, it is necessary to create a conditional functional: . From a mathematical point of view, parameter ϵ is an indefinite Lagrange multiplier, and physically, it is the eigenvalue of the considered system. The minimization procedure produces the equation Λ(A_{α,n + 1} + A_{α,n − 1}) + GA_{αn}^{2}A_{αn} − Π(A_{α + 1,n} + A_{α − 1,n}) + ϵA_{αn} = 0. After dividing this equation by Λ and introducing the notations,
it is possible to reduce it to a dimensionless form:
The function A_{αn} is complex. Therefore, the common solution of the system (6) has the form A_{αn} = a_{αn} · exp(iγ_{αn}). Amplitude a_{αn} and phase γ_{αn} are real functions of the variables α and n. We confine ourselves to the HamiltonianLagrangian approximation in phase [8]. Due to the stationarity of the solved problem, this approximation has the simplest form: γ_{αn} ≡ kn. If the alphahelical part of the molecule is long enough,^{b} a BornKarman condition gives . Here, is the number of turns in the considered alphahelical region of the protein molecule. It plays the role of the dimensionless length of the helical region of the protein in units of an alphahelix step. Parameter j has the values . Then
and Equation 6 takes the form
Separating real and imaginary parts, we have the following formulae:
The solution of this system is usually determined after transition to continuous approximation. But we will analyze systems (8) and (9) without using the continuous approximation, because we are interested in very short alphahelical regions (10 to 30 turns).
There is only condition a_{α,n + 1} − a_{α,n − 1} = 0 (if not to restrict solutions by using the condition k = 0), which does not depend on any symmetry of the alphahelix: whether it is the symmetry of the model or the symmetry of the real molecule. Viewing of other conditions can appear useful on account of the real structure of the alphahelical region. In the simplest case, it may be reduced to the equation a_{αn} = P_{α}. The system (8) now degenerates in the system of three nonlinear equations:
where the following designations are introduced:
The last, fourth, equation arose out from normalization condition (1). The coefficients P_{α} (α = 0, 1, 2) determine the excitement of each peptide chain as a whole.
The system (10) consists of four nonlinear equations for determining the values P_{0}, P_{1}, and P_{2} and the eigenvalue x. By adding and subtracting the first two equations and some transformation of the third equation, the system (10) can be reduced to the form
This transformation does not affect the solutions of the system.
For the solution, the condition P_{0} + P_{1} = 0 should be used. This condition together with the condition P_{2} = 0 turns into an identity the second and third equations. After some simple transformations, we obtain the antisymmetric excitations:
Using Equations 4, 5, and 11, it is possible to find the energy:
Next, we use the condition P_{0} − P_{1} = 0, which turns into an identity the first equation in (12). After some analysis, we can find two types of excitation:
● Symmetrical
● For these excitations, in analogy to the antisymmetric, it is possible to obtain the energy:
● Asymmetrical
● For these excitations, it is also possible to get energy:
The energies E_{a}(k), E_{c}(k), and E_{н}(k) contain parameters Λ = M_{}/2 and Π = M_{⊥}/2. As it was noted between Equations 2 and 3, the relation between these parameters makes the determination of the physical nature of excitation possible: whether they are electronic or intramolecular. Because one of them (Λ) determines the width of the excited energy bands, and the other (Π) their positions, this is the basis for the experimental analysis of the nature of excitations.
There are a few possibilities else for searching for solutions of the system (12). Preliminary analysis shows that the obtained excitations are peculiar in a more or less degree for both symmetries: whether it is the symmetry of the model or the symmetry of the real molecule. The other solutions of the system (12) need to be analyzed only in the conditions of the maximum account of the real structure of an alphahelix. But the general analysis of this system shows that the solutions of a new quality are not present: all of them belong to the asymmetrical type. However, attention should be paid to the equation a_{α,n + 1} − a_{α,n − 1} = 0, which has led to the requirement a_{αn} = P_{α}. This condition is strong enough and essentially limits the solution: it is a constant in variable n, i.e., does not have the spatial distribution along an alphahelix.
Results and discussion
The analysis of the energetics of the protein excitation
From definitions (13), (14), and (15), it ensues that received excitations are located in accordance with the inequality E_{c}(k) > E_{a}(k) > E_{н}(k). Thus, and . It can be seen that for the alphahelical region of finite length, when the number of turns N_{c} ≠ ∞, the lowest energy is the energy of asymmetric excitation E_{н}. Also, it is visible that energy E_{c} is always strongly separated from energies E_{a} and E_{н}. Even when the number of turns N_{c} ⇒ ∞ and the energies E_{a} and E_{н} practically coincide, the energy E_{c} is separated from E_{a} and E_{н} on a value 3Π = 3M_{⊥}/2. Amide I excitations manifested experimentally are probably E_{c} energy.
It is possible to make the supposition that each of the examined energies executes some, expressly certain, function. For example, the main function of symmetric excitations can be activation of muscle proteins. At the same time, they can activate both membrane and enzymatic proteins that are quite often actually observed in the activation of myosin [911].
Antisymmetric excitation energy is not enough to excite the muscle protein because it lies below the symmetric energy. Activation of membrane proteins can be their main function. At the same time, these excitations are able to activate enzymatic proteins that are also actually observed often enough during activation of membranes [1113].
And, lastly, asymmetrical excitations have only one function  to activate exceptionally enzymatic activity in those cases, when membrane and muscular activities are not needed. That is only for intracellular processes.
Conformational response to the excitation of the alphahelical region of the protein molecule
For the analysis of conformational response of the alphahelix on the considered excitations, it is necessary to appeal again to new equilibrium values of the step of the alphahelix. From definition (3), it is possible to find R_{nα} = R_{0} · (1 − βA_{αn}^{2}), where designation is entered: . If we consistently apply the model of dipole interaction between the peptide groups, then , where, as mentioned above, Δd ~ 0.29 D and d ~ 3.7 D. Therefore, in this dipole model [14], β ~ 10^{−1}. Taking into account the definitions of coefficients A_{αn}, given in (7), it is possible to get following:
1. It is possible to obtain the following formula for symmetric excitations: . That is, all three chains are reduced equally and evenly in the space. Then the length of every peptide chain can be appraised, so
This change is small and, at first glance, has no practical significance. But it will be so only in the classical model of the alphahelix (Figure 2). If we consider, for example, that the peptide chains of myosin themselves form superhelices, then the effect of contraction increases. This is done by changing all characteristics of an alphahelix: the step of the helix, its radius, and the effective number of peptide groups on the turn of the helix. Also, additional selftorsion takes place. The strengthening of the effect of contraction is determined by the mutual torsion of long alphahelical regions of light faction of myosin and their torsion on actin filaments.
2. For antisymmetric excitations, it is possible to obtain , . Respective lengths are as follows:
In this type of excitation, one of the peptide chains does not change (here, it is a chain with the number 2), and two others are reduced up to the value . Such asymmetry is enough for the alphahelix to take a form of the segment of torus instead of cylinder (Figure 3). Application of the simple geometric considerations gives for the radius of curvature R_{k} and angle φ:
and for displacement Δ, it is possible to get such estimation:
Taking into account the numerical values β ~ 10^{−1}, R_{0} = 5.4 Å, and d_{α} = 4.56 Å in (16) gives . For the typical number of turns in many enzymes and membrane squirrel (N_{c} > 10), displacement will have an order Δ > 2 Å. This is consistent with the observed values [11].
3. For asymmetrical excitation, the following values are implemented: , . The corresponding lengths of peptide chains equal
The nature of the distribution of deformation along the peptide chain for this type of excitation is similar to that of the antisymmetric excitation. The only difference is that the chain, which in the previous case has not changed at all, now has shortening stronger than the other two. It is possible to estimate displacement for this case too:
Here, Δ is the displacement for antisymmetric excitations, which is determined by Equation 16. Unlike displacement Δ, displacement Δ^{(н)} ‘directed’ to the opposite side. Executing numerical estimates, it is possible to set that Δ^{(н)} > Δ, if the number of turns in the alphahelix N_{c} ≤ 14, but at N_{c} > 14, we will have Δ^{(н)} < Δ accordingly.
Consequently, asymmetrical excitations demonstrate two very interesting features. First, it has the lowest energy and at diminishment of the number of turns N_{c}, it falls down yet more. Second, a conformational response for this type of excitation is the biggest for N_{c} ≤ 14. This is typical for enzymatic proteins only.
Figure 3. Explanation to estimation of displacement Δ of free (here upper) end of alphahelix for antisymmetric excitations.
Conclusions
The general methods [7,1517] of description of the excited states of the condensed environments were applied to the alphahelix region of a protein molecule. The alphahelix is considered as a nanotube, and excitations of the environment are described as quasiparticles. It is shown that three different types of excitation exist, and each of them is probably used by three different types of protein. The symmetrical type of excitation is used for muscle proteins, the antisymmetric type of excitation is used for membrane proteins, and the asymmetric type of excitation is used for enzymatic proteins. It is possible that some excitations of asymmetrical type exist, which are also used by enzymes. The estimations were done for displacements of the free end of the alphahelix. The obtained displacements are in agreement with experimental data. Therefore, the obtained results can be the basis of the interpretation of the functional properties of proteins characterizing their activity related to their conformational changes [11].
Endnotes
^{a}Offsystem unit of energy: 1 eV = 1.602 × 10^{−19} J.
^{b}For example, in the myosin protein, the helical region has about 200 turns or up to 700 amino acids.
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
AD Suprun. He formulated a scientific problem, analyzed the results and take part in the discussing and formulating the conclusions. LV Shmeleva. She made mathematical calculations, take part in the discussing of the results and conclusions. Both authors read and approved the final manuscript.
References

Pauling L, Corey RB, Hayward R: The structure of protein molecules.

Kendrew JC: The threedimensional structure of protein molecule.
Sci Amer 1961, 205(6):96111. PubMed Abstract  Publisher Full Text

Davydov AS, Suprun AD: Configuration changes and optical properties of αhelical protein molecules.
Ukrainian J Phys 1974, 19(1):4450.
(in Russian)

Yu N, Chirgadze E, Rashevskaya P: Intensity of characteristic vibrations of peptide groups.
Biophysics 1969, 14(4):608614.
(in Russian)

Rick SW, Cachau RE: The nonplanarity of the peptide group: molecular dynamics simulations with a polarizable twostate model for the peptide bond.
J Chem Phys 2000, 112(11):52305241. Publisher Full Text

Suprun AD, Atmazha YB: Quantum excitation of protein αspiral and the problem of protein functionality.

Suprun АD: Dynamic Properties of SingleElectron Nonlinear Excitation of the Crystals. Kyiv: Kyiv University; 2008.
(in Ukrainian)

Suprun AD, Shmeleva LV: Degeneracy effect of dynamical properties of quasiparticles of electronic origin in semiconductor materials.

Engelgardt WA, Lubimova MN: Myosin and adenosine triphosphatase.

Hachikubo Y, Ito K, Schiefelbein J, Manstein DJ, Yamamoto K: Enzymatic activity and motility of recombinant Arabidopsis myosin XI, MYA1.
Plant Cell Physiol 2007, 48(6):886891. PubMed Abstract  Publisher Full Text

Davydov AS: Biology and Quantum Mechanics. Kiev: Naukova Dumka (Scientific Thought); 1979.
in Russian

Skon JC: The influence of some cations on an adenosine triphosphatase from peripheral nerves.

Mouritsen OG, Andresen TL, Halperin A, Hansen PL, Jakobsen AF, Jensen UB, Jensen MO, Jørgensen K, Kaasgaard T, Leidy C, Simonsen AC, Peters GH, Weiss M: Activation of interfacial enzymes at membrane surfaces.
J Phys Condens Matter 2006, 18:12931304. Publisher Full Text

Liwo A, Pincus MR, Wawak RJ, Rackovskyp S, Scheraga HA: Calculation of protein backbone geometry from αcarbon coordinates based on peptidegroup dipole alignment.
Protein Sci 1993, 2:16971714. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Natanzon Y, Brizhik LS, Eremko AA: Dynamics of a selftrapped quasiparticle in a onedimensional molecular lattice with two phonon modes.
Phys Status Solidi B 2007, 244(2):545554. Publisher Full Text

Scott AC: Dynamics of Davydov solitons.
Phys Rev A 1982, 26(1):578595. Publisher Full Text