Predicting structure and stability for RNA complexes with intermolecular loop–loop base-pairing

RNA loop–loop interactions are essential for genomic RNA dimerization and regulation of gene expression. In this article, a statistical mechanics-based computational method that predicts the structures and thermodynamic stabilities of RNA complexes with loop–loop kissing interactions is described. The method accounts for the entropy changes for the formation of loop–loop interactions, which is a notable advancement that other computational models have neglected. Benchmark tests with several experimentally validated systems show that the inclusion of the entropy parameters can indeed improve predictions for RNA complexes. Furthermore, the method can predict not only the native structures of RNA/RNA complexes but also alternative metastable structures. For instance, the model predicts that the SL1 domain of HIV-1 RNA can form two different dimer structures with similar stabilities. The prediction is consistent with experimental observation. In addition, the model predicts two different binding sites for hTR dimerization: One binding site has been experimentally proposed, and the other structure, which has a higher stability, is structurally feasible and needs further experimental validation.     

What’s in a loop?

DNAs and proteins are major classes of biomolecules that differ in many aspects. However, a considerable number of their members also share a common architectural feature that enables the assembly of multi-protein complexes and thereby permits the effective processing of signals: loop structures of substantial sizes. Here we briefly review a few representative examples and suggest a functional classification of different types of loop structures. In proteins, these loops occur in protein regions classified as intrinsically disordered. Studying such loops, their binders and their interactions with other loops should reveal much about cellular information computation and signaling network architectures. It is also expected to provide critical information for synthetic biologists and bioengineers.     

A feedforward loop links Gaucher and Parkinson's diseases?

Mutations in the GBA gene that encodes glucocerebrosidase cause the lysosomal storage disorder Gaucher disease but also increase the risk for Parkinson's disease. Mazzulli et?al. (2011) uncover a possible mechanism to explain this connection: loss of glucocerebrosidase creates a positive feedback loop of reduced lysosomal function and α-synuclein accumulation, ultimately leading to neurodegeneration.     

普氏野马线粒体DNAD—loop区序列多态性研究

目的：了解中国新疆吉木萨尔野马繁殖中心的普氏野马（Equus przewalskii）遗传多样性及其遗传背景。方法：采用PCR产物直接测序法,对15匹普氏野马线粒体DNAD—loop高变区进行测序分析。结果：测定15个个体的线粒体DNAD—loop高变区15464—15866片段序列402bp。检测到12种单倍型,包括37个多态位点,占全部序列的9．2％,其中转换位点2,4个、颠换位点20个、转换位点和颠换并存位点8个、缺失位点3个。A％＋T％含量（56．1％）高于G％＋C％含量（43．9％）,平均A含量为28．4％,T含量为27．7％,C含量为29％,G含为14．9％。单倍型间平均遗传距离为0．030,单倍型多态性（h）为1±0．00116,核苷酸多态性（7c）为2．90％。15匹普氏野马线粒体DNAD—loop高变区之间平均核苷酸变异率为2．48％。结论：研究表明我国新疆吉木萨尔野马繁殖中心的普氏野马线粒体DNAD—loop区序列存在丰富的多态性。     

Predicting loop–helix tertiary structural contacts in RNA pseudoknots

Tertiary interactions between loops and helical stems play critical roles in the biological function of many RNA pseudoknots. However, quantitative predictions for RNA tertiary interactions remain elusive. Here we report a statistical mechanical model for the prediction of noncanonical loop–stem base-pairing interactions in RNA pseudoknots. Central to the model is the evaluation of the conformational entropy for the pseudoknotted folds with defined loop–stem tertiary structural contacts. We develop an RNA virtual bond-based conformational model (Vfold model), which permits a rigorous computation of the conformational entropy for a given fold that contains loop–stem tertiary contacts. With the entropy parameters predicted from the Vfold model and the energy parameters for the tertiary contacts as inserted parameters, we can then predict the RNA folding thermodynamics, from which we can extract the tertiary contact thermodynamic parameters from theory–experimental comparisons. These comparisons reveal a contact enthalpy (ΔH) of −14 kcal/mol and a contact entropy (ΔS) of −38 cal/mol/K for a protonated C⁺•(G–C) base triple at pH 7.0, and (ΔH = −7 kcal/mol, ΔS = −19 cal/mol/K) for an unprotonated base triple. Tests of the model for a series of pseudoknots show good theory–experiment agreement. Based on the extracted energy parameters for the tertiary structural contacts, the model enables predictions for the structure, stability, and folding pathways for RNA pseudoknots with known or postulated loop–stem tertiary contacts from the nucleotide sequence alone.     

