Essential Molecular Function Research Articles

M6A-Atlas v2.0: updated resources for unraveling the N6-methyladenosine (m6A) epitranscriptome among multiple species.

N 6-Methyladenosine (m6A) is one of the most abundant internal chemical modifications on eukaryote mRNA and is involved in numerous essential molecular functions and biological processes. To facilitate the study of this important post-transcriptional modification, we present here m6A-Atlas v2.0, an updated version of m6A-Atlas. It was expanded to include a total of 797 091 reliable m6A sites from 13 high-resolution technologies and two single-cell m6A profiles. Additionally, three methods (exomePeaks2, MACS2 and TRESS) were used to identify >16 million m6A enrichment peaks from 2712 MeRIP-seq experiments covering 651 conditions in 42 species. Quality control results of MeRIP-seq samples were also provided to help users to select reliable peaks. We also estimated the condition-specific quantitative m6A profiles (i.e. differential methylation) under 172 experimental conditions for 19 species. Further, to provide insights into potential functional circuitry, the m6A epitranscriptomics were annotated with various genomic features, interactions with RNA-binding proteins and microRNA, potentially linked splicing events and single nucleotide polymorphisms. The collected m6A sites and their functional annotations can be freely queried and downloaded via a user-friendly graphical interface at: http://rnamd.org/m6a.

Unanticipated Large-Scale Deletion in Fusarium graminearum Genome Using CRISPR/Cas9 and Its Impact on Growth and Virulence

Fusarium graminearum, a filamentous fungus, and causal agent of Fusarium head blight (FHB) in wheat and other cereals, leads to significant economic losses globally. This study aimed to investigate the roles of specific genes in F. graminearum virulence using CRISPR/Cas9-mediated gene deletions. Illumina sequencing was used to characterize the genomic changes due to editing. Unexpectedly, a large-scale deletion of 525,223 base pairs on chromosome 2, comprising over 222 genes, occurred in two isolates. Many of the deleted genes were predicted to be involved in essential molecular functions, such as oxidoreductase activity, transmembrane transporter activity, hydrolase activity, as well as biological processes, such as carbohydrate metabolism and transmembrane transport. Despite the substantial loss of genetic material, the mutant isolate exhibited normal growth rates and virulence on wheat under most conditions. However, growth rates were significantly reduced under high temperatures and on some media. Additionally, wheat inoculation assays using clip dipping, seed inoculation, and head point inoculation methods were performed. No significant differences in virulence were observed, suggesting that these genes were not involved in infection or alternative compensatory pathways, and allow the fungi to maintain pathogenicity despite the extensive genomic deletion.

Abstract A030: Deciphering the micro-RNA landscape of AKT isoforms and its clinical impact in patients with lung adenocarcinoma

Abstract Over the last 25 years, AKT kinase has been thoroughly investigated and characterized as a crucial regulator of cellular homeostasis, affecting a plethora of essential molecular functions. The three identified isoforms of this kinase, namely AKT1, AKT2 and AKT3, often display opposing and antagonizing cellular functions. One of these important cellular functions is the expression and activation of micro-RNAs (miRNAs), an important RNA process that epigenetically regulates the human transcriptome. In the context of human cancer and based on several reports, AKT kinase has been associated with aberrant expression of oncogenic miRNAs, that in turn hijack oncogenic pathways, leading to tumor progression, metastasis and resistance to therapy, especially in lung adenocarcinoma (LUAD). Nevertheless, a comprehensive analysis of the impact of each AKT isoform on the regulation and enrichment of a specific miRNA signature in lung adenocarcinoma has yet to be performed. Here, we were willing to elucidate the miRNA epigenetic landscape and downstream targets, driven by each AKT isoform, and its clinical importance in patients with lung adenocarcinoma. To this extent, we utilized RNA-Seq and Small RNA-Seq data from The Cancer Genome Atlas (TCGA-LUAD) to identify the miRNA signatures related with the expression of each AKT isoform in LUAD. Gene Set Enrichment Analysis from RNA-Seq and miRNA levels analysis from Small RNA-Seq, identified a unique miRNA signature for each AKT isoform, demonstrating their largely opposing effect in miRNA profiling. In parallel analysis, in silico target analysis of the top 50 regulated miRNAs, demonstrated another set of unique targets for each isoform-dependent miRNA signature, which is involved in metabolism, DNA damage, RNA processing, cell survival, immune response and metastasis. More importantly, clinical and survival analysis revealed that the top 50 Differentially regulated miRNAs, ranked by q-value, for each isoform correlated with poor progression free survival in patients with LUAD. Collectively, our atlas of AKT isoform-driven miRNA landscapes implicates AKT as a global effector of miRNA landscape in lung adenocarcinoma. Due to the lack of reliable isoform-specific AKT inhibitors, the identification of the oncogenic miRNA networks regulated by each isoform, acts as an important transnational resource for the design and testing of inhibitors that bypass AKT inhibition toxicity, sparing its impactful physiological functions. Citation Format: George Laliotis, Pauline Lin. Deciphering the micro-RNA landscape of AKT isoforms and its clinical impact in patients with lung adenocarcinoma. [abstract]. In: Proceedings of the AACR Special Conference: Cancer Epigenomics; 2022 Oct 6-8; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2022;82(23 Suppl_2):Abstract nr A030.

