Next Article in Journal
High Entropy Oxide Phases with Perovskite Structure
Next Article in Special Issue
The Self-Adaptation Ability of Zinc Oxide Nanoparticles Enables Reliable Cancer Treatments
Previous Article in Journal
Influences of Nanoparticles Characteristics on the Cellular Responses: The Example of Iron Oxide and Macrophages
Previous Article in Special Issue
Eu-Doped Citrate-Coated Carbonated Apatite Luminescent Nanoprobes for Drug Delivery
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Nanomaterials
Volume 10
Issue 2
10.3390/nano10020265
nanomaterials-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigation of the Electrical Properties of Microtubule Ensembles under Cell-Like Conditions

by
Aarat P. Kalra
1,
Sahil D. Patel
2,
Asadullah F. Bhuiyan
2,
Jordane Preto
1,
Kyle G. Scheuer
2,
Usman Mohammed
3,
John D. Lewis
4,
Vahid Rezania
3,
Karthik Shankar
2,* and
Jack A. Tuszynski
1,4
1
Department of Physics, University of Alberta, 11335 Saskatchewan Dr NW, Edmonton, AB T6G 2M9, Canada
2
Department of Electrical and Computer Engineering, University of Alberta, 9107–116 St, Edmonton, AB T6G 2V4, Canada
3
Department of Physical Sciences, MacEwan University, Edmonton, AB T5J 4S2, Canada
4
Department of Oncology, University of Alberta, Edmonton, AB T6G 1Z2, Canada
*
Author to whom correspondence should be addressed.
Nanomaterials 2020, 10(2), 265; https://doi.org/10.3390/nano10020265
Submission received: 4 December 2019 / Revised: 28 January 2020 / Accepted: 29 January 2020 / Published: 5 February 2020
(This article belongs to the Special Issue Application of Nanoscale Materials for Cancer Diagnostic and Therapy)
Browse Figures
Versions Notes

Abstract

:
Microtubules are hollow cylindrical polymers composed of the highly negatively-charged (~23e), high dipole moment (1750 D) protein α, β- tubulin. While the roles of microtubules in chromosomal segregation, macromolecular transport, and cell migration are relatively well-understood, studies on the electrical properties of microtubules have only recently gained strong interest. Here, we show that while microtubules at physiological concentrations increase solution capacitance, free tubulin has no appreciable effect. Further, we observed a decrease in electrical resistance of solution, with charge transport peaking between 20–60 Hz in the presence of microtubules, consistent with recent findings that microtubules exhibit electric oscillations at such low frequencies. We were able to quantify the capacitance and resistance of the microtubules (MT) network at physiological tubulin concentrations to be 1.27 × 10−5 F and 9.74 × 104 Ω. Our results show that in addition to macromolecular transport, microtubules also act as charge storage devices through counterionic condensation across a broad frequency spectrum. We conclude with a hypothesis of an electrically tunable cytoskeleton where the dielectric properties of tubulin are polymerisation-state dependent.

1. Introduction

Microtubules (MTs) are cylindrical polymers composed of the heterodimers of protein α, β- tubulin that play a variety of well-recognised intracellular roles, such as maintaining the shape and rigidity of the cell, aiding in positioning and stabilisation of the mitotic spindle for allowing chromosomal segregation, acting as ‘rails’ for macromolecular transport and forming cilia and flagella for cell movement. Since the tubulin dimer possesses a high negative electric charge of ~23e and a large intrinsic high dipole moment of approximately 1750 D [1,2], MTs have been implicated in electrically-mediated biological roles [3,4,5,6]. They have been modelled as nanowires capable of enhancing ionic transport [7,8], and simulated to receive and attenuate electrical oscillations [4,9,10,11]. In solution, MTs have been shown to align with applied electric fields [2,12,13,14,15,16]. Recently, MTs have also been modelled as the primary cellular targets for low-intensity (1–2 V), intermediate-frequency (100–300 kHz) electric fields termed TTFields (tumour-treating electric fields) that inhibit cancer cell proliferation, in particular glioma [17,18,19]. Indeed, MTs have been reported to decrease buffer solution resistance [12,13], leading to a conductance peak at frequencies close to the TTField regime [20]. While these studies show that MTs are highly sensitive to external electric fields, answers to the questions ‘How do MTs effect a solution’s capacitance?’ and ‘What is the capacitance of a single MT?’ are still elusive and crucial to the determination of the dielectric properties of living cells. The tubulin concentration in mammalian cells varies in the micromolar range (~10–25 μM) [21,22]. In vitro, polymerizing tubulin at such high concentrations can lead to the formation of entangled networks, confounding quantification of the individual MT response to electric fields. Electro-rotation, di-electrophoresis and impedance spectroscopy are thus performed using low concentrations of tubulin, in the nanomolar regime, to enable robust observation of individual MTs.
MT formation and stability are known to be optimal in buffers with ionic strength between 80 and 100 mM [23,24]. A background of BRB80 (which consists of 80 mM PIPES, 2 mM MgCl2 and 0.5 mM EGTA, containing ~46 mM PIPES 2–, ~36 mM PIPES–, ~68 mM Cl–, ~160 mM K+, and ~2 mM Mg2+ [2]), is thus used to study the dynamics and mechanical properties of MTs. To study their electrical properties however, the usage of such high ionic-strength solutions has historically been problematic because any dielectric attenuation caused by MTs is overwhelmed by the noise and high conductivity from the background. In the low-frequency regime (1–100 kHz), two approaches have thus far been used to estimate the dielectric properties of MTs and tubulin. One is to electrically observe low concentrations of MTs (tubulin concentration in the nanomolar regime) in the presence of low ionic strengths [12,13,20,25,26]. Such studies overlook the intrinsic ionic concentration of mammalian cytosol, which varies between 200 to 500 mM depending on the cell type [27,28]. Another approach to electrically interrogate MTs is to dry them: the conductivity of the buffer is nullified by evaporation, leaving polymeric tubulin behind [29,30]. In a physiological situation however, MTs are solvated by the highly conductive and noisy cytosol.
Here, we report on our efforts overcome the barrier posed by a high ionic strength by performing electrochemical impedance spectroscopy (EIS) on cellular concentrations of tubulin. We have been able to successfully observe differences in impedance using a background of BRB80 itself. Surprisingly, we find that MTs increase the solution capacitance of BRB80 whereas free tubulin does not, implicating a difference in electrical properties based only on the morphology of this protein solute. We also report a ‘reversal’ in the resistive behaviour of MTs compared to BRB80, with a reduction in solution resistance peaking in the 20–60 Hz region, a finding consistent with recent reports showing that polymerised tubulin quasi-resonantly responds to electric oscillations at ~39 Hz [9,10]. Using an equivalent circuit model for MTs, we experimentally determine the capacitance and resistance of the MT network to be 1.27 × 10−5 F and 9.74 × 104 Ω respectively, at physiological concentrations of tubulin. Our values indicate that the polymerisation of tubulin into MTs alters the spatial and temporal charge distribution, altering the electrical properties through charge storage in the cell.

2. Materials and Methods

2.1. Tubulin Reconstitution

Fluorescently labelled tubulin solution was prepared using previously published protocols [20]. Notably, tubulin stock powders that were devoid of MAPs were purchased. Lyophilised unlabelled tubulin powder (Cytoskeleton Inc, Denver, CO, USA; T240) was reconstituted in BRB80 supplemented with 1 mM GTP (guanosine triphosphate; Cytoskeleton Inc, Denver, CO, USA; BST06) and mixed with tubulin labelled with a rhodamine-based ester (Cytoskeleton Inc, Denver, CO, USA; TL590m) in a final labelling ratio of 1:15. Aliquots were snap-frozen and stored at −80 °C.

2.2. MT Polymerisation and Stabilisation

