Efficient siRNA selection using hybridization thermodynamics

Small interfering RNA (siRNA) are widely used to infer gene function. Here, insights in the equilibrium of siRNA-target hybridization are used for selection of efficient siRNA. The accessibilities of siRNA and target mRNA for hybridization, as measured by folding free energy change, are shown to be significantly correlated with efficacy. For this study, a partition function calculation that considers all possible secondary structures is used to predict target site accessibility; a significant improvement over calculations that consider only the predicted lowest free energy structure or a set of low free energy structures. The predicted thermodynamic features, in addition to siRNA sequence features, are used as input for a support vector machine that selects functional siRNA. The method works well for predicting efficient siRNA (efficacy >70%) in a large siRNA data set from Novartis. The positive predictive value (percentage of sites predicted to be efficient for silencing that are) is as high as 87.6%. The sensitivity and specificity are 22.7 and 96.5%, respectively. When tested on data from different sources, the positive predictive value increased 8.1% by adding equilibrium terms to 25 local sequence features. Prediction of hybridization affinity using partition functions is now available in the RNAstructure software package.     

The activity of siRNA in mammalian cells is related to structural target accessibility: a comparison with antisense oligonucleotides

The biological activity of siRNA seems to be influenced by local characteristics of the target RNA, including local RNA folding. Here, we investigated quantitatively the relationship between local target accessibility and the extent of inhibition of the target gene by siRNA. Target accessibility was assessed by a computational approach that had been shown earlier to be consistent with experimental probing of target RNA. Two sites of ICAM-1 mRNA predicted to serve as accessible motifs and one site predicted to adopt an inaccessible structure were chosen to test siRNA constructs for suppression of ICAM-1 gene expression in ECV304 cells. The local target-dependent effectiveness of siRNA was compared with antisense oligonucleotides (asON). The concentration dependency of siRNA-mediated suppression indicates a >1000-fold difference between active siRNAs (IC₅₀ ≈ 0.2–0.5 nM) versus an inactive siRNA (IC₅₀ ≥ 1 µM) which is consistent with the activity pattern of asON when relating target suppression to predicted local target accessibility. The extremely high activity of the siRNA si2B (IC₅₀ = 0.24 nM) indicates that not all siRNAs shown to be active at the usual concentrations of >10–100 nM belong to this highly active species. The observations described here suggest an option to assess target accessibility for siRNA and, thus, support the design of active siRNA constructs. This approach can be automated, work at high throughput and is open to include additional parameters relevant to the biological activity of siRNA.     

Thermodynamic calculations and statistical correlations for oligo-probes design

Optimization of probe design for array-based experiments requires improved predictability of oligonucleotide hybridization behavior. Currently, designing oligonucleotides capable of interacting efficiently and specifically with the relevant target is not a routine procedure. Multiple examples demonstrate that oligonucleotides targeting different regions of the same RNA differ in their hybridization ability. The present work shows how thermodynamic evaluations of oligo-target duplex or oligo self-structure stabilities can facilitate probe design. Statistical analysis of large sets of hybridization data reveals that thermodynamic evaluation of oligonucleotide properties can be used to avoid poor RNA binders. Thermodynamic criteria for the selection of 20 and 21mers, which, with high probability, interact efficiently and specifically with their targets, are suggested. The design of longer oligonucleotides can also be facilitated by the same calculations of ΔG°_T values for oligo-target duplex or oligo self-structure stabilities and similar selection schemes.     

Mechanistic Approach to the Problem of Hybridization Efficiency in Fluorescent In Situ Hybridization

