ddbRNA: detection of conserved secondary structures in multiple alignments

MOTIVATION: Structured non-coding RNAs (ncRNAs) have a very important functional role in the cell. No distinctive general features common to all ncRNA have yet been discovered. This makes it difficult to design computational tools able to detect novel ncRNAs in the genomic sequence. RESULTS: We devised an algorithm able to detect conserved secondary structures in both pairwise and multiple DNA sequence alignments with computational time proportional to the square of the sequence length. We implemented the algorithm for the case of pairwise and three-way alignments and tested it on ncRNAs obtained from public databases. On the test sets, the pairwise algorithm has a specificity greater than 97% with a sensitivity varying from 22.26% for Blast alignments to 56.35% for structural alignments. The three-way algorithm behaves similarly. Our algorithm is able to efficiently detect a conserved secondary structure in multiple alignments.     

Identifying potential tRNA genes in genomic DNA sequences.

We have developed an algorithm that automatically and reproducibly identifies potential tRNA genes in genomic DNA sequences, and we present a general strategy for testing the sensitivity of such algorithms. This algorithm is useful for the flagging and characterization of long genomic sequences that have not been experimentally analyzed for identification of functional regions, and for the scanning of nucleotide sequence databases for errors in the sequences and the functional assignments associated with them. In an exhaustive scan of the GenBank database, 97.5% of the 744 known tRNA genes were correctly identified (true-positives), and 42 previously unidentified sequences were predicted to be tRNAs. A detailed analysis of these latter predictions reveals that 16 of the 42 are very similar to known tRNA genes, and we predict that they do, in fact, code for tRNA, yielding a false-positive rate for the algorithm of 0.003%. The new algorithm and testing strategy are a considerable improvement over any previously described strategies for recognizing tRNA genes, and they allow detections of genes (including introns) embedded in long genomic sequences.     

The role of RuvA octamerization for RuvAB function in vitro and in vivo

RuvA plays an essential role in branch migration of the Holliday junction by RuvAB as part of the RuvABC pathway for processing Holliday junctions in Escherichia coli. Two types of RuvA-Holliday junction complexes have been characterized: 1) complex I containing a single RuvA tetramer and 2) complex II in which the junction is sandwiched between two RuvA tetramers. The functional differences between the two forms are still not clear. To investigate the role of RuvA octamerization, we introduced three amino acid substitutions designed to disrupt the E. coli RuvA tetramer-tetramer interface as identified by structural studies. The mutant RuvA was tetrameric and interacted with both RuvB and junction DNA but, as predicted, formed complex I only at protein concentrations up to 500 nm. We present biochemical and surface plasmon resonance evidence for functional and physical interactions of the mutant RuvA with RuvB and RuvC on synthetic junctions. The mutant RuvA with RuvB showed DNA helicase activity and could support branch migration of synthetic four-way and three-way junctions. However, junction binding and the efficiency of branch migration of four-way junctions were affected. The activity of the RuvA mutant was consistent with a RuvAB complex driven by one RuvB hexamer only and lead us to propose that one RuvA tetramer can only support the activity of one RuvB hexamer. Significantly, the mutant failed to complement the UV sensitivity of E. coli DeltaruvA cells. These results indicate strongly that RuvA octamerization is essential for the full biological activity of RuvABC.     

Blocking gap junctional intercellular communication in myoblasts inhibits myogenin and MRF4 expression

Cells rely heavily on cues from their extracellular environment and other cells to coordinate normal physiological processes, and the exchange of molecules via gap junctions has been suggested as an important avenue for cell-cell communication. Gap junctions are found in virtually all mammalian tissues with the notable exception of adult skeletal muscle. However, since functional gap junctions have been detected during the early stages of muscle development, gap junctional intercellular communication (GJIC) may play an important role in myogenesis. In this study, GJIC in normal L6 myoblasts was inhibited using the known blockers 1-octanal and β-glycyrrhetinic acid (β-GA). Under differentiation promoting conditions, L6 cells fused to form multinucleated myotubes, but when treated with either octanol or β-GA, no fusion was observed. The expression of two muscle regulatory factors (MRFs), myogenin and MRF4, was examined in both the blocked and control cells. As expected, the activation of both the myogenin and MRF4 genes coincided with the onset of differentiation in the control L6 cells. Neither of these genes were turned on in the blocked cells, even when grown under low serum conditions. This inhibition of differentiation by octanol and β-GA was reversible, since the activation of both MRF genes as well as myoblast fusion were observed when the blocking medium was replaced with normal differentiating medium. These results suggest that intercellular communication via gap junctions plays an important role in skeletal muscle development and perhaps in the cell signaling events that trigger the activation of muscle-specific MRF genes. Dev. Genet. 20:133–144. 1997. © 1997 Wiley-Liss, Inc.     

