Molecular Dynamics Simulations of Sonic Hedgehog-Receptor and Inhibitor Complexes and Their Applications for Potential Anticancer Agent Discovery

The sonic hedgehog (Shh) signaling pathway is necessary for a variety of development and differentiation during embryogenesis as well as maintenance and renascence of diverse adult tissues. However, an abnormal activation of the signaling pathway is related to various cancers. In this pathway, the Shh signaling transduction is facilitated by binding of Shh to its receptor protein, Ptch. In this study, we modeled the 3D structure of functionally important key loop peptides of Ptch based on homologous proteins. Using this loop model, the molecular interactions between the structural components present in the pseudo-active site of Shh and key residues of Ptch was investigated in atomic level through molecular dynamics (MD) simulations. For the purpose of developing inhibitor candidates of the Shh signaling pathway, the Shh pseudo-active site of this interface region was selected as a target to block the direct binding between Shh and Ptch. Two different structure-based pharmacophore models were generated considering the key loop of Ptch and known inhibitor-induced conformational changes of the Shh through MD simulations. Finally two hit compounds were retrieved through a series of virtual screening combined with molecular docking simulations and we propose two hit compounds as potential inhibitory lead candidates to block the Shh signaling pathway based on their strong interactions to receptor or inhibitor induced conformations of the Shh.     

Syntheses of aminoalcohol-derived macrocycles leading to a small-molecule binder to and inhibitor of Sonic Hedgehog

We report the synthesis and biological activity of a library of aminoalcohol-derived macrocycles from which robotnikinin (17), a binder to and inhibitor of Sonic Hedgehog, was derived. Using an asymmetric alkylation to set a key stereocenter and an RCM reaction to close the macrocycle, we were able to synthesize compounds for testing. High-throughput screening via small-molecule microarray (SMM) technology led to the discovery of a compound capable of binding ShhN. Follow-up chemistry led to a library of macrocycles with enhanced biological activity relative to the original hit compounds. Differences in ring size and stereochemistry, leading to alterations in the mode of binding, may account for differences in the degree of biological activity. These compounds are the first ones reported that inhibit Shh signaling at the ShhN level.     

Hedgehog Pathway Antagonist 5E1 Binds Hedgehog at the Pseudo-active Site

Proper hedgehog (Hh) signaling is crucial for embryogenesis and tissue regeneration. Dysregulation of this pathway is associated with several types of cancer. The monoclonal antibody 5E1 is a Hh pathway inhibitor that has been extensively used to elucidate vertebrate Hh biology due to its ability to block binding of the three mammalian Hh homologs to the receptor, Patched1 (Ptc1). Here, we engineered a murine:human chimeric 5E1 (ch5E1) with similar Hh-binding properties to the original murine antibody. Using biochemical, biophysical, and x-ray crystallographic studies, we show that, like the regulatory receptors Cdon and Hedgehog-interacting protein (Hhip), ch5E1 binding to Sonic hedgehog (Shh) is enhanced by calcium ions. In the presence of calcium and zinc ions, the ch5E1 binding affinity increases 10–20-fold to tighter than 1 nm primarily because of a decrease in the dissociation rate. The co-crystal structure of Shh bound to the Fab fragment of ch5E1 reveals that 5E1 binds at the pseudo-active site groove of Shh with an epitope that largely overlaps with the binding site of its natural receptor antagonist Hhip. Unlike Hhip, the side chains of 5E1 do not directly coordinate the Zn²⁺ cation in the pseudo-active site, despite the modest zinc-dependent increase in 5E1 affinity for Shh. Furthermore, to our knowledge, the ch5E1 Fab-Shh complex represents the first structure of an inhibitor antibody bound to a metalloprotease fold.     

Exploration of Virtual Candidates for Human HMG-CoA Reductase Inhibitors Using Pharmacophore Modeling and Molecular Dynamics Simulations

