Closing loop base pairs in RNA loop–loop complexes: structural behavior, interaction energy and solvation analysis through molecular dynamics simulations

Nanosecond molecular dynamics using the Ewald summation method have been performed to elucidate the structural and energetic role of the closing base pair in loop–loop RNA duplexes neutralized by Mg²⁺ counterions in aqueous phases. Mismatches GA, CU and Watson–Crick GC base pairs have been considered for closing the loop of an RNA in complementary interaction with HIV-1 TAR. The simulations reveal that the mismatch GA base, mediated by a water molecule, leads to a complex that presents the best compromise between flexibility and energetic contributions. The mismatch CU base pair, in spite of the presence of an inserted water molecule, is too short to achieve a tight interaction at the closing-loop junction and seems to force TAR to reorganize upon binding. An energetic analysis has allowed us to quantify the strength of the interactions of the closing and the loop–loop pairs throughout the simulations. Although the water-mediated GA closing base pair presents an interaction energy similar to that found on fully geometry-optimized structure, the water-mediated CU closing base pair energy interaction reaches less than half the optimal value.     

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.     

Gaucher disease glucocerebrosidase and α-synuclein form a bidirectional pathogenic loop in synucleinopathies

Parkinson's disease (PD), an adult neurodegenerative disorder, has been clinically linked to the lysosomal storage disorder Gaucher disease (GD), but the mechanistic connection is not known. Here, we show that functional loss of GD-linked glucocerebrosidase (GCase) in primary cultures or human iPS neurons compromises lysosomal protein degradation, causes accumulation of α-synuclein (α-syn), and results in neurotoxicity through aggregation-dependent mechanisms. Glucosylceramide (GlcCer), the GCase substrate, directly influenced amyloid formation of purified α-syn by stabilizing soluble oligomeric intermediates. We further demonstrate that α-syn inhibits the lysosomal activity of normal GCase in neurons and idiopathic PD brain, suggesting that GCase depletion contributes to the pathogenesis of sporadic synucleinopathies. These findings suggest that the bidirectional effect of α-syn and GCase forms a positive feedback loop that may lead to a self-propagating disease. Therefore, improved targeting of GCase to lysosomes may represent a specific therapeutic approach for PD and other synucleinopathies.     

An autocrine sphingosine-1-phosphate signaling loop enhances NF-κB-activation and survival

Background  

Sphingosine-1-phosphate (S1P) is a bioactive lipid that regulates a multitude of cellular functions, including cell proliferation, survival, migration and angiogenesis. S1P mediates its effects either by signaling through G protein-coupled receptors (GPCRs) or through an intracellular mode of action. In this study, we have investigated the mechanism behind S1P-induced survival signalling.