DMFpred: Predicting protein disorder molecular functions based on protein cubic language model.

Intrinsically disordered proteins and regions (IDP/IDRs) are widespread in living organisms and perform various essential molecular functions. These functions are summarized as six general categories, including entropic chain, assembler, scavenger, effector, display site, and chaperone. The alteration of IDP functions is responsible for many human diseases. Therefore, identifying the function of disordered proteins is helpful for the studies of drug target discovery and rational drug design. Experimental identification of the molecular functions of IDP in the wet lab is an expensive and laborious procedure that is not applicable on a large scale. Some computational methods have been proposed and mainly focus on predicting the entropic chain function of IDRs, while the computational predictive methods for the remaining five important categories of disordered molecular functions are desired. Motivated by the growing numbers of experimental annotated functional sequences and the need to expand the coverage of disordered protein function predictors, we proposed DMFpred for disordered molecular functions prediction, covering disordered assembler, scavenger, effector, display site and chaperone. DMFpred employs the Protein Cubic Language Model (PCLM), which incorporates three protein language models for characterizing sequences, structural and functional features of proteins, and attention-based alignment for understanding the relationship among three captured features and generating a joint representation of proteins. The PCLM was pre-trained with large-scaled IDR sequences and fine-tuned with functional annotation sequences for molecular function prediction. The predictive performance evaluation on five categories of functional and multi-functional residues suggested that DMFpred provides high-quality predictions. The web-server of DMFpred can be freely accessed from http://bliulab.net/DMFpred/.

Melatonin: Regulation of Biomolecular Condensates in Neurodegenerative Disorders.

Biomolecular condensates are membraneless organelles (MLOs) that form dynamic, chemically distinct subcellular compartments organizing macromolecules such as proteins, RNA, and DNA in unicellular prokaryotic bacteria and complex eukaryotic cells. Separated from surrounding environments, MLOs in the nucleoplasm, cytoplasm, and mitochondria assemble by liquid–liquid phase separation (LLPS) into transient, non-static, liquid-like droplets that regulate essential molecular functions. LLPS is primarily controlled by post-translational modifications (PTMs) that fine-tune the balance between attractive and repulsive charge states and/or binding motifs of proteins. Aberrant phase separation due to dysregulated membrane lipid rafts and/or PTMs, as well as the absence of adequate hydrotropic small molecules such as ATP, or the presence of specific RNA proteins can cause pathological protein aggregation in neurodegenerative disorders. Melatonin may exert a dominant influence over phase separation in biomolecular condensates by optimizing membrane and MLO interdependent reactions through stabilizing lipid raft domains, reducing line tension, and maintaining negative membrane curvature and fluidity. As a potent antioxidant, melatonin protects cardiolipin and other membrane lipids from peroxidation cascades, supporting protein trafficking, signaling, ion channel activities, and ATPase functionality during condensate coacervation or dissolution. Melatonin may even control condensate LLPS through PTM and balance mRNA- and RNA-binding protein composition by regulating N6-methyladenosine (m6A) modifications. There is currently a lack of pharmaceuticals targeting neurodegenerative disorders via the regulation of phase separation. The potential of melatonin in the modulation of biomolecular condensate in the attenuation of aberrant condensate aggregation in neurodegenerative disorders is discussed in this review.

