Structure of the D142N mutant of the family 18 chitinase ChiB from Serratia marcescens and its complex with allosamidin

Catalysis by ChiB, a family 18 chitinase from Serratia marcescens, involves a conformational change of Asp142 which is part of a characteristic D₁₄₀XD₁₄₂XE₁₄₄ sequence motif. In the free enzyme Asp142 points towards Asp140, whereas it rotates towards the catalytic acid, Glu144, upon ligand binding. Mutation of Asp142 to Asn reduced k_cat and affinity for allosamidin, a competitive inhibitor. The X-ray structure of the D142N mutant showed that Asn142 points towards Glu144 in the absence of a ligand. The active site also showed other structural adjustments (Tyr10, Ser93) that had previously been observed in the wild-type enzyme upon substrate binding. The X-ray structure of a complex of D142N with allosamidin, a pseudotrisaccharide competitive inhibitor, was essentially identical to that of the wild-type enzyme in complex with the same compound. Thus, the reduced allosamidin affinity in the mutant is not caused by structural changes but solely by the loss of electrostatic interactions with Asp142. The importance of electrostatics was further confirmed by the pH dependence of catalysis and allosamidin inhibition. The pH-dependent apparent affinities for allosamidin were not correlated with k_cat, indicating that it is probably better to view the inhibitor as a mimic of the oxazolinium ion reaction intermediate than as a transition state analogue.     

Molecular Modeling-Based Evaluation of hTLR10 and Identification of Potential Ligands in Toll-Like Receptor Signaling

Toll-like receptors (TLRs) are pattern recognition receptors that recognize pathogens based on distinct molecular signatures. The human (h)TLR1, 2, 6 and 10 belong to the hTLR1 subfamilies, which are localized in the extracellular regions and activated in response to diverse ligand molecules. Due to the unavailability of the hTLR10 crystal structure, the understanding of its homo and heterodimerization with hTLR2 and hTLR1 and the ligand responsible for its activation is limited. To improve our understanding of the TLR10 receptor-ligand interaction, we used homology modeling to construct a three dimensional (3D) structure of hTLR10 and refined the model through molecular dynamics (MD) simulations. We utilized the optimized structures for the molecular docking in order to identify the potential site of interactions between the homo and heterodimer (hTLR10/2 and hTLR10/1). The docked complexes were then used for interaction with ligands (Pam₃CSK₄ and PamCysPamSK₄) using MOE-Dock and ASEDock. Our docking studies have shown the binding orientations of hTLR10 heterodimer to be similar with other TLR2 family members. However, the binding orientation of hTLR10 homodimer is different from the heterodimer due to the presence of negative charged surfaces at the LRR11-14, thereby providing a specific cavity for ligand binding. Moreover, the multiple protein-ligand docking approach revealed that Pam₃CSK₄ might be the ligand for the hTLR10/2 complex and PamCysPamSK_4, a di-acylated peptide, might activate hTLR10/1 hetero and hTLR10 homodimer. Therefore, the current modeled complexes can be a useful tool for further experimental studies on TLR biology.     

Homology modeling and substrate binding study of Nudix hydrolase Ndx1 from Thermos thermophilus HB8

With homology modeling techniques, molecular mechanics, and molecular dynamics methods, a 3D structure model of Ndx1 is created and refined. This model is further assessed by Profile-3D and ProStat, which confirm that the refined model is reliable. With this model, a flexible docking study is performed and the result indicates that Glu46, Arg88, and Glu90 are three important determinant residues in binding, as they have strong hydrogen bonding interactions and electrostatic interactions with Ap6A. In addition, we further find that three residues, Ser38, Leu39 and Glu46, coordinate enzyme-bound Mg2+ ions in complex N-A. The Glu46 is consistent with the experimental results by Iwai et al., and the other four residues mentioned above may also play vital roles in catalysis of Ndx1.     

Homology modelling and molecular dynamics study of human fatty acid amide hydrolase

