How Does Stathmin Destabilize Microtubules? A Root of Consciousness and Alzheimer's disease
Majid Monajjemi1*, Afshar Alihosseini1 and Farnoush Naghsh2
1Department of Chemical engineering, Islamic Azad University, Iran
2Department of Chemistry, Islamic Azad University, Iran
Submission: October 31, 2017; Published: November 16, 2017
*Corresponding author: Majid Monajjemi, Department of Chemical engineering, Islamic Azad University, Iran, Email: m_monajjemi@srbiau.ac.ir
How to cite this article: Majid M, Afshar A, Farnoush N. How Does Stathmin Destabilize Microtubules? A Root of Consciousness and Alzheimer’s disease. Glob J Oto 2017; 11(4): 555819. DOI: 10.19080/GJO.2017.11.555819
Abstract
Alzheimer's disease due to dementia and the memory brain that might be related to the microtubules disorder interrupt daily life. Alzheimer's is a brain disease that causes a slow decline in Consciousness, memory, thinking and reasoning skills. By this research we exhibited the important rules of stathmin destabilize microtubules as a root of this problem in view point of chemical molecular engineering.
Introduction
Microtubule (MTs) is a long protein filament that forms of a dynamic cytoskeletal within a morphological change. MTs have multiple functions in the cellular processes including a peculiar biophysical setting in the internal environment. MTs grow and shrink continually in the living cells qua the harmony between these processes is vital for normal cell function [1-5]. Microtubules are tubular polymers (as its name points)and its diameter is about 24 nm in which the heterodimers are assembled head to tail in a polar fashion (This polarity is reflected by the distinction between the so- called plus and minus ends of proto-filaments).The “building brick" from which they are formed called alpha and betα-tubulins that are proteins and occurs in solution as a dimer of two similar subunits [4,5]. Their assembly is in part determined by the Consciousness of free tubulins in the cytoplasm.
Although most of MTs combinations are largely deprived from the cytoplasmic influence, 2 nm2 lateral pores and 200 nm2 entrances at its ends are connected to cytoplasm. Recently, spherical particles have been found in the microtubules which the particles existence varied between cell types while the neuronal cells have the most particles [6-9]. In 1984, Burton exhibited the particles could be voided from the lumen quickly via reassembly or disassembly of intracellular microtubule [10]. Although identification of components inside of luminal has been difficult, during of the next 30 years observation of these particles in the microtubules lumen has culminated using vitreous cryo- electron microscopy [11]. These subunits polymerize end to end for formation proto-filaments similar to a hollow tubes and their dynamic instability is controlled by numerous compounds. Proteins that destabilize microtubules have been identified recently (such as stathmin family proteins and Colchicines) [12-14].
These proteins increase microtubule turnover in cells, contributing to rapid reorganization of the microtubule cytoskeleton. Stathmin family was initially identified as a protein in response to extracellular signals, over-expressed in highly proliferative breast cancers and malignant ovarian cancers [15,16] . Stathmin is a cellular microtubule inhibitor, which forms a tight complex with two αβ-tubulin dimers, In the references [15,16] the structural information has shown how the stathmin protein family binds tubulin dimers and those studies suggest that phosphorylation occurs in a localized fashion, resulting in decreased microtubule destabilizing activity near microtubule polymer [17]. One of the important stathmin family proteins is RB3 which shares with other members the stathmin-like domain (SLD). Stathmin is a disordered protein, and their activities are down-regulated by multiple phosphorylation. Stathmin is derived from the Greek word “stathmos" and it can also be translated as “terminal or stop", which this translation loosely fits the microtubule destabilizing for the terminating or stopping microtubule growth [18].
Stathmin is a soluble, cytoplasmic protein which enforces an important function in regulating rapid microtubule remaking of the cytoskeleton in response to the cell's requirements. At low Consciousness of free tubulin in the cytoplasm, the growth rate at the microtubule ends is languid and results in an increased rate of disassembly (de-polymerization) [19].The intracellular Consciousness of stathmin varies considerably among different cell types, ranging from 0.005% up to 0.5% of the total cell protein [20], and in vertebrates, is expressed in cells with the potential to proliferate as well as in neurons [20,21]. Proteins related to stathmin are also expressed in the nervous system and include SCLIP, SCG10, RB3 (and two its splice variants RB3'& RB3") [21-23]. Non-polymerized tubulins exists as a hetero dimer of “α-tubulin" & “ β-tubulin" with binding sites for two molecules of guanosine triphosphate (GTP), one exchangeable and the other not. Via electron microscopy (EM) it can be seen that each protofilament consists of globular 4nm subunits and it is possible for refined tubulin to assemble with a range of diameters containing between 9 and16 proto-filaments Scheme 1.
Only two atomic structures of αβ-tubulin are currently known which the first one is GDβ-tubulin in straight anti parallel proto-filaments of two-dimensional zinc sheets stabilized by Taxol (Figure 1) that has been employed to construct pseudo-atomic models of microtubules. The second is the curved structure of two head to tail GDβ-tubulin dimers in complex with the stathmin-like domain of RB3 (RB3-SLD) [15-18]. The structure of tubulin in microtubules is related to the straight zinc-sheet structure, since both are made of similar proto-filaments, although the different lateral contacts cause modifications. Monomers in adjacent proto-filaments arrange a set of shallow helices, which for 13-protofilament microtubule 3 shallow helices run in parallel and for a smaller or larger number of proto-filaments there may be 2 or 4 shallow helices. This difference provides some flexibility in the bonds between adjacent hetero-dimers (at least in the direction running around the microtubule). MT-associated motor proteins, such as dynein and kinesin can be run for long distances along a microtubule when there are 13 proto-filaments (without switching lanes) [24,25]. In more or fewer than 13 proto-filaments, the structure must rotate somewhat, so that the proto-filaments twisted slowly around the microtubule axis [26].
The atomic structure of tubulin proto-filaments is discovered firstly from electron crystallography of zinc-induced 2D sheets [27]. Ravelli et al. [28], extracted important information about conformational change by X-ray crystallography. Each monomer has a pair of spherical domains set which the larger domain, containing the N-terminal half of the polypeptide, has the same fold as a “Rossmann fold" [29]. There is a binding site for the guanosine nucleotide on the plus end area of this domain that contact is made with the next subunit in the proto-filament. A binding site for “Taxol" is placed on the second domain of αβ-tubulin, whichalso makes contact with the core helix, on the opposite side from its contact with the nucleotide base (Figure 1). The C-terminal end of each tubulin polypeptide makes two long helices which those residues (from C-terminal) would be suitable for iso form recognition by tubulin binding proteins.
As it has exhibited in the Figure 2, some long coil loops from the globular domains are complexed in lateral contacts between the proto-filaments in a microtubule. In addition there is common agreement that the “M-loops" of one proto-filament make contact with the guanosine triphosphate (GTPase) domains and the M-loops make different contacts with the adjacent proto-filament (Figure 1). In some positions the proto-filament forms a ring via bending at all of the interfaces between monomers. The bending of protofilament was first seen by electron microscopy [30] and has been confirmed in co crystals of tubulin-sequestering protein stathmin by Gigant et al. [31]. The bending occurs even in the absence of destabilizing agents such as stathmin or colchicines. The core helix is a likely means of communication from the top to the bottom of the β-subunit via a cooperative mechanism. Whole microtubules can twist and bend without coming under great strain or snap and this is clear from images of fluorescently labeled microtubules growing [32].
Although various stathmin activities have been reported, the mechanism by which stathmin affects microtubule dynamics is still a subject of argument. In this work theoretically, we found in vitro that the effects of stathmin can be concluded from its ability to form the strong bonded and non-bonded interaction via a ternary complex including α and β subunit tubulins, αβ hetero dimer of tubulins and one stathmin (Figure 2). When a microtubule bends, individual proto-filaments are bent in a variety of directions. It is likely, therefore, that there are multiple “bent" states for dimers and proto-filaments and there is no certitude that the curved conformations induced by different agents of disassembly are identical [33]. Many efforts have been made to characterize structurally the three major forms of microtubules, hetero dimers, and curved Proto-filament.
The overall structure of the tubulin-stathmin complex has been suggested via scanning transmission electron microscopy combined with digital image processing [34]. Formation of the tubulins-stathmin complex introduced a new path to stabilize tubulin in solution for crystallization. Benoit Gigant et al. [35] has been reported the 4 A° X-ray structure of a complex of GDβ-tubulin with the “Escherichia coli" explicit stathmin-like domain of RB3 (RB3-SLD), a stathmin family protein. This shape defines the 3D structure of the complex, the arrangement of tubulins within this complex, and shows that RB3 contacts tubulin through a 91-residue α-helix. This structure helps to interpreter a clear model for tubulin -stathmin family combination in the non-bonded attraction and provides a structure of GDβ-tubulin which is compared with the structure of tubulin in microtubules.
Stathmin-Microtubules Interactions
Stathmin interacts with two molecules of dimeric α, β-tubulin to form a tight ternary complex which one mole of stathmin binds to two moles of tubulin dimers through the SLD [36]. Tubulin is able to switch between a curved structure in complex with the stathmin- like domain of the RB3 protein and a straight microtubule-like structure. The proto-filament curvature and disassembly has done via GTP hydrolysis. The first opinion for GTP was thought that it would allosterically induce a straight conformation of tubulin subunits thorough-paced of microtubule assembly while the GDP would induce a curved conformation approving disassembly. The microtubule-associated proteins (MAPs) of tubulin-GDP protofilament spirals have indicated different intra- and inter-dimer curvatures [37]. However, this type of allosteric mechanism is challenged by the findings of curved structures of GTP-bound γ-tubulin [38].
This evidence led to a proposition that the free GTβ-tubulin dimer is curved similarly to tubulin rings and is driven into the straight conformation by the microtubule (opposed to what was previously thought). So the GTP γ-phosphate only lowers the unfavorable free energy difference between the curved and the straight form [39]. Consequently, in the allosteric model GTP binding would induce a straighter conformation pre-structured in solution for lateral interactions [40] whereas in the lattice model [39] αβ- tubulin adopts a microtubule incompatible, curved conformation independent of the nucleotide state. It should me important to discuss of the GTP hydrolysis mechanism for destabilizing the microtubule lattice. To provide some answers to this question, Alushin et al., [41] exhibited a structural study for comparing high-resolution cryo-EM reconstructions of GMPCPP microtubules and GDP microtubules. It shows that GTP hydrolysis induces a compression at the linear interface between dimers, immediately over the exchangeable nucleotide-binding site. This compression is amalgamated by conformational changes in α-tubulin. In contrast, lateral contacts between α & β tubulins were basically unchanged in the different nucleotide states.
These understanding suggest that GTP hydrolysis introduces strain into the lattice, but how this strain affects the strength of longitudinal and lateral bonds to destabilize the microtubule remains unknown. In this works 16 main sections have considered for investigation of bonded and non-bonded interaction between Stathmin with microtubules. Although the bonded situations and the GTP hydrolysis essentially are important for the proto-filament curvature and disassembly, the no bonded interaction helps to this processing in other side. More over this kind of non-bonded interaction provides the dynamic behavior for this phenomenon. It now seems clear that changes in the curvature of αβ-tubulin are fundamental to microtubule dynamics and the regulatory activities of microtubule-associated proteins (MAPs) during microtubule polymerization [42].
Results and Discussion
In this work it has been described a novel method of the step- function models for the microtubules within electron density profile in the composition of the stathmin. In this work we have focused on the electron density of the systems when each parts of 16 sections in stathmin has interaction one by one with α and β tubulin (Figures 1-4). The electron density has been defined as
Where η_iis occupation number of orbital i, φ is orbital wave function, is basis function. C is
coefficient matrix, the element of ith row jth column corresponds to the expansion coefficient of
orbital j respect to basis function i.
Atomic unit for electron density can be explicitly written as
e/Bohr3. Bader [43] found that the regions which have large electron localization must have large
magnitudes of Fermi-hole integration. However, the Fermi hole is a six-dimension function and
thus difficult to be studied visually. Becke and Edgecombe noted that spherically averaged likespin conditional pair probability has direct correlation with the Fermi hole and then suggested
electron localization function (ELF) [44].
Savin et al. have reinterpreted ELF in the view of kinetic energy, [45] which makes ELF also meaningful for Kohn-Sham DFT wave- function or even post-HF wave-function. They indicated that D(r) reveals the by Pauli repulsion, while D0(r) can be considered as Thomas-Fermi kinetic energy density. Since D0(r) is introduced into ELF as reference, what the ELF reveals is actually a relative localization. ELF is within the range of [0,1].
A large ELF value means that electrons are greatly localized, indicating that there is a covalent bond, a lone pair or inner shells of the atom involved. ELF has been widely used for a wide variety of systems, such as organic and inorganic small molecules, atomic crystals, coordination compounds, clusters, and for different problems, such as the revealing atomic shell structure, classification of chemical bonding, verification of charge-shift bond, studying aromaticity. Notice that there is a deficiency of ELF, sometimes with r going beyond from molecular boundary, D(r) decreases faster than D0(r) and then ELF reaches 1 (completely localized). To overcome the problem, Multi wfn automatically adds a minimal value 10-5 to D(r). This treatment almost does not affect the ELF value in interesting regions [46]. In which the actual kinetic energy term in D(r) is replaced by Kirzhnits type second-order gradient expansion that is
so that ELF is totally independent from wave-function, and then can be used to analyze electron density from X-ray diffraction data. Of course Tsirelson's ELF can also be used to analyze electron density from quantum chemistry calculation, but is not as good as the ELF defined by Becke owing to the approximation introduced in kinetic energy term; however, qualitative conclusions can still be recovered in general. Localized orbital locator (LOL) is another function for locating high localization regions likewise ELF, defined by Schmider and Becke in the paper [47].
D0 (r) for spin-polarized system and close-shell system are defined in the same way as in ELF (Table 1).
LOL has similar expression compared to ELF. Actually, the chemically significant regions that highlighted by LOL and ELF are generally qualitative comparable, while Jacobsen pointed out that LOL conveys more decisive and clearer picture than ELF, [48]. Obviously LOL can be interpreted in kinetic energy way as for ELF; however LOL can also be interpreted in view of localized orbital. Small (large) LOL value usually appears in boundary (inner) region of localized orbitals because the gradient of orbital wave-function is large (small) in this area. The value range of LOL is identical to ELF, namely [0,1]. The Total electrostatic potential (ESP) measures the electrostatic interaction between a unit point charges placed at r and the system of interest. A positive (negative) value implies that current position is dominated by nuclear (electronic) charges. Molecular electrostatic potential (ESP) has been widely used for prediction of nucleophilic and electrophilic sites for a long time. It is also valuable in studying hydrogen bonds, halogen bonds, molecular recognitions and the intermolecular interaction of aromatics. Moreover, based on statistical analysis, Murray and coworkers found a set of functions called GIPF [49-55], which connects ESP in molecular surface and macroscopic properties (Table 2).
Computational Details
We have simulated a part of microtubule systems including αβ -tubulin hetero dimer and stathmin through QM/MM simulation using Monte Carlo method. Each system was composed of 16 sections of stathmin molecules including interaction with tubulin [56]. Thermodynamic averages for molecular properties were determined from Monte Carlo methods, as can minimum-energy structures [57]. At finite temperature, clusters have finite vapor pressures, and particular cluster sizes are typically unstable to evaporation. Introducing a constraining potential enables one to define clusters of desired sizes. Composed of 16 sections of stathmin molecules were carried out with the simulation microtubules. The pressure was maintained by a variant of the extended system formalism, the Langevin Piston algorithm, which reduces oscillations in the cell parameters. The temperature was maintained at 300 K, well is the body temperature and identical to the relevant experiments.
Configurations of αβ -tubulin heterodimer and stathmin consistent with a mean field were generated by Monte Carlo (MC) simulation, with field values adjusted to obtain agreement with experimental order parameters [58]. In this investigation, differences in force field are illustrated by comparing the calculated energy by using force fields AMBER and OPLS. Furthermore Hyper- Chem professional release 7.01 is used for the calculations. The final parameterization of stathmin was computed using self-consistent field calculations in order to find the optimal starting geometry, as well as the partial charges. We employed density functional theory with the van der Waals density functional to model the exchange- correlation energies of αβ -tubulin heterodimer. All optimization of 16 section of stathmin monomer were performed by Gauessian and GAMESS-US package [59]. We have mainly focused on getting the results from DFT methods such as m062x, m06-L, and m06 for the αβ -tubulin heterodimers. The m062x, m06-L and m06- HF are rather new DFT functional with a good correspondence in non-bonded calculations between tubulin heterodimers and are useful for the energies of distance between two fragments in phospholipids [60].
For non-covalent interactions, the B3LYP method is unable to describe van der Waals [61,62] microtubules systems by medium- range interactions such as the interactions of two tubulins. So the ONIOM methods including 3 levels of 1-high calculation (H), 2-medium calculation (M), and 3-low calculation (L) have been performed in our study for calculating the non-bonded interactions between tubulins. The ab-initio and DFT methods are used for the model system of the ONIOM layers and the semi empirical methods of Pm6 (including pseudo=lanl2) and Pm3MM are used for the medium and low layers, respectively. The semi empirical methods have been used in order to treat the non-bonded interactions between two tubulins.B3LYP and the most other popular functional are insufficient to illustrate the exchange and correlation energy for distant non-bonded medium-range systems correctly. Moreover, some recent studies have shown that inaccuracy for the medium- range exchange energies leads to large systematic errors in the prediction of molecular properties [63,64].
Geometry optimizations and electronic structure calculations have been carried out using the m06 (DFT) functional. This approach is based on an iterative solution of the Kohn-Sham equation [65] of the density functional theory in a plane-wave set with the projector- augmented wave pseudo-potentials. The Perdew-Burke-Ernzerhof (PBE) [66] exchange-correlation (XC) functional of the generalized gradient approximation (GGA) is adopted. The optimizations of the lattice constants and the atomic coordinates are made by the minimization of the total energy. The charge transfer and electrostatic potential-derived charge were also calculated using the Merz-Kollman-Singh [67], chelp [68], or chelpG [69]. The electron density (Both of Gradient norm & Laplacian), value of orbital wave-function, electron spin density, electrostatic potential from nuclear atomic charges, electron localization function (ELF), localized orbital locator (LOL defined by Becke & Tsirelson), total electrostatic potential (ESP), as well as the exchange-correlation density, correlation hole and correlation factor, and the average local ionization energy using the Multifunctional Wave-function Analyzer have also been calculated in this study.
We used Multi wfn software to draw the contour line map. The solid lines indicate positive regions, while the dash lines indicate negative regions. We have plotted the contour line corresponding to vdW surface (electron density=0.001 a.u., which is defined by R. F. W Bader). This is useful to analyze distribution of electrostatic potential on vdW surface. Such a contour line has also been plotted in gradient line and vector field map by the same option. Color, label size and line style of the contour line can be changed based on models. The relief map was used to present the height value at every point. If the values are too large, they will be truncated in the graph. Therefore, it can be chosen to scale the data with a factor to avoid truncation. The graph is shown on interactive interface. Shaded surface map and shaded surface map with projection are used in our representation of height value at each situation.
Conclusion
By this work it has been exhibited the mechanism of microtubules are related to tubulins interaction with stathmin. Each section of 16 parts of stathmin has specific rules for nonbonded interaction between tubulins and stathmin. These subunits polymerize end to end for formation proto -filaments similar to a hollow tubes and their dynamic instability is controlled by numerous compounds. The bending occurs even in the absence of de-stabilizing agents such as stathmin or colchicines. The core helix is a likely means of communication from the top to the bottom of the β-subunit via a cooperative mechanism.
References
- Lowe J, Amos LA (1998) Crystal structure of the bacterial cell-division protein FtsZ. Nature 391(6663): 203-206.
- Nicholas K Gonatas, Elliott Robbins (1965) The homology of spindle and tubules and neuro-tubules in chick embryo retina. Protoplasma 59(3-4): 377-391.
- Clement MJ, Jourdain I, Lachkar S, Savarin P, Gigant B, et al. (2005) N-terminal stathmin-like peptides bind tubulin and impede microtubule assembly. Biochemistry 44(44): 14616-14125.
- Bassot, Martoja (1966) Histological and Ultrastructural data on microtubules cytoplasmic ejaculatory channel insectsorthopteres. Z Zelforsch 74: 145-181.
- Sahu S, Ghosh S, Ghosh B, Aswani K, Hirata K, et al. (2013) Atomic water channel controlling remarkable properties of a single brain microtubule: Correlating single protein to its supramolecular assembly. Biosensors and Bioelectronics 47: 141-148.
- Dusti (1978) Microtubules. Heidelberg, Germany.
- Echandia (1968) Dense core microtubules in neurons and gliocytes of the toad Bufo arenarum Hensel. Am J Anat 122: 157-168.
- Peters (1968) The small pyramidal neuron of the rat cerebral cortex The axon hillock and initil segment. J Cell Biol 39(3): 604-619.
- Stanley (1972) Fine structure of normal spermatic differentiation in Drosophila melanogaster. J Ultrastruc Res 41(5): 433-466.
- Burton (1984) Luminal material in microtubules of frog olfactory axons: Structure and distributin. J Cell Biol 99(2): 520-528.
- Garvalov (2006) Luminal particles within cellular microtubules. J Cell Biol 174(6): 759-765.
- Panda D, Daijo JE, Jordan MA, Wilson L (1995) Kinetic stabilization of microtubule dynamics at steady state in vitro by substoichiometric concentrations of tubulin-colchicine complex. Biochemistry 34: 99219929.
- Belmont LD, Mitchison TJ (1996) Identification of a protein that interacts with tubulin dimers and increases the catastrophe rate of microtubules. Cell 84(4): 623-631.
- Charbaut E (2001) Stathmin family proteins display specific molecular and tubulin binding properties. J Biol Chem 276(19): 16146-16154.
- Hanash SM, Strahler JR, Kuick R, Chu EH, Nichols D (1988) Identificationof a polypeptide associated with the malignant phenotype inacute leukemia. J BiolChem 263: 12813-12815.
- Curmi PA, Nogues C, Lachkar S, Carelle N, Gonthier MP, et al. (2000) Overexpressioin of stathmin in breast carcinomas points out to highly proliferative tumours. Br J Cancer 82(1): 142-150.
- Price DK, Ball JR, Bahrani Mostafavi Z, Vachris JC, Kaufman JS, et al. (2000) The phosphoprotein Op18/stathmin is differentially expressed in ovarian cancer. Cancer Invest 18(8): 722-730.
- Belmont LD, Mitchison TJ (1996) Identification of a protein that interacts with tubulin dimers and increases the catastrophe rate of microtubules. Cell 84: 623-631.
- Cassimeris L (2002) The oncoprotein 18/stathmin family of microtubule destabilizers. Curr Opin Cell Biol 14(1): 18-24.
- Brattsand G, Roos G, Marklund U, Ueda H, Landberg G, et al. (1993) Quantitative analysis of the expression and regulation of an activation- regulated phosphoprotein (oncoprotein18) in normal and neoplastic cells. Leukemia 7: 569-579.
- Rowlands DC, Williams A, Jones NA, Guest SS, Reynolds GM, et al. (1995) Stathmin expression is a feature of proliferating cells of most if not all, cell lineages. Lab Invest 72(1): 100-113.
- Ozon S, Maucuer A, Sobel A (1997) The stathmin family Molecular and biological characterization of novel mammalian proteins expressed in the nervous system. Eur J Biochem 248(3): 794-806.
- Charbaut E, Curmi PA, Ozon S, Lachkar S, Redeker V, et al. (2001) Stathmin family proteins display specific molecular and tubulin binding properties. J BiolChem 276: 16146-16154.
- Hirokawa N (1998) Kinesin and dynein super family proteins and the mechanism of organelle transport. Science 279(5350): 519-526.
- Vale RD (2003) The molecular motor toolbox for intracellular transport. Cell 112(4): 467-480.
- Wade RH, Chretien D, Job D (1990) Characterization of microtubule protofilament numbers How does the surface lattice accommodate? J Mol Biol 212: 775-786.
- Nogales E, Wolf S, Downing KH (1998) Structure of the alpha beta tubulin dimer by electron crystallography. Nature 391: 199-203.
- Ravelli RBG, Gigant B, Curmi PA, Jourdain I, Lachkar S, et al. (2004) Insight into tubulin regulation from a complex with colchicine and a stathmin-like domain. Nature 428(6979): 198-202.
- Nogales E, Downing KH, Amos LA, Lowe J (1998) Tubulin and FtsZ form a distinct family of GTPases. Nat Struct Biol 5: 451-458.
- Watts NR, Cheng N, West W, Steven AC, Sackett DL (2002) The cryptophycin-tubulin ring structure indicates two points of curvature in the tubulin dimer. Biochemistry 41: 12662-12669.
- Gigant B, Curmi PA, Martin Barbey C, Charbaut E, Lachkar S, et al. (2000) The 4 A X-ray structure of a tubulin: stathmin like domain complex. Cell 102(6): 809-816.
- Rodionov V, Nadezhdina E, Borisy G (1999) Centrosomal control of microtubule dynamics. Proc Natl Acad Sci 96(1): 115-120.
- Kraulis PJ (1991) Molscript a program to produce both detailed and schematic plots of protein structures. J Appl Crystallogr 24: 946-950.
- Steinmetz MO, Kammerer RA, Jahnke W, Goldie KN, Lustig A, et al. (2000) Op18/stathmin caps a kinked proto-filament-like tubulin tetramer. EMBO J 19(4): 572-580.
- Benoit Gigant, Patrick A Curmi, Carole Martin Barbey, Elodie Charbaut, Sylvie Lachkar, et al. (2000) The 4 A° X-Ray Structure of a Tubulin: Stathm in like Domain Complex. Cell 102(6): 809-816.
- Sobel A (1991) Stathmin: a relay phosphoprotein for multiple signal transductions? Trends Biochem Sci 16(8): 301-305.
- Wang HW, Nogales E (2005) Nature 435: 911-915.
- Aldaz H, Rice LM, Stearns T, Agard DA (2005) Nature 435: 523-527.
- Buey RM, Di'az JF, Andreu JM (2006) Biochemistry 45: 5933-5938.
- Nogales E, Wang HW (2006) Curr Opin Struct Biol 16: 221-229.
- Alushin GM, GC Lander, EH Kellogg, R Zhang, D Baker, et al. (2014) High- resolution microtubule structures reveal the structural transitions in αβ-tubulin upon GTP hydrolysis. Cell 157(5): 1117-1129.
- Gary J, Brouhard, Luke M Rice J Cell Biol 207(3): 323-334.
- RFW Bader (1990) atoms in Molecule: A quantum Theory.
- Becke, Edgecombe J Chem Phys 92: 5397.
- Savin, Angew. Chem Int Ed Engl 31: 187.
- Tsirelson, Stash Chem Phys Lett 351: 142.
- Schmider, Becke J Mol Struct (THEOCHEM) 527: 51.
- Jacobsen Can J Chem 86: 695.
- Aslangul, Adv Quantum Chem 6: 93.
- Parr J, Phys Chem A 109: 3957.
- Noorizadeh M, Shakerzadeh, Phys. Chem. Chem. Phys 12: 4742.
- Murray, J. Mol. Struct. (THEOCHEM) 307.
- Tieleman DP, Marrink SJ, Berendsen HJ (1997) A computer perspective of membranes: Molecular dynamics studies of lipid bilayer systems. Biochimica et Biophysica Acta 1331(3): 235-270.
- JW Gibbs (1902) Elementary Principles in Statistical Mechanics, Yale University Press, New Haven, CT, USA.
- Feller SE, Pastor RW, Rojnickarin A, Bogusz S, Brooks BR (1996) J Phys Chem 100: 17011.
- Schmidt MW, Baldridge KK, Boatz JA, Elbert ST, Gordon MS (2004) General atomic and molecular electronic structure system. Journal of Computational Chemistry 14(11): 1347-1363.
- Yan Zhao, Donald G Truhlar (2008) The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals. Theoretical Chemistry Accounts 120(1-3): 215-241.
- Zhao Y, Truhlar DG (2008) Density Functionals with Broad Applicability in Chemistry. Accounts of Chemical Research 41(2): 157-167.
- Monajjemi M, Jafari Azan M, Mollaamin F (2013) Density Functional Theory Study on B30N20 Nanocage in Structural Properties and Thermochemical Outlook. Fullerenes Nanotubes and Carbon Nanostructures 21(6): 503-515.
- Grimme S (2006) Seemingly Simple Stereo electronic Effects in Alkane Isomers and the Implications for Kohn-Sham Density Functional Theory. Angewandte Chemie International Edition 45: 4460-4464.
- Monajjemi M, Seyed Hosseini M, Molaamin F (2013) Theoretical Study of Boron Nitride Nanotubes with Armchair Forms. Fullerenes Nanotubes and Carbon Nanostructures 21(5): 381-393.
- Schreiner PR, Fokin AA, Pascal RA, De Meijere A (2006) Many Density Functional Theory Approaches Fail to give Reliable Large Hydrocarbon Isomer Energy Differences. Organic Letters 8: 3635-3638.
- Zhao Y, Truhlar DG (2006) A Density Functional Theory that Accounts for Medium-Range Correlation Energies in Organic Chemistry. Organic Letters 8: 5753-5755.
- Kohn W, Sham LJ (1965) Self-Consistent Equations Including Exchange and Correlation Effects. Phys. Rev 140(4A): 1133-1138.
- Perdew JP, Burke K, Ernzerhof (1996) Generalized Gradient Approximation Made Simple. Physical Review Letters 77(18): 38653868.
- Brent H Besler, Kenneth M Merz, Peter A Kollman (1990) Atomic charges derived from semi empirical methods. Journal of Computational Chemistry 11(4): 431-439.
- Lu T, Chen F (2011) Calculation of Molecular Orbital Composition. Acta Chimica Sinica 69(20): 2393-2406.
- Lu T, Chen F (2012) Quantitative analysis of molecular surface based on improved Marching Tetrahedra algorithm. Journal of Molecular Graphics and Modelling 38: 314-323.
- Lu T, Chen F (2012) Multiwfn: A Multifunctional Wave function Analyzer. Journal of Computational Chemistry 33(5): 580-592.