MT polymerisation was performed by incubating 45.45 μM tubulin aliquots in a 37 °C water bath for 30 minutes. BRB80 solution was heated alongside tubulin during the first 15 minutes of polymerisation. Subsequently, BRB80 was incubated at room temperature, and paclitaxel solution (Cytoskeleton Inc, Denver, CO, USA; TXD01; 2 mM stock) was thawed at room temperature alongside it. After 30 minutes of tubulin polymerisation brought to completion, 100 μL of BRB80 was added to 5 μL of 2 mM paclitaxel (BRB80T). For preparing 0.222, 2.225 and 22.225 μM MTs the above solution was added in different volumes to polymerised tubulin, as shown in Table 1. For preparing the BRB80T background for impedance measurements, 45 μL of this solution was added to 45 μL of BRB80.
For tubulin stabilisation, 2 μL of 5 mM colchicine stock solution (Sigma-Aldrich, St. Louis, MO, USA; C9754; 5 mM in DMSO) was added to 100 μL BRB80 (BRB80C). Subsequently, a similar solution to BRB80T was prepared using colchicine. For preparing 0.222, 2.225 and 22.225 μM free tubulin solutions, the above solution was added in different volumes to free tubulin solutions, as shown in Table 1. For preparing BRB80C, 45 μL of this solution was added to 45 μL of BRB80.

2.3. Fluorescence Imaging of MTs

Imaging was performed on a Zeiss Examiner.Z1 microscope using a Hamamatsu (Hamamatsu City, Japan) EMCCD C9100 camera, a Zeiss (Oberkochen, Germany) plan-Apochromat 1.4 NA 63x lens. After pipetting MT solution (2–5 μL) onto a glass slide (VWR 48382-173) a coverslip (VWR 48393-070) was placed on the solution, allowing it to spread. The microscope used an EXFO X-Cite 120 fluorescence source and excitation and emission filters of 535 nm and 610 nm, respectively. Exposure times between 50 ms and 300 ms were used for imaging to validate the presence of MTs.

2.4. Electrode Design and Device Construction

Each ‘plate’ in the parallel-plate contact device was formed by FTO (fluorine-doped tin oxide)-coated glass slides (Sigma Aldrich, St. Louis, MO, USA; 735140). The slides were cleaved to dimensions of 1.5 mm × 10 mm × 50 mm for the upper contact and 1.5 mm × 27 mm × 50 mm for the lower contact. The cleaving dimensions were set using 3D printed devices that were placed as holders (The Shack, University of Alberta; Figure S1 in Supplementary Materials). The slides were ultrasonicated and subjected to reactive ion etching (RIE) using a 5-minute exposure to oxygen plasma (Oxford Instruments, Abingdon, UK; NGP80) to remove surface particulate matter. A 70-μm thick double-sided tape was used as a spacer, which formed a chamber of dimensions 3 mm × 1.25 cm × 70 μm. The top electrode was placed using a separate 3D-printed holder device (Figure S1). Once the device was constructed using the above protocol, solution was perfused into the chamber using a pipette and a filter paper for suction. We used flat copper electrode clips in a three-electrode configuration to connect to our capacitor device. The counter electrode was connected to the lower electrode, and the working and reference electrodes were connected to the top electrode of our device.

2.5. Impedance Measurements

Experiments were conducted using Electrochemical Impedance Spectroscopy (EIS) on a Zahner Zennium impedance analyser. The parallel-plate contact device was placed into the 3D-printed holder for stabilisation (Figure S1). The contacts from the machine were connected to the parallel-plate device using flat-faced copper alligator clips. A three-electrode configuration was used: The counter electrode was attached to the lower contact of the parallel-plate device, whereas the working electrode was attached to the upper contact with the reference electrode orthogonally clipped onto the clip of the working electrode. Within the Thales Z3.04 environment, the potentiostat mode was ON; the stabilisation delay was set to 1 s, the rest potential drift tolerance was set to 250 μV, Vrms was set to 5 mV. Solutions were perfused into the experimental chamber using a micropipette tip at one opening, and a filter paper at the other opening for suction, similar to protocols used for Total Internal Reflection Fluorescence (TIRF) microscopy [31]. The frequency range of the EIS measurement was set from 4 MHz to 1 Hz and data were subsequently collected.

2.6. Data Analysis

MT and tubulin samples were analysed using data from five to seven days of experiments. Each day consisted of three to seven solutions for each concentration being tested, with one frequency sweep per solution. Readings of each sweep were saved as a csv file, and next sample was loaded by solution exchange method. Water was run as the first solution for each day of experiments. BRB80T was run prior to MT solutions, and BRB80C were run prior to the free tubulin containing solutions. MT- and free tubulin-containing solutions were run on separate days, in increasing order of concentration. MATLAB (The Mathworks; Natick, MA, USA) scripts were used for data analysis. Fitting to the real and imaginary components of impedance was performed using the function lsqnonlin. Initial guess values for the MT network resistance and capacitance were 105 F and 10−5 Ω, respectively, based on visual inspection of raw data. The initial guess values for the nominal series resistor, RH, were set at 1.78, 0.6 and 0.4 Ω with tubulin concentrations of 0.222, 2.222 and 22.222 μM, respectively. The 95% confidence intervals were determined using the function nlparci. Error propagation was performed assuming no relationship between various days of data collection.

3. Results

3.1. Validation of Parallel-Plate Contact Device to Measure Dielectric Properties of Physiologically Relevant Ionic Solutions

To determine the differences in the dielectric properties of solution caused by the presence of MTs, we first aimed to create an electrode geometry that would be experimentally robust and easily modelled. We fabricated an FTO-coated parallel-plate contact device (Figure 1a), which allowed EIS using a solution-exchange method.
We started by performing EIS on electrolytes found in the cytosol and observed that the imaginary component of impedance became less negative as a function of applied input frequency (Figure 1b). The total impedance of our system was given by:
Z = rc + rs/(1 + (rs ωC)2) + j(ωLc – (rs2 ωC)/(1 + (rs ωC)2)),
Here, Z is the impedance, ω is the angular frequency (given by 2πf where f is the input voltage frequency), C is the system capacitance, Lc is the cable inductance, rs and rc are the solution and cable resistances respectively. We also observed a decrease in the real component of impedance as a function of input frequency (Figure 1c). Such a trend is expected from Warburg impedance [32,33] and is in accordance with the equation:
Zcomplex = (Aω)/√ω + (Aω)/(j√ω),
Here, Zcomplex is the complex impedance and Aω is the Warburg coefficient. Our circuit simplifies to the equation below if we ignore the effect of cable inductance ωLc, at frequencies below 105 Hz:
Z = rc – j/ωC,
Our results using various electrolytes emulated previous data [34,35,36] and validated the experimental setup for further analysis.

3.2. The Effect of Microtubule Networks on Solution Capacitance at Physiologically Relevant Conditions

We reconstituted and polymerised fluorescent tubulin from a stock of 45.45 μM tubulin solution (Materials and Methods). MTs were stabilised using 50 μM paclitaxel [37,38] and imaged using an epi-fluorescence microscope. On diluting MT concentration across three orders of magnitude (0.222, 2.222 and 22.225 μM tubulin), we observed that while individual MTs at low concentrations were separated by large distances, those at cellular concentrations formed enmeshed networks reported previously (Figure 2a–c) [39]. Such interconnected MT networks are utilised by molecular motors for long-range macromolecular transport [40,41]. Here, their presence demonstrated successful MT polymerisation for electrical characterisation.
We performed EIS on BRB80, BRB80T (BRB80 supplemented with 50 μM paclitaxel; background for all MT-containing solutions), and MT-containing solutions in increasing order of concentration (Figure 3a,b). We subtracted impedance values obtained for BRB80T alone from those in MT-containing solutions to determine the MT contribution to impedance. Our results showed that with an increasing MT concentration, the value of imaginary impedance became more negative, resulting in positive impedance differences (Figure 3c–f). This effect was greatest at the cell-like 22.225 µM tubulin concentration, with impedance differences lowering in magnitude with increasing input frequency (Figure 5a). Experiments with unpolymerised tubulin at the same concentrations were performed using the identical procedure, but using BRB80C (BRB80 was supplemented with 50 μM colchicine) as a background, to prevent MT nucleation [42,43]. Results with unpolymerised tubulin did not show an appreciable deviation from zero at any concentration (Figure S5a). The above results suggest that polymerisation of tubulin into MTs alters their ensemble electrical properties, increasing the solution’s capacitance on forming MTs and their networks. An increase in the solution’s capacitance because of MTs has previously been modelled, [7,44,45] indicating an increase in charge storage as free tubulin polymerises.

3.3. The Effect of Microtubule Networks on Solution Resistance at Physiologically Relevant Conditions

