RNA Hairpin Research Articles

OS10.7.A Activation of the CCR8-ACP5 axis by human microglia/macrophage derived CCL18 promotes glioma growth

Abstract Background Glioblastoma multiforme is a highly malignant primary brain tumor with an average survival of 14 months and very limited therapeutic options. Glioma associated microglia and macrophages (GAMs) foster tumor growth by releasing several cytokines, which have only partly been identified. Here, we studied the chemokine (C-C motif) ligand 18 (CCL18), a chemokine which is only expressed in human, but not rodent GAMs, in a novel ex-vivo brain slice model including transplantation of human induced pluripotent stem cells (iPSC) derived human microglia (iMGL) and human glioma cells in to murine brain slices, which had been depleted of intrinsic murine microglia before. Material and Methods After establishing the humanized ex-vivo brain slice model, we performed immunohistochemical analysis (IHC) of growth and invasiveness, qrtPCR on glioma cells isolated by magnetic-activated cell sorting (MACS), functional assays measuring invasiveness, proliferation, migration and colony formation of glioma cells in vitro and in slice experiments. Corresponding studies on tumor growth and invasiveness were performed after treatment with a CCL18 neutralizing antibody, a CCR8 neutralizing antibodies and knockdown of CCR8, ACP5 (Acid Phosphatase 5) and PITPNM3 with small interfering RNA (siRNA) and short hairpin RNA (shRNA). QrtPCR, IHC and Westernblot analysis were performed on primary glioma specimens. We also conducted bioinformatic analyses, based on the TCGA GBM, GLIOVIS and GEPIA databases. Results We observed that CCL18 was highly expressed in GAMs, whereas CCR8 was only expressed in glioma cells. We identified the chemokine (C-C motif) receptor 8 (CCR8) as a functional receptor for CCL18 and ACP5 as an important down-stream signaling component in glioma cells. Activation of the CCL18/CCR8/ACP5 signaling pathway in human glioblastoma was associated with enhanced tumor growth and invasiveness. Conclusion GAMs derived CCL18 promoted glioma growth by activation of the CCR8/ACP5 axis in human glioma cells and therefore is a potential therapeutic target.

Extension of the CHARMM Classical Drude Polarizable Force Field to N- and O-Linked Glycopeptides and Glycoproteins.

Molecular dynamic simulations are an effective tool to study complex molecular systems and are contingent upon the availability of an accurate and reliable molecular mechanics force field. The Drude polarizable force field, which allows for the explicit treatment of electronic polarization in a computationally efficient fashion, has been shown to reproduce experimental properties that were difficult or impossible to reproduce with the CHARMM additive force field, including peptide folding cooperativity, RNA hairpin structures, and DNA base flipping. Glycoproteins are essential components of glycoconjugate vaccines, antibodies, and many pharmaceutically important molecules, and an accurate polarizable force field that includes compatibility between the protein and carbohydrate aspect of the force field is essential to study these types of systems. In this work, we present an extension of the Drude polarizable force field to glycoproteins, including both N- and O-linked species. Parameter optimization focused on the dihedral terms using a reweighting protocol targeting NMR solution J-coupling data for model glycopeptides. Validation of the model include eight model glycopeptides and four glycoproteins with multiple N- and O-linked glycosylations. The new glycoprotein carbohydrate force field can be used in conjunction with the remainder of Drude polarizable force field through a variety of MD simulation programs including GROMACS, OPENMM, NAMD, and CHARMM and may be accessed through the Drude Prepper module in the CHARMM-GUI.

In transcription antitermination by Qλ, NusA induces refolding of Qλ to form a nozzle that extends the RNA polymerase RNA-exit channel

