Clinical trials, progression-speed differentiating features and swiftness rule of the innovative targets of first-in-class drugs

An integrated in silico approach for the identification of novel potential drug target and chimeric vaccine against Neisseria meningitides strain 331401 serogroup X by subtractive genomics and reverse vaccinology

Neisseria meningitidis, commonly known as the meningococcus, leads to substantial illness and death among children and young adults globally, revealing as either epidemic or sporadic meningitis and/or septicemia. In this study, we have designed a novel peptide-based chimeric vaccine candidate against the N. meningitidis strain 331,401 serogroup X. Through rigorous analysis of subtractive genomics, two essential cytoplasmic proteins, namely UPI000012E8E0(UDP-3-O-acyl-GlcNAc deacetylase) and UPI0000ECF4A9(UDP-N-acetylglucosamine acyltransferase) emerged as potential drug targets. Additionally, using reverse vaccinology, the outer membrane protein UPI0001F4D537 (Membrane fusion protein MtrC) identified by subcellular localization and recognized for its known indispensable role in bacterial survival was identified as a novel chimeric vaccine target. Following a careful comparison of MHC-I, MHC-II, T-cell, and B-cell epitopes, three epitopes derived from UPI0001F4D537 were linked with three types of linkers-GGGS, EAAAK, and the essential PADRE-for vaccine construction. This resulted in eight distinct vaccine models (V1-V8). Among them V1 model was selected as the final vaccine construct. It exhibits exceptional immunogenicity, safety, and enhanced antigenicity, with 97.7 % of its residues in the Ramachandran plot's most favored region. Subsequently, the vaccine structure was docked with the TLR4/MD2 complex and six different HLA allele receptors using the HADDOCK server. The docking resulted in the lowest HADDOCK score of 39.3 ± 9.0 for TLR/MD2. Immune stimulation showed a strong immune response, including antibodies creation and the activation of B-cells, T Cytotoxic cells, T Helper cells, Natural Killer cells, and interleukins. Furthermore, the vaccine construct was successfully expressed in the Escherichia coli system by reverse transcription, optimization, and ligation in the pET-28a (+) vector for the expression study. The current study proposes V1 construct has the potential to elicit both cellular and humoral responses, crucial for the developing an epitope-based vaccine against N. meningitidis strain 331,401 serogroup X.

Discovery of Potent Selective HDAC6 Inhibitors with 5-Phenyl-1H-indole Fragment: Virtual Screening, Rational Design, and Biological Evaluation

Among the HDACs family, histone deacetylase 6 (HDAC6) has attracted extensive attention due to its unique structure and biological functions. Numerous studies have shown that compared with broad-spectrum HDACs inhibitors, selective HDAC6 inhibitors exert ideal efficacy in tumor treatment with insignificant toxic and side effects, demonstrating promising clinical application prospect. Herein, we carried out rational drug design by integrating a deep learning model, molecular docking, and molecular dynamics simulation technology to construct a virtual screening process. The designed derivatives with 5-phenyl-1H-indole fragment as Cap showed desirable cytotoxicity to the various tumor cell lines, all of which were within 15 μM (ranging from 0.35 to 14.87 μM), among which compound 5i had the best antiproliferative activities against HL-60 (IC50 = 0.35 ± 0.07 μM) and arrested HL-60 cells in the G0/G1 phase. In addition, 5i exhibited better isotype selective inhibitory activities due to the potent potency against HDAC6 (IC50 = 5.16 ± 0.25 nM) and the reduced inhibitory activities against HDAC1 (selective index ≈ 124), which was further verified by immunoblotting results. Moreover, the representative binding conformation of 5i on HDAC6 was revealed and the key residues contributing 5i's binding were also identified via decomposition free-energy analysis. The discovery of lead compound 5i also indicates that virtual screening is still a beneficial tool in drug discovery and can provide more molecular skeletons with research potential for drug design, which is worthy of widespread application.

Identification of new targets for glioblastoma therapy based on a DNA expression microarray

This study provides a comprehensive perspective on the deregulated pathways and impaired biological functions prevalent in human glioblastoma (GBM). In order to characterize differences in gene expression between individuals diagnosed with GBM and healthy brain tissue, we have designed and manufactured a specific, custom DNA microarray. The results obtained from differential gene expression analysis were validated by RT-qPCR. The datasets obtained from the analysis of common differential expressed genes in our cohort of patients were used to generate protein-protein interaction networks of functionally enriched genes and their biological functions. This network analysis, let us to identify 16 genes that exhibited either up-regulation (CDK4, MYC, FOXM1, FN1, E2F7, HDAC1, TNC, LAMC1, EIF4EBP1 and ITGB3) or down-regulation (PRKACB, MEF2C, CAMK2B, MAPK3, MAP2K1 and PENK) in all GBM patients. Further investigation of these genes and enriched pathways uncovered in this investigation promises to serve as a foundational step in advancing our comprehension of the molecular mechanisms underpinning GBM pathogenesis. Consequently, the present work emphasizes the critical role that the unveiled molecular pathways likely play in shaping innovative therapeutic approaches for GBM management. We finally proposed in this study a list of compounds that target hub of GBM-related genes, some of which are already in clinical use, underscoring the potential of those genes as targets for GBM treatment.

Chaotic RIME optimization algorithm with adaptive mutualism for feature selection problems

