Protein folding and tRNA biology

Polypeptides can fold into tertiary structures while they are synthesized by the ribosome. In addition to the amino acid sequence, protein folding is determined by several factors within the cell. Among others, the folding pathway of a nascent polypeptide can be affected by transient interactions with other proteins, ligands, or the ribosome, as well as by the translocation through membrane pores. Particularly, the translation machinery and the population of tRNA under different physiological or adaptive responses can dramatically affect protein folding. This review summarizes the scientific evidence describing the role of translation kinetics and tRNA populations on protein folding and addresses current efforts to better understand tRNA biology. It is organized into three main parts, which are focused on: (i) protein folding in the cellular context; (ii) tRNA biology and the complexity of the tRNA population; and (iii) available methods and technical challenges in the characterization of tRNA pools. In this manner, this work illustrates the ways by which functional properties of proteins may be modulated by cellular tRNA populations.     

Mg2+-induced tRNA folding

Mg(2+)-induced folding of yeast tRNA(Phe) was examined at low ionic strength in steady-state and kinetic experiments. By using fluorescent labels attached to tRNA, four conformational transitions were revealed when the Mg(2+) concentration was gradually increased. The last two transitions were not accompanied by changes in the number of base pairs. The observed transitions were attributed to Mg(2+) binding to four distinct types of sites. The first two types are strong sites with K(diss) of 4 and 16 microM. The sites of the third and fourth types are weak with a K(diss) of 2 and 20 mM. Accordingly, the Mg(2+)-binding sites previously classified as "strong" and "weak" can be further subdivided into two subtypes each. Fluorescent transition I is likely to correspond to Mg(2+) binding to a unique strong site selective for Mg(2+); binding to this site causes only minor A(260) change. The transition at 2 mM Mg(2+) is accompanied by substantial conformational changes revealed by probing with ribonucleases T1 and V1 and likely enhances stacking of the tRNA bases. Fast and slow kinetic phases of tRNA refolding were observed. Time-resolved monitoring of Mg(2+) binding to tRNA suggested that the slow kinetic phase was caused by a misfolded tRNA structure formed in the absence of Mg(2+). Our results suggest that, similarly to large RNAs, Mg(2+)-induced tRNA folding exhibits parallel folding pathways and the existence of kinetically trapped intermediates stabilized by Mg(2+). A multistep scheme for Mg(2+)-induced tRNA folding is discussed.     

Outlining folding nuclei in globular proteins

Our theoretical approach for prediction of folding/unfolding nuclei in three-dimensional protein structures is based on a search for free energy saddle points on networks of protein unfolding pathways. Under some approximations, this search is performed rapidly by dynamic programming and results in prediction of Phi values, which can be compared with those found experimentally. In this study, we optimize some details of the model (specifically, hydrogen atoms are taken into account in addition to heavy atoms), and compare the theoretically obtained and experimental Phi values (which characterize involvement of residues in folding nuclei) for all 17 proteins, where Phi values are now known for many residues. We show that the model provides good Phi value predictions for proteins whose structures have been determined by X-ray analysis (the average correlation coefficient is 0.65), with a more limited success for proteins whose structures have been determined by NMR techniques only (the average correlation coefficient is 0.34), and that the transition state free energies computed from the same model are in a good anticorrelation with logarithms of experimentally measured folding rates at mid-transition (the correlation coefficient is -0.73).     

Computer simulation of tRNA secondary structure folding

Computer simulation results of folding linear RNA moleculesinto secondaty structures are presented. The structure is formedby two interacting processes: the RNA molecular chain growth(beginning from an initial length, L_o), and the structuring(secondary structure sequential growth in the region of theexisting molecular chain, based on the local free energy minimizationby sequential addition of elementary substruc tures-stems).It was found that the final secondary structure formation isgreatly influenced by the ‘structuring period’ T(the ratio of the molecular chain growth rate to the structuringrate), and the direction of RNA synthesis. The computer simulationhas been performed for 219 and 906 tRNA genes from two publishedcatalogues, on the whale two-dimensional domain (T,L₀) parameters,by using four known free-energy models. Minimwn stem lengthand molecular chain growth direction have been also varied Thecalculated secondary structures have been compared to the naturaltRNA structures given in the catalogues, and the region of bestcoincidence for the model parameters has been determined. Ithas been proved that, on average, >86% of the paired basesof natural tRNA structures appear in the folding simulation.     

Kinetics of folding nuclei formation in proteins

