Structural and Dynamical Basis of VP35-RBD Inhibition by Marine Fungi Compounds to Combat Marburg Virus Infection

The Marburg virus (MBV), a deadly pathogen, poses a serious threat to world health due to the lack of effective treatments, calling for an immediate search for targeted and efficient treatments. In this study, we focused on compounds originating from marine fungi in order to identify possible inhibitory compounds against the Marburg virus (MBV) VP35-RNA binding domain (VP35-RBD) using a computational approach. We started with a virtual screening procedure using the Lipinski filter as a guide. Based on their docking scores, 42 potential candidates were found. Four of these compounds-CMNPD17596, CMNPD22144, CMNPD25994, and CMNPD17598-as well as myricetin, the control compound, were chosen for re-docking analysis. Re-docking revealed that these particular compounds had a higher affinity for MBV VP35-RBD in comparison to the control. Analyzing the chemical interactions revealed unique binding properties for every compound, identified by a range of Pi-cation interactions and hydrogen bond types. We were able to learn more about the dynamic behaviors and stability of the protein-ligand complexes through a 200-nanosecond molecular dynamics simulation, as demonstrated by the compounds' consistent RMSD and RMSF values. The multidimensional nature of the data was clarified by the application of principal component analysis, which suggested stable conformations in the complexes with little modification. Further insight into the energy profiles and stability states of these complexes was also obtained by an examination of the free energy landscape. Our findings underscore the effectiveness of computational strategies in identifying and analyzing potential inhibitors for MBV VP35-RBD, offering promising paths for further experimental investigations and possible therapeutic development against the MBV.

Mechanistic and thermodynamic characterization of antiviral inhibitors targeting nucleocapsid N-terminal domain of SARS-CoV-2

The nucleocapsid (N) protein of SARS-CoV-2 plays a pivotal role in encapsulating the viral genome. Developing antiviral treatments for SARS-CoV-2 is imperative due to the diminishing immunity of the available vaccines. This study targets the RNA-binding site located in the N-terminal domain (NTD) of the N-protein to identify the potential antiviral molecules against SARS-CoV-2. A structure-based repurposing approach identified the twelve high-affinity molecules from FDA-approved drugs, natural products, and the LOPAC1280 compound libraries that precisely bind to the RNA binding site within the NTD. The interaction of these potential antiviral agents with the purified NTD protein was thermodynamically characterized using isothermal titration calorimetry (ITC). A fluorescence-based plate assay to assess the RNA binding inhibitory activity of small molecules against the NTD has been employed, and the selected compounds exhibited significant RNA binding inhibition with calculated IC50 values ranging from 8.8 μM to 15.7 μM. Furthermore, the antiviral efficacy of these compounds was evaluated using in vitro cell-based assays targeting the replication of SARS-CoV-2. Remarkably, two compounds, Telmisartan and BMS-189453, displayed potential antiviral activity against SARS-CoV-2, with EC50 values of approximately 1.02 μM and 0.98 μM, and a notable selective index of >98 and > 102, respectively. This study gives valuable insight into developing therapeutic interventions against SARS-CoV-2 by targeting the N-protein, a significant effort given the global public health concern posed due to the virus re-emergence and long COVID-19 disease.

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.

In silico study predicts a key role of RNA-binding domains 3 and 4 in nucleolin-miRNA interactions

RNA binding proteins (RBPs) regulate many important cellular processes through their interactions with RNA molecules. RBPs are critical for posttranscriptional mechanisms keeping gene regulation in a fine equilibrium. Conversely, dysregulation of RBPs and RNA metabolism pathways is an established hallmark of tumorigenesis. Human nucleolin (NCL) is a multifunctional RBP that interacts with different types of RNA molecules, in part through its four RNA binding domains (RBDs). Particularly, NCL interacts directly with microRNAs (miRNAs) and is involved in their aberrant processing linked with many cancers, including breast cancer. Nonetheless, molecular details of the NCL-miRNA interaction remain obscure. In this study, we used an in silico approach to characterize how NCL targets miRNAs and whether this specificity is imposed by a definite RBD-interface. Here, we present structural models of NCL-RBDs and miRNAs, as well as predict scenarios of NCL-miRNA interactions generated using docking algorithms. Our study suggests a predominant role of NCL RBDs 3 and 4 (RBD3-4) in miRNA binding. We provide detailed analyses of specific motifs/residues at the NCL-substrate interface in both these RBDs and miRNAs. Finally, we propose that the evolutionary emergence of more than two RBDs in NCL in higher organisms coincides with its additional role/s in miRNA processing. Our study shows that RBD3-4 display sequence/structural determinants to specifically recognize miRNA precursor molecules. Moreover, the insights from this study can ultimately support the design of novel antineoplastic drugs aimed at regulating NCL-dependent biological pathways with a causal role in tumorigenesis.

