Microtubules as One-Dimensional Crystals: Is Crystal-Like Structure the Key to the Information Processing of Living Systems?

: Each tubulin protein molecule on the cylindrical surface of a microtubule, a fundamental element of the cytoskeleton, acts as a unit cell of a crystal sensor. Electromagnetic sensing enables the 2D surface of microtubule to act as a crystal or a collective electromagnetic signal processing system. We propose a model in which each tubulin dimer acts as the period of a one-dimensional crystal with effective electrical impedance related to its molecular structure. Based on the mathematical crystal theory with one-dimensional translational symmetry, we simulated the electrical transport properties of the signal across the microtubule length and compared it to our single microtubule experimental results. The agreement between theory and experiment suggests that one of the most essential components of any Eukaryotic cell acts as a one-dimensional crystal. we will consider electromagnetic resonance-induced water and ionic ﬂows as ion current waves traveling in a periodic structure. The MT or neuronic crystal period will consist of the tubulin dimer unit (alpha- and beta-tubulin together). The reason for this is because they fold differently, creating two intertwined helicoidal structures. We can model each tubulin monomer as two perpendicular hollow cylinders [23] as depicted in Figure 2. Every monomer has a series resistance RSI (I = β ) and parallel resistance R pi and to keep the model as simple as possible neither capacitance nor inductance effects are taken into account. one described before. The only difference was that it contained only 6 cavities with the same order n n n n n n n . We checked the transmission properties of this shorter MT version and it retained the mean transmission peaks (not shown). The equivalent refractive indices for each monomer were the same as before. A plane wave with an amplitude of 35KV/m excited the MT at an angle of 18 degrees. This value is inferred from the experimental values reported [20] that is 3.5 V/100 microns. In Figure 5 we observe the simulated MT structure.


