Mechanism of AGT-Mediated Repair of dG-dC Cross-Links in the Drug Resistance to Chloroethylnitrosoureas: Molecular Docking, MD Simulation, and ONIOM (QM/MM) Investigation.

Chloroethylnitrosoureas (CENUs) are important chemotherapies applied in the treatment of cancer. They exert anticancer activity by inducing DNA interstrand cross-links (ICLs) via the formation of two O6-alkylguanine intermediates, O6-chloroethylguanine (O6-ClEtG) and N1,O6-ethanoguanine (N1,O6-EtG). However, O6-alkylguanine-DNA alkyltransferase (AGT), a DNA-repair enzyme, can restore the O6-alkylguanine damages and thereby obstruct the formation of ICLs (dG-dC cross-link). In this study, the inhibitory mechanism of ICL formation was investigated to elucidate the drug resistance of CENUs mediated by AGT in detail. Based on the structures of the substrate-enzyme complexes obtained from docking and MD simulations, two ONIOM (QM/MM) models with different sizes of the QM region were constructed. The model with a larger QM region, which included the substrate (O6-ClEtG or N1,O6-EtG), a water molecule, and five residues (Tyr114, Cys145, His146, Lys165, and Glu172) in the active pocket of AGT, accurately described the repairing reaction and generated the results coinciding with the experimental outcomes. The repair process consists of two sequential steps: hydrogen transfer to form a thiolate anion on Cys145 and alkyl transfer from the O6 site of guanine (the rate-limiting step). The repair of N1,O6-EtG was more favorable than that of O6-ClEtG from both kinetics and thermodynamics aspects. Moreover, the comparison of the repairing process with the formation of dG-dC cross-link and the inhibition of AGT by O6-benzylguanine (O6-BG) showed that the presence of AGT could effectively interrupt the formation of ICLs leading to drug resistance, and the inhibition of AGT by O6-BG that was energetically more favorable than the repair of O6-ClEtG could not prevent the repair of N1,O6-EtG. Therefore, it is necessary to completely eliminate AGT activity before CENUs medication to enhance the chemotherapeutic effectiveness. This work provides reasonable explanations for the supposed mechanism of AGT-mediated drug resistance of CENUs and will assist in the development of novel CENU chemotherapies and their medication strategies.

The rat acute oral toxicity of trifluoromethyl compounds (TFMs): a computational toxicology study combining the 2D-QSTR, read-across and consensus modeling methods.

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.

A multi-functional hypoxia/esterase dual stimulus responsive and hyaluronic acid-based nanomicelle for targeting delivery of chloroethylnitrosouea.

Carmustine (BCNU), a vital type of chloroethylnitrosourea (CENU), inhibits tumor cells growth by inducing DNA damage at O6 position of guanine and eventually forming dG-dC interstrand cross-links (ICLs). However, the clinical application of BCNU is hindered to some extent by the absence of tumor selectivity, poor stability and O6-alkylguanine-DNA alkyltransferase (AGT) mediated drug resistance. In recent years, tumor microenvironment has been widely utilized for advanced drug delivery. In the light of the features of tumor microenvironment, we constructed a multifunctional hypoxia/esterase-degradable nanomicelle with AGT inhibitory activity named HACB NPs for tumor-targeting BCNU delivery and tumor sensitization. HACB NPs was self-assembled from hyaluronic acid azobenzene AGT inhibitor conjugates, in which O6-BG analog acted as an AGT inhibitor, azobenzene acted as a hypoxia-responsive linker and carboxylate ester bond acted as both an esterase-sensitive switch and a connector with hyaluronic acid (HA). The obtained HACB NPs possessed good stability, favorable biosafety and hypoxia/esterase-responsive drug-releasing ability. BCNU-loaded HACB/BCNU NPs exhibited superior cytotoxicity and apoptosis-inducing ability toward the human uterine cervix carcinoma HeLa cells compared with traditional combined medication of BCNU plus O6-BG. In vivo studies further demonstrated that after a selective accumulation in the tumor site, the micelles could respond to hypoxic tumor tissue for rapid drug release to an effective therapeutic dosage. Thus, this multifunctional stimulus-responsive nanocarrier could be a new promising strategy to enhance the anticancer efficacy and reduce the side effects of BCNU and other CENUs.

