The Diamagnetic Susceptibility of the Tubulin Dimer
An approximate value of the diamagnetic anisotropy of the tubulin dimer, , has been determined assuming axial symmetry and that only the α-helices and β-sheets contribute to the anisotropy. Two approaches have been utilized: (a) using the value for the for an α-helical peptide bond given by Pauling (1979) and (b) using the previously determined anisotropy of fibrinogen as a calibration standard. The JT−2 obtained from these measurements are similar to within 20%. Although Cotton-Mouton measurements alone cannot be used to estimate directly, the value we measured, T−2cm2mg−1, is consistent with the above estimate for . The method utilized for the determination of the tubulin dimer diamagnetic susceptibility is applicable to other proteins and macromolecular assemblies as well.
Microtubules, MT, are elongated macromolecular structures composed of protofilaments, of tubulin dimers assembled in long hollow tubes. The tubulin dimer is kidney shaped and has an approximate molecular weight of 55 kDa. The assembled microtubule is on average composed of 13 protofilaments which make up a hollow tube with a maximum diameter of nm . The length of this assembly in vivo is variable but can reach several microns.
The aligning effects of magnetic fields on macromolecular microtubules and other rigid fibrillar molecular assemblies are well known [1–3]. The driving force for magnetic orientation is the diamagnetic anisotropy of the tubulin subunits (dimers) combined with the large shape anisotropy and stiffness of the microtubules. Diamagnetic anisotropy originates from the anisotropic nature of chemical bonds and is more pronounced in resonance structures such as aromatic groups, peptide bonds, or double and triple carbon bonds .
A single peptide bond has a very weak diamagnetic anisotropy but when many of these bonds are linked together with a fixed and uniform orientation, as in -helices and -sheets, a relatively strong overall anisotropy can result .
The anisotropy of a single tubulin dimer is rather small, and so far the only effects of magnetic fields that have been observed are in concentrated assembled microtubule solutions [1, 3, 6] and microtubule containing structures [7–10]. Even for these systems, high magnetic fields (>8 Tesla) were required, although Vassilev et al.  reported high orientation in a significantly lower magnetic field for an individual microtubule in a diluted solution. An assessment of the value of the tubulin dimer anisotropy would be interesting from a fundamental point of view.
The molecular structure of the tubulin dimer has been elucidated by cryoelectron microscopy , and, in principle, this information could be used to determine the diamagnetic anisotropy by vectorially adding all the diamagnetically anisotropic components. However, assuming the aromatic amino acids have no overall preferred orientation, it should be possible to calculate the value of the diamagnetic anisotropy, to a good approximation, by considering only the contributions from the -helices and -sheets. Validation of this method by a direct cross correlation with experimentally obtained data is not possible since there is no method to measure the diamagnetic anisotropy directly. Therefore an indirect method has to be used.
In order to obtain experimental information about how a dilute solution of dimers responds to an applied magnetic field, Cotton-Mouton experiments can be performed . For axial symmetric molecules and assuming a relatively low degree of orientation (<5% full alignment), the magnetically induced birefringence, , increases linearly with the square of the applied field, . The Cotton-Mouton constant, , is obtained from the slope; thus , where is the wavelength of the laser used in the birefringence experiments. The depends on the number of particles per unit volume and the product of their optical, , and diamagnetic, , anisotropies; thus where is the Boltzman constant and is the absolute temperature of the solution.
depends on the internal structure of the molecule while has potentially two components, one due to shape, form anisotropy, and the other, intrinsic anisotropy, due to internal structure.
Previous attempts to perform Cotton-Mouton experiments on tubulin dimers have been unsuccessful due to the low signal strength. Therefore, an improved protocol had to be designed with special care being taken to optimize the optical components of the equipment and stabilize the temperature of the solution. Experiments were performed on both dimers and double dimers, formed by the addition of Mg2+ ions.
2. Materials and Methods
2.1. Magnetic Birefringence and Cotton-Mouton Experiments
Magnetic birefringence  and Cotton-Mouton experiments, where the birefringence of the solution is measured as function of a steadily increasing magnetic field , were performed in the Grenoble High Magnetic Field Laboratory (Grenoble, France) using the optical equipment mounted on magnet M2. This magnet has a horizontal bore and can reach fields of upward to 17 Tesla . For both experiments, the same optical equipment was used. In the Cotton-Mouton experiments the sample was maintained at constant temperature (4°C) and the magnetic field was ramped up and down at a constant rate. All experiments were repeated a minimum of 6 times.
The effects that we were able to measure were rather weak and therefore care had to be taken to reduce the background contributions to the birefringence signal from the windows of the sample holder and from the solvent. The latter condition considerably reduces the possible choice of buffer solutions since one should avoid all buffers that cause optical Schlieren effects, that is, all buffer solutions with increased viscosity. In these experiments we found that only a tris(hydroxylmethyl)aminomethane (Tris) buffer did not introduce an unacceptable background.
The samples were placed in a thermostated holder and loaded in the magnet after which the magnetic field was ramped to 17 T.
For the Cotton-Mouton experiments analytical ultracentrifugation was performed on the samples in order to determine the particle weight/size distribution. For these experiments, the tubulin is prepared in 20 mM Tris buffer with 0.1 mM guanosine triphosphate (GTP) and at a pH = 7.5. The samples were diluted to 2, 4, and 8 mg/mL. To one aliquot, 4 mM Mg2+ was added in order to induce the formation of oligomers and none to the other.
The sedimentation velocities were determined using a Beckman Optima XL-I ultracentrifuge (Beckman Coulter, Brea, USA) equipped with an interference optical detection system that allowed us to monitor the sedimentation of tubulin at high (0.1 millimolar) concentrations of nucleotide. Samples were studied at a speed of 40,000 rpm and 20°C by using an An50Ti eight-hole rotor and double-sector centerpieces (Beckman Coulter, USA). Differential sedimentation-coefficient distributions, , were calculated by least-squares boundary modelling of sedimentation-velocity data, using the program SEDFIT .
The samples that were prepared without the addition of Mg2+ contained mainly (%) tubulin dimers, that is, the normal unit in which tubulin appears in vitro. A small amount of protein consisted of double or triple dimers (% and %, resp.). This was independent of the protein concentration.
The samples to which Mg2+ was added had a more complex and concentration-dependent composition. At the lower concentration, 30 wt% of the material consisted of aggregates of several dimers. The remaining 70 wt% was composed of dimers and double dimers. At a concentration of 2 mg/mL the sample consisted of equal fractions of dimers and double dimers. The fraction of double dimers increased regularly as function of concentration such that at 8 mg/mL the fraction of material that was not assembled in larger aggregates (90 wt%) was solely composed of double dimers.
3. Results and Discussion
3.1. Cotton-Mouton Experiments on Dimers and Oligomers
In Figure 1, we show the Cotton-Mouton results obtained on a dimer solution of 4 mg/mL. This was the lowest concentration from which reproducible results could be obtained.
The dimer concentrations used were sufficiently dilute so that interdimer interactions could be ignored. This assumption is confirmed by the observation that the Cotton-Mouton constant normalised to the dimer concentration is independent of concentration (see Figure 2).
From these measurements we determined the value of the Cotton-Mouton constant for the dimer normalized to concentration to be . This value is 11 times smaller than that of fibrinogen  but when normalized to molecular weight, , this difference is reduced by a third to 3.7 ( fibrinogen = tubulin dimer). This comparison shows that tubulin dimers are significantly less anisotropic that fibrinogen molecules.
In order to have a cross check on the reliability of the data, oligomers were formed by adding MgCl2 and GTP. Mass determination via ultracentrifuge showed that the oligomers predominantly consist of double-dimers plus single dimers with a small proportion of larger oligomers (see Figure 2). The relative proportion of the different species was found to vary with concentration, and this is probably the reason for the nonlinear increase in the Cotton-Mouton constant as the concentration rises (see Figure 2).
The average Cotton-Mouton constant of the oligomer solution with the largest proportion of double dimers is about twice that obtained from single dimer solutions. This was expected as an approximate doubling of the Cotton-Mouton constant is the maximum that can be expected when dimers are transformed into double dimers. The anisotropy of double dimers depends on the relative orientation of the constituent single dimers and is maximal when the long axes of the single dimers are parallel to each other as in this configuration their anisotropies add together to a good approximation. Any other relative orientation of single dimers would result in lower anisotropy. It is expected that the long axes will indeed be nearly parallel since this is the normal arrangement when the dimers are incorporated in the protofilaments and the assembled microtubules. However, a small deviation angle between the long axes is feasible since the bond between dimers is not completely rigid.
Unfortunately, the magnetic fields available today are nowhere near strong enough to give rise to complete orientation. This limits further data interpretation since without the saturation value in birefringence it is not feasible to determine the optical anisotropy and hence the diamagnetic anisotropy from the CM constant . One has to revert to comparisons with known materials.
3.2. Calculations of the Diamagnetic Moment
By idealizing the tubulin dimer as being axially symmetric, it is possible to obtain an approximate estimate of its diamagnetic anisotropy, , by summing the contributions from the -helices and -sheets. In this we assume that the other potential sources of diamagnetic anisotropy, principally the aromatic amino acids, have no net preferred orientation.
Firstly, each -helix and -pleated sheet was identified, the number of peptide bonds in each, , was totaled, and the angle between the tubulin symmetry axis and the long axis of these secondary structural elements was estimated. The orientation factors with could thus be obtained for each group. From this information was estimated using the following equation adapted from Torbet and Maret : and are the diamagnetic anisotropies of a single peptide bond in either an -helical or a -pleated sheet conformation. Due to the difference in conformation, is only 25% of (3) so -helices are in general expected to make a greater contribution to the total than -pleated sheets and we can simplify (2) to In this way we estimate with -sheets contributing only 15% to the total anisotropy. Unfortunately the value of is not accurately known. Two different estimates for have been published, Pauling  gives while Worcester’s  estimate is 60% greater. The former value is probably more reliable as it is closer to experimental estimates [18, 19] and is consistent with magnetic birefringence measurements made on two filamentous phages . Thus, using Pauling’s value we obtain .
The value of can also be estimated using fibrinogen for calibration as follows. The diamagnetic anisotropy can be calculated using the Cotton-Mouton constant and the birefringence at complete orientation or saturation orientation, , using the relationship  . This equation cannot be exploited with tubulin as we do not have completely aligned samples. However, by putting the published values  for the Cotton-Mouton constant of fibrinogen and the saturation birefringence of fibrin into the latter equation we obtain for the of a single fibrinogen molecule. Again assuming aromatic residues make no net contribution, this anisotropy is due to the axially aligned -helices constituting about 30% of the molecule (i.e., 930 residues). The of fibrinogen is thus equal to as shown above. By comparison with fibrinogen, (i.e., ). This value is 20% larger than that calculated above using the Pauling value for .
The of fibrinogen is 3.7 times larger than that of normalized to molecular weight. As reported above the for fibrinogen is also 3.7 times larger than that of tubulin. This supports our estimate for ; however, this conclusion assumes that the diamagnetic and optical anisotropies are linked by the same proportionality in fibrinogen and tubulin. This is probably approximately true for the intrinsic component of the optical anisotropy because, like the diamagnetic anisotropy, it depends on the anisotropic mobility of electrons in the molecule . But the optical anisotropy can also have a form component which might be relatively different for fibrinogen and tubulin. While the Cotton-Mouton measurements are supportive of our estimate for , they do not constitute conclusive evidence.
can be used to obtain an estimate of the minimum number of tubulin dimers, acting cooperatively, required to attain a highly oriented system. It is known  that for a diamagnetically anisotropic object to attain better than 80% maximum orientation in a magnetic field, , then . In a very strong magnetic field of 10 Tesla at 20°C the minimum number of tubulin dimers required, , is estimated to be in excess of ( daltons) using the lower value for calculated above. As diamagnetic anisotropy is additive to a good approximation, for dimers arranged in parallel. Thus, if a single microtubule is undergoing orientation, in the absence of interaction with its neighbors and assuming a dimer length of 8 nm along the protofilament direction and an average of 13 dimer protofilaments in the tubulin wall, it would need to have a length in excess of , or alternatively a number of smaller microtubules orienting cooperatively could give rise to the same result. It should be pointed out that in concentrated solutions of rigid molecules like microtubules  The magnetic field might not be the only cause of alignment. Above a certain concentration, phase separation and subsequent formation of oriented domains could occur . However, the directors of these domains will still display random orientation. The magnetic field will force these directors to become aligned with the magnetic field. Further research on this has been done but falls outside the scope of this paper.
We have been able to measure the Cotton-Mouton constant of tubulin dimers with a reasonable accuracy. The absence of a reliable value for the optical polarizability prevents a direct determination of the diamagnetic susceptibility but by comparing the Cotton-Mouton constants of the tubulin dimer with fibrinogen, for which the susceptibility is known, we can make a reasonable estimate. This estimate of the susceptibility corresponds well with the calculated contributions of the -helices and -sheets to the diamagnetic susceptibility. If the crystallographic structure is known these, calculations are relatively simple and can be carried out for other proteins as well. The cross correlations that can be made between the experimental results and the simple calculations validate the calculation method for the tubulin dimer and show that this method can be used as an initial assessment of the diamagnetic susceptibility of proteins for which no other data is available. The method that we have used to determine the tubulin dimer diamagnetic susceptibility can be used, for other, proteins and macromolecular assemblies for which this information is not readily available. This can become relevant in for instance medically applied Magnetic Resonance Imaging (MRI) where the applied fields keep increasing.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
The Grenoble High Magnetic Field Laboratory is gratefully acknowledged for allowing the authors access to their facilities. Professor P. Christianen (Nijmegen High Magnetic Field Laboratory) is thanked for his suggestions which the authors feel have improved the paper considerably.
W. Bras, G. P. Diakun, J. F. Díaz et al., “The susceptibility of pure tubulin to high magnetic fields: a magnetic birefringence and X-ray fiber diffraction study,” Biophysical Journal, vol. 74, no. 3, pp. 1509–1521, 1998.View at: Google Scholar
P. M. Vassilev, R. T. Dronzine, M. P. Vassileva, and G. A. Georgiev, “Parallel arrays of microtubules formed in electric and magnetic fields,” Bioscience Reports, vol. 2, no. 12, pp. 1025–1029, 1982.View at: Google Scholar
J. Torbet, “Using magnetic orientation to study structure and assembly,” Trends in Biochemical Sciences C, vol. 12, pp. 327–330, 1987.View at: Google Scholar
G. Maret and K. Dransfeld, Strong and Ultrastrong Magnetic Fields and their Applications, Springer, 1985.
D. L. Worcester, “Structural origins of diamagnetic anisotropy in proteins,” Proceedings of the National Academy of Sciences of the United States of America, vol. 75, no. 11, pp. 5475–5477, 1978.View at: Google Scholar
U. Raviv, D. J. Needleman, K. K. Ewert, and C. R. Safinya, “Hierarchical bionanotubes formed by the self assembly of microtubules with cationic membranes or polypeptides,” Journal of Applied Crystallography, vol. 40, no. 1, pp. s83–s87, 2007.View at: Publisher Site | Google Scholar
J. M. Denegre, J. M. Valles Jr., K. Lin, W. B. Jordan, and K. L. Mowry, “Cleavage planes in frog eggs are altered by strong magnetic fields,” Proceedings of the National Academy of Sciences of the United States of America, vol. 95, no. 25, pp. 14729–14732, 1998.View at: Publisher Site | Google Scholar
R. Emura, N. Ashida, T. Higashi, and T. Takeuchi, “Orientation of bull sperms in static magnetic fields,” Bioelectromagnetics, vol. 22, no. 1, pp. 60–65, 2001.View at: Publisher Site | Google Scholar
J. M. Valles Jr., “Model of magnetic field-induced mitotic apparatus reorientation in frog eggs,” Biophysical Journal, vol. 82, no. 3, pp. 1260–1265, 2002.View at: Google Scholar
K. Guevorkian and J. M. Valles Jr., “Aligning Paramecium caudatum with static magnetic fields,” Biophysical Journal, vol. 90, no. 8, pp. 3004–3011, 2006.View at: Publisher Site | Google Scholar
E. Nogales, M. Whittaker, R. A. Milligan, and K. H. Downing, “High-resolution model of the microtubule,” Cell, vol. 96, no. 1, pp. 79–88, 1999.View at: Publisher Site | Google Scholar
P. G. Vulfson, Molecular Magnetochemistry, Gordon and Breach, 1998.
P. Rub and G. Maret, “New 18-T resistive magnet with radial bores,” IEEE Transactions on Magnetics, vol. 30, no. 4, pp. 2158–2161, 1994.View at: Google Scholar
P. Schuck and P. Rossmanith, “Determination of the sedimentation coefficient distribution by least-squares boundary modeling,” Biopolymers, vol. 54, no. 5, pp. 328–341, 2000.View at: Publisher Site | Google Scholar
J. M. Freyssinet, J. Torbet, G. Hudry Clergeon, and G. Maret, “Fibrinogen and fibrin structure and fibrin formation measured by using magnetic orientation,” Proceedings of the National Academy of Sciences of the United States of America I, vol. 80, no. 6, pp. 1616–1620, 1983.View at: Google Scholar
J. Torbet and G. Maret, “High-field magnetic birefringence study of the structure of rodlike phages Pf1 and fd in solution,” Biopolymers, vol. 20, no. 12, pp. 2657–2669, 1981.View at: Google Scholar
L. Pauling, “Diamagnetic anisotropy of the peptide group,” Proceedings of the National Academy of Sciences of the United States of America, vol. 76, no. 5, pp. 2293–2294, 1979.View at: Google Scholar
K. Tohyama and N. Miyata, “Diamagnetic susceptibility and its anisotropy of poly-γ-ethyl-L- glutamate,” Journal of the Physical Society of Japan, vol. 34, no. 6, article 1699, 1973.View at: Google Scholar
N. S. Murthy, J. R. Knox, and E. T. Samulski, “Order parameter measurements in polypeptide liquid crystals,” The Journal of Chemical Physics, vol. 65, no. 11, pp. 4835–4839, 1976.View at: Google Scholar
R. Oldenbourg and T. Ruiz, “Birefringence of macromolecules. Wiener's theory revisited, with applications to DNA and tobacco mosaic virus,” Biophysical Journal, vol. 56, no. 1, pp. 195–205, 1989.View at: Google Scholar
J. Torbet, “Internal structural anisotropy of spherical viruses studied with magnetic birefringence,” The EMBO Journal, vol. 2, no. 1, pp. 63–66, 1983.View at: Google Scholar
M. G. L. van den Heuvel, S. Bolhuis, and C. Dekker, “Persistence length measurements from stochastic single-microtubule trajectories,” Nano Letters, vol. 7, no. 10, pp. 3138–3144, 2007.View at: Publisher Site | Google Scholar
L. Onsager, “The effects of shape on the interaction of colloidal particles,” Annals of the New York Academy of Sciences: Molecular Interaction, vol. 51, pp. 627–659, 1949.View at: Publisher Site | Google Scholar