Tertiary Structural Elements Research Articles

Identification of novel candidate predisposing genes in familial nonmedullary thyroid carcinoma implicating DNA damage repair pathways.

The genetic basis of nonsyndromic familial nonmedullary thyroid carcinoma (FNMTC) is still poorly understood, as the susceptibility genes identified so far only account for a small percentage of the genetic burden. Recently, germline mutations in DNA repair-related genes have been reported in cases with thyroid cancer. In order to clarify the genetic basis of FNMTC, 94 genes involved in hereditary cancer predisposition, including DNA repair genes, were analyzed in 48 probands from FNMTC families, through targeted next-generation sequencing (NGS). Genetic variants were selected upon bioinformatics analysis and in silico studies. Structural modeling and network analysis were also performed. In silico results of NGS data unveiled likely pathogenic germline variants in 15 families with FNMTC, in genes encoding proteins involved in DNA repair (ATM, CHEK2, ERCC2, BRCA2, ERCC4, FANCA, FANCD2, FANCF, and PALB2) and in the DICER1, FLCN, PTCH1, BUB1B, and RHBDF2 genes. Structural modeling predicted that most missense variants resulted in the disruption of networks of interactions between residues, with implications for local secondary and tertiary structure elements. Functional annotation and network analyses showed that the involved DNA repair proteins functionally interact with each other, within the same DNA repair pathway and across different pathways. MAPK activation was a common event in tumor progression. This study supports that rare germline variants in DNA repair genes may be accountable for FNMTC susceptibility, with potential future utility in patients' clinical management, and reinforces the relevance of DICER1 in disease etiology.

NMR and Single-Molecule FRET Insights into Fast Protein Motions and Their Relation to Function.

Proteins often undergo large-scale conformational transitions, in which secondary and tertiary structure elements (loops, helices, and domains) change their structures or their positions with respect to each other. Simple considerations suggest that such dynamics should be relatively fast, but the functional cycles of many proteins are often relatively slow. Sophisticated experimental methods are starting to tackle this dichotomy and shed light on the contribution of large-scale conformational dynamics to protein function. In this review, we focus on the contribution of single-molecule Förster resonance energy transfer and nuclear magnetic resonance (NMR) spectroscopies to the study of conformational dynamics. We briefly describe the state of the art in each of these techniques and then point out their similarities and differences, as well as the relative strengths and weaknesses of each. Several case studies, in which the connection between fast conformational dynamics and slower function has been demonstrated, are then introduced and discussed. These examples include both enzymes and large protein machines, some of which have been studied by both NMR and fluorescence spectroscopies.

Watching ion-driven kinetics of ribozyme folding and misfolding caused by energetic and topological frustration one molecule at a time.

Folding of ribozymes into well-defined tertiary structures usually requires divalent cations. How Mg2+ ions direct the folding kinetics has been a long-standing unsolved problem because experiments cannot detect the positions and dynamics of ions. To address this problem, we used molecular simulations to dissect the folding kinetics of the Azoarcus ribozyme by monitoring the path each molecule takes to reach the folded state. We quantitatively establish that Mg2+ binding to specific sites, coupled with counter-ion release of monovalent cations, stimulate the formation of secondary and tertiary structures, leading to diverse pathways that include direct rapid folding and trapping in misfolded structures. In some molecules, key tertiary structural elements form when Mg2+ ions bind to specific RNA sites at the earliest stages of the folding, leading to specific collapse and rapid folding. In others, the formation of non-native base pairs, whose rearrangement is needed to reach the folded state, is the rate-limiting step. Escape from energetic traps, driven by thermal fluctuations, occurs readily. In contrast, the transition to the native state from long-lived topologically trapped native-like metastable states is extremely slow. Specific collapse and formation of energetically or topologically frustrated states occur early in the assembly process.

Expression of membrane beta-barrel protein in E. coli at low temperatures: Structure of Yersinia pseudotuberculosis OmpF porin inclusion bodies

