Ligand-Based Molecular Modeling Study on a Chemically Diverse Series of Cholecystokinin-B/Gastrin Receptor Antagonists: Generation of Predictive Model

Pharmacophore-based virtual screening of ZINC database, molecular modeling and designing new derivatives as potential HDAC6 inhibitors

To date, many HDAC6 inhibitors have been identified and developed but none is clinically approved as of now. Through this study, we aim to obtain novel HDAC6 selective inhibitors and provide new insights into the detailed structural design of potential HDAC6 inhibitors. A HypoGen-based 3D QSAR HDAC6 pharmacophore was built and used as a query model to screen approximately 8 million ZINC database compounds. First, the ZINC Database was filtered using ADMET, followed by pharmacophore-based library screening. Using fit value and estimated activity cutoffs, a final set of 54 ZINC hits was obtained that were further investigated using molecular docking with the crystal structure of human histone deacetylase 6 catalytic domain 2 in complex with Trichostatin A (PDB ID: 5EDU). Through detailed in silico screening of the ZINC database, we shortlisted three hits as the lead molecules for designing novel HDAC6 inhibitors with better efficacy. Docking with 5EDU, followed by ADMET and TOPKAT analysis of modified ZINC hits provided 9 novel potential HDAC6 inhibitors that possess better docking scores and 2D interactions as compared to the control ZINC hit molecules. Finally, a 50 ns MD analysis run followed by Protein-Ligand Interaction Energy (PLIE) analysis of the top scored hits provided a novel molecule N1 that showed promisingly similar results to that of Ricolinostat (a known HDAC6 inhibitor). The comparable result of the designed hits to established HDAC6 inhibitors suggests that these compounds might prove to be successful HDAC6 inhibitors in future. Designed novel hits that might act as good HDAC6 inhibitors derived from ZINC database using combined molecular docking and modeling approaches.

CCK2R antagonists: from SAR to clinical trials

The widespread involvement of the cholecystokinin-2/gastrin receptor (CCK₂R) in multiple (patho)physiological processes has propelled extensive searches for nonpeptide small-molecule CCK₂R antagonists. For the past three decades, considerable research has yielded numerous chemically heterogeneous compounds. None of these entered into the clinic, mainly because of inadequate biological effects. However, it appears that the ultimate goal of a clinically useful CCK₂R antagonist is now just around the corner, with the most promising compounds, netazepide and nastorazepide, now in Phase II clinical trials. Here, we illustrate the structure-activity relationships (SARs) of stablished CCK₂R antagonists of various structural classes, and the most recent proof-of-concept studies where new applicabilities of CCK₂R antagonists as visualizing agents are presented.

Identification of potent cholecystokinin-B receptor antagonists: synthesis, molecular modeling and anti-cancer activity against pancreatic cancer cells

Advanced malignant stages of pancreatic cancer have poor prognosis and very few treatment strategies are available. Pancreatic cancer is known to possess unique growth-related receptors that when activated, stimulate tumour proliferation. Gastrin and its related peptide cholecystokinin (CCK) are also significantly involved in the growth of this cancer type as well as other malignancies through activation of the cholecystokinin-B receptor (CCK-BR). New treatment strategies with CCK-BR antagonists are being suggested that suppress the growth promoting effects of gastrin. In this paper, we report the development of two series of quinazolinone derivatives incorporating hydrazinecarbothioamide (compounds 3a-g) and the hydrazino group (compounds 4a-e) as linkers for developing CCK-BR antagonists. The affinities of the compounds were determined using docking into the CCK-BR homology modeled structure. The compounds were tested for in vitro CCK-BR binding and gastric acid secretion in an isolated lumen-perfused mouse stomach assay. The compounds exhibited CCK-BR binding activity (IC₅₀) in the range of 0.2-975 nM and showed good gastric acid secretion inhibitory activity. Molecular modeling of the compounds was done and pharmacophore mapping results showed good prediction of in vitro activity which correlated well with the experimental antagonistic activity. The compounds were further tested for their cytotoxicity on CCK-BR expressing pancreatic cancer cells. The results of the study provided two potent CCK-BR antagonists which also possess good to moderate growth inhibitory activities against pancreatic cancer cells.

Identification of Structurally Diverse Antimicrobials Through Sequential Application of Pharmacophore Modeling, Virtual Screening, Molecular Docking and In Vitro Microbiological Assay