Next, we investigated the differences between MTs and tubulin in the real component of impedance (solution resistance). Previous studies using nanomolar tubulin concentrations and low ionic strengths (1–12 mM) have indicated that MTs enhance charge-transport in solutions [13,20,46]. To evaluate if this observation held true at physiologically relevant tubulin concentrations and at higher ionic strengths, we also analysed the real component of impedance. Addition of both MTs and tubulin generally led to an increase in solution resistance (Figure 4a–f), with MTs having a higher resistance at low frequencies (1–20 Hz) compared to unpolymerised tubulin. Unexpectedly, a ‘reversal’ of this behaviour was observed at higher frequencies as MTs began to lower the solution resistance compared to tubulin (Figure 5b). The reversal took place gradually between 10 and 300 Hz (Figure 6a–d), with a peak between 20 and 60 Hz (Figure 6e). Interestingly, within this range, we also found that the addition of MTs lowered solution resistance compared to background buffer BRB80T.
Such a reversal in resistance between microtubules and tubulin has not been reported before. Because the extent of this reversal decreased with decreasing concentration, this result also displays the utility of our ‘cell-like’ approach. Our results are consistent with predictions of an increase in solution conductance at ~39 Hz [9,10], which have been hypothesised to arise from oscillatory ionic movement across the MT lattice through nanopores formed between adjacent tubulin dimers (Figure 9a).
It is worth noting that this region falls within the gamma frequency regime (20–60 Hz), implicating such quasi-resonant phenomena as a possible explanation for the source of low frequency intraneuronal electrical oscillations. No such reversal was observed for the corresponding frequency range in the imaginary impedance values (Figure S4).

3.4. The Microtubule Network as an RC Circuit in Parallel

Our next aim is to quantify the resistance and capacitance of the microtubule network. The slope of approximately negative unity on the impedance difference curve suggested that the microtubule network resulted in the addition of a capacitive element to the solution. We examined several combinations but a parallel RC (resistor-capacitor) circuit to represent the entire MT network provided the best fit to observed curves.
We modelled the impedance caused by external circuit elements and BRB80T as Zo and Zs respectively, as shown in Figure 7. The net impedance of the background BRB80T was thus given by:
Zbuffer = Z0 + Zs,
Denoting the impedance, resistance and capacitance of the entire MT network by ZMT, RMT and CMT respectively, the impedance for the circuit with MTs is given by:
ZMT+buffer = Z0 + Zs + RH + ZMT
where,
1/ZMT = 1/RMT + jωCMT
Additionally, the impedance differences between solutions with and without MTs are given by:
ΔZ = ZMT+buffer – Zbuffer = RH + ZMT,
where
ZMT=RMT/(1 + (ωCMTRMT)2) – j(ωCMTRMT2)/(1 + (ωCMTRMT)2),
We subsequently fit experimental impedance difference curves shown in Figure 3 and Figure 4 to real and absolute value of imaginary parts of ΔZ using RH, RMT and CMT as our fit parameters. Here, RH is a resistance ascribed to the nominal fraction of unpolymerised tubulin present in MT containing solutions. The fitted curves are displayed in Figure 8 and the optimal fit parameters are listed in Table 2 (see Materials and Methods for details).

4. Discussion

Our measurements using a parallel plate contact device reveal interesting electrical properties of MTs at physiological concentrations. Unlike studies exposing MT-containing solutions to non-uniform electric fields [12,13,14,20], our device allowed robust quantification of electrical impedance in the presence of spatially uniform electric fields. Our results show that the addition of the MT network mimics a parallel RC element placed in series with the high-ionic strength solution, with a nonlinear dependence on MT number. Unpolymerized tubulin did not alter the capacitance significantly, indicating changes in electrical properties of tubulin as it polymerizes.

4.1. The Physical Underpinnings of An Increased Capacitance

An increase in capacitance arises from dense counterion condensation on the MT surface. This has been extensively predicted and simulated to arise from a variety of sources [7,8,44,47,48]. First, the negative charge of the tubulin dimer attracts counterions in solution, leading to the presence of a double layer and depletion region outside the microtubule surface [7,8,46,48]. The charge distribution in the MT protein wall is also highly non-uniform, with the outer surface containing approximately four times the charge compared to the inner surface [47] (Figure 9c). This asymmetry between the inner and outer electrostatic potentials serves to enhance capacitance and is responsible for the abnormally large dipole moment of the tubulin dimer [1]. The asymmetry also manifests through C-terminal ‘tails’ composed of 10–12 amino-acids, that can extend 4–5 nm outwards from each tubulin monomer. These slender C-termini tails are highly negative, containing about 50% of the charge of the tubulin dimer [49]. As they stretch outwards into the solution in a pH and ionic strength-dependant manner, they increase the effective area of the tubulin dimer and significantly contribute to the overall MT capacitance [7,8].
Coherent oscillations of these C-terminal tails are modelled to generate solitonic pulses of mobile charge along the outer surface of an MT, creating ionic currents along its length [7,44,50]. Ions from the bulk solution are also modelled to be pumped into the hollow MT lumen through nanopores in its wall, resulting in charge accumulation inside the cylindrical MT over time [45]. A recent study using molecular dynamics simulations showed that the permeability of the MT lumen was significantly higher for Na+ and K+ as opposed to Ca2+, allowing for free movement of selective ions into the MT lumen across its porous surface [47]. To the best of our knowledge, our findings are the first to experimentally quantify this resistance encountered by charge flow across the MT cross section. These results implicate not only ionic movement along the microtubule axis, but also across and inside it, enhancing the modelled roles of MTs as complex subcellular nanowires.
Manning’s theory of polyelectrolyte solutions predicts the conditions for ionic condensation on charged polymer surfaces provided a sufficiently high linear charge density is present on these surfaces creating an ionic concentration depletion area surrounding them [51]. The sum total of the charges on polymer surfaces and the associated counterions decreases to values dependent on the valence of the counterions and the Bjerrum length, which is the distance from the polymer surface at which the Coulomb energy of the screened surface charges equals the thermal energy. The double layer of surface charges and counterions separated by the Bjerrum length can be viewed as having capacitor-like properties. Although the Manning theory was originally developed for such polyelectrolytes as DNA, it was also applied to actin filaments [52] and MTs [53]. For actin filaments, its application explained the observed lossless transmission of electric pulses along the filament lengths. In the case of MTs, it provided a plausible explanation of unusual amplification of injected electrical signals that propagated along these nanowires. The calculated Bjerrum length for MTs was found to be approximately 6.7 × 10−10 m [8,50]. Both actin filaments and MTs have been represented in these models by cable equations with effective real and imaginary impedance due to the viscosity of the solution-resisting ionic flows and the capacitive properties of the ionic double layers around the filaments, respectively [52,53]. The capacitance for a single ring of an MT including C-termini was calculated to be approximately 1.3 × 10−15 F [8]. When extended to 20 µm, (representative of the length of a single MT for our measurements), the predicted value would be C = 3 × 10−12 F, although an experimental confirmation of this prediction is not directly available through our measurements or in any previous work. We note the relatively weak dependence of network capacitance on MT concentration, and assign it to the random spatial locations and directional orientations of MTs in our solution. Indeed, the conductivity of randomly distributed RC networks has been shown to scale weakly with the number of elements in the network [54]. Additionally, qualitative similarities can be found in the models of random resistor and capacitor networks with a frequency-dependent crossover for both conductance and impedance in these networks due to percolation-type conduction [55]. We intend to develop a quantitative model for our experimental observations in a subsequent publication.

4.2. Implications for the Cell

