Evaluating the suitability of RNA intervention mechanism exerted by some flavonoid molecules against dengue virus MTase RNA capping site: a molecular docking,molecular dynamics simulation,and binding free energy study

Abstract

Dengue virus (DENV) is one of the most dangerous mosquito-borne human pathogens known to the mankind. Currently, no vaccines or standard therapy is avaliable to treate DENV infection. This makes the drug development against DENV more significant and challenging. The MTase domain of DENV RNA RdRp NS5 is a promising drug target, because this domain hosts the RNA capping process of DENV RNA to escape from human immune system. In the present study, we have analysed the RNA intervention mechanism exerted by flavoniod molecules against NS5 MTase RNA capping site by using molecular docking, molecular dynamics simulation and the binding free energy calculations. The results from the docking analysis confirmed that the RNA intervention mecanism is exerted by the quercetagetin (QGN) molecule with all necessary intermolecular interactions and high binding affinity. Notably, QGN forms strong hydrogen bonding interactions with Asn18, Leu20 and Ser150 residues and π???π stacking interaction with Phe25 residue. The apo and QGN bound NS5 MTase and QGN-NS5 MTase complex were used for MD simulation. The results of MD simulation reveal that the RMSD and RMSF values of QGN-MTase complex have increased on comparing the apo protein due to the effect of ligand binding. The binding free energy calulation includes prediction of total binding free energy of ligand-protein complex and per-residue free energy decomposition. The QGN binding to NS5 MTase affects it’s native motion, this result is found from Principal component analysis.

Communicated by Ramaswamy H. Sarma     

A computational study of binding between 3-(4-fluorophenyl)-N-((4-fluorophenyl)sulphonyl)acrylamide and tubulin

3-(4-Fluorophenyl)-N-((4-fluorophenyl)sulphonyl)acrylamide (FFSA) is a potential tubulin polymerisation inhibitor. In this article, a theoretical study of the binding between FFSA and tubulin in colchicine site was carried out by molecular docking, molecular dynamics (MD) simulation and binding free energy calculations. The docking calculations preliminarily indicate that there are three possible binding modes 1, 2 and 3; MD simulations and binding free energy calculations identify that binding mode 2 is the most favourable, with the lowest binding free energy of ? 29.54 kcal/mol. Moreover, our valuable results for the binding are as follows: the inhibitor FFSA is suitably located at the colchicine site of tubulin, where it not only interacts with residues Leu248β, Lys254β, Leu255β, Lys352β, Met259β and Val181a by hydrophilic interaction, but also interacts with Val181α and Thr179α by hydrogen bond interaction. These two factors are both essential for FFSA strongly binding to tubulin. These theoretical results help understanding the action mechanism and designing new compounds with higher affinity to tubulin.     

Study on the binding mode of a pyrrolotriazin derivative with JAK2 by docking and MD simulation

The pyrrolotriazin derivative 2-(4-(4-((7-(3-(N-methylmethylsulfonamido)phenyl)pyrrolo [2,1-f][1,2,4]triazin-2-yl)amino)phenyl)piperidin-1-yl)acetamide (PPA) is a potential Janus kinase 2 (JAK2) inhibitor. The binding mode between PPA and JAK2 was investigated by using a combined method of docking, molecular dynamics (MD) simulation and binding free-energy calculation. The docking calculations preliminarily indicated that there were two possible binding modes 1 and 2; MD simulations and binding free-energy calculations identified that binding mode 1 was more stable and favourable, with the lower MM-PBSA binding free energy of ?34.00?±?0.17?kcal/mol. Moreover, some valuable binding information is revealed as follows: the inhibitor PPA is suitably located at the ATP-binding site of JAK2 and the hydrophobic interaction plays an essential role. PPA not only interacts with residues Leu855, Val863, Ala880, Tyr931, Leu932 and Leu983 via hydrophobic interaction but also interacts with Ser936 and Asp994 by hydrogen bonds. These two factors are advantageous for PPA to strongly bind to JAK2. These results help to understand the action mechanisms and designing new compounds with a higher affinity to JAK2.     

Negative allosteric modulators of cannabinoid receptor 2: protein modeling,binding site identification and molecular dynamics simulations in the presence of an orthosteric agonist

Abstract

