Bioequivalence in dogs of a meloxicam formulation administered as a transmucosal oral mist with an orally administered pioneer suspension product

A mucosal mist formulation of meloxicam, administered as a spray into the mouth (test article), was compared for bioequivalence to a pioneer meloxicam suspension for oral administration (reference article). Pharmacokinetic profiles and average bioequivalence were investigated in 20 dogs. The study design comprised a two-period, two-sequence, two-treatment cross-over design, with maximum concentration (C(max)) and area under plasma concentration-time curve to last sampling time (AUC(last)) used as pivotal bioequivalence variables. Bioequivalence of the products was confirmed, based on relative ratios of geometric mean concentrations (and 90% confidence intervals within the range 0.80-1.25) for C(max) of 101.9 (97.99-106.0) and for AUC(last) of 97.24 (94.44-100.1). The initial absorption of meloxicam was more rapid for the test article, despite virtually identical C(max) values for the two products. Mean elimination half-lives were 29.6 h (test article) and 30.0 h (reference article). The meloxicam plasma concentration-time profiles were considered in relation to published data on the inhibition of the cyclooxygenase-1 (COX-1) and COX-2 isoenzymes by meloxicam.

On the evolutionary origin of the plant mitochondrion and its genome

Higher plants occupy very different positions in the mitochondrial and nuclear lineages of global phylogenetic trees based on conserved regions of small subunit (SSU) and large subunit (LSU) rRNA sequences. In the nuclear subtree, plants branch off late, at a position reflecting a massive radiation of the major multicellular (and some unicellular) groups; in the mitochondrial subtree, in contrast, plants branch off early, near the point of connection between the mitochondrial and eubacterial lineages. Moreover, in the nuclear lineage, plants branch together with the unicellular green alga Chlamydomonas reinhardtii, whereas in the mitochondrial lineage (in both SSU and LSU trees), metaphytes and chlorophyte branch separately. Statistical evaluation indicates that the anomalous branching position of higher plants in the mitochondrial lineage is not a treeing artifact attributable to the relatively rapid rate of sequence divergence of non-plant mitochondrial rRNA sequences. In considering alternative biological explanations for these results, we are led to propose that the rRNA genes in plant mitochondria may be of more recent evolutionary origin than the rRNA genes in other mitochondria. This proposal has implications for monophyletic vs. polyphyletic scenarios of mitochondrial origin and is consistent with other evidence indicating that plant mtDNA is an evolutionary mosaic.

Cleavage of mitochondria-like transfer RNAs expressed in Escherichia coli

Mitochondrial (mt) transfer RNAs (tRNAs) often harbor unusual structural features causing their secondary structure to differ from the conventional cloverleaf. tRNAs designed with such irregularities, termed mt-like tRNAs, are active in Escherichia coli as suppressors of reporter genes, although they display low steady-state levels. Characterization of fragments produced during mt-like tRNA processing in vitro and in vivo suggests that these RNAs are not fully processed at their 5' ends and are cleaved internally. These abnormal processing events may account for the low levels of mature mt-like RNAs in vivo and are most likely related to defective processing by RNase P.

Identification and characterization of a Plasmodium falciparum RNA polymerase gene with similarity to mitochondrial RNA polymerases

Nearly all mitochondrial RNA polymerase genes identified to date are encoded in the nucleus and have similarities to T3 and T7 bacteriophage RNA polymerases. Some chloroplast genes are also transcribed by T3/T7 phage-like RNA polymerases, raising the possibility that the apicomplexan parasites, which have both a mitochondrion and a plastid, might have two such genes. As part of an investigation of Plasmodium falciparum organelle transcription, we initiated a search for T3/T7 bacteriophage-like RNA polymerase genes. We employed degenerate primers based on highly conserved plant, animal and fungal mitochondrial RNA polymerase sequences to amplify corresponding P. falciparum sequences by polymerase chain reaction (PCR). Less well-conserved flanking sequences were obtained by inverse PCR. The resulting sequence predicts a 1503 amino acid open reading frame with similarity to other T3/T7 phage-like RNA polymerases. Essential amino acids that have been identified in T7 mutant analyses are conserved in the P. falciparum RNA polymerase gene. Comparison of the sequence with preliminary data from the P. falciparum genome sequencing project revealed strain heterogeneity within two regions of the gene. The amino-terminal predicted amino acid sequence of the RNA polymerase gene has similarities to mitochondrial targeting sequences. Taken together, these points suggest that we have identified the P. falciparum mitochondrial RNA polymerase gene.

Hammerhead-mediated processing of satellite pDo500 family transcripts from Dolichopoda cave crickets

