INTEDE: interactome of drug-metabolizing enzymes

Integrative Analysis of Histone Acetylation Regulated CYP4F12 in Esophageal Cancer Development

Current therapeutic strategies for esophageal cancer (EC) patients have yielded limited improvements in survival rates. Recent research has highlighted the influence of drug metabolism enzymes on both drug response and EC development. Our study aims to identify specific drug metabolism enzymes regulated by histone acetylation and to elucidate its molecular and clinical features. CYP4F12 exhibited a notable upregulation subsequent to trichostatin A treatment as evidenced by RNA sequencing analysis conducted on the KYSE-150 cell line. The change in gene expression was associated with increased acetylation level of histone 3 K18 and K27 in the promoter. The regulation was dependent on p300. In silicon analysis of both The Cancer Genome Atlas esophageal carcinoma and GSE53624 dataset suggested a critical role of CYP4F12 in EC development, because CYP4F12 was downregulated in tumor tissues and predicted better disease-free survival. Gene ontology analysis has uncovered a robust correlation between CYP4F12 and processes related to cell migration, as well as its involvement in cytosine-mediated immune activities. Further investigation into the relationship between immune cells and CYP4F12 expression has indicated an increased level of B cell infiltration in samples with high CYP4F12 expression. CYP4F12 was also negatively correlated with the expression of inhibitory checkpoints. An accurate predictive nomogram model was established combining with clinical factors and CYP4F12 expression. In conclusion, CYP4F12 was crucial in EC development, and targeting CYP4F12 may improve the therapeutic efficacy of current treatment in EC patients. SIGNIFICANCE STATEMENT: CYP4F12 expression was downregulated in esophageal cancer (EC) patients and could be induced by trichostatin A. During EC development, CYP4F12 was linked to reduced cell migration and increased infiltration of B cells. CYP4F12 also is a biomarker as prognostic predictors and therapeutic guide in EC patients.

Application of artificial intelligence in drug design: A review

Artificial intelligence (AI) is a field of computer science that involves acquiring information, developing rule bases, and mimicking human behaviour. The fundamental concept behind AI is to create intelligent computer systems that can operate with minimal human intervention or without any intervention at all. These rule-based systems are developed using various machine learning and deep learning models, enabling them to solve complex problems. AI is integrated with these models to learn, understand, and analyse provided data. The rapid advancement of Artificial Intelligence (AI) is reshaping numerous industries, with the pharmaceutical sector experiencing a notable transformation. AI is increasingly being employed to automate, optimize, and personalize various facets of the pharmaceutical industry, particularly in pharmacological research. Traditional drug development methods areknown for being time-consuming, expensive, and less efficient, often taking around a decade and costing billions of dollars. The integration of artificial intelligence (AI) techniques addresses these challenges by enabling the examination of compounds with desired properties from a vast pool of input drugs. Furthermore, it plays a crucial role in drug screening by predicting toxicity, bioactivity, ADME properties (absorption, distribution, metabolism, and excretion), physicochemical properties, and more. AI enhances the drug design process by improving the efficiency and accuracy of predicting drug behaviour, interactions, and properties. These approaches further significantly improve the precision of drug discovery processes and decrease clinical trial costs leading to the development of more effective drugs.

Pan-cancer single-cell landscape of drug-metabolizing enzyme genes

OBJECTIVE

Varied expression of drug-metabolizing enzymes (DME) genes dictates the intensity and duration of drug response in cancer treatment. This study aimed to investigate the transcriptional profile of DMEs in tumor microenvironment (TME) at single-cell level and their impact on individual responses to anticancer therapy.

METHODS

Over 1.3 million cells from 481 normal/tumor samples across 9 solid cancer types were integrated to profile changes in the expression of DME genes. A ridge regression model based on the PRISM database was constructed to predict the influence of DME gene expression on drug sensitivity.

RESULTS

Distinct expression patterns of DME genes were revealed at single-cell resolution across different cancer types. Several DME genes were highly enriched in epithelial cells (e.g. GPX2, TST and CYP3A5) or different TME components (e.g. CYP4F3 in monocytes). Particularly, GPX2 and TST were differentially expressed in epithelial cells from tumor samples compared to those from normal samples. Utilizing the PRISM database, we found that elevated expression of GPX2, CYP3A5 and reduced expression of TST was linked to enhanced sensitivity of particular chemo-drugs (e.g. gemcitabine, daunorubicin, dasatinib, vincristine, paclitaxel and oxaliplatin).

CONCLUSION

Our findings underscore the varied expression pattern of DME genes in cancer cells and TME components, highlighting their potential as biomarkers for selecting appropriate chemotherapy agents.

Exploratory drug discovery in breast cancer patients: A multimodal deep learning approach to identify novel drug candidates targeting RTK signaling

Breast cancer, a highly formidable and diverse malignancy predominantly affecting women globally, poses a significant threat due to its intricate genetic variability, rendering it challenging to diagnose accurately. Various therapies such as immunotherapy, radiotherapy, and diverse chemotherapy approaches like drug repurposing and combination therapy are widely used depending on cancer subtype and metastasis severity. Our study revolves around an innovative drug discovery strategy targeting potential drug candidates specific to RTK signalling, a prominently targeted receptor class in cancer. To accomplish this, we have developed a multimodal deep neural network (MM-DNN) based QSAR model integrating omics datasets to elucidate genomic, proteomic expression data, and drug responses, validated rigorously. The results showcase an R2 value of 0.917 and an RMSE value of 0.312, affirming the model's commendable predictive capabilities. Structural analogs of drug molecules specific to RTK signalling were sourced from the PubChem database, followed by meticulous screening to eliminate dissimilar compounds. Leveraging the MM-DNN-based QSAR model, we predicted the biological activity of these molecules, subsequently clustering them into three distinct groups. Feature importance analysis was performed. Consequently, we successfully identified prime drug candidates tailored for each potential downstream regulatory protein within the RTK signalling pathway. This method makes the early stages of drug development faster by removing inactive compounds, providing a hopeful path in combating breast cancer.

DDID: a comprehensive resource for visualization and analysis of diet-drug interactions

