Pharmacophore definition and three-dimensional quantitative structure-activity relationship study on structurally diverse prostacyclin receptor agonists

Synthesis, Biological Evaluation and Molecular Docking Studies of 5-indolylmethylen-4-oxo-2-thioxothiazolidine Derivatives

Background: Infectious diseases represent a significant global strain on public health security and impact on socio-economic stability all over the world. The increasing resistance to the current antimicrobial treatment has resulted in the crucial need for the discovery and development of novel entities for the infectious treatment with different modes of action that could target both sensitive and resistant strains. Methods: Compounds were synthesized using the classical organic chemistry methods. Prediction of biological activity spectra was carried out using PASS and PASS-based web applications. Pharmacophore modeling in LigandScout software was used for quantitative modeling of the antibacterial activity. Antimicrobial activity was evaluated using the microdilution method. AutoDock 4.2® software was used to elucidate probable bacterial and fungal molecular targets of the studied compounds. Results: All compounds exhibited better antibacterial potency than ampicillin against all bacteria tested. Three compounds were tested against resistant strains MRSA, P.aeruginosa and E.coli and were found to be more potent than MRSA than reference drugs. All compounds demonstrated a higher degree of antifungal activity than the reference drugs bifonazole (6–17-fold) and ketoconazole (13–52-fold). Three of the most active compounds could be considered for further development of the new, more potent antimicrobial agents. Conclusion: Compounds 5b (Z)-3-(3-hydroxyphenyl)-5-((1-methyl-1H-indol-3-yl)methylene)-2-thioxothiazolidin-4-one and 5g (Z)-3-[5-(1H-Indol-3-ylmethylene)-4-oxo-2-thioxo-thiazolidin-3-yl]-benzoic acid as well as 5h (Z)-3-(5-((5-methoxy-1H-indol-3-yl)methylene)-4-oxo-2-thioxothiazolidin-3-yl)benzoic acid can be considered as lead compounds for further development of more potent and safe antibacterial and antifungal agents.

Pharmacophore modeling coupled with molecular dynamic simulation approach to identify new leads for meprin-β metalloprotease

Human meprin beta metalloprotease, a small subgroup of the astacin family, is a potent drug target for the treatment of several disorders such as fibrosis, neurodegenerative disease in particular Alzheimer and inflammatory bowel diseases. In this study, a ligand-based pharmacophore approach has been used for the selection of potentially active compounds to understand the inhibitory activities of meprin-β by using the sulfonamide scaffold based inhibitors. Using this dataset, a pharmacophore model (Hypo1) was selected on the basis of a highest correlation coefficient (0.959), lowest total cost (105.89) and lowest root mean square deviation (1.31 Å) values. All the pharmacophore hypotheses generated from the candidate inhibitors comprised four features: two hydrogen-bond acceptor, one hydrogen-bond donor and one zinc binder feature. The best validated pharmacophore model (Hypo1) was used for virtual screening of compounds from several databases. The selective hit compounds were filtered by drug likeness property, acceptable ADMET profile, molecular docking and DFT study. Molecular dynamic simulations with the final 10 hit compounds revealed that a large number of non-covalent interactions were formed with the active site and specificity sub-pockets of the meprin beta metalloprotease. This study assists in the development of the new lead molecules as well as gives a better understanding of their interaction with meprin-β.

Tubulin inhibitors: pharmacophore modeling, virtual screening and molecular docking

AIM

To construct a quantitative pharmacophore model of tubulin inhibitors and to discovery new leads with potent antitumor activities.

METHODS

Ligand-based pharmacophore modeling was used to identify the chemical features responsible for inhibiting tubulin polymerization. A set of 26 training compounds was used to generate hypothetical pharmacophores using the HypoGen algorithm. The structures were further validated using the test set, Fischer randomization method, leave-one-out method and a decoy set, and the best model was chosen to screen the Specs database. Hit compounds were subjected to molecular docking study using a Molecular Operating Environment (MOE) software and to biological evaluation in vitro.

RESULTS

Hypo1 was demonstrated to be the best pharmacophore model that exhibited the highest correlation coefficient (0.9582), largest cost difference (70.905) and lowest RMSD value (0.6977). Hypo1 consisted of one hydrogen-bond acceptor, a hydrogen-bond donor, a hydrophobic feature, a ring aromatic feature and three excluded volumes. Hypo1 was validated with four different methods and had a goodness-of-hit score of 0.81. When Hypo1 was used in virtual screening of the Specs database, 952 drug-like compounds were revealed. After docking into the colchicine-binding site of tubulin, 5 drug-like compounds with the required interaction with the critical amino acid residues and the binding free energies < -4 kcal/mol were selected as representative leads. Compounds 1 and 3 exhibited inhibitory activity against MCF-7 human breast cancer cells in vitro.

