Probing the anticancer mechanism of prospective herbal drug Withaferin A on mammals: a case study on human and bovine proteasomes
© Grover et al. 2010
Published: 2 December 2010
Skip to main content
© Grover et al. 2010
Published: 2 December 2010
The UPP (ubiquitin proteasome pathway) is the major proteolytic system in the cytosol and nucleus of all eukaryotic cells which regulates cellular events, including mitotis, differentiation, signal transduction, apoptosis, and inflammation. UPP controls activation of the transcriptional factor NF-κB (nuclear factor κB), which is a regulatory protein playing central role in a variety of cellular processes including immune and inflammatory responses, apoptosis, and cellular proliferation. Since the primary interaction of proteasomes occurs with endogenous proteins, the signalling action of transcription factor NF-κB can be blocked by inhibition of proteasomes. A great variety of natural and synthetic chemical compounds classified as peptide aldehydes, peptide boronates, nonpeptide inhibitors, peptide vinyl sulfones and epoxyketones are now widely used as research tools for probing their potential to inhibit proteolytic activities of different proteasomes and to investigate the underlying inhibition mechanisms. The present work reports a bio-computational study carried out with the aim of exploring the proteasome inhibition capability of WA (withaferin A), a steroidal lactone, by understanding the binding mode of WA as a ligand into the mammalian proteasomes (X-ray crystal structure of Bos taurus 20S proteasome and multiple template homology modelled structure of 20S proteasome of Homo sapiens) using molecular docking and molecular dynamics simulation studies.
One possible mode of action which is proposed here for WA to act as a proteasome inhibitor is by suppression of the proteolytic activity which depends on the N-terminal threonine (Thr1) residue hydroxyl group. Docking studies carried out with herbal ligand WA into the structures of bovine and human proteasomes substantiate that WA has the ability to inhibit activity of mammalian 20S proteasomes by blocking the nucleophilic function of N-terminal Thr1. Results from molecular dynamics simulations in water show that the trajectories of both the native human 20S proteasome and the proteasome complexed with WA are stable over a considerably long time period of 4 ns suggesting the dynamic structural stability of human 20S proteasome/WA complex.
Inhibition of proteasomal activity are promising ways to retard or block degradation of specific proteins to correct diverse pathologies. Though quite a number of selective and efficient proteasomal inhibitors exist nowadays, their toxic side effects limit their potential in possible disease treatment. Thus there is an indispensable need for exploration of novel natural products as antitumor drug candidates. The present work supports the mammalian proteasomes inhibiting activity of WA along with elucidation of its possible mode of action. Since WA is a small herbal molecule, it is expected to provide one of the modest modes of inhibition along with added favours of ease in oral administration and decreased immunogenicity. The molecular docking results suggest that WA can inhibit the mammalian proteasomes irreversibly and with a high rate through acylation of the N-terminal Thr1 of the β-5 subunit.
Ubiquitin is a small 76 amino acid protein conserved in all eukaryotic cells with a molecular weight of 8.6 kDa. When polyubiquitin is attached to target proteins, tagged proteins are selected for destruction by cytoplasmic organelles called proteasomes . The mammalian 20S proteasome is characterized by a cylindrical shaped quaternary structure consisting of four heptameric stacked rings, α7β7β7α7, with 7 distinct α-type and 7 distinct β -type subunits, and with C2 symmetry similar to those of the yeast . Within this multienzymatic proteasome system, a proteolytic multifunctional complex 26S proteasome is involved, which consists of a 19S regulatory particle and a 20S core particle [3, 4]. The two outer α rings complex with the two 19S regulatory particles, forming a narrow channel through which only denatured proteins can pass . The catalytic chamber is formed by the two inner rings, each of which contains three well-characterized proteolytic activities. In particular, three active subunits β1, β2, and β5 are responsible for the three major peptidase activities: the peptidylglutamil-hydrolase like, trypsin-like and chymotrypsin-like (ChT-L) activities respectively, as investigated by mutational and crystallographic studies [6, 7]. The core particle is responsible for degradation of the proteins in a progressive manner, generating peptides of 3-25 amino acids in length .
The enzymatic activity of β-subunits is associated with the N-terminal threonine residues, which act as nucleophiles in the hydrolysis reaction catalyzing the cleavage of peptides through nucleophilic attack. Thus these proteasomes are classified as members of the Ntn (N-terminal nucleophilic) hydrolases group [1, 9]. In the eukaryotic proteasome, out of the seven different β-subunit precursors which are processed during particle maturation by autolysis, only β1, β2 and β5 subunits can be activated by the autolytic process, with the release of the amino-terminal Thr1 functioning as the nucleophile [10–12].
UPP is the major proteolytic system in the cytosol and nucleus of all eukaryotic cells [13, 14] which regulates cellular events including mitotis, differentiation, signal transduction, apoptosis, and inflammation . UPP controls activation of the transcriptional factor NF-κB (nuclear factor κB), which is a regulatory protein playing central role in a variety of cellular processes, including immune and inflammatory responses, apoptosis, and cellular proliferation. Since the primary interaction of proteasomes occurs with endogenous proteins, the signalling action of transcription factor NF-κB can be blocked by inhibition of proteasomes, thus inhibiting the completion of the cell cycle and mitotic proliferation of cancerous cells and ultimately leading to cell death. It has been suggested that proteasomal activity is essential for tumour cell proliferation and development of drug resistance. Therefore, the development of specific inhibitors of proteasome mediated degradation pathway is now of considerable interest in the drug discovery research for cancer therapy and prevention . A great variety of natural and synthetic chemical compounds classified as peptide aldehydes, peptide boronates, nonpeptide inhibitors, peptide vinyl sulfones, and epoxyketones are now widely used as research tools for studying their ability to inhibit proteolytic activity of various proteasomes from diverse origins.
Currently two proteasome inhibitors, bortezomib and NPI-0052, have been stated in clinical trials [17, 18]. However, undesirable side effects such as fatigue, nausea, vomiting, peripheral neuropathy, anaemia, diarrhoea, and constipation have also been reported for these drugs . Therefore, there has been an intensive drive to develop new proteasome inhibitors especially those of natural origin having little or no side effects.
The amino acid sequence of the β-5 subunit of h20S (human 20S proteasome) (GenBank: AAH57840.1) comprising of 263 amino acid residues was retrieved from NCBI. Since the crystal structure of h20S is not available, we used Modeller9v7  to perform the homology modelling. The best possible templates were obtained using the build_profile python script. We built three dimensional structure of h20S using multiple template comparative homology modelling based on the following X-ray crystal structures obtained from PDB (Protein Data Bank): b20S (bovine 20S proteasome) at 2.75 Å resolution (PDB:1IRU), yeast 20S proteasome at 2.4 Å resolution (PDB:1RYP) and 20S proteasome from Archaeoglobus fulgidus (PDB:1J2Q). These template proteins were chosen based on a significant sequence similarity of h20S with these proteins in addition to their satisfactory crystallographic resolution. The alignment files of the targets and templates were then prepared using the align2D_mult python script (Sequence alignment module in MODELLER). These alignment files along with the x-ray crystal structures of the templates were used to generate the three dimensional structure models using the model_mult python script of MODELLER.
Discovery Studio (Version 1.7, Accelrys Software Inc.) was then used to execute energy minimization and to perform stereochemical quality checks to arrive at the best possible three dimensional structure of the protein. The force field applied was CHARMm and the energy minimization algorithm used was Conjugate Gradient with an RMS gradient of 0.1 using a maximum of 2000 steps. This resulted in model structures with considerably favourable potential energies. Furthermore, the variability among the models was used to evaluate the reliability of the modelling. The qualities of these models were analyzed by PROCHECKv3.4 .
Binding Site analysis module of Discovery Studio was used to identify the putative binding pockets and protein ligand binding sites in the energy minimized three-dimensional structures of b20S and h20S.
The energy minimized crystal structure of Bos taurus b20S (PDB: 1IRU) and the modelled structure of h20S were used to carry out molecular dockings. For ease in dockings, the subunits which are not directly interacting with the β-5 domain of the protein were removed of b20s crystal structure, and only the domains K, L, M, X & Y were retained for further analysis. The ligand molecule Withaferin-A [PubChem:265237] was retrieved from NCBI-PubChem Compound database .
AutoDock 4.0 suite was used as molecular-docking tool in order to carry out the docking simulations . AutoDock 4.0 was launched in a Cygwin interface in the Windows operating system. Docking logs were analyzed in the graphical user interface of ADT (Auto dock Tools) . Water molecules were cleaned off from the protein crystal structure before docking. H-atoms were added to these target proteins for correct ionization and tautomeric states of amino acid residues and the non-polar hydrogens were then merged up. Kollman united atom charges and solvation parameters were assigned to the proteins. Gasteiger charge was assigned to the ligand and then nonpolar hydrogens were merged. Rigid roots were also assigned to the ligand and five bonds were made “active” or rotatable. The modified structures so obtained: x-ray crystal structure of bovine, modelled three dimensional structure of human 20S proteasome, and the structure of ligand WA accounting the flexibility of its bonds, were converted to PDBQT format in ADT, as required in AutoDock calculations. The Lamarckian Genetic Algorithm was used with a population size of 150 dockings. Five million energy evaluations were used in the docking experiments. All other parameters, e.g. crossover rate and mutation rate, were run with default settings. The grid size for specifying the search space was set at 40 × 30 × 30 centered on Thr1 of the β5 subunit with a default grid point spacing of 0.375 Å. Energy scoring function of AutoDock 4 is based upon the calculation of pair-wise atomic terms including evaluations for different secondary interactions, dispersion/repulsion, hydrogen bonding, electrostatics, and desolvation . Pre-calculated grid maps, which store grids of interaction energy based on the interaction of the ligand atom probes with receptor target, were obtained using AutoGrid. 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. The results are clustered into bins of similar conformations according to the cluster root mean square deviation (rmsd) and orientation.
The docking results obtained using AutoDock were also confirmed using ParDOCK , 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.
The AMBER v.10 package  was used to prepare the protein and the ligand files as well as for the MD (Molecular Dynamics) simulations. The binding complex of h20S/WA obtained using ParDOCK and the free protein simulated in this study were neutralized by adding appropriate number of chloride counterions and were solvated in a octahedron box of TIP4P water with a 10 Å distance between the protein surface and the box boundary . The partial atomic charges for the ligand were obtained using “antechamber”  module of Amber. The energy minimization and MD simulations of h20S and its complex with WA were carried out with the aid of the PMEMD module of the AMBER 10 program. First of all, the simulated binding complex was effected with a 2500 step minimization using the steepest descent algorithm followed by a 1000 step minimization using conjugate gradient to remove bad steric contacts. Topology and parameter files for the protein were generated using “ff03” and for the drug using “gaff” based on the atom types of the force field model developed by Cornell et al . Then the system was equilibrated beginning with the protein atom restrained simulations having 150 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 = 10 to 300 K in 58 small intervals of 2.5 ps each for a total period of 145 ps. The entire system was then equilibrated at 300 K for 100 ps before a sufficiently long MD simulation (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  and constant pressure P=1 atm with isotropic molecule-based scaling. The SHAKE algorithm  was applied to fix all covalent bonds containing hydrogen atoms. We used a time step of 2 fs and a nonbond-interaction cut-off radius of 10 A°. The Particle Mesh Ewald (PME) method  was used to treat long-range electrostatic interactions. The coordinates of the trajectory was sampled every 1 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 Indian Institute of Technology Delhi.
Properties of the docked conformations
Thr1 proximity to
Ester group of lactone ring
Ester group of lactone ring
Total internal energy
DOPE scores of the homology modelled structures
It has been shown using kinetic analysis and X-ray diffraction studies that the ester bond of specific proteasome inhibitor lactacystin covalently modifies the N-terminal threonine of the β5-subunit, which is critical for proteasome inhibition [1, 47–49]. Thus it is quite probable that lactacystin-like reaction occurs with WA also as it contains an internal ester bond in a δ-valero lactone ring (Figure 1B). A similar kind of cleavage of the lactone ring by serine protease has been reported in which the 3-benzyl-2-oxetanone, a β-lactone has been found to be a slowly hydrolyzed substrate of α-chymotrypsin . Other herbal ligands like polyphenols  especially green tea polyphenols like EGCG and its analogs, genistein etc. [52, 53], which have an ester bond susceptible to nucleophilic attack by Thr1 have also been reported to possess proteasome inhibition activity. Our results obtained from docking of WA into bovine and human 20S proteasome structures substantiate the proposed inhibition mechanism.
Since proteasomes play an essential role in the turnover of cellular proteins, modulation or inhibition of proteasomal activity are thus promising ways to retard or block degradation of specific proteins in order to correct diverse pathologies. Though nowadays there exist quite a number of selective and efficient proteasomal inhibitors, the toxic side effects of these compounds strongly limit their potential in possible disease treatment. Thus there is an indispensable need for exploration of novel natural products as anti-cancer drug candidates. The study conducted here makes use of molecular docking and molecular dynamics simulation approaches, which include the search in space for the energetically most favorable conformation of a protein-ligand complex and the scoring of the resulting geometries with respect to binding energy, to analyze the proteasome inhibitory potential of WA and to investigate the underlying inhibitory mechanism. We have obtained significant results delineating the mammalian proteasomes’ inhibitory activity of WA alongwith elucidation of its possible mode of action. Since WA is a small herbal molecule, it is expected to provide one of the modest modes of inhibition alogwith added favors of ease in oral administration and decreased immunogenicity. Conclusively it is strongly suggested here that WA is a potent proteasome inhibitor and should be looked forward for further clinical investigations as a possible proteasome inhibitory drug candidate.
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 Genomics Volume 11 Supplement 4, 2010: Ninth International Conference on Bioinformatics (InCoB2010): Computational Biology. The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2164/11?issue=S4.
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.