In fluorescent in situ hybridization (FISH), the efficiency of hybridization between the DNA probe and the rRNA has been related to the accessibility of the rRNA when ribosome content and cell permeability are not limiting. Published rRNA accessibility maps show that probe brightness is sensitive to the organism being hybridized and the exact location of the target site and, hence, it is highly unpredictable based on accessibility only. In this study, a model of FISH based on the thermodynamics of nucleic acid hybridization was developed. The model provides a mechanistic approach to calculate the affinity of the probe to the target site, which is defined as the overall Gibbs free energy change (ΔG°_overall) for a reaction scheme involving the DNA-rRNA and intramolecular DNA and rRNA interactions that take place during FISH. Probe data sets for the published accessibility maps and experiments targeting localized regions in the 16S rRNA of Escherichia coli were used to demonstrate that ΔG°_overall is a strong predictor of hybridization efficiency and superior to conventional estimates based on the dissociation temperature of the DNA/rRNA duplex. The use of the proposed model also allowed the development of mechanistic approaches to increase probe brightness, even in seemingly inaccessible regions of the 16S rRNA. Finally, a threshold ΔG°_overall of −13.0 kcal/mol was proposed as a goal in the design of FISH probes to maximize hybridization efficiency without compromising specificity.     

Thermodynamic characterization of naturally occurring RNA tetraloops

Although tetraloops are one of the most frequently occurring secondary structure motifs in RNA, less than one-third of the 30 most frequently occurring RNA tetraloops have been thermodynamically characterized. Therefore, 24 stem–loop sequences containing common tetraloops were optically melted, and the thermodynamic parameters ΔH°, ΔS°, ΔG°_37, and T_M for each stem–loop were determined. These new experimental values, on average, are 0.7 kcal/mol different from the values predicted for these tetraloops using the model proposed by Vecenie CJ, Morrow CV, Zyra A, Serra MJ. 2006. Biochemistry 45: 1400–1407. The data for the 24 tetraloops reported here were then combined with the data for 28 tetraloops that were published previously. A new model, independent of terminal mismatch data, was derived to predict the free energy contribution of previously unmeasured tetraloops. The average absolute difference between the measured values and the values predicted using this proposed model is 0.4 kcal/mol. This new experimental data and updated predictive model allow for more accurate calculations of the free energy of RNA stem–loops containing tetraloops and, furthermore, should allow for improved prediction of secondary structure from sequence. It was also shown that tetraloops within the sequence 5′-GCCNNNNGGC-3′ are, on average, 0.6 kcal/mol more stable than the same tetraloop within the sequence 5′-GGCNNNNGCC-3′. More systemic studies are required to determine the full extent of non-nearest-neighbor effects on tetraloop stability.     

Local RNA target structure influences siRNA efficacy: systematic analysis of intentionally designed binding regions

Contradictory reports in the literature have emphasised either the sequence of small interfering RNAs (siRNA) or the structure of their target molecules to be the major determinant of the efficiency of RNA interference (RNAi) approaches. In the present study, we analyse systematically the contributions of these parameters to siRNA activity by using deliberately designed mRNA constructs. The siRNA target sites were included in well-defined structural elements rendering them either highly accessible or completely involved in stable base-pairing. Furthermore, complementary sequence elements and various hairpins with different stem lengths and designs were used as target sites. Only one of the strands of the siRNA duplex was found to be capable of silencing via its respective target site, indicating that thermodynamic characteristics intrinsic to the siRNA strands are a basic determinant of siRNA activity. A significant obstruction of gene silencing by the same siRNA, however, was observed to be caused by structural features of the substrate RNA. Bioinformatic analysis of the mRNA structures suggests a direct correlation between the extent of gene-knockdown and the local free energy in the target region. Our findings indicate that, although a favourable siRNA sequence is a necessary prerequisite for efficient RNAi, complex target structures may limit the applicability even of carefully chosen siRNAs.     

siRNA作用效果的靶点依赖性

小分子双链RNA(siRNA)可以高效、特异地沉默目的基因表达 ,为基因功能研究及基因治疗提供了新工具。近年来研究表明针对基因mRNA不同位点随机设计的siRNA在作用效果上存在差异 ,siRNA作用效果有序列依赖性 ,而且与其在基因mRNA上的结合部位的高级结构有关 ,与反义核酸发挥作用的靶点依赖性类似。这一性质对设计高效siRNA为基因功能和基因治疗研究提供指导作用。     

Determination of the secondary structure of group II bulge loops using the fluorescent probe 2-aminopurine

