Structure Based Drug Design and Molecular Docking Studies of Anticancer Molecules Paclitaxel, Etoposide and Topotecan using Novel Ligands.

BACKGROUND: Tubulin is the biochemical target for several clinically used anticancer drugs as it helps in the formation of mitotic spindle during mitosis stage of cell division. Many of the anti-cancer drugs are known to interact with tubulin and microtubules including some plant alkaloids, such as paclitaxel, etoposide and topotecan. In silico drug design of these molecules were performed prior to testing these drugs in vitro. In silico drug design of these anti-cancer drugs becomes a challenge due to the complex structure of target protein. This challenge was overcome by predicting the structure of the target protein (tubulin) by homology modeling. METHODS: In this study, computer aided drug designing approach was applied to predict the suitable docking site in target protein and the interaction of tubulin protein with paclitaxel, etoposide and topotecan was explored by molecular docking using Schrödinger software. Docking score and glide energy were determined with ligands to validate their anticancer properties. RESULTS: The results indicate that etoposide is the best drug for tubulin with a docking score of - 4.916 and glide energy of -46.470 kcal/mol compared to paclitaxel and topotecan. CONCLUSION: The testing of these drugs in silico provides an alternate to in vitro testing of these molecules on cancer cell lines which is a time and cost intensive process. The in silico study of parameters, such as docking score and glide energy, will help pharmacists in developing new molecules as targets for cancers in a time and cost-effective manner.

Thermodynamic effects of multiple protein conformations on stability and DNA binding.

The side-chain conformations of amino acids in the hydrophobic core are important for protein folding and function. A previous NMR study has shown that a mutant protein of transcriptional activator c-Myb, I155L/I181L R3, has multiple conformations and increased fluctuation in comparison with the wild type. To elucidate the quantitative correlation of structural fluctuation with stability and function, we analyzed the thermodynamic effects of I155L and I181L mutations, using R2R3 that encompasses the minimum specific DNA-binding region. Circular dichroism and differential scanning calorimetry measurements showed that the mutation of I155L had little effect on stability, while the I181L mutation significantly destabilized the protein. It is noteworthy that the decreased stability resulting from the I181L mutation was mainly due to decreased enthalpy change, which is partially compensated by decreased entropy change. Isothermal titration calorimetry measurements showed that the specific DNA-binding affinity was decreased owing to the I181L mutation, which was due to decreased binding entropy change. Entropy in the folded state, which corresponds to the DNA-free state, increases due to the I181L mutation because of the increased conformational fluctuation observed in I155L/I181L mutant of R2R3 by CLEANEX-PM NMR analysis, which in turn results in decreased folding entropy and DNA-binding entropy changes.

Conformational control of the binding of the transactivation domain of the MLL protein and c-Myb to the KIX domain of CREB.

The KIX domain of CBP is a transcriptional coactivator. Concomitant binding to the activation domain of proto-oncogene protein c-Myb and the transactivation domain of the trithorax group protein mixed lineage leukemia (MLL) transcription factor lead to the biologically active ternary MLL∶KIX∶c-Myb complex which plays a role in Pol II-mediated transcription. The binding of the activation domain of MLL to KIX enhances c-Myb binding. Here we carried out molecular dynamics (MD) simulations for the MLL∶KIX∶c-Myb ternary complex, its binary components and KIX with the goal of providing a mechanistic explanation for the experimental observations. The dynamic behavior revealed that the MLL binding site is allosterically coupled to the c-Myb binding site. MLL binding redistributes the conformational ensemble of KIX, leading to higher populations of states which favor c-Myb binding. The key element in the allosteric communication pathways is the KIX loop, which acts as a control mechanism to enhance subsequent binding events. We tested this conclusion by in silico mutations of loop residues in the KIX∶MLL complex and by comparing wild type and mutant dynamics through MD simulations. The loop assumed MLL binding conformation similar to that observed in the KIX∶c-Myb state which disfavors the allosteric network. The coupling with c-Myb binding site faded, abolishing the positive cooperativity observed in the presence of MLL. Our major conclusion is that by eliciting a loop-mediated allosteric switch between the different states following the binding events, transcriptional activation can be regulated. The KIX system presents an example how nature makes use of conformational control in higher level regulation of transcriptional activity and thus cellular events.