Design,synthesis and evaluation of chalcone Mannich base derivatives as multifunctional agents for the potential treatment of Alzheimer’s disease

A series of chalcone Mannich base derivatives were designed, synthesized and evaluated as multifunctional agents for the treatment of Alzheimer’s disease based on the multi-target directed ligands design strategy. In vitro assays demonstrated that most of the derivatives exerted potent selective inhibitory potency on AChE with good multifunctional properties. Among them, representative compound 7c exhibited moderate inhibitory potency for EeAChE (IC₅₀ = 0.44 μM) and MAO-B inhibition (IC₅₀ = 1.21 μM), good inhibitory effect on self-induced Aβ₁₋₄₂ aggregation (55.0%, at 25 μM), biometal chelating property, moderate antioxidant activity with a value 1.93-fold of Trolox. Moreover, both kinetic analysis of AChE inhibition and molecular modeling study revealed that 7c showed a mixed-type inhibition, binding simultaneously to CAS and PAS of AChE. In addition, 7c also displayed high BBB permeability. These properties indicated 7c may be a promising multifunctional agent for the treatment of AD.     

Insight into the binding modes and inhibition mechanisms of adamantyl-based 1,3-disubstituted urea inhibitors in the active site of the human soluble epoxide hydrolase

Soluble epoxide hydrolase (sEH) is a promising new target for treating hypertension and inflammation. Considerable efforts have been devoted to develop novel inhibitors. In this study, the binding modes and interaction mechanisms of a series of adamantyl-based 1,3-disubstituted urea inhibitors were investigated by molecular docking, molecular dynamics simulations, binding free energy calculations, and binding energy decomposition analysis. Based on binding affinity, the most favorable binding mode was determined for each inhibitor. The calculation results indicate that the total binding free energy (ΔG_TOT, the sum of enthalpy ΔG_MM-GB/SA, and entropy ?TΔS) presents a good correlation with the experimental inhibitory activity (IC₅₀, r²?=?.99). The van der Waals energy contributes most to the total binding free energy (ΔG_TOT). A detailed discussion on the interactions between inhibitors and those residues located in the active pocket is made based on hydrogen bond and binding modes analysis. According to binding energy decomposition, the residues Asp333 and Trp334 contribute the most to binding free energy in all systems. Furthermore, Hip523 plays a major role in determining this class of inhibitor-binding orientations. Combined with the results of hydrogen bond analysis and binding free energy, we believe that the conserved hydrogen bonds play a role only in anchoring the inhibitors to the exact site for binding and the number of hydrogen bonds may not directly relate to the binding free energy. The results we obtained will provide valuable information for the design of high potency sEH inhibitors.     

Structural and chemical basis for enhanced affinity to a series of mycobacterial thymidine monophosphate kinase inhibitors: fragment-based QSAR and QM/MM docking studies

Tuberculosis (TB) still remains one of the most deadly infectious diseases. Mycobacterium tuberculosis thymidine monophosphate kinase (TMPKmt) has emerged as an attractive molecular target for the design of a novel class of anti-TB agents since blocking it will affect the pathways involved in DNA replication. Aiming at shedding some light on structural and chemical features that are important for the affinity of thymidine derivatives to TMPKmt, we have employed a special fragment-based method to develop robust quantitative structure-activity relationship models for a large and chemically diverse series of thymidine-based analogues. Significant statistical parameters (r ² ?=?0.94, q ² ?=?0.76, r ² _pred ?=?0.89) were obtained, indicating the reliability of the hologram QSAR model in predicting the biological activity of untested compounds. The 2D model was then used to predict the potency of an external test set, and the predicted values obtained from the HQSAR model were in good agreement with the experimental results. We have accordingly designed novel TMPKmt inhibitors by utilizing the fragments proposed by HQSAR analysis and predicted with good activity in the developed models. The new designed compounds also presented drug-like characteristics based on Lipinski’s rule of 5. The generated molecular recognition patterns gathered from the HQSAR analysis combined with quantum mechanics/molecular mechanics (QM/MM) docking studies, provided important insights into the chemical and structural basis involved in the molecular recognition process of this series of thymidine analogues and should be useful for the design of new potent anti-TB agents.     

Phosphonic and Phosphinic Acid Derivatives as Novel Tyrosinase Inhibitors: Kinetic Studies and Molecular Docking