This work reports the discovery and functional characterization of catalytically active hammerhead motifs within satellite DNA of the pDo500 family from several DOLICHOPODA: cave cricket species. We show that in vitro transcribed RNA of some members of this satellite DNA family do self-cleave in vitro. This self-cleavage activity is correlated with the efficient in vivo processing of long primary transcripts into monomer-sized RNA. The high sequence conservation of the satellite pDo500 DNA family among genetically isolated DOLICHOPODA: schiavazzii populations, as well as other DOLICHOPODA: species, along with the fact that satellite members are actively transcribed in vivo suggests that the hammerhead-encoding satellite transcripts are under selective pressure, perhaps because they fulfil an important physiological role or function. Remarkably, this is the third example of hammerhead ribozyme structures associated with transcribed repetitive DNA sequences from animals. The possibility that such an association may not be purely coincidental is discussed.

Distribution of hammerhead and hammerhead-like RNA motifs through the GenBank

Hammerhead ribozymes previously were found in satellite RNAs from plant viroids and in repetitive DNA from certain species of newts and schistosomes. To determine if this catalytic RNA motif has a wider distribution, we decided to scrutinize the GenBank database for RNAs that contain hammerhead or hammerhead-like motifs. The search shows a widespread distribution of this kind of RNA motif in different sequences suggesting that they might have a more general role in RNA biology. The frequency of the hammerhead motif is half of that expected from a random distribution, but this fact comes from the low CpG representation in vertebrate sequences and the bias of the GenBank for those sequences. Intriguing motifs include those found in several families of repetitive sequences, in the satellite RNA from the carrot red leaf luteovirus, in plant viruses like the spinach latent virus and the elm mottle virus, in animal viruses like the hepatitis E virus and the caprine encephalitis virus, and in mRNAs such as those coding for cytochrome P450 oxidoreductase in the rat and the hamster.

The distribution of RNA motifs in natural sequences

Functional analysis of genome sequences has largely ignored RNA genes and their structures. We introduce here the notion of 'ribonomics' to describe the search for the distribution of and eventually the determination of the physiological roles of these RNA structures found in the sequence databases. The utility of this approach is illustrated here by the identification in the GenBank database of RNA motifs having known binding or chemical activity. The frequency of these motifs indicates that most have originated from evolutionary drift and are selectively neutral. On the other hand, their distribution among species and their location within genes suggest that the destiny of these motifs may be more elaborate. For example, the hammerhead motif has a skewed organismal presence, is phylogenetically stable and recent work on a schistosome version confirms its in vivo biological activity. The under-representation of the valine-binding motif and the Rev-binding element in GenBank hints at a detrimental effect on cell growth or viability. Data on the presence and the location of these motifs may provide critical guidance in the design of experiments directed towards the understanding and the manipulation of RNA complexes and activities in vivo.

The recognition of a noncanonical RNA base pair by a zinc finger protein

BACKGROUND

The zinc finger (ZF) is the most abundant nucleic-acid-interacting protein motif. Although the interaction of ZFs with DNA is reasonably well understood, little is known about the RNA-binding mechanism. We investigated RNA binding to ZFs using the Zif268-DNA complex as a model system. Zif268 contains three DNA-binding ZFs; each independently binds a 3 base pair (bp) subsite within a 9 bp recognition sequence.

RESULTS

We constructed a library of phage-displayed ZFs by randomizing the alpha helix of the Zif268 central finger. Successful selection of an RNA binder required a noncanonical base pair in the middle of the RNA triplet. Binding of the Zif268 variant to an RNA duplex containing a G.A mismatch (rG.A) is specific for RNA and is dependent on the conformation of the mismatched middle base pair. Modeling and NMR analyses revealed that the rG.A pair adopts a head-to-head configuration that counterbalances the effect of S-puckered riboses in the backbone. We propose that the structure of the rG.A duplex is similar to the DNA in the original Zif268-DNA complex.

CONCLUSIONS

It is possible to change the specificity of a ZF from DNA to RNA. The ZF motif can use similar mechanisms in binding both types of nucleic acids. Our strategy allowed us to rationalize the interactions that are possible between a ZF and its RNA substrate. This same strategy can be used to assess the binding specificity of ZFs or other protein motifs for noncanconical RNA base pairs, and should permit the design of proteins that bind specific RNA structures.

A small nucleolar RNA:ribozyme hybrid cleaves a nucleolar RNA target in vivo with near-perfect efficiency

A hammerhead ribozyme has been localized to the yeast nucleolus by using the U3 small nucleolar RNA as a carrier. The hybrid small nucleolar RNA:ribozyme, designated a "snorbozyme," is metabolically stable and cleaves a target U3 RNA with nearly 100% efficiency in vivo. This is the most efficient in vivo cleavage reported for a trans-acting ribozyme. A key advantage of the model substrate featured is that a stable, trimmed cleavage product accumulates. This property allows accurate kinetic measurements of authentic cleavage in vivo. The system offers new avenues for developing effective ribozymes for research and therapeutic applications.

tRNA nucleotide 47: an evolutionary enigma

