Pharmacophore modelling,molecular docking and virtual screening for EGFR (HER 1) tyrosine kinase inhibitors

A pharmacophore model has been developed using diverse classes of epidermal growth factor receptor (EGFR) tyrosine kinase (TK) inhibitors useful in the treatment of human tumours. Among the top 10 generated hypotheses, the second hypothesis, with one hydrogen bond acceptor, one ring aromatic and three hydrophobic features, was found to be the best on the basis of Cat Scramble validation as well as test set prediction (r _training?=?0.89, r _test?=?0.82). The model also maps well to the external test set molecules as well as clinically active molecules and corroborates the docking studies. Finally, 10 hits were identified as potential leads after virtual screening of ZINC database for EGFR TK inhibition. The study may facilitate the designing and discovery of novel EGFR TK inhibitors.     

A comparative QSAR analysis and molecular docking studies of phenyl piperidine derivatives as potent dual NK1R antagonists/serotonin transporter (SERT) inhibitors

Depression is a critical mood disorder that affects millions of patients. Available therapeutic antidepressant agents are associated with several undesirable side effects. Recently, it has been shown that Neurokinin 1 receptor (NK₁R) antagonists can potentiate the antidepressant effects of serotonin-selective reuptake inhibitors (SSRIs). In this study, a series of phenyl piperidine derivatives as potent dual NK₁R antagonists/serotonin transporter (SERT) inhibitors were applied to quantitative structure–activity relationship (QSAR) analysis. A collection of chemometrics methods such as multiple linear regression (MLR), factor analysis–based multiple linear regression (FA-MLR), principal component regression (PCR), and partial least squared combined with genetic algorithm for variable selection (GA-PLS) were applied to make relations between structural characteristics and NK₁R antagonism/SERT inhibitory of these compounds. The best multiple linear regression equation was obtained from GA-PLS and MLR for NK₁R and SERT, respectively. Based on the resulted model, an in silico-screening study was also conducted and new potent lead compounds based on new structural patterns were designed for both targets. Molecular docking studies of these compounds on both targets were also conducted and encouraging results were acquired. There was a good correlation between QSAR and docking results. The results obtained from validated docking studies indicate that the important amino acids inside the active site of the cavity that are responsible for essential interactions are Glu33, Asp395 and Arg26 for SERT and Ala30, Lys7, Asp31, Phe5 and Tyr82 for NK₁R receptors.     

Discrimination of Scrophulariae Radix according to geographical origin and determination of active constituents by near infrared spectroscopy (NIRS)

The potential of near infrared reflectance spectroscopy (NIR) was investigated for its ability to non-destructively discriminate the geographic origins of Scrophularia spp., Andong, Uisung and China. Application of principal component analysis to NIR spectra leads to a clear separation of Andong sample from the others. Moreover, the contents of two neuroprotective constituents of Scrophularia spp., 8-O-(E-p-methoxycinnamoyl)-harpagide (HG), and E-p-methoxycinnamic acid (MCA), were determined by HPLC-DAD. Partial least squares (PLS) regression of NIR spectra combined with these analytical reference data yield the development of calibration models for the contents of the two constituents. The correlation coefficients of prediction models for HG and MCA were > 0.87. These outcomes indicated that the NIRS could be useful for the discrimination of Scrophularia spp.     

Near-infrared spectroscopy and multivariate calibration for the quantitative determination of certain properties in the petrochemical industry

Near-infrared (NIR) spectroscopy in conjunction with chemometric techniques allows on-line monitoring in real time, which can be of considerable use in industry. If it is to be correctly used in industrial applications, generally some basic considerations need to be taken into account, although this does not always apply. This study discusses some of the considerations that would help evaluate the possibility of applying multivariate calibration in combination with NIR to properties of industrial interest. Examples of these considerations are whether there is a relation between the NIR spectrum and the property of interest, what the calibration constraints are and how a sample-specific error of prediction can be quantified. Various strategies for maintaining a multivariate model after it has been installed are also presented and discussed.     

