Search

Peroxiredoxins (Prx) have received considerable attention during recent years. This study demonstrates that two typical Pseudomonas-derived 2-Cys Prx proteins, PpPrx and PaPrx can alternatively function as a peroxidase and chaperone. The amino acid sequences of these two Prx proteins exhibit 93% homology, but PpPrx possesses an additional cysteine residue, Cys112, instead of the alanine found in PaPrx. PpPrx predominates with a high molecular weight (HMW) complex and chaperone activity, whereas PaPrx has mainly low molecular weight (LMW) structures and peroxidase activity. Mass spectrometry and structural analyses showed the involvement of Cys112 in the formation of an inter-disulfide bond, the instability of LMW structures, the formation of HMW complexes, and increased hydrophobicity leading to functional switching of Prx proteins between peroxidase and chaperone.

Background and aims The 2-Cys peroxiredoxin (Prx) A protein of Arabidopsis thaliana performs the dual functions of a peroxidase and a molecular chaperone depending on its conformation and the metabolic conditions. However, the precise mechanism responsible for the functional switching of 2-Cys Prx A is poorly known. This study examines various serine-to-cysteine substitutions on α-helix regions of 2-Cys Prx A in Arabidopsis mutants and the effects they have on the dual function of the protein. Methods Various mutants of 2-Cys Prx A were generated by replacing serine (Ser) with cysteine (Cys) at different locations by site-directed mutagenesis. The mutants were then over-expressed in Escherichia coli. The purified protein was further analysed by size exclusion chromatography, polyacrylamide gel electrophoresis, circular dichroism spectroscopy and transmission electron microscopy (TEM) and image analysis. Peroxidase activity, molecular chaperone activity and hydrophobicity of the proteins were also determined. Molecular modelling analysis was performed in order to demonstrate the relationship between mutation positions and switching of 2-Cys Prx A activity. Key results Replacement of Ser(150) with Cys(150) led to a marked increase in holdase chaperone and peroxidase activities of 2-Cys Prx A, which was associated with a change in the structure of an important domain of the protein. Molecular modelling demonstrated the relationship between mutation positions and the switching of 2-Cys Prx A activity. Examination of the α2 helix, dimer-dimer interface and C-term loop indicated that the peroxidase function is associated with a fully folded α2 helix and easy formation of a stable reduced decamer, while a more flexible C-term loop makes the chaperone function less likely. Conclusions Substitution of Cys for Ser at amino acid location 150 of the α-helix of 2-Cys Prx A regulates/enhances the dual enzymatic functions of the 2-Cys Prx A protein. If confirmed in planta, this leads to the potential for it to be used to maximize the functional utility of 2-Cys Prx A protein for improved metabolic functions and stress resistance in plants.

Recent evidence suggests that aldo-keto reductase family 1 B10 (AKR1B10) may be a potential diagnostic or prognostic marker of human tumors, and that AKR1B10 inhibitors offer a promising choice for treatment of many types of human cancers. The aim of this study was to identify novel chemical scaffolds of AKR1B10 inhibitors using in silico approaches. The 3D QSAR pharmacophore models were generated using HypoGen. A validated pharmacophore model was selected for virtual screening of 4 chemical databases. The best mapped compounds were assessed for their drug-like properties. The binding orientations of the resulting compounds were predicted by molecular docking. Density functional theory calculations were carried out using B3LYP. The stability of the protein-ligand complexes and the final binding modes of the hit compounds were analyzed using 10 ns molecular dynamics (MD) simulations. The best pharmacophore model (Hypo 1) showed the highest correlation coefficient (0.979), lowest total cost (102.89) and least RMSD value (0.59). Hypo 1 consisted of one hydrogen-bond acceptor, one hydrogen-bond donor, one ring aromatic and one hydrophobic feature. This model was validated by Fischer's randomization and 40 test set compounds. Virtual screening of chemical databases and the docking studies resulted in 30 representative compounds. Frontier orbital analysis confirmed that only 3 compounds had sufficiently low energy band gaps. MD simulations revealed the binding modes of the 3 hit compounds: all of them showed a large number of hydrogen bonds and hydrophobic interactions with the active site and specificity pocket residues of AKR1B10. Three compounds with new structural scaffolds have been identified, which have stronger binding affinities for AKR1B10 than known inhibitors.