The RIME optimization algorithm is a newly developed physics-based optimization algorithm used for solving optimization problems. The RIME algorithm proved high-performing in various fields and domains, providing a high-performance solution. Nevertheless, like many swarm-based optimization algorithms, RIME suffers from many limitations, including the exploration-exploitation balance not being well balanced. In addition, the likelihood of falling into local optimal solutions is high, and the convergence speed still needs some work. Hence, there is room for enhancement in the search mechanism so that various search agents can discover new solutions. The authors suggest an adaptive chaotic version of the RIME algorithm named ACRIME, which incorporates four main improvements, including an intelligent population initialization using chaotic maps, a novel adaptive modified Symbiotic Organism Search (SOS) mutualism phase, a novel mixed mutation strategy, and the utilization of restart strategy. The main goal of these improvements is to improve the variety of the population, achieve a better balance between exploration and exploitation, and improve RIME's local and global search abilities. The study assesses the effectiveness of ACRIME by using the standard benchmark functions of the CEC2005 and CEC2019 benchmarks. The proposed ACRIME is also applied as a feature selection to fourteen various datasets to test its applicability to real-world problems. Besides, the ACRIME algorithm is applied to the COVID-19 classification real problem to test its applicability and performance further. The suggested algorithm is compared to other sophisticated classical and advanced metaheuristics, and its performance is assessed using statistical tests such as Wilcoxon rank-sum and Friedman rank tests. The study demonstrates that ACRIME exhibits a high level of competitiveness and often outperforms competing algorithms. It discovers the optimal subset of features, enhancing the accuracy of classification and minimizing the number of features employed. This study primarily focuses on enhancing the equilibrium between exploration and exploitation, extending the scope of local search.

PT-Finder: A multi-modal neural network approach to target identification

Efficient target identification for bioactive compounds, including novel synthetic analogs, is crucial for accelerating the drug discovery pipeline. However, the process of target identification presents significant challenges and is often expensive, which in turn can hinder the drug discovery efforts. To address these challenges machine learning applications have arisen as a promising approach for predicting the targets for novel chemical compounds. These methods allow the exploration of ligand-target interactions, uncovering of biochemical mechanisms, and the investigation of drug repurposing. Typically, the current target identification tools rely on assessing ligand structural similarities. Herein, a multi-modal neural network model was built using a library of proteins, their respective sequences, and active inhibitors. Subsequent validations showed the model to possess accuracy of 82 % and MPRAUC of 0.80. Leveraging the trained model, we developed PT-Finder (Protein Target Finder), a user-friendly offline application that is capable of predicting the target proteins for hundreds of compounds within a few seconds. This combination of offline operation, speed, and accuracy positions PT-Finder as a powerful tool to accelerate drug discovery workflows. PT-Finder and its source codes have been made freely accessible for download at https://github.com/PT-Finder/PT-Finder.

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/.

Computer-aided investigation of Traditional Chinese Medicine mechanisms: A case study of San-Ao decoction in asthma treatment

The San-Ao Decoction (SAD) is a well-known Traditional Chinese Medicine (TCM) formula used to alleviate respiratory symptoms, including asthma. However, its precise mechanisms of action have remained largely unknown. In this study, we utilized computer-aided approaches to explore these mechanisms. Firstly, we conducted a comprehensive analysis of the chemical composition of SAD, which allowed us to identify the 28 main ingredients. Then, we employed computer simulations to investigate the potential active ingredients of SAD and the corresponding binding sites of transient receptor potential vanilloid 1 (TRPV1). The simulations revealed that D509 and D647 were the potential binding sites for TRPV1. Notably, molecular dynamics (MD) studies indicated that site D509 may function as an allosteric site of TRPV1. Furthermore, to validate the computer-aided predictions, we performed experimental studies, including in vitro and in vivo assays. The results of these experiments confirmed the predictions made by our computational models, providing further evidence for the mechanisms of action of San-Ao Decoction in asthma treatment. Our findings demonstrated that: i) D509 and D647 of TRPV1 are the key binding sites for the main ingredients of SAD; ii) SAD or its main ingredients significantly reduce the influx of Ca2+ through TRPV1, following the TCM principle of "Jun, Chen, Zuo, Shi"; iii) SAD shows efficiency in comprehensive in vivo validation. In conclusion, our computer-aided investigation of San-Ao Decoction in asthma treatment has provided valuable insights into the therapeutic mechanisms of this TCM formula. The combination of computational analysis and experimental validation has proven effective in enhancing our understanding of TCM and may pave the way for future discoveries in the field.

RVvictor: Virus RNA-directed molecular interactions for RNA virus infection

RNA viruses are major human pathogens that cause seasonal epidemics and occasional pandemic outbreaks. Due to the nature of their RNA genomes, it is anticipated that virus's RNA interacts with host protein (INTPRO), messenger RNA (INTmRNA), and non-coding RNA (INTncRNA) to perform their particular functions during their transcription and replication. In other words, thus, it is urgently needed to have such valuable data on virus RNA-directed molecular interactions (especially INTPROs), which are highly anticipated to attract broad research interests in the fields of RNA virus translation and replication. In this study, a new database was constructed to describe the virus RNA-directed interaction (INTPRO, INTmRNA, INTncRNA) for RNA virus (RVvictor). This database is unique in a) unambiguously characterizing the interactions between viruses RNAs and host proteins, b) providing, for the first time, the most systematic RNA-directed interaction data resources in providing clues to understand the molecular mechanisms of RNA viruses' translation, and replication, and c) in RVvictor, comprehensive enrichment analysis is conducted for each virus RNA based on its associated target genes/proteins, and the enrichment results were explicitly illustrated using various graphs. We found significant enrichment of a suite of pathways related to infection, translation, and replication, e.g., HIV infection, coronavirus disease, regulation of viral genome replication, and so on. Due to the devastating and persistent threat posed by the RNA virus, RVvictor constructed, for the first time, a possible network of cross-talk in RNA-directed interaction, which may ultimately explain the pathogenicity of RNA virus infection. The knowledge base might help develop new anti-viral therapeutic targets in the future. It's now free and publicly accessible at: https://idrblab.org/rvvictor/.

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/.

Predicting drug-induced liver injury using graph attention mechanism and molecular fingerprints