Molecular docking studies on tetrahydroimidazo-[4,5,1-jk][1,4]-benzodiazepinone (TIBO) derivatives as HIV-1 NNRT inhibitors

At present, chemotherapy seems to be the main weapon in the arsenal of remedies for the ongoing crusade against AIDS. The mode of binding of the TIBO family of inhibitors has been of interest because these compounds do not fit the two-hinged-ring model as generally observed in the NNRTIs. Flexible docking simulations were performed with a series of 53 TIBO derivatives as NNRTIs. Binding preferences as well as the structural and energetic factors associated with them were studied. A good correlation (r ² = 0.849, q ² = 0.843) was observed between the biological activity and binding affinity of the compounds which suggest that the identified binding conformations of these inhibitors are reliable. Further screening of PubChem database yielded novel scaffolds. Our studies suggest that modifications to the TIBO group of inhibitors might enhance their binding efficacy and hence, potentially, their therapeutic utility.     

In-silico profiling,design, molecular docking computation,and drug kinetic model evaluation of novel curcumin derivatives as potential anticancer agents

The alarming trend of leukemia cell lines that are multidrug-resistant has prompted scientists to scramble for effective new anticancer treatments. Therefore, it remains an intriguing scientific task to optimize curcumin by trying to introduce molecular alteration to its vital structure to improve the biological effect against the P388 cell line or get around resistance phenomena. Regardless of the wide range of medications that are now being studied, Prednisone remains the most important and efficient part of chemotherapy that the WHO recommends. This article discusses the QSAR-oriented model and in silico assessment of some potent curcumin derivatives' anticancer activity against the P388 cell line. The solidity and propensity for prediction of the model were ensured by using stringent validation procedures. The activity of these derivatives was shown to be unrelated to lipophilicity, while shorter N-N distances and short substituents result in quite bioactive molecules. This information was used to design potent molecules that demonstrate good quality as per the assessment based on the Lead-Like Soft rule is acceptable for drug-like molecules. Also, molecule d2 does not possess any toxic effects risk alerts, suggesting drug-adherent conduct. While Prednisone the reference drug has a toxicity risk alert in red, suggesting non-adherent conduct for Prednisone. Hence, the novel molecules are promising anticancer agents.     

In silico design novel (5-imidazol-2-yl-4-phenylpyrimidin-2-yl)[2-(2-pyridylamino)ethyl]amine derivatives as inhibitors for glycogen synthase kinase 3 based on 3D-QSAR,molecular docking and molecular dynamics simulation

Glycogen Synthase Kinase 3 (GSK-3) is a member of cellular kinase with various functions, such as glucose regulation, cellular differentiation, neuronal function and cell apoptosis. It has been proved as an important therapeutic target in type 2 diabetes mellitus and Alzheimer's disease. To better understand their structure–activity relationships and mechanism of action, an integrated computational study, including three dimensional quantitative structure-activity relationship (3D-QSAR), molecular docking, and molecular dynamics (MD), was performed on 79 (5-Imidazol-2-yl-4-phenylpyrimidin-2-yl)[2-(2-pyridylamino)ethyl]amine GSK-3 inhibitors. In this paper, we constructed 3D-QSAR using comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA) method. The results showed that the CoMFA model (q ² = 0.743, r² = 0.980) and the CoMSIA model (q² = 0.813, r² = 0.976) had stable and reliable predictive ability. The electrostatic and H-bond donor fields play important roles in the models. The contour maps of the model visually showed the relationship between the activity of compounds and their three-dimensional structure. Molecular docking was used to identify the key amino acid residues at the active site of GSK-3 and explore its binding mode with ligands. Based on 3D-QSAR models, contour maps and the binding feature between GSK-3 and inhibitor, we designed 10 novel compounds with good potential activity and ADME/T profile. Molecular dynamics simulation results validated that Ile62, Val70 and Lys85 located in the active site play a key role for GSK-3 complexed with inhibitors. These results might provide important information for designing GSK-3 inhibitors with high activity.     

