Dependence among sites in RNA evolution

Although probabilistic models of genotype (e.g., DNA sequence) evolution have been greatly elaborated, less attention has been paid to the effect of phenotype on the evolution of the genotype. Here we propose an evolutionary model and a Bayesian inference procedure that are aimed at filling this gap. In the model, RNA secondary structure links genotype and phenotype by treating the approximate free energy of a sequence folded into a secondary structure as a surrogate for fitness. The underlying idea is that a nucleotide substitution resulting in a more stable secondary structure should have a higher rate than a substitution that yields a less stable secondary structure. This free energy approach incorporates evolutionary dependencies among sequence positions beyond those that are reflected simply by jointly modeling change at paired positions in an RNA helix. Although there is not a formal requirement with this approach that secondary structure be known and nearly invariant over evolutionary time, computational considerations make these assumptions attractive and they have been adopted in a software program that permits statistical analysis of multiple homologous sequences that are related via a known phylogenetic tree topology. Analyses of 5S ribosomal RNA sequences are presented to illustrate and quantify the strong impact that RNA secondary structure has on substitution rates. Analyses on simulated sequences show that the new inference procedure has reasonable statistical properties. Potential applications of this procedure, including improved ancestral sequence inference and location of functionally interesting sites, are discussed.     

Prediction of protein secondary structure and active sites using the alignment of homologous sequences

The prediction of protein secondary structure (alpha-helices, beta-sheets and coil) is improved by 9% to 66% using the information available from a family of homologous sequences. The approach is based both on averaging the Garnier et al. (1978) secondary structure propensities for aligned residues and on the observation that insertions and high sequence variability tend to occur in loop regions between secondary structures. Accordingly, an algorithm first aligns a family of sequences and a value for the extent of sequence conservation at each position is obtained. This value modifies a Garnier et al. prediction on the averaged sequence to yield the improved prediction. In addition, from the sequence conservation and the predicted secondary structure, many active site regions of enzymes can be located (26 out of 43) with limited over-prediction (8 extra). The entire algorithm is fully automatic and is applicable to all structural classes of globular proteins.     

Comparison between the toluidine blue stain and the Feulgen reaction for evaluation of rabbit sperm chromatin condensation and their relationship with sperm morphology

Sperm chromatin alteration is an important feature that can affect fertility of the male rabbit. This study compared toluidine blue staining with Feulgen reaction (as methods for evaluating chromatin alteration) and investigated the relationship between sperm morphology and chromatin alteration. Seven hundred rabbit ejaculates of animals with unknown fertility were used. Primary and secondary morphological sperm abnormalities were evaluated in semen smears with phase-contrast microscopy. Chromatin alterations were evaluated in semen smears stained with toluidine blue (pH 4.0 and 5.0) and with the Feulgen reaction. While the three methods were equally efficacious for identification of chromatin alterations, toluidine blue staining was more appropriate to characterize the intensity of chromatin alterations. The correlation between primary sperm defects and chromatin alteration was high and positive, suggesting that sperm chromatin structure affected sperm head morphology. The correlation between secondary sperm defects and chromatin alteration was also positive, but lower. The final chromatin compaction occurs in the epididymus, where secondary sperm defects originate. Therefore, the causes of secondary sperm defects could also intervene with final chromatin compaction. In summary, the toluidine blue stain was an effective means of evaluating the sperm chromatin alteration in rabbit spermatozoa.     

The purple to blue transition of bacteriorhodopsin is accompanied by a loss of the hexagonal lattice and a conformational change