Lambdoid bacteriophage Q proteins are transcription antipausing and antitermination factors that enable RNA polymerase (RNAP) to read through pause and termination sites. Q proteins load onto RNAP engaged in promoter-proximal pausing at a Q binding element (QBE) and adjacent sigma-dependent pause element to yield a Q-loading complex, and they translocate with RNAP as a pausing-deficient, termination-deficient Q-loaded complex. In previous work, we showed that the Q protein of bacteriophage 21 (Q21) functions by forming a nozzle that narrows and extends the RNAP RNA-exit channel, preventing formation of pause and termination RNA hairpins. Here, we report atomic structures of four states on the pathway of antitermination by the Q protein of bacteriophage λ (Qλ), a Q protein that shows no sequence similarity to Q21 and that, unlike Q21, requires the transcription elongation factor NusA for efficient antipausing and antitermination. We report structures of Qλ, the Qλ-QBE complex, the NusA-free pre-engaged Qλ-loading complex, and the NusA-containing engaged Qλ-loading complex. The results show that Qλ, like Q21, forms a nozzle that narrows and extends the RNAP RNA-exit channel, preventing formation of RNA hairpins. However, the results show that Qλ has no three-dimensional structural similarity to Q21, employs a different mechanism of QBE recognition than Q21, and employs a more complex process for loading onto RNAP than Q21, involving recruitment of Qλ to form a pre-engaged loading complex, followed by NusA-facilitated refolding of Qλ to form an engaged loading complex. The results establish that Qλ and Q21 are not structural homologs and are solely functional analogs.

Synthesis of a Peptidoyl RNA Hairpin via a Combination of Solid-Phase and Template-Directed Chain Assembly.

Peptidoyl RNAs are the products of ribosome‐free, single‐nucleotide translation. They contain a peptide in the backbone of the oligoribonucleotide and are interesting from a synthetic and a bioorganic point of view. A synthesis of a stabilized version of peptidoyl RNA, with an amide bond between the C‐terminus of a peptide and a 3′‐amino‐2′,3′‐dideoxynucleoside in the RNA chain was developed. The preferred synthetic route used an N‐Teoc‐protected aminonucleoside support and involved a solution‐phase coupling of the amino‐terminal oligonucleotide to a dipeptido dinucleotide. Exploratory UV‐melting and NMR analysis of the hairpin 5′‐UUGGCGAAAGCdC‐LeuLeu‐AA‐3′ indicated that the peptide‐linked RNA segments do not fold in a cooperative fashion. The synthetic access to doubly RNA‐linked peptides on a scale sufficient for structural biology opens the door to the exploration of their structural and biochemical properties.

Identification of an ERN1 target site within EGFP mRNA.

EGFP (enhanced green fluorescent protein) is one of the most common tools used in life sciences, including research focusing on proteostasis. Here we report that ERN1 (endoplasmic reticulum to nucleus signaling 1), which is upregulated by UPR (unfolded protein response), targets an RNA hairpin loop motif in EGFP mRNA. A silent mutation introduced into EGFP mRNA abolished the ERN1‐dependent mRNA decay. Therefore, experiments that employ EGFP as a reporter gene in studies that involve upregulation of the UPR pathway should be interpreted carefully, and a mutant devoid of the ERN1 target motif may be more suitable for such studies.

Structural basis of AlpA-dependent transcription antitermination.

AlpA positively regulates a programmed cell death pathway linked to the virulence of Pseudomonas aeruginosa by recognizing an AlpA binding element within the promoter, then binding RNA polymerase directly and allowing it to bypass an intrinsic terminator positioned downstream. Here, we report the single-particle cryo-electron microscopy structures of both an AlpA-loading complex and an AlpA-loaded complex. These structures indicate that the C-terminal helix-turn-helix motif of AlpA binds to the AlpA binding element and that the N-terminal segment of AlpA forms a narrow ring inside the RNA exit channel. AlpA was also revealed to render RNAP resistant to termination signals by prohibiting RNA hairpin formation in the RNA exit channel. Structural analysis predicted that AlpA, 21Q, λQ and 82Q share the same mechanism of transcription antitermination.

Analysis of luciferase dsRNA production during baculovirus infection of Hi5 cells: RNA hairpins expressed by very late promoters do not trigger gene silencing.