Does conformational free energy distinguish loop conformations in proteins?

Limitations in protein homology modeling often arise from the inability to adequately model loops. In this paper we focus on the selection of loop conformations. We present a complete computational treatment that allows the screening of loop conformations to identify those that best fit a molecular model. The stability of a loop in a protein is evaluated via computations of conformational free energies in solution, i.e., the free energy difference between the reference structure and the modeled one. A thermodynamic cycle is used for calculation of the conformational free energy, in which the total free energy of the reference state (i.e., gas phase) is the CHARMm potential energy. The electrostatic contribution of the solvation free energy is obtained from solving the finite-difference Poisson-Boltzmann equation. The nonpolar contribution is based on a surface area-based expression. We applied this computational scheme to a simple but well-characterized system, the antibody hypervariable loop (complementarity-determining region, CDR). Instead of creating loop conformations, we generated a database of loops extracted from high-resolution crystal structures of proteins, which display geometrical similarities with antibody CDRs. We inserted loops from our database into a framework of an antibody; then we calculated the conformational free energies of each loop. Results show that we successfully identified loops with a "reference-like" CDR geometry, with the lowest conformational free energy in gas phase only. Surprisingly, the solvation energy term plays a confusing role, sometimes discriminating "reference-like" CDR geometry and many times allowing "non-reference-like" conformations to have the lowest conformational free energies (for short loops). Most "reference-like" loop conformations are separated from others by a gap in the gas phase conformational free energy scale. Naturally, loops from antibody molecules are found to be the best models for long CDRs (> or = 6 residues), mainly because of a better packing of backbone atoms into the framework of the antibody model.     

mRNA secondary structure optimization using a correlated stem–loop prediction

Secondary structure of messenger RNA plays an important role in the bio-synthesis of proteins. Its negative impact on translation can reduce the yield of protein by slowing or blocking the initiation and movement of ribosomes along the mRNA, becoming a major factor in the regulation of gene expression. Several algorithms can predict the formation of secondary structures by calculating the minimum free energy of RNA sequences, or perform the inverse process of obtaining an RNA sequence for a given structure. However, there is still no approach to redesign an mRNA to achieve minimal secondary structure without affecting the amino acid sequence. Here we present the first strategy to optimize mRNA secondary structures, to increase (or decrease) the minimum free energy of a nucleotide sequence, without changing its resulting polypeptide, in a time-efficient manner, through a simplistic approximation to hairpin formation. Our data show that this approach can efficiently increase the minimum free energy by >40%, strongly reducing the strength of secondary structures. Applications of this technique range from multi-objective optimization of genes by controlling minimum free energy together with CAI and other gene expression variables, to optimization of secondary structures at the genomic level.     

Advantages of dynamic “closed loop” stable isotope flux phenotyping over static “open loop” clamps in detecting silent genetic and dietary phenotypes

In vivo insulin sensitivity can be assessed using “open loop” clamp or “closed loop” methods. Open loop clamp methods are static, and fix plasma glucose independently from plasma insulin. Closed loop methods are dynamic, and assess glucose disposal in response to a stable isotope labeled glucose tolerance test. Using PPARα^−/− mice, open and closed loop assessments of insulin sensitivity/glucose disposal were compared. Indirect calorimetry done for the assessment of diurnal substrate utilization/metabolic flexibility showed that chow fed PPARα^−/− mice had increased glucose utilization during the light (starved) cycle. Euglycemic clamps showed no differences in insulin stimulated glucose disposal, whether for chow or high fat diets, but did show differences in basal glucose clearance for chow fed PPARα^−/− versus SV129J-wt mice. In contrast, the dynamic stable isotope labeled glucose tolerance tests reveal enhanced glucose disposal for PPARα^−/− versus SV129J-wt, for chow and high fat diets. Area under the curve for plasma labeled and unlabeled glucose for PPARα^−/− was ≈1.7-fold lower, P < 0.01 during the stable isotope labeled glucose tolerance test for both diets. Area under the curve for plasma insulin was 5-fold less for the chow fed SV129J-wt (P < 0.01) but showed no difference on a high fat diet (0.30 ± 0.1 for SV129J-wt vs. 0.13 ± 0.10 for PPARα^−/−, P = 0.28). This study demonstrates that dynamic stable isotope labeled glucose tolerance test can assess “silent” metabolic phenotypes, not detectable by the static, “open loop”, euglycemic or hyperglycemic clamps. Both open loop and closed loop methods may describe different aspects of metabolic inflexibility and insulin sensitivity.     