Dihydrofolate reductase enzyme has been an attractive biological target for the design and development of antimicrobials. Considering this, we have attempted to identify novel dihydrofolate reductase inhibitors through our well-defined in silico and in vitro work flow. An accurate and predictive pharmacophore model comprising of one hydrogen bond acceptor, two hydrophobic and one ring aromatic was developed and utilized as a query to search the National Cancer Institute and Maybridge database leading to retrieval of various compounds which were filtered on the basis of estimated activity, fit value and Lipinski's violation. Selected hits NSC3423, KM09759, NSC391, NSC2091 and HTS00630 were subjected to docking studies which resulted into visualization of potential interaction capabilities of hits in line to pharmacophoric features. The identified hits were evaluated for in vitro antimicrobial potential, and the results revealed that among all the five hits, NSC3423 is the most potent compound with activity against E. coli, P. aeruginosa, S. aureus, B. substilis, A. niger and F. oxysporum. On the other hand, KM09759, NSC391, NSC2091 and HTS00630 showed varying degree of activities against gram-positive, gram-negative and fungal strains.

Identification of novel HIV 1--protease inhibitors: application of ligand and structure based pharmacophore mapping and virtual screening

A combined ligand and structure-based drug design approach provides a synergistic advantage over either methods performed individually. Present work bestows a good assembly of ligand and structure-based pharmacophore generation concept. Ligand-oriented study was accomplished by employing the HypoGen module of Catalyst in which we have translated the experimental findings into 3-D pharmacophore models by identifying key features (four point pharmacophore) necessary for interaction of the inhibitors with the active site of HIV-1 protease enzyme using a training set of 33 compounds belonging to the cyclic cyanoguanidines and cyclic urea derivatives. The most predictive pharmacophore model (hypothesis 1), consisting of four features, namely, two hydrogen bond acceptors and two hydrophobic, showed a correlation (r) of 0.90 and a root mean square of 0.71 and cost difference of 56.59 bits between null cost and fixed cost. The model was validated using CatScramble technique, internal and external test set prediction. In the second phase of our study, a structure-based five feature pharmacophore hypothesis was generated which signifies the importance of hydrogen bond donor, hydrogen bond acceptors and hydrophobic interaction between the HIV-1 protease enzyme and its inhibitors. This work has taken a significant step towards the full integration of ligand and structure-based drug design methodologies as pharmacophoric features retrieved from structure-based strategy complemented the features from ligand-based study hence proving the accuracy of the developed models. The ligand-based pharmacophore model was used in virtual screening of Maybridge and NCI compound database resulting in the identification of four structurally diverse druggable compounds with nM activities.

Toward the identification of a reliable 3D QSAR pharmacophore model for the CCK2 receptor antagonism

The present study describes application of computational approaches to identify a validated and reliable 3D QSAR pharmacophore model for the CCK-2R antagonism through integrated ligand and structure based studies using anthranilic sulfonamide and 1,3,4-benzotriazepine based CCK-2R antagonists. The best hypothesis consisted five features viz. two aliphatic hydrophobic, one aromatic hydrophobic, one H-bond acceptor, and one ring aromatic feature with an excellent correlation for 34 training set (r²(training) = 0.83) and 58 test set compounds (r²(test) = 0.74). This model was validated through F-test and docking studies at the active site of the plausible CCK-2R where the 99% significance and well corroboration with the pharmacophore model respectively describes the model's reliability. The model also predicts well to other known clinically effective CCK-2R antagonists. Therefore, the developed model may useful in finding new scaffolds that may aid in design and develop new chemical entities (NCEs) as potent CCK-2R antagonists before their synthesis.

Targeting the fatty acid biosynthesis enzyme, beta-ketoacyl-acyl carrier protein synthase III (PfKASIII), in the identification of novel antimalarial agents