Introduction
In the past, several theoretical studies were completed by physicists on the vibrational characteristics of bio-living systems, which are currently becoming more attractive research domains in the present. State-of-the-art tools are used to detect vibrational aspects and principle authentication. Traditional theoretical methods are not sufficient to understand the process of biological systems [1,2]. The cooperative interaction between individual components of the developed dynamic system of the biological system changes in a specific fashion [3]. Biological materials show complex vibrational patterns. In other words, solid structures have a crystalline structure and the translational symmetry/spatial periodicity confirms them as an integrated mechanical system, which is well known in solid-state physics or lattice theory. The lattice, known as the base, can be found by the duplication of building elements [4]. The electrons interact with the atomic lattice with atomic potentials. The electronic bandgap is generated when an electronic state does not exist for a certain energy range [5].
Consequently, a wide range of man-made systems is inspired by the design and organization of such biological systems, which include some of their properties. Whereby we have photonic crystals (crystals in which phonons interact with the periodic variation of elasticity and mass of structure elements) and phononic crystals (crystals in which photons are affected by the variation of the dielectric constant of structure) [6,7].
However, there are many examples of photonic crystals [8] in nature, such as the Chrysinaresplendens beetle in which 120 layers produce color, commonly known as onedimensional photonic crystals [9]. More, a two-dimensional photonic crystal film is found in the peacock Pavomuticus, which is composed of more than ten-unit cells in thickness, and that layer has a highly reflective nature and provides a saturated color. The interactive wavelength is determined by the lattice parameters of the layers [10]. Neurofilament and microfilaments are not photonic crystals, although they form the basic features of neurons, they carry out information processing during the cell division process and serve as carriers of vesicles and several components in the biological system. Crystalline structures are found in various forms in nature. Bressloff PC & Cowan JD [11] present a model of pattern formation in the primary visual cortex (V1) of the brain that considers a crystalline structure. In order to understand biological systems more broadly, an interdisciplinary approach based on mathematical and physical sciences can make a significant contribution as biological crystal structures are not limited to photonic domains. The propagation of electromagnetic quanta through a biological lattice leads to a wide range of applications. For example, quantized energy propagation can have applications in wireless electromagnetic fields in medical treatment. The microtubule is the closest candidate related to it.
Microtubules in living cells are organized as networks (where they regularly shrink and grow) and in neurons (where they form highly complex networks with neurons and microfilaments, which generate large-scale communication pathways). To demonstrate the function of cells, microtubules are interconnected through multiple proteins (Figure 1a). The outer and inner diameter voids of the microtubule cylinder are approximately 25 nm and 15 nm respectively and they are found in animals from 200 nm to 25 µm [12]. The length of the dimer unit is 8 nm distributed over a Gaussian distribution with a 2 nm standard deviation and 13 protofilaments are present in the cylindrical surface of the microtubule [13]. The protofilaments are tightly bound from the inside and have weak ties between them, forming 2D sheets that are wrapped with a tube through the nucleation process [14]. Each monomer of the lattice of the microtubule has H12-a small C-terminal helix, which emits after a sequence of amino acids that extend from the microtubule surface.
Consequently, a wide range of man-made systems is inspired by the design and organization of such biological systems, which include some of their properties. Whereby we have photonic crystals (crystals in which phonons interact with the periodic variation of elasticity and mass of structure elements) and phononic crystals (crystals in which photons are affected by the variation of the dielectric constant of structure) [6,7].
However, there are many examples of photonic crystals [8] in nature, such as the Chrysinaresplendens beetle in which 120 layers produce color, commonly known as one-dimensional photonic crystals [9]. More, a two-dimensional photonic crystal film is found in the peacock Pavomuticus, which is composed of more than ten-unit cells in thickness, and that layer has a highly reflective nature and provides a saturated color. The interactive wavelength is determined by the lattice parameters of the layers [10]. Neurofilament and microfilaments are not photonic crystals, although they form the basic features of neurons, they carry out information processing during the cell division process and serve as carriers of vesicles and several components in the biological system. Crystalline structures are found in various forms in nature. Bressloff PC & Cowan JD [11] present a model of pattern formation in the primary visual cortex (V1) of the brain that considers a crystalline structure. In order to understand biological systems more broadly, an interdisciplinary approach based on mathematical and physical sciences can make a significant contribution as biological crystal structures are not limited to photonic domains. The propagation of electromagnetic quanta through a biological lattice leads to a wide range of applications. For example, quantized energy propagation can have applications in wireless electromagnetic fields in medical treatment. The microtubule is the closest candidate related to it.
Microtubules in living cells are organized as networks (where they regularly shrink and grow) and in neurons (where they form highly complex networks with neurons and microfilaments, which generate large-scale communication pathways). To demonstrate the function of cells, microtubules are interconnected through multiple proteins ( Figure  1a). The outer and inner diameter voids of the microtubule cylinder are approximately 25 nm and 15 nm respectively and they are found in animals from 200 nm to 25 μm [12]. The length of the dimer unit is 8 nm distributed over a Gaussian distribution with a 2 nm standard deviation and 13 protofilaments are present in the cylindrical surface of the microtubule [13]. The protofilaments are tightly bound from the inside and have weak ties between them, forming 2D sheets that are wrapped with a tube through the nucleation process [14]. Each monomer of the lattice of the microtubule has H12-a small C-terminal helix, which emits after a sequence of amino acids that extend from the microtubule surface.  This approximate sequence is known as tubulin tail (TT). TTs have a length of about 4-5 nm when they become completely uncontaminated, protrude from the microtubule (MT) surface into the slurry, and form the rough surface of the microtubule filaments. The TTs are key for MT interactions with related proteins. In addition, TTs from a betatubulin monomer could increase or decrease H12 length and consequently the beta-TT length [15,16].This property would be of importance as we model the transport properties in microtubules because TT's flexibility could significantly act on the ionic current flow across the microtubule surface. MTs are predominantly negatively charged, chiefly on their outward surface. Thus, each MT attracts a positive counter of ions near its surface, while negative ions of the cytosol are repelled in such a way that form depletion layers around the microtubule. The thickness of the depletion layer is known as the Bjerrum length l B , which is the length from the surface of the microtubule where the thermal energy is equal to the coulomb energy of the screened surface charges [14].This length greatly determines MT's electric properties. Additionally, the final spatial arrangement of MT protofilaments can be seen as a tilted multilayer structure with a unit cell of 8 nm (Figure 1b). In this way, we can make explicit the method by which wave-like activity emerges from neuronal activity [17][18][19][20][21][22][23][24][25].
Researchers believe that the dissemination of information in neural networks occurs at synaptic junctions, while experimental evidence suggests that microtubules may function as RAM, which demands correction more accurately than at a synaptic junction. It is also noteworthy that a sequence of protein molecules has to undergo a conformational chain in a particular order to open the ion channel. Therefore, vibration is key to the conventional understanding of membranes and ionic-based information processing. Although the terms of neural oscillations [26][27][28] in brain activity are mostly preferred with electrical impulses through membrane molecules, and the oscillatory term of the brain refers to the rhythm of electrical activity [29,30] in the central nervous system, CNS, and this can form the microtubule from the transmission system of the electromagnetic quanta [31,32]. Poznanski and Cacha have considered "protein polarization as a mechanism by way of passive dendrites becoming active by treating the dipolar core as a conductor of current based on the assumption that interactions between individual moving charges can be treated as a homogeneous intracellular medium of constant conductivity and permittivity in space and time.". As can be seen in Figure 1, the MT can be represented as a structure whose permittivity and permeability parameters are changing periodically. Sahu S, et al. proposed that the "water channel and protein molecule together control the emergent properties of the material" [21]. These facts play an important role in MT modeling as electromagnetic quanta processing crystals. A theoretical model is suggested for magnetoencephalography (MEG) and electroencephalography (EEG) that is based on the activity of cortical networks and the transmission line model [20]. The theory of brain waves is defined theoretically based on synaptic action properties that can infer many silent features and some experimental findings using the characteristic of the forward and backward propagating wave.
We speculate that brain waves can be thought of as Bloch waves that may provide important mechanisms for consistent large-scale integration between brain regions [17,18,24]. In this work, we study, experimentally and theoretically, electrical transport properties within tubulin-based microtubules (MT). Using the above MTs properties and symmetry we model them as a periodic crystal with a defect and compare the experimental transmission with our model predicted transmission spectrum.