X-ray diffraction measurements show that in contrast to the purple membrane, the bacteriorhodopsin molecules are not organized in a hexagonal lattice in the deionized blue membrane. Addition of Ca2+ restores both the purple color and the normal (63 A) hexagonal protein lattice. In the blue state, the circular dichroism spectrum in the visible has the typical exciton features indicating that a trimeric structure is retained. Time-resolved linear dichroism measurements show that the blue patch rotates in aqueous suspension with a mean correlation time of 11 ms and provide no evidence for rotational mobility of bacteriorhodopsin within the membrane. The circular dichroism spectra of the blue and the Ca2+-regenerated purple state in the far-UV are different, indicating a small change in secondary structure. The thermal stability of the blue membrane is much smaller than that of the purple membrane. At pH 5.0, the irreversible denaturation transition of the blue form has a midpoint at 61 degrees C. The photocycle of the blue membrane (lambda ex 590 nm) has an L intermediate around 540 nm whose decay is slowed down into the millisecond time range (5 ms). Light-dark adaptation in the blue membrane is rapid with an exponential decay time of 38 s at 25 degrees C. The purple to blue transition apparently involves a conformational change in the protein leading to a change in the aggregation state from a highly ordered and stable hexagonal lattice to a disordered array of thermally more labile trimers.(ABSTRACT TRUNCATED AT 250 WORDS)     

Structure and reactivity of type-I copper sites

Summary Biological electron transfer is not well understood. The question is addressed in this contribution with reference to the so-called blue copper proteins, each of which has a single copper atom at its active centre. The redox activity (as probed by the electron self exchange reaction) of the Cu centre seems not to be affected. The electron self exchange reaction is known to proceed through His-117, and the hydrophobic patch is most important in the formation of the azurin/azurin encounter complex. Ph effects have not been observed on the three-dimensional structure ofA. denitrificans azurin, which may indicate that if present at all these have no direct physiological implications. Mutants are in process of construction.     

Analysis of catalytic residues in enzyme active sites

We present an analysis of the residues directly involved in catalysis in 178 enzyme active sites. Specific criteria were derived to define a catalytic residue, and used to create a catalytic residue dataset, which was then analysed in terms of properties including secondary structure, solvent accessibility, flexibility, conservation, quaternary structure and function. The results indicate the dominance of a small set of amino acid residues in catalysis and give a picture of a general active site environment. It is hoped that this information will provide a better understanding of the molecular mechanisms involved in catalysis and a heuristic basis for predicting catalytic residues in enzymes of unknown function.     

Conformational properties of azurin in solution as determined from resolution-enhanced Fourier-transform infrared spectra

Infrared spectra of the blue copper protein azurin and of apoazurin from P. fluorescens were obtained in aqueous solution. Using resolution enhancement procedures, a number of component bands were identified in the region of the amide I mode, and these bands were assigned to various components of protein secondary structure. A quantitative analysis of these infrared spectra indicates that the secondary structure of P. fluorescens azurin in solution is very similar to those determined previously by X-ray diffraction for the crystals of azurins from other bacterial species. The major components of this structure are beta strands and turns. Infrared spectra also evince a remarkable thermal stability of the native azurin. A significant unfolding of the protein could only be detected at temperatures above approximately 76 degrees C. While the secondary structure of apoazurin is practically indistinguishable from that of the native protein at room temperature, the thermal stability of the apo form is significantly reduced.     

Blue flickers of hope: secondary structure, dynamics, and putative dimerization interface of the blue-light receptor YtvA from Bacillus subtilis

Bacillus subtilis is capable of responding to various kinds of extracellular, potentially harmful stimuli via a stress response pathway, which involves a signal transduction and integration hub, the stressosome, and finally leads to activation of σ(B). One of the different signals initiating the underlying phosphorylation cascade is blue light. While it is known that the bacterial photoreceptor YtvA is responsible for blue light detection, the intramolecular activation mechanism and the structure of this multidomain protein are unknown. Using solution NMR spectroscopy, we have obtained a near complete backbone assignment of the full-length protein. More importantly, we report relaxation data and data on the solvent accessibility of full-length YtvA in the dark state which are interpreted with respect to secondary structure, the mobility, and the quaternary structure of the protein. Finally, we show that YtvA adopts an elongated domain orientation with LOV-LOV and STAS-STAS interactions on either side.     

Characterization of metal ion-binding sites in bacteriorhodopsin

