RNA-RNA interactions between SARS-CoV-2 and host benefit viral development and evolution during COVID-19 infection

OncoSexome: the landscape of sex-based differences in oncologic diseases

The NIH policy on sex as biological variable (SABV) emphasized the importance of sex-based differences in precision oncology. Over 50% of clinically actionable oncology genes are sex-biased, indicating differences in drug efficacy. Research has identified sex differences in non-reproductive cancers, highlighting the need for comprehensive sex-based cancer data. We therefore developed OncoSexome, a multidimensional knowledge base describing sex-based differences in cancer (https://idrblab.org/OncoSexome/) across four key topics: antineoplastic drugs and responses (SDR), oncology-related biomarkers (SBM), risk factors (SRF) and microbial landscape (SML). SDR covers sex-based differences in 2051 anticancer drugs; SBM describes 12 551 sex-differential biomarkers; SRF illustrates 350 sex-dependent risk factors; SML demonstrates 1386 microbes with sex-differential abundances associated with cancer development. OncoSexome is unique in illuminating multifaceted influences of biological sex on cancer, providing both external and endogenous contributors to cancer development and describing sex-based differences for the broadest oncological classes. Given the increasing global research interest in sex-based differences, OncoSexome is expected to impact future precision oncology practices significantly.

SYNBIP 2.0: epitopes mapping, sequence expansion and scaffolds discovery for synthetic binding protein innovation

Synthetic binding proteins (SBPs) represent a pivotal class of artificially engineered proteins, meticulously crafted to exhibit targeted binding properties and specific functions. Here, the SYNBIP database, a comprehensive resource for SBPs, has been significantly updated. These enhancements include (i) featuring 3D structures of 899 SBP-target complexes to illustrate the binding epitopes of SBPs, (ii) using the structures of SBPs in the monomer or complex forms with target proteins, their sequence space has been expanded five times to 12 025 by integrating a structure-based protein generation framework and a protein property prediction tool, (iii) offering detailed information on 78 473 newly identified SBP-like scaffolds from the RCSB Protein Data Bank, and an additional 16 401 555 ones from the AlphaFold Protein Structure Database, and (iv) the database is regularly updated, incorporating 153 new SBPs. Furthermore, the structural models of all SBPs have been enhanced through the application of the AlphaFold2, with their clinical statuses concurrently refreshed. Additionally, the design methods employed for each SBP are now prominently featured in the database. In sum, SYNBIP 2.0 is designed to provide researchers with essential SBP data, facilitating their innovation in research, diagnosis and therapy. SYNBIP 2.0 is now freely accessible at https://idrblab.org/synbip/.

OrgXenomics: an integrated proteomic knowledge base for patient-derived organoid and xenograft

Patient-derived models (PDMs, particularly organoids and xenografts) are irreplaceable tools for precision medicine, from target development to lead identification, then to preclinical evaluation, and finally to clinical decision-making. So far, PDM-based proteomics has emerged to be one of the cutting-edge directions and massive data have been accumulated. However, such PDM-based proteomic data have not been provided by any of the available databases, and proteomics profiles of all proteins in proteomic study are also completely absent from existing databases. Herein, an integrated database named 'OrgXenomics' was thus developed to provide the proteomic data for PDMs, which was unique in (a) explicitly describing the establishment detail for a wide array of models, (b) systematically providing the proteomic profiles (expression/function/interaction) for all proteins in studied proteomic analysis and (c) comprehensively giving the raw data for diverse organoid/xenograft-based proteomic studies of various diseases. Our OrgXenomics was expected to server as one good complement to existing proteomic databases, and had great implication for the practice of precision medicine, which could be accessed at: https://idrblab.org/orgxenomics/.

SubCELL: the landscape of subcellular compartment-specific molecular interactions

The subcellular compartment-specific molecular interactions (SCSIs) are the building blocks for most molecular functions, biological processes and disease pathogeneses. Extensive experiments have therefore been conducted to accumulate the valuable information of SCSIs, but none of the available databases has been constructed to describe those data. In this study, a novel knowledge base SubCELL is thus introduced to depict the landscape of SCSIs among DNAs/RNAs/proteins. This database is UNIQUE in (a) providing, for the first time, the experimentally-identified SCSIs, (b) systematically illustrating a large number of SCSIs inferred based on well-established method and (c) collecting experimentally-determined subcellular locations for the DNAs/RNAs/proteins of diverse species. Given the essential physiological/pathological role of SCSIs, the SubCELL is highly expected to have great implications for modern molecular biological study, which can be freely accessed with no login requirement at: https://idrblab.org/subcell/.

MolBiC: the cell-based landscape illustrating molecular bioactivities

The measurement of cell-based molecular bioactivity (CMB) is critical for almost every step of drug development. With the booming application of AI in biomedicine, it is essential to have the CMB data to promote the learning of cell-based patterns for guiding modern drug discovery, but no database providing such information has been constructed yet. In this study, we introduce MolBiC, a knowledge base designed to describe valuable data on molecular bioactivity measured within a cellular context. MolBiC features 550 093 experimentally validated CMBs, encompassing 321 086 molecules and 2666 targets across 988 cell lines. Our MolBiC database is unique in describing the valuable data of CMB, which meets the critical demands for CMB-based big data promoting the learning of cell-based molecular/pharmaceutical pattern in drug discovery and development. MolBiC is now freely accessible without any login requirement at: https://idrblab.org/MolBiC/.

SARS-CoV-2: Prediction of critical ionic amino acid mutations

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), that caused coronavirus disease 2019 (COVID-19), has been studied thoroughly, and several variants are revealed across the world with their corresponding mutations. Studies and vaccines development focus on the genetic mutations of the S protein due to its vital role in allowing the virus attach and fuse with the membrane of a host cell. In this perspective, we study the effects of all ionic amino acid mutations of the SARS-CoV-2 viral spike protein S1 when bound to Antibody CC12.1 within the SARS-CoV-2:CC12.1 complex model. Binding free energy calculations between SARS-CoV-2 and antibody CC12.1 are based on the Analysis of Electrostatic Similarities of Proteins (AESOP) framework, where the electrostatic potentials are calculated using Adaptive Poisson-Boltzmann Solver (APBS). The atomic radii and charges that feed into the APBS calculations are calculated using the PDB2PQR software. Our results are the first to propose in silico potential life-threatening mutations of SARS-CoV-2 beyond the present mutations found in the five common variants worldwide. We find each of the following mutations: K378A, R408A, K424A, R454A, R457A, K458A, and K462A, to play significant roles in the binding to Antibody CC12.1, since they are turned into strong inhibitors on both chains of the S1 protein, whereas the mutations D405A, D420A, and D427A, show to play important roles in this binding, as they are turned into mild inhibitors on both chains of the S1 protein.