Eleven RNA hairpins containing 2-aminopurine (2-AP) in either base-paired or single nucleotide bulge loop positions were optically melted in 1 M NaCl; and, the thermodynamic parameters ΔH°, ΔS°, ΔG°₃₇, and T_M for each hairpin were determined. Substitution of 2-AP for an A (adenosine) at a bulge position (where either the 2-AP or A is the bulge) in the stem of a hairpin, does not affect the stability of the hairpin. For group II bulge loops such as AA/U, where there is ambiguity as to which of the A residues is paired with the U, hairpins with 2-AP substituted for either the 5′ or 3′ position in the hairpin stem have similar stability. Fluorescent melts were performed to monitor the environment of the 2-AP. When the 2-AP was located distal to the hairpin loop on either the 5′ or 3′ side of the hairpin stem, the change in fluorescent intensity upon heating was indicative of an unpaired nucleotide. A database of phylogenetically determined RNA secondary structures was examined to explore the presence of naturally occurring bulge loops embedded within a hairpin stem. The distribution of bulge loops is discussed and related to the stability of hairpin structures.     

A study on the fundamental factors determining the efficacy of siRNAs with high C/G contents

Although there are many reports about the efficacy of siRNAs, it is not clear whether those siRNAs with high C/G contents can be used to silence their target mRNAs efficiently. In this study, we investigated the structure and function of a group of siRNAs with high C/G contents. The results showed that single siRNAs against the Calpain, Otoferlin and Her2 mRNAs could induce different silencing effects on their targets, suggesting that the accessibility to target sequences influences the efficacy of siRNA. Unexpectedly, a single siRNA could target its cognate sequence in the 3’UTR of EEF1D or the 5’UTR of hTRF2 or CDC6. Their interaction induced different modes of gene silencing. Furthermore, the introduction of mutations into the 3’ end of the passenger strand showed that the position and number of mutated nucleotides could exert some influence on the efficacy of siRNA. However, these mutations did not completely block the passenger strand from exerting its RNAi effect. Interestingly, our findings also indicated that the target mRNA might play essential roles in maintaining or discarding the guide strand in RISCs. Thus, the conclusion could be drawn that favorable siRNA sequences, accessible target structures and the fast cleavage mode are necessary and sufficient prerequisites for efficient RNAi.     

Adsorption Properties of Tetracycline onto Graphene Oxide: Equilibrium,Kinetic and Thermodynamic Studies

Graphene oxide (GO) nanoparticle is a high potential effective absorbent. Tetracycline (TC) is a broad-spectrum antibiotic produced, indicated for use against many bacterial infections. In the present research, a systematic study of the adsorption and release process of tetracycline on GO was performed by varying pH, sorption time and temperature. The results of our studies showed that tetracycline strongly loads on the GO surface via π–π interaction and cation–π bonding. Investigation of TC adsorption kinetics showed that the equilibrium was reached within 15 min following the pseudo-second-order model with observed rate constants of k₂ = 0.2742–0.5362 g/mg min (at different temperatures). The sorption data has interpreted by the Langmuir model with the maximum adsorption of 323 mg/g (298 K). The mean energy of adsorption was determined 1.83 kJ/mol (298 K) based on the Dubinin–Radushkevich (D–R) adsorption isotherm. Moreover, the thermodynamic parameters such as ΔH°, ΔS° and ΔG° values for the adsorption were estimated which indicated the endothermic and spontaneous nature of the sorption process. The electrochemistry approved an ideal reaction for the adsorption under electrodic process. Simulation of GO and TC was done by LAMMPS. Force studies in z direction showed that tetracycline comes close to GO sheet by C₈ direction. Then it goes far and turns and again comes close from amine group to the GO sheet.     

siRNA-mediated off-target gene silencing triggered by a 7 nt complementation

A growing body of evidence suggests that siRNA could generate off-target effects through different mechanisms. However, the full impact of off-target gene regulation on phenotypic induction and accordingly on data interpretation in the context of large-scale siRNA library screen has not been reported. Here we report on off-target gene silencing effects observed in a large-scale knockdown experiment designed to identify novel regulators of the HIF-1 pathway. All of the three ‘top hits’ from our screen have been demonstrated to result from off-target gene silencing. Two of the three ‘siRNA hits’ were found to directly trigger down-regulation of hif-1α mRNA through a 7 nt motif, AGGCAGT, that is present in both the hif-1α mRNA and the siRNAs. Further analysis revealed that the generation of off-target gene silencing via this 7 nt motif depends on the characteristics of the target mRNA, including the sequence context surrounding the complementary region, the position of the complementary region in the mRNA and the copy number of the complementary region. Interestingly, the off-target siRNA against hif-1α was also shown to trigger mRNA degradation with high probability of other genes that possess multiple copies of the AGGCAGT motif in the 3′-untranslated region. Lessons learned from this study will be a valuable asset to aid in designing siRNAs with more stringent target selectivity and improving ‘hits-follow-up’ strategies for future large-scale knockdown experiments.     

