Assessing the reliability of RNA folding using statistical mechanics

We have analyzed the base-pairing probability distributions of 16 S and 16 S-like, and 23 S and 23 S-like ribosomal RNAs of Archaea, Bacteria, chloroplasts, mitochondria and Eukarya, as predicted by the partition function approach for RNA folding introduced by McCaskill. A quantitative analysis of the reliability of RNA folding is done by comparing the base-pairing probability distributions with the structures predicted by comparative sequence analysis (comparative structures). We distinguish two factors that show a relationship to the reliability of RNA minimum free energy structure. The first factor is the dominance of one particular base-pair or the absence of base-pairing for a given base within the base-pairing probability distribution (BPPD). We characterize the BPPD per base, including the probability of not base-pairing, by its Shannon entropy (S). The S value indicates the uncertainty about the base-pairing of a base: low S values result from BPPDs that are strongly dominated by a single base-pair or by the absence of base-pairing. We show that bases with low S values have a relatively high probability that their minimum free energy (MFE) structure corresponds to the comparative structure. The BPPDs of prokaryotes that live at high temperatures (thermophilic Archaea and Bacteria) have, calculated at 37 degrees C, lower S values than the BPPDs of prokaryotes that live at lower temperatures (mesophilic and psychrophilic Archaea and Bacteria). This reflects an adaptation of the ribosomal RNAs to the environmental temperature. A second factor that is important to consider with regard to the reliability of MFE structure folding is a variable degree of applicability of the thermodynamic model of RNA folding for different groups of RNAs. Here we show that among the bases that show low S values, the Archaea and Bacteria have similar, high probabilities (0.96 and 0.94 in 16 S and 0.93 and 0.91 in 23 S, respectively) that the MFE structure corresponds to the comparative structure. These probabilities are lower in the chloroplasts (16 S 0.91, 23 S 0.79), mitochondria (16 S-like 0.89, 23 S-like 0.69) and Eukarya (18 S 0.81, 28 S 0.86).     

Prediction of RNA secondary structure based on helical regions distribution

MOTIVATION: RNAs play an important role in many biological processes and knowing their structure is important in understanding their function. Due to difficulties in the experimental determination of RNA secondary structure, the methods of theoretical prediction for known sequences are often used. Although many different algorithms for such predictions have been developed, this problem has not yet been solved. It is thus necessary to develop new methods for predicting RNA secondary structure. The most-used at present is Zuker's algorithm which can be used to determine the minimum free energy secondary structure. However many RNA secondary structures verified by experiments are not consistent with the minimum free energy secondary structures. In order to solve this problem, a method used to search a group of secondary structures whose free energy is close to the global minimum free energy was developed by Zuker in 1989. When considering a group of secondary structures, if there is no experimental data, we cannot tell which one is better than the others. This case also occurs in combinatorial and heuristic methods. These two kinds of methods have several weaknesses. Here we show how the central limit theorem can be used to solve these problems. RESULTS: An algorithm for predicting RNA secondary structure based on helical regions distribution is presented, which can be used to find the most probable secondary structure for a given RNA sequence. It consists of three steps. First, list all possible helical regions. Second, according to central limit theorem, estimate the occurrence probability of every helical region based on the Monte Carlo simulation. Third, add the helical region with the biggest probability to the current structure and eliminate the helical regions incompatible with the current structure. The above processes can be repeated until no more helical regions can be added. Take the current structure as the final RNA secondary structure. In order to demonstrate the confidence of the program, a test on three RNA sequences: tRNAPhe, Pre-tRNATyr, and Tetrahymena ribosomal RNA intervening sequence, is performed. AVAILABILITY: The program is written in Turbo Pascal 7.0. The source code is available upon request. CONTACT: Wujj@nic.bmi.ac.cn or Liwj@mail.bmi.ac.cn     

Kinetic and thermodynamic study of the bacteriorhodopsin photocycle over a wide pH range

The photocycle of bacteriorhodopsin and its thermodynamic parameters were studied in the pH range of 4.5-9. Measurements were performed at five different wavelengths (410, 500, 570, 610, and 650 nm), in the time interval 300 ns to 0.5 s, at six temperatures between 5 and 30 degreesC. Data were fitted to different photocycle models. The sequential model with reversible reactions gave a good fit, and the linear character of the Eyring plots was fulfilled. The parallel model with unidirectional reactions gave a poor fit, and the Eyring plot of the rate constants did not follow the expected linear behavior. When a parallel model with reversible reactions, which has twice as many free parameters as the sequential model, was considered, the quality of the fit did not improve and the Eyring plots were not linear. The sequential model was used to determine the thermodynamic activation parameters (activation enthalpy, entropy, and free energy) of the transitions and the free energy levels of the intermediates. pH dependence of the parameters revealed details of the transitions between the intermediates: the transitions M1 to M2 and N to O disclosed a large entropy increase, which could be interpreted as a loosening of the protein structure. The pH dependence of the energy levels explains the disappearance of intermediate O at high pH. A hypothesis is proposed to interpret the relation between the observed pKa of the photocycle energetics and the role of several amino acids in the protein.     

