This Article Abstract Full Text (PDF) All Versions of this Article: biophysj.106.089987v1 92/2/394    most recent Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Drabik, P. Articles by Kovalenko, A. Search for Related Content PubMed PubMed Citation Articles by Drabik, P. Articles by Kovalenko, A.

Microtubule Stability Studied by Three-Dimensional Molecular Theory of Solvation

Piotr Drabik ^*, Sergey Gusarov ^* and Andriy Kovalenko ^* 

^* National Institute for Nanotechnology, National Research Council of Canada, Edmonton, Alberta, Canada; and Department of Mechanical Engineering, University of Alberta, Edmonton, Alberta, Canada

Correspondence: Address reprint requests to Andriy Kovalenko, National Institute for Nanotechnology, 11421 Saskatchewan Dr., Edmonton, AB T6G 2M9, Canada. E-mail: andriy.kovalenko{at}nrc-cnrc.gc.ca.

	ABSTRACT

TOP ABSTRACT INTRODUCTION METHODS MOLECULAR THEORY OF SOLVATION... RESULTS AND DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
We study microtubular supramolecular architectures of tubulin dimers self-assembling into linear protofilaments, in turn forming a closed tube, which is an important component of the cytoskeleton. We identify the protofilament arrangements with the lowest free energy using molecular dynamics to optimize tubulin conformations. We then use the three-dimensional molecular theory of solvation to obtain the hydration structure of protofilaments built of optimized tubulins and the solvent-mediated effective potential between them. The latter theoretical method, based on first principles of statistical mechanics, is capable of predicting the structure and thermodynamics of solvation of supramolecular architectures. We obtained a set of profiles of the potential of mean force between protofilaments in a periodic two-dimensional sheet in aqueous solution. The profiles were calculated for a number of amino acid sequences, tubulin conformations, and spatial arrangements of protofilaments. The results indicate that the effective interaction between protofilaments in aqueous solution depends little on the isotypes studied; however, it strongly depends on the M loop conformation of ß-tubulin. Based on the analysis of the potential of mean force between adjacent protofilaments, we found the optimal arrangement of protofilaments, which is in good agreement with other studies. We also decomposed the potential of mean force into its energetic and entropic components, and found that both are considerable in the free-energy balance for the stabilized protofilament arrangements.

	INTRODUCTION

TOP ABSTRACT INTRODUCTION METHODS MOLECULAR THEORY OF SOLVATION... RESULTS AND DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
All living cells utilize a complicated dense jungle of protein polymers known as the cytoskeleton. This is a three-dimensional fine network of three kinds of filament structures (actin filaments, intermediate filaments, and microtubules), thoroughly coordinated and subtly structurally and functionally linked, and performing as "bones and muscles" of the cell. Microtubules, polymers of tubulin dimers, are considered to be involved in various tasks or cellular features, including cell morphology, intracellular transport (e.g., endoplasmic reticulum-Golgi transport), centralization of nucleus, chromosome segregation during cell division, chromosome motility after DNA damage, cell stiffness control, pathogen infections, memory, and consciousness (1–12). Responsible for different complex tasks, microtubules undergo equally complex posttranslational modifications and unique interactions with multiple families of satellite proteins (13–20). They are very dynamical systems, growing and shrinking in a guanosine-triphosphate (GTP)-hydrolysis-dependent manner (21–26).

The dynamic nature of microtubules makes them susceptible to pharmacological agents. Compounds that perturb microtubule dynamics are currently some of the most effective drugs to treat medical conditions (27–29). Tubulin-assembly modulators represent an important class of antitumor drugs, as they have been proven to be an effective tool for cancer chemotherapy (30,31). Numerous tubulin ligands with antimitotic properties and anticancer potential have been discovered recently (28,32,33). Apart from pharmacology, microtubules can be used in bionanotechnology to form nanowires (when plated with metals) and scaffolds for nanofabrication. Another potential nanotechnology application involves the molecular transport machinery of microtubules and kinesin proteins that can be integrated in kinesin-powered microdevices (34–36). This opens up new possibilities for the development of guided nanotransport mechanisms in nanodevices.

Microtubules are composed of dimers of protein tubulin. The heterodimer consists of one -tubulin and one ß-tubulin subunit (each 40 Å in diameter), which are homologous. The dimers are aligned head-to-tail in rows called protofilaments. A variable number of protofilaments form a hollow tube, that is, a microtubule. The external and internal diameters of the microtubule are 240 and 150 Å; in vitro, the length can reach up to 10–50 µm. Polymerization is a polar process that reflects the polarity of the microtubule. Tubulin polymerizes more quickly from the plus end, which is terminated by the ß-subunit. The other end, growing more slowly, is known as the minus end, and is terminated by the -subunit.

In vivo, microtubule cylinders usually have 13 protofilaments, though the number may be different in particular situations. In vitro, it is possible for purified tubulin to assemble with a fairly wide range of diameters containing 9–18 protofilaments. This variation reveals that there is some flexibility in the bonds between adjacent protofilaments, at least in the direction involved in the curvature of the microtubule wall. In vivo, the occurrence of microtubules with a number of protofilaments other than 13 appears to be determined by isoforms with specific changes in the amino acid sequences of - and ß-tubulin, as, for example, in the case of some specialized 15-protofilament microtubules in neurons of the nematode Caenorhabditis elegans (37,38).

There are three published structures of pure tubulin dimers, i.e., those not complexed with other proteins: Protein Data Bank (pdb) code 1TUB, described by Nogales and Downing (39); pdb code 1JFF, by Lowe et al. (40); and pdb code 1TVK, by Nettles et al. (41). All of these are of relatively poor quality, with the best resolution, 2.9 Å, for the 1TVK structure. Tubulin, an 55-kDa protein, has an /ß fold where a core of two ß-sheets is surrounded by -helices. Each monomer is composed of three functional domains, the nucleotide-binding (containing the Rossman fold), drug-binding, and motor-protein-binding domains (Fig. 1).