Drug-induced liver injury (DILI) is a significant issue in drug development and clinical treatment due to its potential to cause liver dysfunction or damage, which, in severe cases, can lead to liver failure or even fatality. DILI has numerous pathogenic factors, many of which remain incompletely understood. Consequently, it is imperative to devise methodologies and tools for anticipatory assessment of DILI risk in the initial phases of drug development. In this study, we present DMFPGA, a novel deep learning predictive model designed to predict DILI. To provide a comprehensive description of molecular properties, we employ a multi-head graph attention mechanism to extract features from the molecular graphs, representing characteristics at the level of compound nodes. Additionally, we combine multiple fingerprints of molecules to capture features at the molecular level of compounds. The fusion of molecular fingerprints and graph features can more fully express the properties of compounds. Subsequently, we employ a fully connected neural network to classify compounds as either DILI-positive or DILI-negative. To rigorously evaluate DMFPGA's performance, we conduct a 5-fold cross-validation experiment. The obtained results demonstrate the superiority of our method over four existing state-of-the-art computational approaches, exhibiting an average AUC of 0.935 and an average ACC of 0.934. We believe that DMFPGA is helpful for early-stage DILI prediction and assessment in drug development.

Mapping synthetic binding proteins epitopes on diverse protein targets by protein structure prediction and protein-protein docking

Synthetic binding proteins (SBPs) are a class of artificial proteins engineered from privileged protein scaffolds, which can form highly specific molecular recognition interfaces with a variety of targets. Due to the characteristics of small size, high stability, and good tissue permeability, SBPs have important applications in biomedical research, disease diagnosis and treatment. However, knowledge of SBPs epitopes on the structures of target proteins is still limited, which hinder the development of novel SBPs. In this study, based on the currently available information of SBPs and their targets, 96 pairs of interacting proteins referring to 96 representative SBPs and 80 different targets, were systemically investigated using the state-of-the-art computational modeling techniques including AlphaFold2 protein structure prediction and Rosetta protein-protein docking. As a result, 71 out of the 96 pairs were successfully docked, of which 18, 33, and 20 pairs were defined as models with high, medium, and acceptable quality, respectively. In addition, the interface information was analyzed to decipher the interaction types driven SBPs and targets recognition. Overall, this work not only provides important structural information for understanding the mechanism of action of other SBPs with same protein scaffold, but also for aiding the rational protein engineering and to design of novel SBPs with biomedical applications.

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.

Whale optimization with random contraction and Rosenbrock method for COVID-19 disease prediction

Coronavirus Disease 2019 (COVID-19), instigated by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has hugely impacted global public health. To identify and intervene in critically ill patients early, this paper proposes an efficient, intelligent prediction model based on the machine learning approach, which combines the improved whale optimization algorithm (RRWOA) with the k-nearest neighbor (KNN) classifier. In order to improve the problem that WOA is prone to fall into local optimum, an improved version named RRWOA is proposed based on the random contraction strategy (RCS) and the Rosenbrock method. To verify the capability of the proposed algorithm, RRWOA is tested against nine classical metaheuristics, nine advanced metaheuristics, and seven well-known WOA variants based on 30 IEEE CEC2014 competition functions, respectively. The experimental results in mean, standard deviation, the Friedman test, and the Wilcoxon signed-rank test are considered, proving that RRWOA won first place on 18, 24, and 25 test functions, respectively. In addition, a binary version of the algorithm, called BRRWOA, is developed for feature selection problems. An efficient prediction model based on BRRWOA and KNN classifier is proposed and compared with seven existing binary metaheuristics based on 15 datasets of UCI repositories. The experimental results show that the proposed algorithm obtains the smallest fitness value in eleven datasets and can solve combinatorial optimization problems, indicating that it still performs well in discrete cases. More importantly, the model was compared with five other algorithms on the COVID-19 dataset. The experiment outcomes demonstrate that the model offers a scientific framework to support clinical diagnostic decision-making. Therefore, RRWOA is an effectively improved optimizer with efficient value.

ANPELA: Significantly Enhanced Quantification Tool for Cytometry-Based Single-Cell Proteomics