Towards understanding post-COVID-19 condition: A systematic meta-analysis of transcriptomic alterations with sex-specific insights

BACKGROUND

Post COVID-19 Condition (PCC), characterized by lingering symptoms post-acute COVID-19, poses clinical challenges, highlighting the need to understand its underlying molecular mechanisms. This meta-analysis aims to shed light on the transcriptomic landscapes and sex-specific molecular dynamics intrinsic to PCC.

METHODS

A systematic review identified three studies suitable for comprehensive meta-analysis, encompassing 135 samples (57 PCC subjects and 78 recovered subjects). We performed meta-analysis on differential gene expression, a gene set enrichment analysis of Reactome pathways, and weighted gene co-expression network analysis (WGCNA). We performed a drug and disease enrichment analysis and also assessed sex-specific differences in expression patterns.

KEY FINDINGS

A clear difference was observed in the transcriptomic profiles of PCC subjects, with 530 differentially expressed genes (DEGs) identified. Enrichment analysis revealed that the altered pathways were predominantly implicated in cell cycle processes, immune dysregulation and histone modifications. Antioxidant compounds such as hesperitin were predominantly linked to the hub genes of the DEGs. Sex-specific analyses highlighted disparities in DEGs and altered pathways in male and female PCC patients, revealing a difference in the expression of ribosomal proteins. PCC in men was mostly linked to neuro-cardiovascular disorders, while women exhibited more diverse disorders, with a high index of respiratory conditions.

CONCLUSION

Our study reveals the intricate molecular processes underlying PCC, highlighting that the differences in molecular dynamics between males and females could be key to understanding and effectively managing the varied symptomatology of this condition.

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

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

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.

SUnet: A multi-organ segmentation network based on multiple attention

Organ segmentation in abdominal or thoracic computed tomography (CT) images plays a crucial role in medical diagnosis as it enables doctors to locate and evaluate organ abnormalities quickly, thereby guiding surgical planning, and aiding treatment decision-making. This paper proposes a novel and efficient medical image segmentation method called SUnet for multi-organ segmentation in the abdomen and thorax. SUnet is a fully attention-based neural network. Firstly, an efficient spatial reduction attention (ESRA) module is introduced not only to extract image features better, but also to reduce overall model parameters, and to alleviate overfitting. Secondly, SUnet's multiple attention-based feature fusion module enables effective cross-scale feature integration. Additionally, an enhanced attention gate (EAG) module is considered by using grouped convolution and residual connections, providing richer semantic features. We evaluate the performance of the proposed model on synapse multiple organ segmentation dataset and automated cardiac diagnostic challenge dataset. SUnet achieves an average Dice of 84.29% and 92.25% on these two datasets, respectively, outperforming other models of similar complexity and size, and achieving state-of-the-art results.

