Sequence-dependent base-stacking stabilities guide tRNA folding energy landscapes

New substrates and determinants for tRNA recognition of RNA methyltransferase DNMT2/TRDMT1

Methylation is a common post-transcriptional modification of tRNAs, particularly in the anticodon loop region. The cytosine 38 (C38) in tRNAs, such as tRNAAsp-GUC, tRNAGly-GCC, tRNAVal-AAC, and tRNAGlu-CUC, can be methylated by human DNMT2/TRDMT1 and some homologs found in bacteria, plants, and animals. However, the substrate properties and recognition mechanism of DNMT2/TRDMT1 remain to be explored. Here, taking into consideration common features of the four known substrate tRNAs, we investigated methylation activities of DNMT2/TRDMT1 on the tRNAGly-GCC truncation and point mutants, and conformational changes of mutants. The results demonstrated that human DNMT2/TRDMT1 preferred substrate tRNAGly-GCC in vitro. L-shaped conformation of classical tRNA could be favourable for DNMT2/TRDMT1 activity. The complete sequence and structure of tRNA were dispensable for DNMT2/TRDMT1 activity, whereas T-arm was indispensable to this activity. G19, U20, and A21 in D-loop were identified as the important bases for DNMT2/TRDMT1 activity, while G53, C56, A58, and C61 in T-loop were found as the critical bases. The conserved CUXXCAC sequence in the anticodon loop was confirmed to be the most critical determinant, and it could stabilize C38-flipping to promote C38 methylation. Based on these tRNA properties, new substrates, tRNAVal-CAC and tRNAGln-CUG, were discovered in vitro. Moreover, a single nucleotide substitute, U32C, could convert non-substrate tRNAAla-AGC into a substrate for DNMT2/TRDMT1. Altogether, our findings imply that DNMT2/TRDMT1 relies on a delicate network involving both the primary sequence and tertiary structure of tRNA for substrate recognition.

Posttranscriptional modifications at the 37th position in the anticodon stem-loop of tRNA: structural insights from MD simulations

Transfer RNA (tRNA) is the most diversely modified RNA. Although the strictly conserved purine position 37 in the anticodon stem-loop undergoes modifications that are phylogenetically distributed, we do not yet fully understand the roles of these modifications. Therefore, molecular dynamics simulations are used to provide molecular-level details for how such modifications impact the structure and function of tRNA. A focus is placed on three hypermodified base families that include the parent i6A, t6A, and yW modifications, as well as derivatives. Our data reveal that the hypermodifications exhibit significant conformational flexibility in tRNA, which can be modulated by additional chemical functionalization. Although the overall structure of the tRNA anticodon stem remains intact regardless of the modification considered, the anticodon loop must rearrange to accommodate the bulky, dynamic hypermodifications, which includes changes in the nucleotide glycosidic and backbone conformations, and enhanced or completely new nucleobase-nucleobase interactions compared to unmodified tRNA or tRNA containing smaller (m1G) modifications at the 37th position. Importantly, the extent of the changes in the anticodon loop is influenced by the addition of small functional groups to parent modifications, implying each substituent can further fine-tune tRNA structure. Although the dominant conformation of the ASL is achieved in different ways for each modification, the molecular features of all modified tRNA drive the ASL domain to adopt the functional open-loop conformation. Importantly, the impact of the hypermodifications is preserved in different sequence contexts. These findings highlight the likely role of regulating mRNA structure and translation.

RNA modification enzyme TruB is a tRNA chaperone

Cellular RNAs are chemically modified by many RNA modification enzymes; however, often the functions of modifications remain unclear, such as for pseudouridine formation in the tRNA TΨC arm by the bacterial tRNA pseudouridine synthase TruB. Here we test the hypothesis that RNA modification enzymes also act as RNA chaperones. Using TruB as a model, we demonstrate that TruB folds tRNA independent of its catalytic activity, thus increasing the fraction of tRNA that can be aminoacylated. By rapid kinetic stopped-flow analysis, we identified the molecular mechanism of TruB's RNA chaperone activity: TruB binds and unfolds both misfolded and folded tRNAs thereby providing misfolded tRNAs a second chance at folding. Previously, it has been shown that a catalytically inactive TruB variant has no phenotype when expressed in an Escherichia coli truB KO strain [Gutgsell N, et al. (2000) RNA 6(12):1870-1881]. However, here we uncover that E. coli strains expressing a TruB variant impaired in tRNA binding and in in vitro tRNA folding cannot compete with WT E. coli. Consequently, the tRNA chaperone activity of TruB is critical for bacterial fitness. In conclusion, we prove the tRNA chaperone activity of the pseudouridine synthase TruB, reveal its molecular mechanism, and demonstrate its importance for cellular fitness. We discuss the likelihood that other RNA modification enzymes are also RNA chaperones.

