Delineating the interaction of combretastatin A-4 with αβ tubulin isotypes present in drug resistant human lung carcinoma using a molecular modeling approach

Abstract

Tubulin isotypes are known to regulate microtubule dynamic instability and contribute to the development of drug resistance in certain types of cancers. Combretastatin-A4 (CA-4) has a potent anti-mitotic, vascular disrupting and anti-angiogenic activity. It binds at the interface of αβ tubulin heterodimers and inhibits microtubules assembly. Interestingly, the CA-4 resistant human lung carcinoma shows alteration of βI and βIII isotype levels, a higher expression of βI tubulin isotype and a decreased expression of βIII tubulin isotypes has been reported in drug resistant cell lines. However, the origin of CA-4 resistance in lung carcinoma is not well understood. Here, we investigate the interaction and binding affinities of αβI, αβIIb, αβIII and αβIVa tubulin isotypes with CA-4, employing molecular modeling approaches. Sequence analysis shows that variations in residue composition at the CA-4 binding pocket of βI, βIII and βIVa tubulin isotypes when compared to template βIIb isotype. Molecular docking result shows that the CA-4 prefers ‘cis’ conformation in all αβ-tubulin isotypes. Molecular dynamics simulation reveal role of H7 helix, T7 loop and H8 helix of β-tubulin in lower binding affinity of αβI and αβIII isotypes for CA-4. The order of binding energy for CA-4 is αβIIb?>?αβIVa?>?αβI?>?αβIII. This suggest that drug resistance is induced in human lung carcinoma cells by altering the expression of β-tubulin isotypes namely βI and βIII which show lowest binding affinities. Our present study can help in designing potential CA-4 analogs against drug-resistant cancer cells showing altered expression of tubulin isotypes. Abbreviations: CA-4 combretastatin-A4

MD molecular dynamics

RMSD root mean square deviation

DSSP dictionary of secondary structure of proteins

VMD visual molecular dynamics

Communicated by Ramaswamy H. Sarma     

Increased levels of tyrosinated alpha-, beta(III)-, and beta(IV)-tubulin isotypes in paclitaxel-resistant MCF-7 breast cancer cells

Paclitaxel (PTX), the diterpene alkaloid, is a potent anti-cancer drug and is routinely used for the treatment of breast and ovarian cancers. The cellular targets of PTX are microtubules, which are composed of alpha- and beta-tubulin. Development of PTX resistance in patients has been a major problem associated with cancer chemotherapy. In an effort to get insight into this phenomenon of drug resistance, a PTX-resistant cell line from MCF-7 breast cancer cells has been generated. Western analysis of the cell extracts revealed that the resistant cells contain 2-fold higher amount of tyrosinated alpha-tubulin than those of the wild-type MCF-7 cells. Similar analyses of beta-tubulin with the isotype-specific monoclonal antibodies demonstrated that the PTX-resistant cells contain 2.5-fold higher amounts of beta(III) and 1.5-fold higher amount of beta(IV)-tubulin, while no difference was observed in the level of beta(I) isotype. These results demonstrate for the first time that PTX resistance is associated with an increase in the level of tyrosinated alpha-tubulin.     

Molecular dynamics simulation and density functional theory studies on the active pocket for the binding of paclitaxel to tubulin

Paclitaxel (PTX) is used to treat various cancers, but it also causes serious side effects and resistance. To better design similar compounds with less toxicity and more activity against drug-resistant tumors, it is important to clearly understand the PTX-binding pocket formed by the key residues of active sites on β-tubulin. Using a docking method, molecular dynamics (MD) simulation and density functional theory (DFT), we identified some residues (such as Arg278, Asp26, Asp226, Glu22, Glu27, His229, Arg369, Lys218, Ser277 and Thr276) on β-tubulin that are the active sites responsible for interaction with PTX. Another two residues, Leu371 and Gly279, also likely serve as active sites. Most of these sites contact with the “southern hemisphere” of PTX; only one key residue interacts with the “northern hemisphere” of PTX. These key residues can be divided into four groups, which serve as active compositions in the formation of an active pocket for PTX binding to β-tubulin. This active binding pocket enables a very strong interaction (the strength is predicted to be in the range of −327.8 to −365.7 kJ mol⁻¹) between β-tubulin and PTX, with various orientated conformations. This strong interaction means that PTX possesses a high level of activity against cancer cells, a result that is in good agreement with the clinical mechanism of PTX. The described PTX pocket and key active residues will be applied to probe the mechanism of tumor cells resistant to PTX, and to design novel analogs with superior properties.     

