RNA-binding Domain Research Articles

Gadd45β is critical for regulation of type I interferon signaling by facilitating G3BP-mediated stress granule formation

Stress granules (SGs) constitute a signaling hub that plays a critical role in type I interferon responses. Here, we report that growth arrest and DNA damage-inducible beta (Gadd45β) act as a positive regulator of SG-mediated interferon signaling by targeting G3BP upon RNA virus infection. Gadd45β deficiency markedly impairs SG formation and SG-mediated activation of interferon signaling invitro. Gadd45β knockout mice are highly susceptible to RNA virus infection, and their ability to produce interferon and cytokines is severely impaired. Specifically, Gadd45β interacts with the RNA-binding domain of G3BP, leading to conformational expansion of G3BP1 via dissolution of its autoinhibitory electrostatic intramolecular interaction. The acidic loop 1- and RNA-binding properties of Gadd45β markedly increase the conformational expansion and RNA-binding affinity of the G3BP1-Gadd45β complex, thereby promoting assembly of SGs. These findings suggest a role for Gadd45β as a component and critical regulator of G3BP1-mediated SG formation, which facilitates RLR-mediated interferon signaling.

Molecular Docking and Dynamics Identify Potential Drugs to be Repurposed as SARS-CoV-2 Inhibitors

The novel coronavirus disease 19 (COVID-19) has resulted in an estimated 20 million excess deaths and the recent resurgence of COVID-19 in China is predicted to result in up to 1 million deaths over the next few months. With vaccines being ineffective in the case of immunocompromised patients, it is important to continue our quest for safe, effective and affordable drugs that will be available to all countries. Drug repurposing is one of the strategies being explored in this context. Recently, out of the 7817 drugs approved worldwide, 214 candidates were systematically down-selected using a combination of 11 filters including FDA/TGA approval status, assay data against SARS-CoV-2, pharmacokinetic, pharmacodynamic and toxicity profiles. These down-selected drugs were subjected in this study to virtual screening against various SARS-CoV-2 targets followed by molecular dynamics studies of the best scoring ligands against each target. The chosen molecular targets were spike receptor binding domain, nucleocapsid protein RNA binding domain and key nonstructural proteins 3, 5 and 12–14. Four drugs approved for other indications — alendronate, cromolyn, natamycin and treprostinil — look sufficiently promising from our in-silico studies to warrant further in-vitro and in-vivo investigations as appropriate to ascertain their extent of antiviral activities.

Split Luciferase-Fragment Reconstitution for Unveiling RNA Localization and Dynamics in Live Cells.

The intracellular distribution and dynamics of RNAs play pivotal roles in various physiological phenomena. The ability to monitor the amount and localization of endogenous RNAs in living cells allows for elucidating the mechanisms of various intracellular events. Protein-based fluorescent RNA probes are now widely used to visualize and analyze RNAs in living cells. However, continuously monitoring the temporal changes in RNA localization and dynamics in living cells is challenging. In this study, we developed a bioluminescent probe for spatiotemporal monitoring of RNAs in living cells by using a split-luciferase reconstitution technique. The probe consists of split fragments of a bioluminescent protein, NanoLuc, connected with RNA-binding protein domains generated from a custom-made mutation of a PUM-HD. The probe showed rapid luminescence intensity changes in response to an increase or decrease in the amount of a target RNA in vitro. In live-cell imaging, temporal alteration of the intracellular distribution of endogenous β-actin mRNA was visualized in response to extracellular stimulation. Furthermore, the application of the probe to the visualization of the specific localization of β-actin mRNA in primary hippocampal neurons was conducted. These results demonstrate the capability of the bioluminescent RNA probe to monitor the changes in localization, dynamics, and the amount of target RNA in living cells.

Assembly of SARS-CoV-2 ribonucleosomes by truncated N∗ variant of the nucleocapsid protein

