Conduction pathways in microtubules, biological quantum computation, and consciousness.

Technological computation is entering the quantum realm, focusing attention on biomolecular information processing systems such as proteins, as presaged by the work of Michael Conrad. Protein conformational dynamics and pharmacological evidence suggest that protein conformational states-fundamental information units ('bits') in biological systems-are governed by quantum events, and are thus perhaps akin to quantum bits ('qubits') as utilized in quantum computation. 'Real time' dynamic activities within cells are regulated by the cell cytoskeleton, particularly microtubules (MTs) which are cylindrical lattice polymers of the protein tubulin. Recent evidence shows signaling, communication and conductivity in MTs, and theoretical models have predicted both classical and quantum information processing in MTs. In this paper we show conduction pathways for electron mobility and possible quantum tunneling and superconductivity among aromatic amino acids in tubulins. The pathways within tubulin match helical patterns in the microtubule lattice structure, which lend themselves to topological quantum effects resistant to decoherence. The Penrose-Hameroff 'Orch OR' model of consciousness is reviewed as an example of the possible utility of quantum computation in MTs.     

Molecular modeling study on the differential microtubule-stabilizing effect in singly- and doubly-bonded complexes with peloruside A and paclitaxel

Microtubules (MT) are dynamic cytoskeletal components that play a crucial role in cell division. Disrupting MT dynamics by MT stabilizers is a widely employed strategy to control cell proliferation in cancer therapy. Most MT stabilizers bind to the taxol (TX) site located at the luminal interface between protofilaments, except laulimalide and peloruside A (PLA), which bind to an interfacial pocket on outer MT surface. Cryo-electron microscopy MTs reconstructions have shown differential structural effects on the MT lattice in singly- and doubly-bonded complexes with PLA, TX, and PLA/TX, as PLA is able to revert the lattice heterogeneity induced by TX association leading to more regular MT assemblies. In this work, fully-atomistic molecular dynamics simulations were employed to examine the single and double association of MT stabilizers to reduced MT models in the search for structural and energetic evidence that could be related to the differential regularization and stabilization effects exerted by PLA and TX on the MT lattice. Our results revealed that the double association of PLA/TX (a) strengthens the lateral contact between tubulin dimers compared to singly-bonded complexes, (b) favors a more parallel arrangement between tubulin dimers, and (c) induces a larger restriction in the interdimeric conformational motion increasing the probability of finding structures consistent with 13-protofilaments arrangements. These results and are valuable to increase understanding about the molecular mechanism of action of MT stabilizers, and could account for an overstabilization of MTs in doubly-bonded complexes compared to singly-bonded systems.     

Template-free 13-protofilament microtubule-MAP assembly visualized at 8 A resolution

Microtubule-associated proteins (MAPs) are essential for regulating and organizing cellular microtubules (MTs). However, our mechanistic understanding of MAP function is limited by a lack of detailed structural information. Using cryo-electron microscopy and single particle algorithms, we solved the 8 Å structure of doublecortin (DCX)-stabilized MTs. Because of DCX’s unusual ability to specifically nucleate and stabilize 13-protofilament MTs, our reconstruction provides unprecedented insight into the structure of MTs with an in vivo architecture, and in the absence of a stabilizing drug. DCX specifically recognizes the corner of four tubulin dimers, a binding mode ideally suited to stabilizing both lateral and longitudinal lattice contacts. A striking consequence of this is that DCX does not bind the MT seam. DCX binding on the MT surface indirectly stabilizes conserved tubulin–tubulin lateral contacts in the MT lumen, operating independently of the nucleotide bound to tubulin. DCX’s exquisite binding selectivity uncovers important insights into regulation of cellular MTs.     

Regulation of microtubule dynamics by TOG-domain proteins XMAP215/Dis1 and CLASP

The molecular mechanisms by which microtubule-associated proteins (MAPs) regulate the dynamic properties of microtubules (MTs) are still poorly understood. We review recent advances in our understanding of two conserved families of MAPs, the XMAP215/Dis1 and CLASP family of proteins. In vivo and in vitro studies show that XMAP215 proteins act as microtubule polymerases at MT plus ends to accelerate MT assembly, and CLASP proteins promote MT rescue and suppress MT catastrophe events. These are structurally related proteins that use conserved TOG domains to recruit tubulin dimers to MTs. We discuss models for how these proteins might use these individual tubulin dimers to regulate dynamic behavior of MT plus ends.     