Secondary mutations correct fitness defects in Toxoplasma gondii with dinitroaniline resistance mutations

Dinitroanilines (oryzalin, trifluralin, ethafluralin) disrupt microtubules in protozoa but not in vertebrate cells, causing selective death of intracellular Toxoplasma gondii parasites without affecting host cells. Parasites containing α1-tubulin point mutations are dinitroaniline resistant but show increased rates of aberrant replication relative to wild-type parasites. T. gondii parasites bearing the F52Y mutation were previously demonstrated to spontaneously acquire two intragenic mutations that decrease both resistance levels and replication defects. Parasites bearing the G142S mutation are largely dependent on oryzalin for viable growth in culture. We isolated 46 T. gondii lines that have suppressed microtubule defects associated with the G142S or the F52Y mutations by acquiring secondary mutations. These compensatory mutations were α1-tubulin pseudorevertants or extragenic suppressors (the majority alter the β1-tubulin gene). Many secondary mutations were located in tubulin domains that suggest that they function by destabilizing microtubules. Most strikingly, we identified seven novel mutations that localize to an eight-amino-acid insert that stabilizes the α1-tubulin M loop, including one (P364R) that acts as a compensatory mutation in both F52Y and G142S lines. These lines have reduced dinitroaniline resistance but most perform better than parental lines in competition assays, indicating that there is a trade-off between resistance and replication fitness.     

Effects of amino-acid substitutions in β tubulin on benomyl sensitivity and microtubule functions inCoprinus cinereus

The sensitivity of the homobasidiomyceteCoprinus cinereus to the benzimidazole fungicide benomyl allowed us to isolate β-tubulin mutants as strains resistant to benomyl. To understand the molecular basis for the interaction between benomyl and β tubulin and for cellular defects in the β-tubulin mutants, we first analyzed the wild-type β1-tubulin gene (benA) ofC. cinereus, revealing thatbenA contains eight introns and encodes a 445 amino-acid protein. We then characterized 16 β1-tubulin mutants. The 16 mutations involved 11 different amino-acid substitutions at 10 different residues in β1 tubulin. The mutated residues were widely distributed along the primary sequence of β1 tubulin, from residue 3 in the N-terminal domain to residue 350 in the intermediate domain, but half of them appeared to be close to the αβ intradimer interface in an atomic model determined by electron crystallography. The benomyl resistant strain BEN 193, which exhibits clear heat sensitivity for hyphal growth and defects in various cellular processes, had a novel mutation, i.e., the Leu to Phe substitution at residue 350. Benomyl resistance and the heat sensitivity in BEN 193 were suppressed by additional amino-acid substitutions at various residues in β1 tubulin, suggesting that conformational changes of β1 tubulin are involved in the alterations. The DDBJ/GeneBank/EMBL accession number for the sequence reported in this paper is AB000116.     

Comparison of the EGFR resistance mutation profiles generated by EGFR-targeted tyrosine kinase inhibitors and the impact of drug combinations

Recent clinical data indicates that the emergence of mutant drug-resistant kinase alleles may be particularly relevant for targeted kinase inhibitors. In order to explore how different classes of targeted therapies impact upon resistance mutations, we performed EGFR (epidermal-growth-factor receptor) resistance mutation screens with erlotinib, lapatinib and CI-1033. Distinct mutation spectra were generated with each inhibitor and were reflective of their respective mechanisms of action. Lapatinib yielded the widest variety of mutations, whereas mutational variability was lower in the erlotinib and CI-1033 screens. Lapatinib was uniquely sensitive to mutations of residues located deep within the selectivity pocket, whereas mutation of either Gly(796) or Cys(797) resulted in a dramatic loss of CI-1033 potency. The clinically observed T790M mutation was common to all inhibitors, but occurred with varying frequencies. Importantly, the presence of C797S with T790M in the same EGFR allele conferred complete resistance to erlotinib, lapatinib and CI-1033. The combination of erlotinib and CI-1033 effectively reduced the number of drug-resistant clones, suggesting a possible clinical strategy to overcome drug resistance. Interestingly, our results also indicate that co-expression of ErbB2 (v-erb-b2 erythroblastic leukaemia viral oncogene homologue 2) has an impact upon the EGFR resistance mutations obtained, suggesting that ErbB2 may play an active role in the acquisition of drug-resistant mutations.     