3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGR) is a rate-controlling enzyme in the mevalonate pathway which involved in biosynthesis of cholesterol and other isoprenoids. This enzyme catalyzes the conversion of HMG-CoA to mevalonate and is regarded as a drug target to treat hypercholesterolemia. In this study, ten qualitative pharmacophore models were generated based on chemical features in active inhibitors of HMGR. The generated models were validated using a test set. In a validation process, the best hypothesis was selected based on the statistical parameters and used for virtual screening of chemical databases to find novel lead candidates. The screened compounds were sorted by applying drug-like properties. The compounds that satisfied all drug-like properties were used for molecular docking study to identify their binding conformations at active site of HMGR. The final hit compounds were selected based on docking score and binding orientation. The HMGR structures in complex with the hit compounds were subjected to 10 ns molecular dynamics simulations to refine the binding orientation as well as to check the stability of the hits. After simulation, binding modes including hydrogen bonding patterns and molecular interactions with the active site residues were analyzed. In conclusion, four hit compounds with new structural scaffold were suggested as novel and potent HMGR inhibitors.     

Sonic hedgehog promotes proliferation and differentiation of adult muscle cells: Involvement of MAPK/ERK and PI3K/Akt pathways

Sonic hedgehog (Shh) has been reported to act as a mitogen and survival factor for muscle satellite cells. However, its role in their differentiation remains ambiguous. Here, we provide evidence that Shh promotes the proliferation and differentiation of primary cultures of chicken adult myoblasts (also termed satellite cells) and mouse myogenic C2 cells. These effects are reversed by cyclopamine, a specific chemical inhibitor of the Shh pathway. In addition, we show that Shh and its downstream molecules are expressed in adult myoblast cultures and localize adjacent to Pax7 in muscle sections. These gene expressions are regulated during postnatal muscle growth in chicks. Most importantly, we report that Shh induces MAPK/ERK and phosphoinositide 3-kinase (PI3K)-dependent Akt phosphorylation and that activation of both signaling pathways is essential for Shh's signaling in muscle cells. However, the effect of Shh on Akt phosphorylation is more robust than that on MAPK/ERK, and data suggest that Shh influences these pathways in a manner similar to IGF-I. By exploiting specific chemical inhibitors of the MAPK/ERK and PI3K/Akt signaling pathways, UO126 and Ly294002, respectively, we demonstrate that Shh-induced Akt phosphorylation, but not that of MAPK/ERK, is required for its promotive effects on muscle cell proliferation and differentiation. Taken together, we suggest that Shh acts in an autocrinic manner in adult myoblasts, and provide first evidence of a role for PI3K/Akt in Shh signaling during myoblast differentiation.     

Crystallographic analysis at 3.0-A resolution of the binding to human thrombin of four active site-directed inhibitors

The mode of binding of four active-site directed inhibitors to human thrombin has been determined by x-ray crystallographic analysis. The inhibitors studied are benzamidine, PPACK, NAPAP, and MD-805, of which the last three are compounds evolved specifically to inhibit thrombin. Crystal structures were determined in the presence of both the inhibitor and the undecapeptide [des-amino Asp55]hirudin(55-65) which binds distant from the active site. Despite having significantly different chemical structures, NAPAP and MD-805 bind to thrombin in a very similar "inhibitor binding mode" which is not that expected by direct analogy with the binding of substrate. Both inhibitors bind to thrombin in a similar way as to trypsin, but thrombin has an extra loop, the "Tyr-Pro-Pro-Trp loop," not present in trypsin, which gives further binding interactions and is seen to move somewhat to accommodate binding of the different inhibitors. The fact that NAPAP and MD-805 require different stereochemistry for potent inhibition is demonstrated, and its structural basis clarified. The wealth of data on analogs and variants of these lead compounds is shown to be compatible with this inhibitor binding mode.     

Virtual screening of naphthoquinone analogs for potent inhibitors against the cancer-signaling PI3K/AKT/mTOR pathway