The recombinant OmpF porin of Yersinia pseudotuberculosis as a model of transmembrane protein of the β-barrel structural family was used to study low growth temperature effect on the structure of the produced inclusion bodies (IBs). This porin showed a very low expression level in E. coli at a growth temperature below optimal 37 °C. The introduction of a N-terminal hexahistidine tag into the mature porin molecule significantly increased the biosynthesis of the protein at low cultivation temperatures. The recombinant His-tagged porin (rOmpF-His) was expressed in E. coli at 30 and 18 °C as inclusion bodies (IB-30 and IB-18). The properties and structural organization of IBs, as well as the structure of rOmpF-His solubilized from the IBs with urea and SDS, were studied using turbidimetry, electron microscopy, dynamic light scattering, optical spectroscopy, and amyloid-specific dyes. IB-18, in comparison with IB-30, has a higher solubility in denaturants, suggesting a difference between IBs in the conformation of the associated polypeptide chains. The spectroscopic analysis revealed that rOmpF-His IBs have a high content of secondary structure with a tertiary-structure elements, including a native-like conformation, the proportion of which in IB-18 is higher than in IB-30. Solubilization of the porin from IBs is accompanied by a modification of its secondary structure. The studied IBs also contain amyloid-like structures. The results obtained in this study expand our knowledge of the structural organization of IBs formed by proteins of different structural classes and also have a contribution into the new approaches development of producing functionally active recombinant membrane proteins.

Biotechnological potential and initial characterization of two novel sesquiterpene synthases from Basidiomycota Coniophora puteana for heterologous production of \u03b4-cadinol

BackgroundTerpene synthases are versatile catalysts in all domains of life, catalyzing the formation of an enormous variety of different terpenoid secondary metabolites. Due to their diverse bioactive properties, terpenoids are of great interest as innovative ingredients in pharmaceutical and cosmetic applications. Recent advances in genome sequencing have led to the discovery of numerous terpene synthases, in particular in Basidiomycota like the wood rotting fungus Coniophora puteana, which further enhances the scope for the manufacture of terpenes for industrial purposes.ResultsIn this study we describe the identification of two novel (+)-δ-cadinol synthases from C. puteana, Copu5 and Copu9. The sesquiterpene (+)-δ-cadinol was previously shown to exhibit cytotoxic activity therefore having an application as possible, new, and sustainably sourced anti-tumor agent. In an Escherichia coli strain, optimized for sesquiterpene production, titers of 225 mg l−1 and 395 mg l−1, respectively, could be achieved. Remarkably, both enzymes share the same product profile thereby representing the first two terpene synthases from Basidiomycota with identical product profiles. We solved the crystal structure of Copu9 in its closed conformation, for the first time providing molecular details of sesquiterpene synthase from Basidiomycota. Based on the Copu9 structure, we conducted structure-based mutagenesis of amino acid residues lining the active site, thereby altering the product profile. Interestingly, the mutagenesis study also revealed that despite the conserved product profiles of Copu5 and Copu9 different conformational changes may accompany the catalytic cycle of the two enzymes. This observation suggests that the involvement of tertiary structure elements in the reaction mechanism(s) employed by terpene synthases may be more complex than commonly expected.ConclusionThe presented product selectivity and titers of Copu5 and Copu9 may pave the way towards a sustainable, biotechnological production of the potentially new bioactive (+)-δ-cadinol. Furthermore, Copu5 and Copu9 may serve as model systems for further mechanistic studies of terpenoid catalysis.Graphical

Macromolecular Crowding Induces a Binding Competent Transient Structure in Intrinsically Disordered Gab1

Intrinsically disordered proteins (IDPs) are an important class of proteins which lack tertiary structure elements. Their dynamic properties can depend on reversible post-translational modifications and the complex cellular milieu, which provides a crowded environment. Both influences the thermodynamic stability and folding of globular proteins as well as the conformational plasticity of IDPs. Here we investigate the intrinsically disordered C-terminal region (amino acids 613–694) of human Grb2-associated binding protein 1 (Gab1), which binds to the disease-relevant Src homolog region 2 (SH2) domain-containing protein tyrosine phosphatase SHP2 (PTPN11). This binding is mediated by phosphorylation at Tyr 627 and Tyr 659 in Gab1. We characterize induced structure in Gab1613–694 and binding to SHP2 by NMR, CD and ITC under non-crowding and crowding conditions, employing chemical and biological crowding agents and compare the results of the non-phosphorylated and tyrosine phosphorylated C-terminal Gab1 fragment. Our results show that under crowding conditions pre-structured motifs in two distinct regions of Gab1 are formed whereas phosphorylation has no impact on the dynamics and IDP character. These structured regions are identical to the binding regions towards SHP2. Therefore, biological crowders could induce some SHP2 binding capacity. Our results therefore indicate that high concentrations of macromolecules stabilize the preformed or excited binding state in the C-terminal Gab1 region and foster the binding to the SH2 tandem motif of SHP2, even in the absence of tyrosine phosphorylation.