1H, 13C and 15N resonance assignments and solution structures of the two RRM domains of Matrin-3

Matrin-3 is a multifunctional protein that binds to both DNA and RNA. Its DNA-binding activity is linked to the formation of the nuclear matrix and transcriptional regulation, while its RNA-binding activity is linked to mRNA metabolism including splicing, transport, stabilization, and degradation. Correspondingly, Matrin-3 has two zinc finger domains for DNA binding and two consecutive RNA recognition motif (RRM) domains for RNA binding. Matrin-3 has been reported to cause amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) when its disordered region contains pathogenic mutations. Simultaneously, it has been shown that the RNA-binding activity of Matrin-3 mediated by its RRM domains, affects the formation of insoluble cytoplasmic granules, which are related to the pathogenic mechanism of ALS/FTD. Thus, the effect of the RRM domains on the phase separation of condensed protein/RNA mixtures has to be clarified for a comprehensive understanding of ALS/FTD. Here, we report the ¹H, ¹⁵N, and ¹³C resonance assignments of the two RNA binding domains and their solution structures. The resonance assignments and the solution structures obtained in this work will contribute to the elucidation of the molecular basis of Matrin-3 in the pathogenic mechanism of ALS and/or FTD.

Assessing nucleic acid binding activity of four dinoflagellate cold shock domain proteins from Symbiodinium kawagutii and Lingulodinium polyedra

Background

Dinoflagellates have a generally large number of genes but only a small percentage of these are annotated as transcription factors. Cold shock domain (CSD) containing proteins (CSPs) account for roughly 60% of these. CSDs are not prevalent in other eukaryotic lineages, perhaps suggesting a lineage-specific expansion of this type of transcription factors in dinoflagellates, but there is little experimental data to support a role for dinoflagellate CSPs as transcription factors. Here we evaluate the hypothesis that dinoflagellate CSPs can act as transcription factors by binding double-stranded DNA in a sequence dependent manner.

Results

We find that both electrophoretic mobility shift assay (EMSA) competition experiments and selection and amplification binding (SAAB) assays indicate binding is not sequence specific for four different CSPs from two dinoflagellate species. Competition experiments indicate all four CSPs bind to RNA better than double-stranded DNA.

Conclusion

Dinoflagellate CSPs do not share the nucleic acid binding properties expected for them to function as bone fide transcription factors. We conclude the transcription factor complement of dinoflagellates is even smaller than previously thought suggesting that dinoflagellates have a reduced dependance on transcriptional control compared to other eukaryotes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12860-021-00368-4.

Optimal RNA binding by Egalitarian, a Dynein cargo adaptor, is critical for maintaining oocyte fate in Drosophila

The Dynein motor is responsible for the localization of numerous mRNAs within Drosophila oocytes and embryos. The RNA binding protein, Egalitarian (Egl), is thought to link these various RNA cargoes with Dynein. Although numerous studies have shown that Egl is able to specifically associate with these RNAs, the nature of these interactions has remained elusive. Egl contains a central RNA binding domain that shares limited homology with an exonuclease, yet Egl binds to RNA without degrading it. Mutations have been identified within Egl that disrupt its association with its protein interaction partners, BicaudalD (BicD) and Dynein light chain (Dlc), but no mutants have been described that are specifically defective for RNA binding. In this report, we identified a series of positively charged residues within Egl that are required for RNA binding. Using corresponding RNA binding mutants, we demonstrate that specific RNA cargoes are more reliant on maximal Egl RNA biding activity for their correct localization in comparison to others. We also demonstrate that specification and maintenance of oocyte fate requires maximal Egl RNA binding activity. Even a subtle reduction in Egl's RNA binding activity completely disrupts this process. Our results show that efficient RNA localization at the earliest stages of oogenesis is required for specification of the oocyte and restriction of meiosis to a single cell.