一种新的氨基酸描述子及其在肽QSAR中的应用

从天然氨基酸的25个结构与拓扑变量中经主成分分析得到一种新的氨基酸描述子——VSTV （principal component scores vector of structural and topological variables）.应用该描述子对以下3个体系,即血管紧张素转化酶抑制剂（2肽）、抗菌18肽和促凝血酶原激酶抑制剂（6～12肽）进行分子结构参数化表达,并在此基础上通过偏最小二乘回归(PLSR)建立定量构效关系(QSAR)模型,取得了优于文献的结果.模型的复相关系数（R2）和交互检验复相关系数（Q2）分别为0.789, 0.767; 0.996, 0.879; 0.981, 0.480.     

New salen-type manganese(III) Schiff base complexes derived from meso-1,2-diphenyl-1,2-ethylenediamine: In vitro anticancer activity,mechanism of action,and molecular docking studies

Four new manganese(III) Schiff base complexes (1–4) were synthesized and characterized. The complexes have general formula [MnClL^x] in which L represents a Schiff base ligand derived from condensation of meso-1,2-diphenyl-1,2-ethylenediamine with salicylaldehyde or its 3-OMe-, 5-Br-, or 5-OMe-derivatives (x = 1–4, respectively). The crystal structure of [MnClL¹] (1) was characterized by X-ray crystallography. The in vitro anticancer activity of these complexes was evaluated by MTT and apoptosis assays against human breast (MCF-7) and liver (Hep G2) cancer cells. The complexes exhibited considerable antiproliferative activity against both cell lines (IC₅₀ = 10.8–21.02 μM) comparable to cis-platin, except 4 (MCF-7). The highest activity was found for 1 with IC₅₀ values of 13.62 μM (MCF-7) and 10.8 μM (Hep G2). Flow cytometry experiments showed that 1 induced apoptosis on MCF-7 tumor cell line. Docking simulations using AUTODOCK were also carried out. The results showed that all complexes fitted into the minor groove region of DNA.     

Virtual screening of B-Raf kinase inhibitors: A combination of pharmacophore modelling,molecular docking, 3D-QSAR model and binding free energy calculation studies

B-Raf kinase has been identified as an important target in recent cancer treatment. In order to discover structurally diverse and novel B-Raf inhibitors (BRIs), a virtual screening of BRIs against ZINC database was performed by using a combination of pharmacophore modelling, molecular docking, 3D-QSAR model and binding free energy (ΔG_bind) calculation studies in this work. After the virtual screening, six promising hit compounds were obtained, which were then tested for inhibitory activities of A375 cell lines. In the result, five hit compounds show good biological activities (IC₅₀ < 50 μM). The present method of virtual screening can be applied to find structurally diverse inhibitors, and the obtained five structurally diverse compounds are expected to develop novel BRIs.     

Synthesis,anticancer activity and docking studies of pyrazoline and pyrimidine derivatives as potential epidermal growth factor receptor (EGFR) inhibitors

A search for anticancer agents has prompted the design and synthesis of new chalcone, pyrazoline and pyrimidine derivatives as potential epidermal growth factor receptor (EGFR) kinase inhibitors. These derivatives’ binding affinities were predicted by AutoDock, which showed that chalcone, pyrazoline and pyrimidine derivatives as EGFR-kinase inhibitors have good binding energies, ranging from ?10.91 to ?7.32 kcal/mol. These compounds were synthesized and characterized using elemental analysis (CHN analysis) and spectroscopic techniques (FTIR and NMR). Among the pyrazoline derivatives, 4Aiii has revealed a superior in vitro activity, inhibiting the EGFR kinase even at a low concentration of 0.19 μM compared to the pyrimidine derivative, 5Bii. In contrast, the cytotoxic effect of these derivatives was studied against hormonal and non-hormonal breast cancer cell lines. Most of the pyrazoline derivatives were able to express their cytotoxic effect efficiently against hormonal breast cancer but only one pyrimidine derivative managed to express its activity against hormonal breast cancer.     

