Explore the vast range of titles available.

To better understand the mechanism of in vivo toxicity of N-nitroso compounds (NNCs), the toxicity data of 80 NNCs related to their rat acute oral toxicity data (50% lethal dose concentration, LD50) were used to establish quantitative structure-activity relationship (QSAR) and classification models. Quantum chemistry methods calculated descriptors and Dragon descriptors were combined to describe the molecular information of all compounds. Genetic algorithm (GA) and multiple linear regression (MLR) analyses were combined to develop QSAR models. Fingerprints and machine learning methods were used to establish classification models. The quality and predictive performance of all established models were evaluated by internal and external validation techniques. The best GA-MLR-based QSAR model containing eight molecular descriptors was obtained with Q2loo = 0.7533, R2 = 0.8071, Q2ext = 0.7041 and R2ext = 0.7195. The results derived from QSAR studies showed that the acute oral toxicity of NNCs mainly depends on three factors, namely, the polarizability, the ionization potential (IP) and the presence/absence and frequency of C–O bond. For classification studies, the best model was obtained using the MACCS keys fingerprint combined with artificial neural network (ANN) algorithm. The classification models suggested that several representative substructures, including nitrile, hetero N nonbasic, alkylchloride and amine-containing fragments are main contributors for the high toxicity of NNCs. Overall, the developed QSAR and classification models of the rat acute oral toxicity of NNCs showed satisfying predictive abilities. The results provide an insight into the understanding of the toxicity mechanism of NNCs in vivo, which might be used for a preliminary assessment of NNCs toxicity to mammals.

O6-methylguanine-DNA methyltransferase (MGMT), a unique DNA repair enzyme, can confer resistance to DNA anticancer alkylating agents that modify the O6-position of guanine. Thus, inhibition of MGMT activity in tumors has a great interest for cancer researchers because it can significantly improve the anticancer efficacy of such alkylating agents. In this study, we performed a quantitative structure activity relationship (QSAR) and classification study based on a total of 134 base analogs related to their ED50 values (50% inhibitory concentration) against MGMT. Molecular information of all compounds were described by quantum chemical descriptors and Dragon descriptors. Genetic algorithm (GA) and multiple linear regression (MLR) analysis were combined to develop QSAR models. Classification models were generated by seven machine-learning methods based on six types of molecular fingerprints. Performances of all developed models were assessed by internal and external validation techniques. The best QSAR model was obtained with Q2Loo = 0.83, R2 = 0.87, Q2ext = 0.67, and R2ext = 0.69 based on 84 compounds. The results from QSAR studies indicated topological charge indices, polarizability, ionization potential (IP), and number of primary aromatic amines are main contributors for MGMT inhibition of base analogs. For classification studies, the accuracies of 10-fold cross-validation ranged from 0.750 to 0.885 for top ten models. The range of accuracy for the external test set ranged from 0.800 to 0.880 except for PubChem-Tree model, suggesting a satisfactory predictive ability. Three models (Ext-SVM, Ext-Tree and Graph-RF) showed high and reliable predictive accuracy for both training and external test sets. In addition, several representative substructures for characterizing MGMT inhibitors were identified by information gain and substructure frequency analysis method. Our studies might be useful for further study to design and rapidly identify potential MGMT inhibitors.

Dual-specific tyrosine phosphorylation regulated kinase 1 (DYRK1A) has been regarded as a potential therapeutic target of neurodegenerative diseases, and considerable progress has been made in the discovery of DYRK1A inhibitors. Identification of pharmacophoric fragments provides valuable information for structure- and fragment-based design of potent and selective DYRK1A inhibitors. In this study, seven machine learning methods along with five molecular fingerprints were employed to develop qualitative classification models of DYRK1A inhibitors, which were evaluated by cross-validation, test set, and external validation set with four performance indicators of predictive classification accuracy (CA), the area under receiver operating characteristic (AUC), Matthews correlation coefficient (MCC), and balanced accuracy (BA). The PubChem fingerprint-support vector machine model (CA = 0.909, AUC = 0.933, MCC = 0.717, BA = 0.855) and PubChem fingerprint along with the artificial neural model (CA = 0.862, AUC = 0.911, MCC = 0.705, BA = 0.870) were considered as the optimal modes for training set and test set, respectively. A hybrid data balancing method SMOTETL, a combination of synthetic minority over-sampling technique (SMOTE) and Tomek link (TL) algorithms, was applied to explore the impact of balanced learning on the performance of models. Based on the frequency analysis and information gain, pharmacophoric fragments related to DYRK1A inhibition were also identified. All the results will provide theoretical supports and clues for the screening and design of novel DYRK1A inhibitors.