Process Simulation and Control Optimization of a Blast Furnace Using Classical Thermodynamics Combined to a Direct Search Algorithm

Several numerical approaches have been proposed in the literature to simulate the behavior of modern blast furnaces: finite volume methods, data-mining models, heat and mass balance models, and classical thermodynamic simulations. Despite this, there is actually no efficient method for evaluating quickly optimal operating parameters of a blast furnace as a function of the iron ore composition, which takes into account all potential chemical reactions that could occur in the system. In the current study, we propose a global simulation strategy of a blast furnace, the 5-unit process simulation. It is based on classical thermodynamic calculations coupled to a direct search algorithm to optimize process parameters. These parameters include the minimum required metallurgical coke consumption as well as the optimal blast chemical composition and the total charge that simultaneously satisfy the overall heat and mass balances of the system. Moreover, a Gibbs free energy function for metallurgical coke is parameterized in the current study and used to fine-tune the simulation of the blast furnace. Optimal operating conditions and predicted output stream properties calculated by the proposed thermodynamic simulation strategy are compared with reference data found in the literature and have proven the validity and high precision of this simulation.     

Stability of the TAR RNA and its complex with arginine

Most stable secondary structures and their stabilization energies of the TAR RNA with +1 to +104 nucleotide-sequence region were calculated at different temperatures by using thermodynamic parameters for RNA structure prediction. The most stable secondary structure has one bulge and one loop within the region of +20 to +40 nucleotide sequence, and its stabilization energy at 37 degrees C was -46.3 kcal mol-1. The interaction of a TAR bulge oligomer (TARBO) with arginine (Arg) which was in the binding site of a Tat protein was also investigated by CD measurements. The addition of Arg did not affect the CD spectrum of TARBO. The result was different from that of the RNA oligomer with both bulge and loop.     

Prevalence of penile human papillomavirus DNA in husbands of women with and without cervical neoplasia: a study in Spain and Colombia

Quinidine is an active antimalarial compound extracted from the bark of Cinchona trees. The activity differences among structurally related molecules appears to depend on the absolute stereochemistry of some functional groups, a result that stimulated a detailed conformational analysis of these molecules of biological interest. In the present study, the potential energy surface (PES) for the antimalarial agent quinidine (C20H24O2N2) has been comprehensively investigated using the molecular mechanics (MM) and quantum mechanical semiempirical AM1 and PM3 methods. Six distinct minimum energy conformations were located on the multidimensional PES and also characterized as true minima through harmonic frequency analysis. The relative stabilities and thermodynamic properties are reported. The coexistence of different conformers is discussed for the first time in the literature based on the transition state (TS) structures located on the PES for the quinidine molecule. The theoretical results reported in the present study are in agreement with the experimental proposal, based on NMR data, that there are two conformations existing in solution for the quinidine molecule.     

Efficacy of clean v sterile surgical prep kits

Frameshift and readthrough sites within retroviral messenger RNAs are often followed by nucleotide sequences that have the potential to form pseudoknot structures. In the work presented here, NMR methods were used to characterize the base-pairings and structural features of the RNA pseudoknot downstream of the gag-pro frameshift site of simian retrovirus type-1 (SRV-1) and a functional mutant of the SRV-1 pseudoknot. Evidence is presented that these pseudoknots contain two A-form helical stems of six base-pairs each, connected by two loops, in a classic H-type pseudoknot topology. A particularly interesting feature is that the shorter of the two connecting loops, loop 1, consists of only a single adenosine nucleotide that spans the major groove of stem 2. In this respect, the frameshift-associated pseudoknots are structurally similar to the pseudoknot within the gene 32 mRNA of bacteriophage T2, previously characterized by NMR methods. Despite having similar nucleotide sequences, the solvent exchange rates of the imino protons at the junction of the helical stems in the wild-type and mutant frameshifting pseudoknots differ from each other and from the bacteriophage T2 pseudoknot. The implications of this finding are discussed.     

Prediction of RNA Secondary Structure Based on Particle Swarm Optimization