View larger version (66K): [in this window] [in a new window]   FIGURE 1  Ribbon diagram of the conformation of tubulin dimer and the schematic orientation of tubulin dimers in two adjacent protofilaments. Vertical white lines denote the long axis of the microtubule. The position of the M loop is indicated by arrows. The figure was generated with PyMOL (91).

The tubulin polymerization process involves two types of contacts between tubulin subunits: head-to-tail binding of dimers, resulting in protofilaments that run along the length of the microtubule; and lateral interactions between parallel protofilaments, which complete the microtubule wall. Lateral interaction surfaces are more electrostatic and less hydrophobic than the longitudinal contacts. The longitudinal contacts along protofilaments appear to be much stronger than those between adjacent protofilaments (43). Structural studies have established the critically important role of two tubulin loops, the M loop (S7-H9 loop) and the N loop (H1–S2 loop), in microtubule self-assembly. These loops coordinate lateral interactions between protofilaments to build a microtubule. Studies with docking the tubulin dimer structure into low-resolution images of microtubules have established that M loops interact with N loops of laterally adjacent subunits (44). Moreover, organisms with increased microtubule stability (such as arctic fish) have two amino acid substitutions in the M loop (45). The "stickiness" of the M-loop side of the protofilament appears to be largely responsible for the polymorphic nature of tubulin polymers (46).

Microtubules constitute one of the most intriguing questions of modern biophysics. Although much progress has been made, there are still many unresolved issues concerning their molecular structure, architecture, and dynamic behavior. Due to their interesting and important features, microtubules are attracting increased attention from researchers (42–44,47–52). Various aspects of the structure, dynamics, self-assembly, and stability of microtubules have been studied using theoretical and computational approaches. However, most of these studies included no atomic representation of microtubules (7,8,10,53–67). Limited exceptions include an all-atom study of 90 dimers, extremely CPU-intense (700 processors with an aggregate speed of 1 TFlops), where electrostatics was treated by applying the Poisson-Boltzmann equation (68). Another study published by the same group addressed a very small portion of a microtubule (five dimers only) by a combination of the Poisson-Boltzmann electrostatics and the surface-area term (69). The apolar contribution was obtained by assuming that each Ų of buried solvent-accessible area contributed 11 cal/mole to the binding energy. The authors found that the B-lattice is the most favorable configuration of protofilaments and that lateral bonds are significantly weaker than longitudinal bonds along protofilaments.

In this study, we treat sets of protofilaments in the all-atom representation using a method of integral equation theory of molecular liquids and solutions, namely, the three-dimensional reference interaction site model (3D-RISM) (70). In the search for factors determining the microtubule architecture for a sequence of different tubulin isotypes and conformations of the crucial ß-tubulin M loop at the protein-protein interface, the 3D-RISM theory yields potentials of mean force between adjacent protofilaments as a function of the relative orientation of protofilaments and of the vertical and horizontal distances between them.

	METHODS

TOP ABSTRACT INTRODUCTION METHODS MOLECULAR THEORY OF SOLVATION... RESULTS AND DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
The coordinates of the heavy atoms for the bovine tubulin dimer were obtained from the Protein Data Bank, entry 1TVK, as determined by Nettles et al. (41) via electron crystallography. All the cofactors visible in the structure (epothilone, GTP/GDP) were removed for the sake of simplicity. The missing loop (residues 35–60) not visible in the -tubulin monomer was built using the coordinates from the ß-tubulin structure as a guide, with subsequent energy minimization. Hydrogen atoms were added in the AMBER package (71). As the resolution of the experimental structure was relatively low (2.9 Å), the original tubulin dimer was subjected to molecular dynamics for 3 ns using the AMBER force field with the generalized Born model of implicit solvation (72). During molecular dynamics, all hydrogen-containing bonds were constrained using the SHAKE algorithm (73).

The final snapshot from the well equilibrated trajectory of the original bovine dimer was mutated (via in silico substitutions) to generate five different ß-tubulin isotypes (ß_I–ß_V). When building different tubulin isotypes, we decided to keep the monomer unchanged and vary only the ß-monomer sequence, as the variability of sequence is observed mostly in ß-tubulin. A guide for mutations was performed by sequence alignment in the Clustal W program accessible on-line (74). Clustal W is a general-purpose, multiple-sequence alignment program for DNA or proteins. It produces biologically meaningful multiple-sequence alignments of divergent sequences. It calculates the best match for the selected sequences, and lines them up so that identities, similarities, and differences can be seen. Each isotype dimer was again equilibrated in 1 ns of molecular dynamics in the AMBER force field with implicit solvent.

As a next step, a series of four different conformers was generated for the M loop in ß_V-tubulin, as this loop is regarded as the main player in the interprotofilament interactions. The ß_V isotype was chosen as a random example, with no particular reason to give preference to any other isotype. In high-temperature simulation, the M-loop residues were allowed to move while the rest of the protein was frozen. Tubulin was heated in constrained simulated annealing from 0 to 1000 K in 1 ps, and snapshots were then generated every 5 ps in a 40-ps thermal stabilization. Of the resulting eight snapshots, the four most structurally divergent ones were selected by visual inspection (snapshots 1, 3, 6, and 7), cooled down for 50 ps, and then used for further calculations of the potentials of mean force (see below). In the simulated annealing, the following general constraints were used: positional constraints for all C atoms in tubulin, except for M-loop residues; and angular constraints for all peptide bonds and all improper dihedrals to retain proper chirality.