Effects of terbuthylazine on myocardial oxidative stress and ferroptosis via Nrf2/HO-1 signaling pathway in broilers.

Terbuthylazine (TBA) is one of the most commonly used and effective herbicides. However, due to its affinity for soil organic matter and water solubility, TBA can lead to biological health concerns. This study exposed broilers to TBA (0 mg/kg bw, 0.4 mg/kg bw, 4 mg/kg bw) for 28 days. The results showed significant pathological damage in broiler myocardial tissue, such as widening of the interstitial space, rupture of muscle fibers, and deposition of myocardial collagen fibers. In addition, Under the 0.4 mg/kg bw TBA exposure, myocardial oxidative stress was observed in broilers, which was accompanied by the activation of Nrf2/HO-1 pathway and the increased protein and mRNA levels of NQO1, NOX2 and SOD2 antioxidant enzymes. However, Nrf2/HO-1 protein and mRNA levels were reversed at 4 mg/kg bw TBA exposure. Meanwhile, the Nrf2/HO-1 mediated antioxidant defense was impaired. In contrast with the low dose, the protein and gene expression levels of NQO1, NOX2, and SOD2 were reduced in 4 mg/kg bw TBA group. The expression of GPX4 and SLC7A11 was significantly downregulated at both protein and mRNA levels. Beyond that, ACSL4 expression was significantly up-regulated, and the protein result was consistent with the mRNA expression, demonstrating the occurrence of ferroptosis. In general, TBA exposure activated the Nrf2/HO-1 pathway, resulting in ferroptosis. This study links ferroptosis to the Nrf2/HO-1 pathway, providing new insights into the potential role of TBA in myocardial toxicity.

Topological cavity laser with valley edge states.

Topological edge states (ES) arise at the boundary between spatial domains with diverse topological properties in photonic crystals, which can transmit unidirectionally to suppress the backscattering and robustly to be immune to defects and disorders. In addition, optical devices with arbitrary geometries of cavities, such as lasers, are expected to be designed on the basis of ES. Herein, we first propose a topological cavity laser based on a honeycomb lattice of ring holes with the bearded interface in two-dimensional (2D) all-dielectric valley photonic crystals (VPhCs) at telecommunication wavelengths. Specifically, we construct a topological cavity using topological valley edge states (VES) and further study the lasing action of the optically pumped cavity with high-quality factors. Our findings could provide opportunities for practical applications of VES-based lasers as ultra-small light sources with the topological protection.

BMP/Smad Pathway Is Involved in Lithium Carbonate-Induced Neural-Tube Defects in Mice and Neural Stem Cells.

Neural-tube defects (NTDs) are one type of the most serious birth defects. Studies have shown that inositol deficiency is closely related to the occurrence of NTDs. Bone morphogenetic protein (BMP)-mediated Smad signaling pathways have been implicated in neurogenesis and neural-tube closure. However, the role of the BMP/Smad pathway in inositol-deficiency-induced NTDs remains unclear. Inositol-deficiency models in C57 mice and mouse neural stem cells (mNSCs) were induced with Li2CO3 treatment or inositol withdrawal. The role of the BMP/Smad pathway in the regulation of cell proliferation and the development of NTDs was determined utilizing qRT-PCR, HE staining, Western blot, immunostaining, MTT assay, EdU staining, and flow cytometry. The intraperitoneal injection of Li2CO3 at Embryonic Day 7.5 induced the occurrence of NTDs. The mRNA levels of Bmp2, Bmp4, Smad1, Smad5, Smad8 and Runx2, the phosphorylation of Smad1/5/8, and the nuclear translocation of Runx2 were significantly increased in NTD embryonic brain tissues and mNSCs exposed to Li2CO3 or an inositol-free medium, which were suppressed by BMP receptor selective inhibitor LDN-193189. The Li2CO3-induced phosphorylation of Smad1/5/8 was inhibited by inositol supplementation. Cell proliferation was significantly promoted by Li2CO3 exposure or the absence of inositol in mNSCs, which was reversed by LDN-193189. These results suggest that the activation of the BMP/Smad signaling pathway might play an important role in the development of NTDs induced by maternal Li2CO3 exposure via inositol deficiency.

