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.     

Rational design,synthesis, and verification of affinity ligands to a protein surface cleft

The structure-based design, synthesis, and screening of a glucuronic acid scaffold library of affinity ligands directed toward the catalytic cleft on porcine pancreas alpha-amylase are presented. The design was based on the simulated docking to the enzyme active site of 53 aryl glycosides from the Available Chemicals Directory (ACD) selected by in silico screening. Twenty-three compounds were selected for synthesis and screened in solution for binding toward alpha-amylase using nuclear magnetic resonance techniques. The designed molecules include a handle outside of the binding site to allow their attachment to various surfaces with minimal loss of binding activity. After initial screening in solution, one affinity ligand was selected, immobilized to Sepharose (Amersham Biosciences), and evaluated as a chromatographic probe. A column packed with ligand-coupled Sepharose specifically retained the enzyme, which could be eluted by a known inhibitor.     

In silico multi-filter screening approaches for developing novel beta-secretase inhibitors

A large database of chemical structures was screened for potential inhibitors of β-secretase was carried out using in silico multi-filter techniques. Substructure screening, computer-aided ligand docking, binding free energy calculations, and partial interaction energy analyses were performed successively to identify chemical compounds which could serve as different scaffolds from known β-secretase inhibitors for future drug design. We showed that our in silico multi-filter screening retrieved all known inhibitors from the compound database investigated, which suggests that the other compounds identified as inhibitors by this computerized screening process are potential β-secretase inhibitors.     

Identification of novel inhibitors of HCV RNA-dependent RNA polymerase by pharmacophore-based virtual screening and in vitro evaluation

Hepatitis C virus (HCV) is the major etiological agent of non-A, non-B hepatitis where no effective treatment is available. The HCV NS5B with RNA-dependent RNA polymerase (RdRp) activity is a key target for the treatment of HCV infection. Here we report novel NS5B polymerase inhibitors identified by virtual screening and in vitro evaluation of their inhibitory activities. On the basis of a newly identified binding pocket of NS5B, distinct from the nucleotide binding site but highly conserved among various HCV isolates, we performed virtual screening of compounds that fit this binding pocket from the available chemical database of 3.5 million compounds. The inhibitory activities of the in silico selected 119 compounds were estimated with in vitro RdRp assay. Three compounds with IC50 values of about 20 μM were identified, and their kinetic analyses suggest that these compounds are noncompetitive inhibitors with respect to the ribonucleotide substrate. Furthermore, the single-point mutations of the conserved residues in the binding pocket of NS5B resulted in the significant decrease of the RdRp activity, indicating that the binding pocket presented here might be important for the therapeutic intervention of HCV. These novel inhibitors would be useful for the development of effective anti-HCV agents.     

Discovery of potential pancreatic cholesterol esterase inhibitors using pharmacophore modelling, virtual screening, and optimization studies

Pancreatic cholesterol esterase (CEase) is a serine hydrolase involved in the hydrolysis of variety of lipids and transport of free cholesterol. In this study, pharmacophore hypotheses based on known inhibitors were generated using common feature pharmacophore generation protocol available in Discovery Studio program. The best pharmacophore model containing two hydrogen bond acceptor and three hydrophobic features was selected and validated. It was further used in screening three diverse chemical databases. Hit compounds were subjected to drug-likeness and molecular docking studies. Four hits, namely SEW00846, NCI0040784, GK03167, and CD10645, were selected based on the GOLD fitness score and interaction with active site amino acids. All hit compounds were further optimized to improve their binding in the active site. The optimized compounds were found to have improved binding at the active site. Strongly binding optimized hits at the active site can act as virtual leads in potent CEase inhibitor designing.     

Homology modeling,molecular docking and spectra assay studies of sterol 14α-demethylase from <Emphasis Type="Italic">Penicillium digitatum</Emphasis>

Sterol 14α-demethylase from Penicillium digitatum (PdCYP51) is a prime target of antifungal drugs for citrus disease in plants. To design novel antifungal compounds, a homology model of PdCYP51 was constructed using the recently reported crystal structure of human CYP51 as the template. Molecular docking was performed to investigate the interaction of four commercial fungicides with the modeled enzyme. The side chain of these compounds interplayed with PdCYP51 mainly through hydrophobic and van der Waals interactions. Biochemical spectra analysis of inhibitors combined with PdCYP51 are also compatible with the docking results. This is the first molecular modeling for PdCYP51 based on the eukaryotic crystal structure of CYP51. The structural information and binding site mapping of PdCYP51 for different inhibitors obtained from this study could aid in screening and designing new antifungal compounds targeting this enzyme.     