The importance of fatty acids to the human malaria parasite, Plasmodium falciparum, and differences due to a type I fatty acid synthesis (FAS) pathway in the parasite, make it an attractive drug target. In the present study, we developed and a utilized a pharmacophore to select compounds for testing against PfKASIII, the initiating enzyme of FAS. This effort identified several PfKASIII inhibitors that grouped into various chemical classes of sulfides, sulfonamides, and sulfonyls. Approximately 60% of the submicromolar inhibitors of PfKASIII inhibited in vitro growth of the malaria parasite. These compounds inhibited both drug sensitive and resistant parasites and testing against a mammalian cell line revealed an encouraging in vitro therapeutic index for the most active compounds. Docking studies into the active site of PfKASIII suggest a potential binding mode that exploits amino acid residues at the mouth of the substrate tunnel.

Rationalizing the activities of diverse cholecystokinin 2 receptor antagonists using molecular field points

Cholecystokinin 2 receptor antagonists encompass a wide range of structures. This makes them unsuitable candidates for existing 3D-QSAR methods and has led us to develop an alternative approach to account for their observed biological activities. A diverse set of 21 antagonists was subjected to a novel molecular field-based similarity analysis. The hypothesis is that compounds with similar field patterns will bind at the same target site regardless of their underlying structure. This initial report demonstrates a linear correlation between ligand similarity and biological activity for this challenging data set. A model generated with three molecules was used to predict the activity of 18 test compounds, with different chemotypes, with a root-mean-square error of 0.68 pKB units. The ability to automatically derive a molecular alignment without knowledge of the protein structure represents an improvement over existing pharmacophore methods and makes the method particularly suitable for scaffold-hopping.

Partial Agonism, Neutral Antagonism, and Inverse Agonism at the Human Wild-Type and Constitutively Active Cholecystokinin-2 Receptors

Cholecystokinin receptor-2 (CCK2R) is a G protein receptor that regulates a number of physiological functions. Activation of CCK2R and/or expression of a constitutively active CCK2R variant may contribute to human diseases, including digestive cancers. Search for antagonists of the CCK2R has been an important challenge during the last few years, leading to discovery of a set of chemically distinct compounds. However, several early-discovered antagonists turned out to be partial agonists. In this context, we carried out pharmacological characterization of six CCK2R antagonists using COS-7 cells expressing the human CCK2R or a CCK2R mutant having a robust constitutive activity on inositol phosphates production, and we investigated the molecular mechanisms which, at a CCK2R binding site, account for these features. Results indicated that three compounds, 3R(+)-N-(2,3-dihydro-1-methyl-2-oxo-5-phenyl-1H-1,4-benzodiazepin-3-yl)-N'-(3-methylphenyl)urea (L365,260), 4-{[2-[[3-(lH-indol-3-yl)-2-methyl-1-oxo-2-[[[1.7.7-trimethyl-bicyclo[2.2.1]hept-2-yl)-oxy]carbonyl]amino]propyl]amino]-1-phenylethyl]amino-4-oxo-[lS-la.2[S*(S*)]4a]}-butanoate N-methyl-D-glucamine (PD135,158), and (R)-1-naphthalenepropanoic acid, b-[2-[[2-(8-azaspiro-[4.5]dec-8-ylcarbonyl)-4,6-dimethylphenyl]amino]-2-oxoethyl] (CR2945), were partial agonists; one molecule, 1-[(R)-2,3-dihydro-1-(2,3-dihydro-1-(2-methylphenacyl)-2-oxo-5-phenyl-1H-1,4-benzodiazepin-3-yl]-3-(3-methylphenyl)urea (YM022), was a neutral antagonist; and two compounds, N-(+)-[1-(adamant-1-ylmethyl)-2,4-dioxo-5-phenyl2,3,4,5-tetrahydro-1H-1,5-benzodiazepin-3-yl]-N'-phenylurea (GV150,013X) and ([(N-[methoxy-3 phenyl] N-[N-methyl N-phenyl carbamoylmethyl], carbomoyl-methyl)-3 ureido]-3-phenyl)2-propionic acid (RPR101,048), were inverse agonists. Furthermore, target- and pharmacophore-based docking of ligands followed by molecular dynamic simulation experiments resulted in consistent motion of aromatic residues belonging to a network presumably important for activation, thus providing the first structural explanations for the different pharmacological profiles of tested compounds. This study confirms that several referenced so-called antagonists are in fact partial agonists, and because of this undesired activity, we suggest that newly generated molecules should be preferred to efficiently block CCK2R-related physiological effects. Furthermore, data on the structural basis for the different pharmacological features of CCK2R ligands will serve to further clarify CCK2R mechanism of activation.