Diet-drug interactions (DDIs) are pivotal in drug discovery and pharmacovigilance. DDIs can modify the systemic bioavailability/pharmacokinetics of drugs, posing a threat to public health and patient safety. Therefore, it is crucial to establish a platform to reveal the correlation between diets and drugs. Accordingly, we have established a publicly accessible online platform, known as Diet-Drug Interactions Database (DDID, https://bddg.hznu.edu.cn/ddid/), to systematically detail the correlation and corresponding mechanisms of DDIs. The platform comprises 1338 foods/herbs, encompassing flora and fauna, alongside 1516 widely used drugs and 23 950 interaction records. All interactions are meticulously scrutinized and segmented into five categories, thereby resulting in evaluations (positive, negative, no effect, harmful and possible). Besides, cross-linkages between foods/herbs, drugs and other databases are furnished. In conclusion, DDID is a useful resource for comprehending the correlation between foods, herbs and drugs and holds a promise to enhance drug utilization and research on drug combinations.

FERREG: ferroptosis-based regulation of disease occurrence, progression and therapeutic response

Ferroptosis is a non-apoptotic, iron-dependent regulatory form of cell death characterized by the accumulation of intracellular reactive oxygen species. In recent years, a large and growing body of literature has investigated ferroptosis. Since ferroptosis is associated with various physiological activities and regulated by a variety of cellular metabolism and mitochondrial activity, ferroptosis has been closely related to the occurrence and development of many diseases, including cancer, aging, neurodegenerative diseases, ischemia-reperfusion injury and other pathological cell death. The regulation of ferroptosis mainly focuses on three pathways: system Xc-/GPX4 axis, lipid peroxidation and iron metabolism. The genes involved in these processes were divided into driver, suppressor and marker. Importantly, small molecules or drugs that mediate the expression of these genes are often good treatments in the clinic. Herein, a newly developed database, named 'FERREG', is documented to (i) providing the data of ferroptosis-related regulation of diseases occurrence, progression and drug response; (ii) explicitly describing the molecular mechanisms underlying each regulation; and (iii) fully referencing the collected data by cross-linking them to available databases. Collectively, FERREG contains 51 targets, 718 regulators, 445 ferroptosis-related drugs and 158 ferroptosis-related disease responses. FERREG can be accessed at https://idrblab.org/ferreg/.

Gut Microbiome Integration in Drug Discovery and Development of Small Molecules

Human microbiomes, particularly in the gut, could have a major impact on the efficacy and toxicity of drugs. However, gut microbial metabolism is often neglected in the drug discovery and development process. Medicen, a Paris-based human health innovation cluster, has gathered more than 30 international leading experts from pharma, academia, biotech, clinical research organizations, and regulatory science to develop proposals to facilitate the integration of microbiome science into drug discovery and development. Seven subteams were formed to cover the complementary expertise areas of 1) pharma experience and case studies, 2) in silico microbiome-drug interaction, 3) in vitro microbial stability screening, 4) gut fermentation models, 5) animal models, 6) microbiome integration in clinical and regulatory aspects, and 7) microbiome ecosystems and models. Each expert team produced a state-of-the-art report of their respective field highlighting existing microbiome-related tools at every stage of drug discovery and development. The most critical limitations are the growing, but still limited, drug-microbiome interaction data to produce predictive models and the lack of agreed-upon standards despite recent progress. In this paper we will report on and share proposals covering 1) how microbiome tools can support moving a compound from drug discovery to clinical proof-of-concept studies and alert early on potential undesired properties stemming from microbiome-induced drug metabolism and 2) how microbiome data can be generated and integrated in pharmacokinetic models that are predictive of the human situation. Examples of drugs metabolized by the microbiome will be discussed in detail to support recommendations from the working group. SIGNIFICANCE STATEMENT: Gut microbial metabolism is often neglected in the drug discovery and development process despite growing evidence of drugs' efficacy and safety impacted by their interaction with the microbiome. This paper will detail existing microbiome-related tools covering every stage of drug discovery and development, current progress, and limitations, as well as recommendations to integrate them into the drug discovery and development process.

Is fragment-based graph a better graph-based molecular representation for drug design? A comparison study of graph-based models

Graph Neural Networks (GNNs) have gained significant traction in various sectors of AI-driven drug design. Over recent years, the integration of fragmentation concepts into GNNs has emerged as a potent strategy to augment the efficacy of molecular generative models. Nonetheless, challenges such as symmetry breaking and potential misrepresentation of intricate cycles and undefined functional groups raise questions about the superiority of fragment-based graph representation over traditional methods. In our research, we undertook a rigorous evaluation, contrasting the predictive prowess of eight models-developed using deep learning algorithms-across 12 benchmark datasets that span a range of properties. These models encompass established methods like GCN, AttentiveFP, and D-MPNN, as well as innovative fragment-based representation techniques. Our results indicate that fragment-based methodologies, notably PharmHGT, significantly improve model performance and interpretability, particularly in scenarios characterized by limited data availability. However, in situations with extensive training, fragment-based molecular graph representations may not necessarily eclipse traditional methods. In summation, we posit that the integration of fragmentation, as an avant-garde technique in drug design, harbors considerable promise for the future of AI-enhanced drug design.

VARIDT 3.0: the phenotypic and regulatory variability of drug transporter

The phenotypic and regulatory variability of drug transporter (DT) are vital for the understanding of drug responses, drug-drug interactions, multidrug resistances, and so on. The ADME property of a drug is collectively determined by multiple types of variability, such as: microbiota influence (MBI), transcriptional regulation (TSR), epigenetics regulation (EGR), exogenous modulation (EGM) and post-translational modification (PTM). However, no database has yet been available to comprehensively describe these valuable variabilities of DTs. In this study, a major update of VARIDT was therefore conducted, which gave 2072 MBIs, 10 610 TSRs, 46 748 EGRs, 12 209 EGMs and 10 255 PTMs. These variability data were closely related to the transportation of 585 approved and 301 clinical trial drugs for treating 572 diseases. Moreover, the majority of the DTs in this database were found with multiple variabilities, which allowed a collective consideration in determining the ADME properties of a drug. All in all, VARIDT 3.0 is expected to be a popular data repository that could become an essential complement to existing pharmaceutical databases, and is freely accessible without any login requirement at: https://idrblab.org/varidt/.

TTD: Therapeutic Target Database describing target druggability information

Target discovery is one of the essential steps in modern drug development, and the identification of promising targets is fundamental for developing first-in-class drug. A variety of methods have emerged for target assessment based on druggability analysis, which refers to the likelihood of a target being effectively modulated by drug-like agents. In the therapeutic target database (TTD), nine categories of established druggability characteristics were thus collected for 426 successful, 1014 clinical trial, 212 preclinical/patented, and 1479 literature-reported targets via systematic review. These characteristic categories were classified into three distinct perspectives: molecular interaction/regulation, human system profile and cell-based expression variation. With the rapid progression of technology and concerted effort in drug discovery, TTD and other databases were highly expected to facilitate the explorations of druggability characteristics for the discovery and validation of innovative drug target. TTD is now freely accessible at: https://idrblab.org/ttd/.

TheMarker: a comprehensive database of therapeutic biomarkers

Distinct from the traditional diagnostic/prognostic biomarker (adopted as the indicator of disease state/process), the therapeutic biomarker (ThMAR) has emerged to be very crucial in the clinical development and clinical practice of all therapies. There are five types of ThMAR that have been found to play indispensable roles in various stages of drug discovery, such as: Pharmacodynamic Biomarker essential for guaranteeing the pharmacological effects of a therapy, Safety Biomarker critical for assessing the extent or likelihood of therapy-induced toxicity, Monitoring Biomarker indispensable for guiding clinical management by serially measuring patients' status, Predictive Biomarker crucial for maximizing the clinical outcome of a therapy for specific individuals, and Surrogate Endpoint fundamental for accelerating the approval of a therapy. However, these data of ThMARs has not been comprehensively described by any of the existing databases. Herein, a database, named 'TheMarker', was therefore constructed to (a) systematically offer all five types of ThMAR used at different stages of drug development, (b) comprehensively describe ThMAR information for the largest number of drugs among available databases, (c) extensively cover the widest disease classes by not just focusing on anticancer therapies. These data in TheMarker are expected to have great implication and significant impact on drug discovery and clinical practice, and it is freely accessible without any login requirement at: https://idrblab.org/themarker.

ADCdb: the database of antibody-drug conjugates

Antibody-drug conjugates (ADCs) are a class of innovative biopharmaceutical drugs, which, via their antibody (mAb) component, deliver and release their potent warhead (a.k.a. payload) at the disease site, thereby simultaneously improving the efficacy of delivered therapy and reducing its off-target toxicity. To design ADCs of promising efficacy, it is crucial to have the critical data of pharma-information and biological activities for each ADC. However, no such database has been constructed yet. In this study, a database named ADCdb focusing on providing ADC information (especially its pharma-information and biological activities) from multiple perspectives was thus developed. Particularly, a total of 6572 ADCs (359 approved by FDA or in clinical trial pipeline, 501 in preclinical test, 819 with in-vivo testing data, 1868 with cell line/target testing data, 3025 without in-vivo/cell line/target testing data) together with their explicit pharma-information was collected and provided. Moreover, a total of 9171 literature-reported activities were discovered, which were identified from diverse clinical trial pipelines, model organisms, patient/cell-derived xenograft models, etc. Due to the significance of ADCs and their relevant data, this new database was expected to attract broad interests from diverse research fields of current biopharmaceutical drug discovery. The ADCdb is now publicly accessible at: https://idrblab.org/adcdb/.