PDCD2 functions as an evolutionarily conserved chaperone dedicated for the 40S ribosomal protein uS5 (RPS2).

PDCD2 is an evolutionarily conserved protein with previously characterized homologs in Drosophila (zfrp8) and budding yeast (Tsr4). Although mammalian PDCD2 is essential for cell proliferation and embryonic development, the function of PDCD2 that underlies its fundamental cellular role has remained unclear. Here, we used quantitative proteomics approaches to define the protein-protein interaction network of human PDCD2. Our data revealed that PDCD2 specifically interacts with the 40S ribosomal protein uS5 (RPS2) and that the PDCD2-uS5 complex is assembled co-translationally. Loss of PDCD2 expression leads to defects in the synthesis of the small ribosomal subunit that phenocopy a uS5 deficiency. Notably, we show that PDCD2 is important for the accumulation of soluble uS5 protein as well as its incorporation into 40S ribosomal subunit. Our findings support that the essential molecular function of PDCD2 is to act as a dedicated ribosomal protein chaperone that recognizes uS5 co-translationally in the cytoplasm and accompanies uS5 to ribosome assembly sites in the nucleus. As most dedicated ribosomal protein chaperones have been identified in yeast, our study reveals that similar mechanisms exist in human cells to assist ribosomal proteins coordinate their folding, nuclear import and assembly in pre-ribosomal particles.

Protein Network Analysis of the Serum and Their Functional Implication in Patients Subjected to Traumatic Brain Injury.

Traumatic brain injury (TBI) often leads to severe neurobehavioral impairment, but the underlying molecular mechanism remains to be elucidated. Here, we collected the sera from 23 patients (aged from 19 to 81 years old, third day after TBI as TBI-third group) subjected to TBI from The First Hospital of Kunming City, and the sera from 22 healthy donors (aged from 18 to 81 years old and as control group). Then, three samples from TBI-third group and three samples from control group were subjected to the protein microarray detection, and bioinformatics analysis. Then, enzyme-linked immunosorbent assay (ELISA) was used to verify significantly altered protein levels. Results showed that, when compared with the control group, all significantly differentially expressed proteins [DEPs, P < 0.05, FDR < 0.05, fold change (FC) > 2] contained 172 molecules in the TBI-third group, in which 65 proteins were upregulated, while 107 proteins were downregulated. The biological processes of these DEPs, mostly happened in the extracellular region and the extracellular region parts, are mainly involved in the regulation of cellular process, signaling and signal transduction, cell communication, response to stimuli, the immune system process and multicellular organismal development. Moreover, the essential molecular functions of them are cytokine activity, growth factor activity and morphogen activity. Additionally, the most significant pathways are enriched in cytokine–cytokine receptor interaction and PI3K-Akt signaling pathways among downregulated proteins, and pathways in cancer and cytokine–cytokine receptor interaction among upregulated proteins. Of these, we focused on the NGF, NT-3, IGF-2, HGF, NPY, CRP, MMP-9, and ICAM-2 with a high number of interactors in Protein–Protein Interaction (PPI) Network indicated by bioinformatics report. Furthermore, using ELISA test, we confirmed that all increase in the levels of NGF, NT-3, IGF-2, HGF, NPY, CRP, MMP-9, and ICAM-2 in the serum from TBI patients. Together, we determined the screened protein expressional profiles in serum for TBI patients, in which the cross-network between inflammatory factors and growth factors may play a crucial role in TBI damage and repair. Our findings could contribute to indication for the diagnosis and treatment of TBI in future translational medicine and clinical practice.

Essential molecular functions associated with the circular code evolution

