Analysis of a three-dimensional structure of human acidic mammalian chitinase obtained by homology modeling and ligand binding studies

The three-dimensional (3D) model of the human acidic mammalian chitinase (hAMCase) was constructed based on the crystal structure of the human chitotriosidase (EC 3.2.1.44, PDB code 1HKK) by using InsightII/Homology module. With the aid of molecular mechanics and molecular dynamics methods, the last refined model was obtained and further assessed by Profile-3D and Procheck, which confirms that the refined model is reliable. Furthermore, the docking results of the ligands (allosamidin and NAG₂) into the active site of hAMCase indicate that allosamidin is a more preferred ligand than NAG₂, and that Glu119 forms hydrogen bond with allosamidin, which is in good agreement with the experimental results. From the docking studies, we also suggest that Trp10, Glu49, Asp192, and Glu276 in hAMCase are four important determinant residues in binding as they have strong van-der-Waals and electrostatic interactions with the ligand, respectively.     

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.     

Homology model of the CDK1/cyclin B complex

We describe a refined homology model of a CDK1/cyclin B complex that was previously used for the structure-based optimization of the Paullone class of inhibitors. The preliminary model was formed from the homologous regions of the deposited CDK2/cyclin A crystal structure. Further refinement of the CDK1/cyclin B complex was accomplished using molecular mechanics and hydropathic analysis with a protocol of constraints and local geometry searches. For the most part, our CKD1/cyclin B homology model is very similar to the high resolution CDK2/cyclin A crystal structure regarding secondary and tertiary features. However, minor discrepancies between the two kinase structures suggest the possibility that ligand design may be specifically tuned for either CDK1 or CDK2. Our examination of the CDK1/cyclin B model includes a comparison with the CDK2/cyclin A crystal structure in the PSTAIRE interface region, connecting portions to the ATP binding domain, as well as the ATP binding site itself.     

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.     

Homology Model of the CDK1/cyclin B Complex

Abstract

We describe a refined homology model of a CDK1/cyclin B complex that was previously used for the structure-based optimization of the Paullone class of inhibitors. The preliminary model was formed from the homologous regions of the deposited CDK2/cyclin A crystal structure. Further refinement of the CDK1/cyclin B complex was accomplished using molecular mechanics and hydropathic analysis with a protocol of constraints and local geometry searches. For the most part, our CKD1/cyclin B homology model is very similar to the high resolution CDK2/cyclin A crystal structure regarding secondary and tertiary features. However, minor discrepancies between the two kinase structures suggest the possibility that ligand design may be specifically tuned for either CDK1 or CDK2. Our examination of the CDK1/cyclin B model includes a comparison with the CDK2/cyclin A crystal structure in the PSTAIRE interface region, connecting portions to the ATP binding domain, as well as the ATP binding site itself.     

Activation mechanism of CDK2: role of cyclin binding versus phosphorylation

Activation of the cyclin-dependent kinases is a two-step process involving cyclin binding followed by phosphorylation at a conserved threonine residue within the kinase activation loop. In this study, we describe the separate roles of cyclin A binding versus phosphorylation in the overall activation mechanism of CDK2. Interaction of CDK2 with cyclin A results in a partially active complex that is moderately defective in the binding of the protein substrate, but not ATP, and severely defective in both phosphoryl group transfer and turnover. Alternatively, phosphorylation of the CDK2 monomer also results in a partially activated species, but one that is severely (> or = 480-fold) defective in substrate binding exclusively. Catalytic turnover in the phosphorylated CDK2 monomer is largely unimpaired (approximately 8-fold lower). Our data support a model for the activation of CDK2 in vivo, in which interaction of unphosphorylated CDK2 with cyclin A serves to configure the active site for ground-state binding of both ATP and the protein substrate, and further aligns ATP in the transition state for phosphoryl transfer. Optimizing the alignment of protein substrates in the phosphoryl transfer reaction is the principal role of phosphorylation at Thr(160).     

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.     

Proteome-wide identification of cellular targets affected by bisindolylmaleimide-type protein kinase C inhibitors