The PI3K/AKT/mTOR pathway is one of the most commonly disrupted signaling pathways that plays a role in the development and pathogenicity of multiple cancers. Therefore, the critical proteins of this pathway have been targeted for anticancer therapy. The scientific community has increasingly been realizing the anti-cancer therapeutic potential of naphthoquinone analogs. These compounds constitute a major class of diverse sets of plant metabolites, which include various natural products and synthetic compounds with proven anticancer activity. The current study involved structural computational biology approaches to explore compounds from a diverse pool of naphthoquinone analogs that can inhibit key cancer-signaling proteins phosphoinositide 3-kinase (PI3K), protein kinase B, PKB (AKT), and mammalian target of rapamycin (mTOR). The novel compound identified commonly among the top 10 dock score lists of PI3K, AKT, and mTOR was selected for further study and proposed as a potential inhibitor of the 3 cancer-signaling proteins and an anticancer agent. Further, to check the docking accuracy and potential of the compound, post docking analyses, namely, binding comparison with the native ligand, the role of the interacting residue role in binding, predicted binding energy and dissociation constant calculations, etc., were performed. All these measures showed good-quality binding, and thus provide weight to our prediction of the novel compound as a pan PI3K/AKT/mTOR inhibitor and an anticancer agent. Finally, to compare the binding and similarity in the active sites of the 3 protein kinases, a ligand-based active site alignment was performed and analyzed. Thus, the study proposed a novel naphthoquinone analog as a potential anticancer drug, and provided comparative structural insight into its binding to the 3 protein kinases.     

Novel virtual lead identification in the discovery of hematopoietic cell kinase (HCK) inhibitors: application of 3D QSAR and molecular dynamics simulation

High level of hematopoietic cell kinase (Hck) is associated with drug resistance in chronic myeloid leukemia. Additionally, Hck activity has also been connected with the pathogenesis of HIV-1 and chronic obstructive pulmonary disease. In this study, three-dimensional (3D) QSAR pharmacophore models were generated for Hck based on experimentally known inhibitors. A best pharmacophore model, Hypo1, was developed with high correlation coefficient (0.975), Low RMS deviation (0.60) and large cost difference (49.31), containing three ring aromatic and one hydrophobic aliphatic feature. It was further validated by the test set (r?=?0.96) and Fisher’s randomization method (95%). Hypo 1 was used as a 3D query for screening the chemical databases, and the hits were further screened by applying Lipinski’s rule of five and ADMET properties. Selected hit compounds were subjected to molecular docking to identify binding conformations in the active site. Finally, the appropriate binding modes of final hit compounds were revealed by molecular dynamics (MD) simulations and free energy calculation studies. Hence, we propose the final three hit compounds as virtual candidates for Hck inhibitors.     

The hedgehog receptor patched is involved in cholesterol transport

Background

Sonic hedgehog (Shh) signaling plays a crucial role in growth and patterning during embryonic development, and also in stem cell maintenance and tissue regeneration in adults. Aberrant Shh pathway activation is involved in the development of many tumors, and one of the most affected Shh signaling steps found in these tumors is the regulation of the signaling receptor Smoothened by the Shh receptor Patched. In the present work, we investigated Patched activity and the mechanism by which Patched inhibits Smoothened.

Methodology/Principal Findings

Using the well-known Shh-responding cell line of mouse fibroblasts NIH 3T3, we first observed that enhancement of the intracellular cholesterol concentration induces Smoothened enrichment in the plasma membrane, which is a crucial step for the signaling activation. We found that binding of Shh protein to its receptor Patched, which involves Patched internalization, increases the intracellular concentration of cholesterol and decreases the efflux of a fluorescent cholesterol derivative (BODIPY-cholesterol) from these cells. Treatment of fibroblasts with cyclopamine, an antagonist of Shh signaling, inhibits Patched expression and reduces BODIPY-cholesterol efflux, while treatment with the Shh pathway agonist SAG enhances Patched protein expression and BODIPY-cholesterol efflux. We also show that over-expression of human Patched in the yeast S. cerevisiae results in a significant boost of BODIPY-cholesterol efflux. Furthermore, we demonstrate that purified Patched binds to cholesterol, and that the interaction of Shh with Patched inhibits the binding of Patched to cholesterol.

Conclusion/Significance