Using the periodic boundary conditions of a supercell in the 3D-RISM calculation (see below), we built two-dimensional sheets of protofilaments. The protofilaments were incrementally translated along the long axis to vary the offset of dimers (Fig. 2). We started with the zero-offset, a configuration in which the dimers would form a ring in the microtubule, and slid the protofilaments in 2-Å increments until a second ring conformation was obtained at 80 Å offset in direction h. Simultaneously, the distance between the long axes of protofilaments (direction r) was varied in the range 54–63 Å, with a step of 1 Å; the minimal distance between protofilaments with no severe van der Waals collisions between adjacent strands is 54 Å. As a result, we obtained 10 x 40 nodes, giving 400 cases of spatial arrangement (Fig. 2). Counting five different isotypes, with and four conformers for one of them, this totaled 3200 values of the potential of mean force. (We chose to present only the most representative results, i.e., the sections passing through the minima in both the h and r directions).

View larger version (35K): [in this window] [in a new window]   FIGURE 2  Schematic representation of spatial arrangements between two adjacent protofilaments, with the r and h directions indicated as discussed in the text (-tubulins and ß-tubulins are represented as white and dark ribbons, respectively). The figure was generated with PyMOL (91).

	MOLECULAR THEORY OF SOLVATION FOR EFFECTIVE POTENTIALS BETWEEN PROTOFILAMENTS

TOP ABSTRACT INTRODUCTION METHODS MOLECULAR THEORY OF SOLVATION... RESULTS AND DISCUSSION ACKNOWLEDGEMENTS REFERENCES

To obtain the structure and thermodynamics of hydration of protofilaments, we employed the 3D-RISM integral equation complemented with the Kovalenko-Hirata (KH) closure approximation (70,77). The 3D-RISM-KH integral equation theory of molecular solvation explicitly accounts for the effect of chemical specificities of molecular and ionic species in electrolyte solution on the solvation structure and thermodynamics, in particular, association effects and formation of supramolecules. For instance, the 3D-RISM-KH theory has been proven to be efficient and uniquely capable of treating and predicting self-assembly and conformational stability of such complex organic supramolecular architectures as synthetic organic rosette nanotubes in solution (78). This method provides a detailed microscopic insight into the organization of solvent molecules in the solvation shell and their role in thermodynamic stability and optimal arrangement of the fragments that self-assemble into a supramolecule in solution.

The theory yields the solvation structure in the form of a three-dimensional distribution of normalized probability density, the three-dimensional distribution function of solvent site , showing site density enhancement (g > 1) or depletion (g > 1) around the solute macromolecule relative to the average site number density in the solution bulk. The 3D-RISM integral equation is written as (70,77,79,80)

(1)

where is the three-dimensional total correlation function of site related to the three-dimensional distribution function as , and is the three-dimensional direct correlation function which has the asymptotics of the solute-solvent site interaction potential: ; is the site-site susceptibility of pure solvent; superscripts u and v denote solute and solvent, respectively; and the site subscript indices and enumerate all sites on all sorts of solvent molecules. The convolution in Eq. 1 is calculated by using the three-dimensional Fast Fourier Transform technique, with special analytical treatment of the electrostatic asymptotics of all the correlation functions (70,82). The other relation between the three-dimensional total and direct correlation functions, complementing the 3D-RISM integral equation (Eq. 1) is the 3D-KH closure approximation (70,82),

(2)

where is the three-dimensional intermolecular interaction potential between the whole solute and solvent site , which is specified by the molecular force field, and k_BT is the Boltzmann constant times the solution temperature. This relation combines, in a nontrivial manner, the advantages of two known closures: the mean spherical approximation (75) is applied to high association peaks and long-range electrostatic tails for and the hypernetted chain approximation (75) is used in the core repulsion and other depletion regions for , with the switching occurring smoothly, at , or equally, at .

The site-site susceptibility of solvent breaks up into the intra- and intermolecular terms,

(3)

where the intramolecular correlation function, or intramolecular matrix, represents the geometry of solvent molecules with the site-site separations l (z-matrix) specified by the molecular force field, and is the radial total correlation function between sites and in bulk solvent. In advance of the 3D-RISM-KH calculation, the site-site correlation functions of bulk solvent are obtained from the dielectrically consistent RISM theory (81) coupled with the KH closure (DRISM-KH) (70,82). The bulk solvent susceptibility (Eq. 3) is then input into the 3D-RISM integral equation (Eq. 1).

For a given arrangement of the solute supramolecule immersed in multicomponent solvent, the solvation free energy is obtained from the 3D-RISM-KH integral equations (Eqs. 1 and 2) in a closed analytical form in terms of the three-dimensional correlation functions (70,77). The potential of mean force (PMF) between protofilaments can then be obtained as a difference between the solvation free energy of the protofilaments in the supramolecule and those pulled apart (70,82). However, this approach is less practical, as it requires calculation of the three-dimensional correlation functions and the solvation free energy of the supramolecule for a sequence of protofilament arrangements to map the PMF landscape. Instead, in this work we use the definition of the PMF between two molecules in solution in terms of the 6D distribution function,

(4)

where r₁₂ and ₂ are the position and orientation of solute molecule 2 with respect to the coordinate system fixed at solute molecule 1. (The orientation of molecule 1 is thus fixed, too.) The distribution function between molecules 1 and 2 at a given orientation ₂ can be conveniently obtained from the 3D-RISM-HNC integral equation for solute-solute correlations at infinite dilution (70,83). The PMF (Eq. 4) is thus obtained in the form

(5)

where is the direct interaction potential between molecules 1 and 2 as a function of coordinates r₁₂ and ₂, and the latter term gives the solvent-mediated part of the PMF. The three-dimensional direct and total correlation functions and of solvent site around solute molecules u₁ and u₂ are obtained from the 3D-RISM-KH integral equations (Eqs. 1 and 2). The three-dimensional correlation function at a particular orientation ₂ of solute molecule 2 is nothing but rotated by angle ₂, and can be obtained by rendering. Equation 5 yields the PMF between protofilaments 1 and 2 and all their images with the supercell periodicity, specified on a three-dimensional grid of their relative coordinates r₁₂ at a given orientation ₂ of protofilament 2. Notice that the 3D-RISM equations operate with a system at constant density, and so the PMF obtained corresponds to the NVT ensemble.