Witnessing the structural evolution of an RNA enzyme.

An RNA polymerase ribozyme that has been the subject of extensive directed evolution efforts has attained the ability to synthesize complex functional RNAs, including a full-length copy of its own evolutionary ancestor. During the course of evolution, the catalytic core of the ribozyme has undergone a major structural rearrangement, resulting in a novel tertiary structural element that lies in close proximity to the active site. Through a combination of site-directed mutagenesis, structural probing, and deep sequencing analysis, the trajectory of evolution was seen to involve the progressive stabilization of the new structure, which provides the basis for improved catalytic activity of the ribozyme. Multiple paths to the new structure were explored by the evolving population, converging upon a common solution. Tertiary structural remodeling of RNA is known to occur in nature, as evidenced by the phylogenetic analysis of extant organisms, but this type of structural innovation had not previously been observed in an experimental setting. Despite prior speculation that the catalytic core of the ribozyme had become trapped in a narrow local fitness optimum, the evolving population has broken through to a new fitness locale, raising the possibility that further improvement of polymerase activity may be achievable.

Conservation in the Iron Responsive Element Family.

Iron responsive elements (IREs) are mRNA stem-loop targets for translational control by the two iron regulatory proteins IRP1 and IRP2. They are found in the untranslated regions (UTRs) of genes that code for proteins involved in iron metabolism. There are ten “classic” IRE types that define the conserved secondary and tertiary structure elements necessary for proper IRP binding, and there are 83 published “IRE-like” sequences, most of which depart from the established IRE model. Here are structurally-guided discussions regarding the essential features of an IRE and what is important for IRE family membership.

Apolar chemical environments compact unfolded RNAs and can promote folding

It is well documented that the structure, and thus function, of nucleic acids depends on the chemical environment surrounding them, which often includes potential proteinaceous binding partners. The nonpolar amino acid side chains of these proteins will invariably alter the polarity of the local chemical environment around the nucleic acid. However, we are only beginning to understand how environmental polarity generally influences the structural and energetic properties of RNA folding. Here, we use a series of aqueous-organic cosolvent mixtures to systematically modulate the solvent polarity around two different RNA folding constructs that can form either secondary or tertiary structural elements. Using single-molecule Förster resonance energy transfer spectroscopy to simultaneously monitor the structural and energetic properties of these RNAs, we show that the unfolded conformations of both model RNAs become more compact in apolar environments characterized by dielectric constants less than that of pure water. In the case of tertiary structure formation, this compaction also gives rise to more energetically favorable folding. We propose that these physical changes arise from an enhanced accumulation of counterions in the low dielectric environment surrounding the unfolded RNA.

The effect of trehalose on the thermodynamic stability and emulsification of soybean 11S globulin in the molten globule state

The molten globule (MG) state of soybean glycinin (11S) is partially folded at pH 1.0, so various polyhydroxyl osmotic agents are often used to maintain the structural stability of the MG state under acidic conditions. In this study trehalose was used as a penetrant, and intrinsic fluorescence spectroscopy and circular dichroism were performed to investigate the effects of a gradient concentration of trehalose on the structure of soybean 11S. The MG soybean 11S maintained a more α-helix structure (94.7%) and a larger hydrophobic area when treated with a high concentration of trehalose (1.3 mol/L) than when untreated, and retained a typical compact MG structure despite the lack of tertiary structural elements. A linear extrapolation model showed that the Gibbs free energy of the high-concentration trehalose-induced compact MG state of soybean 11S was reduced by 3 times, compared with the untreated globulin. This improvement of thermodynamic stability indicated that trehalose stabilized this structure by forming a stable trehalose–soybean 11S MG state system. With the addition of gradient trehalose, the surface tension of the solution was significantly reduced and the viscosity of the solution was increased. The emulsifying properties were optimized when the trehalose concentration was 0.5 mol/L, the emulsifying activity (EAI) and emulsifying stability (ESI) were 76.67 m2/g and 76.67%. The improvement of emulsifying property of 11S MG-trehalose system was further proved by observing the microstructure of emulsion by Cryo-SEM. This study advances the understanding of the 11S MG-trehalose system and provides a theoretical basis to produce proteins with excellent emulsifying properties.

Sensing of Proteins by ICD Response of Iron(II) Clathrochelates Functionalized by Carboxyalkylsulfide Groups.