INTEDE 2.0: the metabolic roadmap of drugs

The metabolic roadmap of drugs (MRD) is a comprehensive atlas for understanding the stepwise and sequential metabolism of certain drug in living organisms. It plays a vital role in lead optimization, personalized medication, and ADMET research. The MRD consists of three main components: (i) the sequential catalyses of drug and its metabolites by different drug-metabolizing enzymes (DMEs), (ii) a comprehensive collection of metabolic reactions along the entire MRD and (iii) a systematic description on efficacy & toxicity for all metabolites of a studied drug. However, there is no database available for describing the comprehensive metabolic roadmaps of drugs. Therefore, in this study, a major update of INTEDE was conducted, which provided the stepwise & sequential metabolic roadmaps for a total of 4701 drugs, and a total of 22 165 metabolic reactions containing 1088 DMEs and 18 882 drug metabolites. Additionally, the INTEDE 2.0 labeled the pharmacological properties (pharmacological activity or toxicity) of metabolites and provided their structural information. Furthermore, 3717 drug metabolism relationships were supplemented (from 7338 to 11 055). All in all, INTEDE 2.0 is highly expected to attract broad interests from related research community and serve as an essential supplement to existing pharmaceutical/biological/chemical databases. INTEDE 2.0 can now be accessible freely without any login requirement at: http://idrblab.org/intede/.

Enzyme Databases in the Era of Omics and Artificial Intelligence

Enzyme research is important for the development of various scientific fields such as medicine and biotechnology. Enzyme databases facilitate this research by providing a wide range of information relevant to research planning and data analysis. Over the years, various databases that cover different aspects of enzyme biology (e.g., kinetic parameters, enzyme occurrence, and reaction mechanisms) have been developed. Most of the databases are curated manually, which improves reliability of the information; however, such curation cannot keep pace with the exponential growth in published data. Lack of data standardization is another obstacle for data extraction and analysis. Improving machine readability of databases is especially important in the light of recent advances in deep learning algorithms that require big training datasets. This review provides information regarding the current state of enzyme databases, especially in relation to the ever-increasing amount of generated research data and recent advancements in artificial intelligence algorithms. Furthermore, it describes several enzyme databases, providing the reader with necessary information for their use.

The Important Role of Transporter Structures in Drug Disposition, Efficacy, and Toxicity

The ATP-binding cassette (ABC) and solute carrier (SLC) transporters are critical determinants of drug disposition, clinical efficacy, and toxicity as they specifically mediate the influx and efflux of various substrates and drugs. ABC transporters can modulate the pharmacokinetics of many drugs via mediating the translocation of drugs across biologic membranes. SLC transporters are important drug targets involved in the uptake of a broad range of compounds across the membrane. However, high-resolution experimental structures have been reported for a very limited number of transporters, which limits the study of their physiologic functions. In this review, we collected structural information on ABC and SLC transporters and described the application of computational methods in structure prediction. Taking P-glycoprotein (ABCB1) and serotonin transporter (SLC6A4) as examples, we assessed the pivotal role of structure in transport mechanisms, details of ligand-receptor interactions, drug selectivity, the molecular mechanisms of drug-drug interactions, and differences caused by genetic polymorphisms. The data collected contributes toward safer and more effective pharmacological treatments. SIGNIFICANCE STATEMENT: The experimental structure of ATP-binding cassette and solute carrier transporters was collected, and the application of computational methods in structure prediction was described. P-glycoprotein and serotonin transporter were used as examples to reveal the pivotal role of structure in transport mechanisms, drug selectivity, the molecular mechanisms of drug-drug interactions, and differences caused by genetic polymorphisms.

Rm-LR: A long-range-based deep learning model for predicting multiple types of RNA modifications

Recent research has highlighted the pivotal role of RNA post-transcriptional modifications in the regulation of RNA expression and function. Accurate identification of RNA modification sites is important for understanding RNA function. In this study, we propose a novel RNA modification prediction method, namely Rm-LR, which leverages a long-range-based deep learning approach to accurately predict multiple types of RNA modifications using RNA sequences only. Rm-LR incorporates two large-scale RNA language pre-trained models to capture discriminative sequential information and learn local important features, which are subsequently integrated through a bilinear attention network. Rm-LR supports a total of ten RNA modification types (m6A, m1A, m5C, m5U, m6Am, Ψ, Am, Cm, Gm, and Um) and significantly outperforms the state-of-the-art methods in terms of predictive capability on benchmark datasets. Experimental results show the effectiveness and superiority of Rm-LR in prediction of various RNA modifications, demonstrating the strong adaptability and robustness of our proposed model. We demonstrate that RNA language pretrained models enable to learn dense biological sequential representations from large-scale long-range RNA corpus, and meanwhile enhance the interpretability of the models. This work contributes to the development of accurate and reliable computational models for RNA modification prediction, providing insights into the complex landscape of RNA modifications.

DCABM-TCM: A Database of Constituents Absorbed into the Blood and Metabolites of Traditional Chinese Medicine

Traditional Chinese medicine (TCM) not only maintains the health of Asian people but also provides a great resource of active natural products for modern drug development. Herein, we developed a Database of Constituents Absorbed into the Blood and Metabolites of TCM (DCABM-TCM), the first database systematically collecting blood constituents of TCM prescriptions and herbs, including prototypes and metabolites experimentally detected in the blood, together with the corresponding detailed detection conditions through manual literature mining. The DCABM-TCM has collected 1816 blood constituents with chemical structures of 192 prescriptions and 194 herbs and integrated their related annotations, including physicochemical, absorption, distribution, metabolism, excretion, and toxicity properties, and associated targets, pathways, and diseases. Furthermore, the DCABM-TCM supported two blood constituent-based analysis functions, the network pharmacology analysis for TCM molecular mechanism elucidation, and the target/pathway/disease-based screening of candidate blood constituents, herbs, or prescriptions for TCM-based drug discovery. The DCABM-TCM is freely accessible at http://bionet.ncpsb.org.cn/dcabm-tcm/. The DCABM-TCM will contribute to the elucidation of effective constituents and molecular mechanism of TCMs and the discovery of TCM-derived drug-like compounds that are both bioactive and bioavailable.

Discovery of selective HDAC6 inhibitors based on a multi-layer virtual screening strategy

The abnormal enhancement of histone deacetylase 6 (HDAC6) has been demonstrated to be closely related to the occurrence and development of various malignant tumors, attracting extensive attention as a promising target for cancer therapy. Currently, only limited selective HDAC6 inhibitors have entered clinical trials, making the rapid discovery of selective HDAC6 inhibitors with safety profiles particularly urgent. In this study, a multi-layer virtual screening workflow was established, and the representative compounds screened were biologically evaluated in combination with enzyme inhibitory and anti-tumor cell proliferation experiments. The experimental results showed that the screened compounds L-25, L-32, L-45 and L-81 exhibited nanomolar inhibitory activity against HDAC6, and exerted a certain degree of anti-proliferative activities against tumor cells, especially the cytotoxicity of L-45 to A375 (IC₅₀ = 11.23 ± 1.27 μM) and the cytotoxicity of L-81 against HCT-116 (IC₅₀ = 12.25 ± 1.13 μM). Additionally, the molecular mechanisms underlying the subtype selective inhibitory activities of the selected compounds were further elucidated using computational approaches, and the hotspot residues on HDAC6 contributing to the ligands' binding were identified. In summary, this study established a multi-layer screening scheme to quickly and effectively screen out hit compounds with enzyme inhibitory activity and anti-tumor cell proliferation, providing novel scaffolds for the subsequent anti-tumor drug design based on HDAC6 target.

MecDDI: Clarified Drug-Drug Interaction Mechanism Facilitating Rational Drug Use and Potential Drug-Drug Interaction Prediction