1H, 13C and 15N resonance assignment of the YTH domain of YTHDC2

In humans, YTH (YT521-B homology) domain containing protein 2 (YTHDC2) plays a crucial role in the phase-shift from mitosis to meiosis. YTH domains bind to methylated adenosine nucleotides such as m⁶A. In a phylogenic tree, the YTH domain of YTHDC2 (YTH2) and that of the YTH containing protein YTHDC1 (YTH1) belong to the same sub-group. However, the binding affinity of m⁶A differs between these proteins. Here, we report ¹H, ¹³C and ¹⁵N resonance assignment of YTH2 and its solution structure to examine the difference of the structural architecture and the dynamic properties of YTH1 and YTH2. YTH2 adopts a β1-α1-β2-α2-β3-β4-β5-α3-β6-α4 topology, which was also observed in YTH1. However, the β4-β5 loops of YTH1 and YTH2 are distinct in length and amino acid composition. Our data revealed that, unlike in YTH1, the structure of m⁶A-binding pocket of YTH2 formed by the β4-β5 loop is stabilized by electrostatic interaction. This assignment and the structural information for YTH2 will provide the insight on the further functional research of YTHDC2.

Identification and characterization of a silent mutation in RNA binding domain of N protein coding gene from SARS-CoV-2

Objective

This study describes the occurrence of a silent mutation in the RNA binding domain of nucleocapsid phosphoprotein (N protein) coding gene from SARS-CoV-2 that may consequence to a missense mutation by onset of another single nucleotide mutation.

Results

In the DNA sequence isolated from severe acute respiratory syndrome (SARS-CoV-2) in Iran, a coding sequence for the RNA binding domain of N protein was detected. The comparison of Chinese and Iranian DNA sequences displayed that a thymine (T) was mutated to cytosine (C), so “TTG” from China was changed to “CTG” in Iran. Both DNA sequences from Iran and China have been encoded for leucine. In addition, the second T in “CTG” in the DNA or uracil (U) in “CUG” in the RNA sequences from Iran can be mutated to another C by a missense mutation resulting from thymine DNA glycosylase (TDG) of human and base excision repair mechanism to produce “CCG” encoding for proline, which consequently may increase the affinity of the RNA binding domain of N protein to viral RNA and improve the transcription rate, pathogenicity, evasion from human immunity system, spreading in the human body, and risk of human-to-human transmission rate of SARS-CoV-2.

Computational drug repurposing study of the RNA binding domain of SARS-CoV-2 nucleocapsid protein with antiviral agents

The recent outbreak of coronavirus disease (COVID‐19) in China caused by the severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) has led to worldwide human infections and deaths. The nucleocapsid (N) protein of coronaviruses (CoVs) is a multifunctional RNA binding protein necessary for viral RNA replication and transcription. Therefore, it is a potential antiviral drug target, serving multiple critical functions during the viral life cycle. This study addresses the potential to repurpose antiviral compounds approved or in development for treating human CoV induced infections against SARS‐CoV‐2 N. For this purpose, we used the docking methodology to better understand the inhibitory mechanism of this protein with the existing 34 antiviral compounds. The results of this analysis indicate that rapamycin, saracatinib, camostat, trametinib, and nafamostat were the top hit compounds with binding energy (−11.87, −10.40, −9.85, −9.45, −9.35 kcal/mol, respectively). This analysis also showed that the most common residues that interact with the compounds are Phe66, Arg68, Gly69, Tyr123, Ile131, Trp132, Val133, and Ala134. Subsequently, protein‐ligand complex stability was examined with molecular dynamics simulations for these five compounds, which showed the best binding affinity. According to the results of this study, the interaction between these compounds and crucial residues of the target protein were maintained. These results suggest that these residues are potential drug targeting sites for the SARS‐CoV‐2 N protein. This study information will contribute to the development of novel compounds for further in vitro and in vivo studies of SARS‐CoV‐2, as well as possible new drug repurposing strategies to treat COVID‐19 disease.