Our work, which utilizes cell-like tubulin and ionic concentrations for the first time, indicates a cellular role for microtubules as wires that store charge. Neuronal environments where MTs are spontaneously nucleated from free tubulin, such as growth cones, experience large capacitance changes over short bursts of time. This ability significantly impacts the action potentials that are known to depend strongly on the local charge distributions [56]. Additionally, ionic movement across the MT wall enhances their roles as attenuators of local cation distributions. In nonneuronal environments, transient ionic currents around a MT during mitosis could impact MT dynamics and potentially influence the chromosome segregation. Specifically, Ca2+ ion storage/flow about an MT triggers its depolymerisation, whereas waves of Mg2+ or lowering in the local pH (increasing H+) leads to MT stabilisation [57,58]. The attraction of Zn2+ or Mn2+ ions in the vicinity leads to the formation of two-dimensional tubulin polymers [59,60]. Properties of the cytoplasm such as polarisability and relative permittivity get severely attenuated because of the presence of MTs in the vicinity. Because of the polymerisation state of tubulin-altering solution capacitance, our findings implicate a temporal evolution of capacitance and ionic flows as the ratio of MTs to free unpolymerised tubulin changes [61,62,63]. MT lattice defects, which occur when a tubulin dimer is missing in an MT wall [64,65], cause a large ionic flux to develop at the defect site. Spatiotemporal charge distribution shifts are also critical at the MT end, where fluxes form because of sudden changes due to the polymerisation/depolymerisation of the MT. Free/polymerised tubulin hence regulates local and global electrical properties, creating spatially dynamic gradients of charge storage and flux. We envision a cytoskeleton that, in addition to transporting macromolecules, stores and transports ionic signals and electrical information across the cytoplasm (Figure 9a,b).
Our findings can be coupled with a vast array of bio-nanodevices that utilises MTs and MAPs (microtubule-associated proteins) for construction of bio-nanotransporters and bio-actuators [66,67,68,69,70]. Under specific conditions, MAP-MT systems are capable of repositioning macromolecules [71,72], directionally transporting microtubules [15,73] and even drive their movement within zero-mode waveguides [74] and inorganic nanotubes [75]. Storage of electrical charge and its transport along MTs can be coupled to such cutting-edge mechanical MAP-based devices to develop a wide range of nano-actuators and nano-sensors.
When compared to cells, the rates of MT nucleation and polymerisation are significantly lower in BRB80. This difference can be attributed to the absence of MAPs and macromolecular crowding [76,77]. Mammalian cells contain high concentrations of K+ ions (140–300 mM) [27,28], which, in addition to MAPs and molecular crowding agents, will be included in a future study to attain physiological equivalence. We also note that the effect of PTMs (post-translational modifications) on the electrical properties of microtubules has not yet been explored. PTMs involve the addition of residues such as phosphate and glutamate that locally influence counterionic condensation around the outer microtubule surface.
We are in the process of performing DC (direct-current) measurements, determine the contribution of MTs to impedance relaxation time and evaluate the voltage dependence of capacitance on MT-containing solutions. Interestingly, this aspect has been discussed previously: the inductance of a single protofilament is calculated to be <1 fH [8]. Further investigation is required to experimentally confirm these predictions.

5. Conclusions

We used EIS to compare the complex impedance of MT- and tubulin-containing solutions. A physiologically relevant, high ionic strength buffer (BRB80) created a high noise, low impedance background, which was countered through the use of physiological concentrations of tubulin. While the presence of MTs increased solution capacitance, unpolymerised tubulin did not have any appreciable effect. In a study that is the first of its kind to the best of our knowledge, we determined the capacitance and resistance of the MT network at physiological tubulin concentrations to be 1.27 × 10−5 F and 9.74 × 104 Ω. These values correspond to an effective resistance per unit volume of 3.71 × 1010 Ω/L and effective capacitance per unit volume of 7.65 F/L. We envision a dual electrical role for MTs in the cell, that of charge storage devices across a broad frequency spectrum (acting as storage locations for ions), and of charge transporters (bionanowires) in the frequency region between 20 and 60 Hz. Our findings also indicate that the electrical properties of tubulin dimers change as they polymerise, revealing the potential impact of MT nucleation and polymerisation on the cellular charge distribution. Our work shows that by storing charge and attenuating local ion distributions, microtubules play a crucial role in governing the bioelectric properties of the cell.

Supplementary Materials

The following are available online at https://www.mdpi.com/2079-4991/10/2/265/s1. Figure S1. 3D printed holders used to fabricate and align the parallel-plate contact device. (a) Top view (left) and (b) side view (right) of holder for the parallel plate device used to perform impedance measurements. (c) Top view (left) and side view (right) of slider used to position the double-sided tape exactly to fabricate the device. (d) Top view (left) and side view (right) of the holder used to position the upper contact precisely on the lower contact. Figure S2. Validation of parallel-plate contact device using 0.5 mM electrolytic solutions. (a) Imaginary component of impedance for electrolytic solutions at 0.5 mM and de-ionised water. (b) Real component of impedance for electrolytic solutions at 100 mM and de-ionised water. Data display average values collected between 15 and 21 times. Error bars represent standard deviation. Figure S3. Example of microtubule and tubulin subtraction with backgrounds to display typical impedance. Figure S4. No ‘reversal’ in the resistive behaviour of microtubules is observed between 10 and 100 Hz. Graphs showing differences in the real component of impedance as a function of decreasing input AC frequency at total tubulin concentrations of (a) 22.225 µM, (b) 2.222 µM, (c) 0.222 µM, (d) comparison of the effect of paclitaxel and colchicine on impedance. Figure S5. One-sample t-tests were performed to determine if the impedance difference values were significantly above zero. This was carried out using the t-test function (‘ttest’) within MATLAB. Graphs showing the variation of obtained p-values for the imaginary components of impedance in (a) tubulin and (b) MT-containing solutions. Graphs showing the variation of obtained p-values for the real components of impedance in (a) tubulin and (b) MT-containing solutions.

Author Contributions

Conceptualization, methodology, software, formal analysis, validation, investigation, data curation, writing—original draft, visualisation, project administration, funding acquisition, A.P.K. Conceptualization, methodology, investigation, validation, writing—original draft, data curation, S.D.P. Methodology, software, formal analysis, validation, investigation, data curation, A.F.B. Conceptualization, investigation, J.P. Software, validation, formal analysis, K.G.S. Methodology, investigation, U.M. Resources, project administration, funding acquisition, J.D.L. Conceptualization, software, formal analysis, validation, investigation, writing—review and editing, visualization, V.R. Conceptualization, methodology, validation, data curation, writing—review and editing, resources, supervision, project administration, funding acquisition, K.S. Conceptualization, methodology, validation, formal analysis, data curation, writing—review and editing, resources, supervision, project administration, funding acquisition, J.A.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Novocure LLC, CMC Microsystems, voucher number 6067 and NSERC, grant number 06630. The APC was waived by the publisher for this article.

Acknowledgments