A previous analysis of tRNA sequences suggested a correlation between the absence of a nucleotide at position 47 (nt 47) in the extra loop and the presence of a U13:G22 base pair in the D-stem. We have evaluated the significance of this correlation by determining the in vivo activity of tRNAs containing either a C13:G22 or a U13:G22 pair in tRNA molecules with or without nt 47. Although this correlation might reflect some malfunction of tRNAs lacking nt 47, but containing the C13:G22, assays of the in vivo suppressor activity showed that this tRNA is actually more active than the tRNA with the features found in the database, i.e., a U13:G22 base pair and no nt 47. Moreover, analogous constructs with a GGC anticodon permitted the growth of an Escherichia coli strain deleted for tRNA(Ala)GGC genes equally well. On the other hand, long-term growth experiments with competing E. coli strains harboring the tRNA lacking nt 47, either with the C13:G22 or the U13:G22 base pair demonstrated that the U13:G22 tRNA overtook the C13:G22 strain even when the starting proportion of strains favored the C13:G22 strain. Thus, the preference for the U13:G22 tRNA lacking nt 47 in the sequence database is most likely due to factors that come into play during extended growth or latency rather than to the ability of the tRNA to engage in protein synthesis.

Schistosome satellite DNA encodes active hammerhead ribozymes

Using a computer program designed to search for RNA structural motifs in sequence databases, we have found a hammerhead ribozyme domain encoded in the Smalpha repetitive DNA of Schistosoma mansoni. Transcripts of these repeats are expressed as long multimeric precursor RNAs that cleave in vitro and in vivo into unit-length fragments. This RNA domain is able to engage in both cis and trans cleavage typical of the hammerhead ribozyme. Further computer analysis of S. mansoni DNA identified a potential trans cleavage site in the gene coding for a synaptobrevin-like protein, and RNA transcribed from this gene was efficiently cleaved by the Smalpha ribozyme in vitro. Similar families of repeats containing the hammerhead domain were found in the closely related Schistosoma haematobium and Schistosomatium douthitti species but were not present in Schistosoma japonicum or Heterobilharzia americana, suggesting that the hammerhead domain was not acquired from a common schistosome ancestor.

Modeling active RNA structures using the intersection of conformational space: application to the lead-activated ribozyme

The Pb2+ cleavage of a specific phosphodiester bond in yeast tRNA(Phe) is the classical model of metal-assisted RNA catalysis. In vitro selection experiments have identified a tRNA(Phe) variant, the leadzyme, that is very active in cleavage by Pb2+. We present here a three-dimensional modeling protocol that was used to propose a structure for this ribozyme, and is based on the computation of the intersection of conformational space of sequence variants and the use of chemical modification data. Sequence and secondary structure data were used in a first round of computer modeling that allowed identification of conformations compatible with all known leadzyme variants. Common conformations were then tested experimentally by evaluating the activity of analogues containing modified nucleotides in the catalytic core. These experiments led to a new structural hypothesis that was tested in a second round of computer modeling. The resulting proposal for the active conformation of the leadzyme is consistent with all known structural data. The final model suggests an in-line SN2 attack mechanism and predicts two Pb2+ binding sites. The protocol presented here is generally applicable in modeling RNAs whenever the catalytic or binding activity of structural analogues is known.

SnoRNAs as tools for RNA cleavage and modification

Eukaryotic small nucleolar RNAs (snoRNAs) influence rRNA synthesis at two stages: nucleolytic processing and selection of nucleotides to be ribose methylated (Nm) or converted to pseudouridine (psi). The two modification functions and some processing activities involve direct base pairing of snoRNA with rRNA. In addition to rRNA-targeting sequences, snoRNA function depends on the presence of conserved box elements involved in snoRNA synthesis and localization. The present investigation is directed at using snoRNAs as tools for two purposes: 1) introducing nucleotide modifications into novel sites in rRNA and other snoRNAs, and: 2) targeting nucleolar RNAs for destruction using snoRNA:ribozyme chimers ('snorbozymes'). Early results demonstrate that snoRNAs can be used for both applications.

Amber suppression in Escherichia coli by unusual mitochondria-like transfer RNAs

The "cloverleaf" base-pairing pattern was established as the structural paradigm of active tRNA species some 30 years ago. Nevertheless, this pattern does not accommodate the folding of certain mitochondrial tRNAs. For these recalcitrant tRNAs, we have proposed structures having from 5 to 10 base pairs in the anticodon stem rather than the canonical 6. The absence of these types of tRNAs in cytoplasmic translation systems, however, raises the possibility that they may not be bona fide alternate folding patterns for active tRNA molecules. For this reason, we have designed new tRNA genes based on our model of unusual mitochondrial tRNAs, having 7, 8, 9, and 10 base pairs in the anticodon stem with other modifications to the D-stem and connector regions. We show here that these synthetic genes produce tRNAs that actively suppress amber codons in vivo.

Genome structure and gene content in protist mitochondrial DNAs