A pharmacoinformatic approach on Cannabinoid receptor 2 (CB2) and different small molecules: Homology modelling,molecular docking,MD simulations,drug designing and ADME analysis

CB2 receptor belongs to the family of G-protein coupled receptors (GPCRs), which extensively controls a range of pointer transduction. CB2 plays an essential role in the immune system. It also associates in the pathology of different ailment conditions. In this scenario, the synthetic drugs are inducing side effects to the human beings after the drug use. Therefore, this study is seeking novel alternate drug molecules with least side effects than conventional drugs. The alternative drug molecules were chosen from the natural sources. These molecules were selected from cyanobacteria with the help of earlier research findings. The target and ligand molecules were obtained from recognized databases. The bioactive molecules are selected from various cyanobacterial species, which are selected by their biological and pharmacological properties, after, which we incorporated to the crucial findings such as homology modelling, molecular docking, MD simulations along with absorption, distribution, metabolism, and excretion (ADME) analysis. Initially, the homology modelling was performed to frame the target from unknown sequences of CB2, which revealed 44% of similarities and 66% of identities with the A2A receptor. Subsequently, the CB2 protein molecule has docked with already known and prepared bioactive molecules, agonists and antagonist complex. In the present study, the agonists (5) and antagonist (1) were also taken for comparing the results with natural molecules. At the end of the docking analysis, the cyanobacterial molecules and an antagonist TNC-201 are revealed better docking scores with well binding contacts than the agonists. Especially, the usneoidone shows better results than other cyanobacterial molecules, and it is very close docking scores with that of TCN-201. Therefore, the usneoidone has incorporated to MD simulation with Cannabinoid receptors 2 (CB2). In MD simulations, the complex (CB2 and usneoidone) reveals better stability in 30 ns. Based on the computational outcome, we concluded that usneoidone is an effectual and appropriate drug candidate for activating CB2 receptors and it will be serving as a better component for the complications of CB2. Moreover, these computational approaches can be motivated to discover novel drug candidates in the pharmacological and healthcare sectors.     

Chemometrics to chemical modeling: novel molecular distance-edge vector (λ) in alkanes and retention index of gas chromatography

Systematic studies are made on application of chemometrics to chemical modeling and/or molecular modeling as well as the regularity of retention index for gas chromatography (GC). A set of novel molecular graph theoretical parameters, called the molecular distance-edge (MDE) vector (λ), is proposed and found to be excellently correlated to retention index of GC for alkanes. The MDE parameters were tested by the multiple linear regression (MLR) estimation and prediction of the retention index of GC, and the results obtained are satisfactory.     

Adulteration of diesel/biodiesel blends by vegetable oil as determined by Fourier transform (FT) near infrared spectrometry and FT-Raman spectroscopy

In this work it has been shown that the routine ASTM methods (ASTM 4052, ASTM D 445, ASTM D 4737, ASTM D 93, and ASTM D 86) recommended by the ANP (the Brazilian National Agency for Petroleum, Natural Gas and Biofuels) to determine the quality of diesel/biodiesel blends are not suitable to prevent the adulteration of B2 or B5 blends with vegetable oils. Considering the previous and actual problems with fuel adulterations in Brazil, we have investigated the application of vibrational spectroscopy (Fourier transform (FT) near infrared spectrometry and FT-Raman) to identify adulterations of B2 and B5 blends with vegetable oils. Partial least square regression (PLS), principal component regression (PCR), and artificial neural network (ANN) calibration models were designed and their relative performances were evaluated by external validation using the F-test. The PCR, PLS, and ANN calibration models based on the Fourier transform (FT) near infrared spectrometry and FT-Raman spectroscopy were designed using 120 samples. Other 62 samples were used in the validation and external validation, for a total of 182 samples. The results have shown that among the designed calibration models, the ANN/FT-Raman presented the best accuracy (0.028%, w/w) for samples used in the external validation.     