A task-specific encoding algorithm for RNAs and RNA-associated interactions based on convolutional autoencoder

RNAs play essential roles in diverse physiological and pathological processes by interacting with other molecules (RNA/protein/compound), and various computational methods are available for identifying these interactions. However, the encoding features provided by existing methods are limited and the existing tools does not offer an effective way to integrate the interacting partners. In this study, a task-specific encoding algorithm for RNAs and RNA-associated interactions was therefore developed. This new algorithm was unique in (a) realizing comprehensive RNA feature encoding by introducing a great many of novel features and (b) enabling task-specific integration of interacting partners using convolutional autoencoder-directed feature embedding. Compared with existing methods/tools, this novel algorithm demonstrated superior performances in diverse benchmark testing studies. This algorithm together with its source code could be readily accessed by all user at: https://idrblab.org/corain/ and https://github.com/idrblab/corain/.

scAAGA: Single cell data analysis framework using asymmetric autoencoder with gene attention

In recent years, single-cell RNA sequencing (scRNA-seq) has emerged as a powerful technique for investigating cellular heterogeneity and structure. However, analyzing scRNA-seq data remains challenging, especially in the context of COVID-19 research. Single-cell clustering is a key step in analyzing scRNA-seq data, and deep learning methods have shown great potential in this area. In this work, we propose a novel scRNA-seq analysis framework called scAAGA. Specifically, we utilize an asymmetric autoencoder with a gene attention module to learn important gene features adaptively from scRNA-seq data, with the aim of improving the clustering effect. We apply scAAGA to COVID-19 peripheral blood mononuclear cell (PBMC) scRNA-seq data and compare its performance with state-of-the-art methods. Our results consistently demonstrate that scAAGA outperforms existing methods in terms of adjusted rand index (ARI), normalized mutual information (NMI), and adjusted mutual information (AMI) scores, achieving improvements ranging from 2.8% to 27.8% in NMI scores. Additionally, we discuss a data augmentation technology to expand the datasets and improve the accuracy of scAAGA. Overall, scAAGA presents a robust tool for scRNA-seq data analysis, enhancing the accuracy and reliability of clustering results in COVID-19 research.

In-depth characterization and identification of translatable lncRNAs

Long non-coding RNAs (LncRNAs) are non-protein coding transcripts more than 200 nucleotides in length. Deep sequencing technologies have unveiled lncRNAs can harbor translatable short open reading frames (sORFs). Yet the regulatory mechanisms governing lncRNA translation events remain poorly understood. Here, we exhaustively detected the sequence, functional element, and structure features relevant to lncRNA translation in human. Extensive identification and analysis reveal that translatable lncRNAs contain richer protein-coding related sequence features, cap-dependent and cap-independent translation initiation mechanisms, and more stable secondary structures, as compared to untranslatable lncRNAs. These findings strongly support lncRNAs serve as a repository for the production of new small peptides. Based on the feature fusion affecting translation and the extreme gradient boosting (XGBoost) algorithm, we developed the first computational tool that dedicated for predicting translatable lncRNAs, named TransLncPred. Benchmark experimental results show that our method outperforms several state-of-the-art RNA coding potential prediction tools on the same training and testing datasets. The 100-time 10-fold cross-validation tests also demonstrate that regulatory element-derived features, especially N7-methylguanosine (m7G) and internal ribosome entry site (IRES), contribute to the improvement in predictive performance.

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.

Computational Prediction of RNA-RNA Interactions between Small RNA Tracks from Betacoronavirus Nonstructural Protein 3 and Neurotrophin Genes during Infection of an Epithelial Lung Cancer Cell Line: Potential Role of Novel Small Regulatory RNA