Crystal Model
We suppose that there is a flow of water and ions from one end of the MTs to the other, produced by the concentration gradient and by inherent electrical fields generated by the charge of the protein spread aligned to every one of the 13 protofilaments. On top of that, the globe-like configuration of tubulin dimers engenders tiny orifices that can be regarded as channels connecting the outer and the inner surfaces of a microtubule.
In these straits, we should expect further escape of water molecules, which would be capable of rectifying the currents along with the depletions of the Bjerrum layer around a microtubule [23]. Moreover, if the water channel is released, the microtubule loses all its unique vibrational properties, and, thus, the water channel controls all optical and electronic properties of the microtubule [20]. Therefore, our work is based on the theory used for electromagnetic radiation in periodic layered media, but, in the present case, we will consider electromagnetic resonance-induced water and ionic flows as ion current waves traveling in a periodic structure. The MT or neuronic crystal period will consist of the tubulin dimer unit (alpha-and beta-tubulin together). The reason for this is because they fold differently, creating two intertwined helicoidal structures. We can model each tubulin monomer as two perpendicular hollow cylinders [23] as depicted in Figure 2. Every monomer has a series resistance RSI(I = β) and parallel resistance R pi and to keep the model as simple as possible neither capacitance nor inductance effects are taken into account.
tronic properties of the microtubule [20]. Therefore, our work is based on the theory used for electromagnetic radiation in periodic layered media, but, in the present case, we will consider electromagnetic resonance-induced water and ionic flows as ion current waves traveling in a periodic structure. The MT or neuronic crystal period will consist of the tubulin dimer unit (alpha-and beta-tubulin together). The reason for this is because they fold differently, creating two intertwined helicoidal structures. We can model each tubulin monomer as two perpendicular hollow cylinders [23] as depicted in Figure 2. Every monomer has a series resistance RSI(I = β) and parallel resistance Rpi and to keep the model as simple as possible neither capacitance nor inductance effects are taken into account.
The parameter f x  0,1   allows theβ-TTs to shrink.  Both resistances are given by the following: The parameter f x ∈ (0, 1] allows the β-TTs to shrink. l B = 6.7 × 10 −10 m, l c = 4.5 × 10 −9 m, R = 12.5nm, and σ = (0.15 ± 0.01)S/m. The monomers' lengths are l α = 4nm and l β = 4nm. The period length p is then l α + l β = 8nm. The tilted angle θ 0 is approximately 18 degrees (Tan −1 (8nm/25nm)). Herein we use a simplified transmission model like the one developed in [31]. Figure 3 shows the transmission line where Z i and Y i are the impedances and admittances of one tubulin dimer unit. Such a transmission line can be simplified like the one shown in Figure 4, where we observe a unit cell consisting of the alpha-MT unit (effective impedance Z 0α ) and beta-MT unit (effective impedance Z 0β ). We can repeat this unit infinitely. The expressions for Z 0α and Z 0β are as follows: We can then identify: fective impedance Z0β). We can repeat this unit infinitely. The expressions for Z 0 and Z 0 are as follows: We can then identify: We can, thus, define the equivalent refractive indices n i for each monomer as  is the surrounding medium refractive index (water). Therefore, the matrices associated We can then identify: Z   R s , Y   1/ R p and Z   R s , Y   1/ R p .
We can, thus, define the equivalent refractive indices n i for each monomer as  is the surrounding medium refractive index (water). Therefore, the matrices associated with each alpha-or beta units are given as follows: We can, thus, define the equivalent refractive indices n i for each monomer as n i = Z Re f /Z 0i , where Z Re f is a reference impedance and sub-index i means the i-th layer (α or β). The transmission line propagation constant is given by Γ i = jζ i = jωn i cos(θ i )l i /v, where ω shows the angular frequency, l i equals the alpha-or beta units' length, v is the propagation velocity of ions. In physiological conditions like an action potential propagation, the velocity of electric fields regulating the vertical transmission of ions lies between 0.1-10 m/s [33]. Note that ions or electrons move with negligible drift velocity (0.001-0.00001 m/s) when current flows even in metals; hence, we should not confuse field flow with the ionic transmission. Herein, we obtained better results when v = 2.65 m/s. To take into account the tilted angle θ 0 , we use Snell's law n 0 sin(θ 0 ) = n i sin(θ i ) and calculates the angle θ i . The parameter n 0 is the surrounding medium refractive index (water). Therefore, the matrices associated with each alpha-or beta units are given as follows: The transfer matrix in a unit-cell is given as follows: Now, if we want to simulate an MT with N-units the final matrix should as follows: We have the following condition involving the ABCD matrix (where A = D in this case due to the symmetry): where K is the propagation constant also known as the Bloch wave number [34] for the periodic structure and p is the period of the structure. From this equation, we can calculate the band structure as Kp = cos −1 (A) (10) When K is real, we have the transmission in the line and when K is pure imaginary, the transmission is forbidden. To obtain the transmission line band diagram, we have to calculate matrix M k (k = 1,2 . . . N) for each frequency ω, multiply all the matrices shown in Equation (8) to obtain M, then use the element A to calculate the propagation constant K.
Finally, we calculate the electrical transport as We used the Fourier transform technique for the synthesis of the optical thin film to obtain the best theoretical fit. We followed the method developed in [35], where further details can be found. The final design consisted of 376 layers where each layer length was 4 nm (1504 nm in total). The final multilayer sequence was as follows: 12 layers with refractive index n Low = 1.07402 then 36 cavities, the order n Low n High n Low n High n Low n d2 n Low n High n Low n High where, n High = 7.3 and n d1 = 3.3, finally 4 layers with series, n Low n d2 n Low n d2 where n d2 = 2 is added. The lowest refractive index values (n Low ) represent α-tubulin monomers on the contrary n High , n d1 , and n d2 constitute β-tubulin monomers. The whole structure was embedded in the water n 0 = 1.33.