A comprehensive collection of annotations to interpret sequence variation in human mitochondrial transfer RNAs

BACKGROUND

The abundance of biological data characterizing the genomics era is contributing to a comprehensive understanding of human mitochondrial genetics. Nevertheless, many aspects are still unclear, specifically about the variability of the 22 human mitochondrial transfer RNA (tRNA) genes and their involvement in diseases. The complex enrichment and isolation of tRNAs in vitro leads to an incomplete knowledge of their post-transcriptional modifications and three-dimensional folding, essential for correct tRNA functioning. An accurate annotation of mitochondrial tRNA variants would be definitely useful and appreciated by mitochondrial researchers and clinicians since the most of bioinformatics tools for variant annotation and prioritization available so far cannot shed light on the functional role of tRNA variations.

RESULTS

To this aim, we updated our MToolBox pipeline for mitochondrial DNA analysis of high throughput and Sanger sequencing data by integrating tRNA variant annotations in order to identify and characterize relevant variants not only in protein coding regions, but also in tRNA genes. The annotation step in the pipeline now provides detailed information for variants mapping onto the 22 mitochondrial tRNAs. For each mt-tRNA position along the entire genome, the relative tRNA numbering, tRNA type, cloverleaf secondary domains (loops and stems), mature nucleotide and interactions in the three-dimensional folding were reported. Moreover, pathogenicity predictions for tRNA and rRNA variants were retrieved from the literature and integrated within the annotations provided by MToolBox, both in the stand-alone version and web-based tool at the Mitochondrial Disease Sequence Data Resource (MSeqDR) website. All the information available in the annotation step of MToolBox were exploited to generate custom tracks which can be displayed in the GBrowse instance at MSeqDR website.

CONCLUSIONS

To the best of our knowledge, specific data regarding mitochondrial variants in tRNA genes were introduced for the first time in a tool for mitochondrial genome analysis, supporting the interpretation of genetic variants in specific genomic contexts.

Effects of heme binding on myoglobin folding: Coarse grained molecular simulations

Many proteins contain cofactors, such as heme, ATP and metal ions. Binding of cofactors is not only essential for their biological functions, but also can reshape the intrinsic energy landscape of protein molecules and modulate the folding and stability. However, the molecular mechanism of cofactor coupled protein folding is not well understood. In this work, we study the cofactor coupled folding of myoglobin, which is a typical cofactor (heme) containing protein, by performing molecular dynamics simulations with a structure-based protein model developed based on the energy landscape theory. We showed that the heme binding increases the stability of the myoglobin. More importantly, the heme binding tends to increase the protein folding cooperativity, and switch the folding process from a “three-state” mechanism to a “two-state” mechanism. We also showed that the folding pathways of the myoglobin can be modulated by the heme binding. By performing comparative simulations, we revealed that the above effects of heme binding are resulted from the heme induced folding of F-helix, which is otherwise unstructured at apo state, and the heme mediated contacting interactions around the heme binding site. The simulation results are consistent with available experimental data, and provide insights into the molecular mechanism of the effects of cofactor binding on the protein folding and stability.

Searching target sites on DNA by proteins: Role of DNA dynamics under confinement

