3D feasibility of 2D RNA-RNA interaction paths by stepwise folding simulations.

The structure of an RNA, and even more so its interactions with other RNAs, provide valuable information about its function. Secondary structure-based tools for RNA-RNA interaction predictions provide a quick way to identify possible interaction targets and structures. However, these tools ignore the effect of steric hindrance on the tertiary (3D) structure level, and do not consider whether a suitable folding pathway exists to form the interaction. As a consequence, these tools often predict interactions that are unrealistically long and could be formed (in three dimensions) only by going through highly entangled intermediates. Here, we present a computational pipeline to assess whether a proposed secondary (2D) structure interaction is sterically feasible and reachable along a plausible folding pathway. To this end, we simulate the folding of a series of 3D structures along a given 2D folding path. To avoid the complexity of large-scale atomic resolution simulations, our pipeline uses coarse-grained 3D modeling and breaks up the folding path into small steps, each corresponding to the extension of the interaction by 1 or 2 bp. We apply our pipeline to analyze RNA-RNA interaction formation for three selected RNA-RNA complexes. We find that kissing hairpins, in contrast to interactions in the exterior loop, are difficult to extend and tend to get stuck at an interaction length of 6 bp. Our tool, including source code, documentation, and sample data, is available at www.github.com/irenekb/RRI-3D.

Local RNA folding revisited.

Most of the functional RNA elements located within large transcripts are local. Local folding therefore serves a practically useful approximation to global structure prediction. Due to the sensitivity of RNA secondary structure prediction to the exact definition of sequence ends, accuracy can be increased by averaging local structure predictions over multiple, overlapping sequence windows. These averages can be computed efficiently by dynamic programming. Here we revisit the local folding problem, present a concise mathematical formalization that generalizes previous approaches and show that correct Boltzmann samples can be obtained by local stochastic backtracing in McCaskill's algorithms but not from local folding recursions. Corresponding new features are implemented in the ViennaRNA package to improve the support of local folding. Applications include the computation of maximum expected accuracy structures from RNAplfold data and a mutual information measure to quantify the sensitivity of individual sequence positions.

Insights into the secondary and tertiary structure of the Bovine Viral Diarrhea Virus Internal Ribosome Entry Site.

The internal ribosome entry site (IRES) RNA of bovine viral diarrhoea virus (BVDV), an economically significant Pestivirus, is required for the cap-independent translation of viral genomic RNA. Thus, it is essential for viral replication and pathogenesis. We applied a combination of high-throughput biochemical RNA structure probing (SHAPE-MaP) and in silico modelling approaches to gain insight into the secondary and tertiary structures of BVDV IRES RNA. Our study demonstrated that BVDV IRES RNA in solution forms a modular architecture composed of three distinct structural domains (I-III). Two regions within domain III are represented in tertiary interactions to form an H-type pseudoknot. Computational modelling of the pseudoknot motif provided a fine-grained picture of the tertiary structure and local arrangement of helices in the BVDV IRES. Furthermore, comparative genomics and consensus structure predictions revealed that the pseudoknot is evolutionarily conserved among many Pestivirus species. These studies provide detailed insight into the structural arrangement of BVDV IRES RNA H-type pseudoknot and encompassing motifs that likely contribute to the optimal functionality of viral cap-independent translation element.

3D based on 2D: Calculating helix angles and stacking patterns using forgi 2.0, an RNA Python library centered on secondary structure elements.

We present forgi, a Python library to analyze the tertiary structure of RNA secondary structure elements. Our representation of an RNA molecule is centered on secondary structure elements (stems, bulges and loops). By fitting a cylinder to the helix axis, these elements are carried over into a coarse-grained 3D structure representation. Integration with Biopython allows for handling of all-atom 3D information. forgi can deal with a variety of file formats including dotbracket strings, PDB and MMCIF files. We can handle modified residues, missing residues, cofold and multifold structures as well as nucleotide numbers starting at arbitrary positions. We apply this library to the study of stacking helices in junctions and pseudo knots and investigate how far stacking helices in solved experimental structures can divert from coaxial geometries.

Long term follow up of patients after allogeneic stem cell transplantation and transfusion of HSV-TK transduced T-cells.

Allogeneic stem cell transplantation (allo-HSCT) is one of the curative treatments for hematologic malignancies, but is hampered by severe complications, such as acute or chronic graft-versus-host-disease (aGvHD; cGvHD) and infections. CD34-selection of stem cells reduces the risk of aGvHD, but also leads to increased infectious complications and relapse. Thus, we studied the safety, efficacy, and feasibility of transfer of gene modified donor T-cells shortly after allo-HSCT in two clinical trials between 2002 and 2007 and here we compare the results to unmodified donor leukocyte infusion (DLI). The aim of these trials was to provide patients with the protection of T-cells after T-cell-depleted allo-HSCT in the matched or mismatched donor setting with an option to delete transduced T-cells, if severe aGvHD occurred within the trial period. Donor-T-cells were transduced with the replication-deficient retrovirus SFCMM-3, expressing HSV-TK and the truncated ΔLNGFR for selection of transduced cells. Transduced cells were transfused either after day +60 (matched donors) or on day +42 (haploidentical donors). Nine patients were included in the first trial (MHH; 2002 until 2007), two were included in TK007 (2005-2009) and six serves as a control group for outcome after haploidentical transplantation without HSV-TK-transduced DLI. Three patients developed acute GvHD, two had grade I of the skin, one had aGvHD on day +131 (post-HSCT; +89 post-HSV-TK DLI) grade II, which was successfully controlled by ganciclovir (GCV). Donor chimerism was stabilized after transfusion of the transduced cells in all patients treated. Functionality of HSV-TK gene expressing T-cells was shown by loss of bcr-able gene expression as well as by control of cytomegalovirus-reactivation. To date, six patients have relapsed and died, two after a second hematopoietic stem cell transplantation without T-cell depletion or administration of unmodified T-cells. Eleven patients (seven post-HSV-TK DLI) are alive and well to date.