Recognition of elements of protein tertiary structure is crucial for biotechnological and biomedical tasks; this makes the development of optical sensors for certain protein surface elements important. Herein, we demonstrated the ability of iron(II) clathrochelates (1–3) functionalized with mono-, di- and hexa-carboxyalkylsulfide to induce selective circular dichroism (CD) response upon binding to globular proteins. Thus, inherently CD-silent clathrochelates revealed selective inducing of CD spectra when binding to human serum albumin (HSA) (1, 2), beta-lactoglobuline (2) and bovine serum albumin (BSA) (3). Hence, functionalization of iron(II) clathrochelates with the carboxyalkylsulfide group appears to be a promising tool for the design of CD-probes sensitive to certain surface elements of proteins tertiary structure. Additionally, interaction of 1–3 with proteins was also studied by isothermal titration calorimetry, protein fluorescence quenching, electrospray ionization mass spectrometry (ESI-MS) and computer simulations. Formation of both 1:1 and 1:2 assemblies of HSA with 1–3 was evidenced by ESI-MS. A protein fluorescence quenching study suggests that 3 binds with both BSA and HSA via the sites close to Trp residues. Molecular docking calculations indicate that for both BSA and HSA, binding of 3 to Site I and to an “additional site” is more favorable energetically than binding to Site II.

How fast are the motions of tertiary-structure elements in proteins?

Protein motions occur on multiple time and distance scales. Large-scale motions of protein tertiary-structure elements, i.e., domains, are particularly intriguing as they are essential for the catalytic activity of many enzymes and for the functional cycles of protein machines and motors. Theoretical estimates suggest that domain motions should be very fast, occurring on the nanosecond or microsecond time scales. Indeed, free-energy barriers for domain motions are likely to involve salt bridges, which can break in microseconds. Experimental methods that can directly probe domain motions on fast time scales have appeared only in recent years. This Perspective discusses briefly some of these techniques, including nuclear magnetic resonance and single-molecule fluorescence spectroscopies. We introduce a few recent studies that demonstrate ultrafast domain motions and discuss their potential roles. Particularly surprising is the observation of tertiary-structure element dynamics that are much faster than the functional cycles in some protein machines. These swift motions can be rationalized on a case-by-case basis. For example, fast domain closure in multi-substrate enzymes may be utilized to optimize relative substrate orientation. Whether a large mismatch in time scales of conformational dynamics vs functional cycles is a general design principle in proteins remains to be determined.

Inclusion Bodies of Recombinant OmpF Porin from Yersinia pseudotuberculosis: Properties and Structural Characterization.

Mature pore-forming OmpF protein from the outer membrane of Yersinia pseudotuberculosis was expressed in Escherichia coli in the form of inclusion bodies (IBs) under different cultivation conditions. The properties and structural organization of the IBs as well as the structure of the recombinant porin (rOmpF) solubilized from the IBs were investigated using electron microscopy, dynamic light scattering, optical spectroscopy, and specific hydrophobic dyes. The size, shape, and stability of the IBs under denaturing solutions were determined. It was found that the IBs were readily soluble in SDS and more resistant to urea. Dissolution of the IBs in both denaturing agents led to formation of a heterogeneous in size population of oligomeric particles. The IBs contained an intermediate form of the rOmpF with native-like secondary structure and elements of tertiary structure, which was able to penetrate a lipid bilayer and adopt a functionally active conformation. There were no significant differences in the properties and structure between the examined IBs formed at different concentrations of the inducer (IPTG). However, the content of amyloids in the IBs increased with increasing concentration of the inducer. These results contribute to the development of new approaches for the production of active proteins from IBs, as well as biologically and functionally active IBs.

NMR resonance assignments for the GTP-binding RNA aptamer 9-12 in complex with GTP.

Ligand binding RNAs such as artificially created RNA-aptamers are structurally highly diverse. Therefore, they represent important model systems for investigating RNA-folding, RNA-dynamics and the molecular recognition of chemically very different ligands, ranging from small molecules to whole cells. High-resolution structures of RNA-aptamers in complex with their cognate ligands often reveal unexpected tertiary structure elements. Recent studies on different classes of aptamers binding the nucleotide triphosphate GTP as a ligand showed that these systems not only differ widely in binding affinity but also in their ligand binding modes and structural complexity. We initiated the NMR-based structure determination of the high-affinity binding GTP-aptamer 9-12 in order to gain further insights into the diversity of ligand binding modes and structural variability of those aptamers. Here, we report 1H, 13C and 15N resonance assignments for the GTP 9-12-aptamer bound to GTP as the prerequisite for the structure determination by solution NMR.