The baculovirus expression vector system (BEVS) has become an important platform for the expression of recombinant proteins and is especially useful for the production of large protein complexes such as virus-like particles (VLPs). An important application for VLPs is their use as vehicles for targeted delivery of drugs or toxins which requires the development of methods for efficient loading with the intended cargo. Our research intends to employ the BEVS for the production of VLPs for the delivery of insecticidal dsRNA molecules to targeted insect pests (as "dsRNA-VLPs"). A convenient strategy would be the co-expression of long dsRNAs with viral capsid proteins and their simultaneous encapsulation during VLP assembly but the capacity of the BEVS for the production of long dsRNA has not been assessed so far. In this study, the efficiency of production of long RNA hairpins targeting the luciferase gene ("dsLuc") by the polyhedrin promoter during baculovirus infection was evaluated. However, RNAi reporter assays could not detect significant amounts of dsLuc in Hi5 cells infected with recombinant baculovirus, even in the presence of co-expressed dsRNA-binding protein B2-GFP or the employment of the MS2-MCP system. Nevertheless, dot blot analyses using anti-dsRNA antibody revealed that baculovirus-mediated expression of B2-GFP resulted in significant increases in dsRNA levels in infected cells that may correspond to hybridized complementary viral transcripts. Using B2-GFP as a genetically encoded sensor, dsRNA foci were detected in the nuclei that partially co-localized with DAPI staining, consistent with their localization at the virogenic stroma. Co-localization experiments with the baculovirus proteins vp39, Ac93, ODV-E25 and gp64 indicated limited overlap between B2-GFP and the ring zone compartment where assembly of nucleocapsids and virions occurs. Stability experiments showed that exogenous dsRNA is resistant to degradation in extracts of non-infected and infected Hi5 cells and it is proposed that strong unwinding activity at the virogenic stroma in the infected nuclei may neutralize the annealing of complementary RNA strands and block the production of long dsRNAs. Because the strong stability of exogenous dsRNA, transfection can be explored as an alternative method for delivery of cargo for dsRNA-VLPs during their assembly in baculovirus-infected Hi5 cells.

Identification of a guanine-specific pocket in the protein N of SARS-CoV-2

The SARS-CoV-2 nucleocapsid protein (N) is responsible for RNA binding. Here we report the crystal structure of the C-terminal domain (NCTD) in open and closed conformations and in complex with guanine triphosphate, GTP. The crystal structure and biochemical studies reveal a specific interaction between the guanine, a nucleotide enriched in the packaging signals regions of coronaviruses, and a highly conserved tryptophan residue (W330). In addition, EMSA assays with SARS-CoV-2 derived RNA hairpin loops from a putative viral packaging sequence showed the preference interaction of the N-CTD to RNA oligonucleotides containing G and the loss of the specificity in the mutant W330A. Here we propose that this interaction may facilitate the viral assembly process. In summary, we have identified a specific guanine-binding pocket in the N protein that may be used to design viral assembly inhibitors.

Investigating the role of RNA structures in transcriptional pausing using in vitro assays and in silico analyses

ABSTRACT Transcriptional pausing occurs across the bacterial genome but the importance of this mechanism is still poorly understood. Only few pauses were observed during the previous decades, leaving an important gap in understanding transcription mechanisms. Using the well-known Escherichia coli hisL and trpL pause sites as models, we describe here the relation of pause sites with upstream RNA structures suspected to stabilize pausing. We find that the transcription factor NusA influences the pause half-life at leuL, pheL and thrL pause sites. Using a mutagenesis approach, we observe that transcriptional pausing is affected in all tested pause sites, suggesting that the upstream RNA sequence is important for transcriptional pausing. Compensatory mutations assessing the presence of RNA hairpins did not yield clear conclusions, indicating that complex RNA structures or transcriptional features may be playing a role in pausing. Moreover, using a bioinformatic approach, we explored the relation between a DNA consensus sequence important for pausing and putative hairpins among thousands of pause sites in E. coli. We identified 2125 sites presenting hairpin-dependent transcriptional pausing without consensus sequence, suggesting that this mechanism is widespread across E. coli. This study paves the way to understand the role of RNA structures in transcriptional pausing.