Selective activation of the cannabinoid receptor subtype 2 (CB2) shows promise for treating pain, inflammation, multiple sclerosis, cancer, ischemic/reperfusion injury and osteoporosis. Target selectivity and off-target side effects are two major limiting factors for orthosteric ligands, and therefore, the search for allosteric modulators (AMs) is a widely used drug discovery approach. To date, only a limited number of negative CB2 AMs have been identified, possessing only micromolar activity at best, and the CB2 receptor’s allosteric site(s) are not well characterized. Herein, we used computational approaches including receptor modeling, site mapping, docking, molecular dynamics (MD) simulations and binding free energy calculations to predict, characterize and validate allosteric sites within the complex of the CB2 receptor with bound orthosteric agonist CP55,940. After docking of known negative CB2 allosteric modulators (NAMs), dihydro-gambogic acid (DHGA) and trans-β-caryophyllene (TBC) (note that TBC also shows agonist activity), at the predicted allosteric sites, the best total complex with CB2, CP55,940 and NAM was embedded into a hydrated lipid bilayer and subjected to a 200 ns MD simulation. The presence of an AM affected the CB2–CP55,940 complex, altering the relative positioning of the toggle switch residues and promoting a strong π–π interaction between Phe117^3.36 and Trp258^6.48. Binding of either TBC or DHGA to a putative allosteric pocket directly adjacent to the orthosteric ligand reduced the binding free energy of CP55,940, which is consistent with the expected effect of a negative AM. The identified allosteric sites present immense scope for the discovery of novel classes of CB2 AMs.     

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     

Structural insight into Mycobacterium tuberculosis maltosyl transferase inhibitors: pharmacophore-based virtual screening,docking, and molecular dynamics simulations

Pharmacophore-based virtual screening, subsequent docking, and molecular dynamics (MD) simulations have been done to identify potential inhibitors of maltosyl transferase of Mycobacterium tuberculosis (mtb GlgE). Ligand and structure-based pharmacophore models representing its primary binding site (pbs) and unique secondary binding site 2 (sbs2), respectively, were constructed based on the three dimensional structure of mtb GlgE. These pharmacophore models were further used for screening of ZINC and antituberculosis compounds database (ATD). Virtually screened molecules satisfying Lipinski’s rule of five were then analyzed using docking studies and have identified 23 molecules with better binding affinity than its natural substrate, maltose. Four top scoring ligands from ZINC and ATD that either binds to pbs or sbs2 have been subjected to 10 ns each MD simulations and binding free energy calculations. Results of these studies have confirmed stable protein ligand binding. Results reported in the article are likely to be helpful in antitubercular therapeutic development research.     

Analysis of cyclosporin A and a set of analogs as inhibitors of a T. cruzi cyclophilin by docking and molecular dynamics

Cyclophilins (CyPs) are enzymes involved in protein folding. In Trypanosoma cruzi (T. cruzi), the most abundantly expressed CyP is the isoform TcCyP19. It has been shown that TcCyP19 is inhibited by the immunosuppressive drug cyclosporin A (CsA) and analogs, which also proved to have potent trypanosomicidal activity in vitro. In this work, we continue and expand a previous study on the molecular interactions of CsA, and a set of analogs modeled in complexes with TcCyP19. The modeled complexes were used to evaluate binding free energies by molecular dynamics (MD), applying the Linear Interaction Energy (LIE) method. In addition, putative binding sites were identified by molecular docking. In our analysis, the binding free energy calculations did not correlate with experimental data. The heterogeneity of the non-bonded energies and the variation in the pattern of hydrogen bonds suggest that the systems may not be suitable for the application of the LIE method. Further, the docking calculations identified two other putative binding sites with comparable scoring energies to the active site, a fact that may also explain the lack of correlation found. Kinetic experiments are needed to confirm or reject the multiple binding sites hypothesis. In the meantime, MD simulations at the alternative sites, employing other methods to compute binding free energies, might be successful at finding good correlations with the experimental data.     

An explorative study on Staphylococcus aureus MurE inhibitor: induced fit docking,binding free energy calculation,and molecular dynamics

