Acrylamide-gel electrophorograms by mechanical fractionation: radioactive adenovirus proteins.

A mechanical fractionator was developed to produce electrophorograms by extrusion of polyacrylamide gels through a narrow orifice in a continuous, sequential stream. The system permits separation of uniform fractions free of zone distortion. An electrophorogram of radioactive type-2 adenovirus proteins so fractionated gave a pattern in excellent agreement with the pattern obtained by laborious manual sectioning and in agreement with the pattern obtained on a replicate gel stained with Coomassie brilliant blue R250. The adenovirus particle yielded about ten resolvable protein components in unequal amounts. Like picornaviruses, these icosahedral animal viruses have multiple protein components in the viral coat.

Properties of a normal mouse cell DNA sequence (sarc) homologous to the src sequence of Moloney sarcoma virus.

A 15.0-kilobase (kb) Eco RI DNA fragment from normal mouse Balb/c genomic DNA that contains sequences (sarc) homologous to the acquired cell sequences (src) of Moloney sarcoma virus (MSV) has been cloned in phage lambda. The sarc region (1.2 to 1.3 kb) of the 15.0-kb cell fragment is indistinguishable from the src region of two isolates of MSV as judged by heteroduplex and restriction endonuclease analyses. The cellular sequences flanking sarc show no homology to other MSV sequences. Whereas cloned subgenomic portions of MSV that contain src transformed NIH-3T3 cells in vitro, the cloned sarc fragment is inactive.

Data mining of imperfect double-stranded RNA in 3' untranslated regions of eukaryotic mRNAs.

Recent developments in the study of RNA silencing indicate that double-stranded RNA (dsRNA) can be used in eukaryotes to block expression of a corresponding cellular gene. There is also a large class of small non-coding RNAs having potential to form a distinct, stable stem-loop in numbers of eukaryotic genomes. We had reported that a large imperfect dsRNA structure with hundreds of base-pairs (bp) in the 3' untranslated region (3' UTR) of cytotoxic ribonuclease was correlated with the translation suppression. In this study, we search for such dsRNAs in a 3' UTR database. The occurrence rate of large dsRNA in 3' UTRs ranges from 0.01% in plant to 0.30% in vertebrate mRNAs. However, small imperfect dsRNAs of ~ 30 bp are much more prevalent than large ones. The small dsRNAs are statistically very significant and uniquely well-ordered. Most of them have the conserved structural features of pre-miRNAs. Our data mining of the dsRNAs in the 3' UTR database can be used to explore RNA-based regulation of gene expression.

Prediction of common secondary structures of RNAs: a genetic algorithm approach.

In this study we apply a genetic algorithm to a set of RNA sequences to find common RNA secondary structures. Our method is a three-step procedure. At the first stage of the procedure for each sequence, a genetic algorithm is used to optimize the structures in a population to a certain degree of stability. In this step, the free energy of a structure is the fitness criterion for the algorithm. Next, for each structure, we define a measure of structural conservation with respect to those in other sequences. We use this measure in a genetic algorithm to improve the structural similarity among sequences for the structures in the population of a sequence. Finally, we select those structures satisfying certain conditions of structural stability and similarity as predicted common structures for a set of RNA sequences. We have obtained satisfactory results from a set of tRNA, 5S rRNA, rev response elements (RRE) of HIV-1 and RRE of HIV-2/SIV, respectively.

A gender-specific mRNA encoding a cytotoxic ribonuclease contains a 3' UTR of unusual length and structure.

A cDNA (2855 nt) encoding a putative cytotoxic ribonuclease (rapLR1) related to the antitumor protein onconase was cloned from a library derived from the liver of gravid female amphibian Rana pipiens. The cDNA was mainly comprised (83%) of 3' untranslated region (UTR). Secondary structure analysis predicted two unusual folding regions (UFRs) in the RNA 3' UTR. Two of these regions (711-1442 and 1877-2130 nt) contained remarkable, stalk-like, stem-loop structures greater than 38 and 12 standard deviations more stable than by chance, respectively. Secondary structure modeling demonstrated similar structures in the 3' UTRs of other species at low frequencies (0.01-0.3%). The size of the rapLR1 cDNA corresponded to the major hybridizing RNA cross-reactive with a genomic clone encoding onconase (3.6 kb). The transcript was found only in liver mRNA from female frogs. In contrast, immunoreactive onconase protein was detected only in oocytes. Deletion of the 3' UTR facilitated the in vitro translation of the rapLR1 cDNA. Taken together these results suggest that these unusual UFRs may affect mRNA metabolism and/or translation.