A circular code is a set of trinucleotides allowing the reading frames in genes to be retrieved locally, i.e. anywhere in genes and in particular without start codons, and automatically with a window of few nucleotides. In 1996, a common circular code, called X, was identified in large populations of eukaryotic and prokaryotic genes. Hence, it is believed to be an ancestral structural property of genes. A new computational approach based on comparative genomics is developed to identify essential molecular functions associated with circular codes. It is based on a quantitative and sensitive statistical method (FPTF) to identify three permuted trinucleotide sets in the three frames of genes, a flower automaton algorithm to determine if a trinucleotide set is a circular code or not, and an integrated Gene Ontology and Taxonomy (iGOT) database. By carrying out automatic circular code analyses on a huge number of gene populations where each population is associated with a particular molecular function, it identifies 266 gene populations having circular codes close to X. Surprisingly, their molecular functions include 98% of those covered by the essential genes of the DEG database (Database of Essential Genes). Furthermore, three trinucleotides GTG, AAG and GCG, replacing three trinucleotides of the code X and called “evolutionary” trinucleotides, significantly occur in these 266 gene populations. Finally, a new method developed to analyse and quantify the stability of a set of trinucleotides demonstrates that these evolutionary trinucleotides are associated with a significant increase of the stability of the common circular code X. Indeed, its stability increases from the 1502th rank to the 16th rank after the replacement of the three evolutionary trinucleotides among 9920 possible trinucleotide replacement sets.

Dynamic Behavior of an Intrinsically Unstructured Linker Domain Is Conserved in the Face of Negligible Amino Acid Sequence Conservation

Proteins or regions of proteins that do not form compact globular structures are classified as intrinsically unstructured proteins (IUPs). IUPs are common in nature and have essential molecular functions, but even a limited understanding of the evolution of their dynamic behavior is lacking. The primary objective of this work was to test the evolutionary conservation of dynamic behavior for a particular class of IUPs that form intrinsically unstructured linker domains (IULD) that tether flanking folded domains. This objective was accomplished by measuring the backbone flexibility of several IULD homologues using nuclear magnetic resonance (NMR) spectroscopy. The backbone flexibility of five IULDs, representing three kingdoms, was measured and analyzed. Two IULDs from animals, one IULD from fungi, and two IULDs from plants showed similar levels of backbone flexibility that were consistent with the absence of a compact globular structure. In contrast, the amino acid sequences of the IULDs from these three taxa showed no significant similarity. To investigate how the dynamic behavior of the IULDs could be conserved in the absence of detectable sequence conservation, evolutionary rate studies were performed on a set of nine mammalian IULDs. The results of this analysis showed that many sites in the IULD are evolving neutrally, suggesting that dynamic behavior can be maintained in the absence of natural selection. This work represents the first experimental test of the evolutionary conservation of dynamic behavior and demonstrates that amino acid sequence conservation is not required for the conservation of dynamic behavior and presumably molecular function.

Mathematical Models of Host Plant Infection by Helper-Dependent Virus Complexes: Why Are Helper Viruses Always Avirulent?

ABSTRACT Interactions between viruses in plants are common, and some viruses depend on such interactions for their survival. Frequently, a virus lacks some essential molecular function that another provides. In "helper-dependent" virus complexes, the helper virus is transmitted independently by a vector, whereas the dependent virus depends on molecular agents associated with the helper virus for transmission by a vector. A general mathematical model was developed of the dynamics of host plant infection by a helper-dependent virus complex. Four categories of host plants were considered: healthy, infected with helper virus alone, infected with dependent virus alone, and infected with both viruses. New planting of the host crop was constrained by a maximum abundance due to limitation of the cropping area. The ratio of infection rate to host loss rate due to infection is proposed as an important epidemiological quantity, A, that can be used as a measure of the mutual adaptation of the virus and host. A number of alternative equilibria of host infection could occur and were determined exclusively by parameter values; it was informative to display their distribution in the parameter plane: (1/A)(helper) versus (1/A)(dependent). A simple analysis of the distribution of the final equilibria illustrated that the dependent virus could affect the survival of the helper virus, so facilitation between the two can be reciprocal. The distribution of the final equilibria also indicated that a well-adapted helper virus increases the opportunity for a dependent virus to evolve and survive, and the model, therefore, explains why infection with a helper virus usually causes no or little damage to plants, whereas infection with a dependent virus or mixed infection with both often causes very severe damage.