SEARCH · Engineering Papers

Results for “protein function prediction”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 415 records · Page 23

All sources NASA NTRS DOE OSTI

DOE OSTI2012

Sequence-Specific Targeting of Dosage Compensation in Drosophila Favors an Active Chromatin Context

The Drosophila MSL complex mediates dosage compensation by increasing transcription of the single X chromosome in males approximately two-fold. This is accomplished through recognition of the X chromosome and subsequent acetylation of histone H4K16 on X-linked genes. Initial binding to the X is thought to occur at ‘‘entry sites’’ that contain a consensus sequence motif (‘‘MSL recognition element’’ or MRE). However, this motif is only ~2 fold enriched on X, and only a fraction of the motifs on X are initially targeted. Here we ask whether chromatin context could distinguish between utilized and non-utilized copies of the motif, by comparing their relative enrichment for histone modifications and chromosomal proteins mapped in the modENCODE project. Through a comparative analysis of the chromatin features in male S2 cells (which contain MSL complex) and female Kc cells (which lack the complex), we find that the presence of active chromatin modifications, together with an elevated local GC content in the surrounding sequences, has strong predictive value for functional MSL entry sites, independent of MSL binding. We tested these sites for function in Kc cells by RNAi knockdown of Sxl, resulting in induction of MSL complex. We show that ectopic MSL expression in Kc cells leads to H4K16 acetylation around these sites and a relative increase in X chromosome transcription. Collectively, our results support a model in which a pre-existing active chromatin environment, coincident with H3K36me3, contributes to MSL entry site selection. The consequences of MSL targeting of the male X chromosome include increase in nucleosome lability, enrichment for H4K16 acetylation and JIL-1 kinase, and depletion of linker histone H1 on active X-linked genes. Our analysis can serve as a model for identifying chromatin and local sequence features that may contribute to selection of functional protein binding sites in the genome.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2014

Comparative genomics and evolution of regulons of the LacI-family transcription factors