DNA-binding proteins (DBPs) rapidly search and specifically bind to their target sites on genomic DNA in order to trigger many cellular regulatory processes. It has been suggested that the facilitation of search dynamics is achieved by combining 3D diffusion with one-dimensional sliding and hopping dynamics of interacting proteins. Although, recent studies have advanced the knowledge of molecular determinants that affect one-dimensional search efficiency, the role of DNA molecule is poorly understood. In this study, by using coarse-grained simulations, we propose that dynamics of DNA molecule and its degree of confinement due to cellular crowding concertedly regulate its groove geometry and modulate the inter-communication with DBPs. Under weak confinement, DNA dynamics promotes many short, rotation-decoupled sliding events interspersed by hopping dynamics. While this results in faster 1D diffusion, associated probability of missing targets by jumping over them increases. In contrast, strong confinement favours rotation-coupled sliding to locate targets but lacks structural flexibility to achieve desired specificity. By testing under physiological crowding, our study provides a plausible mechanism on how DNA molecule may help in maintaining an optimal balance between fast hopping and rotation-coupled sliding dynamics, to locate target sites rapidly and form specific complexes precisely.

MD Simulations of tRNA and Aminoacyl-tRNA Synthetases: Dynamics, Folding, Binding, and Allostery

While tRNA and aminoacyl-tRNA synthetases are classes of biomolecules that have been extensively studied for decades, the finer details of how they carry out their fundamental biological functions in protein synthesis remain a challenge. Recent molecular dynamics (MD) simulations are verifying experimental observations and providing new insight that cannot be addressed from experiments alone. Throughout the review, we briefly discuss important historical events to provide a context for how far the field has progressed over the past few decades. We then review the background of tRNA molecules, aminoacyl-tRNA synthetases, and current state of the art MD simulation techniques for those who may be unfamiliar with any of those fields. Recent MD simulations of tRNA dynamics and folding and of aminoacyl-tRNA synthetase dynamics and mechanistic characterizations are discussed. We highlight the recent successes and discuss how important questions can be addressed using current MD simulations techniques. We also outline several natural next steps for computational studies of AARS:tRNA complexes.

Search by proteins for their DNA target site: 1. The effect of DNA conformation on protein sliding

The recognition of DNA-binding proteins (DBPs) to their specific site often precedes by a search technique in which proteins slide, hop along the DNA contour or perform inter-segment transfer and 3D diffusion to dissociate and re-associate to distant DNA sites. In this study, we demonstrated that the strength and nature of the non-specific electrostatic interactions, which govern the search dynamics of DBPs, are strongly correlated with the conformation of the DNA. We tuned two structural parameters, namely curvature and the extent of helical twisting in circular DNA. These two factors are mutually independent of each other and can modulate the electrostatic potential through changing the geometry of the circular DNA conformation. The search dynamics for DBPs on circular DNA is therefore markedly different compared with linear B-DNA. Our results suggest that, for a given DBP, the rotation-coupled sliding dynamics is precluded in highly curved DNA (as well as for over-twisted DNA) because of the large electrostatic energy barrier between the inside and outside of the DNA molecule. Under such circumstances, proteins prefer to hop in order to explore interior DNA sites. The change in the balance between sliding and hopping propensities as a function of DNA curvature or twisting may result in different search efficiency and speed.

Topological constraints are major determinants of tRNA tertiary structure and dynamics and provide basis for tertiary folding cooperativity

Recent studies have shown that basic steric and connectivity constraints encoded at the secondary structure level are key determinants of 3D structure and dynamics in simple two-way RNA junctions. However, the role of these topological constraints in higher order RNA junctions remains poorly understood. Here, we use a specialized coarse-grained molecular dynamics model to directly probe the thermodynamic contributions of topological constraints in defining the 3D architecture and dynamics of transfer RNA (tRNA). Topological constraints alone restrict tRNA's allowed conformational space by over an order of magnitude and strongly discriminate against formation of non-native tertiary contacts, providing a sequence independent source of folding specificity. Topological constraints also give rise to long-range correlations between the relative orientation of tRNA's helices, which in turn provides a mechanism for encoding thermodynamic cooperativity between distinct tertiary interactions. These aspects of topological constraints make it such that only several tertiary interactions are needed to confine tRNA to its native global structure and specify functionally important 3D dynamics. We further show that topological constraints are conserved across tRNA's different naturally occurring secondary structures. Taken together, our results emphasize the central role of secondary-structure-encoded topological constraints in defining RNA 3D structure, dynamics and folding.