Modeling the Colchicum autumnale Tubulin and a Comparison of Its Interaction with Colchicine to Human Tubulin.

Tubulin is the target for many small-molecule natural compounds, which alter microtubules dynamics, and lead to cell cycle arrest and apoptosis. One of these compounds is colchicine, a plant alkaloid produced by Colchicum autumnale. While C. autumnale produces a potent cytotoxin, colchicine, and expresses its target protein, it is immune to colchicine's cytotoxic action and the mechanism of this resistance is hitherto unknown. In the present paper, the molecular mechanisms responsible for colchicine resistance in C. autumnale are investigated and compared to human tubulin. To this end, homology models for C. autumnale α-ß tubulin heterodimer are created and molecular dynamics (MD) simulations together with molecular mechanics Poisson-Boltzmann calculations (MM/PBSA) are performed to determine colchicine's binding affinity for tubulin. Using our molecular approach, it is shown that the colchicine-binding site in C. autumnale tubulin contains a small number of amino acid substitutions compared to human tubulin. However, these substitutions induce significant reduction in the binding affinity for tubulin, and subsequently fewer conformational changes in its structure result. It is suggested that such small conformational changes are insufficient to profoundly disrupt microtubule dynamics, which explains the high resistance to colchicine by C. autumnale.

Detailed Per-residue Energetic Analysis Explains the Driving Force for Microtubule Disassembly.

Microtubules are long filamentous hollow cylinders whose surfaces form lattice structures of αß-tubulin heterodimers. They perform multiple physiological roles in eukaryotic cells and are targets for therapeutic interventions. In our study, we carried out all-atom molecular dynamics simulations for arbitrarily long microtubules that have either GDP or GTP molecules in the E-site of ß-tubulin. A detailed energy balance of the MM/GBSA inter-dimer interaction energy per residue contributing to the overall lateral and longitudinal structural stability was performed. The obtained results identified the key residues and tubulin domains according to their energetic contributions. They also identified the molecular forces that drive microtubule disassembly. At the tip of the plus end of the microtubule, the uneven distribution of longitudinal interaction energies within a protofilament generates a torque that bends tubulin outwardly with respect to the cylinder's axis causing disassembly. In the presence of GTP, this torque is opposed by lateral interactions that prevent outward curling, thus stabilizing the whole microtubule. Once GTP hydrolysis reaches the tip of the microtubule (lateral cap), lateral interactions become much weaker, allowing tubulin dimers to bend outwards, causing disassembly. The role of magnesium in the process of outward curling has also been demonstrated. This study also showed that the microtubule seam is the most energetically labile inter-dimer interface and could serve as a trigger point for disassembly. Based on a detailed balance of the energetic contributions per amino acid residue in the microtubule, numerous other analyses could be performed to give additional insights into the properties of microtubule dynamic instability.

Chemical synthesis, pharmacological evaluation and in silico analysis of new 2,3,3a,4,5,6-hexahydrocyclopenta[c]pyrazole derivatives as potential anti-mitotic agents.

We have synthesized new, biologically active mono- and di-substituted 2,3,3a,4,5,6-hexahydrocyclopenta[c]pyrazole derivatives bearing electron withdrawing groups and electron donating groups. These derivative structures were characterized by their spectral and analytical data. The newly synthesized hexahydropyrazole analogues were evaluated for their in vitro anticancer activity against breast and lung cancer cell lines using a cytotoxicity bioassay. To understand their mechanism of action, tubulin binding assays were performed which pointed to their binding to microtubules in a mode similar to but not identical to colchicine, as evidenced by their KD value evaluation. Computational docking studies also suggested binding near the colchicine binding site on tubulin. These results were further confirmed by colchicine-binding assays on the most active compounds, which indicated that they bound to tubulin near but not at the colchicine site. The moderate cytotoxic effects of these compounds may be due to the presence of electron donating groups on the para-position of the phenyl ring, along with the hexahydropyrazole core nucleus. The observed anti-cancer activity based on inhibition of microtubule formation may be helpful in designing more potent compounds with a hexahydropyrazole moiety.

