Additional Secondary Structure Research Articles

Integrated positron emission particle tracking (PEPT) and X-ray computed tomography (CT) imaging of flow phenomena in twisted tape swirl flow

A combined positron emission particle tracking (PEPT) and X-ray computed tomography (CT) technique is presented, and its utility is demonstrated through investigation of flow in a pipe with twisted tape swirl insert with varying flow conditions (diameter-based Reynolds numbers 16,300–63,300). A description of this technique is given, as well as data handling practices used to relate geometric information captured by CT to fluid flow data gathered via PEPT. It is found that the CT component is readily capable of capturing the stainless steel insert geometry in this present system, but the use of combined plastic and metal materials leads to artifacts in imaging of the plastic surface. Nonetheless, CT data are related to PEPT flow measurements, and average velocity fields are calculated via a pseudo-framing and interpolation scheme and used to visualize and interrogate key flow phenomena within the system. Radial velocity profiles of the mean flow characteristics are seen to collapse to a nearly common form across all flow conditions considered. Helical vortices are seen propagating through the flow field, generated by bypass flow around the gap between the insert and pipe wall, with additional coherent secondary flow structures seen in the higher Reynolds number cases. These findings enhance the understanding of the mixing mechanisms in these swirl flows and encourage the continued development of PEPT-CT methodologies for 3D flow measurements in optically inaccessible systems.

Spatial analysis of gold occurrences in the Zuru Schist Belt, northwest Nigeria

ABSTRACT The Zuru Schist Belt in Northwest Nigeria is one of the auriferous schist belts occurring in the western half of Nigeria. The gold mineralization largely hosted within Neoproterozoic metasedimentary rocks is considered orogenic. In an attempt to understand the pattern of occurrence of gold mineralization, their spatial distribution and association with regional structures, spatial analysis including point pattern, fractal, Fry, and distance distribution analyses were carried out to determine and quantify spatial pattern, spatial association, and regional structural control. Point pattern analysis reveals that the gold occurrences in the Zuru Schist Belt have a clustered pattern. Fractal analysis shows that the structures that controlled the distribution of gold occurrences can be divided into two categories: those that operated at a spatial scale of below 8 km, and those that operated at around 8 km and more. Fry analysis shows that the gold occurrences in the Zuru Schist Belt are aligned in a NNE-NE direction. This is parallel with the NNE-NE striking faults, establishing these faults as the main controlling structures. Fry plot of orientations of translated pairs within an 8 km distance shows that in addition to NNE-NE striking faults, folds with ESE trending axial trace and NW striking faults also exert some control on these known gold occurrences. The influence of these structures on the gold occurrences is further established by the result of distance distribution analysis. The results of spatial analysis using point pattern, fractal, Fry, and distance distribution analyses established the NNE-NE striking faults as the major structures that controlled the distribution of gold occurrences in the belt, and NW striking faults as important structures at a local scale. In addition to NW striking faults, this study further identified folds with axial trace oriented in ESE as additional secondary control structures for gold mineralization in the Zuru Schist Belt.

The structural landscape and diversity of Pyricularia oryzae MAX effectors revisited.

Magnaporthe AVRs and ToxB-like (MAX) effectors constitute a family of secreted virulence proteins in the fungus Pyricularia oryzae (syn. Magnaporthe oryzae), which causes blast disease on numerous cereals and grasses. In spite of high sequence divergence, MAX effectors share a common fold characterized by a ß-sandwich core stabilized by a conserved disulfide bond. In this study, we investigated the structural landscape and diversity within the MAX effector repertoire of P. oryzae. Combining experimental protein structure determination and in silico structure modeling we validated the presence of the conserved MAX effector core domain in 77 out of 94 groups of orthologs (OG) identified in a previous population genomic study. Four novel MAX effector structures determined by NMR were in remarkably good agreement with AlphaFold2 (AF2) predictions. Based on the comparison of the AF2-generated 3D models we propose a classification of the MAX effectors superfamily in 20 structural groups that vary in the canonical MAX fold, disulfide bond patterns, and additional secondary structures in N- and C-terminal extensions. About one-third of the MAX family members remain singletons, without strong structural relationship to other MAX effectors. Analysis of the surface properties of the AF2 MAX models also highlights the high variability within the MAX family at the structural level, potentially reflecting the wide diversity of their virulence functions and host targets.

