Arabidopsis COP1 suppresses root hair development by targeting type I ACS proteins for ubiquitination and degradation.

Root hairs (RHs) are an innovation of vascular plants whose development is coordinated by endogenous and environmental cues, such as ethylene and light conditions. However, the potential crosstalk between ethylene and light conditions in RH development is unclear. We report that Arabidopsis constitutive photomorphogenic 1 (COP1) integrates ethylene and light signaling to mediate RH development. Darkness suppresses RH development largely through COP1. COP1 inhibits both cell fate determination of trichoblast and tip growth of RHs based on pharmacological, genetic, and physiological analyses. Indeed, COP1 interacts with and catalyzes the ubiquitination of ACS2 and ACS6. COP1- or darkness-promoted proteasome-dependent degradation of ACS2/6 leads to a low ethylene level in underground tissues. The negative role of COP1 in RH development by downregulating ethylene signaling may be coordinated with the positive role of COP1 in hypocotyl elongation by upregulating ethylene signaling, providing an evolutionary advantage for seedling fitness.

Regulatory roles of miRNA-530 in the post-transcriptional regulation of NF-κB signaling pathway through targeted modulation of IκBα in Sebastesschlegelii.

MicroRNAs (miRNAs) are a crucial type of non-coding RNAs involved in post-transcriptional regulation. The playing essential regulatory roles in the NF-κB signaling pathway and modulate the host immune response to diverse pathogens by targeting IκBα. However, the regulatory mechanism of miRNAs in relation with IκBα in Sebastes schlegelii remains unclear. In our study, we identified two copies of IkBα gene in black rockfish (Sebastes schlegelii), namely IkBα1 and IkBα2. Moreover, we have discovered that miRNA-530 can activate the NF-κB signaling pathway by inhibiting the expression of IκBα, thereby inducing the inflammatory response. This project comprehensively investigated the interactive regulatory roles of miRNA-530 in the NF-κB signaling pathway at both cellular and in vivo levels, while also elucidating the regulatory relationships between miRNA-530 and IκBα. In conclusion, our research confirmed that miRNA-530 can target the 3'UTR region of IκBα, resulting in a decrease in the expression of IκBα at the post-transcriptional level and inhibiting its translation. The findings contribute to the understanding of the regulatory network of non-coding RNA in teleosts and its subsequent regulation of the NF-κB signaling pathway by miRNAs.

Identification of Key Ubiquitination Sites Involved in the Proteasomal Degradation of AtACS7 in Arabidopsis.

The gaseous hormone ethylene plays pivotal roles in plant growth and development. The rate-limiting enzyme of ethylene biosynthesis in seed plants is 1-aminocyclopropane-1-carboxylic acid (ACC) synthase (ACS). ACS proteins are encoded by a multigene family and the expression of ACS genes is highly regulated, especially at a post-translational level. AtACS7, the only type III ACS in Arabidopsis, is degraded in a 26S proteasome-dependent pathway. Here, by using liquid chromatography-mass spectrometry/mass spectrometry (LC-MS/MS) analysis, two lysine residues of AtACS7, lys285 (K285) and lys366 (K366), were revealed to be ubiquitin-modified in young, light-grown Arabidopsis seedlings but not in etiolated seedlings. Deubiquitylation-mimicking mutations of these residues significantly increased the stability of the AtACS7K285RK366R mutant protein in cell-free degradation assays. All results suggest that K285 and K366 are the major ubiquitination sites on AtACS7, providing deeper insights into the post-translational regulation of AtACS7 in Arabidopsis.

Comprehensive Review of Drug-Drug Interaction Prediction Based on Machine Learning: Current Status, Challenges, and Opportunities.

Detecting drug-drug interactions (DDIs) is an essential step in drug development and drug administration. Given the shortcomings of current experimental methods, the machine learning (ML) approach has become a reliable alternative, attracting extensive attention from the academic and industrial fields. With the rapid development of computational science and the growing popularity of cross-disciplinary research, a large number of DDI prediction studies based on ML methods have been published in recent years. To give an insight into the current situation and future direction of DDI prediction research, we systemically review these studies from three aspects: (1) the classic DDI databases, mainly including databases of drugs, side effects, and DDI information; (2) commonly used drug attributes, which focus on chemical, biological, and phenotypic attributes for representing drugs; (3) popular ML approaches, such as shallow learning-based, deep learning-based, recommender system-based, and knowledge graph-based methods for DDI detection. For each section, related studies are described, summarized, and compared, respectively. In the end, we conclude the research status of DDI prediction based on ML methods and point out the existing issues, future challenges, potential opportunities, and subsequent research direction.