Analysis of the strength of interfacial hydrogen bonds between tubulin dimers using quantum theory of atoms in molecules.

Microtubules are key structural elements that, among numerous biological functions, maintain the cytoskeleton of the cell and have a major role in cell division, which makes them important cancer chemotherapy targets. Understanding the energy balance that brings tubulin dimers, the building blocks of microtubules, together to form a microtubule is especially important for revealing the mechanism of their dynamic instability. Several studies have been conducted to estimate various contributions to the free energy of microtubule formation. However, the hydrogen-bond contribution was not studied before as a separate component. In this work, we use concepts such as the quantum theory of atoms in molecules to estimate the per-residue strength of hydrogen bonds contributing to the overall stability that brings subunits together in pair of tubulin heterodimers, across both the longitudinal and lateral interfaces. Our study shows that hydrogen bonding plays a major role in the stability of tubulin systems. Several residues that are crucial to the binding of vinca alkaloids are shown to be strongly involved in longitudinal microtubule stabilization. This indicates a direct relation between the binding of these agents and the effect on the interfacial hydrogen-bonding network, and explains the mechanism of their action. Lateral contacts showed much higher stability than longitudinal ones (-462 ± 70 vs. -392 ± 59 kJ/mol), which suggests a dramatic lateral stabilization effect of the GTP cap in the ß-subunit. The role of the M-loop in lateral stability in absence of taxol was shown to be minor. The B-lattice lateral hydrogen bonds are shown to be comparable in strength to the A-lattice ones (-462 ± 70 vs. -472 ± 46 kJ/mol). These findings establish the importance of hydrogen bonds to the stability of tubulin systems.

Molecular dynamics simulation dataset of a microtubule ring in electric field.

We present molecular dynamics (MD) trajectories of a single ring of B-lattice microtubule ring consisting of 13 tubulin heterodimers. The data contain trajectories of this molecular system ran under various conditions (two temperature values, three ionic strength values, three values of electric field (including no field), and four electric field orientations). Our data enable us to analyze the effects of the electric field on microtubule under a variety of conditions. This data set was a basis of our in silico discovery, which demonstrates that the electric field can open microtubule lattice [1].

Tubulin response to intense nanosecond-scale electric field in molecular dynamics simulation.

Intense pulsed electric fields are known to act at the cell membrane level and are already being exploited in biomedical and biotechnological applications. However, it is not clear if electric pulses within biomedically-attainable parameters could directly influence intra-cellular components such as cytoskeletal proteins. If so, a molecular mechanism of action could be uncovered for therapeutic applications of such electric fields. To help clarify this question, we first identified that a tubulin heterodimer is a natural biological target for intense electric fields due to its exceptional electric properties and crucial roles played in cell division. Using molecular dynamics simulations, we then demonstrated that an intense - yet experimentally attainable - electric field of nanosecond duration can affect the bß-tubulin's C-terminus conformations and also influence local electrostatic properties at the GTPase as well as the binding sites of major tubulin drugs site. Our results suggest that intense nanosecond electric pulses could be used for physical modulation of microtubule dynamics. Since a nanosecond pulsed electric field can penetrate the tissues and cellular membranes due to its broadband spectrum, our results are also potentially significant for the development of new therapeutic protocols.

Anesthetics act in quantum channels in brain microtubules to prevent consciousness.

The mechanism by which anesthetic gases selectively prevent consciousness and memory (sparing non-conscious brain functions) remains unknown. At the turn of the 20(th) century Meyer and Overton showed that potency of structurally dissimilar anesthetic gas molecules correlated precisely over many orders of magnitude with one factor, solubility in a non-polar, 'hydrophobic' medium akin to olive oil. In the 1980s Franks and Lieb showed anesthetics acted in such a medium within proteins, suggesting post-synaptic membrane receptors. But anesthetic studies on such proteins yielded only confusing results. In recent years Eckenhoff and colleagues have found anesthetic action in microtubules, cytoskeletal polymers of the protein tubulin inside brain neurons. 'Quantum mobility' in microtubules has been proposed to mediate consciousness. Through molecular modeling we have previously shown: (1) olive oil-like non-polar, hydrophobic quantum mobility pathways ('quantum channels') of tryptophan rings in tubulin, (2) binding of anesthetic gas molecules in these channels, and (3) capabilities for π-electron resonant energy transfer, or exciton hopping, among tryptophan aromatic rings in quantum channels, similar to photosynthesis protein quantum coherence. Here, we show anesthetic molecules can impair π-resonance energy transfer and exciton hopping in tubulin quantum channels, and thus account for selective action of anesthetics on consciousness and memory.