Finite Element Analysis
It was possible to perform full-wave finite element simulations (FEMLAB 3.1, COM-SOL AB, Burlington, MA, USA) on a shorter MT version (Length of 304 nm). The monomers' lengths were 4 nm each. The final layer sequence was similar to the one described before. The only difference was that it contained only 6 cavities with the same order n Low n High n Low n High n Low n D n Low n High n Low n High . We checked the transmission properties of this shorter MT version and it retained the mean transmission peaks (not shown). The equivalent refractive indices for each monomer were the same as before. A plane wave with an amplitude of 35KV/m excited the MT at an angle of 18 degrees. This value is inferred from the experimental values reported [20] that is 3.5 V/100 microns. In Figure 5 we observe the simulated MT structure.  [20] that is 3.5 V/100 microns. In Figure 5 we observe the simulated MT structure. The mesh was created with isosceles triangles with a base averaged length of 0.5 nm and height of 0.25 nm and 10 −6 was the numerical tolerance. The simulation results are the normal component of the electric field and in-plane magnetic field components. We used TE polarization configuration to solve the Helmholtz equation given as .
(12) The mesh was created with isosceles triangles with a base averaged length of 0.5 nm and height of 0.25 nm and 10 −6 was the numerical tolerance. The simulation results are the normal component of the electric field and in-plane magnetic field components. We used TE polarization configuration to solve the Helmholtz equation given as where ω-the angular frequency; n-the refractive index and the propagation velocity of ions, which has a value of 2.65 m/s.