Undesirable ER stress induced by bavachin contributed to follicular atresia in zebrafish ovary.

Fructus psoraleae (FP) is a commonly used herb with potential reproductive toxicity. Bavachin (BV), one of essential active ingredients of FP, was found to exhibit estrogenic activity, but its effect on female reproductive system remains unknown. In this study, the impact of BV on the female zebrafish reproductive system and underlying molecular mechanism were determined in vivo and ex vivo. The results showed that BV could accumulate in zebrafish ovary, leading to obvious follicular atresia and increase in gonadal index and vitellogenin content. Endoplasmic reticulum (ER) swelling and hypertrophy were observed in the BV-treated zebrafish ovary, accompanied by an increase in the expressions of ER stress and unfolded protein response (UPR) related genes, namely atf6, ire-1α and xbp1s. In the ex vivo study, BV was found to decrease the survival rate and maturation rate of oocytes, while increasing the expression of Ca2+. Additionally, BV led to an elevation in the level of estrogen receptor ESR1 and the expressions of genes involved in ER stress and UPR, including atf6, ire-1α, xbp1s, chop and perk. Moreover, molecular docking revealed that BV could directly bind to immunoglobulin heavy chain binding protein (BiP) and estrogen receptor 1 (ESR1). Besides, the alterations induced by BV could be partially reversed by fulvestrant (FULV) and 4-phenylbutyric acid (4-PBA), respectively. Thus, long-termed BV-containing medicine treatment could generate reproductive toxicity in female zebrafish by causing follicular atresia through BiP- and ESR-mediated ER stress and UPR, providing a potential target for the prevention of reproductive toxicity caused by BV.

A Scd1-mediated metabolic alteration participates in liver responses to low-dose bavachin.

Hepatotoxicity induced by bioactive constituents in traditional Chinese medicines or herbs, such as bavachin (BV) in Fructus Psoraleae, has a prolonged latency to overt drug-induced liver injury in the clinic. Several studies have described BV-induced liver damage and underlying toxicity mechanisms, but little attention has been paid to the deciphering of organisms or cellular responses to BV at no-observed-adverse-effect level, and the underlying molecular mechanisms and specific indicators are also lacking during the asymptomatic phase, making it much harder for early recognition of hepatotoxicity. Here, we treated mice with BV for 7 days and did not detect any abnormalities in biochemical tests, but found subtle steatosis in BV-treated hepatocytes. We then profiled the gene expression of hepatocytes and non-parenchymal cells at single-cell resolution and discovered three types of hepatocyte subsets in the BV-treated liver. Among these, the hepa3 subtype suffered from a vast alteration in lipid metabolism, which was characterized by enhanced expression of apolipoproteins, carboxylesterases, and stearoyl-CoA desaturase 1 (Scd1). In particular, increased Scd1 promoted monounsaturated fatty acids (MUFAs) synthesis and was considered to be related to BV-induced steatosis and polyunsaturated fatty acids (PUFAs) generation, which participates in the initiation of ferroptosis. Additionally, we demonstrated that multiple intrinsic transcription factors, including Srebf1 and Hnf4a, and extrinsic signals from niche cells may regulate the above-mentioned molecular events in BV-treated hepatocytes. Collectively, our study deciphered the features of hepatocytes in response to BV insult, decoded the underlying molecular mechanisms, and suggested that Scd1 could be a hub molecule for the prediction of hepatotoxicity at an early stage.

Generation of Herbicide-Resistant Soybean by Base Editing.

Weeds cause the largest yield loss in soybean production. The development of herbicide-resistant soybean germplasm is of great significance for weed control and yield improvement. In this study, we used the cytosine base editor (BE3) to develop novel herbicide-resistant soybean. We have successfully introduced base substitutions in GmAHAS3 and GmAHAS4 and obtained a heritable transgene-free soybean with homozygous P180S mutation in GmAHAS4. The GmAHAS4 P180S mutants have apparent resistance to chlorsulfuron, flucarbazone-sodium, and flumetsulam. In particular, the resistance to chlorsulfuron was more than 100 times that of with wild type TL-1. The agronomic performance of the GmAHAS4 P180S mutants showed no significant differences to TL-1 under natural growth conditions. In addition, we developed allele-specific PCR markers for the GmAHAS4 P180S mutants, which can easily discriminate homozygous, heterozygous mutants, and wild-type plants. This study demonstrates a feasible and effective way to generate herbicide-resistant soybean by using CRISPR/Cas9-mediated base editing.