Anisotropic elastic network modeling of entire microtubules

Microtubules are supramolecular structures that make up the cytoskeleton and strongly affect the mechanical properties of the cell. Within the cytoskeleton filaments, the microtubule (MT) exhibits by far the highest bending stiffness. Bending stiffness depends on the mechanical properties and intermolecular interactions of the tubulin dimers (the MT building blocks). Computational molecular modeling has the potential for obtaining quantitative insights into this area. However, to our knowledge, standard molecular modeling techniques, such as molecular dynamics (MD) and normal mode analysis (NMA), are not yet able to simulate large molecular structures like the MTs; in fact, their possibilities are normally limited to much smaller protein complexes. In this work, we developed a multiscale approach by merging the modeling contribution from MD and NMA. In particular, MD simulations were used to refine the molecular conformation and arrangement of the tubulin dimers inside the MT lattice. Subsequently, NMA was used to investigate the vibrational properties of MTs modeled as an elastic network. The coarse-grain model here developed can describe systems of hundreds of interacting tubulin monomers (corresponding to up to 1,000,000 atoms). In particular, we were able to simulate coarse-grain models of entire MTs, with lengths up to 350 nm. A quantitative mechanical investigation was performed; from the bending and stretching modes, we estimated MT macroscopic properties such as bending stiffness, Young modulus, and persistence length, thus allowing a direct comparison with experimental data.     

Information processing mechanisms in microtubules at physiological temperature: Model predictions for experimental tests

Both direct and indirect experimental evidence has shown signaling, communication and conductivity in microtubules (MTs). Theoretical models have predicted that MTs can be potentially used for both classical and quantum information processing although controversies arose in regard to physiological temperature effects on these capabilities. In this paper, MTs have been studied using well-established principles of classical statistical physics as applied to information processing, information storage and signal propagation. To investigate the existence of information processing in MTs we used cellular automata (CA) models with neighbor rules based on the electrostatic properties of the molecular structure of tubulin, and both synchronous and asynchronous updating methods. We obtained a phase diagram of possible dynamic behaviors in MTs that depend on the values of characteristic physical parameters that can be experimentally verified.     

Cytoskeletal signaling: is memory encoded in microtubule lattices by CaMKII phosphorylation?

Memory is attributed to strengthened synaptic connections among particular brain neurons, yet synaptic membrane components are transient, whereas memories can endure. This suggests synaptic information is encoded and 'hard-wired' elsewhere, e.g. at molecular levels within the post-synaptic neuron. In long-term potentiation (LTP), a cellular and molecular model for memory, post-synaptic calcium ion (Ca2?) flux activates the hexagonal Ca2?-calmodulin dependent kinase II (CaMKII), a dodacameric holoenzyme containing 2 hexagonal sets of 6 kinase domains. Each kinase domain can either phosphorylate substrate proteins, or not (i.e. encoding one bit). Thus each set of extended CaMKII kinases can potentially encode synaptic Ca2? information via phosphorylation as ordered arrays of binary 'bits'. Candidate sites for CaMKII phosphorylation-encoded molecular memory include microtubules (MTs), cylindrical organelles whose surfaces represent a regular lattice with a pattern of hexagonal polymers of the protein tubulin. Using molecular mechanics modeling and electrostatic profiling, we find that spatial dimensions and geometry of the extended CaMKII kinase domains precisely match those of MT hexagonal lattices. This suggests sets of six CaMKII kinase domains phosphorylate hexagonal MT lattice neighborhoods collectively, e.g. conveying synaptic information as ordered arrays of six "bits", and thus "bytes", with 64 to 5,281 possible bit states per CaMKII-MT byte. Signaling and encoding in MTs and other cytoskeletal structures offer rapid, robust solid-state information processing which may reflect a general code for MT-based memory and information processing within neurons and other eukaryotic cells.     

Modulating Microtubules: A Molecular Perspective on the Effects of Tail Modifications