Stability and structure of three-way DNA junctions containing unpaired nucleotides.

Non-paired nucleotides stabilize the formation of three-way helical DNA junctions. Two or more unpaired nucleotides located in the junction region enable oligomers ten to fifteen nucleotides long to assemble, forming conformationally homogeneous junctions, as judged by native gel electrophoresis. The unpaired bases can be present on the same strand or on two different strands. Up to five extra bases on one strand have been tested and found to produce stable junctions. The formation of stable structures is favored by the presence of a divalent cation such as magnesium and by high monovalent salt concentration. The order-disorder transition of representative three-way junctions was monitored optically in the ultraviolet and analyzed to quantify thermodynamically the stabilization provided by unpaired bases in the junction region. We report the first measurements of the thermodynamics of adding an unpaired nucleotide to a nucleic acid three-way junction. We find that delta delta G degrees (37 degrees C) = +0.5 kcal/mol for increasing the number of unpaired adenosines from two to three. Three-way junctions having reporter arms 40 base-pairs long were also prepared. Each of the three reporter arms contained a unique restriction site 15 base-pairs from the junction. Asymmetric complexes produced by selectively cleaving each arm were analyzed on native gels. Cleavage of the double helical arm opposite the strand having the two extra adenosines resulted in a complex that migrated more slowly than complexes produced by cleavage at either of the other two arms. It is likely that the strand containing the unpaired adenosines is kinked at an acute angle, forming a Y-shaped, rather than a T-shaped junction.     

Structures of bulged three-way DNA junctions.

We have studied a series of three-way DNA junctions containing unpaired bases on one strand at the branch-point of the junctions. The global conformation of the arms of the junctions has been analysed by means of polyacrylamide gel electrophoresis, as a function of conditions. We find that in the absence of added metal ions, all the results for all the junctions can be accounted for by extended structures, with the largest angle being that between the arms defined by the strand containing the extra bases. Upon addition of magnesium (II) or hexamine cobalt (III) ions, the electrophoretic patterns change markedly, indicative of ion-dependent folding transitions for some of the junctions. For the junction lacking the unpaired bases, the three inter-arm angles appear to be quite similar, suggesting an extended structure. However, the addition of unpaired bases permits the three-way junction to adopt a significantly different structure, in which one angle becomes smaller than the other two. These species also exhibit marked protection against osmium addition to thymine bases at the point of strand exchange. These results are consistent with a model in which two of the helical arms undergo coaxial stacking in the presence of magnesium ions, with the third arm defining an angle that depends upon the number of unpaired bases.     

Metabolic pathway reconstruction of eugenol to vanillin bioconversion in Aspergillus niger

Identification of missing genes or proteins participating in the metabolic pathways as enzymes are of great interest. One such class of pathway is involved in the eugenol to vanillin bioconversion. Our goal is to develop an integral approach for identifying the topology of a reference or known pathway in other organism. We successfully identify the missing enzymes and then reconstruct the vanillin biosynthetic pathway in Aspergillus niger. The procedure combines enzyme sequence similarity searched through BLAST homology search and orthologs detection through COG & KEGG databases. Conservation of protein domains and motifs was searched through CDD, PFAM & PROSITE databases. Predictions regarding how proteins act in pathway were validated experimentally and also compared with reported data. The bioconversion of vanillin was screened on UV-TLC plates and later confirmed through GC and GC-MS techniques. We applied a procedure for identifying missing enzymes on the basis of conserved functional motifs and later reconstruct the metabolic pathway in target organism. Using the vanillin biosynthetic pathway of Pseudomonas fluorescens as a case study, we indicate how this approach can be used to reconstruct the reference pathway in A. niger and later results were experimentally validated through chromatography and spectroscopy techniques.     

