对羟苯基丙酮酸双氧化酶抑制剂的研究进展

本文对对羟苯基丙酮酸双氧化酶抑制剂(HPPD)的发现、作用机制、结构与活性、合成方法及主要新品种加以概述。     

HPPD抑制剂类除草剂作用机制和研究进展

介绍了对羟苯基丙酮酸双氧化酶(HPPD)在生物体内的作用机制以及HPPD抑制剂类除草剂的主要品种和最新研究进展。     

新型乙酰胆碱酯酶抑制剂的设计与虚拟筛选

运用分子对接软件GOLD,研究了延胡索类生物碱corydaline与乙酰胆碱酯酶的结合模型,以coryaaline为先导化合物,设计与虚拟筛选一系列四氢异喹啉衍生物。预测活性做好的化合物1A环与Trp84发生π-π堆积;质子化的氮原子与Phe330发生阳离子-π作用;B环与Trp279发生π-π作用;1位的羟甲基、4位甲氧基、C环上的甲氧基分别与Tyr334、Tyr21、Tyrl30形成氢键作用。     

Identification of New KRAS G12D Inhibitors through Computer-Aided Drug Discovery Methods

Owing to several mutations, the oncogene Kirsten rat sarcoma 2 viral oncogene homolog (KRAS) is activated in the majority of cancers, and targeting it has been pharmacologically challenging. In this study, using an in silico approach comprised of pharmacophore modeling, molecular docking, and molecular dynamics simulations, potential KRAS G12D inhibitors were investigated. A ligand-based common feature pharmacophore model was generated to identify the framework necessary for effective KRAS inhibition. The chemical features in the selected pharmacophore model comprised two hydrogen bond donors, one hydrogen bond acceptor, two aromatic rings and one hydrophobic feature. This model was used for screening in excess of 214,000 compounds from InterBioScreen (IBS) and ZINC databases. Eighteen compounds from the IBS and ten from the ZINC database mapped onto the pharmacophore model and were subjected to molecular docking. Molecular docking results highlighted a higher affinity of four hit compounds towards KRAS G12D in comparison to the reference inhibitor, BI-2852. Sequential molecular dynamics (MD) simulation studies revealed all four hit compounds them possess higher KRAS G12D binding free energy and demonstrate stable polar interaction with key residues. Further, Principal Component Analysis (PCA) analysis of the hit compounds in complex with KRAS G12D also indicated stability. Overall, the research undertaken provides strong support for further in vitro testing of these newly identified KRAS G12D inhibitors, particularly Hit1 and Hit2.     

Structural Determination of Three Different Series of Compounds as Hsp90 Inhibitors Using 3D-QSAR Modeling, Molecular Docking and Molecular Dynamics Methods

Hsp90 is involved in correcting, folding, maturation and activation of a diverse array of client proteins; it has also been implicated in the treatment of cancer in recent years. In this work, comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA), molecular docking and molecular dynamics were performed on three different series of Hsp90 inhibitors to build 3D-QSAR models, which were based on the ligand-based or receptor-based methods. The optimum 3D-QSAR models exhibited reasonable statistical characteristics with averaging internal q(2) > 0.60 and external r(2) (pred) > 0.66 for Benzamide tetrahydro-4H-carbazol-4-one analogs (BT), AT13387 derivatives (AT) and Dihydroxylphenyl amides (DA). The results revealed that steric effects contributed the most to the BT model, whereas H-bonding was more important to AT, and electrostatic, hydrophobic, H-bond donor almost contributed equally to the DA model. The docking analysis showed that Asp93, Tyr139 and Thr184 in Hsp90 are important for the three series of inhibitors. Molecular dynamics simulation (MD) further indicated that the conformation derived from docking is basically consistent with the average structure extracted from MD simulation. These results not only lead to a better understanding of interactions between these inhibitors and Hsp90 receptor but also provide useful information for the design of new inhibitors with a specific activity.     

Novel Monocyclam Derivatives as HIV Entry Inhibitors: Design,Synthesis, Anti‐HIV Evaluation,and Their Interaction with the CXCR4 Co‐receptor

