HSPs (Heat shock proteins) are highly conserved ubiquitous proteins among species which are involved in maintaining appropriate folding and conformation of other proteins and are thus referred to as molecular chaperones. Hsp90 (Heat-shock protein 90 kDa) is one of a group of molecular chaperones responsible for managing protein folding and quality control in cell environment. However it is also involved in the maturation and stabilization of a wide range of oncogenic client proteins which are crucial for oncogenesis and malignant progression. Hsp90 requires a series of co-chaperones to assemble into a super-chaperone complex for its function. These co-chaperones bind and leave the complex at various stages to regulate the chaperoning process. Arresting the chaperone cycle at these stages by targeting different co-chaperone/Hsp90 interactions seems to be quite a viable alternative and is likely to achieve similar consequences as that of Hsp90 direct inhibition with added favors of high specificity and reduced side effect profile. The study conducted here is an attempt to explore the potential of Withania somnifera’s major constituent WA (Withaferin A) in attenuating the Hsp90/Cdc37 chaperone/co-chaperone interactions for enhanced tumor arresting activity and to elucidate the underlying mode of action using computational approaches.
Results
Formation of active Hsp90/Cdc37 complex is one of the essential steps for facilitation of chaperone client interaction, non-assembly of which can lead to prevention of the chaperone-client association resulting in apoptosis of tumor cells. From our flexible docking analysis of WA into active Hsp90/Cdc37 complex in which key interfacing residues of the complex were kept flexible, disruption of the active association complex can be discerned. While docking of WA into segregated Hsp90 leaves the interface residues untouched. Thus the molecular docking analysis of WA into Hsp90 and active Hsp90/Cdc37 complex conducted in this study provides significant evidence in support of the proposed mechanism of chaperone assembly suppression by inhibition or disruption of active Hsp90/Cdc37 complex formation being accounted by non-assembly of the catalytically active Hsp90/Cdc37 complex. Results from the molecular dynamics simulations in water show that the trajectories of the protein complexed with ligand WA are stable over a considerably long time period of 4 ns, with the energies of the complex being lowered in comparison to the un-docked association complex, suggesting the thermodynamic stability of WA complexed Hsp90/Cdc37.
Conclusions
The molecular chaperone Hsp90 has been a promising target for cancer therapy. Cancer is a disease marked by genetic instability. Thus specific inhibition of individual proteins or signalling pathways holds a great potential for subversion of this genetic plasticity of cancers. This study is a step forward in this direction. Our computational analysis provided a rationalization to the ability of naturally occurring WA to alter the chaperone signalling pathway. The large value of binding energy involved in binding of WA to the active Hsp90/Cdc37 complex consolidates the thermodynamic stability of the binding. Our docking results obtained substantiate the hypothesis that WA has the potential to inhibit the association of chaperone (Hsp90) to its co-chaperone (Cdc37) by disrupting the stability of attachment of Hsp90 to Cdc37. Conclusively our results strongly suggest that withaferin A is a potent anticancer agent as ascertained by its potent Hsp90-client modulating capability.
Background
Heat shock proteins (HSPs) are highly conserved ubiquitous proteins among species which are involved in maintaining appropriate folding and conformation of other proteins and are thus referred to as molecular chaperones [1]. Hsp90 (90 kDa heat-shock protein) is one of a group of molecular chaperones responsible for managing protein folding and quality control in cell environment [2]. It is one of the most abundant proteins in the eukaryotic cells comprising 1-2% of total proteins under non-stress conditions [3]. Although Hsp90 is involved in assisting protein folding and preventing aggregation of non-native proteins, it also plays a key role under normal conditions in regulating the stability and activation state of more than 200 ‘client’ proteins molecules including transmembrane tyrosine kinasestastable signaling proteins (Akt, Raf-1 and IKK), mutated signaling proteins (p53, Kit, Flt3 and v-src), chimeric signaling proteins (NPM–ALK, Bcr–Abl), steroid receptors (androgen, estrogen and progesterone receptors) and cell-cycle regulators (cdk4, cdk6) [4]. Recent systems-biology studies indicate that Hsp90 is a major interaction node, regulating a very diverse set of cellular functions [5–7].
Hsp90 as a therapeutic target in cancer
Hsp90 is involved in the maturation and stabilization of a wide range of oncogenic client proteins which are crucial for oncogenesis and malignant progression [8], the cancer cells being particularly dependent on proper Hsp90 function [9, 10]. Many of these client proteins are mutated and/or overexpressed in cancers [11]. Moreover, the harsh environmental conditions found in tumors such as hypoxia, low pH, and bad nutritional status may tend to destabilize proteins, making them even more dependent on Hsp90 activity [8, 12, 13]. It has been observed that Hsp90 is constitutively expressed at manifold levels in tumor cells compared to their normal counterparts [14], suggesting that Hsp90 is critically important for tumor cell growth and and/or survival and its inhibition would help check the proliferation of cancer cells. Inhibition of Hsp90 would lead to the proteasomal degradation of a large number of oncogenic client proteins, offering an additional favour of combinatorial impact on multiple oncogenic pathways and amelioration of all the cancer attributing traits, including proliferation, evasion of apoptosis, immortalization, invasion, angiogenesis and metastasis [9]. In addition, this combinatorial action markedly reduces the opportunities for cancer cells to develop resistance to Hsp90 inhibition [15]. The association of Hsp90 with a plethora of pathways including signal transduction, cell-cycle control, transcriptional regulation and many others [16–18] marks it as target for pharmacological modulation in a number of diseases, particularly for cancer, where the chaperoning of mutated and overexpressed oncoproteins is critical [19]. This has led to the development of Hsp90 inhibitors for treatment of human cancers [4, 20].
Development of Hsp90 inhibitors
Chaperoning of Hsp90 client proteins is regulated through a dynamic cycle driven by ATP binding to Hsp90 and subsequent hydrolysis [21]. Thus blocking of ATP binding has developed as a single well known mechanism for development of Hsp90 inhibitors. The currently established natural product inhibitors of Hsp90 like radicicol, benzoquinone ansamycins like geldanamycin and its derivative 17-allylamino geldanamycin, which are all based on the same ATP blocking criterion, cause the catalytic cycle of Hsp90 to arrest in the ADP-bound conformation, leading to the inactivation of chaperone activity, premature ubiquitination and of client proteins [22, 23]. The proteasomal degradation of the client proteins finally result in depletion of oncoproteins, cell cycle arrest and consequent apoptosis [24]. Though blockade of ATP binding seems to be the most direct and the simplest way to inhibit Hsp90 function, the intrinsic non-selectivity of this method for Hsp90 clientele limits its further application [25]. Of all the inhibitors employing ATP blockage mechanism for Hsp90 inhibition, none has received the Food and Drug Administration approval [26], providing space for exploration of alternative inhibition mechanisms.
Blocking of co-chaperone/chaperone interactions
It has been well reported that Hsp90 requires a series of co-chaperones to assemble into a super-chaperone complex for its function. These co-chaperones bind and leave the complex at various stages to regulate the chaperoning process [21]. Arresting the chaperone cycle at these stages by targeting different co-chaperone/Hsp90 interactions seems to be quite a viable alternative and is likely to achieve similar consequences as that of Hsp90 direct inhibition [26, 27]. Cdc37 (Cell division cycle protein 37), one of the essential co-chaperone of Hsp90, plays a specialized and indispensable role in the maturation of a number of kinase clients, whose functions are commonly dysregulated in cancer, such as receptor tyrosine kinases like epidermal growth factor receptor (EGFR), non-receptor tyrosine kinases Src and lymphocyte specific protein tyrosine kinase (Lck), and intracellular serine/threonine kinases such as Raf-1and Cdk4 [28]. Cdc37 facilitates the maturation of these kinase clients by acting as an adaptor, loading these kinases onto the Hsp90 complex [29–31]. Knocking out of Cdc37 using RNAi has been shown to attenuate the kinase clients’ association to Hsp90 in human colon cancer cells, leading to reduced cell proliferation. Furthermore, combining Cdc37 silencing with the HSP90 inhibitor 17-AAG induced more extensive and sustained depletion of kinase clients and potentiated cell cycle arrest and apoptosis. These results support an essential role for Cdc37 in concert with Hsp90 in maintaining oncogenic protein kinase clients and endorse the therapeutic potential of targeting Cdc37/Hsp90 in cancer [32].
Targeting Cdc37/Hsp90 interaction represents a potential alternative to direct Hsp90 inhibition that may offer greater specificity (due to Cdc37 elevated expression in cancer) and an improved side effect profile [27]. Additionally, Cdc37 may be an attractive target in tumor types such as androgen independent prostate cancers that currently lack highly effective therapies. Recently, Celastrol, a quinone methide triterpene from Tripterygium wilfordii , has been shown to exhibit anti-pancreatic cancer activity both in vitro and in vivo through the disruption of Cdc37/Hsp90 association and the subsequent degradation of Hsp90 client proteins [26].
Withania somnifera, also known as “ashwagandha”, “Queen of Ayurveda”, “Indian ginseng”, and “winter cherry”, has been an important herb in the Ayurvedic and indigenous medical systems for more than 3,000 years [33]. Its roots have been used as herb remedy to treat a variety of ailments and to promote general wellness. It has received much attention in recent years due to the presence of a large number of alkaloids and steroidal lactones known as withanolides [34, 35]. The structure of basic skeleton of withanolides is shown in Figure 1. The principal withanolide present in the plant is WA which is known to have anti-inflammatory [36], antitumor [37], antibacterial [38], antioxidant [39], anticonvulsive [40], and immunosuppressive properties [41]. It has the potential to increase tumor sensitization to radiation and chemotherapy while reducing some of the most common side effects of these conventional therapies [42]. Most recently, WA has been shown to potentiate apoptosis of tumor cells by inhibiting Hsp90 in human pancreatic cancer cells [43], protect against UV-induced skin cancer [44] and enhance neurite regeneration and memory [45, 46]. Thus, many studies have been reported depicting the effect of WA on putting a check on human cancers, but the mechanism behind this effect is still eluding the researchers. The study conducted here is an attempt to explore the potential of Withania somnifera’s major constituent WA in attenuating the Hsp90/Cdc37 chaperone/co-chaperone interactions for enhanced tumor arresting activity and to elucidate the underlying mode of action using computational approaches.
Figure 1
Structures of withanolides. (A) Withaferin A falls under the family of compounds known as withanolides which are a group of naturally occurring C28- steroidal lactones built on an intact or rearranged ergostane framework, in which C-22 and C- 26 are appropriately oxidized to form a six-membered lactone ring. The basic structure is designated as the withanolide skeleton defined as a 22-hydroxyergostan-26-oic acid-26,22-lactone. (B) Structure of withaferin A
Methods
Ligand and Receptors
The crystal structure of the Hsp90/Cdc37 association domain [PDB: 2K5B] was obtained from the Protein Data Bank (PDB) [47]. The crystal structure contained many missing atoms which were supplemented by the repairCommands module of AutoDock. Before docking, the protein crystal structure was cleaned by removing the water molecules. H-atoms were added to these target proteins for correct ionization and tautomeric states of amino acid residues. The modified structure so obtained was used for all the flexible docking studies while only Hsp90 [PDB: 3HZ5] was used for performing semi-flexible dockings. The ligand molecule withaferin A [PubChem:265237] was retrieved from NCBI-PubChem Compound database [48]. The structure of this compound is shown in Figure 1. The energy of the ligand molecule and receptors were minimized in Steepest Descent and Conjugate Gradient methods using Amber 11 [49], the most comprehensive suite for carrying molecular dynamics simulations.
Structural aspects of Cdc37-Hsp90 complex
Hsp90 is composed of three domains: an N-terminal domain of ~25-kDa, a middle domain of ~35-kDa, and a C-terminal domain of ~10-kDa [50]. The ATP-binding site of Hsp90 is located in the N-terminal domain. The 44.5-kDa Cdc37 protein can be dissected into three domains [51, 52] : an N-terminal domain (residues 1-127 (Cdc37N), 15.5 kDa), a middle domain (residues 147- 276 (Cdc37M) 16 kDa), and a C-terminal domain (residues 283-378 (Cdc37C), 10.5 kDa). The middle domain Cdc37M is highly resistant to proteolytic digestion and was found to be the most stable domain of Cdc37 [53]. The Hsp90 binding site is located in the middle domain of Cdc37 which associates with the Nterminal domain of Hsp90 [54, 55]. The structural features of the receptor macromolecule [PDB: 2K5B] have been described in detail elsewhere [47] by the depositors of the crystal structure to the Protein Data Bank. Briefly, the receptor is an association complex of Cdc37M and Hsp90N in respectively two chains. Hsp90N density extends from residues Glu 14 to Glu 223 in chain A and the Cdc37M peptide density extends from residues Lys 148 to His 276 in chain B. The Cdc37 peptides are almost helical except for few unwound stretches ranging from 164–168. For Hsp90, residues 117, 121, 123, 124, 125, 126, and 129 have been defined as active. For Cdc37M, residues 160, 161, 164, 165, 166, 167, 168, 193, 202, 204, 205, and 208 have been defined as Active. The Cdc37M residues that appear to be on the interface region of the association complex include Lys-160, Met-164, Arg-166, Arg-167, Trp-168, Leu-205, Glu-207, Gln-208, and Ile-214. Electrostatic interactions are observed between the residues Arg-166, Arg-167, and Asp-170 of Cdc37M and Gln-133 of Hsp90N. Arg-166 and Arg-167 residues of Cdc37 peptide form intermolecular hydrogen bonds with residue Gln-133 of Hsp90N. The hydrophobic core involves residues Met-164, Trp-193, Ala-204, and Leu-205 of Cdc37M and Ala-117, Ala-121, Ala-124, Ala-126, Met-130, and Phe-134 of Hsp90N. The side chain of Glu-120 of human Hsp90 forms a hydrogen bond with the Lys-202 side chain of Cdc37. An important hydrogen bond is formed between Glu-47 of Hsp90 and Arg-167 of Cdc37, which is considered to be the key interaction that arrests the ATPase action of Hsp90 [54]. It was also observed that the amino moiety of Gln 208 of Cdc37 is rather close to the side-chain hydroxyl group of the Ser-113 residue of Hsp90, with which it may possibly form an intermolecular hydrogen bond.
Semi-flexible docking
AutoDock 4.0 suite was used as molecular-docking tool in order to carry out the docking simulations [56]. AutoDock has been found to be able to locate docking modes that are consistent with X-ray crystal structures [57, 58]. AutoDock helps to simulate interactions between substrates or drug candidates as ligands and their macromolecular receptors of known three dimensional structures, allowing ligand flexibility described to a full extent elsewhere [56]. In our docking simulations we first used only the human HSP90 X-ray crystal structure [PDB: 3HZ5] for performing semi-flexible docking, with the ligand WA made flexible while keeping the receptor macromolecule being rigid. Flexibility of the ligand helps it explore six spatial degrees of freedom for rotation and translation and an arbitrary number of torsional degrees of freedom. A random perturbation to each is applied at each time step, and the interaction energy was evaluated for the new location and conformation [59].
Flexible docking
One of the novel features of AutoDock 4.0 allows side chains in the protein as well as in the ligand to be flexible. Thus flexibility of receptor molecule was also exploited in docking studies by making use of AutoDock flexres scripts. It is worth assumable that the incoming ligand arriving at the binding sites of the two chains would try to make its own associations with the residues in order to minimize the energy of the system. To facilitate this notion, the key residues of Hsp90 and Cdc37 chains, majority of which form H-bonds with corresponding residues in the two chains as reported [47], were made flexible in order to observe their mode of interactions with the ligand and to account for the subsequent rearrangements. The peptide bonds of the amino acids chosen to be flexible were kept “inactive” or non-rotatable. The Graphical User Interface program "AutoDock Tools" was used to prepare, run, and analyze the docking simulations. Protein was prepared for docking simulations by assigning of Kollman united atom charges, solvation parameters and polar hydrogens to the receptor PDB file. Water molecules were removed from its PDB file to make it a free receptor. Since ligands are not peptides, Gasteiger charge was assigned and then nonpolar hydrogens were merged. AutoDock assigns the rigid roots to the ligand automatically saving time as compared to manual picking. Five bonds in the ligand were made “active” or rotatable. Atomic solvation parameters were assigned to the receptor using default parameters.
Grid design
Auto-Tors, an auxiliary program using Interactive queries to define rotatable torsion angles, is used to assign all rotatable dihedrals in the ligands and to remove non-polar hydrogen atoms, uniting their partial charges with their bonded carbon atoms. AutoGrid, which calculates grids of interaction energy based on the interaction of the ligand atom probes with receptor target, is used to obtain the grid maps required prior to docking. Each probe consists of an atom type present in the ligand being docked. The pre-calculated grid maps store the potential energy arising from the interaction with the macromolecule. The user defined three dimensional grid must surround the region of interest in the macromolecule, and the ligand was limited to this search space during docking. In the present study, the location and dimensions of the grid box are chosen such that it incorporates the flexible residues which are involved in binding of Cdc37 to Hsp90 for the formation of active complex. The energy scoring grid was prepared as a 36, 30 and 40 A° (x, y, and z) cube for semi-flexible docking and 60, 60, and 60 A° (x, y, and z) cube for flexible docking. The spacing between grid points was 0.375 angstroms.
The Genetic Algorithm
The Lamarckian Genetic Algorithm (LGA) was chosen to search for the best conformers. Using mathematical concepts designed to simulate the conditions influencing biological evolution, genetic algorithms are able to search conformational space by "mutating" a ligand in order to find its lowest energy conformation in the "environment" of a fixed protein. Searches driven by this energy funnelling have been shown to provide a good indication as to the optimum protein-ligand interactions and therefore conveying the structure most likely to be found in vivo[60]. The default parameters for the Lamarckian genetic algorithm [56] were used as the search protocol except for the maximum number of energy evaluations, which were changed to 2.5 million. During the docking process, a maximum of 10 conformers was considered for each compound. The population size was set to 150 and the individuals were initialized randomly. Maximum number of generations was kept as 15000, maximum number of top individuals that automatically survived set to 1 with a mutation rate of 0.02 and a crossover rate of 0.8. Step sizes for translations, quaternions and torsions were kept same as defaults. An efficient and accurate energy assessment of the ligand conformations is as important to the success of a docking simulation as the power of the search algorithm. AutoDock uses a variation on the '95 force field [49] with terms empirically determined by linear regression analysis from a set of protein-ligand complexes with known binding constants [56, 61]. Gibbs Free energy (ΔG) is calculated as a sum of six energy terms of dispersion/repulsion, hydrogen bonding, electrostatic interactions, deviation from covalent geometry, internal ligand torsional constraints, and desolvation effects.
Selection and representation of docking modes
AutoDock reports the best docking solution (lowest docked free energy) for each GA run and also performs a cluster analysis in which the total number of clusters and the rank of each docking mode (cluster rank) are reported. Docking modes were selected on the basis of two criteria: ligand’s proximity to the Asp 166 of middle Cdc37 chain and extent of the ligand’s associations with the Cdc37-Hsp90 interacting residues. For a 10 GA run there would be up to 10 total docking modes from which the lowest energy-docking mode was chosen that met the above two criteria. All the AutoDock docking runs were performed in Intel Core 2 Duo P8400 CPU @ 2.26 GHz of Sony origin, with 3 GB DDR RAM. AutoDock 4.0 was compiled and run under Windows VISTA operating system. The output from AutoDock and all modeling studies as well as images were rendered with PyMOL [62] and Accelerys ViewerLite 5.0. ViewerLite was used to calculate the distances of hydrogen bonds as measured between the hydrogen and its assumed binding partner.
Confirmation of the docking results
The docking results obtained using AutoDock were also confirmed using ParDOCK [63] , which is an all atom energy based monte carlo docking protocol. Docking using ParDOCK requires a reference complex (target protein bound to a reference ligand) and a candidate molecule along with specific mention of the centre of mass of the cavity on which the ligand is to be docked.
MD simulations in water
The AMBER v.11 package [64] was used to prepare the protein and the ligand files as well as for the Molecular Dynamics (MD) simulations. The binding complexes of Hsp90/Cdc37/WA and Hsp90/WA obtained using ParDOCK, and the un-docked proteins simulated in this study were neutralized by adding appropriate number of sodium counter-ions and were solvated in a octahedron box of TIP4PEW water with a 10 Å distance between the protein surface and the box boundary [65]. The partial atomic charges for the ligand were obtained using “antechamber” [66] module of Amber.
The energy minimization and MD simulations of Hsp90, Hsp90/Cdc37 and their complex with WA were carried out with the aid of the SANDER module of the AMBER 11 program. First of all, the simulated binding complex was effected with a 750 step minimization using the steepest descent algorithm followed by a 250 step minimization using conjugate gradient to remove bad steric contacts. Topology and parameter files for the protein were generated using “ff99SB” and for the drug using “gaff” based on the atom types of the force field model developed by Cornell et al [49]. Then the system was equilibrated beginning with the protein atom restrained simulations having 200 ps equilibration dynamics of the solvent molecules at 300 K and a harmonic potential with a 10 kcal/mol restraint force. Next step involved the equilibration of the solute molecules with a fixed configuration of the solvent molecules in which the system was slowly heated from T = 0 to 300 K in 60 intervals each involving heating for a 5 K increase in 2.5 ps followed by an equal time duration equilibration step. The entire system was then equilibrated at 300 K for 200 ps before a sufficiently long MD simulation (for 4 ns) at room temperature. The MD simulations were performed with a periodic boundary condition in the NPT ensemble at T=298.15 K with Berendsen temperature coupling [67] and constant pressure P=1 atm with isotropic molecule-based scaling . The SHAKE algorithm [68] was applied to fix all covalent bonds containing hydrogen atoms. We used a time step of 2 fs and a nonbond-interaction cutoff radius of 10 A°. The Particle Mesh Ewald (PME) method [69] was used to treat long-range electrostatic interactions. The coordinates of the trajectory was sampled every 2.5 ps for analysis of the energy stabilization and RMSD values of the protein as well as that of the complex. MD simulations were performed on a 320 processors SUN Microsystems clusters at Supercomputing Facility (SCFBio) at IIT Delhi.
Results and discussion
Semi-flexible docking of WA into H90
One possible mode of action which is proposed here for WA to act as a chaperone function inhibitor is by restriction/disintegration of the complex between Hsp90 and Cdc37. In order to explore the possibility of non-formation of the complex, we first carried out molecular docking studies with only the Hsp90 protein crystal structure. Before docking to be carried out, the structures of receptor macromolecules were minimized using Sander module of Amber. Macromolecular receptors minimized using Steepest Descent and Conjugate Gradient methods had comparatively lower potential energy values than those of the initial ones and are thus utilized further for carrying out the docking studies. Figure 2A shows the docked ligand WA into the selective HSP90 receptor. WA gets buried inside the pocket of HSP90 as depicted by mesh representation in Figure 2B. For this particular configuration the binding energy of WA with Hsp90 is -9.10 Kcal/mol with an Inhibition constant of 214.73 nM (Table 1). The binding of WA to Hsp90 is characterized by H-bonding between a terminal hydroxyl group of the ligand and the side chain carboxyl group of Asp 102. The other terminal hydroxyl group of WA is also involved in a weak H-bond with Asp 54 of Hsp90 (Figure 3A). The lengths of the two H-bonds formed are 3.78 and 5.71 Å. In the present docked structure, WA is also forming van der waals interactions with Leu 48, Asn 51, Asp 54, Ala 55, Leu 107, Ala 111, Val 136 and Phe 138 of Hsp90 (Figure 3B). These non-covalent interactions help stabilize the binding of the ligand with the macromolecule by lowering the energy. However, no interaction of WA with the proposed Cdc37 interacting residues could be observed, as is evident from Figure 4. These finding suggests that binding of WA to Hsp90 could only provide a naive hinderance to the attachment of Cdc37 to Hsp90, owing to the presence of yet available Cdc37 interacting residues of Hsp 90 upon WA binding. However, WA binding to Hsp90 may result in the formation of a deformed complex as the Hsp90/Cdc37 complex would not get assembled to its catalytically active form.
Figure 2
Docking representations of Withaferin A. (A) Ligand docked into the Hsp90 receptor cavity (B) Docked Ligand inside the pocket of Hsp90 receptor mesh
Table 1
Energies obtained after docking of withaferin A into Hsp90
Property
Quantity
Binding Energy
-9.10 Kcal/mol
Inhibition constant
214.73 nM
Intermolecular energy
-10.36 Kcal/mol
Total internal energy
-0.11 Kcal/mol
Figure 3
Interactions of docked withaferin A with Hsp90 receptor. (A) H-Bond interactions of the docked ligand with Hsp90 residues. (B) Docked withaferin A forming vdw interactions with the hydrophobic residues of Hsp90.
Figure 4
Interactions of docked withaferin A with Hsp90 receptor. The docked withaferin A does not interact with any of the Cdc37 interacting residues of Hsp90.
Flexible docking of WA into active Hsp90/Cdc37 complex
Docking of WA performed on the catalytically active human Hsp90/Cdc37 complex provided interesting results. The binding of WA to the association complex provides significant evidence in support of the proposed mechanism of chaperone/co-chaperone activation suppression by inhibition or disruption of active Hsp90/Cdc37 complex. Two clusters were obtained for Genetic Algorithm run from generation of ten models. Large negative binding energies were obtained for both the clusters as is evident from Table 2. The lowest binding energy (-13.95 Kcal/mol) was much lower than that obtained obtained from binding of WA to the segregated Hsp90 receptor. As shown in Figure 5, the hydroxyl group of WA disrupts Ser 133:Gln 208 H-bond present in the active complex by itself hydrogen bonding with both Ser 133 . One of the ring carbonyl group of the ligand is involved in H-bonding with Gln 208 also, which has been reported as one of the critical residues in complex formation [47]. The same ring lactone group forms H-bond with Arg 166 residue of Cdc 37. Arg 166 and Arg 167 residues of Cdc37 are involved in hydrogen bonding with Gln 133 of Hsp90 in the original functional complex. Thus from our flexible docking results, it is clear that the binding of WA to the interface of Hsp90/Cdc37 structurally disrupts the active complex. The absence of active chaperone/co-chaperone complex would prevent catalysis of the conformational maturation of the Hsp90 client proteins. This would ultimately lead to degradation of Hsp90 clients thus arresting their nefarious acts.
Table 2
Clustering results obtained from docking of withaferin A into Hsp90/Cdc37 complex
Receptor
No. of AutoDockclustersa, b
Cluster rankb
No. of structuresin the cluster
Lowest binding energy of cluster
Energy range within cluster
1
3
-13.95
-13.95 to -10.18
Hsp90/Cdc37 complex
2 (10)
2
7
-13.58
-13.58 to -10.11
a Number of GA runs are shown in parentheses
b Clustering is done with RMS tolerance of 5.0 Å
Figure 5
H-bonding interactions present in the docked complex of withaferin A with active Hsp90/Cdc37 complex using flexible docking. Polar hydroxyl terminal group of withaferin A captivates amino group of Gln 208 of Cdc37 to form H-bond with itself, thus disrupting the earlier present Gln 208 - Ser 113H-bond. WA also disrupts the H-bond between Gln 133-Arg 166 by forming H-bond with Arg 166 of Cdc37. WA also forms H-bonds with Asp 54 and Asp 106 of Hsp90.
The results obtained from docking of WA into the active complex, with a majority of earlier reported H-bond and hydrophobic interaction forming residues kept flexible, clearly show that WA blocks the intermolecular interactions between Hsp90 and Cdc37 at the residues which are significantly involved in formation of the active complex. The large value of binding energy (-13.95 Kcal/mol) involved in binding of WA to the complex consolidates the thermodynamic stability of the binding. These results substantiate the hypothesis that WA possess the potential to disarray the active complex by disrupting the stability of attachment of Hsp90 and Cdc37 chains, being accounted by hydrophobic and H-bond interactions.
MD simulations in water
The Hsp90/Cdc37/WA protein-drug binding complex with the binding energy of -9.68 kcal/mol obtained using ParDOCK was used for carrying out MD simulations. After the MD simulations, we calculated RMSDs between Cα of Hsp90/ Cdc37/WA complex trajectories recorded every 2.5 ps and Cα of Hsp90/Cdc37 association complex PDB structure. The RMSDs for the trajectories of Hsp90 complexed with WA were also calculated using the protein’s X-ray structure as a reference. The results in Figure 6A show that the RMSD trajectory of Hsp90/Cdc37/WA complex was less than 2.5 Å for the entire simulation and achieves a stationary phase in the later length of the simulation. It was noted that Hsp90/Cdc37/WA (complex) has an additional entity WA owing to which the RMSD of the complex is higher than the Hsp90/Cdc37 (protein) in the initial stage (0-700 ps) of simulation. However, as the simulation progresses, the ligand in association with the receptor protein tries to make interactions with the solvent molecules and thus stabilizes the associated structure; thus resulting in lowering of the RMSD trajectory in comparison to the protein. The trajectory for the total energy was also consistent throughout the simulation length. It was also found that the energy of the complex (blue) was lower than that of the protein (red) alone throughout the length of the simulation (Figure 6B), suggesting the stability of our docked protein-drug system. MD simulations were also carried out for Hsp90 in complex with WA. The consistency of RMSD and energy profiles of the complex of WA with Hsp90 (Figure 7) suggests the stability of WA even in complex with Hsp90 alone. However a larger binding energy obtained for Hsp90/Cdc37/WA complex along with the possibility of Hsp90/Cdc37 active complex formation even in WA complexed Hsp90 substantially merits the WA binding to the association complex in comparison to segregated Hsp90 binding. This supports the hypothesis that WA is quite a probable and a worth small molecule for targeting/inhibiting the chaperone/co-chaperone association complex. The simulation length used in this study was long enough to allow rearrangement of side chains of the native as well as the drug complexed receptors to find their most stable binding mode. In conclusion, the present MD simulations made clear the dynamic structural stability of Hsp90/Cdc37 in complex with the drug WA, together with the inhibitory mechanism. These results would be valuable for further designing non-covalent type inhibitors with high specificity and potent activity by alternatively targeting the Hsp90 signalling pathway in contrast to the classic Hsp90 inhibition involving blockage of ATP binding.
Figure 6
(A) Plot of root mean square deviation (RMSD) of Cα of Hsp90/Cdc37 (protein) and Hsp90/Cdc37/WA (complex). RMSDs were calculated using the initial structures as templates. For protein (red) the reference is the PDB structure and for complex (blue) the reference is the initial docked structure. The trajectories were captured every 2.5 ps until the simulation time reached 4000 ps. (B) Plot of total energy of Hsp90/Cdc37 and Hsp90/Cdc37/WA (complex) The energy trajectories of both the protein (red) and the complex (blue) are stable over the entire length of simulation time.
Figure 7
(A) Plot of root mean square deviation (RMSD) of Cα of Hsp90 (protein) and Hsp90/ WA (complex). RMSDs were calculated using the initial structures as templates. For protein (red) the reference is the PDB structure and for complex (blue) the reference is the initial docked structure. The trajectories were captured every 2.5 ps until the simulation time reached 3200 ps. (B) Plot of total energy of Hsp90 and Hsp90/WA (complex) The energy trajectories of both the protein (red) and the complex (blue) are stable over the entire length of simulation time.
Conclusions
The molecular chaperone Hsp90 has been a promising target for cancer therapy. Cancer is a disease marked by genetic instability. Thus specific inhibition of individual proteins or signalling pathways holds a great potential for subversion of this genetic plasticity of cancers. This study is a step forward in this direction. Our computational analysis provided a rationalization to the ability of naturally occurring WA to alter the chaperone signalling pathway. The large value of binding energy involved in binding of WA to the active Hsp90/Cdc37 complex consolidates the thermodynamic stability of the binding. Our docking results obtained substantiate the hypothesis that WA has the potential to inhibit the association of chaperone (Hsp90) to its co-chaperone (Cdc37) by disrupting the stability of attachment of Hsp90 to Cdc37. Conclusively our results strongly suggest that withaferin A is a potent anticancer agent as ascertained by its potent Hsp90-client modulating capability.
Declarations
Acknowledgements
We acknowledge the support and encouragement of Prof. B. Jayaram, Coordinator, Supercomputing Facility for Bioinformatics and Computational Biology (SCFBio) at Indian Institute of Technology Delhi. We also thank Sajeev Chacko of School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi for helpful discussions. Research in the laboratory of DS is supported by grants from Department of Biotechnology (DBT) and Department of Information Technology (DIT), Government of India, New Delhi, India.
This article has been published as part of BMC Bioinformatics Volume 12 Supplement 1, 2011: Selected articles from the Ninth Asia Pacific Bioinformatics Conference (APBC 2011). The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2105/12?issue=S1.
Authors’ Affiliations
(1)
Department of Biochemical Engineering and Biotechnology, Indian Institute of Technology (IIT) Delhi
(2)
Supercomputing Facility for Bioinformatics and Computational Biology, Indian Institute of Technology (IIT) Delhi
References
Hartl FU: Molecular chaperones in cellular protein folding. Nature 1996, 381(6583):571–580. 10.1038/381571a0View ArticlePubMedGoogle Scholar
Young JC, Agashe VR, Siegers K, Hartl FU: Pathways of chaperone-mediated protein folding in the cytosol. Nat Rev Mol Cell Bio 2004, 5(10):781–791. 10.1038/nrm1492View ArticleGoogle Scholar
Welch WJ: Heat-Shock Proteins Functioning as Molecular Chaperones - Their Roles in Normal and Stressed Cells. Philos T Roy Soc B 1993, 339(1289):327–333. 10.1098/rstb.1993.0031View ArticleGoogle Scholar
Kamal A, Boehm MF, Burrows FJ: Therapeutic and diagnostic implications of Hsp90 activation. Trends Mol Med 2004, 10(6):283–290. 10.1016/j.molmed.2004.04.006View ArticlePubMedGoogle Scholar
Kamal A, Thao L, Sensintaffar J, Zhang L, Boehm MF, Fritz LC, Burrows FJ: A high-affinity conformation of Hsp90 confers tumor selectivity on Hsp90 inhibitors. Clin Cancer Res 2003, 9(16):6126s-6126s.Google Scholar
Falsone SF, Gesslbauer B, Tirk F, Piccinini AM, Kungl AJ: A proteomic snapshot of the human heat shock protein 90 interactome. Febs Lett 2005, 579(28):6350–6354. 10.1016/j.febslet.2005.10.020View ArticlePubMedGoogle Scholar
Zhao R, Davey M, Hsu YC, Kaplanek P, Tong A, Parsons AB, Krogan N, Cagney G, Mai D, Greenblatt J, et al.: Navigating the chaperone network: an integrative map of physical, genetic, and chemical-genetic interactions mediated by the yeast Hsp90 chaperone system. Febs J 2005, 272: 349–349.Google Scholar
Pearl LH, Prodromou C, Workman P: The Hsp90 molecular chaperone: an open and shut case for treatment. Biochem J 2008, 410: 439–453. 10.1042/BJ20071640View ArticlePubMedGoogle Scholar
Didelot C, Lanneau D, Brunet M, Joly AL, De Thonel A, Chiosis G, Garrido C: Anti-cancer therapeutic approaches based on intracellular and extracellular heat shock proteins. Curr Med Chem 2007, 14(27):2839–2847. 10.2174/092986707782360079View ArticlePubMedGoogle Scholar
Solit DB, Chiosis G: Development and application of Hsp90 inhibitors. Drug Discov Today 2008, 13(1–2):38–43. 10.1016/j.drudis.2007.10.007View ArticlePubMedGoogle Scholar
Gress TM, Mullerpillasch F, Weber C, Lerch MM, Friess H, Buchler M, Beger HG, Adler G: Differential Expression of Heat-Shock Proteins in Pancreatic-Carcinoma. Cancer Res 1994, 54(2):547–551.PubMedGoogle Scholar
Ferrarini M, Heltai S, Zocchi MR, Rugarli C: Unusual Expression and Localization of Heat-Shock Proteins in Human Tumor-Cells. Int J Cancer 1992, 51(4):613–619. 10.1002/ijc.2910510418View ArticlePubMedGoogle Scholar
Workman P: Combinatorial attack on multistep oncogenesis by inhibiting the Hsp90 molecular chaperone. Cancer Lett 2004, 206(2):149–157. 10.1016/j.canlet.2003.08.032View ArticlePubMedGoogle Scholar
McClellan AJ, Xia Y, Deutschbauer AM, Davis RW, Gerstein M, Frydman J: Diverse cellular functions of the Hsp90 molecular chaperone uncovered using systems approaches. Cell 2007, 131(1):121–135. 10.1016/j.cell.2007.07.036View ArticlePubMedGoogle Scholar
Isaacs JS, Xu WP, Neckers L: Heat shock protein 90 as a molecular target for cancer therapeutics. Cancer Cell 2003, 3(3):213–217. 10.1016/S1535-6108(03)00029-1View ArticlePubMedGoogle Scholar
Pratt WB, Toft DO: Regulation of signaling protein function and trafficking by the hsp90/hsp70-based chaperone machinery. Exp Biol Med 2003, 228(2):111–133.Google Scholar
Whitesell L, Lindquist SL: HSP90 and the chaperoning of cancer. Nat Rev Cancer 2005, 5(10):761–772. 10.1038/nrc1716View ArticlePubMedGoogle Scholar
Powers MV, Workman P: Inhibitors of the heat shock response: Biology and pharmacology. Febs Lett 2007, 581(19):3758–3769. 10.1016/j.febslet.2007.05.040View ArticlePubMedGoogle Scholar
Neckers L: Development of small molecule Hsp90 inhibitors: Utilizing both forward and reverse chemical genomics for drug identification. Curr Med Chem 2003, 10(9):733–739. 10.2174/0929867033457818View ArticlePubMedGoogle Scholar
Whitesell L, Mimnaugh EG, Decosta B, Myers CE, Neckers LM: Inhibition of Heat-Shock Protein Hsp90-Pp60(V-Src) Heteroprotein Complex-Formation by Benzoquinone Ansamycins - Essential Role for Stress Proteins in Oncogenic Transformation. P Natl Acad Sci USA 1994, 91(18):8324–8328. 10.1073/pnas.91.18.8324View ArticleGoogle Scholar
Sharma SV, Agatsuma T, Nakano H: Targeting of the protein chaperone, HSP90, by the transformation suppressing agent, radicicol. Oncogene 1998, 16(20):2639–2645. 10.1038/sj.onc.1201790View ArticlePubMedGoogle Scholar
Hostein I, Robertson D, DiStefano F, Workman P, Clarke PA: Inhibition of signal transduction by the Hsp90 inhibitor 17-allylamino-17-demethoxygeldanamycin results in cytostasis and apoptosis. Cancer Res 2001, 61(10):4003–4009.PubMedGoogle Scholar
Li YY, Zhang T, Schwartz SJ, Sun DX: New developments in Hsp90 inhibitors as anti-cancer therapeutics: Mechanisms, clinical perspective and more potential. Drug Resist Update 2009, 12(1–2):17–27. 10.1016/j.drup.2008.12.002View ArticleGoogle Scholar
Zhang T, Hamza A, Cao XH, Wang B, Yu SW, Zhan CG, Sun DX: A novel Hsp90 inhibitor to disrupt Hsp90/Cdc37 complex against pancreatic cancer cells. Mol Cancer Ther 2008, 7(1):162–170. 10.1158/1535-7163.MCT-07-0484View ArticlePubMedGoogle Scholar
Gray PJ, Stevenson MA, Calderwood SK: Targeting Cdc37 inhibits multiple signaling pathways and induces growth arrest in prostate cancer cells. Cancer Res 2007, 67(24):11942–11950. 10.1158/0008-5472.CAN-07-3162View ArticlePubMedGoogle Scholar
Gray PJ, Prince T, Cheng J, Stevenson MA, Calderwood SK: Targeting the oncogene and kinome chaperone CDC37. Nat Rev Cancer 2008, 8(7):491–495. 10.1038/nrc2420PubMed CentralView ArticlePubMedGoogle Scholar
Pearl LH: Hsp90 and Cdc37 - a chaperone cancer conspiracy. Curr Opin Genet Dev 2005, 15(1):55–61. 10.1016/j.gde.2004.12.011View ArticlePubMedGoogle Scholar
Vaughan CK, Gohlke U, Sobott F, Good VM, Ali MMU, Prodromou C, Robinson CV, Saibil HR, Pearl LH: Structure of an Hsp90-Cdc37-Cdk4 complex. Mol Cell 2006, 23(5):697–707. 10.1016/j.molcel.2006.07.016View ArticlePubMedGoogle Scholar
Silverstein AM, Grammatikakis N, Cochran BH, Chinkers M, Pratt WB: P50(cdc37) binds directly to the catalytic domain of Raf as well as to a site on hsp90 that is topologically adjacent to the tetratricopeptide repeat binding site. J Biol Chem 1998, 273(32):20090–20095. 10.1074/jbc.273.32.20090View ArticlePubMedGoogle Scholar
Smith JR, Clarke PA, de Billy E, Workman P: Silencing the cochaperone CDC37 destabilizes kinase clients and sensitizes cancer cells to HSP90 inhibitors. Oncogene 2009, 28(2):157–169. 10.1038/onc.2008.380PubMed CentralView ArticlePubMedGoogle Scholar
Winters M: Ancient medicine, modern use: Withania somnifera and its potential role in integrative oncology. Alternative Medicine Review 2006, 11(4):269–277.PubMedGoogle Scholar
Matsuda H, Murakami T, Kishi A, Yoshikawa M: Structures of withanosides I, II, III, IV, V, VI, and VII, new withanolide glycosides, from the roots of Indian Withania somnifera Dunal. and inhibitory activity for tachyphylaxis to clonidine in isolated guinea-pig ileum. Bioorgan Med Chem 2001, 9(6):1499–1507. 10.1016/S0968-0896(01)00024-4View ArticleGoogle Scholar
Ray A, Gupta M: Withasteroids, a growing group of naturally occurring steroidal lactones. In Progress in the chemistry of natural organic products. Volume 63. Edited by: W, Kerby G, Moore R, Steglich W, Tamm C. New York: Springer-Verlag; 1994:1–106.Google Scholar
Mishra L, Singh B, Dagenias S: Scientific basis for the therapeutic use of Withania somnifera (ashwagandha): a review. Altern Med Rev 2000, 5: 334–336.PubMedGoogle Scholar
Owais M, Sharad KS, Shehbaz A, Saleemuddin M: Antibacterial efficacy of Withania somnifera (ashwagandha) an indigenous medicinal plant against experimental murine salmonellosis. Phytomedicine 2005, 12(3):229–235. 10.1016/j.phymed.2003.07.012View ArticlePubMedGoogle Scholar
Bhattacharya A, Ghosal S, Bhattacharya SK: Anti-oxidant effect of Withania somnifera glycowithanolides in chronic footshock stress-induced perturbations of oxidative free radical scavenging enzymes and lipid peroxidation in rat frontal cortex and striatum. J Ethnopharmacol 2001, 74(1):1–6. 10.1016/S0378-8741(00)00309-3View ArticlePubMedGoogle Scholar
Kulkarni SK, George B, Mathur R: Protective effect of Withania somnifera root extract on electrographic activity in a lithium-pilocarpine model of status epilepticus. Phytotherapy Research 1998, 12(6):451–453. Publisher Full Text 10.1002/(SICI)1099-1573(199809)12:6<451::AID-PTR328>3.0.CO;2-CView ArticleGoogle Scholar
Furmanowa M, Gajdzis-Kuls D, Ruszkowska J, Czarnocki Z, Obidoska G, Sadowska A, Rani R, Upadhyay SN: In vitro propagation of Withania somnifera and isolation of withanolides with immunosuppressive activity. Planta Med 2001, 67(2):146–149. 10.1055/s-2001-11494View ArticlePubMedGoogle Scholar
Sharada AC, Solomon FE, Devi PU, Udupa N, Srinivasan KK: Antitumor and radiosensitizing effects of withaferin a on mouse Ehrlich ascites carcinoma in vivo. Acta Oncol 1996, 35(1):95–100. 10.3109/02841869609098486View ArticlePubMedGoogle Scholar
Yu YK, Hamza A, Zhang T, Gu MC, Zou P, Newman B, Li YY, Gunatilaka AAL, Zhan CG, Sun DX: Withaferin A targets heat shock protein 90 in pancreatic cancer cells. Biochem Pharmacol 2010, 79(4):542–551. 10.1016/j.bcp.2009.09.017PubMed CentralView ArticlePubMedGoogle Scholar
Mathur S, Kaur P, Sharma M, Katyal A, Singh B, Tiwari M, Chandra R: The treatment of skin carcinoma, induced by UVB radiation, using 1-oxo-5 beta,6 beta-epoxy-witha-2-enolide, isolated from the roots of Withania somnifera, in a rat model. Phytomedicine 2004, 11(5):452–460. 10.1016/j.phymed.2003.05.004View ArticlePubMedGoogle Scholar
Kuboyama T, Tohda C, Komatsu K: Withanoside IV and its active metabolite, sominone, attenuate A beta(25–35)-induced neurodegeneration. Eur J Neurosci 2006, 23(6):1417–1426. 10.1111/j.1460-9568.2006.04664.xView ArticlePubMedGoogle Scholar
Tohda C, Kuboyama T, Komatsu K: Search for natural products related to regeneration of the neuronal network. Neurosignals 2005, 23(1–2):34–45. 10.1159/000085384View ArticleGoogle Scholar
Sreeramulu S, Jonker HRA, Langer T, Richter C, Lancaster CRD, Schwalbe H: The Human Cdc37.Hsp90 Complex Studied by Heteronuclear NMR Spectroscopy. J Biol Chem 2009, 284(6):3885–3896. 10.1074/jbc.M806715200View ArticlePubMedGoogle Scholar
Cornell WD, Cieplak P, Bayly CI, Gould IR, Merz KM, Ferguson DM, Spellmeyer DC, Fox T, Caldwell JW, Kollman PA: A second generation force field for the simulation of proteins, nucleic acids, and organic molecules. J Am Chem Soc 1996, 118(9):2309–2309. 10.1021/ja955032eView ArticleGoogle Scholar
Stebbins CE, Russo AA, Schneider C, Rosen N, Hartl FU, Pavletich NP: Crystal structure of an Hsp90-geldanamycin complex: Targeting of a protein chaperone by an antitumor agent. Cell 1997, 89(2):239–250. 10.1016/S0092-8674(00)80203-2View ArticlePubMedGoogle Scholar
MacLean M, Picard D: Cdc37 goes beyond Hsp90 and kinases. Cell Stress Chaperon 2003, 8(2):114–119. Publisher Full Text 10.1379/1466-1268(2003)008<0114:CGBHAK>2.0.CO;2View ArticleGoogle Scholar
Shao JY, Irwin A, Hartson SD, Malts RL: Functional dissection of Cdc37: Characterization of domain structure and amino acid residues critical for protein kinase binding. Biochemistry-Us 2003, 42(43):12577–12588. 10.1021/bi035138jView ArticleGoogle Scholar
Zhang W, Hirshberg M, McLaughlin SH, Lazar GA, Grossmann JG, Nielsen PR, Sobott F, Robinson CV, Jackson SE, Laue ED: Biochemical and structural studies of the interaction of Cdc37 with Hsp90. J Mol Biol 2004, 340(4):891–907. 10.1016/j.jmb.2004.05.007View ArticlePubMedGoogle Scholar
Roe SM, Ali MMU, Meyer P, Vaughan CK, Panaretou B, Piper PW, Prodromou C, Pearl LH: The mechanism of Hsp90 regulation by the protein kinase-specific cochaperone p50(cdc37). Cell 2004, 116(1):87–98. 10.1016/S0092-8674(03)01027-4View ArticlePubMedGoogle Scholar
Prodromou C, Panaretou B, Chohan S, Siligardi G, O'Brien R, Ladbury JE, Roe SM, Piper PW, Pearl LH: The ATPase cycle of Hsp90 drives a molecular 'clamp' via transient dimerization of the N-terminal domains. Embo J 2000, 19(16):4383–4392. 10.1093/emboj/19.16.4383PubMed CentralView ArticlePubMedGoogle Scholar
Morris GM, Goodsell DS, Halliday RS, Huey R, Hart WE, Belew RK, Olson AJ: Automated docking using a Lamarckian genetic algorithm and an empirical binding free energy function. J Comput Chem 1998, 19(14):1639–1662. Publisher Full Text 10.1002/(SICI)1096-987X(19981115)19:14<1639::AID-JCC10>3.0.CO;2-BView ArticleGoogle Scholar
Dym O, Xenarios I, Ke HM, Colicelli J: Molecular docking of competitive phosphodiesterase inhibitors. Mol Pharmacol 2002, 61(1):20–25. 10.1124/mol.61.1.20View ArticlePubMedGoogle Scholar
Rao MS, Olson AJ: Modelling of Factor Xa-inhibitor complexes: A computational flexible docking approach. Proteins 1999, 34(2):173–183. 10.1002/(SICI)1097-0134(19990201)34:2<173::AID-PROT3>3.0.CO;2-FView ArticlePubMedGoogle Scholar
Goodsell DS, Morris GM, Olson AJ: Automated docking of flexible ligands: Applications of AutoDock. J Mol Recognit 1996, 9(1):1–5. Publisher Full Text 10.1002/(SICI)1099-1352(199601)9:1<1::AID-JMR241>3.0.CO;2-6View ArticlePubMedGoogle Scholar
Hetenyi C, van der Spoel D: Efficient docking of peptides to proteins without prior knowledge of the binding site. Protein Sci 2002, 11(7):1729–1737. 10.1110/ps.0202302PubMed CentralView ArticlePubMedGoogle Scholar
Morris GM, Goodsell DS, Huey R, Olson AJ: Distributed automated docking of flexible ligands to proteins: Parallel applications of AutoDock 2.4. J Comput Aid Mol Des 1996, 10(4):293–304. 10.1007/BF00124499View ArticleGoogle Scholar
DeLano W: The PyMOL Molecular Graphics System 2002. San Carlos, CA: DeLano Scientific; 2002.Google Scholar
Gupta A, Gandhimathi A, Sharma P, Jayaram B: ParDOCK: An all atom energy based Monte Carlo docking protocol for protein-ligand complexes. Protein Peptide Lett 2007, 14(7):632–646. 10.2174/092986607781483831View ArticleGoogle Scholar
Case DA DT, Cheatham TE, Simmerling CL, Wang J, Duke RE, Luo R, Walker RC, Zhang W, Merz KM, Roberts B, Wang B, Hayik S, Roitberg A, Seabra G, Kolossváry I, Wong IF, Paesani F, Vanicek J, Wu X, Brozell SR, Steinbrecher T, Gohlke H, Cai Q, Ye X, Wang J, Hsieh MJ, Cui G, Roe DR, Mathews DH, Seetin MG, Sagui C, Babin V, Luchko T, Gusarov S, Kovalenko A, Kollman PA: AMBER 11. San Francisco: University of California; 2010.Google Scholar
Jorgensen WL, Chandrasekhar J, Madura JD, Impey RW, Klein ML: Comparison of Simple Potential Functions for Simulating Liquid Water. J Chem Phys 1983, 79(2):926–935. 10.1063/1.445869View ArticleGoogle Scholar
Jakalian A, Bush BL, Jack DB, Bayly CI: Fast, efficient generation of high-quality atomic Charges. AM1-BCC model: I. Method. J Comput Chem 2000, 21(2):132–146. Publisher Full Text 10.1002/(SICI)1096-987X(20000130)21:2<132::AID-JCC5>3.0.CO;2-PView ArticleGoogle Scholar
Berendsen HJC, Postma JPM, Vangunsteren WF, Dinola A, Haak JR: Molecular-Dynamics with Coupling to an External Bath. J Chem Phys 1984, 81(8):3684–3690. 10.1063/1.448118View ArticleGoogle Scholar
Ryckaert JP, Ciccotti G, Berendsen HJC: Numerical-Integration of Cartesian Equations of Motion of a System with Constraints - Molecular-Dynamics of N-Alkanes. J Comput Phys 1977, 23(3):327–341. 10.1016/0021-9991(77)90098-5View ArticleGoogle Scholar
Essmann U, Perera L, Berkowitz ML, Darden T, Lee H, Pedersen LG: A Smooth Particle Mesh Ewald Method. J Chem Phys 1995, 103(19):8577–8593. 10.1063/1.470117View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