In contrast to a number of studies on the humanization of non-human antibodies, the reshaping of a non-human antibody into a chicken antibody has never been attempted. Therefore, nothing is known about the animal species-dependent compatibility of the framework regions (FRs) that sustain the appropriate conformation of the complementarity-determining regions (CDRs). In this study, we attempted the reshaping of the variable domains of the mouse catalytic anti-nucleic acid antibody 3D8 (m3D8) into the FRs of a chicken antibody (“chickenization”) by CDR grafting, which is a common method for the humanization of antibodies. CDRs of the acceptor chicken antibody that showed a high homology to the FRs of m3D8 were replaced with those of m3D8, resulting in the chickenized antibody (ck3D8). ck3D8 retained the biochemical properties (DNA binding, DNA hydrolysis, and cellular internalizing activities) and three-dimensional structure of m3D8 and showed reduced immunogenicity in chickens. Our study demonstrates that CDR grafting can be applied to the chickenization of a mouse antibody, probably due to the interspecies compatibility of the FRs.

The modular auxiliary activity (AA) family of proteins is believed to cause amorphogenesis in addition to oxidative cleavage of crystalline cellulose although the supporting evidence is limited. HcAA10-2 is a modular AA10 family protein (58 kDa) composed of a AA10 module and a family two carbohydrate binding module (CBM2), joined by a long stretch of 222 amino acids of unknown function. The protein was expressed in Escherichia coli and purified to homogeneity. Scanning electron microscopy and X-ray diffraction analysis of Avicel treated with HcAA10-2 provided evidence for the disruption of the cellulose microfibrils ("amorphogenesis") and reduction of the crystallinity index, resulting in a twofold increase of cellulase adsorption on the polysaccharide surface. HcAA10-2 exhibited weak endoglucanase-like activity toward soluble cellulose and cello-oligosaccharides with an optimum at pH 6.5 and 45 °C. HcAA10-2 catalyzed oxidative cleavage of crystalline cellulose released native and oxidized cello-oligosaccharides in the presence of copper and an electron donor such as ascorbic acid. Multiple sequence alignment indicated that His1, His109, and Phe197 in the AA10 module formed the conserved copper-binding site. The reducing sugar released from Avicel by the endoglucanase Cel5 and Celluclast accompanying HcAA10-2 was increased by four- and sixfold, respectively. Moreover, HcAA10-2 and Celluclast acted synergistically on pretreated wheat straw biomass resulting in a threefold increase in reducing sugar than Celluclast alone. Taken together, these results suggest that HcAA10-2 is a novel multifunctional modular AA10 protein possessing amorphogenesis, weak endoglucanase, and oxidative cleavage activities useful for efficient degradation of crystalline cellulose.

Grb2 is an adapter protein involved in the signal transduction and cell communication. The Grb2 is responsible for initiation of kinase signaling by Ras activation which leads to the modification in transcription. Ligand based pharmacophore approach was applied to built the suitable pharmacophore model for Grb2. The best pharmacophore model was selected based on the statistical values and then validated by Fischer's randomization method and test set. Hypo1 was selected as a best pharmacophore model based on its statistical values like high cost difference (182.22), lowest RMSD (1.273), and total cost (80.68). It contains four chemical features, one hydrogen bond acceptor (HBA), two hydrophobic (HY), and one ring aromatic (RA). Fischer's randomization results also shows that Hypo1 have a 95% significant level. The correlation coefficient of test set was 0.97 which was close to the training set value (0.94). Thus Hypo1 was used for virtual screening to find the potent inhibitors from various chemical databases. The screened compounds were filtered by Lipinski's rule of five, ADMET and subjected to molecular docking studies. Totally, 11 compounds were selected as a best potent leads from docking studies based on the consensus scoring function and critical interactions with the amino acids in Grb2 active site.

Hsp33, a prokaryotic molecular chaperone, exerts holdase activity in response to oxidative stress. In this study, the stepwise conformational change of Hsp33 upon oxidation was monitored by NMR. In order to overcome its high molecular weight (33 kDa as a monomer and 66 kDa as a dimer), spectra were simplified using a selectively ( 15 N)His-labeled protein. All of the eight histidines were observed in the TROSY spectrum of the reduced Hsp33. Among them, three peaks showed dramatic resonance shifts dependent on the stepwise oxidation, indicating a remarkable conformational change. The results suggest that unfolding of the linker domain is associated with dimerization, but not entire region of the linker domain is unfolded.