Whether RNA-RNA interactions of cytoplasmic RNA viruses, such as Betacoronavirus, might end in the biogenesis of putative virus-derived small RNAs as miRNA-like molecules has been controversial. Even more, whether RNA-RNA interactions of wild animal viruses may act as virus-derived small RNAs is unknown. Here, we address these issues in four ways. First, we use conserved RNA structures undergoing negative selection in the genomes of SARS-CoV, MERS-CoV, and SARS-CoV-2 circulating in different bat species, intermediate animals, and human hosts. Second, a systematic literature review was conducted to identify Betacoronavirus-targeting hsa-miRNAs involved in lung cell infection. Third, we employed sophisticated long-range RNA-RNA interactions to refine the seed sequence homology of hsa-miRNAs with conserved RNA structures. Fourth, we used high-throughput RNA sequencing of a Betacoronavirus-infected epithelial lung cancer cell line (Calu-3) to validate the results. We proposed nine potential virus-derived small RNAs: two vsRNAs in SARS-CoV (Bats: SB-vsRNA-ORF1a-3p; SB-vsRNA-S-5p), one vsRNA in MERS-CoV (Bats: MB-vsRNA-ORF1b-3p), and six vsRNAs in SARS-CoV-2 (Bats: S2B-vsRNA-ORF1a-5p; intermediate animals: S2I-vsRNA-ORF1a-5p; and humans: S2H-vsRNA-ORF1a-5p, S2H-vsRNA-ORF1a-3p, S2H-vsRNA-ORF1b-3p, S2H-vsRNA-ORF3a-3p), mainly encoded by nonstructural protein 3. Notably, Betacoronavirus-derived small RNAs targeted 74 differentially expressed genes in infected human cells, of which 55 upregulate the molecular mechanisms underlying acute respiratory distress syndrome (ARDS), and the 19 downregulated genes might be implicated in neurotrophin signaling impairment. These results reveal a novel small RNA-based regulatory mechanism involved in neuropathogenesis that must be further studied to validate its therapeutic use.

In Silico Identification of Potential Inhibitors of the SARS-CoV-2 Main Protease among a PubChem Database of Avian Infectious Bronchitis Virus 3CLPro Inhibitors

Remarkable structural homologies between the main proteases of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and the avian infectious bronchitis virus (IBV) were revealed by comparative amino-acid sequence and 3D structural alignment. Assessing whether reported IBV 3CLPro inhibitors could also interact with SARS-CoV-2 has been undertaken in silico using a PubChem BioAssay database of 388 compounds active on the avian infectious bronchitis virus 3C-like protease. Docking studies of this database on the SARS-CoV-2 protease resulted in the identification of four covalent inhibitors targeting the catalytic cysteine residue and five non-covalent inhibitors for which the binding was further investigated by molecular dynamics (MD) simulations. Predictive ADMET calculations on the nine compounds suggest promising pharmacokinetic properties.

City-scale model for COVID-19 epidemiology with mobility and social activities represented by a set of hidden Markov models

BACKGROUND AND OBJECTIVE

SARS-CoV-2 emerged by the end of 2019 and became a global pandemic due to its rapid spread. Various outbreaks of the disease in different parts of the world have been studied, and epidemiological analyses of these outbreaks have been useful for developing models with the aim of tracking and predicting the spread of epidemics. In this paper, an agent-based model that predicts the local daily evolution of the number of people hospitalized in intensive care due to COVID-19 is presented.

METHODS

An agent-based model has been developed, taking into consideration the most relevant characteristics of the geography and climate of a mid-size city, its population and pathology statistics, and its social customs and mobility, including the state of public transportation. In addition to these inputs, the different phases of isolation and social distancing are also taken into account. By means of a set of hidden Markov models, the system captures and reproduces virus transmission associated with the stochastic nature of people's mobility and activities in the city. The spread of the virus in the host is also simulated by following the stages of the disease and by considering the existence of comorbidities and the proportion of asymptomatic carriers.

RESULTS

As a case study, the model was applied to Paraná city (Entre Ríos, Argentina) in the second half of 2020. The model adequately predicts the daily evolution of people hospitalized in intensive care due to COVID-19. This adequacy is reflected by the fact that the prediction of the model (including its dispersion), as with the data reported in the field, never exceeded 90% of the capacity of beds installed in the city. In addition, other epidemiological variables of interest, with discrimination by age range, were also adequately reproduced, such as the number of deaths, reported cases, and asymptomatic individuals.

CONCLUSIONS

The model can be used to predict the most likely evolution of the number of cases and hospital bed occupancy in the short term. By adjusting the model to match the data on hospitalizations in intensive care units and deaths due to COVID-19, it is possible to analyze the impact of isolation and social distancing measures on the disease spread dynamics. In addition, it allows for simulating combinations of characteristics that would lead to a potential collapse in the health system due to lack of infrastructure as well as predicting the impact of social events or increases in people's mobility.

Comprehensive analysis of lncRNA-mRNA co-expression networks in HPV-driven cervical cancer reveals the pivotal function of LINC00511-PGK1 in tumorigenesis

BACKGROUND