In Silico Design of Dual-Binding Site Anti-Cholinesterase Phytochemical Heterodimers as Treatment Options for Alzheimer’s Disease

The number of patients with neurodegenerative diseases, particularly Alzheimer’s disease (AD), continues to grow yearly. Cholinesterase inhibitors (ChEIs) represent the first-line symptomatic drug treatment for mild-to-moderate AD; however, there is an unmet need to produce ChEIs with improved efficacy and reduced side effects. Herein, phytochemicals with reported anti-acetylcholinesterase (AChE) activity were ranked in silico for their anti-AChE potential. Ligands with a similar or higher binding affinity to AChE than galantamine were then selected for the design of novel dual-binding site heterodimeric drugs. In silico molecular docking of heterodimers with the target enzymes, AChE and butyrylcholinesterase (BuChE), were performed, and anti-cholinesterase binding affinities were compared with donepezil. Drug-likeliness properties and toxicity of the heterodimers were assessed using the SwissADME and ProTox-II webservers. Nine phytochemicals displayed similar or higher binding affinities to AChE than galantamine: sanguinarine > huperzine A > chelerythrine > yohimbine > berberine > berberastine > naringenin > akuammicine > carvone. Eleven heterodimeric ligands were designed with phytochemicals separated by four- or five-carbon alkyl-linkers. All heterodimers were theoretically potent AChE and BuChE dual-binding site inhibitors, with the highest affinity achieved with huperzine-4C-naringenin, which displayed 34% and 26% improved affinity to AChE and BuChE, respectively, then the potent ChEI drug, donepezil. Computational pharmacokinetic and pharmacodynamic screening suggested that phytochemical heterodimers would display useful gastrointestinal absorption and with relatively low predicted toxicity. Collectively, the present study suggests that phytochemicals could be garnered for the provision of novel ChEIs with enhanced drug efficacy and low toxicity.     

Distinct effects of two HIV-1 capsid assembly inhibitor families that bind the same site within the N-terminal domain of the viral CA protein

The emergence of resistance to existing classes of antiretroviral drugs necessitates finding new HIV-1 targets for drug discovery. The viral capsid (CA) protein represents one such potential new target. CA is sufficient to form mature HIV-1 capsids in vitro, and extensive structure-function and mutational analyses of CA have shown that the proper assembly, morphology, and stability of the mature capsid core are essential for the infectivity of HIV-1 virions. Here we describe the development of an in vitro capsid assembly assay based on the association of CA-NC subunits on immobilized oligonucleotides. This assay was used to screen a compound library, yielding several different families of compounds that inhibited capsid assembly. Optimization of two chemical series, termed the benzodiazepines (BD) and the benzimidazoles (BM), resulted in compounds with potent antiviral activity against wild-type and drug-resistant HIV-1. Nuclear magnetic resonance (NMR) spectroscopic and X-ray crystallographic analyses showed that both series of inhibitors bound to the N-terminal domain of CA. These inhibitors induce the formation of a pocket that overlaps with the binding site for the previously reported CAP inhibitors but is expanded significantly by these new, more potent CA inhibitors. Virus release and electron microscopic (EM) studies showed that the BD compounds prevented virion release, whereas the BM compounds inhibited the formation of the mature capsid. Passage of virus in the presence of the inhibitors selected for resistance mutations that mapped to highly conserved residues surrounding the inhibitor binding pocket, but also to the C-terminal domain of CA. The resistance mutations selected by the two series differed, consistent with differences in their interactions within the pocket, and most also impaired virus replicative capacity. Resistance mutations had two modes of action, either directly impacting inhibitor binding affinity or apparently increasing the overall stability of the viral capsid without affecting inhibitor binding. These studies demonstrate that CA is a viable antiviral target and demonstrate that inhibitors that bind within the same site on CA can have distinct binding modes and mechanisms of action.     

Identification of Putative Plant-Based ALR-2 Inhibitors to Treat Diabetic Peripheral Neuropathy

Diabetic peripheral neuropathy (DPN) is a common diabetes complication (DM). Aldose reductase -2 (ALR-2) is an oxidoreductase enzyme that is most extensively studied therapeutic target for diabetes-related complications that can be inhibited by epalrestat, which has severe adverse effects; hence the discovery of potent natural inhibitors is desired. In response, a pharmacophore model based on the properties of eplarestat was generated. The specified pharmacophore model searched the NuBBE_DB database of natural compounds for prospective lead candidates. To assess the drug-likeness and ADMET profile of the compounds, a series of in silico filtering procedures were applied. The compounds were then put through molecular docking and interaction analysis. In comparison to the reference drug, four compounds showed increased binding affinity and demonstrated critical residue interactions with greater stability and specificity. As a result, we have identified four potent inhibitors: ZINC000002895847, ZINC000002566593, ZINC000012447255, and ZINC000065074786, that could be used as pharmacological niches to develop novel ALR-2 inhibitors.     

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.     