Crystal structure of SARS-CoV-2 nucleocapsid protein RNA binding domain reveals potential unique drug targeting sites

The outbreak of coronavirus disease (COVID-19) caused by SARS-CoV-2 virus continually lead to worldwide human infections and deaths. Currently, there is no specific viral protein-targeted therapeutics. Viral nucleocapsid protein is a potential antiviral drug target, serving multiple critical functions during the viral life cycle. However, the structural information of SARS-CoV-2 nucleocapsid protein remains unclear. Herein, we have determined the 2.7 Å crystal structure of the N-terminal RNA binding domain of SARS-CoV-2 nucleocapsid protein. Although the overall structure is similar as other reported coronavirus nucleocapsid protein N-terminal domain, the surface electrostatic potential characteristics between them are distinct. Further comparison with mild virus type HCoV-OC43 equivalent domain demonstrates a unique potential RNA binding pocket alongside the β-sheet core. Complemented by in vitro binding studies, our data provide several atomic resolution features of SARS-CoV-2 nucleocapsid protein N-terminal domain, guiding the design of novel antiviral agents specific targeting to SARS-CoV-2.

Nascent Pre-rRNA Sorting via Phase Separation Drives the Assembly of Dense Fibrillar Components in the Human Nucleolus

Fibrillar centers (FCs) and dense fibrillar components (DFCs) are essential morphologically distinct sub-regions of mammalian cell nucleoli for rDNA transcription and pre-rRNA processing. Here, we report that a human nucleolus consists of several dozen FC/DFC units, each containing 2-3 transcriptionally active rDNAs at the FC/DFC border. Pre-rRNA processing factors, such as fibrillarin (FBL), form 18-24 clusters that further assemble into the DFC surrounding the FC. Mechanistically, the 5' end of nascent 47S pre-rRNA binds co-transcriptionally to the RNA-binding domain of FBL. FBL diffuses to the DFC, where local self-association via its glycine- and arginine-rich (GAR) domain forms phase-separated clusters to immobilize FBL-interacting pre-rRNA, thus promoting directional traffic of nascent pre-rRNA while facilitating pre-rRNA processing and DFC formation. These results unveil FC/DFC ultrastructures in nucleoli and suggest a conceptual framework for considering nascent RNA sorting using multivalent interactions of their binding proteins.

Molecular insights into the specific recognition between the RNA binding domain qRRM2 of hnRNP F and G-tract RNA: A molecular dynamics study

Heterogeneous nuclear ribonucleoprotein F (hnRNP F) controls the expression of various genes through regulating the alternative splicing of pre-mRNAs in the nucleus. It uses three quasi-RNA recognition motifs (qRRMs) to recognize G-tract RNA which contains at least three consecutive guanines. The structures containing qRRMs of hnRNP F in complex with G-tract RNA have been determined by nuclear magnetic resonance (NMR) spectroscopy, shedding light on the recognition mechanism of qRRMs with G-tract RNA. However, knowledge of the recognition details is still lacking. To investigate how qRRMs specifically bind with G-tract RNA and how the mutations of any guanine to an adenine in the G-tract affect the binding, molecular dynamics simulations with binding free energy analysis were performed based on the NMR structure of qRRM2 in complex with G-tract RNA. Simulation results demonstrate that qRRM2 binds strongly with G-tract RNA, but any mutation of the G-tract leads to a drastic reduction of the binding free energy. Further comparisons of the energetic components reveal that van der Waals and non-polar interactions play essential roles in the binding between qRRM2 and G-tract RNA, but the interactions are weakened by the effect of RNA mutations. Structural and dynamical analyses indicate that when qRRM2 binds with G-tract RNA, both qRRM2 and G-tract maintain stabilized structures and dynamics; however, the stability is disrupted by the mutations of the G-tract. These results provide novel insights into the recognition mechanism of qRRM2 with G-tract RNA that are not elucidated by the NMR technique.