Mutations of the β‐tubulin gene fragments from carbendazim‐resistant isolates of Pestalotiopsis sp. causing strawberry leaf blight in Chiang Mai,Thailand

Screening of field isolates of Pestalotiopsis sp. from strawberry leaf blight in Thailand identified 56 carbendazim‐resistant isolates. Of 56 isolates, 39 highly resistant (HR) phenotypes grew well on PDA amended with carbendazim even at ≥500 mg/L. Isolates with carbendazim‐resistant phenotype had a conspicuous mutation at particular sites in the beta‐tubulin (β‐tubulin) gene sequence. A β‐tubulin encoding gene from this pathogen was cloned and sequenced. Analysis of the β‐tubulin gene in highly resistant (HR) isolates showed a substitution at codon 79 and 102 from serine (AGT) to lysine (AAA) and valine (GTA) to alanine (GCA), respectively. The detection of such point mutations in the β‐tubulin gene allows the rapid screening to detect carbendazim‐resistant isolates in the field.     

Paclitaxel resistance in cells with reduced beta-tubulin

We previously described the isolation of colcemid resistant Chinese hamster ovary cell lines containing alpha- and beta-tubulin mutations that increase microtubule assembly and stability. By analyzing colcemid sensitive revertants from one of the beta-tubulin mutants, we now find that loss or inactivation of the mutant allele represents the most common mechanism of reversion. Consistent with this loss, the revertants have 35% less tubulin at steady state, no evidence for the presence of a mutant polypeptide, and a normal extent of tubulin polymerization. In addition to the loss of colcemid resistance, the revertant cells exhibit increased resistance to paclitaxel relative to wild-type cells. This paclitaxel resistance can be suppressed by transfecting the revertant cells with a cDNA for wild-type beta-tubulin, indicating that the reduction in tubulin in the revertant cells is responsible for the resistance phenotype. We propose that reducing tubulin levels may represent a novel mechanism of paclitaxel resistance.     

Induced transposition of Ds by a stable Ac in crosses of transgenic tobacco plants

Summary Rhizoxin and ansamitocin P-3 (a maytansinoid compound), potent inhibitors of mammalian brain tubulin assembly, inhibit growth of a variety of fungi including Aspergillus nidulans. Mutants of A. nidulans, benA10 which is a benomyl resistant -tubulin gene mutant and tubAl which is a benomyl supersensitive a-tubulin gene mutant, were both sensitive to rhizoxin and ansamitocin P-3 to the same extent as wild-type strains. We isolated 18 rhizoxin resistant mutants of A. nidulans. All of these mutants were cross-resistant to ansamitocin P-3, but not to benzimidazole antimitotic drugs. These mutants mapped to two loci, rhiA and rhiB, and all of those with high resistance mapped to rhiA. The fact that the protein extracts of rhiA mutants lost rhizoxin binding affinity and that rhiA was closely linked to benA, the major -tubulin gene in A. nidulans, indicated that rhiA must be a structural gene for -tubulin and that rhiA mutants are a new class of -tubulin gene mutants. All of this suggested that, in A. nidulans, these antimitotic drugs bind to -tubulin, and that rhizoxin and ansamitocin P-3 share the same binding site but the site does not overlap with the benzimidazole binding site. Protein extracts from a rhiB mutant retained rhizoxin binding affinity, therefore this rhizoxin resistance mechanism should not be a tubulin mediated process.     

Mutations in the β-tubulin binding site for peloruside A confer resistance by targeting a cleft significant in side chain binding