DNA-binding transcription factors (TFs) are essential components of transcriptional regulatory networks in bacteria. LacI-family TFs (LacI-TFs) are broadly distributed among certain lineages of bacteria. The majority of characterized LacI-TFs sense sugar effectors and regulate carbohydrate utilization genes. The comparative genomics approaches enable in silico identification of TF-binding sites and regulon reconstruction. To study the function and evolution of LacI-TFs, we performed genomics-based reconstruction and comparative analysis of their regulons. For over 1300 LacI-TFs from over 270 bacterial genomes, we predicted their cognate DNA-binding motifs and identified target genes. Using the genome context and metabolic subsystem analyses of reconstructed regulons, we tentatively assigned functional roles and predicted candidate effectors for 78 and 67% of the analyzed LacI-TFs, respectively. Nearly 90% of the studied LacI-TFs are local regulators of sugar utilization pathways, whereas the remaining 125 global regulators control large and diverse sets of metabolic genes. The global LacI-TFs include the previously known regulators CcpA in Firmicutes, FruR in Enterobacteria, and PurR in Gammaproteobacteria, as well as the three novel regulators—GluR, GapR, and PckR—that are predicted to control the central carbohydrate metabolism in three lineages of Alphaproteobacteria. Phylogenetic analysis of regulators combined with the reconstructed regulons provides a model of evolutionary diversification of the LacI protein family. The obtained genomic collection of in silico reconstructed LacI-TF regulons in bacteria is available in the RegPrecise database (http://regprecise.lbl.gov). It provides a framework for future structural and functional classification of the LacI protein family and identification of molecular determinants of the DNA and ligand specificity. The inferred regulons can be also used for functional gene annotation and reconstruction of sugar catabolic networks in diverse bacterial lineages.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS1995

Finding the global minimum: a fuzzy end elimination implementation

The 'fuzzy end elimination theorem' (FEE) is a mathematically proven theorem that identifies rotameric states in proteins which are incompatible with the global minimum energy conformation. While implementing the FEE we noticed two different aspects that directly affected the final results at convergence. First, the identification of a single dead-ending rotameric state can trigger a 'domino effect' that initiates the identification of additional rotameric states which become dead-ending. A recursive check for dead-ending rotameric states is therefore necessary every time a dead-ending rotameric state is identified. It is shown that, if the recursive check is omitted, it is possible to miss the identification of some dead-ending rotameric states causing a premature termination of the elimination process. Second, we examined the effects of removing dead-ending rotameric states from further considerations at different moments of time. Two different methods of rotameric state removal were examined for an order dependence. In one case, each rotamer found to be incompatible with the global minimum energy conformation was removed immediately following its identification. In the other, dead-ending rotamers were marked for deletion but retained during the search, so that they influenced the evaluation of other rotameric states. When the search was completed, all marked rotamers were removed simultaneously. In addition, to expand further the usefulness of the FEE, a novel method is presented that allows for further reduction in the remaining set of conformations at the FEE convergence. In this method, called a tree-based search, each dead-ending pair of rotamers which does not lead to the direct removal of either rotameric state is used to reduce significantly the number of remaining conformations. In the future this method can also be expanded to triplet and quadruplet sets of rotameric states. We tested our implementation of the FEE by exhaustively searching ten protein segments and found that the FEE identified the global minimum every time. For each segment, the global minimum was exhaustively searched in two different environments: (i) the segments were extracted from the protein and exhaustively searched in the absence of the surrounding residues; (ii) the segments were exhaustively searched in the presence of the remaining residues fixed at crystal structure conformations. We also evaluated the performance of the method for accurately predicting side chain conformations. We examined the influence of factors such as type and accuracy of backbone template used, and the restrictions imposed by the choice of potential function, parameterization and rotamer database. Conclusions are drawn on these results and future prospects are given.

NASA Program Exobiology↗

DOE OSTI2021

Efficient production of oxidized terpenoids via engineering fusion proteins of terpene synthase and cytochrome P450

The functionalization of terpenes using cytochrome P450 enzymes is a versatile route to the production of useful derivatives that can be further converted to value-added products. Many terpenes are hydrophobic and volatile making their availability as a substrate for P450 enzymes significantly limited during microbial production. In this study, we developed a strategy to improve the accessibility of terpene molecules for the P450 reaction by linking terpene synthase and P450 together. As a model system, fusion proteins of 1,8-cineole synthase (CS) and P450 cin were investigated and it showed an improved hydroxylation of the monoterpenoid 1,8-cineole up to 5.4-fold. Structural analysis of the CS-P450 cin fusion proteins by SEC-SAXS indicated a dimer formation with preferred orientations of the active sites of the two domains. We also applied the enzyme fusion strategy to the oxidation of a sesquiterpene epi-isozizaene and the fusion enzymes significantly improved albaflavenol production in engineered E. coli. From the analysis of positive and negative examples of the fusion strategy, we proposed key factors in structure-based prediction and evaluation of fusion enzymes. Developing fusion enzymes for terpene synthase and P450 presents an efficient strategy toward oxidation of hydrophobic terpene compounds. This strategy could be widely applicable to improve the biosynthetic titer of the functionalized products from hydrophobic terpene intermediates.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2017

Target highlights from the first post-PSI CASP experiment (CASP12, May–August 2016)

In this paper, the functional and biological significance of the selected CASP12 targets are described by the authors of the structures. Finally, the crystallographers discuss the most interesting structural features of the target proteins and assess whether these features were correctly reproduced in the predictions submitted to the CASP12 experiment.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2018

Diacylglycerol kinase γ predicts prognosis and functions as a tumor suppressor by negatively regulating glucose transporter 1 in hepatocellular carcinoma

Highlights: • DGKγ expression was downregulated in hepatocellular carcinoma (HCC). • Lower DGKγ expression in HCC tumor tissues was an independent prognostic factor for poor post-surgical overall survival. • Cell proliferation and migration of HCC cells were suppressed by DGKγ ectopic expression. • DGKγ overexpression led to down regulation of GLUT1 protein level and AMPK activity. Diacylglycerol kinases (DGK) are a family of enzymes catalyzing the transformation of diacylglycerol into phosphatidic acid, which have been recognized as key regulators in cell signaling pathways. The role of DGKγ in human malignancies has seldom been studied. In this study, we investigated the role of DGKγ in hepatocellular carcinoma (HCC). We found that DGKγ was down-regulated in HCC tumor tissues and cell lines as compared to that in non-tumor tissues. The prognostic value of DGKγ expression was evaluated by Cox regression and Kaplan-Meier analyses. Lower DGKγ expression in tumor tissues was an independent prognostic factor for poor post-surgical overall survival. By using HDACs inhibitors treatment and ChIP-PCR, we discovered that histone H3 and H4 deacetylation mainly contributed to the downregulation of DGKγ expression. Functional studies revealed that ectopic expression of DGKγ inhibited cell proliferation and cell migration in HCC cells. Mechanism studies showed that DGKγ overexpression led to down regulation of GLUT1 protein level and AMPK activity, which result in glucose uptake suppression as well as lactate and ATP production declination. The decrease of GLUT1 level could be partially rescued by treatments with either DGK inhibitor and lysosome inhibitor, indicating DGKγ may down-regulate GLUT1 through its kinase activity and lysosome degradation process. Together, this study demonstrated that DGKγ plays a tumor suppressor role in HCC by negatively regulating GLUT1. DGKγ could be a novel prognostic indicator and therapeutic target for HCC.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2017

Heterologous Gene Expression of N -Terminally Truncated Variants of LipPks1 Suggests a Functionally Critical Structural Motif in the N -terminus of Modular Polyketide Synthase

Streptomyces genomes have a high G + C content and typically use an ATG or GTG codon to initiate protein synthesis. Although gene-finding tools perform well in low GC genomes, it is known that the accuracy in predicting a translational start site (TSS) is much less for high GC genomes. LipPks1 is a Streptomyces-derived, well-characterized modular polyketide synthase (PKS). Here, using this enzyme as a model, we experimentally investigated the effects of alternative TSSs using a heterologous host, Streptomyces venezuelae. One of the TSSs employed boosted the protein level by 59-fold and the product yield by 23-fold compared to the originally annotated start codon. Interestingly, a structural model of the PKS indicated the presence of a structural motif in the N-terminus, which may explain the observed different protein levels together with a proline and arginine-rich sequence that may inhibit translational initiation. Lastly, this structure was also found in six other modular PKSs that utilize noncarboxylated starter substrates, which may guide the selection of optimal TSSs in conjunction with start-codon prediction software.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2017

The POTRA domains of Toc75 exhibit chaperone-like function to facilitate import into chloroplasts

Protein trafficking across membranes is an essential function in cells; however, the exact mechanism for how this occurs is not well understood. In the endosymbionts, mitochondria and chloroplasts, the vast majority of proteins are synthesized in the cytoplasm as preproteins and then imported into the organelles via specialized machineries. In chloroplasts, protein import is accomplished by the TOC (translocon on the outer chloroplast membrane) and TIC (translocon on the inner chloroplast membrane) machineries in the outer and inner envelope membranes, respectively. TOC mediates initial recognition of preproteins at the outer membrane and includes a core membrane channel, Toc75, and two receptor proteins, Toc33/34 and Toc159, each containing GTPase domains that control preprotein binding and translocation. Toc75 is predicted to have a β-barrel fold consisting of an N-terminal intermembrane space (IMS) domain and a C-terminal 16-stranded β-barrel domain. Here we report the crystal structure of the N-terminal IMS domain of Toc75 from Arabidopsis thaliana , revealing three tandem polypeptide transport-associated (POTRA) domains, with POTRA2 containing an additional elongated helix not observed previously in other POTRA domains. Functional studies show an interaction with the preprotein, preSSU, which is mediated through POTRA2-3. POTRA2-3 also was found to have chaperone-like activity in an insulin aggregation assay, which we propose facilitates preprotein import. Furthermore, our data suggest a model in which the POTRA domains serve as a binding site for the preprotein as it emerges from the Toc75 channel and provide a chaperone-like activity to prevent misfolding or aggregation as the preprotein traverses the intermembrane space.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2014

Chlamydia trachomatis protein CT009 is a structural and functional homolog to the key morphogenesis component RodZ and interacts with division septal plane localized MreB

Cell division in Chlamydiae is poorly understood as apparent homologs to most conserved bacterial cell division proteins are lacking and presence of elongation (rod shape) associated proteins indicate non-canonical mechanisms may be employed. The rod-shape determining protein MreB has been proposed as playing a unique role in chlamydial cell division. In other organisms, MreB is part of an elongation complex that requires RodZ for proper function. A recent study reported that the protein encoded by ORF CT009 interacts with MreB despite low sequence similarity to RodZ. The studies in this paper expand on those observations through protein structure, mutagenesis and cellular localization analyses. Structural analysis indicated that CT009 shares high level of structural similarity to RodZ, revealing the conserved orientation of two residues critical for MreB interaction. Substitutions eliminated MreB protein interaction and partial complementation provided by CT009 in RodZ deficient Escherichia coli. Cellular localization analysis of CT009 showed uniform membrane staining in Chlamydia. This was in contrast to the localization of MreB, which was restricted to predicted septal planes. Finally, MreB localization to septal planes provides direct experimental observation for the role of MreB in cell division and supports the hypothesis that it serves as a functional replacement for FtsZ in Chlamydia.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2018

Modulating Enzyme Activity by Altering Protein Dynamics with Solvent

Optimal enzyme activity depends on a number of factors, including structure and dynamics. The role of enzyme structure is well recognized; however, the linkage between protein dynamics and enzyme activity has given rise to a contentious debate. We have developed an approach that uses an aqueous mixture of organic solvent to control the functionally relevant enzyme dynamics (without changing the structure), which in turn modulates the enzyme activity. Using this approach, we predicted that the hydride transfer reaction catalyzed by the enzyme dihydrofolate reductase (DHFR) from Escherichia coli in aqueous mixtures of isopropanol (IPA) with water will decrease by ~3 fold at 20% (v/v) IPA concentration. Stopped-flow kinetic measurements find that the pH-independent k hydride rate decreases by 2.2 fold. X-ray crystallographic enzyme structures show no noticeable differences, while computational studies indicate that the transition state and electrostatic effects were identical for water and mixed solvent conditions; quasi-elastic neutron scattering studies show that the dynamical enzyme motions are suppressed. Our approach provides a unique avenue to modulating enzyme activity through changes in enzyme dynamics. Further it provides vital insights that show the altered motions of DHFR cause significant changes in the enzyme's ability to access its functionally relevant conformational substates, explaining the decreased k hydride rate. This approach has important implications for obtaining fundamental insights into the role of rate-limiting dynamics in catalysis and as well as for enzyme engineering.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2017

Human alpha beta hydrolase domain containing protein 11 and its yeast homolog are lipid hydrolases

Mammalian alpha/beta hydrolase domain (ABHD) family of proteins have emerged as key regulators of lipid metabolism and are found to be associated with human diseases. Human α/β-hydrolase domain containing protein 11 (ABHD11) has recently been predicted as a potential biomarker for human lung adenocarcinoma. In silico analyses of the ABHD11 protein sequence revealed the presence of a conserved lipase motif GXSXG. However, the role of ABHD11 in lipid metabolism is not known. To understand the biological function of ABHD11, we heterologously expressed the human ABHD11 in budding yeast, Saccharomyces cerevisiae. In vivo [{sup 14}C]acetate labeling of cellular lipids in yeast cells overexpressing ABHD11 showed a decrease in triacylglycerol content. Overexpression of ABHD11 also alters the molecular species of triacylglycerol in yeast. Similar activity was observed in its yeast homolog, Ygr031w. The role of the conserved lipase motif in the hydrolase activity was proven by the mutation of all conserved amino acid residues of GXSXG motif. Collectively, our results demonstrate that human ABHD11 and its yeast homolog YGR031W have a pivotal role in the lipid metabolism. - Highlights: • Overexpression of ABHD11 protein and its yeast homolog Ygr031w cause a reduction in triacylglycerol levels in yeast. • The reduction in triacylglycerol is due to the presence of lipase motif GXSXG. • Overexpression of ABHD11 and Ygr031w alters the molecular species of triacylglycerol.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2019

Evolution and comparative genomics of the most common Trichoderma species

Background: The growing importance of the ubiquitous fungal genus Trichoderma (Hypocreales, Ascomycota) requires understanding of its biology and evolution. Many Trichoderma species are used as biofertilizers and biofungicides and T. reesei is the model organism for industrial production of cellulolytic enzymes. In addition, some highly opportunistic species devastate mushroom farms and can become pathogens of humans. A comparative analysis of the first three whole genomes revealed mycoparasitism as the innate feature of Trichoderma. However, the evolution of these traits is not yet understood. Results: We selected 12 most commonly occurring Trichoderma species and studied the evolution of their genome sequences. Trichoderma evolved in the time of the Cretaceous-Palaeogene extinction event 66 (±15) mya, but the formation of extant sections (Longibrachiatum, Trichoderma) or clades (Harzianum/Virens) happened in Oligocene. The evolution of the Harzianum clade and section Trichoderma was accompanied by significant gene gain, but the ancestor of section Longibrachiatum experienced rapid gene loss. The highest number of genes gained encoded ankyrins, HET domain proteins and transcription factors. We also identified the Trichoderma core genome, completely curated its annotation, investigated several gene families in detail and compared the results to those of other fungi. Eighty percent of those genes for which a function could be predicted were also found in other fungi, but only 67% of those without a predictable function. Conclusions: Our study presents a time scaled pattern of genome evolution in 12 Trichoderma species from three phylogenetically distant clades/sections and a comprehensive analysis of their genes. The data offer insights in the evolution of a mycoparasite towards a generalist.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2013

Hydrogen Regulation and Global Responses to Electron, Carbon and Nitrogen Sources of Methanococcus Maripaludis

Methanogens catalyze the critical, methane-producing step (called methanogenesis) in the anaerobic decomposition of organic matter. This project has generated the first predictive model of global gene regulation of methanogenesis in a hydrogenotrophic methanogen, Methanococcus maripaludis. We generated a comprehensive list of genes (protein-coding and non-coding) for M. maripaludis through integrated analysis of the transcriptome structure and a newly constructed Peptide Atlas. The environment and gene-regulatory influence network (EGRIN) model of the strain was constructed from a compendium of transcriptome data that was collected over 58 different steady-state and time course experiments that were performed in chemostats or batch cultures, under a spectrum of environmental perturbations that modulated methanogenesis. Analyses of the EGRIN model have revealed novel components of methanogenesis that included at least three additional protein-coding genes of previously unknown function as well as one non-coding RNA. We discovered that at least five regulatory mechanisms act in a combinatorial scheme to inter-coordinate key steps of methanogenesis with different processes such as motility, ATP biosynthesis, and carbon assimilation. Through a combination of genetic and environmental perturbation experiments we have validated the EGRIN-predicted role of two novel transcription factors in the regulation of phosphate-dependent repression of formate dehydrogenase a key enzyme in the methanogenesis pathway. The EGRIN model demonstrates regulatory affiliations within methanogenesis as well as between methanogenesis and other cellular functions. In addition, we have published an analysis of transcriptome architecture in M. maripaludis and an analysis of the effects of H 2 and formate on growth yield and regulation of methanogenesis in M. maripaludis.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2015

Structural insights into the assembly of the histone deacetylase-associated Sin3L/Rpd3L corepressor complex

Acetylation is correlated with chromatin decondensation and transcriptional activation, but its regulation by histone deacetylase (HDAC)-bearing corepressor complexes is poorly understood. Here, we describe the mechanism of assembly of the mammalian Sin3L/Rpd3L complex facilitated by Sds3, a conserved subunit deemed critical for proper assembly. Sds3 engages a globular, helical region of the HDAC interaction domain (HID) of the scaffolding protein Sin3A through a bipartite motif comprising a helix and an adjacent extended segment. Sds3 dimerizes through not only one of the predicted coiled-coil motifs but also, the segment preceding it, forming an ~150-Å-long antiparallel dimer. Contrary to previous findings in yeast, Sin3A rather than Sds3 functions in recruiting HDAC1 into the complex by engaging the latter through a highly conserved segment adjacent to the helical HID subdomain. In the resulting model for the ternary complex, the two copies of the HDACs are situated distally and dynamically because of a natively unstructured linker connecting the dimerization domain and the Sin3A interaction domain of Sds3; these features contrast with the static organization described previously for the NuRD (nucleosome remodeling and deacetylase) complex. Finally, the Sds3 linker features several conserved basic residues that could potentially maintain the complex on chromatin by nonspecific interactions with DNA after initial recruitment by sequence-specific DNA-binding repressors.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Systematic identification of transcriptional activation domains from non-transcription factor proteins in plants and yeast

Transcription factors can promote gene expression through activation domains. Whole-genome screens have systematically mapped activation domains in transcription factors but not in non-transcription factor proteins (e.g., chromatin regulators and coactivators). To fill this knowledge gap, we employed the activation domain predictor PADDLE to analyze the proteomes of Arabidopsis thaliana and Saccharomyces cerevisiae. We screened 18,000 predicted activation domains from >800 non-transcription factor genes in both species, confirming that 89% of candidate proteins contain active fragments. Our work enables the annotation of hundreds of nuclear proteins as putative coactivators, many of which have never been ascribed any function in plants. Analysis of peptide sequence compositions reveals how the distribution of key amino acids dictates activity. Finally, we validated short, "universal" activation domains with comparable performance to state-of-the-art activation domains used for genome engineering. Our approach enables the genome-wide discovery and annotation of activation domains that can function across diverse eukaryotes.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2013

Assembly-Driven Community Genomics of a Hypersaline Microbial Ecosystem

Microbial populations inhabiting a natural hypersaline lake ecosystem in Lake Tyrrell, Victoria, Australia, have been characterized using deep metagenomic sampling, iterative de novo assembly, and multidimensional phylogenetic binning. Composite genomes representing habitat-specific microbial populations were reconstructed for eleven different archaea and one bacterium, comprising between 0.6 and 14.1% of the planktonic community. Eight of the eleven archaeal genomes were from microbial species without previously cultured representatives. These new genomes provide habitat-specific reference sequences enabling detailed, lineage-specific compartmentalization of predicted functional capabilities and cellular properties associated with both dominant and less abundant community members, including organisms previously known only by their 16S rRNA sequences. Together, these data provide a comprehensive, culture-independent genomic blueprint for ecosystem-wide analysis of protein functions, population structure, and lifestyles of co-existing, co-evolving microbial groups within the same natural habitat. The ‘‘assembly-driven’’ community genomic approach demonstrated in this study advances our ability to push beyond single gene investigations, and promotes genome-scale reconstructions as a tangible goal in the quest to define the metabolic, ecological, and evolutionary dynamics that underpin environmental microbial diversity.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2014

New In Situ Imaging and Measurement Technologies for Biological Systems (Brochure)

The U.S. Department of Energy’s (DOE) Office of Biological and Environmental Research (BER) supports fundamental research to advance a predictive understanding of complex biological and environmental systems relevant to DOE’s missions in energy and the environment. Starting with the genetic potential encoded by organisms’ genomes, BER scientists seek to define the principles guiding translation of the genetic code into functional proteins and the metabolic and regulatory networks underlying plant and microbial systems. Concurrent with this research is a need for enabling technologies to place understanding of gene expression, regulation, and function into the spatiotemporal context of whole-cell environments. Imaging and measurement technologies that can resolve multiple key metabolic processes over time within or among cells will act as a crucial bridge toward linking molecular-scale information to whole-cell, systems-level understanding. This brochure outlines the pilot projects at the DOE national laboratories and the investigators leading these efforts. This initial description of BER’s new effort serves as an informational resource for the evolving imaging and measurement technology program.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2019

Identifying RNA contacts from SHAPE-MaP by partial correlation analysis

In a recent paper Siegfried et al. published a new sequence-based structural RNA assay that utilizes mutational profiling to detect base pairing (MaP). Output from MaP provides information about both pairing (via reactivities) and contact (via correlations). Reactivities can be coupled to partition function folding models for structural inference, while correlations can reveal pairs of sites that may be in structural proximity. The possibility for inference of 3D contacts via MaP suggests a novel approach to structural prediction for RNA analogous to covariance structural prediction for proteins. We explore this approach and show that partial correlation analysis outperforms na\"ive correlation analysis. Our results should be applicable to a wide range of high-throughput sequencing based RNA structural assays that are under development.

Tambe, Akshay↗