CONCLUSION

Hypo1 is a quantitative pharmacophore model for tubulin inhibitors, which not only provides a better understanding of their interaction with tubulin, but also assists in discovering new potential leads with antitumor activities.

Pharmacophore modeling and virtual screening for the discovery of new type 4 cAMP phosphodiesterase (PDE4) inhibitors

Type 4 cAMP phosphodiesterase (PDE4) inhibitors show a broad spectrum of anti-inflammatory effects in almost all kinds of inflamed cells, by an increase in cAMP levels which is a pivotal second messenger responsible for various biological processes. These inhibitors are now considered as the potential drugs for treatment of chronic inflammatory diseases. However, some recently marketed inhibitors e.g., roflumilast, have shown adverse effects such as nausea and emesis, thus restricting its use. In order to identify novel PDE4 inhibitors with improved therapeutic indexes, a highly correlating (r = 0.963930) pharmacophore model (Hypo1) was established on the basis of known PDE4 inhibitors. Validated Hypo1 was used in database screening to identify chemical with required pharmacophoric features. These compounds are further screened by using the rule of five, ADMET and molecular docking. Finally, twelve hits which showed good results with respect to following properties such as estimated activity, calculated drug-like properties and scores were proposed as potential leads to inhibit the PDE4 activity. Therefore, this study will not only assist in the development of new potent hits for PDE4 inhibitors, but also give a better understanding of their interaction with PDE4. On a wider scope, this will be helpful for the rational design of novel potent enzyme inhibitors.

1-Aryl-3-(1H-imidazol-1-yl)propan-1-ol esters: synthesis, anti-Candida potential and molecular modeling studies

BACKGROUND

An increased incidence of fungal infections, both invasive and superficial, has been witnessed over the last two decades. Candida species seem to be the main etiology of nosocomial fungal infections worldwide with Candida albicans, which is commensal in healthy individuals, accounting for the majority of invasive Candida infections with about 30-40% of mortality.

RESULTS

New aromatic and heterocyclic esters 5a-k of 1-aryl-3-(1H-imidazol-1-yl)propan-1-ols 4a-d were successfully synthesized and evaluated for their anti-Candida potential. Compound 5a emerged as the most active congener among the newly synthesized compounds 5a-k with MIC value of 0.0833 μmol/mL as compared with fluconazole (MIC value >1.6325 μmol/mL). Additionally, molecular modeling studies were conducted on a set of anti-Candida albicans compounds.

CONCLUSION

The newly synthesized esters 5a-k showed more potent anti-Candida activities than fluconazole. Compounds 7 and 8 revealed significant anti-Candida albicans activity and were able to effectively satisfy the proposed pharmacophore geometry, using the energy accessible conformers (Econf < 20 kcal/mol).

Exploration of the antiplatelet activity profile of betulinic acid on human platelets

Betulinic acid, a natural pentacyclic triterpene acid, presents a diverse mode of biological actions including antiretroviral, antibacterial, antimalarial, and anti-inflammatory activities. The potency of betulinic acid as an inhibitor of human platelet activation was evaluated, and its antiplatelet profile against in vitro platelet aggregation, induced by several platelet agonists (adenosine diphosphate, thrombin receptor activator peptide-14, and arachidonic acid), was explored. Flow cytometric analysis was performed to examine the effect of betulinic acid on P-selectin membrane expression and PAC-1 binding to activated platelets. Betulinic acid potently inhibits platelet aggregation and also reduced PAC-1 binding and the membrane expression of P-selectin. Principal component analysis was used to screen, on the chemical property space, for potential common pharmacophores of betulinic acid with approved antithrombotic drugs. A common pharmacophore was defined between the NMR-derived structure of betulinic acid and prostacyclin agonists (PGI2), and the importance of its carboxylate group in its antiplatelet activity was determined. The present results indicate that betulinic acid has potential use as an antithrombotic compound and suggest that the mechanism underlying the antiplatelet effects of betulinic acid is similar to that of the PGI2 receptor agonists, a hypothesis that deserves further investigation.

Molecular mechanisms regulating the vascular prostacyclin pathways and their adaptation during pregnancy and in the newborn