Ferritinophagy activation and sideroflexin1-dependent mitochondrial iron overload contribute to patulin-induced cardiac inflammation and fibrosis.

Patulin (PAT) is a mycotoxin that is commonly present throughout the ecosystem where fungi grow and mainly contaminates food, soil, and water. PAT was found to be cardiotoxic in previous studies. However, the detailed mechanism has not been fully elucidated. The present study aimed to explore the role and underlying mechanism of ferroptosis in PAT-induced cardiac injury. Here, we confirmed in vivo and in vitro that ferroptosis is involved in PAT-induced myocardial inflammation and fibrosis. Mice exposed to PAT (1 and 2 mg/kg body weight/day for 14 days) exhibited myocardial inflammation and fibrosis along with disrupted iron homeostasis, elevated lipid peroxidation, depletion of glutathione peroxidase 4, and abnormal mitochondrial morphology. When primary neonatal rat cardiomyocytes (NRCMs) and H9c2 cells were exposed to PAT, ferroptosis was initiated in a dose-dependent manner, and this process could be significantly attenuated by ferrostatin-1. Mechanistically, we found that nuclear receptor coactivator (NCOA) 4, a master regulator of ferritinophagy, bound to and degraded ferritin in response to PAT treatment, thereby releasing large amounts of ferrous iron and further leading to sideroflexin (SFXN) 1-dependent mitochondrial iron overload. Conversely, knockdown of NCOA4 or SFXN1 with small interfering RNAs could effectively ameliorate ferroptotic cell death, cellular or mitochondrial iron overload and lipid peroxides accumulation. Furthermore, myocardial inflammation and fibrosis in PAT-exposed mice was alleviated by the mitochondrial iron chelator deferiprone. Overall, our findings underscore that ferritinophagy activation and SFXN1-dependent mitochondrial iron overload play critical roles in PAT-induced myocardial ferroptosis and consequent cardiotoxicity.

Certification of non-classicality in all links of a photonic star network without assuming quantum mechanics.

Networks composed of independent sources of entangled particles that connect distant users are a rapidly developing quantum technology and an increasingly promising test-bed for fundamental physics. Here we address the certification of their post-classical properties through demonstrations of full network nonlocality. Full network nonlocality goes beyond standard nonlocality in networks by falsifying any model in which at least one source is classical, even if all the other sources are limited only by the no-signaling principle. We report on the observation of full network nonlocality in a star-shaped network featuring three independent sources of photonic qubits and joint three-qubit entanglement-swapping measurements. Our results demonstrate that experimental observation of full network nonlocality beyond the bilocal scenario is possible with current technology.

A Genome-Wide Association Study to Identify Novel Candidate Genes Related to Low-Nitrogen Tolerance in Cucumber (Cucumis sativus L.).

Cucumber is one of the most important vegetables, and nitrogen is essential for the growth and fruit production of cucumbers. It is crucial to develop cultivars with nitrogen limitation tolerance or high nitrogen efficiency for green and efficient development in cucumber industry. To reveal the genetic basis of cucumber response to nitrogen starvation, a genome-wide association study (GWAS) was conducted on a collection of a genetically diverse population of cucumber (Cucumis sativus L.) comprising 88 inbred and DH accessions including the North China type, the Eurasian type, the Japanese and South China type mixed subtype, and the South China type subtype. Phenotypic evaluation of six traits under control (14 mM) and treatment (3.5 mM) N conditions depicted the presence of broad natural variation in the studied population. The GWAS results showed that there were significant differences in the population for nitrogen limitation treatment. Nine significant loci were identified corresponding to six LD blocks, three of which overlapped. Sixteen genes were selected by GO annotation associated with nitrogen. Five low-nitrogen stress tolerance genes were finally identified by gene haplotype analysis: CsaV3_3G003630 (CsNRPD1), CsaV3_3G002970 (CsNRT1.1), CsaV3_4G030260 (CsSnRK2.5), CsaV3_4G026940, and CsaV3_3G011820 (CsNPF5.2). Taken together, the experimental data and identification of candidate genes presented in this study offer valuable insights and serve as a useful reference for the genetic enhancement of nitrogen limitation tolerance in cucumbers.

