Advancing Adverse Drug Reaction Prediction with Deep Chemical Language Model for Drug Safety Evaluation.

The accurate prediction of adverse drug reactions (ADRs) is essential for comprehensive drug safety evaluation. Pre-trained deep chemical language models have emerged as powerful tools capable of automatically learning molecular structural features from large-scale datasets, showing promising capabilities for the downstream prediction of molecular properties. However, the performance of pre-trained chemical language models in predicting ADRs, especially idiosyncratic ADRs induced by marketed drugs, remains largely unexplored. In this study, we propose MoLFormer-XL, a pre-trained model for encoding molecular features from canonical SMILES, in conjunction with a CNN-based model to predict drug-induced QT interval prolongation (DIQT), drug-induced teratogenicity (DIT), and drug-induced rhabdomyolysis (DIR). Our results demonstrate that the proposed model outperforms conventional models applied in previous studies for predicting DIQT, DIT, and DIR. Notably, an analysis of the learned linear attention maps highlights amines, alcohol, ethers, and aromatic halogen compounds as strongly associated with the three types of ADRs. These findings hold promise for enhancing drug discovery pipelines and reducing the drug attrition rate due to safety concerns.

A pilot study on identifying gene signatures as markers for predicting patient response to antiseizure medications.

The majority of the biomarkers were associated with the diagnosis of epilepsy and few of them can be applied to predict the response to antiseizure medications (ASMs). In this study, we identified 26 significantly up-regulated genes and 32 down-regulated genes by comparing the gene expression profiles of patients with epilepsy that responded to valproate with those without applying any ASM. The results of gene set enrichment analysis indicated that the ferroptosis pathway was significantly impacted (p = 0.0087) in patients who responded to valproate. Interestingly, the gene NCOA4 in this pathway exhibited significantly different expression levels between the two groups, indicating that NCOA4 could serve as a potential biomarker to better understand the mechanism of valproate resistance. In addition, six up-regulated genes SF3A2, HMGN2, PABPN1, SSBP3, EFTUD2, and CREB3L2 as well as six down-regulated genes ZFP36L1, ACRC, SUB1, CALM2, TLK1, and STX2 also showed significantly different expression patterns between the two groups. Moreover, based on the gene expression profiles of the patients with the treatment of valproate, carbamazepine, and phenytoin, we proposed a strategy for predicting the response to the ASMs by using the Connectivity Map scoring method. Our findings could be helpful for better understanding the mechanisms of drug resistance of ASMs and improving the clinical treatment of epilepsy.

Unraveling Structural Alerts in Marketed Drugs for Improving Adverse Outcome Pathway Framework of Drug-Induced QT Prolongation.

In pharmaceutical treatment, many non-cardiac drugs carry the risk of prolonging the QT interval, which can lead to fatal cardiac complications such as torsades de points (TdP). Although the unexpected blockade of ion channels has been widely considered to be one of the main reasons for affecting the repolarization phase of the cardiac action potential and leading to QT interval prolongation, the lack of knowledge regarding chemical structures in drugs that may induce the prolongation of the QT interval remains a barrier to further understanding the underlying mechanism and developing an effective prediction strategy. In this study, we thoroughly investigated the differences in chemical structures between QT-prolonging drugs and drugs with no drug-induced QT prolongation (DIQT) concerns, based on the Drug-Induced QT Prolongation Atlas (DIQTA) dataset. Three categories of structural alerts (SAs), namely amines, ethers, and aromatic compounds, appeared in large quantities in QT-prolonging drugs, but rarely in drugs with no DIQT concerns, indicating a close association between SAs and the risk of DIQT. Moreover, using the molecular descriptors associated with these three categories of SAs as features, the structure-activity relationship (SAR) model for predicting the high risk of inducing QT interval prolongation of marketed drugs achieved recall rates of 72.5% and 80.0% for the DIQTA dataset and the FDA Adverse Event Reporting System (FAERS) dataset, respectively. Our findings may promote a better understanding of the mechanism of DIQT and facilitate research on cardiac adverse drug reactions in drug development.

