Size, shape, and flexibility of RNA structures

Determination of sizes and flexibilities of RNA molecules is important in understanding the nature of packing in folded structures and in elucidating interactions between RNA and DNA or proteins. Using the coordinates of the structures of RNA in the Protein Data Bank we find that the size of the folded RNA structures, measured using the radius of gyration, $R_G$, follows the Flory scaling law, namely, $R_G =5.5 N^{1/3}$ \AA where N is the number of nucleotides. The shape of RNA molecules is characterized by the asphericity $\Delta$ and the shape $S$ parameters that are computed using the eigenvalues of the moment of inertia tensor. From the distribution of $\Delta$, we find that a large fraction of folded RNA structures are aspherical and the distribution of $S$ values shows that RNA molecules are prolate ($S>0$). The flexibility of folded structures is characterized by the persistence length $l_p$. By fitting the distance distribution function $P(r)$ to the worm-like chain model we extracted the persistence length $l_p$. We find that $l_p\approx 1.5 N^{0.33}$ \AA. The dependence of $l_p$ on $N$ implies the average length of helices should increases as the size of RNA grows. We also analyze packing in the structures of ribosomes (30S, 50S, and 70S) in terms of $R_G$, $\Delta$, $S$, and $l_p$. The 70S and the 50S subunits are more spherical compared to most RNA molecules. The globularity in 50S is due to the presence of an unusually large number (compared to 30S subunit) of small helices that are stitched together by bulges and loops. Comparison of the shapes of the intact 70S ribosome and the constituent particles suggests that folding of the individual molecules might occur prior to assembly.Comment: 28 pages, 8 figures, J. Chem. Phys. in pres

The ribosome structure controls and directs mRNA entry, translocation and exit dynamics

The protein-synthesizing ribosome undergoes large motions to effect the translocation of tRNAs and mRNA; here the domain motions of this system are explored with a coarse-grained elastic network model using normal mode analysis. Crystal structures are used to construct various model systems of the 70S complex with/without tRNA, elongation factor Tu and the ribosomal proteins. Computed motions reveal the well-known ratchet-like rotational motion of the large subunits, as well as the head rotation of the small subunit and the high flexibility of the L1 and L7/L12 stalks, even in the absence of ribosomal proteins. This result indicates that these experimentally observed motions during translocation are inherently controlled by the ribosomal shape and only partially dependent upon GTP hydrolysis. Normal mode analysis further reveals the mobility of A- and P-tRNAs to increase in the absence of the E-tRNA. In addition, the dynamics of the E-tRNA is affected by the absence of the ribosomal protein L1. The mRNA in the entrance tunnel interacts directly with helicase proteins S3 and S4, which constrain the mRNA in a clamp-like fashion, as well as with protein S5, which likely orients the mRNA to ensure correct translation. The ribosomal proteins S7, S11 and S18 may also be involved in assuring translation fidelity by constraining the mRNA at the exit site of the channel. The mRNA also interacts with the 16S 3’ end forming the Shine-Dalgarno complex at the initiation step; the 3’ end may act as a ‘hook’ to reel in the mRNA to facilitate its exit

Probing complex RNA structures by mechanical force

RNA secondary structures of increasing complexity are probed combining single molecule stretching experiments and stochastic unfolding/refolding simulations. We find that force-induced unfolding pathways cannot usually be interpretated by solely invoking successive openings of native helices. Indeed, typical force-extension responses of complex RNA molecules are largely shaped by stretching-induced, long-lived intermediates including non-native helices. This is first shown for a set of generic structural motifs found in larger RNA structures, and then for Escherichia coli's 1540-base long 16S ribosomal RNA, which exhibits a surprisingly well-structured and reproducible unfolding pathway under mechanical stretching. Using out-of-equilibrium stochastic simulations, we demonstrate that these experimental results reflect the slow relaxation of RNA structural rearrangements. Hence, micromanipulations of single RNA molecules probe both their native structures and long-lived intermediates, so-called "kinetic traps", thereby capturing -at the single molecular level- the hallmark of RNA folding/unfolding dynamics.Comment: 9 pages, 9 figure