On-line analysis of gap junctions reveals more efficient electrical than dye coupling between islet cells

Pancreatic beta-cells constitute a well-communicating multicellular network that permits a coordinated and synchronized signal transmission within the islet of Langerhans that is necessary for proper insulin release. Gap junctions are the molecular keys that mediate functional cellular connections, which are responsible for electrical and metabolic coupling in the majority of cell types. Although the role of gap junctions in beta-cell electrical coupling is well documented, metabolic communication is still a matter of discussion. Here, we have addressed this issue by use of a fluorescence recovery after photobleaching (FRAP) approach. This technique has been validated as a reliable and noninvasive approach to monitor functional gap junctions in real time. We show that control pancreatic islet cells did not exchange a gap junction-permeant molecule in either clustered cells or intact islets of Langerhans under conditions that allowed cell-to-cell exchange of current-carrying ions. Conversely, we have detected that the same probe was extensively transferred between islet cells of transgenic mice expressing connexin 32 (Cx32) that have enhanced junctional coupling properties. The results indicate that the electrical coupling of native islet cells is more extensive than dye communication. Dye-coupling domains in islet cells appear more restricted than previously inferred with other methods.     

The three-way DNA junction is a Y-shaped molecule in which there is no helix-helix stacking.

We have studied the structure of a number of three-way DNA junctions that were closely related in sequence to four-way junctions studied previously. We observe that the electrophoretic mobility of the species derived by selective shortening of one arm of a junction are very similar whichever arm is shortened, and that this remains so whether or not magnesium is present in the buffer. This suggests that the angles subtended between the arms of the three-way junctions are similar. All thymine bases located immediately at the junction are reactive to osmium tetroxide, indicating that out-of-plane attack is not prevented by helix-helix stacking, and this is also independent of the presence or absence of metal cations. The results suggest that the three-way junction cannot undergo an ion-induced conformational folding involving helical stacking, but remains fixed in a Y-shaped extended conformation. Thus the three- and four-way junctions are quite different in character in the presence of cations.     

Structure and dynamics of three-way DNA junctions: atomic force microscopy studies

We have used atomic force microscopy (AFM) to study the conformation of three-way DNA junctions, intermediates of DNA replication and recombination. Immobile three-way junctions with one hairpin arm (50, 27, 18 and 7 bp long) and two relatively long linear arms were obtained by annealing two partially homologous restriction fragments. Fragments containing inverted repeats of specific length formed hairpins after denaturation. Three-way junctions were obtained by annealing one strand of a fragment from a parental plasmid with one strand of an inverted repeat-containing fragment, purified from gels, and examined by AFM. The molecules are clearly seen as three-armed molecules with one short arm and two flexible long arms. The AFM analysis revealed two important features of three-way DNA junctions. First, three-way junctions are very dynamic structures. This conclusion is supported by a high variability of the inter-arm angle detected on dried samples. The mobility of the junctions was observed directly by imaging the samples in liquid (AFM in situ). Second, measurements of the angle between the arms led to the conclusion that three-way junctions are not flat, but rather pyramid-like. Non-flatness of the junction should be taken into account in analysis of the AFM data.     

Curation is forever: comparative genomics approaches to functional annotation

The process of deriving accurate annotations for genes in newly sequenced genomes has up to now been based largely on projection of results in one organism to genes in other organisms based on sequence similarity. We are now entering a period in which metabolic reconstructions will play an increasingly significant role. These reconstructions will support a gradual refinement of estimates of function, leading to more consistent and precise understanding. They will become key components in establishing which functional roles have not yet been connected to genes and which aspects of an organism's physiology can be understood in terms of existing functional assignments. As each refinement is made, we incrementally approach the goal of establishing a qualitative understanding of an organism in terms of its inventory of genes.     

The role of the vgf gene and VGF-derived peptides in nutrition and metabolism