ANPELA is widely used for quantifying traditional bulk proteomic data. Recently, there is a clear shift from bulk proteomics to the single-cell ones (SCP), for which powerful cytometry techniques demonstrate the fantastic capacity of capturing cellular heterogeneity that is completely overlooked by traditional bulk profiling. However, the in-depth and high-quality quantification of SCP data is still challenging and severely affected by the large numbers of quantification workflows and extreme performance dependence on the studied datasets. In other words, the proper selection of well-performing workflow(s) for any studied dataset is elusory, and it is urgently needed to have a significantly enhanced and accelerated tool to address this issue. However, no such tool is developed yet. Herein, ANPELA is therefore updated to its 2.0 version (https://idrblab.org/anpela/), which is unique in providing the most comprehensive set of quantification alternatives (>1000 workflows) among all existing tools, enabling systematic performance evaluation from multiple perspectives based on machine learning, and identifying the optimal workflow(s) using overall performance ranking together with the parallel computation. Extensive validation on different benchmark datasets and representative application scenarios suggest the great application potential of ANPELA in current SCP research for gaining more accurate and reliable biological insights.

Methodology to identify a gene expression signature by merging microarray datasets

A vast number of microarray datasets have been produced as a way to identify differentially expressed genes and gene expression signatures. A better understanding of these biological processes can help in the diagnosis and prognosis of diseases, as well as in the therapeutic response to drugs. However, most of the available datasets are composed of a reduced number of samples, leading to low statistical, predictive and generalization power. One way to overcome this problem is by merging several microarray datasets into a single dataset, which is typically a challenging task. Statistical methods or supervised machine learning algorithms are usually used to determine gene expression signatures. Nevertheless, statistical methods require an arbitrary threshold to be defined, and supervised machine learning methods can be ineffective when applied to high-dimensional datasets like microarrays. We propose a methodology to identify gene expression signatures by merging microarray datasets. This methodology uses statistical methods to obtain several sets of differentially expressed genes and uses supervised machine learning algorithms to select the gene expression signature. This methodology was validated using two distinct research applications: one using heart failure and the other using autism spectrum disorder microarray datasets. For the first, we obtained a gene expression signature composed of 117 genes, with a classification accuracy of approximately 98%. For the second use case, we obtained a gene expression signature composed of 79 genes, with a classification accuracy of approximately 82%. This methodology was implemented in R language and is available, under the MIT licence, at https://github.com/bioinformatics-ua/MicroGES.

scMAGS: Marker gene selection from scRNA-seq data for spatial transcriptomics studies

Single-Cell RNA sequencing (scRNA-seq) has provided unprecedented opportunities for exploring gene expression and thus uncovering regulatory relationships between genes at the single-cell level. However, scRNA-seq relies on isolating cells from tissues. Therefore, the spatial context of the regulatory processes is lost. A recent technological innovation, spatial transcriptomics, allows for the measurement of gene expression while preserving spatial information. An initial step in the spatial transcriptomic analysis is to identify the cell type, which requires a careful selection of cell-specific marker genes. For this purpose, currently, scRNA-seq data is used to select a limited number of marker genes from among all genes that distinguish cell types from each other. This study proposes scMAGS (single-cell MArker Gene Selection), a novel method for marker gene selection from scRNA-seq data for spatial transcriptomics studies. scMAGS uses a filtering step in which the candidate genes are identified before the marker gene selection step. For the selection of marker genes, cluster validity indices, the Silhouette index, or the Calinski-Harabasz index (for large datasets) are utilized. Experimental results showed that, in comparison to the existing methods, scMAGS is scalable, fast, and accurate. Even for large datasets with millions of cells, scMAGS could find the required number of marker genes in a reasonable amount of time with fewer memory requirements. scMAGS is made freely available at https://github.com/doganlab/scmags and can be downloaded from the Python Package Directory (PyPI) software repository with the command pip install scmags.

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.

Investigating cardiotoxicity related with hERG channel blockers using molecular fingerprints and graph attention mechanism

Human ether-a-go-go-related gene (hERG) channel blockade by small molecules is a big concern during drug development in the pharmaceutical industry. Failure or inhibition of hERG channel activity caused by drug molecules can lead to prolonging QT interval, which will result in serious cardiotoxicity. Thus, evaluating the hERG blocking activity of all these small molecular compounds is technically challenging, and the relevant procedures are expensive and time-consuming. In this study, we develop a novel deep learning predictive model named DMFGAM for predicting hERG blockers. In order to characterize the molecule more comprehensively, we first consider the fusion of multiple molecular fingerprint features to characterize its final molecular fingerprint features. Then, we use the multi-head attention mechanism to extract the molecular graph features. Both molecular fingerprint features and molecular graph features are fused as the final features of the compounds to make the feature expression of compounds more comprehensive. Finally, the molecules are classified into hERG blockers or hERG non-blockers through the fully connected neural network. We conduct 5-fold cross-validation experiment to evaluate the performance of DMFGAM, and verify the robustness of DMFGAM on external validation datasets. We believe DMFGAM can serve as a powerful tool to predict hERG channel blockers in the early stages of drug discovery and development.

Prediction of lung cancer metastasis by gene expression

Tumor metastasis is the main cause of death in cancer patients. Early prediction of tumor metastasis can allow for timely intervention. At present, research on tumor metastasis mainly focuses on manual diagnosis by imaging or diagnosis by computational methods. With the deterioration of the tumor, gene expression levels in blood change greatly. It is feasible to measure the transcripts of key genes to predict whether cancer will metastasize. Therefore, in this paper, we obtained gene expression data from 226 patients from TCGA. These data included 239,322 transcripts. Background screening and LASSO analysis were used to select 31 transcripts as features. Finally, a deep neural network (DNN) was used to determine whether or not lung cancer would metastasize. We compared our methods with several other methods and found that our method achieved the best precision. In addition, in a previous study, we identified 7 genes that play a vital role in lung cancer. We added those gene transcripts into the DNN and found that the AUC and AUPR of the model were increased.

A novel strategy for designing the magic shotguns for distantly related target pairs

Due to its promising capacity in improving drug efficacy, polypharmacology has emerged to be a new theme in the drug discovery of complex disease. In the process of novel multi-target drugs (MTDs) discovery, in silico strategies come to be quite essential for the advantage of high throughput and low cost. However, current researchers mostly aim at typical closely related target pairs. Because of the intricate pathogenesis networks of complex diseases, many distantly related targets are found to play crucial role in synergistic treatment. Therefore, an innovational method to develop drugs which could simultaneously target distantly related target pairs is of utmost importance. At the same time, reducing the false discovery rate in the design of MTDs remains to be the daunting technological difficulty. In this research, effective small molecule clustering in the positive dataset, together with a putative negative dataset generation strategy, was adopted in the process of model constructions. Through comprehensive assessment on 10 target pairs with hierarchical similarity-levels, the proposed strategy turned out to reduce the false discovery rate successfully. Constructed model types with much smaller numbers of inhibitor molecules gained considerable yields and showed better false-hit controllability than before. To further evaluate the generalization ability, an in-depth assessment of high-throughput virtual screening on ChEMBL database was conducted. As a result, this novel strategy could hierarchically improve the enrichment factors for each target pair (especially for those distantly related/unrelated target pairs), corresponding to target pair similarity-levels.

bSRWPSO-FKNN: A boosted PSO with fuzzy K-nearest neighbor classifier for predicting atopic dermatitis disease

Introduction

Atopic dermatitis (AD) is an allergic disease with extreme itching that bothers patients. However, diagnosing AD depends on clinicians' subjective judgment, which may be missed or misdiagnosed sometimes.

Methods

This paper establishes a medical prediction model for the first time on the basis of the enhanced particle swarm optimization (SRWPSO) algorithm and the fuzzy K-nearest neighbor (FKNN), called bSRWPSO-FKNN, which is practiced on a dataset related to patients with AD. In SRWPSO, the Sobol sequence is introduced into particle swarm optimization (PSO) to make the particle distribution of the initial population more uniform, thus improving the population's diversity and traversal. At the same time, this study also adds a random replacement strategy and adaptive weight strategy to the population updating process of PSO to overcome the shortcomings of poor convergence accuracy and easily fall into the local optimum of PSO. In bSRWPSO-FKNN, the core of which is to optimize the classification performance of FKNN through binary SRWPSO.

Results

To prove that the study has scientific significance, this paper first successfully demonstrates the core advantages of SRWPSO in well-known algorithms through benchmark function validation experiments. Secondly, this article demonstrates that the bSRWPSO-FKNN has practical medical significance and effectiveness through nine public and medical datasets.

Discussion

The 10 times 10-fold cross-validation experiments demonstrate that bSRWPSO-FKNN can pick up the key features of AD, including the content of lymphocytes (LY), Cat dander, Milk, Dermatophagoides Pteronyssinus/Farinae, Ragweed, Cod, and Total IgE. Therefore, the established bSRWPSO-FKNN method practically aids in the diagnosis of AD.

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.

Eres-UNet++: Liver CT image segmentation based on high-efficiency channel attention and Res-UNet+

Computerized tomography (CT) is of great significance for the localization and diagnosis of liver cancer. Many scholars have recently applied deep learning methods to segment CT images of liver and liver tumors. Unlike natural images, medical image segmentation is usually more challenging due to its nature. Aiming at the problem of blurry boundaries and complex gradients of liver tumor images, a deep supervision network based on the combination of high-efficiency channel attention and Res-UNet++ (ECA residual UNet++) is proposed for liver CT image segmentation, enabling fully automated end-to-end segmentation of the network. In this paper, the UNet++ structure is selected as the baseline. The residual block feature encoder based on context awareness enhances the feature extraction ability and solves the problem of deep network degradation. The introduction of an efficient attention module combines the depth of the feature map with spatial information to alleviate the uneven sample distribution impact; Use DiceLoss to replace the cross-entropy loss function to optimize network parameters. The liver and liver tumor segmentation accuracy on the LITS dataset was 95.8% and 89.3%, respectively. The results show that compared with other algorithms, the method proposed in this paper achieves a good segmentation performance, which has specific reference significance for computer-assisted diagnosis and treatment to attain fine segmentation of liver and liver tumors.

Tuberculous pleural effusion prediction using ant colony optimizer with grade-based search assisted support vector machine

Introduction

Although tuberculous pleural effusion (TBPE) is simply an inflammatory response of the pleura caused by tuberculosis infection, it can lead to pleural adhesions and cause sequelae of pleural thickening, which may severely affect the mobility of the chest cavity.

Methods

In this study, we propose bGACO-SVM, a model with good diagnostic power, for the adjunctive diagnosis of TBPE. The model is based on an enhanced continuous ant colony optimization (ACOR) with grade-based search technique (GACO) and support vector machine (SVM) for wrapped feature selection. In GACO, grade-based search greatly improves the convergence performance of the algorithm and the ability to avoid getting trapped in local optimization, which improves the classification capability of bGACO-SVM.

Results

To test the performance of GACO, this work conducts comparative experiments between GACO and nine basic algorithms and nine state-of-the-art variants as well. Although the proposed GACO does not offer much advantage in terms of time complexity, the experimental results strongly demonstrate the core advantages of GACO. The accuracy of bGACO-predictive SVM was evaluated using existing datasets from the UCI and TBPE datasets.

Discussion

In the TBPE dataset trial, 147 TBPE patients were evaluated using the created bGACO-SVM model, showing that the bGACO-SVM method is an effective technique for accurately predicting TBPE.

An efficient multi-threshold image segmentation for skin cancer using boosting whale optimizer

Due to the terrible manifestations of skin cancer, it seriously disturbs the quality of life status and health of patients, so we needs treatment plans to detect it early and avoid it causing more harm to patients. Medical disease image threshold segmentation technique can well extract the region of interest and effectively assist in disease recognition. Moreover, in multi-threshold image segmentation, the selection of the threshold set determines the image segmentation quality. Among the common threshold selection methods, the selection based on metaheuristic algorithm has the advantages of simplicity, easy implementation and avoidable local optimization. However, different algorithms have different performances for different medical disease images. For example, the Whale Optimization Algorithm (WOA) does not give a satisfactory performance for thresholding skin cancer images. We propose an improved WOA (LCWOA) in which the Levy operator and chaotic random mutation strategy are introduced to enhance the ability of the algorithm to jump out of the local optimum and to explore the search space. Comparing with different existing WOA variants on the CEC2014 function set, our proposed and improved algorithm improves the efficiency of the search. Experimental results show that our method outperforms the extant WOA variants in terms of optimization performances, improving the convergence accuracy and velocity. The method is also applied to solve the threshold selection in the skin cancer image segmentation problem, and LCWOA also gives excellent performance in obtaining optimal segmentation results.

In silico analysis of a novel protein in CAR T- cell therapy for the treatment of hematologic cancer through molecular modelling, docking, and dynamics approach

Cellular therapy has emerged as a key tool in the treatment of hematological malignancies. An advanced cell therapy known as chimeric antigen receptor T cell therapy (CAR T-cell therapy) has been approved by the United States Food and Drug Administration (FDA) as KYMRIAH by Novartis and YESCARTA by Gilead/Kite pharma in the year 2017. A chimeric receptor is composed of an extracellular antigen recognition site along with some co-stimulating and signalling domains. On the whole, it turns out to be one of the most potent receptors on T cells targeting a specific type of cancer cell based on its antigenic marker. CD19 CAR T-cell therapy is the first clinically approved therapy for lymphoma with remarkable results in complete remission of B cell lymphoblastic leukemia up to 90%. The high rate of effectiveness of the CAR T-cell therapy against B-ALL justifies the investigation and application of this therapy for fatal diseases like all types of hematological malignancies. The most critical aspect of chimeric receptor therapy is designing and building an artificial receptor that is specific to a given type of cancer. For this reason, the in silico technique is an appropriate model to investigate the integrity and effectiveness of the engineered chimeric receptor prior to commencing in vitro experiments followed by clinical trials. This computerized experimental study aids in predicting the molecular mechanism of chimeric protein and how it interacts with both ligands. We have anticipated various features of the chimeric protein in terms of qualitative analysis (structure, protein modelling, physiological properties) and functional analysis (antigenicity, allergenicity, its receptor-ligand binding ability, involving signalling pathways). Furthermore, the reliability and validation of the binding mode of the chimeric protein against receptors were performed through a complex molecular dynamics simulation for a 100 ns timeframe in an aqueous environment. The obtained simulation study showed that CD30 was a better approachable marker as compared to CD20 due to its better binding energy score and also binding conformations stability.

Multi-threshold image segmentation for melanoma based on Kapur’s entropy using enhanced ant colony optimization

Melanoma is a malignant tumor formed by the cancerous transformation of melanocytes, and its medical images contain much information. However, the percentage of the critical information in the image is small, and the noise is non-uniformly distributed. We propose a new multi-threshold image segmentation model based on the two-dimensional histogram approach to the above problem. We present an enhanced ant colony optimization for continuous domains (EACOR) in the proposed model based on the soft besiege and chase strategies. Further, EACOR is combined with two-dimensional Kapur’s entropy to search for the optimal thresholds. An experiment on the IEEE CEC2014 benchmark function was conducted to measure the reliable global search capability of the EACOR algorithm in the proposed model. Moreover, we have also conducted several sets of experiments to test the validity of the image segmentation model proposed in this paper. The experimental results show that the segmented images from the proposed model outperform the comparison method in several evaluation metrics. Ultimately, the model proposed in this paper can provide high-quality samples for subsequent analysis of melanoma pathology images.

GCN-GENE: A novel method for prediction of coronary heart disease-related genes

Coronary heart disease is the most common heart disease, it can induce myocardial infarction, and the cause of the disease has a lot to do with life and eating habits. The results of a large number of epidemiological studies at home and abroad show that the incidence of coronary heart disease has an obvious familial tendency. However, little is known about the genetic factors of coronary heart disease. Although genome-wide association analysis and gene knockout experiments have found some genes related to coronary heart disease, there are still a large number of genes potentially related to coronary heart disease that have not been discovered. If it is confirmed by biological experimental means, the time and money cost is too high. Therefore, it is urgent to identify genes related to coronary heart disease on a large scale by computational means, so as to conduct targeted biological experimental verification. This paper proposes a deep learning method based on biological networks for the identification of coronary heart disease-related genes. We constructed gene interaction networks and extracted gene expression levels in different tissues as features. Through the association information and expression characteristics between genes, we constructed a model of coronary heart disease-related genes. Through cross-validation, we found that our proposed GCN-GENE that has AUC as 0.75 and AUPR as 0.78, which is more accurate than other methods and is a reliable method for predicting coronary heart disease-related genes.

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/.

CovInter: interaction data between coronavirus RNAs and host proteins

Coronavirus has brought about three massive outbreaks in the past two decades. Each step of its life cycle invariably depends on the interactions among virus and host molecules. The interaction between virus RNA and host protein (IVRHP) is unique compared to other virus-host molecular interactions and represents not only an attempt by viruses to promote their translation/replication, but also the host's endeavor to combat viral pathogenicity. In other words, there is an urgent need to develop a database for providing such IVRHP data. In this study, a new database was therefore constructed to describe the interactions between coronavirus RNAs and host proteins (CovInter). This database is unique in (a) unambiguously characterizing the interactions between virus RNA and host protein, (b) comprehensively providing experimentally validated biological function for hundreds of host proteins key in viral infection and (c) systematically quantifying the differential expression patterns (before and after infection) of these key proteins. Given the devastating and persistent threat of coronaviruses, CovInter is highly expected to fill the gap in the whole process of the 'molecular arms race' between viruses and their hosts, which will then aid in the discovery of new antiviral therapies. It's now free and publicly accessible at: https://idrblab.org/covinter/.

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/.

A novel multi-class classification model for schizophrenia, bipolar disorder and healthy controls using comprehensive transcriptomic data

Two common psychiatric disorders, schizophrenia (SCZ) and bipolar disorder (BP), confer lifelong disability and collectively affect 2% of the world population. Because the diagnosis of psychiatry is based only on symptoms, developing more effective methods for the diagnosis of psychiatric disorders is a major international public health priority. Furthermore, SCZ and BP overlap considerably in terms of symptoms and risk genes. Therefore, the clarity of the underlying etiology and pathology remains lacking for these two disorders. Although many studies have been conducted, a classification model with higher accuracy and consistency was found to still be necessary for accurate diagnoses of SCZ and BP. In this study, a comprehensive dataset was combined from five independent transcriptomic studies. This dataset comprised 120 patients with SCZ, 101 patients with BP, and 149 healthy subjects. The partial least squares discriminant analysis (PLS-DA) method was applied to identify the gene signature among multiple groups, and 341 differentially expressed genes (DEGs) were identified. Then, the disease relevance of these DEGs was systematically performed, including (α) the great disease relevance of the identified signature, (β) the hub genes of the protein-protein interaction network playing a key role in psychiatric disorders, and (γ) gene ontology terms and enriched pathways playing a key role in psychiatric disorders. Finally, a popular multi-class classifier, support vector machine (SVM), was applied to construct a novel multi-class classification model using the identified signature for SCZ and BP. Using the independent test sets, the classification capacity of this multi-class model was assessed, which showed this model had a strong classification ability.

Mutational Slime Mould Algorithm for Gene Selection

A large volume of high-dimensional genetic data has been produced in modern medicine and biology fields. Data-driven decision-making is particularly crucial to clinical practice and relevant procedures. However, high-dimensional data in these fields increase the processing complexity and scale. Identifying representative genes and reducing the data’s dimensions is often challenging. The purpose of gene selection is to eliminate irrelevant or redundant features to reduce the computational cost and improve classification accuracy. The wrapper gene selection model is based on a feature set, which can reduce the number of features and improve classification accuracy. This paper proposes a wrapper gene selection method based on the slime mould algorithm (SMA) to solve this problem. SMA is a new algorithm with a lot of application space in the feature selection field. This paper improves the original SMA by combining the Cauchy mutation mechanism with the crossover mutation strategy based on differential evolution (DE). Then, the transfer function converts the continuous optimizer into a binary version to solve the gene selection problem. Firstly, the continuous version of the method, ISMA, is tested on 33 classical continuous optimization problems. Then, the effect of the discrete version, or BISMA, was thoroughly studied by comparing it with other gene selection methods on 14 gene expression datasets. Experimental results show that the continuous version of the algorithm achieves an optimal balance between local exploitation and global search capabilities, and the discrete version of the algorithm has the highest accuracy when selecting the least number of genes.

Uncovering the dynamic effects of DEX treatment on lung cancer by integrating bioinformatic inference and multiscale modeling of scRNA-seq and proteomics data

Lung cancer is one of the leading causes of cancer-related death, with a five-year survival rate of 18%. It is a priority for us to understand the underlying mechanisms affecting lung cancer therapeutics' implementation and effectiveness. In this study, we combine the power of Bioinformatics and Systems Biology to comprehensively uncover functional and signaling pathways of drug treatment using bioinformatics inference and multiscale modeling of both scRNA-seq data and proteomics data. Based on a time series of lung adenocarcinoma derived A549 cells after DEX treatment, we first identified the differentially expressed genes (DEGs) in those lung cancer cells. Through the interrogation of regulatory network of those DEGs, we identified key hub genes including TGFβ, MYC, and SMAD3 varied underlie DEX treatment. Further gene set enrichment analysis revealed the TGFβ signaling pathway as the top enriched term. Those genes involved in the TGFβ pathway and their crosstalk with the ERBB pathway presented a strong survival prognosis in clinical lung cancer samples. With the basis of biological validation and literature-based curation, a multiscale model of tumor regulation centered on both TGFβ-induced and ERBB-amplified signaling pathways was developed to characterize the dynamic effects of DEX therapy on lung cancer cells. Our simulation results were well matched to available data of SMAD2, FOXO3, TGFβ1, and TGFβR1 over the time course. Moreover, we provided predictions of different doses to illustrate the trend and therapeutic potential of DEX treatment. The innovative and cross-disciplinary approach can be further applied to other computational studies in tumorigenesis and oncotherapy. We released the approach as a user-friendly tool named BIMM (Bioinformatic Inference and Multiscale Modeling), with all the key features available at https://github.com/chenm19/BIMM.

Investigating the structure-activity relationship of marine polycyclic batzelladine alkaloids as promising inhibitors for SARS-CoV-2 main protease (Mpro)

Over a span of two years ago, since the emergence of the first case of the novel coronavirus (SARS-CoV-2) in China, the pandemic has crossed borders causing serious health emergencies, immense economic crisis and impacting the daily life worldwide. Despite the discovery of numerous forms of precautionary vaccines along with other recently approved orally available drugs, yet effective antiviral therapeutics are necessarily needed to hunt this virus and its variants. Historically, naturally occurring chemicals have always been considered the primary source of beneficial medications. Considering the SARS-CoV-2 main protease (M^pro) as the duplicate key element of the viral cycle and its main target, in this paper, an extensive virtual screening for a focused chemical library of 15 batzelladine marine alkaloids, was virtually examined against SARS-CoV-2 main protease (M^pro) using an integrated set of modern computational tools including molecular docking (MDock), molecule dynamic (MD) simulations and structure-activity relationships (SARs) as well. The molecular docking predictions had disclosed four promising compounds including batzelladines H–I (8–9) and batzelladines F-G (6–7), respectively according to their prominent ligand-protein energy scores and relevant binding affinities with the (M^pro) pocket residues. The best two chemical hits, batzelladines H–I (8–9) were further investigated thermodynamically though studying their MD simulations at 100 ns, where they showed excellent stability within the accommodated (M^pro) pocket. Moreover, SARs studies imply the crucial roles of the fused tricyclic guanidinic moieties, its degree of unsaturation, position of the N–OH functionality and the length of the side chain as a spacer linking between two active sites, which disclosed fundamental structural and pharmacophoric features for efficient protein-ligand interaction. Such interesting findings are greatly highlighting further in vitro/vivo examinations regarding those marine natural products (MNPs) and their synthetic equivalents as promising antivirals.

Identification of chromosomal instability-associated genes as hepatocellular carcinoma progression-related biomarkers to guide clinical diagnosis, prognosis and therapy

Hepatocellular carcinoma (HCC) is a type of cancer characterized by high heterogeneity and a complex multistep progression process. Significantly-altered biomarkers for HCC need to be identified. Differentially expressed genes and weighted gene co-expression network analyses were used to identify progression-related biomarkers. LASSO-Cox regression and random forest algorithms were used to construct the progression-related prognosis (PRP) score. Three chromosomal instability-associated genes (KIF20A, TOP2A, and TTK) have been identified as progression-related biomarkers. The robustness of the PRP scores were validated using four independent cohorts. Immune status was observed using the single-sample gene set enrichment analysis (ssGSEA). Comprehensive analysis showed that the patients with high PRP score had wider genomic alterations, more malignant phenotypes, and were in a state of immunosuppression. The diagnostic models constructed via logistic regression based on the three genes showed satisfactory performances in distinguishing HCC from cirrhotic tissues or dysplastic nodules. The nomogram combining PRP scores with clinical factors had a better performance in predicting prognosis than the tumor node metastasis classification (TNM) system. We further confirmed that KIF20A, TOP2A, and TTK were highly expressed in HCC tissues than in cirrhotic tissues. Downregulation of all three genes aggravated chromosomal instabilities in HCC and suppressed HCC cells viability both in vitro and in vivo. Overall, our study highlights the important roles of chromosomal instability-associated genes during the progression of HCC and their potential clinical diagnosis and prognostic value and provides promising new ideas for developing therapeutic strategies to improve the outcomes of HCC patients.

Predicting Entrepreneurial Intention of Students: Kernel Extreme Learning Machine with Boosted Crow Search Algorithm

College students are the group with the most entrepreneurial vitality and potential. How to cultivate their entrepreneurial and innovative ability is one of the important and urgent issues facing this current social development. This paper proposes a reliable, intelligent prediction model of entrepreneurial intentions, providing theoretical support for guiding college students’ positive entrepreneurial intentions. The model mainly uses the improved crow search algorithm (CSA) to optimize the kernel extreme learning machine (KELM) model with feature selection (FS), namely CSA-KELM-FS, to study entrepreneurial intention. To obtain the best fitting model and key features, the gradient search rule, local escaping operator, and levy flight mutation (GLL) mechanism are introduced to enhance the CSA (GLLCSA), and FS is used to extract the key features. To verify the performance of the proposed GLLCSA, it is compared with eight other state-of-the-art methods. Further, the GLLCSA-KELM-FS model and five other machine learning methods have been used to predict the entrepreneurial intentions of 842 students from the Wenzhou Vocational College in Zhejiang, China, in the past five years. The results show that the proposed model can correctly predict the students’ entrepreneurial intention with an accuracy rate of 93.2% and excellent stability. According to the prediction results of the proposed model, the key factors affecting the student’s entrepreneurial intention are mainly the major studied, campus innovation, entrepreneurship practice experience, and positive personality. Therefore, the proposed GLLCSA-KELM-FS is expected to be an effective tool for predicting students’ entrepreneurial intentions.

Multilevel threshold image segmentation for COVID-19 chest radiography: A framework using horizontal and vertical multiverse optimization

COVID-19 is currently raging worldwide, with more patients being diagnosed every day. It usually is diagnosed by examining pathological photographs of the patient's lungs. There is a lot of detailed and essential information on chest radiographs, but manual processing is not as efficient or accurate. As a result, how efficiently analyzing and processing chest radiography of COVID-19 patients is an important research direction to promote COVID-19 diagnosis. To improve the processing efficiency of COVID-19 chest films, a multilevel thresholding image segmentation (MTIS) method based on an enhanced multiverse optimizer (CCMVO) is proposed. CCMVO is improved from the original Multi-Verse Optimizer by introducing horizontal and vertical search mechanisms. It has a more assertive global search ability and can jump out of the local optimum in optimization. The CCMVO-based MTIS method can obtain higher quality segmentation results than HHO, SCA, and other forms and is less prone to stagnation during the segmentation process. To verify the performance of the proposed CCMVO algorithm, CCMVO is first compared with DE, MVO, and other algorithms by 30 benchmark functions; then, the proposed CCMVO is applied to image segmentation of COVID-19 chest radiography; finally, this paper verifies that the combination of MTIS and CCMVO is very successful with good segmentation results by using the Feature Similarity Index (FSIM), the Peak Signal to Noise Ratio (PSNR), and the Structural Similarity Index (SSIM). Therefore, this research can provide an effective segmentation method for a medical organization to process COVID-19 chest radiography and then help doctors diagnose coronavirus pneumonia (COVID-19).

Curcumin inhibits spike protein of new SARS-CoV-2 variant of concern (VOC) Omicron, an in silico study

BACKGROUND

Omicron (B.1.1.529), a variant of SARS-CoV-2 is currently spreading globally as a dominant strain. Due to multiple mutations at its Spike protein, including 15 amino acid substitutions at the receptor binding domain (RBD), Omicron is a variant of concern (VOC) and capable of escaping vaccine generated immunity. So far, no specific treatment regime is suggested for this VOC.

METHODS

The three-dimensional structure of the Spike RBD domain of Omicron variant was constructed by incorporating 15 amino acid substitutions to the Native Spike (S) structure and structural changes were compared that of the Native S. Seven phytochemicals namely Allicin, Capsaicin, Cinnamaldehyde, Curcumin, Gingerol, Piperine, and Zingeberene were docked with Omicron S protein and Omicron S-hACE2 complex. Further, molecular dynamic simulation was performed between Crcumin and Omicron S protein to evaluate the structural stability of the complex in the physiological environment and compared with that of the control drug Chloroquine.

RESULTS

Curcumin, among seven phytochemicals, was found to have the most substantial inhibitory potential with Omicron S protein. Further, it was found that curcumin could disrupt the Omicron S-hACE2 complex. The molecular dynamic simulation demonstrated that Curcumin could form a stable structure with Omicron S in the physiological environment.

CONCLUSION

To conclude, Curcumin can be considered as a potential therapeutic agent against the highly infectious Omicron variant of SARS-CoV-2.