Direct visualization of the native structure of viroid RNAs at single-molecule resolution by atomic force microscopy

ABSTRACTViroids are small infectious, non-protein-coding circular RNAs that replicate independently and, in some cases, incite diseases in plants. They are classified into two families: Pospiviroidae, composed of species that have a central conserved region (CCR) and replicate in the cell nucleus, and Avsunviroidae, containing species that lack a CCR and whose multimeric replicative intermediates of either polarity generated in plastids self-cleave through hammerhead ribozymes. The compact, rod-like or branched, secondary structures of viroid RNAs have been predicted by RNA folding algorithms and further examined using different in vitro and in vivo experimental techniques. However, direct data about their native tertiary structure remain scarce. Here we have applied atomic force microscopy (AFM) to image at single-molecule resolution different variant RNAs of three representative viroids: potato spindle tuber viroid (PSTVd, family Pospiviroidae), peach latent mosaic viroid and eggplant latent viroid (PLMVd and ELVd, family Avsunviroidae). Our results provide a direct visualization of their native, three-dimensional conformations at 0 and 4 mM Mg2+ and highlight the role that some elements of tertiary structure play in their stabilization. The AFM images show that addition of 4 mM Mg2+ to the folding buffer results in a size contraction in PSTVd and ELVd, as well as in PLMVd when the kissing-loop interaction that stabilizes its 3D structure is preserved.

Synthetic bPNAs as allosteric triggers of hammerhead ribozyme catalysis.

The biochemistry and structural biology of the hammerhead ribozyme (HHR) have been well elucidated. The secondary and tertiary structural elements that enable sugar-phosphate bond scission to be catalyzed by this RNA are clearly understood. We have taken advantage of this knowledge base to test the extent to which synthetic molecules, may be used to trigger structure in secondary structure and tertiary interactions and thereby control HHR catalysis. These molecules belong to a family of molecules we generally call "bPNAs" based on our work on bifacial peptide nucleic acid (bPNA). This family of molecules displays the "bifacial" heterocycle melamine, which acts as a base-triple upon capturing two equivalents of thymine or uracil. Loosely structured internal oligouridylate bulges of 4-20 nucleotides can be restructured as triplex hybrid stems upon binding bPNAs. As such, a duplex stem element can be replaced with a bPNA triplex hybrid stem; similarly, a tertiary loop-stem interaction can be replaced with a loop-bPNA-stem complex. The ability to control RNA structure-function facilitates elucidation of these critical aspects of RNA recognition. In this chapter, we discuss how bPNAs are prepared and applied to study structure-function turn on in the hammerhead ribozyme system.

Physicochemical characterization of E. coli-derived human serum albumin and its comparison with the human plasma counterpart reveals it as a promising biosimilar

Human serum albumin one of the most demanded proteins possess an array of clinical and biotechnological applications. Currently, the prime source for HSA production is the human blood which possesses the risk of pathogen contamination and is limited. Thus, there exists an indispensable need to promote non-animal derived HSA production. In the present work, we have exploited the opportunity and promoted the preparation of pathogen-free rHSA from the E. coli host which is blessed with numerous advantages like scalability, cost-effectiveness etc. Upon overcoming the difficulties to produce functional rHSA in E. coli, through engineering the biological system of protein folding in the cell, the E. coli-derived rHSA has been purified to homogeneity. Its detailed physicochemical characterization has been performed, by monitoring its conformational properties, secondary and tertiary structure elements, surface properties, ligand binding properties, stability issues etc. These parameters of the recombinant protein have been compared with the naturally occurring protein from the human source. The outcome of the comparison reveals that the recombinant protein resembles exactly the same as the natural one. Hence, we propose and promote that the E. coli-derived rHSA is an ideal biosimilar for human blood plasma-derived serum albumin.

Single-Molecule Fluorescence Reveals Commonalities and Distinctions among Natural and in Vitro-Selected RNA Tertiary Motifs in a Multistep Folding Pathway.