Peloruside A is a microtubule-stabilizing macrolide that binds to beta tubulin at a site distinct from the taxol site. The site was previously identified by H-D exchange mapping and molecular docking as a region close to the outer surface of the microtubule and confined in a cavity surrounded by a continuous loop of protein folded so as to center on Y340. We have isolated a series of peloruside A-resistant lines of the human ovarian carcinoma cell line A2780(1A9) to better characterize this binding site and the consequences of altering residues in it. Four resistant lines (Pel A-D) are described with single-base mutations in class I β-tubulin that result in the following substitutions: R306H, Y340S, N337D, and A296S in various combinations. The mutations are localized to peptides previously identified by Hydrogen-Deuterium exchange mapping, and center on a cleft in which the drug side chain appears to dock. The Pel lines are 10-15-fold resistant to peloruside A and show cross resistance to laulimalide but not to any other microtubule stabilizers. They show no cross-sensitivity to any microtubule destabilizers, nor to two drugs with targets unrelated to microtubules. Peloruside A induces G2/M arrest in the Pel cell lines at concentrations 10-15 times that required in the parental line. The cells show notable changes in morphology compared to the parental line.     

An <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis> tubulin mutant with conditional root-skewing phenotype

Low doses of microtubule-interacting drugs cause wild-type Arabidopsis thaliana seedling roots to twist in a left-handed helical direction. We here report molecular characterization of an A. thaliana tubulin mutant whose roots twist in a right-handed direction and have shallow left-handed cortical microtubule arrays when challenged with low doses of microtubule drugs. In the absence of the drug, growth and development of the mutant was apparently normal. In this conditional twisting mutant, Cys213 of α-tubulin6 was exchanged with Tyr. The mutant tubulin was incorporated into the microtubule polymer with wild-type tubulins, and thus acted as a dominant-negative mutation. These results suggest that compromised microtubules in wild-type and mutant roots are qualitatively distinct and affect skewing direction differently.     

Double mutation in Eleusine indica α-tubulin increases the resistance of transgenic maize calli to dinitroaniline and phosphorothioamidate herbicides

The repeated use of dinitroaniline herbicides on the cotton and soybean fields of the southern United States has resulted in the appearance of resistant biotypes of one of the world's worst weeds, Eleusine indica. Two biotypes have been characterized, a highly resistant (R) biotype and an intermediate resistant (I) biotype. In both cases the resistance has been attributed to a mutation in α-tubulin, a component of the α/β tubulin dimer that is the major constituent of microtubules. We show here that the I-biotype mutation, like the R-biotype mutation shown in earlier work, can confer dinitroaniline resistance on transgenic maize calli. The level of resistance obtained is the same as that for E. indica I- or R-biotype seedlings. The combined I- and R-biotype mutations increase the herbicide tolerance of transgenic maize calli by a value close to the summation of the maximum herbicide tolerances of calli harbouring the single mutations. These data, taken together with the position of the two different mutations within the atomic structure of the α/β tubulin dimer, imply that each mutation is likely to exert its effect by a different mechanism. These mechanisms may involve increasing the stability of microtubules against the depolymerizing effects of the herbicide or changing the conformation of the α/β dimer so that herbicide binding is less effective, or a combination of both possibilities.     

Involvement of a particular species of beta-tubulin (beta 3) in conidial development in Aspergillus nidulans

Strains of Aspergillus containing the benA22 mutation are resistant to benomyl for vegetative growth but do not produce conidia. To test whether conidiation involved an additional benomyl-sensitive tubulin (i.e., was mediated by a tubulin other than the tubulins coded for by the benA locus), a collection of mutants was produced that formed conidia in the presence of benomyl, i.e., were conidiation-resistant (CR-) mutants. We analyzed the tubulins of these CR- mutants using two-dimensional gel electrophoresis and found that the mutants lacked one species of beta-tubulin (designated beta 3). We have examined two of these mutants in detail. In crosses with strains containing wild-type tubulins, we found that the absence of the beta 3-tubulin co-segregated perfectly with the CR- phenotype. In diploids containing both the benA22 and CR- mutations, we found that the CR- phenotype was recessive and that beta 3-tubulin was present on two-dimensional gels of tubulins prepared from these diploids. In another set of crosses, these two CR- strains and seven others were first made auxotrophic for uridine and then crossed against strains that had homologously integrated a plasmid containing an incomplete internal fragment of the beta 3-tubulin gene and the pyr4 gene of Neurospora crassa (which confers uridine prototrophy on transformants). If the CR- phenotype were produced by a mutation in a gene distinct from the structural gene for beta 3-tubulin (designated the tubC gene), then crossing over should have produced some CR+ segregants among the uridine auxotrophic progeny of the second cross. All of the uridine auxotrophs from this type of cross, however, showed the CR- phenotype, suggesting that the mutation in these strains is at or closely linked to the tubC locus. The most obvious explanation of these results is that beta 3-tubulin is ordinarily used during conidiation and the presence of this species of beta-tubulin renders conidiation sensitive to benomyl. In the CR- mutants, beta 3-tubulin is absent, and in the presence of the benA22 mutation the benomyl-resistant beta 1-and/or beta 2-tubulin substitutes for beta 3 to make conidiation benomyl resistant. We discuss these results and give two models to explain the interactions between these beta-tubulin species.     