The CXCR4 receptor has been shown to interact with the human immunodeficiency virus (HIV) envelope glycoprotein gp120, leading to fusion of viral and cell membranes. Therefore, ligands that can attach to this receptor represent an important class of therapeutic agents against HIV, thus inhibiting the first step in the cycle of viral infection: the virus–cell entry/fusion. Herein we describe the in silico design, synthesis, and biological evaluation of novel monocyclam derivatives as HIV entry inhibitors. In vitro activity testing of these compounds in cell cultures against HIV strains revealed EC₅₀ values in the low micromolar range without cytotoxicity at the concentrations tested. Docking and molecular dynamics simulations were performed to predict the binding interactions between CXCR4 and the novel monocyclam derivatives. A binding mode of these compounds is proposed which is consistent with the main existing site‐directed mutagenesis data on the CXCR4 co‐receptor. Moreover, molecular modeling comparisons were performed between these novel monocyclams, previously reported non‐cyclam compounds from which the monocyclams are derived, and the well‐known AMD3100 bicyclam CXCR4 inhibitors. Our results suggest that these three structurally diverse CXCR4 inhibitors bind to overlapping but not identical amino acid residues in the transmembrane regions of the receptor.     

Exploration of N-Arylsulfonyl-indole-2-carboxamide Derivatives as Novel Fructose-1,6-bisphosphatase Inhibitors by Molecular Simulation

A series of N-arylsulfonyl-indole-2-carboxamide derivatives have been identified as potent fructose-1,6-bisphosphatase (FBPase) inhibitors (FBPIs) with excellent selectivity for the potential therapy of type II diabetes mellitus. To explore the structure–activity relationships (SARs) and the mechanisms of action of these FBPIs, a systematic computational study was performed in the present study, including three-dimensional quantitative structure–activity relationship (3D-QSAR) modeling, pharmacophore modeling, molecular dynamics (MD), and virtual screening. The constructed 3D-QSAR models exhibited good predictive ability with reasonable parameters using comparative molecular field analysis (q² = 0.709, R² = 0.979, r_pre² = 0.932) and comparative molecular similarity indices analysis (q² = 0.716, R² = 0.978, r_pre² = 0.890). Twelve hit compounds were obtained by virtual screening using the best pharmacophore model in combination with molecular dockings. Three compounds with relatively higher docking scores and better ADME properties were then selected for further studies by docking and MD analyses. The docking results revealed that the amino acid residues Met18, Gly21, Gly26, Leu30, and Thr31 at the binding site were of great importance for the effective bindings of these FBPIs. The MD results indicated that the screened compounds VS01 and VS02 could bind with FBPase stably as its cognate ligand in dynamic conditions. This work identified several potential FBPIs by modeling studies and might provide important insights into developing novel FBPIs.     

新型乙酰胆碱酯酶抑制剂的分子对接研究

构建了乙酰胆碱酯酶抑制剂药效团模型,并基于此模型进行新型乙酰胆碱酯酶抑制剂的分子设计。用MVD软件进行分子对接,预测设计分子与乙酰胆碱酯酶的结合活性。通过计算发现,基于此药效团设计的小分子9与14乙酰胆碱酯酶具有较好的结合活性。     

A Computational Workflow for the Identification of Novel Fragments Acting as Inhibitors of the Activity of Protein Kinase CK1δ

Fragment-Based Drug Discovery (FBDD) has become, in recent years, a consolidated approach in the drug discovery process, leading to several drug candidates under investigation in clinical trials and some approved drugs. Among these successful applications of the FBDD approach, kinases represent a class of targets where this strategy has demonstrated its real potential with the approved kinase inhibitor Vemurafenib. In the Kinase family, protein kinase CK1 isoform δ (CK1δ) has become a promising target in the treatment of different neurodegenerative diseases such as Alzheimer’s disease, Parkinson’s disease, and amyotrophic lateral sclerosis. In the present work, we set up and applied a computational workflow for the identification of putative fragment binders in large virtual databases. To validate the method, the selected compounds were tested in vitro to assess the CK1δ inhibition.     