Understanding and modulating cyclin-dependent kinase inhibitor specificity: molecular modeling and biochemical evaluation of pyrazolopyrimidinones as CDK2/cyclin A and CDK4/cyclin D1 inhibitors

Summary Cyclin-dependent kinases (CDKs) play a key role in regulating the cell cycle. The cyclins, their activating agents, and endogenous CDK inhibitors are frequently mutated in human cancers, making CDKs interesting targets for cancer chemotherapy. Our aim is the discovery of selective CDK4/cyclin D1 inhibitors. An ATP-competitive pyrazolopyrimidinone CDK inhibitor was identified by HTS and docked into a CDK4 homology model. The resulting binding model was consistent with available SAR and was validated by a subsequent CDK2/inhibitor crystal structure. An iterative cycle of chemistry and modeling led to a 70-fold improvement in potency. Small substituent changes resulted in large CDK4/CDK2 selectivity changes. The modeling revealed that selectivity is largely due to hydrogen-bonded interactions with only two kinase residues. This demonstrates that small differences between enzymes can efficiently be exploited in the design of selective inhibitors.     

p-TSA-catalyzed one-pot synthesis and docking studies of some 5H-indeno[1,2-b]quinoline-9,11(6H,10H)-dione derivatives as anticonvulsant agents

The present work describes a facile, one-pot three component environment friendly, green synthesis of a series of 5-(4-methoxyphenyl)-7,7-dimethyl-10-phenyl-7,8-dihydro-5H-indeno[1,2-b]quinoline-9,11(6H,10H)-dione derivatives 8(a-n). 1,3-indanedione, aryl-aldehyde and enaminone was thoroughly ground in the presence of catalytic amount of p-toluene sulfonic acid (p-TSA) to give the titled compounds in good yields. All the synthesized derivatives were evaluated for their anticonvulsant activity using the maximal electroshock (MES) method with phenytoin as a standard drug along with their neurotoxicity effect. Derivatives 8b, 8e and 8k exhibited significant anticonvulsant activity (P<0.001). The neurotoxicity study clearly revealed that all the tested compounds are non-toxic at a dose of 40 mg/kg. The molecular modeling studies also predicted good binding interactions of most active molecules with the serotonin 5-HT_2A receptor. Therefore, it can be safely concluded that synthesized derivatives 8(a-n) would represent useful leads for further investigation in the development of a new class of anticonvulsant agents.     

Improving the robustness of a partial least squares (PLS) model based on pure component selectivity analysis and range optimization: case study for the analysis of an etching solution containing hydrogen peroxide

Pure component selectivity analysis (PCSA) was successfully utilized to enhance the robustness of a partial least squares (PLS) model by examining the selectivity of a given component to other components. The samples used in this study were composed of NH₄OH, H₂O₂ and H₂O, a popular etchant solution in the electronic industry. Corresponding near-infrared (NIR) spectra (9000-7500 cm⁻¹) were used to build PLS models. The selective determination of H₂O₂ without influences from NH₄OH and H₂O was a key issue since its molecular structure is similar to that of H₂O and NH₄OH also has a hydroxyl functional group. The best spectral ranges for the determination of NH₄OH and H₂O₂ were found with the use of moving window PLS (MW-PLS) and corresponding selectivity was examined by pure component selectivity analysis. The PLS calibration for NH₄OH was free from interferences from the other components due to the presence of its unique NH absorption bands. Since the spectral variation from H₂O₂ was broadly overlapping and much less distinct than that from NH₄OH, the selectivity and prediction performance for the H₂O₂ calibration were sensitively varied depending on the spectral ranges and number of factors used. PCSA, based on the comparison between regression vectors from PLS and the net analyte signal (NAS), was an effective method to prevent over-fitting of the H₂O₂ calibration. A robust H₂O₂ calibration model with minimal interferences from other components was developed. PCSA should be included as a standard method in PLS calibrations where prediction error only is the usual measure of performance.     