Quantitative Structure-Activity Relationship (QSAR) Model for the Severity Prediction of Drug-Induced Rhabdomyolysis by Using Random Forest.

Drug-induced rhabdomyolysis (DIR) is a rare and potentially life-threatening muscle injury that is characterized by low incidence and high risk. To our best knowledge, the performance of the current predictive models for the early detection of DIR is suboptimal because of the scarcity and dispersion of DIR cases. Therefore, on the basis of the curated drug information from the Drug-Induced Rhabdomyolysis Atlas (DIRA) database, we proposed a random forest (RF) model to predict the DIR severity of the marketed drugs. Compared with the state-of-art methods, our proposed model outperformed extreme gradient boosting, support vector machine, and logistic regression in distinguishing the Most-DIR concern drugs from the No-DIR concern drugs (Matthews correlation coefficient (MCC) and recall rate of our model were 0.46 and 0.81, respectively). Our model was subsequently applied to predicting the potentially serious DIR for 1402 drugs, which were reported to cause DIR by the postmarketing DIR surveillance data in the FDA Spontaneous Adverse Events Reporting System (FAERS). As a result, 62.7% (94) of drugs ranked in the top 150 drugs with the Most-DIR concerns in FAERS can be identified by our model. The top four drugs (odds ratio >30) including acepromazine, rapacuronium, oxyphenbutazone, and naringenin were correctly predicted by our model. In conclusion, the RF model can well predict the Most-DIR concern drug only based on the chemical structure information and can be a facilitated tool for early DIR detection.

Uncovering the prognostic gene signatures for the improvement of risk stratification in cancers by using deep learning algorithm coupled with wavelet transform.

BACKGROUND: The aim of gene expression-based clinical modelling in tumorigenesis is not only to accurately predict the clinical endpoints, but also to reveal the genome characteristics for downstream analysis for the purpose of understanding the mechanisms of cancers. Most of the conventional machine learning methods involved a gene filtering step, in which tens of thousands of genes were firstly filtered based on the gene expression levels by a statistical method with an arbitrary cutoff. Although gene filtering procedure helps to reduce the feature dimension and avoid overfitting, there is a risk that some pathogenic genes important to the disease will be ignored. RESULTS: In this study, we proposed a novel deep learning approach by combining a convolutional neural network with stationary wavelet transform (SWT-CNN) for stratifying cancer patients and predicting their clinical outcomes without gene filtering based on tumor genomic profiles. The proposed SWT-CNN overperformed the state-of-art algorithms, including support vector machine (SVM) and logistic regression (LR), and produced comparable prediction performance to random forest (RF). Furthermore, for all the cancer types, we firstly proposed a method to weight the genes with the scores, which took advantage of the representative features in the hidden layer of convolutional neural network, and then selected the prognostic genes for the Cox proportional-hazards regression. The results showed that risk stratifications can be effectively improved by using the identified prognostic genes as feature, indicating that the representative features generated by SWT-CNN can well correlate the genes with prognostic risk in cancers and be helpful for selecting the prognostic gene signatures. CONCLUSIONS: Our results indicated that gene expression-based SWT-CNN model can be an excellent tool for stratifying the prognostic risk for cancer patients. In addition, the representative features of SWT-CNN were validated to be useful for evaluating the importance of the genes in the risk stratification and can be further used to identify the prognostic gene signatures.

Mutation status coupled with RNA-sequencing data can efficiently identify important non-significantly mutated genes serving as diagnostic biomarkers of endometrial cancer.