Delta-like factor 1 negatively regulates angiogenesis as a target gene of miR-126-5p after indirect revascularization surgery in patients with moyamoya disease

Promoting endothelial cell (EC) proliferation and angiogenesis after indirect revascularization surgery is crucial for decreasing the stroke rate in moyamoya patients. However, the role of delta-like factor 1 (DLK1) in regulating EC proliferation in chronically ischaemic brains and the specific mechanisms remain unclear. Therefore, we compared the expression levels of DLK1 in the dura mater (DM) tissues of patients with moyamoya disease and patients with aneurysms, and dual luciferase reporter and RNA binding protein immunoprecipitation assays were conducted to determine whether DLK1 is a target gene of miR-126-5p. The effect of DLK1 on EC proliferation and the interaction between DLK1 and miR-126-5p were explored in vitro. Then, we established an animal model of two-vessel occlusion together with encephalo-myo-synangiosis (2VO+EMS). The temporalis muscles (TMs) of the animals were transfected with DLK1 lentiviral RNA (Lv-DLK1) and short hairpin RNA (sh-DLK1) to compare the DLK1 expression, angiogenesis (CD31 count), and numbers of vacuoles and impaired tight junctions in the ECs of TM-covered chronically ischaemic brains as well as the cognitive improvement in each group. DLK1 expression was lower in the DM tissues of moyamoya patients than in those of the aneurysm patients, and DLK1 was identified as a target gene of miR-126-5p. In vitro, DLK1 inhibited EC proliferation, migration and angiogenesis and exerted effects opposite those of miR-126-5p. In 2VO+EMS rats, compared to the control transfection, the TM transfection of Lv-DLK1 induced significantly higher DLK1 expression and worse angiogenesis in the TM-covered ischaemic brain as well as less extensive cognitive improvement, while the transfection of sh-DLK1 into the TM had the opposite effects. In summary, DLK1, a target gene of miR-126-5p, negatively regulates angiogenesis in chronically ischaemic brains and is expected to be a new target to improve the efficacy of indirect revascularization surgery and the prognosis of moyamoya patients.

Alternative splicing of Apoptosis Stimulating Protein of TP53-2 (ASPP2) results in an oncogenic isoform promoting migration and therapy resistance in soft tissue sarcoma (STS)

BackgroundMetastatic soft tissue sarcoma (STS) are a heterogeneous group of malignancies which are not curable with chemotherapy alone. Therefore, understanding the molecular mechanisms of sarcomagenesis and therapy resistance remains a critical clinical need. ASPP2 is a tumor suppressor, that functions through both p53-dependent and p53-independent mechanisms. We recently described a dominant-negative ASPP2 isoform (ASPP2κ), that is overexpressed in human leukemias to promote therapy resistance. However, ASPP2κ has never been studied in STS. Materials and methodsExpression of ASPP2κ was quantified in human rhabdomyosarcoma tumors using immunohistochemistry and qRT-PCR from formalin-fixed paraffin-embedded (FFPE) and snap-frozen tissue. To study the functional role of ASPP2κ in rhabdomyosarcoma, isogenic cell lines were generated by lentiviral transduction with short RNA hairpins to silence ASPP2κ expression. These engineered cell lines were used to assess the consequences of ASPP2κ silencing on cellular proliferation, migration and sensitivity to damage-induced apoptosis. Statistical analyses were performed using Student’s t-test and 2-way ANOVA.ResultsWe found elevated ASPP2κ mRNA in different soft tissue sarcoma cell lines, representing five different sarcoma sub-entities. We found that ASSP2κ mRNA expression levels were induced in these cell lines by cell-stress. Importantly, we found that the median ASPP2κ expression level was higher in human rhabdomyosarcoma in comparison to a pool of tumor-free tissue. Moreover, ASPP2κ levels were elevated in patient tumor samples versus adjacent tumor-free tissue within individual patients. Using isogenic cell line models with silenced ASPP2κ expression, we found that suppression of ASPP2κ enhanced chemotherapy-induced apoptosis and attenuated cellular proliferation.ConclusionDetection of oncogenic ASPP2κ in human sarcoma provides new insights into sarcoma tumor biology. Our data supports the notion that ASPP2κ promotes sarcomagenesis and resistance to therapy. These observations provide the rationale for further evaluation of ASPP2κ as an oncogenic driver as well as a prognostic tool and potential therapeutic target in STS.

