2D and 3D Quantitative Structure-Activity Relationship Studies on a Series ofbis-Pyridinium Compounds as Choline Kinase Inhibitors

Studies on the IC50 of Metabolically Stable 1-(3,3-diphenylpropyl)- piperidinyl Amides and Ureas as Human CCR5 Receptor Antagonists Based on QSAR

Background:

In this study, modulators of human Chemotactic cytokine receptor 5 (CCR5) were described using a quantitative structure-activity relationship model (QSAR). This model was based on the molecule’s chemical structure.

Methods::

All 56 compounds of CCR5 receptor antagonists were randomly separated into two sets, 43 were reserved for training and the other 13 for testing. In the course of this study, molecular models were drawn using ChemDraw software. By means of Hyperchem software as well as optimized with both AM1 (semi-empirical self-consistent-field molecular orbital) and MM+ (molecular mechanics plus force field), molecular models were described through numerous descriptors using CODESSA software.

Results:

Linear models were obtained by Heuristic Method (HM) software and nonlinear models were obtained using APS software with optimal descriptor combinations used to build linear QSAR models, involving a group of selected descriptors. As a result, values of the above two different sets were shown to result from 0.82 in testing and 0.86 in training in HM while 0.83 in testing and 0.88 in training in Gene Expression Programming (GEP).

Conclusion:

From this method, the activity of molecules could be predicted, and the molecular structure could be changed to alter its IC50, avoiding the testing of large numbers of compounds.

A QSAR and molecular modelling study towards new lead finding: polypharmacological approach to Mycobacterium tuberculosis

Developing effective inhibitors against Mycobacterium tuberculosis (Mtb) is a challenging task, primarily due to the emergence of resistant strains. In this study, we have proposed and implemented an in silico guided polypharmacological approach, which is expected to be effective against resistant strains by simultaneously inhibiting several potential Mtb drug targets. A combination of pharmacophore and QSAR based virtual screening strategy taking three key targets such as InhA (enoyl-acyl-carrier-protein reductase), GlmU (N-acetyl-glucosamine-1-phosphate uridyltransferase) and DapB (dihydrodipicolinate reductase) have resulted in initial 784 hits from Asinex database of 435,000 compounds. These hits were further subjected to docking with 33 Mtb druggable targets. About 110 potential polypharmacological hits were taken by integrating the aforementioned screening protocols. Further screening was conducted by taking various parameters and properties such as cell permeability, drug-likeness, drug-induced phospholipidosisand structural alerts. A consensus analysis has yielded 59 potential hits that pass through all the filters and can be prioritized for effective drug-resistant tuberculosis. This study proposes about nine potential hits which are expected to be promising molecules, having not only drug-like properties, but also being effective against multiple Mtb targets.

Modeling the permeability of drug-like molecules through the cell wall of Mycobacterium tuberculosis: an analogue based approach

The emergence of drug resistant strains of Mycobacterium Tuberculosis (Mtb) accentuates the urgent need for the development of novel antitubercular drugs. The major causes of drug resistance are genetic mutations, the influx-efflux transporter system, and the complex cell wall system of Mtb, which can function as permeability barriers. The driving force for permeability of small molecules through a biological system depends on various physicochemical factors. To understand the permeability of small molecules and subsequent cell inhibition, we have developed predictive QSAR models based on reported enzyme-based (IC₅₀) and cell-based (MIC) Mtb inhibitory data. The compounds that are highly active in enzyme-based assays and have significant variation in cell-based assays are assumed to possess the required permeability through the Mtb cell wall. The obtained models suggest the importance of molecular connectivity, lipophilicity (log P, size, shape), electrotopology (relative atomic mass, polarizability, electronegativity, ionization potential, atomic charges, van der Waals volume, hybridization, hydrogen bond acceptors/donors, number of fused rings) and functional groups (hydroxyl groups, primary and secondary amines) of a molecule in determining both its inhibitory potency and Mtb cell permeability. The models were validated with known Mtb inhibitors (9804) collected from the ChEMBL database, Mtb drugs (27) and clinical candidates (5). Further, these validated models were used to screen and prioritize a large database of compounds, including Zinc (152 128), Asinex (435 215), DrugBank (6531) and antimicrobial compounds (1324). The results obtained from 2D-QSAR analysis could improve our understanding towards Mtb cell permeability, which may aid in the rational design of novel potent molecules for tuberculosis (TB).