BACKGROUND: Endometrial cancers (ECs) are one of the most common types of malignant tumor in females. Substantial efforts had been made to identify significantly mutated genes (SMGs) in ECs and use them as biomarkers for the classification of histological subtypes and the prediction of clinical outcomes. However, the impact of non-significantly mutated genes (non-SMGs), which may also play important roles in the prognosis of EC patients, has not been extensively studied. Therefore, it is essential for the discovery of biomarkers in ECs to further investigate the non-SMGs that were highly associated with clinical outcomes. RESULTS: For the 9681 non-SMGs reported by the mutation annotation pipeline, there were 1053, 1273 and 395 non-SMGs differentially expressed between the patient groups divided by the clinical endpoints of histological grade, histological type as well as the International Federation of Gynecology and Obstetrics (FIGO) stage of ECs, respectively. In the gene set enrichment analysis, the cancer-related pathways, namely neuroactive ligand-receptor interaction signaling pathway, cAMP signaling pathway and calcium signaling pathway, were significantly enriched with the differentially expressed non-SMGs for all the three endpoints. We further identified 23, 19 and 24 non-SMGs, which were highly associated with histological grade, histological type and FIGO stage, respectively, from the differentially expressed non-SMGs by using the variable combination population analysis (VCPA) approach and found that 69.6% (16/23), 78.9% (15/19) and 66.7% (16/24) of the identified non-SMGs had been previously reported to be correlated with cancers. In addition, the averaged areas under the receiver operating characteristic curve (AUCs) achieved by the predictive models with identified non-SMGs as predictors in predicting histological type, histological grade, and FIGO stage were 0.993, 0.961 and 0.832, respectively, which were superior to those achieved by the models with SMGs as features (averaged AUCs = 0.928, 0.864 and 0.535, resp.). CONCLUSIONS: Besides the SMGs, the non-SMGs reported in the mutation annotation analysis may also involve the crucial genes that were highly associated with clinical outcomes. Combining the mutation status with the gene expression profiles can efficiently identify the cancer-related non-SMGs as predictors for cancer prognostic prediction and provide more supplemental candidates for the discovery of biomarkers.

Expression dynamics and relations with nearby genes of rat transposable elements across 11 organs, 4 developmental stages and both sexes.

BACKGROUND: TEs pervade mammalian genomes. However, compared with mice, fewer studies have focused on the TE expression patterns in rat, particularly the comparisons across different organs, developmental stages and sexes. In addition, TEs can influence the expression of nearby genes. The temporal and spatial influences of TEs remain unclear yet. RESULTS: To evaluate the TEs transcription patterns, we profiled their transcript levels in 11 organs for both sexes across four developmental stages of rat. The results show that most short interspersed elements (SINEs) are commonly expressed in all conditions, which are also the major TE types with commonly expression patterns. In contrast, long terminal repeats (LTRs) are more likely to exhibit specific expression patterns. The expression tendency of TEs and genes are similar in most cases. For example, few specific genes and TEs are in the liver, muscle and heart. However, TEs perform superior over genes on classing organ, which imply their higher organ specificity than genes. By associating the TEs with the closest genes in genome, we find their expression levels are correlated, independent of their distance in some cases. CONCLUSIONS: TEs sex-dependently associate with nearest genes. A gene would be associated with more than one TE. Our works can help to functionally annotate the genome and further understand the role of TEs in gene regulation.

Label-free blood serum detection by using surface-enhanced Raman spectroscopy and support vector machine for the preoperative diagnosis of parotid gland tumors.