A novel method for the prediction of RNA secondary structure was proposed based on the particle swarm optimization(PSO).PSO is known to be effective in solving many different types of optimization problems and known for being able to approximate the global optimal results in the solution space.We designed an efficient objective function according to the minimum free energy,the number of selected stems and the average length of selected stems.We calculated how many legal stems there were in the sequence,and selected some of them to obtain an optimal result using PSO in the right of the objective function.A method based on the improved particle swarm optimization(IPSO)was proposed to predict RNA secondary structure,which consisted of three stages.The first stage was applied to encoding the source sequences,and to exploring all the legal stems.Then,a set of encoded stems were created in order to prepare input data for the second stage.In the second stage,IPSO was responsible for structure selection.At last,the optimal result was obtained from the secondary structures selected via IPSO.Nine sequences from the comparative RNA website were selected for the evaluation of the proposed method.Compared with other six methods,the proposed method decreased the complexity and enhanced the sensitivity and specificity on the basis of the experiment results.     

Specific site selection in RNA resulting from a combination of nonspecific secondary structure and -CCR- boxes: initiation of minus strand synthesis by turnip yellow mosaic virus RNA-dependent RNA polymerase

A turnip yellow mosaic virus RNA-dependent RNA polymerase activity was used to study the template requirements for in vitro minus strand synthesis, which is initiated specifically opposite the 3'-CCA that terminates the 3'-tRNA-like structure. A deletion survey confirmed earlier results suggesting the absence of minus strand promoter elements upstream of the pseudoknotted acceptor stem and 3'-terminus. Reiteration of this 27-nt domain provided two competing initiation sites. By varying the added downstream element, it was shown that the pseudoknotted domain could be functionally replaced by various simple stem/loops, although with some decrease in activity. The addition of varying numbers of consecutive -CCA- triplets to the 3' end of the tRNA-like structure resulted in accurate initiation from each added triplet. A similar spectrum of initiations occurred with an unstructured RNA consisting of 12 consecutive -CCA- triplets and no additional viral sequence. Substitution mutations revealed no influence on minus strand synthesis of the identity of the nucleotide immediately upstream of a -CC- initiation site, but a preference for a purine immediately downstream. The introduction of secondary structure into the linear template showed that the usage of potential -CCR- initiation sites is influenced by nonspecific secondary structure. We conclude that specificity arises from the requirement that a -CCR- sequence be sterically accessible. This mechanism is only applicable to interactions that do not involve RNA unwinding during site selection, but may be used commonly in positive strand RNA virus replication and be applicable to other RNA-protein interactions.     

N 2-methylguanosine is iso-energetic with guanosine in RNA duplexes and GNRA tetraloops

Modified nucleotides are resource-intensive alternatives to the four nucleotides that constitute the bulk of natural RNAs. Yet, even in cases where modifications are highly conserved, their functions are difficult to identify. One possible function might be to modulate the stability of RNA structures. To investigate this possibility for N 2-methylguanosine (m2G), which is present in a wide variety of RNAs, we have determined the thermodynamic consequences of substituting m2G for G in G-C Watson-Crick pairs and G@U wobble pairs within RNA duplexes. The m2G substitution is iso-energetic with G in all cases, except for aninternal m2G@U pair, where it has a modest (0.3 kcal/mol) stabilizing effect. We have also examined theconsequences of replacing G by m2G, and A by N 6, N 6-dimethyladenosine (m26A) in the helix 45 tetraloop of 16S rRNA, which would otherwise be a standard GNRA tetraloop. This loop is a conserved, hypermethylated region of the ribosome where methylation appears to modulate activity. m26A substitution destabilizes the tetraloop, presumably because it prevents the formation of the G@A sheared pair it would otherwise contain. m2G substitution has no effect on tetraloop stability. Together, these results suggest that m2G is equally stable as either the s-cis or s-trans rotamer. The lack of a significant effect on secondary structural stability in these systems suggests that m2G is introduced into naturally occurring RNAs for reasons other than modulation of duplex stability.     

Thermodynamic analysis of an RNA combinatorial library contained in a short hairpin

Prediction of nucleic acid structure from sequence requires thermodynamic parameters for a variety of motifs, many of which are complex and consist of a large number of possible sequence combinations. Here we report an experimental approach for identifying the stable and unstable members of an RNA combinatorial library. Short model RNA hairpins consisting of 13 base pairs (bp) flanked by primer binding sites are constructed and separated according to their relative thermodynamic stabilities using temperature gradient gel electrophoresis (TGGE). Partially denaturing TGGE is carried out with potassium chloride, sodium chloride, or magnesium chloride salts in the gel. The TMs of model hairpins can be tuned by adjusting the concentration of urea in the gel while maintaining the correct order of stabilities for the hairpins. Mixtures of RNAs differing by a single Watson-Crick base pair are resolved according to their relative thermodynamic stabilities, as are mixtures of GC or AU base pair transversions differing in DeltaG degrees37 by only 0.3-0.5 kcal/mol. In addition, a simple combinatorial library with one position of randomization opposite a guanosine is prepared and separated into its four members by parallel and perpendicular TGGE. The order of thermodynamic stabilities for the library determined by TGGE is shown to be the same when assayed by UV-melting experiments. Analysis of the thermodynamics of folding of combinatorial libraries is general and may be applied to a wide variety of complex nucleic acid secondary and tertiary motifs in order to identify the stable and unstable members.     