The acetylcholinesterase (AChE) is important to terminate acetylcholine-mediated neurotransmission at cholinergic synapses. The pivotal role of AChE in apoptosome formation through the interactions with cytochrome c (Cyt c) was demonstrated in recent study. In order to investigate the proper binding conformation between the human AChE (hAChE) and human Cyt c (hCyt c), macro-molecular docking simulation was performed using DOT 2.0 program. The hCyt c was bound to peripheral anionic site (PAS) on hAChE and binding mode of the docked conformation was very similar to the reported crystal structure of the AChE and fasciculin-II (Fas-II) complex. Two 10 ns molecular dynamics (MD) simulations were carried out to refine the binding mode of docked structure and to observe the differences of the binding conformations between the absent (Apo) and presence (Holo) of heme group. The key hydrogen bonding residues between hAChE and hCyt c proteins were found in Apo and Holo systems, as well as each Tyr341 and Trp286 residue of hAChE was participated in cation-pi (π) interactions with Lys79 of hCyt c in Apo and Holo systems, respectively. From the present study, although the final structures of the Apo and Holo systems have similar binding pattern, several differences were investigated in flexibilities, interface interactions, and interface accessible surface areas. Based on these results, we were able to predict the reasonable binding conformation which is indispensable for apoptosome formation.

The behavior of peptide or protein solutes in saline aqueous solution is a fundamental topic in physical chemistry. Addition of ions can strongly alter the thermodynamic and physical properties of peptide molecules in solution. In order to study the effects of added ionic salts on protein conformation and dynamics, we have used the molecular dynamics (MD) simulations to investigate the behavior of Staphylococcus aureus Hfq protein under two different ionic concentrations: 0.1 M NaCl and 1.0 M NaCl in presence and absence of RNA (a hepta-oligoribonucleotide AU5G). Hfq, a global regulator of gene expression is highly conserved and abundant RNA-binding protein. It is already reported that in vivo the increase of ionic strength results in a drastic reduction of Hfq affinity for Qβ RNA and reduces the tendency of aggregation of Escherichia coli host factor hexamers. Our results revealed the crucial role of 0.1 M NaCl Hfq system on the bases with strong hydrogen bonding interactions and by stabilizing the aromatic stacking of Tyr42 residue of the adjacent subunits/monomers with the adenine and uridine nucleobases. An increase in RNA pore diameter and weakened compactness of the Hfq-RNA complex was clearly observed in 1.0 M NaCl Hfq system with bound RNA. Aggregation of monomers in Hfq and the interaction of Hfq with RNA are greatly affected due to the presence of high ionic strength. Higher the ionic concentration, weaker is the aggregation and interaction. Our results were compatible with the experimental data and this is the first theoretical report for the experimental study done in 1980 by Uhlenbeck group for the present system.

Substituted imidazolidin-2-ones deduced as potential inhibitors of IleRS by docking simulations were synthesized from an aziridine-2-carboxaldehyde. Reductive amination of an aziridine-2-carboxaldehyde with dipeptides for the sub- stituents at N1 and followed by aziridine-ring expansion with triphosgene afforded 4-chloromethylimidazolidin-2-ones whose chloride were further manipulated towards phenylurea, pyrimidin-2-yl-urea or benzenesulfonamide at C4.

Human histamine H 1 receptor (HHR1) is a G protein-coupled receptor and a primary target for antiallergic therapy. Here, the ligand-based three-dimensional pharmacophore models were built from a set of known HHR1 inverse agonists using HypoGen module of CATALYST software. All ten generated pharmacophore models consist of five essential features: hydrogen bond acceptor, ring aromatic, positive ionizable and two hydrophobic functions. Best model had a correlation coefficient of 0.854 for training set compounds and it was validated with an external test set with a high correlation value of 0.925. Using this model Maybridge database containing 60,000 compounds was screened for potential leads. A rigorous screening for drug-like compounds unveiled RH01692 and SPB00834, two novel molecules for HHR1 with good CATALYST fit and estimated activity values. The new lead molecules were docked into the active site of constructed HHR 1 homology model based on recently crystallized squid rhodopsin as template. Both the hit compounds were found to have critical interactions with Glu177, Phe432 and other important amino acids. The interpretations of this study may effectively be deployed in designing of novel HHR1 inverse agonists.

