PubMed Health⌕ Search

Biomedical subjects

A B Jacobson

Publications and source records attributed to A B Jacobson.

At least 19 recordsLinked to original sources

Structural plasticity in RNA and its role in the regulation of protein translation in coliphage Q beta.

We have analyzed both conformational and functional changes caused by two large cis-acting deletions (delta 159 and delta 549) located within the read-through domain, a 850 nucleotide hairpin, in coliphage Q beta genomic RNA. Studies in vivo show that co-translational regulation of the viral coat and replicase genes has been uncoupled in viral genomes carrying deletion delta 159. Translational regulation is restored in deletion delta 549, a naturally evolved pseudorevertant. Structural analysis by computer modeling shows that structural features within the read-through domain of delta 159 RNA are less well determined than they are in the read-through domain of wild-type RNA, whereas predicted structure in the read-through domain of evolved pseudorevertant delta 549 is unusually well determined. Structural analysis by electron microscopy of the genomic RNAs shows that several long range helices at the base of the read-through domain, that suppress translational initiation of the viral replicase gene in the wild-type genome, have been destabilized in delta 159 RNA. In addition, the structure of local hairpins within the read-through region is more variable in delta 159 RNA than in wild-type RNA. Stable RNA secondary structure is restored in the read-through domain of delta 549 RNA. Our analyses suggest that structure throughout the read-through domain affects the regulation of viral replicase expression by altering the likelihood that long-range interactions at the base of the domain will form. We discuss possible kinetic and equilibrium models that can explain this effect, and argue that observed changes in structural plasticity within the read-through domain of the mutant genomes are key in understanding the process. During the course of these studies, we became aware of the importance of the information contained in the energy dot plot produced by the RNA secondary structure prediction program mfold. As a result, we have improved the graphical representation of this information through the use of color annotation in the predicted optimal folding. The method is presented here for the first time.

Allolevivirus↗

Using reliability information to annotate RNA secondary structures.