A.P.K. acknowledges Najia Mahdi for helping him with reactive ion etching. A.F.B. acknowledges personnel at The Shack at the University of Alberta and particularly Tristan Stark for teaching him 3D printing. A.F.B. and S.D.P. acknowledge Liam McRae for showing them how to use the impedance analyser. The authors acknowledge Philip Winter, Sheng Zeng, Dayal Pyari Srivastava, Ze’ev Bomzon, Frank Hegmann, Piyush Kar and Kris Carlson for useful discussions.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Mershin, A.; Kolomenski, A.A.; Schuessler, H.A.; Nanopoulos, D.V. Tubulin dipole moment, dielectric constant and quantum behavior: Computer simulations, experimental results and suggestions. Biosystems 2004, 77, 73–85. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Van den Heuvel, M.; De Graaff, M.; Lemay, S.; Dekker, C. Electrophoresis of individual microtubules in microchannels. Proc. Natl. Acad. Sci. USA 2007, 104, 7770–7775. [Google Scholar] [CrossRef] [Green Version]
  3. Timmons, J.J.; Preto, J.; Tuszynski, J.A.; Wong, E.T. Tubulin’s response to external electric fields by molecular dynamics simulations. PLoS ONE 2018, 13, e0202141. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Preto, J.; Pettini, M.; Tuszynski, J.A. Possible role of electrodynamic interactions in long-distance biomolecular recognition. Phys. Rev. E 2015, 91, 052710. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Barvitenko, N.; Lawen, A.; Aslam, M.; Pantaleo, A.; Saldanha, C.; Skverchinskaya, E.; Regolini, M.; Tuszynski, J.A. Integration of intracellular signaling: Biological analogues of wires, processors and memories organized by a centrosome 3D reference system. Biosystems 2018, 173, 191–206. [Google Scholar] [CrossRef]
  6. Friesen, D.E.; Craddock, T.J.; Kalra, A.P.; Tuszynski, J.A. Biological wires, communication systems, and implications for disease. Biosystems 2015, 127, 14–27. [Google Scholar] [CrossRef]
  7. Sekulić, D.L.; Satarić, B.M.; Tuszynski, J.A.; Satarić, M.V. Nonlinear ionic pulses along microtubules. Eur. Phys. J. E 2011, 34, 49. [Google Scholar] [CrossRef]
  8. Satarić, M.; Ilić, D.; Ralević, N.; Tuszynski, J.A. A nonlinear model of ionic wave propagation along microtubules. Eur. Biophys. J. 2009, 38, 637–647. [Google Scholar] [CrossRef]
  9. del Rocío Cantero, M.; Etchegoyen, C.V.; Perez, P.L.; Scarinci, N.; Cantiello, H.F. Bundles of Brain Microtubules Generate Electrical Oscillations. Sci. Rep. 2018, 8, 11899. [Google Scholar] [CrossRef]
  10. del Rocío Cantero, M.; Perez, P.L.; Smoler, M.; Etchegoyen, C.V.; Cantiello, H.F. Electrical oscillations in two-dimensional microtubular structures. Sci. Rep. 2016, 6, 27143. [Google Scholar] [CrossRef] [PubMed]
  11. Cifra, M.; Pokorný, J.; Havelka, D.; Kučera, O. Electric field generated by axial longitudinal vibration modes of microtubule. BioSystems 2010, 100, 122–131. [Google Scholar] [CrossRef]
  12. Uppalapati, M.; Huang, Y.M.; Jackson, T.N.; Hancock, W.O. Microtubule alignment and manipulation using AC electrokinetics. Small 2008, 4, 1371–1381. [Google Scholar] [CrossRef]
  13. Minoura, I.; Muto, E. Dielectric measurement of individual microtubules using the electroorientation method. Biophys. J. 2006, 90, 3739–3748. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Stracke, R.; Böhm, K.; Wollweber, L.; Tuszynski, J.; Unger, E. Analysis of the migration behaviour of single microtubules in electric fields. Biochem. Biophys. Res. Commun. 2002, 293, 602–609. [Google Scholar] [CrossRef]
  15. Isozaki, N.; Ando, S.; Nakahara, T.; Shintaku, H.; Kotera, H.; Meyhöfer, E.; Yokokawa, R. Control of microtubule trajectory within an electric field by altering surface charge density. Sci. Rep. 2015, 5, 7669. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Dujovne, I.; van den Heuvel, M.; Shen, Y.; de Graaff, M.; Dekker, C. Velocity modulation of microtubules in electric fields. Nano Lett. 2008, 8, 4217–4220. [Google Scholar] [CrossRef] [PubMed]
  17. Kirson, E.D.; Gurvich, Z.; Schneiderman, R.; Dekel, E.; Itzhaki, A.; Wasserman, Y.; Schatzberger, R.; Palti, Y. Disruption of cancer cell replication by alternating electric fields. Cancer Res. 2004, 64, 3288–3295. [Google Scholar] [CrossRef] [Green Version]
  18. Giladi, M.; Schneiderman, R.S.; Voloshin, T.; Porat, Y.; Munster, M.; Blat, R.; Sherbo, S.; Bomzon, Z.; Urman, N.; Itzhaki, A. Mitotic spindle disruption by alternating electric fields leads to improper chromosome segregation and mitotic catastrophe in cancer cells. Sci. Rep. 2015, 5, 18046. [Google Scholar] [CrossRef]
  19. Tuszynski, J.; Wenger, C.; Friesen, D.; Preto, J. An overview of sub-cellular mechanisms involved in the action of TTFields. Int. J. Environ. Res. Public Health 2016, 13, 1128. [Google Scholar] [CrossRef] [Green Version]
  20. Santelices, I.B.; Friesen, D.E.; Bell, C.; Hough, C.M.; Xiao, J.; Kalra, A.; Kar, P.; Freedman, H.; Rezania, V.; Lewis, J.D. Response to Alternating Electric Fields of Tubulin Dimers and Microtubule Ensembles in Electrolytic Solutions. Sci. Rep. 2017, 7, 9594. [Google Scholar] [CrossRef] [Green Version]
  21. Shelden, E.; Wadsworth, P. Observation and quantification of individual microtubule behavior in vivo: Microtubule dynamics are cell-type specific. J. Cell Biol. 1993, 120, 935–945. [Google Scholar] [CrossRef] [PubMed]
  22. Van de Water, L.; Olmsted, J. The quantitation of tubulin in neuroblastoma cells by radioimmunoassay. J. Biol. Chem. 1980, 255, 10744–10751. [Google Scholar] [PubMed]
  23. Olmsted, J.; Borisy, G. Ionic and nucleotide requirements for microtubule polymerization in vitro. Biochemistry 1975, 14, 2996–3005. [Google Scholar] [CrossRef] [PubMed]
  24. Weisenberg, R.C. Microtubule formation in vitro in solutions containing low calcium concentrations. Science 1972, 177, 1104–1105. [Google Scholar] [CrossRef]
  25. Böhm, K.J.; Mavromatos, N.E.; Michette, A.; Stracke, R.; Unger, E. Movement and alignment of microtubules in electric fields and electric-dipole-moment estimates. Electromagn. Biol. Med. 2005, 24, 319–330. [Google Scholar] [CrossRef]
  26. Umnov, M.; Palusinski, O.; Deymier, P.; Guzman, R.; Hoying, J.; Barnaby, H.; Yang, Y.; Raghavan, S. Experimental evaluation of electrical conductivity of microtubules. J. Mater. Sci. 2007, 42, 373–378. [Google Scholar] [CrossRef]
  27. van Eunen, K.; Bouwman, J.; Daran-Lapujade, P.; Postmus, J.; Canelas, A.B.; Mensonides, F.I.; Orij, R.; Tuzun, I.; van den Brink, J.; Smits, G.J. Measuring enzyme activities under standardized in vivo-like conditions for systems biology. Febs J. 2010, 277, 749–760. [Google Scholar] [CrossRef] [Green Version]
  28. Theillet, F.-X.; Binolfi, A.; Frembgen-Kesner, T.; Hingorani, K.; Sarkar, M.; Kyne, C.; Li, C.; Crowley, P.B.; Gierasch, L.; Pielak, G.J. Physicochemical properties of cells and their effects on intrinsically disordered proteins (IDPs). Chem. Rev. 2014, 114, 6661–6714. [Google Scholar] [CrossRef]
  29. Fritzsche, W.; Böhm, K.; Unger, E.; Köhler, J.M. Making electrical contact to single molecules. Nanotechnology 1998, 9, 177. [Google Scholar] [CrossRef]
  30. Fritzsche, W.; Köhler, J.M.; Böhm, K.J.; Unger, E.; Wagner, T.; Kirsch, R.; Mertig, M.; Pompe, W. Wiring of metallized microtubules by electron beam-induced structuring. Nanotechnology 1999, 10, 331. [Google Scholar] [CrossRef]
  31. Gell, C.; Bormuth, V.; Brouhard, G.J.; Cohen, D.N.; Diez, S.; Friel, C.T.; Helenius, J.; Nitzsche, B.; Petzold, H.; Ribbe, J. Microtubule dynamics reconstituted in vitro and imaged by single-molecule fluorescence microscopy. In Methods in Cell Biology; Elsevier: New York, NY, USA, 2010; Volume 95, pp. 221–245. [Google Scholar]
  32. Taylor, S.; Gileadi, E. Physical interpretation of the Warburg impedance. Corrosion 1995, 51, 664–671. [Google Scholar] [CrossRef]
  33. Ramos-Barrado, J.; Galan Montenegro, P.; Cambón, C.C. A generalized Warburg impedance for a nonvanishing relaxation process. J. Chem. Phys. 1996, 105, 2813–2815. [Google Scholar] [CrossRef] [Green Version]
  34. Sanabria, H.; Miller, J.H., Jr. Relaxation processes due to the electrode-electrolyte interface in ionic solutions. Phys. Rev. E 2006, 74, 051505. [Google Scholar] [CrossRef] [PubMed]
  35. Lima, L.F.; Vieira, A.L.; Mukai, H.; Andrade, C.M.; Fernandes, P.R. Electric impedance of aqueous KCl and NaCl solutions: Salt concentration dependence on components of the equivalent electric circuit. J. Mol. Liq. 2017, 241, 530–539. [Google Scholar] [CrossRef]
  36. Lenzi, E.; Fernandes, P.; Petrucci, T.; Mukai, H.; Ribeiro, H.; Lenzi, M.; Gonçalves, G. Anomalous diffusion and electrical response of ionic solutions. Int. J. Electrochem. Sci 2013, 8, 2849–2862. [Google Scholar]
  37. Wilson, L.; Miller, H.P.; Farrell, K.W.; Snyder, K.B.; Thompson, W.C.; Purich, D.L. Taxol stabilization of microtubules in vitro: Dynamics of tubulin addition and loss at opposite microtubule ends. Biochemistry 1985, 24, 5254–5262. [Google Scholar] [CrossRef]
  38. Arnal, I.; Wade, R.H. How does taxol stabilize microtubules? Curr. Biol. 1995, 5, 900–908. [Google Scholar] [CrossRef] [Green Version]
  39. Mikhaylova, M.; Cloin, B.M.; Finan, K.; Van Den Berg, R.; Teeuw, J.; Kijanka, M.M.; Sokolowski, M.; Katrukha, E.A.; Maidorn, M.; Opazo, F. Resolving bundled microtubules using anti-tubulin nanobodies. Nat. Commun. 2015, 6, 7933. [Google Scholar] [CrossRef]
  40. Herold, C.; Leduc, C.; Stock, R.; Diez, S.; Schwille, P. Long-Range Transport of Giant Vesicles along Microtubule Networks. ChemPhysChem 2012, 13, 1001–1006. [Google Scholar] [CrossRef]
  41. Kulić, I.M.; Brown, A.E.; Kim, H.; Kural, C.; Blehm, B.; Selvin, P.R.; Nelson, P.C.; Gelfand, V.I. The role of microtubule movement in bidirectional organelle transport. Proc. Natl. Acad. Sci. USA 2008, 105, 10011–10016. [Google Scholar] [CrossRef] [Green Version]
  42. Margolis, R.L.; Wilson, L. Addition of colchicine-tubulin complex to microtubule ends: The mechanism of substoichiometric colchicine poisoning. Proc. Natl. Acad. Sci. USA 1977, 74, 3466–3470. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Olmsted, J.; Borisy, G. Characterization of microtubule assembly in porcine brain extracts by viscometry. Biochemistry 1973, 12, 4282–4289. [Google Scholar] [CrossRef] [PubMed]
  44. Satarić, M.; Sekulić, D.; Živanov, M. Solitonic ionic currents along microtubules. J. Comput. Theor. Nanosci. 2010, 7, 2281–2290. [Google Scholar] [CrossRef]
  45. Freedman, H.; Rezania, V.; Priel, A.; Carpenter, E.; Noskov, S.Y.; Tuszynski, J.A. Model of ionic currents through microtubule nanopores and the lumen. Phys. Rev. E 2010, 81, 051912. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Sanabria, H.; Miller, J.H., Jr.; Mershin, A.; Luduena, R.F.; Kolomenski, A.A.; Schuessler, H.A.; Nanopoulos, D.V. Impedance spectroscopy of α-β tubulin heterodimer suspensions. Biophys. J. 2006, 90, 4644–4650. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Shen, C.; Guo, W. Ion Permeability of a Microtubule in Neuron Environment. J. Phys. Chem. Lett. 2018, 9, 2009–2014. [Google Scholar] [CrossRef]
  48. Priel, A.; Tuszyński, J. A nonlinear cable-like model of amplified ionic wave propagation along microtubules. Epl (Europhys. Lett.) 2008, 83, 68004. [Google Scholar] [CrossRef]
  49. Tuszyński, J.; Brown, J.; Crawford, E.; Carpenter, E.; Nip, M.; Dixon, J.; Satarić, M. Molecular dynamics simulations of tubulin structure and calculations of electrostatic properties of microtubules. Math. Comput. Model. 2005, 41, 1055–1070. [Google Scholar] [CrossRef]
  50. Sekulić, D.L.; Satarić, M.V. Microtubule as nanobioelectronic nonlinear circuit. Serb. J. Electr. Eng. 2012, 9, 107–119. [Google Scholar] [CrossRef]
  51. Manning, G.S. Limiting laws and counterion condensation in polyelectrolyte solutions I. Colligative properties. J. Chem. Phys. 1969, 51, 924–933. [Google Scholar] [CrossRef]
  52. Tuszyński, J.; Portet, S.; Dixon, J.; Luxford, C.; Cantiello, H. Ionic wave propagation along actin filaments. Biophys. J. 2004, 86, 1890–1903. [Google Scholar] [CrossRef] [Green Version]
  53. Priel, A.; Ramos, A.J.; Tuszynski, J.A.; Cantiello, H.F. A biopolymer transistor: Electrical amplification by microtubules. Biophys. J. 2006, 90, 4639–4643. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Aouaichia, M.; McCullen, N.; Bowen, C.R.; Almond, D.P.; Budd, C.; Bouamrane, R. Understanding the anomalous frequency responses of composite materials using very large random resistor-capacitor networks. Eur. Phys. J. B 2017, 90, 39. [Google Scholar] [CrossRef] [Green Version]
  55. Panteny, S.; Stevens, R.; Bowen, C. The frequency dependent permittivity and AC conductivity of random electrical networks. Ferroelectrics 2005, 319, 199–208. [Google Scholar] [CrossRef]
  56. Parri, H.R.; Gould, T.M.; Crunelli, V. Spontaneous astrocytic Ca 2+ oscillations in situ drive NMDAR-mediated neuronal excitation. Nat. Neurosci. 2001, 4, 803. [Google Scholar] [CrossRef]
  57. Regula, C.; Berlin, R. Microtubule assembly and disassembly at alkaline pH. J. Cell Biol. 1981, 89, 45–53. [Google Scholar] [CrossRef] [Green Version]
  58. Borisy, G.; Marcum, J.; Olmsted, J.; Murphy, D.; Johnson, K. Purification of tubulin and associated high molecular weight proteins from porcine brain and characterization of microtubule assembly in vitro. Ann. New York Acad. Sci. USA 1975, 253, 107–132. [Google Scholar] [CrossRef]
  59. Larsson, H.; Wallin, M.; Edström, A. Induction of a sheet polymer of tubulin by Zn2+. Exp. Cell Res. 1976, 100, 104–110. [Google Scholar] [CrossRef]
  60. Nicholson, W.V.; Lee, M.; Downing, K.H.; Nogales, E. Cryo-electron microscopy of GDP-tubulin rings. Cell Biochem. Biophys. 1999, 31, 175–183. [Google Scholar] [CrossRef]
  61. Zhai, Y.; Kronebusch, P.J.; Simon, P.M.; Borisy, G.G. Microtubule dynamics at the G2/M transition: Abrupt breakdown of cytoplasmic microtubules at nuclear envelope breakdown and implications for spindle morphogenesis. J. Cell Biol. 1996, 135, 201–214. [Google Scholar] [CrossRef] [Green Version]
  62. Zhai, Y.; Borisy, G.G. Quantitative determination of the proportion of microtubule polymer present during the mitosis-interphase transition. J. Cell Sci. 1994, 107, 881–890. [Google Scholar] [PubMed]
  63. Kline-Smith, S.L.; Walczak, C.E. Mitotic spindle assembly and chromosome segregation: Refocusing on microtubule dynamics. Mol. Cell 2004, 15, 317–327. [Google Scholar] [CrossRef] [PubMed]
  64. Akhmanova, A.; Dogterom, M. Kinesins lead aging microtubules to catastrophe. Cell 2011, 147, 966–968. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Bugiel, M.; Mitra, A.; Girardo, S.; Diez, S.; Schäffer, E. Measuring microtubule supertwist and defects by three-dimensional-force-clamp tracking of single kinesin-1 motors. Nano Lett. 2018, 18, 1290–1295. [Google Scholar] [CrossRef]
  66. Diez, S.; Reuther, C.; Dinu, C.; Seidel, R.; Mertig, M.; Pompe, W.; Howard, J. Stretching and transporting DNA molecules using motor proteins. Nano Lett. 2003, 3, 1251–1254. [Google Scholar] [CrossRef]
  67. Hess, H.; Saper, G. Engineering with biomolecular motors. Acc. Chem. Res. 2018, 51, 3015–3022. [Google Scholar] [CrossRef]
  68. Song, W.; Zhu, J.; Kong, W.; Möhwald, H.; Li, J. Different Microtubule Structures Assembled by Kinesin Motors. Langmuir 2018, 34, 9768–9773. [Google Scholar] [CrossRef]
  69. Tas, R.P.; Chen, C.-Y.; Katrukha, E.A.; Vleugel, M.; Kok, M.; Dogterom, M.; Akhmanova, A.; Kapitein, L.C. Guided by Light: Optical Control of Microtubule Gliding Assays. Nano Lett. 2018, 18, 7524–7528. [Google Scholar] [CrossRef] [Green Version]
  70. Fujimoto, K.; Kitamura, M.; Yokokawa, M.; Kanno, I.; Kotera, H.; Yokokawa, R. Colocalization of quantum dots by reactive molecules carried by motor proteins on polarized microtubule arrays. Acs Nano 2012, 7, 447–455. [Google Scholar] [CrossRef]
  71. Adrian, M.; Nijenhuis, W.; Hoogstraaten, R.I.; Willems, J.; Kapitein, L.C. A phytochrome-derived photoswitch for intracellular transport. Acs Synth. Biol. 2017, 6, 1248–1256. [Google Scholar] [CrossRef]
  72. Lam, A.T.-C.; Tsitkov, S.; Zhang, Y.; Hess, H. Reversibly bound kinesin-1 motor proteins propelling microtubules demonstrate dynamic recruitment of active building blocks. Nano Lett. 2018, 18, 1530–1534. [Google Scholar] [CrossRef] [PubMed]
  73. Yokokawa, R.; Takeuchi, S.; Kon, T.; Nishiura, M.; Sutoh, K.; Fujita, H. Unidirectional transport of kinesin-coated beads on microtubules oriented in a microfluidic device. Nano Lett. 2004, 4, 2265–2270. [Google Scholar] [CrossRef]
  74. Fujimoto, K.; Morita, Y.; Iino, R.; Tomishige, M.; Shintaku, H.; Kotera, H.; Yokokawa, R. Simultaneous Observation of Kinesin-Driven Microtubule Motility and Binding of Adenosine Triphosphate Using Linear Zero-Mode Waveguides. Acs Nano 2018, 12, 11975–11985. [Google Scholar] [CrossRef] [PubMed]
  75. Li, J.; Jia, Y.; Dong, W.; Feng, X.; Fei, J.; Li, J. Transporting a Tube in a Tube. Nano Lett. 2014, 14, 6160–6164. [Google Scholar] [CrossRef]
  76. Wieczorek, M.; Chaaban, S.; Brouhard, G.J. Macromolecular crowding pushes catalyzed microtubule growth to near the theoretical limit. Cell. Mol. Bioeng. 2013, 6, 383–392. [Google Scholar] [CrossRef]
  77. Sozański, K.; Ruhnow, F.; Wiśniewska, A.; Tabaka, M.; Diez, S.; Hołyst, R. Small crowders slow down kinesin-1 stepping by hindering motor domain diffusion. Phys. Rev. Lett. 2015, 115, 218102. [Google Scholar] [CrossRef]