ProtNAff: protein-bound Nucleic Acid filters and fragment libraries.

Atomistic models of nucleic acids (NA) fragments can be used to model the 3D structures of specific protein-NA interactions and address the problem of great NA flexibility, especially in their single-stranded regions. One way to obtain relevant NA fragments is to extract them from existing 3D structures corresponding to the targeted context (e.g. specific 2D structures, protein families, sequences) and to learn from them. Several databases exist for specific NA 3D motifs, especially in RNA, but none can handle the variety of possible contexts. This article presents protNAff (protein-bound Nucleic Acids filters and fragments), a new pipeline for the conception of searchable databases on the 2D and 3D structures of protein-bound NA, the selection of context-specific (regions of) NA structures by combinations of filters, and the creation of context-specific NA fragment libraries. The strength of this pipeline is its modularity, allowing users to adapt it to many specific modeling problems. As examples, the pipeline is applied to the quantitative analysis of (i) the sequence-specificity of trinucleotide conformations, (ii) the conformational diversity of RNA at several levels of resolution, (iii) the effect of protein binding on RNA local conformations and (iv) the protein-binding propensity of RNA hairpin loops of various lengths. The source code is freely available for download at URL https://github.com/isaureCdB/protNAff. The database and the trinucleotide fragment library are downloadable at URL https://zenodo.org/record/6483823#.YmbVhFxByV4. Supplementary data are available at Bioinformatics online.

Solution Structure of NPSL2, A Regulatory Element in the oncomiR-1 RNA

The miR-17 ∼ 92a polycistron, also known as oncomiR-1, is commonly overexpressed in multiple cancers and has several oncogenic properties. OncomiR-1 encodes six constituent microRNAs (miRs), each enzymatically processed with different efficiencies. However, the structural mechanism that regulates this differential processing remains unclear. Chemical probing of oncomiR-1 revealed that the Drosha cleavage sites of pri-miR-92a are sequestered in a four-way junction. NPSL2, an independent stem loop element, is positioned just upstream of pri-miR-92a and sequesters a crucial part of the sequence that constitutes the basal helix of pri-miR-92a. Disruption of the NPSL2 hairpin structure could promote the formation of a pri-miR-92a structure that is primed for processing by Drosha. Thus, NPSL2 is predicted to function as a structural switch, regulating pri-miR-92a processing. Here, we determined the solution structure of NPSL2 using solution NMR spectroscopy. This is the first high-resolution structure of an oncomiR-1 element. NPSL2 adopts a hairpin structure with a large, but highly structured, apical and internal loops. The 10-bp apical loop contains a pH-sensitive A+·C mismatch. Additionally, several adenosines within the apical and internal loops have elevated pKa values. The protonation of these adenosines can stabilize the NPSL2 structure through electrostatic interactions. Our study provides fundamental insights into the secondary and tertiary structure of an important RNA hairpin proposed to regulate miR biogenesis.

Engineered anti-viral shRNAs are more effective than lhRNAs in transgenic Aedes aegypti