The three-dimensional (3D) model of the human fatty acid amide hydrolase (hFAAH) was constructed based on the crystal structure of the rat FAAH (PDB code 1MT5) in complex with a substrate using Modeller9v2 program. With the aid of molecular mechanics and molecular dynamics method, the last model was obtained and further assessed by Profile-3D, Prosa2003 and Procheck, which confirms that the refined model is reliable. Furthermore, the docking results of propofol and its structural analogue into the active site of hFAAH indicate that 2,6-di-sec-butyl phenol is a more preferred ligand than others, which is in good agreement with the experimental results. From the docking studies, we also suggest that Phe192, Ile238, Thr377, Leu380, Phe381, Phe388 and Leu404 in the hFAAH are seven important determinant residues in binding as they have strong van der Waal interactions with the ligand.     

Structure of the D142N mutant of the family 18 chitinase ChiB from Serratia marcescens and its complex with allosamidin

Catalysis by ChiB, a family 18 chitinase from Serratia marcescens, involves a conformational change of Asp142 which is part of a characteristic D(140)XD(142)XE(144) sequence motif. In the free enzyme Asp142 points towards Asp140, whereas it rotates towards the catalytic acid, Glu144, upon ligand binding. Mutation of Asp142 to Asn reduced k(cat) and affinity for allosamidin, a competitive inhibitor. The X-ray structure of the D142N mutant showed that Asn142 points towards Glu144 in the absence of a ligand. The active site also showed other structural adjustments (Tyr10, Ser93) that had previously been observed in the wild-type enzyme upon substrate binding. The X-ray structure of a complex of D142N with allosamidin, a pseudotrisaccharide competitive inhibitor, was essentially identical to that of the wild-type enzyme in complex with the same compound. Thus, the reduced allosamidin affinity in the mutant is not caused by structural changes but solely by the loss of electrostatic interactions with Asp142. The importance of electrostatics was further confirmed by the pH dependence of catalysis and allosamidin inhibition. The pH-dependent apparent affinities for allosamidin were not correlated with k(cat), indicating that it is probably better to view the inhibitor as a mimic of the oxazolinium ion reaction intermediate than as a transition state analogue.     

Homology modeling and docking study of cyclin-dependent kinase (CDK) 10

In order to understand the mechanisms of ligand binding and the interaction between the ligand and the cyclin-dependent kinase 10 (CDK10), a three-dimensional (3D) model of the CDK10 is generated based on the crystal structure of the cyclin-dependent kinase 2 (CDK2) (PDB code 1AQ1) by using InsightII/Homology module. With the aid of the molecular mechanics and molecular dynamics methods, the last refined model is obtained and is further assessed by PROFILE-3D and PROSTAT, which show that the refined model is reliable. With this model, a flexible docking study is performed and the results indicate that the Lys39 and Asp94 form hydrogen bonds and have strong nonbonding interaction with adenosine 5'-triphosphate (ATP). From the docking studies, we also suggest that the Leu141, Tyr21, and Val24 in CDK10 are three important determinant residues in binding as they have strong nonbonding interaction with ATP. The hydrogen bonding interactions also play an important role for the stability of the complex. Our results may be helpful for further experimental investigations.     

Pharmacophore modeling,multiple docking,and molecular dynamics studies on Wee1 kinase inhibitors

Wee1-like protein kinase (Wee1) is a tyrosine kinase that regulates the G₂ checkpoint and prevents entry into mitosis in response to DNA damage. Based on a series of signaling pathways initiated by Wee1, Wee1 has been recognized as a potential target for cancer therapy. To discover potent Wee1 inhibitors with novel scaffolds, ligand-based pharmacophore model has been built based on 101 known Wee1 inhibitors. Then the best pharmacophore model, AADRRR.340, with good partial least square (PLS) statistics (R²?=?0.9212, Q²?=?0.7457), was selected and validated. The validated model was used as a three-dimensional (3D) search query for databases virtual screening. The filtered molecules were further analyzed and refined by Lipinski’s rule of 5, multiple docking procedures (high throughput virtual screening (HTVS), standard precision (SP), genetic optimization for ligand docking (GOLD), extra precision (XP), and unique quantum polarized ligand docking (QPLD)); absorption, distribution, metabolism, excretion, and toxicity (ADMET) screening; and the Prime/molecular mechanics generalized born surface area (MM-GBSA) method binding free energy calculations. Eight leads were identified as potential Wee1 inhibitors, and a 50?ns molecular dynamics (MD) simulation was carried out for top four inhibitors to predict the stability of ligand–protein complex. Molecular mechanics Poisson–Boltzmann surface area (MM-PBSA) based on MD simulation and the energy contribution per residue to the binding energy were calculated. In the end, three hits with good stabilization and affinity to protein were identified.