Real-time quantification of microRNAs by stem–loop RT–PCR

A novel microRNA (miRNA) quantification method has been developed using stem–loop RT followed by TaqMan PCR analysis. Stem–loop RT primers are better than conventional ones in terms of RT efficiency and specificity. TaqMan miRNA assays are specific for mature miRNAs and discriminate among related miRNAs that differ by as little as one nucleotide. Furthermore, they are not affected by genomic DNA contamination. Precise quantification is achieved routinely with as little as 25 pg of total RNA for most miRNAs. In fact, the high sensitivity, specificity and precision of this method allows for direct analysis of a single cell without nucleic acid purification. Like standard TaqMan gene expression assays, TaqMan miRNA assays exhibit a dynamic range of seven orders of magnitude. Quantification of five miRNAs in seven mouse tissues showed variation from less than 10 to more than 30000 copies per cell. This method enables fast, accurate and sensitive miRNA expression profiling and can identify and monitor potential biomarkers specific to tissues or diseases. Stem–loop RT–PCR can be used for the quantification of other small RNA molecules such as short interfering RNAs (siRNAs). Furthermore, the concept of stem–loop RT primer design could be applied in small RNA cloning and multiplex assays for better specificity and efficiency.     

loop序列对双分子三螺旋DNA稳定性的影响

运用分子力学方法研究了 2 0个双分子三螺旋DNA ,其嘧啶链序列是 5′- dTTCTTTTC- L1TTTL5 -CTTTTCTT -3′,划线的 5个核苷酸组成loop环 ,L1和L5可以是任意的核苷酸 ;嘌呤链的序列是 5′- GAAAAGAA 3′和 5′- AAGAAAAG -3′ ;2条链方向相反 .对 2 0个不同loop序列双分子三螺旋DNA稳定性的研究结果表明 ,5′- loop三螺旋比相应的 3′- loop三螺旋更稳定 ,嘌呤比嘧啶与相邻碱基的堆积作用大 .2 0个三螺旋DNA的相对稳定性主要是由loop环上L1和L5碱基组成不同决定的     

ILOOP – a web application for two-channel microarray interwoven loop design

Microarray technology is widely applied to address complex scientific questions. However, there remain fundamental issues on how to design experiments to ensure that the resulting data enables robust statistical analysis. Interwoven loop design has several advantages over other designs. However it suffers in the complexity of design. We have implemented an online web application which allows users to find optimal loop designs for two-color microarray experiments. Given a number of conditions (such as treatments or time points) and replicates, the application will find the best possible design of the experiment and output experimental parameters. It is freely available from http://mcbc.usm.edu/iloop.     

CDR-H3 loop ensemble in solution – conformational selection upon antibody binding

ABSTRACT

We analyzed pairs of protein-binding, peptide-binding and hapten-binding antibodies crystallized as complex and in the absence of the antigen with and without conformational differences upon binding in the complementarity-determining region (CDR)-H3 loop. Here, we introduce a molecular dynamics-based approach to capture a diverse conformational ensemble of the CDR-H3 loop in solution. The results clearly indicate that the inherently flexible CDR-H3 loop indeed needs to be characterized as a conformational ensemble. The conformational changes of the CDR-H3 loop in all antibodies investigated follow the paradigm of conformation selection, because we observe the experimentally determined binding competent conformation without the presence of the antigen within the ensemble of pre-existing conformational states in solution before binding. We also demonstrate for several examples that the conformation observed in the antibody crystal structure without antigen present is actually selected to bind the carboxyterminal tail region of the antigen-binding fragment (Fab). Thus, special care must be taken when characterizing antibody CDR-H3 loops by Fab X-ray structures, and the possibility that pre-existing conformations are present should always be considered.     

A c-Myc–SIRT1 feedback loop regulates cell growth and transformation

