Secondary structure maps of ribosomal RNA and DNA: I. Processing of Xenopus laevis ribosomal RNA and structure of single-stranded ribosomal DNA

Precursor and mature ribosomal RNA molecules from Xenopus laevis were examined by electron microscopy. A reproducible arrangement of hairpin loops was observed in these molecules. Maps based on this secondary structure were used to determine the arrangement of sequences in precursor RNA molecules and to identify the position of mature rRNAs within the precursors. A processing scheme was derived in which the 40 S rRNA is cleaved to 38 S RNA, which then yields 34 S plus 18 S RNA. The 34 S RNA is processed to 30 S, and finally to 28 S rRNA. The pathway is analogous to that of L-cell rRNA but differs from HeLa rRNA in that no 20 S rRNA intermediate was found. X. laevis 40 S rRNA (M_r = 2.7 × 10⁶) is much smaller than HeLa or L-cell 45 8 rRNA (M_r = 4.7 × 10⁶), but the arrangement of mature rRNA sequences in all precursors is very similar. Experiments with ascites cell 3′-exonuclease show that the 28 S region is located at or close to the 5′-end of the 40 S rRNA.Secondary structure maps were obtained also for single-stranded molecules of ribosomal DNA. The region in the DNA coding for the 40 S rRNA could be identified by its regular structure, which closely resembles that of the RNA. Regions corresponding to the 40 S RNA gene alternate with non-transcribed spacer regions along strands of rDNA. The latter have a large amount of irregular secondary structure and vary in length between different repeating units. A detailed map of the rDNA repeating unit was derived from these experiments.Optical melting studies are presented, showing that rRNAs with a high (G + C) content exhibit significant hypochromicity in the formamide/urea-containing solution that was used for spreading.     

A knowledge-based architecture for protein sequence analysis and structure prediction

Methods for analyzing the amino-acid sequence of a protein for the purposes of predicting its three-dimensional structure were systematically analyzed using knowledge engineering techniques. The resulting entities (data) and relations (processing methods and constraints) have been represented within a generalized dependency network consisting of 29 nodes and over 100 links. It is argued that such a representation meets the requirements of knowledge-based systems in molecular biology. This network is used as the architecture for a prototype knowledge-based system that simulates logically the processes used in protein structure prediction. Although developed specifically for applications in protein structure prediction, the network architecture provides a strategy for tackling the general problem of orchestrating and integrating the diverse sources of knowledge that are characteristic of many areas of science.     

APOBEC3G cytosine deamination hotspots are defined by both sequence context and single-stranded DNA secondary structure

Apolipoprotein B mRNA-editing, enzyme-catalytic, polypeptide-like 3G (i.e., APOBEC3G or A3G) is an evolutionarily conserved cytosine deaminase that potently restricts human immunodeficiency virus type 1 (HIV-1), retrotransposons and other viruses. A3G has a nucleotide target site specificity for cytosine dinucleotides, though only certain cytosine dinucleotides are ‘hotspots’ for cytosine deamination, and others experience little or no editing by A3G. The factors that define these critical A3G hotspots are not fully understood. To investigate how A3G hotspots are defined, we used an in vitro fluorescence resonance energy transfer-based oligonucleotide assay to probe the site specificity of A3G. Our findings strongly suggest that the target single-stranded DNA (ssDNA) secondary structure as well as the bases directly 3′ and 5′ of the cytosine dinucleotide are critically important A3G recognition. For instance, A3G cannot readily deaminate a cytosine dinucleotide in ssDNA stem structures or in nucleotide base loops composed of three bases. Single-stranded nucleotide loops up to seven bases in length were poor targets for A3G activity unless cytosine residues flanked the cytosine dinucleotide. Furthermore, we observed that A3G favors adenines, cytosines and thymines flanking the cytosine dinucleotide target in unstructured regions of ssDNA. Low cytosine deaminase activity was detected when guanines flanked the cytosine dinucleotide. Taken together, our findings provide the first demonstration that A3G cytosine deamination hotspots are defined by both the sequence context of the cytosine dinucleotide target as well as the ssDNA secondary structure. This knowledge can be used to better trace the origins of mutations to A3G activity, and illuminate its impact on processes such as HIV-1 genetic variation.     