Silent information regulator 2 (Sir2) or sirtuins are NAD(+)-dependent deacetylases, which hydrolyze the acetyl-lysine residues. In mammals, sirtuins are classified into seven different classes (SIRT1-7). SIRT1 was reported to be involved in age related disorders like obesity, metabolic syndrome, type II diabetes mellitus and Parkinson’s disease. Activation of SIRT1 is one of the promising approaches to treat these age related diseases. In this study, we have used HipHop module of CATALYST to identify a series of pharmacophore models to screen SIRT1 enhancing molecules. Three molecules from Sirtris Pharmaceuticals were selected as training set and 607 sirtuin activator molecules were used as test set. Five different hypotheses were developed and then validated using the training set and the test set. The results showed that the best pharmacophore model has four features, ring aromatic, positive ionization and two hydrogen-bond acceptors. The best hypothesis from our study, Hypo2, screened high number of active molecules from the test set. Thus, we suggest that this four feature pharmacophore model could be helpful to screen novel SIRT1 activator molecules. Hypo2-virtual screening against Maybridge database reveals seven molecules, which contains all the critical features. Moreover, two new scaffolds were identified from this study. These scaffolds may be a potent lead for the SIRT1 activation.

Adenosine kinase (AK) is a ubiquitous intracellular enzyme, which catalyzes the phosphorylation of adenosine (ADO) to adenosine monophosphate (AMP). AK inhibitors have therapeutic potential as analgesic and antiinflammatory agents. A chemical feature based pharmacophore model has been generated from known AK inhibitors (26 training set compounds) by HypoGen module implemented in CATALYST software. The top ranked hypothesis (Hypo1) contained four features of two hydrogen-bond acceptors (HBA) and two hydrophobic aromatics (Z). Hypo1 was validated by 124 test set molecules with a correlation coefficient of 0.905 between experimental and estimated activity. It was also validated by CatScramble method. Thus, the Hypo1 was exploited for searching new lead compounds over 238,819 chemical compounds in NCI database and then the selected compounds were screened based on restriction estimated activity and Lipinski's rules to evaluate their drug-like properties. Finally we could obtain 72 new lead candidates and the two best compound structures from them were posted.

The activation process of the redox-regulated chaperone heat shock protein 33 (Hsp33) is constituted by the oxidation-induced unfolding of the C-terminal zinc-binding domain and concomitant oligomerization of the N-terminal core domain. Herein, the semi-empirical solution structure of Escherichia coli Hsp33 in the reduced, inactive form was generated through conformational space annealing calculations, utilizing minimalistic NMR data and multiple homology restraints. The various conformations of oxidized Hsp33 and some mutant forms were also investigated in solution. Interestingly, a specific region concentrated around the interdomain linker stretch and its interacting counterparts, the N-terminal β-strand 1 and α-helix 1, hardly showed up as signals in the NMR measurements. The NMR spectra of an Hsp33 derivative with a six-residue deletion in the disordered N-terminus implied a plausible conformational exchange associated with the identified region, and the corresponding exchange rate appeared slower than that of the wild type. Subsequent mutations that destroyed the structure of the β1 or α1 elements resulted in the formation of a reduced but active monomer, without the unfolding of the zinc-binding domain. Collectively, structural insights into the inactive and active conformations, including wild-type and mutant proteins, suggest that the dynamic interactions of the N-terminal segments with their contacting counterpart, the interdomain linker stretch, in the reduced, inactive state are the structural determinants regulating the activation process of the post-translationally regulated chaperone, Hsp33.

AtTDX, a thioredoxin-like plant-specific protein present in Arabidospis is a thermo-stable and multi-functional enzyme. This enzyme is known to act as a thioredoxin and as a molecular chaperone depending upon its oligomeric status. The present study examines the effects of γ-irradiation on the structural and functional changes of AtTDX. Holdase chaperone activity of AtTDX was increased and reached a maximum at 10 kGy of γ-irradiation and declined subsequently in a dose-dependent manner, together with no effect on foldase chaperone activity. However, thioredoxin activity decreased gradually with increasing irradiation. Electrophoresis and size exclusion chromatography analysis showed that AtTDX had a tendency to form high molecular weight (HMW) complexes after γ-irradiation and γ-ray-induced HMW complexes were tightly associated with a holdase chaperone activity. The hydrophobicity of AtTDX increased with an increase in irradiation dose till 20 kGy and thereafter decreased further. Analysis of the secondary structures of AtTDX using far UV-circular dichroism spectra revealed that the irradiation remarkably increased the exposure of β-sheets and random coils with a dramatic decrease in α-helices and turn elements in a dose-dependent manner. The data of the present study suggest that γ-irradiation may be a useful tool for increasing holdase chaperone activity without adversely affecting foldase chaperone activity of thioredoxin-like proteins.