Our results suggest that Patched may contribute to cholesterol efflux from cells, and to modulation of the intracellular cholesterol concentration. This activity is likely responsible for the inhibition of the enrichment of Smoothened in the plasma membrane, which is an important step in Shh pathway activation.     

Enhanced potency of human Sonic hedgehog by hydrophobic modification

Post-translational modifications of the developmental signaling protein Sonic hedgehog (Shh) by a long-chain fatty acid at the N-terminus and cholesterol at the C-terminus greatly activate the protein in a cell-based signaling assay. To investigate the structural determinants of this activation phenomenon, hydrophobic and hydrophilic moieties have been introduced by chemical and mutagenic methods to the soluble N-terminal signaling domain of Shh and tested in both in vitro and in vivo assays. A wide variety of hydrophobic modifications increased the potency of Shh when added at the N-terminus of the protein, ranging from long-chain fatty acids to hydrophobic amino acids, with EC(50) values from 99 nM for the unmodified protein to 0.6 nM for the myristoylated form. The N-myristoylated Shh was as active as the natural form having both N- and C-terminal modifications. The degree of activation appears to correlate with the hydrophobicity of the modification rather than any specific chemical feature of the adduct; moreover, substitution with hydrophilic moieties decreased activity. Hydrophobic modifications at the C-terminus of Shh resulted in only a 2-3-fold increase in activity, and no activation was found with hydrophobic modification at other surface positions. The N-terminal modifications did not appear to alter the binding affinity of the Shh protein for the transfected receptor protein, Patched, and had no apparent effect on structure as measured by circular dichroism, thermal denaturation, and size determination. Activation of Desert Hh through modification of its N-terminus was also observed, suggesting that this is a common feature of Hh proteins.     

Sonic Hedgehog signaling impairs ionizing radiation-induced checkpoint activation and induces genomic instability

The Sonic Hedgehog (Shh) pathway plays important roles in embryogenesis, stem cell maintenance, tissue repair, and tumorigenesis. Haploinsufficiency of Patched-1, a gene that encodes a repressor of the Shh pathway, dysregulates the Shh pathway and increases genomic instability and the development of spontaneous and ionizing radiation (IR)–induced tumors by an unknown mechanism. Here we show that Ptc1^+/− mice have a defect in the IR-induced activation of the ATR–Chk1 checkpoint signaling pathway. Likewise, transient expression of Gli1, a downstream target of Shh signaling, disrupts Chk1 activation in human cells by preventing the interaction of Chk1 with Claspin, a Chk1 adaptor protein that is required for Chk1 activation. These results suggest that inappropriate Shh pathway activation promotes tumorigenesis by disabling a key signaling pathway that helps maintain genomic stability and inhibits tumorigenesis.     

Proliferation of Murine Midbrain Neural Stem Cells Depends upon an Endogenous Sonic Hedgehog (Shh) Source

The Sonic Hedgehog (Shh) pathway is responsible for critical patterning events early in development and for regulating the delicate balance between proliferation and differentiation in the developing and adult vertebrate brain. Currently, our knowledge of the potential role of Shh in regulating neural stem cells (NSC) is largely derived from analyses of the mammalian forebrain, but for dorsal midbrain development it is mostly unknown. For a detailed understanding of the role of Shh pathway for midbrain development in vivo, we took advantage of mouse embryos with cell autonomously activated Hedgehog (Hh) signaling in a conditional Patched 1 (Ptc1) mutant mouse model. This animal model shows an extensive embryonic tectal hypertrophy as a result of Hh pathway activation. In order to reveal the cellular and molecular origin of this in vivo phenotype, we established a novel culture system to evaluate neurospheres (nsps) viability, proliferation and differentiation. By recreating the three-dimensional (3-D) microenvironment we highlight the pivotal role of endogenous Shh in maintaining the stem cell potential of tectal radial glial cells (RGC) and progenitors by modulating their Ptc1 expression. We demonstrate that during late embryogenesis Shh enhances proliferation of NSC, whereas blockage of endogenous Shh signaling using cyclopamine, a potent Hh pathway inhibitor, produces the opposite effect. We propose that canonical Shh signaling plays a central role in the control of NSC behavior in the developing dorsal midbrain by acting as a niche factor by partially mediating the response of NSC to epidermal growth factor (EGF) and fibroblast growth factor (FGF) signaling. We conclude that endogenous Shh signaling is a critical mechanism regulating the proliferation of stem cell lineages in the embryonic dorsal tissue.     