Sequence signals in eukaryotic upstream regions.

In eukaryotes, as in prokaryotes, evidence is accumulating showing that transcription factors recognize and bind to certain promoter elements. Sequence and chromatin structure perturbation specify the transcription initiation site and govern its efficiency. Two oligomers have been implicated in these processes: the TATAAATA and CCAAT. In the present work, all mammalian, non-mammalian vertebrate and invertebrate sequences accumulated in the database have been aligned by their mRNA start positions and scanned for recurrences of the 32 complementary triplets. The more significant signals are summarized here. In particular, TAT/ATA recurs very frequently further upstream, at -275, in addition to the 'classical' -40 position. Downstream their level is very low. The CAAT box complementary triplet components are not among the more striking signals. Closer examination of the -275 region indicates that to a large extent the signal is due to both the ATAT and the TATA quartets. Comparison of the frequencies of these quartets at -275 with the CAAT quartet at -80 suggests that the former signal is twice the strength of the latter. The oligomers' distribution charts support notions that several components are involved in recognition, making the regulatory regions more robust and less sensitive to mutations.

Enhancer elements share local homologous twist-angle variations with a helical periodicity.

Several experiments have shown that some enhancers can be exchanged between different genomes. The transferred enhancers were functional (cf. Levinson, B., Khoury, G., Vande Woude, G. and Gruss, P. (1982) Nature 295, 568-572). This argues that these exchanged fragments are recognized as enhancers and possess some common characteristics which other sequences lack. Extensive comparisons of enhancers yielded only very limited nucleotide sequence homology, which appears to be insufficient for enhancer recognition. We suggest that the enhancers located and sequenced to date have recurring, periodic homologous twist-angle (tg) patterns. This helical periodicity and the symmetric nature of the repeating twist-angle features present a recurring spatial geometry. It also offers a possible explanation of the fact that inverted enhancers are still functional. Regions of large twist, roll or main-chain torsion angle delta deviations from regular B-DNA may facilitate enhancer recognition especially when distant from promoter elements. Tissue specificity may be encoded in additional sequence or structural features.

DNAase I hypersensitive sites may be correlated with genomic regions of large structural variation.

Helical-twist, roll and torsion-angle variations calculated by the Calladine (1982)-Dickerson (1983) rules were scanned along several nucleotide sequences for which DNAase I cleavage data are available. It has been shown that for short synthetic oligomers DNAase I cuts preferentially at positions of high helical twist (Dickerson & Drew, 1981; Lomonossoff et al., 1981). Our calculations indicate that DNAase I sensitive and hypersensitive sites in chromatin are correlated with regions of successive, large, helical-twist angle variations from regular B-DNA. In many cases these regions exhibit large variations in base-pair roll and backbone torsion angles as well. It has been suggested that DNAase I cuts in the vicinity of cruciforms. However, it was recently demonstrated by Courey & Wang (1983) and Gellert et al. (1983) that such cruciform formation in a negatively supercoiled DNA is kinetically forbidden under physiological conditions. We thus propose that clustering of large twist-angle (and/or roll and backbone torsion angle) variations may be among the conformational features recognized by the enzyme. Specific cuts can then preferentially occur at base-pair steps with high helical twists.

Molecular substructure of a viral receptor-recognition protein. The gp17 tail-fiber of bacteriophage T7.

The bacteriophage T7 tail complex consists of a conical tail-tube surrounded by six kinked tail-fibers, which are oligomers of the viral protein gp17 (Mr 61,400). We have derived a molecular model for the tail-fiber by integrating secondary structure predictions with ultrastructural information obtained by correlation averaging of electron micrographs of negatively stained tail complexes. This model has been further refined by high-resolution scanning transmission electron microscopy of purified fibers, both negatively stained and unstained. Mass measurements made from the latter images establish that the fiber is a trimer of gp17. The proximal half-fiber is a uniform rod, about 2.0 nm in diameter and 16.4 nm long, which we infer to be a triple-stranded coiled-coil, containing three copies of an alpha-helical domain of about 117 residues, starting at Phe151. The distal half-fiber is 15.5 nm long, and is made up of four globules, 3.1 to 4.8 nm in diameter, in rigid linear array: it contains the carboxy-terminal halves (residues approximately 268 to 553) of the constituent gp17 chains, arranged with 3-fold symmetry around its long axis. The amino-terminal domains (residues 1 to 149) link the fiber to the tail-tube. We conclude that the three gp17 chains are quasi-equivalent in the proximal half-fiber, equivalent in the distal half-fiber, and non-equivalent in the kink region that separates the two half-fibers: such localized non-equivalence may represent a general mechanism for the formation of kinked joints in segmented homo-oligomeric proteins.