Secondary structure prediction and determination of proteins--a review

The rapid increase in sequence data in combination with a greater understanding of the forces regulating protein structure has been the impetus for an upsurge in the development of theoretical prediction methods. These methods have afforded protein chemists the ability to identify and quantify the various secondary structures along the protein chain. Concurrently, various physico-chemical techniques have been developed such as nuclear Overhauser enhancement n.m.r. and laser Raman spectroscopy. In addition, traditional methods such as infrared and circular dichroism spectroscopy have been refined. Although both predictive and physico-chemical techniques are limited in the types of secondary structure they are capable of determining, they have provided valuable information with regards to protein folding and topology in the absence of X-ray data, and have formed the basis for the development of improved methods for secondary structure determination. This paper reviews some of the predictive and physico-chemical methods presently used to determine protein secondary structure.     

Structural modeling of sequence specificity by an autoantibody against single-stranded DNA

11F8 is a sequence-specific pathogenic anti-single-stranded (ss)DNA autoantibody isolated from a lupus prone mouse. Site-directed mutagenesis of 11F8 has shown that six binding site residues (R31VH, W33VH, L97VH, R98VH, Y100VH, and Y32VL) contribute 80% of the free energy for complex formation. Mutagenesis results along with intermolecular distances obtained from fluorescence resonance energy transfer were implemented here as restraints to model docking between 11F8 and the sequence-specific ssDNA. The model of the complex suggests that aromatic stacking and two sets of bidentate hydrogen bonds between binding site arginine residues (R31VH and R96VH) and loop nucleotides provide the molecular basis for high affinity and specificity. In part, 11F8 utilizes the same ssDNA binding motif of Y32VL, H91VL, and an aromatic residue in the third complementarity-determining region to recognize thymine-rich sequences as do two anti-ssDNA autoantibodies crystallized in complex with thymine. R31SVH is a dominant somatic mutation found in the J558 germline sequence that is implicated in 11F8 sequence specificity. A model of the mutant R31S11F8.ssDNA complex suggests that different interface contacts occur when serine replaces arginine 31 at the binding site. The modeled contacts between the R31S11F8 mutant and thymine are closely related to those observed in other anti-ssDNA binding antibodies, while we find additional contacts between 11F8 and ssDNA that involve amino acids not utilized by the other antibodies. These data-driven 11F8.ssDNA models provide testable hypotheses concerning interactions that mediate sequence specificity in 11F8 and the effects of somatic mutation on ssDNA recognition.     

Secondary structure prediction of interacting RNA molecules

Computational tools for prediction of the secondary structure of two or more interacting nucleic acid molecules are useful for understanding mechanisms for ribozyme function, determining the affinity of an oligonucleotide primer to its target, and designing good antisense oligonucleotides, novel ribozymes, DNA code words, or nanostructures. Here, we introduce new algorithms for prediction of the minimum free energy pseudoknot-free secondary structure of two or more nucleic acid molecules, and for prediction of alternative low-energy (sub-optimal) secondary structures for two nucleic acid molecules. We provide a comprehensive analysis of our predictions against secondary structures of interacting RNA molecules drawn from the literature. Analysis of our tools on 17 sequences of up to 200 nucleotides that do not form pseudoknots shows that they have 79% accuracy, on average, for the minimum free energy predictions. When the best of 100 sub-optimal foldings is taken, the average accuracy increases to 91%. The accuracy decreases as the sequences increase in length and as the number of pseudoknots and tertiary interactions increases. Our algorithms extend the free energy minimization algorithm of Zuker and Stiegler for secondary structure prediction, and the sub-optimal folding algorithm by Wuchty et al. Implementations of our algorithms are freely available in the package MultiRNAFold.     

Secondary structure prediction of seed storage proteins