Although the collection of completely sequenced mitochondrial genomes is expanding rapidly, only recently has a phylogenetically broad representation of mtDNA sequences from protists (mostly unicellular eukaryotes) become available. This review surveys the 23 complete protist mtDNA sequences that have been determined to date, commenting on such aspects as mitochondrial genome structure, gene content, ribosomal RNA, introns, transfer RNAs and the genetic code and phylogenetic implications. We also illustrate the utility of a comparative genomics approach to gene identification by providing evidence that orfB in plant and protist mtDNAs is the homolog of atp8 , the gene in animal and fungal mtDNA that encodes subunit 8 of the F0portion of mitochondrial ATP synthase. Although several protist mtDNAs, like those of animals and most fungi, are seen to be highly derived, others appear to be have retained a number of features of the ancestral, proto-mitochondrial genome. Some of these ancestral features are also shared with plant mtDNA, although the latter have evidently expanded considerably in size, if not in gene content, in the course of evolution. Comparative analysis of protist mtDNAs is providing a new perspective on mtDNA evolution: how the original mitochondrial genome was organized, what genes it contained, and in what ways it must have changed in different eukaryotic phyla.

Modeling RNA-ligand interactions: the Rev-binding element RNA-aminoglycoside complex

An approach to the modeling of ligand-RNA complexes has been developed by combining three-dimensional structure-activity relationship (3D-SAR) computations with a docking protocol. The ability of 3D-SAR to predict bound conformations of flexible ligands was first assessed by attempting to reconstruct the known, bound conformations of phenyloxazolines complexed with human rhinovirus 14 (HRV14) RNA. Subsequently, the same 3D-SAR analysis was applied to the identification of bound conformations of aminoglycosides which associate with the Rev-binding element (RBE) RNA. Bound conformations were identified by parsing ligand conformational data sets with pharmacophores determined by the 3D-SAR analysis. These "bioactive" structures were docked to the receptor RNA, and optimization of the complex was undertaken by extensive searching of ligand conformational space coupled with molecular dynamics computations. The similarity between the bound conformations of the ligand from the 3D-SAR analysis and those found in the docking protocol suggests that this methodology is valid for the prediction of bound ligand conformations and the modeling of the structure of the ligand-RNA complexes.

On the origin of protein synthesis factors: a gene duplication/fusion model

Sequence similarity has given rise to the proposal that IF-2, EF-G, and EF-Tu are related through a common ancestor. We evaluate this proposition and whether the relationship can be extended to other factors of protein synthesis. Analysis of amino acid sequence similarity gives statistical support for an evolutionary affiliation among IF-1, IF-2, IF-3, EF-Tu, EF-Ts, and EF-G and suggests further that this association is a result of gene duplication/fusion events. In support of this mechanism, the three-dimensional structures of IF-3, EF-Tu, and EF-G display a predictable domain structure and overall conformational similarity. The model that we propose consists of three consecutives duplication/fusion events which would have taken place before the divergence of the three superkingdoms: eubacteria, archaea, and eukaryotes. The root of this protein superfamily tree would be an ancestor of the modern IF-1 gene sequence. The repeated fundamental motif of this protein superfamily is a small RNA binding domain composed of two alpha-helices packed along side of an antiparallel beta-sheet.

On the evolution of the single-subunit RNA polymerases

Many eukaryotic nuclear genomes as well as mitochondrial plasmids contain genes displaying evident sequence similarity to those encoding the single-subunit RNA polymerase (ssRNAP) of bacteriophage T7 and its relatives. We have collected and aligned these ssRNAP sequences and have constructed unrooted phylogenetic trees that demonstrate the separation of ssRNAPs into three well-defined and nonoverlapping clusters (phage-encoded, nucleus-encoded, and plasmid-encoded). Our analyses indicate that these three subfamiles of T7-like RNAPs shared a common ancestor; however, the order in which the groups diverged cannot be inferred from available data. On the basis of structural similarities and mutational data, we suggest that the ancestral ssRNAP gene may have arisen via duplication and divergence of a DNA polymerase or reverse transcriptase gene. Considering the current phylogenetic distribution of ssRNAP sequences, we further suggest that the origin of the ancestral ssRNAP gene closely paralleled in time the introduction of mitochondria into eukaryotic cells through a eubacterial endosymbiosis.

PolI-driven integrative expression vectors for yeast

A novel expression vector for yeast has been constructed from the regulatory elements present in the polI promoter and the enhancer/termination region (E/T) of rDNA. Under some conditions, this promoter/vector combination produces small RNAs such as the hammerhead RNA sequence at levels comparable to polII- and polIII-dependent systems. No stable transcription product can be demonstrated with this vector when the enhancer/termination sequence is less than 100 nucleotides downstream from the promoter. On the other hand, high expression of a stable, hammerhead RNA molecule can be obtained from this vector by inserting a 400-bp fragment containing the ADH1 transcription termination region upstream of the E/T. RNAs produced by this vector are polyadenylated and multiple copies of this plasmid can be stably integrated into the yeast chromosome.

An ancestral mitochondrial DNA resembling a eubacterial genome in miniature