The Aedes aegypti mosquito is a major vector of chikungunya virus (CHIKV), which has no licensed vaccine. Engineered mosquitoes expressing long RNA hairpins (lhRNA) or small RNAs against selected arboviruses have been developed to limit virus replication, but their silencing efficiency has not been compared. We developed lhRNA and short hairpin RNA (shRNA) arrays against CHIKV non-structural protein nsP2. We used a Tet-response element (TRE) and a tTA transactivator to control expression of lhRNAs (TRE-lhRNA) and shRNAs (TRE-shRNA). Constitutive expression in Aag2 cells was assessed with a PUb-tTA driverand midgut specific expression in transgenic mosquitoes was assessed using Carboxypeptidase A (AeCPA)-tTA as driver. In vitrointerference ability was determined with a CHIKV split replication system, and a synthetic luciferase reporter with mRNA containing the targeted CHIKV sequence (CHIK-FF). In vivo interference was tested by inserting the TRE-lhRNA and TRE-shRNA constructs into a AeCPA-tTA line so both constructs expressed from the same locus, and a TRE reporter line that expressed an AmCyan reporter with an N terminal CHIKV fusion (CHIK-AmC). In Aag2 cells, shRNAs were more effective than lhRNA in silencing both the CHIKV split replication system (99.7%, S. D. ± 0.47 and 72.8%, S. D. ± 9.50, respectively, P<0.001) and CHIK-FF (98.8%, S. D. ± 0.71 and 50.9%, S. D. ± 7.92, respectively, P<0.05). Similar results were observed in transgenic mosquitoes when comparing AmC expression in the midgut. This study demonstrates that, in mosquitoes, effectively chosen shRNAs can induce greater interference of the desired viral target than the corresponding lhRNA.

Site-specific photolabile roadblocks for the study of transcription elongation in biologically complex systems

Transcriptional pausing is crucial for the timely expression of genetic information. Biochemical methods quantify the half-life of paused RNA polymerase (RNAP) by monitoring restarting complexes across time. However, this approach may produce apparent half-lives that are longer than true pause escape rates in biological contexts where multiple consecutive pause sites are present. We show here that the 6-nitropiperonyloxymethyl (NPOM) photolabile group provides an approach to monitor transcriptional pausing in biological systems containing multiple pause sites. We validate our approach using the well-studied his pause and show that an upstream RNA sequence modulates the pause half-life. NPOM was also used to study a transcriptional region within the Escherichia coli thiC riboswitch containing multiple consecutive pause sites. We find that an RNA hairpin structure located upstream to the region affects the half-life of the 5′ most proximal pause site—but not of the 3′ pause site—in contrast to results obtained using conventional approaches not preventing asynchronous transcription. Our results show that NPOM is a powerful tool to study transcription elongation dynamics within biologically complex systems.

Tetraphenylethene Derivatives Modulate the RNA Hairpin-G-Quadruplex Conformational Equilibria in Proto-oncogenes.

RNA G-quadruplex (GQs) sequences in 5'-UTRs of certain proto-oncogenes co-localize with hairpin (Hp) forming sequences resulting in intramolecular Hp-GQ conformational equilibria, which is suggested to regulate cancer development and progression. Thus, regulation of Hp-GQ equilibria with small molecules is an attractive but less explored therapeutic approach. Herein, two tetraphenylethene (TPE) derivatives, TPE-Py and TPE-MePy, were synthesized and their effect on Hp-GQ equilibrium was explored. FRET, CD and molecular docking experiments suggest that cationic TPE-MePy shifts the Hp-GQ equilibrium significantly towards the GQ conformer mainly through π-π stacking and van der Waals interactions. In the presence of TPE-MePy, the observed rate constant values for first and second folding steps were increased up to 14.6 and 2.6-fold, respectively. The FRET melting assay showed a strong stabilizing ability of TPE-MePy (ΔTm=4.36 °C). Notably, the unmethylated derivative TPE-Py did not alter the Hp-GQ equilibrium. Subsequently, luciferase assay analysis demonstrated that the TPE-MePy derivatives suppressed the translation efficiency by ∼5.7-fold by shifting the Hp-GQ equilibrium toward GQ conformers in the 5'-UTR of TRF2. Our data suggests that HpGQ equilibria could be selectively targeted with small molecules to modulate translation for therapy.

Accelerating the Ensemble Convergence of RNA Hairpin Simulations with a Replica Exchange Structure Reservoir.