Staphylococcus aureus MurE enzyme catalyzes the addition of l-lysine as third residue of the peptidoglycan peptide moiety. Due to the high substrate specificity and its ubiquitous nature among bacteria, MurE enzyme is considered as one of the potential target for the development of new therapeutic agents. In the present work, induced fit docking (IFD), binding free energy calculation, and molecular dynamics (MD) simulation were carried out to elucidate the inhibition potential of 2-thioxothiazolidin-4-one based inhibitor 1 against S. aureus MurE enzyme. The inhibitor 1 formed majority of hydrogen bonds with the central domain residues Asn151, Thr152, Ser180, Arg187, and Lys219. Binding free-energy calculation by MM-GBSA approach showed that van der Waals (ΔG_vdW, ?57.30?kcal/mol) and electrostatic solvation (ΔG_solv, ?36.86?kcal/mol) energy terms are major contributors for the inhibitor binding. Further, 30-ns MD simulation was performed to validate the stability of ligand–protein complex and also to get structural insight into mode of binding. Based on the IFD and MD simulation results, we designed four new compounds D1–D4 with promising binding affinity for the S. aureus MurE enzyme. The designed compounds were subjected to the extra-precision docking and binding free energy was calculated for complexes. Further, a 30-ns MD simulation was performed for D1/4C13 complex.     

Elucidation of the anticancer potential and tubulin isotype-specific interactions of β-sitosterol

Beta-sitosterol (β-SITO), a phytosterol present in many edible vegetables, has been reported to possess antineoplastic properties and cancer treatment potential. We have shown previously that it binds at a unique site (the ‘SITO-site’) compared to the colchicine binding site at the interface of α- and β-tubulin. In this study, we investigated the anticancer efficacy of β-SITO against invasive breast carcinoma using MCF-7 cells. Since ‘isotypes’ of β-tubulin show tissue-specific expression and many are associated with cancer drug resistance, using computer-assisted docking and atomistic molecular dynamic simulations, we also examined its binding interactions to all known isotypes of β-tubulin in αβ-tubulin dimer. β-SITO inhibited MCF-7 cell viability by up to 50%, compared to vehicle-treated control cells. Indicating its antimetastatic potential, the phytosterol strongly inhibited cell migration. Immunofluorescence imaging of β-SITO-treated MCF-7 cells exhibited disruption of the microtubules and chromosome organization. Far-UV circular dichroism spectra indicated loss of helical stability in tubulin when bound to β-SITO. Docking and MD simulation studies, combined with MM-PBSA and MM-GBSA calculations revealed that β-SITO preferentially binds with specific β-tubulin isotypes (β_II and β_III) in the αβ-tubulin dimer. Both these β-tubulin isotypes have been implicated in drug resistance against tubulin-targeted chemotherapeutics. Our data show the tubulin-targeted anticancer potential of β-SITO, and its potential clinical utility against β_II and β_III isotype-overexpressing neoplasms.     

Understanding microscopic binding of macrophage migration inhibitory factor with phenolic hydrazones by molecular docking, molecular dynamics simulations and free energy calculations

Macrophage migration inhibitory factor (MIF), an immunoregulatory protein, is a potential target for a number of inflammatory diseases. In the current work, the interactions between MIF and a series of phenolic hydrazones were studied by molecular docking, molecular dynamics (MD) simulations, binding free energy calculations, and binding energy decomposition analysis to determine the structural requirement for achieving favorable biological activity of phenolic hydrazones. First, molecular docking was used to predict the binding modes of inhibitors in the binding site of MIF. The good correlation between the predicted docking scores and the experimental activities shows that the binding conformations of the inhibitors in the active site of MIF are well predicted. Moreover, our results suggest that the flexibility of MIF is essential in ligand binding process. Then, MD simulations and MM/GBSA free energy calculations were employed to determine the dynamic binding process and compare the binding modes of the inhibitors with different activities. The predicted binding free energies given by MM/GBSA are not well correlated with the experimental activities for the two subsets of the inhibitors; however, for each subset, a good correlation between the predicted binding free energies and the experimental activities is achieved. The MM/GBSA free energy decomposition analysis highlights the importance of hydrophobic residues for the MIF binding of the studied inhibitors. Based on the essential factors for MIF-inhibitor interactions derived from the theoretical predictions, some derivatives were designed and the higher inhibitory activities of several candidates were confirmed by molecular docking studies. The structural insights obtained from our study are useful for designing potent inhibitors of MIF.     

Searching the conformational complexity and binding properties of HDAC6 through docking and molecular dynamic simulations