Switchable tetraplex elements in the heterogeneous nuclear ribonucleoprotein K promoter: micro-environment dictated structural transitions of G/C rich elements

We have elucidated the hnRNP K promoter as a hotspot for tetraplex-based molecular switches receptive to micro-environmental stimuli. We have characterised the structural features of four tetraplex-forming loci and identified them as binding sites of transcription factors. These segments form either G-quadruplex or i-motif structures, the structural dynamicity of which has been studied in depth via several biophysical techniques. The tetraplexes display high dynamicity and are influenced by both pH and KCl concentrations in vitro. The loci complementary to these sequences form additional non-canonical secondary structures. In the cellular context, the most eminent observation of this study is the binding of hnRNP K to the i-motif forming sequences in its own promoter. We are the first to report a probable transcriptional autoregulatory function of hnRNP K in coordination with higher-order DNA structures. Herein, we also report the positive interaction of the endogenous tetraplexes with Sp1, a well-known transcriptional regulator. Treatment with tetraplex-specific small molecule ligands further uncovered G-quadruplexes’ functioning as repressors and i-motifs as activators in this context. Together, our findings strongly indicate the critical regulatory role of the identified tetraplex elements in the hnRNP K promoter. Communicated by Ramaswamy H. Sarma

The Secondary Structure of Potato Spindle Tuber Viroid Determines Its Infectivity in Nicotiana benthamiana.

The function of RNAs is determined by their structure. However, studying the relationship between RNA structure and function often requires altering RNA sequences to modify the structures, which leads to the neglect of the importance of RNA sequences themselves. In our research, we utilized potato spindle tuber viroid (PSTVd), a circular-form non-coding infectious RNA, as a model with which to investigate the role of a specific rod-like structure in RNA function. By generating linear RNA transcripts with different start sites, we established 12 PSTVd forms with different secondary structures while maintaining the same sequence. The RNA secondary structures were predicted using the mfold tool and validated through native PAGE gel electrophoresis after in vitro RNA folding. Analysis using plant infection assays revealed that the formation of a correct rod-like structure is crucial for the successful infection of PSTVd. Interestingly, the inability of PSTVd forms with non-rod-like structures to infect plants could be partially compensated by increasing the amount of linear viroid RNA transcripts, suggesting the existence of additional RNA secondary structures, such as the correct rod-like structure, alongside the dominant structure in the RNA inoculum of these forms. Our study demonstrates the critical role of RNA secondary structures in determining the function of infectious RNAs.

Structural basis of a bi-functional malonyl-CoA reductase (MCR) from the photosynthetic green non-sulfur bacterium Roseiflexus castenholzii.