Microtubules (MTs), an essential component of the eukaryotic cytoskeleton, are a lattice of polymerized tubulin dimers and are crucial for various cellular processes. The genetic and chemical diversity of tubulin and their disordered tails gives rise to a “tubulin code”. The functional role of tubulin post-translational modifications (PTMs), which contribute to the chemical diversity of the tubulin code, is gradually being unraveled. However, variation in the length and spatial organization of tubulin poly-modifications leads to an enormous combinatorial PTM space, which is difficult to study experimentally. Hence, the impact of the combinatorial tubulin PTM space on the biophysical properties of tubulin tails and their interactions with other proteins remains elusive.Here, we combine all-atom and coarse-grained molecular dynamics simulations to elucidate the biophysical implications of the large combinatorial tubulin PTM space in the context of an MT lattice. We find that tail–body interactions are more dominant in the tubulin dimer than in an MT lattice, and are more significant for the tails of α compared with β tubulin. In addition, polyglutamylation, but not polyglycylation, expands the dimensions of the tubulin tails. Polyglutamylation also leads to a decrease in the diffusion rate of MT-associated protein EB1 on MTs, while polyglycylation often increases diffusion rate. These observations are generally not sensitive to the organization of the polymodifications. The effect of PTMs on MT charge density and tail dynamics are also discussed. Overall, this study presents a molecular quantification of the biophysical properties of tubulin tails and their polymodifications, and provides predictions on the functional importance of tubulin PTMs.     

Tubulin Dimers Oligomerize before Their Incorporation into Microtubules

In the presence of GTP, purified dimers of α- and β-tubulin will interact longitudinally and laterally to self-assemble into microtubules (MTs). This property provides a powerful in vitro experimental system to describe MT dynamic behavior at the micrometer scale and to study effects and functioning of a large variety of microtubule associated proteins (MAPs). Despite the plethora of such data produced, the molecular mechanisms of MT assembly remain disputed. Electron microscopy (EM) studies suggested that tubulin dimers interact longitudinally to form short oligomers which form a tube by lateral interaction and which contribute to MT elongation. This idea is however challenged: Based on estimated association constants it was proposed that single dimers represent the major fraction of free tubulin. This view was recently supported by measurements suggesting that MTs elongate by addition of single tubulin dimers. To solve this discrepancy, we performed a direct measurement of the longitudinal interaction energy for tubulin dimers. We quantified the size distribution of tubulin oligomers using EM and fluorescence correlation spectroscopy (FCS). From the distribution we derived the longitudinal interaction energy in the presence of GDP and the non-hydrolysable GTP analog GMPCPP. Our data suggest that MT elongation and nucleation involves interactions of short tubulin oligomers rather than dimers. Our approach provides a solid experimental framework to better understand the role of MAPs in MT nucleation and growth.     

Arp2/3 phosphorylation kickstarts cells

The kinesin-13 motor protein family members drive the removal of tubulin from microtubules (MTs) to promote MT turnover. A point mutation of the kinesin-13 family member mitotic centromere-associated kinesin/Kif2C (E491A) isolates the tubulin-removal conformation of the motor, and appears distinct from all previously described kinesin-13 conformations derived from nucleotide analogues. The E491A mutant removes tubulin dimers from stabilized MTs stoichiometrically in adenosine triphosphate (ATP) but is unable to efficiently release from detached tubulin dimers to recycle catalytically. Only in adenosine diphosphate (ADP) can the mutant catalytically remove tubulin dimers from stabilized MTs because the affinity of the mutant for detached tubulin dimers in ADP is low relative to lattice-bound tubulin. Thus, the motor can regenerate for further cycles of disassembly. Using the mutant, we show that release of tubulin by kinesin-13 motors occurs at the transition state for ATP hydrolysis, which illustrates a significant divergence in their coupling to ATP turnover relative to motile kinesins.     

Motor-dependent microtubule disassembly driven by tubulin tyrosination