The nucleocapsid (N) protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) compacts the RNA genome into viral ribonucleoprotein (vRNP) complexes within virions. Assembly of vRNPs is inhibited by phosphorylation of the N protein serine/arginine (SR) region. Several SARS-CoV-2 variants of concern carry N protein mutations that reduce phosphorylation and enhance the efficiency of viral packaging. Variants of the dominant B.1.1 viral lineage also encode a truncated N protein, termed N∗ or Δ(1-209), that mediates genome packaging despite lacking the N-terminal RNA-binding domain and SR region. Here, we use mass photometry and negative stain electron microscopy to show that purified Δ(1-209) and viral RNA assemble into vRNPs that are remarkably similar in size and shape to those formed with full-length N protein. We show that assembly of Δ(1-209) vRNPs requires the leucine-rich helix of the central disordered region and that this helix promotes N protein oligomerization. We also find that fusion of a phosphomimetic SR region to Δ(1-209) inhibits RNA binding and vRNP assembly. Our results provide new insights into the mechanisms by which RNA binding promotes N protein self-association and vRNP assembly, and how this process is modulated by phosphorylation.

The Saccharomyces cerevisiae SR protein Npl3 interacts with hyperphosphorylated CTD of RNA Polymerase II

The catalytic subunit of RNA Polymerase II contains a highly conserved carboxy terminal domain (CTD) composed of multiple tandem heptad sequence Tyr1Ser2Pro3Thr4Ser5Pro6Ser7. The non-proline residues in CTD undergo posttranslational modifications, with Ser5 phosphorylation (Ser5P) predominating at the start of the transcription cycle and Ser2P at the end, while other phosphorylation levels are high all throughout. The differentially phosphorylated CTD is recognized by regulatory proteins, helpful during mRNA transcription and export. One such protein Npl3 is composed of two RNA binding domains and a C-terminus RGG/SR domain. The Ser411 of Npl3 is reported to make direct contact with Ser2P of CTD for its recruitment and function, while the Npl3 lacking of C-terminal 25 amino acids (Npl3Δ389–414) showed no apparent defects in mRNA synthesis. Here, we report that the RNA binding domains of Npl3 are separate folding units and interact also with the CTD. The interaction between Npl3 and CTD appears to involve not just Ser2P, but also the Ser5P and Ser7P. The Arg126 of the first RNA binding domain interacts with Ser2P whereas the Arg235 of the second RNA binding domain interacts with either Ser7P or Ser5P of another heptad. The finding provides new insight of Npl3 function for mRNA transcription.

ADAR1 links R-loop homeostasis to ATR activation in replication stress response.

Unscheduled R-loops are a major source of replication stress and DNA damage. R-loop-induced replication defects are sensed and suppressed by ATR kinase, whereas it is not known whether R-loop itself is actively involved in ATR activation and, if so, how this is achieved. Here, we report that the nuclear form of RNA-editing enzyme ADAR1 promotes ATR activation and resolves genome-wide R-loops, a process that requires its double-stranded RNA-binding domains. Mechanistically, ADAR1 interacts with TOPBP1 and facilitates its loading on perturbed replication forks by enhancing the association of TOPBP1 with RAD9 of the 9-1-1 complex. When replication is inhibited, DNA-RNA hybrid competes with TOPBP1 for ADAR1 binding to promote the translocation of ADAR1 from damaged fork to accumulate at R-loop region. There, ADAR1 recruits RNA helicases DHX9 and DDX21 to unwind R-loops, simultaneously allowing TOPBP1 to stimulate ATR more efficiently. Collectively, we propose that the tempo-spatially regulated assembly of ADAR1-nucleated protein complexes link R-loop clearance and ATR activation, while R-loops crosstalk with blocked replication forks by transposing ADAR1 to finetune ATR activity and safeguard the genome.

Phytocompound screening, antioxidant activity and molecular docking studies of pomegranate seed: a preventive approach for SARS-CoV-2 pathogenesis