BACKGROUND: It is difficult for the parotid gland neoplasms to make an accurate preoperative diagnosis due to the restriction of biopsy in the parotid gland neoplasms. The aim of this study is to apply the surface-enhanced Raman spectroscopy (SERS) method for the blood serum biochemical detection and use the support vector machine for the analysis in order to develop a simple but accurate blood serum detection for preoperative diagnosis of the parotid gland neoplasms. METHODS: The blood serums were collected from four groups: the patients with pleomorphic adenoma, the patients with Warthin's tumor, the patients with mucoepidermoid carcinoma and the volunteers without parotid gland neoplasms. Au nanoparticles (Au NPs) were mixed with the blood serum as the SERS active nanosensor to enhance the Raman scattering signals produced by the various biochemical materials and high quality SERS spectrum were obtained by using the Raman microscope system. Then the support vector machine was utilized to analyze the differences of the SERS spectrum from the blood serum of different groups and established a diagnostic model to discriminate the different groups. RESULTS: It was demonstrated that there were different intensities of SERS peaks assigned to various biochemical changes in the blood serum between the parotid gland tumor groups and normal control group. Compared with the SERS spectra of the normal serums, the intensities of peaks assigned to nucleic acids and proteins increased in the SERS spectra of the parotid gland tumor serums, which manifested the differences of the biochemical metabolites in the serum from the patients with parotid gland tumors. When the leave-one-sample-out method was used, support vector machine (SVM) played an outstanding performance in the classification of the SERS spectra with the high accuracy (84.1 % ~ 88.3 %), sensitivity (82.2 % ~ 97.4 %) and specificity (73.7 % ~ 86.7 %). Though the accuracy, sensitivity and specificity decreased in the leave-one-patient-out cross validation, the mucoepidermoid carcinoma was still easier to diagnose than other tumors. DISCUSSION: The specific molecular differences of parotid gland tumors and normal serums were significantly demonstrated through the comparison between the various SERS spectra.But compared with the serum SERS spectra reported in the other studies, some differences exist between the spectra in this study and the ones reported in the lietratures. These differences may result from the various nano-particles, the different preparation of serum and equipment parameters, and we could need a further research to find an exact explanation.Based on the SERS spectra of the serum samples, SVM have shown a giant potential to diagnose the parotid gland tumors in our preliminary study. However, different cross validaiton methods could effect the accuracy and a further study involing a great number of samples should be needed. CONCLUSIONS: This exploratory research demonstrated the great potential of SERS combined with SVM into a non-invasive clinical diagnostic method for preoperative diagnosis of parotid gland tumors. And the internal relation between the spectra and patients should be established in the further study.

Genome-wide association study combined with biological context can reveal more disease-related SNPs altering microRNA target seed sites.

BACKGROUND: Emerging studies demonstrate that single nucleotide polymorphisms (SNPs) resided in the microRNA recognition element seed sites (MRESSs) in 3'UTR of mRNAs are putative biomarkers for human diseases and cancers. However, exhaustively experimental validation for the causality of MRESS SNPs is impractical. Therefore bioinformatics have been introduced to predict causal MRESS SNPs. Genome-wide association study (GWAS) provides a way to detect susceptibility of millions of SNPs simultaneously by taking linkage disequilibrium (LD) into account, but the multiple-testing corrections implemented to suppress false positive rate always sacrificed the sensitivity. In our study, we proposed a method to identify candidate causal MRESS SNPs from 12 GWAS datasets without performing multiple-testing corrections. Alternatively, we used biological context to ensure credibility of the selected SNPs. RESULTS: In 11 out of the 12 GWAS datasets, MRESS SNPs were over-represented in SNPs with p-value ≤ 0.05 (odds ratio (OR) ranged from 1.1 to 2.4). Moreover, host genes of susceptible MRESS SNPs in each of the 11 GWAS dataset shared biological context with reported causal genes. There were 286 MRESS SNPs identified by our method, while only 13 SNPs were identified by multiple-testing corrections with a given threshold of 1 × 10-5, which is a common cutoff used in GWAS. 27 out of the 286 candidate SNPs have been reported to be deleterious while only 2 out of 13 multiple-testing corrected SNPs were documented in PubMed. MicroRNA-mRNA interactions affected by the 286 candidate SNPs were likely to present negatively correlated expression. These SNPs introduced greater alternation of binding free energy than other MRESS SNPs, especially when grouping by haplotypes (4210 vs. 4105 cal/mol by mean, 9781 vs. 8521 cal/mol by mean, respectively). CONCLUSIONS: MRESS SNPs are promising disease biomarkers in multiple GWAS datasets. The method of integrating GWAS p-value and biological context is stable and effective for selecting candidate causal MRESS SNPs, it reduces the loss of sensitivity compared to multiple-testing corrections. The 286 candidate causal MRESS SNPs provide researchers a credible source to initialize their design of experimental validations in the future.

Extend the benchmarking indel set by manual review using the individual cell line sequencing data from the Sequencing Quality Control 2 (SEQC2) project.