Communicated by Ramaswamy H. Sarma     

3D-QSAR pharmacophore modelling,virtual screening and docking studies for lead discovery of a novel scaffold for VEGFR 2 inhibitors: Design,synthesis and biological evaluation

A series of novel 6,7-dihydro-5H-cyclopenta[d]pyrimidine derivatives was successfully designed, synthesized and evaluated as a new chemical scaffold with vascular endothelial growth factor receptor (VEGFR 2) inhibitory activity. Compounds 6c and 6b showed enzyme inhibition of 97% and 87% at 10 µM, respectively, and exhibited potent dose-related VEGFR 2 inhibition with IC₅₀ values of 0.85 µM and 2.26 µM, respectively. The design of the 6,7-dihydro-5H-cyclopenta[d]pyrimidine scaffold was implemented via consecutive molecular modelling protocols prior to the synthesis and biological evaluation of the derivatives. First, sorafenib was docked in the binding site of VEGFR 2 to study its binding orientation and affinity, followed by the generation of a valid 3D QSAR pharmacophore model for use in the virtual screening of different 3D databases. Structures with promising pharmacophore-based virtual screening results were refined using molecular docking studies in the binding site of VEGFR 2. A novel scaffold was designed by incorporating the results of the pharmacophore model generation and molecular docking studies. The new scaffold showed hydrophobic interactions with the kinase front pocket that may be attributed to increasing residence time in VEGFR 2, which is a key success factor for ligand optimization in drug discovery. Different derivatives of the novel scaffold were validated using docking studies and pharmacophore mapping, where they exhibited promising results as VEGFR 2 inhibitors to be synthesized and biologically evaluated. 6,7-dihydro-5H-cyclopenta[d]pyrimidine is a new scaffold that can be further optimized for the synthesis of promising VEGFR 2 inhibitors.     

Homology modeling and PAPS ligand (cofactor) binding study of bovine phenol sulfotransferase

In order to understand the mechanisms of ligand binding and the interaction between the ligand and the bovine phenol sulfotransferase, (bSULT1A1, EC 2.8.2.1) a three-dimensional (3D) model of the bSULT1A1 is generated based on the crystal structure of the estrogen sulfotransferase (PDB code 1AQU) by using the InsightII/Homology module. With the aid of the molecular mechanics and molecular dynamics methods, the final refined model is obtained and is further assessed by Profile-3D and ProStat, which show that the refined model is reliable. With this model, a flexible docking study is performed and the results indicate that 3-phosphoadenosine-5- phosphosulfate (PAPS) is a more preferred ligand than coenzyme A (CoA), and that His108 forms hydrogen bond with PAPS, which is in good agreement with the experimental results. From these docking studies, we also suggest that Phe255, Phe24 and Tyr169 in bSULT1A1 are three important determinant residues in binding as they have strong van-der-Waals contacts with the ligand. The hydrogen–bonding interactions also play an important role for the stability of the complex. Our results may be helpful for further experimental investigations.Figure The final 3D-structure of bSULT1A1. The structure is obtained by energy minimizing an average conformation over the last 100 ps of MD simulation. The -helix is represented in red and the -sheet in yellow.     

Interaction Between Umami Peptide and Taste Receptor T1R1/T1R3

The umami taste receptor is a heterodimer composed of two members of the T1R taste receptor family: T1R1 and T1R3. The homology models of the ligand binding domains of the human umami receptor have been constructed based on crystallographic structures of the taste receptor of the central nervous system. Furthermore, the molecular simulations of the ligand binding domain show that the likely conformation was that T1R1 protein exists in the closed conformation, and T1R3 in the open conformation in the heterodimer. The molecular docking study of T1R1 and T1R3 in complex with four peptides, including Lys–Gly–Asp–Glu–Ser–Leu–Leu–Ala, Ser–Glu–Glu, G1u–Ser, and Asp–Glu–Ser, displayed that the amino acid residue of SER146 and Glu277 in T1R3 may play great roles in the synergism of umami taste. This docking result further validated the robustness of the model. In the paper, binding of umami peptide and the T1R1/T1R3 receptor was first described and the interaction is the base of umami activity theory.     