A dozen of phosphonic and phosphinic acid derivatives containing pyridine moiety were synthesized and its inhibitory activity toward mushroom tyrosinase was investigated. Moreover, molecular docking of these compounds to the active site of the enzyme was performed. All the compounds ( 1 – 10 ) demonstrated the inhibitory effect with the IC₅₀ and inhibition constants ranging millimolar concentrations. The obtained results indicate that the compounds show different types of inhibition (competitive, noncompetitive, mixed), but all of them are reversible inhibitors. The obtained outcomes allowed to make the structure–activity relationship (SAR) analysis. Compound 4 ([(benzylamino)(pyridin‐2‐yl)methyl]phenylphosphinic acid) revealed the lowest IC₅₀ value of 0.3 mm and inhibitory constant of K_i 0.076 mm , with noncompetitive type and reversible mechanism of inhibition. According to SAR analysis, introducing bulky phenyl moieties to phosphonic and amino groups plays an important role in the inhibitory potency on activity of mushroom tyrosinase and could be useful in design and development of a new class of potent organophosphorus inhibitors of tyrosinase. Combined results of molecular docking and SAR analysis can be helpful in designing novel tyrosinase inhibitors of desired properties. They may have broad application in food industry and cosmetology.     

Molecular modeling studies,synthesis and biological evaluation of dabigatran analogues as thrombin inhibitors

In this work, 48 thrombin inhibitors based on the structural scaffold of dabigatran were analyzed using a combination of molecular modeling techniques. We generated three-dimensional quantitative structure–activity relationship (3D-QSAR) models based on three alignments for both comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA) to highlight the structural requirements for thrombin protein inhibition. In addition to the 3D-QSAR study, Topomer CoMFA model also was established with a higher leave-one-out cross-validation q² and a non-cross-validation r², which suggest that the three models have good predictive ability. The results indicated that the steric, hydrophobic and electrostatic fields play key roles in QSAR model. Furthermore, we employed molecular docking and re-docking simulation explored the binding relationship of the ligand and the receptor protein in detail. Molecular docking simulations identified several key interactions that were also indicated through 3D-QSAR analysis. On the basis of the obtained results, two compounds were designed and predicted by three models, the biological evaluation in vitro (IC₅₀) demonstrated that these molecular models were effective for the development of novel potent thrombin inhibitors.     

CoMFA studies and in vitro evaluation of some 3-substituted benzylthio quinolinium salts as anticryptococcal agents

The 3-dimensional quantitative structure–activity relationship (3D-QSAR) molecular modeling technique or comparative molecular field analysis (CoMFA) has been used to design analogs of the natural product cryptolepine (1). Twenty-three compounds with their in vitro biological activities (IC₅₀ values) against Crytococcus neoformans were used to generate the training set database of compounds for the CoMFA studies. The cross-validated q², noncross-validated r², and partial least squares (PLS) analysis results were used to predict the biological activity of 11 newly designed test set compounds. The best CoMFA model produced a q² of 0.815 and an r² of 0.976 indicating high statistical significance as a predictive model. The steric and electrostatic contributions from the contour map were interpreted from the color-coded contour plots generated from the PLS model and the active structural components for potency against C. neoformans were determined and validated in the test set compounds. The 3-substituted benzylthio quinolinium salts (4) that make up the test set were synthesized and evaluated based on the predicted activity from the CoMFA model and the results produced a good correlation between the predicted and experimental activity (R = 0.82). Thus, CoMFA has served as an effective tool to aid the design of new analogs and in this case, it has aided the identification of compounds equipotent with amphotericin B, the gold standard in antifungal drug design.     

An atypical residue in the pore of Varroa destructor GABA-activated RDL receptors affects picrotoxin block and thymol modulation