Mitochondria, organelles specialized in energy conservation reactions in eukaryotic cells, have evolved from eubacteria-like endosymbionts whose closest known relatives are the rickettsial group of alpha-proteobacteria. Because characterized mitochondrial genomes vary markedly in structure, it has been impossible to infer from them the initial form of the proto-mitochondrial genome. This would require the identification of minimally derived mitochondrial DNAs that better reflect the ancestral state. Here we describe such a primitive mitochondrial genome, in the freshwater protozoon Reclinomonas americana. This protist displays ultrastructural characteristics that ally it with the retortamonads, a protozoan group that lacks mitochondria. R. americana mtDNA (69,034 base pairs) contains the largest collection of genes (97) so far identified in any mtDNA, including genes for 5S ribosomal RNA, the RNA component of RNase P, and at least 18 proteins not previously known to be encoded in mitochondria. Most surprising are four genes specifying a multisubunit, eubacterial-type RNA polymerase. Features of gene content together with eubacterial characteristics of genome organization and expression not found before in mitochondrial genomes indicate that R. americana mtDNA more closely resembles the ancestral proto-mitochondrial genome than any other mtDNA investigated to date.

Efficient hammerhead ribozyme and antisense RNA targeting in a slow ribosome Escherichia coli mutant

We have evaluated inhibition of the plasmid-born chloramphenicol acetyl transferase gene (CAT) by the hammerhead ribozyme and antisense RNA in Escherichia coli where the translation and transcription rates have been modified. Whereas neither antisense nor the hammerhead had an inhibitory effect on CAT activity in wild-type E. coli, both reduced the level of the messenger RNA and the activity of the CAT gene by almost 60% in a slow ribosome mutant. Streptomycin, which increases the speed of translation in this mutant strain, restored full CAT activity. The level of CAT activity expressed from a T7 RNA polymerase promoter was not affected by the presence of either antisense RNA or the hammerhead ribozyme. When the target gene was expressed from a chromosomal locus in wild-type E. coli, both antisense RNA and the hammerhead ribozyme showed some inhibitory activity, but the level of inhibition was significantly increased in the slow ribosome strain. This bacterial system offers a unique entry to the study of cellular factors which mediate the activity of ribozymes in vivo.

Effect of structural modifications on the activity of the leadzyme

The structure/function properties of functional groups in the leadzyme have been studied by assaying the activity of analog ribozymes generated by the systematic substitution of modified nucleotides in the internal loop region of the ribozyme. Guanosine analogs introduced at positions 4 and 7 occupied by guanosine in the wild-type molecule severely diminished cleavage. The substitution of deoxycytidine for cytidine at the cleavage site completely eliminated the activity of the leadzyme, as expected if the adjacent 2'-OH were the nucleophile in the cleavage reaction. On the other hand, substitution of an abasic nucleotide for adenosine at position 8 did not affect the activity of the ribozyme. An analysis of the activity of these analogs gives rise to the proposal of a triple-base pair motif implicating C1, G4, and G7.

Structural rules and conformational compensations in the tRNA L-form

The mitochondrial tRNAs (mtRNA) of five distinct, secondary structure types have been identified in the tRNA sequence compilation, and the three-dimensional modeling for representative sequences of these types has been carried out using a new criterion for the lengths of the helical domains and connector regions in a full-sized tRNA conformation. This criterion has been derived from the analysis of the known structures of cytosolic tRNAs and states that in the tRNA structure nucleotide 59 of the T-loop should stack onto Domain 1. To ensure this, Domain 1 must have 12 layers of stacked nucleotides, and in the case of a deletion of a base-pair in the T-stem, an additional 13th layer is required. Although a number of mitochondrial tRNAs harbored deficiencies in this criterion and, therefore, could not be modeled directly, this disability could be corrected and modeling accomplished by invoking structural compensations derived from one of two unusual aspects of these tRNAs. One class of these tRNAs contained an unpaired nucleotide in their anticodon stem, and their three-dimensional structure was successfully modeled when the unpaired nucleotide was intercalated into the helical domain of the stem. The second class contained more than the required number of nucleotides connecting the tRNA helical domains; the conformational flexibility of these nucleotides allowed them to take the place of the absent layers. The conformational compensation that we report rationalizes disparate features of these tRNAs and suggests that the stacking of nucleotide 59 on Domain 1 is an essential feature of the three-dimensional L-form of tRNA.

Predicting RNA structures: the model of the RNA element binding Rev meets the NMR structure

BACKGROUND

How accurate are the predictions of RNA three-dimensional structures? Assessing this accuracy requires the detailed comparison of the prediction with the experimentally determined structure. Previously, sequence variation in RNA aptamers that bind the Rev protein was used to infer a three-dimensional model of the Rev-binding element (RBE) RNA. Although much of this model has been substantiated by subsequent experimental data, its validity remains to be determined by confronting it with the structure determined by NMR spectroscopy.

RESULTS

A series of different criteria such as geometric parameters (root mean square deviation, interproton distances, torsions and puckering), helicoidal parameters (base pairing and base stacking) and stability considerations (conformational energies) have been evaluated to identify common and distinguishing structural characteristics of the model and the NMR structure.