Accurate indel calling plays an important role in precision medicine. A benchmarking indel set is essential for thoroughly evaluating the indel calling performance of bioinformatics pipelines. A reference sample with a set of known-positive variants was developed in the FDA-led Sequencing Quality Control Phase 2 (SEQC2) project, but the known indels in the known-positive set were limited. This project sought to provide an enriched set of known indels that would be more translationally relevant by focusing on additional cancer related regions. A thorough manual review process completed by 42 reviewers, two advisors, and a judging panel of three researchers significantly enriched the known indel set by an additional 516 indels. The extended benchmarking indel set has a large range of variant allele frequencies (VAFs), with 87% of them having a VAF below 20% in reference Sample A. The reference Sample A and the indel set can be used for comprehensive benchmarking of indel calling across a wider range of VAF values in the lower range. Indel length was also variable, but the majority were under 10 base pairs (bps). Most of the indels were within coding regions, with the remainder in the gene regulatory regions. Although high confidence can be derived from the robust study design and meticulous human review, this extensive indel set has not undergone orthogonal validation. The extended benchmarking indel set, along with the indels in the previously published known-positive set, was the truth set used to benchmark indel calling pipelines in a community challenge hosted on the precisionFDA platform. This benchmarking indel set and reference samples can be utilized for a comprehensive evaluation of indel calling pipelines. Additionally, the insights and solutions obtained during the manual review process can aid in improving the performance of these pipelines.

Investigating the concordance of Gene Ontology terms reveals the intra- and inter-platform reproducibility of enrichment analysis.

BACKGROUND: Reliability and Reproducibility of differentially expressed genes (DEGs) are essential for the biological interpretation of microarray data. The microarray quality control (MAQC) project launched by US Food and Drug Administration (FDA) elucidated that the lists of DEGs generated by intra- and inter-platform comparisons can reach a high level of concordance, which mainly depended on the statistical criteria used for ranking and selecting DEGs. Generally, it will produce reproducible lists of DEGs when combining fold change ranking with a non-stringent p-value cutoff. For further interpretation of the gene expression data, statistical methods of gene enrichment analysis provide powerful tools for associating the DEGs with prior biological knowledge, e.g. Gene Ontology (GO) terms and pathways, and are widely used in genome-wide research. Although the DEG lists generated from the same compared conditions proved to be reliable, the reproducible enrichment results are still crucial to the discovery of the underlying molecular mechanism differentiating the two conditions. Therefore, it is important to know whether the enrichment results are still reproducible, when using the lists of DEGs generated by different statistic criteria from inter-laboratory and cross-platform comparisons. In our study, we used the MAQC data sets for systematically accessing the intra- and inter-platform concordance of GO terms enriched by Gene Set Enrichment Analysis (GSEA) and LRpath. RESULTS: In intra-platform comparisons, the overlapped percentage of enriched GO terms was as high as ~80% when the inputted lists of DEGs were generated by fold change ranking and Significance Analysis of Microarrays (SAM), whereas the percentages decreased about 20% when generating the lists of DEGs by using fold change ranking and t-test, or by using SAM and t-test. Similar results were found in inter-platform comparisons. CONCLUSIONS: Our results demonstrated that the lists of DEGs in a high level of concordance can ensure the high concordance of enrichment results. Importantly, based on the lists of DEGs generated by a straightforward method of combining fold change ranking with a non-stringent p-value cutoff, enrichment analysis will produce reproducible enriched GO terms for the biological interpretation.

Transcriptional profiling of Chinese medicinal formula Si-Wu-Tang on breast cancer cells reveals phytoestrogenic activity.

