Microtubules are significant therapeutic targets for the treatment of cancer, where suppression of microtubule dynamicity by drugs such as paclitaxel forms the basis of clinical efficacy. Peloruside A, a macrolide isolated from New Zealand marine sponge Mycale hentscheli, is a microtubule-stabilizing agent that synergizes with taxoid drugs through a unique site and is an attractive lead compound in the development of combination therapies. We report here unique allosteric properties of microtubule stabilization via peloruside A and present a structural model of the peloruside-binding site. Using a strategy involving comparative hydrogen-deuterium exchange mass spectrometry of different microtubule-stabilizing agents, we suggest that taxoid-site ligands epothilone A and docetaxel stabilize microtubules primarily through improved longitudinal interactions centered on the interdimer interface, with no observable contributions from lateral interactions between protofilaments. The mode by which peloruside A achieves microtubule stabilization also involves the interdimer interface, but includes contributions from the alpha/beta-tubulin intradimer interface and protofilament contacts, both in the form of destabilizations. Using data-directed molecular docking simulations, we propose that peloruside A binds within a pocket on the exterior of beta-tubulin at a previously unknown ligand site, rather than on alpha-tubulin as suggested in earlier studies.

RASSF1A is a tumor suppressor protein involved in death receptor-dependent apoptosis utilizing the Bax-interacting protein MOAP-1 (previously referred to as MAP-1). However, the dynamics of death receptor recruitment of RASSF1A and MOAP-1 are still not understood. We have now detailed recruitment to death receptors (tumor necrosis factor receptor 1 [TNF-R1] and TRAIL-R1/DR4) and identified domains of RASSF1A and MOAP-1 that are required for death receptor interaction. Upon TNF-␣ stimulation, the C-terminal region of MOAP-1 associated with the death domain of TNF-R1; subsequently, RASSF1A was recruited to MOAP-1/TNF-R1 complexes. Prior to recruitment to TNF-R1/MOAP-1 complexes, RASSF1A homodimerization was lost. RASSF1A associated with the TNF-R1/MOAP-1 or TRAIL-R1/MOAP-1 complex via its N-terminal cysteine-rich (C1) domain containing a potential zinc finger binding motif. Importantly, TNF-R1 association domains on both MOAP-1 and RASSF1A were essential for death receptor-dependent apoptosis. The association of RASSF1A and MOAP-1 with death receptors involves an ordered recruitment to receptor complexes to promote cell death and inhibit tumor formation.

We present a new methodology for computing solvation free energy, which is based upon the reference interaction site model (RISM)/hypernetted chain (HNC) solvation free energy expression, but which substitutes radial distribution functions taken from simulations for those calculated by simultaneous solution of the RISM and HNC equations. Consequently, solvation free energy can be obtained from a single molecular dynamics or Monte Carlo simulation. Here we describe in detail the coupled RISM/simulation approach, and offer some error analysis. Finally we give the results of its application to a set of small test molecules in aqueous solution. The success shown in some of our results demonstrates that the coupled RISM/simulation approach is worth considering further as a potentially useful tool in studies of solvated systems, such as aqueous molecular biosystems.

No abstract

The cytoskeleton is essential to cell morphology, cargo trafficking, and cell division. As the neuronal cytoskeleton is extremely complex, it is no wonder that a startling number of neurodegenerative disorders (including but not limited to Alzheimer’s disease, Parkinson’s disease and Huntington’s disease) share the common feature of a dysfunctional neuronal cytoskeleton. Recently, concern has been raised about a possible link between anesthesia, post-operative cognitive dysfunction, and the exacerbation of neurodegenerative disorders. Experimental investigations suggest that anesthetics bind to and affect cytoskeletal microtubules, and that anesthesia-related cognitive dysfunction involves microtubule instability, hyper-phosphorylation of the microtubule-associated protein tau, and tau separation from microtubules. However, exact mechanisms are yet to be identified. In this paper the interaction of anesthetics with the microtubule subunit protein tubulin is investigated using computer-modeling methods. Homology modeling, molecular dynamics simulations and surface geometry techniques were used to determine putative binding sites for volatile anesthetics on tubulin. This was followed by free energy based docking calculations for halothane (2-bromo-2-chloro-1,1,1-trifluoroethane) on the tubulin body, and C-terminal regions for specific tubulin isotypes. Locations of the putative binding sites, halothane binding energies and the relation to cytoskeleton function are reported in this paper.

