Selection of allosteric dnazymes that can sense phenylalanine by expression-SELEX.

Aptamers are ligand-binding RNA or DNA molecules and have been widely examined as biosensors, diagnostic tools, and therapeutic agents. The application of aptamers as biosensors commonly requires an expression platform to produce a signal to report the aptamer-ligand binding event. Traditionally, aptamer selection and expression platform integration are two independent steps and the aptamer selection requires the immobilization of either the aptamer or the ligand. These drawbacks can be easily overcome through the selection of allosteric DNAzymes (aptazymes). Herein, we used the technique of Expression-SELEX developed in our laboratory to select for aptazymes that can be specifically activated by low concentrations of l-phenylalanine. We chose a previous DNA-cleaving DNAzyme known as II-R1 as the expression platform for its low cleavage rate and used stringent selection conditions to drive the selection of high-performance aptazyme candidates. Three aptazymes were chosen for detailed characterization and these DNAzymes were found to exhibit a dissociation constant for l-phenylalanine as low as 4.8 µM, a catalytic rate constant improvement as high as 20 000-fold in the presence of l-phenylalanine, and the ability to discriminate against closely related l-phenylalanine analogs including d-phenylalanine. This work has established the Expression-SELEX as an effective SELEX method to enrich high-quality ligand-responsive aptazymes.

Integration of an Expression Platform in the SELEX Cycle to Select DNA Aptamer Binding to a Disease Biomarker.

Aptamers can be developed for biosensors, diagnostic tools, and therapeutic reagents. These applications usually require a fusion of aptamers and expression platforms. However, the fusion process is usually time-consuming and laborious. In this study, we integrated the deoxyribozyme (I-R3) as an expression platform in the SELEX cycle (called Expression-SELEX) to select aptazymes that can sense diverse molecules. We used the Maple syrup urine disease (MSUD) biomarker L-allo-isoleucine to test the selection model. After five rounds of screening, the cleavage products were sufficiently enriched to be visualized on polyacrylamide gel electrophoresis (PAGE) gel. Through high-throughput sequencing analysis, several candidates were identified. One such candidate, IR3-I-DNA, binds L-allo-isoleucine with a dissociation constant (K D) of 0.57 mM. When the ligand was present, the cleavage fraction of IR3-I-DNA increased from 0.3 to 0.5, and its K obs value improved from 1.38 min-1 to 1.97 min-1. Our selection approach can also be applied to produce aptazymes that can bind to variable ligands and be used more directly as biosensors.

The application of riboswitch sequencing for human gut bacterial classification and identification.

Bacteria, especially gut bacteria play important roles in human health and diseases. The classification of many bacterial genera by the 16S ribosomal RNA (rRNA) has failed due to its low inter-species resolution. Given the wide distribution of riboswitches in bacteria, they may help 16S rRNA differentiate closely related species. We found that among 28 groups of species that could not be distinguished by 16S rRNA, eight of them could be separated by the TPP riboswitch and other riboswitches. Moreover, the species in the 16S rRNA database and these riboswitch databases overlap, therefore, using riboswitch databases can help 16S rRNA better identify species. In addition, we used Klenow DNA polymerase and a pair of short primers to facilitate the library construction of TPP riboswitches for sequencing. The sequencing results showed that the TPP riboswitch could detect the major phyla similar to those detected by 16S rRNA. Therefore, the TPP riboswitch and other riboswitch classes could potentially be applied to gut bacteria classification.

Identification of 11 candidate structured noncoding RNA motifs in humans by comparative genomics.