We have investigated the effects of the binding of various metal ions to cation-free bacteriorhodopsin ("blue membrane"). The following have been measured: shift of the absorption maximum from 603 to 558 nm (blue to purple transition), binding isotherms, the release of H+ upon binding, and the decay of the deprotonated intermediate of the photocycle, M412. We find that all cations of the lanthanide series, as well as the alkali and alkali earth metals earlier investigated, are able to bring about the absorption shift, whereas Hg2+ and Pt4+ are not. Sigmoidal spectroscopic titration curves and nonsigmoidal binding curves suggest that there are two high affinity sites for cations in bacteriorhodopsin. Binding to the site with the second highest affinity is responsible for the absorption shift. Divalent cation binding to blue membrane causes release of about six protons, whereas higher numbers of protons are released by trivalent cations, suggesting that the shift of absorption maximum involves proton release from carboxyl group(s). The metal ion bound to this site must be surrounded by carboxyl oxygen atoms acting together as a multidentate ligand with a specific geometry because multivalent ions are effective only when capable of octahedral coordination. Lanthanide ions dramatically inhibit M412 decay at pH above 6.3, an effect probably due to binding to lipid phosphoryl groups.     

A possible explanation for the multiple polyadenylation sites in transcripts coding for a winged-bean leghemoglobin

Five different copy DNA clones coding for the same leghemoglobin were isolated from a winged-bean (Psophocarpus tetragonolobus L.) nodule library. Although identical in sequence, they each possess a different side of polyadenylation located 93–128 nucleotides downstream of two overlapping AAUAAA putative signal sequences. By analysis of the untranslated 3′ ends, a potential mRNA secondary structure can be predicted which could explain the observed polyadenylation heterogeneity. The structure is a size-variable hairpin, creating a net topological distance of 25–27 nucleotides between the canonical signal sequence and the different polyadenylation sites observed. We suggest that this type of variable secondary structure could be one among other causes that determines the apparent flexibility of plant polyadenylation. It could also confer particular properties to the mRNA in relation to stability, translation efficiency and-or nuclear export.     

Native-like tertiary structure in the Mucor miehei lipase molten globule state obtained at low pH

Studies on the acid-induced denaturation of Mucor miehei lipase (E.C. 3.1.1.3) were performed by circular dichroism (CD) spectroscopy, fluorescence emission spectroscopy and binding of hydrophobic dye, 1-anilino 8-naphthalenesulfonic acid (ANS). Acid denaturation of the lipase showed loss of secondary structure and alterations in the tertiary structure in the pH range 4 to 2 and 7 to 2 respectively, suggesting that the lipase exists as an acid-unfolded state approximately pH 2.0. A further decrease in pH (from 2.0 to 1.0) resulted in a second transition, which corresponded to the formation of both secondary and tertiary structures. The acid unfolded state at around pH 2.0 has been characterized by significant loss of secondary structure and a small increase in fluorescence intensity with a blue shift of 2 nm, indicating shift of tryptophan residues to less polar environment. Interestingly, the lipase at pH 1.0 exhibits characteristics of molten globule, such as enhanced binding of hydrophobic dye (ANS), native-like secondary structure and slightly altered tryptophanyl environments. That the molten globule of the lipase at pH 1.0 also possesses native-like tertiary structure is an interesting observation made for this lipase.     

Statistical prediction of single-stranded regions in RNA secondary structure and application to predicting effective antisense target sites and beyond

Single-stranded regions in RNA secondary structure are important for RNA–RNA and RNA–protein interactions. We present a probability profile approach for the prediction of these regions based on a statistical algorithm for sampling RNA secondary structures. For the prediction of phylogenetically-determined single-stranded regions in secondary structures of representative RNA sequences, the probability profile offers substantial improvement over the minimum free energy structure. In designing antisense oligonucleotides, a practical problem is how to select a secondary structure for the target mRNA from the optimal structure(s) and many suboptimal structures with similar free energies. By summarizing the information from a statistical sample of probable secondary structures in a single plot, the probability profile not only presents a solution to this dilemma, but also reveals ‘well-determined’ single-stranded regions through the assignment of probabilities as measures of confidence in predictions. In antisense application to the rabbit β-globin mRNA, a significant correlation between hybridization potential predicted by the probability profile and the degree of inhibition of in vitro translation suggests that the probability profile approach is valuable for the identification of effective antisense target sites. Coupling computational design with DNA–RNA array technique provides a rational, efficient framework for antisense oligonucleotide screening. This framework has the potential for high-throughput applications to functional genomics and drug target validation.     