The expanding universe of ribonucleoproteins: of novel RNA-binding proteins and unconventional interactions

Post-transcriptional regulation of gene expression plays a critical role in almost all cellular processes. Regulation occurs mostly by RNA-binding proteins (RBPs) that recognise RNA elements and form ribonucleoproteins (RNPs) to control RNA metabolism from synthesis to decay. Recently, the repertoire of RBPs was significantly expanded owing to methodological advances such as RNA interactome capture. The newly identified RNA binders are involved in diverse biological processes and belong to a broad spectrum of protein families, many of them exhibiting enzymatic activities. This suggests the existence of an extensive crosstalk between RNA biology and other, in principle unrelated, cell functions such as intermediary metabolism. Unexpectedly, hundreds of new RBPs do not contain identifiable RNA-binding domains (RBDs), raising the question of how they interact with RNA. Despite the many functions that have been attributed to RNA, our understanding of RNPs is still mostly governed by a rather protein-centric view, leading to the idea that proteins have evolved to bind to and regulate RNA and not vice versa. However, RNPs formed by an RNA-driven interaction mechanism (RNA-determined RNPs) are abundant and offer an alternative explanation for the surprising lack of classical RBDs in many RNA-interacting proteins. Moreover, RNAs can act as scaffolds to orchestrate and organise protein networks and directly control their activity, suggesting that nucleic acids might play an important regulatory role in many cellular processes, including metabolism.

Characterization of Two Dinoflagellate Cold Shock Domain Proteins

Dinoflagellate transcriptomes contain cold shock domain proteins as the major component of the proteins annotated as transcription factors. We show here that the major family of cold shock domain proteins in the dinoflagellate Lingulodinium do not bind specific sequences, suggesting that transcriptional control is not a predominant mechanism for regulating gene expression in this group of protists.

Roughly two-thirds of the proteins annotated as transcription factors in dinoflagellate transcriptomes are cold shock domain-containing proteins (CSPs), an uncommon condition in eukaryotic organisms. However, no functional analysis has ever been reported for a dinoflagellate CSP, and so it is not known if they do in fact act as transcription factors. We describe here some of the properties of two CSPs from the dinoflagellate Lingulodinium polyedrum, LpCSP1 and LpCSP2, which contain a glycine-rich C-terminal domain and an N-terminal cold shock domain phylogenetically related to those in bacteria. However, neither of the two LpCSPs act like the bacterial CSP, since they do not functionally complement the Escherichia coli quadruple cold shock domain protein mutant BX04, and cold shock does not induce LpCSP1 and LpCSP2 to detectable levels, based on two-dimensional gel electrophoresis. Both CSPs bind to RNA and single-stranded DNA in a nonspecific manner in electrophoretic mobility shift assays, and both proteins also bind double-stranded DNA nonspecifically, albeit more weakly. These CSPs are thus unlikely to act alone as sequence-specific transcription factors.

IMPORTANCE Dinoflagellate transcriptomes contain cold shock domain proteins as the major component of the proteins annotated as transcription factors. We show here that the major family of cold shock domain proteins in the dinoflagellate Lingulodinium do not bind specific sequences, suggesting that transcriptional control is not a predominant mechanism for regulating gene expression in this group of protists.

Structure-specific nucleic acid recognition by L-motifs and their diverse roles in expression and regulation of the genome

The high-mobility group (HMG) domain containing proteins regulate transcription, DNA replication and recombination. They adopt L-shaped folds and are structure-specific DNA binding motifs. Here, I define the L-motif super-family that consists of DNA-binding HMG-box proteins and the L-motif of the histone mRNA binding domain of stem-loop binding protein (SLBP). The SLBP L-motif and HMG-box domains adopt similar L-shaped folds with three α-helices and two or three small hydrophobic cores that stabilize the overall fold, but have very different and distinct modes of nucleic acid recognition. A comparison of the structure, dynamics, protein-protein and nucleic acid interactions, and regulation by PTMs of the SLBP and the HMG-box L-motifs reveals the versatile and diverse modes by which L-motifs utilize their surfaces for structure-specific recognition of nucleic acids to regulate gene expression.