A global hazard to public health has been generated by the coronavirus infection 2019 (COVID-19), which is spreading quickly. Pomegranate is a strong source of antioxidants and has demonstrated a number of pharmacological characteristics. This work was aimed to analyze the phytochemicals present in ethanolic pomegranate seed extract (PSE) and their in vitro antioxidant potential and further in-silico evaluation for antiviral potential against crystal structure of two nucleocapsid proteins i.e., N-terminal RNA binding domain (NRBD) and C-terminal Domain (CTD) of SARS-CoV-2. The bioactive components from ethanolic extract of PSE were assessed by gas chromatography-mass spectroscopy (GC–MS). Free radical scavenging activity of PSE was determined using DPPH dye. Molecular docking was executed through the Glide module of Maestro software. Lipinski’s 5 rule was applied for drug-likeness characteristics using cheminformatics Molinspiration software while OSIRIS Data Warrior V5.5.0 was used to predict possible toxicological characteristics of components. Thirty-two phytocomponents was detected in PSE by GC–MS technique. Free radical scavenging assay revealed the high antioxidant capacity of PSE. Docking analysis showed that twenty phytocomponents from PSE exhibited good binding affinity (Docking score ≥ − 1.0 kcal/mol) towards NRBD and CTD nucleocapsid protein. This result increases the possibility that the top 20 hits could prevent the spread of SARS-CoV-2 by concentrating on both nucleocapsid proteins. Moreover, molecular dynamics (MD) simulation using GROMACS was used to check their binding efficacy and internal dynamics of top complexes with the lowest docking scores. The metrics root mean square deviation (RMSD), root mean square fluctuation (RMSF), intermolecular hydrogen bonding (H-bonds) and radius of gyration (Rg) revealed that the lead phytochemicals form an energetically stable complex with the target protein. Majority of the phytoconstituents exhibited drug-likeness with non-tumorigenic properties. Thus, the PSE phytoconstituents could be useful source of drug or nutraceutical development in SARS-CoV-2 pathogenesis.

Chemoproteomic capture of RNA binding activity in living cells

Proteomic methods for RNA interactome capture (RIC) rely principally on crosslinking native or labeled cellular RNA to enrich and investigate RNA-binding protein (RBP) composition and function in cells. The ability to measure RBP activity at individual binding sites by RIC, however, has been more challenging due to the heterogenous nature of peptide adducts derived from the RNA-protein crosslinked site. Here, we present an orthogonal strategy that utilizes clickable electrophilic purines to directly quantify protein-RNA interactions on proteins through photoaffinity competition with 4-thiouridine (4SU)-labeled RNA in cells. Our photo-activatable-competition and chemoproteomic enrichment (PACCE) method facilitated detection of >5500 cysteine sites across ~3000 proteins displaying RNA-sensitive alterations in probe binding. Importantly, PACCE enabled functional profiling of canonical RNA-binding domains as well as discovery of moonlighting RNA binding activity in the human proteome. Collectively, we present a chemoproteomic platform for global quantification of protein-RNA binding activity in living cells.

A transcriptional network required for bradyzoite development in Toxoplasma gondii is dispensable for recrudescent disease

Identification of regulators of Toxoplasma gondii bradyzoite development and cyst formation is the most direct way to address the importance of parasite development in long-term persistence and reactivation of this parasite. Here we show that a T. gondii gene (named Regulator of Cystogenesis 1; ROCY1) is sufficient for T. gondii bradyzoite formation in vitro and in vivo. ROCY1 encodes an RNA binding protein that has a preference for 3’ regulatory regions of hundreds of T. gondii transcripts, and its RNA-binding domains are required to mediate bradyzoite development. Female mice infected with ΔROCY1 parasites have reduced (>90%) cyst burden. While viable parasites can be cultivated from brain tissue for up to 6 months post-infection, chronic brain-resident ΔROCY1 parasites have reduced oral infectivity compared to wild type. Despite clear defects in bradyzoite formation and oral infectivity, ΔROCY1 parasites were able to reactivate with similar timing and magnitude as wild type parasites for up to 5 months post-infection. Therefore while ROCY1 is a critical regulator of the bradyzoite developmental pathway, it is not required for parasite reactivation, raising new questions about the persisting life stage responsible for causing recrudescent disease.