Mounting evidence suggests that noncoding RNAs (lncRNAs) were involved in various human cancers. However, the role of these lncRNAs in HPV-driven cervical cancer (CC) has not been extensively studied. Considering that HR-HPV infections contribute to cervical carcinogenesis by regulating the expression of lncRNAs, miRNAs and mRNAs, we aim to systematically analyze lncRNAs and mRNAs expression profile to identify novel lncRNAs-mRNAs co-expression networks and explore their potential impact on tumorigenesis in HPV-driven CC.

METHODS

LncRNA/mRNA microarray technology was utilized to identify the differentially expressed lncRNAs (DElncRNAs) and mRNAs (DEmRNAs) in HPV-16 and HPV-18 cervical carcinogenesis compared to normal cervical tissues. Venn diagram and weighted gene co-expression network analysis (WGCNA) were used to identify the hub DElncRNAs/DEmRNAs that were both significantly correlated with HPV-16 and HPV-18 CC patients. LncRNA-mRNA correlation analysis and functional enrichment pathway analysis were performed on these key DElncRNAs/DEmRNAs in HPV-16 and HPV-18 CC patients to explore their mutual mechanism in HPV-driven CC. A lncRNA-mRNA co-expression score (CES) model was established and validated by using the Cox regression method. Afterward, the clinicopathological characteristics were analyzed between CES-high and CES-low groups. In vitro, functional experiments were performed to evaluate the role of LINC00511 and PGK1 in cell proliferation, migration and invasion in CC cells. To understand whether LINC00511 play as an oncogenic role partially via modulating the expression of PGK1, rescue assays were used.

RESULTS

We identified 81 lncRNAs and 211 mRNAs that were commonly differentially expressed in HPV-16 and HPV-18 CC tissues compared to normal tissues. The results of lncRNA-mRNA correlation analysis and functional enrichment pathway analysis showed that the LINC00511-PGK1 co-expression network may make an important contribution to HPV-mediated tumorigenesis and be closely associated with metabolism-related mechanisms. Combined with clinical survival data, the prognostic lncRNA-mRNA co-expression score (CES) model based on LINC00511 and PGK1 could precisely predict patients' overall survival (OS). CES-high patients had a worse prognosis than CES-low patients and the enriched pathways and potential targets of applicable drugs were explored in CES-high patients. In vitro experiments confirmed the oncogenic functions of LINC00511 and PGK1 in the progression of CC, and revealed that LINC00511 functions in an oncogenic role in CC cells partially via modulating the expression of PGK1.

CONCLUSIONS

Together, these data identify co-expression modules that provide valuable information to understand the pathogenesis of HPV-mediated tumorigenesis, which highlights the pivotal function of the LINC00511-PGK1 co-expression network in cervical carcinogenesis. Furthermore, our CES model has a reliable predicting ability that could stratify CC patients into low- and high-risk groups of poor survival. This study provides a bioinformatics method to screen prognostic biomarkers which leads to lncRNA-mRNA co-expression network identification and construction for patients' survival prediction and potential drug applications in other cancers.

A novel five-gene metabolism-related risk signature for predicting prognosis and immune infiltration in endometrial cancer: A TCGA data mining

BACKGROUND

Metabolism dysfunction can affect the biological behavior of tumor cells and result in carcinogenesis and the development of various cancers. However, few thoughtful studies focus on the predictive value and efficacy of immunotherapy of metabolism-related gene signatures in endometrial cancer (EC). This research aims to construct a predictive metabolism-related gene signature in EC with prognostic and therapeutic implications.

METHODS

We downloaded the RNA profile and clinical data of 503 EC patients and screened out different expressions of metabolism-related genes with prognosis influence of EC from The Cancer Genome Atlas (TCGA) database. We first established a metabolism-related genes model using univariate and multivariate Cox regression and Lasso regression analysis. To internally validate the predictive model, 503 samples (entire set) were randomly assigned into the test set and the train set. Then, we applied the receiver operating characteristic (ROC) curve to confirm our previous predictive model and depicted a nomogram integrating the risk score and the clinicopathological feature. We employed a gene set enrichment analysis (GSEA) to explore the biological processes and pathways of the model. Afterward, we used ESTIMATE to evaluate the TME. Also, we adopted CIBERSORT and ssGSEA to estimate the fraction of immune infiltrating cells and immune function. At last, we investigated the relationship between the predictive model and immune checkpoint genes.

RESULTS