Homology modeling and molecular dynamics study on N-acetylneuraminate lyase

With homology modeling techniques, molecular mechanics and molecular dynamics methods, a 3D structure model of N-acetylneuraminate lyase from human (hNAL, EC 4.1.3.3) was created and refined. This model was further assessed by Profile-3D and PROCHECK, which confirms that the refined model is reliable. Furthermore, the docking results of the substrates (sialic acid and KDO) into the active site of hNAL indicate that hNAL can cleave the sialic acid and KDO. Thr51 and Tyr143 may be the key amino acids residues as they have strong hydrogen bonding interactions with the substrates, which is in good agreement with the experimental results by Izard et al. (Structure 2:361–369. doi:10.1016/S0969-2126(00)00038-1 (1994)). From the docking studies, we also suggest that Asp176 and Ser218 only form hydrogen bonds with sialic acid, therefore, they may help sialic acid interact with hNAL steadly.     

Homology modeling and atomic level binding study of Leishmania MAPK with inhibitors

The current therapy for leishmaniasis is not sufficient and it has two severe drawbacks, host-toxicity and drug resistance. The substantial knowledge of parasite biology is not yet translating into novel drugs for leishmaniasis. Based on this observation, a 3D structural model of Leishmania mitogen-activated protein kinase (MAPK) homologue has been developed, for the first time, by homology modeling and molecular dynamics simulation techniques. The model provided clear insight in its structure features, i.e. ATP binding pocket, phosphorylation lip, and common docking site. Sequence-structure homology recognition identified Leishmania CRK3 (LCRK3) as a distant member of the MAPK superfamily. Multiple sequence alignment and 3D structure model provided the putative ATP binding pocket of Leishmania with respect to human ERK2 and LCRK3. This analysis was helpful in identifying the binding sites and molecular function of the Leishmania specific MAPK homologue. Molecular docking study was performed on this 3D structural model, using different classes of competitive ATP inhibitors of LCRK3, to check whether they exhibit affinity and could be identified as Leishmania MAPK specific inhibitors. It is well known that MAP kinases are extracellular signal regulated kinases ERK1 and ERK2, which are components of the Ras-MAPK signal transduction pathway which is complexed with HDAC4 protein, and their inhibition is of significant therapeutic interest in cancer biology. In order to understand the mechanism of action, docking of indirubin class of molecules to the active site of histone deacetylase 4 (HDAC4) protein is performed, and the binding affinity of the protein-ligand interaction was computed. The new structural insights obtained from this study are all consistent with the available experimental data, suggesting that the homology model of the Leishmania MAPK and its ligand interaction modes are reasonable. Further the comparative molecular electrostatic potential and cavity depth analysis of Leishmania MAPK and human ERK2 suggested several important differences in its ATP binding pocket. Such differences could be exploited in the future for designing Leishmania specific MAPK inhibitors.     

Molecular docking and CoMFA studies of thiazoloquin(az)olin(on)es as CD38 inhibitors: determination of inhibitory mechanism,pharmacophore interactions,and design of new inhibitors

In this research, molecular docking and 3D-QSAR studies were carried out on a series of 79 thiazoloquin(az)olin(on)es as CD38 inhibitors. Based on docking results, four interactions including hydrogen bonding with main chain of GLU-226 (H-M-GLU-226), Van der Waals interactions with side chain of TRP-125 (V-S-TRP-125), TRP-189 (V-S-TRP-189), and THR-221 (V-S-THR-221) were considered as pharmacological interactions. Active conformation of each ligand was extracted from docking studies and was used for carrying out 3D-QSAR modeling. Comparative molecular field analysis (CoMFA) was performed on CD38 inhibitory activities of these compounds on human and mouse. We developed CoMFA models with five components as optimum models for both data-sets. For human data-set, a model with high predictive power was developed. R², RMSE, and F-test values for training set of this model were .94, .24, and 179.58, respectively, and R² and RMSE for its test set were .92 and .32, respectively. The q² and RMSE values for leave-one-out cross validation test on training set were .78 and .46, respectively, that demonstrate created model is robust. Based on extracted steric and electrostatic contour maps for this model, three inhibitors with pIC₅₀ larger than 8.85 were designed.     