Bisindolylmaleimide compounds such as GF109203X are potent inhibitors of protein kinase C (PKC) activity. Although bisindolylmaleimides are not entirely selective for PKC and are known to inhibit a few other protein kinases, these reagents have been extensively used to study the functional roles of PKC family enzymes in cellular signal transduction for more than a decade. Here, we establish a proteomics approach to gain further insights into the cellular effects of this compound class. Functional immobilization of suitable bisindolylmaleimide analogues in combination with the specific purification of cellular binding proteins by affinity chromatography led to the identification of several known and previously unknown enzyme targets. Subsequent in vitro binding and activity assays confirmed the protein kinases Ste20-related kinase and cyclin-dependent kinase 2 (CDK2) and the non-protein kinases adenosine kinase and quinone reductase type 2 as novel targets of bisindolylmaleimide inhibitors. As observed specifically for CDK2, minor chemical variation of the ligand by immobilizing the closely related bisindolylmaleimides III, VIII, and X dramatically affected target binding. These observed changes in affinity correlated with both the measured IC(50) values for in vitro CDK2 inhibition and results from molecular docking into the CDK2 crystal structure. Moreover, the conditions for affinity purification could be adapted in a way that immobilized bisindolylmaleimide III selectively interacted with either PKC alpha or ribosomal S6 protein kinase 1 only after activation of these kinases. Thus, we have established an efficient technique for the rapid identification of cellular bisindolylmaleimide targets and further demonstrate the comparative selectivity profiling of closely related kinase inhibitors within a cellular proteome.     

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.     

Discovery of cyclin-dependent kinase inhibitor, CR229, using structurebased drug screening

To generate new scaffold candidates as highly selective and potent cyclin-dependent kinase (CDK) inhibitors, structure-based drug screening was performed utilizing 3D pharmacophore conformations of known potent inhibitors. As a result, CR229 (6-bromo-2,3,4,9-tetrahydro-carbolin-1-one) was generated as the hit-compound. A computational docking study using the X-ray crystallographic structure of CDK2 in complex with CR229 was evaluated. This predicted binding mode study of CR229 with CDK2 demonstrated that CR229 interacted effectively with the Leu83 and Glu81 residues in the ATP-binding pocket of CDK2 for the possible hydrogen bond formation. Furthermore, biochemical studies on inhibitory effects of CR229 on various kinases in the human cervical cancer HeLa cells demonstrated that CR229 was a potent inhibitor of CDK2 (IC50: 3 microM), CDK1 (IC50: 4.9 microM), and CDK4 (IC50: 3 microM), yet had much less inhibitory effect (IC50: >20 microM) on other kinases, such as casein kinase 2-1 (CK2- alpha1), protein kinase A (PKA), and protein kinase C (PKC). Accordingly, these data demonstrate that CR229 is a potent CDK inhibitor with anticancer efficacy.     

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.     

A flavonoid gossypin binds to cyclin-dependent kinase 2

Flavonoids have low toxicity and mild activity. In order to find flavonoids showing cyclin-dependent kinase 2 (CDK2) binding effects, 347 flavonoid derivatives were docked into the crystal structure of the CDK2. The docking study showed that gossypin has a good conformational match with CDK2, which was confirmed by the binding affinity assay using NMR experiments.     

Targeting the cyclin-binding groove site to inhibit the catalytic activity of CDK2/cyclin A complex using p27KIP1-derived peptidomimetic inhibitors

Functionally activated cyclin-dependent kinase 2 (CDK2)/cyclin A complex has been validated as an interesting therapeutic target to develop the efficient antineoplastic drug based on the cell cycle arrest. Cyclin A binds to CDK2 and activates the kinases as well as recruits the substrate and inhibitors using a hydrophobic cyclin-binding groove (CBG). Blocking the cyclin substrate recruitment on CBG is an alternative approach to override the specificity hurdle of the currently available ATP site targeting CDK2 inhibitors. Greater understanding of the interaction of CDK2/cyclin A complex with p27 (negative regulator) reveals that the Leu-Phe-Gly (LFG) motif region of p27 binds with the CBG site of cyclin A to arrest the malignant cell proliferation that induces apoptosis. In the present study, Replacement with Partial Ligand Alternatives through Computational Enrichment (REPLACE) drug design strategies have been applied to acquire LFG peptide-derived peptidomimetics library. The peptidomimetics function is equivalent with respect to substrate p27 protein fashion but does not act as an ATP antagonist. The combined approach of molecular docking, molecular dynamics (MD), and molecular electrostatic potential and ADME/T prediction were carried out to evaluate the peptidomimetics. Resultant interaction and electrostatic potential maps suggested that smaller substituent is desirable at the position of phenyl ring to interact with Trp217, Arg250, and Gln254 residues in the active site. The best docked poses were refined by the MD simulations which resulted in conformational changes. After equilibration, the structure of the peptidomimetic and receptor complex was stable. The results revealed that the various substrate protein-derived peptidomimetics could serve as perfect leads against CDK2 protein.

Electronic supplementary material

The online version of this article (doi:10.1007/s12154-014-0124-y) contains supplementary material, which is available to authorized users.     

An insight into the inhibitory selectivity of 4-(Pyrazol- 4-yl)-pyrimidines to CDK4 over CDK2