Biochemical and structural characterization of the RT domain of Leishmania sp. telomerase reverse transcriptase

BackgroundThe Leishmania genus comprises parasites that cause leishmaniasis, a neglected disease spread worldwide. Leishmania sp. telomeres are composed of TTAGGG repeats maintained by telomerase. In most eukaryotes, the enzyme minimal complex contains the TER (telomerase RNA) and the TERT (telomerase reverse transcriptase) components. The TERT holds the enzyme catalytic core and is formed by four structural and functional domains (TEN, Telomerase Essential N-terminal; TRBD, Telomerase RNA Binding Domain; RT, the reverse transcriptase domain and CTE, C-Terminal Extension domain). Methods and resultsAmino acid sequence alignments, protein structure prediction analysis, and protein: nucleic acid interaction assays were used to show that the Leishmania major RT domain preserves the canonical structural elements found in higher eukaryotes, including the canonical motifs and the aspartic acid residues that stabilize the Mg2+ ion cofactor. Furthermore, amino acid substitutions specific to the Leishmania genus and partial conservation of the residues involved with nucleic acid interactions are shown. The purified recombinant Leishmania RT protein is biochemically active and interacts with the G-rich telomeric strand and the TER template sequence. ConclusionOur results highlight that the telomerase catalysis mechanism is conserved in a pathogen of medical importance despite the structural peculiarities present in the parasite's RT domain.

The characterization of RNA-binding proteins and RNA metabolism-related proteins in fungal extracellular vesicles.

RNA-binding proteins (RBPs) are essential for regulating RNA metabolism, stability, and translation within cells. Recent studies have shown that RBPs are not restricted to intracellular functions and can be found in extracellular vesicles (EVs) in different mammalian cells. EVs released by fungi contain a variety of proteins involved in RNA metabolism. These include RNA helicases, which play essential roles in RNA synthesis, folding, and degradation. Aminoacyl-tRNA synthetases, responsible for acetylating tRNA molecules, are also enriched in EVs, suggesting a possible link between these enzymes and tRNA fragments detected in EVs. Proteins with canonical RNA-binding domains interact with proteins and RNA, such as the RNA Recognition Motif (RRM), Zinc finger, and hnRNP K-homology (KH) domains. Polyadenylate-binding protein (PABP) plays a critical role in the regulation of gene expression by binding the poly(A) tail of messenger RNA (mRNA) and facilitating its translation, stability, and localization, making it a key factor in post-transcriptional control of gene expression. The presence of proteins related to the RNA life cycle in EVs from different fungal species suggests a conserved mechanism of EV cargo packing. Various models have been proposed for selecting RNA molecules for release into EVs. Still, the actual loading processes are unknown, and further molecular characterization of these proteins may provide insight into the mechanism of RNA sorting into EVs. This work reviews the current knowledge of RBPs and proteins related to RNA metabolism in EVs derived from distinct fungi species, and presents an analysis of proteomic datasets through GO term and orthology analysis, Our investigation identified orthologous proteins in fungal EVs on different fungal species.

Suppression of hnRNP A1 binding to HK1 RNA leads to glycolytic dysfunction in Alzheimer's disease models.