Spatial transcriptomes and microbiota reveal immune mechanism that respond to pathogen infection in the posterior intestine of Sebastes schlegelii.

The intestine is a site of immune cell priming at birth. Therefore, spatial transcriptomes were performed to define how the transcriptomic landscape was spatially organized in the posterior intestine of Sebastes schlegelii following Edwardsiella piscicida infection. In the healthy condition, we identified a previously unappreciated molecular regionalization of the posterior intestine. Following bacterial infection, most immune-related genes were identified in mucosa layer. Moreover, investigation of immune-related genes and genes in immune-related KEGG pathways based on spatial transcriptomes shed light on which sections of these genes are in the posterior intestine. Meanwhile, the high expression of genes related to regeneration also indicated that the posterior intestine was responding to the invasion of pathogens by constantly proliferating new cells. In addition, the increasing microbiota communities indicated that these bacteria maintained posterior intestine integrity and shaped the mucosal immune system. Taken together, spatial transcriptomes and microbiota compositions have significant implications for understanding the immune mechanism that responds to E. piscicida infection in the posterior intestine of S. schlegelii, which also provides a theoretical basis for the spatial distribution of immune genes and changes in bacterial flora in other teleosts in the process of resisting pathogens.

Domain Swapping between AtACS7 and PpACL1 Results in Chimeric ACS-like Proteins with ACS or Cß-S Lyase Single Enzymatic Activity.

The gaseous hormone ethylene plays a pivotal role in plant growth and development. In seed plants, the key rate-limiting enzyme that controls ethylene biosynthesis is ACC synthase (ACS). ACS has, for a long time, been believed to be a single-activity enzyme until we recently discovered that it also possesses Cß-S lyase (CSL) activity. This discovery raises fundamental questions regarding the biological significance of the dual enzymatic activities of ACS. To address these issues, it is highly necessary to obtain ACS mutants with either ACS or CSL single activity. Here, domain swapping between Arabidopsis AtACS7 and moss CSL PpACL1 were performed. Enzymatic activity assays of the constructed chimeras revealed that, R10, which was produced by replacing AtACS7 box 6 with that of PpACL1, lost ACS but retained CSL activity, whereas R12 generated by box 4 substitution lost CSL and only had ACS activity. The activities of both chimeric proteins were compared with previously obtained single-activity mutants including R6, AtACS7Q98A, and AtACS7D245N. All the results provided new insights into the key residues required for ACS and CSL activities of AtACS7 and laid an important foundation for further in-depth study of the biological functions of its dual enzymatic activities.

Pinelliae rhizoma alleviated acute lung injury induced by lipopolysaccharide via suppressing endoplasmic reticulum stress-mediated NLRP3 inflammasome.

Pinelliae rhizoma (PR), one kind of commonly-used Chinese herbs, is generally prescribed to treat various respiratory diseases, including acute lung injury (ALI). However, the accurate bioactive ingredients of PR and the underlying pharmacological mechanism have both not been fully elucidated. Therefore, this study aimed to identify the bioactive ingredients that could alleviate lipopolysaccharide (LPS)-induced ALI and explore the possible mechanism involved. Our results confirmed that LPS infection indeed caused acute inflammatory damage in mice lung, accompanying with the enhancement of IL-1ß contents and the activation of the NLRP3 inflammasome in lung tissue and macrophagocyte, all of which were remarkably ameliorated by PR treatment. Next, mechanistically, LPS was found to trigger endoplasmic reticulum (ER) stress and downstream cellular calcium ions (Ca2+) release via activating Bip/ATF4/CHOP signaling pathway. Like PR, 4-PBA (a specific inhibitor of ER stress) not only obviously reversed Bip/ATF4/CHOP-mediated ER stress, but also significantly attenuated LPS-induced activation of the NLRP3 inflammasome. Furthermore, the bioactive ingredients of PR, which generated the anti-inflammatory effects, were screened by metabolomics and network pharmacology. In vitro experiments showed that chrysin, dihydrocapsaicin, and 7,8-dihydroxyflavone (7,8-DHF) notably suppressed LPS-induced ER stress and following NLRP3 inflammasome activation. In conclusion, our findings suggested that PR alleviated LPS-induced ALI by inhibiting ER stress-mediated NLRP3 inflammasome activation, which is mainly relevant with these three bioactive ingredients. This study provided a theoretical basis for the clinical application of PR to treat ALI, and these bioactive ingredients of PR would be promising therapeutic drugs for the treatment of ALI.