Histone deacetylases (HDACs) are a family of proteins involved in the deacetylation of histones and other non-histones substrates. HDAC6 belongs to class II and shares similar biological functions with others of its class. Nevertheless, its three-dimensional structure that involves the catalytic site remains unknown for exploring the ligand recognition properties. Therefore, in this contribution, homology modeling, 100-ns-long Molecular Dynamics (MD) simulation and docking calculations were combined to explore the conformational complexity and binding properties of the catalytic domain 2 from HDAC6 (DD2-HDAC6), for which activity and affinity toward five different ligands have been reported. Clustering analysis allowed identifying the most populated conformers present during the MD simulation, which were used as starting models to perform docking calculations with five DD2-HDAC6 inhibitors: Cay10603 (CAY), Rocilinostat (RCT), Tubastatin A (TBA), Tubacin (TBC), and Nexturastat (NXT), and then were also submitted to 100-ns-long MD simulations. Docking calculations revealed that the five inhibitors bind at the DD2-HDAC6 binding site with the lowest binding free energy, the same binding mode is maintained along the 100-ns-long MD simulations. Overall, our results provide structural information about the molecular flexibility of apo and holo DD2-HDAC6 states as well as insight of the map of interactions between DD2-HDAC6 and five well-known DD2-HDAC6 inhibitors allowing structural details to guide the drug design. Finally, we highlight the importance of combining different theoretical approaches to provide suitable structural models for structure-based drug design.     

Homology modeling,docking and structure-based virtual screening for new inhibitor identification of Klebsiella pneumoniae heptosyltransferase-III

Abstract

Klebsiella pneumoniae (K. pneumoniae) is a Gram-negative opportunistic pathogen commonly associated with hospital-acquired infections that are often resistant even to antibiotics. Heptosyltransferase (HEP) belongs to the family of glycosyltransferase-B (GT-B) and plays an important in the synthesis of lipopolysaccharides (LPS) essential for the formation of bacterial cell membrane. HEP-III participates in the transfer of heptose sugar to the outer surface of bacteria to synthesize LPS. LPS truncation increases the bacterial sensitivity to hydrophobic antibiotics and detergents, making the HEP as a novel drug target. In the present study, we report the 3D homology model of K. pneumoniae HEP-III and its structure validation. Active site was identified based on similarities with known structures using Dali server, and structure-based pharmacophore model was developed for the active site substrate ADP. The generated pharmacophore model was used as a 3D search query for virtual screening of the ASINEX database. The hit compounds were further filtered based on fit value, molecular docking, docking scores, molecular dynamics (MD) simulations of HEP-III complexed with hit molecules, followed by binding free energy calculations using Molecular Mechanics-Poisson–Boltzmann Surface Area (MM-PBSA). The insights obtained in this work provide the rationale for design of novel inhibitors targeting K. pneumoniae HEP-III and the mechanistic aspects of their binding.

Communicated by Ramaswamy H. Sarma     

Computational investigations of gram-negative bacteria phosphopantetheine adenylyltransferase inhibitors using 3D-QSAR,molecular docking and molecular dynamic simulations

Abstract

Phosphopantetheine adenylyltransferase (PPAT) has been recognized as a promising target to develop novel antimicrobial agents, which is a hexameric enzyme that catalyzes the penultimate step in coenzyme A biosynthesis. In this work, molecular modeling study was performed with a series of PPAT inhibitors using molecular docking, three-dimensional qualitative structure-activity relationship (3D-QSAR) and molecular dynamic (MD) simulations to reveal the structural determinants for their bioactivities. Molecular docking study was applied to understand the binding mode of PPAT with its inhibitors. Subsequently, 3D-QSAR model was constructed to find the features required for different substituents on the scaffolds. For the best comparative molecular field analysis (CoMFA) model, the Q² and R² values of which were calculated as 0.702 and 0.989, while they were calculated as 0.767 and 0.983 for the best comparative molecular similarity index analysis model. The statistical data verified the significance and accuracy of our 3D-QSAR models. Furthermore, MD simulations were carried out to evaluate the stability of the receptor–ligand contacts in physiological conditions, and the results were consistent with molecular docking studies and 3D-QSAR contour map analysis. Binding free energy was calculated with molecular mechanics generalized born surface area approach, the result of which coincided well with bioactivities and demonstrated that van der Waals accounted for the largest portion. Overall, our study provided a valuable insight for further research work on the recognition of potent PPAT inhibitors.

