Assembly and analysis of eukaryotic Argonaute-RNA complexes in microRNA-target recognition

Frontiers and Challenges of Computing ncRNAs Biogenesis, Function and Modulation

Non-coding RNAs (ncRNAs), generated from nonprotein coding DNA sequences, constitute 98-99% of the human genome. Non-coding RNAs encompass diverse functional classes, including microRNAs, small interfering RNAs, PIWI-interacting RNAs, small nuclear RNAs, small nucleolar RNAs, and long non-coding RNAs. With critical involvement in gene expression and regulation across various biological and physiopathological contexts, such as neuronal disorders, immune responses, cardiovascular diseases, and cancer, non-coding RNAs are emerging as disease biomarkers and therapeutic targets. In this review, after providing an overview of non-coding RNAs' role in cell homeostasis, we illustrate the potential and the challenges of state-of-the-art computational methods exploited to study non-coding RNAs biogenesis, function, and modulation. This can be done by directly targeting them with small molecules or by altering their expression by targeting the cellular engines underlying their biosynthesis. Drawing from applications, also taken from our work, we showcase the significance and role of computer simulations in uncovering fundamental facets of ncRNA mechanisms and modulation. This information may set the basis to advance gene modulation tools and therapeutic strategies to address unmet medical needs.

A spontaneous thermo-sensitive female sterility mutation in rice enables fully mechanized hybrid breeding

Male sterility enables hybrid crop breeding to increase yields and has been extensively studied. But thermo-sensitive female sterility, which is an ideal property that may enable full mechanization in hybrid rice breeding, has rarely been investigated due to the absence of such germplasm. Here we identify the spontaneous thermo-sensitive female sterility 1 (tfs1) mutation that confers complete sterility under regular/high temperature and partial fertility under low temperature as a point mutation in ARGONAUTE7 (AGO7). AGO7 associates with miR390 to form an RNA-Induced Silencing Complex (RISC), which triggers the biogenesis of small interfering RNAs (siRNAs) from TRANS-ACTING3 (TAS3) loci by recruiting SUPPRESSOR OF GENE SILENCING (SGS3) and RNA-DEPENDENT RNA POLYMERASE6 (RDR6) to TAS3 transcripts. These siRNAs are known as tasiR-ARFs as they act in trans to repress auxin response factor genes. The mutant TFS1 (mTFS1) protein is compromised in its ability to load the miR390/miR390* duplex and eject miR390* during RISC formation. Furthermore, tasiR-ARF levels are reduced in tfs1 due to the deficiency in RDR6 but not SGS3 recruitment by mTFS1 RISC under regular/high temperature, while low temperature partially restores mTFS1 function in RDR6 recruitment and tasiR-ARF biogenesis. A miR390 mutant also exhibits female sterility, suggesting that female fertility is controlled by the miR390-AGO7 module. Notably, the tfs1 allele introduced into various elite rice cultivars endows thermo-sensitive female sterility. Moreover, field trials confirm the utility of tfs1 as a restorer line in fully mechanized hybrid rice breeding.

Different N-Glycosylation Sites Reduce the Activity of Recombinant DSPAα2

Bat plasminogen activators α2 (DSPAα2) has extremely high medicinal value as a powerful natural thrombolytic protein. However, wild-type DSPAα2 has two N-glycosylation sites (N185 and N398) and its non-human classes of high-mannose-type N-glycans may cause immune responses in vivo. By mutating the N-glycosylation sites, we aimed to study the effect of its N-glycan chain on plasminogen activation, fibrin sensitivity, and to observe the physicochemical properties of DSPAα2. A logical structure design was performed in this study. Four single mutants and one double mutant were constructed and expressed in Pichia pastoris. When the N398 site was eliminated, the plasminogen activator in the mutants had their activities reduced to ~40%. When the N185 site was inactivated, there was a weak decrease in the plasminogen activation of its mutant, while the fibrin sensitivity significantly decreased by ~10-fold. Neither N-glycosylation nor deglycosylation mutations changed the pH resistance or heat resistance of DSPAα2. This study confirms that N-glycosylation affects the biochemical function of DSPAα2, which provides a reference for subsequent applications of DSPAα2.

Structural Modeling of the SARS-CoV-2 Spike/Human ACE2 Complex Interface can Identify High-Affinity Variants Associated with Increased Transmissibility