Novel benzothiazine-piperazine derivatives by peptide-coupling as potential anti-proliferative agents

In an attempt to develop potential and selective anti-proliferative agents, a series of novel benzothiazine-piperazine derivatives 8a-i and 10a-g were synthesized by coupling of 2H-1,4-benzothiazin-3(4H)-one with various amines 7a-i and 9a-g in excellent yields and evaluated for their in vitro anti-proliferative activity against four cancer cell lines, HeLa (cervical), MIAPACA (pancreatic), MDA-MB-231 (breast) and IMR32 (neuroblastoma). In vitro inhibitory activity indicated that compounds 8a, 8d, 8g, 10a, 10b, 10e, 10f were found to be good anti-proliferative agents. Among them the derivatives 8g, 10e and 10f were found to be the most active members exhibiting remarkable growth inhibitory activity. Molecular docking was undertaken to investigate the probable binding mode and key active site interactions in HDAC8 and EHMT2 proteins. The docking results are complementary to the experimental results.

Choline kinase alpha-Putting the ChoK-hold on tumor metabolism

It is well established that lipid metabolism is drastically altered during tumor development and response to therapy. Choline kinase alpha (ChoKα) is a key mediator of these changes, as it represents the first committed step in the Kennedy pathway of phosphatidylcholine biosynthesis and ChoKα expression is upregulated in many human cancers. ChoKα activity is associated with drug resistant, metastatic, and malignant phenotypes, and represents a robust biomarker and therapeutic target in cancer. Effective ChoKα inhibitors have been developed and have recently entered clinical trials. ChoKα's clinical relevance was, until recently, attributed solely to its production of second messenger intermediates of phospholipid synthesis. The recent discovery of a non-catalytic scaffolding function of ChoKα may link growth receptor signaling to lipid biogenesis and requires a reinterpretation of the design and validation of ChoKα inhibitors. Advances in positron emission tomography, magnetic resonance spectroscopy, and optical imaging methods now allow for a comprehensive understanding of ChoKα expression and activity in vivo. We will review the current understanding of ChoKα metabolism, its role in tumor biology and the development and validation of targeted therapies and companion diagnostics for this important regulatory enzyme. This comes at a critical time as ChoKα-targeting programs receive more clinical interest.

Impact of cation-π interactions on the cell voltage of carbon nanotube-based Li batteries

Carbon nanotube (CNT)-based Li batteries have attracted wide attention because of their high capacity, high cyclability and high energy density and are believed to be one of the most promising electrochemical energy storage systems. In CNT-based Li batteries, the main interaction between the Li(+) ions and the CNT is the cation-π interaction. However, up to now, it is still not clear how this interaction affects the storage characteristics of CNT-based Li batteries. Here, using density functional theory (DFT) calculations, we report a highly favorable impact of cation-π interactions on the cell voltage of CNT-based Li batteries. Considering both Li(+)-π interaction and Li-π interaction, we show that cell voltage enhances with the increase of the CNT diameter. In addition, when the Li(+) ion adsorbs on the external wall, the cell voltage is larger than that when it adsorbs on the internal wall. This suggests that CNTs with a large diameter and a low array density are more advantageous to enhance storage performance of CNT-based Li batteries. Compared with Li(+) ions on the (4,4) CNT internal wall, the cell voltage of Li(+) on the (10,10) CNT external wall is 0.55 V higher, which indicates an improvement of about 38%. These results will be helpful for the design of more efficient CNT-based Li batteries.

Luotonin-A based quinazolinones cause apoptosis and senescence via HDAC inhibition and activation of tumor suppressor proteins in HeLa cells

A series of novel quinazolinone hybrids were synthesized by employing click chemistry and evaluated for anti-proliferative activities against MCF-7, HeLa and K562 cell lines. Among these cell lines, HeLa cells were found to respond effectively to these quinazolinone hybrids with IC50 values ranging from 5.94 to 16.45 μM. Some of the hybrids (4q, 4r, 4e, 4k, 4t, 4w) with promising anti-cancer activity were further investigated for their effects on the cell cycle distribution. FACS analysis revealed the G1 cell cycle arrest nature of these hybrids. Further to assess the senescence inducing ability of these compounds, a senescence associated β-gal assay was performed. The senescence inducing nature of these compounds was supported by the effect of hybrid (4q) on p16 promoter activity, the marker for senescence. Moreover, cells treated with most effective compound (4q) show up-regulation of p53, p21 and down-regulation of HDAC-1, HDAC-2, HDAC-5 and EZH2 mRNA levels. Docking results suggest that, the triazole nitrogen showed Zn(+2) mediated interactions with the histidine residue of HDACs.