We first constructed a predictive model based on five metabolism-related genes (INPP5K, PLPP2, MBOAT2, DDC, and ITPKA). This model showed the ability to predict EC patients' prognosis accurately and performed well in the train set, test set, and entire set. Then we confirmed the predictive signature was a novel independent prognostic factor in EC patients. In addition, we drew and validated a nomogram to precisely predict the survival rate of EC patients at 1-, 3-, and 5-years (ROC_1-year = 0.714, ROC_3-year = 0.750, ROC_5-year = 0.767). Furthermore, GSEA unveiled that the cell cycle, certain malignant tumors, and cell metabolism were the main biological functions enriched in this identified model. We found the five metabolism-related genes signature was associated with the immune infiltrating cells and immune functions. Most importantly, it was linked with specific immune checkpoints (PD-1, CTLA4, and CD40) that could predict immunotherapy's clinical response.

CONCLUSION

The metabolism-related genes signature (INPP5K, PLPP2, MBOAT2, DDC, and ITPKA) is a valuable index for predicting the survival outcomes and efficacy of immunotherapy for EC in clinical settings.

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.

Computational approach investigation bioactive molecules from Saussurea Costus plant as SARS-CoV-2 main protease inhibitors using reverse docking, molecular dynamics simulation, and pharmacokinetic ADMET parameters

SARS-COV-2 virus causes (COVID-19) disease; it has become a global pandemic since 2019 and has negatively affected all aspects of human life. Scientists have made great efforts to find a reliable cure, vaccine, or treatment for this emerging disease. Efforts have been directed towards using medicinal plants as alternative medicines, as the active chemical compounds in them have been discovered as potential antiviral or anti-inflammatory agents. In this research, the potential of Saussurea costus (S. Costus) or QUST Al Hindi chemical consistent as potential antiviral agents was investigated by using computational methods such as Reverse Docking, ADMET, and Molecular Dynamics with different proteases COVID-19 such as PDB: 2GZ9; 6LU7; 7AOL, 6Y2E, 6Y84. The results of Reverse Docking the complex between 6LU7 proteases and Cynaropicrin compound being the best complex, as the same result, is achieved by molecular dynamics. Also, the toxicity testing result from ADMET method proved that the complex is the least toxic and the safest possible drug. In addition, 6LU7-Cynaropicrin complex obeyed Lipinski rule; it formed ≤5 H-bond donors and ≤10 H bond acceptors, MW < 500 Daltons, and octanol/water partition coefficient <5.

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 Putative long-range RNA-RNA interaction between ORF8 and Spike of SARS-CoV-2

SARS-CoV-2 has affected people worldwide as the causative agent of COVID-19. The virus is related to the highly lethal SARS-CoV-1 responsible for the 2002-2003 SARS outbreak in Asia. Research is ongoing to understand why both viruses have different spreading capacities and mortality rates. Like other beta coronaviruses, RNA-RNA interactions occur between different parts of the viral genomic RNA, resulting in discontinuous transcription and production of various sub-genomic RNAs. These sub-genomic RNAs are then translated into other viral proteins. In this work, we performed a comparative analysis for novel long-range RNA-RNA interactions that may involve the Spike region. Comparing in-silico fragment-based predictions between reference sequences of SARS-CoV-1 and SARS-CoV-2 revealed several predictions amongst which a thermodynamically stable long-range RNA-RNA interaction between (23660-23703 Spike) and (28025-28060 ORF8) unique to SARS-CoV-2 was observed. The patterns of sequence variation using data gathered worldwide further supported the predicted stability of the sub-interacting region (23679-23690 Spike) and (28031-28042 ORF8). Such RNA-RNA interactions can potentially impact viral life cycle including sub-genomic RNA production rates.

GC-MS profiling of Bauhinia variegata major phytoconstituents with computational identification of potential lead inhibitors of SARS-CoV-2 Mpro