Secondary structure in heterogeneous nuclear RNA: involvement of regions from repeated DNA sites

Heterogeneous nuclear RNA was found to contain regions of secondary structure based on a relative resistance to nuclease treatment compared with mRNA or poliovirus RNA and a shift in density toward double-stranded RNA early in the course of nuclease digestion. The regions involved in this secondary structure are enriched for RNA segments transcribed from repeated sites in the DNA. Thus, to maximize hybridization to repetitive sites heterogeneous nuclear RNA molecules must be both denatured and fragmented. Some of the self-complementary regions in heterogeneous nuclear RNA are released by alkali denaturation and fragmentation below 1500 nucleotides but maximum release is not achieved until fragmentation below 500 nucleotides. These results indicate that these self-complementary regions (“loops” plus “stems”) are mainly below 500 nucleotides in length.     

Characterizing the microenvironment surrounding protein sites.

Sites are microenvironments within a biomolecular structure, distinguished by their structural or functional role. A site can be defined by a three-dimensional location and a local neighborhood around this location in which the structure or function exists. We have developed a computer system to facilitate structural analysis (both qualitative and quantitative) of biomolecular sites. Our system automatically examines the spatial distributions of biophysical and biochemical properties, and reports those regions within a site where the distribution of these properties differs significantly from control nonsites. The properties range from simple atom-based characteristics such as charge to polypeptide-based characteristics such as type of secondary structure. Our analysis of sites uses non-sites as controls, providing a baseline for the quantitative assessment of the significance of the features that are uncovered. In this paper, we use radial distributions of properties to study three well-known sites (the binding sites for calcium, the milieu of disulfide bridges, and the serine protease active site). We demonstrate that the system automatically finds many of the previously described features of these sites and augments these features with some new details. In some cases, we cannot confirm the statistical significance of previously reported features. Our results demonstrate that analysis of protein structure is sensitive to assumptions about background distributions, and that these distributions should be considered explicitly during structural analyses.     

奥利亚罗非鱼Piwi基因的克隆与生物信息学分析

Piwi基因是影响鱼类生殖细胞发育的重要因子。通过设计特异引物,以奥利亚罗非鱼的精巢RNA为模板,通过PCR扩增和克隆测序,获得了两种基因序列,分别长2 571 bp和3 204 bp。通过生物信息学分析,表明这两个片段与尼罗罗非鱼的Piwil-1和Piwil-2相似性都达到99%。通过翻译,Piwil-1和Piwil-2的蛋白序列具有典型的PAZ和Piwi结构域。因此,这两个基因为奥利亚罗非鱼的Piwil-1和Piwil-2基因。另外,生物信息学分析也表明Piwil-1和Piwil-2在一级结构和二级结构上十分相似,说明在奥利亚罗非鱼基因组中这两个基因分化时间则较短。本研究为后续研究鱼类生殖机理和罗非鱼进化研究奠定了基础。     

Spectral and physical characterization of the inverted terminal repeat DNA structure from adenoassociated virus 2

An oligodeoxynucleotide (ODN) that includes elements found in secondary structures at the 5'- and 3'- ends of adenoassociated virus 2 virion DNA was synthesized by ligation of three overlapping ODNs. The most stable secondary structure was calculated to be branched, with a 61 bp duplex stem, terminating in a three-way junction with 9 bp arms. The electrophoretic mobility of the ODN is slower than expected for normal duplex DNA of the same size, suggesting a bent or branched conformation. CD spectra indicate that the ITR structure is largely B form DNA, although there is a slight blue shift compared to the spectra of the isolated stem and loop elements. Thermal melting experiments indicate that the hairpin is significantly more stable than the isolated stem and loop elements. Singular value decomposition of UV spectra obtained as a function of temperature indicates that four species contribute to changes in the spectra upon denaturation, indicating that the melting is not a simple two-state process. Characterization of the branched ODN by differential scanning calorimetry permits estimation of the enthalpy of melting by a model-free analysis, yielding DeltaHcal= 614 kcal mol-1. This value agrees with the enthalpy computed for the most stable secondary structure.     