Multiple isoforms of Arabidopsis thaliana h-type thioredoxins (AtTrx-hs) have distinct structural and functional specificities. AtTrx-h3 acts as both a disulfide reductase and as a molecular chaperone. We prepared five representative AtTrx-hs and compared their protein structures and disulfide reductase and molecular chaperone activities. AtTrx-h2 with an N-terminal extension exhibited distinct functional properties with respect to other AtTrx-hs. AtTrx-h2 formed low-molecular-mass structures and exhibited only disulfide reductase activity, whereas the other AtTrx-h isoforms formed high-molecular-mass complexes and displayed both disulfide reductase and molecular chaperone activities. The domains that determine the unique structural and functional properties of each AtTrx-hs protein were determined by constructing a domain-swap between the N- and C-terminal regions of AtTrx-h2 and AtTrx-h3 (designated AtTrx-h-2N3C and AtTrx-h-3N2C respectively), an N-terminal deletion mutant of AtTrx-h2 [AtTrx-h2-N(∆19)] and site-directed mutagenesis of AtTrx-h3. AtTrx-h2-N(∆19) and AtTrx-h-3N2C exhibited similar properties to those of AtTrx-h2, but AtTrx-h-2N3C behaved more like AtTrx-h3, suggesting that the structural and functional specificities of AtTrx-hs are determined by their C-terminal regions. Hydrophobicity profiling and molecular modelling revealed that Ala100 and Ala106 in AtTrx-h3 play critical roles in its structural and functional regulation. When these two residues in AtTrx-h3 were replaced with lysine, AtTrx-h3 functioned like AtTrx-h2. The chaperone function of AtTrx-hs conferred enhanced heat-shock-resistance on a thermosensitive trx1/2-null yeast mutant.

Human chymase catalyzes the hydrolysis of peptide bonds. Three chymase inhibitors with very similar chemical structures but highly different inhibitory profiles towards the hydrolase function of chymase were selected with the aim of elucidating the origin of disparities in their biological activities. As a substrate (angiotensin-I) bound crystal structure is not available, molecular docking was performed to dock the substrate into the active site. Molecular dynamics simulations of chymase complexes with inhibitors and substrate were performed to calculate the binding orientation of inhibitors and substrate as well as to characterize conformational changes in the active site. The results elucidate details of the 3D chymase structure as well as the importance of K40 in hydrolase function. Binding mode analysis showed that substitution of a heavier Cl atom at the phenyl ring of most active inhibitor produced a great deal of variation in its orientation causing the phosphinate group to interact strongly with residue K40. Dynamics simulations revealed the conformational variation in region of V36-F41 upon substrate and inhibitor binding induced a shift in the location of K40 thus changing its interactions with them. Chymase complexes with the most active compound and substrate were used for development of a hybrid pharmacophore model which was applied in databases screening. Finally, hits which bound well at the active site, exhibited key interactions and favorable electronic properties were identified as possible inhibitors for chymase. This study not only elucidates inhibitory mechanism of chymase inhibitors but also provides key structural insights which will aid in the rational design of novel potent inhibitors of the enzyme. In general, the strategy applied in the current study could be a promising computational approach and may be generally applicable to drug design for other enzymes.

2-Cys peroxiredoxins (Prxs) play important roles in the protection of chloroplast proteins from oxidative damage. Arabidopsis NADPH-dependent thioredoxin reductase isotype C (AtNTRC) was identified as efficient electron donor for chloroplastic 2-Cys Prx-A. There are three isotypes (A, B, and C) of thioredoxin reductase (TrxR) in Arabidopsis. AtNTRA contains only TrxR domain, but AtNTRC consists of N-terminal TrxR and C-terminal thioredoxin (Trx) domains. AtNTRC has various oligomer structures, and Trx domain is important for chaperone activity. Our previous experimental study has reported that the hybrid protein (AtNTRA-(Trx-D)), which was a fusion of AtNTRA and Trx domain from AtNTRC, has formed variety of structures and shown strong chaperone activity. But, electron transfer mechanism was not detected at all. To find out the reason of this problem with structural basis, we performed two different molecular dynamics (MD) simulations on AtNTRC and AtNTRA-(Trx-D) proteins with same cofactors such as NADPH and flavin adenine dinucleotide (FAD) for 50 ns. Structural difference has found from superimposition of two structures that were taken relatively close to average structure. The main reason that AtNTRA-(Trx-D) cannot transfer the electron from TrxR domain to Trx domain is due to the difference of key catalytic residues in active site. The long distance between TrxR C153 and disulfide bond of Trx C387-C390 has been observed in AtNTRA-(Trx-D) because of following reasons: i) unstable and unfavorable interaction of the linker region, ii) shifted Trx domain, and iii) different or weak interface interaction of Trx domains. This study is one of the good examples for understanding the relationship between structure formation and reaction activity in hybrid protein. In addition, this study would be helpful for further study on the mechanism of electron transfer reaction in NADPH-dependent thioredoxin reductase proteins.