Malonyl-CoA reductase (MCR) is a NADPH-dependent bi-functional enzyme that performs alcohol dehydrogenase and aldehyde dehydrogenase (CoA-acylating) activities in the N- and C-terminal fragments, respectively. It catalyzes the two-step reduction of malonyl-CoA to 3-hydroxypropionate (3-HP), a key reaction in the autotrophic CO2 fixation cycles of Chloroflexaceae green non-sulfur bacteria and the archaea Crenarchaeota. However, the structural basis underlying substrate selection, coordination, and the subsequent catalytic reactions of full-length MCR is largely unknown. For the first time, we here determined the structure of full-length MCR from the photosynthetic green non-sulfur bacterium Roseiflexus castenholzii (RfxMCR) at 3.35 Å resolution. Furthermore, we determined the crystal structures of the N- and C-terminal fragments bound with reaction intermediates NADP+ and malonate semialdehyde (MSA) at 2.0 Å and 2.3 Å, respectively, and elucidated the catalytic mechanisms using a combination of molecular dynamics simulations and enzymatic analyses. Full-length RfxMCR was a homodimer of two cross-interlocked subunits, each containing four tandemly arranged short-chain dehydrogenase/reductase (SDR) domains. Only the catalytic domains SDR1 and SDR3 incorporated additional secondary structures that changed with NADP+-MSA binding. The substrate, malonyl-CoA, was immobilized in the substrate-binding pocket of SDR3 through coordination with Arg1164 and Arg799 of SDR4 and the extra domain, respectively. Malonyl-CoA was successively reduced through protonation by the Tyr743-Arg746 pair in SDR3 and the catalytic triad (Thr165-Tyr178-Lys182) in SDR1 after nucleophilic attack from NADPH hydrides. IMPORTANCE The bi-functional MCR catalyzes NADPH-dependent reduction of malonyl-CoA to 3-HP, an important metabolic intermediate and platform chemical, from biomass. The individual MCR-N and MCR-C fragments, which contain the alcohol dehydrogenase and aldehyde dehydrogenase (CoA-acylating) activities, respectively, have previously been structurally investigated and reconstructed into a malonyl-CoA pathway for the biosynthetic production of 3-HP. However, no structural information for full-length MCR has been available to illustrate the catalytic mechanism of this enzyme, which greatly limits our capacity to increase the 3-HP yield of recombinant strains. Here, we report the cryo-electron microscopy structure of full-length MCR for the first time and elucidate the mechanisms underlying substrate selection, coordination, and catalysis in the bi-functional MCR. These findings provide a structural and mechanistic basis for enzyme engineering and biosynthetic applications of the 3-HP carbon fixation pathways.

Physics-based approach to extend a de novo TIM barrel with rationally designed helix-loop-helix motifs.

Computational protein design promises the ability to build tailor-made proteins de novo. While a range of de novo proteins have been constructed so far, the majority of these designs have idealized topologies that lack larger cavities which are necessary for the incorporation of small molecule binding sites or enzymatic functions. One attractive target for enzyme design is the TIM-barrel fold, due to its ubiquity in nature and capability to host versatile functions. With the successful de novo design of a 4-fold symmetric TIM barrel, sTIM11, an idealized, minimalistic scaffold was created. In this work, we attempted to extend this de novo TIM barrel by incorporating a helix-loop-helix motif into its βα-loops by applying a physics-based modular design approach using Rosetta. Further diversification was performed by exploiting the symmetry of the scaffold to integrate two helix-loop-helix motifs into the scaffold. Analysis with AlphaFold2 and biochemical characterization demonstrate the formation of additional α-helical secondary structure elements supporting the successful extension as intended.

Targeted Sequencing of Cytokine-Induced PI3K-Related Genes in Ulcerative Colitis, Colorectal Cancer and Colitis-Associated Cancer.

Chronic relapsing inflammatory bowel disease is strongly linked to an increased risk of colitis-associated cancer (CAC). One of the well-known inflammatory carcinogenesis pathways, phosphatidylinositol 3-kinase (PI3K), was identified to be a crucial mechanism in long-standing ulcerative colitis (UC). The goal of this study was to identify somatic variants in the cytokine-induced PI3K-related genes in UC, colorectal cancer (CRC) and CAC. Thirty biopsies (n = 8 long-standing UC, n = 11 CRC, n = 8 paired normal colorectal mucosa and n = 3 CAC) were subjected to targeted sequencing on 13 PI3K-related genes using Illumina sequencing and the SureSelectXT Target Enrichment System. The Genome Analysis Toolkit was used to analyze variants, while ANNOVAR was employed to detect annotations. There were 5116 intronic, 355 exonic, 172 untranslated region (UTR) and 59 noncoding intronic variations detected across all samples. Apart from a very small number of frameshifts, the distribution of missense and synonymous variants was almost equal. We discovered changed levels of IL23R, IL12Rß1, IL12Rß2, TYK2, JAK2 and OSMR in more than 50% of the samples. The IL23R variant in the UTR region, rs10889677, was identified to be a possible variant that might potentially connect CAC with UC and CRC. Additional secondary structure prediction using RNAfold revealed that mutant structures were more unstable than wildtype structures. Further functional research on the potential variants is, therefore, highly recommended since it may provide insight on the relationship between inflammation and cancer risk in the cytokine-induced PI3K pathway.