Communicated by Ramaswamy H. Sarma     

Substrate docking and molecular dynamic simulation for prediction of fungal enzymes from Trichoderma species-assisted extraction of nanocellulose from oil palm leaves

Abstract

Fungi of the Trichoderma species are valued industrial enzymes in support of the ‘zero-waste’ technology to convert agro-industrial biomass into valuable products, i.e. nanocellulose (NC). In this study, an in silico approach using substrate docking and molecular dynamic (MD) simulation was used to predict the order of which the multilayers of cellulosic polymers, i.e. lignin, hemicellulose and cellulose in oil palm leaves (OPL) are degraded by fungal enzymes, endocellulase and exocellulase. The study aimed to establish the catalytic tendencies of the enzymes to optimally degrade the cellulosic components of OPL for high yield production of NC. Energy minimized endocellulase and exocellulase models revealed satisfactory scores of PROCHECK (90.0% and 91.2%), Verify3D (97.23% and 98.85%) and ERRAT (95.24% and 91.00%) assessments. Active site prediction by blind docking, COACH meta-server and multiple sequence alignment indicated the catalytic triads for endocellulase and exocellulase were Ser116–His205–Glu249 and Ser382–Arg124–Asp385, respectively. Binding energy of endocellulase docked with hemicellulose (?6.0 ? kcal mol^?1) was the most favourable followed by lignin (?5.6 ? kcal mol^?1) and cellulose (?4.4 ? kcal mol^?1). Exocellulase, contrarily, bonded favorably with lignin (?8.7 ? kcal mol^?1), closely followed by cellulose (?8.5 ? kcal mol^?1) and hemicellulose (?8.4 ? kcal mol^?1). MDs simulations showed that interactions of complexes, endocellulase–hemicellulose and the exocellulase–cellulose being the most stable. Thus, the findings of the study successfully identified the specific actions of sugar-acting enzymes for NC production.

Communicated by Ramaswamy H. Sarma     

Exploration of potential RSK2 inhibitors by pharmacophore modelling,structure-based 3D-QSAR,molecular docking study and molecular dynamics simulation

Abstract

The p90 ribosomal s6 kinase 2 (RSK2) is a promising target because of its over expression and activation in human cancer cells and tissues. Over the last few years, significant efforts have been made in order to develop RSK2 inhibitors to treat myeloma, prostatic cancer, skin cancer and etc., but with limited success so far. In this paper, pharmacophore modelling, molecular docking study and molecular dynamics (MD) simulation have been performed to explore the novel inhibitors of RSK2. Pharmacophore models were developed by 95 molecules having pIC₅₀ ranging from 4.577 to 9.000. The pharmacophore model includes one hydrogen bond acceptor (A), one hydrogen bond donor (D), one hydrophobic feature (H) and one aromatic ring (R). It is the best pharmacophore hypothesis that has the highest correlation coefficient (R² = 0.91) and cross validation coefficient (Q² = 0.71) at 5 component PLS factor. It was evaluated using enrichment analysis and the best model was used for virtual screening. The constraints used in this study were docking score, ADME properties, binding free energy estimates and IFD Score to screen the database. Ultimately, 12 hits were identified as potent and novel RSK2 inhibitors. A 15 ns molecular dynamics (MD) simulation was further employed to validate the reliability of the docking results.     

Probing the binding mechanism of mercaptoguanine derivatives as inhibitors of HPPK by docking and molecular dynamics simulations

6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) is a promising antimicrobial target involved in the folate biosynthesis pathway. Although, the results from crystallographic studies of HPPK have attracted a great interest in the design of novel HPPK inhibitors, the mechanism of action of HPPK due to inhibitor binding remains questionable. Recently, mercaptoguanine derivatives were reported to inhibit the pyrophosphoryl transfer mechanism of Staphylococcus aureus HPPK (SaHPPK). The present study is an attempt to understand the SaHPPK-inhibitors binding mechanism and to highlight the key residues that possibly involve in the complex formation. To decipher these questions, we used the state-of-the-art advanced insilico approach such as molecular docking, molecular dynamics (MD), molecular mechanics-generalized Born surface area approach. Domain cross correlation and principle component analysis were applied to the snapshots obtained from MD revealed that the compounds with high binding affinity stabilize the conformational dynamics of SaHPPK. The binding free energy estimation showed that the van der Waals and electrostatic interactions played a vital role for the binding mechanism. Additionally, the predicted binding free energy was in good agreement with the experimental values (R² = .78). Moreover, the free energy decomposition on per-residue confirms the key residues that significantly contribute to the complex formation. These results are expected to be useful for rational design of novel SaHPPK inhibitors.     