In cells, stable microtubules (MTs) are covalently modified by a carboxypeptidase, which removes the C-terminal Tyr residue of α-tubulin. The significance of this selective detyrosination of MTs is not understood. In this study, we report that tubulin detyrosination in fibroblasts inhibits MT disassembly. This inhibition is relieved by overexpression of the depolymerizing motor mitotic centromere-associated kinesin (MCAK). Conversely, suppression of MCAK expression prevents disassembly of normal tyrosinated MTs in fibroblasts. Detyrosination of MTs suppresses the activity of MCAK in vitro, apparently as the result of a decreased affinity of the adenosine diphosphate (ADP)–inorganic phosphate- and ADP-bound forms of MCAK for the MT lattice. Detyrosination also impairs MT disassembly in neurons and inhibits the activity of the neuronal depolymerizing motor KIF2A in vitro. These results indicate that MT depolymerizing motors are directly inhibited by the detyrosination of tubulin, resulting in the stabilization of cellular MTs. Detyrosination of transiently stabilized MTs may give rise to persistent subpopulations of disassembly-resistant polymers to sustain subcellular cytoskeletal differentiation.     

Structural studies of N-terminal mutants of Connexin 32 using (1)H NMR spectroscopy

Microtubules (MTs) control cell replication, material transport and motion in eukaryotic cells, but MT role in several pathologies is still unknown. These functions are related to the MT physico-chemical properties and MT formation mode starting from tubulin molecules. This study describes a new method, based on the computer aided analysis of the electron paramagnetic resonance (EPR) spectra of selected spin probes to obtain structural and dynamical information on tubulins and MTs and the kinetics of MTs formation promoted by guanosine-5'-triphosphate (GTP). It was found that tubulin and MTs avoid radical quenching caused by ethylene glycol tetraacetic acid (EGTA). MT formation showed different kinetics as a function of tubulin concentration. At 5 mg/mL of tubulin, MTs were formed in 8 min. These results are also useful for getting information on MT-drug interactions.     

Intrinsic bending of microtubule protofilaments

The complex polymerization dynamics of the microtubule (MT) plus end are closely linked to the hydrolysis of the GTP nucleotide bound to the β-tubulin. The destabilization is thought to be associated with the conformational change of the tubulin dimers from the straight conformation in the MT lattice to a curved conformation. It remains under debate whether this transformation is directly related to the nucleotide state, or a consequence of the longitudinal or lateral contacts in the MT lattice. Here, we present large-scale atomistic simulations of short tubulin protofilaments with both nucleotide states, starting from both extreme conformations. Our simulations indicate that both interdimer and intradimer contacts in both GDP and GTP-bound tubulin dimers and protofilaments in solution bend. There are no observable differences between the mesoscopic properties of the contacts in GTP and GDP-bound tubulin or the intradime and interdimer interfaces.     

Relating nucleotide‐dependent conformational changes in free tubulin dimers to tubulin assembly

The complex dynamic behavior of microtubules (MTs) is believed to be primarily due to the αβ‐tubulin dimer architecture and its intrinsic GTPase activity. Hence, a detailed knowledge of the conformational variations of isolated α‐GTP‐β‐GTP‐ and α‐GTP‐β‐GDP‐tubulin dimers in solution and their implications to interdimer interactions and stability is directly relevant to understand the MT dynamics. An attempt has been made here by combining molecular dynamics (MD) simulations and protein–protein docking studies that unravels key structural features of tubulin dimer in different nucleotide states and correlates their association to tubulin assembly. Results from simulations suggest that tubulin dimers and oligomers attain curved conformations in both GTP and GDP states. Results also indicate that the tubulin C‐terminal domain and the nucleotide state are closely linked. Protein–protein docking in combination with MD simulations suggest that the GTP‐tubulin dimers engage in relatively stronger interdimer interactions even though the interdimer interfaces are bent in both GTP and GDP tubulin complexes, providing valuable insights on in vitro finding that GTP‐tubulin is a better assembly candidate than GDP‐tubulin during the MT nucleation and elongation processes. © 2012 Wiley Periodicals, Inc. Biopolymers 99: 282–291, 2013.     

Quantum Information Processing in the Wall of Cytoskeletal Microtubules

Microtubules (MT) are composed of 13 protofilaments, each of which is a series of two-state tubulin dimers. In the MT wall, these dimers can be pictured as “lattice” sites similar to crystal lattices. Based on the pseudo-spin model, two different location states of the mobile electron in each dimer are proposed. Accordingly, the MT wall is described as an anisotropic two-dimensional (2D) pseudo-spin system considering a periodic triangular “lattice”. Because three different “spin-spin” interactions in each cell exist periodically in the whole MT wall, the system may be shown to be an array of three types of two-pseudo-spin-state dimers. For the above-mentioned condition, the processing of quantum information is presented by using the scheme developed by Lloyd.     