Identification of Pharmacophoric Fragments of DYRK1A Inhibitors Using Machine Learning Classification Models.

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.

Chemometric QSAR modeling of acute oral toxicity of Polycyclic Aromatic Hydrocarbons (PAHs) to rat using simple 2D descriptors and interspecies toxicity modeling with mouse.

The information of the acute oral toxicity for most polycyclic aromatic hydrocarbons (PAHs) in mammals are lacking due to limited experimental resources, leading to a need to develop reliable in silico methods to evaluate the toxicity endpoint. In this study, we developed the quantitative structure-activity relationship (QSAR) models by genetic algorithm (GA) and multiple linear regression (MLR) for the rat acute oral toxicity (LD50) of PAHs following the strict validation principles of QSAR modeling recommended by OECD. The best QSAR model comprised eight simple 2D descriptors with definite physicochemical meaning, which showed that maximum atom-type electrotopological state, van der Waals surface area, mean atomic van der Waals volume, and total number of bonds are main influencing factors for the toxicity endpoint. A true external set (554 compounds) without rat acute oral toxicity values, and 22 limit test compounds, were firstly predicted along with reliability assessment. We also compared our proposed model with the OPERA predictions and recently published literature to prove the prediction reliability. Furthermore, the interspecies toxicity (iST) models of PAHs between rat and mouse were also established, validated and employed for filling data gap. Overall, our developed models should be applicable to new or untested or not yet synthesized PAHs falling within the applicability domain (AD) of the models for rapid acute oral toxicity prediction, thus being important for environmental or personal exposure risk assessment under regulatory frameworks.

Quantitative Structure-Activity Relationship (QSAR) Studies on the Toxic Effects of Nitroaromatic Compounds (NACs): A Systematic Review.

Nitroaromatic compounds (NACs) are ubiquitous in the environment due to their extensive industrial applications. The recalcitrance of NACs causes their arduous degradation, subsequently bringing about potential threats to human health and environmental safety. The problem of how to effectively predict the toxicity of NACs has drawn public concern over time. Quantitative structure-activity relationship (QSAR) is introduced as a cost-effective tool to quantitatively predict the toxicity of toxicants. Both OECD (Organization for Economic Co-operation and Development) and REACH (Registration, Evaluation and Authorization of Chemicals) legislation have promoted the use of QSAR as it can significantly reduce living animal testing. Although numerous QSAR studies have been conducted to evaluate the toxicity of NACs, systematic reviews related to the QSAR modeling of NACs toxicity are less reported. The purpose of this review is to provide a thorough summary of recent QSAR studies on the toxic effects of NACs according to the corresponding classes of toxic response endpoints.

Synergistic Effect between Human Papillomavirus 18 and 4-(Methylnitrosamino)-1-(3-pyridyl)-1-butanone on Malignant Transformation of Immortalized SHEE Cells.

4-(Methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK) is an important tobacco-specific nitrosamine (TSNA) that induces malignant tumors in rodents. High-risk human papillomavirus (hr-HPV) infection is an important cause of several human cancers. Epidemiological evidence has shown that HPV cooperatively induces carcinogenesis with tobacco smoke. In the present study, the synergistic carcinogenesis of NNK and HPV18 was investigated. Immortalized human esophageal epithelial SHEE cells containing the HPV18 E6E7 gene were constructed by lentiviral transfection. SHEE-E6E7 cells were exposed to NNK along with SHEE-V cells without HPV18 E6E7 as a negative control. The cooperation of NNK and HPV was examined by wound-healing, transwell, and colony-forming assays. The results showed that NNK exposure promoted the migration, invasion, and proliferation abilities of both SHEE-E6E7 and SHEE-V cells; however, the changes in these phenotypic features were remarkably stronger in SHEE-E6E7 cells than those in SHEE-V cells. Our findings indicate that NNK promotes malignant transformation of human esophageal epithelial cells and suggest a synergistic carcinogenesis with the HPV18 E6E7 oncogene. As reported previously, the formation of pyridyloxybutylated DNA adducts is a crucial step in NNK-mediated carcinogenesis. In order to clarify the influence of HPV on the formation of NNK-induced DNA adducts, the amounts of 4-hydroxy-1-(3-pyridyl)-1-butanone (HPB)-releasing DNA adducts were determined using high-performance liquid chromatography-electrospray ionization-tandem mass spectrometry. We observed that the levels of HPB-releasing adducts in SHEE-E6E7 cells were significantly higher (p < 0.01) than those of SHEE-V cells, which was in line with results of the phenotypic assays. In conclusion, this study provides direct evidence that NNK and HPV18 exhibit a synergistic effect on formation of DNA adducts, resulting in malignant transformation of esophageal epithelial cells. Such knowledge on the interaction between infection and smoking habits in the development of cancers informs cancer-prevention strategies. Further studies to delineate the molecular mechanism and to identify specific intervention targets are worthwhile.