To obtain the PMF for each system, the 3D-RISM-KH equations were solved on a three-dimensional grid in a rectangular supercell of the following sizes set to accommodate the two-dimensional periodic sheet of protofilaments: 1) the current value of separation r between protofilaments in the sheet; 2) 128 Å (commensurate with the distance between tubulin dimers in a protofilament) in direction h; and 3) the same value 128 Å. The grid resolution of 0.5 Å we used was fine enough to obtain the results without significant numerical errors. The SPC/E water model (84) was used, and the water solvent was at ambient thermodynamic conditions of T = 298 K and density = 0.997 g/cm³. The Amber parameter set (amber99) was used for the interaction site charges and Lennard-Jones parameters of proteins (85).

	RESULTS AND DISCUSSION

TOP ABSTRACT INTRODUCTION METHODS MOLECULAR THEORY OF SOLVATION... RESULTS AND DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
The results of sequence alignment for ß-tubulin are shown in Fig. 3. It is very clear that there is almost no variability in the M-loop sequence (residues 279–287), whereas most of the differences among the isotypes are clustered in the C-terminal region. This region plays a major role in regulating microtubule assembly (86). The majority of differences between the known isotypes localize within the last 15 residues of the sequences. The PMF between protofilaments as a function of the offset h is presented in Fig. 4. Only the ß_V isotype is presented, as all the isotypes display almost identical plots. It might be surprising that there are no differences in the PMF among isotypes. The main difference between the isotype sequences is within the C-terminus, where the homology is the lowest. These C-termini tails are exposed on the outside surface of the microtubule and are very flexible in molecular dynamics runs. The difference in the C-terminus sequence could possibly affect the dynamical process of microtubule assembly or translate into different spatial arrangements of protofilaments in a microtubular cylinder. However, their conformations do not affect the PMFs between protofilaments arranged in a 2D sheet. It can be seen clearly that there is one sharp global minimum at approximately h = 6 Å. It corresponds very well with the value found by Sept et al. (69) in their CPU-time-consuming approach to a single system. This global minimum corresponds to the so-called B-lattice. It is generally established that a B-lattice is the predominant arrangement of tubulin monomers in the microtubule lattice (87), where lateral contacts are made between homologous subunits, that is, - and ß-ß. Our results confirm that.

View larger version (34K): [in this window] [in a new window]   FIGURE 3  Results of the sequence alignment performed in Clustal W for five ß-tubulin isotypes. An asterisk (*) means that the residues or nucleotides in the column are identical in all sequences in the alignment; a colon (:) means that conserved substitutions have been observed; and a period (.) means that semiconserved substitutions are observed.

View larger version (7K): [in this window] [in a new window]   FIGURE 4  Profile of the potential of mean force (PMF) between adjacent protofilaments of ßV-tubulin, conformer 1, as a function of their offset h at separation r = 54 Å.

View larger version (11K): [in this window] [in a new window]   FIGURE 5  Global minimum region of the potential of mean force (PMF) between adjacent protofilaments as a function of the offset h at separation r = 54 Å. All conformers of ßV-tubulin.

View larger version (7K): [in this window] [in a new window]   FIGURE 6  Potential of mean force (PMF) between adjacent protofilaments of ßV-tubulin, conformer 1, as a function of the separation r at a the offset h = 6 Å in the B-lattice.

View larger version (10K): [in this window] [in a new window]   FIGURE 7  Global minimum region of the potential of mean force (PMF) between adjacent protofilaments as a function of the separation r at the offset h = 6 Å in the B-lattice. All conformers of ßV-tubulin.

(6)

In this calculation, we chose the temperature difference value to be T = 10 K. It should be noted that the intramolecular entropy of tubulin is not included in the decomposition (6) giving the thermodynamics of solvation only.

We performed the decomposition (6) for one selected system built of ß_v-tubulin, conformer 1. Fig. 8 shows the PMF decomposition into the energetic and entropic terms E(h) and –TS(h) as functions of the offset h between the protofilaments. The entropic contribution clearly forms two minima corresponding to lattices A and B. It is interesting to note that the second minimum is pronounced much more clearly than in the full PMF profile. Investigating the lateral interactions, Sept et al. (69) decomposed the free energy into the electrostatic and apolar contributions. The latter term can be seen as hydrophobic, and was calculated assuming that each Ų of solvent-accessible surface area (SASA) that is buried contributes 11 cal/mol to the binding energy. The hydrophobic effect is generally known to be entropically rather than enthalpically driven, which gives a rationale for comparing this apolar term to our entropic contribution to the PMF. Both the studies localize two distinct minima, and both locate them at approximately the same offset h. However, somewhat contrary to our results, Sept et al. (69) see both minima as equally deep, arguing that because the two tubulin monomers are structurally homologous, the A- and B-lattices have similar amounts of buried surface area. As distinct, our results show that the minimum associated with the B-lattice is the main one, being deeper by 80 kcal/mol than the second one for the A-lattice. It is difficult to speculate about the exact nature of the observed difference, but one obvious explanation might be that we calculate a pure entropic term, whereas Sept at al. (69) deal with an "entropy-like" term within the empiric SASA approach. Another explanation could be that the PMF decomposition from 3D-RISM is simply inaccurate. Some problems with the 1D-RISM-HNC theory for hydrophobic interactions were reported previously (90). In that study, the enthalpy and entropy were found to make comparable contributions to the attraction of nonpolar solutes in water, and that observation disagreed with the conventional description of an entropy-driven hydrophobic interaction (90). However, in our case of the 3D-RISM-KH theory, we actually do not observe comparable contributions of energetic and entropic parts, i.e., the enntropic term is about twice as large in magnitude and opposite in sign. Nevertheless, the changes in these terms at the positions around the second minimum are quite comparable, implying that they are both important in formation of the second minimum. Our results indicate that the entropic contribution around the offset h corresponding to the A-lattice has a well defined minimum that is strongly suppressed in the PMF by the maximum of the energetic term (see Fig. 8). This makes the second minimum much shallower and wider, unlike the arrangement corresponding to the B-lattice, where both of them have a well pronounced minimum.

