COVID-19 is a global health emergency that causes serious concerns. A global effort is underway to identify drugs for the treatment of COVID-19. One possible solution to the present problem is to develop drugs that can inhibit SARS-CoV-2 main protease (M^pro), a coronavirus protein that been considered as one among many drug targets. In this work, lactoferrin from Bos taurus L. was in silico hydrolyzed. The bioactivity, water solubility, and ADMET properties of the generated peptides were predicted using various online tools. The molecular interactions between M^pro and the peptides were studied using molecular docking and molecular dynamic simulation. The results demonstrated that peptide GSRY was predicted to have better physicochemical properties, and the value of ‘-C DOCKER interaction energy’ between peptide GSRY and M^pro was 80.8505 kcal/mol. The interaction between the peptide GSRY and the native ligand N3 co-crystallized with M^pro had overlapped amino acids, i.e., HIS163, GlY143, GLU166, GLN189 and MET165. Molecular dynamic simulation revealed that M^pro/GSRY complexes were stable. Collectively, the peptide GSRY may be a potential candidate drug against M^pro of SARS-CoV-2.

Adenosine-to-inosine (A-to-I) RNA editing, constituting nearly 90% of all RNA editing events in humans, has been reported to contribute to the tumorigenesis in diverse cancers. However, the comprehensive map for functional A-to-I RNA editing events in cancers is still insufficient. To fill this gap, we systematically and intensively analyzed multiple tumorigenic mechanisms of A-to-I RNA editing events in samples across 33 cancer types from The Cancer Genome Atlas. For individual candidate among ∼ 1,500,000 quantified RNA editing events, we performed diverse types of downstream functional annotations. Finally, we identified 24,236 potentially functional A-to-I RNA editing events, including the cases in APOL1, IGFBP3, GRIA2, BLCAP, and miR-589-3p. These events might play crucial roles in the scenarios of tumorigenesis, due to their tumor-related editing frequencies or probable effects on altered expression profiles, protein functions, splicing patterns, and microRNA regulations of tumor genes. Our functional A-to-I RNA editing events (https://ccsm.uth.edu/CAeditome/) will help better understand the cancer pathology from the A-to-I RNA editing aspect.

Parkinson’s disease (PD) is characterized by dopaminergic neurodegeneration and an abnormal accumulation of α-synuclein aggregates. A number of genetic factors have been shown to increase the risk of PD. Exploring the underlying molecular mechanisms that mediate PD’s transcriptomic diversity can help us understand neurodegenerative pathogenesis. In this study, we identified 9897 A-to-I RNA editing events associated with 6286 genes across 372 PD patients. Of them, 72 RNA editing events altered miRNA binding sites and this may directly affect miRNA regulations of their host genes. However, RNA editing effects on the miRNA regulation of genes are more complex. They can (1) abolish existing miRNA binding sites, which allows miRNAs to regulate other genes; (2) create new miRNA binding sites that may sequester miRNAs from regulating other genes; or (3) occur in the miRNA seed regions and change their targets. The first two processes are also referred to as miRNA competitive binding. In our study, we found 8 RNA editing events that may alter the expression of 1146 other genes via miRNA competition. We also found one RNA editing event that modified a miRNA seed region, which was predicted to disturb the regulation of four genes. Considering the PD-related functions of the affected genes, 25 A-to-I RNA editing biomarkers for PD are proposed, including the 3 editing events in the EIF2AK2, APOL6, and miR-4477b seed regions. These biomarkers may alter the miRNA regulation of 133 PD-related genes. All these analyses reveal the potential mechanisms and regulations of RNA editing in PD pathogenesis.

Objective: Preeclampsia (PE) is a serious condition in pregnant women and hence an important topic in obstetrics. The current research aimed to recognize the potential and significant immune-related diagnostic biomarkers for PE.Methods: From the Gene Expression Omnibus (GEO) data sets, three public gene expression profiles (GSE24129, GSE54618, and GSE60438) from the placental samples of PE and normotensive pregnancy were downloaded. Differentially expressed genes (DEGs) were selected and determined among 73 PE and 85 normotensive control pregnancy samples. The DEGs were used for Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), Disease Ontology (DO) enrichment analysis, and Gene Set Enrichment Analysis (GSEA). The candidate biomarkers were identified by the least absolute shrinkage and selection operator (LASSO) and support vector machine recursive feature elimination (SVM-RFE) analysis. The receiver operating characteristic curve (ROC) was applied to evaluate diagnostic ability. For further confirmation, the expression levels and diagnostic value of biomarkers in PE were verified in the GSE75010 data set (80 PE and 77 controls) and validated by qRT-RCR, Western blot, and immunohistochemistry (IHC). The CIBERSORT algorithm was used to calculate the compositional patterns of 22 types of immune cells in PE.Results: In total, 15 DEGs were recognized. The GO and KEGG analyses revealed that the DEGs were enriched in the steroid metabolic process, receptor ligand activity, GnRH secretion, and neuroactive ligand−receptor interaction. The recognized DEGs were primarily implicated in cell-type benign neoplasm, kidney failure, infertility, and PE. Gene sets related to hormone activity, glycosylation, multicellular organism process, and response to BMP were activated in PE. The LEP gene was distinguished as a diagnostic biomarker of PE (AUC = 0.712) and further certified in the GSE75010 data set (AUC = 0.850). The high expression of LEP was associated with PE in clinical samples. In addition, the analysis of the immune microenvironment showed that gamma delta T cells, memory B cells, M0 macrophages, and regulatory T cells were positively correlated with LEP expression (P < 0.05).Conclusion:LEP expression can be considered to be a diagnostic biomarker of PE and can offer a novel perspective for future studies regarding the occurrence and molecular mechanisms of PE.