Comparative structure analysis of vertebrate ribonuclease P RNA

Ribonuclease P cleaves 5'-precursor sequences from pre-tRNAs. All cellular RNase P holoenzymes contain homologous RNA elements; the eucaryal RNase P RNA, in contrast to the bacterial RNA, is catalytically inactive in the absence of the protein component(s). To understand the function of eucaryal RNase P RNA, knowledge of its structure is needed. Considerable effort has been devoted to comparative studies of the structure of this RNA from diverse organisms, including eucaryotes, primarily fungi, but also a limited set of vertebrates. The substantial differences in the sequences and structures of the vertebrate RNAs from those of other organisms have made it difficult to align the vertebrate sequences, thus limiting comparative studies. To expand our understanding of the structure of diverse RNase P RNAs, we have isolated by PCR and sequenced 13 partial RNase P RNA genes from 11 additional vertebrate taxa representing most extant major vertebrate lineages. Based on a recently proposed structure of the core elements of RNase P RNA, we aligned the sequences and propose a minimum consensus secondary structure for the vertebrate RNase P RNA.     

Ricin A-chain: kinetics, mechanism, and RNA stem-loop inhibitors

Ricin A-chain (RTA) catalyzes the depurination of a single adenine at position 4324 of 28S rRNA in a N-ribohydrolase reaction. The mechanism and specificity for RTA are examined using RNA stem-loop structures of 10-18 nucleotides which contain the required substrate motif, a GAGA tetraloop. At the optimal pH near 4.0, the preferred substrate is a 14-base stem-loop RNA which is hydrolyzed at 219 min-1 with a kcat/Km of 4.5 x 10(5) M-1 s-1 under conditions of steady-state catalysis. Smaller or larger stem-loop RNAs have lower kcat values, but all have Km values of approximately 5 microM. Both the 10- and 18-base substrates have kcat/Km near 10(4) M-1 s-1. Covalent cross-linking of the stem has a small effect on the kinetic parameters. Stem-loop DNA (10 bases) of the same sequence is also a substrate with a kcat/Km of 0.1 that for RNA. Chemical mechanisms for enzymatic RNA depurination reactions include leaving group activation, stabilization of a ribooxocarbenium transition state, a covalent enzyme-ribosyl intermediate, and ionization of the 2'-hydroxyl. A stem-loop RNA with p-nitrophenyl O-riboside at the depurination site is not a substrate, but binds tightly to the enzyme (Ki = 0.34 microM), consistent with a catalytic mechanism of leaving group activation. The substrate activity of stem-loop DNA eliminates ionization of the 2'-hydroxyl as a mechanism. Incorporation of the C-riboside formycin A at the depurination site provides an increased pKa of the adenine analogue at N7. Binding of this analogue (Ki = 9.4 microM) is weaker than substrate which indicates that the altered pKa at this position is not an important feature of transition state recognition. Stem-loop RNA with phenyliminoribitol at the depurination site increases the affinity substantially (Ki = 0.18 microM). The results are consistent with catalysis occurring by leaving group protonation at ring position(s) other than N7 leading to a ribooxocarbenium ion transition state. Small stem-loop RNAs have been identified with substrate activity within an order of magnitude of that reported for intact ribosomes.     

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.     

Automatic detection of conserved RNA structure elements in complete RNA virus genomes

We propose a new method for detecting conserved RNA secondary structures in a family of related RNA sequences. Our method is based on a combination of thermodynamic structure prediction and phylogenetic comparison. In contrast to purely phylogenetic methods, our algorithm can be used for small data sets of approximately 10 sequences, efficiently exploiting the information contained in the sequence variability. The procedure constructs a prediction only for those parts of sequences that are consistent with a single conserved structure. Our implementation produces reasonable consensus structures without user interference. As an example we have analysed the complete HIV-1 and hepatitis C virus (HCV) genomes as well as the small segment of hantavirus. Our method confirms the known structures in HIV-1 and predicts previously unknown conserved RNA secondary structures in HCV.     

M-Si-B合金非晶形成能力的CALPHAD模式评估

Evaluation of the amorphous-forming ability of M-Si-B ternary alloys using CALPHAD approach@長谷部光弘$日本九州工业大学 @OHTANI Hiroshi$Department of Materials Science and Engineering, Kyushu Institute of Technology @HASEBE Mitsuhiro$Department of Materials Science and Engineering, Kyushu Institute of Technology…