BACKGROUND: Only 1.5% of the human genome encodes proteins, while large part of the remaining encodes noncoding RNAs (ncRNA). Many ncRNAs form structures and perform many important functions. Accurately identifying structured ncRNAs in the human genome and discovering their biological functions remain a major challenge. RESULTS: Here, we have established a pipeline (CM-line) with the following features for analyzing the large genomes of humans and other animals. First, we selected species with larger genetic distances to facilitate the discovery of covariations and compatible mutations. Second, we used CMfinder, which can generate useful alignments even with low sequence conservation. Third, we removed repetitive sequences and known structured ncRNAs to reduce the workload of CMfinder. Fourth, we used Infernal to find more representatives and refine the structure. We reported 11 classes of structured ncRNA candidates with significant covariations in humans. Functional analysis showed that these ncRNAs may have variable functions. Some may regulate circadian clock genes through poly (A) signals (PAS); some may regulate the elongation factor (EEF1A) and the T-cell receptor signaling pathway by cooperating with RNA binding proteins. CONCLUSIONS: By searching for important features of RNA structure from large genomes, the CM-line has revealed the existence of a variety of novel structured ncRNAs. Functional analysis suggests that some newly discovered ncRNA motifs may have biological functions. The pipeline we have established for the discovery of structured ncRNAs and the identification of their functions can also be applied to analyze other large genomes.

DeepPSP: A Global-Local Information-Based Deep Neural Network for the Prediction of Protein Phosphorylation Sites.

Identification of phosphorylation sites is an important step in the function study and drug design of proteins. In recent years, there have been increasing applications of the computational method in the identification of phosphorylation sites because of its low cost and high speed. Most of the currently available methods focus on using local information around potential phosphorylation sites for prediction and do not take the global information of the protein sequence into consideration. Here, we demonstrated that the global information of protein sequences may be also critical for phosphorylation site prediction. In this paper, a new deep neural network model, called DeepPSP, was proposed for the prediction of protein phosphorylation sites. In the DeepPSP model, two parallel modules were introduced to extract both local and global features from protein sequences. Two squeeze-and-excitation blocks and one bidirectional long short-term memory block were introduced into each module to capture effective representations of the sequences. Comparative studies were carried out to evaluate the performance of DeepPSP, and four other prediction methods using public data sets The F1-score, area under receiver operating characteristic curves (AUROC), and area under precision-recall curves (AUPRC) of DeepPSP were found to be 0.4819, 0.82, and 0.50, respectively, for S/T general site prediction and 0.4206, 0.73, and 0.39, respectively, for Y general site prediction. Compared with the MusiteDeep method, the F1-score, AUROC, and AUPRC of DeepPSP were found to increase by 8.6, 2.5, and 8.7%, respectively, for S/T general site prediction and by 20.6, 5.8, and 18.2%, respectively, for Y general site prediction. Among the tested methods, the developed DeepPSP method was also found to produce best results for different kinase-specific site predictions including CDK, mitogen-activated protein kinase, CAMK, AGC, and CMGC. Taken together, the developed DeepPSP method may offer a more accurate phosphorylation site prediction by including global information. It may serve as an alternative model with better performance and interpretability for protein phosphorylation site prediction.

Identification of 15 candidate structured noncoding RNA motifs in fungi by comparative genomics.

BACKGROUND: With the development of rapid and inexpensive DNA sequencing, the genome sequences of more than 100 fungal species have been made available. This dataset provides an excellent resource for comparative genomics analyses, which can be used to discover genetic elements, including noncoding RNAs (ncRNAs). Bioinformatics tools similar to those used to uncover novel ncRNAs in bacteria, likewise, should be useful for searching fungal genomic sequences, and the relative ease of genetic experiments with some model fungal species could facilitate experimental validation studies. RESULTS: We have adapted a bioinformatics pipeline for discovering bacterial ncRNAs to systematically analyze many fungal genomes. This comparative genomics pipeline integrates information on conserved RNA sequence and structural features with alternative splicing information to reveal fungal RNA motifs that are candidate regulatory domains, or that might have other possible functions. A total of 15 prominent classes of structured ncRNA candidates were identified, including variant HDV self-cleaving ribozyme representatives, atypical snoRNA candidates, and possible structured antisense RNA motifs. Candidate regulatory motifs were also found associated with genes for ribosomal proteins, S-adenosylmethionine decarboxylase (SDC), amidase, and HexA protein involved in Woronin body formation. We experimentally confirm that the variant HDV ribozymes undergo rapid self-cleavage, and we demonstrate that the SDC RNA motif reduces the expression of SAM decarboxylase by translational repression. Furthermore, we provide evidence that several other motifs discovered in this study are likely to be functional ncRNA elements. CONCLUSIONS: Systematic screening of fungal genomes using a computational discovery pipeline has revealed the existence of a variety of novel structured ncRNAs. Genome contexts and similarities to known ncRNA motifs provide strong evidence for the biological and biochemical functions of some newly found ncRNA motifs. Although initial examinations of several motifs provide evidence for their likely functions, other motifs will require more in-depth analysis to reveal their functions.