Computational insights into the binding pattern of mitochondrial calcium uniporter inhibitor through homology modeling,molecular dynamics simulation,binding free energy prediction and density functional theory calculation

Abstract

The mitochondrial calcium uniporter (MCU) is the critical protein of the inner mitochondrial membrane that is the primary mediator for calcium uptake into the mitochondrial matrix. Herein we built the optimal homology model of human MCU which was refined through all-atom molecular dynamics simulation. Then, the binding mode of known inhibitor was predicted through molecular docking method, along with molecular dynamics simulation and binding free energy calculation to verify the docking result and stability of the protein-inhibitor complex. Finally, density functional theory (DFT) calculation enhanced our understanding of the molecular interaction of MCU inhibitor. Our research would provide a deeper insight into the interactions between human MCU and its inhibitor, which boosts to develop novel therapy against MCU related disease.

Communicated by Ramaswamy H. Sarma     

Discovery of potential sclerostin inhibitors from plants with loop2 region of sclerostin inhibition by interacting with residues outside Pro-Asn-Ala-Ile-Gly motif

Abstract

Sclerostin, an antagonist of the Wnt/β-catenin signaling pathway, was discovered as a potential therapeutic target for stimulating bone formation in osteoporosis. In this study, molecular docking was employed to predict the binding of 29 herbal compounds, which were reported as bone formation stimulators, to the loop2 region of sclerostin. Then, the 50 ns molecular dynamics (MD) simulation of the complexes between sclerostin and the top 10 hits obtained from molecular docking were carried out. Root mean square deviations (RMSDs) analysis of MD trajectories pointed out that all ligands-complexes remain stable throughout the duration of MD simulations. In addition, the molecular mechanics/generalized born surface area (MM/GBSA) binding free energy and energy decomposition analyses were determined. The results here suggested that baicalin is the most promising inhibitor of sclerostin. Interestingly, baicalin binds to sclerostin via the hydrophobic interaction with the amino acid residues on loop2 region but outside the Pro-Asn-Ala-Ile-Gly (PNAIG) motif, particularly the Arg-Gly-Lys-Trp-Trp-Arg (RGKWWR) motif. This finding could be a novel strategy for developing new sclerostin inhibitors in the future.

Communicated by Ramaswamy H. Sarma     

Structural insight into the mechanism of epothilone A bound to beta-tubulin and its mutants at Arg282Gln and Thr274Ile

Epothilone A (EpoA) is under investigation as an antitumor agent. To provide better understanding of the activity of EpoA against cancers, by theoretical studies such as using docking method, molecular dynamics simulation and density functional theory calculations, we identify several key residues located on β-tubulin as the active sites to establish an active pocket responsible for interaction with EpoA. Eight residues (Arg276, Asp224, Asp26, His227, Glu27, Glu22, Thr274, and Met363) are identified as the active sites to form the active pocket on β-tubulin. The interaction energy is predicted to be -121.3?kJ/mol between EpoA and β-tubulin. In the mutant of β-tubulin at Thr274Ile, three residues (Arg359, Glu27, and His227) are identified as the active sites for the binding of EpoA. In the mutant of β-tubulin at Arg282Gln, three residues (Arg276, Lys19, and His227) serve as the active sites. The interaction energy is reduced to -77.2?kJ/mol between EpoA and Arg282Gln mutant and to -50.2?kJ/mol between EpoA and Thr274Ile mutant. The strong interaction with β-tubulin is significant to EpoA's activity against cancer cells. When β-tubulin is mutated either at Arg282Gln or at Thr274Ile, the decreased strength of interaction explains the activity reduced for EpoA. Therefore, this work shows that the structural basis of the active pocket plays an important role in regulating the activity for EpoA with a Taxol-like mechanism of action to be promoted as an antitumor agent.