Detection of initiation sites in protein folding of the four helix bundle ACBP by chemical shift analysis

A simple alternative method for obtaining "random coil" chemical shifts by intrinsic referencing using the protein's own peptide sequence is presented. These intrinsic random coil backbone shifts were then used to calculate secondary chemical shifts, that provide important information on the residual secondary structure elements in the acid-denatured state of an acyl-coenzyme A binding protein. This method reveals a clear correlation between the carbon secondary chemical shifts and the amide secondary chemical shifts 3-5 residues away in the primary sequence. These findings strongly suggest transient formation of short helix-like segments, and identify unique sequence segments important for protein folding.     

Prediction of the secondary structure and functional sites of major histocompatibility complex molecules.

A method for the theoretical prediction of the antigenic determinants and the antigen-interactive receptor sites of immunological proteins from their primary structure would constitute a useful tool for their study. Such a method developed in this laboratory uses hydrophilicity, accessibility, flexibility, and recognition profiles, together with the predicted secondary structure (alpha-helices, beta-sheets, and turns). The secondary structure is determined by a modification of the method of Lim (1974), as described below. A study of human and mouse class I and class II major histocompatibility complex (MHC) antigens, central to the regulation of immune responses and to the phenomenon of graft rejection, was carried out using the above method. Comparison of the predictions with some of the available experimental and theoretical information supports the validity and usefulness of the approach.     

An algorithm for predicting protein-protein interaction sites: Abnormally exposed amino acid residues and secondary structure elements

Multiprotein systems mediate most regulatory processes in living organisms. Although the structures of the individual proteins are often defined, less is known of the structures of multiprotein systems. Computational methods for predicting interfaces, using evolutionary conservation and/or physicochemical data, have been developed. Here we consider the use of solvent accessibility, residue propensity, and hydrophobicity, in conjunction with secondary structure data, as prediction parameters. We analyze the influence of residue type and secondary structure on solvent accessibility and define a measure of "relative exposedness." Clustering abnormally high scoring residues provides a basis for predicting interaction sites. The analysis is extended to investigate abnormally exposed secondary structure elements, particularly beta-sheet strands. We show that surface-exposed beta-strands lacking protective features are more likely to be found at protein-protein interfaces, allowing us to create an algorithm with approximately 68% and approximately 75% accuracy in differentiating between interacting and edge strands in isolated beta-strands and beta-sheet strands, respectively. These methods of identifying abnormally exposed surface regions are combined in an algorithm, which, on a data set of 77 unbound and disjoint (single chain extracted from complex) structures, predicts 79% of the protein-protein interfaces correctly. If enzyme-inhibitor complexes, where the inhibitor mimics a nonprotein substrate, are excluded, the accuracy increases to 85%.     

Optimal mutation sites for PRE data collection and membrane protein structure prediction

Nuclear magnetic resonance paramagnetic relaxation enhancement (PRE) measures long-range distances to isotopically labeled residues, providing useful constraints for protein structure prediction. The method usually requires labor-intensive conjugation of nitroxide labels to multiple locations on the protein, one at a time. Here a computational procedure, based on protein sequence and simple secondary structure models, is presented to facilitate optimal placement of a minimum number of labels needed to determine the correct topology of?a helical transmembrane protein. Tests on DsbB (four helices) using just one label lead to correct topology predictions in four of five cases, with the predicted structures <6 ? to the native structure. Benchmark results using simulated PRE data show that we can generally predict the correct topology for five and six to seven helices using two and three labels, respectively, with an average success rate of 76% and structures of similar precision. The results show promise in facilitating experimentally constrained structure prediction of membrane proteins.