O6-methylguanine-DNA methyltransferase (MGMT) constitutes an important cellular mechanism for repairing potentially cytotoxic DNA damage induced by guanine O6-alkylating agents and can render cells highly resistant to certain cancer chemotherapeutic drugs. A wide variety of potential MGMT inactivators have been designed and synthesized for the purpose of overcoming MGMT-mediated tumor resistance. We determined the inactivation potency of these compounds against human recombinant MGMT using [3H]-methylated-DNA-based MGMT inactivation assays and calculated the IC50 values. Using the results of 370 compounds, we performed quantitative structure–activity relationship (QSAR) modeling to identify the correlation between the chemical structure and MGMT-inactivating ability. Modeling was based on subdividing the sorted pIC50 values or on chemical structures or was random. A total of nine molecular descriptors were presented in the model equation, in which the mechanistic interpretation indicated that the status of nitrogen atoms, aliphatic primary amino groups, the presence of O-S at topological distance 3, the presence of Al-O-Ar/Ar-O-Ar/R..O..R/R-O-C=X, the ionization potential and hydrogen bond donors are the main factors responsible for inactivation ability. The final model was of high internal robustness, goodness of fit and prediction ability (R2pr = 0.7474, Q2Fn = 0.7375–0.7437, CCCpr = 0.8530). After the best splitting model was decided, we established the full model based on the entire set of compounds using the same descriptor combination. We also used a similarity-based read-across technique to further improve the external predictive ability of the model (R2pr = 0.7528, Q2Fn = 0.7387–0.7449, CCCpr = 0.8560). The prediction quality of 66 true external compounds was checked using the “Prediction Reliability Indicator” tool. In summary, we defined key structural features associated with MGMT inactivation, thus allowing for the design of MGMT inactivators that might improve clinical outcomes in cancer treatment.

DNA repair enzyme O6-methylguanine-DNA methyltransferase (MGMT), which plays an important role in inducing drug resistance against alkylating agents that modify the O6 position of guanine in DNA, is an attractive target for anti-tumor chemotherapy. A series of MGMT inhibitors have been synthesized over the past decades to improve the chemotherapeutic effects of O6-alkylating agents. In the present study, we performed a three-dimensional quantitative structure activity relationship (3D-QSAR) study on 97 guanine derivatives as MGMT inhibitors using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) methods. Three different alignment methods (ligand-based, DFT optimization-based and docking-based alignment) were employed to develop reliable 3D-QSAR models. Statistical parameters derived from the models using the above three alignment methods showed that the ligand-based CoMFA (Qcv2 = 0.672 and Rncv2 = 0.997) and CoMSIA (Qcv2 = 0.703 and Rncv2 = 0.946) models were better than the other two alignment methods-based CoMFA and CoMSIA models. The two ligand-based models were further confirmed by an external test-set validation and a Y-randomization examination. The ligand-based CoMFA model (Qext2 = 0.691, Rpred2 = 0.738 and slope k = 0.91) was observed with acceptable external test-set validation values rather than the CoMSIA model (Qext2 = 0.307, Rpred2 = 0.4 and slope k = 0.719). Docking studies were carried out to predict the binding modes of the inhibitors with MGMT. The results indicated that the obtained binding interactions were consistent with the 3D contour maps. Overall, the combined results of the 3D-QSAR and the docking obtained in this study provide an insight into the understanding of the interactions between guanine derivatives and MGMT protein, which will assist in designing novel MGMT inhibitors with desired activity.