Discovery of potential inhibitors targeting the kinase domain of polynucleotide kinase/phosphatase (PNKP): Homology modeling,virtual screening based on multiple conformations,and molecular dynamics simulation

In recent years, the level of interest has been increased in developing the DNA-repair inhibitors, to enhance the cytotoxic effects in the treatment of cancers. Polynucleotide kinase/phosphatase (PNKP) is a critical human DNA repair enzyme that repairs DNA strand breaks by catalyzing the restoration of 5’-phosphate and 3’-hydroxyl termini that are required for subsequent processing by DNA ligases and polymerases. PNKP is the only protein that repairs the 3′-hydroxyl group and 5′-phosphate group, which depicts PNKP as a potential therapeutic target. Besides, PNKP is the only DNA-repair enzyme that contains the 5′-kinase activity, therefore, targeting this kinase domain would motivate the development of novel PNKP-specific inhibitors. However, there are neither crystal structures of human PNKP nor the kinase inhibitors reported so far. Thus, in this present study, a sequential molecular docking-based virtual screening with multiple PNKP conformations integrating homology modeling, molecular dynamics simulation, and binding free energy calculation was developed to discover novel PNKP kinase inhibitors, and the top-scored molecule was finally submitted to molecular dynamics simulation to reveal the binding mechanism between the inhibitor and PNKP. Taken together, the current study could provide some guidance for the molecular docking based-virtual screening of novel PNKP kinase inhibitors.     

Synthesis,crystal structure and activity evaluation of novel 3,4‐dihydro‐1‐benzoxepin‐5(2H)‐one derivatives as protein–tyrosine kinase (PTK) inhibitors

Four new 3,4‐dihydro‐1‐benzoxepin‐5(2H )‐one derivatives, namely (E )‐4‐(5‐bromo‐2‐hydroxybenzylidene)‐6,8‐dimethoxy‐3,4‐dihydrobenzo[b ]oxepin‐5(2H )‐one, ( 7 ), (E )‐4‐[(E )‐3‐(5‐bromo‐2‐hydroxyphenyl)allylidene]‐6,8‐dimethoxy‐3,4‐dihydrobenzo[b ]oxepin‐5(2H )‐one, ( 8 ), (E )‐4‐(5‐bromo‐2‐hydroxybenzylidene)‐6‐hydroxy‐8‐methoxy‐3,4‐dihydrobenzo[b ]oxepin‐5(2H )‐one, C₁₈H₁₅BrO₅, ( 9 ), and (E )‐4‐[(E )‐3‐(5‐bromo‐2‐hydroxyphenyl)allylidene]‐6‐hydroxy‐8‐methoxy‐3,4‐dihydrobenzo[b ]oxepin‐5(2H )‐one, ( 10 ), have been synthesized and characterized by FT–IR, NMR and MS. The structure of ( 9 ) was confirmed by single‐crystal X‐ray diffraction. Crystal structure analysis shows that molecules of ( 9 ) are connected into a one‐dimensional chain in the [010] direction through classical hydrogen bonds and these chains are further extended into a three‐dimensional network via C—H…O interactions. The inhibitory activities of these compounds against protein–tyrosine kinases (PTKs) show that 6‐hydroxy‐substituted compounds ( 9 ) and ( 10 ) are more effective for inhibiting ErbB1 and ErbB2 than are 6‐methoxy‐substituted compounds ( 7 ) and ( 8 ). This may be because ( 9 ) and ( 10 ) could effectively bind to the active pockets of the protein through intermolecular interactions.