GABA-activated RDL receptors are the insect equivalent of mammalian GABA_A receptors, and play a vital role in neurotransmission and insecticide action. Here we clone the pore lining M2 region of the Varroa mite RDL receptor and show that it has 4 atypical residues when compared to M2 regions of most other insects, including bees, which are the major host of Varroa mites. We create mutant Drosophila RDL receptors containing these substitutions and characterise their effects on function. Using two electrode voltage clamp electrophysiology we show that one substitution (T6′M) ablates picrotoxin inhibition and increases the potency of GABA. This mutation also alters the effect of thymol, which enhances both insect and mammalian GABA responses, and is widely used as a miticide. Thymol decreases the GABA EC₅₀ of WT receptors, enhancing responses, but in T6′M-containing receptors it is inhibitory. The other 3 atypical residues have no major effects on either the GABA EC₅₀, the picrotoxin potency or the effect of thymol. In conclusion we show that the RDL 6′ residue is important for channel block, activation and modulation, and understanding its function also has the potential to prove useful in the design of Varroa-specific insecticidal agents.     

Molecular modelling studies on d-annulated benzazepinones as VEGF-R2 kinase inhibitors using docking and 3D-QSAR

Chemotypes comprising the d-annulated 1,3-dihydro-2H-1-benzazepin-2-one scaffold derived from the paullone structure were found to be potent vascular endothelial growth factor receptor 2 (VEGF-R2) kinase inhibitors. Three-dimensional quantitative structure-activity relationship (3D-QSAR) and docking studies were performed on a series of d-annulated benzazepinones with VEGF-R2 kinase inhibition activities. The comparative molecular field analysis and comparative molecular similarity indices analysis models using 32 molecules in the training set gave r²_cv values of 0.811 and 0.769, r² values of 0.962 and 0.953, respectively. 3D contour maps generated from the two models revealed that the electron-withdrawing groups at R₁ and the bulky, electron-withdrawing as well as hydrogen bond donor groups at R₂ position are favourable; the bulky, hydrogen bond acceptor substituent at R₃ and the minor groups at R₄ position may benefit the potency. We have designed a series of novel VEGF-R2 inhibitors by utilizing the SAR results revealed in the present study, which were predicted with excellent potencies in the developed models. The results may aid in designing of potential VEGF-R2 inhibitors with better activities.     

Design,synthesis, and evaluation of 2-piperidone derivatives for the inhibition of β-amyloid aggregation and inflammation mediated neurotoxicity

A series of novel multipotent 2-piperidone derivatives were designed, synthesized and biologically evaluated as chemical agents for the treatment of Alzheimer’s disease (AD). The results showed that most of the target compounds displayed significant potency to inhibit Aβ_1–42 self-aggregation. Among them, compound 7q exhibited the best inhibition of Aβ_1–42 self-aggregation (59.11% at 20 μM) in a concentration-dependent manner. Additionally, the compounds 6b, 7p and 7q as representatives were found to present anti-inflammation properties in lipopolysaccharide (LPS)-induced microglial BV-2 cells. They could effectively suppress the production of pro-inflammatory cytokines such as TNF-α, IL-1β and IL-6. Meanwhile, compound 7q could prevent the neuronal cell SH-SY5Y death by LPS-stimulated microglia cell activation mediated neurotoxicity. The molecular modeling studies demonstrated that compounds matched the pharmacophore well and had good predicted physicochemical properties and estimated IC₅₀ values. Moreover, compound 7q exerted a good binding to the active site of myeloid differentiation factor 88 (MyD88) through the docking analysis and could interfere with its homodimerization or heterodimerization. Consequently, these compounds emerged as promising candidates for further development of novel multifunctional agents for AD treatment.     

Binding conformations and QSAR of CA-4 analogs as tubulin inhibitors

A theoretical study on the binding conformations and the quantitative structure–activity relationship (QSAR) of combretastatin A4 (CA-4) analogs as inhibitors toward tubulin has been carried out using docking analysis and comparative molecular field analysis (CoMFA). The appropriate binding orientations and conformations of these compounds interacting with tubulin were revealed by the docking study; and a 3D-QSAR model showing significant statistical quality and satisfactory predictive ability was established, in which the correlation coefficient (R²) and cross-validation coefficient (q²) were 0.955 and 0.66, respectively. The same model was further applied to predict the pIC₅₀ values for 16 congeneric compounds as external test set, and the predictive correlation coefficient R²_pred reached 0.883. Other tests on additional validations further confirmed the satisfactory predictive power of the model. In this work, it was very interesting to find that the 3D topology structure of the active site of tubulin from the docking analysis was in good agreement with the 3D-QSAR model from CoMFA for this series of compounds. Some key structural factors of the compounds responsible for cytotoxicity were reasonably presented. These theoretical results can offer useful references for understanding the action mechanism and directing the molecular design of this kind of inhibitor with improved activity.     