Identification and Biological Evaluation of CK2 Allosteric Fragments through Structure-Based Virtual Screening.

Casein kinase II (CK2) is considered as an attractive cancer therapeutic target, and recent efforts have been made to develop its ATP-competitive inhibitors. However, achieving selectivity with respect to related kinases remains challenging due to the highly conserved ATP-binding pocket of kinases. Allosteric inhibitors, by targeting the much more diversified allosteric site relative to the highly conserved ATP-binding pocket, might be a promising strategy with the enhanced selectivity and reduced toxicity than ATP-competitive inhibitors. The previous studies have highlighted the traditional serendipitousity of discovering allosteric inhibitors owing to the complicate allosteric modulation. In this current study, we identified the novel allosteric inhibitors of CK2α by combing structure-based virtual screening and biological evaluation methods. The structure-based pharmacophore model was built based on the crystal structure of CK2α-compound 15 complex. The ChemBridge fragment library was searched by evaluating the fit values of these molecules with the optimized pharmacophore model, as well as the binding affinity of the CK2α-ligand complexes predicted by Alloscore web server. Six hits forming the holistic interaction mechanism with the αD pocket were retained after pharmacophore- and Alloscore-based screening for biological test. Compound 3 was found to be the most potent non-ATP competitive CK2α inhibitor (IC50 = 13.0 µM) with the anti-proliferative activity on A549 cancer cells (IC50 = 23.1 µM). Our results provide new clues for further development of CK2 allosteric inhibitors as anti-cancer hits.

Structure-based identification of novel CK2 inhibitors with a linear 2-propenone scaffold as anti-cancer agents.

Protein kinase CK2 has emerged as an attractive cancer therapeutic target. Previous studies have highlighted the challenge of optimizing CK2 ATP-competitive inhibitors that have low druggability due to their polycyclic ring scaffolds. Therefore the development of novel inhibitors with non-polycyclic scaffolds emerges as a promising strategy for drug discovery targeting CK2. In this current study, based on the similar predicted binding poses of the linear 2-propenone scaffold of isoliquiritigenin with that of the polycyclic inhibitor CX-4945, a series of 2-propenone derivatives containing an amine-substituted five-membered heterocycle and a benzoic acid were designed, synthesized and evaluated for their in vitro CK2 inhibition and anti-cancer activity. Compound 8b was found to be the most potent CK2 inhibitor (IC50 = 0.6 µM) with the anti-proliferative activity on HepG2 cancer cells (IC50 = 14 µM), compared to the activity of isoliquiritigenin (IC50 = 17 µM and 51 µM, respectively). Molecular docking was performed to understand the binding modes of the newly designed 2-propenone derivatives with CK2. Compound 8b formed the most favorable network of hydrogen bonds with both the hinge region and positive area. Our results indicate that CK2 derivatives with a linear 2-propenone scaffold are promising candidates for anti-cancer drug discovery.

Mass Spectrometric Quantitation of Pyridyloxobutyl DNA Phosphate Adducts in Rats Chronically Treated with N'-Nitrosonornicotine.