Pharmacophore and Molecular Docking Guided 3D-QSAR Study of Bacterial Enoyl-ACP Reductase (FabI) Inhibitors

Enoyl acyl carrier protein (ACP) reductase (FabI) is a potential target for the development of antibacterial agents. Three-dimensional quantitative structure-activity relationships (3D-QSAR) for substituted formamides series of FabI inhibitors were investigated using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) techniques. Pharmacophore and molecular docking methods were used for construction of the molecular alignments. A training set of 36 compounds was performed to create the 3D-QSAR models and their external predictivity was proven using a test set of 11 compounds. Graphical interpretation of the results revealed important structural features of the formamides related to the active site of FabI. The results may be exploited for further optimization of the design of new potent FabI inhibitors.     

AKT Inhibitors: The Road Ahead to Computational Modeling-Guided Discovery

AKT, is a serine/threonine protein kinase comprising three isoforms—namely: AKT1, AKT2 and AKT3, whose inhibitors have been recognized as promising therapeutic targets for various human disorders, especially cancer. In this work, we report a systematic evaluation of multi-target Quantitative Structure-Activity Relationship (mt-QSAR) models to probe AKT’ inhibitory activity, based on different feature selection algorithms and machine learning tools. The best predictive linear and non-linear mt-QSAR models were found by the genetic algorithm-based linear discriminant analysis (GA-LDA) and gradient boosting (Xgboost) techniques, respectively, using a dataset containing 5523 inhibitors of the AKT isoforms assayed under various experimental conditions. The linear model highlighted the key structural attributes responsible for higher inhibitory activity whereas the non-linear model displayed an overall accuracy higher than 90%. Both these predictive models, generated through internal and external validation methods, were then used for screening the Asinex kinase inhibitor library to identify the most potential virtual hits as pan-AKT inhibitors. The virtual hits identified were then filtered by stepwise analyses based on reverse pharmacophore-mapping based prediction. Finally, results of molecular dynamics simulations were used to estimate the theoretical binding affinity of the selected virtual hits towards the three isoforms of enzyme AKT. Our computational findings thus provide important guidelines to facilitate the discovery of novel AKT inhibitors.     

A Computational Study on Thiourea Analogs as Potent MK-2 Inhibitors

Mitogen-activated protein kinase-activated protein kinase 2 (MK-2) has been identified as a drug target for the treatment of inflammatory diseases. Currently, a series of thiourea analogs as potent MK-2 inhibitors were studied using comprehensive computational methods by 3D-QSAR, molecular docking and molecular dynamics simulations for a further improvement in activities. The optimal 3D models exhibit high statistical significance of the results, especially for the CoMFA results with r(2) (ncv), q(2) values of 0.974, 0.536 for the internal validation, and r(2) (pred), r(2) (m) values of 0.910, 0.723 for the external validation and Roy's index, respectively. In addition, more rigorous validation criteria suggested by Tropsha were also employed to check the built models. Graphic representation of the results, as contoured 3D coefficient plots, also provides a clue to the reasonable modification of molecules: (i) The substituent with a bulky size and electron-rich group at the C5 position of the pyrazine ring is required to enhance the potency; (ii) The H-bond acceptor group in the C3 position of the pyrazine ring is likely to be helpful to increase MK-2 inhibition; (iii) The small and electropositive substituent as a hydrogen bond donor of the C2 position in the oxazolone ring is favored; In addition, several important amino acid residues were also identified as playing an important role in MK-2 inhibition. The agreement between 3D-QSAR, molecular docking and molecular dynamics simulations also proves the rationality of the developed models. These results, we hope, may be helpful in designing novel and potential MK-2 inhibitors.     

Rational Design,Synthesis, and DNA Binding Properties of Novel Sequence‐Selective Peptidyl Congeners of Ametantrone