Structure-based rational design of novel hit compounds for pyruvate dehydrogenase multienzyme complex E1 components from Escherichia coli

Pyruvate dehydrogenase multienzyme complex (PDHc) E1 component plays a pivotal role in cellular metabolism to convert the product of glycolysis (pyruvate) to acetyl-CoA, and has been reported as a potential target for anti-microbial and herbicide. In present study, based on the thiamin diphosphate (ThDP) site, four novel hit compounds with high inhibitory activity against the PDHc-E1 from Escherichia coli were firstly designed by using structure-based molecular docking methods. As expected, among four compounds, the compound 3a is the best inhibitor by far, with IC₅₀ value of 6.88 μM against PDHc-E1 from E. coli. To elucidate the interaction mechanism between the active site of PDHc-E1 and its inhibitor, the docking-based molecular dynamics simulation (MD) and MD-based ab initio fragment molecular orbital (FMO) calculations were also further performed. The positive results indicated that all modeling strategies presented in the current study most like to be an encouraging way in design of novel lead compounds with structural diversity for PDHc-E1 in the future.     

Chemical kinomics: a powerful strategy for target deconvolution

Kinomics is an emerging and promising approach for deciphering kinomes. Chemical kinomics is a discipline of chemical genomics that is also referred to as "chemogenomics", which is derived from chemistry and biology. Chemical kinomics has become a powerful approach to decipher complicated phosphorylation-based cellular signaling networks with the aid of small molecules that modulate kinase functions. Moreover, chemical kinomics has played a pivotal role in the field of kinase drug discovery as it enables identification of new molecular targets of small molecule kinase modulators and/or exploitation of novel functions of known kinases and has also provided novel chemical entities as hit/lead compounds. In this short review, contemporary chemical kinomics technologies such as activity-based protein profiling, T7 kinasetagged phages, kinobeads, three-hybrid systems, fluorescenttagged kinase binding assays, and chemical genomic profiling are discussed along with a novel allosteric Bcr-Abl kinase inhibitor (GNF-2/GNF-5) as a successful application of chemical kinomics approaches.     

An NMR and MD study of complexes of bacteriophage lambda lysozyme with tetra- and hexa-N-acetylchitohexaose

The X-ray structure of lysozyme from bacteriophage lambda (λ lysozyme) in complex with the inhibitor hexa-N-acetylchitohexaose (NAG6) (PDB: 3D3D) has been reported previously showing sugar units from two molecules of NAG6 bound in the active site. One NAG6 is bound with four sugar units in the ABCD sites and the other with two sugar units in the E′F′ sites potentially representing the cleavage reaction products; each NAG6 cross links two neighboring λ lysozyme molecules. Here we use NMR and MD simulations to study the interaction of λ lysozyme with the inhibitors NAG4 and NAG6 in solution. This allows us to study the interactions within the complex prior to cleavage of the polysaccharide. ¹H^N and ¹⁵N chemical shifts of λ lysozyme resonances were followed during NAG4/NAG6 titrations. The chemical shift changes were similar in the two titrations, consistent with sugars binding to the cleft between the upper and lower domains; the NMR data show no evidence for simultaneous binding of a NAG6 to two λ lysozyme molecules. Six 150 ns MD simulations of λ lysozyme in complex with NAG4 or NAG6 were performed starting from different conformations. The simulations with both NAG4 and NAG6 show stable binding of sugars across the D/E active site providing low energy models for the enzyme-inhibitor complexes. The MD simulations identify different binding subsites for the 5th and 6th sugars consistent with the NMR data. The structural information gained from the NMR experiments and MD simulations have been used to model the enzyme-peptidoglycan complex.