Severe acute respiratory syndrome coronavirus 2 was originally identified in Wuhan city of China in December 2019 and it spread rapidly throughout the globe, causing a threat to human life. Since targeted therapies are deficient, scientists all over the world have an opportunity to develop novel drug therapies to combat COVID-19. After the declaration of a global medical emergency, it was established that the Food and Drug Administration (FDA) could permit the use of emergency testing, treatments, and vaccines to decrease suffering, and loss of life, and restore the nation's health and security. The FDA has approved the use of remdesivir and its analogs as an antiviral medication, to treat COVID-19. The primary protease of SARS-CoV-2, which has the potential to regulate coronavirus proliferation, has been a viable target for the discovery of medicines against SARS-CoV-2. The present research deals with the in silico technique to screen phytocompounds from a traditional medicinal plant, Bauhinia variegata for potential inhibitors of the SARS-CoV-2 main protease. Dried leaves of the plant B. variegata were used to prepare aqueous and methanol extract and the constituents were analyzed using the GC-MS technique. A total of 57 compounds were retrieved from the aqueous and methanol extract analysis. Among these, three lead compounds (2,5 dimethyl 1-H Pyrrole, 2,3 diphenyl cyclopropyl methyl phenyl sulphoxide, and Benzonitrile m phenethyl) were shown to have the highest binding affinity (−5.719 to −5.580 kcal/mol) towards SARS-CoV-2 M^pro. The post MD simulation results also revealed the favorable confirmation and stability of the selected lead compounds with M^pro as per trajectory analysis. The Prime MM/GBSA binding free energy supports this finding, the top lead compound 2,3 diphenyl cyclopropyl methyl phenyl sulphoxide showed high binding free energy (−64.377 ± 5.24 kcal/mol) towards M^pro which reflects the binding stability of the molecule with M^pro. The binding free energy of the complexes was strongly influenced by His, Gln, and Glu residues. All of the molecules chosen are found to have strong pharmacokinetic characteristics and show drug-likeness properties. The lead compounds present acute toxicity (LD50) values ranging from 670 mg/kg to 2500 mg/kg; with toxicity classifications of 4 and 5 classes. Thus, these compounds could behave as probable lead candidates for treatment against SARS-CoV-2. However further in vitro and in vivo studies are required for the development of medication against SARS-CoV-2.

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

In silico exploration of disulfide derivatives of Ferula foetida oleo-gum (Covexir®) as promising therapeutics against SARS-CoV-2

Although vaccines have been significantly successful against coronavirus, due to the high rate of the Omicron variant spread many researchers are trying to find efficient drugs against COVID-19. Herein, we conducted a computational study to investigate the binding mechanism of four potential inhibitors (including disulfide derivatives isolated from Ferula foetida) to SARS-CoV-2 main protease. Our findings revealed that the disulfides mainly interacted with HIS41, MET49, CYS145, HIS64, MET165, and GLN189 residues of SARS-CoV-2 main protease. The binding free energy decomposition results also showed that the van der Waals (vdW) energy plays the main role in the interaction of HIS41, MET49, CYS145, HIS64, MET165, and GLN189 residues with the inhibitors. Furthermore, it is found that the Z-isomer derivatives have a stronger interaction with SARS-CoV-2, and the strongest interaction belongs to the (Z)-1-(1-(methylthio)propyl)-2-(prop-1-enyl)disulfane (ΔG = -18.672 kcal/mol). The quantum mechanical calculations demonstrated that the second-order perturbation stabilization energy and the electron density values for MET49-ligand interactions are higher than the other residue-ligand complexes. This finding confirms the stronger interaction of this residue with the ligands.

Most frequently harboured missense variants of hACE2 across different populations exhibit varying patterns of binding interaction with spike glycoproteins of emerging SARS-CoV-2 of different lineages

Since the emergence of SARS-CoV-2 at Wuhan in the Hubei province of China in 2019, the virus has accumulated various mutations, giving rise to many variants. According to the combinations of mutations acquired, these variants are classified into lineages and greatly differ in infectivity and transmissibility. In 2021 alone, a variant of interest (VoI) Mu (B.1.621), as well as, variants of concern (VoC) Delta (B.1.617.2) and Omicron (BA.1, BA.2) and later in 2022, BA.4, BA.5, and BA.2.12.1 have emerged. Since then, the world has seen prominent surges in the rate of infection during short periods of time. However, not all populations have suffered equally, which suggests a possible role of host genetic factors. Here, we investigated the strength of binding of the spike glycoprotein receptor-binding domain (RBD) of the SARS-CoV-2 variants: Mu, Delta, Delta Plus (AY.1), Omicron sub-variants BA.1, BA.2, BA.4, BA.5, and BA.2.12.1 with the human angiotensin-converting enzyme 2 (hACE2) missense variants prevalent in major populations. In this purpose, molecular docking analysis, as well as, molecular dynamics simulation was performed of the above-mentioned SARS-CoV-2 RBD variants with the hACE2 containing the single amino acid substitutions prevalent in African (E37K), Latin American (F40L), non-Finnish European (D355 N), and South Asian (P84T) populations, in order to predict the effects of the lineage-defining mutations of the viral variants on receptor binding. The effects of these mutations on protein stability were also explored. The protein-protein docking and molecular dynamics simulation analyses have revealed variable strength of attachment and exhibited altered interactions in the case of different hACE2-RBD complexes. In vitro studies are warranted to confirm these findings which may enable early prediction regarding the risk of transmissibility of newly emerging variants across different populations in the future.

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

De Novo design of potential inhibitors against SARS-CoV-2 Mpro