Drug-drug interactions (DDIs) are a major concern in clinical practice and have been recognized as one of the key threats to public health. To address such a critical threat, many studies have been conducted to clarify the mechanism underlying each DDI, based on which alternative therapeutic strategies are successfully proposed. Moreover, artificial intelligence-based models for predicting DDIs, especially multilabel classification models, are highly dependent on a reliable DDI data set with clear mechanistic information. These successes highlight the imminent necessity to have a platform providing mechanistic clarifications for a large number of existing DDIs. However, no such platform is available yet. In this study, a platform entitled "MecDDI" was therefore introduced to systematically clarify the mechanisms underlying the existing DDIs. This platform is unique in (a) clarifying the mechanisms underlying over 1,78,000 DDIs by explicit descriptions and graphic illustrations and (b) providing a systematic classification for all collected DDIs based on the clarified mechanisms. Due to the long-lasting threats of DDIs to public health, MecDDI could offer medical scientists a clear clarification of DDI mechanisms, support healthcare professionals to identify alternative therapeutics, and prepare data for algorithm scientists to predict new DDIs. MecDDI is now expected as an indispensable complement to the available pharmaceutical platforms and is freely accessible at: https://idrblab.org/mecddi/.

Computational drug repurposing by exploiting large-scale gene expression data: Strategy, methods and applications

De novo drug development is an extremely complex, time-consuming and costly task. Urgent needs for therapies of various diseases have greatly accelerated searches for more effective drug development methods. Luckily, drug repurposing provides a new and effective perspective on disease treatment. Rapidly increased large-scale transcriptome data paints a detailed prospect of gene expression during disease onset and thus has received wide attention in the field of computational drug repurposing. However, how to efficiently mine transcriptome data and identify new indications for old drugs remains a critical challenge. This review discussed the irreplaceable role of transcriptome data in computational drug repurposing and summarized some representative databases, tools and strategies. More importantly, it proposed a practical guideline through establishing the correspondence between three gene expression data types and five strategies, which would facilitate researchers to adopt appropriate strategies to deeply mine large-scale transcriptome data and discover more effective therapies.

A systematic review of state-of-the-art strategies for machine learning-based protein function prediction

New drug discovery is inseparable from the discovery of drug targets, and the vast majority of the known targets are proteins. At the same time, proteins are essential structural and functional elements of living cells necessary for the maintenance of all forms of life. Therefore, protein functions have become the focus of many pharmacological and biological studies. Traditional experimental techniques are no longer adequate for rapidly growing annotation of protein sequences, and approaches to protein function prediction using computational methods have emerged and flourished. A significant trend has been to use machine learning to achieve this goal. In this review, approaches to protein function prediction based on the sequence, structure, protein-protein interaction (PPI) networks, and fusion of multi-information sources are discussed. The current status of research on protein function prediction using machine learning is considered, and existing challenges and prominent breakthroughs are discussed to provide ideas and methods for future studies.

Label-free proteome quantification and evaluation

The label-free quantification (LFQ) has emerged as an exceptional technique in proteomics owing to its broad proteome coverage, great dynamic ranges and enhanced analytical reproducibility. Due to the extreme difficulty lying in an in-depth quantification, the LFQ chains incorporating a variety of transformation, pretreatment and imputation methods are required and constructed. However, it remains challenging to determine the well-performing chain, owing to its strong dependence on the studied data and the diverse possibility of integrated chains. In this study, an R package EVALFQ was therefore constructed to enable a performance evaluation on >3000 LFQ chains. This package is unique in (a) automatically evaluating the performance using multiple criteria, (b) exploring the quantification accuracy based on spiking proteins and (c) discovering the well-performing chains by comprehensive assessment. All in all, because of its superiority in assessing from multiple perspectives and scanning among over 3000 chains, this package is expected to attract broad interests from the fields of proteomic quantification. The package is available at https://github.com/idrblab/EVALFQ.

In vitro and in silico perspectives to explain anticancer activity of a novel syringic acid analog ((4-(1H-1, 3-benzodiazol-2-yl)-2, 6-dimethoxy phenol)) through apoptosis activation and NFkB inhibition in K562 leukemia cells

Syringic acid (SA) is an active carcinogenesis inhibitor; however, the low bioavailability and unstable functional groups hinder its activity. Here, a chemically synthesized novel SA analog (SA10) is evaluated for its anticancer activity using in-vitro and in-silico studies. K562 cell line study revealed that SA10 had shown a higher rate of inhibition (IC50 = 50.40 μg/mL) than its parental compound, SA (IC50 = 96.92 μg/mL), at 50 μM concentration. The inhibition ratio was also been evaluated by checking the expression level of NFkB and Bcl-2 and showing that SA10 has two-fold increase in the inhibitory mechanism than SA. This result demonstrates that SA10 acts as an NFkB inhibitor and an apoptosis inducer. Further, molecular docking and simulation have been performed to get insights into the possible inhibitory mechanism of SA and SA10 on NFkB at the atomistic level. The molecular docking results exemplify that both SA and SA10 bind to the active site of NFkB, thereby interfering with the association between DNA and NFkB. SA10 exhibits a more robust binding affinity than SA and is firmly docked well into the interior of the NFkB, as confirmed by MM-PBSA calculations. In a nutshell, the Benzimidazole scaffold containing SA10 has shown more NFkB inhibitory activity in K562 cells than SA, which could be helpful as an ideal therapeutic NFkB inhibitor for treating cancers.

A systematic review on the state-of-the-art strategies for protein representation

The study of drug-target protein interaction is a key step in drug research. In recent years, machine learning techniques have become attractive for research, including drug research, due to their automated nature, predictive power, and expected efficiency. Protein representation is a key step in the study of drug-target protein interaction by machine learning, which plays a fundamental role in the ultimate accomplishment of accurate research. With the progress of machine learning, protein representation methods have gradually attracted attention and have consequently developed rapidly. Therefore, in this review, we systematically classify current protein representation methods, comprehensively review them, and discuss the latest advances of interest. According to the information extraction methods and information sources, these representation methods are generally divided into structure and sequence-based representation methods. Each primary class can be further divided into specific subcategories. As for the particular representation methods involve both traditional and the latest approaches. This review contains a comprehensive assessment of the various methods which researchers can use as a reference for their specific protein-related research requirements, including drug research.

Transcription Factors and ncRNAs Associated with CYP3A Expression in Human Liver and Small Intestine Assessed with Weighted Gene Co-Expression Network Analysis

CYP3A4, CYP3A5, and CYP3A7, which are located in a multigene locus (CYP3A), play crucial roles in drug metabolism. To understand the highly variable hepatic expression of CYP3As, regulatory network analyses have focused on transcription factors (TFs). Since long non-coding RNAs (lncRNAs) likely contribute to such networks, we assessed the regulatory effects of both TFs and lncRNAs on CYP3A expression in the human liver and small intestine, main organs of CYP3A expression. Using weighted gene co-expression network analysis (WGCNA) of GTEx v8 RNA expression data and multiple stepwise regression analysis, we constructed TF-lncRNA-CYP3A co-expression networks. Multiple lncRNAs and TFs displayed robust associations with CYP3A expression that differed between liver and small intestines (LINC02499, HNF4A-AS1, AC027682.6, LOC102724153, and RP11-503C24.6), indicating that lncRNAs contribute to variance in CYP3A expression in both organs. Of these, HNF4A-AS1 had been experimentally demonstrated to affect CYP3A expression. Incorporating ncRNAs into CYP3A expression regulatory network revealed additional candidate TFs associated with CYP3A expression. These results serve as a guide for experimental studies on lncRNA-TF regulation of CYP3A expression in the liver and small intestines.

LargeMetabo: an out-of-the-box tool for processing and analyzing large-scale metabolomic data