The metastatic non-small cell lung cancer (NSCLC) is one of the cancers with high incidence, poor survival, and limited treatment. Epithelial-mesenchymal transition (EMT) is the first step by which an early tumor converts to an invasive one. Studying the underlying mechanisms of EMT can help the understanding of cancer metastasis and improve the treatment. In this study, 1013 NSCLC patients and 123 NSCLC cell lines are deeply analyzed for the potential roles of alternative polyadenylation (APA) in the EMT process. A trend of shorter 3′-UTRs (three prime untranslated region) is discovered in the mesenchymal samples. The identification of EMT-related APA events highlights the proximal poly(A) selection of CARM1. It is a pathological biomarker of mesenchymal tumor and cancer metastasis through losing miRNA binding to upregulate the EMT inducer of CARM1 and releasing miRNAs to downregulate the EMT inhibitor of RBM47. The crucial role of this APA event in EMT also guides its effect on drug responses. The patients with shorter 3′-UTR of CARM1 are more benefit from chemotherapy drugs, especially cisplatin. A stratification of NSCLC patients based on this APA event is useful for chemotherapy design in future clinics.

A-to-I RNA editing diversifies human transcriptome to confer its functional effects on the downstream genes or regulations, potentially involving in neurodegenerative pathogenesis. Its variabilities are attributed to multiple regulators, including the key factor of genetic variants. To comprehensively investigate the potentials of neurodegenerative disease-susceptibility variants from the view of A-to-I RNA editing, we analyzed matched genetic and transcriptomic data of 1596 samples across nine brain tissues and whole blood from two large consortiums, Accelerating Medicines Partnership-Alzheimer’s Disease and Parkinson’s Progression Markers Initiative. The large-scale and genome-wide identification of 95 198 RNA editing quantitative trait loci revealed the preferred genetic effects on adjacent editing events. Furthermore, to explore the underlying mechanisms of the genetic controls of A-to-I RNA editing, several top RNA-binding proteins were pointed out, such as EIF4A3, U2AF2, NOP58, FBL, NOP56 and DHX9, since their regulations on multiple RNA-editing events were probably interfered by these genetic variants. Moreover, these variants may also contribute to the variability of other molecular phenotypes associated with RNA editing, including the functions of 3 proteins, expressions of 277 genes and splicing of 449 events. All the analyses results shown in NeuroEdQTL (https://relab.xidian.edu.cn/NeuroEdQTL/) constituted a unique resource for the understanding of neurodegenerative pathogenesis from genotypes to phenotypes related to A-to-I RNA editing.