Natural and synthetic compounds characterized by an anthraquinone nucleus represent an important class of anti‐neoplastic agents, the mechanism of action of which is related to intercalation into DNA. Ametantrone (AM) is a synthetic 9,10‐anthracenedione bearing two (hydroxyethylamino)ethylamino residues at positions 1 and 4; along with other anthraquinones and anthracyclines, it shares a polycyclic intercalating moiety and charged side chains that stabilize DNA binding. All these drugs elicit adverse side effects, which represent a challenge for antitumor chemotherapy. In the present work the structure of AM was augmented with appropriate groups that target well‐defined base pairs in the major groove. These should endow AM with DNA sequence selectivity. We describe the rationale for the synthesis and the evaluation of activity of a new series of compounds in which the planar anthraquinone is conjugated at positions 1 and 4 through the side chains of AM or other bioisosteric linkers to appropriate dipeptides. The designed novel AM derivatives were shown to selectively stabilize two oligonucleotide duplexes that both have a palindromic GC‐rich hexanucleotide core, but their stabilizing effects on a random DNA sequence was negligible. In the case of the most effective compound, the 1,4‐bis‐[Gly‐(L ‐Lys)] derivative of AM, the experimental results confirm the predictions of earlier theoretical computations. In contrast, AM had equal stabilizing effects on all three sequences and showed no preferential binding. This novel peptide derivative can be classified as a strong binder regarding the sequences that it selectively targets, possibly opening the exploitation of less cytotoxic conjugates of AM to the targeted treatment of oncological and viral diseases.     

Generation of Non-Nucleotide CD73 Inhibitors Using a Molecular Docking and 3D-QSAR Approach

Radiotherapy and chemotherapy are conventional cancer treatments. Around 60% of all patients who are diagnosed with cancer receive radio- or chemotherapy in combination with surgery during their disease. Only a few patients respond to the blockage of immune checkpoints alone, or in combination therapy, because their tumours might not be immunogenic. Under these circumstances, an increasing level of extracellular adenosine via the activation of ecto-5’-nucleotidase (CD73) and consequent adenosine receptor signalling is a typical mechanism that tumours use to evade immune surveillance. CD73 is responsible for the conversion of adenosine monophosphate to adenosine. CD73 is overexpressed in various tumour types. Hence, targetting CD73’s signalling is important for the reversal of adenosine-facilitated immune suppression. In this study, we selected a potent series of the non-nucleotide small molecule inhibitors of CD73. Molecular docking studies were performed in order to examine the binding mode of the inhibitors inside the active site of CD73 and 3D-QSAR was used to study the structure–activity relationship. The obtained CoMFA (q² = 0.844, ONC = 5, r² = 0.947) and CoMSIA (q² = 0.804, ONC = 4, r² = 0.954) models showed reasonable statistical values. The 3D-QSAR contour map analysis revealed useful structural characteristics that were needed to modify non-nucleotide small molecule inhibitors. We used the structural information from the overall docking and 3D-QSAR results to design new, potent CD73 non-nucleotide inhibitors. The newly designed CD73 inhibitors exhibited higher activity (predicted pIC₅₀) than the most active compound of all of the derivatives that were selected for this study. Further experimental studies are needed in order to validate the new CD73 inhibitors.     

3D-QSAR Studies on Barbituric Acid Derivatives as Urease Inhibitors and the Effect of Charges on the Quality of a Model

Urease enzyme (EC 3.5.1.5) has been determined as a virulence factor in pathogenic microorganisms that are accountable for the development of different diseases in humans and animals. In continuance of our earlier study on the helicobacter pylori urease inhibition by barbituric acid derivatives, 3D-QSAR (three dimensional quantitative structural activity relationship) advance studies were performed by Comparative Molecular Field Analysis (CoMFA) and Comparative Molecular Similarity Indices Analysis (CoMSIA) methods. Different partial charges were calculated to examine their consequences on the predictive ability of the developed models. The finest developed model for CoMFA and CoMSIA were achieved by using MMFF94 charges. The developed CoMFA model gives significant results with cross-validation (q²) value of 0.597 and correlation coefficients (r²) of 0.897. Moreover, five different fields i.e., steric, electrostatic, and hydrophobic, H-bond acceptor and H-bond donors were used to produce a CoMSIA model, with q² and r² of 0.602 and 0.98, respectively. The generated models were further validated by using an external test set. Both models display good predictive power with r²_pred ≥ 0.8. The analysis of obtained CoMFA and CoMSIA contour maps provided detailed insight for the promising modification of the barbituric acid derivatives with an enhanced biological activity.     