View larger version (10K): [in this window] [in a new window]   FIGURE 8  Decomposition of the potential of mean force (PMF) between protofilaments into the energetic and entropic contributions, W(h) = E(h) – TS(h), as functions of the offset h at the separation r = 54 Å. ßv-tubulin, conformer 1.

View larger version (8K): [in this window] [in a new window]   FIGURE 9  Decomposition of the potential of mean force (PMF) between the protofilaments into the energetic and entropic contributions, W(r) = E(r) – TS(r), as functions of the separation r at offset h = 6 Å corresponding to the B-lattice. ßv-tubulin, conformer 1.

To summarize, the 3D-RISM method proved to be useful and capable of reproducing the thermodynamically preferable self-organizing supramolecular architecture, which has been suggested by experimental results and obtained in computationally intense molecular simulations. With reasonable computational effort much smaller than molecular simulations the 3D-RISM method allows one to predict, from the first principles, the three-dimensional solvation structure and solvation thermodynamics of a supramolecule in molecular solvent. This gives access to the PMFs between protofilaments in solution, yielding predictions for the conformational stability and optimal arrangement of protofilaments in a supramolecule. These characteristics are very difficult, and mostly not feasible, to obtain from molecular simulations. Being cost-effective, the 3D-RISM method allows one to scan a large number of systems.

We can conclude that our results for the potential of mean force between adjacent protofilaments in aqueous solution indicate that the interactions between protofilaments are independent of isotype studied. However, they strongly depend on the M-loop conformation of ß-tubulin, suggesting that the microtubular architecture may be more sensitive to the M-loop conformation than to the isotype sequence. These data were obtained for the first time that we know of. Based on the analysis of the potential of mean force between adjacent protofilaments, we also found the optimal arrangement of protofilaments, which is in good agreement with other theoretical studies, as well as experiment.

	ACKNOWLEDGEMENTS

TOP ABSTRACT INTRODUCTION METHODS MOLECULAR THEORY OF SOLVATION... RESULTS AND DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
We are grateful to Prof. Jack Tuszynski for stimulating discussions. We thank Mr. Tyler Luchko for preliminary MD simulations for tubulin conformation.

We gratefully acknowledge support from the Alberta Prion Research Institute, Prion Project G224200003, and from the National Research Council of Canada. P.D. and S.G. were National Science and Engineering Research Council (NSERC) Visiting Fellows. The computations were supported by the Centre of Excellence in Integrated Nanotools at the University of Alberta.

Submitted on May 24, 2006; accepted for publication September 28, 2006.

	REFERENCES

TOP ABSTRACT INTRODUCTION METHODS MOLECULAR THEORY OF SOLVATION... RESULTS AND DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
1. Vaughan, K. T. 2005. Microtubule plus ends, motors, and traffic of Golgi membranes. Biochim. Biophys. Acta. 1744:316–324.[Medline]

2. Becker, B. E., and L. Cassimeris. 2005. Cytoskeleton: microtubules born on the run. Curr. Biol. 15:R551–R554.[CrossRef][Medline]

3. Maly, I. V., and G. G. Borisy. 2002. Self-organization of treadmilling microtubules into a polar array. Trends Cell Biol. 12:462–465.[CrossRef][Medline]

4. Henry, T., J.-P. Gorvel, and S. Meresse. 2006. Molecular motors hijacking by intracellular pathogens. Cell. Microbiol. 8:23–32.[CrossRef][Medline]

5. Stamenovic, D. 2005. Microtubules may harden or soften cells, depending on the extent of cell distension. J. Biomech. 38:1728–1732.[CrossRef][Medline]

6. Crawford-Young, S. J. 2006. Effects of microgravity on cell cytoskeleton and embryogenesis. Int. J. Dev. Biol. 50:183–191.[CrossRef][Medline]

7. Portet, S., J. A. Tuszynski, J. M. Dixon, and M. V. Sataric. 2003. Models of spatial and orientational self-organization of microtubules under the influence of gravitational fields. Phys. Rev. E. 68:021903.[CrossRef]

8. Faber, J., R. Portugala, and L. P. Rosa. 2006. Information processing in brain microtubules. Biosystems. 83:1–9.[Medline]

9. Hameroff, S. R., and R. Penrose. 1998. Conscious events as orchestrated space time selections. In Explaining Consciousness: The Hard Problem. J. Shear, editor. MIT Press, Cambridge, MA.

10. Tuszynski, J. A., B. Trpisova, D. Sept, and M. V. Sataric. 1997. The enigma of microtubules and their self-organization behavior in the cytoskeleton. Biosystems. 42:153–175.[CrossRef][Medline]

11. Tuszynski, J. A., J. A. Brown, and P. Hawrylak. 1998. Dielectric polarization, electric conduction, information processing and quantum computation in microtubules. Are they plausible? The Royal Society. 356:1897–1926.

12. Tuszynski, J. A., S. R. Hameroff, M. V. Sataric, B. Trpisova, and M. L. A. Nip. 1995. Ferroelectric behavior in microtubule dipole lattices: implications for information processing, signaling and assembly/disassembly. J. Theor. Biol. 174:371–380.[CrossRef]

13. Bloom, K. 2004. Microtubule composition: Cryptography of dynamic polymers. Proc. Natl. Acad. Sci. USA. 101:6839–6840.[Free Full Text]