Structural insights into ring-building motif domains involved in bacterial sporulation

Components of specialized secretion systems, which span the inner and outer membranes in Gram-negative bacteria, include ring-forming proteins whose oligomerization was proposed to be promoted by domains called RBM for “Ring-Building Motifs”. During spore formation in Gram-positive bacteria, a transport system called the SpoIIIA-SpoIIQ complex also assembles in the double membrane that surrounds the forespore following its endocytosis by the mother cell. The presence of RBM domains in some of the SpoIIIA proteins led to the hypothesis that they would assemble into rings connecting the two membranes and form a conduit between the mother cell and forespore. Among them, SpoIIIAG forms homo-oligomeric rings in vitro but the oligomerization of other RBM-containing SpoIIIA proteins, including SpoIIIAH, remains to be demonstrated. In this work, we identified RBM domains in the YhcN/YlaJ family of proteins that are not related to the SpoIIIA-SpoIIQ complex. We solved the crystal structure of YhcN from Bacillus subtilis, which confirmed the presence of a RBM fold, flanked by additional secondary structures. As the protein did not show any oligomerization ability in vitro, we investigated the structural determinants of ring formation in SpoIIIAG, SpoIIIAH and YhcN. We showed that in vitro, the conserved core of RBM domains alone is not sufficient for oligomerization while the β-barrel forming region in SpoIIIAG forms rings on its own. This work suggests that some RBMs might indeed participate in the assembly of homomeric rings but others might have evolved toward other functions.

Impact of C-Terminal Chemistry on Self-Assembled Morphology of Guanosine Containing Nucleopeptides.

Herein, we report the design and characterization of guanosine-containing self-assembling nucleopeptides that form nanosheets and nanofibers. Through spectroscopy and microscopy analysis, we propose that the peptide component of the nucleopeptide drives the assembly into β-sheet structures with hydrogen-bonded guanosine forming additional secondary structures cooperatively within the peptide framework. Interestingly, the distinct supramolecular morphologies are driven not by metal cation responsiveness common to guanine-based materials, but by the C-terminal peptide chemistry. This work highlights the structural diversity of self-assembling nucleopeptides and will help advance the development of applications for these supramolecular guanosine-containing nucleopeptides.

Contributions of the C-terminal domain to poly(A)-specific ribonuclease (PARN) stability and self-association

Poly(A)-specific ribonuclease (PARN) catalyzes the degradation of mRNA poly(A) tail to regulate translation efficiency and mRNA decay in higher eukaryotic cells. The full-length PARN is a multi-domain protein containing the catalytic nuclease domain, the R3H domain, the RRM domain and the C-terminal intrinsically unstructured domain (CTD). The roles of the three well-structured RNA-binding domains have been extensively studied, while little is known about CTD. In this research, the impact of CTD on PARN stability and aggregatory potency was studied by comparing the thermal inactivation and denaturation behaviors of full-length PARN with two N-terminal fragments lacking CTD. Our results showed that K+ induced additional regular secondary structures and enhanced PARN stability against heat-induced inactivation, unfolding and aggregation. CTD prevented PARN from thermal inactivation but promoted thermal aggregation to initiate at a temperature much lower than that required for inactivation and unfolding. Blue-shift of Trp fluorescence during thermal transitions suggested that heat treatment induced rearrangements of domain organizations. CTD amplified the stabilizing effect of K+, implying the roles of CTD was mainly achieved by electrostatic interactions. These results suggested that CTD might dynamically interact with the main body of the molecule and release of CTD promoted self-association via electrostatic interactions.

Functional morpho-anatomy of water-gap complexes in physically dormant seed