Structural Determinants of CX-4945 Derivatives as Protein Kinase CK2 Inhibitors: A Computational Study

Protein kinase CK2, also known as casein kinase-2, is involved in a broad range of physiological events including cell growth, proliferation and suppression of apoptosis which are related to human cancers. A series of compounds were identified as CK2 inhibitors and their inhibitory activities varied depending on their structures. In order to explore the structure-activity correlation of CX-4945 derivatives as inhibitors of CK2, in the present study, a set of ligand- and receptor-based 3D-QSAR models were developed employing Comparative Molecular Field Analysis (CoMFA) and Comparative Molecular Similarity Index Analysis (CoMSIA). The optimum CoMFA (R(cv) (2) = 0.618, R(pred) (2) = 0.892) and CoMSIA (R(cv) (2) = 0.681, R(pred) (2) = 0.843) models exhibited reasonable statistical characteristics for CX-4945 derivatives. The results indicated that electrostatic effects contributed the most to both CoMFA and CoMSIA models. The combination of docking analysis and molecular dynamics (MD) simulation showed that Leu45, Lys68, Glu81, Val116, Asp175 and Trp176 of CK2 which formed several direct or water-bridged H-bonds with CX-4945 are crucial for CX-4945 derivatives recognition to CK2. These results can offer useful theoretical references for designing more potent CK2 inhibitors.     

Tetrahydroindoles as Multipurpose Screening Compounds and Novel Sirtuin Inhibitors

Indoles are privileged structures in medicinal and bioorganic chemistry that are particularly well suited to serve as platforms for diversity. Among many other therapeutic areas, the indole scaffold has been used to design aromatic compounds useful to interfere with enzymes engaged in the regulation of substrate acylation status, such as sirtuins. However, the planarity of the indole ring is not necessarily optimal for all target enzymes, especially when functionalization with aromatic side chains is required. Replacement of flat scaffolds by nonplanar molecular cores dominated by sp³ hybridization is a common strategy to avoid the disadvantages associated with poor solubility and high promiscuity, while covering less-well-explored areas of chemical space. Thus, we synthesized fragment-like tetrahydroindoles suitable for fragment-based drug discovery as well as a well-characterized small library intended as multipurpose screening compounds. For proof of principle, these compounds were screened against sirtuins 1–3, enzymes known to be addressable by indoles. We found that 2,6,6-trimethyl-4-oxo-4,5,6,7-tetrahydro-1H-indole-3-carboxamides are potent and selective SIRT2 inhibitors. Compound 16 t displayed an IC₅₀ value of 0.98 μm and could serve as exquisite starting point for hit-to-lead profiling.     

Ligand-Based Virtual Screening and Molecular Docking of Benzimidazoles as Potential Inhibitors of Triosephosphate Isomerase Identified New Trypanocidal Agents

Trypanosoma cruzi (T. cruzi) is a parasite that affects humans and other mammals. T. cruzi depends on glycolysis as a source of adenosine triphosphate (ATP) supply, and triosephosphate isomerase (TIM) plays a key role in this metabolic pathway. This enzyme is an attractive target for the design of new trypanocidal drugs. In this study, a ligand-based virtual screening (LBVS) from the ZINC15 database using benzimidazole as a scaffold was accomplished. Later, a molecular docking on the interface of T. cruzi TIM (TcTIM) was performed and the compounds were grouped by interaction profiles. Subsequently, a selection of compounds was made based on cost and availability for in vitro evaluation against blood trypomastigotes. Finally, the compounds were analyzed by molecular dynamics simulation, and physicochemical and pharmacokinetic properties were determined using SwissADME software. A total of 1604 molecules were obtained as potential TcTIM inhibitors. BP2 and BP5 showed trypanocidal activity with half-maximal lytic concentration (LC₅₀) values of 155.86 and 226.30 µM, respectively. Molecular docking and molecular dynamics simulation analyzes showed a favorable docking score of BP5 compound on TcTIM. Additionally, BP5 showed a low docking score (−5.9 Kcal/mol) on human TIM compared to the control ligand (−7.2 Kcal/mol). Both compounds BP2 and BP5 showed good physicochemical and pharmacokinetic properties as new anti-T. cruzi agents.     