The impact of the ravages of COVID-19 on people's lives is obvious, and the development of novel potential inhibitors against SARS-CoV-2 main protease (Mpro), which has been validated as a potential target for drug design, is urgently needed. This study developed a model named MproI-GEN, which can be used for the de novo design of potential Mpro inhibitors (MproIs) based on deep learning. The model was mainly composed of long-short term memory modules, and the last layer was re-trained with transfer learning. The validity (0.9248), novelty (0.9668), and uniqueness (0.0652) of the designed potential MproI library (PMproIL) were evaluated, and the results showed that MproI-GEN could be used to design structurally novel and reasonable molecules. Additionally, PMproIL was filtered based on machine learning models and molecular docking. After filtering, the potential MproIs were verified with molecular dynamics simulations to evaluate the binding stability levels of these MproIs and SARS-CoV-2 Mpro, thereby illustrating the inhibitory effects of the potential MproIs against Mpro. Two potential MproIs were proposed in this study. This study provides not only new possibilities for the development of COVID-19 drugs but also a complete pipeline for the discovery of novel lead compounds.

MutCov: A pipeline for evaluating the effect of mutations in spike protein on infectivity and antigenicity of SARS-CoV-2

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), causing an outbreak of coronavirus disease 2019 (COVID-19), is a major threat to public health worldwide. Previous studies have shown that the spike protein of SARS-CoV-2 determines viral infectivity and major antigenicity. However, the spike protein has been undergoing various mutations, which bring a great challenge to the prevention and treatment of COVID-19. Here we present the MutCov, a pipeline for evaluating the effect of mutations in spike protein on infectivity and antigenicity of SARS-CoV-2 by calculating the binding free energy between spike protein and angiotensin-converting enzyme 2 (ACE2) or neutralizing monoclonal antibody (mAb). The predicted infectivity and antigenicity were highly consistent with biologically experimental results, and demonstrated that the MutCov achieved good prediction performance. In conclusion, the MutCov is of high importance for systematically evaluating the effect of novel mutations and improving the prevention and treatment of COVID-19. The source code and installation instruction of MutCov are freely available at http://jianglab.org.cn/MutCov.

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.

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.

Computational studies on the interaction of SARS-CoV-2 Omicron SGp RBD with human receptor ACE2, limonin and glycyrrhizic acid

On November 24, 2021, the SARS-CoV-2 Omicron variant (B.1.1.529) was first identified in South Africa. The World Health Organization (WHO) declared the Omicron as a variant of concern (VoC) because of the unexpected and large numbers of mutations occurred in the genome, higher viral transmission and immune evasions. The present study was performed to explore the interactions of SARS-CoV-2 spike glycoprotein receptor-binding domain (SGp RBD) of the three variants (Omicron, Delta, and WT) with the receptor hACE2. The structural changes occurred in Omicron due to the mutations at key positions improved the ability to mediate SARS-CoV-2 viral infection compared to other VoCs. The phytochemicals limonin and glycyrrhizic acid were docked with the SGp RBD of the variants WT, Delta and Omicron. The computed dock score revealed that limonin and glycyrrhizic acid binds effectively at the SGp RBD of all three variants, and showed almost similar binding affinity at the binding interface of ACE2. Therefore, despite the multiple mutations occurred in Omicron and its viral transmission is comparatively high, the computed binding affinity of the phytochemicals limonin and glycyrrhizic acid supported that the traditional medicines can be useful in formulating adjuvant therapies to fight against the SARS-CoV-2 Omicron.

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.

Multi-omics integration in the age of million single-cell data

An explosion in single-cell technologies has revealed a previously underappreciated heterogeneity of cell types and novel cell-state associations with sex, disease, development and other processes. Starting with transcriptome analyses, single-cell techniques have extended to multi-omics approaches and now enable the simultaneous measurement of data modalities and spatial cellular context. Data are now available for millions of cells, for whole-genome measurements and for multiple modalities. Although analyses of such multimodal datasets have the potential to provide new insights into biological processes that cannot be inferred with a single mode of assay, the integration of very large, complex, multimodal data into biological models and mechanisms represents a considerable challenge. An understanding of the principles of data integration and visualization methods is required to determine what methods are best applied to a particular single-cell dataset. Each class of method has advantages and pitfalls in terms of its ability to achieve various biological goals, including cell-type classification, regulatory network modelling and biological process inference. In choosing a data integration strategy, consideration must be given to whether the multi-omics data are matched (that is, measured on the same cell) or unmatched (that is, measured on different cells) and, more importantly, the overall modelling and visualization goals of the integrated analysis.