Analysis of the structure and dynamics of human serum albumin

A computational approach to designing a peptide-based ligand for the purification of human serum albumin (HSA) was undertaken using molecular docking and molecular dynamics (MD) simulation. A three-step procedure was performed to design a specific ligand for HSA. Based on the candidate pocket structure of HSA (warfarin binding site), a peptide library was built. These peptides were then docked into the pocket of HSA using the GOLD program. The GOLDscore values were used to determine the affinity of peptides for HSA. Consequently, the dipeptide Trp–Trp, which shows a high GOLDscore value, was selected and linked to a spacer arm of Lys[CO(CH₂)₅NH] on the surface of ECH-lysine sepharose 4 gel. For further evaluation, the Autodock Vina program was used to dock the linked compound into the pocket of HSA. The docking simulation was performed to obtain a first guess of the binding structure of the spacer–Trp–Trp–HSA complex and subsequently analyzed by MD simulations to assess the reliability of the docking results. These MD simulations indicated that the ligand–HSA complex remains stable, and water molecules can bridge between the ligand and the protein by hydrogen bonds. Finally, absorption spectroscopic studies were performed to illustrate the appropriateness of the binding affinity of the designed ligand toward HSA. These studies demonstrate that the designed dipeptide can bind preferentially to the warfarin binding site. Graphical Abstract

Three-step computational approach to the design of a dipeptide ligand for human serum albumin purification exploiting structure-based docking and molecular dynamics simulation     

Can molecular dynamics simulations improve the structural accuracy and virtual screening performance of GPCR models?

The determination of G protein-coupled receptor (GPCR) structures at atomic resolution has improved understanding of cellular signaling and will accelerate the development of new drug candidates. However, experimental structures still remain unavailable for a majority of the GPCR family. GPCR structures and their interactions with ligands can also be modelled computationally, but such predictions have limited accuracy. In this work, we explored if molecular dynamics (MD) simulations could be used to refine the accuracy of in silico models of receptor-ligand complexes that were submitted to a community-wide assessment of GPCR structure prediction (GPCR Dock). Two simulation protocols were used to refine 30 models of the D₃ dopamine receptor (D₃R) in complex with an antagonist. Close to 60 μs of simulation time was generated and the resulting MD refined models were compared to a D₃R crystal structure. In the MD simulations, the receptor models generally drifted further away from the crystal structure conformation. However, MD refinement was able to improve the accuracy of the ligand binding mode. The best refinement protocol improved agreement with the experimentally observed ligand binding mode for a majority of the models. Receptor structures with improved virtual screening performance, which was assessed by molecular docking of ligands and decoys, could also be identified among the MD refined models. Application of weak restraints to the transmembrane helixes in the MD simulations further improved predictions of the ligand binding mode and second extracellular loop. These results provide guidelines for application of MD refinement in prediction of GPCR-ligand complexes and directions for further method development.     

Docking and SAR studies of calystegines: Binding orientation and influence on pharmacological chaperone effects for Gaucher’s disease

We report on the identification of the required configuration and binding orientation of nor-tropane alkaloid calystegines against β-glucocerebrosidase. Calystegine B₂ is a potent competitive inhibitor of human lysosomal β-glucocerebrosidase with K_i value of 3.3 μM. A molecular docking study revealed that calystegine B₂ had a favorable van der Waals interactions (Phe128, Trp179, and Phe246) and the hydrogen bonding (Glu235, Glu340, Asp127, Trp179, Asn234, Trp381 and Asn396) was similar to that of isofagomine. All calystegine isomers bound into the same active site as calystegine B₂ and the essential hydrogen bonds formed to Asp127, Glu235 and Glu340 were maintained. However, their binding orientations were obviously different. Calystegine A₃ bound to β-glucocerebrosidase with the same orientations as calystegine B₂ (Type 1), while calystegine B₃ and B₄ had different binding orientations (Type 2). It is noteworthy that Type 1 orientated calystegines B₂ and A₃ effectively stabilized β-glucocerebrosidase, and consequently increased intracellular β-glucocerebrosidase activities in N370S fibroblasts, while Type 2 orientated calystegines B₃ and B₄ could not keep the enzyme activity. These results clearly indicate that the binding orientations of calystegines are changed by the configuration of the hydroxyl groups on the nor-tropane ring and the suitable binding orientation is a requirement for achieving a strong affinity to β-glucocerebrosidase.     