Identification and Functional Analysis of Key Autophosphorylation Residues of Arabidopsis Senescence Associated Receptor-like Kinase.

Reversible protein phosphorylation mediated by protein kinases and phosphatases plays important roles in the regulation of leaf senescence. We previously reported that the senescence-associated leucine-rich repeat receptor-like kinase AtSARK autophosphorylates on both serine/threonine and tyrosine residues and functions as a positive regulator of Arabidopsis leaf senescence; the senescence-suppressed protein phosphatase SSPP interacts with and dephosphorylates the cytoplasmic domain of AtSARK, thereby negatively regulating leaf senescence. Here, 27 autophosphorylation residues of AtSARK were revealed by mass spectrometry analysis, and six of them, including two Ser, two Thr, and two Tyr residues, were further found to be important for the biological functions of AtSARK. All site-directed mutations of these six residues that resulted in decreased autophosphorylation level of AtSARK could significantly inhibit AtSARK-induced leaf senescence. In addition, mutations mimicking the dephosphorylation form of Ser384 (S384A) or the phosphorylation form of Tyr413 (Y413E) substantially reduced the interaction between AtSARK and SSPP. All results suggest that autophosphorylation of AtSARK is essential for its functions in promoting leaf senescence. The possible roles of S384 and Y413 residues in fine-tuning the interaction between AtSARK and SSPP are discussed herein.

N7 -SSPP fusion gene improves salt stress tolerance in transgenic Arabidopsis and soybean through ROS scavenging.

Considerable signal crosstalk exists in the regulatory network of senescence and stress response. Numerous senescence-associated genes are also involved in plant stress tolerance. However, the underlying mechanisms and application potential of these genes in stress-tolerant crop breeding remain poorly explored. We found that overexpression of SENESCENCE-SUPPRESSED PROTEIN PHOSPHATASE (SSPP), a negative regulator of leaf senescence, significantly improved plant salt tolerance by increasing reactive oxygen species (ROS) scavenging in both Arabidopsis and soybean. However, overexpression of SSPP severely suppressed normal plant growth, limiting its direct use in agriculture. We previously revealed that the N-terminal 1-14 residues of ACS7 (termed 'N7 ') negatively regulated its protein stability through the ubiquitin/proteasome pathway, and the N7 -mediated protein degradation was suppressed by environmental and senescence signals. To avoid the adverse effects of SSPP, the N7 element was fused to the N-terminus of SSPP. We demonstrated that N7 -SSPP fusion gene effectively rescued SSPP-induced growth suppression but maintained enhanced salt tolerance in Arabidopsis and soybean. Particularly, N7 -SSPP enhanced tolerance to long-term salt stress and increased seed yield in soybean. These results suggest that N7 -SSPP overcomes the disadvantages of SSPP on plant growth inhibition and effectively improves salt tolerance through enhanced ROS scavenging, providing an effective strategy of using posttranslational regulatory element for salt-tolerant crop breeding.

Ginsenoside Rb1 Ameliorated Bavachin-Induced Renal Fibrosis via Suppressing Bip/eIF2α/CHOP Signaling-Mediated EMT.