Sonic hedgehog (Shh) plays an important role in the activation of Shh signaling pathway that regulates preservation and rebirth of adult tissues. An abnormal activation of this pathway has been identified in hyperplasia and various tumorigenesis. Hence the inhibition of this pathway using a Shh inhibitor might be an efficient way to treat a wide range of malignancies. This study was done in order to develop a lead chemical candidate that has an inhibitory function in the Shh signaling pathway. We have generated common feature pharmacophore models using three-dimensional (3D) structural information of robotnikinin, an inhibitor of the Shh signaling pathway, and its analogs. These models have been validated with fit values of robotnikinin and its analogs, and the best model was used as a 3D structural query to screen chemical databases. The hit compounds resulted from the screening docked into a proposed binding site of the Shh named pseudo-active site. Molecular dynamics (MD) simulations were performed to investigate detailed binding modes and molecular interactions between the hit compounds and functional residues of the pseudo-active site. The results of the MD simulation analyses revealed that the hit compounds can bind the pseudo-active site with high affinity than robotnikinin. As a result of this study, a candidate inhibitor (GK 03795) was selected as a potential lead to be employed in future Shh inhibitor design.

Bacterial Hfq is a highly conserved thermostable protein of about 10 kDa. The Hfq protein was discovered in 1968 as an E. coli host factor that was essential for replication of the bacteriophage Q. It is now clear that Hfq has many important physiological roles. In E. coli, Hfq mutants show a multiple stress response related phenotypes. Hfq is now known to regulate the translation of two major stress transcription factors RpoS and RpoE in Enterobacteria and mediates its plieotrophic effects through several mechanisms. It interacts with regulatory sRNA and facilitates their antisense interaction with their targets. It also acts independently to modulate mRNA decay and in addition acts as a repressor of mRNA translation. Recent paper from Arluison et al. [9] provided the first evidence indicating that Hfq is an ATP-binding protein. They determined a plausible ATP-binding site in Hfq and tested Hfq’s ATP-binding affinity and stoichiometry. Experimental data suggest that the ATP-binding by the Hfq–RNA complex results in its significant destabilization of the protein and the result also proves important role of Tyr25 that flanks the cleft and stabilizes the adenine portion of ATP, possibly via aromatic stacking. In our study, the ATP molecule was docked into the predicted binding cleft using GOLD docking software. The binding nature of ATP and its effect on Hfq–RNA complex was studied using molecular dynamics simulations. Importance of Tyr25 residue was monitored and revealed using mutational study on the modeled systems. Our data and the corresponding results point to one of Hfq functional structural consequences due to ATP binding and Tyr25Ala mutation.

Histone deacetylases (HDACs) have emerged as effective therapeutic targets in the treatment of various diseases including cancers as these enzymes directly involved in the epigenetic regulation of genes. However the development of isoform-selective HDAC inhibitors has been a challenge till date since all HDAC enzymes possess conserved tunnel-like active site. In this study, using molecular dynamics simulation we have analyzed the behavior of tunnels present in HDAC8, 10, and 11 enzymes of class I, II, and IV, respectively. We have identified the equivalent tunnel forming amino acids in these three isoforms and found that they are very much conserved with subtle differences to be utilized in selective inhibitor development. One amino acid, methionine of HDAC8, among six tunnel forming residues is different in isoforms of other classes (glutamic acid (E) in HDAC10 and leucine (L) in HDAC 11) based on which mutations were introduced in HDAC11, the less studied HDAC isoform, to observe the effects of this change. The HDAC8-like (L268M) mutation in the tunnel forming residues has almost maintained the deep and narrow tunnel as present in HDAC8 whereas HDAC10-like (L268E) mutation has changed the tunnel wider and shallow as observed in HDAC10. These results explained the importance of the single change in the tunnel formation in different isoforms. The observations from this study can be utilized in the development of isoform-selective HDAC inhibitors.