When a protein folds or unfolds, it passes through many half-folded microstates. Only a few of them can accumulate and be seen experimentally, and this happens only when the folding (or unfolding) occurs far from the point of thermodynamic equilibrium between the native and denatured states. The universal features of folding, though, are observed in the vicinity of the equilibrium point. Here the "two-state" transition proceeds without any accumulation of metastable intermediates, and only the transition state ("folding nucleus") is outlined by its key influence on the folding/unfolding kinetics. This review covers recent experimental and theoretical studies of folding nuclei.     

Role of divalent ions in folding of tRNA.

The native structure of tRNA is not achieved in low salt (4.5 mM Na+, 25 degrees C), but can be restored by addition of divalent ions. We have explored the structure of the central region in Escherichia coli tRNAfMet by absorption and emission spectroscopy of 4-thiouracil, and the structure of the anticodon loop in yeast tRNAPhe by fluorescence of the 'Y' base, versus the number of manganese ions bound to tRNA, which was derived from electron spin resonance. The fluorescence of the reduced 8-13 photoproduct (in which 4-thiouracil at position 8 is crosslinked to cytosine at position 13) was also analysed. In low salt (e.g. 4.5 mM Na+), the region of 4-thiouracil is affected strongly as the first eight Mn2+ bind to tRNA, whereas the fluorescence of the 'Y' base is affected only after four Mn2+ are bound. Considering the structural similarities of the two tRNAs, this suggests that the reorganisation brought about by divalent ions starts in the central region, the anticodon loop being affected later. The binding of divalent ions to each region starts together with its restructuration. Monovalent ions can substitute for divalent ions in this process, a 15 mM sodium concentration being equivalent to the binding of the first five Mn2+. If divalent ions are then added, even the first ones distribute themselves between both the central and the anticodon region. Alternatively, the renaturation may be achieved by monovalent ions only, implying that no sites exist whose occupancy by divalent ions is crucial for the native structure. These observations suggest that the role and means of divalent ion binding to tRNA are largely explainable in terms of a simple maganese-phosphate binding supplemented by electrostatic interaction with distant phosphates.     

Cation-induced conformational changes in tRNA molecules

The uv absorption spectra and melting profiles of an initially ion-free solution of E. coli unfractionated tRNA are significantly modified by the addition of either Na⁺, Mg²⁺, or Mn²⁺ or of other first-series transition-metal ions such as Ni²⁺, Co²⁺, and Zn²⁺. The main effect of the addition of all monovalent or divalent cations examined is an increase of the ordered and stacking stabilized tRNA structure, as revealed by a drop in the absorption near 260 nm, as well as in the 4-TU absorption region. Sharp differences have, however, been detected in the 290–305-nm range in the presence of the various ions studied. When transition-metal ions were added to a tRNA solution, an absorption peak appeared at 294 nm. This effect is interpreted as a perturbation of the electronic structure of the bases due to direct binding of metal ions to the bases. An analysis of the variation in the spectrum as a function of metal concentration and of the thermal melting reversibility in the presence of various metal ions supports the conclusion that while all ions investigated are involved in binding to the phosphate groups of tRNA, transition-metal ions are also able to bind directly to the bases.     

Structural features of protein folding nuclei

A crucial event of protein folding is the formation of a folding nucleus. We demonstrate the presence of a considerable coincidence between the location of folding nuclei and the location of so-called "root structural motifs", which have unique overall folds and handedness. In the case of proteins with a single root structural motif, the involvement in the formation of a folding nucleus is in average significantly higher for amino acids residues that are in root structural motifs, compared to residues in other parts of the protein. The tests carried out revealed that the observed difference is statistically reliable. Thus, a structural feature that corresponds to the protein folding nucleus is now found.     

On prediction of folding nuclei in globular proteins

The approach described in this paper on the prediction of folding nuclei in globular proteins with known three dimensional structures is based on a search of the lowest saddle points through the barrier separating the unfolded state from the native structure on the free-energy landscape of protein chain. This search is performed by a dynamic programming method. Comparison of theoretical results with experimental data on the folding nuclei of two dozen of proteins shows that our model provides good phi value predictions for proteins whose structures have been determined by X-ray analysis, with a less limited success for proteins whose structures have been determined by NMR techniques only. Consideration of a full ensemble of transition states results in more successful prediction than consideration of only the transition states with the minimal free energy. In conclusion we have predicted the localization of folding nuclei for three dimensional protein structures for which kinetics of folding is studied now but the localization of folding nuclei is still unknown.     

