pH modulates the TGF‐β ligands binding to the receptors: a computational analysis

In spite of showing high sequence similarity and forming structurally similar ternary complex in vitro, the in vivo role of TGF‐β1 and TGF‐β3 ligands suggests against their functional redundancy and necessitates the importance for the study of the specificity of these ligands. A comparative computational analysis of binary and ternary complexes of these two ligands shows that anchor residues of ligand and receptor at TGF‐β:TβR2 interface are similar in both complexes. However, the potential anchor residues of TGF‐β at TGF‐β:TβR1 interface are different, Tyr50 and Lys51 in TGF‐β3 complex and Lys60 and Tyr6 in TGF‐β1 complex. Pro55 and Asp57 of TβRI may act as anchor residues in complexes of both ligands along with Ile54 for TGF‐β3 complex and Val61 for TGF‐β1 complex. Arg58 of TβR1 acts as a potential hot residue for TGF‐β3 ternary complex but not for TGF‐β1 ternary complex formation whereas Pro55 and Phe60 may act as hot residues for both complexes. The Delphi analysis of the pH dependence of the binding energy indicates that pH has a remarkable effect on the binding energy of TβR2 to the open form of TGF‐β3. Lowering of pH from 7 to 4 favors binding of the open form of TGF‐β3 to TβR2. Now, apart from the residues at pH 7, residues Arg25, Lys31 and Arg94 of TGF‐β3 and Asp118 and Glu119 of TβR2 also contribute significantly to the binding energy. Contrary to the binding energy of TβR2 to TGF‐β3/TGF‐β1, TβR1 shows appreciable pH dependence for its binding in ternary complex of TGF‐β3/TGF‐β1. In TGF‐β3 ternary complex, the TβR1 electrostatic interaction energy disfavors complex formation at pH 7 while it is favored at pH 4. Copyright © 2014 John Wiley & Sons, Ltd.     

The structure of the C-terminal domain of the largest editosome interaction protein and its role in promoting RNA binding by RNA-editing ligase L2

Trypanosomatids, such as the sleeping sickness parasite Trypanosoma brucei, contain a ～ 20S RNA-editing complex, also called the editosome, which is required for U-insertion/deletion editing of mitochondrial mRNAs. The editosome contains a core of 12 proteins including the large interaction protein A1, the small interaction protein A6, and the editing RNA ligase L2. Using biochemical and structural data, we identified distinct domains of T. brucei A1 which specifically recognize A6 and L2. We provide evidence that an N-terminal domain of A1 interacts with the C-terminal domain of L2. The C-terminal domain of A1 appears to be required for the interaction with A6 and also plays a key role in RNA binding by the RNA-editing ligase L2 in trans. Three crystal structures of the C-terminal domain of A1 have been elucidated, each in complex with a nanobody as a crystallization chaperone. These structures permitted the identification of putative dsRNA recognition sites. Mutational analysis of conserved residues of the C-terminal domain identified Arg703, Arg731 and Arg734 as key requirements for RNA binding. The data show that the editing RNA ligase activity is modulated by a novel mechanism, i.e. by the trans-acting RNA binding C-terminal domain of A1.     

Structure-based and ligand-based drug design for microsomal prostaglandin E synthase-1 inhibitors