Scaffold and cell line based approaches for QSAR studies on anticancer agents

Importance of 2D QSAR in drug discovery, lower number of descriptors containing models shows best statistical parameters, number of involved scaffolds in models affects the statistical values.

FDA approved drugs complexed to their targets: evaluating pose prediction accuracy of docking protocols

Efficient drug discovery programs can be designed by utilizing existing pools of knowledge from the already approved drugs. This can be achieved in one way by repositioning of drugs approved for some indications to newer indications. Complex of drug to its target gives fundamental insight into molecular recognition and a clear understanding of putative binding site. Five popular docking protocols, Glide, Gold, FlexX, Cdocker and LigandFit have been evaluated on a dataset of 199 FDA approved drug-target complexes for their accuracy in predicting the experimental pose. Performance for all the protocols is assessed at default settings, with root mean square deviation (RMSD) between the experimental ligand pose and the docked pose of less than 2.0 Å as the success criteria in predicting the pose. Glide (38.7 %) is found to be the most accurate in top ranked pose and Cdocker (58.8 %) in top RMSD pose. Ligand flexibility is a major bottleneck in failure of docking protocols to correctly predict the pose. Resolution of the crystal structure shows an inverse relationship with the performance of docking protocol. All the protocols perform optimally when a balanced type of hydrophilic and hydrophobic interaction or dominant hydrophilic interaction exists. Overall in 16 different target classes, hydrophobic interactions dominate in the binding site and maximum success is achieved for all the docking protocols in nuclear hormone receptor class while performance for the rest of the classes varied based on individual protocol.

3D-QSAR STUDIES OF SUBSTITUTED 4-ARYL/HETEROARYL-4H-CHROMENES AS APOPTOSIS INDUCERS USING CoMFA AND CoMSIA

Three-dimensional (3D) quantitative structure–activity relationships (QSARs) of 36 apoptosis inducers, substituted 4-aryl/heteroaryl-4H-chromenes with anticancer activity against human breast cancer cell lines T47D, have been studied by using methods of comparative molecular field analysis (CoMFA) and comparative molecular similarity analysis (CoMSIA). The established 3D-QSAR models in training set show not only significant statistical quality, but also predictive ability, with high correlation coefficient (R2) values and cross-validation coefficient (q2) values: CoMFA (R2, q2: 0.944, 0.747), CoMSIA (R2, q2: 0.944, 0.704). Moreover, the predictive abilities of the CoMFA and CoMSIA models were further confirmed by a test set, giving the predictive correlation coefficients ([Formula: see text] values) of 0.845 and 0.851, respectively. Based on the CoMFA and CoMSIA contour map analyses, some key factors responsible for anticancer activity of this series of compounds have been found as follows: the steric interaction plays a decisive role in determining the anticancer activities of these compounds; bulky groups as substituent R 1 are not tolerated; in addition to a steric moderation, higher degree of electropositivity and hydrophobicity on the terminal alkyl of substituent R 2 might be favorable to the activity; the substituent R 3 should be hydrophobic; bulky and strong electron withdrawing groups for the substituent R 4 are not advantageous to the activity; simultaneously introducing large electronegative atoms as hydrogen-acceptors to the first atoms of the substituents R 5 and R 6 may increase the activity, but substituents R 5 and R 6 with a linking group – OCH 2 O – may decrease the activity. Such results can offer some useful theoretical references for understanding the action mechanism, designing more potent derivatives, and predicting their activities prior to synthesis.

2D/3D-QSAR STUDY ON ANALOGUES OF 2-METHOXYESTRADIOL WITH ANTICANCER ACTIVITY