To investigate the mechanism of RNA-binding protein hnRNP A1 in mouse hippocampal neurons (HT22) on glycolysis. RIP and CLIP-qPCR were performed by HT22 in vitro to observe the mechanism of hnRNP A1 regulating the expression of key proteins in glycolysis. The RNA binding domain of hnRNP A1 protein in HT22 was inhibited by VPC-80051, and the effect of hnRNP A1 on glycolysis of HT22 was observed. Lentivirus overexpression of hnRNP A1 was used to observe the effect of overexpression of hnRNP A1 on glycolysis of Aβ25-35-injured HT22. The expression of hnRNP A1 in brain tissues of wild-type mice and triple-transgenic (APP/PS1/Tau) AD mice at different ages was studied by Western blot assay. The results of RIP experiment showed that hnRNP A1 and HK1 mRNA were significantly bound. The results of CLIP-qPCR showed that hnRNP A1 directly bound to the 2605-2821 region of HK1 mRNA. hnRNP A1 inhibitor can down-regulate the expression of HK1 mRNA and HK1 protein in HT22 cells. Overexpression of hnRNP A1 can significantly reduce the toxic effect of Aβ25-35 on neurons via the hnRNP A1/HK1/ pyruvate pathway. In addition, inhibition of hnRNP A1 binding to amyloid precursor protein (APP) RNA was found to increase Aβ expression, while Aβ25-35 also down-regulated hnRNP A1 expression by enhancing phosphorylation of p38 MAPK in HT22. They interact to form bidirectional regulation, further down-regulating the expression of hnRNP A1, and ultimately aggravating glycolytic dysfunction. Protein immunoblotting showed that hnRNP A1 decreased with age in mouse brain tissue, and the decrease was greater in AD mice, suggesting that the decrease of hnRNP A1 may be a predisposed factor in the pathogenesis of AD.

SARS-CoV-2 nucleocapsid protein inhibits the PKR-mediated integrated stress response through RNA-binding domain N2b.

The nucleocapsid protein N of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) enwraps and condenses the viral genome for packaging but is also an antagonist of the innate antiviral defense. It suppresses the integrated stress response (ISR), purportedly by interacting with stress granule (SG) assembly factors G3BP1 and 2, and inhibits type I interferon responses. To elucidate its mode of action, we systematically deleted and over-expressed distinct regions and domains. We show that N via domain N2b blocks PKR-mediated ISR activation, as measured by suppression of ISR-induced translational arrest and SG formation. N2b mutations that prevent dsRNA binding abrogate these activities also when introduced in the intact N protein. Substitutions reported to block post-translation modifications of N or its interaction with G3BP1/2 did not have a detectable additive effect. In an encephalomyocarditis virus-based infection model, N2b - but not a derivative defective in RNA binding-prevented PKR activation, inhibited β-interferon expression and promoted virus replication. Apparently, SARS-CoV-2 N inhibits innate immunity by sequestering dsRNA to prevent activation of PKR and RIG-I-like receptors. Similar observations were made for the N protein of human coronavirus 229E, suggesting that this may be a general trait conserved among members of other orthocoronavirus (sub)genera.

Cell-type specific regulator RBPMS switches alternative splicing via higher-order oligomerization and heterotypic interactions with other splicing regulators.

Alternative pre-mRNA splicing decisions are regulated by RNA binding proteins (RBPs) that can activate or repress regulated splice sites. Repressive RBPs typically harness multivalent interactions to bind stably to target RNAs. Multivalency can be achieved by homomeric oligomerization and heteromeric interactions with other RBPs, often mediated by intrinsically disordered regions (IDRs), and by possessing multiple RNA binding domains. Cell-specific splicing decisions often involve the action of widely expressed RBPs, which are able to bind multivalently around target exons, but without effect in the absence of a cell-specific regulator. To address how cell-specific regulators can collaborate with constitutive RBPs in alternative splicing regulation, we used the smooth-muscle specific regulator RBPMS. Recombinant RBPMS is sufficient to confer smooth muscle cell specific alternative splicing of Tpm1 exon 3 in cell-free assays by preventing assembly of ATP-dependent splicing complexes. This activity depends upon a C-terminal IDR that facilitates dynamic higher-order self-assembly, cooperative binding to multivalent RNA and interactions with widely expressed splicing co-regulators, including MBNL1 and RBFOX2, allowing cooperative assembly of stable cell-specific regulatory complexes.

Examination of Factors Affecting Site-Directed RNA Editing by the MS2-ADAR1 Deaminase System.