Large-scale metabolomics is a powerful technique that has attracted widespread attention in biomedical studies focused on identifying biomarkers and interpreting the mechanisms of complex diseases. Despite a rapid increase in the number of large-scale metabolomic studies, the analysis of metabolomic data remains a key challenge. Specifically, diverse unwanted variations and batch effects in processing many samples have a substantial impact on identifying true biological markers, and it is a daunting challenge to annotate a plethora of peaks as metabolites in untargeted mass spectrometry-based metabolomics. Therefore, the development of an out-of-the-box tool is urgently needed to realize data integration and to accurately annotate metabolites with enhanced functions. In this study, the LargeMetabo package based on R code was developed for processing and analyzing large-scale metabolomic data. This package is unique because it is capable of (1) integrating multiple analytical experiments to effectively boost the power of statistical analysis; (2) selecting the appropriate biomarker identification method by intelligent assessment for large-scale metabolic data and (3) providing metabolite annotation and enrichment analysis based on an enhanced metabolite database. The LargeMetabo package can facilitate flexibility and reproducibility in large-scale metabolomics. The package is freely available from https://github.com/LargeMetabo/LargeMetabo.

DRESIS: the first comprehensive landscape of drug resistance information

Widespread drug resistance has become the key issue in global healthcare. Extensive efforts have been made to reveal not only diverse diseases experiencing drug resistance, but also the six distinct types of molecular mechanisms underlying this resistance. A database that describes a comprehensive list of diseases with drug resistance (not just cancers/infections) and all types of resistance mechanisms is now urgently needed. However, no such database has been available to date. In this study, a comprehensive database describing drug resistance information named 'DRESIS' was therefore developed. It was introduced to (i) systematically provide, for the first time, all existing types of molecular mechanisms underlying drug resistance, (ii) extensively cover the widest range of diseases among all existing databases and (iii) explicitly describe the clinically/experimentally verified resistance data for the largest number of drugs. Since drug resistance has become an ever-increasing clinical issue, DRESIS is expected to have great implications for future new drug discovery and clinical treatment optimization. It is now publicly accessible without any login requirement at: https://idrblab.org/dresis/.

DrugMAP: molecular atlas and pharma-information of all drugs

The efficacy and safety of drugs are widely known to be determined by their interactions with multiple molecules of pharmacological importance, and it is therefore essential to systematically depict the molecular atlas and pharma-information of studied drugs. However, our understanding of such information is neither comprehensive nor precise, which necessitates the construction of a new database providing a network containing a large number of drugs and their interacting molecules. Here, a new database describing the molecular atlas and pharma-information of drugs (DrugMAP) was therefore constructed. It provides a comprehensive list of interacting molecules for >30 000 drugs/drug candidates, gives the differential expression patterns for >5000 interacting molecules among different disease sites, ADME (absorption, distribution, metabolism and excretion)-relevant organs and physiological tissues, and weaves a comprehensive and precise network containing >200 000 interactions among drugs and molecules. With the great efforts made to clarify the complex mechanism underlying drug pharmacokinetics and pharmacodynamics and rapidly emerging interests in artificial intelligence (AI)-based network analyses, DrugMAP is expected to become an indispensable supplement to existing databases to facilitate drug discovery. It is now fully and freely accessible at: https://idrblab.org/drugmap/.

Molecular dynamics simulations highlight the altered binding landscape at the spike-ACE2 interface between the Delta and Omicron variants compared to the SARS-CoV-2 original strain

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) B.1.1.529 variant (Omicron), represents a significant deviation in genetic makeup and function compared to previous variants. Following the BA.1 sublineage, the BA.2 and BA.3 Omicron subvariants became dominant, and currently the BA.4 and BA.5, which are quite distinct variants, have emerged. Using molecular dynamics simulations, we investigated the binding characteristics of the Delta and Omicron (BA.1) variants in comparison to wild-type (WT) at the interface of the spike protein receptor binding domain (RBD) and human angiotensin converting enzyme-2 (ACE2) ectodomain. The primary aim was to compare our molecular modelling systems with previously published observations, to determine the robustness of our approach for rapid prediction of emerging future variants. Delta and Omicron were found to bind to ACE2 with similar affinities (−39.4 and −43.3 kcal/mol, respectively) and stronger than WT (−33.5 kcal/mol). In line with previously published observations, the energy contributions of the non-mutated residues at the interface were largely retained between WT and the variants, with F456, F486, and Y489 having the strongest energy contributions to ACE2 binding. Further, residues N440K, Q498R, and N501Y were predicted to be energetically favourable in Omicron. In contrast to Omicron, which had the E484A and K417N mutations, intermolecular bonds were detected for the residue pairs E484:K31 and K417:D30 in WT and Delta, in accordance with previously published findings. Overall, our simplified molecular modelling approach represents a step towards predictive model systems for rapidly analysing arising variants of concern.

M6AREG: m6A-centered regulation of disease development and drug response

As the most prevalent internal modification in eukaryotic RNAs, N6-methyladenosine (m6A) has been discovered to play an essential role in cellular proliferation, metabolic homeostasis, embryonic development, etc. With the rapid accumulation of research interest in m6A, its crucial roles in the regulations of disease development and drug response are gaining more and more attention. Thus, a database offering such valuable data on m6A-centered regulation is greatly needed; however, no such database is as yet available. Herein, a new database named 'M6AREG' is developed to (i) systematically cover, for the first time, data on the effects of m6A-centered regulation on both disease development and drug response, (ii) explicitly describe the molecular mechanism underlying each type of regulation and (iii) fully reference the collected data by cross-linking to existing databases. Since the accumulated data are valuable for researchers in diverse disciplines (such as pathology and pathophysiology, clinical laboratory diagnostics, medicinal biochemistry and drug design), M6AREG is expected to have many implications for the future conduct of m6A-based regulation studies. It is currently accessible by all users at: https://idrblab.org/m6areg/.

DeepPSE: Prediction of polypharmacy side effects by fusing deep representation of drug pairs and attention mechanism

Polypharmacy (multiple use of drugs) is an effective strategy for combating complex or co-existing diseases. However, a major consequence of polypharmacy is a higher risk of adverse side effects due to drug-drug interactions, which are rare and observed in relatively small clinical testing. Thus, identification of polypharmacy side effects remains challenging. Here, we propose a deep learning-based method, DeepPSE, to predict polypharmacy side effects in an end-to-end way. DeepPSE is composed of two main modules. First, multiple types of neural networks are constructed and fused to learn the deep representation of a drug pair. Second, the encoder block of transformer that includes self-attention mechanism is built to get latent features, which are further fed into the fully connected layer to predict polypharmacy side effects of drug pairs. Further, DeepPSE is compared with five baseline or state-of-the-art methods on a benchmark dataset of 964 types of polypharmacy side effects across 63473 drug pairs. Experimental results demonstrate that DeepPSE achieves better performance than that of all five methods. The source codes and data are available at https://github.com/ShenggengLin/DeepPSE.

Integrate thermostabilized fusion protein apocytochrome b562RIL and N-glycosylation mutations: A novel approach to heterologous expression of human UDP-glucuronosyltransferase (UGT) 2B7

Human UDP-glucuronosyltransferase (UGT) 2B7 is a crucial phase II metabolic enzyme that transfers glucuronic acid from UDP-glucuronic acid (UDPGA) to endobiotic and xenobiotic substrates. Biophysical and biochemical investigations of UGT2B7 are hampered by the challenge of the integral membrane protein purification. This study focused on the expression and purification of recombinant UGT2B7 by optimizing the insertion sites for the thermostabilized fusion protein apocytochrome b₅₆₂RIL (BRIL) and various mutations to improve the protein yields and homogeneity. Preparation of the recombinant proteins with high purity accelerated the measurement of pharmacokinetic parameters of UGT2B7. The dissociation constants (K_D) of two classical substrates (zidovudine and androsterone) and two inhibitors (schisanhenol and hesperetin) of UGT2B7 were determined using the surface plasmon resonance spectroscopy (SPR) for the first time. Using negative-staining transmission electron microscopy (TEM), UGT2B7 protein particles were characterized, which could be useful for further exploring its three-dimensional structure. The methods described in this study could be broadly applied to other UGTs and are expected to provide the basis for the exploration of metabolic enzyme kinetics, the mechanisms of drug metabolisms and drug interactions, changes in pharmacokinetics, and pharmacodynamics studies in vitro.