Designing and Testing of Novel Taxanes to Probe the Highly Complex Mechanisms by Which Taxanes Bind to Microtubules and Cause Cytotoxicity to Cancer Cells.

Our previous work identified an intermediate binding site for taxanes in the microtubule nanopore. The goal of this study was to test derivatives of paclitaxel designed to bind to this intermediate site differentially depending on the isotype of ß-tubulin. Since ß-tubulin isotypes have tissue-dependent expression--specifically, the ßIII isotype is very abundant in aggressive tumors and much less common in normal tissues--this is expected to lead to tubulin targeted drugs that are more efficacious and have less side effects. Seven derivatives of paclitaxel were designed and four of these were amenable for synthesis in sufficient purity and yield for further testing in breast cancer model cell lines. None of the derivatives studied were superior to currently used taxanes, however computer simulations provided insights into the activity of the derivatives. Our results suggest that neither binding to the intermediate binding site nor the final binding site is sufficient to explain the activities of the derivative taxanes studied. These findings highlight the need to iteratively improve on the design of taxanes based on their activity in model systems. Knowledge gained on the ability of the engineered drugs to bind to targets and bring about activity in a predictable manner is a step towards personalizing therapies.

Similarity-based virtual screening for microtubule stabilizers reveals novel antimitotic scaffold.

Microtubules are among the most studied and best characterized cancer targets identified to date. Many microtubule stabilizers have been introduced so far that work by disrupting the dynamic instability of microtubules causing mitotic block and apoptosis. However, most of these molecules, especially taxol and epothilone, suffer absorption, toxicity and/or resistance problems. Here we employ a novel similarity-based virtual screening approach in the hope of finding other microtubule stabilizers that perform better and have lower toxicity and resistance. Epothilones, discodermolide, eleutherobin and sarcodictyin A have been found to compete with taxanes for the ß-tubulin binding site, which suggests common chemical features qualifying for that. Our approach was based on similarity screening against all these compounds and other microtubule stabilizers, followed by virtual screening against the taxol binding site. Some novel hits were found, together with a novel highly rigid molecular scaffold. After visual manipulations, redocking and rescoring of this novel scaffold, its affinity dramatically increased in a promising trend, which qualifies for biological testing.

The DNA region of phage BF23 encoding receptor binding protein and receptor blocking lipoprotein lacks homology to the corresponding region of closely related phage T5.

Analysis of the DNA region upstream of the BF23 hrs gene revealed a genetic organisation similar to that of closely related phage T5. A gene encoding a lipoprotein (llp(BF23)) is located directly upstream of the gene encoding the receptor binding protein (hrs) but is transcribed in opposite direction. The gene is followed by four open reading frames transcribed in the same direction. The llp (BF23) gene product does not show similarity to the corresponding T5 Llp(T5) protein, however, like Llp(T5) does for FhuA it blocks the BtuB receptor for BF23 infection. While no similarity between BF23 and T5 was observed for the DNA region encoding Llp and the receptor binding protein, the flanking regions were highly similar. Based on our results we conclude that a genetic module, the receptor-binding/receptor-blocking module, exists in phages BF23 and T5. Due to exclusion of homologous recombination within this module, it is hereditary only as an intact module: separation of the receptor-binding gene from the receptor-blocking gene, which apparently results in reduced fitness of the phage due to inactivation of progeny phage by active receptor proteins in the outer membranes of lysed cells, is thus effectively prevented.