The COVID-19 pandemic has triggered concerns about the emergence of more infectious and pathogenic viral strains. As a public health measure, efficient screening methods are needed to determine the functional effects of new sequence variants. Here we show that structural modeling of SARS-CoV-2 Spike protein binding to the human ACE2 receptor, the first step in host-cell entry, predicts many novel variant combinations with enhanced binding affinities. By focusing on natural variants at the Spike-hACE2 interface and assessing over 700 mutant complexes, our analysis reveals that high-affinity Spike mutations (including N440K, S443A, G476S, E484R, G502P) tend to cluster near known human ACE2 recognition sites (K31 and K353). These Spike regions are structurally flexible, allowing certain mutations to optimize interface interaction energies. Although most human ACE2 variants tend to weaken binding affinity, they can interact with Spike mutations to generate high-affinity double mutant complexes, suggesting variation in individual susceptibility to infection. Applying structural analysis to highly transmissible variants, we find that circulating point mutations S477N, E484K and N501Y form high-affinity complexes (~40% more than wild-type). By combining predicted affinities and available antibody escape data, we show that fast-spreading viral variants exploit combinatorial mutations possessing both enhanced affinity and antibody resistance, including S477N/E484K, E484K/N501Y and K417T/E484K/N501Y. Thus, three-dimensional modeling of the Spike/hACE2 complex predicts changes in structure and binding affinity that correlate with transmissibility and therefore can help inform future intervention strategies.

Germline AGO2 mutations impair RNA interference and human neurological development

ARGONAUTE-2 and associated miRNAs form the RNA-induced silencing complex (RISC), which targets mRNAs for translational silencing and degradation as part of the RNA interference pathway. Despite the essential nature of this process for cellular function, there is little information on the role of RISC components in human development and organ function. We identify 13 heterozygous mutations in AGO2 in 21 patients affected by disturbances in neurological development. Each of the identified single amino acid mutations result in impaired shRNA-mediated silencing. We observe either impaired RISC formation or increased binding of AGO2 to mRNA targets as mutation specific functional consequences. The latter is supported by decreased phosphorylation of a C-terminal serine cluster involved in mRNA target release, increased formation of dendritic P-bodies in neurons and global transcriptome alterations in patient-derived primary fibroblasts. Our data emphasize the importance of gene expression regulation through the dynamic AGO2-RNA association for human neuronal development.

miR-122 and Ago interactions with the HCV genome alter the structure of the viral 5' terminus

Hepatitis C virus (HCV) is a positive-sense RNA virus that interacts with the liver-specific microRNA, miR-122. miR-122 binds to two sites in the 5' untranslated region (UTR) and this interaction promotes HCV RNA accumulation, although the precise role of miR-122 in the HCV life cycle remains unclear. Using biophysical analyses and Selective 2' Hydroxyl Acylation analyzed by Primer Extension (SHAPE) we investigated miR-122 interactions with the 5' UTR. Our data suggests that miR-122 binding results in alteration of nucleotides 1-117 to suppress an alternative secondary structure and promote functional internal ribosomal entry site (IRES) formation. Furthermore, we demonstrate that two hAgo2:miR-122 complexes are able to bind to the HCV 5' terminus simultaneously and SHAPE analyses revealed further alterations to the structure of the 5' UTR to accommodate these complexes. Finally, we present a computational model of the hAgo2:miR-122:HCV RNA complex at the 5' terminus of the viral genome as well as hAgo2:miR-122 interactions with the IRES-40S complex that suggest hAgo2 is likely to form additional interactions with SLII which may further stabilize the HCV IRES. Taken together, our results support a model whereby hAgo2:miR-122 complexes alter the structure of the viral 5' terminus and promote formation of the HCV IRES.

The Role of Tertiary Structure in MicroRNA Target Recognition

Translational repression and degradation of transcripts by microRNAs (miRNAs) is mediated by a ribonucleoprotein complex called the miRNA-induced silencing complex (miRISC, or RISC). Advances in experimental determination of RISC structures have enabled detailed analysis and modeling of known miRNA targets, yet a full appreciation of the structural factors influencing target recognition remains a challenge, primarily because target recognition involves a combination of RNA-RNA and RNA-protein interactions that can vary greatly among different miRNA-target pairs. In this chapter, we review progress toward understanding the role of tertiary structure in miRNA target recognition using computational approaches to assemble RISC complexes at known targets and physics-based methods for computing target interactions. Using this framework to examine RISC structures and dynamics, we describe how the conformational flexibility of Argonautes plays an important role in accommodating the diversity of miRNA-target duplexes formed at canonical and noncanonical target sites. We then discuss applications of tertiary structure-based approaches to emerging topics, including the structural effects of SNPs in miRNA targets and cooperative interactions involving Argonaute-Argonaute complexes. We conclude by assessing the prospects for genome-scale modeling of RISC structures and modeling of higher-order Argonaute complexes associated with miRNA biogenesis, mRNA regulation, and other functions.