Two-dimensional (2D) and three-dimensional (3D) quantitative structure–activity relationships (QSARs) of 23 analogs of 2-Methoxyestradiol with anticancer activity (expressed as p GI50) against MCF-7 human breast cancer cells have been studied by using a combined method of the DFT, MM2 and statistics for 2D, as well as the comparative molecular field analysis (CoMFA) for 3D. The established 2D-QSAR model in training set shows not only significant statistical quality, but also predictive ability, with the square of adjusted correlation coefficient [Formula: see text] and the square of the cross-validation coefficient (q2= 0.779). The same model was further applied to predict p GI50values of the four compounds in the test set, and the resulting [Formula: see text] being as high as 0.827, further confirms that this 2D-QSAR model has high predictive ability for this kind of compound. The 3D-QSAR model also shows good correlative and predictive capabilities in terms of R2(0.927) and q2(0.786) obtained from CoMFA model. The results that 2D- and 3D-QSAR analyses accord with each other, suggest that the electrostatic interaction plays a decisive role in determining the anticancer activity of the studied compounds, and that increasing the negative charge of substituent R2and the positive charge of substituents linking to C17as well as decreasing the size of substituent R1are advantageous to improving the cytotoxicity. Such results can offer some useful theoretical references for directing the molecular design and understanding the action mechanism of this kind of compound with anticancer activity.

Analogue-based approaches in anti-cancer compound modelling: the relevance of QSAR models

BACKGROUND

QSAR is among the most extensively used computational methodology for analogue-based design. The application of various descriptor classes like quantum chemical, molecular mechanics, conceptual density functional theory (DFT)- and docking-based descriptors for predicting anti-cancer activity is well known. Although in vitro assay for anti-cancer activity is available against many different cell lines, most of the computational studies are carried out targeting insufficient number of cell lines. Hence, statistically robust and extensive QSAR studies against 29 different cancer cell lines and its comparative account, has been carried out.

RESULTS

The predictive models were built for 266 compounds with experimental data against 29 different cancer cell lines, employing independent and least number of descriptors. Robust statistical analysis shows a high correlation, cross-validation coefficient values, and provides a range of QSAR equations. Comparative performance of each class of descriptors was carried out and the effect of number of descriptors (1-10) on statistical parameters was tested. Charge-based descriptors were found in 20 out of 39 models (approx. 50%), valency-based descriptor in 14 (approx. 36%) and bond order-based descriptor in 11 (approx. 28%) in comparison to other descriptors. The use of conceptual DFT descriptors does not improve the statistical quality of the models in most cases.

CONCLUSION

Analysis is done with various models where the number of descriptors is increased from 1 to 10; it is interesting to note that in most cases 3 descriptor-based models are adequate. The study reveals that quantum chemical descriptors are the most important class of descriptors in modelling these series of compounds followed by electrostatic, constitutional, geometrical, topological and conceptual DFT descriptors. Cell lines in nasopharyngeal (2) cancer average R2 = 0.90 followed by cell lines in melanoma cancer (4) with average R2 = 0.81 gave the best statistical values.

Design of 1-arylsulfamido-2-alkylpiperazine derivatives as secreted PLA2 inhibitors

Structure and analog based analysis of 3D-QSAR, CoMFA and CoMSIA, along with different docking protocols were used to evaluate the structure activity relationship of 26 analogues of 1-aryl sulfamido-2-alkyl piperazines to model the activities of group I and II secreted phospholipases A(2) (sPLA(2)s) and probe into the chemical space and nature of receptor--ligand interactions. The best CoMFA model yields cross-validated (q(2)) and conventional correlation coefficients (r(2)) of 0.703 and 0.962 respectively whereas CoMSIA model yields q(2) and r(2) values of 0.408 and 0.922 respectively, followed by docking analysis using FlexX and GOLD methodologies on the X-ray structure of human and bovine PLA(2)s. A comparative study was made to find out the differences in the active site residues of both PLA(2)s. The information enunciated from the analysis of CoMFA and CoMSIA maps and docking results were analyzed and employed in the design of 29 new ligands using molecules 4, 21, 22 from the initial set as templates. New ligands for group I and II secreted phospholipases A(2) (sPLA(2)s) have been thus designed based on the 32 analogues of 1-aryl sulfamido-2-alkyl piperazine with a cursory note on its synthetic feasibility. Molecular modeling studies indicate that the newly designed ligands are expected to show high affinity and experimental efforts in this direction is highly rewarding.

Antibacterial activity and QSAR of chalcones against biofilm-producing bacteria isolated from marine waters