Decades of study of the RNA folding problem have revealed that diverse and complex structured RNAs are built from a common set of recurring structural motifs, leading to the perspective that a generalizable model of RNA folding may be developed from understanding of the folding properties of individual structural motifs. We used single-molecule fluorescence to dissect the kinetic and thermodynamic properties of a set of variants of a common tertiary structural motif, the tetraloop/tetraloop-receptor (TL/TLR). Our results revealed a multistep TL/TLR folding pathway in which preorganization of the ubiquitous AA-platform submotif precedes the formation of the docking transition state and tertiary A-minor hydrogen bond interactions form after the docking transition state. Differences in ion dependences between TL/TLR variants indicated the occurrence of sequence-dependent conformational rearrangements prior to and after the formation of the docking transition state. Nevertheless, varying the junction connecting the TL/TLR produced a common kinetic and ionic effect for all variants, suggesting that the global conformational search and compaction electrostatics are energetically independent from the formation of the tertiary motif contacts. We also found that in vitro-selected variants, despite their similar stability at high Mg2+ concentrations, are considerably less stable than natural variants under near-physiological ionic conditions, and the occurrence of the TL/TLR sequence variants in Nature correlates with their thermodynamic stability in isolation. Overall, our findings are consistent with modular but complex energetic properties of RNA structural motifs and will aid in the eventual quantitative description of RNA folding from its secondary and tertiary structural elements.

Identification and characterization of a thermostable endolytic β-agarase Aga2 from a newly isolated marine agarolytic bacteria Cellulophaga omnivescoria W5C

Research on the enzymatic breakdown of seaweed-derived agar has recently gained attention due to the progress in green technologies for marine biomass utilization. The enzymes known as agarases catalyze the cleavage of glycosidic bonds within the polysaccharide. In this study, a new β-agarase, Aga2, was identified from Cellulophaga omnivescoria W5C. Aga2 is one of four putative agarases from the W5C genome, and it belongs to the glycoside hydrolase 16 family. It was shown to be exclusive to the Cellulophaga genus. Agarase activity assays showed that Aga2 is an endolytic-type β-agarase that produces tetrameric and hexameric neoagaro-oligosaccharides, with optimum activity at 45°C and pH 8.0. Zinc ions slightly enhanced its activity while manganese ions had inhibitory effects even at very low concentrations. Aga2 has a Km of 2.59mgmL−1 and Vmax of 275.48Umg−1. The Kcat is 1.73×102s−1, while the Kcat/Km is 8.04×106s−1M−1. Aga2 also showed good thermostability at 45°C and above, and retained >90% of its activity after repeated freeze-thaw cycles. Bioinformatic analysis of its amino acid sequence revealed that intrinsic properties of the protein (e.g. presence of certain dipeptides and the relative volume occupied by aliphatic amino acids) and tertiary structural elements (e.g. presence of salt bridges, hydrophobic interactions and H-bonding) contributed to its thermostability.

Recognizing RNA structural motifs in HT-SELEX data for ribosomal protein S15

BackgroundProteins recognize many different aspects of RNA ranging from single stranded regions to discrete secondary or tertiary structures. High-throughput sequencing (HTS) of in vitro selected populations offers a large scale method to study RNA-proteins interactions. However, most existing analysis methods require that the binding motifs are enriched in the population relative to earlier rounds, and that motifs are found in a loop or single stranded region of the potential RNA secondary structure. Such methods do not generalize to all RNA-protein interaction as some RNA binding proteins specifically recognize more complex structures such as double stranded RNA.ResultsIn this study, we use HT-SELEX derived populations to study the landscape of RNAs that interact with Geobacillus kaustophilus ribosomal protein S15. Our data show high sequence and structure diversity and proved intractable to existing methods. Conventional programs identified some sequence motifs, but these are found in less than 5-10% of the total sequence pool. Therefore, we developed a novel framework to analyze HT-SELEX data. Our process accounts for both sequence and structure components by abstracting the overall secondary structure into smaller substructures composed of a single base-pair stack, which allows us to leverage existing approaches already used in k-mer analysis to identify enriched motifs. By focusing on secondary structure motifs composed of specific two base-pair stacks, we identified significantly enriched or depleted structure motifs relative to earlier rounds.ConclusionsDiscrete substructures are likely to be important to RNA-protein interactions, but they are difficult to elucidate. Substructures can help make highly diverse sequence data more tractable. The structure motifs provide limited accuracy in predicting enrichment suggesting that G. kaustophilus S15 can either recognize many different secondary structure motifs or some aspects of the interaction are not captured by the analysis. This highlights the importance of considering secondary and tertiary structure elements and their role in RNA-protein interactions.