BACKGROUND: Si-Wu-Tang (SWT), comprising the combination of four herbs, Paeoniae, Angelicae, Chuanxiong and Rehmanniae, is one of the most popular traditional oriental medicines for women's diseases. In our previous study, the microarray gene expression profiles of SWT on breast cancer cell line MCF-7 were found similar to the effect of ß-estradiol (E2) on MCF-7 cells in the Connectivity Map database. METHODS: Further data analysis was conducted to find the main similarities and differences between the effects of SWT and E2 on MCF-7 gene expression. The cell proliferation assay on MCF-7 (ER-positive) and MDA-MB-231 (ER-negative) cells were used to examine such estrogenic activity. The estrogenic potency of SWT was further confirmed by estrogen-responsive element (ERE) luciferase reporter assay in MCF-7 cells. RESULTS: Many estrogen regulated genes strongly up-regulated by E2 were similarly up-regulated by SWT, e.g., GREB1, PGR and EGR3. Of interest with regard to safety of SWT, the oncogenes MYBL1 and RET were strongly induced by E2 but not by SWT. Quantitative RT-PCR analysis revealed a highly concordant expression change in selected genes with data obtained by microarrays. Further supporting SWT's estrogenic activity, in MCF-7 but not in MDA-MB-231 cells, SWT stimulated cell growth at lower concentrations (< 3.0 mg/ml), while at high concentrations, it inhibits the growth of both cell lines. The growth inhibitory potency of SWT was significantly higher in MDA-MB-231 than in MCF-7 cells. The SWT-induced cell growth of MCF-7 could be blocked by addition of the estrogen receptor antagonist tamoxifen. In addition, SWT was able to activate the ERE activity at lower concentrations. The herbal components Angelicae, Chuanxiong and Rehmanniae at lower concentrations (< 3.0 mg/ml) also showed growth-inducing and ERE-activating activity in MCF-7 cells. CONCLUSIONS: These results revealed a new mechanism to support the clinical use of SWT for estrogen related diseases and possibly for cancer prevention. This study also demonstrated the feasibility of using microarray transcriptional profiling to discover phytoestrogenic components that are present in natural products.

Differentiation of epilepsy and psychogenic nonepileptic events based on body fluid characteristics.

OBJECTIVE: Differential diagnosis between epileptic seizures and psychogenic nonepileptic events (PNEEs) is a worldwide problem for neurologists. The present study aims to identify important characteristics from body fluid tests and develop diagnostic models based on them. METHODS: This is a register-based observational study in patients with a diagnosis of epilepsy or PNEEs at West China Hospital of Sichuan University. Data from body fluid tests between 2009 and 2019 were used as a training set. We constructed models with a random forest approach in eight training subsets divided by sex and categories of tests, including electrolyte, blood cell, metabolism, and urine tests. Then, we collected data prospectively from patients between 2020 and 2022 to validate our models and calculated the relative importance of characteristics in robust models. Selected characteristics were finally analyzed with multiple logistic regression to establish nomograms. RESULTS: A total of 388 patients, including 218 with epilepsy and 170 with PNEEs, were studied. The AUROCs of random forest models of electrolyte and urine tests in the validation phase achieved 80.0% and 79.0%, respectively. Carbon dioxide combining power, anion gap, potassium, calcium, and chlorine in electrolyte tests and specific gravity, pH, and conductivity in urine tests were selected for the logistic regression analysis. C (ROC) of the electrolyte and urine diagnostic nomograms achieved 0.79 and 0.85, respectively. SIGNIFICANCE: The application of routine indicators of serum and urine may help in the more accurate identification of epileptic and PNEEs.

Wavelet transform analysis of NMR structure ensembles to reveal internal fluctuations of enzymes.

Internal motions and flexibility are essential for biological functions in proteins. To assess the internal fluctuations and conformational flexibility of proteins, reliable computational methods are needed. In this study, wavelet transformation was used to filter out the noise and facilitate investigating the internal positional fluctuations of enzymes within nuclear magnetic resonance (NMR) structure ensembles. Moreover, potential active sites were identified by combining with positional fluctuation score, sequence conservation, and solvent accessible surface area. Among the total 107 catalytic residues in 44 examined enzymes, 69 residues were identified correctly. Our results suggest that wavelet transform analysis of structure ensemble is applicable to extract essential fluctuation information of proteins; furthermore, analysis of positional fluctuations is helpful for the identification of catalytic residues.

Drug-induced QT Prolongation Atlas (DIQTA) for enhancing cardiotoxicity management.

Drug-induced prolongation of the QT interval is common in a variety of pharmaceutical treatments and can lead to serious clinical outcomes. Although substantial efforts have been made to prevent drug-induced QT interval prolongation, the lack of a centralized data source remains the main obstacle to further study of the underlying mechanism and the development of effective prediction strategies. To fill this gap, we propose a schema for stratifying the risk of marketed QT prolonging drugs based on US Food and Drug Administration (FDA)-approved drug labeling and developed a Drug-Induced QT Prolongation Atlas (DIQTA). Potential application of DIQTA was shown by precision dosing in off-label use and therapeutic strategy optimization, as well as the facilitation of artificial intelligence (AI)-based modeling in predictive toxicity.