Biofouling in the marine environment is a major problem. In this study, three marine organisms, namely Bacillus flexus (LD1), Pseudomonas fluorescens (MD3) and Vibrio natriegens (MD6), were isolated from biofilms formed on polymer and metal surfaces immersed in ocean water. Phylogenetic analysis of these three organisms indicated that they were good model systems for studying marine biofouling. The in vitro antifouling activity of 47 synthesized chalcone derivatives was investigated by estimating the minimum inhibitory concentration against these organisms using a twofold dilution technique. Compounds C-5, C-16, C-24, C-33, C-34 and C-37 were found to be the most active. In the majority of the cases it was found that these active compounds had hydroxyl substitutions. A quantitative structure-activity relationship (QSAR) was developed after dividing the total data into training and test sets. The statistical measures r(2), [image omitted] (>0.6) q(2) (>0.5) and the F-ratio were found to be satisfactory. Spatial, structural and electronic descriptors were found to be predominantly affecting the antibiofouling activity of these compounds. Among the spatial descriptors, Jurs descriptors showed their contribution in all the three antibacterial QSARs.

Comparison of the cellular and biochemical properties of Plasmodium falciparum choline and ethanolamine kinases

The proliferation of the malaria-causing parasite Plasmodium falciparum within the erythrocyte is concomitant with massive phosphatidylcholine and phosphatidylethanolamine biosynthesis. Based on pharmacological and genetic data, de novo biosynthesis pathways of both phospholipids appear to be essential for parasite survival. The present study characterizes PfCK (P. falciparum choline kinase) and PfEK (P. falciparum ethanolamine kinase), which catalyse the first enzymatic steps of these essential metabolic pathways. Recombinant PfCK and PfEK were expressed as His6-tagged fusion proteins from overexpressing Escherichia coli strains, then purified to homogeneity and characterized. Using murine polyclonal antibodies against recombinant kinases, PfCK and PfEK were shown to be localized within the parasite cytoplasm. Protein expression levels increased during erythrocytic development. PfCK and PfEK appeared to be specific to their respective substrates and followed Michaelis-Menten kinetics. The Km value of PfCK for choline was 135.3+/-15.5 microM. PfCK was also able to phosphorylate ethanolamine with a very low affinity. PfEK was found to be an ethanolamine-specific kinase (Km=475.7+/-80.2 microM for ethanolamine). The quaternary ammonium compound hemicholinium-3 and an ethanolamine analogue, 2-amino-1-butanol, selectively inhibited PfCK or PfEK. In contrast, the bis-thiazolium compound T3, which was designed as a choline analogue and is currently in clinical trials for antimalarial treatment, affected PfCK and PfEK activities similarly. Inhibition exerted by T3 was competitive for both PfCK and PfEK and correlated with the impairment of cellular phosphatidylcholine biosynthesis. Comparative analyses of sequences and structures for both kinase types gave insights into their specific inhibition profiles and into the dual capacity of T3 to inhibit both PfCK and PfEK.

Potential choline kinase inhibitors: a molecular modeling study of bis-quinolinium compounds

Choline kinase (ChoK) is reported to involve in cell signaling pathways and cell growth by regulating the intermediate, phosphocholine (PCho), which is the first step to biosynthesis a membrane phospholipid, phosphatidylcholine. The PCho levels are overexpressed due to elevated activation of the protein under carcinogenesis conditions. ChoK has thus evolved as a novel target for various cancers and a range of compounds has been reported in this course as potent ChoK inhibitors. However, not much information is known about the binding site of the inhibitors. Therefore, we ventured to unravel the possible binding site of 39 bis-quinolinium inhibitors from which the structural requirement for better protein-ligand complex was delved. Molecular docking and 3D-QSAR studies namely comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were performed on the series. The knowledge of the active site was obtained from the site id search and molcad surface calculations of Sybyl, which was further considered for docking studies. In 3D-QSAR, the best predictions were obtained from the model where 29 compounds were considered in the training set and remaining 10 in the test set. The best CoMFA statistics were obtained with r(2) of 0.99 and q(2) of 0.81 while, CoMSIA was resulted with r(2) of 0.98 and q(2) of 0.77. A comparative analysis was done with the resulted 3D-QSAR maps and the docked poses by overlaying the maps on the active site residues. Since, there is no reported ligand co-crystallized structure of ChoK the present study provides valuable clues on the binding conformation of the ligand and its interactions with the active site.