DNA multiplex hybridization on microarrays and thermodynamic stability in solution: a direct comparison

Hybridization intensities of 30 distinct short duplex DNAs measured on spotted microarrays, were directly compared with thermodynamic stabilities measured in solution. DNA sequences were designed to promote formation of perfect match, or hybrid duplexes containing tandem mismatches. Thermodynamic parameters ΔH°, ΔS° and ΔG° of melting transitions in solution were evaluated directly using differential scanning calorimetry. Quantitative comparison with results from 63 multiplex microarray hybridization experiments provided a linear relationship for perfect match and most mismatch duplexes. Examination of outliers suggests that both duplex length and relative position of tandem mismatches could be important factors contributing to observed deviations from linearity. A detailed comparison of measured thermodynamic parameters with those calculated using the nearest-neighbor model was performed. Analysis revealed the nearest-neighbor model generally predicts mismatch duplexes to be less stable than experimentally observed. Results also show the relative stability of a tandem mismatch is highly dependent on the identity of the flanking Watson–Crick (w/c) base pairs. Thus, specifying the stability contribution of a tandem mismatch requires consideration of the sequence identity of at least four base pair units (tandem mismatch and flanking w/c base pairs). These observations underscore the need for rigorous evaluation of thermodynamic parameters describing tandem mismatch stability.     

Correlation of exhaled breath temperature with bronchial blood flow in asthma

In asthma elevated rates of exhaled breath temperature changes (Δe°T) and bronchial blood flow (Q_aw) may be due to increased vascularity of the airway mucosa as a result of inflammation.We investigated the relationship of Δe°T with Q_aw and airway inflammation as assessed by exhaled nitric oxide (NO). We also studied the anti-inflammatory and vasoactive effects of inhaled corticosteroid and β₂-agonist.Δe°T was confirmed to be elevated (7.27 ± 0.6 Δ°C/s) in 19 asthmatic subjects (mean age ± SEM, 40 ± 6 yr; 6 male, FEV₁ 74 ± 6 % predicted) compared to 16 normal volunteers (4.23 ± 0.41 Δ°C/s, p < 0.01) (30 ± 2 yr) and was significantly increased after salbutamol inhalation in normal subjects (7.8 ± 0.6 Δ°C/ s, p < 0.05) but not in asthmatic patients. Q_aw, measured using an acetylene dilution method was also elevated in patients with asthma compared to normal subjects (49.47 ± 2.06 and 31.56 ± 1.6 μl/ml/min p < 0.01) and correlated with exhaled NO (r = 0.57, p < 0.05) and Δe°T (r = 0.525, p < 0.05). In asthma patients, Q_aw was reduced 30 minutes after the inhalation of budesonide 400 μg (21.0 ± 2.3 μl/ml/min, p < 0.05) but was not affected by salbutamol.Δe°T correlates with Q_aw and exhaled NO in asthmatic patients and therefore may reflect airway inflammation, as confirmed by the rapid response to steroids.     

Determination of base and backbone contributions to the thermodynamics of premelting and melting transitions in B DNA