Computational Studies of Ligand-Receptor Interactions in Bitter Taste Receptors

Phenylthiocarbamide tastes intensely bitter to some individuals, but others find it completely tasteless. Recently, it was suggested that phenylthiocarbamide elicits bitter taste by interacting with a human G protein-coupled receptor (hTAS2R38) encoded by the PTC gene. The phenylthiocarbamide nontaster trait was linked to three single nucleotide polymorphisms occurring in the PTC gene. Using the crystal structure of bovine rhodopsin as template, we generated the 3D structure of hTAS2R38 bitter taste receptor. We were able to map on the receptor structure the amino acids affected by the genetic polymorphisms and to propose molecular functions for two of them that explained the emergence of the nontaster trait. We used molecular docking simulations to find that phenylthiocarbamide exhibited a higher affinity for the target receptor than the structurally similar molecule 6-n-propylthiouracil, in line with recent experimental studies. A 3D model was constructed for the hTAS2R16 bitter taste receptor as well, by applying the same protocol. We found that the recently published experimental ligand binding affinity data for this receptor correlated well with the binding scores obtained from our molecular docking calculations.     

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

Abstract

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

Communicated by Ramaswamy H. Sarma     

Study of a ligand complexed with Cdk2/Cdk4 by computer simulation

Cyclin-dependent kinases (Cdks) play important roles in the regulation of the cell cycle. Their inhibitors have entered clinical trials to treat cancer. Very recently, Davis et al. (Nat Struct Biol 9:745–749, 2002) have found a ligand NU6102, which has a high affinity with cyclin-dependent kinase 2 (K _i =6 nM) but a low affinity with cyclin-dependent kinase 4 (K _i =1,600 nM). To understand the selectivity, we use homology modeling, molecular docking, molecular dynamics and free-energy calculations to analyze the interactions. A rational 3D model of the Cdk4–NU6102 complex is built. Asp86 is a key residue that recognizes NU6102 more effectively with Cdk2 rather than Cdk4. Good binding free energies are obtained. Energetic analysis reveals that van der Waals interaction and nonpolar contributions to solvent are favorable in the formation of complexes and the sulfonamide group of the ligand plays a crucial role for binding selectivity between Cdk2 and Cdk4. Figure Two-dimensional representative for the interacting model of NU6102 complexed with the Cdk4 from a predicted structure by LIGPLOT.      

Molecular modeling and docking characterization of CzR1, a CC-NBS-LRR R-gene from Curcuma zedoaria Loeb. that confers resistance to Pythium aphanidermatum

Plant NBS-LRR R-genes recognizes several pathogen associated molecular patterns (PAMPs) and limit pathogen infection through a multifaceted defense response. CzR1, a coiled-coil-nucleotide-binding-site-leucine-rich repeat R-gene isolated from Curcuma zedoaria L exhibit constitutive resistance to different strains of P. aphanidermatum. Majority of the necrotrophic oomycetes are characterized by the presence of carbohydrate PAMPs β-glucans in their cell walls which intercat with R-genes. In the present study, we predicted the 3D (three dimensional) structure of CzR1 based on homology modeling using the homology module of Prime through the Maestro interface of Schrodinger package ver 2.5. The docking investigation of CzR1 with β-glucan using the Glide software suggests that six amino acid residues, Ser186, Glu187, Ser263, Asp264, Asp355 and Tyr425 act as catalytic residues and are involved in hydrogen bonding with ligand β-(1,3)-D-Glucan. The calculated distance between the carboxylic oxygen atoms of Glu187–Asp355 pair is well within the distance of 5Å suggesting a positive glucanase activity of CzR1. Elucidation of these molecular characteristics will help in in silico screening and understanding the structural basis of ligand binding to CzR1 protein and pave new ways towards a broad spectrum rhizome rot resistance development in the cultivated turmeric.