Adenosine deaminases acting on RNA (ADARs) have double-stranded RNA binding domains and a deaminase domain (DD). We used the MS2 system and specific guide RNAs to direct ADAR1-DD to target adenosines in the mRNA encoding-enhanced green fluorescence protein. Using this system in transfected HEK-293 cells, we evaluated the effects of changing the length and position of the guide RNA on the efficiency of conversion of amber (TAG) and ochre (TAA) stop codons to tryptophan (TGG) in the target. Guide RNAs of 19, 21 and 23 nt were positioned upstream and downstream of the MS2-RNA, providing a total of six guide RNAs. The upstream guide RNAs were more functionally effective than the downstream guide RNAs, with the following hierarchy of efficiency: 21 nt > 23 nt > 19 nt. The highest editing efficiency was 16.6%. Off-target editing was not detected in the guide RNA complementary region but was detected 50 nt downstream of the target. The editing efficiency was proportional to the amount of transfected deaminase but inversely proportional to the amount of the transfected guide RNA. Our results suggest that specific RNA editing requires precise optimization of the ratio of enzyme, guide RNA, and target RNA.

RNA-binding domain of SARS-CoV2 nucleocapsid: MD simulation study of the effect of the proline substitutions P67S and P80R on the structure of the protein

The nucleocapsid component of SARS-CoV2 is involved in the viral genome packaging. GammaP.1(Brazil) and the 20 C-US(USA) variants had a high frequency of the P80R and P67S mutations respectively in the RNA-binding domain of the nucleocapsid. Since RNA-binding domain participates in the electrostatic interactions with the viral genome, the study of the effects of proline substitutions on the flexibility of the protein will be meaningful. It evinced that the trajectory of the wildtype and mutants was stable during the simulation and exhibited distinct changes in the flexibility of the protein. Moreover, the beta-hairpin loop region of the protein structures exhibited high amplitude fluctuations and dominant motions. Additionally, modulations were detected in the drug binding site. Besides, the extent of correlation and anti-correlation motions involving the protruding region, helix, and the other RNA binding sites differed between the wildtype and mutants. The secondary structure analysis disclosed the variation in the occurrence pattern of the secondary structure elements between the proteins. Protein-ssRNA interaction analysis was also done to detect the amino acid contacts with ssRNA. R44, R59, and Y61 residues of the wildtype and P80R mutant exhibited different duration contacts with the ssRNA. It was also noticed that R44, R59, and Y61 of the wildtype and P80R formed hydrogen bonds with the ssRNA. However in P67S, residues T43, R44, R45, R40, R59, and R41 displayed contacts and formed hydrogen bonds with ssRNA. Binding free energy was also calculated and was lowest for P67S than wildtype andP80R. Thus, proline substitutions influence the structure of the RNA-binding domain and may modulate viral genome packaging besides the host-immune response. Communicated by Ramaswamy H. Sarma

Structural basis for specific RNA recognition by the alternative splicing factor RBM5

The RNA-binding motif protein RBM5 belongs to a family of multi-domain RNA binding proteins that regulate alternative splicing of genes important for apoptosis and cell proliferation and have been implicated in cancer. RBM5 harbors structural modules for RNA recognition, such as RRM domains and a Zn finger, and protein-protein interactions such as an OCRE domain. Here, we characterize binding of the RBM5 RRM1-ZnF1-RRM2 domains to cis-regulatory RNA elements. A structure of the RRM1-ZnF1 region in complex with RNA shows how the tandem domains cooperate to sandwich target RNA and specifically recognize a GG dinucleotide in a non-canonical fashion. While the RRM1-ZnF1 domains act as a single structural module, RRM2 is connected by a flexible linker and tumbles independently. However, all three domains participate in RNA binding and adopt a closed architecture upon RNA binding. Our data highlight how cooperativity and conformational modularity of multiple RNA binding domains enable the recognition of distinct RNA motifs, thereby contributing to the regulation of alternative splicing. Remarkably, we observe surprising differences in coupling of the RNA binding domains between the closely related homologs RBM5 and RBM10.

HKDC1 reprograms lipid metabolism to enhance gastric cancer metastasis and cisplatin resistance via forming a ribonucleoprotein complex

