New compounds identified through in silico approaches reduce the α-synuclein expression by inhibiting prolyl oligopeptidase in vitro.

Prolyl oligopeptidase (POP) is a serine protease that is responsible for the maturation and degradation of short neuropeptides and peptide hormones. The inhibition of POP has been demonstrated in the treatment of α-synucleinopathies and several neurological conditions. Therefore, ligand-based and structure-based pharmacophore models were generated and validated in order to identify potent POP inhibitors. Pharmacophore-based and docking-based virtual screening of a drug-like database resulted in 20 compounds. The in vitro POP assays indicated that the top scoring compounds obtained from virtual screening, Hit 1 and Hit 2 inhibit POP activity at a wide range of concentrations from 0.1 to 10 µM. Moreover, treatment of the hit compounds significantly reduced the α-synuclein expression in SH-SY5Y human neuroblastoma cells, that is implicated in Parkinson's disease. Binding modes of Hit 1 and Hit 2 compounds were explored through molecular dynamics simulations. A detailed investigation of the binding interactions revealed that the hit compounds exhibited hydrogen bond interactions with important active site residues and greater electrostatic and hydrophobic interactions compared to those of the reference inhibitors. Finally, our findings indicated the potential of the identified compounds for the treatment of synucleinopathies and CNS related disorders.

Novel chemical scaffolds of the tumor marker AKR1B10 inhibitors discovered by 3D QSAR pharmacophore modeling.

AIM: Recent evidence suggests that aldo-keto reductase family 1 B10 (AKR1B10) may be a potential diagnostic or prognostic marker of human tumors, and that AKR1B10 inhibitors offer a promising choice for treatment of many types of human cancers. The aim of this study was to identify novel chemical scaffolds of AKR1B10 inhibitors using in silico approaches. METHODS: The 3D QSAR pharmacophore models were generated using HypoGen. A validated pharmacophore model was selected for virtual screening of 4 chemical databases. The best mapped compounds were assessed for their drug-like properties. The binding orientations of the resulting compounds were predicted by molecular docking. Density functional theory calculations were carried out using B3LYP. The stability of the protein-ligand complexes and the final binding modes of the hit compounds were analyzed using 10 ns molecular dynamics (MD) simulations. RESULTS: The best pharmacophore model (Hypo 1) showed the highest correlation coefficient (0.979), lowest total cost (102.89) and least RMSD value (0.59). Hypo 1 consisted of one hydrogen-bond acceptor, one hydrogen-bond donor, one ring aromatic and one hydrophobic feature. This model was validated by Fischer's randomization and 40 test set compounds. Virtual screening of chemical databases and the docking studies resulted in 30 representative compounds. Frontier orbital analysis confirmed that only 3 compounds had sufficiently low energy band gaps. MD simulations revealed the binding modes of the 3 hit compounds: all of them showed a large number of hydrogen bonds and hydrophobic interactions with the active site and specificity pocket residues of AKR1B10. CONCLUSION: Three compounds with new structural scaffolds have been identified, which have stronger binding affinities for AKR1B10 than known inhibitors.

Finding off-targets, biological pathways, and target diseases for chymase inhibitors via structure-based systems biology approach.