Identification of small-molecule inhibitors of the steroidogenic acute regulatory protein (STARD1) by structure-based design

A homology model of the steroidogenic acute regulatory protein (STAR)-related lipid transfer (START) domain of STARD1 was built, and the cholesterol binding site was identified. Structure-based design studies were performed to identify small molecule inhibitors of the START domain. The lead compounds were selected based on cAMP-induced, but not 22R-hydroxycholesterol-supported, inhibition of steroid synthesis by 50% at 10 μM. The results obtained by molecular docking & dynamics show a good correlation between bioactivity, docking scores and calculated binding energies of ligand-protein complexes. The best active compounds will be optimized further and used to develop potential drugs to control excessive steroid formation.     

Pharmacophore modeling,resistant mutant isolation,docking, and MM-PBSA analysis: Combined experimental/computer-assisted approaches to identify new inhibitors of the bovine viral diarrhea virus (BVDV)

Starting from a series of our new 2-phenylbenzimidazole derivatives, shown to be selectively and potently active against the bovine viral diarrhea virus (BVDV), we developed a hierarchical combined experimental/molecular modeling strategy to explore the drug leads for the BVDV RNA-dependent RNA-polymerase. Accordingly, a successful 3D pharmacophore model was developed, characterized by distinct chemical features that may be responsible for the activity of the inhibitors. BVDV mutants resistant to lead compounds in our series were then isolated, and the mutant residues on the viral molecular target, the RNA-dependent RNA-polymerase, were identified. Docking procedures upon pharmacophoric constraints and mutational data were carried out, and the binding affinity of all active compounds for the RdRp were estimated. Given the excellent agreement between in silico and in vitro data, this procedure is currently being employed in the design a new series of more selective and potent BVDV inhibitors.     

Characterization of isoflavonoids as inhibitors of β-hydroxyacyl-acyl carrier protein dehydratase (FabZ) from Moraxella catarrhalis: Kinetics,spectroscopic, thermodynamics and in silico studies

Backgroud

β-hydroxyacyl-acyl carrier protein dehydratase (FabZ) is an essential component of type II fatty acid biosynthesis (FAS II) pathway in bacteria. It performs dehydration of β-hydroxyacyl-ACP to trans-2-acyl-ACP in the elongation cycle of the FAS II pathway. FabZ is ubiquitously expressed and has uniform distribution, which makes FabZ an excellent target for developing novel drugs against pathogenic bacteria.

Structure modeling, ligand binding, and binding affinity calculation (LR-MM-PBSA) of human heparanase for inhibition and drug design