AbstractPhysical dormancy (PY) occurs in at least 18 angiosperm plant families and is caused by water-impermeable palisade cells in seed (or fruit) coats. Breaking of PY involves disruption or dislodgement of water-gap structures causing the seeds/fruits to become water permeable (non-dormant). The water-gap region is a morphologically distinct area of the seed or fruit coat that forms a water-gap complex. The location, anatomy, morphology and origin of water-gaps can differ between and even within families and genera. Water-gap structures sense environmental conditions that allow seeds with PY to become permeable just prior to the commencement of conditions favourable for germination and plant establishment. There are three basic water-gap morpho-anatomies characterized by the way the water-gap opens: Type-I, Type-II and Type-III. In Type-I water-gaps, specific kinds of cells pull apart to form a surface opening, while in Type-II a portion of the surface structure is pulled away from adjacent cells, opening the water-gap. Type-III is the least common type and has a circular, plug-like structure that is dislodged, whereby water entry occurs. In addition, water-gap complexes are either simple or compound, depending on whether only a single primary water-gap structure is involved in dormancy release or an additional secondary water-gap structure opens, permitting water entry.

Stem-Loop RNA Hairpins in Giant Viruses: Invading rRNA-Like Repeats and a Template Free RNA.

We examine the hypothesis that de novo template-free RNAs still form spontaneously, as they did at the origins of life, invade modern genomes, contribute new genetic material. Previously, analyses of RNA secondary structures suggested that some RNAs resembling ancestral (t)RNAs formed recently de novo, other parasitic sequences cluster with rRNAs. Here positive control analyses of additional RNA secondary structures confirm ancestral and de novo statuses of RNA grouped according to secondary structure. Viroids with branched stems resemble de novo RNAs, rod-shaped viroids resemble rRNA secondary structures, independently of GC contents. 5′ UTR leading regions of West Nile and Dengue flavivirid viruses resemble de novo and rRNA structures, respectively. An RNA homologous with Megavirus, Dengue and West Nile genomes, copperhead snake microsatellites and levant cotton repeats, not templated by Mimivirus' genome, persists throughout Mimivirus' infection. Its secondary structure clusters with candidate de novo RNAs. The saltatory phyletic distribution and secondary structure of Mimivirus' peculiar RNA suggest occasional template-free polymerization of this sequence, rather than noncanonical transcriptions (swinger polymerization, posttranscriptional editing).

Abstract LB-048: Investigating mechanism of action of oligonucleotides homologous to the telomere overhang in melanoma