The comparison of partial primary structure of seed storage proteins leads to show homologies inside of each considered family (Legume seed legumins and cereal prolamins). Predicted secondary structures deduced from the presently known sequences also exhibit considerable homologies, which implies a severe conservatism of these proteins. Short repetitive segments of sequence of 5-20 residues are frequently occurring and give rise to the prediction either of beta-structure (or alpha-helix) bonded by beta-turns or of successive beta-turns. The latter conformation, which would be able to form a helicoidal arrangement, could contribute to a maximal packing of the protein molecules inside of the subcellular organelles (protein bodies) within which they are confined. As the only known function of seed storage proteins is to provide amino acids to the embryo, it is suggested that their ability to occupy a minimal volume is actually a reasonable explanation of their extreme conservatism in the course of evolution.     

Extraordinary stable structure of short single-stranded DNA fragments containing a specific base sequence: d(GCGAAAGC).

Short single-stranded DNA fragments carrying a GCGAAAGC sequence were found to move unexpectedly faster than other fragments of the same length in electrophoresis on a polyacrylamide gel containing a denaturing agent. The fragments were noted to have a stable structure even in 7M urea solution, but the stability cannot be explained simply on the basis of base pair formation alone. Physical characterization of the GCGAAAGC fragment indicated that it takes a hairpin-like structure in spite of the short chain length with only two G-C base pairs, comprised of GCG and AAAGC subsegments, each possessing a helical configuration independent of the others. Some biological implications of this unusual structure are discussed.     

PROTINFO: Secondary and tertiary protein structure prediction

Information about the secondary and tertiary structure of a protein sequence can greatly assist biologists in the generation and testing of hypotheses, as well as design of experiments. The PROTINFO server enables users to submit a protein sequence and request a prediction of the three-dimensional (tertiary) structure based on comparative modeling, fold generation and de novo methods developed by the authors. In addition, users can submit NMR chemical shift data and request protein secondary structure assignment that is based on using neural networks to combine the chemical shifts with secondary structure predictions. The server is available at http://protinfo.compbio.washington.edu.     

Protein local structure prediction from sequence

A basis set of protein canonical fragments, or centroids, represents the range of local structure found in globular proteins. We develop a methodology to predict centroids from the amino acid sequence. The predictor gives the probability of each centroid in the basis set, at each loci along the backbone. The predictor selects the best-fit centroid at about 40% of the loci. The predicted probabilities are accurate and can be used to judge the confidence of each centroid prediction. For example, when filtering out centroids with <0.50 probability, the predictor is 65% accurate, although such high-probability centroids occur at only 28% of the loci. Centroids with high probability can be interpreted as segments that are highly influenced by the amino acid sequence, whereas centroids with low probability can be interpreted as segments that are more likely influenced by tertiary contacts. Low-resolution, starting point structures, can be generated by fitting the predicted centroids together.     

Secondary structure prediction for aligned RNA sequences

Most functional RNA molecules have characteristic secondary structures that are highly conserved in evolution. Here we present a method for computing the consensus structure of a set aligned RNA sequences taking into account both thermodynamic stability and sequence covariation. Comparison with phylogenetic structures of rRNAs shows that a reliability of prediction of more than 80% is achieved for only five related sequences. As an application we show that the Early Noduline mRNA contains significant secondary structure that is supported by sequence covariation.     

Secondary structure prediction with support vector machines

MOTIVATION: A new method that uses support vector machines (SVMs) to predict protein secondary structure is described and evaluated. The study is designed to develop a reliable prediction method using an alternative technique and to investigate the applicability of SVMs to this type of bioinformatics problem. METHODS: Binary SVMs are trained to discriminate between two structural classes. The binary classifiers are combined in several ways to predict multi-class secondary structure. RESULTS: The average three-state prediction accuracy per protein (Q(3)) is estimated by cross-validation to be 77.07 +/- 0.26% with a segment overlap (Sov) score of 73.32 +/- 0.39%. The SVM performs similarly to the 'state-of-the-art' PSIPRED prediction method on a non-homologous test set of 121 proteins despite being trained on substantially fewer examples. A simple consensus of the SVM, PSIPRED and PROFsec achieves significantly higher prediction accuracy than the individual methods.     