Localization and molecular modelling of the membrane-inserted domain of the ninth component of human complement and perforin.

Upon interaction with the membrane-bound C5b-8 complex, the ninth component of complement (C9) unfolds and inserts into the membrane of cells on which surface complement has been activated. Consequently C9 oligomerization occurs and transmembrane channels of varying sizes are formed. The domain of the unfolded protein interacting with the cell membrane has so far not been identified since, unlike many integral membrane proteins, the C9 sequence does not contain a continuous stretch of hydrophobic amino acids. We studied the interaction of C9 with the lipid bilayer using the membrane-restricted photoaffinity label 3-(trifluoromethyl)-3-(m[125I]iodophenyl)diazirine (125I-TID). C9 was assembled on liposomes and after photoactivation, several labeled and non-labeled peptides, obtained by chemical and enzymatic cleavage or the 125I-TID-labeled C9, were analyzed. The segment from 176 to 345 was identified as the region containing the membrane-interacting structure. By means of secondary structure predictions, we identified two amphipathic alpha-helices (292-308 and 313-333) separated by a turn (309-312). Based on these results, we constructed a molecular model for the membrane-spanning region of C9. By analogy, we also constructed a model for this domain in perforin/cytolysin, a pore-forming protein found in the cytoplasmic granules of cytotoxic T-lymphocytes.

Distinguishing between sequential and nonsequentially folded proteins: implications for folding and misfolding.

We describe here an algorithm for distinguishing sequential from nonsequentially folding proteins. Several experiments have recently suggested that most of the proteins that are synthesized in the eukaryotic cell may fold sequentially. This proposed folding mechanism in vivo is particularly advantageous to the organism. In the absence of chaperones, the probability that a sequentially folding protein will misfold is reduced significantly. The problem we address here is devising a procedure that would differentiate between the two types of folding patterns. Footprints of sequential folding may be found in structures where consecutive fragments of the chain interact with each other. In such cases, the folding complexity may be viewed as being lower. On the other hand, higher folding complexity suggests that at least a portion of the polypeptide backbone folds back upon itself to form three-dimensional (3D) interactions with noncontiguous portion(s) of the chain. Hence, we look at the mechanism of folding of the molecule via analysis of its complexity, that is, through the 3D interactions formed by contiguous segments on the polypeptide chain. To computationally splice the structure into consecutively interacting fragments, we either cut it into compact hydrophobic folding units or into a set of hypothetical, transient, highly populated, contiguous fragments ("building blocks" of the structure). In sequential folding, successive building blocks interact with each other from the amino to the carboxy terminus of the polypeptide chain. Consequently, the results of the parsing differentiate between sequentially vs. nonsequentially folded chains. The automated assessment of the folding complexity provides insight into both the likelihood of misfolding and the kinetic folding rate of the given protein. In terms of the funnel free energy landscape theory, a protein that truly follows the mechanism of sequential folding, in principle, encounters smoother free energy barriers. A simple sequentially folded protein should, therefore, be less error prone and fold faster than a protein with a complex folding pattern.

Difference in sizes of human compared to murine alpha-subunits of the glycoprotein hormones arises by four-codon gene deletion or insertion.