CONCLUSIONS

The detailed comparison of the two structures reveals striking structural similarities at both the global and local level that validate the RNA modeling approach that we have used. Analysis of the structural differences and the precision of the model suggest that the limitations of the method are related to the amount of structural information available for modeling.

Stimulation of mitotic recombination upon transcription from the yeast GAL1 promoter but not from other RNA polymerase I, II and III promoters

Homologous recombination in Saccharomyces cerevisiae and other organisms can be stimulated by transcription. Consistent with this, we find that recombination of a chromosomal ade1 allele with a plasmid-borne ADE1 ORF under the control of the GAL1 promoter increased from 6.1x10(-6) to 1.7x10(-4) when transcription of the plasmid locus was induced by growing the cells in the presence of galactose. Recombination could also be stimulated by over-expressing the Gal4 transcription factor in the presence of the GAL1-ADE1 plasmid, while culturing the cells in dextrose medium. However, when transcription of the same ORF was driven from the highly active promoters of the rDNA (RNA polymerase I), and ADH1 (RNA polymerase II) genes, only background levels of recombination (5-10x10(-6)) were observed, irrespective of the carbon source. Recombination was found to involve integration of the whole plasmid and to depend on RAD51, RAD52 and RAD54. The results indicate that increased accessibility of transcriptionally active chromatin is not sufficient to cause increased rates of this kind of reciprocal exchange.

Cell cycle arrest promotes trans-hammerhead ribozyme action in yeast

A hammerhead ribozyme designed to cleave the yeast ADE1 mRNA has been expressed in yeast under the control of a galactose-inducible promoter. RNA prepared from the galactose-induced yeast cultures possesses an activity that cleaves ADE1 mRNA in vitro. However, in spite of high expression levels of the ribozyme, no cleavage activity could be demonstrated in vivo. On the other hand, when the yeast cells expressing hammerhead RNA were treated with the alpha-factor mating pheromone, the level of ADE1 mRNA was reduced by 50%. Similar reductions were observed when this strain was cultured in the presence of lithium acetate or in nitrogen-free medium. Moreover, control experiments in which disabled hammerhead genes were expressed showed no such reductions. Extension of the length of the flanking recognition arms of the ribozyme from a total of 10 to 16 or 24 nucleotides diminished the inhibitory effect of the ribozyme. These data suggest that ribozymes are able to cleave a trans-RNA target in yeast.

The puzzling origin of the genetic code

Recent results add to the mystery of the origin of the genetic code. In spite of early doubts, RNA can discriminate between hydrophobic amino acids under certain contexts. Moreover, codon reassignment, which has taken place in several organisms and mitochondria, is not a random process. Finally, phylogenies of some aminoacyl-tRNA synthetases suggest that the entire code was not completely assigned at the time of the divergence of bacteria from nucleated cells.

Sequences homologous to yeast mitochondrial and bacteriophage T3 and T7 RNA polymerases are widespread throughout the eukaryotic lineage

Although mitochondria and chloroplasts are considered to be descendants of eubacteria-like endo- symbionts, the mitochondrial RNA polymerase of yeast is a nucleus-encoded, single-subunit enzyme homologous to bacteriophage T3 and T7 RNA polymerases, rather than a multi-component, eubacterial-type alpha 2 beta beta' enzyme, as encoded in chloroplast DNA. To broaden our knowledge of the mitochondrial transcriptional apparatus, we have used a polymerase chain reaction (PCR) approach designed to amplify an internal portion of phage T3/T7-like RNA polymerase genes. Using this strategy, we have recovered sequences homologous to yeast mitochondrial and phage T3/T7 RNA polymerases from a phylogenetically broad range of multicellular and unicellular eukaryotes. These organisms display diverse patterns of mitochondrial genome organization and expression, and include species that separated from the main eukaryotic line early in the evolution of this lineage. In certain cases, we can deduce that PCR-amplified sequences, some of which contain small introns, are localized in nuclear DNA. We infer that the T3/T7-like RNA polymerase sequences reported here are likely derived from genes encoding the mitochondrial RNA polymerase in the organisms in which they occur, suggesting a phage T3/T7-like RNA polymerase was recruited to act in transcription in the mitochondrion at an early stage in the evolution of this organelle.

A docking and modelling strategy for peptide-RNA complexes: applications to BIV Tat-TAR and HIV Rev-RBE

BACKGROUND

In spite of the great interest in the interaction between RNAs and proteins, no general protocol for modelling these complexes is presently available. This methodological vacuum is particularly acute because the structure of few such complexes is known.

RESULTS

A general strategy for docking and modelling RNA-protein complexes has been developed. The docking procedure involves minimizing electrostatic and van der Waals' interaction energies of conformationally rigid structures during docking. After docking, libraries of amino acid sidechain conformations are searched to obtain the best interactions between the peptide and the RNA. Using this method, we have reproduced the structure of a bovine immunodeficiency virus (BIV) Tat peptide bound to BIV TAR RNA and have developed a model for the structure of the arginine-rich HIV-1 Rev peptide (Rev34-50) interacting with the Rev-binding element (RBE).