Protein folding on the ribosome studied using NMR spectroscopy

NMR spectroscopy is a powerful tool for the investigation of protein folding and misfolding, providing a characterization of molecular structure, dynamics and exchange processes, across a very wide range of timescales and with near atomic resolution. In recent years NMR methods have also been developed to study protein folding as it might occur within the cell, in a de novo manner, by observing the folding of nascent polypeptides in the process of emerging from the ribosome during synthesis. Despite the 2.3 MDa molecular weight of the bacterial 70S ribosome, many nascent polypeptides, and some ribosomal proteins, have sufficient local flexibility that sharp resonances may be observed in solution-state NMR spectra. In providing information on dynamic regions of the structure, NMR spectroscopy is therefore highly complementary to alternative methods such as X-ray crystallography and cryo-electron microscopy, which have successfully characterized the rigid core of the ribosome particle. However, the low working concentrations and limited sample stability associated with ribosome-nascent chain complexes means that such studies still present significant technical challenges to the NMR spectroscopist. This review will discuss the progress that has been made in this area, surveying all NMR studies that have been published to date, and with a particular focus on strategies for improving experimental sensitivity

Probing the Assembly of the Ribosome: Insights from Computational Studies on Ribosomal Proteins

Ribosomes are complex cellular machines that synthesize new proteins in the cell. The accurate and efficient assembly of ribosomal proteins (r-proteins) and ribosomal RNA (rRNA) to form a functional ribosome is important for cell growth, metabolic reactions, and other cellular processes. Ribosomal assembly has been an active research topic for many years because understanding the assembly mechansims can provide insight into protein/RNA recognitions that are important in many other cellular processes, as well as help optimize the development of antibacterial therapeutics. Experimental and computational sutdies thus far have greatly improved our understanding of assembly, yet many questions remain unanswered regarding the complex behaviors of r-proteins and rRNA during the process. To further understand ribosome assembly, we have computationally studied the sequences, structures, and dynamic properties of r-proteins from the 30S subunit and their relationships to RNA binding. We discuss the statistically greater amount of positively charged residues in r-proteins compared to other housekeeping proteins and observe a high level of charged interactions between r-proteins and rRNA in the assembled structure. We also detect a significant correlation between the overall flexibility of a protein and the number of contact points it makes with its rRNA binding site. Protein residues contacting with rRNA are observed to be more mobile in solution when compared to the non-contacting residues. We also describe common modes of structural dynamics, revealing likely conformational changes the proteins make prior to binding, how they relate to possible binding mechanisms used during the assembly and to the location of the protein in the fully assembled ribosome

A Computational Investigation on the Connection between Dynamics Properties of Ribosomal Proteins and Ribosome Assembly

Assembly of the ribosome from its protein and RNA constituents has been studied extensively over the past 50 years, and experimental evidence suggests that prokaryotic ribosomal proteins undergo conformational changes during assembly. However, to date, no studies have attempted to elucidate these conformational changes. The present work utilizes computational methods to analyze protein dynamics and to investigate the linkage between dynamics and binding of these proteins during the assembly of the ribosome. Ribosomal proteins are known to be positively charged and we find the percentage of positive residues in r-proteins to be about twice that of the average protein: Lys+Arg is 18.7% for E. coli and 21.2% for T. thermophilus. Also, positive residues constitute a large proportion of RNA contacting residues: 39% for E. coli and 46% for T. thermophilus. This affirms the known importance of charge-charge interactions in the assembly of the ribosome. We studied the dynamics of three primary proteins from E. coli and T. thermophilus 30S subunits that bind early in the assembly (S15, S17, and S20) with atomic molecular dynamic simulations, followed by a study of all r-proteins using elastic network models. Molecular dynamics simulations show that solvent-exposed proteins (S15 and S17) tend to adopt more stable solution conformations than an RNA-embedded protein (S20). We also find protein residues that contact the 16S rRNA are generally more mobile in comparison with the other residues. This is because there is a larger proportion of contacting residues located in flexible loop regions. By the use of elastic network models, which are computationally more efficient, we show that this trend holds for most of the 30S r-proteins