Off-target binding connotes the binding of a small molecule of therapeutic significance to a protein target in addition to the primary target for which it was proposed. Progressively such off-targeting is emerging to be regular practice to reveal side effects. Chymase is an enzyme of hydrolase class that catalyzes hydrolysis of peptide bonds. A link between heart failure and chymase is ascribed, and a chymase inhibitor is in clinical phase II for treatment of heart failure. However, the underlying mechanisms of the off-target effects of human chymase inhibitors are still unclear. Here, we develop a robust computational strategy that is applicable to any enzyme system and that allows the prediction of drug effects on biological processes. Putative off-targets for chymase inhibitors were identified through various structural and functional similarity analyses along with molecular docking studies. Finally, literature survey was performed to incorporate these off-targets into biological pathways and to establish links between pathways and particular adverse effects. Off-targets of chymase inhibitors are linked to various biological pathways such as classical and lectin pathways of complement system, intrinsic and extrinsic pathways of coagulation cascade, and fibrinolytic system. Tissue kallikreins, granzyme M, neutrophil elastase, and mesotrypsin are also identified as off-targets. These off-targets and their associated pathways are elucidated for the effects of inflammation, cancer, hemorrhage, thrombosis, and central nervous system diseases (Alzheimer's disease). Prospectively, our approach is helpful not only to better understand the mechanisms of chymase inhibitors but also for drug repurposing exercises to find novel uses for these inhibitors.

Characterization of modular bifunctional processive endoglucanase Cel5 from Hahella chejuensis KCTC 2396.

Cel5 from marine Hahella chejuensis is composed of glycoside hydrolase family-5 (GH5) catalytic domain (CD) and two carbohydrate binding modules (CBM6-2). The enzyme was expressed in Escherichia coli and purified to homogeneity. The optimum endoglucanase and xylanase activities of recombinant Cel5 were observed at 65 °C, pH 6.5 and 55 °C, pH 5.5, respectively. It exhibited K m of 1.8 and 7.1 mg/ml for carboxymethyl cellulose and birchwood xylan, respectively. The addition of Ca(2+) greatly improved thermostability and endoglucanase activity of Cel5. The Cel5 retained 90 % of its endoglucanase activity after 24 h incubation in presence of 5 M concentration of NaCl. Recombinant Cel5 showed production of cellobiose after hydrolysis of cellulosic substrates (soluble/insoluble) and methylglucuronic acid substituted xylooligosaccharides after hydrolysis of glucuronoxylans by endo-wise cleavage. These results indicated that Cel5 as bifunctional enzyme having both processive endoglucanase and xylanase activities. The multidomain structure of Cel5 is clearly distinguished from the GH5 bifunctional glycoside hydrolases characterized to date, which are single domain enzymes. Sequence analysis and homology modeling suggested presence of two conserved binding sites with different substrate specificities in CBM6-2 and a single catalytic site in CD. Residues Glu132 and Glu219 were identified as key catalytic amino acids by sequence alignment and further verified by using site directed mutagenesis. CBM6-2 plays vital role in catalytic activity and thermostability of Cel5. The bifunctional activities and multiple substrate specificities of Cel5 can be utilized for efficient hydrolysis of cellulose and hemicellulose into soluble sugars.

A combination of receptor-based pharmacophore modeling & QM techniques for identification of human chymase inhibitors.

Inhibition of chymase is likely to divulge therapeutic ways for the treatment of cardiovascular diseases, and fibrotic disorders. To find novel and potent chymase inhibitors and to provide a new idea for drug design, we used both ligand-based and structure-based methods to perform the virtual screening(VS) of commercially available databases. Different pharmacophore models generated from various crystal structures of enzyme may depict diverse inhibitor binding modes. Therefore, multiple pharmacophore-based approach is applied in this study. X-ray crystallographic data of chymase in complex with different inhibitors were used to generate four structure-based pharmacophore models. One ligand-based pharmacophore model was also developed from experimentally known inhibitors. After successful validation, all pharmacophore models were employed in database screening to retrieve hits with novel chemical scaffolds. Drug-like hit compounds were subjected to molecular docking using GOLD and AutoDock. Finally four structurally diverse compounds with high GOLD score and binding affinity for several crystal structures of chymase were selected as final hits. Identification of final hits by three different pharmacophore models necessitates the use of multiple pharmacophore-based approach in VS process. Quantum mechanical calculation is also conducted for analysis of electrostatic characteristics of compounds which illustrates their significant role in driving the inhibitor to adopt a suitable bioactive conformation oriented in the active site of enzyme. In general, this study is used as example to illustrate how multiple pharmacophore approach can be useful in identifying structurally diverse hits which may bind to all possible bioactive conformations available in the active site of enzyme. The strategy used in the current study could be appropriate to design drugs for other enzymes as well.

Development of predictive quantitative structure-activity relationship model and its application in the discovery of human leukotriene A4 hydrolase inhibitors.

BACKGROUND: Human LTA4H catalyzes the conversion of LTA4 to LTB4 and plays a key role in innate immune responses. Inhibition of this enzyme can be a valid method in the treatment of inflammatory response exhibited through LTB4. RESULTS & DISCUSSION: The quantitative structure-activity relationship (QSAR) models were developed using genetic function approximation and validated. A training set of 26 diverse compounds and their molecular descriptors were used to develop highly correlating QSAR models. A six-descriptor model explaining the biological activity of the training and test sets with correlation values of 0.846 and 0.502, respectively, was selected as the best model and used in a database screening of drug-like Maybridge database followed by molecular docking. CONCLUSION: Based on the predicted potent inhibitory activities, expected binding mode and molecular interactions at the active site of hLTA4H final leads were selected as to be utilized in designing future hLTA4H inhibitors.

Molecular modeling study on tunnel behavior in different histone deacetylase isoforms.

Histone deacetylases (HDACs) have emerged as effective therapeutic targets in the treatment of various diseases including cancers as these enzymes directly involved in the epigenetic regulation of genes. However the development of isoform-selective HDAC inhibitors has been a challenge till date since all HDAC enzymes possess conserved tunnel-like active site. In this study, using molecular dynamics simulation we have analyzed the behavior of tunnels present in HDAC8, 10, and 11 enzymes of class I, II, and IV, respectively. We have identified the equivalent tunnel forming amino acids in these three isoforms and found that they are very much conserved with subtle differences to be utilized in selective inhibitor development. One amino acid, methionine of HDAC8, among six tunnel forming residues is different in isoforms of other classes (glutamic acid (E) in HDAC10 and leucine (L) in HDAC 11) based on which mutations were introduced in HDAC11, the less studied HDAC isoform, to observe the effects of this change. The HDAC8-like (L268M) mutation in the tunnel forming residues has almost maintained the deep and narrow tunnel as present in HDAC8 whereas HDAC10-like (L268E) mutation has changed the tunnel wider and shallow as observed in HDAC10. These results explained the importance of the single change in the tunnel formation in different isoforms. The observations from this study can be utilized in the development of isoform-selective HDAC inhibitors.