The nephrotoxicity of Fructus Psoraleae, an effective traditional Chinese medicine for vitiligo treatment, has been reported. As one of the main toxic components in Fructus Psoraleae, bavachin (BV) was considered to be related to Fructus Psoraleae-caused adverse outcomes, but the direct evidence and molecular mechanism underlying BV-induced nephrotoxicity are not well elucidated. Therefore, this study was designed to confirm whether BV would cause toxic effects on the kidney and explore the possible mode of action. Our results demonstrated that days' treatment with 0.5 µM BV indeed caused obvious renal fibrosis in the zebrafish kidney. The obvious E- to N-cadherin switch and the expressions of proteins promoting epithelial-mesenchymal transition (EMT) were observed in BV-treated human renal tubular epithelial and zebrafish kidneys. In addition, elevated reactive oxygen species (ROS) levels and Bip/eIF2α/CHOP-mediated endoplasmic reticulum (ER) stress and the unfolded protein response (UPR) were caused by BV, both of which could be reversed by ROS scavenger N-acetyl-L-cysteine (NAC). Also, blocking ER stress-caused cytoplasmic Ca2+ overload with 4-PBA notably alleviated BV-induced alterations in key molecular events related to EMT and renal fibrosis. Furthermore, of the natural compounds subjected to screening, ginsenoside Rb1 significantly downregulated BV-induced ER stress by inhibiting ROS generation and following the activation of Bip/eIF2α/CHOP signaling in HK2 cells. Subsequently, BV-triggered EMT and renal fibrosis were both ameliorated by ginsenoside Rb1. In summary, our findings suggested that BV-induced ROS promoted the appearance of EMT and renal fibrosis mainly via Bip/eIF2α/CHOP-mediated ER stress. This ER stress-related toxic pathway might be a potential intervention target for BV-caused renal fibrosis, and ginsenoside Rb1 would be a promising drug against BV- or Fructus Psoraleae-induced nephrotoxicity.

Abnormal BMI in Male and/or Female Partners Are Deleterious for Embryonic Development and Pregnancy Outcome During ART Process: A Retrospective Study.

Objective: To investigate the effect of BMI in male and/or female partners on embryo development and clinical pregnancy outcome during ART. Methods: Data of 11,130 cycles between January 2018 and December 2020 were retrospectively analyzed. They were divided into Group A, B, C, and D based on couples' BMI values, also were divided into Group I, II, III and IV in IVF cycles and Group i, ii, iii, and iv in ICSI cycles. After grouping, inter-group indicators comparisons and logistic regression analysis were performed. Results: In IVF cycles, CPR in Group I and Group III were higher than Group IV. In Group III, it was higher than Group II. The AR in Group IV was higher, but the LBR was lower than Group I, Group II, and Group III. Logistic regression analysis results suggested that AR in Group IV was higher than that in Group I in IVF cycles, whereas LBR was lower.In ICSI cycles, high-quality embryo rate in Group i and Group ii were both higher than that in Group iii and Group iv. The CPR in Group i was higher than Group ii and Group iv, and in Group iii was higher than Group ii and Group iv. The AR in Group i was lower than Group iii and Group iv, and AR in Group ii was lower than Group iv. LBR, in Group I it was higher than Group ii, Group iii, and Group iv. Logistic regression analysis results suggested CPR in Group ii was significantly lower than that in Group i. AR in Group iii was considerably higher than that in Group i. LBR in Group ii and Group iv were significantly lower than that in Group i. Conclusion: Female higher BMI was not conducive to the formation of high-quality embryos in ICSI cycle. Female and/or male BMI affected AR and LBR more than CPR not only in IVF cycles, but also in ICSI cycles.

Machine learning to predict metabolic drug interactions related to cytochrome P450 isozymes.

Drug-drug interaction (DDI) often causes serious adverse reactions and thus results in inestimable economic and social loss. Currently, comprehensive DDI evaluation has become a major challenge in pharmaceutical research due to the time-consuming and costly process of the experimental assessment and it is of high necessity to develop effective in silico methods to predict and evaluate DDIs accurately and efficiently. In this study, based on a large number of substrates and inhibitors related to five important CYP450 isozymes (CYP1A2, CYP2C9, CYP2C19, CYP2D6 and CYP3A4), a series of high-performance predictive models for metabolic DDIs were constructed by two machine learning methods (random forest and XGBoost) and 4 different types of descriptors (MOE_2D, CATS, ECFP4 and MACCS). To reduce the uncertainty of individual models, the consensus method was applied to yield more reliable predictions. A series of evaluations illustrated that the consensus models were more reliable and robust for the DDI predictions of new drug combination. For the internal validation, the whole prediction accuracy and AUC value of the DDI models were around 0.8 and 0.9, respectively. When it was applied to the external datasets, the model accuracy was 0.793 and 0.795 for multi-level validation and external validation, respectively. Furthermore, we also compared our model with some recently published tools and then applied the final model to predict FDA-approved drugs and proposed 54,013 possible drug pairs with potential DDIs. In summary, we developed a powerful DDI predictive model from the perspective of the CYP450 enzyme family and it will help a lot in the future drug development and clinical pharmacy research.