1 School of Biotechnology and Bioinformatics, DY Patil University, CBD Belapur, Navi Mumbai 400614, India
2 Advanced Centre for Treatment, Research and Education in Cancer, Kharghar, Navi Mumbai 410210, India
Human β tubulin displays resistance to drugs like Taxol and Vinblastine due to amino acids substitutions within and outside the drug binding site.
This study focuses on the effect of amino acid substitutions outside the drug binding site on drug resistance. Amino acid substitution like R306C (mut2) is associated with Taxol resistance and D197N (mut1) and K350N (mut3) are associated with Vinblastine resistance. However, the mechanism of resistance has not been understood yet. This study has attempted to investigate the mechanism of resistance.
SWISSMODEL server was used to model the wild and the mutant β subunits which were later considered for protein-protein and protein-ligand docking using HADDOCK and AutoDock 1.5.6 software respectively. Dimer mutants were generated using Swisspdbviewer. POCASA 1.1 server was used to calculate the overall effect of substitution on pocket volume and the effect of substitution on domain mobility was explored using GROMACS software.
From sequence perspective, amino acid replacement in all three positions viz. D197N (mut1), R306C (mut2) and K350N (mut3) were found to have a deleterious effect on the stability of the protein. This study was further confirmed through structural analysis. Change in hydrogen bonding pattern was observed within the site of substitution in modeled mut1 and mut3 which is known to be specifically involved in Vinblastine interaction. In mut2 associated with Taxol binding, the hydrogen bonding pattern remained unaltered. All three mutants showed better protein-protein (β-β) interactions compared to the wild-type. Pocket size analysis in β subunit revealed that Taxol binding site increased in size after substitution in mut2 compared to the wild-type. However, the size of the Vinblastine binding site in the dimer interface remained the same before and after the substitution in wild and the mutants. Wild-type (β monomer and αβ dimer) associated with Taxol and Vinblastine, respectively showed better drug interaction compared to their mutants.
This study throws light on the mechanism of drug resistance due to amino acid substitutions outside the drug binding site. It was found that amino acid substitution outside the drug site enhanced protein-protein interaction between the β-β subunits.
Keywords: R306C, D197N, Drug resistance, K350N, Taxol, Tubulin mutations, Vinblastine.
open-access license: This is an open access article distributed under the terms of the Creative Commons Attribution 4.0 International Public License (CC-BY 4.0), a copy of which is available at: (https://creativecommons.org/licenses/by/4.0/legalcode). This license permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
* Address correspondence to this author at the School of Biotechnology and Bioinformatics, DY Patil University, Selvaa Kumar C, CBD Belapur, Navi Mumbai 400614, India; Tel: 912227563600; Email: email@example.com
Microtubules play a key role in cell division, mobility, architecture, intracellular transportation and chromosomal segregation [1Chang P, Stearns T. δ-tubulin and ϵ-tubulin: two new human centrosomal tubulins reveal new aspects of centrosome structure and function. Nat Cell Biol 2000; 2(1): 30-5.[http://dx.doi.org/10.1038/71350] [PMID: 10620804] ]. The fundamental unit of a microtubule is a tubulin dimer made of α and β subunits [2Williams RC Jr, Shah C, Sackett D. Separation of tubulin isoforms by isoelectric focusing in immobilized pH gradient gels. Anal Biochem 1999; 275(2): 265-7.[http://dx.doi.org/10.1006/abio.1999.4326] [PMID: 10552916] ]. Longitudinal interactions occur between β-α and α-β monomer of tubulin classified as intra- and inter-dimer interactions, respectively. As per literature review, intra-dimer interface (non-exchangeable) is rigid compared to the inter-dimer interface with exchangeable GDP [3Caplow M, Ruhlen RL, Shanks J. The free energy for hydrolysis of a microtubule-bound nucleotide triphosphate is near zero: all of the free energy for hydrolysis is stored in the microtubule lattice. J Cell Biol 1994; 127(3): 779-88.[http://dx.doi.org/10.1083/jcb.127.3.779] [PMID: 7962059] , 4Mickey B, Howard J. Rigidity of microtubules is increased by stabilizing agents. J Cell Biol 1995; 130(4): 909-17.[http://dx.doi.org/10.1083/jcb.130.4.909] [PMID: 7642706] ]. There are specific loop regions namely H1-S2 and M loops which exhibit lateral electrostatic interactions between H1-S2 loop of protofilament1 with M loop of adjacent protofilament2. Their overall binding affinity stabilizes the microtubule [4Mickey B, Howard J. Rigidity of microtubules is increased by stabilizing agents. J Cell Biol 1995; 130(4): 909-17.[http://dx.doi.org/10.1083/jcb.130.4.909] [PMID: 7642706] , 5Nogales E, Whittaker M, Milligan RA, Downing KH. High-resolution model of the microtubule. Cell 1999; 96(1): 79-88.[http://dx.doi.org/10.1016/S0092-8674(00)80961-7] [PMID: 9989499] ].
Microtubules are also considered as potential targets for anticancer drugs as their disruption results in cell cycle arrest associated with apoptosis [6Shing JC, Choi JW, Chapman R, et al. A novel synthetic 1,3-phenyl bis-thiourea compound targets microtubule polymerization to cause cancer cell death. Cancer Biol Ther 2014; 15(7): 895-905.[http://dx.doi.org/10.4161/cbt.28881] [PMID: 24755487] ]. Many anticancer drugs like Taxol, Vinblastine, Vincristine and Colchicine cause cell cycle arrest [7Thorpe MF, Hespenheide BM, Yang Y, Kuhn LA. Flexibility and critical hydrogen bonds in cytochrome c. Pac Symp Biocomput 2000; 191-202.[PMID: 10902168] , 8Kuppens IE. Current state of the art of new tubulin inhibitors in the clinic. Curr Clin Pharmacol 2006; 1(1): 57-70.[http://dx.doi.org/10.2174/157488406775268200] [PMID: 18666378] ]. It has been reported that amino acid substitutions have led to drug resistance in β tubulin. These substitutions could be within Taxol binding site (A231T, F270V, T274I, R282N and A364T) or outside the drug binding pockets (R306) [9Berrieman HK, Lind MJ, Cawkwell L. Do β-tubulin mutations have a role in resistance to chemotherapy? Lancet Oncol 2004; 5(3): 158-64.[http://dx.doi.org/10.1016/S1470-2045(04)01411-1] [PMID: 15003198] , 10Verma K, Ramanathan K. Investigation of paclitaxel resistant R306C mutation in β-tubulin—A computational approach. J Cell Biochem 2015; 116(7): 1318-24.[http://dx.doi.org/10.1002/jcb.25087] [PMID: 25735511] ]. For Vinca alkaloid binding site, substitution within the binding pocket is S172A and substitution observed outside are D197N, K350N which causes altered drug binding affinity (Fig. 1). Given the high incidence of tubulin mutation as a cause of drug resistance in these studies, it is imperative to investigate the role of different amino acid substitutions in drug resistance [6Shing JC, Choi JW, Chapman R, et al. A novel synthetic 1,3-phenyl bis-thiourea compound targets microtubule polymerization to cause cancer cell death. Cancer Biol Ther 2014; 15(7): 895-905.[http://dx.doi.org/10.4161/cbt.28881] [PMID: 24755487] , 9Berrieman HK, Lind MJ, Cawkwell L. Do β-tubulin mutations have a role in resistance to chemotherapy? Lancet Oncol 2004; 5(3): 158-64.[http://dx.doi.org/10.1016/S1470-2045(04)01411-1] [PMID: 15003198] , 11Poruchynsky MS, Kim J-H, Nogales E, et al. Tumor cells resistant to a microtubule-depolymerizing hemiasterlin analogue, HTI-286, have mutations in α- or β-tubulin and increased microtubule stability. Biochemistry 2004; 43(44): 13944-54.[http://dx.doi.org/10.1021/bi049300+] [PMID: 15518543] -13Giannakakou P, Gussio R, Nogales E, et al. A common pharmacophore for epothilone and taxanes: molecular basis for drug resistance conferred by tubulin mutations in human cancer cells. Proc Natl Acad Sci USA 2000; 97(6): 2904-9.[http://dx.doi.org/10.1073/pnas.040546297] [PMID: 10688884] ]. Existing literature review investigated residues within the drug binding site [14Gökmen-Polar Y, Escuin D, Walls CD, et al. β-Tubulin mutations are associated with resistance to 2-methoxyestradiol in MDA-MB-435 cancer cells. Cancer Res 2005; 65(20): 9406-14.[http://dx.doi.org/10.1158/0008-5472.CAN-05-0088] [PMID: 16230404] -17Verma K, Ramanathan K. Exploring the impact of F270V mutation in the β-tubulin (Bos Taurus) structure and its function: a computational perspective. Biotechnol Lett 2015; 37(5): 1003-11.[http://dx.doi.org/10.1007/s10529-015-1765-9] [PMID: 25631779] ]. However, not much has been reported for the substitutions outside the drug binding site. In the present study, the effect of these three mutations D197N, R306C and K350N outside the drug binding pocket and their role in drug resistance will be understood from docking and dynamics perspective.
Human beta tubulin with mutations outside the drug binding sites (R306C, D197N, K350N are highlighted in wireframe format.
2. MATERIALS AND METHODS
2.1. Amino Acid Composition Analysis
β tubulin of Homo sapiens was downloaded from Uniprot database (Accession number: Q9BVA1) [18Bairoch A, Boeckmann B, Ferro S, Gasteiger E. Swiss-Prot: juggling between evolution and stability. Brief Bioinform 2004; 5(1): 39-55.[http://dx.doi.org/10.1093/bib/5.1.39] [PMID: 15153305] ]. To understand the effect of substitutions at positions 197, 306 and 350 we substituted each of them individually to generate D197N (mut1), R306C (mut2) and K350N (mut3). These mutants protein sequences were submitted to Amino Acid substitution (AAS) tools like PANTHER (Protein Analysis Through Evolutionary Relationships) [19Thomas PD, Campbell MJ, Kejariwal A, et al. PANTHER: a library of protein families and subfamilies indexed by function. Genome Res 2003; 13(9): 2129-41.[http://dx.doi.org/10.1101/gr.772403] [PMID: 12952881] ] and Polyphen 2 (Polymorphism Phenotyping v2) [20Adzhubei IA, Schmidt S, Peshkin L, et al. A method and server for predicting damaging missense mutations. Nat Methods 2010; 7(4): 248-9.[http://dx.doi.org/10.1038/nmeth0410-248] [PMID: 20354512] ] to study the effect of amino acid substitutions on structural stability. Consensus reports of these two softwares were considered for further analysis.
2.2. Homology Modeling and Mutant Generation
The protein sequence of human β tubulin was considered for homology modeling using SWISSMODEL server [21Guex N, Peitsch MC. SWISS-MODEL and the Swiss-PdbViewer: an environment for comparative protein modeling. Electrophoresis 1997; 18(15): 2714-23.[http://dx.doi.org/10.1002/elps.1150181505] [PMID: 9504803] ] to model the wild type. Of the listed templates, the 3D structure with better query coverage and amino acid identity was selected. Similar step was followed for modeling mut2 wherein Arg at 306th position was substituted by Cys. Regarding Vinblastine docking, instead of modeling and generating dimer, here we downloaded the 3D structure of Bos taurus (PDB id: 4I50) [22Prota AE, Bargsten K, Zurwerra D, et al. Molecular mechanism of action of microtubule-stabilizing anticancer agents. Science 2013; 339(6119): 587-90.[http://dx.doi.org/10.1126/science.1230582] [PMID: 23287720] ] from Protein Data Bank [23Berman HM, Westbrook J, Feng Z, et al. The protein data bank. Nucleic Acids Res 2000; 28(1): 235-42.[http://dx.doi.org/10.1093/nar/28.1.235] [PMID: 10592235] ] which showed 99% amino acid identity with Homo sapiens. This has four subunits of which αβ dimer was extracted using Swisspdbviewer [21Guex N, Peitsch MC. SWISS-MODEL and the Swiss-PdbViewer: an environment for comparative protein modeling. Electrophoresis 1997; 18(15): 2714-23.[http://dx.doi.org/10.1002/elps.1150181505] [PMID: 9504803] ]. Furthermore, mut1 and mut3 were generated through the substitution of Asp and Lys with Asn at position 197 and 350, respectively. All the modelled structures were energy minimized using Swisspdbviewer software and were further validated using an PROSA (Protein Structure Analysis Software) [24Sippl MJ. Recognition of errors in three-dimensional structures of proteins. Proteins 1993; 17(4): 355-62.[http://dx.doi.org/10.1002/prot.340170404] [PMID: 8108378] ].
2.3. Hydrogen Bonding Pattern
Hydrogen bonding patterns were separately investigated for the wild and the mutants using CHIMERA software [25Pettersen EF, Goddard TD, Huang CC, et al. UCSF Chimera--a visualization system for exploratory research and analysis. J Comput Chem 2004; 25(13): 1605-12.[http://dx.doi.org/10.1002/jcc.20084] [PMID: 15264254] ].
2.4. Protein-Protein Docking
β subunit of the modeled wild and the mutants (mut1, mut2 and mut3) were considered for protein-protein docking (β-β) using HADDOCK server [26de Vries SJ, van Dijk M, Bonvin AM. The HADDOCK web server for data-driven biomolecular docking. Nat Protoc 2010; 5(5): 883-97.[http://dx.doi.org/10.1038/nprot.2010.32] [PMID: 20431534] ] to investigate their lateral binding affinity. The entire docking process consisted of three stages. They were (a) rigid-body energy minimization (b) a semi-flexible refinement torsion angle space and (c) a final refinement in explicit solvent. After each stage, structures were scored and ranked accordingly, and further selected for the next step. Moreover, HADDOCK score is the summation of van der Waals, electrostatic, desolvation and restraint violation energies with a buried surface area. From the webpage of HADDOCK, Easy interface module was selected to upload wild tubulins as molecule1 and molecule2. Regarding the active residues, sub-segments critically involved in lateral interactions like H1-S2 loop (residues 26-63) and M loop (residues 272-288) were considered for molecule1 and molecule2 respectively. Similar steps were followed for the mutant structures.
2.5. Pocket Identification
POCASA 1.1 software [27Yu J, Zhou Y, Tanaka I, Yao M. Roll: a new algorithm for the detection of protein pockets and cavities with a rolling probe sphere. Bioinformatics 2010; 26(1): 46-52.[http://dx.doi.org/10.1093/bioinformatics/btp599] [PMID: 19846440] ] was mainly used for measuring the ligand binding pocket size in monomer (for Taxol) and αβ dimer interface (for Vinblastine) of wild and mutants. Here, the probe radius and the grid size was 2 Å and 1.0Å respectively.
2.6. Protein-Ligand Docking
Modelled wild and the mutants were considered for docking using AUTODOCK tool 1.5.6 [28Morris GM, Huey R, Lindstrom W, et al. AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility. J Comput Chem 2009; 30(16): 2785-91.[http://dx.doi.org/10.1002/jcc.21256] [PMID: 19399780] ]. Drugs like Taxol and Vinblastine were obtained from 1JFF [29Löwe J, Li H, Downing KH, Nogales E. Refined structure of α β-tubulin at 3.5 A resolution. J Mol Biol 2001; 313(5): 1045-57.[http://dx.doi.org/10.1006/jmbi.2001.5077] [PMID: 11700061] ] and 1Z2B [30Gigant B, Wang C, Ravelli RB, et al. Structural basis for the regulation of tubulin by vinblastine. Nature 2005; 435(7041): 519-22.[http://dx.doi.org/10.1038/nature03566] [PMID: 15917812] ] which were separated from the receptor using Swisspdbviewer. Kollman charges were added to the modelled wild and the mutant receptors which were saved as .pdbqt files. Gasteiger charges were added to the ligands and saved as .pdbqt files. Firstly, wild type β subunit and Taxol were considered for docking. For the same, Gridbox was constructed near Taxol binding site (active site as reported in the crystal structure of 1JFF) wherein the number of grid points for x,y and z direction are 40,28 and 48, respectively. The grid centre for x,y and z coordinates are 6.855,-0.611 and -26.64, respectively. Same Gridbox value was considered for mut2 protein. Secondly, wild αβ dimer and Vinblastine were considered for docking. Basically, Vinca alkaloid prefers to bind to the dimer interface of the α-β subunits. The Gridbox is constructed near the inter-dimer interface (active sites as reported in the crystal structure of 1Z2B) wherein the number of grid point for x,y and z direction are 32,52 and 36, respectively. The grid centre for x, y and z coordinates are -0.413, 41.812 and 18.333, respectively. Same Gridbox values were considered for mut1 and mut3 proteins associated with Vinblastine docking. All docked poses were visualized using Discover Studio Visualizer software [31Discovery Studio Modeling Environment, Release xx 2007.].
2.7. Protein Simulation
The coordinates of wild and mutants (monomer and dimer) were considered for molecular simulation using GROMACS 4.5.5 [32Hess B, Kutzner C, van der Spoel D, Lindahl E. GROMACS 4: algorithms for highly efficient, load-balanced, and scalable molecular simulation. J Chem Theory Comput 2008; 4(3): 435-47.[http://dx.doi.org/10.1021/ct700301q] [PMID: 26620784] ] with implementation of OPLS-AA/L force field [33Jorgensen WL, Tirado-Rives J. The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. J Am Chem Soc 1988; 110(6): 1657-66.[http://dx.doi.org/10.1021/ja00214a001] [PMID: 27557051] ]. These systems were solvated using TIP3P water modeling a cubic box with periodic boundary conditions [34Jorgensen WL, Chandrasekhar J, Madura JD, Impey RW, Klein ML. Comparison of simple potential functions for simulating liquid water. J Chem Phys 1983; 79: 926-35.[http://dx.doi.org/10.1063/1.445869] ]. Additionally, counter ions were added to neutralize the system. These systems were first energy minimized using steepest descent algorithm with a tolerance of 1000 kJ/mol/nm. Long-range electrostatic interactions were calculated using Particle Mesh Ewald (PME) summation with 1 nm cutoffs for Columbic interactions, and Van der Waal interactions were calculated with a distance cut-off of 1.4 nm. Later on, systems were equilibrated by applying positional restraints on the structure using NVT followed by NPT ensemble for 100 ps each. Temperature of 300 K was coupled using Berendsen thermostat [35Berendsen HJ, Postma JV, van Gunsteren WF, DiNola A, Haak J. Molecular dynamics with coupling to an external bath. J Chem Phys 1984; 81: 3684-90.[http://dx.doi.org/10.1063/1.448118] ] with pressure at 1 bar, coupled by Parrinello-Rahman algorithm [36Car R, Parrinello M. Unified approach for molecular dynamics and density-functional theory. Phys Rev Lett 1985; 55(22): 2471-4.[http://dx.doi.org/10.1103/PhysRevLett.55.2471] [PMID: 10032153] ]. The equilibrated systems were then subjected to 15 ns of production run with time-step integration of 2 fs. The trajectories were saved at every 2 ps and analyzed using analysis tools from GROMACS 4.5.5. Basic analyses like Root Mean Square Deviation (RMSD) and Root Mean Square Fluctuations (RMSF) were calculated individually for wild and the mutants. RMSD measure the average distance between the atoms (usually the backbone atoms) of superimposed proteins [37Coutsias EA, Seok C, Dill KA. Using quaternions to calculate RMSD. J Comput Chem 2004; 25(15): 1849-57.[http://dx.doi.org/10.1002/jcc.20110] [PMID: 15376254] ]. The RMSF is one of the best methods of comparing dynamics. This calculation finds the degree of movement of each Cα around its average position, i.e. parts of the protein that are highly flexible will have a large RMSF value, while portions that are contained will result in a low RMSF [38Mitra A, Sept D. Taxol allosterically alters the dynamics of the tubulin dimer and increases the flexibility of microtubules. Biophys J 2008; 95(7): 3252-8.[http://dx.doi.org/10.1529/biophysj.108.133884] [PMID: 18621813] ].
3.1. Amino Acid Substitution Analysis
Three β tubulin mutants mut1(D197N), mut2(R306C) and mut3(K350N) were generated through amino acid substitutions tools and as per their consensus report all three substitutions could probably damage the protein stability Table (1). However, this needs further confirmation from a structural perspective. For the same, the structures were modelled in our next step.
Table 1 Summary of amino acid substitution tools based analysis of beta tubulin using PANTHER and Polyphen2.
3.2. Homology Modeling and Mutant Generation
SWISSMODEL server was considered for the modeling of human wild β subunit and mut2 which was later considered for Taxol binding. This model was based on the template of 4I50 of Bos Taurus with query coverage of 99%. For Vinblastine docking, αβ dimer itself was extracted from PDB id 4I50. Both wild and the mutant models showcased good local model quality below 0 for a window size of 40 residues (Supplementary material-S1 (a-d)).
3.3. Hydrogen Bonding Pattern
Hydrogen bonding patterns revealed that in wild β subunit, salt bridges were observed between Asp197 and Arg158 of helix4 wherein OD1 and OD2 of Asp197 interacted with the hydrogen atom of Arg158. This interaction is basically associated with the anchoring of COOH-terminal end of helix 4 during lateral protofilament interaction as per literature review [39Gökmen-Polar Y, Escuin D, Walls CD, et al. β-Tubulin mutations are associated with resistance to 2-methoxyestradiol in MDA-MB-435 cancer cells. Cancer Res 2005; 65(20): 9406-14.[http://dx.doi.org/10.1158/0008-5472.CAN-05-0088] [PMID: 16230404] ]. Furthermore, hydrogen atom of Asp197 (renumbered as 199 in wild-type after modeling) interacts with the oxygen and hydrogen atom of Met166 and Arg158. Thus, in total, there are five interactions observed between Asp199 and the other two amino acids (Met166, Arg158). After the substitution, mut1 displays only three hydrogen bonds between Met166 and Arg158. In Arg306 (wild) and mut2 the hydrogen bonding pattern remains the same. With respect to Lys350 (wild), the hydrogen bonds were observed between Thr353 and Ala317. However, in mut3 Asn350 interacts with Ala315 and Val313. With substitution, the contact with Thr353 ceases. Overall, mut1 and mut3 showcases fewer hydrogen bonds which are otherwise unscathed in mut2 (Fig. 2a-f).
Hydrogen bonding pattern analysis before and after the substitution of amino acids in human β tubulin (a) wild type197, (b) wildtype 306 (c) wildtype 350 (d)mut1(e)mut2(f)mut3.
3.4. Protein-Protein Docking
Wild and the mutant models were considered for protein-protein docking to calculate their overall binding affinity and the protofilament interactions before and after mutations. In wild beta subunits, the binding affinity score is -76.9 +/-8.4 with a buried surface of 1643.5 +/-64.4. Conversely, the binding affinity of mut1, mut2 and mut3 were -112.6 +/-4.0,-106.0 +/-19.1,-92 +/-4.2, respectively. Moreover, the buried surface of mut1, mut2 and mut3 were 1685.4 +/-96.4, 2108.3 +/-155.5 and 1607.0 +/-57.3 respectively. Thus, here a considerable increase in the binding affinity between mut1, mut2 and mut3 was observed compared to their wild counterpart. Similarly, a steady rise in the buried surface area was observed in mut1 and mut2 with an exception in mut3 (Fig. 3a-d). Thus, mutants showed better lateral interactions which could assist in better protofilament interactions during microtubule formation.
Protein-protein docking analysis of the wild and mutants (a) wild type (b)mut1(c)mut2 and (d)mut3.
3.5. Pocket Identification
As per the POCASA report, in wildtype β subunit the pocket size associated with Taxol interaction is 242Å3 which increases to 251Å3 in mut2. Thus, an increase in pocket size was observed after the substitution of arginine with cysteine. In Vinblastine binding dimer interface, the wild and the mutant dimmers displayed a pocket size of 268Å3 which remained unaltered after the substitution in D197N and K350N.
3.6. Protein-Ligand Docking
In protein-ligand docking, drug Taxol was docked with modeled wild and mut2. In wild type, the binding affinity is -8.1 kcal/mol wherein Taxol interacts with His229, Thr276 and Gln281. In mut2, the binding affinity is –6.4 kcal/mol and a single interaction with Thr276 was observed after the substitution of Arg with Cys. Here, we assume that an increase in pocket size volume after substitution could be the reason behind the absence two interactions (Figs. 4a, b). Furthermore, Vinblastine was docked with wild and the generated mutant dimers. In wild-type, the binding affinity is -7.9 kcal/mol and the drug interacts with Asn 329. In mut1 and mut3 the binding affinity is -7.4 kcal/mol and -6.8 kcal/mol respectively with no interaction between the receptor and the drug. The overall change in hydrogen bonding pattern within the site of substitution could be the reason behind the reduced binding affinity (Fig. 4c).
Protein-ligand docking analysis of the wild and the mutants (a) wild type with Taxol (b)mut2 with Taxol (c) wild type with Vinblastine.
3.7. Molecular Simulation
To estimate the effect of mutation on structural stability/flexibility of the tubulins with respect to the wild type, the RMSD and RMSF values were calculated individually. Here, we monitored the structural stability of the wild and the mutants throughout the simulation (Fig. 5a-d). In mut1, mut2 and mut3 there is a dramatic increase in the RMSD compared to the wild type. As per the RMSF report, in mut1 and mut2 the subsegment of D197N and R306C are flexible compared to the wild type. However, in mut3 the subsegment K350N shows an increase in rigidness compared to the wild type after substitution with asparagine (Fig. 6a-d).
Root Mean Square fluctuation analysis in wild and mutants (a) wild type (b)mut1(c)mut2 (d)mut3.
It is generally observed that substitutions within the drug binding site of β tubulin lead to resistance. Since there is a paucity of information about substitution outside the drug binding site, in the present study the role of amino acid substitutions outside the drug binding site in β tubulin from Homo sapiens has been investigated. Firstly, from sequence perspective, all three substitutions (D197N, R306C and K350N) were found to be probably damaging the protein stability. Secondly, hydrogen bonding network analysis of the wild and the mutant structures reveals a decrease in bonds in mut1 and mut3 compared to the wild-type. In particular, salt bridge formation between Asp197 and Arg158 is lost in mut1 after substitution which increases the flexibility of this region as observed during simulation. Mut2 also shows higher flexibility irrespective of no change in their bonding pattern. In mut3, substitution reduces the flexibility of this subsegment. Thirdly, protein-protein interaction analysis through docking confirms better interactions among mut1-mut1, mut2-mut2 and mut3-mut3 compared to the wildtype. Thus, a better protofilament interaction is inevitable after amino acid substitutions which could further enhance their domain stability. Fourthly, protein-ligand docking analysis reveals lower Vinblastine binding affinity in mut1 and mut3 compared to the wild-type due to the loss of salt bridge formation and higher domain flexibility which could ultimately reduce the protein-ligand interaction. Even in mut2, higher flexibility associated with larger pocket size ultimately hampered with the Taxol binding. This was evident from a single residue binding instead of three. In mut3, the domain rigidness plays a crucial role in reduced Vinblastine binding within the pocket. Thus, all three mutants have reduced drug binding and improved protein-protein interactions after substitutions which could certainly assist in better protofilament interactions much needed for their stability. Thus to conclude, even mutations outside the drug binding site have critical role in protein-protein and protein–ligand interaction which may play pivotal role during drug resistance.
In this study, we investigated the effect of amino acid substitutions outside the drug binding site on drug resistance. And it was observed that amino acid substitution outside the drug site enhanced protein-protein interaction between the β-β subunits.
Supplementary material is available on the publishers Website along with the published article.
No Animals/Humans were used for studies that are base of this research.
CONSENT FOR PUBLICATION
CONFLICT OF INTEREST
The authors declare no conflict of interest, financial or otherwise.
Authors would like to acknowledge the Department of Biotechnology (DBT), Government of India, sponsored Distributed Information Sub Centre (SubDIC) of Biotechnology Information System (BTIS) Network at ACTREC where analysis of MD trajectories and docking studies were carried out. Authors also thank Super-computing facility (Utkarsh) of Computer division, BARC, for providing access to molecular dynamics simulation of structures. This Project was not funded by any external funding agencies.
Chang P, Stearns T. δ-tubulin and ϵ-tubulin: two new human centrosomal tubulins reveal new aspects of centrosome structure and function. Nat Cell Biol 2000; 2(1): 30-5.[http://dx.doi.org/10.1038/71350] [PMID: 10620804]
Caplow M, Ruhlen RL, Shanks J. The free energy for hydrolysis of a microtubule-bound nucleotide triphosphate is near zero: all of the free energy for hydrolysis is stored in the microtubule lattice. J Cell Biol 1994; 127(3): 779-88.[http://dx.doi.org/10.1083/jcb.127.3.779] [PMID: 7962059]
Shing JC, Choi JW, Chapman R, et al. A novel synthetic 1,3-phenyl bis-thiourea compound targets microtubule polymerization to cause cancer cell death. Cancer Biol Ther 2014; 15(7): 895-905.[http://dx.doi.org/10.4161/cbt.28881] [PMID: 24755487]
Thorpe MF, Hespenheide BM, Yang Y, Kuhn LA. Flexibility and critical hydrogen bonds in cytochrome c. Pac Symp Biocomput 2000; 191-202.[PMID: 10902168]
Poruchynsky MS, Kim J-H, Nogales E, et al. Tumor cells resistant to a microtubule-depolymerizing hemiasterlin analogue, HTI-286, have mutations in α- or β-tubulin and increased microtubule stability. Biochemistry 2004; 43(44): 13944-54.[http://dx.doi.org/10.1021/bi049300+] [PMID: 15518543]
Giannakakou P, Gussio R, Nogales E, et al. A common pharmacophore for epothilone and taxanes: molecular basis for drug resistance conferred by tubulin mutations in human cancer cells. Proc Natl Acad Sci USA 2000; 97(6): 2904-9.[http://dx.doi.org/10.1073/pnas.040546297] [PMID: 10688884]
Navarrete KR, Alderete JB, Jiménez VA. Structural basis for drug resistance conferred by β-tubulin mutations: a molecular modeling study on native and mutated tubulin complexes with epothilone B. J Biomol Struct Dyn 2015; 33(12): 2530-40.[http://dx.doi.org/10.1080/07391102.2015.1063455] [PMID: 26081685]
Discovery Studio Modeling Environment, Release xx 2007.
Hess B, Kutzner C, van der Spoel D, Lindahl E. GROMACS 4: algorithms for highly efficient, load-balanced, and scalable molecular simulation. J Chem Theory Comput 2008; 4(3): 435-47.[http://dx.doi.org/10.1021/ct700301q] [PMID: 26620784]
Jorgensen WL, Tirado-Rives J. The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. J Am Chem Soc 1988; 110(6): 1657-66.[http://dx.doi.org/10.1021/ja00214a001] [PMID: 27557051]
Jorgensen WL, Chandrasekhar J, Madura JD, Impey RW, Klein ML. Comparison of simple potential functions for simulating liquid water. J Chem Phys 1983; 79: 926-35.[http://dx.doi.org/10.1063/1.445869]
Berendsen HJ, Postma JV, van Gunsteren WF, DiNola A, Haak J. Molecular dynamics with coupling to an external bath. J Chem Phys 1984; 81: 3684-90.[http://dx.doi.org/10.1063/1.448118]