The yjdF riboswitch candidate regulates gene expression by binding diverse azaaromatic compounds.

The yjdF motif RNA is an orphan riboswitch candidate that almost exclusively associates with the yjdF protein-coding gene in many bacteria. The function of the YjdF protein is unknown, which has made speculation regarding the natural ligand for this putative riboswitch unusually challenging. By using a structure-probing assay for ligand binding, we found that a surprisingly broad diversity of nitrogen-containing aromatic heterocycles, or "azaaromatics," trigger near-identical changes in the structures adopted by representative yjdF motif RNAs. Regions of the RNA that undergo ligand-induced structural modulation reside primarily in portions of the putative aptamer region that are highly conserved in nucleotide sequence, as is typical for riboswitches. Some azaaromatic molecules are bound by the RNA with nanomolar dissociation constants, and a subset of these ligands activate riboswitch-mediated gene expression in cells. Furthermore, genetic elements most commonly adjacent to the yjdF motif RNA or to the yjdF protein-coding region are homologous to protein regulators implicated in mitigating the toxic effects of diverse phenolic acids or polycyclic compounds. Although the precise type of natural ligand sensed by yjdF motif RNAs remains unknown, our findings suggest that this riboswitch class might serve as part of a genetic response system to toxic or signaling compounds with chemical structures similar to azaaromatics.

Biochemical analysis of pistol self-cleaving ribozymes.

Pistol RNAs are members of a distinct class of self-cleaving ribozymes that was recently discovered by using a bioinformatics search strategy. Several hundred pistol ribozymes share a consensus sequence including 10 highly conserved nucleotides and many other modestly conserved nucleotides associated with specific secondary structure features, including three base-paired stems and a pseudoknot. A representative pistol ribozyme from the bacterium Lysinibacillus sphaericus was found to promote RNA strand scission with a rate constant of ∼10 min(-1) under physiological Mg(2+) and pH conditions. The reaction proceeds via the nucleophilic attack of a 2'-oxygen atom on the adjacent phosphorus center, and thus adheres to the same general catalytic mechanism of internal phosphoester transfer as found with all other classes of natural self-cleaving ribozymes discovered to date. Analyses of the kinetic characteristics and the metal ion requirements of the cleavage reaction reveal that members of this ribozyme class likely use several catalytic strategies to promote the rapid cleavage of RNA.

Biochemical analysis of hatchet self-cleaving ribozymes.

Hatchet RNAs are members of a novel self-cleaving ribozyme class that was recently discovered by using a bioinformatics search strategy. The consensus sequence and secondary structure of this class includes 13 highly conserved and numerous other modestly conserved nucleotides interspersed among bulges linking four base-paired substructures. A representative hatchet ribozyme from a metagenomic source requires divalent ions such as Mg(2+) to promote RNA strand scission with a maximum rate constant of ∼4 min(-1). As with all other small self-cleaving ribozymes discovered to date, hatchet ribozymes employ a general mechanism for catalysis involving the nucleophilic attack of a ribose 2'-oxygen atom on an adjacent phosphorus center. Kinetic characteristics of the reaction demonstrate that members of this ribozyme class have an essential requirement for divalent metal ions and that they might have a complex active site that employs multiple catalytic strategies to accelerate RNA cleavage by internal phosphoester transfer.

New classes of self-cleaving ribozymes revealed by comparative genomics analysis.

Enzymes made of RNA catalyze reactions that are essential for protein synthesis and RNA processing. However, such natural ribozymes are exceedingly rare, as evidenced by the fact that the discovery rate for new classes has dropped to one per decade from about one per year during the 1980s. Indeed, only 11 distinct ribozyme classes have been experimentally validated to date. Recently, we recognized that self-cleaving ribozymes frequently associate with certain types of genes from bacteria. Herein we exploited this association to identify divergent architectures for two previously known ribozyme classes and to discover additional noncoding RNA motifs that are self-cleaving RNA candidates. We identified three new self-cleaving classes, which we named twister sister, pistol and hatchet, from this collection, suggesting that even more ribozymes remain hidden in modern cells.