CONCLUSIONS

The resulting model of the Rev34-50-RBE complex predicts that although no single arginine sidechain is responsible for complex formation, residues Arg2, Arg5 and Arg11 are more important for binding than the other arginine residues in the peptide. One model is supported by binding measurements performed on wild-type and mutant RBE molecules with the peptide.

Using SAR and QSAR analysis to model the activity and structure of the quinolone-DNA complex

A set of 78 quinolone derivatives were used in a structure-activity study to identify structural features correlating with antibacterial activity. Distinct combinations of functional properties were identified for Gram-negative and Gram-positive bacteria. 3-D Quantitative structure-activity relationship (QSAR) studies identified specific hydrophobic, topologic and electronic properties of the molecules for both in vitro and in vivo activities. From these results, a three-dimensional model of a DNA-quinolone complex was built using molecular modeling techniques. It was based on the intercalation of quinolone into the double helix of DNA. We conclude that the intercalation model is consistent with most available data on the structure of the quinolone complex. This predicted structure is stabilized by the binding of magnesium ion with the sp2 oxygens present in quinolone, a phosphate and a purine base of the DNA. Substituents R1 and R7 are predicted to make hydrophobic interactions in the major and minor groove of DNA, respectively. R7 could also form hydrogen bonds with amino groups of guanines and the aspartic acid residue at position 87 in DNA gyrase subunit A.

A correlation between N2-dimethylguanosine presence and alternate tRNA conformers

Even though the evolutionary conservation of the cloverleaf model is strongly suggestive of powerful constraints on the secondary structure of functional tRNAs, some mitochondrial tRNAs cannot be folded into this form. From the optimal base pairing pattern of these recalcitrant tRNAs, structural correlations between the length of the anticodon stem and the lengths of connector regions between the two helical domains, formed by the coaxial stacking of the anticodon and D-stems and the acceptor and T-stems, have been derived and used to scan the tRNA and tRNA gene database. We show here that some cytosolic tRNA gene sequences that are compatible with the cloverleaf model can also be folded into patterns proposed for the unusual mitochondrial tRNAs. Furthermore, the ability to be folded into these atypical structures correlates in the mature RNA sequences with the presence of dimethylguanosine, whose role may be to prevent the unusual mitochondrial tRNA pattern folding.

An oligodeoxyribonucleotide that supports catalytic activity in the hammerhead ribozyme domain

A study of the activity of deoxyribonucleotide-substituted analogs of the hammerhead domain of RNA catalysis has led to the design of a 14mer oligomer composed entirely of deoxyribonucleotides that promotes the cleavage of an RNA substrate. Characterization of this reaction with sequence variants and mixed DNA/RNA oligomers shows that, although the all-deoxyribonucleotide oligomer is less efficient in catalysis, the DNA/substrate complex shares many of the properties of the all-RNA hammerhead domain such as multiple turnover kinetics and dependence on Mg2+ concentration. On the other hand, the values of kinetic parameters distinguish the DNA oligomer from the all-RNA oligomer. In addition, an analog of the oligomer having a single ribonucleotide in a strongly conserved position of the hammerhead domain is associated with more efficient catalysis than the all-RNA oligomer.

Structural and thermodynamic properties of DNA uncover different evolutionary histories

We propose an index of DNA homogeneity (IDH) based on a binary distribution model that quantifies structural and thermodynamic aggregates present in DNA primary structures. Extensive analysis of sequence databases with the IDH uncovers significant constraints on DNA sequence other than those derived from codon usage or protein function. This index clearly distinguishes between organisms of different evolutive origins and places them in disjoint domains of DNA sequence space.

A hammerhead ribozyme inhibits ADE1 gene expression in yeast

To study factors that affect in vivo ribozyme (Rz) activity, a model system has been devised in Saccharomyces cerevisiae based on the inhibition of ADE1 gene expression. This gene was chosen because Rz action can be evaluated visually by the Red phenotype produced when the activity of the gene product is inhibited. Different plasmid constructs allowed the expression of the Rz either in cis or in trans with respect to ADE1. Rz-related inhibition of ADE1 expression was correlated with a Red phenotype and a diminution of ADE1 mRNA levels only when the Rz gene was linked 5' to ADE1. The presence of the expected 3' cleavage fragment was demonstrated using a technique combining RNA ligation and PCR. This yeast system and detection technique are suited to the investigation of general factors affecting Rz-catalyzed inhibition of gene expression under in vivo conditions.

A three-dimensional model of the Rev-binding element of HIV-1 derived from analyses of aptamers

Coordinated variations in the sequence of the Rev-binding element of HIV-1, identified by in vitro genetic selections, have been used as distance and conformational constraints for molecular modelling. Three-dimensional models of the wild-type Rev-binding element and several, evolved RNA ligands (aptamers) have been constructed. These models demonstrate that non-Watson-Crick pairings open the major groove allowing access of an alpha-helical peptide from Rev, and explain why some selected RNA sequences can bind Rev more tightly than others.