A number of heuristic descriptors have been developed previously in conjunction with the mfold package that describe the propensity of individual bases to participate in base pairs and whether or not a predicted helix is "well-determined." They were developed for the "energy dot plot" output of mfold. Two descriptors, P-num and H-num, are used to measure the level of promiscuity in the association of any given nucleotide or helix with alternative complementary pairs. The third descriptor, S-num, measures the propensity of bases to be single-stranded. In the current work, we describe a series of programs that were developed in order to annotate individual structures with "well-definedness" information. We use color annotation to present the information. The programs can annotate PostScript files that are created by the mfold package or the PostScript secondary structure plots produced by the Weiser and Noller program XRNA (Weiser B, Noller HF, 1995, XRNA: Auto-interactive program for modeling RNA, The Center for Molecular Biology of RNA, Santa Cruz, California: University of California; Internet: ftp://fangio.ucsc.edu/pub/XRNA). In addition, these programs can annotate ss files that serve as input to XRNA. The annotation package can also handle structure comparison with a reference structure. This feature can be used to compare predicted structure with a phylogenetically deduced model, to compare two different predicted foldings, and to identify conformational changes that are predicted between wild-type and mutant RNAs. We provide several examples of application. Predicted structures of two RNase P RNAs were colored with P-num information and further annotated with comparative information. The comparative model of a 16S rRNA was annotated with P-num information from mfold and with base pair probabilities obtained from the Vienna RNA folding package. Further annotation adds comparisons with the optimal foldings obtained from mfold and the Vienna package, respectively. The results of all of these analyses are discussed in the context of the reliability of structure prediction.

Computer Simulation↗

cis-acting elements within an RNA coliphage genome: fold as you please, but fold you must!!

Using an in vivo complementation system, we conducted a mutational analysis of the bacteriophage Q beta readthrough cistron. In the Q beta cDNA-containing plasmid, pQ beta m100, we constructed six defined Q beta deletion cDNA genomes, each missing between 86 and 447 nucleotides from within the readthrough cistron. These deletion plasmids were introduced into host cells that are constitutively supplied with Q beta readthrough protein from the plasmid pQ beta RT. Under these conditions, all six deletion genomes spontaneously generated phage particles, each exhibiting a characteristic plaque phenotype and virus forming potential. Isolated readthrough-defective phage particles were subsequently used to infect host cells that carried helper readthrough protein. Passaged viruses yielded both larger plaques and higher titers, compared with those of the parent phages. Sequence analysis revealed that the genomes of the passaged viruses had deleted additional regions of readthrough RNA sequence. We discuss the possibilities that (1) the disruption of a well-defined structural domain in Q beta RNA was selectively disadvantageous to phage infection, and that (2) the evolved viral populations were selected by virtue of their ability to restore critical integrity of short and/or long-range nucleotide interactions within this region of Q beta RNA.

Base Sequence↗

"Well-determined" regions in RNA secondary structure prediction: analysis of small subunit ribosomal RNA.

Recent structural analyses of genomic RNAs from RNA coliphages suggest that both well-determined base paired helices and well-determined structural domains that are identified by "energy dot plot" analysis using the RNA folding package mfold, are likely to be predicted correctly. To test these observations with another group of large RNAs, we have analyzed 15 ribosomal RNAs. Published secondary structure models that were derived by comparative sequence analysis were used to evaluate the predicted structures. Both the optimal predicted fold and the predicted "energy dot plot" of each sequence were examined. Each prediction was obtained from a single computer run on an entire ribosomal RNA sequence. All predicted base pairs in optimal foldings were examined for agreement with proven base pairs in the comparative models. Our analyses show that the overall correspondence between the predicted and comparative models varied for different RNAs and ranges from a low of 27% to high of 70%, with a mean value of 49%. The correspondence improves to a mean value of 81% when the analysis is limited to well-determined helices. In addition to well-determined helices, large well-determined structural domains can be observed in "energy dot plots" of some 16S ribosomal RNAs. The predicted domains correspond closely with structural domains that are found by the comparative method in the same RNAs. Our analyses also show that measuring the agreement between predicted and comparative secondary structure models underestimates the reliability of structural prediction by mfold.

Base Sequence↗

A two-dimensional model at the nucleotide level for the central hairpin of coliphage Q beta RNA.

We describe the construction and testing of a structural model at the nucleotide level for conformation CH of the central hairpin of genomic RNA from coliphage Q beta. The model was developed with the computer program MFOLD using both optimal and suboptimal predictions. Structural information obtained by electron microscopic analysis of Kleinschmidt spreadings of Q beta RNA was used to guide the modeling. The model was tested in solution with three enzymatic probes: RNase T1, RNase T2, and RNase V1, as well as four chemical probes: dimethylsulfate, diethylpyrocarbonate, kethoxal and 1-cyclohexyl-3-(2-morpholinoethyl)carbodiimide metho-p-toluene sulfonate (CMCT). The structural analyses in solution are consistent with the predicted structural model. The model is also supported by comparative structural analysis with the related coliphage SP. The model provides a structural basis for published biochemical and genetic studies implicating large, long-range structural features in the co-regulation of viral coat and replicase expression. In addition, we show that the read-through region of the viral protein A1 forms a separate structural domain, and we suggest that it functions as a nucleation site that participates in the folding and refolding of the molecule during replication and translation. In addition to the central hairpin, we have analyzed the structure of the viral coat initiation region. Our studies show that the entire region consists of small local hairpins and that 26 nucleotides immediately surrounding the coat initiation codon are single-stranded.

Base Sequence↗

Structural analysis by energy dot plot of a large mRNA.

We have predicted the secondary structure of the entire 4217 nucleotide sequence of the genomic RNA of coliphage Q beta in one computer run using the computer program MFOLD that computes RNA structures within any prescribed increment of the computed minimum free energy. The results are presented in the form of an "energy dot plot" that shows both an optimal folding as well as the superposition of all base-pairs that can form in slightly suboptimal foldings. The plot reveals five large, well-determined, independent structural domains that cover approximately 50% of the viral genome. The predicted structural domains are consistent with and provide support for five large structural domains identified previously by quantitative electron microscopy in Q beta RNA. The dot plot also contains cluttered regions that indicate large numbers of alternative foldings within or between segments of an RNA molecule. These reflect the impossibility of accurate structure prediction and/or the biological reality of more than one folding. Weaker, long range structures, that are observed by electron microscopy in two alternate competing conformations, are located in the regions of the Q beta sequence that correspond to cluttered regions of the dot plot. The potential biological significance of these secondary structures is discussed.

Coliphages↗

Secondary structure of coliphage Q beta RNA. Analysis by electron microscopy.

The secondary structure of genomic RNA from the coliphage Q beta has been examined by electron microscopy in the presence of varying concentrations of spermidine using the Kleinschmidt spreading technique. The size and position of structural features that cover 70% of the viral genome have been mapped. The structural features that are visualized by electron microscopy in Q beta RNA are large. They range in size from 170 to 1600 nucleotides. A loop containing approximately 450 nucleotides is located at the 5' end of the RNA. It includes the initiation region for the viral maturation protein. A large hairpin containing approximately 1600 nucleotides is located in the center of the molecule. It is multibranched and includes most of the viral coat gene, the readthrough region of the A1 gene, and approximately one third of the viral replicase gene. Within the central hairpin, the initiation region for the viral replicase gene pairs with a region within the distal third of the viral coat gene. This structure may participate in the regulation of translational initiation of the viral replicase gene. Two structural variants of the central hairpin were observed. One of them brings the internal S and M viral replicase binding regions into juxtaposition. These observations suggest that the central hairpin may also participate in the regulation of translation of the viral coat gene. The secondary structures that are observed in Q beta RNA differ significantly from structures that we described previously in the genomic RNA of coliphage MS2 but are similar to structures we observed by electron microscopy in the related group B coliphage SP.

Coliphages↗

A computer method for finding common base paired helices in aligned sequences: application to the analysis of random sequences.

We describe a new computer program that identifies conserved secondary structures in aligned nucleotide sequences of related single-stranded RNAs. The program employs a series of hash tables to identify and sort common base paired helices that are located in identical positions in more than one sequence. The program gives information on the total number of base paired helices that are conserved between related sequences and provides detailed information about common helices that have a minimum of one or more compensating base changes. The program is useful in the analysis of large biological sequences. We have used it to examine the number and type of complementary segments (potential base paired helices) that can be found in common among related random sequences similar in base composition to 16S rRNA from Escherichia coli. Two types of random sequences were analyzed. One set consisted of sequences that were independent but they had the same mononucleotide composition as the 16S rRNA. The second set contained sequences that were 80% similar to one another. Different results were obtained in the analysis of these two types of random sequences. When 5 sequences that were 80% similar to one another were analyzed, significant numbers of potential helices with two or more independent base changes were observed. When 5 independent sequences were analyzed, no potential helices were found in common. The results of the analyses with random sequences were compared with the number and type of helices found in the phylogenetic model of the secondary structure of 16S ribosomal RNA. Many more helices are conserved among the ribosomal sequences than are found in common among similar random sequences. In addition, conserved helices in the 16S rRNAs are, on the average, longer than the complementary segments that are found in comparable random sequences. The significance of these results and their application in the analysis of long non-ribosomal nucleotide sequences is discussed.

Base Composition↗

Analysis of the complete nucleotide sequence of the group IV RNA coliphage SP.

We report the nucleotide sequence of the Group IV RNA bacteriophage SP. The entire sequence is 4276 nucleotides long. Four cistrons have been identified by comparison with the related Group III phage Q beta. The maturation protein contains 449 amino acids, the coat protein contains 131 amino acids, the read-through protein contains 330 amino acids and the replicase beta-subunit contains 575 amino acids. SP is 59 nucleotides longer than Q beta. We have analyzed both sequence and structural conservation between SP and Q beta and shown that the sequences for the coat and central region of the replicase are strongly conserved between the two genomes. We also show that the S and M replicase binding sites of Q beta are strongly conserved in SP. Interestingly, the base composition of SP and Q beta differ significantly from one another, and most of the differences can be accounted for by a strong preponderance of U in the third position of each codon of Q beta relative to SP. We also compare conserved hairpins associated with potential coat protein and replicase binding sites.

Base Sequence↗

The complete nucleotide sequence of the group II RNA coliphage GA.

The complete nucleotide sequence of the RNA coliphage GA, a group II phage, is presented. The entire genome comprises 3466 bases. Three large open reading frames were identified, which correspond to the maturation protein gene (390 amino acids), the coat protein gene (129 amino acids) and the replicase beta-subunit protein gene (531 amino acids). In addition, untranslated regions occur at the 5' (135 bases) and 3' (122 bases) ends of the molecule. Two intercistronic untranslated regions occur between the cistrons for the maturation and coat proteins, and between the coat and beta-subunit proteins. We have compared the nucleotide sequence of GA RNA with the published sequence of MS2 RNA, and show that they are related. The comparative structures of two important regulatory regions are presented; the coat protein binding site which is involved in translational repression of the replicase beta-subunit protein gene, and a hairpin in a region proximal to the lysis protein gene.

Base Sequence↗

Effect of spermidine on the conformation of bacteriophage MS2 RNA. Electron microscopy and computer modeling.

The structure of single-stranded RNA from the bacteriophage MS2 has been examined by electron microscopy in the presence of the polyamine spermidine. The molecules are found in two alternate conformations. The first of these can be characterized as a cruciform structure composed of three large loops approximately 500 to 700 nucleotides in size. The interior of the molecule has extensive base-paired regions which connect distant regions of the molecule; the farthest being 2500 nucleotides apart. In the second conformation, the molecules appear rod-like. Two of the large loops disappear, and these regions form, instead, extensive long-range helices. Computer modeling has been employed to explore the base-pairing potential of the sequence of bacteriophage MS2 RNA. Double-stranded regions identified by electron microscopy are shown to occur in local G + C-rich stretches of the RNA. Detailed models have been calculated for two regions of long-range contact. One of these includes the ribosome-binding site for the viral coat protein gene. The results are discussed in the context of the known role of RNA structure in the regulation of viral gene expression.

Bacteriophages↗

Some simple computational methods to improve the folding of large RNAs.

Computational methods are described which increase the efficiency of the RNA folding algorithm described by Zuker and Stiegler. Bit addressing has been used to reduce the memory requirements from 2NxN to NxN/2. The order in which the nucleotide sequence is examined internally has been altered, and some additional short arrays which carry temporary information have been introduced. These changes optimize the management of the large data arrays generated by the algorithm. The methods were developed for use with a UNIVAC 1100/82 computer. They are, however, easily adaptable to other computers; especially those with virtual memory capabilities. The analysis of sequences up to 1000 nucleotides long are relatively routine, and larger searches are also feasible. Some limitations and applications of the algorithm are also discussed.

Base Sequence↗

Small changes in free energy assignments for unpaired bases do not affect predicted secondary structures in single stranded RNA.

We present extensive calculations of the secondary structure of mRNA which point to its insensitivity to small changes in the free energy assignments of single stranded regions. Truncating the free energies of hairpin loops, bulges, internal loops and multibranched junctions to two significant digits yields structures nearly identical to those generated using three digit values. The results show that one can safely use truncated values in RNA folding calculations. The implementation of these results enabled us to carry out secondary structure calculations on 2600 nucleotides in a single computer run.

Base Composition↗

Secondary structure model for the complete simian virus 50 late precursor mRNA.

Structures for all sequences containing less than 1790 nucleotides in the 2600 nucleotide late region of the SV40 virus have been computed and saved on magnetic tape. Previously the longest sequence whose secondary structure was calculated in a single computer run contained 950 nucleotides. In the past, analysis of long molecules required numerous repeated, partially overlapping computations on much shorter segments. The structure obtained for the late half of the SV40 is Y-shaped with two unequal arms. It has 52 short hairpins. Two long range interactions between nucleotides near 650 and 1350 and between 1450 and 2450 appear to play an important role. The first is within the 16S intron; the second is in the 3' exon. The 5' and 3' ends of the molecule are close to each other and are found in the major elongated stem in the vicinity of the fork.

Base Sequence↗

Fast algorithm for predicting the secondary structure of single-stranded RNA.

A computer method is presented for finding the most stable secondary structures in long single-stranded RNAs. It is 1-2 orders of magnitude faster than existing codes. The time required for its application increases as N3 for a chain N nucleotides long. As many as 1000 nucleotides can be searched in a single run. The approach is systematic and builds an optimal structure in a straightforward inductive procedure based on an exact mathematical algorithm. Two simple half-matrices are constructed and the best folded form is read directly from the second matrix by a simple back-tracking procedure. The program utilizes published values for base-pairing energies to compute one structure with the lowest free energy.

Base Sequence↗

Studies on secondary structure of single-stranded RNA from bacteriophage MS2 by electron microscopy.

A method allowing the demonstration and study by electron microscopy of secondary structure of viral RNA has been developed. Single-stranded RNA from the bacteriophage MS2 has been analyzed in the electron microscope in the presence of various concentrations of MgCl2. Depending on the salt concentration, the molecules display one to three large open loops which range in size from 10 to 20% of the total RNA length, and smaller closed loops which are approximately 3-5% of the total RNA length. Within one spreading, the conformation of the molecules is variable. However, the average complexity of the molecules increases with increasing salt, and individual loops which are infrequent at low salt increase in frequency with increasing salt. By analyzing the manner in which the individual loop appeared, it was possible to show that all molecules could be described by one basic pattern of secondary structure formation.

Coliphages↗