Discovery of Novel HIV Protease Inhibitors Using Modern Computational Techniques

The human immunodeficiency virus type 1 (HIV-1) has continued to be a global concern. With the new HIV incidence, the emergence of multi-drug resistance and the untoward side effects of currently used anti-HIV drugs, there is an urgent need to discover more efficient anti-HIV drugs. Modern computational tools have played vital roles in facilitating the drug discovery process. This research focuses on a pharmacophore-based similarity search to screen 111,566,735 unique compounds in the PubChem database to discover novel HIV-1 protease inhibitors (PIs). We used an in silico approach involving a 3D-similarity search, physicochemical and ADMET evaluations, HIV protease-inhibitor prediction (IC₅₀/percent inhibition), rigid receptor–molecular docking studies, binding free energy calculations and molecular dynamics (MD) simulations. The 10 FDA-approved HIV PIs (saquinavir, lopinavir, ritonavir, amprenavir, fosamprenavir, atazanavir, nelfinavir, darunavir, tipranavir and indinavir) were used as reference. The in silico analysis revealed that fourteen out of the twenty-eight selected optimized hit molecules were within the acceptable range of all the parameters investigated. The hit molecules demonstrated significant binding affinity to the HIV protease (PR) when compared to the reference drugs. The important amino acid residues involved in hydrogen bonding and п-п stacked interactions include ASP25, GLY27, ASP29, ASP30 and ILE50. These interactions help to stabilize the optimized hit molecules in the active binding site of the HIV-1 PR (PDB ID: 2Q5K). HPS/002 and HPS/004 have been found to be most promising in terms of IC₅₀/percent inhibition (90.15%) of HIV-1 PR, in addition to their drug metabolism and safety profile. These hit candidates should be investigated further as possible HIV-1 PIs with improved efficacy and low toxicity through in vitro experiments and clinical trial investigations.     

Combined 3D-QSAR, Molecular Docking and Molecular Dynamics Study on Derivatives of Peptide Epoxyketone and Tyropeptin-Boronic Acid as Inhibitors Against the β5 Subunit of Human 20S Proteasome

An abnormal ubiquitin-proteasome is found in many human diseases, especially in cancer, and has received extensive attention as a promising therapeutic target in recent years. In this work, several in silico models have been built with two classes of proteasome inhibitors (PIs) by using 3D-QSAR, homology modeling, molecular docking and molecular dynamics (MD) simulations. The study resulted in two types of satisfactory 3D-QSAR models, i.e., the CoMFA model (Q(2) = 0.462, R(2) (pred) = 0.820) for epoxyketone inhibitors (EPK) and the CoMSIA model (Q(2) = 0.622, R(2) (pred) = 0.821) for tyropeptin-boronic acid derivatives (TBA). From the contour maps, some key structural factors responsible for the activity of these two series of PIs are revealed. For EPK inhibitors, the N-cap part should have higher electropositivity; a large substituent such as a benzene ring is favored at the C6-position. In terms of TBA inhibitors, hydrophobic substituents with a larger size anisole group are preferential at the C8-position; higher electropositive substituents like a naphthalene group at the C3-position can enhance the activity of the drug by providing hydrogen bond interaction with the protein target. Molecular docking disclosed that residues Thr60, Thr80, Gly106 and Ser189 play a pivotal role in maintaining the drug-target interactions, which are consistent with the contour maps. MD simulations further indicated that the binding modes of each conformation derived from docking is stable and in accord with the corresponding structure extracted from MD simulation overall. These results can offer useful theoretical references for designing more potent PIs.