Revealing selectivity mechanism of cyclin-dependent kinases (CDK) and their inhibitors is an important issue to develop potential anticancer drugs. The substituted 4-(Pyrazol-4-yl)-pyrimidines are potent inhibitors of CDK4 but not of the highly homologous CDK2. In order to reveal the inhibitory selectivity of these inhibitors to CDK4 over CDK2, we select one of substituted 4-(Pyrazol-4-yl)-pyrimidines as a representative (marked as A1 hereunder) and perform molecular docking, molecular dynamics simulations and binding free energy analysis for CDK4/A1 and CDK2/A1, respectively. The electrostatic and van der Waals (vdW) interactions of the A1 inhibitor with CDK4/CDK2 are discussed. The computed binding free energies based on the MM-PBSA method are consistent with experimental bioactivity ranking of A1 inhibitor to CDK4/CDK2. On the other hand, the conformational characteristics of CDK2 and CDK4 induced by A1 inhibitor are analysed and revealed. Results demonstrate that the vdW interactions considerably contribute to binding of CDK4/CDK2 with A1 inhibitor and are similar in size. The hydrogen bonding between A1 inhibitor and CDK4/CDK2 is considerably favourable to the binding, in which the hydrogen bond between the NH group of the pyrazole group of A1 and the residue Asp158 of CDK4 plays a crucial role in inhibitory selectivity of A1 inhibitor to CDK4 over CDK2. The electrostatic interaction energy differences between the corresponding residues of CDK4/A1 and CDK2/A1 confirm the above inference. The conformational changes of CDK2 and CDK4 induced by A1 inhibitor influence the selectivity of A1 inhibitor to CDK4/CDK2.     

Regulation of the cell cycle in response to inhibition of mitochondrial generated energy

Cell cycle control is regulated through the temporal action of both cyclin-dependent kinases and cyclin binding partners. Previously, we have demonstrated that low doses of oligomycin result in a cell cycle arrest of HL-60 cells in G(1) [S. Sweet, G. Singh, Accumulation of human promyelocytic leukemic (HL-60) cells at two energetic cell cycle checkpoints, Cancer Res. 55 (1995) 5164-5167]. In this study, we provide the molecular mechanisms for the observed G(1) arrest following mitochondrial ATPase inhibition. Protein expression of cyclin E and CDK2, the kinase activity of complexed cyclin E/CDK2, and protein expression of p16, p21, and p27 were all unaffected by oligomycin administration. While CDK4 levels were unchanged following oligomycin treatment, a dramatic reduction in cyclin D(1) was observed. Moreover, increased amounts of hypo-phosphorylated retinoblastoma protein (Rbp) and Rbp bound E2F were observed following mitochondrial ATP synthase inhibition. These data provide further evidence that surveillance of available energy occurs during G(1) and ATP deprivation results in cell cycle arrest via a reduction in cyclin D.     

Acetylation of cyclin-dependent kinase 5 is mediated by GCN5

Cyclin-dependent kinase 5 (CDK5), a member of atypical serine/threonine cyclin-dependent kinase family, plays a crucial role in pathophysiology of neurodegenerative disorders. Its kinase activity and substrate specificity are regulated by several independent pathways including binding with its activator, phosphorylation and S-nitrosylation. In the present study, we report that acetylation of CDK5 comprises an additional posttranslational modification within the cells. Among many candidates, we confirmed that its acetylation is enhanced by GCN5, a member of the GCN5-related N-acetyl-transferase family of histone acetyltransferase. Co-immunoprecipitation assay and fluorescent localization study indicated that GCN5 physically interacts with CDK5 and they are co-localized at the specific nuclear foci. Furthermore, liquid chromatography in conjunction with a mass spectrometry indicated that CDK5 is acetylated at Lys33 residue of ATP binding domain. Considering this lysine site is conserved among a wide range of species and other related cyclin-dependent kinases, therefore, we speculate that acetylation may alter the kinase activity of CDK5 via affecting efficacy of ATP coordination.     

Molecular insights of benzodipyrazole as CDK2 inhibitors: combined molecular docking,molecular dynamics,and 3D QSAR studies

Abstract

Benzodipyrazoles have been previously evaluated for their in vitro CDK2 inhibitory activity. In the current investigation, we identified a six-feature common pharmacophore model (AADDRR.33) which is predicted to be responsible for CDK2 inhibition. An efficient 3D QSAR (r^2?=?0.98 and q^2?=?0.82) model was also constructed by employing PLS regression analysis. From the molecular docking studies, we examined the binding patterns of compound 7aa with the target protein and also calculated the binding energy using MM-GBSA calculations. Three hydrogen bonds with Lys 33, Glu 81, and Leu 83 are conserved even after 1000?ps run in a molecular dynamics simulation. We identified the slight displacement in bond lengths and the conformational changes occurred during the dynamics. The results also elucidated the protein residue–ligand interaction fractions which clearly explained the involvement of non-H-bond interactions.     

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.