The tobacco-specific carcinogens N'-nitrosonornicotine (NNN) and 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK) require metabolic activation to exert their carcinogenicity. NNN and NNK are metabolized to the same reactive diazonium ions, which alkylate DNA forming pyridyloxobutyl (POB) DNA base and phosphate adducts. We have characterized the formation of both POB DNA base and phosphate adducts in NNK-treated rats and the formation of POB DNA base adducts in NNN-treated rats. However, POB DNA phosphate adducts in NNN-treated rats are still uncharacterized. In this study, we quantified the levels of POB DNA phosphate adducts in tissues of rats chronically treated with ( S)-NNN or ( R)-NNN for 10, 30, 50, and 70 weeks during a carcinogenicity study. The highest amounts of POB DNA phosphate adducts were observed in the esophagus of the ( S)-NNN-treated rats, with a maximum level of 5400 ± 317 fmol/mg DNA at 50 weeks. The abundance of POB DNA phosphate adducts in the esophagus was consistent with the results of the carcinogenicity study showing that the esophagus was the primary site of tumor formation from treatment with ( S)-NNN. Compared to the ( R)-NNN group, the levels of POB DNA phosphate adducts were higher in the oral mucosa, esophagus, and liver, while lower in the nasal mucosa of the ( S)-NNN-treated rats. Among 10 combinations of all isomers of POB DNA phosphate adducts, Ap(POB)C and combinations with thymidine predominated across all the rat tissues examined. In the primary target tissue, esophageal mucosa, Ap(POB)C accounted for ∼20% of total phosphate adducts in the ( S)-NNN treatment group throughout the 70 weeks, with levels ranging from 780 ± 194 to 1010 ± 700 fmol/mg DNA. The results of this study showed that POB DNA phosphate adducts were present in high levels and persisted in target tissues of rats chronically treated with ( S)- or ( R)-NNN. These results improve our understanding of DNA damage during NNN-induced carcinogenesis. The predominant POB DNA phosphate isomers observed, such as Ap(POB)C, may serve as biomarkers for monitoring chronic exposure of tobacco-specific nitrosamines in humans.

Prediction on the mutagenicity of nitroaromatic compounds using quantum chemistry descriptors based QSAR and machine learning derived classification methods.

Nitroaromatic compounds (NACs) are an important type of environmental organic pollutants. However, it is lack of sufficient information relating to their potential adverse effects on human health and the environment due to the limited resources. Thus, using in silico technologies to assess their potential hazardous effects is urgent and promising. In this study, quantitative structure activity relationship (QSAR) and classification models were constructed using a set of NACs based on their mutagenicity against Salmonella typhimurium TA100 strain. For QSAR studies, DRAGON descriptors together with quantum chemistry descriptors were calculated for characterizing the detailed molecular information. Based on genetic algorithm (GA) and multiple linear regression (MLR) analyses, we screened descriptors and developed QSAR models. For classification studies, seven machine learning methods along with six molecular fingerprints were applied to develop qualitative classification models. The goodness of fitting, reliability, robustness and predictive performance of all developed models were measured by rigorous statistical validation criteria, then the best QSAR and classification models were chosen. Moreover, the QSAR models with quantum chemistry descriptors were compared to that without quantum chemistry descriptors and previously reported models. Notably, we also obtained some specific molecular properties or privileged substructures responsible for the high mutagenicity of NACs. Overall, the developed QSAR and classification models can be utilized as potential tools for rapidly predicting the mutagenicity of new or untested NACs for environmental hazard assessment and regulatory purposes, and may provide insights into the in vivo toxicity mechanisms of NACs and related compounds.

Reductive Activity and Mechanism of Hypoxia- Targeted AGT Inhibitors: An Experimental and Theoretical Investigation.

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.

QSAR and Classification Study on Prediction of Acute Oral Toxicity of N-Nitroso Compounds.

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 Q²loo = 0.7533, R² = 0.8071, Q²ext = 0.7041 and R²ext = 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.

Insights into the Impact of Linker Flexibility and Fragment Ionization on the Design of CK2 Allosteric Inhibitors: Comparative Molecular Dynamics Simulation Studies.

