Next Article in Journal
The Effect of Size on Ag Nanosphere Toxicity in Macrophage Cell Models and Lung Epithelial Cell Lines Is Dependent on Particle Dissolution
Next Article in Special Issue
Potential Activity of Fevicordin-A from Phaleria macrocarpa (Scheff) Boerl. Seeds as Estrogen Receptor Antagonist Based on Cytotoxicity and Molecular Modelling Studies
Previous Article in Journal
Effects of RNAi-Mediated Knockdown of Histone Methyltransferases on the Sex-Specific mRNA Expression of Imp in the Silkworm Bombyx mori
Previous Article in Special Issue
Towards Personalized Medicine Mediated by in Vitro Virus-Based Interactome Approaches
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 15
Issue 4
10.3390/ijms15046797
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling and Docking Studies on Novel Mutants (K71L and T204V) of the ATPase Domain of Human Heat Shock 70 kDa Protein 1

by
Asita Elengoe
*,
Mohammed Abu Naser
† and
Salehhuddin Hamdan
†
Faculty of Bioscience and Medical Engineering, Universiti Teknologi Malaysia, Skudai, Johor 81310, Malaysia
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Int. J. Mol. Sci. 2014, 15(4), 6797-6814; https://doi.org/10.3390/ijms15046797
Submission received: 25 January 2014 / Revised: 3 April 2014 / Accepted: 4 April 2014 / Published: 22 April 2014
(This article belongs to the Special Issue Proteins and Protein-Ligand Interactions)
Browse Figures
Versions Notes

Abstract

:
The purpose of exploring protein interactions between human adenovirus and heat shock protein 70 is to exploit a potentially synergistic interaction to enhance anti-tumoral efficacy and decrease toxicity in cancer treatment. However, the protein interaction of Hsp70 with E1A32 kDa of human adenovirus serotype 5 remains to be elucidated. In this study, two residues of ATPase domain of human heat shock 70 kDa protein 1 (PDB: 1 HJO) were mutated. 3D mutant models (K71L and T204V) using PyMol software were then constructed. The structures were evaluated by PROCHECK, ProQ, ERRAT, Verify 3D and ProSA modules. All evidence suggests that all protein models are acceptable and of good quality. The E1A32 kDa motif was retrieved from UniProt (P03255), as well as subjected to docking interaction with NBD, K71L and T204V, using the Autodock 4.2 program. The best lowest binding energy value of −9.09 kcal/mol was selected for novel T204V. Moreover, the protein-ligand complex structures were validated by RMSD, RMSF, hydrogen bonds and salt bridge analysis. This revealed that the T204V-E1A32 kDa motif complex was the most stable among all three complex structures. This study provides information about the interaction between Hsp70 and the E1A32 kDa motif, which emphasizes future perspectives to design rational drugs and vaccines in cancer therapy.

1. Introduction

The current study explores the interaction between human heat shock protein 70 (Hsp70) and adenovirus. Hyperthermia has been explored as an anti-cancer agent for many decades. While the treatment effects of hyperthermia as a single agent are limited, its ability to potentiate the effects of standard chemo radiotherapies has generated lasting interest. Despite the fact that hyperthermia in combination with adenoviral therapy has shown their effectiveness in vitro, currently no clinical trials for this combination treatment is underway. Although these oncolytic adenoviruses are promising as anticancer agents, clinical experiences show that these agents alone failed to generate sustained clinical responses or to cause complete tumor regressions. This is because heterogeneity or indeed lack of expression of receptors (coxsackie adenovirus receptor, CAR) and co-receptors (integrin αvβ3 and αvβ5 classes) in tumors could be implicated in the poor efficiency of infectivity by adenovirus. In addition, many tumors cells fail to support adenovirus replication because of its replication deficiency. Traditionally, elevated temperature as a part of the febrile response of humans was thought to inhibit viral replication. However, several investigators suggest that hyperthermia may in fact enhance viral replication, particularly in tumor cells [1]. The E1A gene products of adenovirus are responsible for activation of the Hsp70 and may affect Hsp70 levels during the cell cycle [2]. Hsp70 enhances viral proteins import and colocalizes in the nucleus with E1A. Bacterial DNAJ and DNAK, which are important for bacteriophage DNA replication, may depend on Hsp70 induction, as documented by Wickner et al. in 1992 [3]. Hsp40 and Hsp70 induction promotes production of viral proteins for avian adenovirus CELO [4]. Thus, the purpose of investigating the protein interaction between adenovirus and Hsp70 is to exploit a potentially synergistic interaction to enhance anti-tumoral efficacy and decrease toxicity in cancer therapy; but the molecular interaction between Hsp70 with E1A32 kDa of human adenovirus serotype 5 remains to be elucidated. E1A32 kDa protein consists of 289 amino acids. According to experimental studies, mutation of residues 154, 157, 171 and 174 from cysteine to serine causes a loss of E1A transactivation. In addition, mutation of residue 115 (leucine) to alanine causes complete loss of interaction with host ZMYND11, while mutation of leucine (residue 122) to isoleucine abolishes binding to UBE2I.
Hsp70 and other members of the Hsp70 protein family are described as molecular chaperones which help in the non-covalent folding or unfolding and the assembly or disassembly of other protein structures [5–8]. Human Hsp70 contains 640 amino acids and has two major domains. The 44 kDa, 380 amino acid N-terminal domain binds and hydrolyzes ATP, whereas the C-terminal domain is required for binding peptides and folding non-native polypeptides. Two functionally relevant subdomains can be differentiated within the substrate binding domain (SBD), an 18 kDa peptide-binding domain and a 10 kDa C-terminal domain that contains the Glu-Glu-Val-Asp (EEVD) regulatory motif [9]. It has recently been shown that the EEVD motif of the human Hsp70 molecular chaperone regulates ATP hydrolysis [9], and that it interacts with substrates and the co-chaperone, HDJ-1 [10]. Allosteric interaction with the SBD and the interactions with co-chaperones such as Hsp40 and nucleotide exchange factors all critically depend on the conformation of the nucleotide binding domain (NBD), thus affecting the chaperone function of Hsp70.
Since the structure of Hsp70 is yet to be determined experimentally, development of computational three-dimensional model structures of the protein from the sequence and molecular dynamic simulations on the model structure might help in understanding the dynamics of conformational changes and hence the functional mechanism of the NBDs of Hsp70. Understanding the conformational dynamics of the NBD is the key to understanding how the ATPase motor drives the periodic client binding and releasing of the Hsp70 machine. In this study, K71 and T204 of ATPase domain of human heat shock 70 kDa protein 1 (PDB: 1 HJO) involved in playing a vital role in catalytic activity were mutated. Then, we created 3D mutant structures using PyMol software. Furthermore, the docking program (Autodock Version 4.2) has been used to predict the preferred sites of interaction between the NBD, K71L and T204V mutant models with E1A32 kDa motif.