杨树木质素合成酶CCR基因的序列分析及蛋白结构预测

利用RT—PCR从欧美杨107次生木质部中克隆出-961bp的CCR基因片段。通过生物信息学软件对该序列的核苷酸序列、拟翻译的氨基酸序列的疏水性、残基带电量及表面暴露区、蛋白质二级结构、亚细胞定位及三维结构等进行了初步分析预测。结果表明该CCR基因含一个编码301个氨基酸的完整开放阅读框,其成熟蛋白为亲水性的、主要存在于细胞膜,具有大多数植物CCR蛋白普遍存在的KNWYCYGK的保守性基序,其二级结构中共包含12个α螺旋,20个β折叠,11个卷曲,并构建了其三维结构图。     

Secondary structure prediction of 11 mammalian growth hormones

The secondary structure of 11 mammalian growth hormones has been predicted by combining five different methods. Three long helical regions located around residues 20, 120, and 170 constitute the most prominent common feature in the species studied. The strong amphiphilic character of these helices suggests that they can play an important role in protein folding or stability.     

Secondary structure prediction of human salivary proline-rich proteins

Conformations associated with secondary structure in human salivary proline-rich proteins A (PRPA), C (PRPC), P-D and P-E were predicted by analysis of their respective hydrophobicity profiles by computer programming. Structurally, PRPA and PRPC would present a globular head and a tail that consists of type 3(10) polyproline helices. P-D and P-E would be fibrilar molecules with helical zones of the polyproline 3(10) type. Alternatively for PRPA and PRPC, the head and tail would form one globular domain with the tail folding upon itself at places where random coils occur.     

Secondary structure prediction and unrefined tertiary structure prediction for cyclin A,B, and D

We present heuristic-based predictions of the secondary and tertiary structures of the cyclins A, B, and D, representatives of the cyclin superfamily. The list of suggested constraints for tertiary structure assembly was left unrefined in order to submit this report before an announced crystal structure for cyclin A becomes available. To predict these constraints, a master sequence alignment over 270 positions of cyclin types A, B, and D was adjusted based on individual secondary structure predictions for each type. We used new heuristics for predicting aromatic residues at protein-protein interfaces and to identify sequentially distinct regions in the protein chain that cluster in the folded structure. The boundaries of two conjectured domains in the cyclin fold were predicted based on experimental data in the literature. The domain that is important for interaction of the cyclins with cyclin-dependent kinases (CDKs) is predicted to contain six helices; the second domain in the consensus model contains both helices and a β-sheet that is formed by sequentially distant regions in the protein chain. A plausible phosphorylation site is identified. This work represents a blinded test of the method for prediction of secondary and, to a lesser extent, tertiary structure from a set of homologous protein sequences. Evaluation of our predictions will become possible with the publication of the announced crystal structure.     

Thermostable alanine racemase from Bacillus stearothermophilus: DNA and protein sequence determination and secondary structure prediction

The nucleotide sequence of the alanine racemase (EC 5.1.1.1) gene from a thermophile, Bacillus stearothermophilus, was determined by the dideoxy chain termination method with universal and synthetic site-specific primers. The amino acid sequence of the enzyme predicted from the nucleotide sequence was confirmed by peptide sequence information derived from the N-terminal amino acid residues and several tryptic fragments. The alanine racemase gene consists of 1158 base pairs encoding a protein of 386 amino acid residues; the molecular weight of the apoenzyme is estimated as 43,341. The racemase gene of B. stearothermophilus has a closely similar size (1158 vs 1167 base pairs) to that of the gene of a mesophile, B. subtilis, but shows a higher preference for codons ending in G or C. A comparison of the amino acid sequence with those of Bacillus subtilis and Salmonella typhimurium dadB and alr enzymes revealed overall sequence homologies of 31-54%, including an identical octapeptide bearing the pyridoxal 5'-phosphate binding site. Although the residues common in the four racemases are not continuously arrayed, these constitute distinct domains and their hydropathy profiles are very similar. The secondary structure of B. stearothermophilus alanine racemase was predicted from the results obtained by theoretical analysis and circular dichroism measurement.