14. Luduena, R. F. 1998. The multiple forms of tubulin: different gene products and covalent modifications. Int. Rev. Cytol. 178:207–275.[Medline]

15. Bertrand, S., I. Barthelemy, M. A. Oliva, J. L. Carrascosa, J. M. Andreu, and J. M. Valpuesta. 2005. Folding, stability and polymerization properties of FtsZ chimeras with inserted tubulin loops involved in the interaction with the cytosolic chaperonin CCT and in microtubule formation. J. Mol. Biol. 346:319–330.[CrossRef][Medline]

16. Lewis, S. A., G. Tian, I. Vainberg, and N. J. Cowan. 1996. Chaperonin-mediated folding of actin and tubulin. J. Cell. Biol. 132:1–4.[Free Full Text]

17. Szymanski, D. 2002. Tubulin folding cofactors: half a dozen for a dimmer. Curr. Biol. 12:R767–R769.[CrossRef][Medline]

18. Lewis, S. A., G. Tian, and N. J. Cowan. 1997. The - and ß-tubulin folding pathways. Trends Cell Biol. 7:479–484.[CrossRef][Medline]

19. Surrey, T., F. Nedelec, S. Leibler, and E. Karsenti. 2001. Physical properties determining self-organization of motors and microtubules. Science. 292:1167–1171.[Abstract/Free Full Text]

20. Downing, K. H. 2000. Structural basis for the interaction of tubulin with proteins and drugs that affect microtubule dynamics. Annu. Rev. Cell Dev. Biol. 16:89–111.[CrossRef][Medline]

21. Carvalho, P., J. S. Tirnauer, and D. Pellman. 2003. Surfing on microtubule ends. Trends Cell Biol. 13:229–237.[CrossRef][Medline]

22. Dammermann, A., A. Desai, and K. Oegema. 2003. The minus end in sight. Curr. Biol. 13:R614–R624.[CrossRef][Medline]

23. Howard, J., and A. A. Hyman. 2003. Dynamics and mechanics of the microtubule plus end. Nature. 422:753–758.[CrossRef][Medline]

24. Krebs, A., K. N. Goldie, and A. Hoenger. 2005. Structural rearrangements in tubulin following microtubule formation. EMBO Rep. 6:227–232.[CrossRef][Medline]

25. Nogales, E., and H.-W. Wang. 2006. Structural intermediates in microtubule assembly and disassembly: how and why? Curr. Opin. Cell Biol. 18:1–6.[CrossRef]

26. Wang, H.-W., and E. Nogales. 2005. Nucleotide-dependent bending flexibility of tubulin regulates microtubule assembly. Nature. 435:911–915.[CrossRef][Medline]

27. Honore, S., E. Pasquier, and D. Braguer. 2005. Understanding microtubule dynamics for improved cancer therapy. Cell. Mol. Life Sci. 62:3039–3056.[CrossRef][Medline]

28. Jordan, A., J. A. Hadfield, N. J. Lawrence, and A. T. McGown. 1998. Tubulin as a target for anticancer drugs: agents which interact with the mitotic spindle. Med. Res. Rev. 18:259–296.[CrossRef][Medline]

29. Morrissette, N. S., A. Mitra, D. Sept, and L. D. Sibley. 2004. Dinitroanilines bind -tubulin to disrupt microtubules. Mol. Biol. Cell. 4:1960–1968.

30. Buey, R. M., I. Barasoain, E. Jackson, A. Meyer, P. Giannakakou, I. Paterson, S. Mooberry, J. M. Andreu, and J. F. Diaz. 2005. Microtubule interactions with chemically diverse stabilizing agents: thermodynamics of binding to the paclitaxel site predicts cytotoxicity. Chem. Biol. 12:1269–1279.[CrossRef][Medline]

31. Wilson, L., and M. A. Jordan. 2004. New microtubule/tubulin-targeted anticancer drugs and novel chemotherapeutic strategies. J. Chemother. 16:83–85.

32. Hamel, E. 1996. Antimitotic natural products and their interactions with tubulin. Med. Res. Rev. 16:207–231.[CrossRef][Medline]

33. Shi, Q., K. Chen, S. L. Morris-Natschke, and K. H. Lee. 1998. Recent progress in the development of tubulin inhibitors as antimitotic antitumor agents. Curr. Pharm. Des. 4:219–248.[Medline]

34. van den Heuvel, M. G. L., C. T. Butcher, S. G. Lemay, S. Diez, and C. Dekker. 2005. Electrical docking of microtubules for kinesin-driven motility in nanostructures. Nano Lett. 5:235–241.[CrossRef][Medline]

35. Limberis, L., and R. J. Stewart. 2000. Toward kinesis-powered microdevices. Nanotechnology. 11:1–5.

36. Romet-Lemonne, G., M. VanDuijn, and M. Dogterom. 2005. Three-dimensional control of protein patterning in microfabricated devices. Nano Lett. 5:2350–2354.[CrossRef][Medline]

37. Fukushige, T., Z. K. Siddiqui, M. Chou, J. G. Culotti, C. B. Gogonea, S. S. Siddiqui, and M. Hamelin. 1999. MEC-12, an -tubulin required for touch sensitivity in C. elegans. J. Cell Sci. 112:395–403.[Abstract]

38. Savage, C., M. Hamelin, J. G. Culotti, A. Coulson, D. G. Albertson, and M. Chalfie. 1989. Mec-7 is a ß-tubulin gene required for the production of 15-protofilament microtubules in Caenorhabditis elegans. Genes Dev. 3:870–881.[Abstract/Free Full Text]

39. Nogales, E., and K.A. Downing. 1998. Tubulin and microtubule structure. Curr. Opin. Cell. Biol. 10:16–22.[CrossRef][Medline]

40. Lowe, J., H. Li, K. H. Downing, and E. Nogales. 2001. Refined structure of ß-tubulin at 3.5 Å resolution. J. Mol. Biol. 313:1045–1057.[CrossRef][Medline]