Abstract Telomeres have become a popular target for novel cancer therapeutics because activation of telomerase is associated with development of immortality in nearly 90% of all cancers. It has been shown that administration of an oligonucleotide homologous to the telomere overhang called T-oligo (5’-GTTAGGGTTAG-3’) induces potent DNA damage responses (DDRs) in several cancers, including melanoma, and causes minimal toxicity to normal cells. Our study focuses on T-oligo which has the guanine rich character of the telomere 3’ overhang. The folding and stabilization of G-quadruplexes (G4) which are G-rich structures at the end of telomeres may be involved in regulation of telomere length and cell cycle regulation. To date, T-oligo has not been shown to form G4, although evidence suggests it has a G4-mediated mechanism of action. In this study, we investigate T-oligo's ability to form G4 structures in vitro. We hypothesize the additional secondary structure of G4 configured T oligo may confer increased stability, resistance to nucleases in serum and, initiate DDRs. G4 formation is a potential target since G4 presence can cause telomere uncapping leading to DDRs. In our studies, G4 T-oligo was formed by heating various concentrations of T-oligo in KCl solutions to 90°C and incubating at 60°C for 24-36 hours. G4 T-oligo formation was verified using native PAGE and NMR spectroscopy before subjecting it to DNase1 degradation and analyzing with native PAGE. We also analyzed the uptake and anti-proliferative ability of G4 T-oligo using FACS analysis and MTT assay respectively. The change in expression of TRF2 and JNK activity was performed using western blot analysis since T-oligo may induce activation of SAPK/JNK pathway. We further investigated the interaction between T-oligo and TRF2 using immunofluorescence (IF) with FITC labeled T-oligo in MM-AN melanoma cells. Densitometric analysis after nuclease digestion demonstrated a 52.6% decrease in T-oligo and a 20% decrease in G4 T-oligo as well as a 96% increase of single stranded T-oligo within the G4 T-oligo sample after 90 minutes, indicating nuclease resistance of G4-T-oligo. FACS analysis revealed G4 T-oligo and T-oligo had similar cellular uptake but G4 T-oligo was shown to be less effective in inhibiting cell proliferation suggesting the reduced anti-proliferative activity of G4 T-oligo is not likely attributed to cellular uptake. Western blot results showed upregulation of both p-JNK and total JNK by 4.0- and 2.0-fold respectively after 24 hour treatment. JNK has previously been demonstrated to be activated by G4 stabilizing ligands, providing more evidence that T-oligo may induce telomere dysfunction and upregulation of JNK activity. IF showed 80% co-localization of T-oligo and TRF2 in the nucleus suggesting T oligo may induce dissociation of TRF2 from the telomere shelterin complex causing DDRs. In conclusion, formation of G4 species by T-oligo may induce DDRs and future studies may provide insights into its mechanism of action and suggests new avenues for research. Citation Format: Mark Frakes, Luke Wojdyla, Marko Ivancich, Perlina Fortinberry, Pooja Vinay, Priyanka Pandey, Gagan Chhabra, Neelu Puri. Investigating mechanism of action of oligonucleotides homologous to the telomere overhang in melanoma. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr LB-048.

Molecular docking of Glycine max and Medicago truncatula ureases with urea; bioinformatics approaches.

Urease (EC 3.5.1.5) is a nickel-dependent metalloenzyme catalyzing the hydrolysis of urea into ammonia and carbon dioxide. It is present in many bacteria, fungi, yeasts and plants. Most species, with few exceptions, use nickel metalloenzyme urease to hydrolyze urea, which is one of the commonly used nitrogen fertilizer in plant growth thus its enzymatic hydrolysis possesses vital importance in agricultural practices. Considering the essentiality and importance of urea and urease activity in most plants, this study aimed to comparatively investigate the ureases of two important legume species such as Glycine max (soybean) and Medicago truncatula (barrel medic) from Fabaceae family. With additional plant species, primary and secondary structures of 37 plant ureases were comparatively analyzed using various bioinformatics tools. A structure based phylogeny was constructed using predicted 3D models of G. max and M. truncatula, whose crystallographic structures are not available, along with three additional solved urease structures from Canavalia ensiformis (PDB: 4GY7), Bacillus pasteurii (PDB: 4UBP) and Klebsiella aerogenes (PDB: 1FWJ). In addition, urease structures of these species were docked with urea to analyze the binding affinities, interacting amino acids and atom distances in urease-urea complexes. Furthermore, mutable amino acids which could potentially affect the protein active site, stability and flexibility as well as overall protein stability were analyzed in urease structures of G. max and M. truncatula. Plant ureases demonstrated similar physico-chemical properties with 833-878 amino acid residues and 89.39-90.91 kDa molecular weight with mainly acidic (5.15-6.10 pI) nature. Four protein domain structures such as urease gamma, urease beta, urease alpha and amidohydro 1 characterized the plant ureases. Secondary structure of plant ureases also demonstrated conserved protein architecture, with predominantly α-helix and random coil structures. In structure-based phylogeny, plant ureases from G. max, M. truncatula and C. ensiformis were clearly diverged from bacterial ureases of B. pasteurii and K. aerogenes. Glu, Thr, His and Gly were commonly found as interacting residues in most urease-urea docking complexes while Glu was available in all docked structures. Besides, Ala and Arg residues, which are reported in active-site architecture of plant and bacterial ureases were present in G. max urea-urease complex but not present in others. Moreover, Arg435 and Arg437 in M. truncatula and G. max, respectively were identified as highly mutable hotspot residues residing in amidohydro 1 domain of enzyme. In addition, a number of stabilizing residues were predicted upon mutation of these hotspot residues however Cys and Thr made strong implications since they were also found in codon-aligned sequences as substitutions of hotspot residues. Comparative analyses of primary sequence and secondary structure in 37 different plants demonstrated quite conserved natures of ureases in plant kingdom. Structure-based phylogeny indicated the presence of a possible prokaryote-eukaryote split and implicated the subjection of bacterial ureases to heavy selection in prokaryotic evolution compared to plants. Urea-urease docking complexes suggested that different species could share common interacting residues as well as may have some other uncommon residues at species-dependent way. In silico mutation analyses identified mutable amino acids, which were predicted to reside in catalytic site of enzyme therefore mutagenesis at these sites seemed to have adverse effects on enzyme efficiency or function. This study findings will become valuable preliminary resource for future studies to further understand the primary, secondary and tertiary structures of urease sequences in plants as well as it will provide insights about various binding features of urea-urease complexes.