Human heparanase is an endo-beta-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 (alpha/beta)(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.     

Nanomolar binding affinity of quinine-based antimalarial compounds by the cocaine-binding aptamer

An unusual feature of the cocaine-binding aptamer is that it binds quinine much tighter than the ligand it was selected for, cocaine. Here we expand the repertoire of ligands that this aptamer binds to include the quinine-based antimalarial compounds amodiaquine, mefloquine, chloroquine and primaquine. Using isothermal titration calorimetry (ITC) we show that amodiaquine is bound by the cocaine-binding aptamer with an affinity of (7?±?4) nM, one of the tightest aptamer-small molecule affinities currently known. Amodiaquine, mefloquine and chloroquine binding are driven by both a favorable entropy and enthalpy of binding, while primaquine, quinine and cocaine binding are enthalpy driven with unfavorable binding entropy. Using nuclear magnetic resonance (NMR) and ITC methods we show that these ligands compete for the same binding sites in the aptamer. Our identification of such a tight binding ligand for this aptamer should prove useful in developing new biosensor techniques and applications using the cocaine-binding aptamer as a model system.     

Serotonin uptake inhibitors differentially modulate high affinity imipramine dissociation in human platelet membranes

The influence of selective serotonin uptake inhibitors on the dissociation rate of 3H-imipramine from its high affinity binding site in human platelet membranes was studied. Of the uptake inhibitors tested, one group of compounds attenuated dissociation, another group accelerated it, while a third group had little effect on the dissociation process. Drugs unrelated to the serotonin uptake system were ineffective. Removal of sodium ions markedly increased imipramine dissociation. Dose response curves of the active compounds indicated that micromolar concentrations were required to exert an effect on imipramine dissociation. These results can be adequately explained by an allosteric model which includes effector binding sites and distinct conformational states of the high affinity imipramine binding site.     

Discovery of a novel family of CDK inhibitors with the program LIDAEUS: structural basis for ligand-induced disordering of the activation loop

A family of 4-heteroaryl-2-amino-pyrimidine CDK2 inhibitor lead compounds was discovered with the new database-mining program LIDAEUS through in silico screening. Four compounds with IC(50) values ranging from 17 to 0.9 microM were selected for X-ray crystal analysis. Two distinct binding modes are observed, one of which resembles the hydrogen bonding pattern of bound ATP. In the second binding mode, the ligands trigger a conformational change in the activation T loop by inducing movement of Lys(33) and Asp(145) side chains. The family of molecules discovered provides an excellent starting point for the design and synthesis of tight binding inhibitors, which may lead to a new class of antiproliferative drugs.     

Drug screening of α-amylase inhibitors as candidates for treating diabetes

In the present study, the identification of potential α-amylase inhibitors is explored as a potential strategy for treating type-2 diabetes mellitus. A computationally driven approach using molecular docking was employed to search for new α-amylase inhibitors. The interactions of potential drugs with the enzyme's active site were investigated and compared with the contacts established by acarbose (a reference drug for α-amylase inhibition) in the crystallographic structure 1B2Y. For this active site characterization, both molecular docking and molecular dynamics simulations were performed, and the residues involved in the α-amylase–acarbose complex were considered to analyse the potential drug's interaction with the enzyme. Two potential α-amylase inhibitors (AN-153I105594 and AN-153I104845) have been selected following this computational strategy. Both compounds established a large number of interactions with key binding site α-amylase amino acids and obtained a comparable docking score concerning the reference drug (acarbose). Aiming to further analyse candidates' properties, their ADME (absorption, distribution, metabolism, excretion) parameters, druglikeness, organ toxicity, toxicological endpoints and median lethal dose (LD₅₀) were estimated. Overall estimations are promising for both candidates, and in silico toxicity predictions suggest that a low toxicity should be expected.     

Towards more effective acetylcholinesterase inhibitors: a comprehensive modelling study based on human acetylcholinesterase protein–drug complex

Abstract

Alzheimer's disease is a progressive neurodegenerative disorder and as the exact cause of the disease remains unknown, it still has no cure to date. Due to the fact that, until recently, there has been no crystal structure of human AChE in complex with drugs, researchers have had to use mainly Torpedo californica homologues which were later reported to have significantly different binding sites. In this study, an energy-based pharmacophore model that has the advantages of both ligand- and structure-based approaches was generated using hAChE crystal structures in complex with drugs. This model was validated utilizing several commonly used statistical measures such as the enrichment factor, BEDROC, RIE and AUAC. A huge database consisting of around 7 million compounds with approximately 150,000,000 conformations from 8 vendors was used in virtually screening workflow, in which Lipinski’s filter and basic and precise docking algorithms were utilized. A rigorous MM-GBSA binding affinity calculation was also applied to accurately predict relative free energy that produced 361 hits. The selected top ranked 15 compounds were shown to have extra intermolecular interactions with hAChE which is an indication of a more stable complex and high binding affinity. The e-pharmacophore model and overall results obtained might be used for further experimental studies in designing the next generation of hAChE inhibitors.

Communicated by Ramaswamy H. Sarma     

Profiling the interaction of 1‐phenylbenzimidazoles to cyclooxygenases

In the design of 1‐phenylbenzimidazoles as model cyclooxygenase (COX) inhibitors, docking to a series of crystallographic COX structures was performed to evaluate their potential for high‐affinity binding and to reproduce the interaction profile of well‐known COX inhibitors. The effect of ligand‐specific induced fit on the calculations was also studied. To quantitatively compare the pattern of interactions of model compounds to the profile of several cocrystallized COX inhibitors, a geometric parameter, denominated ligand‐receptor contact distance (LRCD), was developed. The interaction profile of several model complexes showed similarity to the profile of COX complexes with inhibitors such as iodosuprofen, iodoindomethacin, diclofenac, and flurbiprofen. Shaping of high‐affinity binding sites upon ligand‐specific induced fit mostly determined both the affinity and the binding mode of the ligands in the docking calculations. The results suggest potential of 1‐phenylbenzimidazole derivatives as COX inhibitors on the basis of their predicted affinity and interaction profile to COX enzymes. The analyses also provided insights into the role of induced fit in COX enzymes. While inhibitors produce different local structural changes at the COX ligand binding site, induced fit allows inhibitors in diverse chemical classes to share characteristic interaction patterns that ensure key contacts to be achieved. Different interaction patterns may also be associated with different inhibitory mechanisms.