Raman Spectroscopy: A Potential Diagnostic Tool for Oral Diseases.

Oral diseases impose a major health burden worldwide and have a profound effect on general health. Dental caries, periodontal diseases, and oral cancers are the most common oral health conditions. Their occurrence and development are related to oral microbes, and effective measures for their prevention and the promotion of oral health are urgently needed. Raman spectroscopy detects molecular vibration information by collecting inelastic scattering light, allowing a "fingerprint" of a sample to be acquired. It provides the advantages of rapid, sensitive, accurate, and minimally invasive detection as well as minimal interference from water in the "fingerprint region." Owing to these characteristics, Raman spectroscopy has been used in medical detection in various fields to assist diagnosis and evaluate prognosis, such as detecting and differentiating between bacteria or between neoplastic and normal brain tissues. Many oral diseases are related to oral microbial dysbiosis, and their lesions differ from normal tissues in essential components. The colonization of keystone pathogens, such as Porphyromonas gingivalis, resulting in microbial dysbiosis in subgingival plaque, is the main cause of periodontitis. Moreover, the components in gingival crevicular fluid, such as infiltrating inflammatory cells and tissue degradation products, are markedly different between individuals with and without periodontitis. Regarding dental caries, the compositions of decayed teeth are transformed, accompanied by an increase in acid-producing bacteria. In oral cancers, the compositions and structures of lesions and normal tissues are different. Thus, the changes in bacteria and the components of saliva and tissue can be used in examinations as special markers for these oral diseases, and Raman spectroscopy has been acknowledged as a promising measure for detecting these markers. This review summarizes and discusses key research and remaining problems in this area. Based on this, suggestions for further study are proposed.

Structure-Activity Relationship (SAR) Model for Predicting Teratogenic Risk of Antiseizure Medications in Pregnancy by Using Support Vector Machine.

Teratogenicity is one of the main concerns in clinical medications of pregnant women. Prescription of antiseizure medications (ASMs) in women with epilepsy during pregnancy may cause teratogenic effects on the fetus. Although large scale epilepsy pregnancy registries played an important role in evaluating the teratogenic risk of ASMs, for most ASMs, especially the newly approved ones, the potential teratogenic risk cannot be effectively assessed due to the lack of evidence. In this study, the analyses are performed on any medication, with a focus on ASMs. We curated a list containing the drugs with potential teratogenicity based on the US Food and Drug Administration (FDA)-approved drug labeling, and established a support vector machine (SVM) model for detecting drugs with high teratogenic risk. The model was validated by using the post-marketing surveillance data from US FDA Spontaneous Adverse Events Reporting System (FAERS) and applied to the prediction of potential teratogenic risk of ASMs. Our results showed that our proposed model outperformed the state-of-art approaches, including logistic regression (LR), random forest (RF) and extreme gradient boosting (XGBoost), when detecting the high teratogenic risk of drugs (MCC and recall rate were 0.312 and 0.851, respectively). Among 196 drugs with teratogenic potential reported by FAERS, 136 (69.4%) drugs were correctly predicted. For the eight commonly used ASMs, 4 of them were predicted as high teratogenic risk drugs, including topiramate, phenobarbital, valproate and phenytoin (predicted probabilities of teratogenic risk were 0.69, 0.60 0.59, and 0.56, respectively), which were consistent with the statement in FDA-approved drug labeling and the high reported prevalence of teratogenicity in epilepsy pregnancy registries. In addition, the structural alerts in ASMs that related to the genotoxic carcinogenicity and mutagenicity, idiosyncratic adverse reaction, potential electrophilic agents and endocrine disruption were identified and discussed. Our findings can be a good complementary for the teratogenic risk assessment in drug development and facilitate the determination of pharmacological therapies during pregnancy.

Predicting disease-associated substitution of a single amino acid by analyzing residue interactions.