41. Nettles, J. H., H. Li, B. Cornett, J. M. Krahn, J. P. Snyder, and K. H. Downing. 2004. The binding mode of epothilone A on ,ß-tubulin by electron crystallography. Science. 305:866–869.[Abstract/Free Full Text]

42. Nogales, E. 2001. Structural insights into microtubule function. Annu. Rev. Biophys. Biomol. Struct. 30:397–420.[CrossRef][Medline]

43. Nogales, E., M. Whittaker, R. A. Milligan, and K. H. Downing. 1999. High-resolution model of the microtubule. Cell. 96:79–88.[CrossRef][Medline]

44. Li, H., D. DeRosier, W. Nicholson, E. Nogales, and K. H. Downing. 2002. Microtubule structure at 8 Å resolution. Structure. 10:1317–1328.[Medline]

45. Detrich 3rd, H. W., S. K. Parker, R. C. Williams, Jr., E. Nogales, and K. H. Downing. 2000. Cold adaptation of microtubule assembly and dynamics. Structural interpretation of primary sequence changes present in the - and ß-tubulins of Antarctic fishes. J. Biol. Chem. 275:37038–37047.[Abstract/Free Full Text]

46. Amos, L. A., and J. Lowe. 1999. How Taxol stabilises microtubule structure. Chem. Biol. 6:R65–R69.[CrossRef][Medline]

47. Amos, L. A. 2004. Microtubule structure and its stabilisation. Org. Biomol. Chem. 2:2153–2160.[CrossRef][Medline]

48. Insinna, E. M., P. Zaborski, and J. A. Tuszynski. 1996. Electrodynamics of microtubular motors: the building blocks of a new model. Biosystems. 39:187–226.[CrossRef][Medline]

49. Meurer-Grob, P., J. Kasparian, and R. H. Wade. 2001. Microtubule structure at improved resolution. Biochemistry. 40:8000–8008.[CrossRef][Medline]

50. Nogales, E., and H.-W. Wang. 2006. Structural mechanisms underlying nucleotide-dependent self-assembly of tubulin and its relatives. Curr. Opin. Struct. Biol. 16:221–229.[CrossRef][Medline]

51. Tuszynski, J. A., J. A. Brown, and D. Sept. 2003. Models of the collective behavior of proteins in cells: tubulin, actin and motor proteins. J. Biol. Phys. 29:401–428.[CrossRef]

52. Ravelli, R. B. G., B. Gigant, P. A. Curmi, I. Jourdain, S. Lachkar, A. Sobel, and M. Knossow. 2004. Insight into tubulin regulation from a complex with colchicines and a stathmin-like domain. Nature. 428:198–202.[CrossRef][Medline]

53. Bolterauer, H., H. J. Limbach, and J. A. Tuszynski. 1999. Microtubules: strange polymers inside the cell. Bioelectrochem. Bioenerg. 48:285–295.[CrossRef][Medline]

54. Cytrynbaum, E. N., V. Rodionov, and A. Mogilner. 2004. Computational model of dynein-dependent self-organization of microtubule asters. J. Cell Sci. 117:1381–1397.[Abstract/Free Full Text]

55. Hunyadi, V., D. Chretien, and I. M. Janosi. 2005. Mechanical stress induced mechanism of microtubule catastrophes. J. Mol. Biol. 348:927–938.[CrossRef][Medline]

56. Janulevicius, A., J. van Pelt, and A. van Ooyen. 2006. Compartment volume influences microtubule dynamic instability: a model study. Biophys. J. 90:788–798.[Abstract/Free Full Text]

57. Martin, S. R., M. J. Schilstra, and P. M. Bayley. 1993. Dynamic instability of microtubules: Monte Carlo simulation and application to different types of microtubule lattice. Biophys. J. 65:578–596.[Abstract/Free Full Text]

58. Mershin, A., A. A. Kolomenski, H. A. Schuessler, and D. V. Nanopoulos. 2004. Tubulin dipole moment, dielectric constant and quantum behavior: computer simulations, experimental results and suggestions. Biosystems. 77:73–85.[CrossRef][Medline]

59. Mishra, P. K., A. Kunwar, S. Mukherji, and D. Chowdhury. 2005. Dynamic instability of microtubules: effect of catastrophe-suppressing drugs. Phys. Rev. E. 72:051914.[CrossRef]

60. Mitra, A., and D. Sept. 2004. Localization of the antimitotic peptide and depsipeptide binding site on ß-tubulin. Biochemistry. 43:13955–13962.[CrossRef][Medline]

61. Portet, S., J. A. Tuszynski, C. W. V. Hogue, and J. M. Dixon. 2005. Elastic vibrations in seamless microtubules. Eur. Biophys. J. 34:912–920.[CrossRef][Medline]

62. Sataric, M. V., and J. A. Tuszynski. 2003. Relationship between the nonlinear ferroelectric and liquid crystal models for microtubules. Phys. Rev. E. 67:011901.[CrossRef]

63. Son, J., G. Orkoulas, and A. B. Kolomeisky. 2005. Monte Carlo simulations of rigid biopolymer growth processes. J. Chem. Phys. 123:124902.[CrossRef][Medline]

64. Tuszynski, J. A., T. Luchko, S. Portet, and J. M. Dixon. 2005. Anisotropic elastic properties of microtubules. Eur. Phys. J. E. 17:29–35.

65. Tuszynski, J. A., E. J. Carpenter, J. T. Huzil, W. Malinski, T. Luchko, and R. F. Luduena. 2006. The evolution of the structure of tubulin and its potential consequences for the role and function of microtubules in cells and embryos. Int. J. Dev. Biol. 50:341–358.[CrossRef][Medline]