2. Results and Discussion

2.1. In-Silico Mutagenesis

Residues 71 and 204 of the NBD protein mutated from lysine (basic polar) to leucine (non-polar) and threonine (polar) to valine (non-polar) respectively. These mutations do not change the overall electronic nature of the side chains. The classification of amino acid chemical properties was based on the research done by Biro et al. (2003) [11]. Proteins that contain changes in residues may have some effects on the overall structure or function of the protein. Therefore, all mutants were chosen for the subsequent in-silico based modeling.
The conserved LYS71 is a catalytically important residue that affects ATP hydrolysis [12]. The proposed mechanism of ATP hydrolysis suggested that the role of LYS71 in accepting a proton from the hydroxide ion or water molecule involved is in-line with a nucleophilic attack [12–15]. The inorganic phosphate group (Pi) is coordinated by a salt bridge with LYS71, hydrogen bonds to THR13 and THR204 and interacts directly with a calcium ion. A water molecule mediates additional interactions with the protein’s main chain at positions 202, 203 and 204. The Pi-binding site is on the protein face opposite the highly conserved GLY32 loop that has been implicated in the binding of nucleotide release factor (GrpE) to the ATPase domain of Hsp40 (DNAK) [16]. Therefore, there are potential channels for Pi exit to the protein surface. However, release of the inorganic phosphate group has been implicated in the conformational transition of Hsp70 molecular chaperone [17]. Phospho-threonine was postulated as an intermediate of ATP hydrolysis. In addition, ATPase activity of Hsp70 initiates viral DNA replication. This has been demonstrated for bacterial DNAJ which stimulates ATPase activity of Hsp70 to start DNA replication of SV40 [18]. Thus, mutational study of these two residues which are important for ATPase activity was carried out.

2.2. Physiochemical Characterization

The computed pI value for NBD, K71L and T204V (pI < 7) indicated their acidic character. The extinction coefficient was calculated as 20,525–20,625 M−1·cm−1 based on the molar extinction coefficient of TYR, TRP and CYS residues. This measure indicates how much light is absorbed by a protein at a particular wave length. On the basis of instability index, Expasy’s ProtParam [19] classified NBD, K71L and T204V proteins as stable (instability index < 40). Instability index relies upon the occurrence of certain dipeptides along the length of the enzyme. The aliphatic index is defined as the relative volume of a protein that is occupied by an aliphatic side chain. An increase in the aliphatic index increases the thermo stability of globular proteins. The very high aliphatic index of all NBD and mutant proteins infers that these proteins may be stable for a wide range of temperatures. The very low grand average of hydropathicity (GRAVY) index (a negative value GRAVY) of NBD and all the mutant proteins infers that these proteins could result in a better interaction with water (hydrophilic in nature) (Table 1). The secondary structure indicates whether a given amino acid lies in a helix, strand or coil. The results from the SOPMA server [20] revealed that alpha helix dominated among secondary structure elements followed by random coils, extended strand and beta turns (Table 2).

2.3. Model Simulation and Evaluation