Energy homeostasis is a complex physiological function coordinated at multiple levels. The issue of genetic regulation of nutrition and metabolism is attracting increasing interest and new energy homeostasis-regulatory genes are continuously identified. Among these genes, vgf is gaining increasing interest following two observations: (1) VGF-/- mice have a lean and hypermetabolic phenotype; (2) the first VGF-derived peptide involved in energy homeostasis, named TLQP-21, has been identified. The aim of this review will be to discuss the role of the vgf gene and VGF derived peptides in metabolic and nutritional functions. In particular we will: (1) provide a brief overview on the central systems regulating energy homeostasis and nutrition particularly focusing on the melanocortin system; (2) introduce the structure and molecular characteristic of vgf; (3) describe the phenotype of VGF deficient mice; (4) present recent data on the metabolic role of VGF-derived peptides, particularly focusing on one peptide named TLQP-21.     

Coenzyme Q - Biosynthesis and functions

In addition to its role as a component of the mitochondrial respiratory chain and our only lipid-soluble antioxidant synthesized endogenously, in recent years coenzyme Q (CoQ) has been found to have an increasing number of other important functions required for normal metabolic processes. A number of genetic mutations that reduce CoQ biosynthesis are associated with serious functional disturbances that can be eliminated by dietary administration of this lipid, making CoQ deficiencies the only mitochondrial diseases which can be successfully treated at present. In connection with certain other diseases associated with excessive oxidative stress, the level of CoQ is elevated as a protective response. Aging, certain experimental conditions and several human diseases reduce this level, resulting in serious metabolic disturbances. Since dietary uptake of this lipid is limited, up-regulation of its biosynthetic pathway is of considerable clinical interest. One approach for this purpose is administration of epoxidated all-trans polyisoprenoids, which enhance both CoQ biosynthesis and levels in experimental systems.     

Alterations of metabolic genes and metabolites in cancer

Altered metabolic regulation has long been observed in human cancer and broadly used in the clinic for tumor detection. Two recent findings--the direct regulation of metabolic enzymes by frequently mutated cancer genes and frequent mutations of several metabolic enzymes themselves in cancer--have renewed interest in cancer metabolism. Supporting a causative role of altered metabolic enzymes in tumorigenesis, abnormal levels of several metabolites have been found to play a direct role in cancer development. The alteration of metabolic genes and metabolites offer not only new biomarkers for diagnosis and prognosis, but also potential new targets for cancer therapy.     

Sequence-dependent folding of DNA three-way junctions

Three-way DNA junctions can adopt several different conformers, which differ in the coaxial stacking of the arms. These structural variants are often dominated by one conformer, which is determined by the DNA sequence. In this study we have compared several three-way DNA junctions in order to assess how the arrangement of bases around the branch point affects the conformer distribution. The results show that rearranging the different arms, while retaining their base sequences, can affect the conformer distribution. In some instances this generates a structure that appears to contain parallel coaxially stacked helices rather than the usual anti-parallel arrangement. Although the conformer equilibrium can be affected by the order of purines and pyrimidines around the branch point, this is not sufficient to predict the conformer distribution. We find that the folding of three-way junctions can be separated into two groups of dinucleotide steps. These two groups show distinctive stacking properties in B-DNA, suggesting there is a correlation between B-DNA stacking and coaxial stacking in DNA junctions.     

Evidence for widespread genomic methylation in the migratory locust, Locusta migratoria (Orthoptera: Acrididae)

The importance of DNA methylation in mammalian and plant systems is well established. In recent years there has been renewed interest in DNA methylation in insects. Accumulating evidence, both from mammals and insects, points towards an emerging role for DNA methylation in the regulation of phenotypic plasticity. The migratory locust (Locusta migratoria) is a model organism for the study of phenotypic plasticity. Despite this, there is little information available about the degree to which the genome is methylated in this species and genes encoding methylation machinery have not been previously identified. We therefore undertook an initial investigation to establish the presence of a functional DNA methylation system in L. migratoria. We found that the migratory locust possesses genes that putatively encode methylation machinery (DNA methyltransferases and a methyl-binding domain protein) and exhibits genomic methylation, some of which appears to be localised to repetitive regions of the genome. We have also identified a distinct group of genes within the L. migratoria genome that appear to have been historically methylated and show some possible functional differentiation. These results will facilitate more detailed research into the functional significance of DNA methylation in locusts.