In previous papers of this series the temperature-dependent Raman spectra of poly(dA)·poly(dT) and poly(dA–dT)·poly(dA–dT) were used to characterize structurally the melting and premelting transitions in DNAs containing consecutive A·T and alternating A·T/T·A base pairs. Here, we describe procedures for obtaining thermodynamic parameters from the Raman data. The method exploits base-specific and backbone-specific Raman markers to determine separate thermodynamic contributions of A, T and deoxyribosyl-phosphate moieties to premelting and melting transitions. Key findings include the following: (i) Both poly(dA)·poly(dT) and poly(dA–dT)· poly(dA–dT) exhibit robust premelting transitions, due predominantly to backbone conformational changes. (ii) The significant van’t Hoff premelting enthalpies of poly(dA)·poly(dT) [ΔH_vH^pm = 18.0 ± 1.6 kcal·mol^–1 (kilocalories per mole cooperative unit)] and poly(dA–dT)·poly(dA–dT) (ΔH_vH^pm = 13.4 ± 2.5 kcal·mol^–1) differ by an amount (~4.6 kcal·mol^–1) estimated as the contribution from three-centered inter-base hydrogen bonding in (dA)_n·(dT)_n tracts. (iii) The overall stacking free energy of poly(dA)· poly(dT) [–6.88 kcal·mol_bp^–1 (kilocalories per mole base pair)] is greater than that of poly(dA–dT)· poly(dA–dT) (–6.31 kcal·mol_bp^–1). (iv) The difference between stacking free energies of A and T is significant in poly(dA)·poly(dT) (ΔΔG_st = 0.8 ± 0.3 kcal· mol_bp^–1), but marginal in poly(dA–dT)·poly(dA–dT) (ΔΔG_st = 0.3 ± 0.3 kcal·mol_bp^–1). (v) In poly(dA)· poly(dT), the van’t Hoff parameters for melting of A (ΔH_vH^A = 407 ± 23 kcal·mol^–1, ΔS_vH^A = 1166 ± 67 cal·°K^–1·mol^–1, ΔG_vH(25°C)^A = 60.0 ± 3.2 kcal·mol^–1) are clearly distinguished from those of T (ΔH_vH^T = 185 ± 38 kcal·mol^–1, ΔS_vH^T = 516 ± 109 cal·°K^–1·mol^–1, ΔG_vH(25°C)^T = 27.1 ± 5.5 kcal·mol^–1). (vi) Similar relative differences are observed in poly(dA–dT)· poly(dA–dT) (ΔH_vH^A = 333 ± 54 kcal·mol^–1, ΔS_vH^A = 961 ± 157 cal·°K^–1·mol^–1, ΔG_vH(25°C)^A = 45.0 ± 7.6 kcal· mol^–1; ΔH_vH^T = 213 ± 30 kcal·mol^–1, ΔS_vH^T = 617 ± 86 cal·°K^–1·mol^–1, ΔG_vH(25°C)^T = 29.3 ± 4.9 kcal·mol^–1). The methodology employed here distinguishes thermodynamic contributions of base stacking, base pairing and backbone conformational ordering in the molecular mechanism of double-helical B DNA formation.     

A chemical synthesis of LNA-2,6-diaminopurine riboside, and the influence of 2'-O-methyl-2,6-diaminopurine and LNA-2,6-diaminopurine ribosides on the thermodynamic properties of 2'-O-methyl RNA/RNA heteroduplexes

Modified nucleotides are useful tools to study the structures, biological functions and chemical and thermodynamic stabilities of nucleic acids. Derivatives of 2,6-diaminopurine riboside (D) are one type of modified nucleotide. The presence of an additional amino group at position 2 relative to adenine results in formation of a third hydrogen bond when interacting with uridine. New method for chemical synthesis of protected 3′-O-phosphoramidite of LNA-2,6-diaminopurine riboside is described. The derivatives of 2′-O-methyl-2,6-diaminopurine and LNA-2,6-diaminopurine ribosides were used to prepare complete 2′-O-methyl RNA and LNA-2′-O-methyl RNA chimeric oligonucleotides to pair with RNA oligonucleotides. Thermodynamic stabilities of these duplexes demonstrated that replacement of a single internal 2′-O-methyladenosine with 2′-O-methyl-2,6-diaminopurine riboside (D^M) or LNA-2,6-diaminopurine riboside (D^L) increases the thermodynamic stability (ΔΔG°₃₇) on average by 0.9 and 2.3 kcal/mol, respectively. Moreover, the results fit a nearest neighbor model for predicting duplex stability at 37°C. D-A and D-G but not D-C mismatches formed by D^M or D^L generally destabilize 2′-O-methyl RNA/RNA and LNA-2′-O-methyl RNA/RNA duplexes relative to the same type of mismatches formed by 2′-O-methyladenosine and LNA-adenosine, respectively. The enhanced thermodynamic stability of fully complementary duplexes and decreased thermodynamic stability of some mismatched duplexes are useful for many RNA studies, including those involving microarrays.     