Recent Molecular Genetic Explorations of Caenorhabditis elegans MicroRNAs

MicroRNAs are small, noncoding RNAs that regulate gene expression at the post-transcriptional level in essentially all aspects of Caenorhabditis elegans biology. More than 140 genes that encode microRNAs in C. elegans regulate development, behavior, metabolism, and responses to physiological and environmental changes. Genetic analysis of C. elegans microRNA genes continues to enhance our fundamental understanding of how microRNAs are integrated into broader gene regulatory networks to control diverse biological processes, including growth, cell division, cell fate determination, behavior, longevity, and stress responses. As many of these microRNA sequences and the related processing machinery are conserved over nearly a billion years of animal phylogeny, the assignment of their functions via worm genetics may inform the functions of their orthologs in other animals, including humans. In vivo investigations are especially important for microRNAs because in silico extrapolation of their functions using mRNA target prediction programs can easily assign microRNAs to incorrect genetic pathways. At this mezzanine level of microRNA bioinformatic sophistication, genetic analysis continues to be the gold standard for pathway assignments.

A non-canonical site reveals the cooperative mechanisms of microRNA-mediated silencing

Although strong evidence supports the importance of their cooperative interactions, microRNA (miRNA)-binding sites are still largely investigated as functionally independent regulatory units. Here, a survey of alternative 3΄UTR isoforms implicates a non-canonical seedless site in cooperative miRNA-mediated silencing. While required for target mRNA deadenylation and silencing, this site is not sufficient on its own to physically recruit miRISC. Instead, it relies on facilitating interactions with a nearby canonical seed-pairing site to recruit the Argonaute complexes. We further show that cooperation between miRNA target sites is necessary for silencing in vivo in the C. elegans embryo, and for the recruitment of the Ccr4-Not effector complex. Using a structural model of cooperating miRISCs, we identified allosteric determinants of cooperative miRNA-mediated silencing that are required for both embryonic and larval miRNA functions. Our results delineate multiple cooperative mechanisms in miRNA-mediated silencing and further support the consideration of target site cooperation as a fundamental characteristic of miRNA function.

MicroRNA-34 directly targets pair-rule genes and cytoskeleton component in the honey bee

MicroRNAs (miRNAs) are key regulators of developmental processes, such as cell fate determination and differentiation. Previous studies showed Dicer knockdown in honeybee embryos disrupt the processing of functional mature miRNAs and impairs embryo patterning. Here we investigated the expression profiles of miRNAs in honeybee embryogenesis and the role of the highly conserved miR-34-5p in the regulation of genes involved in insect segmentation. A total of 221 miRNAs were expressed in honey bee embryogenesis among which 97 mature miRNA sequences have not been observed before. Interestingly, we observed a switch in dominance between the 5-prime and 3-prime arm of some miRNAs in different embryonic stages; however, most miRNAs present one dominant arm across all stages of embryogenesis. Our genome-wide analysis of putative miRNA-target networks and functional pathways indicates miR-34-5p is one of the most conserved and connected miRNAs associated with the regulation of genes involved in embryonic patterning and development. In addition, we experimentally validated that miR-34-5p directly interacts to regulatory elements in the 3'-untranslated regions of pair-rule (even-skipped, hairy, fushi-tarazu transcription factor 1) and cytoskeleton (actin5C) genes. Our study suggests that miR-34-5p may regulate the expression of pair-rule and cytoskeleton genes during early development and control insect segmentation.

TargetLink, a new method for identifying the endogenous target set of a specific microRNA in intact living cells

MicroRNAs are small non-coding RNAs acting as posttranscriptional repressors of gene expression. Identifying mRNA targets of a given miRNA remains an outstanding challenge in the field. We have developed a new experimental approach, TargetLink, that applied locked nucleic acid (LNA) as the affinity probe to enrich target genes of a specific microRNA in intact cells. TargetLink also consists a rigorous and systematic data analysis pipeline to identify target genes by comparing LNA-enriched sequences between experimental and control samples. Using miR-21 as a test microRNA, we identified 12 target genes of miR-21 in a human colorectal cancer cell by this approach. The majority of the identified targets interacted with miR-21 via imperfect seed pairing. Target validation confirmed that miR-21 repressed the expression of the identified targets. The cellular abundance of the identified miR-21 target transcripts varied over a wide range, with some targets expressed at a rather low level, confirming that both abundant and rare transcripts are susceptible to regulation by microRNAs, and that TargetLink is an efficient approach for identifying the target set of a specific microRNA in intact cells. C20orf111, one of the novel targets identified by TargetLink, was found to reside in the nuclear speckle and to be reliably repressed by miR-21 through the interaction at its coding sequence.