Eukaryotic resistance to fluoride toxicity mediated by a widespread family of fluoride export proteins.

Fluorine is an abundant element and is toxic to organisms from bacteria to humans, but the mechanisms by which eukaryotes resist fluoride toxicity are unknown. The Escherichia coli gene crcB was recently shown to be regulated by a fluoride-responsive riboswitch, implicating it in fluoride response. There are >8,000 crcB homologs across all domains of life, indicating that it has an important role in biology. Here we demonstrate that eukaryotic homologs [renamed FEX (fluoride exporter)] function in fluoride export. FEX KOs in three eukaryotic model organisms, Neurospora crassa, Saccharomyces cerevisiae, and Candida albicans, are highly sensitized to fluoride (>200-fold) but not to other halides. Some of these KO strains are unable to grow in fluoride concentrations found in tap water. Using the radioactive isotope of fluoride, (18)F, we developed an assay to measure the intracellular fluoride concentration and show that the FEX deletion strains accumulate fluoride in excess of the external concentration, providing direct evidence of FEX function in fluoride efflux. In addition, they are more sensitive to lower pH in the presence of fluoride. These results demonstrate that eukaryotic FEX genes encode a previously unrecognized class of fluoride exporter necessary for survival in standard environmental conditions.

Eukaryotic TPP riboswitch regulation of alternative splicing involving long-distance base pairing.

Thiamin pyrophosphate (TPP) riboswitches are found in organisms from all three domains of life. Examples in bacteria commonly repress gene expression by terminating transcription or by blocking ribosome binding, whereas most eukaryotic TPP riboswitches are predicted to regulate gene expression by modulating RNA splicing. Given the widespread distribution of eukaryotic TPP riboswitches and the diversity of their locations in precursor messenger RNAs (pre-mRNAs), we sought to examine the mechanism of alternative splicing regulation by a fungal TPP riboswitch from Neurospora crassa, which is mostly located in a large intron separating protein-coding exons. Our data reveal that this riboswitch uses a long-distance (∼530-nt separation) base-pairing interaction to regulate alternative splicing. Specifically, a portion of the TPP-binding aptamer can form a base-paired structure with a conserved sequence element (α) located near a 5' splice site, which greatly increases use of this 5' splice site and promotes gene expression. Comparative sequence analyses indicate that many fungal species carry a TPP riboswitch with similar intron architecture, and therefore the homologous genes in these fungi are likely to use the same mechanism. Our findings expand the scope of genetic control mechanisms relying on long-range RNA interactions to include riboswitches.

Fluoride enhances the activity of fungicides that destabilize cell membranes.

Fluoride has long been known to inhibit bacterial and fungal cell growth most likely by blocking the functions of key metabolic enzymes. In this study, we demonstrate that antifungal compounds that disrupt cell membrane integrity exhibit improved ability to inhibit cell growth when used with millimolar concentrations of fluoride. Specifically, antifungal compounds of the polyene class and an antifungal peptide exhibit synergy with fluoride to inhibit the growth of various fungal species, including Candida albicans. Our results demonstrate that certain compounds can be found that increase the cellular uptake of fluoride, and provide new opportunities for creating antimicrobial compounds whose functions are enhanced when combined with otherwise sub-inhibitory concentrations of small ions.

A new mutation affecting FRQ-less rhythms in the circadian system of Neurospora crassa.