Multiple vortex structures in the wake of a rectangular winglet in ground effect

Patterns of vorticity in the wake of a single rectangular winglet (vortex generator) embedded in a turbulent boundary layer have been studied using Stereoscopic Particle Image Velocimetry (SPIV). The winglet was mounted normally to a flat surface with an angle to the oncoming flow. A parametric study varying the winglet height (constant aspect ratio) and angle has shown, contrary to the common classical single tip-vortex conception, that the wake generally consists of a complex system of multiple vortex structures. The primary vortex has previously been discovered to contain a direct coupling between the axial and the rotational flow. In the current work, even the longitudinal secondary structures detected from measured streamwise vorticity display similar behavior. A regime map depicting the observed stable far wake states of the multiple vortices as a function of winglet height and angle reveals complex patterns of the flow topologies not only with the primary tip vortex, but with the additional secondary structures as well. A bifurcation diagram shows distinct regimes of the various secondary structures as well as how the primary vortex is in some cases significantly affected by their presence. These data should serve as inspiration in the process of generating longitudinal vortices for enhancement of heat and mass transfer in industrial devices since the multiple vortex regimes can help improve the conditions for these exchanges. Further, these results point to a weakness in existing inviscid models not accounting for the possibility of multiple vortical structures in the wake.

Conserved Structural Domains in FoxD4L1, a Neural Forkhead Box Transcription Factor, Are Required to Repress or Activate Target Genes

FoxD4L1 is a forkhead transcription factor that expands the neural ectoderm by down-regulating genes that promote the onset of neural differentiation and up-regulating genes that maintain proliferative neural precursors in an immature state. We previously demonstrated that binding of Grg4 to an Eh-1 motif enhances the ability of FoxD4L1 to down-regulate target neural genes but does not account for all of its repressive activity. Herein we analyzed the protein sequence for additional interaction motifs and secondary structure. Eight conserved motifs were identified in the C-terminal region of fish and frog proteins. Extending the analysis to mammals identified a high scoring motif downstream of the Eh-1 domain that contains a tryptophan residue implicated in protein-protein interactions. In addition, secondary structure prediction programs predicted an α-helical structure overlapping with amphibian-specific Motif 6 in Xenopus, and similarly located α-helical structures in other vertebrate FoxD proteins. We tested functionality of this site by inducing a glutamine-to-proline substitution expected to break the predicted α-helical structure; this significantly reduced FoxD4L1’s ability to repress zic3 and irx1. Because this mutation does not interfere with Grg4 binding, these results demonstrate that at least two regions, the Eh-1 motif and a more C-terminal predicted α-helical/Motif 6 site, additively contribute to repression. In the N-terminal region we previously identified a 14 amino acid motif that is required for the up-regulation of target genes. Secondary structure prediction programs predicted a short β-strand separating two acidic domains. Mutant constructs show that the β-strand itself is not required for transcriptional activation. Instead, activation depends upon a glycine residue that is predicted to provide sufficient flexibility to bring the two acidic domains into close proximity. These results identify conserved predicted motifs with secondary structures that enable FoxD4L1 to carry out its essential functions as both a transcriptional repressor and activator of neural genes.