RNA is a key participant in many biological processes, but studies of RNA using computer simulations lag behind those of proteins, largely due to less-developed force fields and the slow dynamics of RNA. Generating converged RNA ensembles for force field development and other studies remains a challenge. In this study, we explore the ability of replica exchange molecular dynamics to obtain well-converged conformational ensembles for two RNA hairpin systems in an implicit solvent. Even for these small model systems, standard REMD remains computationally costly, but coupling to a pre-generated structure library using the reservoir REMD approach provides a dramatic acceleration of ensemble convergence for both model systems. Such precise ensembles could facilitate RNA force field development and validation and applications of simulation to more complex RNA systems. The advantages and remaining challenges of applying R-REMD to RNA are investigated in detail.

Nucleolin’s Role in the Biogenesis of Pri‐miRNA 16 by the Microprocessor Complex: An In silico Approach

MicroRNAs (miRNAs) are short, non‐coding strands of RNA found in eukaryotes. They play a role in gene silencing through inhibiting translation or degrading mRNA. Primary miRNA (pri‐miRNA) processing is the initial step in producing a mature microRNA molecule. The microprocessor complex (MP) proteins, DROSHA and DGCR8, cleaves the pri‐miRNA molecule to a shorter hairpin structure. This cleaved RNA hairpin, precursor miRNA, is sent to the cytoplasm for further processing into the mature miRNA. Previous research has uncovered that the expression and cellular localization of the stress‐responsive protein, nucleolin (NCL) affects the processing of several miRNAs, including microRNA 16 (miR16). Nucleolin is an RNA binding protein and is primarily located in the nucleolus but can also be localized in various cellular compartments. When NCL is present in the nucleus, the MP complex works more efficiently in cleaving pri‐miR 16‐1 to its pre‐miRNA form resulting in an upregulation of miR 16‐1. It has been proposed that nucleolin binds to pri‐miR 16‐1 and stabilizes its confirmation for cleavage by the MP complex. However, the mechanistic and structural details on how nucleolin interacts with pri‐miRNA 16‐1 or the MP are still unclear. It is also unknown if NCL is similarly involved in miR 16‐2 biogenesis. NCL and miRNA 16 are also intricately involved in the regulation of the expression of bcl2, an anti‐apoptotic and oncogenic protein. When NCL is in the cytoplasm, it stabilizes the 3’ untranslated region of bcl2 mRNA leading to an increased expression of the protein. Unlike nucleolin, miR 16 represses the translation of bcl2 mRNA. At high levels of miR 16 when NCL is localized to the nucleus, bcl2 protein expression is low. When cellular stress occurs, NCL returns to the cytoplasm to stabilize bcl2 and this also results in less efficiency in pri‐miR 16 processing, resulting in an overall increase in bcl2 protein. It is evident that NCL balances the levels of bcl2 and miR16 through its cellular localization. Since the molecular mechanism underlying NCL’s role in the MP‐based processing of pri‐miR16‐1 remains largely unknown, in this study, we used computational approaches including structural alignments, protein‐protein, and protein‐RNA docking analyses to determine how pri‐miR 16‐1, nucleolin, and the MP complex interact. In addition, we have investigated if pri‐miR 16‐2 interacts with the NCL and the MP complex in a similar fashion. We show in our predictive models that the interaction interface of pri‐miR 16‐1 and pri‐miR 16‐2 with nucleolin and the MP complex nucleolin is distinct. For pri‐miR 16‐1, the interaction interface involves NCL and DGCR8 B and C, but in the case of pri‐miR 16‐2, it involves NCL and DROSHA. Additionally, in the presence of nucleolin, the hydrogen bond interactions are reduced between pri‐miR 16‐2 and the MP complex. Although the number of interactions between the different components of the MP complex is the same in both scenarios, the specific residues that form the interaction interface differ, suggesting distinct mechanisms. Since the levels of miR16 expression can lead to apoptosis of cancerous cells and prevent cell proliferation, the results of our study provide important clues in understanding how the NCL‐miR‐MP complex can be targeted in drug design for cancer therapy.