We are using the fungus Neurospora crassa as a model organism to study the circadian system of eukaryotes. Although the FRQ/WCC feedback loop is said to be central to the circadian system in Neurospora, rhythms can still be seen under many conditions in FRQ-less (frq knockout) strains. To try to identify components of the FRQ-less oscillator (FLO), we carried out a mutagenesis screen in a FRQ-less strain and selected colonies with altered conidiation (spore-formation) rhythms. A mutation we named UV90 affects rhythmicity in both FRQ-less and FRQ-sufficient strains. The UV90 mutation affects FRQ-less rhythms in two conditions: the free-running long-period rhythm in choline-depleted chol-1 strains becomes arrhythmic, and the heat-entrained rhythm in the frq(10) knockout is severely altered. In a FRQ-sufficient background, the UV90 mutation causes damping of the free-running conidiation rhythm, reduction of the amplitude of the FRQ protein rhythm, and increased phase-resetting responses to both light and heat pulses, consistent with a decreased amplitude of the circadian oscillator. The UV90 mutation also has small but significant effects on the period of the conidiation rhythm and on growth rate. The wild-type UV90 gene product appears to be required for a functional FLO and for sustained, high-amplitude rhythms in FRQ-sufficient conditions. The UV90 gene product may therefore be a good candidate for a component of the FRQ-less oscillator. These results support a model of the Neurospora circadian system in which the FRQ/WCC feedback loop mutually interacts with a single FLO in an integrated circadian system.

Effects of prd circadian clock mutations on FRQ-less rhythms in Neurospora.

Rhythmic conidiation (spore formation) in Neurospora crassa provides a model system for investigating the molecular mechanisms of circadian rhythmicity. A feedback loop involving the frq, wc-1, and wc-2 gene products (FRQ/ WCC) is an important component of the mechanism; however, rhythmic conidiation can still be observed when these gene products are absent. The nature of the oscillator(s) that drives this FRQ-less rhythmicity (FLO) is an important question in Neurospora circadian biology. We have looked for interactions between FRQ/WCC and FLO by assaying the effects on FRQ-less rhythms of mutations known to affect the period in the presence of FRQ. We assayed 4 prd mutations (prd-1, prd-2, prd-3, and prd-4) under 2 conditions in frq(null) strains: long-period free-running rhythms in chol-1 strains grown without choline, and heat-entrainable rhythms in choline-sufficient conditions. We found effects of all 4 mutations on both types of FRQ-less rhythms. The greatest effects were seen with prd-1 and prd-2, which abolished free-running rhythms in the chol-1; frq(10) backgrounds and significantly affected entrained peak timing under heat-entrainment conditions in frq( 10) backgrounds. The prd-3 and prd-4 mutations had more subtle effects on period and stability of free-running rhythms in the chol-1; frq(10) backgrounds and had little effect on peak timing under heat-entrainment conditions in frq(10) backgrounds. These results, along with previously published evidence for effects of prd mutations on other FRQ-less rhythms, suggest that either there are common components shared between the FRQ/WCC oscillator and several FRQ-less oscillators or that there is a single oscillator driving all conidiation rhythms. We favor a model of the Neurospora circadian system in which a single FRQ-less oscillator drives conidiation and interacts with the FRQ/WCC feedback loop; the output or amplitude of the FRQ-less oscillator can be affected by many gene products and metabolic conditions that reveal FRQ-less rhythmicity. We propose that prd-1 and prd-2 are good candidates for components of the FRQ-less oscillator and that prd-3 and prd-4 act on the system mainly through effects on FRQ/WCC.

Rhythmic conidiation in constant light in vivid mutants of Neurospora crassa.

In Neurospora crassa, a circadian rhythm of conidiation (asexual spore formation) can be seen on the surface of agar media. This rhythm has a period of 22 hr in constant darkness (D/D). Under constant illumination (L/L), no rhythm is visible and cultures show constant conidiation. However, here we report that strains with a mutation in the vivid (vvd) gene, previously shown to code for the photoreceptor involved in photo-adaptation, exhibit conidiation rhythms in L/L as well as in D/D. The period of the rhythm of vvd strains ranges between 6 and 21 hr in L/L, depending upon the intensity of the light, the carbon source, and the presence of other mutations. Temperature compensation of the period also depends on light intensity. Dark pulses given in L/L shift the phase of the rhythm. Shifts from L/L to D/D show unexpected after effects; i.e., the short period of a vvd strain in L/L gradually lengthens over 2-3 days in D/D. The rhythm in L/L requires the white collar (wc-1) gene, but not the frequency (frq) gene. FRQ protein shows no rhythm in L/L in a vvd strain. The conidiation rhythm in L/L in vvd is therefore driven by a FRQ-less oscillator (FLO).