Nanomaterials 10 00265 g001 550
Figure 1. A parallel-plate contact device to measure the impedance properties of microtubules (MTs) compared to tubulin. The operation of the parallel plate device showing (a) top view (left) and side view (right). The upper and lower contacts, double-sided tape and solution are labelled in green, grey and blue, respectively. (b) Imaginary component of impedance for electrolytic solutions at 100 mM and de-ionised water. (c) Real component of impedance for electrolytic solutions at 100 mM and de-ionised water. Data display average values collected between 15 and 21 times. Error bars represent standard deviation.
Figure 1. A parallel-plate contact device to measure the impedance properties of microtubules (MTs) compared to tubulin. The operation of the parallel plate device showing (a) top view (left) and side view (right). The upper and lower contacts, double-sided tape and solution are labelled in green, grey and blue, respectively. (b) Imaginary component of impedance for electrolytic solutions at 100 mM and de-ionised water. (c) Real component of impedance for electrolytic solutions at 100 mM and de-ionised water. Data display average values collected between 15 and 21 times. Error bars represent standard deviation.
Nanomaterials 10 00265 g001
Nanomaterials 10 00265 g002 550
Figure 2. Microtubule imaging at different tubulin concentrations. Polymerisation was performed using 45 µM tubulin, and MTs were stabilised with 50 µM paclitaxel, and subsequently diluted to a final concentration of (a) 0.222 µM tubulin (b) 2.222 µM tubulin (c) 22.225 µM tubulin, respectively. Scale bars represent 10 μm.
Figure 2. Microtubule imaging at different tubulin concentrations. Polymerisation was performed using 45 µM tubulin, and MTs were stabilised with 50 µM paclitaxel, and subsequently diluted to a final concentration of (a) 0.222 µM tubulin (b) 2.222 µM tubulin (c) 22.225 µM tubulin, respectively. Scale bars represent 10 μm.
Nanomaterials 10 00265 g002
Nanomaterials 10 00265 g003 550
Figure 3. Examples of raw values of imaginary component of impedance in solutions containing (a) MTs and (b) tubulin, for the purpose of displaying typical impedance values. Mean differences in the imaginary component of impedance as a function of input AC (alternating current) frequency at total tubulin concentrations of (c) 22.225 µM (n = 22 experiments for tubulin, n = 21 for MTs), (d) 2.222 µM (n = 35 experiments for tubulin, n = 49 for MTs) (e) 0.222 µM (n = 35 experiments for tubulin, n = 49 for MTs), (f) comparison of the effects of paclitaxel (BRB80T) and colchicine (BRB80C, n = 49 experiments for BRB80T, n = 35 for BRB80C, n = 84 experiments for BRB80). Error-bars represent standard deviation.
Figure 3. Examples of raw values of imaginary component of impedance in solutions containing (a) MTs and (b) tubulin, for the purpose of displaying typical impedance values. Mean differences in the imaginary component of impedance as a function of input AC (alternating current) frequency at total tubulin concentrations of (c) 22.225 µM (n = 22 experiments for tubulin, n = 21 for MTs), (d) 2.222 µM (n = 35 experiments for tubulin, n = 49 for MTs) (e) 0.222 µM (n = 35 experiments for tubulin, n = 49 for MTs), (f) comparison of the effects of paclitaxel (BRB80T) and colchicine (BRB80C, n = 49 experiments for BRB80T, n = 35 for BRB80C, n = 84 experiments for BRB80). Error-bars represent standard deviation.
Nanomaterials 10 00265 g003
Nanomaterials 10 00265 g004 550
Figure 4. Examples of raw values of real component of impedance in solutions containing (a) MTs and (b) tubulin, for the purpose of displaying typical impedance values. Mean differences in the real component of impedance as a function of input AC frequency at total tubulin concentrations of (c) 22.225 µM (n = 22 experiments for tubulin, n = 21 for MTs), (d) 2.222 µM (n = 35 experiments for tubulin, n = 49 for MTs) (e) 0.222 µM (n = 35 experiments for tubulin, n = 49 for MTs), (f) comparison of the effects of paclitaxel (BRB80T) and colchicine (BRB80C, n = 49 experiments for BRB80T, n = 35 for BRB80C, n = 84 experiments for BRB80). Error-bars represent standard deviation.
Figure 4. Examples of raw values of real component of impedance in solutions containing (a) MTs and (b) tubulin, for the purpose of displaying typical impedance values. Mean differences in the real component of impedance as a function of input AC frequency at total tubulin concentrations of (c) 22.225 µM (n = 22 experiments for tubulin, n = 21 for MTs), (d) 2.222 µM (n = 35 experiments for tubulin, n = 49 for MTs) (e) 0.222 µM (n = 35 experiments for tubulin, n = 49 for MTs), (f) comparison of the effects of paclitaxel (BRB80T) and colchicine (BRB80C, n = 49 experiments for BRB80T, n = 35 for BRB80C, n = 84 experiments for BRB80). Error-bars represent standard deviation.
Nanomaterials 10 00265 g004
Nanomaterials 10 00265 g005 550
Figure 5. Graphs showing differences in the (a) imaginary from Figure 3 and (b) real component of impedance from Figure 4 as a function of tubulin concentration at input AC frequencies of 1 Hz, 10 Hz, 100 Hz, 1 kHz, 10 kHz and 86 kHz. Graphs display average values. Error-bars represent standard deviation.
Figure 5. Graphs showing differences in the (a) imaginary from Figure 3 and (b) real component of impedance from Figure 4 as a function of tubulin concentration at input AC frequencies of 1 Hz, 10 Hz, 100 Hz, 1 kHz, 10 kHz and 86 kHz. Graphs display average values. Error-bars represent standard deviation.
Nanomaterials 10 00265 g005
Nanomaterials 10 00265 g006 550
Figure 6. Zoomed in view of the mean differences in the real component of impedance as a function of decreasing input AC frequency at total tubulin concentrations of (a) 22.225 µM, (b) 2.222 µM, (c) 0.222 µM, (d) comparison of the effect of paclitaxel and colchicine on impedance. (e) A logarithmic plot obtained by translating the graphs (a), (b) and (c) upwards. The translation is performed by adding (1+minimum MT solution resistance) to the resistance of each MT concentration. A resistance reversal between 20–60 Hz is observed, with a peak at 39 Hz for the 22.225 µM concentration. Error-bars represent standard deviation.
Figure 6. Zoomed in view of the mean differences in the real component of impedance as a function of decreasing input AC frequency at total tubulin concentrations of (a) 22.225 µM, (b) 2.222 µM, (c) 0.222 µM, (d) comparison of the effect of paclitaxel and colchicine on impedance. (e) A logarithmic plot obtained by translating the graphs (a), (b) and (c) upwards. The translation is performed by adding (1+minimum MT solution resistance) to the resistance of each MT concentration. A resistance reversal between 20–60 Hz is observed, with a peak at 39 Hz for the 22.225 µM concentration. Error-bars represent standard deviation.
Nanomaterials 10 00265 g006
Nanomaterials 10 00265 g007 550
Figure 7. The equivalent electrical circuit model representing the microtubule network as a parallel RC circuit, with network resistance RMT and capacitance CMT. The external element has impedance Z0, while solution has impedance Zs. RH is the small constant resistance that is ascribed to small fraction of unpolymerised tubulin that is present in MT containing solutions.
Figure 7. The equivalent electrical circuit model representing the microtubule network as a parallel RC circuit, with network resistance RMT and capacitance CMT. The external element has impedance Z0, while solution has impedance Zs. RH is the small constant resistance that is ascribed to small fraction of unpolymerised tubulin that is present in MT containing solutions.
Nanomaterials 10 00265 g007
Nanomaterials 10 00265 g008 550
Figure 8. Mean differences of (a) imaginary and (b) real impedance curves for 0.222 µM, 2.222 µM and 22.222 µM, are fitted with the model described in Equation (6) and Figure 7. Fit parameters and confidence intervals are displayed in Table 2. The region between 1–100 Hz was not fit because of the negative differences in resistance from background BRB80T solutions.
Figure 8. Mean differences of (a) imaginary and (b) real impedance curves for 0.222 µM, 2.222 µM and 22.222 µM, are fitted with the model described in Equation (6) and Figure 7. Fit parameters and confidence intervals are displayed in Table 2. The region between 1–100 Hz was not fit because of the negative differences in resistance from background BRB80T solutions.
Nanomaterials 10 00265 g008
Nanomaterials 10 00265 g009 550
Figure 9. Schematic of charge transport along and across an MT. (a) A representation of charge flow across the MT cross section through nanopores present between adjacent protofilaments. (b) A representation of charge flow through both inner and outer modes along an MT. Arrows depict charge flow via both mechanisms, enabling MT charge storage across a broad spectrum of frequencies, and charge transport at low AC frequencies in the cell. (c) Side view (left) and top view (right) of the tubulin dimer, displaying distribution of electrostatic potential at different locations. The negatively charged C-termini face towards the solution and contains ~50% of the total negative charge on a tubulin dimer.
Figure 9. Schematic of charge transport along and across an MT. (a) A representation of charge flow across the MT cross section through nanopores present between adjacent protofilaments. (b) A representation of charge flow through both inner and outer modes along an MT. Arrows depict charge flow via both mechanisms, enabling MT charge storage across a broad spectrum of frequencies, and charge transport at low AC frequencies in the cell. (c) Side view (left) and top view (right) of the tubulin dimer, displaying distribution of electrostatic potential at different locations. The negatively charged C-termini face towards the solution and contains ~50% of the total negative charge on a tubulin dimer.
Nanomaterials 10 00265 g009
Table
Table 1. Volumes of tubulin and buffer solution (BRB80T or BRB80C) used to stabilise microtubules (BRB80T) or free tubulin (BRB80C) in solution.
Table 1. Volumes of tubulin and buffer solution (BRB80T or BRB80C) used to stabilise microtubules (BRB80T) or free tubulin (BRB80C) in solution.
Tubulin Concentration (μM)Volume of BRB80T or BRB80C (μL)Tubulin Volume (μL)
0.22299.50.5
2.222955
22.22555
Table
Table 2. Fit parameters attained by fitting the real and imaginary components of impedance to Equation (7). Fit parameters represent effective capacitance CMT, and resistance RMT introduced into the solution through the addition of the MT network at different concentrations. RH is the small constant resistance that is ascribed to small fraction of unpolymerised tubulin that is present in MT-containing solutions. γR and γC describe the effective resistance and capacitance per unit volume introduced by the microtubule network. δ RMT, δ CMT and δ RH correspond to 95% confidence intervals for the fit parameters. The values δ γR and δ γC correspond to the uncertainties in the resistance and capacitance per unit volume. Corresponding graphs are displayed in Figure 8.
Table 2. Fit parameters attained by fitting the real and imaginary components of impedance to Equation (7). Fit parameters represent effective capacitance CMT, and resistance RMT introduced into the solution through the addition of the MT network at different concentrations. RH is the small constant resistance that is ascribed to small fraction of unpolymerised tubulin that is present in MT-containing solutions. γR and γC describe the effective resistance and capacitance per unit volume introduced by the microtubule network. δ RMT, δ CMT and δ RH correspond to 95% confidence intervals for the fit parameters. The values δ γR and δ γC correspond to the uncertainties in the resistance and capacitance per unit volume. Corresponding graphs are displayed in Figure 8.
[Tub] (µM).CMT (F).δCMT (F)RMT (Ω)δRMT (Ω)RH (Ω)δRH (Ω)γR
(Ω/L)		γC
(F/L)		δγR
(Ω/L)		δγC
(F/L)
22.2221.27 × 10−51.48 × 10−79.74 × 1041.18 × 1042.1240.613.71 × 10107.654.49 × 1090.056
2.2221.25 × 10−51.67 × 10−71.00 × 1051.40 × 1040.6134.793.81 × 10104.765.33 × 1090.063
0.2222.01 × 10−53.38 × 10−79.97 × 1042.82 × 1040.4131.953.80 × 10104.831.07 × 10100.12