Newly assembled microtubules are concentrated in the proximal and distal regions of growing axons

We have investigated the sites of microtubule (MT) assembly in neurons during axon growth by taking advantage of the relationship between the proportion of tyrosinated alpha-tubulin (tyr-tubulin) in MTs and their age. Specifically, young (newly assembled) MTs contain more tyr-tubulin than older (more long-lived) MTs. To quantify the relative proportion of tyr-tubulin in MTs, cultured rat sympathetic neurons were permeabilized under conditions that stabilize existing MTs and remove unassembled tubulin. The MTs were then double-stained with antibodies to tyr-tubulin (as a measure of the amount of tyr-tubulin in MTs) and to beta-tubulin (as a measure of total MT mass), using immunofluorescence procedures. Cells were imaged with a cooled charge-coupled device camera and the relative proportion of tyr-tubulin in the MTs was quantified by computing the ratio of the tyr-tubulin fluorescence to the beta-tubulin fluorescence using a novel application of digital image processing and analysis techniques. The amount of tyr-tubulin in the MTs was highest in the cell body and at the growth cone; peak ratios in these two regions were approximately 10-fold higher than for the axon shaft. Moving out from the cell body into the axon, the tyr-tubulin content declined over an average distance of 40 microns to reach a constant low value within the axon shaft and then rose again more distally, over an average distance of 110 microns, to reach a peak at the growth cone (average axon length = 358 microns). These observations indicate that newly assembled MTs are concentrated in the proximal and distal regions of growing axons and therefore that the cell body and growth cone are the most active sites of MT assembly dynamics in neurons that are actively extending axons.     

Quantitative analysis of microtubule transport in growing nerve processes

In neurons, tubulin is synthesized primarily in the cell body, whereas the molecular machinery for neurite extension and elaboration of microtubule (MT) array is localized to the growth cone region. This unique functional and biochemical compartmentalization of neuronal cells requires transport mechanisms for the delivery of newly synthesized tubulin and other cytoplasmic components from the cell body to the growing axon. According to the polymer transport model, tubulin is transported along the axon as a polymer. Because the majority of axonal MTs are stationary at any given moment, it has been assumed that only a small fraction of MTs translocates along the axon by saltatory movement reminiscent of the fast axonal transport. Such intermittent "stop and go" MT transport has been difficult to detect or to exclude by using direct video microscopy methods. In this study, we measured the translocation of MT plus ends in the axonal shaft by expressing GFP-EB1 in Xenopus embryo neurons in culture. Formal quantitative analysis of MT assembly/disassembly indicated that none of the MTs in the axonal shaft were rapidly transported. Our results suggest that transport of axonal MTs is not required for delivery of newly synthesized tubulin to the growing nerve processes.     

Fission yeast Alp14 is a dose-dependent plus end-tracking microtubule polymerase

XMAP215/Dis1 proteins are conserved tubulin-binding TOG-domain proteins that regulate microtubule (MT) plus-end dynamics. Here we show that Alp14, a XMAP215 orthologue in fission yeast, Schizosaccharomyces pombe, has properties of a MT polymerase. In vivo, Alp14 localizes to growing MT plus ends in a manner independent of Mal3 (EB1). alp14-null mutants display short interphase MTs with twofold slower assembly rate and frequent pauses. Alp14 is a homodimer that binds a single tubulin dimer. In vitro, purified Alp14 molecules track growing MT plus ends and accelerate MT assembly threefold. TOG-domain mutants demonstrate that tubulin binding is critical for function and plus end localization. Overexpression of Alp14 or only its TOG domains causes complete MT loss in vivo, and high Alp14 concentration inhibits MT assembly in vitro. These inhibitory effects may arise from Alp14 sequestration of tubulin and effects on the MT. Our studies suggest that Alp14 regulates the polymerization state of tubulin by cycling between a tubulin dimer-bound cytoplasmic state and a MT polymerase state that promotes rapid MT assembly.