As essential modulators of transcription and translation, RNA-binding proteins (RBPs) are frequently dysregulated in cancer. Bioinformatics study reveals that the RNA-binding protein hexokinase domain component 1 (HKDC1) is overexpressed in gastric cancer (GC). As HKDC1 plays a role in lipid homeostasis in the liver and glucose metabolism in certain cancers, the exact mechanism of action of HKDC1 in GC remains largely unknown. Upregulation of HKDC1 correlates with chemoresistance and poor prognosis in GC patients. HKDC1 enhances invasion, migration and resistance to cisplatin (CDDP) in GC cells in vitro and in vivo. Comprehensive transcriptomic sequencing and metabolomic analysis reveal that HKDC1 mediates abnormal lipid metabolism in GC cells. Herein, we identify a number of HKDC1–binding endogenous RNAs in GC cells, including protein kinase, DNA-activated, catalytic subunit (PRKDC) mRNA. We further validate that PRKDC is a crucial downstream effector of HKDC1 induced-GC tumorigenesis depends on lipid metabolism. Interestingly, G3BP1, a well-known oncoprotein, can be bound by HKDC1. HKDC1 cooperates with G3BP1 to enhance the stability of PRKDC transcript. Our results reveal a novel HKDC1/G3BP1-PRKDC regulatory axis that induces GC metastasis and chemoresistance via reprogramming lipid metabolism, which may provide an effective therapeutic strategy for a subset of GC with HKDC1 overexpression.

An Evidence of Drug Repurposing for COVID-19 Pandemic Based on <i>In silico</i> Investigation from Phenolic Derivatives of Silybum Marianum Against SARS-Cov-2 Proteins

The outbreak of coronavirus disease-2019 (COVID-19) had a striking impact on the worldwide healthcare system within a very short period. The availability of a large number of clinical data on SARS-CoV-2, conventional precautionary majors, and treatment strategies with the existing therapeutic antiviral drug molecules also fails to control progression and disease transmission among the population. Hence, we implemented pharmacoinformatics approaches to facilitate the drug discovery by repurposing naturally available therapeutic molecules as an effective intervention. The major phenolic derivatives of <em>Silybum marianum </em>(Milk thistle) have been identified and investigated for ADME (Absorption, Distribution, Metabolism and Excretion)/tox properties. Co-crystallized structure of three major proteins (i.e., main protease, RNA binding domain of nucleocapsid phosphoprotein and Spike receptor binding domain) from SARS-CoV-2 investigated with molecular docking (MD) interaction with the phenolic compounds from milk thistle. Furthermore, a 100 ns MD simulation was performed with silibinin molecule based on ADMET and MD interaction. Being less toxic in ADME, a good MD interaction and stability of silibinin molecule across the MD simulation trajectories with targeted proteins explicate that silibinin molecule can be a promising drug candidate against the main protease and will be helpful to cease the enzymatic activity in viral replication and transcription.

HydRA: Deep-learning models for predicting RNA-binding capacity from protein interaction association context and protein sequence

RNA-binding proteins (RBPs) control RNA metabolism to orchestrate gene expression and, when dysfunctional, underlie human diseases. Proteome-wide discovery efforts predict thousands of RBP candidates, many of which lack canonical RNA-binding domains (RBDs). Here, we present a hybrid ensemble RBP classifier (HydRA), which leverages information from both intermolecular protein interactions and internal protein sequence patterns to predict RNA-binding capacity with unparalleled specificity and sensitivity using support vector machines (SVMs), convolutional neural networks (CNNs), and Transformer-based protein language models. Occlusion mapping by HydRA robustly detects known RBDs and predicts hundreds of uncharacterized RNA-binding associated domains. Enhanced CLIP (eCLIP) for HydRA-predicted RBP candidates reveals transcriptome-wide RNA targets and confirms RNA-binding activity for HydRA-predicted RNA-binding associated domains. HydRA accelerates construction of a comprehensive RBP catalog and expands the diversity of RNA-binding associated domains.