Prostacyclin (PGI(2)) is a member of the prostanoid group of eicosanoids that regulate homeostasis, hemostasis, smooth muscle function and inflammation. Prostanoids are derived from arachidonic acid by the sequential actions of phospholipase A(2), cyclooxygenase (COX), and specific prostaglandin (PG) synthases. There are two major COX enzymes, COX1 and COX2, that differ in structure, tissue distribution, subcellular localization, and function. COX1 is largely constitutively expressed, whereas COX2 is induced at sites of inflammation and vascular injury. PGI(2) is produced by endothelial cells and influences many cardiovascular processes. PGI(2) acts mainly on the prostacyclin (IP) receptor, but because of receptor homology, PGI(2) analogs such as iloprost may act on other prostanoid receptors with variable affinities. PGI(2)/IP interaction stimulates G protein-coupled increase in cAMP and protein kinase A, resulting in decreased [Ca(2+)](i), and could also cause inhibition of Rho kinase, leading to vascular smooth muscle relaxation. In addition, PGI(2) intracrine signaling may target nuclear peroxisome proliferator-activated receptors and regulate gene transcription. PGI(2) counteracts the vasoconstrictor and platelet aggregation effects of thromboxane A(2) (TXA(2)), and both prostanoids create an important balance in cardiovascular homeostasis. The PGI(2)/TXA(2) balance is particularly critical in the regulation of maternal and fetal vascular function during pregnancy and in the newborn. A decrease in PGI(2)/TXA(2) ratio in the maternal, fetal, and neonatal circulation may contribute to preeclampsia, intrauterine growth restriction, and persistent pulmonary hypertension of the newborn (PPHN), respectively. On the other hand, increased PGI(2) activity may contribute to patent ductus arteriosus (PDA) and intraventricular hemorrhage in premature newborns. These observations have raised interest in the use of COX inhibitors and PGI(2) analogs in the management of pregnancy-associated and neonatal vascular disorders. The use of aspirin to decrease TXA(2) synthesis has shown little benefit in preeclampsia, whereas indomethacin and ibuprofen are used effectively to close PDA in the premature newborn. PGI(2) analogs have been used effectively in primary pulmonary hypertension in adults and have shown promise in PPHN. Careful examination of PGI(2) metabolism and the complex interplay with other prostanoids will help design specific modulators of the PGI(2)-dependent pathways for the management of pregnancy-related and neonatal vascular disorders.

Schiff bases of indoline-2,3-dione (isatin) with potential antiproliferative activity

Background

Cancer is one of the most dreaded diseases and it is a leading cause of mankind death worldwide. Recent reports documented a remarkable antiproliferative activity of isatin nucleus against various cancer cell lines. The current work describes the antiproliferative activity of Schiff bases of combinatorial mixtures of the isatin derivatives M1-M22 as well as the individual compounds 1-11(A-K) of these combinatorial mixtures.

Results

The designed combinatorial library composed from eleven hydrazides A-K and eleven isatin derivatives 1-11 has been synthesized to formally generate 22 mixtures, M1-M22 of 121 Schiff bases, and their antiproliferative activity against K562 chronic myelogenous leukemia cells was evaluated. The indexed method of analysis of the prepared library was applied to elucidate the active components in the tested mixtures M1-M22. The predictions from the crossing procedure was validated through evaluation of the antiproliferative activity of individual compounds 1-11(A-K) of the library. Individual compounds 1-11(A-K) were also evaluated against the non-tumorigenic MCF-12A cell line to investigate their selectivity. A pharmacophore model was developed to further optimize the antiproliferative activity among this series of compounds.

Conclusions

Variable antiproliferative activity was revealed with the investigated mixtures M1-M22 and the individual compounds 1-11(A-K). Most of the tested mixtures and several individual Schiff bases displayed high potency with IC₅₀ values in the low micromolar range. A considerable selectivity of some individual compounds to the tumorigenic K562 cell line compared with the non-tumorigenic MCF-12A cell line was observed as indicated by their selectivity index (SI).

Development, evaluation and application of 3D QSAR Pharmacophore model in the discovery of potential human renin inhibitors

Background

Renin has become an attractive target in controlling hypertension because of the high specificity towards its only substrate, angiotensinogen. The conversion of angiotensinogen to angiotensin I is the first and rate-limiting step of renin-angiotensin system and thus designing inhibitors to block this step is focused in this study.

Methods

Ligand-based quantitative pharmacophore modeling methodology was used in identifying the important molecular chemical features present in the set of already known active compounds and the missing features from the set of inactive compounds. A training set containing 18 compounds including active and inactive compounds with a substantial degree of diversity was used in developing the pharmacophore models. A test set containing 93 compounds, Fischer randomization, and leave-one-out methods were used in the validation of the pharmacophore model. Database screening was performed using the best pharmacophore model as a 3D structural query. Molecular docking and density functional theory calculations were used to select the hit compounds with strong molecular interactions and favorable electronic features.

Results

The best quantitative pharmacophore model selected was made of one hydrophobic, one hydrogen bond donor, and two hydrogen bond acceptor features with high a correlation value of 0.944. Upon validation using an external test set of 93 compounds, Fischer randomization, and leave-one-out methods, this model was used in database screening to identify chemical compounds containing the identified pharmacophoric features. Molecular docking and density functional theory studies have confirmed that the identified hits possess the essential binding characteristics and electronic properties of potent inhibitors.