In this paper, we performed three 50 ns (50,000 ps) MD simulations of NBD, K71L and T204V to explore and compare the protein internal dynamics. To analyze the global behavior of the studied systems, the root mean square deviations (RMSDs) of the protein backbone with respect to the initial conformation were plotted versus simulation time (Figure 1A). RMSD played an important role in protein stability. The RMSD of the NBD model increased slightly and stabilized at 16,000 ps. K71L RMSD value increased gradually and stabilized at 5000 ps. It decreased slightly and then increased until it stabilized at 12,000 ps. At 40,000 ps, the NBD increased until 45,000 ps and dropped. In T204V’s simulation, the RMSD increased and stabilized at 5000 ps. T204V then increased slightly until it attained a constant level at 21,000 ps; and decreased at 31,000 ps until it stabilized at 35,000 ps. From the RMSD analysis, the T204V structure seemed more stable compared to the K71L and NBD structures due to lower fluctuation of RMSD values from 25,000 to 50,000 ps compared to those of the NBD and K71L. The root mean square fluctuation (RMSF) of the Cα atom of NBD, K71L and T204V as a function of residue number was plotted to evaluate the average fluctuation of each residue during the simulation (Figure 1B). All the residues in mutant models fluctuated around the NBD value which was approximately 0.2 nm. The RMSF of the Cα value of T204V exhibited a higher fluctuation than NBD and K71L at residue 56 (0.88 nm). In addition, NBD and mutant models stabilized at a gyration distance of about 2.11 nm at 20,000 ps (Figure 1C).
After 50,000 ps MD simulation, the geometry of three dimensional protein models was carried out with Ramachandran’s plot calculations using PROCHECK [21]. In the current study, the stereo-chemical evaluation of backbone psi and Phi dihedral angles of the NBD revealed that 81.7%, 15.4%, 2.1% and 0.9% of residues were falling within the most favored regions, additionally allowed regions, generously allowed regions and disallowed regions (ASN33, ASP44 and ASP97) respectively. In general, a score close to 100% implies good stereo-chemical quality of the model [22]. Therefore, these PROCHECK results suggest that the predicted model was of good quality (Figure 2).
The results also show that residues of the NBD model in the most favorable region were more than 80% except K71L and T204V mutants which scored slightly less than 80% for the most favored region (Table 3). However, the stereo chemical quality of the predicted models were found to be satisfactory and a low percentage of residues having phi/psi angles in the outlier region. The analysis explored that no bad contacts and no bad scores for main-chain or side-chain parameters. In spite of that, the overall G-factor values of NBD and mutants were slightly out of range because the values were lower than −0.5 but higher than −1.0. The acceptable values of the G-factor in PROCHECK are between 0 and −0.5, with the best quality models displaying values close to zero [23]. The quality of the protein structures were checked using ProQ [24]. The results show that the predicted LG score (>4: extremely good model) and predicted MaxSub score (>0.5 good model) for all protein models were in an acceptable range of a good model (Table 3).
The protein structures were also validated by other structure verification servers such as Verify 3D and ERRAT to check the quality of the models. ERRAT works by analyzing the statistics of non-bonded interactions between different atom types, with higher scores indicating higher quality [25]. The ERRAT score for T204V was the highest (94.536%). Nevertheless, this score was close with score values of the NBD (91.530%) and K71L (89.607%) proteins (Figure 3). None of the residues were above the 99% cut off of error-value. However, the generally accepted range is >50 for a high quality model [21]. Thus, this analysis revealed that the backbone conformation and non-bonded interactions of the NBD and mutant models fit well within the range of a high quality model.
In the Verify 3D analysis, it was found that none of the amino acids had a negative score (Figure 4). Therefore, the predicted models were compatible with its amino acid sequence. It should be noted that compatibility scores above zero correspond to an acceptable side chain environment [26].
PROSA was used to check three dimensional models of proteins under this study for potential errors [27]. The program displays two characteristics of the input structure: its Z-score and a plot of its residue energies. The Z-score indicates overall model quality and measures the deviation of the total energy of the structure with respect to an energy distribution derived from random conformations. Analysis of the K71L and T204V mutants with PROSA showed a Z-Score of −10.38 and −10.63 respectively, indicating no significant deviation from typical native structures of a similar size as the target protein’s Z-Score which was −10.64 (Figure 5). The quality of the protein folds of NBD, K71L and T204V were also evaluated in terms of energy function of amino acid residues. In general, folding energy of the protein showed minimum value as this accounts for the stability and nativity of the molecules. The energy profiles of mutant models in comparison to that of the X-ray structure of the NBD is presented in Figure 5. The energy profile of the K71L and T204V mutant models were consistent with a reliable conformation based on its similarity to that of the NBD.

2.4. Active Site Identification

Size volume, protein volume of active site (Table 4) and the residues forming pocket were obtained using Q-SiteFinder (Figure 6).

2.5. Molecular Docking

The negative and low value of ΔGbind (−9.09 kcal/mol) indicated strong bonds between T204V and the E1A32 kDa motif, and demonstrated that the protein was in a favorable conformation. Sum of intermolecular energy and torsion energy was the binding energy. Furthermore, the total intermolecular energy of T204V (−12.97 kcal/mol) was found to be lower than the NBD and K71L (−11.93 and −10.64 kcal/mol), stating that the mutant model (T204V) had a better binding affinity than the NBD and K71L in this analysis (Table 5). Hydrogen bonds formed between the compound and the protein usually contribute to the stability of the protein-ligand complexes; a large number of hydrogen bonds form more stable complexes [28,29]. The results show that the NBD, K71L and T204V were stabilized by three, two and four hydrogen bonds with E1A32 kDa motif, respectively (Table 6 and Figure 7). The active residues of the T204V mutant (THR13, THR14 and ARG72) were also involved in the formation of hydrogen bonds, suggesting the protein (T204V) forms a more stable complex than the NBD and K71L (Table 6). Therefore, in the study of protein-ligand binding mechanism, it was revealed that the novel T204V mutant has stronger interaction energy with the E1A32 kDa motif than other protein models.

2.6. Model Simulation and Evaluation of Protein-Ligand Complex

In this study, three 50 ns (50,000 ps) molecular dynamics simulation were run with NBD, K71L and T204V-E1A32 kDa motif complexes. Figure 8A shows the RMSD value of the protein-ligand complex structures over the simulation time. The RMSD value of 0.15 nm for NBD, K71L and T204V-E1A32 kDa motif complex structures were less deviated until 50 ps from their starting structure. NBD-E1A32 kDa motif complex reached stabilization at 35,000 ps (0.17 nm) while the RMSD attained a stable value of 0.20 nm at 35,000 ps for the K71L-E1A32 kDa complex. The T204V-E1A32 kDa complex attained 0.15 nm of RMSD backbone at 35,000 ps during simulation. From these results, it can be concluded that T204V-E1A32 kDa complex deviated less compared to the NBD and K71L-E1A32 kDa complexes indicating that the T204V-E1A32 kDa complex was more stable than the other two complex structures. The RMSF values of carbon alpha for each amino acid residue were obtained from the trajectory data of NBD, K71L and T204V-E1A32 kDa motif complexes shown in Figure 8B. In RMSF analysis, all the residues in the protein model fluctuated between 0.05 and 0.20 nm throughout the simulation period. The NBD-E1A32 kDa motif complex exhibited a high fluctuation up to 0.37 nm at residue 86.
History independent hydrogen bond autocorrelation function was calculated between the protein and the ligand (Figure 8C). The independent autocorrelation function measured the probability of a hydrogen bond present (broken or reformed between a time interval allowed) at time (t), given that it was present at time zero. The hydrogen bond analysis revealed that the life-time of a hydrogen bond between protein and ligand is highest for T204V followed by K71L and NBD; the results imply that the association between the protein and the ligand T204V is stronger.
Salt bridges formed between the amino acid side chains at positive ions in NBD, K71L and T204V, and negative ions in the E1A32 kDa motif. The salt bridge is important for stabilizing the protein’s structure. The presence of salt-bridges was a proof of the close proximity in the structure [30]. Salt-bridges occurring between the NBD, K71L and T204V-E1A32 kDa motif were calculated (Figure 8D). All three complexes attained the stable distance of 2.70, 2.60 and 2.40 nm in the entire simulation period. This suggested that the salt bridge with the shortest distance stabilizes the protein the most.

3. Experimental Section

3.1. Target Sequence

The tertiary structure of the ATPase domain of human heat shock 70 kDa protein 1 was publicly available. The complete amino acid sequence of ATPase domain of human heat shock 70 kDa protein 1, which consists of 380 amino acids, was retrieved from the RCSB Protein Databank (PDB: 1 HJO). The 3D structure of 1 HJO protein viewed using PyMol software [31].

3.2. In-Silico Mutagenesis

The functional amino acid residues (K71 and T204) were mutated computationally. These residues played a crucial role in catalytic activity and stabilization of the protein structure. PyMol software was used to alter the amino acid residues.

3.3. Physiochemical Characterization

Protein structure analysis was performed using the Expasy’s ProtParam Proteomics server [19]. The secondary structures prediction was carried out using Self Optimized Prediction Method from Alignment (SOPMA) [20].

3.4. Model Simulation and Evaluation

The Gromacs package 4.6.3 and the GROMOS 53a6 force field were used to perform MD simulation [32]. In the MD simulations, the NBD and the mutants were simulated using the GROMACS MD simulation software to examine the structural stability at a temperature of 300 K (27 °C). The protein models were solvated in a box of explicit simple point charge (SPC) water molecules and simulated using periodic boundary conditions (PBC) and particle mesh Ewald (PME) summation to take into account the long range electrostatic interactions. One Na+ ion was added to neutralise the total charge of the system for NBD and T204V proteins, meanwhile two sodium ions were added for the K71L protein simulation box. 749, 1004 and 961 steps of steepest descent energy minimization were carried out for NBD, K71L and T204V proteins, respectively. After energy minimization, the system was equilibrated at a constant temperature and pressure for 50 ps. The equilibrated structures were then subjected to molecular dynamic simulations for 50 ns (50,000 ps), LINCS constraint algorithm and 2-fs time step were set to run the simulation. All of the resulting trajectories were analysed using GROMACS utilities. RMSD and RMSF relative to the initial structure were calculated. The stereochemical quality and accuracy of the predicted models were assessed by PROCHECK program [21]. Validation of generated models was further performed by ProQ [24], ERRAT [25], Verify 3D [26] and ProSA [27] programs.

3.5. Active Site Identification

The binding sites of the protein were identified using Q-SiteFinder [33,34].

3.6. Homology Modeling

The three dimensional model of E1A32 kDa of human adenovirus serotype 5 was not available in the protein database at this time. The complete amino acid sequence of E1A32 kDa was retrieved from UniProtKB (accession number: P03255). BLASTP against the RCSB Protein Databank was carried out to find a suitable template for homology modeling [34]. Crystal structure of (PDB ID: 2 KJE) was selected as a template based on maximum identity with high positives and lower gaps percentage. The percentage of query coverage, sequence identity, positive and gap between the template and target protein were 13%, 100%, 100% and 0% respectively. The three dimensional structure of E1A32 kDa was built using EasyModeller 2.1 software [35], the Graphical User Interface (GUI) of Modeller 9.10, [36] and the model was then viewed using PyMol software. The three dimensional model of the E1A32 kDa motif (PNLVP) was created using the built three-dimensional model of E1A32 kDa as a template. The same homology modeling and 50 ns (50,000 ps) MD simulation approaches were performed before docking with NBD, K71L and T204V.

3.7. Molecular Docking

To understand the molecular interactions between E1A32 kDa motif of human adenovirus serotype 5 motif (PNLVP) and the NBD, flexible small molecule rigid protein docking was performed using Autodock Version 4.2 [37]. In the protein, non-polar hydrogen atoms were merged and total Kollman and Gasteiger charge was added to the protein. It was made sure that there were no non-bonded atoms in the protein. Kollman and Gasteiger partial charges were also assigned to the ligand and all torsions were allowed to rotate during docking. NBD and ligand were converted from PDB format to PDBQT format. Residues of the active site were also specified. A grid box was used around the active site to cover the entire protein-binding site and to allow ligands to move freely; and affinity maps of NBD, K71L and T204V (74 × 88 × 108, 70 × 60 × 70, 60 × 70 × 95 containing total grid points of 727,575, 307,501, 411,445 respectively) were calculated using AutoGrid. One hundred Lamarckian Genetic Algorithm (LGA) runs with default parameter settings were performed. Docking was reclustered for 0.5, 1.0 and 2.0 tolerances. The largest docked conformations were clustered at RMS of 1.0 nm and played ranked according to the native Autodock scoring function. The best conformation with the lowest docked energy was chosen from the docking search. The interactions of complex NBD protein-ligand conformations including hydrogen bonds and bond lengths were analyzed. The same docking simulation approach was performed with the single point mutants of NBD (K71L and T204V).

3.8. Molecular Dynamics Simulation of Protein-Ligand Complex

The docked complexes of E1A32 kDa motif with NBD and mutants (K71L and T204V) were used as a starting point for MD simulation. The GROMACS package 4.6.3 adopting the GROMOS53a6 force field parameter was carried out to run MD simulation [32]. The protein structures were solvated in cubic box 0.9 nm, using periodic boundary conditions and the SPC (simple point charge) water model. The ligand topology file was generated using the PRODRG server to include heteroatom due to limitations of GROMACS to parameterize the heteroatom group in PDB file [38]. The total charge of the system was neutralized by adding one sodium ion around the molecule for NBD and T204V proteins whereas two sodium ions were added for K71L protein. 993, 945 and 1023 steps of steepest descent energy minimization were carried out for NBD, K71L and T204V proteins respectively. The system was then equilibrated at a constant temperature (303 K) and pressure for 50 ps. Finally, the equilibrated structures were subjected to molecular dynamic simulations for 50 ns (50000 ps) with a LINCS algorithm 2-fs time step. The non-bonded list was generated using an atom-based cut-off of 10 Å. The long range electrostatic interactions were handled by the particle-mesh Ewald algorithm. The trajectory snapshots were taken for structural analysis at every pico-second. RMSD, RMSF, H-bonds and salt bridge formed between the protein and ligand in the docked complex during the simulation were analyzed through Gromacs utilities g_rmsd, g_rmsf, g_hbond and g_salt respectively.

4. Conclusions

In this study, we developed 3D mutant models (K71L and T204V), using NBD as a template. The structures were evaluated using RMSD, RMSF analysis, PROCHECK, ProQ, ERRAT, Verify 3D and ProSA programs. All the evidence suggests that the geometric quality of the backbone conformation, energy profile, residue interaction and contact of the predicted 3D structures of K71L and T204V mutants were well within the limits of reliable structures. The protein-protein docking study was elucidated for the purpose of finding protein interaction between Hsp70 and E1A32 kDa of human adenovirus serotype 5 motif; which might help in exploiting a potentially synergistic interaction to enhance anti-tumoral efficacy and decrease toxicity in cancer treatment. The interaction energy of docking between Hsp70 and E1A32 kDa motif was calculated and analyzed using the Autodock 4.2 programme. The results show that the novel T204V mutant was found to have the lowest binding energy (−9.09 kcal/mol) and intermolecular energy (−12.97 kcal/mol) among other protein models. Moreover, the protein-ligand complex structures were validated by RMSD, RMSF, hydrogen bonds and salt bridge analysis. It proved that T204V had strong bonds with E1A32 kDa motif and the complex structure was stable. Therefore, further biochemical and in vivo investigations of in-silico interpretations of this protein structure will be studied for the development of new therapy for the efficient cancer treatment.

Acknowledgments

This work was supported by research funding from Faculty of Bioscience and Medical Engineering, Universiti Teknologi Malaysia.

Conflicts of Interest

The authors declare no conflict of interest.
  • Author ContributionsAE designed the research, performed the experiments, analysed data and wrote the paper. MAN and SH contributed critical comments and suggestions. All authors read and approved the final manuscript.

References

  1. Hightower, L.E. Heat shock, stress proteins, chaperones, and proteotoxicity. Cell 1991, 66, 191–197. [Google Scholar]
  2. Gething, M.J.; Sambrook, J. Protein folding in the cell. Nature 1992, 355, 33–45. [Google Scholar]
  3. Netzer, W.J.; Hartl, F.U. Protein folding in the cytosol: Chaperon independent and -independent mechanisms. Trends Biochem. Sci 1998, 23, 68–73. [Google Scholar]
  4. Freeman, B.C.; Yamamoto, K.R. Disassembly of transcriptional regulatory complexes by molecular chaperones. Science 2002, 296, 2232–2235. [Google Scholar]
  5. Thorne, S.H.; Brooks, G.; Lee, Y.L.; Au, T.; Eng, L.F.; Reid, T. Effects of febrile temperature on adenoviral infection and replication: Implications for viral therapy of cancer. J. Virol 2005, 79, 581–591. [Google Scholar]
  6. Kao, H.T.; Capasso, O.; Heintz, N.; Nevins, J.R. Cell cycle control of the human HSP70 gene: Implications for the role of a cellular E1A-like function. Mol. Cell. Biol 1985, 5, 628–633. [Google Scholar]
  7. Wickner, S.; Skowyra, D.; Hoskins, J.; Mckenney, K. DnaJ, DnaK, and GrpE heat shock proteins are required in oriP1 DNA replication solely at the RepA monomerization step. Proc. Natl. Acad. Sci. USA 1992, 89, 10345–10349. [Google Scholar]
  8. Glotzer, J.B.; Saltik, M.; Chiocca, S.; Michou, A.I.; Moseley, P.; cotton, M. Activation of heat-shock response by an adenovirus is essential for virus replication. Nature 2000, 407, 207–211. [Google Scholar]
  9. Freeman, B.C.; Myers, M.P.; Schumacher, R.; Morimoto, R.I. Identification of a regulatory motif in Hsp70 that affects ATPase activity, substrate binding and interaction with HDJ-1. EMBO J 1995, 14, 2281–2292. [Google Scholar]
  10. Freeman, B.C.; Morimoto, R.I. The human cytosolic molecular chaperones Hsp90, Hsp70 (Hsc70), and HDJ-1 have distinct roles in recognition of a non-native protein and protein refolding. EMBO J 1996, 15, 2969–2979. [Google Scholar]
  11. Biro, J.C.; Benyo, B.; Sansom, C.; Szlavecz, A.; Fordos, G.; Micsik, T.; Benyo, Z. A common periodic table of codons and amino acids. Biochem. Biophys. Res. Commun 2003, 306, 408–415. [Google Scholar]
  12. O’Brien, M.C.; Flaherty, K.M.; McKay, D.B. Lysine 71 of the chaperone protein Hsc70 is essential for ATP hydrolysis. J. Biol. Chem 1996, 271, 15874–15878. [Google Scholar]
  13. Flaherty, K.M.; Wilbanks, S.M.; DeLuca-Flaherty, C.; McKay, D.B. Structural basis of the 70 kilodalton heat shock cognate protein ATP hydrolytic activity. J. Biol. Chem 1994, 269, 12899–12907. [Google Scholar]
  14. Wilbanks, S.M.; McKay, D.B. How potassium affects the activity of the molecular chaperone Hsc70. J. Biol. Chem 1995, 270, 2251–2257. [Google Scholar]
  15. McCarty, J.S.; Walker, G.C. DnaK as a thermometer: Threonine-199 is site of autophosphorylation and is critical for ATPase activity. Proc. Natl. Acad. Sci. USA 1991, 88, 9513–9517. [Google Scholar]
  16. Buchberger, A.; Schroder, H.; Buttner, M.; Valencia, A.; Bukau, B. A conserved loop in the ATPase domain of the DnaK chaperone is essential for stable binding of GrpE. Nat. Struct. Biol 1994, 1, 95–101. [Google Scholar]
  17. Wawrzynow, A.; Banecki, B.; Wall, D.; Liberek, K.; Georgopoulos, C. ATP hydrolysis is required for the DnaJ-dependent activation of DnaK chaperone for binding to both native and denatured protein substrates. J. Biol. Chem 1995, 270, 19307–19311. [Google Scholar]
  18. Campbell, K.S.; Mullane, K.P.; Aksoy, I.A.; Stubdal, H.; Zalvide, J.; Pipas, J.M.; Silver, P.A.; Roberts, T.M.; Schaffhausen, B.S.; DeCaprio, J.A. DnaJ/hsp40 chaperone domain of SV40 large T antigen promotes efficient viral DNA replication. Genes Dev 1997, 11, 1098–1110. [Google Scholar]
  19. Gasteiger, E.; Hoogland, C.; Gattiker, A.; Duvaud, S.; Wilkins, M.R.; Appel, R.D.; Bairoch, A. Protein identification and analysis tools on the ExPASy server. In The Proteomics Protocols Handbook; Walker, J.M., Ed.; Humana Press: Totowa, NJ, USA, 2005. [Google Scholar]
  20. Geourjon, C.; Deleage, G. SOPMA: Significant improvements in protein secondary structure prediction by consensus prediction from multiple alignments. Comput. Appl. Biosci 1995, 11, 681–684. [Google Scholar]
  21. Laskowski, R.A.; MacArthur, M.W.; Moss, D.S.; Thornton, J.M. PROCHECK: A program to check the stereochemical quality of protein structures. J. Appl. Cryst 1993, 26, 283–291. [Google Scholar]
  22. Reddy, Ch.S.; Vijayasarathy, K.; Srinivas, E.; Sastry, G.M.; Sastry, G.N. Homology modeling of membrane proteins: A critical assessment. Comput. Biol. Chem 2006, 30, 120–126. [Google Scholar]
  23. Gundampati, R.K.; Chikati, R.; Kumari, M.; Sharma, A.; Pratyush, D.D.; Jagannadham, M.V.; Kumar, C.S.; Das, M.D. Protein-protein docking on molecular models of Aspergillus niger RNase and human actin: Novel target for anticancer therapeutics. J. Mol. Model 2012, 18, 653–662. [Google Scholar]
  24. Wallner, B.; Elofsson, A. Can correct protein models be identified? Protein Sci 2003, 12, 1073–1086. [Google Scholar]
  25. Colovos, C.; Yeates, T.O. Verification of protein structures: Patterns of non-bonded atomic interactions. Protein Sci 1993, 2, 1511–1519. [Google Scholar]
  26. Eisenberg, D.; Luthy, R.; Bowie, J.U. VERIFY3D: Assessment of protein models with three-dimensional profiles. Methods Enzymol 1997, 277, 396–404. [Google Scholar]
  27. Wiederstein, M.; Sippl, M. ProSA-web: Interactive web service for the recognition of errors in three-dimensional structures of proteins. Nucleic Acids Res 2007, 35, W407–W410. [Google Scholar]
  28. Steiner, T.; Koellner, G. Hydrogen bonds with p–acceptors in proteins: Frequencies and role in stabilizing local 3–D structures. J. Mol. Biol 2001, 305, 535–557. [Google Scholar]
  29. Priya, George; Doss, C.; Nagasundaram, N. Investigating the structural impacts of I64T and P311S mutations in APE1-DNA complex: A Molecular Dynamics Approach. PLoS One 2012, 7, 1–11. [Google Scholar]
  30. Weiss, M.S.; Brandl, M.; Sühnel, J.; Pal, D.; Hilgenfeld, R. More hydrogen bonds for the (structural) biologist. Trends Biochem. Sci 2001, 26, 521–523. [Google Scholar]
  31. Delano, W.L. The PyMOL Molecular Graphics System, 2002. 10 September 2013. Available online: http://www.pymol.org.
  32. Van Der Spoel, D.; Lindahl, E.; Hess, B.; Groenhof, G.; Mark, A.E.; Berendsen, H.J.C. GROMACS: Fast, flexible, and free. J. Comput. Chem 2005, 26, 1701–1718. [Google Scholar]
  33. Laurie, A.T.; Jackson, R.M. Q-SiteFinder: An energy-based method for the prediction of protein-ligand binding sites. Bioinformatics 2005, 21, 1908–1916. [Google Scholar]
  34. Burgoyne, N.J.; Jackson, R.M. Predicting protein interaction sites: Binding hot-spots in protein-protein and protein-ligand interfaces. Bioinformatics 2006, 22, 1335–1342. [Google Scholar]
  35. Kuntal, B.K.; Aparoy, P.; Reddanna, P. EasyModeller: A graphical interface to MODELLER. BMC Res. Notes 2010, 3, 226. [Google Scholar]
  36. Fiser, A.; Sali, A. Modeller: Generation and refinement of homology-based protein structure models. Methods Enzymol 2003, 374, 461–491. [Google Scholar]
  37. Sanner, M.F. Python: A programming language for software integration and development. J. Mol. Graph. Model 1999, 17, 57–61. [Google Scholar]
  38. Schuttelkopf, A.W.; van Aalten, D.M. PRODRG—A tool for high-throughput crystallography of protein-ligand complexes. Acta Crystallogr. D Biol. Crystallogr 2004, 60, 1355–1363. [Google Scholar]
Ijms 15 06797f1 1024
Figure 1. Dynamic changes of the NBD upon mutations. (A) Root mean square deviations (RMSD) of NBD (black), K71L (yellow) and T204V (brown); (B) Backbone atomic fluctuations (RMSF) of NBD (black), K71L (red) and T204V (green); (C) Radius gyration of NBD (black), K71L (red) and T204V (green).
Figure 1. Dynamic changes of the NBD upon mutations. (A) Root mean square deviations (RMSD) of NBD (black), K71L (yellow) and T204V (brown); (B) Backbone atomic fluctuations (RMSF) of NBD (black), K71L (red) and T204V (green); (C) Radius gyration of NBD (black), K71L (red) and T204V (green).
Ijms 15 06797f1
Ijms 15 06797f2 1024
Figure 2. Ramachandran plots generated via PROCHECK for (A) NBD protein; (B) K71L and (C) T204V mutants. PROCHECK shows that the residues in most favored (red), additionally allowed (yellow), generously allowed (pale yellow) and disallowed regions (white color).
Figure 2. Ramachandran plots generated via PROCHECK for (A) NBD protein; (B) K71L and (C) T204V mutants. PROCHECK shows that the residues in most favored (red), additionally allowed (yellow), generously allowed (pale yellow) and disallowed regions (white color).
Ijms 15 06797f2
Ijms 15 06797f3a 1024Ijms 15 06797f3b 1024
Figure 3. ERRAT plots for (A) NBD; (B) K71L and (C) T204V mutants. Black bars identify the misfolded region located distantly from the active site, gray bars demonstrate the error region between 95% and 99%, and white bars indicate the region with a lower error rate for protein folding. * On the error axis, two lines are drawn to indicate the confidence with which it is possible to reject regions that exceed that error value. ** Expressed as the percentage of the protein for which the calculated error value falls below the 95% rejection limit. Good high resolution structures generally produce values around 95% or higher. For lower resolutions (2.5 to 3 A) the average overall quality factor is around 91%.
Figure 3. ERRAT plots for (A) NBD; (B) K71L and (C) T204V mutants. Black bars identify the misfolded region located distantly from the active site, gray bars demonstrate the error region between 95% and 99%, and white bars indicate the region with a lower error rate for protein folding. * On the error axis, two lines are drawn to indicate the confidence with which it is possible to reject regions that exceed that error value. ** Expressed as the percentage of the protein for which the calculated error value falls below the 95% rejection limit. Good high resolution structures generally produce values around 95% or higher. For lower resolutions (2.5 to 3 A) the average overall quality factor is around 91%.
Ijms 15 06797f3aIjms 15 06797f3b
Ijms 15 06797f4 1024
Figure 4. Verify 3D plots for (A) NBD protein; (B) K71L and (C) T204V mutants. Each residue was assigned a structural class based on its location and environment (alpha, beta, loop, polar, nonpolar, etc.). A collection of good structures was used as a reference to obtain a score for each of the 20 amino acids in this structural class. The scores of a sliding 21-residue window (from −10 to +10) were added and plotted for individual residues.
Figure 4. Verify 3D plots for (A) NBD protein; (B) K71L and (C) T204V mutants. Each residue was assigned a structural class based on its location and environment (alpha, beta, loop, polar, nonpolar, etc.). A collection of good structures was used as a reference to obtain a score for each of the 20 amino acids in this structural class. The scores of a sliding 21-residue window (from −10 to +10) were added and plotted for individual residues.
Ijms 15 06797f4
Ijms 15 06797f5 1024
Figure 5. Protein quality scores for (A) NBD protein; (B) K71L and (C) T204V mutants generated through ProSA web server. The results generated display the Z-scores which indicate the overall model quality and energy plots which indicate the local model quality. PROSA-web Z-scores of all protein chains in PDB are determined by X-ray crystallography (light blue) and NMR spectroscopy (dark blue) with respect to their length. The Z-score of protein models were present in the range represented by the large black dot.
Figure 5. Protein quality scores for (A) NBD protein; (B) K71L and (C) T204V mutants generated through ProSA web server. The results generated display the Z-scores which indicate the overall model quality and energy plots which indicate the local model quality. PROSA-web Z-scores of all protein chains in PDB are determined by X-ray crystallography (light blue) and NMR spectroscopy (dark blue) with respect to their length. The Z-score of protein models were present in the range represented by the large black dot.
Ijms 15 06797f5
Ijms 15 06797f6 1024
Figure 6. Projection of the predicted active site for (A) NBD protein; (B) K71L and (C) T204V mutants; obtained using Q-SiteFinder web server (shown as red colour).
Figure 6. Projection of the predicted active site for (A) NBD protein; (B) K71L and (C) T204V mutants; obtained using Q-SiteFinder web server (shown as red colour).
Ijms 15 06797f6
Ijms 15 06797f7a 1024Ijms 15 06797f7b 1024
Figure 7. Docking of (A) NBD protein; (B) K71L and (C) T204V mutants with the E1A32 kDa motif. Hydrogen bonds are shown by a green line with its distance (Å).
Figure 7. Docking of (A) NBD protein; (B) K71L and (C) T204V mutants with the E1A32 kDa motif. Hydrogen bonds are shown by a green line with its distance (Å).
Ijms 15 06797f7aIjms 15 06797f7b
Ijms 15 06797f8 1024
Figure 8. Dynamic changes of the NBD, K71L and T204V-E1A32 kDa motif complex structures. (A) Root mean square deviations (RMSD) of NBD (black), K71L (yellow) and T204V (brown); (B) Backbone atomic fluctuations (RMSF) of NBD (black), K71L (red) and T204V (green); (C) Hydrogen bond autocorrelation of NBD (black), K71L (red) and T204V (green); (D) Salt bridge of NBD (black), K71L (red) and T204V (green).
Figure 8. Dynamic changes of the NBD, K71L and T204V-E1A32 kDa motif complex structures. (A) Root mean square deviations (RMSD) of NBD (black), K71L (yellow) and T204V (brown); (B) Backbone atomic fluctuations (RMSF) of NBD (black), K71L (red) and T204V (green); (C) Hydrogen bond autocorrelation of NBD (black), K71L (red) and T204V (green); (D) Salt bridge of NBD (black), K71L (red) and T204V (green).
Ijms 15 06797f8
Table
Table 1. The physiochemical character of the NBD protein, K71L and T204V mutants as predicted by Expasy’s Prot-Param program.
Table 1. The physiochemical character of the NBD protein, K71L and T204V mutants as predicted by Expasy’s Prot-Param program.
ProteinLengthM.wt (Daltons)pI−R+RExtinction coefficient (M−1·cm−1)Instability indexAliphatic indexGRAVY
NBD38041,827.46.69504920,52535.0988.32−0.274
K71L38041,812.46.38504820,52534.7489.34−0.253
T204V38041,825.56.69504920,52534.9989.08−0.261
Table
Table 2. Secondary structures of the NBD protein, K71L and T204V mutants.
Table 2. Secondary structures of the NBD protein, K71L and T204V mutants.
Secondary structureAlpha helix(Hh)Extended strand (Ee)Beta turn (Tt)Random coil (Cc)
NBD42.8919.747.6329.74
K71L44.2118.688.1628.95
T204V44.4718.956.8429.74
Table
Table 3. Validation of NBD, K71L and T204V mutants using PROCHECK and ProQ.
Table 3. Validation of NBD, K71L and T204V mutants using PROCHECK and ProQ.
StructureRamachandran plot statistics (%)Goodness factorProQ
Most favouredAdditionally allowedGenerously allowedDisallowedDihedral anglesCovalent forcesOverall averageLG ScoreMax–sub
NBD81.715.42.10.9−0.61 *−0.95 *−0.66 *5.7070.451
K71L79.918.60.90.6−0.59 *−1.02 **−0.67 *5.4970.425
T204V78.719.21.20.9−0.61 *−0.91 *−0.64 *5.8620.424
*Lower than −0.5-unusual;
**Lower than −1.0-highly unusual.
Table
Table 4. Predicted active sites of the NBD protein, K71L and T204V mutants.
Table 4. Predicted active sites of the NBD protein, K71L and T204V mutants.
ProteinSite volume (cubic Å)Protein volume (cubic Å)Residues that forming pocket
NBD43434357ASP10,LEU11,GLY12,THR13,THR14,TYR15,PHE68,ASP69,LYS71,ARG72,TRP90,THR145,VAL146,PRO147,ALA148,GLU175,PRO176,ILE197,PHE198,ASP199,GLY201,GLY202,GLY203,THR204,ASP206,VAL207,SER208,THR222,VAL337,VAL369
K71L49634418ASP10,LEU11,GLY12,THR13,PRO14,TYR15,CYS17,ARG72,VAL146,PRO147,ALA148,TYR149,GLU175,PRO176,ALA179,ILE197,PHE198,ASP199,LEU200,GLY201,GLY202,GLY203,THR204,ASP206,VAL207,SER208,THR222,ALA223,GLY224,LYS271,ARG272,VAL337,GLY338,GLY339,GLY340,ALA368,VAL369,ALA370
T204V53234224ASP10,LEU11,GLY12,THR13,THR14,TYR15,SER16,CYS17,LYS71,ARG72,VAL82,THR145,PRO147,ALA148,TYR149,PHE150,GLU175,ASP199,LEU200,GLY201,GLY202,GLY203,VAL204,PHE205,ASP206,ARG272,VAL337,GLY338,GLY339,PRO365,ASP366,GLU367,ALA368,VAL369,ALA370
Table
Table 5. Docking results of NBD protein, K71L and T204V mutants with the E1A32 kDa motif.
Table 5. Docking results of NBD protein, K71L and T204V mutants with the E1A32 kDa motif.
ProteinNBDK71LT204V
Binding energy (kcal/mol)−8.05−6.76−9.09
kI (μM)1.2611.040.22
Intermolecular Energy (kcal/mol)−11.93−10.64−12.97
Internal energy (kcal/mol)−2.45−2.55−1.49
Torsion energy (kcal/mol)3.883.883.88
Unbounded Extended energy (kcal/mol)−2.45−2.55−1.49
Cluster RMS0.000.000.00
Reference RMS86.0972.9376.96
Table
Table 6. Hydrogen bonds interaction studies of the NBD protein, K71L and T204V mutants with the E1A32 kDa motif.
Table 6. Hydrogen bonds interaction studies of the NBD protein, K71L and T204V mutants with the E1A32 kDa motif.
ProteinDonor atomAcceptor atomDistance (Å)
NBDARG72:HH21ASN2:OD11.913
TRP90:HE1ASN2:O2.006
ARG264:HH11PRO5:O,OXT2.000
K71LGLU231:HNPRO1:O1.867
ARG261:HEPRO5:O1.899
T204VTHR14:HNPRO5:O2.065
ARG72:HEPRO:O2.199
ARG72:HH21ASN2:O1.594
THR13:HNVAL4:O1.922

Share and Cite

MDPI and ACS Style

Elengoe, A.; Naser, M.A.; Hamdan, S. Modeling and Docking Studies on Novel Mutants (K71L and T204V) of the ATPase Domain of Human Heat Shock 70 kDa Protein 1. Int. J. Mol. Sci. 2014, 15, 6797-6814. https://doi.org/10.3390/ijms15046797

AMA Style

Elengoe A, Naser MA, Hamdan S. Modeling and Docking Studies on Novel Mutants (K71L and T204V) of the ATPase Domain of Human Heat Shock 70 kDa Protein 1. International Journal of Molecular Sciences. 2014; 15(4):6797-6814. https://doi.org/10.3390/ijms15046797

Chicago/Turabian Style

Elengoe, Asita, Mohammed Abu Naser, and Salehhuddin Hamdan. 2014. "Modeling and Docking Studies on Novel Mutants (K71L and T204V) of the ATPase Domain of Human Heat Shock 70 kDa Protein 1" International Journal of Molecular Sciences 15, no. 4: 6797-6814. https://doi.org/10.3390/ijms15046797

APA Style

Elengoe, A., Naser, M. A., & Hamdan, S. (2014). Modeling and Docking Studies on Novel Mutants (K71L and T204V) of the ATPase Domain of Human Heat Shock 70 kDa Protein 1. International Journal of Molecular Sciences, 15(4), 6797-6814. https://doi.org/10.3390/ijms15046797

Article Metrics

Back to TopTop