Molecular docking and 3D-QSAR studies on β-phenylalanine derivatives as dipeptidyl peptidase IV inhibitors

Three-dimensional quantitative structure-activity relationship (3D-QSAR) and molecular docking studies were carried out to explore the binding of 73 inhibitors to dipeptidyl peptidase IV (DPP-IV), and to construct highly predictive 3D-QSAR models using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). The negative logarithm of IC₅₀ (pIC₅₀) was used as the biological activity in the 3D-QSAR study. The CoMFA model was developed by steric and electrostatic field methods, and leave-one-out cross-validated partial least squares analysis yielded a cross-validated value (r_cv² {\hbox{r}}_{{\rm{cv}}}^{\rm{2}} ) of 0.759. Three CoMSIA models developed by different combinations of steric, electrostatic, hydrophobic and hydrogen-bond fields yielded significant r_cv² {\hbox{r}}_{{\rm{cv}}}^{\rm{2}} values of 0.750, 0.708 and 0.694, respectively. The CoMFA and CoMSIA models were validated by a structurally diversified test set of 18 compounds. All of the test compounds were predicted accurately using these models. The mean and standard deviation of prediction errors were within 0.33 and 0.26 for all models. Analysis of CoMFA and CoMSIA contour maps helped identify the structural requirements of inhibitors, with implications for the design of the next generation of DPP-IV inhibitors for the treatment of type 2 diabetes.     

Induced fit docking,pharmacophore modeling,and molecular dynamic simulations on thiazolidinedione derivatives to explore key interactions with Tyr48 in polyol pathway

AT₁ receptor is an interesting biological target involved in several important diseases, such as blood hypertension and cardiovascular pathologies. In this study we investigated the main electrostatic and steric features of a series of AT₁ antagonists related to hypertensive activity using structure and ligand-based strategies (docking and CoMFA). The generated 3D model had good internal and external consistency and was used to predict the potency of an external test set. The predicted values of pIC₅₀ are in good agreement with the experimental results of biological activity, indicating that the 3D model can be used to predict the biological property of untested compounds. The electrostatic and steric CoMFA maps showed molecular recognition patterns, which were analyzed with structure-based molecular modeling studies (docking). The most and the least potent compounds docked into the AT₁ binding site were subjected to molecular dynamics simulations with the aim to verify the stability and the flexibility of the ligand-receptor interactions. These results provided valuable insights on the electronic/structural requirements to design novel AT₁ antagonists.     

Discovery of arylamide-5-anilinoquinazoline-8-nitro derivatives as VEGFR-2 kinase inhibitors: Synthesis,in vitro biological evaluation and molecular docking

Herein, we embarked on a structural optimization campaign aiming at the discovery of novel anticancer agents with our previously reported XL-6f as a lead compound. A library of 23 compounds has been synthesized based on the highly conserved active site of VEGFR-2. Several title compounds exhibited selective inhibitory activities against VEGFR-2, which also displayed selective anti-proliferation potency against HepG2 cell. All synthesized compounds were evaluated for anti-angiogenesis capability. Compound 7o showed the most potent anti-angiogenesis ability, the efficient cytotoxic activities (in vitro against HUVEC and HepG2 cell lines with IC₅₀ values of 0.58 and 0.23 µM, respectively). The molecular docking analysis revealed 7o is a Type-II inhibitor of VEGFR-2 kinase. In general, these results indicated these arylamide-5-anilinoquinazoline-8-nitro derivatives are promising inhibitors of VEGFR-2 for the potential treatment of anti-angiogenesis.     

Design,synthesis and biological evaluation of novel oseltamivir derivatives as potent neuraminidase inhibitors