Dynamics and Driving Forces of Macromolecular Complexes

Many functions in living cells are governed by macromolecular complexes. To fully describe the underlying mechanisms, they have to be understood at atomic level. The present study combines data obtained by X-ray crystallography and cryo-electron microscopy (cryo-EM) with molecular dynamics (MD) simulations. Two functions of macromolecular complexes, the downregulation of neurotransmitter release by the SNARE protein complex under oxidative stress and the translocation of transfer RNAs (tRNAs) through the ribosome during protein synthesis, were investigated. First, the hypothesis that oxidation of two cysteines on linker of the SNARE protein SNAP-25B and consequent disulfide bond formation shortens this linker sufficiently to hinder complex formation was tested. For this purpose, MD simulations of the SNARE complex with and without the disulfide bond were compared. Disulfide bond formation lead to conformational changes of the linker and of three central hydrophobic layers necessary to form the SNARE complex. Previously, mutations of residues contributing to these layers have been shown to reduce neurotransmitter release, suggesting that the stability of these layers is crucial for complex formation. The results from the simulations agree with the hypothesis that disulfide bond formation leads to a destabilization of the SNARE complex thus rendering it dysfunctional. This mechanism is interpreted as a chemomechanical regulation to shut down neurotransmitter release under oxidative stress, which has been linked to neurodegenerative diseases. In a second part I investigated the ribosome, where after peptide bond formation, bound tRNAs translocate by more than 7 nm to adjacent binding sites (A and P to P and E), accompanied by large-scale conformational motions (L1-stalk, intersubunit rotation) of the ribosome. By combining existing cryo-EM reconstructions of translocation intermediates with high resolution crystal structures, we obtained 13 near-atomic resolution structures. Subsequently, MD simulations of were carried out for each intermediate state. The obtained dynamics within these states allowed to estimate transition rates between states for motions of the L1-stalk, tRNAs and intersubunit rotations. Rapid motions of the L1-stalk and the small (30S) subunit on sub-microsecond timescales were revealed, whereas tRNA motions were seen to be rate-limiting for most transitions. By calculating the interaction free energy between L1-stalk and tRNA, molecular forces were derived showing that the L1-stalk is actively pulling the tRNA from P to E binding site, thereby overcoming barriers for the tRNA motion. Further, ribosomal proteins L5 and L16 guide the tRNAs by 'sliding' and 'stepping' mechanisms involving key protein-tRNA contacts. This explains how tRNA binding affinity is kept sufficiently constant to allow rapid translocation despite large-scale displacements. Translocation is accompanied by rotations of the 30S ribosomal subunit of more than 20 degrees relative to the large (50S) subunit. For each translocation intermediate, the affinity of the two subunits with each other must be finely tuned enabling such conformational flexibility while maintaining integrity of the ribosomal complex. Analyzing the trajectories at residue level reveals two classes of intersubunit contact interactions: i) persistent residue contacts which are independent of 30S rotation and primarily located close to the axis of rotation. ii) contacts that are formed and ruptured depending on the rotation angle, seen mainly on the periphery. Strikingly, also these rotation specific contacts substantially contribute to the overall stability of the ribosomal assembly and are expected to maintain a constant interaction energy with low barriers for rotation. The simulations reveal that upon removal of tRNAs peripheral contacts are weakened and, in turn, intersubunit rotation angles decrease, in agreement with cryo-EM analysis of tRNA depleted ribosomes. The identified mechanisms for lowering free energy barriers and for fine-tuning affinities might have developed similarly in other macromolecular complexes