Conclusion

A quantitative pharmacophore model of predictive ability was developed with essential molecular features of a potent renin inhibitor. Using this pharmacophore model, two potential inhibitory leads were identified to be used in designing novel and future renin inhibitors as antihypertensive drugs.

Potent BACE-1 inhibitor design using pharmacophore modeling, in silico screening and molecular docking studies

BACKGROUND

Beta-site amyloid precursor protein cleaving enzyme (BACE-1) is a single-membrane protein belongs to the aspartyl protease class of catabolic enzymes. This enzyme involved in the processing of the amyloid precursor protein (APP). The cleavage of APP by BACE-1 is the rate-limiting step in the amyloid cascade leading to the production of two peptide fragments Aβ40 and Aβ42. Among two peptide fragments Aβ42 is the primary species thought to be responsible for the neurotoxicity and amyloid plaque formation that lead to memory and cognitive defects in Alzheimer's disease (AD). AD is a ravaging neurodegenerative disorder for which no disease-modifying treatment is currently available. Inhibition of BACE-1 is expected to stop amyloid plaque formation and emerged as an interesting and attractive therapeutic target for AD.

METHODS

Ligand-based computational approach was used to identify the molecular chemical features required for the inhibition of BACE-1 enzyme. A training set of 20 compounds with known experimental activity was used to generate pharmacophore hypotheses using 3D QSAR Pharmacophore Generation module available in Discovery studio. The hypothesis was validated by four different methods and the best hypothesis was utilized in database screening of four chemical databases like Maybridge, Chembridge, NCI and Asinex. The retrieved hit compounds were subjected to molecular docking study using GOLD 4.1 program.

RESULTS

Among ten generated pharmacophore hypotheses, Hypo 1 was chosen as best pharmacophore hypothesis. Hypo 1 consists of one hydrogen bond donor, one positive ionizable, one ring aromatic and two hydrophobic features with high correlation coefficient of 0.977, highest cost difference of 121.98 bits and lowest RMSD value of 0.804. Hypo 1 was validated using Fischer randomization method, test set with a correlation coefficient of 0.917, leave-one-out method and decoy set with a goodness of hit score of 0.76. The validated Hypo 1 was used as a 3D query in database screening and retrieved 773 compounds with the estimated activity value <100 nM. These hits were docked into the active site of BACE-1 and further refined based on molecular interactions with the essential amino acids and good GOLD fitness score.

CONCLUSION

The best pharmacophore hypothesis, Hypo 1, with high predictive ability contains chemical features required for the effective inhibition of BACE-1. Using Hypo 1, we have identified two compounds with diverse chemical scaffolds as potential virtual leads which, as such or upon further optimization, can be used in the designing of new BACE-1 inhibitors.

Allosteric site on muscarinic acetylcholine receptors: identification of two amino acids in the muscarinic M2 receptor that account entirely for the M2/M5 subtype selectivities of some structurally diverse allosteric ligands in N-methylscopolamine-occupied receptors

Two epitopes have been identified recently to be responsible for the high-affinity binding of alkane-bisammonium and caracurine V type allosteric ligands to N-methylscopolamine (NMS)-occupied M2 muscarinic acetylcholine receptors, relative to M5 receptors: the amino acid M2-Thr423 at the top of transmembrane region (TM) 7 and an epitope comprising the second extracellular loop (o2) of the M2 receptor including the flanking regions of TM4 and TM5. We aimed to find out whether a single amino acid could account for the contribution of this epitope to binding affinity. Allosteric interactions were investigated in wild-type and mutant receptors in which the orthosteric binding site was occupied by [3H]NMS (5 mM Na,K,Pi buffer, pH 7.4, 23 degrees C). Using M2/M5 chimeric and point-mutated receptors, the relevant epitope was narrowed down to M2-Tyr177. A double point-mutated M2 receptor in which both M2-Tyr177 and M2-Thr423 were replaced by the corresponding amino acids of M5 revealed that these two amino acids account entirely for the (approximately 100-fold) M2/M5 selectivity of the alkane-bisammonium and the caracurine V type allosteric ligands. At NMS-free M2 receptors, the caracurine V derivative also displayed approximately 100-fold M2/M5 selectivity, but the double point mutation reduced the M2 affinity by only approximately 10-fold; thus, additional epitopes may influence selectivity for the free receptors. A three-dimensional model of the M2 receptor was used to simulate allosteric agent docking to NMS-occupied receptors. M2-Tyr177 and M2-Thr423 seem to be located near the junction of the allosteric and the orthosteric areas of the M2 receptor ligand binding cavity.