Protein kinase is a novel therapeutic target for human diseases. The off-target and side effects of ATP-competitive inhibitors preclude them from the clinically relevant drugs. The compounds targeting the druggable allosteric sites outside the highly conversed ATP binding pocket have been identified as promising alternatives to overcome current barriers of ATP-competitive inhibitors. By simultaneously interacting with the αD region (new allosteric site) and sub-ATP binding pocket, the attractive compound CAM4066 was named as allosteric inhibitor of CK2α. It has been demonstrated that the rigid linker and non-ionizable substituted fragment resulted in significant decreased inhibitory activities of compounds. The molecular dynamics simulations and energy analysis revealed that the appropriate coupling between the linker and pharmacophore fragments were essential for binding of CAM4066 with CK2α. The lower flexible linker of compound 21 lost the capability of coupling fragments A and B to αD region and positive area, respectively, whereas the methyl benzoate of fragment B induced the re-orientated Pre-CAM4066 with the inappropriate polar interactions. Most importantly, the match between the optimized linker and pharmacophore fragments is the challenging work of fragment-linking based drug design. These results provide rational clues to further structural modification and development of highly potent allosteric inhibitors of CK2.

In Silico Prediction of O6-Methylguanine-DNA Methyltransferase Inhibitory Potency of Base Analogs with QSAR and Machine Learning Methods.

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 Q²Loo = 0.83, R² = 0.87, Q²ext = 0.67, and R²ext = 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.

Exploring the Pivotal Role of the CK2 Hinge Region Sub-Pocket in Binding with Tricyclic Quinolone Analogues by Computational Analysis.

Protein kinase CK2 has been considered as an attractive therapeutic target of cancer therapy. The tricyclic quinoline compound CX-4945 is the first representative of CK2 inhibitors used in human clinical trials. The binding of non-2,6-naphtyridine substituted compounds 27e (IC50 > 500 nM) and 27h (IC50 > 1000 nM) to CK2 is abolished. However, the unbinding mechanisms due to the key pharmacophore group replacement of compounds 27e and 27h are unveiled. In the present work, combined computational analysis was performed to investigate the underlying structural basis of the low-affinity of two systems. As indicated in the results, the loss of hydrogen bonds between the non-2,6-naphtyridine and the hinge region destroyed the proper recognition of the two complexes. Besides, the allosteric mechanisms between the deviated ligands and the changed regions (G-loop, C-loop and ß4/ß5 loop) are proposed. Furthermore, energetic analysis was evaluated by detailed energy calculation and residue-based energy decomposition. More importantly, the summary of known polar pharmacophore groups elucidates the pivotal roles of hinge region sub-pocket in the binding of CK2 inhibitors. These results provide rational clues to the fragment-based design of more potent CK2 inhibitors.

Synthesis and antitumor activity evaluation of a novel combi-nitrosourea prodrug: Designed to release a DNA cross-linking agent and an inhibitor of O(6)-alkylguanine-DNA alkyltransferase.

The drug resistance of CENUs induced by O(6)-alkylguanine-DNA alkyltransferase (AGT), which repairs the O(6)-alkylated guanine and subsequently inhibits the formation of dG-dC cross-links, hinders the application of CENU chemotherapies. Therefore, the discovery of CENU analogs with AGT inhibiting activity is a promising approach leading to novel CENU chemotherapies with high therapeutic index. In this study, a new combi-nitrosourea prodrug 3-(3-(((2-amino-9H-purin-6-yl)oxy)methyl)benzyl)-1-(2-chloroethyl)-1-nitrosourea (6), designed to release a DNA cross-linking agent and an inhibitor of AGT, was synthesized and evaluated for its antitumor activity and ability to induce DNA interstrand cross-links (ICLs). The results indicated that 6 exhibited higher cytotoxicity against mer(+) glioma cells compared with ACNU, BCNU, and their respective combinations with O(6)-benzylguanine (O(6)-BG). Quantifications of dG-dC cross-links induced by 6 were performed using HPLC-ESI-MS/MS. Higher levels of dG-dC cross-link were observed in 6-treated human glioma SF763 cells (mer(+)), whereas lower levels of dG-dC cross-link were observed in 6-treated calf thymus DNA, when compared with the groups treated with BCNU and ACNU. The results suggested that the superiority of 6 might result from the AGT inhibitory moiety, which specifically functions in cells with AGT activity. Molecular docking studies indicated that five hydrogen bonds were formed between the O(6)-BG analogs released from 6 and the five residues in the active pocket of AGT, which provided a reasonable explanation for the higher AGT-inhibitory activity of 6 than O(6)-BG.