The electric field
→ E z is oriented in the z direction (perpendicular to monomers' layers).
All external boundaries fulfilled a matched boundary condition given by where k water = ωn 0 /v is the water wave number and → n is a unitary vector normal to external boundaries. This boundary condition permits the inclusion of water as a surrounding medium (n 0 ).
The transverse field components act as current and voltage on a transmission line model and are known as transverse equivalent networks [36]. The correspondence relationship can be understood as (i = TM or TE): The unit vectors are chosen thatn z ×n T = ±n P , wheren P is a unit vector in the propagation direction,n T is a unit vector perpendicular to propagation direction andn z is a unit vector perpendicular to the propagation plane. The function Ψ t (y, z) is known as the transverse wave function. For instance, in the case of plane waves traveling in the ±x direction, it is written as Ψ t (y, z) = e −i(k y y+k z z) . Therefore, magnetic field modes should provide insights into the current distributions along the MTs. The simulated magnetic modes are given by the normalized summation of the real part of both magnetic field components parallel to multilayers. That is H x + H y / H x + H y .

Discussion
We detected experimentally the microtubule's electrical transport properties and following various protocols we measured the transmission spectra across the single microtubule. For over a decade we have published a series of experimental measurements where transmissions of various carriers have been measured using quantum tunneling [21,22,29,31]. Figure 6a shows MT's electrical transport bands for a spectral bandwidth spanning from 10 KHz to 250 MHz. The MT length was 1500 nm. We simulated the power and transport band of MT as shown in Figure 6a from the electrical parameters obtained from the model created. Such a model consists of 376 layers, where each layer has the length l α = l β = 4nm (Figure 6b). The model fits very well the experimental transmission spectra. Only peak values of the order of 20% at normalized frequencies bigger than 1 are not present in the theoretical model. Nevertheless, we do not know if the experimental peaks in that region are meaningful or just noise. The electrical impedance Z 0α has only one value of 2.68885GΩ calculated with the help of Equation (2). The other Z 0β electrical impedances were calculated from the effective refractive indices using the ratio between them. That is, since we know the value of the ratios n High /n Low , n d1 /n Low , and n d2 /n Low and, by definition, they are equal to Z 0α /Z 0β High , Z 0α /Z 0β d 1 , and Z 0α /Z 0β d2 , respectively. Thus, Z 0β High = 0.395234GΩ, Z 0β d1 = 0.87498GΩ, and Z 0β d2 = 1.43725GΩ. These values are obtained when the parameter f x equals 0.02063, 0.105, and 0.285 respectively.
The built model has two open parameters as the length of beta-TT and the propagation velocity of ions. We found the finest value of beta-TT length are 0.09 nm, 0.47 nm, and 1.28 nm and 2.65 m/s for the propagation velocity. Although, the model is simple and could speculate the whole electrical band structure (see Figure 6). These results highlight the fact that molecular configurational outcomes are significant in microtubules because these outcomes could modulate energy transportation. Indeed, simply changing β-TTs' lengths drastically modifies the electrical impedances to create the electrical band structure.
Other results that support the idea of standing wave formation are presented in [21,37] where scanning tunneling images (STM) of a single microtubule were presented. After an AC signal is pumped over the microtubule at a particular frequency, the STM images of the induced current were recorded. From the images, one can clearly observe different ion current configurations that might be thought of as ion current modes similar to electromagnetic modes in a photonic crystal. Along the same lines, the appearance of these ion current modes only happens when water is presented. That is if the water is out of the MT core, the STM images cannot be taken, and the MT disintegrates as the researchers scan it. This also reinforces the idea that water is necessary to create ion waves that later could form standing waves within the MT. From the finite element analysis (Figure 7) we can observe that simulation describes fairly well the experimental current distribution for different applied AC voltage frequencies up to 113 MHz (Figure 7a-c). Nonetheless, for the frequency of 185 MHz the experimental pattern is replicated only in two layers (shown by asterisks). obtained from the model created. Such a model consists of 376 layers, where eac has the length l   l   4nm (Figure 6b). The model fits very well the exper transmission spectra. Only peak values of the order of 20% at normalized freq bigger than 1 are not present in the theoretical model. Nevertheless, we do not k the experimental peaks in that region are meaningful or just noise. The electri pedance Z 0 has only one value of 2.68885G calculated with the help of Eq (2). The other Z 0 electrical impedances were calculated from the effective re indices using the ratio between them. That is, since we know the value of th n High / n Low , n d1 / n Low , and n d 2 / n Low and, by definition, they are eq Z 0 / Z 0 High , Z 0 / Z 0 d 1 , and Z 0 / Z 0 d 2 ,respectively. Thus, Z 0 High  0.39523 Z 0 d1  0.87498G , and Z 0 d 2  1.43725G . These values are obtained wh parameter f x equals 0.02063, 0.105 , and 0.285 respectively. Microtubule's STM images are taken while pumping materials using a GHz antenna. We have described the experimental details in earlier works [21,22,29,37].

Conclusions
From this analysis, we can conclude that MT is a one-dimensional structure with linear translational symmetry, but it would be interesting to use a helicoidal symmetry to make our results more realistic. Another important restriction of our model is that capacitance and inductance effects were not taken into account, which is equivalent to ne-

Conclusions
From this analysis, we can conclude that MT is a one-dimensional structure with linear translational symmetry, but it would be interesting to use a helicoidal symmetry to make our results more realistic. Another important restriction of our model is that capacitance and inductance effects were not taken into account, which is equivalent to neglect energy emission or absorption by MT's dimers. These effects can modulate electrical transmission peaks' amplitudes and bandwidths. In future work, the model can be improved by taking into account both aspects.

Microtubules Preparation
In the experiment, porcine brain neuron-derived microtubules were used. The microtubules were reorganized in our laboratory and we purchased tubulin proteins with all associated tools from Cytoskeleton Inc. (Denver, CO, USA). The −80 • C temperature freezer was used to preserve tubulin protein. To polymerize tubulin into a microtubule of 6.5 µm length, a 160 µL microtubule cushion buffer is mixed with 10 µL of 100 Mm GTP solution and 830 µL general tubulin buffer. Such a solution is placed in an ice bath for 10 min. 200 µL solution of the mixture was added to 1mg general tubulin solution. We put it in an ice bath again for the same period. After that, such a solution is kept at 35-37 • C for 40 min. To achieve stabilization of the microtubule, paclitaxel is decomposed into anhydrous DMSO and mixed with the solution. The entire mixture is again kept at 37 • C for 10 min. The length of the microtubule varied 4-20 µm. The composition of the microtubule cushion buffer and general tubulin buffer solutions are given as 1 mM MgCl 2 , 60% v/v glycerol, 1 mm EGTA, 80 mm PIPES pH 6.8 and 2 mM MgCl 2 , 80 mM PIPES pH 7, 1 mM EGTA, respectively.
To prepare the film, the prepared solution is diluted several times in a paclitaxel DMSO solution and microtubule cushion buffer, after which, it is dropped onto a Si substrate, and interconnected electrodes are placed in an ice bath (−20 • C) for the entire night. For the proper alignment of the microtubule, an electric field is applied to the Si substrate. A filter paper is used to remove the accesses solution of microtubules on the substrate which is the appropriate method of study of AFM and STM. The Si substrate is then dissolved in tubulin solution and N2 is dried in the same manner [21].
The reformation of α and β tubulin into αβ heterodimer has the dimensions 46 × 80 × 65 Å as confirmed via UHV-STM and AFM [21,22,29,37]. The solution containing the microtubule is dropped on the pre-grown electrodes and then the chip is cleaned from large ion deposits by dielectrophoresis. The dielectric property of the microtubule was measured with the electro-orientation method [19]. The analysis of the dielectric properties of the microtubules reports a large charge density because the reported value is much greater than in an organic semiconductor, but smaller than in a metal, the microtubule cannot be classified as either.

Protein Preparation
The −80 • C temperature is used in protein preservation. Protein orientation is performed when the protein solution is dropped into a small gap between electrodes under a 4V DC electric field. The nanomolar solution is precipitated on a highly oriented pyrolytic graphite, HOPG substrate. The resonant AC frequencies are used to scan the protein and its formation, otherwise, they go away. The solution of tubulin solution is dropped at the junction of the 4 probe Au electrodes. A constant value of a 0.5 V electric field is reliably kept on each pair of electrodes.

Electrical Transport Measurements
MT transmittance was measured with a four-probe circuit [21], there are some special frequencies where the microtubule conducts the AC signal with almost no AC resistance (<<<0.04 MΩ) at those frequencies.
The details can be found in [21]. Both spectra were measured between 10 kHz to 20 GHz. It has been suggested that the microtubule resonance frequencies significantly change the lattice geometry of the microtubule. Every single resonance frequency is associated with particular lattice geometry. Using a specially designed scanning tunneling microscope, the cylindrical surface has been imaged live. Several fundamentally different experiments support each other in the finding that the microtubule is an insulator; however, several interesting properties of microtubule emerge as soon as it is triggered by resonance frequencies. Using several versions of the tubulin proteins, it has been demonstrated that electromagnetic resonance has a profound effect, even in their self-assembly.