The sizes of the human and subhuman alpha-subunits of the glycoprotein hormones differ by four amino acids (hCG alpha, 92 amino acids; murine, equine, bovine, and ovine alpha, 96 amino acids). The shortening of the human alpha-subunit has been attributed to posttranslational proteolysis. We have recently determined the nucleotide sequences of the mRNAs encoding the precursors of the alpha-subunit of mouse TSH and rat gonadotropins using recombinant DNA techniques. In this report, we have compared these nucleotide sequences and their deduced amino acid sequences with those of the pre- alpha-subunit of hCG (hCG pre-alpha) and the gene encoding the human alpha-subunit. We show that the difference in size of four amino acids between the apoproteins of the murine and human alpha-subunits is a result of a deletion, or insertion, of four codons close to the site of the second intron in the human gene. In addition, the sequence of four codons absent in the mRNA encoding hCG pre-alpha is similar to a region in the 3' end of this intron. These findings indicate that the differences in size of these alpha-subunits originate at the gene level rather than at the posttranslational level at which proteolytic processing occurs.

Unusual folding regions and ribosome landing pad within hepatitis C virus and pestivirus RNAs.

A statistically significant folding region is identified in the 5' untranslated region (5'-UTR) of hepatitis C virus (HCV), bovine viral diarrhea virus and hog cholera virus. This unusual folding region (UFR) detected in HCV encompasses 199 nucleotides (nt) and coincides with the reported internal ribosome entry site or ribosome landing pad (RLP), as determined by the 5' and 3' deletions [Tsukiyama-Kohara et al., J. Virol. 66 (1992) 1476-1483]. The RNA structure predicted in the UFR of HCV consists of a large stem-loop and a pseudoknot. The proposed structural model is consistent with RNase sensitivity studies [Brown et al., Nucleic Acids Res. 20 (1992) 5041-5045]. Moreover, the structure is highly conserved among these divergent HCV and pestivirus RNAs. The covariation of paired bases in the helical regions offers support for the proposed structural models. The pseudoknot predicted in these UFR shares a similar structural feature to those proposed in the RLP of cardioviruses, aphthoviruses and hepatitis A virus. Based on the common structural motif, a putative base-pairing model between HCV RNA and 18S rRNA, as well as pestiviral RNAs and 18S rRNA are suggested. Intriguingly, the proposed base-pairing models in this study are comparable to those proposed in picornaviruses in terms of their folded shape and location of the predicted complementary sequences between viral RNAs and 18S rRNA. Taken together, we suggest that the common base-pairing model between the UFR detected in the 5'-UTR of pestivirus and HCV and 18S rRNA have a general function in the internal initiation of cap-independent translation.

Identification of unusual RNA folding patterns encoded by bacteriophage T4 gene 60.

A 50-nucleotide (nt) untranslated region (coding gap sequence) that interrupts the amino acid coding sequence in T4 gene 60, plus an additional 5 nt upstream and another 3 nt downstream from the gap sequence, shows unusual folding patterns according to RNA structure prediction. A predicted highly stable and significant hairpin structure in the 5' half of the gap sequence and a plausible tertiary structural element computed in the 3' part of the gap sequence seem significant by statistical tests on the wild-type (wt) sequence. This feature is absent in insertion, deletion and substitution variants of the gap sequence, in which template activities are markedly lower than that of the wt. The proposed feature is consistent with currently available data showing that the translational bypass of the coding gap is correlated with a stop codon involved in a stem-loop structure folded in the gap sequence. We suggest that the role of this segment in 'ribosomal bypass' of a portion of the mRNA sequence is a property of its special folded structure.

HIV-1 reverse transcriptase: structure predictions for the polymerase domain.

Reverse transcriptase (RT) plays an essential role in the life cycle of the human immunodeficiency viruses (HIV). A better understanding of this enzyme, and its two catalytic functions, the DNA polymerase and the RNase H, could lead to the development of new drugs that would specifically block HIV replication. The available genetic, sequence, biochemical, and immunological data on the reverse transcriptase of HIV-1 constrain the possible structure of the DNA polymerase domain. The purpose of this review is to correlate the data and to discuss, in light of that data, a model for the structure of the polymerase domain. In this model, the polymerase domain is approximately 50 to 60 A in diameter with a 20 A opening to accommodate the nucleic acid duplex. The most evolutionarily conserved region of RT (amino acids 20-190 of HIV-1 RT) is proposed to form the inner surface of the 20 A opening to which the nucleic acid hemiduplex is bound.

Short-strong hydrogen bonds and a low barrier transition state for the proton transfer reaction in RNase A catalysis: a quantum chemical study.