Share and Cite

MDPI and ACS Style

Kalra, A.P.; Patel, S.D.; Bhuiyan, A.F.; Preto, J.; Scheuer, K.G.; Mohammed, U.; Lewis, J.D.; Rezania, V.; Shankar, K.; Tuszynski, J.A. Investigation of the Electrical Properties of Microtubule Ensembles under Cell-Like Conditions. Nanomaterials 2020, 10, 265. https://doi.org/10.3390/nano10020265

AMA Style

Kalra AP, Patel SD, Bhuiyan AF, Preto J, Scheuer KG, Mohammed U, Lewis JD, Rezania V, Shankar K, Tuszynski JA. Investigation of the Electrical Properties of Microtubule Ensembles under Cell-Like Conditions. Nanomaterials. 2020; 10(2):265. https://doi.org/10.3390/nano10020265

Chicago/Turabian Style

Kalra, Aarat P., Sahil D. Patel, Asadullah F. Bhuiyan, Jordane Preto, Kyle G. Scheuer, Usman Mohammed, John D. Lewis, Vahid Rezania, Karthik Shankar, and Jack A. Tuszynski. 2020. "Investigation of the Electrical Properties of Microtubule Ensembles under Cell-Like Conditions" Nanomaterials 10, no. 2: 265. https://doi.org/10.3390/nano10020265

APA Style

Kalra, A. P., Patel, S. D., Bhuiyan, A. F., Preto, J., Scheuer, K. G., Mohammed, U., Lewis, J. D., Rezania, V., Shankar, K., & Tuszynski, J. A. (2020). Investigation of the Electrical Properties of Microtubule Ensembles under Cell-Like Conditions. Nanomaterials, 10(2), 265. https://doi.org/10.3390/nano10020265

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