Neuraminidase (NA) is one of the particular potential targets for novel antiviral therapy. In this work, a series of neuraminidase inhibitors with the cyclohexene scaffold were studied based upon the combination of 3D-QSAR, molecular docking, and molecular dynamics techniques. The results indicate that the built 3D-QSAR models yield reliable statistical information: the correlation coefficient (r²) and cross-validation coefficient (q²) of CoMFA (comparative molecular field analysis) are 0.992 and 0.819; the r² and q² of CoMSIA (comparative molecular similarity analysis) are 0.992 and 0.863, respectively. Molecular docking and MD simulations were conducted to confirm the detailed binding mode of enzyme-inhibitor system. The new NA inhibitors had been designed, synthesized, and their inhibitory activities against group-1 neuraminidase were determined. One agent displayed excellent neuraminidase inhibition, with IC₅₀ value of 39.6?μM against NA, while IC₅₀ value for oseltamivir is 61.1?μM. This compound may be further investigated for the treatment of infection by the new type influenza virus.     

Probing the structural requirements of A-type Aurora kinase inhibitors using 3D-QSAR and molecular docking analysis

Aurora-A, the most widely studied isoform of Aurora kinase overexpressed aberrantly in a wide variety of tumors, has been implicated in early mitotic entry, degradation of natural tumor suppressor p53 and centrosome maturation and separation; hence, potent inhibitors of Aurora-A may be therapeutically useful drugs in the treatment of various forms of cancer. Here, we report an in silico study on a group of 220 reported Aurora-A inhibitors with six different substructures. Three-dimensional quantitative structure–activity relationship (3D-QSAR) studies were carried out using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) techniques on this series of molecules. The resultant optimum 3D-QSAR models exhibited an r _cv² value of 0.404-0.582 and their predictive ability was validated using an independent test set, ending in r _pred² 0.512-0.985. In addition, docking studies were employed to explore these protein–inhibitor interactions at the molecular level. The results of 3D-QSAR and docking analyses validated each other, and the key structural requirements affecting Aurora-A inhibitory activities, and the influential amino acids involved were identified. To the best of our knowledge, this is the first report on 3D-QSAR modeling of Aurora-A inhibitors, and the results can be used to accurately predict the binding affinity of related analogues and also facilitate the rational design of novel inhibitors with more potent biological activities.     

Structure-activity relationship investigation of Phe-Arg mimetic region of human glutaminyl cyclase inhibitors

Glutamyl cyclase (QC) is a promising therapeutic target because of its involvement in the pathogenesis of Alzheimer’s disease. In this study, we developed novel QC inhibitors that contain 3-aminoalkyloxy-4-methoxyphenyl and 4-aminoalkyloxyphenyl groups to replace the previously developed pharmacophore. Several potent inhibitors were identified, showing IC₅₀ values in a low nanomolar range, and were further studied for in vitro toxicity and in vivo activity. Among these, inhibitors 51 and 53 displayed the most potent Aβ_N3pE?40-lowering effects in in vivo acute model with reasonable BBB penetration, without showing cytotoxicity and hERG inhibition. The molecular modeling analysis of 53 indicated that the salt bridge interaction and the hydrogen bonding in the active site provided a high potency. Given the potent activity and favorable BBB penetration with low cytotoxicity, we believe that compound 53 may serve as a potential candidate for anti-Alzheimer’s agents.     

Discovery of highly potent V600E-B-RAF kinase inhibitors: Molecular modeling study

A series of 20 triarylpyrazole derivatives containing amide or urea linker have been synthesized. Their in vitro antiproliferative activity against NCI-60 cancer cell lines panel has been reported. Upon investigating the mechanism of action at molecular level, compound 1e showed selectivity and potency against V600E-B-RAF (IC₅₀?=?390?nM). Herein, we decided to investigate the potency of the other nineteen target compounds against V600E-B-RAF. This led to discovery of several more potent compounds against that kinase. The IC₅₀ values of compounds 1g–i and 2f–i were within the range of 7–47?nM. Among them, the diarylurea compound 1i was the most potent (IC₅₀?=?7?nM). Results of docking and molecular dynamic analysis suggested the presence of consistent binding mode among our compound series with type-IIA class of inhibition pattern. Subsequently, the contribution of structural features to bioactivity were explored by means of QSAR analysis, where such effort led to the development of predictive QSAR model with significant statistical parameters (R²?=?0.912, F?=?38.64, Q²_LOO?=?0.834, Q²_LMO?=?0.816, s?=?0.334). Furthermore, pharmacophoric features existed among our compound series were investigated employing molecular interaction field (MIF) analysis, which led to the development of partial least squares model consisted of four latent variables (4LV-PLS) with statistical parameters of (R²acc.?=?0.98, Q²acc.?=?0.81).