Evidence for unbenignant nature of glucose as a replacement for water in purple membranes

The net angle (θ_α) between the seven helical segments of the bacteriorhodopsin (bR) polypeptide and the normal to the membrane plane of the purple membrane (PM) is ~0° when determined by oriented far-ultraviolet (UV) circular dichroism (OCD) and midinfrared linear dichroism (IRLD). However, θ_α is ~11° when determined by high-resolution electron cryo-microscopy and electron diffraction (EMD). The spectral studies are made with fresh hydrated PM films at ambient temperature, whereas diffraction studies are made with aged glucose-embedded PM at -120 to -268°. The current study presents oriented far-UV OCD results of hydrated PM films embedded with glucose, which can best be interpreted as a change in the magnitude of θ_α (Δθ_α) from 0 to 23° as a consequence of glucose embedment. Possible alternative explanations contrary to this conclusion are discussed and ruled out. Therefore, it is suggested that a θ_α of ~11° as determined by the EMD method may not be an intrinsic structural characteristic of the native PM but an induced one. The differences in the Δθ_α value due to glucose embedment as determined by the two different approaches (23 vs. 11°) may be attributed to the drastic differences in the experimental conditions used, especially temperature. It is expected that at extremely low temperatures protein dynamics would be highly restricted and Δθ_α relatively curtailed. It is concluded that glucose may not be as benign to biological structures as has been assumed in the past.     

RNA Silencing Is Resistant to Low-Temperature in Grapevine

RNA silencing is a natural defence mechanism against viruses in plants, and transgenes expressing viral RNA-derived sequences were previously shown to confer silencing-based enhanced resistance against the cognate virus in several species. However, RNA silencing was shown to dysfunction at low temperatures in several species, questioning the relevance of this strategy in perennial plants such as grapevines, which are often exposed to low temperatures during the winter season. Here, we show that inverted-repeat (IR) constructs trigger a highly efficient silencing reaction in all somatic tissues in grapevines. Similarly to other plant species, IR-derived siRNAs trigger production of secondary transitive siRNAs. However, and in sharp contrast to other species tested to date where RNA silencing is hindered at low temperature, this process remained active in grapevine cultivated at 4°C. Consistently, siRNA levels remained steady in grapevines cultivated between 26°C and 4°C, whereas they are severely decreased in Arabidopsis grown at 15°C and almost undetectable at 4°C. Altogether, these results demonstrate that RNA silencing operates in grapevine in a conserved manner but is resistant to far lower temperatures than ever described in other species.     

Designer proteins that competitively inhibit Gαq by targeting its effector site

During signal transduction, the G protein, Gα_q, binds and activates phospholipase C-β isozymes. Several diseases have been shown to manifest upon constitutively activating mutation of Gα_q, such as uveal melanoma. Therefore, methods are needed to directly inhibit Gα_q. Previously, we demonstrated that a peptide derived from a helix-turn-helix (HTH) region of PLC-β3 (residues 852–878) binds Gα_q with low micromolar affinity and inhibits Gα_q by competing with full-length PLC-β isozymes for binding. Since the HTH peptide is unstructured in the absence of Gα_q, we hypothesized that embedding the HTH in a folded protein might stabilize the binding-competent conformation and further improve the potency of inhibition. Using the molecular modeling software Rosetta, we searched the Protein Data Bank for proteins with similar HTH structures near their surface. The candidate proteins were computationally docked against Gα_q, and their surfaces were redesigned to stabilize this interaction. We then used yeast surface display to affinity mature the designs. The most potent design bound Gα_q/i with high affinity in vitro (K_D = 18 nM) and inhibited activation of PLC-β isozymes in HEK293 cells. We anticipate that our genetically encoded inhibitor will help interrogate the role of Gα_q in healthy and disease model systems. Our work demonstrates that grafting interaction motifs into folded proteins is a powerful approach for generating inhibitors of protein–protein interactions.