Tumor formation and growth depend on various biological metabolism processes that are distinctly different with normal tissues. Abnormal energy metabolism is one of the typical characteristics of tumors. It has been proven that most tumor cells highly rely on aerobic glycolysis to obtain energy rather than mitochondrial oxidative phosphorylation (OXPHOS) even in the presence of oxygen, a phenomenon called “Warburg effect”. Thus, inhibition of aerobic glycolysis becomes an attractive strategy to specifically kill tumor cells, while normal cells remain unaffected. In recent years, a small molecule alkylating agent, 3-bromopyruvate (3-BrPA), being an effective glycolytic inhibitor, has shown great potential as a promising antitumor drug. Not only it targets glycolysis process, but also inhibits mitochondrial OXPHOS in tumor cells. Excellent antitumor effects of 3-BrPA were observed in cultured cells and tumor-bearing animal models. In this review, we described the energy metabolic pathways of tumor cells, mechanism of action and cellular targets of 3-BrPA, antitumor effects, and the underlying mechanism of 3-BrPA alone or in combination with other antitumor drugs (e.g., cisplatin, doxorubicin, daunorubicin, 5-fluorouracil, etc.) in vitro and in vivo. In addition, few human case studies of 3-BrPA were also involved. Finally, the novel chemotherapeutic strategies of 3-BrPA, including wafer, liposomal nanoparticle, aerosol, and conjugate formulations, were also discussed for future clinical application.

All areas of the modern society are affected by fluorine chemistry. In particular, fluorine plays an important role in medical, pharmaceutical and agrochemical sciences. Amongst various fluoro-organic compounds, trifluoromethyl (CF3) group is valuable in applications such as pharmaceuticals, agrochemicals and industrial chemicals. In the present study, following the strict OECD modelling principles, a quantitative structure–toxicity relationship (QSTR) modelling for the rat acute oral toxicity of trifluoromethyl compounds (TFMs) was established by genetic algorithm-multiple linear regression (GA-MLR) approach. All developed models were evaluated by various state-of-the-art validation metrics and the OECD principles. The best QSTR model included nine easily interpretable 2D molecular descriptors with clear physical and chemical significance. The mechanistic interpretation showed that the atom-type electro-topological state indices, molecular connectivity, ionization potential, lipophilicity and some autocorrelation coefficients are the main factors contributing to the acute oral toxicity of TFMs against rats. To validate that the selected 2D descriptors can effectively characterize the toxicity, we performed the chemical read-across analysis. We also compared the best QSTR model with public OPERA tool to demonstrate the reliability of the predictions. To further improve the prediction range of the QSTR model, we performed the consensus modelling. Finally, the optimum QSTR model was utilized to predict a true external set containing many untested/unknown TFMs for the first time. Overall, the developed model contributes to a more comprehensive safety assessment approach for novel CF3-containing pharmaceuticals or chemicals, reducing unnecessary chemical synthesis whilst saving the development cost of new drugs.

N’-nitrosonornicotine (NNN) is one of the tobacco-specific nitrosamines (TSNAs) that exists widely in smoke and smokeless tobacco products. NNN can induce tumors in various laboratory animal models and has been identified by International Agency for Research on Cancer (IARC) as a human carcinogen. Metabolic activation of NNN is primarily initiated by cytochrome P450 enzymes (CYP450s) via 2′-hydroxylation or 5′-hydroxylation. Subsequently, the hydroxylating intermediates undergo spontaneous decomposition to generate diazohydroxides, which can be further converted to alkyldiazonium ions, followed by attacking DNA to form various DNA damages, such as pyridyloxobutyl (POB)-DNA adducts and pyridyl-N-pyrrolidinyl (py-py)-DNA adducts. If not repaired correctly, these lesions would lead to tumor formation. In the present study, we performed density functional theory (DFT) computations and molecular docking studies to understand the mechanism of metabolic activation and carcinogenesis of NNN. DFT calculations were performed to explore the 2′- or 5′- hydroxylation reaction of (R)-NNN and (S)-NNN. The results indicated that NNN catalyzed by the ferric porphyrin (Compound I, Cpd I) at the active center of CYP450 included two steps, hydrogen abstraction and rebound reactions. The free energy barriers of the 2′- and 5′-hydroxylation of NNN are 9.82/8.44 kcal/mol (R/S) and 7.99/9.19 kcal/mol (R/S), respectively, suggesting that the 2′-(S) and 5′-(R) pathways have a slight advantage. The free energy barriers of the decomposition occurred at the 2′-position and 5′-position of NNN are 18.04/18.02 kcal/mol (R/S) and 18.33/19.53 kcal/mol (R/S), respectively. Moreover, we calculated the alkylation reactions occurred at ten DNA base sites induced by the 2′-hydroxylation product of NNN, generating the free energy barriers ranging from 0.86 to 4.72 kcal/mol, which indicated that these reactions occurred easily. The docking study showed that (S)-NNN had better affinity with CYP450s than that of (R)-NNN, which was consistent with the experimental results. Overall, the combined results of the DFT calculations and the docking obtained in this study provide an insight into the understanding of the carcinogenesis of NNN and other TSNAs.