Statistical torsion angle potential energy functions for protein structure modeling: A bicubic interpolation approach

A set of grid type knowledge-based energy functions is introduced for ϕ-χ₁ , ψ-χ₁ , ϕ-ψ, and χ₁ -χ₂ torsion angle combinations. Boltzmann distribution is assumed for the torsion angle populations from protein X-ray structures, and the functions are named as statistical torsion angle potential energy functions. The grid points around periodic boundaries are duplicated to force periodicity, and the remedy relieves the derivative discontinuity problem. The devised functions rapidly improve the quality of model structures. The potential bias in the functions and the usefulness of additional secondary structure information are also investigated. The proposed guiding functions are expected to facilitate protein structure modeling, such as protein structure prediction, protein design, and structure refinement.

MicroRNA discovery by similarity search to a database of RNA-seq profiles

In silico generated search for microRNAs (miRNAs) has been driven by methods compiling structural features of the miRNA precursor hairpin, as well as to some degree combining this with the analysis of RNA-seq profiles for which the miRNA typically leave the drosha/dicer fingerprint of 1–2 ~22 nt blocks of reads corresponding to the mature and star miRNA. In complement to the previous methods, we present a study where we systematically exploit these patterns of read profiles. We created two datasets comprised of 2540 and 4795 read profiles obtained after preprocessing short RNA-seq data from miRBase and ENCODE, respectively. Out of 4795 ENCODE read profiles, 1361 are annotated as non-coding RNAs (ncRNAs) and of which 285 are further annotated as miRNAs. Using deepBlockAlign (dba), we align ncRNA read profiles from ENCODE against the miRBase read profiles (cleaned for “self-matches”) and are able to separate ENCODE miRNAs from the other ncRNAs by a Matthews Correlation Coefficient (MCC) of 0.8 and obtain an area under the curve of 0.93. Based on the dba score cut-off of 0.7 at which we observed the maximum MCC of 0.8, we predict 523 novel miRNA candidates. An additional RNA secondary structure analysis reveal that 42 of the candidates overlap with predicted conserved secondary structure. Further analysis reveal that the 523 miRNA candidates are located in genomic regions with MAF block (UCSC) fragmentation and poor sequence conservation, which in part might explain why they have been overlooked in previous efforts. We further analyzed known human and mouse miRNA read profiles and found two distinct classes; the first containing two blocks and the second containing >2 blocks of reads. Also the latter class holds read profiles that have less well defined arrangement of reads in comparison to the former class. On comparison of miRNA read profiles from plants and animals, we observed kingdom specific read profiles that are distinct in terms of both length and distribution of reads within the read profiles to each other. All the data, as well as a server to search miRBase read profiles by uploading a BED file, is available at http://rth.dk/resources/mirdba.

RING: networking interacting residues, evolutionary information and energetics in protein structures

Residue interaction networks (RINs) have been used in the literature to describe the protein 3D structure as a graph where nodes represent residues and edges physico-chemical interactions, e.g. hydrogen bonds or van-der-Waals contacts. Topological network parameters can be calculated over RINs and have been correlated with various aspects of protein structure and function. Here we present a novel web server, RING, to construct physico-chemically valid RINs interactively from PDB files for subsequent visualization in the Cytoscape platform. The additional structure-based parameters secondary structure, solvent accessibility and experimental uncertainty can be combined with information regarding residue conservation, mutual information and residue-based energy scoring functions. Different visualization styles are provided to facilitate visualization and standard plugins can be used to calculate topological parameters in Cytoscape. A sample use case analyzing the active site of glutathione peroxidase is presented. The RING server, supplementary methods, examples and tutorials are available for non-commercial use at URL: http://protein.bio.unipd.it/ring/.