Structure modeling, ligand binding, and binding affinity calculation (LR‐MM‐PBSA) of human heparanase for inhibition and drug design
- 13 September 2006
- journal article
- Published by Wiley in Proteins-Structure Function and Bioinformatics
- Vol. 65 (3) , 580-592
- https://doi.org/10.1002/prot.21065
Abstract
Human heparanase is an endo‐β‐D‐glycosidase that cleaves heparan sulphate (HS) chains in the extracellular matrix and basement membrane. It is known that the cleavage of HS by heparanase results in cell invasion and metastasis of cancer. Therefore, heparanase is considered an important target for cancer drug development. The three‐dimensional structure of heparanase would be useful in the rational design of inhibitors targeted to the enzyme; however, the three‐dimensional structure has not yet been determined. In our effort to design inhibitors, we developed a three‐dimensional structure of heparanase using a homology‐modeling approach. The homology‐built structure is consistent to previous bioinformatics and site‐mutation experimental results. The heparanase features a (α/β)8 TIM‐barrel fold with two glutamate residues (Glu225 and Glu343) located in the active‐site cleft. This feature supports the putative mechanism of proton donor and nucleophilic sites. Docking simulations yielded 41 complex structures, which indicate that the bound inhibitor could block ligand binding into the catalytic site. A free energy of binding model was established for 25 heparanase inhibitors with a training set of 25 heparanase inhibitors using the linear response MM‐PBSA approach (LR‐MM‐PBSA). The correlation between calculated and experimental activity was 0.79 and the reliability of the model was validated with leave‐one‐out cross‐validation method. Its predictive capability was further validated using a test set of 16 inhibitors similar to the training set of inhibitors. The correlation between the predicted and observed activities is significantly improved by the protein “induced‐fit” that accounts for the flexibility of the receptor. These interaction and pharmacophore elements provide a unique insight to the rational design of new ligands targeted to the enzyme. Proteins 2006.Keywords
This publication has 79 references indexed in Scilit:
- Homology modeling provides insights into the binding mode of the PAAD/DAPIN/pyrin domain, a fourth member of the CARD/DD/DED domain familyProtein Science, 2003
- Rapid grid‐based construction of the molecular surface and the use of induced surface charge to calculate reaction field energies: Applications to the molecular systems and geometric objectsJournal of Computational Chemistry, 2001
- Gapped BLAST and PSI-BLAST: a new generation of protein database search programsNucleic Acids Research, 1997
- AMBER, a package of computer programs for applying molecular mechanics, normal mode analysis, molecular dynamics and free energy calculations to simulate the structural and energetic properties of moleculesComputer Physics Communications, 1995
- Proteoglycans on endothelial cells present adhesion-inducing cytokines to leukocytesImmunology Today, 1993
- A method for predicting protein structure from sequenceCurrent Biology, 1993
- Assessment of protein models with three-dimensional profilesNature, 1992
- Molecular dynamics with coupling to an external bathThe Journal of Chemical Physics, 1984
- Comparison of simple potential functions for simulating liquid waterThe Journal of Chemical Physics, 1983
- A general method applicable to the search for similarities in the amino acid sequence of two proteinsJournal of Molecular Biology, 1970