In-silico screening and in-vitro assay show the antiviral effect of Indomethacin against SARS-CoV-2

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused the worldwide spread of coronavirus disease 19 (COVID-19), and till now, it has caused death to more than 6.2 million people. Although various vaccines and drug candidates are being tested globally with limited to moderate success, a comprehensive therapeutic cure is yet to be achieved. In this study, we applied computational drug repurposing methods complemented with the analyses of the already existing gene expression data to find better therapeutics in treatment and recovery. Primarily, we identified the most crucial proteins of SARS-CoV-2 and host human cells responsible for viral infection and host response. An in-silico screening of the existing drugs was performed against the crucial proteins for SARS-CoV-2 infection, and a few existing drugs were shortlisted. Further, we analyzed the gene expression data of SARS-CoV-2 in human lung epithelial cells and investigated the molecules that can reverse the cellular mRNA expression profiles in the diseased state. LINCS L1000 and Comparative Toxicogenomics Database (CTD) were utilized to obtain two sets of compounds that can be used to counter SARS-CoV-2 infection from the gene expression perspective. Indomethacin, a nonsteroidal anti-inflammatory drug (NSAID), and Vitamin-A were found in two sets of compounds, and in the in-silico screening of existing drugs to treat SARS-CoV-2. Our in-silico findings on Indomethacin were further successfully validated by in-vitro testing in Vero CCL-81 cells with an IC₅₀ of 12 μM. Along with these findings, we briefly discuss the possible roles of Indomethacin and Vitamin-A to counter the SARS-CoV-2 infection in humans.

REGLIV: Molecular regulation data of diverse living systems facilitating current multiomics research

Multiomics is a powerful technique in molecular biology that facilitates the identification of new associations among different molecules (genes, proteins & metabolites). It has attracted tremendous research interest from the scientists worldwide and has led to an explosive number of published studies. Most of these studies are based on the regulation data provided in available databases. Therefore, it is essential to have molecular regulation data that are strictly validated in the living systems of various cell lines and in vivo models. However, no database has been developed yet to provide comprehensive molecular regulation information validated by living systems. Herein, a new database, Molecular Regulation Data of Living System Facilitating Multiomics Study (REGLIV) is introduced to describe various types of molecular regulation tested by the living systems. (1) A total of 2996 regulations describe the changes in 1109 metabolites triggered by alterations in 284 genes or proteins, and (2) 1179 regulations describe the variations in 926 proteins induced by 125 endogenous metabolites. Overall, REGLIV is unique in (a) providing the molecular regulation of a clearly defined regulatory direction other than simple correlation, (b) focusing on molecular regulations that are validated in a living system not simply in an in vitro test, and (c) describing the disease/tissue/species specific property underlying each regulation. Therefore, REGLIV has important implications for the future practice of not only multiomics, but also other fields relevant to molecular regulation. REGLIV is freely accessible at: https://idrblab.org/regliv/.

Exploration and validation of Taraxacum mongolicum anti-cancer effect

Taraxacum mongolicum gained a lot of concern and was applied in 93 formulas in China due to its fame as a traditional Chinese medicine. The earliest recorded application of Taraxacum mongolicum was traced back to the Han dynasty. Generations of doctors boosted the usage and enriched the pharmacological mechanism. Clinical application of the Taraxacum mongolicum is flourishing as it treats multiple diseases. This study aims to explore the anti-cancer effect, retrieve the active ingredients and screen the key targets of Taraxacum mongolicum in cancer therapy. We collected and evaluated 10 key active compounds to investigate the anti-cancer effect via 69 significant targets and a variety of biological processes and pathways. Gene Ontology (GO) enrichment analysis uncovered targets associated with protein phosphorylation, cell proliferation and apoptotic processes via regulation of kinases, ATP and enzyme binding activities. Half of the top 20 enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were directly involved in cancer. Based on standard selection criteria, seven hub targets were obtained. These targets functioned through distinct patterns and pathways in realizing the anti-cancer effect. Molecular docking was conducted to validate the potential combination between compounds and hub targets to explore the pharmacological mechanism of key compounds in Taraxacum mongolicum against cancer. In summary, our findings indicate that the famous and widely used Chinese herb, Taraxacum mongolicum, shows good anti-cancer effect through its active compounds, targeted genes, and multiple involved biological processes. The results may provide a theoretical basis for subsequent experimental validation and drug development of Taraxacum mongolicum extract against cancer.

ConSIG: consistent discovery of molecular signature from OMIC data

The discovery of proper molecular signature from OMIC data is indispensable for determining biological state, physiological condition, disease etiology, and therapeutic response. However, the identified signature is reported to be highly inconsistent, and there is little overlap among the signatures identified from different biological datasets. Such inconsistency raises doubts about the reliability of reported signatures and significantly hampers its biological and clinical applications. Herein, an online tool, ConSIG, was constructed to realize consistent discovery of gene/protein signature from any uploaded transcriptomic/proteomic data. This tool is unique in a) integrating a novel strategy capable of significantly enhancing the consistency of signature discovery, b) determining the optimal signature by collective assessment, and c) confirming the biological relevance by enriching the disease/gene ontology. With the increasingly accumulated concerns about signature consistency and biological relevance, this online tool is expected to be used as an essential complement to other existing tools for OMIC-based signature discovery. ConSIG is freely accessible to all users without login requirement at https://idrblab.org/consig/.

Computational approach to decode the mechanism of curcuminoids against neuropathic pain

BACKGROUND

Curcumin (CUR), demethoxycurcumin (DMC) and bisdemethoxycurcumin (BDMC) are the main components of turmeric that commonly used to treat neuropathic pain (NP). However, the mechanism of the therapy is not sufficiently clarified. Herein, network pharmacology, molecular docking and molecular dynamics (MD) approaches were used to investigate the mechanism of curcuminoids for NP treatment.

METHODS

Active targets of curcuminoids were obtained from the Swiss Target database, and NP-related targets were retrieved from GeneCards, OMIM, Drugbank and TTD databases. A protein-protein interaction (PPI) network was built to screen the core targets. Furthermore, DAVID was used for GO and KEGG pathway enrichment analyses. Interactions between potential targets and curcuminoids were assessed by molecular docking and the MD simulations were run for 100ns to validate the docking results on the top six complexes.

RESULTS

CUR, DMC, and BDMC had 100, 99 and 100 targets respectively. After overlapping with NP there were 33, 33 and 31 targets respectively. PPI network analysis of TOP 10 core targets, TNF, GSK3β were common targets of curcuminoids. Molecular docking and MD results indicated that curcuminoids bind strongly with the core targets. The GO and KEGG showed that curcuminoids regulated nitrogen metabolism, the serotonergic synapse and ErbB signaling pathway to alleviate NP. Furthermore, specific targets in these three compounds were also analysed at the same time.

CONCLUSIONS

This study systematically explored and compared the anti-NP mechanism of curcuminoids, providing a novel perspective for their utilization.

A multidrug efflux protein in Mycobacterium tuberculosis; tap as a potential drug target for drug repurposing