Microsomal prostaglandin E synthase-1 (mPGES-1) has been regarded as an attractive drug for inflammation-related diseases. In search of new mPGES-1 inhibitors, we performed virtual screening using our traditional Chinese medicine and natural products database (http://tcm.cmu.edu.tw/) and constructed comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) using a training set of 30 experimentally tested mPGES-1 inhibitors. The CoMFA and CoMSIA models derived were statistically significant with cross-validated coefficient values of 0.808 for CoMFA and 0.829 for CoMSIA and non-cross-validated coefficient values of 0.829 for CoMFA and 0.980 for CoMSIA. Docking and de novo evolution design gave three top derivatives, 2-O-caffeoyl tartaric acid-Evo_2, glucogallin-Evo_1 and 3-O-feruloylquinic acid-Evo_7 that have higher binding affinities than the control, glutathione. These three derivatives have interactions with Arg70, Arg73, Arg110, Arg126 and Arg38, which all are mPGES-1 key active site residues. In addition, these derivatives fit well into the CoMFA and CoMSIA models, with hydrophobic, hydrophilic and electropositive substructures mapped onto corresponding contour plots. Hence, we suggest that these three de novo compounds could be a starting basis for new mPGES-1 inhibitors.     

Identification of hotspot regions of MurB oxidoreductase enzyme using homology modeling,molecular dynamics and molecular docking techniques

Despite the availability of effective chemotherapy and a moderately protective vaccine, new anti-tuberculosis agents are urgently needed to decrease the global incidence of tuberculosis (TB) disease. The MurB gene belongs to the bacterial cell wall biosynthesis pathway and is an essential drug target in Mycobacterium tuberculosis (Mtb) that has no mammalian counterparts. Here, we present an integrated approach involving homology modeling, molecular dynamics and molecular docking studies on Mtb-MurB oxidoreductase enzyme. A homology model of Mtb-MurB enzyme was built for the first time in order to carry out structure-based inhibitor design. The accuracy of the model was validated using different techniques. The molecular docking study on this enzyme was undertaken using different classes of well known MurB inhibitors. Estimation of binding free energy by docking analysis indicated the importance of Tyr155, Arg156, Ser237, Asn241 and His304 residues within the Mtb-MurB binding pocket. Our computational analysis is in good agreement with experimental results of site-directed mutagenesis. The present study should therefore play a guiding role in the experimental design of Mtb-MurB inhibitors for in vitro/in vivo analysis.     

Evidence for subsites in the galectins involved in sugar binding at the nonreducing end of the central galactose of oligosaccharide ligands: sequence analysis, homology modeling and mutagenesis studies of hamster galectin-3

A model of the carbohydrate recognition domain CRD, residues 111-245, of hamster galectin-3 has been made using homology modeling and dynamics minimization methods. The model is based on the known x-ray structures of bovine galectin-1 and human galectin-2. The oligosaccharides NeuNAc-alpha2,3-Gal-beta1,4-Glc and GalNAc-alpha1, 3- [Fuc-alpha1,2]-Gal-beta1,4-Glc, known to be specific high-affinity ligands for galectin-3, as well as lactose recognized by all galectins were docked in the galectin-3 CRD model structure and a minimized binding conformation found in each case. These studies indicate a putative extended carbohydrate-binding subsite in the hamster galectin- 3 involving Arg139, Glu230, and Ser232 for NeuNAc-alpha2,3-; Arg139 and Glu160 for fucose-alpha1,2-; and Arg139 and Ile141 for GalNAc-alpha1,3- substituents on the primary galactose. Each of these positions is variable within the whole galectin family. Two of these residues, Arg139 and Ser232, were selected for mutagenesis to probe their importance in this newly identified putative subsite. Residue 139 adopts main-chain dihedral angles characteristic of an isolated bridge structural feature, while residue 232 is the C-terminal residue of beta- strand-11, and is followed immediately by an inverse gamma-turn. A systematic series of mutant proteins have been prepared to represent the residue variation present in the aligned sequences of galectins-1, - 2, and -3. Minimized docked models were generated for each mutant in complex with NeuNAc-alpha2,3-Gal-beta1,4-Glc, GalNAc-alpha1, 3-[Fuc- alpha1,2]-Gal-beta1,4- Glc, and Gal-beta1,4-Glc. Correlation of the computed protein-carbohydrate interaction energies for each lectin- oligosaccharide pair with the experimentally determined binding affinities for fetuin and asialofetuin or the relative potencies of lactose and sialyllactose in inhibiting binding to asiolofetuin is consistent with the postulated key importance of Arg139 in recognition of the extended sialylated ligand.      

Origins of the specificity of inhibitor P218 toward wild-type and mutant PfDHFR: a molecular dynamics analysis

Molecular dynamics simulations were performed to evaluate the origin of the antimalarial effect of the lead compound P218. The simulations of the ligand in the cavities of wild-type, mutant Plasmodium falciparum Dihydrofolate Reductase (PfDHFR) and the human DHFR revealed the differences in the atomic-level interactions and also provided explanation for the specificity of this ligand toward PfDHFR. The binding free energy estimation using Molecular Mechanics Poisson-Boltzmann Surface Area method revealed that P218 has higher binding affinity (~ ?30 to ?35 kcal/mol) toward PfDHFR (both in wild-type and mutant forms) than human DHFR (~ ?22 kcal/mol), corroborating the experimental observations. Intermolecular hydrogen bonding analysis of the trajectories showed that P218 formed two stable hydrogen bonds with human DHFR (Ile7 and Glu30), wild-type and double-mutant PfDHFR’s (Asp54 and Arg122), while it formed three stable hydrogen bonds with quadruple-mutant PfDHFR (Asp54, Arg59, and Arg122). Additionally, P218 binding in PfDHFR is stabilized by hydrogen bonds with residues Ile14 and Ile164. It was found that mutant residues do not reduce the binding affinity of P218 to PfDHFR, in contrast, Cys59Arg mutation strongly favors inhibitor binding to quadruple-mutant PfDHFR. The atomistic-level details explored in this work will be highly useful for the design of non-resistant novel PfDHFR inhibitors as antimalarial agents.     

Structure-function analysis of the kinase-CPD domain of yeast tRNA ligase (Trl1) and requirements for complementation of tRNA splicing by a plant Trl1 homolog

Trl1 is an essential 827 amino acid enzyme that executes the end-healing and end-sealing steps of tRNA splicing in Saccharomyces cerevisiae. Trl1 consists of two domains—an N-terminal ligase component and a C-terminal 5′-kinase/2′,3′-cyclic phosphodiesterase (CPD) component—that can function in tRNA splicing in vivo when expressed as separate polypeptides. To understand the structural requirements for the kinase-CPD domain, we performed an alanine scan of 30 amino acids that are conserved in Trl1 homologs from other fungi. We thereby identified four residues (Arg463, His515, Thr675 and Glu741) as essential for activity in vivo. Structure–function relationships at these positions, and at four essential or conditionally essential residues defined previously (Asp425, Arg511, His673 and His777), were clarified by introducing conservative substitutions. Biochemical analysis showed that lethal mutations of Asp425, Arg463, Arg511 and His515 in the kinase module abolished polynucleotide kinase activity in vitro. We report that a recently cloned 1104 amino acid Arabidopsis RNA ligase functions in lieu of yeast Trl1 in vivo and identify essential side chains in the ligase, kinase and CPD modules of the plant enzyme. The plant ligase, like yeast Trl1 but unlike T4 RNA ligase 1, requires a 2′-PO₄ end for tRNA splicing in vivo.     

Insights into the binding specificity of wild type and mutated wheat germ agglutinin towards Neu5Acα(2‐3)Gal: a study by in silico mutations and molecular dynamics simulations

Wheat germ agglutinin (WGA) is a plant lectin, which specifically recognizes the sugars NeuNAc and GlcNAc. Mutated WGA with enhanced binding specificity can be used as biomarkers for cancer. In silico mutations are performed at the active site of WGA to enhance the binding specificity towards sialylglycans, and molecular dynamics simulations of 20 ns are carried out for wild type and mutated WGAs (WGA1, WGA2, and WGA3) in complex with sialylgalactose to examine the change in binding specificity. MD simulations reveal the change in binding specificity of wild type and mutated WGAs towards sialylgalactose and bound conformational flexibility of sialylgalactose. The mutated polar amino acid residues Asn114 (S114N), Lys118 (G118K), and Arg118 (G118R) make direct and water mediated hydrogen bonds and hydrophobic interactions with sialylgalactose. An analysis of possible hydrogen bonds, hydrophobic interactions, total pair wise interaction energy between active site residues and sialylgalactose and MM‐PBSA free energy calculation reveals the plausible binding modes and the role of water in stabilizing different binding modes. An interesting observation is that the binding specificity of mutated WGAs (cyborg lectin) towards sialylgalactose is found to be higher in double point mutation (WGA3). One of the substituted residues Arg118 plays a crucial role in sugar binding. Based on the interactions and energy calculations, it is concluded that the order of binding specificity of WGAs towards sialylgalactose is WGA3 > WGA1 > WGA2 > WGA. On comparing with the wild type, double point mutated WGA (WGA3) exhibits increased specificity towards sialylgalactose, and thus, it can be effectively used in targeted drug delivery and as biological cell marker in cancer therapeutics. Copyright © 2014 John Wiley & Sons, Ltd.     

Identification of ligand effector binding sites in transmembrane regions of the human G protein-coupled C3a receptor

The human C3a anaphylatoxin receptor (C3aR) is a G protein-coupled receptor (GPCR) composed of seven transmembrane alpha-helices connected by hydrophilic loops. Previous studies of chimeric C3aR/C5aR and loop deletions in C3aR demonstrated that the large extracellular loop2 plays an important role in noneffector ligand binding; however, the effector binding site for C3a has not been identified. In this study, selected charged residues in the transmembrane regions of C3aR were replaced by Ala using site-directed mutagenesis, and mutant receptors were stably expressed in the RBL-2H3 cell line. Ligand binding studies demonstrated that R161A (helix IV), R340A (helix V), and D417A (helix VII) showed no binding activity, although full expression of these receptors was established by flow cytometric analysis. C3a induced very weak intracellular calcium flux in cells expressing these three mutant receptors. H81A (helix II) and K96A (helix III) showed decreased ligand binding activity. The calcium flux induced by C3a in H81A and K96A cells was also consistently reduced. These findings suggest that the charged transmembrane residues Arg161, Arg340, and Asp417 in C3aR are essential for ligand effector binding and/or signal coupling, and that residues His81 and Lys96 may contribute less directly to the overall free energy of ligand binding. These transmembrane residues in C3aR identify specific molecular contacts for ligand interactions that account for C3a-induced receptor activation.     

Binding studies with mutants of Zif268. Contribution of individual side chains to binding affinity and specificity in the Zif268 zinc finger-DNA complex.

The Zif268 zinc finger-DNA complex has served as a model system for understanding how Cys2His2 type zinc fingers recognize DNA. Structural studies of the Zif268-DNA complex revealed that residues at four positions in the alpha helix of each zinc finger play key roles in recognition, but there has been no information about the precise contributions of individual residues. Here we report the results of binding studies involving five mutants of Zif268 that have changes in the base-contacting residues of finger one. These studies let us evaluate the contributions that Arg18 (position -1 of the alpha helix), Asp20 (position 2), Glu21 (position 3), and Arg24 (position 6) make to the overall energy of DNA binding. Our results confirm the important role played by these arginines. By comparing the affinities of the wild type and mutant peptides for various sites, we also prove that Asp20 and Glu21 play important roles in determining binding site specificity.     

Influence of cation-pi interactions on RNA-binding proteins

The energy contribution due to cation-π interactions has been computed for 37 RNA binding proteins. The contribution of these cation-π interacting residues in sequential separation, secondary structure involvement, solvent accessibility, and stabilization centers has been evaluated. Sequential separation of the cation-π involving residues show that, long range contacts predominates in all the proteins studied. Lys and Arg prefers to be in helical structures. Of the cation-π interacting residues, Arg and Lys were in the exposed regions and the aromatic residues (Phe, Tyr and Trp) were in the buried and partially buried regions in the protein structures. Stabilization centers for these proteins showed that all the five residues found in cation-π interactions are important in locating one or more of such centers. On the whole, the results presented in this work will be very useful for further investigations on the specificity and selectivity of RNA binding proteins and also for their structural studies.     

Molecular dynamics simulation studies of GSK-3β ATP competitive inhibitors: understanding the factors contributing to selectivity

Glycogen synthase kinase-3 is a constitutively acting, multifunctional serine threonine kinase, the role of which has been implicated in several physiological pathways and has emerged as a promising target for the treatment of type-II diabetes and Alzheimer’s disease. In order to provide a detailed understanding of the origin of selectivity determinants of ATP competitive inhibitors, molecular dynamics simulations in combination with MM-PBSA binding energy calculations were performed using crystal structures of GSK-3β and CDK-2 in complex with 12 ATP competitive inhibitors. Analysis of energy contributions indicate that electrostatic interaction energy dictates the selectivity of ATP competitive inhibitors against CDK-2. Key interactions as well as residues that potentially make a major contribution to the binding free energy were identified at the ATP binding site. This analysis stresses the need for the inhibitors to interact with Lys85, Thr138, and Arg141 in the binding site of GSK-3β to show selectivity. The residue-wise energy decomposition analysis further suggested the additional role of Gln185 in determining the selectivity of maleimides. The results obtained in this study can be utilized to design new selective GSK-3 ATP competitive inhibitors.     

Molecular dynamics (MD) simulations and binding free energy calculation studies between inhibitors and type II dehydroquinase (DHQ2)

Type II dehydroquinase (DHQ2) is the third enzyme of the shikimic acid pathway, and it has been the effective target for tuberculosis (TB). So far, developing multiple potent inhibitors of the DHQ2 of Mycobacterium tuberculosis (DHQ2-Mt) has been considered to be the new therapy to TB. Molecular dynamics simulations followed by molecular mechanics-generalised Born surface area were carried out to calculate the free binding energy and to determine the affinity ability of the four chosen inhibitor molecules, L1, L2, L3 and L4. Energy decomposition analyses show the electrostatic interaction and van der Waals interaction of the ligands to every residue of the DHQ2-Mt. The results suggest that some important residues have different interactions with the four ligands, such as Arg19 and Tyr24. These interactions may have an effect on the ligand binding affinity. The binding affinity of monosubstituted inhibitor is higher than that of disubstituted inhibitor, due to some important interactions with the DHQ2-Mt residues. These computational works will be significant to the theoretical research in the future.     

The systematic modeling studies and free energy calculations of the phenazine compounds as anti-tuberculosis agents

Abstract

Phenazine compounds have good activity against Mycobacterium tuberculosis (MTB). Based on the reported activities that were obtained in MTB H37Rv, a three-dimensional quantitative structure–activity relationship (3D-QSAR) model was built to design novel compounds against MTB. A fivefold cross-validation method and external validation were used to analyze the accuracy of forecasting. The model has a cross-validation coefficient q²=0.7 and a non-cross-validation coefficient r² = 0.903, indicating that the model has good predictive possibility. The design of anti-pneumococcus MTB compounds was guided by the obtained 3D-QSAR model, and several compounds with better activity were obtained. To test the activity of these compounds, molecular docking, molecular dynamics simulation, and post-simulation analysis of the already reported drug targets in MTB were carried out. Among the total 15 drug targets, only three targets (Rv2361c, Rv2965c, and Rv3048c) were selected based on the docking results. Initial results reported that these compounds possessed good inhibition activity for Rv2361c. The top nine complexes of Rv2361 ligands were only subjected to MD simulation which resulted in a stable dynamics of the structures and showed a residual fluctuation in inhibitors binding pocket. Free energy reported that overall, the derivatives hold strong energy against the protein target. Energetic contribution results showed that residues, Asp76, Arg80, Asn124, Arg127, Arg244, and Arg250, play a major role in total energy. Systems biology approach validates shortlisted drug effect on the entire system which might be useful to predict potential drug in wet lab as well.

Communicated by Ramaswamy H. Sarma     

Mutagenic analysis of conserved arginine residues in and around the novel sulfate binding pocket of the human Theta class glutathione transferase T2-2.

The human Theta class glutathione transferase GSTT2-2 has a novel sulfatase activity that is not dependent on the presence of a conserved hydrogen bond donor in the active site. Initial homology modeling and the crystallographic studies have identified three conserved Arg residues that contribute to the formation of (Arg107 and Arg239), and entry to (Arg242), a sulfate binding pocket. These residues have been individually mutated to Ala to investigate their potential role in substrate binding and catalysis. The mutation of Arg107 had a significant detrimental effect on the sulfatase reaction, while the Arg242 mutation caused only a small reduction in sulfatase activity. Surprisingly, the Arg239 had an increased activity resulting from a reduction in stability. Thus, Arg239 appears to play a role in maintaining the architecture of the active site. Electrostatic calculations performed on the wild-type and mutant forms of the enzyme are in good agreement with the experimental results. These findings, along with docking studies, suggest that prior to conjugation, the location of 1-menaphthyl sulfate, a model substrate for the sulfatase reaction, is approximately midway between the position ultimately occupied by the naphthalene ring of 1-menaphthylglutathione and the free sulfate. It is further proposed that the Arg residues in and around the sulfate binding pocket have a role in electrostatic substrate recognition.     

Mutagenesis of the AT1 receptor reveals different binding modes of angiotensin II and [Sar1]-angiotensin II

Homology modeling of the structure of the AT1 receptor, based on the high resolution rhodopsin crystal structure, indicated that it is unlikely that the binding of AngII to AT1 involves simultaneously all the receptor's residues reported in the literature to participate in this process. Site-directed mutagenesis using Ala substitution of charged residues Lys20, Arg23, Glu91 and Arg93 was performed to evaluate the participation of their side-chains in ligand binding and in triggering the cell's response. A comparative analysis by competition binding and functional assays using angiotensin II and the analog [Sar1]-angiotensin II suggests an important role for Arg23 of AT1 receptor in binding of the natural agonist. It is discussed whether some receptor's residues participate directly in the binding with AngII or whether they are part of a regulatory site.     

The maintenance of high affinity plasminogen binding by group A streptococcal plasminogen-binding M-like protein is mediated by arginine and histidine residues within the a1 and a2 repeat domains

Subversion of the plasminogen activation system is implicated in the virulence of group A streptococci (GAS). GAS displays receptors for the human zymogen plasminogen on the cell surface, one of which is the plasminogen-binding group A streptococcal M-like protein (PAM). The plasminogen binding domain of PAM is highly variable, and this variation has been linked to host selective immune pressure. Site-directed mutagenesis of full-length PAM protein from an invasive GAS isolate was undertaken to assess the contribution of residues in the a1 and a2 repeat domains to plasminogen binding function. Mutagenesis to alanine of key plasminogen binding lysine residues in the a1 and a2 repeats (Lys98 and Lys111) did not abrogate plasminogen binding by PAM nor did additional mutagenesis of Arg101 and His102 and Glu104, which have previously been implicated in plasminogen binding. Plasminogen binding was only abolished with the additional mutagenesis of Arg114 and His115 to alanine. Furthermore, mutagenesis of both arginine (Arg101 and Arg114) and histidine (His102 and His115) residues abolished interaction with plasminogen despite the presence of Lys98 and Lys111 in the binding repeats. This study shows for the first time that residues Arg101, Arg114, His102, and His115 in both the a1 and a2 repeat domains of PAM can mediate high affinity plasminogen binding. These data suggest that highly conserved arginine and histidine residues may compensate for variation elsewhere in the a1 and a2 plasminogen binding repeats, and may explain the maintenance of high affinity plasminogen binding by naturally occurring variants of PAM.     

Molecular dynamics and free energy studies on the Drosophila melanogaster and Leptinotarsa decemlineata ecdysone receptor complexed with agonists: Mechanism for binding and selectivity

The ecdysone receptor is a nuclear hormone receptor that plays a pivotal role in the insect metamorphosis and development. To address the molecular mechanisms of binding and selectivity, the interactions of two typical agonists Ponasterone A and 20-Hydroxyecdysone with Drosophila melanogaster (DME) and Leptinotarsa decemlineata ecdysone (LDE) receptors were investigated by homology modeling, molecular docking, molecular dynamic simulation, and thermodynamic analysis. We discover that 1) the L5-loop, L11-loop, and H12 helix for DME, L7-loop, and L11-loop for LDE are more flexible, which affect the global dynamics of the ligand-binding pocket, thus facilitating the ligand recognition of ecdysone receptor; 2) several key residues (Thr55/Thr37, Phe109/Phe91, Arg95/Arg77, Arg99/Arg81, Phe108/Leu90, and Ala110/Val92) are responsible for the binding of the proteins; 3) the binding-free energy is mainly contributed by the van der Waals forces as well as the electrostatic interactions of ligand and receptor; 4) the computed binding-free energy difference between DME-C1 and LDE-C1 is –4.65 kcal/mol, explains that C1 can form many more interactions with the DME; 5) residues Phe108/Leu90 and Ala110/Val92 have relatively position and orientation difference in the two receptors, accounting most likely for the ligand selectivity of ecdysone receptor from different orders of insects. This study underscores the expectation that different insect pests should be able to discriminate among compounds from different as yet undiscovered compounds, and the results firstly show a structural and functional relay between the agonists and receptors (DME and LDE), which can provide an avenue for the development of target-specific insecticides.

Communicated by Ramaswamy H. Sarma     

In silico modeling and hydrogen peroxide binding study of rice catalase

Homology modeling of the catalase, CatC cloned and sequenced from rice (Oryza sativa L., cv Ratna an Indica cultivar) has been performed based on the crystal structure of the catalase CatF (PDB code 1m7s) by using the software MODELLER. With the aid of molecular mechanics and molecular dynamics methods, the final model is obtained and is further assessed by PROCHECK and VERIFY - 3D graph, which show that the final refined model is reliable. With this model, a flexible docking study with the hydrogen peroxide, the substrate for catalase, is performed and the results indicate that Arg310, Asp343 and Arg346 in catalase are three important determinant residues in binding as they have strong hydrogen bonding contacts with the substrate. These hydrogen-bonding interactions play an important role for the stability of the complex. Our results may be helpful for further experimental investigations.