The protein deacetylase SIRT1 has been implicated in a variety of cellular functions, including development, cellular stress responses, and metabolism. Increasing evidence suggests that similar to its counterpart, Sir2, in yeast, Caenorhabditis elegans, and Drosophila melanogaster, SIRT1 may function to regulate life span in mammals. However, SIRT1''s role in cancer is unclear. During our investigation of SIRT1, we found that c-Myc binds to the SIRT1 promoter and induces SIRT1 expression. However, SIRT1 interacts with and deacetylates c-Myc, resulting in decreased c-Myc stability. As a consequence, c-Myc''s transformational capability is compromised in the presence of SIRT1. Overall, our experiments identify a c-Myc–SIRT1 feedback loop in the regulation of c-Myc activity and cellular transformation, supporting/suggesting a role of SIRT1 in tumor suppression.     

A stem–loop structure directs oskar mRNA to microtubule minus ends

mRNA transport coupled with translational control underlies the intracellular localization of many proteins in eukaryotic cells. This is exemplified in Drosophila, where oskar mRNA transport and translation at the posterior pole of the oocyte direct posterior patterning of the embryo. oskar localization is a multistep process. Within the oocyte, a spliced oskar localization element (SOLE) targets oskar mRNA for plus end-directed transport by kinesin-1 to the posterior pole. However, the signals mediating the initial minus end-directed, dynein-dependent transport of the mRNA from nurse cells into the oocyte have remained unknown. Here, we show that a 67-nt stem–loop in the oskar 3′ UTR promotes oskar mRNA delivery to the developing oocyte and that it shares functional features with the fs(1)K10 oocyte localization signal. Thus, two independent cis-acting signals, the oocyte entry signal (OES) and the SOLE, mediate sequential dynein- and kinesin-dependent phases of oskar mRNA transport during oogenesis. The OES also promotes apical localization of injected RNAs in blastoderm stage embryos, another dynein-mediated process. Similarly, when ectopically expressed in polarized cells of the follicular epithelium or salivary glands, reporter RNAs bearing the oskar OES are apically enriched, demonstrating that this element promotes mRNA localization independently of cell type. Our work sheds new light on how oskar mRNA is trafficked during oogenesis and the RNA features that mediate minus end-directed transport.     

Favorable 2′-substitution in the loop region of a thrombin-binding DNA aptamer

Simple 2′-OMe-chemical modification in the loop region of the 15mer G-rich DNA sequence GGTTGGTGTGGTTGG is reported. The G-quadruplex structure of this thrombin-binding aptamer (TBA), is stabilized by single modifications (T?→?2′-OMe-U), depending on the position of the modification. The structural stability also renders significantly increased inhibition of thrombin-induced fibrin polymerization, a process closely associated with blood-clotting.     

Is the V3 loop involved in HIV binding to CD4?

The entry of the human immunodeficiency virus into cells requires the interaction of the viral envelope glycoprotein gp120 with CD4 and a chemokine receptor. The gp120 binding site has been previously mapped to the Ig-CDR2-like region of CD4 first domain. A second area of this domain (Ig-CDR3-like region) is involved in gp120-CD4 interactions, but its gp120 counterpart remained so far unknown. Using a photoaffinity labeling experiment, we demonstrate that a peptide, mapping the (307-330)m region of HIV-MN-gp120 V3 loop, binds a sequence including a part of the Ig-CDR3-like region. These results may contribute to explain the complex mechanism of human immunodeficiency virus penetration, helping the development of new therapeutic agents.     

Loop–loop interaction in an adenine-sensing riboswitch: A molecular dynamics study

Riboswitches are mRNA-based molecules capable of controlling the expression of genes. They undergo conformational changes upon ligand binding, and as a result, they inhibit or promote the expression of the associated gene. The close connection between structural rearrangement and function makes a detailed knowledge of the molecular interactions an important step to understand the riboswitch mechanism and efficiency. We have performed all-atom molecular dynamics simulations of the adenine-sensing add A-riboswitch to study the breaking of the kissing loop, one key tertiary element in the aptamer structure. We investigated the aptamer domain of the add A-riboswitch in complex with its cognate ligand and in the absence of the ligand. The opening of the hairpins was simulated using umbrella sampling using the distance between two loops as the reaction coordinate. A two-step process was observed in all the simulated systems. First, a general loss of stacking and hydrogen bond interactions is seen. The last interactions that break are the two base pairs G37-C61 and G38-C60, but the break does not affect the energy profile, indicating their pivotal role in the tertiary structure formation but not in the structure stabilization. The junction area is partially organized before the kissing loop formation and residue A24 anchors together the loop helices. Moreover, when the distance between the loops is increased, one of the hairpins showed more flexibility by changing its orientation in the structure, while the other conserved its coaxial arrangement with the rest of the structure.