Tuberculosis (TB) is a serious communicative disease caused by Mycobacterium tuberculosis. Although there are vaccines and drugs available to treat the disease, they are not efficient, moreover, multidrug-resistant TB (MDR-TB) become a major hurdle in its therapy. These MDR strains utilize the multidrug efflux pump as a decisive weapon to fight against antitubercular drugs. Tap membrane protein was observed as a crucial multidrug efflux pump in M. tuberculosis and its critical implication in MDR-MTB development makes it an effective drug target. In the present study, we have utilized various in silico approaches to predict the applicability of FDA-approved ion channel inhibitors and blockers as therapeutic leads against Tuberculosis by targeting multidrug efflux protein; Tap in MTB. Tap protein structure is predicted by Phyre2 server followed by model refinement, validation, physio-chemical catheterization and target prediction. Further, the interaction between Tap protein and ligands were analysed by molecular docking and MD simulation run of 100 ns. Based on implication and compatibility, 18 FDA-approved ion channel inhibitors and blockers are selected as a ligand against the Tap protein and eventually observed five ligands; Glimepiride, Flecainide, Flupiritine, Nimodipine and Amlodipine as promising compounds which have exhibited the significant stable interaction with Tap protein and are proposed to modulate or interfere with its activity. These compounds illustrated the substantial docking score and total binding enthalpy more than -7 kcal/mol and -42 kcal/mol respectively which implies that the selected FDA-approved compounds can spontaneously interact with the Tap protein to modulate its function. This study proposed Tap protein as a prominent drug target in MTB and investigated compounds that show considerable interaction with the Tap protein as potential therapeutic molecules. These interactions may lead to modulating or inhibit the activity of drug efflux protein thereby making MTB susceptible to antitubercular drugs.

Biological activities of drug inactive ingredients

In a drug formulation (DFM), the major components by mass are not Active Pharmaceutical Ingredient (API) but rather Drug Inactive Ingredients (DIGs). DIGs can reach much higher concentrations than that achieved by API, which raises great concerns about their clinical toxicities. Therefore, the biological activities of DIG on physiologically relevant target are widely demanded by both clinical investigation and pharmaceutical industry. However, such activity data are not available in any existing pharmaceutical knowledge base, and their potentials in predicting the DIG-target interaction have not been evaluated yet. In this study, the comprehensive assessment and analysis on the biological activities of DIGs were therefore conducted. First, the largest number of DIGs and DFMs were systematically curated and confirmed based on all drugs approved by US Food and Drug Administration. Second, comprehensive activities for both DIGs and DFMs were provided for the first time to pharmaceutical community. Third, the biological targets of each DIG and formulation were fully referenced to available databases that described their pharmaceutical/biological characteristics. Finally, a variety of popular artificial intelligence techniques were used to assess the predictive potential of DIGs' activity data, which was the first evaluation on the possibility to predict DIG's activity. As the activities of DIGs are critical for current pharmaceutical studies, this work is expected to have significant implications for the future practice of drug discovery and precision medicine.

Mosquitocidal efficacy of embelin and its derivatives against Aedes aegypti L. and Culex quinquefasciatus Say. (Diptera: Culicidae) and computational analysis of acetylcholinesterase 1 (AChE1) inhibition

Embelin was isolated from the chloroform extract of Embelia ribes (Burm.f.) fruits; its derivative compounds 6-bromoembelin and vilangin were prepared, and they were evaluated for mosquitocidal activities against the third instar larvae and pupae of Aedes aegypti L. and Culex quinquefasciatus Say. (Diptera: Culicidae). The concentrations used were 0.5, 1.0, 1.5, and 2.0 ppm. Embelin recorded LC₅₀ values of 5.79 and 5.54 ppm against the larvae of Ae. aegypti and Cx. quinquefasciatus, respectively. Similarly, the LC₅₀ values of embelin were 10.23 and 6.93 ppm against the pupae of Ae. aegypti and Cx. quinquefasciatus, respectively. Of the two derivatives tested, vilangin showed the highest larvicidal activity with LC₅₀ values of 1.38 and 1.28 ppm against the larvae of Ae. aegypti and Cx. quinquefasciatus, respectively. Similarly, the LC₅₀ values of vilangin were 1.60 and 1.43 ppm against the pupae of Ae. aegypti and Cx. quinquefasciatus, respectively. The LC₅₀ values of 6-bromoembelin were 3.30 and 2.83 ppm against the larvae and 4.40 and 4.30 ppm against the pupae of Ae. aegypti and Cx. quinquefasciatus, respectively. The histopathological results displayed significant damage on cuboidal cells of the midgut (CU) in vilangin treated larvae of Ae. aegypti and Cx. quinquefasciatus at a concentration of 2.0 ppm. Similarly, peritrophic membrane (PM) was completely impaired in vilangin-treated larvae of Cx. quinquefasciatus and midgut content (MC) was very low in vilangin-treated larvae of Cx. quinquefasciatus. In addition, molecular docking and molecular dynamics studies demonstrated the efficacy of vilangin on the inhibition of acetylcholinesterase (AChE1) in Ae. aegypti and Cx. quinquefasciatus. The present results suggest that vilangin could be used to develop a natural active product against mosquito larvae.

PregTox: A Resource of Knowledge about Drug Fetal Toxicity

Background

It is of vital importance to determine the safety of drugs. Pregnant women, as a special group, need to evaluate the effects of drugs on pregnant women as well as the fetus. The use of drugs during pregnancy may be subject to fetal toxicity, thus affecting the development of the fetus or even leading to stillbirth. The U.S. Food and Drug Administration (FDA) issued a toxicity rating for drugs used during pregnancy in 1979. These toxicity ratings are denoted by the letters A, B, C, D, and X. However, the query of drug pregnancy category has yet to be well established as electronic service.

Results

Here, we presented PregTox, a publicly accessible resource for pregnancy category information of 1114 drugs. The PregTox database also included chemical structures, important physico-chemical properties, protein targets, and relevant signaling pathways. An advantage of the database is multiple search options which allow systematic analyses. In a case study, we demonstrated that a set of chemical descriptors could effectively discriminate high-risk drugs from others (area under ROC curve reached 0.81).

Conclusions

PregTox can serve as a unique drug safety data source for drug development and pharmacological research.

RGCNCDA: Relational graph convolutional network improves circRNA-disease association prediction by incorporating microRNAs

Recently, a large number of studies have indicated that circRNAs with covalently closed loops play important roles in biological processes and have potential as diagnostic biomarkers. Therefore, research on the circRNA-disease relationship is helpful in disease diagnosis and treatment. However, traditional biological verification methods require considerable labor and time costs. In this paper, we propose a new computational method (RGCNCDA) to predict circRNA-disease associations based on relational graph convolutional networks (R-GCNs). The method first integrates the circRNA similarity network, miRNA similarity network, disease similarity network and association networks among them to construct a global heterogeneous network. Then, it employs the random walk with restart (RWR) and principal component analysis (PCA) models to learn low-dimensional and high-order information from the global heterogeneous network as the topological features. Finally, a prediction model based on an R-GCN encoder and a DistMult decoder is built to predict the potential disease-associated circRNA. The predicted results demonstrate that RGCNCDA performs significantly better than the other six state-of-the-art methods in a 5-fold cross validation. Furthermore, the case study illustrates that RGCNCDA can effectively discover potential circRNA-disease associations.

PFmulDL: a novel strategy enabling multi-class and multi-label protein function annotation by integrating diverse deep learning methods

Bioinformatic annotation of protein function is essential but extremely sophisticated, which asks for extensive efforts to develop effective prediction method. However, the existing methods tend to amplify the representativeness of the families with large number of proteins by misclassifying the proteins in the families with small number of proteins. That is to say, the ability of the existing methods to annotate proteins in the 'rare classes' remains limited. Herein, a new protein function annotation strategy, PFmulDL, integrating multiple deep learning methods, was thus constructed. First, the recurrent neural network was integrated, for the first time, with the convolutional neural network to facilitate the function annotation. Second, a transfer learning method was introduced to the model construction for further improving the prediction performances. Third, based on the latest data of Gene Ontology, the newly constructed model could annotate the largest number of protein families comparing with the existing methods. Finally, this newly constructed model was found capable of significantly elevating the prediction performance for the 'rare classes' without sacrificing that for the 'major classes'. All in all, due to the emerging requirements on improving the prediction performance for the proteins in 'rare classes', this new strategy would become an essential complement to the existing methods for protein function prediction. All the models and source codes are freely available and open to all users at: https://github.com/idrblab/PFmulDL.