<Emphasis Type="Italic">Agrobacterium</Emphasis>-mediated transformation of flax with a mutant tubulin gene responsible for resistance to dinitroaniline herbicides

Agrobacterium tumefaciens was used to transform fiber flax with the pBITUBA8 plasmid carrying the mutant α-tubulin gene imparting resistance to dinitroaniline herbicides and the nptII selective marker gene imparting resistance to kanamycin. The transformants were selected in parallel on media containing kanamycin and trifluralin (a dinitroaniline herbicide). The transgenic nature of the resultant regenerants resistant to dinitroaniline herbicides was confirmed by means of Southern blotting and polymerase chain reaction (PCR) analysis using specific probes for the ntpII gene and the gene of α-tubulin.     

Alcohol Dehydrogenase Mutants of Chinese Hamster Somatic Cells Resistant to Allyl Alcohol

Alcohol dehydrogenase (alcohol: NAD oxidoreductase, E.C. 1.1.1.1.) mutants of Chinese hamster somatic cells were isolated as resistant to allyl alcohol (ALL^R). The ALL^R phenotypes of the mutant clones were reproducible with high fidelity and stable over long intervals of growth in the absence of the selecting drug. Several mutants, Adh-1, Adh-2, Adh-9 and Adh-13, resistant to allyl alcohol were characterized. They have between 15 and 40% of the alcohol dehydrogenase activity of the wild-type cell lines. Cell-cell hybridization experiments using Adh-1 and wild-type Chinese hamster cells indicate that resistance to allyl alcohol is recessive to the wild-type allele. This phenotype is therefore a useful marker to analyze gene segregation of somatic cell mutations and to study the expression of the genes involved in the metabolism of ethanol in mammalian cells.     

Molecular determinants for α-tubulin methylation by SETD2

Post-translational modifications to tubulin are important for many microtubule-based functions inside cells. It was recently shown that methylation of tubulin by the histone methyltransferase SETD2 occurs on mitotic spindle microtubules during cell division, with its absence resulting in mitotic defects. However, the catalytic mechanism of methyl addition to tubulin is unclear. We used a truncated version of human wild type SETD2 (tSETD2) containing the catalytic SET and C-terminal Set2–Rpb1–interacting (SRI) domains to investigate the biochemical mechanism of tubulin methylation. We found that recombinant tSETD2 had a higher activity toward tubulin dimers than polymerized microtubules. Using recombinant single-isotype tubulin, we demonstrated that methylation was restricted to lysine 40 of α-tubulin. We then introduced pathogenic mutations into tSETD2 to probe the recognition of histone and tubulin substrates. A mutation in the catalytic domain (R1625C) allowed tSETD2 to bind to tubulin but not methylate it, whereas a mutation in the SRI domain (R2510H) caused loss of both tubulin binding and methylation. Further investigation of the role of the SRI domain in substrate binding found that mutations within this region had differential effects on the ability of tSETD2 to bind to tubulin versus the binding partner RNA polymerase II for methylating histones in vivo, suggesting distinct mechanisms for tubulin and histone methylation by SETD2. Finally, we found that substrate recognition also requires the negatively charged C-terminal tail of α-tubulin. Together, this study provides a framework for understanding how SETD2 serves as a dual methyltransferase for both histone and tubulin methylation.