The Role of Vts1/Smaug in Adaptive Heritable Gene Expression

The ancient RNA Binding Protein (RBP) Vts1 has been identified to self‐assemble along with the Smaug regulator to produce enhancing non‐amyloid prions in mRNA decay. Seen in Drosophila, the sterile‐alpha‐motif (SAM) domain of Vts and the Smaug regulator acts as an embryonic regulator for maternal transcript degradation through cytoplasmic deadenylase in progeny development. Specifically, of Saccharomyces cerevisiae, the Vts/Smaug regulator forms condensates that self‐template to create prions functioning towards epigenetic heritable inheritance. The RBP Vts1 SAM domain recognizes in RNA hairpins Smaug recognition elements, SREs, initiating the process of degradation and self‐assembly. The structure of RBPs are composed of an ordered RNA binding domain, and an unventured intrinsically disordered region (IDR). At the quaternary level structure, the Vts1 protein matches as a hexamer of 489 kDa dominant peak with an oligomer of more than eight monomers, and the IDR is around 490 kDa as well while the RBD‐Vts1 is around 75 kDa. The condensed version of Vts1 fuse into structures of 10 micrometers due to the IDR driving self‐association to create oligomers. The RNA binding SAM domain and IDR have induced transcription degradation and post‐transcriptional regulation in a new mechanism, prompting adaptive gene expression for developmental decisions in organisms over evolutionary time periods. Though many RBP prion condensates are not heritable and inactivate protein function, the SMAUG+ prions of the condensed Vts1/Smaug protein in S. Cerevisiae instead respond to environmental stress and produce new emerging phenotypes that are actually heritable over a large biological time‐scale. Furthermore, SMAUG+ conversion is reversible and an enhancer of RBP Vts1 function, producing adaptive gene expression systems in fluctuating environments. In various yeast models, the SMAUG+ allows yeast to decide nutrient replenishment after starvation which creates a new adaptive advantage for organisms. The new cellular decision making regulated by the prion is in two forms of rapid mitosis or meiosis creating a non‐stress state; SMAUG+ has created an adaptive memory that can influence developmental and cellular decisions in addition to plant growth in limited nutrients. Such advancements in heritable prion inheritance can uncover the function of intrinsically disordered sequences. The Walton SMART team has designed a 3D model with MSOE Center for Biomolecular Modeling to investigate the SAM domain of Vts1 further and examine the role of SMAUG+ prions in adaptive heritable gene expression.

Thermodynamics of unfolding mechanisms of mouse mammary tumor virus pseudoknot from a coarse-grained loop-entropy model.

Pseudoknotted RNA molecules play important biological roles that depend on their folded structure. To understand the underlying principles that determine their thermodynamics and folding/unfolding mechanisms, we carried out a study on a variant of the mouse mammary tumor virus pseudoknotted RNA (VPK), a widely studied model system for RNA pseudoknots. Our method is based on a coarse-grained discrete-state model and the algorithm of PK3D (pseudoknot structure predictor in three-dimensional space), with RNA loops explicitly constructed and their conformational entropic effects incorporated. Our loop entropy calculations are validated by accurately capturing previously measured melting temperatures of RNA hairpins with varying loop lengths. For each of the hairpins that constitutes the VPK, we identified alternative conformations that are more stable than the hairpin structures at low temperatures and predicted their populations at different temperatures. Our predictions were validated by thermodynamic experiments on these hairpins. We further computed the heat capacity profiles of VPK, which are in excellent agreement with available experimental data. Notably, our model provides detailed information on the unfolding mechanisms of pseudoknotted RNA. Analysis of the distribution of base-pairing probability of VPK reveals a cooperative unfolding mechanism instead of a simple sequential unfolding of first one stem and then the other. Specifically, we find a simultaneous "loosening" of both stems as the temperature is raised, whereby both stems become partially melted and co-exist during the unfolding process.