Memristors represent the fourth electrical circuit element complementing resistors, capacitors and inductors. Hallmarks of memristive behavior include pinched and frequency-dependent I-V hysteresis loops and most importantly a functional dependence of the magnetic flux passing through an ideal memristor on its electrical charge. Microtubules (MTs), cylindrical protein polymers composed of tubulin dimers are key components of the cytoskeleton. They have been shown to increase solution's ionic conductance and re-orient in the presence of electric fields. It has been hypothesized that MTs also possess intrinsic capacitive and inductive properties, leading to transistor-like behavior. Here, we show a theoretical basis and experimental support for the assertion that MTs under specific circumstances behave consistently with the definition of a memristor. Their biophysical properties lead to pinched hysteretic current-voltage dependence as well a classic dependence of magnetic flux on electric charge. Based on the information about the structure of MTs we provide an estimate of their memristance. We discuss its significance for biology, especially neuroscience, and potential for nanotechnology applications. MemristorsThe term memristor is the contraction of memory and resistor and it was first proposed in 1971 as the fourth element of the electric circuits 1 . A memristor is defined as a two-terminal passive circuit element that provides a functional relation between electric charge and magnetic flux 1,2 . The first physical realization of a memristor was achieved in 2008 2,3 and it has held a promise of nanoelectronics beyond Moore's law 4 , although this realization has been both difficult and controversial 5 . One of the possible breakthrough applications of memristors is neuromorphic computing 6 . Memristance refers to a property of the memristor that is analogous to resistance but it also depends on the history of applied voltage or injected current, unlike in other electrical circuit elements. When the electrical charge flows in one direction, the resistance of some memristors increases while it decreases when the charge flows in the opposite direction or vice versa. If the applied voltage is turned off, the memristor retains the last resistance value that it exhibited. This history dependence of memristance is expressed via a self-crossing or pinched I-V loop, which is frequency dependent 3,6 , and whose lobe area tends to zero as the frequency tends to infinity.A memristor is said to be charge-controlled if the relation between flux ϕ and charge q is: ϕ = ϕ (q). Conversely, it is said to be flux-controlled if q = q(ϕ). The voltage v of a charge-controlled memristor obeys a linear relationship with the current i(t) representing a charge-dependent Ohm's law such that: v(t) M(q) i(t)(1) = where memristance is defined as:

Tubulin, the primary subunit of microtubules, is remarkable for the variety of small molecules to which it binds. Many of these are very useful or promising agents in cancer chemotherapy. One of the most useful of these is paclitaxel. The tubulin molecule is itself an alpha/beta heterodimer, both alpha- and beta-tubulin monomers existing as multiple isotypes. Despite the success of paclitaxel as an anticancer drug, resistance often occurs in cancer cells and has been associated with variations in tubulin isotype expression, most notably with the increased expression of betaIII-tubulin. Paclitaxel is thought to reach its binding site on beta-tubulin by diffusion through nanopores in the microtubule wall. It has been suggested that a transitional step in this process may be the binding of paclitaxel to an intermediate site within a nanopore, from which it moves directly to its binding site in the microtubule interior facing the lumen. To test this hypothesis, we have computationally docked paclitaxel within a microtubule nanopore and simulated its passage to the intermediate binding site. Targeted molecular dynamics was then used to test the hypothesis that paclitaxel utilizes the H6/H7 loop as a hinge to move directly from this intermediate binding site to its final position in the luminal binding site. We observed that this motion appears to be stabilized by the formation of a hydrogen bond involving serine 275 in beta-tubulin isotypes I, IIa, IIb, IVa, IVb, V, VII, and VIII. Interestingly, this residue is replaced by alanine in the betaIII and VI isotypes. This observation raises the possibility that the observed isotype difference in paclitaxel binding may be a kinetic effect arising from the isotype difference at this residue. We are now able to suggest derivatives of paclitaxel that may reverse the isotype-specificity or lead to an alternate stabilizing hydrogen-bond interaction with tubulin, thus increasing the rate of passage to the luminal binding site and hopefully offering a therapeutic advantage in paclitaxel resistant cases.

Tubulin, an α/β heterodimer, has had most of its 3D structure analyzed; however, the carboxy (C)-termini remain elusive. Importantly, the C-termini play critical roles in regulating microtubule structure and function. They are sites of most of the post-translational modifications of tubulin and interaction sites with molecular motors and microtubule-associated proteins. Simulated annealing was used in our molecular dynamics modeling to predict the interactions of the C-terminal tails with the tubulin dimer. We examined differences in their flexibility, interactions with the body of tubulin, and the existence of structural motifs. We found that the α-tubulin tail interacts with the H11 helix of β-tubulin, and the β-tubulin tail interacts with the H11 helix of α-tubulin. Tail domains and H10/B9 loops interact with each other and compete for interactions with positively-charged residues of the H11 helix on the neighboring monomer. In a simulation in which α-tubulin's H10/B9 loop switches on sub-nanosecond intervals between interactions with the C-terminal tail of α-tubulin and the H11 helix of β-tubulin, the intermediate domain of α-tubulin showed more fluctuations compared to those in the other simulations, indicating that tail domains may cause shifts in the position of this domain. This suggests that C-termini may affect the conformation of the tubulin dimer which may explain their essential function in microtubule formation and effects on ligand binding to microtubules. Our modeling also provides evidence for a disordered-helical/helical double-state system of the T3/H3 region of the microtubule, which could be linked to depolymerization following GTP hydrolysis.