On the optimal folding of protein molecules

The process of formation of a globular structure by a long molecular chain has been examined. In this process, various regions of the chain interact with one another. We classify the contacts thus formed as “correct” and “erroneous” ones. The correct contacts are those characteristic of the final native globular structure. All other contacts can be treated as erroneous. It is demonstrated that globule formation may proceed actually without formation and subsequent decay of erroneous contacts. Our model permits avoiding examination of numerous erroneous variants inasmuch as the regions of the chain that form correct contacts enter “long-range” interactions that at the same time can be highly selective. The existence of interactions of this kind facilitates the mutual approach and interaction of just those regions of the chain that yield correct contacts. Based on database analysis, it is shown that the model is valid not only for abstract structures but also for real polypeptide chains capable of forming protein globules and helical fibrils.     

Loop folds in proteins and evolutionary conservation of folding nuclei

We show that loops of close contacts involving hydrophobic residues are important in protein folding. Contrary to Berezovsky Berezovsky and Trifonov (J Biomol Struct Dyn 20, 5-6, 2002) the loops important in protein folding usually are much larger in size than 23-31 residues, being instead comparable to the size of the protein for single domain proteins. Additionally what is important are not single loop contacts, but a highly interconnected network of such loop contacts, which provides extra stability to a protein fold and which leads to their conservation in evolution.     

Prediction of protein folding pathways.

Recent 1H nuclear magnetic resonance (n.m.r.) hydrogen exchange experiments on five different proteins have delineated the secondary structures formed in trapped, partially folded intermediates. The early forming structural elements are identifiable through a technique described in this work to predict folding pathways. The method assumes that the sequential selection of structural fragments such as alpha-helices and beta-strands involved in the folding process is founded upon the maximal burial of solvent accessible surface from both the formation of internal structure and substructure association. The substructural elements were defined objectively by major changes in main-chain direction. The predicted folding pathways are in complete correspondence with the n.m.r. results in that the formed structural fragments found in the folding intermediates are those predicted earliest in the pathways. The technique was also applied to proteins of known tertiary structure and with fold similar to one of the five proteins examined by 1H n.m.r. The pathways for these structures also showed general consistency with the n.m.r. observations, suggesting conservation of a secondary structural framework or molten globule about which folding nucleates and proceeds.     

Concatemerization of tRNA molecules in the presence of trivaline derivative

The interaction of tRNA with trivaline dansyl hydrazide trifluoroacetate (DHTV) has been studied. The shape of curves of fluorimetric titration of tRNA with DHTV and vice versa can be explained only by formation of DHTV dimers on tRNA molecules, and subsequent association of DHTV-saturated tRNA molecules with each other. The ability of tRNA molecules to form concatemers in solution in the presence of DHTV has been demonstrated by electron microscopy. Electron microscopy of the tRNA-DHTV complexes stained with uranyl acetate revealed flexible rods 6-7 nm thick and up to several micrometers long.     

Protein folding: looping from hydrophobic nuclei.

Protein structure can be viewed as a compact linear array of nearly standard size closed loops of 25-30 amino acid residues (Berezovsky et al., FEBS Letters 2000; 466: 283-286) irrespective of details of secondary structure. The end-to-end contacts in the loops are likely to be hydrophobic, which is a testable hypothesis. This notion could be verified by direct comparison of the loop maps with Kyte and Doolittle hydropathicity plots. This analysis reveals that most of the ends of the loops are hydrophobic, indeed. The same conclusion is reached on the basis of positional autocorrelation analysis of protein sequences of 23 fully sequenced bacterial genomes. Hydrophobic residues valine, alanine, glycine, leucine, and isoleucine appear preferentially at the 25-30 residues distance one from another. These observations open a new perspective in the understanding of protein structure and folding: a consecutive looping of the polypeptide chain with the loops ending primarily at hydrophobic nuclei.     

Prediction of folding mechanism for circular-permuted proteins

Recent theoretical and experimental studies have suggested that real proteins have sequences with sufficiently small energetic frustration that topological effects are central in determining the folding mechanism. A particularly interesting and challenging framework for exploring and testing the viability of these energetically unfrustrated models is the study of circular-permuted proteins. Here we present the results of the application of a topology-based model to the study of circular permuted SH3 and CI2, in comparison with the available experimental results. The folding mechanism of the permuted proteins emerging from our simulations is in very good agreement with the experimental observations. The differences between the folding mechanisms of the permuted and wild-type proteins seem then to be strongly related to the change in the native state topology.