A novel mTOR-associated gene signature for predicting prognosis and evaluating tumor immune microenvironment in lung adenocarcinoma

BACKGROUND

The mechanistic target of rapamycin (mTOR) was proven to have great impact on apoptosis, cell proliferation, autophagy, and many other fundamental cellular processes; moreover, it closely correlates with tumor occurrence and development. However, few studies have constructed signatures based on mTOR-associated genes to assess multiple indicators of prognosis in lung adenocarcinoma (LUAD) patients.

METHODS

mTOR-associated gene sets, whole mRNA expression matrices, and clinical information of LUAD patients in training and validation cohorts were obtained from multiple public databases. Multiple methods were used to screen candidate genes, construct signatures, validate internally and externally, and conduct further studies: differentially expressed gene analysis, LASSO Cox regression analysis, Cox regression analysis, risk factor analysis, nomogram analysis, functional enrichment analysis, analyses in tumor immune microenvironment, and therapy.

RESULTS

A prognostic signature containing 8 genes (LDHA, SLA, WNT7A, PLK1, CCT6A, BTG2, TXNRD1, and DDIT4) was constructed. It performed well in both internal and external validation. Subsequent analysis found that the prognostic signature was of great significance in evaluating the tumor immune microenvironment and could guide the treatment of patients with LUAD to a certain extent.

CONCLUSION

The constructed mTOR-associated gene signature accurately predicted the prognostic pattern of patients with LUAD and is expected to be extremely useful in guiding LUAD therapy.

HKAM-MKM: A hybrid kernel alignment maximization-based multiple kernel model for identifying DNA-binding proteins

The identification of DNA-binding proteins (DBPs) has always been a hot issue in the field of sequence classification. However, considering that the experimental identification method is very resource-intensive, the construction of a computational prediction model is worthwhile. This study developed and evaluated a hybrid kernel alignment maximization-based multiple kernel model (HKAM-MKM) for predicting DBPs. First, we collected two datasets and performed feature extraction on the sequences to obtain six feature groups, and then constructed the corresponding kernels. To ensure the effective utilisation of the base kernel and avoid ignoring the difference between the sample and its neighbours, we proposed local kernel alignment to calculate the kernel between the sample and its neighbours, with each sample as the centre. We combined the global and local kernel alignments to develop a hybrid kernel alignment model, and balance the relationship between the two through parameters. By maximising the hybrid kernel alignment value, we obtained the weight of each kernel and then linearly combined the kernels in the form of weights. Finally, the fused kernel was input into a support vector machine for training and prediction. Finally, in the independent test sets PDB186 and PDB2272, we obtained the highest Matthew's correlation coefficient (MCC) (0.768 and 0.5962, respectively) and the highest accuracy (87.1% and 78.43%, respectively), which were superior to the other predictors. Therefore, HKAM-MKM is an efficient prediction tool for DBPs.

TMPpred: A support vector machine-based thermophilic protein identifier

MOTIVATION

The thermostability of proteins will cause them to break the temperature binding and play more functions. Using machine learning, we explored the mechanism of and reasons for protein thermostability characteristics.

RESULTS

Different from other methods that only pursue the performance of models, we aim to find important features so as to provide a powerful reference for in vitro experiments. We transformed this problem into a binary classification problem, that is, the distinction between thermophilic proteins and nonthermophilic proteins. Using support vector machine-based model construction and analysis, we inferred that Gly, Ala, Ser and Thr may be the most important components at the residue level that determine the thermal stability of proteins. It is also noteworthy that our proposed model obtains an Sn of 0.892, an Sp of 0.857, an ACC of 0.87566 and an AUC of 0.874. To facilitate other researchers, we wrapped our model and deployed it as a web server, which is accessible at http://112.124.26.17:7000/TMPpred/index.html.

POSREG: proteomic signature discovered by simultaneously optimizing its reproducibility and generalizability

Mass spectrometry-based proteomic technique has become indispensable in current exploration of complex and dynamic biological processes. Instrument development has largely ensured the effective production of proteomic data, which necessitates commensurate advances in statistical framework to discover the optimal proteomic signature. Current framework mainly emphasizes the generalizability of the identified signature in predicting the independent data but neglects the reproducibility among signatures identified from independently repeated trials on different sub-dataset. These problems seriously restricted the wide application of the proteomic technique in molecular biology and other related directions. Thus, it is crucial to enable the generalizable and reproducible discovery of the proteomic signature with the subsequent indication of phenotype association. However, no such tool has been developed and available yet. Herein, an online tool, POSREG, was therefore constructed to identify the optimal signature for a set of proteomic data. It works by (i) identifying the proteomic signature of good reproducibility and aggregating them to ensemble feature ranking by ensemble learning, (ii) assessing the generalizability of ensemble feature ranking to acquire the optimal signature and (iii) indicating the phenotype association of discovered signature. POSREG is unique in its capacity of discovering the proteomic signature by simultaneously optimizing its reproducibility and generalizability. It is now accessible free of charge without any registration or login requirement at https://idrblab.org/posreg/.

Biomedical data, computational methods and tools for evaluating disease-disease associations

In recent decades, exploring potential relationships between diseases has been an active research field. With the rapid accumulation of disease-related biomedical data, a lot of computational methods and tools/platforms have been developed to reveal intrinsic relationship between diseases, which can provide useful insights to the study of complex diseases, e.g. understanding molecular mechanisms of diseases and discovering new treatment of diseases. Human complex diseases involve both external phenotypic abnormalities and complex internal molecular mechanisms in organisms. Computational methods with different types of biomedical data from phenotype to genotype can evaluate disease-disease associations at different levels, providing a comprehensive perspective for understanding diseases. In this review, available biomedical data and databases for evaluating disease-disease associations are first summarized. Then, existing computational methods for disease-disease associations are reviewed and classified into five groups in terms of the usages of biomedical data, including disease semantic-based, phenotype-based, function-based, representation learning-based and text mining-based methods. Further, we summarize software tools/platforms for computation and analysis of disease-disease associations. Finally, we give a discussion and summary on the research of disease-disease associations. This review provides a systematic overview for current disease association research, which could promote the development and applications of computational methods and tools/platforms for disease-disease associations.

Molecular Mechanism for the Allosteric Inhibition of the Human Serotonin Transporter by Antidepressant Escitalopram

Human serotine transporter (hSERT) is one of the most influential drug targets, and its allosteric modulators (e.g., escitalopram) have emerged to be the next-generation medication for psychiatric disorders. However, the molecular mechanism underlying the allosteric modulation of hSERT is still elusive. Here, the simulation strategies of conventional (cMD) and steered (SMD) molecular dynamics were applied to investigate this molecular mechanism from distinct perspectives. First, cMD simulations revealed that escitalopram's binding to hSERT's allosteric site simultaneously enhanced its binding to the orthosteric site. Then, SMD simulation identified that the occupation of hSERT's allosteric site by escitalopram could also block its dissociation from the orthosteric site. Finally, by comparing the simulated structures of two hSERT-escitalopram complexes with and without allosteric modulation, a new conformational coupling between an extracellular (Arg104-Glu494) and an intracellular (Lys490-Glu494) salt bridge was identified. In summary, this study explored the mechanism underlying the allosteric modulation of hSERT by collectively applying two MD simulation strategies, which could facilitate our understanding of the allosteric modulations of not only hSERT but also other clinically important therapeutic targets.