66. VanBuren, V., D. J. Odde, and L. Cassimeris. 2002. Estimates of lateral and longitudinal bond energies within the microtubule lattice. Proc. Natl. Acad. Sci. USA. 99:6035–6040.[Abstract/Free Full Text]

67. VanBuren, V., L. Cassimeris, and D. J. Odde. 2005. Mechanochemical model of microtubule structure and self-assembly kinetics. Biophys. J. 89:2911–2926.[Abstract/Free Full Text]

68. Baker, N. A., D. Sept, S. Joseph, M. J. Holst, and J. A. McCammon. 2001. Electrostatics of nanosystems: application to microtubules and the ribosome. Proc. Natl. Acad. Sci. USA. 98:10037–10041.[Abstract/Free Full Text]

69. Sept, D., N. A. Baker, and J. A. McCammon. 2003. The physical basis of microtubule structure and stability. Protein Sci. 12:2257–2261.[Abstract/Free Full Text]

70. Kovalenko, A. 2003. Three-dimensional RISM theory for molecular liquids and solid-liquid interfaces. In Molecular Theory of Solvation. Understanding Chemical Reactivity. F. Hirata and P. G. Mezey, editors. Kluwer Academic Publishers, Dordrecht, The Netherlands. 169–275.

71. Pearlman, D. A., D. A. Case, J. W. Caldwell, W. S. Ross, T. E. Cheatham III, S. DeBolt, D. Ferguson, G. Seibel, and P. A. Kollman. 1995. AMBER: a package of computer programs for applying molecular mechanics, normal mode analysis, molecular dynamics and free energy calculations to simulate the structural and energetic properties of molecules. Comput. Phys. Commun. 91:1–41.

72. Still, W. C., A. Tempczyk, R. C. Hawley, and T. Hendrickson. 1990. Semianalytical treatment of solvation for molecular mechanics and dynamics. J. Am. Chem. Soc. 112:6127–6129.[CrossRef]

73. Ryckaert, J.-P., G. Ciccotti, and H. J. C. Berendsen. 1977. Numerical integration of the Cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanes. J. Comp. Phys. 23:327–341.[CrossRef]

74. Higgins, D., J. Thompson, T. Gibson, J. D. Thompson, D. G. Higgins, and T. J. Gibson. 1994. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 22:4673–4680.[Abstract/Free Full Text]

75. Hansen, J. P., and I. R. McDonald. 1986. Theory of Simple Liquids, 2nd Ed. Academic Press, London, UK.

76. Hirata, F., editor. 2003. Molecular Theory of Solvation, Series, Vol. 24: Understanding Chemical Reactivity. P.G. Mezey, editor. Kluwer Academic, Dordrecht, The Netherlands.

77. Kovalenko, A., and F. Hirata. 1999. Self-consistent description of a metal-water interface by the Kohn-Sham density functional theory and the three-dimensional reference interaction site model. J. Chem. Phys. 110:10095–10112.[CrossRef]

78. Moralez, J. G., J. Raez, T. Yamazaki, R. K. Motkuri, A. Kovalenko, and H. Fenniri. 2005. Helical rosette nanotubes with tunable stability and hierarchy. J. Am. Chem. Soc. 127:8307–8309.[CrossRef]

79. Beglov, D., and B. J. Roux. 1997. An integral equation to describe the solvation of polar molecules in liquid water. J. Phys. Chem. B. 101:7821–7826.[CrossRef]

80. Kovalenko, A., and F. Hirata. 1998. Three-dimensional density profiles of water in contact with a solute of arbitrary shape: a RISM approach. Chem. Phys. Lett. 290:237–244.[CrossRef]

81. Perkyns, J. S., and B. M. Pettitt. 1992. A site-site theory for finite concentration saline solutions. J. Chem. Phys. 97:7656–7666.[CrossRef]

82. Kovalenko, A., and F. Hirata. 2000. Potentials of mean force of simple ions in ambient aqueous solution. I. Three-dimensional reference interaction site model approach. J. Chem. Phys. 112:10391–10402.[CrossRef]

83. Kovalenko, A., and F. Hirata. 1999. Potential of mean force between two molecular ions in a polar molecular solvent: a study by the three-dimensional reference interaction site model. J. Phys. Chem. B. 103:7942–7957.[CrossRef]

84. Berendsen, H. J. C., J. R. Grigera, and T. P. Straatsma. 1987. The missing term in effective pair potentials. J. Phys. Chem. 91:6269–6271.[CrossRef]

85. Wang, J., P. Cieplak, and P. A. Kollman. 2000. How well does a restrained electrostatic potential (RESP) model perform in calculating conformational energies of organic and biological molecules? J. Comput. Chem. 21:1049–1074.[CrossRef]

86. Luduena, R. F. 1993. Are tubulin isotypes functionally significant. Mol. Biol. Cell. 4:445–457.[Medline]

87. Song, Y.-H., and E. Mandelkow. 1993. Recombinant kinesin motor domain binds to ß-tubulin and decorates microtubules with a B surface lattice. Proc. Natl. Acad. Sci. USA. 99:1671–1675.

88. Pettitt, B. M., and P. J. Rossky. 1986. Alkali halides in water: ion-solvent correlations and ion-ion potentials of mean force at infinite dilution. J. Chem. Phys. 84:5836–5844.[CrossRef]

89. Choudhury, N., and B. M. Pettitt. 2006. Enthalpy-entropy contributions to the potential of mean force of nanoscopic hydrophobic solutes. J. Phys. Chem. B. 110:8459–8463.[Medline]

90. Yu, H.-A., B. Roux, and M. Karplus. 1990. Solvation thermodynamics: an approach from analytic temperature derivatives. J. Chem. Phys. 92:5020–5033.[CrossRef]

91. Delano, W. L. 2002. The PyMOL Molecular Graphics System. DeLano Scientific, San Carlos, CA. http://www.pymol.org/.

This article has been cited by other articles: