Mechanical Models of Microtubules

Microtubules are the major part of the cytoskeleton. They are involved in nuclear and cell division and serve as a network for motor proteins. The first model that describes nonlinear dynamics of microtubules was introduced in 1993. Three nonlinear models are described in this chapter. They are longitudinal U-model, representing an improved version of the first model, radial φ-model and new general model. Also, two mathematical procedures are explained. These are continuum and semi-discrete approximations. Continuum approximation yields to either kink-type or bell-type solitons, while semidiscrete one predicts localized modulated waves moving along microtubules. Some possible improvements and suggestions for future research are discussed.


Introduction
A cell is defined as eukaryotic if it has a membrane-bound nucleus.Such cells are generally larger and much more sophisticated than prokaryotic ones.Microtubules (MTs) are the basic components of cytoskeleton existing in eukaryotes [1].They are long structures that spread between a nucleus and a cell membrane.MTs play an essential role in the shaping and the maintenance of cells and are involved in cell division.Also, they represent a network for motor proteins.These proteins move with a velocity of 0:1 À 2μm=s [2] carrying a certain cargo such as mitochondrion.
All eukaryotic cells produce two kinds of tubulin proteins.These are α and β tubulins, or monomers, and they spontaneously arrange head to tail forming biologically functional subunit that we call a heterodimer, or a dimer for short.When intracellular conditions favor assembly, the dimers assemble into long structures called protofilaments (PFs).Microtubules are usually formed of 13 PFs, as shown in Figure 1.
Hence, MTs are long cylindrical polymers whose lengths vary from a few hundred nanometers up to meters in long nerve axons [4].Each dimer is an electric dipole whose mass and length are m ¼ 1:8 Â 10 À22 kg and l ¼ 8 nm, respectively [5].The component of its electric dipole moment in the direction of PF is p ¼ 337Debye ¼ 1:13 Â 10 -27 Cm [6].Consequently, MT as a whole appears to be a giant dipole with negatively charged end coinciding with biologically positive end (more active) and vice versa.This is the reason why an intrinsic electric field exists within MT.
MTs in non-neuronal cells are unstable structures.They exhibit dynamic instability behavior existing in phases of elongation (polymerization) or rapid shortening (depolymerization).This size fluctuation has been called as dynamic instability [7,8].Notice that the shrinkage rate is bigger than the growth rate (see Ref. [9] and references therein).MTs grow steadily at positive end, corresponding to the β-subunit, and then shrink rapidly by loss of tubulin dimers at the negative end, corresponding to the α-monomer.Many anticancer drugs, for example, taxol (paclitaxel), prevent growth and shrinkage of MTs and thus prevent cell proliferation [10].
MTs existing in neuronal cells are stable and, consequently, neurons, once formed, do not divide [4].This stability is crucial as there are evidences that neuronal MTs are responsible for processing, storage and transduction of biological information in a brain [4,11].
It was mentioned that MTs represent the traffic road for motor proteins.Some more information can be found in Ref. [9] and in an exhaustive review paper [12].It suffices now to state that the cellular motors with dimensions of less than 100 nm convert chemical energy into useful work.These small machines have the fundamental role of dissipation in biological systems, which has been confirmed by both the theoretical and the experimental investigations [13].The molecular motors dissipate continuously and operate as irreversible systems [13].
It is clear that any molecular motor, to start moving, should obtain a certain signal.One of the promising dynamical mechanisms for intracellular signaling is solitary waves, which is explained in this chapter.

Mechanical models
MTs, as well as all biological systems, are nonlinear in nature.Strong covalent chemical bonds are usually modeled by linear "springs", while weak chemical interactions, existing in all biological Figure 1.A tubulin dimer, a protofilament and a microtubule [3].
Complexity in Biological and Physical Systems -Bifurcations, Solitons and Fractals systems, are modeled by nonlinear "springs".This means that expressions for energy of biological systems require nonlinear terms, which brings about nonlinear partial differential equations (PDEs) explaining nonlinear dynamics of these systems.This is the topic of the present chapter.We will see that, in case of MTs, the solutions of these nonlinear PDEs are solitary waves.
The word soliton was introduced in 1965 to designate solitary waves describing the propagation of excitations in continuous media with nonlinearity and dispersion [14].The first qualitative description of solitary waves dates back to 1834 when hydrodynamic engineer John Scott Russell observed them on a surface in a shallow channel [15].The wave was so stable that the engineer followed it about 1 or 2 miles.From then, there has been tremendous interest for various kinds of solitons in many branches of physics [15][16][17][18][19].In this chapter, the terms soliton and solitary waves are treated as synonyms, which is commonly accepted in literature.
Solitons are localized waves possessing some interesting properties.The most important is their stability in a sense that they conserve their shape and energy after mutual interaction.In other words, they can pass through one another without annihilation.This was experimentally observed in neurons [20].
To model complex MT dynamics, we should introduce some simplifications.To the best of the author's knowledge, all the models introduced so far have only one degree of freedom per dimer.Hence, for the models explained in this chapter, elementary subunits of PFs are dimers and they perform either longitudinal or angular oscillations and the appropriate models can be called as longitudinal or angular (radial), respectively.
The longitudinal contacts along PFs are much stronger than those between adjacent PFs [21,22], which allows us to construct a simplified Hamiltonian of MT, which is, practically, Hamiltonian for a single PF only.However, the influence of the neighboring PFs is taken into consideration through the electric field.Namely, each dimer exists in the electric field coming from the dimers belonging to all PFs.Also, the nearest neighbor approximation is assumed.

U-model
The first model that describes nonlinear dynamics of MTs is a longitudinal one.It was introduced in 1993 by Satarić et al. [23].According to the model, the dimers perform angular oscillations but a coordinate u, describing the dimer's displacement, is a projection of the top of the dimer on the direction of PF.Therefore, the displacements are radial but the used coordinate is longitudinal.There is a real longitudinal model assuming longitudinal displacements of the dimers that we call as Z-model [24].Both U-and Z-models bring about equal crucial differential equations and the latter one will not be studied here.
Somewhat improved and more general version of the first nonlinear model is what we call as U-model [25] and this will be explained in the following paragraphs.Both models are based on the fact mentioned above that the dimers are electric dipoles and that the whole MT can be regarded as ferroelectric [23,26], which means that the interaction between a single dimer and its surrounding can be modeled by W-potential [23,27].This yields to the following Hamiltonian for MT [23,25,28] where dot means the first derivative with respect to time while the integer n determines the position of the considered dimer in PF.The first term obviously represents a kinetic energy of the dimer of mass m.The second one is interaction between the neighboring dimers belonging to the same PF in the nearest neighboring approximation and k is an intradimer stiffness parameter.The next two terms represent the W-potential energy mentioned earlier, where the parameters A and B should be determined or, at least, estimated and are assumed to be positive.We should point out that the double-well potential is rather common in physics [27,29,30].The very last term is coming from the fact that the dimer is the electric dipole existing in the field of all other dimers, where Q > 0 represents the excess charge within the dipole and E > 0 is internal electric field.The last three terms together can be regarded as unsymmetrical W-potential.
Our final goal is the function u n t ðÞ , describing nonlinear dynamics of MT.This function is a solution of so-called dynamical equation of motion, which can be obtained from Eq. (1).To derive it, we introduce generalized coordinates q n and p n defined as q n ¼ u n and p n ¼ mdu n =dt.Using well-known Hamilton's equations of motion dp n =dt ¼ÀdH=dq n and dq n =dt ¼ dH=dp n , we obtain the following discrete differential equation that should be solved The last term is a viscosity force with γ being a viscosity coefficient [23].Therefore, nonlinear dynamics of MTs has been described by Eq. (2).Obviously, nonlinearity is coming from the fourth degree term in the W-potential.
It was explained earlier that we used some approximations to derive Eq. ( 2).However, we need one more to solve it.We now explain two mathematical methods for solving this equation.Practically, these two approaches are two approximations.They are continuum and semidiscrete approximations.We will see that the different mathematical procedures yield to different solutions.Therefore, the function u n t ðÞdepends not only on the physical system but also on the used mathematical method.
Let us explain the continuum approximation first.A question if MTs are discrete or continuum systems was studied in Ref. [31], where it was shown that the continuum approximation is valid.The continuum approximation means a transition u n t ðÞ!ux ; t ðÞ , which allows a series expan- w h e r el is the dimer'sl e n g t he x p l a i n e d earlier.In fact, PF can be seen as one-dimensional crystal with l being a period of the lattice.This straightforwardly brings about the following continuum dynamical equation of motion This is PDE that cannot be easily solved.Hopefully, this equation can be transformed into an ordinary differential equation (ODE).It is well known that, for a given wave equation, a traveling wave u ξ ðÞ is a solution which depends upon x and t only through a unified variable ξ as ξ ¼ κx À ωt, where κ and ω are constants.If we substitute the variables x and t by ξ we straightforwardly transform Eq. ( 3) into the following ODE where u 0 du=dξ and Eq. ( 4) becomes the appropriate one in Ref. [23] for α u ¼À1.Therefore, the U-model is more general than its predecessor introduced in Ref. [23].It is crucial that the parameter α u can be determined together with the function Ψ for known or estimated σ and r u .This is because α u has very important physical meaning.The first term in Eq. ( 3) is the inertial term and it is coming from the kinetic energy in Hamiltonian (1), while the second one is the elastic one.Therefore, positive α u means that the inertial term is bigger than the elastic one and vice versa.
Eq. ( 4) has already been solved using different mathematical procedures like standard procedure [23,27,29,30] and method of factorization [31,32].There exists a group of procedures where the function Ψ is represented as a serious expansion over other known function like The function Φ is usually known and we plug Ψ into Eq.( 4) and determine the coefficients A k .A common example for Φ is a solution of Riccati equation, which is either tangent or tangent hyperbolic.As only the latter function may have physical meaning, we call the method as tangent hyperbolic function method (THFM) [25,[33][34][35] and extended or modified extended THFM [36].The function Φ can also be one of Jacobian elliptic functions [37] and, even, unknown [38].
where A 0 , A k and B k are coefficients that should be determined and Φ 0 represents the first derivative.In general, the function Φ ¼ Φ ξ ðÞis known and represents a solution of a certain ODE of lower order than the equation that should be solved.A commonly used example is the Riccati equation [40] To determine the positive integer N in Eq. ( 6), we should plug Ψ ¼ c=ξ p , c ¼ const, into Eq.( 4) and concentrate our attention on the leading terms [42].One can easily show that N ¼ 1 for Eq. ( 4) as the leading terms are proportional to ξ À pþ2 ðÞ and ξ À3p .
Our next step is determination of the parameters A 0 , A 1 , B 1 , a, b and α u .According to Eqs. (6) for N ¼ 1 and ( 7), we obtain the expressions for Ψ 0 , Ψ 00 and Ψ 3 as required by Eq. ( 4), which yields to the following expression: Obviously, this is satisfied if all the coefficients are simultaneously equal to zero.This brings about a system of seven equations, which can be obtained using Mathematica or similar software [39].One of them can be written as indicating two possible relationships between the parameters A 1 and B 1 .Hence, there are a few cases to be studied.They are as follows [39]: It is obvious that the first case represents nothing but a simpler method called extended tanhfunction method.The system mentioned earlier brings about [39] 8A The final result is [39] where A 0i is the following three real solutions of the first of Eqs. ( 11) Of course, these three solutions exist for σ < σ 0 .The case σ > σ 0 was discussed in Ref. [25].
All the three solutions are shown in Figure 2 for σ ≈ 0:9σ 0 and r u ¼ 1.Of course, these solutions reproduce previously known results [25].Figure 2 shows that the solutions of Eq. ( 4) are kink and antikink solitons.More detailed analysis of their physical meaning is given in Ref. [25].
It was shown [42] that the second case is equal to the first one indicating that the value of a is irrelevant if B 1 ¼ 0. In other words, we could have assumed a simpler version of the Riccati equation neglecting the term 2aΦ.
The third case is more interesting.It turns out that, instead of the three lines in Figure 2, that is, the three solutions, we obtain infinitely many lines corresponding to each of them [39].However, they represent three groups of parallel lines, which means that all these solutions are only shifted functions and, consequently, have equal physical meaning.Therefore, this case does not bring about any physically new result.
Case 4 is suggested by Eq. (10).The system of seven equations, mentioned earlier, gives the first and the last term in Eq. ( 11) as well as The final expression for Ψ is This case yields to a new solution, which was not obtained using less general mathematical methods.However, it may be interesting from mathematical point of view only as Ψ diverges for ξ ¼ 0. Case 5 is a simplified version of SEM, explained in Ref. [42].The mentioned system brings about r u ¼ 0 as well as where a notation A 0 has been introduced to distinguish this parameter from A 0 , used in the previous cases.It is interesting to compare the polynomials for A 0 and A 0 , existing in Eqs.(11) and (17).We can see that.
which means that the values for A 0i are given by Eqs. ( 13), ( 14) and (18).We can easily show that the final solution for Ψ is [39] Obviously, this function cannot diverge for any value of aξ but only for À1 < K < 1.Also, K should be real and these two requirements eliminate Ψ 2 and Ψ 3 [39], which means that Ψ and A 0 in Eq. ( 19) are Ψ 1 and A 01 .
The function Ψ 1 ξ ðÞis shown in Figure 3 for a ¼ 0:1 and for two values of the parameter σ.W e notice very interesting result that is a bell-type soliton!This certainly demonstrates the advantage of SEM method over the less general ones.
It is important to study the physical meanings of the parameters a and σ.Eq. (19) indicates that solitonic width is inversely proportional to a jjand that a does not affect maximum of the wave.Figure 3 shows that the amplitude of Ψ 1 is a decreasing function with respect to σ.
Finally, the last case gives the solution Therefore, all the cases are explained and we can see that the continuum approximation yields to both kink solitons and bell-type solitons.The latter may exist only if viscosity is neglected.
It was stated earlier that the coordinate u was the projection of the top of the dimer on the direction of MT.A patient reader may ask how u can be negative when this is the projection.This question is answered in Ref. [28].
It was mentioned earlier that there are two approximations that can be used to solve Eq. ( 2).Now we get back to Eq. ( 2) and study semi-discrete one [15,28,43].A mathematical basis for the method is a multiple-scale method or a derivative-expansion method [16,44].We assume small oscillations u n t ðÞ ¼ε Φ n t ðÞ , ε << 1 (21) which straightforwardly transforms Eq. ( 2) into According to the semi-discrete approximation, we look for wave solution which is a modulated wave, that is [28,45] where ω is the optical frequency of the linear approximation, q =2π/λ > 0 is the wave number, cc represents complex conjugate terms and the function F 0 is real.Of course, l is the dimer's length, as mentioned earlier.The function F is continuous and represents an envelope, while exp(iθ n ), including discreteness, is a carrier component.Notice that the parameter ε exists in Mechanical Models of Microtubules http://dx.doi.org/10.5772/intechopen.71181 the function F, but does not in exp(iθ n ).This is because the frequency of the carrier wave is much higher than the frequency of the envelope and we need two time scales, t and ε t, for those two functions.The same holds for the coordinate scales.
To simplify the problem, a continuum limit nl !z should be introduced as well as new transformations Z ¼ εz and T ¼ ε t.This allows a series expansion of F ξ ðÞ , that is where indexes Z and ZZ denote the first and the second derivative with respect to Z. Hence, the function Φ n t ðÞbecomes where * stands for complex conjugate and F FZ ; T ðÞ .All this allows us to obtain the expressions existing in Eq. ( 22), such as Φ nþ1 þ Φ nÀ1 À 2Φ n , Φ n and Φ n 3 , and Eq. ( 22) becomes [28] This crucial expression represents a starting point for a series of important expressions.They can be obtained equating the coefficients for the various harmonics.For example, equating the coefficients for e iθ and neglecting all the terms with ε 2 and ε 3 one obtains the following expressions for the dispersion relation ω ¼ ω q ðÞand the group velocity dω=dq: Also, the coefficients for e i0 ¼ 1 and e i2θ , respectively give [28] which yields to Eqs. ( 28) and ( 29) and new coordinates S and τ, defined as S ¼ Z À V g T and τ ¼ ε T, allows us to simplify Eq. ( 27).An explanation for why the parameter ε exists in the time scaling but is absent in the space scaling is given in Refs.[45,46].If we consider the terms for e iθ again we obtain the well-known nonlinear Schrödinger equation (NLSE) for the function where the dispersion coefficient P and the coefficient of nonlinearity Q are Even though Eq. ( 31) is PDE, its solution exists.This well-known solution, existing for PQ > 0, is [15,47,48] FS ; τ ðÞ ¼ A e sech S À u e τ L e exp iu e S À u c τ ðÞ 2P ,u e > 2u c (33) where parameters u e and u c represent envelope and carrier component velocities, while the amplitude A e and the soliton width L e have the forms It is very difficult to deal with the parameters u e and u c as u e > 2u c is completely unprecise statement.However, u c =u e < 0:5 seems to be more practical.Hence, new parameters U e and η have been introduced as U e ¼ εu e , η ¼ u c =u e and 0 ≤ η < 0:5 [45].Finally, we can easily obtain the expression for the longitudinal displacement of the dimer at the position n where One more parameter can be eliminated using the idea of coherent mode [49].This mode means that the envelope and the carrier wave velocities are equal.It follows from Eq. ( 35) that V e ¼ Ω=Θ, which yields to the function U e η ðÞ .This means that the wave u n t ðÞis the one phase function, preserving its shape in time.
To plot the function u n t ðÞor, equivalently, U n t ðÞwe should know or estimate the values of a couple of the parameters.Of course, if 2D plot is chosen, U n t ðÞcan be presented as either a Mechanical Models of Microtubules http://dx.doi.org/10.5772/intechopen.71181function of t at a certain position n or as a function of n for chosen t.Very detailed analysis of the parameter selection was done in Ref. [28].One example for q ¼ 2π=Nl is shown in Figure 4. Obviously, this is a localized modulated wave usually called as breather.We can see that its width is about 200 nm, which means that it covers about 25 dimers.
As a conclusion, we can state that the two mathematical procedures bring about even three results, that is, three different solitons.These are kinks, bell-type solitons and breathers.They may be signals for the motor proteins to start moving, as explained in Introduction.
Obviously, viscosity has been neglected.This will be explained in the following section, within the φ-model.

φ-model
A weak point of the U-model is the last term in Eq. (1).A scalar product À p ! Á E ! ¼ÀQdE cos φ n would be better choice for the potential energy, where d is the distance between the centers of positive and negative charges within the dipole.This potential indicates the angle as a coordinate instead of the projection u and the Hamiltonian for the radial model, which we call as φ-model, is [50,51]  Complexity in Biological and Physical Systems -Bifurcations, Solitons and Fractals where I is a moment of inertia of the dimer at the position n.Notice that the W-potential does not exist in Eq. ( 38) even though the terms including φ n 2 and φ n 4 appear as a result of a series expansion of the cosine function.Instead of viscosity force introduced in the previous section, we introduce viscosity momentum M v ¼ÀΓ _ φ, where Γ is the viscosity coefficient [51][52][53].Following the procedure explained earlier, we obtain where Like above, the solutions are kink solitons [50].
It is very interesting to compare the expressions for α u and α φ , given by Eqs. ( 5) and (40).They can be written as and where v ¼ ω=κ is the soliton velocity, while c u and c φ are corresponding sound velocities.According to Eq. ( 11), we can see that the U-model predicts c u > v as A is positive.On the other hand, α φ ¼ 2r φ 2 =9 > 0 [50] means that, according to the φ-model, the kink belongs to the class of supersonic solitons.We will return to this issue in the next section.Now, we switch to the semi-discrete approximation within the φ-model to solve the dynamical equation of motion, which is [51] εI where φ ¼ εΦ has been used.Of course, Eq. ( 43) is analog to Eq. (22).Following the procedure explained in the previous section, we straightforwardly obtain F 0 ¼ 0 and F 2 ξ ðÞ¼0, as well as where ω 0 is the lowest frequency of the oscillations [59].Also, we easily obtain NLSE (31), where Mechanical Models of Microtubules http://dx.doi.org/10.5772/intechopen.71181 The final solution φ n t ðÞis the same as U n t ðÞexcept that P and Q are different.Therefore, both the U-and the φ-models predict the breather waves moving through MT.
Finally, viscosity should be introduced in the semi-discrete approximation [51].Due to viscosity momentum M v ¼ÀΓ _ Φ n, the final result φ n t ðÞincludes the expected exponential term e Àβ t , where β ¼ Γ=2I [51].

General model of MTs
It was mentioned earlier that the weak point of the U-model is the last term in Eq. ( 1).Also, it is better to use the radial coordinate φ than the longitudinal one as we assume angular oscillations of the dimers.The scalar product À p ! Á E ! ¼ÀQdE cos φ n , existing in the φ-model, solved these problems but the W-potential has been missing.In fact, a series expansion of cos φ n gives φ n 2 and φ n 4 terms but with opposite signs from those in the U-model.These two terms are, practically, a potential that looks like W in a mirror having only one minimum surrounded by two maxima and, due to its shape, can be called as M-potential [54].This potential brings about α φ > 0, which is disputable result.
Therefore, we want to solve the mentioned problem regarding the U-model but to keep the Wpotential, the coordinate φ and, probably, Iω 2 < kl 2 κ 2 , that is, α < 0. This suggests the following Hamiltonian where A > 0, B > 0 and φ n has the same meaning as in the φ-model.Let us call the model as general one (GM).The procedure mentioned earlier brings about where, of course, φ φξ ðÞ .If we consider Eqs. ( 3) and ( 47), we can see that the last two terms in Eq. ( 47) may be the first derivatives of either W-or M-potential, depending on the sign of the terms in the brackets.However, these brackets may have different signs or can be zero.Therefore, the possible cases are: All of them are studied in Ref. [54] and they will be explained here briefly.
Case 1 straightforwardly yields to where and the final solution is [54] Eq. ( 50) holds for both positive and negative r 1 .Therefore, φξ ðÞrepresents kink soliton if r 1 > 0 and antikink one for negative r 1 , which is shown in Figure 5.
One of the advantages of the GM over the φ-model is the value of amplitude.Namely, the amplitude of the kink soliton, according to the φ-model, is ffiffi ffi 6 p , coming from Eq. ( 40).This is unrealistic, too big value.Instead of ffiffi ffi 6 p , the appropriate factor, existing in the GM, is K, given by Eq. ( 49).
If viscosity is neglected, the GM brings about Case 2 straightforwardly yields to and to the final results and It is obvious that these results do not have physical meaning as φ 2 is complex, while φ 20 may diverge.
The remaining Case 4 linearizes Eq. ( 47) and will not be studied here.
Therefore, the GM yields to the kink solitons as the previous two models do.However, this is the radial model and the problems with both the last term in Eq. ( 1) and the huge amplitude in the case of the φ-model have been solved.We should study one more issue.It was mentioned earlier that the U-model predicts the subsonic kink soliton, while the φ-model predicts the supersonic wave.How about the GM?It was shown that Case 1 yields to the solutions having physical sense and that α 1 > 0. According to Eq. ( 49), we easily reach the final conclusion: a.If A > pE and B > pE=6 , then r 1 < 0 and the function φ x; t ðÞ is subsonic soliton, kink for the positive K in Eq. ( 50) and antikink otherwise.
b.If A < pE and B < pE=6 , then r 1 > 0 and the function φ x; t ðÞ is supersonic soliton, antikink for the positive K in Eq. ( 50) and kink otherwise.
All this certainly suggests the advantages of the GM with respect to the previous two.Complexity in Biological and Physical Systems -Bifurcations, Solitons and Fractals

Conclusion and future research
In this chapter, the three models describing nonlinear dynamics of MTs are shown.The first one, the U-model, is the improved version of the first nonlinear model and it predicts subsonic kink solitons moving along MT.The second one, the radial φ-model, predicts the supersonic kinks.Finally, the GM is explained.This is the radial model which yields both possibilities regarding the kink's speed.If we assume that the kink soliton is subsonic wave then we know the minimum value of the parameter A, that is, A > pE, as explained earlier.
Two mathematical procedures are explained, continuum and semi-discrete approximations.It is very interesting that the final result depends not only on the physical system but on the mathematical methods as well.These solutions are the kink soliton and the breather.The question is which one, if any, really moves along MTs.This is not known in the moment and cannot be without experimental results.
It was demonstrated that the GM is better than the previous two models.However, this does not mean that it should not be improved.For example, there has been an attempt to improve the model introducing Morse potential instead of the harmonic one [55].The harmonic potential energy assumes that attractive and repulsive forces are equal.Morse potential is not symmetric and is good for both strong and weak interactions.
In this chapter, the dimers are considered as elementary units.However, their structure is more complicated and they include tubulin tales (TTs).Consequently, nonlinear dynamics of TTs should also be studied and some results already exist [56,57].
The W-potential has two minima which means that it assumes existence of the two angles between the dimer and the direction of PF around which the dimer oscillates.One of the future tasks should be measuring these angles.First of all, such experiment would check if the Wpotential is correct or not.If it is, then our knowledge of their values would improve the theory a lot.
One of the future research goals should be two-component model.This may mean that we should construct the model assuming two degrees of freedom.However, one of these degrees can be an internal one, which means that oscillations if monomers within the dimer should be taken into consideration.Notice that the two-component model may be the one studying electro-acoustic wave excitations [58].
Finally, we should bear in mind the cytological and medical applications of the research explained in this chapter [59][60][61].