Fitting the structurally diverse animal mitochondrial tRNAs(Ser) to common three-dimensional constraints

We propose three-dimensional models for animal mitochondrial (amt) tRNAs lacking the D-domain based on consideration of universal constraints on tRNA to maintain functionality. The available tRNA sequences are classified into two groups, and distinct models are proposed for both classes derived from common structural features. The distance between the anticodon and the acceptor stem is comparable in the models and corresponds to that observed in conventional tRNAs. This fact averts the problem of how a shorter mitochondrial tRNA could function within the context of a protein synthesis machinery suited to full-sized tRNAs. In the models, the angle which defines the relationship between the helical domains composed of the acceptor/T-stem and the anticodon/D-stem is greater than in conventional tRNAs. These structures resemble more a "boomerang" than an "L". However, even in the boomerang model, the inner surface of tRNA would be sufficiently uncluttered to avoid steric clashes when two tRNA molecules cohabit the ribosome.

Catalytic activity of an RNA domain derived from the U6-U4 RNA complex

U6 RNA contains two regions that are essential for proper splicing of nuclear precursor messenger RNA (pre-mRNA). A comparison of putative secondary structures of the U6-U4 RNA complexes from different phyla revealed a conserved domain that is similar to the catalytic hammerhead RNA motif. Although no catalytic activity was detected in the mammalian U6-U4 RNA complexes, two nucleotide changes in U6 RNA and one in U4 RNA conferred cleavage activity to the complex. Furthermore, the highly conserved domain of the wild-type complex, without the accompanying flanking regions, cleaved an RNA substrate and exhibited other characteristics of the hammerhead ribozyme. The possible involvement of this structure in pre-mRNA splicing is also discussed.

Reproducing the three-dimensional structure of a tRNA molecule from structural constraints

The three-dimensional structure of yeast tRNA(Phe) was reproduced at atomic resolution with the automated RNA modeling program MC-SYM, which is based on a constraint-satisfaction algorithm. Structural constraints used in the modeling were derived from the secondary structure, four tertiary base pairs, and other information available prior to the determination of the x-ray crystal structure of the tRNA. The program generated 26 solutions (models), all of which had the familiar "L" form of tRNA and root-mean-square deviations from the crystal structure in the range of 3.1-3.8 A. The interaction between uridine-8 and adenosine-14 was crucial in the modeling procedure, since only this among the tertiary pairs is necessary and sufficient to reproduce the L form of tRNA. Other tertiary interactions were critical in reducing the number of solutions proposed by the program.

Modeling the three-dimensional structure of RNA using discrete nucleotide conformational sets

The flexibility about seven torsion angles in nucleotides constitutes a severe obstacle to computer modeling of RNA. The computational feasibility of RNA conformational searches can be enhanced by assigning to each nucleotide a set of discrete conformations. In this work, four types of discrete conformational sets for the atomic representation of nucleotide structures were defined and evaluated. These sets, comprising between 10 and 30 conformations, were tested for their ability to reproduce known RNA structures and to generate structures responding to new specifications. Conformational searches were performed with the MC-SYM program, which allows for the generation of all structures satisfying a predetermined set of three-dimensional constraints in a given discrete space. Results with known hairpin loop structures show that root-mean-square deviations of about 1.5 A for backbone atoms and about 2.0 A for all atoms between the modeled and X-ray crystal structures can be expected. The conformational set that gives the most faithful representation of test structures is based on the classification of nucleotide conformations derived from a structural database. Representative conformations are selected from each class that adequately sample variations in backbone direction, sugar pucker and base orientation. With this conformational set, most of the important features of test hairpin structures are reproduced with fidelity, indicating that biologically useful models can be constructed from the combination of discrete nucleotide conformations and an algorithm that rapidly and systematically scans the pre-defined conformational space.

Modeling the three-dimensional structure of RNA

The limited number of RNA structures determined by X-ray crystallography and NMR spectroscopy compels the use of experimental and theoretical methods that are less precise to obtain information on RNA conformation. RNA flexibility, a consequence of rotational freedom about seven intra- and internucleotide bonds, is unfortunately of such magnitude that these alternate techniques fall short of providing sufficient information to build robust tertiary structures. Various RNA modeling methods, described herein, permit the organization of this structural data to the form of three-dimensional structures. Interactive computer graphics techniques, for example, have generated several useful models. Also, conventional computer algorithms involving the minimization of empirical energy functions, previously limited to small molecules, are giving way to methods able to handle much larger molecules. Modified distance geometry and molecular mechanics algorithms, using simplified "pseudoatom" representations, can generate structures consistent with input data. A constraint satisfaction algorithm combined with discrete representations of nucleotide conformations systematically explores poorly defined regions of a molecule yielding all-atom representations, but requires enough structural constraints to avoid a computational explosion.