There is growing evidence that some enzymes catalyze reactions through the formation of short-strong hydrogen bonds as first suggested by Gerlt and Gassman. Support comes from several experimental and quantum chemical studies that include correlation energies on model systems. In the present study, the process of proton transfer between hydroxyl and imidazole groups, a model of the crucial step in the hydrolysis of RNA by the enzymes of the RNase A family, is investigated at the quantum mechanical level of density functional theory and perturbation theory at the MP2 level. The model focuses on the nature of the formation of a complex between the important residues of the protein and the hydroxyl group of the substrate. We have also investigated different configurations of the ground state that are important in the proton transfer reaction. The nature of bonding between the catalytic unit of the enzyme and the substrate in the model is investigated by Bader's atoms in molecule theory. The contributions of solvation and vibrational energies corresponding to the reactant, the transition state and the product configurations are also evaluated. Furthermore, the effect of protein environment is investigated by considering the catalytic unit surrounded by complete proteins--RNase A and Angiogenin. The results, in general, indicate the formation of a short-strong hydrogen bond and the formation of a low barrier transition state for the proton transfer model of the enzyme.

An interactive dot matrix system for locating potentially significant features in nucleic acid molecules.

An interactive computer system using a dot matrix approach has been developed and used to determine potentially significant features due to distortions in the B-DNA helix as a result of variations of purine and pyrimidine patterns. Sequences were compared using matrices which were generated using the Calladine-Dickerson rules (C.R. Calladine, J. Mol. Biol. 161, 343-352, 1982 and R.E. Dickerson, J. Mol. Biol. 166, 419-441, 1983). Having control over various parameters to enhance different aspects of the visual appearance of these matrices was helpful in discovering patterns that were not known a priori. Specifically, it was found that a pattern of alternating doublets of purines and pyrimidines appear to exist in regulatory regions. This event is shown to be beyond probabilistic expectation.

Ion-RNA interactions in the RNA pseudoknot of a ribosomal frameshifting site: molecular modeling studies.

The three-dimensional (3-D) structure of a RNA pseudoknot that causes the efficient ribosomal frameshifting in the gag-pro region of mouse mammary tumor virus (MMTV) has been determined recently by nuclear magnetic resonance (NMR) studies. But since the structure refinement in the studies did not use metal ions and waters, it is not clear how metal ions participate in the stabilization of the pseudoknot, and what kind of ion-RNA interactions dominate in the tertiary contacts for the RNA pseudoknotting. Based on the reported structure data of the pseudoknot VPK of MMTV, we gradually refined the structure by restrained molecular dynamics (MD) using NMR distance restraints. Restrained MD simulation of the RNA pseudoknot was performed with sodium ions and water molecules. Our results are in good agreement with known NMR data and delineate the importance of the metal ion coordination in the stability of the pseudoknot. In the non-coaxially stacking pseudoknot, stem 1 (S1), stem 2 (S2), and the intervening A14 involves unconventional stacking of base pairs coordinated by Na+ and/or bridging water molecules. A6 and G7 of loop L1 make a perfect base stacking in the major groove and are further stabilized by coordinated Na+ ions and water molecules. The first 4-nucleotide (nt) ACUC of loop L2 form a sharp turn and the following 4-nt AAAA cross the minor groove of S1 and are steadied by interactions with the nucleotides of S , bridging water molecules and coordinated Na+ ions. Our studies suggest that the metal ion plays a crucial role in the RNA pseudoknotting of VPK. In the stacking interior of S1 and S2, the Na+ ion is positioned in the major groove and interacts directly with the carbonyl group O6 of G28 and carbonyl group O4 of U13 in the wobble base pair U13:G28. The ion-RNA interactions in MMTV VPK not only stabilize the RNA pseudoknot but also modify the electrostatic properties of the nucleotides at the critical parts of the pseudoknot VPK.

Speeding up the dynamic algorithm for planar RNA folding.

The simplest dynamic algorithm for planar RNA folding searches for the maximum number of base pairs. The algorithm uses O(n3) steps. The more general case, where different weights (energies) are assigned to stacked base pairs and to the various types of single-stranded region topologies, requires a considerably longer computation time because of the partial backtracking involved. Limiting the loop size reduces the running time back to O(n3). Reduction in the number of steps in the calculations of the various RNA topologies has recently been suggested, thereby improving the time behavior. Here we show how a "jumping" procedure can be used to speed up the computation, not only for the maximal number of base pairs algorithm, but for the minimal energy algorithm as well.