Protein kinase CK2 is a validated target for cancer therapy. Many natural products have shown inhibitory activity against CK2 as potential anti-cancer drug candidates. A compatible quantitative structure-activity relationship (QSAR) model of natural products is necessary to identify the structural determinants related to their biological activities and provides valuable clues for the discovery of natural leads as anticancer drugs. In this study, genetic algorithm (GA) and multiple linear regression (MLR) methods, combined with preferred molecular descriptors, were employed to build QSAR models of CK2 natural product inhibitors. The best model, composed of eight molecular descriptors, yielded Q 2 Loo  = 0.7914 and R 2  = 0.8220 for the training set and Q 2 ext  = 0.7921 and R 2 ext  = 0.7998 for the test set, indicating the model’s robust reliability and high predictability. As a proof of concept, a true external test set, distinct from the training and test sets, was synthesized and tested in vitro to verify the predictive ability of this model. The predicted pIC 50 values of 13 compounds showed less than 30% relative error (including 10 compounds with relative errors less than 20%), further validating the predictive performance of this model. And compound M18, M24, and M26 were identified as potential CK2 inhibitors with the predicted pIC 50 values of 11.29, 8.79, and 12.03 respectively. Furthermore, the underlying structural mechanisms through which key molecular descriptors influenced their inhibitory activities against CK2 were elucidated. All these results provide valuable information for the discovery of CK2 inhibitors. Graphical Abstract

O6-alkylguanine-DNA alkyltransferase (AGT) is the main cause of tumor cell resistance to DNA-alkylating agents, so it is valuable to design tumor-targeted AGT inhibitors with hypoxia activation. Based on the existing benchmark inhibitor O6-benzylguanine (O6-BG), four derivatives with hypoxia-reduced potential and their corresponding reduction products were synthesized. A reductase system consisting of glucose/glucose oxidase, xanthine/xanthine oxidase, and catalase were constructed, and the reduction products of the hypoxia-activated prodrugs under normoxic and hypoxic conditions were determined by high-performance liquid chromatography electrospray ionization tandem mass spectrometry (HPLC-ESI-MS/MS). The results showed that the reduction products produced under hypoxic conditions were significantly higher than that under normoxic condition. The amount of the reduction product yielded from ANBP (2-nitro-6-(3-amino) benzyloxypurine) under hypoxic conditions was the highest, followed by AMNBP (2-nitro-6-(3-aminomethyl)benzyloxypurine), 2-NBP (2-nitro-6-benzyloxypurine), and 3-NBG (O6-(3-nitro)benzylguanine). It should be noted that although the levels of the reduction products of 2-NBP and 3-NBG were lower than those of ANBP and AMNBP, their maximal hypoxic/normoxic ratios were higher than those of the other two prodrugs. Meanwhile, we also investigated the single electron reduction mechanism of the hypoxia-activated prodrugs using density functional theory (DFT) calculations. As a result, the reduction of the nitro group to the nitroso was proven to be a rate-limiting step. Moreover, the 2-nitro group of purine ring was more ready to be reduced than the 3-nitro group of benzyl. The energy barriers of the rate-limiting steps were 34–37 kcal/mol. The interactions between these prodrugs and nitroreductase were explored via molecular docking study, and ANBP was observed to have the highest affinity to nitroreductase, followed by AMNBP, 2-NBP, and 3-NBG. Interestingly, the theoretical results were generally in a good agreement with the experimental results. Finally, molecular docking and molecular dynamics simulations were performed to predict the AGT-inhibitory activity of the four prodrugs and their reduction products. In summary, simultaneous consideration of reduction potential and hypoxic selectivity is necessary to ensure that such prodrugs have good hypoxic tumor targeting. This study provides insights into the hypoxia-activated mechanism of nitro-substituted prodrugs as AGT inhibitors, which may contribute to reasonable design and development of novel tumor-targeted AGT inhibitors.