Holly Freedman

A Unique Mode of Microtubule Stabilization Induced by Peloruside A

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.

Dynamics of RASSF1A/MOAP-1 Association with Death Receptors†

ABSTRACT 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.

Identification and Characterization of an Intermediate Taxol Binding Site Within Microtubule Nanopores and a Mechanism for Tubulin Isotype Binding Selectivity

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.

Computational Predictions of Volatile Anesthetic Interactions with the Microtubule Cytoskeleton: Implications for Side Effects of General Anesthesia

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.

MICROTUBULE IONIC CONDUCTION AND ITS IMPLICATIONS FOR HIGHER COGNITIVE FUNCTIONS

The neuronal cytoskeleton has been hypothesized to play a role in higher cognitive functions including learning, memory and consciousness. Experimental evidence suggests that both microtubules and actin filaments act as biological electrical wires that can transmit and amplify electric signals via the flow of condensed ion clouds. The potential transmission of electrical signals via the cytoskeleton is of extreme importance to the electrical activity of neurons in general. In this regard, the unique structure, geometry and electrostatics of microtubules are discussed with the expected impact on their specific functions within the neuron. Electric circuit models of ionic flow along microtubules are discussed in the context of experimental data, and the specific importance of both the tubulin C-terminal tail regions, and the nano-pore openings lining the microtubule wall is elucidated. Overall, these recent results suggest that ions, condensed around the surface of the major filaments of the cytoskeleton, flow along and through microtubules in the presence of potential differences, thus acting as transmission lines propagating intracellular signals in a given cell. The significance of this conductance to the functioning of the electrically active neuron, and to higher cognitive function is also discussed.

Molecular dynamics modeling of tubulin C‐terminal tail interactions with the microtubule surface

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 α‐tubulins 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. Proteins 2011;

Model of ionic currents through microtubule nanopores and the lumen.

It has been suggested that microtubules and other cytoskeletal filaments may act as electrical transmission lines. An electrical circuit model of the microtubule is constructed incorporating features of its cylindrical structure with nanopores in its walls. This model is used to study how ionic conductance along the lumen is affected by flux through the nanopores, both with and without an external potential applied across its two ends. Based on the results of Brownian dynamics simulations, the nanopores were found to have asymmetric inner and outer conductances, manifested as nonlinear IV curves. Our simulations indicate that a combination of this asymmetry and an internal voltage source arising from the motion of the C-terminal tails causes cations to be pumped across the microtubule wall and propagate in both directions down the microtubule through the lumen, returning to the bulk solution through its open ends. This effect is demonstrated to add directly to the longitudinal current through the lumen resulting from an external voltage source applied across the two ends of the microtubule. The predicted persistent currents directed through the microtubule wall and along the lumen could be significant in directing the dissipation of weak, endogenous potential gradients toward one end of the microtubule within the cellular environment.

Identification of tubulin drug binding sites and prediction of relative differences in binding affinities to tubulin isotypes using digital signal processing

Microtubules are involved in numerous cellular processes including chromosome segregation during mitosis and, as a result, their constituent protein, tubulin, has become a successful target of several chemotherapeutic drugs. In general, these drugs bind indiscriminately to tubulin within both cancerous and healthy cells, resulting in unwanted side effects. However, differences between beta-tubulin isotypes expressed in a wide range of cell types may aid in the development of anti-tubulin drugs having increased specificity for only certain types of cells. Here, we describe a digital signal processing (DSP) method that is capable of predicting hot spots for the tubulin family of proteins as well as determining relative differences in binding affinities to these hot spots based only on the primary sequence of 10 human tubulin isotypes. Due to the fact that several drug binding sites have already been characterized within beta-tubulin, we are able to correlate hot spots with the binding sites for known chemotherapy drugs. We have also verified the accuracy of this method using the correlation between the binding affinities of characterized drugs and the tubulin isotypes. Additionally, the DSP method enables the rapid estimation of relative differences in binding affinities within the binding sites of tubulin isotypes that are yet to be experimentally determined.

Explicitly Solvated Ligand Contribution to Continuum Solvation Models for Binding Free Energies: Selectivity of Theophylline Binding to an RNA Aptamer

To provide more accurate computational estimates of binding free energies in solution from molecular dynamics (MD) simulations, a separate solvation contribution for the binding ligand is determined from a linear response treatment. We use explicit water coordinates for this term and combine with MM-PBSA (molecular mechanics, Poisson-Boltzmann, and surface area contributions) in a new approach (MM-PB/LRA-SA). To assess this method, application to the binding between theophylline and its derivatives to an RNA aptamer was performed and compared with experimental binding affinities. Explicitly solvated MD trajectories were generated with the same parameter set used in the previous work by Gouda et al., who compared the relative binding of these molecules by both the MM-PBSA and thermodynamic integration methods. Substituting the linear response term for the ligand in the MM-PB/LRA-SA approach led to an improvement upon MM-PBSA when compared with experimental and thermodynamic integration results at approximately twice the computational cost. The balance between accuracy and computational expense achieved using this method suggests potential advantages in applying it in the virtual drug-screening process.

Computational Prediction of the Heterodimeric and Higher-Order Structure of gpE1/gpE2 Envelope Glycoproteins Encoded by Hepatitis C Virus

ABSTRACT Despite the recent success of newly developed direct-acting antivirals against hepatitis C, the disease continues to be a global health threat due to the lack of diagnosis of most carriers and the high cost of treatment. The heterodimer formed by glycoproteins E1 and E2 within the hepatitis C virus (HCV) lipid envelope is a potential vaccine candidate and antiviral target. While the structure of E1/E2 has not yet been resolved, partial crystal structures of the E1 and E2 ectodomains have been determined. The unresolved parts of the structure are within the realm of what can be modeled with current computational modeling tools. Furthermore, a variety of additional experimental data is available to support computational predictions of E1/E2 structure, such as data from antibody binding studies, cryo-electron microscopy (cryo-EM), mutational analyses, peptide binding analysis, linker-scanning mutagenesis, and nuclear magnetic resonance (NMR) studies. In accordance with these rich experimental data, we have built an in silico model of the full-length E1/E2 heterodimer. Our model supports that E1/E2 assembles into a trimer, which was previously suggested from a study by Falson and coworkers (P. Falson, B. Bartosch, K. Alsaleh, B. A. Tews, A. Loquet, Y. Ciczora, L. Riva, C. Montigny, C. Montpellier, G. Duverlie, E. I. Pecheur, M. le Maire, F. L. Cosset, J. Dubuisson, and F. Penin, J. Virol. 89:10333–10346, 2015, https://doi.org/10.1128/JVI.00991-15 ). Size exclusion chromatography and Western blotting data obtained by using purified recombinant E1/E2 support our hypothesis. Our model suggests that during virus assembly, the trimer of E1/E2 may be further assembled into a pentamer, with 12 pentamers comprising a single HCV virion. We anticipate that this new model will provide a useful framework for HCV envelope structure and the development of antiviral strategies. IMPORTANCE One hundred fifty million people have been estimated to be infected with hepatitis C virus, and many more are at risk for infection. A better understanding of the structure of the HCV envelope, which is responsible for attachment and fusion, could aid in the development of a vaccine and/or new treatments for this disease. We draw upon computational techniques to predict a full-length model of the E1/E2 heterodimer based on the partial crystal structures of the envelope glycoproteins E1 and E2. E1/E2 has been widely studied experimentally, and this provides valuable data, which has assisted us in our modeling. Our proposed structure is used to suggest the organization of the HCV envelope. We also present new experimental data from size exclusion chromatography that support our computational prediction of a trimeric oligomeric state of E1/E2.