BACKGROUND: The rapid accumulation of data on non-synonymous single nucleotide polymorphisms (nsSNPs, also called SAPs) should allow us to further our understanding of the underlying disease-associated mechanisms. Here, we use complex networks to study the role of an amino acid in both local and global structures and determine the extent to which disease-associated and polymorphic SAPs differ in terms of their interactions to other residues. RESULTS: We found that SAPs can be well characterized by network topological features. Mutations are probably disease-associated when they occur at a site with a high centrality value and/or high degree value in a protein structure network. We also discovered that study of the neighboring residues around a mutation site can help to determine whether the mutation is disease-related or not. We compiled a dataset from the Swiss-Prot variant pages and constructed a model to predict disease-associated SAPs based on the random forest algorithm. The values of total accuracy and MCC were 83.0% and 0.64, respectively, as determined by 5-fold cross-validation. With an independent dataset, our model achieved a total accuracy of 80.8% and MCC of 0.59, respectively. CONCLUSIONS: The satisfactory performance suggests that network topological features can be used as quantification measures to determine the importance of a site on a protein, and this approach can complement existing methods for prediction of disease-associated SAPs. Moreover, the use of this method in SAP studies would help to determine the underlying linkage between SAPs and diseases through extensive investigation of mutual interactions between residues.

Comparative Analysis for the Performance of Long-Read-Based Structural Variation Detection Pipelines in Tandem Repeat Regions.

There has been growing recognition of the vital links between structural variations (SVs) and diverse diseases. Research suggests that, with much longer DNA fragments and abundant contextual information, long-read technologies have advantages in SV detection even in complex repetitive regions. So far, several pipelines for calling SVs from long-read sequencing data have been proposed and used in human genome research. However, the performance of these pipelines is still lack of deep exploration and adequate comparison. In this study, we comprehensively evaluated the performance of three commonly used long-read SV detection pipelines, namely PBSV, Sniffles and PBHoney, especially the performance on detecting the SVs in tandem repeat regions (TRRs). Evaluated by using a robust benchmark for germline SV detection as the gold standard, we thoroughly estimated the precision, recall and F1 score of insertions and deletions detected by the pipelines. Our results revealed that all these pipelines clearly exhibited better performance outside TRRs than that in TRRs. The F1 scores of Sniffles in and outside TRRs were 0.60 and 0.76, respectively. The performance of PBSV was similar to that of Sniffles, and was generally higher than that of PBHoney. In conclusion, our findings can be benefit for choosing the appropriate pipelines in real practice and are good complementary to the application of long-read sequencing technologies in the research of rare diseases.

Prediction of Synergistic Drug Combinations for Prostate Cancer by Transcriptomic and Network Characteristics.

Prostate cancer (PRAD) is a major cause of cancer-related deaths. Current monotherapies show limited efficacy due to often rapidly emerging resistance. Combination therapies could provide an alternative solution to address this problem with enhanced therapeutic effect, reduced cytotoxicity, and delayed the appearance of drug resistance. However, it is prohibitively cost and labor-intensive for the experimental approaches to pick out synergistic combinations from the millions of possibilities. Thus, it is highly desired to explore other efficient strategies to assist experimental researches. Inspired by the challenge, we construct the transcriptomics-based and network-based prediction models to quickly screen the potential drug combination for Prostate cancer, and further assess their performance by in vitro assays. The transcriptomics-based method screens nine possible combinations. However, the network-based method gives discrepancies for at least three drug pairs. Further experimental results indicate the dose-dependent effects of the three docetaxel-containing combinations, and confirm the synergistic effects of the other six combinations predicted by the transcriptomics-based model. For the network-based predictions, in vitro tests give opposite results to the two combinations (i.e. mitoxantrone-cyproheptadine and cabazitaxel-cyproheptadine). Namely, the transcriptomics-based method outperforms the network-based one for the specific disease like Prostate cancer, which provide guideline for selection of the computational methods in the drug combination screening. More importantly, six combinations (the three mitoxantrone-containing and the three cabazitaxel-containing combinations) are found to be promising candidates to synergistically conquer Prostate cancer.